Objective National Spanish studies show that prevalence of cervical human papillomavirus (HPV) infection in the female population is increasingly frequent, with an overall estimate of 14% in women aged 18–65 years. The objective of this study is to know the prevalence and distribution of HPV types in the female population of the Canary Islands prior to the introduction of HPV vaccines and to investigate the associated clinical and sociodemographic factors.
Methods Based on the Primary Health Care database, a sample of adult women (aged 18–65 years) of Gran Canaria (GC) and Tenerife (TF) stratified into nine age groups was carried out between 2002 and 2007. Women were contacted by postal letter and telephone call and were visited in their primary care centre. A clinical-epidemiological survey was completed and cervical samples were taken for cytological study and HPV detection. HPV prevalence and its 95% CI were estimated, and multivariate analyses were performed using logistic regression to identify factors associated with the infection.
Results 6010 women participated in the study, 3847 from GC and 2163 from TF. The overall prevalence of HPV infection was 13.6% (CI 12.8%–14.5%) and 11.1% (CI 10.3%–11.9%) for high-risk types. The most frequent HPV type was 16 followed by types 51, 53, 31, 42 and 59. HPV types included in the nonavalent vaccine were detected in 54.1% of infected women. Factors associated with an increased risk of infection were: young ages (18–29 years), the number of sexual partners throughout life, not being married, being a smoker, and having had previous cervical lesions or genital warts.
Conclusions It is confirmed that prevalence of HPV infection in the female population of the Canary Islands is high, but similar to that of Spain, HPV 16 being the most frequent genotype. The determinants of infection are consistent with those of other populations.
- gynaecological oncology
- infectious diseases
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
Strengths and limitations of this study
This is the first prevalence study of human papillomavirus infection in Canary Islands.
The study design is population-based, including the main healthcare centres of the participant regions.
Cytological and molecular samples were analysed in the same laboratory by the same staff, using highly sensitive and partially automated techniques that ensured consistency, homogeneity and reproducibility of diagnostic methods.
Study recruitment time was extensive, from 3 to 6 years depending on the region.
Characteristics of the study participants could be different over time.
Cervical cancer is the fourth most common female cancer worldwide and the second most frequent among young women aged 15–44 years, with an estimated 569 847 new cases in 2018.1 In Spain, cervical cancer is the 15th most frequent cancer in women (4th in women aged 15–44 years), with an estimated 1942 new cases in 2018.1 In the Canary Islands autonomous community, 356 new cases were diagnosed in 2008–2011, with a crude rate of 10.1 cases per 100 000 women,2 one of the highest incidence rates in Spain.3
Human papillomavirus (HPV) is a necessary but not sufficient cause of cervical cancer.4 More than 200 HPV genotypes are currently known, epidemiologically classified into low-oncogenic risk (LR-HPV) and high-oncogenic risk (HR-HPV) types.5 HR-HPV types include 16 and 18 genotypes, present in >70% of cervical cancer cases6 and included in the three prophylactic HPV vaccines currently commercialised.7 8
No robust estimations of HPV infection prevalence are available for the Canary Islands, which hinders comparisons with the rest of Spain. Changes in Spanish women’s sexual behaviour in the last decades have leads to increased HPV infection rates (up to 14% in women aged 18–65 years, 29% of them in women younger than 25 years).9 Baseline prevalence estimations of HPV infection and the genotype distributions are essential to monitor the impact of HPV-vaccination campaigns. Therefore, the goal of this study was to estimate the prevalence and distribution of HPV types in the female population of the Canary Islands before introducing HPV vaccination, as well as to study the clinical and sociodemographic factors associated with HPV infection.
The study was conducted between 2002 and 2007 on a sample of women aged 18–65 years living in any of the two most populated Canarian Islands: Gran Canaria and Tenerife. Participants were randomly selected from the regional Health Administration databases, stratified and selected with a probability proportional to the different healthcare areas on both islands. Selected women were stratified into nine age groups (18–24, 25–29, 30–34, 35–39, 40–44, 45–49, 50–54, 55–59 and 60–65 years). The initial sample included 2276 women. For each age group, four reserve groups were obtained to supply the absences or refusals to participate. Participants were contacted by letter and a subsequent telephone call. A visit to the nearest healthcare centre was scheduled. A total of 24 345 letters were sent, 15 577 in Gran Canaria and 8768 in Tenerife, of which 23.3% agreed to participate. Women who did not attend the visit were recalled by phone to schedule another visit. Subsequently, a group of 934 women from Gran Canaria asked to participate in the study (volunteers) of which 665 finally attended the arranged appointment.
Patient and public involvement
No patients or the public were involved in the design, conduct, reporting or dissemination of this study.
Participants were asked to fulfil an informed consent form and to complete a clinical and epidemiological questionnaire (adapted from International Agency for Research on Cancer (IARC) surveys). For cytological collection, the wooden Ayre spatula and endocervical brush (cytobrush) were used, stained with the Papanicolaou technique and the cytological diagnosis was made by a single pathologist according to the criteria of the Bethesda system. For the molecular study, a sterile cotton-tipped polystyrene swab without culture medium (Deltalab, Spain) was used. The obtained cell pellet was subjected to enzymatic digestion with stirring for 2 hours at 55°C with proteinase K following the inactivation of the process with incubation for 10 min at 90°C and subsequent centrifugation, obtaining DNA from the sample supernatant.
To detect HPV infection, two separated PCR were conducted: one using My09/My11 consensus primer and the other using Gp5+/Gp6+ consensus primer. DNA quality was evaluated by PCR testing for the β-globin gene. Samples that were negative for both HPV DNA and β-globin were excluded from the final analysis. Samples showing positive results for any of the HPV PCR reactions or any cytological alteration (atypical squamous cells of undetermined significance (ASCUS) or higher) were genotyped using the Linear Array HPV Genotyping Test (CE-IVD; Roche Diagnostics) or the INNO-LIPA HPV Genotyping Extra Amp kit (Immnogenetics (now Fujirebio Europe), Belgium).
Descriptive analysis of sociodemographic variables was conducted, globally and stratified according to the study subpopulation (ie, selected participants from Gran Canaria, volunteers from Gran Canaria, selected participants from Tenerife). Estimated HPV infection prevalence and genotype distribution and corresponding 95% CIs (Confidence Interval) were calculated as the number of HPV-positive women among the total number of women tested for each age group, study subpopulation and cytological outcome (normal, abnormal). For each genotype, estimated prevalences were calculated independently including the presence of a given type either as a single type or in combination with others (multiple infections). Multivariate analysis was conducted using basic and adjusted logistic regression models in order to assess potential risk factors associated with infections by any HPV type and by HR types. Variables were introduced one by one into a basic regression model adjusted for age group and subpopulation. Variables showing statistically significant association (p value <0.05) were kept as adjustment variables in the final model. Statistical analysis was carried out with the R software (R Development Core Team, 2005, http://www.r-project.org).
Table 1 shows the characteristics of the study population. A total of 6091 women were included: 3212 selected from the general Gran Canaria population (52.7%), 665 volunteers from Gran Canaria (10.9%) and 2214 selected from Tenerife (36.3%). Up to 8.4% of participants were not born in Spain and came mostly from Latin American countries (5.4%); participants’ mean age was 40.7 years; 64.4% were married at recruitment; 77.5% had been pregnant at least once and the mean number of children was 2.2. Regarding cytology screening, 53.7% of subjects had undergone more than five cytological studies in their lives, while 3.5% of them had never undergone one. Regarding HPV infection-related epidemiological factors, 56.5% of subjects were non-smokers and 28.5% were smokers at recruitment; 53.9% of subjects had only one sexual partner and 47.3% were younger than 19 years at first sexual intercourse. Demographic characteristics were slightly different between both islands: education level, proportion of smokers and number of sexual partners were statistically higher in Tenerife than in Gran Canaria. Regarding the subgroup of Gran Canaria volunteers, they were younger, with a high level of education, more divorced or separated, ex-smokers and with more previous cervical pap smears compared with the general population of the island (see online supplementary table 1).
Prevalence of cervical HPV infection
For the prevalence study, 6010 women were included in the analysis after excluding 81 women due to poor DNA quality in their samples. Prevalence of any-type HPV infection was 13.6% (95% CI 12.8 to 14.5) while the prevalence of HR-HPV infection was 11.1% (95% CI 10.3 to 11.9) (table 2). The youngest age group (18–24 years) showed the highest prevalence with 26.7% of any-type HPV infection (95% CI 23.1 to 30.4). Prevalence progressively decreased with increasing age, although the two oldest groups (55–65 years) showed a slightly non-significant increase compared with the immediately younger group (figure 1).
Although volunteers from Gran Canaria showed higher prevalence of any-type HPV infection than participants from the general population from both Gran Canaria (14.5%, 95% CI 11.8 to 17.2 vs 12.7%, 95% CI 11.6 to 13.9; see online supplementary table 2), the difference was not statistically significant. A comparison between the two populations from Gran Canaria (general population and volunteers) and the population from Tenerife showed statistically significant differences in HR-HPV infection prevalence (10.6%, 95% CI 9.6 to 11.6 vs 12.1%, 95% CI 10.7 to 13.4, p=0.002; see online supplementary table 2).
Table 3 shows the distribution of the most frequent HPV genotypes. Single-type HPV infection was detected in 6% of subjects and multiple infections in 7.2% (corresponding to 43.8% and 52.8% of all HPV-positive women, respectively). Among HR-HPV types, type 16 was the most frequent one found in 27.8% of positive women (including both single and multiple HPV types), followed by types 51 (13.7%), 53 (13.3%), 59 (9.9%), 31 (8.5%), 52 (7.7%) and 18 (6.1%).
Among LR-HPV types, type 42 was the most common one (9.3%). In an analysis combining the genotypes included in the HPV vaccines, 31.8% of HPV-positive women were infected by types 16 and/or 18 while the percentage increased to 36.2% when types 6 and/or 11 were added and to 54.1% when the nine types included in the nonavalent vaccine were considered. Figure 1 and online supplementary table 3 show the genotype distribution per age group.
Cytopathological study and cervical HPV infection
The cytological study yielded 317 pathological findings (5.3%) with 69.1% (95% CI 64.0 to 74.2) of HPV positivity vs 5693 non-pathological cytologies (94.7%) with 10.6% (95% CI 9.8 to 11.4) of HPV positivity (see online supplementary table 4), 214 cases of ASCUS were detected (3.6%) with 60.7% of HPV positivity, 91 cases of low-grade squamous intraepithelial lesions (LSIL) (1.5%) with 86.8% of HPV positivity and 12 cases of high-grade squamous intraepithelial lesions or worse (HSIL+) (0.2%) with 83.3% of HPV positivity. Genotype 16 was the most frequently type found in these cytological alterations. Multiple infections were more frequent in women with LSIL or HSIL+ as compared with ASCUS (see online supplementary table 5).
Cervical HPV infection and associated risk factors
Considering all cases of cervical HPV infection (LR-HPV and HR-HPV) and according to the final adjusted model, the following statistically significant variables were detected in the association with HPV infection: younger ages (18–29 years, with a significant lineal trend), not married, smokers, more than one sexual partner (statistically significant trend), history of cervical alterations or genital warts and practising coitus interruptus (table 4). When only cases of HR-HPV cervical infection were considered, the same variables showed statistical significance except for practising coitus interruptus (see online supplementary table 6).
Prevalence of cervical HPV infection
The prevalence of cervical HPV infection (LR-HPV and HR-HPV) in the whole studied population was 13.6% and 11.1% for HR-HPV. HPV prevalence in Spain reported in other published studies ranges from 2.7% to 17.5%.9–15 Two published studies were population-based: one by de Sanjosé et al10 with a random sample of 973 women from the metropolitan area of Barcelona reporting an HPV prevalence of 3.4% (95% CI 2.3 to 4.5), which is rather lower than ours, and one by García et al15 conducted in Castilla y León and reporting 9.6% of HPV prevalence, closer to ours. Differences between both studies could be explained by changes in sexual behaviour in the Spanish population in recent years, with lower age at first sexual intercourse and more sexual partners.16
Non population-based studies include CLEOPATRE (Prevalence and Genotype Distribution of Human Papillomavirus Infection of the Cervix in Spain),9 a study conducted in 17 Autonomous Communities in Spain, using the HC2 test and reporting 14.3% (95% CI 13.1 to 15.5) of HPV prevalence and 12.2% (95% CI 11.1 to 13.4) of HR-HPV infection, both results were similar to ours.
Studies conducted in other European countries reported varied results, with diverse populations and different HPV testing methods. In a review of 18 European studies conducted in 14 countries using the HPV-test as first screening (HC2 or PCR) the HR-HPV prevalence, standardised by age, ranged from 1.7% in Spain to 12.5% in Belgium.17 Bruni et al in a meta-analysis including 1 million women worldwide with normal cytological findings observed 8.8% global adjusted HPV prevalence in Southern Europe, 9% in Western Europe and 10% in Northern Europe.18 Studies conducted among women from different European screening programmes showed HPV prevalences ranging from 6.4% in Germany,19 8.8% in Italy,20 13.7% in France,21 15.2% in Belgium,22 19.4% in Portugal23 to 26.4% in a population-based study in Denmark.24
Prevalence of cervical HPV infection per age group
As expected, the highest HPV prevalence found in our study was observed in women aged 18–24 years (26.7%), an age group potentially associated with a higher number of sexual partners. This finding was also observed in previous Spanish and European studies.9 16 17 In our study, after this first peak in women <25 years, the prevalence declines in older ages, although a slightly, not significant, increased was observed in women older than 55 years. This second peak in older women was also reported by other authors.17 18 20–22 Such a bimodal pattern could be due to changes in the sexual behaviour or the reactivation of latent viral infections,25 HPV types and their variants in such infections, individual susceptibility or regional differences in the screening programmes.18
HPV 16 was the most prevalent genotype in our population, present in 27.8% of positive samples. This prevalence was similar to that reported in other studies in Spain,10 14 although higher than the 16.9% found in the CLEOPATRE study.9 After HPV 16, the most frequent types in decreasing order were: HPV 51, 53, 59, 31 and 52. Our results are similar to most studies conducted in Spain9–11 14 and other European countries.19 21–24
Many studies have reported the percentage of multiple infections9 12 13 15 18–21 23 24 ranging from 8.1% in Spain13 to 54.3% in Denmark.24 The one from Denmark was similar to ours (52.8%), although it included a higher percentage of infections by more HPV types. This finding could be explained by the use of a HPV detection technique (hybridisation technology) with a high sensitivity for detecting multiple infections.
A total of 31.8% of HPV-positive women (4.3% of the total population) were infected by types 16 and/or 18, which were included in the bivalent vaccine. Regarding HPV types included in the quadrivalent vaccine (HPV 6, 11, 16 and 18), at least one of them was found in 36.2% of women (4.9% of the total population). This prevalence increased up to 54.1% with the addition of HPV types 31/33/45/52/58, included in the nonavalent vaccine. Such proportions were higher than those reported in Denmark24 (27.7%) and in the CLEOPATRE study (22.1% in Spain9 and 32.6% in Portugal23). These data illustrate the degree of protection offered by HPV vaccines; one out of three HPV-infected women would have been protected by the bivalent or the quadrivalent vaccine and one out of two women would have been protected by the nonavalent one. However, the frequency of HPV types 51, 53, 59, frequently found in our population, indicate the need to continue the cytological screening population.
Cytopathological study and cervical HPV infection
Cytological alterations found in our study (5.3%) were similar to those observed in other studies, both in Spain9 10 14 and Europe,19 20 22–24 ranging between 1.6% and 7%. The HPV prevalence increased with lesion severity (60.7% in women with ASCUS; 86.8% in women with LSIL and 83.3% in women with HSIL+). This finding was in agreement with other published studies.9 10 12 19 21–24 The HPV prevalence in normal cytologies was 10.6%, similar to that reported by Bruni et al18 in our geographical area (8.8%), although lower than that reported in most studies.9 21–24
Risk factors and cervical HPV infection
Age consistently appears as a risk factor for HPV infection, both in our study and other published ones,14 20 26 27 directly associated with younger women’s sexual behaviour as compared with older ones.
Number of sexual partners in life extensively appears10 11 14 26–28 as a risk factor for HPV infection and was the factor with the largest impact in our study. As in our study, most authors failed to find a relationship with age at first intercourse.10 26 27 This later parameter seems to influence number of sexual partners but does not seem to be an independent risk factor for HPV infection.
In our analysis, not being married (divorced, widow or single) was a statistically significant risk factor for HPV infection, as was also reported in other studies.10 20 26 This finding could be explained by the sexual behaviour of not married women, who may probably have more sexual partners.
Coitus interruptus was the only contraception-related practice found to be associated with higher risk of any-type HPV infection, both in the basic and the adjusted models, although such an association was not found for HR-HPV types. This factor might possibly be linked to younger groups, where other risk increasing factors coexist.
Smoking was a risk factor for HPV infection in our population, in accordance with data reported by other authors,26 27 29 although not by others.10 Quitting smoking has been considered to potentially revert infection risk.29 In order to explain for the relationship between smoking and increased risk of HPV infection, it has been postulated that tobacco and its metabolites may alter the immune system of the cervical epithelium, thus reducing the number of CD4 lymphocytes and Langerhans cells29 and impairing the activity of natural killer cells.
The presence of genital warts and previous cervical alterations was associated with higher risk in our population, as well as in other studies,26 which is not surprising since both events are directly related.
Country of origin especially African ones, appeared as a risk factor for HPV infection in our basic model, although not in our adjusted model. Earlier published Spanish studies showed higher HPV infection risk in women born out of Spain,10 11 26 probably due to differences in the sexual behaviour of men and women.
Regarding parity and HPV infection risk, similar to other authors,20 we found some protective effect in women with one or two births in our basic model for any-type HPV, although not for the adjusted model or for HR-HPV types, a finding also reported by some authors.10 26 27 In a meta-analysis published by the IARC,30 a slight risk increase in nulliparous women (younger and more sexually active) as compared with women who have been pregnant was described.
The relationship between taking oral contraceptives (OC) and the risk of HPV infection is controversial. In our population, a slightly increased risk was found for women taking OC in the basic model though not in the adjusted model, a finding also described in other studies.10 20 26 27 30
Infection by other sexually transmitted diseases analysed in our population increased the risk in the basic model but not in the adjusted model (data not shown), consistent with other published studies.26 27
Some authors have reported no association between using condoms and increased risk of HPV infection14 20 27 28; some even reported some protective effect.26 In our study, like with other contraceptive methods we failed to find an association with HPV infection (data not shown). The evidence is controversial regarding the association between HPV infection and level of education.26 27 31
Strengths and weaknesses
One of the main strengths of our study was our population-based design, which covered the main healthcare centres on the islands and recruited potential participants from an official source, ensuring a random sample. Additionally, the fact that all cytological and molecular studies were conducted in the same laboratory, by the same technical and medical staff, using highly sensitive and partially automated analytic systems ensured consistency, homogeneity and reproducibility of diagnostic methods.
The prolonged recruitment time was a weakness of this study. Three years were needed for Tenerife and 6 years for Gran Canaria, although 2 years had been originally planned. Potential variations over time could have influenced the sociodemographic characteristics of the population. Thus, the characteristics of participants recruited at the beginning of the recruitment period could have been different from those of women recruited by the end.
This study provides population-based references for the prevalence of HPV infection in the Canary Islands, which enables future assessment of the impact of HPV vaccination campaigns. The prevalence of HPV infection in the female population of Gran Canaria and Tenerife was high, although similar to that of previous studies conducted in Spain, with genotype HPV 16 being the most frequent one. These results support the potential benefits of HPV vaccines in terms of reducing infection as well as the consequent development of HPV-related lesions, including cancer.
The authors would like to thank all patients for their participation in the study. The authors would like to thank the colleagues and the study staff for their commitment to data collection and to Dr Jorge Luis Doreste from University of Las Palmas de Gran Canaria for statistical analysis for simple size of the study. The authors would like to thank Tenesoya Alamo, Tanausú de la Cruz and Maria Dolores Navarro for technical assistance. The authors would also like to thank Sanofi Pasteur MSD (Merck Sharp & Dome), Roche Diagnostics SL and Fujirebio Ibérica SL for receiving donations for analytical kits.
Collaborators HPV Canary Study Group. Gran Canaria Team: Diana Alemán, Mónica Almeida, Ana María Arencibia, María Isabel Armas, Guillermina Batista, Victoria Bernal, Francisca Bernaldo de Quirós, Sili Bolaños, Dolores Casaña, Luisa Celedón, Isabel Cruz, Elisa Díaz, Inocencia Duarte, Felisa Expósito, Carmelo Felipe, Carlos Galván, María José García, María Isabel García, Vanesa García, Virginia García, Elena Giménez, Teresa Godoy, Catalina Gómez, Lucia González, Luisa Gutiérrez, Mónica Hernández, Delia Herrera, Laura Herrera, Rosario Laseca, Carmen Marrero, Ofelia Marrero, Noa Mateos, Olivia Medina, Josefa Mendoza, Lucia Montesdeoca, Rosa Monzón, Cristina Morales, Mercedes Morales, M Dolores Navarro, María Ángeles Nieto, Noelia Pérez, Yurena Pérez, Antonio Ramos, Antonio Rico, Margarita Roldán, Esther Salamanca, Rosario Sánchez, Raquel Santana, Elvira Santos, Antonia Solanes, Elisabeth Soutto, Dulce Suarez, María Jesús Suárez, María Ángeles Tadeo, Virgen Valdés, Gabriela Valido, Iralla Vega, Maria del Pino Vega. Tenerife Team: María Angeles Afonso, Elisa Baena, María Pilar Baz, José de Armas, Alicia de la Puerta, Josefina García, María Asuncion González, Célida González, María Teresa Hernández, Josefa Limiñana, Carmen Rosa León, Fernando Marín, Emma Manrique, José Roberto Negrín, Rosa Olavarrieta, Verónica Perera, Concepción Sabater, Candelaria Sosvilla.
Contributors MA designed the study, performed HPV diagnostic molecular methods, data analysis, interpretation of data and drafted the manuscript. ER performed statistical analysis of data, designed the figures and drafted the manuscript. MP performed cytopathological diagnosis. MS performed HPV diagnostic molecular methods. MAS designed and supervised a base data and processed the experimental data. AT, BV, LA, RH, HPV Canary Study Group received the patients and took cervical samples. MdCC and ARdP were involved in planning and supervised the management of cervical a molecular samples. AL, JLT, OA, VB, NM, SC, AQ treated patients with cytological and molecular disorders. LB, SdS and ES aided in interpreting the results and worked on the manuscript. All authors read and approved the final manuscript.
Funding This research has been funded by: Fondo de Investigación Sanitaria (Instituto de Salud Carlos III), grant FIS 00/714; Fundación Canaria de Investigación y Salud, grant FUNCIS 00/14 and FUNCIS 02/19 and Fundación Amurga.
Disclaimer The funders had neither any involvement in the study design, collection, analysis, interpretation of the data, writing of the report nor in the decision to submit the paper for publication.
Competing interests None declared.
Patient consent for publication Not required.
Ethics approval This study was favourably evaluated by the Ethics and Clinical Trial Committee of the hospital Complejo Hospitalario Universitario Insular Materno Infantil.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data are available on reasonable request. The database obtained from this study is kept under the supervision of the authors (MA and ER) in an anonymised form. These data will be shared in a raw form by emailing to firstname.lastname@example.org.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.