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

Development and validation of an all-in-one rabies virus Bat-Clade genomic sequencing and host identification protocol

Rabies virus (RABV), remains a significant public health concern, with bat-maintained lineages accounting for all currently documented cases in Brazil. Despite the availability of pharmacological prophylaxis for humans and animals, the high genetic diversity of RABV in diverse natural bat hosts and continued circulation in multiple animals pose challenges for effective surveillance. Here, we developed and validated a novel, rapidly deployable amplicon-based sequencing approach for RABV genomic surveillance. This "all-in-one" protocol integrates whole RABV genome sequencing with host species identification through COI gene amplification and sequencing, addressing the challenges posed by RABV's high genetic diversity and complex transmission dynamics. We assessed the protocol's effectiveness by sequencing 25 near-complete RABV genomes from host species across four distinct families (Bovidae, Equidae, Felidae, and Microchiroptera) obtained from the Rabies Control and Surveillance Program from the Rio Grande do Sul State, Southern Brazil. The method achieved an average genome coverage of 91.4 % at a minimum 5x read depth, with a mean depth coverage of 816x across sequenced genomes. The results demonstrated significant Bat-Clade sublineage diversity, which was classified using the MADDOG RABV lineage system. The protocol successfully identified three bat species (Tadarida brasiliensis, Desmodus rotundus, and Myotis nigricans) among the samples, highlighting its capability for precise host identification. This study presents a powerful tool for high-resolution evaluation of RABV genomic features and host identification, enabling more targeted animal and human health interventions. This new approach has the potential to enhance RABV surveillance capabilities, contributing to more effective rabies control strategies within a One Health framework.

Integrated Community-Based Reporting and Field Diagnostics for Improved Rabies Surveillance in Rural Laikipia, Kenya

Rabies vaccination in domestic dog populations has increased globally in a bid to protect human health. Surveillance efforts, however, are inconsistent in endemic regions such as in sub-Saharan Africa, due to fragmented reporting and limited diagnostic capacity for suspected cases, limiting successful monitoring and evaluation of vaccination campaigns. Here, we conducted a pilot study aiming to strengthen rabies surveillance by combining community-based surveillance with field-based diagnostic testing in pastoral and agro-pastoral communities in central Kenya; communities which are frequently marginalised from health systems. During the 6-month pilot study, there were 14 alerts of suspected rabid dogs in the community, of which eight were tested and five diagnostically confirmed as rabid. Two positive samples processed successfully for whole genome sequencing indicated that the rabies variant circulating in central Kenya during the study period belonged to the Africa 1b subclade, which is similar to variants identified in eastern Kenya and Tanzania, suggesting regional transmission. This pilot study indicates that rabies continues to circulate in the region and that community-based surveillance, when combined with enhanced diagnostic testing, can help alleviate underreporting and guide vaccination campaigns.

Genomic characterization of a dog-mediated rabies outbreak in El Pedregal, Arequipa, Peru

BACKGROUND

Rabies, a re-emerging zoonosis with the highest known human case fatality rate, has been largely absent from Peru, except for endemic circulation in the Puno region on the Bolivian border and re-emergence in Arequipa City in 2015, where it has persisted. In 2021, an outbreak occurred in the rapidly expanding city of El Pedregal near Arequipa, followed by more cases in 2022 after nearly a year of epidemiological silence. While currently under control, questions persist regarding the origin of the El Pedregal outbreak and implications for maintaining rabies control in Peru.

METHODS

We sequenced 25 dog rabies virus (RABV) genomes from the El Pedregal outbreak (n=11) and Arequipa City (n=14) from 2021-2023 using Nanopore sequencing in Peru. Historical genomes from Puno (n=4, 2010-2012) and Arequipa (n=5, 2015-2019), were sequenced using an Illumina approach in the UK. In total, 34 RABV genomes were generated, including archived and newly obtained samples. The genomes were analyzed phylogenetically to understand the outbreak's context and origins.

RESULTS

Phylogenomic analysis identified two genetic clusters in El Pedregal: 2021 cases stemmed from a single introduction unrelated to Arequipa cases, while the 2022 sequence suggested a new introduction from Arequipa rather than persistence. In relation to canine RABV diversity in Latin America, all new sequences belonged to the new minor clade, Cosmopolitan Am5, sharing relatives from Bolivia, Argentina, and Brazil.

CONCLUSION

