Practical considerations for measuring the effective reproductive number, Rt

Utilizing disease transmission and response capacities to optimize covid-19 control in Malaysia

OBJECTIVES

Public Health Social Measures (PHSM) such as movement restriction movement needed to be adjusted accordingly during the COVID-19 pandemic to ensure low disease transmission alongside adequate health system capacities based on the COVID-19 situational matrix proposed by the World Health Organization (WHO). This paper aims to develop a mechanism to determine the COVID-19 situational matrix to adjust movement restriction intensity for the control of COVID-19 in Malaysia.

METHODS

Several epidemiological indicators were selected based on the WHO PHSM interim guidance report and validated individually and in several combinations to estimate the community transmission level (CT) and health system response capacity (RC) variables. Correlation analysis between CT and RC with COVID-19 cases was performed to determine the most appropriate CT and RC variables. Subsequently, the CT and RC variables were combined to form a composite COVID-19 situational matrix (SL). The SL matrix was validated using correlation analysis with COVID-19 case trends. Subsequently, an automated web-based system that generated daily CT, RC, and SL was developed.

RESULTS

CT and RC variables were estimated using case incidence and hospitalization rate; Hospital bed capacity and COVID-19 ICU occupancy respectively. The estimated CT and RC were strongly correlated [ρ = 0.806 (95% CI 0.752, 0.848); and ρ = 0.814 (95% CI 0.778, 0.839), p < 0.001] with the COVID-19 cases. The estimated SL was strongly correlated with COVID-19 cases (ρ = 0.845, p < 0.001) and responded well to the various COVID-19 case trends during the pandemic. SL changes occurred earlier during the increase of cases but slower during the decrease, indicating a conservative response. The automated web-based system developed produced daily real-time CT, RC, and SL for the COVID-19 pandemic.

CONCLUSIONS

The indicators selected and combinations formed were able to generate validated daily CT and RC levels for Malaysia. Subsequently, the CT and RC levels were able to provide accurate and sensitive information for the estimation of SL which provided valuable evidence on the progression of the pandemic and movement restriction adjustment for the control of Malaysia.

An Efficient Approach to Nowcasting the Time-varying Reproduction Number

Estimating the instantaneous reproduction number () in near real time is crucial for monitoring and responding to epidemic outbreaks on a daily basis. However, such estimates often suffer from bias due to reporting delays inherent in surveillance systems. We propose a fast and flexible Bayesian methodology to overcome this challenge by estimating while taking into account reporting delays. Furthermore, the method naturally takes into account the uncertainty associated with the nowcasting of cases to get a valid uncertainty estimation of the nowcasted reproduction number. We evaluate the proposed methodology through a simulation study and apply it to COVID-19 incidence data in Belgium.

Impacts of human mobility on the citywide transmission dynamics of 18 respiratory viruses in pre- and post-COVID-19 pandemic years

Many studies have used mobile device location data to model SARS-CoV-2 dynamics, yet relationships between mobility behavior and endemic respiratory pathogens are less understood. We studied the effects of population mobility on the transmission of 17 endemic viruses and SARS-CoV-2 in Seattle over a 4-year period, 2018-2022. Before 2020, visits to schools and daycares, within-city mixing, and visitor inflow preceded or coincided with seasonal outbreaks of endemic viruses. Pathogen circulation dropped substantially after the initiation of COVID-19 stay-at-home orders in March 2020. During this period, mobility was a positive, leading indicator of transmission of all endemic viruses and lagging and negatively correlated with SARS-CoV-2 activity. Mobility was briefly predictive of SARS-CoV-2 transmission when restrictions relaxed but associations weakened in subsequent waves. The rebound of endemic viruses was heterogeneously timed but exhibited stronger, longer-lasting relationships with mobility than SARS-CoV-2. Overall, mobility is most predictive of respiratory virus transmission during periods of dramatic behavioral change and at the beginning of epidemic waves.

A simulation-based approach for estimating the time-dependent reproduction number from temporally aggregated disease incidence time series data

Tracking pathogen transmissibility during infectious disease outbreaks is essential for assessing the effectiveness of public health measures and planning future control strategies. A key measure of transmissibility is the time-dependent reproduction number, which has been estimated in real-time during outbreaks of a range of pathogens from disease incidence time series data. While commonly used approaches for estimating the time-dependent reproduction number can be reliable when disease incidence is recorded frequently, such incidence data are often aggregated temporally (for example, numbers of cases may be reported weekly rather than daily). As we show, commonly used methods for estimating transmissibility can be unreliable when the timescale of transmission is shorter than the timescale of data recording. To address this, here we develop a simulation-based approach involving Approximate Bayesian Computation for estimating the time-dependent reproduction number from temporally aggregated disease incidence time series data. We first use a simulated dataset representative of a situation in which daily disease incidence data are unavailable and only weekly summary values are reported, demonstrating that our method provides accurate estimates of the time-dependent reproduction number under such circumstances. We then apply our method to two outbreak datasets consisting of weekly influenza case numbers in 2019-20 and 2022-23 in Wales (in the United Kingdom). Our simple-to-use approach will allow accurate estimates of time-dependent reproduction numbers to be obtained from temporally aggregated data during future infectious disease outbreaks.

SARS-CoV-2 Transmission in Alberta, British Columbia, and Ontario, Canada, January 2020-January 2022

We estimated COVID-19 transmission potential and case burden by variant type in Alberta, British Columbia, and Ontario, Canada, during January 23, 2020-January 27, 2022; we also estimated the effectiveness of public health interventions to reduce transmission. We estimated time-varying reproduction number (Rt) over 7-day sliding windows and nonoverlapping time-windows determined by timing of policy changes. We calculated incidence rate ratios (IRRs) for each variant and compared rates to determine differences in burden among provinces. Rt corresponding with emergence of the Delta variant increased in all 3 provinces; British Columbia had the largest increase, 43.85% (95% credible interval [CrI] 40.71%-46.84%). Across the study period, IRR was highest for Omicron (8.74 [95% CrI 8.71-8.77]) and burden highest in Alberta (IRR 1.80 [95% CrI 1.79-1.81]). Initiating public health interventions was associated with lower Rt and relaxing restrictions and emergence of new variants associated with increases in Rt.

Respiratory disease contact patterns in the US are stable but heterogeneous

Background

Contact plays a critical role in infectious disease transmission. Characterizing heterogeneity in contact patterns across individuals, time, and space is necessary to inform accurate estimates of transmission risk, particularly to explain superspreading, predict age differences in vulnerability, and inform social distancing policies. Current respiratory disease models often rely on data from the 2008 POLYMOD study conducted in Europe, which is now outdated and potentially unrepresentative of behavior in the US. We seek to understand the variation in contact patterns across spatial scales and demographic and social classifications, whether there is seasonality to contact patterns, and what social behavior looks like at baseline in the absence of an ongoing pandemic.

Methods

