BACKGROUND

Temporal associations between Kawasaki disease (KD) and childhood vaccines have been reported. Limited data on KD following 13-valent pneumococcal conjugate (PCV13) and rotavirus vaccines are available.

METHODS

We conducted a self-controlled risk interval study using Vaccine Safety Datalink electronic health record data to investigate the risk of KD following PCV13 and rotavirus vaccines in children <2 years of age who were born from 2006 to 2017. All hospitalized KD cases identified by International Classification of Diseases diagnosis codes that fell within predefined risk (days 1–28 postvaccination) and control (days 29–56 for doses 1 and 2, and days 43–70 for doses 3 and 4) intervals were confirmed by manual chart review.

RESULTS

During the study period, 655 cases of KD were identified by International Classification of Diseases codes. Of these, 97 chart-confirmed cases were within risk or control intervals. In analyses, the age-adjusted relative risk for KD following any dose of PCV13 was 0.75 (95% confidence interval, 0.47–1.21). Similarly, the age-adjusted relative risk for KD following any dose of rotavirus vaccine was 0.66 (95% CI, 0.40–1.09). Overall, there was no evidence of an elevated risk of KD following PCV13 or rotavirus vaccines by dose. In addition, no statistically significant temporal clustering of KD cases was identified during days 1 to 70 postvaccination.

CONCLUSIONS

PCV13 and rotavirus vaccination were not associated with an increased risk of KD in children <2 years of age. Our findings provide additional evidence for the overall safety of PCV13 and rotavirus vaccines.

What’s Known on This Subject:

Temporal associations between Kawasaki disease and childhood vaccines have been reported. Although the evidence supports the safety of routine childhood vaccinations, there are limited data on Kawasaki disease following more recently licensed vaccines, including 13-valent pneumococcal conjugate and rotavirus vaccines.

What This Study Adds:

In self-controlled risk analyses, 13-valent pneumococcal conjugate and rotavirus vaccines were not associated with an increased risk of Kawasaki disease in children <2 years of age. Our findings provide additional evidence for the overall safety of these childhood vaccinations.

Kawasaki disease (KD) is a rare, potentially life-threatening disease of unknown etiology in children that is associated with vascular inflammation, most notably affecting the coronary arteries. Although the etiology of KD remains uncertain, a variety of theories have been proposed based upon pathologic, demographic, and epidemiologic data.1  One plausible theory is that KD may occur in genetically susceptible individuals because of an immunologic response triggered by any of a variety of infectious and/or environmental agents.1  The peak incidence of KD is at less than 2 years of age, at a time when many childhood vaccines are given. As a result, case reports noting the temporal associations between KD and childhood vaccines have been published.25  Because of the clinical significance of the disease and public concern, KD is considered by the Centers for Disease Control and Prevention (CDC) and Vaccine Safety Datalink (VSD) investigators as a high priority outcome to investigate as a potential vaccine-associated adverse event.6 

Prior studies of routine childhood immunizations (specifically, diphtheria-tetanus-acellular pertussis, 7-valent pneumococcal conjugate vaccine [PCV7], hepatitis B, Haemophilus influenzae type b, measles, mumps, and rubella, varicella, influenza, and poliovirus), including postmarketing surveillance using large databases, have not identified an association between KD and vaccination.710  In 2007, however, the US Food and Drug Administration revised the RotaTeq (a pentavalent rotavirus vaccine [RV5]) Prescribing Information to include information about KD cases identified in the large Phase 3 prelicensure trials and subsequently reported to the Vaccine Adverse Event Reporting System.11  Rotarix (a monovalent rotavirus vaccine [RV1]) Prescribing Information also includes information on an imbalance of KD cases between the RV1 recipients and the placebo recipients in prelicensure clinical trials.12  Furthermore, in a prior VSD study, the cumulative risk of KD within 28 days of PCV13 vaccination was compared with that of PCV7 vaccination. A nonstatistically significant increased risk of KD (relative risk [RR] 1.94, 95% CI: 0.79–4.86) was observed and the authors concluded that the possible association between PCV13 and KD deserved further investigation.13  In this study, we evaluated the association between KD and rotavirus and PCV13 vaccinations in the large VSD cohort using a self-controlled risk interval design (SCRI).

