key: cord-0000414-782iew86 authors: Kapoor, A; Simmonds, P.; Slikas, B.; Li, Linlin; Bodhidatta, L.; Sethabutr, O.; Triki, H.; Bahri, Olfa; Oderinde, B.; Baba, M.; Bukbuk, D.; Besser, J.; Bartkus, J.; Delwart, E. title: Human bocaviruses are highly diverse, dispersed, recombination prone, and prevalent enteric infections date: 2010-06-01 journal: The Journal of Infectious Diseases DOI: 10.1086/652416 sha: 97082b55ee23b20688f82dac1ba4f03e6d793c85 doc_id: 414 cord_uid: 782iew86 A new species of parvovirus tentatively named human bocavirus 4 (HBoV4) was genetically characterized. Among 641 feces samples from children and adults the most commonly detected bocaviruses species were HBoV2>HBoV3>HBoV4>HBoV1 with HBoV2 prevalence of 21% and 26% in Nigerian and Tunisian children. HBoV3 and HBoV4 species combined were found in 12/192 cases of non-polio acute flaccid paralysis (AFP) from Tunisia and Nigeria and 0/96 healthy Tunisian contacts (p=0.01). Evidence of extensive recombination at the NP1 and VP1 gene boundary between and within species was found. The multiple species and high degree of genetic diversity seen among the human bocaviruses found in feces relative to the highly homogeneous HBoV1 suggest that this world-wide distributed respiratory pathogen may have recently evolved from an enteric bocavirus, perhaps after acquiring an expanded tropism favoring the respiratory track. Elucidating the possible role of the newly identified enteric bocaviruses in human diseases including AFP and diarrhea will require further epidemiological studies. B19 in 1976, several other human parvoviruses have recently been genetically characterized. PARV4 was found in the blood of a febrile adult intravenous drug user [5] , HBoV1 in the nasopharyngeal secretion of a child with respiratory problems [6] , HBoV2 in the stool of children with non-polio acute flaccid paralysis (AFP) [7] , and HBoV3 in the stool of Australian children with diarrhea [8] . In contrast to PARV4 or B19, bocaviruses contain a third open reading frame of unknown function [9, 10] . The first bocavirus identified was in cows [10] , and the name of the genus is derived from its first known hosts (bovine-canine). Animal bocaviruses can cause both respiratory and gastrointestinal diseases, as well as embryonic and fetal death [11] . HBoV1 infection has been linked with mild-to-severe, primarily lower respiratory tract infections in children, frequently in association with other viral infections [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] . HBoV1 has also been detected in stool samples, although association with diarrhea appears weaker than with respiratory symptoms [26] [27] [28] [29] [30] [31] . HBoV1 strains show a very low degree of genetic variability worldwide [32, 33] . In this study, pan-bocavirus polymerase chain reaction (PCR) primers were designed and used on human stool samples from several countries. All 3 recently identified bocaviruses plus a fourth species, HBoV4, were identified. A high degree of genetic diversity, relative to that seen for HBoV1, was seen among human bocaviruses in feces. Stool samples were collected as part of previous clinical studies and were anonymized. This study was approved by the University of California San Francisco Committee on Human Research. Samples from Nigeria and Tunisia were collected as part of the World Health Organization's poliovirus eradication program from children with non-polio AFP between the ages of 4 months and 15 years. Stool samples from healthy contacts of case patients from Tunisia were matched for age. Stool samples from Nepal were obtained from adult travelers and resident expatriates with diarrhea with no known pathogens detected by standard microbiologic testing for enteric bacteria; enzyme immunoassay for rotavirus, adenovirus, astrovirus, Giardia and Cryptosporidium; and reverse transcription PCR for norovirus. Stool samples from healthy, asymptomatic control persons were collected from the same population. Stool samples from the Minnesota Department of Health are from individuals with diarrhea and healthy individuals, matched for age and residential area code. PCR amplification of bocaviruses. We used nested PCR targeting the VP1/2 region of both HBoV1 and HBoV2 (nucleotide positions 3233-3808, numbered here and subsequently with use of the HBoV2 prototype sequence; GenBank accession number FJ170278). Nucleic acids (both DNA and RNA) were extracted (QIAamp Viral RNA mini kit) from 140 mL of clarified stool supernatant and were eluted into 60 mL of water. First-round PCR primers were AK-VP-F1 (5 -CGCCGTGGCT-CCTGCTCT-3 ) and AK-VP-R1 (5 -TGTTCGCCATCACAAA-AGATGTG-3 ) and second-round primers were AK-VP-F2 (5 -GGCTCCTGCTCTAGGAAATAAAGAG-3 ) and AK-VP-R2 (5 -CCTGCTGTTAGGTCGTTGTTGTATGT-3 ). PCR reactions contained 2.5 U of Taq DNA polymerase (NEB) in 1.1ϫ Thermopol reaction buffer with MgCl (2.0 mmol/L), 20 pmol/ L (each) of forward and reverse primers, and 2.5 mL of nucleic acids (first round) or 1 mL of the first-round PCR product (second round) as template in a 50 mL total volume. Firstround conditions were 10 cycles at 95ЊC for 35 s, 58ЊC for 1 min, and 72ЊC for 1 min, with a decrease of 0.5ЊC in annealing temperature each cycle; 30 cycles at 95ЊC for 30 s, 54ЊC for 45 s, and 72ЊC for 45 s; and a final extension at 72ЊC for 10 min. Similar conditions were used for the second PCR round, except that the initial annealing temperatures were 60ЊC and 58ЊC in the first and second group of PCR cycles, respectively. Amplicons of the appropriate size, as detected by agarose gel electrophoresis, were directly sequenced. The PCR products whose sequences are reported here produced unambiguous dideoxysequencing electrophoregram peaks indicating the predominance (190%) of a single bocavirus variant. The sensitivity of the pan-bocavirus nested PCR was determined using dilution of plasmids containing the target sequences from each of the 4 bocavirus species and was estimated at 10-100 plasmid copies for each species. Complete genome sequencing. Nearly complete genomes were amplified using PCR primers designed from alignments of HBoV1 and HBoV2 genomes and were then directly sequenced by primer walking [7] . The terminal sequences were acquired by a modified protocol for rapid amplification of complementary DNA ends [5] . The terminal sequences are incomplete because of extensive hairpin structures that prevented extensions to the viral 5 and 3 extremities. Distance measurements and phylogenetic analyses. Phylogenetic relationships were evaluated with use of Mega 4.1 (http:// www.megasoftware.net/mega41.html) [34] . Neighbor-joining trees were inferred using a matrix of pair-wise maximum-likelihood distances computed from a nucleotide alignment including the genomes obtained in this study and in GenBank, plus a matrix of PAM distances computed from the inferred amino acid alignment. Recombination analyses. Similarity values based on Jukes-Cantor corrected nucleotide distances between full-length sequences were calculated using the program SequenceDist in the Simmonic2005 v1.6 Sequence Editor Package [35] . To assess similarity across the genomes, sequence scans were performed using a fragment length of 300 bases and an increment of 9 bases between fragments. For sequence comparisons with HBoV1 and HBoV2, a mean pair-wise distance was computed using a set of 14 HBoV1 and 7 HBoV2A sequences. Nucleotide sequence accession numbers. The near-full genomes and partial VP1 gene sequences have been deposited in GenBank under accession numbers FJ973558-FJ973563 and GQ506558-GQ506661. Pan-bocavirus PCR primers were designed that could anneal to both HBoV1 and HBoV2, amplifying a ∼576-nucleotide fragment of the VP1 capsid gene. Nucleic acids from human stools collected from Nigeria, Tunisia, Nepal, and the United States were then analyzed using this nested PCR (Table 1) . Of the 641 samples tested, 101 (16%) had confirmed positive results by PCR sequencing, with the highest prevalence in Tunisian children with AFP (33%). To determine the phylogenetic relationship of these strains, the sequences were aligned with available sequences of HBoV1, HBoV2, and HBoV3. Only 4 of 101 strains (from Nigeria and Tunisia) grouped with HBoV1. The remaining strains shared a more recent common ancestry with HBoV2, although the sequences clustered into 4 distinct genetic lineages, labeled HBoV2A, HBoV2B, HBoV3, HBoV4 ( Figure 1 ). The HBoV2A clade included all published HBoV2 genotypes [7] , plus 4 new strains from Nigeria and Tunisia and the recently reported W153 strain from Australia [8] . Because of the high degree of genetic diversity observed and to reduce splitting of HBoV2 into a multitude of genotypes, all strains in that cluster were reclassified into a new, now more diverse genotype A (HBoV2A). The nearly full genome of a new HBoV2A variant was sequenced (TU-C-114-06) and showed very high protein identity (198%) to the Pakistani HBoV2 prototype PK5510 (FJ170278). The HBoV2B cluster included 76 of the 101 bocavirus strains reported here ( Figure 1 ). Two nearly full HBoV2B genome sequences demonstrated a pair-wise amino acid divergence in VP1 of 0.45% (Table 2 ) (NI-213 and NI-327; GenBank accession numbers FJ973560 and FJ973559). Greater divergence was observed when these strains were compared with the VP1 of HBoV2A strains (average divergence, 3.9% [range, 3.3%-4.5%] ( Table 2 ). The nucleotide and amino acid distances at the other loci are shown in Table 2 . The HBoV3 cluster included 11 strains, and the genomes of 2 representative strains, NI-374 and TU-A-210-07, were sequenced (GenBank accession numbers FJ973563 and FJ973562). When compared with HBoV2, the NS1 region of both strains showed an average amino acid divergence of 26% in the NS1 region and 9% in the VP1 region. Compared with HBoV1, the HBoV3 strains showed an average amino acid divergence of 12% in the NS1 region and 20% in the VP1 region ( Table 2 ). All 11 strains were classified as members of the HBoV3 species, recently described by Arthur et al [8] . A distantly related variant of HBoV3 was also recently identified in a sewage sample in the United States (HBoV3B-CA-1-C1; Figure 1 ) [36] . To confirm whether this variant represents a second genotype of HBoV3 will require full genome sequencing. The HBoV4 cluster included 6 strains. The complete genome of 1 representative strain (NI-385) was obtained (GenBank accession number FJ973561). The NS1 protein of HBoV4 (NI-385) showed an average of 11% divergence with HBoV2 and 25%-27% divergence with HBoV1 or HBoV3 (Table 2) . For the VP1 protein, the relative divergences were reversed, with an average divergence of 8.5% with HBoV3, 9.5% with HBoV2, and 19% with HBoV1. NI-385, therefore, also appears to be a recombinant, with the 5 genes NS1 and NP1 closely related to HBoV2 (particularly genotype A) and the 3 VP1 gene slightly more similar to HBoV3 than HBoV2 (Table 2 ). According to the International Committee on Taxonomy of Viruses species demarcation criteria in the genus Bocavirus, members of different species must show 15% divergence in their nonstructural gene nucleotide sequences [37] . The genetic distance of NI-385 to its closest relatives in the NS1 gene (HBoV2A) was 10.8% (range, 6.8%-12.6%), indicating that NI-385 qualifies, pending International Committee on Taxonomy of Viruses review, as the prototype of a fourth HBoV species (HBoV4). A distant VP1 (partial) variant of HBoV4 was also detected in the United States (US-MN-964-05; Figure 1 ). To determine whether the latter variant represents a second HBoV4 genotype will require full genome sequencing. Nearly complete genomes and phylogenetic analysis of new bocavirus species. In a manner similar to HBoV1 and 2, all the new genomes of HBoV2, 3, and 4 encoded 3 large open reading frames (ORF). The left ORF encodes the nonstructural protein (NS), the middle ORF encodes NP1, and the right ORF encodes overlapping VP1/VP2 capsid proteins. Conserved motifs associated with rolling circle replication, helicase, and ATPase were identified within the NS protein. NP1 is a highly phosphorylated protein of currently undetermined function [38] ; NP1 differed in length between species, ranging from 214 to 219 amino acids. Situated within the VP1-unique (VP1u) region, the phospholipase A2 motifs required for parvovirus infectivity were found in all 6 genomes, together with the calcium-binding loop and catalytic residues. Further evidence of recombination in human bocaviruses. To further determine the relationship between members of the Bocavirus genus, phylogenetic analyses of NS1, NP1, and VP1/ VP2 were performed, with use of both nucleotide sequences and deduced protein sequences (Figure 2 ). The NS1 and NP1 genes of HBoV3 clustered with HBoV1, whereas their VP1/2 gene clustered with HBoV2. The incongruence in phylogenetic association between loci provided further evidence that HBoV3 originated from a recombination event bringing together the NS1/NP1 gene of HBoV1 and the VP1/2 gene of HBoV2 [8] . The likely recombinant origin of HBoV4, clustering with HBoV2 in the NS1/NP1 but with HBoV3 in the VP1, is also shown (Figure 2) . A scan of sequence divergence between complete genome sequences further supported the hypothesis of past recombi-nation between HBoV1 and 2 in the generation of HBoV3 and recombination of HBoV2 and HBoV3 in the generation of HBoV4, with both recombination points near the NP1 and VP1 junction (Figure 3 ). When different HBoV2 variants were similarly analyzed for recombination, intraspecies HBoV2 recombinants were also detected (data not shown) [7] . Diversity among respiratory HBoV1 and enteric HBoV2-4. We compared the intraspecies diversity of HBoV1 with that of HBoV2 by use of the partial VP1 sequence data generated with the pan-bocavirus PCR primers available for the greatest number of HBoV2, 3, 4 variants ( Figure 4A and 4B) . A very low average pair-wise difference was seen for HBoV1 collected worldwide ( Figure 4A ). HBoV2, including both genotypes, was more diversified than HBoV1, although HBoV2B alone showed a low level of diversity comparable to that of HBoV1 (HBoV2B generated the large low divergence peak in Figure 4A ). The homogeneity of HBoV1 and HBoV2B can also be visualized in the small branch lengths in Figure 1 . When the pair-wise dis- tances of all the enteric (non-HBoV1) sequences were plotted ( Figure 4C ), the distribution was much larger than that for HBoV1. The low level of intraspecies diversity of HBoV1 relative to the other species is also reflected in Table 2 . The genomic organization of HBoV species are remarkably similar to those of animal bocaviruses, except that all HBoV NS ORFs encode a shorter NS2 protein (630-650 amino acids), compared with the longer NS1 of animal bocaviruses (716-726 amino acids) ( Figure 5A ). We noticed in all 4 HBoV species the presence of a stretch of encoded amino acids similar to the C-terminus of the longer NS1 of animal bocaviruses overlapping the NP ORF, but in a different frame ( Figure 5A) . Genomes of all HBoV species were aligned to determine the presence of conserved potential RNA splicing signals near the end of the smaller NS2 ORF and the putative second exon encoding the C-terminal region of NS1 (Figure 5B). The putative NS1 resulting from such a spliced transcript encoded a 750-780 amino acid protein with a carboxy terminus that showed significant similarity to that of the canine and bovine bocaviruses NS1 ( Figure 5C ) [39, 40] . A study using RT-PCR for the detection of HBoV1 viral transcript in human lung epithelial cells did not detect the NS splicing proposed here [41] . On the other hand, the proposed NS1 RNA splicing and NS1 protein expression itself were detected using Northern blots and NS1 C-termini specific sera, respectively, in 293 and human epithelial cells transfected with plasmids expressing HBoV1 transcripts (Jianming Qiu, personal communication). Most PCR-positive stool samples contained HBoV2B (76 of 101), making this genotype the most commonly detected enteric human bocavirus (Table 1) . HBoV3 was identified in 11 stool samples and HBoV4 in 6, making them the second and third most common enteric human bocaviruses, respectively, in the regions analyzed here. HBoV1 was detected in only 4 of the 101 positive samples. Stool samples from patients with AFP in Nigeria and Tunisia both showed a high prevalence of HBoV2A+B (21%-26%). Because both patients with AFP and healthy contacts from Tunisia showed a comparable prevalence of HBoV2, no asso- ciation was observed between AFP and HBoV2 excretion. The prevalence of HBoV2 in adults with diarrhea from the United States and Nepal were also compared with those in healthy matched subjects. No associations were observed between HBoV2 shedding and diarrhea. Although the numbers of HBoV3 and HBoV4 detected were relatively small, it was noticed that the 8 HBoV3-and 4 HBoV4-positive samples were found among 192 samples from patients with AFP, whereas these 2 species were not detected among 96 healthy Tunisian matched contacts. Comparing the combined Nigerian and Tunisian patients with AFP infected with HBoV3 and HBoV4 to Tunisian controls yielded a 2-tailed Fisher's exact P value of .01, whereas comparing only the Tunisian patients with AFP with Tunisian controls gave a P value of .059. HBoV3 was also found in 1 patient with diarrhea from the United States and from 1 healthy person each from the United States and Nepal. HBoV4 was also found in a patient with diarrhea and 1 healthy adult from the United States. Further testing will be needed to confirm this trend of an association of HBoV3/4 with AFP in children or of the association of HBoV2 with diarrhea [8] . We report on a previously uncharacterized species of human bocavirus which we tentatively named HBoV4. HBoV1, HBoV2, and HBoV3 were also detected using a pan-PCR approach. HBoV2 has been detected in the stool of Pakistani children [7] and HBoV3 in stool samples of Australian children [8] , and both were recognized as recombinant viruses. The newly reported Australian HBoV3 (EU918736) is closely related to a Tunisian strain (TA-210-07) (Figure 2) . A highly prevalent genotype of HBoV2 (HBoV2B), together with partial genomic support for second genotypes of HBoV3 and HBoV4, were also identified. The availability of novel bocavirus genomes will allow the design of species-specific PCR or microarray oligonucleotides for their detection and for the disease association studies that are now required for the 3 recently characterized enteric human bocavirus species (HBoV2, HBoV3, and HBoV4). Based on the phylogenetic clustering observed for a large number of partial VP1 sequences ( Figure 1 ) and the distances among full genomes (Table 2) , we propose for future classification that HBoV strains showing 18% protein and 110% nucleotide difference in the complete VP1 gene should be considered different species, whereas those showing 11.5% protein and 15% nucleotide difference should be considered different genotypes. This VP1-based classification retains the 4 proposed human bocavirus species. The VP1 locus was selected because it is likely to strongly influence tissue tropism and potentially pathogenesis [42] . HBoV1 is primarily, although not exclusively [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] , a respiratory virus. We show here a higher prevalence of HBoV2 (particularly HBoV2B) than HBoV1 in stool samples from children and adults of different countries. A recent study testing for HBoV1 and HBoV2 DNA with use of species-specific nested PCRs failed to detect any HBoV2 in 16500 respiratory secretion samples from Edinburgh and Bangkok, whereas HBoV1 was found in 3% and 14% of these respiratory samples, respectively [43] . Another study found 5 (2%) of 212 nasopharyngeal samples from Korean children with acute lower respiratory track infections to be positive for HBoV2 DNA, but unexpectedly, no HBoV1 DNA was found [44] . Analyzing for HBoV1 and HBoV2 in both respiratory secretion and stool samples collected from the same individuals will be required to confirm whether the tropism of HBoV2 favors the digestive track and is distinct from that of the largely respiratory HBoV1. Extensive evidence for recombination was observed through full genome analyses, including the likely recombinant origin of HBoV3 and HBoV4 and the high level of intraspecies recombination between HBoV2 variants [7, 8] . The high prev-alence of bocavirus infection does provide the opportunity for coinfections, the first step in generating recombinant viruses. Indeed, a HBoV3 and HBoV4 coinfection was detected based on the pattern of mixed bases in 1 directly sequenced PCR product (data not shown). The frequent detection of HBoVs in stool samples from both healthy children and adults supports the likelihood of long periods of viral shedding and/or frequent reinfections. Determination of whether symptoms such as diarrhea are associated with acute infection in subsequently healthy viral shedders will require quantitative viral load measurements, analysis of longitudinally collected samples, and serological assays. Whether prior infection provides any protection against reinfection with the same or different genotype or species is also unknown. In both diarrhea sample sets analyzed here, consisting of samples from adults, no association with HBoV2 shedding was detected. If HBoV2 causes diarrhea, it may do so in only a small subset of infected children, possibly those without passively transferred maternal antibodies or without protective immune responses from prior infections. Coinfections with other enteric viruses may also exacerbate symptoms. Given the very large number of childhood infections (viral shedding prevalence of 120% in some developing countries), even low virulence would translate into a large disease burden. A trend of an association of HBoV3/HBoV4 detection with AFP was detected. The small numbers of cases will require independent confirmation. The neurological damages caused by bovine and canine bocaviruses in their animal hosts provide a precedent for infant nervous system pathogenicity [11] . The totality of HBoV1 sequences collected worldwide by multiple groups show very low protein and nucleotide sequence diversity (Table 2 and Figure 4A ) [35, 36] . In contrast, this single study found substantial diversity among HBoV2 and HBoV3, a fourth species (HBoV4), and extensive viral recombination (Figures 3, 4B, and 4C ). Assuming comparable rates of evolution, the genetically homogeneous and largely respiratory HBoV1, therefore, appears to be the more recently evolved species relative to the more diverse HBoV2, 3, and 4 found predominantly in feces. We propose that HBoV1 evolved from an enteric bocavirus ancestor that acquired, through mutation and/or recombination, enhanced respiratory tract tropism. Single-stranded DNA parvoviruses have been shown to have a mutation rate approaching that of RNA viruses, and recombination among animal parvoviruses has also been reported [45] [46] [47] . Parvoviruses have also demonstrated the capacity to rapidly expand their host species tropism, resulting in a recent pandemic in dogs [42, 48, 49] . A recent study showed that HBoV1 could replicate in differentiated human airway epithelial cells [41] . Whether HBoV2, 3, and 4 show an in vitro tropism and in vivo distribution that is more biased towards digestive track cells will require further study. Duration of immunity for canine and feline vaccines: a review Review of companion animal viral diseases and immunoprophylaxis Forty years of canine vaccination New DNA viruses identified in patients with acute viral infection syndrome Cloning of a human parvovirus by molecular screening of respiratory tract samples A newly identified bocavirus species in human stool A novel bocavirus associated with acute gastroenteritis in Australian children Human bocavirus Complete nucleotide sequence and genome organization of bovine parvovirus Animal bocaviruses: a brief review The association of newly identified respiratory viruses with lower respiratory tract infections in Korean children Human bocavirus infection Human Bocavirus in children Severe pneumonia and human bocavirus in adult Detection of human bocavirus in Japanese children with lower respiratory tract infections Evidence of human coronavirus HKU1 and human bocavirus in Australian children Human bocavirus: developing evidence for pathogenicity Human bocavirus infection in young children in the United States: molecular epidemiological profile and clinical characteristics of a newly emerging respiratory virus Epidemiological profile and clinical associations of human bocavirus and other human parvoviruses Frequent detection of human rhinoviruses, paramyxoviruses, coronaviruses, and bocavirus during acute respiratory tract infections Human bocavirus: clinical significance and implications Serodiagnosis of human bocavirus infection Human bocavirus: passenger or pathogen in acute respiratory tract infections? Emerging respiratory agents: new viruses for old diseases? Human bocavirus infection in children with gastroenteritis Human bocavirus in children hospitalized for acute gastroenteritis: a case-control study Clinical and molecular epidemiology of human bocavirus in respiratory and fecal samples from children in Hong Kong Detection of human bocavirus in children hospitalized because of acute gastroenteritis Human bocavirus, a respiratory and enteric virus Human bocavirus infection in children hospitalized with acute gastroenteritis in China Complete coding sequences and phylogenetic analysis of human bocavirus (HBoV) Genotyping of human bocavirus using a restriction length polymorphism MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences Structural constraints on RNA virus evolution Frequent detection of highly diverse variants of cardiovirus, cosavirus, bocavirus, and circovirus in sewage samples collected in the United States Virus taxonomy. Eighth report of the International Committee on Taxonomy of Viruses Virally coded noncapsid protein associated with bovine parvovirus infection Molecular characterization of infectious clones of the minute virus of canines reveals unique features of bocaviruses The transcription profile of the bocavirus bovine parvovirus is unlike those of previously characterized parvoviruses Human bocavirus can be cultured in differentiated human airway epithelial cells The origins of new pandemic viruses: the acquisition of new host ranges by canine parvovirus and influenza A viruses Absence of detectable replication of the human bocavirus species 2 in the respiratory tract despite frequent detection in faecal samples Human bocavirus 2 in children Rates of evolutionary change in viruses: patterns and determinants Rapid sequence change and geographical spread of human parvovirus B19; comparison of B19 evolution in acute and persistent infections Comparative analysis reveals frequent recombination in the parvoviruses Phylogenetic analysis reveals the emergence, evolution and dispersal of carnivore parvoviruses High rate of viral evolution associated with the emergence of carnivore parvovirus We thank Dr Michael P. Busch and Blood Systems Research Institute for sustained support and Hope Biswas for statistical assistance.