We analyze spatiotemporal non-household contact patterns across 11 million survey responses from June 2020 - April 2021 post-stratified on age and gender to correct for sample representation. To characterize spatiotemporal heterogeneity in respiratory contact patterns at the county-week scale, we use generalized additive models. In the absence of pre-pandemic data on contact in the US, we also use a regression approach to produce baseline contact estimates to fill this gap.

Findings

Although contact patterns varied over time during the pandemic, contact is relatively stable after controlling for disease. We find that the mean number of non-household contacts is spatially heterogeneous regardless of disease. There is additional heterogeneity across age, gender, race/ethnicity, and contact setting, with mean contact decreasing with age and lower in women. The contacts of white individuals and contacts at work or social events change the most under increased national incidence.

Interpretation

We develop the first county-level estimates of non-pandemic contact rates for the US that can fill critical gaps in parameterizing disease models. Our results identify that spatiotemporal, demographic, and social heterogeneity in contact patterns is highly structured, informing the risk landscape of respiratory disease transmission in the US.

MODELS: a six-step framework for developing an infectious disease model

Since the COVID-19 pandemic began, a plethora of modeling studies related to COVID-19 have been released. While some models stand out due to their innovative approaches, others are flawed in their methodology. To assist novices, frontline healthcare workers, and public health policymakers in navigating the complex landscape of these models, we introduced a structured framework named MODELS. This framework is designed to detail the essential steps and considerations for creating a dependable epidemic model, offering direction to researchers engaged in epidemic modeling endeavors.

Estimating the reproduction number and transmission heterogeneity from the size distribution of clusters of identical pathogen sequences

Quantifying transmission intensity and heterogeneity is crucial to ascertain the threat posed by infectious diseases and inform the design of interventions. Methods that jointly estimate the reproduction number R and the dispersion parameter k have however mainly remained limited to the analysis of epidemiological clusters or contact tracing data, whose collection often proves difficult. Here, we show that clusters of identical sequences are imprinted by the pathogen offspring distribution, and we derive an analytical formula for the distribution of the size of these clusters. We develop and evaluate an inference framework to jointly estimate the reproduction number and the dispersion parameter from the size distribution of clusters of identical sequences. We then illustrate its application across a range of epidemiological situations. Finally, we develop a hypothesis testing framework relying on clusters of identical sequences to determine whether a given pathogen genetic subpopulation is associated with increased or reduced transmissibility. Our work provides tools to estimate the reproduction number and transmission heterogeneity from pathogen sequences without building a phylogenetic tree, thus making it easily scalable to large pathogen genome datasets.

Assessing the Time Evolution of COVID-19 Effective Reproduction Number in Brazil

In this paper, we use a Bayesian method to estimate the effective reproduction number ( R ( t ) ), in the context of monitoring the time evolution of the COVID-19 pandemic in Brazil at different geographic levels. The focus of this study is to investigate the similarities between the trends in the evolution of such indicators at different subnational levels with the trends observed nationally. The underlying question addressed is whether national surveillance of such variables is enough to provide a picture of the epidemic evolution in the country or if it may hide important localized trends. This is particularly relevant in the scenario where health authorities use information obtained from such indicators in the design of non-pharmaceutical intervention policies to control the epidemic. A comparison between R ( t ) estimates and the moving average (MA) of daily reported infections is also presented, which is another commonly monitored variable. The analysis carried out in this paper is based on the data of confirmed infected cases provided by a public repository. The correlations between the time series of R ( t ) and MA in different geographic levels are assessed. Comparing national with subnational trends, higher degrees of correlation are found for the time series of R ( t ) estimates, compared to the MA time series. Nevertheless, differences between national and subnational trends are observed for both indicators, suggesting that local epidemiological surveillance would be more suitable as an input to the design of non-pharmaceutical intervention policies in Brazil, particularly for the least populated states.

Incorporating testing volume into estimation of effective reproduction number dynamics

Branching process inspired models are widely used to estimate the effective reproduction number-a useful summary statistic describing an infectious disease outbreak-using counts of new cases. Case data is a real-time indicator of changes in the reproduction number, but is challenging to work with because cases fluctuate due to factors unrelated to the number of new infections. We develop a new model that incorporates the number of diagnostic tests as a surveillance model covariate. Using simulated data and data from the SARS-CoV-2 pandemic in California, we demonstrate that incorporating tests leads to improved performance over the state of the art.

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.

AH, Fernandez S, Huang AT, Vatanaprasan P, Saethang T, Luk-in S, Storer RJ, Ounjai P, Devanga Ragupathi NK, Kanthawee P, Ngamwongsatit N, Badavath VN, Thuptimdang W, Leelahavanichkul A, Kanjanabuch T, Miyanaga K, Cui L, Nanbo A, Shibuya K, Kupwiwat R, Sano D, Furukawa T, Sei K, Higgins PG, Kicic A, Singer AC, Chatsuwan T, Trowsdale S, Abe S, Ishikawa H, Amarasiri M, Modchang C, Wannigama DL. Exploring indoor and outdoor dust as a potential tool for detection and monitoring of COVID-19 transmission

This study investigated the potential of using SARS-CoV-2 viral concentrations in dust as an additional surveillance tool for early detection and monitoring of COVID-19 transmission. Dust samples were collected from 8 public locations in 16 districts of Bangkok, Thailand, from June to August 2021. SARS-CoV-2 RNA concentrations in dust were quantified, and their correlation with community case incidence was assessed. Our findings revealed a positive correlation between viral concentrations detected in dust and the relative risk of COVID-19. The highest risk was observed with no delay (0-day lag), and this risk gradually decreased as the lag time increased. We observed an overall decline in viral concentrations in public places during lockdown, closely associated with reduced human mobility. The effective reproduction number for COVID-19 transmission remained above one throughout the study period, suggesting that transmission may persist in locations beyond public areas even after the lockdown measures were in place.

Underdetected dispersal and extensive local transmission drove the 2022 mpox epidemic

The World Health Organization declared mpox a public health emergency of international concern in July 2022. To investigate global mpox transmission and population-level changes associated with controlling spread, we built phylogeographic and phylodynamic models to analyze MPXV genomes from five global regions together with air traffic and epidemiological data. Our models reveal community transmission prior to detection, changes in case reporting throughout the epidemic, and a large degree of transmission heterogeneity. We find that viral introductions played a limited role in prolonging spread after initial dissemination, suggesting that travel bans would have had only a minor impact. We find that mpox transmission in North America began declining before more than 10% of high-risk individuals in the USA had vaccine-induced immunity. Our findings highlight the importance of broader routine specimen screening surveillance for emerging infectious diseases and of joint integration of genomic and epidemiological information for early outbreak control.

Impact of public health policy and mobility change on transmission potential of severe acute respiratory syndrome coronavirus 2 in Rhode Island, March 2020 - November 2021