The VSD is a collaboration between the CDC and 9 integrated healthcare organizations that is used to conduct surveillance and epidemiologic research on vaccine safety.14  The VSD population encompasses approximately 12 million persons per year. VSD electronic data includes patient demographics, inpatient and outpatient encounters, and immunizations and can be supplemented with medical record review. Large cohorts of children can be followed while enrolled in a VSD site, and these data can be used to examine potential associations between vaccines and rare adverse events. For the current study, data were collected from 8 sites: Kaiser Permanente Northern California (CA), Kaiser Permanente Southern California (CA), Kaiser Permanente Washington (WA), Kaiser Permanente Colorado (CO), Denver Health (CO), Marshfield Clinic (WI), Kaiser Permanente Northwest (OR), and Health Partners Institute (MN). The study population consisted of children <24 months of age who were born between January 1, 2006 and December 31, 2017 (RV5, RV1, and PCV13 were introduced in the United States in 2006, 2008, and 2010, respectively). For inclusion in the study, children were required to have continuous enrollment from ≤42 days of life through either the index KD event or 23 months of age.

We conducted a SCRI analysis to compare the occurrence of KD between a risk interval 1 to 28 days postvaccination and an equal length 28-day comparison (control) interval following the risk interval (Fig 1).1517  SCRI analysis was chosen for the analytic method since it can inherently account for several important potential confounders including sex, VSD site, and season.18  The day of vaccination was day 0. The control intervals were days 29 through 56 for doses 1 and 2 (RV1, RV5, and PCV13) and days 43 to 70 for dose 3 (RV5 and PCV13) and dose 4 (PCV13) based on prior literature.13,19  The latter control interval (43–70 days postvaccination) was preferable to avoid misclassification of the outcome timing caused by uncertainty about the true period of the risk associated with vaccines.19  However, it is recommended by the CDC that PCV13 be given at 2, 4, 6, and 12 to 15 months of age, RV5 be given at 2, 4, and 6 months of age, and RV1 be given at 2 and 4 months of age.20  As such, the control interval for doses 1 and 2 was needed to fit into a 2-month period to avoid overlapped intervals with a subsequent dose. Vaccines were identified using CVX codes for PCV13, RV1, and RV5.21  We examined KD cases following any PCV13 vaccination, any rotavirus vaccinations, simultaneous PCV13 and rotavirus vaccinations, PCV13 vaccination without simultaneous rotavirus vaccination, and rotavirus vaccination without simultaneous PCV13 vaccination. Chart-confirmed KD cases were assigned to risk or control intervals based on the timing of their adjudicated symptom onset. Cases were excluded if the risk or control interval of 1 dose overlapped the risk or control interval of a subsequent dose.

FIGURE 1

Self-controlled risk interval design to compare the occurrence of Kawasaki disease between a risk interval and a control interval. Dose 1 and 2: RV1, RV5, and PCV13; dose 3: RV5 and PCV13; dose 4: PCV13.

FIGURE 1

Self-controlled risk interval design to compare the occurrence of Kawasaki disease between a risk interval and a control interval. Dose 1 and 2: RV1, RV5, and PCV13; dose 3: RV5 and PCV13; dose 4: PCV13.

Close modal

Cases of KD were identified by the International Classification of Diseases, Ninth Revision (ICD-9) code 446.1 (Kawasaki disease) and ICD-10 code M30.3 (acute febrile mucocutaneous lymph node syndrome) assigned in the inpatient setting. The diagnosis was limited to the inpatient setting for specificity as true KD cases are severe and require close monitoring and treatment with high-dose intravenous immunoglobulin (IVIG), which cannot be administered in outpatient settings. For children with multiple KD medical encounters, only the first inpatient visit was included. We conducted manual chart review of cases identified by ICD codes for KD that fell into the risk and control windows to confirm the diagnosis. Because KD symptom onset may precede the ICD code-based diagnosis date, it is possible that some cases with onset during the risk or control windows may be missed. To account for this, KD cases identified by ICD codes that occurred during a risk or control interval as well as cases identified during expanded intervals (ie, 70 days following dose 1 or 2 and 84 days following dose 3 or 4) were also selected for manual chart review (Supplemental Fig 4). We calculated the confirmation rates of this code-based algorithm for identification of incident KD.

