Phenotypic assessment of antiviral activity for spiro-annulated oxepanes and azepenes.

Evolutionary potential of viruses can result in outbreaks of well-known viruses and emergence of novel ones. Pharmacological methods of intervening the reproduction of various less popular, but not less important viruses are not available, as well as the spectrum of antiviral activity for most known compounds. In the framework of chemical biology paradigm, characterization of antiviral activity spectrum of new compounds allows to extend the antiviral chemical space and provides new important structure-activity relationships for data-driven drug discovery. Here we present a primary assessment of antiviral activity of spiro-annulated derivatives of seven-membered heterocycles, oxepane and azepane, in phenotypic assays against viruses with different genomes, virion structures, and genome realization schemes: orthoflavivirus (tick-borne encephalitis virus, TBEV), enteroviruses (poliovirus, enterovirus A71, echovirus 30), adenovirus (human adenovirus C5), hantavirus (Puumala virus). Hit compounds inhibited reproduction of adenovirus C5, the only DNA virus in the studied set, in the yield reduction assay, and did not inhibit reproduction of RNA viruses.

m6A Methylation in Regulation of Antiviral Innate Immunity.

The epitranscriptomic modification m6A is a prevalent RNA modification that plays a crucial role in the regulation of various aspects of RNA metabolism. It has been found to be involved in a wide range of physiological processes and disease states. Of particular interest is the role of m6A machinery and modifications in viral infections, serving as an evolutionary marker for distinguishing between self and non-self entities. In this review article, we present a comprehensive overview of the epitranscriptomic modification m6A and its implications for the interplay between viruses and their host, focusing on immune responses and viral replication. We outline future research directions that highlight the role of m6A in viral nucleic acid recognition, initiation of antiviral immune responses, and modulation of antiviral signaling pathways. Additionally, we discuss the potential of m6A as a prognostic biomarker and a target for therapeutic interventions in viral infections.

The menace of severe adverse events and deaths associated with viral gene therapy and its potential solution.

Over the past decade, in vivo gene replacement therapy has significantly advanced, resulting in market approval of numerous therapeutics predominantly relying on adeno-associated viral vectors (AAV). While viral vectors have undeniably addressed several critical healthcare challenges, their clinical application has unveiled a range of limitations and safety concerns. This review highlights the emerging challenges in the field of gene therapy. At first, we discuss both the role of biological barriers in viral gene therapy with a focus on AAVs, and review current landscape of in vivo human gene therapy. We delineate advantages and disadvantages of AAVs as gene delivery vehicles, mostly from the safety perspective (hepatotoxicity, cardiotoxicity, neurotoxicity, inflammatory responses etc.), and outline the mechanisms of adverse events in response to AAV. Contribution of every aspect of AAV vectors (genomic structure, capsid proteins) and host responses to injected AAV is considered and substantiated by basic, translational and clinical studies. The updated evaluation of recent AAV clinical trials and current medical experience clearly shows the risks of AAVs that sometimes overshadow the hopes for curing a hereditary disease. At last, a set of established and new molecular and nanotechnology tools and approaches are provided as potential solutions for mitigating or eliminating side effects. The increasing number of severe adverse reactions and, sadly deaths, demands decisive actions to resolve the issue of immune responses and extremely high doses of viral vectors used for gene therapy. In response to these challenges, various strategies are under development, including approaches aimed at augmenting characteristics of viral vectors and others focused on creating secure and efficacious non-viral vectors. This comprehensive review offers an overarching perspective on the present state of gene therapy utilizing both viral and non-viral vectors.

Non-Polio Enteroviruses Isolated by Acute Flaccid Paralysis Surveillance Laboratories in the Russian Federation in 1998-2021: Distinct Epidemiological Features of Types.