To study the SARS-CoV-2 transmission potential in Rhode Island (RI) and its association with policy changes and mobility changes, the time-varying reproduction number, Rt, was estimated. The daily incident case counts (16 March 2020, through 30 November 2021) were bootstrapped within a 15-day sliding window and multiplied by Poisson-distributed multipliers (λ = 4, sensitivity analysis: 11) to generate 1000 estimated infection counts, to which EpiEstim was applied to generate Rt time series. The median Rt percentage change when policies changed was estimated. The time lag correlations were assessed between the 7-day moving average of the relative changes in Google mobility data in the first 90 days, and Rt and estimated infection count, respectively. There were three major pandemic waves in RI in 2020-2021: spring 2020, winter 2020-2021 and fall-winter 2021. The median Rt fluctuated within the range of 0.5-2 from April 2020 to November 2021. Mask mandate (18 April 2020) was associated with a decrease in Rt (-25.99%, 95% CrI: -37.42%, -14.30%). Termination of mask mandates on 6 July 2021 was associated with an increase in Rt (36.74%, 95% CrI: 27.20%, 49.13%). Positive correlations were found between changes in grocery and pharmacy, Rt retail and recreation, transit, and workplace visits, for both Rt and estimated infection count, respectively. Negative correlations were found between changes in residential area visits for both Rt and estimated infection count, respectively. Public health policies enacted in RI were associated with changes in the pandemic trajectory. This ecological study provides further evidence of how non-pharmaceutical interventions and vaccination slowed COVID-19 transmission in RI.

Comparing the Performance of Three Computational Methods for Estimating the Effective Reproduction Number

The effective reproduction number ( R t ) is one of the most important epidemiological parameters, providing suggestions for monitoring the development trend of diseases and also for adjusting the prevention and control policies. However, a few studies have focused on the performance of some common computational methods for Rt. The purpose of this article is to compare the performance of three computational methods for Rt: the time-dependent (TD) method, the new time-varying (NT) method, and the sequential Bayesian (SB) method. Four evaluation methods-accuracy, correlation coefficient, similarity based on trend, and dynamic time warping distance-were used to compare the effectiveness of three computational methods for Rt under different time lags and time windows. The results showed that the NT method was a better choice for real-time monitoring and analysis of the epidemic in the middle and late stages of the infectious disease. The TD method could reflect the change of the number of cases stably and accurately, and was more suitable for monitoring the change of Rt during the whole process of the epidemic outbreak. When the data were relatively stable, the SB method could also provide a reliable estimate for Rt, while the error would increase when the fluctuation in the number of cases increased. The results would provide suggestions for selecting appropriate Rt estimation methods and making policy adjustments more timely and effectively according to the change of Rt.

Wastewater-based reproduction numbers and projections of COVID-19 cases in three areas in Japan, November 2021 to December 2022

BackgroundWastewater surveillance has expanded globally as a means to monitor spread of infectious diseases. An inherent challenge is substantial noise and bias in wastewater data because of the sampling and quantification process, limiting the applicability of wastewater surveillance as a monitoring tool.AimTo present an analytical framework for capturing the growth trend of circulating infections from wastewater data and conducting scenario analyses to guide policy decisions.MethodsWe developed a mathematical model for translating the observed SARS-CoV-2 viral load in wastewater into effective reproduction numbers. We used an extended Kalman filter to infer underlying transmissions by smoothing out observational noise. We also illustrated the impact of different countermeasures such as expanded vaccinations and non-pharmaceutical interventions on the projected number of cases using three study areas in Japan during 2021-22 as an example.ResultsObserved notified cases were matched with the range of cases estimated by our approach with wastewater data only, across different study areas and virus quantification methods, especially when the disease prevalence was high. Estimated reproduction numbers derived from wastewater data were consistent with notification-based reproduction numbers. Our projections showed that a 10-20% increase in vaccination coverage or a 10% reduction in contact rate may suffice to initiate a declining trend in study areas.ConclusionOur study demonstrates how wastewater data can be used to track reproduction numbers and perform scenario modelling to inform policy decisions. The proposed framework complements conventional clinical surveillance, especially when reliable and timely epidemiological data are not available.

Estimates of early outbreak-specific SARS-CoV-2 epidemiological parameters from genomic data

We estimate the basic reproductive number and case counts for 15 distinct Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreaks, distributed across 11 populations (10 countries and one cruise ship), based solely on phylodynamic analyses of genomic data. Our results indicate that, prior to significant public health interventions, the reproductive numbers for 10 (out of 15) of these outbreaks are similar, with median posterior estimates ranging between 1.4 and 2.8. These estimates provide a view which is complementary to that provided by those based on traditional line listing data. The genomic-based view is arguably less susceptible to biases resulting from differences in testing protocols, testing intensity, and import of cases into the community of interest. In the analyses reported here, the genomic data primarily provide information regarding which samples belong to a particular outbreak. We observe that once these outbreaks are identified, the sampling dates carry the majority of the information regarding the reproductive number. Finally, we provide genome-based estimates of the cumulative number of infections for each outbreak. For 7 out of 11 of the populations studied, the number of confirmed cases is much bigger than the cumulative number of infections estimated from the sequence data, a possible explanation being the presence of unsequenced outbreaks in these populations.

Near-term forecasting of Covid-19 cases and hospitalisations in Aotearoa New Zealand

Near-term forecasting of infectious disease incidence and consequent demand for acute healthcare services can support capacity planning and public health responses. Despite well-developed scenario modelling to support the Covid-19 response, Aotearoa New Zealand lacks advanced infectious disease forecasting capacity. We develop a model using Aotearoa New Zealand's unique Covid-19 data streams to predict reported Covid-19 cases, hospital admissions and hospital occupancy. The method combines a semi-mechanistic model for disease transmission to predict cases with Gaussian process regression models to predict the fraction of reported cases that will require hospital treatment. We evaluate forecast performance against out-of-sample data over the period from 2 October 2022 to 23 July 2023. Our results show that forecast performance is reasonably good over a 1-3 week time horizon, although generally deteriorates as the time horizon is lengthened. The model has been operationalised to provide weekly national and regional forecasts in real-time. This study is an important step towards development of more sophisticated situational awareness and infectious disease forecasting tools in Aotearoa New Zealand.

Accuracy of Inferences About the Reproductive Number and Superspreading Potential of SARS-CoV-2 with Incomplete Contact Tracing Data

The basic reproductive number (R0) and superspreading potential ( k ) are key epidemiological parameters that inform our understanding of a disease's transmission. Often these values are estimated using the data obtained from contact tracing studies. Here we performed a simulation study to understand how incomplete data due to preferential contact tracing impacted the accuracy and inferences about the transmission of SARS-CoV-2. Our results indicate that as the number of positive contacts traced decreases, our estimates of R0 tend to decrease and our estimates of ktend to increase. Notably, when there are large amounts of positive contacts missed in the tracing process, we can conclude that there is no indication of superspreading even if we know there is. The results of this study highlight the need for a unified public health response to transmissible diseases.

Substantial decline of phasic dopamine signaling in senescent male rats does not impact dopamine-dependent Pavlovian conditioning

