Elucidating the phylodynamics of endemic rabies virus in eastern Africa using whole-genome sequencing

Dispersion and diversification of Lyssavirus rabies transmitted from haematophagous bats Desmodus rotundus: a phylogeographical study

Rabies is worldwide zoonosis caused by Lyssavirus rabies (RABV) a RNA negative sense virus with low level of fidelity during replication cycle. Nucleoprotein of RABV is the most conserved between all five proteins of the virus and is the most used gene for phylogenetic and phylogeographic studies. Despite of rabies been very important in Public Health concern, it demands continuous prophylactic care for herbivores with economic interest, such as cattle and horses. The main transmitter of RABV for these animals in Brazil is the hematophagous bats Desmodus rotundus. The aim of this study was to determine the dispersion over time and space of RABV transmitted by D. rotundus. Samples of RABV from the State of São Paulo (SP), Southeast Brazil isolated from the central nervous system (CNS) of cattle, were submitted to RNA extraction, RT-PCR, sequencing and phylogeographic analyzes with BEAST (Bayesian Evolutionary Analysis Sampling Trees) v 2.5 software. Was possible to identify high rate of diversification in starts sublineages of RABV what are correlated with a behavior of D. rotundus, the main transmitter of rabies to cattle. This study also highlights the importance of continuous monitoring of genetic lineages of RABV in Brazil.

Integrating full and partial genome sequences to decipher the global spread of canine rabies virus

Despite the rapid growth in viral genome sequencing, statistical methods face challenges in handling historical viral endemic diseases with large amounts of underutilized partial sequence data. We propose a phylogenetic pipeline that harnesses both full and partial viral genome sequences to investigate historical pathogen spread between countries. Its application to rabies virus (RABV) yields precise dating and confident estimates of its geographic dispersal. By using full genomes and partial sequences, we reduce both geographic and genetic biases that often hinder studies that focus on specific genes. Our pipeline reveals an emergence of the present canine-mediated RABV between years 1301 and 1403 and reveals regional introductions over a 700-year period. This geographic reconstruction enables us to locate episodes of human-mediated introductions of RABV and examine the role that European colonization played in its spread. Our approach enables phylogeographic analysis of large and genetically diverse data sets for many viral pathogens.

From Field Tests to Molecular Tools-Evaluating Diagnostic Tests to Improve Rabies Surveillance in Namibia

Rabies is endemic in Namibia and is present both in wildlife carnivores and domestic free-roaming dogs. The disease thus represents a challenge for public human and veterinary disease control. Namibia has implemented a national strategic plan to control rabies and the country's activities are supported by international organizations. To this end, rabies diagnosis at the Central Veterinary Laboratory (CVL) was improved in the frame of a World Organization for Animal Health (WOAH) laboratory twinning program: from practical sampling techniques and the use of lateral flow devices to a novel universal and discriminatory quantitative real-time Reverse transcription polymerase chain reaction (RT-qPCR), which easily identify dog-associated rabies viruses. The procedures applied and the results can be used as a template to improve rabies laboratory diagnosis.

Identifying high-risk areas for dog-mediated rabies using Bayesian spatial regression

Despite ongoing control efforts, rabies remains an endemic zoonotic disease in many countries. Determining high-risk areas and the space-time patterns of rabies spread, as it relates to epidemiologically important factors, can support policymakers and program managers alike to develop evidence-based targeted surveillance and control programs. In this One Health approach which selected Thailand as the example site, the location-based risk of contracting dog-mediated rabies by both human and animal populations was quantified using a Bayesian spatial regression model. Specifically, a conditional autoregressive (CAR) Bayesian zero-inflated Poisson (ZIP) regression was fitted to the reported human and animal rabies case counts of each district, from the 2012–2017 period. The human population was used as an offset. The epidemiologically important factors hypothesized as risk modifiers and therefore tested as predictors included: number of dog bites/attacks, the population of dogs and cats, number of Buddhist temples, garbage dumps, animal vaccination, post-exposure prophylaxis, poverty, and shared administrative borders. Disparate sources of data were used to improve the estimated associations and predictions. Model performance was assessed using cross-validation. Results suggested that accounting for the association between human and animal rabies with number of dog bites/attacks, number of owned and un-owned dogs; shared country borders, number of Buddhist temples, poverty levels, and accounting for spatial dependence between districts, may help to predict the risk districts for dog-mediated rabies in Thailand. The fitted values of the spatial regression were mapped to illustrate the risk of dog-mediated rabies. The cross-validation indicated an adequate performance of the spatial regression model (AUC = 0.81), suggesting that had this spatial regression approach been used to identify districts at risk in 2015, the cases reported in 2016/17 would have been predicted with model sensitivity and specificity of 0.71 and 0.80, respectively. While active surveillance is ideal, this approach of using multiple data sources to improve risk estimation may inform current rabies surveillance and control efforts including determining rabies-free zones, and the roll-out of human post-exposure prophylaxis and anti-rabies vaccines for animals in determining high-risk areas.

•

Bayesian spatial regression was used to quantify location-based risk of dog-mediated rabies

•

Available and publicly accessible data from disparate sources were gathered

•

Risk was estimated using the association between

•

Risk estimates were compared over time to determine the prediction ability

•

Study suggests while active surveillance is ideal, using multiple data sources may improve risk estimation

Epidemiological Interface of Sylvatic and Dog Rabies in the North West Province of South Africa

Rabies is a viral zoonosis that causes an estimated 59,000 preventable human fatalities every year. While more than 120 countries remain endemic for dog-mediated rabies, the burden is the highest in Africa and Asia where 99% of human rabies cases are caused by domestic dogs. One such rabies-endemic country is South Africa where an estimated 42 preventable human deaths occur every year. Although canine rabies had been well described for most of the provinces in South Africa, the epidemiology of rabies within the North West Province had not been well defined prior to this investigation. As such, the aim of this study was to use nucleotide sequence analyses to characterise the extant molecular epidemiology of rabies in the North West Province of South Africa—with specific focus on the interface between dogs and sylvatic species. To this end, Rabies lyssavirus isolates originating from the North West Province were subjected to molecular epidemiological analyses relying on the Bayesian Markov Chain Monte Carlo methodology on two distinct gene regions, viz. the G-L intergenic region and partial nucleoprotein gene. Our results provided strong evidence in support of an endemic cycle of canine rabies in the East of the province, and three independent endemic cycles of sylvatic rabies spread throughout the province. Furthermore, evidence of specific events of virus spill-over between co-habiting sylvatic species and domestic dogs was found. These results suggest that the elimination of canine-mediated rabies from the province will rely not only on eliminating the disease from the dog populations, but also from the co-habiting sylvatic populations using oral rabies vaccination campaigns.

Making genomic surveillance deliver: A lineage classification and nomenclature system to inform rabies elimination

