Search

Abstract The SARS-CoV-2 BA.2.86 lineage, and its sublineage JN.1 in particular, achieved widespread transmission in the US during winter 2023–24. However, this surge in infections was not accompanied by COVID-19 hospitalizations and mortality commensurate with prior waves. To understand shifts in COVID-19 epidemiology associated with JN.1 emergence, we compared characteristics and clinical outcomes of time-matched cases infected with BA.2.86 lineages (predominantly representing JN.1) versus co-circulating XBB-derived lineages in December, 2023 and January, 2024. Cases infected with BA.2.86 lineages received greater numbers of COVID-19 vaccine doses, including XBB.1.5-targeted boosters, in comparison to cases infected with XBB-derived lineages. Additionally, cases infected with BA.2.86 lineages experienced greater numbers of documented prior SARS-CoV-2 infections. Cases infected with BA.2.86 lineages also experienced lower risk of progression to severe clinical outcomes requiring emergency department consultations or hospital admission. Sensitivity analyses suggested under-ascertainment of prior infections could not explain this apparent attenuation of severity. Our findings implicate escape from immunity acquired from prior vaccination or infection in the emergence of the JN.1 lineage and suggest infections with this lineage are less likely to experience clinically-severe disease. Monitoring of immune escape and clinical severity in emerging SARS-CoV-2 variants remains a priority to inform responses.

ABSTRACT After introducing pneumococcal conjugate vaccines (PCVs), serotype replacement occurred in Streptococcus pneumoniae. Predicting which pneumococcal strains will become common in carriage after vaccination can enhance vaccine design, public health interventions, and understanding of pneumococcal evolution. Invasive pneumococcal isolates were collected during 1998–2018 by the Active Bacterial Core surveillance (ABCs). Carriage data from Massachusetts (MA) and Southwest United States were used to calculate weights. Using pre-vaccine data, serotype-specific inverse-invasiveness weights were defined as the ratio of the proportion of the serotype in carriage to the proportion in invasive data. Genomic data were processed under bioinformatic pipelines to define genetically similar sequence clusters (i.e., strains), and accessory genes (COGs) present in 5–95% of isolates. Weights were applied to adjust observed strain proportions and COG frequencies. The negative frequency-dependent selection (NFDS) model predicted strain proportions by calculating the post-vaccine strain composition in the weighted invasive disease population that would best match pre-vaccine COG frequencies. Inverse-invasiveness weighting increased the correlation of COG frequencies between invasive and carriage data in linear or logit scale for pre-vaccine, post-PCV7, and post-PCV13; and between different epochs in the invasive data. Weighting the invasive data significantly improved the NFDS model’s accuracy in predicting strain proportions in the carriage population in the post-PCV13 epoch, with the adjusted R2 increasing from 0.254 before weighting to 0.545 after weighting. The weighting system adjusted invasive disease data to better represent the pneumococcal carriage population, allowing the NFDS mechanism to predict strain proportions in carriage in the post-PCV13 epoch. Our methods enrich the value of genomic sequences from invasive disease surveillance.IMPORTANCEStreptococcus pneumoniae, a common colonizer in the human nasopharynx, can cause invasive diseases including pneumonia, bacteremia, and meningitis mostly in children under 5 years or older adults. The PCV7 was introduced in 2000 in the United States within the pediatric population to prevent disease and reduce deaths, followed by PCV13 in 2010, PCV15 in 2022, and PCV20 in 2023. After the removal of vaccine serotypes, the prevalence of carriage remained stable as the vacated pediatric ecological niche was filled with certain non-vaccine serotypes. Predicting which pneumococcal clones, and which serotypes, will be most successful in colonization after vaccination can enhance vaccine design and public health interventions, while also improving our understanding of pneumococcal evolution. While carriage data, which are collected from the pneumococcal population that is competing to colonize and transmit, are most directly relevant to evolutionary studies, invasive disease data are often more plentiful. Previously, evolutionary models based on negative frequency-dependent selection (NFDS) on the accessory genome were shown to predict which non-vaccine strains and serotypes were most successful in colonization following the introduction of PCV7. Here, we show that an inverse-invasiveness weighting system applied to invasive disease surveillance data allows the NFDS model to predict strain proportions in the projected carriage population in the post-PCV13/pre-PCV15 and pre-PCV20 epoch. The significance of our research lies in using a sample of invasive disease surveillance data to extend the use of NFDS as an evolutionary mechanism to predict post-PCV13 population dynamics. This has shown that we can correct for biased sampling that arises from differences in virulence and can enrich the value of genomic data from disease surveillance and advance our understanding of how NFDS impacts carriage population dynamics after both PCV7 and PCV13 vaccination.

Background: Nursing home (NH) residents are at high risk of COVID-19 from exposure to infected staff and other residents. Understanding SARS-CoV-2 viral RNA kinetics in residents and staff can guide testing, isolation, and return to work recommendations. We sought to determine the duration of antigen test and polymerase chain reaction (PCR) positivity in a cohort of NH residents and staff. Methods: We prospectively collected data on SARS-CoV-2 viral kinetics from April 2023 through November 2023. Staff and residents could enroll prospectively or upon a positive test (identified through routine clinical testing, screening, or outbreak response testing). Participating facilities performed routine clinical testing; asymptomatic testing of contacts was performed within 48 hours if an outbreak or known exposure occurred and upon (re-) admission. Enrolled participants who tested positive for SARS-CoV-2 were re-tested daily for 14 days with both nasal antigen and nasal PCR tests. All PCR tests were run by a central lab with the same assay. We conducted a Kaplan-Meier survival analysis on time to first negative test restricted to participants who initially tested positive (day zero) and had at least one test ≥10 days after initially testing positive with the same test type; a participant could contribute to both antigen and PCR survival curves. We compared survival curves for staff and residents using the log-rank test. Results: Twenty-four nursing homes in eight states participated; 587 participants (275 residents, 312 staff) enrolled in the evaluation, participants were only tested through routine clinical or outbreak response testing. Seventy-two participants tested positive for antigen; of these, 63 tested PCR-positive. Residents were antigen- and PCR-positive longer than staff (Figure 1), but this finding is only statistically significant (p=0.006) for duration of PCR positivity. Five days after the first positive test, 56% of 50 residents and 59% of 22 staff remained antigen-positive; 91% of 44 residents and 79% of 19 staff were PCR-positive. Ten days after the first positive test, 22% of 50 residents and 5% of 22 staff remained antigen-positive; 61% of 44 residents and 21% of 19 staff remained PCR-positive. Conclusions: Most NH residents and staff with SARS-CoV-2 remained antigen- or PCR-positive 5 days after the initial positive test; however, differences between staff and resident test positivity were noted at 10 days. These data can inform recommendations for testing, duration of NH resident isolation, and return to work guidance for staff. Additional viral culture data may strengthen these conclusions.

