Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Prevalence, incidence and correlates of low risk HPV infection and anogenital warts in a cohort of women living with HIV in Burkina Faso and South Africa

  • Admire Chikandiwa ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Data curation, Formal analysis, Investigation, Project administration, Writing – original draft

    achikandiwa@wrhi.ac.za

    Affiliation Wits Reproductive Health and HIV Institute, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa

  • Helen Kelly ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Formal analysis, Project administration, Writing – review & editing

    Affiliation Clinical Research Department, London School of Hygiene & Tropical Medicine, London, United Kingdom

  • Bernard Sawadogo ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Data curation, Investigation, Project administration

    Affiliation Centre de Recherche Internationale pour la Santé, Université de Ouagadougou, Ouagadougou, Burkina Faso

  • Jean Ngou ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Investigation

    Affiliation INSERM, EFS, University of Montpellier, Montpellier, France

  • Pedro T. Pisa ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Formal analysis, Writing – review & editing

    Affiliation Wits Reproductive Health and HIV Institute, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa

  • Lorna Gibson ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Data curation, Formal analysis

    Affiliation Clinical Research Department, London School of Hygiene & Tropical Medicine, London, United Kingdom

  • Marie-Noelle Didelot ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Formal analysis, Investigation, Methodology

    Affiliation INSERM, EFS, University of Montpellier, Montpellier, France

  • Nicolas Meda ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Conceptualization, Methodology, Resources, Supervision

    Affiliation Centre de Recherche Internationale pour la Santé, Université de Ouagadougou, Ouagadougou, Burkina Faso

  • Helen A. Weiss ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Conceptualization, Methodology, Supervision, Validation

    Affiliation Clinical Research Department, London School of Hygiene & Tropical Medicine, London, United Kingdom

  • Michel Segondy ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Conceptualization, Resources, Supervision

    Affiliation INSERM, EFS, University of Montpellier, Montpellier, France

  • Philippe Mayaud ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing

    Affiliations Wits Reproductive Health and HIV Institute, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa, Clinical Research Department, London School of Hygiene & Tropical Medicine, London, United Kingdom

  • Sinead Delany-Moretlwe ,

    Contributed equally to this work with: Admire Chikandiwa, Helen Kelly, Bernard Sawadogo, Jean Ngou, Pedro T. Pisa, Lorna Gibson, Marie-Noelle Didelot, Nicolas Meda, Helen A. Weiss, Michel Segondy, Philippe Mayaud, Sinead Delany-Moretlwe

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Resources, Supervision, Writing – review & editing

    Affiliation Wits Reproductive Health and HIV Institute, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa

  • on behalf of the HARP Study Group

    Membership of the HARP Study Group is provided in the Acknowledgments.

Abstract

Objective

To report the prevalence and incidence of low-risk human papillomavirus infection (LR-HPV) and anogenital warts (AGW) among women living with HIV (WLHIV) in Burkina Faso (BF) and South Africa (SA), and to explore HIV-related factors associated with these outcomes.

Methods

We enrolled 1238 WLHIV (BF = 615; SA = 623) aged 25–50 years and followed them at three time points (6, 12 and 16 months) after enrolment. Presence of AGW was assessed during gynaecological examination. Cervico-vaginal swabs for enrolment and month 16 follow-up visits were tested for HPV infection by Inno-LiPA® genotyping. Logistic regression was used to assess risk factors for prevalent infection or AGW. Cox regression was used to assess risk factors for incident AGW.

Results

Women in SA were more likely than those in BF to have prevalent LR-HPV infection (BF: 27.1% vs. SA: 40.9%; p<0.001) and incident LR-HPV infection (BF: 25.8% vs. SA: 31.6%, p = 0.05). Prevalence of persistent LR-HPV was similar in the two countries (BF: 33.3% vs. SA: 30.4%; p = 0.54), as were prevalence and incidence of AGW (Prevalence: BF: 7.5% vs. SA: 5.7%; p = 0.21; Incidence: BF: 2.47 vs. SA: 2.33 per 100 person-years; p = 0.41). HPV6 was associated with incident AGW (BF: adjusted Hazard Ratio (aHR) = 4.88; 95%CI: 1.36–17.45; SA: aHR = 5.02; 95%CI: 1.40–17.99). Prevalent LR-HPV (BF: adjusted Odds Ratio [aOR = 1.86]; 95%CI: 1.01–3.41; SA: aOR = 1.75; 95%CI: 0.88–3.48); persistent LR-HPV (BF: aOR = 1.92; 95%CI: 0.44–8.44; SA: aOR = 2.81; 95%CI: 1.07–7.41) and prevalent AGW (BF: aOR = 1.53; 95%CI: 0.61–3.87; SA: aOR = 4.11; 95%CI: 1.20–14.10) were each associated with low CD4+ counts (i.e. <200 vs. >500 cells/μL). Duration of ART and HIV plasma viral load were not associated with any LR-HPV infection or AGW outcomes.