Genomic insights into the El Pedregal outbreak revealed multiple introductions over a 2-year window. Eco-epidemiological conditions, including migratory worker patterns, suggest human-mediated movement drove introductions. Despite outbreak containment, El Pedregal remains at risk of dog-mediated rabies due to ongoing circulation in Arequipa, Puno, and Bolivia. Human-mediated movement of dogs presents a major risk for rabies re-emergence in Peru, jeopardizing regional dog-mediated rabies control. Additional sequence data is needed for comprehensive phylogenetic analyses.

Craft: A Machine Learning Approach to Dengue Subtyping

Motivation

The dengue virus poses a major global health threat, with nearly 390 million infections annually. A recently proposed hierarchical dengue nomenclature system enhances spatial resolution by defining major and minor lineages within genotypes, aiding efforts to track viral evolution. While current subtyping tools - Genome Detective, GLUE, and NextClade - rely on computationally intensive sequence alignment and phylogenetic inference, machine learning presents a promising alternative for achieving accurate and rapid classification.

Results

We present Craft (Chaos Random Forest), a machine learning framework for dengue subtyping. We demonstrate that Craft is capable of faster classification speeds while matching or surpassing the accuracy of existing tools. Craft achieves 99.5% accuracy on a hold-out test set and processes over 140 000 sequences per minute. Notably, Craft maintains remarkably high accuracy even when classifying sequence segments as short as 700 nucleotides.

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

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

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

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

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

Background

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

Methods

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

Results

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

Conclusions

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

Examining the molecular clock hypothesis for the contemporary evolution of the rabies virus

The molecular clock hypothesis assumes that mutations accumulate on an organism's genome at a constant rate over time, but this assumption does not always hold true. While modelling approaches exist to accommodate deviations from a strict molecular clock, assumptions about rate variation may not fully represent the underlying evolutionary processes. There is considerable variability in rabies virus (RABV) incubation periods, ranging from days to over a year, during which viral replication may be reduced. This prompts the question of whether modelling RABV on a per infection generation basis might be more appropriate. We investigate how variable incubation periods affect root-to-tip divergence under per-unit time and per-generation models of mutation. Additionally, we assess how well these models represent root-to-tip divergence in time-stamped RABV sequences. We find that at low substitution rates (<1 substitution per genome per generation) divergence patterns between these models are difficult to distinguish, while above this threshold differences become apparent across a range of sampling rates. Using a Tanzanian RABV dataset, we calculate the mean substitution rate to be 0.17 substitutions per genome per generation. At RABV's substitution rate, the per-generation substitution model is unlikely to represent rabies evolution substantially differently than the molecular clock model when examining contemporary outbreaks; over enough generations for any divergence to accumulate, extreme incubation periods average out. However, measuring substitution rates per-generation holds potential in applications such as inferring transmission trees and predicting lineage emergence.

A new lineage nomenclature to aid genomic surveillance of dengue virus

Dengue virus (DENV) is currently causing epidemics of unprecedented scope in endemic settings and expanding to new geographical areas. It is therefore critical to track this virus using genomic surveillance. However, the complex patterns of viral genomic diversity make it challenging to use the existing genotype classification system. Here, we propose adding 2 sub-genotypic levels of virus classification, named major and minor lineages. These lineages have high thresholds for phylogenetic distance and clade size, rendering them stable between phylogenetic studies. We present assignment tools to show that the proposed lineages are useful for regional, national, and subnational discussions of relevant DENV diversity. Moreover, the proposed lineages are robust to classification using partial genome sequences. We provide a standardized neutral descriptor of DENV diversity with which we can identify and track lineages of potential epidemiological and/or clinical importance. Information about our lineage system, including methods to assign lineages to sequence data and propose new lineages, can be found at: dengue-lineages.org.

Genomic Characterization of a Dog-Mediated Rabies Outbreak in El Pedregal, Arequipa, Peru

Background

Rabies, a re-emerging zoonosis with the highest known human case fatality rate, has been largely absent from Peru, except for endemic circulation in the Puno region on the Bolivian border and re-emergence in Arequipa City in 2015, where it has persisted. In 2021, an outbreak occurred in the rapidly expanding city of El Pedregal near Arequipa, followed by more cases in 2022 after nearly a year of epidemiological silence. While currently under control, questions persist regarding the origin of the El Pedregal outbreak and implications for maintaining rabies control in Peru.

Methods

