key: cord-0902589-328b4t5c
authors: Aryal, Manita; Liu, Guangliang
title: Porcine Bocavirus: A 10-Year History since Its Discovery
date: 2021-04-28
journal: Virol Sin
DOI: 10.1007/s12250-021-00365-z
sha: ed169161aec7ad41a56a0674141118508285050a
doc_id: 902589
cord_uid: 328b4t5c

Porcine bocavirus (PBoV) is a single-stranded DNA virus, belongs to the genus Bocaparvovirus of family Parvoviridae. It was discovered along with porcine circovirus 2 (PCV 2) and torque tenovirus (TTV) in the lymph nodes of pigs suffering from postweaning multisystemic wasting syndrome (PMWS) in Sweden in 2009. PBoV has been reported throughout the world, mostly in weaning piglets, and has a broad range of tissue tropism. Since PBoV is prevalent in healthy as well as clinically infected pigs and is mostly associated with coinfection with other viruses, the pathogenic nature of PBoV is still unclear. Currently, there are no cell lines available for the study of PBoV, and animal model experiments have not been described. This review summarizes the current state of knowledge about PBoV, including the epidemiology, evolution analysis, detection methods, pathogenesis and public health concerns.

Bocaviruses (Bo: bovine; ca: canine) were first reported in the 1960s (Blomström et al. 2009; Zhou et al. 2014; Lanave et al. 2015) with possible recombination between bovine parvovirus and canine minute virus (Allander et al. 2005; Lau et al. 2011) . Bocavirus is an emerging pathogen and has been reported in canines (Bodewes et al. 2014) , bovines (Mitra et al. 2016) , humans (Allander et al. 2005) , pigs (Blomström et al. 2009 ), gorillas (Nze-Nkogue et al. 2017) , chimpanzees (Brožová et al. 2016) , California sea lion , rabbits (Lanave et al. 2015) , Himalayan marmots (Ao et al. 2017) , camels (Woo et al. 2017) , wild boars (Cadar et al. 2011) , bats (Lau et al. 2017) , minks (Wang et al. 2017) , house shrews (Xiong et al. 2018 ) and rodents (Lau et al. 2016; Xiong et al. 2018; Zhang C et al. 2018) .

Porcine bocavirus (PBoV), the new member of Bocaparvovirus genus, was discovered in Sweden in 2009 in pigs from finishing herds suffering from postweaning multisystemic wasting syndrome (PMWS) (Blomström et al. 2009 ). Within a year, it was identified in healthy pigs (Shan et al. 2011a) , diarrheic pigs (Shan et al. 2011b ) and the pigs suffering from respiratory problems (Cheung et al. 2010; Zhai et al. 2010) . It has been reported in Asia (Zhai et al. 2010) , Europe (Cadar et al. 2011) , Africa (Blomstrom et al. 2013; Ndze et al. 2013) and North America (Schirtzinger et al. 2015) . Although this virus has been found for 10 years, there still many questions remain unresolved. This review describes the genome structure, classification, detection methods, molecular epidemiology and public health importance of this virus.

TTV-2 and PCV-2 but had no detectable level of a boca-like virus. However, the pig that suffered from PMWS had TTV-1, TTV-2, PCV-2, and boca-like parvovirus (Blomström et al. 2009 ). The prevalence rate of porcine boca-like parvovirus was 46% in pigs without PMWS and 88% in pigs with PMWS. The prevalence was almost twice as high in clinically affected pigs as in unaffected pigs (Blomstrom et al. 2010) . The complete genome of PBoV was first determined in China in 2010 (PBoV1/2CHN) from fecal samples (Cheng et al. 2010) .

Later, PBoV was identified in archived samples (collected in 2006) from Croatia (gene accession number: KF206161.1, KF206162.1) and China (Cheng et al. 2010) . Similarly, bocaviruses were also identified in samples from the wild boar population in Romania collected in 2006 Romania collected in /2007 Romania collected in and 2010 Romania collected in /2011 . The prevalence rate of bocaviruses was found to increase with time [from 9.14% (43/470) in 2006/2007 to 17.74% (66/ 372) in 2010/2011] (Cadar et al. 2011) . These studies show that bocaviruses circulated in the swine population independently before they were discovered in 2009 (Zhou et al. 2014) .

