BACKGROUND AND OBJECTIVES:

Human papillomavirus is the most common sexually transmitted infection in the United States and causes certain anogenital and oropharyngeal cancers. The 9-valent human papillomavirus vaccine (9vHPV) provides protection against additional types not included in the quadrivalent vaccine. We conducted near real-time vaccine safety surveillance for 24 months after the vaccine became available in the Vaccine Safety Datalink.

METHODS:

Immunizations and adverse events were extracted weekly from October 2015 to October 2017 from standardized data files for persons 9 to 26 years old at 6 Vaccine Safety Datalink sites. Prespecified adverse events included anaphylaxis, allergic reaction, appendicitis, Guillain-Barré syndrome, chronic inflammatory demyelinating polyneuropathy, injection site reaction, pancreatitis, seizure, stroke, syncope, and venous thromboembolism. The observed and expected numbers of events after 9vHPV were compared weekly by using sequential methods. Both historical and concurrent comparison groups were used to identify statistical signals for adverse events. Unexpected signals were investigated by medical record review and/or additional analyses.

RESULTS:

During 105 weeks of surveillance, 838 991 doses of 9vHPV were administered. We identified unexpected statistical signals for 4 adverse events: appendicitis among boys 9 to 17 years old after dose 3; pancreatitis among men 18 to 26 years old; and allergic reactions among girls 9 to 17 years old and women 18 to 26 years old after dose 2. On further evaluation, which included medical record review, temporal scan analysis, and additional epidemiological analyses, we did not confirm signals for any adverse events.

CONCLUSIONS:

After 2 years of near real-time surveillance of 9vHPV and several prespecified adverse events, no new safety concerns were identified.

What’s Known on This Subject:

Human papillomavirus (HPV) is a common sexually transmitted infection that can cause cancer and other illnesses. Clinical trials of the 9-valent HPV vaccine have demonstrated efficacy and safety, but population-based studies are needed to evaluate its safety profile.

What This Study Adds:

In this postlicensure study, we documented ∼839 000 9-valent HPV doses administered from 2015 to 2017 among persons 9 to 26 years old in the Vaccine Safety Datalink; no new safety concerns were identified. With these findings, we add to the safety profile established in prelicensure clinical trials.

Human papillomavirus (HPV) is the most common sexually transmitted infection in the United States.1  Although most of the 14 million infections that occur each year in the United States are transient, HPV is a known cause of certain anogenital and oropharyngeal cancers and accounts for 3% of all cancers among women and 2% of all cancers among men.1,2  The Advisory Committee on Immunization Practices (ACIP) first recommended the quadrivalent vaccine (4-valent human papillomavirus vaccine [4vHPV], Gardasil; Merck and Co) in 2006 for routine vaccination of 11- to 12-year-old girls with catchup through age 26 years3 ; recommendations for boys followed in 2011.4  In February 2015, the ACIP recommended 3 doses (0, 1–2, and 6 months) of the recently licensed 9-valent human papillomavirus vaccine (9vHPV) (Gardasil 9; Merck and Co) be routinely administered to boys and girls starting at the age of 11 to 12 years. Vaccination was also recommended for females 13 to 26 years and males 13 to 21 years, if not previously vaccinated, and through age 26 for men in specific risk groups.5  The ACIP changed its recommendation in October 2016 from 3 doses to 2 doses (0, 6–12 months) for those who begin the series before the age of 15 years.6,7  In June 2019, the ACIP recommended shared clinical decision-making for HPV vaccination in individuals aged 27 to 45 years but did not extend the recommended catch-up age group beyond age 26 years.8  Compared with the quadrivalent vaccine (4vHPV), 9vHPV has more than twice the amount of an aluminum-based adjuvant and protects against an additional 5 HPV types.5 

Although prelicensure clinical trials of the 9-valent vaccine did not detect any important safety concerns, these studies were not powered to detect uncommon adverse events; population-based monitoring of 9vHPV is warranted.9,10  Several months after the ACIP recommendation for 9vHPV, the Vaccine Safety Datalink (VSD) initiated near real-time surveillance to assess the risks of prespecified adverse events after receipt of 9vHPV.

The VSD was established in 1990 as a collaborative project between several integrated health care organizations and the Centers for Disease Control and Prevention to monitor vaccine safety.11  Six integrated health care organizations in the VSD contributed data for this surveillance: Kaiser Permanente Northern California, Southern California, Colorado, Oregon, and Washington and Marshfield Clinic (Marshfield, Wisconsin). The study population consisted of persons 9 to 26 years old during the study period (October 4, 2015 through October 7, 2017). Standardized electronic data files were generated weekly and contained information on demographics, vaccinations, and medical encounters. Data were aggregated to create prospective cohorts that were followed for up to 180 days from the date of vaccination.

This study was approved by institutional review boards at all participating sites.

Adverse events were prespecified on the basis of reports from clinical trials, the Vaccine Adverse Event Report System, or other published investigations including a VSD safety study of 4vHPV.12  Adverse events were identified by using International Classification of Diseases, 10th Revision (ICD-10) codes assigned during inpatient, outpatient, or emergency department (ED) encounters (Supplemental Table 4); only diagnoses assigned in event-specific, postvaccination exposure windows and settings were included in the analysis (Table 1). A 10-week lag was applied before each weekly analysis to permit administrative corrections to the electronic data, enhance data completeness, and ensure that sufficient time had passed to cover postvaccination risk windows.13 

TABLE 1

Characteristics of Prespecified Adverse Events

