Human papillomavirus infection and seroprevalence among female university students in Mexico

ABSTRACT

Human papillomavirus (HPV) is one of the most common causes of sexually transmitted diseases, and the main etiology of cervical cancer. This study was aimed to assess type-specific cervical HPV prevalence and their association with HPV-specific antibodies in a cohort of female university students. HPV genotyping was performed by amplifying and sequencing a fragment of the L1 protein. A BLAST search was performed to identify HPV types. HPV-specific IgG antibodies were measured by ELISA in serum samples. A total of 129 women participated, with an average age of 21.75 years. The prevalence of vaginal HPV infection was 74.42%. The most predominant high-risk HPV types were 18 (13.95%), 31 (10.85%), and 16 (9.3%). We found that early age at coitarche and a higher number of sexual partners were significantly associated with a high prevalence of HPV infection. In addition to sexual behavior, we observed that the presence of serum-specific IgG antibodies against HPV can impact the prevalence of the virus. Seropositivity to HPV-16 and HPV-18 was associated with a lower prevalence of HPV-16, but not for other HPV types. Of note, there was a lower proportion of HPV-specific seropositivity in women who had the presence of the same HPV type in a cervical specimen, suggesting an immunoregulatory mechanism associated with the viral infection. In conclusion, the prevalence of HPV in university women was higher than expected and it was associated with early age of sexual debut, an increasing number of sexual partners, and a low proportion of HPV seropositivity.

KEYWORDS:

• Human papilloma virus
• prophylactic HPV vaccination
• cervical cancer
• México
• seroprevalence

Introduction

Human Papillomavirus (HPV) infection is among the most common sexually transmitted infections. Papillomaviruses are small non-enveloped viruses with around 8 kb of circular double-stranded DNA genomes including a non-coding regulatory long control region (LCR), and eight open reading frames (ORF): L1 and L2 that encodes capsid proteins; E1, E2, E4, E5, E6, and E7 that encodes proteins involved in replication, transcription, and transformation.¹ HPV types are established when the DNA sequence of the L1 differs by at least 10%. More than 170 types of HPV have been identified; approximately 40 of them infect the epithelium in the genital tract, the mucosa of the upper respiratory tract, and the skin developing epidermodysplasia verruciformis.^2,3 HPV infection has a central role in common dermatologic and sexually transmitted diseases as well as in one of the most frequent cancers worldwide.^4,5 The risk of being infected by HPV at least once in a lifetime among both men and women is calculated as 50%.⁶ However, in the majority of cases, HPV infection is transient or asymptomatic and is resolved spontaneously.⁷ According to their oncogenic potential, HPV are classified as low-risk (wart-causing) and high-risk (cancer-causing) viruses.⁷ Persistent infection with high-risk types is associated with precancerous and cancerous lesions.⁸ High-risk HPV types 16, 18, 31, 33, and 35 have been associated with cervical, vulvar, vaginal, penile, anal, and oropharyngeal cancers and with pre-cancer lesions.^9,10 HPV types 16 and 18 are the most common high-risk types worldwide and are considered to be responsible for more than 70% of all cervical cancer cases.^11,12 In Mexico, previous studies have identified HPV types 16, 18, 31, 45, 51, 58, and 59 as the most prevalent in cervical samples.^13–16

Prophylactic HPV vaccination was introduced in 2006 with the development of two vaccines: the bivalent vaccine (Cervarix) containing HPV 16 and 18 antigens, and the quadrivalent vaccine (Gardasil) containing HPV 6, 11, 16, and 18 antigens. Recently, a nonavalent vaccine (Gardasil 9) was approved, protecting against infection from HPV types 6, 11, 16, 18, 31, 33, 45, 52, and 58. These vaccines are made with L1 viral proteins that assemble into morphologically similar particles to the virus but do not contain its genome, the so-called virus-like particles (VLP). Different studies have demonstrated the efficiency of HPV vaccination in preventing infection with the types of HPV they target when given before initial exposure to the virus. High-risk HPV-16 and 18 are virtually preventable by vaccination.^17,18 In Mexico, vaccination against HPV was first introduced in 2008 with low coverage to girls aged 12–16 years using a 0–6 month schedule; 1 year later, it changed to an extended dosing schedule that targets girls aged 9–12 years for the first 2 doses, applied 6 months apart, followed by a third dose 60 months later.¹⁹ The vaccine was included in the national vaccine program until 2012.⁹ The coverage has been increased over time; according to the last reported data in 2018, about 1 million doses were applied in Mexico.²⁰ However, this number is still very low considering the whole population of 125 million persons, from who 5.7% are females between 9 and 14 years old.²¹

Despite the demonstrated benefits of HPV immunization, there are still many obstacles to implement the vaccination schedules in high- and low-income countries. Those barriers are multifactorial and include limitations in costs, infrastructure, and even social stigma.²² Due to these difficulties in less developed countries the prevalence of HPV infections remains high. Mexico is within the regions with a high rate of HPV infection.^9,12 Several studies in different regions have shown the highest prevalence of HPV in women younger than 25 years of age and college women because of their risky sexual behavior, lack of knowledge of HPV infection, low rate of vaccination, and other cultural factors.^23–25 Because of the need for data that support the actual status of the prevalence of HPV infection in our country as well as the levels of specific antibodies in individuals at high risk of HPV infection the current study was designed to identify HPV prevalence, genotypes, and seropositivity to help establish appropriate preventive strategies for a high-risk population.

