Introduction

Immune thrombocytopenic purpura (ITP) is a haemorrhagic disorder characterised by thrombocytopenia, a purpuric rash, normal bone marrow and the absence of signs of other identifiable causes of thrombocytopenia [1,2,3,4,5]. Previous studies have shown that ITP risk increases after measles, rubella, chickenpox and influenza infection [2, 6,7,8]. Recently, live measles–mumps–rubella (MMR) [8,9,10,11] and varicella [12] vaccines have been reported to potentially increase the risk for ITP. Studies have also suggested that inactivated hepatitis B [13] and diphtheria–tetanus–acellular pertussis (DTaP) [12] vaccines may increase ITP risk.

Evidence suggests that ITP risk after vaccination increases through the same mechanism as that by which microbial infections induce antiplatelet autoantibodies [14]. Because vaccines are designed to induce protective immunity by mimicking infections in the human body, both live and inactivated vaccinations can theoretically trigger the development of ITP [15]. Researchers have investigated associations between multiple vaccinations and ITP risk, but findings regarding inactivated vaccinations were inconsistent [16,17,18,19,20]. However, because many vaccines are administered to children in infancy, an epidemiological study for vaccine safety conducted some time after a single vaccine dose is difficult to design. Therefore, more studies are needed to clarify the ITP risk associated with live and inactivated vaccines.

A combined vaccine is defined as a single product in which equivalent component vaccines are administered as a single vaccine to prevent more than one disease or to protect against multiple strains of infectious microbes [21]. For example, MMR, diphtheria–pertussis–tetanus and multivalent pneumococcal conjugate vaccines have been licensed for large-scale supply. The nature of combined vaccines means it is difficult to separately measure the adverse reaction rate of each component.

Simultaneous vaccination, which is commonly conducted for children in Europe [22], North America [23, 24] and Asia [25], is defined as administering more than one vaccine at different anatomic sites during the same clinic visit, without combining these vaccines in the same syringe [21]. Early research reported that parents were concerned about the efficacy of simultaneous immunisation [26]. However, experimental and clinical epidemiological studies provided a scientific basis for the efficacy of this practice [27]. Recently, the focus of interest in simultaneous vaccination has moved from questions of efficacy to questions of safety [28,29,30]. If a vaccination has adverse reaction risks, these risks may accumulate in simultaneous vaccination [31]. However, no research has compared the risks when administering vaccines simultaneously versus separately in the same population.

Studies conducted in the 1990s indicated rates of adverse reactions were similar between simultaneous and separate vaccine administration [32,33,34]. Correspondingly, the Centers for Disease Control and Prevention issued a guideline recommending that parents use simultaneous vaccination for children younger than 2 years to avoid missing the appropriate vaccination timing [21]. However, some studies suggested that rates of several adverse reactions may increase with simultaneous vaccination [27, 35]. Therefore, the risk associated with simultaneous vaccine administration needs to be investigated more comprehensively.

In this case–control study, we aimed to determine the ITP risk after live, inactivated and simultaneous vaccination in Japan. We also examined the risk associated with simultaneous vaccination among subjects aged < 2 years, referred to in this study as infants [36], because this age group are frequently immunised using this method.

Materials and methods

Enrolment

We requested physicians from paediatrics and internal medicine departments at seven university hospitals and 18 regional centre hospitals in Japan to recruit inpatients and outpatients for this study. Recruitment occurred from 1 October 2015 to 27 March 2017. During this period, participating physicians enrolled in all new cases with ITP that attended their hospital and met the inclusion criteria for this study. These physicians filled in the study questionnaire with the requested data for the case and control patients.

Vaccines and other measurements

We defined exposure as vaccination within 28 days before the onset of ITP. To measure this exposure, participating physicians who treated case (ITP) and control (other diseases) patients completed questionnaires covering retrospective information on vaccination history and other characteristics. In Japan, vaccination history for infants and preschool children is recorded in the Maternal and Child Health Handbook, which is provided by the municipality. Vaccination history for other children and adults was captured in a medical interview conducted by physicians or obtained from medical records. Additional information was obtained from hospital medical records. In Japan, the Immunization Law defines ITP as an adverse event following immunisation that occurs within 28 days after vaccine administration [37]. Based on this law, we considered ≤ 28 days after administration as the duration in which vaccinated people could potentially develop ITP as a result of vaccination.

