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.

Whole Genome Sequencing and Phylogenetic Analysis of Rabies Viruses from Bats in Connecticut, USA, 2018-2019

We performed whole genome sequencing and genetic characterization of rabies viruses (RABV) detected in bats submitted to the Connecticut Veterinary Medical Diagnostic Laboratory (CVMDL) during 2018–2019. Among 88 bats submitted to CVMDL, six brain samples (6.8%, 95% confidence interval: 1.6% to 12.1%) tested positive by direct fluorescent antibody test. RABVs were detected in big brown bats (Eptesicus fuscus, n = 4), a hoary bat (Lasiurus cinereus, n = 1), and an unidentified bat species (n = 1). Complete coding sequences of four out of six detected RABVs were obtained. In phylogenetic analysis, the RABVs (18-62, 18-4347, and 19-2274) from big brown bats belong to the bats EF-E1 clade, clustering with RABVs detected from the same bat species in Pennsylvania and New Jersey. The bat RABV (19-2898) detected from the migratory hoary bat belongs to the bats LC clade, clustering with the eleven viruses detected from the same species in Arizona, Washington, Idaho, and Tennessee. The approach used in this study generated novel data regarding genetic relationships of RABV variants, including their reservoirs, and their spatial origin and it would be useful as reference data for future investigations on RABV in North America. Continued surveillance and genome sequencing of bat RABV would be needed to monitor virus evolution and transmission, and to assess the emergence of genetic mutations that may be relevant for public health.

A Novel Terrestrial Rabies Virus Lineage Occurring in South America: Origin, Diversification, and Evidence of Contact between Wild and Domestic Cycles

The rabies virus (RABV) is characterized by a history dominated by host shifts within and among bats and carnivores. One of the main outcomes of long-term RABV maintenance in dogs was the establishment of variants in a wide variety of mesocarnivores. In this study, we present the most comprehensive phylogenetic and phylogeographic analysis, contributing to a better understanding of the origins, diversification, and the role of different host species in the evolution and diffusion of a dog-related variant endemic of South America. A total of 237 complete Nucleoprotein gene sequences were studied, corresponding to wild and domestic species, performing selection analyses, ancestral states reconstructions, and recombination analyses. This variant originated in Brazil and disseminated through Argentina and Paraguay, where a previously unknown lineage was found. A single host shift was identified in the phylogeny, from dog to the crab-eating fox (Cerdocyon thous) in the Northeast of Brazil. Although this process occurred in a background of purifying selection, there is evidence of adaptive evolution -or selection of sub-consensus sequences- in internal branches after the host shift. The interaction of domestic and wild cycles persisted after host switching, as revealed by spillover and putative recombination events.

Molecular Epidemiology of Rabies in Wild Canidae in Tunisia

Rabies is a viral zoonosis that is transmissible to humans via domestic and wild animals. There are two epidemiological cycles for rabies, the urban and the sylvatic cycles. In an attempt to study the epidemiological role of wild canidae in rabies transmission, the present study aimed to analyze the genetic characteristics of virus isolates and confirm prior suggestions that rabies is maintained through a dog reservoir in Tunisia. Virus strains isolated from wild canidae were subject to viral sequencing, and Bayesian phylogenetic analysis was performed using Beast2 software. Essentially, the virus strains isolated from wild canidae belonged to the Africa-1 clade, which clearly diverges from fox-related strains. Our study also demonstrated that genetic characteristics of the virus isolates were not as distinct as could be expected if a wild reservoir had already existed. On the contrary, the geographic landscape is responsible for the genetic diversity of the virus. The landscape itself could have also acted as a natural barrier to the spread of the virus.

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.

Determination of freedom-from-rabies for small Indian mongoose populations in the United States Virgin Islands, 2019-2020

Mongooses, a nonnative species, are a known reservoir of rabies virus in the Caribbean region. A cross-sectional study of mongooses at 41 field sites on the US Virgin Islands of St. Croix, St. John, and St. Thomas captured 312 mongooses (32% capture rate). We determined the absence of rabies virus by antigen testing and rabies virus exposure by antibody testing in mongoose populations on all three islands. USVI is the first Caribbean state to determine freedom-from-rabies for its mongoose populations with a scientifically-led robust cross-sectional study. Ongoing surveillance activities will determine if other domestic and wildlife populations in USVI are rabies-free.

Mongooses in the Caribbean region are known to carry rabies and infect and kill humans with the deadly virus. While many countries in the Caribbean region assume they are rabies-free, there is a lack of scientifically based surveillance. We prospectively designed a statistically valid study to determine that mongoose populations in the US Virgin Islands are rabies-free and pose no risk of transmission to human or animal populations. We are pursuing further surveillance efforts to establish rabies freedom in all wildlife and domestic animal population in the US Virgin Islands.

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.

Modeling Mongoose Rabies in the Caribbean: A Model-Guided Fieldwork Approach to Identify Research Priorities

We applied the model-guided fieldwork framework to the Caribbean mongoose rabies system by parametrizing a spatially-explicit, individual-based model, and by performing an uncertainty analysis designed to identify parameters for which additional empirical data are most needed. Our analysis revealed important variation in output variables characterizing rabies dynamics, namely rabies persistence, exposure level, spatiotemporal distribution, and prevalence. Among epidemiological parameters, rabies transmission rate was the most influential, followed by rabies mortality and location, and size of the initial infection. The most influential landscape parameters included habitat-specific carrying capacities, landscape heterogeneity, and the level of resistance to dispersal associated with topography. Movement variables, including juvenile dispersal, adult fine-scale movement distances, and home range size, as well as life history traits such as age of independence, birth seasonality, and age- and sex-specific mortality were other important drivers of rabies dynamics. We discuss results in the context of mongoose ecology and its influence on disease transmission dynamics. Finally, we suggest empirical approaches and study design specificities that would provide optimal contributing data addressing the knowledge gaps identified by our approach, and would increase our potential to use epidemiological models to guide mongoose rabies control and management in the Caribbean.

Characterizing and Evaluating the Zoonotic Potential of Novel Viruses Discovered in Vampire Bats

The contemporary surge in metagenomic sequencing has transformed knowledge of viral diversity in wildlife. However, evaluating which newly discovered viruses pose sufficient risk of infecting humans to merit detailed laboratory characterization and surveillance remains largely speculative. Machine learning algorithms have been developed to address this imbalance by ranking the relative likelihood of human infection based on viral genome sequences, but are not yet routinely applied to viruses at the time of their discovery. Here, we characterized viral genomes detected through metagenomic sequencing of feces and saliva from common vampire bats (Desmodus rotundus) and used these data as a case study in evaluating zoonotic potential using molecular sequencing data. Of 58 detected viral families, including 17 which infect mammals, the only known zoonosis detected was rabies virus; however, additional genomes were detected from the families Hepeviridae, Coronaviridae, Reoviridae, Astroviridae and Picornaviridae, all of which contain human-infecting species. In phylogenetic analyses, novel vampire bat viruses most frequently grouped with other bat viruses that are not currently known to infect humans. In agreement, machine learning models built from only phylogenetic information ranked all novel viruses similarly, yielding little insight into zoonotic potential. In contrast, genome composition-based machine learning models estimated different levels of zoonotic potential, even for closely related viruses, categorizing one out of four detected hepeviruses and two out of three picornaviruses as having high priority for further research. We highlight the value of evaluating zoonotic potential beyond ad hoc consideration of phylogeny and provide surveillance recommendations for novel viruses in a wildlife host which has frequent contact with humans and domestic animals.

