Early introductions and transmission of SARS-CoV-2 variant B.1.1.7 in the United States.

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2,500 COVID-19 cases associated with this variant have been detected in the United States (US) since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight that the primary ports of entry for B.1.1.7 in the US were in New York, California, and Florida. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid- to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

Projected decline in European bumblebee populations in the twenty-first century.

Habitat degradation and climate change are globally acting as pivotal drivers of wildlife collapse, with mounting evidence that this erosion of biodiversity will accelerate in the following decades1-3. Here, we quantify the past, present and future ecological suitability of Europe for bumblebees, a threatened group of pollinators ranked among the highest contributors to crop production value in the northern hemisphere4-8. We demonstrate coherent declines of bumblebee populations since 1900 over most of Europe and identify future large-scale range contractions and species extirpations under all future climate and land use change scenarios. Around 38-76% of studied European bumblebee species currently classified as 'Least Concern' are projected to undergo losses of at least 30% of ecologically suitable territory by 2061-2080 compared to 2000-2014. All scenarios highlight that parts of Scandinavia will become potential refugia for European bumblebees; it is however uncertain whether these areas will remain clear of additional anthropogenic stressors not accounted for in present models. Our results underline the critical role of global change mitigation policies as effective levers to protect bumblebees from manmade transformation of the biosphere.

Considerable escape of SARS-CoV-2 Omicron to antibody neutralization.

The SARS-CoV-2 Omicron variant was first identified in November 2021 in Botswana and South Africa1-3. It has since spread to many countries and is expected to rapidly become dominant worldwide. The lineage is characterized by the presence of around 32 mutations in spike-located mostly in the N-terminal domain and the receptor-binding domain-that may enhance viral fitness and enable antibody evasion. Here we isolated an infectious Omicron virus in Belgium from a traveller returning from Egypt. We examined its sensitivity to nine monoclonal antibodies that have been clinically approved or are in development4, and to antibodies present in 115 serum samples from COVID-19 vaccine recipients or individuals who have recovered from COVID-19. Omicron was completely or partially resistant to neutralization by all monoclonal antibodies tested. Sera from recipients of the Pfizer or AstraZeneca vaccine, sampled five months after complete vaccination, barely inhibited Omicron. Sera from COVID-19-convalescent patients collected 6 or 12 months after symptoms displayed low or no neutralizing activity against Omicron. Administration of a booster Pfizer dose as well as vaccination of previously infected individuals generated an anti-Omicron neutralizing response, with titres 6-fold to 23-fold lower against Omicron compared with those against Delta. Thus, Omicron escapes most therapeutic monoclonal antibodies and, to a large extent, vaccine-elicited antibodies. However, Omicron is neutralized by antibodies generated by a booster vaccine dose.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant.

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

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.

Phylogeographic reconstruction of the emergence and spread of Powassan virus in the northeastern United States.

Powassan virus is an emerging tick-borne virus of concern for public health, but very little is known about its transmission patterns and ecology. Here, we expanded the genomic dataset by sequencing 279 Powassan viruses isolated from Ixodes scapularis ticks from the northeastern United States. Our phylogeographic reconstructions revealed that Powassan virus lineage II was likely introduced or emerged from a relict population in the Northeast between 1940 and 1975. Sequences strongly clustered by sampling location, suggesting a highly focal geographical distribution. Our analyses further indicated that Powassan virus lineage II emerged in the northeastern United States mostly following a south-to-north pattern, with a weighted lineage dispersal velocity of ~3 km/y. Since the emergence in the Northeast, we found an overall increase in the effective population size of Powassan virus lineage II, but with growth stagnating during recent years. The cascading effect of population expansion of white-tailed deer and I. scapularis populations likely facilitated the emergence of Powassan virus in the northeastern United States.

Variant-specific introduction and dispersal dynamics of SARS-CoV-2 in New York City - from Alpha to Omicron.