The live vaccines investigated were the bacillus Calmette–Guérin (BCG), rotavirus, varicella, mumps and measles–rubella (MR) vaccines (in Japan, MR and mumps vaccines are administered separately). Inactivated vaccines were the influenza, diphtheria–pertussis–tetanus–polio (DPT-IPV), hepatitis B, 13-valent pneumococcal conjugate, Haemophilus influenzae type b (Hib) and Japanese encephalitis vaccines. Manufacturers of vaccines administered in Japan are Daiichi Sankyo Company Ltd.; Takeda Pharmaceutical Company Ltd.; KM Biologics Company Ltd.; Biken Company Ltd.; Denka Seiken Company Ltd.; Japan BCG Laboratory; Merck Sharp & Dohme Corp.; Sanofi S.A.; Pfizer Inc.; and GlaxoSmithKline K.K.

Cases and controls

Following Japanese and American guidelines [38,39,40,41], cases with ITP were identified when patients met all of the following conditions. (1) Peripheral blood platelet count ≤ 100,000/µL [41,42,43]. (2) Without anaemia unless the patient was bleeding as a result of ITP. (3) Without deformation of red or white blood cells. (4) Without aplastic anaemia, myelodysplastic syndrome, leukaemia, malignant lymphoma, paroxysmal nocturnal haemoglobinuria, systemic lupus erythematosus, disseminated intravascular coagulation, bone marrow metastasis, myelofibrosis, thrombotic thrombocytopenic purpura, haemolytic uraemic syndrome, liver cirrhosis, hypersplenism and megaloblastic anaemia. (5) Without Bernard–Soulier, Wiskott–Aldrich, May–Hegglin and Kasabach–Merrit syndromes. (6) Condition was not pseudo-thrombocytopenia due to ethylenediaminetetraacetic acid. (7) Condition was not thrombocytopenia caused by pharmacological agents, radiation or measles. Because this study included patients with new onset ITP, all ITP cases were considered to be as acute [42].

Each patient with ITP (case) was matched with one or two control patients [29, 30]. Participating physicians matched controls with case patients by the institution, timing of hospital visit (within a 1-month difference), sex and age. In the matched study design, we permitted overlapping use of case patient datum to be paired with two different control patient data; within an institution, one control datum was paired with two different case patient data. We requested that participating physicians reduced the difference in timing of hospital visit between case and control patients to within 2 weeks where possible. Age was matched as age in months for case patients aged under 1 year, and within 2 years for case patients aged 1–17 years. For case patients aged ≥ 18 years, age was matched within 2 years where possible, with a maximum difference of 5 years. All participating patients were of Asian ethnicity. Sex, age at hospital visit and consulted department were measured; height and body weight were not measured.

Statistical analysis

To exploit the matched study design, we calculated McNemar’s odds ratios (ORs) and their 95% confidence intervals (CIs) [44] for the primary outcome. McNemar’s OR is suitable for matched case–control design. Briefly, in this calculation method, vaccination histories for case–control pairs were grouped into four profiles: (1) vaccination ( +) for a patient with immune thrombocytopenic purpura (case) and vaccination ( +) for a patient with other diseases (control); (2) vaccination ( +) for case and vaccination ( −) for control; (3) vaccination ( −) for case and vaccination ( +) for control; and (4) vaccination ( −) for case and vaccination ( −) for control. In calculating McNemar’s ORs, only two of these profiles (2 and 3) were used based on an assumption of the binary distribution. For example, if the number of the case–control vaccination profile 2 equalled that of profile 3, McNemar’s OR = 1 (null hypothesis) [45]. If the former number of profile 2 doubled the latter number of profile 3, McNemar’s OR = 2. This calculation method originated when McNemar’s OR was developed to evaluate the results of a matched case–control study. Subsequently, conditional logistic regression was developed to further adjust for covariates.

We also calculated McNemar’s ORs stratified by age group. Because vaccination for measles, rubella and varicella are known risk factors for ITP onset [8,9,10,11,12], we calculated those ORs as the reference values. Tests for p values were based on a binary distribution with the null hypothesis being the same pair of numbers between the two vaccination history profiles, which is usually used in calculating NcNemar’s OR [44].