The availability of pathogen sequence data and use of genomic surveillance is rapidly increasing. Genomic tools and classification systems need updating to reflect this. Here, rabies virus is used as an example to showcase the potential value of updated genomic tools to enhance surveillance to better understand epidemiological dynamics and improve disease control. Previous studies have described the evolutionary history of rabies virus, however the resulting taxonomy lacks the definition necessary to identify incursions, lineage turnover and transmission routes at high resolution. Here we propose a lineage classification system based on the dynamic nomenclature used for SARS-CoV-2, defining a lineage by phylogenetic methods for tracking virus spread and comparing sequences across geographic areas. We demonstrate this system through application to the globally distributed Cosmopolitan clade of rabies virus, defining 96 total lineages within the clade, beyond the 22 previously reported. We further show how integration of this tool with a new rabies virus sequence data resource (RABV-GLUE) enables rapid application, for example, highlighting lineage dynamics relevant to control and elimination programmes, such as identifying importations and their sources, as well as areas of persistence and routes of virus movement, including transboundary incursions. This system and the tools developed should be useful for coordinating and targeting control programmes and monitoring progress as countries work towards eliminating dog-mediated rabies, as well as having potential for broader application to the surveillance of other viruses.

Long-Read Sequencing and Hybrid Assembly for Genomic Analysis of Clinical Brucella melitensis Isolates

Brucella melitensis is a key etiological agent of brucellosis and has been increasingly subject to characterization using sequencing methodologies. This study aimed to investigate and compare short-read, long-read, and hybrid assemblies of B. melitensis. Eighteen B. melitensis isolates from Southern Israel were sequenced using Illumina and the Oxford Nanopore (ONP) MinION, and hybrid assemblies were generated with ONP long reads scaffolded on Illumina short reads. Short reads were assembled with INNUca with SPADes, long reads and hybrid with dragonflye. Abricate with the virulence factor database (VFDB) and in silico PCR (for the genes BetB, BPE275, BSPB, manA, mviN, omp19, perA, PrpA, VceC, and ureI) were used for identifying virulence genes, and a total of 61 virulence genes were identified in short-read, long-read, and hybrid assemblies of all 18 isolates. The phylogenetic analysis using long-read assemblies revealed several inconsistencies in cluster assignment as compared to using hybrid and short-read assemblies. Overall, hybrid assembly provided the most comprehensive data, and stand-alone short-read sequencing provided comparable data to stand-alone long-read sequencing regarding virulence genes. For genomic epidemiology studies, stand-alone ONP sequencing may require further refinement in order to be useful in endemic settings.

Evaluating the Impact of Anthropogenic Factors on the Dissemination of Contemporary Cosmopolitan, Arctic, and Arctic-like Rabies Viruses

Rabies is a globally prevalent viral zoonosis that causes 59,000 deaths per year and has important economic consequences. Most virus spread is associated with the migration of its primary hosts. Anthropogenic dissemination, mainly via the transportation of rabid dogs, shaped virus ecology a few hundred years ago and is responsible for several current outbreaks. A systematic analysis of aberrant long-distance events in the steppe and Arctic-like groups of rabies virus was performed using statistical (Bayesian) phylogeography and plots of genetic vs. geographic distances. The two approaches produced similar results but had some significant differences and complemented each other. No phylogeographic analysis could be performed for the Arctic group because polar foxes transfer the virus across the whole circumpolar region at high velocity, and there was no correlation between genetic and geographic distances in this virus group. In the Arctic-like group and the steppe subgroup of the cosmopolitan group, a significant number of known sequences (15-20%) was associated with rapid long-distance transfers, which mainly occurred within Eurasia. Some of these events have been described previously, while others have not been documented. Most of the recent long-distance transfers apparently did not result in establishing the introduced virus, but a few had important implications for the phylogeographic history of rabies. Thus, human-mediated long-distance transmission of the rabies virus remains a significant threat that needs to be addressed.

Reservoir dynamics of rabies in south-east Tanzania and the roles of cross-species transmission and domestic dog vaccination

Understanding the role of different species in the transmission of multi-host pathogens, such as rabies virus, is vital for effective control strategies. Across most of sub-Saharan Africa domestic dogs Canis familiaris are considered the reservoir for rabies, but the role of wildlife has been long debated. Here we explore the multi-host transmission dynamics of rabies across south-east Tanzania.Between January 2011 and July 2019, data on probable rabies cases were collected in the regions of Lindi and Mtwara. Hospital records of animal-bite patients presenting to healthcare facilities were used as sentinels for animal contact tracing. The timing, location and species of probable rabid animals were used to reconstruct transmission trees to infer who infected whom and the relative frequencies of within- and between-species transmission.During the study, 688 probable human rabies exposures were identified, resulting in 47 deaths. Of these exposures, 389 were from domestic dogs (56.5%) and 262 from jackals (38.1%). Over the same period, 549 probable animal rabies cases were traced: 303 in domestic dogs (55.2%) and 221 in jackals (40.3%), with the remainder in domestic cats and other wildlife species.Although dog-to-dog transmission was most commonly inferred (40.5% of transmission events), a third of inferred events involved wildlife-to-wildlife transmission (32.6%), and evidence suggested some sustained transmission chains within jackal populations.A steady decline in probable rabies cases in both humans and animals coincided with the implementation of widespread domestic dog vaccination during the first 6 years of the study. Following the lapse of this program, dog rabies cases began to increase in one of the northernmost districts. Synthesis and applications. In south-east Tanzania, despite a relatively high incidence of rabies in wildlife and evidence of wildlife-to-wildlife transmission, domestic dogs remain essential to the reservoir of infection. Continued dog vaccination alongside improved surveillance would allow a fuller understanding of the role of wildlife in maintaining transmission in this area. Nonetheless, dog vaccination clearly suppressed rabies in both domestic dog and wildlife populations, reducing both public health and conservation risks and, if sustained, has potential to eliminate rabies from this region.

Use of partial N-gene sequences as a tool to monitor progress on rabies control and elimination efforts in Ethiopia

Ethiopia is one of the African countries most affected by rabies. A coarse catalog of rabies viruses (RABV) was created as a benchmark to assess the impact of control and elimination activities. We evaluated a 726 bp amplicon at the end of the N-gene to infer viral lineages in circulation using maximum likelihood and Bayesian methods for phylogenetic reconstruction. We sequenced 228 brain samples from wild and domestic animals collected in five Ethiopian regions during 2010-2017. Results identified co-circulating RABV lineages that are causing recurrent spillover infections into wildlife and domestic animals. We found no evidence of importation of RABVs from other African countries or vaccine-induced cases in the area studied. A divergent RABV lineage might be involved in an independent rabies cycle in jackals. This investigation provides a feasible approach to assess rabies control and elimination efforts in resource-limited countries.

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.

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.

