Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Immunol., 25 August 2025

Sec. Vaccines and Molecular Therapeutics

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1621939

Risk of developing juvenile idiopathic arthritis after quadrivalent HPV vaccination: a retrospective cohort study using the TriNetX U.S. Network

Wen-Yu WuWen-Yu Wu1Hsin-Hua Chen,,,,&#x;Hsin-Hua Chen2,3,4,5,6†Ming-Chin Tsai&#x;Ming-Chin Tsai1†Yung-Chieh Huang,,&#x;Yung-Chieh Huang1,4,7†Jun-Peng Chen&#x;Jun-Peng Chen8†Lin-Shien Fu,,*Lin-Shien Fu1,4,7*
  • 1Department of Pediatrics, Taichung Veterans General Hospital, Taichung, Taiwan
  • 2Division of Allergy, Immunology and Rheumatology, Department of Medicine, Taichung Veterans General Hospital, Taichung, Taiwan
  • 3Division of Clinical Informatics, Center for Quality Management, Taichung Veterans General Hospital, Taichung, Taiwan
  • 4Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
  • 5Program in Translational Medicine and Rong Hsing Research Center for Translational Medicine, National Chung Hsing University, Taichung, Taiwan
  • 6Big Data Center, National Chung Hsing University, Taichung, Taiwan
  • 7Department of Pediatrics, National Yang Ming Chiao Tung University, Taipei, Taiwan
  • 8Department of Research, Taichung Veterans General Hospital, Taichung, Taiwan

Introduction: Human papillomavirus (HPV) infection has been implicated in autoimmune processes, yet concerns remain about the potential autoimmune risks of HPV vaccination. Juvenile idiopathic arthritis (JIA) is a chronic autoimmune condition that typically manifests in childhood. The relationship between HPV vaccination and the development of JIA remains uncertain.

Methods: We conducted a retrospective cohort study using data from the TriNetX U.S. Collaborative Network. Females aged 9–13 years were included. Three analyses were performed: (1) comparing HPV4-vaccinated vs. unvaccinated matched cohorts; (2) a stricter comparison excluding subjects with positive ANA; (3) comparing single vs. multiple HPV4 doses. Propensity score matching and Cox proportional hazards models were used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs).

Results: In Analysis 1 (n=55,257 pairs) and Analysis 2 (n=53,827 pairs), the HPV4-vaccinated groups showed significantly reduced rates of JIA from 12 to 36 months post-vaccination (HR range: 0.33–0.52). No difference in JIA risk was observed between single and multiple doses in Analysis 3 (n=20,822 pairs). Early-onset JIA (<6 months after HPV4 vaccine) showed inconsistent trends, with only limited protective signals.

Conclusions: Our findings suggest that HPV4 vaccination is not associated with an increased risk of JIA. On the contrary, vaccination may confer a long-term protective effect against new-onset JIA, observable from 6–12 months and lasting for at least 3 years. These findings support the safety and possible immunoregulatory benefit of HPV4 in adolescents.

1 Introduction

Human papillomaviruses (HPVs) are double-stranded DNA viruses associated with a wide spectrum of human diseases, including cervical cancer and autoimmune conditions (1, 2). Prophylactic HPV vaccination, particularly the quadrivalent HPV4 vaccine, has demonstrated robust efficacy in preventing HPV-related infections and malignancies (3, 4). Despite its favorable safety profile, concerns remain regarding its potential role in triggering autoimmune disorders such as juvenile idiopathic arthritis (JIA).

JIA is the most prevalent autoimmune rheumatic disease in children, characterized by persistent joint inflammation without a known etiology. It is believed to result from complex interactions between genetic susceptibility and environmental triggers, including viral infections (5). Molecular mimicry and immune dysregulation have been proposed as key mechanisms by which pathogens may initiate or exacerbate autoimmunity (6, 7).