Material and methods

Study population

One hundred twenty-nine female students who attended gynecological consultation at the Clinica Universitaria de Salud Integral Iztacala (CUSI) from March 2017 to April 2018 were included. All the subjects that provide informed consent were eligible for participation. This study was conducted with the approval of the Institutional Review Board of Facultad de Estudios Superiores Iztacala (FESI). The clinical data for the study was extracted from the clinical history of each individual.

Preparation of DNA from cervical cell swabs

Porous foam-tipped swabs (Epicenter Biotechnologies, Madison, WI, USA) were used for the collection of cervical exfoliated epithelial cells from patients. Cervical DNA samples were extracted with DNA Extraction Kits (Epicenter Biotechnologies). Samples were processed according to the manufacturer’s instructions.

DNA genotyping

Viral DNA was amplified with polymerase-chain reaction in the L1 region of the HPV genome, with human beta-actin as the reference gene. The sequences of primers used in this were three consensus primer pairs, L1C1 ⁄ L1C2 (Fw: 5ʹCGTAAACGTTTTCCCTATTTTTTT 3’, Rev 5ʹTACCCTAAATACTCTGTATTG 3’), GP5+ ⁄ 6+ (5’ TTTGTTACTGTGGTAGATACTAC3’, Rev 5’ GAAAAATAAACTGTAAATCATATTC 3’, and actin (5’ GCACAGGGACATAATAATGG 3’, Rev 5’ CGTCCAAAAGGAAACTGATC 3’). Polymerase chain reaction amplification was done in a 50 μL reaction mixture. Thermal cycling conditions were as follows: denaturation at 95°C for 10 minutes followed by 40 cycles of denaturation at 95°C for 45 seconds, annealing at 50°C for 55 seconds, and extension at 72°C for 1 minute. The amplification program was followed by a final extension step at 72°C for ten minutes. The PCR products were electrophoretically migrated on 1.5% agarose gels and the bands corresponding to the specific amplicons were cut off from gels. Then DNA was extracted using the WizardRSV Gel and PCR Clean-Up System. Every product was used as a template to conduct a sequencing reaction with the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Austin, TX, USA), according to the manufacturer’s procedures. The sequences were analyzed using the ABI PRISM 310 system (Applied Biosystems). The sequence data obtained by automated DNA sequencing were analyzed using the BLAST search (http://www.ncbi.nlm.nih.gov/ BLAST/, accessed in 2017–2018) for HPV genotypes. DNA sequencing identified the most dominant genotype in a given specimen, and it served as a confirmatory assay.

HPV-negative detected samples were confirmed by qPCR using the consensus primer set GP5+/GP6+. All qPCR was carried out using an StepOne Real-Time PCR System (Applied Biosystems) with Maxima SYBR Green/ROX qPCR Kit (Thermo Scientific, USA). Reactions comprised 1× SYBR Green mix, 50 nM primer pairs, and 5 μl of template. The following cycling conditions were employed: initial denaturation 95°C 15 minutes, followed by 40 cycles of 95°C 15 seconds, 46°C 1 minute, 72°C 1 minute. A final extension of 72°C 10 minutes, and melting curve of 95°C 15 seconds, 60°C 1 minute, 95°C 15 seconds transition were incorporated.

HPV ELISA

Antibodies against HPV were measured by a VLP-based direct enzyme-linked immunoabsorbent assay (ELISA) base on the recommendations of the World Health Organization.²⁶ Maxisorp Nunc-Immuno plates (Thermo Fisher Scientific, MA, USA) were coated with 1 µg/ml of either purified recombinant HPV-16 (Abcam, Cat. ab119880) or HPV-18 (Abcam, Cat. ab119881) L1 proteins in phosphate-buffered saline (PBS), pH 7.2 at 4°C overnight. The plates were blocked with ELISA Assay Diluent (Biolegend, Cat. 421203) at 4°C overnight. Plates were incubated with serum samples diluted 1:1000 in ELISA Assay Diluent for 2 hrs at 4°C. After washing the wells 5 times with PBS containing 0.1% Tween 20 (Sigma-Aldrich), HRP-conjugated goat anti-human IgG (Abcam, Cat. ab6858) diluted 1:7500 was added and incubated at room temperature for 1 hr. Following an additional washing cycle, TMB ELISA Substrate (Abcam, Cat. ab171527) was added and incubated at room temperature for 15 minutes, followed by the addition of 1 M H₂SO₄. The absorbance was measured at 490 nm. To control the inter, and intra-assay variation it was included a positive control serum, identified during the optimization process of the immunoassay, diluted in the same way as the specimens in every ELISA plate. For data analysis, it was first calculated the net mean OD value for each sample by subtracting the mean blank OD. Then a normalized absorbance ratio (NAR) was obtained by dividing the net mean OD value of each sample by the net mean OD from the positive control serum included in the same plate, similar to the approach reported by Ramanakumar et al.²⁷ aimed to reduce the effect of inter-, and intra-assay variations. The seropositivity cut points were determined by three standard deviations above the mean values obtained from a pool of HPV-negative samples.

Heatmap analysis

The heatmap analysis of seropositivity for HPV16 and HPV18 versus vaccination and infections status was made with the free online software Heatmapper (http://www.heatmapper.ca).²⁸

Statistical analysis

The statistical analysis was performed with the GraphPad Prism 7.0 statistics package software. The study data were expressed in terms of arithmetic mean, standard deviation, number, and percentage. Differences between prevalence were analyzed by Fisher’s exact test. P-values less than 0.05 were considered statistically significant (*P < .05; **P < .01; ***P < .001).

Results

HPV prevalence

A total of 129 female university students were recruited voluntarily. The analyzed group consisted of students from the following careers: Biology (25%), Optometry (12%), Medicine (18.4%), Psychology (30.2%), and Odontology (14.4%). The mean age of the study population was 21.75 ± 2.86 with a range of 17 to 35 years. The clinical characteristics of the patients are grouped in Table 1 (supplementary material). The overall incidence of HPV infection in cervical samples was 74.4% (96 students) (Figure 1a). Among the HPV-positive women, 70 were positive for a high-risk HPV type, accounting for 72.9% (70/96) of the HPV infection, and 26 students were positive for low-risk HPV types. The different HPV genotypes assessed had different prevalence. The frequencies of infection with low-risk HPV types 6 and 11 were 16.3% (21 students) and 3.9% (5 students), respectively (Figure 1b). Concerning genotype distribution of high-risk types viruses, the higher prevalence was for HPV-18 (14%). HPV-31 was the second most frequent genotype (11%), followed by HPV-16 (9%), HPV-51 (9%), and HPV-59 (5%). Moreover, we also detect the presence of the HPV genotypes 39, 45, 53, 66, 68, and 90 in very low frequencies (data not shown). Twenty-three HPV negative samples were further evaluated by qPCR. As expected, even the higher sensitivity of the fluorescent qPCR method, the 23 samples analyzed were confirmed as HPV-negative (data not shown). We did not find multiple HPV genotypes in any sample.

Table 1. Clinical characteristics of the participants. IUD = Intrauterine device

Sample	HPV genotype	HPV risk	Age	Sexual debut age	Number of lifetime sexual partners	Vaccine doses	Contraceptive method	Medical consultation
M 1	18	High	18	15	2	0	Condom	Gynecological evaluation
M 2	18	High	19	17	8	0	Condom	IUD
M 3	51	High	21	17	10	0	Oral contraceptive	Vaginal infection and condylomatosis
M 4	Negative	Negative	22	15	3	0	Oral contraceptive	Subdermal contraceptive
M 5	31	High	22	14	5	0	Condom	Gynecological evaluation
M 6	Negative	Negative	20	17	4	2	Condom	Gynecological evaluation
M 7	6	Low	22	15	3	0	IUD	Cervical screening
M 8	6	Low	24	16	4	0	Condom	Cervical screening
M 9	16	High	19	16	Unknown	0	Condom	Cervical screening
M 10	18	High	19	17	4	2	None	Vaginal infection
M 11	51	High	25	17	10	0	IUD and Condom	Cervical screening
M 12	53	Low	19	16	6	0	None	Vaginal infection
M 13	6	Low	19	14	3	1	Condom and contraceptive injection	Cervical screening
M 14	6	Low	22	20	12	0	Condom	Cervical screening
M 15	68	High	21	18	2	0	Condom	Gynecological evaluation
M 16	18	High	19	16	Unknown	0	Subdermal contraceptive and Condom	Vaginal infection
M 17	Negative	Negative	24	20	2	1	Condom	Cervical screening
M 18	Negative	Negative	21	14	4	2	Condom	Cervical screening
M 19	18	High	22	18	7	1	Condom	Condylomatosis
M 20	51	High	20	19	Unknown	0	None	Vaginal infection
M 21	31	High	20	17	1	0	Condom	Gynecological evaluation
M 22	Negative	Negative	24	17	2	0	Condom	Vaginal infection
M 23	Negative	Negative	23	16	4	0	Condom	Cervical screening
M 24	16	High	23	18	6	0	Condom	Cervical screening
M 25	31	High	24	19	3	2	Condom	Vaginal infection
M 26	Negative	Negative	19	17	2	1	Condom	IUD
M 27	59	High	19	17	6	2	Condom	Vaginal infection
M 28	16	High	21	16	6	0	IUD and Condom	Cervical screening
M 29	59	High	19	15	1	2	Condom	Gynecological evaluation
M 30	Negative	Negative	19	17	2	0	Condom	Vaginal infection
M 31	51	High	22	20	4	0	None	Cervical screening
M 32	51	High	18	17	6	3	Condom	Vaginal infection
M 33	66	?	19	17	1	1	Subdermal contraceptive	Cervical screening
M 34	51	High	19	16	5	0	None	Estudio de papanicolaou
M 35	Negative	Negative	18	18	1	0	Condom	Gynecological evaluation
M 36	16	High	21	20	2	1	Condom	Vaginal infection
M 37	18	High	19	17	1	3	Condom	Gynecological evaluation
M 38	31	High	24	19	4	0	Condom	Cervical screening
M 39	31	High	20		1	0	Condom	Gynecological evaluation
M 40	6	Low	20	18	1	0	Condom	Gynecological evaluation
M 41	16	High	23	21	1	0	Condom	Gynecological evaluationpor lesión en mama
M 42	31	High	20	16	4	0	Oral contraceptive and Condom	Condylomatosis
M 43	51	High	20	16	5	2	Condom	Condylomatosis
M 44	18	High	20	16	5	0	Condom	Vaginal infection
M 45	6	Low	21	19	4	2	Condom	Gynecological evaluation
M 46	Negative	Negative	22	17	6	0	Condom	Cervical screening
M 47	18	High	19	16	6	1	Subdermal contraceptive	Condylomatosis
M 48	16	High	19	18	2	0	Condom	Vaginal infection
M 49	31	High	24	19	2	0	IUD	Gynecological evaluation
M 50	Negative	Negative	23	16	3	2	None	Cervical screening
M 51	Negative	Negative	23	21	1	0	Condom	Condylomatosis
M 52	6	Low	22	18	3	0	None	IUD
M 53	16	High	21	18	2	0	Condom	Condylomatosis
M 54	Negative	Negative	19	18	1	1	None	Vaginal infection
M 55	31	High	17	16	2	0	Condom	Cervical screening
M 56	31	High	21	20	2	0	Condom	IUD
M 57	Negative	Negative	21	15	3	1	Condom	Cervical screening
M 58	51	High	18	16	2	2	Condom	Cervical screening
M 59	51	High	22	17	4	0	Condom	Condylomatosis
M 60	Negative	Negative	24	20	2	0	None	Cervical screening
M 61	Negative	Negative	35	18	8	0	Condom	Gynecological evaluation
M 62	18	High	25	15	15	0	IUD	Cervical screening
M 63	Negative	Negative	24	20	2	0	Condom	IUD
M 64	59	High	23	15	7	0	Condom	Cervical screening
M 65	31	High	19	14	29	0	IUD	Cervical screening
M 66	18	High	22	17	3	0	Condom	Gynecological evaluation
M 67	11	Low	23	15	5	0	Condom	Gynecological evaluation
M 68	6	Low	22	16	5	0	Condom	Cervical screening
M 69	Negative	Negative	19	18	3	0	Condom	Subdermal contraceptive
M 70	Negative	Negative	20	21	1	0	Condom	Vaginal infection
M 71	59	High	18	17	2	2	Condom	Condylomatosis
M 72	31	High	25	16	10	1	Subdermal contraceptive and Condom	Cervical screening
M 73	39	High	22	17	1	0	Condom	Vaginal infection
M 74	16	High	24	15	3	0	Oral contraceptive	Gynecological evaluation
M 75	Negative	Negative	21	20	1	2	Condom	Pregnancy test
M 76	18	High	25	18	4	0	Condom	Vaginal infection
M 77	45	High	23	22	1	0	Condom	Cervical screening
M 78	18	High	22	17	9	0	IUD	Vaginal infection
M 79	18	High	26	15	5	0	Condom	Gynecological evaluation
M 80	6	Low	20	14	8	0	Condom	IUD
M 81	6	Low	21	16	5	1	Condom	Gynecological evaluation
M 82	11	Low	22	16	10	0	Condom	Cervical screening
M 83	Negative	Negative	26	18	10	0	IUD	Cervical screening
M 84	16	High	19	17	1	1	Condom	Vaginal infection
M 85	16	High	20	15	3	2	Condom	Vaginal infection
M 86	18	High	28	17	1	0	Condon	Gynecological evaluation
M 87	31	High	28	24	5	0	Condom	Vaginal infection
M 88	6	Low	25	17	2	0	Contraceptive injection	Subdermal contraceptive
M 89	6	Low	19	16	1	0	Condom	Cervical screening
M 90	16	High	23	15	7	1	IUD	Cervical screening
M 91	18	High	21	14	5	0	IUD	Cervical screening
M 92	6	Low	18	15	7	0	None	Condylomatosis
M 93	Negative	Negative	23	16	3	0	Condom	IUD
M 94	6	Low	24	19	2	0	Oral contraceptive and Condom	Gynecological evaluation
M 95	Negative	Negative	25	21	2	0	Oral contraceptive	Estudio de papanicolaou
M 96	Negative	Negative	28	21	6	0	Oral contraceptive and Condom	Cervical screening
M 97	11	Low	19	15	2	1	Condom	Gynecological evaluation
M 98	59	High	24	18	2	0	Condon	Cervical screening
M 99	59	High	24	17	1	0	Condom	Cervical screening
M 100	51	High	23	17	4	0	Condom	Condylomatosis
M 101	Negative	Negative	29	19	20	3	Condom	Cervical screening
M 102	51	High	20	18	2	2	Condom	Cervical screening
M 103	11	Low	18	15	1	0	Condom	Condylomatosis
M 104	Negative	Negative	21	20	3	0	Condom	Gynecological evaluation
M 105	Negative	Negative	22	16	3	0	IUD	Gynecological evaluation
M 106	Negative	Negative	27	16	1	0	Oral contraceptive	Cervical screening
M 107	11	Low	17	16	1	2	Condom	Gynecological evaluation
M 108	6	Low	27	17	4	0	Oral contraceptive	Gynecological evaluation
M 109	18	High	24	22	3	1	Condom	Cervical screening
M 110	6	Low	20	15	10	2	Condom	Condylomatosis
M 111	31	High	20	17	2	1	Oral contraceptive and Condom	Subdermal contraceptive
M 112	6	Low	23	21	6	0	Condom	Gynecological evaluation
M 113	59	High	25	17	5	0	Condom	Gynecological evaluation
M 114	Negative	Negative					None	Gynecological evaluation
M 115	45	High	20	16	3	2	IUD	Gynecological evaluation
M 116	90	?	26	17	5	0	None	Vaginal infection
M 117	16	High	25	15	11	1	Subdermal contraceptive	Cervical screening
M 118	6	Low	25	17	2	0	Condom	Cervical screening
M 119	6	Low	24	17	5	0	Subdermal contraceptive and Condom	Cervical screening
M 120	51	High	24	16	1	0	Condom	Vaginal infection
M 121	18	High	20	19	1	2	Condom	Subdermal contraceptive
M 123	Negative	Negative	21	17	3	0	Oral contraceptive	Cervical screening
M125	Negative	Negative	20	18	2	0	Condom	Gynecological evaluation
M126	Negative		20	14	4	0	Condom	Cervical screening
M127	Negative	Negative	21	18	2	0	Condom	Cervical screening
M128	31	High	19	16	2	0	Condom	Gynecological evaluation
M129	18	High	19	17	3	3	Subdermal contraceptive	Gynecological evaluation
M124	6	ND	19	16	2	0	Condom	Cervical screening
M122	6					0	None	Gynecological evaluation

Figure 1. The overall prevalence of HPV genotypes in female university students. a) The infection prevalence is expressed as the number of individual positives for each genotype. From 129 samples, 74.41% (96 students) were positive for HPV (gray bar). b) Comparison of the most common HPV genotypes detected in all samples. The prevalence of most common high-risk HPV was HPV-18 13.95%; HPV-31 10.85%; HPV-16 9.3%, HPV-51 9.3%; and HPV-59 5.4%. c) Age distribution of sexual debut among students. The median age of coitarche for HPV- was 17.5 and for HPV+ girls was 17 years, which were used as a cutoff value for stratification of the girls by age. Fisher’s exact test showed that girls over 17 years old have a higher risk of HPV infection (P = .016). d) Comparison of the distribution of sexual partners among HPV- and HPV+ girls showed a statistically significant association between HPV infection and an increased number of lifetime sexual partners (Fisher’s exact test, P = .038). The cutoff value was the median number of lifetime sexual partners that was 3 for HPV- girls and 3.5 for HPV+ girls.