Normal aging is associated with cognitive decline which impacts financial decision making. One of the underlying features of decision making is probability estimation, in which nucleus accumbens dopamine signaling has been implicated. Here we used fast-scan cyclic voltammetry to probe for age differences in dopamine signaling, and pharmacological manipulation to test for age differences in the dopamine dependence of Pavlovian conditioning. We found differences in phasic dopamine signaling to reward delivery, and unconditioned and conditioned stimuli, but no difference in conditioned approach between adult and senescent groups. In addition, we found that dopamine receptor antagonism with flupenthixol (225 μg/kg, i.p.) partially inhibited conditioned approach in the adult group, whereas it completely blocked conditioned approach in the senescent group. Further increase in concentration to 300 μg/kg, i.p. resulted in complete inhibition of conditioned approach behavior in both age groups. Therefore, while phasic dopamine signaling in the nucleus accumbens of senescent animals is greatly diminished in concentration, these animals maintain dopamine dependent Pavlovian conditioning.

Estimating the instantaneous reproduction number (Rt) by using particle filter

Background

Monitoring the transmission of coronavirus disease 2019 (COVID-19) requires accurate estimation of the effective reproduction number (R t ). However, existing methods for calculating R t may yield biased estimates if important real-world factors, such as delays in confirmation, pre-symptomatic transmissions, or imperfect data observation, are not considered.

Method

To include real-world factors, we expanded the susceptible-exposed-infectious-recovered (SEIR) model by incorporating pre-symptomatic (P) and asymptomatic (A) states, creating the SEPIAR model. By utilizing both stochastic and deterministic versions of the model, and incorporating predetermined time series of R t , we generated simulated datasets that simulate real-world challenges in estimating R t . We then compared the performance of our proposed particle filtering method for estimating R t with the existing EpiEstim approach based on renewal equations.

Results

The particle filtering method accurately estimated R t even in the presence of data with delays, pre-symptomatic transmission, and imperfect observation. When evaluating via the root mean square error (RMSE) metric, the performance of the particle filtering method was better in general and was comparable to the EpiEstim approach if perfectly deconvolved infection time series were provided, and substantially better when R t exhibited short-term fluctuations and the data was right truncated.

Conclusions

The SEPIAR model, in conjunction with the particle filtering method, offers a reliable tool for predicting the transmission trend of COVID-19 and assessing the impact of intervention strategies. This approach enables enhanced monitoring of COVID-19 transmission and can inform public health policies aimed at controlling the spread of the disease.

Estimating time-dependent contact: a multi-strain epidemiological model of SARS-CoV-2 on the island of Ireland

Mathematical modelling plays a key role in understanding and predicting the epidemiological dynamics of infectious diseases. We construct a flexible discrete-time model that incorporates multiple viral strains with different transmissibilities to estimate the changing patterns of human contact that generates new infections. Using a Bayesian approach, we fit the model to longitudinal data on hospitalisation with COVID-19 from the Republic of Ireland and Northern Ireland during the first year of the pandemic. We describe the estimated change in human contact in the context of government-mandated non-pharmaceutical interventions in the two jurisdictions on the island of Ireland. We take advantage of the fitted model to conduct counterfactual analyses exploring the impact of lockdown timing and introducing a novel, more transmissible variant. We found substantial differences in human contact between the two jurisdictions during periods of varied restriction easing and December holidays. Our counterfactual analyses reveal that implementing lockdowns earlier would have decreased subsequent hospitalisation substantially in most, but not all cases, and that an introduction of a more transmissible variant - without necessarily being more severe - can cause a large impact on the health care burden.

Estimating effective reproduction number revisited

Accurately estimating the effective reproduction number is crucial for characterizing the transmissibility of infectious diseases to optimize interventions and responses during epidemic outbreaks. In this study, we improve the estimation of the effective reproduction number through two main approaches. First, we derive a discrete model to represent a time series of case counts and propose an estimation method based on this framework. We also conduct numerical experiments to demonstrate the effectiveness of the proposed discretization scheme. By doing so, we enhance the accuracy of approximating the underlying epidemic process compared to previous methods, even when the counting period is similar to the mean generation time of an infectious disease. Second, we employ a negative binomial distribution to model the variability of count data to accommodate overdispersion. Specifically, given that observed incidence counts follow a negative binomial distribution, the posterior distribution of secondary infections is obtained as a Dirichlet multinomial distribution. With this formulation, we establish posterior uncertainty bounds for the effective reproduction number. Finally, we demonstrate the effectiveness of the proposed method using incidence data from the COVID-19 pandemic.

A mathematical model of Marburg virus disease outbreaks and the potential role of vaccination in control

BACKGROUND

Marburg virus disease is an acute haemorrhagic fever caused by Marburg virus. Marburg virus is zoonotic, maintained in nature in Egyptian fruit bats, with occasional spillover infections into humans and nonhuman primates. Although rare, sporadic cases and outbreaks occur in Africa, usually associated with exposure to bats in mines or caves, and sometimes with secondary human-to-human transmission. Outbreaks outside of Africa have also occurred due to importation of infected monkeys. Although all previous Marburg virus disease outbreaks have been brought under control without vaccination, there is nevertheless the potential for large outbreaks when implementation of public health measures is not possible or breaks down. Vaccines could thus be an important additional tool, and development of several candidate vaccines is under way.

METHODS

We developed a branching process model of Marburg virus transmission and investigated the potential effects of several prophylactic and reactive vaccination strategies in settings driven primarily by multiple spillover events as well as human-to-human transmission. Linelist data from the 15 outbreaks up until 2022, as well as an Approximate Bayesian Computational framework, were used to inform the model parameters.

RESULTS

Our results show a low basic reproduction number which varied across outbreaks, from 0.5 [95% CI 0.05-1.8] to 1.2 [95% CI 1.0-1.9] but a high case fatality ratio. Of six vaccination strategies explored, the two prophylactic strategies (mass and targeted vaccination of high-risk groups), as well as a combination of ring and targeted vaccination, were generally most effective, with a probability of potential outbreaks being terminated within 1 year of 0.90 (95% CI 0.90-0.91), 0.89 (95% CI 0.88-0.90), and 0.88 (95% CI 0.87-0.89) compared with 0.68 (0.67-0.69) for no vaccination, especially if the outbreak is driven by zoonotic spillovers and the vaccination campaign initiated as soon as possible after onset of the first case.

CONCLUSIONS

Our study shows that various vaccination strategies can be effective in helping to control outbreaks of MVD, with the best approach varying with the particular epidemiologic circumstances of each outbreak.

Chronicling the 3-year evolution of the COVID-19 pandemic: analysis of disease management, characteristics of major variants, and impacts on pathogenicity