Development of Monoclonal Antibodies for Detection of Conserved and Variable Epitopes of Large Protein of Rabies Virus

Rabies virus (RABV) causes fatal neurological encephalitis and results in approximately 6000 human death cases worldwide every year. The large (L) protein of RABV, possessing conserved domains, is considered as the target for detection. In this study, three monoclonal antibodies (mAbs), designated as 3F3, 3A6 and L-C, against L protein were generated by using the recombinant truncated L protein (aa 1431-1754) and the epitopes were also identified using a series of overlapping truncated polypeptides for testing the reactivity of mAbs with different RABV strains. The ¹⁴⁷⁹EIFSIP¹⁴⁸⁴, ¹⁶⁵⁹RALSK¹⁶⁶³ and ¹⁷²⁴VFNSL¹⁷²⁸ were identified as the minimal linear epitopes recognized by mAbs 3F3, 3A6 and L-C, respectively. Amino acid alignment showed epitope ¹⁷²⁴VFNSL¹⁷²⁸ recognized by mAb L-C is completely conserved among RABV strains, indicating that mAb L-C could be used to detect all of the RABV strains. Epitope ¹⁴⁷⁹EIFSIP¹⁴⁸⁴ is highly conserved among RABV strains except for a P1484S substitution in a China I sub-lineage strain of Asian lineage, which eliminated the reactivity of the epitope with mAb 3F3. However, the epitope ¹⁶⁵⁹RALSK¹⁶⁶³ was only completely conserved in the Africa-2 and Indian lineages, and a single A1660T substitution, mainly appeared in strains of the China I belonging to Asian lineage and a Cosmopolitan lineage strain, still retained the reactivity of the epitope with mAb 3A6. While both A1660T and K1663R substitutions in a China I lineage strain, single K1663R/Q substitution in some China II strains of Asian lineage and some Arctic-like lineage strains and R1659Q mutation in a strain of Africa-3 lineage eliminated the reactivity of the epitope with mAb 3A6, suggesting mAb 3A6 could be used for differentiation of variable epitopes of some strains in different lineages. Thus, variability and conservation of the three epitopes of L protein showed the reactive difference of mAbs among RABV strains of different lineages. These results may facilitate future studies in development of detection methods for RABV infection, the structure and function of RABV L protein.

Evaluation of Two Real-Time, TaqMan Reverse Transcription-PCR Assays for Detection of Rabies Virus in Circulating Variants from Argentina: Influence of Sequence Variation

In rabies diagnosis, it is essential to count on a rapid test to give a quick response. The combined sensitivity and robustness of the TaqMan RT-PCR assays (qRT-PCR) have made these methods a valuable alternative for rabies virus (RABV) detection. We conducted a study to compare the applicability of two widely used qRT-PCR assays targeting the nucleoprotein gene (LysGT1 assay) and leader sequences (LN34 qRT-PCR assay) of RABV genomes, in all variants circulating in Argentina. A total of 44 samples obtained from bats, dogs, cattle, and horses, that were previously tested for rabies by FAT and conventional RT-PCR, were used in the study. All variants were successfully detected by the pan-lyssavirus LN34 qRT-PCR assay. The LysGT1 assay failed to detect three bat-related variants. We further sequenced the region targeted by LysGT1 and demonstrated that the presence of three or more mismatches with respect to the primers and probe sequences precludes viral detection. We conclude that the LysGT1 assay is prone to yield variant-dependent false-negative test results, and in consequence, the LN34 assay would ensure more effective detection of RABV in Argentina.

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.

Quantifying and mapping the burden of human and animal rabies in Iraq

Rabies was first reported in ancient Iraqi civilizations, yet it remains a poorly quantified and important public health threat in the region. Efforts to control rabies in Iraq including dog population control, and vaccination of livestock and dogs, have increased since 2010. Officially reported data on human rabies, dog bites, and animal rabies cases between 2012 and 2017 are analysed here to assess the effect of existing control efforts, to inform future strategies, and to highlight gaps in surveillance and reporting. The results of molecular characterization of 32 viruses from animal cases from throughout Iraq are presented, to improve the understanding of rabies dynamics in the animal reservoir. Although annual numbers of reported human cases were lower in the period between 2012 and 2017 than prior to 2010, human cases continue. There was a distinct gender and age bias among human cases with nine cases in males for every one female and twice as many cases in children than adults. Spatial clustering analysis and phylogenetic evidence suggests rabies is endemic throughout the country, with no regional variation in risk, but better surveillance and reporting is required to underpin control strategies.

Genetic and Phenotypic Characterization of a Rabies Virus Strain Isolated from a Dog in Tokyo, Japan in the 1940s

The rabies virus strain Komatsugawa (Koma), which was isolated from a dog in Tokyo in the 1940s before eradication of rabies in Japan in 1957, is known as the only existent Japanese field strain (street strain). Although this strain potentially provides a useful model to study rabies pathogenesis, little is known about its genetic and phenotypic properties. Notably, this strain underwent serial passages in rodents after isolation, indicating the possibility that it may have lost biological characteristics as a street strain. In this study, to evaluate the utility of the Koma strain for studying rabies pathogenesis, we examined the genetic properties and in vitro and in vivo phenotypes. Genome-wide genetic analyses showed that, consistent with previous findings from partial sequence analyses, the Koma strain is closely related to a Russian street strain within the Arctic-related phylogenetic clade. Phenotypic examinations in vitro revealed that the Koma strain and the representative street strains are less neurotropic than the laboratory strains. Examination by using a mouse model demonstrated that the Koma strain and the street strains are more neuroinvasive than the laboratory strains. These findings indicate that the Koma strain retains phenotypes similar to those of street strains, and is therefore useful for studying rabies pathogenesis.

A long-distance translocation initiated an outbreak of raccoon rabies in Hamilton, Ontario, Canada

Despite proactive measures to prevent raccoon rabies entering Canada from the United States, several incursions of this disease have occurred. The largest outbreak, first reported in December 2015 in the city of Hamilton, Ontario, has resulted in the reporting of 449 animal cases as of December 31, 2018. Initial phylogenetic studies on the index case suggested that this outbreak was not due to local cross-border spread from the Niagara region of the United States where raccoon rabies has persisted for several years. Phylogenetic analysis of whole genome sequences of a viral collection from the Hamilton area and several US states indicates that a long-distance translocation of a diseased animal from southeastern New York State was responsible for this incursion. The role of the skunk as a potential secondary host supporting persistence and / or spread of the virus is also examined.