The average age of sexual debut among students was 17.1 ± 1.9 years old. Among the HPV negative and positive cases, the average age at coitarche was 17.8 ± 2, and 16.9 ± 1.9, respectively. To analyze whether there was an association between the age of sexual debut and genital HPV infection, we stratified the population of girls into two age groups, the cutoff value for age distribution was the median age that was 17.5 for the HPV- group and 17 for HPV+ girls. We found a statistically significant association between early age at coitarche and the presence of HPV infection, the frequency of infection at early age (14–17 years) of sexual debut was 82% in comparison with 60% when girls were above 17 years old (Figure 1c, Fisher’s exact test P = .016). Regarding the number of lifetime sexual partners, the mean was 4.16 ± 3.85 for the whole group. For HPV-negative girls, the mean number of sexual partners was 3.6 ± 3.6, and for the HPV+ group was 4.4 ± 3.9. There was a significant association between HPV infection and the number of lifetime sexual partners (Figure 1d, Fisher’s exact test P = .038). The cutoff value for the distribution was the median number of sexual partners that was 3 for HPV- girls and 3.5 for HPV+ girls.

Seroprevalence and HPV infection

Of the 129 individuals enrolled, 38 reported having received at least one vaccination against papillomavirus, while 86 were not vaccinated, and 5 have no record. The girls without vaccination records were excluded from the analysis. The percentage of vaccination was 30.6. Among vaccinated women 17 had received one dose; the same number had received two doses, and 4 had three doses. Sexual behavior analysis indicates that the median age of sexual debut and the median number of lifetime sexual partners were not significantly different between vaccinated and unvaccinated girls (T-test P = .245 and P = .822, respectively). In addition, vaccine effectiveness (VE) was analyzed against incidence of infections as follows: VE = (Risk among unvaccinated group – Risk among vaccinated group)/Risk among unvaccinated group. In our study, the vaccine effectiveness against incident HPV 16/18 infection was 58.1%; nevertheless, a cross-protective HPV infection was not observed since the vaccine effectiveness was 4.9% for all HPV genotypes infection.

