Human papillomavirus prevalence and genotype distribution among young women and men in Maputo city, Mozambique.

OBJECTIVES: Human papillomavirus (HPV) is a well-known cause of cervical cancer, the second most frequent cancer in female African populations. This study aimed at determining the prevalence of HPV infections and the genotype distribution in young adults aged 18-24, in Maputo city, Mozambique, and to assess the suitability of commercially available HPV vaccines.
METHODS: This cross-sectional study was conducted between 2009 and 2011 at a youth clinic in Maputo Central Hospital. Cervical and urethral samples were obtained from 236 women and 176 men, respectively. Demographic and behavioural data were collected using structured questionnaires. HPV genotyping was performed for 35 different high, probably or possibly high-risk and low-risk HPV types using the CLART Human Papillomavirus 2.
RESULTS: HPV prevalence was 168/412 (40.8%; 95% CI 36.0 to 45.5) and was significantly higher in women than in men (63.6%vs10.2%). HPV52 was the most frequent type found in women, followed by HPV35, -16,-53, -58,-6 and -51. In men, HPV51 ranked the highest, followed by HPV6, -11,-52, -59 and -70. HIV infection and sexual debut before 18 years of age were associated with multiple HPV infections (OR 3.03; 95% CI 1.49 to 6.25 and OR 6.03; 95% CI 1.73 to 21.02, respectively). Women had a significantly higher HPV infection prevalence than men (p<0.001). The 9-valent HPV vaccine would cover 36.8% of the high-risk genotypes circulating in women in this study, compared with 26.3% and 15.8% coverage by the bivalent and quadrivalent vaccines, respectively.
CONCLUSION: This study confirmed the high burden of HPV infections in young women in Maputo city, Mozambique. The HPV prevalence was associated with high-risk sexual behaviour. Sex education and sexually transmitted infection prevention interventions should be intensified in Mozambique. Only a proportion of the high-risk HPV genotypes (37%) were covered by currently available vaccines. © Article author(s) (or their employer(s) unless otherwise stated in the text of the article) 2017. All rights reserved. No commercial use is permitted unless otherwise expressly granted.

This is the first study to describe the human papillomavirus (HPV) prevalence and genotype distribution in young adults (aged 18–24) in Mozambique.

This study provides insights on the HPV types circulating in young populations in Maputo city, which is an important public health information particularly for discussions on HPV vaccine introduction in the country.

The absence of an older cohort (aged ≥25 years) for comparison of circulating HPV genotypes was a limitation of this study.

The lower HPV prevalence in the male population may have been related to the site of sample collection, therefore additional studies in male population in Mozambique, using different anatomic collection sites should be considered.

There is a high burden of human papillomavirus (HPV) infection in young women in Maputo city, Mozambique.

High-risk sexual behaviour was demonstrated to be associated with high prevalence of HPV infection.

There is a need to intensify sex education and prevention of sexually transmitted infections in young adults.

Currently available HPV vaccines would only cover a proportion (37%) of the prevalent high-risk genotypes.

Introduction

Globally, cervical cancer (CC) is the fourth most common cancer in women and the second most common cancer in female African populations.1 More than 75% of cases occur in less developed areas. East Africa is the most affected region with an age-standardised incidence rate of 42.7 per 100 000 women. Mozambique has the second highest incidence of CC after Malawi.2

Human papillomavirus (HPV) infection has been implicated in virtually all CC cases. To date, 201 different HPV types have been identified.3 More than 40 of these types can infect the anogenital tract and have been classified as low-risk (lrHPV6, −11, –40, −42, –43, −44, –54, −61, –62, −71, –72, −81, –83, −84 and −89), probably or possibly high-risk (phr-HPV26, −30, –34, −53,– 66, −67, –68, −69, –70, −73, –82, −85 and −97) and high-risk HPV types (hr-HPV16, −18, –31, −33, –35, −39, –45, −51, –52, −56, –58 and −59) depending on their ability to lead to malignant progression.4–6

Worldwide, HPV16 has been responsible for more than half (54.4%) of the CC cases followed by HPV18 (16.5%), HPV58 (5.1%), HPV33 (4.7%), HPV45 (4.4%), HPV31 (3.6%), HVP52 (3.4%) and HPV35 (1.9%).7 However, only a small proportion of infected women develop the disease. Persistent infection with hr-HPV types and integration of HPV–DNA in the host cells are usually required for malignant transformation of the cervix. The time between infection and initial precancer lesions is around 7–10 years. Cervical cancer is more frequently diagnosed in middle-aged women (40–64 years), thus suggesting that infections are taking place early in life. Globally, around 20% of all cervical cancers are being diagnosed between the ages of 15–39 years.8 This suggests that a proportion of infections may be occurring even earlier.