The aim of medical record review was to confirm whether subjects had “complete KD” or “incomplete KD,” and to determine the onset date of symptoms. Verified complete KD cases met the CDC epidemiologic case definition for KD (Supplemental Table 3).8,22  Although the CDC epidemiologic case definition establishes the diagnosis of complete KD, these criteria do not identify all children with the illness and will leave out incomplete KD. As long term cardiac outcomes in children with incomplete KD are comparable to those of children with complete KD, we also used the American Heart Association KD scientific statement to identify incomplete KD.1,23  Verified incomplete KD was defined as a case has prolonged unexplained fever for ≥5 days (or fever until the date of administration of IVIG if it was given before the fifth day of fever), and 2 or 3 of the 5 principal clinical findings as described in Supplemental Table 3, and compatible laboratory or positive echocardiographic findings.1  Only chart-confirmed complete KD and incomplete KD were included in the final analyses.

Only cases of KD identified by ICD codes that were confirmed by manual chart review were analyzed. Descriptive analyses were conducted on the study cohort, including by age, sex, race and ethnicity, and season of year when the case occurred (winter = January–March, spring = April–June, summer = July–September, fall = October–December). The rates of chart-confirmed KD with onset in the risk interval were compared with the rates of chart-confirmed KD in the control interval using conditional Poisson regression. As age is a time-varying confounder, we adjusted for the differential risk of KD according to age in each interval using offset terms calculated based on data on incidence of KD by age within the VSD cohort. Additionally, analyses specific to dose number partially accounted for the confounding effect of age.24  As the risk and control intervals are narrow (28 days for each interval) and the incidence rates of KD in the United States were stable during the study period, calendar year of vaccination was not included in the models.25  Finally, a temporal scan analysis was conducted to detect any temporal clustering of chart-confirmed KD cases within 70 days after PCV13 or rotavirus vaccines.26,27  All statistical tests were 2-sided, and a P value of <.05 considered statistically significant. All data were analyzed using SAS version 9.4 (2019 SAS Institute, Cary, NC). The human subjects review board at each participating site and the CDC approved the study and determined informed consent was not required.

A total of 2 508 207 doses of PCV13 and 2 048 186 doses of rotavirus vaccines were administered to children <24 months of age who were born between January 1, 2006 and December 31, 2017. A total of 655 cases of KD were ascertained by ICD codes. Of these, 511 cases were excluded because they did not occur within the expanded intervals for manual chart review. Of 144 cases that underwent manual chart review, 122 cases were confirmed by manual medical record review, and the symptom onset of 97 cases were in the prespecified risk or control intervals (Fig 2). The chart confirmation rate was 85% overall and comprised 94 cases of complete KD (65%) and 28 incomplete KD (19%). An additional 13 cases (9%) did not meet the criteria for complete and incomplete KD and 9 (6%) were ruled out or lacked the necessary information to make a case determination.

FIGURE 2

Flowchart of Kawasaki disease cases confirmed by manual chart review, Vaccine Safety Datalink, 2006 to 2017.

FIGURE 2

Flowchart of Kawasaki disease cases confirmed by manual chart review, Vaccine Safety Datalink, 2006 to 2017.

Close modal

Demographic characteristics of chart-confirmed cases are shown in Table 1. The mean age of children with chart-confirmed KD was 7.1 months. In addition, 63% were male and 40% were Asian. Overall, there was no evidence of an elevated risk of KD following PCV13 or rotavirus vaccines by dose in SCRI analyses using chart-confirmed KD cases (Table 2). In SCRI analyses of chart-confirmed KD cases receiving any PCV13 vaccine with or without simultaneous rotavirus vaccines, there were 29 cases in the risk interval and 41 cases in the control interval (the age-adjusted RR, 0.75; 95% confidence interval [CI], 0.47–1.21). Similarly, in SCRI analyses using chart-confirmed KD cases who received any rotavirus vaccine with or without simultaneous PCV13, there were 24 cases in the risk interval and 41 cases in the control interval (the age-adjusted RR, 0.66; 95% CI, 0.40–1.09). In the temporal scan analysis using the symptom onset date of 97 chart-confirmed KD cases, no statistically significant clusters of KD cases were identified during days 1 to 70 postvaccination (P = .18) (Fig 3).

FIGURE 3

Temporal distribution of chart-confirmed Kawasaki Disease during days 1 to 70 postvaccinations (PCV13 and rotavirus vaccines), Vaccine Safety Datalink, 2006 to 2017.

FIGURE 3

Temporal distribution of chart-confirmed Kawasaki Disease during days 1 to 70 postvaccinations (PCV13 and rotavirus vaccines), Vaccine Safety Datalink, 2006 to 2017.

Close modal
TABLE 1