Identification of a novel Aichivirus D in sheep

A novel kobuvirus was found in diarrheal fecal samples of Tibetan sheep using a viral metagenomics approach, and a full kobuvirus genome was successfully obtained by RT-PCR from a diarrheal fecal sample. The full genomic sequence was 8485 nucleotides (nt) in length with a standard picornavirus genome organization. The novel genome shares 62.9% and 77.8% nt homology with Aichivirus D1 genotype strain 1-22-KoV, and Aichivirus D2 genotype strain 2-44-KoV, respectively. According to the species classification criteria of the International Committee on Taxonomy of Viruses (ICTV), the new kobuvirus belongs to Aichivirus species D. Interestingly, compared with 2 known Aichivirus D genotype strains, the novel Aichivirus D has unique amino acid substitutions in the 5'untranslated region (-UTR), VP0, VP3, and VP1, with a recombination event in the 2C region.These characteristics make the novel Aichivirus D cluster into an independent branch in the phylogenetic tree, suggesting that strain may represent a novel genotype in Aichivirus D. Moreover, the novel Aichivirus D was detected in 9.2% (18/195) of the sheep diarrheal fecal samples from 4 farms in 3 counties of the Qinghai Tibet Plateau in China. In addition, full-length VP0, VP3, and VP1 genes were successfully obtained from 12 samples from 4 farms, and phylogenetic analysis based on these genes revealed a unique evolutionary pattern for this novel Aichivirus D strain. This study identified a novel Aichivirus D that is circulating in sheep in Qinghai Tibet Plateau in China and these findings provide a better understanding of the epidemiologic and genetic evolution of kobuviruses.

Portable Rabies Virus Sequencing in Canine Rabies Endemic Countries Using the Oxford Nanopore MinION

As countries with endemic canine rabies progress towards elimination by 2030, it will become necessary to employ techniques to help plan, monitor, and confirm canine rabies elimination. Sequencing can provide critical information to inform control and vaccination strategies by identifying genetically distinct virus variants that may have different host reservoir species or geographic distributions. However, many rabies testing laboratories lack the resources or expertise for sequencing, especially in remote or rural areas where human rabies deaths are highest. We developed a low-cost, high throughput rabies virus sequencing method using the Oxford Nanopore MinION portable sequencer. A total of 259 sequences were generated from diverse rabies virus isolates in public health laboratories lacking rabies virus sequencing capacity in Guatemala, India, Kenya, and Vietnam. Phylogenetic analysis provided valuable insight into rabies virus diversity and distribution in these countries and identified a new rabies virus lineage in Kenya, the first published canine rabies virus sequence from Guatemala, evidence of rabies spread across an international border in Vietnam, and importation of a rabid dog into a state working to become rabies-free in India. Taken together, our evaluation highlights the MinION's potential for low-cost, high volume sequencing of pathogens in locations with limited resources.

Towards the elimination of dog-mediated rabies: development and application of an evidence-based management tool

BACKGROUND

International organizations advocate for the elimination of dog-mediated rabies, but there is only limited guidance on interpreting surveillance data for managing elimination programmes. With the regional programme in Latin America approaching elimination of dog-mediated rabies, we aimed to develop a tool to evaluate the programme's performance and generate locally-tailored rabies control programme management guidance to overcome remaining obstacles.

METHODS

We developed and validated a robust algorithm to classify progress towards rabies elimination within sub-national administrative units, which we applied to surveillance data from Brazil and Mexico. The method combines criteria that are easy to understand, including logistic regression analysis of case detection time series, assessment of rabies virus variants, and of incursion risk. Subjecting the algorithm to robustness testing, we further employed simulated data sub-sampled at differing levels of case detection to assess the algorithm's performance and sensitivity to surveillance quality.

RESULTS

Our tool demonstrated clear epidemiological transitions in Mexico and Brazil: most states progressed rapidly towards elimination, but a few regressed due to incursions and control lapses. In 2015, dog-mediated rabies continued to circulate in the poorest states, with foci remaining in only 1 of 32 states in Mexico, and 2 of 27 in Brazil, posing incursion risks to the wider region. The classification tool was robust in determining epidemiological status irrespective of most levels of surveillance quality. In endemic settings, surveillance would need to detect less than 2.5% of all circulating cases to result in misclassification, whereas in settings where incursions become the main source of cases the threshold detection level for correct classification should not be less than 5%.

CONCLUSION

Our tool provides guidance on how to progress effectively towards elimination targets and tailor strategies to local epidemiological situations, while revealing insights into rabies dynamics. Post-campaign assessments of dog vaccination coverage in endemic states, and enhanced surveillance to verify and maintain freedom in states threatened by incursions were identified as priorities to catalyze progress towards elimination. Our finding suggests genomic surveillance should become increasingly valuable during the endgame for discriminating circulating variants and pinpointing sources of incursions.

Using affinity propagation clustering for identifying bacterial clades and subclades with whole-genome sequences of Francisella tularensis

By combining a reference-independent SNP analysis and average nucleotide identity (ANI) with affinity propagation clustering (APC), we developed a significantly improved methodology allowing resolving phylogenetic relationships, based on objective criteria. These bioinformatics tools can be used as a general ruler to determine phylogenetic relationships and clustering of bacteria, exemplary done with Francisella (F.) tularensis. Molecular epidemiology of F. tularensis is currently assessed mostly based on laboratory methods and molecular analysis. The high evolutionary stability and the clonal nature makes Francisella ideal for subtyping with single nucleotide polymorphisms (SNPs). Sequencing and real-time PCR can be used to validate the SNP analysis. We investigate whole-genome sequences of 155 F. tularensis subsp. holarctica isolates. Phylogenetic testing was based on SNPs and average nucleotide identity (ANI) as reference independent, alignment-free methods taking small-scale and large-scale differences within the genomes into account. Especially the whole genome SNP analysis with kSNP3.0 allowed deciphering quite subtle signals of systematic differences in molecular variation. Affinity propagation clustering (APC) resulted in three clusters showing the known clades B.4, B.6, and B.12. These data correlated with the results of real-time PCR assays targeting canSNPs loci. Additionally, we detected two subtle sub-clusters. SplitsTree was used with standard-setting using the aligned SNPs from Parsnps. Together APC, HierBAPS, and SplitsTree enabled us to generate hypotheses about epidemiologic relationships between bacterial clusters and describing the distribution of isolates. Our data indicate that the choice of the typing technique can increase our understanding of the pathogenesis and transmission of diseases with the eventual for prevention. This is opening perspectives to be applied to other bacterial species. The data provide evidence that Germany might be the collision zone where the clade B.12, also known as the East European clade, overlaps with the clade B.6, also known as the Iberian clade. Described methods allow generating a new, more detailed perspective for F. tularensis subsp. holarctica phylogeny. These results may encourage to determine phylogenetic relationships and clustering of other bacteria the same way.

