Rapid incidence estimation from SARS-CoV-2 genomes reveals decreased case detection in Europe during summer 2020

Genomic Epidemiology for Estimating Pathogen Burden in a Population

The role of genomics in outbreak response and pathogen surveillance has expanded and ushered in the age of pathogen intelligence. Genomic surveillance enables detection and monitoring of novel pathogens; case clusters; and markers of virulence, antimicrobial resistance, and immune escape. We can leverage pathogen genomic diversity to estimate total pathogen burden in populations and environments, which was previously challenging because of unreliable data. Pathogen genomics might allow pathogen burdens to be estimated by sequencing even a small percentage of cases. Deeper genomic epidemiology analyses require multidisciplinary collaboration to ensure accurate and actionable real-time pathogen intelligence.

A novel approach to forecasting reproduction numbers of spatiotemporal stochastic epidemic spread using a PDE-based model and real-time infection data

The COVID-19 pandemic highlighted the need for improved epidemic spread forecasting, a critical precursor for developing optimal control measures for spread mitigation. Well-recognized shortcomings in computing basic and effective reproduction numbers ( R 0 , R e )-fundamental metrics for forecasting-underscore the need for new methods for estimating them from available data. We present a novel computational framework for estimating reproduction numbers from empirical spread data. The framework is derived from a mechanistic, spatiotemporal, Partial Differential Equation (PDE) model of epidemic spread utilizing mathematical results from PDE epidemic models. Forecasts of spatiotemporal effective reproduction number R e using the framework are found to be in excellent agreement with COVID-19 spread trends for Hamilton County, Ohio, USA, for three distinct periods. Furthermore, the forecasts are shown to align with corresponding reproduction numbers computed independently using the Wallinga-Teunis and Cori retrospective methods used in epidemiology. In summary, the results establish the validity of the framework and indicate applicability to future epidemics-especially for regions such as counties and for timeframes extending in weeks-even during dynamic phases when obtainable real-time infection spread data will likely be sparse.

SARS-CoV-2 evolution on a dynamic immune landscape

Since the onset of the pandemic, many SARS-CoV-2 variants have emerged, exhibiting substantial evolution in the virus' spike protein1, the main target of neutralizing antibodies2. A plausible hypothesis proposes that the virus evolves to evade antibody-mediated neutralization (vaccine- or infection-induced) to maximize its ability to infect an immunologically experienced population1,3. Because viral infection induces neutralizing antibodies, viral evolution may thus navigate on a dynamic immune landscape that is shaped by local infection history. Here we developed a comprehensive mechanistic model, incorporating deep mutational scanning data4,5, antibody pharmacokinetics and regional genomic surveillance data, to predict the variant-specific relative number of susceptible individuals over time. We show that this quantity precisely matched historical variant dynamics, predicted future variant dynamics and explained global differences in variant dynamics. Our work strongly suggests that the ongoing pandemic continues to shape variant-specific population immunity, which determines a variant's ability to transmit, thus defining variant fitness. The model can be applied to any region by utilizing local genomic surveillance data, allows contextualizing risk assessment of variants and provides information for vaccine design.

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 hemagglutinin 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 1 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 at https://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.

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.

Evolution of Sequence and Structure of SARS-CoV-2 Spike Protein: A Dynamic Perspective

Novel coronavirus (SARS-CoV-2) enters its host cell through a surface spike protein. The viral spike protein has undergone several modifications/mutations at the genomic level, through which it modulated its structure-function and passed through several variants of concern. Recent advances in high-resolution structure determination and multiscale imaging techniques, cost-effective next-generation sequencing, and development of new computational methods (including information theory, statistical methods, machine learning, and many other artificial intelligence-based techniques) have hugely contributed to the characterization of sequence, structure, function of spike proteins, and its different variants to understand viral pathogenesis, evolutions, and transmission. Laying on the foundation of the sequence-structure-function paradigm, this review summarizes not only the important findings on structure/function but also the structural dynamics of different spike components, highlighting the effects of mutations on them. As dynamic fluctuations of three-dimensional spike structure often provide important clues for functional modulation, quantifying time-dependent fluctuations of mutational events over spike structure and its genetic/amino acidic sequence helps identify alarming functional transitions having implications for enhanced fusogenicity and pathogenicity of the virus. Although these dynamic events are more difficult to capture than quantifying a static, average property, this review encompasses those challenging aspects of characterizing the evolutionary dynamics of spike sequence and structure and their implications for functions.

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

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

