- Research article
- Open Open Peer Review
Human papillomavirus infection in Rwanda at the moment of implementation of a national HPV vaccination programme
BMC Infectious Diseasesvolume 16, Article number: 225 (2016)
Cervical cancer is the most common female cancer in Rwanda that, in 2011, became the first African country to implement a national vaccination programme against human papillomavirus (HPV).
To provide a robust baseline for future evaluations of vaccine effectiveness, cervical cell specimens were obtained from 2508 women aged 18–69 years from the general population in Kigali, Rwanda, during 2013/14. 20 % of women were HIV-positive. Samples were used for liquid-based cytology and HPV testing (44 types) with GP5+/6+ PCR.
HPV prevalence was 34 %, being highest (54 %) in women ≤19 years and decreasing to 20 % at age ≥50. Prevalence of high risk (HR) HPV and cytological abnormalities was 22 and 11 % respectively (including 2 % with high-grade squamous intraepithelial lesions, HSIL) decreasing with age. Age-standardised prevalence of HR HPV was 22 % (or 19 % among HIV-negative women), and HPV16 was the most common type. Prevalence of HPV and cytological abnormalities were significantly higher in HIV-positive than HIV-negative women, and the difference increased with age. Other significant risk factors for HPV positivity in multivariate analyses were high lifetime number of sexual partners, receiving cash for sex, and being a farmer. 40 % of women with HSIL were infected with HPV16/18 and there was no significant difference between HIV-positive and HIV-negative women.
This study confirms Rwanda to be a setting of high prevalence of HPV and cervical disease that is worsened by HIV. These data will serve as a robust baseline for future evaluations of HPV vaccine programme effectiveness.
Cervical cancer represents by far the most common cancer among females in Rwanda, and has a high age-standardised incidence rate of approximately 42 cases per 100,000 women per year  that is typical of many sub-Saharan African countries. In response to this high burden, Rwanda became the first African country to initiate a national vaccination programme against human papillomavirus (HPV), the necessary cause of cervical cancer. In 2011, over 92,000 girls in primary school grade six (approximately 12 years old) were vaccinated with quadrivalent (HPV16/18/6/11) HPV vaccine. The three-dose vaccination coverage was estimated at 93 % in the target population . During 2012 and 2013, a catch-up vaccination targeted girls also in secondary school grade three (approximately 15 years old) whereas from 2014, vaccination has reverted to targeting primary school grade six only.
However, the impact of HPV vaccination on cervical cancer will likely not be observed for 20 years at least, so earlier evidence of vaccine programme effectiveness is important to encourage planners in Rwanda and other African countries to sustain their programmes. To this end, the Rwandan Ministry of Health, in collaboration with the International Agency for Research on Cancer (IARC), initiated a series of studies to evaluate the early impact of HPV vaccination on HPV prevalence.
A pre-vaccination study of the prevalence of, and risk factors for, HPV infection, was initiated in 2011 according to the standard protocol of the IARC HPV prevalence surveys . Cervicovaginal samples were collected from women aged 18–69 years living in the capital, Kigali, with a particular effort to oversample women ≤25 years of age in whom the impact of vaccination on HPV prevalence should be seen the earliest. The influence of HIV positivity, which is known to be an important determinant of HPV prevalence and may affect the future effectiveness of HPV vaccination, was also evaluated.
Between July 2013 and May 2014, a survey was conducted by the Ministry of Health of Rwanda in collaboration with the IARC, Lyon, France. The study aim was to enrol 2500 women from the general population using an age-stratified approach, namely 1000 women aged below 25 years old, 500 women aged 25 to 29 years old, 200 women in each 5-year age group between 30–34 and 45–49 years, and 200 women aged ≥50 years. All mentally and physically competent women were eligible for the study, regardless of their marital status, with exception of those who were known to be pregnant.
All participants were residents of Nyarugenge District, Kigali and were recruited in two ways: 1) A population-based invitation of women aged 18–69 years residing in pre-defined areas of Nyarugenge District was made through community workers and community meetings. Women were invited to attend the Reproductive Health Department of the district hospital (Muhima). Participation rates are difficult to estimate due to the uncertainty on the number of women who had been reached by the information campaign. 2) An additional opportunistic approach was used to recruit women aged ≤29 years who were spontaneously attending Muhima hospital or eight other health centers in the Nyarugenge District (Muhima, Butamwa, Rugarama, Mwendo, Corum, Gitega, Kabusunzu, Biryogo). Principal reasons for consultation were family planning, childhood vaccination, post-natal care services and HIV-related services. Study procedures were undertaken in the centre where the woman was consulting and very few refusals were recorded.
Following signature of an informed consent form, a structured risk questionnaire was administered to all study participants, and data transferred electronically to IARC via a specifically designed on-line webpage/application.
A cytobrush (Cervex-Brush, Rovers Medical Devices B.V., The Netherlands) was used for the collection of exfoliated cervical cells from the endocervix and ectocervix for all women. The brush containing cellular material was placed in a vial containing PreservCyt medium (Cytyc, Boxbourough, MA, USA) and stored at +4C until shipment.
Women aged ≥30 years old also underwent screening by visual inspection with acetic acid (VIA). Women with visible abnormalities of the cervix, and/or who were positive for high risk (HR) HPV infection (for which results became available at a later date, see below) were to be recalled for a second VIA and treatment according to national protocols (including cryotherapy, thermo coagulation or loop electrosurgical excision of the cervix transition zone). Biopsies were to be taken from all HPV-positive women, either VIA-directed to a visible lesion or randomly from 12 o’clock of the transformation zone in the absence of a visible lesion, prior to any treatment. Histological confirmation of cervical tissue was performed at the Department of Pathology at the Kigali Teaching University Hospital. Histology results are the subject of ongoing quality control and, apart from invasive cervical cancer, are not reported here. Among women ≤29 years, no VIA was performed, nor were HR HPV-positive women recalled for follow-up.
HPV testing and liquid based cytology (LBC)
HPV testing was performed on exfoliated cervical cells in the Department of Pathology at the VU University Medical Center, Amsterdam. DNA was extracted from the PreservCyt sample using magnetic beads (Macherey-Nagel) on a robotic system (Hamilton Robotics) according to the manufacturer’s instructions. Beta-Globin PCR analysis was conducted first to confirm the presence of human DNA in all specimens . The presence of HPV DNA was determined by conducting a general primer GP5+/6+ − mediated PCR, which permits the detection of a broad spectrum of genital HPV types . HPV positivity was assessed by hybridization of PCR products in an enzyme immunoassay with two oligoprobe cocktails that, together, detect 44 mucosal HPV types. Subsequent HPV genotyping was conducted by reverse-line blot hybridization of GP5+/6+ PCR products . HPV types considered HR types for this analysis were HPV16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68 ; possible HR types HPV26, 30, 34, 53, 66, 67, 69, 70, 73, 82, and 85; low-risk (LR) types HPV6, 11, 32, 40, 42, 43, 44, 54, 55, 57, 61, 64, 71, 72, 81, 83, 84, 86, 89, and 90.
Liquid-based cytology slides were prepared using a Thin Prep 3000 processor (Cytyc-Hologic), and were read at the Department of Pathology at the Vrije University Medical Center, Amsterdam, The Netherlands. Cytological diagnoses were reported according to CISOE-A standards  and were translated into Bethesda 2001 terminology system as 1) normal, 2) atypical squamous cells of undetermined significance; atypical squamous cells cannot exclude high-grade lesion; atypical glandular cells of undetermined significance; low-grade squamous intraepithelial lesions (ASCUS/ASC-H/AGUS/LSIL), 3) high-grade squamous intraepithelial lesions (HSIL), or 4) invasive cancer.
Age-standardised HPV prevalence was computed using the world standardised population to allow comparisons of HPV prevalence with other IARC HPV surveys . Adjusted prevalence ratios (PR) for HPV positivity and corresponding 95 % confidence intervals (CI) were computed using two binomial regression models with a log link, the first adjusted for age group (<24; 25–34, 35–44 and ≥45 years) only, and a second model that adjusted additionally for lifetime sexual partners (1, 2 or ≥3) and HIV status (positive, negative, unknown). Risk trends were assessed by considering categories as continuous variables. Effect modification of HIV status on the association between HPV16/18 positivity and cytology was tested using the log-likelihood ratio test.
General female population
Of 2511 participants providing cervical cell samples, 2508 had valid HPV results and were included in the following analyses, including 1233 and 1275 recruited at Muhima Hospital and health centres, respectively. Ninety nine percent of women reported no HPV vaccination and 98 % no history of previous cervical cancer screening. Among the 2391 women with a valid result on LBC, 255 (11 %) had an abnormal cytological diagnosis, including 210 (9 %) with ASCUS/ASC-H/AGUS/LSIL and 45 (2 %) with HSIL.
Overall HPV prevalence was 34 %, being 29 % among women with normal cytology (Table 1). The prevalence of HPV was 73 and 96 % in ASCUS/ASC-H/AGUS/LSIL and HSIL, respectively. In total, 480 (19 %) women had single-type and 377 (15 %) had multiple-type infections. Crude and age-standardised prevalence of HR HPV types was 22 %. Prevalence of possible HR HPV and LR HPV was 10 and 15 % respectively. HPV16 was the most common HPV type in women with normal cytology (4 %), ASCUS/ASC-H/AGUS/LSIL (11 %), and HSIL (29 %). HPV52 and HPV35 were the next most common HR types in women with normal cytology and ASCUS/ASC-H/AGUS/LSIL, compared to HPV18 and HPV58 in HSIL. Four invasive cervical cancers were found in the study population (2 HPV16, 1 HPV33, 1 HPV negative), and two additional invasive cervical cancers among participants without a valid sample for HPV and cytology.
Figure 1 shows the age-specific prevalence of HPV and cytological abnormalities, with HPV-positive women classified hierarchically into (1) HPV16 and/or 18, (2) other HR types and (3) LR types only. Overall HPV prevalence was highest (54 %) among women ≤19 years, and decreased with age, down to 20 % among women aged ≥50 years. Prevalence of cytological abnormalities showed no clear trend by age group.
A total of 490 (20 %) study participants self-reported to be HIV-positive, and 1414 (56 %) to be HIV-negative. An additional 604 (24 %) who reported not to know their HIV-status are hereafter considered as HIV-negative, based on similar risk factor profiles and sensitivity analyses after the exclusion of women with unknown HIV status. Eighty-four percent of HIV-positive women were receiving combined anti-retroviral therapy (cART) (data not shown). Of note, HIV-positive women represented 11 % of the study population below age 25 but 34 % in those age 45 or older.
Table 2 shows the relationship between HPV positivity and selected women’s characteristics. In the age-adjusted model, variables significantly associated with HPV positivity were age group, occupation, tobacco use, marital status, number of lifetime sexual partners, history of receiving cash for sex, the presence of anogenital warts, HIV positivity, number of pregnancies and contraceptive use in the last year. Additional adjustment for the two strongest risk factors, lifetime number of sexual partners and HIV, attenuated most PRs for HPV positivity, and the remaining significant associations included age group (PR for ≥45 vs ≤24 years = 0.45; 95 % CI: 0.37–0.55), marital status (PR for single versus married = 1.16; 95 % CI: 1.02–1.30), being a farmer (PR for farmers versus housewives = 0.78; 95 % CI: 0.66–0.92), number of lifetime sexual partners (PR for ≥3 vs 1 = 1.53; 95 % CI: 1.33–1.77), HIV positivity (PR = 1.54; 95 % CI: 1.36–1.75), number of pregnancies (PR for ≥4 vs 1 = 0.78; 95 % CI: 0.63–0.97), and use of contraceptive methods in the last year (PR for any = 0.82; 95 % CI: 0.74–0.91). Among contraceptive methods, condom did not reach statistical significance in multi-variate analysis. No significant associations with HPV positivity were seen for body mass index (5 % underweight, 7 % obese), age at menarche (median = 14 years), age at first sexual intercourse (median = 18 years), nor difference in age with partner at first sexual intercourse (median = 5 years younger) (data not shown).
In addition to total HPV prevalence, HR HPV prevalence was also significantly higher among HIV-positive (32 %) than in HIV-negative (20 %) women. Figure 2 shows prevalence of HR HPV infection, and of cytological abnormalities among HR HPV-positive women, stratified by age group and HIV positivity. PRs for HR HPV prevalence by HIV-positivity increased with increasing age from 1.6 (95 % CI: 1.3–2.0) in women below age 25 up to 4.0 (95 % CI: 2.3–7.0) in women aged 45 or older (Fig. 2). The percentage of HR HPV-positive that had cytological abnormalities was significantly higher in HIV-positive than HIV-negative women below age 25 (PR = 1.7, 95 % CI: 1.2–2.6) but similar in older age groups.
The proportion of women who were positive for HPV16/18 is shown in Table 3, by cytological diagnosis and HIV status. Overall, the proportion of HPV16/18-positive women increased from 5.5 % in normal cytology up to 40 % in HSIL with no statistically significant heterogeneity by HIV status (p for heterogeneity = 0.33) (Table 3).
Of 976 women that underwent VIA, only 5 % were classified as positive and only one of the 22 HSILs detected in this group of women was VIA-positive (data not shown).
We describe the epidemiology of HPV infection in the adult female population of Rwanda soon after the country successfully embarked upon a national programme of HPV vaccination. Data represent that of an almost entirely vaccine- and screening-naive population and disclose a high prevalence of HPV infection (34 %), particularly among women ≤19 years (54 %). Findings can be robustly compared with other IARC HPV surveys that were conducted according to a similar protocol. In this respect, the age-standardised HR HPV prevalence of 22 % in Rwanda (or 19 % in HIV-negative women) ranks below that of Guinea (31 %)  and Mongolia (25 %) , is similar to the prevalence detected in Nigeria (18 %)  and Vanuatu (19 %) , but remains more elevated than that found in many other areas known to be at high-risk for cervical cancer, including that in India (12 %)  and Colombia (10 %) . The high prevalence of HSIL and cervical cancer in the study population is consistent with recent estimates of high cervical cancer incidence in Rwanda .
HPV16 was by far the most frequently detected type in Rwanda, increasing from 4 % of normal cytology up to 29 % of HSIL. This agrees with a meta-analysis in which the proportion of HPV16-positive HSIL in Africa was 30 %, and lower than in other world regions . Vaccine types HPV16 and/or 18 accounting for 42 % of HSIL, and a relatively high importance of HPV35 and HPV58 in Rwanda are also consistent with previous African studies . The importance of HPV16 and 18 increased with the severity of cytological lesions similarly in HIV-negative and HIV-positive, and we did not observe any statistically significant under-representation of HPV16/18 in HIV-positive HSIL, as has been suggested by meta-analyses .
HIV infection was associated with higher prevalence of HPV, confirming previous findings from Rwanda [18–20]; sub-Saharan Africa [21–25] and elsewhere . This association persisted after adjustment for lifetime number of sexual partners, confirming an independent effect of HIV-infection beyond that of confounding by sexual behaviour. Indeed, lower CD4+ cell counts (a specific measure of HIV-related immunosuppression) are known to be negatively associated with the prevalence [26–29], persistence , and cumulative incidence  of HPV infection.
In agreement with our present findings, some studies in sub-Saharan African populations have reported declining HPV prevalence with age [21, 32–34], whereas others have showed HPV prevalence remaining equally high [10, 12, 35–37], or even increasing  in middle and old age. Even if it declined, HR HPV prevalence remained 10–15 % in women age 35 or older in Rwanda, suggesting that a large proportion of women might test positive for HPV-based screening and require diagnostic follow-up and treatment. This problem particularly applies to HIV-positive women, among whom HR HPV prevalence decreased less with age than among HIV-negative women.
Risk factors, other than HIV and age group, associated with cervical HPV infection were mostly proxies of sexual behaviour, as identified in previous HPV prevalence surveys. Lifetime number of sexual partners, being single, and history of receiving cash for sex (albeit rarely reported), were positively associated with HPV prevalence in multivariate analyses. Multiparity and recent use of contraceptives were also significantly and negatively associated with HPV positivity. Lower prevalence in farmers than other occupational groups suggest the possibility of an urban–rural gradient in the frequency of the infection. Residual confounding by sexual behaviour on the association with occupation, multiparity, and contraceptive use is difficult to rule out [39, 40].
Although only 2 % of study participants reported a prior Pap smear, the current study coincided with the initiation of screening and treatment programmes based on VIA, and/or testing for HPV DNA with careHPV (QIAGEN, Hilden, Germany) with use of VIA to triage HPV-positives for treatment . A potential warning about the performance of VIA in this setting comes from the finding that, despite very high HR HPV prevalence, VIA detected abnormalities in only 5 % of women in this high-risk population, and missed the vast majority of HSIL. This mirrors findings from other IARC HPV surveys [10, 42].
This study has strengths and weaknesses. Strengths include the relatively large number of women, especially young women, and the use of high-quality HPV testing and liquid-based cytology. HPV infection is asymptomatic and self-selection of women by level of HPV risk is unlikely to have occurred. The large number of HIV-positive women in our study will allow the opportunity to estimate HPV vaccine effectiveness in this high-risk group separately but also raises the question of whether our study population is representative of the general female population in Kigali, and Rwanda as a whole. HIV-positive women are over-represented in our survey in comparison to the general population . In particular, HIV-positive women (in the vast majority under cART treatment) were clearly more willing to accept screening invitations, or to be spontaneously seeking medical care, than HIV-negative women and represented approximately one third of study women age 35 or older. Had overall age-standardised HR HPV prevalence been estimated among HIV-negative women only, however, findings would have not been very different: 19 % instead of 22 %.
The findings of this study confirm Rwanda to be a setting of high prevalence of HPV and cervical disease that is worsened by HIV co-infection. These data will serve as a robust baseline for evaluations of HPV vaccine programme effectiveness in future cohorts of vaccinated girls.
Ethics approval and consent to participate
The present study had the approval of the Health Research Committee and Rwanda National Ethical Committee of the Rwandan Ministry of Health and the IARC Ethics Committee.
Informed consent was obtained from all study participants.
Consent for publication
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article.
atypical squamous cells of undetermined significance, atypical squamous cells cannot exclude high-grade lesion or atypical glandular cells of undetermined significance, low-grade squamous intraepithelial lesions
combined anti-retroviral therapy
high-grade squamous intraepithelial lesions
International Agency for Research on Cancer
visual inspection with acetic acid
Ferlay J, Soerjomataram I, Ervik M, Dikshit R, Eser S, Mathers C, Rebelo M, Parkin DM, Forman D, and Bray F. GLOBOCAN 2012 v1.0, Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11 [Internet]. Lyon, France: International Agency for Research on Cancer; http://globocan.iarc.fr. Accessed 20 June 15.
Binagwaho A, Wagner CM, Gatera M, Karema C, Nutt CT, Ngabo F. Achieving high coverage in Rwanda’s national human papillomavirus vaccination programme. Bull World Health Organ. 2012;90:623–8.
Clifford GM, Gallus S, Herrero R, Munoz N, Snijders PJ, Vaccarella S, Anh PT, Ferreccio C, Hieu NT, Matos E, Molano M, Rajkumar R, Ronco G, de Sanjosé S, Shin HR, Sukvirach S, Thomas JO, Tunsakul S, Meijer CJ, Franceschi S, and the IARC HPV Prevalence Surveys Study Group. Worldwide distribution of human papillomavirus types in cytologically normal women in the International Agency for Research on Cancer HPV prevalence surveys: a pooled analysis. Lancet. 2005;366:991–8.
Roda Husman AM, Snijders PJ, Stel HV, van den Brule AJ, Meijer CJ, Walboomers JM. Processing of long-stored archival cervical smears for human papillomavirus detection by the polymerase chain reaction. Br J Cancer. 1995;72:412–7.
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–7.
van den Brule AJ, Pol R, Fransen-Daalmeijer N, Schouls LM, Meijer CJ, Snijders PJ. GP5+/6+ PCR followed by reverse line blot analysis enables rapid and high-throughput identification of human papillomavirus genotypes. J Clin Microbiol. 2002;40:779–87.
IARC. Human papillomaviruses. IARC Monogr Eval Carcinog Risks Hum. 2007;90:1–636. http://monographs.iarc.fr/ENG/Monographs/vol90/index.php. Accessed 20 June 15.
Bulk S, van Kemenade FJ, Rozendaal L, Meijer CJ. The Dutch CISOE-A framework for cytology reporting increases efficacy of screening upon standardisation since 1996. J Clin Pathol. 2004;57:388–93.
Crosbie EJ, Einstein MH, Franceschi S, Kitchener HC. Human papillomavirus and cervical cancer. Lancet. 2013;382:889–99.
Keita N, Clifford GM, Koulibaly M, Douno K, Kabba I, Haba M, Sylla BS, van Kemenade FJ, Snijders PJ, Meijer CJ, Franceschi S. HPV infection in women with and without cervical cancer in Conakry, Guinea. Br J Cancer. 2009;101:202–8.
Dondog B, Clifford GM, Vaccarella S, Waterboer T, Unurjargal D, Avirmed D, Enkhtuya S, Kommoss F, Wentzensen N, Snijders PJ, Meijer CJ, Franceschi S, Pawlita M. Human papillomavirus infection in Ulaanbaatar, Mongolia: a population-based study. Cancer Epidemiol Biomarkers Prev. 2008;17:1731–8.
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–45.
Aruhuri B, Tarivonda L, Tenet V, Sinha R, Snijders PJ, Clifford G, Pang J, McAdam M, Meijer CJ, Frazer IH, Franceschi S. Prevalence of cervical Human Papillomavirus (HPV) infection in Vanuatu. Cancer Prev Res (Phila). 2012;5:746–53.
Franceschi S, Rajkumar R, Snijders PJF, Arslan A, Mahé C, Plummer M, Sankaranarayanan R, Cherian J, Meijer CJLM, Weiderpass E. Papillomavirus infection in rural women in southern India. Br J Cancer. 2005;92:601–6.
Molano M, Posso H, Weiderpass E, van den Brule AJ, Ronderos M, Franceschi S, Meijer CJ, Arslan A, Munoz N, and the HPV Study Group. Prevalence and determinants of HPV infection among Colombian women with normal cytology. Br J Cancer. 2002;87:324–33.
Guan P, Howell-Jones R, Li N, Bruni L, de Sanjosé S, Franceschi S, Clifford GM. Human papillomavirus types in 115,789 HPV-positive women: a meta-analysis from cervical infection to cancer. Int J Cancer. 2012;131:2349–59.
Clifford GM, Goncalves MA, Franceschi S, for the HPV and HIV Study Group. Human papillomavirus types among women infected with HIV: a meta-analysis. AIDS. 2006;20:2337–44.
Sinayobye J, Sklar M, Hoover DR, Shi Q, Dusingize JC, Cohen M, Mutimura E, Asiimwe-Kateera B, Castle PE, Strickler H, Anastos K. Prevalence and risk factors for High-Risk Human Papillomavirus (hrHPV) infection among HIV-infected and Uninfected Rwandan women: implications for hrHPV-based screening in Rwanda. Infect Agent Cancer. 2014;9:40.
Singh DK, Anastos K, Hoover DR, Burk RD, Shi Q, Ngendahayo L, Mutimura E, Cajigas A, Bigirimani V, Cai X, Rwamwejo J, Vuolo M, Cohen M, Castle PE. Human papillomavirus infection and cervical cytology in HIV-infected and HIV-uninfected Rwandan women. J Infect Dis. 2009;199:1851–61.
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.
Womack SD, Chirenje ZM, Gaffikin L, Blumenthal PD, McGrath JA, Chipato T, Ngwalle S, Munjoma M, Shah KV. HPV-based cervical cancer screening in a population at high risk for HIV infection. Int J Cancer. 2000;85:206–10.
Safaeian M, Kiddugavu M, Gravitt PE, Gange SJ, Ssekasanvu J, Murokora D, Sklar M, Serwadda D, Wawer MJ, Shah KV, Gray R. Determinants of incidence and clearance of high-risk human papillomavirus infections in rural Rakai, Uganda. Cancer Epidemiol Biomarkers Prev. 2008;17:1300–7.
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.
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.
Dartell M, Rasch V, Munk C, Kahesa C, Mwaiselage J, Iftner T, Kjaer SK. Risk factors for high-risk human papillomavirus detection among HIV-negative and HIV-positive women from Tanzania. Sex Transm Dis. 2013;40:737–43.
Chaturvedi AK, Dumestre J, Gaffga AM, Mire KM, Clark RA, Braly PS, Dunlap K, Beckel TE, Hammons AF, Kissinger PJ, Hagensee ME. Prevalence of human papillomavirus genotypes in women from three clinical settings. J Med Virol. 2005;75:105–13.
Strickler HD, Palefsky JM, Shah KV, Anastos K, Klein RS, Minkoff H, Duerr A, Massad LS, Celentano DD, Hall C, Fazzari M, Cu-Uvin S, Bacon M, Schuman P, Levine AM, Durante AJ, Gange S, Melnick S, Burk RD. Human papillomavirus type 16 and immune status in human immunodeficiency virus-seropositive women. J Natl Cancer Inst. 2003;95:1062–71.
Strickler HD, Burk RD, Fazzari M, Anastos K, Minkoff H, Massad LS, Hall C, Bacon M, Levine AM, Watts DH, Silverberg MJ, Xue X, Schlecht NF, Melnick S, Palefsky JM. Natural history and possible reactivation of human papillomavirus in human immunodeficiency virus-positive women. J Natl Cancer Inst. 2005;97:577–86.
Harris TG, Burk RD, Palefsky JM, Massad LS, Bang JY, Anastos K, Minkoff H, Hall CB, Bacon MC, Levine AM, Watts DH, Silverberg MJ, Xue X, Melnick SL, Strickler HD. Incidence of cervical squamous intraepithelial lesions associated with HIV serostatus, CD4 cell counts, and human papillomavirus test results. JAMA. 2005;293:1471–6.
Rowhani-Rahbar A, Hawes SE, Sow PS, Toure P, Feng Q, Dem A, Dembele B, Critchlow CW, N'Doye I, Kiviat NB. The impact of HIV status and type on the clearance of human papillomavirus infection among Senegalese women. J Infect Dis. 2007;196:887–94.
Massad LS, Xie X, Burk R, Keller MJ, Minkoff H, D’Souza G, Watts DH, Palefsky J, Young M, Levine AM, Cohen M, Strickler HD. Long-term cumulative detection of human papillomavirus among HIV seropositive women. AIDS. 2014;28:2601–8.
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–9.
Castellsagué 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–30.
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–42.
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–7.
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–25.
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–76.
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–9.
Vaccarella S, Herrero R, Dai M, Snijders PJ, Meijer CJ, Thomas JO, Hoang Anh PT, Ferreccio C, Matos E, Posso H, de SS, Shin HR, Sukvirach S, Lazcano-Ponce E, Ronco G, Rajkumar R, Qiao YL, Munoz N, Franceschi S. Reproductive factors, oral contraceptive use, and human papillomavirus infection: pooled analysis of the IARC HPV prevalence surveys. Cancer Epidemiol Biomarkers Prev. 2006;15:2148–53.
Vaccarella S, Franceschi S, Herrero R, Munoz N, Snijders PJF, Clifford GM, Smith JS, Lazcano-Ponce E, Sukvirach S, Shin HR, de Sanjose S, Molano M, Matos E, Ferreccio C, Anh PTH, Thomas JO, Meijer CJLM, and the IARC HPV Prevalence Surveys Study Group. Sexual behavior, condom use and HPV: pooled analysis of the International Agency for Research on Cancer HPV Prevalence Surveys. Cancer Epidemiol Biomarkers Prev. 2006;15:326–33.
Binagwaho A, Ngabo F, Wagner CM, Mugeni C, Gatera M, Nutt CT, Nsanzimana S. Integration of comprehensive women’s health programmes into health systems: cervical cancer prevention, care and control in Rwanda. Bull World Health Organ. 2013;91:697–703.
Li N, Shi JF, Franceschi S, Zhang WH, Dai M, Liu B, Zhang YZ, Li LK, Wu RF, De VH, Plummer M, Qiao YL, Clifford G. Different cervical cancer screening approaches in a Chinese multicentre study. Br J Cancer. 2009;100:532–7.
Rwanda Biomedical Center. Rwanda Global AIDS Response Progress Report 2014. Rwanda, Kigali: Republic of Rwanda, Ministry of Health; http://www.unaids.org/sites/default/files/country/documents/RWA_narrative_report_2014.pdf. Accessed 20 May 15
We would like to thank Marcelline Mahirwe and Evodie Mudaheranwa, the staff of Muhima Hospital, the staff in the health centres of Muhima, Butamwa, Rugarama, Mwendo, Corum, Gitega, Kabusunzu and Biryogo for their contribution to the field work in Rwanda. We also acknowledge the help of Dr Venerand Bigirimana, Pathologist at Kigali University Hospital for histological confirmation of biopsies, and Veronique Chabanis and Susan Gamon at IARC for administrative assistance.
The support for this project came from the Rwandan Ministry of Health, the International Agency for Research on Cancer and grants from the Bill & Melinda Gates Foundation, USA (grant numbers 35537 and OPP1053353). The Bill & Melinda Gates Foundation did not have any role in the writing process of this manuscript.
P.J.F. Snijders has Honoraria from Speakers Bureau from Roche, Gen-Probe, Qiagen, Abbott and Seegene and has a minority stake in Self-screen. No competing interests were disclosed by the other authors.
GC and SF conceived the study protocol. FN, MG and AG were involved in the development of the study methodology and infrastructure in Rwanda. MCU supervised the study protocol and data entry in Rwanda. IB, FL and VT were responsible for designing the electronic supports for data acquisition and/or statistical analyses. PS and AU were responsible for HPV testing and cytological reading, respectively. GC wrote the first draft of the manuscript. All authors were involved in the interpretation of the analyses and gave input to the final manuscript. All authors read and approved the final manuscript.