Medical SettingPostvaccination Risk Window, dFirst Episode in What Period?Primary Comparison Groupa
Uncommon or rare adverse event     
 Anaphylaxis ED, inpatient 0–2 First in 42 d Historic 
 Appendicitis ED, inpatient 1–42 First in 42 d Historic 
 Pancreatitis ED, inpatient 1–42 First ever Historic 
 GBS Outpatient, ED, inpatient 1–42 First in 42 d Historic 
 CIDP Outpatient, ED, inpatient 1–180 First ever Historic 
 Seizure ED, inpatient 0–42 First in 42 d and first ever Concurrent 
 Stroke ED, inpatient 0–42 First in 42 d Historic 
 VTE Outpatient, ED, inpatient 1–42 First in 1 y Historic 
Common adverse event     
 Allergic reactions Outpatient, ED, inpatient 0–2 for ED and inpatient,1–2 for outpatient First in 42 d Concurrent 
 Injection site reactions Outpatient, ED, inpatient 1–6 First in 42 d Concurrent 
 Nonspecific reactions Outpatient, ED, inpatient 0–6 First in 42 d Concurrent 
 Syncope Outpatient, ED, inpatient Day 0 First in 2 d Concurrent 
Medical SettingPostvaccination Risk Window, dFirst Episode in What Period?Primary Comparison Groupa
Uncommon or rare adverse event     
 Anaphylaxis ED, inpatient 0–2 First in 42 d Historic 
 Appendicitis ED, inpatient 1–42 First in 42 d Historic 
 Pancreatitis ED, inpatient 1–42 First ever Historic 
 GBS Outpatient, ED, inpatient 1–42 First in 42 d Historic 
 CIDP Outpatient, ED, inpatient 1–180 First ever Historic 
 Seizure ED, inpatient 0–42 First in 42 d and first ever Concurrent 
 Stroke ED, inpatient 0–42 First in 42 d Historic 
 VTE Outpatient, ED, inpatient 1–42 First in 1 y Historic 
Common adverse event     
 Allergic reactions Outpatient, ED, inpatient 0–2 for ED and inpatient,1–2 for outpatient First in 42 d Concurrent 
 Injection site reactions Outpatient, ED, inpatient 1–6 First in 42 d Concurrent 
 Nonspecific reactions Outpatient, ED, inpatient 0–6 First in 42 d Concurrent 
 Syncope Outpatient, ED, inpatient Day 0 First in 2 d Concurrent 
a

Concurrent comparison group: visits during which comparator vaccines were administered, with analytic strata defined by age (in 1-y increments), site, sex, and week of the vaccination visit. Two historic comparison groups (2007–2014) were used: general VSD population and vaccinated VSD population, vaccinated with comparator vaccines.

The adverse events included both common and uncommon diseases. Primary analyses for uncommon diseases used large historical comparison groups, whereas common diseases were analyzed with a smaller concurrent comparison group; both are described below. To be consistent with the earlier study of 4vHPV,12  we considered the following adverse events to be uncommon in the age groups of interest: anaphylaxis, appendicitis, Guillain-Barré syndrome (GBS), pancreatitis, seizures, stroke, and venous thromboembolism (VTE). Chronic inflammatory demyelinating polyneuropathy (CIDP), not evaluated previously, was also considered an uncommon event. Injection site reactions, allergic reactions (evaluated separately in the outpatient, ED, and inpatient settings), syncope, and nonspecific reactions (eg, adverse effect of viral vaccines) were considered common (Table 1).14 

Individuals were considered exposed if electronic medical records indicated receipt of 9vHPV during the study period. Doses given within 42 days of a previous dose in the same person were excluded from the analysis to prevent overlapping risk windows.

Comparison groups were age-comparable persons observed in either historical or concurrent periods of time relative to 9vHPV uptake who were unexposed to 9vHPV vaccine. Person-time rates based on the general population provide increased stability for rare adverse events but may introduce bias because of systematic differences in exposed and unexposed groups. In contrast, data based on comparator vaccines, both historical and concurrent, were derived from smaller populations with risk estimates that are less stable, but exposed and unexposed groups are more likely to be similar. For the analysis of common adverse events, we used a concurrent comparison group. For uncommon events, the primary analytic method used 2 complementary historical comparison groups, but we also evaluated uncommon events (eg, pancreatitis and seizures) using the concurrent comparison group because the analytic infrastructure for uncommon events was under development at the start of the study.

Data for both historical comparison groups were derived from persons 9 to 26 years of age from 2007 through 2014. One group was the general VSD population that was used to estimate (via modeling) general background person-time rates by sex, site, and single-year of age; there was no exclusion based on 4vHPV vaccination. The rates were multiplied by the observed number of 9vHPV doses to produce expected counts, which were prorated to the length of the postvaccination risk window (Table 1). For informational and comparison purposes, the rates summarized across subgroups (but not prorated to the length of the postvaccination risk window) are provided in Supplemental Table 5. The other historical group consisted of persons with visits at which comparator vaccines routinely given to this age group (tetanus, diphtheria; tetanus, diphtheria, acellular pertussis; meningococcal conjugate; hepatitis A; varicella) were administered. Visits in which 4vHPV was administered were not included in this group. The number of events observed in postvaccination risk windows for comparator vaccine visits and the number of comparator vaccine doses administered in the historical period were incorporated into the analyses.

The concurrent comparison group was defined in the same manner as the historical comparator vaccine group, except that the visits among the VSD population occurred during the study period. Across analytic strata defined by age, site, sex, and week of the vaccination visit, exposed versus unexposed comparisons were performed.

We used the Rapid Cycle Analysis (RCA) methodology and near real-time data to compare adverse event rates in a recently vaccinated group with rates from an unvaccinated group.15  The VSD has investigated multiple vaccines using the RCA approach, including 4vHPV.12,1521 

