- Research article
- Open Access
- Open Peer Review
Analysis of influenza transmission in the households of primary and junior high school students during the 2012–13 influenza season in Odate, Japan
BMC Infectious Diseasesvolume 15, Article number: 282 (2015)
Households are one of the major settings of influenza transmission in the community and transmission is frequently initiated by school-aged children. We surveyed households with primary school (PS) and/ or junior high school (JH) children for the 2012–13 influenza season in Odate, Japan then characterized the epidemiology of influenza household transmission as well as estimated the serial intervals.
We delivered a self-reported questionnaire survey to households with PS and/or JH school children in Odate City, Japan. Influenza A (H3N2) virus predominantly circulated during the 2012–13 influenza season. We investigated the epidemiological characteristics of within-household transmission and calculated the serial intervals (SI). SIs were drew by a non-parametric model and compared with parametric models by the Akaike Information Criterion. The covariable contributions were investigated by the accelerated failure model.
Household influenza transmission was identified in 255 out of 363 household respondents. Primary school (PS) children accounted for 45.1 % of primary cases, and disease transmission was most commonly observed between PS children and parents, followed by transmission from PS children to siblings. In primary cases of PS or JH children, younger age and longer absence from school were significantly associated with household transmission events. The mean SI was estimated as 2.8 days (95 % confidence interval 2.6-3.0 days) in the lognormal model. The estimated acceleration factors revealed that while secondary school age and the absence duration > 7 days were associated with shorter and longer SIs, respectively, antiviral prescriptions for primary cases made no contribution.
High frequencies of household transmission from primary school with shorter SI were found. These findings contribute to the development of future mitigation strategies against influenza transmission in Japan.
Influenza is an acute viral respiratory disease, which is usually self-limiting but can lead to severe complications and reach pandemic proportions. The influenza pandemic of 2009 affected an estimated 24 % of the global population  with over 100,000 deaths worldwide . Various non-pharmaceutical measures, such as hand hygiene and school closures, have been implemented, not only during pandemics but also during the periods of seasonal influenza . To optimise these preventive measures, we must understand the dynamics of influenza transmission. Households provide excellent environments for influenza transmission, as contact among household members is exceptionally high . Although school-aged children play an important role in influenza transmission , the rate of household transmission is affected by factors such as the structure and the size of the household , pre-existing immunity  and the household environment. The situation in Japan is rather unique partly because class or grade dismissal is commonly implemented during seasonal influenza epidemics in schools and partly because the proportion of the elderly has reached up to 23 % of the total population in 2010 . Although several studies have investigated the epidemiology of influenza transmission in households during the 2009 H1N1 influenza pandemic [9, 10] and the influenza A(H2N2) of 1957 pandemic  in Japan, studies on seasonal influenza have been rarely reported.
Household studies usually measure the secondary infection risk; typically the secondary attack rate (SAR). This proportional parameter enables us to estimate the number of subsequent cases and evaluate the risk ratio, especially by age group. Another important parameter is the serial interval (SI), defined as the time between the onset of specific signs and the symptoms of a primary case and the onset of a secondary case [12, 13]. This parameter characterises the infectiousness profile and determines the rate of epidemic growth. By analysing these variables, we can elucidate (at least partially) the influenza transmission dynamics .
We conducted an enhanced influenza surveillance in Odate City, Akita Prefecture, Japan, during the 2011–12 season and found nearly 95 % of patients came to health facilities within 2 days from their onsets . We also conducted a cross-sectional survey among households with primary school (PS) or junior high (JH) school children during that season. This study aimed to characterise the influenza epidemiology in the households of a rural city in Japan and to estimate the SI in this setting.
Study site and design
Odate City is located in Akita Prefecture, northeast of Honshu, Japan. Among its population of 78,946, 15.6 and 31.7 % were registered as under 20 years old and over 64, respectively, in the national 2010 census. A questionnaire survey, designed to obtain information on household influenza transmission, was administered to households with PS or JH children during the 2012–13 influenza season. The questionnaire was handed out to students two times, in January and in March, and collected altogether at the end of March 2013. A total of 5,225 questionnaires were distributed at each time of the survey. Respondents were instructed to record influenza episodes for each household member. The questionnaire requested the demographic information of the influenza cases, students’ grade, the date of fever onset, the status of influenza vaccination, any medical consultation during the clinical course of the disease, duration of school absence and the administration of antiviral drugs. We did not collect the total size of respondents’ households.
Influenza surveillance during the 2012–2013 season
Influenza surveillance has been undertaken since the 2011–12 season in the city . Throughout the surveillance, an influenza case was suspected if the subject developed fever with cough, runny nose or congested nose. All suspected cases were tested using a commercial rapid test kit. From December 2012 to March 2013, 1,806 suspected influenza cases tested positive, 99.4 % of which were positive for influenza A. The number of PS and JH children cases detected through the surveillance was 352 while the number reported in the survey was 384. The daily numbers of cases reported from the surveillance and the survey were plotted and the correlation was examined by the cross-correlation test. Unfortunately, we could not totally link cases between the surveillance and the survey. In addition, as per a surveillance report of national infectious diseases, 38 out of 40 influenza A strains in Akita Prefecture during that season were influenza A(H3N2); the remaining two were influenza A(H1N1)pdm09 . Only three strains of influenza B were isolated in the prefecture, all belonging to the influenza B Yamagata lineage.
In the survey, an influenza case was defined as those who reported influenza episodes with a history of medical consultation. A primary case is defined as the first influenza case within a household, whereas a secondary case is any influenza case following the primary case in the same household. Episodes whose intervals were either zero or more than seven days after the onset of the primary case were excluded from the analysis. Family members were categorised into seven groups: preschool (PreS), PS, JH school, high school (HS), father, mother and grandparents/other adults. The antiviral treatment history, whether subjects had received Oseltamivir, Zanamivir, Laninamivir or Peramivir, was also recorded. Then we characterised the epidemiological factors of primary PS and JH cases between households with and without household transmission events. These variables were the month of onset, the age of primary case, antiviral treatment, pre-season influenza vaccination and duration of school absence of primary cases.
The SI was defined as the interval between the onset of symptoms in a primary and secondary case. Generally SI data constitute a form of time-to-event data and thus it was necessary to consider the censoring time. Unlike in previous studies [13, 17], the questionnaire was distributed through PS or JH children twice in the study season, potentially minimising the data truncation. The SI estimation model was applied from a Hong Kong study . Briefly, the SIs were visualised by a non-parametric model called the Kaplan–Meier estimator, which was compared with the parametric Weibull, gamma and lognormal distributions. Ninety-five percent (95 %) confidence interval (CI) for each model were calculated by using a parametric bootstrap approach with 1,000 resamples. The best-fit parametric model was decided by the Akaike Information Criterion (AIC). Households in which the primary case was a PS or JH student were subjected to multi-variable analysis using an accelerated failure time model to evaluate the contribution of covariates to the SI. The model included the age groups of primary and secondary cases, antiviral prescription, pre-season influenza vaccination and absence duration of primary cases. The transformed regression coefficients derived from the model were interpreted as the acceleration factor (AF). This factor was considered a multiplicative increase (if >1) or decrease (if <1) of the median SI relative to the reference variable.
Categorical data was analysed either by Wilcoxon test or Kurskal-Wallis test and continuous variables were analysed by the chi-square test. The level of statistical significance was determined at p <0.05. All statistical analyses were conducted by R 3.1.0 .
All participants were provided with an instruction document stating our research purpose and describing the survey method, along with an informed consent form. Participants completed the questionnaire on agreement and returned that with the written consent. The entire study design was approved by the Ethics Committee of Tohoku University Graduate School of Medicine (ID 2011–268).
Among the 2,930 identified eligible households, 363 responded to the survey (12.4 %). Of the responding households, 356 (98.1 %) reported at least one influenza case, and 589 influenza cases were recorded in total. The daily number of surveyed influenza cases was highly correlated with the number of influenza cases detected in an enhanced influenza surveillance of the city  (Fig. 1) (correlation coefficient = 0.877). Forty-three percent of the surveyed households reported more than two influenza cases (Table 1). Throughout the study period, 93 households reported two cases and 4 households reported five cases.
Junior high (JH) and PS students accounted for the majority of cases (71.1 %), but a substantial number of cases developed in mothers, pre-school children and fathers (9.2 %, 6.5 % and 6.1 % respectively). The age distribution of cases reflects their status category (Table 2). Parents, grandparents and other adults accounted for 7.3 % of primary cases and 42.2 % of secondary cases (p < 0.0001). More secondary cases than primary cases were also reported among PreS children.
Among those households reporting influenza cases, 255 reported transmission episodes within the household (Additional file 1: Table S1). Household transmissions from PS, JH and PreS children accounted for 45.1 %, 17.2 % and 13.3 % of all episodes, respectively. Transmissions from PS to mother were most frequently observed (30 episodes), followed by transmissions from PS to father (20 episodes), and from PS to JH or PS to PS (18 episodes) (Fig. 2). The grade distributions were not significantly different between primary and secondary PS cases (p = 0.58), although 38 % of the total episodes developed in grade 1 and grade 2 students. Within households, JH cases were most commonly transmitted to PS children, followed by mothers and grandparents/other adults. Stratified by grades, no significant distribution differences appeared among episodes (p = 0.75), but 60 % of the total episodes developed in JH grade 1 students. Transmissions between siblings, such as PreS to PS and JH to PreS, were also observed. Although the number was small, all of the HS cases were primary cases in household transmissions, except one transmission event from PS to HS. Among parent–parent transmissions, the father was more frequently the primary source than the mother (Fig. 2a, Additional file 1: Table S1).
Table 3 summarises the epidemiological factors observed in primary PS and JH cases between households with and without transmission events. Among the explanatory factors, younger age and longer absence interval in primary cases were significantly associated with household transmission events while the percentages of antiviral drugs prescription as well as influenza vaccination prior to the season were not significantly different.
The distribution of SIs was constructed for the 298 household transmissions (Fig. 2b) and fitted with parametric models (Fig. 3). These parametric models yielded a mean SI of 2.8 days (95 % CI 2.6-3.0), 3.1 days (95 % CI 2.9-3.3), and 2.8 days (95%CI 2.6-3.0) in Weibull, gamma, and lognormal distributions, respectively. The lowest AIC value was obtained in the lognormal model; hence, this model was used in further analysis. We then performed multi-variable analysis among households whose primary cases were PS or JH students. Specifically, we evaluated the covariate contributions to the SI and the obtained AF using accelerated failure time models. The results are presented in Table 4. Longer SI was more closely associated with JH cases than with PS cases. A significantly shorter SI was associated with secondary cases in both PS and JH children, using fathers as reference. The SI was slightly extended in secondary PreS cases, although the result is not statistically significant. Both antiviral prescription and influenza vaccination shortened the SI in primary cases but not to a statistically significant extent. On the other hand, the duration of absence > 7 days in primary cases was significantly associated with a longer SI.
Numerous studies have investigated household influenza transmission during periods of seasonal influenza [19, 20] and during the 2009 influenza A(H1N1) pandemic [21–24] and determined that households were a major source of influenza transmission in communities [25, 26]. This study was based on a questionnaire survey administered to households through primary school or junior high school children. From the results, we characterised the household transmission dynamics during the 2012–13 influenza season in Odate City, Japan. During this season, approximately 40 % of responding households reported ≥ 1 member affected by the disease.
Children in pre-school (PreS), primary school (PS), junior high school (JH) and high school (HS) collectively initiated 78 % of household transmissions, of which 80 % were sourced from PS or JH children. We found more household transmission events in younger grades of PS or JH children compared with older counterparts and the regression analysis revealed the odds of household influenza transmission declined with age. One reason is that younger children transmit influenza more efficiently due to closer contacts with other family members . Another reason is that the spread of influenza was initiated in those grades and consequently introduced to the household through those grade groups  as the probability of infection in different grades was smaller than that within the same grade .
In the present study, secondary cases were most common among the age groups of 10–19 and 20–49 years. Similar trends were reported in other studies. The exception was a U.S. study, in which secondary cases were most common in children aged below 9 years . We determined that 30 and 38 % of cases were transmitted from children to parents and from sibling to sibling, respectively. Intense contact between same age groups and between children and parents has been documented in studies of social mixing patterns [4, 29]. These data suggest that household transmission can be prevented by decreasing the contact intensity through measures such as home quarantine . In addition, longer school absence duration was associated with longer SIs. In Japan, school children with influenza are required by law to stay home until 5 days after the onset and 2 days after the resolution of fever. We did not obtain the timing of symptom resolution; however, the aforementioned association implies that prolonging the time to disease resolution increases the opportunity for virus exposure for a longer time. Carers of infected children need to adopt personal hygiene procedures to minimise disease transmission.
The percentages of reported vaccination as well as antiviral administration of the primary case were not significantly different between households with and without transmission events. Previous studies determined that vaccination exerts no significant effect on secondary infection risk [5, 19]. Similarly, seasonal influenza studies have determined that antiviral treatment imposes no significant reduction in secondary infection risk among household contacts [31, 32]. On the other hand, influenza A(H1N1)pdm09 outbreaks have been significantly  or non-significantly  reduced by antivirals. As previously reported , nearly 95 % of the medical consultations in the Odate population were made within two days of symptom onset; thus, we expected an early antiviral administration could have a protective effect on secondary transmission. However, the extent of protection was insufficient to prevent the transmission. The mean SI of 2.8 days in our study supports a scenario that an exposure to influenza virus is likely to occur prior to antiviral administration. A shorter SI was associated with antiviral administration on a primary case without statistically significant. This suggests the possibility to inhibit the exposure of influenza in relatively late timing by the drug. The infectious period for influenza has not been elucidated, but cohort studies have revealed that viral shedding peaks around the time of symptom appearance  and reduces over time . The data of antiviral administration timing can give an opportunity of further analysis of its effect on the SI.
The SI depends on the setting as well as on the causative agent. As per previous studies, the mean SI for influenza A(H3N2) ranges between 1.9 and 3.3 days inclusive [17, 20, 22, 35] and the mean generation time was computed as 3.1 days . Our estimated SI, as well as mean generation time, was comparable to that in studies in Thailand  and Hong Kong  while another two studies [20, 35] estimated a shorter SI. These differences could be attributed to the right censoring of data  or the strict inclusion criteria of secondary contacts among households  in the other studies, and to the variation of the sizes of the enrolled households among studies . However, the compatible estimation of SI in our study adds to the evidence that an influenza virus transmission to the next generation occurs within 3 days in household settings.
The AF in secondary cases with fathers as reference indicated that SIs are reduced among children and those shorter SIs were statistically significant among PS and JH children. These findings, which are partially compatible with a previous study , suggest that these school age groups can transmit more efficiently because of frequent contact or physical proximity to the primary case in households.
Our study has several limitations. First, since the study questionnaires were distributed through PS and JH students, households without children in these groups were excluded from the analysis. Therefore, we could not investigate the epidemiology of household transmission in general populations. In fact, we conducted an intensive study on households with PS and JH children. Second, we largely lacked information on influenza-free households which consequently led to the low response rate of the survey. Third, we may have underestimated the number of household transmission episodes since some households did not respond to the survey in spite of the instruction. Also, asymptomatic cases can alter the number of secondary influenza episodes in a population; indeed, one study found a substantial proportion of asymptomatic influenza A(H3N2) infections in households . Fourth, we did not obtain the sizes of respondents’ households and thus we could not calculate the SAR as mentioned. Fifth, we relied on households to self-report their medical consultation histories without confirming a laboratory-tested positive influenza result. In Japan, most patients who visit medical facilities with suspected influenza are routinely examined with a commercial rapid test kit . We included the history of medical consultation in our case definition to increase the sensitivity. Sixth, because we collected the questionnaires at the end of March, we missed any influenza cases that occurred after the study period. Finally, we excluded cases whose onsets coincided with the primary onset, potentially missing some transmissions.
Despite these limitations, we analysed a substantial number of household transmission events in a rural city of Japan and identified school children as the major initial source of infection and estimated a next generation episode within a mean of 2.8 days during influenza A (H3N2) epidemics in the 2012–13 season. Frequent acquirers of familial influenza infection were parents and siblings. A household is considered to provide an opportunity to transmit influenza across different generations. Our findings contribute to the development of future mitigation strategies against influenza transmission in Japan.
Secondary attack rate
Junior high school
Akaike information criterion
Van Kerkhove MD, Hirve S, Koukounari A, Mounts AW. Estimating age-specific cumulative incidence for the 2009 influenza pandemic: a meta-analysis of A(H1N1)pdm09 serological studies from 19 countries. Influenza Other Respir Viruses. 2013;7(5):872–86.
Simonsen L, Spreeuwenberg P, Lustig R, Taylor RJ, Fleming DM, Kroneman M, et al. Global mortality estimates for the 2009 Influenza Pandemic from the GLaMOR project: a modeling study. PLoS Med. 2013;10(11):e1001558.
World Health Organization Writing G, Bell D, Nicoll A, Fukuda K, Horby P, Monto A, et al. Non-pharmaceutical interventions for pandemic influenza, national and community measures. Emerg Infect Dis. 2006;12(1):88–94.
Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008;5(3):e74.
Azman AS, Stark JH, Althouse BM, Vukotich Jr CJ, Stebbins S, Burke DS, et al. Household transmission of influenza A and B in a school-based study of non-pharmaceutical interventions. Epidemics. 2013;5(4):181–6.
Marathe A, Lewis B, Chen J, Eubank S. Sensitivity of household transmission to household contact structure and size. PLoS One. 2011;6(8):e22461.
Klick B, Nishiura H, Ng S, Fang VJ, Leung GM, Peiris JS, et al. Transmissibility of seasonal and pandemic influenza in a cohort of households in Hong Kong in 2009. Epidemiology. 2011;22(6):793–6.
Muramatsu N, Akiyama H. Japan: super-aging society preparing for the future. Gerontologist. 2011;51(4):425–32.
Hirotsu N, Wada K, Oshitani H. Risk factors of household transmission of pandemic (H1N1) 2009 among patients treated with antivirals: a prospective study at a primary clinic in Japan. PLoS One. 2012;7(2):e31519.
Komiya N, Gu Y, Kamiya H, Yahata Y, Yasui Y, Taniguchi K, et al. Household transmission of pandemic 2009 influenza A (H1N1) virus in Osaka, Japan in May 2009. J Infect. 2010;61(4):284–8.
Nishiura H, Chowell G. Household and community transmission of the Asian influenza A (H2N2) and influenza B viruses in 1957 and 1961. Southeast Asian J Trop Med Public Health. 2007;38(6):1075–83.
Fine PE. The interval between successive cases of an infectious disease. Am J Epidemiol. 2003;158(11):1039–47.
Cowling BJ, Fang VJ, Riley S, Malik Peiris JS, Leung GM. Estimation of the serial interval of influenza. Epidemiology. 2009;20(3):344–7.
Keeling MJ, Danon L. Mathematical modelling of infectious diseases. Br Med Bull. 2009;92:33–42.
Huo X, Kamigaki T, Mimura S, Takahashi Y, Oshitani H. Analysis of medical consultation interval between the symptom onset and consultation observed in multiple medical facilities in Odate city, Japan, 2011/2012 and 2012/2013 seasons. J Infect Chemother. 2014;20(6):370–4.
National Institute of Infectious Diseases. Influenza 2012/13 season. In: Infectious Agents Surveillance Report.2013. http://www.nih.go.jp/niid/ja/flu-m/flu-iasrtpc/4100-tpc405-j.html (in Japanese). Accessed 18 Sep 2014
Levy JW, Cowling BJ, Simmerman JM, Olsen SJ, Fang VJ, Suntarattiwong P, et al. The serial intervals of seasonal and pandemic influenza viruses in households in Bangkok. Thailand Am J Epidemiol. 2013;177(12):1443–51.
R Development Core Team. R: A language and environment for statistical computing. In: editor. R Development Core Team. Vienna, Austria: R Foundation for Statistical Computing; 2014.
Petrie JG, Ohmit SE, Cowling BJ, Johnson E, Cross RT, Malosh RE, et al. Influenza transmission in a cohort of households with children: 2010–2011. PLoS One. 2013;8(9):e75339.
Viboud C, Boelle PY, Cauchemez S, Lavenu A, Valleron AJ, Flahault A, et al. Risk factors of influenza transmission in households. Br J Gen Pract. 2004;54(506):684–9.
Kim CY, Breiman RF, Cosmas L, Audi A, Aura B, Bigogo G, et al. Secondary household transmission of 2009 pandemic influenza A (H1N1) virus among an urban and rural population in Kenya, 2009–2010. PLoS One. 2012;7(6):e38166.
Cowling BJ, Chan KH, Fang VJ, Lau LL, So HC, Fung RO, et al. Comparative epidemiology of pandemic and seasonal influenza A in households. N Engl J Med. 2010;362(23):2175–84.
Loustalot F, Silk BJ, Gaither A, Shim T, Lamias M, Dawood F, et al. Household transmission of 2009 pandemic influenza A (H1N1) and nonpharmaceutical interventions among households of high school students in San Antonio, Texas. Clin Infect Dis. 2011;52 Suppl 1:S146–153.
Cauchemez S, Donnelly CA, Reed C, Ghani AC, Fraser C, Kent CK, et al. Household transmission of 2009 pandemic influenza A (H1N1) virus in the United States. N Engl J Med. 2009;361(27):2619–27.
Longini Jr IM, Koopman JS, Monto AS, Fox JP. Estimating household and community transmission parameters for influenza. Am J Epidemiol. 1982;115(5):736–51.
Monto AS. Studies of the community and family: acute respiratory illness and infection. Epidemiol Rev. 1994;16(2):351–73.
Kucharski AJ, Kwok KO, Wei VW, Cowling BJ, Read JM, Lessler J, et al. The contribution of social behaviour to the transmission of influenza A in a human population. PLoS Pathog. 2014;10(6):e1004206.
Cauchemez S, Bhattarai A, Marchbanks TL, Fagan RP, Ostroff S, Ferguson NM, et al. Role of social networks in shaping disease transmission during a community outbreak of 2009 H1N1 pandemic influenza. Proc Natl Acad Sci U S A. 2011;108(7):2825–30.
Horby P, Pham QT, Hens N, Nguyen TT, Le QM, Dang DT, et al. Social contact patterns in Vietnam and implications for the control of infectious diseases. PLoS One. 2011;6(2):e16965.
Zhang D, Liu W, Yang P, Zhang Y, Li X, Germ KE, et al. Factors associated with household transmission of pandemic (H1N1) 2009 among self-quarantined patients in Beijing. China PLoS One. 2013;8(10):e77873.
Ng S, Cowling BJ, Fang VJ, Chan KH, Ip DK, Cheng CK, et al. Effects of oseltamivir treatment on duration of clinical illness and viral shedding and household transmission of influenza virus. Clin Infect Dis. 2010;50(5):707–14.
Carrat F, Duval X, Tubach F, Mosnier A, Van der Werf S, Tibi A, et al. Effect of oseltamivir, zanamivir or oseltamivir-zanamivir combination treatments on transmission of influenza in households. Antivir Ther. 2012;17(6):1085–90.
Goldstein E, Cowling BJ, O’Hagan JJ, Danon L, Fang VJ, Hagy A, et al. Oseltamivir for treatment and prevention of pandemic influenza A/H1N1 virus infection in households, Milwaukee, 2009. BMC Infect Dis. 2010;10:211.
Loeb M, Singh PK, Fox J, Russell ML, Pabbaraju K, Zarra D, et al. Longitudinal study of influenza molecular viral shedding in Hutterite communities. J Infect Dis. 2012;206(7):1078–84.
Suess T, Remschmidt C, Schink SB, Schweiger B, Heider A, Milde J, et al. Comparison of shedding characteristics of seasonal influenza virus (sub)types and influenza A(H1N1)pdm09; Germany, 2007–2011. PLoS One. 2012;7(12):e51653.
Carrat F, Vergu E, Ferguson NM, Lemaitre M, Cauchemez S, Leach S, et al. Time lines of infection and disease in human influenza: a review of volunteer challenge studies. Am J Epidemiol. 2008;167(7):775–85.
Black AJ, Ross JV. Estimating a Markovian epidemic model using household serial interval data from the early phase of an epidemic. PLoS One. 2013;8(8):e73420.
We would like to thank the students, students’ parents and teachers in schools who collaborated for the questionnaire survey. We also appreciate the assistance of the Bureau of Education, Odate City. The study was financially supported by the Health and Labour Sciences Research Grants for the Research on Emerging and Re-emerging Infectious Diseases, Ministry of Health, Labour and Welfare, Japan and JSPS KAKENHI Grant Number 24249041.
The authors declare that they have no competing interests.
Conceived and Designed the study: TK, SM, YT, HO. Coordinated to obtain the data: YT. Data analysed: TK. Drafted the manuscript: TK, SM. All authors read and approved the final manuscript.
Number of household transmission events and mean intervals sorted by from–to transmissions.