We sequenced 25 dog rabies virus (RABV) genomes from the El Pedregal outbreak (n=11) and Arequipa City (n=14) from 2021-2023 using Nanopore sequencing in Peru. Historical genomes from Puno (n=4, 2010-2012) and Arequipa (n=5, 2015-2019), were sequenced using an Illumina approach in the UK. In total, 34 RABV genomes were analyzed, including archived and newly obtained samples. The genomes were analyzed phylogenetically to understand the outbreak's context and origins.

Results

Phylogenomic analysis identified two genetic clusters in El Pedregal: 2021 cases stemmed from a single introduction unrelated to Arequipa cases, while the 2022 sequence suggested a new introduction from Arequipa rather than persistence. In relation to canine RABV diversity in Latin America, all new sequences belonged to a new minor clade, Cosmopolitan Am5, sharing relatives from Bolivia, Argentina, and Brazil.

Conclusion

Genomic insights into the El Pedregal outbreak revealed multiple introductions over a 2-year window. Eco-epidemiological conditions, including migratory worker patterns, suggest human-mediated movement drove introductions. Despite outbreak containment, El Pedregal remains at risk of dog-mediated rabies due to ongoing circulation in Arequipa, Puno, and Bolivia. Human-mediated movement of dogs presents a major risk for rabies re-emergence in Peru, jeopardizing regional dog-mediated rabies control. Additional sequence data is needed for comprehensive phylogenetic analyses.

Geographic Distribution of Rabies Virus and Genomic Sequence Alignment of Wild and Vaccine Strains, Kenya

Rabies, a viral disease that causes lethal encephalitis, kills ≈59,000 persons worldwide annually, despite availability of effective countermeasures. Rabies is endemic in Kenya and is mainly transmitted to humans through bites from rabid domestic dogs. We analyzed 164 brain stems collected from rabid animals in western and eastern Kenya and evaluated the phylogenetic relationships of rabies virus (RABV) from the 2 regions. We also analyzed RABV genomes for potential amino acid changes in the vaccine antigenic sites of nucleoprotein and glycoprotein compared with RABV vaccine strains commonly used in Kenya. We found that RABV genomes from eastern Kenya overwhelmingly clustered with the Africa-1b subclade and RABV from western Kenya clustered with Africa-1a. We noted minimal amino acid variances between the wild and vaccine virus strains. These data confirm minimal viral migration between the 2 regions and that rabies endemicity is the result of limited vaccine coverage rather than limited efficacy.

A new lineage nomenclature to aid genomic surveillance of dengue virus

Dengue virus (DENV) is currently causing epidemics of unprecedented scope in endemic settings and expanding to new geographical areas. It is therefore critical to track this virus using genomic surveillance. However, the complex patterns of viral genomic diversity make it challenging to use the existing genotype classification system. Here we propose adding two sub-genotypic levels of virus classification, named major and minor lineages. These lineages have high thresholds for phylogenetic distance and clade size, rendering them stable between phylogenetic studies. We present an assignment tool to show that the proposed lineages are useful for regional, national and sub-national discussions of relevant DENV diversity. Moreover, the proposed lineages are robust to classification using partial genome sequences. We provide a standardized neutral descriptor of DENV diversity with which we can identify and track lineages of potential epidemiological and/or clinical importance. Information about our lineage system, including methods to assign lineages to sequence data and propose new lineages, can be found at: dengue-lineages.org.

Genetic analyses of rabies virus glycoprotein and nucleoprotein gene sequences reveal the emergence of multiple lineages in animals in Arkhangai province, a central region of Mongolia

Genetic characterizations of rabies viruses circulating in carnivore and non-carnivore animals were investigated for the first time in Arkhangai province, a central region of Mongolia. Also, glycoprotein gene of the rabies virus was sequenced for the first time in Mongolia. The nucleotide sequences of the glycoprotein and nucleoprotein genes were analysed, revealing the presence of multiple lineages in this area. Of particular concern are the lineages identified in carnivores, which might emerge to spread throughout Mongolia, further facilitating transboundary transmission to neighbouring countries, including China and Russia.

Analyzing the Evolution and Host Adaptation of the Rabies Virus from the Perspective of Codon Usage Bias

