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.

A subpopulation of the arctic fox lineage of rabies virus has circulated extensively in red fox populations of Ontario, Canada, between the 1960s and 1990s. An intensive wildlife rabies control program, in which field operations were initiated in 1989, resulted in elimination of the disease in eastern Ontario. However in southwestern Ontario, as numbers of rabid foxes declined the proportion of skunks confirmed to be infected with this rabies virus variant increased and concerted control efforts targeting this species were employed to eliminate the disease. Since 2012 no cases due to this viral variant were reported in southwestern Ontario until 2015 when a single case of rabies due to the arctic fox variant was reported in a bovine. Several additional cases have been documented subsequently. Since routine antigenic typing cannot discriminate between the variants which previously circulated in Ontario and those from northern Canada it was unknown whether these recent cases were the result of a new introduction of this variant or a continuation of the previous enzootic. To explore the origins of this new outbreak whole genome sequences of a collection of 128 rabies viruses recovered from Ontario between the 1990s to the present were compared with those representative of variants circulating in the Canadian north. Phylogenetic analysis shows that the variant responsible for current cases in southwestern Ontario has evolved from those variants known to circulate in Ontario previously and is not due to a new introduction from northern regions. Thus despite ongoing passive surveillance the persistence of wildlife rabies went undetected in the study area for almost three years. The apparent adaptation of this rabies virus variant to the skunk host provided the opportunity to explore coding changes in the viral genome which might be associated with this host shift. Several such changes were identified including a subset for which the operation of positive selection was supported. The location of a small number of these amino acid substitutions in or close to protein motifs of functional importance suggests that some of them may have played a role in this host shift. Rabies, a serious disease which almost invariably results in death once clinical signs appear, is caused by the rabies virus. The arctic fox lineage of rabies virus persists in fox populations in northern Canada and from the 1960s onwards a sub-population of this virus has circulated in red foxes in the province of Ontario, Canada. By 2000 a provincial wildlife control program attempting to control fox rabies in southern Ontario had eliminated the disease in most of the targeted regions except for a focus of rabies in certain counties in southwestern Ontario which involved mostly skunks. Modified control efforts resulted in apparent rabies elimination with no cases reported in southern Ontario after May 2012 until a rabid bovine was identified in December 2015. Additional cases, mostly in skunks and livestock, have been reported subsequently. To discriminate between a re-introduction of the virus from the north and re-emergence of the Ontario viral variant, we sequenced and compared the genome of multiple viruses from both areas. The sequence data indicate that sub-types of the Ontario variant which eluded control efforts are responsible for these recent cases. These data also permit an examination of potential genetic changes in the virus that may facilitate its circulation in this new wildlife host, the skunk.

Rabies spreads in both Arctic (Vulpes lagopus) and red foxes (Vulpes vulpes) throughout the Canadian Arctic but limited wildlife disease surveillance, due to the extensive landmass of the Canadian north and its small widely scattered human population, undermines our knowledge of disease transmission patterns. This study has explored genetic population structure in both the rabies virus and its fox hosts to better understand factors that impact rabies spread. Phylogenetic analysis of 278 samples of the Arctic lineage of rabies virus recovered over 40 years identified four sub-lineages, A1 to A4. The A1 lineage has been restricted to southern regions of the Canadian province of Ontario. The A2 lineage, which predominates in Siberia, has also spread to northern Alaska while the A4 lineage was recovered from southern Alaska only. The A3 sub-lineage, which was also found in northern Alaska, has been responsible for virtually all cases across northern Canada and Greenland, where it further differentiated into 18 groups which have systematically evolved from a common predecessor since 1975. In areas of Arctic and red fox sympatry, viral groups appear to circulate in both hosts, but both mitochondrial DNA control region sequences and 9-locus microsatellite genotypes revealed contrasting phylogeographic patterns for the two fox species. Among 157 Arctic foxes, 33 mitochondrial control region haplotypes were identified but little genetic structure differentiating localities was detected. Among 162 red foxes, 18 control region haplotypes delineated three groups which discriminated among the Churchill region of Manitoba, northern Quebec and Labrador populations, and the coastal Labrador locality of Cartwright. Microsatellite analyses demonstrated some genetic heterogeneity among sampling localities of Arctic foxes but no obvious pattern, while two or three clusters of red foxes suggested some admixture between the Churchill and Quebec-Labrador regions but uniqueness of the Cartwright group. The limited population structure of Arctic foxes is consistent with the rapid spread of rabies virus subtypes throughout the north, while red fox population substructure suggests that disease spread in this host moves most readily down certain independent corridors such as the northeastern coast of Canada and the central interior. Interestingly the evidence suggests that these red fox populations have limited capacity to maintain the virus over the long term, but they may contribute to viral persistence in areas of red and Arctic fox sympatry.