Phylodynamic analysis of foot-and-mouth disease virus evolution in Mar Chiquita, Argentina.

Foot-and-mouth disease is a highly contagious disease affecting cloven-hoofed animals, resulting in considerable economic losses. Its causal agent is foot-and-mouth disease virus (FMDV), a picornavirus. Due to its error-prone replication and rapid evolution, the transmission and evolutionary dynamics of FMDV can be studied using genomic epidemiological approaches. To analyze FMDV evolution and identify possible transmission routes in an Argentinean region, field samples that tested positive for FMDV by PCR were obtained from 21 farms located in the Mar Chiquita district. Whole FMDV genome sequences were obtained by PCR amplification in seven fragments and sequencing using the Sanger technique. The genome sequences obtained from these samples were then analyzed using phylogenetic, phylogeographic, and evolutionary approaches. Three local transmission clusters were detected among the sampled viruses. The dataset was analyzed using Bayesian phylodynamic methods with appropriate coalescent and relaxed molecular clock models. The estimated mean viral evolutionary rate was 1.17 × 10- 2 substitutions/site/year. No significant differences in the rate of viral evolution were observed between farms with vaccinated animals and those with unvaccinated animals. The most recent common ancestor of the sampled sequences was dated to approximately one month before the first reported case in the outbreak. Virus transmission started in the south of the district and later dispersed to the west, and finally arrived in the east. Different transmission routes among the studied herds, such as non-replicating vectors and close contact contagion (i.e., aerosols), may be responsible for viral spread.

Synchrony of Bird Migration with Global Dispersal of Avian Influenza Reveals Exposed Bird Orders.

Highly pathogenic avian influenza virus (HPAIV) A H5, particularly clade 2.3.4.4, has caused worldwide outbreaks in domestic poultry, occasional spillover to humans, and increasing deaths of diverse species of wild birds since 2014. Wild bird migration is currently acknowledged as an important ecological process contributing to the global dispersal of HPAIV H5. However, this mechanism has not been quantified using bird movement data from different species, and the timing and location of exposure of different species is unclear. We sought to explore these questions through phylodynamic analyses based on empirical data of bird movement tracking and virus genome sequences of clade 2.3.4.4 and 2.3.2.1. First, we demonstrate that seasonal bird migration can explain salient features of the global dispersal of clade 2.3.4.4. Second, we detect synchrony between the seasonality of bird annual cycle phases and virus lineage movements. We reveal the differing exposed bird orders at geographical origins and destinations of HPAIV H5 clade 2.3.4.4 lineage movements, including relatively under-discussed orders. Our study provides a phylodynamic framework that links the bird movement ecology and genomic epidemiology of avian influenza; it highlights the importance of integrating bird behavior and life history in avian influenza studies.

Characterization of the genomic sequence of a circo-like virus and of three chaphamaparvoviruses detected in mute swan (Cygnus olor).

We report the complete genomes of four ssDNA viruses: a circular replication-associated protein-encoding single-stranded DNA virus belonging to a clade previously detected only in mammals, and three chaphamaparvoviruses, which were detected by viromic surveillance of mute swan (Cygnus olor) fecal samples from the United Kingdom.

Ancient chicken remains reveal the origins of virulence in Marek's disease virus.

The pronounced growth in livestock populations since the 1950s has altered the epidemiological and evolutionary trajectory of their associated pathogens. For example, Marek's disease virus (MDV), which causes lymphoid tumors in chickens, has experienced a marked increase in virulence over the past century. Today, MDV infections kill >90% of unvaccinated birds, and controlling it costs more than US$1 billion annually. By sequencing MDV genomes derived from archeological chickens, we demonstrate that it has been circulating for at least 1000 years. We functionally tested the Meq oncogene, one of 49 viral genes positively selected in modern strains, demonstrating that ancient MDV was likely incapable of driving tumor formation. Our results demonstrate the power of ancient DNA approaches to trace the molecular basis of virulence in economically relevant pathogens.

SARS-CoV-2 evolution in the Omicron era.

Since SARS-CoV-2 BA.5 (Omicron) emerged and spread in 2022, Omicron lineages have markedly diversified. Here we review the evolutionary trajectories and processes that underpin the emergence of these lineages, and identify the most prevalent sublineages. We discuss the potential origins of second-generation BA.2 lineages. Simple and complex recombination, antigenic drift and convergent evolution have enabled SARS-CoV-2 to accumulate mutations that alter its antigenicity. We also discuss the potential evolutionary trajectories of SARS-CoV-2 in the future.