Secondary data for global health digitalisation

Substantial opportunities for global health intelligence and research arise from the combined and optimised use of secondary data within data ecosystems. Secondary data are information being used for purposes other than those intended when they were collected. These data can be gathered from sources on the verge of widespread use such as the internet, wearables, mobile phone apps, electronic health records, or genome sequencing. To utilise their full potential, we offer guidance by outlining available sources and approaches for the processing of secondary data. Furthermore, in addition to indicators for the regulatory and ethical evaluation of strategies for the best use of secondary data, we also propose criteria for assessing reusability. This overview supports more precise and effective policy decision making leading to earlier detection and better prevention of emerging health threats than is currently the case.

Sequence diversity of hepatitis D virus in Mongolia

Introduction

The Hepatitis Delta Virus (HDV) is a defective, single-stranded RNA virusoid encoding for a single protein, the Hepatitis Delta Antigen (HDAg), which requires the hepatitis B virus (HBV) envelope protein (HBsAg) for its transmission. Currently, hepatitis D is the most aggressive form of viral hepatitis and treatment options are limited. Worldwide 12 million people are chronically infected with HDV being at high risk for progression to cirrhosis and development of liver cancer.

Objectives

Although it is well established that Mongolia is the country with the highest prevalence of HDV infections, the information on the molecular epidemiology and factors contributing to HDV sequence diversity are largely unclear. The aim of the study was to characterize the sequence diversity of HDV in rural areas from Mongolia and to determine the extent of HLA class I-associated selection pressure.

Patients and methods

From the HepMongolia cohort from rural areas in Mongolia, 451 HBsAg-positive individuals were selected and anti-HDV, HDV-RNA and the sequence of the large HDAg was determined. For all individuals the HLA class I locus was genotyped. Residues under selection pressure in the presence of individual HLA class I types were identified with the recently published analysis tool HAMdetector.

Results

Of 431 HBsAg positive patients, 281 were anti-HDV positive (65%), and HDV-RNA could be detected in 207 of 281 (74%) of patients. The complete large HDAg was successfully sequenced from 131 samples. Phylogenetic analysis revealed that all Mongolian HDV isolates belong to genotype 1, however, they separate into several different clusters without clear regional association. In turn, from phylogeny there is strong evidence for recent local transmission events. Importantly, we found multiple residues with strong support for HLA class I-associated selection pressure consistent with a functional CD8⁺ T cell response directed against HDV.

Conclusion

HDV isolates from Mongolia are highly diverse. The molecular epidemiology suggests circulation of multiple subtypes and provides evidence for ongoing recent transmissions.

Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic

Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.

Women in the European Virus Bioinformatics Center

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

Weekly Nowcasting of New COVID-19 Cases Using Past Viral Load Measurements

The rapid spread of the coronavirus disease COVID-19 has imposed clinical and financial burdens on hospitals and governments attempting to provide patients with medical care and implement disease-controlling policies. The transmissibility of the disease was shown to be correlated with the patient’s viral load, which can be measured during testing using the cycle threshold (Ct). Previous models have utilized Ct to forecast the trajectory of the spread, which can provide valuable information to better allocate resources and change policies. However, these models combined other variables specific to medical institutions or came in the form of compartmental models that rely on epidemiological assumptions, all of which could impose prediction uncertainties. In this study, we overcome these limitations using data-driven modeling that utilizes Ct and previous number of cases, two institution-independent variables. We collected three groups of patients (n = 6296, n = 3228, and n = 12,096) from different time periods to train, validate, and independently validate the models. We used three machine learning algorithms and three deep learning algorithms that can model the temporal dynamic behavior of the number of cases. The endpoint was 7-week forward number of cases, and the prediction was evaluated using mean square error (MSE). The sequence-to-sequence model showed the best prediction during validation (MSE = 0.025), while polynomial regression (OLS) and support vector machine regression (SVR) had better performance during independent validation (MSE = 0.1596, and MSE = 0.16754, respectively), which exhibited better generalizability of the latter. The OLS and SVR models were used on a dataset from an external institution and showed promise in predicting COVID-19 incidences across institutions. These models may support clinical and logistic decision-making after prospective validation.

