Quantifying the information in noisy epidemic curves

Machine learning mathematical models for incidence estimation during pandemics

Accurate estimates of the incidence of infectious diseases are key for the control of epidemics. However, healthcare systems are often unable to test the population exhaustively, especially when asymptomatic and paucisymptomatic cases are widespread; this leads to significant and systematic under-reporting of the real incidence. Here, we propose a machine learning approach to estimate the incidence of a pandemic in real-time, using reported cases and the overall test rate. In particular, we use Bayesian symbolic regression to automatically learn the closed-form mathematical models that most parsimoniously describe incidence. We develop and validate our models using COVID-19 incidence values for nine different countries, confirming their ability to accurately predict daily incidence. Remarkably, despite the differences in epidemic trajectories and dynamics across countries, we find that a single model for all countries offers a more parsimonious description and is more predictive of actual incidence compared to separate models for each country. Our results show the potential to accurately model incidence in real-time using closed-form mathematical models, providing a valuable tool for public health decision-makers.

Jointly estimating epidemiological dynamics of Covid-19 from case and wastewater data in Aotearoa New Zealand

BACKGROUND

Timely and informed public health responses to infectious diseases such as COVID-19 necessitate reliable information about infection dynamics. The case ascertainment rate (CAR), the proportion of infections that are reported as cases, is typically much less than one and varies with testing practices and behaviours, making reported cases unreliable as the sole source of data. The concentration of viral RNA in wastewater samples provides an alternate measure of infection prevalence that is not affected by clinical testing, healthcare-seeking behaviour or access to care.

METHODS

We construct a state-space model with observed data of levels of SARS-CoV-2 in wastewater and reported case incidence and estimate the hidden states of the effective reproduction number, R, and CAR using sequential Monte Carlo methods.

RESULTS

We analyse data from 1 January 2022 to 31 March 2023 from Aotearoa New Zealand. Our model estimates that R peaks at 2.76 (95% CrI 2.20, 3.83) around 18 February 2022 and the CAR peaks around 12 March 2022. We calculate that New Zealand's second Omicron wave in July 2022 is similar in size to the first, despite fewer reported cases. We estimate that the CAR in the BA.5 Omicron wave in July 2022 is approximately 50% lower than in the BA.1/BA.2 Omicron wave in March 2022.

CONCLUSIONS

Estimating R, CAR, and cumulative number of infections provides useful information for planning public health responses and understanding the state of immunity in the population. This model is a useful disease surveillance tool, improving situational awareness of infectious disease dynamics in real-time.

Host behaviour driven by awareness of infection risk amplifies the chance of superspreading events

We demonstrate that heterogeneity in the perceived risks associated with infection within host populations amplifies chances of superspreading during the crucial early stages of epidemics. Under this behavioural model, individuals less concerned about dangers from infection are more likely to be infected and attend larger sized (riskier) events, where we assume event sizes remain unchanged. For directly transmitted diseases such as COVID-19, this leads to infections being introduced at rates above the population prevalence to those events most conducive to superspreading. We develop an interpretable, computational framework for evaluating within-event risks and derive a small-scale reproduction number measuring how the infections generated at an event depend on transmission heterogeneities and numbers of introductions. This generalizes previous frameworks and quantifies how event-scale patterns and population-level characteristics relate. As event duration and size grow, our reproduction number converges to the basic reproduction number. We illustrate that even moderate levels of heterogeneity in the perceived risks of infection substantially increase the likelihood of disproportionately large clusters of infections occurring at larger events, despite fixed overall disease prevalence. We show why collecting data linking host behaviour and event attendance is essential for accurately assessing the risks posed by invading pathogens in emerging stages of outbreaks.

Generative Bayesian modeling to nowcast the effective reproduction number from line list data with missing symptom onset dates

The time-varying effective reproduction number Rt is a widely used indicator of transmission dynamics during infectious disease outbreaks. Timely estimates of Rt can be obtained from reported cases counted by their date of symptom onset, which is generally closer to the time of infection than the date of report. Case counts by date of symptom onset are typically obtained from line list data, however these data can have missing information and are subject to right truncation. Previous methods have addressed these problems independently by first imputing missing onset dates, then adjusting truncated case counts, and finally estimating the effective reproduction number. This stepwise approach makes it difficult to propagate uncertainty and can introduce subtle biases during real-time estimation due to the continued impact of assumptions made in previous steps. In this work, we integrate imputation, truncation adjustment, and Rt estimation into a single generative Bayesian model, allowing direct joint inference of case counts and Rt from line list data with missing symptom onset dates. We then use this framework to compare the performance of nowcasting approaches with different stepwise and generative components on synthetic line list data for multiple outbreak scenarios and across different epidemic phases. We find that under reporting delays realistic for hospitalization data (50% of reports delayed by more than a week), intermediate smoothing, as is common practice in stepwise approaches, can bias nowcasts of case counts and Rt, which is avoided in a joint generative approach due to shared regularization of all model components. On incomplete line list data, a fully generative approach enables the quantification of uncertainty due to missing onset dates without the need for an initial multiple imputation step. In a real-world comparison using hospitalization line list data from the COVID-19 pandemic in Switzerland, we observe the same qualitative differences between approaches. The generative modeling components developed in this work have been integrated and further extended in the R package epinowcast, providing a flexible and interpretable tool for real-time surveillance.