Convergent trends and spatiotemporal patterns of Aedes-borne arboviruses in Mexico and Central America.

BACKGROUND: Aedes-borne arboviruses cause both seasonal epidemics and emerging outbreaks with a significant impact on global health. These viruses share mosquito vector species, often infecting the same host population within overlapping geographic regions. Thus, comparative analyses of the virus evolutionary and epidemiological dynamics across spatial and temporal scales could reveal convergent trends. METHODOLOGY/PRINCIPAL FINDINGS: Focusing on Mexico as a case study, we generated novel chikungunya and dengue (CHIKV, DENV-1 and DENV-2) virus genomes from an epidemiological surveillance-derived historical sample collection, and analysed them together with longitudinally-collected genome and epidemiological data from the Americas. Aedes-borne arboviruses endemically circulating within the country were found to be introduced multiple times from lineages predominantly sampled from the Caribbean and Central America. For CHIKV, at least thirteen introductions were inferred over a year, with six of these leading to persistent transmission chains. For both DENV-1 and DENV-2, at least seven introductions were inferred over a decade. CONCLUSIONS/SIGNIFICANCE: Our results suggest that CHIKV, DENV-1 and DENV-2 in Mexico share evolutionary and epidemiological trajectories. The southwest region of the country was determined to be the most likely location for viral introductions from abroad, with a subsequent spread into the Pacific coast towards the north of Mexico. Virus diffusion patterns observed across the country are likely driven by multiple factors, including mobility linked to human migration from Central towards North America. Considering Mexico's geographic positioning displaying a high human mobility across borders, our results prompt the need to better understand the role of anthropogenic factors in the transmission dynamics of Aedes-borne arboviruses, particularly linked to land-based human migration.

Correction: Implementation of Genomic Surveillance of SARS-CoV-2 in the Caribbean: Lessons learned for sustainability in resource-limited settings.

[This corrects the article DOI: 10.1371/journal.pgph.0001455.].

Genetic diversity, recombination and cross-species transmission of a waterbird gammacoronavirus in the wild.

Viruses emerging from wildlife can cause outbreaks in humans and domesticated animals. Predicting the emergence of future pathogens and mitigating their impacts requires an understanding of what shapes virus diversity and dynamics in wildlife reservoirs. In order to better understand coronavirus ecology in wild species, we sampled birds within a coastal freshwater lagoon habitat across 5 years, focussing on a large population of mute swans (Cygnus olor) and the diverse species that they interact with. We discovered and characterised the full genome of a divergent gammacoronavirus belonging to the Goose coronavirus CB17 species. We investigated the genetic diversity and dynamics of this gammacoronavirus using untargeted metagenomic sequencing of 223 faecal samples from swans of known age and sex, and RT-PCR screening of 1632 additional bird samples. The virus circulated persistently within the bird community; virus prevalence in mute swans exhibited seasonal variations, but did not change with swan age-class or epidemiological year. One whole genome was fully characterised, and revealed that the virus originated from a recombination event involving an undescribed gammacoronavirus species. Multiple lineages of this gammacoronavirus co-circulated within our study population. Viruses from this species have recently been detected in aquatic birds from both the Anatidae and Rallidae families, implying that host species habitat sharing may be important in shaping virus host range. As the host range of the Goose coronavirus CB17 species is not limited to geese, we propose that this species name should be updated to 'Waterbird gammacoronavirus 1'. Non-invasive sampling of bird coronaviruses may provide a tractable model system for understanding the evolutionary and cross-species dynamics of coronaviruses.

Dispersal patterns and influence of air travel during the global expansion of SARS-CoV-2 variants of concern.

The Alpha, Beta, and Gamma SARS-CoV-2 variants of concern (VOCs) co-circulated globally during 2020 and 2021, fueling waves of infections. They were displaced by Delta during a third wave worldwide in 2021, which, in turn, was displaced by Omicron in late 2021. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of VOCs worldwide. We find that source-sink dynamics varied substantially by VOC and identify countries that acted as global and regional hubs of dissemination. We demonstrate the declining role of presumed origin countries of VOCs in their global dispersal, estimating that India contributed <15% of Delta exports and South Africa <1%-2% of Omicron dispersal. We estimate that >80 countries had received introductions of Omicron within 100 days of its emergence, associated with accelerated passenger air travel and higher transmissibility. Our study highlights the rapid dispersal of highly transmissible variants, with implications for genomic surveillance along the hierarchical airline network.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico.