Conclusion

LR-HPV infection and AGW are common in WLHIV in sub-Saharan Africa. Type-specific HPV vaccines and effective ART with immunological reconstitution could reduce the burden of AGW in this population.

Introduction

Anogenital human papillomavirus (HPV) infection is the most frequent sexually transmitted infection (STI) globally.[1] Low Risk (LR) HPV types 6 and 11 are the predominant causes of anogenital warts (AGW), a source of psychosocial distress,[2] and physical discomfort (including bleeding and itching).[3] Approximately 25% of AGW spontaneously regress,[3] but recurrence is common, resulting in high medical costs from repeated treatment.[4]

There is little literature on the epidemiology of LR-HPV infection and AGW in Sub-Saharan Africa. Regional variations in LR-HPV prevalence varies from 2% in Western Africa, 4% in Southern Africa and 8% in Eastern Africa.[5] Similary, the prevalence of AGW ranges from 4–11% in Western, 2–14% in Southern and 3–11% in Eastern African regions.[6] The natural history of HPV infection is altered by HIV infection, and HPV prevalence is higher where HIV is more common or among individuals with HIV infection.[7] People living with HIV who have AGW tend to experience florid and prolonged clinical manifestations of AGW because of their impaired immune response.[8] Low CD4+ counts (<200 cells/μL) are associated with an increased risk of AGW.[9, 10] The impact of antiretroviral therapy (ART) on the prevalence and incidence of AGW is unclear, with conflicting findings from different studies. Some report that ART reduces incidence,[10] while other studies have shown that it does not reduce incidence.[7, 11]

Three vaccines against HPV have been licensed,[12] two of which (quadrivalent and nonavalent) have activity against HPV types 6 and 11 that cause 90% of AGW.[13, 14] These vaccines have been steadily introduced into national immunization programmes in both high and low and middle income countries (LMIC), although the pace of introduction has been slower in LMIC.[15] Data from countries which have included the quadrivalent vaccine (against HPV6/11/16/18 types) in their vaccination programme show significant declines in the number of cases of AGW among vaccinated women and some herd-immunity effect in unvaccinated heterosexual men.[16] Decisions about which vaccine to include in a vaccination programme are driven by cost-effectiveness considerations. Since a high-burden of AGW is associated with significant costs to the health system in LMIC and this is worsened by high HIV prevalence, data on the prevalence, incidence and correlates of LR-HPV and AGW is important for informing vaccine programme decision-making and resource allocation.[12] To address this question, we analysed data from a cohort of women living with HIV (WLHIV) who participated in the HPV in Africa Research Partnership (HARP) study of cervical cancer screening in Burkina Faso and South Africa.

Materials and methods

Study design, population and sample collection

We conducted a prospective cohort study among WLHIV aged 25–50 years and residing in Ouagadougou, Burkina Faso (BF) and Johannesburg, South Africa (SA) between December 2011 and October 2012. A full description of the study design and procedures has been published elsewhere.[17] Briefly, enrolled participants were followed-up every 6 months up to 18 months. Data on clinical, socio-demographic and behavioural characteristics were collected by interviewer-administered questionnaire at each visit. The presence of AGW was assessed at enrolment and months 6, 12 and endline by trained nurse midwives during standardised pelvic examination. Cervical samples were collected using a Digene cervical sampler (Qiagen, Courtaboeuf, France) for HPV-DNA testing and genotyping, a swab from the ecto/endocervix to detect cervical sexually transmitted infections (STIs) by molecular methods and a vaginal smear to diagnose bacterial vaginosis and Candida albicans by Gram stain at both enrolment and endline visits. Blood samples were collected at enrolment to confirm HIV-1 serostatus, perform Herpes Simplex Virus-2 (HSV-2) and syphilis serologies as well as obtain baseline HIV-1 plasma viral load (PVL) measurement and CD4+ T-lymphocytes counts. The monitoring of CD4+ counts was done at subsequent 6-monthly visits.

Laboratory methods

HIV-1 serostatus was diagnosed using two rapid tests according to national guidelines.[18, 19] Testing for CD4+ cells was performed using FACScount (Becton-Dickinson, NJ). Plasma HIV-1 RNA was assessed using real-time PCR (Abbott RT HIV-1) in BF and COBAS Taqman (Roche Diagnostics) in SA, with a lower limit of detection of 40 copies/ml in both countries. Laboratories subscribed to international external quality assessment schemes such as the UK-NEQAS for CD4+ counts[20] and QCMD for HIV-1 PVL testing.[21]