Previous reports have shown a high prevalence of HPV infections in younger populations9 10 and reports in southern Africa demonstrated a high prevalence of HPV infections among women of 30 years or younger.11 12 Limited information is available about the HPV epidemiology in Mozambique particularly in men and young populations in urban settings. Castellsagué et al have shown that the HPV prevalence is high (75%) in women aged 14–61 years in rural southern Mozambique.13 Cervical cancer is the most frequent cancer in Mozambican women.14 Others have described the contribution of HPV infection in the development of cervical cancers in Mozambique.15–18

By the end of 2016, HPV vaccination had only been introduced in 65 countries worldwide, seven being in the African continent. In Mozambique, HPV vaccine has yet not been introduced in the national expanded programme on immunisation. Currently, 23 African countries, including Mozambique, are conducting or have recently finalised their HPV vaccine pilot projects.1

The following three HPV vaccines are available: Gardasil (which confers protection against HPV6, −11, –16 and −18 and cross-protection against HPV31); Cervarix (which prevents against HPV16 and −18 and provides high and moderate cross-protection against HPV31, −45 and −33, respectively); and Gardasil 9 (which protects against HPV6, −11, –16, −18, –31, −33, –45, −52 and −58). Gardasil 9 is a second-generation 9-valent vaccine that can theoretically prevent 87% of all hr-HPV infections in Africa.19 Nonetheless, the HPV distribution varies between African regions, and the identification of country-specific profiles and their adequacy to the currently available vaccines remains of great importance when deciding on vaccine introduction.

In the present study, we aimed to (1) determine the prevalence and distribution of HPV infections in sexually active young adults and (2) determine the suitability of the current HPV vaccines in the context of the Mozambican epidemic.

Subjects and methods

This report adheres to the Strengthening The Reporting of OBservational Studies in Epidemiology (STROBE) guidelines for the reporting of observational studies.

Study design and population

This was a cross-sectional ancillary study of an HIV incidence study conducted at a youth clinic in Maputo city, Mozambique, from August 2009 to October 2011.20 In the parent study, briefly, 1380 males and females aged 18–24 years were screened for HIV, syphilis and hepatitis B virus. HIV-seronegative subjects were enrolled and followed for 1 year with quarterly visits to the clinic for assessment of HIV serostatus. Baseline demographic and behavioural data were collected using structured questionnaires. In the present study, samples from male and female subjects were collected at one time point, at screening or at one of the follow-up visits. We aimed at sampling 500 individuals with equal gender distribution (250 females and 250 males) based on the possibility of recruiting male subjects at the youth clinic, which mainly caters to women. Cervical samples were collected in female participants via speculoscopy, using a Rovers Viba-Brush (Rovers Medical Devices BV, Oss, the Netherlands). Urethral samples were collected from male participants by gently inserting a cotton swab approximately 2–4 cm in the urethral meatus and rotating it in one direction. Both the brushes and swabs were immersed in 5 mL of SurePath cell-preservation solution (TriPath Imaging, Burlington, North Carolina, USA), transported to the study laboratory within the same day of collection and stored at +4°C to +8°C for 3 months and then at −80°C.

Laboratory testing

DNA was extracted from an initial sample volume of 1.5 mL using a QIAamp DNA Mini Kit (Qiagen, GmbH, Hilden, Germany) according to the manufacturer’s instructions and eluted to a final volume of 200 µl.

HPV genotyping was performed using the CLART Human Papillomavirus 2 (Genomica, Madrid, Spain), a low-density microarray platform based on PCR amplification of a 450 bp fragment within an HPV L1 highly conserved region from 35 different HPV types (HPV6, −11, –16, −18, –26, −31, –33, −35, –39, −40, –42, −43, –44, −45, –51, −52, –53, −54, –56, −58, –59, −61, –62, −66, –68, −70, –71, −72, –73, −81, –82, −83, –84, −85, and −89) and a human gene control (CTFR). Individual genotyping results were analysed on the Clinical Array Reader (Genomica, Madrid, Spain). Adequacy of samples was assessed by amplification of the CTFR. Samples with undetectable DNA were rerun and the second result was considered final. The CLART Human Papillomavirus 2 assay complies with the European Union (EU) safety, health or environmental requirements, with the EU legislation and with the European In-Vitro Diagnostic Devices Directive. This assay has been shown to have a similar performance as that of other well-established HPV screening assays, such as the Food and Drug Administration-approved Hybrid Capture 2 test.21 22

DNA extraction was performed at the Instituto Nacional de Saúde laboratories in Maputo, Mozambique. DNA samples were then shipped to the Molecular Biology laboratories at the Department of Laboratory Medicine, Faculty of Medicine and Health, Örebro University Hospital, Örebro, Sweden for HPV genotyping.

HIV testing was performed as previously described on site in Maputo.20

Statistical analysis

Demographic and risk behavioural data were entered into a MySQL database V.5.1 (MySQL AB, 2008) and HPV laboratory results in a Microsoft Office Excel 2010 spreadsheet. For processing and analysis, data were imported into Stata V.14 (StataCorp. 2015. Stata: Release 14. Statistical Software. College Station, TX: StataCorp LP).

