Using phylogeographic approaches to analyse the dispersal history, velocity and direction of viral lineages - Application to rabies virus spread in Iran

Typing Arctic and Africa-2 clades of rabies virus using clade-defining single nucleotide polymorphisms

Rabies is a deadly neurotropic, zoonotic disease with a mortality rate of 100% after symptoms appear. Rabies virus (RABV) is the primary cause of rabies disease in humans, and it mainly spreads via dog bites in developing countries. Over the course of RABV evolution, multiple RABV variants, called clades, have emerged. However, our understanding of these clades is limited, as the only method to identify a clade is sequencing, followed by phylogeny. In this study, we have developed a rapid, PCR-based method for typing two RABV clades. We utilised highly conserved amino acid changes specific to the Arctic and Africa-2 clades of the rabies virus (RABV). A single nucleotide substitution from adenine to thymine at position 178 within the nucleoprotein gene was found to be clade-specific in the Arctic clade. Similarly, adenine at position 638 is a distinctive marker for the Africa-2 clade. The assay demonstrated high specificity and offers the added benefit of PCR-based amplification, enabling virus detection even when viral titers are low. The assay can identify the Arctic clade and Africa-2 clade without sequencing and is highly specific and sensitive. Furthermore, this method can be adapted to detect other RABV clades and a wide range of viruses.

Combining genomics and epidemiology to investigate a zoonotic outbreak of rabies in Romblon Province, Philippines

Rabies is a viral zoonosis that kills thousands of people annually in low- and middle-income countries across Africa and Asia where domestic dogs are the reservoir. 'Zero by 30', the global strategy to end dog-mediated human rabies, promotes a One Health approach underpinned by mass dog vaccination, post-exposure vaccination of bite victims, robust surveillance and community engagement. Using Integrated Bite Case Management (IBCM) and whole genome sequencing (WGS), we enhanced rabies surveillance to detect an outbreak in a formerly rabies-free island province in the Philippines. We inferred that the outbreak was seeded by at least three independent human-mediated introductions that were identified as coming from neighbouring rabies-endemic provinces. Considerable local transmission went undetected, and two human deaths occurred within 6 months of outbreak detection. Suspension of routine dog vaccination due to COVID-19 restrictions likely facilitated rabies spread from these introductions. Emergency response, consisting of awareness measures, and ring vaccination, were performed, but swifter and more widespread implementation is needed to contain and eliminate the outbreak and to secure rabies freedom. We conclude that strengthened surveillance making use of new tools such as IBCM, WGS, and rapid diagnostic tests can support One Health in action and progress towards the 'Zero by 30' goal.

How fast are viruses spreading in the wild?

Genomic data collected from viral outbreaks can be exploited to reconstruct the dispersal history of viral lineages in a two-dimensional space using continuous phylogeographic inference. These spatially explicit reconstructions can subsequently be used to estimate dispersal metrics that can be informative of the dispersal dynamics and the capacity to spread among hosts. Heterogeneous sampling efforts of genomic sequences can however impact the accuracy of phylogeographic dispersal metrics. While the impact of spatial sampling bias on the outcomes of continuous phylogeographic inference has previously been explored, the impact of sampling intensity (i.e., sampling size) when aiming to characterise dispersal patterns through continuous phylogeographic reconstructions has not yet been thoroughly evaluated. In our study, we use simulations to evaluate the robustness of 3 dispersal metrics - a lineage dispersal velocity, a diffusion coefficient, and an isolation-by-distance (IBD) signal metric - to the sampling intensity. Our results reveal that both the diffusion coefficient and IBD signal metrics appear to be the most robust to the number of samples considered for the phylogeographic reconstruction. We then use these 2 dispersal metrics to compare the dispersal pattern and capacity of various viruses spreading in animal populations. Our comparative analysis reveals a broad range of IBD patterns and diffusion coefficients mostly reflecting the dispersal capacity of the main infected host species but also, in some cases, the likely signature of rapid and/or long-distance dispersal events driven by human-mediated movements through animal trade. Overall, our study provides key recommendations for the use of lineage dispersal metrics to consider in future studies and illustrates their application to compare the spread of viruses in various settings.

Identification of clade-defining single nucleotide polymorphisms for improved rabies virus surveillance

Background