HSV-2 serology was performed using the Kalon® gG2 ELISA (Kalon Diagnostics, UK) and syphilis serology by a combination of a Treponema pallidum haemagglutination (TPHA) and rapid plasma reagin (RPR; BioMérieux, Lyon, France in BF and Immutrep carbon antigen RPR, Omega Diagnostics in SA). Neisseria gonorrhoeae, Chlamydia trachomatis, Mycoplasma genitalium and Trichomonas vaginalis were detected using nucleic acid amplification tests, the Sacace simplex assays (Sacace, Como, Italy) in BF and the APTIMA Combo (Gen-Probe, San Diego, CA) in SA. The Nugent’s score was used for vaginal flora reading of Gram-stained vaginal smears, with diagnosis of bacterial vaginosis made for scores ≥7, and examined for presence of Candida.

HPV DNA testing and genotyping was performed at the virology laboratory of University of Montpellier (UM) using INNO-LiPA HPV Genotyping Extra Assay (Innogenetics, Courtaboeuf, France). The assay is based on the amplification of a 65-bp fragment in the L1 gene and detects 28 HPV types, including all the high-risk or probable high-risk types (16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, 82) as well as the following LR HPV genotypes (6, 11, 40, 43, 44, 54, 69/71, 70,74).[22]

Statistical analysis

Any LR-HPV type infection was defined as detection of 1 of the 9 genotypes described above. LR-HPV type-specific persistence was defined as being positive for the same type at enrolment and endline. Type-specific LR-HPV incidence was defined as the proportion of women who were negative for a specific LR type at enrolment and positive for that type at endline. Prevalent AGW were defined as the presence of clinically-defined condyloma acuminata at the enrolment visit. Incident AGW were defined as the first documented occurrence of condyloma acuminata in a woman previously documented not to have AGW. Time to incident AGW was calculated from date of enrolment visit to the date of first AGW detection.

Descriptive statistics were used to summarize the prevalence of any LR-HPV infection and AGW. For bivariate analysis, socio-demographic, sexual behaviour, presence of other laboratory detected STI (as part of the screening process at enrolment), clinical signs and HIV-related risk factors were chosen for assessment a priori as known correlates for HPV and AGW acquisition. As LR-HPV prevalence was common, associations with exposure variables were estimated from prevalence ratios obtained from logistic regression using marginal standardization to estimate prevalence ratios, and the delta method to estimate 95% confidence intervals (CIs). Associations between LR-HPV persistence and exposure variables were estimated with generalized estimating equations to account for multiple LR-HPV infection and multiple infection states (persistence and clearance).[23] To explore associations of any LR-HPV and AGW with HIV-related factors, pre-specified analyses included stratification by site, ART use, duration on ART (≤ or >2 years), HIV-1 viral suppression (< or ≥1000 copies/ml) and CD4+ cell counts at enrolment. Logistic regression (Model 1) was used to identify independent associations between any prevalent LR-HPV or AGW and HIV-related factors, after adjustment for potential confounders (non HIV-related factors) that were associated with the outcomes in bivariate analysis at p<0.10.[24] A second logistic regression model (Model 2) adjusted for CD4+ cell count at enrolment in addition to factors adjusted for in Model 1 was run to explore associations with ART duration and HIV-1 viral suppression. Cox regression was performed to analyse predictors of incident AGW, with and without adjustment for potential confounders using the same model-building approach as described above under logistic regression. Data were analysed using Stata version 14 (Stata Statistical Software, College Station. TX: Stata Corporation).

Ethics statement

Study protocols were approved by the Ministry of Health in Burkina Faso (no. 2012-12-089), the University of the Witwatersrand in South Africa (no. 110707), and the London School of Hygiene and Tropical Medicine (no. 7400). All participants provided written informed consent.

Results

Study population

Of the 1473 women screened, 1238 (615 in BF and 623 in SA) were enrolled. The median age of participants was 36 (interquartile range [IQR], 31–42) years in BF and 34 (IQR, 30–40) years in SA. At enrolment, two-thirds of the participants were on ART (BF: 68.6% in BF; 65.2% in SA) and this is because enrolment was stratified in the ratio of 2 ART users: 1 non-ART user. In BF, the median CD4+ count was 417 (IQR, 315–606) cells/μL among ART-naive participants and 446 (IQR, 309–600) cells/μL among those on ART. In SA, the median CD4+ count was 448 (IQR, 353–614) cells/μL among ART-naive and 420 (IQR, 279–567) cells/μL among those on ART. A more detailed description of study participants has been published elsewhere.[17]