Over 200 different SARS-CoV-2 lineages have been observed in Mexico by November 2021. To investigate lineage replacement dynamics, we applied a phylodynamic approach and explored the evolutionary trajectories of five dominant lineages that circulated during the first year of local transmission. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in Mexico. Lineages B.1.1.222 and B.1.1.519 exhibited similar dynamics, constituting clades that likely originated in Mexico and persisted for >12 months. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. For the largest B.1.617.2 clades, we further explored viral lineage movements across Mexico. Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

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.

Identification of Evolutionary Trajectories Shared across Human Betacoronaviruses.

Comparing the evolution of distantly related viruses can provide insights into common adaptive processes related to shared ecological niches. Phylogenetic approaches, coupled with other molecular evolution tools, can help identify mutations informative on adaptation, although the structural contextualization of these to functional sites of proteins may help gain insight into their biological properties. Two zoonotic betacoronaviruses capable of sustained human-to-human transmission have caused pandemics in recent times (SARS-CoV-1 and SARS-CoV-2), although a third virus (MERS-CoV) is responsible for sporadic outbreaks linked to animal infections. Moreover, two other betacoronaviruses have circulated endemically in humans for decades (HKU1 and OC43). To search for evidence of adaptive convergence between established and emerging betacoronaviruses capable of sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1, and SARS-CoV-2), we developed a methodological pipeline to classify shared nonsynonymous mutations as putatively denoting homoplasy (repeated mutations that do not share direct common ancestry) or stepwise evolution (sequential mutations leading towards a novel genotype). In parallel, we look for evidence of positive selection and draw upon protein structure data to identify potential biological implications. We find 30 candidate mutations, from which 4 (codon sites 18121 [nsp14/residue 28], 21623 [spike/21], 21635 [spike/25], and 23948 [spike/796]; SARS-CoV-2 genome numbering) further display evolution under positive selection and proximity to functional protein regions. Our findings shed light on potential mechanisms underlying betacoronavirus adaptation to the human host and pinpoint common mutational pathways that may occur during establishment of human endemicity.

Implementation of genomic surveillance of SARS-CoV-2 in the Caribbean: Lessons learned for sustainability in resource-limited settings.

The COVID-19 pandemic highlighted the importance of global genomic surveillance to monitor the emergence and spread of SARS-CoV-2 variants and inform public health decision-making. Until December 2020 there was minimal capacity for viral genomic surveillance in most Caribbean countries. To overcome this constraint, the COVID-19: Infectious disease Molecular epidemiology for PAthogen Control & Tracking (COVID-19 IMPACT) project was implemented to establish rapid SARS-CoV-2 whole genome nanopore sequencing at The University of the West Indies (UWI) in Trinidad and Tobago (T&T) and provide needed SARS-CoV-2 sequencing services for T&T and other Caribbean Public Health Agency Member States (CMS). Using the Oxford Nanopore Technologies MinION sequencing platform and ARTIC network sequencing protocols and bioinformatics pipeline, a total of 3610 SARS-CoV-2 positive RNA samples, received from 17 CMS, were sequenced in-situ during the period December 5th 2020 to December 31st 2021. Ninety-one Pango lineages, including those of five variants of concern (VOC), were identified. Genetic analysis revealed at least 260 introductions to the CMS from other global regions. For each of the 17 CMS, the percentage of reported COVID-19 cases sequenced by the COVID-19 IMPACT laboratory ranged from 0·02% to 3·80% (median = 1·12%). Sequences submitted to GISAID by our study represented 73·3% of all SARS-CoV-2 sequences from the 17 CMS available on the database up to December 31st 2021. Increased staffing, process and infrastructural improvement over the course of the project helped reduce turnaround times for reporting to originating institutions and sequence uploads to GISAID. Insights from our genomic surveillance network in the Caribbean region directly influenced non-pharmaceutical countermeasures in the CMS countries. However, limited availability of associated surveillance and clinical data made it challenging to contextualise the observed SARS-CoV-2 diversity and evolution, highlighting the need for development of infrastructure for collecting and integrating genomic sequencing data and sample-associated metadata.

Bidirectional Movement of Emerging H5N8 Avian Influenza Viruses Between Europe and Asia via Migratory Birds Since Early 2020.

Migratory birds play a critical role in the rapid spread of highly pathogenic avian influenza (HPAI) H5N8 virus clade 2.3.4.4 across Eurasia. Elucidating the timing and pattern of virus transmission is essential therefore for understanding the spatial dissemination of these viruses. In this study, we surveyed >27,000 wild birds in China, tracked the year-round migration patterns of 20 bird species across China since 2006, and generated new HPAI H5N8 virus genomic data. Using this new data set, we investigated the seasonal transmission dynamics of HPAI H5N8 viruses across Eurasia. We found that introductions of HPAI H5N8 viruses to different Eurasian regions were associated with the seasonal migration of wild birds. Moreover, we report a backflow of HPAI H5N8 virus lineages from Europe to Asia, suggesting that Europe acts as both a source and a sink in the global HPAI virus transmission network.

