INTRODUCTION
Each year, about 30 000 Americans die from influenza. During the 2009 H1N1 pandemic, children and young adults were about 10 times more likely to die than during seasonal influenza [Reference Morens, Taubenberger and Fauci1]. However, most fatalities were attributed to secondary bacterial infection predominately caused by Streptococcus pneumoniae (Sp) [Reference Morens, Taubenberger and Fauci1, Reference Weinberger2] and, to a lesser extent, Staphylococcus aureus (Sa), rather than influenza [Reference Mertz3]. Both species are common commensals of the nasopharynx and oropharynx [Reference Mertz3, Reference Nilsson and Ripa4]; Sp and influenza have well-characterized interactions, Sa and influenza less so. The duration and frequency of pneumococcal colonization varies by age and exposure to young children, the group at highest risk (30–80%) of carriage [Reference Peerbooms5], but all individuals are likely colonized with Sp at some point [Reference Bogaert, De Groot and Hermans6]. Pneumococcal colonization rates are lower in older children and adults (8–15%) [Reference Cardozo7, Reference Rawlings, Higgins and Han8], and children seem to carry Sp for longer periods than adults [Reference Hill9]. Sa carriage in the nares is also common (on average 32·4% of the US population), although neonates and youths aged between 7 and 19 years have the highest rates [Reference Mainous10]. Carriage of Sa in the nares and oropharynx are overall similar, and throat swabs may even be more sensitive than nasal swabs for detecting Sa carriage [Reference Mertz11].
Animal models [Reference Kukavica-Ibrulj12, Reference McCullers13], epidemiological studies [Reference Nilsson and Ripa4, Reference Wolter14] and previous work from our group [Reference Morens, Taubenberger and Fauci1] indicate that the transmission and pathogenicity of Sp is strongly affected by influenza infection. Influenza increases susceptibility and severity of pneumococcal invasive disease and facilitates pneumococcal colonization and growth [Reference Nilsson and Ripa4, Reference Kukavica-Ibrulj12–Reference Shrestha16]. Sp colonization may also modify risk of upper respiratory viral infection and viral infection may facilitate transmission of Sp [Reference Rawlings, Higgins and Han8, Reference Mertz11, Reference Diavatopoulos17]. Sa colonization is also increased by concurrent influenza infection [Reference van den Bergh18]. Sa was a major cause of mortality in each of the recent influenza pandemics – surpassing the mortality attributable to Sp in the 1957 and 1968 pandemics [Reference McCullers19]. One possible explanation for the role of Sp and Sa in post-viral pneumonia is that both bacteria are better able to colonize the lower respiratory tract than other commensals, and compete for the same ecological niches. Further, in vitro [Reference Regev-Yochay20] and epidemiological studies [Reference Regev-Yochay21–Reference Bogaert23] suggest a negative interaction between Sp and Sa carriage in the nasopharynx. However, this negative interaction seems to be specific to certain pneumococcal serotypes, and may not be sustained with viral infection [Reference van den Bergh18, Reference Regev-Yochay22, Reference Rodrigues24]. Additionally, although carriage in the two sites may be concordant, no study has specifically analysed co-colonization in the oropharynx, which may be more likely to be a source of lower-respiratory infections than the nasopharynx.
Despite being a common complication of influenza infection, the overall incidence of pneumonia is relatively low in the general population (~1/1000 cases of influenza). This has made it difficult to study the interactions of influenza, Sp, and Sa in human populations. Hence, much of the understanding of these complex co-colonization dynamics comes from animal models. However, since bacterial pneumonia requires aspiration of the infecting bacteria into the lungs, the same mechanisms that ultimately lead to pneumonia may also lead to increases in colonization rates and colonizing bacterial load in the oropharynx. High Sp colonization density is associated with invasive pneumococcal pneumonia [Reference Wolter14]. In a study conducted among college students, Sp colonization, but not Sa, was higher in those with viral acute respiratory illness (ARI) than in healthy contacts [Reference Davis25]. Similarly, a Portuguese study conducted in 585 daycare children aged <3 years found the prevalence of Sp, but not Sa, increased with the number of respiratory symptoms [Reference Rodrigues24].
Therefore, we chose to study the carriage rates and relative abundances of Sp and Sa in the oropharynx as an intermediate outcome towards the understanding of the increased risk of invasive pneumococcal or staphylococcal infections associated with ARI. We hypothesized that (i) the negative interaction between Sp and Sa may not be sustained during viral infection, and (ii) that in the presence of virus – particularly influenza, both Sp and Sa overgrow. Because rates of colonization by both bacteria are age-dependent, we took advantage of an ongoing prospective cohort study of ARI conducted in households including young children and their parents, to estimate the rate of carriage and relative abundance of Sp and Sa in these different age groups by detection of influenza or other respiratory viruses. We also conducted an exploratory analysis that assessed whether Sp colonization was associated with repeated ARI.
MATERIAL AND METHODS
Study population
We prospectively studied the subset of participants in The Household Influenza Vaccine Effectiveness (HIVE) study aged ⩾2 years that developed ARI. Between October 2013 and May 2014, ARI was actively monitored in 290 southeastern Michigan households with at least two children, as described previously [Reference Monto26]. Households were contacted weekly to identify members with ARI defined by the presence of an illness of <7 days duration with ⩾2 of the following symptoms: cough, fever or feverishness, nasal congestion, chills, headache, body aches, or sore throat. Specimens were collected only from those meeting the ARI case definition. Participants could be included on more than one visit if another ARI occurred. The HIVE study collected sociodemographic data, and medical record/state registry documentation of both the 2013–2014 influenza and pneumococcal vaccines [Prevnar7® (PCV7; Wyeth Pharmaceuticals, USA) and Prevnar13® (PCV13; Pfizer Inc., USA) and/or the polysaccharide pneumococcal vaccine, Pneumovax23® (PPV23; Merck, USA)]. Participants were considered fully immunized for influenza if they had been vaccinated during the current influenza season before the onset of ARI and for pneumococcus if they had received three PCV doses in the first year of life or one after age 12 months for PCV7 or PCV13, and one dose in the last 10 years for PPV23 [Reference Loughlin27]. Health system medical records were reviewed to document high-risk health conditions for complications of influenza, following the Advisory Committee on Immunization Practices (ACIP) [28]. The HIVE study, and this study protocol, were approved by the University of Michigan Medical School Institutional Review Board.
Sample collection and conservation
Participants with ARI had a throat and nasal swab collected which were combined together and conserved in 1 ml of universal transport media (UTM, Copan, USA) at −80 °C for viral testing. A second throat swab was collected and conserved separately in 1 ml UTM for bacterial testing.
Detection of Sp and Sa
Total bacterial DNA was extracted from the throat sample using a commercially available protocol (Qiagen Dneasy, USA). The presence of Sp and Sa was assessed by quantitative real-time PCR (qPCR). For Sp, the lytA gene was targeted using lytA-CDC primers, as described previously [Reference Carvalho29]. For Sa, the well-conserved and specific nuc gene was amplified [Reference Brakstad, Aasbakk and Maeland30]. Quantitative PCR reactions were carried out in a total volume of 20 µl, using a commercially available SYBR Green PCR master mix (Ssofast EvaGreen Supermix, Bio-Rad, USA), 200 nm primers and 5 µl DNA. The lytA-CDC gene fragment was obtained from PCR from a reference Sp strain (S. pneumoniae ATCC 49619). Assay conditions were optimized to reduce the amount of primer dimer in order to maximize the bacteria-specific signal in the qPCR assay. The limit of detection was 103 genomic copies/ml. All assays were conducted in triplicate and the mean values used for analysis. A no-template control and Sp-positive (ATCC 49619) and Sa-positive (MRSA 1268 strain) DNA controls were included in every run. Negative samples were defined as those with cycle threshold values >40. Serotypes for the Sp-positive samples were determined by multiplex PCR, as described previously [Reference Pai, Gertz and Beall31]. Isolates were typed for all serotypes found in PCV7 and PCV13 and additional common clinically relevant serotypes (11A, 15A, 15B, 16F, 20, 22F, 33F, 35B). Positive controls were used for all serotypes tested, from our laboratory collection and the Active Bacterial Core surveillance [ABCs/Emerging Infections Programs (EIP)] network.
Viral testing
Combined nasal and throat swab specimens were tested for influenza (in real-time) and 11 other common respiratory viruses (after conservation at −80 °C) by real-time reverse transcription–PCR (RT–PCR) on an ABI 7500 RT-PCR platform (Life Technologies, USA). RNA and DNA/RNA was extracted for influenza and additional respiratory virus testing using the Qiagen QIAamp Viral RNA Mini kit and MinElute Virus Vacuum kit (Qiagen, Germany), respectively. The Super-Script III Platinum One-Step Quantitative RT–PCR system and the AgPath-ID One-Step RT–PCR system (ThermoFisher Scientific, USA), plus primers and probes developed by the Centers for Disease Control and Prevention Influenza and Viral, Gastroenteritis, and Respiratory Viruses Divisions were used for detection of influenza A and B, respiratory syncytial viruses (RSVs), human metapneumoviruses (HMPVs), parainfluenza viruses 1–3, rhinoviruses, adenoviruses, and four coronaviruses (229E, OC43, HKU1, NL63), respectively, as described previously [Reference Ohmit32].
Laboratory testing was performed in the investigators' microbiology and respiratory virus laboratories at the University of Michigan School of Public Health.
Statistical analysis
Patients' characteristics were compared between adults and children, and Sp and Sa carriers and non-carriers using χ 2 test or Fisher's exact test for categorical variables and Student's t test or Wilcoxon rank-sum test for continuous variables. We used similar age groups as other studies in order to be able to compare our results with previously published works [Reference Regev-Yochay21, Reference Bogaert23]. In order to account for the longitudinal data (multiple observations per subject) and household correlation (multiple subjects per household) we used generalized estimating equation (GEE) models with an exchangeable correlation structure and robust standard errors to estimate the association between Sp detection and the following variables: age group, detection of influenza or other respiratory viruses, and Sa carriage. Since detection of Sp was a relatively common outcome, log-binomial regression models were fit to estimate the relative risk (RR). Daycare or school was not included in multivariate models because it was strongly associated with age group. We also fitted a logistic regression model to determine whether visiting the clinic more than twice with ARI during the season was independently associated with carriage of Sp or Sa on the first visit, age group or participants' daytime activity. Both models were tested for the whole population, and separately for children.
We evaluated co-colonization between respiratory virus detection and Sa carriage at any point during surveillance by describing departures from additivity and multiplicativity on the risk scale. Departures from additivity were evaluated by first calculating the RR of Sp in those with a virus detected compared to those with no virus detected, stratified by Sa carriage. These estimates were compared to the expected joint effect by estimating the relative excess risk due to interaction (RERI) and 95% confidence intervals (CIs), calculated using the delta method [Reference Hosmer and Lemeshow33]. A negative RERI value implies antagonism [Reference Knol34]. To evaluate departures from multiplicativity we included a product term in multivariable regression models. For all analyses a P value <0·05 was considered to be statistically significant. Multivariable regression models were performed using SAS software v. 9.3 (SAS Institute Inc., USA), RERI estimates and 95% CIs were estimated using the epiR package in R v. 3.1.3 (R Foundation, Austria).
RESULTS
This analysis included 338 participants (152 adults, 186 children) with ARI residing in 143 different households. Overall, 119/338 (35·2%, 95% CI 30·1–40·4) participants carried Sp on at least one visit, 52 participants carried Sa (13·4%, 95% CI 11·9–19·6) and 25 (7·4%, 95% CI 5·1–10·7) carried both, with young children representing the majority of both Sp and Sa carriers (Table 1). Sp carriage was the highest (70%, 95% CI 48–85) for the 2–3 years age group and Sa carriage was highest (22%, 95% CI 7–22·7) for the 7–12 years age group (Fig. 1). Characteristics of subjects by Sp and/or Sa carriage are described in Table 2. At least one of the 13 viruses tested was detected in 224 participants (66·3%, 95% CI 61–71·1), including 31 (9%, 95% CI 6·3–12·7) with influenza A/H1N1 and four (1·2%, 95% CI 0·3–3·1) with influenza B/Yamagata. Of the 35 participants testing positive for influenza, 14 (40%, 95% CI 23·8–57·8) carried Sp and six (17%, 95% CI 6·2–32·3) carried Sa, significantly more than those carrying Sp, 105/303 (34·6%, 95% CI 29·3–40·3), and Sa, 22/303 (7·2%, 95% CI 4·6–10·7), in the influenza-negative group (P = 0·03 and P = 0·04 for Sp and Sa, respectively) (Fig. 2a ). Sp, but not Sa, was carried significantly more often by participants with any virus detected than those with no virus detected (39·2% vs. 27·9%, P = 0·03 for Sp and 11·6% vs. 14%, P = 0·6 for Sa) (Fig. 2b ).
ARI, Acute respiratory illness.
* As defined by the Advisory Committee on Immunization Practices (ACIP) [28].
† Including home schooling for children.
‡ Pneumococcal vaccination was considered complete if the participant had received three PCV doses in the first year of life or at least one dose after age 12 months for PCV7 or PCV13, and at least one immunization in the last 10 years for PPV23 [Reference Loughlin27].
§ One child received both PCV7 and PPV23 vaccinations.
* As defined by the Advisory Committee on Immunization Practices (ACIP) [28].
† Pneumococcal vaccination was considered complete according to recommendations [Reference Loughlin27].
In children, both Sa carriage and influenza infection suggested an increase in risk of Sp carriage (RR 1·35, 95% CI 0·96–1·91; RR 1·47, 95% CI 1·01–2·14, respectively). However, when both Sa and influenza were present, Sp carriage tended to be less frequent than the expected additive effect of both Sa and influenza infection (observed RR 1·14, 95% CI 0·42–3·09 for an expected additive RR of 1·6) (Table 3). A similar departure from the expected joint effect was observed when all tested viruses were taken into account (for Sa alone: RR 1·76, 95% CI 0·81–3·80; for viruses alone: RR 1·31, 95% CI 0·99–1·75, observed RR for both Sa and virus 1·8, 95% CI 1·06–3·05, for an expected RR of 2·28). The negative RERI point estimate suggests that Sp carriage in children with ARI is less frequent than expected when Sa is also present, although this result is not statistically significant. A similar trend was not observed in adults, who were less likely than children to carry both Sp or Sa and to be infected by a virus. In the log-binomial adjusted regression model, only younger age was independently associated with Sp carriage (Table 4).
Sp+, Subjects carrying S. pneumoniae; RR, relative risk; CI, confidence interval; RERI, relative excess risk due to interaction.
* Any virus includes all 13 virus tested for (influenza A/B, adenovirus; coronaviruses 229E, HKU1, NL63, and OC43; human metapneumovirus, parainfluenza 1, 2, 3; rhinovirus; respiratory syncytial virus).
† The RERI estimates departures from additivity on the relative risk scale. A negative value implies antagonism, the null value is 0. RERI measures antagonism between virus detection (influenza or any virus) and Sa carriage in regards to Sp carriage and not antagonism between Sa and Sp.
* Influenza A (pH1N1 and H3N2) and influenza B Yamagata viruses were detected by RT–PCR.
† The following 11 viruses were detected by RT–PCR: adenovirus; coronaviruses 229E, HKU1, NL63, and OC43; human metapneumovirus; parainfluenza 1, 2, 3; rhinovirus; and respiratory syncytial virus.
‡ Adjusted models included age group, Sa carriage, influenza detection, other respiratory virus detection, and product terms for the interaction between influenza and Sa and other respiratory viruses and Sa.
All 338 participants were included in the model, with no missing data. The participants' daytime activities (being in daycare, at school or at work) were not independent from age, they were therefore removed from the final regression model.
Sa-positive children were more likely than Sa-negative children to be infected with any virus (29/36, 80% vs. 107/186, 57·5%, P < 0·001). However, after adjustment for age and Sp carriage, this effect was minimal [odds ratio (OR) 1·3, 95% CI 0·75–2·18, P = 0·8]. Sp bacterial load was also higher in subjects with virus detected than those without (3·65 log genomic copies/ml vs. 3.3 log genomic copies/ml, P = 0·05). By contrast, Sa bacterial loads (3·91 vs. 3·98 log genomic copies/ml, P = 0·5) did not differ between virus-positive and virus-negative participants.
For 110 participants (54 adults, 56 children), ⩾2 ARI episodes occurred during the influenza season. Of these, 45 (41%) carried Sp (14 adults, 31 children) and nine (8·2%) carried Sa (three adults, six children) at least once. The probability of returning for a second ARI episode was significantly higher in adults carrying Sp on their first visit (39% for Sp carriers vs. 20% for non-Sp carriers, P < 0·001), but not for children (27% and 33%, respectively) (Fig. 3). After adjusting using logistic regression for age, household and daytime activity for children (staying at home, being at school, or daycare), Sp carriers were twice as likely as non-carriers (OR 2·14, 95% CI 1·07–4·26, P = 0·04) to return for a second visit. By contrast, Sa carriage was not associated with number of ARI episodes.
Sp and Sa carriage in households
Participating households had a median of 4·5 subjects (range 4–8), with at least one adult and two children per household, but household size did not vary by carriage of Sp or Sa. Ninety (53·8%) households accounted for two-thirds (62·7%) of participants in which Sp was detected; children in these households were of younger age than households with no Sp (average age 6·7 vs. 9·7 years for Sp-negative households, P < 0·001). By contrast, the average age of children in the 45 (27%) households with at least one Sa carrier [141 (41·7%) participants] was not different from those where no Sa was detected (7·7 vs. 8·1 years, P = 0·4).
Sp serotypes
Serotypes were identified for 66 (45%) of the 148 Sp-positive samples (Supplementary Table S1). One adult (3·2%) and nine children (10·3%) carried more than one serotype concurrently. Only one adult (3%) and five children (3%) carried a serotype included in PCV7. Serotypes 35B (three adults, 11 children), 16 F (two adults, eight children) and 33 F (three adults, four children) were the most frequently detected. Thirteen subjects (four adults, nine children) carried the same serotype consistently during 2–4 visits, with an average follow-up time of 64 days (minimum 8 days between two visits, maximum 164). In three households, the same serotypes were carried by all participants: in one household, three children carried serotype 19A on the same date, in another household, three children carried serotype 35B and all tested positive for influenza B on the same date, and in a third household both a child and his parents carried 19A on the same date.
DISCUSSION
Our household study of 338 adults and children with ARI provides some of the first estimates of the prevalence of Sp and Sa carriage in the oropharynx in individuals with acute respiratory illness. We identified three major insights into the co-colonization dynamics of Sp, Sa and ARI. First, influenza and other common respiratory viruses were positively associated with Sp carriage and an association with Sa was observed only in subjects infected by influenza. Second, we did not observe a negative interaction between Sa and Sp colonization in the oropharynx. Third, those who carried Sp on their initial visit were more likely than non-carriers to return with subsequent ARI episodes, independently of age or household. Consistent with recent findings, we found a decrease in carriage of PCV vaccine serotypes in immunized children and non-immunized adults, but no decrease in overall carriage rates, indicating a switch in the serotype distribution since the implementation of PCV vaccination and efficient herd immunity [Reference Loughlin27, Reference van Hoek35].
Overall carriage rates of Sp in children were similar to those found in the literature, while adults' rates were higher, probably because they lived in households with young children [Reference Peerbooms5, Reference van den Bergh18]. However, Sa rates were significantly lower than those found in healthy populations [Reference Mainous10, Reference van den Bergh18, Reference Regev-Yochay21] where the prevalence based on culture-dependent methods is close to 37% [Reference Mertz11], with Sa carriage decreasing with age [Reference Graham, Lin and Larson36]. We are not the first to note a lower Sa prevalence with respiratory symptoms: in a Portuguese daycare study Sa prevalence was (13·9%), and Sa carriage decreased with severity of rhinitis symptoms. Similarly, in a previous work from our group, Sa recovery decreased with days since onset of ARI symptoms [Reference Davis25]. The understanding of this phenomenon is currently unclear, although many molecular pathways, such as increased mucus production, may explain the decrease in Sa detection in the oropharynx during ARI [Reference Shuter, Hatcher and Lowy37].
When the first PCV7 vaccination was approved and implemented, concern arose that decreasing Sp nasopharyngeal colonization might allow overgrowth of Sa and possibly be responsible for an increase in subsequent Sa infections [Reference Regev-Yochay21, Reference Bogaert23]. A 2004 Israeli study of 790 healthy children showed that Sp nasopharyngeal carriage, specifically of PCV7 vaccine-type strains, was negatively associated with Sa nasal carriage in children, even after controlling for age [Reference van den Bergh18]. As the authors argue, this apparent negative correlation may be due to specific bacterial antagonisms or could be a consequence of confounding effects [Reference Regev-Yochay22, Reference Reiss-Mandel and Regev-Yochay38]. Yet, although serotype replacement has largely occurred since PCV vaccination [Reference Loughlin27, Reference van Hoek35], there has been no corresponding increase in incidence of Sa infections in the post-vaccine era. Later works, including a recent randomized controlled trial, have confirmed that the implementation of PCV vaccinations had no impact in terms of Sa carriage despite driving Sp serotype replacement [Reference Lewnard39].
Previous epidemiological and experimental findings of an antagonism between Sa and Sp co-colonization in the face of viral infection focused on the nares and the nasopharynx [Reference Regev-Yochay20–Reference Bogaert23]. Although our sample size is small, and we did not have an ARI-negative control group, we did not observe an antagonism between Sa and Sp in the oropharynx. However, in children, who have the highest colonization rates of both Sa and Sp, the effect of viral infection on Sp colonization was lower, although not significantly, in the presence of Sa. In all age groups Sa carriage was lower than reported in the healthy US population, despite the use of a very sensitive detection method [Reference Mainous10, Reference Graham, Lin and Larson36].
In our sample, Sp colonization occurred more frequently in the oropharynx of participants infected by influenza and/or one of 11 other viruses than in those with no virus found, whereas Sa colonization was more frequent only in participants infected by influenza. Our inability to detect a viral infection may be due to the actual lack of viral infection, a high viral detection limit, or an infection with a virus not tested here. We screened for the 13 viruses most frequently associated with ARI in both adults and children and used a highly sensitive detection system [Reference Ohmit32, Reference Sakthivel40]; yet we cannot rule out that subjects may have been infected with another virus or with very low titres of one of the viruses screened for here, which constitutes a limitation to this work.
Our findings are consistent with a previous study from our group where Sp, but not Sa nor Haemophilus influenzae colonization was associated with detection of a respiratory virus in college students [Reference Davis25], and the Portuguese daycare study which showed an increase in Sp colonization, but not of Sa or H. influenzae, in children with respiratory symptoms [Reference Rodrigues24]. The association of Sp colonization and viral detection may be the consequence of Sp proliferation during viral infection, as found in animal models [Reference McCullers13, Reference Short41]. Importantly, higher Sp colonization densities are independently associated with invasive pneumococcal infection [Reference Wolter14]. Other studies have found a positive association between Sp and other respiratory viruses, such as HMPV and RSV [Reference Kukavica-Ibrulj12, Reference Rodrigues24, Reference Stark42, Reference Weinberger43]. This suggests specific, multifactorial, immunological consequences of viral infection that lead to increasing Sp growth following ARI. These may include aberrant innate immune responses to bacterial invaders in the immediate post viral state, uncontrolled pro- and anti-inflammatory cytokine production and diminished epithelial and mucosal defences [Reference Short41].
Viral infection may also increase susceptibility to Sp acquisition, as suggested by the mechanistic transmission model previously published by our group [Reference Shrestha16], and a 2014 study which found that influenza and parainfluenza ARIs appeared to facilitate pneumococcal acquisition in young children [Reference Grijalva44]. Although the mechanisms remain unclear, knowing that Sp colonization is more frequent in the oropharynx during viral disease including influenza is an important step towards the understanding of the pathophysiology of pneumococcal infection following viral infection.
The co-colonization dynamics between Sp and respiratory viruses may work in both directions. Indeed, we found that participants carrying Sp on their first visit were more likely to return with another ARI episode, independently of age and of household (with the limit of the small number of households involved). This finding is consistent with results from animal models that found that Sp infection may increase influenza titres and viral release from infected cells [Reference Short41], and in vitro data suggesting that exposure to Sp can increase susceptibility to HMPV infection [Reference Verkaik45]. Unfortunately, we did not have samples collected before the onset of ARI, which would have allowed us to characterize the impact of Sp colonization on the risk of ARI. We plan to specifically address this question in future studies.
In summary, our study complements the growing literature regarding the complex dynamics of co-colonization between Sp and Sa and respiratory viral illnesses in animal and in vitro models by adding results of co-colonization in the oropharynx from a household study in adults and children with ARI. Future studies should test for Sp carriage prior to the onset of ARI, so that the risk of ARI given Sp colonization can be estimated. Understanding these complex relationships is an essential step towards helping clinicians prevent pneumococcal and staphylococcal infections in the aftermath of acute respiratory viral infections.
SUPPLEMENTARY MATERIAL
For supplementary material accompanying this paper visit http://dx.doi.org/10.1017/S0950268816001473.
ACKNOWLEDGEMENTS
This work was made possible thanks to an Advancing Science through Pfizer-Investigator Research (ASPIRE) Award 2013 awarded to Betsy Foxman. The Household Influenza Vaccine Effectiveness (HIVE) study is supported by the Centers for Disease Control and Prevention (U01 IP000474) and the National Institute of Allergy and Infectious Diseases (R01 AI097150).
V.dL. would like to thank the Philippe Foundation and the Assistance-Publique Hôpitaux de Paris for financial support.
DECLARATION OF INTEREST
None.