# American Institute of Mathematical Sciences

October  2013, 18(8): 2151-2174. doi: 10.3934/dcdsb.2013.18.2151

## Qualitative analysis of an age- and sex-structured vaccination model for human papillomavirus

 1 Department of Applied Mathematics and Sciences, Khalifa University of Science, Technology and Research, PO Box 127788, Abu Dhabi 2 Department of Mathematics, University of Manitoba, Winnipeg, Manitoba, R3T 2N2 3 Merck Research Laboratories, UG1C-60, PO Box 1000, North Wales, PA 19454-1099

Received  December 2011 Revised  February 2013 Published  July 2013

A new model for the transmission dynamics of human pappilomavirus (HPV) is designed and analysed. The model, which stratifies the total population in terms of age and gender, incorporates an imperfect anti-HPV vaccine with some therapeutic benefits. Rigorous qualitative analysis of the resulting age-structured model, which takes the form of a deterministic system of non-linear partial differential equations with separable transmission coefficients, shows that the disease-free equilibrium of the model is locally-asymptotically stable whenever the effective reproduction number (denoted by $\mathcal{R}_v$) is less than unity. It is shown to be globally-asymptotically stable if certain additional conditions hold. Furthermore, it is shown that the model has at least one endemic equilibrium when $\mathcal{R}_v$ exceeds unity. Hence, the effective control of HPV spread in a community, using a vaccine, is governed by the threshold quantity $\mathcal{R}_v$ (the use of the vaccine will lead to effective disease control or elimination only if it reduces the threshold quantity to a value less than unity; and the use of such vaccine will not lead to effective disease control if it fails to make the threshold quantity to be less than unity).
Citation: Tufail Malik, Abba Gumel, Elamin H. Elbasha. Qualitative analysis of an age- and sex-structured vaccination model for human papillomavirus. Discrete & Continuous Dynamical Systems - B, 2013, 18 (8) : 2151-2174. doi: 10.3934/dcdsb.2013.18.2151
##### References:
 [1] M. Al-arydah and R. J. Smith, An age-structured model of human papillomavirus vaccination,, Mathematics and Computers in Simulation, 82 (2011), 629. doi: 10.1016/j.matcom.2011.10.006. Google Scholar [2] R. Barnabas and G. Garnett, The potential public health impact of vaccines against human papillomavirus,, in, (2004), 61. Google Scholar [3] C. Castillo-Chavez and Z. Feng, Global stability of an age-structure model for TB and its applications to optimal vaccination strategies,, Math. Biosci., 151 (1998), 135. doi: 10.1016/S0025-5564(98)10016-0. Google Scholar [4] "Genital HPV Infection - CDC Fact Sheet,", Centers for Disease Control and Prevention (CDC)., , (). Google Scholar [5] H. W. Chesson, D. U. Ekwueme, M. Saraiya and L. E. Markowitz, Cost-effectiveness of human papillomavirus vaccination in the United States,, Emerg. Infect. Dis., 14 (2008), 244. doi: 10.3201/eid1402.070499. Google Scholar [6] M. E. Cruickshank, L. Sharp, G. Chambers, L. Smart and G. Murray, Persistent infection with human papillomavirus following the successful treatment of high grade cervical intraepithelial neoplasia,, BJOG, 109 (2002), 579. Google Scholar [7] O. Diekmann and J. A. P. Heesterbeek, "Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation,", John Wiley & Sons Ltd., (2000). Google Scholar [8] E. H. Elbasha, E. J. Dasback and R. P. Insinga, Model for assessing human papillomavirus vaccination strategies,, Emerg. Infect. Dis., 13 (2007), 28. doi: 10.3201/eid1301.060438. Google Scholar [9] E. H. Elbasha, Global stability of equilibria in a two-sex HPV vaccination model,, Bull. Math. Biol., 70 (2008), 894. doi: 10.1007/s11538-007-9283-0. Google Scholar [10] E. H. Elbasha, E. J. Dasbach and R. P. Insinga, A multi-type HPV transmission model,, Bull. Math. Biol., 70 (2008), 2126. doi: 10.1007/s11538-008-9338-x. Google Scholar [11] E. H. Elbasha and E. J. Dasbach, Impact of vaccinating boys and men against HPV in the United States,, Vaccine, 28 (2010), 6858. doi: 10.1016/j.vaccine.2010.08.030. Google Scholar [12] A. Ferenczy, Persistent human papillomavirus infection and cervical neoplasia,, The Lancet Oncology, 3 (2002), 11. doi: 10.1016/S1470-2045(01)00617-9. Google Scholar [13] D. Greenhalgh, Threshold and stability results for an epidemic model with an age-structured meeting rate,, IMA J. Math. Appl. Med. Biol., 5 (1988), 81. doi: 10.1093/imammb/5.2.81. Google Scholar [14] H. W. Hethcote, Age-structured epdiemiology models and expressions for $R_0$,, in, (2009), 91. doi: 10.1142/9789812834836_0003. Google Scholar [15] J. Hughes, G. Garnett and L. Koutsky, The theoretical population level impact of a prophylactic human papillomavirus vaccine,, Epidemiology, 13 (2002), 631. Google Scholar [16] H. Inaba, Threshold and stability results for an age-structured epidemic model,, J. Math. Biol., 28 (1990), 411. doi: 10.1007/BF00178326. Google Scholar [17] W. Kaplan, "Advanced Calculus,", Addison-Wesley Mathematics Series, (1984). Google Scholar [18] J. J. Kim and S. J. Goldie, Cost effectiveness analysis of including boys in a human papillomavirus vaccination programme in the United States,, BMJ, 339 (2009), 1. doi: 10.1136/bmj.b3884. Google Scholar [19] C. J. N. Lacey, C. M. Lowndes and K. V. Shah, Burden and management of non-cancerous HPV-related conditions: HPV-6/11 disease,, Vaccine, 24 (2006), 35. doi: 10.1016/j.vaccine.2006.06.015. Google Scholar [20] X. Li, G. Gupur and G. Zhu, Threshold and stability results for an age-structured SEIR epidemic model,, Comput. Math. Appl., 42 (2001), 883. doi: 10.1016/S0898-1221(01)00206-1. Google Scholar [21] X. Li and J. Liu, Stability of an age-structured epidemiological model for hepatitis C,, J. Appl. Math. Comput., 27 (2008), 159. doi: 10.1007/s12190-008-0060-5. Google Scholar [22] X. Li, J. Liu and M. Martcheva, An age-structured two-strain epidemic model with super-infection,, Math. Biosci. Engg., 7 (2009), 123. doi: 10.3934/mbe.2010.7.123. Google Scholar [23] A. T. Newall et al., Cost-effectiveness analyses of human papillomavirus vaccination,, The Lancet Infectious Diseases, 7 (2007), 289. Google Scholar [24] D. M. Parkin and F. Bray, The burden of HPV-related cancers,, Vaccine, 23 (2006), 11. doi: 10.1016/j.vaccine.2006.05.111. Google Scholar [25] D. M. Parkin, The global health burden of infection-associated cancers in the year 2002,, Int. J. Cancer, 118 (2006), 3030. doi: 10.1002/ijc.21731. Google Scholar [26] M. A. Safi et al., Qualitative analysis of an age-structred SEIR epidemic model with treatment,, Appl. Math. Comput., 219 (2013), 10627. doi: 10.1016/j.amc.2013.03.126. Google Scholar [27] N. F. Schlecht et al., Persistent human papillomavirus infection as a predictor of cervical intraepithelial neoplasia,, JAMA, 286 (2001), 3106. Google Scholar [28] L. L. Villa, R. L. R. Costa, C. A. Petta et al., Prophylactic quadrivalent human papillomavirus (types 6, 11, 16, and 18) L1 virus-like particle vaccine in young women: A randomised double-blind placebo-controlled multicentre phase II efficacy trial,, Lancet Oncol., 6 (2005), 271. doi: 10.1016/S1470-2045(05)70101-7. Google Scholar [29] L. Zou, S. Ruan and W. Zhang, An age-structured model for the transmission dynamics of Hepatitis B,, SIAM J. Appl. Math., 70 (2010), 3121. doi: 10.1137/090777645. Google Scholar