Prevalence and correlates of LR-HPV infection

Almost all participants (n = 1215; 98.1%) had valid HPV genotyping results [BF = 594 (96.6%); SA = 621 (99.7%)]. The prevalence of any LR-HPV was lower in BF (27.1%) than in SA (40.9%) (p<0.001), as were the prevalence of the following LR types: 11 (1.5% vs. 5.3%, p<0.001), 44 (7.7% vs. 11.6%, p = 0.02), 70 (4.0% vs. 8.9%, p = 0.001) and 74 (6.2% vs. 10.0%, p = 0.02). Type 6 prevalence was similar in both countries (BF 5.7% vs. SA 5.3%). Type 44 was the most prevalent LR-HPV type, for both countries (Fig 1A). In both countries, LR-HPV prevalence was lower among women who have had pregnancies compared to those who have been never pregnant (BF: adjusted Prevalence Ratio (aPR) = 0.54; 95%CI: 0.29–0.99 vs. SA: aPR = 0.62; 95%CI: 0.44–0.89). In BF alcohol use (aPR = 1.41; 95%CI: 1.01–1.97) and self-reported consistent condom use (aPR = 1.93; 95%CI: 1.07–3.47) were associated with higher LR-HPV prevalence. In SA LR-HPV prevalence was higher among women without a regular sexual partner compared to those with a regular partner (aPR = 2.39; 95%CI: 1.37–4.19) (Table 1).

thumbnail
Fig 1. Type specific LR-HPV infection.

(A) Prevalence at enrolment. (B) Incidence over 16 months. (C) Persistence over 16 months. *Significantly different between the two countries (i.e. p-value <0.05). **Persistence as a proportion of those positive for that specific type at enrolment.

https://doi.org/10.1371/journal.pone.0196018.g001

thumbnail
Table 1. The prevalence and correlates of LR-HPVa infection among 594 WLHIV in BF and 621 in SA.

https://doi.org/10.1371/journal.pone.0196018.t001

Associations between prevalent LR-HPV infection and HIV-related factors

Of the 1215 women with valid HPV results, 32 (2.6%) reported that they were not on ART at enrolment despite undetectable HIV-1 PVL (BF = 24; SA = 8). Given the uncertainty of their ART status, these women were excluded from HIV-related factors analyses. Further analyses were conducted on the remaining 1183 participants (BF = 570; SA = 613) to determine the effect of HIV-related factors on any LR-HPV prevalence (Table 2). In both countries, LR-HPV prevalence was higher among women with a low CD4+ count (<200 cells/μL vs. >500 cells/μL; Model 1 BF: adjusted Odds Ratio (aOR) = 1.37; 95%CI: 1.13–1.65 vs. Model 1 SA: aOR = 1.15; 95%CI: 0.99–1.33) (Table 2). Among ART users, a low CD4+ count was also associated with a higher risk of prevalent LR-HPV infection (Model 1 BF: aOR = 1.41; 95%CI: 0.89–2.24 vs. Model 1 SA: aOR = 1.58; 95%CI: 1.06–2.34). The duration on ART (i.e. >2 years vs. ≤2 years on ART) and HIV-1 PVL at enrolment were not significantly associated with prevalent LR-HPV infection (Table 2).

thumbnail
Table 2. Associations between prevalent LR-HPVa and HIV-related factors among 570 WLHIV in BF and 613 in SA.

https://doi.org/10.1371/journal.pone.0196018.t002

Incidence and persistence of LR-HPV infection and effect of HIV-related factors

Overall, 963/1238 (77.8%) WLHIV attended the endline visit at a median 16 (IQR, 15.6–16.8) months and the majority (922) had HPV-genotyping results at both enrolment and endline (BF = 476; SA = 446). Of these, 264 (28.7%) had incident LR-HPV of any type during follow-up. The incidence of any LR-HPV infection was lower in BF compared to SA (BF: 25.8% vs. SA: 31.6%, p = 0.05). This was also true for most individual types (11, 54 and 74) with the exception of types 44 (BF: 4.7% vs. SA: 0.9%, p = 0.01) and 69/71 (BF: 8.5% vs. SA: 4.5%, p = 0.02) which were significantly higher in BF compared to SA (Fig 1B).

