Earliest infections predict the age distribution of seasonal influenza A cases ============================================================================== * Philip Arevalo * Huong Q. McLean * Edward A. Belongia * Sarah Cobey ## Abstract Seasonal variation in the age distribution of influenza A cases suggests that factors other than age shape susceptibility to medically attended infection. We ask whether these differences can be partly explained by protection conferred by childhood influenza infection, which has lasting impacts on immune responses to influenza and protection against new influenza A subtypes (phenomena known as original antigenic sin and immune imprinting). Fitting a statistical model to data from studies of influenza vaccine effectiveness (VE), we find that primary infection appears to reduce the risk of medically attended infection with that subtype throughout life. This effect is stronger for H1N1 compared to H3N2. Additionally, we find evidence that VE varies with both age and birth year, suggesting that VE is sensitive to early exposures. Our findings may improve estimates of age-specific risk and VE in similarly vaccinated populations and thus improve forecasting and vaccination strategies to combat seasonal influenza. ## Introduction Seasonal influenza is a serious public health concern, resulting in approximately 100,000-600,000 hospitalizations and 5000-27,000 deaths per year in the United States despite extensive annual vaccination campaigns (Reed et al., 2015). The rapid evolution of the virus to escape preexisting immunity contributes to the relatively high incidence of influenza, including in previously infected older children and adults. How susceptibility arises and changes over time in the host population has been difficult to quantify. A pathogen’s rate of antigenic evolution should affect the mean age of the hosts it infects, and differences in the rate of antigenic evolution have been proposed to explain differences in the age distributions of the two subtypes of influenza A. Compared to H3N2, H1N1 disproportionately infects children (Gagnon et al., 2018b; Caini et al., 2018; Khiabanian et al., 2009). It also evolves antigenically more slowly (Bedford et al., 2015). Thus, compared to H3N2, H1N1 is slower to escape immunity in individuals who have experienced prior infection (namely older children and adults), making them less susceptible to reinfection (Bedford et al., 2015; Beauté et al., 2015; Caini et al., 2018; Khiabanian et al., 2009). H3N2, in contrast, exhibits well known changes in antigenic phenotype that are expected to drive cases toward adults (Smith et al., 2004; Cobey and Koelle, 2008). Under this simple model, hosts previously infected with a subtype face equal risk of reinfection (on challenge) with an antigenic variant of that subtype. The age distributions of influenza cases in exceptional circumstances—pandemics and spillovers of avian influenza—have shown unexpected variation that suggests important effects of prior infection. Excess mortality in some adult cohorts during the 1918 and 2009 H1N1 pandemics correlates with childhood infection with other subtypes (Gagnon et al., 2013; Worobey et al., 2014; Gagnon et al., 2018a). In the post-2009 pandemic period, excess mortality and hospitalization were observed among cohorts first exposed to H2N2 or H3N2 during H1N1pdm-dominated seasons (Budd et al., 2019). Similarly, the subtypes circulating in childhood predict individuals’ susceptibility to severe zoonotic infections with avian H5N1 and H7N9, regardless of later exposure to other seasonal subtypes (Gostic et al., 2016). These patterns suggest that early influenza infections, and not prior infection per se, strongly shape susceptibility. Early infections might also affect the protection conferred by influenza vaccination. Foundational work on the theory of original antigenic sin demonstrated that an individual’s immune response to influenza vaccination is biased toward antigens similar to those encountered in childhood (Davenport and Hennessy, 1956, 1957). In some cases, this may result in a narrow antibody response focused on a single epitope (Davis et al., 2018). This phenomenon has been suggested to explain an unexpected decrease in vaccine effectiveness (VE) in the middle-aged in the 2015-2016 influenza season (Skowronski et al., 2017b; Flannery et al., 2018). More generally, it has been hypothesized that biases in immune memory can arise from both past infections and vaccinations and lead to variation in VE that is sensitive to the precise history of exposures (Smith et al., 1999; Skowronski et al., 2017a). To measure the effect of early exposures on medically attended infection risk and VE, we fitted statistical models to 3493 PCR-confirmed influenza cases identified through seasonal studies of influenza VE from the 2007-2008 to 2017-2018 seasons in the Marshfield Epidemiologic Study Area (MESA) in Marshfield, Wisconsin (Belongia et al., 2009, 2011; Griffin et al., 2011; Treanor et al., 2012; Ohmit et al., 2014; McLean et al., 2014; Gaglani et al., 2016; Zimmerman et al., 2016; Jackson et al., 2017; Flannery et al., 2018, Figure 1-Supplement 1). Each influenza season, individuals in a defined community cohort were recruited and tested for influenza when seeking outpatient care for acute respiratory infection. Eligibility was restricted to individuals >6 months of age living in MESA who received routine care from the Marshfield Clinic and who presented in an outpatient setting. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/03/29/19001875/F1) Figure 1. Summary of models tested. Ten different models result from considering different combinations of HA imprinting and VE. We also tested one additional model excluding the effects of N2 and HA imprinting (Materials and Methods: “Modeling approach”). We sought to explain the variation in the age distribution of these cases by subtype and over time. Our model predicted the relative number of cases of influenza in each birth year each season as a function of the age structure of the population, age-specific differences in the risk of medically attended influenza A infection, early influenza infection, and vaccination. Despite the extensive antigenic evolution in both subtypes over the study period, we found strong evidence of protection from the subtype to which a birth cohort was likely first infected (the imprinting subtype) and variation in VE by birth cohort. ## Materials and Methods ### Study cohort Cases of PCR-confirmed, medically attended influenza were identified from annual community cohorts based on residency in MESA. MESA is a contiguous geographic area surrounding Marshfield, Wisconsin, where nearly all 61,000 residents receive outpatient and inpatient care from the Marshfield Clinic Health System (Kieke et al., 2015). For each influenza season from 2007-2008 through 2017-2018, we identified MESA residents >6 months of age who received routine care from the Marshfield Clinic. These individuals were eligible for recruitment into that season’s VE study if they sought care for acute respiratory infection. Trained research coordinators recruited patients during clinical encounters in primary care departments, including urgent care, pediatrics, combined internal medicine and pediatrics, internal medicine, and family practice. Patients were enrolled on weekdays, evenings, and weekends when clinical services were provided. Research staff used an electronic appointment system to screen the chief complaints for respiratory or febrile illness. Patients were then approached in-person to assess eligibility based on specific respiratory symptoms and duration of illness. The proportion of patients with medically attended acute respiratory infection (MAARI) who were screened for enrollment varied by season and was largely determined by the volume of patients each day and staffing capacity. Only symptoms and illness duration were used to determine eligibility among those patients who were in the predefined cohort. Patients were also assessed for the presence of medical conditions that put them at high risk for complications from influenza infection, as defined by the Advisory Committee on Immunization Practice (Smith et al., 2006). These conditions included cardiovascular disease, diabetes, pulmonary disease, cancer, kidney disease, liver disease, blood disorders, immunosuppressive disorders, metabolic disorders, and neurological/musculoskeletal disorders. We considered subjects vaccinated if they received that season’s influenza vaccine >14 days before enrollment. For the 2009-2010 season, we only considered receipt of the 2009 monovalent vaccine. The Marshfield Clinic generally does not capture MAARI in nursing facilities with dedicated medical staff, causing undersampling of the oldest age groups. We adjusted for this (Appendix 1: “Age-specific rates of approachment, enrollment, and nursing home residence”). Each season, recruitment began when influenza activity was detected in the community and usually continued for 12-15 weeks. Symptom eligibility criteria varied by season but included fever/feverishness or cough during most seasons. We retroactively standardized symptom eligibility criteria to only require cough as a symptom. Individuals with illness duration >7 days or presenting in an inpatient (hospital) setting were excluded. After obtaining informed consent, a mid-turbinate swab was obtained for influenza detection. RT-PCR was performed using CDC primers and probes to identify influenza cases, including type and subtype. ### Calculating differences in the age distribution between seasons We defined the age distribution of each season as the number of cases of the dominant (more common) subtype in each of nine age groups (0-4 year-olds, 5-9 year-olds, 10-14 year-olds, 15-19 year-olds, 20-29 year-olds, 30-39 year-olds, 40-49 year-olds, 50-64 year-olds, and >65 years old). We excluded the subdominant subtype in each season due to concerns that short-term interference between the subtypes (Laurie et al., 2015; Goldstein et al., 2011) would affect the age distribution of the rarer subtype. The G-test of independence was used to measure differences in seasons’ age distributions. ### Calculating relative risk To evaluate relative infection risk in different age groups, we measured their relative risk of infection in the first versus second half of each season. This risk is a combination of the chance of infection, conditional on infection (susceptibility), and the rate of contact with infected people. Attack rates should be higher in populations that experience more risk, and therefore these populations should be infected earlier in the epidemic (Worby et al., 2015). To calculate relative risk we used an approach similar to Worby et al., 2015. We defined the midpoint of each season as the week in which the cumulative number of cases of the dominant subtype among all people exceeded half the total for that season. Weeks before and after this point were assigned to the first and second half of the season, respectively. We assigned each case to one of the five age groups used by Worby et al., 2015 (0-4 year-olds, 5-17 year-olds, 18-49 year-olds, 50-64 year olds, and >65 years old). For each age group *g*, we defined relative risk as ![Formula][1] where *C*first,*t*,*g* and *C*second,*t*,*g* are the fraction of cases of the dominant subtype during influenza season *t* that occurred during the first or second half of the season, respectively. A relative risk >1 indicates that cases in an age group were more likely to occur during the first half of the season. ### Calculating imprinting probabilities We hypothesized that the subtype of a person’s first influenza A infection affects their future susceptibility to that subtype. Testing this hypothesis requires knowing the probability that a person’s primary influenza A infection was with a particular subtype. To calculate these probabilities, we emulated the approach of Gostic et al., 2016, which assumes these probabilities are determined by a person’s year of birth and subsequent exposure to each subtype. First, we calculated the probability that an individual born in year *y* received their first influenza A exposure in influenza season *t*. Assuming a constant per-season rate of infection *i*, the probability of infection in one season (i.e., the attack rate) is given by ![Formula][2] By assuming that the average probability that a naive individual is infected in a single season is 0.28 (Bodewes et al., 2011; Gostic et al., 2016), we calculated the expected per-season infection rate (*i*) as ![Formula][3] However, because the intensity of epidemics varies between seasons (*I**t*, Appendix 1: “Seasonal intensity”) and the fraction of the epidemic experienced by a person depends on their birth year (*y**y*,*t*, Appendix 1: “Fraction of season experienced”), we considered the time-varying per-season infection rate, ![Formula][4] Therefore, the probability that a naive individual born in year *y* is infected in season *t* is ![Formula][5] We used *a**y*,*t* to calculate the fraction of a birth cohort *y* that received their first influenza A infection in season *t*. Let *U**y*,*t* be the fraction of people born in year *y* who were unexposed at the beginning of season *t* (Appendix 1: “Calculating the fraction unexposed”). The probability that a person born in year *y* has their first infection in season *t* is ![Formula][6] We calculated *m**s*,*t*,*y*, the probability that a person born in year *y* had their first influenza A infection with subtype *s* in season *t*, by multiplying *a**y*,*t**U**y*,*t* by the frequency of subtype *s* in season *t, l**s*,*t* (Figure 3-Supplement 1), ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F2.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/03/29/19001875/F2) Figure 3. Imprinting probabilities by age group across seasons. Each panel shows the imprinting probabilities of an age group from the 2007-2008 season through the 2017-2018 season. The color of each bar corresponds to the imprinting subtype or naive individuals, who have not yet been infected. ![Formula][7] ### Modeling approach We aimed to predict *p**s*,*t*,*y*,*v*, the fraction of cases of subtype *s* in season *t* among people born in year *y* with vaccination status *v*. Our models assume that this is proportional to a combination of the following factors: 1. *Demography*. The age distribution of our study cohort is not static over the study period. All models adjusted for the changing fractions of the population in each birth cohort and season (Figure 1-Supplement 2; Mathematical expressions for model components: “Demography”). 2. *Age-specific effects*. We considered that age itself may be associated with differences in medically attended influenza A infection risk stemming from differences in susceptibility and/or rates of contact with infectious people. Additionally, we expect that age groups may intrinsically differ in their healthcare-seeking behaviors. These factors are inseparable in our data, and all models represent their combined effects with a static age-specific parameter shared by both subtypes that describes the risk of age-specific medically attended influenza A infection (Mathematical expressions for model components: “Age-specific factors”). We assumed no intrinsic differences in the age-specific virulence of the two subtypes. These age-specific parameters were fitted. We also adjusted for other potential sources of age-specific bias, including age-specific differences in study approachment and enrollment rates (Appendix 1: “Age-specific rates of approachment, enrollment, and nursing home residence”). 3. *Imprinting*. We tested several hypotheses of how primary exposures could affect the risk of medically attended infection with H1N1 and H3N2. In each version, we estimated fractional reductions in risk of medically attended H1N1 and H3N2 infection due to primary (i.e., imprinting) exposure to the same type: * Subtype-specific imprinting: Influenza has two main antigens, hemagglutinin (HA) and neuraminadase (NA). Imprinting could in theory derive from responses to either or both antigens. Because H1N1 is the only seasonal subtype of influenza with N1, we cannot separate the effects of initial N1 exposure from initial H1 exposure. However, since N2 appears in both H3N2 and H2N2 viruses, we can estimate protection against H3N2 infection from initial N2 exposure separately from protection from initial H3 exposure (Mathematical expressions for model components: “HA subtype imprinting” and “N2 imprinting”). * Group-level imprinting: Influenza A viruses fall into two groups (I and II) corresponding to the two phylogenetic clades of HA. Gostic et al., 2016 found that primary infection by a virus belonging to one group protected against severe infection by another subtype in the same group. If group-level imprinting were influential, we would see primary infection with H2N2 conferring protection against H1N1, another group I virus, as well as H1N1 protecting against H1N1, and H3N2 against H3N2. We considered a separate class of models that assumes group-level protection instead of subtype-specific protection (Mathematical expressions for model components: “HA group imprinting”). 4. *Vaccination*. Approximately 45% of the MESA population was vaccinated against influenza each year (Figure 1-Supplement 3; Appendix 1: “Vaccination coverage”). We estimated cases in vaccinated and unvaccinated individuals of each birth year separately. Naively, we expect that vaccinated individuals should seek medical attention for acute respiratory infection proportionally to the fraction of their cohort vaccinated that season. However, vaccinated individuals may seek medical attention for acute respiratory infection more frequently than non-vaccinees due to correlations between the decision to vaccinate, healthcare-seeking behavior, and underlying medical conditions (Jackson et al., 2005a,b; Belongia et al., 2009). Indeed, we generally observed higher rates of high-risk medical conditions among vaccinated people compared to unvaccinated people (Figure 1-Supplement 4). We attempted to adjust for this by calculating the fraction of vaccinated people among those who had MAARI and tested negative for influenza (i.e., the test-negative controls, “Mathematical expressions for model components: Vaccination”). We found that the vaccinated fraction exceeds vaccination coverage for most age groups, suggesting vaccinated individuals are overrepresented among cases for reasons unrelated to influenza (Figure 1-Supplement 5). We also assumed that vaccination is not perfectly effective, and defined VE as the fractional reduction in cases expected in vaccinated compared to unvaccinated individuals after controlling for the effects described above. We estimated subtype-specific VE under five scenarios: (i) constant across age groups and seasons; (ii) constant across age groups but season-specific; (iii) age-specific but constant across seasons; (iv) imprinting-specific; and (v) birth-cohort-specific. We assumed that vaccination affects risk only in the current season, i.e., vaccination in a prior season confers no residual protection (Mathematical expressions for model components: “Vaccination”; Ohmit et al., 2014, 2015; Jackson et al., 2017; Skowronski et al., 2016; Pebody et al., 2013; McLean et al., 2018). We defined models as specific combinations of the above factors. We tested a set of 10 models by pairing each of the possible implementations of HA imprinting with each implementation of VE (Figure 1). Demography, age-specific effects, and N2 imprinting were included in all these models. To test whether more complex models truly improved model fit, we also tested a simple model with constant VE and no effect of imprinting. We evaluated these 11 models by maximum likelihood and compared their performance using the corrected Akaike information criterion (cAIC, “Model likelihood”) and leave-one-out cross-validation. ### Mathematical expressions for model components #### Demography We expect that the fraction of cases in each birth cohort should be proportional to the underlying demographic birth year distribution of the population. To calculate the demographic birth year distribution, we used MESA-specific data on the age distribution for each season (Kieke et al., 2015). Because people >90 years old were grouped into a single age class, we estimated the number of people in each age >90 years old by assuming a geometric decline in population with age. We converted the age distribution for each season into a distribution by birth year by assigning people of a specific age into the two possible birth years of that age (Appendix 1: “Birth year distribution of the study population”). Therefore, ![Formula][8] where *D**t*,*y* is the fraction of the population in season *t* who were born in year *y*. #### Age-specific factors We modeled intrinsically age-specific differences in medically attended influenza A infection risk and healthcare-seeking behavior by using parameters that represent the relative risk of medically attended influenza A infection in each age group. These parameters combine the effects of underlying age-specific differences in influenza A medically attended infection risk as well as age-specific differences in healthcare-seeking behavior. We considered the same age groups as before (0-4 year-olds, 5-9 year-olds, 10-14 year-olds, 15-19 year-olds, 20-29 year-olds, 30-39 year-olds, 40-49 year-olds, 50-64 year-olds, and >65 years old). We chose 20-29 year-olds as our reference age group. All age groups *g* aside from 20-29 year-olds had an associated parameter (*A**g*) that scaled their risk of medically attended influenza A infection relative to 20-29 year-olds. These parameters can take on any positive value. Since our models describe the distribution of cases by birth year and not by age, we mapped the age-group-specific parameters (*A**g*) to birth cohorts in each season *t* (*A**t*,*y*). We considered that each birth cohort has two possible ages in each season (*a*1 and *a*2). Let *G*(*a*) be a function that specifies the age group *g* of a given age *a*. Then *A**t*,*y*, the age-specific relative risk in season *t* of medically attended influenza A infection for a person born in year *y*, is ![Formula][9] where *f**a*1,*t*,*y* and *f**a*2,*t*,*y* are the fractions of birth cohort *y* who are age *a*1 or *a*2 in influenza season *t* (Appendix 1: “Fraction of birth cohort with specific age”), and *A**G*(*a*1) and *A**G*(*a*2) are the age-group-specific parameters for *a*1 and *a*2. Our models also included age-specific approachment rates ![Graphic][10], enrollment rates ![Graphic][11], and nursing home enrollment (*k**t*,*y*) as covariates, all of which bias the age distribution of medically attended influenza infections (Appendix 1: “Age-specific rates of approachment, enrollment, and nursing home residence”). The combination of estimated age-specific effects and age-specific covariates was modeled as ![Formula][12] #### HA subtype imprinting We considered that imprinting to HA reduces a birth cohort’s risk of future infection from the same HA subtype. Therefore, ![Formula][13] where *h**s* is the strength of HA imprinting for subtype *s* and *m**s*,*t*,*y* is the imprinting probability in season *t* of birth cohort *y* to subtype *s* (“Calculating imprinting probabilities”). #### HA group imprinting We considered that imprinting to HA reduces a birth cohort’s risk of future infection with viruses from the same HA group. Therefore, ![Formula][14] ![Formula][15] where *g*1 is the strength of HA imprinting for group 1 viruses; *g*2 is the strength of HA imprinting for group 2 viruses; and *m*H1N1,*t*,*y*, *m*H2N2,*t*,*y*, and *m*H3N2,*t*,*y* are the imprinting probabilities in season *t* of birth cohort *y* to H1N1, H2N2, and H3N2. #### N2 imprinting We considered that imprinting to N2 reduces a birth cohort’s risk of H3N2 infection. Therefore, ![Formula][16] where *n**m* is the strength of N2 imprinting, and *m*H3N2,*t*,*y* and *m*H2N2,*t*,*y* are the imprinting probabilities of birth cohort *y* in season *t* to H3N2 and H2N2. #### Vaccination We assumed that vaccination decreases the risk of medically attended infection. However, vaccinated individuals may seek healthcare for symptomatic influenza at a different rate than unvaccinated individuals. Moreover, because vaccines are routinely recommended for individuals with underlying health conditions, pre-existing susceptibility to MAARI among vaccinated individuals may also differ from unvaccinated individuals. Let *R**t*,*g* represent the fraction of vaccinated individuals in age group *g* in season *t* that present with MAARI. We use test-negative controls to estimate this as ![Formula][17] where ![Graphic][18] and ![Graphic][19] are the number of vaccinated or unvaccinated individuals born in year *g* presenting with MAARI and testing negative for influenza in season *t*. We converted *R**t*,*g* to *R**t*,*y* (i.e., to a covariate indexed by birth cohort) using the same method described in “Age-specific factors.” We tested five different VE schemes: subtype-specific VE that remained constant across seasons and cohorts (2 parameters), subtype-specific VE that varied between the age groups described above (18 parameters), VE that varied between seasons (12 parameters), VE for each possible imprinting subtype (6 parameters), and birth-cohort-specific VE (18 parameters). These VE parameters (*V*) reduced the probability of medically attended influenza A infection among vaccinated individuals in a birth cohort, i.e., ![Formula][20] ![Formula][21] where *V* depends on the specific implementation of VE used. Constant VE only varies with the infecting subtype, thus ![Formula][22] Season-specific VE varies with subtype and season, thus ![Formula][23] For age-specific VE, we used the same age classes described above for “Age-specific factors” but did not consider a reference age class, so that each age group had an associated VE for each subtype. We used these age-specific VE parameters to calculate the VE against subtype *s* in birth cohort *y* during season *t* using the same procedure described in “Age-specific factors” (Equation 9). Therefore, ![Formula][24] where *v**G*(*a*1),*s* and *v**G*(*a*2),*s* are age-specific VE parameters for *a*1 and *a*2. For imprinting-specific VE, we used the imprinting probabilities for each birth cohort described in “Calculating imprinting probabilities” to scale V such that ![Formula][25] where *v**s*,*z* is the VE among people imprinted to subtype *z* against infection by dominant subtype *s*, and *m**z*,*t*,*y* is the imprinting probability for subtype *z* in season *t* for birth cohort *y*. For birth-cohort-specific VE, we defined nine birth cohorts corresponding to the nine age groups we used for the 2017-2018 season: 1918-1952, 1953-1967, 1968-1977, 1978-1987, 1988-1997, 1998-2002, 2003-2007, 2008-2012, and 2013-2017. Let *Q*(*y*) be the birth cohort of people born in year *y*. Then ![Formula][26] where *v**Q*(*y*),*s* is the VE among people in cohort *Q*(*y*) against infection by dominant subtype *s*. ### Model likelihood Recall that our aim is to predict *p**s*,*t*,*y*,*v*, the fraction of all PCR-confirmed influenza cases of dominant subtype *s* in influenza season *t* among people born in year *y* with vaccination status *v*. These fractions can also be interpreted as multinomial parameters that describe the probability that in season *t*, a medically attended influenza infection of subtype *s* occurs among people born in year *y* with vaccination status *v*. Each model *M* assumes that *p**s*,*t*,*y*,*v* is proportional to a collection of model components *j* described above (demography, age, imprinting, and vaccination). Thus, ![Formula][27] where *p**M*,*s*,*t*,*y*,*v* is a multinomial probability under model *M, ϕ**M*,*j* indicates whether model *M* contains component *j*, and *η**j*,*s*,*t*,*y*,*v* is the mathematical expression for model component *j* given *s, t, y*, and *v* (e.g., for HA subtype imprinting, *η**j*,*s*,*t*,*y*,*v* = 1 − *h**s**m**s*,*t*,*y*). To obtain proper multinomial probabilities, we calculated a normalizing constant for each season *t* such that all probabilities in that season sum to 1. For convenience, let ![Graphic][28] be the unnormalized multinomial probability for model *M*. Then for a specific season *t*, the normalized multinomial probability is ![Formula][29] where *y*max,*t* is the maximum birth year possible for a specific season *t*. To calculate the likelihood of a given model, we used the multinomial probabilities and the observed birth year distribution of cases. Let *n**s*,*t*,*y*,*v* be the number of PCR-confirmed cases of dominant subtype *s* in influenza season *t* among people born in year *y* with vaccination status *v*. The total number of PCR-confirmed cases of dominant subtype *s* in season *t* is ![Formula][30] For models fitted to a restricted set of ages, we limited the cases for each season to the birth cohorts that were guaranteed to meet the age requirements in that season. Then, the likelihood of model *M* in season *t* is given by the multinomial likelihood, ![Formula][31] Finally, the full model likelihood for model *M* over all observed seasons is ![Formula][32] We fitted the model to case data using the L-BFGS-B algorithm implemented in the R package *optimx*. We estimated 95% confidence intervals for parameters of the best-fitting model by evaluating likelihood profiles at 14 evenly spaced points and interpolating the entire profile using a smoothing spline. ### Code and data availability The code and data used to perform the analyses for this project are available at [https://github.com/cobeylab/FluAImprinting](https://github.com/cobeylab/FluAImprinting). ## Results ### The age distribution of cases varies between seasons and subtypes The age distribution of cases varies between subtypes. The relative burden of cases is consistently higher in people >65 years old during H3N2-dominated seasons compared to H1N1-dominated seasons (Figure 2). The age distribution tends to vary more between subtypes than within either over time (Figure 2-Supplement 1, off-diagonal quadrants). This is consistent with recent work showing that the ratios of H3N2 to H1N1 cases differ between age groups (Gagnon et al., 2018b). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F3.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/03/29/19001875/F3) Figure 2. The age distribution of cases. **A**. The age distributions of cases from the 2007-2008 through the 2017-2018 influenza seasons in MESA. Dark lines with open circles indicate the average fraction of cases in each age group. Lighter-colored lines show the age distribution for individual seasons. **B**. The age distribution of cases in H1N1-dominated seasons. **C**. The age distribution of cases in H3N2-dominated seasons. The age distribution also varies within subtypes over time (Figure 2-Supplement 1, diagonal quadrants). The seven H3N2-dominated seasons display three types of age distributions (Figure 2-Supplement 1, clusters of lighter-colored cells in the upper left-hand quadrant), and two correspond to major antigenic clusters (2007-2008, Fonville et al., 2015; 2010-2012, Ann et al., 2012). These differences sometimes coincide with significant shifts in the age distribution between seasons. For instance, the highest fraction of H3N2 cases occurs in 20-29 year olds in the 2007-2008 season, but this age group has the lowest fraction of cases in the next H3N2-dominated season (2010-2011, Figure 2C). In H1N1, the shift from seasonal to pandemic strains is associated with large changes in the age distribution (Figure 2-Supplement 1, lower right-hand quadrant). We found further evidence that age groups differed in their susceptibility across seasons by examining the relative risk of infection during the first versus second half of each epidemic period (Materials and Methods: “Calculating relative risk”). Individuals at greater risk of infection should be infected disproportionately early rather than late in an epidemic (Worby et al., 2015). We confirmed that an age group’s relative risk correlates with the fraction of cases within that age group in the same season (Pearson’s r=0.58, 95% CI 0.38-0.73; Figure 2-Supplement 2A; Appendix 1: “Correlation of relative risk and fraction of cases”). This trend is evident for H1N1 (Pearson’s r=0.73, 95% CI 0.45-0.88; Figure 2-Supplement 2A) and H3N2 seasons separately (Pearson’s r=0.52, 95% CI 0.30-0.69; Figure 2-Supplement 2A). The positive correlation in all seasons is robust to undersampling of cases at the start and end of seasons (Appendix 1: “Sensitivity to sampling effort”, Figure 2-Supplement 2B). This provides supporting evidence that the different numbers of cases in each age group reflect underlying differences in infection risk. Just as the age distribution of cases varies over time, the age groups with high relative risks of infection change over time. If people contact one another similarly from one season to the next, these shifting relative risks imply that age groups’ relative susceptibilities change over time. For instance, 5-17 year olds had the highest relative risk of early infection in the 2008-2009 season, whereas 50-64 year-olds had the highest relative risk in the 2013-2014 season (Figure 2-Supplement 3). Relative risks in MESA vary more than national estimates, which show that 5-17 year-olds had the highest relative risk in all but one season from the 2009 pandemic to 2013-2014 (Worby et al., 2015). These differences may partly be due to the fact that our measurements of relative risk use outpatient visits, whereas the national estimates use hospitalizations. Taken together, these findings suggest that the risk of influenza infection is not a simple function of age alone. Other factors, such as past influenza infections and vaccination, might explain the changing age distributions of cases in time. ### Imprinting probabilities of age groups change over time We hypothesized that variation in the age distribution of cases could be explained by the aging of birth cohorts with similar early exposure histories. This would cause the early exposure history of an age group, and thus potentially its susceptibility, to change in time. To calculate the probability that people in a particular age group had their first influenza A infection with a particular subtype, we adapted the approach from Gostic et al., 2016. Briefly, we calculated the probability that an individual born in a specific year had a primary infection with H1N1, H2N2, or H3N2 using data on relative epidemic sizes and the frequencies of circulating subtypes (Figure 3-Supplement 1; Materials and Methods: “Calculating imprinting probabilities”). As expected, age groups’ early exposures are not static and change over time (Figure 3). Older people nonetheless tend to be imprinted to H1N1 or H2N2, whereas younger people have higher probabilities of imprinting to H3N2. The effects of the 2009 H1N1 pandemic are evident in the three youngest age groups as a transient increase (from 2009 to approximately 2013) in their H1N1 imprinting probability. These imprinting probabilities are relatively well-constrained even after for accounting for uncertainty in epidemic size (Figure 3-Supplement 2; Appendix 1: “Sensitivity to uncertainty in ILI and the frequency of influenza A”). ### Age-specific differences in medically attended influenza A infection risk affect epidemic patterns We fitted models to estimate the underlying effects of age, early infections, and vaccination on the age distributions of cases. As expected, the cases reveal age-specific differences in the risk of medically attended influenza A infection (Figure 4; Figure 4-Supplement 1; Appendix 2: Table 1). This risk is roughly threefold higher among children <4 years old compared to adults 20-29 years old, after adjusting for other effects (Figure 4). The decline in risk through middle age is generally consistent with attack rates estimated from serology (Monto et al., 1985; Bodewes et al., 2011; Wu et al., 2010; Huang et al., 2019) and clinical infections (Wu et al., 2017). We recently observed smaller differences in the attack rates of school-aged children and their parents when estimating infections serologically (Ranjeva et al., 2019). We hypothesize that the attack rates estimated from clinical infections might show larger differences by age due to age-related changes in infection severity and healthcare-seeking behavior. Indeed, rates of healthcare-seeking behavior have been shown to decline with age before rising in adults >65 years old (Biggerstaff et al., 2014; Brooks-Pollock et al., 2011; Van Cauteren et al., 2012), consistent with our results. Finally, the increased risk of medically attended influenza A infection among people >65 years old compared to other adults may be related to the increasing prevalence of high-risk medical conditions with age (Figure 1-Supplement 4). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/03/29/19001875/F4) Figure 4. Estimates of relative age-specific medically attended influenza infection risk. Open circles represent the maximum likelihood estimates of parameters describing age-specific differences in the relative risk of medically attended influenza A infection. Lines show the 95% confidence interval. ### Initial infection confers long-lasting, subtype-specific protection against future infection Our best-fitting model supports subtype-specific imprinting for H1N1 and H3N2 (Figure 5, top row; Appendix 2: Table 1). This model also provides the best predictive power compared to other models in a leave-one-out cross-validation analysis (Figure 5-Supplement 1; Figure 5-Supplement 2; Appendix 1: “Evaluation of predictive power”). The risk of future medically attended infection by H1N1 is reduced by 66% (95% CI 53-77%) in people imprinted to H1N1, whereas the risk of future medically attended infection by H3N2 is reduced by 33% (95% CI 17-46%) in people imprinted to H3N2. We found no evidence of a protective effect from imprinting to N2 (0%, 95% CI 0-7%). These estimates of imprinting protection are insensitive to ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2020/03/29/19001875/F5) Figure 5. Estimates of imprinting strength. Imprinting is more protective against H1N1 infection than H3N2 infection. Open circles represent the maximum likelihood estimates of imprinting parameters from the model including HA subtype imprinting and age-specific VE fitted to the indicated age group (y-axis). Black lines show 95% confidence intervals. * uncertainty in imprinting probabilities due to uncertainty in past epidemic sizes (Figure 3-Supplement 2; Appendix 1: “Sensitivity to uncertainty in ILI and the frequency of influenza A”; Appendix 2: Table 3), * choice of age groups for medically attended influenza A infection risk and VE (Appendix 1: “Sensitivity to age groups”; Appendix 2: Table 4), and * undersampling of influenza cases in some seasons (Figure 5-Supplement 3). In theory, the estimated protective effects of imprinting could be influenced by cross-protection rather than the impact of first infection per se. Because first infections are also recent infections in children, we reasoned that the observed imprinting effects might arise from confounding with recent infections in these ages. Based on an estimated 7-year half-life of homologous protection after H1N1pdm infection in children (Ranjeva et al., 2019) and the fact that most children experience primary influenza A infection by 5 years of age (Bodewes et al., 2011), we reasoned that excluding children <15 years old would diminish the impact of protection from recent infection on our results. When we excluded the youngest age groups, our estimates of H1N1 imprinting protection decreased while H3N2 imprinting protection increased (Figure 5, second row). However, initial infection by H1N1 was still more protective than initial infection by H3N2, both imprinting effects remained positive, and there was no significant change in the values of other estimated parameters (Appendix 2: Table 1 and Table 2). The effects of recent infection should also manifest in the difference between the observed and estimated numbers of cases (i.e., the excess cases, Appendix 1: “Calculating excess cases”), since unlike typical transmission models, our model does not take prior-season infections into account when estimating cases for the current season. More infections in a birth cohort in one season should reduce susceptibility in that birth cohort at the start of the next season. We thus expect that excess cases in one season will be followed by missing cases in the next season dominated by that subtype (i.e., a negative correlation in excess cases). Instead, we observed that excess cases for each birth cohort are weakly positively correlated from season to season, suggesting that immunity from recent infections is not a major driver of temporal variation in the age distribution of cases (Figure 5-Supplement 5). Since older adults have the highest probability of primary infection with H1N1, we also reasoned that older adults might disproportionately drive the strong protection from H1N1 imprinting we observe. People born before 1947 were likely exposed to H1N1 strains that are antigenically similar to the post-pandemic H1N1 strains that comprise most of our H1N1 infection data (Manicassamy et al., 2010; O’Donnell et al., 2012), creating the possibility that strain-specific cross-immunity drives the pattern we attribute to subtype-specific imprinting. These people nearly all fall into the 65 year-old age group in the study period. The study also underenrolled medically attended infections among people in nursing facilities, which would artificially lower the case count in this age group and may affect estimates of imprinting protection. Therefore, we excluded adults >65 years old and refitted our models. Excluding the oldest adults does not significantly change estimated imprinting protection or other parameters (Appendix 2 Table 1 and Table 2). When we exclude both the youngest and oldest age groups, initial infections by H1N1 and H3N2 have similar protective effects (Figure 5, bottom row). This shows that the combined effects of cross-protection in both the youngest and oldest individuals contribute to the signal of imprinting protection we observe, but they are not its sole drivers. ### VE varies by birth cohort in older children and adults The best-fitting model includes age-specific VE (Figure 4-Supplement 1; Appendix 2: Table 2). While serological responses to influenza vaccination are weakest in the young (Englund et al., 2005; Neuzil et al., 2006) and old (Lee et al., 2018; DiazGranados et al., 2014), it is unclear what age-related factors would drive variation in VE in other age groups. We hypothesized that VE in these ages varies with early exposure history, which correlates with birth year, rather than age. To test this hypothesis, we fitted a model with birth-cohort-specific VE to the cases, excluding either children <15 years old or adults >65 years old. We chose birth cohorts that corresponded to the age groups of the original model in 2017-2018 (Materials and Methods: “Vaccination”), keeping the number of parameters the same (e.g., VE in the 20-29 age group became VE in the 1988-1997 birth year cohort). We find that age-specific VE still outperforms all other models after we exclude the oldest age group (>65 years old). In contrast, birth-cohort-specific VE performs better when we exclude children <15 years old (Figure 6-Supplement 1). Estimates of imprinting protection and age-specific risk of medically attended influenza in the birth-cohort-specific VE models are not significantly different from estimates from the best-fitting model fitted to all ages (Appendix 2: Table 1). Taken together, these results suggest that birth-cohort-specific VE best explains the case distribution in older children and adults, who have likely experienced their first influenza infection, whereas age-specific VE best explains cases in younger children, who have less influenza exposure. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2020/03/29/19001875/F6) Figure 6. Estimates of birth-cohort-specific VE. Birth-cohort-specific VE differs significantly between subtypes and birth cohorts. The location of each pie chart represents the H3N2 (x-axis) and H1N1 (y-axis) VE estimates for a birth cohort (indicated by text) obtained from our model fitted to people >15 years old. Pie charts are colored by the probability of first infection by each subtype (i.e., imprinting probability). 95% confidence intervals of the VE estimates are indicated by light grey solid lines. The dashed grey line shows the diagonal where the VE estimate for H1N1 is equal to the VE estimate for H3N2. VE differs between birth cohorts that have similar imprinting by subtype (Figure 6; Appendix 2: Table 5). For example, the 1968-1977 and 1988-1997 cohorts have similar probabilities of primary exposure to H1N1 and H3N2, but they differ substantially in their VE to both subtypes (Figure 6). The 1988-1997 and 1998-2002 cohorts also have similar probabilities of primary exposure to each subtype and have similar H1N1 VEs, but have significantly different H3N2 VEs (Figure 6). Antigenic differences within each subtype might explain this variation. ### Discrepancies partly explained by antigenic evolution The best-fitting model accurately reproduces the age distributions of vaccinated and unvaccinated cases of each subtype, aggregated across seasons (Figure 7A). The only exception is that it underestimates aggregate H1N1 cases in unvaccinated 5-9 year-olds. By examining the differences between predicted and observed cases for each season, we see that this is largely driven by infection during the 2009 H1N1 pandemic (Figure 7B). Such a large antigenic change may have negated any protection from previous infection in 5-9 year-olds and made them particularly susceptible to pandemic infection. ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2020/03/29/19001875/F7) Figure 7. Model predictions compared to observed case counts. **A**. The model including age-specific VE and subtype-specific HA imprinting accurately predicts the overall age distribution of cases across seasons and age groups. Each row depicts the age distribution of cases among unvaccinated (top) and vaccinated (bottom) individuals over all sampled seasons (2007-2008 through 2017-2018). Each column indicates H1N1 cases (left, blue) and H3N2 cases (right, red). Open circles represent observed cases, solid lines represent the predicted number of cases from the best-fitting model, the shaded area represents the 95% prediction interval of the best-fitting model. **B**. Excess cases of dominant subtype for each season. Excess cases are defined as the predicted number of cases from the best-fitting model - observed cases (Appendix 1: “Calculating excess cases”). Each panel shows the excess cases of the dominant subtype for each season for each age group among unvaccinated (dark bars) and vaccinated (light bars) individuals. Grey error bars show the 95% prediction interval. The model underestimates cases in unvaccinated individuals who were 30-39 years old and over 50 years old in the 2013-2014 season (Figure 7B), as indicated by the many excess cases in these age groups in that season. This is further evidence that subtype-specific imprinting cannot explain all age variation. As mentioned before, this season provided one of the first examples that original antigenic sin could affect protection: middle-aged adults had been targeting a familiar site on the pandemic strain that then mutated, rendering them susceptible. Other age groups were effectively blind to these changes, owing to their different exposure histories (Linderman et al., 2014; Huang et al., 2015; Arriola et al., 2014; Dávila et al., 2014; Petrie et al., 2016). ## Discussion The distribution of influenza cases by birth year is consistent with subtype-level imprinting, whereby initial infection with a subtype protects against future medically attended infections by the same subtype. The stronger protective effect observed from primary H1N1 infection compared to primary H3N2 infection may be caused by stronger cross-protective responses to conserved epitopes in the more slowly evolving H1N1 (Bedford et al., 2015). This is in line with previous work showing that protection conferred by H1N1 infection lasts longer than protection conferred by H3N2 infection (Ranjeva et al., 2019). Another recent study found stronger imprinting protection from primary H1N1 compared to primary H3N2 infection (Gostic et al., 2019). Subtype-specific protection observed in seasonal influenza is narrower than the previously reported HA-group-level imprinting protection against avian influenza (Gostic et al., 2016), but in both cases, the protection correlates strongly with primary infection rather than any prior exposure. Examining cases of seasonal influenza over a 20-year period in Arizona, Gostic et al., 2019 find evidence of imprinting protection not only from HA but also NA, which we do not. We speculate that this discrepancy may be due to increasing vaccination coverage over time in middle-aged adults. During the period of the Arizona study (1993-1994 through 2014-2015), vaccination coverage in U.S. adults increased most rapidly in this age group (NHIS, 2009), which corresponds to the H2N2-imprinted cohorts near the end of the study. Without adjustment for vaccination, the apparently increased protection in the middle aged might resemble N2 imprinting. Accounting for vaccination in the MESA population, including the relatively stable vaccination coverage in each age group over time (Figure 1-Supplement 3), suggests imprinting protection is driven by HA. In contrast to the clear role of the imprinting subtype in protection against medically attended infection, the model implicates the imprinting strain or other attributes of early exposure history in VE. We expect that people born around the same time were likely exposed to similar strains, not just subtypes, of influenza A early in life, and our results support the idea that biases in immune memory from these early exposures (i.e., original antigenic sin; Davenport and Hennessy, 1957; Francis, 1960; Groth and Webster, 1966) influence VE. Specifically, we observe that our model is consistent with previous suggestions of birth-cohort-specific VE. The model with birth-cohort-specific VE better estimates cases in vaccinated 50-64 year-olds (born 1953-1967) in the 2015-2016 season than the model with age-specific VE, as indicated by the fewer excess cases predicted in that age group and an improved fit of 1.1 log-likelihood units (Figure 6-Supplement 2; Appendix 1: “Calculating excess cases”). Reduced VE in this group during the 2015-2016 season has been attributed to the exacerbation of antigenic mismatch by the vaccine in adults whose antibody responses were focused on a non-protective site (Skowronski et al., 2017b; Flannery et al., 2018). The improved performance of birth-cohort-specific VE relative to age-specific VE suggests other seasons and age groups where original antigenic sin might have influenced VE, such as 20-29 year-olds in the 2007-2008 influenza season. Although seasonal estimates of VE routinely stratify by age, shifts in VE from one season to the next might thus be easier to interpret in light of infection history (e.g., Skowronski et al., 2017b; Flannery et al., 2018). The results suggest this effect may be subtle, i.e., influenced by strains’ specific identities rather than merely their subtype. Our model cannot distinguish between the possibility that the precise identity of the imprinting strain primarily determines later VE, or if individuals’ responses to vaccination are shaped by a particular succession of exposures, which would be common to others in the same birth cohort. Regardless, variation in VE between birth cohorts appears substantial and presents a challenge for vaccination strategies (Erbelding et al., 2018). The use of different influenza vaccines in MESA during this period is unlikely to affect the results. Most people enrolled in the study received the standard-dose inactivated influenza vaccine (IIV-SD) (Figure 1-Supplement 7). However, between 9-26% of vaccinated children <18 years old received the live attenuated influenza vaccine (LAIV) between the 2008-2009 and 2015-2016 seasons (Figure 1-Supplement 7B). A separate study of LAIV VE in the United States found that LAIV and IIV-SD recipients who were repeat vaccinees (as most children were) had similar VE, and thus we do not expect that LAIV receipt should affect VE estimates (McLean et al., 2018). Similarly, 1-15% of adults >65 years old received the high-dose inactivated influenza vaccine (IIV-HD) between 2009-2010 and 2017-2018 (Figure 1-Supplement 7C). This vaccine is 20% more effective than IIV-SD (Lee et al., 2018). Therefore, the changing ratio of IIV-HD to IIV-SD recipients over time might bias results toward cohort-specific VE in models that include people >65 years old. However, when we fitted to cases between 15-64 years old, we found that cohort-specific VE still performed best. Thus, we conclude that changes in IIV-HD coverage do not substantially influence results. Potential methodological biases and the vaccination history of the study population nonetheless suggest caution in interpreting VE estimates. Selection and misclassification biases can arise when using influenza test-negative controls to control for differences in healthcare-seeking behavior (Lewnard et al., 2018; Sullivan et al., 2016). Because we also use test-negative controls to set our null expectation for the distribution of cases among birth cohorts, our VE estimates are subject to these biases as well. Moreover, since 45% of the study population is vaccinated, and most participants are frequent vaccinees (Figure 1-Supplement 6), we are limited in our ability to generalize the VE results to populations with much lower vaccination coverage and/or a shorter history of vaccination. Frequent vaccination has been associated with reduced VE (McLean et al., 2014; Saito et al., 2018; Skowronski et al., 2016). Therefore, the model may underestimate VE in less vaccinated populations. Underestimation of VE could also occur if unvaccinated people are protected by vaccination in the preceding season. Inference might also be distorted if vaccination has large indirect effects, which our model does not consider. Finally, our analysis is worth repeating in a larger population to reduce stochastic influences. We observed an unusually high H1N1 VE in the 1998-2002 birth cohort. Because we restricted cases in this analysis to people >15 years old, this VE estimate included data from only the 2013-2014 and 2015-2016 influenza seasons. No H1N1 cases among vaccinated or unvaccinated individuals were observed in this birth cohort in those seasons, which led to the high VE. This might have been due to particular epidemic dynamics in MESA. Incorporating differences in susceptibility based on early exposures might improve methods to forecast influenza seasons. The analysis of the relative risk of infection during the first half of each season shows more variation in the susceptible age groups from season to season than previously estimated (Worby et al., 2015). While the smaller sample sizes in MESA introduce uncertainty, the correlation between the relative risk and total fraction of cases indicates that the age groups driving epidemics indeed change from season to season. Because the contact structure of the population is probably constant over influenza seasons, variation in the driving age group may be determined by fluctuating susceptibility, which is partly determined by early infections. Therefore, incorporating information on early exposure history into epidemic models may allow for more accurate identification of at-risk populations and fine-scale epidemic timing. While the rate of antigenic evolution affects the rate at which different populations become susceptible to infection, we propose that the heterogeneity in susceptibility observed here may also drive antigenic evolution. Heterogeneity in susceptibility implies that influenza viruses face different selective pressures in groups with different exposure histories (Cobey and Hensley, 2017; Nakajima et al., 2000). Recent research consistent with this hypothesis has shown that sera isolated from different individuals can select for distinct escape mutants (Lee et al., 2019). More careful study of how immune memory to influenza evolves from infection and vaccination might improve understanding of influenza’s evolution. ## Ethics Human subjects: Study procedures for the VE study were approved by the IRB at the Marshfield Clinic Research Institute. Informed consent was obtained from all participants at the time of enrollment into the VE study. This analysis was subsequently approved by the IRB with a waiver of informed consent. The analysis of data was approved by the University of Chicago IRB. ## Data Availability Data for this study are available at [www.github.com/cobeylab/FluAImprinting](https://www.github.com/cobeylab/FluAImprinting) [https://github.com/cobeylab/FluAImprinting](https://github.com/cobeylab/FluAImprinting) ## Acknowledgments We thank Jennifer King and Carla Rottscheit for their assistance in providing the data for this study and Rohan Dandavati for compiling historical data on subtype frequencies and ILI. We thank Marcos Vieira and Kangchon Kim for their assistance in calculating imprinting probabilities. This work was completed with computational resources provided by the University of Chicago’s Research Computing Center. Funding for this project was provided by the National Institutes of Health (NIH), Department of Health and Human Services, under grant DP2AI117921 (to SC), CEIRS Contract No. HHSN272201400005C (to SC), and NRSA Fellowship F32AI145177-01 (to PA). HQM receives research support from Seqirus unrelated to this work. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. ## Appendix 1: Supplementary methods ### Vaccination coverage Seasonal influenza vaccination coverage for MESA was collected by age in the 2007-2008 through 2017-2018 seasons using a regional immunization registry (Irving et al., 2009). Monovalent vaccination coverage for the 2009-2010 season was obtained by directly measuring monovalent vaccination coverage in enrolled individuals and fitting a smoothing spline to the data (Figure 1-Supplement 3). We also calculated the fraction of people who received different vaccination formulations, and found that most people received IIV-SD (Figure 1-Supplement 7). ### Correlation of relative risk and fraction of cases To assess whether an age group’s relative risk correlates with the fraction of cases of that age group in the same season, we performed a rank correlation analysis. For each season, we ranked each age group based on its relative risk and the fraction of cases within that age group. If age groups were tied in either relative risk or fraction of cases, we assigned them the average rank they spanned. We then calculated the Pearson’s correlation coefficient for these two rankings. A positive correlation indicates that an age group with a large relative risk compared to other age groups will also make up a large proportion of cases compared to other age groups. ### Seasonal intensity We defined the intensity of an influenza season as the product of the mean fraction of patients with influenza-like illness (ILI) and the percentage of specimens testing positive for influenza A that season, ![Formula][33] where ILI*t* is the mean fraction of all patients with ILI in season *t* adjusted for differences in state population size (CDC, 2018), *F**t* is the number of respiratory specimens testing positive for influenza A in season *t*, and *N**t* is the total number of respiratory specimens tested in season *t*. For seasons 1997-1998 through 2017-2018, these data were obtained from the U.S. Outpatient Influenza-like Illness Surveillance Network (ILINet) and the World Health Organization/National Respiratory and Enteric Virus Surveillance System (WHO/NREVSS) Collaborating Labs (CDC, 2018). For seasons 1976-1977 through 1996-1997 when seasonal ILI data were not available, we assumed that the mean ILI was equal to the mean of mean ILI for seasons 1997-1998 through 2017-2018. We obtained data on *F**t* and *N**t* for these seasons from Thompson et al., 2003. We then normalized the intensity of each season by dividing *I**t* by the mean of *I**t* from the 1976-1977 through 2017-2018 seasons. For all seasons before 1976-1977, for which no seasonal intensity data were available, we assumed that the intensity of influenza A equalled the mean intensity of seasons 1976-1977 through 2017-2018. ### Fraction of season experienced We defined the fraction of a given influenza season *f**w*,*t* occurring in week *w* of season *t* as ![Formula][34] where ILI*w*,*t* is the weighted fraction of all patients with ILI in week *w* of season *t, F**w*,*t* is the number of respiratory specimens testing positive for influenza A in week *w* of season *t*, and *N**w*,*t* is the number of specimens tested in week *w* of season *t*. ![Graphic][35] is the product of ILI and the fraction of positive influenza A specimens summed over all weeks of the influenza season *t*, where *w* is the first week of the season and *w**f* is the final week of the season. We defined the start of the influenza season as week 40 of the calendar year, which usually falls at the beginning of October. For seasons before 1997-1998, where weekly data is unavailable, we assumed that the fraction of the influenza season experienced in week *w* was ![Formula][36] where ![Graphic][37] is the mean fraction of the influenza season experienced at week *w* for all seasons after 1997-1998. We used *f**w*,*t* to calculate the fraction of an influenza season experienced by an individual born in year *y*. We assumed that people born in year *y* were born uniformly throughout the year. We also assumed that due to maternal immunity, infants did not experience immunizing exposure to influenza until they were at least 180 days old. Let *p**y*,*w*,*t* be the proportion of individuals born in year *y* that are over 180 days old in week *w* of season *t* and *y**y*,*t* be the fraction of individuals born in year *y* exposed to influenza season *t*. Then ![Formula][38] ### Calculating the fraction unexposed When calculating imprinting probabilities, we used an iterative approach to calculate *U**y*,*t* the fraction of people in birth cohort *y* who were unexposed at the start of season *t*. First, we assumed that in the first year of life (i.e., when *t* = *y*), the entire population was unexposed. For seasons where *t > y*, the fraction unexposed depends on the fraction unexposed at the start of the previous season (*U**y*,*t*−1) and the attack rate in the previous season (*a**y*,*t*−1). Thus, ![Formula][39] ### Birth year distribution of the study population In order to convert the demographic age distribution to a birth year distribution, we assumed that people were born uniformly throughout the year. We defined a breakpoint date prior to the start of the enrollment period based on when the the 6 month-old age limit cutoff was set (e.g., if the breakpoint date was Ocotober 1, then infants had to be 6 months old by that date to be eligible for enrollment). We used this date to calculate the fraction of people of age *a* in season *t* who were born in year *t* – *y* ![Graphic][40] or year ![Graphic][41]. A fraction ![Graphic][42] of the total population of age *a* in season *t* was assigned to birth year *t* − *y* and ![Graphic][43] to *t* − *y* − 1. Breakpoint dates ranged from September 1 through January 1 with the exception of the pandemic season which had a breakpoint date of May 1, 2009. The start of the enrollment period ranged from December to January with the exception of the 2009 pandemic season, when enrollment began in May 2009. ### Fraction of birth cohort with specific age When converting an age-specific parameter to a birth-cohort-specific parameter as in Materials and Methods “Age-specific factors”, we considered that each birth cohort had two possible ages (*a*1 and *a*2) in a given season *t*. We assumed that people were born uniformly throughout the year and used the same breakpoint dates described above in “Birth year distribution of the study population.” Then, *f* (*a*1, *t, y*), the fraction of people born in year *y* who were age *a*1 in season *t*, is the fraction of people born in year *y* who were born on a date prior to the breakpoint date for season *t*. Finally, *f* (*a*2, *t, y*), the fraction of people born in year *y* who were age *a*2 in season *t*, is 1 - *f* (*a*1, *t, y*). ### Age-specific rates of approachment, enrollment, and nursing home residence The relative rates at which different age groups were approached for study enrollment (the approachment rate, *x*approach) varied between seasons. Similarly, the relative rates at which different age groups enrolled in the study after being approached (the enrollment rate, *x*enroll) also varied between seasons. Enrollment rates also varied between vaccinated and unvaccinated individuals. We defined the approachment rate of an age group *g* in season *t* as ![Formula][44] where *N*approached,*t*,*g* is the number of people in age group *g* during season *t* who were approached for enrollment, and *N*MAARI,*t*,*g* is the total number of people in the MESA cohort who presented with MAARI regardless of whether they were approached for enrollment. We defined the enrollment rate of age group *g* in season *t* with vaccination status *v* as ![Formula][45] where *N*enrolled,*t*,*g*,*v* is the number of people in age group *g* with vaccination status *v* who enrolled in the study in season *t*, and *N*approached,*t*,*g*,*v* is the number of people in age group *g* with vaccination status *v* who were approached for enrollment in season *t*. Due to differences in data collection for the 2007-2008 and 2008-2009 seasons, complete vaccination records for eligible unenrolled individuals were not available, so we assumed that the enrollment rates by age group and vaccination status in those seasons were equal to the mean enrollment rate for each age group and vaccination status across all other seasons. We normalized *x*approach,*t*,*g* by the value of *x*approach,*t*,*g* for the reference age group (i.e., 20-29 year-olds) in each season. Similarly, we normalized *x*enroll,,*t*,*g*,*v* to the value of *x*enroll,,*t*,*g*,*v* for unvaccinated members of the reference age group for each season. This yielded the relative approachment and enrollment rates ![Graphic][46] and ![Graphic][47]. We converted both ![Graphic][48] and ![Graphic][49] to birth-year specific covariates (i.e. covariates by *y* instead of *g*) using the same procedure described in Materials and Methods: “Age-specific factors” (Equation 9). Finally, the study did not enroll residents of skilled nursing facilities with dedicated medical staff. To account for this, we estimated the proportion of the population in nursing facilities within the study area. We obtained the total number of beds in nursing facilities within MESA in 2018 from the Wisconsin Department of Health Services (WDHS, 2018). We assumed that the total number of beds did not change between 2007-2008 and 2017-2018. We also used data from the Centers for Medicare and Medicaid Services (CMS, 2015) to calculate the percent of beds occupied in Wisconsin nursing facilities by age for 2011 through 2014 and the fraction of people in a nursing facility by age group. We used a smoothing spline to obtain the fraction of people of a given age in a nursing facility. For seasons before 2010-2011 and after 2013-2014, we assumed that the fraction of people of a given age in a nursing facility was the average value for 2011-2014. Given the total population of the study area by age and season, we calculated the fraction of people in a given age *a* and season *t* who are in nursing facilities (*k**t*,*a*). We converted this to a covariate by birth year (*k**t*,*y*) using the same procedure described in Materials and Methods: “Age-specific factors” (Equation 9). ### Evaluation of predictive power To evaluate the predictive power of each model, we performed leave-one-out cross-validation. We excluded data from each season and fitted our models to the remaining seasons. Because our goal was to evaluate how well our models predict seasonal epidemics, we excluded the 2009 pandemic season from all cross-validation analyses. We also did not test seasonal VE models with cross-validation since estimation of seasonal VE requires data from the excluded season. Let *n**s*,*t*,*y*,*v* be the number of observed cases of subtype *s* in season *t* among people born in year *y* with vaccination status *v*. Let ![Graphic][50] be the multinomial probability of a case of subtype *s* in season *t*− among people born in year *y* with vaccination status *v* under model *M* fitted to all seasons except *t*−. Let ![Graphic][51] be the total number of cases of subtype *s* in season *t*−. Then, the predicted number of cases of subtype *s* in season *t*− among people born in year *y* with vaccination status *v* under model *M* fitted to all seasons except *t*− is ![Formula][52] The sum squared prediction error for model *M* in season *t*− is given by ![Formula][53] where *y*max,*t* is the maximum possible birth year in season *t*. We evaluated each model *M* by its mean-squared prediction error across all excluded seasons *t*−. Let − be the set of all seasons left out and *X* be the size of *T*−. Then the mean-squared prediction error for model *M* is ![Formula][54] ### Sensitivity to uncertainty in ILI and the frequency of influenza A Because of the lack of ILI data prior to the 1997-1998 season and the lack of data on the frequency of influenza A prior to the 1976-1977 season, we used simulated datasets to test the robustness of our results. We randomly assigned ILI values from the 1997-1998 through 2017-2018 seasons to every season which did not have a measured ILI value. Similarly, we randomly assigned values of the frequency of influenza A from the 1976-1977 through 2017-2018 seasons to every season which did not have a measured value for the frequency of influenza A. We created 10000 simulated datasets using this procedure and recalculated imprinting probabilities for each dataset (Figure 3-Supplement 2). In the period of H1N1 and H3N2 co-circulation, the maximum H1N1 imprinting probability for a particular birth cohort corresponds to the minimum H3N2 imprinting probability for that cohort and vice-versa. Therefore, to generate datasets representing the upper and lower bounds of imprinting probabilities, we assigned imprinting probabilities from the simulation with either the lowest or highest H1N1 imprinting probability to each birth cohort in each season. We then fitted our models to these two datasets and evaluated model fit using cAIC. ### Sensitivity to age groups To test whether our models were sensitive to our choice of age groups, we fit revised versions of all our models with different age groups: * 0-4 years, 5-17 years, 18-49 years, 50-64 years, and >65 years * 0-4 years, 5-17 years, 18-64 years, and >65 years These models with alternate age groupings were fitted to case data to determine whether our findings on the strength of protection from initial H1N1 and H3N2 infection were altered compared to fits using the higher-resolution age grouping described above (Appendix 2: Table 4). ### Sensitivity to sampling effort Sampling effort was not even across seasons, and analysis of the number of influenza cases per sampling day suggested that a significant number of cases may have been missed at the beginning or end of a specific seasons (Figure 5-Supplement 3). As our analysis of relative risk indicates, different age groups are more susceptible during different points in the influenza season, and therefore missing data from the beginning or end of a season could introduce bias in the observed age distribution of cases. To adjust for this, we simulated cases for seasons which did not have sufficient sampling of the start or end of the epidemic period. We considered a season sufficiently sampled if the sampling period spanned the start and end of the epidemic. We expect that the start and end of the epidemic have few cases per sampling day, and we therefore defined sufficiently sampled seasons as seasons where * the number of cases per sampling day in the first week of the enrollment period is <1 and * the number of cases per sampling day in the last week of the enrollment period is <1. To extrapolate the start of a season, we linearly regressed the number of cases of the dominant subtype per sampling day for each week of the first half of the season and identified the week of the season where the number of cases per sampling day fell below 1 (*t*). For each week from *t* to the first week of the enrollment period, we used the regression of cases per sampling day to calculate the number of cases we expected to see in each week. Summing these yields the total number of unsampled cases at the beginning of the season. We used a similar approach to extrapolate the number of unsampled cases at the end of a season by instead regressing cases per sampling day for each week of the latter half of the season. We did not extrapolate cases for the 2010-2011 season for this analysis since the observed number of cases per sampling day did not follow a typical epidemic curve. We stochastically assigned a birth year and vaccination status to these cases according to a multinomial distribution. The success probabilities of this distribution were set using the age distribution of cases of the dominant subtype from the first two weeks of the enrollment period (if extrapolating the beginning of a season) or the last two weeks of the enrollment period (if extrapolating the end of a season). Specifically, we calculated the distribution of observed cases in the first or last two weeks of the enrollment period among nine age groups (Materials and Methods: “Age-specific factors”) with their associated vaccination status. We then assumed that cases were uniformly distributed among all birth years contained in an age group. This yielded a set of probabilities describing the probability of infection given birth year and vaccination status in a specific season. We sampled from these multinomial distributions 1000 times to obtain augmented datasets that combined observed and extrapolated cases. For each replicate simulation, we calculated the age distribution of cases for the entire season as well as the relative risk of each age group in the first versus the latter half of the season (Figure 2-Supplement 2B). We also fitted the best-fitting model to 100 of these datasets (excluding the 2010-2011 season) and recorded the estimated imprinting strength for both H1N1 and H3N2 for each fit (Figure 5-Supplement 4). ### Calculating excess cases We defined excess cases for a given birth cohort or age group as the number of observed cases for that birth cohort or age group minus the number of predicted cases for that age group. Predictions were obtained by multiplying the multinomial probabilities produced by the model by the total number of cases of the dominant subtype in each season. A 95% prediction interval was obtained by simulating 10000 datasets using the multinomial probabilities from a specific model (Figure 6-Supplement 2, Figure 7). To test whether recent infection might be confounding our estimates, we calculated the correlation between excess cases in each birth cohort in each season with excess cases of the same birth cohort in the next season with the same dominant subtype (Figure 5-Supplement 5). ## Appendix 2: Supplementary tables and figures View this table: [Appendix 2 Table 1.](http://medrxiv.org/content/early/2020/03/29/19001875/T1) Appendix 2 Table 1. Estimates of parameters shared by the age-specific VE and birth-cohort-specific VE models. View this table: [Appendix 2 Table 2.](http://medrxiv.org/content/early/2020/03/29/19001875/T2) Appendix 2 Table 2. Estimates of age-specific VE parameters in models fitted to different age groups. View this table: [Appendix 2 Table 3.](http://medrxiv.org/content/early/2020/03/29/19001875/T3) Appendix 2 Table 3. Estimates of imprinting protection fitted to datasets representing upper and lower bounds of imprinting probabilities. View this table: [Appendix 2 Table 4.](http://medrxiv.org/content/early/2020/03/29/19001875/T4) Appendix 2 Table 4. Estimates of imprinting protection for models with different age groups. View this table: [Appendix 2 Table 5.](http://medrxiv.org/content/early/2020/03/29/19001875/T5) Appendix 2 Table 5. Estimates for VE from model with birth-cohort-specific VE fitted to people >15 years old. ![Figure 1–Supplement 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F8.medium.gif) [Figure 1–Supplement 1.](http://medrxiv.org/content/early/2020/03/29/19001875/F8) Figure 1–Supplement 1. Sample collection and final study population. Flowchart of sample collection (**A**) and final study population stratified by season, age, test status, and vaccination status (**B**). “Test-positive” is defined as testing positive for the dominant circulating influenza A subtype in that season. ![Figure 1–Supplement 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F9.medium.gif) [Figure 1–Supplement 2.](http://medrxiv.org/content/early/2020/03/29/19001875/F9) Figure 1–Supplement 2. Birth year distribution of population. Each panel shows the population distribution of all individuals in the study area who met the age criteria for study enrollment. People under 6 months old at the start of the sampling period in a season were not eligible to participate. ![Figure 1–Supplement 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F10.medium.gif) [Figure 1–Supplement 3.](http://medrxiv.org/content/early/2020/03/29/19001875/F10) Figure 1–Supplement 3. Vaccination coverage. We estimated monovalent vaccination coverage in 2009-2010 by measuring vaccination coverage among enrolled people and fitting a smoothing spline to the data (solid line). ![Figure 1–Supplement 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F11.medium.gif) [Figure 1–Supplement 4.](http://medrxiv.org/content/early/2020/03/29/19001875/F11) Figure 1–Supplement 4. Age distribution of high-risk medical status. High-risk medical status (Materials and Methods, “Study cohort”) varies with age and vaccination status but stays relatively consistent across seasons. Each plot shows the fraction of enrolled people who had a high-risk medical condition for each season stratified by age, vaccination status, and test status. High-risk medical condition data was not collected for the 2009 pandemic season. ![Figure 1–Supplement 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F12.medium.gif) [Figure 1–Supplement 5.](http://medrxiv.org/content/early/2020/03/29/19001875/F12) Figure 1–Supplement 5. Rate of MAARI in vaccinated and unvaccinated controls. Vaccinated individuals seek health-care for MAARI at a higher rate than predicted by vaccination coverage. We measured the fraction of vaccinated people among all who presented with MAARI and tested negative for influenza ![Graphic][55] Materials and Methods: “Vaccination”). This is plotted against vaccination coverage by season for different age groups. The dashed grey line shows where *R* and vaccination coverage are equal. Vaccination coverage for the 2009-2010 season uses monovalent vaccination coverage estimated directly from all individuals with MAARI. We do not show the 2009 pandemic season because the monovalent vaccine was not distributed until the second wave of the pandemic. ![Figure 1–Supplement 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F13.medium.gif) [Figure 1–Supplement 6.](http://medrxiv.org/content/early/2020/03/29/19001875/F13) Figure 1–Supplement 6. Repeat vaccination by age group and season. Each bar shows the fraction of individuals who were vaccinated in that season who also received at least one influenza vaccination in the previous two seasons. ![Figure 1–Supplement 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F14.medium.gif) [Figure 1–Supplement 7.](http://medrxiv.org/content/early/2020/03/29/19001875/F14) Figure 1–Supplement 7. Vaccine type receive. Most vaccinated study participants received the inactivated influenza vaccine. The fraction of vaccinated people who received the standard-dose inactivated influenza vaccine (IIV-SD), the high-dose inactivated influenza vaccine (IIV-HD), or the live attenuated influenza vaccine (LAIV) is shown for all participants (**A**), children < 18 years old (**B**), and adults >65 years old (**C**). ![Figure 2–Supplement 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F15.medium.gif) [Figure 2–Supplement 1.](http://medrxiv.org/content/early/2020/03/29/19001875/F15) Figure 2–Supplement 1. Statistical analysis of age distribution of cases. Seasons differ in their age distributions. The color intensity of each cell shows the observed G-test statistic, which measures how much the age distributions of two seasons differ from the null expectation that they are drawn from the same distribution (Materials and Methods: “Calculating differences in the age distribution between seasons.”). The text in each cell shows the Bonferroni-corrected p-value for each G-test. The dominant subtype of each season is indicated by the label color. ![Figure 2–Supplement 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F16.medium.gif) [Figure 2–Supplement 2.](http://medrxiv.org/content/early/2020/03/29/19001875/F16) Figure 2–Supplement 2. Correlation of relative risk and fraction of cases within an age group. **A**. Each point shows the rank of an age group’s relative risk of infection during the first half compared to the second half of an epidemic period (x-axis) and the rank of the fraction of cases belonging to that age group in the same epidemic period (y-axis) (Appendix 1: “Correlation of relative risk and fraction of cases”). Points are colored by the dominant subtype of the season and x-axis values are offset to facilitate visualization. Points with the same x and y values overlap and are indicated by darker shading. **B**. To account for potential undersampling of cases at the beginning and end of specific seasons, we simulated 1000 replicate epidemics (Appendix 1 : “Sensitivity to sampling effort”) and calculated the same correlation as in panel **A**. The range is indicated by a vertical line and the median by a square. ![Figure 2–Supplement 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F17.medium.gif) [Figure 2–Supplement 3.](http://medrxiv.org/content/early/2020/03/29/19001875/F17) Figure 2–Supplement 3. Relative risk among different age groups across seasons. Each panel shows the relative risk of infection in the first versus the second half of an epidemic for different age groups in each season (Materials and Methods: “Calculating relative risk”). Relative risk greater than 1 (indicated by the grey dashed line) means that an age group was more likely to be infected at during the first rather than second half of an epidemic. Age groups with no cases in the latter half of a season are indicated by asterisks and no bar. The dominant subtype of each subtype is indicated by the bar color. 95% binomial confidence intervals are indicated by grey vertical lines. Bars with asterisks over them indicate that the 95% confidence interval includes the scenario where all cases occur in the first half of the season. ![Figure 3–Supplement 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F18.medium.gif) [Figure 3–Supplement 1.](http://medrxiv.org/content/early/2020/03/29/19001875/F18) Figure 3–Supplement 1. Intensity and subtype frequencies of influenza **A**. The intensity (top panel) and subtype frequencies (bottom panel) of influenza A seasons in the United States. Intensity is measured as the product of influenza-like illness (ILI) and the fraction of respiratory specimens testing positive for influenza A in national surveillance data (Appendix 1: “Seasonal intensity”). This is normalized to the average intensity value between 1977 and 2017-2018. Seasons before 1977 where United States ILI surveillance data are unavailable are assumed to have an intensity score of 1 (i.e., the average score over all other seasons). Subtype frequencies were obtained from national surveillance data before the 2007-2008 season and directly from the MESA studies afterwards. ![Figure 3–Supplement 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F19.medium.gif) [Figure 3–Supplement 2.](http://medrxiv.org/content/early/2020/03/29/19001875/F19) Figure 3–Supplement 2. Imprinting probabilities with random sampling of seasonal intensity. Uncertainty in ILI and the frequency of A have a small impact on imprinting probabilities. We simulated 10000 datasets to represent the range of possible epidemic sizes for seasons where we did not have data on either ILI or the frequency of influenza A (Appendix 1: “Sensitivity to uncertainty in ILI and the frequency of influenza A”). The vertical dashed line shows the point at which data on ILI and the frequency of influenza A are available while the vertical dotted line shows the point at which data on only the frequency of A is available. The median imprinting probabilities for those simulations is shown as a solid line with the maximum and minimum imprinting probabilities shown by the bounds of the shaded area. ![Figure 4–Supplement 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F20.medium.gif) [Figure 4–Supplement 1.](http://medrxiv.org/content/early/2020/03/29/19001875/F20) Figure 4–Supplement 1. Ranking of models fitted to all ages. The best-fitting model includes age-specific risk of medically attended influenza A infection, HA subtype imprinting, and age-specific VE. The 11 main models are shown as rows with colored squares indicating whether that model included parameters indicated by the columns. Orange squares indicate covariates that were not estimated. Light green squares mean that a given estimated parameter was supported. Dark green squares mean that the model did not support the inclusion of the parameters indicated by the column (i.e., the CI includes 0). Models are sorted by their cAIC relative to the best-fitting model. ![Figure 5–Supplement 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F21.medium.gif) [Figure 5–Supplement 1.](http://medrxiv.org/content/early/2020/03/29/19001875/F21) Figure 5–Supplement 1. Ranking of models by predictive power. The model which best-predicts excluded seasons includes HA subtype imprinting and age-specific VE. Models are shown as rows with colored squares indicating whether that model included parameters indicated by the columns. Orange squares indicate covariates that were not estimated. Light green squares mean that a given estimated parameter was supported. Dark green squares mean that the model did not support the inclusion of the parameters indicated by the column (i.e., the CI includes 0). Models are sorted by their MSE in predicting excluded seasons (Appendix 1: “Evaluation of predictive power”). ![Figure 5–Supplement 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F22.medium.gif) [Figure 5–Supplement 2.](http://medrxiv.org/content/early/2020/03/29/19001875/F22) Figure 5–Supplement 2. Model performance on excluded seasons. Each panel shows the number of observed and predicted cases by birth year among unvaccinated (**A**) and vaccinated (**B**) study participants. Predictions and 95% prediction intervals were generated by fitting the model including age-specific risk of medically attended influenza A infection, HA subtype imprinting, and age-specific VE fitted to all seasons except the season in the panel (Appendix 1: “Evaluation of predictive power”). ![Figure 5–Supplement 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F23.medium.gif) [Figure 5–Supplement 3.](http://medrxiv.org/content/early/2020/03/29/19001875/F23) Figure 5–Supplement 3. Cases per sampling day. Each panel shows the number of cases per sampling day (green circles). We extrapolated cases at the start and end of the season (orange dashed line) if the observed number of cases per day exceeded 1 (black line) at the start and end of that season (Appendix 1: “Sensitivity to sampling effort”). ![Figure 5–Supplement 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F24.medium.gif) [Figure 5–Supplement 4.](http://medrxiv.org/content/early/2020/03/29/19001875/F24) Figure 5–Supplement 4. Estimates of imprinting protection with added simulated cases. We fitted the model including HA subtype imprinting and age-specific VE to simulated cases in seasons where the enrollment period does not fully overlap the epidemic period and recorded the maximum likelihood estimates for H1N1 and H3N2 imprinting protection (Appendix 1: “Sensitivity to sampling effort”). The distributions of these values are shown as violin plots and the medians are shown as squares. Estimates of imprinting protection from the best-fitting model without simulated data with a 95% confidence interval are shown as circles with error bars. ![Figure 5–Supplement 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F25.medium.gif) [Figure 5–Supplement 5.](http://medrxiv.org/content/early/2020/03/29/19001875/F25) Figure 5–Supplement 5. Correlation of excess cases between seasons. We tested whether excess cases in each birth cohort were negatively correlated with excess cases in the same birth cohort in the next season of the same subtype (Appendix 1: “Calculating excess cases”). We find a weak positive correlation for cases of H1N1 (Spearman’s *p*=0.12, 95% CI 0.02-0.22) and H3N2 (Spearman’s *p*=0.05, 95% CI -0.03-0.14). ![Figure 6–Supplement 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F26.medium.gif) [Figure 6–Supplement 1.](http://medrxiv.org/content/early/2020/03/29/19001875/F26) Figure 6–Supplement 1. Ranking of models fitted to people >15 years old. A model including age-specific risk of medically attended influenza A infection, HA subtype imprinting, and birth-cohort-specific VE best fits cases of people >15 years old. The 11 main models are shown as rows with colored squares indicating whether that model uses parameters indicated by the columns. Orange squares indicate covariates that were not estimated. Light green squares mean that a given estimated parameter was supported. Dark green squares mean that the model did not support the inclusion of the parameters indicated by the column (i.e., the CI includes 0). Models are sorted by their cAIC relative to the best-fitting model. ![Figure 6–Supplement 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/29/19001875/F27.medium.gif) [Figure 6–Supplement 2.](http://medrxiv.org/content/early/2020/03/29/19001875/F27) Figure 6–Supplement 2. Excess cases for models using birth-cohort-specific VE and age-specific VE. The birth-cohort-specific VE model predicts observed cases better than the age-specific VE model for people >15 years old. Bars show the excess cases in vaccinated individuals relative to the birth-cohort-specific VE model (dark colors) and the age-specific VE model (light colors) for age groups >15 years old. Colors indicate the dominant subtype of a given season. 95% prediction intervals are shown as grey error bars. * Received July 9, 2019. * Revision received March 24, 2020. * Accepted March 29, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. Ann J, Papenburg J, Bouhy X, Rheaume C, Hamelin ME, Boivin G. Molecular and antigenic evolution of human influenza A/H3N2 viruses in Quebec, Canada, 2009-2011. J Clin Virol. 2012; 53(1):88–92. doi: 10.1016/j.jcv.2011.09.016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jcv.2011.09.016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22036039&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 2. Arriola CS, Brammer L, Epperson S, Blanton L, Kniss K, Mustaquim D, Steffens C, Dhara R, Leon M, Perez A, Chaves SS, Katz J, Wallis T, Villanueva J, Xu X, Abd Elal AI, Gubareva L, Cox N, Finelli L, Bresee J, et al. Update: influenza activity - United States, September 29, 2013-February 8, 2014. MMWR Morbidity and mortality weekly report. 2014; 63(7):148–154. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24553198&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 3. Beauté J, Zucs P, Korsun N, Bragstad K, Enouf V, Kossyvakis A, Griškevicius A, Olinger CM, Meijer A, Guiomar R, Prosenc K, Staroiiová E, Delgado C, Brytting M, Broberg E, European Influenza Surveillance N. Age-specific differences in influenza virus type and sub-type distribution in the 2012/2013 season in 12 European countries. Epidemiology and infection. 2015; 143(14):2950–2958. doi: 10.1017/S0950268814003422. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268814003422&link_type=DOI) 4. Bedford T, Riley S, Barr IG, Broor S, Chadha M, Cox NJ, Daniels RS, Gunasekaran CP, Hurt AC, Kelso A, Klimov A, Lewis NS, Li X, McCauley JW, Odagiri T, Potdar V, Rambaut A, Shu Y, Skepner E, Smith DJ, et al. Global circulation patterns of seasonal influenza viruses vary with antigenic drift. Nature. 2015; 523:217. doi: 10.1038/nature14460. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature14460&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26053121&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 5. Belongia EA, Kieke BA, Donahue JG, Coleman LA, Irving SA, Meece JK, Vandermause M, Lindstrom S, Gargiullo P, Shay DK. Influenza vaccine effectiveness in Wisconsin during the 2007-08 season: comparison of interim and final results. Vaccine. 2011; 29(38):6558–63. doi: 10.1016/j.vaccine.2011.07.002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2011.07.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21767593&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295148800021&link_type=ISI) 6. Belongia EA, Kieke BA, Donahue JG, Greenlee RT, Balish A, Foust A, Lindstrom S, Shay DK. Effectiveness of Inactivated Influenza Vaccines Varied Substantially with Antigenic Match from the 2004–2005 Season to the 2006–2007 Season. The Journal of Infectious Diseases. 2009; 199(2):159–167. doi: 10.1086/595861. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/595861&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19086915&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262086500002&link_type=ISI) 7. Biggerstaff M, Jhung MA, Reed C, Fry AM, Balluz L, Finelli L. Influenza-like illness, the time to seek healthcare, and influenza antiviral receipt during the 2010–2011 influenza season—United States. The Journal of infectious diseases. 2014; 210(4):535–544. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiu224&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24731959&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 8. Bodewes R, de Mutsert G, van der Klis FRM, Ventresca M, Wilks S, Smith DJ, Koopmans M, Fouchier RAM, Osterhaus ADME, Rimmelzwaan GF. Prevalence of Antibodies against Seasonal Influenza A and B Viruses in Children in Netherlands. Clinical and Vaccine Immunology. 2011; 18(3):469–476. doi: 10.1128/cvi.00396-10. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/cvi.00396-10&link_type=DOI) 9. Brooks-Pollock E, Tilston N, Edmunds WJ, Eames KTD. Using an online survey of healthcare-seeking behaviour to estimate the magnitude and severity of the 2009 H1N1v influenza epidemic in England. BMC Infectious Diseases. 2011; 11(1):68. doi: 10.1186/1471-2334-11-68. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2334-11-68&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21410965&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 10. Budd AP, Beacham L, Smith CB, Garten RJ, Reed C, Kniss K, Mustaquim D, Ahmad FB, Cummings CN, Garg S, Levine MZ, Fry AM, Brammer L. Birth Cohort Effects in Influenza Surveillance Data: Evidence that First Influenza Infection Affects Later Influenza-Associated Illness. The Journal of Infectious Diseases. 2019; doi: 10.1093/infdis/jiz201. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiz201&link_type=DOI) 11. Caini S, Spreeuwenberg P, Kusznierz GF, Rudi JM, Owen R, Pennington K, Wangchuk S, Gyeltshen S, Ferreira de Almeida WA, Pessanha Henriques CM, Njouom R, Vernet MA, Fasce RA, Andrade W, Yu H, Feng L, Yang J, Peng Z, Lara J, Bruno A, et al. Distribution of influenza virus types by age using case-based global surveillance data from twenty-nine countries, 1999-2014. BMC infectious diseases. 2018; 18(1):269–269. doi: 10.1186/s12879-018-3181-y. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-018-3181-y&link_type=DOI) 12. CDC, FluView National, Regional, and State Level Outpatient Illness and Viral Surveillance. Accessed 23-October-2018; 2018. [https://gis.cdc.gov/grasp/fluview/fluportaldashboard.html](https://gis.cdc.gov/grasp/fluview/fluportaldashboard.html). 13. CMS, Nursing Home Compendium 2015 Edition. Accessed 30-September-2019; 2015.[https://www.cms.gov/Medicare/Provider-Enrollment-and-Certification/CertificationandComplianc/downloads/nursinghomedatacompendium\_508-2015.pdf](https://www.cms.gov/Medicare/Provider-Enrollment-and-Certification/CertificationandComplianc/downloads/nursinghomedatacompendium_508-2015.pdf). 14. Cobey S, Koelle K. Capturing escape in infectious disease dynamics. Trends Ecol Evol. 2008; 23(10):572–7. doi: 10.1016/j.tree.2008.06.008. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tree.2008.06.008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18715671&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000260144800009&link_type=ISI) 15. Cobey S, Hensley SE. Immune history and influenza virus susceptibility. Current Opinion in Virology. 2017; 22:105–111. doi: [https://doi.org/10.1016/j.coviro.2016.12.004](https://doi.org/10.1016/j.coviro.2016.12.004). 16. Davenport FM, Hennessy AV. A serologic recapitulation of past experiences with influenza A; antibody response to monovalent vaccine, vol. 104; 1956. doi: 10.1084/jem.104.1.85. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjg6IjEwNC8xLzg1IjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMjkvMTkwMDE4NzUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 17. Davenport FM, Hennessy AV. Predetermination by infection and by vaccination of antibody response to influenza virus vaccines. J Exp Med. 1957; 106(6):835–50. doi: 10.1084/jem.106.6.835. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjk6IjEwNi82LzgzNSI7czo0OiJhdG9tIjtzOjM5OiIvbWVkcnhpdi9lYXJseS8yMDIwLzAzLzI5LzE5MDAxODc1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 18. Davis AKF, McCormick K, Gumina ME, Petrie JG, Martin ET, Xue KS, Bloom JD, Monto AS, Bushman FD, Hensley SE. Sera from Individuals with Narrowly Focused Influenza Virus Antibodies Rapidly Select Viral Escape Mutations In Ovo. Journal of Virology. 2018; 92(19). doi: 10.1128/JVI.00859-18. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjE1OiI5Mi8xOS9lMDA4NTktMTgiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMy8yOS8xOTAwMTg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 19. DiazGranados CA, Dunning AJ, Kimmel M, Kirby D, Treanor J, Collins A, Pollak R, Christoff J, Earl J, Landolfi V. Efficacy of high-dose versus standard-dose influenza vaccine in older adults. New England Journal of Medicine. 2014; 371(7):635–645. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1315727&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25119609&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340819600009&link_type=ISI) 20. Dávila J, Chowell G, Borja-Aburto VH, Viboud C, Grajales Muñiz C, Miller M. Substantial Morbidity and Mortality Associated with Pandemic A/H1N1 Influenza in Mexico, Winter 2013-2014: Gradual Age Shift and Severity. PLoS currents. 2014; 6:ecur- rents.outbreaks.a855a92f19db1d90ca955f5e908d6631. doi: 10.1371/currents.outbreaks.a855a92f19db1d90ca955f5e908d6631. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/currents.outbreaks.a855a92f19db1d90ca955f5e908d6631&link_type=DOI) 21. Englund JA, Walter EB, Fairchok MP, Monto AS, Neuzil KM. A comparison of 2 influenza vaccine schedules in 6- to 23-month-old children. Pediatrics. 2005; 115(4):1039–47. doi: 10.1542/peds.2004-2373. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBlZGlhdHJpY3MiO3M6NToicmVzaWQiO3M6MTA6IjExNS80LzEwMzkiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMy8yOS8xOTAwMTg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 22. Erbelding EJ, Post DJ, Stemmy EJ, Roberts PC, Augustine AD, Ferguson S, Paules CI, Graham BS, Fauci AS. A Universal Influenza Vaccine: The Strategic Plan for the National Institute of Allergy and Infectious Diseases. The Journal of Infectious Diseases. 2018 02; 218(3):347–354. doi: 10.1093/infdis/jiy103. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiy103&link_type=DOI) 23. Flannery B, Smith C, Garten RJ, Levine MZ, Chung JR, Jackson ML, Jackson LA, Monto AS, Martin ET, Belongia EA, McLean HQ, Gaglani M, Murthy K, Zimmerman R, Nowalk MP, Griffin MR, Keipp Talbot H, Treanor JJ, Wentworth DE, Fry AM. Influence of Birth Cohort on Effectiveness of 2015-2016 Influenza Vaccine Against Medically Attended Illness Due to 2009 Pandemic Influenza A(H1N1) Virus in the United States. J Infect Dis. 2018; 218(2):189–196. doi: 10.1093/infdis/jix634. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jix634&link_type=DOI) 24. Fonville JM, Fraaij PLA, de Mutsert G, Wilks SH, van Beek R, Fouchier RAM, Rimmelzwaan GF. Antigenic Maps of Influenza A(H3N2) Produced With Human Antisera Obtained After Primary Infection. The Journal of Infectious Diseases. 2015; 213(1):31–38. doi: 10.1093/infdis/jiv367. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiv367&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26142433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 25. Francis T. On the doctrine of original antigenic sin. Proceedings of the American Philosophical Society. 1960; 104(6):572–578. 26. Gaglani M, Pruszynski J, Murthy K, Clipper L, Robertson A, Reis M, Chung JR, Piedra PA, Avadhanula V, Nowalk MP, Zimmerman RK, Jackson ML, Jackson LA, Petrie JG, Ohmit SE, Monto AS, McLean HQ, Belongia EA, Fry AM, Flannery B. Influenza Vaccine Effectiveness Against 2009 Pandemic Influenza A(H1N1) Virus Differed by Vaccine Type During 2013-2014 in the United States. J Infect Dis. 2016; 213(10):1546–56. doi: 10.1093/infdis/jiv577. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiv577&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26743842&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 27. Gagnon A, Acosta E, Hallman S, Bourbeau R, Dillon LY, Ouellette N, Earn DJD, Herring DA, Inwood K, Madrenas J, Miller MS. Pandemic Paradox: Early Life H2N2 Pandemic Influenza Infection Enhanced Susceptibility to Death during the 2009 H1N1 Pandemic. mBio. 2018; 9(1):e02091–17. doi: 10.1128/mBio.02091-17. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.02091-17&link_type=DOI) 28. Gagnon A, Acosta E, Miller MS. Reporting and evaluating influenza virus surveillance data: An argument for incidence by single year of age. Vaccine. 2018; 36(42):6249–6252. doi: [https://doi.org/10.1016/j.vaccine.2018.08.077](https://doi.org/10.1016/j.vaccine.2018.08.077). 29. Gagnon A, Miller MS, Hallman SA, Bourbeau R, Herring DA, Earn DJD, Madrenas J. Age-Specific Mortality During the 1918 Influenza Pandemic: Unravelling the Mystery of High Young Adult Mortality. PLOS ONE. 2013; 8(8):e69586. doi: 10.1371/journal.pone.0069586. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0069586&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23940526&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 30. Goldstein E, Cobey S, Takahashi S, Miller JC, Lipsitch M. Predicting the Epidemic Sizes of Influenza A/H1N1, A/H3N2, and B: A Statistical Method. PLOS Medicine. 2011; 8(7):e1001051. doi: 10.1371/journal.pmed.1001051. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1001051&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21750666&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 31. Gostic KM, Ambrose M, Worobey M, Lloyd-Smith JO. Potent protection against H5N1 and H7N9 influenza via childhood hemagglutinin imprinting. Science. 2016; 354(6313):722–726. doi: 10.1126/science.aag1322. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTQvNjMxMy83MjIiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMy8yOS8xOTAwMTg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 32. Gostic KM, Bridge R, Brady S, Viboud C, Worobey M, Lloyd-Smith JO. Childhood immune imprinting to influenza A shapes birth year- specific risk during seasonal H1N1 and H3N2 epidemics. PLOS Pathogens. 2019 12; 15(12):1–20. doi: 10.1371/journal.ppat.1008109. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1007529&link_type=DOI) 33. Griffin MR, Monto AS, Belongia EA, Treanor JJ, Chen Q, Chen J, Talbot HK, Ohmit SE, Coleman LA, Lofthus G, Petrie JG, Meece JK, Hall CB, Williams JV, Gargiullo P, Berman L, Shay DK. Effectiveness of non-adjuvanted pandemic influenza A vaccines for preventing pandemic influenza acute respiratory illness visits in 4 U.S. communities. PLoS One. 2011; 6(8):e23085. doi: 10.1371/journal.pone.0023085. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0023085&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21857999&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 34. Groth SFdS, Webster R. Disquisitions on original antigenic sin: I. Evidence in man. Journal of Experimental Medicine. 1966; 124(3):331–345. [Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjk6IjEyNC8zLzMzMSI7czo0OiJhdG9tIjtzOjM5OiIvbWVkcnhpdi9lYXJseS8yMDIwLzAzLzI5LzE5MDAxODc1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 35. Huang KYA, Rijal P, Schimanski L, Powell TJ, Lin TY, McCauley JW, Daniels RS, Townsend AR. Focused antibody response to influenza linked to antigenic drift. The Journal of clinical investigation. 2015; 125(7):2631–2645. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI81104&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 36. Huang QS, Bandaranayake D, Wood T, Newbern EC, Seeds R, Ralston J, Waite B, Bissielo A, Prasad N, Todd A, Jelley L, Gunn W, McNicholas A, Metz T, Lawrence S, Collis E, Retter A, Wong SS, Webby R, Bocacao J, et al. Risk Factors and Attack Rates of Seasonal Influenza Infection: Results of the Southern Hemisphere Influenza and Vaccine Effectiveness Research and Surveillance (SHIVERS) Seroepidemiologic Cohort Study. J Infect Dis. 2019; 219(3):347–357. doi: 10.1093/infdis/jiy443. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiy443&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30016464&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 37. Irving SA, Donahue JG, Shay DK, Ellis-Coyle TL, Belongia EA. Evaluation of self-reported and registry-based influenza vaccination status in a Wisconsin cohort. Vaccine. 2009; 27(47):6546–9. doi: 10.1016/j.vaccine.2009.08.050. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2009.08.050&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19729083&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272056300008&link_type=ISI) 38. Jackson LA, Jackson ML, Nelson JC, Neuzil KM, Weiss NS. Evidence of bias in estimates of influenza vaccine effectiveness in seniors. International Journal of Epidemiology. 2005 12; 35(2):337–344. doi: 10.1093/ije/dyi274. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyi274&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16368725&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236817900028&link_type=ISI) 39. Jackson LA, Nelson JC, Benson P, Neuzil KM, Reid RJ, Psaty BM, Heckbert SR, Larson EB, Weiss NS. Functional status is a confounder of the association of influenza vaccine and risk of all cause mortality in seniors. International Journal of Epidemiology. 2005 12; 35(2):345–352. doi: 10.1093/ije/dyi275. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyi275&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16368724&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236817900029&link_type=ISI) 40. Jackson ML, Chung JR, Jackson LA, Phillips CH, Benoit J, Monto AS, Martin ET, Belongia EA, McLean HQ, Gaglani M, Murthy K, Zimmerman R, Nowalk MP, Fry AM, Flannery B. Influenza Vaccine Effectiveness in the United States during the 2015-2016 Season. N Engl J Med. 2017; 377(6):534–543. doi: 10.1056/NEJMoa1700153. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1700153&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28792867&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 41. Khiabanian H, Farrell GM, St George K, Rabadan R. Differences in patient age distribution between influenza A subtypes. PloS one. 2009; 4(8):e6832–e6832. doi: 10.1371/journal.pone.0006832. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0006832&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19718262&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 42. Kieke AL, Kieke BA, Kopitzke SL, McClure DL, Belongia EA, VanWormer JJ, Greenlee RT. Validation of Health Event Capture in the Marshfield Epidemiologic Study Area. Clinical Medicine & Research. 2015; 13(3-4):103–111. doi: 10.3121/cmr.2014.1246. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNsaW5tZWRyZXMiO3M6NToicmVzaWQiO3M6MTA6IjEzLzMtNC8xMDMiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMy8yOS8xOTAwMTg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 43. Laurie KL, Guarnaccia TA, Carolan LA, Yan AWC, Aban M, Petrie S, Cao P, Heffernan JM, McVernon J, Mosse J, Kelso A, McCaw JM, Barr IG. Interval Between Infections and Viral Hierarchy Are Determinants of Viral Interference Following Influenza Virus Infection in a Ferret Model. The Journal of infectious diseases. 2015; 212(11):1701–1710. doi: 10.1093/infdis/jiv260. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiv260&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25943206&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 44. Lee JKH, Lam GKL, Shin T, Kim J, Krishnan A, Greenberg DP, Chit A. Efficacy and effectiveness of high-dose versus standard-dose influenza vaccination for older adults: a systematic review and meta-analysis. Expert Review of Vaccines. 2018; 17(5):435–443. doi: 10.1080/14760584.2018.1471989. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/14760584.2018.1471989&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29715054&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 45. Lee JM, Eguia R, Zost SJ, Choudhary S, Wilson PC, Bedford T, Stevens-Ayers T, Boeckh M, Hurt AC, Lakdawala SS, Hensley SE, Bloom JD. Mapping person-to-person variation in viral mutations that escape polyclonal serum targeting influenza hemagglutinin. 2019 aug; 8:e49324. doi: 10.7554/eLife.49324. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.49324&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31452511&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 46. Lewnard JA, Tedijanto C, Cowling BJ, Lipsitch M. Measurement of Vaccine Direct Effects Under the Test-Negative Design. American Journal of Epidemiology. 2018 08; 187(12):2686–2697. doi: 10.1093/aje/kwy163. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwy163&link_type=DOI) 47. Linderman SL, Chambers BS, Zost SJ, Parkhouse K, Li Y, Herrmann C, Ellebedy AH, Carter DM, Andrews SF, Zheng NY, Huang M, Huang Y, Strauss D, Shaz BH, Hodinka RL, Reyes-Terán G, Ross TM, Wilson PC, Ahmed R, Bloom JD, et al. Potential antigenic explanation for atypical H1N1 infections among middle-aged adults during the 2013–2014 influenza season. Proceedings of the National Academy of Sciences. 2014; 111(44):15798. doi: 10.1073/pnas.1409171111. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTExLzQ0LzE1Nzk4IjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMjkvMTkwMDE4NzUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 48. Manicassamy B, Medina RA, Hai R, Tsibane T, Stertz S, Nistal-Villán E, Palese P, Basler CF, García-Sastre A. Protection of mice against lethal challenge with 2009 H1N1 influenza A virus by 1918-like and classical swine H1N1 based vaccines. PLoS pathogens. 2010; 6(1):e1000745. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1000745&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20126449&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 49. McLean HQ, Caspard H, Griffin MR, Gaglani M, Peters TR, Poehling KA, Ambrose CS, Belongia EA. Association of Prior Vaccination With Influenza Vaccine Effectiveness in Children Receiving Live Attenuated or Inactivated Vaccine. JAMA Network Open. 2018 10; 1(6):e183742–e183742. doi: 10.1001/jamanetworkopen.2018.3742. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2018.3742&link_type=DOI) 50. McLean HQ, Thompson MG, Sundaram ME, Meece JK, McClure DL, Friedrich TC, Belongia EA. Impact of repeated vaccination on vaccine effectiveness against influenza A (H3N2) and B during 8 seasons. Clinical Infectious Diseases. 2014; 59(10):1375–1385. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciu680&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25270645&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 51. Monto AS, Koopman JS, Longini J I M. Tecumseh study of illness. XIII. Influenza infection and disease, 1976-1981. Am J Epidemiol. 1985; 121(6):811–22. doi: 10.1093/oxfordjournals.aje.a114052. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/oxford-journals.aje.a114052&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=4014174&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1985AHV9600005&link_type=ISI) 52. Nakajima S, Nobusawa E, Nakajima K. Variation in Response among Individuals to Antigenic Sites on the HA Protein of Human Influenza Virus May Be Responsible for the Emergence of Drift Strains in the Human Population. Virology. 2000; 274(1):220–231. doi: [https://doi.org/10.1006/viro.2000.0453](https://doi.org/10.1006/viro.2000.0453). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/viro.2000.0453&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10936103&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000088957600023&link_type=ISI) 53. Neuzil KM, Jackson LA, Nelson J, Klimov A, Cox N, Bridges CB, Dunn J, DeStefano F, Shay D. Immunogenicity and Reactogenicity of 1 versus 2 Doses of Trivalent Inactivated Influenza Vaccine in Vaccine-Naive 5–8-Year-Old Children. The Journal of Infectious Diseases. 2006; 194(8):1032–1039. doi: 10.1086/507309. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/507309&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16991077&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000240548500003&link_type=ISI) 54. NHIS, Influenza Vaccination Coverage Trends 1989 – 2008. Accessed 11-March-2020; 2009. [https://www.cdc.gov/flu/pdf/professionals/nhis89\_08fluvaxtrendtab.pdf](https://www.cdc.gov/flu/pdf/professionals/nhis89_08fluvaxtrendtab.pdf). 55. O’Donnell CD, Wright A, Vogel LN, Wei CJ, Nabel GJ, Subbarao K. Effect of Priming with H1N1 Influenza Viruses of Variable Antigenic Distances on Challenge with 2009 Pandemic H1N1 Virus. Journal of Virology. 2012; 86(16):8625–8633. doi: 10.1128/jvi.00147-12. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjEwOiI4Ni8xNi84NjI1IjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMjkvMTkwMDE4NzUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 56. Ohmit SE, Thompson MG, Petrie JG, Thaker SN, Jackson ML, Belongia EA, Zimmerman RK, Gaglani M, Lamerato L, Spencer SM, Jackson L, Meece JK, Nowalk MP, Song J, Zervos M, Cheng PY, Rinaldo CR, Clipper L, Shay DK, Piedra P, et al. Influenza vaccine effectiveness in the 2011-2012 season: protection against each circulating virus and the effect of prior vaccination on estimates. Clin Infect Dis. 2014; 58(3):319–27. doi: 10.1093/cid/cit736. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cit736&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24235265&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 57. Ohmit SE, Petrie JG, Malosh RE, Johnson E, Truscon R, Aaron B, Martens C, Cheng C, Fry AM, Monto AS. Substantial Influenza Vaccine Effectiveness in Households With Children During the 2013–2014 Influenza Season, When 2009 Pandemic Influenza A(H1N1) Virus Predominated. The Journal of Infectious Diseases. 2015 11; 213(8):1229–1236. doi: 10.1093/infdis/jiv563. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiv563&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26597255&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 58. Pebody RG, Andrews N, Fleming DM, Mcmenamin J, Cottrell S, Smyth B, Durnall H, Robertson C, Carman W, Ellis J, et al. Age-specific vaccine effectiveness of seasonal 2010/2011 and pandemic influenza A(H1N1) 2009 vaccines in preventing influenza in the United Kingdom. Epidemiology and Infection. 2013; 141(3):620–630. doi: 10.1017/S0950268812001148. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268812001148&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000315374600025&link_type=ISI) 59. Petrie JG, Parkhouse K, Ohmit SE, Malosh RE, Monto AS, Hensley SE. Antibodies Against the Current Influenza A(H1N1) Vaccine Strain Do Not Protect Some Individuals From Infection With Contemporary Circulating Influenza A(H1N1) Virus Strains. The Journal of infectious diseases. 2016; 214(12):1947–1951. doi: 10.1093/infdis/jiw479. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiw479&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27923954&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 60. Ranjeva S, Subramanian R, Fang VJ, Leung GM, Ip DKM, Perera Rapm, Peiris JSM, Cowling BJ, Cobey S. Age-specific differences in the dynamics of protective immunity to influenza. Nature Communications. 2019; 10(1):1660. doi: 10.1038/s41467-019-09652-6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-09652-6&link_type=DOI) 61. Reed C, Chaves SS, Daily Kirley P, Emerson R, Aragon D, Hancock EB, Butler L, Baumbach J, Hollick G, Bennett NM, Laidler MR, Thomas A, Meltzer MI, Finelli L. Estimating influenza disease burden from population-based surveillance data in the United States. PLoS One. 2015; 10(3):e0118369. doi: 10.1371/journal.pone.0118369. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0118369&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25738736&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 62. Saito N, Komori K, Suzuki M, Kishikawa T, Yasaka T, Ariyoshi K. Dose-Dependent Negative Effects of Prior Multiple Vaccinations Against Influenza A and Influenza B Among Schoolchildren: A Study of Kamigoto Island in Japan During the 2011–2012, 2012–2013, and 2013–2014 Influenza Seasons. Clinical Infectious Diseases. 2018; 67(6):897–904. 63. Skowronski DM, Chambers C, De Serres G, Sabaiduc S, Winter AL, Dickinson JA, Gubbay JB, Fonseca K, Drews SJ, Charest H, Martineau C, Krajden M, Petric M, Bastien N, Li Y, Smith DJ. Serial Vaccination and the Antigenic Distance Hypothesis: Effects on Influenza Vaccine Effectiveness During A(H3N2) Epidemics in Canada, 2010-2011 to 2014-2015. J Infect Dis. 2017; 215(7):1059–1099. doi: 10.1093/infdis/jix074. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jix074&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28180277&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 64. Skowronski DM, Chambers C, Sabaiduc S, De Serres G, Winter AL, Dickinson JA, Gubbay JB, Drews SJ, Martineau C, Charest H, Krajden M, Bastien N, Li Y. Beyond Antigenic Match: Possible Agent-Host and Immuno-epidemiological Influences on Influenza Vaccine Effectiveness During the 2015-2016 Season in Canada. J Infect Dis. 2017; 216(12):1487–1500. doi: 10.1093/infdis/jix526. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jix526&link_type=DOI) 65. Skowronski DM, Chambers C, Sabaiduc S, De Serres G, Winter AL, Dickinson JA, Krajden M, Gubbay JB, Drews SJ, Martineau C, et al. A perfect storm: impact of genomic variation and serial vaccination on low influenza vaccine effectiveness during the 2014–2015 season. Clinical Infectious Diseases. 2016; 63(1):21–32. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciw176&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27025838&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 66. 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(5682):371–6. doi: 10.1126/science.1097211. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMDUvNTY4Mi8zNzEiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMy8yOS8xOTAwMTg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 67. Smith DJ, Forrest S, Ackley DH, Perelson AS. Variable efficacy of repeated annual influenza vaccination. Proceedings of the National Academy of Sciences. 1999; 96(24):14001–14006. doi: 10.1073/pnas.96.24.14001. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiOTYvMjQvMTQwMDEiO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMy8yOS8xOTAwMTg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 68. Smith NM, Bresee JS, Shay DK, Uyeki TM, Cox NJ, Strikas RA. Prevention and Control of Influenza: recommendations of the Advisory Committee on Immunization Practices (ACIP). MMWR Recomm Rep. 2006; 55(Rr-10):1–42. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16498385&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 69. Sullivan SG, Tchetgen Tchetgen EJ, Cowling BJ. Theoretical Basis of the Test-Negative Study Design for Assessment of Influenza Vaccine Effectiveness. American journal of epidemiology. 2016; 184(5):345–353. doi: 10.1093/aje/kww064. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kww064&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27587721&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 70. Thompson WW, Shay DK, Weintraub E, Brammer L, Cox N, Anderson LJ, Fukuda K. Mortality Associated With Influenza and Respiratory Syncytial Virus in the United States. JAMA. 2003; 289(2):179–186. doi: 10.1001/jama.289.2.179. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.289.2.179&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12517228&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180226400028&link_type=ISI) 71. Treanor JJ, Talbot HK, Ohmit SE, Coleman LA, Thompson MG, Cheng PY, Petrie JG, Lofthus G, Meece JK, Williams JV, Berman L, Breese Hall C, Monto AS, Griffin MR, Belongia E, Shay DK. Effectiveness of seasonal influenza vaccines in the United States during a season with circulation of all three vaccine strains. Clin Infect Dis. 2012; 55(7):951–9. doi: 10.1093/cid/cis574. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cis574&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22843783&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 72. Van Cauteren D, Vaux S, de Valk H, Le Strat Y, Vaillant V, Lévy-Bruhl D. Burden of influenza, healthcare seeking behaviour and hygiene measures during the A(H1N1)2009 pandemic in France: a population based study. BMC Public Health. 2012; 12(1):947. doi: 10.1186/1471-2458-12-947. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2458-12-947&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23127166&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 73. WDHS, Nursing Home Directory. Accessed 30-September-2019; 2018. [https://www.dhs.wisconsin.gov/guide/nursing-home.htm](https://www.dhs.wisconsin.gov/guide/nursing-home.htm). 74. Worby CJ, Chaves SS, Wallinga J, Lipsitch M, Finelli L, Goldstein E. On the relative role of different age groups in influenza epidemics. Epidemics. 2015; 13:10–16. doi: 10.1016/j.epidem.2015.04.003. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2015.04.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26097505&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) 75. Worobey M, Han GZ, Rambaut A. Genesis and pathogenesis of the 1918 pandemic H1N1 influenza A virus. Proceedings of the National Academy of Sciences of the United States of America. 2014; 111(22):8107–8112. doi: 10.1073/pnas.1324197111. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTExLzIyLzgxMDciO3M6NDoiYXRvbSI7czozOToiL21lZHJ4aXYvZWFybHkvMjAyMC8wMy8yOS8xOTAwMTg3NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 76. Wu JT, Ma ES, Lee CK, Chu DK, Ho PL, Shen AL, Ho A, Hung IF, Riley S, Ho LM, Lin CK, Tsang T, Lo SV, Lau YL, Leung GM, Cowling BJ, Malik Peiris JS. The infection attack rate and severity of 2009 pandemic H1N1 influenza in Hong Kong. Clin Infect Dis. 2010; 51(10):1184–91. doi: 10.1086/656740. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/656740&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20964521&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283331800012&link_type=ISI) 77. Wu S, L Vana, Wang L, McDonald SA, Pan Y, Duan W, Zhang L, Sun Y, Zhang Y, Zhang X, Pilot E, Krafft T, W Vdh, Mab VDS, Yang P, Wang Q. Estimated incidence and number of outpatient visits for seasonal influenza in 2015-2016 in Beijing, China. Epidemiology and infection. 2017; 145(16):3334–3344. doi: 10.1017/s0950268817002369. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/s0950268817002369&link_type=DOI) 78. Zimmerman RK, Nowalk MP, Chung J, Jackson ML, Jackson LA, Petrie JG, Monto AS, McLean HQ, Belongia EA, Gaglani M, Murthy K, Fry AM, Flannery B, Investigators ftUFV. 2014–2015 Influenza Vaccine Effectiveness in the United States by Vaccine Type. Clinical Infectious Diseases. 2016; 63(12):1564–1573. doi: 10.1093/cid/ciw635. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciw635&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27702768&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F29%2F19001875.atom) [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif [7]: /embed/graphic-9.gif [8]: /embed/graphic-10.gif [9]: /embed/graphic-11.gif [10]: /embed/inline-graphic-1.gif [11]: /embed/inline-graphic-2.gif [12]: /embed/graphic-12.gif [13]: /embed/graphic-13.gif [14]: /embed/graphic-14.gif [15]: /embed/graphic-15.gif [16]: /embed/graphic-16.gif [17]: /embed/graphic-17.gif [18]: /embed/inline-graphic-3.gif [19]: /embed/inline-graphic-4.gif [20]: /embed/graphic-18.gif [21]: /embed/graphic-19.gif [22]: /embed/graphic-20.gif [23]: /embed/graphic-21.gif [24]: /embed/graphic-22.gif [25]: /embed/graphic-23.gif [26]: /embed/graphic-24.gif [27]: /embed/graphic-25.gif [28]: /embed/inline-graphic-5.gif [29]: /embed/graphic-26.gif [30]: /embed/graphic-27.gif [31]: /embed/graphic-28.gif [32]: /embed/graphic-29.gif [33]: /embed/graphic-35.gif [34]: /embed/graphic-36.gif [35]: /embed/inline-graphic-6.gif [36]: /embed/graphic-37.gif [37]: /embed/inline-graphic-7.gif [38]: /embed/graphic-38.gif [39]: /embed/graphic-39.gif [40]: /embed/inline-graphic-8.gif [41]: /embed/inline-graphic-9.gif [42]: /embed/inline-graphic-10.gif [43]: /embed/inline-graphic-11.gif [44]: /embed/graphic-40.gif [45]: /embed/graphic-41.gif [46]: /embed/inline-graphic-12.gif [47]: /embed/inline-graphic-13.gif [48]: /embed/inline-graphic-14.gif [49]: /embed/inline-graphic-15.gif [50]: /embed/inline-graphic-16.gif [51]: /embed/inline-graphic-17.gif [52]: /embed/graphic-42.gif [53]: /embed/graphic-43.gif [54]: /embed/graphic-44.gif [55]: F12/embed/inline-graphic-18.gif