Next, we compare the prevalence of the main HPV types found in our cohort in unvaccinated girls, and those with one, two or three doses of HPV vaccine. We found a downward trend in the prevalence of infection as the number of doses of the vaccine increased for HPV6, HPV16 and HPV31 (Figure 2a). In contrast, for HPV18 there were not differences observed for one and two doses, but there was an increase in the prevalence in girls with three doses of the vaccine. For HPV51 there was a decrease of the prevalence in girls with one dose of the vaccine, but not differences were observed in the group with two doses, and similarly to HPV18 in girls with three doses there was an increase in the prevalence.

Figure 2. HPV-vaccination decrease the prevalence of some HPV types. a) The effect of HPV vaccination was evaluated in girls who received one, two or three doses of the vaccine in comparison with non-vaccinated girls. The percentages of infections were for HPV6 16% (12/75 girls), 12% (2/17), 12% (2/17) and 0% for unvaccinated, one dose, two doses and three doses of the vaccine, respectively. For HPV16 were 9% (7/75),18% (3/17), 6% (1/17), and 0%; for HPV18 were 12% (9/75), 18% (3/17), 12% (2/17), and 50% (2/4); for HPV31 were 15% (11/75), 12% (2/17), 6% (1/17), and 0%; and for HPV51 were 8% (6/75), 0%, 12% (2/17), and 25% (1/4). Data were analyzed by Fisher’s exact test. P-values less than 0.05 were considered statistically significant (*P < .05; **P < .01; ***P < .001). b) Percentage of seropositive individuals for HPV-16-specific IgG antibodies (left graph), and IgG anti-HPV-18 antibodies (right graph) in girls based on their vaccination status. Forty-three percent (32 of 75) of unvaccinated girls were seropositive for HPV16 and 41% (31/75) for HPV18. Forty-four percent (7/16) of vaccinated girls were seropositive for HPV16 after one dose of the HPV vaccine, 56% (9/16) for two doses, and 75% (3/4) for three doses. For HPV18 38% (6/16) of vaccinated girls were seropositive for one dose of the vaccine, while 56% (9/16) and 75% (3/4) were seropositive after two and three doses, respectively. c) Prevalence (percentage) of specific HPV infection in girls that were seropositive for HPV16 (left graph) or HPV18 antibodies (right graph). White bars represent unvaccinated girls and gray bar vaccinated girls. Data were analyzed by Fisher’s exact test. P-values less than 0.05 were considered statistically significant (*P < .05; **P < .01; ***P < .001).