Previous studies have yielded conflicting findings on the association between HPV4 vaccination and the development of autoimmune arthritis. A U.S.-based analysis of the Vaccine Adverse Event Reporting System (VAERS) reported an elevated odds ratio for rheumatoid arthritis following HPV4 vaccination (8). Similarly, a Colombian study employing inverse probability of treatment weighting also observed increased arthritis risk (9). However, both relied on self-reported adverse events, potentially introducing reporting bias. Conversely, a retrospective study using electronic medical records from Kaiser Permanente found no significant increase in JIA incidence within 180 days post-vaccination (10).

To address these inconsistencies, we conducted a large-scale retrospective cohort study using the TriNetX U.S. Collaborative Network. This federated platform offers access to de-identified electronic health records across multiple healthcare organizations, enabling robust longitudinal analyses (11). We aimed to evaluate whether HPV4 vaccination is associated with an increased or decreased risk of new-onset JIA in adolescent females. Three complementary analyses were performed to examine the effect of vaccination status, antinuclear antibody (ANA) positivity, and HPV4 dose number on subsequent JIA risk.

2 Materials and methods

2.1 Data sources

Data were obtained from the TriNetX U.S. Collaborative Network, comprising 57 healthcare organizations with integrated electronic health records (EHRs). The study population included females aged 9 to 13 years between January 1, 2010 and December 31, 2015 (for Analyses 1 and 2) or up to December 31, 2017 (for Analysis 3). Three analyses were conducted as follows:

● Analysis 1: Compared HPV4-vaccinated individuals (n = 55,257) to an unvaccinated control cohort (n = 1,323,995). Individuals with a prior diagnosis of JIA (ICD-10: M08.xx), or prescription of methotrexate or antirheumatic drugs (Supplementary Tables 1, 2), were excluded.

● Analysis 2: Built on the criteria in Analysis 1, with additional exclusion of individuals with positive ANA serology (ICD-10: R76.0 or R68.89). This analysis involved 53,827 vaccinated individuals and 1,399,261 unvaccinated controls.

● Analysis 3: Focused on HPV4-vaccinated individuals only, comparing those who received a single dose (n = 10,411) versus multiple doses (n = 10,411). The extended study period accounted for the 6–12 month dosing interval.

TriNetX’s integrated propensity score matching algorithm was used to match cohorts 1:1 based on age at index event (first vaccination or clinical visit), race (White, Black, Asian), ICD-10 diagnoses related to the musculoskeletal system (M00–M99), and prior NSAID use. Only individuals identified as White, Black or African American, or Asian were included in race-based matching, due to data completeness and sample size considerations.

2.2 Study design

As shown in Figures 1A, B, matched cohorts were followed for up to 36 months post-index event. Incident JIA was defined by ICD-10 code M08.xx, with or without concurrent use of disease-modifying antirheumatic drugs (DMARDs) or NSAIDs. Cases diagnosed within 7 days post-index event were excluded to avoid capturing pre-existing or unrelated conditions.

Figure 1
Flowchart depicting inclusion and exclusion criteria for a study on HPV4 vaccination in females aged nine to thirteen between 2010 and 2015. Two parallel paths are illustrated: one excludes those with previous HPV vaccination and specific drug codes, leading to a group matched by age, race, diagnosis, and medication with 55,257 individuals. The other path differentiates between one-dose (21,311 individuals) and more than one-dose groups (29,746 individuals), ending with matching criteria resulting in 20,822 individuals in both groups.

Figure 1. (A) Flowchart of cohort construction for Analyses 1 and 2. HPV4-vaccinated and unvaccinated individuals (n = 55,257 per group) were selected from the TriNetX U.S. Collaborative Network and matched 1:1 using propensity scores based on age at index event, race, prior musculoskeletal diagnoses, and NSAID use. Follow-up lasted up to 36 months post-index. (B) Flowchart of cohort construction for Analysis 3. HPV4-vaccinated individuals were stratified into single-dose and multiple-dose groups (n = 10,411 per group) using 1:1 propensity score matching with identical covariates. Follow-up continued for up to 36 months post-index.

The number of new-onset JIA cases was assessed at intervals of 42 days, 3, 6, 12, 18, 24, 30, and 36 months post-event. Medication records were reviewed to confirm JIA diagnoses based on prescription of NSAIDs, methotrexate, or other DMARDs.