We also calculated common ORs to exploit all available data for unmatched case and control patients. Data for those patients whose paired partner was excluded from the matched data analysis were included in the unmatched analysis. Common ORs were the secondary outcome. We aimed to measure confounders of history of infection with helicobacter pylori [42, 46], other viruses and bacteria, but only obtained these data for a subset of patients. Therefore, in estimating common ORs, we only adjusted for sex [10, 47, 48], and calculated ORs stratified by age group. The age groups were < 2, 2–5, 6–17 and ≥ 18 years, based on the vaccination schedule for children, adults and older adults recommended by the Japan Pediatric Society and the Infectious Disease Surveillance Center (English versions) [25, 49].

Because the risk associated with simultaneous vaccination could not be assessed in the matched analysis because of the small sample size, this risk was only assessed in the unmatched analysis. We calculated sex-adjusted common ORs for ITP occurrence in the simultaneous administration of two or more vaccines. We performed all statistical analyses using SAS version 9.4 (SAS Institute, Cary, NC, USA). All reported p values were two-sided, and p < 0.05 was considered significant.

Results

Participants

In total, 114 pairs were included in the matched analyses, and 175 patients in the unmatched analyses. Most patients with ITP were enrolled from paediatrics departments. In the ITP case and control groups, varied infectious diseases were reported 1 month before the hospital visit. There were no underlying diseases involving an immunocompromised status, but there was an epilepsy case in the ITP case group and an asthma case in the control group. Table 1 shows the characteristics of the matched case and control patients. In total, 64 participants (36.5%) were under 2 years of age, 52 (29.7%) were aged 2–5 years and 76 (43.4%) were male. For the 16 cases on whom data was provided by their physicians, the mean (standard deviation) of duration from vaccination to ITP onset was 17.4 (11.1) days. Among case patients, there were 18 cases of purpuras, three of bleeding in the oral cavity and difficulty stopping bleeding, three of hypermenorrhoea and one of haematochezia. Diagnoses recorded for hospital visits among control patients included 16 respiratory tract infections, nine digestive tract infections, three urinary tract infections, nine allergic diseases and four malignant neoplasms. There were 16 positive results for any helicobacter pylori test in the case group, but no positive results in the control group. Figure 1 presents a histogram of the ages of case patients. Most case patients with ITP were aged < 1–14 years.

Table 1 Characteristics of matched immune thrombocytopenic purpura case and control patients in this case–control study
Fig. 1
figure 1

Histogram of ages among case patients with immune thrombocytopenic purpura

Pairwise analyses

Table 2 shows the McNemar’s ORs, stratified by live/inactivated vaccine and age group. For all age groups, the McNemar’s OR for ITP occurrence was 1.80 (95% CI 0.54–6.84, p = 0.64) for all vaccines; 1.00 (95% CI 0.19–5.37, p = 0.50) for live vaccines; 2.33 (95% CI 0.53–14.0, p = 0.70) for inactivated vaccines; and 2.00 (95% CI 0.03–157, p = 1.00) for the MR vaccine. The OR for the varicella vaccine could not be calculated because of an unbalanced number combination in the four vaccination history profiles of case–control pairs. In infants aged under 2 years, the McNemar’s ORs for ITP occurrence in cases for all vaccines, live vaccines and inactivated vaccines were 1.50 (95% CI 0.17–18.0, p = 0.50), 2.00 (95% CI 0.29–22.1, p = 0.67) and 1.00 (95% CI 0.01–78.5, p = 0.50), respectively. Supplementary Tables 1–3 present cross-tabulations of the four profiles of matched pairs for all, live and inactivated vaccines.

Table 2 McNemar’s odds ratios (95% confidence intervals) for immune thrombocytopenic purpura occurrence after vaccination by age group

Non-pairwise analyses

Table 3 shows the sex-adjusted common ORs stratified by live/inactivated vaccine and age group. The common OR for ITP occurrence was 1.39 (95% CI 0.61–3.14, p = 0.50) for all vaccines, 2.18 (95% CI 0.53–9.40, p = 0.34) for live vaccines, 1.28 (95% CI 0.52–3.08, p = 0.69) for inactivated vaccines and 3.55 (95% CI 0.18–216, p = 0.60) for the MR vaccine. The OR for the varicella vaccine could not be calculated because of the small number of ITP cases in this exposure. In infants, the common ORs for ITP occurrence in cases for all vaccines, live vaccines and inactivated vaccines were 1.43 (95% CI 0.43–4.68, p = 0.68), 4.26 (95% CI 0.80–29.3, p = 0.10) and 0.99 (95% CI 0.26–3.55, p = 1.00), respectively. The common OR for inactivated vaccines was significantly high among children aged 6–17 years (OR 13.6, 95% CI 1.23–747, p = 0.028). Supplementary Tables 4–21 present the cross-tabulations of individual vaccines for the unmatched data. The numbers of case and control patients differed between pairwise and non-pairwise analyses.