Estimation of the Time-Varying Effective Reproductive Number of COVID-19 Based on Multivariate Time Series of Severe Health Outcomes

The time-varying effective reproduction number (Rt at time t) measures the transmissibility of SARS-CoV-2 and is conventionally based on daily case counts, which may suffer from time-varying ascertainment. We analyzed Rt estimates from case counts and severe COVID-19 (intensive care unit admissions, severe or critical cases, and mortality) across 2022 in Hong Kong's fifth and sixth waves of infection. Within the fifth wave, the severe disease-based Rt (3.5) was significantly higher than the case-based Rt (2.4) but not in the sixth wave. During periods with fluctuating underreporting, data based on severe diseases may provide more reliable Rt estimates.

Angular reproduction numbers improve estimates of transmissibility when disease generation times are misspecified or time-varying

We introduce the angular reproduction number Ω, which measures time-varying changes in epidemic transmissibility resulting from variations in both the effective reproduction number R, and generation time distribution w. Predominant approaches for tracking pathogen spread infer either R or the epidemic growth rate r. However, R is biased by mismatches between the assumed and true w, while r is difficult to interpret in terms of the individual-level branching process underpinning transmission. R and r may also disagree on the relative transmissibility of epidemics or variants (i.e. rA > rB does not imply RA > RB for variants A and B). We find that Ω responds meaningfully to mismatches and time-variations in w while mostly maintaining the interpretability of R. We prove that Ω > 1 implies R > 1 and that Ω agrees with r on the relative transmissibility of pathogens. Estimating Ω is no more difficult than inferring R, uses existing software, and requires no generation time measurements. These advantages come at the expense of selecting one free parameter. We propose Ω as complementary statistic to R and r that improves transmissibility estimates when w is misspecified or time-varying and better reflects the impact of interventions, when those interventions concurrently change R and w or alter the relative risk of co-circulating pathogens.

Risk averse reproduction numbers improve resurgence detection

The effective reproduction number R is a prominent statistic for inferring the transmissibility of infectious diseases and effectiveness of interventions. R purportedly provides an easy-to-interpret threshold for deducing whether an epidemic will grow (R>1) or decline (R<1). We posit that this interpretation can be misleading and statistically overconfident when applied to infections accumulated from groups featuring heterogeneous dynamics. These groups may be delineated by geography, infectiousness or sociodemographic factors. In these settings, R implicitly weights the dynamics of the groups by their number of circulating infections. We find that this weighting can cause delayed detection of outbreak resurgence and premature signalling of epidemic control because it underrepresents the risks from highly transmissible groups. Applying E-optimal experimental design theory, we develop a weighting algorithm to minimise these issues, yielding the risk averse reproduction number E. Using simulations, analytic approaches and real-world COVID-19 data stratified at the city and district level, we show that E meaningfully summarises transmission dynamics across groups, balancing bias from the averaging underlying R with variance from directly using local group estimates. An E>1generates timely resurgence signals (upweighting risky groups), while an E<1ensures local outbreaks are under control. We propose E as an alternative to R for informing policy and assessing transmissibility at large scales (e.g., state-wide or nationally), where R is commonly computed but well-mixed or homogeneity assumptions break down.

Resurgence of Omicron BA.2 in SARS-CoV-2 infection-naive Hong Kong

Hong Kong experienced a surge of Omicron BA.2 infections in early 2022, resulting in one of the highest per-capita death rates of COVID-19. The outbreak occurred in a dense population with low immunity towards natural SARS-CoV-2 infection, high vaccine hesitancy in vulnerable populations, comprehensive disease surveillance and the capacity for stringent public health and social measures (PHSMs). By analyzing genome sequences and epidemiological data, we reconstructed the epidemic trajectory of BA.2 wave and found that the initial BA.2 community transmission emerged from cross-infection within hotel quarantine. The rapid implementation of PHSMs suppressed early epidemic growth but the effective reproduction number (Re) increased again during the Spring festival in early February and remained around 1 until early April. Independent estimates of point prevalence and incidence using phylodynamics also showed extensive superspreading at this time, which likely contributed to the rapid expansion of the epidemic. Discordant inferences based on genomic and epidemiological data underscore the need for research to improve near real-time epidemic growth estimates by combining multiple disparate data sources to better inform outbreak response policy.

Effectiveness assessment of non-pharmaceutical interventions: lessons learned from the COVID-19 pandemic

Effectiveness of non-pharmaceutical interventions (NPIs), such as school closures and stay-at-home orders, during the COVID-19 pandemic has been assessed in many studies. Such assessments can inform public health policies and contribute to evidence-based choices of NPIs during subsequent waves or future epidemics. However, methodological issues and no standardised assessment practices have restricted the practical value of the existing evidence. Here, we present and discuss lessons learned from the COVID-19 pandemic and make recommendations for standardising and improving assessment, data collection, and modelling. These recommendations could contribute to reliable and policy-relevant assessments of the effectiveness of NPIs during future epidemics.