Hamiltonian Monte Carlo sampling to estimate past population dynamics using the skygrid coalescent model in a Bayesian phylogenetics framework

Untangling lineage introductions, persistence and transmission drivers of HP-PRRSV sublineage 8.7

Despite a rapid expansion of Porcine reproductive and respiratory syndrome virus (PRRSV) sublineage 8.7 over recent years, very little is known about the patterns of virus evolution, dispersal, and the factors influencing this dispersal. Relying on a national PRRSV surveillance project established over 20 years ago, we expand the available genomic data of sublineage 8.7 from China. We perform independent interlineage and intralineage recombination analyses for the entire study period, which showed a heterogeneous recombination pattern. A series of Bayesian phylogeographic analyses uncover the role of Guangdong as an important infection hub within Asia. The spatial spread of PRRSV is highly linked with a composite of human activities and the heterogeneous provincial distribution of the swine industry, largely propelled by the smaller-scale Chinese rural farming systems in the past years. We sequence all four available modified live vaccines (MLVs) and perform genomic analyses with publicly available data, of which our results suggest a key "leaky" period spanning 2011-2017 with two concurrent amino acid mutations in ORF1a 957 and ORF2 250. Overall, our study provides an in-depth overview of the evolution, transmission dynamics, and potential leaky status of HP-PRRS MLVs, providing critical insights into new MLV development.

Emergence of the B.1.214.2 SARS-CoV-2 lineage with an Omicron-like spike insertion and a unique upper airway immune signature

We investigate the emergence, mutation profile, and dissemination of SARS-CoV-2 lineage B.1.214.2, first identified in Belgium in January 2021. This variant, featuring a 3-amino acid insertion in the spike protein similar to the Omicron variant, was speculated to enhance transmissibility or immune evasion. Initially detected in international travelers, it substantially transmitted in Central Africa, Belgium, Switzerland, and France, peaking in April 2021. Our travel-aware phylogeographic analysis, incorporating travel history, estimated the origin to the Republic of the Congo, with primary European entry through France and Belgium, and multiple smaller introductions during the epidemic. We correlate its spread with human travel patterns and air passenger data. Further, upon reviewing national reports of SARS-CoV-2 outbreaks in Belgian nursing homes, we found this strain caused moderately severe outcomes (8.7% case fatality ratio). A distinct nasopharyngeal immune response was observed in elderly patients, characterized by 80% unique signatures, higher B- and T-cell activation, increased type I IFN signaling, and reduced NK, Th17, and complement system activation, compared to similar outbreaks. This unique immune response may explain the variant's epidemiological behavior and underscores the need for nasal vaccine strategies against emerging variants.

Travel surveillance uncovers dengue virus dynamics and introductions in the Caribbean

Dengue is the most prevalent mosquito-borne viral disease in humans, and cases are continuing to rise globally. In particular, islands in the Caribbean have experienced more frequent outbreaks, and all four dengue virus (DENV) serotypes have been reported in the region, leading to hyperendemicity and increased rates of severe disease. However, there is significant variability regarding virus surveillance and reporting between islands, making it difficult to obtain an accurate understanding of the epidemiological patterns in the Caribbean. To investigate this, we used travel surveillance and genomic epidemiology to reconstruct outbreak dynamics, DENV serotype turnover, and patterns of spread within the region from 2009-2022. We uncovered two recent DENV-3 introductions from Asia, one of which resulted in a large outbreak in Cuba, which was previously under-reported. We also show that while outbreaks can be synchronized between islands, they are often caused by different serotypes. Our study highlights the importance of surveillance of infected travelers to provide a snapshot of local introductions and transmission in areas with limited local surveillance and suggests that the recent DENV-3 introductions may pose a major public health threat in the region.

Scalable gradients enable Hamiltonian Monte Carlo sampling for phylodynamic inference under episodic birth-death-sampling models