Announced on December 31, 2019, the novel coronavirus arising in Wuhan City, Hubei Province resulted in millions of cases and lives lost. Following intense tracking, coronavirus disease 2019 (COVID-19) was declared a pandemic by the World Health Organization (WHO) in 2020. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified as the cause of COVID-19 and the continuous evolution of the virus has given rise to several variants. In this review, a comprehensive analysis of the response to the pandemic over the first three-year period is provided, focusing on disease management, development of vaccines and therapeutics, and identification of variants. The transmissibility and pathogenicity of SARS-CoV-2 variants including Alpha, Beta, Gamma, Delta, and Omicron are compared. The binding characteristics of the SARS-CoV-2 spike protein to the angiotensin-converting enzyme 2 (ACE2) receptor and reproduction numbers are evaluated. The effects of major variants on disease severity, hospitalisation, and case-fatality rates are outlined. In addition to the spike protein, open reading frames mutations are investigated. We also compare the pathogenicity of SARS-CoV-2 with SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). Overall, this study highlights the strengths and weaknesses of the global response to the pandemic, as well as the importance of prevention and preparedness. Monitoring the evolution of SARS-CoV-2 is critical in identifying and potentially predicting the health outcomes of concerning variants as they emerge. The ultimate goal would be a position in which existing vaccines and therapeutics could be adapted to suit new variants in as close to real-time as possible.

Why are different estimates of the effective reproductive number so different? A case study on COVID-19 in Germany

The effective reproductive number Rt has taken a central role in the scientific, political, and public discussion during the COVID-19 pandemic, with numerous real-time estimates of this quantity routinely published. Disagreement between estimates can be substantial and may lead to confusion among decision-makers and the general public. In this work, we compare different estimates of the national-level effective reproductive number of COVID-19 in Germany in 2020 and 2021. We consider the agreement between estimates from the same method but published at different time points (within-method agreement) as well as retrospective agreement across eight different approaches (between-method agreement). Concerning the former, estimates from some methods are very stable over time and hardly subject to revisions, while others display considerable fluctuations. To evaluate between-method agreement, we reproduce the estimates generated by different groups using a variety of statistical approaches, standardizing analytical choices to assess how they contribute to the observed disagreement. These analytical choices include the data source, data pre-processing, assumed generation time distribution, statistical tuning parameters, and various delay distributions. We find that in practice, these auxiliary choices in the estimation of Rt may affect results at least as strongly as the selection of the statistical approach. They should thus be communicated transparently along with the estimates.

Transmissibility of severe acute respiratory syndrome coronavirus 2 among household contacts of coronavirus disease 2019-positive patients: A community-based study in India

BACKGROUND

This study identified the risk factors for severe acute respiratory syndrome coronavirus 2 infection among household contacts of index patients and determined the incubation period (IP), serial interval, and estimates of secondary infection rate in Kerala, India.

METHODS

We conducted a cohort study in three districts of Kerala among the inhabitants of households of reverse transcriptase polymerase chain reaction-positive coronavirus disease 2019 patients between January and July 2021. About 147 index patients and 362 household contacts were followed up for 28 days to determine reverse transcriptase polymerase chain reaction positivity and the presence of total antibodies against SARS-CoV-2 on days 1, 7, 14, and 28.

RESULTS

The mean IP, serial interval, and generation time were 1.6, 3, and 3.9 days, respectively. The secondary infection rate at 14 days was 43.0%. According to multivariable regression analysis persons who worked outside the home were protected (adjusted odds ratio [aOR], 0.45; 95% confidence interval [CI], 0.24-0.85), whereas those who had kissed the coronavirus disease 2019-positive patients during illness were more than twice at risk of infection (aOR, 2.23; 95% CI, 1.01-5.2) than those who had not kissed the patients. Sharing a toilet with the index patient increased the risk by more than twice (aOR, 2.5; 95% CI, 1.42-4.64) than not sharing a toilet. However, the contacts who reported using masks (aOR, 2.5; 95% CI, 1.4-4.4) were at a higher risk of infection in household settings.

CONCLUSIONS

Household settings have a high secondary infection rate and the changing transmissibility dynamics such as IP, serial interval should be considered in the prevention and control of SARS-CoV-2.

Monitoring the reproductive number of COVID-19 in France: Comparative estimates from three datasets

BACKGROUND

The effective reproduction number (Rt) quantifies the average number of secondary cases caused by one person with an infectious disease. Near-real-time monitoring of Rt during an outbreak is a major indicator used to monitor changes in disease transmission and assess the effectiveness of interventions. The estimation of Rt usually requires the identification of infected cases in the population, which can prove challenging with the available data, especially when asymptomatic people or with mild symptoms are not usually screened. The purpose of this study was to perform sensitivity analysis of Rt estimates for COVID-19 surveillance in France based on three data sources with different sensitivities and specificities for identifying infected cases.

METHODS

We applied a statistical method developed by Cori et al. to estimate Rt using (1) confirmed cases identified from positive virological tests in the population, (2) suspected cases recorded by a national network of emergency departments, and (3) COVID-19 hospital admissions recorded by a national administrative system to manage hospital organization.

RESULTS

Rt estimates in France from May 27, 2020, to August 12, 2022, showed similar temporal trends regardless of the dataset. Estimates based on the daily number of confirmed cases provided an earlier signal than the two other sources, with an average lag of 3 and 6 days for estimates based on emergency department visits and hospital admissions, respectively.

CONCLUSION

The COVID-19 experience confirmed that monitoring temporal changes in Rt was a key indicator to help the public health authorities control the outbreak in real time. However, gaining access to data on all infected people in the population in order to estimate Rt is not straightforward in practice. As this analysis has shown, the opportunity to use more readily available data to estimate Rt trends, provided that it is highly correlated with the spread of infection, provides a practical solution for monitoring the COVID-19 pandemic and indeed any other epidemic.

Phylodynamic approaches to studying avian influenza virus

Avian influenza viruses can cause severe disease in domestic and wild birds and are a pandemic threat. Phylodynamics is the study of how epidemiological, evolutionary, and immunological processes can interact to shape viral phylogenies. This review summarizes how phylodynamic methods have and could contribute to the study of avian influenza viruses. Specifically, we assess how phylodynamics can be used to examine viral spread within and between wild or domestic bird populations at various geographical scales, identify factors associated with virus dispersal, and determine the order and timing of virus lineage movement between geographic regions or poultry production systems. We discuss factors that can complicate the interpretation of phylodynamic results and identify how future methodological developments could contribute to improved control of the virus.

The limits of human mobility traces to predict the spread of COVID-19: A transfer entropy approach

Mobile phone data have been widely used to model the spread of COVID-19; however, quantifying and comparing their predictive value across different settings is challenging. Their quality is affected by various factors and their relationship with epidemiological indicators varies over time. Here, we adopt a model-free approach based on transfer entropy to quantify the relationship between mobile phone-derived mobility metrics and COVID-19 cases and deaths in more than 200 European subnational regions. Using multiple data sources over a one-year period, we found that past knowledge of mobility does not systematically provide statistically significant information on COVID-19 spread. Our approach allows us to determine the best metric for predicting disease incidence in a particular location, at different spatial scales. Additionally, we identify geographic and demographic factors, such as users' coverage and commuting patterns, that explain the (non)observed relationship between mobility and epidemic patterns. Our work provides epidemiologists and public health officials with a general-not limited to COVID-19-framework to evaluate the usefulness of human mobility data in responding to epidemics.

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.