2.3 Ethical approval

The study protocol was approved by the Institutional Review Board of Taichung Veterans General Hospital (Approval No. CE23480C). Given the use of de-identified data, informed consent was waived. All data accessed through the TriNetX platform were de-identified in compliance with the Health Insurance Portability and Accountability Act (HIPAA) Privacy Rule, and no patient-level identifiers were available to the investigators.

2.4 Statistical analysis

To control for potential confounding, propensity score matching (1:1) was applied to balance baseline covariates between cohorts, including age at index event, race, prior musculoskeletal or connective tissue diagnoses (ICD-10 M00–M99), and NSAID use. Time-to-event analysis was performed using the Kaplan–Meier method based on daily intervals. Censoring was applied at the last recorded encounter. Hazard ratios (HRs) with 95% confidence intervals (CIs), log-rank tests, and proportionality assessments were generated through the TriNetX survival analysis module.

Due to the platform’s design, individual-level Schoenfeld residuals were unavailable; therefore, the proportional hazards assumption was indirectly assessed using cumulative incidence curves derived from daily JIA-free probabilities over a 36-month follow-up period. To account for multiple comparisons across time points, p-values were adjusted using the False Discovery Rate (FDR) method proposed by Benjamini and Hochberg. All statistical analyses were conducted using the TriNetX platform and R software (version 3.5; survival package). A two-tailed p-value < 0.05 was considered statistically significant.

3 Results

3.1 Cohort characteristics

The demographic and baseline clinical characteristics of the matched cohorts for all three analyses are presented in Table 1. In Analysis 1, 55,257 participants were included in both the HPV4-vaccinated and unvaccinated cohorts, with a mean age of 11.7 ± 1.0 years. The racial distribution was 51.0–51.1% White, 24.9% Black or African American, and 2.7% Asian. Baseline prevalence of musculoskeletal/connective tissue disorders was 18.9%, and NSAID use was noted in 15.8% of participants.

Table 1
www.frontiersin.org

Table 1. Baseline characteristics of study subjects after matching.

In Analysis 2, the matched vaccinated and unvaccinated groups each consisted of 53,827 individuals with similar age and race distributions. ANA-positive individuals were excluded. In Analysis 3, 20,822 HPV4-vaccinated individuals were divided into single-dose and multiple-dose cohorts (10,411 each). Demographics and clinical variables were well balanced following propensity score matching.

3.2 Incidence of new-onset JIA

The cumulative incidence and odds ratios (ORs) of new-onset JIA over 36 months are summarized in Tables 2A–C. JIA was defined either solely by ICD-10 M08.xx diagnosis (criterion 1) or by M08.xx with concurrent DMARD or NSAID usage (criterion 2). To account for multiple comparisons across time points, the False Discovery Rate (FDR) method adjusted p-values were also shown in Table 2.

www.frontiersin.org

.

In Analysis 1, there was no significant difference in JIA incidence at 42 days post-vaccination (criterion 1: OR = 0.71, 95% CI = 0.32–1.60; criterion 2: OR = 20.99, 95% CI = 1.23–358.19). However, significantly lower JIA rates were observed in the HPV4 cohort at 3 and 6 months by criterion 1 (OR = 0.45 and 0.33, respectively), but not by criterion 2. From 12 to 36 months, both criteria consistently showed a significantly reduced JIA incidence in the vaccinated cohort (p < 0.05 for all time points).

Analysis 2 showed a similar trend. No significant differences were observed at 42 days and 3 months, but a significant reduction in JIA incidence appeared from 6 months onward (criterion 1: OR = 0.33 at 6 months; p < 0.01). This protective effect remained consistent through 36 months (Tables 2B, Figure 2).

Figure 2
Bar chart showing rates of vaccination and non-vaccination over various time periods. Blue and hatched bars denote Analysis 1, red and crosshatched bars denote Analysis 2. Rates increase over time, with non-vaccination generally higher than vaccination for most intervals.

