Next Article in Journal
Correction: Narongkiatikhun et al. Immunogenicity and Safety of Homologous and Heterologous Prime-Boost of CoronaVac® and ChAdOx1 nCoV-19 among Hemodialysis Patients: An Observational Prospective Cohort Study. Vaccines 2023, 11, 175
Next Article in Special Issue
The Impact of Informational Intervention on HPV Vaccination Intention among Heterosexual Men
Previous Article in Journal
Long Term Follow-Up Study of a Randomized, Open-Label, Uncontrolled, Phase I/II Study to Assess the Safety and Immunogenicity of Intramuscular and Intradermal Doses of COVID-19 DNA Vaccine (AG0302-COVID19)
Previous Article in Special Issue
Acceptability of HPV Vaccines: A Qualitative Systematic Review and Meta-Summary
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Vaccines
Volume 11
Issue 10
10.3390/vaccines11101536
vaccines-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Understanding Prognostic Factors for Human Papillomavirus Vaccination: A Rural Community Case–Control Study

by
Lara Colomé-Ceballos
1,2,*,
Josep Lluís Clua-Espuny
3,4,*,
Concepción Ceballos-García
5,
Josep Clua-Queralt
3,
Maria Jesús Pla-Farnós
6 and
José Fernández-Sáez
4,7
1
Sexual and Reproductive Attention, SAP Terres de l’Ebre, Catalonian Health Institute, 43500 Catalonia, Spain
2
Gynaecology Department, Hospital Verge de la Cinta de Tortosa, Catalonian Health Institute, 43500 Catalonia, Spain
3
EAP Tortosa Est. Primary Care, SAP Terres de l’Ebre, Catalonian Health Institute, 43500 Catalonia, Spain
4
Unitat de Suport a la Recerca Terres de l’Ebre, Foundation University Institute for Primary Health Care Research Jordi Gol i Gurina (IDIAPJGol), 43500 Tortosa, Spain
5
Santa Creu de Jesús Hospital, Salut Tortosa, 43590 Catalonia, Spain
6
Gynaecology Department, Hospital Universitari de Bellvitge, Carrer de la Feixa Llarga, s/n, L’Hospitalet de Llobregat, 08907 Barcelona, Spain
7
Faculty of Nursing, Terres de l’Ebre Campus, Rovira i Virgili University, 43500 Tortosa, Spain
*
Authors to whom correspondence should be addressed.
Vaccines 2023, 11(10), 1536; https://doi.org/10.3390/vaccines11101536
Submission received: 28 July 2023 / Revised: 18 September 2023 / Accepted: 25 September 2023 / Published: 28 September 2023
(This article belongs to the Special Issue Vaccination Progress in COVID-19 and HPV)
Browse Figures
Review Reports Versions Notes

Abstract

:
HPV vaccination coverage rates can vary depending on several factors. The main objective of this study is to identify possible independent prognostic factors that have an impact on HPV vaccination in a rural community, specifically related to sexual and reproductive health. A case–control, retrospective, community-based study was carried out on women aged 15 to 40 in the primary health centers of Southern Catalonia’s Terres de l’Ebre region, Spain, from 1 January 2020 to 31 December 2022. A random sample of 520 women with an average age of 29.3 (SD 7.8) years old was included in the study. Independent prognostic factors: age OR 0.680 (95% CI: 0.635–0.729, p < 0.001), immigrant origin OR 0.215 (95% CI: 0.109–0.422, p < 0.001), and HPV PCR OR 7.402 (95% CI: 2.504–21.880, p < 0.001). The variables that showed a barrier effect for HPV vaccination were age (OR 0.680, 95% CI 0.635–0.729, p < 0.001), and immigrant origin (OR 0.215, 95% CI 0.109–0.422, p < 0.001). The variable that showed a facilitating effect for HPV vaccination was HPV PCR (OR 7.402, 95% CI 2.504–21.880, p < 0.001).

1. Introduction

The human papillomavirus (HPV) vaccine is a critical component of the World Health Organization’s (WHO) 2030 Objectives program. This initiative places a strong emphasis on achieving high vaccination rates among women, with a specific target of attaining a 90% vaccination coverage rate among 15-year-old females [1]. The program aims to strengthen immunization efforts in primary healthcare [2]. However, its implementation shows substantial variability worldwide, with a global estimated coverage rate of only about 15% [3,4].
Regions undergoing demographic changes, such as the one in which this study was conducted, especially due to immigration phenomena, rurality, and socioeconomic characteristics [5,6,7,8,9], should not only evaluate the current status of the WHO program’s implementation but also explore ways to make its execution easier within the context of these demographic shifts. Additionally, they should address changes in the HPV guidelines and how to identify and address potential risk subgroups that may arise as a result [10,11,12,13].
On one hand, it is known that vaccination coverage rates vary due to numerous factors [14]: income level [3,4,15], accessibility [15,16,17], educational attainment [18,19], immigration status [20,21], smoking habits [22,23], and self-perception [24]. Moreover, most factors described are particularly evident in rural communities, resistant to cultural changes and exposed to deficits in basic resources such as the supply, monitoring, and recording of vaccinations of public interest and associated with lower vaccination coverage compared to urban regions. Understanding of these specific factors may help to devise targeted strategies for enhancing HPV vaccination coverage in rural locale communities [25,26,27,28].
On other hand, facilitating factors have been described for vaccination such as practicing safe sex and undergoing regular sexual and reproductive health screenings [29]. Likewise, women who have had Pap smears within the past year are more likely to have received the HPV vaccine [30]; and adolescent females who have received the HPV vaccine are more likely to undergo health examinations, including tests for sexually transmitted infections [31]. In addition, women who have received the HPV vaccine are more inclined to practice safer sexual behaviors when compared to their unvaccinated counterparts [32]. The recommended population screening for women aged 30–35 and older prioritizes HPV detection as the preferred test over conventional cytology. However, since this age group coincides with lower rates of HPV vaccination, it may be faced with new barriers to follow-up vaccination. The recent inclusion of boys in systematic vaccination programs [33] is expected to have a pivotal role in altering the landscape of HPV vaccination [34].
Effectively addressing these challenges necessitates a multifaceted approach that acknowledges the wider context of rural communities, takes into account their unique obstacles, and caters to their specific needs. The primary objective of this study is to ascertain potential autonomous prognostic determinants that impact HPV vaccination within a rural community, with a specific focus on aspects related to sexual and reproductive health.

2. Materials and Methods

2.1. Study Design

A retrospective, case–control, community-based study was carried out on women aged 15 to 40 in the primary health centers of Southern Catalonia’s Terres de l’Ebre region, Spain, with data available up to 31 December 2022.

2.2. Study Scope

The Terres de l’Ebre region (Appendix A) comprises four districts, 52 municipalities, and has a total population of 191,791 inhabitants. The population density is 54.5 inhabitants per square kilometer, which is significantly lower than the average population density in Catalonia. This region is also recognized as part of the demographic phenomenon referred to as “emptied Spain” [5,6,7]. The region demonstrates a progressive aging trend due to a negative migratory balance and a low birth rate, resulting in a higher aging index (162.7) compared to Catalonia (127.1) and Spain (118.43) [8,9,35]. Additionally, the average income per inhabitant is 77.4%, which is lower than the overall average for Catalonia (100%) [10]. A previous article [11] has already described the study population and provided information about their demographics, age at vaccination, vaccination rate, socioeconomic status, and other relevant factors.

