Phylodynamics reveals the role of human travel and contact tracing in controlling the first wave of COVID-19 in four island nations

Deep Learning and Likelihood Approaches for Viral Phylogeography Converge on the Same Answers Whether the Inference Model Is Right or Wrong

Analysis of phylogenetic trees has become an essential tool in epidemiology. Likelihood-based methods fit models to phylogenies to draw inferences about the phylodynamics and history of viral transmission. However, these methods are often computationally expensive, which limits the complexity and realism of phylodynamic models and makes them ill-suited for informing policy decisions in real-time during rapidly developing outbreaks. Likelihood-free methods using deep learning are pushing the boundaries of inference beyond these constraints. In this paper, we extend, compare, and contrast a recently developed deep learning method for likelihood-free inference from trees. We trained multiple deep neural networks using phylogenies from simulated outbreaks that spread among 5 locations and found they achieve close to the same levels of accuracy as Bayesian inference under the true simulation model. We compared robustness to model misspecification of a trained neural network to that of a Bayesian method. We found that both models had comparable performance, converging on similar biases. We also implemented a method of uncertainty quantification called conformalized quantile regression that we demonstrate has similar patterns of sensitivity to model misspecification as Bayesian highest posterior density (HPD) and greatly overlap with HPDs, but have lower precision (more conservative). Finally, we trained and tested a neural network against phylogeographic data from a recent study of the SARS-Cov-2 pandemic in Europe and obtained similar estimates of region-specific epidemiological parameters and the location of the common ancestor in Europe. Along with being as accurate and robust as likelihood-based methods, our trained neural networks are on average over 3 orders of magnitude faster after training. Our results support the notion that neural networks can be trained with simulated data to accurately mimic the good and bad statistical properties of the likelihood functions of generative phylogenetic models.

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.

Genomic surveillance reveals dynamic shifts in the connectivity of COVID-19 epidemics

The maturation of genomic surveillance in the past decade has enabled tracking of the emergence and spread of epidemics at an unprecedented level. During the COVID-19 pandemic, for example, genomic data revealed that local epidemics varied considerably in the frequency of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage importation and persistence, likely due to a combination of COVID-19 restrictions and changing connectivity. Here, we show that local COVID-19 epidemics are driven by regional transmission, including across international boundaries, but can become increasingly connected to distant locations following the relaxation of public health interventions. By integrating genomic, mobility, and epidemiological data, we find abundant transmission occurring between both adjacent and distant locations, supported by dynamic mobility patterns. We find that changing connectivity significantly influences local COVID-19 incidence. Our findings demonstrate a complex meaning of "local" when investigating connected epidemics and emphasize the importance of collaborative interventions for pandemic prevention and mitigation.

PhyloJunction: a computational framework for simulating, developing, and teaching evolutionary models

We introduce PhyloJunction, a computational framework designed to facilitate the prototyping, testing, and characterization of evolutionary models. PhyloJunction is distributed as an open-source Python library that can be used to implement a variety of models, through its flexible graphical modeling architecture and dedicated model specification language. Model design and use are exposed to users via command-line and graphical interfaces, which integrate the steps of simulating, summarizing, and visualizing data. This paper describes the features of PhyloJunction - which include, but are not limited to, a general implementation of a popular family of phylogenetic diversification models - and, moving forward, how it may be expanded to not only include new models, but to also become a platform for conducting and teaching statistical learning.

Leveraging insect-specific viruses to elucidate mosquito population structure and dynamics

Several aspects of mosquito ecology that are important for vectored disease transmission and control have been difficult to measure at epidemiologically important scales in the field. In particular, the ability to describe mosquito population structure and movement rates has been hindered by difficulty in quantifying fine-scale genetic variation among populations. The mosquito virome represents a possible avenue for quantifying population structure and movement rates across multiple spatial scales. Mosquito viromes contain a diversity of viruses, including several insect-specific viruses (ISVs) and "core" viruses that have high prevalence across populations. To date, virome studies have focused on viral discovery and have only recently begun examining viral ecology. While nonpathogenic ISVs may be of little public health relevance themselves, they provide a possible route for quantifying mosquito population structure and dynamics. For example, vertically transmitted viruses could behave as a rapidly evolving extension of the host's genome. It should be possible to apply established analytical methods to appropriate viral phylogenies and incidence data to generate novel approaches for estimating mosquito population structure and dispersal over epidemiologically relevant timescales. By studying the virome through the lens of spatial and genomic epidemiology, it may be possible to investigate otherwise cryptic aspects of mosquito ecology. A better understanding of mosquito population structure and dynamics are key for understanding mosquito-borne disease ecology and methods based on ISVs could provide a powerful tool for informing mosquito control programs.

