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

Emergence of Arctic-like Rabies Lineage in India

Progenitors of Arctic-like rabies viruses, which now circulate extensively in India, may have been responsible for the emergence of the Arctic rabies lineage

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

Genetic Tracking of the Raccoon Variant of Rabies Virus in Eastern North America

AbstractTo gain insight into the incursion of the raccoon variant of rabies into the raccoon population in three Canadian provinces, a collection of 192 isolates of the raccoon rabies virus (RRV) strain was acquired from across its North American range and was genetically characterized. A 516-nucleotide segment of the non-coding region between the G and L protein open reading frames, corresponding to the most variable region of the rabies virus genome, was sequenced. This analysis identified 119 different sequences, and phylogenetic analysis of the dataset supports the documented history of RRV spread. Three distinct geographically restricted RRV lineages were identified. Lineage 1 was found in Florida, Alabama and Georgia and appears to form the ancestral lineage of the raccoon variant of rabies. Lineage 2, represented by just two isolates, was found only in Florida, while the third lineage appears broadly distributed throughout the rest of the eastern United States and eastern Canada. In New York State, two distinct spatially segregated variants were identified; the one occupying the western and northern portions of the state was responsible for an incursion of raccoon rabies into the Canadian province of Ontario. Isolates from New Brunswick and Quebec form distinct, separate clusters, consistent with their independent origins from neighboring areas of the United States. The data are consistent with localized northward incursion into these three separate areas with no evidence of east–west viral movement between the three Canadian provinces

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

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

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

Influence of genomic variations on glanders serodiagnostic antigens using integrative genomic and transcriptomic approaches

Glanders is a highly contagious and life-threatening zoonotic disease caused by Burkholderia mallei (B. mallei). Without an effective vaccine or treatment, early diagnosis has been regarded as the most effective method to prevent glanders transmission. Currently, the diagnosis of glanders is heavily reliant on serological tests. However, given that markedly different host immune responses can be elicited by genetically different strains of the same bacterial species, infection by B. mallei, whose genome is unstable and plastic, may result in various immune responses. This variability can make the serodiagnosis of glanders challenging. Therefore, there is a need for a comprehensive understanding and assessment of how B. mallei genomic variations impact the appropriateness of specific target antigens for glanders serodiagnosis. In this study, we investigated how genomic variations in the B. mallei genome affect gene content (gene presence/absence) and expression, with a special focus on antigens used or potentially used in serodiagnosis. In all the genome sequences of B. mallei isolates available in NCBI’s RefSeq database (accessed in July 2023) and in-house sequenced samples, extensive small and large variations were observed when compared to the type strain ATCC 23344. Further pan-genome analysis of those assemblies revealed variations of gene content among all available genomes of B. mallei. Specifically, differences in gene content ranging from 31 to 715 genes with an average of 334 gene presence-absence variations were found in strains with complete or chromosome-level genome assemblies, using the ATCC 23344 strain as a reference. The affected genes included some encoded proteins used as serodiagnostic antigens, which were lost due mainly to structural variations. Additionally, a transcriptomic analysis was performed using the type strain ATCC 23344 and strain Zagreb which has been widely utilized to produce glanders antigens. In total, 388 significant differentially expressed genes were identified between these two strains, including genes related to bacterial pathogenesis and virulence, some of which were associated with genomic variations, particularly structural variations. To our knowledge, this is the first comprehensive study to uncover the impacts of genetic variations of B. mallei on its gene content and expression. These differences would have significant impacts on host innate and adaptive immunity, including antibody production, during infection. This study provides novel insights into B. mallei genetic variants, knowledge which will help to improve glanders serodiagnosis

Early Cenozoic denudation of central west Britain in response to transient and permanent uplift above a mantle plume

Upwelling mantle plumes beneath continental crust are predicted to produce difficult to quantify, modest uplift and denudation. The contribution of permanent and transient components to the uplift is also difficult to distinguish. A pulse of denudation in Britain in the Early Paleogene has been linked, although with some controversy, with the arrival of the proto-Iceland mantle plume. In this contribution we show that combining apatite and zircon (U-Th-Sm)/He and apatite fission track analyses from central west Britain with numerical modeling clearly identifies a pulse of early Cenozoic denudation. The data indicate that rock uplift and denudation were centered on the northern East Irish Sea Basin and 1.0–2.4 km of rocks were removed during the latest Cretaceous-early Paleogene. Uplift and erosion appears to have started a few million years before the earliest magmatism in the region. The regional denudation pattern mirrors the distribution of low-density magmatic rocks that has been imaged in the deep crust. However, the injection of the underplating melt is not enough to account for the total denudation. An additional regional uplift of at least 300 m is required, which is consistent with a transient thermal effect from the hot mantle plume. The rapid exhumation event ceased by ~40 Ma and the data do not require significant Neogene exhumation