HPV vaccine-related thyroid adverse events: temporal patterns and reporting trends

Background: Human papillomavirus (HPV) vaccination is essential for cervical cancer prevention, but concerns about thyroid-related adverse events (AEs) have emerged.

Methods: This pharmacovigilance study aimed to assess potential associations between HPV vaccination and thyroid disorders using spontaneous report data. Reports from the U.S. Vaccine Adverse Event Reporting System through 31 December 2024, were analyzed. Disproportionality analyses (PRR, ROR, BCPNN, MGPS) were performed, with subgroup analyses by gender, age, and vaccine type. A Weibull shape parameter model assessed the temporal risk pattern.

Results: Among 60,840 HPV vaccine-related reports, 13 thyroid-associated AEs showed positive signals. Hypothyroidism (ROR = 11.65) and autoimmune thyroiditis (ROR = 4.26) were the strongest signals. Most cases occurred in females under 65 years. HPV-4 was linked to 72.8% of thyroid AEs. The cumulative reporting rate reached 69.1% within 180 days.

Conclusion: Our findings suggest a potential association between HPV vaccination and thyroid disorders, notably hypothyroidism and autoimmune thyroiditis. Continued pharmacovigilance and further mechanistic investigations are needed.

1 Introduction

Human papillomavirus (HPV) is a double-stranded DNA virus that infects epithelial cells and evades immune surveillance by modulating host immune responses, thereby promoting oncogenesis (Steinbach and Riemer, 2018; Yousefi et al., 2021). Clinical manifestations encompass benign lesions caused by low-risk carcinogenic HPV subtypes (e.g., genital, anal, oral, or laryngeal warts) and malignancies associated with high-risk carcinogenic subtypes (e.g., vulvovaginal, oropharyngeal, cervical, and head/neck cancers) (Sharma et al., 2023; Yousaf et al., 2024).

HPV vaccination is currently one of the most effective preventive measures against cervical cancer. According to World Health Organization (WHO) data, 133 countries/regions have incorporated HPV vaccines into their National Immunization Programs (NIPs) (Yousaf et al., 2024). Three prophylactic HPV vaccines are commercially available (2vHPV, 4vHPV, 9vHPV), all containing the major viral capsid protein L1 derived from specific HPV subtypes. These proteins are produced via recombinant DNA technology and self-assemble into virus-like particles (VLPs). VLPs lack viral DNA and live infectious particles, rendering them non-infectious and non-carcinogenic. The vaccines utilize L1 protein as target antigens, inducing humoral immunity to generate virus-neutralizing antibodies with 10- to 100-fold higher titers than natural infection, thereby effectively neutralizing virions (Yousefi et al., 2021).

Although the vaccines are generally safe, rare adverse events have been reported, including autoimmune conditions such as premature ovarian insufficiency (POI), Guillain–Barré syndrome, and complex regional pain syndrome (Tsai et al., 2023). Emerging studies have suggested potential associations between HPV vaccination and thyroid disorders (Willame et al., 2020). However, findings are inconsistent, with some studies suggesting increased risk, while others report null or even protective associations (Grimaldi-Bensouda et al., 2017; Hviid et al., 2018; Rosillon et al., 2020; Yoon et al., 2021).

The Vaccine Adverse Event Reporting System (VAERS), a national spontaneous surveillance system, accepts reports of potential vaccine-related adverse events (AEs) from healthcare professionals, vaccine recipients, manufacturers, and other reporters (Vaccine Adverse Events Reporting System. VAERS data use guide). Trained medical coders systematically review all reports, classifying signs and symptoms using Medical Dictionary for Regulatory Activities (MedDRA), a clinically validated international standardized terminology (Shimabukuro et al., 2015). Based on this, data mining analysis of adverse reactions using the VAERS database may provide a preliminary validation strategy for the study of the relationship between HPV vaccines and AEs.

This study aimed to systematically analyze thyroid-related adverse events reported after HPV vaccination in VAERS, focusing on their frequency and temporal patterns, to provide evidence supporting future research and public health risk assessment.

2 Materials and methods

VAERS was established in 1990 as a national spontaneous reporting system for post-vaccination adverse events. It is jointly administered by the Centers for Disease Control and Prevention (CDC) and the Food and Drug Administration (FDA) in the United States (Sever et al., 2004; Moro et al., 2023). Reports are submitted to VAERS by patients, clinicians, regulators, and others. It is imperative that submissions accurately reflect the information regarding the administered vaccine, the individual who received the vaccine, and the adverse events encountered by the vaccinated person (Chen et al., 2024). Reports are categorised as either serious or non-serious. Incidents of death, life-threatening illness, hospitalisation, or prolongation of existing hospitalisation, or permanent disability must be classified as serious, as defined by the Code of Federal Regulations (Myers et al., 2020). The standardisation of symptoms reported in VAERS is carried out by third-party experts according to the standard International Medical Dictionary of Terms (MedDRA) (Zhou et al., 2003), a dictionary of medical terms used for the reporting and monitoring of adverse drug reactions. The dictionary is organised into five levels, from broad to specific, and includes Lowest Level Terminology (LLT), Preferred Terminology (PT), Higher Level Terminology (HLT), Higher Level Group Terminology (HLGT), and 27 System Organ Class (SOC). This classification and standardisation have been shown to improve the efficiency of data analysis and signal detection (Bousquet et al., 2005). In this study, an initial set of thyroid-related Preferred Terms (PTs) was derived from relevant narrow scope Structured MedDRA Queries (SMQs). This set was subsequently refined by two independent clinical endocrinologists to ensure clinical relevance and alignment with the study objective.