LinguaPhylo: A probabilistic model specification language for reproducible phylogenetic analyses

Phylogenetic models have become increasingly complex, and phylogenetic data sets have expanded in both size and richness. However, current inference tools lack a model specification language that can concisely describe a complete phylogenetic analysis while remaining independent of implementation details. We introduce a new lightweight and concise model specification language, 'LPhy', which is designed to be both human and machine-readable. A graphical user interface accompanies 'LPhy', allowing users to build models, simulate data, and create natural language narratives describing the models. These narratives can serve as the foundation for manuscript method sections. Additionally, we present a command-line interface for converting LPhy-specified models into analysis specification files (in XML format) compatible with the BEAST2 software platform. Collectively, these tools aim to enhance the clarity of descriptions and reporting of probabilistic models in phylogenetic studies, ultimately promoting reproducibility of results.

Genomic epidemiology of SARS-CoV-2 variants during the first two years of the pandemic in Colombia

BACKGROUND

The emergence of highly transmissible SARS-CoV-2 variants has led to surges in cases and the need for global genomic surveillance. While some variants rapidly spread worldwide, other variants only persist nationally. There is a need for more fine-scale analysis to understand transmission dynamics at a country scale. For instance, the Mu variant of interest, also known as lineage B.1.621, was first detected in Colombia and was responsible for a large local wave but only a few sporadic cases elsewhere.

METHODS

To better understand the epidemiology of SARS-Cov-2 variants in Colombia, we used 14,049 complete SARS-CoV-2 genomes from the 32 states of Colombia. We performed Bayesian phylodynamic analyses to estimate the time of variants' introduction, their respective effective reproductive number, and effective population size, and the impact of disease control measures.

RESULTS

Here, we detect a total of 188 SARS-CoV-2 Pango lineages circulating in Colombia since the pandemic's start. We show that the effective reproduction number oscillated drastically throughout the first two years of the pandemic, with Mu showing the highest transmissibility (Re and growth rate estimation).

CONCLUSIONS

Our results reinforce that genomic surveillance programs are essential for countries to make evidence-driven interventions toward the emergence and circulation of novel SARS-CoV-2 variants.

COVID-19 impact on EuroTravNet infectious diseases sentinel surveillance in Europe

BACKGROUND

The COVID-19 pandemic resulted in a sharp decline of post-travel patient encounters at the European sentinel surveillance network (EuroTravNet) of travellers' health. We report on the impact of COVID-19 on travel-related infectious diseases as recorded by EuroTravNet clinics.

METHODS

Travelers who presented between January 1, 2019 and September 30, 2021 were included. Comparisons were made between the pre-pandemic period (14 months from January 1, 2019 to February 29, 2020); and the pandemic period (19 months from March 1, 2020 to September 30, 2021).

RESULTS

Of the 15,124 visits to the network during the 33-month observation period, 10,941 (72%) were during the pre-pandemic period, and 4183 (28%) during the pandemic period. Average monthly visits declined from 782/month (pre-COVID-19 era) to 220/month (COVID-19 pandemic era). Among non-migrants the top-10 countries of exposure changed after onset of the COVID-19 pandemic; destinations such as Italy and Austria, where COVID-19 exposure peaked in the first months, replaced typical travel destinations in Asia (Thailand, Indonesia, India). There was a small decline in migrant patients reported, with little change in the top countries of exposure (Bolivia, Mali). The three top diagnoses with the largest overall decreases in relative frequency were acute gastroenteritis (-5.3%), rabies post-exposure prophylaxis (-2.8%), and dengue (-2.6%). Apart from COVID-19 (which rose from 0.1% to 12.7%), the three top diagnoses with the largest overall relative frequency increase were schistosomiasis (+4.9%), strongyloidiasis (+2.7%), and latent tuberculosis (+2.4%).