A third of WLHIV (310/922) were defined as having persistent LR-HPV infection; the total number of persistent infections was 396 (BF = 156; SA = 240). The frequency of persistent infection of any LR-HPV at endline was similar in both countries (BF: 33.3% vs. SA: 30.4%; p = 0.54). Persistence for individual types was largely similar between the two countries (Fig 1C). A low CD4+ count (i.e. <200 cells/μL vs. >500 cells/μL) was associated with a higher likelihood of having persistent LR-HPV infection. The association was stronger among ART users (Model 1 BF: aOR = 3.21; 95%CI: 0.90–11.40 vs. Model 1 SA: aOR = 3.33; 95%CI: 0.96–11.55). The duration on ART at enrolment and HIV-1 PVL were not significantly associated with persistent LR HPV infection (S1 Table).

Prevalent AGW and effect of prevalent LR-HPV infection and HIV-related factors

Overall, 81 (6.7%) WLHIV had AGW detected at enrolment; the prevalence did not vary by country (BF: 7.5% vs. SA: 5.7%; p = 0.21). In univariate analysis (data not shown), prevalent AGW were associated with age, number of partners in the last 3 months before enrolment, presence of genital ulcers for both countries and being HSV-2 seropositive in BF or having bacterial vaginosis in SA. Infection with any LR-HPV type, multiple LR-HPV infections, and infection with types 6 (Model 2 BF aOR = 3.14; 95%CI: 1.11–9.01 and SA aOR = 1.79; 95%CI: 0.49–6.51) and 11 (Model 2 BF aOR = 20.12; 95%CI: 4.06–99.71 and SA aOR = 7.90; 95%CI: 3.11–20.09) were significantly associated with prevalent AGW (Table 3). In BF, HPV type 44 infection (Model 2 aOR = 3.59; 95%CI: 1.53–8.39), was strongly associated with prevalent AGW. Similary in SA, HPV type 74 infection (Model 2 aOR = 3.28; 95%CI: 1.35–7.98) was also significantly associated with prevalent AGW. AGW prevalence was higher among those with low CD4+ counts (i.e. <200 cells/μL vs. >500 cells/μL), particularly in SA among all women (aOR = 4.11; 95%CI: 1.20–14.10) and among ART users only (aOR = 3.93; 95%CI: 1.01–15.28). AGW prevalence was higher among those who had been on ART for two years or less compared to those on ART for more than two years at enrolment in both BF (Model 1 aOR = 1.78; 95%CI: 0.87–3.64) and SA (Model 1 aOR = 1.88; 95%CI: 0.80–4.38), although the association was less pronounced when adjusted for CD4+ count (Table 3).

thumbnail
Table 3. Effect of HIV-related factors and LR-HPV on AGW prevalence at enrolment among 572 WLHIV in BF and 613 in SA.

https://doi.org/10.1371/journal.pone.0196018.t003

AGW incidence and association with LR HPV infection and HIV-related factors

A total of 32 new AGW events (BF = 16; SA = 16) were recorded for the 992 (BF = 468; SA = 524) WLHIV who were followed up for a total of 1390 person-years. The overall AGW incidence was 2.30 per 100 person-years (BF: 2.47 vs. SA: 2.33 per 100 person-years, p = 0.41). In univariate analysis (data not shown), incident AGW were associated with age for both countries and number of partners in the last 3 months before enrolment in BF or having mycoplasma genitalium in SA. Infection with any LR-HPV at enrolment was associated with having an incident AGW; this association was statistically significant in BF (Model 1 adjusted Hazard Ratio [aHR] = 5.78; 95%CI: 2.13–16.22) (Table 4). Prevalent type 6 infection at enrolment was associated with an almost 5-fold higher odds of incident AGW in both BF (Model 1 aHR = 4.88; 95%CI: 1.36–17.45) and SA (Model 1 aHR = 5.02; 95%CI: 1.40–17.99). In BF, prevalent HPV 74 infection was also significantly associated with incident AGW (Model 1 aHR = 6.02; 95%CI: 1.90–19.04). The incidence of AGW was higher in women with lower CD4+ counts, with an inverse dose-response relationship in SA; but the duration on ART at enrolment and HIV-1 PVL enrolment were not significantly associated with AGW incidence, even after adjustment for CD4+ count (Model 2).

thumbnail
Table 4. Effect of HIV-related factors and LR-HPV infection on incident AGW among 468 WLHIV in BF and 524 in SA.

https://doi.org/10.1371/journal.pone.0196018.t004

Discussion

We found that the prevalence and incidence of LR-HPV infection in this population were high and consistent with other studies of WLHIV in sub-Saharan Africa.[7, 25, 26] We also found that infections with types 6 and 11 were strongly associated with both prevalent and incident (type 6 only) AGW, which confirms the literature.[7, 27] The prevalence of AGW in both countries is similar to what has been reported among HIV-1 seropositive female sex workers in Burkina Faso, [7] but the AGW incidence (BF = 2.47 and SA = 2.33 per 100 person-years) was almost 3-fold lower than the 7.14 per 100 person-years reported by Low et al.,[7] possibly owing to the repeated occupational exposure and lower proportion of WLHIV on ART (10%) in that population.