Rabies is an ancient disease that remains endemic in many countries. It causes many human deaths annually, predominantly in resource-poor countries. Over evolutionary timelines, several rabies virus (RABV) genotypes have stabilised, forming distinct clades. Extensive studies have been conducted on the origin, occurrence and spread of RABV clades. Single nucleotide polymorphisms (SNPs) distribution across the RABV genome and its clades remains largely unknown, highlighting the need for comprehensive whole-genome analyses.

Methods

We accessed whole genome sequences for RABV from public databases and identified SNPs across the whole genome sequences. Then, we annotated these SNPs using an R script, and these SNPs were categorised into different categories; universal, clade-specific, and clade-defining, based on the frequency of occurrence.

Results

In this study, we present the SNPs occurring in the RABV based on whole genome sequences belonging to 8 clades isolated from 7 different host species likely to harbour dog-related rabies. We classified mutations into several classes based on their location within the genome and assessed the effect of SNP mutations on the viral glycoprotein.

Conclusions

The clade-defining mutations have implications for targeted surveillance and classification of clades. Additionally, we investigated the effects of these mutations on the Glycoprotein of the virus. Our findings contribute to expanding knowledge about RABV clade diversity and evolution, which has significant implications for effectively tracking and combatting RABV transmission.

Dynamics and ecology of a multi-stage expansion of Oropouche virus in Brazil

In March 2024, the Pan American Health Organization (PAHO) issued an alert in response to a rapid increase in Oropouche fever cases across South America. Brazil has been particularly affected, reporting a novel reassortant lineage of the Oropouche virus (OROV) and expansion to previously non-endemic areas beyond the Amazon Basin. Utilising phylogeographic approaches, we reveal a multi-scale expansion process with both short and long-distance dispersal events, and diffusion velocities in line with human-mediated jumps. We identify forest cover, banana and cocoa cultivation, temperature, and human population density as key environmental factors associated with OROV range expansion. Using ecological niche modelling, we show that OROV circulated in areas of enhanced ecological suitability immediately preceding its explosive epidemic expansion in the Amazon. This likely resulted from the virus being introduced into simultaneously densely populated and environmentally favourable regions in the Amazon, such as Manaus, leading to an amplified epidemic and spread beyond the Amazon. Our study provides valuable insights into the dispersal and ecological dynamics of OROV, highlighting the role of human mobility in colonisation of new areas, and raising concern over high viral suitability along the Brazilian coast.

Interrogating Genomes and Geography to Unravel Multiyear Vesicular Stomatitis Epizootics

We conducted an integrative analysis to elucidate the spatial epidemiological patterns of the Vesicular Stomatitis New Jersey virus (VSNJV) during the 2014-15 epizootic cycle in the United States (US). Using georeferenced VSNJV genomics data, confirmed vesicular stomatitis (VS) disease cases from surveillance, and a suite of environmental factors, our study assessed environmental and phylogenetic similarity to compare VS cases reported in 2014 and 2015. Despite uncertainties from incomplete virus sampling and cross-scale spatial processes, patterns suggested multiple independent re-invasion events concurrent with potential viral overwintering between sequential seasons. Our findings pointed to a geographically defined southern virus pool at the US-Mexico interface as the source of VSNJV invasions and overwintering sites. Phylodynamic analysis demonstrated an increase in virus diversity before a rise in case numbers and a pronounced reduction in virus diversity during the winter season, indicative of a genetic bottleneck and a significant narrowing of virus variation between the summer outbreak seasons. Environment-vector interactions underscored the central role of meta-population dynamics in driving disease spread. These insights emphasize the necessity for location- and time-specific management practices, including rapid response, movement restrictions, vector control, and other targeted interventions.

How fast are viruses spreading in the wild?

Genomic data collected from viral outbreaks can be exploited to reconstruct the dispersal history of viral lineages in a two-dimensional space using continuous phylogeographic inference. These spatially explicit reconstructions can subsequently be used to estimate dispersal metrics allowing to unveil the dispersal dynamics and evaluate the capacity to spread among hosts. Heterogeneous sampling intensity of genomic sequences can however impact the accuracy of dispersal insights gained through phylogeographic inference. In our study, we implement a simulation framework to evaluate the robustness of three dispersal metrics - a lineage dispersal velocity, a diffusion coefficient, and an isolation-by-distance signal metric - to the sampling effort. Our results reveal that both the diffusion coefficient and isolation-by-distance signal metrics appear to be robust to the number of samples considered for the phylogeographic reconstruction. We then use these two dispersal metrics to compare the dispersal pattern and capacity of various viruses spreading in animal populations. Our comparative analysis reveals a broad range of isolation-by-distance patterns and diffusion coefficients mostly reflecting the dispersal capacity of the main infected host species but also, in some cases, the likely signature of rapid and/or long-distance dispersal events driven by human-mediated movements through animal trade. Overall, our study provides key recommendations for the lineage dispersal metrics to consider in future studies and illustrates their application to compare the spread of viruses in various settings.

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.