Since the latter part of 2020, SARS-CoV-2 evolution has been characterised by the emergence of viral variants associated with distinct biological characteristics. While the main research focus has centred on the ability of new variants to increase in frequency and impact the effective reproductive number of the virus, less attention has been placed on their relative ability to establish transmission chains and to spread through a geographic area. Here, we describe a phylogeographic approach to estimate and compare the introduction and dispersal dynamics of the main SARS-CoV-2 variants - Alpha, Iota, Delta, and Omicron - that circulated in the New York City area between 2020 and 2022. Notably, our results indicate that Delta had a lower ability to establish sustained transmission chains in the NYC area and that Omicron (BA.1) was the variant fastest to disseminate across the study area. The analytical approach presented here complements non-spatially-explicit analytical approaches that seek a better understanding of the epidemiological differences that exist among successive SARS-CoV-2 variants of concern.

Early Genomic Surveillance and Phylogeographic Analysis of Getah Virus, a Reemerging Arbovirus, in Livestock in China.

Getah virus (GETV) mainly causes disease in livestock and may pose an epidemic risk due to its expanding host range and the potential of long-distance dispersal through animal trade. Here, we used metagenomic next-generation sequencing (mNGS) to identify GETV as the pathogen responsible for reemerging swine disease in China and subsequently estimated key epidemiological parameters using phylodynamic and spatially-explicit phylogeographic approaches. The GETV isolates were able to replicate in a variety of cell lines, including human cells, and showed high pathogenicity in a mouse model, suggesting the potential for more mammal hosts. We obtained 16 complete genomes and 79 E2 gene sequences from viral strains collected in China from 2016 to 2021 through large-scale surveillance among livestock, pets, and mosquitoes. Our phylogenetic analysis revealed that three major GETV lineages are responsible for the current epidemic in livestock in China. We identified three potential positively selected sites and mutations of interest in E2, which may impact the transmissibility and pathogenicity of the virus. Phylodynamic inference of the GETV demographic dynamics identified an association between livestock meat consumption and the evolution of viral genetic diversity. Finally, phylogeographic reconstruction of GETV dispersal indicated that the sampled lineages have preferentially circulated within areas associated with relatively higher mean annual temperature and pig population density. Our results highlight the importance of continuous surveillance of GETV among livestock in southern Chinese regions associated with relatively high temperatures. IMPORTANCE Although livestock is known to be the primary reservoir of Getah virus (GETV) in Asian countries, where identification is largely based on serology, the evolutionary history and spatial epidemiology of GETV in these regions remain largely unknown. Through our sequencing efforts, we provided robust support for lineage delineation of GETV and identified three major lineages that are responsible for the current epidemic in livestock in China. We further analyzed genomic and epidemiological data to reconstruct the recent demographic and dispersal history of GETV in domestic animals in China and to explore the impact of environmental factors on its genetic diversity and its diffusion. Notably, except for livestock meat consumption, other pig-related factors such as the evolution of live pig transport and pork production do not show a significant association with the evolution of viral genetic diversity, pointing out that further studies should investigate the potential contribution of other host species to the GETV outbreak. Our analysis of GETV demonstrates the need for wider animal species surveillance and provides a baseline for future studies of the molecular epidemiology and early warning of emerging arboviruses in China.

Complete genome reconstruction of the global and European regional dispersal history of the lumpy skin disease virus.

IMPORTANCE: Lumpy skin disease virus (LSDV) has a complex epidemiology involving multiple strains, recombination, and vaccination. Its DNA genome provides limited genetic variation to trace outbreaks in space and time. Sequencing of LSDV whole genomes has also been patchy at global and regional scales. Here, we provide the first fine-grained whole genome sequence sampling of a constrained LSDV outbreak (southeastern Europe, 2015-2017), which we analyze along with global publicly available genomes. We formally evaluate the past occurrence of recombination events as well as the temporal signal that is required for calibrating molecular clock models and subsequently conduct a time-calibrated spatially explicit phylogeographic reconstruction. Our study further illustrates the importance of accounting for recombination events before reconstructing global and regional dynamics of DNA viruses. More LSDV whole genomes from endemic areas are needed to obtain a comprehensive understanding of global LSDV dispersal dynamics.