Distinct variants of rabies virus are harbored by several wildlife species across North America. The variant associated with the raccoon host is especially problematic given its rapid spread throughout the eastern seaboard of the United States over the second half of the 20^th century and the high population density of this species in urban habitats that raises the potential for disease spread to both companion animals and their owners. Accordingly, extensive efforts to control raccoon rabies have been undertaken at many jurisdictional levels. This study uses current methods of genome analysis to explore the origins of a raccoon rabies outbreak in and around the city of Hamilton, Ontario, Canada, located close to the border with western New York state. In contrast to previous Canadian outbreaks of raccoon rabies, which were due to cross-border spread of the virus, this study demonstrates that the incursion into the Hamilton area was the result of a long-distance animal translocation from south-eastern New York. The implications of this event for strategies to control raccoon rabies are discussed.

Molecular Detection of Rabies Lyssaviruses from Dogs in Southeastern Nigeria: Evidence of TransboundaryTransmission of Rabies in West Africa

Despite being the first country to register confirmed cases of Mokola and Lagos bat lyssaviruses (two very distant lyssaviruses), knowledge gaps, particularly on the molecular epidemiology of lyssaviruses, still exist in Nigeria. A total of 278 specimens were collected from dogs in southeastern Nigeria between October 2015 and July 2016, and 23 (8.3%) of these tested positive for lyssaviruses with the direct fluorescent antibody test (DFA). The lyssaviruses were genetically characterized by amplifying the highly conserved nucleoprotein (N) gene of the rabies lyssaviruses (RABVs) of the viral genome. Phylogenetic analyses of the nucleotide sequences showed that all the RABV sequences in this study were of the Africa-2 lineage. Our results demonstrated that transboundary transmission of rabies lyssavirus is a key event, given that one of the RABV sequences (MN196576) clustered with rabies variants from neighboring Niger Republic. Furthermore, three RABVs from dogs from Anambra State clustered separately forming a novel and distinct group. Our results demonstrated that transboundary transmission of RABLVs is a key driver in the spread of rabies in West Africa. In order for the successful control of this zoonotic disease, a multinational stepwise surveillance and elimination of rabies in Africa by 2030 is probably the solution for regional elimination.

Population- and Variant-Based Genome Analyses of Viruses from Vaccine-Derived Rabies Cases Demonstrate Product Specific Clusters and Unique Patterns

Rabies in wildlife has been successfully controlled in parts of Europe and North America using oral rabies vaccination, i.e., the distribution of baits containing live-attenuated virus strains. Occasionally, these vaccines caused vaccine virus-induced rabies cases. To elucidate the mechanisms of genetic selection and the effect of viral populations on these rabies cases, a next generation sequencing approach as well as comprehensive data analyses of the genetic diversity of Street Alabama Dufferin (SAD) and ERA vaccine virus strains and vaccine-induced rabies cases from Canada and several European countries were conducted. As a result, twelve newly generated sets of sequencing data from Canada and Poland were added to a pool of previously investigated samples. While the population-based analysis showed a segregation of viruses of ERA vaccine-induced rabies cases from those of SAD Bern original (SAD Bern_orig)-derived rabies cases, the in-depth variant analysis revealed three distinct combinations of selected variants for the ERA vaccine-induced cases, suggesting the presence of multiple replication-competent haplotypes in the investigated ERA-BHK21 vaccine. Our findings demonstrate the potential of a deep sequencing approach in combination with comprehensive analyses on the consensus, population, and variant level.

Origins of the arctic fox variant rabies viruses responsible for recent cases of the disease in southern Ontario

A subpopulation of the arctic fox lineage of rabies virus has circulated extensively in red fox populations of Ontario, Canada, between the 1960s and 1990s. An intensive wildlife rabies control program, in which field operations were initiated in 1989, resulted in elimination of the disease in eastern Ontario. However in southwestern Ontario, as numbers of rabid foxes declined the proportion of skunks confirmed to be infected with this rabies virus variant increased and concerted control efforts targeting this species were employed to eliminate the disease. Since 2012 no cases due to this viral variant were reported in southwestern Ontario until 2015 when a single case of rabies due to the arctic fox variant was reported in a bovine. Several additional cases have been documented subsequently. Since routine antigenic typing cannot discriminate between the variants which previously circulated in Ontario and those from northern Canada it was unknown whether these recent cases were the result of a new introduction of this variant or a continuation of the previous enzootic. To explore the origins of this new outbreak whole genome sequences of a collection of 128 rabies viruses recovered from Ontario between the 1990s to the present were compared with those representative of variants circulating in the Canadian north. Phylogenetic analysis shows that the variant responsible for current cases in southwestern Ontario has evolved from those variants known to circulate in Ontario previously and is not due to a new introduction from northern regions. Thus despite ongoing passive surveillance the persistence of wildlife rabies went undetected in the study area for almost three years. The apparent adaptation of this rabies virus variant to the skunk host provided the opportunity to explore coding changes in the viral genome which might be associated with this host shift. Several such changes were identified including a subset for which the operation of positive selection was supported. The location of a small number of these amino acid substitutions in or close to protein motifs of functional importance suggests that some of them may have played a role in this host shift.

Rabies, a serious disease which almost invariably results in death once clinical signs appear, is caused by the rabies virus. The arctic fox lineage of rabies virus persists in fox populations in northern Canada and from the 1960s onwards a sub-population of this virus has circulated in red foxes in the province of Ontario, Canada. By 2000 a provincial wildlife control program attempting to control fox rabies in southern Ontario had eliminated the disease in most of the targeted regions except for a focus of rabies in certain counties in southwestern Ontario which involved mostly skunks. Modified control efforts resulted in apparent rabies elimination with no cases reported in southern Ontario after May 2012 until a rabid bovine was identified in December 2015. Additional cases, mostly in skunks and livestock, have been reported subsequently. To discriminate between a re-introduction of the virus from the north and re-emergence of the Ontario viral variant, we sequenced and compared the genome of multiple viruses from both areas. The sequence data indicate that sub-types of the Ontario variant which eluded control efforts are responsible for these recent cases. These data also permit an examination of potential genetic changes in the virus that may facilitate its circulation in this new wildlife host, the skunk.

Rapid identification of the raccoon rabies virus variant using a real-time reverse-transcriptase polymerase chain reaction

The raccoon-associated variant of rabies virus (RRV) is enzootic throughout the eastern seaboard of the United States with frequent incursions into Canada. Many wildlife management agencies are actively engaged in control programmes targeting elimination of this disease and rapid identification of raccoon rabies cases is crucial to the success of these operations. This report documents the development of a reverse transcriptase real-time PCR (RT-qPCR) that specifically identifies this rabies virus variant (RRV RT-qPCR) and which can be readily multiplexed with a generic rabies virus RT-qPCR for use as a typing tool. Using a large collection of rabies virus samples representative of the variants circulating around the world, but with a focus on those occurring in the Americas, the RRV RT-qPCR was 100% sensitive and 99.31% specific. To further apply these assays for diagnostic purposes, addition of an RT-qPCR targeting the host β-actin mRNA, which serves as an internal amplification control, in a triplex format was shown to yield highly comparable results using a subset of our viral collection. Use of these assays for early and accurate identification of this viral variant will help to optimize the utilization of resources required for control of this disease.

The reemergence of human rabies and emergence of an Indian subcontinent lineage in Tibet, China