show all references

##### References:
 [1] M. Al-arydah and R. J. Smith, An age-structured model of human papillomavirus vaccination,, Mathematics and Computers in Simulation, 82 (2011), 629. doi: 10.1016/j.matcom.2011.10.006. Google Scholar [2] R. Barnabas and G. Garnett, The potential public health impact of vaccines against human papillomavirus,, in, (2004), 61. Google Scholar [3] C. Castillo-Chavez and Z. Feng, Global stability of an age-structure model for TB and its applications to optimal vaccination strategies,, Math. Biosci., 151 (1998), 135. doi: 10.1016/S0025-5564(98)10016-0. Google Scholar [4] "Genital HPV Infection - CDC Fact Sheet,", Centers for Disease Control and Prevention (CDC)., , (). Google Scholar [5] H. W. Chesson, D. U. Ekwueme, M. Saraiya and L. E. Markowitz, Cost-effectiveness of human papillomavirus vaccination in the United States,, Emerg. Infect. Dis., 14 (2008), 244. doi: 10.3201/eid1402.070499. Google Scholar [6] M. E. Cruickshank, L. Sharp, G. Chambers, L. Smart and G. Murray, Persistent infection with human papillomavirus following the successful treatment of high grade cervical intraepithelial neoplasia,, BJOG, 109 (2002), 579. Google Scholar [7] O. Diekmann and J. A. P. Heesterbeek, "Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation,", John Wiley & Sons Ltd., (2000). Google Scholar [8] E. H. Elbasha, E. J. Dasback and R. P. Insinga, Model for assessing human papillomavirus vaccination strategies,, Emerg. Infect. Dis., 13 (2007), 28. doi: 10.3201/eid1301.060438. Google Scholar [9] E. H. Elbasha, Global stability of equilibria in a two-sex HPV vaccination model,, Bull. Math. Biol., 70 (2008), 894. doi: 10.1007/s11538-007-9283-0. Google Scholar [10] E. H. Elbasha, E. J. Dasbach and R. P. Insinga, A multi-type HPV transmission model,, Bull. Math. Biol., 70 (2008), 2126. doi: 10.1007/s11538-008-9338-x. Google Scholar [11] E. H. Elbasha and E. J. Dasbach, Impact of vaccinating boys and men against HPV in the United States,, Vaccine, 28 (2010), 6858. doi: 10.1016/j.vaccine.2010.08.030. Google Scholar [12] A. Ferenczy, Persistent human papillomavirus infection and cervical neoplasia,, The Lancet Oncology, 3 (2002), 11. doi: 10.1016/S1470-2045(01)00617-9. Google Scholar [13] D. Greenhalgh, Threshold and stability results for an epidemic model with an age-structured meeting rate,, IMA J. Math. Appl. Med. Biol., 5 (1988), 81. doi: 10.1093/imammb/5.2.81. Google Scholar [14] H. W. Hethcote, Age-structured epdiemiology models and expressions for $R_0$,, in, (2009), 91. doi: 10.1142/9789812834836_0003. Google Scholar [15] J. Hughes, G. Garnett and L. Koutsky, The theoretical population level impact of a prophylactic human papillomavirus vaccine,, Epidemiology, 13 (2002), 631. Google Scholar [16] H. Inaba, Threshold and stability results for an age-structured epidemic model,, J. Math. Biol., 28 (1990), 411. doi: 10.1007/BF00178326. Google Scholar [17] W. Kaplan, "Advanced Calculus,", Addison-Wesley Mathematics Series, (1984). Google Scholar [18] J. J. Kim and S. J. Goldie, Cost effectiveness analysis of including boys in a human papillomavirus vaccination programme in the United States,, BMJ, 339 (2009), 1. doi: 10.1136/bmj.b3884. Google Scholar [19] C. J. N. Lacey, C. M. Lowndes and K. V. Shah, Burden and management of non-cancerous HPV-related conditions: HPV-6/11 disease,, Vaccine, 24 (2006), 35. doi: 10.1016/j.vaccine.2006.06.015. Google Scholar [20] X. Li, G. Gupur and G. Zhu, Threshold and stability results for an age-structured SEIR epidemic model,, Comput. Math. Appl., 42 (2001), 883. doi: 10.1016/S0898-1221(01)00206-1. Google Scholar [21] X. Li and J. Liu, Stability of an age-structured epidemiological model for hepatitis C,, J. Appl. Math. Comput., 27 (2008), 159. doi: 10.1007/s12190-008-0060-5. Google Scholar [22] X. Li, J. Liu and M. Martcheva, An age-structured two-strain epidemic model with super-infection,, Math. Biosci. Engg., 7 (2009), 123. doi: 10.3934/mbe.2010.7.123. Google Scholar [23] A. T. Newall et al., Cost-effectiveness analyses of human papillomavirus vaccination,, The Lancet Infectious Diseases, 7 (2007), 289. Google Scholar [24] D. M. Parkin and F. Bray, The burden of HPV-related cancers,, Vaccine, 23 (2006), 11. doi: 10.1016/j.vaccine.2006.05.111. Google Scholar [25] D. M. Parkin, The global health burden of infection-associated cancers in the year 2002,, Int. J. Cancer, 118 (2006), 3030. doi: 10.1002/ijc.21731. Google Scholar [26] M. A. Safi et al., Qualitative analysis of an age-structred SEIR epidemic model with treatment,, Appl. Math. Comput., 219 (2013), 10627. doi: 10.1016/j.amc.2013.03.126. Google Scholar [27] N. F. Schlecht et al., Persistent human papillomavirus infection as a predictor of cervical intraepithelial neoplasia,, JAMA, 286 (2001), 3106. Google Scholar [28] L. L. Villa, R. L. R. Costa, C. A. Petta et al., Prophylactic quadrivalent human papillomavirus (types 6, 11, 16, and 18) L1 virus-like particle vaccine in young women: A randomised double-blind placebo-controlled multicentre phase II efficacy trial,, Lancet Oncol., 6 (2005), 271. doi: 10.1016/S1470-2045(05)70101-7. Google Scholar [29] L. Zou, S. Ruan and W. Zhang, An age-structured model for the transmission dynamics of Hepatitis B,, SIAM J. Appl. Math., 70 (2010), 3121. doi: 10.1137/090777645. Google Scholar
 [1] C. Connell McCluskey. Global stability for an SEI epidemiological model with continuous age-structure in the exposed and infectious classes. Mathematical Biosciences & Engineering, 2012, 9 (4) : 819-841. doi: 10.3934/mbe.2012.9.819 [2] Andrea Franceschetti, Andrea Pugliese, Dimitri Breda. Multiple endemic states in age-structured $SIR$ epidemic models. Mathematical Biosciences & Engineering, 2012, 9 (3) : 577-599. doi: 10.3934/mbe.2012.9.577 [3] Dimitri Breda, Stefano Maset, Rossana Vermiglio. Numerical recipes for investigating endemic equilibria of age-structured SIR epidemics. Discrete & Continuous Dynamical Systems - A, 2012, 32 (8) : 2675-2699. doi: 10.3934/dcds.2012.32.2675 [4] Tufail Malik, Jody Reimer, Abba Gumel, Elamin H. Elbasha, Salaheddin Mahmud. The impact of an imperfect vaccine and pap cytology screening on the transmission of human papillomavirus and occurrence of associated cervical dysplasia and cancer. Mathematical Biosciences & Engineering, 2013, 10 (4) : 1173-1205. doi: 10.3934/mbe.2013.10.1173 [5] Najat Ziyadi. A male-female mathematical model of human papillomavirus (HPV) in African American population. Mathematical Biosciences & Engineering, 2017, 14 (1) : 339-358. doi: 10.3934/mbe.2017022 [6] Tom Burr, Gerardo Chowell. The reproduction number $R_t$ in structured and nonstructured populations. Mathematical Biosciences & Engineering, 2009, 6 (2) : 239-259. doi: 10.3934/mbe.2009.6.239 [7] Hui Cao, Yicang Zhou. The basic reproduction number of discrete SIR and SEIS models with periodic parameters. Discrete & Continuous Dynamical Systems - B, 2013, 18 (1) : 37-56. doi: 10.3934/dcdsb.2013.18.37 [8] Svend Christensen, Preben Klarskov Hansen, Guozheng Qi, Jihuai Wang. The mathematical method of studying the reproduction structure of weeds and its application to Bromus sterilis. Discrete & Continuous Dynamical Systems - B, 2004, 4 (3) : 777-788. doi: 10.3934/dcdsb.2004.4.777 [9] Ariel Cintrón-Arias, Carlos Castillo-Chávez, Luís M. A. Bettencourt, Alun L. Lloyd, H. T. Banks. The estimation of the effective reproductive number from disease outbreak data. Mathematical Biosciences & Engineering, 2009, 6 (2) : 261-282. doi: 10.3934/mbe.2009.6.261 [10] Nicolas Bacaër, Xamxinur Abdurahman, Jianli Ye, Pierre Auger. On the basic reproduction number $R_0$ in sexual activity models for HIV/AIDS epidemics: Example from Yunnan, China. Mathematical Biosciences & Engineering, 2007, 4 (4) : 595-607. doi: 10.3934/mbe.2007.4.595 [11] Ling Xue, Caterina Scoglio. Network-level reproduction number and extinction threshold for vector-borne diseases. Mathematical Biosciences & Engineering, 2015, 12 (3) : 565-584. doi: 10.3934/mbe.2015.12.565 [12] Gerardo Chowell, Catherine E. Ammon, Nicolas W. Hengartner, James M. Hyman. Estimating the reproduction number from the initial phase of the Spanish flu pandemic waves in Geneva, Switzerland. Mathematical Biosciences & Engineering, 2007, 4 (3) : 457-470. doi: 10.3934/mbe.2007.4.457 [13] Geni Gupur, Xue-Zhi Li. Global stability of an age-structured SIRS epidemic model with vaccination. Discrete & Continuous Dynamical Systems - B, 2004, 4 (3) : 643-652. doi: 10.3934/dcdsb.2004.4.643 [14] Diána H. Knipl, Gergely Röst. Modelling the strategies for age specific vaccination scheduling during influenza pandemic outbreaks. Mathematical Biosciences & Engineering, 2011, 8 (1) : 123-139. doi: 10.3934/mbe.2011.8.123 [15] Jacek Banasiak, Eddy Kimba Phongi, MirosŁaw Lachowicz. A singularly perturbed SIS model with age structure. Mathematical Biosciences & Engineering, 2013, 10 (3) : 499-521. doi: 10.3934/mbe.2013.10.499 [16] Gerardo Chowell, R. Fuentes, A. Olea, X. Aguilera, H. Nesse, J. M. Hyman. The basic reproduction number $R_0$ and effectiveness of reactive interventions during dengue epidemics: The 2002 dengue outbreak in Easter Island, Chile. Mathematical Biosciences & Engineering, 2013, 10 (5&6) : 1455-1474. doi: 10.3934/mbe.2013.10.1455 [17] Sophia R.-J. Jang. Discrete host-parasitoid models with Allee effects and age structure in the host. Mathematical Biosciences & Engineering, 2010, 7 (1) : 67-81. doi: 10.3934/mbe.2010.7.67 [18] Suxia Zhang, Xiaxia Xu. A mathematical model for hepatitis B with infection-age structure. Discrete & Continuous Dynamical Systems - B, 2016, 21 (4) : 1329-1346. doi: 10.3934/dcdsb.2016.21.1329 [19] Toshikazu Kuniya, Jinliang Wang, Hisashi Inaba. A multi-group SIR epidemic model with age structure. Discrete & Continuous Dynamical Systems - B, 2016, 21 (10) : 3515-3550. doi: 10.3934/dcdsb.2016109 [20] C. Connell McCluskey. Global stability for an $SEI$ model of infectious disease with age structure and immigration of infecteds. Mathematical Biosciences & Engineering, 2016, 13 (2) : 381-400. doi: 10.3934/mbe.2015008

2018 Impact Factor: 1.008