To estimate the relative risk (RR) for prespecified adverse events, we used the Poisson-based maximized sequential probability ratio test (MaxSPRT) for analyses in which the historical general VSD population comparison group was used,22  the conditional maximized sequential probability ratio test (CMaxSPRT) for analyses in which the historical vaccinated VSD population was used,23  and the exact sequential analysis (ESA) for analyses in which the concurrent comparison group was used.12,24  We conducted the analyses in the overall study population as well as in subgroups defined by age group and sex. Analyses for more-common outcomes (eg, syncope) were dose specific, whereas others (eg, GBS) pertain to 1 or more doses (ie, any dose). Because data were analyzed on a weekly basis, sequential methods were used to maintain an overall 1-sided type I error rate of 0.05 across the multiple tests performed for each adverse event, subgroup, and statistical method combination. Sequential methods periodically compare the number of cases of each adverse event among the exposed with the number of cases observed or expected among a comparison group unexposed to 9vHPV.12,2224  A preliminary statistical signal was generated if a test statistic from an analysis exceeded a predetermined threshold. A more-detailed description of the analytical methods used in this study is provided in the Supplemental Information.

We conducted the first analytic run for ESA in week 25, which encompassed data from October 4, 2015 (week 1), to March 20, 2016. We conducted the first analytic run for MaxSPRT and CMaxSPRT in week 53, which encompassed data from October 4, 2015, to October 2, 2016. All subsequent analyses were conducted at weekly intervals. For reporting purposes, signals detected on the first analytical runs were arbitrarily designated as having occurred in weeks 25 (ESA) and 53 (MaxSPRT and CMaxSPRT), respectively. The analyses were delayed to allow for development of the analytic infrastructure.

In general, follow-up of a signal entailed 1 or more additional evaluations to determine if the signal was an indicator of increased risk. These included data quality assessments, use of temporal scan statistics,25  medical record reviews to validate the adverse event, and a self-controlled risk interval (SCRI) analysis.26  No follow-up investigations were conducted for syncope, injection site reactions, and nonspecific reactions because associations were expected based on 9vHPV clinical trials27  and clinical or observational studies of 4vHPV.28,29 

All analyses were conducted by using SAS version 9.4 (SAS Institute, Inc, Cary, NC).

Over the 2-year study period, we observed 838 991 doses of 9vHPV vaccine administered among the VSD population. Of 638 947 (76.2%) doses administered among individuals aged 9 to 17 years, 47.6% of doses were administered among girls. Of 200 044 (23.8%) doses administered among persons aged 18 to 26 years, 64.4% of doses were administered among women. Doses administered to 9- to 17-year-olds increased sharply during August of each year; the counts for 18- to 26-year-olds are relatively constant during the study period (Fig 1).

FIGURE 1

Histogram of weekly counts of 9vHPV vaccinations administered in the VSD population from the week starting on October 4, 2015, through the week starting on October 1, 2017, by age and sex.

FIGURE 1

Histogram of weekly counts of 9vHPV vaccinations administered in the VSD population from the week starting on October 4, 2015, through the week starting on October 1, 2017, by age and sex.

Close modal

No signals were observed in the MaxSPRT or CMaxSPRT analyses for anaphylaxis, appendicitis, CIDP, GBS, seizure, stroke, or VTE (Tables 1, 2, and 3, and Suplemental Tables 4 and 5). In week 71, the MaxSPRT analysis identified a statistical signal for pancreatitis among men 18 to 26 years old after any dose. There were 8 exposed cases of pancreatitis in this subgroup; the RR was 3.1 (P < .05, Table 2). No pancreatitis statistical signal was observed in the CMaxSPRT analysis (RR = 1.9, P > .05). Pancreatitis was also evaluated using ESA; the RR was elevated but not statistically significant (RR = 4.7, P = .47).

TABLE 2

Summary of Uncommon Adverse Events in Selected Subgroups Evaluated With Historical Comparison Groups Using MaxSPRT and CMaxSPRT

AEType of SPRT AnalysisaSubgroup, Sex, Age Group in yNo. Observed VaccinationsNo. Observed AENo. Comparator Vaccinations, Historical PeriodNo. AEs in the Historical PeriodNo. Expected AEsRRTest StatisticCritical ValuebSignal
Anaphylaxis Maximized Male, 9–17 334 381 — — 0.9 2.7 No 
 CMax Male, 9–17 334 381 1 053 642 — 2.6 No 
Appendicitis Maximized Male, 9–17 78 885 33 — — 25.6 1.3 0.97 3.5 No 
 CMax Male, 9–17 78 885 33 1 053 642 312 — 1.4 1.6 3.8 No 
GBS Maximized Female, 18–26 128 645 — — 0.6 2.6 No 
 CMax Female, 18–26 128 645 431 401 — 2.7 No 
Pancreatitis Maximized Male, 18–26 51 944 — — 2.6 3.1 3.7 2.9 Yes 
 CMax Male, 18–26 51 944 349 966 29 — 1.9 1.1 3.4 No 
Seizures Maximized Female, 9–17 304 384 44 — — 115.2 0.4 3.8 No 
 CMax Female, 9–17 304 384 44 698 263 105 — 1.0 3.7 No 
Stroke Maximized Female, 18–26 128 645 — — 4.9 1.0 2.9 No 
 CMax Female, 18–26 128 645 431 401 11 — 1.5 0.3 2.9 No 
VTE Maximized Female, 18–26 128 645 — — 9.8 0.8 3.1 No 
 CMax Female, 18–26 128 645 431 401 54 — 0.5 3.4 No 
AEType of SPRT AnalysisaSubgroup, Sex, Age Group in yNo. Observed VaccinationsNo. Observed AENo. Comparator Vaccinations, Historical PeriodNo. AEs in the Historical PeriodNo. Expected AEsRRTest StatisticCritical ValuebSignal
Anaphylaxis Maximized Male, 9–17 334 381 — — 0.9 2.7 No 
 CMax Male, 9–17 334 381 1 053 642 — 2.6 No 