Coordinated surveillance, vaccination and public information efforts have brought the Chinese rabies epizootic under control, but significant numbers of fatalities are still reported annually with some cases occurring in previously rabies free regions. Tibet has remained virtually rabies free for 16 years, but since 2015 one human rabies case has been reported each year. To better understand the origins of these cases, we sequenced three human samples and an additional sample isolated from a dog in 2012. Three genomes were sequenced from brain samples: human case 1 (reported in 2015), human case 3 (2017), and the 2012 dog case. For human case 2 (2016), the rabies N gene was sequenced from a limited saliva sample. Phylogenetic analysis shows that Case 1 (CXZ1501H) and the dog case (CXZ1201D) belong to China IV lineage (equivalent to Arctic-like-2 in global rabies), suggesting an association with a wildlife spillover event. However, Case 2 (CXZ1601H) is placed within the dominant lineage China I, and was most similar with recent strains from neighboring Yunnan province, indicating the current epizootic has finally reached Tibet. Most surprisingly however, was the finding that Case 3 (CXZ1704H) is distinct from other Chinese isolates. This isolate is placed in the Indian Subcontinent clade, similar to recent Nepal strains, indicating that cross-border transmission is a new source for rabies infections. Thus, the complex mixture of the rabies epizootic in Tibet represents a major new challenge for Tibet and national rabies control.

Processes Underlying Rabies Virus Incursions across US-Canada Border as Revealed by Whole-Genome Phylogeography

Disease control programs aim to constrain and reduce the spread of infection. Human disease interventions such as wildlife vaccination play a major role in determining the limits of a pathogen’s spatial distribution. Over the past few decades, a raccoon-specific variant of rabies virus (RRV) has invaded large areas of eastern North America. Although expansion into Canada has been largely prevented through vaccination along the US border, several outbreaks have occurred in Canada. Applying phylogeographic approaches to 289 RRV whole-genome sequences derived from isolates collected in Canada and adjacent US states, we examined the processes underlying these outbreaks. RRV incursions were attributable predominantly to systematic virus leakage of local strains across areas along the border where vaccination has been conducted but also to single stochastic events such as long-distance translocations. These results demonstrate the utility of phylogeographic analysis of pathogen genomes for understanding transboundary outbreaks.

Evaluation and Phylogenetic Analysis of Regular Rabies Virus Vaccine Strains

BACKGROUND

This study aimed to evaluate Rabies virus vaccine strains. The obtained results may be helpful for vaccine producers and researchers to compare the strains with wild type and other vaccine strains and select the correct strain to challenge their products.

METHODS

Fourteen rabies virus vaccine strains were compared with each other. The full genomes of the selected strains were taken from the GenBank and the N, P and G genes were labeled. The major and minor antigenic sites of these sequences were identified and contrasted with each other. The identity matrix was designed for rabies virus full genome, N and G genes. In addition, the phylogenetic tree was drawn based on rabies virus N gene for deep analysis.

RESULTS

Although there were no significant differences between antigenic sites in N, P, and G genes, there were noticeable differences for full genome identity matrix and this significant difference can also be observed in N and G identity matrix. In the phylogenetic tree, the Iranian sequences were distant from currently applied vaccine strains.

CONCLUSION

It is necessary to pay attention to the results shown in phylogenetic tree because they warn us about distance between the Iranian sequences and current strains used in applied vaccines. In addition, the obtained results help vaccine producers to choose a correct strain to challenge their product and evaluate their vaccine potency.

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 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.

The history of rabies in the Western Hemisphere

Before the introduction of control programs in the 20th century, rabies in domestic dogs occurred throughout the Western Hemisphere. However, historical records and phylogenetic analysis of multiple virus isolates indicate that, before the arrival of the first European colonizers, rabies virus was likely present only in bats and skunks. Canine rabies was either rare or absent among domestic dogs of Native Americans, and first arrived when many new dog breeds were imported during the period of European colonization. The introduction of the cosmopolitan dog rabies lyssavirus variant and the marked expansion of the dog population provided ideal conditions for the flourishing of enzootic canine rabies. The shift of dog-maintained viruses into gray foxes, coyotes, skunks and other wild mesocarnivores throughout the Americas and to mongooses in the Caribbean has augmented the risk of human rabies exposures and has complicated control efforts. At the same time, the continued presence of bat rabies poses novel challenges in the absolute elimination of canine and human rabies. This article compiles existing historical and phylogenetic evidence of the origins and subsequent dynamics of rabies in the Western Hemisphere, from the era preceding the arrival of the first European colonizers through the present day. A companion article reviews the current status of canine rabies control throughout the Western Hemisphere and steps that will be required to achieve and maintain its complete elimination (Velasco-Villa et al., 2017).

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.

Molecular Characterization of Canine Rabies Virus, Mali, 2006-2013

We genetically characterized 32 canine rabies viruses isolated in Mali during 2006–2013 and identified 3 subgroups that belonged to the Africa 2 lineage. We also detected subgroup F rabies virus. This information should be useful for development of mass vaccination campaigns for dogs and eventual large-scale control programs in this country.

Spatio-temporal Analysis of the Genetic Diversity of Arctic Rabies Viruses and Their Reservoir Hosts in Greenland

There has been limited knowledge on spatio-temporal epidemiology of zoonotic arctic fox rabies among countries bordering the Arctic, in particular Greenland. Previous molecular epidemiological studies have suggested the occurrence of one particular arctic rabies virus (RABV) lineage (arctic-3), but have been limited by a low number of available samples preventing in-depth high resolution phylogenetic analysis of RABVs at that time. However, an improved knowledge of the evolution, at a molecular level, of the circulating RABVs and a better understanding of the historical perspective of the disease in Greenland is necessary for better direct control measures on the island. These issues have been addressed by investigating the spatio-temporal genetic diversity of arctic RABVs and their reservoir host, the arctic fox, in Greenland using both full and partial genome sequences. Using a unique set of 79 arctic RABV full genome sequences from Greenland, Canada, USA (Alaska) and Russia obtained between 1977 and 2014, a description of the historic context in relation to the genetic diversity of currently circulating RABV in Greenland and neighboring Canadian Northern territories has been provided. The phylogenetic analysis confirmed delineation into four major arctic RABV lineages (arctic 1–4) with viruses from Greenland exclusively grouping into the circumpolar arctic-3 lineage. High resolution analysis enabled distinction of seven geographically distinct subclades (3.I – 3.VII) with two subclades containing viruses from both Greenland and Canada. By combining analysis of full length RABV genome sequences and host derived sequences encoding mitochondrial proteins obtained simultaneously from brain tissues of 49 arctic foxes, the interaction of viruses and their hosts was explored in detail. Such an approach can serve as a blueprint for analysis of infectious disease dynamics and virus-host interdependencies. The results showed a fine-scale spatial population structure in Greenland arctic foxes based on mitochondrial sequences, but provided no evidence for independent isolated evolutionary development of RABV in different arctic fox lineages. These data are invaluable to support future initiatives for arctic fox rabies control and elimination in Greenland.

