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

Associations of Human Papillomavirus (HPV) genotypes with high-grade cervical neoplasia (CIN2+) in a cohort of women living with HIV in Burkina Faso and South Africa

  • Helen A. Kelly ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    helen.kelly@lshtm.ac.uk

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

  • Jean Ngou ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliation INSERM U1058 and University Hospital (CHRU), Montpellier, France

  • Admire Chikandiwa ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

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

  • Bernard Sawadogo ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliation Centre de Recherches Internationales en Santé, University of Ouagadougou, Ouagadougou, Burkina Faso

  • Clare Gilham ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

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

  • Tanvier Omar ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliation National Health Laboratory Services, Johannesburg, South Africa

  • Olga Lompo ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliation Centre de Recherches Internationales en Santé, University of Ouagadougou, Ouagadougou, Burkina Faso

  • Sylviane Doutre ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliation INSERM U1058 and University Hospital (CHRU), Montpellier, France

  • Nicolas Meda ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliation Centre de Recherches Internationales en Santé, University of Ouagadougou, Ouagadougou, Burkina Faso

  • Helen A. Weiss ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliations Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, London, United Kingdom, MRC Tropical Epidemiology Group, London School of Hygiene and Tropical Medicine, London, United Kingdom

  • Sinead Delany-Moretlwe ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliation National Health Laboratory Services, Johannesburg, South Africa

  • Michel Segondy ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

    Affiliation INSERM U1058 and University Hospital (CHRU), Montpellier, France

  • Philippe Mayaud ,

    Contributed equally to this work with: Helen A. Kelly, Jean Ngou, Admire Chikandiwa, Bernard Sawadogo, Clare Gilham, Tanvier Omar, Olga Lompo, Sylviane Doutre, Nicolas Meda, Helen A. Weiss, Sinead Delany-Moretlwe, Michel Segondy, Philippe Mayaud

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

  • for the HARP Study Group

    The complete membership of the author group can be found in the Acknowledgments section.

Abstract

Objective

To describe associations of high-risk human papillomavirus (HR-HPV) with high-grade cervical intraepithelial neoplasia (CIN2+) in women living with HIV (WLHIV) in Burkina Faso (BF) and South Africa (SA).

Methods

Prospective cohort of WLHIV attending HIV outpatient clinics and treatment centres. Recruitment was stratified by ART status. Cervical HPV genotyping using INNO-LiPA and histological assessment of 4-quadrant cervical biopsies at enrolment and 16 months later.

Results

Among women with CIN2+ at baseline, the prevalence of any HR-HPV genotypes included in the bi/quadrivalent (HPV16/18) or nonavalent (HPV16/18/31/35/45/52/58) HPV vaccines ranged from 37% to 90%. HPV58 was most strongly associated with CIN2+ (aOR = 5.40, 95%CI: 2.77–10.53). At 16-months follow-up, persistence of any HR-HPV was strongly associated with incident CIN2+ (aOR = 7.90, 95%CI: 3.11–20.07), as was persistence of HPV16/18 (aOR = 5.25, 95%CI: 2.14–12.91) and the additional HR types in the nonavalent vaccine (aOR = 3.23, 95%CI: 1.23–8.54).

Conclusion

HR-HPV persistence is very common among African WLHIV and is linked to incident CIN2+. HPV vaccines could prevent between 37–90% of CIN2+ among African WLHIV.

Introduction

Women living with HIV (WLHIV) have a higher prevalence of genital high-risk human papillomavirus (HR-HPV) infection than the general population [1], are more likely to be infected with multiple HR types [2, 3] and have greater persistence of infection [4] and risk of cervical intraepithelial neoplasia (CIN) progression [5]. WLHIV have been shown to be more commonly infected with types other than HPV16 or 18 [6] and their high-grade cytological lesions are frequently attributed to types other than HPV16/18 [3].

As WLHIV are living longer due to increased availability of antiretroviral therapy (ART), many retain a high risk of infection with HPV and high risk of progressing to cervical and other genital cancers. Primary prevention of HPV infection through vaccination could reduce the burden of infection and disease and on screening and treatment services. Current bivalent and quadrivalent HPV vaccines target two HR types (HPV16 and 18) responsible for about 70% of cervical cancers [7], whereas the nonavalent vaccine which protects against a wider range of HR-HPV types (HPV16/18/31/33/45/52/58), is estimated to prevent up to 90% of cervical cancers in women from the general population [7, 8]. These potential benefits have seldom been estimated among WLHIV, particularly in sub-Saharan Africa. An improved understanding of HPV type distribution associated with histological lesions in this population is needed to guide HPV vaccine programme decisions.