Birth-death models play a key role in phylodynamic analysis for their interpretation in terms of key epidemiological parameters. In particular, models with piecewise-constant rates varying at different epochs in time, to which we refer as episodic birth-death-sampling (EBDS) models, are valuable for their reflection of changing transmission dynamics over time. A challenge, however, that persists with current time-varying model inference procedures is their lack of computational efficiency. This limitation hinders the full utilization of these models in large-scale phylodynamic analyses, especially when dealing with high-dimensional parameter vectors that exhibit strong correlations. We present here a linear-time algorithm to compute the gradient of the birth-death model sampling density with respect to all time-varying parameters, and we implement this algorithm within a gradient-based Hamiltonian Monte Carlo (HMC) sampler to alleviate the computational burden of conducting inference under a wide variety of structures of, as well as priors for, EBDS processes. We assess this approach using three different real world data examples, including the HIV epidemic in Odesa, Ukraine, seasonal influenza A/H3N2 virus dynamics in New York state, America, and Ebola outbreak in West Africa. HMC sampling exhibits a substantial efficiency boost, delivering a 10- to 200-fold increase in minimum effective sample size per unit-time, in comparison to a Metropolis-Hastings-based approach. Additionally, we show the robustness of our implementation in both allowing for flexible prior choices and in modeling the transmission dynamics of various pathogens by accurately capturing the changing trend of viral effective reproductive number.

Dispersion patterns of SARS-CoV-2 variants Gamma, Lambda and Mu in Latin America and the Caribbean

Latin America and Caribbean (LAC) regions were an important epicenter of the COVID-19 pandemic and SARS-CoV-2 evolution. Through the COVID-19 Genomic Surveillance Regional Network (COVIGEN), LAC countries produced an important number of genomic sequencing data that made possible an enhanced SARS-CoV-2 genomic surveillance capacity in the Americas, paving the way for characterization of emerging variants and helping to guide the public health response. In this study we analyzed approximately 300,000 SARS-CoV-2 sequences generated between February 2020 and March 2022 by multiple genomic surveillance efforts in LAC and reconstructed the diffusion patterns of the main variants of concern (VOCs) and of interest (VOIs) possibly originated in the Region. Our phylogenetic analysis revealed that the spread of variants Gamma, Lambda and Mu reflects human mobility patterns due to variations of international air passenger transportation and gradual lifting of social distance measures previously implemented in countries. Our results highlight the potential of genetic data to reconstruct viral spread and unveil preferential routes of viral migrations that are shaped by human mobility patterns.

Travel surveillance uncovers dengue virus dynamics and introductions in the Caribbean

Dengue is the most prevalent mosquito-borne viral disease in humans, and cases are continuing to rise globally. In particular, islands in the Caribbean have experienced more frequent outbreaks, and all four dengue virus (DENV) serotypes have been reported in the region, leading to hyperendemicity and increased rates of severe disease. However, there is significant variability regarding virus surveillance and reporting between islands, making it difficult to obtain an accurate understanding of the epidemiological patterns in the Caribbean. To investigate this, we used travel surveillance and genomic epidemiology to reconstruct outbreak dynamics, DENV serotype turnover, and patterns of spread within the region from 2009-2022. We uncovered two recent DENV-3 introductions from Asia, one of which resulted in a large outbreak in Cuba, which was previously under-reported. We also show that while outbreaks can be synchronized between islands, they are often caused by different serotypes. Our study highlights the importance of surveillance of infected travelers to provide a snapshot of local introductions and transmission in areas with limited local surveillance and suggests that the recent DENV-3 introductions may pose a major public health threat in the region.

Scalable gradients enable Hamiltonian Monte Carlo sampling for phylodynamic inference under episodic birth-death-sampling models

Birth-death models play a key role in phylodynamic analysis for their interpretation in terms of key epidemiological parameters. In particular, models with piecewise-constant rates varying at different epochs in time, to which we refer as episodic birth-death-sampling (EBDS) models, are valuable for their reflection of changing transmission dynamics over time. A challenge, however, that persists with current time-varying model inference procedures is their lack of computational efficiency. This limitation hinders the full utilization of these models in large-scale phylodynamic analyses, especially when dealing with high-dimensional parameter vectors that exhibit strong correlations. We present here a linear-time algorithm to compute the gradient of the birth-death model sampling density with respect to all time-varying parameters, and we implement this algorithm within a gradient-based Hamiltonian Monte Carlo (HMC) sampler to alleviate the computational burden of conducting inference under a wide variety of structures of, as well as priors for, EBDS processes. We assess this approach using three different real world data examples, including the HIV epidemic in Odesa, Ukraine, seasonal influenza A/H3N2 virus dynamics in New York state, America, and Ebola outbreak in West Africa. HMC sampling exhibits a substantial efficiency boost, delivering a 10- to 200-fold increase in minimum effective sample size per unit-time, in comparison to a Metropolis-Hastings-based approach. Additionally, we show the robustness of our implementation in both allowing for flexible prior choices and in modeling the transmission dynamics of various pathogens by accurately capturing the changing trend of viral effective reproductive number.