Rabies virus (RABV) is a highly pathogenic virus that causes a fatal disease in humans and other mammals, but the mechanism of its evolution, spread, and spillover remains unknown. In this study, we analyzed the codon usage pattern of 2,018 RABV full-length genome sequences from 79 countries collected between 1931 and 2021 to provide an insight into its molecular evolution and unravel its unknown host-adapted pattern. We found that RABV exhibited a weak codon usage bias, with a preference for the codons ending in A (28.10 ± 0.01) or U (26.43 ± 0.02). Moreover, natural selection plays a major role in shaping the codon usage bias of the RABV. Notably, nearly half of the 18 codons in the virus were best matched to the hosts' most abundant isoacceptor tRNAs, which might account for the wide range of RABV hosts. Furthermore, significant differences were observed in the codon usage patterns of RABV for different host species, suggesting that codon usage bias may be influenced by host-specific factors. In conclusion, our study reveals codon usage patterns of RABV that may help in the development of control strategies and effective vaccines and therapies against this deadly virus.

Evolution and distribution of rabies viruses from a panorama view

Rabies kills more than 59,000 people annually, mainly in developing countries. Previous studies on the evolution and distribution of rabies viruses (RABVs) were scattered. Here, we explore the evolution and distribution of this deadly virus from a novel panorama view. Multiple bioinformatic software tools were employed to analyze the phylogenetic diversity, evolution, spatiotemporal, and distribution of RABVs. The analyses were based on 1,202 qualified full-length genomes of RABVs and numerous literatures. Of the 10 distinct phylogenetic clades of RABV that we identified, more frequent intra- and inter-clade recombination occurs in the sequences of Asian-SEA, Arctic, and Cosmopolitan clades isolated from China, while according to existing sequence information, RABV might originate from bats (posterior probability, PP = 0.75, PP = 0.60 inferred from N and L genes, separately) in North America (PP = 0.57, PP = 0.62 inferred from N and L genes, separately). Due to the difference in evolutionary rate of N (2.22 × 10-4 subs/site/year, 95% HPD 1.99-2.47 × 10-4 subs/site/year) and L genes (1.67 × 10-4 subs/site/year, 95% HPD 1.59-1.74 × 10-4 subs/site/year), the root age was 1,406.6 (95% HPD 1,291.2-1,518.2) and 1,122.7 (95% HPD 1,052.4-1,193.9) inferred from N and L genes, separately. Among other findings, Mephitidae plays an important role in the interspecific transmission and communication of RABV, which we found tends to spread to populations genetically proximate to the host. We also identified amino acids under positive selection in different genes of different clades as well as single nucleotide variation sites important for different lineages. IMPORTANCE Rabies virus is widely distributed all over the world, and wild animals are its largest potential reservoir. Our study offers a panorama view about evolution and distribution of rabies viruses and emphasizes the need to monitor the transmission dynamics of animal rabies.

Uncovering the endemic circulation of rabies in Cambodia

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

Genomic Surveillance of Rabies Virus in Georgian Canines

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

Emerging Genomic Trends on Rabies Virus in Davao Region, Philippines, 2018-2021

Rabies, caused by the rabies virus (RABV), remains a significant public health issue in the Philippines despite efforts to control it. To eliminate rabies by 2030, effective surveillance strategies are crucial. In this study, we examined RABV evolution and phylodynamics in the Davao Region using genome sequences from Davao City and nearby provinces. We adapted the RABV ARTIC Protocol for Oxford Nanopore High-Throughput Sequencing to optimize workflow efficiency under limited resources. Comparing new virus samples collected from June 2019 to June 2021 (n = 38) with baseline samples from June 2018 to May 2019 (n = 49), new sub-clades were observed in the phylogenetic tree, suggesting divergence from older variants that were previously undetected. Most of the new viruses belonged to the Asian SEA4_A1.1.1 lineage, but new (SEA4_B1 and SEA4_B1.1) and emerging (SEA4_B1.1_E1) lineages that have never been reported in the Philippines were also identified. The baseline study reported phylogeographic clustering of RABV isolates from the same areas. However, this pattern was disrupted in the current biosurveillance, with variants detected in areas outside the original cluster. Furthermore, our findings revealed significant transmission routes between Davao City and neighboring provinces, contrasting with the predominantly intra-city transmission observed in the baseline study. These results underscore the need for ongoing and timely genomic surveillance to monitor genetic diversity changes and the emergence of novel strains, as well as to track alterations in transmission pathways. Implementing cost-effective next-generation sequencing workflows will facilitate the integration of genomic surveillance into rabies control programs, particularly in resource-limited settings. Collaborations between different sectors can empower local laboratories and experts in genomic technologies and analysis.

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

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

Molecular Analysis of Rabies Virus Using RNA Extracted from Used Lateral Flow Devices

