Variability of high risk HPV genotypes among HIV infected women in Mwanza, Tanzania- the need for evaluation of current vaccine effectiveness in developing countries

Background

High risk (HR) human papilloma Virus (HPV) genotypes have been associated with cervical cancer. In Tanzania there is a limited data on the epidemiology of HPV and genotypes distribution among HIV infected women. Here we document varieties of HPV genotypes associated with cervical squamous intraepithelial lesions (SIL) among HIV- infected women at Bugando Medical Centre, Mwanza-Tanzania.

Methods

A cross sectional hospital based study involving HIV infected women was conducted between August and October, 2014. Exfoliated cells from ectocervix and endocervix were collected using cytobrush. HPV genotypes were detected using polymerase chain reaction (PCR) followed by sequencing using specific primers targeting broad range of HPV types. Cytology was done to establish squamous intraepithelial lesions. Log binomial regression analysis was done to establish risk ratios (RR) associated with HPV infection using STATA version 11.

Results

A total of 255 HIV infected women with mean age 39.2 ± 9.1 years were enrolled in the study. HPV DNA was detected in 138/255 (54.1 %, 95 % CI: 47-60) of HIV infected women. Twenty six genotypes were detected in various combinations; of these 17(65.3 %) were of HR genotypes. HR genotypes were detected in 124(48.6 %) of HIV infected women. Common HR genotypes detected were HPV-52(26), HPV-58(21), HPV-35(20) and HPV-16(14). The risk of being HPV positive was significantly higher among women with CD4 counts <100 (RR: 1.20, 95 % CI: 1.05-1.35, P = 0.006) and women with SIL (RR: 1.37, 95 % CI: 1.11-1.68, P = 0.005)

Conclusion

Significant proportion of HIV infected women with low CD4 counts have various grades of cervical SIL associated with varieties of uncommon HR genotypes. There is a need to evaluate the effectiveness of the current vaccine in preventing cervical cancer in developing countries where HIV is endemic.

Human papillomavirus (HPV) is a very common sexually transmitted infection which is acquired through body fluids and skin-to-skin contact [1]. It has been known as the causal agent of cervical cancer [2]. In sub-Saharan Africa, the incidence of cervical cancer has reached up to 40.7 cases per 100,000 women and it is the second cause of death among cancers affecting women [3, 4]. In HIV infected women the risk of acquiring HPV is even higher due to immune suppression caused by HIV infection which has the same epidemiological pattern as HPV in sub-Saharan Africa [5, 6]. In most of these cases infection with multiple genotypes has been found to predominate [5, 7–13].

Using molecular techniques, different HPV genotypes have been identified and categorized as high risk (HR) genotypes that are considered to be carcinogenic. These include HPV- 16, 18, 31, 33, 35, 95, 45, 51, 52, 56, 58, 59, 68, 73 and 82; and are considered to cause approximately 95 % of cervical cancers [14]. Human papilloma virus genotypes 6, 11, 40, 42, 43, 44, 54, 61, 70, 72, 81, and cp 6108 have been classified as low risk (LR); majority of these have been found to cause genital warts [15]. Vaccines for HPV are available in some countries and have been found to be effective in reducing infection and its associated consequences [16].

As reported earlier in Africa, women particularly those who are HIV infected are under-represented in the global estimate of HPV genotypes [17]. This necessitates the need to emphasize HPV screening programs for understanding its epidemiology and effective control. Here we document varieties of HPV genotypes associated with cervical squamous intraepithelial lesions (SIL) and associated factors among HIV- infected women at Bugando Medical Centre, Mwanza-Tanzania underscoring the importance of focused control strategies to prevent cervical cancer in developing countries.

Study design, site and population

A cross sectional hospital based study involving 255 HIV infected women was conducted between August and October, 2014 at Bugando Medical Centre (BMC) HIV care and treatment clinic (CTC). The clinic attends patients who are on antiretroviral therapy (ART) and those who have not started ART. The CTC provides routine Papanicolaou (PAP) smear and gynaecological services with no routine screening of HPV. An average of 80 revisiting HIV infected women are seen per day [18].