Complete Coding Sequences of 23 South African Domestic and Wildlife Rabies Viruses

South African rabies viruses originating from dogs and jackals (canid viruses) are closely related and highlight cross-species transmission events between the two canine species. Rabies due to the canid lyssavirus variant is a significant public health matter in this country. The complete coding sequences of 23 canid lyssaviruses from South Africa are reported here.

South African rabies viruses originating from dogs and jackals (canid viruses) are closely related and highlight cross-species transmission events between the two canine species. Rabies due to the canid lyssavirus variant is a significant public health matter in this country. The complete coding sequences of 23 canid lyssaviruses from South Africa are reported here.

Molecular study of rabies virus in slaughtered dogs in Billiri and Kaltungo local government areas of Gombe state, Nigeria

Rabies is one of the most dreadful diseases and a major viral zoonosis which has been shown to cause an almost 100% fatality rate in infected victims. It is characterized by acute progressive encephalitis in mammals. This study determined the genotypic characteristics of rabies virus in dogs slaughtered for human consumption based on sequence of a fragment of nucleoprotein gene. Brain tissues were collected from 50 dogs slaughtered in Billiri and Kaltungo Local Government Areas of Gombe State, Nigeria. Direct fluorescent antibody test (DFAT) was used to screen for the presence of rabies virus antigen. Viral RNA isolated from DFAT positive brain tissues were subjected to the reverse transcription polymerase chain reaction (RT-PCR) followed by sequencing of the amplicons. Maximum Likelihood (ML) was used to construct a phylogenetic tree for sequences obtained with 1000 bootstrap replicates. The DFAT detected rabies antigen in 3 (6%) of the 50 dog brain tissues, from which 1 (2%) was positive by RT-PCR. ML phylogeny approach of the nucleotide sequences inferred members as originating lyssavirus genus and dog species. Essentially, MK234794 in this study displayed 99.3% sequence similarity with other related rabies viruses in the Africa 2 cluster (Nigeria, Cameroon, Chad and Niger). Interestingly, MK234794 showed no cluster relation with the Africa 1a, 1b, 3 and Africa 4 clades, respectively. This indicates there is in-country and trans-boundary circulation of the rabies viruses with no co-circulation between the Africa lineages, especially as dogs are continuously being traded due to consumption of dog meat in West Africa. This finding has given additional insight into the molecular epidemiology of rabies virus in Nigeria, therefore providing more baseline information for future design of rabies control programs in the country.

The Expectations and Challenges of Wildlife Disease Research in the Era of Genomics: Forecasting with a Horizon Scan-like Exercise

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented "Big Data" tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural "Genomics of Disease in Wildlife" workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed 5 themes considered to be the most pressing to the application of genomics in wildlife disease research: 1) "Improving communication," 2) "Methodological and analytical advancements," 3) "Translation into practice," 4) "Integrating landscape ecology and genomics," and 5) "Emerging new questions." Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

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.

Evaluation of Direct Metagenomics and Target Enriched Approaches for High-throughput Sequencing of Field Rabies Viruses

INTRODUCTION

High-throughput sequencing (HTS) identifies random viral fragments in environmental samples metagenomically. High reliability gains it broad application in virus evolution, host-virus interaction, and pathogenicity studies. Deep sequencing of field samples with content of host genetic material and bacteria often produces insufficient data for metagenomics and must be preceded by target enrichment. The main goal of the study was the evaluation of HTS for complete genome sequencing of field-case rabies viruses (RABVs).

MATERIAL AND METHODS

The material was 23 RABVs isolated mainly from red foxes and one European bat lyssavirus-1 isolate propagated in neuroblastoma cells. Three methods of RNA isolation were tested for the direct metagenomics and RABV-enriched approaches. Deep sequencing was performed with a MiSeq sequencer (Illumina) and reagent v3 kit. Bioinformatics data were evaluated by Kraken and Centrifuge software and de novo assembly was done with metaSPAdes.

RESULTS

Testing RNA extraction procedures revealed the deep sequencing scope superiority of the combined TRIzol/column method. This HTS methodology made it possible to obtain complete genomes of all the RABV isolates collected in the field. Significantly greater rates of RABV genome coverages (over 5,900) were obtained with RABV enrichment. Direct metagenomic studies sequenced the full length of 6 out of 16 RABV isolates with a medium coverage between 1 and 71.

CONCLUSION

Direct metagenomics gives the most realistic illustration of the field sample microbiome, but with low coverage. For deep characterisation of viruses, e.g. for spatial and temporal phylogeography during outbreaks, target enrichment is recommended as it covers sequences much more completely.

Phylodynamics reveals extinction-recolonization dynamics underpin apparently endemic vampire bat rabies in Costa Rica

Variation in disease incidence in wildlife is often assumed to reflect environmental or demographic changes acting on an endemic pathogen. However, apparent endemicity might instead arise from spatial processes that are challenging to identify from traditional data sources including time series and field studies. Here, we analysed longitudinal sequence data collected from rabies virus outbreaks over 14 years in Costa Rica, a Central American country that has recorded continuous vampire bat-transmitted rabies outbreaks in humans and livestock since 1985. We identified five phylogenetically distinct lineages which shared most recent common ancestors with viruses from North and South America. Bayesian phylogeographic reconstructions supported bidirectional viral dispersals involving countries to the north and south of Costa Rica at different time points. Within Costa Rica, viruses showed little contemporaneous spatial overlap and no lineage was detected across all years of surveillance. Statistical models suggested that lineage disappearances were more likely to be explained by viral extinctions than undetected viral circulation. Our results highlight the importance of international viral dispersal for shaping the burden of rabies in Costa Rica, suggest a Central American corridor of rabies virus invasions between continents, and show that apparent disease endemicity may arise through recurrent pathogen extinctions and reinvasions which can be readily detected in relatively small datasets by joining phylodynamic and modelling approaches.

A hundred years of rabies in Kenya and the strategy for eliminating dog-mediated rabies by 2030