More than 100 types of non-polio enteroviruses (NPEVs) are ubiquitous in the human population and cause a variety of symptoms ranging from very mild to meningitis and acute flaccid paralysis (AFP). Much of the information regarding diverse pathogenic properties of NPEVs comes from the surveillance of poliovirus, which also yields NPEV. The analysis of 265 NPEV isolations from 10,433 AFP cases over 24 years of surveillance and more than 2500 NPEV findings in patients without severe neurological lesions suggests that types EV-A71, E13, and E25 were significantly associated with AFP. EV-A71 was also significantly more common among AFP patients who had fever at the onset and residual paralysis compared to all AFP cases. In addition, a significant disparity was noticed between types that were common in humans (CV-A2, CVA9, EV-A71, E9, and E30) or in sewage (CVA7, E3, E7, E11, E12, and E19). Therefore, there is significant evidence of non-polio viruses being implicated in severe neurological lesions, but further multicenter studies using uniform methodology are needed for a definitive conclusion.

Transient and tunable CRISPRa regulation of APOBEC/AID genes for targeting hepatitis B virus.

APOBEC/AID cytidine deaminases play an important role in innate immunity and antiviral defenses and were shown to suppress hepatitis B virus (HBV) replication by deaminating and destroying the major form of HBV genome, covalently closed circular DNA (cccDNA), without toxicity to the infected cells. However, developing anti-HBV therapeutics based on APOBEC/AID is complicated by the lack of tools for activating and controlling their expression. Here, we developed a CRISPR-activation-based approach (CRISPRa) to induce APOBEC/AID transient overexpression (>4-800,000-fold increase in mRNA levels). Using this new strategy, we were able to control APOBEC/AID expression and monitor their effects on HBV replication, mutation, and cellular toxicity. CRISPRa prominently reduced HBV replication (∼90%-99% decline of viral intermediates), deaminated and destroyed cccDNA, but induced mutagenesis in cancer-related genes. By coupling CRISPRa with attenuated sgRNA technology, we demonstrate that APOBEC/AID activation can be precisely controlled, eliminating off-site mutagenesis in virus-containing cells while preserving prominent antiviral activity. This study untangles the differences in the effects of physiologically expressed APOBEC/AID on HBV replication and cellular genome, provides insights into the molecular mechanisms of HBV cccDNA mutagenesis, repair, and degradation, and, finally, presents a strategy for a tunable control of APOBEC/AID expression and for suppressing HBV replication without toxicity.

Depleting hepatitis B virus relaxed circular DNA is necessary for resolution of infection by CRISPR-Cas9.

CRISPR-Cas9 systems can directly target the hepatitis B virus (HBV) major genomic form, covalently closed circular DNA (cccDNA), for decay and demonstrate remarkable anti-HBV activity. Here, we demonstrate that CRISPR-Cas9-mediated inactivation of HBV cccDNA, frequently regarded as the "holy grail" of viral persistence, is not sufficient for curing infection. Instead, HBV replication rapidly rebounds because of de novo formation of HBV cccDNA from its precursor, HBV relaxed circular DNA (rcDNA). However, depleting HBV rcDNA before CRISPR-Cas9 ribonucleoprotein (RNP) delivery prevents viral rebound and promotes resolution of HBV infection. These findings provide the groundwork for developing approaches for a virological cure of HBV infection by a single dose of short-lived CRISPR-Cas9 RNPs. Blocking cccDNA replenishment and re-establishment from rcDNA conversion is critical for completely clearing the virus from infected cells by site-specific nucleases. The latter can be achieved by widely used reverse transcriptase inhibitors.

Genomic analysis of Leishmania turanica strains from different regions of Central Asia.

The evolution in Leishmania is governed by the opposite forces of clonality and sexual reproduction, with vicariance being an important factor. As such, Leishmania spp. populations may be monospecific or mixed. Leishmania turanica in Central Asia is a good model to compare these two types. In most areas, populations of L. turanica are mixed with L. gerbilli and L. major. Notably, co-infection with L. turanica in great gerbils helps L. major to withstand a break in the transmission cycle. Conversely, the populations of L. turanica in Mongolia are monospecific and geographically isolated. In this work, we compare genomes of several well-characterized strains of L. turanica originated from monospecific and mixed populations in Central Asia in order to shed light on genetic factors, which may drive evolution of these parasites in different settings. Our results illustrate that evolutionary differences between mixed and monospecific populations of L. turanica are not dramatic. On the level of large-scale genomic rearrangements, we confirmed that different genomic loci and different types of rearrangements may differentiate strains originated from mixed and monospecific populations, with genome translocations being the most prominent example. Our data suggests that L. turanica has a significantly higher level of chromosomal copy number variation between the strains compared to its sister species L. major with only one supernumerary chromosome. This suggests that L. turanica (in contrast to L. major) is in the active phase of evolutionary adaptation.