Next to dog-mediated rabies, wildlife rabies continues to pose a public health problem, particularly in the northern hemisphere. Control of this zoonosis at the animal source has been proven the most efficient route to reduction of human rabies burden. Successful elimination of red fox-mediated rabies in Western Europe and parts of North America has demonstrated the viability of wildlife rabies control strategies. In some regions, the epidemiology of wildlife rabies is well understood; this is not the case for arctic rabies, particularly in Greenland. Previous molecular epidemiological studies demonstrated the occurrence of one particular arctic rabies virus (RABV) lineage (arctic-3) but were limited by low sample numbers and limited sequence length so as to preclude generation of high resolution phylogenetic analysis. Here, a unique set comprised of 79 complete genome sequences of RABVs from Greenland, Canada, USA (Alaska) and Russia collected over the past four decades was analysed. The use of next generation sequencing (NGS) allowed simultaneous determination of host derived sequences encoding mitochondrial proteins from the same brain tissue of 49 arctic foxes. These sequence data combined with geographical and temporal information permit the study of the genetic diversity and evolution of circulating RABVs in Greenland against the background of reservoir host genetics. The results reveal the existence of a single arctic RABV lineage (arctic-3) in Greenland, which has evolved into multiple distinct variants. These analyses provide an improved knowledge of the evolution of the circulating viruses at the molecular level and a better understanding of the historical perspective of the disease in Greenland compared to other parts of the Arctic. This knowledge will support policy on rabies control in mammalian wildlife reservoirs.

A Systematic Review of Human Bat Rabies Virus Variant Cases: Evaluating Unprotected Physical Contact with Claws and Teeth in Support of Accurate Risk Assessments

In the United States and Canada, the most recent documented cases of rabies have been attributed to bat rabies viruses (RABV). We undertook this systematic review in an effort to summarize and enhance understanding of the risk of infection for individuals who have been potentially exposed to a suspect or confirmed rabid bat. United States rabies surveillance summaries documented a total of 41 human bat-rabies virus variant verified non-transplant cases between 1990 and 2015. All cases were fatal. Seven (17.1%) of 41 cases reported a bite from a bat. Ten (24.3%) cases had unprotected physical contact (UPC); these included seven cases that had a bat land or crawl on them (contact with claws) and one case that touched a bat's teeth. Seven (17.1%) cases had probable UPC. Insectivorous bat teeth are extremely sharp and highly efficient for predation upon arthropod prey. Bats also have sharp claws on the end of their thumbs and feet. One of the most common bat RABV variants has an ability to replicate in non-neural cells. Questioning individuals about unprotected contact with bat teeth and claws (including a bat landing or crawling on a person) may help identify additional exposures.

Surveillance of Canine Rabies in the Central African Republic: Impact on Human Health and Molecular Epidemiology

BACKGROUND

Although rabies represents an important public health threat, it is still a neglected disease in Asia and Africa where it causes tens of thousands of deaths annually despite available human and animal vaccines. In the Central African Republic (CAR), an endemic country for rabies, this disease remains poorly investigated.

METHODS

To evaluate the extent of the threat that rabies poses in the CAR, we analyzed data for 2012 from the National Reference Laboratory for Rabies, where laboratory confirmation was performed by immunofluorescence and PCR for both animal and human suspected cases, and data from the only anti-rabies dispensary of the country and only place where post-exposure prophylaxis (PEP) is available. Both are located in Bangui, the capital of the CAR. For positive samples, a portion of the N gene was amplified and sequenced to determine the molecular epidemiology of circulating strains.

RESULTS

In 2012, 966 exposed persons visited the anti-rabies dispensary and 632 received a post-exposure rabies vaccination. More than 90% of the exposed persons were from Bangui and its suburbs and almost 60% of them were under 15-years of age. No rabies-related human death was confirmed. Of the 82 samples from suspected rabid dogs tested, 69 were confirmed positive. Most of the rabid dogs were owned although unvaccinated. There was a strong spatiotemporal correlation within Bangui and within the country between reported human exposures and detection of rabid dogs (P<0.001). Phylogenetic analysis indicated that three variants belonging to Africa I and II lineages actively circulated in 2012.

CONCLUSIONS

These data indicate that canine rabies was endemic in the CAR in 2012 and had a detrimental impact on human health as shown by the hundreds of exposed persons who received PEP. Implementation of effective public health interventions including mass dog vaccination and improvement of the surveillance and the access to PEP are urgently needed in this country.

Rabies in the Baltic States: Decoding a Process of Control and Elimination

Rabies is a fatal zoonosis that still causes nearly 70, 000 human deaths every year. In Europe, the oral rabies vaccination (ORV) of red foxes (Vulpes vulpes) was developed in the late 1970s and has demonstrated its effectiveness in the eradication of the disease in Western and some Central European countries. Following the accession of the three Baltic countries--Estonia, Latvia and Lithuania--to the European Union in 2004, subsequent financial support has allowed the implementation of regular ORV campaigns since 2005-2006. This paper reviews ten years of surveillance efforts and ORV campaigns in these countries resulting in the near eradication of the disease. The various factors that may have influenced the results of vaccination monitoring were assessed using generalized linear models (GLMs) on bait uptake and on herd immunity. As shown in previous studies, juveniles had lower bait uptake level than adults. For the first time, raccoon dogs (Nyctereutes procyonoides) were shown to have significantly lower bait uptake proportion compared with red foxes. This result suggests potentially altered ORV effectiveness in this invasive species compared to the red foxes. An extensive phylogenetic analysis demonstrated that the North-East European (NEE) rabies phylogroup is endemic in all three Baltic countries. Although successive oral vaccination campaigns have substantially reduced the number of detected rabies cases, sporadic detection of the C lineage (European part of Russian phylogroup) underlines the risk of reintroduction via westward spread from bordering countries. Vaccine induced cases were also reported for the first time in non-target species (Martes martes and Meles meles).

Bioecological Drivers of Rabies Virus Circulation in a Neotropical Bat Community

INTRODUCTION

In addition to the commonly accepted importance of the vampire bat in the maintenance and transmission of the rabies virus (RABV) in South America, RABV infection of other species is widely evidenced, challenging their role in the viral cycle.

METHODOLOGY / PRINCIPLES FINDINGS

To identify the bioecological drivers of RABV circulation in neotropical bat communities, we conducted a molecular and serological survey on almost 1,000 bats from 30 species, and a 4-year longitudinal survey in two colonies of vampire bats in French Guiana. RABV was molecularly detected in a common vampire and in a frugivorous bat. The sequences corresponded to haematophagous bat-related strains and were close to viruses circulating in the Brazilian Amazon region. Species' seroprevalence ranged from 0 to 20%, and the risk of seropositivity was higher in bats with a haematophagous diet, living in monospecific colonies and in dense forests. The longitudinal survey showed substantial temporal fluctuations, with individual waves of seroconversions and waning immunity. The high prevalences observed in bat communities, in most habitats and in species that do not share the same microhabitats and bioecological patterns, the temporal variations, and a rather short period of detectable antibodies as observed in recaptured vampires suggest (i) frequent exposure of animals, (ii) an ability of the infected host to control and eliminate the virus, (iii) more relaxed modes of exposure between bats than the commonly assumed infection via direct contact with saliva of infected animals, all of which should be further investigated.