Appendicitis Maximized Male, 9–17 78 885 33 — — 25.6 1.3 0.97 3.5 No 
 CMax Male, 9–17 78 885 33 1 053 642 312 — 1.4 1.6 3.8 No 
GBS Maximized Female, 18–26 128 645 — — 0.6 2.6 No 
 CMax Female, 18–26 128 645 431 401 — 2.7 No 
Pancreatitis Maximized Male, 18–26 51 944 — — 2.6 3.1 3.7 2.9 Yes 
 CMax Male, 18–26 51 944 349 966 29 — 1.9 1.1 3.4 No 
Seizures Maximized Female, 9–17 304 384 44 — — 115.2 0.4 3.8 No 
 CMax Female, 9–17 304 384 44 698 263 105 — 1.0 3.7 No 
Stroke Maximized Female, 18–26 128 645 — — 4.9 1.0 2.9 No 
 CMax Female, 18–26 128 645 431 401 11 — 1.5 0.3 2.9 No 
VTE Maximized Female, 18–26 128 645 — — 9.8 0.8 3.1 No 
 CMax Female, 18–26 128 645 431 401 54 — 0.5 3.4 No 

Results were extracted from the final report (week starting October 1, 2017) for all events except for pancreatitis, which comes from the week when it first signaled (week 71). CIDP was not included because there were no events during the study period. AE, adverse event; CMax, conditional maximized; SPRT, sequential probability ratio test; —, not applicable.

a

MaxSPRT: analyses were conducted in combination with the general VSD population historic comparison group (2007–2014). CMaxSPRT: analyses were conducted in combination with the historic VSD population vaccinated with comparator vaccines (2007–2014).

b

Critical values are threshold values of the test statistic above, in which the null hypothesis would be rejected.

Medical record review determined that 7 of the 8 pancreatitis cases were either not incident (n = 2) or were attributed to causes other than vaccination (n = 5). The one confirmed case had no known risk factors for pancreatic disease. Consequently, the pancreatitis statistical signal was classified as a false-positive.

A total of 12 statistical signals were identified for 5 types of adverse events in the ESA analysis: appendicitis, allergic reaction, injection site reaction, syncope, and nonspecific reactions (Table 3). No signals were detected for anaphylaxis, CIDP, GBS, pancreatitis, seizure, stroke, or VTE.

TABLE 3

Summary of Common Adverse Events Using ESA in Subgroups Evaluated With the Concurrent Comparison Group

Adverse EventSubgroup, Sex, Age Group in y, 9vHPV DoseWeek When First SignaledNo. 9vHPV Vaccinations in SubgroupNo. Total CasesaNo. Exposed CasesaRRP
Allergic reactionb Female, 9–17, ED or inpatient, any 82 242 726 33 26 2.7 .04 
 Female, 9–17, outpatient, any No signalc 242 726 86 60 0.8 .75 
 Female, 9–17, ED or inpatient, dose 1 94 109 896 26 17 2.8 .04 
 Female, 9–17, outpatient, dose 1 No signald 109 896 82 50 1.2 .28 
 Female, 18–26, outpatient, dose 2 86 33 118 38 15 1.9 .04 
 Female, 18–26, ED or inpatient, dose 2 No signale 33 118 0.4 .92 
Appendicitisf Male, 9–17, dose 3 84 73 122 50 30 2.1 .03 
 Male, 9–17, any No signalg 271 679 103 81 1.5 .09 
 Male, 9–17, dose 1 No signalg 106 741 47 25 1.4 .23 
 Male, 9–17, dose 2 No signalg 91 156 47 26 1.5 .23 
Injection site reaction Male, 9–17, dose 3 26 23 409 29 18 2.5 .03 
Nonspecific reaction Male, 18–26, dose 3 25 4054 95.0 .04 
 Male, 18–26, dose 1 34 13 228 14 11.1 .04 
 Female, 18–26, dose 1 50 26 711 71 34 1.7 .03 
 Female, 18–26, any 105 128 806 215 126 1.3 .04 
Syncope Female, 18–26, any 25 28 234 98 67 1.8 .007 
 Female, 18–26, dose 1 25 12 245 65 35 2.0 .004 
 Female, 18–26, dose 2 31 10 924 60 25 1.7 .04 
Adverse EventSubgroup, Sex, Age Group in y, 9vHPV DoseWeek When First SignaledNo. 9vHPV Vaccinations in SubgroupNo. Total CasesaNo. Exposed CasesaRRP
Allergic reactionb Female, 9–17, ED or inpatient, any 82 242 726 33 26 2.7 .04 
 Female, 9–17, outpatient, any No signalc 242 726 86 60 0.8 .75 
 Female, 9–17, ED or inpatient, dose 1 94 109 896 26 17 2.8 .04 
 Female, 9–17, outpatient, dose 1 No signald 109 896 82 50 1.2 .28 
 Female, 18–26, outpatient, dose 2 86 33 118 38 15 1.9 .04 
 Female, 18–26, ED or inpatient, dose 2 No signale 33 118 0.4 .92 
Appendicitisf Male, 9–17, dose 3 84 73 122 50 30 2.1 .03 
 Male, 9–17, any No signalg 271 679 103 81 1.5 .09 
 Male, 9–17, dose 1 No signalg 106 741 47 25 1.4 .23 
 Male, 9–17, dose 2 No signalg 91 156 47 26 1.5 .23 
Injection site reaction Male, 9–17, dose 3 26 23 409 29 18 2.5 .03 
Nonspecific reaction Male, 18–26, dose 3 25 4054 95.0 .04 
 Male, 18–26, dose 1 34 13 228 14 11.1 .04 
 Female, 18–26, dose 1 50 26 711 71 34 1.7 .03 
 Female, 18–26, any 105 128 806 215 126 1.3 .04 
