Metagenomic Next-Generation Sequencing Reveal Presence of a Novel Ungulate Bocaparvovirus in Alpacas

Viruses belonging to the genus Bocaparvovirus (BoV) are a genetically diverse group of DNA viruses known to cause respiratory, enteric, and neurological diseases in animals, including humans. An intestinal sample from an alpaca (Vicugna pacos) herd with reoccurring diarrhea and respiratory disease was submitted for next-generation sequencing, revealing the presence of a BoV strain. The alpaca BoV strain (AlBoV) had a 58.58% whole genome nucleotide percent identity to a camel BoV from Dubai, belonging to a tentative ungulate BoV 8 species (UBoV8). Recombination events were lacking with other UBoV strains. The AlBoV genome was comprised of the NS1, NP1, and VP1 proteins. The NS1 protein had the highest amino acid percent identity range (57.89–67.85%) to the members of UBoV8, which was below the 85% cut-off set by the International Committee on Taxonomy of Viruses. The low NS1 amino acid identity suggests that AlBoV is a tentative new species. The whole genome, NS1, NP1, and VP1 phylogenetic trees illustrated distinct branching of AlBoV, sharing a common ancestor with UBoV8. Walker loop and Phospholipase A2 (PLA2) motifs that are vital for virus infectivity were identified in NS1 and VP1 proteins, respectively. Our study reports a novel BoV strain in an alpaca intestinal sample and highlights the need for additional BoV research.