CONCLUSIONS / SIGNIFICANCE

We hypothesize that RABV circulation in French Guiana is mainly maintained in the pristine forest habitats that may provide sufficient food resources to allow vampire bats, the main prevalent species, to survive and RABV to be propagated. However, on the forest edge and in disturbed areas, human activities may induce more insidious effects such as defaunation. One of the ecological consequences is the disappearance of resources for tertiary or secondary consumers. Populations of vampires may then shift to alternative resources such as cattle, domestic animals and humans. Therefore, a good forest status, allowing both a dilution effect in highly rich bat communities and the maintenance of large populations of medium-sized and large mammals used as prey by vampires, should prevent their migration to anthropized areas.

Rabies Cases in the West of China Have Two Distinct Origins

In China, rabies remains an ongoing threat to public health. Although control efforts have been effective in reducing the number of annual cases, the virus continues to spread into new areas. Tibet, Qinghai, Gansu and Ningxia in western China have, until recently, reported only a handful of events. However, since 2011, there have been increasing numbers of cases recorded in these areas. In this study, we report the collection and analysis of samples collected from these regions. We find that cases originate from two different sources. Strains collected from Gansu and Ningxia are closely related to the primary lineage associated with the current epizootic, whereas those from Tibet and Qinghai are related to the Arctic-like-2 lineage that is most commonly associated with wildlife cases in China. Thus, it appears that while the epizootic is beginning to encroach into Gansu and Ningxia, Tibet and Qinghai a significant number of rabies cases originate from wildlife.

Overall, the number of annual cases of human rabies reported in China has been decreasing since 2007. However, some Western provinces, where few human cases have been reported in recent years, are beginning to see increasing incidence of rabies. Specifically, Ningxia, Qinghai and Gansu began to report human cases respectively from 2011, 2012 and 2013, while Tibet had its first laboratory confirmed dog rabies case in 2012. Consequently, as part of the national rabies surveillance program, we collected specimens from biting dogs or human saliva from suspected rabies cases in these areas, and after sequencing positive samples, performed a phylogenetic analysis based on the nucleoprotein gene complete sequences. Our results indicate that while Ningxia and Gansu rabies strains are very close to the lineage associated with most cases in mainland China, Tibet and Qinghai strains belong to the global Arctic-like-2 clade, which is typically associated with wild life in China and neighboring countries. Thus, it appears that, rabies reemergence in the west of China has two distinct origins, and Tibet and Qinghai rabies cases were isolated events rather than an indication of the ongoing epizootic in China.

Molecular epidemiology of reemergent rabies in Yunnan Province, southwestern China

This province is a focal point for spread of rabies between Southeast Asia and China.

Yunnan Province in China borders 3 countries (Vietnam, Laos, and Myanmar) in Southeast Asia. In the 1980s, a large-scale rabies epidemic occurred in this province, which subsided by the late 1990s. However, 3 human cases of rabies in 2000 indicated reemergence of the disease in 1 county. In 2012, rabies was detected in 77 counties; 663 persons died of rabies during this new epidemic. Fifty two rabies virus strains obtained during 2008–2012 were identified and analyzed phylogenetically by sequencing the nucleoprotein gene. Of the 4 clades identified, clades YN-A and YN-C were closely related to strains from neighboring provinces, and clade YN-B was closely related to strains from Southeast Asia, but formed a distinct branch. Rabies virus diversity might be attributed to dog movements among counties, provinces, and neighboring countries. These findings suggest that Yunnan Province is a focal point for spread of rabies between Southeast Asia and China.

Genetic characterisation of the rabies virus vaccine strains used for oral immunization of foxes in Poland to estimate the effectiveness of vaccination

The main reservoir of rabies virus in Poland has been the red fox. To control rabies in wildlife, oral immunization of foxes was introduced in 1993. The vaccine is effective when it confers immunity against the virus circulating in the environment. To assess the above issue, a study of the molecular characteristics of 570-bp fragments of the N and G genes of vaccine strains SAD B19 and SAD Bern against street virus strains was performed. The results confirmed the similarity of the vaccine strains and rabies virus strains circulating in the environment and also demonstrate the genetic stability of vaccine strains that have been distributed in Poland for 20 years.

The phylogeography of rabies in Grenada, West Indies, and implications for control

In Grenada, West Indies, rabies is endemic, and is thought to be maintained in a wildlife host, the small Indian mongoose (Herpestes auropunctatus) with occasional spillover into other hosts. Therefore, the present study was undertaken to improve understanding of rabies epidemiology in Grenada and to inform rabies control policy. Mongooses were trapped island-wide between April 2011 and March 2013 and examined for the presence of Rabies virus (RABV) antigen using the direct fluorescent antibody test (dFAT) and PCR, and for serum neutralizing antibodies (SNA) using the fluorescent antibody virus neutralization test (FAVN). An additional cohort of brain samples from clinical rabies suspects submitted between April 2011 and March 2014 were also investigated for the presence of virus. Two of the 171 (1.7%) live-trapped mongooses were RABV positive by FAT and PCR, and 20 (11.7%) had SNAs. Rabies was diagnosed in 31 of the submitted animals with suspicious clinical signs: 16 mongooses, 12 dogs, 2 cats and 1 goat. Our investigation has revealed that rabies infection spread from the northeast to the southwest of Grenada within the study period. Phylogenetic analysis revealed that the viruses from Grenada formed a monophyletic clade within the cosmopolitan lineage with a common ancestor predicted to have occurred recently (6-23 years ago), and are distinct from those found in Cuba and Puerto Rico, where mongoose rabies is also endemic. These data suggest that it is likely that this specific strain of RABV was imported from European regions rather than the Americas. These data contribute essential information for any potential rabies control program in Grenada and demonstrate the importance of a sound evidence base for planning interventions.

Bat rabies in France: a 24-year retrospective epidemiological study

Since bat rabies surveillance was first implemented in France in 1989, 48 autochthonous rabies cases without human contamination have been reported using routine diagnosis methods. In this retrospective study, data on bats submitted for rabies testing were analysed in order to better understand the epidemiology of EBLV-1 in bats in France and to investigate some epidemiological trends. Of the 3176 bats submitted for rabies diagnosis from 1989 to 2013, 1.96% (48/2447 analysed) were diagnosed positive. Among the twelve recognised virus species within the Lyssavirus genus, two species were isolated in France. 47 positive bats were morphologically identified as Eptesicus serotinus and were shown to be infected by both the EBLV-1a and the EBLV-1b lineages. Isolation of BBLV in Myotis nattereri was reported once in the north-east of France in 2012. The phylogenetic characterisation of all 47 French EBLV-1 isolates sampled between 1989 and 2013 and the French BBLV sample against 21 referenced partial nucleoprotein sequences confirmed the low genetic diversity of EBLV-1 despite its extensive geographical range. Statistical analysis performed on the serotine bat data collected from 1989 to 2013 showed seasonal variation of rabies occurrence with a significantly higher proportion of positive samples detected during the autumn compared to the spring and the summer period (34% of positive bats detected in autumn, 15% in summer, 13% in spring and 12% in winter). In this study, we have provided the details of the geographical distribution of EBLV-1a in the south-west of France and the north-south division of EBLV-1b with its subdivisions into three phylogenetic groups: group B1 in the north-west, group B2 in the centre and group B3 in the north-east of France.