ABSTRACT When respiratory viruses co-circulate in a population, individuals may be infected with multiple pathogens and experience possible virus–virus interactions, where concurrent or recent prior infection with one virus affects the infection process of another virus. While experimental studies have provided convincing evidence for within-host mechanisms of virus–virus interactions, evaluating evidence for viral interference or potentiation using population-level data has proven more difficult. Recent studies have quantified the prevalence of co-detections using populations drawn from clinical settings. Here, we focus on selection bias issues associated with this study design. We provide a quantitative account of the conditions under which selection bias arises in these studies, review previous attempts to address this bias, and propose unbiased study designs with sample size estimates needed to ascertain viral interference. We show that selection bias is expected in cross-sectional co-detection prevalence studies conducted in clinical settings, except under a strict set of assumptions regarding the relative probabilities of being included in a study limited to individuals with clinical disease under different viral states. Population-wide studies that collect samples from participants irrespective of their clinical status would meanwhile require large sample sizes to be sufficiently powered to detect viral interference, suggesting that a study’s timing, inclusion criteria, and the expected magnitude of interference are instrumental in determining feasibility.

The COVID-19 pandemic has highlighted the need to upgrade systems for infectious disease surveillance and forecasting and modeling of the spread of infection, both of which inform evidence-based public health guidance and policies. Here, we discuss requirements for an effective surveillance system to support decision making during a pandemic, drawing on the lessons of COVID-19 in the U.S., while looking to jurisdictions in the U.S. and beyond to learn lessons about the value of specific data types. In this report, we define the range of decisions for which surveillance data are required, the data elements needed to inform these decisions and to calibrate inputs and outputs of transmission-dynamic models, and the types of data needed to inform decisions by state, territorial, local, and tribal health authorities. We define actions needed to ensure that such data will be available and consider the contribution of such efforts to improving health equity.

Public health surveillance for pathogens presents an optimization problem: we require enough sampling to identify intervention-triggering shifts in pathogen epidemiology, such as new introductions or sudden increases in prevalence, but not so much that costs due to surveillance itself outweigh those from pathogen-associated illness. To determine this optimal sampling frequency, we developed a general mathematical model for the introduction of a new pathogen that, once introduced, increases in prevalence exponentially. Given the relative cost of infection vs. sampling, we derived equations for the expected combined cost per unit time of disease burden and surveillance for a specified sampling frequency, and thus the sampling frequency for which the expected total cost per unit time is lowest.

Abstract Host immune responses are a key source of selective pressure driving pathogen evolution. Emergence of many SARS-CoV-2 lineages has been associated with enhancements in their ability to evade population immunity resulting from both vaccination and infection. Here we show diverging trends of escape from vaccine-derived and infection-derived immunity for the emerging XBB/XBB.1.5 Omicron lineage. Among 31,739 patients tested in ambulatory settings in Southern California from December, 2022 to February, 2023, adjusted odds of prior receipt of 2, 3, 4, and ≥5 COVID-19 vaccine doses were 10% (95% confidence interval: 1–18%), 11% (3–19%), 13% (3–21%), and 25% (15–34%) lower, respectively, among cases infected with XBB/XBB.1.5 than among cases infected with other co-circulating lineages. Similarly, prior vaccination was associated with greater point estimates of protection against progression to hospitalization among cases with XBB/XBB.1.5 than among non-XBB/XBB.1.5 cases (70% [30–87%] and 48% [7–71%], respectively, for recipients of ≥4 doses). In contrast, cases infected with XBB/XBB.1.5 had 17% (11–24%) and 40% (19–65%) higher adjusted odds of having experienced 1 and ≥2 prior documented infections, respectively, including with pre-Omicron variants. As immunity acquired from SARS-CoV-2 infection becomes increasingly widespread, fitness costs associated with enhanced vaccine sensitivity in XBB/XBB.1.5 may be offset by increased ability to evade infection-derived host responses.

The COVID-19 epidemic in the United States has been characterized by two stark disparities. COVID-19 burden has been unequally distributed among racial and ethnic groups and at the same time the mortality rates have been sharply higher among older age groups. These disparities have led some to suggest that inequalities could be reduced by vaccinating front-line workers before vaccinating older individuals, as older individuals in the US are disproportionately Non-Hispanic White. We compare the performance of two distribution policies, one allocating vaccines to front-line workers and another to older individuals aged 65-74-year-old. We estimate both the number of lives saved and the number of years of life saved under each of the policies, overall and in every race/ethnicity groups, in the United States and every state. We show that prioritizing COVID-19 vaccines for 65-74-year-olds saves both more lives and more years of life than allocating vaccines front-line workers in each racial/ethnic group, in the United States as a whole and in nearly every state. When evaluating fairness of vaccine allocation policies, the overall benefit to impact of each population subgroup should be considered, not only the proportion of doses that is distributed to each subgroup. Further work can identify prioritization schemes that perform better on multiple equity metrics.

Social gatherings can be an important locus of transmission for many pathogens including SARS-CoV-2. During an outbreak, restricting the size of these gatherings is one of several non-pharmaceutical interventions available to policy-makers to reduce transmission. Often these restrictions take the form of prohibitions on gatherings above a certain size. While it is generally agreed that such restrictions reduce contacts, the specific size threshold separating “allowed” from “prohibited” gatherings often does not have a clear scientific basis, which leads to dramatic differences in guidance across location and time. Building on the observation that gathering size distributions are often heavy-tailed, we develop a theoretical model of transmission during gatherings and their contribution to general disease dynamics. We find that a key, but often overlooked, determinant of the optimal threshold is the distribution of gathering sizes. Using data on pre-pandemic contact patterns from several sources as well as empirical estimates of transmission parameters for SARS-CoV-2, we apply our model to better understand the relationship between restriction threshold and reduction in cases. We find that, under reasonable transmission parameter ranges, restrictions may have to be set quite low to have any demonstrable effect on cases due to relative frequency of smaller gatherings. We compare our conceptual model with observed changes in reported contacts during lockdown in March of 2020.

Sparse testing early in the SARS-CoV-2 pandemic hinders estimation of the dates and origins of initial case importations. Here, the authors show that the main source of cases imported from China shifted from Wuhan to other Chinese cities by mid-February, especially for African locations.

Abstract Background Staphylococcus aureus is a leading cause of healthcare- and community-associated infections and can be difficult to treat due to antimicrobial resistance. About 30% of individuals carry S. aureus asymptomatically in their nares, a risk factor for later infection, and interactions with other species in the nasal microbiome likely modulate its carriage. It is thus important to identify ecological or functional genetic elements within the maternal or infant nasal microbiomes that influence S. aureus acquisition and retention in early life. Results We recruited 36 mother-infant pairs and profiled a subset of monthly longitudinal nasal samples from the first year after birth using shotgun metagenomic sequencing. The infant nasal microbiome is highly variable, particularly within the first 2 months. It is weakly influenced by maternal nasal microbiome composition, but primarily shaped by developmental and external factors, such as daycare. Infants display distinctive patterns of S. aureus carriage, positively associated with Acinetobacter species, Streptococcus parasanguinis, Streptococcus salivarius, and Veillonella species and inversely associated with maternal Dolosigranulum pigrum. Furthermore, we identify a gene family, likely acting as a taxonomic marker for an unclassified species, that is significantly anti-correlated with S. aureus in infants and mothers. In gene content-based strain profiling, infant S. aureus strains are more similar to maternal strains. Conclusions This improved understanding of S. aureus colonization is an important first step toward the development of novel, ecological therapies for controlling S. aureus carriage.

