 Research article
 Open Access
 Published:
Is the impact of childhood influenza vaccination less than expected: a transmission modelling study
BMC Infectious Diseases volume 17, Article number: 258 (2017)
Abstract
Background
To reduce the burden of severe influenza, most industrialized countries target specific riskgroups with influenza vaccines, e.g. the elderly or individuals with comorbidities. Since children are the main spreaders, some countries have recently implemented childhood vaccination programs to reduce overall virus transmission and thereby influenza disease in the whole population. The introduction of childhood vaccination programs was often supported by modelling studies that predicted substantial incidence reductions. We developed a mathematical transmission model to examine the potential impact of childhood influenza vaccination in Germany, while also challenging established modelling assumptions.
Methods
We developed an agestratified SEIRtype transmission model to reproduce the epidemic influenza seasons between 2003/04 and 2013/14. The model was built upon German population counts, contact patterns, and vaccination history and was fitted to seasonal data on influenzaattributable medically attended acute respiratory infections (IMAARI) and strain distribution using Bayesian methods. As novelties we (i) implemented a stratified model structure enabling seasonal variability and (ii) deviated from the commonly assumed massactionprinciple by employing a phenomenological transmission rate.
Results
According to the model, by vaccinating primarily the elderly over ten seasons 4 million (95% prediction interval: 3.84 – 4.19) IMAARI were prevented which corresponds to an 8.6% (8.3% – 8.9%) reduction compared to a novaccination scenario and a numberneededtovaccinate (NNV) to prevent one IMAARI of 37.1 (35.5 – 38.7). Additional vaccination of 210 yearold children at 40% coverage would have led to an overall IMAARI reduction of 17.8% (17.1 – 18.7%) mostly due to indirect effects with a NNV of 20.7 (19.6 – 21.6). When employing the traditional massactionprinciple, the model predicted a more than 3fold higher IMAARI reduction (55.6%) due to childhood vaccination.
Conclusion
In Germany, the introduction of routine childhood influenza vaccination could considerably reduce IMAARI among all agegroups and improve the NNV. However, the predicted impact is much lower compared to previous studies, which is primarily caused by our phenomenological approach to modelling influenza virus transmission.
Background
Around the globe, annual influenza epidemics result in considerable morbidity and mortality in all parts of the population. Children have the highest attack rate [1, 2], but the elderly and people with comorbidities are at particular risk for severe influenza disease and influenzaassociated complications [3]. Therefore, Germany and many other industrialized countries recommend annual vaccination of these atrisk groups focusing on the direct benefits of the vaccine [4].
However, since the effectiveness of influenza vaccines is limited in these atrisk groups [5, 6] and because children are considered the primary spreaders of influenza, childhood vaccination programs might reduce influenza disease burden not only in the vaccinated but also in the nonvaccinated population through indirect effects [7]. In the UK, for example, routine vaccination of all children aged two to 17 years is recommended in addition to the vaccination of riskgroups. In the US, annual vaccination is recommended for all people above 6 months of age [8]. According to impact results, the UK childhood vaccination program caused considerable indirect effects also among adults [9]. Field studies from the USA suggested that indirect effects are very pronounced in closed communities [10]. However, the indirect impact of school vaccination programs on US county level turned out to be low [11, 12].
The availability on new live attenuated influenza vaccines (LAIV) also supported the introduction of childhood vaccination, as LAIV promised higher vaccine efficacy (VE) among children and easier administration compared to inactivated vaccines [13]. Recent results though question the beneficial effectiveness of LAIV [14], which led to withdrawal of LAIV recommendation in the US [15].
To predict the potential impact of including a new vaccine or implementing a new strategy in the national immunization program, many national decision makers utilize mathematical transmission modelling prior to introduction. For influenza, every transmission model published as of today predicted considerable populationwide case reductions due to childhood immunization programs, which contrasts the heterogeneous impact of vaccination observed in the field [16–19].
To explore explanations for the partial discrepancy between modelled and observed childhood vaccination impact, we propose a new modelling approach of stratification by season and subtype to handle variability of observed influenza epidemics. Moreover, we challenge established assumptions regarding the modelling of contact pattern and force of infection. Using Bayesian methods, our employed transmission model is calibrated through estimated numbers of medicallyattended influenza cases in Germany and virological data on the yearly distribution of influenza subtypes [20].
The objective of our work was to develop a transmission model that is capable of reproducing the influenza epidemics in Germany observed over the past 10 years. Using the model we estimated the number of medicallyattended influenza cases that could have been prevented by a childhood vaccination strategy in Germany, which provides necessary evidence to support immunization decision making.
Methods
We developed an agestratified SEIRtype transmission model to reproduce the epidemic influenza seasons between 2003/04 and 2013/14 excluding the pandemic season 2009/10. The model was based on German population counts and vaccination history and fitted to seasonal data on influenza excess consultations and strain distribution as described below. Here, we give an overview on the most relevant model aspects whereas full details are provided in the supplemental material (see Additional file 1).
Epidemiological data
Influenza surveillance and virological data
The German working group on influenza (Arbeitsgemeinschaft Influenza, AGI) operates a syndromic surveillance system using a GP sentinel network [20]. On average 600 participating practitioners provide weekly reports on the frequency of medicallyattended acute respiratory infections (MAARI) for five age groups: 04, 514, 1534, 3559, and 60+ years of age.
The season and agegroup specific number of influenzaattributable MAARI (IMAARI) are estimated through a time series approach based on the reported MAARI incidence [20].
The National Influenza Reference Center provided data of the virological surveillance of the AGI including seasonal distribution over the influenza types A(H1N1), A(H3N2) and B. Here, A(H1N1) refers either to the pandemic variant A(H1N1)pdm09 or to the prepandemic variant A(H1N1)prepan, depending on season. Since the season 2008/09 Bpositive samples are further divided into BYamagata and BVictoria yielding overall four different subtypes. Although inaccurate, from here on we use subtype to refer to the two A subtypes and the two Blineages. For seasons prior to 2008/09 with missing lineage information we identified the circulating Blineages from the annual AGI reports [21]. For the epidemic season 2005/06 that exhibited considerable cocirculation with 90% BVictoria and 10% BYamagata, we randomly assigned each Bpositive sample to one of the two lineages according to the Blineage distribution given in the AGI report [22].
Vaccination data
In Germany during the modelled time horizon, influenza vaccination was primarily recommended for atrisk groups, i.e. people aged 60+ years and individuals with comorbidities [23].
To reproduce the past vaccination history, we obtained vaccine coverage rates for each of the included seasons from health insurance claims data as described by Rieck et al. [24]. The seasonal coverage was highest among elderly over 60 years with on average 44.2 % between 2003/2004 and 2013/2014. The mean coverage among children (019 years) and adults under 60 years (2059 years) was 5.4% and 9.7%, respectively.
For the majority of modelled seasons only trivalent inactivated vaccines (TIV) from various manufacturers were administered in Germany, although in recent years also quadrivalent vaccines (QIV) and for children LAIV became available. Since market share of QIV and LAIV was small in Germany (data obtained from Insight Health GmbH & Co. KG) and information on the type/brand of administered influenza vaccine is not included in the insurance claims data, we assumed one generic trivalent vaccine to be applied in the model.
Because influenza viruses change dynamically over time (antigenic drift) and the vaccine is reformulated every season according to recommendations from the World Health Organization (WHO), the vaccine effectiveness (VE) may differ from season to season, by subtype, and by age. Thus, the VE in our model also distinguished between subtypes, seasons, and three age groups: 014, 1559, and 60+ years. For VE estimates of past seasons we relied on data from the European IMOVE network (Influenza  Monitoring Vaccine Effectiveness), Cochrane reviews, and AGI reports [25–32]. For each season since 2008/2009, we specified the model VE based on IMOVE that assessed country, subtype, and agespecific VE estimates [27–32]. We prioritized IMOVE estimates that were adjusted for study site, chronic condition, and other influential factors. Where IMOVE results were not available, we used estimates from Cochrane reviews that measured VE in both adults and children for wellmatched and poorlymatched seasons [25, 26]. Information on seasonal matching was taken from the annual AGI reports.
Although TIV contains only one Blineage component, some clinical studies detected a crossprotection for the Blineage that was not included in the vaccine, which was measured at 60% compared to the VE against the included Blineage [33, 34]. Thus, for the Blineage not included in TIV according to WHO recommendation we assumed this 60% crossprotection as already done in previous models [35]. Combining the observed seasonal Blineage mix with the overall Blineage VE estimates then yields VE estimates specific to each Blineage.
Transmission model
Model structure
The season and subtype specific transmission dynamics are captured through an agestructured SEIR model introduced by Vynnycky et al. and revisited elsewhere [16–19, 35–37]. The model divides each agegroup into susceptibles (S), latently infected but not yet infectious individuals (E), infectious (I), and recovered individuals (R). These are further split into vaccinated and unvaccinated individuals which yields in total eight agestratified model compartments.
During an epidemic season and for each subtype, susceptible people might acquire infection (S → E), become infectious (E → I), and upon recovery (I → R) become immune for the remainder of the season. These dynamics are formulated as a system of ordinary differential equations.
In each season only a fraction of the population starts as susceptible since initial immunity could be inherited from prior infections or established through vaccination. Here, we assume an allornothing vaccination effect [38]. Waning of vaccine protection during the season is not taken into account, but the protection is assumed to vanish completely at the end of each season.
We pursue a stratified modelling approach such that each subtype and season is embedded within a global parametric framework [39]. Hence, although the transmission dynamics of each subtype within each season are modelled separately, these single models share not only the same structure but also certain parameter values that are likely unaffected by season and subtype, e.g. contact patterns and consultation rates. Conversely, those transmission aspects that are presumably specific to either season or subtype, and thus are causing the observed seasonal variability, are allowed to differ accordingly. See Table 1 for a list of all model parameters and their corresponding stratification.
Susceptible Population
The susceptible fraction of the model population may differ with respect to each subtype and season, since varying susceptibility is assumed to be the primary cause for variability within seasonal influenza epidemics. Within our model, the susceptible fraction for each subtypeseason is controlled by two aspects: (i) an agespecific susceptibility σ _{ a } accounting for increased subtypespecific immunity within the two agegroups 1559 and 60+ years, and (ii) a seasonspecific factor φ that additionally shrinks the subtypespecific susceptible population fraction due to immunity gained within recently past seasons.
Disease transmission
For each agegroup i the force of infection λ _{ i } results from the age dependent contact behavior, the prevalence of infection, its transmissibility, and a constant risk λ _{ o } to import influenza infection from outside of the German population.
In contrast to implementing the massactionprinciple that assumes spatial homogeneity and is commonly used for pathogen transmission modelling, we apply a phenomenological transmission rate that enables a dampened growth of the force of infection for an increasing prevalence subject to the power parameter ρ ≤ 1 [40, 41]. This transmission rate is motivated by effects which lead to a declining effective reproduction rate, such as reactive behavior changes during an epidemic or a potential spatial clustering of the infection resulting from the nonhomogenous spread within a population [42]. In our model, the infectious pressure λ _{ i } applicable to susceptibles in age group i therefore increases concavely, i.e. at a decreasing marginal rate, with the infections prevalence in each age group j (given by the number of infectious people I _{ j }(t) divided by the population size N(t)_{ j }, j = 1 , … , n _{ a }). This concaveness is more pronounced for small values of ρ [43].
The contact rates \( {\beta}_{i, j}^{(eff)} \) denote the average number of effective contacts from individuals in age group i with individuals in age group j that are – according to the social contact hypothesis – proportional to the number of social contacts between these age group [44]. We estimated social contact frequencies based on the German part of the POLYMOD survey using a spline regression approach [44, 45]. We restricted the data to contacts of physical nature or of at least 15 minutes duration as these were found to be a good proxy for contacts likely leading to transmission of seasonal influenza [36, 46].
However, the POLYMOD study surveyed the behavior of primarily healthy people. A contact survey conducted throughout the pandemic season in 2009/2010 measured the social activity of symptomatically ill people and once again after their recovery [47]. Based on this data it was shown that ill people have a higher number of household contacts whereas school and work contacts are less frequent [48]. We utilized these results to estimate a contact matrix β ^{(sick)} subject to symptomatic illness. The effective contact matrix was then modelled to be a weighted combination of the two contact matrices for healthy and sick people.
The weighting parameter m > 1 secures that the contact matrix subject to ill people has a bigger impact due to the high probability for developing symptoms, but also because symptomatic people have a higher infectiousness [48].
The parameter R _{ e } controls the transmissibility of the infection whereas z(t) reflects the intraseasonal variation of the transmissibility as one oscillation over the season.
The parameter t _{ z } and t _{ s } control the timely shift of peak transmission with respect to subtype and season, respectively, and thus may vary accordingly. The magnitude of the seasonal oscillation is controlled by δ.
An overview on each parameters interpretation and potential stratification by season or subtype is given in Table 1.
Parameter estimation
Parameters that are well supported through data such as demographics and vaccination history are kept fixed within the model. Model parameters based on weak evidence were estimated within a Bayesian framework using the available disease burden data (see Table 1).
Prior distributions and plausible ranges for all estimated parameters were defined based on literature if available. Additionally, the likelihood function specified below measures the plausibility of the two data sets, IMAARI data and virological data, subject to the model.
The number of IMAARI \( {D}_{t, a}^{(s)} \) per week t and agegroup a determines for each season s the magnitude of the influenza wave as predicted by the model that is the aggregation of the four subtype specific waves. Thus, let \( {X}_{t, a}^{\left( z, s\right)}\left(\boldsymbol{\theta} \right) \) denote the model predicted number of influenza cases due to subtype z and subject to the parameter vector θ. The predicted number of IMAARI is then given by
Where 0.67 gives the probability for developing symptoms [49] and \( {p}_a^{(c)} \) refers to the age specific probability for seeking medical treatment, i.e. the consultation rate due to influenza. The probability for developing symptoms was chosen to be constant in order to secure identifiability of the model. The number of IMAARI estimated by AGI is then assumed to be negativebinomially distributed with expectation \( {Y}_{t, a}^{(s)}\left(\boldsymbol{\theta} \right) \) and dispersion d _{ t }, i.e.
The subtype distribution is governed through the virological data containing the weekly number of positive tests \( {P}_{t, a}^{\left( z, s\right)} \) for each subtype z. We constructed the likelihood for observing the subtype distribution \( {\boldsymbol{P}}_{t, a}^{(s)}=\left({P}_{t, a}^{\left( AH1 N1, s\right)},{P}_{t, a}^{\left( AH3 N2, s\right)},{P}_{t, a}^{\left( B Yam,\kern0.5em s\right)},{P}_{t, a}^{\left( B Vic, s\right)}\right) \) by assuming that each influenza case has the same probability of leading to a positive virological test and identification that leads to a Dirichletmultinomial likelihood.
The posterior distribution is then derived as the product of the prior and the two likelihood functions L(D θ) and L(P θ) assuming conditional independence. A sample from the posterior distribution is obtained by applying adaptive Markov chain Monte Carlo sampling as done previously [16, 50].
Investigated model scenarios
We developed alternative model scenarios to investigate the sensitivity of the respective results caused by different assumptions. Alternative assumptions included assuming no Blineage crossprotection (S1), switching from phenomenological transmission to massactiontransmission (S2), applying the contact matrix from only healthy individuals (S3), and disabling any indirect effect (S4). The fit of each model was measured through its respective marginal likelihood of the data [51].
Vaccination impact analysis
To assess the impact of the past vaccination program (historic vaccination) on the number of prevented IMAARI we simulated a scenario without vaccination. The potential impact of routine childhood vaccination was estimated by simulating a scenario with increased vaccination coverage of 40% among all children from 2 to 10 years additional to the historic vaccination. Furthermore, we varied the target group to include children up to 17 years or only up to 6 years and also examined alternative childhood coverage rates of 20% or 60%. The potential effect of QIV was examined by assuming 100% Blineage cross protection. Although the beneficial effect of LAIV is uncertain, in order to investigate the potential impact of an overall more effective influenza vaccine we assumed LAIV having a 50% higher VE among children age 2 to 6 years compared to TIV. We also investigated a scenario of increased coverage among elderly over 60 years. The impact of each scenario was measured as the IMAARI reduction in comparison to the scenario without vaccination.
Results
Model fitting
The ten reproduced influenza waves are shown in Fig. 1. The IMAARI curves including their magnitude, peak time, and age distribution are reflected by the modelpredicted waves. The distribution and timely occurrence of the single subtypes are well reproduced, with the only exception being the 2003/04 season for which the model predicts the circulation of also the BYamagata lineage. Corresponding parameter estimates are displayed in Table 1.
The estimated susceptible population fractions given in Fig. 2 vary considerably by season and subtype. For all subtypes the susceptible fraction among children was estimated to be much larger compared to older age groups. Moreover, the overall susceptible population fraction with respect to one single subtype, especially AH1N1, can vary considerably for subsequent seasons. An overall smaller population fraction was susceptible to the Blineages as compared to influenza A. Comparing prepandemic and pandemic A(H1N1), a smaller fraction of people over 60 years was susceptible to the pandemic strain A(H1N1)pdm09, whereas in the age group 15 to 60 years a larger fraction is susceptible to A(H1N1)pdm09.
Childhood vaccination impact
Across the ten modelled seasons the historic vaccination program targeting only atrisk groups has prevented 8.6% of all IMAARI that would have occurred in a scenario without vaccination as predicted by the model (Table 2). An additional vaccination uptake to 40% among 2 to 10 year old children could have reduced the influenza disease burden by overall 17.8% compared to a scenario without vaccination, which corresponds to an average 828,000 prevented IMAARI annually. Within the childhood vaccination scenario, the relative reduction per season ranged from 10.6% to 27.2% (Fig. 3). It was most pronounced among elderly with a mean IMAARI reduction of 31.0%, followed by children with 25.1% and 18.5% prevented IMAARI in the age groups <5 years and 514 years, respectively. IMAARI incidence among adults was predicted to be reduced by 12.1% and 14.5% for the age groups 15 to 34 years and 35 to 60 years, respectively. Extending the hypothetical childhood vaccination to children aged 2 to 17 years increased the overall IMAARI reduction to 22.9% whereas restricting childhood vaccination to children aged 2 to 6 years yields a reduction of 14.3% compared to a no vaccination scenario (Table 2). Administration of live attenuated vaccines for children or utilization of quadrivalent vaccines within the base childhood vaccination scenario (210 years, 40% coverage) would have led to an additional reduction of 3.8% or 2.3%, respectively. Increasing the coverage among elderly over 60 years to 100% would reduce IMAARI by 12.1% compared to no vaccination.
Vaccination scenario effectiveness
The additional implementation of childhood vaccination yields large indirect effects within the whole population, e.g. reducing IMAARI among elderly by an additional 11% compared to the historic vaccination scenario (Fig. 3). As a consequence childhood vaccination scenarios always reduced the number needed to vaccinate to prevent one IMAARI (Table 2). Thus, the historic vaccination scenario is least efficient at preventing IMAARI as it requires the most vaccine doses to prevent one IMAARI (Fig. 4). A vaccine doses allocation analyses yields that — under a restriction of only one million available vaccinations per year — the largest overall IMAARI reduction can be achieved by targeting the age group 24 years (5% reduction) whereas targeting any one age group over 60 years reduced IMAARI by less than 0.4% (Fig. 5).
Sensitivity analyses
Investigating different modelling scenarios, the results were found to be comparable for models assuming a different contact pattern or no Blineage cross protection (Table 3). The strongest influence on the model had the assumption of massactiontransmission, which resulted in a 55.6% relative reduction of IMAARI.
The highest marginal likelihood of the data was given by the base model and the model assuming a contact pattern according to POLYMOD. The massactiontransmission model yielded a considerably lower marginal likelihood (Table 3 which corresponds to a visibly worse model fit as displayed in Additional file 1: Figure S6.
Discussion
We developed a dynamic SEIRtype transmission model to reconstruct the seasonal influenza epidemics between the years 2003/04 and 2013/14. Each epidemic season was represented through its estimated number of IMAARI and its subtype distribution. When examining hypothetical childhood vaccination scenarios, the fitted model predicted considerable reduction in IMAARI across all agegroups, which go beyond the expected direct effects. The implementation of a phenomenological transmission rate had a major impact on the predicted effects of augmented vaccination scenarios, since massactiontransmission models suggested even more drastic incidence reductions.
The observed influenza waves are well reflected by the fitted model; in particular the variability in the magnitude of the epidemic seasons and the circulating subtypes (Fig. 1). The modelpredicted epidemics do not capture each nuance of the weekly IMAARI numbers. However, it should be noted that the weekly IMAARI are already estimated numbers such that the true IMAARI curve might be less fluctuate than suggested by the data [20]. This applies in particular for the weeks prior to beginning and after the end of each epidemic season. In those weeks the data yields zero IMAARI whereas the true number is certainly positive. Nevertheless, the advantage of relying on the estimated IMAARI is that they represent only influenzaattributable medical consultations. Thus – unlike e.g. ILI incidence data – the IMAARI are better capable to capture the true influenza epidemics, which directly improves the validity of the fitted model.
Applying a partial stratification of the model parameters offers the advantage that all observed variability in influenza epidemics is explained only by transmission aspects for which variation by season and subtype is epidemiologically plausible. On the contrary, transmission factors that likely stay consistent over the years – such as contact patterns or the consultation rate, i.e. the proportion of sick people seeking medical advice – are enforced to do so. Thus, the stratification approach provides robust estimates for such constant quantities, e.g. the consultation rate of children under 5 years was estimated to be 45.3% over all seasons and subtypes (Table 1). In comparison, by modelling each influenza season separately Baguelin et al. obtained heavily diverse ascertainment probabilities, i.e. probabilities for an influenza case leading to notification. For children under 15 years these ranged between 0.1% and 2% depending on season and subtype [16]. This diversity in the estimates is difficult to explain. Moreover, if the ascertainment probabilities were over or underestimated it also implies that other parameters might be falsely estimated since the unobserved influenza epidemics may be actually higher or lower, respectively. Secondly, our stratified model framework provides distinct estimates for those parameters where heterogeneity is expected, such as the distinct susceptibility profiles with respect to different subtypes and seasons (Fig. 2). These are in line with available data on population immunity against A(H1N1)pdm09 following the pandemic season 2009/10 in the UK, indicating a very high immunity among elderly and low immunity among children [52]. However, in Germany postpandemic antibody titres against A(H1N1)pdm09 were comparably low among elderly, although in the subsequent season the lowest infection rates were also detected among elderly [53]. This inconsistency might be caused from preexisting immunity not measurable through crossreactive antibodies [53], such that our estimated susceptibilities provide a rough assessment of the actual preexisting immunity. On the other side, model stratification by season in contrast to utilizing a continuous model as done by Vynnycky et al. and others [17, 18, 37] is necessary to obtain additional insight into the potential range of vaccination impact that might vary by season e.g. due to different VEs. Thus, the model predictions suggest that a strong influenza season is associated with a small relative vaccination impact (compare Figs. 1 and 3), presumably because antigenic virus drift may lead to high population susceptibility and vaccine mismatch at the same time. This association was also observed by more realistic networkbased models [54].
As a novel approach our model also incorporated contact data from people with ILI symptoms into the transmission process whereas most other models utilize POLYMOD data that have investigated contact patterns of healthy people only [44, 47]. Although the number of school contacts is reduced in the adjusted contact matrix, children still remain the group with the highest contact frequency and even have a higher number of household contacts such that transmission between children and adults is even more pronounced [48]. Childhood vaccination impact was thus even slightly larger compared to a model utilizing solely the POLYMOD data (S3) (Table 3). Also note that the POLYMOD based model was more likely from a statistical standpoint (Table 3). The worse fit of the model incorporating behavioural changes due to illness might originate from the respective contact data being collected in the UK during the pandemic season and thus these data are not necessarily applicable for seasonal influenza transmission in Germany. However, for the eventual selection of our base analysis we preferred a more realistic model accounting for adjusted contact behaviour over a slightly improved model fit.
Impact of Childhood vaccination
Investigating a childhood vaccination scenario targeting 210 year old children at 40% coverage in addition to the observed vaccination uptake in the past, our model predicts an overall reduction in IMAARI of 17.8% that is mostly caused by reductions among the vaccination target group of children and elderly themselves. Compared to the expected direct effects of 11.0% decrease in IMAARI, routine childhood vaccination is predicted to causes considerable positive indirect effects. However, throughout all other influenza transmission models the predicted impact of childhood vaccination is much more pronounced [16–19, 37], e.g. Rose et al. estimated a 40% reduction in symptomatic cases through increased vaccination coverage among children and utilization of LAIV, whereas Baguelin et al. suggest a 55% incidence reduction in a scenario with 50% vaccination coverage among 116 year old children.
Those deviating results can be explained by the assumption of a phenomenological transmission rate. Our alternative model scenario (S2) employing the massactionprinciple – like all previous influenza models – also predicts a 56.5% reduction in IMAARI when vaccinating 40% of 210 yearold children (Table 3). In the case of phenomenological transmission (base model) the indirect vaccination effects are much less pronounced, because small reductions in the prevalence of infection have an even smaller effect on the resulting force of infection, whereas when assuming massactiontransmission these effects are proportional. From a statistical perspective, the phenomenological model is more valid as indicated by its higher marginal likelihood of the data compared to the massactiontransmission model, which yields a considerably worse fit especially for weak waves such as those in the seasons 2007/08 or 2010/11 (see Additional file 1: Figure S5). Furthermore, the phenomenological modelling approach has a transmissiondynamic justification, since influenza prevalence is often spatially clustered across Germany [21]. This disagrees with the massactionprinciple assuming a homogeneous distribution of infectious people in the population.
Expanding the vaccination program to 11 to 17 year old children leads to an additional 5.1% decrease in IMAARI (Table 2), whereas increasing the vaccination uptake among children aged 2 to 10 years to 60% leads to an additional 5.2% reduction. Thus, allocating resources to the agegroup with the highest contact rate is worth considering. Analogously, administration of a vaccine providing higher VE among children — such as LAIV in our model example — also enhances the vaccination impact. Again note, that the assumption of LAIV granting improved protection must be challenged in the face of new evidence from the US [15]. An additional benefit would also have been provided by quadrivalent vaccines, since in seasons 2005/06, 2007/08 and 2008/09 we observed a considerable circulation of the Blineage that was missing in the vaccine. Due to indirect effects routine childhood influenza vaccination can prevent more IMAARI compared to a vaccination strategy among elderly even with 100% uptake. When taking required vaccine doses into account (Fig. 4), childhood vaccination is much more efficient compared to targeted vaccination of highrisk groups as it yields much lower NNVs to prevent one IMAARI. This results primarily from indirect effects affecting the overall population, which are most pronounced when targeting younger age groups (Figs. 3, 5).
Limitations
Although the model fit to the IMAARI and virological data is overall satisfactory, there is one model mismatch in 2003/04 where the virological data detected A(H3N2) as the only circulating subtype whereas the model additionally predicts circulation of BYamagata. This is presumably a result of the IMAARI data forming a double peak which was interpreted as two overlapping waves by the model. In fact, this double peak might have originated from a reduced midseason transmission activity during school vacation in February 2004, which cannot be captured by the model, or maybe from inaccuracies in the estimation of the IMAARI.
Although our stratification approach allows for the necessary seasonal heterogeneity when modelling influenza transmission, our model is missing some explicit linkage between subsequent seasons, e.g. acquired immunity could be carried over into the following seasons as implemented by Goeyvaerts et al. [36].
We modelled vaccine administration being completed with beginning of each season, which does not account for possible late administration, e.g. in February or March, or for the delay in developing vaccineinduced immunity. This approach thus slightly overestimates the impact of vaccination in our model. However, in the past most vaccines were in fact administered much prior to the beginning of influenza season such that we found this model simplification to have marginal effects [24].
We assumed a fixed probability of 0.67 for developing symptoms for all age groups and subtypes although this probability likely varies accordingly. However, Carrat et al. [49] did not detect any significant differences at least with respect to subtype and estimating this probability as an age and subtypestratified model parameter would have led to identifiability problems since it functions in a similar way as especially the probability for seeking medical care.
Finally, the focus of our model predictions was on the number of IMAARI by age groups. Thus, our predictions do not account for different hospitalization or mortality rates which could provide additional insight especially when comparing the overall benefits of vaccination programs targeting different age groups. Moreover, our model did not explicitly account for individuals with an underlying disease, that might suffer from more severe influenza and among which the vaccine might be less effective. However, the highest proportion of influenzaassociated complications and mortality is usually seen in young children and the elderly, the same agegroups where our model also predicted the strongest relative reduction in IMAARI (Fig. 3).
Conclusion
The present study provides an expansion of existing modelling studies investigating the impact of vaccinating the main spreaders of influenza – children. We employed a stratified model approach that simultaneously addressed seasonal variability but also consistent patterns within influenza transmission dynamics. Among the revisited methodological approaches, we found that a switch from the established massactionprinciple to a phenomenological approach has considerable impact on the predicted effects of childhood vaccination programs. Thus, more insight on the accurate modelling of transmission rates in compartment models or, alternatively, the application of more realistic individualbased models as in Eichner et al. are required in the future [35]. Even with our more conservative approach, our model predicted that childhood vaccination could considerably reduce influenza infections on the ambulant level with a much lower numberneededtovaccinate than the currently implemented vaccination strategy.
Abbreviations
 AGI:

German working group of influenza
 MAARI:

medically attended acute respiratory infections
 IMAARI:

influenza attributable MAARI
 TIV:

trivalent influenza vaccine
 QIV:

quadrivalent influenza vaccine
 LAIV:

live attenuated influenza vaccine
 WHO:

World Health Organization
 IMOVE:

Influenza  Monitoring Vaccine Effectiveness
 VE:

vaccine effectiveness
 NNV:

number needed to vaccinate
References
 1.
Neuzil KM, Hohlbein C, Zhu Y. Illness among schoolchildren during influenza season: Effect on school absenteeism, parental absenteeism from work, and secondary illness in families. Arch Pediatr Adolesc Med. 2002;156(10):986–91.
 2.
Monto AS, Sullivan KM. Acute respiratory illness in the community. Frequency of illness and the agents involved. Epidemiol Infect. 1993;110(1):145–60.
 3.
WHO. Influenza vaccines—WHO position paper. Wkly Epidemiol Rec. 2005;33:279–87.
 4.
ECDC. ECDC Vaccination schedule. [cited 2016; Available from: http://vaccineschedule.ecdc.europa.eu/Pages/Scheduler.aspx.
 5.
Michiels B, et al. A systematic review of the evidence on the effectiveness and risks of inactivated influenza vaccines in different target groups. Vaccine. 2011;29(49):9159–70.
 6.
Jefferson, T., et al., Vaccines for preventing influenza in healthy adults, in Cochrane Database Syst Rev. 2010.
 7.
Reichert TA, et al. The Japanese Experience with Vaccinating Schoolchildren against Influenza. N Engl J Med. 2001;344(12):889–96.
 8.
Grohskopf LA, et al. Prevention and Control of Influenza with Vaccines: Recommendations of the Advisory Committee on Immunization Practices. Morb Mortal Wkly Rep. 2015;64(30):818–25.
 9.
Pebody, R.G., et al., Uptake and impact of vaccinating school age children against influenza during a season with circulation of drifted influenza A and B strains, England, 2014/15. Euro Surveill, 2015. 20(39).
 10.
Loeb M, et al. Effect of influenza vaccination of children on infection rates in hutterite communities: A randomized trial. JAMA. 2010;303(10):943–50.
 11.
Spier R, et al. 3rd Vaccine Global Congress, Singapore 2009. Assessing Herd Immunity in the Elderly Following the Vaccination of School Children with Live Attenuated Trivalent Influenza Vaccine (LAIV): A CountyLevel Analysis. Proc Vaccinol. 2010;2(1):92–100.
 12.
Keipp Talbot H, et al. Influenza in Older Adults: Impact of Vaccination of School Children. Vaccine. 2009;27(13):1923–7.
 13.
Rhorer J, et al. Efficacy of live attenuated influenza vaccine in children: A metaanalysis of nine randomized clinical trials. Vaccine. 2009;27(7):1101–10.
 14.
Flannery B. LAIV vs IIV effectiveness: Summary of evidence since 2009. In: Presented at: Centers for Disease Control and Prevention Advisory Committee on Immunization Practices (ACIP) Meeting; 2016 Jun 22–23; Atlanta, GA; 2016.
 15.
Centers for Disease Control and Prevention. ACIP votes down use of LAIV for 20162017 flu season. 2016.
 16.
Baguelin M, et al. Assessing Optimal Target Populations for Influenza Vaccination Programmes: An Evidence Synthesis and Modelling Study. PLoS Med. 2013;10(10):e1001527.
 17.
Vynnycky E, et al. Estimating the impact of childhood influenza vaccination programmes in England and Wales. Vaccine. 2008;26(41):5321–30.
 18.
Rose MA, et al. The epidemiological impact of childhood influenza vaccination using liveattenuated influenza vaccine (LAIV) in Germany: predictions of a simulation study. BMC Infect Dis. 2014;14:4040.
 19.
Meeyai A, et al. Seasonal Influenza Vaccination for Children in Thailand: A CostEffectiveness Analysis. PLoS Med. 2015;12(5):e1001829.
 20.
an der Heiden M, et al. Estimates of Excess Medically Attended Acute Respiratory Infections in Periods of Seasonal and Pandemic Influenza in Germany from 2001/02 to 2010/11. PLoS One. 2013;8(7):e64593.
 21.
Buda S, et al. Bericht zur in Deutschland Saison 2013/14. Berlin: Robert KochInstitut; 2014.
 22.
Buchholz, U., A. Grüber, and B. Schweiger. Abschlussbericht der Influenzasaison. Robert KochInstitut
 23.
STIKO. Statement of the German Standing Committee on Vaccination at the RKI: Recommendations of the Standing Committee on Vaccination (STIKO) at the Robert Koch Institute. Epidemiol Bull. 2015;34:327–62.
 24.
Rieck T, et al. Vaccination coverage among children in Germany estimated by analysis of health insurance claims data. Hum Vaccin Immunother. 2014;10(2):476–84.
 25.
Jefferson, T., et al., Vaccines for preventing influenza in healthy adults. Cochrane Database Syst Rev, 2010 (7).
 26.
Jefferson, T., et al., Vaccines for preventing influenza in healthy children. Cochrane Database Syst Rev, 2008(2).
 27.
Kissling, E., et al., Influenza vaccine effectiveness estimates in Europe in a season with three influenza type/subtypes circulating: the IMOVE multicentre casecontrol study, influenza season 2012/13. Euro Surveill. 2014, 19(6).
 28.
Kissling E, et al. IMOVE MultiCentre Case Control Study 201011: Overall and Stratified Estimates of Influenza Vaccine Effectiveness in Europe. PLoS One. 2011;6(11):e27622.
 29.
Kissling, E., et al., "IMOVE" towards monitoring seasonal and pandemic influenza vaccine effectiveness: lessons learnt from a pilot multicentric casecontrol study in Europe, 20089. Euro Surveill, 2009. 14(44).
 30.
Kissling, E., et al., Low and decreasing vaccine effectiveness against influenza A(H3) in 2011/12 among vaccination target groups in Europe: results from the IMOVE multicentre casecontrol study. Euro Surveill. 18(5).
 31.
Valenciano M, et al. Estimates of Pandemic Influenza Vaccine Effectiveness in Europe, 2009–2010: Results of Influenza Monitoring Vaccine Effectiveness in Europe (IMOVE) Multicentre CaseControl Study. PLoS Med. 2011;8(1):e1000388.
 32.
Valenciano M, et al. The European IMOVE Multicentre 2013–2014. CaseControl Study. Homogeneous moderate influenza vaccine effectiveness against A(H1N1)pdm09 and heterogenous results by country against A(H3N2). Vaccine. 2015;33(24):2813–22.
 33.
Langley JM, et al. Immunogenicity and safety of an inactivated quadrivalent influenza vaccine candidate: a phase III randomized controlled trial in children. J Infect Dis. 2013;208(4):544–53.
 34.
Tricco AC, et al. Comparing influenza vaccine efficacy against mismatched and matched strains: a systematic review and metaanalysis. BMC Med. 2013;11(1):1–19.
 35.
Eichner M, et al. 4Flu  an individual based simulation tool to study the effects of quadrivalent vaccination on seasonal influenza in Germany. BMC Infect Dis. 2014;14:365365.
 36.
Goeyvaerts N, et al. Estimating dynamic transmission model parameters for seasonal influenza by fitting to age and seasonspecific influenzalike illness incidence. Epidemics. 2015;13:1–9.
 37.
Pitman RJ, White LJ, Sculpher M. Estimating the clinical impact of introducing paediatric influenza vaccination in England and Wales. Vaccine. 2012;30(6):1208–24.
 38.
Ultsch B, et al. Methods for Health Economic Evaluation of Vaccines and Immunization Decision Frameworks: A Consensus Framework from a European Vaccine Economics Community. PharmacoEconomics. 2016;34:227–44.
 39.
Gelman A, Hill J. Data analysis using regression and multilevel/hierarchical models. Cambridge: Cambridge University Press; 2006.
 40.
Finkenstädt BF, Grenfell BT. Time series modelling of childhood diseases: a dynamical systems approach. J R Stat Soc: Ser C: Appl Stat. 2000;49(2):187–205.
 41.
McCallum H, Barlow N, Hone J. How should pathogen transmission be modelled? Trends Ecol Evol. 2001;16(6):295–300.
 42.
Chowell G, et al. Characterizing the reproduction number of epidemics with early subexponential growth dynamics. J R Soc Interface. 2016;13
 43.
Kong, L., et al., Modeling Heterogeneity in Direct Infectious Disease Transmission in a Compartmental Model. Int J Environ Res Public Health, 2016. 13(3).
 44.
Mossong J, et al. Social Contacts and Mixing Patterns Relevant to the Spread of Infectious Diseases. PLoS Med. 2008;5(3):e74.
 45.
Goeyvaerts N, et al. Estimating Infectious Disease Parameters from Data on Social Contacts and Serological Status. J Royal Stat Soc Ser C Appl Stat. 2010;59(2):255–77.
 46.
Hermes J, et al. Lack of evidence for presymptomatic transmission of pandemic influenza virus A(H1N1) 2009 in an outbreak among teenagers; Germany, 2009. Influenza Other Respir Viruses. 2011;5(6):e499–503.
 47.
Eames K, et al. The impact of illness and the impact of school closure on social contact patterns. Health Technol Assess. 2010:14(34).
 48.
Van Kerckhove K, et al. The Impact of Illness on Social Networks: Implications for Transmission and Control of Influenza. Am J Epidemiol. 2013;178(11):1655–62.
 49.
Carrat F, 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.
 50.
Weidemann F, et al. Modelling the epidemiological impact of rotavirus vaccination in Germany – A Bayesian approach. Vaccine. 2014;32(40):5250–7.
 51.
Kass RE, Raftery AE. Bayes Factors. J Am Stat Assoc. 1995;90(430):773–95.
 52.
Hardelid, P., et al., Assessment of baseline agespecific antibody prevalence and incidence of infection to novel influenza A H1N1 2009. Health Technol Assess, 2011. 14(55).
 53.
Dudareva S, et al. Prevalence of Antibodies to 2009 Pandemic Influenza A (H1N1) Virus in German Adult Population in Pre and PostPandemic Period. PLoS One. 2011;6(6):e21340.
 54.
Basta NE, et al. Strategies for Pandemic and Seasonal Influenza Vaccination of Schoolchildren in the United States. Am J Epidemiol. 2009;170(6):679–86.
 55.
Bayer, C., et al., Internetbased syndromic monitoring of acute respiratory illness in the general population of Germany, weeks 35/2011 to 34/2012. Euro Surveill., 2014 19(4).
Acknowledgements
We would like to thank Brunhilde Schweiger for providing virological data, Joël Mossong and Kim Van Kerckhove for providing data on contact frequencies, and Matthias an der Heiden for advice on IMAARI data and model design.
Funding
The work of FW was financially supported by the accompanying surveillance research program of the Robert Koch Institute.
Availability of data and materials
Data on population counts, vaccine coverage rates and vaccine effectiveness are provided within the supplemental material. Data on estimated IMAARI and virological distribution are displayed in the supplemental material and are available from the authors upon request. RCode used for data analyses is also available from the corresponding author upon request.
Authors’ Contribution
FW, CR and OW conceived and designed the study. SB and UB processed the data. FW conducted all modelling work. FW drafted the manuscript. CR, SB, UB, BU and OW helped to revise the model and the manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable
Ethics approval and consent to participate
Not applicable
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Additional file
Additional file 1:
Supplementary information. (PDF 2648 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Weidemann, F., Remschmidt, C., Buda, S. et al. Is the impact of childhood influenza vaccination less than expected: a transmission modelling study. BMC Infect Dis 17, 258 (2017). https://doi.org/10.1186/s1287901723446
Received:
Accepted:
Published:
Keywords
 Influenza
 childhood vaccination
 transmission model
 NNV
 Bayesian inference