All analyses were stratified by gender. Descriptive statistics were employed as follows: frequencies for categorical variables, and means and SD for quantitative variables. We used list-wise deletion as the missing values were negligible. Here, we report the prevalence of HPV infections and HPV genotypes per gender and HIV status. Bivariate logistic analysis between sociodemographic and sexual behavioural characteristics and presence of any HPV infection was conducted. Age of individuals and age of sexual debut and characteristics whose p value on bivariate logistic regression analysis was below 0.25 were included in the multivariable logistic regression. Additionally, to identify characteristics that would contribute to multiple HPV genotype infections a multinomial logistic regression (MNLR) was used. For this analysis, males and females are not stratified given the sample constraints. We first grouped the participants as follows: (a) subjects without an HPV infection; (b) subjects with one HPV infection (monoinfection); and (c) subjects with at least two HPV infections (multiple infections). The MNLR per level of the factor under consideration provided two ORs (one for monoinfected versus not-infected subjects and another for the subjects with multiple-infections versus not-infected subjects). We report both the unadjusted and adjusted ORs and their respective 95% CI. For adjustment, we included factors whose p values through the global likelihood ratio test were below 0.25 in the MNLR bivariate analysis; gender and age were forced factors to include in the model. The significance level was set at 5%.

Ethics statement

Ethical approval was granted by the National Health Bioethics Committee of Mozambique (Ref. 148/CNBS of 8 May 2009 and Ref. 18/CNBS/11). Testing performed in Sweden was approved by the Regional Ethical Review Board Uppsala. Study investigators followed the Good Clinical Practice-International Conference on Harmonisation (GCP-ICH) guidelines. Subjects signed an informed consent form prior to any study activities.

Results

In total, 489 participants (263 females and 226 males) were enrolled, but samples from 77 (15.7%) subjects had undetectable DNA levels and were excluded from the data analysis, being 27/263 (10.3%) from females and 50/226 (22.1%) from males.

Demographic and behavioural characteristics

Table 1 shows the socio demographic and sexual behavioural characteristics of study participants by sex. Of the 412 subjects with valid HPV results, 236 (57.3%) were females, and 176 (42.7%) were males. The mean age of study participants was 21.1 years (SD: 1.71) with males being slightly older than females (21.5 vs 20.8). Almost all participants were students (94.2%) and all had some formal education. Male subjects were more educated than females with more than half having secondary or higher educational degrees (56.8% vs 41.9%). Nearly all participants were single (98.8%). The median age at sexual debut was 17 years (IQR: 15–18) being 16 years (IQR: 14–18) for males and 17 years (IQR: 16–18) for females. The majority of subjects (89.3%) had more than one sexual partner in life (97.2% of males and 83.5% of females) and 24% reported two or more sexual partners in the last 6 months (37.5% of males and 14% of females). Around one-fourth of the study population reported having at least one episode of STI in life and this was seen more frequently in women than in men (34.3% vs 13.6%). Only 66.0% of subjects reported having used a condom at the time of their last sexual intercourse with females reporting a lower frequency compared with males (59.7% vs 74.4%). HIV infection was diagnosed in 21 (5.1%) participants, of which 17 (81.0%) were females.

Table 1

Sociodemographic and sexual behavioural characteristics of study participants by sex

	Total		Male		Female
	n	%	n	%	n	%
Total	412	100.0	176	100.0	236	100.0
Age (years)
Mean (SD)	21.1 (1.71)		21.5 (1.67)		20.8 (1.68)
Marital status
Single	407	98.8	175	99.4	232	98.3
Married/cohabitating	5	1.2	1	0.6	4	1.7
Education
Secondary or less	213	51.7	76	43.2	137	58.1
Technical training	109	26.5	50	28.4	59	25.0
University grade	90	21.8	50	28.4	40	16.9
Occupation
Student	388	94.2	163	92.6	225	95.3
Employed	23	5.6	12	6.8	11	4.7
Missing	1	0.2	1	0.6	0	0.0
Religion
Christian	375	91.0	162	92.0	213	90.3
Other	37	9.0	14	8.0	23	9.7
Age at sexual debut
Less than 18	285	69.2	115	65.3	170	72.0
18 or more	122	29.6	58	33.0	64	27.1
Missing	5	1.2	3	1.7	2	0.8
Number of sex partners in life
0–1	44	10.7	5	2.8	39	16.5
>1	368	89.3	171	97.2	197	83.5
Number of sex partners in the last 6 months
0–1	313	76.0	110	62.5	203	86.0
>1	99	24.0	66	37.5	33	14.0
Condom use in the last sexual intercourse
No	139	33.7	45	25.6	94	39.8
Yes	272	66.0	131	74.4	141	59.7
Missing	1	0.2	0	0.0	1	0.4
Had a STI in life
No	305	74.0	150	85.2	155	65.7
Yes	105	25.5	24	13.6	81	34.3
Missing	2	0.5	2	1.1	0	0.0
HIV infection
Positive	21	5.1	4	2.3	17	7.2
Negative	391	94.9	172	97.7	219	92.8

STI, sexually transmitted infection.

HPV prevalence and genotype distributions

