INTRODUCTION
Kawasaki disease (KD), which is the most common cause of paediatric vasculitis in developed countries, was first reported in Japan in the 1960s [Reference Kawasaki1]. Following three large epidemics of KD, which occurred in 1979, 1982 and 1986, the incidence of KD has been steadily increasing in Japan [Reference Nakamura2]. Although the aetiology of KD has not been identified, KD is assumed to be infectious based on the presence of seasonality [Reference Burns3] and spatio-temporal clustering during epidemics [Reference Nakamura, Yanagawa and Kawasaki4]. In addition, KD shares many clinical manifestations (e.g. fever, conjunctivitis, rash, cervical adenopathy) with other viral infections [Reference Shingadia, Bose and Booy5]. The rarity of KD in infants aged <6 months indicates the presence of an effective passive immunity [Reference Tsuchida6]. The recurrence of KD is rare [Reference Hirata, Nakamura and Yanagawa7], which suggests that KD may be caused by a limited number of – if not a single – aetiological agents. In recent years, a number of large-scale epidemiological studies have enhanced our understanding of the genetic susceptibility to KD [Reference Burgner8–Reference Onouchi10]. By contrast, however, the aetiological agent(s) of KD has not been identified.
KD can cause a coronary artery lesion and a subsequent fatal myocardial infarction in 1–2% of patients [Reference Kato11, Reference Kato12], unless treated promptly with high-dose intravenous immunoglobulin (IVIG) [Reference Furusho13]; however, 40% of Japanese KD patients who are treated with IVIG therapy still develop transient coronary artery lesions [Reference Ogata14]. IVIG therapy is not only expensive but also harbours the potential for several adverse effects. The identification of the aetiological agent for KD is therefore critical to the development of less expensive and more specific therapies. In the present study, we compared KD with other paediatric infectious diseases with known aetiologies in order to characterize the likely aetiology of KD. A comparison of the epidemiological profiles of multiple diseases can shed light on the determinants of transmission of the individual diseases [Reference Amaku15]. Here, we focused on two epidemiological characteristics, the geographical distribution of the force of infection (FOI) and the periodicity, for KD and each of the seven other paediatric infections.
METHODS
KD data
We obtained the KD data from the Nationwide Surveillance of Kawasaki Disease (NSKD), which includes almost all of the reported cases of KD stratified by age across all 47 prefectures of Japan [16] since 1979 [Reference Nagao17].
Paediatric infectious disease data
The data for nine infectious diseases stratified by age have been reported weekly to the National Epidemiological Surveillance of Infectious Diseases (NESID) from approximately 3000 sentinel paediatric clinics since 1999 (e.g. 2978 clinics in 2000 and 3028 clinics in 2010) [18]. Our analysis excluded infectious gastroenteritis, which is caused by multiple heterogeneous pathogens, and mumps, for which the information about vaccine coverage cannot be obtained. Therefore, we analysed the data for seven paediatric infectious diseases (Table 1). Because the allocation of sentinel clinics is not precisely proportional to the population size, the absolute incidence of disease cannot be estimated from these data. Instead, we estimated the mean patient age from the patient age distribution (see Supplementary Fig. S1).
Table 1. Kawasaki disease and paediatric infectious diseases compared in the present study and their mean patient ages (years) in Japan measured between 2000 and 2010
* Synonym for roseola infantum.
Estimation of the mean patient age
The mean patient age of each infectious disease, including KD, varies by season [Reference Pitzer19]. We pooled the data reported between 2000 and 2010 from both the NSKD and the NESID to circumvent any interference from this seasonal variation. We estimated the crude mean patient age for KD and the seven other paediatric infectious diseases using the following equation [Reference Nagao17]:
$$cA_{ij} \, = \,\displaystyle{{\sum\limits_{k = 1}^{11} {\sum\limits_{m = 2000}^{2010} {\left( {{M}_k \times {N}_{ijkm}} \right)}}} \over {\sum\limits_{k = 1}^{11} {\sum\limits_{ m = 2000}^{2010} {\left( {{N}_{ijkm}} \right)}}}}, $$
where cA ij denotes the crude mean age for the ith disease in the jth prefecture, M k represents the midpoint age for the kth age group, and N ijkm represents the number of patients with the ith disease from the kth age group in the jth prefecture in year m. The definitions of the 11 age groups (Supplementary Fig. S1) were consistent across the diseases and throughout the study period.
To eliminate interference due to the demographic structure of the data, the mean patient age for each disease was adjusted [Reference Nagao17] using the following equation:
$$aA_{ij} \, = \,\displaystyle{{\sum\limits_{k = 1}^{ 11} {\sum\limits_{m = 2000}^{2010} {\left( {{M}_k \times {N}_{ijkm} \times {S}_k { /C}_{\,jkm}} \right)}}} \over {\sum\limits_{k = 1}^{ 11} {\sum\limits_{m = 2000}^{2010} {\left( {{N}_{ijkm} \times {S}_k {/C}_{\,jkm}} \right)}}}}, $$
where aA ij represents the adjusted mean patient age for the ith disease in the jth prefecture, S k represents the proportion of the population in the kth age group in the standard population structure (i.e. Japan's population structure in 2000), and C jkm represents the proportion of the population in the kth age group in the jth prefecture in year m.
In addition, we estimated the national-level mean patient ages using the following equations:
$${\rm national\_}cA_i = \displaystyle{{\sum\limits_{\,j = 1}^{47} \sum\limits_{k = 1}^{11} \sum\limits_{m = 2000}^{2010} {(M_{k} \times} {N}_{ijkm} )} \over {\sum\limits_{{\,j} = 1}^{47} \sum\limits_{{k} = 1}^{11} \sum\limits_{{m} = 2000}^{2010} {({N}_{ijkm} )}}}, $$
$${\rm national\_}aA_i = \displaystyle{{\sum\limits_{j = 1}^{47} \sum\limits_{k = 1}^{11} \sum\limits_{m = 2000}^{2010} {(M_{k} \times} { N}_{ijkm} \times S_{k} /C_{jkm} )} \over {\sum\limits_{{j} = 1}^{47} \sum\limits_{{k} = 1}^{11} \sum\limits_{{m} = 2000}^{2010} {({N}_{ijkm} \times S_{k} /C_{jkm} )}}}, $$
where national_cAi and national_aAi represent the crude and adjusted mean patient ages of the ith disease at the national level, respectively.
Estimation of the FOI
The distribution of the FOI often identifies the risk factors for transmission [Reference Nagao20] and may also reveal the mode of transmission. In the present study, a surrogate measure of the FOI was estimated as the inverse of the mean patient age as shown in the following equation:
$$A = 1/\lambda, $$
where A represents the mean patient age and λ represents the FOI for an infectious agent that manifests as an acute illness and confers lifelong immunity [Reference Anderson and May21]. The FOI estimated by equation (5) has been shown to be a good approximation of more elaborate mathematical techniques used to compute the FOI [Reference Yu22]. The FOI and the mean patient age have been shown to be negatively correlated even in a disease caused by multiple viral strains [Reference Nagao and Koelle23].
Geographical analysis
The geographical distribution of the mean patient age was similar between KD and the other paediatric infectious diseases (Fig. 1). To identify the determinants of this shared geographical distribution, we estimated the correlation between the adjusted mean patient ages and climatic/socioeconomic variables. The following climate data was obtained from the University Corporation of Atmospheric Research [24] as described previously [Reference Nagao17]: mean temperature (°C), precipitation (mm/day), average vapour pressure (mmHg), and the average vapour pressure deficit (which represents aridity, mmHg). Thirty-seven socioeconomic variables were obtained from the Ministry of Internal Affairs and Communications of Japan [25], including seven demographic variables, eight education variables, eight health variables, seven infrastructure variables, four variables related to the standard of living, and three landscape variables (Supplementary Table S1). The values of these variables, which were surveyed at least three times between 2000 and 2010, were averaged for each prefecture. Mitsubishi Tanabe Pharma (Osaka, Japan), the sole retailer of varicella vaccine in Japan, provided the uptake rate of the varicella vaccine at the prefecture level for each year between 2000 and 2010.
Fig. 1. Similarity in the geographical distribution of the crude mean patient age for Kawasaki disease and seven paediatric infectious diseases. The prefectures are categorized from red to blue in ascending order of the crude mean patient age. The distributions of the adjusted mean patient ages were similar to the results presented here. Only data from the main part of each prefecture is presented; minor islands were omitted in this and subsequent figures. Okinawa, which is 650 km away from the main island, is shown close to the main island in this representation. HFMD, Hand, foot and mouth disease; PCF, pharyngoconjunctival fever; GAS, group A streptococcus.
A rank correlation analysis was used to screen the climatic/socioeconomic variables that exhibited a statistically significant (P < 0·05) rank correlation with the adjusted mean patient age of all the diseases. We applied multivariate regression analyses to these preselected variables. The adjusted mean patient age was regressed against the explanatory variables using a conventional multivariate linear regression analysis. Then, to account for spatial autocorrelation, we applied a spatial multivariate regression analysis [Reference Anselin and Hudak26] using the following equation:
$${\bi Y} = {\rho} {\bi UY} + {\bi X}{\bi \beta} + {\rm const}.,$$
where Y represents the dependent variable vector, X represents the independent variable matrix, ρ represents the spatial autoregressive parameter, β denotes the regression coefficient vector, and U represents the spatial adjacency matrix. Each element of U was assigned a value of 1 if the two prefectures were geographically adjacent or connected by a bridge or tunnel, and 0 if the two prefectures were not adjacent or connected. The use of the crude mean patient age (cAij) instead of the adjusted mean patient age (aAij) did not result in qualitatively different results (data not shown). We used Stata v. 11.1 (StataCorp., USA) for all statistical analyses.
Estimation of the periodicity
We also compared the periodicity between KD and other paediatric infectious diseases. The super-annual periodicity (T) in the number of cases was determined by the following equation:
$$T \cong \,2{\rm \pi} [(D + D^{\prime})A]^{{1\over2}}, $$
where D and D’ represent the incubation and infectious periods, respectively, and A represents the mean patient age of the infectious disease [Reference Anderson and May21]. This relationship holds for persistent infections (e.g. Plasmodium falciparum) as well as diseases caused by complex interactions between multiple strains (e.g. dengue haemorrhagic fever) [Reference Hay27]. We tested whether this equation applies to paediatric infectious diseases, and then applied this equation to predict the nature of the aetiological agent for KD.
Time-series analysis
The monthly number of cases of KD in Japan were available from the NSKD. We converted the weekly number of the seven paediatric infectious disease cases, which were reported from the sentinel paediatric clinics to the NESID between 2000 and 2010, into the monthly number of cases. The monthly number of KD cases and the seven other paediatric diseases were transformed as the proportion to the maximum number of cases recorded during the study period, and were subsequently analysed by the wavelet method [Reference Grenfell, Bjornstad and Kappey28] using the Morlet wavelet [Reference Torrence and Compo29]; we followed the standard wavelet procedure used in many recent epidemiological studies (e.g. [Reference Nagao and Koelle23]). We used Mathematica 9.0 (Wolfram, USA) for this analysis.
Estimation of the probability of KD manifestation
We expressed the proportion of individuals who are naive to the aetiological agent for KD as V. The rate of change of V at the age of t years can be expressed as:
$${\rm d}V = -{\rm \lambda} V{\rm d}t,$$
where λ represents the FOI. Equation (5) holds for a persistent infection if the incubation period is not very long. Therefore, we solved equation (8) using equation (5) and a condition of V(0) = 1 into the following equation:
$$V = \exp ( - {\rm \lambda} t)\, = \exp( - t/A),$$
where A represents the mean patient age. Based on this equation and epidemiological and demographic data from Japan, we estimated the probability of developing KD when infected with its aetiological agent.
RESULTS
Geographical analysis
The distribution of the mean patient age was similar between KD and the other paediatric infectious diseases (Fig. 1); there was a statistically significant rank correlation between all pairs of these diseases with respect to the prefectural mean patient age (Supplementary Table S2). The rank correlations between the adjusted mean patient ages and the climatic/socioeconomic variables were estimated (Supplementary Table S3). This screening process revealed eight variables that were significantly correlated with the adjusted mean patient ages of all the diseases (Table 2). Four of these variables were related to healthcare, and hence were mutually correlated. We excluded non-unique variables to avoid multi-collinearity effects (see uniqueness in Table 2), and included the varicella vaccine uptake rate as a covariate. The final statistical model for each disease was determined by a step-wise addition and elimination of these variables (Table 3).
Table 2. Variables which exhibited statistically significant rank correlations with the adjusted mean patient age of Kawasaki disease and paediatric infectious diseases (n = 47)
HFMD, Hand, foot and mouth disease; PCF, pharyngoconjunctival fever; GAS, group A streptococcus.
* Uniqueness was estimated by factor analysis. A variable of smaller ‘uniqueness’ can be largely explained by a linear combination of the other variables.
Table 3. Conventional and spatial regression models to explain adjusted mean patient ages of Kawasaki disease and paediatric infectious diseases
HFMD, Hand, foot and mouth disease; PCF, pharyngoconjunctival fever; GAS, group A streptococcus; n.s., not significant.
* Although five variables (total fertility rate, pupils per class, health insurance, waste process, and chickenpox vaccination) were incorporated into the multivariate analysis, only total fertility rate, health insurance, and chickenpox vaccination remained as statistically significant contributor(s).
The health insurance paid per insuree remained as a significant negative contributor to the adjusted mean patient age for most of the diseases. The total fertility rate (TFR) remained as a significant negative contributor to the adjusted mean patient age for all the diseases except group A streptococcus (GAS). These results suggest that the availability/utilization of health services, and more importantly, the TFR (Fig. 2) are critical factors in the transmission of KD and paediatric infectious diseases. Notably, Okinawa prefecture, which has the largest TFR and lowest health insurance paid per insuree, emerged as an outlier for some diseases (Supplementary Figs S2, S3).
Fig. 2. Distribution of socioeconomic factors that were correlated with the mean patient age for each prefecture. The prefectures were categorized from red to blue in descending order for (a) the total fertility rate (TFR) and (b) the health insurance paid per insuree.
Seasonality
We compared the seasonality across diseases. We presented our data starting in April to reflect the beginning of the Japanese academic school year (Fig. 3). Hand, foot, and mouth disease (HFMD) and herpangina exhibited a unimodal distribution with a sharp peak during the summer (July). GAS exhibited trimodal seasonality with three major peaks in June, December, and March. KD, chickenpox, pharyngoconjunctival fever (PCF), and erythema infectiosum exhibited a bimodal distribution with two major peaks, one in June and the other in December/January, as well as a subtle minor peak in March.
Fig. 3. Seasonality of Kawasaki disease and other paediatric infectious diseases. The proportion of the number of cases in a month to the annual number of cases was estimated for each year between 2000 and 2010 and averaged across these 11 years. The dashed lines indicate the standard deviation. The x-axis starts in April to reflect the beginning of the Japanese academic school year. HFMD, Hand, foot and mouth disease; PCF, pharyngoconjunctival fever; GAS, group A streptococcus.
Super-annual periodicity
All the diseases examined in the present study, except KD and exanthema subitum, exhibited substantial fluctuations in the number of cases from year to year (Fig. 4). The wavelet analysis revealed that all the diseases except KD and exanthema subitum exhibited some coloured area, which implies that super-annual periodicity was detected in the corresponding period (Fig. 5). We obtained the incubation period (D) and infectious period (D’) of each disease from Richardson et al. [Reference Richardson30] (Table 4) and used the adjusted mean patient age (A) (Table 1) to test whether the super-annual periodicity could be explained by equation (7). Most of the super-annual periodicities estimated by the wavelet analysis were consistent with the predicted values (Table 4). We could not judge whether the prediction for GAS was successful because D + D’ was ambiguous for this disease. We also predicted D + D’ for KD using equation (7). Because T was longer than 10 years after 1987 (Fig. 6) and A is 2·4 years, the D + D’ for KD was estimated to be larger than 1·1 years. Since D for KD is approximately 10 days [Reference Fujita31], D’ for KD therefore is likely to be longer than 1 year. We will report the hypothesis elsewhere that the epidemics in 1979, 1982, and 1986 were driven by dynamics unique to these years (Y. Nagao, unpublished data ).
Fig. 4. Annual time series of Kawasaki disease (KD) and other paediatric infectious diseases in Japan between 2000 and 2010. The annual number of cases of KD and other paediatric infectious diseases was expressed as a proportion to the maximum number of reported cases. PCF, pharyngoconjunctival fever; GAS, group A streptococcus.
Fig. 5. Super-annual periodicities of KD and other paediatric infectious diseases as detected by the wavelet analysis. The blue lines indicate the 95% confidence limit of the wavelet analysis. The area outside the blue lines is unreliable due to the edge effect. PCF, pharyngoconjunctival fever; GAS, group A streptococcus.
Fig. 6. Monthly number of Kawasaki disease (KD) cases between 1979 and 2010. The monthly number of KD cases recorded between 1979 and 2010 in Japan are expressed as the proportion to the maximum number of cases (a). The results of the wavelet analysis applied to this time series are plotted in (b).
Table 4. The lengths of predicted and actual super-annual periodicities (T)
HFMD, Hand, foot and mouth disease; PCF, pharyngoconjunctival fever; GAS, group A streptococcus; d, days; wk, weeks; mo., months; yr, years.
Probability of KD manifestation
Between 2006 and 2010, the average annual number of KD cases in Japan was 7699 for children aged between 0 and 2 years and 10 115 for children aged between 0 and 4 years. During the same period, the average number of live births was 1 082 997 per year and the death rate for children aged <5 years was 3667 per year. From these numbers, we estimated that approximately 0·7% and 0·9% of Japanese children had been diagnosed with KD by their third and fifth birthdays, respectively. By substituting A in equation (9) with the mean age of KD (i.e. 2·4 years), the proportion of individuals who experienced infection(s) with the aetiological agent for KD (1 – V) was computed as 71% and 88% at t = 3 and t = 5 years, respectively. Therefore, only about 1% of infections with the KD agent(s) manifested with symptoms of KD.
DISCUSSION
The geographical distribution of the FOI of KD and other paediatric infectious diseases was similar across Japan. We identified the availability/utilization of health services and the TFR as the main determinants of this shared geographical structure. Particularly, the TFR remained as a statistically significant risk factor to the FOI of KD and all the viral diseases, after accounting for the possible confounders and the spatial autocorrelation.
Indeed, there have been numerous reports that health services facilitate the transmission of paediatric infectious diseases [Reference Johnston and Burke32, Reference Buettcher and Heininger33]. In addition, the causal relationship between the TFR and the FOI is intuitive: a higher TFR implies a larger number of siblings within a household, which increases the probability of sibling-to-sibling infection. Previous studies have also identified the importance of sibling-to-sibling infections for KD [Reference Fujita31, Reference Tsai34].
The finding that sibling-to-sibling infection and, to a lesser extent, hospital infection emerged as important routes of transmission for KD and other paediatric infectious diseases in this study suggests that KD shares a mode of transmission with other paediatric infectious diseases, i.e. close contact with an infected individual. KD also exhibited a bimodal seasonality similar to other viral infections such as chickenpox and erythema infectiosum, each of which is caused by a single agent. Therefore, we cannot assert that KD is caused by multiple agents based solely on its multimodal seasonality. The seasonality of a paediatric infectious disease is largely determined by climate and the school term [Reference Brisson35, Reference Grassly and Fraser36]. Indeed, the troughs in the seasonality of GAS occurred in April, September, and January, immediately following school/kindergarten holidays in the spring (late March to early April), the summer (late July to late August), and the winter (late December to early January), respectively. This finding suggests that these troughs were due to the separation of children during the holidays. This finding also suggests that the bimodal seasonalities of KD, chickenpox, PCF, and erythema infectiosum were affected by the school term. Our finding that KD shares a similar seasonality pattern to other paediatric diseases suggests that the transmission routes for these diseases are common. Taken together, these results suggest that KD is not likely to be transmitted through a novel route such as a transcontinental wind [Reference Rodo37, Reference Rodo38].
One notable outlier in our regression analysis was Okinawa prefecture, a tropical prefecture composed of 160 islands located 650 km from Japan's mainland. Host genetic heterogeneity is an unlikely explanation for this observation, because the genetic structure of the Okinawan population is similar to other populations within Japan [Reference Yamaguchi-Kabata39]. The tropical climate of Okinawa may contribute to this observation [Reference Burns40]. Alternatively, due to the remoteness and multi-island structure of Okinawa, the transient extinction of a pathogen may be more frequent than on the mainland of Japan [Reference Nodelijk41].
One major limitation of the present study was that the data about the paediatric infectious diseases was derived from sentinel paediatric clinics. The allocation of sentinel clinics is not precisely proportional to the population size. However, no factor was suggested as a possible source of bias in estimating the mean patient age from the sentinel data. Therefore, we believe that our analyses were not affected significantly by the use of sentinel data.
Another possible refutation against the present analysis is that our statistical analysis may have been flawed by ecological confounding. However, we believe that the identified contributions of fertility rate and hospital utilization to the FOI are not statistical artefacts. Four of the eight variables, which were preselected from a large number of variables during the initial screening process, were related to the utilization of health services (Table 2). This result strongly suggests that health services play an important role in determining the FOI. In addition, the varicella vaccine uptake rate remained as a statistically significant predictor of risk reduction on the FOI of chickenpox only (Table 3). This result validates our analytical procedure, because vaccination is assumed to decrease the FOI [Reference Anderson and May21, Reference Laurenson42]. Moreover, the fertility rate (or birth rate) is one of the most important factors that determine the epidemiology of paediatric infectious diseases [Reference Grenfell, Bjornstad and Kappey28].
Our mathematical assumptions may seem to contain excessive approximations. For example, equation (5) presumes a fixed longevity (or Type I survival). This, however, is a well-founded approximation. First, Type I survival is more suitable ‘for humans, especially in developed countries’ than Type II survival, which assumes age-independent mortality [Reference Anderson and May21]. Second, regardless of whether Type I or II survival is assumed, the mean patient age can effectively approximate the inverse of the FOI, if the mortality is very small as in Japan. Equation (8), which has been used in epidemiological studies [Reference Thai43], presumes continuous susceptibility to a given disease from birth and complete endemicity of the disease agent. Therefore, our prediction based upon equation (8) should be tested in future serological surveys once the aetiological agent for KD is identified.
Given that any infectious disease has finite values for the incubation and infectious periods (i.e. D + D’), equation (7) predicts that any infectious disease will have a finite super-annual periodicity (T). However, among the diseases examined in this study, only KD and exanthema subitum exhibited the least evidence for super-annual periodicities. This finding suggests that the super-annual periodicities of these diseases are very long, and hence the sum of the incubation period and infectious period of these diseases is large [equation (7)]. Consistent with this finding, the viral agent for exanthema subitum persists across the lifespan [Reference Suga44]. Therefore, it would not be surprising if the aetiological agent for KD also persists for a long period. Our finding is consistent with a previous report that the agent (or at least one of the agents) for KD is capable of persistent infection [Reference Pitzer19].
The low mean patient age, and therefore the strong FOI of KD, led to an estimation that only about 1% of infections with the KD agent(s) manifest with symptoms of KD. Our finding that KD may persist over a long period and cause frequent asymptomatic infections predicts that the aetiological agent(s) for KD exists in KD patients as well as in many healthy individuals. Therefore, an epidemiological study to determine whether a particular microbe is a causative agent for KD will require a large sample size to detect subtle difference in the prevalence of the microbe between cases and controls. In addition, the possibility that sibling-to-sibling infection is an important route of transmission for KD implies that the number of siblings and the rank of the patient among them should be considered as covariates in epidemiological analyses.
SUPPLEMENTARY MATERIAL
For supplementary material accompanying this paper visit http://dx.doi.org/10.1017/S0950268815001223.
ACKNOWLEDGEMENTS
We are grateful to Yoshikazu Nakamura and Mayumi Yashiro for providing the KD data. We are also grateful to Maurizio Pisati and Junji Shiraishi, who assisted in the spatial analysis. We thank Tetsuya Sano for helpful comments on the manuscript.
This research was supported in part by the Aihara Project, the FIRST programme from the Japan Society for the Promotion of Science, which was initiated by the Council for Science, Technology, and Innovation.
DECLARATION OF INTEREST
None.
Open access