Background: Rabies causes an estimated 59,000 human deaths annually. In Kenya, rabies was first reported in a dog in 1912, with the first human case reported in 1928. Here we examine retrospective rabies data in Kenya for the period 1912 - 2017 and describe the spatial and temporal patterns of rabies occurrence in the country. Additionally, we detail Kenya's strategy for the elimination of dog-mediated human rabies by 2030. Methods: Data on submitted samples and confirmed cases in humans, domestic animals and wildlife were obtained from Kenya's Directorate of Veterinary Services. These data were associated with the geographical regions where the samples originated, and temporal and spatial trends examined. Results: Between 1912 and the mid 1970's, rabies spread across Kenya gradually, with fewer than 50 cases reported per year and less than half of the 47 counties affected. Following an outbreak in the mid 1970's, rabies spread rapidly to more than 85% of counties, with a 4 fold increase in the percent positivity of samples submitted and number of confirmed rabies cases. Since 1958, 7,584 samples from domestic animals (93%), wildlife (5%), and humans (2%) were tested. Over two-thirds of all rabies cases came from six counties, all in close proximity to veterinary diagnostic laboratories, highlighting a limitation of passive surveillance. Conclusions: Compulsory annual dog vaccinations between 1950's and the early 1970's slowed rabies spread. The rapid spread with peak rabies cases in the 1980's coincided with implementation of structural adjustment programs privatizing the veterinary sector leading to breakdown of rabies control programs. To eliminate human deaths from rabies by 2030, Kenya is implementing a 15-year step-wise strategy based on three pillars: a) mass dog vaccination, b) provision of post-exposure prophylaxis and public awareness and c) improved surveillance for rabies in dogs and humans with prompt responses to rabies outbreaks.

Epidemiological aspects of the persistent transmission of rabies during an outbreak (2010 - 2017) in Harare, Zimbabwe

Canine-mediated human rabies is endemic to the entire African continent, where the disease burden is often highest in rural communities of resource-limited countries. In this study, we analysed an animal rabies outbreak, which had persisted since 2010 in the predominantly metropolitan capital city of Zimbabwe, Harare. As rabies is considered to disproportionally affect rural communities, the persistence of urban rabies in this metropolitan setting is of interest. In order to gain an improved understanding of the epidemiology of the outbreak under investigation, we utilised both routine surveillance data that had been collected during the first eight years of the outbreak (2010-2017), as well as molecular epidemiological analyses relying on the Bayesian Markov Chain Monte Carlo methodology. This approach allowed us to characterize virus transmission by identifying specific suburbs within the city limits where persistent disease transmission took place, while also confirming that immunologically naïve dogs were the most likely reservoir species in and around the city. In addition to gaining an improved local understanding of the outbreak, we are also able to infer that rabies was likely introduced to the city in 2010 when a rabid animal was moved from the north-east of Zimbabwe into Harare-resulting in an epizootic event. The work presented here not only showcased the value of combining conventional and molecular epidemiological data, but also highlighted the importance of maintaining rabies vaccination coverage and continued public awareness in urban areas where the risk appears to be low. By educating the general population on rabies and relying on owners to bring their companion animals to strategically placed vaccination points, the control and elimination of rabies from Harare may be feasible.

Transmission dynamics of re-emerging rabies in domestic dogs of rural China

Despite ongoing efforts to control transmission, rabies prevention remains a challenge in many developing countries, especially in rural areas of China where re-emerging rabies is under-reported due to a lack of sustained animal surveillance. By taking advantage of detailed genomic and epidemiological data for the re-emerging rabies outbreak in Yunnan Province, China, collected between 1999 and 2015, we reconstruct the demographic and dispersal history of domestic dog rabies virus (RABV) as well as the dynamics of dog-to-dog and dog-to-human transmission. Phylogeographic analyses reveal a lower diffusion coefficient than previously estimated for dog RABV dissemination in northern Africa. Furthermore, epidemiological analyses reveal transmission rates between dogs, as well as between dogs and humans, lower than estimates for Africa. Finally, we show that reconstructed epidemic history of RABV among dogs and the dynamics of rabid dogs are consistent with the recorded human rabies cases. This work illustrates the benefits of combining phylogeographic and epidemic modelling approaches for uncovering the spatiotemporal dynamics of zoonotic diseases, with both approaches providing estimates of key epidemiological parameters.

A hundred years of rabies in Kenya and the strategy for eliminating dog-mediated rabies by 2030

Background: Rabies causes an estimated 59,000 human deaths annually. In Kenya, rabies was first reported in a dog in 1912, with the first human case reported in 1928. Here we examine retrospective rabies data in Kenya for the period 1912 – 2017 and describe the spatial and temporal patterns of rabies occurrence in the country. Additionally, we detail Kenya’s strategy for the elimination of dog-mediated human rabies by 2030. Methods: Data on submitted samples and confirmed cases in humans, domestic animals and wildlife were obtained from Kenya’s Directorate of Veterinary Services. These data were associated with the geographical regions where the samples originated, and temporal and spatial trends examined. Results: Between 1912 and the mid 1970’s, rabies spread across Kenya gradually, with fewer than 50 cases reported per year and less than half of the 47 counties affected. Following an outbreak in the mid 1970’s, rabies spread rapidly to more than 85% of counties, with a 4 fold increase in the percent positivity of samples submitted and number of confirmed rabies cases. Since 1958, 7,584 samples from domestic animals (93%), wildlife (5%), and humans (2%) were tested. Over two-thirds of all rabies cases came from six counties, all in close proximity to veterinary diagnostic laboratories, highlighting a limitation of passive surveillance. Conclusions: Compulsory annual dog vaccinations between 1950’s and the early 1970’s slowed rabies spread. The rapid spread with peak rabies cases in the 1980’s coincided with implementation of structural adjustment programs privatizing the veterinary sector leading to breakdown of rabies control programs. To eliminate human deaths from rabies by 2030, Kenya is implementing a 15-year step-wise strategy based on three pillars: a) mass dog vaccination, b) provision of post-exposure prophylaxis and public awareness and c) improved surveillance for rabies in dogs and humans with prompt responses to rabies outbreaks.

Rabies elimination research: juxtaposing optimism, pragmatism and realism

More than 100 years of research has now been conducted into the prevention, control and elimination of rabies with safe and highly efficacious vaccines developed for use in human and animal populations. Domestic dogs are a major reservoir for rabies, and although considerable advances have been made towards the elimination and control of canine rabies in many parts of the world, the disease continues to kill tens of thousands of people every year in Africa and Asia. Policy efforts are now being directed towards a global target of zero human deaths from dog-mediated rabies by 2030 and the global elimination of canine rabies. Here we demonstrate how research provides a cause for optimism as to the feasibility of these goals through strategies based around mass dog vaccination. We summarize some of the pragmatic insights generated from rabies epidemiology and dog ecology research that can improve the design of dog vaccination strategies in low- and middle-income countries and which should encourage implementation without further delay. We also highlight the need for realism in reaching the feasible, although technically more difficult and longer-term goal of global elimination of canine rabies. Finally, we discuss how research on rabies has broader relevance to the control and elimination of a suite of diseases of current concern to human and animal health, providing an exemplar of the value of a ‘One Health’ approach.

Geography but not alternative host species explain the spread of raccoon rabies virus in Vermont