The overall HPV prevalence was 40.8% (168 infected subjects; 95% CI 36.0 to 45.5) and was significantly higher in women than in men (150/236 [63.6%, 95% CI 57.1 to 69.7] vs 18/176 [10.2%, 95% CI 6.2 to 15.7], p<0.001). The HPV prevalence among HIV-infected subjects (14/21, 66.7%) was 69% higher than in the HIV-negative subjects (154/391, 39.4%; p=0.013). Figure 1 and table 2 show the distribution of HPV genotypes in male and female participants. Of the 412 enrolled participants, 115 (27.9%), 54 (13.1%) and 89 (21.6%) were infected with one or more high-risk, probably high-risk and low-risk HPV types, respectively. High-risk and phr-HPV types were present in 124/150 (82.7%, 95% CI 75.6 to 88.4) and 14/18 (77.8%, 95% CI 52.4 to 93.6) of HPV-infected women and men, respectively (p=0.533). Overall, HPV52 was the most frequent type found (9.1%), followed by HPV35, −6, –16, −53, –58 and −51. Altogether, these HPV types accounted for 49.9% of the infections in women and men. In female participants, the most common genotypes were similar to the overall distribution (HPV52, −35, –16, −53, –58, −6, and −51) and accounted for half (50.3%) of the infections in women (figure 1A). In contrast, male participants’ HPV infections were less diverse with HPV51 being the most frequent type found, followed by HPV6, −11, –52, −59 and −70 in equal proportions. These six types were responsible for 52% of the infections in men (figure 1B). Hr-HPV genotypes were identified in 11/21 (52.4%) of HIV-infected subjects. Three hr-HPVs (HPV35, −52 and −58) accounted for 60% of all hr-HPV infections in this population. There were no differences in the number of subjects infected with hr-HPVs and phr-HPVs in the HIV-infected (p=0.552) and uninfected populations (p=1.000). In contrast, lr-HPVs were more frequently observed in the HIV-infected subjects (p=0.022).

Figure 1

HPV genotype distribution in female (A) and male (B) participants. For calculation of proportions, the denominator used was the number of HPV-infected participants per gender. HPV, human papillomavirus.

Table 2

HPV prevalence and HPV genotype distribution by sex

	Total	Gender		HIV serostatus		HPV genotype coverage
		Male	Female	Negative	Positive	Gardasil	Cervarix	Gardasil 9
Total subjects, n	412	176	236	391	21
Number of HPV infections, n	351	25	326	305	46
Subjects with at least 1 HPV infection*, n (%)	168 (40.8)	18 (10.2)	150 (63.6)	154 (39.4)	14 (66.7)
Subjects with hr-HPV infections†, n (%)	115 (27.9)	10 (5.7)	105 (44.5)	104 (26.6)	11 (52.4)
16	24 (5.8)	1 (0.6)	23 (9.7)	22 (5.6)	2 (9.5)	++++	++++	++++
18	11 (2.7)	0 (0.0)	11 (4.7)	10 (2.6)	1 (4.8)	++++	++++	++++
31	6 (1.5)	1 (0.6)	5 (2.1)	5 (1.3)	1 (4.8)	++++ρ	++++§	++++
33	4 (1.0)	0 (0.0)	4 (1.7)	2 (0.5)	2 (9.5)	++++ρ	++++ρ	++++
35	27 (6.6)	1 (0.6)	26 (11.0)	22 (5.6)	5 (23.8)
39	4 (1.0)	0 (0.0)	4 (1.7)	4 (1.0)	0 (0.0)
45	2 (0.5)	0 (0.0)	2 (0.8)	1 (0.3)	1 (4.8)		++++§	++++
51	20 (4.9)	3 (1.7)	17 (7.2)	19 (4.9)	1 (4.8)
52	32 (7.8)	2 (1.1)	30 (12.7)	28 (7.2)	4 (19.0)			++++
56	4 (1.0)	1 (0.6)	3 (1.3)	4 (1.0)	0 (0.0)
58	24 (5.8)	1 (0.6)	23 (9.7)	21 (5.4)	3 (14.3)			++++
59	13 (3.2)	2 (1.1)	11 (4.7)	13 (3.3)	0 (0.0)
Total number of infections‡	171 (48.7)	12 (48.0)	159 (48.8)	151 (49.5)	20 (43.5)
Subjects with probably or possibly hr-HPV infections, n (%)	54 (13.1)	5 (2.8)	49 (20.8)	50 (12.8)	4 (19.0)
26	1 (0.2)	0 (0.0)	1 (0.4)	1 (0.3)	0 (0.0)
53	24 (5.8)	1 (0.6)	23 (9.7)	22 (5.6)	2 (9.5)
66	10 (2.4)	1 (0.6)	9 (3.8)	9 (2.3)	1 (4.8)
68	3 (0.7)	0 (0.0)	3 (1.3)	3 (0.8)	0 (0.0)
70	12 (2.9)	2 (1.1)	10 (4.2)	11 (2.8)	1 (4.8)
73	2 (0.5)	0 (0.0)	2 (0.8)	2 (0.5)	0 (0.0)
82	6 (1.5)	1 (0.6)	5 (2.1)	6 (1.5)	0 (0.0)
85	0 (0.0)	0 (0.0)	0 (0.0)	0 (0.0)	0 (0.0)
Total number of infections‡	58 (16.5)	5 (20.0)	53 (16.3)	54 (17.7)	4 (8.7)
Subjects with lr-HPV infections, n (%)	89 (21.6)	7 (4.0)	82 (34.7)	79 (20.2)	12 (57.1)
6	24 (5.8)	2 (1.1)	22 (9.3)	21 (5.4)	3 (14.3)	++++		++++
11	8 (1.9)	2 (1.1)	6 (2.5)	6 (1.5)	2 (9.5)	++++		++++
40	6 (1.5)	0 (0.0)	6 (2.5)	5 (1.3)	1 (4.8)
42	11 (2.7)	0 (0.0)	11 (4.7)	10 (2.6)	1 (4.8)
43	1 (0.2)	1 (0.6)	0 (0.0)	1 (0.3)	0 (0.0)
44	3 (0.7)	0 (0.0)	3 (1.3)	3 (0.8)	0 (0.0)
54	7 (1.7)	0 (0.0)	7 (3.0)	4 (1.0)	3 (14.3)
61	13 (3.2)	0 (0.0)	13 (5.5)	11 (2.8)	2 (9.5)
62	14 (3.4)	1 (0.6)	13 (5.5)	11 (2.8)	3 (14.3)
71	2 (0.5)	0 (0.0)	2 (0.8)	2 (0.5)	0 (0.0)
72	2 (0.5)	0 (0.0)	2 (0.8)	1 (0.3)	1 (4.8)
81	10 (2.4)	0 (0.0)	10 (4.2)	7 (1.8)	3 (14.3)
83	8 (1.9)	1 (0.6)	7 (3.0)	7 (1.8)	1 (4.8)
84	12 (2.9)	1 (0.6)	11 (4.7)	10 (2.6)	2 (9.5)
89	1 (0.2)	0 (0.0)	1 (0.4)	1 (0.3)	0 (0.0)
Total number of infections‡	122 (34.8)	8 (32.0)	114 (35.0)	100 (32.8)	22 (47.8)