CONCLUSIONS

A marked COVID-19 pandemic-induced decline in global travel activities is reflected in reduced travel-related infectious diseases sentinel surveillance reporting.

Swiss public health measures associated with reduced SARS-CoV-2 transmission using genome data

Genome sequences from evolving infectious pathogens allow quantification of case introductions and local transmission dynamics. We sequenced 11,357 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from Switzerland in 2020-the sixth largest effort globally. Using a representative subset of these data, we estimated viral introductions to Switzerland and their persistence over the course of 2020. We contrasted these estimates with simple null models representing the absence of certain public health measures. We show that Switzerland's border closures decoupled case introductions from incidence in neighboring countries. Under a simple model, we estimate an 86 to 98% reduction in introductions during Switzerland's strictest border closures. Furthermore, the Swiss 2020 partial lockdown roughly halved the time for sampled introductions to die out. Last, we quantified local transmission dynamics once introductions into Switzerland occurred using a phylodynamic model. We found that transmission slowed 35 to 63% upon outbreak detection in summer 2020 but not in fall. This finding may indicate successful contact tracing over summer before overburdening in fall. The study highlights the added value of genome sequencing data for understanding transmission dynamics.

Genome Structure, Life Cycle, and Taxonomy of Coronaviruses and the Evolution of SARS-CoV-2

Coronaviruses have a broad host range and exhibit high zoonotic potential. In this chapter, we describe their genomic organization in terms of encoded proteins and provide an introduction to the peculiar discontinuous transcription mechanism. Further, we present evolutionary conserved genomic RNA secondary structure features, which are involved in the complex replication mechanism. With a focus on computational methods, we review the emergence of SARS-CoV-2 starting with the 2019 strains. In that context, we also discuss the debated hypothesis of whether SARS-CoV-2 was created in a laboratory. We focus on the molecular evolution and the epidemiological dynamics of this recently emerged pathogen and we explain how variants of concern are detected and characterised. COVID-19, the disease caused by SARS-CoV-2, can spread through different transmission routes and also depends on a number of risk factors. We describe how current computational models of viral epidemiology, or more specifically, phylodynamics, have facilitated and will continue to enable a better understanding of the epidemic dynamics of SARS-CoV-2.

Tracing the international arrivals of SARS-CoV-2 Omicron variants after Aotearoa New Zealand reopened its border

In the second quarter of 2022, there was a global surge of emergent SARS-CoV-2 lineages that had a distinct growth advantage over then-dominant Omicron BA.1 and BA.2 lineages. By generating 10,403 Omicron genomes, we show that Aotearoa New Zealand observed an influx of these immune-evasive variants (BA.2.12.1, BA.4, and BA.5) through the border. This is explained by the return to significant levels of international travel following the border's reopening in March 2022. We estimate one Omicron transmission event from the border to the community for every ~5,000 passenger arrivals at the current levels of travel and restriction. Although most of these introductions did not instigate any detected onward transmission, a small minority triggered large outbreaks. Genomic surveillance at the border provides a lens on the rate at which new variants might gain a foothold and trigger new waves of infection.

A systematic review of COVID-19 transport policies and mitigation strategies around the globe

This paper reports a Scopus-based systematic literature review of a wide variety of transportation policies and mitigation strategies that have been conducted around the world to minimize COVID-19 contagion risk in transportation systems. The review offers a representative coverage of countries across all continents of the planet, as well as among representative climate regions - as weather is an important factor to consider. The readership interested in policies and mitigation strategies is expected to involve a wide range of actors, each involving a particular application context; hence, the literature is also characterized by key attributes such as: transportation mode; actor (users, operators, government, industry); jurisdiction (national, provincial, city, neighborhood); and area of application (planning, regulation, operations, research, incentives). An in-depth analysis of the surveyed literature is then reported, focusing first on condensing the literature into 151 distinct policies and strategies, which are subsequently categorized into 25 broad categories that are discussed at length. The compendium and discussion of strategies and policies reported not only provide comprehensive guidelines to inform various courses of action for decision-makers, planners, and social communicators, but also emphasize on future work and the potential of some of these strategies to be the precursors of meaningful, more sustainable behavioral changes in future mobility patterns.

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.

