Next-generation sequencing approach to investigate genome variability of Parapoxvirus in Canadian muskoxen (Ovibos moschatus)

In 2016, the first orf virus, a double-stranded DNA (dsDNA) virus of the genus parapoxvirus, from a muskox was isolated on Victoria Island, Nunavut (NU), Canada. We used deep sequencing on DNA extracted from orf virus -positive tissues from wild muskoxen from locations on Victoria Island and the adjacent mainland. Orf virus sequence reads derived from four samples were nearly identical. The consensus sequences generated from pooled reads of MxOV comprises of a large contiguous sequence (contig) of 131,759 bp and a smaller right terminal contig of 3552 bp, containing all coding sequences identified as Parapoxvirus. Individual gene comparisons reveal that MxOV shares genetic characteristics with reference strains from both sheep and goat origin. Recombination analysis using Bootscan, MAXCHI, GENECONV, CHIMAERA, SISCAN, and RDP algorithms within the RDP4 software predicted recombination events in two virulence factors, and a large 3000 bp segment of the MxOV genome. Partial B2L nucleotide sequences from strains around the world and other North American isolates were compared to MxOV using MUSCLE alignments and RAxML phylogenetic trees. MxOV was identical to our previously characterized isolate, and shared similarity with orf virus isolated from sheep and goats. The phylogenetic grouping of partial B2L nucleotide sequences did not follow the sample geographic distribution. More full genomes of orf virus, or at least full B2L gene squences, in wildlife are needed especially in North America to better understand the epidemiology of the disease in muskoxen.