key: cord-0727582-ekwcrzv0 authors: Gunn, Lynda; Collins, Patrick James; Fanning, Séamus; McKillen, John; Morgan, John; Staines, Anthony; O'Shea, Helen title: Detection and characterisation of novel bocavirus (genus Bocaparvovirus) and gastroenteritis viruses from asymptomatic pigs in Ireland date: 2015-06-09 journal: Infect Ecol Epidemiol DOI: 10.3402/iee.v5.27270 sha: 07241235ee553a5b165ed19b2091101128eca805 doc_id: 727582 cord_uid: ekwcrzv0 BACKGROUND: Livestock animals have been the assumed source of several human epidemics in recent years, for example, influenza H1N1, rotavirus G8/G9, and MERS-CoV. Surveillance of novel viruses in animals is essential to evaluate the risk to human and animal health and to determine any economic impact, for example, failure to thrive. There is a paucity of data regarding detection and characterisation of gastroenteritis viruses, particularly novel viruses, in porcines in Ireland. Recently, a number of small novel porcine DNA viruses have emerged globally, for example, torque teno sus virus, porcine bocavirus, and parvoviruses 2 & 4, and little is known about the biology and potential pathogenicity of these viruses. Bocaparvovirus is a genetically distinct group of viruses which has been recently detected in humans and animals. METHODS: In this study, the presence of gastroenteritis viruses (rotavirus A, porcine circovirus, adenovirus, and porcine bocavirus) was investigated in a selection of archived faecal samples from asymptomatic piglets from a commercial farm in Ireland. A total of 104 specimens were pooled and screened using conventional molecular techniques (PCR and RT-PCR), a subset of specimens (n=44) were then examined individually. Viral diversity was then investigated using statistical and phylogenetic techniques. RESULTS: Initial screening showed a high prevalence of PBoV in this farm, with the formation of three distinct groups in phylogenetic analysis. Other viruses were also investigated in this study with the first report of PCV, PAdV and lineage I G5 RVA in Ireland. Some specimens contained >1 virus, with statistical analysis indicating a strong correlation for mixed infections of PBoV and PAdV on this farm. CONCLUSION: Investigating the diversity of circulating enteric viruses on Irish porcine farms is important to improve the prophylactic tools available and to facilitate the early detection of changes in circulating viruses. and infertility (SMEDI) in pigs (PPV). Bocaviruses (BoV) are a recently described group of DNA viruses; members are characterised based on their genetic homology to bovine parvovirus and canine minute virus. BoVs were first described in humans (3) and have also been detected in canine, feline, and porcine hosts (4) . They are a genetically diverse group, distinct from the other members of Parvovirinae as they contain three major open reading frame (ORFs), coding for four genes (NS1, NP1, and VP1/2) (5, 6) . The presence of ORF3 is a near exclusive trait of bocavirus species (shared with PPV4), coding for NP1 protein, which is essential in viral propagation (7, 8) . PBoV was first detected by Blomströ m et al. (9) , in piglets suffering from post-weaning multi-systemic failure with porcine circovirus 2 (PCV2) aetiology. There have been many studies to characterise this group of viruses in pigs, and, to date, five putative species of PBoV species have been described, but classification of these isolates is poorly defined (10) . The ninth ICTV report (International Committee of Taxonomy of Viruses) defined a new bocavirus species as any isolate with less than 95% identity in non-structural gene sequences (11) . However, with the surge of sequence data for this virus family, classification of isolates has become confused. For consistency, PBoV isolates will be discussed as the groups as outlined by Yang et al. and Xiao et al. (10, 12) , throughout this report. In addition, recently proposed rationalisation and extension of the taxonomy of the family Parvoviridae (www.ictv online.org/virusTaxonomy.asp) is also used (13) . The pathogenicity of PBoV remains unclear, with viral isolates being detected in both diseased and healthy animals (14Á16). In addition, viral DNA can be detected in a range of tissues, including lymph nodes (17) , serum (18) , lung (19) , saliva (20) , and faecal matter (21) . Furthermore, PBoV can be found in a high percentage of mixed infections leading to the question as to whether it is a true pathogen or opportunistic in nature, requiring activation through co-infection (12, 22Á24) . Most studies have used molecular techniques to detect and characterise PBoV (17Á19, 21, 25Á27), while some limited progress has been made using cell culture and immunological techniques for group 3 PBoV (28, 29) . Other gastroenteritis viruses such as group A rotavirus (RVA), porcine adenovirus (PAdV), and PCV are economically important viruses. PCV and PAdV have been associated with gastroenteritis, respiratory disease, or post-weaning multi-systemic wasting (PWMS) in piglets. RVA in pigs is attributed to severe gastrointestinal disease and accounts for 89% of diarrhoea in commercial pig farms (30, 31) ; in addition RVA in pigs has been recognised to be both enzootic and epizootic, and many of the human RVA gene segments share a common origin with porcine RVA (32, 33) . In this study, the presence of bocavirus and other gastroenteritis viruses was investigated in archived faecal samples from a commercial pig farm in Ireland. Results highlight the diversity of gastroenteritis viruses present, also the diversity of BoV isolated from just one farm in Ireland and the detection of more than one BoV group in one host. In addition, other viruses such as adenovirus (PAdV), PCV1, and rotavirus A (RVA) were also detected and characterised. This is the first report of group 1 and 2 PBoV, PCV, and PAdV in Ireland. Sample preparation and extraction Archived faecal specimens (n0104) which had been collected from a commercial farm in Ireland in 2007 and stored at 48C were examined. Samples were taken from asymptomatic suckling and weaning piglets, aged from 2 to 6 weeks. Specimens were pooled prior to extraction and analysis (n026). Nucleic acids were extracted from pooled specimens using Qiaamp DNA stool mini kit (Qiagen), following manufacturer's instructions. Viral nucleic acid was eluted in 100 ml of elution buffer and stored at (208C prior to analysis. For analysis of RVA in pooled specimens, nucleic acid was extracted using phenol-chloroform, with ethanol precipitation (30) . Viral RNA was re-suspended in 100 ml of DEPC H 2 0 and stored at (208C. PBoV was detected using three sets of primers targeting the VP1/2 region of PBoV; primers from Shan et al. (34) were used to detect PBoV groups 1 and 2, while primers designed during this study were used to detect PBoV group 3 (Forward Á 5? Á AAA TTG CGC CTG CGC TCA ACG Á 3? and Reverse Á 5? Á TGC GTC CAA GGA AAG GCG TG Á 3?). To characterise isolates, primers which amplify the NP1 gene were selected; previously published primers (14) were used for group 1 PBoV; those reported by Choi et al. (20) were used to characterise group 2 PBoV, while primers designed in this group were used for group 3 PBoV (Forward Á 5? Á CAT ACT ACT ACC AGC GAC GG Á 3? and Reverse Á 5? Á CTC GTC GGC TTT ATT GAC AGG Á 3?). The PCR reactions were carried out in 50 ml volumes using the following reagents: 10 ml 5) Go Taq Flexi buffer † , 3 ml of 25 mM MgCl 2 , 8 ml 1.25 mM dNTP mix, 1 ml of each primer (50 pmol), 0.5 ml of 5 U/ml of Taq Polymerase (Promega, USA), and 5 ml of template. Selected positive pooled samples were subsequently re-extracted using individual specimens, to further characterise and determine the genetic diversity of gastroenteritis viruses, in particular BoV, therein. Specimens were also tested for the presence of other DNA viruses, namely adenovirus (35) and PCV types 1 and 2 (36), using the conditions described above. The presence of RVA was also investigated using primers specific to VP7 and VP4 genes (30, 37) . A one-step RT-PCR reaction was used for RVAusing a similar reaction mix to that described previously: 10 ml 5) Go Taq Flexi buffer † , 3 ml of 25 mM MgCl 2 , 8 ml 1.25 mM dNTP mix, 1 ml of each primer (50 pmol), 0.5 ml of 5 U/ml Taq Polymerase (Promega, USA), 0.2 ml of 10 U/ml AMV-RT (Promega, USA), and 3 ml of template. All reactions were carried out on a Biometra T3000 thermocycler. Amplified products were visualised on a 1.5% agarose gel stained with ethidium bromide using a UV transilluminator. Sequencing and data analysis Selected positive samples were cleaned using Roche High Pure PCR clean kit (Roche) and sequenced using a commercial service (GATC, Germany). For PBoV, seven specimens were sequenced; initially the VP1 gene was used in detection, while the NP1 gene was selected for sequencing as it is unique to BoV species. For group 1 PBoV, one specimen was sequenced ( Â1,040 bp) as this group typically displays little genetic diversity; likewise, two specimens were sequenced for group 2 PBoV (Â820 bp) and four specimens were sequenced for group 3 PBoV (Â1,000 bp), as these groups are more genetically diverse. For other gastroenteritis viruses (PCV, PAdV, and RVA), capsid genes were targeted, as these genes are more useful for detection purposes. Representative specimens were sequenced for PCV1 ( Â900 bp) and PAdV (Â610 bp), as only one PCV1 positive was detected and PAdV3 has limited genetic diversity. For RVA, two specimens were sequenced to determine the G and P types, using VP7 (Â1,070 bp) and VP4 (Â876 bp) RT-PCR products, respectively, to examine and compare diversity of RVA to previous studies in Ireland. Resulting sequence data were subsequently analysed and edited using BioEdit v7.0.9.0 (38) and online BLAST tool (www.blast.ncbi.nlm.nih.gov/Blast.cgi). Alignment of BoV isolates was performed using MUSCLE, and the resulting alignment was used for phylogenetic analysis. The phylogeny for BoV was performed using maximum likelihood criterion as implemented in MetaPIGA v3.1 (39) . The dataset was evaluated for saturation and ambiguous sequences; furthermore, excessively gapped regions were trimmed using the TrimAl (40) tool available in MetaPIGA. Bayesian Information Criterion was used to select the best model for the dataset; GTR model with rate heterogeneity and 1,000 bootstrap replicates was implemented and values less than 70% were collapsed. Analysis of other viral isolates was performed using Clustal W alignment as implemented in BioEdit; the subsequent phylogenies were constructed using Neighbourjoining method as executed in Mega 6 (41), using maximum composite likelihood model and 1,000 bootstrap replicates. Taxa are displayed with accession number, virus and genotype, strain name, country of origin, and year of isolation, provided the data are available on GenBank. Isolates from this study are highlighted with a filled circle (m). In this study, 104 archived faecal specimens from a pig farm in Ireland were investigated for the presence of bocavirus and other selected gastroenteritis viruses. Initially, specimens were pooled (n026, pools of 4) prior to analysis and results of initial screening were as follows: (Tables 1 and 2 ). Subsequently, individual positive specimens were re-analysed and further characterised from Table 1 . Results of initial screening of pooled samples have single PBoV infections. The frequency of mixed infections, including RVA, PAdV, group 2 or 3 PBoV, in individual specimens is represented using a mosaic plot (Fig. 1) , correlation disparity is measured using Pearson residuals. The phylogeny for PBoV isolates (Fig. 2 ) reveals three major clades, similar to those described in Refs. (10) and (12) . Group 1 and 3 PBoV form distinct monophyletic groups, in comparison to group 2 PBoV, which forms a paraphyletic group within a large monophyly, including other animal and human bocaviruses. Group 2 PBoV 5 were more related to other animal BoV than to HBoV. Subgroups for group 3 PBoV have been proposed in previous publications (10, 12) , but, depending on the gene being investigated, these subgroups do not consistently cluster together. Phylogenies for other viruses were also determined; Figure 3 indicates that the PCV1 isolated in this study is closely related to the Hungarian wild boar isolate WB-H-8, although the corresponding bootstrap value does not support this relationship (B70%). Figure 4 reveals that PAdV from this study is related to PAdV 3, similar to the reference strain IAF and wild-type isolates from Spain. Two RVA isolates were sequenced to determine the G and P types, these isolates were genotyped as G5P [13] and G9P [6] (Figs. 5 and 6 ). Phylogenetic analysis revealed greater genetic diversity of G5P [13] than in the G9P [6] strain. For G9P [6] , the VP7 gene tree shows that G9 has a strong relationship with the Irish porcine RVA from a previous study; clustering within lineage VI, the VP4 gene tree indicates that P [6] is similar to other porcine and human strains within lineage I. For G5P [13] , the tree shows more marked genetic diversity within the VP7 gene tree, G5 clusters within lineage I, a different lineage to previously described Irish isolates, and P [13] also displays greater diversity to previously isolated Irish strains. In this study, archived porcine faecal specimens were investigated for the presence of viruses proven or implicated in having a role in gastroenteritis, including the novel DNA virus, PBoV. Other porcine viruses (rotavirus A, porcine circovirus, adenovirus) were also investigated as there is still ongoing debate as to whether PBoV is a true pathogen or an opportunistic infection (22, 42) . In other animals, there is increasing evidence that BoV is a true pathogen found in respiratory and gastrointestinal disease cases, but the pathology is still unknown (18, 20, 42, 43) . In this study, 48% of PBoVpositive specimens contained more than one group of PBoV, while 48% were positive for other viruses. The most commonly detected additional virus was AdV (48%), followed by PCV and RVA, which were each detected in 4% of PBoV-positive specimens. A mosaic plot representing the frequency of mixed infections including RVA, PAdV, group 2 or 3 PBoV was constructed (Fig. 1) , Pearson residuals suggest that there is a strong correlation between PBoV and PAdV, but further large-scale investigations are necessary to confirm this relationship. As only a subset of the total number of specimens were examined individually using PCR/RT-PCR, these figures are not completely representative but are indicative of the wide diversity of gastrointestinal viruses detected in the samples. Since the discovery of PBoV in 2009, five putative species have been identified through the use of ICTV guidelines (B95% homology) (11) . Members of genus parvovirus have been suggested to have mutation rates similar to that of RNA viruses (44, 45) ; in addition, BoV have been noted to undergo a high incidence of genetic recombination (5, 21) . Considering these factors, a B95% homology value was too low for accurate identification of new species. Recent release of ICTV guidelines stipulate new Parvoviridae species are identified by B85% homology of the NS1 gene. In addition, the nomenclature of taxonomic groups has also been changed and this has also been incorporated in this manuscript. Typically, three major monophyletic groups can be identified as outlined in previous reports 10, 12); group 1 contains isolates from species 1, including original isolate Swebo; group 2 consists mainly of species 2 isolates, and group 3 contains species 3Á5; in addition, group 3 contains viruses from the same species that do not cluster together. In the BoV phylogeny displayed here (Fig. 2) , three major monophyletic groups can be identified; two of these groups are exclusively PBoV groups 1 and 3, containing no other human or animal isolates. Group 2 PBoV, however, is part of a large monophyletic group which contains other animal and human BoV; in addition, PBoV in this group are more closely related to other animal BoV than to HBoV (Fig. 2) . In general, Irish isolates cluster together with isolates previously described globally, displaying a high degree of diversity. For group 1 PBoV, only one isolate was sequenced; as seen from previous studies, there is very low genetic diversity in this group. CIT-B1.61 has 98.4Á99.8% sequence identity to other isolates in group 1 PBoV. Group identity to each other; overall, isolates within this group share 87.8Á96% identity. There was greater diversity observed in group 3 PBoV in comparison. Irish strains share 73.3Á83.3% sequence identity; overall, group 3 PBoV has 73Á98.5% identity to global isolates within this group. Group 3 PBoV has been previously described in Northern Ireland (28) and isolates CIT-B3.98 and CIT-B3.57 appear to be more closely related to these isolates; isolate CIT-B3.98 has 95.4% identity to the Irish isolate 64-1 (Accession no. JF512472), while CIT-B3.57 has 82.2% identity to isolate F41 (Accession no. JF512473). Many of the PBoV-positive specimens described here contain either more than one group of PBoV or other potential pathogens. There is ongoing debate as to whether BoV is a true pathogen or a passenger in which co-infection is required for replication (12, 22, 23) . The original report of PBoV (group 1) was in the presence of piglets that died from PWMS, with PCV2 aetiology (9). Our study describes PCV1 with group 3 PBoV (CIT-20). Although PCV1 is not considered pathogenic, there is some limited data which suggest that PCV1 may have potential to cause disease (46) . RVA and PBoV have been reported previously in an investigation into co-infection and porcine diarrhoeal disease (24) , but the isolates were not genotyped. In our study, one specimen (CIT-43) was positive for both RVA (G9P [6] ) and group 3 PBoV. RVA isolated from this study was for the most part similar to isolates previously isolated from Ireland (Figs. 5 and 6) Lynda Gunn et al. (47) , except for isolate CIT-53, which is the first report of lineage 1 G5 rotavirus in Ireland. To our knowledge, this is the first report of group 1 and group 2 PBoV in Ireland; additionally, PBoV in the presence PAdV. The majority of PBoV and PAdV 3 co-infections were with group 3 PBoV (44%). The issue of diarrhoeal coinfections in pigs has been investigated previously by Zhang et al. (24) , who reported that 75.1% of clinically diseased animals had co-infections with up to four pathogens. In our study, PBoV was isolated from asymptomatic animals, which has also been noted in previous studies (14Á16). The 2012 All-island Surveillance Report in Ireland on a wide variety of animals including bovines, porcines, equines, and ovines reported that enteritis was among the leading causes of mortality in pigs (48) . Currently, there is no vaccine licensed for the prevention of viral enteritis in piglets in Ireland and given the broad genetic heterogeneity of the strains detected in this study (and previous studies in Ireland), it remains unclear whether this diversity may pose a challenge for future prophylaxis programmes for the prevention of enteritis in piglets. This study, along with other previous studies (49) , indicates that a wide variety of enteric viruses are circulating in piglets in Ireland; therefore, it will be important to establish surveillance for such viral pathogens, including porcine bocavirus, circovirus, and adenovirus, to gather information on the genetic diversity of these viruses and to plan more adequate measures of control for porcine enteritis. A clear understanding of viral diversity is still hampered by the fact that enteric viruses are not included routinely in the diagnostic tests of porcine enteritis, and the information is limited. It is likely that porcine enteritis has a multi-factorial nature, with viruses and other micro-organisms acting in synergy, yet under the influence of other factors, such as environmental and social stress due to overcrowding and/or incorrect management. Pigs have been an assumed reservoir for previous human epidemics such as influenza H1N1 and rotavirus G9 infections. Bocavirus/Bocaparvovirus is a recently described virus genus which can be detected in humans and animals. Although animal BoVs are relatively distinct from HBoV, PBoV does appear to share greater genetic homology with HBoV, based on NS1 gene (50) . PBoV also appears to be much more genetically diverse in comparison to other animal and human BoV. To date, there have been no reported cases of BoV zoonosis. However, monitoring and surveillance of pathogens which can be isolated from humans and animals is important to understand the risks of zoonosis, approaches to vaccination, and implications for both human and animal health. RVA/Pig-wt/Can/F7P4-A Neighbour-joining phylogeny of rotavirus A VP4 gene for genotypes Confronting zoonoses, linking human and veterinary medicine Viral host jumps: moving toward a predictive framework Cloning of a human parvovirus by molecular screening of respiratory tract samples Animal bocaviruses: a brief review Human Bocaviruses are highly diverse, dispersed, recombination prone, and prevalent in enteric infections A newly identified Bocavirus species in human stool The nonstructural protein NP1 of human bocavirus 1 induces cell cycle arrest and apoptosis in Hela cells Human bocavirus NP1 inhibits IFN-b production by blocking association of IFN regulatory factor 3 with IFNB promoter Detection of a novel porcine boca-like virus in the background of porcine circovirus type 2 induced postweaning multisystemic wasting syndrome Genome characterization of a novel porcine bocavirus Virus taxonomy: classification and nomenclature of viruses: ninth report of the International Committee on Taxonomy of Viruses Molecular evolutionary genetic analysis of emerging parvoviruses identified in pigs Rationalization and extension of the taxonomy of the family Parvoviridae. International Committee on Taxonomy of Viruses 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 Emerging and re-emerging swine viruses Detection of novel porcine bocaviruses in fecal samples of asymptomatic pigs in Cameroon Genetic characterisation of a porcine bocavirus detected in domestic pigs in Uganda High prevalence of a novel porcine bocavirus in weanling piglets with respiratory tract symptoms in China Detection, prevalence and analysis of emerging porcine parvovirus infections Molecular detection and genetic analysis of porcine bocavirus in Korean domestic swine herds 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 Human bocavirus: passenger or pathogen in acute respiratory tract infections? Human bocavirus: lessons learned to date Occurrence and investigation of enteric viral infections in pigs with diarrhea in China Identification and nearly full-length genome characterization of Novel Porcine Bocaviruses Complete genome sequence of a novel species of porcine bocavirus, PBoV5 Complete sequence and phylogenetic analysis of a porcine bocavirus strain swBoV CH437 Isolation in cell cultures and initial characterisation of two novel bocavirus species from swine in Northern Ireland Production, characterisation and applications of monoclonal antibodies to two novel porcine bocaviruses from swine in Northern Ireland Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens Evaluation of rotavirus infection and diarrhea in Iowa commercial pigs based on an epidemiologic study of a population represented by diagnostic laboratory cases Zoonotic aspects of rotaviruses Full genome-based classification of rotaviruses reveals a common origin between human wa-like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains Genomic characterization and high prevalence of bocaviruses in swine Detection of bovine and porcine adenoviruses for tracing the source of fecal contamination The emergence of a new strain of porcine circovirus-2 in Ontario and Quebec swine and its association with severe porcine circovirus associated disease Á Identification of group A rotavirus gene 4 types by polymerase chain reaction BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT MetaPIGA v2.0: maximum likelihood large phylogeny estimation using the metapopulation genetic algorithm and other stochastic heuristics trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses MEGA6: molecular evolutionary genetics analysis version 6.0 Human bocavirus: increasing evidence for virulence Human bocavirus: still more questions than answers Rates of evolutionary change in viruses: patterns and determinants High rate of viral evolution associated with the emergence of carnivore parvovirus Outcome of experimental porcine circovirus type 1 infections in mid-gestational porcine foetuses Detection and characterisation of group A rotavirus in asymptomatic piglets in southern Ireland All-island Animal Disease Surveillance Report Molecular characterization of equine rotavirus in Ireland Complete coding sequences and phylogenetic analysis of porcine bocavirus The authors have not received any funding or benefits from industry or elsewhere to conduct this study.