Hydroxychloroquine Enhances Cytotoxic Properties of Extracellular Vesicles and Extracellular Vesicle-Mimetic Nanovesicles Loaded with Chemotherapeutics.

Because of their high biocompatibility, biological barrier negotiation, and functionalization properties, biological nanoparticles have been actively investigated for many medical applications. Biological nanoparticles, including natural extracellular vesicles (EVs) and synthetic extracellular vesicle-mimetic nanovesicles (EMNVs), represent novel drug delivery vehicles that can accommodate different payloads. In this study, we investigated the physical, biological, and delivery properties of EVs and EMNVs and analyzed their ability to deliver the chemotherapeutic drug doxorubicin. EMNVs and EVs exhibit similar properties, but EMNVs are more effectively internalized, while EVs show higher intracellular doxorubicin release activity. In addition, these nanotherapeutics were investigated in combination with the FDA-approved drug hydroxychloroquine (HCQ). We demonstrate that HCQ-induced lysosome destabilization and could significantly increase nanoparticle internalization, doxorubicin release, and cytotoxicity. Altogether, these data demonstrate that, from the delivery standpoint in vitro, the internalization of EMNVs and EVs and their payload release were slightly different and both nanotherapeutics had comparable cytotoxic performance. However, the synthesis of EMNVs was significantly faster and cost-effective. In addition, we highlight the benefits of combining biological nanoparticles with the lysosome-destabilizing agent HCQ that increased both the internalization and the cytotoxic properties of the particles.

Patterns and Temporal Dynamics of Natural Recombination in Noroviruses.

Noroviruses infect a wide range of mammals and are the major cause of gastroenteritis in humans. Recombination at the junction of ORF1 encoding nonstructural proteins and ORF2 encoding major capsid protein VP1 is a well-known feature of noroviruses. Using all available complete norovirus sequences, we systematically analyzed patterns of natural recombination in the genus Norovirus both throughout the genome and across the genogroups. Recombination events between nonstructural (ORF1) and structural genomic regions (ORF2 and ORF3) were found in all analyzed genogroups of noroviruses, although recombination was most prominent between members of GII, the most common genogroup that infects humans. The half-life times of recombinant forms (clades without evidence of recombination) of human GI and GII noroviruses were 10.4 and 8.4-11.3 years, respectively. There was evidence of many recent recombination events, and most noroviruses that differed by more than 18% of nucleotide sequence were recombinant relative to each other. However, there were no distinct recombination events between viruses that differed by over 42% in ORF2/3, consistent with the absence of systematic recombination between different genogroups. The few inter-genogroup recombination events most likely occurred between ancient viruses before they diverged into contemporary genogroups. The recombination events within ORF1 or between ORF2/3 were generally rare. Thus, noroviruses routinely exchange full structural and nonstructural blocks of the genome, providing a modular evolution.

Technological aspects of manufacturing and analytical control of biological nanoparticles.

Extracellular vesicles (EVs) are cell-derived biological nanoparticles that gained great interest for drug delivery. EVs have numerous advantages compared to synthetic nanoparticles, such as ideal biocompatibility, safety, ability to cross biological barriers and surface modification via genetic or chemical methods. On the other hand, the translation and the study of these carriers resulted difficult, mostly because of significant issues in up-scaling, synthesis and impractical methods of quality control. However, current manufacturing advances enable EV packaging with any therapeutic cargo, including DNA, RNA (for RNA vaccines and RNA therapeutics), proteins, peptides, RNA-protein complexes (including gene-editing complexes) and small molecules drugs. To date, an array of new and upgraded technologies have been introduced, substantially improving EV production, isolation, characterization and standardization. The used-to-be "gold standards" of EV manufacturing are now outdated, and the state-of-art requires extensive revision. This review re-evaluates the pipeline for EV industrial production and provides a critical overview of the modern technologies required for their synthesis and characterization.