We conducted a large prospective study of cervical cancer screening in a cohort of women living with HIV-1 in Burkina Faso and South Africa (HARP—HPV in Africa Research Partnership). In this paper, we describe (i) the prevalence, persistence, incidence and genotype distribution of HPV, and (ii) their association with prevalent and incident CIN grade 2 and higher (CIN2+).

Materials and methods

Study population

Participants were recruited from the HIV outpatient clinic of the University Teaching Hospital of Ouagadougou, Burkina Faso (BF), and HIV treatment centres from inner city Johannesburg South Africa (SA), from December 2011 to October 2012. Inclusion criteria were being HIV-1 seropositive, aged 25–50 years and resident in the city. Exclusion criteria were history of prior treatment for cervical cancer, previous hysterectomy, and being pregnant or less than 8 weeks postpartum. Enrolment was stratified in a 2:1 ratio of ART-users:ART-naïve. Written informed consent was obtained at the screening visit when eligibility for the study was assessed. Participants were followed-up every 6 months for CD4+ T-lymphocytes cell count and up to Month 18 (endline) when procedures similar to baseline were repeated.

Specimen collection

A venous blood sample was collected to confirm HIV-1 serostatus if necessary, and to obtain CD4+ cell counts. Cervical samples were collected by a nurse-midwife including a Digene cervical sampler at enrolment (baseline) and careHPV cervical sampler at endline (both Qiagen, Gaithesburg, MD) for HPV-DNA testing and genotyping. All women were screened for abnormal cervical lesions using cervical smear for cytology by Papanicolaou staining and visual inspection using acetic acid and Lugol’s iodine (VIA/VILI), and referred for colposcopy performed by trained colposcopists. Systematic 4-quadrant cervical biopsy, including directed biopsy of any suspicious lesions, was performed for participants who had abnormalities detected on cytology, VIA/VILI or colposcopy, or who were HR-HPV DNA positive. The same genital sampling and examination procedures were repeated at a scheduled follow-up visit at Month 18.

Laboratory testing

HIV infection was diagnosed according to national guidelines for each country [9, 10]. Measurement of CD4+ cell count was performed using FACScount (Becton-Dickinson, NJ) in both countries. HPV-DNA genotyping of all samples was performed at the virology laboratory of the Gui de Chauliac Teaching Hospital, University of Montpellier, France, using the INNO-LiPA HPV genotyping Extra® assay (Fujirebio, Courtaboeuf, France) with external quality assurance [11]. HR-HPV types were defined using the current International Agency for Research on Cancer (IARC) classification [12] to include ‘carcinogenic to humans’ (HPV16/18/31/33/35/39/45/51/52/56/58/59) and ‘probable carcinogenic’ (HPV68). The ‘possible carcinogenic’ types (HPV26/53/66/69/70/73/82) and other known low-risk types (HPV6/11/40/43/44/54/71/74) were considered as low-risk (LR-HPV).

Cervical biopsies were processed at the local pathology laboratories in Ouagaoudou and Johannesburg and read using the 3-tier CIN classification system [13]. Histology was classified as “negative” (≤CIN1) or “positive” (CIN2+) based on the highest reading across all findings from the 4-quadrant biopsies and endocervical curettage, if indicated. All histological slides from women with a local diagnosis of CIN2+ and approximately 10% of slides from women with normal or CIN1 histological findings were reviewed by the HARP Endpoint Committee of five pathologists for final classification, after baselineand endline rounds.

Statistical analysis

HR-HPV genotype-specific persistence was defined as being positive for the same HR-HPV type at baseline and endline visits. Any type clearance was defined as being positive for a specific HR type at baseline and negative for that type at endline visit. Cumulative incidence of any HR-HPV was defined as the proportion of women who were negative for a specific type at baseline and positive for that type at endline visit.

For associations of prevalent and incident CIN2+ with individual HR-HPV types, logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CI). Multivariable analyses adjusted for site and the socio-demographic and behavioural factors which were independently associated in univariate analyses (p<0.10) with CIN2+ for each country, as described in [14], were conducted. Data were analysed using Stata version 14 (Stata Statistical Software, College Station. TX: Stata Corporation).

This study received ethical approval from 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 in the UK (no. 7400).

Results

Study population