Molecular characterization of cryptically circulating rabies virus from ferret badgers, Taiwan

After the last reported cases of rabies in a human in 1959 and a nonhuman animal in 1961, Taiwan was considered free from rabies. However, during 2012-2013, an outbreak occurred among ferret badgers in Taiwan. To examine the origin of this virus strain, we sequenced 3 complete genomes and acquired multiple rabies virus (RABV) nucleoprotein and glycoprotein sequences. Phylogeographic analyses demonstrated that the RABV affecting the Taiwan ferret badgers (RABV-TWFB) is a distinct lineage within the group of lineages from Asia and that it has been differentiated from its closest lineages, China I (including isolates from Chinese ferret badgers) and the Philippines, 158-210 years ago. The most recent common ancestor of RABV-TWFB originated 91-113 years ago. Our findings indicate that RABV could be cryptically circulating in the environment. An understanding of the underlying mechanism might shed light on the complex interaction between RABV and its host.

Epidemiology and molecular diversity of rabies viruses in Bulgaria

A health emergency situation occurred in Bulgaria in 2007 when positive rabies cases were notified in Sofia district in the central-western part of the country, suggesting a southward spread of the disease for the first time in the last 10 years. Phylogenetic analysis on 49 isolates sampled between 2009 and 2011 showed, for the first time, evidence of the existence of NEE and D clustered lineages in Bulgaria. Their geographical distribution clearly reveals the permeability of natural barriers, as already suggested by the disease spread that occurred across the Balkan mountain range in 2007. The monitoring and passive surveillance programmes conducted since the first 2009 oral vaccination campaign, the spatio-temporal evolution of the disease in the country since 2007, and the need for further investigation of the role of jackals in virus dispersion are discussed.

Molecular characterization of three ferret badger (Melogale moschata) rabies virus isolates from Jiangxi province, China

Ferret badger (FB) rabies viruses JX09-17(fb), JX09-18 and JX10-37 were isolated from three different regions in Jiangxi province, China, in 2009 and 2010. The complete nucleotide sequence identity between these three isolates was 87-93 %. Compared with the other Chinese rabies virus isolates and vaccine strains, 101 substitutions (53 in JX10-37, 23 in JX09-17(fb) and 25 in JX09-18) in the five structural proteins were observed, and 47 of these substitutions (27 in JX10-37, 14 in JX09-17(fb) and 6 in JX09-18) were unique among lyssaviruses. Amino acid substitutions of S231 and Q333 were noted respectively in the G protein antigenic site I of JX10-37 and site III in JX09-17(fb). Phylogenetic analysis showed that JX09-17(fb) is rooted within the China I lineage, JX09-18 is in China II, and JX10-37 is independent. Evolutionary analysis and comparative sequence data indicate that isolate JX10-37 is a variant virus that diverged from canine rabies viruses around 1933 (range 1886-1963).

Molecular epidemiology of rabies virus in Poland

The paper describes a phylogenetic study of 58 Polish isolates of rabies virus collected between 1992 and 2010. Sequences of the nucleoprotein (N) and glycoprotein (G) genes approximately 600 bp long were compared with reference sequences (GenBank) of European rabies viruses from neighbouring countries. The study confirmed a very high level of homology (94.4–100 %) of the Polish rabies virus strains irrespective of the date of isolation. Two variants of rabies virus: NEE (Northeastern Europe variant) and CE (Central Europe variant), depending on the geographical place of isolation, were circulating in Poland from 1992 to 2010. The Polish rabies virus isolates showed high similarity to European RABV strains, especially those collected in Ukraine and Romania. They were clearly different from vaccine strains SAD B19 and SAD Bern, which have been used for oral vaccination of foxes against rabies in Poland since 1993.

Ecology and geography of transmission of two bat-borne rabies lineages in Chile

Rabies was known to humans as a disease thousands of years ago. In America, insectivorous bats are natural reservoirs of rabies virus. The bat species Tadarida brasiliensis and Lasiurus cinereus, with their respective, host-specific rabies virus variants AgV4 and AgV6, are the principal rabies reservoirs in Chile. However, little is known about the roles of bat species in the ecology and geographic distribution of the virus. This contribution aims to address a series of questions regarding the ecology of rabies transmission in Chile. Analyzing records from 1985–2011 at the Instituto de Salud Pública de Chile (ISP) and using ecological niche modeling, we address these questions to help in understanding rabies-bat ecological dynamics in South America. We found ecological niche identity between both hosts and both viral variants, indicating that niches of all actors in the system are undifferentiated, although the viruses do not necessarily occupy the full geographic distributions of their hosts. Bat species and rabies viruses share similar niches, and our models had significant predictive power even across unsampled regions; results thus suggest that outbreaks may occur under consistent, stable, and predictable circumstances.

The situation of rabies in America has been changing: rabies in dogs has decreased considerably, but bats are increasingly documented as natural reservoirs of other rabies variants. A significant gap exists in understanding of bat-borne rabies in Latin America. We identified bat species known to be connected with enzootic rabies with different antigenic variants in Chile, and compiled large-scale data sets by which to test for ecological niche differences among virus lineages and bat hosts. Our results begin to characterize important ecological factors affecting rabies distribution; modeling rabies in Chile allows comparisons across different latitudes and diverse landscapes. We found that rabies virus strains are found in similar environments, regardless of the bat host involved. This research improves understanding of bat-borne rabies dynamics, and important step towards preventing and controlling this and other emergent diseases linked to bats.

The phylogeography and spatiotemporal spread of south-central skunk rabies virus

The south-central skunk rabies virus (SCSK) is the most broadly distributed terrestrial viral lineage in North America. Skunk rabies has not been efficiently targeted by oral vaccination campaigns and represents a natural system of pathogen invasion, yielding insights to rabies emergence. In the present study we reconstructed spatiotemporal spread of SCSK in the whole territory of its circulation using a combination of Bayesian methods. The analysis based on 241 glycoprotein gene sequences demonstrated that SCSK is much more divergent phylogenetically than was appreciated previously. According to our analyses the SCSK originated in the territory of Texas ~170 years ago, and spread geographically during the following decades. The wavefront velocity in the northward direction was significantly greater than in the eastward and westward directions. Rivers (except the Mississippi River and Rio Grande River) did not constitute significant barriers for epizootic spread, in contrast to deserts and mountains. The mean dispersal rate of skunk rabies was lower than that of the raccoon and fox rabies. Viral lineages circulate in their areas with limited evidence of geographic spread during decades. However, spatiotemporal reconstruction shows that after a long period of stability the dispersal rate and wavefront velocity of SCSK are increasing. Our results indicate that there is a need to develop control measures for SCSK, and suggest how such measure can be implemented most efficiently. Our approach can be extrapolated to other rabies reservoirs and used as a tool for investigation of epizootic patterns and planning interventions towards disease elimination.

Recent emergence and spread of an Arctic-related phylogenetic lineage of rabies virus in Nepal