2.3. Patients/Subjects of Study

The study’s target population encompassed 24,415 women aged 15 to 40 residing in the Terres de l’Ebre region, constituting 94.7% of the total number of women in that age group registered within the territory.
For this case–control study, a randomized sample size of 520 participants between the ages of 15 and 40 years, who are registered in the territory and have an active medical history in any of the healthcare centers within it, was calculated. The study aimed to include 182 cases of vaccinated patients and 338 controls of unvaccinated patients. The calculations were based on the following parameters: alpha risk of 0.05, beta risk 0.15, and a bilateral contrast to detect a minimum odds ratio of 0.5. A rate of 0.3 exposed in the control group is assumed, and a loss of 10% in follow-up has been estimated. The POISSON approach has been used (https://www.imim.es/ofertadeserveis/software-public/granmo/, accessed on 23 January 2023). The STROBE checklist was used for the review of items included.

2.4. Observation Period

All the data pertain to the clinical history of the subjects up to 31 December 2022.

2.5. Inclusion Criteria

The inclusion criteria were as follows:
−
Availability of clinical record at any of the health centers in the territory (HC3).
−
Women who have been residents for at least five years in one of the four regions of the Terres de l’Ebre territory (Baix Ebre, Montsià, Terra Alta, and Ribera d’Ebre), regardless of their place of birth.

2.6. Exclusion Criteria

The exclusion criteria were as follows:
−
Patients without active or open clinical history and/or with insufficient data or inaccessibility to said data.
−
Disease with vital prognosis of less than one year.

2.7. Variables

2.7.1. Dependent Variable: HPV Vaccination Status

The vaccination status was classified as “correctly vaccinated” if the individual received either 2 doses before their 15th birthday or 2 or 3 doses after their 15th birthday, as well as for immunocompromised patients [36,37,38,39]. The patients were classified as non-vaccinated if they had never received a dose of HPV vaccination and incompletely vaccinated if they had received one dose.
The participants in this study were categorized into two groups based on their vaccination status, as per the criteria outlined by the World Health Organization (WHO) and the Centers for Disease Control and Prevention (CDC) [36,37]. Since 2017, all individuals who received vaccines subsidized by the public health system in the study region were administered the nonavalent vaccine (9vHPV), irrespective of the type of access to the vaccine. Patients under 15 years of age receive two doses, while those of 15 years or more receive three doses. Before 2017, the financed vaccine was the quadrivalent (4vHPV), and all patients received three doses until 2014, two doses from 2014 to 2017 if they were less than 15 years of age, and three doses in the rest of the patients. In the case of the patients vaccinated after medical counselling, the health personnel decide what type of vaccine is administered to the patient [40,41]. The number of doses administered, the type of vaccine, and the date of administration were included in the study.

2.7.2. Independent Variables

-
Demographic variables: age, ethnic origin, immigrant origin, and educational level.
Age was divided into intervals in Table 1 as a categorical variable with ordinal properties, as a binary categorical variable in Table 3, and as a continuous variable in the rest of the calculations. Immigrant origin was defined as individuals who were born in a different country but have moved to Spain to live and work. Ethnic origin was defined based on the country of origin as well as race, and educational level is defined as illiterate (inability to read or write), mandatory primary/secondary schooling, or higher education (beyond secondary high school level, such as bachelor’s, master’s, and doctoral degrees). These data were recorded in the patients’ clinical history.
−
Variables related to cervical cancer screening: first registered cervical cytology, result of the first registered cervical cytology, first registered cervical PCR HPV test, and result of the first registered HPV PCR test.
−
Variables related to sexually transmitted diseases: first registered PCR sexually transmitted disease (STD) test in the vagina or cervix, result of the first registered PCR STD test in the vagina or cervix, and history of any diagnosis of an STD or HIV in the medical record. The STDs that are included in the test available at our laboratory are: Chlamydia trachomatis, Neisseria gonorrhoeae, Trichomonas vaginalis, Ureaplasma spp., and Mycoplasma spp.
−
The variables related to sexual and reproductive health included factors such as childbirth, breastfeeding, voluntary termination of pregnancy or abortion, contraception methods, and the use of condoms.
−
Other clinical variables: smoking and immune status.
The variables have been categorized based on the individual’s vaccination status, captured at the time of HPV vaccination, and the results have stratified into different age groups or periods.

2.8. Data Collection and Information Sources

All study participants belong to the Catalan Institute of Health (ICS), which manages the area through 11 Primary Care Centers (EAPs). Data were obtained retrospectively through an anonymized database provided by the Service of Technology and Communication. The data are recorded using the International Classification of Diseases (ICD-10) and sent to the principal investigator in an anonymized format. Data were supervised and analyzed according to the General Data Protection Regulation of Spain/Europe as of 1 February 2017. The study was conducted in accordance with the most relevant standards on data handling, experimental context with patients, ethics, and data protection and privacy, following Directive 95/46/EC (protection of individuals with regard to the processing of personal data and on the free movement of such data). All collected data were compiled and stored in a dedicated repository specifically created for this study. The repository was then handed over to the principal investigator for analysis. Prior to conducting the study, the research protocol underwent thorough ethical evaluation and received approval from the Ethics Committee of the Jordi Gol Primary Care Research Institute, with the registration number 21/064-P.
The datasets utilized for this project include:
  • The Health Plan of the Terres de l’Ebre Healthcare Region 2021–2025 [42]: This strategic document outlines the goals, priorities, and actions that will guide healthcare services in the Terres de l’Ebre region of Catalonia, Spain, from 2021 to 2025.
  • The HC3 Patient Episode Dataset for Catalonia (CatSalut, Department of Health): This dataset contains demographic and clinical information on inpatient and outpatient care in Catalan hospitals. It provides details on patient episodes of care, including hospital admissions, outpatient visits, and diagnostic tests.
  • The shared clinical information database of the 11 primary care teams managed by the Catalan Health Institute: This database includes clinical data, symptoms, tests, diagnoses, comorbidities, prescribed medication, referrals to secondary and tertiary care, and the alive/dead status of 97.7% of the residents in the territory.
  • Pharmacological variables collected from the Integrated System of Electronic Prescription (SIRE): This dataset captures information on prescribed medications.
  • Data provided by the Statistics Institute of Catalonia: This dataset includes information on gross household income per capita relative to the Catalonia average (100%), population density per square kilometer, and aging index compared to Catalonia (100%) [7,8,10].
Data on these factors was automatically collected whenever possible or manually collected otherwise.

2.9. Statistical Analysis

The vaccination coverage rate was calculated by dividing the number of patients with complete vaccination by the total target population. Data are presented using frequencies and percentages for categorical variables and using means and standard deviations for continuous variables, according to region, vaccination status, and age. To detect differences between the two groups, the chi-square test for categorical variables and the Mann–Whitney U test for continuous variables were used.
The study employed binary logistic regression using the backward stepwise method of Wald to identify the independent predictive variables associated with HPV vaccination. The backward stepwise method of Wald was used to eliminate variables that did not significantly contribute to the model while identifying those that were significantly associated with vaccination. The results were presented in the form of odds ratios (ORs) and their 95% confidence intervals. The ORs enabled the determination of the probability of vaccination for each predictive variable while considering the effect of other variables in the model. The variables analyzed included age, smoking, ethnic origin, immigrant origin, offspring, breastfeeding, contraception, cervical cytology, HPV PCR test, and STD PCR test. The association between variables was calculated using the OR, with a statistical significance level of p < 0.05 established, and 95% CI. The association between variables in the logistic regression model was evaluated using the area under the curve (AUC) of the receiver operating characteristic (ROC) curve. By utilizing the kappa concordance methodology, researchers assessed the agreement between the outcomes derived from the multivariate regression analysis and the real vaccination status against HPV. Specificity was employed to determine the proportion of true negative results (non-vaccinated individuals correctly identified as non-vaccinated) among all individuals who were not vaccinated, while sensitivity measured the proportion of true positive results (vaccinated individuals correctly identified as vaccinated) among all individuals who were vaccinated. All analyses were performed using the IBM SPSS Statistics 20 statistical package.

3. Results

A randomized sample of 520 cases was included, with an average age of 29.3 (SD 7.8 years old). HPV vaccination displayed a significant association with various factors, including regular Pap smear screenings (p < 0.001), contraception usage (p = 0.008), parenthood status (p < 0.001), and the practice of breastfeeding (p < 0.001). In the logistic regression analysis, the age (OR 0.680, 95% CI 0.635–0.729, p < 0.001), immigrant origin (OR 0.215, 95% CI 0.109–0.422, p < 0.001), and HPV PCR (OR 7.402, 95% CI 2.504–21.880, p < 0.001) emerged as independent facilitators of HPV vaccination.

3.1. Basal Characteristics

Statistically significant differences between vaccinated vs. unvaccinated groups are shown for average age in Table 1 and race, immigrant origin, and educational level in Table 2. In total, 27.3% (95% CI 23.3–31.2) of unvaccinated women have an immigrant background and differences by ethnic origins (p = 0.044) and immigrant status (p < 0.001) are shown in relation to vaccination status.
A total of 175 (33.6%) individuals were classified as correctly vaccinated. Among the vaccinated patients, the average number of doses administered was 2.54 (SD 0.59). Moreover, the probability of successfully finishing the entire vaccination schedule once initiated exhibited a 25-fold increase (95% CI 11.74–53.21; p < 0.001).
Significant statistical differences (Table 3) were observed between the screening tests performed and vaccination status, as well as the results of the cervical cancer screening tests and the mean ages of the different groups (p < 0.001), except for the mean age in the different results of the HPV PCR test (p = 0.125).
The average age of individuals without a recorded STD PCR test was 29.4 years old (SD 7.8), while the average age of those with at least one test recorded was 28.9 years old (SD 7.5) (p = 0.793). No statistically significant differences were found in the mean age of individuals who had undergone an STD PCR test or not or in the various outcomes of the test (Table 4).
Significant disparities (Table 5) were observed between the vaccinated and unvaccinated groups across several variables: mean age (21.8 (SD 5.2) vs. 33.1 (SD 5.9), p < 0.001); use of contraception (p = 0.008); parenthood (p < 0.001); breastfeeding practice (p < 0.001); and regular Pap smear screenings (p < 0.001).
Related to other variables, no significant differences were found in terms of immune status alterations (p = 0.170) as shown in Appendix B. In relation to tobacco use, 15.0% (95% CI 11.83–18.16) of the participants were registered as smokers with a statistically significant difference (p < 0.001) as shown in Appendix C. Non-vaccinated smoking patients were significantly older than non-vaccinated non-smoking patients (35.5 (SD 4.3) vs. 33.0 (SD 5.7), p = 0.002).

3.2. Prognostic Factors of HPV Vaccination

The initial model included the following variables: age, cervical cytology test, HPV PCR test, STD PCR test, offspring, breastfeeding, contraception, smoking, immune status, ethnic origin, and immigrant origin. The results (Table 6) showed age and immigrant origin as significant barriers to HPV vaccination. Conversely, cytology and HPV PCR are facilitators for vaccination.
Regarding the ROC curve analysis (Figure 1), the variable age exhibited an area under the curve (AUC) of 0.917 (95% CI 0.89–0.94), indicating strong discriminatory power; immigrant origin, 0.575 (95% CI 0.52–0.62); cervical cytology performed, 0.633 (95% CI 0.52–0.62); and HPV PCR performed, 0.498 (95% CI: 0.44–0.55).
The kappa concordance methodology showed good results, with a Kappa index of 0.755, sensibility of 0.81 (95% CI 0.75–0.86), specificity of 0.93 (95% CI 0.91–0.96), positive predictive value of 0.87 (95% CI 0.82–0.92), and negative predictive value of 0.90 (95% CI 0.87–0.93). Figure 2 shows the probability distribution of the forecast from the logistic regression model.
Eventually, three target risk groups were identified related to challenges in accessing HPV vaccination: women over 25 years old, women of immigrant origin, and women without regular cervical cancer screening, and the probability of successfully finishing the entire vaccination schedule, once initiated, exhibited a 25-fold increase (95% CI 11.74–53.21; p < 0.001). These findings underscore the significance of considering age, immigrant status, and accessibility to cervical cancer screening when addressing the issue of HPV vaccination access.

4. Discussion

A randomized sample of 520 cases was included, with an average age of 29.3 (SD 7.8) years. HPV vaccination displayed a significant association with various factors, including regular Pap smear screenings (p < 0.001), contraception usage (p = 0.008), parenthood status (p < 0.001), and the practice of breastfeeding (p < 0.001). In the logistic regression analysis, the age (OR 0.680, 95% CI 0.635–0.729, p < 0.001), immigrant origin (OR 0.215, 95% CI 0.109–0.422, p < 0.001), and HPV PCR (OR 7.402, 95% CI 2.504–21.880, p < 0.001) emerged as independent facilitators of HPV vaccination.
The characteristics of rurality, low population density, demographic aging, and below-average income levels in the country [5,6,7,8,9,10,35] may be negatively linked to the administration of HPV vaccine doses [43,44,45]. However, our results differ from significant associations with Hispanic ethnicity, county poverty, household characteristics, and uninsured rates, which could be attributed to variations in geographic and time frames [46,47,48].
The results regarding age can be explained by the school calendar factor. In Spain, the introduction of HPV vaccines into vaccination schedules began in 2008. Consequently, only individuals born after 1996 were integrated into school vaccination programs, enabling them to receive publicly funded HPV vaccination while attending school [49]. This creates a distinction from patients born before 1995. If these patients choose to get vaccinated, they require medical advice and must purchase the vaccine using their own funds. Altogether, these factors have become barriers to accessibility for HPV vaccination for immigrant women [17,18,19,20,21]. The study area has a higher percentage of immigrant population compared to the averages documented for Catalonia [50]. Lack of knowledge about the vaccination schedule also poses a barrier for them to access the catch-up option for funded vaccination. These results are similar to those of other studies [17,18,19] that illustrate noteworthy disparities in the initiation of HPV vaccination, influenced by factors such as ethnic background [51,52,53,54,55], parental attitudes [56], limited awareness [57,58], and financial constraints [59].
The findings related to cervical cytology or Pap smear closely paralleled established evidence [11] while also shedding light on certain distinctions, such as a reduced screening coverage and an elevated proportion of abnormal results in contrast to Spain’s national data [60]. Comparable patterns were evident in HPV PCR testing, with the younger vaccinated cohorts displaying elevated rates of positive infections. Moreover, given that PCR testing was identified as a contributing factor to vaccination, ironically, healthcare professionals should view the infection as an opportunity to increase vaccine coverage rates (VCRs) and emphasize the accessibility of healthcare professionals in this regard. These findings provide additional validation for the evidence-based suggestion to implement HPV PCR screening across the general population, as endorsed by scientific associations [61]. In 2019, the Ministry of Health, Consumption, and Welfare modified the Royal Decree on common services of the National Health System, through which the autonomous communities were gradually required to replace cytological screening with HPV testing and opportunistic screening with population-based screening [13].
Individuals who engage in risky sexual behavior, particularly in relation to unwanted pregnancies, are also at a heightened risk of contracting STDs and HPV [62]. Furthermore, there is an observed correlation between cervical dysplasia and sexually transmitted infections, particularly among younger patients [63]. The results suggest that there may be an underdiagnosis of these diseases. Individuals who face barriers to accessing healthcare services may encounter difficulties in obtaining screening for sexually transmitted diseases and cervical cancer. Consequently, they may miss out on vital primary prevention measures such as vaccination, as well as important secondary prevention measures such as regular screening for cervical cancer.
In order to improve protection against HPV-related diseases, it is recommended to complete the full vaccination schedule, which typically consists of two or three doses. However, the study highlights that receiving the first dose of vaccination increases the likelihood of completing the full schedule by 25 times. Studies [64,65,66] have demonstrated that individuals who initiate the HPV vaccination series are more likely to complete the entire schedule if they receive a recommendation from their healthcare provider and have a reminder system in place. Furthermore, these coverage results still have room for improvement by reducing the number of doses, which appears to enhance adherence to the regimen [67]. The WHO has published a document that enables better access to vaccination by reducing vaccine doses. This way, it could be used as a single-dose vaccine in individuals aged 9 to 20 years [68].
Regarding the remaining variables, “offspring” and “breastfeeding,” no previous evidence has been found. Further research is needed on this topic. As for immune status, modifications in the screening interval should be considered based on the woman’s medical history and her risk factors [69], due to their heightened vulnerability to HPV infection [41]. The role of smoking is debated in numerous studies [22,23,24], but our study did not find smoking to be a prognostic factor for HPV vaccination.
There are some limitations, such as the lack of data on the male population and potential limitations due to a registration bias in some variables in addition to the data sources selected. In Spain, it was in the year 2008 that HPV vaccines began to be introduced into the vaccination schedules of the autonomous communities, and it was in the year 2023 that vaccination was included for 12-year-old male adolescents. The results suggest a possible relationship between HPV vaccination and sexual and reproductive health behaviors but do not imply causation. Additional research would be necessary to determine how changes in vaccination coverage correlate with accessibility to sexual and reproductive health services and their potential interplay with the age and gender of immigrant individuals, especially among newcomers and due to the new indication of the HPV vaccine for males.

5. Conclusions

The study identified as barriers age (OR 0.68, 95% CI 0.63–0.72, p < 0.001) and immigrant origin (OR 0.21, 95% CI 0.10–0.42, p < 0.001), which were found to decrease the likelihood of HPV vaccination, and HPV PCR (OR 7.40, 95% CI 2.50–21.88, p < 0.001) as a facilitator for HPV vaccination. Also, the accessibility to the initial vaccine greatly enhances the likelihood of completing the full vaccination schedule, increasing it by a remarkable 25 times. These findings highlight opportunities to improve HPV vaccination rates.

Author Contributions

Conceptualization, L.C.-C. and J.L.C.-E.; methodology, L.C.-C., J.F.-S., M.J.P.-F. and J.L.C.-E.; software, C.C.-G., J.C.-Q. and J.F.-S.; validation, L.C.-C., C.C.-G., J.C.-Q., J.F.-S., M.J.P.-F. and J.L.C.-E.; formal analysis, L.C.-C., J.F.-S. and J.L.C.-E.; investigation, L.C.-C., C.C.-G., J.C.-Q., J.F.-S., M.J.P.-F. and J.L.C.-E.; resources, C.C.-G. and J.C.-Q.; data curation, L.C.-C. and J.F.-S.; writing—original draft preparation, L.C.-C., C.C.-G., J.C.-Q., J.F.-S., M.J.P.-F. and J.L.C.-E.; writing—review and editing, L.C.-C., C.C.-G., J.C.-Q., J.F.-S., M.J.P.-F. and J.L.C.-E.; visualization, L.C.-C. and J.F.-S.; supervision, L.C.-C., C.C.-G., J.C.-Q., J.F.-S., M.J.P.-F. and J.L.C.-E.; project administration, L.C.-C. and J.L.C.-E.; funding acquisition, L.C.-C. and C.C.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of IDIAP Jordi Gol University Institute of Primary Care Research (protocol code, 21/064-P; date of approval, 24 March 2021).

Informed Consent Statement

Patient consent was waived prior to the inclusion of medical data since formal consent is not required for this type of study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author (L.C.-C.) upon reasonable request.

Acknowledgments

The authors would like to thank Tom Foy for his technical support and Staci L Sudenga for her invaluable mentorship.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Demographic Characteristics of the Region

Vaccines 11 01536 i001

Appendix B. Immune Status by Vaccination Status

Immune StatusNon-VaccinatedVaccinatedAll
n (%)n (%)n (%)
Alterations5 (1.5)6 (3.3)11 (2.1)p = 0.170
No Alteration333 (98.5)176 (96.7)509 (97.8)

Appendix C. Tobacco Use by Vaccination Status

Vaccination Status
Age GroupsNon-VaccinatedVaccinatedAll
n (%)n (%)n (%)
15–190 (0.0)1 (1.2)1 (0.1)p < 0.001
20–240 (0.0)10 (15.2)10 (1.9)
25–297 (13.7)3 (17.6)10 (1.9)
30–3414 (12.8)3 (23.1)17 (3.2)
>3538 (24.4)2 (33.3)40 (7.6)
All59 (11.3)19 (3.6)78 (15.0)

References

  1. World Health Organization. Cervical Cancer. Available online: https://www.who.int/news-room/fact-sheets/detail/cervical-cancer (accessed on 23 January 2023).
  2. World Health Organization. Reaching 2030 Cervical Cancer Elimination Targets—New WHO Recommendations for Screening and Treatment of Cervical Pre-Cancer. Available online: https://www.who.int/news-room/events/detail/2021/07/06/default-calendar/reaching-2030-cervical-cancer-elimination-targets (accessed on 25 June 2023).
  3. Bruni, L.; Saura-Lázaro, A.; Montoliu, A.; Brotons, M.; Alemany, L.; Diallo, M.S.; Afsar, O.Z.; LaMontagne, D.S.; Mosina, L.; Contreras, M.; et al. HPV vaccination introduction worldwide and WHO and UNICEF estimates of national HPV immunization coverage 2010–2019. Prev. Med. 2021, 144, 106399. [Google Scholar] [CrossRef] [PubMed]
  4. Braeye, T.; Emborg, H.-D.; Llorente-García, A.; Huerta, C.; Martín-Merino, E.; Duarte-Salles, T.; Danieli, G.; Tramontan, L.; Weibel, D.; McGee, C.; et al. Age-specific vaccination coverage estimates for influenza, human papillomavirus and measles containing vaccines from seven population-based healthcare databases from four EU countries—The ADVANCE project. Vaccine 2020, 38, 3243–3254. [Google Scholar] [CrossRef] [PubMed]
  5. El Proceso de Vaciado de la España Interior No se Detiene, Advierten los Expertos. Available online: https://www.elperiodicodearagon.com/sociedad/2023/01/09/proceso-vaciado-espana-interior-detiene-80864275.html (accessed on 18 April 2023).
  6. La España Vaciada: La Decadencia de la Agricultura y la Aparición de Nuevos Partidos. Available online: https://diario16.com/la-espana-vaciada-la-decadencia-de-la-agricultura-y-la-aparicion-de-nuevos-partidos/ (accessed on 18 April 2023).
  7. Idescat. Anuario Estadístico de Cataluña. Densidad de Población. Comarcas y Aran, Ámbitos y Provincias. Available online: https://www.idescat.cat/pub/?id=aec&n=249&lang=es (accessed on 18 April 2023).
  8. Idescat. Indicadors Demogràfics i de Territori. Estructura per Edats, Envelliment i Dependència. Comarques i Aran. Available online: http://www.idescat.cat/pub/?id=inddt&n=915&by=com (accessed on 18 April 2023).
  9. Agència de Qualitat i Avaluació Sanitàries de Catalunya, Generalitat de Catalunya. Què Destaquem de la Regió Sanitària de Camp de Tarragona? Observatori del Sistema de Salut de Catalunya. Generalitat de Catalunya. Available online: https://observatorisalut.gencat.cat/web/.content/minisite/observatorisalut/ossc_central_resultats/Presentacions/Presentacio-RSTarragona-i-RSTE_def.pdf (accessed on 18 April 2023).
  10. Idescat. Anuari Estadístic De Catalunya. Renda Familiar Disponible Bruta. Índex. Comarques i Aran, i Àmbits. Available online: http://www.idescat.cat/pub/?id=aec&n=941 (accessed on 18 April 2023).
  11. Colomé-Ceballos, L.; Clua-Espuny, J.L.; Fernández-Sáez, J.; Ceballos-García, C.; Andrés-Cubells, N.; Pla-Farnós, M.J. HPV Vaccination Coverage Rate in a Rural Area: An Observational, Retrospective, and Cohort Study. Vaccines 2022, 10, 1274. [Google Scholar] [CrossRef] [PubMed]
  12. Protocol de les Activitats per al Cribratge del Càncer de Coll Uterí a l’Atenció Primària. Available online: https://canalsalut.gencat.cat/web/.content/_A-Z/C/cancer/recursos_prof/deteccio_precoc/protocol_cribratge_cancer_coll_uter.pdf (accessed on 23 April 2023).
  13. Ministry of Health, Consumption and Social Welfare. Order SCB/480/2019, of April 26, Amending Annexes I, II, and IV of Royal Decree 1030/2006, of September 15, Which Establishes the Common Services Portfolio of the National Health System and the Procedure for Its Updating. Official State Gazette Number 101, Dated 27 April 2019, pp. 43018–43028. Available online: https://www.boe.es/boe/dias/2019/04/27/pdfs/BOE-A-2019-6277.pdf (accessed on 20 June 2023).
  14. Agència de Salut Pública de Catalunya (ASPCAT). Available online: https://salutpublica.gencat.cat/ca/ambits/promocio_salut/determinants_de_salut/ (accessed on 23 January 2023).
  15. Spayne, J.; Hesketh, T. Estimate of global human papillomavirus vaccination coverage: Analysis of country-level indicators. BMJ Open 2021, 11, e052016. [Google Scholar] [CrossRef] [PubMed]
  16. Gallagher, K.; LaMontagne, D.; Watson-Jones, D. Status of HPV vaccine introduction and barriers to country uptake. Vaccine 2018, 36 Pt A, 4761–4767. [Google Scholar] [CrossRef]
  17. Vu, M.; King, A.R.; Jang, H.M.; Bednarczyk, R.A. Practice-, provider- and patient-level facilitators of and barriers to HPV vaccine promotion and uptake in Georgia: A qualitative study of healthcare providers’ perspectives. Health Educ. Res. 2020, 35, 512–523. [Google Scholar] [CrossRef]
  18. Robles, C.; Bruni, L.; Acera, A.; Riera, J.C.; Prats, L.; Poljak, M.; Mlakar, J.; Valenčak, A.O.; Eriksson, T.; Lehtinen, M.; et al. Determinants of Human Papillomavirus Vaccine Uptake by Adult Women Attending Cervical Cancer Screening in 9 European Countries. Am. J. Prev. Med. 2021, 60, 478–487. [Google Scholar] [CrossRef]
  19. Wemrell, M.; Vicente, R.P.; Merlo, J. Mapping sociodemographic and geographical differences in human papillomavirus non-vaccination among young girls in Sweden. Scand. J. Public Health 2022, 51, 288–295. [Google Scholar] [CrossRef]
  20. Groene, E.A.; Horvath, K.J.; Yared, N.; Mohammed, I.; Muscoplat, M.; Kuramoto, S.; Richter, T.; Kulasingam, S. Missed Opportunities for Human Papillomavirus Vaccination by Parental Nativity, Minnesota, 2015–2018. Public Health Rep. 2021, 137, 867–877. [Google Scholar] [CrossRef]
  21. Cofie, L.E.; Hirth, J.M.; Guo, F.; Berenson, A.B.; Markides, K.; Wong, R. HPV Vaccination Among Foreign-Born Women: Examining the National Health Interview Survey 2013–2015. Am. J. Prev. Med. 2018, 54, 20–27. [Google Scholar] [CrossRef]
  22. Newman, P.A.; Logie, C.H.; Doukas, N.; Asakura, K. HPV vaccine acceptability among men: A systematic review and meta-analysis. Sex. Transm. Infect. 2013, 89, 568–574. [Google Scholar] [CrossRef]
  23. Manhart, L.E.; Burgess-Hull, A.J.; Fleming, C.B.; Bailey, J.A.; Haggerty, K.P.; Catalano, R.F. HPV vaccination among a community sample of young adult women. Vaccine 2011, 29, 5238–5244. [Google Scholar] [CrossRef]
  24. Nielsen, A.; Munk, C.; Liaw, K.-L.; Kjær, S.K. Awareness of human papillomavirus in 23,000 Danish men from the general male population. Eur. J. Cancer Prev. 2009, 18, 236–239. [Google Scholar] [CrossRef]
  25. Gunn, R.; Ferrara, L.K.; Dickinson, C.; Stock, I.; Griffith-Weprin, J.; Wiser, A.; Hatch, B.; Fagnan, L.; Carney, P.A.; Davis, M.M. Human Papillomavirus Immunization in Rural Primary Care. Am. J. Prev. Med. 2020, 59, 377–385. [Google Scholar] [CrossRef]
  26. Patel, D.A.; Zochowski, M.K.; Peterman, S.; Dempsey, A.F.; Bernstein, D.I. Disparities in HPV vaccine awareness, knowledge, and receipt among rural compared to urban women in the United States, 2007–2016. Vaccine 2018, 36, 4998–5004. [Google Scholar]
  27. Mohammed, K.A.; Subramaniam, D.S.; Geneus, C.J.; Henderson, E.R.; Dean, C.A.; Subramaniam, D.P.; Burroughs, T.E. Rural-urban differences in human papillomavirus knowledge and awareness among US adults. Prev. Med. 2018, 109, 39–43. [Google Scholar] [CrossRef]
  28. Goessl, C.L.; Christianson, B.; Hanson, K.E.; Polter, E.J.; Olson, S.C.; Boyce, T.G.; Dunn, D.; Williams, C.L.; Belongia, E.A.; McLean, H.Q.; et al. Human papillomavirus vaccine beliefs and practice characteristics in rural and urban adolescent care providers. BMC Public Health 2022, 22, 1322. [Google Scholar] [CrossRef]
  29. McRee, A.-L.; Gilkey, M.B.; Dempsey, A.F. HPV Vaccine Hesitancy: Findings from a Statewide Survey of Health Care Providers. J. Pediatr. Health Care 2014, 28, 541–549. [Google Scholar] [CrossRef]
  30. Dempsey, A.F.; Butchart, A.; Singer, D.; Clark, S.; Davis, M. Factors Associated with Parental Intentions for Male Human Papillomavirus Vaccination: Results of a National Survey. Sex. Transm. Dis. 2011, 38, 769–776. [Google Scholar] [CrossRef]
  31. Jena, A.B.; Goldman, D.P.; Seabury, S.A. Incidence of Sexually Transmitted Infections After Human Papillomavirus Vaccination Among Adolescent Females. JAMA Intern. Med. 2015, 175, 617–623. [Google Scholar] [CrossRef]
  32. Kester, L.M.; Zimet, G.D.; Fortenberry, J.D.; Kahn, J.A.; Shew, M.L. A national study of HPV vaccination of adolescent girls: Rates, predictors, and reasons for non-vaccination. Matern. Child Health J. 2012, 17, 879–885. [Google Scholar] [CrossRef]
  33. Els Nois També Rebran la Vacuna del Papilloma Humà. Available online: https://web.gencat.cat/ca/actualitat/detall/Els-nois-tambe-rebran-la-vacuna-del-papilloma-huma (accessed on 23 April 2023).
  34. Álvarez García, F.J.; Cilleruelo Ortega, M.J.; Álvarez Aldeán, J.; Garcés-Sánchez, M.; García Sánchez, N.; Garrote Llanos, E.; Hernández Merino, Á.; Iofrío de Arce, A.; Merino Moína, M.; Montesdeoca Melián, A.; et al. Immunisation schedule of the Spanish Association of Paediatrics: 2020 recommendations. Ann. Pediatr. Engl. Ed. 2020, 92, 52.e1–52.e10. [Google Scholar] [CrossRef]
  35. Generalitat de Catalunya. Projeccions de Població Principals Resultats 2013–2051. 2008. Available online: https://www.idescat.cat/serveis/biblioteca/docs/cat/pp2021-2041pr.pdf (accessed on 14 April 2023).
  36. WHO Updates Recommendations on HPV Vaccination Schedule. Available online: https://www.who.int/news/item/20-12-2022-WHO-updates-recommendations-on-HPV-vaccination-schedule (accessed on 18 April 2023).
  37. HPV Vaccination Recommendations. Available online: https://www.cdc.gov/vaccines/vpd/hpv/hcp/recommendations.html (accessed on 18 April 2023).
  38. World Health Organization. Strategic Advisory Group of Experts on Immunization. Available online: https://www.who.int/groups/strategic-advisory-group-of-experts-on-immunization (accessed on 18 April 2023).
  39. Pan American Health Organization. One-Dose Human Papillomavirus (HPV) Vaccine Offers Solid Protection against Cervical Cancer. Available online: https://www.paho.org/en/news/11-4-2022-one-dose-human-papillomavirus-hpv-vaccine-offers-solid-protection-against-cervical (accessed on 18 April 2023).
  40. Departament de Salut. Manual de Vacunacions de Catalunya; Agència de Salut Pública de Catalunya: Barcelona, Spain, 2020. [Google Scholar]
  41. Campins, M.T.A.; Alemany, L.; Bayas, J.M.; Borruel, N.; Castellsagué, X.; Curran, A.; Diaz Heredia, C.; Martinez, X.; Moraga llop, F.A. AEPCC-Guía: Vacunación Selectiva Frente al Virus del Papiloma Humano en Poblaciones de Riesgo Elevado; Publicaciones AEPCC: Madrid, Spain, 2016; pp. 1–46. [Google Scholar]
  42. Pla de Salut de la Regió Sanitària Terres de l’Ebre 2021–2025. Available online: https://scientiasalut.gencat.cat/handle/11351/7964 (accessed on 2 May 2023).
  43. Henry, K.A.; Stroup, A.M.; Warner, E.L.; Kepka, D. Geographic Factors and Human Papillomavirus (HPV) Vaccination Initiation among Adolescent Girls in the United States. Cancer Epidemiol. Biomark. Prev. 2016, 25, 309–317. [Google Scholar] [CrossRef] [PubMed]
  44. Henry, K.A.; Swiecki-Sikora, A.L.; Stroup, A.M.; Warner, E.L.; Kepka, D. Area-based socioeconomic factors and Human Papillomavirus (HPV) vaccination among teen boys in the United States. BMC Public Health 2017, 18, 19. [Google Scholar] [CrossRef] [PubMed]
  45. Sull, M.; Eavey, J.; Papadouka, V.; Mandell, R.; Hansen, M.A.; Zucker, J.R. Adolescent Vaccine Co-administration and Coverage in New York City: 2007–2013. Pediatrics 2014, 134, e1576–e1583. [Google Scholar] [CrossRef]
  46. Eberth, J.M.; Hossain, M.; Tiro, J.A.; Zhang, X.; Holt, J.B.; Vernon, S.W. Human Papillomavirus Vaccine Coverage Among Females Aged 11 to 17 in Texas Counties: An Application of Multilevel, Small Area Estimation. Women’s Health Issues 2013, 23, e131–e141. [Google Scholar] [CrossRef]
  47. Pruitt, S.L.; Schootman, M. Geographic Disparity, Area Poverty, and Human Papillomavirus Vaccination. Am. J. Prev. Med. 2010, 38, 525–533. [Google Scholar] [CrossRef]
  48. Elenwo, C.; Batioja, K.; Davis, T.; Greiner, B.; Markey, C.; Hartwell, M. Associations of Maternal Age, Education, and Marital Status with HPV Vaccine Uptake and Hesitancy among United States Youth: A Cross-Sectional Analysis of the 2020 National Immunization Survey. J. Pediatr. Adolesc. Gynecol. 2023, 36, 273–279. [Google Scholar] [CrossRef]
  49. Manual de Vacunacions de Catalunya. 2020. Available online: https://salutpublica.gencat.cat/ca/ambits/promocio_salut/vacunacions/Manual-de-vacunacions/Versio-desglossada-per-capitols/ (accessed on 25 June 2023).
  50. Idescat. Anuario Estadístico de Cataluña. Población Extranjera. Comarcas y Aran, Ámbitos y Provincias. Available online: https://www.idescat.cat/poblacioestrangera/?b=4 (accessed on 23 April 2023).
  51. Wilson, L.; Rubens-Augustson, T.; Murphy, M.; Jardine, C.; Crowcroft, N.; Hui, C.; Wilson, K. Barriers to immunization among newcomers: A systematic review. Vaccine 2018, 36, 1055–1062. [Google Scholar] [CrossRef]
  52. Cofie, L.E.; Tailor, H.D.; Lee, M.H.; Xu, L. HPV vaccination uptake among foreign-born Blacks in the US: Insights from the National Health Interview Survey 2013–2017. Cancer Causes Control 2022, 33, 583–591. [Google Scholar] [CrossRef]
  53. Netfa, F.; Tashani, M.; Booy, R.; King, C.; Rashid, H.; Skinner, S.R. Knowledge, Attitudes and Perceptions of Immigrant Parents Towards Human Papillomavirus (HPV) Vaccination: A Systematic Review. Trop. Med. Infect. Dis. 2020, 5, 58. [Google Scholar] [CrossRef] [PubMed]
  54. Ratnasamy, P.; Chagpar, A.B. HPV vaccination and factors influencing vaccine uptake among people of Indian ancestry living in the United States. Epidemiol. Infect. 2022, 150, e152. [Google Scholar] [CrossRef] [PubMed]
  55. White, A.A.; Neelon, B.; Martin, R.H.; Roberts, J.R.; Korte, J.E.; Williams, E.M.; Cartmell, K.B. Predicting HPV vaccination among Tdap vaccinated adolescents in Georgia at the county level. Ann. Epidemiol. 2022, 70, 74–78. [Google Scholar] [CrossRef] [PubMed]
  56. Warner, E.L.; Ding, Q.; Pappas, L.M.; Henry, K.; Kepka, D. White, affluent, educated parents are least likely to choose HPV vaccination for their children: A cross-sectional study of the National Immunization Study—Teen. BMC Pediatr. 2017, 17, 200. [Google Scholar] [CrossRef]
  57. Javanbakht, M.; Stahlman, S.; Walker, S.; Gottlieb, S.; Markowitz, L.; Liddon, N.; Plant, A.; Guerry, S. Provider perceptions of barriers and facilitators of HPV vaccination in a high-risk community. Vaccine 2012, 30, 4511–4516. [Google Scholar] [CrossRef]
  58. Allen, E.M.; Lee, H.Y.; Pratt, R.; Vang, H.; Desai, J.R.; Dube, A.; Lightfoot, E. Facilitators and Barriers of Cervical Cancer Screening and Human Papilloma Virus Vaccination Among Somali Refugee Women in the United States: A Qualitative Analysis. J. Transcult. Nurs. 2018, 30, 55–63. [Google Scholar] [CrossRef]
  59. Stratoberdha, D.; Gobis, B.; Ziemczonek, A.; Yuen, J.; Giang, A.; Zed, P.J. Barriers to adult vaccination in Canada: A qualitative systematic review. Can. Pharm. J. 2022, 155, 206–218. [Google Scholar] [CrossRef]
  60. Torné Bladé, A.; Del Pino Saladrigues, M.; Cusidó Gimferrer, M.; Alameda Quitllet, F.; Andía Ortiz, D.; Castellsagué Piqué, X.; Cortés Bordoy, J.; Granados Carreño, R. Guia del Cribado de Cáncer de Cuello en España, 1st ed.; Publicaciones AEPCC: Madrid, Spain, 2014. [Google Scholar]
  61. Torné, A.; Andía, D.; Bruni, L.; Centeno, C.; Coronado, P.; Cruz Quílez, J.; de la Fuente, J.; de Sanjosé, S.; Granados, R.; Ibáñez, R.; et al. AEPCC-Guía: Prevención Secundaria del Cáncer de Cuello de Útero, 1st ed.; Publicaciones AEPCC: Madrid, Spain, 2022. [Google Scholar]
  62. D’Souza, G.; Wentz, A.; Kluz, N.; Zhang, Y.; Sugar, E.; Youngfellow, R.M.; Guo, Y.; Xiao, W.; Gillison, M.L. Sex Differences in Risk Factors and Natural History of Oral Human Papillomavirus Infection. J. Infect. Dis. 2016, 213, 1893–1896. [Google Scholar] [CrossRef]
  63. Adhikari, I.; Eriksson, T.; Harjula, K.; Hokkanen, M.; Apter, D.; Nieminen, P.; Luostarinen, T.; Lehtinen, M. Association of Chlamydia trachomatis infection with cervical atypia in adolescent women with short-term or long-term use of oral contraceptives: A longitudinal study in HPV vaccinated women. BMJ Open 2022, 12, e056824. [Google Scholar] [CrossRef]
  64. Drolet, M.; Bénard, É.; Boily, M.-C.; Ali, H.; Baandrup, L.; Bauer, H.; Beddows, S.; Brisson, J.; Brotherton, J.M.L.; Cummings, T.; et al. Population-level impact and herd effects following human papillomavirus vaccination programmes: A systematic review and meta-analysis. Lancet Infect. Dis. 2015, 15, 565–580. [Google Scholar] [CrossRef]
  65. Drolet, M.; Bénard, É.; Pérez, N.; Brisson, M.; HPV Vaccination Impact Study Group. Population-level impact and herd effects following the introduction of human papillomavirus vaccination programmes: Updated systematic review and meta-analysis. Lancet 2019, 394, 497–509. [Google Scholar] [CrossRef] [PubMed]
  66. Gallagher, K.E.; Kadokura, E.; Eckert, L.O.; Miyake, S.; Mounier-Jack, S.; Aldea, M.; Ross, D.A.; Watson-Jones, D. Factors influencing completion of multi-dose vaccine schedules in adolescents: A systematic review. BMC Public Health 2016, 16, 172. [Google Scholar] [CrossRef] [PubMed]
  67. Sankaranarayanan, R.; Joshi, S.; Muwonge, R.; Esmy, P.O.; Basu, P.; Prabhu, P.; Bhatla, N.; Nene, B.M.; Shaw, J.; Poli, U.R.R.; et al. Can a single dose of human papillomavirus (HPV) vaccine prevent cervical cancer? Early findings from an Indian study. Vaccine 2018, 36, 4783–4791. [Google Scholar] [CrossRef] [PubMed]
  68. International Agency for Research on Cancer. Protection from a Single Dose of HPV Vaccine: A Major Public Health Impact from IARC Studies of Vaccine Efficacy (IARC Evidence Summary Briefs, No. 4); IARC: Lyon, France, 2023; Available online: https://www.iarc.who.int/evidence-summary-briefs-series/ (accessed on 20 June 2023).
  69. SAGE. Strategic Advisory Group of Experts (SAGE) Working Group on Potential Contribution of HPV Vaccines and Immunization Towards Cervical Cancer Elimination. 2022. Available online: https://cdn.who.int/media/docs/default-source/immunization/position_paper_documents/human-papillomavirus-(hpv)/hpv-background-document--report-march-2022.pdf?sfvrsn=b600e252_1 (accessed on 20 June 2023).
Vaccines 11 01536 g001
Figure 1. ROC curve of the logistic regression statistical model.
Figure 1. ROC curve of the logistic regression statistical model.
Vaccines 11 01536 g001
Vaccines 11 01536 g002
Figure 2. Probability distribution of the forecast from the logistic regression model.
Figure 2. Probability distribution of the forecast from the logistic regression model.
Vaccines 11 01536 g002
Table 1. Vaccination coverage rate by age groups.
Table 1. Vaccination coverage rate by age groups.
Age GroupsNon-VaccinatedVaccinatedAllp
n%n%n%
15–1967.08093.08616.5<0.001
20–241619.56680.58215.8
25–295175.01725.06813.1
30–3410989.31310.712223.5
>3515696.363.716231.2
All33865.018235.0520100.0
Average age
(SD 1) 2		33.1 (5.9)21.8 (5.2)29.3 (7.8)<0.001
1 Standard deviation; 2 in years old.
Table 2. Demographic variables by vaccination status.
Table 2. Demographic variables by vaccination status.
Non-VaccinatedVaccinatedAllp = 0.044
Ethnic originn (%)n (%)n (%)
European277 (62.7)165 (37.3)442 (85.0)
African35 (83.3)7 (16.7)42 (8.1)
American19 (73.1)7 (26.9)26 (5.0)
Asian7 (70.0)3 (30.0)10 (1.9)
Non-VaccinatedVaccinatedAllp < 0.001
Immigrant originn (%)n (%)n (%)
Non-Immigrant origin228 (43.8)150 (28.8)378 (72.7)
Immigrant origin110 (21.1)32 (6.1)142 (27.3)
Non-VaccinatedVaccinatedAllp = 0.081
Educationn (%)n (%)n (%)
No education21 (91.3)2 (8.7)23 (4.4)
Compulsory education90 (73.2)33 (26.8)123 (23.7)
Higher education66 (82.5)14 (17.5)80 (15.4)
Table 3. Cervical cancer screening test by vaccination status.
Table 3. Cervical cancer screening test by vaccination status.
Cervical Cytology TestNon-VaccinatedVaccinatedAllp < 0.001
n (%)n (%)n (%)
Performed140 (41.4)28 (15.3)168 (32.3)
Not Performed198 (58.6)154 (84.6)352 (67.7)
Results of Cervical Cytology TestNon-vaccinatedVaccinatedAllp < 0.001
n (%)n (%)n (%) *
Negative136 (40.2)19 (10.4)155 (92.3)
Positive4 (1.2)9 (4.9)13 (7.7)
HPV 1 PCR 2 testNon-vaccinatedVaccinatedAllp = 0.885
n (%)n (%)n (%)
Performed23 (6.8)13 (7.1)36 (6.9)
Not Performed315 (93.2)169 (92.9)484 (93.1)
HPV 1 PCR 2 testNon-vaccinatedVaccinatedAllp = 0.047
n (%)n (%)n (%) *
Negative21 (87.5)7 (58.3)28 (77.8)
Positive3 (12.5)5 (41.7)8 (22.2)
* From the total of those performed; 1 human papillomavirus; 2 polymerase chain reaction.
Table 4. Sexual transmitted disease screening test by vaccination status.
Table 4. Sexual transmitted disease screening test by vaccination status.
STD 1 PCR 2 testNon-VaccinatedVaccinatedAllp = 0.760
n (%)n (%)n (%)
Performed20 (5.9)12 (6.6)32 (6.2)
Not Performed318 (94.1)170 (93.4)488 (93.8)
Results of STD 1 PCR 2 testNon-vaccinatedVaccinatedAllp = 0.555
n (%)n (%)n (%) *
Negative16 (76.2)8 (66.7)24 (72.7)
Positive5 (23.8)4 (33.3)9 (27.3)
1 Sexually transmitted disease; 2 polymerase chain reaction; * from the total of those performed.
Table 5. Variables of sexual and reproductive health by vaccination status.
Table 5. Variables of sexual and reproductive health by vaccination status.
OffspringNon-VaccinatedVaccinatedAllp < 0.001
n (%)n (%)n (%)
No Offspring166 (50.3)164 (49.7)330 (63.4)
Offspring171 (90.5)19 (9.5)189 (36.3)
BreastfeedingNon-VaccinatedVaccinatedAllp < 0.001
n (%)n (%)n (%)
No Breastfeeding227 (56.8)173 (43.2)400 (76.9)
Breastfeeding111 (92.5)9 (7.5)120 (23.1)
AbortionNon-VaccinatedVaccinatedAllp = 0.729
n (%)n (%)n (%)
No Abortion326 (65.1)175 (34.9)502 (96.5)
Abortion11 (61.1)7 (38.9)18 (3.5)
ContraceptionNon-VaccinatedVaccinatedAllp = 0.008
n (%)n (%)n (%)
No Contraception40 (80.0)10 (20.0)50 (9.6)
Contraception76 (58.9)53 (41.1)129 (24.8)
Condom UseNon-VaccinatedVaccinatedAllp = 0.110
n (%)n (%)n (%)
No Condom Use136 (72.3)52 (27.7)188 (36.1)
Condom Use50 (62.5)30 (37.5)80 (15.4)
Table 6. Logistic regression analysis.
Table 6. Logistic regression analysis.
VariablepOR95% CI
Lower LimitUpper Limit
Age<0.0010.6740.6250.726
Cervical cytology test performed0.0632.3150.9555.613
HPV PCR performed0.0017.3372.35222.888
STD PCR performed0.3370.5690.1801.797
Offspring0.5161.4200.4934.086
Breastfeeding0.2430.4770.1381.653
Smoking0.4621.3980.5723.415
Contraception0.4460.6030.1652.210
Immune status0.1423.9660.62925.005
Ethnic origin0.3880.5230.1202.280
Immigrant origin0.0110.2980.1180.754
Results statistical model of logistic regression
Age<0.0010.6800.6350.729
Cervical cytology test performed0.0782.1350.9174.970
HPV PCR performed<0.0017.4022.50421.880
Immigrant origin<0.0010.2150.1090.422
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Colomé-Ceballos, L.; Clua-Espuny, J.L.; Ceballos-García, C.; Clua-Queralt, J.; Pla-Farnós, M.J.; Fernández-Sáez, J. Understanding Prognostic Factors for Human Papillomavirus Vaccination: A Rural Community Case–Control Study. Vaccines 2023, 11, 1536. https://doi.org/10.3390/vaccines11101536

AMA Style

Colomé-Ceballos L, Clua-Espuny JL, Ceballos-García C, Clua-Queralt J, Pla-Farnós MJ, Fernández-Sáez J. Understanding Prognostic Factors for Human Papillomavirus Vaccination: A Rural Community Case–Control Study. Vaccines. 2023; 11(10):1536. https://doi.org/10.3390/vaccines11101536

Chicago/Turabian Style

Colomé-Ceballos, Lara, Josep Lluís Clua-Espuny, Concepción Ceballos-García, Josep Clua-Queralt, Maria Jesús Pla-Farnós, and José Fernández-Sáez. 2023. "Understanding Prognostic Factors for Human Papillomavirus Vaccination: A Rural Community Case–Control Study" Vaccines 11, no. 10: 1536. https://doi.org/10.3390/vaccines11101536

APA Style

Colomé-Ceballos, L., Clua-Espuny, J. L., Ceballos-García, C., Clua-Queralt, J., Pla-Farnós, M. J., & Fernández-Sáez, J. (2023). Understanding Prognostic Factors for Human Papillomavirus Vaccination: A Rural Community Case–Control Study. Vaccines, 11(10), 1536. https://doi.org/10.3390/vaccines11101536

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop