key: cord-0790330-1xo4uc8e
authors: Kumar, Deepak; Chaudhary, Suman; Lu, Nanyan; Duff, Michael; Heffel, Mathew; McKinney, Caroline A.; Bedenice, Daniela; Marthaler, Douglas
title: Metagenomic Next-Generation Sequencing Reveal Presence of a Novel Ungulate Bocaparvovirus in Alpacas
date: 2019-07-31
journal: Viruses
DOI: 10.3390/v11080701
sha: 8b1ef7f0f343414329cefb9a014bffcf3f92f117
doc_id: 790330
cord_uid: 1xo4uc8e

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.

Bocaparvoviruses (BoVs) belong to the genus Bocaparvovirus and are emerging pathogens of the Parvoviridae family. BoVs are nonenveloped, single-stranded DNA viruses with an icosahedral symmetry and were originally named according to their first identified members, bovine parvovirus (BPV) and minute virus of canine (MVC) [1] . In the past few years, novel BoVs have been identified in a variety of animals, including bats [2] , camels [3] , gorillas [4] , marmots [5] , pigs [6] , and rodents [7] . BoVs are comprised of 21 species, including carnivore BoV 1-6, chiropteran BoV 1-4, lagomorph BoV 1, pinniped BoV 1 and 2, primate BoV 1 and 2, and ungulate BoV (UBoV) 1-6. A few new UBoVs have been identified in dromedary camels (tentatively UBoV7 and UBoV8) [3] but have yet to be classified by the International Committee on Taxonomy of Viruses (ICTV).

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] .

Alpaca (Vicugna pacos) are domesticated members of the new world camelids closely related to llama (Lama glama), guanaco (Lama guanicoe), and vicuna (Vicugna vicugna) [13] . Over the past couple of decades, alpacas have gained significant popularity as pets, show animals, and fiber animals in the United States, with a total of 264,587 alpacas registered in the US as of May 2019 [14] . A variety of viruses have been identified in alpacas, including bovine viral diarrhea virus, coronavirus, adenovirus, equine viral arteritis virus, rotavirus, rabies, bluetongue virus, foot-and-mouth disease virus, bovine respiratory syncytial virus, influenza A virus, bovine papillomavirus, vesicular stomatitis virus, parainfluenza-3 virus, West Nile virus, and equine herpesvirus [12, [15] [16] [17] [18] [19] [20] . However, BoVs have yet to be reported in alpacas.

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. 

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

Genome Sequence of a Porcine Bocavirus Detected in Feces of Domestic Minks in China

Identification and Interspecies Transmission of a Novel Bocaparvovirus among Different Bat Species in China

Two novel dromedary camel bocaparvoviruses from dromedaries in the Middle East with unique genomic features

Identification and characterization of a new Bocavirus species in Gorillas

Two novel bocaparvovirus species identified in wild Himalayan marmots

Detection and characterisation of novel bocavirus (genus Bocaparvovirus) and gastroenteritis viruses from asymptomatic pigs in Ireland

Identification and characterization of a novel rodent bocavirus from different rodent species in China

Identification and nearly full-lengthg characterization of novel porcine Bocaviruses

Bocavirus infection in otherwise healthy children with respiratory disease

Structure of an enteric pathogen, bovine parvovirus

Recovery and characterization of a minute virus of canines

Discovery of a Bovine Enterovirus in Alpaca

Alpaca Owners Association Inc Homepage

Identification and characterization of a novel alpaca respiratory coronavirus most closely related to the human coronavirus 229E

Variation in Bluetongue virus real-time reverse transcription polymerase chain reaction assay results in blood samples of sheep, cattle, and alpaca

Viral diseases of new world camelids

Analysis of the genome sequence of an alpaca coronavirus

Biosecurity and biocontainment in alpaca operations

Foot-and-mouth disease in camelids: A review

Bovine rhinitis viruses are common in U.S. cattle with bovine respiratory disease

Simultaneous rapid sequencing of multiple RNA virus genomes

A porcine enterovirus G associated with enteric disease contains a novel papain-like cysteine protease

Trimmomatic: A flexible trimmer for Illumina sequence data

Sickle-A Windowed Adaptive Trimming Tool for FASTQ Files Using Quality

Ultrafast metagenomic sequence classification using exact alignments

MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability

Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data

RDP4: Detection and analysis of recombination patterns in virus genomes

High rate of viral evolution associated with the emergence of carnivore parvovirus

NP1 protein of the bocaparvovirus minute virus of canines controls access to the viral capsid genes via its role in RNA processing

A Viral Phospholipase A2 Is Required for Parvovirus Infectivity

The determinants for the enzyme activity of human parvovirus b19 phospholipase a2 (pla2) and its influence on cultured cells

The metagenomic approach and causality in virology

The authors would like to thank the Kansas State Veterinary Diagnostic Laboratory and Tufts University for supporting the project.

The authors declare no conflict of interest.