Dispersal history of SARS-CoV-2 in Galicia, Spain.

The dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron-BA.1 variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the United States became increasingly significant. The number of detected introductions varied from 96 and 101 for Alpha and Delta to 39 for Omicron-BA.1. Most of these introductions left a low number of descendants (<10), suggesting a limited impact on the evolution of the pandemic in Galicia. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.

Combined Phylogeographic Analyses and Epidemiologic Contact Tracing to Characterize Atypically Pathogenic Avian Influenza (H3N1) Epidemic, Belgium, 2019.

The high economic impact and zoonotic potential of avian influenza call for detailed investigations of dispersal dynamics of epidemics. We integrated phylogeographic and epidemiologic analyses to investigate the dynamics of a low pathogenicity avian influenza (H3N1) epidemic that occurred in Belgium during 2019. Virus genomes from 104 clinical samples originating from 85% of affected farms were sequenced. A spatially explicit phylogeographic analysis confirmed a dominating northeast to southwest dispersal direction and a long-distance dispersal event linked to direct live animal transportation between farms. Spatiotemporal clustering, transport, and social contacts strongly correlated with the phylogeographic pattern of the epidemic. We detected only a limited association between wind direction and direction of viral lineage dispersal. Our results highlight the multifactorial nature of avian influenza epidemics and illustrate the use of genomic analyses of virus dispersal to complement epidemiologic and environmental data, improve knowledge of avian influenza epidemiologic dynamics, and enhance control strategies.

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.

Dispersal dynamics of SARS-CoV-2 lineages during the first epidemic wave in New York City.

During the first phase of the COVID-19 epidemic, New York City rapidly became the epicenter of the pandemic in the United States. While molecular phylogenetic analyses have previously highlighted multiple introductions and a period of cryptic community transmission within New York City, little is known about the circulation of SARS-CoV-2 within and among its boroughs. We here perform phylogeographic investigations to gain insights into the circulation of viral lineages during the first months of the New York City outbreak. Our analyses describe the dispersal dynamics of viral lineages at the state and city levels, illustrating that peripheral samples likely correspond to distinct dispersal events originating from the main metropolitan city areas. In line with the high prevalence recorded in this area, our results highlight the relatively important role of the borough of Queens as a transmission hub associated with higher local circulation and dispersal of viral lineages toward the surrounding boroughs.

Uncovering the endemic circulation of rabies in Cambodia.

In epidemiology, endemicity characterizes sustained pathogen circulation in a geographical area, which involves a circulation that is not being maintained by external introductions. Because it could potentially shape the design of public health interventions, there is an interest in fully uncovering the endemic pattern of a disease. Here, we use a phylogeographic approach to investigate the endemic signature of rabies virus (RABV) circulation in Cambodia. Cambodia is located in one of the most affected regions by rabies in the world, but RABV circulation between and within Southeast Asian countries remains understudied. Our analyses are based on a new comprehensive data set of 199 RABV genomes collected between 2014 and 2017 as well as previously published Southeast Asian RABV sequences. We show that most Cambodian sequences belong to a distinct clade that has been circulating almost exclusively in Cambodia. Our results thus point towards rabies circulation in Cambodia that does not rely on external introductions. We further characterize within-Cambodia RABV circulation by estimating lineage dispersal metrics that appear to be similar to other settings, and by performing landscape phylogeographic analyses to investigate environmental factors impacting the dispersal dynamic of viral lineages. The latter analyses do not lead to the identification of environmental variables that would be associated with the heterogeneity of viral lineage dispersal velocities, which calls for a better understanding of local dog ecology and further investigations of the potential drivers of RABV spread in the region. Overall, our study illustrates how phylogeographic investigations can be performed to assess and characterize viral endemicity in a context of relatively limited data.

