Projections of the impact of transmission-reducing pediatric influenza vaccine in South Africa using high resolution immunologic and infection data

preprint OA: closed
Full text 83,403 characters Β· extracted from preprint-html Β· click to expand
Projections of the impact of transmission-reducing pediatric influenza vaccine in South Africa using high resolution immunologic and infection data | medRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-P4HH5NV'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search Projections of the impact of transmission-reducing pediatric influenza vaccine in South Africa using high resolution immunologic and infection data View ORCID Profile Jackie Kleynhans , CΓ©cile Viboud , Molly Sauter , Stefano Tempia , Nicole Wolter , Jocelyn Moyes , Anne von Gottberg , Lorens Maake , Alexandra Moerdyk , Cheryl Cohen , Kaiyuan Sun doi: https://doi.org/10.1101/2025.04.04.25325228 Jackie Kleynhans 1 Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service , Johannesburg, South Africa 2 School of Public Health, Faculty of Health Sciences, University of the Witwatersrand , Johannesburg, South Africa Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jackie Kleynhans For correspondence: jackiel{at}nicd.ac.za kaiyuan.sun{at}nih.gov CΓ©cile Viboud 3 Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health , Bethesda, Maryland, United States of America Find this author on Google Scholar Find this author on PubMed Search for this author on this site Molly Sauter 3 Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health , Bethesda, Maryland, United States of America Find this author on Google Scholar Find this author on PubMed Search for this author on this site Stefano Tempia 1 Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service , Johannesburg, South Africa 2 School of Public Health, Faculty of Health Sciences, University of the Witwatersrand , Johannesburg, South Africa Find this author on Google Scholar Find this author on PubMed Search for this author on this site Nicole Wolter 1 Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service , Johannesburg, South Africa 4 School of Pathology, Faculty of Health Sciences, University of the Witwatersrand , Johannesburg, South Africa Find this author on Google Scholar Find this author on PubMed Search for this author on this site Jocelyn Moyes 1 Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service , Johannesburg, South Africa Find this author on Google Scholar Find this author on PubMed Search for this author on this site Anne von Gottberg 1 Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service , Johannesburg, South Africa 4 School of Pathology, Faculty of Health Sciences, University of the Witwatersrand , Johannesburg, South Africa Find this author on Google Scholar Find this author on PubMed Search for this author on this site Lorens Maake 1 Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service , Johannesburg, South Africa Find this author on Google Scholar Find this author on PubMed Search for this author on this site Alexandra Moerdyk 1 Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service , Johannesburg, South Africa Find this author on Google Scholar Find this author on PubMed Search for this author on this site Cheryl Cohen 1 Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service , Johannesburg, South Africa 2 School of Public Health, Faculty of Health Sciences, University of the Witwatersrand , Johannesburg, South Africa Find this author on Google Scholar Find this author on PubMed Search for this author on this site Kaiyuan Sun 3 Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health , Bethesda, Maryland, United States of America Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: jackiel{at}nicd.ac.za kaiyuan.sun{at}nih.gov Abstract Full Text Info/History Metrics Data/Code Preview PDF Abstract Background Traditional influenza vaccination programs focus on preventing severe illness through the vaccine’s direct effect. However, vaccinating children may substantially reduce overall burden by lowering transmission, especially in high infection settings and when effective transmission reducing vaccines become available. Methods To evaluate the impact of a transmission-reducing vaccination campaign, we developed a Susceptible-Latent-Infectious-Recovered (SLIR) compartmental transmission model calibrated with high-resolution virological and immunological data sourced from a multiyear household cohort study in urban and rural South Africa. The mode explicitly accounted for variations in susceptibility, infectiousness, and viral shedding duration by age group, pre-season hemagglutination inhibition antibody titer and influenza subtype/lineage. Using this framework, we simulated pediatric vaccination scenarios targeting children aged 6 months – 5 years or 6-12 years, varying both vaccine coverage and the vaccine’s effectiveness in reducing susceptibility and transmission. Scenarios were compared to no-vaccination counterfactuals to project reductions in infections, illness episodes, hospitalizations, and deaths in South Africa. Results Our model projects that vaccinating 50% of children aged 6 months - 5 years or 6-12 years with a vaccine that is 50% effective in reducing both susceptibility and viral shedding could lower the overall infection incidence by 58-60% and 48%-51%, respectively. Illness episodes would be reduced by 20-60%, hospitalizations by 26-64%, and deaths by 21-59% across all age groups. The target age group among children that resulted in the highest population-wide vaccination benefits varied with the circulating lineage, season, and urbanization level. Conclusions Our work demonstrates that pediatric immunization with transmission-reducing vaccines holds significant potential for lowering community influenza burden, particularly in high transmission settings. Our modeling framework offers a practical translational tool for anticipating the impact of new vaccines in low and middle-income settings. Introduction Influenza remains a significant cause of morbidity and mortality globally despite well-established influenza vaccine programs. Influenza-associated lower respiratory tract infections are globally responsible for 124 hospitalizations and 1.9 deaths per 100,000 population annually, with higher rates (156 hospitalizations and 2.7 deaths per 100,000 population) in sub-Saharan Africa [ 1 ]. In South Africa, where the influenza burden pyramid has been well-characterized, an estimated 10.7 million episodes of influenza-associated illness occur annually, 1.2% of these are severe and 0.1% result in death [ 2 ], amounting to an annual cost of US$270.5 million [ 3 ]. While children younger than 13 years in this setting have very high attack rates (51-66 per 100 person seasons) [ 4 ], severe illness and mortality is concentrated at the extremes of age, and among people living with HIV and other underlying conditions [ 5 ]. Seasonal influenza vaccination reduces the risk of severe influenza-associated illness and death in vaccinated individuals [ 6 ]. Globally, vaccination programs have focused on high risk groups, including pregnant woman, children aged 6-59 months, the elderly and those with comorbidities [ 7 ]. Although forty-five percent of South Africa’s population falls into high-risk groups targeted for influenza vaccination [ 8 ], the public national influenza program procures only a limited number of vaccines, prioritizing individuals in these groups and healthcare workers. Yet, seasonal influenza coverage in South Africa remains low, estimated at 4.6% of this target population in 2018 [ 9 ]. An alternative influenza control strategy could be considered, which would harness the transmission-reducing effects of vaccines. By prioritizing vaccination in age groups with high attack rates and significant roles in viral transmission, this approach would aim to reduce influenza circulation at the community level. High-risk individuals who remain unvaccinated benefit indirectly from this policy through reduced exposure to influenza – the so-called indirect effect. To date, however, this transmission-focused strategy has been primarily explored through mathematical modeling [ 10 – 13 ], with real-world assessments mostly focused in high income countries such as Canada, Japan, Russia, the United States of America and the United Kingdom [ 14 – 20 ]. While a school-based influenza vaccination program has been implemented in England [ 19 , 20 ], this approach has not yet been implemented in any low and middle-income countries, where influenza burden is the highest [ 1 ], and potential indirect impacts may be greater. Two major challenges impede the practical implementation of influenza immunization programs designed to harness the vaccine’s indirect protective effects. First, current vaccines are primarily designed to prevent severe outcomes, and their effects on transmission are rarely assessed in clinical trials. Second, vaccination strategies relying on indirect vaccine effects should target the population groups that contribute the most to influenza transmission, regardless of clinical severity. Although children have been identified as major transmitters by a number of symptom-based surveillance and epidemiological studies [ 4 , 21 , 22 ], these studies often overlook asymptomatic and subclinical infections that also drive transmission. As a result, the potential benefits of pediatric vaccination may be under- or overestimated depending on the age group in which much asymptomatic infections occur, and the findings may not generalize to regions with different population demographics and influenza dynamics. In this study, we introduce a comprehensive modeling framework to assess the potential impact of transmission-reducing vaccines targeting children. By leveraging detailed infection, serological and contact data from South African household cohorts, our model integrates the roles of asymptomatic transmission, variations in susceptibility and infectiousness by age and pre-existing HAI titers, and mixing patterns [ 4 , 23 , 24 ]. By systematically simulating a range of vaccine efficacy profiles and immunization coverages against different influenza subtypes/strains, we identify strategies that could markedly reduce the overall burden of influenza by targeting the high-transmitting pediatric population. Our approach not only addresses critical gaps in the current assessments of vaccine impact on transmission but also provides actionable tools for optimizing vaccination strategies. Materials and Methods To evaluate the potential impact of transmission-reducing influenza vaccines, we developed a transmission and vaccine-impact modeling framework calibrated to detailed infection and serological data from longitudinal household cohorts in South Africa. First, we describe the underlying cohort data. Next, we introduce the structure, parameterization and calibration of our transmission model, along with assumptions underpinning our vaccine intervention simulations. Finally, we outline our methodology to project the amount of disease burden averted by implementing a pediatric vaccination program in South Africa. Household cohort study for model calibration We used data from a prospective household cohort study where 50 households were randomly selected annually from each of an urban and a rural community and intensively followed up for six to ten months spanning the 2016-2018 influenza seasons in South Africa (PHIRST study [ 4 ]). Briefly, participants were visited twice a week for the collection of nasopharyngeal specimens for PCR testing of influenza viruses A and B irrespective of symptoms, with typing to A(H1N1)pdm09, A(H3N2), B/Victoria and B/Yamagata for all positive samples. In addition, pre-season serum was collected to test for presence of influenza antibodies for each subtype/lineage, using hemagglutination inhibition (HAI) assays with viruses that circulated during the period (A/South Africa/2517/2016 A(H1N1)pdm09, A/Singapore/INFIMH-16-0019/2016 A(H3N2), B/South Africa/R3037/2016 (B/Victoria), and B/South Africa/R05631/2017 (B/Yamagata)). In 2018, a contact survey was performed to collect age-based contact rates in the cohort [ 25 ]. Model We constructed a deterministic Susceptible-Latent-Infectious-Recovered (SLIR) transmission model, where individuals move between susceptible (S), latent (L), infectious (I) and recovered (R) compartments. The model allowed for both susceptibility to infection and infectiousness to vary by age and influenza antibody titers (hereafter referred to as titers), with age- and titer-specific susceptibility and infectiousness calibrated on our prior analysis of shedding dynamics and infection risk in the South African household cohort [ 24 ]. The compartmental model was further stratified by age (<6 months, 6 months – 5 years, 6-12 years, 13-18 years, 19-44 years, 45-64 years and β‰₯65 years), and pre-season titer for each subtype/lineage (<1:40, 1:40-1:80, β‰₯1:160). The impact of age and antibody titers on transmissibility was modeled by assigning different latent periods, infectious periods, and relative infectiousness depending on age and titer, following prior analysis ([ 24 ] and Supplementary Excel File 1). Within each age and titer group, the average latent period was defined as the sum of half the duration of the viral shedding proliferation phase and the interval from infection to the first detectable PCR positivity [ 26 ]; the average infectious period was half of the shedding duration; and the relative infectiousness was approximated by the ratios of peak viral loads ( Fig. S1 ). Compared to the original estimates in Sauter et al. [ 24 ], we updated the age categories and discretized titer as a categorical variable to accommodate the parameterization of the SLIR model used in this paper. Adjusted susceptibility estimates are provided in Supplementary Excel File 1. Download figure Open in new tab Fig. S1: Definition of latent period, infectious period and infectiousness. Estimated within age and pre-season HAI titer. Separate models were built for each cohort study site (rural: Agincourt, urban: Klerksdorp), season (2016, 2017, 2018) and subtype/lineage (A(H1N1)pdm09, A(H3N2), B/Victoria and B/Yamagata). In some of our analyses, we projected influenza epidemics separately for each subtype/lineage, and in others we combined projections for all subtypes/lineages co-circulating the same year to reflect the β€œtotal” annual influenza burden (by summing infections). Population size ( N ) was set at 10,000,000, with the age distribution reflecting the age distribution of the PHIRST cohort for the respective sites ( Figures 1A ). The site-specific mixing patterns ( Figures 1B ) among different age groups were based on a cross-sectional contact survey done in the 2018 cohort [ 25 ]. Within each age group, the proportion of individuals with each antibody titer ( h ) was based on the cohort site- and year-specific pre-season titers for the respective subtype/lineage ( Figures 1C ). Download figure Open in new tab Fig. 1. Model inputs (A-G) and calibration targets (H-I) for A(H1N1)pdm09, urban site, 2018. A) Model population age structure (m – months, y - years). B) Cumulative contact rates between age groups at home, school, work and within the community, weighted by calibrated omega (matrix weight). C) Proportion of population with pre-season hemagglutination inhibition assay (HAI) titers for A(H1N1)pdm09 within each age group. D) Susceptibility (Ξ±, risk of infection), E) infectiousness (Ξ², minimum Ct value as proxy for viral load), F) infectious period (Β΅, Β½ duration of shedding) and G) latent period (Οƒ, infection to infectiousness) by age and pre-season HAI titers. H) Observed and modelled incidence per 100 population (with 95% confidence intervals) by age group and I) weekly observed and modelled A(H1N1)pdm09 incidence used for model calibration. We did not consider antibody (immunity) waning during the season, or cross-protection of antibodies across subtypes/lineages. To project the potential impact of vaccination, the model was further broken down by vaccination status (vaccinated or not at the start of the season), in addition to age and antibody titers. This allowed us to model the potential impact of vaccination on susceptibility against infection as well as shedding kinetics and infectiousness. The equations of this SLIR model were: Where 𝑁 π‘”β„Žπ‘– is the number of individuals with vaccination status g , pre-season titer h in age group represent the number of individuals in the Susceptible, Latent, Infectious and Recovered compartments at time point t. are the rate at which individuals move between compartments. 𝛼 π‘”β„Žπ‘— is the susceptibility, 𝛽 π‘”β„Žπ‘— is the infectiousness, is the daily contact rate between individuals in age group i with age group j in contact setting k (e.g. household, school, workplace, or community). πœ” π‘˜ is the relative weight of contacts in setting k , as per-contact transmission risk is setting dependent. 𝛾 𝑑 is the baseline transmissibility at time point t . 𝜎 π‘”β„Žπ‘– is the 1/latency period, πœ‡ π‘”β„Žπ‘– the recovery rate (1/duration of shedding). All individuals in 𝑅 compartment behave the same, irrespective of their vaccine, pre-season titer level, and age, with no transition out of 𝑅 compartment. 𝑉𝐸 𝑠 and 𝑉𝐸 𝑖 represent the vaccine effectiveness on susceptibility and infectiousness (viral load), respectively. 𝑉𝐸 𝑖 and 𝑉𝐸 𝑑 collectively is the vaccine effectiveness against transmission. In the current models, 𝑉𝐸 𝑖 was set to 0 and 𝑉𝐸 𝑑 is considered the VE against transmission, although the framework would allow 𝑉𝐸 𝑖 to be adjusted. To initialize the model, we set the initial number of individuals infected , where is the proportion of individuals initially infected (the same across stratification) at the beginning of each influenza season, which is a free parameter to be estimated during model calibration. The initial number of susceptible was , and the initial number in the latent and recovered compartments were 0. 𝑑 0 represents the starting time of the epidemic season. A description of model parameters is provided in Table S1 and actual parameter values used for individual models are provided in Supplementary Excel file 1. Calibrated model parameters included the initial proportion infected , the time-varying transmissibility (𝛾 𝑑 , unique to each strain and year), and a contact matrix weight (πœ”, shared across strains and years). All of the model parameters were based on estimates directly derived from the PHIRST cohort, as previously described in Sauter et al. [ 24 ] and Kleynhans et al. [ 25 ]. View this table: View inline View popup Download powerpoint Table S1: Parameters used in SLIR transmission and vaccination model. Model calibration Models were calibrated to match site-, season-, and subtype-specific infection data from PHIRST based on PCR-confirmed influenza. We calibrated the model against 1) weekly all-age incidence of infections (where infection was assessed by twice weekly reverse transcription real-time polymerase chain reaction (rRT-PCR) tests on nasopharyngeal swabs irrespective of symptoms [ 4 ]) and 2) age-specific cumulative seasonal infection attack rate ( Table S2 ). View this table: View inline View popup Download powerpoint Table S2: Influenza incidence by age 2016-2018, urban and rural site, South Africa. NC (not considered) demotes subtype/lineage with <2/100 population incidence overall and not considered in modeling framework. For each subtype/lineage model for each site and year, we calculated the squared residual of the weekly incidence (π‘…π‘’π‘ π‘–π‘‘π‘’π‘Žπ‘™ 2 π‘€π‘˜_𝑖𝑛𝑐 ) and age-based and overall attack rates (π‘…π‘’π‘ π‘–π‘‘π‘’π‘Žπ‘™ 2 π‘Žπ‘”π‘’_π‘œπ‘£π‘’π‘Ÿπ‘Žπ‘™π‘™_𝑖𝑛𝑐 ) between the model and observed (PHIRST) data. Model error π‘…π‘’π‘ π‘–π‘‘π‘’π‘Žπ‘™ 2 subtype/lineage was calculated as the product of π‘…π‘’π‘ π‘–π‘‘π‘’π‘Žπ‘™ 2 π‘€π‘˜_𝑖𝑛𝑐 and π‘…π‘’π‘ π‘–π‘‘π‘’π‘Žπ‘™ 2 π‘Žπ‘”π‘’_π‘œπ‘£π‘’π‘Ÿπ‘Žπ‘™π‘™_𝑖𝑛𝑐 . We then calculated an overall observational error as the product of π‘…π‘’π‘ π‘–π‘‘π‘’π‘Žπ‘™ 2 subtype/lineage across all lineages and years, which we minimized. We only modelled and considered subtypes/lineages with an attack rate above 2/100 population per season in the combined models ( Table S2 ). For calibration purposes, we set vaccine coverage at 0 for all ages since vaccination was not routinely used in the cohort study population (<1% [ 4 ]). Modeling vaccine impact on infections To assess the potential impact of different coverage levels, age groups targeted, and vaccine profiles, we ran year-, site- and subtype/lineage-specific models for different vaccination scenarios, exploring different levels of vaccine effectiveness for 𝑉𝐸 𝑠 and 𝑉𝐸 𝑑 and vaccine coverage in 6 months – 5 years, 6-12 years or combined 6 months – 12 years. The calibrated model, with vaccine coverage set to 0, represented the counterfactual no-vaccination base model. Vaccine impact was calculated as the absolute (𝐼𝑛𝑐𝑖𝑑𝑒𝑛𝑐𝑒 π‘π‘Žπ‘ π‘’ π‘šπ‘œπ‘‘π‘’π‘™ βˆ’ 𝐼𝑛𝑐𝑖𝑑𝑒𝑛𝑐𝑒 π‘‰π‘Žπ‘π‘π‘–π‘›π‘’ π‘šπ‘œπ‘‘π‘’π‘™ ) and percentage reduction in incidence by age group and overall between the model with (vaccine model) and without (base model) vaccination. We also calculated infections averted per dose as To calculate the overall influenza incidence and vaccine impact, we created β€˜typical’ influenza seasons, through combining the infections from the respective subtypes/lineages reported to co-circulate in a given site and year and divided by the population size. These seasonal incidence estimates were used in the vaccine impact calculations. Projections of vaccine reductions in disease burden To project the overall impact of different vaccination strategies on influenza burden, we used age-based multipliers developed for the South African population to derive the number of (symptomatic) illness episodes, hospitalizations and deaths based on projected infections by our model ( Table S3 ). To estimate the number of symptomatic illness episodes, we used the symptomatic fraction reported in the PHIRST cohort [ 4 ]. To estimate hospitalizations and deaths, we used the proportion of infections resulting in hospitalizations or deaths from a disease burden study in South Africa [ 2 ]. We adapted the age breakdowns of the disease burden study [ 2 ] to our transmission model ( Table S3 ). View this table: View inline View popup Download powerpoint Table S3: Age-specific proportion of infectious resulting in illness episodes [ 4 ], severe episodes or death [ 2 ]. To account for inter-annual fluctuations in influenza disease burden, we averaged illness outcomes over the three study years to obtain annual rates. This was done for the counterfactual and for each vaccine coverage scenario, with a VE set at 50% against both susceptibility and duration of shedding for illustration. We assumed that the vaccine had no additional impact on illness severity beyond the effect on preventing infection and reducing transmission. In other words, we ignored the potential effect of the vaccine on reducing the risk of developing severe illness among breakthrough infections. The percentage reduction of each vaccine scenario and illness outcome was calculated as , with illness outcome representing symptomatic illness, hospitalizations or deaths. Model simulation and analyses were performed using Python v3.9.12. Results In Fig. 1 we summarize the model inputs and calibration targets for the model. The age distribution in the rural and urban sites differed in the South African household cohorts, with a higher proportion of children aged 6 months – 18 years at the rural site ( Fig. S2 A-B ). These differences likely drove the differences in contact patterns between the two sites, with much higher contact rates in the rural site among children and teenagers ( Fig. S2 C-J ). Pre-season antibody titers differed across age groups, years and subtypes/lineages ( Fig. S3 - 4 ). Susceptibility ( Fig. 1D , S2), infectiousness ( Fig. 1E , S6), infectious period ( Fig. 1F , S7) and latent period ( Fig. 1G , S8) estimates followed a decreasing gradient with the increase in age and titer. These factors likely contributed to the year-to-year and site variations in subtypes/lineages circulating, with differences in incidence rates observed by season and site ( Fig. 1H-I , S9-10). Our models performed well to estimate age-based incidence within the 95% confidence intervals of observed data, and to track weekly epidemic curves each season ( Fig. 1H-I , S9-10). Download figure Open in new tab Fig. S2: Age distribution at the A) rural and B) urban sites. Average number of contacts in the household (C, G), school (D, H), work (E, I) and community (F, J) at the rural (C-F) and urban (G-J) sites, 2018, South Africa Download figure Open in new tab Fig. S3: Proportion of population with pre-season hemagglutination inhibition assay (HAI) titers for A(H1N1)pdm09, A(H3N2), B/Victoria and B/Yamagata within each age group and season, rural site, 2016-2018, South Africa Download figure Open in new tab Fig. S4: Proportion of population with pre-season hemagglutination inhibition assay (HAI) titers for A(H1N1)pdm09, A(H3N2), B/Victoria and B/Yamagata within each age group and season, urban site, 2016-2018, South Africa In Fig. 2 , we illustrate the impact of vaccine modeling results for 2018 at the urban site with 50% of children aged 6-12 years vaccinated, with additional results for other vaccine scenarios, sites and years shown in Fig. S9 - 26 . Our model predicted that reductions in influenza incidence varied according to vaccine impact on susceptibility and viral shedding, with both mechanisms contributing comparably ( Fig. 2A ). For an exemplar season marked by co-circulation of A(H1N1)pdm09 and B/Victoria (as in 2018), we observed a 42% reduction in influenza incidence in all ages with a VE of 50% against duration of shedding and no effect on susceptibility, compared to a 53% reduction with 50% VE against susceptibility and no effect on shedding. Assuming a 50% VE against both duration of shedding and susceptibility resulted in a 67% reduction in influenza incidence. We used this VE assumption for the remaining analyses in the paper. Download figure Open in new tab Fig. S5: Susceptibility (Ξ±, risk of infection) by age and pre-season hemagglutination inhibition assay (HAI) titers for A(H1N1)pdm09, A(H3N2), B/Victoria and B/Yamagata, both sites, 2016-2018, South Africa Download figure Open in new tab Fig. S6: Infectiousness (Ξ², minimum Ct value as proxy for viral load) by age and pre-season hemagglutination inhibition assay (HAI) titers for A(H1N1)pdm09, A(H3N2), B/Victoria and B/Yamagata, both sites, 2016-2018, South Africa Download figure Open in new tab Fig. S7: Infectious period (Β΅, Β½ duration of shedding) by age and pre-season hemagglutination inhibition assay (HAI) titers for A(H1N1)pdm09, A(H3N2), B/Victoria and B/Yamagata, both sites, 2016-2018, South Africa Download figure Open in new tab Fig. S8: Latent period (Οƒ, infection to infectiousness) by age and pre-season hemagglutination inhibition assay (HAI) titers for A(H1N1)pdm09, A(H3N2), B/Victoria and B/Yamagata, both sites, 2016-2018, South Africa Download figure Open in new tab Fig. S9: Observed and modelled incidence per 100 population (with 95% confidence intervals) by age group and weekly observed and modelled incidence used for model calibration for each season (2016 (A-D), 2017 (E-H), 2018 (I-L) and subtype/lineage, rural site, South Africa. Download figure Open in new tab Fig. 2. Vaccine model results for vaccinating children aged 6-12 years, urban site, 2018 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. The weekly incidence in vaccination scenarios initially followed the same trajectory as the no-vaccination scenario, diverging primarily at the peakβ€”where the magnitude was reduced, though its timing remained consistent ( Fig. 2B ). The overall percentage reduction depended on age groups targeted for vaccination, with greatest benefits observed when the age group with highest incidence was immunized. Notably, the age group with highest incidence varied based on the circulating subtypes/lineages, year and site. For example, during the 2018 urban season, children aged 6-12 years experienced an attack rate of 62 per 100 persons, and vaccinating this age group reduced the overall incidence by 67% ( Fig. 2C ). In contrast, the urban 2017 season, characterized by circulation of A(H3N2) and B/Yamagata, recorded extremely high attack rates in children aged 13-18 years (98 per 100 persons, which includes multiple infections), 2 to 3 times higher than in children aged 6 months - 12 years or 6-12 years. In this season, we find more modest incidence reduction in all ages when vaccinating children aged 6 months – 5 years or 6 – 12 years, at 9% and 34% respectively ( Fig. S14 and S15 ). As expected, higher vaccine coverage led to greater reductions in incidence, with a rapid decline in incidence as coverage increased up to 50%, followed by more gradual benefits beyond that level ( Fig. 2D ). Additional model results for all seasons, vaccination strategies and both sites are presented in S11-S27. Download figure Open in new tab Fig. S10: Observed and modelled incidence per 100 population (with 95% confidence intervals) by age group and weekly observed and modelled incidence per 100 population used for model calibration for each season (2016 (A-D), 2017 (E-H), 2018 (I-L) and subtype/lineage, urban site, South Africa. Download figure Open in new tab Fig. S11: Vaccine model results for vaccinating children aged 6 months - 5 years, urban site, 2016 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S12: Vaccine model results for vaccinating children aged 6 - 12 years, urban site, 2016 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S13: Vaccine model results for vaccinating children aged 6 months - 12 years, urban site, 2016 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S14: Vaccine model results for vaccinating children aged 6 months - 5 years, urban site, 2017 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S15: Vaccine model results for vaccinating children aged 6 - 12 years, urban site, 2017 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S16: Vaccine model results for vaccinating children aged 6 months - 12 years, urban site, 2017 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S17: Vaccine model results for vaccinating children aged 6 months - 5 years, urban site, 2018 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S18: Vaccine model results for vaccinating children aged 6 - 12 years, urban site, 2018 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S19: Vaccine model results for vaccinating children aged 6 months - 5 years, rural site, 2016 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S20: Vaccine model results for vaccinating children aged 6 - 12 years, rural site, 2016 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S21: Vaccine model results for vaccinating children aged 6 months - 12 years, rural site, 2016 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S22: Vaccine model results for vaccinating children aged 6 months - 5 years, rural site, 2017 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S23: Vaccine model results for vaccinating children aged 6 - 12 years, rural site, 2017 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S24: Vaccine model results for vaccinating children aged 6 months - 12 years, rural site, 2017 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S25: Vaccine model results for vaccinating children aged 6 months - 5 years, rural site, 2018 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S26: Vaccine model results for vaccinating children aged 6 - 12 years, rural site, 2018 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. Download figure Open in new tab Fig. S27: Vaccine model results for vaccinating children aged 6 months - 12 years, rural site, 2018 season. A) Percentage reduction in incidence at differing levels of vaccine effectiveness (VE) against susceptibility (x axis) and duration of shedding (y axis) when vaccinating 50% of the target age group. Red block denotes VE chosen (50%) for subsequent analysis. B) Modelled weekly incidence per 100 population of individual influenza subtype/lineage in base model (no vaccination) and vaccine model. C) Overall incidence of influenza per 100 population subtype/lineage by age in base and vaccine model. Both B and C based on 50% coverage in target age group with vaccine with VE of 50% against susceptibility and duration of shedding. D) Incidence of influenza per 100 population subtype/lineage (left and middle) and combined season (right) at varying vaccine coverage levels in the target age group with vaccine with VE of 50% against susceptibility and duration of shedding. To better understand the potential range of vaccine benefits, we summarized infection reductions by targeted age group, year and site ( Fig. 3 ). We assumed 50% vaccine coverage using a vaccine with 50% VE in reducing both susceptibility to infection and the duration of viral shedding. Projected vaccine benefits varied substantially by year and study site, ranging from a 9% reduction in incidence in all ages when vaccinating children aged 6 months – 5 years in the urban site in 2017, to a 93% reduction when vaccinating children aged 6 months – 12 years in the 2018 season in the rural site ( Fig. 3 ). This correspond to a range of 0.5 to 2 infections averted for every dose administered. When only one of the two age groups was targeted, we saw the greatest percentage reduction when vaccinating children aged 6-12 years in the rural site in the 2018 season when A(H1N1)pdm09 and B/Victoria were circulating, with a 69% reduction in incidence and 2.5 infections averted per dose administered ( Fig. 3B ). In both sites, we observed the lowest reductions in incidence when A(H3N2) and B/Yamagata were circulating, and highest when A(H1N1)pdm09 and B/Victoria were circulating. There were also differences in the level of reduction between the two sites even when the same subtypes/lineages were circulating ( Fig. 3A, B ). For example, in 2018, although both sites experienced A(H1N1)pdm09 and B/Victoria, the urban site benefited more from vaccinating children aged 6–12 years, whereas the rural site saw greater effects from vaccinating children aged 6 months – 5 years. Download figure Open in new tab Fig. 3. Modelled influenza incidence per 100 population (bars, primary axis) in the base model (no vaccination) and three scenarios of vaccinating 50% of the target age group (6 months to 5 years, 6-12 years and 6 months – 12 years) with a vaccine with 50% effectiveness against susceptibility and duration of shedding; infections averted per dose (dots, secondary axis); reference for one infection averted per dose (grey dotted line, the secondary axis); percentage reduction (text at edge of bar), at A) urban and B) rural site, 2016-2018 seasons, South Africa Given the heterogeneity of influenza circulation and vaccine effects across different seasons, we summarized the putative benefits of a vaccination program on infections, symptomatic illness, hospitalizations and deaths by averaging results over 3 seasons ( Fig. 4 ). Across all illness outcomes and for both metrics of percentage reduction and burden averted per dose administered, the strategy to vaccinate children aged 6 months – 5 years fared better in the rural site, while the strategy to vaccinate children aged 6-12 years fared better in the urban site ( Fig. 4B ). Depending on the outcome considered, the influenza disease burden could be reduced by 58-64% in the rural site, and 43-51% in the urban site. For both sites, irrespective of strategy, at least 1.2 and 0.6 infections and illness episodes could be averted per dose, corresponding to at least 8 hospitalizations and 0.4 deaths averted per 1000 doses administered. Download figure Open in new tab Fig. 4. Annual average projected burden of influenza and reductions due to vaccinations. A) infections, B) illness episodes (symptomatic illness), C) hospitalizations and D) deaths per 100 person seasons (bars, primary axis). Results are shown for the base model (no vaccination) and three scenarios of vaccinating 50% of the target age group (6 months to 5 years, 6-12 years and 6 months – 12 years) with a vaccine with 50% effectiveness against infection and duration of shedding. The graph displays annual event rates (left y-axis, infections, illness episodes, hospitalizations or deaths), and the number of averted events per vaccine dose (dots, secondary axis). Percentage reduction is also shown (text at edge of bar). Data are based on model projections for urban and rural sites, averaged over the 2016-2018 seasons, South Africa Discussion By modeling high-resolution data on infection, serology, and contacts from a South African household cohort study, we evaluated the potential benefits of pediatric influenza vaccination programs in mitigating community-wide transmission and disease burden. Our modeling analysis indicates that even with moderate vaccination coverage (50%) and modest vaccine efficacy against both infection and transmission (50% each, respectively), a transmission-reducing immunization campaign targeting children could yield substantial benefits. Specifically, such a strategy may reduce infections, illness episodes, hospitalizations, and deaths by 21–87% in an average season, depending on circulating strain, geographical settings, and targeted age groups ( Figure 4 ). Even without accounting for the effects of vaccine strain match and the timing of influenza vaccination relative to the influenza season, we observed considerable differences in the overall reduction of influenza incidence across seasons and sites. These variations were likely driven by a combination of factors, including the circulating subtype/lineage and age-based attack rates for the specific season, as well as the population age structure and contact patterns. Here we only considered how variability in influenza attack rates may affect the population-level benefits of vaccination, but additional variation in VE have been reported between seasons and populations due to vaccine mismatch and immunologic factors [ 27 ]. Our data highlights how the optimal pediatric influenza vaccine program may vary across seasons and sites, driven by different mean ages at infection by season and geography. Our cohort study in South Africa revealed very high attack rates of influenza infection in the population (44 infections per 100 person seasons), with rates as high as 67 infections per 100 person seasons in children younger than 5 years [ 4 ]. Children also shed virus for longer and at higher viral loads [ 4 , 28 – 30 ], which translates to higher transmission potential [ 4 , 21 , 30 – 32 ]. Taken together, these factors make a pediatric influenza vaccine program particularly attractive for South Africa, as demonstrated by our model projections. Previous modeling studies have similarly demonstrated promising direct and indirect benefits of pediatric vaccination. For example, modeling of next-generation influenza vaccines in Kenya demonstrated both direct and indirect effects of pediatric programs, with a 29% to 66% reduction in infections in all ages depending on the product used, achieving cost-effectiveness [ 12 ]. These results align well with estimates from our models. A model from the Netherlands also showed pediatric influenza vaccination to have large indirect effects, with almost half of the symptomatic cases and almost all deaths averted arising from the non-vaccinated age groups [ 10 ]. A United Kingdom (UK) model showed lower reductions on infections but higher reductions in deaths when vaccinating 50% of children aged 6 months to 4 years compared to our model (0.5 infections averted per dose and 1.8 deaths averted per 1,000 doses in UK model vs 1.2 infections averted per dose and 0.4 deaths averted per dose in our model) [ 11 ]. Differences between the two studies are likely due to the difference in age structure between the two settings, the UK having less children (hence less infections and presumably less indirect effects) and a greater population aged β‰₯65 years (hence more deaths per capita) [ 33 ]. Further evidence for the impact of pediatric influenza vaccine programs also comes from observational studies of large-scale vaccine rollouts. The UK extended their influenza vaccination program from the elderly to children aged 6 months to 17 years starting in 2013. Results from comparing primary school pilot to non-pilot areas showed a 59% reduction in general practitioner influenza-like illness consultations in the 2014/15 season [ 19 ], although no significant reductions were observed in the 2015/16 season [ 20 ]. Estimation of vaccine effects from observational ecological studies of time trends is notoriously difficult, and impact estimates may differ considerably between observational studies and transmission models. A US study of school-based vaccination programs revealed notable differences between the two approaches, with observational approaches estimating much higher indirect effects (33%) than mathematical models (2%). These discrepancies between empirical and modeling studies highlight the need to integrate real-world data with mathematical approaches to improve benefits estimates and address potential biases [ 13 ]. The assumption of a vaccine that can reduce 50% of onward transmission by shortening the duration of shedding is not based on a currently available product. For mRNA COVID-19 vaccines, it has been demonstrated that vaccination provides indirect protection to close contacts [ 34 – 36 ], as well as reduces onward transmission from breakthrough infections [ 35 ]. Data are scarce on the indirect effects of vaccinating children for influenza. A community-based randomized control trial investigating the direct and indirect effects of inactivated influenza vaccine (IIV) showed only a 20% reduction in infection probability for household members if a single child aged 5-17 years in the household is vaccinated [ 37 ]. On the other hand, vaccinating children in daycare with IIV resulted in 42% fewer respiratory illnesses in household contacts [ 15 ]. A cluster randomized control trial in which children aged 3-15 years were vaccinated IIV showed a 61% protection against influenza infection in non-vaccinated community members [ 14 ]. Specifically, due to the induction of mucosal immunity, live attenuated influenza vaccine (LAIV) has greater potential to reduce onward transmission. A cluster randomized trial of elementary schools where schools were allocated to receive either IIV or LAIV showed a 64% greater protection in household members against influenza infection when a child was vaccinated with LAIV compared to IIV [ 16 ]. With new influenza vaccine products on the horizon, including >200 universal vaccine candidates on different virus, nucleic acid, virus-like particles (VLP), non-VLP nanoparticles, recombinant protein, and virus vector platforms [ 38 ], it is important to understand the key targets and endpoints that future vaccines should aim for to achieve the greatest impact on transmission and overall disease burden. Our model framework is adaptable to various VE estimates, including reductions in susceptibility, duration and peak of viral shedding, making it a versatile tool to assess and optimize new influenza vaccine products as they become available. A key strength of our model stems from robust parameterization using comprehensive data from multi-year cohorts, providing detailed characterization of influenza transmission dynamics, including pre-season antibody titers, symptomatic fraction, duration of shedding, risk of infection and transmission [ 4 , 24 ], as well as contact rates within these communities [ 25 ]. These data enabled us to significantly reduce the number of parameters that need to be estimated, while reproducing the overall temporal and age dynamics of influenza. In addition, use of a detailed model allowed us to simulate different mechanisms of vaccine protection. A few limitations are worth noting. Firstly, our model was focused on the dynamics of influenza within a single season, and we did not consider waning immunity and viral evolution that may affect multi-year influenza dynamics. We implicitly accounted for inter-annual fluctuations in population immunity by initializing each season based on pre-existing titers to circulating strains in the cohort. We also estimated subtype/lineage-specific transmission rates each season, allowing for differences in force-of-infection between seasons and subtype/lineage. Our model was already highly complex and it would have been difficult to explicitly account for additional mechanisms such as influenza inter-annual variability, especially given the unknown duration of protection following infection. Secondly, while our model was calibrated to data from two South African communities, and may not be generalizable to other countries, we believe our results provide reasonable estimates for communities with similar population structures and contact patterns. Thirdly, recognizing that HAI titers are only one aspect of the immune response to influenza, we did not incorporate the boosting of HAI titers following vaccination or its potential impact on vaccine efficacy. Instead, we applied a pre-specified vaccine efficacy, defined as a percentage reduction in both shedding duration and susceptibility. Lastly, our analysis primarily focused on the evaluation of a vaccine campaign’s indirect effects. We thus ignored the direct effects of vaccination in further reducing mild to severe disease among vaccine breakthrough infections, particularly because the campaign was targeted at children, who generally have a lower risk of severe illness. As a result, our projections of vaccine benefits are likely conservative. In conclusion, our study provides compelling evidence that pediatric influenza immunization programs employing transmission-reducing vaccines can substantially lower the community burden of influenza, particularly in settings with high attack rates among children. Our model offers a valuable tool for determining the minimum efficacy requirement for targeted direct and indirect impact on influenza burden, as well as for evaluating vaccine impact with established efficacy profiles. From a translational perspective, it can inform clinical trial design by assisting in defining target product profiles to achieve a predetermined population impact, and by guiding the development of post-licensure immunization programs to maximize benefits during large-scale rollouts. Data Availability All data produced in the present work are contained in the manuscript. Model code will be made available after publication in peer-review journal. Author contributions Conception and design of the study JK, CV, ST, CC, KS Data collection and laboratory processing JK, ST, NW, JM, AvG, LM, AM, CC Analysis and interpretation JK, CV, MS, ST, NW, JM, AvG, LM, AM, CC, KS Drafting or critical review of the article JK, CV, MS, ST, NW, JM, AvG, LM, AM, CC, KS Data reporting Model code is available at https://github.com/crdm-nicd/phirst_fluvax_transm_model . [ Will be made available upon publication ] Supplementary files Supplementary Excel file 1: Model parameter values. Supplementary tables and figures. Acknowledgement This work was supported by the National Institute for Communicable Diseases of the National Health Laboratory Service and the US Centers for Disease Control and Prevention (cooperative agreement number 5U51IP000155), and the in-house research program of the Division of International Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health. The findings and conclusions in this paper are those of the authors and do not necessarily represent the official position of the NICD, NIH, US CDC. CC has received grant support from Sanofi Pasteur, US CDC, Welcome Trust, Taskforce for Global Health and Bill & Melinda Gates Foundation. AvG has received grant support from US CDC, ASLM/Africa CDC and WHO Afro. NW reports receiving grants from Sanofi Pasteur, US CDC and the Gates Foundation. JM has received grant support from Sanofi Pasteur and PATH. Other authors declare no competing interests. References 1. ↡ Troeger CE , Blacker BF , Khalil IA , Zimsen SRM , Albertson SB , Abate D , et al. Mortality, morbidity, and hospitalisations due to influenza lower respiratory tract infections, 2017: an analysis for the Global Burden of Disease Study 2017 . The Lancet Respiratory Medicine . 2019 ;7(1):69-89. doi: 10.1016/S2213-2600(18)30496-X . OpenUrl CrossRef PubMed 2. ↡ Tempia S , Walaza S , Moyes J , Cohen AL , McMorrow ML , Treurnicht FK , et al. Quantifying How Different Clinical Presentations, Levels of Severity, and Healthcare Attendance Shape the Burden of Influenza-associated Illness: A Modeling Study From South Africa. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America . 2019 ; 69 ( 6 ): 1036 – 48 . doi: 10.1093/cid/ciy1017 . PubMed PMID: 30508065 . OpenUrl CrossRef PubMed 3. ↡ Tempia S , Moyes J , Cohen AL , Walaza S , Edoka I , McMorrow ML , et al. Health and economic burden of influenza-associated illness in South Africa, 2013-2015 . Influenza Other Respir Viruses . 2019 . Epub 2019/06/13. doi: 10.1111/irv.12650 . PubMed PMID: 31187609 . OpenUrl CrossRef PubMed 4. ↡ Cohen C , Kleynhans J , Moyes J , McMorrow ML , Treurnicht FK , Hellferscee O , et al. Asymptomatic transmission and high community burden of seasonal influenza in an urban and a rural community in South Africa, 2017-18 (PHIRST): a population cohort study . The Lancet Global Health . 2021 ;9(6):e863-e74. Epub 2021/05/22. doi: 10.1016/S2214-109X(21)00141-8 . PubMed PMID: 34019838 ; PubMed Central PMCID: PMCPMC8262603 . OpenUrl CrossRef PubMed 5. ↡ Tempia S , Walaza S , Moyes J , Cohen AL , von Mollendorf C , Treurnicht FK , et al. Risk Factors for Influenza-Associated Severe Acute Respiratory Illness Hospitalization in South Africa, 2012–2015 . Open Forum Infect Dis . 2017 ;4(1):ofw262-ofw. doi: 10.1093/ofid/ofw262 . OpenUrl CrossRef 6. ↡ Ferdinands JM , Thompson MG , Blanton L , Spencer S , Grant L , Fry AM . Does influenza vaccination attenuate the severity of breakthrough infections? A narrative review and recommendations for further research . Vaccine . 2021 ; 39 ( 28 ): 3678 – 95 . doi: 10.1016/j.vaccine.2021.05.011 . PubMed PMID: 2540098393 . OpenUrl CrossRef PubMed 7. ↡ World Health Organization . Vaccine against influenza WHO position paper - November 2012 . Wkly Epidemiol Rec . 2012 ; 47 ( 87 ): 461 – 76 . Epub 476. OpenUrl 8. ↡ Tempia S , Walaza S , Moyes J , McMorrow ML , Cohen AL , Edoka I , et al. Influenza disease burden among potential target risk groups for immunization in South Africa, 2013–2015 . Vaccine . 2020 ;38(27):4288-97. doi: 10.1016/j.vaccine.2020.04.045 . PubMed PMID: 2425693001 . OpenUrl CrossRef PubMed 9. ↡ Fraser H , Tombe-Mdewa W , Kohli-Lynch C , Hofman K , Tempia S , McMorrow M , et al. Costs of seasonal influenza vaccination in South Africa . Influenza Other Respi Viruses . 2022 ; 16 ( 5 ): 873 – 80 . doi: 10.1111/irv.12987 . OpenUrl CrossRef 10. ↡ de Boer PT , Nagy L , Dolk FCK , Wilschut JC , Pitman R , Postma MJ . Cost-Effectiveness of Pediatric Influenza Vaccination in The Netherlands . Value Health . 2021 ; 24 ( 1 ): 19 – 31 . doi: 10.1016/j.jval.2020.10.011 . OpenUrl CrossRef PubMed 11. ↡ Baguelin M , Flasche S , Camacho A , Demiris N , Miller E , Edmunds WJ . Assessing Optimal Target Populations for Influenza Vaccination Programmes: An Evidence Synthesis and Modelling Study . PLoS Med . 2013 ; 10 ( 10 ): e1001527 . doi: 10.1371/journal.pmed.1001527 . OpenUrl CrossRef PubMed 12. ↡ Waterlow NR , Radhakrishnan S , Dawa J , van Leeuwen E , Procter SR , Lambach P , et al. Potential health and economic impact of paediatric vaccination using next-generation influenza vaccines in Kenya: a modelling study . BMC Med . 2023 ; 21 ( 1 ): 106 . doi: 10.1186/s12916-023-02830-w . OpenUrl CrossRef PubMed 13. ↡ Arinaminpathy N , Reed C , Biggerstaff M , Nguyen AT , Athni TS , Arnold BF , et al. Estimating community-wide indirect effects of influenza vaccination: triangulation using mathematical models and bias analysis . Am J Epidemiol . 2024 . doi: 10.1093/aje/kwae365 . OpenUrl CrossRef 14. ↡ Loeb M , Russell ML , Moss L , Fonseca K , Fox J , Earn DJD , et al. Effect of Influenza Vaccination of Children on Infection Rates in Hutterite Communities: A Randomized Trial . JAMA . 2010 ; 303 ( 10 ): 943 – 50 . doi: 10.1001/jama.2010.250 . OpenUrl CrossRef PubMed Web of Science 15. ↡ Hurwitz ES , Haber M , Chang A , Shope T , Teo S , Ginsberg M , et al. Effectiveness of Influenza Vaccination of Day Care Children in Reducing Influenza-Related Morbidity Among Household Contacts . JAMA . 2000 ; 284 ( 13 ): 1677 – 82 . doi: 10.1001/jama.284.13.1677 . OpenUrl CrossRef PubMed Web of Science 16. ↡ Kwong JC , Pereira JA , Quach S , Pellizzari R , Dusome E , Russell ML , et al. Randomized evaluation of live attenuated vs. inactivated influenza vaccines in schools (RELATIVES) cluster randomized trial: Pilot results from a household surveillance study to assess direct and indirect protection from influenza vaccination . Vaccine . 2015 ; 33 ( 38 ): 4910 – 5 . doi: 10.1016/j.vaccine.2015.07.044 . OpenUrl CrossRef PubMed 17. Rudenko LG , Slepushkin AN , Monto AS , Kendal AP , Grigorieva EP , Burtseva EP , et al. Efficacy of Live Attenuated and Inactivated Influenza Vaccines in Schoolchildren and Their Unvaccinated Contacts in Novgorod, Russia . The Journal of Infectious Diseases . 1993 ; 168 ( 4 ): 881 – 7 . doi: 10.1093/infdis/168.4.881 . OpenUrl CrossRef PubMed Web of Science 18. Matsuoka T , Sato T , Akita T , Yanagida J , Ohge H , Kuwabara M , et al. High Vaccination Coverage among Children during Influenza A(H1N1)pdm09 as a Potential Factor of Herd Immunity . Int J Environ Res Public Health . 2016 ; 13 ( 10 ). Epub 2016/10/21. doi: 10.3390/ijerph13101017 . PubMed PMID: 27763532 ; PubMed Central PMCID: PMCPMC5086756 . OpenUrl CrossRef PubMed 19. ↡ Pebody RG , Green HK , Andrews N , Boddington NL , Zhao H , Yonova I , et al. Uptake and impact of vaccinating school age children against influenza during a season with circulation of drifted influenza A and B strains , England , 2014 /15. Euro Surveill. 2015;20(39). doi: 10.2807/1560-7917.ES.2015.20.39.30029 . PubMed PMID: 26537222 . OpenUrl CrossRef PubMed 20. ↡ Pebody RG , Sinnathamby MA , Warburton F , Andrews N , Boddington NL , Zhao H , et al. Uptake and impact of vaccinating primary school-age children against influenza: experiences of a live attenuated influenza vaccine programme , England , 2015 /16. Euro Surveill. 2018;23(25). Epub 2018/06/28. doi: 10.2807/1560-7917.Es.2018.23.25.1700496 . PubMed PMID: 29945698 ; PubMed Central PMCID: PMCPMC6152241 . OpenUrl CrossRef PubMed 21. ↡ Nayak J , Hoy G , Gordon A . Influenza in Children . Cold Spring Harb Perspect Med . 2021 ; 11 ( 1 ). doi: 10.1101/cshperspect.a038430 . OpenUrl Abstract / FREE Full Text 22. ↡ Kamigaki T , Mimura S , Takahashi Y , Oshitani H . Analysis of influenza transmission in the households of primary and junior high school students during the 2012-13 influenza season in Odate, Japan . BMC Infect Dis . 2015 ; 15 : 282 -. doi: 10.1186/s12879-015-1007-8 . PubMed PMID: 26201316 . OpenUrl CrossRef PubMed 23. ↡ Cohen C , McMorrow ML , Martinson NA , Kahn K , Treurnicht FK , Moyes J , et al. Cohort profile: A Prospective Household cohort study of Influenza, Respiratory syncytial virus and other respiratory pathogens community burden and Transmission dynamics in South Africa, 2016–2018 . Influenza Other Respi Viruses . 2021 ;15(6):789-803. Epub 2021/07/24. doi: 10.1111/irv.12881 . PubMed PMID: 34296810 ; PubMed Central PMCID: PMCPMC8542945 . OpenUrl CrossRef PubMed 24. ↡ Sauter MK , Kleynhans J , Moyes J , McMorrow ML , Treurnicht FK , Hellferscee O , et al. The roles of pre-season immunity, age, viral shedding, and community exposures in shaping influenza household transmission dynamics . medRxiv . 2025 :2025.03.25.25324622. doi: 10.1101/2025.03.25.25324622 . OpenUrl Abstract / FREE Full Text 25. ↡ Kleynhans J , Tempia S , McMorrow ML , von Gottberg A , Martinson NA , Kahn K , et al. A cross- sectional study measuring contact patterns using diaries in an urban and a rural community in South Africa, 2018 . BMC Public Health . 2021 ;21:1055. Epub 2021/06/04. doi: 10.1186/s12889-021-11136-6 . PubMed PMID: 34078327 ; PubMed Central PMCID: PMCPMC8172361 . OpenUrl CrossRef PubMed 26. ↡ Canini L , Carrat F . Population Modeling of Influenza A/H1N1 Virus Kinetics and Symptom Dynamics . J Virol . 2011 ; 85 ( 6 ): 2764 – 70 . doi : doi: 10.1128/jvi.01318-10 . OpenUrl Abstract / FREE Full Text 27. ↡ Guo J , Chen X , Guo Y , Liu M , Li P , Tao Y , et al. Real-world effectiveness of seasonal influenza vaccination and age as effect modifier: A systematic review, meta-analysis and meta-regression of test-negative design studies . Vaccine . 2024 ; 42 ( 8 ): 1883 – 91 . doi: 10.1016/j.vaccine.2024.02.059 . PubMed PMID: 2956720062 . OpenUrl CrossRef PubMed 28. ↡ Ng S , Lopez R , Kuan G , Gresh L , Balmaseda A , Harris E , et al. The Timeline of Influenza Virus Shedding in Children and Adults in a Household Transmission Study of Influenza in Managua, Nicaragua . The Pediatric Infectious Disease Journal . 2016 ; 35 ( 5 ): 583 – 6 . doi: 10.1097/inf.0000000000001083 . PubMed PMID: 00006454 -201605000-00029. OpenUrl CrossRef PubMed 29. Roosenhoff R , Reed V , Kenwright A , Schutten M , Boucher CA , Monto A , et al. Viral Kinetics and Resistance Development in Children Treated with Neuraminidase Inhibitors: The Influenza Resistance Information Study (IRIS) . Clin Infect Dis . 2020 ; 71 ( 5 ): 1186 – 94 . doi: 10.1093/cid/ciz939 . OpenUrl CrossRef 30. ↡ Tsang TK , Cowling BJ , Fang VJ , Chan K-H , Ip DKM , Leung GM , et al. Influenza A Virus Shedding and Infectivity in Households . The Journal of Infectious Diseases . 2015 ; 212 ( 9 ): 1420 – 8 . Epub 2015/04/18. doi: 10.1093/infdis/jiv225 . PubMed PMID: 25883385 ; PubMed Central PMCID: PMCPMC4601913 . OpenUrl CrossRef PubMed 31. 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 – 6 . doi: 10.1016/j.epidem.2015.04.003 . OpenUrl CrossRef PubMed 32. ↡ Gordon A , Tsang TK , Cowling BJ , Kuan G , Ojeda S , Sanchez N , et al. Influenza Transmission Dynamics in Urban Households, Managua, Nicaragua, 2012-2014 . Emerg Infect Dis . 2018 ;24(10):1882-8. Epub 2018/09/19. doi: 10.3201/eid2410.161258 . PubMed PMID: 30226161 ; PubMed Central PMCID: PMCPMC6154158 . OpenUrl CrossRef PubMed 33. ↡ 33. Office of National Statistics. Population estimates for the UK, England, Wales, Scotland, and Northern Ireland: mid-2022 . 2025 21 January 2025]]. Available from: https://www.ons.gov.uk/peoplepopulationandcommunity/populationandmigration/populationestimates/bulletins/annualmidyearpopulationestimates/mid2022 . 34. ↡ Tan ST , RodrΓ­guez-Barraquer I , Kwan AT , Blumberg S , Park HJ , Hutchinson J , et al. Strength and durability of indirect protection against SARS-CoV-2 infection through vaccine and infection- acquired immunity . Nature communications . 2025 ; 16 ( 1 ): 1090 . doi: 10.1038/s41467-024-55029-9 . OpenUrl CrossRef PubMed 35. ↡ Hayek S , Shaham G , Ben-Shlomo Y , Kepten E , Dagan N , Nevo D , et al. Indirect protection of children from SARS-CoV-2 infection through parental vaccination . Science . 2022 ; 375 (6585): 1155 -9. doi : doi: 10.1126/science.abm3087 . OpenUrl CrossRef 36. ↡ Salo J , HΓ€gg M , Kortelainen M , Leino T , Saxell T , Siikanen M , et al. The indirect effect of mRNA-based COVID-19 vaccination on healthcare workers’ unvaccinated household members . Nature communications . 2022 ; 13 ( 1 ): 1162 . doi: 10.1038/s41467-022-28825-4 . OpenUrl CrossRef PubMed 37. ↡ Tsang TK , Fang VJ , Ip DKM , Perera RAPM , So HC , Leung GM , et al. Indirect protection from vaccinating children against influenza in households . Nature communications . 2019 ; 10 ( 1 ): 106 . doi: 10.1038/s41467-018-08036-6 . OpenUrl CrossRef PubMed 38. ↡ University of Minnesota. Universal Influenza Vaccine Technology Landscape . 2024 11 October 2024]]. Available from: https://ivr.cidrap.umn.edu/universal-influenza-vaccine-technology-landscape . View the discussion thread. Back to top Previous Next Posted April 04, 2025. Download PDF Data/Code Email Thank you for your interest in spreading the word about medRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Projections of the impact of transmission-reducing pediatric influenza vaccine in South Africa using high resolution immunologic and infection data Message Subject (Your Name) has forwarded a page to you from medRxiv Message Body (Your Name) thought you would like to see this page from the medRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Projections of the impact of transmission-reducing pediatric influenza vaccine in South Africa using high resolution immunologic and infection data Jackie Kleynhans , CΓ©cile Viboud , Molly Sauter , Stefano Tempia , Nicole Wolter , Jocelyn Moyes , Anne von Gottberg , Lorens Maake , Alexandra Moerdyk , Cheryl Cohen , Kaiyuan Sun medRxiv 2025.04.04.25325228; doi: https://doi.org/10.1101/2025.04.04.25325228 Share This Article: Copy Citation Tools Projections of the impact of transmission-reducing pediatric influenza vaccine in South Africa using high resolution immunologic and infection data Jackie Kleynhans , CΓ©cile Viboud , Molly Sauter , Stefano Tempia , Nicole Wolter , Jocelyn Moyes , Anne von Gottberg , Lorens Maake , Alexandra Moerdyk , Cheryl Cohen , Kaiyuan Sun medRxiv 2025.04.04.25325228; doi: https://doi.org/10.1101/2025.04.04.25325228 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Infectious Diseases (except HIV/AIDS) Subject Areas All Articles Addiction Medicine (569) Allergy and Immunology (863) Anesthesia (300) Cardiovascular Medicine (4442) Dentistry and Oral Medicine (444) Dermatology (383) Emergency Medicine (609) Endocrinology (including Diabetes Mellitus and Metabolic Disease) (1510) Epidemiology (15230) Forensic Medicine (30) Gastroenterology (1126) Genetic and Genomic Medicine (6608) Geriatric Medicine (668) Health Economics (998) Health Informatics (4542) Health Policy (1370) Health Systems and Quality Improvement (1613) Hematology (543) HIV/AIDS (1266) Infectious Diseases (except HIV/AIDS) (15923) Intensive Care and Critical Care Medicine (1103) Medical Education (623) Medical Ethics (147) Nephrology (668) Neurology (6607) Nursing (346) Nutrition (999) Obstetrics and Gynecology (1146) Occupational and Environmental Health (957) Oncology (3336) Ophthalmology (974) Orthopedics (369) Otolaryngology (420) Pain Medicine (436) Palliative Medicine (130) Pathology (664) Pediatrics (1693) Pharmacology and Therapeutics (692) Primary Care Research (712) Psychiatry and Clinical Psychology (5448) Public and Global Health (9237) Radiology and Imaging (2201) Rehabilitation Medicine and Physical Therapy (1370) Respiratory Medicine (1196) Rheumatology (596) Sexual and Reproductive Health (714) Sports Medicine (530) Surgery (712) Toxicology (99) Transplantation (289) Urology (265) (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'a01843234ddb58d3',t:'MTc3OTc1MTEzNw=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();

Text is read by the "Ask this paper" AI Q&A widget below. Extraction quality varies by source β€” PMC NXML preserves structure cleanly, OA-HTML may include some navigation residue, and OA-PDF can have broken hyphenation. The publisher copy (via DOI) is the canonical version.

My notes (saved in your browser only)

βš™ Ask this paper AI returns verbatim quotes from the full text Β· source: preprint-html β“˜

Answers must be backed by verbatim quotes from this paper's full text. Hallucinated quotes are dropped automatically; if no verbatim passage answers the question, we say so. How this works

Citation neighborhood (no data yet)

We don't have any in-corpus citations linked to this paper yet. This is a recent paper (2025) β€” citers typically take a year or two to land, and the OpenAlex reference graph may still be filling in.

Source provenance

europepmc
last seen: 2026-05-20T01:45:00.602351+00:00