Syncope Female, 18–26, any 25 28 234 98 67 1.8 .007 
 Female, 18–26, dose 1 25 12 245 65 35 2.0 .004 
 Female, 18–26, dose 2 31 10 924 60 25 1.7 .04 

Results were extracted from the report for the week when the adverse event first signaled. NS, no signal.

a

Cases in a specific subgroup are only counted for analytic strata with ≥1 case (either exposed or not), ≥1 9vHPV vaccine, and ≥1 comparator vaccine, in which analytic strata are defined by age (in 1-y increments), site, sex, and week of the vaccination visit.

b

Diagnoses were made in the ED or inpatient setting or in the outpatient setting.

c

No signal was detected for this subgroup. Data were extracted from the report for the week when allergic reaction signaled for girls 9 to 17 y old with any dose in the ED or inpatient setting.

d

No signal was detected for this subgroup. Data were extracted from the report for the week when allergic reaction signaled for girls 9 to 17 y old with dose 1 in the ED or inpatient setting.

e

No signal was detected for this subgroup. Data were extracted from the report for the week when allergic reaction signaled for women 18 to 26 y old with dose 2 in the outpatient setting.

f

Appendicitis was classified as an uncommon adverse event in this study, but a statistical signal was detected with ESA, and it is therefore included here.

g

No signal was detected for this subgroup. Data were extracted from the report for the week when appendicitis signaled for boys 9 to 17 y old with dose 3.

Appendicitis

A statistical signal for appendicitis was observed in week 84 among boys 9 to 17 years old after the third dose of 9vHPV (RR = 2.1, P = .03). Medical record reviews were performed for patients vaccinated with 9vHPV, and all 30 were confirmed to be acute appendicitis within 42 days postvaccination. A temporal scan analysis revealed no statistically significant clustering within the 42-day postvaccination risk interval (Fig 2); the P value for the various scan widths ranged from .78 to .98. RRs for appendicitis ranged from 1.4 to 1.5 among other dose-specific subgroups, but none were statistically significant (Table 3). There were no statistical signals for appendicitis in either the MaxSPRT or CMaxSPRT analyses.

FIGURE 2

Distribution of days to onset of appendicitis in the 42-day risk window after administration of 9vHPV vaccine among boys 9 to 17 years old.

FIGURE 2

Distribution of days to onset of appendicitis in the 42-day risk window after administration of 9vHPV vaccine among boys 9 to 17 years old.

Close modal

SCRI analysis of the appendicitis signal was conducted to further assess its validity. The 1- to 42-day risk interval was compared with the control interval of 43 to 84 days postvaccination within the same person. Of the 30 cases that were identified and reviewed from the control interval, 2 cases were not confirmed and 4 were reclassified because onset occurred within the 1- to 42-day risk interval. The rate ratio of appendicitis after 9vHPV was 1.4 (95% confidence interval = 0.8–2.6). The appendicitis statistical signal was classified as a false-positive because 2 of 3 sequential analytic methods did not signal, there was no temporal clustering, and SCRI analysis failed to confirm the association.

Allergic Reaction

We separately assessed allergic reactions occurring in the outpatient setting and in ED and inpatient settings. The first of 3 statistical signals was observed in the ED or inpatient setting after any 9vHPV dose among girls 9 to 17 years old (RR = 2.7, P = .04). Medical record review identified a possible vaccine-related allergic reaction in 8 (31%) of 26 cases. The most common reasons for nonconfirmation of the diagnosis included injection site reaction (n = 7) and miscoding (n = 4). There was no analogous signal in the outpatient setting (RR = 0.85, P = .75).

The second statistical signal for allergic reactions followed the second dose of 9vHPV among women 18 to 26 years old in the outpatient setting (RR = 1.9, P = .04). Medical records were reviewed for 14 of 15 vaccinated cases, and 6 (43%) were determined as possibly vaccine related. There was no analogous signal for allergic reaction in the ED or inpatient setting (RR = 0.4, P = .75).

The third allergic reaction statistical signal was in the ED or inpatient setting after the first dose of 9vHPV among girls 9 to 17 years old (RR = 2.8, P = .04). This signal included 17 vaccinated cases. Medical record reviews were not conducted for this subgroup given the low predictive value of allergic reaction diagnosis codes demonstrated in the 2 previous allergic reaction statistical signals, including the first signal, which was also in the ED or inpatient setting and in the same age and sex group. There was no analogous signal in the outpatient setting (RR = 1.2, P = .28).

The allergic reaction statistical signals were classified as false-positives because signals were only observed in ESA analyses; discordant results were observed in outpatient, ED, and inpatient settings; and most cases were not confirmed by record review.

Syncope, Injection Site Reaction, and Nonspecific Reaction

There was a signal for syncope in multiple subgroups of women 18 to 26 years old; the RRs were ≤2.0 in each of these subgroups (Table 3). There was one statistical signal for injection site reaction among boys 9 to 17 years old; the RR was 2.5 (P = .03). There was a signal for nonspecific reactions among several subgroups. Men 18 to 26 years old signaled after dose 3; the RR was 95.0, but there were only 3 cases, 2 of which were exposed. The risk estimates dropped sharply over the subsequent weeks, reaching 2.5 with 6 exposed cases in week 79. The statistical signal among men 18 to 26 years old after dose 1 behaved similarly. No follow-up investigations were conducted for syncope, injection site reaction, or nonspecific reactions because they were expected based on clinical trials of 9vHPV, clinical experience with 4vHPV,27,30  and because the diagnoses were unlikely to indicate a serious adverse event.

After 2 years of surveillance and nearly 839 000 administered doses of 9vHPV in the VSD population, we did not identify any new safety concerns from a group of prespecified adverse events. RCA methodology allows rapid, near real-time assessment to identify potential safety concerns, but it is based on unconfirmed electronic diagnosis coded outcomes. In this framework, false-positive statistical signals are expected and further investigation is conducted to determine if a statistical signal represents a valid safety concern. During the surveillance period, there were several statistical signals, but they were either expected on the basis of prelicensure trials or classified as false-positives after further investigation.