Table 3 Sex-adjusted odds ratios (95% confidence intervals) for immune thrombocytopenic purpura occurrence after vaccination by age group

Simultaneous vaccination

For simultaneous vaccination, we investigated children aged 2 months to 11 years: six children had two vaccines, four had three vaccines, two had four vaccines and one had five vaccines. Table 4 shows the sex-adjusted common ORs for ITP occurrence in the simultaneous administration of two or more vaccines. The ORs were 1.52 (95% CI 0.45–5.21, p = 0.71) in children aged < 12 years and 1.95 (95% CI 0.44–8.72, p = 0.28) in children aged < 2 years. The ORs for the simultaneous administration of inactivated vaccines were 1.83 (95% CI 0.44–7.59, p = 0.40) in children of all ages and 1.41 (95% CI 0.29–6.94, p = 0.67) in children aged < 2 years. The ORs for the simultaneous administration of live and inactivated vaccines were 1.36 (95% CI 0.29–6.30, p = 0.69) in children of all ages and 2.85 (95% CI 0.43–18.9, p = 0.28) in children aged < 2 years.

Table 4 Sex-adjusted odds ratios (95% confidence intervals) for immune thrombocytopenic purpura occurrence in simultaneous administration of vaccines among children aged < 12 years

Discussion

Statement of principal findings

The McNemar’s and common ORs for ITP occurrence for all, live and inactivated vaccines ranged from almost null to more than 1 among people of all ages, children aged < 18 years and children aged < 2 years, but these results were not statistically significant. The common ORs for simultaneous administration of two or more vaccines were greater than 1 for all, live and inactivated vaccines; these results were not statistically significant.

Results in the context of previous reports

Previous research suggests that people may develop ITP as an immune reaction to vaccination in combination with a genetic predisposition [4, 50,51,52,53,54]. Because vaccines are designed to mimic real infections and trigger immune reactions, clinicians have been concerned that vaccination—particularly the administration of live vaccines—may cause autoimmune disease [50]. MMR-associated acute ITP has been frequently investigated in studies on this topic [9, 51,52,53]. In 1994, the United States Institute of Medicine acknowledged that evidence had established a causal relationship between MMR vaccination and thrombocytopenia [55]. In our study, the common OR for MR vaccination was 3.55, although this result was not statistically significant (Supplementary Table 15). In summary, the present study did not find many significant effects, but the McNemar’s and common ORs estimated for live vaccination ranged from null to high across age groups.

ITP risk associated with inactivated vaccination has not been sufficiently investigated in previous studies. A post-licensure retrospective study in the United States did not find any ITP cases among 2-month-old infants vaccinated with the DTaP-IPV/Hib or other DTaP-containing vaccine from 1 October 2008 to 31 July 2010 [17]. A Canadian surveillance study (1992–2001) identified eight ITP cases after DTP or DTaP ± IPV ± Hib vaccinations, three after hepatitis B vaccination and two after influenza vaccination in 12 paediatric hospitals [18]. An Italian post-licensure prospective study involving a large number of older adults conducted from 2006 to 2009 reported 13 cases of ITP after seasonal influenza vaccinations [19]. A study using a Taiwanese national database found an incidence rate ratio for ITP of 1.09 (95% CI 0.65–1.85) following pandemic influenza (H1N1) 2009 monovalent vaccination without adjuvant versus non-vaccination among people aged 0–42 years, and of 1.05 (95% CI 0.41–2.85) among children [20]. The United States Vaccine Adverse Event Reporting System database reported that from 1990 to 2008, influenza vaccination had a proportional reporting ratio (an index of identifying unexpectedly frequent reports) for thrombocytopenia occurrence of 2.10 (17 cases, lower bound of 90% CI of the empirical Bayesian geometric mean < 2.0) [16]. However, there is a lack of consistent evidence about the ITP risk associated with inactivated vaccines.

In our study, the estimated McNemar’s ORs for inactivated vaccines varied from under 1 to over 1, and were not statistically significant in the examined age groups. Calculation of McNemar’s ORs used discordant case–control pairs in terms of vaccination history; therefore, calculating McNemar’s ORs required a much larger sample size than would have been required to calculate common ORs [56]. Because of insufficient power, we were unable to clarify whether ITP risk accompanied inactivated vaccination.