The search was conducted on the database provided by VAERS for data with an initial cut-off date of 31 December 2024. The raw data files, including VAERSDATA, VAERSVAX, and VAERSSYMPTOMS, were downloaded in CSV format and subsequently collated and merged. The report was subject to screening in accordance with two conditions: firstly, the presence of HPV vaccination (HPV2, HPV4, HPV9, HPVX). Secondly, the presence of thyroid-related AE was taken into consideration. The selection statistics collated baseline information (state/region, age, gender, time from vaccination to adverse event, serious vs. non-serious, etc.), vaccine information (type, dose, number of times vaccinated, etc.), clinical outcomes, and annual uptake frequency statistics; and time to induction was statistically grouped.

Multiple pharmacovigilance signal detection algorithms were applied in order to validate potential positive signals for thyroid-associated adverse events following HPV vaccination. These included the proportional reporting ratio (PRR), the reporting reporting odds ratio (ROR), the Bayesian confidence propagation neural network (BCPNN), and the Multi-item Gamma Poisson Shrinker (MGPS) algorithm. The underlying principles of these methodologies are predicated upon a 2x2 columnar table (wherein ‘a' denotes target ADE occurrence, ‘b' denotes target ADE non-occurrence, ‘c' denotes target vaccination, and ‘d' denotes non-vaccination or vaccination with other vaccines). Supplementary Table S1 provides detailed formulas for this analytical method and the thresholds for positive signals. According to Toshiyuki Sakaeda et al., when at least one of the four indicators met the aforementioned criteria, the adverse event was classified as drug-related (Sakaeda et al., 2013) and could be further reviewed. All analyses were performed using R (version 4.4.3) in RStudio.

3 Results

3.1 General report characteristics

Figure 1 illustrates the data processing and analysis flow. As of 31 December 2024, the VAERS database contained 2,647,548 case reports and 11,222,196 individual adverse events. Among these, 78,641 reports were identified as HPV vaccine-related. After deduplication, 60,840 unique cases involving individuals who received only HPV vaccines were retained for analysis. The distribution of vaccine types was as follows: HPV2 (7.7%), HPV4 (61.3%), HPV9 (27.3%), and unspecified type (HPVX, 3.7%).

Figure 1

Figure 1. Flowchart for Selecting HPV Vaccine-Related Adverse Events from the VAERS Database. This flowchart outlines the process for identifying adverse events (AEs) related to the HPV vaccine using data from the Vaccine Adverse Event Reporting System (VAERS). It illustrates the total number of reports, the removal of duplicate records, and the methodologies employed for signal detection, including Reporting Odds Ratio (ROR), Proportional Reporting Ratio (PRR), Bayesian Confidence Propagation Neural Network (BCPNN), and Multi-item Gamma Poisson Shrinker (MGPS). Finally, it details the PT-level screening focused on identifying thyroid-related AEs through various analytical approaches.

Demographic analysis revealed a predominance of female recipients (69.4%), with the most frequent age group being under 18 years (40.5%), followed by those aged 18–65 years (22.8%). A total of 10,349 reports (17.0%) were classified as serious, 65.1% of which were associated with HPV4, and 90.6% occurred in females. Among all HPV-related reports, 579 deaths were recorded, accounting for 1.0% of the total (Table 1). Furthermore, there were 21,844 adverse event reports indicating that they occurred within 30 days after vaccination, which suggests a temporal clustering phenomenon.

Table 1

Table 1. Clinical characteristics of reports with HPV vaccine from the VAERS database.

3.2 Signal detection of thyroid-related adverse events

To identify potential associations between HPV vaccination and thyroid-related AEs, four validated signal detection algorithms were applied: ROR, PRR, BCPNN, and MGPS.

Analysis was conducted across 27 SOCs, as defined by MedDRA terminology (Figure 2). Among these, 15 SOCs demonstrated statistically significant positive signals. Notably, four SOCs were identified as having strong signals—detected by at least three of the four algorithms—namely: (1) Pregnancy, puerperium and perinatal conditions; (2) Benign, malignant and unspecified neoplasms (including cysts and polyps); (3) Social circumstances; and (4) Endocrine system disorders.

Figure 2

Figure 2. Adverse Event Reports by SOC Following HPV Vaccination. The Fig classifies statistics of adverse events (aes) after HPV vaccination according to SOC. The ordinate lists SOC categories (e.g., neurological disorders, gastrointestinal disorders) and corresponding AE terms. The values on the abscissa represent the number of cases.

At the PT level, 13 thyroid-related AEs exhibited statistically significant positive signals (Supplementary Table S2; Figures 3, 4). The most prominent signals were observed for autoimmune thyroiditis (ROR = 4.26, 95% CI: 3.48–5.23; EBGM05 = 3.32) and autoimmune hypothyroidism (ROR = 11.65, 95% CI: 5.00–27.15; EBGM05 = 4.52), with 103 and 82 reported cases, respectively. Additional positive signals included: thyroid cyst (ROR = 5.89), general thyroid disease (ROR = 2.76), hypothyroidism (ROR = 2.46), goiter (ROR = 1.93), Graves’ disease (ROR = 1.47), Hashimoto’s thyrotoxicosis (ROR = 38.28), hypothyroid goiter (ROR = 38.28), chronic thyroiditis (ROR = 12.76), and euthyroid sick syndrome (ROR = 6.38). Importantly, significant signals were also detected for thyroid malignancies, including thyroid cancer (ROR = 7.22) and papillary thyroid carcinoma (ROR = 4.25). Among these, autoimmune thyroiditis, autoimmune hypothyroidism, and thyroid cancer are categorized as important medical events (IMEs).