Genotype covered by the vaccine.

Genotype highly covered through vaccine crossprotection.

Genotype moderately covered through vaccine crossprotection.

*The denominator is total subjects.

†The denominator is total subjects.

‡The denominator is number of HPV infections.

HPV, human papillomavirus; lr-HPV, low-risk HPV; hr-HPV, high-risk HPV.

Vaccine-matched HPV genotypes

Table 2 shows the distribution of HPV genotypes in women and men, and their correspondence to the genotypes present in the current HPV vaccines. Genotypes associated with vaccine crossprotection were included in the analysis and considered as covered by the vaccine. The total number of circulating HPV genotypes matching the vaccine genotypes were 5/33 (15.2%; Gardasil: HPV6,−11, –16, −18 and −31), 5/33 (15.2%; Cervarix: HPV16, −18, –31, −33 and −45) and 9/33 (27.3%; Gardasil 9: HPV6, −11, –16, −18, –31, −33, –45, −52 and −58) for the female population and 4/18 (22.2%; Gardasil: HPV6, −11, –16 and −31), 2/18 (11.1%; Cervarix: HPV6 and −31) and 6/18 (33.3%; Gardasil 9: HPV6, −11, –16, −31, –52 and −58) for the male population. The three vaccines can cover 3/19 (15.8%; Gardasil: HPV16, −18 and −31), 5/19 (26.3%; Cervarix: HPV16, −18, –31, −33 and −45) and 7/19 (36.8%; Gardasil 9: HPV16, −18, –31, −33, –45, −52 and −58) hr-HPVs circulating in women and 2/12 (16.7%; Gardasil: HPV16 and −31), 2/12 (16.7%; Cervarix: HPV16 and −31) and 4/12 (33.3%; Gardasil 9: HPV16, −31, –52 and −58) in men.

Factors associated with HPV infection

Overall, women had a significantly higher HPV infection prevalence than men (p<0.001). Although the odds of being infected with HPV reduced by 80% for each year of age, in the multivariate analysis this was not significant. Sexual debut before the age of 18 years, history of STI and infection with HIV were significantly associated with the presence of HPV infection (p=0.008, p<0.001, p=0.013, respectively) in the univariate analysis. When stratifying by gender (table 3A,B), the univariate analysis show that women who initiated sexual activity before the age of 18 were significantly more at risk of having a HPV infection (p=0.041), and a marginally significance association was seen in women who reported to have had more than two sexual partners in the last 6 months (p=0.055), but this was not seen in the multivariate analysis. No significant associations were seen for male participants both in the univariate and multivariate analysis, but a marginally significant association was seen for a reported history of STI and HPV infection.

Table 3