Scalable Bayesian Divergence Time Estimation With Ratio Transformations

Divergence time estimation is crucial to provide temporal signals for dating biologically important events from species divergence to viral transmissions in space and time. With the advent of high-throughput sequencing, recent Bayesian phylogenetic studies have analyzed hundreds to thousands of sequences. Such large-scale analyses challenge divergence time reconstruction by requiring inference on highly correlated internal node heights that often become computationally infeasible. To overcome this limitation, we explore a ratio transformation that maps the original $N-1$ internal node heights into a space of one height parameter and $N-2$ ratio parameters. To make the analyses scalable, we develop a collection of linear-time algorithms to compute the gradient and Jacobian-associated terms of the log-likelihood with respect to these ratios. We then apply Hamiltonian Monte Carlo sampling with the ratio transform in a Bayesian framework to learn the divergence times in 4 pathogenic viruses (West Nile virus, rabies virus, Lassa virus, and Ebola virus) and the coralline red algae. Our method both resolves a mixing issue in the West Nile virus example and improves inference efficiency by at least 5-fold for the Lassa and rabies virus examples as well as for the algae example. Our method now also makes it computationally feasible to incorporate mixed-effects molecular clock models for the Ebola virus example, confirms the findings from the original study, and reveals clearer multimodal distributions of the divergence times of some clades of interest.

Genomic assessment of invasion dynamics of SARS-CoV-2 Omicron BA.1

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) now arise in the context of heterogeneous human connectivity and population immunity. Through a large-scale phylodynamic analysis of 115,622 Omicron BA.1 genomes, we identified >6,000 introductions of the antigenically distinct VOC into England and analyzed their local transmission and dispersal history. We find that six of the eight largest English Omicron lineages were already transmitting when Omicron was first reported in southern Africa (22 November 2021). Multiple datasets show that importation of Omicron continued despite subsequent restrictions on travel from southern Africa as a result of export from well-connected secondary locations. Initiation and dispersal of Omicron transmission lineages in England was a two-stage process that can be explained by models of the country's human geography and hierarchical travel network. Our results enable a comparison of the processes that drive the invasion of Omicron and other VOCs across multiple spatial scales.

Dynamics of eastern equine encephalitis virus during the 2019 outbreak in the Northeast United States

Eastern equine encephalitis virus (EEEV) causes a rare but severe disease in horses and humans and is maintained in an enzootic transmission cycle between songbirds and Culiseta melanura mosquitoes. In 2019, the largest EEEV outbreak in the United States for more than 50 years occurred, centered in the Northeast. To explore the dynamics of the outbreak, we sequenced 80 isolates of EEEV and combined them with existing genomic data. We found that, similar to previous years, cases were driven by multiple independent but short-lived virus introductions into the Northeast from Florida. Once in the Northeast, we found that Massachusetts was important for regional spread. We found no evidence of any changes in viral, human, or bird factors which would explain the increase in cases in 2019, although the ecology of EEEV is complex and further data is required to explore these in more detail. By using detailed mosquito surveillance data collected by Massachusetts and Connecticut, however, we found that the abundance of Cs. melanura was exceptionally high in 2019, as was the EEEV infection rate. We employed these mosquito data to build a negative binomial regression model and applied it to estimate early season risks of human or horse cases. We found that the month of first detection of EEEV in mosquito surveillance data and vector index (abundance multiplied by infection rate) were predictive of cases later in the season. We therefore highlight the importance of mosquito surveillance programs as an integral part of public health and disease control.

Dynamics of Eastern equine encephalitis virus during the 2019 outbreak in the Northeast United States

Eastern equine encephalitis virus (EEEV) causes a rare but severe disease in horses and humans, and is maintained in an enzootic transmission cycle between songbirds and Culiseta melanura mosquitoes. In 2019, the largest EEEV outbreak in the United States for more than 50 years occurred, centered in the Northeast. To explore the dynamics of the outbreak, we sequenced 80 isolates of EEEV and combined them with existing genomic data. We found that, like previous years, cases were driven by frequent short-lived virus introductions into the Northeast from Florida. Once in the Northeast, we found that Massachusetts was important for regional spread. We found no evidence of any changes in viral, human, or bird factors which would explain the increase in cases in 2019. By using detailed mosquito surveillance data collected by Massachusetts and Connecticut, however, we found that the abundance of Cs. melanura was exceptionally high in 2019, as was the EEEV infection rate. We employed these mosquito data to build a negative binomial regression model and applied it to estimate early season risks of human or horse cases. We found that the month of first detection of EEEV in mosquito surveillance data and vector index (abundance multiplied by infection rate) were predictive of cases later in the season. We therefore highlight the importance of mosquito surveillance programs as an integral part of public health and disease control.

