The importance of the generation interval in investigating dynamics and control of new SARS-CoV-2 variants

Early mutational signatures and transmissibility of SARS-CoV-2 Gamma and Lambda variants in Chile

Genomic surveillance (GS) programmes were crucial in identifying and quantifying the mutating patterns of SARS-CoV-2 during the COVID-19 pandemic. In this work, we develop a Bayesian framework to quantify the relative transmissibility of different variants tailored for regions with limited GS. We use it to study the relative transmissibility of SARS-CoV-2 variants in Chile. Among the 3443 SARS-CoV-2 genomes collected between January and June 2021, where sampling was designed to be representative, the Gamma (P.1), Lambda (C.37), Alpha (B.1.1.7), B.1.1.348, and B.1.1 lineages were predominant. We found that Lambda and Gamma variants' reproduction numbers were 5% (95% CI: [1%, 14%]) and 16% (95% CI: [11%, 21%]) larger than Alpha's, respectively. Besides, we observed a systematic mutation enrichment in the Spike gene for all circulating variants, which strongly correlated with variants' transmissibility during the studied period (r = 0.93, p-value = 0.025). We also characterised the mutational signatures of local samples and their evolution over time and with the progress of vaccination, comparing them with those of samples collected in other regions worldwide. Altogether, our work provides a reliable method for quantifying variant transmissibility under subsampling and emphasises the importance of continuous genomic surveillance.

Concepts and methods for predicting viral evolution

The seasonal human influenza virus undergoes rapid evolution, leading to significant changes in circulating viral strains from year to year. These changes are typically driven by adaptive mutations, particularly in the antigenic epitopes, the regions of the viral surface protein haemagglutinin targeted by human antibodies. Here we describe a consistent set of methods for data-driven predictive analysis of viral evolution. Our pipeline integrates four types of data: (1) sequence data of viral isolates collected on a worldwide scale, (2) epidemiological data on incidences, (3) antigenic characterization of circulating viruses, and (4) intrinsic viral phenotypes. From the combined analysis of these data, we obtain estimates of relative fitness for circulating strains and predictions of clade frequencies for periods of up to one year. Furthermore, we obtain comparative estimates of protection against future viral populations for candidate vaccine strains, providing a basis for pre-emptive vaccine strain selection. Continuously updated predictions obtained from the prediction pipeline for influenza and SARS-CoV-2 are available on the website previr.app .

Concepts and methods for predicting viral evolution

The seasonal human influenza virus undergoes rapid evolution, leading to significant changes in circulating viral strains from year to year. These changes are typically driven by adaptive mutations, particularly in the antigenic epitopes, the regions of the viral surface protein haemagglutinin targeted by human antibodies. Here we describe a consistent set of methods for data-driven predictive analysis of viral evolution. Our pipeline integrates four types of data: (1) sequence data of viral isolates collected on a worldwide scale, (2) epidemiological data on incidences, (3) antigenic characterization of circulating viruses, and (4) intrinsic viral phenotypes. From the combined analysis of these data, we obtain estimates of relative fitness for circulating strains and predictions of clade frequencies for periods of up to one year. Furthermore, we obtain comparative estimates of protection against future viral populations for candidate vaccine strains, providing a basis for pre-emptive vaccine strain selection. Continuously updated predictions obtained from the prediction pipeline for influenza and SARS-CoV-2 are available on the website previr.app.

The impact of quality-adjusted life years on evaluating COVID-19 mitigation strategies: lessons from age-specific vaccination roll-out and variants of concern in Belgium (2020-2022)

BACKGROUND

When formulating and evaluating COVID-19 vaccination strategies, an emphasis has been placed on preventing severe disease that overburdens healthcare systems and leads to mortality. However, more conventional outcomes such as quality-adjusted life years (QALYs) and inequality indicators are warranted as additional information for policymakers.

METHODS

We adopted a mathematical transmission model to describe the infectious disease dynamics of SARS-COV-2, including disease mortality and morbidity, and to evaluate (non)pharmaceutical interventions. Therefore, we considered temporal immunity levels, together with the distinct transmissibility of variants of concern (VOCs) and their corresponding vaccine effectiveness. We included both general and age-specific characteristics related to SARS-CoV-2 vaccination. Our scenario study is informed by data from Belgium, focusing on the period from August 2021 until February 2022, when vaccination for children aged 5-11 years was initially not yet licensed and first booster doses were administered to adults. More specifically, we investigated the potential impact of an earlier vaccination programme for children and increased or reduced historical adult booster dose uptake.