Genomic Surveillance of Rabies Virus in Georgian Canines

Rabies is a fatal zoonosis that is considered a re-emerging infectious disease. Although rabies remains endemic in canines throughout much of the world, vaccination programs have essentially eliminated dog rabies in the Americas and much of Europe. However, despite the goal of eliminating dog rabies in the European Union by 2020, sporadic cases of dog rabies still occur in Eastern Europe, including Georgia. To assess the genetic diversity of the strains recently circulating in Georgia, we sequenced seventy-eight RABV-positive samples from the brain tissues of rabid dogs and jackals using Illumina short-read sequencing of total RNA shotgun libraries. Seventy-seven RABV genomes were successfully assembled and annotated, with seventy-four of them reaching the coding-complete status. Phylogenetic analyses of the nucleoprotein (N) and attachment glycoprotein (G) genes placed all the assembled genomes into the Cosmopolitan clade, consistent with the Georgian origin of the samples. An amino acid alignment of the G glycoprotein ectodomain identified twelve different sequences for this domain among the samples. Only one of the ectodomain groups contained a residue change in an antigenic site, an R264H change in the G5 antigenic site. Three isolates were cultured, and these were found to be efficiently neutralized by the human monoclonal antibody A6. Overall, our data show that recently circulating RABV isolates from Georgian canines are predominantly closely related phylogroup I viruses of the Cosmopolitan clade. Current human rabies vaccines should offer protection against infection by Georgian canine RABVs. The genomes have been deposited in GenBank (accessions: OQ603609-OQ603685).

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.

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.

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.

The Comparison of Full G and N Gene Sequences From Turkish Rabies Virus Field Strains

The rabies infection is a zoonotic viral disease in humans and is spread by both wild and domestic carnivores. This study aimed to molecularly characterize the field strains of the rabies virus circulating in Turkey between 2013 and 2020. Brain samples obtained from 16 infected animals (8 cattle, one donkey, three foxes, three dogs, and one marten) were tested. Full nucleoprotein (N) and glycoprotein (G) gene sequences were used to determine the genetic and antigenic characteristics of the rabies virus field strains. The phylogenetic analyses revealed that the 16 field strains identified in Turkey belonged to the Cosmopolitan lineage.

The Importance of Accurate Host Species Identification in the Framework of Rabies Surveillance, Control and Elimination

Accurate host identification is paramount to understand disease epidemiology and to apply appropriate control measures. This is especially important for multi-host pathogens such as the rabies virus, a major and almost invariably fatal zoonosis that has mobilized unanimous engagement at an international level towards the final goal of zero human deaths due to canine rabies. Currently, diagnostic laboratories implement a standardized identification using taxonomic keys. However, this method is challenged by high and undiscovered biodiversity, decomposition of carcasses and subjective misevaluation, as has been attested to by findings from a cohort of 242 archived specimens collected across Sub-Saharan Africa and submitted for rabies diagnosis. We applied two simple and cheap methods targeting the Cytochrome b and Cytochrome c oxidase subunit I to confirm the initial classification. We therefore suggest prioritizing a standardized protocol that includes, as a first step, the implementation of taxonomic keys at a family or subfamily level, followed by the molecular characterization of the host species.

On the Use of Phylogeographic Inference to Infer the Dispersal History of Rabies Virus: A Review Study

Rabies is a neglected zoonotic disease which is caused by negative strand RNA-viruses belonging to the genus Lyssavirus. Within this genus, rabies viruses circulate in a diverse set of mammalian reservoir hosts, is present worldwide, and is almost always fatal in non-vaccinated humans. Approximately 59,000 people are still estimated to die from rabies each year, leading to a global initiative to work towards the goal of zero human deaths from dog-mediated rabies by 2030, requiring scientific efforts from different research fields. The past decade has seen a much increased use of phylogeographic and phylodynamic analyses to study the evolution and spread of rabies virus. We here review published studies in these research areas, making a distinction between the geographic resolution associated with the available sequence data. We pay special attention to environmental factors that these studies found to be relevant to the spread of rabies virus. Importantly, we highlight a knowledge gap in terms of applying these methods when all required data were available but not fully exploited. We conclude with an overview of recent methodological developments that have yet to be applied in phylogeographic and phylodynamic analyses of rabies virus.

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.

