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

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.

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.