In North America, the raccoon-associated variant of rabies virus (RRV) is of special concern, given its relatively rapid spread throughout the eastern USA and its potential public health impact due to high raccoon host densities in urban areas. Northward expansion of this epizootic included an outbreak in the Canadian province of Quebec in 2006-2009 due to trans-border spread from the State of Vermont. To inform a more proactive approach to future control efforts, this study uses phylogenetic analyses to explore the role of geography and alternative carnivore hosts in the dynamics of RRV spread within Vermont. Specifically, we sought to examine whether striped skunks, a species frequently infected by RRV, could be part of the maintenance host community. Whole genome sequencing of 160 RRV samples from Vermont and neighbouring US states were used for fine-scale phylogeographic analyses. Results, together with the complete surveillance record of raccoon rabies since its entry into Vermont in 1994, document incursions by two distinct viral lineages and identify topographical features of the landscape which have significantly influenced viral spread, resulting in a complex distribution pattern of viral variants throughout the state. Results of phylogenetic cluster analysis and discrete state reconstruction contained some evidence of skunk-to-skunk and skunk-to-raccoon transmission but overall failed to support a role for skunks as alternative maintenance hosts.

The spread and evolution of rabies virus: conquering new frontiers

Rabies is a lethal zoonotic disease that is caused by lyssaviruses, most often rabies virus. Despite control efforts, sporadic outbreaks in wildlife populations are largely unpredictable, underscoring our incomplete knowledge of what governs viral transmission and spread in reservoir hosts. Furthermore, the evolutionary history of rabies virus and related lyssaviruses remains largely unclear. Robust surveillance efforts combined with diagnostics and disease modelling are now providing insights into the epidemiology and evolution of rabies virus. The immune status of the host, the nature of exposure and strain differences all clearly influence infection and transmission dynamics. In this Review, we focus on rabies virus infections in the wildlife and synthesize current knowledge in the rapidly advancing fields of rabies virus epidemiology and evolution, and advocate for multidisciplinary approaches to advance our understanding of this disease.

Recombination in the rabies virus and other lyssaviruses

Recombination is a common event in RNA viruses; however, in the rabies virus there have been only a few reports of isolated recombination events. Comprehensive analysis found traces of recent recombination events within Arctic, Arctic-like and Africa 1b rabies virus groups, as well as recombination between distinct lyssaviruses. Recombination breakpoints were not linked to gene boundaries and could be detected all over the genome. However, there was no evidence that recombination is an important factor in the genetic variability of the rabies virus. It is therefore likely that recombination in the rabies virus is limited by ecological factors (e.g., rare co-circulation of distinguishable lineages and a narrow window for productive coinfection in most carnivore hosts), rather than molecular barriers (e.g., incompatibility of genome fragments).

Landscape attributes governing local transmission of an endemic zoonosis: Rabies virus in domestic dogs

Landscape heterogeneity plays an important role in disease spread and persistence, but quantifying landscape influences and their scale dependence is challenging. Studies have focused on how environmental features or global transport networks influence pathogen invasion and spread, but their influence on local transmission dynamics that underpin the persistence of endemic diseases remains unexplored. Bayesian phylogeographic frameworks that incorporate spatial heterogeneities are promising tools for analysing linked epidemiological, environmental and genetic data. Here, we extend these methodological approaches to decipher the relative contribution and scale-dependent effects of landscape influences on the transmission of endemic rabies virus in Serengeti district, Tanzania (area ~4,900 km2 ). Utilizing detailed epidemiological data and 152 complete viral genomes collected between 2004 and 2013, we show that the localized presence of dogs but not their density is the most important determinant of diffusion, implying that culling will be ineffective for rabies control. Rivers and roads acted as barriers and facilitators to viral spread, respectively, and vaccination impeded diffusion despite variable annual coverage. Notably, we found that landscape effects were scale-dependent: rivers were barriers and roads facilitators on larger scales, whereas the distribution of dogs was important for rabies dispersal across multiple scales. This nuanced understanding of the spatial processes that underpin rabies transmission can be exploited for targeted control at the scale where it will have the greatest impact. Moreover, this research demonstrates how current phylogeographic frameworks can be adapted to improve our understanding of endemic disease dynamics at different spatial scales.

Defining objective clusters for rabies virus sequences using affinity propagation clustering

Rabies is caused by lyssaviruses, and is one of the oldest known zoonoses. In recent years, more than 21,000 nucleotide sequences of rabies viruses (RABV), from the prototype species rabies lyssavirus, have been deposited in public databases. Subsequent phylogenetic analyses in combination with metadata suggest geographic distributions of RABV. However, these analyses somewhat experience technical difficulties in defining verifiable criteria for cluster allocations in phylogenetic trees inviting for a more rational approach. Therefore, we applied a relatively new mathematical clustering algorythm named ‘affinity propagation clustering’ (AP) to propose a standardized sub-species classification utilizing full-genome RABV sequences. Because AP has the advantage that it is computationally fast and works for any meaningful measure of similarity between data samples, it has previously been applied successfully in bioinformatics, for analysis of microarray and gene expression data, however, cluster analysis of sequences is still in its infancy. Existing (516) and original (46) full genome RABV sequences were used to demonstrate the application of AP for RABV clustering. On a global scale, AP proposed four clusters, i.e. New World cluster, Arctic/Arctic-like, Cosmopolitan, and Asian as previously assigned by phylogenetic studies. By combining AP with established phylogenetic analyses, it is possible to resolve phylogenetic relationships between verifiably determined clusters and sequences. This workflow will be useful in confirming cluster distributions in a uniform transparent manner, not only for RABV, but also for other comparative sequence analyses.

Rabies is one of the oldest known zoonoses, caused by lyssaviruses. In recent years, more than 21,000 nucleotide sequences for rabies viruses (RABV) have been deposited in public databases. In this study, a novel mathematical approach called affinity propagation (AP) clustering, a highly powerful tool, to verifiably divide full genome RABV sequences into genetic clusters, was used. A panel of existing and novel RABV full genome sequences was used to demonstrate the application of AP for RABV clustering. Using a combination of AP with established phylogenetic analyses is useful in resolving phylogenetic relationships between more objectively determined clusters and sequences. This workflow will help to substantiate a transparent cluster distribution, not only for RABV, but also for other comparative sequence analyses.

Molecular Epidemiology and Evolution of European Bat Lyssavirus 2