Virus genomes reveal factors that spread and sustained the Ebola epidemic.

The 2013-2016 West African epidemic caused by the Ebola virus was of unprecedented magnitude, duration and impact. Here we reconstruct the dispersal, proliferation and decline of Ebola virus throughout the region by analysing 1,610 Ebola virus genomes, which represent over 5% of the known cases. We test the association of geography, climate and demography with viral movement among administrative regions, inferring a classic 'gravity' model, with intense dispersal between larger and closer populations. Despite attenuation of international dispersal after border closures, cross-border transmission had already sown the seeds for an international epidemic, rendering these measures ineffective at curbing the epidemic. We address why the epidemic did not spread into neighbouring countries, showing that these countries were susceptible to substantial outbreaks but at lower risk of introductions. Finally, we reveal that this large epidemic was a heterogeneous and spatially dissociated collection of transmission clusters of varying size, duration and connectivity. These insights will help to inform interventions in future epidemics.

Investigation of an international water polo tournament in Czechia as a potential source for early introduction of the SARS-CoV-2 Omicron variant into Belgium, Switzerland and Germany, November 2021.

BackgroundThe earliest recognised infections by the SARS-CoV-2 Omicron variant (Pango lineage B.1.1.529) in Belgium and Switzerland suggested a connection to an international water polo tournament, held 12-14 November 2021 in Brno, Czechia.AimTo study the arrival and subsequent spread of the Omicron variant in Belgium and Switzerland, and understand the overall importance of this international sporting event on the number of infections in the two countries.MethodsWe performed intensive forward and backward contact tracing in both countries, supplemented by phylogenetic investigations using virus sequences of the suspected infection chain archived in public databases.ResultsThrough contact tracing, we identified two and one infected athletes of the Belgian and Swiss water polo teams, respectively, and subsequently also three athletes from Germany. In Belgium and Switzerland, four and three secondary infections, and three and one confirmed tertiary infections were identified. Phylogenetic investigation demonstrated that this sporting event played a role as the source of infection, but without a direct link with infections from South Africa and not as a superspreading event; the virus was found to already be circulating at that time in the countries involved.ConclusionThe SARS-CoV-2 Omicron variant started to circulate in Europe several weeks before its identification in South Africa on 24 November 2021. Accordingly, it can be assumed that travel restrictions are usually implemented too late to prevent the spread of newly detected SARS-CoV-2 variants to other regions. Phylogenetic analysis may modify the perception of an apparently clear result of intensive contact tracing.

Relax, Keep Walking - A Practical Guide to Continuous Phylogeographic Inference with BEAST.

Spatially explicit phylogeographic analyses can be performed with an inference framework that employs relaxed random walks to reconstruct phylogenetic dispersal histories in continuous space. This core model was first implemented 10 years ago and has opened up new opportunities in the field of phylodynamics, allowing researchers to map and analyze the spatial dissemination of rapidly evolving pathogens. We here provide a detailed and step-by-step guide on how to set up, run, and interpret continuous phylogeographic analyses using the programs BEAUti, BEAST, Tracer, and TreeAnnotator.

A Phylodynamic Workflow to Rapidly Gain Insights into the Dispersal History and Dynamics of SARS-CoV-2 Lineages.

Since the start of the COVID-19 pandemic, an unprecedented number of genomic sequences of SARS-CoV-2 have been generated and shared with the scientific community. The unparalleled volume of available genetic data presents a unique opportunity to gain real-time insights into the virus transmission during the pandemic, but also a daunting computational hurdle if analyzed with gold-standard phylogeographic approaches. To tackle this practical limitation, we here describe and apply a rapid analytical pipeline to analyze the spatiotemporal dispersal history and dynamics of SARS-CoV-2 lineages. As a proof of concept, we focus on the Belgian epidemic, which has had one of the highest spatial densities of available SARS-CoV-2 genomes. Our pipeline has the potential to be quickly applied to other countries or regions, with key benefits in complementing epidemiological analyses in assessing the impact of intervention measures or their progressive easement.