Characteristics of Patients <2 Years of Age Diagnosed with Kawasaki Disease, Vaccine Safety Datalink, 2006 to 2017

Demographic CharacteristicsChart-confirmed Cases, (n = 97)
Age (month), mean (SD) 7.1 (4) 
Sex, n (%)  
 Male 61 (63) 
 Female 36 (37) 
Race and ethnicity, n (%)  
 Non-Hispanic white 25 (26) 
 Non-Hispanic Black 3 (3) 
 Hispanic 22 (23) 
 American Indian or Alaska Native 0 (0) 
 Asian 39 (40) 
 Pacific Islander 2 (2) 
 Multiple races 1 (1) 
 Other 2 (2) 
 Unknown 3 (3) 
Season of year at diagnosis, n (%)  
 Spring 23 (24) 
 Summer 26 (27) 
 Autumn 25 (26) 
 Winter 23 (24) 
Demographic CharacteristicsChart-confirmed Cases, (n = 97)
Age (month), mean (SD) 7.1 (4) 
Sex, n (%)  
 Male 61 (63) 
 Female 36 (37) 
Race and ethnicity, n (%)  
 Non-Hispanic white 25 (26) 
 Non-Hispanic Black 3 (3) 
 Hispanic 22 (23) 
 American Indian or Alaska Native 0 (0) 
 Asian 39 (40) 
 Pacific Islander 2 (2) 
 Multiple races 1 (1) 
 Other 2 (2) 
 Unknown 3 (3) 
Season of year at diagnosis, n (%)  
 Spring 23 (24) 
 Summer 26 (27) 
 Autumn 25 (26) 
 Winter 23 (24) 
TABLE 2

Age-adjusted Relative Risks (RR) of Chart-confirmed Kawasaki Disease Following PCV13 and Rotavirus Vaccines by Dose, Vaccine Safety Datalink, 2006 to 2017

Cases in Risk IntervalCases in Control IntervalRR95%CIP
Any PCV13a      
 All 29 41 0.75 0.47–1.21 .24 
 Dose 1 10 0.40 0.11–1.47 .17 
 Dose 2 11 1.24 0.51–2.99 .63 
 Dose 3 10 10 1.02 0.42–2.45 .97 
 Dose 4 12 0.41 0.14–1.17 .09 
Any rotavirus vaccinesb      
 All 24 41 0.66 0.40–1.09 .11 
 Dose 1 17 0.40 0.15–1.09 .07 
 Dose 2 15 15 1.02 0.50–2.09 .96 
 Dose 3 0.45 0.14–1.48 .19 
Simultaneous PCV13 and rotavirus vaccines      
 All 16 22 0.81 0.43–1.55 .53 
 Dose 1 10 0.40 0.11–1.47 .17 
 Dose 2 11 1.24 0.51–2.99 .63 
 Dose 3 0.68 0.11–4.09 .68 
PCV13 without simultaneous rotavirus vaccines      
 All 13 19 0.68 0.34–1.39 .29 
 Dose 1 NA NA NA 
 Dose 2 NA NA NA 
 Dose 3 1.16 0.42–3.20 .77 
 Dose 4 12 0.41 0.14–1.17 .09 
Rotavirus vaccines without simultaneous PCV13      
 All 19 0.48 0.21–1.09 .08 
 Dose 1 0.40 0.08–1.95 .26 
 Dose 2 0.69 0.19–2.43 .56 
 Dose 3 0.34 0.07–1.69 .19 
Cases in Risk IntervalCases in Control IntervalRR95%CIP
Any PCV13a      
 All 29 41 0.75 0.47–1.21 .24 
 Dose 1 10 0.40 0.11–1.47 .17 
 Dose 2 11 1.24 0.51–2.99 .63 
 Dose 3 10 10 1.02 0.42–2.45 .97 
 Dose 4 12 0.41 0.14–1.17 .09 
Any rotavirus vaccinesb      
 All 24 41 0.66 0.40–1.09 .11 
 Dose 1 17 0.40 0.15–1.09 .07 
 Dose 2 15 15 1.02 0.50–2.09 .96 
 Dose 3 0.45 0.14–1.48 .19 
Simultaneous PCV13 and rotavirus vaccines      
 All 16 22 0.81 0.43–1.55 .53 
 Dose 1 10 0.40 0.11–1.47 .17 
 Dose 2 11 1.24 0.51–2.99 .63 
 Dose 3 0.68 0.11–4.09 .68 