Rabies is a zoonotic disease that is endemic in many parts of the developing world, especially in Africa and Asia. However its epidemiology remains largely unappreciated in much of these regions, such as in Nepal, where limited information is available about the spatiotemporal dynamics of the main etiological agent, the rabies virus (RABV). In this study, we describe for the first time the phylogenetic diversity and evolution of RABV circulating in Nepal, as well as their geographical relationships within the broader region. A total of 24 new isolates obtained from Nepal and collected from 2003 to 2011 were full-length sequenced for both the nucleoprotein and the glycoprotein genes, and analysed using neighbour-joining and maximum-likelihood phylogenetic methods with representative viruses from all over the world, including new related RABV strains from neighbouring or more distant countries (Afghanistan, Greenland, Iran, Russia and USA). Despite Nepal's limited land surface and its particular geographical position within the Indian subcontinent, our study revealed the presence of a surprising wide genetic diversity of RABV, with the co-existence of three different phylogenetic groups: an Indian subcontinent clade and two different Arctic-like sub-clades within the Arctic-related clade. This observation suggests at least two independent episodes of rabies introduction from neighbouring countries. In addition, specific phylogenetic and temporal evolution analysis of viruses within the Arctic-related clade has identified a new recently emerged RABV lineage we named as the Arctic-like 3 (AL-3) sub-clade that is already widely spread in Nepal.

Rabies is endemic in most Asian countries and represents a serious public health issue, with an estimated 31,000 people dying each year of this disease. The majority of human cases are transmitted by domestic dogs, which act as the principal reservoir host and vector. However, molecular epidemiology and evolutionary dynamics of the main etiological agent, the rabies virus (RABV), remains largely unappreciated in some regions such as in Nepal. Based on a subset of 24 new Nepalese isolates collected from 2003 to 2011 and representative RABV strains at a global scale, phylogenetic analysis based on the complete nucleoprotein and glycoprotein genes sequences revealed the presence of a surprising wide genetic diversity of RABV circulating in this country. The presence of three different co-existing phylogenetic groups was identified: an Indian subcontinent clade and two different Arctic-like sub-clades within the Arctic-related clade, namely Arctic-like (AL)-1, lineage a (AL-1a), and AL-3. Among these clusters, the AL-3 sub-clade appears as the major Nepalese phylogroup which emerged relatively recently in this country, within the last 30 years. These data has raised some concerns about the exchange of RABV between different countries, and provided key elements for implementation of effective control measures of rabies in Nepal.

Molecular Epidemiology of Rabies in Vietnam

The present study was done to determine the molecular epidemiology of rabies virus (RV) in Vietnam. The nucleoprotein (N) and glycoprotein (G) genes of RVs were amplified from the brains of ten rabid dogs of Ho Chi Minh City, Vietnam. The nucleotide sequences of these genes were compared with those of other Asian strains to find the possible relationship among them. Phylogenetic analysis revealed that the Asian N gene segregated into three main branches, namely South-East Asia 1 (SEA 1), South-East Asia 2 (SEA 2) and Indian subcontinent (ISC) genotypes. The SEA 1 genotype comprised RVs from Malaysia, Vietnam and Thailand. The SEA 2 genotype contained strains from the Philippines, and the ISC genotype comprised strains from Sri Lanka and India. Phylogenetically G genes of RVs from Vietnam and Thailand were clustered together. Our study suggests that Vietnamese and Thai RVs are closely related and might have originated from a common ancestor.

Genetic Typing of Feline Rabies Virus Isolated in Greater Bangkok, Thailand

To study the molecular epidemiology of rabies virus that is prevalent among cats in greater Bangkok, Thailand, a total of 17 rabies virus isolates from cats were characterized and compared with 120 rabies virus isolates from dogs. Analyses were performed on the genetic polymorphism in the rabies virus nucleoprotein (N) gene. Rabies virus N gene of isolates was amplified by reverse transcriptionpolymerase chain reaction. The diversity of N gene was revealed by the restriction fragment length polymorphism (RFLP) method. The rabies virus isolates from cats could be classified into 5 types, designated as Dd I-Hf I, Dd II-Hf II, Dd III-Hf I, Dd IV-Hf I, and Dd IV-Hf III. Type Dd I-Hf I was encountered more frequently than the others. It was apparent that no less than five rabies virus types presented in the areas of Bangkok. Moreover, all five RFLP patterns were typical of those which had been observed in dogs. Our findings suggest that there had been viral transmission between the dogs and the cats.

Evolutionary history and phylogeography of rabies viruses associated with outbreaks in Trinidad

Bat rabies is an emerging disease of public health significance in the Americas. The Caribbean island of Trinidad experiences periodic outbreaks within the livestock population. We performed molecular characterisation of Trinidad rabies virus (RABV) and used a Bayesian phylogeographic approach to investigate the extent to which outbreaks are a result of in situ evolution versus importation of virus from the nearby South American mainland. Trinidadian RABV sequences were confirmed as bat variant and clustered with Desmodus rotundus (vampire bat) related sequences. They fell into two largely temporally defined lineages designated Trinidad I and II. The Trinidad I lineage which included sequences from 1997-2000 (all but two of which were from the northeast of the island) was most closely related to RABV from Ecuador (2005, 2007), French Guiana (1990) and Venezuela (1993, 1994). Trinidad II comprised sequences from the southwest of the island, which clustered into two groups: Trinidad IIa, which included one sequence each from 2000 and 2007, and Trinidad IIb including all 2010 sequences. The Trinidad II sequences were most closely related to sequences from Brazil (1999, 2004) and Uruguay (2007, 2008). Phylogeographic analyses support three separate RABV introductions from the mainland from which each of the three Trinidadian lineages arose. The estimated dates for the introductions and subsequent lineage expansions suggest periods of in situ evolution within Trinidad following each introduction. These data also indicate co-circulation of Trinidad lineage I and IIa during 2000. In light of these findings and the likely vampire bat origin of Trinidadian RABV, further studies should be conducted to investigate the relationship between RABV spatiotemporal dynamics and vampire bat population ecology, in particular any movement between the mainland and Trinidad.

Genetic strain modification of a live rabies virus vaccine widely used in Europe for wildlife oral vaccination

In Europe, the main reservoir and vector of rabies has been the red fox (Vulpes vulpes). Oral immunization of foxes with live vaccines, using attenuated rabies strains (SAD B19, SAD Bern), apathogenic mutants of an attenuated strain (SAG2) and the vaccinia-rabies glycoprotein recombinant virus vaccine (V-RG), has been shown to be the most effective method for the control and elimination of rabies. Among all vaccines currently used for wildlife oral vaccination, one vaccine (marketed as SAD Bern strain) has been widely used in Europe since 1992 with the distribution of 17million of baits in 2011. Because of the potential environmental safety risk of a live virus which could revert to virulence, the full genome sequencing of this vaccine was undertaken and the sequence was characterized and compared with those of referenced rabies viruses. The vaccine showed higher similarity to the strains belonging to the SAD B19 vaccine virus strains than to the SAD Bern vaccines. This study is the first one reporting on virus strain identity changes in this attenuated vaccine.