Skip to main content
  • Research article
  • Open access
  • Published:

The epidemiology of human papillomavirus infection in HIV-positive and HIV-negative high-risk women in Kigali, Rwanda



The prevalence, incidence and persistence of human papillomavirus (HPV) types in sub-Saharan Africa are not well established. The objectives of the current study are to describe (predictors of) the epidemiology of HPV among high-risk women in Kigali, Rwanda.


HIV-negative, high-risk women were seen quarterly for one year, and once in Year 2. HIV serostatus, clinical, and behavioral information were assessed at each visit, HPV types at Month 6 and Year 2, and other sexually transmitted infections (STI) at selected visits. HPV prevalence was also assessed in HIV-positive, high-risk women.


Prevalence of any HPV was 47.0% in HIV-negative women (median age 25 years) compared to 72.2% in HIV-positive women (median age 27 years; OR 2.9, 95% CI 1.9-4.6). Among HIV-negative women, cumulative incidence of high-risk (HR)-HPV was 28.0% and persistence 32.0% after a mean period of 16.6 and 16.9 months, respectively. Prior Chlamydia trachomatis and Neisseria gonorrhoeae infection, concurrent low-risk (LR)-HPV infection and incident HSV-2 were associated with HR-HPV prevalence among HIV-negative women; prior C. trachomatis infection and co-infection with LR-HPV and HPV16-related HPV types with HR-HPV acquisition. HPV16-related types were the most prevalent and persistent.


High HPV prevalence, incidence and persistence were found among high-risk women in Kigali. HPV52 had the highest incidence; and, together with HPV33 and HPV58, were strongly associated with acquisition of other HR-HPV types in HIV-negative women.

Peer Review reports


Human papillomavirus (HPV) is the most common sexually transmitted infection (STIs) worldwide [1]. Prevalence is highest among sexually active women below 25-35 years, usually declining thereafter [2, 3]. Worldwide age-adjusted HPV prevalence in women with normal cytology is estimated at 12%, ranging from 5% in Northern America to 34% in East Africa% [3]. Age and high-risk sexual behavior (such as high life-time number of sex partners and recent acquisition of a new sex partner) are well-recognized risk factors for acquiring HPV. Moreover, higher HPV prevalence rates have been documented among women at high risk for other STIs and among HIV-positive women [4, 5]. Most infections are cleared within two years but a minority of infections becomes persistent [6]. Persistent infection with high-risk (HR), oncogenic HPV genotypes is causally associated with the development of cervical cancer, the most common malignancy among women in Sub-Saharan Africa [7].

An increasing number of articles report on the type-specific HPV prevalence among women in Sub-Saharan Africa, but incidence and persistence data remain scarce. HPV dynamics, and predictors of incidence and persistence, are best studied in prospectively collected data. This information will aid in the design of primary or secondary prevention programs.

As part of secondary data-analysis of an HIV incidence study, we evaluated HPV prevalence, incidence, persistence, and their predictors, among HIV-negative high-risk women in Kigali, Rwanda. In addition, HPV prevalence was assessed among HIV-positive high-risk women cross-sectionally.


Study population

A cross-sectional HIV prevalence survey was conducted among high-risk women in Kigali, Rwanda. Subsequently, women testing HIV-negative during the survey, and who were 18 years or older, at high risk for HIV (defined as having exchanged sex for money at least once in the last month and/or currently having sex with multiple partners plus having sex at least twice per week) and were not pregnant at enrolment or planning a pregnancy within the next year were invited to participate in a prospective cohort study to assess HIV incidence (2006-2009) [8].

The target sample size for the HIV-negative cohort was 400 and this was based on the primary HIV incidence objective. HPV prevalence, incidence and persistence were among secondary objectives of the study. Enrolled women were followed quarterly for one year, and were invited for one additional visit during year 2 ("Year 2 visit").

Women testing HIV-positive during the survey were excluded from cohort participation, but were invited for one additional study visit during year 2 ("Year 2 visit") during which HPV prevalence was assessed.

Women were recruited via community meetings in three districts in Kigali. Written informed consent was obtained, including parental or guardian consent for women who were 18-20 years old. The study was conducted at Projet Ubuzima, an international non-governmental organization that operates a research clinic and laboratory, in Kigali, Rwanda.

Study procedures for HIV-negative participants

HIV-negative participants were tested for HIV at every study visit, for Neisseria gonorrhoeae, C. trachomatis, Treponema pallidum (syphilis), Trichomonas vaginalis, and bacterial vaginosis (BV) at enrolment, Month 6, and Month 12, and for BV also in Year 2. HPV typing and cervical cytology were done at Month 6 and Year 2. HSV-2 serology testing was performed at enrolment, at Month 12 for women testing negative at enrolment, and in Year 2 for women testing negative at Month 12. Standardized questionnaires were administered at each study visit.

Study procedures for HIV-positive participants