Understanding how antibiotic use drives resistance is crucial for guiding effective strategies to limit the spread of resistance, but the use-resistance relationship across pathogens and antibiotics remains unclear. We applied sinusoidal models to evaluate the seasonal use-resistance relationship across 3 species (Staphylococcus aureus, Escherichia coli, and Klebsiella pneumoniae) and 5 antibiotic classes (penicillins, macrolides, quinolones, tetracyclines, and nitrofurans) in Boston, Massachusetts. Outpatient use of all 5 classes and resistance in inpatient and outpatient isolates in 9 of 15 species-antibiotic combinations showed statistically significant amplitudes of seasonality (false discovery rate (FDR) < 0.05). While seasonal peaks in use varied by class, resistance in all 9 species-antibiotic combinations peaked in the winter and spring. The correlations between seasonal use and resistance thus varied widely, with resistance to all antibiotic classes being most positively correlated with use of the winter peaking classes (penicillins and macrolides). These findings challenge the simple model of antibiotic use independently selecting for resistance and suggest that stewardship strategies will not be equally effective across all species and antibiotics. Rather, seasonal selection for resistance across multiple antibiotic classes may be dominated by use of the most highly prescribed antibiotic classes, penicillins and macrolides.

When responding to infectious disease outbreaks, rapid and accurate estimation of the epidemic trajectory is critical. However, two common data collection problems affect the reliability of the epidemiological data in real time: missing information on the time of first symptoms, and retrospective revision of historical information, including right censoring. Here, we propose an approach to construct epidemic curves in near real time that addresses these two challenges by 1) imputation of dates of symptom onset for reported cases using a dynamically-estimated "backward" reporting delay conditional distribution, and 2) adjustment for right censoring using the NobBS software package to nowcast cases by date of symptom onset. This process allows us to obtain an approximation of the time-varying reproduction number (Rt) in real time. We apply this approach to characterize the early SARS-CoV-2 outbreak in two Spanish regions between March and April 2020. We evaluate how these real-time estimates compare with more complete epidemiological data that became available later. We explore the impact of the different assumptions on the estimates, and compare our estimates with those obtained from commonly used surveillance approaches. Our framework can help improve accuracy, quantify uncertainty, and evaluate frequently unstated assumptions when recovering the epidemic curves from limited data obtained from public health systems in other locations.

Background: Mathematical models predict that community-wide access to HIV testing-and-treatment can rapidly and substantially reduce new HIV infections. Yet several large universal test-and-treat HIV prevention trials in high-prevalence epidemics demonstrated variable reduction in population-level incidence. Methods: To elucidate patterns of HIV spread in universal test-and-treat trials, we quantified the contribution of geographic-location, gender, age, and randomized-HIV-intervention to HIV transmissions in the 30-community Ya Tsie trial in Botswana. We sequenced HIV viral whole genomes from 5114 trial participants among the 30 trial communities. Results: Deep-sequence phylogenetic analysis revealed that most inferred HIV transmissions within the trial occurred within the same or between neighboring communities, and between similarly aged partners. Transmissions into intervention communities from control communities were more common than the reverse post-baseline (30% [12.2 – 56.7] vs. 3% [0.1 – 27.3]) than at baseline (7% [1.5 – 25.3] vs. 5% [0.9 – 22.9]) compatible with a benefit from treatment-as-prevention. Conclusions: Our findings suggest that population mobility patterns are fundamental to HIV transmission dynamics and to the impact of HIV control strategies. Funding: This study was supported by the National Institute of General Medical Sciences (U54GM088558), the Fogarty International Center (FIC) of the U.S. National Institutes of Health (D43 TW009610), and the President’s Emergency Plan for AIDS Relief through the Centers for Disease Control and Prevention (CDC) (Cooperative agreements U01 GH000447 and U2G GH001911).

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection exported from mainland China could lead to self-sustained outbreaks in other countries. By February 2020, several countries were reporting imported SARS-CoV-2 cases. To contain the virus, early detection of imported SARS-CoV-2 cases is critical. We used air travel volume estimates from Wuhan, China, to international destinations and a generalized linear regression model to identify locations that could have undetected imported cases. Our model can be adjusted to account for exportation of cases from other locations as the virus spreads and more information on importations and transmission becomes available. Early detection and appropriate control measures can reduce the risk for transmission in all locations.

Abstract Background Herd protection through interruption of transmission has contributed greatly to the impact of pneumococcal conjugate vaccines (PCVs) and may enable the use of cost-saving reduced dose schedules. To aid PCV age targeting to achieve herd protection, we estimated which population age groups contribute most to vaccine serotype (VT) pneumococcal transmission. Methods We used transmission dynamic models to mirror pre-PCV epidemiology in England and Wales, Finland, Kilifi in Kenya and Nha Trang in Vietnam where data on carriage prevalence in infants, pre-school and school-aged children and adults as well as social contact patterns was available. We used Markov Chain Monte Carlo methods to fit the models and then extracted the per capita and population-based contribution of different age groups to VT transmission. Results We estimated that in all settings,

Abstract Background Antibiotic use contributes to the rates of sepsis and the associated mortality, particularly through lack of clearance of resistant infections following antibiotic treatment. At the same time, there is limited information on the effects of prescribing of some antibiotics vs. others on subsequent sepsis and sepsis-related mortality. Methods We used a multivariable mixed-effects model to relate state-specific rates of outpatient prescribing overall for oral fluoroquinolones, penicillins, macrolides, and cephalosporins between 2014 and 2015 to state-specific rates of mortality with sepsis (ICD-10 codes A40–41 present as either underlying or contributing causes of death on a death certificate) in different age groups of US adults between 2014 and 2015, adjusting for additional covariates and random effects associated with the ten US Health and Human Services (HHS) regions. Results Increase in the rate of prescribing of oral penicillins by 1 annual dose per 1000 state residents was associated with increases in annual rates of mortality with sepsis of 0.95 (95% CI (0.02,1.88)) per 100,000 persons aged 75-84y, and of 2.97 (0.72,5.22) per 100,000 persons aged 85 + y. Additionally, the percent of individuals aged 50-64y lacking health insurance, as well as the percent of individuals aged 65-84y who are African-American were associated with rates of mortality with sepsis in the corresponding age groups. Conclusions Our results suggest that prescribing of penicillins is associated with rates of mortality with sepsis in older US adults. Those results, as well as the related epidemiological data suggest that replacement of certain antibiotics, particularly penicillins in the treatment of different syndromes should be considered with the aim of reducing the rates of severe outcomes, including mortality related to bacterial infections.