Figure 2. Rate of new onset juvenile arthritis defined by M08.xx plus drug use at 42 days, 3, 6, 9, 12, 24, 30, 36 months after the event (visit or HPV4 vaccination) in Analysis 1 and Analysis 2.

In Analysis 3, comparing single versus multiple HPV4 doses, no significant differences in JIA incidence were detected across all time points (Table 2C).

3.3 Survival and risk analyses

The cumulative incidence approach based on daily JIA-free probabilities allowed for time-resolved comparisons without assuming a constant hazard ratio over the entire follow-up period. The Cumulative probability curves comparing Cohorts 1 and 2 defined by criteria 1 for Analysis 1 and 2 are presented in Figures 3A, B (log-rank p < 0.0001). In Analysis 3, there was no difference in cumulative JIA incidence between single- and multiple-dose recipients (Figure 3, log-rank p = 0.717).

Figure 3
Three cumulative incidence curve charts compare different cohorts over 36 months. The first and second charts show two cohorts with Cohort 1 (blue) and Cohort 2 (green), highlighting higher cumulative probabilities for Cohort 2. Both indicate statistical significance with a p-value less than 0.001. The third chart compares two dosing regimens for JIA, with 1 dose (blue) and 2 or 3 doses (green), showing similar cumulative probabilities and a non-significant p-value of 0.7117.

Figure 3. Cumulative probability curves for new onset juvenile idiopathic arthritis. (A) the curve derived from Analysis 1, (B) the curve derived from Analysis 2. (C) the curve derived from Analysis 3.

Hazard ratios for JIA development are shown in Figure 4 and Table 3. HPV4 vaccination was associated with reduced hazard across all post-vaccination intervals in Analyses 1 and 2. No significant difference in HRs was found in Analysis 3.

Figure 4
Forest plot comparing hazard ratios over time intervals from 8 days to 3 months, 6 months, 12 months, 18 months, 24 months, 30 months, and 3 years. The plot includes three analyses represented by green, purple, and orange lines with dots indicating hazard ratios and horizontal lines showing 95% confidence intervals. A vertical dashed red line at hazard ratio of 1 marks the reference value for comparison.

Figure 4. Forest Plots showed hazard ratios of new onset JIA in A. Analysis 1, B. Analysis 2 and C. Analysis 3.

Table 3
www.frontiersin.org

Table 3. Hazard ratio of new onset JIA in 3 analyses.

4 Discussion

This large-scale, multicenter retrospective cohort study found no evidence of increased risk of juvenile idiopathic arthritis (JIA) following quadrivalent HPV vaccination. In contrast, our findings consistently demonstrated a protective association between HPV4 vaccination and new-onset JIA, particularly from 12 months to 36 months post-vaccination. These results held true across two distinct analyses, including one that excluded ANA-positive individuals, thus strengthening the validity of the observed protective trend.

The absence of early-onset autoimmune flares (i.e., within 42 days) aligns with current immunological understanding, which suggests that an adaptive autoimmune response generally takes longer to manifest (12). The 42-day threshold is also the standard time frame used by the VAERS system for defining vaccine-associated arthritis (13). Therefore, our observation that JIA onset was significantly lower in the vaccinated cohort beginning at 6 to 12 months post-vaccination is biologically plausible and statistically robust.

Our findings corroborate prior research suggesting a potential autoimmune link between HPV infection and rheumatic diseases, including RA and systemic lupus erythematosus (14, 15). One tertiary center study showed that 48.2% of HPV DNA+ patients had autoantibody or autoimmune disease (16). The mechanism for their condition was generally considered to be via molecular mimicry and cross reactions between self-antigens and viral proteins (17). HPV-47, in particular, contains sequences with molecular mimicry to citrullinated autoantigens involved in RA pathogenesis (18). Moreover, HPV infection has been associated with altered immune cell activation, including Th17 cell expansion and elevated IL-17 levels, which may potentiate local and systemic inflammation (19).

One possible explanation for our findings is that HPV4 vaccination may reduce the systemic immune burden associated with persistent HPV infection, thereby lowering the probability of autoimmune activation in genetically predisposed individuals. This hypothesis is supported by recent NHANES data showing that HPV4 vaccination was associated with reduced RA incidence among adults aged 18–59 years (20).