Figure 3

Figure 3. Volcano plot of thyroid-related adverse events post-HPV vaccination. This volcano plot illustrates the association between thyroid-related adverse events and HPV vaccination, with the log of ROR plotted against the negative log of adjusted p-values. Each point represents a specific thyroid condition, with color gradients from blue to red indicating the logarithm of the number of cases reported (log (Number of Cases)). A dashed vertical line denotes an ROR of 1, while a horizontal dashed line indicates a false discovery rate (FDR) adjusted p-value threshold of 0.05.

Figure 4

Figure 4. Forest Plot of Thyroid-Related Adverse Events Following HPV Vaccination. This forest plot shows the relative odds ratios (RORs) and 95% confidence intervals (CIs) for various thyroid-related adverse events reported after HPV vaccination. Each square represents the ROR for a particular condition, and the horizontal black line represents the corresponding 95% ci. The dashed line at ROR 1 represents the null hypothesis, indicating no association. Data suggest that risk estimates for some thyroid conditions increase after vaccination.

Subgroup analyses by sex and age (Tables 2, 3) revealed distinct patterns. The vast majority of thyroid-related AEs were reported in female recipients, including autoimmune thyroiditis, autoimmune hypothyroidism, chronic thyroiditis, Hashimoto’s thyrotoxicosis, thyroid cysts, goiter, thyroid cancer, and papillary thyroid carcinoma. Autoimmune thyroiditis was the most frequently reported AE (n = 100). Among male recipients, positive signals were limited to autoimmune hypothyroidism.

Table 2

Table 2. Gender-stratified detection of thyroid adverse event signals linked to HPV vaccination.

Table 3

Table 3. Age-stratified detection of thyroid adverse event signals linked to HPV vaccination.

No thyroid-related AEs were reported in individuals over the age of 65. All relevant reports originated from patients aged 65 years or younger, with a notable concentration among adolescents under 18 years of age. In this group, commonly reported thyroid-related AEs included autoimmune hypothyroidism (n = 66), thyroid cysts, autoimmune thyroiditis, hypothyroidism, goiter, thyroid tumors, and hyperthyroidism. Autoimmune thyroiditis remained the most frequently reported AE (n = 74). Several previously unreported AEs, such as goiter, hyperthyroidism (n = 29), and thyroid tumors, were also observed. These findings suggest a potential association between HPV vaccination and thyroid-related AEs, particularly in women and younger populations.

In addition, subgroup analysis, stratified by HPV vaccine types, revealed that HPV-4 was associated with the highest absolute number of thyroid-related adverse event (AE) reports (n = 187, 73.0% of total thyroid-related AEs), followed by HPV-9 (n = 54, 21.1%) and HPV-2 (n = 15, 5.9%) (Supplementary Table S3). Autoimmune thyroiditis was the predominant AE subtype across all vaccine groups, accounting for 66.7% (10/15) of AEs in HPV-2 recipients, 35.8% (67/187) in HPV-4 recipients, and 40.7% (22/54) in HPV-9 recipients. Hypothyroidism exhibited a higher reporting proportion in the HPV-4 cohort (32.6%, 61/187), with no comparable signal observed in the HPV-2 or HPV-9 subgroups.

3.3 Time-to-onset (TTO) analysis

Reports of thyroid-related AEs following HPV vaccination were primarily concentrated within the first 30 days post-vaccination, accounting for 45.08% of cases. Although the frequency declined over time, a notable proportion (16.06%) of reports still occurred beyond 360 days (Figure 5). This pattern reflects a bimodal distribution, characterized by a peak in early-onset reporting and a secondary rise in the long-term post-vaccination period.

Figure 5

Figure 5. Temporal distribution of Thyroid Adverse Events Following HPV Vaccination. (A) Survival curve demonstrating cumulative incidence of thyroid AEs during 360-day follow-up. Vertical axis: Cumulative incidence (%), horizontal axis: Days post-vaccination. (B) Butterfly plot showing cumulative incidence of thyroid adverse events (AEs) stratified by post-vaccination intervals. Left panel (blue): Cumulative incidence (%) across sequential time windows Right panel (orange): Absolute case counts per interval.

After excluding reports with missing, inaccurate, or unknown onset dates, 68 cases of autoimmune thyroiditis and 4 cases of thyroid cancer were included in the TTO analysis. For autoimmune thyroiditis, the median time to onset was 44 days (interquartile range: 12–248 days). The cumulative reporting rate reached 35.3% within 30 days, increased to 69.1% within 180 days, and remained at 20.6% 1 year after vaccination. Due to the small number of thyroid cancer cases (n = 4), statistical comparison was not performed; the median onset time for these cases was 185.5 days.