PCV13 without simultaneous rotavirus vaccines      
 All 13 19 0.68 0.34–1.39 .29 
 Dose 1 NA NA NA 
 Dose 2 NA NA NA 
 Dose 3 1.16 0.42–3.20 .77 
 Dose 4 12 0.41 0.14–1.17 .09 
Rotavirus vaccines without simultaneous PCV13      
 All 19 0.48 0.21–1.09 .08 
 Dose 1 0.40 0.08–1.95 .26 
 Dose 2 0.69 0.19–2.43 .56 
 Dose 3 0.34 0.07–1.69 .19 

PCV13, 13-valent pneumococcal conjugate vaccine; NA, not applicable.

a

Any PCV13 with or without rotavirus vaccines.

b

Any Rotavirus vaccines with or without PCV13 vaccines.

In this study using US post-licensure vaccine safety monitoring surveillance data supplemented with standardized medical record review, we found no evidence of an increased KD risk following PCV13 or rotavirus vaccination in children under 2 years of age. As consistent with prior studies, males and Asian were more common in chart-confirmed KD cases in this study.1,24  Although the rate of incomplete KD (23%) is higher than prior reports, this is likely because our KD cases are limited to children less than 2 years of age (the median age, 7.1 months), and incomplete KD is more common in infants.1,28,29 

Although a prior study showed a nonstatistically significant increased point estimate risk of KD following PCV13, we did not identify a similar trend.13  Whereas the prior study used the VSD data between 2010 and 2012 (the number of vaccine doses of PCV13 included was 599 229), we expanded the study period through 2017, which increased the population size and resulted in over 2.5 million doses of PCV13 being included in the current study. As such, we were able to identify more KD cases following PCV13 vaccination than the prior study (70 cases vs 12 cases). Furthermore, we conducted a SCRI analysis with age-adjustment to account for important potential time-invariant and time-varying confounders, which reduces the risk of bias in our results. In fact, our results are consistent with the recent large studies on KD following PCV13 conducted in the Postlicensure Rapid Immunization Safety Monitoring program in the United States and in the Hospital Episodes Statistics system in the United Kingdom.19,30  Compared with these prior studies, we used manual chart review with stringent criteria for all KD cases identified by ICD codes.

In contrast to a previous ICD code-based study in Taiwan showing a small increased risk of KD during the third week post dose 2 of RV5 (incidence rate ratio 2.33, 95% CI 1.35–4.00), and fourth week post dose 1 of RV1 (incidence rate ratio 1.98, 95% CI 1.16–3.40), we did not detect any risks for KD among rotavirus vaccine recipients.31  Although this may be, in part, caused by the different risk intervals analyzed, it is more likely because of the standardized chart review process we used, which made misclassification of case status less likely. As KD symptom onset could precede physician diagnosis, the ICD code-based results of the study from Taiwan could have been influenced by misclassification of the symptom onset date in relation to the risk or control intervals. In addition, our findings on the safety of rotavirus vaccines are consistent with the results of other small or ICD code-based studies.3235 

Our study has several important strengths. First, we conducted manual chart reviews for all KD cases identified by ICD codes. This is particularly critical as diagnosing KD can be challenging and can only be made based on clinical criteria, which may pose a risk of misclassification of the outcome without chart-confirmation. Second, we used the standardized case definitions based on the CDC and American Heart Association criteria, which included supplemental laboratory and echocardiographic findings. Therefore, our criteria for incomplete KD were more stringent than prior studies and provide more accurate classification of the outcome. Third, we used a SCRI analysis with age-adjustment using a large VSD cohort with comprehensive vaccination and follow up information.36 

This study is subject to limitations. First, we identified cases of KD only in the inpatient setting and we did not include KD cases which did not require hospitalization. However, omission of KD cases would have been rare given that patients with a KD diagnosis generally require inpatient care for high-dose IVIG treatment, which is infused for a much longer time (eg, 10 hours or more) than for other indications and requires close monitoring during and after the infusion. Second, another potential limitation of this study is bias from misclassification of the exposure. It is possible that some children in the study received vaccines outside their VSD health care system, which were not captured in the study databases. We expect this to be rare as we required continuous enrollment in the system up until the time of the diagnosis of KD. Third, although our sample size was relatively large, this study was not powered to conduct subanalyses for risk assessment by outcome type (complete versus incomplete KD) or vaccine brand type (RV1 versus RV5). In addition, the interpretation for the secondary findings by vaccine dose may be limited because of the relatively small number of cases for each dose. Finally, we did not adjust for simultaneous vaccines, such as diphtheria-tetanus-acellular pertussis. However, the majority of vaccines in this age group are given together with other vaccines in accordance with the Advisory Committee on Immunization Practices schedule, and these other vaccines have been previously shown to confer no increased risk of KD.8 