The increased rate of infection for HPV18 after vaccination was unexpected, but it could be related to a decay in the levels of specific antibodies along the time according to some reports. For that reason, we decided to evaluate the presence of specific serum IgG antibodies against HPV16 and HPV18. The seropositivity observed between unvaccinated and vaccinated girls showed an evident impact of HPV vaccination from two doses of the vaccine, and it increased with the third dose (Figure 2b). The prevalence of HPV infection was affected by the presence of specific antibodies against the virus; the analysis with specific HPV types showed similar frequencies of infection in seronegative girls: HPV6, HPV16, and HPV18 have a prevalence of 14% (9/63 girls), and for HPV 31 and HPV51 the prevalence was of 13% (8/63). In contrast, for seropositive girls there was a significant reduction for HPV16 with only 2 of 51 girls infected (4%), for HPV 31 and HPV51 there was a slight decrease in the prevalence, but it was not statistically significant (Figure 2c); no effect was observed for other HPV types. For HPV-18 seropositivity, even though there were no significant changes in the prevalence of HPV-18 infection, it was observed a slight decrease, similarly for HPV6, HPV31, and HPV51 in seropositive girls (right graph in Figure 2c).

Of note, we observed in HPV-18 seropositive girls a significant decrease in the prevalence of HPV16 infection (Figure 2c), which could probably be due to the presence of specific antibodies for HPV-16 in the same individuals. For that reason, we perform a heatmap that shows the seroprevalence for both types of antibodies in all the individuals analyzed together with the infection and vaccination status (Figure 3). There were 34 girls that showed seropositivity for both HPV16 and HPV18, while 17 girls were seropositive only for HPV16, and 15 girls seropositive only for HPV18. Interestingly, the double seropositivity can be acquired by both natural infection according to data displayed in the rectangles 1 and 2, or by vaccination as observed in the rectangles 3 and 4 (Figure 3).