Genomic Surveillance of Yellow Fever Virus Epizootic in São Paulo, Brazil, 2016 - 2018.

São Paulo, a densely inhabited state in southeast Brazil that contains the fourth most populated city in the world, recently experienced its largest yellow fever virus (YFV) outbreak in decades. YFV does not normally circulate extensively in São Paulo, so most people were unvaccinated when the outbreak began. Surveillance in non-human primates (NHPs) is important for determining the magnitude and geographic extent of an epizootic, thereby helping to evaluate the risk of YFV spillover to humans. Data from infected NHPs can give more accurate insights into YFV spread than when using data from human cases alone. To contextualise human cases, identify epizootic foci and uncover the rate and direction of YFV spread in São Paulo, we generated and analysed virus genomic data and epizootic case data from NHPs in São Paulo. We report the occurrence of three spatiotemporally distinct phases of the outbreak in São Paulo prior to February 2018. We generated 51 new virus genomes from YFV positive cases identified in 23 different municipalities in São Paulo, mostly sampled from NHPs between October 2016 and January 2018. Although we observe substantial heterogeneity in lineage dispersal velocities between phylogenetic branches, continuous phylogeographic analyses of generated YFV genomes suggest that YFV lineages spread in São Paulo at a mean rate of approximately 1km per day during all phases of the outbreak. Viral lineages from the first epizootic phase in northern São Paulo subsequently dispersed towards the south of the state to cause the second and third epizootic phases there. This alters our understanding of how YFV was introduced into the densely populated south of São Paulo state. Our results shed light on the sylvatic transmission of YFV in highly fragmented forested regions in São Paulo state and highlight the importance of continued surveillance of zoonotic pathogens in sentinel species.

Genome-associations of extended-spectrum ß-lactamase producing (ESBL) or AmpC producing E. coli in small and medium pig farms from Khon Kaen province, Thailand.

Thailand is undergoing rapid intensification of livestock production where small subsistence farms and medium sized commercial farms coexist. In medium farms, antimicrobials are prescribed by a veterinarian, whereas in small farms antimicrobial use remains largely unsupervised. The impact of these differences as well as other farming practices on the emergence and composition of antimicrobial resistance genes (ARGs) remains largely unknown. We analyzed 363 genomes of extended-spectrum ß-lactamase producing (ESBL) and/or AmpC producing Escherichia coli recovered from humans and pigs at small and medium farms from the Khon Kaen province, Thailand. We tested for genome-wide associations to identify links between ARGs, host, and farm size. Pig isolates from small farms were associated with mcr and qnr genes conferring resistance to colistin and fluoroquinolones, respectively. In contrast, pig isolates from medium farms were associated with ARGs conferring resistance to drugs commonly used on medium farms (i.e., streptomycin). ESBL plasmids from small farms co-carried ARGs conferring resistance to critically important antimicrobials more frequently compared to plasmid from medium farms. Frequent ARG combinations included blaCTX-M-55 + qnrS1 (29.8% vs 17.5% in small and medium farms, respectively), blaCTX-M-55 + qnrS1 + mcr-3.19 (5% vs 0%), blaCTX-M-14 + qnrS1 (9.3% vs 6.2%), and blaCTX-M-14 + qnrS1 + mcr-1.1 (3.1% vs 0%). The co-location on plasmids of ARGs conferring resistance to critically important antimicrobials as defined by the World Health Organization is concerning, and actions to curb their spread are urgently needed. Legislation on limiting antimicrobial sales and initiatives to better inform farmers and veterinarians on appropriate antimicrobial usage and farm biosecurity could help reduce antimicrobial use on farms.