Phylogeographic analysis of foot-and-mouth disease virus serotype O dispersal and associated drivers in East Africa

The continued endemicity of foot and mouth disease virus (FMDV) in East Africa has significant implications for livestock production and poverty reduction, yet its complex epidemiology in endemic settings remains poorly understood. Identifying FMDV dispersal routes and drivers of transmission is key to improved control strategies. Environmental heterogeneity and anthropogenic drivers (e.g., demand for animal products) can impact viral spread by influencing host movements. Here, we utilized FMDV serotype O VP1 genetic sequences and corresponding spatiotemporal data in order to (i) infer the recent dispersal history, and (II) investigate the impact of external factors (cattle density, human population density, proximity to livestock markets, and drought) on dispersal velocity, location, and direction of FMDV serotype O in East Africa. We identified statistical evidence of long-distance transmission events, and we found that FMDV serotype O tends to remain circulating in areas of high cattle density, high human population density, and in close proximity to livestock markets. The latter two findings highlight the influence of anthropogenic factors on FMDV serotype O spread in this region. These findings contribute to the understanding of FMDV epidemiology in East Africa and can help guide improved control measures.

Whole-genome sequencing and phylogenetic analysis of rabies viruses from Jordan

Human fatalities caused by rabies are rarely reported in Jordan; however, domestic animals are more likely to fall victim to rabies compared to wild animals, at least this is the case in Jordan due to the presence of canine rabies. In this study, twelve brain samples from domestic and wild animals suspected of being infected with rabies virus from different regions of Jordan were collected during 2019. Seven of them tested positive using the fluorescent antibody test and real-time SYBR RT-PCR assay. Five specimens were from stray dogs and two from foxes. The whole genome sequences were obtained from the positive samples. Sequence analysis showed that one dog virus from Al Quwaysimah city located in Amman governorate, was closely related to an Israeli strain belonging to a Cosmopolitan ME1a clade. The genomes of the remaining six viruses (four from dogs and two from foxes) collected from different areas of Jordan were genetically-related to each other and clustered together with sequences from Iran and Turkey; all belong to Cosmopolitan ME2 clade. These sequences were analyzed with six other Jordanian rabies virus nucleoprotein (N) gene sequences available in the public database, five of them belong to ME1a clade and one belongs to ME1b clade. Rabies virus whole genome data is scarce across the Middle East. This study provides a better understanding of the molecular epidemiology of rabies virus in the region.

In this study, we performed whole genome sequencing (WGS) for rabies virus (RABV) isolates from seven samples, five of which were of stray dogs, and the other two were from foxes. Specimens were collected from animals across Jordan, including Balqa, Amman, Irbid, Tafilah, and Madaba governorates. Six out of the seven isolates were belonging to the Cosmopolitan ME2 clade, which related to the Iranian and Turkish sequences. This is not the case previously, where the majority of the Jordanian isolates belong to Cosmopolitan ME1a clade and closely related to the sequences from Israel. This shift might be due to the applied regulations across borders between Jordan and Israel. Besides the growth in travel and trade movement between Jordan and Turkey, where the latter is a border country with Iran. These collected data, where such studies are not common in the Middle East countries, will enhance our understanding of the RABV evolution and epidemiology in the region for rapid and effective response for rabies virus outbreaks.

Development of a Dissemination Platform for Spatiotemporal and Phylogenetic Analysis of Avian Infectious Bronchitis Virus