The Weibull shape parameter (β) for thyroid-related AEs, calculated using the WSP model, was 0.54 (95% CI: 0.485–0.596) (Table 4), indicating that the reporting risk was highest during the early post-vaccination period and declined thereafter.

Table 4

Table 4. Time to onset of HPV vaccine-associated adverse thyroid reaction signals was analysed using the Weibull distribution test.

4 Discussion

Although HPV vaccination plays a critical role in cervical cancer prevention, concerns remain regarding its potential association with thyroid disorders. This study systematically analyzed thyroid-related AEs following HPV vaccination using the VAERS database. A total of 13 thyroid-related AEs demonstrated statistically significant positive signals, with autoimmune thyroid disease (AITD) and hypothyroidism being the most frequently reported.

The potential association between HPV vaccination and AITD remains a subject of ongoing debate. In our study, positive signals for AITD onset were observed in the context of HPV vaccination. This finding echoes a large cohort study that showed HPV vaccination was statistically significantly associated with an increased risk of autoimmune diseases such as Hashimoto’s thyroiditis (HT) (Hviid et al., 2018). A meta-analysis further supported this association, indicating a 1.5-fold increased risk of AITD after vaccination (Rosillon et al., 2020). However, other studies have reported no significant elevation in thyroiditis risk within 2 years post-vaccination (Castillo-Cano et al., 2022), and a Korean cohort study found no increased risk of HT in either cohort or self-controlled case series analyses (adjusted rate ratio = 1.24; 95% CI: 0.78–1.90; adjusted relative risk = 1.28; 95% CI: 0.53–3.08) (Yoon et al., 2021). Interestingly, Grimaldi-Bensouda et al. observed a decreased risk of AITD in vaccinated girls compared to unvaccinated individuals over a 24-month period (OR = 0.35; 95% CI: 0.13–0.92) (Grimaldi-Bensouda et al., 2017). Our study highlighted a notable signal in the safety data: reports of hypothyroidism and autoimmune thyroiditis following HPV vaccination that we believe merit closer scrutiny and further investigation. Notably, three preferred terms—autoimmune thyroiditis, autoimmune hypothyroidism, and thyroid carcinoma—were classified as IMEs. These findings suggest that individualized vaccination strategies, such as avoiding periods of hormonal fluctuation (e.g., adolescence), or implementing baseline screening and post-vaccination monitoring of thyroid function (e.g., TSH, anti-TPO antibodies), especially in individuals presenting with fatigue or unexplained weight changes, may warrant further evaluation to mitigate potential risks.

The underlying pathogenesis may involve adjuvant-mediated immune activation or molecular mimicry. Immunostimulatory adjuvants used in some vaccines have been implicated in triggering autoimmune phenomena (Perricone et al., 2016; Cohen Tervaert et al., 2023). Cervarix^® employs the AS04 adjuvant system. Studies demonstrate enhanced immunogenicity of AS04-adjuvanted vaccines compared to aluminum-only formulation (Garçon et al., 2011b). AS04 potentiates CD4⁺ T-cell proliferation and functional activation, thereby promoting memory B-cell differentiation (Garçon et al., 2011a). Notably, CD4⁺ T-cell hyperactivation constitutes a key pathogenic mechanism in Hashimoto’s thyroiditis development (Pyzik et al., 2015). Another hypothesis involves molecular mimicry and cross-reactivity between HPV antigens and thyroid autoantigens. Clinical observations have reported the presence of anti-thyroid peroxidase (TPO) antibodies in some vaccinated individuals (Colafrancesco et al., 2013), and experimental models have demonstrated that administration of thyroglobulin (Tg) or TPO can induce autoimmune thyroiditis in mice, suggesting these autoantigens play a critical role in HT pathogenesis (Pyzik et al., 2015).

Subgroup analysis stratified by HPV vaccine type revealed that 4vHPV exhibited the highest absolute number of thyroid-related adverse events reports (n = 187, accounting for 72.8% of total thyroid-related AEs), followed by 9vHPV (n = 54, 21.0%) and 2vHPV (n = 16, 6.2%). The distinct immunogenic and reactogenic profiles observed among the three HPV vaccines may arise from variations in their manufacturing platforms and adjuvant compositions (Leung et al., 2015). 2vHPV utilizes a baculovirus-mediated expression system in *Trichoplusia ni* Rix4446 cells to generate virus-like particles, which are combined with the AS04 adjuvant system (Garçon et al., 2011a). Mechanistic studies indicate that AS04 potentiates both humoral and cellular immunity through transient cytokine induction at the injection site (Giannini et al., 2006). This process enhances antigen-presenting cell activation and facilitates efficient antigen presentation to CD4⁺ T lymphocytes (Garçon et al., 2011a). In contrast,4vHPV and nonavalent 9vHPV vaccines are produced in the yeast Saccharomyces cerevisiae, coupled with amorphous aluminum hydroxyphosphate sulfate (AAHS) as an adjuvant (Garçon et al., 2007). AAHS exhibits superior binding affinity to L1 VLPs relative to conventional aluminum salts (Caulfield et al., 2007). However, there has not been any direct comparison between AS04 and AAHS adjuvants using identical HPV antigens. Autoimmune thyroiditis emerged as the predominant AE subtype across all vaccine groups, with potential pathogenic mechanisms possibly involving adjuvant-mediated immune activation or molecular mimicry, consistent with our overall study. Furthermore, hypothyroidism demonstrated a higher reporting proportion in the 4vHPV cohort (32.6%, 61/187), while no similar signal was observed in 2vHPV or 9vHPV subgroups. This potential association between hypothyroidism and 4vHPV vaccination may represent a novel finding. However, considering the study’s limitations and inherent complexity, we recommend conducting additional targeted investigations to validate and expand these observations, thereby enabling a comprehensive assessment of HPV vaccination’s impact on thyroid function.

