Revisiting HPV infection pattern among urban Indonesian women in general population and its implication on health burden: A cross-sectional analysis from Indonesian Noncommunicable Disease Research 2016

Abstract

Objective: To identify circulating HPV types among urban Indonesian women and their specific co-infection patterns in bid to curb HPV infection in the general population and minimize its complications.

Methods: Urban Indonesian women from general population were selected as sample framework. Sample size and distribution across regions were determined by the Indonesian Bureau of Statistics (Badan Pusat Statistik, BPS), which represented the national level. Up to 35 408 cervical swab specimens were collected from August to September 2016 in 34 Indonesian provinces, categorized into six regions based on the development criteria set by the Ministry of National Development Planning (Badan Perencanaan Pembangunan Nasional, BAPPENAS). From all 1874 samples identified as HPV-positive, hybrid capture was implemented to evaluate type-specific HPV. This study analyzed descriptive data to determine the core-cluster of HPV combination. Co-occurrence HPV network was assessed using ‘qgraph’ package version 1.6.3 and computed in R version 3.6.3. Two-HPV association was analyzed in logistic regression using bias-reduction generalized linear model (brglm2) package version 0.5.1 adjusted by age and six main Indonesian regions.

Results: The logistic regression analysis demonstrated that HPV type 52 had rare relationship despite its common co-occurrence, cementing its role in single HPV infection. HPV type 16 and 18 tended to form infection cluster and were strongly associated with other types.

Conclusions: HPV type 52 was the most frequent HPV type among urban Indonesian women and accounted for most single infection cases. Concurrently, HPV 16 and HPV 18 accounted for most multiple infection cases and had strong tendency to attract other types, which may add further complications. However, due to lack of cytology and histological examination and information for other potential determinants, further in-depth studies are necessary to confirm whether these infection patterns truly connect to certain clinical outcomes.

Keywords

Human papilloma virus / HPV / Infection pattern / Urban women / Indonesia

References

[1] Arbyn M, Weiderpass E, de Sanjosé S, Saraiya M, Ferlay J, Bray F. Estimates of incidence and mortality of cervical cancer in 2018: A worldwide analysis. Lancet Global Health 2019; 8(2): e191-e203.

[2] Torres-Poveda K, Ruiz-Fraga I, Madrid-Marina V, Chavez M, Richardson V. High risk HPV infection prevalence and associated cofactors: A population-based study in female ISSSTE beneficiaries attending the HPV screening and early detection of cervical cancer program. BMC Cancer 2019; 19(1): 1-12.

[3] Handayani S, Dewi RM, Wibowo HA, Dany F, Kipuw NL, Adelina R, et al. Cervical cancer risk factor: HPV infection among Indonesian urban women Proceedings of the 4th International Symposium on Health Research.. USA: Atlantis Press; 2019, p. 452-455.

[4] Schmitt M, Depuydt C, Benoy I, Bogers J, Antoine J, Arbyn M, et al. Multiple human papillomavirus infections with high viral loads are associated with cervical lesions but do not differentiate grades of cervical abnormalities. J Clin Microbiol 2013; 51(5): 1458-1464.

[5] Aguilar-lemarroy A, Vallejo-Ruiz V, Cortés-Gutiérrez EI, Salgado-Bernabé ME, Ramos-González NP, Ortega-Cervantes L, et al. Human papillomavirus infections in Mexican women with normal cytology, precancerous lesions, and cervical cancer: Type-specific prevalence and HPV coinfections. J Med Virol 2015; 87(5): 871-884.

[6] Xu H, Lin A, Chen Y, Dong S, Shi W. Prevalence characteristics of cervical human papillomavirus (HPV) genotypes in the Taizhou area, China: A cross-sectional study of 37967 women from the general population. Br Med J Open 2017; 7(6): 1-9.

[7] Vu LTH. High-risk and multiple human papillomavirus infections among married women in Can Tho, Viet Nam. WPSAR 2012; 3(3): 57-62.

[8] Tran LTH, Tran LT, Bui TC, Le DT, Nyitray AG, Markham CM, et al. Risk factors for high-risk and multi-type human papillomavirus infections among women in Ho Chi Minh City, Vietnam: A cross-sectional study. BMC Women's Health 2015; 15(1): 1-9.

[9] Sahiratmadja E, Tobing MDL, Dewayani BM, Hernowo BS, Susanto H. Multiple human papilloma virus infections predominant in squamous cell cervical carcinoma in Bandung. Universa Med 2014; 33(1): 58-64.

[10] Bao H, Chen W, Zhang X, Bi H, Zhao Y, Fang L, et al. Prevalence of high-risk human papillomavirus in cervical intraepithelial neoplasia in the pre-vaccine era-China, 2017-2018. China CDC Weekly 2022; 4(48): 1083-1087.

[11] Kristina SA, Permitasari NPAL. Knowledge, attitudes and barriers towards human papillomavirus (HPV) vaccination in developing economies countries of South-East Asia region: A systematic review. Syst Rev Pharm 2019; 10(1): 81-86.

[12] Gallegos-Bolaños J, Rivera-Domínguez JA, Presno-Bernal JM, Cervantes-Villagrana RD. High prevalence of co-infection between human papillomavirus (HPV) 51 and 52 in Mexican population. BMC Cancer 2017; 17(531): 1-8.