We found that low enrolment CD4+ cell count was associated with a higher prevalence and persistence of LR-HPV as well as a higher prevalence and incidence of AGW. This is consistent with previous studies reporting CD4+ count as a strong predictor of HPV infection and AGW development.[9, 10] A longer duration (>2 years) of ART was not protective per se against prevalent LR-HPV infection, incident LR-HPV infection and incident AGW, confirming other studies findings.[7, 11] This suggests that other co-factors such as adherence to ART, HIV-1 PVL suppression and immunological reconstitution (as measured by increases in CD4+ count) are also important for the control of HPV infection and progression to associated disease.

We also found that, as reported in previous studies, non-HIV related factors such as alcohol use, number of recent sexual partners, changes in vaginal flora and coinfection with other STIs were significant co-factors for LR-HPV infection and AGW.[6, 7] In contrast to reports from previous studies, we did not find significant association between smoking or hormonal contraception with LR-HPV infection or AGW.[6, 28, 29] It was also surprising to note that LR-HPV prevalence was high among women who reported consistent condom use. The reason why infection prevalence was high among women who reported consistent condom use is not known. However, this could be related to social desirability bias which meant that women with high risk sexual behaviours reported high condom use.

Limitations of the study include the absence of HIV-1 PVL monitoring, in particular at endline and the possible subjective nature of clinical examination for AGW detection without histological confirmation leading to possible misclassification. Cervical LR-HPV infection rather than vaginal/vulval LR-HPV infection was measured which would have been more appropriate given the frequent external location of AGW. This can lead to difficulties in ascribing causality. However, cervical HPV infection has been used previously in determining associations with AGW.[7] Despite these limitations, this study had a number of strengths, including its longitudinal design, availability of AGW data at four time points and availability of genotyping data at both enrolment and endline visits. More importantly, the study contributes data on the epidemiology of LR-HPV infection and AGW from two sub-Saharan African countries with different HIV epidemics and this allows the findings to be extended to a range of countries and settings in the region.

Conclusions

This study shows that WLHIV in sub-Saharan Africa are at high risk of LR-HPV infections and AGW. AGW were associated with HPV 6/11 and these findings support the argument for the use of quadri- or nano-valent HPV vaccines for the prevention of AGW among WLHIV. Effective use of ART with immunological reconstitution could also help control the burden of AGW among WLHIV in the region.

Supporting information

S1 Table. Effect of HIV-related factors on LR-HPV persistence, using infections as unit of measure.

Adjusted Odds Ratio (aOR) using generalised estimating equations: aIn BF, associations with LR-HPV were adjusted for alcohol use, condom use and number of pregnancies; bIn SA, associations with LR-HPV were adjusted for number of regular sex partners, number of pregnancies and number of lifetime sex partners; For Model 2: same adjustment as Model 1, in addition to enrolment CD4+ count; c Of the 156 infections in Burkina Faso at enrolment (Table 3), 5 were excluded from analysis as they were ART-naïve but with undetectable viral load and similarly of the 240 infections at enrolment in South Africa (Table 3), 2 were excluded from the analysis; dART use was defined as being on ART at both enrolment and endline; *All factors with p<0.05 are in bold.

https://doi.org/10.1371/journal.pone.0196018.s001

(DOC)

Acknowledgments

The research leading to these results has received funding from the European Commission (EC) 7th Framework Programme under grant agreement No. HEALTH-2010-F2-265396.

Contributing members of the HARP study group included: A. Chikandiwa, E. Cutler, S. Delany-Moretlwe, L. Diep, C. Hwinya, D. A. Lewis, M.P. Magooa, V. Maseko, P. Michelow, B. Muzah, T. Omar, A. Puren (Johannesburg, South Africa); F. Djigma, J. Drabo, O. Goumbri-Lompo, D. Karambiri, N. Meda, B. Sawadogo, J. Simporé, R. Sombie, A. Yonli, Z. Zamane, S Zan (Ouagadougou, Burkina Faso); V. Costes, M.N. Didelot, S. Doutre, N. Leventoux, N. Nagot, J. Ngou, M. Segondy (Montpellier, France); and A. Devine, C. Gilham, L. Gibson, H. Kelly, R. Legood, P. Mayaud, H.A. Weiss (London, UK).