Advancing Precision Vaccinology by Molecular and Genomic Surveillance of Severe Acute Respiratory Syndrome Coronavirus 2 in Germany, 2021

BACKGROUND

Comprehensive pathogen genomic surveillance represents a powerful tool to complement and advance precision vaccinology. The emergence of the Alpha variant in December 2020 and the resulting efforts to track the spread of this and other severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern led to an expansion of genomic sequencing activities in Germany.

METHODS

At Robert Koch Institute (RKI), the German National Institute of Public Health, we established the Integrated Molecular Surveillance for SARS-CoV-2 (IMS-SC2) network to perform SARS-CoV-2 genomic surveillance at the national scale, SARS-CoV-2-positive samples from laboratories distributed across Germany regularly undergo whole-genome sequencing at RKI.

RESULTS

We report analyses of 3623 SARS-CoV-2 genomes collected between December 2020 and December 2021, of which 3282 were randomly sampled. All variants of concern were identified in the sequenced sample set, at ratios equivalent to those in the 100-fold larger German GISAID sequence dataset from the same time period. Phylogenetic analysis confirmed variant assignments. Multiple mutations of concern emerged during the observation period. To model vaccine effectiveness in vitro, we employed authentic-virus neutralization assays, confirming that both the Beta and Zeta variants are capable of immune evasion. The IMS-SC2 sequence dataset facilitated an estimate of the SARS-CoV-2 incidence based on genetic evolution rates. Together with modeled vaccine efficacies, Delta-specific incidence estimation indicated that the German vaccination campaign contributed substantially to a deceleration of the nascent German Delta wave.

CONCLUSIONS

SARS-CoV-2 molecular and genomic surveillance may inform public health policies including vaccination strategies and enable a proactive approach to controlling coronavirus disease 2019 spread as the virus evolves.

When might host heterogeneity drive the evolution of asymptomatic, pandemic coronaviruses?

Controlling many infectious diseases, including SARS-Coronavirus-2 (SARS-CoV-2), requires surveillance followed by isolation, contact-tracing and quarantining. These interventions often begin by identifying symptomatic individuals. However, actively removing pathogen strains causing symptomatic infections may inadvertently select for strains less likely to cause symptomatic infections. Moreover, a pathogen's fitness landscape is structured around a heterogeneous host pool; uneven surveillance efforts and distinct transmission risks across host classes can meaningfully alter selection pressures. Here, we explore this interplay between evolution caused by disease control efforts and the evolutionary consequences of host heterogeneity. Using an evolutionary epidemiology model parameterized for coronaviruses, we show that intense symptoms-driven disease control selects for asymptomatic strains, particularly when these efforts are applied unevenly across host groups. Under these conditions, increasing quarantine efforts have diverging effects. If isolation alone cannot eradicate, intensive quarantine efforts combined with uneven detections of asymptomatic infections (e.g., via neglect of some host classes) can favor the evolution of asymptomatic strains. We further show how, when intervention intensity depends on the prevalence of symptomatic infections, higher removal efforts (and isolating symptomatic cases in particular) more readily select for asymptomatic strains than when these efforts do not depend on prevalence. The selection pressures on pathogens caused by isolation and quarantining likely lie between the extremes of no intervention and thoroughly successful eradication. Thus, analyzing how different public health responses can select for asymptomatic pathogen strains is critical for identifying disease suppression efforts that can effectively manage emerging infectious diseases.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11071-022-07548-7.

Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data

BACKGROUND

As a new epi-center of COVID-19 in Asia and a densely populated developing country, Indonesia is facing unprecedented challenges in public health. SARS-CoV-2 lineage B.1.466.2 was reported to be an indigenous dominant strain in Indonesia (once second only to the Delta variant). However, it remains unclear how this variant evolved and spread within such an archipelagic nation.

METHODS