Although the HPV4 vaccine is not classified as a tolerogenic vaccine, its nature as a non-live virus-like particle (VLP) formulation typically administered in multiple doses may favor regulatory T cell (Treg) induction under low-inflammatory conditions. Previous studies have suggested that repeated antigen exposure from inactivated or subunit vaccines can promote tolerogenic responses by providing persistent but non-replicative stimulation, which facilitates Treg differentiation via IL-10 production, reduced co-stimulation, and antigen presentation in tolerogenic dendritic cell contexts (21). However, in our study, no significant difference in JIA incidence was observed between individuals who received only one dose versus those who completed the full vaccination schedule. This may reflect the relatively small sample size in the single-dose subgroup, limiting the power to detect potential dose-dependent effects. Nonetheless, the possibility that even a single dose of HPV4 vaccine provides sufficient antigen exposure to initiate protective immune modulation cannot be excluded.

Our results are consistent with the known epidemiology of JIA in the United States. Previously reported incidence rates in girls aged 11–15 years range from 13.2 to 30.2 per 100,000 (22). In our study, the observed JIA incidence—based on both diagnostic coding and medication use—fell within or slightly above this range, affirming the external validity of our cohort. Importantly, using a more stringent diagnostic criterion (M08.xx plus drug use) revealed a delayed but more robust protective effect beginning at 6–12 months, likely reflecting improved specificity for true JIA cases. Taken together, the consistency between our observations and prior evidence—including NHANES epidemiologic data, molecular mimicry mechanisms, and the Treg-promoting potential of inactivated vaccines—provides additional support for the biological plausibility of a protective association.

This study has several limitations. First, the TriNetX dataset is predominantly U.S.-based and may not generalize to other populations, particularly Asian cohorts where sample sizes were limited. Second, temporal variation in vaccine administration may confound short-term effects, although this is unlikely to influence long-term outcomes. Third, while the use of both diagnostic codes and medication prescriptions enhances case validation, misclassification remains possible. Not all individuals with an M08.xx code received antirheumatic therapy, suggesting possible overcoding or transient arthropathy. As with all EHR-based studies, the potential for diagnostic misclassification exists; however, our strict case definitions and exclusion of pre-index autoimmune diagnoses or treatments help reduce such risk. Additionally, although extensive matching was performed, residual confounding from unmeasured variables—such as family history, HLA genotype, environmental exposures, and socioeconomic status —cannot be excluded. Although vaccination histories other than HPV were not available for analysis, routine immunizations are limited between ages 9 and 13, with seasonal influenza vaccine being the only widely administered immunization during this period.

In conclusion, this large-scale TriNetX study showed the protective effect which at least one dose of HPV4 vaccination from 12 months to 36 months has on developing JIA. Single or multiple doses of the vaccine showed similar rates of new onset JIA from 42 days to 3 years after HPV vaccination.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics statement

The study protocol was approved by the Institutional Review Board of Taichung Veterans General Hospital (Approval No. CE23480C). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

W-YW: Writing – original draft, Conceptualization. H-HC: Conceptualization, Writing – review & editing. M-CT: Formal analysis, Writing – review & editing. Y-CH: Writing – review & editing. J-PC: Data curation, Writing – review & editing. L-SF: Conceptualization, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

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/fimmu.2025.1621939/full#supplementary-material

Supplementary Table 1 | Drugs in MS109.

Supplementary Table 2 | Diagnosis of ICD 10=M.08.xx.

References

1. Franco EL, Duarte-Franco E, and Ferenczy A. Cervical cancer: Epidemiology, prevention and the role of human papillomavirus infection. Cmaj. (2001) 164:1017–25. doi: 10.1158/1055-9965.EPI-12-1406

PubMed Abstract | Crossref Full Text | Google Scholar

