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.     

An age-structured within-host HIV-1 infection model with virus-to-cell and cell-to-cell transmissions

In this paper, a within-host HIV-1 infection model with virus-to-cell and direct cell-to-cell transmission and explicit age-since-infection structure for infected cells is investigated. It is shown that the model demonstrates a global threshold dynamics, fully described by the basic reproduction number. By analysing the corresponding characteristic equations, the local stability of an infection-free steady state and a chronic-infection steady state of the model is established. By using the persistence theory in infinite dimensional system, the uniform persistence of the system is established when the basic reproduction number is greater than unity. By means of suitable Lyapunov functionals and LaSalle's invariance principle, it is shown that if the basic reproduction number is less than unity, the infection-free steady state is globally asymptotically stable; if the basic reproduction number is greater than unity, the chronic-infection steady state is globally asymptotically stable. Numerical simulations are carried out to illustrate the feasibility of the theoretical results.     

Acquired immune response to oncogenic human papillomavirus associated with prophylactic cervical cancer vaccines

Human papillomavirus (HPV) is a common infection among women and a necessary cause of cervical cancer. Oncogenic HPV types infecting the anogenital tract have the potential to induce natural immunity, but at present we do not clearly understand the natural history of infection in humans and the mechanisms by which the virus can evade the host immune response. Natural acquired immune responses against HPV may be involved in the clearance of infection, but persistent infection with oncogenic virus types leads to the development of precancerous lesions and cancer. B cell responses are important for viral neutralization, but antibody responses in patients with cervical cancer are poor. Prophylactic vaccines targeting oncogenic virus types associated with cervical cancer have the potential to prevent up to 80% of cervical cancers by targeting HPV types 16 and 18. Clinical data show that prophylactic vaccines are effective in inducing antibody responses and in preventing persistent infection with HPV, as well as the subsequent development of high-grade cervical intraepithelial neoplasia. This article reviews the known data regarding natural immune responses to HPV and those developed by prophylactic vaccination.     

Genetic diversity and phylogenetic analysis of HPV 16 & 18 variants isolated from cervical specimens of women in Saudi Arabia

Human papillomavirus (HPV) are well known to be associated with the development of cervical cancer. HPV16 and HPV 18 are known as high-risk types and reported to be predominantly associated with cervical cancer. The prevalence and genetic diversity of HPV have been well documented globally but, in the Kingdom of Saudi Arabia, data on HPV genetic diversity are lacking. In this study, we have analyzed the genetic diversity of both HPV16 and HPV18 based on their L1 gene sequence because L1 gene is a major capsid protein gene and has been utilized to develop a prophylactic vaccine. In January 2011–2012, a total of forty samples from cervical specimens of women in Saudi Arabia were collected. The association of HPV16, HPV18 was detected by polymerase chain reaction, sequenced and submitted to GenBank. The sequences identity matrix and the phylogenetic relationship were analyzed with selected HPVs. The highest sequence identity (99.5%) for HPV16 and (99.3%) for HPV was observed with selected HPVs. The phylogenetic analysis results showed that HPVs from Saudi Arabia formed a closed cluster with African, Asian, East Asian as well as American HPVs distributed into multiple linages from various geographical locations. The results provided the valuable information about genetic diversity, but there is an urgent need to generate full genome sequence information which will provide a clearer picture of the genetic diversity and evolution of HPVs in Saudi Arabia. In conclusion, the generated data will be highly beneficial for developing molecular diagnostic tools, analyzing and correlating the epidemiological data to determine the risk of cervical cancer and finally to develop a vaccine for Saudi Arabian population.     

Potential Benefits of Second-Generation Human Papillomavirus Vaccines

Background

Current prophylactic vaccines against human papillomavirus (HPV) target two oncogenic types (16 and 18) that contribute to 70% of cervical cancer cases worldwide. Our objective was to quantify the range of additional benefits conferred by second-generation HPV prophylactic vaccines that are expected to expand protection to five additional oncogenic types (31, 33, 45, 52 and 58).

Methods

A microsimulation model of HPV and cervical cancer calibrated to epidemiological data from two countries (Kenya and Uganda) was used to estimate reductions in lifetime risk of cervical cancer from the second-generation HPV vaccines. We explored the independent and joint impact of uncertain factors (i.e., distribution of HPV types, co-infection with multiple HPV types, and unidentifiable HPV types in cancer) and vaccine properties (i.e., cross-protection against non-targeted HPV types), compared against currently-available vaccines.