Abstract Background Antibiotics remain the cornerstone of modern medicine. Yet there exists an inherent dilemma in their use: we are able to prevent harm by administering antibiotic treatment as necessary to both humans and animals, but we must be mindful of limiting the spread of resistance and safeguarding the efficacy of antibiotics for current and future generations. Policies that strike the right balance must be informed by a transparent rationale that relies on a robust evidence base. Main text One way to generate the evidence base needed to inform policies for managing antibiotic resistance is by using mathematical models. These models can distil the key drivers of the dynamics of resistance transmission from complex infection and evolutionary processes, as well as predict likely responses to policy change in silico. Here, we ask whether we know enough about antibiotic resistance for mathematical modelling to robustly and effectively inform policy. We consider in turn the challenges associated with capturing antibiotic resistance evolution using mathematical models, and with translating mathematical modelling evidence into policy. Conclusions We suggest that in spite of promising advances, we lack a complete understanding of key principles. From this we advocate for priority areas of future empirical and theoretical research.

Background: SARS-CoV-2 variant Beta (B.1.351) was designated as a Variant of Concern (VoC) after becoming the dominant strain in South Africa and spreading internationally. BNT162b2 showed lower levels of neutralizing antibodies against Beta than against other strains raising concerns about effectiveness of vaccines against infections caused by Beta. We estimated BNT162b2 vaccine effectiveness (VE) against Beta infections in Israel, a country with high vaccine uptake. Methods: The Ministry of Health (MoH) identified Beta cases through mandatory reporting of SARS-CoV-2 cases and whole genome sequencing (WGS) of specimens from vaccination-breakthrough infections, reinfections, arriving international travelers, and a selection of other infected persons. A cohort analysis was conducted of exposure events of contacts of primary Beta cases. WGS was conducted on available PCR-positive specimens collected from contacts. VE estimates with 95% confidence intervals (CIs) against confirmed and probable Beta infections were determined by comparing infection risk between unvaccinated and fully-vaccinated (≥7 days after the second dose) contacts, and between unvaccinated and partially-vaccinated (

Background: BNT162b2 was shown to be 92% effective in preventing COVID-19. Prioritizing vaccine rollout, and achievement of herd immunity depend on SARS-CoV-2 transmission reduction. The vaccine's effect on infectivity is thus a critical priority. Methods: Among all 9650 HCW of a large tertiary medical center in Israel, we calculated the prevalence of positive SARS-CoV-2 qRT-PCR cases with asymptomatic presentation, tested following known or presumed exposure and the infectious subset (N-gene-Ct-value

Effectively designing and evaluating public health responses to the ongoing COVID-19 pandemic requires accurate estimation of the prevalence of COVID-19 across the United States (US). Equipment shortages and varying testing capabilities have however hindered the usefulness of the official reported positive COVID-19 case counts. We introduce four complementary approaches to estimate the cumulative incidence of symptomatic COVID-19 in each state in the US as well as Puerto Rico and the District of Columbia, using a combination of excess influenza-like illness reports, COVID-19 test statistics, COVID-19 mortality reports, and a spatially structured epidemic model. Instead of relying on the estimate from a single data source or method that may be biased, we provide multiple estimates, each relying on different assumptions and data sources. Across our four approaches emerges the consistent conclusion that on April 4, 2020, the estimated case count was 5 to 50 times higher than the official positive test counts across the different states. Nationally, our estimates of COVID-19 symptomatic cases as of April 4 have a likely range of 2.3 to 4.8 million, with possibly as many as 7.6 million cases, up to 25 times greater than the cumulative confirmed cases of about 311,000. Extending our methods to May 16, 2020, we estimate that cumulative symptomatic incidence ranges from 4.9 to 10.1 million, as opposed to 1.5 million positive test counts. The proposed combination of approaches may prove useful in assessing the burden of COVID-19 during resurgences in the US and other countries with comparable surveillance systems.

Background: The impact of variable infection risk by race and ethnicity on the dynamics of SARS-CoV-2 spread is largely unknown. Methods: Here, we fit structured compartmental models to seroprevalence data from New York State and analyze how herd immunity thresholds (HITs), final sizes, and epidemic risk change across groups. Results: A simple model where interactions occur proportionally to contact rates reduced the HIT, but more realistic models of preferential mixing within groups increased the threshold toward the value observed in homogeneous populations. Across all models, the burden of infection fell disproportionately on minority populations: in a model fit to Long Island serosurvey and census data, 81% of Hispanics or Latinos were infected when the HIT was reached compared to 34% of non-Hispanic whites. Conclusions: Our findings, which are meant to be illustrative and not best estimates, demonstrate how racial and ethnic disparities can impact epidemic trajectories and result in unequal distributions of SARS-CoV-2 infection. Funding: K.C.M. was supported by National Science Foundation GRFP grant DGE1745303. Y.H.G. and M.L. were funded by the Morris-Singer Foundation. M.L. was supported by SeroNet cooperative agreement U01 CA261277.

Abstract Background Rates of sepsis/septicemia hospitalization in the US have risen significantly during recent years. Antibiotic resistance and use may contribute to those rates through various mechanisms, including lack of clearance of resistant infections following antibiotic treatment, with some of those infections subsequently devolving into sepsis. At the same time, there is limited information on the effect of prescribing of certain antibiotics vs. others on the rates of septicemia and sepsis-related hospitalizations and mortality. Methods We used multivariable linear regression to relate state-specific rates of outpatient prescribing overall for oral fluoroquinolones, penicillins, macrolides, and cephalosporins between 2011 and 2012 to state-specific rates of septicemia hospitalization (ICD-9 codes 038.xx present anywhere on a discharge diagnosis) in each of the following age groups of adults: (18-49y, 50-64y, 65-74y, 75-84y, 85 + y) reported to the Healthcare Cost and Utilization Project (HCUP) between 2011 and 2012, adjusting for additional covariates, and random effects associated with the ten US Health and Human Services (HHS) regions. Results Increase in the rate of prescribing of oral penicillins by 1 annual dose per 1000 state residents was associated with increases in annual septicemia hospitalization rates of 0.19 (95% CI (0.02,0.37)) per 10,000 persons aged 50-64y, of 0.48(0.12,0.84) per 10,000 persons aged 65-74y, and of 0.81(0.17,1.40) per 10,000 persons aged 74-84y. Increase by 1 in the percent of African Americans among state residents in a given age group was associated with increases in annual septicemia hospitalization rates of 2.3(0.32,4.2) per 10,000 persons aged 75-84y, and of 5.3(1.1,9.5) per 10,000 persons aged over 85y. Average minimal daily temperature was positively associated with septicemia hospitalization rates in persons aged 18-49y, 50-64y, 75-84y and over 85y. Conclusions Our results suggest positive associations between the rates of prescribing for penicillins and the rates of hospitalization with septicemia in US adults aged 50-84y. Further studies are needed to better understand the potential effect of antibiotic replacement in the treatment of various syndromes, including the potential impact of the recent US FDA guidelines on restriction of fluoroquinolone use, as well as the potential effect of changes in the practices for prescribing of penicillins on the rates of sepsis-related hospitalization and mortality.

Abstract Background There is great interest in the use of reduced dosing schedules for pneumococcal conjugate vaccines, a strategy premised on maintaining an acceptable level of protection against disease and carriage of the organism. We asked about the practicality of measuring differential effectiveness against carriage in a population with and without widespread use of the vaccine for infants. Methods We adapted an existing transmission-dynamic, individual-based stochastic model fitted to the prevaccine epidemiology of pneumococcal carriage in the United States, and compared the observed vaccine-type carriage prevalence in different arms of a simulated trial with one, two, or three infant doses plus a 12-month booster. Using these simulations, we calculated vaccine efficacy that would be estimated at different times post-enrollment in the trial and calculated required sample sizes to see a difference in carriage prevalence. Results In a pneumococcal conjugate vaccine (PCV)-naïve population, the difference in vaccine-type (VT) pneumococcal carriage prevalence between trial arms was less than 7% and varied with sampling time. In a population already receiving routine PCV administration, VT pneumococcal prevalence is nearly indistinguishable between trial arms. Relative efficacy of different dosing schedules was strongly dependent on the time between enrollment and sampling, with maximal prevalence differences reached 1–3 years post-enrollment. Moreover, vaccine efficacy estimates were typically slightly higher in trials in communities already receiving vaccination. Despite this, much larger sample sizes—by more than an order of magnitude—are required for a vaccine trial conducted in a population receiving routine PCV administration as compared to in a PCV-naïve population. Conclusions These findings highlight some underappreciated aspects of clinical trials of pneumococcal conjugate vaccines with efficacy endpoints, such as the context- and time-dependence of efficacy estimates. They support the wisdom of conducting comparative dose schedule trials of conjugate vaccine effects on carriage in vaccine-naïve populations.

Estimation of the effective reproductive number Rt is important for detecting changes in disease transmission over time. During the Coronavirus Disease 2019 (COVID-19) pandemic, policy makers and public health officials are using Rt to assess the effectiveness of interventions and to inform policy. However, estimation of Rt from available data presents several challenges, with critical implications for the interpretation of the course of the pandemic. The purpose of this document is to summarize these challenges, illustrate them with examples from synthetic data, and, where possible, make recommendations. For near real-time estimation of Rt, we recommend the approach of Cori and colleagues, which uses data from before time t and empirical estimates of the distribution of time between infections. Methods that require data from after time t, such as Wallinga and Teunis, are conceptually and methodologically less suited for near real-time estimation, but may be appropriate for retrospective analyses of how individuals infected at different time points contributed to the spread. We advise caution when using methods derived from the approach of Bettencourt and Ribeiro, as the resulting Rt estimates may be biased if the underlying structural assumptions are not met. Two key challenges common to all approaches are accurate specification of the generation interval and reconstruction of the time series of new infections from observations occurring long after the moment of transmission. Naive approaches for dealing with observation delays, such as subtracting delays sampled from a distribution, can introduce bias. We provide suggestions for how to mitigate this and other technical challenges and highlight open problems in Rt estimation.

Predicting how pathogen populations will change over time is challenging. Such has been the case with Streptococcus pneumoniae, an important human pathogen, and the pneumococcal conjugate vaccines (PCVs), which target only a fraction of the strains in the population. Here, we use the frequencies of accessory genes to predict changes in the pneumococcal population after vaccination, hypothesizing that these frequencies reflect negative frequency-dependent selection (NFDS) on the gene products. We find that the standardized predicted fitness of a strain, estimated by an NFDS-based model at the time the vaccine is introduced, enables us to predict whether the strain increases or decreases in prevalence following vaccination. Further, we are able to forecast the equilibrium post-vaccine population composition and assess the invasion capacity of emerging lineages. Overall, we provide a method for predicting the impact of an intervention on pneumococcal populations with potential application to other bacterial pathogens in which NFDS is a driving force.

Genotype-based diagnostics for antibiotic resistance represent a promising alternative to empiric therapy, reducing inappropriate antibiotic use. However, because such assays infer resistance based on known genetic markers, their utility will wane with the emergence of novel resistance. Maintenance of these diagnostics will therefore require surveillance to ensure early detection of novel resistance variants, but efficient strategies to do so remain undefined. We evaluate the efficiency of targeted sampling approaches informed by patient and pathogen characteristics in detecting antibiotic resistance and diagnostic escape variants in Neisseria gonorrhoeae, a pathogen associated with a high burden of disease and antibiotic resistance and the development of genotype-based diagnostics. We show that patient characteristic-informed sampling is not a reliable strategy for efficient variant detection. In contrast, sampling informed by pathogen characteristics, such as genomic diversity and genomic background, is significantly more efficient than random sampling in identifying genetic variants associated with resistance and diagnostic escape.

Achieving accurate, real-time estimates of disease activity is challenged by delays in case reporting. "Nowcast" approaches attempt to estimate the complete case counts for a given reporting date, using a time series of case reports that is known to be incomplete due to reporting delays. Modeling the reporting delay distribution is a common feature of nowcast approaches. However, many nowcast approaches ignore a crucial feature of infectious disease transmission-that future cases are intrinsically linked to past reported cases-and are optimized to one or two applications, which may limit generalizability. Here, we present a Bayesian approach, NobBS (Nowcasting by Bayesian Smoothing) capable of producing smooth and accurate nowcasts in multiple disease settings. We test NobBS on dengue in Puerto Rico and influenza-like illness (ILI) in the United States to examine performance and robustness across settings exhibiting a range of common reporting delay characteristics (from stable to time-varying), and compare this approach with a published nowcasting software package while investigating the features of each approach that contribute to good or poor performance. We show that introducing a temporal relationship between cases considerably improves performance when the reporting delay distribution is time-varying, and we identify trade-offs in the role of moving windows to accurately capture changes in the delay. We present software implementing this new approach (R package "NobBS") for widespread application and provide practical guidance on implementation.

ABSTRACT We propose here changes to the U.S. government policy on potential pandemic pathogen (PPP) oversight and implementation, emphasizing transparency of the review process and the content of the review, publication of the review in advance, responsible publication of enhanced PPP research, high-level signoff on approvals of enhanced PPP experiments, and the need for a significant effort to establish a common international approach to enhanced PPP work. We advocate that the U.S. government recommend, and non-U.S. government funders and journals adopt, a set of best practices that would extend important considerations of biosafety and biosecurity to all work on enhanced potential pandemic pathogens regardless of funding source.

The relationship between antibiotic stewardship and population levels of antibiotic resistance remains unclear. In order to better understand shifts in selective pressure due to stewardship, we use publicly available data to estimate the effect of changes in prescribing on exposures to frequently used antibiotics experienced by potentially pathogenic bacteria that are asymptomatically colonizing the microbiome. We quantify this impact under four hypothetical stewardship strategies. In one scenario, we estimate that elimination of all unnecessary outpatient antibiotic use could avert 6% to 48% (IQR: 17% to 31%) of exposures across pairwise combinations of sixteen common antibiotics and nine bacterial pathogens. All scenarios demonstrate that stewardship interventions, facilitated by changes in clinician behavior and improved diagnostics, have the opportunity to broadly reduce antibiotic exposures across a range of potential pathogens. Concurrent approaches, such as vaccines aiming to reduce infection incidence, are needed to further decrease exposures occurring in ‘necessary’ contexts.

The sensitivity of genotype-based diagnostics that predict antimicrobial susceptibility is limited by the extent to which they detect genes and alleles that lead to resistance. As novel resistance variants are expected to emerge, such sensitivity is expected to decline unless the new variants are detected and incorporated into the diagnostic. Here, we present a mathematical framework to define how many diagnostic failures may be expected under varying surveillance regimes and thus quantify the surveillance needed to maintain the sensitivity of genotype-based diagnostics.

Abstract Background Interventions in infectious diseases can have both direct effects on individuals who receive the intervention as well as indirect effects in the population. In addition, intervention combinations can have complex interactions at the population level, which are often difficult to adequately assess with standard study designs and analytical methods. Discussion Herein, we urge the adoption of a new paradigm for the design and interpretation of intervention trials in infectious diseases, particularly with regard to emerging infectious diseases, one that more accurately reflects the dynamics of the transmission process. In an increasingly complex world, simulations can explicitly represent transmission dynamics, which are critical for proper trial design and interpretation. Certain ethical aspects of a trial can also be quantified using simulations. Further, after a trial has been conducted, simulations can be used to explore the possible explanations for the observed effects. Conclusion Much is to be gained through a multidisciplinary approach that builds collaborations among experts in infectious disease dynamics, epidemiology, statistical science, economics, simulation methods, and the conduct of clinical trials.

We investigated the impact of the displacement of sexual activity from adherent recipients of an intervention to others within or outside a trial population on the results from hypothetical trials of different sexual behavior interventions. A short-term model of HIV-prevention interventions that lead to female rejection of male partnership requests showed the impact of displacement expected at the start of a trial. An agent-based model, with sexual mixing and other South African specific demographics, evaluated consequences of displacement for sexual behavior interventions targeting young females in South Africa. This model measured the cumulative incidence among adherent, non-adherent, control and non-enrolled females in a hypothetical trial of HIV prevention. When males made more than one attempt to seek a partnership, interventions reduced short-term HIV infection risk among adherent females, but increased it among non-adherent females as well as controls, non-enrolled (females eligible for the trial but not chosen to participate) and ineligible females (females that did not qualify for the trial due to age). The impact of displacement depends on the intervention and the adherence. In both models, the risk to individuals who are not members of the adherent intervention group will increase with displacement leading to a biased calculation for the effect estimates for the trial. Likewise, intent-to-treat effect estimates become nonlinear functions of the proportion adherent.

Some past epidemics of different influenza subtypes (particularly A/H3N2) in the US saw co-circulation of vaccine-type and variant strains. There is evidence that natural infection with one influenza subtype offers short-term protection against infection with another influenza subtype (henceforth, cross-immunity). This suggests that such cross-immunity for strains within a subtype is expected to be strong. Therefore, while vaccination effective against one strain may reduce transmission of that strain, this may also lead to a reduction of the vaccine-type strain's ability to suppress spread of a variant strain. It remains unclear what the joint effect of vaccination and cross-immunity is for co-circulating influenza strains within a subtype, and what is the potential benefit of a bivalent vaccine that protects against both strains. We simulated co-circulation of vaccine-type and variant strains under a variety of scenarios. In each scenario, we considered the case when the vaccine efficacy against the variant strain is lower than the efficacy against the vaccine-type strain (monovalent vaccine), as well the case when vaccine is equally efficacious against both strains (bivalent vaccine). Administration of a bivalent vaccine results in a significant reduction in the overall incidence of infection compared to administration of a monovalent vaccine, even with lower coverage by the bivalent vaccine. Additionally, we found that with greater cross-immunity, increasing coverage levels for the monovalent vaccine becomes less beneficial, while introducing the bivalent vaccine becomes more beneficial. Our work exhibits the limitations of influenza vaccines that have low efficacy against non-vaccine strains, and demonstrates the benefits of vaccines that offer good protection against multiple influenza strains. The results elucidate the need for guarding against the potential co-circulation of non-vaccine strains for an influenza subtype, at least during select seasons, possibly through inclusion of multiple strains within a subtype (particularly A/H3N2) in a vaccine.

Rapid point-of-care resistance diagnostics (POC-RD) are a key tool in the fight against antibiotic resistance. By tailoring drug choice to infection genotype, doctors can improve treatment efficacy while limiting costs of inappropriate antibiotic prescription. Here, we combine epidemiological theory and data to assess the potential of resistance diagnostics (RD) innovations in a public health context, as a means to limit or even reverse selection for antibiotic resistance. POC-RD can be used to impose a nonbiological fitness cost on resistant strains by enabling diagnostic-informed treatment and targeted interventions that reduce resistant strains' opportunities for transmission. We assess this diagnostic-imposed fitness cost in the context of a spectrum of bacterial population biologies and find that POC-RD have a greater potential against obligate pathogens than opportunistic pathogens already subject to selection under "bystander" antibiotic exposure during asymptomatic carriage (e.g., the pneumococcus). We close by generalizing the notion of RD-informed strategies to incorporate carriage surveillance information and illustrate that coupling transmission-control interventions to the discovery of resistant strains in carriage can potentially select against resistance in a broad range of opportunistic pathogens.

Resistance against different antibiotics appears on the same bacterial strains more often than expected by chance, leading to high frequencies of multidrug resistance. There are multiple explanations for this observation, but these tend to be specific to subsets of antibiotics and/or bacterial species, whereas the trend is pervasive. Here, we consider the question in terms of strain ecology: explaining why resistance to different antibiotics is often seen on the same strain requires an understanding of the competition between strains with different resistance profiles. This work builds on models originally proposed to explain another aspect of strain competition: the stable coexistence of antibiotic sensitivity and resistance observed in a number of bacterial species. We first identify a partial structural similarity in these models: either strain or host population structure stratifies the pathogen population into evolutionarily independent sub-populations and introduces variation in the fitness effect of resistance between these sub-populations, thus creating niches for sensitivity and resistance. We then generalise this unified underlying model to multidrug resistance and show that models with this structure predict high levels of association between resistance to different drugs and high multidrug resistance frequencies. We test predictions from this model in six bacterial datasets and find them to be qualitatively consistent with observed trends. The higher than expected frequencies of multidrug resistance are often interpreted as evidence that these strains are out-competing strains with lower resistance multiplicity. Our work provides an alternative explanation that is compatible with long-term stability in resistance frequencies.

We are writing to reply to the comment by Pouwels et al., 2019 about our recent study (Olesen et al., 2018) on antibiotic use and antibiotic resistance.

BackgroundSeasonal epidemics of bacterial meningitis in the African Meningitis Belt carry a high burden of disease and mortality. Reactive mass vaccination is used as a control measure during epidemics, but the time taken to gain immunity from the vaccine reduces the flexibility and effectiveness of these campaigns. Targeted reactive antibiotic prophylaxis could be used to supplement reactive mass vaccination and further reduce the incidence of meningitis, and the potential effectiveness and efficiency of these strategies should be explored.Methods and findingsData from an outbreak of meningococcal meningitis in Niger, caused primarily by Neisseria meningitidis serogroup C, is used to estimate clustering of meningitis cases at the household and village level. In addition, reactive antibiotic prophylaxis and reactive vaccination strategies are simulated to estimate their potential effectiveness and efficiency, with a focus on the threshold and spatial unit used to declare an epidemic and initiate the intervention. There is village-level clustering of suspected meningitis cases after an epidemic has been declared in a health area. Risk of suspected meningitis among household contacts of a suspected meningitis case is no higher than among members of the same village. Village-wide antibiotic prophylaxis can target subsequent cases in villages: across of range of parameters pertaining to how the intervention is performed, up to 220/672 suspected cases during the season are potentially preventable. On the other hand, household prophylaxis targets very few cases. In general, the village-wide strategy is not very sensitive to the method used to declare an epidemic. Finally, village-wide antibiotic prophylaxis is potentially more efficient than mass vaccination of all individuals at the beginning of the season, and than the equivalent reactive vaccination strategy.ConclusionsVillage-wide antibiotic prophylaxis should be considered and tested further as a response against outbreaks of meningococcal meningitis in the Meningitis Belt, as a supplement to reactive mass vaccination.

Pneumococcal whole cell vaccines (WCVs) could cost-effectively protect against a greater strain diversity than current capsule-based vaccines. Immunoglobulin G (IgG) responses to a WCV were characterised by applying longitudinally-sampled sera, available from 35 adult placebo-controlled phase I trial participants, to a panproteome microarray. Despite individuals maintaining distinctive antibody ‘fingerprints’, responses were consistent across vaccinated cohorts. Seventy-two functionally distinct proteins were associated with WCV-induced increases in IgG binding. These shared characteristics with naturally immunogenic proteins, being enriched for transporters and cell wall metabolism enzymes, likely unusually exposed on the unencapsulated WCV’s surface. Vaccine-induced responses were specific to variants of the diverse PclA, PspC and ZmpB proteins, whereas PspA- and ZmpA-induced antibodies recognised a broader set of alleles. Temporal variation in IgG levels suggested a mixture of anamnestic and novel responses. These reproducible increases in IgG binding to a limited, but functionally diverse, set of conserved proteins indicate WCV could provide species-wide immunity. Clinical trial registration: The trial was registered with ClinicalTrials.gov with Identifier NCT01537185; the results are available from https://clinicaltrials.gov/ct2/show/results/NCT01537185.

Antibiotic use is a primary driver of antibiotic resistance. However, antibiotic use can be distributed in different ways in a population, and the association between the distribution of use and antibiotic resistance has not been explored. Here, we tested the hypothesis that repeated use of antibiotics has a stronger association with population-wide antibiotic resistance than broadly-distributed, low-intensity use. First, we characterized the distribution of outpatient antibiotic use across US states, finding that antibiotic use is uneven and that repeated use of antibiotics makes up a minority of antibiotic use. Second, we compared antibiotic use with resistance for 72 pathogen-antibiotic combinations across states. Finally, having partitioned total use into extensive and intensive margins, we found that intense use had a weaker association with resistance than extensive use. If the use-resistance relationship is causal, these results suggest that reducing total use and selection intensity will require reducing broadly distributed, low-intensity use.

The identification of key “driver” groups in influenza epidemics is of much interest for the implementation of effective public health response strategies, including vaccination programs. However, the relative importance of different age groups in propagating epidemics is uncertain. During a communicable disease outbreak, some groups may be disproportionately represented during the outbreak's ascent due to increased susceptibility and/or contact rates. Such groups or subpopulations can be identified by considering the proportion of cases within the subpopulation occurring before (Bp) and after the epidemic peak (Ap) to calculate the subpopulation's relative risk, RR =Bp/Ap. We estimated RR for several subpopulations (age groups) using data on laboratory-confirmed US influenza hospitalizations during epidemics between 2009 and 2014. Additionally, we simulated various influenza outbreaks in an age-stratified population, relating the RR to the impact of vaccination in each subpopulation on the epidemic's initial effective reproductive number Re(0). We found that children aged 5–17 had the highest estimates of RR during the five largest influenza A outbreaks, though the relative magnitude of RR in this age group compared to other age groups varied, being highest for the 2009 A/H1N1 pandemic. For the 2010–2011 and 2012–2013 influenza B epidemics, adults aged 18–49, and 0–4 year-olds had the highest estimates of RR, respectively. For 83% of simulated epidemics, the group with the highest RR was also the group for which initial distribution of a given quantity of vaccine would result in the largest reduction of Re(0). In the largest 40% of simulated outbreaks, the group with the highest RR and the largest vaccination impact was children 5–17. While the relative importance of different age groups in propagating influenza outbreaks varies, children aged 5–17 play the leading role during the largest influenza A epidemics. Extra vaccination efforts for this group may contribute to reducing the epidemic's impact in the whole community.

For encapsulated bacteria such as Streptococcus pneumoniae, asymptomatic carriage is more common and longer in duration than disease, and hence is often a more convenient endpoint for clinical trials of vaccines against these bacteria. However, using a carriage endpoint entails specific challenges. Carriage is almost always measured as prevalence, whereas the vaccine may act by reducing incidence or duration. Thus, to determine sample size requirements, its impact on prevalence must first be estimated. The relationship between incidence and prevalence (or duration and prevalence) is convex, saturating at 100% prevalence. For this reason, the proportional effect of a vaccine on prevalence is typically less than its proportional effect on incidence or duration. This relationship is further complicated in the presence of multiple pathogen strains. In addition, host immunity to carriage accumulates rapidly with frequent exposures in early years of life, creating potentially complex interactions with the vaccine's effect. We conducted a simulation study to predict the impact of an inactivated whole cell pneumococcal vaccine-believed to reduce carriage duration-on carriage prevalence in different age groups and trial settings. We used an individual-based model of pneumococcal carriage that incorporates relevant immunological processes, both vaccine-induced and naturally acquired. Our simulations showed that for a wide range of vaccine efficacies, sampling time and age at vaccination are important determinants of sample size. There is a window of favorable sampling times during which the required sample size is relatively low, and this window is prolonged with a younger age at vaccination, and in a trial setting with lower transmission intensity. These results illustrate the ability of simulation studies to inform the planning of vaccine trials with carriage endpoints, and the methods we present here can be applied to trials evaluating other pneumococcal vaccine candidates or comparing alternative dosing schedules for the existing conjugate vaccines.

In a Policy Forum, Marc Lipsitch and colleagues discuss trial design issues in infectious disease outbreaks.

In an Essay, Michael Johansson and colleagues advocate the posting of research studies addressing infectious disease outbreaks as preprints.

In the United States, the introduction of the heptavalent pneumococcal conjugate vaccine (PCV) largely eliminated vaccine serotypes (VT); non-vaccine serotypes (NVT) subsequently increased in carriage and disease. Vaccination also disrupts the composition of the pneumococcal pangenome, which includes mobile genetic elements and polymorphic non-capsular antigens important for virulence, transmission, and pneumococcal ecology. Antigenic proteins are of interest for future vaccines; yet, little is known about how the they are affected by PCV use. To investigate the evolutionary impact of vaccination, we assessed recombination, evolution, and pathogen demographic history of 937 pneumococci collected from 1998-2012 among Navajo and White Mountain Apache Native American communities. We analyzed changes in the pneumococcal pangenome, focusing on metabolic loci and 19 polymorphic protein antigens. We found the impact of PCV on the pneumococcal population could be observed in reduced diversity, a smaller pangenome, and changing frequencies of accessory clusters of orthologous groups (COGs). Post-PCV7, diversity rebounded through clonal expansion of NVT lineages and inferred in-migration of two previously unobserved lineages. Accessory COGs frequencies trended toward pre-PCV7 values with increasing time since vaccine introduction. Contemporary frequencies of protein antigen variants are better predicted by pre-PCV7 values (1998-2000) than the preceding period (2006-2008), suggesting balancing selection may have acted in maintaining variant frequencies in this population. Overall, we present the largest genomic analysis of pneumococcal carriage in the United States to date, which includes a snapshot of a true vaccine-naïve community prior to the introduction of PCV7. These data improve our understanding of pneumococcal evolution and emphasize the need to consider pangenome composition when inferring the impact of vaccination and developing future protein-based pneumococcal vaccines.

ABSTRACT Heterogeneity in host susceptibility is a key determinant of infectious disease dynamics but is rarely accounted for in assessment of disease control measures. Understanding how susceptibility is distributed in populations, and how control measures change this distribution, is integral to predicting the course of epidemics with and without interventions. Using multiple experimental and modeling approaches, we show that rainbow trout have relatively homogeneous susceptibility to infection with infectious hematopoietic necrosis virus and that vaccination increases heterogeneity in susceptibility in a nearly all-or-nothing fashion. In a simple transmission model with an R0 of 2, the highly heterogeneous vaccine protection would cause a 35 percentage-point reduction in outbreak size over an intervention inducing homogenous protection at the same mean level. More broadly, these findings provide validation of methodology that can help to reduce biases in predictions of vaccine impact in natural settings and provide insight into how vaccination shapes population susceptibility. IMPORTANCE Differences among individuals influence transmission and spread of infectious diseases as well as the effectiveness of control measures. Control measures, such as vaccines, may provide leaky protection, protecting all hosts to an identical degree, or all-or-nothing protection, protecting some hosts completely while leaving others completely unprotected. This distinction can have a dramatic influence on disease dynamics, yet this distribution of protection is frequently unaccounted for in epidemiological models and estimates of vaccine efficacy. Here, we apply new methodology to experimentally examine host heterogeneity in susceptibility and mode of vaccine action as distinct components influencing disease outcome. Through multiple experiments and new modeling approaches, we show that the distribution of vaccine effects can be robustly estimated. These results offer new experimental and inferential methodology that can improve predictions of vaccine effectiveness and have broad applicability to human, wildlife, and ecosystem health.

BACKGROUND:The 2014-6 West African Ebola epidemic highlights the need for rigorous, rapid clinical trial methods for vaccines. A challenge for trial design is making sample size calculations based on incidence within the trial, total vaccine effect, and intracluster correlation, when these parameters are uncertain in the presence of indirect effects of vaccination. METHODS AND FINDINGS:We present a stochastic, compartmental model for a ring vaccination trial. After identification of an index case, a ring of contacts is recruited and either vaccinated immediately or after 21 days. The primary outcome of the trial is total vaccine effect, counting cases only from a pre-specified window in which the immediate arm is assumed to be fully protected and the delayed arm is not protected. Simulation results are used to calculate necessary sample size and estimated vaccine effect. Under baseline assumptions about vaccine properties, monthly incidence in unvaccinated rings and trial design, a standard sample-size calculation neglecting dynamic effects estimated that 7,100 participants would be needed to achieve 80% power to detect a difference in attack rate between arms, while incorporating dynamic considerations in the model increased the estimate to 8,900. This approach replaces assumptions about parameters at the ring level with assumptions about disease dynamics and vaccine characteristics at the individual level, so within this framework we were able to describe the sensitivity of the trial power and estimated effect to various parameters. We found that both of these quantities are sensitive to properties of the vaccine, to setting-specific parameters over which investigators have little control, and to parameters that are determined by the study design. CONCLUSIONS:Incorporating simulation into the trial design process can improve robustness of sample size calculations. For this specific trial design, vaccine effectiveness depends on properties of the ring vaccination design and on the measurement window, as well as the epidemiologic setting.

The threat of an influenza A virus pandemic stems from continual virus spillovers from reservoir species, a tiny fraction of which spark sustained transmission in humans. To date, no pandemic emergence of a new influenza strain has been preceded by detection of a closely related precursor in an animal or human. Nonetheless, influenza surveillance efforts are expanding, prompting a need for tools to assess the pandemic risk posed by a detected virus. The goal would be to use genetic sequence and/or biological assays of viral traits to identify those non-human influenza viruses with the greatest risk of evolving into pandemic threats, and/or to understand drivers of such evolution, to prioritize pandemic prevention or response measures. We describe such efforts, identify progress and ongoing challenges, and discuss three specific traits of influenza viruses (hemagglutinin receptor binding specificity, hemagglutinin pH of activation, and polymerase complex efficiency) that contribute to pandemic risk.

Animal studies aimed at understanding influenza virus mutations that change host specificity to adapt to replication in mammalian hosts are necessarily limited in sample numbers due to high cost and safety requirements. As a safe, higher-throughput alternative, we explore the possibility of using readily available passage bias data obtained mostly from seasonal H1 and H3 influenza strains that were differentially grown in mammalian (MDCK) and avian cells (eggs). Using a statistical approach over 80,000 influenza hemagglutinin sequences with passage information, we found that passage bias sites are most commonly found in three regions: (i) the globular head domain around the receptor binding site, (ii) the region that undergoes pH-dependent structural changes and (iii) the unstructured N-terminal region harbouring the signal peptide. Passage bias sites were consistent among different passage cell types as well as between influenza A subtypes. We also find epistatic interactions of site pairs supporting the notion of host-specific dependency of mutations on virus genomic background. The sites identified from our large-scale sequence analysis substantially overlap with known host adaptation sites in the WHO H5N1 genetic changes inventory suggesting information from passage bias can provide candidate sites for host specificity changes to aid in risk assessment for emerging strains.