Bat rabies cases in Europe are mainly attributed to two lyssaviruses, namely European Bat Lyssavirus 1 (EBLV-1) and European Bat Lyssavirus 2 (EBLV-2). Prior to the death of a bat worker in Finland in 1985, very few bat rabies cases were reported. Enhanced surveillance in the two subsequent years (1986-1987) identified 263 cases (more than a fifth of all reported cases to date). Between 1977 and 2016, 1183 cases of bat rabies were reported, with the vast majority (>97%) being attributed to EBLV-1. In contrast, there have been only 39 suspected cases of EBLV-2, of which 34 have been confirmed by virus typing and presently restricted to just two bat species; Myotis daubentonii and Myotis dasycneme. The limited number of EBLV-2 cases in Europe prompted the establishment of a network of European reference laboratories to collate all available viruses and data. Despite the relatively low number of EBLV-2 cases, a large amount of anomalous data has been published in the scientific literature, which we have here reviewed and clarified. In this review, 29 EBLV-2 full genome sequences have been analysed to further our understanding of the diversity and molecular evolution of EBLV-2 in Europe. Analysis of the 29 complete EBLV-2 genome sequences clearly corroborated geographical relationships with all EBLV-2 sequences clustering at the country level irrespective of the gene studied. Further geographical clustering was also observed at a local level. There are high levels of homogeneity within the EBLV-2 species with nucleotide identities ranging from 95.5-100% and amino acid identities between 98.7% and 100%, despite the widespread distribution of the isolates both geographically and chronologically. The mean substitution rate for EBLV-2 across the five concatenated genes was 1.65 × 10^-5, and evolutionary clock analysis confirms the slow evolution of EBLV-2 both between and within countries in Europe. This is further supported by the first detailed EBLV-2 intra-roost genomic analysis whereby a relatively high sequence homogeneity was found across the genomes of three EBLV-2 isolates obtained several years apart (2007, 2008, and 2014) from M. daubentonii at the same site (Stokesay Castle, Shropshire, UK).

Rabies

Rabies is a life-threatening neglected tropical disease: tens of thousands of cases are reported annually in endemic countries (mainly in Africa and Asia), although the actual numbers are most likely underestimated. Rabies is a zoonotic disease that is caused by infection with viruses of the Lyssavirus genus, which are transmitted via the saliva of an infected animal. Dogs are the most important reservoir for rabies viruses, and dog bites account for >99% of human cases. The virus first infects peripheral motor neurons, and symptoms occur after the virus reaches the central nervous system. Once clinical disease develops, it is almost certainly fatal. Primary prevention involves dog vaccination campaigns to reduce the virus reservoir. If exposure occurs, timely post-exposure prophylaxis can prevent the progression to clinical disease and involves appropriate wound care, the administration of rabies immunoglobulin and vaccination. A multifaceted approach for human rabies eradication that involves government support, disease awareness, vaccination of at-risk human populations and, most importantly, dog rabies control is necessary to achieve the WHO goal of reducing the number of cases of dog-mediated human rabies to zero by 2030.

Molecular characterization and phylogenetic relatedness of dog-derived Rabies Viruses circulating in Cameroon between 2010 and 2016

Rabies is enzootic among dog populations in some parts of Cameroon and the risk of human rabies is thought to be steadily high in these regions. However, the molecular epidemiology of circulating Rabies Virus (RABV) has been hardly considered in Cameroon as well as in most neighboring central African countries. To address this fundamental gap, 76 nucleoprotein (N) gene sequences of dog-derived RABV were obtained from 100 brain specimens sampled in Cameroon from 2010 to 2016. Studied sequences were subjected to molecular and phylogenetic analyses with reference strains retrieved from databases. The 71 studied Africa-1 isolates displayed 93.5–100% nucleotide (nt) and 98.3–100% amino-acid (aa) identities to each other while, the 5 studied Africa-2 isolates shared 99.4–99.7% sequence similarities at nt and aa levels. Maximum Likelihood based phylogenies inferred from nucleotide sequences confirmed all studied RABV isolates as members of the dog-related species 1 of the Lyssavirus genus. Individual isolates could be unambiguously assigned as either the Africa-1 subclade of the Cosmopolitan clade or the Africa 2 clade. The Africa-1 subclade appeared to be more prevalent and diversified. Indeed, 70 studied isolates segregated into 3 distinct circulating variants within Africa-1a lineage while a unique isolate was strikingly related to the Africa-1b lineage known to be prevalent in the neighboring Central African Republic and eastern Africa. Interestingly, all five Africa-2 isolates fell into the group-E lineage even though they appeared to be loosely related to databases available reference RABV; including those previously documented in Cameroon. This study uncovered the co-circulation of several Africa-1 and Africa-2 lineages in the southern regions of Cameroon. Striking phylogenetic outcasts to the geographic differentiation of RABV variants indicated that importation from close regions or neighboring countries apparently contributes to the sustainment of the enzootic cycle of domestic rabies in Cameroon.

Rabies has been repeatedly reported among dog populations in Cameroon, especially in Yaounde, its capital city. However, the relative rates and genetic variability of Rabies Virus (RABV) variants circulating among dog populations in Cameroon are still to be documented. This study aimed to estimate the frequency and genetic diversity of RABV isolates originating from rabid dogs in the southern regions of Cameroon from 2010 to 2016. Overall, 76 of the 100 dog-derived RABV isolates sampled in Cameroon from 2010 to 2016 were successfully characterized. Our findings revealed that studied isolates belonged to the dog-related species 1 of the Lyssavirus genus, specifically 70 Africa-1a, 1 Africa-1b and 5 Africa-2 group-E lineages. The general phylogenetic pattern suggested an in-country geographic differentiation of the circulating RABV variants. This apparent geographic differentiation was contradicted by striking outcasts indicating importation from close or distant regions. Overall, this study uncovered the co-circulation of several Africa-1 and Africa-2 lineages in some southern regions of Cameroon, thus providing base-line molecular data that would be of interest for future stages of implementation of the rabies surveillance and control plan that is being setup in Cameroon.

Exploiting Genetic Information to Trace Plant Virus Dispersal in Landscapes

During the past decade, knowledge of pathogen life history has greatly benefited from the advent and development of molecular epidemiology. This branch of epidemiology uses information on pathogen variation at the molecular level to gain insights into a pathogen's niche and evolution and to characterize pathogen dispersal within and between host populations. Here, we review molecular epidemiology approaches that have been developed to trace plant virus dispersal in landscapes. In particular, we highlight how virus molecular epidemiology, nourished with powerful sequencing technologies, can provide novel insights at the crossroads between the blooming fields of landscape genetics, phylogeography, and evolutionary epidemiology. We present existing approaches and their limitations and contributions to the understanding of plant virus epidemiology.

First detection of European bat lyssavirus type 2 (EBLV-2) in Norway

BACKGROUND

In Europe, bat rabies is primarily attributed to European bat lyssavirus type 1 (EBLV-1) and European bat lyssavirus type 2 (EBLV-2) which are both strongly host-specific. Approximately thirty cases of infection with EBLV-2 in Daubenton's bats (Myotis daubentonii) and pond bats (M. dasycneme) have been reported. Two human cases of rabies caused by EBLV-2 have also been confirmed during the last thirty years, while natural spill-over to other non-flying mammals has never been reported. Rabies has never been diagnosed in mainland Norway previously.