Seasonal dynamics of the wild rodent faecal virome.

Viral discovery studies in wild animals often rely on cross-sectional surveys at a single time point. As a result, our understanding of the temporal stability of wild animal viromes remains poorly resolved. While studies of single host-virus systems indicate that host and environmental factors influence seasonal virus transmission dynamics, comparable insights for whole viral communities in multiple hosts are lacking. Utilizing noninvasive faecal samples from a long-term wild rodent study, we characterized viral communities of three common European rodent species (Apodemus sylvaticus, A. flavicollis and Myodes glareolus) living in temperate woodland over a single year. Our findings indicate that a substantial fraction of the rodent virome is seasonally transient and associated with vertebrate or bacteria hosts. Further analyses of one of the most common virus families, Picornaviridae, show pronounced temporal changes in viral richness and evenness, which were associated with concurrent and up to ~3-month lags in host density, ambient temperature, rainfall and humidity, suggesting complex feedbacks from the host and environmental factors on virus transmission and shedding in seasonal habitats. Overall, this study emphasizes the importance of understanding the seasonal dynamics of wild animal viromes in order to better predict and mitigate zoonotic risks.

Impact of host age on viral and bacterial communities in a waterbird population.

Wildlife harbour pathogens that can harm human or livestock health and are the source of most emerging infectious diseases. It is rarely considered how changes in wildlife population age-structures or how age-stratified behaviours might alter the level of pathogen detection within a species, or risk of spillover to other species. Micro-organisms that occur in healthy animals can be an important model for understanding and predicting the dynamics of pathogens of greater health concern, which are hard to study in wild populations due to their relative rarity. We therefore used a metagenomic approach to jointly characterise viral and prokaryotic carriage in faeces collected from a healthy wild bird population (Cygnus olor; mute swan) that has been subject to long-term study. Using 223 samples from known individuals allowed us to compare differences in prokaryotic and eukaryotic viral carriage between adults and juveniles at an unprecedented level of detail. We discovered and characterised 77 novel virus species, of which 21% belong putatively to bird-infecting families, and described the core prokaryotic microbiome of C. olor. Whilst no difference in microbiota diversity was observed between juveniles and adult individuals, 50% (4/8) of bird-infecting virus families (picornaviruses, astroviruses, adenoviruses and bornaviruses) and 3.4% (9/267) of prokaryotic families (including Helicobacteraceae, Spirochaetaceae and Flavobacteriaceae families) were differentially abundant and/or prevalent between juveniles and adults. This indicates that perturbations that affect population age-structures of wildlife could alter circulation dynamics and spillover risk of microbes, potentially including pathogens.

Identification of evolutionary trajectories shared across human betacoronaviruses.

Comparing the evolution of distantly related viruses can provide insights into common adaptive processes related to shared ecological niches. Phylogenetic approaches, coupled with other molecular evolution tools, can help identify mutations informative on adaptation, whilst the structural contextualization of these to functional sites of proteins may help gain insight into their biological properties. Two zoonotic betacoronaviruses capable of sustained human-to-human transmission have caused pandemics in recent times (SARS-CoV-1 and SARS-CoV-2), whilst a third virus (MERS-CoV) is responsible for sporadic outbreaks linked to animal infections. Moreover, two other betacoronaviruses have circulated endemically in humans for decades (HKU1 and OC43). To search for evidence of adaptive convergence between established and emerging betacoronaviruses capable of sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1 and SARS-CoV-2), we developed a methodological pipeline to classify shared non-synonymous mutations as putatively denoting homoplasy (repeated mutations that do not share direct common ancestry) or stepwise evolution (sequential mutations leading towards a novel genotype). In parallel, we look for evidence of positive selection, and draw upon protein structure data to identify potential biological implications. We find 30 mutations, with four of these [codon sites 18121 (nsp14/residue 28), 21623 (spike/21), 21635 (spike/25) and 23948 (spike/796); SARS-CoV-2 genome numbering] displaying evolution under positive selection and proximity to functional protein regions. Our findings shed light on potential mechanisms underlying betacoronavirus adaptation to the human host and pinpoint common mutational pathways that may occur during establishment of human endemicity.

