Age-specific prevalence of human papilloma virus infection among Nigerian women
BMC Public Health volume 14, Article number: 656 (2014)
Inconsistent trends in HPV prevalence by age have been described in Africa. We examined the age prevalence pattern and distribution of 37 HPV-DNA types among urban Nigerian women.
The study population was a sample of 278 women who presented to cervical cancer screening programs in Abuja, Nigeria, between April and August 2012. Using a nurse administered questionnaire, information on demographic characteristics and risk factors of cervical cancer was collected and samples of cervical exfoliated cells were obtained from all participants. Roche Linear Array HPV Genotyping Test® was used to characterize prevalent HPV and log-binomial regression models were used to examine the association between potential correlates and the prevalence of HPV infection.
The mean age (SD) of the women enrolled was 38 (8) years. The overall prevalence of HPV was 37%. HPV 35 was the most prevalent HPV type in the study population. Among women age ≤ 30 years, 52% had HPV infection compared to 23% of those women who were older than 45 years (p = 0.006). We observed a significant linear association between age and the prevalence of HPV infections. The prevalence ratio (PR) and 95% confidence interval (CI) was 2.26 (1.17, 4.34) for any HPV infection, 3.83 (1.23, 11.94) for Group 1 HPV (definite carcinogens), and 2.19 (0.99, 4.84) for Group 2a or 2b HPV (probable or possible carcinogens) types, among women aged 18–30 years, compared to women who were older than 45 years.
The prevalence of HPV infection was highest among younger women and decreased steadily with age among this population of urban Nigerian women.
Human papilloma viruses (HPV) are small DNA viruses, of which there are several types associated with benign and malignant conditions of the cervix, penis, vulva, vagina, anus and oropharynx [1, 2]. These viruses are ubiquitous and most women in the world are probably infected with at least one type of HPV during their sexual life giving a point prevalence of 10.1% . The Population Attributable Fraction of HPV associated cancer varies from 14.2% in sub-Saharan Africa and 15.5% in India to 1.2% in Australia/New Zealand and 1.6% in Northern America .
Previous studies have shown that the prevalence of HPV infection rises soon after the onset of sexual activity [4, 5] and attains the highest value in adolescence and early adulthood after which it declines during the second and third decades of life . While there seems to be a clear peak in HPV prevalence in younger women in Asia, Europe and Northern America, with a downward trend as age increases, age-specific HPV prevalence data from Central and South America appears to follow a slight U-shaped curve .
In Africa, the age-specific prevalence of HPV differs across countries. HPV prevalence was highest in younger women and decreased steadily with age in Kenya, Uganda and Zimbabwe [8–10] but generally reached a plateau at approximately 40 years of age in Nigeria and Mozambique [11, 12]. Whereas in Senegal, South Africa and another study from rural Nigeria, the peak prevalence in young women was accompanied by a second but smaller increase in prevalence among older women [13–16]. In Gambia, HPV prevalence remained relatively constant among women aged 15 to 54 years . The reason for these different age prevalence patterns in Africa remains unclear. Possible reasons include variations in sexual behavioral practices, sexual hygiene, co-morbidities, polygamy, male sexual partner behavior and male circumcision [15, 16, 18].
In this study, we examined the age prevalence pattern, distribution and socio-demographic correlates of HPV types among women in an urban, ethnically diverse population in Central Nigeria.
Between April and August 2012, we enrolled women from cervical cancer screening clinics at National Hospital, Abuja and University of Abuja Teaching Hospital, Nigeria. Abuja is the Federal Capital of Nigeria and it is located in the central part of the country. Its population is comprised of individuals from almost all ethnic groups in Nigeria, thus it is the most ethnically diverse city in the country . The study participants were 18 years or older, did not have a history of hysterectomy, not pregnant and had engaged in vaginal sexual intercourse. Interviewer administered questionnaires were used to collect data on socio-demographic characteristics, sexual and reproductive histories. Samples of exfoliated cervical cells were obtained from the cervical os with a cervical brush, which was stored in Thinprep® at -80°C until processing for HPV genotyping.
HPV detection by genotyping
HPV DNA was extracted using Gentra Puregene Cell Kit® according to the manufacturer’s instruction. HPV genotyping was done using Linear Array HPV genotyping test® (Roche Diagnostics), a qualitative in vitro test which utilizes amplification of target DNA by the Polymerase Chain Reaction and nucleic acid hybridization, and detects 37 HPV genotypes - (6, 11, 16, 18, 26, 31, 33, 35, 39, 40, 42, 45, 51, 52, 53, 54, 55, 56, 58, 59, 61, 62, 64, 66, 67, 68, 69, 70, 71, 72, 73 (MM9), 81, 82 (MM4), 83 (MM7), 84 (MM8), IS39, and CP6108. It has been validated and offers a reliable and sensitive approach for HPV DNA detection and genotyping in cervical specimens, using standardized quality-controlled reagents [20–23].
We classified HPV types according to the World Health Organization International Agency for Research on Cancer (IARC) Monographs Working Group assessment of the carcinogenicity of different HPV types . HPV types were classified as follows: definite carcinogens (group 1), HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58 and 59; Probable carcinogens (group 2a), HPV 68; and Possible carcinogens (group 2b), HPV 6, 11, 26, 40, 42, 53, 54, 55, 61, 62, 64, 66, 67, 69, 70, 71, 72, 73, 81, 82, 83, 84, IS39 and CP6108.
Descriptive analyses were performed to characterize the sampled population. Student t-tests were used to assess differences in the distribution of continuous variables between groups while χ2 and Fishers exact tests were used for categorical variables. We created age categories as follows: ≤ 30, 31 – 35, 36 – 40, 41 – 45 and > 45 years and calculated the Prevalence Ratio (PR) and 95% confidence intervals (CI) by age categories and by HPV group classification using univariate and multivariable log-binomial regression models. In separate models, we examined the association between HPV positivity and the following covariates: age, marital status, age at sexual initiation, total number of sexual partners, HIV status, contraceptive use, douching and body-mass index (BMI). To test for a linear trend, we assigned each age category the median value and modeled this variable as a continuous variable. All P-values were 2-sided and the analyses were performed using SAS 9.3 for UNIX statistical software.
The study was conducted according to the Nigerian National Code for Health Research Ethics and the Declaration of Helsinki. Ethical approval to conduct this study was obtained from the IHVN Health Research Ethics Committee and the University of Maryland, Baltimore, Institutional Review Board. Informed consent was obtained from all participants before enrollment in the study.
A total of 278 participants were enrolled in this study. We excluded 3 persons who had missing data on HPV genotype or several demographic variables, data from 275 participants were available for analysis. The mean age (SD) of the participants was 38 (8) and the median age (inter-quartile range) was 37 (31–42) years. Table 1 describes the demographic characteristics of the study participants.
The prevalence of any of the 37 HPV types in the study population was 37% (101/275). Group 1 HPV types (definite carcinogens) were found in 22% (61/275) of the study participants. HPV 35, a group 1 type, was the most prevalent type of HPV found in this study, it was found in 5% (14/275) of the women studied. Other group 1 HPV types found among the women studied were HPV 56, 58, 45 and 59 in decreasing order of prevalence. The prevalence of HPV16 and 18 in this population was low; HPV 16 was found in 2% (6/275) and HPV 18 was found in 2.5% (7/275) of the women studied. The prevalence of group 2a (Probable carcinogens) HPV type in the study participants (HPV68) was 3% (8/275). Group 2b HPV types (Possible carcinogens) were found in 24% (67/275) of the study participants. HPV81 and HPV42 were the most prevalent group 2b types, each was found in 4.7% (13/275) of the women studied.
Among women who had HPV infection, the age group distribution was 32% (33/101) for ≤30 years; 27% (27/101) for 31 – 35 years, 20% (20/101) for age group 36 – 40 years, 13% (13/101) for 41 – 45 years and 8% (8/101) for women older than 45 years (Table 1). The prevalence curves of HPV types for each age group are shown in Figure 1. The prevalence of any HPV infection was significantly reduced with increasing age. The multivariable PR of any HPV infection was 2.26 (1.17 to 4.34) for the lowest age category compared with the highest age category (referent group). This association remained significant for definite carcinogenic and probable/possible carcinogenic HPV types.
We found multiple HPV infections in 12.7% (35/275) of the women. Multiple HPV infections were more prevalent among younger women. Compared to women in the highest age group (>45 years), the PR (95% CI) of multiple HPV infection was 4.38 (1.07, 17.94), 5.00 (1.23, 20.30), 2.15 (0.47, 9.77) and 0.94 (0.16, 5.33), p for linear trend = 0.0004, for age group ≤ 30, 31 – 35, 36 – 40 and 41 – 45, respectively (Table 2).
The PR (95% CI) of any HPV infection was 2.54 (1.69, 3.81) p = < 0.001, comparing women who were HIV positive with those who were HIV negative. In multivariable analysis, we found a significant association between HIV and prevalent HPV infection, the PR (95% CI) of HIV for any HPV infection was 2.27 (1.48 to 3.48). HIV status was also associated with the prevalence of group 1 and group 2 HPV types in multivariable analysis, adjusting for age, marital status, age at sexual initiation, number of sexual partners, use of contraceptives, douching and body mass index. The multivariable adjusted PR (95% CI) was 4.01 (1.96, 8.21) for HIV and group 1 HPV, and 1.87 (1.11 to 3.17) for HIV and group 2 HPV (Table 3).
The multivariable PR (95% CI) of any HPV infection was 0.96 (0.93 to 1.00), comparing women with BMI ≥ 25 kg/m2 to those with BMI 18.5 - < 25 kg/m2. There was no significant association between BMI and group 1 or group 2 HPV infections. The multivariable PR (95% CI) was 0.97 (0.91, 1.02) for BMI and group 1; and 0.97 (0.92, 1.03) for BMI and group 2 HPV infections, comparing women with BMI ≥ 25 kg/m2 to those with BMI 18.5 - < 25 kg/m2.
In this cross-sectional study of urban Nigerian women, the prevalence of HPV infection was 37%. We observed a significant inverse linear association between age and the prevalence of HPV infections. This prevalence was approximately 4 times higher in women aged 18 – 30 years compared to women who were older than 45 years and younger women were more likely to have carcinogenic HPV types.
The prevalence of HPV infection among women in our study population was higher than had been reported in previous studies from Nigeria. Among 932 women surveyed between 1999 and 2000 in southwestern Nigeria, the prevalence of HPV was 26.3% . Our finding of HPV35 as the most prevalent type is similar to findings from southwest Nigeria and southern Mozambique, among women without invasive cervical cancer [12, 25].
Different trends have been reported on the age prevalence of HPV in Africa. The results of some studies showed there was a high peak of HPV prevalence among younger women which decreased with age, but among older women the prevalence either decreased [8, 9], reached a plateau [12, 26] or increased [13, 14]; whereas one study showed a relatively constant prevalence across all ages . Studies in rural and semi-rural southwestern Nigeria showed that the age prevalence of HPV was highest in younger women, declined with age and reached a plateau  or increased  among older women. Although the variation in age prevalence of HPV is well documented and appears to largely reflect differences in sexual behavior across geographical regions , the difference between our result and that of the previous studies from Nigeria may reflect characteristics of the populations surveyed. We studied an urban population of women in Abuja, central Nigeria who presented to hospitals for cervical cancer screening, while these other studies were community surveys in which participants were visited at their homes in southwestern Nigeria and asked to participate.
The results of some studies suggests that HIV infection is a known risk factor for HPV infection . In our study, we found an association between HIV and HPV infection and this association remained statistically significant when we examined the association between HIV and different categories of HPV. The results of our analyses on HIV and high-risk HPV infections among Nigerian women have been previously described . Although the prevalence of HPV infection appeared to be associated with some socio-demographic factors, these results did not reach significant levels in our multivariable analysis.
We found a marginal inverse association between any HPV infection and BMI in univariate analysis but this association was null in multivariable analysis. Prior studies of the association between body weight and HPV infections have yielded conflicting results. Among female sexual workers in Spain, there was an inverse association between high-risk HPV prevalence and BMI. In contrast, a study based on the US National Health and Nutrition Examination Survey (NHANES) showed that BMI was not associated with prevalence of HPV . The relationship between BMI and the risk of HPV infection needs further evaluation.
Our study has several important limitations. First, the sample size was small hence we were unable to examine the prevalence of HPV infection among women age 60 years and above. We are therefore unable to determine whether the prevalence of HPV in this study population would have continued to fall, remained stable or risen. Second, because of the cross sectional design of our study we could not examine temporal trends of potential correlates of HPV and the risk factors for persistent HPV infections. With continued follow up of this population and repeated sampling, we will be able to report the pattern and spectrum of persistent HPV infection in Nigeria in future.
The age prevalence pattern of HPV infection among women in cervical cancer screening programs in Nigeria is similar to that in many developed countries. Replication of this study in other urban African populations is required. Further studies to investigate epidemiological correlates of persistent high risk HPV infection are also needed. These studies will be useful in providing evidence for policies relating to the use of HPV DNA tests for screening, utility and coverage of current HPV vaccines and allocation of resources for HPV infections and cervical cancer prevention and treatment.
Body mass index
Cervical intraepithelial neoplasm
Human immunodeficiency virus
International Agency for Cancer Research
Institute of Human Virology Nigeria
National Health and Nutrition Examination Survey
Polymerase Chain Reaction REDCap Research Electronic Data Capture
Hausen H: Papillomavirus infections - a major cause of human cancers. Biochim Biophys Acta. 1996, 1288: F55-F78.
Forman D, de Martel C, Lacey CJ, Soerjomataram I, Lortet-Tieulent J, Bruni L, Vignat J, Ferlay J, Bray F, Plummer M, Franceschi S: Global burden of human papillomavirus and related diseases. Vaccine. 2012, 30 (Suppl 5): F12-F23.
Baseman JG, Koutsky LA: The epidemiology of human papillomavirus infections. J Clin Virol. 2005, 32 (Suppl 1): S16-S24.
Dunne EF, Unger ER, Sternberg M, McQuillan G, Swan DC, Patel SS, Markowitz LE: Prevalence of HPV infection among females in the United States. JAMA. 2007, 297: 813-819. 10.1001/jama.297.8.813.
Jacobs MV, Walboomers JM, Snijders PJ, Voorhorst FJ, Verheijen RH, Fransen-Daalmeijer N, Meijer CJ: Distribution of 37 mucosotropic HPV types in women with cytologically normal cervical smears: the age-related patterns for high-risk and low-risk types. Int J Cancer. 2000, 87: 221-227. 10.1002/1097-0215(20000715)87:2<221::AID-IJC11>3.0.CO;2-2.
Castellsague X, Schneider A, Kaufmann AM, Bosch FX: HPV vaccination against cervical cancer in women above 25 years of age: key considerations and current perspectives. Gynecol Oncol. 2009, 115: S15-S23. 10.1016/j.ygyno.2009.09.021.
Smith JS, Melendy A, Rana RK, Pimenta JM: Age-specific prevalence of infection with human papillomavirus in females: a global review. J Adolesc Health. 2008, 43: S5-S25. S25 e21-41
Temmerman M, Tyndall MW, Kidula N, Claeys P, Muchiri L, Quint W: Risk factors for human papillomavirus and cervical precancerous lesions, and the role of concurrent HIV-1 infection. Int J Gynaecol Obstet. 1999, 65: 171-181. 10.1016/S0020-7292(99)00043-0.
Serwadda D, Wawer MJ, Shah KV, Sewankambo NK, Daniel R, Li C, Lorincz A, Meehan MP, Wabwire-Mangen F, Gray RH: Use of a hybrid capture assay of self-collected vaginal swabs in rural Uganda for detection of human papillomavirus. J Infect Dis. 1999, 180: 1316-1319. 10.1086/315026.
Womack SD, Chirenje ZM, Blumenthal PD, Gaffikin L, McGrath JA, Chipato T, Ngwalle E, Shah KV: Evaluation of a human papillomavirus assay in cervical screening in Zimbabwe. BJOG. 2000, 107: 33-38. 10.1111/j.1471-0528.2000.tb11576.x.
Thomas JO, Herrero R, Omigbodun AA, Ojemakinde K, Ajayi IO, Fawole A, Oladepo O, Smith JS, Arslan A, Munoz N, Snijders PJ, Meijer CJ, Franceschi S: Prevalence of papillomavirus infection in women in Ibadan, Nigeria: a population-based study. Br J Cancer. 2004, 90: 638-645. 10.1038/sj.bjc.6601515.
Castellsague X, Menendez C, Loscertales MP, Kornegay JR, dos Santos F, Gomez-Olive FX, Lloveras B, Abarca N, Vaz N, Barreto A, Bosch FX, Alonso P: Human papillomavirus genotypes in rural Mozambique. Lancet. 2001, 358: 1429-1430. 10.1016/S0140-6736(01)06523-0.
Xi LF, Toure P, Critchlow CW, Hawes SE, Dembele B, Sow PS, Kiviat NB: Prevalence of specific types of human papillomavirus and cervical squamous intraepithelial lesions in consecutive, previously unscreened, West-African women over 35 years of age. Int J Cancer. 2003, 103: 803-809. 10.1002/ijc.10876.
Kuhn L, Denny L, Pollack A, Lorincz A, Richart RM, Wright TC: Human papillomavirus DNA testing for cervical cancer screening in low-resource settings. J Natl Cancer Inst. 2000, 92: 818-825. 10.1093/jnci/92.10.818.
Clarke MA, Gage JC, Ajenifuja KO, Wentzensen NA, Adepiti AC, Wacholder S, Burk RD, Schiffman M: A population-based, cross-sectional study of age-specific risk factors for high risk human papillomavirus prevalence in rural Nigeria. Infect Agent Cancer. 2011, 6: 12-10.1186/1750-9378-6-12.
Gage JC, Ajenifuja KO, Wentzensen NA, Adepiti AC, Eklund C, Reilly M, Hutchinson M, Wacholder S, Harford J, Soliman AS, Burk RD, Schiffman M: The age-specific prevalence of human papillomavirus and risk of cytologic abnormalities in rural Nigeria: implications for screen-and-treat strategies. Int J Cancer. 2012, 130: 2111-2117. 10.1002/ijc.26211.
Wall SR, Scherf CF, Morison L, Hart KW, West B, Ekpo G, Fiander AN, Man S, Gelder CM, Walraven G, Borysiewicz LK: Cervical human papillomavirus infection and squamous intraepithelial lesions in rural Gambia, West Africa: viral sequence analysis and epidemiology. Br J Cancer. 2005, 93: 1068-1076. 10.1038/sj.bjc.6602736.
Veldhuijzen NJ, Braunstein SL, Vyankandondera J, Ingabire C, Ntirushwa J, Kestelyn E, Tuijn C, Wit FW, Umutoni A, Uwineza M, Crucitti T, van de Wijgert JH: The epidemiology of human papillomavirus infection in HIV-positive and HIV-negative high-risk women in Kigali, Rwanda. BMC Infect Dis. 2011, 11: 333-10.1186/1471-2334-11-333.
NPC: Nigeria Demographic and Health Survey. Book Nigeria Demographic and Health Survey. 2008, City
Castle PE, Sadorra M, Garcia F, Holladay EB, Kornegay J: Pilot study of a commercialized human papillomavirus (HPV) genotyping assay: comparison of HPV risk group to cytology and histology. J Clin Microbiol. 2006, 44: 3915-3917. 10.1128/JCM.01305-06.
Coutlee F, Rouleau D, Petignat P, Ghattas G, Kornegay JR, Schlag P, Boyle S, Hankins C, Vezina S, Cote P, Macleod J, Voyer H, Forest P, Walmsley S, Franco E: Enhanced detection and typing of human papillomavirus (HPV) DNA in anogenital samples with PGMY primers and the Linear array HPV genotyping test. J Clin Microbiol. 2006, 44: 1998-2006. 10.1128/JCM.00104-06.
Stevens MP, Garland SM, Tabrizi SN: Validation of an automated detection platform for use with the roche linear array human papillomavirus genotyping test. J Clin Microbiol. 2008, 46: 3813-3816. 10.1128/JCM.01169-08.
van Hamont D, van Ham MA, Bakkers JM, Massuger LF, Melchers WJ: Evaluation of the SPF10-INNO LiPA human papillomavirus (HPV) genotyping test and the roche linear array HPV genotyping test. J Clin Microbiol. 2006, 44: 3122-3129. 10.1128/JCM.00517-06.
Schiffman M, Clifford G, Buonaguro FM: Classification of weakly carcinogenic human papillomavirus types: addressing the limits of epidemiology at the borderline. Infect Agent Cancer. 2009, 4: 8-10.1186/1750-9378-4-8.
Okolo C, Franceschi S, Adewole I, Thomas JO, Follen M, Snijders PJ, Meijer CJ, Clifford GM: Human papillomavirus infection in women with and without cervical cancer in Ibadan, Nigeria. Infect Agent Cancer. 2010, 5: 24-10.1186/1750-9378-5-24.
De Vuyst H, Steyaert S, Van Renterghem L, Claeys P, Muchiri L, Sitati S, Vansteelandt S, Quint W, Kleter B, Van Marck E, Temmerman M: Distribution of human papillomavirus in a family planning population in nairobi, kenya. Sex Transm Dis. 2003, 30: 137-142. 10.1097/00007435-200302000-00009.
Hanisch RA, Sow PS, Toure M, Dem A, Dembele B, Toure P, Winer RL, Hughes JP, Gottlieb GS, Feng Q, Kiviat NB, Hawes SE: Influence of HIV-1 and/or HIV-2 infection and CD4 count on cervical HPV DNA detection in women from Senegal, West Africa. J Clin Virol. 2013, 58: 696-702. 10.1016/j.jcv.2013.10.012.
Akarolo-Anthony SN, Al-Mujtaba M, Famooto AO, Dareng EO, Olaniyan OB, Offiong R, Wheeler CM, Adebamowo CA: HIV associated high-risk HPV infection among Nigerian women. BMC Infect Dis. 2013, 13: 521-10.1186/1471-2334-13-521.
Wee CC, Huang A, Huskey KW, McCarthy EP: Obesity and the likelihood of sexual behavioral risk factors for HPV and cervical cancer. Obesity (Silver Spring). 2008, 16: 2552-2555. 10.1038/oby.2008.394.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/14/656/prepub
Acknowledgement and funding
This work was supported by the UM-Capacity Development for Research in AIDS Associated Malignancy Grant (NIH/NCI 1D43CA153792-01).
CMW has received through the University of New Mexico, funds from grants and cooperative agreements from the US National institutes of Health related to cervical screening, funds from GSK for HPV vaccine studies and reimbursements for travel related to publication activities and equipment and reagents from Roche Molecular Systems for HPV genotyping.
Other authors report no competing interest.
SNA analyzed the data and drafted the manuscript. AOF, EOD, MA and IL coordinated the study and managed data accrual. OO and RU, site investigators, contributed to the study design, implementation and provided revisions of the manuscript. CMW performed HPV genotyping, contributed to data interpretation and provided revisions of the manuscript. CAA conceived the study, obtained funds, contributed to the study design and provided critical revisions of the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Akarolo-Anthony, S.N., Famooto, A.O., Dareng, E.O. et al. Age-specific prevalence of human papilloma virus infection among Nigerian women. BMC Public Health 14, 656 (2014). https://doi.org/10.1186/1471-2458-14-656