Patient selection, socio-demographic and clinical information

The study included all HIV infected women aged >18 years attending the CTC during the study period while excluding women who underwent wedge resection of the cervix and those with missed CD4 counts within 3 months of specimen collection. Socio-demographic and relevant clinical information collected were number of sexual partners, parity, marital status, age at 1st sexual debut, history of STI, history of genital warts, contraceptive use, and CD4 counts.

Sample collection and Laboratory procedures

Cervical exfoliated cells from the ectocervix and endocervix were obtained using cytobrush; there after the tip of the cytobrush was placed into a transport medium and stored at -10 °C until processing. Papanicolaou (PAP) smear was also taken, processed and interpreted based on the revised 2001 Bethesda system [19].

Detection of HPV DNA

After vigorous vortexing of transport media, cytobrush was removed and 200 μl of sample was used for nucleic acid extraction with the MagNaPure LC DNA Large Volume Kit on a MagNaPure LC instrument (Roche Diagnostics, Mannheim, Germany) with elution volume set at 100 μl. For the detection of HPV DNA, two different end-point PCR assays (GEN-1 and GEN-2) were done as previously described [20]. Both assays were done using primers targeting sequences which are conserved among a broad range of HPV genotypes, including all major genital HPV genotypes (Tables 1 and 2). These primers amplify fragments of approximately 460 bp for GEN-1 and 220 bp for GEN-2. All PCR products were sequenced (Eurofins Genomics, Ebersberg, Germany); for sequencing of the GEN-1 PCR products, only primer CP4 was used, whereas for GEN-2 PCR products both PPF1 and PPR2 primers were used as previously described [20]. After editing the sequence files using Vector NTI Advance 9 software (Life Technologies, Darmstadt, Germany), sequences were subjected to BLAST search (http://blast.ncbi.nlm.nih.gov/Blast.cgi) to identify the HPV genotypes. The efficiency of nucleic acid extraction was controlled by a real-time PCR assay from our routine diagnostics simultaneously targeting human cytomegalovirus (HCMV) and the human albumin gene.

Data analysis

The data were entered in the computer using Microsoft excel software and analysed using STATA version 11(STATA Corp LP, USA). Continuous variables were summarized as mean with standard deviation or as median with interquartile range whilst categorical variables were summarized as proportions. Socioeconomic status (SES) was defined using education, employment and business status of the participant. Log binomial regression for estimation of relative risk ratios was done followed by log multinomial regression analysis for the factors with P value of <0.2. Cross tabulation was done to detect variables with collinearity. Factors found with collinearity were not subjected into log multinomial regression analysis. The chi square test and fisher’s exact test were done to compare the distribution of HPV genotypes and SIL. Women with mixed genotypes were regarded as having HR genotype if HR genotype was present. Level of significance was measured by p-value whereby factors with p-values of less than 0.05 were considered statistically significant.

Sociodemographic and obstetric characteristics

The mean age of enrolled women was 39.2 ± 9.14 years. A total of 111(43.5 %) were married and 46 (18.03 %) and 209(81.7 %) had higher and lower education respectively. Majority of the enrolled women 138(54.1 %) had high socioeconomic status. A total of 158(61.9 %) were residing in rural areas (Table 3).

Prevalence of HPV

Of the 255 enrolled women, 138 (54.1 %, 95 % CI: 47-60) were found to have HPV DNA in their cervical cells. Out of 138 HPV infected women; 124 (48.6 %) were found to be infected with HR genotypes.

Genotypes distribution

In 138 HPV infected women; 26 genotypes were detected in various combinations. A total of 17/138(12.3 %) of HIV infected women had multiple HR genotypes. Of 26 genotypes, 17(65.4 %; 95%CI: 44.7-81.2) were HR while 9 (34.6 %, 95%CI: 14-49) were LR genotypes (p = 0.028). High risk genotypes detected were 16, 18, 31, 33, 34, 35, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73 and 82 whereas LR genotypes detected were 6, 40, 42, 47, 54, 74, 83 and 90. The total frequency of the genotypes detected was 160 (Table 4); of these common HR genotypes were HPV-52(26, 16.3 %), HPV-58 (21, 13.1 %), HPV-35(20, 12.5 %) and HPV-16(14, 8.8 %).

Cervical squamous intraepithelial lesions (SIL) and genotypes distribution

Out of the 255 women; 91(35.6 %) had SIL on cytological studies. Of 91 HIV infected women with SIL, 74 (81.3 %) had low grade SIL whilst 17(18.7 %) had high grade SIL. Eleven (64.7 %) of women with high grade SIL had positive HPV DNA in their cervical cells. Amongst these eleven women none-16 and non-18 genotypes were detected in 9(81.8 %) of them. Out of 124 women infected with HR genotypes, 55 (44.4 %) had any kind of SIL compared to only 31 (26.3 %) of 118 women with no HPV infection (P = 0.003) Fig. 1. Non-significant association between LR genotypes and SIL (5/14 vs. 31/118, P = 0.264) was observed. Of 17 women with multiple high risk genotypes, 1(5.8 %) had high grade SIL, 9(52.9 %) had low grade SIL and 7(41.2 %) had no SIL.

High risk and low risk HPV genotypes in relation to cervical squamous intraepithelial lesions

Full size image

Factors associated with HPV positivity among HIV infected women

On log binomial regression analysis; women with good socioeconomic status, low CD4 counts, use of ARV ≤6 months and presence of SIL significant high risk of HPV infection as shown in Table 4. Only women with low CD4 counts (RR: 1.2, 95 % CI: 1.05-1.35, P = 0.006) and presence of SIL (RR: 1.37, 95 % CI: 1.11-1.68, P = 0.005) remained with significant higher risk of being infected with HPV on log multinomial regression analysis.

Infection with oncogenic or HR human papillomavirus (HPV) genotypes is the major cause of cervical cancer. With the introduction of the quadrivalent HPV vaccine which prevents four HPV types: HPV 16 and 18, as well as HPV 6 and 11, that cause 90 % of genital warts the knowledge concerning the distribution of HPV genotypes has become increasingly public health concern [21]. This study for the first time in Mwanza, Tanzania; documents high proportion of HIV infected women carrying varieties of HR genotypes similarly to previous reports [6, 22, 23]. As in previous studies [24], a significant proportion of women were uncommonly infected with multiple HR genotypes. These findings are in agreement with previous studies which reported uncommon HR genotypes infecting HIV infected women in developing countries [4, 22, 25, 26]. This calls for the need to include uncommon high risk genotypes in the current vaccine production especially for areas where HIV is endemic. The recently introduced a 9 valent (nanovalent) vaccine which includes HPV 52, 58, 45[27]; the commonest genotypes detected in this study is of great advancement especially in sub Saharan Africa where HIV is endemic. The World health organization (WHO) should consider advocating replacing the current quadrivalent vaccine which is available in many centers in developing countries with nanovalent vaccine. One of the interesting findings in this study is the presence of HPV 35 genotype in relative high frequency which has been also found in other African countries [22, 26, 28]. Unexpectedly, the recent introduced nanovalent vaccine does not include HPV 35 genotype emphasing the need to use epidemiological data for future vaccine development.

Variation of HPV genotypes can occur even within the country which might have implication in vaccine efficacy in the same country. As evidenced by the fact that in the current study, HPV-16 was not the commonest genotype among women with high grade SIL while in the previous study in Tanzania among HIV and non HIV infected women HPV-16 was the commonest genotype among women with high grade SIL [12]. There is a possibility that the current existing HPV vaccines might have diverse efficacy especially in developing countries.

The current study confirms the clear association between HPV infection and SIL. As in previous studies [29–31], it was observed that women infected with HR genotypes were more likely to have SIL than non-infected and those infected with LR genotypes. Though none of the women with LR genotypes had high grade SIL; it was observed that 35.7 % of women infected with LR genotypes had low grade SIL, following these women over time to determine cervical changes and acquisition of other genotypes will provide important information regarding the role of LR genotypes and SIL. In addition, women with low CD4 counts (<100) had significantly higher risk of acquiring HPV infection (low and high risk genotypes) as compared to their counterparts. Our findings are similar to what have been reported earlier [11, 32–34]. This is due to immune suppression which poses the risk of opportunistic infections including HPV.

The major limitation of this study is inability to perform HIV viral load which could be used to correlate virological failure and HPV infection. In addition some important information such as tobacco use and information for those who declined were not collected. However the study established HPV prevalence and genotypes distribution among HIV infected women.

In present study, HIV seropositive women with low CD4 counts and various grades of cervical SIL are significantly infected with HR HPV genotypes. In countries where HIV is endemic the governments should consider addition of HPV vaccine in national immunization programme focusing on epidemiology of HPV genotypes. Data on the efficacy of the current vaccine in preventing SIL/cancer and warts are highly needed in developing countries for appropriate policy decision making.

ART:

Anti-Retroviral Therapy

BMC:

Bugando Medical Centre

CD4:

Cluster of Differentiation type 4

CTC:

Care and Treatment Clinics

CUHAS:

Catholic University of Health and Allied Sciences

HAART:

Highly Active Ant-Retroviral Therapy

SIL:

Squamous Intraepithelial Lesions

HIV:

Human Immunodeficiency Virus

HPV:

Human Papilloma Virus

HR:

High Risk

LR:

Low Risk

RR:

Risk ratio

SIL:

Squamous Intraepithelial Lesion

SES:

Social Economic Status

WHO:

World Health Organization

The authors would like to acknowledge the support provided by CTC authority and Institute of clinical and molecular virology, Enlargen University, Schlossgarten for their technical support.

Funding

This study was supported by research grant from CUHAS to FM and Medical Mission Institute of Wuerzburg. Funders had no any role in this study.

Availability of data and materials

All data have been included in the manuscript.

Authors’ contributions

MMM, DM, CK and SEM participated in the design of the work. FM and DM participated in the collection of specimens and clinical data. FM KK, PR and AM performed laboratory analysis of the specimens. MMM, MM and SEM analyzed and interpreted the data. MMM and SEM wrote the first draft of the manuscript. MMM and SEM did critical revision of the manuscript. All authors read and approved the final version of the manuscript.

Authors’ information

None.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

The ethical clearance for the study (CREC/057/2014) was granted by the Joint catholic university of health and allied sciences/Bugando Medical Centre (CUHAS/BMC) research ethics and review committee (CREC). Importance of the study was explained to the participants and written informed consent was sought from each participant prior to enrolment.

Affiliations

1. Department of Obstetrics and Gynecology, Weill Bugando School of Medicine, P.O.Box 1464, Mwanza, Tanzania
   • Fridolin Mujuni
   •  & Dismas Matovelo
2. Department of Microbiology and Immunology, Weill Bugando School of Medicine, P.O.Box 1464, Mwanza, Tanzania
   • Mariam M. Mirambo
   •  & Stephen E. Mshana
3. Department of Pathology, Weill Bugando School of Medicine, P.O.Box 1464, Mwanza, Tanzania
   • Peter Rambau
4. Institute of Clinical and Molecular Virology, Enlargen University, Schlossgarten 4, 91054, Erlangen, Germany
   • Korn Klaus
5. Medical Mission Institute, Salvatorstrasse 7, 97067, Wuerzburg, Germany
   • Muller Andreas
   •  & Christa Kasang
6. Department of Microbiology and Immunology, Muhimbili University of Health and Allied Sciences, P.O. Box 65001, Dar es Salaam, Tanzania
   • Mtebe Majigo

Corresponding authors

Correspondence to Mariam M. Mirambo or Stephen E. Mshana.

Reprints and Permissions