CASE PRESENTATION

In late September 2015, a subadult male Daubenton's bat was found in a poor condition 800 m above sea level in the southern part of Norway. The bat was brought to the national Bat Care Centre where it eventually displayed signs of neurological disease and died after two days. EBLV-2 was detected in brain tissues by polymerase chain reaction (PCR) followed by sequencing of a part of the nucleoprotein gene, and lyssavirus was isolated in neuroblastoma cells.

CONCLUSIONS

The detection of EBLV-2 in a bat in Norway broadens the knowledge on the occurrence of this zoonotic agent. Since Norway is considered free of rabies, adequate information to the general public regarding the possibility of human cases of bat-associated rabies should be given. No extensive surveillance of lyssavirus infections in bats has been conducted in the country, and a passive surveillance network to assess rabies prevalence and bat epidemiology is highly desired.

Application of high-throughput sequencing to whole rabies viral genome characterisation and its use for phylogenetic re-evaluation of a raccoon strain incursion into the province of Ontario

Raccoon rabies remains a serious public health problem throughout much of the eastern seaboard of North America due to the urban nature of the reservoir host and the many challenges inherent in multi-jurisdictional efforts to administer co-ordinated and comprehensive wildlife rabies control programmes. Better understanding of the mechanisms of spread of rabies virus can play a significant role in guiding such control efforts. To facilitate a detailed molecular epidemiological study of raccoon rabies virus movements across eastern North America, we developed a methodology to efficiently determine whole genome sequences of hundreds of viral samples. The workflow combines the generation of a limited number of overlapping amplicons covering the complete viral genome and use of high throughput sequencing technology. The value of this approach is demonstrated through a retrospective phylogenetic analysis of an outbreak of raccoon rabies which occurred in the province of Ontario between 1999 and 2005. As demonstrated by the number of single nucleotide polymorphisms detected, whole genome sequence data were far more effective than single gene sequences in discriminating between samples and this facilitated the generation of more robust and informative phylogenies that yielded insights into the spatio-temporal pattern of viral spread. With minor modification this approach could be applied to other rabies virus variants thereby facilitating greatly improved phylogenetic inference and thus better understanding of the spread of this serious zoonotic disease. Such information will inform the most appropriate strategies for rabies control in wildlife reservoirs.

Lyssaviruses and rabies: current conundrums, concerns, contradictions and controversies

Lyssaviruses are bullet-shaped, single-stranded, negative-sense RNA viruses and the causative agents of the ancient zoonosis rabies. Africa is the likely home to the ancestors of taxa residing within the Genus Lyssavirus, Family Rhabdoviridae. Diverse lyssaviruses are envisioned as co-evolving with bats, as the ultimate reservoirs, over seemingly millions of years. In terms of relative distribution, overt abundance, and resulting progeny, rabies virus is the most successful lyssavirus species today, but for unknown reasons. All mammals are believed to be susceptible to rabies virus infection. Besides reservoirs among the Chiroptera, meso-carnivores also serve as major historical hosts and are represented among the canids, raccoons, skunks, mongooses, and ferret badgers.  Perpetuating as a disease of nature with the mammalian central nervous system as niche, host breadth alone precludes any candidacy for true eradication. Despite having the highest case fatality of any infectious disease and a burden in excess of or comparative to other major zoonoses, rabies remains neglected. Once illness appears, no treatment is proven to prevent death. Paradoxically, vaccines were developed more than a century ago, but the clear majority of human cases are unvaccinated. Tens of millions of people are exposed to suspect rabid animals and tens of thousands succumb annually, primarily children in developing countries, where canine rabies is enzootic. Rather than culling animal populations, one of the most cost-effective strategies to curbing human fatalities is the mass vaccination of dogs. Building on considerable progress to date, several complementary actions are needed in the near future, including a more harmonized approach to viral taxonomy, enhanced de-centralized laboratory-based surveillance, focal pathogen discovery and characterization, applied pathobiological research for therapeutics, improved estimates of canine populations at risk, actual production of required vaccines and related biologics, strategies to maximize prevention but minimize unnecessary human prophylaxis, and a long-term, realistic plan for sustained global program support to achieve success in disease control, prevention, and elimination.

Large-Scale Phylogenomic Analysis Reveals the Complex Evolutionary History of Rabies Virus in Multiple Carnivore Hosts

The natural evolution of rabies virus (RABV) provides a potent example of multiple host shifts and an important opportunity to determine the mechanisms that underpin viral emergence. Using 321 genome sequences spanning an unprecedented diversity of RABV, we compared evolutionary rates and selection pressures in viruses sampled from multiple primary host shifts that occurred on various continents. Two major phylogenetic groups, bat-related RABV and dog-related RABV, experiencing markedly different evolutionary dynamics were identified. While no correlation between time and genetic divergence was found in bat-related RABV, the evolution of dog-related RABV followed a generally clock-like structure, although with a relatively low evolutionary rate. Subsequent molecular clock dating indicated that dog-related RABV likely underwent a rapid global spread following the intensification of intercontinental trade starting in the 15^th century. Strikingly, although dog RABV has jumped to various wildlife species from the order Carnivora, we found no clear evidence that these host-jumping events involved adaptive evolution, with RABV instead characterized by strong purifying selection, suggesting that ecological processes also play an important role in shaping patterns of emergence. However, specific amino acid changes were associated with the parallel emergence of RABV in ferret-badgers in Asia, and some host shifts were associated with increases in evolutionary rate, particularly in the ferret-badger and mongoose, implying that changes in host species can have important impacts on evolutionary dynamics.

Zoonoses account for most recently emerged infectious diseases of humans, although little is known about the evolutionary mechanisms involved in cross-species virus transmission. Understanding the evolutionary patterns and processes that underpin such cross-species transmission is of importance for predicting the spread of zoonotic infections, and hence to their ultimate control. We present a large-scale and detailed reconstruction of the evolutionary history of rabies virus (RABV) in domestic and wildlife animal species. RABV is of particular interest as it is capable of infecting many mammals but, paradoxically, is only maintained in distinct epidemiological cycles associated with animal species from the orders Carnivora and Chiroptera. We show that bat-related RABV and dog-related RABV have experienced very different evolutionary dynamics, and that host jumps are sometimes characterized by significant increases in evolutionary rate. Among Carnivora, the association between RABV and particular host species most likely arose from a combination of the historical human-mediated spread of the virus and jumps into new primary host species. In addition, we show that changes in host species are associated with multiple evolutionary pathways including the occurrence of host-specific parallel evolution. Overall, our data indicate that the establishment of dog-related RABV in new carnivore hosts may only require subtle adaptive evolution.