For statistical description, the spatiotemporal distributions of the B.1.466.2 variant were plotted using the publicly accessible metadata in GISAID. A total of 1302 complete genome sequences of Indonesian B.1.466.2 strains with high coverage were downloaded from the GISAID's EpiCoV database on 28 August 2021. To determine the molecular evolutionary characteristics, we performed a time-scaled phylogenetic analysis using the maximum likelihood algorithm and called the single nucleotide variants taking the Wuhan-Hu-1 sequence as reference. To investigate the spatiotemporal transmission patterns, we estimated two dynamic parameters (effective population size and effective reproduction number) and reconstructed the phylogeography among different islands.

RESULTS

As of the end of August 2021, nearly 85% of the global SARS-CoV-2 lineage B.1.466.2 sequences (including the first one) were obtained from Indonesia. This variant was estimated to account for over 50% of Indonesia's daily infections during the period of March-May 2021. The time-scaled phylogeny suggested that SARS-CoV-2 lineage B.1.466.2 circulating in Indonesia might have originated from Java Island in mid-June 2020 and had evolved into two disproportional and distinct sub-lineages. High-frequency non-synonymous mutations were mostly found in the spike and NSP3; the S-D614G/N439K/P681R co-mutations were identified in its larger sub-lineage. The demographic history was inferred to have experienced four phases, with an exponential growth from October 2020 to February 2021. The effective reproduction number was estimated to have reached its peak (11.18) in late December 2020 and dropped to be less than one after early May 2021. The relevant phylogeography showed that Java and Sumatra might successively act as epi-centers and form a stable transmission loop. Additionally, several long-distance transmission links across seas were revealed.

CONCLUSIONS

SARS-CoV-2 variants circulating in the tropical archipelago may follow unique patterns of evolution and transmission. Continuous, extensive and targeted genomic surveillance is essential.

A Bioinformatics Tool for Predicting Future COVID-19 Waves Based on a Retrospective Analysis of the Second Wave in India: Model Development Study

BACKGROUND

Since the start of the COVID-19 pandemic, health policymakers globally have been attempting to predict an impending wave of COVID-19. India experienced a devastating second wave of COVID-19 in the late first week of May 2021. We retrospectively analyzed the viral genomic sequences and epidemiological data reflecting the emergence and spread of the second wave of COVID-19 in India to construct a prediction model.

OBJECTIVE

We aimed to develop a bioinformatics tool that can predict an impending COVID-19 wave.

METHODS

We analyzed the time series distribution of genomic sequence data for SARS-CoV-2 and correlated it with epidemiological data for new cases and deaths for the corresponding period of the second wave. In addition, we analyzed the phylodynamics of circulating SARS-CoV-2 variants in the Indian population during the study period.

RESULTS

Our prediction analysis showed that the first signs of the arrival of the second wave could be seen by the end of January 2021, about 2 months before its peak in May 2021. By the end of March 2021, it was distinct. B.1.617 lineage variants powered the wave, most notably B.1.617.2 (Delta variant).

CONCLUSIONS

Based on the observations of this study, we propose that genomic surveillance of SARS-CoV-2 variants, complemented with epidemiological data, can be a promising tool to predict impending COVID-19 waves.

Deciphering early-warning signals of SARS-CoV-2 elimination and resurgence from limited data at multiple scales

Inferring the transmission potential of an infectious disease during low-incidence periods following epidemic waves is crucial for preparedness. In such periods, scarce data may hinder existing inference methods, blurring early-warning signals essential for discriminating between the likelihoods of resurgence versus elimination. Advanced insight into whether elevating caseloads (requiring swift community-wide interventions) or local elimination (allowing controls to be relaxed or refocussed on case-importation) might occur can separate decisive from ineffective policy. By generalizing and fusing recent approaches, we propose a novel early-warning framework that maximizes the information extracted from low-incidence data to robustly infer the chances of sustained local transmission or elimination in real time, at any scale of investigation (assuming sufficiently good surveillance). Applying this framework, we decipher hidden disease-transmission signals in prolonged low-incidence COVID-19 data from New Zealand, Hong Kong and Victoria, Australia. We uncover how timely interventions associate with averting resurgent waves, support official elimination declarations and evidence the effectiveness of the rapid, adaptive COVID-19 responses employed in these regions.