Results

Assuming complete uptake of the second-generation vaccine, reductions in lifetime cancer risk were 86.3% in Kenya and 91.8% in Uganda, representing an absolute increase in cervical cancer reduction of 26.1% in Kenya and 17.9% in Uganda, compared with complete uptake of current vaccines. The range of added benefits was 19.6% to 29.1% in Kenya and 14.0% to 19.5% in Uganda, depending on assumptions of cancers attributable to multiple HPV infections and unidentifiable HPV types. These effects were blunted in both countries when assuming vaccine cross-protection with both the current and second-generation vaccines.

Conclusion

Second-generation HPV vaccines that protect against additional oncogenic HPV types have the potential to improve cervical cancer prevention. Co-infection with multiple HPV infections and unidentifiable HPV types can influence vaccine effectiveness, but the magnitude of effect may be moderated by vaccine cross-protective effects. These benefits must be weighed against the cost of the vaccines in future analyses.     

An SIR pairwise epidemic model with infection age and demography

The demography and infection age play an important role in the spread of slowly progressive diseases. To investigate their effects on the disease spreading, we propose a pairwise epidemic model with infection age and demography on dynamic networks. The basic reproduction number of this model is derived. It is proved that there is a disease-free equilibrium which is globally asymptotically stable if the basic reproduction number is less that unity. Besides, sensitivity analysis is performed and shows that increasing the variance in recovery time and decreasing the variance in infection time can effectively control the diseases. The complex interaction between the death rate and equilibrium prevalence suggests that it is imperative to correctly estimate the parameters of demography in order to assess the disease transmission dynamics accurately. Moreover, numerical simulations show that the endemic equilibrium is globally asymptotically stable.     

Global Properties of Infectious Disease Models with Nonlinear Incidence

We consider global properties for the classical SIR, SIRS and SEIR models of infectious diseases, including the models with the vertical transmission, assuming that the horizontal transmission is governed by an unspecified function f(S,I). We construct Lyapunov functions which enable us to find biologically realistic conditions sufficient to ensure existence and uniqueness of a globally asymptotically stable equilibrium state. This state can be either endemic, or infection-free, depending on the value of the basic reproduction number.     

Differential susceptibility epidemic models

We formulate compartmental differential susceptibility (DS) susceptible-infective-removed (SIR) models by dividing the susceptible population into multiple subgroups according to the susceptibility of individuals in each group. We analyze the impact of disease-induced mortality in the situations where the number of contacts per individual is either constant or proportional to the total population. We derive an explicit formula for the reproductive number of infection for each model by investigating the local stability of the infection-free equilibrium. We further prove that the infection-free equilibrium of each model is globally asymptotically stable by qualitative analysis of the dynamics of the model system and by utilizing an appropriately chosen Liapunov function. We show that if the reproductive number is greater than one, then there exists a unique endemic equilibrium for all of the DS models studied in this paper. We prove that the endemic equilibrium is locally asymptotically stable for the models with no disease-induced mortality and the models with contact numbers proportional to the total population. We also provide sufficient conditions for the stability of the endemic equilibrium for other situations. We briefly discuss applications of the DS models to optimal vaccine strategies and the connections between the DS models and predator-prey models with multiple prey populations or host-parasitic interaction models with multiple hosts are also given.This research was partially supported by the Department of Energy under contracts W-7405-ENG-36 and the Applied Mathematical Sciences Program KC-07-01-01.     

Modelling the effects of pre-exposure and post-exposure vaccines in tuberculosis control

Epidemic control strategies alter the spread of the disease in the host population. In this paper, we describe and discuss mathematical models that can be used to explore the potential of pre-exposure and post-exposure vaccines currently under development in the control of tuberculosis. A model with bacille Calmette-Guerin (BCG) vaccination for the susceptibles and treatment for the infectives is first presented. The epidemic thresholds known as the basic reproduction numbers and equilibria for the models are determined and stabilities are investigated. The reproduction numbers for the models are compared to assess the impact of the vaccines currently under development. The centre manifold theory is used to show the existence of backward bifurcation when the associated reproduction number is less than unity and that the unique endemic equilibrium is locally asymptotically stable when the associated reproduction number is greater than unity. From the study we conclude that the pre-exposure vaccine currently under development coupled with chemoprophylaxis for the latently infected and treatment of infectives is more effective when compared to the post-exposure vaccine currently under development for the latently infected coupled with treatment of the infectives.     