Contributions to clinicians and parents

Care should be taken in the interpretation of the ORs reported for inactivated vaccines in this study. In the context of Japan, many inactivated vaccinations are scheduled for children aged < 1 year [25], meaning that infants of this age have a high chance of being immunised with inactivated vaccines. Previous research suggests that among Japanese children, the peak age of ITP onset is < 2 years [48]. The ORs for inactivated vaccination could, therefore, be overestimated in this study because most data used for their estimation included data for infants aged < 1 year, who frequently receive inactivated vaccines and have a relatively high incidence of ITP. Therefore, we interpreted the McNemar’s OR of 1.00 (Table 2) and common OR of 0.99 (Table 3) among children aged < 2 years as indicating that this study did not detect ITP risk in inactivated vaccination.

Vaccination is frequently criticised in the media [29, 30]. However, vaccination has three types of benefits: direct effects (immunity of the vaccinated person), indirect effects (herd protection of non-vaccinated people) and reduction in the risk for developing ITP caused by infection with viruses or bacteria. The third type of benefit could be rephrased as the observation that the incidence of ITP following vaccination is lower than the incidence of ITP after infection with wild viruses [50]. The incidence of ITP development in children immunised with MMR is reportedly around one in 30,000, whereas it is approximately one in 3000 following natural rubella infection and approximately one in 6000 following measles infection. Another report suggested that the incidence of ITP caused by MMR vaccination is 10–20 times lower than the incidence of ITP following natural infection [6]. As is inherent with a single outcome (in this case, the development of ITP), the present findings of ORs greater than 1 do not consider the third type of vaccination benefit. In other words, it is difficult to investigate the long-term effect of vaccination on preventing ITP development from the perspective of the third type of benefit—if it exists—in cohort or case–control studies.

The present study is notable for investigating ITP risk following simultaneous vaccination. Although we did not find statistically significant results, the sex-adjusted ORs were relatively high (OR 1.52 among children aged < 12 years and 1.95 among children aged < 2 years; Table 3). The point estimate of the OR was larger than that of single vaccination (1.38 among children aged < 12 years and 1.21 among children aged < 2 years; Table 3) in children of each age group. Few studies have investigated the risks associated with simultaneous vaccination. A randomised controlled trial involving people receiving simultaneous influenza vaccination reported increased haemagglutinin inhibition titre, seroconversion and seropositivity rates, suggesting that the immune response was boosted [57]. Because vaccine efficacy and immune-related adverse reactions have the same origin [15], the rate of ITP development may increase depending on the administration method. Table 4 suggests that there is potential for a relatively higher risk for ITP with simultaneous vaccination compared with separate vaccination. However, further analyses using large databases are needed to confirm this finding.

Limitations and strengths

There are several limitations inherent to this study. First, the sample size meant the study was underpowered because we adopted a matched design to adjust for several confounders between cases and controls. This sample size also made it difficult to estimate ORs separately by age group. To compensate for this limitation and present additional results, we calculated common ORs using all available data, although adjustment for confounders might have been insufficient for this analysis. Second, we could not adjust for infection history as a confounder. Although we aimed to measure infection history in the questionnaire, we did not obtain sufficient data. Third, because the investigation depended on the voluntary replies of physicians by post, reporting bias might have occurred. Physicians who were interested in potential adverse events might have been more likely to participate in this study, which could have led to an overestimation of the ORs. A national surveillance system of adverse events following vaccinations in Japan is needed in an era when people are concerned about drug-induced diseases. Fourth, each case was matched to one or two controls. This method might have resulted in an imbalance in the weight of samples for the primary outcome. Because several matching factors may restrict selecting potential control patients, we chose a design to counter this limitation.

In terms of study strengths, the detection of ITP cases by physicians was a main advantage of the present study. The use of physician-registered cases enhanced the internal validity of this investigation. This validity would not be possible in a large database study. Second, the data were gathered from hospitals across Japan. The reported ITP risks, therefore, reflect a broad area, suggesting relatively high external validity. Third, the analyses included detailed results stratified by age group and type of vaccine. This information will be informative for researchers and clinicians seeking to assess the potential risks for adverse events for their patients.

Conclusions

Based on the limited data available in this study, high ITP risk was not found following inactivated or simultaneous vaccination. Future investigations of ITP risk should include analyses using large databases.