Revisiting epidemiology of leishmaniasis in central Asia: lessons learnt.

In this work we reviewed historical and recent data on Leishmania spp. infection combining data collected in Turkmenistan, Uzbekistan, Kazakhstan, Kyrgyzstan, Iran, China and Mongolia. We specifically focused on a complex of co-existing species (Leishmania major, Leishmania turanica and Leishmania gerbilli) sharing the same animal reservoirs and vectors. In addition, we analysed the presence of dsRNA viruses in these species and discussed future research directions to identify species-specific traits, which may determine susceptibility of different Leishmania spp. to viral infection.

Robust AAV Genotyping Based on Genetic Distances in Rep Gene That Are Maintained by Ubiquitous Recombination.

Adeno-associated viruses (AAVs) are a convenient tool for gene therapy delivery. According to the current classification, they are divided into the species AAV A and AAV B within the genus Dependoparvovirus. Historically AAVs were also subdivided on the intraspecies level into 13 serotypes, which differ in tissue tropism and targeted gene delivery capacity. Serotype, however, is not a universal taxonomic category, and their assignment is not always robust. Cross-reactivity has been shown, indicating that classification could not rely on the results of serological tests alone. Moreover, since the isolation of AAV4, all subsequent AAVs were subdivided into serotypes based primarily on genetic differences and phylogenetic reconstructions. An increased interest in the use of AAV as a gene delivery tool justifies the need to improve the existing classification. Here, we suggest genotype-based AAV classification below the species level based on the rep gene. A robust threshold was established as 10% nt differences within the 1248 nt genome fragment, with 4 distinct AAV genotypes identified. This distinct sub-species structure is maintained by ubiquitous recombination within, but not between, rep genes of the suggested genotypes.

Cases of Acute Flaccid Paralysis Associated with Coxsackievirus A2: Findings of a 20-Year Surveillance in the Russian Federation.

Surveillance for acute flaccid paralysis syndrome (AFP) in children under 15 is the backbone of the Global Polio Eradication Initiative. Laboratory examination of stool samples from AFP cases allows the detection of, along with polioviruses, a variety of non-polio enteroviruses (NPEV). The etiological significance of these viruses in the occurrence of AFP cases has been definitively established only for enteroviruses A71 and D68. Enterovirus Coxsackie A2 (CVA2) is most often associated with vesicular pharyngitis and hand, foot and mouth disease. Among 7280 AFP cases registered in Russia over 20 years (2001-2020), CVA2 was isolated only from five cases. However, these included three children aged 3 to 4 years, without overt immune deficiency, immunized with 4-5 doses of poliovirus vaccine in accordance with the National Vaccination Schedule. The disease resulted in persistent residual paralysis. Clinical and laboratory data corresponded to poliomyelitis developing during poliovirus infection. These findings are compatible with CVA2 being the cause of AFP. Molecular analysis of CVA2 from these patients and a number of AFP cases in other countries did not reveal association with a specific phylogenetic group, suggesting that virus genetics is unlikely to explain the pathogenic profile. The overall results highlight the value of AFP surveillance not just for polio control but for studies of uncommon AFP agents.

Analyses of Leishmania-LRV Co-Phylogenetic Patterns and Evolutionary Variability of Viral Proteins.

Leishmania spp. are important pathogens causing a vector-borne disease with a broad range of clinical manifestations from self-healing ulcers to the life-threatening visceral forms. Presence of Leishmania RNA virus (LRV) confers survival advantage to these parasites by suppressing anti-leishmanial immunity in the vertebrate host. The two viral species, LRV1 and LRV2 infect species of the subgenera Viannia and Leishmania, respectively. In this work we investigated co-phylogenetic patterns of leishmaniae and their viruses on a small scale (LRV2 in L. major) and demonstrated their predominant coevolution, occasionally broken by intraspecific host switches. Our analysis of the two viral genes, encoding the capsid and RNA-dependent RNA polymerase (RDRP), revealed them to be under the pressure of purifying selection, which was considerably stronger for the former gene across the whole tree. The selective pressure also differs between the LRV clades and correlates with the frequency of interspecific host switches. In addition, using experimental (capsid) and predicted (RDRP) models we demonstrated that the evolutionary variability across the structure is strikingly different in these two viral proteins.