New Phylogenetic Models Incorporating Interval-Specific Dispersal Dynamics Improve Inference of Disease Spread

Phylodynamic methods reveal the spatial and temporal dynamics of viral geographic spread, and have featured prominently in studies of the COVID-19 pandemic. Virtually all such studies are based on phylodynamic models that assume-despite direct and compelling evidence to the contrary-that rates of viral geographic dispersal are constant through time. Here, we: (1) extend phylodynamic models to allow both the average and relative rates of viral dispersal to vary independently between pre-specified time intervals; (2) implement methods to infer the number and timing of viral dispersal events between areas; and (3) develop statistics to assess the absolute fit of discrete-geographic phylodynamic models to empirical datasets. We first validate our new methods using simulations, and then apply them to a SARS-CoV-2 dataset from the early phase of the COVID-19 pandemic. We show that: (1) under simulation, failure to accommodate interval-specific variation in the study data will severely bias parameter estimates; (2) in practice, our interval-specific discrete-geographic phylodynamic models can significantly improve the relative and absolute fit to empirical data; and (3) the increased realism of our interval-specific models provides qualitatively different inferences regarding key aspects of the COVID-19 pandemic-revealing significant temporal variation in global viral dispersal rates, viral dispersal routes, and the number of viral dispersal events between areas-and alters interpretations regarding the efficacy of intervention measures to mitigate the pandemic.

Minimising the use of costly control measures in an epidemic elimination strategy: A simple mathematical model

Countries such as New Zealand, Australia and Taiwan responded to the Covid-19 pandemic with an elimination strategy. This involves a combination of strict border controls with a rapid and effective response to eliminate border-related re-introductions. An important question for decision makers is, when there is a new re-introduction, what is the right threshold at which to implement strict control measures designed to reduce the effective reproduction number below 1. Since it is likely that there will be multiple re-introductions, responding at too low a threshold may mean repeatedly implementing controls unnecessarily for outbreaks that would self-eliminate even without control measures. On the other hand, waiting for too high a threshold to be reached creates a risk that controls will be needed for a longer period of time, or may completely fail to contain the outbreak. Here, we use a highly idealised branching process model of small border-related outbreaks to address this question. We identify important factors that affect the choice of threshold in order to minimise the expect time period for which control measures are in force. We find that the optimal threshold for introducing controls decreases with the effective reproduction number, and increases with overdispersion of the offspring distribution and with the effectiveness of control measures. Our results are not intended as a quantitative decision-making algorithm. However, they may help decision makers understand when a wait-and-see approach is likely to be preferable over an immediate response.

Genomic epidemiology of Delta SARS-CoV-2 during transition from elimination to suppression in Aotearoa New Zealand

New Zealand's COVID-19 elimination strategy heavily relied on the use of genomics to inform contact tracing, linking cases to the border and to clusters during community outbreaks. In August 2021, New Zealand entered its second nationwide lockdown after the detection of a single community case with no immediately apparent epidemiological link to the border. This incursion resulted in the largest outbreak seen in New Zealand caused by the Delta Variant of Concern. Here we generated 3806 high quality SARS-CoV-2 genomes from cases reported in New Zealand between 17 August and 1 December 2021, representing 43% of reported cases. We detected wide geographical spread coupled with undetected community transmission, characterised by the apparent extinction and reappearance of genomically linked clusters. We also identified the emergence, and near replacement, of genomes possessing a 10-nucleotide frameshift deletion that caused the likely truncation of accessory protein ORF7a. By early October, New Zealand moved from an elimination strategy to a suppression strategy and the role of genomics changed markedly from being used to track and trace, towards population-level surveillance.

The End of the Elimination Strategy: Decisive Factors towards Sustainable Management of COVID-19 in New Zealand