Vaccination against multiple HPV types

Vaccines against the most common human papillomavirus (HPV) types are currently under development. Epidemiologic data suggest that the transmission dynamics of different HPV types are not independent. Some studies indicate that interactions among HPV types are synergistic, where infection with one type facilitates concurrent or subsequent infection with another HPV type. Other studies point to antagonistic interference among HPV types. Here we develop a mathematical model to explore how these interactions may either enhance or diminish the effectiveness of vaccination programs designed to reduce the prevalence of the HPV types associated with cervical cancer. We analyze the local stability of the infection-free and boundary equilibria and characterize the conditions leading to a coexistence equilibrium. We also illustrate the results with numerical simulations using realistic model parameters. We show that if interactions among HPV types are synergistic, mass vaccination may reduce the prevalence of types that are not even included in the vaccine.     

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.     

Chancroid transmission dynamics: a mathematical modeling approach

Mathematical models have long been used to better understand disease transmission dynamics and how to effectively control them. Here, a chancroid infection model is presented and analyzed. The disease-free equilibrium is shown to be globally asymptotically stable when the reproduction number is less than unity. High levels of treatment are shown to reduce the reproduction number suggesting that treatment has the potential to control chancroid infections in any given community. This result is also supported by numerical simulations which show a decline in chancroid cases whenever the reproduction number is less than unity.     

Dynamics of an SIQS epidemic model with transport-related infection and exit-entry screenings

Population dispersal, as a common phenomenon in human society, may cause the spreading of many diseases such as influenza, SARS, etc. which are easily transmitted from one region to other regions. Exit and entry screenings at the border are considered as effective ways for controlling the spread of disease. In this paper, the dynamics of an SIQS model are analyzed and the combined effects of transport-related infection enhancing and exit-entry screenings suppressing on disease spread are discussed. The basic reproduction number is computed and proved to be a threshold for disease control. If it is not greater than the unity, the disease free equilibrium is globally asymptotically stable. And there exists an endemic equilibrium which is locally asymptotically stable if the reproduction number is greater than unity. It is shown that the disease is endemic in the sense of permanence if and only if the endemic equilibrium exists. Exit screening and entry screening are shown to be helpful for disease eradication since they can always have the possibility to eradicate the disease endemic led by transport-related infection and furthermore have the possibility to eradicate disease even when the isolated cites are disease endemic.     

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.     

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.     

Global Properties of <Emphasis Type="Italic">SIR</Emphasis> and <Emphasis Type="Italic">SEIR</Emphasis> Epidemic Models with Multiple Parallel Infectious Stages

We consider global properties of compartment SIR and SEIR models of infectious diseases, where there are several parallel infective stages. For instance, such a situation may arise if a fraction of the infected are detected and treated, while the rest of the infected remains undetected and untreated. We assume that the horizontal transmission is governed by the standard bilinear incidence rate. The direct Lyapunov method enables us to prove that the considered models are globally stable: There is always a globally asymptotically stable equilibrium state. Depending on the value of the basic reproduction number R ₀, this state can be either endemic (R ₀>1), or infection-free (R ₀≤1).     

Global Dynamics of an In-host Viral Model with Intracellular Delay

The dynamics of a general in-host model with intracellular delay is studied. The model can describe in vivo infections of HIV-I, HCV, and HBV. It can also be considered as a model for HTLV-I infection. We derive the basic reproduction number R ₀ for the viral infection, and establish that the global dynamics are completely determined by the values of R ₀. If R ₀≤1, the infection-free equilibrium is globally asymptotically stable, and the virus are cleared. If R ₀>1, then the infection persists and the chronic-infection equilibrium is locally asymptotically stable. Furthermore, using the method of Lyapunov functional, we prove that the chronic-infection equilibrium is globally asymptotically stable when R ₀>1. Our results shows that for intercellular delays to generate sustained oscillations in in-host models it is necessary have a logistic mitosis term in target-cell compartments.     

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.     

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.