Tracking the turnover of SARS-CoV-2 VOCs Gamma to Delta in a Brazilian state (Minas Gerais) with a high-vaccination status

The emergence and global dissemination of Severe Acute Respiratory Syndrome virus 2 (SARS-CoV-2) variants of concern (VOCs) have been described as the main factor driving the Coronavirus Disease 2019 pandemic. In Brazil, the Gamma variant dominated the epidemiological scenario during the first period of 2021. Many Brazilian regions detected the Delta variant after its first description and documented its spread. To monitor the introduction and spread of VOC Delta, we performed Polymerase Chain Reaction (PCR) genotyping and genome sequencing in ten regional sentinel units from June to October 2021 in the State of Minas Gerais (MG). We documented the introduction and spread of Delta, comprising 70 per cent of the cases 8 weeks later. Comparing the viral loads of the Gamma and Delta dominance periods, we provide additional evidence that the latter is more transmissible. The spread and dominance of Delta did not culminate in the increase in cases and deaths, suggesting that the vaccination may have restrained the epidemic growth. Analysis of 224 novel Delta genomes revealed that Rio de Janeiro state was the primary source for disseminating this variant in the state of MG. We present the establishment of Delta, providing evidence of its enhanced transmissibility and showing that this variant shift did not aggravate the epidemiological scenario in a high immunity setting.

Molecular Epidemiology of Porcine Parvovirus Type 1 (PPV1) and the Reactivity of Vaccine-Induced Antisera Against Historical and Current PPV1 Strains

Porcine Parvovirus Type 1 (PPV1) contributes to important losses in the swine industry worldwide. During a PPV1 infection, embryos and fetuses are targeted, resulting in stillbirth, mummification, embryonic death, and infertility (SMEDI syndrome). Even though vaccination is common in gilts and sows, strains mainly belonging to the 27a-like group have been spreading in Europe since early 2000s, resulting in SMEDI problems and requiring in-depth studies into the molecular epidemiology and vaccination efficacy of commercial vaccines. Here, we show that PPV1 has evolved since 1855 [1737, 1933] at a rate of 4.71 × 10⁻⁵ nucleotide substitutions per site per year. Extensive sequencing allowed evaluating and reassessing the current PPV1 VP1-based classifications, providing evidence for the existence of four relevant phylogenetic groups. While most European strains belong to the PPV1a (G1) or PPV1b (G2 or 27a-like) group, most Asian and American G2 strains and some European strains were divided into virulent PPV1c (e.g. NADL-8) and attenuated PPV1d (e.g. NADL-2) groups. The increase in the swine population, vaccination degree, and health management (vaccination and biosafety) influenced the spread of PPV1. The reactivity of anti-PPV1 antibodies from sows vaccinated with Porcilis^© Parvo, Eryseng^© Parvo, or ReproCyc^© ParvoFLEX against different PPV1 field strains was the highest upon vaccination with ReproCyc^© ParvoFLEX, followed by Eryseng^© Parvo, and Porcilis^© Parvo. Our findings contribute to the evaluation of the immunogenicity of existing vaccines and support the development of new vaccine candidates. Finally, the potential roles of cluster-specific hallmark amino acids in elevated pathogenicity and viral entry are discussed.

Antibody escape and global spread of SARS-CoV-2 lineage A.27

In spring 2021, an increasing number of infections was observed caused by the hitherto rarely described SARS-CoV-2 variant A.27 in south-west Germany. From December 2020 to June 2021 this lineage has been detected in 31 countries. Phylogeographic analyses of A.27 sequences obtained from national and international databases reveal a global spread of this lineage through multiple introductions from its inferred origin in Western Africa. Variant A.27 is characterized by a mutational pattern in the spike gene that includes the L18F, L452R and N501Y spike amino acid substitutions found in various variants of concern but lacks the globally dominant D614G. Neutralization assays demonstrate an escape of A.27 from convalescent and vaccine-elicited antibody-mediated immunity. Moreover, the therapeutic monoclonal antibody Bamlanivimab and partially the REGN-COV2 cocktail fail to block infection by A.27. Our data emphasize the need for continued global monitoring of novel lineages because of the independent evolution of new escape mutations.

Phylogeography reveals association between swine trade and the spread of porcine epidemic diarrhea virus in China and across the world