Factors associated with any HPV infection: unadjusted and adjusted

	Any HPV infection		Univariate			Multivariate
	Negative (n)	Positive (n)	OR	95% CI	p	OR	95% CI	p
Among male participants
Total observations	158	18				171
Age (years)			0.90*	0.67 to 1.20	0.465	0.92*	0.66 to 1.27	0.605
Education					0.173			0.239
Primary and secondary	66	10	1.00	–		1.00	–
Technical training	44	6	0.90	0.31 to 2.65		1.02	0.31 to 3.34
University grade	48	2	0.28	0.06 to 1.31		0.27	0.05 to 1.32
Occupation					0.813
Student	146	17	1.00	–
Employed	11	1	0.78	0.09 to 6.43
Religion					0.618
Christian	146	16	1.00	–
Other	12	2	1.52	0.31 to 7.41
Age at sexual debut					0.268			0.191
Less than 18	101	14	1.00	–		1.00	–
18 or more	54	4	0.53	0.17 to 1.70		0.43	0.12 to 1.52
Number of sex partners in the last 6 months					0.254
0–1	101	9	1.00	–
>1	57	9	1.77	0.67 to 4.72
Condom use in the last sexual intercourse					0.822
No	40	5	1.00	–
Yes	118	13	0.88	0.30 to 2.63
Had a STI in life					0.097			0.057
No	137	13	1.00	–		1.00	–
Yes	19	5	2.77	0.89 to 8.65		3.21	0.97 to 10.64
HIV infection					0.347
Positive	3	1	3.04	0.30 to 30.86
Negative	155	17	1.00	–
Among female participants
Total observations	86	150				234
Age (years)			0.87*	0.74 to 1.02	0.091	0.88*	0.75 to 1.04	0.139
Education					0.246			0.267
Primary and secondary	51	86	1.00	–		1.00	–
Technical training	17	42	1.47	0.76 to 2.84		1.57	0.80 to 3.11
University grade	18	22	0.72	0.36 to 1.48		0.81	0.39 to 1.68
Occupation					0.521
Student	83	142	1.00	–
Employed	3	8	1.56	0.40 to 6.04
Religion					0.862
Christian	78	135	1.00	–
Other	8	15	1.08	0.44 to 2.67
Age at sexual debut					0.041			0.109
Less than 18	55	115	1.00	–		1.00	–
18 or more	30	34	0.54	0.30 to 0.97		0.61	0.33 to 1.12
Number of sex partners in the last 6 months					0.055			0.150
0–1	79	124	1.00	–		1.00	–
>1	7	26	2.37	0.98 to 5.71		1.94	0.79 to 4.80
Condom use in the last sexual intercourse					1.000
No	34	60	1.00	–
Yes	51	90	1.00	0.58 to 1.72
Had a STI in life					0.665
No	58	97	1.00	–
Yes	28	53	1.13	0.65 to 1.98
HIV infection					0.258
Positive	4	13	1.95	0.61 to 6.17
Negative	82	137	1.00	–

*Associated effect of 1-year increase.

HPV, human papillomavirus.

HPV monoinfection versus multiple infections

Almost half, 75/168 (44.6%), of HPV infected subjects had a single infection, 48/168 (28.6%) had two concurrent HPV infections and 45/168 (26.8%) had three or more HPV coinfections (figure 2). HPV-infected male participants had more frequently a single infection (72.2% of the cases), whereas more than 50% of HPV-infected females had two or more HPV coinfections. In the univariate analysis, several factors were shown to be associated with the presence of two or more (multiple) HPV coinfections. Being a woman, younger age (the odds of being infected with HPV reduced by 80% for each year of age), lower educational degree, sexual debut before the age of 18, reported STI in life and the presence of a HIV infection were factors associated with multiple-HPV infections. Nonetheless, in the multivariate analysis only three factors were demonstrated to be significantly associated, namely, being a woman (p=0.001), having started sexual activity before the age of 18 (p=0.008) and the presence of a HIV infection (p=0.003). Condom use in the last sexual intercourse was not significantly associated with protection against HPV infection (p=0.143) (table 4).

Figure 2

Distribution of HPV infections in female and male participants. The denominator is the total number of infected women and men, respectively. HPV, human papillomavirus.

Table 4

Factors associated with the number (mono or multiple) of HPV infection: unadjusted and adjusted using multinominal logistic regression

	Univariate				p	Multivariate				p
	Monoinfection vs negative		Multiple infections vs negative			Monoinfection vs negative		Multiple infections vs negative
	OR	95% CI	OR	95% CI		OR	95% CI	OR	95% CI