[13] Li Z, Liu F, Cheng S, Shi L, Yan Z, Yang J, et al. Prevalence of HPV infection among 28457 Chinese women in Yunnan Province, southwest China. Sci Rep 2016; 6(1): 1-8.

[14] Tan SC, Ismail MP, Duski DR, Othman NH, Ankathil R. Prevalence and type distribution of human papillomavirus (HPV) in Malaysian women with and without cervical cancer: An updated estimate. Biosci Rep 2018; 38(2): 1-10.

[15] Murdiyarso LS, Kartawinata M, Jenie I, Widjajahakim G, Hidajat H, Sembiring R, et al. Single and multiple high-risk and low-risk human papillomavirus association with cervical lesions of 11224 women in Jakarta. Cancer Causes & Control 2016; 27: 1371-1379.

[16] Firman T, Guilmoto CZ, Jones GW, eds . Demographic patterns of Indonesia's urbanization, 2000-2010: Continuity and change at the macro level. In: Demographic transformation and socio-economic development. UK: Springer; 2016, p. 255-270.

[17] Giannella L, Delli Carpini G, Di Giuseppe J, Bogani G, Gardella B, Monti E, et al. Trend of HPV 16/18 genotypes in cervical intraepithelial neoplasia grade 3: Data for 2007-2018. Infect Drug Resist 2021; 14: 3763-3771.

[18] Salazar KL, Zhou HS, Xu J, Peterson LE, Schwartz MR, Mody DR, et al. Multiple human papilloma virus infections and their impact on the development of high-risk cervical lesions. Acta Cytol 2015; 59(5): 391-398.

[19] Hajia M, Sohrabi A. Possible synergistic interactions among multiple HPV genotypes in women suffering from genital neoplasia. Asian Pacific J Cancer Prev 2018; 19(3): 785-789.

[20] De Brot L, Pellegrini B, Moretti ST, Carraro DM, Soares FA, Rocha RM, et al. Infections with multiple high-risk HPV types are associated with high-grade and persistent low-grade intraepithelial lesions of the cervix. Cancer Cytopathol 2017; 125(2): 138-143.

[21] Smith MA, Tellier PP, Roger M, Coutlée F, Franco EL, Richardson H. Determinants of human papillomavirus coinfections among montreal university students: The influence of behavioral and biologic factors. Cancer Epidemiol Biomarkers Prev 2014; 23(5): 812-822.

[22] Biryukov J, Meyers C. Superinfection exclusion between two high-risk human papillomavirus types during a coinfection. J Virol 2018; 92(8): 1-20.

[23] Sakamoto J, Kamiura S, Okayama K, Okodo M, Shibata T, Osaka Y, et al. Single type infection of human papillomavirus as a cause for high-grade cervical intraepithelial neoplasia and invasive cancer in Japan. Papillomavirus Res 2018; 6: 46-51.

[24] Bee KJ, Gradissimo A, Chen Z, Harari A, Schiffman M, Raine-Bennett T, et al. Genetic and epigenetic variations of HPV 52 in cervical precancer. Int J Mol Sci 2021; 22(12): 1-15.

[25] Senapati R, Nayak B, Kar SK, Dwibedi B. HPV genotypes co-infections associated with cervical carcinoma: Special focus on phylogenetically related and non-vaccine targeted genotypes. PLoS One 2017; 12(11): 1-10.

[26] So KA, Lee IH, Lee KH, Hong SR, Kim YJ, Seo HH, et al. Human papillomavirus genotype-specific risk in cervical carcinogenesis. J Gynecol Oncol 2019; 30(4): 1-9.

[27] Kosen S, Andrijono A, Ocviyanti D, Indriatmi W. The cost-effectiveness of quadrivalent human papillomavirus vaccination in Indonesia. Asian Pacific J Cancer Prev 2017; 18(7): 2011-2017.

[28] Mariani L, Preti M, Cristoforoni P, Stigliano CM, Perino A. Overview of the benefits and potential issues of the nonavalent HPV vaccine. Int J Gynecol Obstet 2017; 136(3): 258-265.

[29] Frati ER, Bianchi S, Amendola A, Colzani D, Petrelli F, Zehender G, Tanzi E. Genetic characterization of variants of HPV-16, HPV-18 and HPV-52 circulating in Italy among general and high-risk populations. Mol Med Rep 2020; 21(2): 894-902.

[30] Wang X, Cheng Z. Cross-sectional studies: Strengths, weaknesses, and recommendations. Chest 2020; 158(1): S65-S71.

[31] Dickson EL, Vogel RI, Bliss RL, Downs LS Jr. Multiple-type human papillomavirus (HPV) infections: A cross-sectional analysis of the prevalence of specific types in 309 000 women referred for HPV testing at the time of cervical cytology. Int J Gynecol Cancer 2013; 23(7): 1295-1302.

[32] Rositch AF, Koshiol J, Hudgens MG, Razzaghi H, Backes DM, Pimenta JM, et al. Patterns of persistent genital human papillomavirus infection among women worldwide: A literature review and meta-analysis. Int J Cancer 2013; 133(6): 1271-1285.

[33] Hoffman SR, Le T, Lockhart A, Sanusi A, Dal Santo L, Davis M, et al. Patterns of persistent HPV infection after treatment for cervical intraepithelial neoplasia (CIN): A systematic review. Int J Cancer 2017; 141(1): 8-23.