Marked gender disparities were observed in this study. The vast majority of thyroid-related AEs occurred in females. Among male recipients, only autoimmune thyroid disorders exhibited stronger signals. Gender-specific vaccination rates may potentially confound the outcomes. This sex bias may reflect known differences in immune responses: females typically mount stronger humoral and cellular immunity, generate higher antibody titers to various antigens (especially T-cell independent), and display more rapid allograft rejection (Chiovato et al., 1993). Estrogen also modulates immune responses and is associated with increased susceptibility to anxiety, depression, and eating disorders—all of which have been linked to thyroid dysfunction and may contribute to increased AE reporting in women (Plaisier et al., 2008). Furthermore, AITD affects women up to four times more frequently than men, which may partially account for the observed imbalance (Hu et al., 2022). In summary, the immune and endocrine mechanistic differences discussed above provide an inherent biological explanation for the observed sex disparities, which constitutes the primary finding of this study. It is noteworthy that confounding factors such as vaccination rates and reporting behaviors may introduce additional potential bias, which could also exert certain influences on the results and still requires further investigation through more rigorously designed prospective studies in the future.

Furthermore, age-stratified analyses revealed heterogeneity in AE distribution, with thyroid-related events predominantly observed in younger vaccinees and no thyroid-related AEs were found in older vaccinees. This may be associated with age-related immunosenescence in older adults, a condition characterized by an impaired capacity to mount effective protective humoral and cellular immune responses against pathogens or vaccines (Crooke et al., 2019). People who can be vaccinated according to the guidelines are individuals aged 9–45 years (Zadeh Mehrizi et al., 2025). While high vaccination rates among the younger population may confound the results, our methodology employing signal intensity calculations effectively mitigates the impact of population size disparities. This phenomenon aligns with AITD epidemiology. AITD is most common between the ages of 30 and 60 (Wrońska et al., 2024). Additionally, depression and anxiety disorders are more likely to occur in middle age and childbearing years (Plaisier et al., 2008), while environmental factors such as stress play significant roles in AITD pathogenesis (Tomer and Huber, 2009; Ajjan and Weetman, 2015). The absence of significant thyroid AE signals in elderly subgroups remains incompletely elucidated pathophysiologically, warranting further investigation.

The time-to-onset analysis revealed a distinct bimodal distribution in thyroid-related adverse events post HPV vaccination. Early-phase risks were predominant, with 45.08% (95% CI: 39.2%–51.1%) of cases emerging within 30 days post-vaccination. A secondary peak comprising 16.06% (95% CI: 12.3%–20.5%) of cases manifested beyond 365 days. Our study identified that thyroid-related adverse events significantly associated with HPV vaccination were predominantly comprised of autoimmune diseases, including autoimmune thyroiditis (ROR = 4.26), autoimmune hypothyroidism (ROR = 11.65), and Hashimoto’s disease (ROR = 38.28), among others. This disease spectrum provides critical insights for interpreting the bimodal temporal distribution—particularly the second delayed peak. We hypothesize that the second peak may be primarily driven by biological mechanisms, whereby vaccination triggers a slowly progressive autoimmune process (such as Hashimoto’s thyroiditis) in genetically or immunologically susceptible individuals (Cohen Tervaert et al., 2023). Hashimoto’s thyroiditis typically exhibit an insidious onset, characterized by chronic and progressive inflammation and tissue damage (Ralli et al., 2020). The progression from immune initiation to the manifestation of typical clinical symptoms often takes months to years, which is highly consistent with the delayed peak observed in this study. On the other hand, this peak may also be partially attributable to reporting bias or diagnostic delay, wherein non-specific thyroid-related symptoms (e.g., fatigue, weight changes) are not promptly recognized and associated with vaccination, or are reported to pharmacovigilance systems long after the adverse event occurrence. The Weibull distribution analysis (shape parameter β = 0.54, 95% CI: 0.485-0.596; β < 1) indicates a monotonically decreasing hazard function over time, consistent with progressive risk reduction. These results suggested that we should be vigilant about the AEs associated with HPV vaccination in the first month, and early recognition of AEs caused by HPV vaccination is important, while also promptly addressing AEs that may arise. These findings need to be confirmed through prospective studies and long-term follow-up.

The strengths of our study are the ability of VAERS to accept reports from a variety of sources (healthcare providers, patients, and others) nationwide, in addition to rapid signal recognition and efficacy in monitoring rare adverse events.

5 Limitations

Notwithstanding, limitations exist. VAERS operates as a passive surveillance system collecting post-marketing adverse event reports from pharmaceutical manufacturers, physicians, regulatory agencies, and global vaccine recipients following FDA-approved administrations. Underreporting and delayed reporting persist despite mandatory AE reporting requirements. Submitted reports occasionally contain incomplete data and lack standardized reporting formats. Collectively, these constraints preclude robust causal inference. Nevertheless, VAERS remains invaluable for hypothesis generation, with subsequent validation recommended through more rigorous databases like the Vaccine Safety Datalink (VSD), which enables statistically stringent causality assessments.