HIV-positive women participating in the Year 2 visit followed the same study procedures as described for the Year 2 visit of the HIV-negative cohort participants.

Women with abnormal cytology results were referred for visual inspection with acetic acid (VIA) and colposcopy by a gynecologist, and were offered appropriate treatment. Women with curable STIs received treatment and were offered partner notification and treatment. HIV-positive women were referred to HIV care. As part of the study procedures women received HIV prevention, safer sex and family planning counseling at each visit.

Human papillomavirus typing

HPV genotyping was performed at the Institute of Tropical Medicine (ITM, Antwerp, Belgium) using the Linear Array HPV Genotyping Test (LA) (Roche Molecular Systems, USA) according to manufacturer's instructions. In brief, cervical cells collected in Liquid Based Cytology media PreservCyt were stored at -80°C until testing. An aliquot of 250 μl was used for DNA extraction by AmpliLute Liquid Media Extraction Kit (Roche Molecular Systems, USA). HPV genotype targets were amplified using standard polymerase chain reaction (PCR). LA applies PGMY09/11 L1 as consensus primer system and co-amplifies human β-globulin as a cellular, extraction and inhibition control. Reverse line blot hybridization and detection identifies 37 high- and low-risk HPV types (6, 11, 16, 18, 26, 31, 33, 35, 39, 40, 42, 45, 51, 53, 54, 55, 56, 58, 59, 61, 62, 64, 66 to 73, 81, 82 subtype IS39, 82, 83, 84, CP108 and a mixed probe for types 33/35/52/58). An in-house Real-Time PCR (RT-PCR) was established at the ITM to identify DNA of HPV genotype 52 in the mixed probe positive samples. Briefly, DNA was prepared using the NucliSens miniMAG (BioMιrieux, Boxtel, the Netherlands). DNA templates were amplified in a qualitative RT-PCR assay. HPV52-F/R and hydrolysis probe HPV52 labeled at the 5' with FAM and at the 3' with BHQ-1 (Eurogentec, Belgium) were used for amplification and detection [9]. Amplification was performed using the Absolute Blue QPCR mix kit (Thermo Scientific, UK), 500 nm of each primer and probe, and 10 μL of DNA template, in a final reaction volume of 30 μL. The RT-PCR was performed on the Rotor-Gene™ 6000 platform (Corbett Life Science, USA) using the following thermal cycling conditions: initial enzyme activation at 95°C for 15 min; followed by 60 cycles of 95°C for 15 sec, 58°C for 45 sec, and 65°C for 45 sec. Detection of the fluorescence signal was recorded using the setting for hydrolysis probes (470 nm, 510 nm) following the extension phase cycle by cycle.

Other laboratory testing

HIV testing was done at Projet Ubuzima according to the national testing guidelines, with two positive rapid tests indicating an HIV diagnosis (First Response, Premier Medical Corporation, Daman, India; and Uni-Gold HIV, Trinity Biotech Plc, Ireland), and a third rapid test (Capillus, Trinidity Biotech Plc, Ireland) as tiebreaker. Positive results were confirmed with Murex HIV Ag/Ab Combination ELISA (Abbott Laboratories, Germany). CD4 cytometry was performed at the National Reference Laboratory in Kigali.

InPouch culture (Biomed Diagnostics, USA) was used for T. vaginalis diagnosis. Syphilis was diagnosed using RPR (Spinreact, Spain) and TPHA (Spinreact, Spain). HerpeSelect-2 ELISA (Focus Technologies, USA) was used for HSV-2 antibody testing. Optical density ≥ 3.5 were considered positive for HSV-2, and OD ≤ 1.1 negative. All tests were performed at Projet Ubuzima. Dry swabs (BD Culturette™ EZ Collection), stored at -80°C were used for Amplicor NG/CT PCR (Roche Molecular Systems, USA) for N. gonorrhoeae and C. trachomatis testing at ITM. Positive results were confirmed with a strand displacement amplification assay (Probetec, Becton Dickinson, USA). External quality control for STI testing at Projet Ubuzima was coordinated by ITM. BV was diagnosed by Gram stain Nugent score. Nugent scoring was performed at Projet Ubuzima; at the Academic Medical Center of the University of Amsterdam, Department of Medical Microbiology; and at the Amsterdam Municipal Health Service, Public Health Laboratory. Conventional cytology was performed at Ghent University Hospital, Belgium, using the Bethesda system for classification [10].

Statistical analysis

Data were double-entered, and analyzed using STATA 9.2 (StatCorp, TX, USA). Baseline characteristics were estimated as percentages with 95% confidence intervals (CI) for binary and categorical variables; median and interquartile range (IQR) for non-normally distributed continuous variables or mean and standard deviations for normally distributed continuous variables.

HPV classification

In these analyses, HPV16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 68, 73, 82 and IS39 were considered HR-HPV types [1113]. All other types were considered LR-HPV. HPV16-related types included HPV16, 31, 33, 35, 52, and 58. 'Any HPV' included HR and/or LR-HPV.

HPV prevalence

HR-HPV prevalence was estimated from the first HPV results in both HIV-negative (Month 6) and HIV-positive women (Year 2). Multivariable logistic regression models for predictors of HR or LR-HPV detection in HIV-positive and HIV-negative women were fitted using manual backward modeling. Variables associated with the outcome with a P-value < 0.10 in univariable models were included in the initial multivariable logistic regression model. Variables were retained in the model based upon Wald statistics (P < 0.05) and/or likelihood ratio tests for the overall fit of the model (P < 0.1). Age was included as a covariate in all models. Pair-wise correlations between the variables in the final models were investigated, and correlated variables were assessed for interaction. The comparison group for HR-HPV analysis included women without HR-HPV, irrespective of LR-HPV detection. Concurrent LR-HPV was included as a covariate. Similarly, the comparison group in LR-HPV analysis included women without LR-HPV, irrespective of HR-HPV, but concurrent HR-HPV was included as a covariate. Age-adjusted odds ratios and 95% confidence intervals (C.I.) were calculated comparing HPV prevalence among HIV-positive (n = 349) and HIV-negative women (n = 126). Women with abnormal cytology results (HIV-negative women n = 16; HIV positive women n = 14) were excluded.

HPV incidence (HIV-negative women)

Type-specific HR-HPV incidence was defined as newly detected type-specific HR-HPV at the second HPV measurement (Year 2) among women free of the type of interest at the first HPV measurement (Month 6). The (cumulative) incidence was calculated over the mean interval between Month 6 and Year 2 visit. All incident type-specific HR-HPV cases were grouped together into one binary incident HR-HPV outcome variable, with women without incident HR-HPV as the referent group (regardless of presence of other prevalent HR-HPV types or prevalent or incident LR-HPV). The association between prevalent HR-HPV types and incidence of other HR-HPV types was assessed. First, univariable associations between each type-specific HR-HPV at Month 6 and subsequent acquisition of other HR-HPV types were evaluated. To increase statistical power, types were then grouped according to their phylogenetic group and within group according to strength of the association (0.1 > P < 0.1). These groups were included in the final incidence model. Sensitivity analyses were performed excluding women with any prevalent HR-HPV and/or LR-HPV at Month 6 from the reference group. These results are not shown because they were similar to those of the incidence model as described above. Of the 349 HIV-negative women with normal cytology results and HPV results at Month 6, 296 women also had HPV results at Year 2 visit and contributed to the incidence analyses. Methods for multivariable logistic regression modeling were identical to those for HPV prevalence.

HPV persistence (HIV-negative women)

HR-HPV persistence was defined as women who tested positive for the same HR-HPV type at both visits during which HPV typing was performed (Month 6 and Year 2). Persistence was calculated over the mean interval between these two visits. All persistent type-specific HR-HPV cases were grouped into one binary persistent HR-HPV outcome variable, irrespective of LR-HPV. The reference group included women who cleared all HR-HPV infections at the second visit that were present at the first visit, but who may however tested positive for new HPV types at the second HPV measurement. Of the 111 HR-HPV positive women with normal cytology at Month 6, 90 women also had HPV results at Year 2 visit and contributed to the persistence analyses. Multivariable logistic regression modeling was conducted as described for HPV prevalence above.

Ethical approvals were obtained from the National Ethics Committee and the National AIDS Control Commission (CNLS) in Rwanda; the Institutional Review Board of Columbia University Medical Center, New York, USA; and the Medical Ethics Committee of the University of Antwerp, Belgium.


Study population characteristics

Eight hundred women participated in the HIV prevalence survey between October 2006 and August 2007 and 192 women tested HIV-positive. Four hundred three (403) HIV-negative, non-pregnant women were enrolled in the cohort. The flow-diagram of the study is given in Figure 1. All women, except one, self-identified as female sex worker (FSW). Median age at survey was 25 (IQR 22-28) for HIV-negative women and 27 (IQR 23-32) for HIV-positive women. Among HIV-negative cohort participants, prevalence of N. gonorrhoeae was 10% (CI 7.3-13.6), C. trachomatis 5% (CI 2.7-7.2), T. vaginalis 17% (CI 13.4-21.2), syphilis 7% (CI 4.7-10.1), and BV 45% (CI 39.4-49.9).

Figure 1
figure 1

Flow diagram - HPV analysis data set (cytologically normal women). Footnote: Abbreviations: HIV = human immunodeficiency virus; HPV = human papillomavirus; LFU = lost to follow-up *in total there were 6 HIV seroconversion between M12 and Year 2, but one participant was already dropped from HPV data analysis set because of abnormal cytology results at M6. HPV genotyping and cervical cytology was performed among HIV-negative women at Month 6 and Year 2 and among HIV-positive women at Year 2 only. Women with abnormal cytology results are excluded from the HPV analysis dataset.