Molecular analysis of rabies virus can provide accurate diagnosis and information on its genetic diversity. The transportation of rabies brain samples from remote areas to a central laboratory is challenging owing to biohazard risks and decomposability. We investigated the utility of used lateral flow devices (LFDs) for subsequent molecular analysis and assessed the necessary storage temperatures. Using RNA extracted from used LFD strips, we performed conventional reverse transcription-PCR (RT-PCR) using an LN34 primer set to amplify short fragments (165 bp) for rabies virus detection and the P1-304 primer set to amplify long fragments of the entire N gene amplicon (1,506 bp) for phylogenetic analysis. Among 71 used LFDs stored in a refrigerator and 64 used LFDs stored at room temperature, the LN34 assay showed high sensitivities (96.2% and 100%, respectively) for the diagnosis of rabies, regardless of the storage temperature. A significant reduction in the sensitivity of rabies diagnosis was observed when using the P1-304 primer set for used LFDs stored at room temperature compared to those stored at refrigeration temperature (20.9% versus 100%; P < 0.05). Subsequent sequencing and phylogenetic analysis were successfully performed using the amplicons generated by the P1-304 RT-PCR assays. Used LFDs are thus promising resources for rabies virus RNA detection and sequence analysis. Virus detection via RT-PCR, amplifying a short fragment, was possible regardless of the storage temperature of the used LFDs. However, refrigerated storage is recommended for RT-PCR amplification of long fragments for phylogenetic analysis.

An advanced sequence clustering and designation workflow reveals the enzootic maintenance of a dominant West Nile virus subclade in Germany

West Nile virus (WNV) is the most widespread arthropod-borne (arbo) virus and the primary cause of arboviral encephalitis globally. Members of WNV species genetically diverged and are classified into different hierarchical groups below species rank. However, the demarcation criteria for allocating WNV sequences into these groups remain individual and inconsistent, and the use of names for different levels of the hierarchical levels is unstructured. In order to have an objective and comprehensible grouping of WNV sequences, we developed an advanced grouping workflow using the 'affinity propagation clustering' algorithm and newly included the 'agglomerative hierarchical clustering' algorithm for the allocation of WNV sequences into different groups below species rank. In addition, we propose to use a fixed set of terms for the hierarchical naming of WNV below species level and a clear decimal numbering system to label the determined groups. For validation, we applied the refined workflow to WNV sequences that have been previously grouped into various lineages, clades, and clusters in other studies. Although our workflow regrouped some WNV sequences, overall, it generally corresponds with previous groupings. We employed our novel approach to the sequences from the WNV circulation in Germany 2020, primarily from WNV-infected birds and horses. Besides two newly defined minor (sub)clusters comprising only three sequences each, Subcluster 2.5.3.4.3c was the predominant WNV sequence group detected in Germany from 2018 to 2020. This predominant subcluster was also associated with at least five human WNV infections in 2019-20. In summary, our analyses imply that the genetic diversity of the WNV population in Germany is shaped by enzootic maintenance of the dominant WNV subcluster accompanied by sporadic incursions of other rare clusters and subclusters. Moreover, we show that our refined approach for sequence grouping yields meaningful results. Although we primarily aimed at a more detailed WNV classification, the presented workflow can also be applied to the objective genotyping of other virus species.

Comparative analysis of genome-encoded viral sequences reveals the evolutionary history of flavivirids (family Flaviviridae)

Flavivirids (family Flaviviridae) are a group of positive-strand ribonucleic acid (RNA) viruses that pose serious risks to human and animal health on a global scale. Here, we use flavivirid-derived deoxyribonucleic acid (DNA) sequences, identified in animal genomes, to reconstruct the long-term evolutionary history of family Flaviviridae. We demonstrate that flavivirids are >100 million years old and show that this timing can be combined with dates inferred from co-phyletic analysis to produce a cohesive overview of their evolution, distribution, and diversity wherein the main flavivirid subgroups originate in early animals and broadly co-diverge with major animal phyla. In addition, we reveal evidence that the 'classical flaviviruses' of vertebrates, most of which are transmitted via blood-feeding arthropod vectors, originally evolved in haematophagous arachnids and later acquired the capacity to be transmitted by insects. Our findings imply that the biological properties of flavivirids have been acquired gradually over the course of animal evolution. Thus, broad-scale comparative analysis will likely reveal fundamental insights into their biology. We therefore published our results via an open, extensible, database (Flavivirid-GLUE), which we constructed to facilitate the wider utilisation of genomic data and evolution-related domain knowledge in flavivirid research.