Estimating the time-varying reproduction number for COVID-19 in South Africa during the first four waves using multiple measures of incidence for public and private sectors across four waves

OBJECTIVES

The aim of this study was to quantify transmission trends in South Africa during the first four waves of the COVID-19 pandemic using estimates of the time-varying reproduction number (R) and to compare the robustness of R estimates based on three different data sources, and using data from public and private sector service providers.

METHODS

R was estimated from March 2020 through April 2022, nationally and by province, based on time series of rt-PCR-confirmed cases, hospitalisations, and hospital-associated deaths, using a method that models daily incidence as a weighted sum of past incidence, as implemented in the R package EpiEstim. R was also estimated separately using public and private sector data.

RESULTS

Nationally, the maximum case-based R following the introduction of lockdown measures was 1.55 (CI: 1.43-1.66), 1.56 (CI: 1.47-1.64), 1.46 (CI: 1.38-1.53) and 3.33 (CI: 2.84-3.97) during the first (Wuhan-Hu), second (Beta), third (Delta), and fourth (Omicron) waves, respectively. Estimates based on the three data sources (cases, hospitalisations, deaths) were generally similar during the first three waves, but higher during the fourth wave for case-based estimates. Public and private sector R estimates were generally similar except during the initial lockdowns and in case-based estimates during the fourth wave.

CONCLUSION

Agreement between R estimates using different data sources during the first three waves suggests that data from any of these sources could be used in the early stages of a future pandemic. The high R estimates for Omicron relative to earlier waves are interesting given a high level of exposure pre-Omicron. The agreement between public and private sector R estimates highlights that clients of the public and private sectors did not experience two separate epidemics, except perhaps to a limited extent during the strictest lockdowns in the first wave.

The effective reproductive number: Modeling and prediction with application to the multi-wave Covid-19 pandemic

Classical compartmental models of infectious disease assume that spread occurs through a homogeneous population. This produces poor fits to real data, because individuals vary in their number of epidemiologically-relevant contacts, and hence in their ability to transmit disease. In particular, network theory suggests that super-spreading events tend to happen more often at the beginning of an epidemic, which is inconsistent with the homogeneity assumption. In this paper we argue that a flexible decay shape for the effective reproductive number (Rt) indexed by the susceptible fraction (St) is a theory-informed modeling choice, which better captures the progression of disease incidence over human populations. This, in turn, produces better retrospective fits, as well as more accurate prospective predictions of observed epidemic curves. We extend this framework to fit multi-wave epidemics, and to accommodate public health restrictions on mobility. We demonstrate the performance of this model by doing a prediction study over two years of the SARS-CoV2 pandemic.

Evaluating the use of social contact data to produce age-specific short-term forecasts of SARS-CoV-2 incidence in England

Mathematical and statistical models can be used to make predictions of how epidemics may progress in the near future and form a central part of outbreak mitigation and control. Renewal equation based models allow inference of epidemiological parameters from historical data and forecast future epidemic dynamics without requiring complex mechanistic assumptions. However, these models typically ignore interaction between age groups, partly due to challenges in parameterising a time varying interaction matrix. Social contact data collected regularly during the COVID-19 epidemic provide a means to inform interaction between age groups in real-time. We developed an age-specific forecasting framework and applied it to two age-stratified time-series: incidence of SARS-CoV-2 infection, estimated from a national infection and antibody prevalence survey; and, reported cases according to the UK national COVID-19 dashboard. Jointly fitting our model to social contact data from the CoMix study, we inferred a time-varying next generation matrix which we used to project infections and cases in the four weeks following each of 29 forecast dates between October 2020 and November 2021. We evaluated the forecasts using proper scoring rules and compared performance with three other models with alternative data and specifications alongside two naive baseline models. Overall, incorporating age interaction improved forecasts of infections and the CoMix-data-informed model was the best performing model at time horizons between two and four weeks. However, this was not true when forecasting cases. We found that age group interaction was most important for predicting cases in children and older adults. The contact-data-informed models performed best during the winter months of 2020-2021, but performed comparatively poorly in other periods. We highlight challenges regarding the incorporation of contact data in forecasting and offer proposals as to how to extend and adapt our approach, which may lead to more successful forecasts in future.

estimateR: an R package to estimate and monitor the effective reproductive number

BACKGROUND

Accurate estimation of the effective reproductive number ([Formula: see text]) of epidemic outbreaks is of central relevance to public health policy and decision making. We present estimateR, an R package for the estimation of the reproductive number through time from delayed observations of infection events. Such delayed observations include confirmed cases, hospitalizations or deaths. The package implements the methodology of Huisman et al. but modularizes the [Formula: see text] estimation procedure to allow easy implementation of new alternatives to the currently available methods. Users can tailor their analyses according to their particular use case by choosing among implemented options.

RESULTS

The estimateR R package allows users to estimate the effective reproductive number of an epidemic outbreak based on observed cases, hospitalization, death or any other type of event documenting past infections, in a fast and timely fashion. We validated the implementation with a simulation study: estimateR yielded estimates comparable to alternative publicly available methods while being around two orders of magnitude faster. We then applied estimateR to empirical case-confirmation incidence data for COVID-19 in nine countries and for dengue fever in Brazil; in parallel, estimateR is already being applied (i) to SARS-CoV-2 measurements in wastewater data and (ii) to study influenza transmission based on wastewater and clinical data in other studies. In summary, this R package provides a fast and flexible implementation to estimate the effective reproductive number for various diseases and datasets.

CONCLUSIONS

The estimateR R package is a modular and extendable tool designed for outbreak surveillance and retrospective outbreak investigation. It extends the method developed for COVID-19 by Huisman et al. and makes it available for a variety of pathogens, outbreak scenarios, and observation types. Estimates obtained with estimateR can be interpreted directly or used to inform more complex epidemic models (e.g. for forecasting) on the value of [Formula: see text].

Estimating the epidemic reproduction number from temporally aggregated incidence data: A statistical modelling approach and software tool

The time-varying reproduction number (Rt) is an important measure of epidemic transmissibility that directly informs policy decisions and the optimisation of control measures. EpiEstim is a widely used opensource software tool that uses case incidence and the serial interval (SI, time between symptoms in a case and their infector) to estimate Rt in real-time. The incidence and the SI distribution must be provided at the same temporal resolution, which can limit the applicability of EpiEstim and other similar methods, e.g. for contexts where the time window of incidence reporting is longer than the mean SI. In the EpiEstim R package, we implement an expectation-maximisation algorithm to reconstruct daily incidence from temporally aggregated data, from which Rt can then be estimated. We assess the validity of our method using an extensive simulation study and apply it to COVID-19 and influenza data. For all datasets, the influence of intra-weekly variability in reported data was mitigated by using aggregated weekly data. Rt estimated on weekly sliding windows using incidence reconstructed from weekly data was strongly correlated with estimates from the original daily data. The simulation study revealed that Rt was well estimated in all scenarios and regardless of the temporal aggregation of the data. In the presence of weekend effects, Rt estimates from reconstructed data were more successful at recovering the true value of Rt than those obtained from reported daily data. These results show that this novel method allows Rt to be successfully recovered from aggregated data using a simple approach with very few data requirements. Additionally, by removing administrative noise when daily incidence data are reconstructed, the accuracy of Rt estimates can be improved.