Figure 3. Heatmap analysis of seroprevalence, infection, and vaccination. Analysis of the levels of specific antibodies against HPV16 and 18, and their association with the presence of HPV infection and vaccination status were analyzed by heatmap. For the antibody levels the scoring star with blue color for negative or low levels and change to red color for higher amounts of antibodies detected by ELISA. For the vaccination status, the blue color represents unvaccinated girls and the increasing intensity in red represents 1, 2, or 3 doses of the HPV vaccine. Finally, for HPV infection blue color represent HPV- girls and light red color represent girls with HPV infection (any HPV). Each horizontal line represents one girl, and the sample number is indicated in the left side of the graph. Rectangles 1 and 2 show unvaccinated double seropositive girls, while rectangles 3 and 4 show vaccinated double seropositive girls.

Finally, to assess whether the presence of HPV infection could affect the immune response against the virus we compare in the group of unvaccinated girls the seroconversion rate between HPV-, HPV+ or specific HPV for the antibodies measured (Figure 4). We found that seroconversion was decreased for anti-HPV16 antibodies when the girls were positive for HPV16, but not for other genotypes (Fisher exact test, P = .027). Similarly, this effect was also observed for HPV18 seroconversion in the presence of HPV18 infection (P = .002), however for HPV18 the inhibitory effect in the seroconversion was also observed with the group HPV+ that include all genotypes of HPV detected.

