BioMed Research International

BioMed Research International / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 2250690 |

James Kinoti Njue, Margaret Muturi, Lucy Kamau, Raphael Lwembe, "Human Papillomavirus Types Associated with Cervical Dysplasia among HIV- and Non-HIV-Infected Women Attending Reproductive Health Clinics in Eastern Kenya", BioMed Research International, vol. 2021, Article ID 2250690, 10 pages, 2021.

Human Papillomavirus Types Associated with Cervical Dysplasia among HIV- and Non-HIV-Infected Women Attending Reproductive Health Clinics in Eastern Kenya

Academic Editor: Jorge F. Quarleri
Received13 May 2021
Accepted03 Sep 2021
Published20 Sep 2021


Background. Human papillomavirus (HPV) causes over 99% of all cervical cancer globally. In 2019, it was responsible for 3286 deaths in Kenya. Data on the epidemiological distribution of HPV genotypes by cervical dysplasia and HIV-infected women which is important in designing prevention strategy monitoring treatment and management of cervical cancer is lacking in Eastern Kenya. Objective. To determine HPV genotype prevalence and their association with cervical dysplasia among HIV-infected (cases) and noninfected (control) women aged 18-48 years seeking reproductive healthcare. Methods. A cervical broom was softly rotated 360 degrees five times to exfoliate cells from the region of the transformation zone, squamocolumnar junction, and endocervical canal for HPV genotyping. Social-demographic and risk factors responsible for HPV acquisition were collected using a questionnaire. Laboratory outcome and questionnaire data statistical relationships were computed using Pearson chi-square test. Results. 317 women (cases: 161 (50.8%), control 156 (49.2%), mean age: 34.3,, range 18-46 years) were recruited from Embu (85/317 (26.8%)), Isiolo (64/317 (20.2%)), Kirinyaga (56/317 (17.7%)), Meru (81/317 (25.6%)), and Tharaka-Nithi (31/317 (9.8%)). The frequency HPV genotypes detected by cervical dysplasia were CIN1 (cases: HPV81 (12/317 (3.8%)), HPV11 (2/317 (0.6%)); control: HPV53 and 66 coinfection (1/317 (0.3%)), CIN2 (cases: HPV11, HPV16, HPV66 ((1/317 (0.3%) each), HPV81 (6/317 (1.9%)), and single case (1/317 (0.3%)) of HPV11 and 66, HPV81 and 44, HPV81 and 88, HPV9 and 53, and HPV16 and 58 coinfection; control: HPV81 (2/317 (0.6%)) and invasive cervical cancer (cases: HPV16 (1/317 (0.3%)) and HPV81 (3/317 (0.9%)); control: HPV16 and 66 (1/317 (0.3%))).Conclusions. There was a higher frequency of both high-risk and low-risk HPV genotypes associated with cervical dysplasia among HIV-infected than HIV-uninfected women seeking reproductive health care. This study provides epidemiological data on the existence of nonvaccine HPV types associated with cervical dysplasia in the region.

1. Introduction

Human papillomavirus (HPV) is primarily responsible for 99.7% of cervical cancers globally [1]. It can be sexually transmitted and causes cervical neoplastic changes leading to cervical cancer and is the second type of cancer in women aged 15-44 years in Kenya [1, 2]. It was responsible for 311,365 (8.2%) annual global mortality, 37,017 (16.9%) in East Africa, and 3250 (12.8%) in Kenya in 2019 [1, 2]. The cervical screening rate in Kenya is only 3.2% for women aged over 18 years [2].

Human papillomaviruses are grouped based on their oncogenicity as Group I: “Carcinogenic” and include HPV16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, and 59; Group 2A: “probably-carcinogenic” which includes HPV26, 53, 66, 67, 70, 73, 82, 30, 34, 69, 85 and 97; Group 2B: “possibly-carcinogenic” which includes HPV6 and 11; and Group 3: “unclassifiable to carcinogenicity” [2, 3].

Oncogenicity of HPV is further increased by immunocompromised status following human immunodeficiency virus (HIV) infection, long-term exposure to hormonal-contraceptives, and autoimmune drugs which favor tumorigenic effects of HPV-genome and weaken the immune system responsible of clearing cancer cells [3, 4]. Other risk factors include early sex debut, long-term inflammation caused by recurrent genital infections, hormonal changes during pregnancy, cervical trauma during labor, and failure to undergo HPV vaccination [35].

Literatures review wide spectra of HPV genotypes associated with cervical dysplasia among HIV-infected women [1, 57]. No such data is available from Eastern Kenya which has implications regarding cervical cancer burden and efficacy of HPV vaccination in the region [6]. Cross-protection between HPV types in the current vaccine and those that are predominant in the population has also not been established. This study therefore aimed at determining HPV types associated with cervical dysplasia among HIV-infected and noninfected women attending reproductive health clinics in the region.

2. Methods

2.1. Study Design and Participants

This cross-sectional study involved 317 women aged 18-46 years who attended reproductive health clinics in referral hospitals of Isiolo, Kirinyaga, Meru, Tharaka-Nithi, and Embu Counties in Eastern Kenya from January to December 2019. A sample size was calculated by using an ICO HPV infection rate (1) of 2.8% (total: 158 cases and 158 controls). Samples were distributed according to reproductive health clinic’s diagnosis 27 VIA/VILLI-positive cases starting July–December 2017 in the region. The stratified technique of sampling was used to obtain five strata based on county of residence then simple random sampling was then used to recruit the required number of a participant in each stratum.

2.2. Inclusion and Exclusion Criteria

Inclusion criteria are as follows: consenting women aged 18-46 years from Eastern Kenya.

Exclusion criteria are as follows: menstruating and pregnant women, mentally incompetent, and those with eroded cervix or history of ablative procedures or medical treatment for cervical disease in the last six months [6, 7].

2.3. Detection of Cervical Dysplasia

A two-step approach was applied where visual inspection with acetic acid solution (VIA) was performed followed by Lugol’s iodine solution (VILLI) which are based on the colors taken up by the cervical transformation zone [8, 9]. The outcome was reported as VIA/VILLI positive or negative.

2.4. Social-Demographic Data Collection

Social-demographic data was collected from participants with precaution after consenting using a questionnaire. The questionnaire focused on data on the residence, age, education level, parity, sexual orientation, and choice of family planning method.

2.5. HIV Determination

HIV serostatus was determined as per national algorithm. Alere Determine®HIV-1/2 test by Abbort Co. was used as baseline screening while the First Response® HIV 1-2-0 card test by Premier Medical Corporation was used as a confirmatory test in case of discrepancies. Uni-Gold™ Recombigen® HIV-1/2 by Trinity Biotech was used as a tie breaker test [10].

2.6. Collection and Storage of Cervical Exfoliated Cell Samples

External genitalia were examined while the participant lays in a lithotomic position. A speculum was rinsed in warm water, lubricated, and used to locate the cervical opening (os) under direct light. The mucus plug in os was removed and wiped to ensure sufficient cells were collected. A cervical broom (Dacron cervical broom; Digene Corporation, Silver Spring, Maryland STM™) was softly rotated 360 degrees five times to exfoliate cells from the region of the transformation zone, squamocolumnar junction, and endocervical canal [9]. Exfoliated cells were spread evenly and fixed immediately on a clean glass slide. The broom bristles were then dipped into aqueous minimum essential media (MEM) and the broom handle snapped so that it remained in the tightly closed vial and stored at 1-4°C [2, 5, 8].

2.7. Cytology

A standardized protocol for cytopathology request form [5, 6], Pap smears staining, and examination [6, 9] was followed for detection of nuclei and cytoplasm cytological changes following HPV infection. Cytopathologists supervised by a pathologist at Embu and Meru Hospitals were required to fill a pathology synoptic reporting form [6] using the Bethesda 2001 guidelines for reporting slides using a binocular microscope. Slides were later transferred to Kenya Medical Research Institute (KEMRI) for examination by pathologist to ensure the quality of cytology results. Pap smear results were classified as normal or abnormal (ASCUS, CIN1, CIN2, CIN3, or ICC). Cervical infections (candidiasis, cervicitis, trachomatis, and bacterial vaginitis) reported in cytology were exempted from the study [9, 11, 12].

2.8. HPV DNA Extraction and Amplification

All samples stored for HPV DNA detection were analyzed by the following procedure: (a)DNA extraction

A 96-well format HighPrep™ Viral-DNA/RNA, MagBio Genomics, Inc. USA/Canada Lysis kit was used. Master Mix was prepared by mixing 240 μL Lysis-buffer, 8 μLRNA carrier, and 280 μL isopropanol sample/volume. Samples stored at 1-4°C were thawed, vortexed for 5 minutes, and then centrifuged at 10000 r/min for 5 minutes to extract cytological material from the brush into MEM. Samples were subjected to the kit protocol to obtain DNA extracts by magnetic bead technique The eluate was stored at -20°C [13, 14]. (b)HPV DNA PCR

HPV detection was achieved by amplifying an L1 portion of the HPV genome that is relatively conserved through L1 consensus nested PCR in the ABI-thermocycler Model 9600; Applied Biosystems®. HPV consensus primary primers PGMY09 (GCACAGGGACATAACAATGG) and PGMY11 (CGTCCCAAAGGAAACTGATC) [15] that target the 450 bp region in the L1 ORF; the genome was included in the first PCR reaction. Additional primer sets targeting the same region of L1, MGP5+(ACGTTGGATGTTTGTTACTGTGGTGGATACTAC) and MGP6+(ACGTTGGATGGAAAAATAAACTGTAAATCATATTCCT) were used to produce shorter amplicon of ~160 bp [15]. 5 μM working stock of each primer was prepared by adding 50 μL biotinylated PGMY09 100 μM to 350 μL nanopure water and 50 μL PGMY11 100 μM primers to 750 μL nanopure water. They were later distributed each 5 μM working stock in 45–90 μL aliquots and stored at -20°C. In the primary PCR, 5 μL of the extract was amplified in a Master Mix containing PCR buffer (1x), 2.0 mM MgCl2, 100 μM dNTPs, 0.13 parts Taq polymerase enzyme and 500 nM MY09, and 500 nM MY11 (forward and reverse primer, respectively). First reaction: 4 minutes at 95°C (initial denaturation) then 30 cycles of 20 seconds at 95°C, 40 seconds at 56°C, and 2 minutes at 72°C, and final extension for 7 minutes at 72°C. In nested PCR, 5 μL of the first PCR product, 2 mM MgCl2, 500 nM GP5+, and 500 nM GP6+ (forward and reverse primer, respectively), 400 μM of dNTPs, and 0.13 units of Taq polymerase enzyme composed the Master Mix. Nested reaction: 4 minutes at 95°C (initial denaturation) then 30 cycles of 20 seconds at 95°C, 40 seconds at 60°C, 40 minutes at 72°C, and final extension for 7 minutes at 72°C. Positive control of CIN2+ and negative control of distilled water were incorporated in all primer cycles [15]. (c)Gel electrophoresis and UV visualization

Tris-Borate-EDTA 10x was prepared by dissolving 162 g Tris-base, 50 g boric acid, and 9.5 g EDTA in nanopure water to 1.0 liter volume (pH 8.8). 5 μL PCR-product in 4% agarose was used in gel electrophoresis where samples that showed a band of 160 bp were considered positive. The positive PCR product was purified using the QIAquick DNA purification kit™ (Qiagen, Germany) before sequencing [7, 13]. (d)HPV DNA sequencing

DNA sequencing was performed in ABI-thermocycler Model 9600 (Applied Biosystems) for 20 reaction cycles of 1 μL positive PCR-product, 1 μL of 5 μM GP6+ primer, 1 μL BigDye® Terminator, 3.5 μL buffer (5x), and 13.5 mL nanopure water according to the protocol supplied by the manufacturer. Centri-Sep column (Princeton Separations, Adelphia, NJ) was used for dye-terminator cleanup; then, the reaction mixture was sequenced in an automated ABI 3130 four-capillary genetic analyzer [15, 16]. (e)HPV genotyping and phylogenetic analysis

CHROMAS software Version 2.4.3 was used to edit the sequences obtained. HPV sequences were blasted in NCBI Sequences with unique divergences from the same HPV type underwent phylogenic analysis, and references were obtained from GenBank. The input file of representative sequences and their references was made and subjected for multiple alignments with CLUSTAL W in MEGA X software [15, 16]. The maximum likelihood method and the Tamura-Nei model were used to infer evolutionary history. Neighbor-Join and BioNJ algorithms were used to construct initial trees for the heuristic search of the matrix of pairwise distances by the maximum composite likelihood (MCL) method by selecting the topology with superior log likelihood value. Eighty-six nucleotide sequences were involved while codon positions were 1st+2nd+3rd [1618].

2.9. Ethics Approval

This study was approved by the KEMRI Scientific Ethical Review Unit (approval number: KEMRI/SERU/CVR/004/3342). Participants were consented orally where all steps and procedures for HIV and cervical exfoliated cell sample collection; analysis and collection of their results were explained. Laboratory and questionnaire data were confidentially stored by the principal investigator.

3. Results

3.1. Distribution of HPV Genotypes among HIV- and Non-HIV-Infected Women ()

Thirteen HPV types detected were low-risk HPV9, HPV11, HPV81, HPV66, HPV87, and HPV88 while high-risk types were HPV16, HPV53, HPV61, HPV45, HPV52, and HPV58. There was significant association between HPV types detected and HIV serostatus (, ) (Figure 1).

3.2. Distribution of HPV Genotypes with Cytology (Levels of Cervical Dysplasia) among HIV- and Non-HIV-Infected Women

A total of 62 (19.2%) HIV-infected women were infected with single HPV type compared to (93.2%) HIV-uninfected women while 73 (23.03%) were infected by multiple HPV types compared with 13 (4.1%) HIV-uninfected women (, ) (Table 1).

Infection type (%)HIV negativeHIV positive value

Single HPV type infection
 Low-risk type67 (21.1)9 (2.8): HPV817 (2.2): HPV11<0.001
51 (16.1): HPV81
 High-risk type4 (1.3)3 (0.1): HPV16
1 (0.3): HPV66
Total (single HPV types)71 (22.4)9 (3.2)62 (19.2)
Multiple HPV types of infection
 Low-risk type2 (0.6)1 (0.3): HPV81, 44
1 (0.3): HPV81, 88
 High-risk type9 (2.8)3 (0.9): HPV16, 664 (1.3): HPV16, 66
1 (0.3): HPV53, 66
1 (0.3): HPV16, 58
 Low- and high-risk type4 (1.2)1 (0.3): HPV11, 16, 53, 61, 811 (0.3): HPV9, 53
1 (0.3): HPV11, 66
1 (0.3): HPV11, 45, 52, 87
Total (multiple HPV types)15 (4.7)4 (1.2)11 (3.5)
Total ()86 (27.12)13 (4.1)73 (23.03)

: the probability at the 0.001 level.
3.3. Social Demographic and HPV-Associated Risk Factors among HIV- and Non-HIV-Infected Women

A total of 317 (mean age: 34.3years, , range 18-46) women were recruited. Overall HPV prevalence was 27.1% ((86/317); cases: 23.2% (73/317) and control: 4.1% (13/317)). There was significant association between HPV infection and age, parity, education level, and reported number of sex partners (Table 2).

Category (%)TotalHIV negative ( (%))HIV positive ( (%)]
HPVHPV statusTotal valueHPV statusTotal value

Embu85 (26.8)27 (8.5)36 (11.4)5 (1.6)41 (13.0)0.85322 (6.9)22 (6.9)44 (13.8)0.383
Isiolo64 (20.2)18 (5.7)35 (11.0)3 (0.9)38 (11.9)11 (3.5)15 (4.7)26 (8.2)
Kirinyaga56 (17.7)12 (3.8)22 (6.9)1 (0.3)23 (7.2)22 (6.9)11 (3.5)33 (10.4)
Meru81 (25.6)20 (6.3)37 (11.7)3 (0.9)40 (12.6)24 (7.6)17 (5.4)41 (13.0)
T.Nithi31 (9.8)9 (2.8)13 (4.1)1 (0.3)14 (4.4)9 (2.8)5)17 (5.3)
Age (mean: 34.3, range 18-46)
<35162 (51.2)60 (18.9)59 (18.6)10 (3.2)69 (21.8)0.01443 (13.6)50 (15.8)93 (29.4)0.009
>35155 (48.8)26 (8.2)84 (26.5)3 (0.9)87 (27.4)45 (14.2)23 (7.3)68 (21.5)
Christian255 (80.4)70 (22.1)108 (34.1)10 (3.1)118 (37)0.60777 (24.3)60 (18.9)137 (43.2)0.237
Muslim62 (19.6)16 (5.0)35 (11.1)3 (0.9)38 (12.0)11 (3.5)13 (4.1)24 (7.6)
Education level
Primary96 (30.3)20 (6.3)40 (12.6)3 (0.9)43 (13.5)0.11536 (11.4)17 (5.4)53 (16.8)0.004
Secondary135 (42.6)39 (12.3)66 (20.8)4 (1.3)70 (22.1)30 (9.5)35 (11.0)65 (19.5)
College67 (21.1)21 (6.6)33 (10.4)4 (1.3)37 (11.7)13 (4.1)17 (5.4)30 (9.5)
University19 (6.0)6 (1.9)4 (1.3)2 (0.6)6 (1.9)9 (2.8)4 (1.3)13 (4.2)
Contraceptive use
Other223 (70.3)66 (20.8)95 (30.0)8 (2.5)103 (32.5)0.46962 (19.6)58 (18.3)120 (37.9)0.131
Hormonal94 (29.7)20 (6.3)48 (15.1)5 (1.6)53 (16.7)26 (8.2)15 (4.7)41 (12.7)
>368 (21.5)16 (5.0)30 (9.5)8 (2.5)38 (12.0)0.00322 (6.9)39 (12.3)41 (39.2)0.001
<3249 (78.5)70 (22.1)113 (35.6)5 (1.6)118 (37.2)66 (20.8)34 (10.7)100 (21.5)
Marital status
Married226 (71.3)55 (17.4)109 (34.4)8 (2.5)117 (36.9)0.24962 (19.6)47 (14.8)109 (34.4)0.299
Separated32 (10.1)8 (2.5)13 (4.1)1 (0.3)14 (4.4)11 (3.5)7 (2.2)18 (5.7)
Single41 (12.9)12 (3.8)17 (5.4)2 (0.6)19 (6.0)12 (3.8)10 (3.2)22 (7.0)
Divorced6 (1.9)3 (0.9)2 (0.6)1 (0.3)3 (0.9)1 (0.3)2 (0.6)3 (0.9)
Widowed12 (3.8)8 (2.5)2 (0.6)1 (0.3)3 (0.9)2 (0.6)7 (2.2)9 (2.8)
Sex partners
1186 (58.7)43 (13.6)91 (28.7)7 (2.2)98 (30.9)0.33952 (16.4)36 (11.4)88 (27.8)0.014
>1131 (41.3)43 (13.6)52 (16.4)6 (1.9)58 (18.3)36 (11.4)37 (11.7)73 (23.1)
Total317 (100.0)86 (27.1)143 (45.1)13 (4.1)156 (49.2)88 (27.8)73 (23.2)161 (50.8)0.001

N: negative; P: positive; : the probability at the 0.001 level; : the probability at the 0.005 level.
3.4. Distribution of HPV Genotypes among HIV- and Non-HIV-Infected Women

HPV type distribution among HIV-infected women was significantly associated by residence (Embu (), Isiolo (), Meru (), and Tharaka-Nithi ()), age (≤30 years (); <30 years ()), parity (≤35 (); >35 ()), family planning method (hormonal (, ); nonhormonal contraceptives () use), and number of sex partners (1 (; >1 ()) (Table 3).

CategoryHIV testHPV infectionotal value
High-risk HPV typesLow-risk HPV types

 EmbuN1 (0.9)1 (0.9)4 (3.8)6 (5.7)0.042
P3 (2.8)1 (0.9)2 (1.9)1 (0.9)2 (1.9)1 (0.9)4 (3.8)1 (0.9)13 (12.8)1 (0.9)29 (27.6)
 IsioloN2 (1.9)3 (2.8)1 (0.9)6 (5.7)0.001
P1 (0.9)3 (2.8)3 (2.8)1 (0.9)11 (10.5)1 (0.9)20 (19.0)
 KirinyagaN1 (0.9)1 (0.9)0.076
P2 (1.9)8 (7.6)10 (9.5)
 MeruN1 (0.9)1 (0.9)1 (0.9)1 (0.9)3 (2.8)7 (6.6)0.001
P2 (1.9)15 (14.3)17 (16.2)
 T.NithiN1 (0.9)1 (0.9)0.048
P2 (1.9)6 (5.7)8 (7.6)
 <35N3 (2.8)1 (0.9)2 (1.9)1 (0.9)5 (4.8)12 (11.4)0.001
P3 (2.8)1 (0.9)2 (1.9)5 (4.8)1 (0.9)8 (7.6)1 (0.9)39 (37.1)1 (0.9)61 (58.1)
 >35N1 (0.9)1 (0.9)1 (0.9)4 (3.8)7 (6.6)0.017
P5 (4.8)2 (1.9)1 (0.9)1 (0.9)1 (0.9)14 (13.3)24 (22.8)
Family planning
 hormonalN1 (0.9)1 (0.9)4 (3.8)6 (5.7)0.047
P1 (0.9)2 (1.9)3 (2.8)1 (0.9)2 (1.9)10 (9.5)19 (18.1)
 OtherN3 (2.8)1 (0.9)2 (1.9)1 (0.9)1 (0.9)5 (4.8)13 (12.8)0.001
P7 (6.6)1 (0.9)1 (0.9)4 (3.8)7 (6.6)1 (0.9)1 (0.9)43 (40.9)1 (0.9)1 (0.9)67 (63.8)
 <3N1 (0.9)1 (0.9)2 (1.9)4 (3.8)0.001
P2 (1.9)1 (0.9)3 (2.8)1 (0.9)8 (7.6)15 (14.3)
 >3N3 (2.8)1 (0.9)2 (1.9)1 (0.9)7 (6.6)14 (13.3)
P6 (5.7)1 (0.9)2 (1.9)1 (0.9)6 (5.7)1 (0.9)6 (5.7)1 (0.9)1 (0.9)45 (42.8)1 (0.9)1 (0.9)72 (68.6)
Number of sex partners
 OneN3 (2.8)3 (2.8)4 (3.8)10 (9.5)0.001
P3 (2.8)1 (0.9)1 (0.9)1 (0.9)5 (4.8)1 (0.9)28 (26.8)1 (0.9)1 (0.9)42 (40.0)
 >oneN1 (0.9)1 (0.9)1 (0.9)5 (4.8)8 (7.6)0.001
P5 (4.8)1 (0.9)6 (5.7)1 (0.9)4 (3.8)1 (0.9)25 (23.8)43 (40.9)
TotalN4 (3.8)3 (2.8)1 (0.9)1 (0.9)9 (8.6)18 (17.1)0.001
P8 (7.6)1 (0.9)1 (0.9)1 (0.9)7 (6.6)1 (0.9)9 (8.6)1 (0.9)1 (0.9)53 (50.5)1 (0.9)1 (0.9)85 (80.9)

N: negative; P: positive; : the probability at the 0.001 level; : the probability at the 0.005 level; T. Nithi: Tharaka-Nithi County.
3.5. Association of Cervical Cytology with Other Clinical Reproductive Health Ailments

A total of 96 (30.3%) HIV-positive women had normal cytology as compared to 143 (45.1%) HIV-negative, whereas 65 (20.5%) HIV-positive women had abnormal cytology results compared with 13 (4.1%) HIV-negative () (Figure 2).

3.6. Association of HPV Genotypes with Cervical Dysplasia among HIV-Infected Women and Noninfected Women

Frequency of common HPV genotypes detected by cervical dysplasia among HIV-infected women were as follows: CIN1 (17 (5.4%)): HPV81 (12 (3.8%)), HPV11 (2 (0.6%)), CIN2 (16 (5.0%)): single case (0.3%) of HPV11, HPV16, and HPV66 and HPV81 (6 (1.9%)), CIN3 (6 (1.8%)): HPV11 (1 (0.3%)), HPV81 (4 (1.2%)) and invasive cancer (5 (1.5%)): HPV16 (1 (0.3%)) and HPV81 (3 (0.9%)), () (Table 4).

HIVHPV type of infectionNormal cytology ( (%))Abnormal cytology ( (%))Total value

NegativeSingle type
HPV813 (0.9)1 (0.3)2 (0.6)2 (0.6)1 (0.3)9 (2.7)0.001
Multiple types
HPV16, 662 (0.6)1 (0.3)3 (0.9)
HPV11, 16, 53, 81, 611 (0.3)1 (0.3)
HPV positive6 (1.9)1 (0.3)2 (0.6)2 (0.6)1 (0.3)1 (0.3)13 (4.1)
HPV negative137 (43.2)6 (1.8)143 (45.1)
Total (HPV positive and negative)143 (45.1)7 (2.1)2 (0.6)2 (0.6)1 (0.3)1 (0.3)156 (49.2)

PositiveSingle type
HPV112 (0.6)1 (0.3)2 (0.6)1 (0.3)1 (0.3)7 (2.2)0.001
HPV161 (0.3)1 (0.3)1 (0.3)3 (0.9)
HPV661 (0.3)1 (0.3)
HPV8121 (6.6)5 (1.6)12 (3.8)6 (1.9)4 (1.2)3 (0.9)51 (16.1)
Multiple types
HPV11, 661 (0.3)1 (0.3)
HPV81, 441 (0.3)1 (0.3)
HPV81, 881 (0.3)1 (0.3)
HPV9, 531 (0.3)1 (0.3)
HPV16, 581 (0.3)1 (0.3)
HPV16, 663 (0.9)1 (0.3)4 (1.2)
HPV66, 531 (0.3)1 (0.3)
HPV11, 45, 52, 871 (0.3)1 (0.3)
HPV positive28 (8.8)7 (2.1)15 (4.7)14 (4.4)5 (1.6)4 (1.2)73
HPV negative68 (21.4)20 (6.3)88 (27.8)
Total (HIV positive and negative)96 (30.3)27 (8.5)15 (4.7)14 (4.4)5 (1.6)4 (1.2)161 (50.8)
Total (HIV positive and negative)239 (75.4)34 (10.7)17 (5.4)16 (5.0)6 (1.8)5 (1.5)317 (100.0)

ASCUS: atypical cells of unknown significance; CIN: cervical intraepithelial neoplasia; ICC: invasive cervical cancer; : the probability at the 0.001 level.
3.7. Distribution of HPV Genotypes among HIV- and Non-HIV-Infected Women with Cervical Dysplasia

Phylogenetic tree of HPV samples marked in red aligned against the representation of the different HPV genotypes and subtypes distributed in different regions worldwide, constructed using MEGA6 neighbor-joining using 1000 bootstrap reference method. Most HPV81 clustered with those cases was detected in Bangkok Morocco and Thailand, HPV66 clustered with those from Kenya and Iran, while HPV66 clustered with cases reported in Tunisia, Morocco, Iran, and India (Figure 3).

4. Discussion

This study established an overall HPV infection rate of 27.12% (cases: 73 (23.03%); control: 13 (4.1%)). It agrees with the overall age-specific HPV infection rate of 27% in neighbouring Nairobi region [1]. HIV infection was significantly associated with HPV infection. A high HPV prevalence was reported in Embu and Meru Counties where HIV was most prevalent. There is a need to establish HPV prevalence in other regions of the country with high HIV burden [19] that could be harbouring more HPV infection rate, hence raising overall national HPV prevalence (40%) [1]. This indicator should be prioritised by public health interventions to reduce cervical cancer morbidity and mortality.

Women aged below 35 years of age had a high rate of mixed HPV genotypes and a significant association between HIV infection and abnormal cytology outcome which agrees with published observations [5]. A possible explanation is that HIV infection may have facilitated cervical dysplasia in young women, and an inverse relationship of high-risk HPV prevalence and age has been described [20].

A high infection rate was also noted among Christians who comprised the majority of participants than Muslims. This study concurs with other studies that have shown high HPV infection rate among secondary school educated women. This is attributed to a higher risky sexual behaviour and low knowledge on HPV and associated risk factors as reported among young women in South Africa [20, 21].

HPV infection rate was high among HIV-infected women with high parity, married, and those using hormonal contraceptive than other categories. This could be associated with high persistence and subsequent low clearance of HPV upon infection among these groups as reported by other studies [5, 22].

There was a significant association between having multiple sex partners and HPV infection. Most of the participants in this category were HPV infected. This could be because having multiple sex partners increased their chance of contracting HPV and HIV. The dominance of HPV infection in normal cytology was detailed high (28.8%) as published [1]. Younger women showed a higher infection rate in normal cytology than their older counterparts. This is because they are more sexually active and highly likely to encounter wider spectra of HPV types [7, 15]. However, they are reported to have a high HPV clearance rate than their older counterparts [23] hence less likely to develop cervical dysplasia.

The foremost common HPV detected in the region was lrHPV81—Alphapapillomavirus 3 (α3) HPV species by both HIV serostatus and cytology outcomes. The predominance of HPV81 (CP8304) has moreover been detailed in women with anomalous cytology in Kenya [6, 7] and Qatar [18]. It is associated with precancerous and cancerous lesions and mostly detected among immunocompromised women [7, 24, 25] as seen in this study.

HPV16 (α9) and HPV66 (α6) were most common high-risk types associated with invasive cancer by HIV serostatus. Their predominance in abnormal cytology among HIV-infected women has been established [7, 15, 24, 25]. HPV66 is classified as oncogenic because of its close phylogenetic relationship with HPV56 (α9). This study and others therefore recommend routine testing for HPV16, HPV52, and HPV66 [20] since these members of (α7) and (α9) dominate neoplastic tissues. Though they were also detected in normal cytology, a study has demonstrated the genetic variability and frequent mutating L1, E6, and E7 genes mostly HPV52 and HPV53 (α7) resulting into intratypic variants with increased oncogenic potential [21].

Single HPV type infection in CIN1+ showed diversity compared with multiple HPV type’s infection in CIN1+ by HIV infection. Diverse infection was also seen among HIV-infected women with multiple HPV infections. Other common lrHPV types detected were HPV11, HPV44 (α10), and HPV88 which do not feature in many studies as potential oncogenic types Here, their association with CIN+ could have been increased by HIV infection as reported in other studies [5, 15, 21].

There were a significantly low number of participants infected with HPV types of the same alpha group in this study. There were also a reduced number of hrHPV52 and 58 and absence of most of other hrHPV types, notably HPV18. This could be attributed to antibodies’ cross-reaction which is achieved within species unlike across species. Antibodies against HPV16 cross-react and protect against HPV52 and HPV58 of the same α9 species, HPV18 cross-react with HPV66, and HPV53 (α7) species and HPV81 will cross-react with HPV61 [26].

This study and literature data in Kenya establish the predominance of mixed HPV types other than those included in the current vaccine [15]. This puts doubt on the effectiveness of the current vaccine in protection and/or cross-protection against hrHPV among women with compromised immunity.

The phylogenetic trees illustrated that single or multiple HPV types infected each participant as shown by 15 participants. Another reason is that HPV genes are replicated by the host replication machinery suggesting that a very low human autosomal–like mutation rate would be operating as established in another study [21].

5. Limitation

The study group is a minimal presentation, and data presented cannot be generalized as exact outcome if all women of reproductive age in Eastern Kenya were sampled. Responses to the questionnaire could have been altered by missing key words like ‘cervical cancer’ in local direct during translation but caution was taken to give accurate translation and meaning.

6. Conclusions

There was a higher frequency of both high-risk and low-risk HPV genotypes associated with cervical dysplasia among HIV-infected than HIV-uninfected women seeking reproductive health care. This study provides epidemiological data on the existence of nonvaccine HPV types associated with cervical dysplasia in Eastern Kenya. Genotypic HPV diversity established has potential impact on new vaccines design and need for increased cervical screening among HIV-infected women.

Data Availability

The datasets are available from the corresponding author on reasonable request.


The funding institution had no role whatsoever in designing the study, sample and data analysis, or writing of the manuscript.

Conflicts of Interest

The authors declare that they have no conflicting interests.

Authors’ Contributions

NJK, MM, LK, and RL designed the study. NJK conducted the survey and laboratory analysis, interpreted the data, and wrote the main manuscript text. All authors reviewed the final manuscript.


Special thanks was given to Zakayo Maingi, Caroline Baliach, and Alex Maiyo for their assistance in data analysis. Special thanks is also given to the National Research Fund for funding this project and Kenyatta University for accepting this project as an academic proposal. We thank all staff of all County Referral Hospitals in Eastern Kenya for enabling us to recruit their clients. Special thanks is given to KEMRI Laboratories for the support in sample analysis. This study was supported by the Kenya National Council for Science Technology and Innovation.


  1. ICO HPV Information Centre, Human Papillomavirus and Related Diseases Report-World, HPV Information Centre, 2019.
  2. J. Ferlay, M. Ervik, F. Lam et al., Global Cancer Observatory: Cancer Today, International Agency for Research on Cancer, Lyon, France, 2018,
  3. D. Stelzle, L. F. Tanaka, K. K. Lee et al., “Estimates of the global burden of cervical cancer associated with HIV,” The Lancet Global Health, vol. 9, no. 2, pp. e161–e169, 2021. View at: Publisher Site | Google Scholar
  4. J. Lei, “HPV vaccination and the risk of invasive cervical cancer,” The New England Journal of Medicine, vol. 383, pp. 1340–1348, 2020. View at: Publisher Site | Google Scholar
  5. K. van Doorslaer, Z. Chen, H. U. Bernard et al., “ICTV virus taxonomy profile: Papillomaviridae,” Journal of General Virology, vol. 99, no. 8, pp. 989-990, 2018. View at: Publisher Site | Google Scholar
  6. H. Mabeya, S. Menon, S. Weyers et al., “Uptake of three doses of HPV vaccine by primary school girls in Eldoret, Kenya; a prospective cohort study in a malaria-endemic setting,” BMC Cancer., vol. 18, no. 1, 2018. View at: Publisher Site | Google Scholar
  7. A. Omire, “Cervical Dysplasia, Infection, and Phylogeny of Human Papillomavirus in HIV- Infected and HIV-Uninfected Women at a Reproductive Health Clinic in Nairobi, Kenya,” BioMed Research International, vol. 2020, 10 pages, 2020. View at: Publisher Site | Google Scholar
  8. E. Krug and C. Varghese, Guide for Establishing a Pathology Laboratory in the Context Of Cancer Control, World Health Organization, Geneva, Switzerland, 2019.
  9. R. Jug and S. M. Bean, “Bethesda system,” 2021, View at: Google Scholar
  10. National AIDS and STI Control Program, Guidleines for HIV Testing Services in Kenya, NASCOP, Nairobi, Kenya, 2020.
  11. Z. Pan, Y. Song, X. Zhe et al., “Screening for HPV infection in exfoliated cervical cells of women from different ethnic groups in Yili, Xinjiang, China,” Scientific Reports, vol. 9, no. 1, p. 3468, 2019. View at: Publisher Site | Google Scholar
  12. X. Huang, C. Li, F. Li, J. Zhao, X. Wan, and K. Wang, “Cervicovaginal microbiota composition correlates with the acquisition of high-risk human papillomavirus types,” International Journal of Cancer, vol. 143, no. 3, pp. 621–634, 2018. View at: Publisher Site | Google Scholar
  13. T. Shibata and M. Nakagawa, “Evaluation of DNA extraction protocols from liquid-based cytology specimens for studying cervical microbiota,” bioRxiv, vol. 10, no. 1101, article 921619, 2020. View at: Publisher Site | Google Scholar
  14. H. Onywera, A. L. Williamson, Z. Z. A. Mbulawa, D. Coetzee, and T. L. Meiring, “The cervical microbiota in reproductive-age South African women with and without human papillomavirus infection,” Papillomavirus Research, vol. 7, pp. 154–163, 2019. View at: Publisher Site | Google Scholar
  15. W. S. Chan, T. L. Chan, C. H. Au et al., “An economical nanopore sequencing assay for human papillomavirus (HPV) genotyping,” Diagnostic Pathology, vol. 15, no. 1, p. 45, 2020. View at: Publisher Site | Google Scholar
  16. S. Menon, A. Wusiman, M. C. Boily et al., “Epidemiology of HPV genotypes among HIV positive women in Kenya: a systematic review and meta-analysis,” PLoS One, vol. 11, no. 10, article e0163965, 2016. View at: Publisher Site | Google Scholar
  17. D. Hawkes, M. H. T. Keung, Y. Huang et al., “Self-collection for cervical screening programs: from research to reality,” Cancers, vol. 12, no. 4, p. 1053, 2020. View at: Publisher Site | Google Scholar
  18. A. A. Elmi, D. Bansal, A. Acharya et al., “Human papillomavirus (HPV) infection: molecular epidemiology, genotyping, Seroprevalence and associated risk factors among Arab women in Qatar,” PLoS One, vol. 12, no. 1, article e0169197, 2017. View at: Publisher Site | Google Scholar
  19. M. Tsehaiu, UNAIDS, Kenya HIV and AIDS Estimates, 2020.
  20. A. Molina-Pineda, M. G. López-Cardona, L. P. Limón-Toledo et al., “High frequency of HPV genotypes 59, 66, 52, 51, 39 and 56 in women from Western Mexico,” BMC Infectious Diseases, vol. 20, no. 1, p. 889, 2020. View at: Publisher Site | Google Scholar
  21. J. Zhang, S. Zhang, M. Wang et al., “Genetic variability in E5, E6, E7 and L1 genes of human papillomavirus type 31,” Molecular Medicine Reports, vol. 17, pp. 5498–5507, 2018. View at: Publisher Site | Google Scholar
  22. M. Akaaboune, B. Kenfack, M. Viviano et al., “Clearance and persistence of the human papillomavirus infection among Cameroonian women,” Womens Health, vol. 14, 2018. View at: Publisher Site | Google Scholar
  23. M. Li, T. Liu, G. Luo et al., “Incidence, persistence and clearance of cervical human papillomavirus among women in Guangdong, China 2007-2018: A retrospective cohort study,” Journal of Infection and Public Health, vol. 14, no. 1, pp. 42–49, 2021. View at: Publisher Site | Google Scholar
  24. G. Kovacevic, V. Milosevic, N. Nikolic et al., “The prevalence of 30 HPV genotypes detected by EUROArray HPV in cervical samples among unvaccinated women from Vojvodina province, Serbia,” PLoS One, vol. 16, no. 4, article e0249134, 2021. View at: Publisher Site | Google Scholar
  25. Z. Z. A. Mbulawa, N. I. Somdyala, S. A. Mabunda, and A.-L. Williamson, “High human papillomavirus prevalence among females attending high school in the Eastern Cape Province of South Africa,” PLoS One, vol. 16, no. 6, article e0253074, 2021. View at: Publisher Site | Google Scholar
  26. H. Kann, M. Lehtinen, T. Eriksson, H.-M. Surcel, J. Dillner, and H. Faust, “Sustained cross-reactive antibody responses after human papillomavirus vaccinations: up to 12 years follow-up in the Finnish Maternity Cohort,” The Journal of Infectious Diseases, vol. 223, no. 11, pp. 1992–2000, 2021. View at: Publisher Site | Google Scholar

Copyright © 2021 James Kinoti Njue et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.