Since 2009, the virus has been reported throughout the world (Fig. 1) , including Sweden (Blomström et al. 2009 ), China (Cheng et al. 2010; Zhai et al. 2010) , the USA (Cheung et al. 2010; Shan et al. 2011b; Jiang et al. 2014 ), Ireland (McKillen et al. 2011 Gunn et al. 2015) , Romania (Cadar et al. 2011) , Hungary (Csagola et al. 2012) , Croatia Keros et al. 2017) , Cameroon (Ndze et al. 2013) , the UK (McMenamy et al. 2013), Thailand (Saekhow and Ikeda 2014) , Korea (Choi et al. 2014) , the Czech Republic, Slovakia (Vlasakova et al. 2014) , Mexico (Schirtzinger et al. 2015) , Germany (Pfankuche et al. 2016) , Japan (Zhang et al. 2016) , Uganda, Kenya (Amimo et al. 2017) , Slovenia (Keros et al. 2017 ), Belgium (Conceicao-Neto et al. 2017 , and Malaysia (Jacob et al. 2018) . Even though PBoV has been identified in many countries (Fig. 1) , its prevalence may differ based on geographical location, pig age, and pig herd management.

In China, the prevalence rate of PBoV in stool samples of piglets was found to be 12.59% (Cheng et al. 2010) . Similarly, the prevalence rate of bocaviruses in clinical samples from one-month-old piglets was reported to be 5.77% (9/156) . Another study reported that the prevalence rate of PBoV was 18.1% (84/463) (Zhou et al. 2018) . Serological analysis of 2025 serum specimens, revealed that 856 (42.3%) of serum specimens were positive for anti-PBoV IgG antibodies (Shi et al. 2019) . Comparative analysis of the prevalence rate of PBoV in the USA and China showed that the prevalence rate of PBoV in the USA (42%) was higher than that in China (11.4%) (Zhang et al. 2015) .

In Croatia, the prevalence rate of PBoV was found to be 100% (10/10) in small commercial farms, 45.45% (5/11) in composite fecal samples after small backyard holdings, 42.1% (24/57) in small commercial pig farms and 27.9% (12/43) in small backyard holdings (Keros et al. 2017) . In Malaysia, out of 103 samples from 17 pigs, 32 samples from 15 pigs were positive for PBoV, revealing a prevalence rate of 90.9% (Jacob et al. 2018) . In Thailand, eighty tonsil samples from pigs were collected and analyzed from Virologica Sinica a slaughterhouse; the prevalence rate was 18% (23/80) (Saekhow and Ikeda 2014) .

The occurrence of bocavirus infection was found to be higher in the spring (72.7%, 16/22) than in the summer (28.9%, 26/90), autumn (38.7%, 12/31) or winter (41.7%, 20/48). Nevertheless, there is no significant difference in PBoV infection rate between seasons (Zhai et al. 2010) .

In a study from China, the prevalence of PBoV was found to be higher in grower pigs than in pigs of other age groups . In another study from Malaysia, the prevalence of bocavirus was higher in weaning and growing pigs than in newborn piglets (Jacob et al. 2018) . Similarly, in samples from Romania, the prevalence of bocavirus infection was found to be higher in 6-to 12-month old animals (77.06%, 84/109) than in 12-to 36-month old animals (22.94%, 25/109). This result indicates that infection is more common in younger pigs than in older pigs (Cadar et al. 2011) . In Korea, the Czech Republic and Slovakia, it has been reported that infection was more common in weaned piglets than in pigs of other age groups (Choi et al. 2014; Vlasakova et al. 2014) . In piglets, maternal antibodies might play an important role in the prevention of bocavirus infection (Vlasakova et al. 2014 ). However, this hypothesis requires further investigation for confirmation.

PBoV, which is a nonenveloped virus, exhibits icosahedral symmetry (Cheng et al. 2010; Zhou et al. 2017) , has a diameter of 26-30 nm (Yang et al. 2012 ) and consists of linear, single-stranded DNA of 4-6 kb (Cheng et al. 2010; Zhou et al. 2017) of either positive or, most often, negative polarity, with most strains being negative-sense (Anne Christin Böhmer et al. 2009 ). As shown in Fig. 2 , the genome of PBoV consists of three open reading frames (ORFs) (Zhou et al. 2018) . ORF-1 encodes the nonstructural protein NS1, and ORF-2 encodes the structural capsid proteins VP1 and VP2 (Zhou et al. 2018) . VP1 and VP2 overlap each other (Zhou et al. 2014) . ORF-3 is regarded as the characteristic feature of bocavirus (Cheng et al. 2010; Shan et al. 2011b ) and encodes the nonstructural protein NP1 (Zhou et al. 2018) .

The ORF-1 region of PBoV resembles bovine parvovirus (BPV), and ORF-2 resembles canine minute virus (CMV) (Lau et al. 2011) . ORF-1 encodes a protein associated with rolling circle replication, helicase, and ATPase activities. ORF-2 encodes a protein associated with phospholipase A2 motifs (the most conserved motif is YXGXG) required for parvovirus infectivity, calcium-binding loop, and catalytic residues (the most conserved motif is HDXXY) (Cheng et al. 2010; Zhou et al. 2017) . The conserved motif YXGXF carried by porcine bocavirus is different from the YXGXG motif of other parvoviruses (Cheng et al. 2010; Zhou et al. 2017) . However, in addition to the HDXXY motif and YXGXF motif other motifs have been reported in PBoV, including the HDQAY motif in strain GD6, the HDLAY motif in strain GD10, and the YLGPF motif in strains GD6, GD10, and GD23 (Zhou et al. 2017) . The function of the protein encoded by ORF-3 is not clear, but in CMV, it has been associated with viral replication (Lau et al. 2011) .

Episomes have been detected in PBoV, and replication takes place by a rolling circle mechanism with a head-totail structure. Hairpin-1 of PBoV has a different position and is conserved (Zhou et al. 2014 ).

According to the International Committee on Taxonomy of Viruses, PBoV was classed to the genera Bocaparvovirus (BoV), subfamily Parvovirinae and family Parvoviridae (Cotmore et al. 2013 (Cotmore et al. , 2019 Zhou et al. 2018) . Previously, PBoV strains were designated as novel strains if they had \ 95% homology of the nonstructural gene (NS1) with existing strain. Because of this high level of required homology, many current species consist of single isolates, creating confusion between taxa and viruses in the literature and hindering the taxonomic structure. According to a new proposal, bocavirus strains are designated as novel strains if they have less than 85% amino acid sequence identity of the NS1 protein. Based on this new criterion, most of the previously identified viruses belong to the same taxa, allowing individual species to contain more diverse viruses. Other existing criteria, such as host, antigenic properties, and genome characteristics, are still used. The PBoV has been placed under the Ungulate bocaparvovirus five type species with seven PBoV variants (Cotmore et al. 2013) .

Originally, PBoV were classified into three (PBoV-likeV, PBoV1/PBoV2, and 6 V/7 V) phylogenetic clades based on nucleotide and amino acid alignments (Shan et al. 2011a ). Further, some other researchers proposed classification methods based on ORF-2 encoding complete VP1 genes using the criteria of [ 40% nucleotide differences considering members of different groups and [ 10% difference from the VP1 nucleotide sequence considering members of different subgroups (Yang et al. 2012) . In another study, PBoV from China and the USA were divided into six groups based on 400 bp VP1 fragment: PBoV1 (group 1), PBoV2 (group 2), PBoV3 (group 3), PBoV4 (group 4), PBoV5 (group 5), and PBoV6V7V (group 6). PBoV1 and PBoV2 were epidemic strains from 2006 to 2011 in China followed by PBoV3 subtypes epidemic from 2010 to 2012 in China and the USA. Further reported that PBoV3C, PBoV5, PBoV3/4 were epidemic viruses in China and the USA (Zhang et al. 2015) . Based on phylogeographical analysis, it was reported that Hong Kong, Hebei Province, and Jiangsu Province were likely the geographical origins of PBoV in China, and that Minnesota and North Carolina were likely the geographical origins of PBoV in the USA (Zhang et al. 2015) . In the recent study, PBoV were classified into three groups based on the complete VP1 gene; group 1 (HQ223038, HQ291308), group 2 (KX017193, HM053694), and group 3 (JF12472, KF025381). Group 3 PBoV were further subdivided into five subgroups (Zhou et al. 2017) .

The relationships among PBoV strains isolated from different countries can be observed from the phylogenetic tree based on the full length sequence (Fig. 3) . The result demonstrate that there are three major groups. Group 2 and 3 forms distinct monophyletic groups in comparison to group 1 which forms a paraphyletic group with many monophyletic groups. Group 2 contains PBoV isolated from China, Kenya and Belgium. The nucleotide sequence of these PBoV were truncated in the NS1 gene. Group 3 contains PBoV isolated from China and European countries. Group 1 contains PBoV isolated from Asia, Europe and the USA. It can be observed that PBoV isolated throughout the world are diverse and tends to cluster with other strains isolated from different countries (Fig. 3) .

Recombination is considered an important event for emerging viruses and has been described in human and porcine bocaviruses (Lau et al. 2011; Zhou et al. 2014 Zhou et al. , 2017 . Recombination has been described in the NS1 gene, NP1 gene or VP1/VP2 gene. A recombination event between two recombination breakpoints near sites 1170 and 1780 of the VP1 sequence was demonstrated in a fecal sample (SH14F/12) containing PBoV4 (Lau et al. 2011) . Similarly, a natural recombination breakpoint was identified in the NS1 gene (nucleotides 925 to 1959), but no genetic recombination was detected (Zhou et al. 2017) . Also no recombination signal was reported for PBoV1 or PBoV2 (Cheng et al. 2010) . Parvoviruses are known to have a high mutation rate of the capsid gene, with 10 -4 nucleotide substitutions per site per year (Karuppannan and Opriessnig 2018) . Recombination events as well as a high mutation rate may lead to better adaptation of bocaviruses as well as host escape.

Bocavirus isolated from other host species show near relationship with porcine bocavirus. The bocavirus isolated from minks shared 87% nucleotide similarity with porcine bocavirus (strain: PBoV-KU14; HQ223038, HQ291308, KJ622366) and formed separate clade with porcine bocavirus (Wang et al. 2017) . Similarly, five full length genomes of PBoV (HD.44, HD.52, MM.57, YN.27, and YN.22) were constructed based on sequences generated from murine rodents and formed separate clade with group four porcine bocavirus (Xiong et al. 2018) . Additionally, bocavirus isolated from bats, rodents and Himalayan marmonts are closely related with porcine bocavirus. The bocavirus isolated from Himalayan marmonts were closely related to porcine bocavirus (PBoV-SX/1-2) (Ao et al. 2017) . The bocavirus isolated from bats were closely related to porcine bocavirus (NC_024453.2, NC_016032.1, HQ223038.1, HQ291308.1) (Lau et al. 2017) . Similarly, the bocavirus isolated from rodents species were related with porcine bocavirus isolate PBoV-KU14 with 92.1%-92.9% amino acid identities in NS1 protein . Therefore these studies clearly reflects that interspecies transmission must have occurred from minks, bats, Himalayan marmonts and rodents to the porcine. However no any experimental evidence of interspecies transmission of bocavirus have been reported till date. Such studies can reflect the wide range of host adaptation of bocavirus as well as genetic variation in the porcine bocavirus.

The random amplification and high-throughput sequencing approach was implemented in Sweden (Blomström et al. 2009 ) to discover the porcine boca-like virus. These methods are gaining popularity for the discovery of novel viruses and to address metagenomic issues (Blomström et al. 2009; Li et al. 2011; Shan et al. 2011b; Yang et al. 2012; Xiong et al. 2018 ).

PCR-based detection methods have been developed to rapidly detect PBoV. Conventional multiplex PCR can be used to detect the three groups of PBoV. The sensitivity of single PCR assays is equal to that of the multiplex PCR assay and is 1.0 9 10 3 genomic copies/lL for PBoV group 1, 4.5 9 10 3 copies/lL for PBoV group 2 and 3.8 9 10 3 copies/lL for PBoV group 3 (Zheng et al. 2016a ). Loopmediated isothermal amplification (LAMP) was developed for rapid, specific, and sensitive detection of porcine bocalike virus. It has been reported that this method is 100 times more sensitive than conventional PCR and has a detection limit of approximately ten copies per reaction. Moreover, no crossreactivity has been reported with common circulating swine viruses such as porcine reproductive and respiratory syndrome virus (PRRSV), porcine circovirus type 2 (PCV-2), porcine parvovirus (PPV), and classical swine fever virus (CSFV) (Li et al. 2013) . TaqMan-based real-time PCR (qPCR) has been used for the detection and quantitation of ungulate bocavirus 2 by targeting the VP1 gene. This method is highly sensitive and specific and can detect the real-time quantification of viral load without post-PCR detection procedures (Zhou et al. 2018 ). An EvaGreen-based multiplex real-time PCR method was developed for simultaneous detection and grouping of the three bocavirus groups (groups 1, 2, and 3) in a single step (Zheng et al. 2016b) . Duplex nano PCR has been used for the simultaneous amplification and Fig. 3 Phylogenetic analysis of bocaviruses. The phylogenetic tree was generated using full length nucleotide sequences from bocaviruses in GenBank using the MEGA 7 software package (neighbor joining method with 1000 bootstrap replicates). The relationship among the porcine bocaviruses from different countries can be observed from the tree. Detailed information related to certain PBoV strains is listed in Table 1. quantification of PBoV and pseudorabies virus (Luo et al. 2015) . SYBR Green-based duplex real-time PCR has been used for detection of porcine epidemic diarrhea virus (PEDV) as well as PBoV genotypes (3, 4, 5), with detection limits of 10 copies/lL for both viruses (Zheng et al. 2020) . Sequence-independent single primer amplification (SISPA) is a primer-initiated method that uses a single primer to amplify nucleic acids of an unknown sequence by sequenceindependent PCR. This method has been used to amplify PBoV (Cheng et al. 2010) . Nested polymerase chain reaction has been used for the identification of PBoV from throat swab, fecal, and serum samples by using two nested PCR primers (Xiong et al. 2018; Zhou et al. 2018 ). In the primer walking method, each specific primer is used along with a degenerate primer to hybridize to all known bocavirus sequences (Zeng et al. 2011) . For this method, other sequences should be available to design specific primers. Nucleic acid sequence-based amplification (NASBA) has been used to determine the polarity of the human bocavirus genome and to analyze the transcription of uncharacterized genes (Anne Christin Böhmer et al. 2009 ).

Immunological methods, such as indirect immunofluorescence assays, have been developed to detect unidentified isolated PBoV in primary porcine kidney cell lines. In 

Many researchers have reported PBoV as a potential emerging pathogen (Blomström et al. 2009; Zhou et al. 2018 ) of the respiratory tract (Zhai et al. 2010) or the gastrointestinal tract (Amimo et al. 2017 ). Diversity of host species is assumed to be the specific feature of emerging viruses . To date, PBoV has not been reported as an enteric pathogen in the literature. However, the co-occurrence of this virus with circulating viruses has been reported since its discovery (Table 1) Zhou et al. 2018) , porcine parvovirus , porcine pseudorabies virus (Luo et al. 2015) and porcine kobuvirus . Clinically, the incidence of coinfection of PBoV with PCV-2 has been reported to be 83.8%, suggesting that PCV-2 enhances the infectivity of PBoV ). Thus, PBoV may not be directly associated with disease and might function as a helper virus for triggering other infectious agents .

Analysis of 12 fecal samples from diarrheic piglets and 24 fecal samples from healthy piglets by Shan et al. revealed that bocavirus (5/12 vs. 3/24) and coronavirus (11/12 vs. 13/24) are more prevalent in diarrheic than in healthy piglets (Shan et al. 2011b ). The researchers concluded that the high rate of coinfection seen in high-density farms provides favorable conditions for recombination and accelerates viral evolution. In Uganda and Kenya, bocavirus was isolated from piglets with no signs of clinical diarrhea (Amimo et al. 2017) . During a diarrhea outbreak in fattening pigs from Belgium, PEDV was detected along with PBoV, recombinant enterovirus-toroviruses, astroviruses, enteroviruses, and picobirnaviruses. The presence of PEVD (Cq [ 30) in low loads along with high loads of PBoV and several enteroviruses suggests that the cause of viral diarrhea was a virus other than PEVD and that bocavirus was the most likely cause (Conceicao-Neto et al. 2017 ).

The pathogenesis of PBoV has not been determined to date. It is thought that, as PBoV is highly prevalent in the swine population and has high genetic diversity, its pathogenesis could be determined by direct evidence from clinical disease (Zhou et al. 2017) . The detection of bocavirus in different tissues represents a wide range of tissue tropism (Lau et al. 2011) . PBoV was detected for the first time in the lymph nodes of weaning pigs suffering from PMWS (Blomström et al. 2009 ). Soon after that, it was detected in fecal samples from piglets (Cheng et al. 2010 ) as well as the respiratory tract (Zhai et al. 2010) . In diarrheal disease, histopathological changes, including microscopic lesions and villous atrophy, are mainly detected in the jejunum, ileum, and duodenum . Bocavirus has been isolated from piglets with encephalomyelitis, demonstrating the pathological role of this virus (Pfankuche et al. 2016) . In a study novel porcine parvovirus, PPV4, similar to porcine bocavirus was identified from the lung lavage of a diseased pig coinfected with PCV-2. Tissue homogenates consisting of lung, lymph node, spleen, and heart were prepared and inoculated into 

two colostrum-deprived pig oronasally and monitored clinically. Both pigs developed respiratory disease and were euthanized. However, due to the coinfecting virus PCV-2, it was not clear if it caused disease on its own or contributed to the clinical disease (Cheung et al. 2010) . Information related to the isolation of PBoV from different tissues is presented in Table 2 . Lymph nodes, spleen, and tonsil have the highest detection rate suggesting that these organs are the sites of active replication and the pathogenesis of PBoV infections may involve the lymphoid tissues and gastrointestinal tract (Jacob et al. 2018) . Besides tissues, bocavirus has been detected from the saliva (Choi et al. 2014 ) and serum (Csagola et al. 2012; Choi et al. 2014; Meng et al. 2018) . In brief, bocavirus has been isolated from different tissues of the pig, but it is still not clear which cells in pigs are responsible for its multiplication and dissemination throughout the body. However, in a study conducted in Thailand in dogs, it was reported that most of the viruses belonging to the Parvoviridae family replicate in mitotically active cells such as intestinal crypt epithelial cells. Analysis of tissue samples from the small intestine revealed nuclear signals for canine bocavirus 2 in enterocytes, mainly those located at the villus tips and crypts. Transmission electron microscopy showed numerous electron-dense icosahedral viral particles measuring approximately 20 nm in diameter that aggregated to form large intranuclear inclusion bodies within apical small intestinal enterocytes (Piewbang et al. 

Key findings References

Lymph nodes The infection rates of PBoV in the PMWS-affected pigs were twice higher than in the non-PMWS affected pigs. The co-infection of PBoV along with the TTSV and PCV-2 might have facilitate the development of PMWS Blomström et al. (2009) The co-existence of two bocavirus strain within the same fecal sample revealing inter and intra host genetic diversity

Lau et al.

The genetic diversity of the circulating bocavirus strains in Xinjiang belong to three subgroups of three different genetic groups Meng et al. (2018) Gastrointestinal tract

The prevalence rate of PBoV was higher in stool samples and these viruses multiply in the intestinal tract of piglets Cheng et al. (2010) Respiratory tract The first evidence of infection of weaning piglets with respiratory tract symptoms representing an emerging virus for swine respiratory tract diseases

Zhai et al.

Nasopharyngeal sample

PBoV were higher in nasopharyngeal samples in deceased pigs than in healthy pigs Lau et al.

Mesenteric lymph nodes

The mesenteric lymph node had the highest detection rate suggesting the pathogenesis of PBoV infection involves the lymphoid tissues Jacob et al.

Inguinal lymph nodes 25% of the organ tested were positive for PBoV Jacob et al.

Spleen 23.5% of the spleen tested were positive for PBoV Jacob et al. . This study provides novel insights into the pathogenicity of canine bocavirus infections, which could also apply to other bocaviruses isolated from different animals. Rats naturally infected with parvovirus are mostly healthy. However, in some cases, growth retardation and fetal loss is observed in pregnant rats experimentally inoculated with parvovirus (Lau et al. 2016) . Cells lines for the isolation of bocavirus have not yet been reported. Different cells such as porcine kidney cells (PK-15), swine testicular cells, porcine alveolar macrophages, monkey kidney cells , and human embryonic kidney epithelial cells (HEK293T) have been used for the propagation of PBoV. However, no success has been reported (Zeng et al. 2011) . During the diarrheal disease outbreak in Belgium, an attempt was made to isolate bocavirus and recombinant enterovirus-torovirus from primary porcine kidney epithelial cells and swine testicular cells. Although a cytopathic effect was observed in porcine kidney epithelial cells after two days, no virus could be detected in passaged cells (Conceicao-Neto et al. 2017) . In another study in Northern Ireland, fecal suspensions and homogenates of the small intestines of 6-weekold piglets exhibited cytopathic effect after four passages in cultured primary porcine kidney cells. PCV1, PCV-2, PPV, porcine enterovirus types 1, 2, and 3, porcine adenovirus, and porcine reovirus were not detected by PCR or RT-PCR; thus, this study was the first to describe the growth of bocavirus in a primary cell line . Because this virus cannot be cultured in cell lines and since it commonly coexists with other circulating swine enteric viruses, the exact role of this virus in disease progression is unclear.

Due to the association of bocavirus with other circulating pathogenic enteric viruses, its nature has always been confusing. Due to the involvement of bocavirus in respiratory tract infection as well as gastrointestinal tract infection, it is still unclear whether bocaviruses are respiratory pathogens or enteric pathogens, especially in humans, pigs and other animals. Due to the lack of an experimental animal model, the pathogenic nature of PBoV has not been investigated.

As for other enteric pathogens, the fecal-oral route may be the main transmission route of PBoV. Diarrheal feces and/or vomitus and other contaminated fomites, such as transport vehicles and feeds, may be significant transmission sources. Since PBoV is a nonenveloped virus, it can survive harsh environment and is also responsible for the widespread distribution of various parvoviruses (Saekhow and Ikeda 2014) . In 2007, human bocaviruses were isolated from 17 of 21 (81%) raw sewage samples from different states in the USA (Blinkova et al. 2009 ). This finding reflects the robust nature of bocaviruses. Once they enter swine herds, they can cause subclinical infection of pigs and circulate in the swine population. However, whether aerosolized PBoV is infectious should be investigated.

Human bocavirus, human parvovirus 4 and parvovirus B19 of the family Parvoviridae are known to be associated with human disease (Martin et al. 2009; Karuppannan and Opriessnig 2018) . Human bocaviruses have been reported worldwide (Rikhotso et al. 2018 ). The first case of human bocavirus was reported in the respiratory tract of a child from Sweden in 2005, four years before the discovery of bocavirus in swines (Allander et al. 2005) . Bocaviruses have been isolated from the respiratory tract (Vicente et al. 2007; Arthur et al. 2009; Martin et al. 2009 ) and gastrointestinal tract (Zhao et al. 2012) , and from cancer patients (Schildgen et al. 2013) . They have also been isolated from the respiratory and gastrointestinal tracts of healthy and clinically infected children. Therefore, in humans, it is still unclear whether bocaviruses are real pathogens, opportunistic pathogens or indwelling flora.

In a case study from Iran published in 2018, PBoV was detected in a three-year-old child with an acute respiratory tract infection (Safamanesh et al. 2018) , implying that PBoV may be evolving as a human pathogen.

The xenotransplantation of pig organs is gaining greater attention (Karuppannan and Opriessnig 2018) . With the increasing concern of the pathogenic role of PBoV, screening, monitoring and risk analysis of ssDNA viruses is recommended to eliminate the risk of transmission during xenotransplantation of pig organs.

In southern China, brown rats from markets or restaurants were infected with the rat bocavirus. Further investigations are required to study the potential for interspecies transmission from rats to humans (Lau et al. 2016) .

Zoonotic diseases have caused devastating effects in human beings. There is increasing concern regarding the broad host range, wide tissue tropism and increasing pathogenic potential of bocaviruses. These studies indicate that bocaviruses have important public health implications.

Bocavirus is an emerging pathogen with a broad host range and wide tissue tropism. In the past decade, many researchers have tried to isolate this virus in cell culture, but success has not been achieved. Thus, an organ coculture system might be beneficial for its propagation. Studies in such a system could be used to determine whether this virus favors the respiratory tract, the gastrointestinal tract or other tissue systems for its propagation.

Many studies from various countries have reported that PBoV is associated with other coinfecting pathogens of swine, including PEDV, PRRSV, CSFV, and rotavirus A. While PBoV was discovered ten years ago, we still do not know if it is a real pathogen or an opportunistic pathogen. Therefore, single virus isolation could be used for the study of animal infection models. Additionally, direct evidence of the pathogenesis of this virus is needed.

Phylogenetic analysis of PBoV isolated from minks and mice has revealed complex evolutionary mechanisms, interspecies transmission and recombination points. Further studies are needed to elucidate interspecies transmission between different host species.

Bocavirus poses a significant public health concern. The clinical case from Iran (Safamanesh et al. 2018 ) is a reminder that it is necessary to understand emerging pathogens such as PBoV because they might have devastating effects in the future. As mentioned above, our significant contribution to the understanding of PBoV will facilitate the prevention and control of this emerging virus and reduce its economic impact.

Cloning of a human parvovirus by molecular screening of respiratory tract samples

First complete genome sequences of porcine bocavirus strains from East Africa

Novel application for isothermal nucleic acid sequence-based amplification (NASBA)

Two novel bocaparvovirus species identified in wild Himalayan marmots

A novel bocavirus associated with acute gastroenteritis in Australian children

Frequent detection of highly diverse variants of cardiovirus, cosavirus, bocavirus, and circovirus in sewage samples collected in the United States

Detection of a novel porcine boca-like virus in the background of porcine circovirus type 2 induced postweaning multisystemic wasting syndrome

Studies of porcine circovirus type 2, porcine boca-like virus and torque teno virus indicate the presence of multiple viral infections in postweaning multisystemic wasting syndrome pigs

Genetic characterisation of a porcine bocavirus detected in domestic pigs in Uganda

Novel canine bocavirus strain associated with severe enteritis in a dog litter

Genetic and phylogenetic characterization of novel bocaparvovirus infecting chimpanzee

Genetic detection and analysis of porcine bocavirus type 1 (PoBoV1) in European wild boar (Sus scrofa)

Emerging novel porcine parvoviruses in Europe: origin, evolution, phylodynamics and phylogeography

Identification and nearly full-length genome characterization of novel porcine bocaviruses

Identification and molecular cloning of a novel porcine parvovirus

Molecular detection and genetic analysis of porcine bocavirus in Korean domestic swine herds

Identification of an enterovirus recombinant with a torovirus-like gene insertion during a diarrhea outbreak in fattening pigs

The family parvoviridae

ICTV virus taxonomy profile: parvoviridae

Detection, prevalence and analysis of emerging porcine parvovirus infections

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

First molecular detection of porcine bocavirus in Malaysia

High prevalence and genetic diversity of porcine bocaviruses in pigs in the USA, and identification of multiple novel porcine bocaviruses

Possible risks posed by singlestranded DNA viruses of pigs associated with xenotransplantation

The silent spread of porcine bocavirus in croatian pigs: Should we be concerned?

Novel bocaparvoviruses in rabbits

Co-existence of multiple strains of two novel porcine bocaviruses in the same pig, a previously undescribed phenomenon in members of the family Parvoviridae, and evidence for inter-and intra-host genetic diversity and recombination

Identification and genomic characterization of a novel rat bocavirus from brown rats in China

Bats host diverse parvoviruses as possible origin of mammalian dependoparvoviruses and source for bat-swine interspecies transmission

The fecal viral flora of California sea lions

Development of a loopmediated isothermal amplification method for rapid detection of porcine boca-like virus

Concurrent infections of pseudorabies virus and porcine bocavirus in China detected by duplex nanoPCR

Detection of bocavirus in saliva of children with and without respiratory illness

Isolation in cell cultures and initial characterisation of two novel bocavirus species from swine in Northern Ireland

Detection of a porcine boca-like virus in combination with porcine circovirus type 2 genotypes and Torque teno sus virus in pigs from postweaning multisystemic wasting syndrome (PMWS)-affected and non-PMWS-affected farms in archival samples from Great Britain

Production, characterisation and applications of monoclonal antibodies to two novel porcine bocaviruses from swine in Northern Ireland

Molecular detection and genetic diversity of porcine bocavirus in piglets in China

Metagenomic characterization of the virome associated with bovine respiratory disease in feedlot cattle identified novel viruses and suggests an etiologic role for influenza D virus

Detection of novel porcine bocaviruses in fecal samples of asymptomatic pigs in Cameroon

Identification and molecular characterization of novel primate bocaparvoviruses from wild western lowland gorillas of Moukalaba-Doudou National Park

Porcine bocavirus infection associated with encephalomyelitis in a pig

Canine bocavirus type 2 infection associated with intestinal lesions

Prevalence of human bocavirus in africa and other developing countries between 2005 and 2016: a potential emerging viral pathogen for diarrhea

Prevalence and genomic characterization of porcine parvoviruses detected in Chiangmai area of Thailand in 2011

Detection of porcine bocavirus from a child with acute respiratory tract infection

The human bocavirus is associated with some lung and colorectal cancers and persists in solid tumors

First identification of porcine parvovirus 6 in North America by viral metagenomic sequencing of serum from pigs infected with porcine reproductive and respiratory syndrome virus

Genomic characterization and high prevalence of bocaviruses in swine

The fecal virome of pigs on a high-density farm

Seroprevalence of porcine bocavirus in pigs in northcentral China using a recombinant-NP1-protein-based indirect ELISA

Human bocavirus, a respiratory and enteric virus

The presence of six potentially pathogenic viruses in pigs suffering from post-weaning multisystemic wasting syndrome

Genome sequence of a porcine bocavirus detected in feces of domestic minks in China

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

Detection and phylogenetic analysis of porcine bocaviruses carried by murine rodents and house shrews in China

Genome characterization of a novel porcine bocavirus

Complete coding sequences and phylogenetic analysis of porcine bocavirus

High prevalence of a novel porcine bocavirus in weanling piglets with respiratory tract symptoms in China

Porcine bocaviruses: genetic analysis and prevalence in Chinese swine population

Occurrence and investigation of enteric viral infections in pigs with diarrhea in China

Evolutionary, epidemiological, demographical, and geographical dissection of porcine bocavirus in China and America

Virus-like particles of porcine bocavirus generated by recombinant baculoviruses can be applied to sero-epidemic studies

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

Co-infection with porcine bocavirus and porcine circovirus 2 affects inflammatory cytokine production and tight junctions of IPEC-J2 cells

Detection of a bocavirus circular genome in fecal specimens from children with acute diarrhea in Beijing

Development and validation of a multiplex conventional PCR assay for simultaneous detection and grouping of porcine bocaviruses

Rapid detection and grouping of porcine bocaviruses by an EvaGreen((R)) based multiplex real-time PCR assay using melting curve analysis

Development of a duplex SYBR Green based real-time PCR assay for detection of porcine epidemic diarrhea virus and porcine bocavirus3/4/5

Porcine bocavirus: achievements in the past five years

Genetic analysis of three porcine bocaparvoviruses and identification of a natural recombinant breakpoint in NS1

A TaqMan-based real-time PCR assay for the detection of ungulate bocaparvovirus 2

Acknowledgements The author would like to acknowledge the sup-