Figure 4. HPV seropositivity according to cervical HPV infection status. Percentages of HPV16 or HPV18 seropositive girls, which were grouped in those negative for HPV (HPV-), positive for any HPV (HPV+), and positive only for HPV-16 (HPV-16+) or HPV-18 (HPV-18+). Data were analyzed by Fisher’s exact test. P-values less than 0.05 were considered statistically significant (*P < .05; **P < .01).

Discussion

Human papillomavirus is an infectious agent of epithelial tissue with high clinical relevance for its association with the generation of cervical carcinoma. HPV is considered to be among the most common sexually transmitted infections. The HPV infection is detected frequently in young women; the epidemiological data indicated that at least 70% of sexually active women are infected at least once during their lifetimes.⁸ In this study, a group of women attending an outpatient clinic for gynecological exams were assessed for molecular detection of HPV in cervical samples, showing that the overall HPV infection rate in the students that were analyzed was 74.4%. In comparison to other studies, the frequency reported in our research was higher than the world prevalence of 10.4%, and some regional rates such as Asia (8%),^12,29 Africa (22.1%), Central America, Mexico (20.4%), and North America (11.3%).^30,31 However, the population in our study is considered at high risk for HPV infection according to several studies in which college young women had a high prevalence. This population has been reported to have a prevalence of 30–50%.^23–25 Moreover, it is important to consider the fact that young women attending clinics often do so because of signs or symptoms of a genital tract infection, which increases the possibility of having a sexually transmitted infectious disease. The high incidence of HPV infection reported in this study is closely related to the different risk factors associated with the participants such as the coitarche age between 14 and 17 years of age (65%), multiple sexual partners (4 on average). Domínguez et al. mention that the main risk factors for HPV infection in women under 25 years of age are both the early sexual intercourse debut and having multiple sexual partners.³² Another reason to which we attribute the high incidence of HPV infection was due to the fact that a significant percentage of the patients attended a gynecological evaluation for both presence of condylomatosis or another infection as observed in supplementary Table 1.

Understanding the HPV prevalence and genotype distribution is important for the planning of diagnostic and preventive strategies toward HPV-related diseases. HPV types 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, and 59 are the most common in the general female population worldwide, accounting for 70% of HPV infections in the presence of normal cytological findings.¹⁰ However, some reports have shown that the prevalence and genotype distribution of HPV varies greatly all over the world.¹² In our cohort, the HPV infection was associated mainly with high-risk viruses (72.9%), which increase the possibility of cervical cancer development. Some studies have reported that HPV-58 is the most frequent genotype in the southern regions of Mexico, with a prevalence of 55.9%.¹⁴ Further, Gallegos-Bolaños et al. in 2017 reported that the HPV-16, 33, 51, and 52 were the most prevalent types in Mexico City³³. In our study, we found that HPV-18 showed the highest prevalence (14%) followed by genotypes 31, 51, and 16. The discrepancy in prevalence with the Gallegos-Bolaños study might be due to population differences because our study included women from Mexico City and the State of Mexico. In line with this hypothesis, it has been observed that depending on the geographical region, some viral types are more frequently found than others. For example, HPV types 33 and 31 are more prevalent in Europe and the USA, types 35 and 45 in Africa, and types 52 and 58 are more frequently observed in Asia.^12,34 Moreover, the age range was also different between both studies, which also can account for the different prevalence observed. Regarding sexual behavior, we found an association of HPV prevalence with the early age of sexual debut and with a higher number of lifetime sexual partners, these sexual behaviors have been previously reported as a risk factor for HPV infection in different populations.^35,36 Taking into account the proportion of girls that had an early sexual debut (65.9%), and the proportion of those who have had more than three sexual partners (46.6%), the population analyzed here can be considered as a high-risk population for HPV infection, which is confirmed by the high prevalence of the virus observed.

In addition to sexual behavior, immunization against HPV is critical to prevent dissemination of the viral infection, especially in young populations. The low vaccination rate (30%) reported in this study could be in part responsible for the high prevalence observed for HPV infection. The low HPV vaccination coverage in our cohort reflects the fact that the HPV vaccine was introduced in Mexico in 2008 with very low coverage, and it was included in the national vaccine program until 2012 with only 52 thousand doses applied.^9,20 Between 2013 and 2018 an average of 1.4 million doses were administered every year,²⁰ this amount covered less than one-third of the female population between 9 and 14 years of age, which was assessed at 6.5 million in the year 2020.²¹

Based on the effectiveness of HPV vaccine in our study, we found that the efficacy of vaccination is lower than other reports. Basu et al. showed that a single dose of HPV vaccine provides a similar protection against persistent infection from HPV 16 and 18 to that provided by two or three doses. The authors mention that adjusted vaccine efficacy against incident HPV 16/18 infection was 63.5% with a single dose with very low variation with two or three doses.³⁷ In our study, the effectiveness was 58.96% against incident HPV 16/18. We could attribute these differences to some sociodemographic characteristics of the participants, for example, in our study, the coitarche age was before 20 years of age and most patients have had multiple sexual partners, which increases the risk to HPV exposure compared to the low-risk exposure of women studied by Basu et al. The sociodemographic differences between the participants could explain the differences in vaccine effectiveness in both populations. It would be interesting to increase the number of participants to assess the vaccine effectiveness, in both, incident and persistent HPV 16/18 infection.