2. Pellegrino P, Carnovale C, Pozzi M, Antoniazzi S, Perrone V, Salvati D, et al. On the relationship between human papilloma virus vaccine and autoimmune diseases. Autoimmun Rev. (2014) 13:736–41. doi: 10.1016/j.autrev.2014.01.054

PubMed Abstract | Crossref Full Text | Google Scholar

3. Kamolratanakul S and Pitisuttithum P. Human papillomavirus vaccine efficacy and effectiveness against cancer. Vaccines (Basel). (2021) 9(12):1413. doi: 10.3390/vaccines9121413

PubMed Abstract | Crossref Full Text | Google Scholar

4. Giuliano AR, Palefsky JM, Goldstone S, Moreira ED Jr., Penny ME, Aranda C, et al. Efficacy of quadrivalent HPV vaccine against HPV Infection and disease in males. N Engl J Med. (2011) 364:401–11. doi: 10.1056/NEJMoa0909537

PubMed Abstract | Crossref Full Text | Google Scholar

5. Hussein HM and Rahal EA. The role of viral infections in the development of autoimmune diseases. Crit Rev Microbiol. (2019) 45:394–412. doi: 10.1080/1040841X.2019.1614904

PubMed Abstract | Crossref Full Text | Google Scholar

6. Rojas M, Restrepo-Jiménez P, Monsalve DM, Pacheco Y, Acosta-Ampudia Y, Ramírez-Santana C, et al. Molecular mimicry and autoimmunity. J Autoimmun. (2018) 95:100–23. doi: 10.1016/j.jaut.2018.10.012

PubMed Abstract | Crossref Full Text | Google Scholar

7. Sogkas G, Atschekzei F, Adriawan IR, Dubrowinskaja N, Witte T, and Schmidt RE. Cellular and molecular mechanisms breaking immune tolerance in inborn errors of immunity. Cell Mol Immunol. (2021) 18:1122–40. doi: 10.1038/s41423-020-00626-z

PubMed Abstract | Crossref Full Text | Google Scholar

8. Geier DA and Geier MR. Quadrivalent human papillomavirus vaccine and autoimmune adverse events: A case-control assessment of the vaccine adverse event reporting system (VAERS) database. Immunol Res. (2017) 65:46–54. doi: 10.1007/s12026-016-8815-9

PubMed Abstract | Crossref Full Text | Google Scholar

9. Maldonado I, Rodríguez Niño N, Valencia CF, Ortiz D, Ayala N, D’Silva L, et al. Evaluation of the safety profile of the quadrivalent vaccine against human papillomavirus in the risk of developing autoimmune, neurological, and hematological diseases in adolescent women in Colombia. Vaccine. (2024) 42:2414–20. doi: 10.1016/j.vaccine.2024.02.085

PubMed Abstract | Crossref Full Text | Google Scholar

10. Chao C, Klein NP, Velicer CM, Sy LS, Slezak JM, Takhar H, et al. Surveillance of autoimmune conditions following routine use of quadrivalent human papillomavirus vaccine. J Intern Med. (2012) 271:193–203. doi: 10.1111/j.1365-2796.2011.02467.x

PubMed Abstract | Crossref Full Text | Google Scholar

11. Topaloglu U and Palchuk MB. Using a federated network of real-world data to optimize clinical trials operations. JCO Clin Cancer Inform. (2018) 2:1–10. doi: 10.1200/CCI.17.00067

PubMed Abstract | Crossref Full Text | Google Scholar

12. Gil-Vila A, Ravichandran N, Selva-O’Callaghan A, Sen P, Nune A, Gaur PS, et al. COVID-19 Vaccination in Autoimmune Diseases (COVAD) study: Vaccine safety in idiopathic inflammatory myopathies. Muscle Nerve. (2022) 66:426–37. doi: 10.1002/mus.27681

PubMed Abstract | Crossref Full Text | Google Scholar

13. VAERS table of reportable events following vaccination(2024). Available online at: https://vaers.hhs.gov/docs/VAERS_Table_of_Reportable_Events_Following_Vaccination.pdf. (Accessed May 6, 2024)

Google Scholar