The ongoing SARS (severe acute respiratory syndrome)-CoV (coronavirus)-2 pandemic has exposed major gaps in our knowledge on the origin, ecology, evolution, and spread of animal coronaviruses. Porcine epidemic diarrhea virus (PEDV) is a member of the genus Alphacoronavirus in the family Coronaviridae that may have originated from bats and leads to significant hazards and widespread epidemics in the swine population. The role of local and global trade of live swine and swine-related products in disseminating PEDV remains unclear, especially in developing countries with complex swine production systems. Here, we undertake an in-depth phylogeographic analysis of PEDV sequence data (including 247 newly sequenced samples) and employ an extension of this inference framework that enables formally testing the contribution of a range of predictor variables to the geographic spread of PEDV. Within China, the provinces of Guangdong and Henan were identified as primary hubs for the spread of PEDV, for which we estimate live swine trade to play a very important role. On a global scale, the United States and China maintain the highest number of PEDV lineages. We estimate that, after an initial introduction out of China, the United States acted as an important source of PEDV introductions into Japan, Korea, China, and Mexico. Live swine trade also explains the dispersal of PEDV on a global scale. Given the increasingly global trade of live swine, our findings have important implications for designing prevention and containment measures to combat a wide range of livestock coronaviruses.

Emergence and spread of SARS-CoV-2 lineage B.1.620 with variant of concern-like mutations and deletions

Distinct SARS-CoV-2 lineages, discovered through various genomic surveillance initiatives, have emerged during the pandemic following unprecedented reductions in worldwide human mobility. We here describe a SARS-CoV-2 lineage - designated B.1.620 - discovered in Lithuania and carrying many mutations and deletions in the spike protein shared with widespread variants of concern (VOCs), including E484K, S477N and deletions HV69Δ, Y144Δ, and LLA241/243Δ. As well as documenting the suite of mutations this lineage carries, we also describe its potential to be resistant to neutralising antibodies, accompanying travel histories for a subset of European cases, evidence of local B.1.620 transmission in Europe with a focus on Lithuania, and significance of its prevalence in Central Africa owing to recent genome sequencing efforts there. We make a case for its likely Central African origin using advanced phylogeographic inference methodologies incorporating recorded travel histories of infected travellers.

Untangling introductions and persistence in COVID-19 resurgence in Europe

After the first wave of SARS-CoV-2 infections in spring 2020, Europe experienced a resurgence of the virus starting in late summer 2020 that was deadlier and more difficult to contain1. Relaxed intervention measures and summer travel have been implicated as drivers of the second wave2. Here we build a phylogeographical model to evaluate how newly introduced lineages, as opposed to the rekindling of persistent lineages, contributed to the resurgence of COVID-19 in Europe. We inform this model using genomic, mobility and epidemiological data from 10 European countries and estimate that in many countries more than half of the lineages circulating in late summer resulted from new introductions since 15 June 2020. The success in onward transmission of newly introduced lineages was negatively associated with the local incidence of COVID-19 during this period. The pervasive spread of variants in summer 2020 highlights the threat of viral dissemination when restrictions are lifted, and this needs to be carefully considered in strategies to control the current spread of variants that are more transmissible and/or evade immunity. Our findings indicate that more effective and coordinated measures are required to contain the spread through cross-border travel even as vaccination is reducing disease burden.

SARS-CoV-2 European resurgence foretold: interplay of introductions and persistence by leveraging genomic and mobility data

Following the first wave of SARS-CoV-2 infections in spring 2020, Europe experienced a resurgence of the virus starting late summer that was deadlier and more difficult to contain. Relaxed intervention measures and summer travel have been implicated as drivers of the second wave. Here, we build a phylogeographic model to evaluate how newly introduced lineages, as opposed to the rekindling of persistent lineages, contributed to the COVID-19 resurgence in Europe. We inform this model using genomic, mobility and epidemiological data from 10 West European countries and estimate that in many countries more than 50% of the lineages circulating in late summer resulted from new introductions since June 15th. The success in onwards transmission of these lineages is predicted by SARS-CoV-2 incidence during this period. Relatively early introductions from Spain into the United Kingdom contributed to the successful spread of the 20A.EU1/B.1.177 variant. The pervasive spread of variants that have not been associated with an advantage in transmissibility highlights the threat of novel variants of concern that emerged more recently and have been disseminated by holiday travel. Our findings indicate that more effective and coordinated measures are required to contain spread through cross-border travel.