Modular Evolution of Coronavirus Genomes.

The viral family Coronaviridae comprises four genera, termed Alpha-, Beta-, Gamma-, and Deltacoronavirus. Recombination events have been described in many coronaviruses infecting humans and other animals. However, formal analysis of the recombination patterns, both in terms of the involved genome regions and the extent of genetic divergence between partners, are scarce. Common methods of recombination detection based on phylogenetic incongruences (e.g., a phylogenetic compatibility matrix) may fail in cases where too many events diminish the phylogenetic signal. Thus, an approach comparing genetic distances in distinct genome regions (pairwise distance deviation matrix) was set up. In alpha, beta, and delta-coronaviruses, a low incidence of recombination between closely related viruses was evident in all genome regions, but it was more extensive between the spike gene and other genome regions. In contrast, avian gammacoronaviruses recombined extensively and exist as a global cloud of genes with poorly corresponding genetic distances in different parts of the genome. Spike, but not other structural proteins, was most commonly exchanged between coronaviruses. Recombination patterns differed between coronavirus genera and corresponded to the modular structure of the spike: recombination traces were more pronounced between spike domains (N-terminal and C-terminal parts of S1 and S2) than within domains. The variability of possible recombination events and their uneven distribution over the genome suggest that compatibility of genes, rather than mechanistic or ecological limitations, shapes recombination patterns in coronaviruses.

The effects of genetic drift and genomic selection on differentiation and local adaptation of the introduced populations of Aedes albopictus in southern Russia.

BACKGROUND: Asian tiger mosquito Aedes albopictus is an arbovirus vector that has spread from its native habitation areal in Southeast Asia throughout North and South Americas, Europe, and Africa. Ae. albopictus was first detected in the Southern Federal District of the Russian Federation in the subtropical town of Sochi in 2011. In subsequent years, this species has been described in the continental areas with more severe climate and lower winter temperatures. METHODS: Genomic analysis of pooled Ae. albopictus samples collected in the mosquito populations in the coastal and continental regions of the Krasnodar Krai was conducted to look for the genetic changes associated with the spread and potential cold adaptation in Ae. albopictus. RESULTS: The results of the phylogenetic analysis based on mitochondrial genomes corresponded well with the hypothesis that Ae. albopictus haplotype A1a2a1 was introduced into the region from a single source. Population analysis revealed the role of dispersal and genetic drift in the local adaptation of the Asian tiger mosquito. The absence of shared haplotypes between the samples and high fixation indices suggest that gene flow between samples was heavily restricted. Mitochondrial and genomic differentiation together with different distances between dispersal routes, natural and anthropogenic barriers and local effective population size reduction could lead to difficulties in local climatic adaptations due to reduced selection effectiveness. We have found genomic regions with selective sweep patterns which can be considered as having been affected by recent selection events. The genes located in these regions participate in neural protection, lipid conservation, and cuticle formation during diapause. These processes were shown to be important for cold adaptation in the previous transcriptomic and proteomic studies. However, the population history and relatively low coverage obtained in the present article could have negatively affect sweep detection.

Citizen science initiative points at childhood BCG vaccination as a risk factor for COVID-19.

Current results do not provide conclusive evidence on the effect of BCG vaccination on COVID-19 alone or in combination with other factors. To address this limitation, in this study we used a citizen science initiative on the COVID-19 pandemic to collect data worldwide during 2 October 2020-30 October 2020 (1,233 individuals) in a structured way for analysing factors and characteristics of affected individuals in relation to BCG vaccination. For the first time, the results of our study suggested that vaccination with BCG may increase the risk for COVID-19 at certain age, particularly in individuals vaccinated at childhood. Childhood BCG vaccination increased the likelihood of being diagnosed with COVID-19 fivefold in COVID-19 low-incidence countries and threefold in high-incidence countries. A reasonable explanation for this effect is the activation of certain innate immunity mechanisms associated with inflammatory reactions. These factors should be considered when analysing the risks associated with this global pandemic.

