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.

Table 1

Sequences of primers and probes for HPV detection and controls

Name

Sequence

Target

HPV GEN1 PCR

CP4

5´-ATG GTA CAR TGG GCA TWT GA-3´

HPV E1 (nt 1942-1961)

CP5

5´-GAG GYT GCA ACC AAA AMT GRC-3´

HPV E1 (nt 2400-2378)

HPV GEN2 PCR

PPF1

5´-AAC AAT GTG TAG ACA TTA TAA ACG AGC-3´

HPVE1 (nt 2082-2108)

PPR2

5´-ATT AAA CTC ATT CCA AAA TAT GA-3´

HPVE1 (nt 2336-2314)

Numbering is according to the sequence of HPV16W12E (GenbankIDNr. AF125673)

Table 2

Socio-demographic & Reproductive characteristics of the study population

Characteristics

Frequency/mean/median

Percentage (%)

aAge(years)

39.2 ± 9.1

bARV duration(years)

3(IQR:1-6)

Residence

Rural

158

62.0

Urban

97

38.0

SES

Low

117

45.9

High

138

54.1

H/STI

NO

213

83.5

YES

42

16.5

Parity

Nullporous

24

9.4

≥ 1 child

231

90.6

Level of education

Lower level

209

81.7

High level

46

18.3

CD4 cells/μl

> 200

134

52.6

200-100

99

38.8

< 100

22

8.6

Age at first intercourse

Below 18

181

71.0

Above 18

74

29.0

Contraceptive use

None

164

64.3

Hormonal

62

24.3

Condom

29

11.4

Marital status

Not married

144

56.5

Married

111

43.5

ARV use

None

24

9.4

≤ 6 months

23

9.0

> 6 months

208

81.6

amean, bmedian

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

Table 3

Genotypes distribution of HPV

HPV genotypes

Frequency

LR/HR

HPV 6

1

LR

HPV 16

14

HR

HPV 18

8

HR

HPV 31

9

HR

HPV 33

4

HR

HPV 34

1

HR

HPV 35

20

HR

HPV 40

2

LR

HPV 42

5

LR

HPV 45

10

HR

HPV 47

1

LR

HPV 51

3

HR

HPV 52

26

HR

HPV 53

5

HR

HPV 54

2

LR

HPV 56

5

LR

HPV 58

21

HR

HPV 59

3

HR

HPV 66

5

HR

HPV 68

2

HR

HPV 73

1

HR

HPV 74

1

LR

HPV 82

4

HR

HPV 83

3

LR

HPV 90

3

LR

HPV 103

1

Unclassified

Total

160

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.

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.

Acknowledgements

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.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.