RESULTS

Through simulations, we demonstrate that increasing vaccine uptake in children aged 5-11 years in August-September 2021 could have led to reduced disease incidence and ICU occupancy, which was an essential indicator for implementing non-pharmaceutical interventions and maintaining healthcare system functionality. However, an enhanced booster dose regimen for adults from November 2021 onward could have resulted in more substantial cumulative QALY gains, particularly through the prevention of elevated levels of infection and disease incidence associated with the emergence of Omicron VOC. In both scenarios, the need for non-pharmaceutical interventions could have decreased, potentially boosting economic activity and mental well-being.

CONCLUSIONS

When calculating the impact of measures to mitigate disease spread in terms of life years lost due to COVID-19 mortality, we highlight the impact of COVID-19 on the health-related quality of life of survivors. Our study underscores that disease-related morbidity could constitute a significant part of the overall health burden. Our quantitative findings depend on the specific setup of the interventions under review, which is open to debate or should be contextualised within future situations.

Differential contagiousness of respiratory disease across the United States

The initial contagiousness of a communicable disease within a given population is quantified by the basic reproduction number, R0. This number depends on both pathogen and population properties. On the basis of compartmental models that reproduce Coronavirus Disease 2019 (COVID-19) surveillance data, we used Bayesian inference and the next-generation matrix approach to estimate region-specific R0 values for 280 of 384 metropolitan statistical areas (MSAs) in the United States (US), which account for 95% of the US population living in urban areas and 82% of the total population. We focused on MSA populations after finding that these populations were more uniformly impacted by COVID-19 than state populations. Our maximum a posteriori (MAP) estimates for R0 range from 1.9 to 7.7 and quantify the relative susceptibilities of regional populations to spread of respiratory diseases. ONE-SENTENCE SUMMARY: Initial contagiousness of Coronavirus Disease 2019 varied over a 4-fold range across urban areas of the United States.

Population immunity predicts evolutionary trajectories of SARS-CoV-2

The large-scale evolution of the SARS-CoV-2 virus has been marked by rapid turnover of genetic clades. New variants show intrinsic changes, notably increased transmissibility, and antigenic changes that reduce cross-immunity induced by previous infections or vaccinations. How this functional variation shapes global evolution has remained unclear. Here, we establish a predictive fitness model for SARS-CoV-2 that integrates antigenic and intrinsic selection. The model is informed by tracking of time-resolved sequence data, epidemiological records, and cross-neutralization data of viral variants. Our inference shows that immune pressure, including contributions of vaccinations and previous infections, has become the dominant force driving the recent evolution of SARS-CoV-2. The fitness model can serve continued surveillance in two ways. First, it successfully predicts the short-term evolution of circulating strains and flags emerging variants likely to displace the previously predominant variant. Second, it predicts likely antigenic profiles of successful escape variants prior to their emergence.

Fitness, growth and transmissibility of SARS-CoV-2 genetic variants

The massive scale of the global SARS-CoV-2 sequencing effort created new opportunities and challenges for understanding SARS-CoV-2 evolution. Rapid detection and assessment of new variants has become one of the principal objectives of genomic surveillance of SARS-CoV-2. Because of the pace and scale of sequencing, new strategies have been developed for characterizing fitness and transmissibility of emerging variants. In this Review, I discuss a wide range of approaches that have been rapidly developed in response to the public health threat posed by emerging variants, ranging from new applications of classic population genetics models to contemporary synthesis of epidemiological models and phylodynamic analysis. Many of these approaches can be adapted to other pathogens and will have increasing relevance as large-scale pathogen sequencing becomes a regular feature of many public health systems.

Improved time-varying reproduction numbers using the generation interval for COVID-19