A hepatitis B virus causes chronic infections in equids worldwide.

Preclinical testing of novel therapeutics for chronic hepatitis B (CHB) requires suitable animal models. Equids host homologs of hepatitis C virus (HCV). Because coinfections of hepatitis B virus (HBV) and HCV occur in humans, we screened 2,917 specimens from equids from five continents for HBV. We discovered a distinct HBV species (Equid HBV, EqHBV) in 3.2% of donkeys and zebras by PCR and antibodies against EqHBV in 5.4% of donkeys and zebras. Molecular, histopathological, and biochemical analyses revealed that infection patterns of EqHBV resembled those of HBV in humans, including hepatotropism, moderate liver damage, evolutionary stasis, and potential horizontal virus transmission. Naturally infected donkeys showed chronic infections resembling CHB with high viral loads of up to 2.6 × 109 mean copies per milliliter serum for >6 mo and weak antibody responses. Antibodies against Equid HCV were codetected in 26.5% of donkeys seropositive for EqHBV, corroborating susceptibility to both hepatitis viruses. Deltavirus pseudotypes carrying EqHBV surface proteins were unable to infect human cells via the HBV receptor NTCP (Na+/taurocholate cotransporting polypeptide), suggesting alternative viral entry mechanisms. Both HBV and EqHBV deltavirus pseudotypes infected primary horse hepatocytes in vitro, supporting a broad host range for EqHBV among equids and suggesting that horses might be suitable for EqHBV and HBV infections in vivo. Evolutionary analyses suggested that EqHBV originated in Africa several thousand years ago, commensurate with the domestication of donkeys. In sum, EqHBV naturally infects diverse equids and mimics HBV infection patterns. Equids provide a unique opportunity for preclinical testing of novel therapeutics for CHB and to investigate HBV/HCV interplay upon coinfection.

At Least Seven Distinct Rotavirus Genotype Constellations in Bats with Evidence of Reassortment and Zoonotic Transmissions.

Bats host many viruses pathogenic to humans, and increasing evidence suggests that rotavirus A (RVA) also belongs to this list. Rotaviruses cause diarrheal disease in many mammals and birds, and their segmented genomes allow them to reassort and increase their genetic diversity. Eighteen out of 2,142 bat fecal samples (0.8%) collected from Europe, Central America, and Africa were PCR-positive for RVA, and 11 of those were fully characterized using viral metagenomics. Upon contrasting their genomes with publicly available data, at least 7 distinct bat RVA genotype constellations (GCs) were identified, which included evidence of reassortments and 6 novel genotypes. Some of these constellations are spread across the world, whereas others appear to be geographically restricted. Our analyses also suggest that several unusual human and equine RVA strains might be of bat RVA origin, based on their phylogenetic clustering, despite various levels of nucleotide sequence identities between them. Although SA11 is one of the most widely used reference strains for RVA research and forms the backbone of a reverse genetics system, its origin remained enigmatic. Remarkably, the majority of the genotypes of SA11-like strains were shared with Gabonese bat RVAs, suggesting a potential common origin. Overall, our findings suggest an underexplored genetic diversity of RVAs in bats, which is likely only the tip of the iceberg. Increasing contact between humans and bat wildlife will further increase the zoonosis risk, which warrants closer attention to these viruses.IMPORTANCE The increased research on bat coronaviruses after severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) allowed the very rapid identification of SARS-CoV-2. This is an excellent example of the importance of knowing viruses harbored by wildlife in general, and bats in particular, for global preparedness against emerging viral pathogens. The current effort to characterize bat rotavirus strains from 3 continents sheds light on the vast genetic diversity of rotaviruses and also hints at a bat origin for several atypical rotaviruses in humans and animals, implying that zoonoses of bat rotaviruses might occur more frequently than currently realized.

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.