The HARP Study Group also wishes to thank its International Scientific Advisory Group (ISAG) constituting of Prof. C. Lacey (Chair, University of York, UK), Prof. Y. Qiao (Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China), Prof. M. Chirenje (University of Harare, Zimbabwe) and Prof. S. de Sanjosé (Institut Catala d’Oncologia, Barcelona, Spain).

The HARP Study Group also wishes to thank Prof. W. Prendiville (International Federation of Cervical Pathology and Colposcopy) and Dr A. Olaitan (Consultant Gynaecological Oncologist at the University College London Hospital) who acted as external clinical advisors, and Dr O. E Clavero (Institut Catala d’Oncologia (Barcelona, Spain) and H. Perocchia (Department of Cytology and Anatomical Pathology, Gui de Chauliac Hospital, Montpellier, France) who acted as independent readers for the Endpoint Committee.

References

  1. 1. Schiffman M, Castle PE. Human papillomavirus: epidemiology and public health. Archives of pathology & laboratory medicine. 2003;127(8):930–4. Epub 2003/07/23. pmid:12873163.
  2. 2. Jeynes C, Chung MC, Challenor R. 'Shame on you'—the psychosocial impact of genital warts. Int J STD AIDS. 2009;20(8):557–60. Epub 2009/07/25. pmid:19625588.
  3. 3. Wiley DJ, Douglas J, Beutner K, Cox T, Fife K, Moscicki AB, et al. External genital warts: diagnosis, treatment, and prevention. Clin Infect Dis. 2002;35(Suppl 2):S210–24. Epub 15 Oct 2002. pmid:12353208.
  4. 4. Insinga RP, Dasbach EJ, Elbasha EH. Assessing the annual economic burden of preventing and treating anogenital human papillomavirus-related disease in the US: analytic framework and review of the literature. Pharmacoeconomics. 2005;23(11):1107–22. pmid:16277547
  5. 5. De Vuyst H, Alemany L, Lacey C, Chibwesha CJ, Sahasrabuddhe V, Banura C, et al. The Burden of Human Papillomavirus Infections and Related Diseases in Sub-Saharan Africa. Vaccine. 2013;31, Supplement 5:F32–F46. http://dx.doi.org/10.1016/j.vaccine.2012.07.092.
  6. 6. Banura C, Mirembe FM, Orem J, Mbonye AK, Kasasa S, Mbidde EK. Prevalence, incidence and risk factors for anogenital warts in Sub Saharan Africa: a systematic review and meta analysis. Infect Agent Cancer. 2013;8(1):27. Epub 2013/07/12. pmid:23842471; PubMed Central PMCID: PMCPMC3712022.
  7. 7. Low AJ, Clayton T, Konate I, Nagot N, Ouedraogo A, Huet C, et al. Genital warts and infection with human immunodeficiency virus in high-risk women in Burkina Faso: a longitudinal study. BMC Infect Dis. 2011;11(20):1471–2334.
  8. 8. Silverberg MJ, Ahdieh L, Munoz N, Anastos K, Burk RD, Cu-Uvin S. The Impact of HIV Infection and Immunodeficiency on Human Papillomavirus Type 6 or 11 Infection and on Genital Warts. Sex Transm Dis. 2002;29:427–35. pmid:12172526
  9. 9. Kavanaugh BE, Odem-Davis K, Jaoko W, Estambale B, Kiarie JN, Masese LN, et al. Prevalence and Correlates of Genital Warts in Kenyan Female Sex Workers. Sexually Transmitted Diseases. 2012;39(11):902–5. 00007435-201211000-00016. pmid:23060082
  10. 10. Massad LS, Silverberg MJ, Springer G, Minkoff H, Hessol N, Palefsky JM, et al. Effect of antiretroviral therapy on the incidence of genital warts and vulvar neoplasia among women with the human immunodeficiency virus. Am J Obstet Gynecol. 2004;190(5):1241–8. pmid:15167825
  11. 11. Dolev JC, Maurer T, Springer G, Glesby MJ, Minkoff H, Connell C, et al. Incidence and risk factors for verrucae in women. Aids. 2008;22. pmid:18525267
  12. 12. WHO. Summary of the WHO Position Paper on Vaccines against Human Papillomavirus (HPV)2017 2 October 2017. Available from: http://www.who.int/immunization/policy/position_papers/pp_hpv_may2017_summary.pdf?ua=1.
  13. 13. Joura EA, Giuliano AR, Iversen O-E, Bouchard C, Mao C, Mehlsen J, et al. A 9-Valent HPV Vaccine against Infection and Intraepithelial Neoplasia in Women. New England Journal of Medicine. 2015;372(8):711–23. pmid:25693011.
  14. 14. Garland SM, Kjaer SK, Munoz N, Block SL, Brown DR, DiNubile MJ, et al. Impact and Effectiveness of the Quadrivalent Human Papillomavirus Vaccine: A Systematic Review of 10 Years of Real-world Experience. Clin Infect Dis. 2016;63(4):519–27. Epub 2016/05/28. pmid:27230391; PubMed Central PMCID: PMCPMC4967609.
  15. 15. LaMontagne DS, Bloem PJN, Brotherton JML, Gallagher KE, Badiane O, Ndiaye C. Progress in HPV vaccination in low- and lower-middle-income countries. Int J Gynaecol Obstet. 2017;138 Suppl 1:7–14. pmid:28691329.
  16. 16. Drolet M, Bénard É, Boily M-C, Ali H, Baandrup L, Bauer H, et al. Population-level impact and herd effects following human papillomavirus vaccination programmes: a systematic review and meta-analysis. The Lancet Infectious Diseases. 2015;15(5):565–80. pmid:25744474
  17. 17. Kelly HA, Sawadogo B, Chikandiwa A, Segondy M, Gilham C, Lompo O, et al. Epidemiology of High-risk Human Papillomavirus and Cervical Lesions in African women living with HIV/AIDS: Effect of Anti-Retroviral Therapy. Aids. 2017;31(2):273–85. Epub 2016/10/19. pmid:27755107.
  18. 18. Department of Health RoSA. National consolidated guidelines for the prevention of mother-to-child transmission of HIV (PMTCT) and the management of HIV in children, adolescents and adults http://www.health.gov.za/2014
  19. 19. VIH/SIDA Mdlscmdlcl. Normes et protocoles de prise en charge medicale des personnes vivant avec le VIH au Burkina Faso http://www.remed.org/recommandations_officielles_du_ministere_de_la_sante_nov2008.pdf2008.
  20. 20. UK-NEQAS. Immune Monitoring Programme (Accredited) 2016. Available from: http://www.ukneqasli.co.uk/eqa-pt-programme-information/flow-cytometry-programmes/immune-monitoring/.
  21. 21. QCMD. Quality Control for Molecular Diagnostics 2016. Available from: http://www.qcmd.org/index.php?pageId=3&pageVersion=EN.
  22. 22. Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F, et al. A review of human carcinogens—Part B: biological agents. Lancet Oncol. 2009;10(4):321–2. Epub 2009/04/08. pmid:19350698.
  23. 23. Xue X, Gange SJ, Zhong Y, Burk RD, Minkoff H, Massad LS, et al. Marginal and mixed-effects models in the analysis of human papillomavirus natural history data. Cancer Epidemiol Biomarkers Prev. 2010;19(1):159–69. Epub 2010/01/09. pmid:20056635; PubMed Central PMCID: PMCPMC2839537.
  24. 24. Bursac Z, Gauss CH, Williams DK, Hosmer DW. Purposeful selection of variables in logistic regression. Source code for biology and medicine. 2008;3:17. Epub 2008/12/18. pmid:19087314; PubMed Central PMCID: PMCPMC2633005.
  25. 25. Mbulawa ZZ, Marais DJ, Johnson LF, Boulle A, Coetzee D, Williamson AL. Influence of human immunodeficiency virus and CD4 count on the prevalence of human papillomavirus in heterosexual couples. J Gen Virol. 2010;91(Pt 12):3023–31. Epub 2010/08/20. pmid:20719990.
  26. 26. Djigma FW, Ouedraogo C, Karou DS, Sagna T, Bisseye C, Zeba M, et al. Prevalence and genotype characterization of human papillomaviruses among HIV-seropositive in Ouagadougou, Burkina Faso. Acta tropica. 2011;117(3):202–6. Epub 2010/12/21. pmid:21167118.
  27. 27. Dunne EF, Nielson CM, Stone KM, Markowitz LE, Giuliano AR. Prevalence of HPV infection among men: A systematic review of the literature. J Infect Dis. 2006;194(8):1044–57. Epub 2006/09/23. pmid:16991079.
  28. 28. Mzarico E, Gomez-Roig MD, Guirado L, Lorente N, Gonzalez-Bosquet E. Relationship between smoking, HPV infection, and risk of Cervical cancer. European journal of gynaecological oncology. 2015;36(6):677–80. Epub 2016/01/19. pmid:26775350.
  29. 29. Ghanem KG, Datta SD, Unger ER, Hagensee M, Shlay JC, Kerndt P, et al. The association of current hormonal contraceptive use with type-specific HPV detection. Sex Transm Infect. 2011;87(5):385–8. pmid:21772042; PubMed Central PMCID: PMCPMC3512558.