14. Liu YT, Tsou HK, Chiou JY, Wang YH, Chou MC, and Wei JC. Association of human papillomavirus infection with risk for rheumatoid arthritis: A population-based cohort study. Ann Rheum Dis. (2019) 78:1734–6. doi: 10.1136/annrheumdis-2019-215931

PubMed Abstract | Crossref Full Text | Google Scholar

15. Shi LH, Huang JY, Liu YZ, Chiou JY, Wu R, and Wei JC. Risk of systemic lupus erythematosus in patients with human papillomavirus infection: A population-based retrospective cohort study. Lupus. (2018) 27:2279–83. doi: 10.1177/0961203318809179

PubMed Abstract | Crossref Full Text | Google Scholar

16. Donmez HG, Tanacan A, Unal C, Fadiloglu E, Onder SC, Portakal O, et al. Human papillomavirus infection and autoimmune disorders: A tertiary center experience. Pathog Dis. (2019) 77(3):ftz028. doi: 10.1093/femspd/ftz028

PubMed Abstract | Crossref Full Text | Google Scholar

17. Smolen JS, Aletaha D, and McInnes IB. Rheumatoid arthritis. Lancet. (2016) 388:2023–38. doi: 10.1016/S0140-6736(16)30173-8

PubMed Abstract | Crossref Full Text | Google Scholar

18. Shi J, Sun X, Zhao Y, Zhao J, and Li Z. Prevalence and significance of antibodies to citrullinated human papilloma virus-47 E2345–362 in rheumatoid arthritis. J Autoimmun. (2008) 31:131–5. doi: 10.1016/j.jaut.2008.04.021

PubMed Abstract | Crossref Full Text | Google Scholar

19. Xue JS, Wang YL, Chen C, Zhu XJ, Zhu H, and Hu Y. Effects of Th17 cells and IL-17 in the progression of cervical carcinogenesis with high-risk human papillomavirus infection. Cancer Med. (2018) 7:297–306. doi: 10.1002/cam4.1279

PubMed Abstract | Crossref Full Text | Google Scholar

20. Yang G, Ren Z, and Wang K. Association between human papillomavirus infection or immunization and risk for rheumatoid arthritis. Front Immunol. (2023) 14:1130217. doi: 10.3389/fimmu.2023.1130217

PubMed Abstract | Crossref Full Text | Google Scholar

21. Prietl B, Treiber G, Chaparro N, Graninger M, and Obermayer-Pietsch B. Harnessing regulatory T cells through vaccination strategies: a path toward immune tolerance. Vaccines (Basel). (2024) 12:992. doi: 10.3390/vaccines12030992

Crossref Full Text | Google Scholar

22. Harrold LR, Salman C, Shoor S, Curtis JR, Asgari MM, Gelfand JM, et al. Incidence and prevalence of juvenile idiopathic arthritis among children in a managed care population, 1996-2009. J Rheumatol. (2013) 40:1218–25. doi: 10.3899/jrheum.120661

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: quadrivalent human papillomavirus vaccine (HPV4 vaccine), juvenile idiopathic arthritis, TriNetX, vaccine safety, cohort study

Citation: Wu W-Y, Chen H-H, Tsai M-C, Huang Y-C, Chen J-P and Fu L-S (2025) Risk of developing juvenile idiopathic arthritis after quadrivalent HPV vaccination: a retrospective cohort study using the TriNetX U.S. Network. Front. Immunol. 16:1621939. doi: 10.3389/fimmu.2025.1621939

Received: 02 May 2025; Accepted: 04 August 2025;
Published: 25 August 2025.

Edited by:

Pei-Hui Wang, Shandong University, China

Reviewed by:

Kai Wang, The Affiliated Huaian No.1 People’s Hospital of Nanjing Medical University, China
Julio Cesar Gonzalez Rodriguez, National Institute of Cancerology (INCAN), Mexico

Copyright © 2025 Wu, Chen, Tsai, Huang, Chen and Fu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Lin-Shien Fu, bGluc2hpZW5mdUB5YWhvby5jb20udHc=

These authors have contributed equally to this work

Disclaimer: 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.