Vaccines including PCV13 and rotavirus vaccines have prevented substantial morbidity and mortality in children and also have other indirect effects, including herd protection, decreasing the use of broad-spectrum antimicrobial agents, and decreasing antimicrobial drug resistance by reducing the circulation of multidrug-resistant bacterial pathogens.37,38  However, public concerns about the safety of the immunization schedule threaten to undercut these substantial benefits. In particular, there is parental concern that the increasing prevalence of allergic diseases or autoimmune diseases may be associated with the increasing number of vaccines administered during early childhood.39  As KD is thought to be an immunologic phenomenon, our findings of the safety of these vaccines at a young age provide additional reassurance to concerned families, healthcare providers, and the public.

Our extensive chart-confirmed SCRI study showed no evidence of an increased risk of KD following PCV13 or rotavirus vaccinations. These findings provide additional safety evidence of PCV13 and rotavirus vaccines to support the current US childhood immunization schedule.

We thank Vaccine Safety Datalink project managers for their contributions to overall project management, data managers for their contributions to overall data management, and medical record abstractors for their contributions to medical record review; and Joseph Y. Abrams, PhD (Centers for Disease Control and Prevention) for providing critical feedback for this study.

Dr Kamidani conceptualized and designed the study, designed the data collection instruments, participated in data collection, participated in the analysis, drafted the initial manuscript, and reviewed and revised the manuscript; Dr Panagiotakopoulos participated in the conceptualization and design of the study, designed the data collection instruments, and reviewed and revised the manuscript; Mr Weintraub and Dr Licata participated in the conceptualization and design of the study, conducted the analysis, and critically reviewed the manuscript; Dr Belay contributed to the interpretation of data and reviewed and revised the manuscript; Drs Daley, Yih, Zerbo, Tseng, DeSilva, Nelson, Williams, Hambidge, and Donahue and Ms Groom participated in the design of the study, contributed to the acquisition and interpretation of data, and reviewed and revised the manuscript; all authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

DISCLAIMER: The manuscript was approved by the Centers for Disease Control and Prevention clearance process. The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention. Use of trade names is for identification only and does not imply endorsement by the Public Health Service or the US Department of Health and Human Services.

FUNDING: This research is supported by the Centers for Disease Control and Prevention, and Emory University Warshaw Fellow Research Grant. The Vaccine Safety Datalink Project is funded by the Centers for Disease Control and Prevention. The funder of this study had no role in the design of the study, data collection, data analysis, data interpretation, or writing of the report.

CONFLICT OF INTERST DISCLOSURES: Dr Kamidani’s institution, Emory University, has received funding from the National Institutes of Health to conduct clinical trials of Moderna and Janssen COVID-19 vaccines and funding from Pfizer to conduct clinical trials of Pfizer-BioNTech coronavirus disease 2019 vaccines. Dr Yih has received research support from GlaxoSmithKline and Pfizer for unrelated studies. Dr Tseng has received research support from Moderna, GlaxoSmithKline, and Janssen for unrelated studies. Dr DeSilva has received research support from National Institutes of Health for unrelated studies. Dr Nelson reported receiving payments from Moderna, Southern California Permanente Medical Group, Elsevier, and Harvard Pilgrim Health Care outside the submitted work. The other authors have indicated they have no potential conflicts of interest to disclose.

CDC

Centers for Disease Control and Prevention

CI

confidence interval

ICD

International Classification of Diseases

IVIG

intravenous immune globulin

KD

Kawasaki disease

PCV

pneumococcal conjugate vaccine

RR

relative risk

SCRI

self-controlled risk interval

VSD

Vaccine Safety Datalink

