 Research article
 Open Access
 Published:
4Flu  an individual based simulation tool to study the effects of quadrivalent vaccination on seasonal influenza in Germany
BMC Infectious Diseases volume 14, Article number: 365 (2014)
Abstract
Background
Influenza vaccines contain Influenza A and B antigens and are adjusted annually to match the characteristics of circulating viruses. In Germany, Influenza B viruses belonged to the B/Yamagata lineage, but since 2001, the antigenically distinct B/Victoria lineage has been cocirculating. Trivalent influenza vaccines (TIV) contain antigens of the two A subtypes A(H3N2) and A(H1N1), yet of only one B lineage, resulting in frequent vaccine mismatches. Since 2012, the WHO has been recommending vaccine strains from both B lineages, paving the way for quadrivalent influenza vaccines (QIV).
Methods
Using an individualbased simulation tool, we simulate the concomitant transmission of four influenza strains, and compare the effects of TIV and QIV on the infection incidence. Individuals are connected in a dynamically evolving agedependent contact network based on the POLYMOD matrix; their agedistribution reproduces German demographic data and predictions. The model considers maternal protection, boosting of existing immunity, loss of immunity, and crossimmunizing events between the B lineages. Calibration to the observed annual infection incidence of 10.6% among young adults yielded a basic reproduction number of 1.575. Vaccinations are performed annually in October and November, whereby coverage depends on the vaccinees’ age, their risk status and previous vaccination status. New drift variants are introduced at random time points, leading to a sudden loss of protective immunity for part of the population and occasionally to reduced vaccine efficacy. Simulations run for 50 years, the first 30 of which are used for initialization. During the final 20 years, individuals receive TIV or QIV, using a mirrored simulation approach.
Results
Using QIV, the mean annual infection incidence can be reduced from 8,943,000 to 8,548,000, i.e. by 395,000 infections, preventing 11.2% of all Influenza B infections which still occur with TIV (95% CI: 10.711.8%). Using a lower B lineage cross protection than the baseline 60%, the number of Influenza B infections increases and the number additionally prevented by QIV can be 5.5 times as high.
Conclusions
Vaccination with TIV substantially reduces the Influenza incidence compared to no vaccination. Depending on the assumed degree of B lineage cross protection, QIV further reduces Influenza B incidence by 1133%.
Background
Before 1985, circulating Influenza B viruses were closely related to the precursor of the subsequent B/Yamagata (B/Yam) lineage [1]. Since 2001, viruses of the distinct B/Victoria (B/Vic) lineage have been detected in clinical samples from various parts of the world. The two Influenza B lineages are antigenically distinct [1–3]; genome sequence analysis revealed that they substantially differ as to the gene encoding the surface hemagglutinin [1, 2] and that they continue to diverge [4, 5]. Viruses of both B lineages now cocirculate together with the Influenza A viruses A(H3N2) and A(H1N1) [6]. The incidence of Influenza B infection varies substantially from year to year [6]. Surveillance data of the seasons from 2001/02 to 2010/11 show that Influenza B is responsible for 1 to 60% (mean 23%) of influenza cases in Europe and for <1 to 44% (mean 24%) in the United States, respectively [6]. Influenza B infections cause severe disease in all age groups [7] and their clinical symptoms and outcomes are similar to those of Influenza A infections [8].
Influenza vaccine design is intended to anticipate which variants will circulate in the subsequent transmission season. The antigenic composition of the vaccine is revised twice a year by the World Health Organization (WHO) and adjusted to the antigenic characteristics of circulating influenza viruses in the Northern and Southern hemisphere. Current trivalent influenza vaccines (TIV) contain antigens of A(H3N2) and A(H1N1), but of only one Influenza B lineage [9]. The efficacy of these influenza vaccines varies with the age of the vaccinees and with the degree at which the vaccines match the circulating influenza strains. In case of a good match, their efficacy can range from 70% to 90% in young adults, but it is lower in very young children and in elderly [9]. A recent metaanalysis of randomized controlled vaccine efficacy trials showed that TIV confers a 52% protection (95% confidence interval (CI): 19 to 72%) against the Blineage not contained in the vaccine, compared to 77% (95% CI: 18 to 94%) for the lineage which matches the vaccine [10]. Cocirculation, fluctuating prevalence, and differing geographical distributions of the two B lineages have challenged global recommendations for vaccine composition [3]. In the US, the Influenza B lineage in TIV failed to match the predominant circulating B lineage in five out of ten seasons between 2001 and 2011 [11]. In 2012, WHO recommended for the first time vaccine strains from both B lineages, paving the way for the development of quadrivalent influenza vaccines (QIV) which contain strains of two B lineages and two A subtypes [12].
A variety of deterministic and stochastic models have been proposed to address various aspects of seasonal [13–16] or pandemic influenza [17–21] or to predict the size and timing of seasonal waves [22–24]. We have decided to use an individualbased stochastic model which allows for maximum flexibility in modeling immunity and vaccination strategies [25] and, thus, circumvents common oversimplifications of deterministic models [26]. We have developed the simulation tool 4Flu to compare the effect of trivalent and quadrivalent vaccination on the incidence of influenza infection in Germany. 4Flu is based on an individualbased stochastic model, which allows for the independent and simultaneous transmission of the four influenza variants A(H1N1), A(H3N2), B/Vic and B/Yam. We use the official German agedistribution for simulation years which lie in the past, and official demographic predictions for the future [27]. Individuals are connected by a dynamically changing contact network which is synchronized with the German POLYMOD contact matrix [28]. The simulation model further allows for the introduction of new drift variants and for mismatched vaccine design against new drift variants. Special care has been taken to provide a realistic immunity model which allows for maternal protection of newborns, for booster effects and waning of existing immunity and for crossimmunizing effects of infections and vaccinations. Vaccination either uses TIV or QIV. Each simulation runs for several decades and reports the annual incidence of influenza infections. Simulations are typically repeated one thousand times to average out random effects.
Methods
We present an individualbased stochastic simulation model for the simultaneous and independent transmission of four influenza strains (A(H1N1), A(H3N2), B/Vic, B/Yam). Each simulation runs for 50 years, with each simulation year starting on July 1st: the first 10 years are used to initialize the contact network, the next 20 years for initializing agedependent infection and immunity patterns in the population, and the final 20 years for comparing vaccines.
Demography
As demographic changes of the population have been shown to have a major impact on the result of transmission models for seasonal influenza [13], we have made sure that the simulated demography is at all times identical to the population structure of Germany. Our simulated population corresponds to one thousandth of the German population (about 80,000 individuals). Every individual has an age and a risk status. For the years from 1993 to 2011, we use the observed age distribution of Germany, for later years, we use the official prediction “variant 1, W2” of Deutsches Statistisches Bundesamt [27]. Births and deaths occur throughout the year and are implemented such that the simulated population exactly matches the official demographic distribution at the end of each simulation year. If the size of a cohort increases from one year to the next, new individuals immigrate (see Appendix 1 for details). Depending on the age of the individuals, they are given a risk status which is used to determine the vaccination coverage and which can also be used to determine the risk of complications upon infection. As individuals age and the demography keeps changing, the individuals’ risk status is updated during the simulation to keep the specified risk fractions in the different age classes constant.
Contact network
For the transmission of the influenza infections, we have constructed a dynamic network which connects the individuals, using the German POLYMOD matrix as summary statistic for these connections. This contact matrix is based on an EU sponsored study in which individuals were asked to supply a contact diary for one day. Contacts were then summarized in 5year age classes [28]. To use this matrix for a population with changing demography, we connect each individual with randomly picked other individuals such that the expected number of these “outgoing” connections is proportionate to the values of the POLYMOD matrix (which depend on the ages of the individuals and their contacts). On average, each individual has a sum of 10 “outgoing” and “incoming” contacts. When individuals age, their contacts gradually change such that the population statistics of the agedependent contact distribution remain in agreement with the POLYMOD matrix at any time (see Appendix 2 for details). In order to initialize the contact structure, we start with the German demography of 1993 and chose random contacts of individuals within 5year age classes. We then run the simulation for 10 years without infection while keeping the age distribution constant, but allowing for births, deaths and immigrations as needed.
Initialization of immunity, infection and vaccination
After the initialization of the contact network, a random fraction of 30% of the population is set to be immune against the Influenza A subtypes and a random fraction of 40% is set to be immune against the Influenza B lineages. To transform this crude initial immunity into a realistic agedependent immunity pattern, transmission of influenza is simulated for 20 years before vaccination strategies are compared. As B/Vic did not circulate in Germany for at least 10 years before it reappeared in the 2001/02 season [29], no B/Vic infections are used in the simulations before the 2001/02 simulation year. During the initialization phase from 1993 to 2012, individuals are vaccinated with TIV, using the B lineage in the vaccine which was actually used in Germany in the given year. In the 20year evaluation period which starts with the simulation year 2013/14, individuals are either vaccinated with TIV or with QIV. As exactly the same individuals are given either one of the two vaccines on the same vaccination days, pairwise comparisons can be made between the TIV and the QIV scenario. A 50/50 random decision is made for each future simulation year to determine which B lineage will be used in TIV.
Natural history and transmission
For the transmission of the infection, it is not important whether a contact is labeled as “incoming” and “outgoing”; any connection of two individuals allows transmitting the infection either way. In order to calculate the probability of infecting a contact, the largest eigenvalue of the next generation matrix is calculated after initialization of the contact network. The transmission probability per contact is then given by dividing the basic reproduction number R _{ 0 } (chosen to be 1.575 for the baseline parameter setting) by the largest eigenvalue [30] (see Appendix 3 for details). The transmission probability is further assumed to be subject to seasonal fluctuation; its maximum around Christmas is 43% higher than the baseline, its minimum occurs on 21 June of each year [31]. Throughout the year, individuals are constantly exposed to a very low rate of infection “from the outside” of the population, irrespective of their age. Infected individuals first pass through a latent period of 2 days and then become infectious for 4 days (individuals below 18 years of age) or 2 days (older individuals), respectively [32]. Transmission of the infection occurs at a random time point during the infectious period.
Naturally acquired immunity
After infection, individuals acquire a temporary immunity which prevents further infections; this immunity is lost continuously at a constant rate. After losing immunity, individuals become fully susceptible again. If an immune individual gets in contact with an infectious one, or if an immune individual is vaccinated, the duration of his or her immunity can be extended by a booster vaccination or infection (see details in Appendix 4) without rendering the individual infectious. As simulation models have shown that crossimmunizing events caused by cocirculating B lineages can play an important role in the transmission dynamics of seasonal influenza [14, 26], we have also incorporated this feature in 4Flu: infection with any one of the two B lineages has a 60% probability to lead to a crossimmunizing event against the other B lineage [10]. We do not consider specific crossprotection between the two Influenza A strains or among Influenza A and B, but examine the possible effect of a shortterm unspecific protection in a sensitivity analysis. Newborn individuals are protected by maternal antibodies against the influenza strains against which their mother is immune (irrespective of whether her immunity was acquired by infection or vaccination). This maternal protection lasts for 2 to 4 months [33, 34] and is assumed to prevent infection. The dynamics of acquiring and losing immunity is schematically depicted in Figure 1.
Drift variants
After an average circulation time of 3.5 years (A(H3N2)) [35, 36] or 7 years (all other influenza viruses), a circulating variant is replaced by a new drift variant in the simulation. For the sake of simplicity, we assume that all variants of a given influenza subtype or lineage which are introduced into the population in a simulation year belong to the same drift variant, i.e. switching to a new drift variant always occurs at the beginning of the simulation year (on July 1st). We further assume that a variant will not be introduced again after it has been replaced by a new drift variant. A percentage of 60% randomly chosen individuals who are immune against the old variant are also immune against the new one, whereas the remaining 40% are assumed to be susceptible against the new drift variant (see Appendix 6 for details) [36].
Vaccination
Vaccinations are performed in October and November of each simulation year [37], whereby a percentage of individuals is vaccinated, using either TIV or QIV. As shown in Table 1, the vaccinated percentage depends on the age of the individuals, on their risk status and on their vaccination status [38]: individuals who were vaccinated in the previous simulation year have twice the chance of being vaccinated again, compared to previously unvaccinated individuals (see Appendix 5 for details). An agedependent percentage of previously susceptible vaccinees become immune after vaccination (vaccine efficacy) [39–41]. In case of TIV vaccination, we also consider the possibility of crossimmunizing effects against the B lineage which is not included in the vaccine: the agedependent vaccine efficacy is multiplied with 0.6 to determine what percentage of vaccinees is protected against the B lineage which is missing in the vaccine. Successful vaccination can also extend the duration of preexisting immunity (“boostervaccination”; see Appendix 5 for details). When a new drift variant has just been introduced, the vaccine may not match this variant (vaccine mismatch). In this case, the agedependent vaccine efficacy is multiplied by 0.6 (calculated from [39]). The decision of whether the vaccine is matched or mismatched to a new drift variant [42] is made identically for the TIV and the QIV scenario. Vaccine efficacy depends on the age of the vaccinee, but is assumed to be independent on the virus subtype or lineage. We use the same vaccine efficacy for TIV and for QIV (Table 1) [39–41], but distinguish between age classes and individual risk status [43]. For convenience, we use the same agedependent values of the vaccination coverage in every simulation year. In the first simulation year, a total of 26.8% of the population is vaccinated. Due to demographic changes (Figure 2) and due to a high vaccination coverage in the elderly (Table 1), the overall vaccination coverage reaches 29.4% at the beginning of the evaluation period and continues increasing until the end of the simulation when it reaches 33.4%. Vaccinationderived immunity is lost much more quickly than naturally acquired immunity [44].
Comparison of TIV and QIV
During the initialization period, only TIV is used as vaccine. During the evaluation period, each single simulation is run as a “TIV” and as a “QIV scenario”: (1) In order to keep the two scenarios closely matched, they use exactly the same contact network and demographic composition of the population. (2) Exactly the same individuals are vaccinated in each scenario, using either TIV or QIV at exactly the same time points. (3) If a drift event or a vaccine mismatch occurs, it occurs in the same simulation year in both scenarios. As a consequence, individuals can only differ in their epidemiologic and immunologic state with respect to the B lineages in the two scenarios. We calculate the cumulative numbers of infections for the 20 year evaluation period, using TIV or QIV and finally report the annual number of infections prevented by QIV (difference “TIVQIV” divided by 20). As simulations are strongly influenced by random numbers, at least 1,000 simulations are run for each parameter setting and arithmetic means of the prevented number of infections are calculated.
Model calibration
After an initial screening process, 1,000 to 2,000 simulation runs are performed for each proposed value of the basic reproduction number, starting with values from 1.5 up to 1.6, using a step size of 0.025. For each simulation, it is recorded, what percentage of individuals in the age class from 16 to 60 years experiences an influenza infection in the 2006/07 transmission period. These results are compared to the observation of Williams et al.[45]. As 2006/07 was labeled as a “normal” influenza season by the Arbeitsgemeinschaft Influenza [42], we used the median of the simulated infections for model calibration.
Sensitivity analyses
The basic reproduction number R _{ 0 } has been shown to be one of the most influential parameters which determine the outcome of influenza models [13]. Different ranges for R _{ 0 } have been proposed in the literature (0.92.1 [46], 1.21.3 [47], 1.21.4 [48], 1.22.3 [49], 1.31.7 [50], 1.31.9 [51], 1.41.6 [52], 1.63.0 [53]). Differences in these ranges may derive from the differing epidemiologic circumstances of the places and seasons in which the estimates were derived, yet they may also be due to different ways of estimating R _{ 0 }[30, 54, 55]. As it has been stressed to address uncertainties in the value of R _{ 0 }[56], we perform a sensitivity analysis in which we screen the parameter range from 1.2 to 2.0, using 1,000 simulations for each value. Further sensitivity analyses address the influence of the duration of naturally acquired immunity and of the degree of B lineage cross protection.
Results
Figure 2 shows how the contact matrix changes over time: Figure 2a shows the state of the contact network before initialization: the 5year age classes of the German POLYMOD matrix are slightly disturbed by the German demography, but are still clearly visible. After about five years of initialization, the contact structure has changed to mostly allocate contacts to individuals of the same age, although (by construction of the algorithm) the summary statistics of the number of “outgoing” contacts per 5year age group remain identical to the POLYMOD matrix. After ten years of network initialization, the contact network has converged to a steady state (Figure 2b). After that time, the immunologic and epidemiologic initializations are started and the demography of the simulated population is allowed to change according to the official German demography. The resulting age distributions of the German population are shown on the right hand sides of the contact matrix visualizations in Figure 2bf.
As the number of children keeps declining and the number of elderly grows in the German population, the frequency of contacts becomes less focused on children and young adults and tends to be more spread out over the ages. In 2012, contacts are most frequent in juveniles and then decline over age (Figure 3b). Taking into consideration the longer duration of infectiousness of children and juveniles and considering the susceptibility status of their contacts, we can calculate how many individuals are infected on average by a single infectious individual of a given age as shown in Figure 3c. The median time until an infected individual passes on the infection (generation period) is 3.3 days in our simulations.
Our model is calibrated to observed data from Germany. As our simulation results vary strongly (coefficient of variation CV = 91.5% for R _{ 0 } = 1.575), we use the median of the simulated incidence when evaluating simulation results: for each value of the basic reproduction number R _{ 0 } between 1.5 and 1.6 (step size 0.025), we ran 1,000 to 2,000 simulations and calculated the median of the infection incidence in the 2006/07 transmission season for young adults. A basic reproduction number of R _{ 0 } = 1.575 yields a median infection incidence of nearly exactly 10.6% as was observed among German health care workers and controls in that season [45]. Figure 4a shows the large fluctuation of annual influenza waves in a single simulation in which TIV is used throughout, whereby all infections with the two A strains and the two B lineages are added up. Figure 4b selectively shows the incidence with B/Yam for the same simulation. Years in which B/Yam is not contained in TIV are shown as dark grey bars. New B/Yam drift variants are introduced in 2015, 2022, 2028 and 2031, shown as vertical dashed lines. In 2028, the B/Yam lineage in the vaccine does not match the newly introduced drift variant, indicated by a light grey background. During two seasons, large B/Yam outbreaks occur, dwarfing the smaller waves in the other seasons. Starting with the year 2013, simulations branch into a TIV scenario and a QIV scenario (Figure 4c). The cumulative numbers of infections are calculated in each of these scenarios and are used to calculate the average annual number of infections prevented by QIV as compared to TIV.
Averaged annual simulation results for the baseline set of parameter values are given in Table 2. As the simulated population is exactly one thousandth of the German population, these results have to be multiplied with 1,000 to obtain estimates for the whole of Germany. The annual number of influenza infections in Germany can be reduced from 8,943,000 to 8,548,000, i.e. by 395,000 infections per year without increasing the vaccination coverage if QIV is used instead of TIV. The prevented infections constitute 4.3% (95% confidence interval (CI): 4.14.5%) of all influenza infections which still occur under TIV vaccination. The simulated fraction of Influenza B infections among all infections varies widely from one season to the next; even if averaged over the 20 years of each simulation, the 95% reference interval for the percentage of Influenza B infections among all influenza infections ranges from 24.6 to 51.3% (median 37.3%) in the TIV scenario, and from 21.3 to 49.1% (median 34.3%) in the QIV scenario. Switching from TIV to QIV prevents on average 11.2% (CI: 10.711.8%) of all Influenza B infections which still occur under TIV vaccination. In 53% of 10,000 simulation years, the B lineage not contained in TIV caused more infections than the B lineage which was contained in TIV.Additionally to lowering the expected number of infections, QIV furthermore reduces the number of seasons with rather extreme Influenza B outbreaks (Figure 5).
The risk that in Germany more than 10 million Influenza B infections occur in a single transmission season is about 13.4% for the TIV scenario and 11.2% for the QIV scenario (Figure 5a), meaning that the expected duration between such extreme events can be increased from 7.5 years (TIV) to 8.9 years (QIV). For 15 million infections per year, the difference is even more pronounced (21.3 years for TIV, 29.2 years for QIV). The risk of strong Influenza B seasons can be reduced by up to 27% if QIV is used (Figure 5b).
We address parameter uncertainty in various sensitivity analyses: the basic reproduction number, R _{ 0 }, is varied from 1.2 to 2.0 (baseline value: 1.575) while the values of the other parameters are kept at their baseline values. As R _{ 0 } summarizes the contagiousness of the disease, the average annual number of infections in Germany grows with increasing R _{ 0 } from 3,354,000 to 12,973,000 per year in the TIV scenario, and from 3,029,000 to 12,587,000 in the QIV scenario (Figure 6ac and Appendix 8); the average number of infections annually prevented by QIV rises from 325,000 infections for R _{ 0 } = 1.2, to about 386,000 infections for R _{ 0 } = 1.5, and, thereafter, remains about constant on a level of 395,000 infections. Due to the rising denominator, the percentage of prevented infections falls with increasing R _{ 0 } from 9.5 to 3.0% (all infections), or from 21.6 to 7.8% (Influenza B infections).
In a further sensitivity analysis, the rate at which immunity is lost is varied such that the duration of naturally acquired immunity lasts on average 4 to 12 years against A(H1N1) and the B lineages, and for 3.2 to 7.1 years against A(H3N2) (see Appendix 4 for details). The average annual number of infections drops with increasing immunity duration from 11,267,000 to 6,735,000 infections per year in the TIV scenario, and from 10,751,000 to 6,453,000 in the QIV scenario (Figure 7ac and Appendix 8). The average number of infections annually prevented by QIV correspondingly drops from 517,000 (4.6% of all infections, or 11.7% of influenza B infections) to 282,000 infections (4.1% of all infections, or 10.7% of influenza B infections).
Following the parameter choice of Vynnycky et al.[31], we perform a sensitivity analysis in which we only increase the duration of naturally acquired Influenza B immunity from the baseline value of 6 years to 12 years while keeping the immunity duration of Influenza A at the baseline value. This reduces the annual number of infections from baseline 8,943,000 to 8,098,000 in the TIV scenario, and from baseline 8,548,000 to 7,814,000 in the QIV scenario. The average number of infections annually prevented by QIV correspondingly drops from baseline 395,000 (4.3% of all infections, or 11.2% of all Influenza B infections) to 284,000 infections (3.5% of all infections, or 10.8% of all Influenza B infections). Increasing the duration of Influenza B immunity shifts the 95% reference intervals of the percentage of Influenza B infections among all influenza infections towards lower values, reaching 16.346.6% (median 30.5%) for the TIV scenario, and 13.743.3% (median 27.5%) for the QIV scenario.
Another highly influential parameter is the degree of B lineage cross protection: in a sensitivity analysis, we increase the percentage of individuals who are immunized against an Influenza B lineage when they are successfully vaccinated or infected with the other B lineage, using values from 0% up to the baseline value of 60% (Appendix 8). The higher the probability of B lineage cross protection is, the fewer Influenza B infections occur (Figure 8ac). Whereas for the baseline value, the median percentage of Influenza B infections among all infections is 37.3% (TIV scenario), it increases to 54.1% if no B lineage cross protection is considered. Accordingly, the number of infections which can annually be prevented by QIV increases from 395,000 (baseline value) up to 2,180,000 infections if lower values are considered. If the degree of B lineage cross protection decreases over time from initially 60% to a lower value which is reached at the end of the simulation, the additional benefit of using QIV increases from the baseline result of 395,000 infections per year (obtained for a constant level of 60% B lineage cross protection) to 559,000 (decreasing the cross protection from 60% to 50%), to 706,000 (40%) and to 853,000 (30%), respectively (average values of 1,000 simulations for each parameter set). Our model does not explicitly consider cross immunity between the A strains or among Influenza A and B, but we have run a sensitivity analysis where we consider a shortterm unspecific immunity which prevents infection with any influenza virus within the first weeks following an infection. As shown in Appendix 8, such an indirect protection, which may last for one to four weeks, has only a marginal effect on the results. If anything, it may slightly increase the difference between the QIV and the TIV scenario.
In order to keep the age distribution of “outgoing contacts” synchronized to the POLYMOD contact matrix despite demographic changes, contacts have to be removed or added daily. Although we keep the number of changes to a minimum, only 0.36% of the childhood contacts of 10 year old children at the beginning of the simulation are still among their contacts at the end of the simulation when they are 50 years old. We have examined the influence of contact turnover on our results in a sensitivity analysis in which we have either increased the contact turnover by replacing at random 10% of all contacts each year, or decreased contact turnover by updating contacts only once a year instead of daily (resulting in larger and longer lasting deviations of the simulated agedependent contact distribution from the POLYMOD matrix). The effect of contact turnover on the number of infections additionally prevented by QIV is marginal (Appendix 8).
Discussion
Our results indicate that annually an average number of 395,000 influenza infections (representing 11.2% of all Influenza B infections) can additionally be prevented without increasing vaccination coverage by switching from TIV to QIV. Although the composition of TIV during the evaluation period of our simulations is determined at random, our finding that in 53% of seasons, the B lineage not contained in TIV causes more infections than the B lineage which actually is contained in the vaccine is practically identical to findings from the US (the Influenza B lineage in TIV failed to match the predominant circulating B lineage in five out of ten seasons between 2001 and 2011) [11]. It had to be expected that this percentage should even be higher if the vaccination coverage was higher, yet this finding frequently is regarded as an indicator of having chosen the wrong B lineage in TIV.
As QIV is designed to additionally prevent Influenza B infections, our results strongly depend on the assumptions made on the transmissibility and on the duration of acquired immunity against Influenza B. Cross immunization after infection with an Influenza B lineage or after successful vaccination with TIV increases Influenza B immunity in the population. As we show in sensitivity analyses (Figure 8ac and Appendix 8), our baseline value of 60% lineage cross protection leads to a highly conservative estimate of the effects of QIV. Using lower cross protection values, the additional benefit of using QIV instead of TIV can be up to 5.5 times as high as the baseline result of preventing 395,000 infections per year (Figure 8ac). As the antigenic overlap of the two B lineages is clearly disadvantageous to their transmission, mutations which increase the genetic difference between the two lineages and, thus, reduce the degree of cross protection should be selected for. To check the effect of this hypothetical antigenic drift, we have continually lowered the degree of B lineage cross protection over time in a sensitivity analysis. If the degree of B lineage cross protection drops from initially 60% to finally 30%, the additional benefit of using QIV increases up to an average of 853,000 infections per year. It must be noted that we restrict crossimmunizing vaccine effects to TIV and omit them when individuals receive QIV, because it has been shown that increasing antigen amount does not enhance the efficacy of influenza vaccines [57, 58] and because there is no convincing evidence that increasing the amount of antigen enhances the efficacy of influenza vaccine [59].
The simulated absolute and relative numbers of Influenza B infections fluctuate strongly from year to year. Even if they are averaged over 20 simulation years, the percentage of Influenza B infections varies from 24.6 to 51.3% (95% reference interval in the TIV scenario). This range is contained in the range of reported values on clinical Influenza B cases, yet the median percentage of Influenza B infections (37.3%) is higher than the reported average of 23% of clinical cases who are caused by Influenza B [6]. This difference between the simulated percentage of Influenza B infections and the reported percentage of Influenza B cases may partly be due to a different fraction of Influenza A and B cases being reported: it has been shown that Influenza A infections more frequently lead to high fever which may lead to more doctoral visits and finally to more case reports [60]. During a mixed A(H3N2) and B season in the Netherlands, the proportion of laboratory confirmed Influenza A was higher in hospital samples than in GP sentinels [61]. This indicates an underestimation of the incidence with types which tend to cause mild illness when laboratory data include hospital samples. With increasing severity, the percentage of strains which cause milder infections should increasingly be underestimated (e.g. from diagnosis of acute respiratory infections, ARI, to diagnosis of influenza like illness, ILI, or hospitalization). Influenza infections are diagnosed by using Flurapid tests. Most of these tests can identify Influenza A and B, but the proportion of selective tests which only test for Influenza A is unknown. This might cause a slight underestimation of the true Influenza B incidence. Vynnycky et al.[31] use in their simulation study different immunity durations for Influenza A (6 years) and Influenza B (12 years). We have addressed this in a sensitivity analysis in which we double the duration of Influenza B immunity (unlike Vynnycky et al. we use the resulting duration of 12 years only for naturally acquired immunity and not for immunity conferred by vaccination). As a consequence, the median percentage of Influenza B infections among all influenza infections drops to 30.5% in the TIV scenario, and the average number of infections prevented by QIV drops to 284,000 per year (preventing 3.5% of all influenza infections, or 10.8% of all Influenza B infections, respectively). If, on the other hand, we assume that the B lineage cross protection is less than 60%, the number of Influenza B infections increases, and, consequently, QIV can prevent more infections (Figure 8ac). Without any B lineage cross protection, we expect to see more Influenza B than A infections (54.1 vs. 45.9%) in the TIV scenario, because TIV only immunizes against one B lineage, but against both A subtypes. In that extreme case, QIV would additionally prevent over two million infections per year. Although the annual incidence of infections strongly depends on R _{ 0 }, the number of prevented infections is nearly constant for R _{ 0 } > 1.5.
The seasonal influenza waves depicted in Figure 4a show a much higher variability than in the countrywide statistics published by the Robert Koch Institute [42]: the coefficient of variation of the simulated annual number of influenza infections is in the order of 100% (TIV scenario). This high variability mainly results from using a population size of only one thousandth of the German population: geographic differences in immunity patterns and transmission reduce the variability of the reported influenza incidence. When adding up the annual incidences of ten randomly picked simulations (to mimic the combined result of ten locations with independent transmission), the coefficient of variation reduces to about 30%. In reality, much more locations ought to be combined, but their epidemiology must not be assumed to be independent.
Various contact networks have been proposed to study the spread of influenza [62–65]. In our simulations, we use a novel contact network which is based on contact information gathered for the German population [28], yet which allows for dynamic changes during the ageing process of the population. As was shown in a sensitivity analysis, the network turnover had only negligible effect on the difference between the TIV and QIV scenario. It has been suggested to use bipartite networks which explicitly consider the place of contact (school, work place, household) [66]. Although the POLYMOD data include some information on the place of the contact [28], many additional assumptions would be required to explicitly consider contact place in our simulations. Some of the contacts which regularly occur at school or work place [67] are covered by the agedependency of contacts in the underlying POLYMOD matrix [28], but our contact network may be too rigid to allow for random encounters (e.g. during public transportation).
In our simulation studies, we make several simplifying assumptions: (1) the basic reproduction number is assumed to be equally high for all four influenza strains; differences between Influenza A and B are only due to the higher drift frequency of A(H3N2), to the missing B lineage in TIV, and to B lineage cross protection. (2) Infected children are as contagious as adults, yet they have a longer infectious period (cf. Table 1) and have more contacts than adults  specifically with other children, resulting in a larger number of secondary infections (Figure 3c). (3) We apply the same reduction factor 0.6 to the vaccine efficacy which crossprotects against the B lineage which is not contained in TIV, and to a vaccine design mismatch; likewise, we assume that only 60% of previously immune individuals are immune against a newly introduced drift variant. Tricco et. al.[10] report for vaccinees from 18 to 49 years of age that the vaccine efficacy against Influenza B is 77% if the lineage is contained in TIV and 52% if it is not contained in TIV, respectively. The ratio 52/77 = 0.675 is similar to our ratio of 0.6, calculated from a Cochrane review for any type of mismatch. Langley et al. show in a vaccination study with children that the seroconversion rate with respect to the B lineage which is not contained in TIV is half as high as the seroconversion rate with respect to the lineage which is contained in the vaccine [59] (cf. Appendix 6). Using a B lineage cross protection of only 50%, leads to a larger than baseline average annual number of infections prevented by QIV (Figure 8c). (4) School holidays have been shown to have a strong impact on the transmission of influenza [62], but we did not want to further complicate our model by changing contact patterns during such time periods, as different school holiday periods are used in different German states and as they change from year to year. We also do not distinguish between week days and weekends in our simulation. Using averaged social contacts throughout the week has been shown not to affect the major outcome of influenza simulation studies [63]. (5) Although vaccination coverage in Germany varies strongly between different locations and may also slightly vary over time [68], we use a constant agedependent vaccination coverage throughout our simulations [38] which is assumed to represent a countrywide average. (6) We do not explicitly consider human travel in our simulations [20, 21, 47], but we apply a very small “external” infection rate which is present throughout the year. (7) Individuals who were vaccinated against A(H1N1) pdm09 in 2009 were more likely to be revaccinated in 2010 with an odds ratio (OR) of 6.02 in The Netherlands [69], yet in our simulations, the OR of revaccination ranges from 2.2 to 6.5 (for details, see Appendix 5). If the reported OR is also representative for nonpandemic years in Germany, too few individuals are revaccinated in our simulations and, thus, the effects of vaccinations may be overestimated. (8) We use the same average immunity duration after booster infections as after infections of previously susceptible individuals as it is not known whether this immunity duration should be larger or smaller than after infection. (10) Although great care has been given to use realistic assumptions on the immunity acquired by infection, the parameters associated with immunity are quite uncertain. We have addressed this uncertainty in various sensitivity analyses, which explore the effect of the duration of immunity and the degree of heterologous immunity between the B lineages. It has been difficult to obtain quantifiable evidence for crossimmunizing effects between the A strains or between Influenza A and B, yet an indirect protection may be caused by an increased production of interferons, tumor necrosis factor alpha and beta and other cytokines as has been reported after infection with influenza and various other viral infections [70–73]. As we did not want to incorporate noninfluenza infections, we have restricted the study on the effects of unspecific protection to a sensitivity analysis which showed only a minor effect of indirect protection (Appendix 8).
Conclusions
4Flu is a novel simulation tool which describes the spread of influenza viruses using a dynamically evolving contact network, reflecting German demography, with different vaccinations strategies. The results indicate that the current vaccination with TIV in Germany, although covering less than 30% of the total population [38], substantially reduces the number of influenza infections compared to no vaccination. Replacing TIV with QIV would further reduce Influenza B incidence by 395,000 to 2,180,000 influenza infections per year, depending on the degree of B lineage cross protection.
Appendix 1 Dynamic changes in demographic structure and risk status
The observed and predicted demographic distributions of Germany are given in 1year cohorts for every year until 2060, together with annual birth rates [27]. Our simulations start with the demographic distribution of 1993. Each individual is given a random birth day; in each simulation year, individuals are added according to the predicted number of births. If the number of individuals of a cohort declines from one year to the next, individuals are chosen to die during the simulation year; if it increases, individuals are introduced as “immigrants”. Birth, death and immigration take place at random times during the year. On each birthday, an individual’s age grows by one year. Individuals are also given a risk status which influences their probability of being vaccinated and which may influence the clinical outcome of an infection. Some individuals have an increased risk when they are born, others acquire a risk status over time. In order to avoid sudden transitions from one age group to the next one, we assume that the percentage with risk changes linearly within each age group, whereby we keep the total percentage in the age group identical to the literature value. Due to births, ageing and death, the risk fractions in the simulations keep drifting away from the specified values given in Table 1. At each birthday, it is checked whether there are currently too many or too few individuals with risk status in the cohort of the individual who is now one year older; if needed, the individual’s risk status is changed in order to reduce this deviation.
Appendix 2 Dynamic changes in contact structure
In order to transmit the infection from one individual to another one, individuals are connected in a dynamic contact network. The average number of contacts depends on the age of the individual and on the ages of the contact persons, and is given by the POLYMOD matrix [28]. Before initializing the contact network, each individual “initiates” on average 5 contacts with others, whereby the age distribution of these contacts is proportional to the numbers given by the POLYMOD matrix. These contacts are classified as “outgoing contacts” for the source individual and as “incoming contacts” for the contacted individual (i.e. on average, every individual has 10 contacts). For each newborn individual, a random person of the agegroup from 20 to 40 years is chosen (denoted as the “mother”) who initiates a contact to the newborn individual. If the mother is connected to other individuals (which usually is the case), the newborn initiates a contact to the individual whose age is closest to that of the mother; thus newborns are initially connected in triangles. Immigrants are incorporated into the existing contact network in the same way as newborns. Due to births, immigrations, ageing and deaths, the contact distribution of “outgoing” contacts keeps drifting away from the values specified by the POLYMOD matrix. In order to correct this, the distribution of “outgoing” contacts is compared at the end of each simulation day with the POLYMOD matrix assuming a population average of 5 “outgoing” contacts per individual. If too many contacts between two age groups exist, superfluous contacts are chosen at random and removed; if too few contacts exist, random individuals are chosen from the required age groups and connected. In a sensitivity analysis this “repair mechanism” is only performed annually instead of daily (resulting in lower network turnover); or an additional random percentage of contacts is removed and replaced by new random contacts (resulting in higher network turnover). In order to let the contact network adjust to the demographic structure of Germany, we run the simulation for 10 years (“contact network initialization”) without changing the demographic structure of Germany before introducing immunity and infections in the simulations.
Appendix 3 Transmission dynamics and basic reproduction number
An infected individual does not necessarily infect all of his or her contacts; the expected number of secondary infections in a completely susceptible population is given by the basic reproduction number R _{ 0 }. In order to calculate the infection probability per contact from the basic reproduction number, the “next generation matrix” is calculated which combines the contact information in the population and the duration of contagiousness [32]. As we start the simulations with 80,976 individuals, the contact network can be described by an 80,976 × 80,976 matrix which is mostly filled with zeros. The next generation matrix is defined by the product of the network matrix and the average duration D _{ a } of contagiousness (which depends on the age a of the infectious individual). The largest eigenvalue E of this matrix is calculated by power iteration (Von Mises iteration). The probability that a given contact of an infectious individual of age a is infected by that individual is then given by p _{ a } = R _{ 0 } D _{ a } /E[30]. During the first 5 years of the initialization of the contact structure, the largest eigenvalue of the next generation changes markedly, but after that it remains stable, indicating that the contact network has converged to a steadystate. As influenza transmission is assumed to vary seasonally, we multiply the transmission probability with a seasonal factor: p _{ a } (1 + 0.43cos(2π(t173)/365)). The seasonality is chosen such that the transmission reaches its maximum around Christmas where it is 43% higher than the average value [31]. When an individual becomes contagious, a random decision is made for each contact whether an infection event occurs. Such an event can either lead to the infection of a susceptible individual or booster existing immunity. The time point of infection is sampled at random from the infectious individual's period of infectiousness, assuming a uniform distribution. Additional to the transmission of infection between individuals, all individuals (irrespective of their age or immunity status) are constantly exposed to a very low infection rate from the “outside” of the population. Transmission of the four influenza variants occurs independently, thus, on very rare occasions, individuals may be infected with more than one variant at the same time.
The basic reproduction number of the model is calibrated such that the median percentage of young adults (16 to 60 years of age) who are infected during the simulation year 2006/07 closely matches the observed value of 10.6% of German health care workers and controls who had a serologically confirmed influenza infection or at least a fourfold rise in influenza antibody titer during the 2006/07 season [45]. The chosen baseline value of the basic reproduction number of 1.575 yields a median infection rate of 10.6% (calculated from 2,000 simulations).
Appendix 4 Acquisition and loss of immunity
Immunity is assumed to be lost at a constant rate, leading to an exponentially distributed duration of immunity. The average duration of immunity is shorter if the immunity is caused by vaccination than when it is caused by infection. Immunity can be boostered by another infection or by a successful vaccination. If a booster event occurs, a new exponentially distributed duration of immunity is calculated for the individual. If the new duration is longer than the existing immunity, the individual’s immunity is extended accordingly. Although we use the same average length of immunity when calculating the new “booster immunity” duration, using the maximum of the existing and the new immunity duration automatically leads to a larger average duration of immunity after a booster infection than after a nonbooster infection. If an individual is infected or successfully vaccinated with a virus of one of the B lineages, a crossimmunizing event can occur which protects the individual also against the other B lineage (this is determined by calculating a random number). Such crossimmunizing events can immunize susceptible individuals or booster existing immunity. In most of our analyses, we keep the percentage of B lineage cross protection constant, but in a sensitivity analysis, we allow it to decrease linearly over time, starting with a value of 60% in 2001 when B/Vic is introduced and reaching the final value of e.g. 30% on the last simulation day.
A single infection with Influenza A(H1N1) is assumed to result in an immunity which lasts on average for 6 years [31]. As we assume that new A(H1N1) drift variants appear on average every 7 years and that 60% of the old immunity is active against new drift variants (see Appendix 6 for details), immunity loss due to drift alone would, therefore, result in an average immunity duration of 7 years/(1–0.6) = 17.5 years. In order to obtain an average immunity duration of 6 years, we assume that immunity is constantly lost at a rate λ while no drift variant is introduced. The average immunity duration D _{ I } is then given by D _{ I } = ((1–0.6)/7 + λ) ^{−1}. The value λ = 1/(9.13 years) which yields an average duration of 6 years is used for all four influenza strains (due to the higher drift frequency, the average immunity duration against A(H3N2) results in about 4.5 years).
The average duration of immunity after vaccination is calculated from a figure given by Ruf et al.[44] in which the seroprotection rate (SPR) (i.e. the fraction of vaccines with a heamagglutination inhibition titer of ≥ 40) is shown for various time periods after vaccination. Values of the percentage of vaccinees who were protected t months after vaccination were read from a figure using Aotsuka’s freeware package “Data Picker” M: (Shareware Version 1.2. http://hp.vector.co.jp/authors/VA019223/). We obtain the parameters of the function p _{ 0 } + (100p _{ 0 } ) _{ * } VE _{ * } e ^{αt/12} by least squares fit to the picked values. The SPR of vaccinees who were immune before vaccination (p _{ 0 }) is estimated to be 28.6% (95% CI from 23.5 to 33.7%), the fraction of individuals who were unprotected before vaccination, yet become protected by vaccination (VE) is estimated as 92.2% (95% CI from 84.9 to 99.6%), the immunity loss rate α is estimated to be 0.553 per year (95% CI from 0.394 to 0.729). The average duration of vaccination derived immunity (in the absence of drift) is given by 1/α = 1.81 years (95% CI from 1.37 to 2.54 years).
Each newborn individual inherits maternal antibodies against all influenza strains or lineages against which his or her mother is immune at the time of birth; maternal antibodies protect the newborn individual against infection for some months. The immunological and epidemiological statuses of immigrants also need to be initialized: this is done by choosing a random individual of similar age (preferably in the year cohort of the immigrant) and transferring his or her immunity status and infection status to the immigrant. In a sensitivity analysis, we additionally consider the possibility that infection causes a shortterm unspecific protection which temporarily prevents infection with any one of the influenza viruses.
Appendix 5 Vaccination coverage
The vaccination coverage depends on the age and the risk status of individuals. Furthermore, people who were vaccinated in the previous year are more likely to be vaccinated again. The baseline revaccination factor f _{ rev } of 2 means that previously vaccinated individuals are twice as likely to be chosen for vaccination; it can be interpreted as the “relative risk” (RR) of revaccination given a previous vaccination. At the beginning of each simulation year, each individual is assigned a “vaccination day”, irrespective of whether this individual will be vaccinated or not. The vaccination decision will be made on that day. If a percentage c _{ a,r } of individuals of age a and risk category r is intended to be vaccinated in a simulation year, we first calculate the expected number of vaccinees c _{ a,r } N _{ a,r }, whereby N _{ a,r } is the number of individuals in the population with age a and risk status r. If there are N ^{+} _{ a,r } previously vaccinated and N ^{−} _{ a,r } previously unvaccinated individuals, we sample (c _{ a,r } N _{ a,r } ) _{ * } (f _{ rev } N ^{+} _{ a,r } )/(N ^{−} _{ a,r } + f _{ rev } N ^{+} _{ a,r } ) previously vaccinated and (c _{ a,r } N _{ a,r } ) _{ * } N ^{−} _{ a,r } /(N ^{−} _{ a,r } + f _{ rev } N ^{+} _{ a,r } ) previously unvaccinated individuals to be vaccinated. During the vaccination period, it is continuously memorized how many individuals of each combination of age group, risk status and previous vaccination status have already been vaccinated and how many of them have already made the “decision” not to be vaccinated. Due to demographic changes (births, ageing, deaths, and immigrations) and due to individual changes in the risk status, individuals may change their age group or risk status during the vaccination period. If this happens, the memorized vaccination statistics is changed accordingly, so that the exact numbers of yet vaccinated and of yet unvaccinated individuals are known at any time during the vaccination period.
For an approximate calculation of the OR of revaccination, we neglect deaths and immigrations and assume that the vaccination coverage in the two subsequent years is the same (i.e. N ^{+} _{ a,r } = c _{ a,r * } N _{ a,r }, and N ^{−} _{ a,r } = (1c _{ a,r } ) _{ * } N _{ a,r } ), the relationship between the revaccination factor and the odds ratio (OR) of being vaccinated again can be approximately calculated as follows:
whereby z = N ^{2} _{ a,r } /(N ^{−} _{ a,r } + f _{ rev* } N ^{+} _{ a,r } ). Using the standard parameter value of f _{ rev } = 2, we get OR = 2/(1c _{ a,r } ). For a vaccination coverage of 20%, the OR is 2.5, for 60%, it is 5.0.
Appendix 6 Drift variants, vaccine mismatch, and cross protection
Smith et al.[35] constructed an antigenic map of Influenza A(H3N3) data from 35 years of surveillance (1968 to 2003) and found 11 antigenically distant clusters; the average antigenic distance between the centers of consecutive clusters was 4.5 units (SD 1.3). As influenza vaccines are updated when there is an antigenic difference of at least 2.0, these clusters can be regarded as representing genetically relevant drift variants. As ten new drift variants occurred in the time period of 35 years, we get an average duration of D = 3.5 years for new A(H3N2) drift variants (95% CI from 2.2 to 6.8). No direct estimates for the duration between consecutive Influenza B drift variants are available, but as they change at a rate which is about two to three times lower than this [3], we assume the average sojourn time of Influenza B drift variants to be 7 years. Judging by the frequency of changes in vaccine composition from 1987 to 2000, the variability of A(H1N1) may be similar to that of Influenza B [3]; we, therefore, use the same average duration of 7 years for Influenza A(H1N1) drift variants. For each one of the two A subtypes and the two B lineages, an independent random number r is calculated at the beginning of each simulation year to determine whether the circulating variant is replaced by a new drift variant. If r < 1/D, a new drift variant is introduced.
The possibility of a vaccine design mismatch is only considered for years in which new drift variants appear. Of course, TIV vaccination can only be affected by an Influenza B mismatch, if the mismatched B lineage is actually contained in the vaccine. For the last 18 years, four mismatch events have been reported (A(H3N2) in 2003 and 2004, and A(H1N1) in 2006 and 2009 [42]). For this time period, we calculate that 18/3.5 = 5.14 new drift variants for A(H3N2) and 18/7 = 2.57 new drift variants for A(H1N1) appeared. For the B lineage which was included in TIV during the 18 years, we also calculate 18/7 = 2.57 new drift variants. As we expect altogether 10.28 drift events and observed 4 mismatch events, we calculate a mismatch probability of 40% per drift. We use this probability for all four influenza variants, including the two B lineages for which no vaccine design mismatch has been reported.
The immunity shared between a previously circulating variant and a new variant is not well known. We infer information on this shared immunity from vaccination studies summarized by Jefferson et al.[39]: wellmatched TIV vaccination protects on average 73% of healthy adults, whereas mismatched TIV vaccination protects only 44%. We take the ratio of 44/73 = 60% also as indicator for the shared immunity between subsequent drift variants. 60% of the individuals who were immune against the previously circulating variant are also immune against the new drift variant, whereas the remaining 40% are susceptible to the new drift variant. For the sake of simplicity, we replace variants by new drift variants at the beginning of a simulation year (i.e. on July 1st). We do not allow cocirculation of the old and the new variant in the simulation, and we assume that a previously circulating variant will never be introduced again.
We furthermore use this estimate of 60% for shared immunity between the two B lineages: if an individual is infected with one B lineage, an additional immunization (or booster event) with respect to the other lineage takes place with a probability of 60%. Likewise, TIV vaccination can cause immunity (or booster events) against the B lineage which is not contained in the vaccine; the agedependent vaccine efficacy is multiplied by 0.6 to obtain the efficacy against the missing lineage. We restrict this crossimmunizing vaccination effect to TIV because there is no convincing evidence that increasing the antigen amount enhances the efficacy of influenza vaccines [57–59]. Tricco et. al.[10] report for vaccinees from 18 to 49 years of age that the vaccine efficacy against Influenza B is 77% if the lineage is contained in TIV and 52% if it is not contained in TIV, respectively. The ratio 52/77 = 0.675 indicates a B lineage cross protection of 67.5%. Langley et al. examine in a vaccination study with children what percentages seroconvert with respect to the B lineage which is contained and which is not contained in TIV. As they use two different TIV vaccines which either contain B/Vic or B/Yam, they obtain four different seroconversion rates with respect to the Influenza B lineages [59] (summarized in Table 3). We have developed a maximum likelihood model based on the assumption that either one of the two TIV vaccines causes a direct seroconversion S with respect to the B lineage contained in the vaccine and an indirect effect S x with respect to the lineage not contained in the vaccine. The joint likelihood L of the observations is then given by $\mathit{L}=\left(\begin{array}{c}\hfill 870\hfill \\ \hfill 359\hfill \end{array}\right){\left(\mathit{Sx}\right)}^{359}{\left(1\mathit{Sx}\right)}^{511}\left(\begin{array}{l}870\\ 622\end{array}\right){\mathit{S}}^{622}{\left(1\mathit{S}\right)}^{248}\left(\begin{array}{l}877\\ 644\end{array}\right){\mathit{S}}^{644}{\left(1\mathit{S}\right)}^{233}\left(\begin{array}{l}877\\ 262\end{array}\right){\left(\mathit{Sx}\right)}^{262}{\left(1\mathit{Sx}\right)}^{615}$, omitting the constant factors, this simplifies to L ∝ (Sx)^{359 + 262}(1 − Sx)^{511 + 615} S ^{622 + 644}(1 − S)^{248 + 233}.
The maximum likelihood estimates are S = 72.5% (95% confidence interval CI: 70.374.5, based on the profile likelihood) and x = 49.1% (CI: 45.752.5), respectively, i.e. the seroconversion rate of the B lineage not contained in TIV is about half as high as that of the lineage which is contained in the vaccine. Having obtained three indicators for the amount of B lineage cross protection (49.1% Langley, 60% Jefferson and 67.5% Tricco, respectively), we have decided to use the value 0.6 as the baseline value in our simulation studies.
Appendix 7 Results of TIV and QIV vaccination
During the 20 years of each simulation’s evaluation period, about half a million TIV or QIV vaccinations are performed in a population of 80,000 individuals. Due to the efficacy of the vaccine, about 38 percent of all vaccinations remain without immunologic effect (Figure 9). As many individuals are already immune when they are vaccinated, nearly a third of the successful vaccinations can only booster existing immunity. On average, 35% of all vaccinees newly acquire immunity against Influenza A. As TIV only contains one B lineage, the effect on Influenza B would only be half of that if no crossimmunizing events occurred. Due to crossimmunizing vaccination effects (light green and light blue areas), the effect of TIV vaccination is less pessimistic (Figure 9, middle). Using QIV, the failure rate drops to 38% (conveyed by an overall vaccine efficacy of 62%).
Appendix 8 Sensitivity analyses
In the following, we present the results of sensitivity analyses: the number of infections with and without vaccination is first accumulated for the 20 year evaluation period of each simulation and divided by 20 to obtain mean numbers of infections per simulation year. These means are used to calculate overall averages and 95% confidence intervals (given in Tables 4, 5, 6, 7 and 8). Pairwise differences between the TIV scenario and the QIV scenario of each simulation allow calculating the number of influenza infections prevented by TIV. The values in the tables are arithmetic means (averaged over 1,000 or 2,000 simulations) and their 95% confidence intervals. Please note that the percentages reported in Tables 4, 5, 6, 7 and 8 are calculated as arithmetic means of the percentages observed in the individual simulations. Calculating the percentages from the mean numbers of infections yield slightly different results: for the baseline parameter setting (R0 = 1.575 in Table 4), the calculation based on the mean values would yield 100%*395/(5597 + 3346) = 4.4% of infections prevented by QIV whereas the average percentage (calculated from 2,000 simulations) yields 4.3%.
Abbreviations
 4Flu:

Name of the individualbased seasonal influenza simulation tool
 TIV:

Trivalent influenza vaccine
 QIV:

Quadrivalent influenza vaccine
 CI:

Confidence interval
 OR:

Odds ration
 R0:

Basic reproduction number
 RR:

Relative risk, WHO, World health organization
 A(H1N1):

Influenza A subtype H1N1
 A(H3N2):

Influenza A subtype H3N2
 B/Yam:

Influenza B, Yamagata lineage
 B/Vic:

Influenza B, Victoria lineage.
References
 1.
Rota PA, Wallis TR, Harmon MW, Rota JS, Kendal AP, Nerome K: Cocirculation of two distinct evolutionary lineages of influenza type B virus since 1983. Virology. 1990, 175: 5968.
 2.
Kanegae Y, Sugita S, Endo A, Ishida M, Senya S, Osako K, Nerome K, Oya A: Evolutionary pattern of the hemagglutinin gene of influenza B viruses isolated in Japan: cocirculating lineages in the same epidemic season. J Virol. 1990, 64: 28602865.
 3.
Hay AJ, Gregory V, Douglas AR, Lin YP: The evolution of human influenza viruses. Philos Trans R Soc Lond B Biol Sci. 2001, 356: 18611870.
 4.
Chen R, Holmes EC: The evolutionary dynamics of human influenza B virus. J Mol Evol. 2008, 66: 655663.
 5.
McCullers JA, Saito T, Iverson AR: Multiple genotypes of influenza B virus circulated between 1979 and 2003. J Virol. 2004, 78: 1281712828.
 6.
Ambrose CS, Levin MJ: The rationale for quadrivalent influenza vaccines. Hum Vaccin Immunother. 2012, 8: 8188.
 7.
Thompson WW, Shay DK, Weintraub E, Brammer L, Bridges CB, Cox NJ, Fukuda K: Influenzaassociated hospitalizations in the United States. JAMA. 2004, 292: 13331340.
 8.
Irving SA, Patel DC, Kieke BA, Donahue JG, Vandermause MF, Shay DK, Belongia EA: Comparison of clinical features and outcomes of medically attended influenza A and influenza B in a defined population over four seasons: 2004–2005 through 2007–2008. Influenza Other Respi Viruses. 2012, 6: 3743.
 9.
World Health Organization: Vaccines against influenza WHO position paper  November 2012. Wkly Epidemiol Rec. 2012, 87: 461476.
 10.
Tricco AC, Chit A, Soobiah C, Hallett D, Meier G, Chen MH, Tashkandi M, Bauch CT, Loeb M: Comparing influenza vaccine efficacy against mismatched and matched strains: a systematic review and metaanalysis. BMC Med. 2013, 11: 153
 11.
Centers for Disease Control and Prevention: Past weekly surveillance reports. 2013, [http://www.cdc.gov/flu/weekly/pastreports.htm]
 12.
World Health Organization: Recommended composition of influenza virus vaccines for use in the 2013 southern hemisphere influenza season. 2012, [http://www.who.int/influenza/vaccines/virus/recommendations/2013_south/en/index.html]
 13.
McDonald SA, van Lier A, Plass D, Kretzschmar ME: The impact of demographic change on the estimated future burden of infectious diseases: examples from hepatitis B and seasonal influenza in the Netherlands. BMC Public Health. 2012, 12: 1046
 14.
Alexander ME, Kobes R: Effects of vaccination and population structure on influenza epidemic spread in the presence of two circulating strains. BMC Public Health. 2011, 11 (Suppl 1): S8
 15.
Edlund S, Kaufman J, Lessler J, Douglas J, Bromberg M, Kaufman Z, Bassal R, Chodick G, Marom R, Shalev V, Mesika Y, Ram R, Leventhal A: Comparing three basic models for seasonal influenza. Epidemics. 2011, 3: 135142.
 16.
Pitman RJ, Nagy LD, Sculpher MJ: Costeffectiveness of childhood influenza vaccination in England and Wales: results from a dynamic transmission model. Vaccine. 2012, 31: 927942.
 17.
Eichner M, Schwehm M, Duerr HP, Brockmann SO: The influenza pandemic preparedness planning tool InfluSim. BMC Infect Dis. 2007, 7: 17
 18.
Chao DL, Halloran ME, Obenchain VJ, Longini IM: FluTE, a publicly available stochastic influenza epidemic simulation model. PLoS Comput Biol. 2010, 6: e1000656
 19.
Halloran ME, Ferguson NM, Eubank S, Longini IM, Cummings DA, Lewis B, Xu S, Fraser C, Vullikanti A, Germann TC, Wagener D, Beckman R, Kadau K, Barrett C, Macken CA, Burke DS, Cooley P: Modeling targeted layered containment of an influenza pandemic in the United States. Proc Natl Acad Sci U S A. 2008, 105: 46394644.
 20.
Kenah E, Chao DL, Matrajt L, Halloran ME, Longini IM: The global transmission and control of influenza. PLoS One. 2011, 6: e19515
 21.
Colizza V, Barrat A, Barthelemy M, Valleron AJ, Vespignani A: Modeling the worldwide spread of pandemic influenza: baseline case and containment interventions. PLoS Med. 2007, 4: e13
 22.
Shaman J, Karspeck A: Forecasting seasonal outbreaks of influenza. Proc Natl Acad Sci U S A. 2012, 109: 2042520430.
 23.
Nishiura H: Realtime forecasting of an epidemic using a discrete time stochastic model: a case study of pandemic influenza (H1N12009). Biomed Eng Online. 2011, 10: 15
 24.
Tizzoni M, Bajardi P, Poletto C, Ramasco JJ, Balcan D, Goncalves B, Perra N, Colizza V, Vespignani A: Realtime numerical forecast of global epidemic spreading: case study of 2009 A/H1N1pdm. BMC Med. 2012, 10: 165
 25.
Elveback LR, Fox JP, Ackerman E, Langworthy A, Boyd M, Gatewood L: An influenza simulation model for immunization studies. Am J Epidemiol. 1976, 103: 152165.
 26.
Minayev P, Ferguson N: Improving the realism of deterministic multistrain models: implications for modelling influenza A. J R Soc Interface. 2009, 6: 509518.
 27.
Deutsches Statistisches Bundesamt: 2009, [https://www.destatis.de/DE/Startseite.html]
 28.
Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, Massari M, Salmaso S, Tomba GS, Wallinga J, Heijne J, SadkowskaTodys M, Rosinka M, Edmunds WJ: Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008, 5: e74
 29.
Biere B, Schweiger B: Molecular analyses of human influenza viruses. Circulation of new variants since 1995/96. Bundesgesundheitsblatt, Gesundheitsforschung, Gesundheitsschutz. 2008, 51: 10501060.
 30.
Britton T, Lindenstrand D: Inhomogeneous epidemics on weighted networks. Math Biosci. 2012, 24: 124131.
 31.
Vynnycky E, Pitman R, Siddiqui R, Gay N, Edmunds WJ: Estimating the impact of childhood influenza vaccination programmes in England and Wales. Vaccine. 2008, 26: 53215330.
 32.
Haas W: Influenza: Prävention, Diagnostik, Therapie und öffentliche Gesundheit. 2009, Munich, Germany: Elsevier GmbH
 33.
Zaman K, Roy E, Arifeen SE, Rahman M, Raqib R, Wilson E, Omer SB, Shahid NS, Breiman RF, Steinhoff MC: Effectiveness of maternal influenza immunization in mothers and infants. N Eng J Med. 2008, 359: 15551564.
 34.
Poehling KA, Szilagyi PG, Staat MA, Snively BM, Payne DC, Bridges CB, Chu SY, Light LS, Prill MM, Finelli L, Griffin MR, Edwards KM, New Vaccine Surveillance Network: Impact of maternal immunization on influenza hospitalizations in infants. Am J Obstetr Gynecol. 2011, 204: S141S148.
 35.
Smith DJ, Lapedes AS, de Jong JC, Bestebroer TM, Rimmelzwaan GF, Osterhaus AD, Fouchier RA: Mapping the antigenic and genetic evolution of influenza virus. Science. 2004, 305: 371376.
 36.
Park AW, Daly JM, Lewis NS, Smith DJ, Wood JL, Grenfell BT: Quantifying the impact of immune escape on transmission dynamics of influenza. Science. 2009, 326: 726728.
 37.
Arbeitsgemeinschaft für Influenza: Bericht zur Epidemiologie der Influenza in Deutschland. 2012, Saison 2011/12, [http://www.rki.de/DE/Content/InfAZ/I/Influenza/PK_AGI_2012_AGISaisonbericht_11_12.html]
 38.
Blank PR, Schwenkglenks M, Szucs TD: Disparities in influenza vaccination coverage rates by target group in five European countries: trends over seven consecutive seasons. Infection. 2009, 37: 390400.
 39.
Jefferson T, Di Pietrantonj C, Rivetti A, Bawazeer GA, AlAnsary LA, Ferroni E: Vaccines for preventing influenza in healthy adults. Cochrane Database Syst Rev. 2010, 7: CD001269
 40.
Jefferson T, Rivetti A, Di Pietrantonj C, Demicheli V, Ferroni E: Vaccines for preventing influenza in healthy children. Cochrane Database Syst Rev. 2012, 8: CD004879
 41.
Jefferson T, Di Pietrantonj C, AlAnsary LA, Ferroni E, Thorning S, Thomas RE: Vaccines for preventing influenza in the elderly. Cochrane Database Syst Rev. 2010, 7: CD004876
 42.
Arbeitsgemeinschaft für Influenza: Saisonberichte. 2013, [http://influenza.rki.de/Saisonbericht.aspx]
 43.
Robert Koch Institute: Influenzapandemieplan. [http://www.rki.de/DE/Content/InfAZ/I/Influenza/Influenzapandemieplan.html]
 44.
Ruf B, Colberg K, Frick M, Preusche A: Immunogenicity and antibody persistence of FluarixTM vs. Fluad(R) vs. Inflexal VTM in the elderly. 2004, Lisbon, Portugal: Poster presented during the First International Conference on Influenza Vaccines for the World
 45.
Williams CJ, Schweiger B, Diner G, Gerlach F, Haaman F, Krause G, Nienhaus A, Buchholz A: Seasonal influenza risk in hospital healthcare workers is more strongly associated with household than occupational exposures: results from a prospective cohort study in Berlin, Germany, 2006/07. BMC Infect Dis. 2010, 10: 8
 46.
Chowell G, Miller MA, Viboud C: Seasonal influenza in the United States, France, and Australia: transmission and prospects for control. Epidemiol Infect. 2008, 136: 852864.
 47.
Smieszek T, Balmer M, Hattendorf J, Axhausen KW, Zinsstag J, Scholz RW: Reconstructing the 2003/2004 H3N2 influenza epidemic in Switzerland with a spatially explicit, individualbased model. BMC Infect Dis. 2011, 11: 115
 48.
Lunelli A, Rizzo C, Puzelli S, Bella A, Montomoli E, Rota MC, Donatelli I, Pugliese A: Understanding the dynamics of seasonal influenza in Italy: incidence, transmissibility and population susceptibility in a 9year period. Influenza Other Respi Viruses. 2013, 7: 286295.
 49.
Boelle PY, Ansart S, Cori A, Valleron AJ: Transmission parameters of the A/H1N1 (2009) influenza virus pandemic: a review. Influenza Other Respi Viruses. 2011, 5: 306316.
 50.
Yang Y, Sugimoto JD, Halloran ME, Basta NE, Chao DL, Matrajt L, Potter G, Kenah E, Longini IM: The transmissibility and control of pandemic influenza A (H1N1) virus. Science. 2009, 326: 729733.
 51.
Tan X, Yuan L, Zhou J, Zheng Y, Yang F: Modeling the initial transmission dynamics of influenza A H1N1 in Guangdong Province, China. IJID Official Public Int Soc Infect Dis. 2013, 17: e479e484.
 52.
Marcelin G, Sandbulte MR, Webby RJ: Contribution of antibody production against neuraminidase to the protection afforded by influenza vaccines. Rev Med Virol. 2012, 22: 267279.
 53.
Truscott J, Fraser C, Cauchemez S, Meeyai A, Hinsley W, Donnelly CA, Ghani A, Ferguson N: Essential epidemiological mechanisms underpinning the transmission dynamics of seasonal influenza. J R Soc Interface. 2012, 9: 304312.
 54.
Dietz K: The estimation of the basic reproduction number for infectious diseases. Stat Methods Med Res. 1993, 2: 2341.
 55.
Goldstein E, Paur K, Fraser C, Kenah E, Wallinga J, Lipsitch M: Reproductive numbers, epidemic spread and control in a community of households. Math Biosci. 2009, 221: 1125.
 56.
Roberts MG: Epidemic models with uncertainty in the reproduction number. J Mathematic Biol. 2013, 66: 14631474.
 57.
Centers for Disease Control and Prevention: Prevention and control of influenza with vaccines: recommendations of the Advisory Committee on Immunization Practices (ACIP)United States, 2012–13 influenza season. MMWR. 2012, 61: 613618.
 58.
Centers for Disease Control and Prevention: Prevention and control of influenza with vaccines: recommendations of the Advisory Committee on Immunization Practices (ACIP), 2010. MMWR Recommend Rep/Centers Dis Control. 2010, 59: 162.
 59.
Langley JM, Carmona Martinez A, Chatterjee A, Halperin SA, McNeil S, Reisinger KS, Aggarwal N, Huang LM, Peng CT, GarciaSicilia J, Salamanca de la Cueva I, Cabanas F, TrevinoGarza C, RodriguezWeber MA, dela OM, Chandrasekaran V, Dewé W, Liu A, Innis BL, Jain VK: Immunogenicity and safety of an inactivated quadrivalent influenza vaccine candidate: a phase III randomized controlled trial in children. J Infect Dis. 2013, 208: 544553.
 60.
Carrat F, Vergu E, Ferguson NM, Lemaitre M, Cauchemez S, Leach S, Valleron AJ: Time lines of infection and disease in human influenza: a review of volunteer challenge studies. Am J Epidemiol. 2008, 167: 775785.
 61.
Claas EC, de Jong JC, Bartelds AI, Bijlsma K, Rothbarth P, de Groot R, Rimmelzwaan GF, Osterhaus AD: Influenza types and patient population. Lancet. 1995, 346: 180
 62.
Eames KT, Tilston NL, BrooksPollock E, Edmunds WJ: Measured dynamic social contact patterns explain the spread of H1N1v influenza. PLoS Comput Biol. 2012, 8: e1002425
 63.
Towers S, Chowell G: Impact of weekday social contact patterns on the modeling of influenza transmission, and determination of the influenza latent period. J Theor Biol. 2012, 312C: 8795.
 64.
Martin G, Marinescu MC, Singh DE, Carretero J: Leveraging social networks for understanding the evolution of epidemics. BMC Syst Biol. 2011, 5 (Suppl 3): S14
 65.
Salathe M, Kazandjieva M, Lee JW, Levis P, Feldman MW, Jones JH: A highresolution human contact network for infectious disease transmission. Proc Natl Acad Sci U S A. 2010, 107: 2202022025.
 66.
Cauchemez S, Bhattarai A, Marchbanks TL, Fagan RP, Ostroff S, Ferguson NM, Swerdlow D, Pennsylvania H1N1 working group: Role of social networks in shaping disease transmission during a community outbreak of 2009 H1N1 pandemic influenza. Proc Natl Acad Sci U S A. 2011, 108: 28252830.
 67.
Glass LM, Glass RJ: Social contact networks for the spread of pandemic influenza in children and teenagers. BMC Public Health. 2008, 8: 61
 68.
Riens B, Mangiapane S, Erhart M, Stillfried D: Analyse regionaler Unterschiede der InfluenzaImpfraten in der Impfsaison 2007/2008. Versorgungsatlas 2011. [http://www.versorgungsatlas.de/themen/versorgungsprozesse/?tab=6&uid=2]
 69.
Tacken MA, Jansen B, Mulder J, Visscher S, Heijnen ML, Campbell SM, Braspenning JC: Pandemic influenza A(H1N1)pdm09 improves vaccination routine in subsequent years: a cohort study from 2009 to 2011. Vaccine. 2013, 31: 900905.
 70.
Seo SH, Webster RG: Tumor necrosis factor alpha exerts powerful antiinfluenza virus Effects in Lung Epithelial Cells. J Virology. 2002, 76: 10711076.
 71.
Haller O, Kochs G, Weber F: The interferon response circuit: induction and suppression by pathogenic viruses. Virology. 2006, 344: 119130.
 72.
Gröndahl B, Ankermann T, von Bismarck P, Rockahr S, Kowalzik F, Gehring S, Meyer C, Knuf M, Puppe W: The 2009 pandemic influenza A(H1N1) coincides with changes in the epidemiology of other viral pathogens causing acute respiratory tract infections in children. Infection. 2014, 42: 303308.
 73.
Goodman AG, Zeng H, Proll SC, Peng X, Cilloniz C, Carter VS, Korth MJ, Tumpey TM, Katze MG: The alpha/beta interferon receptor provides protection against influenza virus replication but is dispensable for inflammatory response signaling. J Virology. 2010, 84: 20272037.
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712334/14/365/prepub
Acknowledgements
The authors would like to thank A. Anastassopoulou, G. Meier and E. Thommes, employed by the GlaxoSmithKline group of company, for their contribution in providing advice and helping with the parameterization and discussion of the model. The authors also thank Dr. Thomas Kochmann, CreativeAgents.net e. K, Wallenhorst, Germany, for valuable suggestions regarding the issue of indirect protection. Funding for this study was provided by GlaxoSmithKline GmbH & Co. KG, Munich, Germany. We acknowledge support by Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of Tübingen University.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
ME is shareholder of Epimos GmbH & Co. KG, which has received research support from the GlaxoSmithKline group of companies and AstraZeneca. MS is shareholder of ExploSYS GmbH, which has received payments from Epimos GmbH & Co. KG for developing the simulation tool 4Flu. JH and RSO are employees of the GlaxoSmithKline group of companies, and have restricted shares, stocks and/or stock options ownership in the GlaxoSmithKline group of companies. BS has received consultancy fees and honoraria for presentations from the GlaxoSmithKline group of companies and Novartis. MK has received consultancy fees and honoraria for advisory activities and presentations from Novartis, the GlaxoSmithKline group of companies and AstraZeneca.
Authors’ contributions
ME has created the transmission model and has written the major part of the manuscript. MS has designed and programmed the simulation tool and has taken part in the model development. JH and RSO participated to the structural design of the model and drafted parts of the introduction and the discussion section of the paper. HU, BS and MK have provided expert input on the epidemiology of seasonal influenza. They also have aided in the parameterisation and discussion of the model and have written parts of the manuscript. All authors had full access to the data and gave final approval before submission.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
12879_2013_3671_MOESM1_ESM.pdf
Authors’ original file for figure 1
12879_2013_3671_MOESM2_ESM.pdf
Authors’ original file for figure 2
12879_2013_3671_MOESM3_ESM.tiff
Authors’ original file for figure 3
12879_2013_3671_MOESM4_ESM.pdf
Authors’ original file for figure 4
12879_2013_3671_MOESM5_ESM.pdf
Authors’ original file for figure 5
12879_2013_3671_MOESM6_ESM.pdf
Authors’ original file for figure 6
12879_2013_3671_MOESM7_ESM.pdf
Authors’ original file for figure 7
12879_2013_3671_MOESM8_ESM.pdf
Authors’ original file for figure 8
12879_2013_3671_MOESM9_ESM.pdf
Authors’ original file for figure 9
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.
About this article
Cite this article
Eichner, M., Schwehm, M., Hain, J. et al. 4Flu  an individual based simulation tool to study the effects of quadrivalent vaccination on seasonal influenza in Germany. BMC Infect Dis 14, 365 (2014). https://doi.org/10.1186/1471233414365
Received:
Accepted:
Published:
Keywords
 Influenza
 Vaccination
 Simulation
 Mathematical model