The impact of HPV vaccination in our cohort was evident from the second dose of the vaccine, in contrast with other studies that reported that one dose can induce seroconversion in most individuals, especially with HPV16.^38,39 A possible explanation is a decay in the levels of the specific antibodies to HPV16 and HPV18 along the time as it has been described in other reports, where the seropositivity decay after 18 and 36 months.^40,41 This possibility may be the main reason in our cohort because the average age of the girls was 22 years, meaning that they have around 10 years post vaccination. According to the study of Sankaranarayanan et al. in 2016, the levels of antibodies against HPV in girls with one dose of the HPV-vaccine can decrease below the limit of seropositivity after 36 months of vaccination.⁴¹

Another factor that affects the seropositivity is the fact that in previous studies they have excluded girls with positive infection at the time of recruitment, this condition in our cohort is unknown and it could affect the development of the immune response according to some reports that showed that HPV infection delays the seroconversion and even prevent it.^42,43 Indeed, in one report the efficacy of the vaccine among individuals with cervical HPV16 or HPV18 DNA detected at enrollment was only 25% compared with 85% for women who were negative for HPV.⁴⁴ To assess this possibility, we compare the seroconversion rate between HPV-, HPV+, HPV16+, or HPV18+ girls (Figure 4). We only performed this analysis with unvaccinated girls to avoid the seroconversion due to vaccination. In agreement with the previous reports mentioned early, we found a decreased rate in the seroconversion when the girls were infected by the same HPV genotype of the antibodies analyzed. This delay or decrease in seroconversion could be due to some immune-regulatory mechanisms associated with the virus such as the induction of anti-inflammatory cytokines like IL-10 in keratinocytes, macrophages, and Langerhans cells.⁴⁵ Moreover, it has also been reported an increased frequency of regulatory T cells in cervical lesions that had persistent HPV16 infection.⁴⁶

Of note, the prevalence of HPV-16 infection was significantly decreased in girls that were HPV-16 seropositive in comparison with those that were negative for HPV-16-specific antibodies, but the prevalence of other HPV types was unaffected. This clearly shows that the induction of specific antibodies to the L1 protein of HPV-16 is protective against the infection with HPV16, but they are not protective against other types of HPV. This is in line with other clinical studies.^47–50 In contrast, the prevalence of HPV-18 infection in seronegative and seropositive girls for HPV-18 specific antibodies was similar between both groups. This might reflect a less effective immunization for HPV18 antigens or shorter protective levels of neutralizing antibodies, in support of this presumption in other reports the levels of antibodies against HPV18 are lower than those observed for HPV16 after vaccination, and they have a faster decay effect along the time.^39–41 Interestingly, we observed that HPV18 seropositive girls had a reduction in the number of HPV16 infections, suggesting the presence of seropositivity for this virus that could be associated with vaccination, which includes both viruses. In a heatmap analysis, we corroborate that 34 of 49 girls that were positive for HPV18-antibodies were also positive for antibodies to HPV16, nevertheless, this double seroconversion is not only associated with vaccination, but also with natural infection. This suggests that unvaccinated girls, that were seropositive for both HPV16 and HPV18 viruses were infected by both viruses or there is a crossed immune response, this last possibility could also explain why HPV18 seropositive girls did not have a significant reduction in the rate of infection because crossed immune response is less effective than direct immune response against the specific HPV.^51–53

This study has some limitations that should be considered. Though there was statistically significant analysis between variables, the number of patients included in the study was low. Only some risk factors related to sexual behavior were recorded during enrollment like the number of lifetime sexual partners and age of sexual debut, but others like the type of sexual activity, use of condoms, number of current sexual partners are missing.

In conclusion, this study shows that the incidence of high-risk HPV infection was higher than that reported previously in some regions of Mexico, being the most common genotypes HPV-18, HPV-31, and HPV-51 for female college students. The prevalence and genotype distribution of HPV provides the basis for designing HPV prevention programs including vaccination. The observation that HPV-specific antibodies are mainly protective for the virus used in immunization is very important for the application of HPV vaccines in each geographic region where the viral types can change, and it suggests that the new nonavalent vaccine would be more beneficial in our population, that have a prevalence of diverse HPV types. Furthermore, the coverage of early vaccination in girls who are naïve to HPV infection is of vital importance to reduce the prevalence of HPV. Our study was focused on a high-risk population for HPV infection; therefore, the data obtained can be used for the design of strategies aimed to target this population.

Consent to participate

Written informed consent was obtained from each patient before blood and tissue donation.

Consent for publication

This work has not been published before; it is not under consideration for publication anywhere else, and its publication has been approved by all co-authors.

Ethics approval

All procedures performed in studies involving human participants were following the ethical standards of the institutional research committee of the FES Iztacala and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

Acknowledgments

We thank Estanislao Antonio Calixto for the technical assistant.

Data availability statement

The data sets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Supplementary material

Supplemental data for this article can be accessed on the publisher’s website at https://doi.org/10.1080/21645515.2022.2028514

Additional information

Funding

This study was supported by the program: “Programa de Apoyo a los Profesores de Carrera (PAPCA)” from the Faculty of Higher Studies Iztacala (FES Iztacala) of the National Autonomous University of Mexico (UNAM). Project number FESI-DIP-PAPCA-2016-10.