Anchoring the mean generation time in the SEIR to mitigate biases in ℜ0 estimates due to uncertainty in the distribution of the epidemiological delays

The basic reproduction number, ℜ 0 , is of paramount importance in the study of infectious disease dynamics. Primarily, ℜ 0 serves as an indicator of the transmission potential of an emerging infectious disease and the effort required to control the invading pathogen. However, its estimates from compartmental models are strongly conditioned by assumptions in the model structure, such as the distributions of the latent and infectious periods (epidemiological delays). To further complicate matters, models with dissimilar delay structures produce equivalent incidence dynamics. Following a simulation study, we reveal that the nature of such equivalency stems from a linear relationship between ℜ 0 and the mean generation time, along with adjustments to other parameters in the model. Leveraging this knowledge, we propose and successfully test an alternative parametrization of the SEIR model that produces accurate ℜ 0 estimates regardless of the distribution of the epidemiological delays, at the expense of biases in other quantities deemed of lesser importance. We further explore this approach's robustness by testing various transmissibility levels, generation times and data fidelity (overdispersion). Finally, we apply the proposed approach to data from the 1918 influenza pandemic. We anticipate that this work will mitigate biases in estimating ℜ 0 .

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.

Dynamics-informed deconvolutional neural networks for super-resolution identification of regime changes in epidemiological time series

The ability to infer the timing and amplitude of perturbations in epidemiological systems from their stochastically spread low-resolution outcomes is crucial for multiple applications. However, the general problem of connecting epidemiological curves with the underlying incidence lacks the highly effective methodology present in other inverse problems, such as super-resolution and dehazing from computer vision. Here, we develop an unsupervised physics-informed convolutional neural network approach in reverse to connect death records with incidence that allows the identification of regime changes at single-day resolution. Applied to COVID-19 data with proper regularization and model-selection criteria, the approach can identify the implementation and removal of lockdowns and other nonpharmaceutical interventions (NPIs) with 0.93-day accuracy over the time span of a year.

Modeling vaccination coverage during the 2022 central Ohio measles outbreak: a cross-sectional study

Background

Of the eight large (>50 cases) US postelimination outbreaks, the first and last occurred in Ohio. Ohio's vaccination registry is incomplete. Community-level immunity gaps threaten more than two decades of measles elimination in the US. We developed a statistical model, VaxEstim, to rapidly estimate the early-phase vaccination coverage and immunity gap in the exposed population during the 2022 Central Ohio outbreak.

Methods

We used reconstructed daily incidence (from publicly available data) and assumptions about the distribution of the serial interval, or the time between symptom onset in successive measles cases, to estimate the effective reproduction number (i.e., the average number of secondary infections caused by an infected individual in a partially immune population). We estimated early-phase measles vaccination coverage by comparing the effective reproduction number to the basic reproduction number (i.e., the average number of secondary infections caused by an infected individual in a fully susceptible population) while accounting for vaccine effectiveness. Finally, we estimated the early-phase immunity gap as the difference between the estimated critical vaccination threshold and vaccination coverage.

Findings

VaxEstim estimated the early-phase vaccination coverage as 53% (95% credible interval, 21%-77%), the critical vaccination threshold as 93%, and the immunity gap as 42% (95% credible interval, 18%-74%).

Interpretation

This study estimates a significant immunity gap in the exposed population during the early phase of the 2022 Central Ohio measles outbreak, suggesting a robust public health response is needed to identify the susceptible community and develop community-specific strategies to close the immunity gap.

Funding

This work was supported in part by the National Institute of General Medical Sciences, National Institutes of Health; the UK Medical Research Council (MRC); the Foreign, Commonwealth and Development Office; the National Institute for Health Research (NIHR) Health Protection Research Unit in Modelling Methodology; Imperial College London, and the London School of Hygiene & Tropical Medicine, Community Jameel; the EDCTP2 programme, supported by the EU; and the Sergei Brin Foundation.

Understanding the role of mask-wearing during COVID-19 on the island of Ireland

Non-pharmaceutical interventions have played a key role in managing the COVID-19 pandemic, but it is challenging to estimate their impacts on disease spread and outcomes. On the island of Ireland, population mobility restrictions were imposed during the first wave, but mask-wearing was not mandated until about six months into the pandemic. We use data on mask-wearing, mobility, and season, over the first year of the pandemic to predict independently the weekly infectious contact estimated by an epidemiological model. Using our models, we make counterfactual predictions of infectious contact, and ensuing hospitalizations, under a hypothetical intervention where 90% of the population wore masks from the beginning of community spread until the dates of the mask mandates. Over periods including the first wave of the pandemic, there were 1601 hospitalizations with COVID-19 in Northern Ireland and 1521 in the Republic of Ireland. Under the counterfactual mask-wearing scenario, we estimate 512 (95% CI 400, 730) and 344 (95% CI 266, 526) hospitalizations in the respective jurisdictions during the same periods. This could be partly due to other factors that were also changing over time.

Trajectories of COVID-19: A longitudinal analysis of many nations and subnational regions

The COVID-19 pandemic is the first to be rapidly and sequentially measured by nation-wide PCR community testing for the presence of the viral RNA at a global scale. We take advantage of the novel "natural experiment" where diverse nations and major subnational regions implemented various policies including social distancing and vaccination at different times with different levels of stringency and adherence. Initially, case numbers expand exponentially with doubling times of ~1-2 weeks. In the nations where interventions were not implemented or perhaps lees effectual, case numbers increased exponentially but then stabilized around 102-to-103 new infections (per km2 built-up area per day). Dynamics under effective interventions were perturbed and infections decayed to low levels. They rebounded concomitantly with the lifting of social distancing policies or pharmaceutical efficacy decline, converging on a stable equilibrium setpoint. Here we deploy a mathematical model which captures this V-shape behavior, incorporating a direct measure of intervention efficacy. Importantly, it allows the derivation of a maximal estimate for the basic reproductive number Ro (mean 1.6-1.8). We were able to test this approach by comparing the approximated "herd immunity" to the vaccination coverage observed that corresponded to rapid declines in community infections during 2021. The estimates reported here agree with the observed phenomena. Moreover, the decay (0.4-0.5) and rebound rates (0.2-0.3) were similar throughout the pandemic and among all the nations and regions studied. Finally, a longitudinal analysis comparing multiple national and regional results provides insights on the underlying epidemiology of SARS-CoV-2 and intervention efficacy, as well as evidence for the existence of an endemic steady state of COVID-19.

Multivariate epidemic count time series model