Infecting large portions of the global poultry populations, the avian infectious bronchitis virus (IBV) remains a major economic burden in North America. With more than 30 serotypes globally distributed, Arkansas, Connecticut, Delaware, Georgia, and Massachusetts are among the most predominant serotypes in the United States. Even though vaccination is widely used, the high mutation rate exhibited by IBV is continuously triggering the emergence of new viral strains and hindering control and prevention measures. For that reason, targeted strategies based on constantly updated information on the IBV circulation are necessary. Here, we sampled IBV-infected farms from one US state and collected and analyzed 65 genetic sequences coming from three different lineages along with the immunization information of each sampled farm. Phylodynamic analyses showed that IBV dispersal velocity was 12.3 km/year. The majority of IBV infections appeared to have derived from the introduction of the Arkansas DPI serotype, and the Arkansas DPI and Georgia 13 were the predominant serotypes. When analyzed against IBV sequences collected across the United States and deposited in the GenBank database, the most likely viral origin of our sequences was from the states of Alabama, Georgia, and Delaware. Information about vaccination showed that the MILDVAC-MASS+ARK vaccine was applied on 26% of the farms. Using a publicly accessible open-source tool for real-time interactive tracking of pathogen spread and evolution, we analyzed the spatiotemporal spread of IBV and developed an online reporting dashboard. Overall, our work demonstrates how the combination of genetic and spatial information could be used to track the spread and evolution of poultry diseases, providing timely information to the industry. Our results could allow producers and veterinarians to monitor in near-real time the current IBV strain circulating, making it more informative, for example, in vaccination-related decisions.

Mathematical modelling and phylodynamics for the study of dog rabies dynamics and control: A scoping review

BACKGROUND

Rabies is a fatal yet vaccine-preventable disease. In the last two decades, domestic dog populations have been shown to constitute the predominant reservoir of rabies in developing countries, causing 99% of human rabies cases. Despite substantial control efforts, dog rabies is still widely endemic and is spreading across previously rabies-free areas. Developing a detailed understanding of dog rabies dynamics and the impact of vaccination is essential to optimize existing control strategies and developing new ones. In this scoping review, we aimed at disentangling the respective contributions of mathematical models and phylodynamic approaches to advancing the understanding of rabies dynamics and control in domestic dog populations. We also addressed the methodological limitations of both approaches and the remaining issues related to studying rabies spread and how this could be applied to rabies control.

METHODOLOGY/PRINCIPAL FINDINGS

We reviewed how mathematical modelling of disease dynamics and phylodynamics have been developed and used to characterize dog rabies dynamics and control. Through a detailed search of the PubMed, Web of Science, and Scopus databases, we identified a total of n = 59 relevant studies using mathematical models (n = 30), phylodynamic inference (n = 22) and interdisciplinary approaches (n = 7). We found that despite often relying on scarce rabies epidemiological data, mathematical models investigated multiple aspects of rabies dynamics and control. These models confirmed the overwhelming efficacy of massive dog vaccination campaigns in all settings and unraveled the role of dog population structure and frequent introductions in dog rabies maintenance. Phylodynamic approaches successfully disentangled the evolutionary and environmental determinants of rabies dispersal and consistently reported support for the role of reintroduction events and human-mediated transportation over long distances in the maintenance of rabies in endemic areas. Potential biases in data collection still need to be properly accounted for in most of these analyses. Finally, interdisciplinary studies were determined to provide the most comprehensive assessments through hypothesis generation and testing. They also represent new avenues, especially concerning the reconstruction of local transmission chains or clusters through data integration.

CONCLUSIONS/SIGNIFICANCE

Despite advances in rabies knowledge, substantial uncertainty remains regarding the mechanisms of local spread, the role of wildlife in dog rabies maintenance, and the impact of community behavior on the efficacy of control strategies including vaccination of dogs. Future integrative approaches that use phylodynamic analyses and mechanistic models within a single framework could take full advantage of not only viral sequences but also additional epidemiological information as well as dog ecology data to refine our understanding of rabies spread and control. This would represent a significant improvement on past studies and a promising opportunity for canine rabies research in the frame of the One Health concept that aims to achieve better public health outcomes through cross-sector collaboration.

Sampling bias and model choice in continuous phylogeography: Getting lost on a random walk

Phylogeographic inference allows reconstruction of past geographical spread of pathogens or living organisms by integrating genetic and geographic data. A popular model in continuous phylogeography—with location data provided in the form of latitude and longitude coordinates—describes spread as a Brownian motion (Brownian Motion Phylogeography, BMP) in continuous space and time, akin to similar models of continuous trait evolution. Here, we show that reconstructions using this model can be strongly affected by sampling biases, such as the lack of sampling from certain areas. As an attempt to reduce the effects of sampling bias on BMP, we consider the addition of sequence-free samples from under-sampled areas. While this approach alleviates the effects of sampling bias, in most scenarios this will not be a viable option due to the need for prior knowledge of an outbreak’s spatial distribution. We therefore consider an alternative model, the spatial Λ-Fleming-Viot process (ΛFV), which has recently gained popularity in population genetics. Despite the ΛFV’s robustness to sampling biases, we find that the different assumptions of the ΛFV and BMP models result in different applicabilities, with the ΛFV being more appropriate for scenarios of endemic spread, and BMP being more appropriate for recent outbreaks or colonizations.