6 Conclusion

This study identified disproportionality signals suggesting a potential association between HPV vaccination and various thyroid-related adverse events, particularly autoimmune thyroiditis and hypothyroidism. The majority of events were reported in females and younger individuals, with most occurring within the first 6 months post-vaccination. Subgroup and time-to-onset analyses further supported temporal clustering and demographic variability in AE reporting patterns. Although the underlying mechanisms remain to be fully elucidated, immune activation and molecular mimicry may contribute to the observed associations.

Given the limitations inherent in passive surveillance systems, these findings do not establish causality. However, they highlight the importance of continued pharmacovigilance and suggest that targeted monitoring of thyroid function in susceptible individuals may be warranted. Further prospective studies and mechanistic investigations are needed to clarify the potential immunological links between HPV vaccination and thyroid autoimmunity.

Data availability statement

All data analyzed in this study were obtained from the publicly available Vaccine Adverse Event Reporting System (VAERS; https://vaers.hhs.gov/). Derived data and analysis code are available from the corresponding author on reasonable request.

Ethics statement

This study used de-identified and publicly available data from the Vaccine Adverse Event Reporting System (VAERS). In accordance with institutional and national policies, ethics approval and informed consent were not required.

Author contributions

LS: Conceptualization, Data curation, Formal Analysis, Supervision, Visualization, Writing – original draft, Writing – review and editing. JD: Data curation, Formal Analysis, Methodology, Software, Supervision, Visualization, Writing – original draft, Writing – review and editing. JF: Data curation, Investigation, Project administration, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing. RX: Data curation, Formal Analysis, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing. SL: Data curation, Formal Analysis, Software, Validation, Writing – original draft, Writing – review and editing. GD: Investigation, Methodology, Project administration, Validation, Writing – original draft, Writing – review and editing. YD: Formal Analysis, Validation, Visualization, Writing – original draft, Writing – review and editing. SH: Conceptualization, Data curation, Validation, Writing – original draft, Writing – review and editing. JL: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Shanxi Provincial Key Research and Development Project (Grant No. 2022048), The First Hospital of Shanxi Medical University.

Acknowledgments

Acknowledgements

This study utilized data from the FDA Adverse Event Reporting System (FAERS), made publicly available by the U.S. Food and Drug Administration. The analyses, interpretations, and conclusions presented in this study are solely those of the authors and do not reflect the views or official positions of the FDA.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2025.1664697/full#supplementary-material

References

Ajjan, R. A., and Weetman, A. P. (2015). The pathogenesis of hashimoto's thyroiditis: further developments in our understanding. Horm. Metab. Res. 47 (10), 702–710. doi:10.1055/s-0035-1548832

PubMed Abstract | CrossRef Full Text | Google Scholar

Bousquet, C., Lagier, G., Lillo-Le Louët, A., Le Beller, C., Venot, A., and Jaulent, M. C. (2005). Appraisal of the MedDRA conceptual structure for describing and grouping adverse drug reactions. Drug Saf. 28 (1), 19–34. doi:10.2165/00002018-200528010-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Castillo-Cano, B., Martín-Pérez, M., Llorente-García, A., Montero-Corominas, D., Comas-Cufí, M., and Martín-Merino, E. (2022). Assessment of thyroiditis risk associated with HPV vaccination among girls aged 9-18 years: a time-varying cohort study. Vaccine 40 (33), 4816–4826. doi:10.1016/j.vaccine.2022.06.060

PubMed Abstract | CrossRef Full Text | Google Scholar

Caulfield, M. J., Shi, L., Wang, S., Wang, B., Tobery, T. W., Mach, H., et al. (2007). Effect of alternative aluminum adjuvants on the absorption and immunogenicity of HPV16 L1 VLPs in mice. Hum. Vaccin. 3 (4), 139–145. doi:10.4161/hv.3.4.4309

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Chen, C., Cao, L., Fang, J., and Xiao, J. (2024). Comparative safety profile of bivalent and original COVID-19 mRNA vaccines regarding myocarditis/pericarditis: a pharmacovigilance study. Int. Immunopharmacol. 133, 112022. doi:10.1016/j.intimp.2024.112022

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiovato, L., Lapi, P., Fiore, E., Tonacchera, M., and Pinchera, A. (1993). Thyroid autoimmunity and female gender. J. Endocrinol. Invest. 16 (5), 384–391. doi:10.1007/bf03348863

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen Tervaert, J. W., Martinez-Lavin, M., Jara, L. J., Halpert, G., Watad, A., Amital, H., et al. (2023). Autoimmune/Inflammatory syndrome induced by adjuvants (ASIA) in 2023. Autoimmun. Rev. 22 (5), 103287. doi:10.1016/j.autrev.2023.103287

PubMed Abstract | CrossRef Full Text | Google Scholar

Colafrancesco, S., Perricone, C., Tomljenovic, L., and Shoenfeld, Y. (2013). Human papilloma virus vaccine and primary ovarian failure: another facet of the autoimmune/inflammatory syndrome induced by adjuvants. Am. J. Reprod. Immunol. 70 (4), 309–316. doi:10.1111/aji.12151

PubMed Abstract | CrossRef Full Text | Google Scholar

Crooke, S. N., Ovsyannikova, I. G., Poland, G. A., and Kennedy, R. B. (2019). Immunosenescence: a systems-level overview of immune cell biology and strategies for improving vaccine responses. Exp. Gerontol. 124, 110632. doi:10.1016/j.exger.2019.110632

PubMed Abstract | CrossRef Full Text | Google Scholar

Garçon, N., Chomez, P., and Van Mechelen, M. (2007). GlaxoSmithKline adjuvant systems in vaccines: concepts, achievements and perspectives. Expert Rev. Vaccines 6 (5), 723–739. doi:10.1586/14760584.6.5.723

PubMed Abstract | CrossRef Full Text | Google Scholar

Garçon, N., Morel, S., Didierlaurent, A., Descamps, D., Wettendorff, M., and Van Mechelen, M. (2011a). Development of an AS04-adjuvanted HPV vaccine with the adjuvant system approach. BioDrugs 25 (4), 217–226. doi:10.2165/11591760-000000000-00000

PubMed Abstract | CrossRef Full Text | Google Scholar

Garçon, N., Segal, L., Tavares, F., and Van Mechelen, M. (2011b). The safety evaluation of adjuvants during vaccine development: the AS04 experience. Vaccine 29 (27), 4453–4459. doi:10.1016/j.vaccine.2011.04.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Giannini, S. L., Hanon, E., Moris, P., Van Mechelen, M., Morel, S., Dessy, F., et al. (2006). Enhanced humoral and memory B cellular immunity using HPV16/18 L1 VLP vaccine formulated with the MPL/aluminium salt combination (AS04) compared to aluminium salt only. Vaccine 24 (33-34), 5937–5949. doi:10.1016/j.vaccine.2006.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Grimaldi-Bensouda, L., Rossignol, M., Koné-Paut, I., Krivitzky, A., Lebrun-Frenay, C., Clet, J., et al. (2017). Risk of autoimmune diseases and human papilloma virus (HPV) vaccines: six years of case-referent surveillance. J. Autoimmun. 79, 84–90. doi:10.1016/j.jaut.2017.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, X., Chen, Y., Shen, Y., Tian, R., Sheng, Y., and Que, H. (2022). Global prevalence and epidemiological trends of hashimoto's thyroiditis in adults: a systematic review and meta-analysis. Front. Public Health 10, 1020709. doi:10.3389/fpubh.2022.1020709

PubMed Abstract | CrossRef Full Text | Google Scholar

Hviid, A., Svanström, H., Scheller, N. M., Grönlund, O., Pasternak, B., and Arnheim-Dahlström, L. (2018). Human papillomavirus vaccination of adult women and risk of autoimmune and neurological diseases. J. Intern. Med. 283 (2), 154–165. doi:10.1111/joim.12694

PubMed Abstract | CrossRef Full Text | Google Scholar

Leung, T. F., Liu, A. P., Lim, F. S., Thollot, F., Oh, H. M., Lee, B. W., et al. (2015). Comparative immunogenicity and safety of human papillomavirus (HPV)-16/18 AS04-adjuvanted vaccine and HPV-6/11/16/18 vaccine administered according to 2- and 3-dose schedules in girls aged 9-14 years: results to month 12 from a randomized trial. Hum. Vaccin. Immunother. 11 (7), 1689–1702. doi:10.1080/21645515.2015.1050570

PubMed Abstract | CrossRef Full Text | Google Scholar

Moro, P. L., Zhang, B., Ennulat, C., Harris, M., McVey, R., Woody, G., et al. (2023). Safety of co-administration of mRNA COVID-19 and seasonal inactivated influenza vaccines in the vaccine adverse event reporting system (VAERS) during July 1, 2021-June 30, 2022. Vaccine 41 (11), 1859–1863. doi:10.1016/j.vaccine.2022.12.069

PubMed Abstract | CrossRef Full Text | Google Scholar

Myers, T. R., McNeil, M. M., Ng, C. S., Li, R., Marquez, P. L., Moro, P. L., et al. (2020). Adverse events following quadrivalent meningococcal diphtheria toxoid conjugate vaccine (menactra®) reported to the vaccine adverse event reporting system (VAERS), 2005-2016. Vaccine 38 (40), 6291–6298. doi:10.1016/j.vaccine.2020.07.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Perricone, C., Colafrancesco, S., Mazor, R. D., Soriano, A., Agmon-Levin, N., and Shoenfeld, Y. (2016). Corrigendum to Autoimmune/inflammatory syndrome induced by adjuvants (ASIA) 2013: unveiling the pathogenic, clinical and diagnostic aspects. J. Autoimmun. 72, 126. doi:10.1016/j.jaut.2016.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Plaisier, I., de Bruijn, J. G., Smit, J. H., de Graaf, R., Ten Have, M., Beekman, A. T., et al. (2008). Work and family roles and the association with depressive and anxiety disorders: differences between men and women. J. Affect Disord. 105 (1-3), 63–72. doi:10.1016/j.jad.2007.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Pyzik, A., Grywalska, E., Matyjaszek-Matuszek, B., and Roliński, J. (2015). Immune disorders in hashimoto's thyroiditis: what do we know so far? J. Immunol. Res. 2015, 979167. doi:10.1155/2015/979167

PubMed Abstract | CrossRef Full Text | Google Scholar

Ralli, M., Angeletti, D., Fiore, M., D'Aguanno, V., Lambiase, A., Artico, M., et al. (2020). Hashimoto's thyroiditis: an update on pathogenic mechanisms, diagnostic protocols, therapeutic strategies, and potential malignant transformation. Autoimmun. Rev. 19 (10), 102649. doi:10.1016/j.autrev.2020.102649

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosillon, D., Willame, C., Tavares Da Silva, F., Guignard, A., Caterina, S., Welby, S., et al. (2020). Meta-analysis of the risk of autoimmune thyroiditis, guillain-barré syndrome, and inflammatory bowel disease following vaccination with AS04-adjuvanted human papillomavirus 16/18 vaccine. Pharmacoepidemiol. Drug Saf. 29 (9), 1159–1167. doi:10.1002/pds.5063

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakaeda, T., Tamon, A., Kadoyama, K., and Okuno, Y. (2013). Data mining of the public version of the FDA adverse event reporting system. Int. J. Med. Sci. 10 (7), 796–803. doi:10.7150/ijms.6048

PubMed Abstract | CrossRef Full Text | Google Scholar

Sever, J. L., Brenner, A. I., Gale, A. D., Lyle, J. M., Moulton, L. H., Ward, B. J., et al. (2004). Safety of anthrax vaccine: an expanded review and evaluation of adverse events reported to the vaccine adverse event reporting system (VAERS). Pharmacoepidemiol. Drug Saf. 13 (12), 825–840. doi:10.1002/pds.936

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, H., Parekh, S., Pujari, P., Shewale, S., Desai, S., Bhatla, N., et al. (2023). Immunogenicity and safety of a new quadrivalent HPV vaccine in girls and boys aged 9-14 years versus an established quadrivalent HPV vaccine in women aged 15-26 years in India: a randomised, active-controlled, multicentre, phase 2/3 trial. Lancet Oncol. 24 (12), 1321–1333. doi:10.1016/s1470-2045(23)00480-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimabukuro, T. T., Nguyen, M., Martin, D., and DeStefano, F. (2015). Safety monitoring in the vaccine adverse event reporting system (VAERS). Vaccine 33 (36), 4398–4405. doi:10.1016/j.vaccine.2015.07.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Steinbach, A., and Riemer, A. B. (2018). Immune evasion mechanisms of human papillomavirus: an update. Int. J. Cancer 142 (2), 224–229. doi:10.1002/ijc.31027

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomer, Y., and Huber, A. (2009). The etiology of autoimmune thyroid disease: a story of genes and environment. J. Autoimmun. 32 (3-4), 231–239. doi:10.1016/j.jaut.2009.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai, S. A., Lu, C. Y., Chen, T. I., Huang, S. P., and Chen, Y. C. (2023). Adverse events from HPV vaccination in Taiwan. Vaccine 41 (49), 7444–7449. doi:10.1016/j.vaccine.2023.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Willame, C., Gadroen, K., Bramer, W., Weibel, D., and Sturkenboom, M. (2020). Systematic review and meta-analysis of postlicensure observational studies on human papillomavirus vaccination and autoimmune and other rare adverse events. Pediatr. Infect. Dis. J. 39 (4), 287–293. doi:10.1097/inf.0000000000002569

PubMed Abstract | CrossRef Full Text | Google Scholar

Wrońska, K., Hałasa, M., and Szczuko, M. (2024). The role of the immune system in the course of hashimoto's thyroiditis: the current state of knowledge. Int. J. Mol. Sci. 25 (13), 6883. doi:10.3390/ijms25136883

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, D., Lee, J. H., Lee, H., and Shin, J. Y. (2021). Association between human papillomavirus vaccination and serious adverse events in South Korean adolescent girls: nationwide cohort study. Bmj 372, m4931. doi:10.1136/bmj.m4931

PubMed Abstract | CrossRef Full Text | Google Scholar

Yousaf, S., Shehzadi, A., Ahmad, M., Asrar, A., Ahmed, I., Iqbal, H. M. N., et al. (2024). Recent advances in HPV biotechnology: understanding host-virus interactions and cancer progression - a review. Int. J. Surg. 110 (12), 8025–8036. doi:10.1097/js9.0000000000002117

PubMed Abstract | CrossRef Full Text | Google Scholar

Yousefi, Z., Aria, H., Ghaedrahmati, F., Bakhtiari, T., Azizi, M., Bastan, R., et al. (2021). An update on human papilloma virus vaccines: history, types, protection, and efficacy. Front. Immunol. 12, 805695. doi:10.3389/fimmu.2021.805695

PubMed Abstract | CrossRef Full Text | Google Scholar

Zadeh Mehrizi, T., Ataei-Pirkooh, A., Eshrati, B., and Ebrahimi Shahmabadi, H. (2025). Investigating factors affecting the effectiveness of gardasil 4, cervarix, and gardasil 9 vaccines considering the WHO regions in females: a systematic review. Cancer Epidemiol. 95, 102759. doi:10.1016/j.canep.2025.102759

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, W., Pool, V., Iskander, J. K., English-Bullard, R., Ball, R., Wise, R. P., et al. (2003). Surveillance for safety after immunization: Vaccine adverse event reporting system (VAERS)--United states, 1991-2001. MMWR Surveill. Summ. 52 (1), 1–24.

PubMed Abstract | Google Scholar