Total observations	412					406
Gender
Male	1.00	–	1.00	–	<0.001	1.00	–	1.00	–	<0.001
Female	8.76	4.56 to 16.83	32.33	12.65 to 82.66		7.23	3.64 to 14.33	28.38	10.78 to 74.70
Age (years)	0.80*	0.69 to 0.94	0.80*	0.69 to 0.92	0.001	0.85*	0.71 to 1.02	0.87*	0.72 to 1.04	0.134
Marital status
Single	1.00	–	1.00	–	0.632
Married/cohabitating	1.64	0.15 to 18.29	2.66	0.37 to 19.16
Education
Secondary or less	1.00	–	1.00	–	0.040	1.00	–	1.00	–	0.246
Technical training	0.98	0.54 to 1.79	0.94	0.54 to 1.65		1.42	0.71 to 2.83	1.77	0.87 to 3.61
University grade	0.41	0.19 to 0.87	0.47	0.24 to 0.91		0.59	0.26 to 1.33	0.84	0.38 to 1.85
Occupation
Student	1.00	–	1.00	–	0.784
Employed	1.17	0.41 to 3.36	0.74	0.24 to 2.29
Religion
Christian	1.00	–	1.00	–	0.764
Other	1.15	0.47 to 2.84	1.35	0.61 to 3.00
Age at sexual debut
Less than 18	1.00	–	1.00	–	0.003	1.00	–	1.00	–	0.008
18 or more	0.82	0.47 to 1.44	0.36	0.20 to 0.67		0.91	0.48 to 1.72	0.33	0.16 to 0.67
Number of sex partners in life
0–1	1.00	–	1.00	–	0.424
>1	0.64	0.29 to 1.43	0.67	0.32 to 1.41
Number of sex partners in the last 6 months
0–1	1.00	–	1.00	–	0.369
>1	0.65	0.34 to 1.23	0.82	0.47 to 1.44
Condom use in the last sexual intercourse
No	1.00	–	1.00	–	0.143	1.00	–	1.00	–	0.423
Yes	0.59	0.35 to 1.00	0.80	0.48 to 1.32		0.81	0.45 to 1.48	1.25	0.68 to 2.31
Had a STI in life
No	1.00	–	1.00	–	0.003	1.00	–	1.00	–	0.530
Yes	2.33	1.32 to 4.12	2.07	1.21 to 3.55		1.44	0.76 to 2.72	1.18	0.62 to 2.24
HIV infection
Positive	0.46	0.06 to 3.79	5.50	2.12 to 14.27	<0.001	0.44	0.05 to 4.04	6.03	1.73 to 21.02	0.003
Negative	1.00	–	1.00	–		1.00	–	1.00	–

*Associated effect of 1-year increase.

HPV, human papillomavirus.

Discussion

This was a cross-sectional study reporting the prevalence of HPV infections and genotype distributions in young adults aged 18–24 in Maputo city, Mozambique. The HPV prevalence in the female population was 63.6% (95% CI 57.1 to 69.7) and was similar to previous reports of 75.9% in women aged 14–61 years in rural southern Mozambique and other neighbouring countries.11–13 23 The HPV prevalence in the male population was 10.2% (95% CI 6.1 to 15.7), which was lower than previous reports for African countries and other regions.2 Although no significant difference was observed in the number of hr-HPV infections between women and men, a high susceptibility of women to HPV infections was demonstrated.

We report a lower prevalence of HPV infections in men than in women, which is contrary to what has been shown in neighbouring countries and other regions.24 Studies assessing different male anatomic collection sites show that HPV detection in urethral samples is suboptimal compared with other sampling sites in the male genitalia, but the sensitivity reported for urethral sampling has varied.25 26 In a study comparing sampling from seven male genital sites, Giuliano et al reported that detection was highest at the penile shaft (48%), while only 10% of HPV infections were detected when using urethral sampling.26 In contrast, the sensitivity of urethral samples for overall HPV detection was 54.2% in a study that explored expressed prostate secretions as an option for detection of HPV infections in Russian men.27 We found 50/226 (22.1%) urethral samples with undetectable DNA and 158/176 (89.8%) of valid samples tested negative. The first could be related to an inadequate sampling technique and the second to the selected anatomic collection site. Altogether, this could have contributed to the lower male HPV prevalence found in our study.

Worldwide, the prevalence of cervical HPV infections in healthy women has been demonstrated to be 11.7%. Nonetheless, in less developed regions, such as East Africa, the prevalence is three times the global figures (33.6%).28 High-risk HPV16, −18, –52, −31, –58, −39, –51 and −56 have been found in women with normal cytology. In this study, we demonstrated that hr-HPV genotype distribution in Mozambican healthy young women differs from the global figures with HPV52, −35, –16, −53, –58 and −51 being the most prevalent. Castellsagué et al have previously described that hr-HPVs 51, –35, −18, –31 and −52 were the most commonly found in women aged 14–61 years with normal cytology in a study conducted in southern Mozambique,13 which is similar to our findings. Furthermore, HPV16 contributed to only 6.8% of female infections in the present study. These findings are concordant with those from Castellsagué et al.13 Globally, HPV16 has been responsible for 22.5% of all infections in women, whereas in sub-Saharan Africa, it seems to have the lowest contribution (13.7%, 11.3% and 11.1% for Southern, Eastern and Western Africa, respectively).28 The peak prevalence of HPV infections in women was shown to occur in early ages (before the age of 30 years).11 12 Some reports have also described a less pronounced second peak prevalence in older women (around their 50s). Reasons for the occurrence of the first HPV peak are related to absence of immune responses to HPV and sexual activity in younger women. The second peak is not consistently reported across the globe, but it is speculated that hormonal changes my play a role. Although the majority of HPV infections in early ages clear spontaneously,29 30 some do become chronic infections and may potentially evolve to malignant transformation of the cervix. In the present study, we described the epidemiology of HPV infections in both women and men that belong to the ‘first HPV peak’ population.

Our report shows that more than 50% of the study population had two or more concomitant HPV infections with more than 50% of female participants having two or more coinfections compared with only 27.8% of male participants. It has been demonstrated by others that women with multiple HPV infections take longer to clear their infections compared with women with a single infection,29 30 thus potentially contributing to the establishment of a persistent infection, a prerequisite for cervical cancer development. Several studies conducted in Mozambique have shown a strong association between HPV16, −33, –35, −45, –18, −31 and −58 with cervical neoplasia.15–18 In another Mozambican study, HPV52 was reported to be the second most common hr-HPV in women with abnormal cytology and the fourth most common in women with invasive cervical cancer.13 In the present study, we demonstrated a 50% match between the hr-HPV types commonly circulating in the female population (HPV35, −16, –52 and −58) and HPVs types commonly present in cervical neoplasia in Maputo, Mozambique.

We confirmed the association between early onset of sexual debut, history of past STI and HIV infection with multiple HPV infections.31 32 These findings reiterate the need to intensify sexual behaviour education and STI prevention in young populations.

This study shows a high prevalence of HPV infections (66.7%) in the HIV-infected population. A previous report showed a prevalence of HPV infections of 56.6% in HIV-infected African women, which is concordant with our findings.33 Other reports from South Africa have demonstrated an even higher HPV prevalence (75%) in the age group of 17–19 years.34 We demonstrated that HPV infections and multiple HPV infections are more frequent in HIV-infected subjects. Similar descriptions have been previously reported.34 35 We have shown that HPV35, −52 and −58 were responsible for 60% of all hr-HPV infections in this population. McDonald et al have reported that hrHPV35, −58, –18, −45, –16 and −52 were the most common genotypes in HIV-infected women with normal cytology in neighbouring South Africa.34 Lastly, other reports in Africa have shown that the most frequent genotypes circulating in HIV-infected women were HPV16, −58, –52, −31, –35 and −18.33 These findings suggest a high homology between the common hr-HPV genotypes in Mozambique versus South Africa and other African regions. We have also described a high prevalence of lr-HPV infections in the HIV-infected subjects. These findings were also demonstrated elsewhere.36

Since 2014, Mozambique has been engaged in the HPV vaccine demonstration project which aimed at accessing the country’s preparedness for the introduction of an HPV vaccine in the national vaccination programme. Information regarding HPV types circulating in young populations is crucial for selection of a suitable vaccine. Currently available HPV vaccines, Gardasil, Cervarix and Gardasil 9 can prevent up to seven hr-HPV infections. Moderate to high crossprotection against HPV types phylogenetically related to HPV16 and −18 have also been reported with the use of Gardasil and Cervarix. Nonetheless, the level of protection has shown to be heterogeneous across the studies and a signal for reduction of vaccine efficacy with extended follow-up has also been demonstrated.37 Our study shows that Gardasil 9 has the highest genotype coverage in both Mozambican women and men, whereas Gardasil has the lowest in women and Cervarix in men. Gardasil 9 can provide protection against HPV52, the most common genotype found in this study. Nonetheless, Gardasil 9 would only cover 37% of the high-risk genotypes found. This has important implications for vaccine-strategic discussions. It should also send signals to vaccine manufacturers concerning future developments of new versions of the HPV vaccines. In the present study, we did not assess HPV vaccination status since no HPV vaccine had been introduced or administered within the National Health System.

This study has some limitations. Power calculation was applied for the HIV incidence in the parent study but not for the HPV ancillary study. We did not include an older cohort (aged ≥25 years) for comparison of circulating HPV genotypes due to funding restrictions. This must be considered in future studies to better describe the epidemiology of HPV infections in Mozambique. In our study, the questionnaire assessing sexual behavioural captured information about the number of sexual partners categorised into two variables only, either having one or more than one partner. This did not enable us to further analyse the influence of the number of sexual partners on HPV status. The low-HPV prevalence in the male population may have been related to the site and technique of sampling and we therefore recommend that additional studies in male populations, in Mozambique, are considered.

Conclusion

This study confirmed the high burden of HPV infections in young women in the Maputo city, Mozambique. The HPV prevalence was demonstrated to be associated with high-risk sexual behaviour, thus confirming the need to intensify sex education and STI prevention interventions. Additional studies involving other anatomic collection sites in men to confirm the comparatively low prevalence of HPV are needed. For vaccine-strategic discussions, it is important to consider that the current HPV vaccines only cover a proportion of prevailing HPV types in this area.