New Zealand has long been praised for the effectiveness of its COVID-19 elimination strategy. It resulted in fewer COVID-19-related deaths, better economic recovery, and less stringent policy measures within its borders compared with other OECD countries, which opted for mitigation or suppression. However, since September 2021, the rising number of infections has not been contained anymore by the contact tracing and self-isolation system in place and the government has shifted towards a policy strategy similar to suppression to manage the crisis. In this case study, we analyse the factors that led the government to switch policy and discuss why elimination became unsustainable to manage the COVID-19 epidemic in New Zealand. Results showed that the socioeconomic and political factors, along with the appearance of new variants and a delayed vaccination program, were accountable for the switch in strategy. This switch allows the country to better adapt to the evolving nature of the disease and to address the social and economic repercussions of the first year of measures. Our conclusion does not disregard elimination as an appropriate initial strategy to contain this pandemic in the absence of a vaccine or treatment, but rather suggests that borders cannot remain closed for long periods of time without creating social, economical, and political issues.

An Efficient Coalescent Epoch Model for Bayesian Phylogenetic Inference

We present a two-headed approach called Bayesian Integrated Coalescent Epoch PlotS (BICEPS) for efficient inference of coalescent epoch models. Firstly, we integrate out population size parameters, and secondly, we introduce a set of more powerful Markov chain Monte Carlo (MCMC) proposals for flexing and stretching trees. Even though population sizes are integrated out and not explicitly sampled through MCMC, we are still able to generate samples from the population size posteriors. This allows demographic reconstruction through time and estimating the timing and magnitude of population bottlenecks and full population histories. Altogether, BICEPS can be considered a more muscular version of the popular Bayesian skyline model. We demonstrate its power and correctness by a well-calibrated simulation study. Furthermore, we demonstrate with an application to SARS-CoV-2 genomic data that some analyses that have trouble converging with the traditional Bayesian skyline prior and standard MCMC proposals can do well with the BICEPS approach. BICEPS is available as open-source package for BEAST 2 under GPL license and has a user-friendly graphical user interface.[Bayesian phylogenetics; BEAST 2; BICEPS; coalescent model.].

The impact of global lineage dynamics, border restrictions, and emergence of the B.1.1.7 lineage on the SARS-CoV-2 epidemic in Norway

As the COVID-19 pandemic swept through an immunologically naïve human population, academics and public health professionals scrambled to establish methods and platforms for genomic surveillance and data sharing. This offered a rare opportunity to study the ecology and evolution of SARS-CoV-2 over the course of the ongoing pandemic. Here, we use population genetic and phylogenetic methodology to characterize the population dynamics of SARS-CoV-2 and reconstruct patterns of virus introductions and local transmission in Norway against this backdrop. The analyses demonstrated that the epidemic in Norway was largely import driven and characterized by the repeated introduction, establishment, and suppression of new transmission lineages. This pattern changed with the arrival of the B.1.1.7 lineage, which was able to establish a stable presence concomitant with the imposition of severe border restrictions.

Progress and challenges in virus genomic epidemiology

Genomic epidemiology, which links pathogen genomes with associated metadata to understand disease transmission, has become a key component of outbreak response. Decreasing costs of genome sequencing and increasing computational power provide opportunities to generate and analyse large viral genomic datasets that aim to uncover the spatial scales of transmission, the demographics contributing to transmission patterns, and to forecast epidemic trends. Emerging sources of genomic data and associated metadata provide new opportunities to further unravel transmission patterns. Key challenges include how to integrate genomic data with metadata from multiple sources, how to generate efficient computational algorithms to cope with large datasets, and how to establish sampling frameworks to enable robust conclusions.

Real-Time Genomics for Tracking Severe Acute Respiratory Syndrome Coronavirus 2 Border Incursions after Virus Elimination, New Zealand

Since severe acute respiratory syndrome coronavirus 2 was first eliminated in New Zealand in May 2020, a total of 13 known coronavirus disease (COVID-19) community outbreaks have occurred, 2 of which led health officials to issue stay-at-home orders. These outbreaks originated at the border via isolating returnees, airline workers, and cargo vessels. Because a public health system was informed by real-time viral genomic sequencing and complete genomes typically were available within 12 hours of community-based positive COVID-19 test results, every outbreak was well-contained. A total of 225 community cases resulted in 3 deaths. Real-time genomics were essential for establishing links between cases when epidemiologic data could not do so and for identifying when concurrent outbreaks had different origins.