The origins and molecular evolution of SARS-CoV-2 lineage B.1.1.7 in the UK.

The first SARS-CoV-2 variant of concern (VOC) to be designated was lineage B.1.1.7, later labelled by the World Health Organization as Alpha. Originating in early autumn but discovered in December 2020, it spread rapidly and caused large waves of infections worldwide. The Alpha variant is notable for being defined by a long ancestral phylogenetic branch with an increased evolutionary rate, along which only two sequences have been sampled. Alpha genomes comprise a well-supported monophyletic clade within which the evolutionary rate is typical of SARS-CoV-2. The Alpha epidemic continued to grow despite the continued restrictions on social mixing across the UK and the imposition of new restrictions, in particular, the English national lockdown in November 2020. While these interventions succeeded in reducing the absolute number of cases, the impact of these non-pharmaceutical interventions was predominantly to drive the decline of the SARS-CoV-2 lineages that preceded Alpha. We investigate the only two sampled sequences that fall on the branch ancestral to Alpha. We find that one is likely to be a true intermediate sequence, providing information about the order of mutational events that led to Alpha. We explore alternate hypotheses that can explain how Alpha acquired a large number of mutations yet remained largely unobserved in a region of high genomic surveillance: an under-sampled geographical location, a non-human animal population, or a chronically infected individual. We conclude that the latter provides the best explanation of the observed behaviour and dynamics of the variant, although the individual need not be immunocompromised, as persistently infected immunocompetent hosts also display a higher within-host rate of evolution. Finally, we compare the ancestral branches and mutation profiles of other VOCs and find that Delta appears to be an outlier both in terms of the genomic locations of its defining mutations and a lack of the rapid evolutionary rate on its ancestral branch. As new variants, such as Omicron, continue to evolve (potentially through similar mechanisms), it remains important to investigate the origins of other variants to identify ways to potentially disrupt their evolution and emergence.

Global Expansion of SARS-CoV-2 Variants of Concern: Dispersal Patterns and Influence of Air Travel.

In many regions of the world, the Alpha, Beta and Gamma SARS-CoV-2 Variants of Concern (VOCs) co-circulated during 2020-21 and fueled waves of infections. During 2021, these variants were almost completely displaced by the Delta variant, causing a third wave of infections worldwide. This phenomenon of global viral lineage displacement was observed again in late 2021, when the Omicron variant disseminated globally. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of SARS-CoV-2 VOCs worldwide. We find that the source-sink dynamics of SARS-CoV-2 varied substantially by VOC, and identify countries that acted as global hubs of variant dissemination, while other countries became regional contributors to the export of specific variants. We demonstrate a declining role of presumed origin countries of VOCs to their global dispersal: we estimate that India contributed <15% of all global exports of Delta to other countries and South Africa <1-2% of all global Omicron exports globally. We further estimate that >80 countries had received introductions of Omicron BA.1 100 days after its inferred date of emergence, compared to just over 25 countries for the Alpha variant. This increased speed of global dissemination was associated with a rebound in air travel volume prior to Omicron emergence in addition to the higher transmissibility of Omicron relative to Alpha. Our study highlights the importance of global and regional hubs in VOC dispersal, and the speed at which highly transmissible variants disseminate through these hubs, even before their detection and characterization through genomic surveillance. Highlights: Global phylogenetic analysis reveals relationship between air travel and speed of dispersal of SARS-CoV-2 variants of concern (VOCs)Omicron VOC spread to 5x more countries within 100 days of its emergence compared to all other VOCsOnward transmission and dissemination of VOCs Delta and Omicron was primarily from secondary hubs rather than initial country of detection during a time of increased global air travelAnalysis highlights highly connected countries identified as major global and regional exporters of VOCs.

Global disparities in SARS-CoV-2 genomic surveillance.

Genomic sequencing is essential to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments, vaccines, and guide public health responses. To investigate the global SARS-CoV-2 genomic surveillance, we used sequences shared via GISAID to estimate the impact of sequencing intensity and turnaround times on variant detection in 189 countries. In the first two years of the pandemic, 78% of high-income countries sequenced >0.5% of their COVID-19 cases, while 42% of low- and middle-income countries reached that mark. Around 25% of the genomes from high income countries were submitted within 21 days, a pattern observed in 5% of the genomes from low- and middle-income countries. We found that sequencing around 0.5% of the cases, with a turnaround time <21 days, could provide a benchmark for SARS-CoV-2 genomic surveillance. Socioeconomic inequalities undermine the global pandemic preparedness, and efforts must be made to support low- and middle-income countries improve their local sequencing capacity.