Phylogeography studies past location and migration using information from current geographic locations of genetic sequences. For example, phylogeography can be used to reconstruct the history of geographical spread of an outbreak using the genetic sequences of the pathogen collected at different times and locations. Here, we investigate the effects of different model assumptions on phylogeographic inference. In particular, we examine the effects of the strategy used to collect samples. We show that sample collection biases can have a strong impact on the quality of phylogeographic reconstruction: geographically biased sampling scheme can be very detrimental for popular continuous phylogeography models. We consider different ways to counter these effects, from utilising alternative phylogeographic models, to the inclusion of partially informative samples (known cases without genetic sequences). While these strategies do alleviate the effects of sampling biases, they also lead to considerable additional computational burden. We also investigate the intrinsic differences of different phylogeographic models, and their effects on reconstructed patterns in different scenarios.

[New aspects of rabies control]

Rabies still causes about 60,000 human deaths per year, mainly in poor populations in Africa and Asia. However, since Louis Pasteur developed the first vaccine 130 years ago, prophylactic measures have been considerably improved and simplified. They now consist of the vaccine combined with purified rabies immunoglobulins of equine or human origin. In general, however, post-exposure prophylaxis protocols are long and expensive. Furthermore, the immunoglobulins used for associated serotherapy are costly and not widely available in developing countries. Approaches have been developed to deal with these two issues that offer hope for a paradigm shift for the benefit of exposed populations. Finally, mass rabies vaccination in dogs, which are the most cost-effective measure for preventing rabies in humans, are difficult to implement and sometimes have moderate effectiveness. The identification and analysis of the epidemiological drivers conditioning the circulation of the virus in dog populations allow a better understanding of the key control points that need to be associated with these campaigns for a better efficacy.

Update on lyssaviruses and rabies: will past progress play as prologue in the near term towards future elimination?

Rabies is an ancient, much-feared, and neglected infectious disease. Caused by pathogens in the family Rhabdoviridae, genus Lyssavirus, and distributed globally, this viral zoonosis results in tens of thousands of human fatalities and millions of exposures annually. All mammals are believed susceptible, but only certain taxa act as reservoirs. Dependence upon direct routing to, replication within, and passage from the central nervous system serves as a basic viral strategy for perpetuation. By a combination of stealth and subversion, lyssaviruses are quintessential neurotropic agents and cause an acute, progressive encephalitis. No treatment exists, so prevention is the key. Although not a disease considered for eradication, something of a modern rebirth has been occurring within the field as of late with regard to detection, prevention, and management as well as applied research. For example, within the past decade, new lyssaviruses have been characterized; sensitive and specific diagnostics have been optimized; pure, potent, safe, and efficacious human biologics have improved human prophylaxis; regional efforts have controlled canine rabies by mass immunization; wildlife rabies has been controlled by oral rabies vaccination over large geographic areas in Europe and North America; and debate has resumed over the controversial topic of therapy. Based upon such progress to date, there are certain expectations for the next 10 years. These include pathogen discovery, to uncover additional lyssaviruses in the Old World; laboratory-based surveillance enhancement by simplified, rapid testing; anti-viral drug appearance, based upon an improved appreciation of viral pathobiology and host response; and improvements to canine rabies elimination regionally throughout Africa, Asia, and the Americas by application of the best technical, organizational, economic, and socio-political practices. Significantly, anticipated Gavi support will enable improved access of human rabies vaccines in lesser developed countries at a national level, with integrated bite management, dose-sparing regimens, and a 1 week vaccination schedule.

Epidemiological hypothesis testing using a phylogeographic and phylodynamic framework

Computational analyses of pathogen genomes are increasingly used to unravel the dispersal history and transmission dynamics of epidemics. Here, we show how to go beyond historical reconstructions and use spatially-explicit phylogeographic and phylodynamic approaches to formally test epidemiological hypotheses. We illustrate our approach by focusing on the West Nile virus (WNV) spread in North America that has substantially impacted public, veterinary, and wildlife health. We apply an analytical workflow to a comprehensive WNV genome collection to test the impact of environmental factors on the dispersal of viral lineages and on viral population genetic diversity through time. We find that WNV lineages tend to disperse faster in areas with higher temperatures and we identify temporal variation in temperature as a main predictor of viral genetic diversity through time. By contrasting inference with simulation, we find no evidence for viral lineages to preferentially circulate within the same migratory bird flyway, suggesting a substantial role for non-migratory birds or mosquito dispersal along the longitudinal gradient.

