Household transmission of respiratory viruses – assessment of viral, individual and household characteristics in a population study of healthy Australian adults

Background Household transmission of influenza-like illness (ILI) may vary with viral and demographic characteristics. We examined the effect of these factors in a population-based sample of adults with ILI. Methods We conducted a prospective cohort study in community-dwelling Australian adults nested within an influenza vaccine effectiveness trial. On presentation with ILI, participants were swabbed for a range of respiratory viruses and asked to return a questionnaire collecting details of household members with or without similar symptoms. We used logistic and Poisson regression to assess the key characteristics of household transmission. Results 258 participants from multi-occupancy households experienced 279 ILI episodes and returned a questionnaire. Of these, 183 were the primary case in the household allowing assessment of factors associated with transmission. Transmission was significantly associated in univariate analyses with female sex (27% vs. 13%, risk ratio (RR) = 2.13 (1.08, 4.21)) and the presence of a child in the house (33% vs. 17%, RR = 1.90 (1.11, 3.26)). The secondary household attack proportion (SHAP) was 0.14, higher if influenza was isolated (RR = 2.1 (1.0, 4.5)). Vaccinated participants who nonetheless became infected with influenza had a higher SHAP (Incidence RR = 5.24 (2.17, 12.6)). Conclusions The increased SHAP in households of vaccinated participants who nonetheless had confirmed influenza infection supports the hypothesis that in years of vaccine mismatch, not only is influenza vaccine less protective for the vaccine recipient, but that the population’s immunity is also lower.