Type-specific HPV prevalence in HIV-negative and HIV-positive women

Among 349 HIV-negative women with no cervical abnormalities, the HPV prevalence at the first HPV measurement visit (Month 6) was 47.0% for any HPV, 31.8% for HR-HPV and 32.4% for LR-HPV. The five most prevalent HR-HPV types among HIV-negative women were HPV16 (5.4%), HPV45 (4.9%), and HPV52 (4.3%) HPV58 (4.0%) and HPV68 (3.7%) (Figure 2). Infection with multiple HR-HPV types was present in 11.7% of HR-HPV positive/HIV-negative women.

Figure 2
figure 2

Distribution of high-risk HPV types in HIV-negative (n = 349) and HIV-positive (n = 126) high-risk women with normal cytology.

A total of 126 HIV-positive women were tested for HPV at the Year 2 visit, had valid HPV results, and did not have cervical abnormalities. HPV prevalence among these women was 72.2% for any HPV, 50.8% for HR-HPV and 54.8% for LR-HPV. HIV-positive compared to HIV-negative women were 2.9 times (CI 1.9-4.6) more likely to have any HPV detected; 2.2 times (CI 1.5-3.4) HR-HPV and 2.5 times (CI 1.7-3.9) LR-HPV. The five most prevalent HR-HPV types among HIV-positive women were HPV52 (13.9%), HPV58 (11.1%), HPV51 (11.1%), HPV16, HPV 45 and HPV35 (each 7.9%). Multiple HR-HPV types were detected in 29.4% of HR-HPV positive/HIV-positive women.

Correlates of prevalent high- and low-risk HPV in HIV-negative women

In univariable analysis, the odds of prevalent HR-HPV at Month 6 were significantly lower among older, higher educated, widowed/divorced, and multiparous (≥ 3 pregnancies) women (Table 1). The odds were significantly higher among women with N. gonorrhoeae or C. trachomatis infection at baseline or among women with concurrent BV or LR-HPV infection. Women experiencing incident HSV-2 infection at Month 6 and women using hormonal contraceptices for longer than 6 months were also more likely to have HR-HPV detected at Month 6 (P < 0.05). Sexual behavior indicators (number of clients, number of vaginal sex acts, condom use, age at sexual debut, years as FSW) were not associated with HR-HPV prevalence among this highly sexually active population. In multivariable analysis, being single, having N. gonorrhoeae or C. trachomatis infection at baseline, concurrent LR-HPV infection, and incident HSV-2 infection remained independent predictors of HR-HPV prevalence among HIV-negative women (Table 1).

Table 1 Correlates of HR-HPV in HIV-negative and HIV-positive high-risk women with normal cytology, Kigali, Rwanda

In univariable analysis, LR-HPV was negatively associated with older age (P < 0.05) (data not shown). Hormonal contraceptive use for longer than 6 months and co-infection with HR-HPV were positively associated with LR-HPV (P < 0.05). Age, and concurrent HR-HPV remained independently associated with LR-HPV in multivariable analysis, and there was a trend in that direction for hormonal contraceptive use for longer than 6 months (P = 0.059; data not shown).

Correlates of prevalent high- and low-risk HPV in HIV-positive women

Among the 126 HIV-positive women with no cervical abnormalities who attended the Year 2 visit and had HPV results, prevalent HR-HPV was associated with concurrent LR-HPV infection in uni- and multivariable models (Table 1) There was some evidence of an increased risk of HR-HPV detection among women never/inconsistently using condoms, although this did not reach statistical significance in multivariable model.

Other exposure indicators, CD4 count, use of combination antiretroviral therapy (cART), and the interaction between CD4 count and cART were not associated with HR-HPV detection.

HR-HPV incidence among HIV-negative women

Of the 296 women who contributed to the HPV incidence analyses, incident HR-HPV types were detected in 28.0% (95% CI 23.0-33.1) after a mean interval of 16.6 months (SD 1.8). The mean duration between the two HPV measurements was not significantly associated with acquisition of type-specific HR-HPV (data not shown). Among HR-HPV positive women at Month 6, 34.4% acquired a new HR-HPV type by Year 2, compared to 24.8% among women who were HR-HPV negative at Month 6 (P = 0.087). Cumulative HR-HPV incidence was 39.2% among LR-HPV positive women at Month 6 compared to 22.1% among LR-HPV negative women (P = 0.002).