Although pancreatitis has not been identified as a safety concern in any prelicensure studies of 9vHPV, we included it as an adverse event because a temporal association between pancreatitis and 4vHPV was described in 2 case reports.31,32  Additionally, in a Vaccine Adverse Event Report System study, researchers reported 9 cases after 4vHPV vaccination,33  although the postvaccination reporting rate was not greater than expected. In our analyses, the pancreatitis statistical signal was not confirmed after further investigation.

Appendicitis was the most frequent serious adverse event (excluding fetal loss) in a pivotal prelicensure trial of 9vHPV, although none of the cases were vaccine related.9,34  Our results are similar to the VSD study of 4vHPV, which also identified a signal for appendicitis, but a causal association was judged unlikely after further analysis.12 

Preliminary statistical signals for allergic reactions were detected in 3 subgroups. However, we concluded that these signals represented false-positives because medical record reviews failed to confirm most cases. Our findings are consistent with prelicensure clinical trials in which serious allergic reactions were rare.35 

We detected preliminary statistical signals for syncope in 18- to 26-year-old women, but not in younger girls, a group in which higher rates of syncope have been reported.36  Our results are in general agreement with both cohort and passive surveillance studies in which researchers have found associations between 4vHPV and vasovagal syncope, particularly among adolescents.29,33  In contrast, the previous VSD study of 4vHPV, in which researchers used analytic methods and comparison groups similar to ours, did not find an association.12 

In previous VSD studies, researchers have assessed the risk of VTE after the 4vHPV vaccine. In the RCA analysis, researchers found a nonstatistically significant elevated risk of VTE (RR = 1.98) among girls 9 to 17 years old after receipt of 4vHPV.12  In a follow-up study, researchers determined that VTE risk among 9- to 26-year-old males and females in the VSD population was not elevated after 4vHPV exposure.37  In our analysis of 9vHPV, we observed 4 VTE cases among girls 9 to 17 years old but no signal.

Our study is subject to a number of limitations. Presumptive cases were identified by using coded diagnoses, but the validity of electronic diagnosis codes varies substantially38 ; diagnoses were validated by medical record review for some statistical signals but not for others. Although sequential methods accounted for weekly testing within a subgroup, they did not account for the number of tests performed each week across subgroups (ie, examining many subgroups increases the likelihood of a false-positive signal). We chose to enhance the sensitivity of our analyses with the understanding that it would require additional investigation to rule out false-positives. Finally, despite the large size of the VSD population, this analysis had limited power to detect signals for rare adverse events, such as GBS.

With this large observational study, we contribute reassuring postlicensure data that will help bolster the safety profile of 9vHPV. We documented nearly 839 000 9vHPV doses administered over 2 years and did not identify any new safety concerns. Although we detected several unexpected potential safety signals, none were confirmed after further evaluation. Our findings are consistent with prelicensure clinical trials, which have determined that 9vHPV, similar to 4vHPV, has a favorable safety profile.9,10,35 

We thank the following individuals for help with this research: Pat Ross and Margarita Magallon (Kaiser Permanente Northern California); Erika L. Kiniry and Jennifer M. Covey (Kaiser Permanente Washington); Kristin Goddard, Stacy Harsh, Stephanie Irving, and Mara Kalter (Kaiser Permanente Northwest); Jo Ann Shoup and Kate Burniece (Kaiser Permanente Colorado); Cheryl Carlson, Denison Ryan, Lina Sy, Claire Jang, Gina Lee, Joy Gelfond, Kerresa Morrissette, Lee Tillman, Nancy Canul, and Sean Anand (Kaiser Permanente Southern California); and Deanna Cole, Tara Johnson, Diane Kohnhorst, Madalyn Palmquist, Suellyn Murray, Becky Pilsner, Carla Rottscheit, Erica Scotty, Samantha Smith, Sandy Strey, and Jane Wesely (Marshfield Clinic Research Institute).

Disclaimer: The findings and conclusion in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention or the US Food and Drug Administration.

Dr Donahue designed the study, provided analytic support, drafted the initial manuscript, and revised subsequent draft manuscripts; Mr Kieke designed the study, conducted the statistical analysis, and reviewed and revised the manuscript; Mr Lewis, Mr Weintraub, and Dr McClure designed the study, provided methodologic and analytic support, and reviewed and revised the manuscript; Ms Vickers and Ms Hanson designed chart abstraction tools, had oversight of data collection, and reviewed and revised the manuscript; Drs Daley, Hechter, Jackson, Klein, Naleway, and Nelson had oversight of study activities at their sites, provided study design and analytic support, and reviewed and revised the manuscript; Ms Gee and Dr DeStefano assisted with study design, interpretation of results, and reviewed and revised the manuscript; Dr Belongia had oversight of all study activities, assisted with study design and interpretation of results, and reviewed and revised the manuscript; and all authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

Results from this study have been presented at the meeting of the Advisory Committee on Immunization Practices (February 21–22, 2018; Atlanta, GA) and at the 2018 Annual Conference on Vaccinology Research (April 23–25, 2018; Bethesda, MD; National Foundation for Infectious Diseases).

FUNDING: Supported by the US Centers for Disease Control and Prevention (contract number 200-2012-53587).

COMPANION PAPER: Companions to this article can be found online at www.pediatrics.org/cgi/doi/10.1542/peds.2019-1791 and www.pediatrics.org/cgi/doi/10.1542/peds.2019-2345.

ACIP

Advisory Committee on Immunization Practices

CIDP

chronic inflammatory demyelinating polyneuropathy

CMaxSPRT

conditional maximized sequential probability ratio test

ED

emergency department

ESA

exact sequential analysis

GBS

Guillain-Barré Syndrome

HPV

human papillomavirus

ICD-10

International Classification of Diseases, 10th Revision

MaxSPRT

maximized sequential probability ratio test

RCA

Rapid Cycle Analysis

RR

relative risk

SCRI

self-controlled risk interval

VSD

Vaccine Safety Datalink

VTE

venous thromboembolism

4vHPV

4-valent human papillomavirus vaccine

9vHPV

9-valent human papillomavirus vaccine

1
Markowitz
LE
,
Dunne
EF
,
Saraiya
M
, et al;
Centers for Disease Control and Prevention (CDC)
.
Human papillomavirus vaccination: recommendations of the Advisory Committee on Immunization Practices (ACIP)
.
MMWR Recomm Rep
.
2014
;
63
(
RR
):
1
30
2
Jemal
A
,
Simard
EP
,
Dorell
C
, et al
.
Annual Report to the Nation on the Status of Cancer, 1975-2009, featuring the burden and trends in human papillomavirus (HPV)-associated cancers and HPV vaccination coverage levels
.
J Natl Cancer Inst
.
2013
;
105
(
3
):
175
201
3
Markowitz
LE
,
Dunne
EF
,
Saraiya
M
, et al;
Centers for Disease Control and Prevention (CDC)
;
Advisory Committee on Immunization Practices (ACIP)
.
Quadrivalent human papillomavirus vaccine: recommendations of the Advisory Committee on Immunization Practices (ACIP)
.
MMWR Recomm Rep
.
2007
;
56
(
RR
):
1
24
4
Centers for Disease Control and Prevention (CDC)
.
Recommendations on the use of quadrivalent human papillomavirus vaccine in males–Advisory Committee on Immunization Practices (ACIP), 2011
.
MMWR Morb Mortal Wkly Rep
.
2011
;
60
(
50
):
1705
1708
5
Petrosky
E
,
Bocchini
JA
 Jr