Background
Improved characterisation of the determinants of household transmission of influenza-like illness (ILI) remains an important public health priority, particularly in light of the past decade's events in which we have witnessed the emergence of severe-acute-respiratory-syndrome (SARS) and the 2009 H1N1 influenza pandemic. The evidence base for pandemic influenza public health interventions such as home-quarantine, provision of antiviral agents for post-exposure prophylaxis, school-closure and vaccination builds upon an appropriate understanding of the patterns and timing of infection within the household unit [1][2][3][4][5][6].
Reflecting the need to improve our understanding of household transmission of ARI, the literature examining factors associated with household transmission of influenza [6,[18][19][20][21] has expanded significantly since the 2009 H1N1 influenza pandemic [22][23][24][25][26][27][28][29][30][31], including a systematic review and meta-analysis [32]. Donnelly et al. estimated the serial interval for all ILI (without laboratory confirmation) from case reports during the 2009 pandemic [25]. Only two studies of which we are aware explicitly consider the impact of virus type on infectiousness. Principi et al. found less onwards transmission to household members from influenza-negative than influenza-positive children presenting to a hospital emergency department [33]. Similarly, in a cohort study of ARI in young children, Lambert et. al. observed significant heterogeneity in the proportion of participants' households in which one or more illness events were observed (ranging from 13% for isolation of hMPV from the child to 61% for isolation of influenza) [34].
Here we report on the household transmission of a range of viruses in a cohort study of healthy, community-dwelling adults reporting symptoms of influenza-like illness (ILI). The study population was sourced from a large, industry sponsored placebocontrolled phase IV efficacy trial of a licensed seasonal trivalent influenza vaccine (Fluvax W , CSL Ltd), conducted prior to the 2009 pandemic between March and November 2008. In a previous article [35] we have described the viral aetiology of ILI in the cohort, and examined the influence of virus type, host and spatiotemporal factors on disease symptomatology.

Subject recruitment and selection
Full details of subject recruitment and selection for the primary phase IV vaccine efficacy trial have been described previously [35]. Briefly, across 23 study sites in Australia and New Zealand, 7544 healthy adults aged ≥18 to <65 years were recruited for a placebocontrolled trial of a licensed trivalent influenza vaccine (Fluvax W , CSL Ltd) in 2008 (Clinicaltrials.gov #NCT00562484). Study participants were randomized to receive either placebo or vaccine in a 1:2 ratio prior to the southern hemisphere 2008 influenza season. From an available pool of 5624 participants from the primary study at 12 study sites, we consider the 581 persons (adults) who experienced at least one ILI episodemeeting the case definition of at least one respiratory symptom (cough, sore throat, runny nose or nasal congestion) and at least one systemic symptom (fever (oral temperature ≥ 37.8°C), feverishness, chills or myalgia) [35], and who provided written informed consent for participation in the nested cohort-study which required contribution of a valid biological sample (COPAN TM dry flocked swab). Samples were tested for a range of respiratory viruses using a combination of multiplexed and uniplexed conventional and real-time polymerase chain reaction (PCR) assays [35]. Of those, 258 were members of multi-occupancy households, allowing investigation of transmission within the household. Using a non-specific and sensitive ILI definition, they reported an episode of ILI on 279 occasions. For each episode they returned a study questionnaire (Additional file 1) detailing respiratory symptoms (see [35] for details), health seeking behaviour (health care provider consultations, hospital admission, time off work), household characteristics (number of adults (≥18 years) and children (<18 years)) and temporally associated symptoms of ILI (if any) in other household members. From herein, we consider the illness episode as the primary unit of analysis. The

Outcome measures and analysis
Virology results were classified into 5 virus groups [35]: influenza (influenza A, influenza B), coronaviruses (OC43, 229E, NL63, HKU1), picornaviruses, other viruses (parainfluenza viruses (1, 2, 3), adenoviruses, human metapneumovirus (hMPV), bocaviruses, RSV and KI and WU polyomaviruses) or none, where none indicates that no 'tested-for' virus was detected in the participant's sample, as opposed to a missing sample or inconclusive result. Study participants' vaccination status (as determined by the primary phase IV trial intervention), physical location (i.e. study site), and socio-demographic characteristics were also recorded.
The relevant outcome measure for this sub-analysis was evidence of transmission within households based on experience of symptomatic illness in at least one other member of the study participant's household. Study participants were asked to complete the diary on the day following cessation of their own symptoms. They recorded the date of onset of symptomatic ILIs in household members from between 14 days prior to the study participant's illness through to the day of diary completion. Note that the primary case in the household may or may not be the study participant.
For household ILI events in which the participant was the primary case (183 of 279), transmission may or may not have occurred in the household and so it is statistically valid to develop univariate and multivariate explanatory models. For participants who reported recurrent ILI episodes during the study in which the same virus was isolated, we exclude all but the first episode. Events with co-introduction, defined as onset of symptoms in the participant and one or more household members on the same day, were also excluded. Following these exclusions, 177 episodes remained for analysis.
For household ILI events in which the participant was not the primary case (95 of 279), while we do know who the introducer was for these events (via the questionnaire data), other household ILI events initiated by that introducer that did not involve the participant are unobserved. That is, any events in which a child (or for that matter, any other adult member of the household) introduced an infection that did not infect the study participant are not captured by the study protocol. This observation necessitates the exclusion of all household ILI events in which the participant was not the primary case from the analyses.
In the one remaining household ILI event, the status of the participant (primary or not) was unknown, so the episode was excluded from the analysis.
For the 177 episodes in which the participant was the primary case, logistic regression models were used to explore associations between host, demographic or virus variables with any observation of within-household transmission (outcome variable = transmission in household for each recorded ILI episode in a study participant). The secondary household attack proportion (SHAP) was calculated as the proportion of potentially exposed household members (assumed susceptible) experiencing illness, averaged over all recorded episodes. We present descriptive statistics for the SHAP and its variation by virus, participant and demographic variables. Poisson regression models were used to assess the influence of virus, participant and demographic variables on the number of secondary cases within a given household, offset against the number of potentially exposed household members (outcome variable = number of secondary cases in household for each recorded ILI episode in a study participant).
Vaccination status of participants was not included in the primary logistic and Poisson statistical analyses due to its known mitigating effect on the likelihood of influenza acquisition [35]. Investigation of the influence of prior immunisation on influenza transmission in 'breakthrough cases' was explored in a secondary analysis by inclusion of a statistical interaction term between vaccination and influenza-identification status.
We make an empirical calculation of the mean time between the onset of symptoms in the primary case and the onset of symptoms in the household contacts (the serial interval), for all household ILI events, events in which the participant was the primary case, and events in which influenza was isolated from the participant's virological sample.
All statistical analyses were conducted in Stata/IC 11.1. Figure 1 reports characteristics of the 258 multioccupancy households in which transmission did and did not occur. 28 study participants reported two or more ILI episodes during the course of the study. For two participants, who both experienced two episodes, picornavirus was isolated on both occasions. We only retain the first episode for each participant. The distribution of household size is dramatically different based on the absence or presence of children within the household ( Figure 2). In households without children, the distribution is left-skewed (mean household size = 3.02, standard deviation (SD) = 1.18, skewness = 1.45), while in households with children there is minimal skew (mean household size = 4.18, SD = 1.11, skewness = 0.103).  Table 1 summarizes the descriptive statistics (and logistic model results) associated with presence or absence of transmission in the household for the 177 household ILI events in which the participant was the primary case.

Results
There is marked co-linearity between the variables 'presence of child in household' , 'age-category' and  'household size'. For example, respondents aged 35 -44 years had significantly greater odds of having a child in the household than those aged 18 -24 years (OR 53.2 (13.0, 217)), while no participant aged more than 55 years lived with a child. The relationship between the household size distribution and presence or absence of children is depicted in Figure 2. We retained 'presence of child in household' in the final multivariate model for transmission due to its strong predictive role, intuitive appeal, presumed causal role in our observed (univariate) association with age-category, and previous research indicating an association between transmission and children [18,22,23,29,36]. In the multivariate model, the observed increased risk of transmission with female sex remains (OR = 2.45 (1.01, 5.93), p = 0.047). Presence of children in the household is both the strongest and most statistically significant factor associated with transmission (OR = 2.63 (1.18, 5.88), p = 0.018).

Factors associated with the number of secondary household infections
Within 258 multi-occupancy households, 177 primaryparticipant introductions gave rise to 54 secondary cases among 391 potentially exposed individuals, a secondary household attack rate (SHAP) of 0.138. Of 102 exposed children, 22  In households in which the participant was female, 41 secondary infections were reported among 238 exposed household members (SHAP = 0.172), compared with 13 secondary cases among 153 contacts in households in which the participant was male (SHAP = 0.085), a riskratio of 2.03 (1.12, 3.66), p = 0.016 (2-sided Fisher's exact). In households with children, 27 secondary infections were reported among 162 exposed household members (SHAP = 0.167), compared with 27 secondary cases among 229 contacts in households without children (SHAP = 0.118), a risk-ratio of 1.41 (0.863, 2.32), p = 0.182.
A multivariate Poisson regression model ( Table 2) was used to consider the influence of virus group and demographic characteristics on the number of reported secondary cases within a given household, offset against the number of potentially exposed household members. In correspondence with the logistic regression model for transmission, we include presence of children in the In a secondary analysis, we considered the influence of prior vaccination on the reported number of secondary household cases among participants testing positive for influenza compared with all other participants. In a Poisson model for secondary attacks including a statistical interaction between influenza detection (true/false) and vaccination status (placebo/vaccine), the IRR for influenza positive cases in those receiving placebo was 1.69 (0.421, 6.80), p = 0.459. The factor increase (interaction term) for the IRR for vaccinated participants was 3.10 (0.608, 15.8), p = 0.174, yielding a net IRR for vaccinated influenza-positive participants relative to vaccinated influenza-negative participants of 5.24 (2.17, 12.6), p < 0.001).

Serial interval
Under the simplifying assumption that the introducer of infection into the household is responsible for all subsequent infections, we may calculate an empiric serial interval, the time from symptom onset in one individual until symptom onset in another. We first consider infections to be related if symptoms are reported within 14 days following onset in the primary case. Across all virus-type isolations, we calculate a mean serial interval of 6.0 days (SD = 3.6) for all household ILI events (where the study participant was the primary case or otherwise), and 5.1 days (SD = 3.2) for the events in which the participant was the primary case. For the five events in which the primary participant had virologically confirmed influenza and transmission occurred, the mean serial interval was 4.5 days (SD = 1.6). If we limited the maximum serial interval to seven days, the mean was reduced to 4.0 days (SD = 1.7) for all household ILI events and 3.9 days (SD = 1.9) for events in which the

Sensitivity analysis
In our main analysis, we made two assumptions that we now subject to a sensitivity analysis. Of the 183 events in which our participant was the primary case, 6 were classified as co-introductions as (at least) one other household member recorded symptoms beginning on the same day. As the latent period for respiratory infections may vary from individual to individual, here we exclude a further 4 episodes in which there was a 1 day interval from onset of symptoms in the study participant to onset of symptoms in another household member. The resulting multivariate models (equivalent to Tables 1  and 2) are materially unchanged, with the expected slight reduction in statistical power (data not shown). A second assumption made was that, for participants who reported multiple ILI episodes during the study period, we only excluded the latter ILI episode where the same respiratory pathogen was isolated on both occasions. However, if we conservatively exclude all ILI episodes except for the first (10 episodes excluded (by virus type: 6 "none", 1 "picornavirus", 2 "influenza", 1 "coronavirus")), again we find no material change in either the logistic or Poisson analyses (data not shown).

Discussion
This study, notable in its consideration of a broad range of respiratory pathogens in addition to influenza, demonstrates that household transmission of ILI is most strongly associated with host and demographic factors: female sex and the presence of children within the household (Tables 1 and 2).
The observation that female sex may be associated with increased transmission in the absence of children (RR = 2.33 (0.919, 5.90), p = 0.059) is novel, perhaps suggesting that females are fundamentally more infectious, and not simply more connected to children (in terms of both their susceptibility compared with males if a child introduces infection, and their infectiousness to children if they are the primary household case). Behavioural differences whilst ill may drive such an observation. Alternatively, mechanisms by which influenza pathogenesis is sex dependent have been investigated [37]; whether or not differences extend to infectiousness and susceptibility is not clear. Barbara et al. have recently identified that the reporting of respiratory symptoms may be linked with risk perception [38] and hence gender [39]. Clearly, we cannot exclude the possibility of gender difference in the reporting of within household transmission.
The association between transmission and the presence of children within the household is consistent with many other studies [18,22,23,29,36]. The logistic and Poisson model findings (Tables 1 and 2) are consistent with an increased susceptibility for children. This is further supported by the observed increased SHAP in children compared to adults (0.216 compared to 0.110, a risk-ratio of 1. 95 (1.19, 3.19). The SHAP in adults did not differ by whether or not their household contained children, suggesting that other 'indirect' effects of children are less likely. As our study design limited the analysis to household events with an adult introducer, we were unable to assess the hypothesis that children may be more infectious than adults.
Our Poisson regression analysis on the number of secondary cases given that the participant was the primary case ( Table 2) indicates that isolation of influenza in the introducer of infection to the household is associated with an increase in the number of secondary cases. We explored this finding more deeply using a statistical interaction model. While somewhat limited by sample size, we found that in placebo recipients identification of influenza was not significantly associated with an increase in the number of secondary cases (IRR = 1.69 (0.421, 6.80), p = 0.459), while in vaccine recipients the IRR (relative to identification of any other virus, including 'none') was 5.24 (2.17, 12.6), p < 0.001. Note that our previous analysis confirms that vaccination is associated with a reduced probability of influenza virus identification [35]. Additionally, 'breakthrough' influenza cases have similar symptoms compared to unvaccinated individuals [35]. We therefore suggest that our finding of increased transmission may be explained by infection with an influenza virus mismatched to the vaccine-strain (known to be in circulation during the year of study [40]), which furthermore may be relatively antigenically novel and to which household members may be expected to have heightened susceptibility. With no virological samples available from other household members and the small number of vaccinated participants who were infected with influenza we are unable to explore this hypothesis further.
Across all virus types isolated and all household ILI events, and assuming that all secondary cases within the household are directly infected by the introducer, we calculate a serial interval of 6.0 days. Restricting to events in which the participant was the primary case and in which influenza was isolated, we calculate a serial interval of 4.5 days. This simple approach, as taken by others [22,29], cannot account for two important factors: community importation and infection of household members by other non-introducing members (i.e. tertiary cases). While others have partially accounted for these effects [25,29,41], a mechanistically-motivated statistical model is required to fully account for such possibilities, for example as introduced by Cauchemez et al. [23] who determined a serial interval for influenza of 2.6 days (SD = 1.3) compared to 2.9 days if calculated directly from empirical observations. With just 5 events, application of these more advanced model-based techniques is not justified for our data.
Of primary interest for this sub-analysis focussed on transmission is the complication introduced by the monitoring and assessment of ILI in an individual rather than a household. Ideally, a protocol such as that suggested by Klick et al. would have been employed [42]. The lack of virological assessment of household secondary cases and the broad nature of the question used to establish the secondary case count in each house also contributes to uncertainty with regards to our assignment of temporally associated ILI to within-household transmission. Both of these limitations were an unavoidable consequence of the nesting of the data-collection protocol within a randomized placebo-controlled trial. Furthermore, due to the requested timing for completion of the questionnaire, we cannot exclude the possibility that late onset of secondary (or tertiary etc.) cases may have been missed, particularly if a participant's experience of symptoms was of short duration. Similarly, because the study protocol and analyses effectively assume that individuals are infectious until the end of their symptoms, any systematic differences (by virus type) in this relationship may influence the results. However the prompt to return the diary upon symptom cessation was in an effort to ensure timely reporting of questionnaire information to minimise recall bias. Conversely, our Poisson model implicitly assumes independence among household members, attributing all household infections to the primary case. More advanced model based methods that account for tertiary (and subsequent) cases and community introduction would be warranted with more complete data sources.
As with all protocols based purely on symptomatic presentation (as opposed to active surveillance for nonclinical signs of infection such as virological or immunological measures [24,42]), we are unable to account for potential sub-clinical infection routes, with potential impact for our assessment of whether or not transmission did occur, the primary case status of our participants and determination of the size of the susceptible pool within a given household. Conversely, taking a nonsimulation approach to analysis, we are unable to discount our estimate for the SHAP due to the effects of community introduction into the household, or account for community introduction and tertiary cases in our estimate for the serial interval [23].
Our study sample had an over-representation of females (166 of 258 (64.3%) individuals for the 279 captured episodes; 105 of 167 (62.9%) individuals for the 177 primary-participant introductions). Furthermore, it should be noted that the study population were originally volunteers in a randomized controlled trial and as such more likely to represent a group who were more concerned with their health than the general population. Eligibility was restricted to healthy adults without recognized risk factors for severe influenza infection.

Conclusions
In the context of a literature focussed on the transmission characteristics of laboratory confirmed influenza, our study is the only one that we know of to systematically explore the relationship between transmission and virus aetiology. The analyses suggest that influenza is more transmissible than other causative agents of ILI, at least when introduced to the household by an adult. Host and demographic factors are also of importance. Further studies combining active surveillance of all household members with specimen collection and testing for a range of respiratory pathogens are warranted to elucidate these relationships.

Additional file
Additional file 1: Illness Visit Questionnaire.
Competing interests PCR has previously served on a scientific advisory board regarding influenza vaccines for CSL Ltd and has received a grant for an investigator initiated epidemiological study of otitis media from GlaxoSmithKline Australia. He has also received travel support for himself and staff employed by the Vaccine Trial Group to attend and present data at scientific meetings from Baxter, GlaxoSmithKline, Sanofi and Pfizer. MDN has received travel grants from Wyeth Australia to present independent research at international meetings, and currently and previously has been the principal investigator for clinical trials sponsored by Abbott, Baxter, CSL, GSK, MedImmune, Merck, Novartis, Sanofi-Pasteur, Wyeth, and Pfizer. ML is an employee of CSL Limited and has an equity interest in the company.
Authors' contributions JMC, PH and JMV conducted the statistical analyses, provided the primary interpretation of the results and wrote the manuscript. PR was principal investigator on the vaccine efficacy trial within which the sub-study was conducted. TN, JMV, TS, MN, SL and PR conceived the sub-study and secured funding for its conduct, in partnership with CSL Limited represented by ML and MG. JMV coordinated conduct of the study at multiple sites and oversaw collation of the questionnaire data. TS, MN and SL oversaw conduct of and reporting of the virological testing at the Queensland Paediatric Infectious Diseases Laboratory. ML was medical monitor for the main vaccine study and a partner investigator on the sub-study, as was MG. All authors contributed to critical revision of the manuscript and have seen and approved the final version of the manuscript.