<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Immunol.</journal-id>
<journal-title>Frontiers in Immunology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Immunol.</abbrev-journal-title>
<issn pub-type="epub">1664-3224</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fimmu.2025.1621939</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Immunology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Risk of developing juvenile idiopathic arthritis after quadrivalent HPV vaccination: a retrospective cohort study using the TriNetX U.S. Network</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Wu</surname>
<given-names>Wen-Yu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3053878/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Chen</surname>
<given-names>Hsin-Hua</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<xref ref-type="aff" rid="aff6">
<sup>6</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/806664/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Tsai</surname>
<given-names>Ming-Chin</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2019800/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Huang</surname>
<given-names>Yung-Chieh</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff7">
<sup>7</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2398302/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Chen</surname>
<given-names>Jun-Peng</given-names>
</name>
<xref ref-type="aff" rid="aff8">
<sup>8</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Fu</surname>
<given-names>Lin-Shien</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff7">
<sup>7</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Pediatrics, Taichung Veterans General Hospital</institution>, <addr-line>Taichung</addr-line>,&#xa0;<country>Taiwan</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Division of Allergy, Immunology and Rheumatology, Department of Medicine, Taichung Veterans General Hospital</institution>, <addr-line>Taichung</addr-line>,&#xa0;<country>Taiwan</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Division of Clinical Informatics, Center for Quality Management, Taichung Veterans General Hospital</institution>, <addr-line>Taichung</addr-line>,&#xa0;<country>Taiwan</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University</institution>, <addr-line>Taichung</addr-line>,&#xa0;<country>Taiwan</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>Program in Translational Medicine and Rong Hsing Research Center for Translational Medicine, National Chung Hsing University</institution>, <addr-line>Taichung</addr-line>,&#xa0;<country>Taiwan</country>
</aff>
<aff id="aff6">
<sup>6</sup>
<institution>Big Data Center, National Chung Hsing University</institution>, <addr-line>Taichung</addr-line>,&#xa0;<country>Taiwan</country>
</aff>
<aff id="aff7">
<sup>7</sup>
<institution>Department of Pediatrics, National Yang Ming Chiao Tung University</institution>, <addr-line>Taipei</addr-line>,&#xa0;<country>Taiwan</country>
</aff>
<aff id="aff8">
<sup>8</sup>
<institution>Department of Research, Taichung Veterans General Hospital</institution>, <addr-line>Taichung</addr-line>,&#xa0;<country>Taiwan</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Pei-Hui Wang, Shandong University, China</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1269902/overview">Kai Wang</ext-link>, The Affiliated Huaian No.1 People&#x2019;s Hospital of Nanjing Medical University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2548914/overview">Julio Cesar Gonzalez Rodriguez</ext-link>, National Institute of Cancerology (INCAN), Mexico</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Lin-Shien Fu, <email xlink:href="mailto:linshienfu@yahoo.com.tw">linshienfu@yahoo.com.tw</email>
</p>
</fn>
<fn fn-type="equal" id="fn003">
<p>&#x2020;These authors have contributed equally to this work</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>25</day>
<month>08</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>16</volume>
<elocation-id>1621939</elocation-id>
<history>
<date date-type="received">
<day>02</day>
<month>05</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>04</day>
<month>08</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2025 Wu, Chen, Tsai, Huang, Chen and Fu.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Wu, Chen, Tsai, Huang, Chen and Fu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>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.</p>
</sec>
<sec>
<title>Methods</title>
<p>We conducted a retrospective cohort study using data from the TriNetX U.S. Collaborative Network. Females aged 9&#x2013;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).</p>
</sec>
<sec>
<title>Results</title>
<p>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&#x2013;0.52). No difference in JIA risk was observed between single and multiple doses in Analysis 3 (n=20,822 pairs). Early-onset JIA (&lt;6 months after HPV4 vaccine) showed inconsistent trends, with only limited protective signals.</p>
</sec>
<sec>
<title>Conclusions</title>
<p>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&#x2013;12 months and lasting for at least 3 years. These findings support the safety and possible immunoregulatory benefit of HPV4 in adolescents.</p>
</sec>
</abstract>
<kwd-group>
<kwd>quadrivalent human papillomavirus vaccine (HPV4 vaccine)</kwd>
<kwd>juvenile idiopathic arthritis</kwd>
<kwd>TriNetX</kwd>
<kwd>vaccine safety</kwd>
<kwd>cohort study</kwd>
</kwd-group>
<counts>
<fig-count count="4"/>
<table-count count="3"/>
<equation-count count="0"/>
<ref-count count="22"/>
<page-count count="10"/>
<word-count count="4414"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Vaccines and Molecular Therapeutics</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Human papillomaviruses (HPVs) are double-stranded DNA viruses associated with a wide spectrum of human diseases, including cervical cancer and autoimmune conditions (<xref ref-type="bibr" rid="B1">1</xref>, <xref ref-type="bibr" rid="B2">2</xref>). Prophylactic HPV vaccination, particularly the quadrivalent HPV4 vaccine, has demonstrated robust efficacy in preventing HPV-related infections and malignancies (<xref ref-type="bibr" rid="B3">3</xref>, <xref ref-type="bibr" rid="B4">4</xref>). Despite its favorable safety profile, concerns remain regarding its potential role in triggering autoimmune disorders such as juvenile idiopathic arthritis (JIA).</p>
<p>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 (<xref ref-type="bibr" rid="B5">5</xref>). Molecular mimicry and immune dysregulation have been proposed as key mechanisms by which pathogens may initiate or exacerbate autoimmunity (<xref ref-type="bibr" rid="B6">6</xref>, <xref ref-type="bibr" rid="B7">7</xref>).</p>
<p>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 (<xref ref-type="bibr" rid="B8">8</xref>). Similarly, a Colombian study employing inverse probability of treatment weighting also observed increased arthritis risk (<xref ref-type="bibr" rid="B9">9</xref>). 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 (<xref ref-type="bibr" rid="B10">10</xref>).</p>
<p>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 (<xref ref-type="bibr" rid="B11">11</xref>). 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.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Data sources</title>
<p>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:</p>
<list list-type="bullet">
<list-item>
<p>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 (<xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary Tables&#xa0;1</bold>
</xref>, <xref ref-type="supplementary-material" rid="SF2">
<bold>2</bold>
</xref>), were excluded.</p>
</list-item>
<list-item>
<p>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.</p>
</list-item>
<list-item>
<p>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&#x2013;12 month dosing interval.</p>
</list-item>
</list>
<p>TriNetX&#x2019;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&#x2013;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.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Study design</title>
<p>As shown in <xref ref-type="fig" rid="f1">
<bold>Figures 1A, B</bold>
</xref>, 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.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>
<bold>(A)</bold> 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. <bold>(B)</bold> 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.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-16-1621939-g001.tif">
<alt-text content-type="machine-generated">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.</alt-text>
</graphic>
</fig>
<p>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.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Ethical approval</title>
<p>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.</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Statistical analysis</title>
<p>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&#x2013;M99), and NSAID use. Time-to-event analysis was performed using the Kaplan&#x2013;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.</p>
<p>Due to the platform&#x2019;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 &lt; 0.05 was considered statistically significant.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Cohort characteristics</title>
<p>The demographic and baseline clinical characteristics of the matched cohorts for all three analyses are presented in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>. In Analysis 1, 55,257 participants were included in both the HPV4-vaccinated and unvaccinated cohorts, with a mean age of 11.7 &#xb1; 1.0 years. The racial distribution was 51.0&#x2013;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.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Baseline characteristics of study subjects after matching.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" rowspan="3" align="left">Variables</th>
<th valign="middle" colspan="4" align="left">Analysis 1</th>
<th valign="middle" colspan="4" align="left">Analysis 2</th>
<th valign="middle" colspan="4" align="left">Analysis 3</th>
</tr>
<tr>
<th valign="middle" align="left">Vaccination group (n= 55,257)</th>
<th valign="middle" align="left">Non-vaccination (n= 55,257)</th>
<th valign="middle" rowspan="2" align="left">p-value</th>
<th valign="middle" rowspan="2" align="left">SMD</th>
<th valign="middle" align="left">Vaccination group (n= 53,827)</th>
<th valign="middle" align="left">Non-vaccination (n= 53,827)</th>
<th valign="middle" rowspan="2" align="left">p-value</th>
<th valign="middle" rowspan="2" align="left">SMD</th>
<th valign="middle" align="left">Vaccination 1/time (n= 20,822)</th>
<th valign="middle" align="left">Vaccination&gt;1/time (n= 20,822)</th>
<th valign="middle" rowspan="2" align="left">p-value</th>
<th valign="middle" rowspan="2" align="left">SMD</th>
</tr>
<tr>
<th valign="middle" align="left">n (%) or mean &#xb1; SD</th>
<th valign="middle" align="left">n (%) or mean &#xb1; SD</th>
<th valign="middle" align="left">n (%) or mean &#xb1; SD</th>
<th valign="middle" align="left">n (%) or mean &#xb1; SD</th>
<th valign="middle" align="left">n (%) or mean &#xb1; SD</th>
<th valign="middle" align="left">n (%) or mean &#xb1; SD</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">Age at index (y/o)</td>
<td valign="middle" align="left">11.7 &#xb1; 1.0</td>
<td valign="middle" align="left">11.7 &#xb1; 1.0</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">11.7 &#xb1; 1.0</td>
<td valign="middle" align="left">11.7 &#xb1; 1.0</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">11.7 &#xb1; 1.0</td>
<td valign="middle" align="left">11.7 &#xb1; 1.0</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
</tr>
<tr>
<th valign="middle" colspan="13" align="left">Race or ethnicity</th>
</tr>
<tr>
<td valign="middle" align="left">White</td>
<td valign="middle" align="left">27,917 (51.0)</td>
<td valign="middle" align="left">28,005 (51.1)</td>
<td valign="middle" align="left">0.6</td>
<td valign="middle" align="left">0.003</td>
<td valign="middle" align="left">27,392 (50.9)</td>
<td valign="middle" align="left">27,470 (51.0)</td>
<td valign="middle" align="left">0.6</td>
<td valign="middle" align="left">0.003</td>
<td valign="middle" align="left">10,254 (49.2)</td>
<td valign="middle" align="left">10,245 (49.2)</td>
<td valign="middle" align="left">0.93</td>
<td valign="middle" align="left">0.001</td>
</tr>
<tr>
<td valign="middle" align="left">Black or African American</td>
<td valign="middle" align="left">13,625 (24.9)</td>
<td valign="middle" align="left">13,625 (24.9)</td>
<td valign="middle" align="left"/>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">13,417 (24.9)</td>
<td valign="middle" align="left">13,417 (24.9)</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">5,485 (26.3)</td>
<td valign="middle" align="left">5,437 (26.1)</td>
<td valign="middle" align="left">0.59</td>
<td valign="middle" align="left">0.005</td>
</tr>
<tr>
<td valign="middle" align="left">Asian</td>
<td valign="middle" align="left">1,468 (2.7)</td>
<td valign="middle" align="left">1,468 (2.7)</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1,437 (2.7)</td>
<td valign="middle" align="left">1,437 (2.7)</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">495 (2.4)</td>
<td valign="middle" align="left">498 (2.4)</td>
<td valign="middle" align="left">0.92</td>
<td valign="middle" align="left">0.001</td>
</tr>
<tr>
<th valign="middle" colspan="13" align="left">Diagnosis</th>
</tr>
<tr>
<td valign="middle" align="left">Disease of the musculoskeletal<break/>system and connective tissue</td>
<td valign="middle" align="left">10,349 (18.9)</td>
<td valign="middle" align="left">10,349 (18.9)</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">9,927 (18.4)</td>
<td valign="middle" align="left">9,927 (18.4)</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">3,472 (16.7)</td>
<td valign="middle" align="left">3,460 (16.6)</td>
<td valign="middle" align="left">0.88</td>
<td valign="middle" align="left">0.002</td>
</tr>
<tr>
<th valign="middle" colspan="13" align="left">Medication</th>
</tr>
<tr>
<td valign="middle" align="left">NSAIs</td>
<td valign="middle" align="left">8,679 (15.8)</td>
<td valign="middle" align="left">8,679 (15.8)</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">8,295 (15.4)</td>
<td valign="middle" align="left">8,295 (15.4)</td>
<td valign="middle" align="left">1</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">2,518 (12.1)</td>
<td valign="middle" align="left">2,555 (12.3)</td>
<td valign="middle" align="left">0.58</td>
<td valign="middle" align="left">0.005</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>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.</p>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Incidence of new-onset JIA</title>
<p>The cumulative incidence and odds ratios (ORs) of new-onset JIA over 36 months are summarized in <xref ref-type="table" rid="T2">
<bold>Tables&#xa0;2A&#x2013;C</bold>
</xref>. 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 <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>.</p>
<table-wrap-group id="T2" position="float">
<table-wrap>
<label>Table&#xa0;2</label>
<caption>
<p>Incidence and odds ratio of new onset JIA and JIA with drug.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" colspan="11" align="left">A. Analysis 1</th>
</tr>
<tr>
<th valign="middle" rowspan="2" align="center">Analysis 1</th>
<th valign="middle" colspan="5" align="left">JIA</th>
<th valign="middle" colspan="5" align="left">JIA with Drug</th>
</tr>
<tr>
<th valign="middle" align="left">Vaccination (%)</th>
<th valign="middle" align="left">Non-vaccination (%)</th>
<th valign="middle" align="left">OR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
<th valign="middle" align="left">FDR p-value</th>
<th valign="middle" align="left">Vaccination (%)</th>
<th valign="middle" align="left">Non-vaccination (%)</th>
<th valign="middle" align="left">OR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
<th valign="middle" align="left">FDR p-value</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">8 days to 42 days</td>
<td valign="middle" align="right">10 (0.018)</td>
<td valign="middle" align="right">14 (0.026)</td>
<td valign="middle" align="right">0.71 (0.32-1.60)</td>
<td valign="middle" align="right">0.4153</td>
<td valign="middle" align="right">0.4153</td>
<td valign="middle" align="right">10 (0.018)</td>
<td valign="middle" align="right">0 (0.000)</td>
<td valign="middle" align="right">20.99(1.23-358.19)</td>
<td valign="middle" align="right">0.0355</td>
<td valign="middle" align="right">0.0473</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 3 months</td>
<td valign="middle" align="right">10 (0.018)</td>
<td valign="middle" align="right">22 (0.040)</td>
<td valign="middle" align="right">0.45 (0.22-0.96)</td>
<td valign="middle" align="right">0.0385</td>
<td valign="middle" align="right">0.0440</td>
<td valign="middle" align="right">9 (0.016)</td>
<td valign="middle" align="right">4 (0.007)</td>
<td valign="middle" align="right">2.25 (0.69-7.30)</td>
<td valign="middle" align="right">0.1776</td>
<td valign="middle" align="right">0.2030</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 6 months</td>
<td valign="middle" align="right">10 (0.018)</td>
<td valign="middle" align="right">30 (0.055)</td>
<td valign="middle" align="right">0.33 (0.16-0.68)</td>
<td valign="middle" align="right">0.0026</td>
<td valign="middle" align="right">0.0035</td>
<td valign="middle" align="right">9 (0.016)</td>
<td valign="middle" align="right">12 (0.022)</td>
<td valign="middle" align="right">0.75 (0.32-1.78)</td>
<td valign="middle" align="right">0.5130</td>
<td valign="middle" align="right">0.5130</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 12 months</td>
<td valign="middle" align="right">10 (0.018)</td>
<td valign="middle" align="right">42 (0.077)</td>
<td valign="middle" align="right">0.24 (0.12-0.47)</td>
<td valign="middle" align="right">&lt;0.0001</td>
<td valign="middle" align="right">0.0003</td>
<td valign="middle" align="right">6 (0.011)</td>
<td valign="middle" align="right">22 (0.040)</td>
<td valign="middle" align="right">0.27 (0.11-0.67)</td>
<td valign="middle" align="right">0.0048</td>
<td valign="middle" align="right">0.0122</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 18 months</td>
<td valign="middle" align="right">19 (0.035)</td>
<td valign="middle" align="right">52 (0.095)</td>
<td valign="middle" align="right">0.37 (0.22-0.62)</td>
<td valign="middle" align="right">0.0002</td>
<td valign="middle" align="right">0.0004</td>
<td valign="middle" align="right">13 (0.024)</td>
<td valign="middle" align="right">32 (0.058)</td>
<td valign="middle" align="right">0.41 (0.21-0.77)</td>
<td valign="middle" align="right">0.0061</td>
<td valign="middle" align="right">0.0122</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 24 months</td>
<td valign="middle" align="right">26 (0.047)</td>
<td valign="middle" align="right">61 (0.111)</td>
<td valign="middle" align="right">0.43 (0.27-0.67)</td>
<td valign="middle" align="right">0.0003</td>
<td valign="middle" align="right">0.0005</td>
<td valign="middle" align="right">19 (0.035)</td>
<td valign="middle" align="right">37 (0.068)</td>
<td valign="middle" align="right">0.51 (0.30-0.89)</td>
<td valign="middle" align="right">0.0180</td>
<td valign="middle" align="right">0.0288</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 30 months</td>
<td valign="middle" align="right">32 (0.058)</td>
<td valign="middle" align="right">79 (0.144)</td>
<td valign="middle" align="right">0.40 (0.27-0.61)</td>
<td valign="middle" align="right">&lt;0.0001</td>
<td valign="middle" align="right">0.0003</td>
<td valign="middle" align="right">26 (0.047)</td>
<td valign="middle" align="right">54 (0.099)</td>
<td valign="middle" align="right">0.48 (0.30-0.77)</td>
<td valign="middle" align="right">0.0022</td>
<td valign="middle" align="right">0.0122</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 3 years</td>
<td valign="middle" align="right">37 (0.068)</td>
<td valign="middle" align="right">83 (0.152)</td>
<td valign="middle" align="right">0.45 (0.30-0.66)</td>
<td valign="middle" align="right">&lt;0.0001</td>
<td valign="middle" align="right">0.0003</td>
<td valign="middle" align="right">31 (0.057)</td>
<td valign="middle" align="right">59 (0.108)</td>
<td valign="middle" align="right">0.52 (0.34-0.81)</td>
<td valign="middle" align="right">0.0037</td>
<td valign="middle" align="right">0.0122</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>JIA, Juvenile idiopathic arthritis.</p>
</fn>
<fn>
<p>JIA with Drug: JIA plus at least one drug use. The drugs included non-salicylate drugs, methotrexate and other anti-rheumatic drugs shown in Supplementary data.</p>
</fn>
<fn>
<p>FDR, False Discovery Rate (FDR) method proposed by Benjamini and Hochberg.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" colspan="11" align="left">B. Analysis 2</th>
</tr>
<tr>
<th valign="middle" rowspan="2" align="left">Analysis 2</th>
<th valign="middle" colspan="5" align="left">JIA</th>
<th valign="middle" colspan="5" align="left">JIA with Drug</th>
</tr>
<tr>
<th valign="middle" align="left">Vaccination (%)</th>
<th valign="middle" align="left">Non-vaccination (%)</th>
<th valign="middle" align="left">OR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
<th valign="middle" align="left">FDR p-value</th>
<th valign="middle" align="left">Vaccination (%)</th>
<th valign="middle" align="left">Non-vaccination (%)</th>
<th valign="middle" align="left">OR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
<th valign="middle" align="left">FDR p-value</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">8 days to 42 days</td>
<td valign="middle" align="right">10 (0.019)</td>
<td valign="middle" align="right">12 (0.022)</td>
<td valign="middle" align="right">0.83 (0.36-1.93)</td>
<td valign="middle" align="right">0.6691</td>
<td valign="middle" align="right">0.6691</td>
<td valign="middle" align="right">10 (0.019)</td>
<td valign="middle" align="right">0 (0.007)</td>
<td valign="middle" align="right">2.45 (0.78-7.97)</td>
<td valign="middle" align="right">0.1217</td>
<td valign="middle" align="right">0.1623</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 3 months</td>
<td valign="middle" align="right">10 (0.019)</td>
<td valign="middle" align="right">21 (0.039)</td>
<td valign="middle" align="right">0.48 (0.22-1.01)</td>
<td valign="middle" align="right">0.0532</td>
<td valign="middle" align="right">0.0608</td>
<td valign="middle" align="right">9 (0.017)</td>
<td valign="middle" align="right">7 (0.013)</td>
<td valign="middle" align="right">1.28 (0.48-3.45)</td>
<td valign="middle" align="right">0.6189</td>
<td valign="middle" align="right">0.7073</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 6 months</td>
<td valign="middle" align="right">10 (0.019)</td>
<td valign="middle" align="right">30 (0.056)</td>
<td valign="middle" align="right">0.33 (0.16-0.68)</td>
<td valign="middle" align="right">0.0026</td>
<td valign="middle" align="right">0.0035</td>
<td valign="middle" align="right">9 (0.017)</td>
<td valign="middle" align="right">9 (0.017)</td>
<td valign="middle" align="right">0.99 (0.40-2.52)</td>
<td valign="middle" align="right">0.9989</td>
<td valign="middle" align="right">0.9989</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 12 months</td>
<td valign="middle" align="right">10 (0.019)</td>
<td valign="middle" align="right">47 (0.087)</td>
<td valign="middle" align="right">0.21 (0.11-0.42)</td>
<td valign="middle" align="right">&lt;0.0001</td>
<td valign="middle" align="right">0.0002</td>
<td valign="middle" align="right">6 (0.011)</td>
<td valign="middle" align="right">22 (0.041)</td>
<td valign="middle" align="right">0.27 (0.11-0.67)</td>
<td valign="middle" align="right">0.0048</td>
<td valign="middle" align="right">0.0096</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 18 months</td>
<td valign="middle" align="right">19 (0.035)</td>
<td valign="middle" align="right">59 (0.110)</td>
<td valign="middle" align="right">0.33 (0.20-0.55)</td>
<td valign="middle" align="right">&lt;0.0001</td>
<td valign="middle" align="right">0.0002</td>
<td valign="middle" align="right">13 (0.024)</td>
<td valign="middle" align="right">32 (0.059)</td>
<td valign="middle" align="right">0.41 (0.21-0.77)</td>
<td valign="middle" align="right">0.0061</td>
<td valign="middle" align="right">0.0098</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 24 months</td>
<td valign="middle" align="right">26 (0.048)</td>
<td valign="middle" align="right">72 (0.134)</td>
<td valign="middle" align="right">0.36 (0.23-0.57)</td>
<td valign="middle" align="right">&lt;0.0001</td>
<td valign="middle" align="right">0.0002</td>
<td valign="middle" align="right">19 (0.035)</td>
<td valign="middle" align="right">43 (0.080)</td>
<td valign="middle" align="right">0.44 (0.26-0.76)</td>
<td valign="middle" align="right">0.003</td>
<td valign="middle" align="right">0.0096</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 30 months</td>
<td valign="middle" align="right">32 (0.059)</td>
<td valign="middle" align="right">83 (0.154)</td>
<td valign="middle" align="right">0.39 (0.26-0.58)</td>
<td valign="middle" align="right">&lt;0.0001</td>
<td valign="middle" align="right">0.0002</td>
<td valign="middle" align="right">26 (0.048)</td>
<td valign="middle" align="right">52 (0.097)</td>
<td valign="middle" align="right">0.50 (0.32-0.80)</td>
<td valign="middle" align="right">0.0039</td>
<td valign="middle" align="right">0.0096</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 3 years</td>
<td valign="middle" align="right">36 (0.067)</td>
<td valign="middle" align="right">91 (0.169)</td>
<td valign="middle" align="right">0.40 (0.27-0.28)</td>
<td valign="middle" align="right">&lt;0.0001</td>
<td valign="middle" align="right">0.0002</td>
<td valign="middle" align="right">30 (0.056)</td>
<td valign="middle" align="right">61 (0.113)</td>
<td valign="middle" align="right">0.49 (0.32-0.76)</td>
<td valign="middle" align="right">0.0014</td>
<td valign="middle" align="right">0.0096</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>JIA, Juvenile idiopathic arthritis.</p>
</fn>
<fn>
<p>JIA with Drug: JIA plus at least one drug use. The drugs included non-salicylate drugs, methotrexate and other anti-rheumatic drugs shown in Supplementary data.</p>
</fn>
<fn>
<p>FDR, False Discovery Rate (FDR) method proposed by Benjamini and Hochberg.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" colspan="11" align="left">C. Analysis 3</th>
</tr>
<tr>
<th valign="middle" rowspan="2" align="left">Analysis 3</th>
<th valign="middle" colspan="5" align="left">JIA</th>
<th valign="middle" colspan="5" align="left">JIA with Drug</th>
</tr>
<tr>
<th valign="middle" align="left">Vaccination 1 dose (%)</th>
<th valign="middle" align="left">Vaccination 2 or 3 doses (%)</th>
<th valign="middle" align="left">OR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
<th valign="middle" align="left">FDR p-value</th>
<th valign="middle" align="left">Vaccination 1 dose (%)</th>
<th valign="middle" align="left">Vaccination 2 or 3 doses (%)</th>
<th valign="middle" align="left">OR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
<th valign="middle" align="left">FDR p-value</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">8 days to 42 days</td>
<td valign="middle" align="right">0 (0.000)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">0.05 (0.003-0.81)</td>
<td valign="middle" align="right">0.0354</td>
<td valign="middle" align="right">0.2832</td>
<td valign="middle" align="right">0 (0.000)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">0.05 (0.003-0.81)</td>
<td valign="middle" align="right">0.0354</td>
<td valign="middle" align="right">0.2832</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 3 months</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">1.00 (0.42-2.40)</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">9 (0.043)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">0.90 (0.37-2.21)</td>
<td valign="middle" align="right">0.8183</td>
<td valign="middle" align="right">0.9352</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 6 months</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">1.00 (0.42-2.40)</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">9 (0.043)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">0.90 (0.37-2.21)</td>
<td valign="middle" align="right">0.8183</td>
<td valign="middle" align="right">0.9352</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 12 months</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">1.00 (0.42-2.40)</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">9 (0.043)</td>
<td valign="middle" align="right">7 (0.034)</td>
<td valign="middle" align="right">1.29 (0.48-3.45)</td>
<td valign="middle" align="right">0.6181</td>
<td valign="middle" align="right">0.9352</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 18 months</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">1.00 (0.42-2.40)</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">7 (0.034)</td>
<td valign="middle" align="right">7 (0.034)</td>
<td valign="middle" align="right">1.00 (0.35-2.85)</td>
<td valign="middle" align="right">0.9998</td>
<td valign="middle" align="right">0.9998</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 24 months</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">13 (0.062)</td>
<td valign="middle" align="right">0.77 (0.34-1.75)</td>
<td valign="middle" align="right">0.5325</td>
<td valign="middle" align="right">0.9997</td>
<td valign="middle" align="right">7 (0.034)</td>
<td valign="middle" align="right">8 (0.038)</td>
<td valign="middle" align="right">0.87 (0.32-2.41)</td>
<td valign="middle" align="right">0.7962</td>
<td valign="middle" align="right">0.9352</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 30 months</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">15 (0.072)</td>
<td valign="middle" align="right">0.67 (0.30-1.48)</td>
<td valign="middle" align="right">0.3203</td>
<td valign="middle" align="right">0.8541</td>
<td valign="middle" align="right">8 (0.038)</td>
<td valign="middle" align="right">10 (0.048)</td>
<td valign="middle" align="right">0.80 (0.32-2.03)</td>
<td valign="middle" align="right">0.6378</td>
<td valign="middle" align="right">0.9352</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 3 years</td>
<td valign="middle" align="right">11 (0.053)</td>
<td valign="middle" align="right">18 (0.086)</td>
<td valign="middle" align="right">0.61 (0.29-1.29)</td>
<td valign="middle" align="right">0.1979</td>
<td valign="middle" align="right">0.7916</td>
<td valign="middle" align="right">8 (0.038)</td>
<td valign="middle" align="right">13 (0.062)</td>
<td valign="middle" align="right">0.62 (0.24-1.48)</td>
<td valign="middle" align="right">0.2797</td>
<td valign="middle" align="right">0.9352</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>JIA, Juvenile idiopathic arthritis.</p>
</fn>
<fn>
<p>JIA with Drug: JIA plus at least one drug use. The drugs included non-salicylate drugs, methotrexate and other anti-rheumatic drugs shown in Supplementary data.</p>
</fn>
<fn>
<p>FDR, False Discovery Rate (FDR) method proposed by Benjamini and Hochberg.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</table-wrap-group>
<p>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&#x2013;1.60; criterion 2: OR = 20.99, 95% CI = 1.23&#x2013;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 &lt; 0.05 for all time points).</p>
<p>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 &lt; 0.01). This protective effect remained consistent through 36 months (<xref ref-type="table" rid="T2">
<bold>Tables&#xa0;2B</bold>
</xref>, <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>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.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-16-1621939-g002.tif">
<alt-text content-type="machine-generated">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.</alt-text>
</graphic>
</fig>
<p>In Analysis 3, comparing single versus multiple HPV4 doses, no significant differences in JIA incidence were detected across all time points (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2C</bold>
</xref>).</p>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Survival and risk analyses</title>
<p>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 <xref ref-type="fig" rid="f3">
<bold>Figures&#xa0;3A, B</bold>
</xref> (log-rank p &lt; 0.0001). In Analysis 3, there was no difference in cumulative JIA incidence between single- and multiple-dose recipients (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>, log-rank p = 0.717).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Cumulative probability curves for new onset juvenile idiopathic arthritis. <bold>(A)</bold> the curve derived from Analysis 1, <bold>(B)</bold> the curve derived from Analysis 2. <bold>(C)</bold> the curve derived from Analysis 3.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-16-1621939-g003.tif">
<alt-text content-type="machine-generated">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.</alt-text>
</graphic>
</fig>
<p>Hazard ratios for JIA development are shown in <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref> and <xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>. 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.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Forest Plots showed hazard ratios of new onset JIA in A. Analysis 1, B. Analysis 2 and C. Analysis 3.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-16-1621939-g004.tif">
<alt-text content-type="machine-generated">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.</alt-text>
</graphic>
</fig>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Hazard ratio of new onset JIA in 3 analyses.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" rowspan="2" align="left">Time</th>
<th valign="middle" colspan="2" align="left">Analysis 1</th>
<th valign="middle" colspan="2" align="left">Analysis 2</th>
<th valign="middle" colspan="2" align="left">Analysis 3</th>
</tr>
<tr>
<th valign="middle" align="left">HR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
<th valign="middle" align="left">HR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
<th valign="middle" align="left">HR (95% CI)</th>
<th valign="middle" align="left">p-value</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">8 days to 3 months</td>
<td valign="middle" align="right">0.120 (0.036-0.400)</td>
<td valign="middle" align="right">0.066</td>
<td valign="middle" align="right">0.125 (0.037-0.420)</td>
<td valign="middle" align="right">0.421</td>
<td valign="middle" align="right">0.529 (0.048-5.833)</td>
<td valign="middle" align="right">0.701</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 6 months</td>
<td valign="middle" align="right">0.087 (0.027-0.286)</td>
<td valign="middle" align="right">0.970</td>
<td valign="middle" align="right">0.087 (0.027-0.285)</td>
<td valign="middle" align="right">0.718</td>
<td valign="middle" align="right">0.529 (0.048-5.833)</td>
<td valign="middle" align="right">0.905</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 12 months</td>
<td valign="middle" align="right">0.206 (0.103-0.410)</td>
<td valign="middle" align="right">0.017</td>
<td valign="middle" align="right">0.183 (0.092-0.362)</td>
<td valign="middle" align="right">0.137</td>
<td valign="middle" align="right">0.425 (0.082-2.191)</td>
<td valign="middle" align="right">0.287</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 18 months</td>
<td valign="middle" align="right">0.313 (0.185-0.530)</td>
<td valign="middle" align="right">0.001</td>
<td valign="middle" align="right">0.280 (0.167-0.470)</td>
<td valign="middle" align="right">0.001</td>
<td valign="middle" align="right">0.713 (0.254-2.002)</td>
<td valign="middle" align="right">0.902</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 24 months</td>
<td valign="middle" align="right">0.363 (0.230-0.575)</td>
<td valign="middle" align="right">0.012</td>
<td valign="middle" align="right">0.306 (0.196-0.480)</td>
<td valign="middle" align="right">0.057</td>
<td valign="middle" align="right">0.660 (0.273-1.592)</td>
<td valign="middle" align="right">0.602</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 30 months</td>
<td valign="middle" align="right">0.343 (0.227-0.517)</td>
<td valign="middle" align="right">0.386</td>
<td valign="middle" align="right">0.325 (0.216-0.489)</td>
<td valign="middle" align="right">0.384</td>
<td valign="middle" align="right">0.644 (0.282-1.472)</td>
<td valign="middle" align="right">0.759</td>
</tr>
<tr>
<td valign="middle" align="left">8 days to 3 years</td>
<td valign="middle" align="right">0.376 (0.255-0.554)</td>
<td valign="middle" align="right">0.095</td>
<td valign="middle" align="right">0.333 (0.226-0.489)</td>
<td valign="middle" align="right">0.124</td>
<td valign="middle" align="right">0.658 (0.311-1.394)</td>
<td valign="middle" align="right">0.795</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>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.</p>
<p>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 (<xref ref-type="bibr" rid="B12">12</xref>). The 42-day threshold is also the standard time frame used by the VAERS system for defining vaccine-associated arthritis (<xref ref-type="bibr" rid="B13">13</xref>). 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.</p>
<p>Our findings corroborate prior research suggesting a potential autoimmune link between HPV infection and rheumatic diseases, including RA and systemic lupus erythematosus (<xref ref-type="bibr" rid="B14">14</xref>, <xref ref-type="bibr" rid="B15">15</xref>). One tertiary center study showed that 48.2% of HPV DNA+ patients had autoantibody or autoimmune disease (<xref ref-type="bibr" rid="B16">16</xref>). The mechanism for their condition was generally considered to be via molecular mimicry and cross reactions between self-antigens and viral proteins (<xref ref-type="bibr" rid="B17">17</xref>). HPV-47, in particular, contains sequences with molecular mimicry to citrullinated autoantigens involved in RA pathogenesis (<xref ref-type="bibr" rid="B18">18</xref>). 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 (<xref ref-type="bibr" rid="B19">19</xref>).</p>
<p>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&#x2013;59 years (<xref ref-type="bibr" rid="B20">20</xref>).</p>
<p>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 (<xref ref-type="bibr" rid="B21">21</xref>). 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.</p>
<p>Our results are consistent with the known epidemiology of JIA in the United States. Previously reported incidence rates in girls aged 11&#x2013;15 years range from 13.2 to 30.2 per 100,000 (<xref ref-type="bibr" rid="B22">22</xref>). In our study, the observed JIA incidence&#x2014;based on both diagnostic coding and medication use&#x2014;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&#x2013;12 months, likely reflecting improved specificity for true JIA cases. Taken together, the consistency between our observations and prior evidence&#x2014;including NHANES epidemiologic data, molecular mimicry mechanisms, and the Treg-promoting potential of inactivated vaccines&#x2014;provides additional support for the biological plausibility of a protective association.</p>
<p>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&#x2014;such as family history, HLA genotype, environmental exposures, and socioeconomic status &#x2014;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.</p>
<p>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.</p>
</sec>
</body>
<back>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>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/<xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary Material</bold>
</xref>.</p>
</sec>
<sec id="s6" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>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&#x2019; legal guardians/next of kin in accordance with the national legislation and institutional requirements.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>W-YW: Writing &#x2013; original draft, Conceptualization. H-HC: Conceptualization, Writing &#x2013; review &amp; editing. M-CT: Formal analysis, Writing &#x2013; review &amp; editing. Y-CH: Writing &#x2013; review &amp; editing. J-PC: Data curation, Writing &#x2013; review &amp; editing. L-SF: Conceptualization, Writing &#x2013; review &amp; editing.</p>
</sec>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research and/or publication of this article.</p>
</sec>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec id="s10" sec-type="ai-statement">
<title>Generative AI statement</title>
<p>The author(s) declare that no Generative AI was used in the creation of this manuscript.</p>
<p>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.</p>
</sec>
<sec id="s11" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>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.</p>
</sec>
<sec id="s12" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fimmu.2025.1621939/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fimmu.2025.1621939/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Table1.docx" id="SF1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document">
<label>Supplementary Table&#xa0;1</label>
<caption>
<p>Drugs in MS109.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Table2.docx" id="SF2" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document">
<label>Supplementary Table&#xa0;2</label>
<caption>
<p>Diagnosis of ICD 10=M.08.xx.</p>
</caption>
</supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Franco</surname> <given-names>EL</given-names>
</name>
<name>
<surname>Duarte-Franco</surname> <given-names>E</given-names>
</name>
<name>
<surname>Ferenczy</surname> <given-names>A</given-names>
</name>
</person-group>. <article-title>Cervical cancer: Epidemiology, prevention and the role of human papillomavirus infection</article-title>. <source>Cmaj</source>. (<year>2001</year>) <volume>164</volume>:<page-range>1017&#x2013;25</page-range>. doi: <pub-id pub-id-type="doi">10.1158/1055-9965.EPI-12-1406</pub-id>, PMID: <pub-id pub-id-type="pmid">23549399</pub-id></citation></ref>
<ref id="B2">
<label>2</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pellegrino</surname> <given-names>P</given-names>
</name>
<name>
<surname>Carnovale</surname> <given-names>C</given-names>
</name>
<name>
<surname>Pozzi</surname> <given-names>M</given-names>
</name>
<name>
<surname>Antoniazzi</surname> <given-names>S</given-names>
</name>
<name>
<surname>Perrone</surname> <given-names>V</given-names>
</name>
<name>
<surname>Salvati</surname> <given-names>D</given-names>
</name>
<etal/>
</person-group>. <article-title>On the relationship between human papilloma virus vaccine and autoimmune diseases</article-title>. <source>Autoimmun Rev</source>. (<year>2014</year>) <volume>13</volume>:<page-range>736&#x2013;41</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.autrev.2014.01.054</pub-id>, PMID: <pub-id pub-id-type="pmid">24468416</pub-id></citation></ref>
<ref id="B3">
<label>3</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kamolratanakul</surname> <given-names>S</given-names>
</name>
<name>
<surname>Pitisuttithum</surname> <given-names>P</given-names>
</name>
</person-group>. <article-title>Human papillomavirus vaccine efficacy and effectiveness against cancer</article-title>. <source>Vaccines (Basel)</source>. (<year>2021</year>) <volume>9</volume>(<issue>12</issue>):<fpage>1413</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/vaccines9121413</pub-id>, PMID: <pub-id pub-id-type="pmid">34960159</pub-id></citation></ref>
<ref id="B4">
<label>4</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Giuliano</surname> <given-names>AR</given-names>
</name>
<name>
<surname>Palefsky</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Goldstone</surname> <given-names>S</given-names>
</name>
<name>
<surname>Moreira</surname> <given-names>ED</given-names>
<suffix>Jr.</suffix>
</name>
<name>
<surname>Penny</surname> <given-names>ME</given-names>
</name>
<name>
<surname>Aranda</surname> <given-names>C</given-names>
</name>
<etal/>
</person-group>. <article-title>Efficacy of quadrivalent HPV vaccine against HPV Infection and disease in males</article-title>. <source>N Engl J Med</source>. (<year>2011</year>) <volume>364</volume>:<page-range>401&#x2013;11</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1056/NEJMoa0909537</pub-id>, PMID: <pub-id pub-id-type="pmid">21288094</pub-id></citation></ref>
<ref id="B5">
<label>5</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hussein</surname> <given-names>HM</given-names>
</name>
<name>
<surname>Rahal</surname> <given-names>EA</given-names>
</name>
</person-group>. <article-title>The role of viral infections in the development of autoimmune diseases</article-title>. <source>Crit Rev Microbiol</source>. (<year>2019</year>) <volume>45</volume>:<fpage>394</fpage>&#x2013;<lpage>412</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/1040841X.2019.1614904</pub-id>, PMID: <pub-id pub-id-type="pmid">31145640</pub-id></citation></ref>
<ref id="B6">
<label>6</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rojas</surname> <given-names>M</given-names>
</name>
<name>
<surname>Restrepo-Jim&#xe9;nez</surname> <given-names>P</given-names>
</name>
<name>
<surname>Monsalve</surname> <given-names>DM</given-names>
</name>
<name>
<surname>Pacheco</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Acosta-Ampudia</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Ram&#xed;rez-Santana</surname> <given-names>C</given-names>
</name>
<etal/>
</person-group>. <article-title>Molecular mimicry and autoimmunity</article-title>. <source>J Autoimmun</source>. (<year>2018</year>) <volume>95</volume>:<page-range>100&#x2013;23</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jaut.2018.10.012</pub-id>, PMID: <pub-id pub-id-type="pmid">30509385</pub-id></citation></ref>
<ref id="B7">
<label>7</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sogkas</surname> <given-names>G</given-names>
</name>
<name>
<surname>Atschekzei</surname> <given-names>F</given-names>
</name>
<name>
<surname>Adriawan</surname> <given-names>IR</given-names>
</name>
<name>
<surname>Dubrowinskaja</surname> <given-names>N</given-names>
</name>
<name>
<surname>Witte</surname> <given-names>T</given-names>
</name>
<name>
<surname>Schmidt</surname> <given-names>RE</given-names>
</name>
</person-group>. <article-title>Cellular and molecular mechanisms breaking immune tolerance in inborn errors of immunity</article-title>. <source>Cell Mol Immunol</source>. (<year>2021</year>) <volume>18</volume>:<page-range>1122&#x2013;40</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41423-020-00626-z</pub-id>, PMID: <pub-id pub-id-type="pmid">33795850</pub-id></citation></ref>
<ref id="B8">
<label>8</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Geier</surname> <given-names>DA</given-names>
</name>
<name>
<surname>Geier</surname> <given-names>MR</given-names>
</name>
</person-group>. <article-title>Quadrivalent human papillomavirus vaccine and autoimmune adverse events: A case-control assessment of the vaccine adverse event reporting system (VAERS) database</article-title>. <source>Immunol Res</source>. (<year>2017</year>) <volume>65</volume>:<fpage>46</fpage>&#x2013;<lpage>54</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12026-016-8815-9</pub-id>, PMID: <pub-id pub-id-type="pmid">27406735</pub-id></citation></ref>
<ref id="B9">
<label>9</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maldonado</surname> <given-names>I</given-names>
</name>
<name>
<surname>Rodr&#xed;guez Ni&#xf1;o</surname> <given-names>N</given-names>
</name>
<name>
<surname>Valencia</surname> <given-names>CF</given-names>
</name>
<name>
<surname>Ortiz</surname> <given-names>D</given-names>
</name>
<name>
<surname>Ayala</surname> <given-names>N</given-names>
</name>
<name>
<surname>D&#x2019;Silva</surname> <given-names>L</given-names>
</name>
<etal/>
</person-group>. <article-title>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</article-title>. <source>Vaccine</source>. (<year>2024</year>) <volume>42</volume>:<page-range>2414&#x2013;20</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.vaccine.2024.02.085</pub-id>, PMID: <pub-id pub-id-type="pmid">38458869</pub-id></citation></ref>
<ref id="B10">
<label>10</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chao</surname> <given-names>C</given-names>
</name>
<name>
<surname>Klein</surname> <given-names>NP</given-names>
</name>
<name>
<surname>Velicer</surname> <given-names>CM</given-names>
</name>
<name>
<surname>Sy</surname> <given-names>LS</given-names>
</name>
<name>
<surname>Slezak</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Takhar</surname> <given-names>H</given-names>
</name>
<etal/>
</person-group>. <article-title>Surveillance of autoimmune conditions following routine use of quadrivalent human papillomavirus vaccine</article-title>. <source>J Intern Med</source>. (<year>2012</year>) <volume>271</volume>:<fpage>193</fpage>&#x2013;<lpage>203</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-2796.2011.02467.x</pub-id>, PMID: <pub-id pub-id-type="pmid">21973261</pub-id></citation></ref>
<ref id="B11">
<label>11</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Topaloglu</surname> <given-names>U</given-names>
</name>
<name>
<surname>Palchuk</surname> <given-names>MB</given-names>
</name>
</person-group>. <article-title>Using a federated network of real-world data to optimize clinical trials operations</article-title>. <source>JCO Clin Cancer Inform</source>. (<year>2018</year>) <volume>2</volume>:<fpage>1</fpage>&#x2013;<lpage>10</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1200/CCI.17.00067</pub-id>, PMID: <pub-id pub-id-type="pmid">30652541</pub-id></citation></ref>
<ref id="B12">
<label>12</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gil-Vila</surname> <given-names>A</given-names>
</name>
<name>
<surname>Ravichandran</surname> <given-names>N</given-names>
</name>
<name>
<surname>Selva-O&#x2019;Callaghan</surname> <given-names>A</given-names>
</name>
<name>
<surname>Sen</surname> <given-names>P</given-names>
</name>
<name>
<surname>Nune</surname> <given-names>A</given-names>
</name>
<name>
<surname>Gaur</surname> <given-names>PS</given-names>
</name>
<etal/>
</person-group>. <article-title>COVID-19 Vaccination in Autoimmune Diseases (COVAD) study: Vaccine safety in idiopathic inflammatory myopathies</article-title>. <source>Muscle Nerve</source>. (<year>2022</year>) <volume>66</volume>:<page-range>426&#x2013;37</page-range>. doi: <pub-id pub-id-type="doi">10.1002/mus.27681</pub-id>, PMID: <pub-id pub-id-type="pmid">35869701</pub-id></citation></ref>
<ref id="B13">
<label>13</label>
<citation citation-type="web">
<article-title>VAERS table of reportable events following vaccination</article-title>(<year>2024</year>). Available online at: <uri xlink:href="https://vaers.hhs.gov/docs/VAERS_Table_of_Reportable_Events_Following_Vaccination.pdf">https://vaers.hhs.gov/docs/VAERS_Table_of_Reportable_Events_Following_Vaccination.pdf</uri>. (Accessed May 6, 2024)</citation></ref>
<ref id="B14">
<label>14</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>YT</given-names>
</name>
<name>
<surname>Tsou</surname> <given-names>HK</given-names>
</name>
<name>
<surname>Chiou</surname> <given-names>JY</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>YH</given-names>
</name>
<name>
<surname>Chou</surname> <given-names>MC</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>JC</given-names>
</name>
</person-group>. <article-title>Association of human papillomavirus infection with risk for rheumatoid arthritis: A population-based cohort study</article-title>. <source>Ann Rheum Dis</source>. (<year>2019</year>) <volume>78</volume>:<page-range>1734&#x2013;6</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1136/annrheumdis-2019-215931</pub-id>, PMID: <pub-id pub-id-type="pmid">31413002</pub-id></citation></ref>
<ref id="B15">
<label>15</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname> <given-names>LH</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>JY</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>YZ</given-names>
</name>
<name>
<surname>Chiou</surname> <given-names>JY</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>R</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>JC</given-names>
</name>
</person-group>. <article-title>Risk of systemic lupus erythematosus in patients with human papillomavirus infection: A population-based retrospective cohort study</article-title>. <source>Lupus</source>. (<year>2018</year>) <volume>27</volume>:<page-range>2279&#x2013;83</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1177/0961203318809179</pub-id>, PMID: <pub-id pub-id-type="pmid">30451639</pub-id></citation></ref>
<ref id="B16">
<label>16</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Donmez</surname> <given-names>HG</given-names>
</name>
<name>
<surname>Tanacan</surname> <given-names>A</given-names>
</name>
<name>
<surname>Unal</surname> <given-names>C</given-names>
</name>
<name>
<surname>Fadiloglu</surname> <given-names>E</given-names>
</name>
<name>
<surname>Onder</surname> <given-names>SC</given-names>
</name>
<name>
<surname>Portakal</surname> <given-names>O</given-names>
</name>
<etal/>
</person-group>. <article-title>Human papillomavirus infection and autoimmune disorders: A tertiary center experience</article-title>. <source>Pathog Dis</source>. (<year>2019</year>) <volume>77</volume>(<issue>3</issue>):<fpage>ftz028</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/femspd/ftz028</pub-id>, PMID: <pub-id pub-id-type="pmid">31034015</pub-id></citation></ref>
<ref id="B17">
<label>17</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Smolen</surname> <given-names>JS</given-names>
</name>
<name>
<surname>Aletaha</surname> <given-names>D</given-names>
</name>
<name>
<surname>McInnes</surname> <given-names>IB</given-names>
</name>
</person-group>. <article-title>Rheumatoid arthritis</article-title>. <source>Lancet</source>. (<year>2016</year>) <volume>388</volume>:<page-range>2023&#x2013;38</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0140-6736(16)30173-8</pub-id>, PMID: <pub-id pub-id-type="pmid">27156434</pub-id></citation></ref>
<ref id="B18">
<label>18</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname> <given-names>J</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>X</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>J</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Z</given-names>
</name>
</person-group>. <article-title>Prevalence and significance of antibodies to citrullinated human papilloma virus-47 E2345&#x2013;362 in rheumatoid arthritis</article-title>. <source>J Autoimmun</source>. (<year>2008</year>) <volume>31</volume>:<page-range>131&#x2013;5</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jaut.2008.04.021</pub-id>, PMID: <pub-id pub-id-type="pmid">18511237</pub-id></citation></ref>
<ref id="B19">
<label>19</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xue</surname> <given-names>JS</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>YL</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>C</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>XJ</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>H</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>Y</given-names>
</name>
</person-group>. <article-title>Effects of Th17 cells and IL-17 in the progression of cervical carcinogenesis with high-risk human papillomavirus infection</article-title>. <source>Cancer Med</source>. (<year>2018</year>) <volume>7</volume>:<fpage>297</fpage>&#x2013;<lpage>306</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/cam4.1279</pub-id>, PMID: <pub-id pub-id-type="pmid">29277958</pub-id></citation></ref>
<ref id="B20">
<label>20</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname> <given-names>G</given-names>
</name>
<name>
<surname>Ren</surname> <given-names>Z</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>K</given-names>
</name>
</person-group>. <article-title>Association between human papillomavirus infection or immunization and risk for rheumatoid arthritis</article-title>. <source>Front Immunol</source>. (<year>2023</year>) <volume>14</volume>:<elocation-id>1130217</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fimmu.2023.1130217</pub-id>, PMID: <pub-id pub-id-type="pmid">37122753</pub-id></citation></ref>
<ref id="B21">
<label>21</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Prietl</surname> <given-names>B</given-names>
</name>
<name>
<surname>Treiber</surname> <given-names>G</given-names>
</name>
<name>
<surname>Chaparro</surname> <given-names>N</given-names>
</name>
<name>
<surname>Graninger</surname> <given-names>M</given-names>
</name>
<name>
<surname>Obermayer-Pietsch</surname> <given-names>B</given-names>
</name>
</person-group>. <article-title>Harnessing regulatory T cells through vaccination strategies: a path toward immune tolerance</article-title>. <source>Vaccines (Basel)</source>. (<year>2024</year>) <volume>12</volume>:<elocation-id>992</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/vaccines12030992</pub-id>
</citation></ref>
<ref id="B22">
<label>22</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Harrold</surname> <given-names>LR</given-names>
</name>
<name>
<surname>Salman</surname> <given-names>C</given-names>
</name>
<name>
<surname>Shoor</surname> <given-names>S</given-names>
</name>
<name>
<surname>Curtis</surname> <given-names>JR</given-names>
</name>
<name>
<surname>Asgari</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Gelfand</surname> <given-names>JM</given-names>
</name>
<etal/>
</person-group>. <article-title>Incidence and prevalence of juvenile idiopathic arthritis among children in a managed care population, 1996-2009</article-title>. <source>J Rheumatol</source>. (<year>2013</year>) <volume>40</volume>:<page-range>1218&#x2013;25</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.3899/jrheum.120661</pub-id>, PMID: <pub-id pub-id-type="pmid">23588938</pub-id></citation></ref>
</ref-list>
</back>
</article>