Estimating key epidemiological parameters, such as incubation period, serial interval (SI), generation interval (GI) and latent period, is essential to quantify the transmissibility and effects of various interventions of COVID-19. These key parameters play a critical role in quantifying the basic reproduction number. With the hard work of epidemiological investigators in South Korea, estimating these key parameters has become possible based on infector-infectee surveillance data of COVID-19 between February 2020 and April 2021. Herein, the mean incubation period was estimated to be 4.9 days (95% CI: 4.2, 5.7) and the mean generation interval was estimated to be 4.3 days (95% CI: 4.2, 4.4). The mean serial interval was estimated to be 4.3, with a standard deviation of 4.2. It is also revealed that the proportion of presymptomatic transmission was ~57%, which indicates the potential risk of transmission before the disease onset. We compared the time-varying reproduction number based on GI and SI and found that the time-varying reproduction number based on GI may result in a larger estimation of Rt, which refers to the COVID-19 transmission potential around the rapid increase of cases. This highlights the importance of considering presymptomatic transmission and generation intervals when estimating the time-varying reproduction number.

Community incidence patterns drive the risk of SARS-CoV-2 outbreaks and alter intervention impacts in a high-risk institutional setting

Optimization of control measures for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in high-risk institutional settings (e.g., prisons, nursing homes, or military bases) depends on how transmission dynamics in the broader community influence outbreak risk locally. We calibrated an individual-based transmission model of a military training camp to the number of RT-PCR positive trainees throughout 2020 and 2021. The predicted number of infected new arrivals closely followed adjusted national incidence and increased early outbreak risk after accounting for vaccination coverage, masking compliance, and virus variants. Outbreak size was strongly correlated with the predicted number of off-base infections among staff during training camp. In addition, off-base infections reduced the impact of arrival screening and masking, while the number of infectious trainees upon arrival reduced the impact of vaccination and staff testing. Our results highlight the importance of outside incidence patterns for modulating risk and the optimal mixture of control measures in institutional settings.

Epidemic modeling with heterogeneity and social diffusion

We propose and analyze a family of epidemiological models that extend the classic Susceptible-Infectious-Recovered/Removed (SIR)-like framework to account for dynamic heterogeneity in infection risk. The family of models takes the form of a system of reaction-diffusion equations given populations structured by heterogeneous susceptibility to infection. These models describe the evolution of population-level macroscopic quantities S, I, R as in the classical case coupled with a microscopic variable f, giving the distribution of individual behavior in terms of exposure to contagion in the population of susceptibles. The reaction terms represent the impact of sculpting the distribution of susceptibles by the infection process. The diffusion and drift terms that appear in a Fokker-Planck type equation represent the impact of behavior change both during and in the absence of an epidemic. We first study the mathematical foundations of this system of reaction-diffusion equations and prove a number of its properties. In particular, we show that the system will converge back to the unique equilibrium distribution after an epidemic outbreak. We then derive a simpler system by seeking self-similar solutions to the reaction-diffusion equations in the case of Gaussian profiles. Notably, these self-similar solutions lead to a system of ordinary differential equations including classic SIR-like compartments and a new feature: the average risk level in the remaining susceptible population. We show that the simplified system exhibits a rich dynamical structure during epidemics, including plateaus, shoulders, rebounds and oscillations. Finally, we offer perspectives and caveats on ways that this family of models can help interpret the non-canonical dynamics of emerging infectious diseases, including COVID-19.

Correlation between times to SARS-CoV-2 symptom onset and secondary transmission undermines epidemic control efforts

Severe acute respiratory coronavirus 2 (SARS-CoV-2) infections have been associated with substantial presymptomatic transmission, which occurs when the generation interval-the time between infection of an individual with a pathogen and transmission of the pathogen to another individual-is shorter than the incubation period-the time between infection and symptom onset. We collected a dataset of 257 SARS-CoV-2 transmission pairs in Japan during 2020 and jointly estimated the mean incubation period of infectors (4.8 days, 95 % CrI: 4.4-5.1 days), mean generation interval to when they infect others (4.3 days, 95 % credible interval [CrI]: 4.0-4.7 days), and the correlation (Kendall's tau: 0.5, 95 % CrI: 0.4-0.6) between these two epidemiological parameters. Our finding of a positive correlation and mean generation interval shorter than the mean infector incubation period indicates ample infectiousness before symptom onset and suggests that reliance on isolation of symptomatic COVID-19 cases as a focal point of control efforts is insufficient to address the challenges posed by SARS-CoV-2 transmission dynamics.