,
Hariri
S
, et al;
Centers for Disease Control and Prevention (CDC)
.
Use of 9-valent human papillomavirus (HPV) vaccine: updated HPV vaccination recommendations of the advisory committee on immunization practices
.
MMWR Morb Mortal Wkly Rep
.
2015
;
64
(
11
):
300
304
6
Iversen
OE
,
Miranda
MJ
,
Ulied
A
, et al
.
Immunogenicity of the 9-valent HPV vaccine using 2-dose regimens in girls and boys vs a 3-dose regimen in women
.
JAMA
.
2016
;
316
(
22
):
2411
2421
7
Meites
E
,
Kempe
A
,
Markowitz
LE
.
Use of a 2-dose schedule for human papillomavirus vaccination - updated recommendations of the advisory committee on immunization practices
.
MMWR Morb Mortal Wkly Rep
.
2016
;
65
(
49
):
1405
1408
8
Meites
E
,
Szilagyi
PG
,
Chesson
HW
, et al
.
Human papillomavirus vaccination for adults: Updated recommendations of the Advisory Committee on Immunization Practices
.
MMWR
.
2019
;
68
(
32
):
698
702
9
Huh
WK
,
Joura
EA
,
Giuliano
AR
, et al
.
Final efficacy, immunogenicity, and safety analyses of a nine-valent human papillomavirus vaccine in women aged 16-26 years: a randomised, double-blind trial
.
Lancet
.
2017
;
390
(
10108
):
2143
2159
10
Garland
SM
,
Cheung
TH
,
McNeill
S
, et al
.
Safety and immunogenicity of a 9-valent HPV vaccine in females 12-26 years of age who previously received the quadrivalent HPV vaccine
.
Vaccine
.
2015
;
33
(
48
):
6855
6864
11
Baggs
J
,
Gee
J
,
Lewis
E
, et al
.
The Vaccine Safety Datalink: a model for monitoring immunization safety
.
Pediatrics
.
2011
;
127
(
suppl 1
):
S45
S53
12
Gee
J
,
Naleway
A
,
Shui
I
, et al
.
Monitoring the safety of quadrivalent human papillomavirus vaccine: findings from the Vaccine Safety Datalink
.
Vaccine
.
2011
;
29
(
46
):
8279
8284
13
Greene
SK
,
Kulldorff
M
,
Yin
R
, et al
.
Near real-time vaccine safety surveillance with partially accrued data
.
Pharmacoepidemiol Drug Saf
.
2011
;
20
(
6
):
583
590
14
Jackson
LA
,
Yu
O
,
Belongia
EA
, et al
.
Frequency of medically attended adverse events following tetanus and diphtheria toxoid vaccine in adolescents and young adults: a Vaccine Safety Datalink study
.
BMC Infect Dis
.
2009
;
9
:
165
15
Lieu
TA
,
Kulldorff
M
,
Davis
RL
, et al;
Vaccine Safety Datalink Rapid Cycle Analysis Team
.
Real-time vaccine safety surveillance for the early detection of adverse events
.
Med Care
.
2007
;
45
(
10
suppl 2
):
S89
S95
16
Weintraub
ES
,
Baggs
J
,
Duffy
J
, et al
.
Risk of intussusception after monovalent rotavirus vaccination
.
N Engl J Med
.
2014
;
370
(
6
):
513
519
17
Lee
GM
,
Greene
SK
,
Weintraub
ES
, et al;
Vaccine Safety Datalink Project
.
H1N1 and seasonal influenza vaccine safety in the vaccine safety datalink project
.
Am J Prev Med
.
2011
;
41
(
2
):
121
128
18
Klein
NP
,
Fireman
B
,
Yih
WK
, et al;
Vaccine Safety Datalink
.
Measles-mumps-rubella-varicella combination vaccine and the risk of febrile seizures
.
Pediatrics
.
2010
;
126
(
1
). Available at: www.pediatrics.org/cgi/content/full/126/1/e1
19
Nelson
JC
,
Yu
O
,
Dominguez-Islas
CP
, et al
.
Adapting group sequential methods to observational postlicensure vaccine safety surveillance: results of a pentavalent combination DTaP-IPV-Hib vaccine safety study
.
Am J Epidemiol
.
2013
;
177
(
2
):
131
141
20
Daley
MF
,
Yih
WK
,
Glanz
JM
, et al
.
Safety of diphtheria, tetanus, acellular pertussis and inactivated poliovirus (DTaP-IPV) vaccine
.
Vaccine
.
2014
;
32
(
25
):
3019
3024
21
Yih
WK
,
Nordin
JD
,
Kulldorff
M
, et al
.
An assessment of the safety of adolescent and adult tetanus-diphtheria-acellular pertussis (Tdap) vaccine, using active surveillance for adverse events in the Vaccine Safety Datalink
.
Vaccine
.
2009
;
27
(
32
):
4257
4262
22
Kulldorff
M
,
Davis
RL
,
Kolczak
M
, et al
.
A maximized sequential probability ratio test of drug and vaccine safety
.
Seq Anal
.
2011
;
30
:
58
78
23
Li
L
,
Kulldorff
M
.
A conditional maximized sequential probability ratio test for pharmacovigilance
.
Stat Med
.
2010
;
29
(
2
):
284
295
24
Lewis
E
,
Fireman
B
,
Klein
NP
,
Baxter
R
.
Exact sequential analysis for vaccine safety surveillance.
In:
Proceedings from the NFID 12th Annual Conference on Vaccine Research
;
April 27–29, 2009
25
McClure
DL
,
Xu
S
,
Weintraub
E
,
Glanz
JM
.
An efficient statistical algorithm for a temporal scan statistic applied to vaccine safety analyses
.
Vaccine
.
2012
;
30
(
27
):
3986
3991
26
Li
R
,
Stewart
B
,
Weintraub
E
.
Evaluating efficiency and statistical power of self-controlled case series and self-controlled risk interval designs in vaccine safety
.
J Biopharm Stat
.
2016
;
26
(
4
):
686
693
27
Joura
EA
,
Giuliano
AR
,
Iversen
OE
, et al;
Broad Spectrum HPV Vaccine Study
.
A 9-valent HPV vaccine against infection and intraepithelial neoplasia in women
.
N Engl J Med
.
2015
;
372
(
8
):
711
723
28
Crawford
NW
,
Clothier
HJ
,
Elia
S
, et al
.
Syncope and seizures following human papillomavirus vaccination: a retrospective case series
.
Med J Aust
.
2011
;
194
(
1
):
16
18
29
Klein
NP
,
Hansen
J
,
Chao
C
, et al
.
Safety of quadrivalent human papillomavirus vaccine administered routinely to females
.
Arch Pediatr Adolesc Med
.
2012
;
166
(
12
):
1140
1148
30
Vichnin
M
,
Bonanni
P
,
Klein
NP
, et al
.
An overview of quadrivalent human papillomavirus vaccine safety: 2006 to 2015
.
Pediatr Infect Dis J
.
2015
;
34
(
9
):
983
991
31
Bizjak
M
,
Bruck
O
,
Praprotnik
S
,
Dahan
S
,
Shoenfeld
Y
.
Pancreatitis after human papillomavirus vaccination: a matter of molecular mimicry
.
Immunol Res
.
2017
;
65
(
1
):
164
167
32
Das
A
,
Chang
D
,
Biankin
AV
,
Merrett
ND
.
Pancreatitis following human papillomavirus vaccination
.
Med J Aust
.
2008
;
189
(
3
):
178
33
Slade
BA
,
Leidel
L
,
Vellozzi
C
, et al
.
Postlicensure safety surveillance for quadrivalent human papillomavirus recombinant vaccine
.
JAMA
.
2009
;
302
(
7
):
750
757
34
US Food and Drug Administration
.
Gardasil 9
.
Silver Spring, MD
:
US FDA Original Application
;
2014
:
1
157
35
Moreira
ED
 Jr
,
Block
SL
,
Ferris
D
, et al
.
Safety profile of the 9-valent HPV vaccine: a combined analysis of 7 phase III clinical trials
.
Pediatrics
.
2016
;
138
(
2
):
e20154387
36
Centers for Disease Control and Prevention (CDC)
.
Syncope after vaccination–United States, January 2005-July 2007
.
MMWR Morb Mortal Wkly Rep
.
2008
;
57
(
17
):
457
460
37
Naleway
AL
,
Crane
B
,
Smith
N
, et al;
Vaccine Safety Datalink
.
Absence of venous thromboembolism risk following quadrivalent human papillomavirus vaccination, Vaccine Safety Datalink, 2008-2011
.
Vaccine
.
2016
;
34
(
1
):
167
171
38
Mullooly
JP
,
Donahue
JG
,
DeStefano
F
,
Baggs
J
,
Eriksen
E
;
VSD Data Quality Working Group
.
Predictive value of ICD-9-CM codes used in vaccine safety research
.
Methods Inf Med
.
2008
;
47
(
4
):
328
335

Competing Interests

POTENTIAL CONFLICT OF INTEREST: Dr Naleway has received funding from Pfizer for an unrelated study; Dr Klein reports research support from Merck (including a 9-valent human papillomavirus vaccine phase 4 postmarketing safety study), Sanofi Pasteur, GlaxoSmithKline, Protein Science (now Sanofi Pasteur), and Pfizer; the other authors have indicated they have no potential conflicts of interest to disclose.

FINANCIAL DISCLOSURE: Other than those listed in Potential Conflict of Interest, the authors have indicated they have no financial relationships relevant to this article to disclose.

Supplementary data