An infectious disease spreads not only over a single population or community but also across multiple and heterogeneous communities. Moreover, its transmissibility varies over time because of various factors such as seasonality and epidemic control, which results in strongly nonstationary behavior. In conventional methods for assessing transmissibility trends or changes, univariate time-varying reproduction numbers are calculated without taking into account transmission across multiple communities. In this paper, we propose a multivariate-count time series model for epidemics. We also propose a statistical method for estimating the transmission of infections across multiple communities and the time-varying reproduction numbers of each community simultaneously from a multivariate time series of case counts. We apply our method to incidence data for the novel coronavirus disease 2019 (COVID-19) pandemic to reveal the spatiotemporal heterogeneity of the epidemic process.

Key epidemiological indicators and spatial autocorrelation patterns across five waves of COVID-19 in Catalonia

This research studies the evolution of COVID-19 crude incident rates, effective reproduction number R(t) and their relationship with incidence spatial autocorrelation patterns in the 19 months following the disease outbreak in Catalonia (Spain). A cross-sectional ecological panel design based on n = 371 health-care geographical units is used. Five general outbreaks are described, systematically preceded by generalized values of R(t) > 1 in the two previous weeks. No clear regularities concerning possible initial focus appear when comparing waves. As for autocorrelation, we identify a wave's baseline pattern in which global Moran's I increases rapidly in the first weeks of the outbreak to descend later. However, some waves significantly depart from the baseline. In the simulations, both baseline pattern and departures can be reproduced when measures aimed at reducing mobility and virus transmissibility are introduced. Spatial autocorrelation is inherently contingent on the outbreak phase and is also substantially modified by external interventions affecting human behavior.

No evidence of sustained nonzoonotic Plasmodium knowlesi transmission in Malaysia from modelling malaria case data

Reported incidence of the zoonotic malaria Plasmodium knowlesi has markedly increased across Southeast Asia and threatens malaria elimination. Nonzoonotic transmission of P. knowlesi has been experimentally demonstrated, but it remains unknown whether nonzoonotic transmission is contributing to increases in P. knowlesi cases. Here, we adapt model-based inference methods to estimate RC, individual case reproductive numbers, for P. knowlesi, P. falciparum and P. vivax human cases in Malaysia from 2012-2020 (n = 32,635). Best fitting models for P. knowlesi showed subcritical transmission (RC < 1) consistent with a large reservoir of unobserved infection sources, indicating P. knowlesi remains a primarily zoonotic infection. In contrast, sustained transmission (RC > 1) was estimated historically for P. falciparum and P. vivax, with declines in RC estimates observed over time consistent with local elimination. Together, this suggests sustained nonzoonotic P. knowlesi transmission is highly unlikely and that new approaches are urgently needed to control spillover risks.

Identification of growing tuberculosis incidence clusters in a region with a decrease in tuberculosis prevalence in Moscow (2000-2019)

Background

The control of tuberculosis (TB) may benefit from a prospective identification of areas where the incidence may increase in addition to the traditionally identified foci of high incidence. We aimed to identify residential areas with growing tuberculosis incidence rates and assess their significance and stability.

Methods

We analysed the changes in TB incidence rates using case data georeferenced with spatial granularity to apartment buildings in the territory of Moscow from 2000 to 2019. We identified sparsely distributed areas with significant increases in the incidence rate inside residential areas. We tested the stability of found growth areas to case underreporting via stochastic modelling.

Results

For 21 350 cases with smear- or culture-positive pulmonary TB among residents from 2000 to 2019, we identified 52 small-scale clusters of growing incidence rate responsible for 1% of all registered cases. We tested clusters of disease growth for underreporting and found them to be relatively unstable to resampling with case drop-out, but their spatial displacement was small. Territories with a stable increase in TB incidence rate were identified and compared to the rest of the city, which is characterised by a significant decrease in incidence.

Conclusions

Identified areas with a tendency for an increase in the TB incidence rate may be important targets for disease control services.

Design of effective outpatient sentinel surveillance for COVID-19 decision-making: a modeling study

BACKGROUND

Decision-makers impose COVID-19 mitigations based on public health indicators such as reported cases, which are sensitive to fluctuations in supply and demand for diagnostic testing, and hospital admissions, which lag infections by up to two weeks. Imposing mitigations too early has unnecessary economic costs while imposing too late leads to uncontrolled epidemics with unnecessary cases and deaths. Sentinel surveillance of recently-symptomatic individuals in outpatient testing sites may overcome biases and lags in conventional indicators, but the minimal outpatient sentinel surveillance system needed for reliable trend estimation remains unknown.

METHODS

We used a stochastic, compartmental transmission model to evaluate the performance of various surveillance indicators at reliably triggering an alarm in response to, but not before, a step increase in transmission of SARS-CoV-2. The surveillance indicators included hospital admissions, hospital occupancy, and sentinel cases with varying levels of sampling effort capturing 5, 10, 20, 50, or 100% of incident mild cases. We tested 3 levels of transmission increase, 3 population sizes, and conditions of either simultaneous transmission increase or lagged increase in the older population. We compared the indicators' performance at triggering alarm soon after, but not prior, to the transmission increase.

RESULTS

Compared to surveillance based on hospital admissions, outpatient sentinel surveillance that captured at least 20% of incident mild cases could trigger an alarm 2 to 5 days earlier for a mild increase in transmission and 6 days earlier for a moderate or strong increase. Sentinel surveillance triggered fewer false alarms and averted more deaths per day spent in mitigation. When transmission increase in older populations lagged the increase in younger populations by 14 days, sentinel surveillance extended its lead time over hospital admissions by an additional 2 days.

CONCLUSIONS

Sentinel surveillance of mild symptomatic cases can provide more timely and reliable information on changes in transmission to inform decision-makers in an epidemic like COVID-19.

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.

A new method for the joint estimation of instantaneous reproductive number and serial interval during epidemics

Although some methods for estimating the instantaneous reproductive number during epidemics have been developed, the existing frameworks usually require information on the distribution of the serial interval and/or additional contact tracing data. However, in the case of outbreaks of emerging infectious diseases with an unknown natural history or undetermined characteristics, the serial interval and/or contact tracing data are often not available, resulting in inaccurate estimates for this quantity. In the present study, a new framework was specifically designed for joint estimates of the instantaneous reproductive number and serial interval. Concretely, a likelihood function for the two quantities was first introduced. Then, the instantaneous reproductive number and the serial interval were modeled parametrically as a function of time using the interpolation method and a known traditional distribution, respectively. Using the Bayesian information criterion and the Markov Chain Monte Carlo method, we ultimately obtained their estimates and distribution. The simulation study revealed that our estimates of the two quantities were consistent with the ground truth. Seven data sets of historical epidemics were considered and further verified the robust performance of our method. Therefore, to some extent, even if we know only the daily incidence, our method can accurately estimate the instantaneous reproductive number and serial interval to provide crucial information for policymakers to design appropriate prevention and control interventions during epidemics.