1
McCrindle
BW
,
Rowley
AH
,
Newburger
JW
, et al;
American Heart Association Rheumatic Fever, Endocarditis, and Kawasaki Disease Committee of the Council on Cardiovascular Disease in the Young
;
Council on Cardiovascular and Stroke Nursing
;
Council on Cardiovascular Surgery and Anesthesia
;
and Council on Epidemiology and Prevention
.
Diagnosis, treatment, and long-term management of Kawasaki disease: a scientific statement for health professionals from the American Heart Association
.
Circulation
.
2017
;
135
(
17
):
e927
e999
2
Yin
S
,
Liubao
P
,
Chongqing
T
,
Xiaomin
W
.
The first case of Kawasaki disease in a 20-month old baby following immunization with rotavirus vaccine and hepatitis A vaccine in China: a case report
.
Hum Vaccin Immunother
.
2015
;
11
(
11
):
2740
2743
3
Shimada
S
,
Watanabe
T
,
Sato
S
.
A patient with Kawasaki disease following influenza vaccinations
.
Pediatr Infect Dis J
.
2015
;
34
(
8
):
913
4
Kraszewska-Glomba
B
,
Kuchar
E
,
Szenborn
L
.
Three episodes of Kawasaki disease including one after the Pneumo 23 vaccine in a child with a family history of Kawasaki disease
.
J Formos Med Assoc
.
2016
;
115
(
10
):
885
886
5
Miron
D
,
Fink
D
,
Hashkes
PJ
.
Kawasaki disease in an infant following immunisation with hepatitis B vaccine
.
Clin Rheumatol
.
2003
;
22
(
6
):
461
463
6
Glanz
JM
,
Newcomer
SR
,
Jackson
ML
, et al
.
White paper on studying the safety of the childhood immunization schedule in the Vaccine Safety Datalink
.
Vaccine
.
2016
;
34
(
Suppl 1
):
A1
A29
7
Hua
W
,
Izurieta
HS
,
Slade
B
, et al
.
Kawasaki disease after vaccination: reports to the vaccine adverse event reporting system 1990-2007
.
Pediatr Infect Dis J
.
2009
;
28
(
11
):
943
947
8
Abrams
JY
,
Weintraub
ES
,
Baggs
JM
, et al
.
Childhood vaccines and Kawasaki disease, Vaccine Safety Datalink, 1996–2006
.
Vaccine
.
2015
;
33
(
2
):
382
387
9
Phuong
LK
,
Bonetto
C
,
Buttery
J
, et al;
Brighton Collaboration Kawasaki Disease (KD) Working Group
.
Kawasaki disease and immunisation: a systematic review
.
Vaccine
.
2017
;
35
(
14
):
1770
1779
10
Hall
GC
,
Tulloh
RM
,
Tulloh
LE
.
The incidence of Kawasaki disease after vaccination within the UK pre-school National Immunisation Programme: an observational THIN database study
.
Pharmacoepidemiol Drug Saf
.
2016
;
25
(
11
):
1331
1336
11
US Food and Drug Administration
.
RotaTeq
.
12
US Food and Drug Administration
.
Rotarix
.
13
Tseng
HF
,
Sy
LS
,
Liu
IL
, et al
.
Postlicensure surveillance for pre-specified adverse events following the 13-valent pneumococcal conjugate vaccine in children
.
Vaccine
.
2013
;
31
(
22
):
2578
2583
14
Centers for Disease Control and Prevention
.
Vaccine safety datalink
.
15
Baker
MA
,
Lieu
TA
,
Li
L
, et al
.
A vaccine study design selection framework for the postlicensure rapid immunization safety monitoring program
.
Am J Epidemiol
.
2015
;
181
(
8
):
608
618
16
Greene
SK
,
Kulldorff
M
,
Lewis
EM
, et al
.
Near real-time surveillance for influenza vaccine safety: proof-of-concept in the Vaccine Safety Datalink Project
.
Am J Epidemiol
.
2010
;
171
(
2
):
177
188
17
Duffy
J
,
Lewis
M
,
Harrington
T
, et al
.
Live attenuated influenza vaccine use and safety in children and adults with asthma
.
Ann Allergy Asthma Immunol
.
2017
;
118
(
4
):
439
444
18
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
19
Baker
MA
,
Baer
B
,
Kulldorff
M
, et al
.
Kawasaki disease and 13-valent pneumococcal conjugate vaccination among young children: a self-controlled risk interval and cohort study with null results
.
PLoS Med
.
2019
;
16
(
7
):
e1002844
20
Centers for Disease Control and Prevention
.
Child and adolescent immunization schedule
. .
21
HL7 FHIR
.
HL7 FHIR
.
22
Centers for Disease Control and Prevention
.
Kawasaki syndrome
.
https://www.cdc.gov/kawasaki/. Accessed October, 30, 2019
23
Cho
MA
,
Choi
YJ
,
Jung
JW
.
Affects of “age at diagnosis” on coronary artery lesions in patients with incomplete kawasaki disease
.
Korean Circ J
.
2010
;
40
(
6
):
283
287
24
Holman
RC
,
Belay
ED
,
Christensen
KY
,
Folkema
AM
,
Steiner
CA
,
Schonberger
LB
.
Hospitalizations for Kawasaki syndrome among children in the United States, 1997-2007
.
Pediatr Infect Dis J
.
2010
;
29
(
6
):
483
488
25
Vasudeva
R
,
Poku
FA
,
Thommana
M
, et al
.
Trends and resource utilization in Kawasaki disease hospitalizations in the United States, 2008-2017
.
Hosp Pediatr
.
2022
;
12
(
3
):
257
266
26
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
27
Kulldorff
M
.
A spatial scan statistic
.
Commun Stat Theory Methods
.
1997
;
26
(
6
):
1481
1496
28
Fukushige
J
,
Takahashi
N
,
Ueda
Y
,
Ueda
K
.
Incidence and clinical features of incomplete Kawasaki disease
.
Acta Paediatr
.
1994
;
83
(
10
):
1057
1060
29
Yeom
JS
,
Woo
HO
,
Park
JS
,
Park
ES
,
Seo
JH
,
Youn
HS
.
Kawasaki disease in infants
.
Korean J Pediatr
.
2013
;
56
(
9
):
377
382
30
Stowe
J
,
Andrews
NJ
,
Turner
PJ
,
Miller
E
.
The risk of Kawasaki disease after pneumococcal conjugate & meningococcal B vaccine in England: a self-controlled case-series analysis
.
Vaccine
.
2020
;
38
(
32
):
4935
4939
31
Huang
WT
,
Juan
YC
,
Liu
CH
,
Yang
YY
,
Chan
KA
.
Intussusception and Kawasaki disease after rotavirus vaccination in Taiwanese infants
.
Vaccine
.
2020
;
38
(
40
):
6299
6303
32
Layton
JB
,
Butler
AM
,
Panozzo
CA
,
Brookhart
MA
.
Rotavirus vaccination and short-term risk of adverse events in US infants
.
Paediatr Perinat Epidemiol
.
2018
;
32
(
5
):
448
457
33
Loughlin
J
,
Mast
TC
,
Doherty
MC
,
Wang
FT
,
Wong
J
,
Seeger
JD
.
Postmarketing evaluation of the short-term safety of the pentavalent rotavirus vaccine
.
Pediatr Infect Dis J
.
2012
;
31
(
3
):
292
296
34
Belongia
EA
,
Irving
SA
,
Shui
IM
, et al;
Vaccine Safety Datalink Investigation Group
.
Real-time surveillance to assess risk of intussusception and other adverse events after pentavalent, bovine-derived rotavirus vaccine
.
Pediatr Infect Dis J
.
2010
;
29
(
1
):
1
5
35
Hoffman
V
,
Abu-Elyazeed
R
,
Enger
C
, et al
.
Safety study of live, oral human rotavirus vaccine: a cohort study in United States health insurance plans
.
Hum Vaccin Immunother
.
2018
;
14
(
7
):
1782
1790
36
McNeil
MM
,
Gee
J
,
Weintraub
ES
, et al
.
The Vaccine Safety Datalink: successes and challenges monitoring vaccine safety
.
Vaccine
.
2014
;
32
(
42
):
5390
5398
37
Moore
MR
,
Link-Gelles
R
,
Schaffner
W
, et al
.
Effect of use of 13-valent pneumococcal conjugate vaccine in children on invasive pneumococcal disease in children and adults in the USA: analysis of multisite, population-based surveillance
.
Lancet Infect Dis
.
2015
;
15
(
3
):
301
309
38
National Vaccine Advisory Committee
.
A call for greater consideration for the role of vaccines in national strategies to combat antibiotic-resistant bacteria: recommendations from the National Vaccine Advisory Committee: approved by the National Vaccine Advisory Committee on June 10, 2015
.
Public Health Rep
.
2016
;
131
(
1
):
11
16
39
Hilton
S
,
Petticrew
M
,
Hunt
K
.
‘Combined vaccines are like a sudden onslaught to the body’s immune system’: parental concerns about vaccine ‘overload’ and ‘immune-vulnerability’
.
Vaccine
.
2006
;
24
(
20
):
4321
4327

Supplementary data