A full description of study participants has been published elsewhere [14]. Of 1479 women screened, 1238 were enrolled (BF: 615; SA: 623). The median age of participants was 36 (interquartile range [IQR], 31–42) years in BF and 34 (IQR, 30–40) years in SA. About half (49.7%) of SA participants had ever had a Pap smear, and a fifth (20.8%) of BF participants had ever had a VIA/VILI examination, the primary cervical cancer screening modality in each country respectively. At enrolment, 422 (68.6%) participants were on ART in BF and 406 (65.2%) in SA. In BF, the median CD4+ count was 417 (IQR, 315–606) cells/mm3 among ART-naive participants and 446 (IQR, 309–600) cells/mm3 among ART users. In SA, the median CD4+ count was 448 (IQR, 353–614) cells/mm3 among ART-naive and 420 (IQR, 279–567) cells/mm3 among ART users.

Prevalence of HPV and association with CIN2+ at enrolment

Of the 1238 participants enrolled, 1215 (98.1%) had valid HPV genotyping results (BF: 96.6%; SA: 99.7%) and 998 (82.1%) were positive for any HPV (BF: 75.3%; SA: 88.7%). The proportion of multiple HR-HPV genotypes was lower in BF (BF: 41.9% vs. SA: 55.2%; p<0.001; Table 1). HPV52 was the most prevalent type in both countries (BF: 20.4%; SA: 24.2%). HPV16/18 were detected in 15.0% of participants in BF vs. 29.0% in SA (p<0.001). Among those without HPV16/18 co-infection, the proportion infected by any of the additional HR types included in the nonavalent (9v) vaccine (HPV31/33/45/52/58) was 30.1% in BF and 36.7% in SA (p<0.001). The prevalence of any of the non-vaccine HR types (HPV35/39/51/56/59/68) was 35.5% in BF and 49.1% in SA (p<0.001) irrespective of co-infection with any of the vaccine types, and was 14.0% in BF and 13.4% in SA among those not co-infected by any vaccine types.

thumbnail
Table 1. HR-HPV infection at baseline and endline follow-up among Women Living with HIV (WLHIV) in Burkina Faso (BF) and South Africa (SA).

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

Among women with both histology and HPV results (BF: 546; SA: 573), the prevalence of HR-HPV was higher among those with higher CIN grades in both settings (BF: 48.4% in CIN normal, 71.8% in CIN1, 100% in CIN2 and 100% in CIN3+, p<0.001; SA: 72.4% in CIN normal, 81.4% in CIN1, 89.5% in CIN2 and 92.5% in CIN3+, p<0.001;) (Fig 1(A) and 1(B)). The prevalence of any HR-HPV genotypes included in the bi/quadrivalent and nonavalent HPV vaccines was 45.2% and 90.3% among women with CIN2+ in BF, respectively, and 37.2% and 79.8% in SA, respectively. In the absence of any vaccine types, the prevalence of any non-vaccine types was 9.7% (3/31) in BF and 10.9% (14/129) in SA.

thumbnail
Fig 1.

(A) HR-HPV prevalence by CIN grade among 546 women living with HIVin Burkina Faso. (B) HR-HPV prevalence by CIN grade among 573 women living with HIV in South Africa.

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

In both countries, HPV58 was the genotype most strongly associated with CIN2 compared to ≤CIN1 cases (adjusted Odds Ratio [aOR] = 5.40, 95%CI: 2.77–10.53), while types most strongly associated with CIN3+ were HPV16 in BF (aOR = 34.56, 95%CI: 5.70–209.49) and HPV58 in SA (aOR = 3.65, 95%CI: 1.42–9.37) (Fig 2(A) and 2(B) and S1 Table). The prevalence of HPV58 was 22.1% and 21.1% among women with CIN2 in BF and SA respectively, and 23.1% and 15.1% among women with CIN3+.

thumbnail
Fig 2.

(A) Association of HR-HPV prevalence with prevalent CIN2 and CIN3+ among 546 women living with HIV in Burkina Faso. (B) Association of HR-HPV prevalence with prevalent CIN2 and CIN3+ among 573 women living with HIV in South Africa.

https://doi.org/10.1371/journal.pone.0174117.g002

The prevalence of CIN2+ was significantly higher among those positive for HPV16/18 compared to those who were negative for both types in BF (17.7% vs. 3.6%; aOR = 7.90, 95%CI: 3.23–19.33) and in SA (29.1% vs. 19.9%, aOR = 1.49, 95%CI: 0.95–2.34), although the latter was not statistically significant (S1 Table). Among 875 women without co-infection of either HPV16 or 18, the prevalence of CIN2+ was higher among women positive for any of the additional HR-HPV types targeted by the nonavalent vaccine (31/33/45/52/58) compared to women negative for all of these types in both countries (BF: 8.8% vs. 1.0%; SA: 25.2% vs. 13.7%; aOR = 2.78, 95%CI: 1.67–4.64 for both countries combined). The prevalence of CIN2+ was similar among women who were positive or negative for any of the non-vaccine HR-HPV types. However, among 498 women without infection with vaccine types (BF: 308, SA: 190), the prevalence of CIN2+ was higher among women positive for any of the non-vaccine types compared to women who were HR-HPV negative (11.1% vs. 3.4%; aOR = 3.35, 95%CI: 1.36–8.26).

Incidence and persistence of HR-HPV

Of the 963 women without CIN2+ at enrolment who returned at follow-up M18 visit (median of actual follow-up 16 months [IQR, 15.6–16.8]), genotyping data at both visits was available for 922 (95.7%) women. The proportion of women with an incident HR-HPV infection was similar in both sites (BF: 47.9% [228/476] vs. SA: 49.3% [220/446]; p = 0.67). Type-specific incident infection was highest for HPV52 in both countries (BF: 12.7%; SA: 16.6%). Incident infection of HPV16/18 was observed in 14.6% of women in BF and 14.7% in SA. In the absence of HPV16/18 co-infection, incident infection of any of the additional HR-HPV types of the nonavalent vaccine was observed in 19.7% of women in BF and 23.9% in SA; incident infection of non-vaccine HR-HPV types in the absence of vaccine type co-infection was observed in 14.0% in BF and 13.5% in SA (Table 1).

Persistence was highest for HPV58 (72.2%) in BF and HPV35 (42.6%) in SA. Forty-nine women (18.2%) in BF and 80 (23.5%) in SA had at least one persistent and one cleared infection during the 16 months follow-up period (Table 1).

Association of incident CIN2+ with HR-HPV persistence

Among 922 women with genotyping data at both visits, concomitant histology results were available for 780 (84.6%) women (BF: 405; SA: 375). Cumulative incidence of CIN2+ was higher among those with HR-HPV persistence compared to those who cleared or who were HR-HPV negative at baseline in both settings (BF: 3.5% vs. 0.3%; SA: 13.4% vs. 2.0%; both countries combined aOR = 7.90, 95%CI: 3.11–20.07, Table 2).

thumbnail
Table 2. Risk of incident CIN2+ according to HR-HPV infection status over 16 months among 405 Women Living with HIV (WLHIV) without CIN2+ at enrolment in Burkina Faso and 375 WLHIV in South Africa.

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

CIN2+ incidence was found to be significantly associated with the persistence of any members of the alpha-9 HPV (HPV16-related) and any alpha-7 HPV (HPV-18 related) genotypes in both countries (Table 2). In particular, CIN2+ incidence was higher among women with persistent HPV16/18 compared to those who cleared either or who were negative for both at baseline (BF: 7.7% vs. 0.8%; SA: 18.2% vs. 4.7%; both sites combined: aOR = 5.25, 95%CI: 2.14–12.91).

Among 721 women without persistent HPV16 or 18, CIN2+ incidence was higher among those with persistence of any of the additional 5 HR-HPV included in the nonavalent (both sites combined: 6.1% vs. 2.0% aOR = 3.23, 95%CI: 1.23–8.54). In the absence of co-persistence by any of the vaccine types (606 women), CIN2+ incidence was higher among women with persistent non-vaccine types compared to those who cleared or who were HR-HPV negative at baseline in both settings (8.1% vs. 1.1%; aOR = 7.87, 95%CI: 2.40–25.81).

Association of HIV-related factors with HR-HPV prevalence at baseline

The prevalence of HR-HPV types was similar among ART users and ART-naive participants at baseline in BF, whilst in SA, the prevalence of HPV16/18 was higher among ART-naïve participants (S2 and S3 Tables). Among ART users in BF, the prevalence of any nonavalent HR-HPV in the absence of HPV16/18 was higher among women with low CD4+ count (<200 cells/mm3) compared to those with high CD4+ count (>500 cells/mm3) and this association was significant among ART users only, however no association was observed for HPV16/18 or the non-vaccine types. In SA, only the prevalence of non-vaccine HR-HPV was higher among ART-taking women with low CD4+ count (S3 Table).

Discussion

Women living with HIV (WLHIV) in Burkina Faso and South Africa have a high prevalence, incidence and persistence of HR-HPV and correspondingly a high prevalence and incidence of cervical neoplasia. The high rates of HR-HPV found in both countries are similar to other studies in Sub Saharan Africa [1, 2, 15, 16]. The increased risk in the prevalence of HR-HPV and CIN2+ among women in South Africa compared to Burkina Faso may be explained by the less well controlled HIV disease (in terms of lower rates of ART adherence and HIV viral suppression) found among women in South Africa, as we reported elsewhere [14], as well as by other cofactors for HR-HPV and CIN2+, such as greater frequency of contraceptive use and smoking and higher prevalence of mucosal STIsobserved in this and other cohort of WLHIV [1725].

While HR-HPV prevalence was higher in South Africa at baseline, HR-HPV persistence was higher among women in Burkina Faso. This may be due to the fact that estimates of HR-HPV persistence were restricted to women without evidence of CIN2+ at baseline, who underwent treatment. As there was a higher prevalence of CIN2+ at baseline in South Africa, their exclusion from the follow-up analysis potentially reduced the overall proportion of women with persistent infection at endline. By contrast, because of the lower prevalence of CIN2+ in Burkina Faso at baseline, fewer women were excluded in the follow-up analysis, which potentially allowed a greater number of persistent HR-HPV cases to accrue during the follow-up period. In addition, we cannot exclude the possibility that the 4-quadrant biopsy may have helped clear HR-HPV from the cervical epithelium, as a greater proportion of women in South Africa underwent a 4-quadrant biopsy.

HPV52 was the type most frequently detected overall in both countries, but HPV58 was the type most strongly associated with cervical lesions, and one of the most likely to persist in both countries. In a meta-analysis describing HPV type distribution among 5578 WLHIV from 20 studies including 5 from African countries, HPV16, 18 and 58 were the most prevalent among those with cytological high-grade squamous intraepithelial lesions (HSIL+) [3, 26]. Additionally, WLHIV with HSIL+ were less likely to harbour HPV16 compared to those with HSIL in the general population, and more likely to be infected with HPV18/33/51/52/58, suggesting that the HPV16/18 based bi/quadrivalent vaccines may have a lower impact on cervical cancer rates in settings where HIV prevalence is high [27]. The high prevalence of types other than HPV16 and HPV18 among WLHIV in this study is similar to what has been reported elsewhere [6]. While we found that up to 45% of CIN2+ in BF and 37% of CIN2+ in SA could be prevented by the HPV16/18-targeting vaccines, up to 90% of CIN2+ in BF and 80% of CIN2+ in SA could be prevented by the nonavalent vaccine, as predicted from studies in the general population [8, 28].

We found that up to 45% of CIN2/3 in BF and 37% of CIN2/3 in SA could be prevented by the HPV16/18-targeting vaccines, and that a further 45% of CIN2/3 in BF and 43% of CIN2/3 in SA could be prevented by the additional 5 types (HPV31/33/45/51/52) contained in the nonavalent vaccine. However, not all CIN2/3 develop into invasive cervical cancer (ICC), and not all HPV types found in CIN2 have the same propensity to evolve towards ICC. Clifford et al, in their systematic review comparing the HPV type distribution in ICC biopsy and cervical specimens of 770 HIV-positive women from 21 studies in 12 African countries [29], report a higher pooled prevalence of HPV16/18 (61.7%) than found in this study, whilst the prevalence of the other nonavalent vaccine types was 11.2%, lower than what we report in our study. These findings indicate that the relative contribution of HPV16/18 increases with increasing lesion grade, whilst the contribution of the additional nonavalent vaccine still represents a significant preventable fraction of ICC.

Our study also shows that the persistence of any of the bi/quadrivalent and nonavalent HR-HPV types is strongly associated with the incidence of cervical lesions among WLHIV, as is the case for HIV uninfected women [30]. Moreover, our study found a significant proportion of women had an incident infection with any of the nonavalent vaccine types (35%) as well as incident infection of multiple HR-HPV types over the 16-month follow-up period. This emphasises the potential impact of vaccination in this population using an HPV vaccine with a broader-range protection potential, given that these types are also likely to persist in this population [31].

A smaller but non-negligible proportion (11%) of women with CIN2+ were exclusively infected by non-vaccine types at enrolment. In addition, the persistence of non-vaccine types in the absence of current vaccine types, was also strongly associated with incident CIN2+, and this was largely driven by persistence of HPV35 and 56. A large retrospective cross-sectional study using histologically conformed ICC collected from 38 countries worldwide reported that 91% of cases were positive for any of the current nonavalent vaccine types in addition to HPV35 [7]. Until further generation vaccines incorporate a wider range of HPV types, cervical screening will remain important among WLHIV where non-vaccine types are prevalent.

ART-naïve women were at increased risk of HPV16 or 18 infection in South Africa at baseline. Among ART users, a decrease in CD4+ count was also associated with an increased risk of any of the nonavalent other than HPV16/18 and the non-vaccine type prevalence in Burkina Faso, and of the non-vaccine types in South Africa. However, there was no reduction in HPV16 or 18 prevalence with increasing CD4+ counts in either country. This is consistent with previous reports suggesting that HPV16 might be more weakly associated with immune suppression than other HR types due to its potential to evade the host immune response, thus explaining its greater contribution to high-grade lesions and cancer in immunocompetent women [32]. However, these associations may also be reflective of the overall increased risk of any HR-HPV prevalence among those with low CD4+, as we have previously reported elsewhere [14].

Given that the bi/quadrivalent and nonavalent types are the most important contributors to cervical lesions in this population, maintaining a stable high CD4+ count is necessary to reduce the risks of HR-HPV persistence and cervical lesion development. This should be achieved by early initiation of ART before CD4+ counts decline and maintaining good HIV-1 virological control thereafter. In addition, HPV vaccination could be offered to these women before CD4+ counts decline.

Limitations of this study include the absence of HPV genotyping from the biopsies, which preclude accurate attribution of CIN to particular genotypes. Furthermore, our definition of cumulative HR-HPV incidence over 16 months is at the same time crude and of limited duration to assess the actual roles of HR-HPV incidence and persistence on CIN2+ incidence. Other studies have used shorter time intervals for defining HR-HPV incidence [33] because of the transient nature of HR-HPV infection. Furthermore, the study could not rule out type-specific clearance and reinfection when estimating persistence during the 16-months interval between HPV testing. The evaluation of HPV at two time points only did not allow precise estimation of the duration of infections. Longer duration of follow-up would have allowed to accrue a larger number of incident CIN2+ cases and more robustly assess the role of incident HR-HPV on CIN development over sufficient follow-up time. A key strength of this study, however, is the availability of a large proportion of women with verified histological endpoints and concurrent genotyping data. In addition, the HARP study was undertaken in two African countries with very different HIV epidemics and possibly different burdens of HPV infection and cervical cancer, but the similarity of data allows the findings to be extended to a range of countries and settings in sub-Saharan Africa.

In conclusion, we confirm that HR-HPV infection and cervical lesions are very common among WLHIV in Africa, and that a broader range of genotypes are potentially associated with CIN2+ development. While currently available bi/quadrivalent vaccines could prevent up to 45% of treatable precursor lesions, the nonavalent vaccine has the potential to prevent up to 90% of cases in WLHIV. HPV vaccination could reduce the incidence of HR-HPV related disease among WLHIV in addition to contributing cost saving to current screening and treatment programmes, but this would require a formal health economic evaluation.

Supporting information

S1 Table. Association of HPV type prevalence with prevalent CIN2 and CIN3+ among 546 women living with HIV in Burkina Faso and 573 in South Africa.

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

(DOCX)

S2 Table. Association of HR-HPV type prevalence with ART and CD4+ count at enrolment among 570 women living with HIV in Burkina Faso.

https://doi.org/10.1371/journal.pone.0174117.s002

(DOCX)

S3 Table. Association of HR-HPV type prevalence with ART and CD4+ count at enrolment among 613 women living with HIV in South Africa.

https://doi.org/10.1371/journal.pone.0174117.s003

(DOCX)

Acknowledgments

Contributing members of the HARP study group included: A. Chikandiwa, E. Cutler, S. Delany-Moretlwe, 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, N. Meda, B. Sawadogo, J. Simporé, A. Yonli, 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); and the members of the HARP Endpoint Committee (EPC) including T. Omar (Johannesburg, South Africa), O. Goumbri-Lompo (Ouagadougou, Burkina Faso), S. Doutre, N. Leventoux, V. Costes (Montpellier, France) and O. Clavero (Institut Catala d’Oncologia, Barcelona, Spain).

Author Contributions

  1. Conceptualization: PM SDM NM HAW MS HK.
  2. Data curation: CG HK HAW.
  3. Formal analysis: HK CG HAW.
  4. Funding acquisition: PM.
  5. Investigation: HK JN AC BS TO OL SD.
  6. Methodology: PM HK SDM NM MS.
  7. Project administration: HK AC BS.
  8. Resources: HK JN AC BS TO OL SD CG.
  9. Software: CG.
  10. Supervision: PM SDM NM HK.
  11. Validation: TO OL SD CG HAW HK.
  12. Visualization: HK PM.
  13. Writing – original draft: HK.
  14. Writing – review & editing: HK JN AC BS CG TO OL SD NM HAW SDM MS PM.

References

  1. 1. McDonald AC, Tergas AI, Kuhn L, Denny L, Wright TC Jr. Distribution of Human Papillomavirus Genotypes among HIV-Positive and HIV-Negative Women in Cape Town, South Africa. Front Oncol. 2014;4:48. pmid:24672770
  2. 2. De Vuyst H, Mugo NR, Chung MH, McKenzie KP, Nyongesa-Malava E, Tenet V, et al. Prevalence and determinants of human papillomavirus infection and cervical lesions in HIV-positive women in Kenya. Br J Cancer. 2012;107(9):1624–30. pmid:23033006
  3. 3. Clifford GM, Goncalves MA, Franceschi S. Human papillomavirus types among women infected with HIV: a meta-analysis. AIDS. 2006;20(18):2337–44. pmid:17117020
  4. 4. Ahdieh L, Klein RS, Burk R, Cu-Uvin S, Schuman P, Duerr A, et al. 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–90. pmid:11517428
  5. 5. Denslow SA, Rositch AF, Firnhaber C, Ting J, Smith JS. Incidence and progression of cervical lesions in women with HIV: a systematic global review. Int J STD AIDS. 2014;25(3):163–77. pmid:24216030
  6. 6. Konopnicki D, Manigart Y, Gilles C, Barlow P, De Marchin J, Feoli F, et al. High-risk human papillomavirus genotypes distribution in a cohort of HIV-positive women living in Europe: epidemiological implication for vaccination against human papillomavirus. Aids. 2016;30(3):425–33. pmid:26765936
  7. 7. de Sanjose S, Quint WGV, Alemany L, Geraets DT, Klaustermeier JE, Lloveras B, et al. Human papillomavirus genotype attribution in invasive cervical cancer: a retrospective cross-sectional worldwide study. The Lancet Oncology. 2010;11(11):1048–56. pmid:20952254
  8. 8. Joura EA, Ault KA, Bosch FX, Brown D, Cuzick J, Ferris D, et al. Attribution of 12 high-risk human papillomavirus genotypes to infection and cervical disease. Cancer Epidemiol Biomarkers Prev. 2014;23(10):1997–2008. pmid:25274978
  9. 9. 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
  10. 10. 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.
  11. 11. Ngou J, Gilham C, Omar T, Goumbri-Lompo O, Doutre S, Michelow P, et al. Comparison of analytical and clinical performances of the digene HC2 HPV DNA assay and the INNO-LiPA HPV genotyping assay for detecting high-risk HPV infection and cervical neoplasia among HIV-positive African women. J Acquir Immune Defic Syndr. 2015;68(2):162–8. pmid:25394189
  12. 12. Bouvard V, Baan R., Straif K., Grosse Y., Secretan B., El Ghissassi F., Benbrahim-Tallaa L., Guha N., Freeman C., Galichet L., Cogliano V. on behalf of the WHO International Agency for Research on Cancer Monograph Working Group. A review of human carcinogens—Part B: biological agents. Lancet Oncology. 2009;10:321–2. pmid:19350698
  13. 13. Clement B, Young R. Atlas of Gynecologic Pathology: Chapter 5 –Invasive squamous carcinoma of the cervix and its precursors. 3rd ed: Elsevier Inc; 2014. p. 100–19.
  14. 14. 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. 2016.
  15. 15. Firnhaber C, Van Le H, Pettifor A, Schulze D, Michelow P, Sanne IM, et al. Association between cervical dysplasia and human papillomavirus in HIV seropositive women from Johannesburg South Africa. Cancer Causes Control. 2010;21(3):433–43. pmid:19949850
  16. 16. Didelot-Rousseau MN, Nagot N, Costes-Martineau V, Valles X, Ouedraogo A, Konate I, et al. Human papillomavirus genotype distribution and cervical squamous intraepithelial lesions among high-risk women with and without HIV-1 infection in Burkina Faso. BrJCancer. 2006;95(3):355–62.
  17. 17. Giuliano AR, Sedjo RL, Roe DJ, Harri R, Baldwi S, Papenfuss MR, et al. Clearance of oncogenic human papillomavirus (HPV) infection: effect of smoking (United States). Cancer Causes Control. 2002;13(9):839–46. pmid:12462549
  18. 18. Lazenby GB, Taylor PT, Badman BS, McHaki E, Korte JE, Soper DE, et al. An association between Trichomonas vaginalis and high-risk human papillomavirus in rural Tanzanian women undergoing cervical cancer screening. Clin Ther. 2014;36(1):38–45. pmid:24417784
  19. 19. Gillet E, Meys JF, Verstraelen H, Bosire C, De Sutter P, Temmerman M, et al. Bacterial vaginosis is associated with uterine cervical human papillomavirus infection: a meta-analysis. BMC Infect Dis. 2011;11:10. pmid:21223574
  20. 20. Castellsague X, Munoz N. Chapter 3: Cofactors in human papillomavirus carcinogenesis—role of parity, oral contraceptives, and tobacco smoking. J Natl Cancer Inst Monogr. 2003(31):20–8. pmid:12807941
  21. 21. Smith JS, Green J, Berrington de Gonzalez A, Appleby P, Peto J, Plummer M, et al. Cervical cancer and use of hormonal contraceptives: a systematic review. Lancet. 2003;361(9364):1159–67. pmid:12686037
  22. 22. Zhu H, Shen Z, Luo H, Zhang W, Zhu X. Chlamydia Trachomatis Infection-Associated Risk of Cervical Cancer: A Meta-Analysis. Medicine (Baltimore). 2016;95(13):e3077.
  23. 23. Smith JS, Munoz N, Herrero R, Eluf-Neto J, Ngelangel C, Franceschi S, et al. Evidence for Chlamydia trachomatis as a human papillomavirus cofactor in the etiology of invasive cervical cancer in Brazil and the Philippines. J Infect Dis. 2002;185(3):324–31. pmid:11807714
  24. 24. Zhang ZF, Begg CB. Is Trichomonas vaginalis a cause of cervical neoplasia? Results from a combined analysis of 24 studies. Int J Epidemiol. 1994;23(4):682–90. pmid:8002180
  25. 25. Smith JS, Herrero R, Bosetti C, Munoz N, Bosch FX, Eluf-Neto J, et al. Herpes simplex virus-2 as a human papillomavirus cofactor in the etiology of invasive cervical cancer. J Natl Cancer Inst. 2002;94(21):1604–13. pmid:12419786
  26. 26. Clifford GM, Smith JS, Aguado T, Franceschi S. Comparison of HPV type distribution in high-grade cervical lesions and cervical cancer: a meta-analysis. Br J Cancer. 2003;89(1):101–5. pmid:12838308
  27. 27. Denny LA, Franceschi S, de Sanjose S, Heard I, Moscicki AB, Palefsky J. Human papillomavirus, human immunodeficiency virus and immunosuppression. Vaccine. 2012;30 Suppl 5:F168–74.
  28. 28. Clifford G, Franceschi S, Diaz M, Munoz N, Villa LL. Chapter 3: HPV type-distribution in women with and without cervical neoplastic diseases. Vaccine. 2006;24 Suppl 3:S3/26–34.
  29. 29. Clifford GM, de Vuyst H, Tenet V, Plummer M, Tully S, Franceschi S. Effect of HIV Infection on Human Papillomavirus Types Causing Invasive Cervical Cancer in Africa. J Acquir Immune Defic Syndr. 2016;73(3):332–9. pmid:27331659
  30. 30. Trottier H, Mahmud SM, Lindsay L, Jenkins D, Quint W, Wieting SL, et al. Persistence of an incident human papillomavirus infection and timing of cervical lesions in previously unexposed young women. Cancer Epidemiol Biomarkers Prev. 2009;18(3):854–62. pmid:19223559
  31. 31. Blitz S, Baxter J, Raboud J, Walmsley S, Rachlis A, Smaill F, et al. Evaluation of HIV and highly active antiretroviral therapy on the natural history of human papillomavirus infection and cervical cytopathologic findings in HIV-positive and high-risk HIV-negative women. J Infect Dis. 2013;208(3):454–62. pmid:23624362
  32. 32. Strickler HD, Palefsky JM, Shah KV, Anastos K, Klein RS, Minkoff H, et al. Human papillomavirus type 16 and immune status in human immunodeficiency virus-seropositive women. JNatlCancer Inst. 2003;95(14):1062–71.
  33. 33. Houlihan CF, Baisley K, Bravo IG, Kapiga S, de Sanjose S, Changalucha J, et al. The Incidence of Human Papillomavirus in Tanzanian Adolescent Girls Before Reported Sexual Debut. J Adolesc Health. 2016;58(3):295–301. pmid:26725717