We propose a fractional-order differential equation model to study the gender-specific transmission dynamics of human papillomavirus (HPV). Due to limited gender-disaggregated HPV data, oropharyngeal cancer incidence was used as a proxy. The model includes a theoretical analysis of disease-free and endemic equilibria and computation of the basic reproduction number, with stability conditions derived for both cases. Structural identifiability was confirmed using differential algebra, and data fitting indicated that the fractional models $ (\alpha = 0.9) $ better captured U.S. dynamics, while classical models $ (\alpha = 1) $ suited Turkish data. Practical identifiability, assessed via Monte Carlo simulations, revealed that most parameters were identifiable, though some (e.g., the male exit rate $ \mu_m $) showed limitations. Numerical simulations demonstrated sensitivity to fractional effects. This work provides a robust framework for modeling the transmission of HPV and offers insights for future epidemiological modeling and intervention planning.
Citation: Veysel Fuat Hatipoğlu, Maia Martcheva. Fractional gender structured model of human papillomavirus (HPV)[J]. Mathematical Biosciences and Engineering, 2026, 23(6): 1708-1742. doi: 10.3934/mbe.2026062
We propose a fractional-order differential equation model to study the gender-specific transmission dynamics of human papillomavirus (HPV). Due to limited gender-disaggregated HPV data, oropharyngeal cancer incidence was used as a proxy. The model includes a theoretical analysis of disease-free and endemic equilibria and computation of the basic reproduction number, with stability conditions derived for both cases. Structural identifiability was confirmed using differential algebra, and data fitting indicated that the fractional models $ (\alpha = 0.9) $ better captured U.S. dynamics, while classical models $ (\alpha = 1) $ suited Turkish data. Practical identifiability, assessed via Monte Carlo simulations, revealed that most parameters were identifiable, though some (e.g., the male exit rate $ \mu_m $) showed limitations. Numerical simulations demonstrated sensitivity to fractional effects. This work provides a robust framework for modeling the transmission of HPV and offers insights for future epidemiological modeling and intervention planning.
| [1] | HPV and Cancer, National Cancer Institute. Avaiable from: https://www.cancer.gov/about-cancer/causes-prevention/risk/infectious-agents/hpv-and-cancer. |
| [2] | Sexually Transmitted Infections. Avaiable from: https://www.cdc.gov/sti/about/about-genital-hpv-infection.html. |
| [3] |
V. Hancer, M. Buyukdogan, I. Bylykbashi, B. Oksuz, M. Acar, Prevalence of human papilloma virus types in Turkish and Albanian women, J. Cytol., 35 (2018), 252–254. https://doi.org/10.4103/JOC.JOC_162_17 doi: 10.4103/JOC.JOC_162_17
|
| [4] |
R. M. Lewis, J. Laprise, J. Gargano, E. Unger, T. Querec, H. Chesson, et al., Estimated prevalence and incidence of disease-associated human papillomavirus types among 15-to 59-year-olds in the United States, Sex. Transm. Dis., 48 (2021), 273–277. https://doi.org/10.1097/OLQ.0000000000001356 doi: 10.1097/OLQ.0000000000001356
|
| [5] | Human Papillomavirus and Related Cancers, Fact Sheet 2023, ICO/IARC Information Centre on HPV and Cancer Fact Sheet 2023. Avaiable from: https://hpvcentre.net/statistics/reports/TUR_FS.pdf. |
| [6] |
A. Bayram, Y. Derici, N. Yilmaz, S. Hanci, N. Agus, M. Sirin, et al., Prevalence of high-risk human papillomavirus in women from Turkey, Clin. Obstet. Gynecol. Reprod. Med., 1 (2015), 84–86. https://doi.org/10.15761/COGRM.1000122 doi: 10.15761/COGRM.1000122
|
| [7] | HPV and Related Cancers, ICO/IARC Information Center on HPV and Cancer, Fact Sheet 2023. Avaiable from: https://hpvcentre.net/statistics/reports/USA_FS.pdf. |
| [8] |
S. Lee, A. Tameru, A mathematical model of human papillomavirus (HPV) in the United States and its impact on cervical cancer, J. Cancer, 3 (2012), 262–268. https://doi.org/10.7150/jca.4161 doi: 10.7150/jca.4161
|
| [9] |
A. Sado, Mathematical modeling of cervical cancer with HPV transmission and vaccination, Sci. J. Appl. Mat. Stat., 7 (2019), 21–25. https://doi.org/10.11648/j.sjams.20190702.13 doi: 10.11648/j.sjams.20190702.13
|
| [10] |
K. Zhang, X. Wang, H. Liu, Y. Ji, Q. Pan, Y. Wei, et al., Mathematical analysis of a human papillomavirus transmission model with vaccination and screening, Math. Biosci. Eng., 17 (2020), 5449–5476. https://doi.org/10.3934/mbe.2020294 doi: 10.3934/mbe.2020294
|
| [11] |
E. J. Dasbach, E. H. Elbasha, R. P. Insinga, Mathematical models for predicting the epidemiologic and economic impact of vaccination against human papillomavirus infection and disease, Epidemiol. Rev., 28 (2006), 88–100. https://doi.org/10.1093/epirev/mxj006 doi: 10.1093/epirev/mxj006
|
| [12] |
T. Malik, J. Reimer, A. Gumel, E. H. Elbasha, S. 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, Math. Biosci. Eng., 10 (2013), 1173–1205. https://doi.org/10.3934/mbe.2013.10.1173 doi: 10.3934/mbe.2013.10.1173
|
| [13] |
H. D. Desta, G. T. Tilahun, T. M. Tolasa, M. G. Geleso, Mathematical model of human papillomavirus (HPV) dynamics with double-dose vaccination and its impact on cervical cancer, Discrete Dyn. Nat. Soc., 2024 (2024), 9971859. https://doi.org/10.1155/ddns/9971859 doi: 10.1155/ddns/9971859
|
| [14] |
C. Liu, H. Liu, X. Zhu, X. Lin, Q. Zhang, Y. Wei, Dynamic analysis of human papillomavirus transmission model under vaccine intervention: A case study of cervical cancer patients from Hungary, Adv. Cont. Discr. Mod., 2024 (2024), 36. https://doi.org/10.1186/s13662-024-03838-z doi: 10.1186/s13662-024-03838-z
|
| [15] |
E. H. Elbasha, Global stability of equilibria in a two-sex HPV vaccination model, Bull. Math. Biol., 70 (2008), 894–909. https://doi.org/10.1007/s11538-007-9283-0 doi: 10.1007/s11538-007-9283-0
|
| [16] |
A. Omame, R. Umana, D. Okuonghae, S. Inyama, Mathematical analysis of a two-sex human papillomavirus (HPV) model, Int. J. Biomath., 11 (2018), 1850092. https://doi.org/10.1142/S1793524518500924 doi: 10.1142/S1793524518500924
|
| [17] |
O. Sharomi, T. Malik, A model to assess the effect of vaccine compliance on human papillomavirus infection and cervical cancer, Appl. Math. Modell., 47 (2017), 528–550. https://doi.org/10.1016/j.apm.2017.03.025 doi: 10.1016/j.apm.2017.03.025
|
| [18] |
S. Gao, M. Martcheva, H. Miao, L. Rong, A dynamic model to assess human papillomavirus vaccination strategies in a heterosexual population combined with men who have sex with men, Bull. Math. Biol., 83 (2021), 1–36. https://doi.org/10.1007/s11538-020-00830-y doi: 10.1007/s11538-020-00830-y
|
| [19] |
S. Gao, M. Martcheva, H. Miao, L. Rong, The impact of vaccination on human papillomavirus infection with disassortative geographical mixing: A two-patch modeling study, J. Math. Biol., 84 (2022), 43. https://doi.org/10.1007/s00285-022-01745-z doi: 10.1007/s00285-022-01745-z
|
| [20] |
S. Gao, M. Martcheva, H. Miao, L. Rong, A two-sex model of human papillomavirus infection: Vaccination strategies and a case study, J. Theoret. Biol., 536 (2022), 111006. https://doi.org/10.1016/j.jtbi.2022.111006 doi: 10.1016/j.jtbi.2022.111006
|
| [21] |
F. Saldaña, J. A. Camacho-Gutiérrez, G. Villavicencio-Pulido, J. X. Velasco-Hernandez, Modeling the transmission dynamics and vaccination strategies for human papillomavirus infection: An optimal control approach, Appl. Math. Modell., 112 (2022), 767–785. https://doi.org/10.1016/j.apm.2022.08.017 doi: 10.1016/j.apm.2022.08.017
|
| [22] |
M. Al-arydah, R. Smith, An age-structured model of human papillomavirus vaccination, Math. Comput. Simul., 82 (2011), 629–652. https://doi.org/10.1016/j.matcom.2011.10.006 doi: 10.1016/j.matcom.2011.10.006
|
| [23] |
Z. Afzal, M. Alshehri, Fractional-order SIR model for ADHD as a neurobiological and genetic disorder, Sci. Rep., 15 (2025), 22992. https://doi.org/10.1038/s41598-025-07646-7 doi: 10.1038/s41598-025-07646-7
|
| [24] |
B. C. Agbata, R. Dervishi, D. F. Agbebaku, E. Cenaj, O. C. Collins, A. U. Ezeafulukwe, et al., A comprehensive analysis of fractional-order model of tuberculosis with treatment intervention, BMC Infect. Dis., 25 (2025), 1070. https://doi.org/10.1186/s12879-025-11303-9 doi: 10.1186/s12879-025-11303-9
|
| [25] |
Y. Chen, F. Liu, Q. Yu, T. Li, Review of fractional epidemic models, Appl. Math. Modell., 97 (2021), 281–307. https://doi.org/10.1016/j.apm.2021.03.044 doi: 10.1016/j.apm.2021.03.044
|
| [26] |
Z. U. A. Zafar, M. T. Hussain, M. Inc, D. Baleanu, B. Almohsen, A. S. Oke, et al., Fractional-order dynamics of human papillomavirus, Results Phys., 34 (2022), 105281. https://doi.org/10.1016/j.rinp.2022.105281 doi: 10.1016/j.rinp.2022.105281
|
| [27] |
U. K. Nwajeri, A. Omame, C. P. Onyenegecha, Analysis of a fractional order model for HPV and CT co-infection, Results Phys., 28 (2021), 104643. https://doi.org/10.1016/j.rinp.2021.104643 doi: 10.1016/j.rinp.2021.104643
|
| [28] |
A. El-Mesady, T. M. Al-shami, H. M. Ali, Optimal control efforts to reduce the transmission of HPV in a fractional-order mathematical model, Bound. Value Probl., 2025 (2025), 1–39. https://doi.org/10.1186/s13661-024-01991-8 doi: 10.1186/s13661-024-01991-8
|
| [29] |
M. C. Bahi, S. Bahramand, R. Jan, S. Boulaaras, H. Ahmad, R. Guefaifia, Fractional view analysis of sexual transmitted human papilloma virus infection for public health, Sci. Rep., 14 (2024), 3048. https://doi.org/10.1038/s41598-024-53696-8 doi: 10.1038/s41598-024-53696-8
|
| [30] |
B. Bajjah, M. Modanli, Finite difference method for infection model of HPV with cervical cancer under Caputo operator, Discrete Dyn. Nat. Soc., 2024 (2024), 2580745. https://doi.org/10.1155/2024/2580745 doi: 10.1155/2024/2580745
|
| [31] |
N. Raza, A. Raza, Y. Chahlaoui, J. Gomez-Aguilar, Numerical analysis of hpv HPV and its association with cervical cancer using Atangana–Baleanu fractional derivative, Modell. Earth Syst. Environ., 11 (2025), 60. https://doi.org/10.1007/s40808-024-02243-5 doi: 10.1007/s40808-024-02243-5
|
| [32] |
P. K. Rajan, M. Kuppusamy, A. Yusuf, A fractional-order modeling of human papillomavirus transmission and cervical cancer, Model. Earth Syst. Environ., 10 (2024), 1337–1357. https://doi.org/10.1007/s40808-023-01843-x doi: 10.1007/s40808-023-01843-x
|
| [33] |
P. K. Rajan, M. Kuppusamy, A fractional order human papillomavirus model with Caputo derivative, J. Anal., 32 (2024), 2135–2156. https://doi.org/10.1007/s41478-023-00641-z doi: 10.1007/s41478-023-00641-z
|
| [34] |
K. Hattaf, A new mixed fractional derivative with applications in computational biology, Computation, 12 (2024), 7. https://doi.org/10.3390/computation12010007 doi: 10.3390/computation12010007
|
| [35] |
K. Hattaf, Useful results for the qualitative analysis of generalized Hattaf mixed fractional differential equations with applications to medicine, Computation, 13 (2025), 167. https://doi.org/10.3390/computation13070167 doi: 10.3390/computation13070167
|
| [36] |
F. Assadiki, K. Hattaf, N. Yousfi, Global stability of fractional partial differential equations applied to the biological system modeling a viral infection with Hattaf time-fractional derivative, Math. Model. Comput., 11 (2024), 430–437. https://doi.org/10.23939/mmc2024.02.430 doi: 10.23939/mmc2024.02.430
|
| [37] |
S. Bhatter, S. Kumawat, S. D. Purohit, D. Suthar, Mathematical modeling of tuberculosis using Caputo fractional derivative: A comparative analysis with real data, Sci. Rep., 15 (2025), 12672. https://doi.org/10.1038/s41598-025-97502-5 doi: 10.1038/s41598-025-97502-5
|
| [38] |
E. Kharazmi, M. Cai, X. Zheng, Z. Zhang, G. Lin, G. E. Karniadakis, Identifiability and predictability of integer and fractional-order epidemiological models using physics-informed neural networks, Nat. Comput. Sci., 1 (2021), 744–753. https://doi.org/10.1038/s43588-021-00158-0 doi: 10.1038/s43588-021-00158-0
|
| [39] |
A. Morán-Torres, N. G. Pazos-Salazar, S. Téllez-Lorenzo, R. Jiménez-Lima, M. Lizano, D. O. Reyes-Hernández, et al., HPV oral and oropharynx infection dynamics in young population, Braz. J. Microbiol., 52 (2021), 1991–2000. https://doi.org/10.1007/s42770-021-00602-3 doi: 10.1007/s42770-021-00602-3
|
| [40] |
S. L. Ranjeva, E. B. Baskerville, V. Dukic, L. L. Villa, E. Lazcano-Ponce, A. R. Giuliano, et al., Recurring infection with ecologically distinct HPV types can explain high prevalence and diversity, Proc. Natl. Acad. Sci. U.S.A., 114 (2017), 13573–13578. https://doi.org/10.1073/pnas.1714712114 doi: 10.1073/pnas.1714712114
|
| [41] |
H. Trottier, S. Ferreira, P. Thomann, M. C. Costa, J. S. Sobrinho, J. C. M. Prado, et al., Human papillomavirus infection and reinfection in adult women: The role of sexual activity and natural immunity, Cancer. Res., 70 (2010), 8569–8577. https://doi.org/10.1158/0008-5472.CAN-10-0621 doi: 10.1158/0008-5472.CAN-10-0621
|
| [42] | Wikipedia Contributors. Caputo Fractional Derivative—Wikipedia, The Free Encyclopedia. Avaiable from: https://en.wikipedia.org/w/index.php?title=Caputo_fractional_derivative&oldid=1274649699. |
| [43] |
H. Jahanshahi, J. M. Munoz-Pacheco, S. Bekiros, N. D. Alotaibi, A fractional-order SIRD model with time-dependent memory indexes for encompassing the multi-fractional characteristics of the COVID-19, Chaos Solitons Fractals, 143 (2021), 110632. https://doi.org/10.1016/j.chaos.2020.110632 doi: 10.1016/j.chaos.2020.110632
|
| [44] |
J. M. Munoz-Pacheco, C. Posadas-Castillo, E. Zambrano-Serrano, The effect of a non-local fractional operator in an asymmetrical glucose-insulin regulatory system: Analysis, synchronization and electronic implementation, Symmetry, 12 (2020), 1395. https://doi.org/10.3390/sym12091395 doi: 10.3390/sym12091395
|
| [45] | J. Gómez-Aguilar, J. Rosales-García, J. Bernal-Alvarado, T. Córdova-Fraga, R. Guzmán-Cabrera, Fractional mechanical oscillators, Rev. Mex. Fis., 58 (2012), 348–352. |
| [46] |
P. Van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
|
| [47] | D. Matignon, Stability results for fractional differential equations with applications to control processing, Comput. Eng. Syst. Appl., 2 (1996), 963–968. |
| [48] |
Z. M. Odibat, Analytic study on linear systems of fractional differential equations, Comput. Math. Appl., 59 (2010), 1171–1183. https://doi.org/10.1016/j.camwa.2009.06.035 doi: 10.1016/j.camwa.2009.06.035
|
| [49] |
C. Vargas-De-León, Volterra-type Lyapunov functions for fractional-order epidemic systems, Commun. Nonlinear Sci. Numer. Simul., 24 (2015), 75–85. https://doi.org/10.1016/j.cnsns.2014.12.013 doi: 10.1016/j.cnsns.2014.12.013
|
| [50] | J. LaSalle, Limiting equations and stability of nonautonomous ordinary differential equations, in The Stability of Dynamical Systems, Society for Industrial and Applied Mathematics, (1976), 57–76. |
| [51] |
J. A. Gallegos, M. A. Duarte-Mermoud, On the Lyapunov theory for fractional order systems, Appl. Math. Comput., 287 (2016), 161–170. https://doi.org/10.1016/j.amc.2016.04.039 doi: 10.1016/j.amc.2016.04.039
|
| [52] |
K. K. Ang, J. Harris, R. Wheeler, R. Weber, D. I. Rosenthal, P. F. Nguyen-Tân, et al., Human papillomavirus and survival of patients with oropharyngeal cancer, N. Engl. J. Med., 363 (2010), 24–35. https://doi.org/10.1056/NEJMoa0912217 doi: 10.1056/NEJMoa0912217
|
| [53] | A. K. Chaturvedi, E. A. Engels, R. M. Pfeiffer, B. Y. Hernandez, W. Xiao, E. Kim, et al., Human papillomavirus and rising oropharyngeal cancer incidence in the United States, J. Clin. Oncol., 29 (2011), 4294–4301. https://doi.org/0.1200/JCO.2011.36.4596 |
| [54] |
S. R. Khode, R. C. Dwivedi, P. Rhys-Evans, R. Kazi, Exploring the link between human papilloma virus and oral and oropharyngeal cancers, J. Cancer. Res. Ther., 10 (2014), 492–498. https://doi.org/10.4103/0973-1482.138213 doi: 10.4103/0973-1482.138213
|
| [55] | E. Kokkinis, N. S. Bastas, I. Mega, C. Tsironis, A. D. Lianou, Association of HPV with oral and oropharyngeal cancer: Current evidence, Maedica, 19 (2024), 801–806. |
| [56] |
M. Lechner, J. Liu, L. Masterson, T. R. Fenton, HPV-associated oropharyngeal cancer: Epidemiology, molecular biology and clinical management, Nat. Rev. Clin. Oncol., 19 (2022), 306–327. https://doi.org/10.1038/s41571-022-00603-7 doi: 10.1038/s41571-022-00603-7
|
| [57] |
H. Damgacioglu, K. Sonawane, Y. Zhu, R. Li, B. A. Balasubramanian, D. R. Lairson, et al., Oropharyngeal cancer incidence and mortality trends in all 50 states in the US, 2001–2017, JAMA Otolaryngol. Head Neck Surg., 148 (2022), 155–165. https://doi.org/10.1001/jamaoto.2021.3567 doi: 10.1001/jamaoto.2021.3567
|
| [58] | M. Ervik, F. Lam, M. Laversanne, M. Colombet, J. Ferlay, A. Miranda-Filho, et al., Global Cancer Observatory: Cancer Over Time. Lyon, France: International Agency for Research on Cancer. Avaiable from: https://gco.iarc.fr/overtime. |
| [59] |
O. Chiş, J. R. Banga, E. Balsa-Canto, GenSSI: A software toolbox for structural identifiability analysis of biological models, Bioinformatics, 27 (2011), 2610–2611. https://doi.org/10.1093/bioinformatics/btr431 doi: 10.1093/bioinformatics/btr431
|
| [60] |
H. Hong, A. Ovchinnikov, G. Pogudin, C. Yap, SIAN: Software for structural identifiability analysis of ODE models, Bioinformatics, 35 (2019), 2873–2874. https://doi.org/10.1093/bioinformatics/bty1069 doi: 10.1093/bioinformatics/bty1069
|
| [61] |
G. Bellu, M. P. Saccomani, S. Audoly, L. D'Angiò, DAISY: A new software tool to test global identifiability of biological and physiological systems, Comput. Methods Programs Biomed., 88 (2007), 52–61. https://doi.org/10.1016/j.cmpb.2007.07.002 doi: 10.1016/j.cmpb.2007.07.002
|
| [62] |
M. Gultekin, S. Dundar, B. Keskinkilic, M. Turkyilmaz, N. Ozgul, K. Yuce, et al., How to triage HPV positive cases: Results of four million females, Gynecol. Oncol., 158 (2020), 105–111. https://doi.org/10.1016/j.ygyno.2020.04.698 doi: 10.1016/j.ygyno.2020.04.698
|
| [63] |
H. Miao, X. Xia, A. S. Perelson, H. Wu, On identifiability of nonlinear ODE models and applications in viral dynamics, SIAM Rev., 53 (2011), 3–39. https://doi.org/10.1137/090757009 doi: 10.1137/090757009
|
| [64] |
O. R. Isik, N. Tuncer, M. Marctheva, Mathematical model of measles in Turkey, J. Biol. Syst., 32 (2024), 941–970. https://doi.org/10.1142/S0218339024500323 doi: 10.1142/S0218339024500323
|
| [65] |
N. Tuncer, M. Marctheva, B. LaBarre, S. Payoute, Structural and practical identifiability analysis of zika epidemiological models, Bull. Math. Biol., 80 (2018), 2209–2241. https://doi.org/10.1007/s11538-018-0453-z doi: 10.1007/s11538-018-0453-z
|
| [66] | K. Diethelm, A. D. Freed, The fracPECE subroutine for the numerical solution of differential equations of fractional order, in Orschung und Wissenschaftliches Rechnen: Beiträge zum Heinz-Billing-Preis 1998, (1999), 57–71. |
| [67] |
R. Garrappa, Numerical solution of fractional differential equations: A survey and a software tutorial. Mathematics, 6 (2018), 16. https://doi.org/10.3390/math6020016 doi: 10.3390/math6020016
|
| [68] | A. Zinihi, Identifying memory effects in epidemics via a fractional SEIRD model and physics-informed neural networks, preprint, arXiv: 2509.22760. https://doi.org/10.48550/arXiv.2509.22760 |
Figures(11) / Tables(14)
Veysel Fuat Hatipoğlu, Maia Martcheva. Fractional gender structured model of human papillomavirus (HPV)[J]. Mathematical Biosciences and Engineering, 2026, 23(6): 1708-1742. doi: 10.3934/mbe.2026062
DownLoad: