https://orcid.org/0000-0002-0246-4429
Figures
Abstract
Objectives
This study aimed to: (1) Estimate HPV prevalence and genotype distribution among female sex workers (FSWs) in Mali and Benin as well as the prevalence of multiple HPV type infections in this group, and (2) Identify potential risk factors associated with high-risk (HR) HPV infections.
Methods
We analyzed baseline data of 665 FSWs aged ≥ 18 years recruited during a prospective cohort of cervical cancer screening in Cotonou (Benin) and Bamako (Mali) from 2017 to 2018. The Linear Array HPV genotyping test was used to identify HPV genotypes. Descriptive statistics and multivariate log-binomial regression were used. Adjusted prevalence ratios (APR) with 95% confidence intervals (95%CI) were estimated to identify risk factors associated with HR-HPV infections.
Results
HPV data were available for 659 FSWs (Benin: 309; Mali: 350). The mean age was 35.0 years (± 10.7) in Benin and 26.8 years (± 7.6) in Mali. The overall HPV prevalence rates were 95.5% in Benin and 81.4% in Mali. About 87.7% and 63.4% of FSWs harbored ≥ 2 HPV types in Benin and Mali, respectively. The top three prevalent HR-HPV among FSWs in Benin were: HPV58 (37.5%), HPV16 (36.6%) and HPV52 (28.8%). Corresponding patterns in Mali were HPV16 (15.7%), HPV51 (14.3%) and HPV52 (12.9%). In Benin, the main factors associated with HR-HPV were vaginal douching (APR = 1.17; 95%CI:1.02–1.34) and gonococcal infection (APR = 1.16; 95%CI:1.04–1.28), while in Mali they were sex work duration ≤ 1 year (APR = 1.35; 95%CI:1.10–1.65) and HIV infection (APR = 1.26; 95%CI: 1.06–1.51).
Citation: Tounkara FK, Téguété I, Guédou FA, Goma-Matsétsé E, Koné A, Béhanzin L, et al. (2020) Human papillomavirus genotype distribution and factors associated among female sex workers in West Africa. PLoS ONE 15(11): e0242711. https://doi.org/10.1371/journal.pone.0242711
Editor: Magdalena Grce, Rudjer Boskovic Institute, CROATIA
Received: June 5, 2020; Accepted: November 9, 2020; Published: November 25, 2020
Copyright: © 2020 Tounkara et al. 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 author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: This work was supported by a grant from the Canadian Institutes of Health Research (grant # FDN-143218, Dr. Michel Alary) which was a principal funder. The Réseau FRQS SIDA-MI, Dr François Coutlée, finance a little part of HPV genotyping analysis. FKT was recipient of a PhD studentship from the “Fonds de la recherche en santé du Québec (FRQS)” (Grant # 35546). The funders had no role in study design, data collection and analysis.
Competing interests: MA declares a grant from the Canadian Institutes of Health Research used to fund the present study (grant # FDN-143218) and grants from the Bill & Melinda Gates Foundation and the Public Health Agency of Canada, not related to this work. All other authors have declared that no competing interests exist.
Introduction
Cervical cancer is the second most frequent cancer among women and the leading cause of cancer-related death in many countries in Sub-Saharan Africa (SSA) [1]. It is well established that early detection via screening and effective precancer treatment can significantly reduce the incidence of invasive cervical cancer as well as cancer-related death [2, 3]. However, due to limited access to infrastructure, lack of financial and technical resources and well-trained physicians, effective organized cervical cancer screening programs do not exist in many SSA’s countries [4–6]. As a consequence, this part of the world bears the highest burden of cervical cancer in terms of incidence and mortality [7].
Human papillomavirus (HPV) infection is the most common sexually transmitted infection (STI) worldwide [8]. According to 2007 data, the overall prevalence of HPV infection among sexually active women with normal cytology is estimated at 10.4% worldwide. The lowest prevalence rates are observed in Asia (8%) and the highest in Africa (22.1%) [9]. To date, more than 100 HPV genotypes have been identified. According to their oncogenic potential, they are divided into high risk (HR), possible or probable high risk (pHR) and low risk (LR) HPV types [10]. Most of cervical HPV infections are transient and resolve spontaneously. However, persistent infections with HR-HPV types are the etiological agents of cervical pre-malignant and malignant lesions [11]. HPV16 and HPV18 are the most oncogenic genotypes responsible for 70% of all cases of cervical cancer worldwide, while other HPV types like HPV31, HPV33, HPV35, HPV45, HPV52, and HPV58 represent an additional 20% of cervical cancer cases [12].
Benin and Mali are two West African countries where cervical cancer is a major public health issue [13]. Recent data in 2018 reported age-standardized incidence and mortality rates of 23.7 and 20.2 per 100,000 in Benin, while in Mali they were estimated at 43.9 and 36.2 per 100,000, respectively. Thus, in West Africa, Benin and Mali ranked 13th and 3rd, respectively, in terms of age-standardized incidence and 12th and 3rd rank according to mortality of cervical cancer [13]. Furthermore, the prevalence of HPV infections in women from the general population with normal cytology is estimated at 26.7% in Benin (2011) [14] and 12% in Mali (2011) [15], but little is known about HPV genotype distribution among these women.
In regions of the world where the HIV epidemic is driven by heterosexual transmission, FSWs and their clients constitute a core group for the spread of HIV/STI to the general population [16]. In addition, because of their high level of sexual activity, FSWs are simultaneously at high risk of contracting HPV infections as well as other STIs, including HIV infection [17]. Like for HIV infection, their clients may spread HPV towards women of the general population after acquiring it from FSWs [18]. A literature review conducted in 2013 reported an overall prevalence of HPV infections of 45.7% (range 2.3% to 100%) among FSWs worldwide [17]. Studies conducted in East and West Africa reported an HPV prevalence of 57.7% in Kenya in 2017 [19]; 26% in Ghana in 2019 [20]; 45.2% in Togo in 2019 [21]; 51.1% in Cote d’Ivoire in 2018 [22]; 66,1% in Burkina Faso in 2006 [23] and 79.8% in Senegal in 2019 [24]. Also, according to the studies conducted in 2017 in Kenya and 2019 in Burkina Faso, high prevalence of multiple HPV type infections has been reported in this group, significantly increasing their risk of cervical pre-cancer and cancer [19, 24]. Understanding determinants of cervical HPV infections in FSWs is essential for developing effective prevention programs. Moreover, with the advent of HPV vaccination programs in African countries, it is essential to better understand the distribution of circulating HPV genotypes in this group. To our knowledge, there is no data published concerning HPV infection among FSWs in Benin and Mali. Therefore, this study aimed to: (1) Estimate HPV prevalence and genotype distribution among FSWs in Mali and Benin as well as the prevalence of multiple HPV type infections in this group, and (2) Identify risk factors associated with HR-HPV infections among FSWs in each country.
Materials and methods
Study design and settings
We analyzed baseline data from a prospective cohort study on cervical cancer screening, HPV and HIV infections. From January 2017 to March 2018, the study was conducted in collaboration with three Non-Governmental Organizations (NGOs) in Bamako (Mali) namely SOUTOURA, DANAYA SO and ARCAD-SIDA. These NGOs are the most active in the field of STI/HIV prevention in Mali, with a recognized leadership and trusted by FSWs. The details of NGOs from Mali are shown in a previous publication [25]. In Benin, we collaborated with “Santé et Développement (SED)”, a NGO working with FSWs for over 20 years with confirmed proficiency in STI/HIV prevention and “SOLIDARITÉ”, the local association of women practicing sex work, which operates a network of peer-educators (PEs).
Study planning, recruitment procedures and study population
For uniform procedures in both countries, we have carefully planned our study. Before starting the activities, we organized a training session targeting physicians at the “Dispensaire des IST (DIST)” in Cotonou and at the ARCAD-SIDA STI clinic in Bamako. These sessions covered theorical and practical aspects of cervical cancer screening, review of questionnaire and practical issues related to data collection procedures in the field.
For field work issues, we recruited four PEs from each participating NGO in Mali, two PEs as well as two facilitators in Benin. Specific training was provided to these PEs and facilitators to familiarize them with the survey procedures. A cervical cancer screening awareness campaign preceded the recruitment of FSWs in the two cities. These activities were performed by trained PEs and senior staff from the participating NGOs in the bars, brothels, hotels, etc. Women interested to participate in the study were invited to come to the DIST in Cotonou (STI clinic offering adapted services to FSWs) or to the ARCAD-SIDA STI clinic specialized in FSW care in Bamako. The service package offered in these two FSW-friendly STI clinics includes small talk discussion for behavior change, peer education, condom use demonstration, as well as free distribution of condoms and lubricants. In addition, they provide STI care using syndromic management as well as HIV testing and treatment for free of charge.
Inclusion and exclusion criteria
A FSW was defined as any woman who receives money or gifts in exchange for sex. To be included in the study participants had to fulfill the three criteria as follows: (a) being a FSW in the study settings for at least 6 months; (b) being referred by one of the PEs of the participating NGOs; and (c) being aged between 18 and 65 years old. All FSWs who had previously been diagnosed with cervical cancer, who had a hysterectomy and those who were pregnant were excluded from the study.
Data collection
At the clinics, the study procedures were explained to each FSW and informed consent was obtained. A questionnaire was administered face-to-face by qualified interviewers to collect data on socio-demographic characteristics, including age, marital status, educational level, country of origin; behavioral and sex work characteristics like the frequency of alcohol consumption, drug consumption, number of cigarette packs smoked per week, age at first sexual intercourse, total number of sexual partners in the last week of work, number of paying clients in the last week of work, having regular partner, frequency of condom use, duration of sex work, place where practicing sex work, oral and anal sexual intercourse, vaginal douching; reproductive health variables as well as family and medical histories, such as parity, number of abortions, number of vaginal and cesarian deliveries, history of cervical cancer, ovarian cancer, vaginal cancer or breast cancer, family history of these cancers; history of STIs in the last six months.
Following the interview, each woman underwent a gynecological examination performed by a well-trained physician. Vaginal and cervical swabs were obtained from each participant for curable STI testing. A cytobrush (Rovers Medical Devices B.V. Oss, The Netherlands) was used to collect cells from the endocervix and ectocervix and stored in tubes containing 3 ml of fresh PBS (phosphate-buffered saline) for HPV DNA testing. Cervical cancer screening was performed using visual inspection methods (VIA/VILI). Finally, a venous blood sample was taken for HIV testing and CD4 count.
Laboratory analyses
In both countries, the same procedures were used for the following tests. Wet mounts of the vaginal swabs were microscopically examined immediately for Trichomonas vaginalis and Candida Albicans. The diagnosis of bacterial vaginosis (BV) was made using Nugent score.
For Neisseria gonorrhoeae and Chlamydia trachomatis detection, we used different tests. In Benin, these infections were detected using the NG/CT Probetec® assay from Becton Dickenson (Cockeysville, MD, USA), while in Mali, these infections were detected using the Abbott Real-Time CT/NG assay. These assays were carried out according to the manufacturers’ instructions. These two tests have similar sensitivity and specificity as reported elsewhere [26, 27]. All FSWs with laboratory diagnosed STIs received appropriate treatment, free of charge.
Both, Mali and Benin have the same national HIV testing algorithms and we thus used the same tests for the detection of HIV antibodies for all the study participants. Alere Determine HIV-1/2 test (Alere Medical Co. Ltd) was used as first line and positive specimens were then confirmed with a rapid and discriminatory test SD Bioline (Giheung-gu, Yongin-si, Korea). Women testing positive for HIV were immediately offered antiretroviral therapy.
Amplification and detection of HPV DNA
HPV testing was performed as part of the study. In both countries, endocervical and ectocervical cells were centrifuged at 3000 rpm for 5 minutes at 4°C for cervical cell concentration. Cell pellets were stored at -20°C until DNA extraction. At the end of the study, all cell pellets were sent to Montreal (Canada) in a laboratory specialized in HPV testing. HPV Genotyping was performed using the Linear Array HPV genotyping test (Roche Molecular Systems, Inc., Laval, Qc, Canada). This test identifies 36 genotypes by hybridization on a linear array of PGMY-generated amplicons with 34 type-specific probes for HPV-6, -11, -16, -18, -26, -31, -33, -35, -39, -40, -42, -45, -51, -53, -54, -55, -56, -58, -59, -61, -62, -64, -66, -67, -68, -69, -70, -71, -72, -73, -81, -83, -84, and -89, two probes for two variants of HPV-82, and one probe that cross-reacts with HPV-33, -35, -52, and -58. Amplification and type detection were performed according to the manufacturer’s recommendations [28]. The cross-reactivity between HPV 52 and HPV-33, 35 or 58 was further tested with a Real-Time PCR assay specific for HPV52 as described elsewhere [29].
Outcome variable
The dependent variable was “HR-HPV” defined as being positive for at least one HR-HPV type.
Statistical analysis
Data were analyzed using SAS version 9.4 (SAS institute, Inc, Cary, North Carolina, USA). We used descriptive statistics to analyze demographic, behavioral and sex work characteristics. The results were presented as percentages, means with standard deviations (± SD) or medians with inter-quartile ranges (IQR). Pearson’s Chi-Square test or Fisher’s exact test were used for categorical variables, the Student’s t test for continuous variables and rank test with median score. We also computed descriptive statistics to estimate HPV prevalence rates. An exact binomial 95% confidence interval (95%CI) was calculated for each prevalence rate. To identify factors associated with HR-HPV infection, we carried out, for each country, univariate and multivariate log-binomial regression models with a robust “sandwich” variance estimator to calculate the adjusted prevalence ratios (APRs) with 95% CI. All variables significant at p ≤ 0.2 (as recommended for variable selection [30]) in univariate analysis or known from the literature as potential risk factors were computed into multivariate log-binomial regression models to build the full model. Manual backwards elimination procedures were applied to remove covariates from the full model if they were neither significant (p-value ≥ 0.05) or changing the effect estimate for the association of HR-HPV and other variables by more than 10%.
Ethical issues
The study was reviewed and approved by the ethics committees of the School of medicine of Bamako, Mali (#no2017/93/CE/FMPOS), and the “CHU de Québec-Université Laval” (#2017–3313), as well as the National Ethics Committee for Health Research in Benin (# no01 du 25 janvier 2017). The objectives, procedures and potential risks related to the study were explained to each woman and written informed consent was obtained before enrolment. Consenting participants signed or apposed their fingerprint on the consent forms. In Mali, participants received 5000 CFA (about US$8.40), while in Benin, they received 3000 CFA (about US$5) for compensation of transportation and the time spent at the clinic. The amount of compensation was slightly different between the two countries because of the socio-economic context (the cost of living is slightly higher in Mali than in Benin). Condoms and lubricants were distributed for free to each woman.
Results
Socio-demographic and sexual behavioral characteristics
A total of 710 FSWs were approached (337 in Benin and 373 in Mali). Of these 20 FSWs were excluded in Mali (one previous cervical cancer and 19 pregnancies) and 25 FSWs in Benin (9 had their menstruations and 16 pregnancies). A total of 665 women were thus included in the study, 312 in Benin and 353 in Mali. FSWs from Benin were older than those from Mali, the mean (± SD) age was 35.0 ± 10.7 versus 26.8 ± 7.6 years respectively (p < .0001; Table 1). FSWs in Mali were mostly never married (69%), whereas they were mostly separated/widowed/divorced in Benin (55.8%). In both countries nearly 40% were uneducated. Mean age at first sexual intercourse and mean age at first paid sex were lower in Mali (p < .0001). Vaginal douching was more frequently practiced in Benin (79.2%) as compared with Mali FSWs (40.3%, p < .0001). Curable STIs (C. trachomatis, N. gonorrhoeae, T. vaginalis) were more common in Mali than in Benin (p < 0.0001). HIV prevalence was 26.3% in Benin versus 20.4% in Mali and the proportions of positive VIA/VILI tests were 20.2% and 10.5%, respectively.
HPV prevalence and type-specific distribution
Among 665 women enrolled, HPV data were not available for 6 FSWs, three in Mali and three in Benin, because of invalid cervical samples, but all other variables were available for these FSWs. The prevalence rates of any HPV among FSWs were 95.5% (95%CI: 92.5–97.5) and 81.4% (95%CI: 77.0–85.4) in Benin and Mali, respectively (Table 2). HR-HPV prevalence was higher among FSWs from Benin than those from Mali (87.1% versus 62.3%; p < .0001). A similar pattern was found with LR-HPV types (77.4% versus 55.4%, respectively; p < .0001). HPV prevalence was 87.5% in HIV-positive FSWs vs 86.9% among HIV-negative FSWs in Benin. These statistics in Mali were 71.8% and 59.9% for HIV-positive and HIV-negative FSWs, respectively.
HPV genotype distribution varied widely between the two countries with HPV58 (37.5%), HPV16 (36.6%) and HPV52 (28.8%) being the top three in Benin, compared to HPV16 (15.7%), HPV51 (14.3%) and HPV52 (12.9%) in Mali. Concerning LR-HPV the top three in Benin were: HPV62 (35.6%), HPV81 (23.6%) and HPV61 (23.0%). This profile in Mali was HPV62 (15.1%), HPV61 (12.9%) and HPV84 (10.6%).
The prevalence of HR-HPV multiple type infections (≥ 2 HR-HPV) among FSWs was high in this study (Table 3). Multiple HR-HPV type infections occurred in 61.8% and 35.4% of FSWs in Benin and Mali respectively (p < 0.0001). A similar profile was observed for LR-HPV types.
To evaluate the potential effectiveness of the HPV vaccines in FSWs, we estimated the overall prevalence rate and multiple infections with HPV types covered by available vaccines (Table 3). This analysis was restricted to HR-HPV positive FSWs. In Benin, about 65.1% of HR-HPV positive FSWs had at least one of the 4 genotypes covered by the Gardasil-4 vaccine versus 43.1% in Mali (p < 0.0001). Regarding the Gardasil-9 vaccine, 91.2% and 76.6% of HR-HPV infected FSWs harbored at least one HPV type prevented by the 9-valent HPV vaccine in Benin and Mali, respectively (p < 0.0001). Finally, none of the study participants had previously received any HPV vaccine at the time of the study.
Risk factors for HR-HPV infections
The country-specific risk factors for HR-HPV types are shown in Table 4. In Benin, the main risk factors for HR-HPV infections were age < 25 years (APR = 1.25; 95%CI: 1.06–1.47), p-trend test = 0.015, vaginal douching before and after sex (APR = 1.17; 95%CI: 1.02–1.34) and gonorrhea infection (APR = 1.16; 95%CI: 1.04–1.28). In Mali, these figures were being home-based FSWs (APR = 2.13; 95%CI: 1.17–3.86), sex work duration ≤ 1 year (APR = 1.35 95%CI: 1.10–1.65) and being HIV positive (APR = 1.26 95%CI: 1.06–1.51). There was an inverse association between high number of clients (≥ 5) compared to a lower number (< 5; p < 0.001).
Discussion
This is the first epidemiological study on HPV infection among FSWs in Benin and Mali. Our findings show a high overall prevalence rate of HPV infection as well as a high rate of HR-HPV types among this group in both countries. Similarly, there was a high rate of multiple HPV type infections among them. We also noted a wide variation of HPV genotype distribution among FSWs between both countries. Furthermore, we identified several factors associated HR-HPV infections that were different between countries.
Prevalence rates of HPV infection of 95.5% in Benin and 81.4% in Mali were several-fold higher than those found among women in the general population (26.7% in Benin and 12% in Mali according to 2011 data) [14, 15]. These results corroborate previous reports [31] and are partly explained by the cumulative exposure of FSWs to high number of sexual partners and other STIs.
Furthermore, compared with other studies among FSWs in West Africa where HR-HPV prevalence rates varied from 32.9% - 72.5% [21, 24], the prevalence rate of HR-HPV infections was higher in our study. There was also a high prevalence of multiple infections (≥ 2 HR-HPV) as reported elsewhere among FSWs in SSA settings [17, 19, 24]. Two hypotheses may explain these differences. First, HIV prevalence among FSWs participating in other HPV studies in West Africa (10.6% to 15.4) [21, 24] is lower than what we found (26.3% among FSWs in Benin and 20.4% in Mali). The strong association between HIV and HPV is well documented and a higher HIV prevalence could thus lead to a higher HPV prevalence [32]. Secondly, the HPV genotyping methods and the number of HPV types detected differ between studies cited here (21 HPV to 28 HPV) [21, 24]. We used a more sensitive assay, detecting up to 36 HPV types.
A striking feature of HPV infection in FSWs is the wide variation of genotypes across countries. HPV58 and HPV16 were the most frequent in Benin, while in Mali, the predominant genotypes were HPV16 and HPV51. These two prevalent genotypes are quite different from those observed in Madagascar [33], but shared the presence of one HR-HPV (HPV16) with observations in Ghana, Senegal and Kenya [19, 20, 24]. Such findings confirm the epidemiological particularity of circulating HPV types in Sub-Saharan Africa [9].
Like elsewhere [20, 21], in Benin, the main potential risk factors of HR-HPV were younger age, gonorrhea as well as vaginal douching. Multiple sexual partners and immature cervix producing inadequate cervical mucus are known factors increasing the likelihood of acquiring HPV infection in young women [34, 35]. Similarly, the associations between HR-HPV and could partly be explained by risky behavior or because of the chronic inflammation caused by this infection, which would facilitate the acquisition of HPV [36, 37]. On the other hand, the relationship between vaginal douching and higher rates of genital infections, including STIs/HPV, is well documented in the general population [38, 39]. Supporting these findings are local immune system disturbance and removal of the cervical mucus protective barrier [40]. In contrast, a study in Cambodia found less HPV infections in FSWs practicing vaginal douching just after sexual intercourse [41]. Such controversies, combined with the rarity of studies exploring the link between HPV infection and vaginal douching among FSWs, despite its common practice among this population in SSA [42], suggest the need for additional studies.
In Mali the highest prevalence of HR-HPV was observed in home-based FSWs. While the relation between HIV infection and sex work place is well documented [43], little is known about the link between HPV infection and sex work place. Differences in number and social categories of sexual partners, prevalence of other STIs, cumulative exposure according to place of sex work could potentially explain our findings. As reported elsewhere [44], shorter sex work duration association with increased HPV prevalence shares a mechanism similar to that of younger age. Our analysis revealed a moderate association between HIV and HR-HPV infections. Supporting this finding is the HIV-induced immunosuppression that can increase the susceptibility to virus acquisition as well as the inability to eliminate HPV infection [45, 46]. Unexpectedly, we noted an inverse association between the high number of sexual partners and HPV infections. This is probably due to the high HIV prevalence among FSWs with < 5 clients as observed in our data publish elsewhere [25].
Limitations
We studied the epidemiology of HPV infections among FSWs from a convenience sample. To minimize selection bias, we used different recruitment approaches either by the PEs or FSW leaders. In addition, self-reported information on risky behaviors such as number of sexual partners, age at first sexual interaction, etc., may be subject to recall and social desirability bias. Indeed, due to stigma or social desirability, this information collected may be underestimated. However, such misclassification bias would be limited in our study since the data were collected by trusted qualified interviewers. Failure to measure certain variables could cause residual confounding in our analyses on risk factors. To deal with this bias, we adjusted for a large number of potential confounding variables based on the literature. Furthermore, the associations we found between several factors and HR-HPV cannot be viewed as causal because of the cross-sectional nature of the data. In addition, although our results reflect current HPV infection status, the lack of information about HPV antibodies in the serum, which is reported to be a marker of past infections [47], constitute another limitation of this study. Finally, our results may not be generalizable to other FSW populations in West Africa because of the different characteristics observed in FSWs.
Implications
Understanding the distribution of HPV genotypes among FSWs is crucial to estimating the burden of HPV infections as well as the future impact of HPV vaccines within this high-risk group. The high prevalence of HR-HPV and other STIs/HIV found among FSWs implies an increased risk of cervical cancer.
Currently, the bivalent and quadrivalent vaccines are approved in both countries, but they are not integrated in the routine immunization program. However, these vaccines are available in private pharmacies. The relative low representativeness of HPV types covered by the quadrivalent vaccine like HPV6, HPV11 and HPV18 may question the potential effectiveness of this vaccine in this group. The presence of HPV16, HPV52 and HPV58 in the top three would favor the nonavalent vaccine for cervical cancer prevention. Since evidences support vaccination against HPV in sexually active women up to the age of 44 years [48], FSW immunization program may be one of the most promising strategies to decrease the HPV burden. However, socio-cultural and financial barriers in most SSA countries may impact the implementation and the effectiveness of programs, at least on the short term. The high rates of positive VIA/VILI found in our study call for an emphasis on secondary prevention, namely cervical cancer screening using HPV molecular testing [3]. An integration of these services to the package of services, as recommended by international professional societies, would be a valuable strategy.
Conclusion
Our study, the first in Benin and Mali, showed a high prevalence of HR-HPV infections among FSWs. These results make FSW a priority group for cervical cancer prevention programs adapted to their context. Additionally, the identification of the predominant HPV genotypes in this population shed the light on the vaccine potentially needed in this specific group.
Acknowledgments
The authors acknowledge the contribution all the staff involved in this project, in particular the staff of the DIST in Benin, staff of the ARCAD-SIDA, SOUTOURA and DANAYA SO in Mali. We thank also the staff of all laboratories involved in the project such as “Unité de Biochimie et de Biologie Moléculaire (UBBM) in Benin, ALGI, “Centre universitaire de recherche clinique, Bamako”, and the pathology unit of Point G Teaching Hospital in Mali, as well as the staff of the gynecology units of the CHU-MEL Teaching Hospital in Benin and the Gabriel Toure Teaching Hospital in Mali. Finally, we thank all women who participated in the study.
References
- 1. Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424. pmid:30207593
- 2. Jansen EEL, Zielonke N, Gini A, et al. Effect of organised cervical cancer screening on cervical cancer mortality in Europe: a systematic review. Eur J Cancer. 2020. pmid:31980322
- 3. Canfell K, Kim JJ, Brisson M, et al. Mortality impact of achieving WHO cervical cancer elimination targets: a comparative modelling analysis in 78 low-income and lower-middle-income countries. Lancet. 2020;395:591–603. pmid:32007142
- 4. Sankaranarayanan R, Budukh AM, Rajkumar R. Effective screening programmes for cervical cancer in low- and middle-income developing countries. Bull World Health Organ. 2001;79:954–62. pmid:11693978
- 5. Gustafsson L, Ponten J, Zack M, et al. International incidence rates of invasive cervical cancer after introduction of cytological screening. CCC. 1997;8:755–63. pmid:9328198
- 6. Nygard JF, Skare GB, Thoresen SO. The cervical cancer screening programme in Norway, 1992–2000: changes in Pap smear coverage and incidence of cervical cancer. J Med Screen. 2002;9:86–91. pmid:12133929
- 7. Parkin DM, Bray F. Chapter 2: The burden of HPV-related cancers. Vaccine. 2006;24 Suppl 3:S3/11–25. pmid:16949997
- 8. Bekkers RL, Massuger LF, Bulten J, et al. Epidemiological and clinical aspects of human papillomavirus detection in the prevention of cervical cancer. Rev Med Virol. 2004;14:95–105. pmid:15027002
- 9. de Sanjose S, Diaz M, Castellsague X, et al. Worldwide prevalence and genotype distribution of cervical human papillomavirus DNA in women with normal cytology: a meta-analysis. Lancet Infect Dis. 2007;7:453–9. pmid:17597569
- 10. Munoz N, Bosch FX, de Sanjose S, et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003;348:518–27. pmid:12571259
- 11. Bosch FX, Burchell AN, Schiffman M, et al. Epidemiology and natural history of human papillomavirus infections and type-specific implications in cervical neoplasia. Vaccine. 2008;26 Suppl 10:K1–16. pmid:18847553
- 12. Clifford G, Franceschi S, Diaz M, et al. Chapter 3: HPV type-distribution in women with and without cervical neoplastic diseases. Vaccine. 2006;24 Suppl 3:S3/26–34. pmid:16950015
- 13.
Bruni L, Albero G, Serrano B, et al. ICO/IARC Information Centre on HPV and Cancer (HPV Information Centre). Human Papillomavirus and Related Diseases in the World. Summary Report 17 June 2019. [Accessed, October 08th, 2019].
- 14. Piras F, Piga M, De Montis A, et al. Prevalence of human papillomavirus infection in women in Benin, West Africa. Virol J. 2011;8:514. pmid:22074103
- 15. Tracy JK, Traore CB, Bakarou K, et al. Risk factors for high-risk human papillomavirus infection in unscreened Malian women. Trop Med Int Health. 2011;16:1432–8. pmid:21749583
- 16. Lowndes CM, Alary M, Meda H, et al. Role of core and bridging groups in the transmission dynamics of HIV and STIs in Cotonou, Benin, West Africa. Sex Transm Infect. 2002;78 Suppl 1:i69–77.
- 17. Soohoo M, Blas M, Byraiah G, et al. Cervical HPV Infection in Female Sex Workers: A Global Perspective. Open AIDS J. 2013;7:58–66. pmid:24511334
- 18. Ghosh I, Ghosh P, Bharti AC, et al. Prevalence of human papillomavirus and co-existent sexually transmitted infections among female sex workers, men having sex with men and injectable drug abusers from eastern India. Asian Pac J Cancer Prev. 2012;13:799–802. pmid:22631651
- 19. Menon S, van den Broeck D, Rossi R, et al. Multiple HPV infections in female sex workers in Western Kenya: implications for prophylactic vaccines within this sub population. Infect Agent Cancer. 2017;12:2. pmid:28070215
- 20. Adams AR, Nortey PA, Dortey BA, et al. Cervical Human Papillomavirus Prevalence, Genotypes, and Associated Risk Factors among Female Sex Workers in Greater Accra, Ghana. J Oncol. 2019;2019:1–7.
- 21. Ferre VM, Ekouevi DK, Gbeasor-Komlanvi FA, et al. Prevalence of human papillomavirus, human immunodeficiency virus and other sexually transmitted infections among female sex workers in Togo: a national cross-sectional survey. Clin Microbiol Infect. 2019;25:1560.e1–.e7. pmid:31051265
- 22. Ouattara A, Yeo A, Blavo-Kouame EB, et al. Humans Papillomavirus (Hpv) Infections in Female Sex Workers in Cote D’ivoire. Am J Cancer Res. 2017;1:1–12.
- 23. Didelot-Rousseau MN, Nagot N, Costes-Martineau V, et al. Human papillomavirus genotype distribution and cervical squamous intraepithelial lesions among high-risk women with and without HIV-1 infection in Burkina Faso. Br J Cancer. 2006;95:355–62. pmid:16832413
- 24. Diop-Ndiaye H, Beiter K, Gheit T, et al. Human Papillomavirus infection in senegalese female sex workers. Papillomavirus Res. 2019;7:97–101. pmid:30771492
- 25. Tounkara FK, Téguété I, Guédou F, et al. Prevalence and Factors Associated with HIV and Sexually Transmitted Infections among Female Sex Workers in Bamako, Mali. Sex Transm Dis. 2020;47:679–85. pmid:32932403
- 26. Fontana C, Favaro M, Cicchetti O, et al. Performance of strand displacement amplification assay in the detection of Chlamydia trachomatis and Neisseria gonorrhoeae. Jpn J Infect Dis. 2005;58:283–8. pmid:16249622
- 27.
Abbott RealTime CT/NG[package insert]. Des Plaines, IL: Abbott Molecular; 2007.
- 28. Coutlee F, Rouleau D, Petignat P, et al. Enhanced detection and typing of human papillomavirus (HPV) DNA in anogenital samples with PGMY primers and the Linear array HPV genotyping test. J Clin Microbiol. 2006;44:1998–2006. pmid:16757590
- 29. Coutlee F, Rouleau D, Ghattas G, et al. Confirmatory real-time PCR assay for human papillomavirus (HPV) type 52 infection in anogenital specimens screened for HPV infection with the linear array HPV genotyping test. J Clin Microbiol. 2007;45:3821–3. pmid:17898159
- 30. Steyerberg EW, Eijkemans MJ, Harrell FE Jr., et al. Prognostic modelling with logistic regression analysis: a comparison of selection and estimation methods in small data sets. Stat Med. 2000;19:1059–79. pmid:10790680
- 31. Gonzalez C, Torres M, Canals J, et al. Higher incidence and persistence of high-risk human papillomavirus infection in female sex workers compared with women attending family planning. Int J Infect Dis. 2011;15:e688–94. pmid:21757383
- 32. Clifford GM, Goncalves MA, Franceschi S. Human papillomavirus types among women infected with HIV: a meta-analysis. AIDS. 2006;20:2337–44. pmid:17117020
- 33. Smith JS, Van Damme K, Randrianjafisamindrakotroka N, et al. Human papillomavirus and cervical neoplasia among female sex workers in Madagascar. Int J Gynecol Cancer. 2010;20:1593–6. pmid:21370602
- 34. Kahn JA, Rosenthal SL, Succop PA, et al. The interval between menarche and age of first sexual intercourse as a risk factor for subsequent HPV infection in adolescent and young adult women. J Pediatr. 2002;141:718–23. pmid:12410205
- 35. McSorley J. Non-HIV sexually transmitted infections. Obstet Gynaecol Reprod Med. 2013;23:180–84.
- 36. de Abreu AL, Malaguti N, Souza RP, et al. Association of human papillomavirus, Neisseria gonorrhoeae and Chlamydia trachomatis co-infections on the risk of high-grade squamous intraepithelial cervical lesion. Am J Cancer Res. 2016;6:1371–83. pmid:27429850
- 37. Lenz JD, Dillard JP. Pathogenesis of Neisseria gonorrhoeae and the Host Defense in Ascending Infections of Human Fallopian Tube. Front Immunol. 2018;9:2710. pmid:30524442
- 38. Sun CA, Hsiung CA, Lai CH, et al. Epidemiologic correlates of cervical human papillomavirus prevalence in women with abnormal Pap smear tests: a Taiwan Cooperative Oncology Group (TCOG) study. J Med Virol. 2005;77:273–81. pmid:16121376
- 39. Moscicki AB, Ma Y, Farhat S, et al. Redetection of cervical human papillomavirus type 16 (HPV16) in women with a history of HPV16. J Infect Dis. 2013;208:403–12. pmid:23599313
- 40. Chu TY, Chang YC, Ding DC. Cervicovaginal secretions protect from human papillomavirus infection: effects of vaginal douching. Taiwan J Obstet Gynecol. 2013;52:241–5. pmid:23915858
- 41. Bui TC, Scheurer ME, Pham VTT, et al. Intravaginal practices and genital human papillomavirus infection among female sex workers in Cambodia. J Med Virol. 2018;90:1765–74. pmid:30016541
- 42. Myer L, Kuhn L, Stein ZA, et al. Intravaginal practices, bacterial vaginosis, and women's susceptibility to HIV infection: epidemiological evidence and biological mechanisms. Lancet Infect Dis. 2005;5:786–94. pmid:16310150
- 43. Pitpitan EV, Kalichman SC, Eaton LA, et al. HIV/STI risk among venue-based female sex workers across the globe: a look back and the way forward. Curr HIV/AIDS Rep. 2013;10:65–78. pmid:23160840
- 44. Sarkar K, Bhattacharya S, Bhattacharyya S, et al. Oncogenic human papilloma virus and cervical pre-cancerous lesions in brothel-based sex workers in India. J Infect Public Health. 2008;1:121–8. pmid:20701853
- 45. Hawes SE, Critchlow CW, Faye Niang MA, et al. Increased risk of high-grade cervical squamous intraepithelial lesions and invasive cervical cancer among African women with human immunodeficiency virus type 1 and 2 infections. J Infect Dis. 2003;188:555–63. pmid:12898443
- 46. Ahdieh L, Klein RS, Burk R, 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:682–90. pmid:11517428
- 47. Wentzensen N, Rodriguez AC, Viscidi R, et al. A competitive serological assay shows naturally acquired immunity to human papillomavirus infections in the Guanacaste Natural History Study. J Infect Dis. 2011;204:94–102. pmid:21628663
- 48. Wheeler CM, Skinner SR, Del Rosario-Raymundo MR, et al. Efficacy, safety, and immunogenicity of the human papillomavirus 16/18 AS04-adjuvanted vaccine in women older than 25 years: 7-year follow-up of the phase 3, double-blind, randomised controlled VIVIANE study. Lancet Infect Dis. 2016;16:1154–68. pmid:27373900