Global Spread of the B5 Subgenotype EV-A71 and the Phylogeographical Analysis of Chinese Migration Events

The subgenotype B5 of EV-A71 is a widely circulating subgenotype that frequently spreads across the globe. Several outbreaks have occurred in nations, such as Malaysia, Thailand, Vietnam, and Japan. Appearing first in Taiwan, China, the subgenotype has been frequently reported in mainland of China even though no outbreaks have been reported so far. The current study reconstructed the migration of the B5 subgenotype of EV-A71 in China via phylogeographical analysis. Furthermore, we investigated its population dynamics in order to draw more credible inferences. Following a dataset cleanup of B5 subgenotype of EV-A71, we detected earlier B5 subgenotypes of EV-A71 sequences that had been circulating in Malaysia and Singapore since the year 2000, which was before the 2003 outbreak that occurred in Sarawak. The Bayesian inference indicated that the most recent common ancestor of B5 subgenotype EV-A71 appeared in September, 1994 (1994.75). With respect to the overall prevalence, geographical reconstruction revealed that the B5 subgenotype EV-A71 originated singly from single-source cluster and subsequently developed several active lineages. Based on a large amount of data that was accumulated, we conclude that the appearance of the B5 subgenotype of EV-A71 in mainland of China was mainly due to multiple migrations from different origins.

Genomic Epidemiology, Evolution, and Transmission Dynamics of Porcine Deltacoronavirus

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has shown once again that coronavirus (CoV) in animals are potential sources for epidemics in humans. Porcine deltacoronavirus (PDCoV) is an emerging enteropathogen of swine with a worldwide distribution. Here, we implemented and described an approach to analyze the epidemiology of PDCoV following its emergence in the pig population. We performed an integrated analysis of full genome sequence data from 21 newly sequenced viruses, along with comprehensive epidemiological surveillance data collected globally over the last 15 years. We found four distinct phylogenetic lineages of PDCoV, which differ in their geographic circulation patterns. Interestingly, we identified more frequent intra- and interlineage recombination and higher virus genetic diversity in the Chinese lineages compared with the USA lineage where pigs are raised in different farming systems and ecological environments. Most recombination breakpoints are located in the ORF1ab gene rather than in genes encoding structural proteins. We also identified five amino acids under positive selection in the spike protein suggesting a role for adaptive evolution. According to structural mapping, three positively selected sites are located in the N-terminal domain of the S1 subunit, which is the most likely involved in binding to a carbohydrate receptor, whereas the other two are located in or near the fusion peptide of the S2 subunit and thus might affect membrane fusion. Finally, our phylogeographic investigations highlighted notable South-North transmission as well as frequent long-distance dispersal events in China that could implicate human-mediated transmission. Our findings provide new insights into the evolution and dispersal of PDCoV that contribute to our understanding of the critical factors involved in CoVs emergence.

The World Rabies Day 2020: Collaborate and Vaccinate

Rabies constantly kills 59,000 people annually, mostly in Asia and Africa. Rabies, which is responsible for 99% of human rabies cases, is totally preventable by standard vaccinations. In 2015, a global call for action was made by the World Health Organization, World Organization for Animal Health, Food and Agriculture Organization of the United Nations, and the Global Alliance for Rabies Control to join forces toward the elimination of dog-transmitted human rabies by the year 2030. All the tools and protocols to reach that target are readily available, and the feasibility of dog rabies elimination has been proven. Countries should drive the changes needed to engage into this global movement. Certainly, countries in the Eastern Mediterranean Region require taking more critical steps to reach the rabies elimination target by 2030. The international awareness campaign of the World Rabies Day is an excellent occasion to assess challenges and opportunities toward rabies elimination.

Porcine reproductive and respiratory syndrome virus dissemination across pig production systems in the United States