Initially, the classification of parvoviruses required the isolation of the virus; however, reporting of the viral sequence containing all the non-structural and structural coding regions is now acceptable provided the genomic, serological, or biological data supports infectious etiology [8]. Most of the members of the Bocaparvovirus genus have been identified using molecular methods and lack isolation in cell culture [3,4,9]. Human BoVs cause severe respiratory and gastrointestinal infections in young children [10]. Bovine parvovirus (BPV) causes gastrointestinal and respiratory symptoms, reproductive Viruses 2019, 11, 701 2 of 7 failure, and conjunctivitis in cattle worldwide [11]. Another important member of the BoV genus, canine minute virus (MVC), causes sub-clinical disease and fetal infections often leading to neonatal respiratory disease or abortions [12]. However, Koch's postulates have yet to be fulfilled to link newly emerging BoVs with the clinical disease in animals [1,3,5].
An alpaca farm in the mid-eastern United States reported recurrent diarrhea and respiratory failure in young alpacas, with a case fatality rate up to 100%. In 2017, an alpha coronavirus was identified as causing clinical disease in two animals, and vaccination was subsequently attempted. However, diarrhea and respiratory distress continued to occur in juvenile animals despite increased biosecurity measures and supportive herd management. In 2018, an intestinal sample from a deceased alpaca was submitted to Kansas State University Veterinary Diagnostic Laboratory for metagenomic next-generation sequencing (NGS) to further evaluate potential causes of disease. The intestinal sample was processed, extracted, and sequenced using previously described methods [21,22]. The raw data was analyzed using a custom bioinformatic pipeline [23]. Reads were trimmed, and the adapter/index sequences were removed using Trimmomatic [24], Sickle [25], and Scythe [26].
A total of 334,052 cleaned reads were classified as eukaryotes (41%), bacteria (28%), viral (7%), and other organisms (4%) by Kraken software, which applies a k-mer search strategy from a sequence database to taxonomically classify reads ( Figure 1) [27]. Kraken revealed a majority of the viral reads (22,170) as BoV (77%); bacteriophages (14%); and miscellaneous viruses composed of retroviruses, bacteria viruses, and unclassified viruses (9%). Reads lacking classification (no hits, n = 67,604) and identified as viral reads (n = 22,170) were de novo assembled into contigs and BLAST (Basic Local Alignment Search Tool) against the National Center for Biotechnology Information (NCBI) database, identifying a contig with a 58.58% nucleotide percent identity to a camel BoV from Dubai (KY640435). A full-length genome (5155 nucleotides) of an alpaca BoV (AlBoV, GenBank number MK014742) had an average read coverage of 2440X. AlBoV was aligned with the 108 complete UBoV genomes from GenBank using Multiple Alignment using Fast Fourier Transform (MAFFT) [28] in Geneious Prime [29]. AlBoV shared a 57.77-58.58% whole genome nucleotide identity to the UBoV8 strains (Table 1). Recombination events were not detected in the UBoV alignment using RDP4 software [30], although single-stranded DNA viruses such as parvoviruses possess a mutation rate similar to single-stranded RNA viruses and a higher mutation rate than double-stranded DNA viruses [31].
AlBoV contained three open reading frames (ORFs), NS1, NP1, and VP1/VP2, which were 2154 bp (411 to 2564 bp), 507 bp (2701 to 3207 bp), and 1395bp (3194 to 4588 bp), respectively. ICTV indicates a new parvovirus species should have less than 85% amino acid identity of the NS1 protein with other parvovirus species. The AlBoV identified in the present study shared the highest NS1 amino acid percent identity (57.89-67.85%) with camel BoVs in UBoV8 (Table 1) and represents a tentative new BoV species, UBoV9. The ancillary protein NP1, which is a unique feature of BoVs, is known to influence RNA processing events by suppressing internal polyadenylation and splicing of an upstream intron [32]. Unlike some of the other BoV sequences, the coding region of NP1 of AlBoV did not overlap with the C-terminal region of NS1. Interestingly, VP1/VP2 gene of AlBoV was the shortest among the identified UBoVs in the GenBank ( Figure 2). Frameshift mutation were lacking in the AlBoV VP1/VP2 sequence, and the nucleotide sequence after the VP1/VP2 stop codon varied among the 108 complete UBoV sequences in the GenBank.   To study the phylogenetic relationship between AlBoV and other UBoVs, whole genome, NS1, NP1, and VP1 phylogenetic trees were created using a maximum likelihood method (phyML), using 500 bootstrap replicates in Geneious Prime. The trees were curated in FigTree (available from http://tree.bio.ed.ac.uk/software/figtree/) and Adobe Illustrator CS6 (Adobe Systems Inc, San Jose, CA, USA). Whole genome and NS1 phylogenetic trees illustrated that AlBoV shared a common ancestor with the UBoV8 species from camels ( Figure 3). All eight species of UBoV (1-8) illustrated clear grouping in phylogenetic trees, which was observed in NP1 and VP1 phylogenetic trees as well.  the AlBoV VP1/VP2 sequence, and the nucleotide sequence after the VP1/VP2 stop codon varied among the 108 complete UBoV sequences in the GenBank. a b   To study the phylogenetic relationship between AlBoV and other UBoVs, whole genome, NS1, NP1, and VP1 phylogenetic trees were created using a maximum likelihood method (phyML), using 500 bootstrap replicates in Geneious Prime. The trees were curated in FigTree (available from http://tree.bio.ed.ac.uk/software/figtree/) and Adobe Illustrator CS6 (Adobe Systems Inc, San Jose, CA, USA). Whole genome and NS1 phylogenetic trees illustrated that AlBoV shared a common ancestor with the UBoV8 species from camels ( Figure 3). All eight species of UBoV (1-8) illustrated clear grouping in phylogenetic trees, which was observed in NP1 and VP1 phylogenetic trees as well. To study the phylogenetic relationship between AlBoV and other UBoVs, whole genome, NS1, NP1, and VP1 phylogenetic trees were created using a maximum likelihood method (phyML), using 500 bootstrap replicates in Geneious Prime. The trees were curated in FigTree (available from http://tree.bio.ed.ac.uk/software/figtree/) and Adobe Illustrator CS6 (Adobe Systems Inc, San Jose, CA, USA). Whole genome and NS1 phylogenetic trees illustrated that AlBoV shared a common ancestor with the UBoV8 species from camels ( Figure 3). All eight species of UBoV (1-8) illustrated clear grouping in phylogenetic trees, which was observed in NP1 and VP1 phylogenetic trees as well. To investigate and identify the presence of virulence attributes, AlBoV was screened for the ATP or GTP-binding Walker loop motif (GPASTGKT) and Phospholipase A2 (PLA2) motif with the calcium-binding loop and phospholipase catalytic residues; GPASTGKT and PLA2 were found in the NS1 and N-terminal of VP1 proteins, respectively (Figure 4). These protein motifs are conserved and are required for parvovirus infectivity. Phospholipase A2 activity, with the calcium-binding loop and phospholipase catalytic residues, is critical for efficient transfer of the viral genome from the late endosomes/lysosomes to the nucleus for the initiation of replication, and hence is considered essential for virus infectivity [33]. Mutations of critical amino acid residues in the VP1 protein of human parvovirus B19 induces a strong reduction in phospholipase A2 activity and virus infectivity [34]. Considering their vital role in parvovirus infectivity, PLA2 inhibitors are also targeted for antiviral drugs against parvovirus-associated diseases. The presence of the Walker loop and Phospholipase A2 motifs suggests that the newly identified alpaca BoV possesses the virulence determinants necessary to cause disease. To investigate and identify the presence of virulence attributes, AlBoV was screened for the ATP or GTP-binding Walker loop motif (GPASTGKT) and Phospholipase A2 (PLA2) motif with the calcium-binding loop and phospholipase catalytic residues; GPASTGKT and PLA2 were found in the NS1 and N-terminal of VP1 proteins, respectively ( Figure 4). These protein motifs are conserved and are required for parvovirus infectivity. Phospholipase A2 activity, with the calcium-binding loop and phospholipase catalytic residues, is critical for efficient transfer of the viral genome from the late endosomes/lysosomes to the nucleus for the initiation of replication, and hence is considered essential for virus infectivity [33]. Mutations of critical amino acid residues in the VP1 protein of human parvovirus B19 induces a strong reduction in phospholipase A2 activity and virus infectivity [34]. Considering their vital role in parvovirus infectivity, PLA2 inhibitors are also targeted for antiviral drugs against parvovirus-associated diseases. The presence of the Walker loop and Phospholipase A2 motifs suggests that the newly identified alpaca BoV possesses the virulence determinants necessary to cause disease.
for virus infectivity [33]. Mutations of critical amino acid residues in the VP1 protein of human parvovirus B19 induces a strong reduction in phospholipase A2 activity and virus infectivity [34]. Considering their vital role in parvovirus infectivity, PLA2 inhibitors are also targeted for antiviral drugs against parvovirus-associated diseases. The presence of the Walker loop and Phospholipase A2 motifs suggests that the newly identified alpaca BoV possesses the virulence determinants necessary to cause disease. A new species of UBoV was identified in an alpaca intestinal sample. Bocaparvoviruses cause respiratory and gastrointestinal infections in humans, bovines, and other animals. However, the causation of clinical disease in this alpaca farm is unclear. Establishing an association between the presence of BoV and clinical disease would require comprehensive PCR testing of the alpaca farms. Moreover, Koch's postulates is required to establish a virus-disease association. A causal association between the presence of BoV and clinical disease is often difficult due to prolonged viral shedding by the host after infection, high prevalence of BoV infection, and high rate of co-infections [35]. Nevertheless, the discovery of a new BoV will help in developing new PCR diagnostics to determine the prevalence of BoV in alpaca herds and also to develop vaccines to prevent clinical disease. Given the high mutation rate of BoVs and increasing domestication of alpacas, identification of a new BoV in alpaca presents a true risk of cross-species transmission to other mammals.

Funding:
The study was partially supported by Marthaler's start up fund.