HPV52 incidence was highest among the HIV-negative high-risk women (Table 2). Women above the age of 30 were less likely to acquire HR-HPV types (Table 3). LR-HPV and HR-HPV types 33 and/or 52 and/or 58 at Month 6,C. trachomatis infection at Month 12, the preceding visit to the second HPV measurement, and concurrent LR-HPV were associated with increased risk of HR-HPV acquisition. Furthermore, women who worked continuously throughout follow-up as a FSW were at increased risk for HR-HPV acquisition. Infection with HPV16-related HPV types 33, 52 and 58 at Month 6, C. trachomatis infection at Month 12, concurrent LR-HPV and self-reported sex work remained independent predictors of HR-HPV incidence in multivariable analysis.

Table 2 Type-specific high-risk HPV prevalence and incidence among HIV-negative high-risk women with normal cytology, Kigali, Rwanda
Table 3 Predictors of HR-HPV incidence and persistence in HIV-negative high-risk women with normal cytology, Kigali, Rwanda

HR-HPV persistence

Ninety HIV-negative women without cervical abnormalities contributed to the HR-HPV persistence analysis. Of them, 32.0% had persistent HR-HPV at Year 2 after a mean interval of 16.9 months (SD 1.7). Mean duration between the two HPV measurements was not associated with clearance or persistence. Persistence was most common for HPV31 and 68 (data not shown). In univariable analyses, only BV at Month 12 was associated with a decreased likelihood of persistent HR-HPV infection (Table 3).


We found high prevalence rates of HR-HPV (50.8% and 31.8%) and LR-HPV (54.8% and 32.4%) in HIV-positive and HIV-negative high-risk Rwandan women, respectively. This is considerably higher than the prevalence rates that we found in women from the general population in Kigali, Rwanda (HR-HPV 19.9%; LR-HPV 20.5%; unpublished data Veldhuijzen et al [14]). Our results are in line with results from another study among HIV-positive women in Kigali that reported a HR-HPV prevalence of 46% [15], as well as studies in other regions that reported HR-HPV prevalence rates ranging from 30-46% in HIV-negative sex workers and 44-73% in HIV-positive sex workers [1517]. In general, increased HPV prevalence among HIV-positive women can be attributed to increased susceptibility, decreased ability to clear infection and frequent reactivation of latent infection associated with immunosuppression [18].

Very few prospective studies on HPV epidemiology have been conducted in Africa. We found a cumulative HR-HPV incidence of 28% and persistence of 32% over a mean interval of 17 months, in HIV-negative high-risk women. In comparison, HR-HPV incidence among a rural Ugandan population was estimated at 7.0 per 100 PY among HIV-negative and 17.3 per 100 PY among HIV-positive women; 66% of these infections were cleared by HIV-negative women compared to 40% by HIV-positive women [19]. HR-HPV incidence among HIV-negative women in Zimbabwe was estimated at 11% and persistence at 37% after 12 months of follow-up [20]. Compared to the other two studies, our study population had higher risk for STIs, including HPV, because of their sexual risk behavior.

HPV52 had the highest incidence rate, the presence of infection by HPV52, 58, and 33 combined at Month 6 was strongly associated with increased acquisition of other HR-HPV types in Year 2. In a previous analysis of the same cohort, the presence of HR-HPV at Month 6 was strongly associated with HIV acquisition, and HPV52 was the most frequent HR-HPV type detected in HIV seroconverters [21]. Infection with multiple HR-HPV types was more common among HIV-positive compared to HIV-negative women. Other studies have also found higher rates of (non-16/18) HR-HPV types in HIV-positive compared to HIV-negative women, as well as higher rates of infections with multiple types [4, 17, 2225].

In our study, prior C. trachomatis and N. gonorrhoeae infection were associated with HR-HPV prevalence in HIV-negative women, and prior C. trachomatis infection with HR-HPV incidence. The literature is mixed, with some studies finding associations between C. trachomatis and HR-HPV [26, 27], and others not [28]. The biological interactions between C. trachomatis infection and HPV infection are not yet fully understood [29, 30]. Inflammatory changes of cervical epithelium, microabrasions of the epithelium and host immune response may facilitate (re) uptake of HPV and/or impede clearance [27, 30]. Most women in our study had either received treatment or spontaneously cleared the C. trachomatis infection by the next visit. Active epithelial regeneration following infection might increase the susceptibility of target cells to HPV. Conversely, insufficient Th1 cells responses (which are required for resolving C. trachomatis infection and clearing HPV) may promote persistent C. trachomatis and persistent HPV infection [3133]. As the majority of C. trachomatis infections in our population were treated or cleared, we were not able to establish an association between persistent C. trachomatis infection and HPV persistence. We also did not find an association between BV and any or HR-HPV, as has been found in other studies [34, 35]. Sexual behavior indicators (number of clients, number of vaginal sex acts, condom use, age at sexual debut, years as FSW) were not associated with HR-HPV prevalence among this highly sexually active population.