Porcine reproductive and respiratory syndrome virus (PRRSV) remains widespread in the North American pig population. Despite improvements in virus characterization, it is unclear whether PRRSV infections are a product of viral circulation within production systems (local) or across production systems (external). Here, we examined the local and external dissemination dynamics of PRRSV and the processes facilitating its spread in three production systems. Overall, PRRSV genetic diversity has declined since 2018, while phylodynamic results support frequent external transmission. We found that PRRSV dissemination predominantly occurred mostly through transmission between farms of different production companies for several months, especially from November until May, a timeframe already established as PRRSV season. Although local PRRSV dissemination occurred mainly through regular pig flow (from sow to nursery and then to finisher farms), an important flux of PRRSV dissemination also occurred in the opposite direction, from finisher to sow and nursery farms, highlighting the importance of downstream farms as sources of the virus. Our results also showed that farms with pig densities of 500 to 1,000 pig/km² and farms located at a range within 0.5 km and 0.7 km from major roads were more likely to be infected by PRRSV, whereas farms at an elevation of 41 to 61 meters and surrounded by denser vegetation were less likely to be infected, indicating their role as dissemination barriers. In conclusion, our results demonstrate that external dissemination was intense, and reinforce the importance of farm proximity on PRRSV spread. Thus, consideration of farm location, geographic characteristics and animal densities across production systems may help to forecast PRRSV collateral dissemination.

Rapid in-country sequencing of whole virus genomes to inform rabies elimination programmes

Genomic surveillance is an important aspect of contemporary disease management but has yet to be used routinely to monitor endemic disease transmission and control in low- and middle-income countries. Rabies is an almost invariably fatal viral disease that causes a large public health and economic burden in Asia and Africa, despite being entirely vaccine preventable. With policy efforts now directed towards achieving a global goal of zero dog-mediated human rabies deaths by 2030, establishing effective surveillance tools is critical. Genomic data can provide important and unique insights into rabies spread and persistence that can direct control efforts. However, capacity for genomic research in low- and middle-income countries is held back by limited laboratory infrastructure, cost, supply chains and other logistical challenges. Here we present and validate an end-to-end workflow to facilitate affordable whole genome sequencing for rabies surveillance utilising nanopore technology. We used this workflow in Kenya, Tanzania and the Philippines to generate rabies virus genomes in two to three days, reducing costs to approximately £60 per genome. This is over half the cost of metagenomic sequencing previously conducted for Tanzanian samples, which involved exporting samples to the UK and a three- to six-month lag time. Ongoing optimization of workflows are likely to reduce these costs further. We also present tools to support routine whole genome sequencing and interpretation for genomic surveillance. Moreover, combined with training workshops to empower scientists in-country, we show that local sequencing capacity can be readily established and sustainable, negating the common misperception that cutting-edge genomic research can only be conducted in high resource laboratories. More generally, we argue that the capacity to harness genomic data is a game-changer for endemic disease surveillance and should precipitate a new wave of researchers from low- and middle-income countries.

Rapid in-country sequencing of whole virus genomes to inform rabies elimination programmes

Genomic surveillance is an important aspect of contemporary disease management but has yet to be used routinely to monitor endemic disease transmission and control in low- and middle-income countries. Rabies is an almost invariably fatal viral disease that causes a large public health and economic burden in Asia and Africa, despite being entirely vaccine preventable. With policy efforts now directed towards achieving a global goal of zero dog-mediated human rabies deaths by 2030, establishing effective surveillance tools is critical. Genomic data can provide important and unique insights into rabies spread and persistence that can direct control efforts. However, capacity for genomic research in low- and middle-income countries is held back by limited laboratory infrastructure, cost, supply chains and other logistical challenges. Here we present and validate an end-to-end workflow to facilitate affordable whole genome sequencing for rabies surveillance utilising nanopore technology. We used this workflow in Kenya, Tanzania and the Philippines to generate rabies virus genomes in two to three days, reducing costs to approximately £60 per genome. This is over half the cost of metagenomic sequencing previously conducted for Tanzanian samples, which involved exporting samples to the UK and a three- to six-month lag time. Ongoing optimization of workflows are likely to reduce these costs further. We also present tools to support routine whole genome sequencing and interpretation for genomic surveillance. Moreover, combined with training workshops to empower scientists in-country, we show that local sequencing capacity can be readily established and sustainable, negating the common misperception that cutting-edge genomic research can only be conducted in high resource laboratories. More generally, we argue that the capacity to harness genomic data is a game-changer for endemic disease surveillance and should precipitate a new wave of researchers from low- and middle-income countries.