There is no consistent definition of HPV persistence in the scientific literature, with some published studies using a duration of 6-12 months, and others using a duration of longer than 12 months [36]. The mean interval between the two HPV measurements in our study was approximately 17 months. With this long duration, the extent of persistence may have been overestimated because we cannot exclude re-activation or re-infection after clearance of the initial infection. In addition, left censoring, by using prevalent HPV infections, may have further inflated our persistence results because long-term infections are often over-represented among prevalent cases. In contrast, the incidence rate may have been underestimated because some women may have acquired new infections but cleared them before the second HPV test was performed. Furthermore, our type-specific prevalence and incidence rates are imprecise due to small sample sizes.


In conclusion, high HPV prevalence, incidence and persistence rates were found in this study among high-risk women in Kigali. HPV52 had the highest incidence; and, together with HPV33 and HPV58, was strongly associated with acquisition of other HR-HPV types in HIV-negative women. HPV52 and 58 were also the most prevalent in HIV-positive women and, as previously shown, in HIV seroconverters. In addition, prior C. trachomatis was associated with HR-HPV incidence.

Women at high-risk for sexually transmitted infections (including HIV) experience a higher incidence and prevalence of HPV infections, including infections with high-risk HPV types. Therefore these women are at greater risk for developing cervical cancer. Cervical cancer prevention programs should incorporate strategies to ensure a high coverage among these high-risk women.


  1. Aral SO, Holmes KK: The epidemiology of STIs and their social and behavioral determinants: industrialized and developing countries. Sexually transmitted Diseases. Edited by: Holmes KK, Sparling PF, Stamm WE, Piot P, Wasserheit JN, Corey L, Cohen MS, Watts DH. 2008, New York: McGraw-Hill, 53-92. 4

    Google Scholar 

  2. 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 (4 Suppl): S5-25. S25 e21-41

    PubMed  Google Scholar 

  3. Bruni L, Diaz M, Castellsague X, Ferrer E, Bosch FX, de Sanjose S: Cervical human papillomavirus prevalence in 5 continents: meta-analysis of 1 million women with normal cytological findings. J Infect Dis. 2010, 202 (12): 1789-99. 10.1086/657321.

    Article  PubMed  Google Scholar 

  4. Clifford GM, Goncalves MA, Franceschi S: Human papillomavirus types among women infected with HIV: a meta-analysis. Aids. 2006, 20 (18): 2337-2344. 10.1097/01.aids.0000253361.63578.14.

    Article  PubMed  Google Scholar 

  5. Juarez-Figueroa LA, Wheeler CM, Uribe-Salas FJ, Conde-Glez CJ, Zampilpa-Mejia LG, Garcia-Cisneros S, Hernandez-Avila M: Human papillomavirus: a highly prevalent sexually transmitted disease agent among female sex workers from Mexico City. Sex Transm Dis. 2001, 28 (3): 125-130. 10.1097/00007435-200103000-00001.

    Article  CAS  PubMed  Google Scholar 

  6. Veldhuijzen NJ, Snijders PJ, Reiss P, Meijer CJ, van de Wijgert JH: Factors affecting transmission of mucosal human papillomavirus. Lancet Infect Dis. 2010, 10 (12): 862-74. 10.1016/S1473-3099(10)70190-0.

    Article  PubMed  Google Scholar 

  7. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM: GLOBOCAN 2008, Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 2008, 10 Lyon, France: International Agency for Research on Cancer

    Google Scholar 

  8. Braunstein Sarah, Ingabire Chantal, Kestelyn Evelyne, Uwizera Aline Umutoni, Mwamarangwe Lambert, Ntirushwa Justin, Nash Denis, Veldhuijzen Nienke, Nel Annalene, Vyankandondera Joseph, van de Wijgert Janneke: High HIV incidence in a cohort of Rwandan female sex workers. Sex Transm Dis. 2011, 38 (5): 384-94.

    Article  Google Scholar 

  9. Stevens MP, Garland SM, Tabrizi SN: Development and validation of a real-time PCR assay specifically detecting human papillomavirus 52 using the Roche LightCycler 480 system. J Virol Methods. 2008, 147 (2): 290-296. 10.1016/j.jviromet.2007.09.018.

    Article  CAS  PubMed  Google Scholar 

  10. Solomon D, Davey D, Kurman R, Moriarty A, O'Connor D, Prey M, Raab S, Sherman M, Wilbur D, Wright T, Young N: The 2001 Bethesda System: terminology for reporting results of cervical cytology. Jama. 2002, 287 (16): 2114-2119. 10.1001/jama.287.16.2114.

    Article  PubMed  Google Scholar 

  11. Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV, Snijders PJ, Meijer CJ: Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003, 348 (6): 518-527. 10.1056/NEJMoa021641.

    Article  PubMed  Google Scholar 

  12. Cogliano V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F: Carcinogenicity of human papillomaviruses. Lancet Oncol. 2005, 6 (4): 204-10.1016/S1470-2045(05)70086-3.

    Article  PubMed  Google Scholar 

  13. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Veldhuijzen NJ, Dhont N, Vyankandondera J, Gasarabwe A, Busasa R, Crucitti T, vdWijgert JHHM: Prevalence and concordance of HPV, HIV and HSV-2 in heterosexual couples in Kigali, Rwanda. Sex Transm Dis - Accepted for publication.

  15. 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 (12): 1851-1861. 10.1086/599123.

    Article  PubMed  PubMed Central  Google Scholar 

  16. del Amo J, Gonzalez C, Belda J, Fernandez E, Martinez R, Gomez I, Torres M, Saiz AG, Ortiz M: Prevalence and risk factors of high-risk human papillomavirus in female sex workers in Spain: differences by geographical origin. J Womens Health (Larchmt). 2009, 18 (12): 2057-2064. 10.1089/jwh.2008.1293.

    Article  Google Scholar 

  17. Luchters SM, Vanden Broeck D, Chersich MF, Nel A, Delva W, Mandaliya K, Depuydt CE, Claeys P, Bogers JP, Temmerman M: Association of HIV infection with distribution and viral load of HPV types in Kenya: a survey with 820 female sex workers. BMC Infect Dis. 2010, 10: 18-10.1186/1471-2334-10-18.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Palefsky J: Human papillomavirus-related disease in people with HIV. Curr Opin HIV AIDS. 2009, 4 (1): 52-56. 10.1097/COH.0b013e32831a7246.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Safaeian M, Kiddugavu M, Gravitt PE, Gange SJ, Ssekasanvu J, Murokora D, Sklar M, Serwadda D, Wawer MJ, Shah KV, 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 (6): 1300-1307. 10.1158/1055-9965.EPI-07-2678.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Fukuchi E, Sawaya GF, Chirenje M, Magure T, Tuveson J, Ma Y, Shiboski S, Da Costa M, Palefsky J, Moscicki AB, Makunike-Mutasa R, Chipato T, Smith-McCune KK: Cervical human papillomavirus incidence and persistence in a cohort of HIV-negative women in Zimbabwe. Sex Transm Dis. 2009, 36 (5): 305-311. 10.1097/OLQ.0b013e318194eb76.

    Article  PubMed  Google Scholar 

  21. Veldhuijzen NJ, Vyankandondera J, van de Wijgert JH: HIV acquisition is associated with prior high-risk human papillomavirus infection among high-risk women in Rwanda. Aids. 2010, 24 (14): 2289-2292. 10.1097/QAD.0b013e32833cbb71.

    Article  PubMed  Google Scholar 

  22. Ahdieh L, Klein RS, Burk R, Cu-Uvin S, Schuman P, Duerr A, Safaeian M, Astemborski J, Daniel R, Shah K: Prevalence, incidence, and type-specific persistence of human papillomavirus in human immunodeficiency virus (HIV)-positive and HIV-negative women. J Infect Dis. 2001, 184 (6): 682-690. 10.1086/323081.

    Article  CAS  PubMed  Google Scholar 

  23. 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 (2): 137-142. 10.1097/00007435-200302000-00009.

    Article  PubMed  Google Scholar 

  24. Mbulawa ZZ, Coetzee D, Marais DJ, Kamupira M, Zwane E, Allan B, Constant D, Moodley JR, Hoffman M, Williamson AL: Genital human papillomavirus prevalence and human papillomavirus concordance in heterosexual couples are positively associated with human immunodeficiency virus coinfection. J Infect Dis. 2009, 199 (10): 1514-1524. 10.1086/598220.

    Article  PubMed  Google Scholar 

  25. McKenzie ND, Kobetz EN, Hnatyszyn J, Twiggs LB, Lucci JA: Women with HIV are more commonly infected with non-16 and -18 high-risk HPV types. Gynecol Oncol. 116 (3): 572-577.

  26. Samoff E, Koumans EH, Markowitz LE, Sternberg M, Sawyer MK, Swan D, Papp JR, Black CM, Unger ER: Association of Chlamydia trachomatis with persistence of high-risk types of human papillomavirus in a cohort of female adolescents. Am J Epidemiol. 2005, 162 (7): 668-675. 10.1093/aje/kwi262.

    Article  PubMed  Google Scholar 

  27. Verhoeven V, Baay M, Weyler J, Avonts D, Lardon F, Van Royen P, Vermorken JB: Concomitant Chlamydia trachomatis and human papilloma virus infection cannot be attributed solely to sexual behaviour. Eur J Clin Microbiol Infect Dis. 2004, 23 (9): 735-737.

    Article  CAS  PubMed  Google Scholar 

  28. Nielsen A, Iftner T, Munk C, Kjaer SK: Acquisition of high-risk human papillomavirus infection in a population-based cohort of Danish women. Sex Transm Dis. 2009, 36 (10): 609-615. 10.1097/OLQ.0b013e3181a96d0e.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Burchell AN, Winer RL, de Sanjose S, Franco EL: Chapter 6: Epidemiology and transmission dynamics of genital HPV infection. Vaccine. 2006, 24 (Suppl 3): S3/52-61.

    Google Scholar 

  30. Simonetti AC, Melo JH, de Souza PR, Bruneska D, de Lima Filho JL: Immunological's host profile for HPV and Chlamydia trachomatis, a cervical cancer cofactor. Microbes Infect. 2009, 11 (4): 435-442. 10.1016/j.micinf.2009.01.004.

    Article  CAS  PubMed  Google Scholar 

  31. Loomis WP, Starnbach MN: T cell responses to Chlamydia trachomatis. Curr Opin Microbiol. 2002, 5 (1): 87-91. 10.1016/S1369-5274(02)00291-6.

    Article  CAS  PubMed  Google Scholar 

  32. Scott M, Stites DP, Moscicki AB: Th1 cytokine patterns in cervical human papillomavirus infection. Clin Diagn Lab Immunol. 1999, 6 (5): 751-755.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Moscicki AB, Schiffman M, Kjaer S, Villa LL: Chapter 5: Updating the natural history of HPV and anogenital cancer. Vaccine. 2006, 24 (Suppl 3): S3/42-51.

    Google Scholar 

  34. Lippman SA, Sucupira MC, Jones HE, Luppi CG, Palefsky J, van de Wijgert JH, Oliveira RL, Diaz RS: Prevalence, distribution and correlates of endocervical human papillomavirus types in Brazilian women. Int J STD AIDS. 2010, 21 (2): 105-109. 10.1258/ijsa.2009.008436.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Watts DH, Fazzari M, Minkoff H, Hillier SL, Sha B, Glesby M, Levine AM, Burk R, Palefsky JM, Moxley M, Ahdieh-Grant L, Strickler HD: Effects of bacterial vaginosis and other genital infections on the natural history of human papillomavirus infection in HIV-1-infected and high-risk HIV-1-uninfected women. J Infect Dis. 2005, 191 (7): 1129-1139. 10.1086/427777.

    Article  PubMed  Google Scholar 

  36. Koshiol J, Lindsay L, Pimenta JM, Poole C, Jenkins D, Smith JS: Persistent human papillomavirus infection and cervical neoplasia: a systematic review and meta-analysis. Am J Epidemiol. 2008, 168 (2): 123-137. 10.1093/aje/kwn036.

    Article  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references


The authors acknowledge the following people and institutions for their contributions to this study: all study participants; the Projet Ubuzima study team, Governing Council, Community Advisory Group, and community mobilizers; laboratory staff at the National Reference Laboratory in Kigali; Institute of Tropical Medicine, Antwerp, Belgium (Said Abdellati and Vicky Cuylaert), Ghent University Hospital (Prof. Marleen Praet), AMC (Anneke Oei), and the Amsterdam Municipal Health Service- (Ineke Linde, Anita Hadisusanto, and Ilse van Meelt); Robert Meester; the clinical affairs team at the International Partnership for Microbicides, Inc.; and our collaborators within the EDCTP-funded project "Preparing for Phase III vaginal microbicide trials in Rwanda and Kenya: Preparedness studies, capacity building, and strengthening of medical referral systems".


The work described in this paper was funded by the European and Developing Countries Clinical Trials Partnership (EDCTP) through a project entitled: "Preparing for Phase III vaginal microbicide trials in Rwanda and Kenya: Preparedness studies, capacity building, and strengthening of medical referral systems" (CT_ct_05_33070_001). EDCTP cannot accept any responsibility for information or views expressed herein. The main HIV incidence study, of which the reported work was a substudy, was also supported by the International Partnership for Microbicides, Inc. NV received a travel grant from the Van Walree fund for an extended stay in Rwanda during implementation of the HPV-related work.

Previous presentations

Preliminary HPV prevalence results were presented during an oral presentation at Eurogin Conference, Nice, France, 12-15 November 2008 (Abstract FC-6a-13). A poster presentation with results on the association between HPV prevalence and concurrent STIs was presented at ISSTDR, London, UK, 28 June - 3 July 2009. HPV incidence, prevalence and persistence results were presented during an oral presentation at the 27th International Papillomavirus Conference in Montreal, 3-8 July 2010.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Nienke J Veldhuijzen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JvdW, JV and SB have been involved in the conception and design; CI, JN, EK, CT, AU, MU, CT, JV, SB and NV in the acquisition of data; and NV, FW and JvdW in data analysis and interpretation. NV took the lead in drafting the manuscript, all authors critically reviewed manuscripts drafts and approved the final version.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Veldhuijzen, N.J., Braunstein, S.L., Vyankandondera, J. et al. The epidemiology of human papillomavirus infection in HIV-positive and HIV-negative high-risk women in Kigali, Rwanda. BMC Infect Dis 11, 333 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: