Next Article in Journal
Generation and Evaluation of Recombinant Thermostable Newcastle Disease Virus Expressing the HA of H9N2 Avian Influenza Virus
Next Article in Special Issue
Next-Generation Sequencing Reveals Four Novel Viruses Associated with Calf Diarrhea
Previous Article in Journal
Retrospective Study of the Upsurge of Enterovirus D68 Clade D1 among Adults (2014–2018)
Previous Article in Special Issue
Whole-Genome Sequence Analysis of Pseudorabies Virus Clinical Isolates from Pigs in China between 2012 and 2017 in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Cause of Diarrhoea and Poor Growth in 8–11-Week-Old Pigs from an Australian Pig Herd Using Metagenomic Sequencing

by
Tarka Raj Bhatta
1,2,
Anthony Chamings
1,2 and
Soren Alexandersen
1,2,3,*
1
Geelong Centre for Emerging Infectious Diseases, Geelong, VIC 3220, Australia
2
School of Medicine, Deakin University, Geelong, VIC 3220, Australia
3
Barwon Health, University Hospital Geelong, Geelong, VIC 3220, Australia
*
Author to whom correspondence should be addressed.
Viruses 2021, 13(8), 1608; https://doi.org/10.3390/v13081608
Submission received: 22 July 2021 / Revised: 6 August 2021 / Accepted: 11 August 2021 / Published: 13 August 2021
(This article belongs to the Special Issue Applications of Next-Generation Sequencing in Virus Discovery)

Abstract

:
Diarrhoea and poor growth among growing pigs is responsible for significant economic losses in pig herds globally and can have a wide range of possible aetiologies. Next generation sequencing (NGS) technologies are useful for the detection and characterisation of diverse groups of viruses and bacteria and can thereby provide a better understanding of complex interactions among microorganisms potentially causing clinical disease. Here, we used a metagenomics approach to identify and characterise the possible pathogens in colon and lung samples from pigs with diarrhoea and poor growth in an Australian pig herd. We identified and characterized a wide diversity of porcine viruses including RNA viruses, in particular several picornaviruses—porcine sapelovirus (PSV), enterovirus G (EV-G), and porcine teschovirus (PTV), and a porcine astrovirus (PAstV). Single stranded DNA viruses were also detected and included parvoviruses like porcine bocavirus (PBoV) and porcine parvovirus 2 (PPV2), porcine parvovirus 7 (PPV7), porcine bufa virus (PBuV), and porcine adeno-associated virus (AAV). We also detected single stranded circular DNA viruses such as porcine circovirus type 2 (PCV2) at very low abundance and torque teno sus viruses (TTSuVk2a and TTSuVk2b). Some of the viruses detected here may have had an evolutionary past including recombination events, which may be of importance and potential involvement in clinical disease in the pigs. In addition, our metagenomics data found evidence of the presence of the bacteria Lawsonia intracellularis, Brachyspira spp., and Campylobacter spp. that may, together with these viruses, have contributed to the development of clinical disease and poor growth.

1. Introduction

Porcine enteric diseases resulting in diarrhoea or intestinal pathology and poor growth performance cause huge economic losses around the world [1,2,3,4,5,6,7,8]. It has been suggested that various etiological agents acting either individually or synergistically including viruses, bacteria, and parasites are responsible for clinical consequences such as diarrhoea and poor growth in pigs [9,10]. Viruses, including picornaviruses such as porcine sapelovirus (PSV), porcine enteroviruses (PEV), and porcine teschoviruses (PTV), have been identified in association with various disorders, including diarrhoea, polioencephalomyelitis, respiratory distress, dullness, skin lesions, pyrexia, and flaccid paralysis [11,12,13,14,15]. Porcine astrovirus (PAstV) has also been associated with gastrointestinal, respiratory, and neurological disorders [16,17,18]. Porcine circovirus type 2 (PCV2) has been associated with a number of diseases/clinical syndromes together called porcine circovirus associated diseases (PCVAD) [19], including postweaning multisystemic wasting syndrome (PMWS) and porcine dermatitis and nephropathy syndrome (PDNS) [20,21]. Some of the clinical sign of PMWS include growth retardation, enlarged lymph nodes, and wasting [22]. It has been shown that co-infection of PCV2 with porcine parvoviruses (PPVs), porcine bocaviruses (PBoVs), and torque teno sus viruses (TTSuVs) can be associated with porcine respiratory diseases complex (PRDC) [23] and PMWS [24]. Bacteria such as Lawsonia intracellularis [25,26,27], Brachyspira spp. [28,29,30], Campylobacter spp. [31], Escherichia coli, and Salmonella spp. [32] have also been found to be associated with diarrhoea and poor growth in pigs. Culture and targeted polymerase chain reaction (PCR) can be performed for the detection of these bacteria [33,34]. Detection of specific porcine viruses known or likely to be present, e.g., swine influenza A virus (swIAV) [35], PSV [36], PEV [37], PTV [38], PAstV [39], PBoV [40,41], PCV2 [42], and many others, can also be performed using targeted PCR assays. However, targeted detection may not be suitable for investigating co-infection with multiple pathogens, or if there is an involvement of novel or variant viruses and other pathogens [43].
Next generation sequencing (NGS) technology is a very useful tool for the metagenomic detection and characterisation of multiple, unexpected, novel, or unidentified viruses from healthy and diseased individuals [44,45,46,47]. Various independent studies have used NGS for the detection of diverse groups of viruses in pig samples including PEV, PAstV, PCV2, atypical porcine pestivirus (APPV), PBoV, TTSuVs, and many others [47,48,49,50]. Similarly, different pathogens were also identified in bat [51], dog [44,52], bird [45,46,53], rodent [54], and human samples [55,56] using NGS technology.
Previously we have used our metagenomics approach for the detection and characterisation of various viruses (both RNA and DNA) in dogs suffering from diarrhoea [44] and in faecal samples from wild birds [45,46] and humans [45]. The study described here focused on the metagenomic detection and characterisation of pathogens from a single pig herd located in South Eastern Australia that was experiencing diarrhoea and poor growth in growing pigs from eight weeks of age onwards.

2. Materials and Methods

2.1. Clinical History and Sample Collection

The affected pig herd was located in South Eastern Australia. The farm produced and recruited gilts internally but did use commercially bought local boar semen as part of their breeding program. There were no recent introductions of pigs to this farm; the last introduction being 10 years prior from a genetic herd in Queensland, Australia. Progeny pigs were not vaccinated against PCV2 or any other pathogens. When visited by the veterinarian in late August 2018, approximately eight to twelve pigs from a total of 52 pigs in a pen had started showing relatively severe signs of gastrointestinal disease/diarrhoea from approximately eight weeks of age. A few pigs had died while some were skinny and not growing well compared to other pigs in the pen. Treatment was initiated, including Tylosin injection (Tylan, Elanco) which resulted in no noticeable improvement. The concentration of Tylosin in the feed was increased from 50 to 200ppm with a slight beneficial effect. The diet was also supplemented with 100ppm Olaquindox as a growth promoter. The pigs were then injected with 10mg/kg Engemycin (Oxytetracycline), a broad spectrum antibiotic active against both Gram-positive and Gram-negative bacteria, which resulted in some clinical improvement in some of the pigs, but the problem persisted to eleven weeks of age. At eleven weeks, (mid-September 2018), two (Pig 45 and Pig 46) of the affected pigs in the pen were euthanised for post-mortem examination by the attending farm veterinarian. Gross findings included mild proliferative lesions in the small and large intestines. The colon was enlarged and there was some inflammation of the colon mucosa along with excessive mucus. Proliferative enteropathy (Lawsonia intracellularis) was suspected, and samples were submitted for histopathology along with culture for Salmonella and Brachyspira spp. at a veterinary diagnostic laboratory. In addition, fresh samples of surplus colon and lung tissues from the two pigs were forwarded overnight to the Geelong Centre for Emerging Infectious Diseases (GCEID) in an Esky with ice packs to perform a metagenomics investigation using NGS. Intestinal and respiratory tissues were chosen as these systems are common sites of virus excretion. Upon receipt, swabs were taken from each of the mucosal epithelial surfaces of the colon and multiple sites within the lung tissues from each pig, and the swabs were placed in universal transport medium (UTM). Aliquots of the tissue suspensions were frozen at minus 80 °C until processed further.

2.2. Enrichment, Nucleic Acid Extraction, cDNA Synthesis and Non-Targeted Amplification

For NGS, each colon and lung sample from each pig were processed separately, and small particle-protected nucleic acids enriched using a previously described protocol optimised in our laboratory with some modifications [44,45,57]. Briefly, the samples were homogenised at 25 Hz for 2 min using a TissueLyzer (Qiagen, Hilden, Germany) followed by centrifugation at 17,000× g for 3 min and then filtered using a 0.8 µm polyether sulfone (PES) spin-column filter. The samples were then ultracentrifuged at 178,000 g for 1 h. Pelleted small particles were then treated with benzonase and micrococcal nuclease for 2 h to enrich for nucleic acids being protected in the particles [44]. The ultracentrifuged and nuclease treated samples were then treated with 200 µM propidium monoazide (PMAxx) for 30 min in a PMA-Lit LED photolysis device to cross-link nucleic acids not protected in particles, as indicated by the manufacturer (Biotium, Fremont, CA, USA) [57]. Finally, total nucleic acids were extracted using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany). For cDNA synthesis and DNA amplification, an initial RNA denaturation step of 95 °C for 3 min followed by snap-cooling in an ethanol–ice bath (minus 20 °C) was carried out before using the SeqPlex RNA Amplification Kit (Sigma, St. Louis, MO, USA), as per the manufacturer’s protocol [44,45]. Bioanalyzer and high-sensitivity DNA chips (Agilent, Waldbronn, Germany) were used to check the quantity and quality of the amplified product.

2.3. Library Preparation, Next Generation Sequencing (NGS), and Detection of Different Virus Sequences

Library preparation was performed using the Ion Plus Fragment Library Kit and IonXpress Barcode Adapters 1-96 Kit (Thermo Fisher Scientific, Waltham, MA, USA). The libraries were quantified using the Ion Library TaqMan Quantitation Kit (Thermo Fisher Scientific, Waltham, MA, USA) and then pooled. Sequencing was performed using Ion 530 chips and an Ion Torrent S5XL System (Thermo Fisher Scientific, Waltham, MA, USA). Each of the pooled individual samples generated approximately 2.8 to 6.9 million reads. Generated sequence reads were compared to a local database of virus sequences downloaded from GenBank (December 2018) using BLASTN and BLASTX [58,59] queries with an e-value cut-off score of 1 × 10−10. Contigs were generated using the AssemblerSPAdes 5.12.0.0 plugin on the Ion Torrent Server and were compared to the virus reference sequences in NCBI GenBank using BLASTN and BLASTX to identify virus sequences. NGS reads were mapped to the contigs matching virus sequences using the TMAP plugin [60] to determine the coverage depth. The TMAP plugin [60] from the Ion Torrent suite was also used for mapping the NGS reads using selected virus sequences from NCBI GenBank that were identified as most similar to the viruses in the pig samples. Consensus sequences were generated based on the number of high-quality reads (initially MAPQ >20) matched to the reference viruses and coverage/abundance of the reads across the reference genomes.
Near full or partial consensus sequences were generated from the mapped sequences using Integrative Genomics Viewer software (IGV) (Broad Institute, MA, USA) [61] with a mapping quality score of 64 or higher and a coverage depth of at least two with further specific details noted in the results below. Generated sequences were compared with the online NCBI GenBank database using BLASTN [58,59], and the closest, and additional related or representative sequences from around the globe, were selected for further analyses. Selected sequences were aligned using Clustal-W [62] and phylogenetic trees, created using the maximum likelihood (ML) method with the best fitting model, as determined by MEGA 7 [63]. The robustness of different phylogenetic nodes was assessed using 1000 bootstrap replicates for nucleotide sequences.

2.4. Similarity Plot (SimPlot) Analysis of Identified Virus Sequences

Similarity plot (SimPlot) was used to look for possible recombination between our virus sequences obtained from the pig samples and related sequences from NCBI GenBank. A 200 nucleotide sliding window at 20 nucleotide intervals and the F84 distance model with maximum likelihood method was used. Percentage identities at each analysis point were plotted on a line chart [64].

2.5. Mapping of Bacterial and Antimicrobial Resistance (AMR) Gene Reads

Although our NGS method enriched for virus particles, the final NGS reads also includes some reads derived from other nucleic acids, e.g., from ribosomes, mitochondria, and bacterial or host genomes [45]. Consequently, mapping of NGS reads to bacterial DNA was performed using CLC Workbench version 20.0 Microbial Genomics Module (Qiagen, Hilden, Germany) with the included “Taxonomic Profiling” tool and the microbial genome database from January 2020. Mapping of reads to the 16S ribosomal RNA microbial databases was also done on the Thermofisher Ion Reporter website, using the setting of at least two identical reads mapping over at least 60 nucleotides (60-60-2). Additionally, a BLASTN query with an e-value cut-off score of 1 × 10−10 was performed against the refseq ribosomal (5S, 16S and 23S) database from the NCBI GenBank genetic sequence database (downloaded September 2019). Finally, CLC and the Microbial Genomics Module were used to map against the CLC antimicrobial resistance (AMR) gene databases using CLC assembled contigs and the tool “Find resistance with nucleotide database”.

3. Results

NGS was conducted on swab samples taken from fresh chilled samples of surplus colon and lung tissues from two pigs experiencing diarrhoea and poor growth at 8–11 weeks of age. Based on the initial BLASTN and BLASTX analyses of the NGS reads and their abundance, both RNA and DNA viruses representing 0.07–1.37% of the total sequenced reads were identified in samples from both pigs (Table 1). RNA viruses identified were picornaviruses: including porcine sapelovirus (PSV), enterovirus G (EV-G), porcine teschovirus (PTV), and a porcine astrovirus (PAstV). Single Stranded DNA viruses identified were parvoviruses: porcine bocavirus (PBoV), porcine parvovirus 2 (PPV2), porcine parvovirus 7 (PPV7), porcine bufa virus (PBuV), porcine adeno-associated virus (AAV), porcine circovirus type 2 (PCV2), and torque teno sus viruses (TTSuVk2a and TTSuVk2b).
Despite employing a virus enrichment method to remove/degrade much of the non-virus nucleic acids, between 0.30–1.29% of the total NGS sequenced reads were mapped to bacterial chromosomal or plasmid DNA and a small number of reads mapped to the corresponding bacterial ribosomal RNA in each of the samples from both pigs (Supplementary Table S1). These reads mapped to a number of different bacterial species, but of particular interest were NGS reads mapping to pathogenic bacteria which are known to cause diarrhoea and poor growth in pigs (e.g., Lawsonia intracellularis, Brachyspira spp., and Campylobacter spp.).
Conventional diagnostic methods (histopathology and bacterial culture) detected Lawsonia intracellularis and Brachyspira spp. The NGS metagenomics sequencing confirmed the presence of these agents and was able to detect sequences from at least 12 viruses (Table 1) that were not tested for by conventional diagnostics methods. Nine of these viruses were of particular interest because of either their relatively high abundance and/or their recognised potential to cause clinical signs consistent with those observed in the herd. These viruses included porcine sapelovirus (PSV), enterovirus G (EV-G), porcine teschovirus (PTV), porcine astrovirus (PAstV), porcine bocavirus (PBoV), porcine parvovirus 2 (PPV2), porcine circovirus type 2 (PCV2), and torque teno sus viruses (TTSuVk2a and TTSuVk2b) (see below).

3.1. Porcine Sapelovirus (PSV)

Sequences belonging to PSV made up between 0.02 and 1% (Table 1) of total sequenced reads from the tissue samples. Four and six partial consensus sequences representing different regions of at least two genetically distinct PSVs were assembled from the mapped reads using a mapping quality of 80–90 and coverage of 2–5638 from the swabs of Pig 45 and Pig 46, respectively (Supplementary Table S2). These assembled sequences each included the full sequence of the virus capsid structural proteins (VP4, VP2, VP3, and VP1), the full sequence of the non-structural virus proteins P2 (2A-2B-2C), and a partial sequence of the non-structural proteins P3 (3A-3B-3C (but not 3D)) of PSV except for sapelovirus sequenced in the Pig 46 lung sample.
Sequence comparison of the VP1 region has been used for the genotyping of several picornaviruses as this is one of the most variable capsid proteins and important in inducing specific host immunity [65]. In the case of sapelovirus, the VP1 sequence is 844 nt long and encodes for a 281 amino acid (aa) long protein. Two VP1 sequences from the colon and lung of Pig 45 were 100% identical with each other and were also 100% identical to the partial VP1 sequence from the lung sample of Pig 46. The VP1 sequences from the colon and lung samples of Pig 45 were only 80.2% identical with the closest PSV sequence (LC508229-PSapV/22-B/Zambia-2018) from Zambia and were only 74.6% identical with the VP1 sequence obtained from the colon sample from Pig 46, indicating that different porcine sapeloviruses were circulating on the farm and were possibly co-infecting Pig 46. The VP1 sequence from Pig 46 colon was 88.9% identical with the closest PSV sequence (MK378907-PSapV/GX04 C1/China-2017) from China. Interestingly, the VP1 sequence from the Pig 46 colon sample has a 21 (ACAGGTTTCTACCCTGCAACC) nucleotide (7 amino acids: TGFYPAT) insertion/repeat towards the end of VP1, while this is absent in the sapelovirus sequences from the colon and lung of Pig 45 and that from the lung sample of Pig 46. The insertion has been reported for some other sapelovirus and other picornavirus sequences as “protein of unknown function (DUF3724)” and is a short domain of approximately 20 aa in length in which a completely conserved tyrosine (Y) amino acid residue which may be functionally important [66].
Pairwise comparison of the 1827 nt long (609 aa) VP4, VP2, and VP3 (VP4-VP3) sequence of the capsid region and the 3216 nt long (1072 aa) P2-P3 (but not 3D) sequence of the non-structural region with the reference sapelovirus sequences from NCBI GenBank was conducted. VP4-VP3 sequences of Pig 45 colon and lung and Pig 46 lung swab were 100% identical with each other. These three sequences were 85.2% identical at the nt level with the closest PSV sequence LC425396-PSV/DeTk2-2/Japan-2015 from Japan and were only 81.7% identical with the VP4-VP3 sequence from the Pig 46 colon swab. KJ821019-PSV/KS04105/KOR-2004 from Korea was the closest sequence with 88.1% identity at nt level with VP4-VP3 sequence from the Pig 46 colon swab. The non-structural P2-P3 region of the PSVs from Pig 45 colon and lung were 100% identical with each other and were only one nucleotide different from the partial P2-P3 sequences of Pig 46 lung. The P2-P3 sequences from Pig 45 colon and lung were 96.6% identical with P2-P3 sequences from Pig 46 colon. All the P2-P2 sequences from this study were 89% to 89.1% identical at nt level with the closest reference, MK378928-PSV/HLJ01 C1/China-2017 from China. The PSV sequences from Pig 45 colon and lung and Pig 46 lung samples were therefore reasonably similar in their non-structural regions to Pig 46 colon sample but had a more different capsid region between nucleotides 400–3000 (Figure 1). The different capsid regions of these two viruses may have evolved through recombination or possibly immune selection of different capsid genotypes/serotypes whereas the non-structural region was less divergent.

3.2. Porcine Enterovirus G (PEV-G)

Sequences corresponding to PEV-G made up 0.02 to 0.35% of the sequenced reads from the colon and lung tissue samples from the two pigs (Table 1). Two identical complete coding sequences of PEV-G (Pig 45 colon; 6993 nt (coverage 2–3448) and Pig 45 lung; 7189 nt (coverage 2-1686)) representing a single porcine enterovirus-G were obtained from Pig 45 at mapping quality of 80. However, 15 partial PEV-G sequences (Supplementary Table S3) representing multiple porcine enteroviruses-G were obtained from the colon sample of Pig 46 and five partial consensus sequences of PEV-G were obtained from the lung sample of Pig 46 (Supplementary Table S3).
As the VP1 sequence determines the genotype of EV-G, and the criteria to distinguish/separate one genotype from another is a variation of VP1 nucleotide sequence of more than 25% [49,67], a pairwise comparison of the VP1 nucleotide sequences was conducted. We obtained four full sequences of the 729 nt long VP1 capsid of PEV-G, one from each of the lung and colon tissue samples of Pig 45 and Pig 46 and these were used for phylogenetic analysis using thirty nine reference sequences belonging to EV-G1 to EV-G20 [49,68,69] from https://www.picornaviridae.com/sg3_ensavirinae/enterovirus/ev-g/ev-g.htm (accessed on 18 March 2021) and from NCBI GenBank. The two VP1 sequences from the colon and lung from Pig 45 (MZ515508 and MZ515509, respectively), and one VP1 sequence from the lung sample of Pig 46 (MZ515525) were 100% identical with each other. These three sequences were only two nucleotides different from one of the VP1 sequences (MZ515510-EV-G-VP1-AUS-2018) from the colon sample of Pig 46. These four VP1 sequences clustered with EV-G1 genotype (Figure 2) and were closest (~80.1–80.4% identical) to KT265971-EV-G1/734141/ThanhBinh-Vietnam-2012 from Vietnam by nucleotide pairwise comparison. The other partial (363 nt long) VP1 sequence (MZ515511-EV-G-VP1-AUS-2018) from the colon of Pig 46 was most similar to the EV-G4 reference genotype and was only 69% identical with the PEV-G VP1 sequence mentioned above (MZ515510-EV-G-VP1-AUS-2018). The EV-G4 like VP1 sequence from Pig 46 colon sample was closest (~81% identical at nt level) to MF113307-EV-G4/F22-1/03-05-GER-2013 from Germany. Six near full length sequences of VP4, VP2, and VP3 (VP4-VP3) capsid region of PEV-G were obtained from Pig 45 colon and lung and Pig 46 colon samples (Supplementary Table S3) and used for phylogenetic analysis. Sequences most similar to genotypes EV-G1 and EV-G4 were again identified, as was one sequence most similar to genotype, EV-G2 and another to genotype EV-G9 (Supplementary Figure S1) both from the Pig 46 colon sample.
The non-structural P2-P3 sequences from the tissues of Pig 45 and lung of Pig 46 were most similar to the reference sequences from the EV-G1 genotype whereas sequences closest to EV-G1 and EV-G2 genotypes were obtained from the Pig 46 colon samples. No P2-P3 sequences belonging to the genotypes EV-G4 and EV-G9 were identified. Despite only obtaining partial sequences of PEV-G from Pig 46, there was clearly a number of genetically different EV’s present in the tissues of this pig at the time of sampling, although some were at very low abundance and therefore unlikely to be actively replicating at the time of sampling.

3.3. Porcine Teschovirus (PTV)

Sequences corresponding to PTV with 0.007 to 0.024% (Table 1) abundance of total sequenced reads were obtained. Four and seven partial consensus sequences representing different regions of porcine teschovirus (PTV) were assembled from the mapped reads using a mapping quality of 80 and coverage of 2–20 from the colon swabs of Pig 45 and Pig 46, respectively (Supplementary Table S4). No sequences of PTV were obtained from either of the pig lung samples.
To genotype the PTV sequences in the sample, we aligned the sequences we identified belonging to the VP1 gene sequence with the reference PTV sequences as it is considered a highly variable region containing neutralizing epitopes for differentiation of serotypes of the PTV [71,72,73]. A 610 nt long partial VP1 sequence (MZ515535-PTV-VP1-AUS-2018) encoding for 203 aa from the colon sample of Pig 46 was used for pairwise comparison and was found to belong to the PTV3 genotype with 82.2% nucleotide identity with the nearest PTV sequence, GQ293230-PTV3/Vir 1925/02/GER-2002 from Germany. Another 292 nt long partial VP1 sequence (MZ515534-PTV-VP1-AUS-2018) encoding 97 aa from the colon sample of Pig 46 (Supplementary Table S4) was found to belong to the PTV2 genotype with around 84.2% nucleotide identity with its nearest sequence AY392533-PTV2/DS 183/93/GER-2003 from Germany. Our two VP1 sequences representing PTV3 and PTV2 genotypes were only around 65% identical with each other at the nucleotide level and indicated the presence of at least two different genotypes/serotypes of PTV present in the colon of Pig 46.

3.4. Porcine Astrovirus (PAstV)

Sequences corresponding to PastV were only identified in the colon sample of Pig 46 and comprised 0.02% (Table 1) of the total sequenced reads in this sample. As a result, we were only able to generate 10 partial sequences of porcine astrovirus (PastV) from the NGS reads of the Pig 46 colon sample. The number of reads were much fewer than for PSV, PEV, PTV, and PBoV (Table 1), and therefore we were only able to generate low coverage consensus sequences of PAstV reads using a mapping quality threshold of 64 and above (Supplementary Table S5).
The ORF2 capsid region of AstV is the most genetically variable region of the genome [74] as it contains neutralizing epitopes that form the capsid protein including the spike domain [75,76]. The 1122nt (374aa) long partial ORF2 sequence (MZ515541-PAstV-ORF2-AUS-2018) was used for pairwise comparison and found to be most similar to the PAstV4 genotype and also most similar (81.2% identical at nt level) with the reference sequence LC201607-PAstV4-JPN-HgTa2-1-2-Japan-2015 from Japan. Another partial ORF2 (MZ515544-PAstV-ORF2-AUS-2018) sequence of 484 nt long was also found to be most similar to PAstV4 genotype and was closest (75.4% identical at nt level) to the reference sequence KU764486-PAstV4-15-12-USA-2015 from the USA. These two capsid sequences shared a 271 nt long overlapping region, however they were only 59.1% identical at nt level indicating that these two sequences belonged to two distinct PoAstV4 in the colon of Pig 46.
One 869 nt long partial ORF1a sequence (MZ515542-PAstV-ORF1a-AUS-2018) was obtained from the Pig 46 colon swab and was most similar (89.3% (776/869) nucleotide identity) to MK460231-PAstV4/CHN/WG-R2-China-2017 from China which belonged to PAstV genotype 4. Another 281 nt long ORF1a sequence (MZ515548-PAstV-ORF1a-AUS-2018) was obtained from the same Pig 46 colon sample which was only 75.1% (211/281) identical to our first ORF1a sequence MZ515542-PAstV-ORF1a-AUS-2018. It was therefore evident that there were at least two different PAstV4 viruses within the colon sample from Pig 46, with two different PAstV4 capsid sequences and two different PAstV4 ORF1a sequences present.

3.5. Porcine Bocavirus (PBoV)

Sequences corresponding to PBoV were identified in both the colon and lung samples from both Pig 45 and 46 and made up 0.002 to 0.225% (Table 1) of the total sequenced reads. Twenty-two partial consensus sequences representing different regions of multiple genetically distinct porcine bocavirus 3 (PBoV3) were assembled from the mapped reads using a mapping quality of 64 and above and coverage of 2–888 from the swabs of Pig 45 and Pig 46 (Supplementary Table S6).
Phylogenetic classification of parvoviruses is performed using a relatively conserved region of NS1 [77,78,79,80]. Two 1278 nt long partial NS1 sequences from the Pig 45 lung (MZ544035-PBoV-NS1-AUS-2018) and Pig 46 colon (MZ544039-PBoV-NS1-AUS-2018) samples were used for pairwise sequence comparison and phylogenetic analysis along with 21 other reference sequences representing different clusters [77,78,79,80] of PBoV3 that were chosen from NCBI GenBank database. The NS1 sequence obtained from the Pig 45 lung sample was most similar to virus sequences belonging to the PBoV3 cluster D [77,78,81,82] (Figure 3) with KF025384-PBoV/MN154-1/USA-2011 as the closest (~97.3% identical at nt level) reference sequence from the USA. However, the NS1 sequence from the colon of Pig 46 was closest (~92.4% identical at nt level) to KF360033-PBoV/swBoV CH437/China-2012 from China (Figure 3), which was proposed as belonging to a novel PBoV3 cluster [83] (here tentatively designated PBoV3 cluster F/?). These two NS1 sequences from the Pig 45 lung and Pig 46 colon swabs were only ~80.1% identical at nucleotide level with each other.
One partial 1193 nt long VP1 sequence from Pig 45 lung sample and two VP1 partial sequences each 1193 nt long from Pig 46 colon sample were created and used for pairwise sequence comparison. It was found that VP1 sequence from Pig 45 lung sample (MZ544034-PBoV-VP1-AUS-2018) was closest (90.6% nucleotide identity) to MN747338-PBoV/GXHX2017-3/China-2017 from China. However, the two VP1 sequences (MZ544041-PBoV-VP1-AUS-2018 and MZ544040-PBoV-VP1-AUS-2018) from Pig 46 colon sample were 79.4% identical at the nucleotide level from each other and were closest (78.5–97% nucleotide identity) to MT321514-PBoV/GDZ1/China-2018 from China. Overall, there are at least two genetically distinct PBoVs circulating within the colon of individual Pig 46 and multiple genetically distinct PBoVs circulating on the farm.

3.6. Porcine Parvovirus 2 (PPV2)

In addition to the porcine bocaviruses described above, other single stranded DNA parvoviruses identified were PPV2, PPV7, PBuV, and AAV. Sequences corresponding to PPV2 comprised 0.0001 to 0.0279% of the total sequenced reads from both Pig 45 and 46 lung samples. Two partial consensus sequences representing different regions of PPV2 were obtained from Pig 45 lung sample whereas one partial sequence representing VP1 region of PPV2 was obtained from Pig 46 lung sample (Supplementary Table S7).
A 1494 nt long NS1 sequence (MZ544051-PPV2-NS1-AUS-2018) from the Pig 45 lung sample was used for pairwise comparison and phylogenetic analysis with 17 other reference sequences from NCBI GenBank. Our NS1 sequence was found to be closest (~95.4% nucleotide identity) to reference MG345014-PPV2/S4/China-2017 from China but made its own separate branch as shown in supplementary Figure S2. From the Pig 45 lung sample, a 2982 nt long partial VP1 sequence (MZ544050-PPV2-VP1-AUS-2018) was obtained and codes for a 993 amino acid long sequence. On pairwise comparison, our VP1 sequence was closest (~93.7% identical at nt level) to reference KC687097-PPV2/64-PL/Poland-2011 from Poland. Moreover, our partial VP1 sequence (PL46-BC-25-PPV2-3622-3913-AUS-2018) from Pig 46 lung sample was 100% identical with our VP1 sequence from the Pig 45 lung sample. Like our NS1 sequence, our VP1 sequence also formed its own separate branch in the phylogenetic analysis (Supplementary Figure S3), indicating a more distant relationship with other known PPV2 references and also the possibility of a distinct genotype.

3.7. Porcine Circovirus Type 2 (PCV2)

Although in very low abundance (only 0.00003 to 0.0007% of total sequenced reads), we were able to assemble four partial consensus sequences of porcine circovirus type 2 (PCV2) from the Pig 46 colon sample. The number of reads was less than for the other viruses described above, but we were interested in this virus because of its recognized ability to interact with other pathogens and its possible role in immunosuppression, therefore potentially making pigs susceptible to other infections [85]. The consensus sequences obtained included 627 nt (coverage 1–6) and 161 nt (coverage 1–2) long segments of the ORF1 region and 451 nt (coverage 1–3) and 100 nt (coverage 1–4) long segments of the ORF2 region. These were assembled from the mapped reads using a mapping quality threshold of 64 and above and have been deposited in the NCBI Sequence Read Archive (SRA) under SRA accession PRJNA745489.
The ORF2 region of the PCV2 genome is the most variable region and is also the most suitable region for performing phylogenetic studies [86]. The 451 nt long partial ORF2 sequence (PC46-BC-26-PCV2-ORF2-AUS-2018) from the Pig 46 colon sample was used for pairwise sequence comparison and phylogenetic analysis along with 18 other reference sequences representing different genotypes of PCV2 chosen from the NCBI GenBank. The 451 nucleotide long ORF2 sequence from this study fell into the PCV2b genotype cluster with only 1 nucleotide different from the nearest sequences from USA, Canada, and Denmark as shown in Supplementary Figure S4. The 100 nucleotide long ORF2 partial sequence (coverage 1-4) also paired with 100% identity with the nearest sequences from USA, Canada, and Denmark (an overall identity of 99.82% if the 451 nt and 100nt ORF2 sequences are considered together) and also fell within the PCV2b genotype cluster indicating our two partial ORF2 sequences may have come from the same virus.
The 627 nt long partial ORF1 sequence (PC46-BC-26-PCV2-ORF1-AUS-2018) was used for pairwise sequence comparison and phylogenetic analysis with 18 other reference sequences as mentioned above. This analysis suggested that the partial ORF1 section of our PCV2 sequence from Pig 46 colon swab fell into the PCV2a genotype with the nearest reference sequences being MK504415-PCV2/KSU-KS-2017-PCV2-38/USA-2017 and KX828217-PCV2/KU-1614/KOR-2016 with 100% nucleotide identity as shown in supplementary Figure S5. In addition, the 161 nt long (coverage 1–2) ORF1 partial sequence also shared 100% identity with the two same reference sequences from USA and Korea, indicating that our two partial ORF1 sequences may have come from the same virus. The finding that the ORF1 region from the Pig 46 colon is of PCV2a genotype while the ORF2 region is of PCV2b genotype may indicate that two different PCV2 genotypes may be circulating together in the colon sample of Pig 46 or alternatively that it is a single PCV2a/b recombinant PCV2 [87].

3.8. Porcine Torque Teno Sus Virus (TTSuV)

Sequences representing two species of genus Kappatorquevirus, namely TTSuVk2a (0.0004–0.0119% abundance of total sequenced NGS reads) and TTSuVk2b (0.0058-0.0249% abundance of total sequenced NGS reads) were obtained from both pigs (Table 1). Nine partial consensus sequences (117–1276 nt long) of TTSuVk2a representing different ORFs (Supplementary Table S8) were assembled from the mapped reads using the mapping quality threshold of 90 from Pig 45 colon and lung and Pig 46 lung samples. One complete (2565 nt long) coding sequence (MZ544061-TTSuVk2b-AUS-2018) of TTSuVk2b was obtained from the Pig 45 lung sample. Four additional partial consensus sequences of TTSuVk2b of 203–1177 nt long representing different ORFs from Pig 45 colon and Pig 46 lung swabs were also obtained (Supplementary Table S9).
After pairwise comparison of our two partial TTSuVk2a sequences, MZ544058-TTSuVk2a-AUS-2018 and MZ544057-TTSuVk2a-AUS-2018, from the Pig 45 lung sample with the reference sequences in the NCBI GenBank, MN272073-TTSuVk2a/TTV_LV06/Brazil-2016 from Brazil was found to be one of the closest sequences with 97% (840/866) and 95.22% (1215/1276) identity at nt level, respectively. All the other partial sequences of TTSuVk2a obtained from Pig 45 colon and Pig 46 lung samples were identical with the TTSuVk2a sequences from the Pig 45 lung sample.
Pairwise comparison of the 2565 nt long TTSuVk2b sequence (MZ544061-TTSuVk2b-AUS-2018) from the Pig 45 lung sample with reference sequences in the NCBI GenBank, identified MT671970-TTSuVk2b/pork147/Brazil-2018 from Brazil as the closest sequence with 99.92% (2551/2553) identity at nucleotide level with only two nucleotides different (A to C at 64, A to G at 923). The nucleotide positions mentioned here were given on the basis of our TTSuVk2b sequence MZ544061-TTSuVk2b-AUS-2018. In addition, our TTSuVk2b sequence had a 12 nucleotides (TAAAAAGACAAC) insertion at our nucleotide position 2517–2528, which codes for an additional four amino acids (QIKR) towards the 3′end of our TTSuVk2b sequence as compared to the sequence from Brazil (MT671970-TTSuVk2b/pork147/Brazil-2018). Other sequences also shared high nucleotide identity with our sequences, with a 1384 nt long partial TTSuVk2b sequence (MN272088) from Brazil in 2016 being 99.78% identical with our TTSuVk2b sequence, while another 707 nt long partial TTSuVk2b sequence (MK378015) from China in 2017 was 99.86% identical with our TTSuVk2b sequence from the lung of Pig 45. All the partial sequences of TTSuVk2b obtained from Pig 45 colon and Pig 46 lung samples were identical with the TTSuVk2b sequence from the lung of Pig 45. We also observed the same 12 nt insertion in TTSuVk2b sequence from the Pig 45 colon sample whereas we did not get any TTSuk2b sequence from the lung sample from Pig 46 from that region. The 12 nt insertion has also been observed previously by other investigators in other TTSuVK2b sequences in NCBI GenBank. One of the closest sequences with the 12 nt insertion was JQ406844-TTSuVk2b/38E05/NZ-2007 from New Zealand, which was 95.32% identical at nt level with our TTSuVk2b sequence. Although the findings of TTSuV1 sequences have been reported in the past from Australia [88], we did not find any previous sequences reported from Australia representing TTSuVk2a and 2b. To the best of our knowledge, our TTSuVk2b sequence is the first report of this genotype from pigs in Australia. Overall, a near identical TTSuVk2b sequence to that found in the Pig 45 lung sample was found in Brazil in 2016 and again in 2018 and was also found in China in 2017.

3.9. Analysis of Reads Mapping to Bacteria and to Antimicrobial Resistance (AMR) Genes

Among the detected bacterial species of interest, Lawsonia intracellularis was the most abundant and constituted around 2% of all the reads mapped to bacterial nucleic acid in the colon and lung of Pig 45 (2.05–2.08%) and 7.2% of the bacterial reads in the colon of Pig 46. Similarly, Brachyspira intermedia constituted around 0.12 to 0.59% of total bacterial reads and Brachyspira hyodysenteriae constituted around 0.05 to 0.11% of total bacterial reads with reads matching nearly 100% identical with the chromosome or specific plasmids of these bacteria. Similarly, very few reads mapping to the corresponding bacterial ribosomal RNA were found in these pig samples. Moreover, DNA reads (0.06–0.5% of total bacterial reads) representing Campylobacter spp. were also obtained (Supplementary Table S1). The high abundance of reads mapped to Lawsonia intracellularis was consistent with the histopathology findings reported independently by a commercial diagnostic laboratory with intracellular organisms positive on a Warthin–Starry stain and tissue changes consistent with Lawsonia intracellularis enteritis observed by a veterinary pathologist. Bacterial culture isolated a non-haemolytic Brachyspira spp. and was unable to isolate any Salmonella spp. This was also consistent with what we observed in the bacterial read mapping, identifying reads belonging to Brachyspira spp. but not to Salmonella spp.
A low number of NGS reads could be assembled into short contigs representing AMR genes, specifically the tetracycline resistance genes tet (W) and tet (Q) at 97–99% identity over 426–1200 nucleotides and the transposon-mediated lincosamide nucleotidyltransferase, lun (C) gene at 99.5% identity over 388 nucleotides from both pigs were detected except Pig 46 lung sample in which AMR genes were not detected. The low abundance of reads and very short contigs assembled of these AMR genes unfortunately did not allow us to assess which bacteria carried these AMR genes, only that these genes were present in one or more species of the microbial community associated with the pigs. Interestingly, the pigs had been treated with drugs within both of these classes of antibiotics in an attempt to control the clinical signs. We found a relatively high abundance of NGS reads for bacterial chromosomal or plasmid DNA for both Lawsonia intracellularis and Brachyspira spp. compared to their respective ribosomal RNAs, which may indicate that these bacteria may have been dead or inactive possibly as a result of the antibiotic treatments given, and therefore were most likely not the species harbouring the resistance genes.

4. Discussion

Porcine enteric diseases like diarrhoea or other intestinal problems and poor growth performance are some of the major problems in pig herds causing huge economic losses around the world [1,2,3,4,5,6,7,8]. In addition, control of enteric disease in pigs is a major reason why antimicrobials are used on farms. Like in our previous study in dogs [44], here NGS was found to be an effective non-targeted technique for the detection and characterisation of a diverse group of viruses and bacteria in pig colon and lung samples. The method was able to identify and quantify a relatively high abundance of the bacterium Lawsonia intracellularis, which was initially suspected by the veterinarian and subsequently identified by the histopathologist and could be consistent with the clinical disease, poor growth, and lesions observed [25,26,27]. Furthermore, it could also detect the presence of several species of Brachyspira and Campylobacter, which may have also been involved in the diarrhoea [28,29,30,31] and also demonstrated that Salmonella was not present in the samples similar to the microbiological testing. Moreover, the NGS method was able to identify several viruses that may have also played a role in the clinical presentation, and enabled the characterisation of these viruses, many for which this study represents the first time they have been genetically characterised in Australian pigs. The method was also able to detect several viruses in very low abundance (e.g., several of the sapelovirus and enterovirus sequences), which suggested that these viruses were not actively replicating at the time of sample collection. These viruses may have infected the pigs at a younger age, and as they were non-enveloped, were likely stable and present in the pigs’ environment.
The pig samples had 0.07–1.37% of the total sequenced reads mapping to different viruses which individually or together (and in combination with the bacteria detected) may have been involved in causing clinical disease in the pigs. Porcine picornaviruses like PSV, PEV-G, and PTV have been isolated from healthy pigs from different countries [50,67,89,90,91,92,93], but were also found to be associated with diseases like diarrhoea, polioencephalomyelitis, respiratory distress, and pneumonia [11,12,15,49,94,95,96,97]. Similarly, PAstVs have been found in both clinically healthy and diarrheic pigs [98,99]. Experiments performed in seven-day-old piglets to study pathogenesis of PAstV have shown clinical manifestations of mild diarrhoea, growth retardation, and damage of the villi of the small intestinal mucosa of piglets after PAstV infection [100]. In our study, the picornaviruses and astroviruses detected showed a high level of genetic diversity when compared with other such viruses identified and sequenced internationally, which may be the result of a high evolution rate together with selection and possibly recombination [101,102,103]. Moreover, multiple genotypes of picornaviruses were detected, particularly for EV-G and PTV, and at least two genetically distinct PoAstV4 in the Pig 46 colon sample were also detected, showing the possibility of future recombination between different EV-Gs genotypes [104], between different PTVs [105], and also between the different PAstV4 genotypes [48], which may give rise to genetically diverse respective virus genotype.
Porcine single stranded DNA viruses like PBoV, PCV2, and TTSuV have also been found in clinically healthy pigs [106,107,108,109]. However, in a study from Germany, PBoV was also found to be associated with coughing, growth retardation, diarrhoea, and encephalomyelitis [110]. Experimental inoculation of six-week-old healthy pigs with PCV2 has been associated with PCV2- systemic disease (PCV2-SD) displaying weight loss [111]. PCV2 in addition to PBoV, PPV, TTSuV, and other viruses have been found to be associated with the development of PMWS, PDNS, and PRRS, causing huge economic losses in the pork industry worldwide [20,40,107,112,113,114]. Interestingly, in this study PCV2 was found, although in low abundance. It is possible that PCV2 infection played some role in the observed diarrhoea and poor growth in conjunction with the other pathogens. It has been shown that pigs affected with PCV2-SD develop severe lymphocyte depletion [115] and hence immunosuppression leading to secondary infection by other pathogens [116]. In this study, lymph nodes were not collected during the initial necropsy to look specifically for lymphoid depletion, however based on our communication with the veterinarian, this farm initiated PCV2 vaccination after this disease event, and subsequent batches of pigs have not suffered from the same performance issues. Similar responses have been observed by other investigators [117,118], and the response perhaps provides some empirical evidence to suggest that PCV2 may have played some role in the clinical disease in this herd before vaccination. Like other studies, parvoviruses in our study showed an intermediate level of genetic diversity [119,120], however the PBoV3 sequences obtained from Pig 46 colon was highly different from those overseas, particularly the VP1 capsid region which may give rise to genetically distinct virus genotypes possibly by recombination events as shown by other studies [77,80].
While most of the virus sequences identified here indicated some divergence from similar viruses overseas, most likely due to geographical isolation of the Australian pig herd and the policy of no imports of live pigs or semen over the last 40 years or so [121,122], surprisingly, we found in the lung sample of Pig 45 a near identical (99.92% and 99.78% identity at nt level) TTSuVk2b sequence to sequences MT671970 and MN272088 reported from Brazil in 2018 and 2016, and 99.86% identical to MK378015 from China in 2017. Similarly, a near identical (99.82% identical at nt level) PCV2b sequence with sequences (DQ220730; Canada-2005, DQ629118; USA-2006, EF565357; Denmark-2007, FJ233907; Canada-2007 and MK504403; USA-2018) reported from Europe and America was obtained in the colon sample from Pig 46, though some studies have suggested a relatively high evolutionary rate among single stranded DNA (ssDNA) viruses like PCV2 [123,124]. Our TTSuVk2b and PCV2b virus sequences were analysed in-depth to check if the obtained virus sequences were a result of laboratory reagent contaminants through, e.g., spin column as shown for another virus [125]. No such TTSuVk2b and PCV2b reads were obtained from human and bird samples processed on the same NGS chip as our pig samples. Moreover, we did not detect any reads of these virus sequences from a nuclease free water negative control that was processed using the same protocol. Interestingly, by retrospectively looking at some NGS results from pig slaughterhouse pooled materials collected by our laboratory from South Eastern Australia in 2017, we were also able to detect the presence of a similar but not identical TTSuVK2b virus, albeit at low levels in the sample. Interestingly, this slaughterhouse sample had a relatively high abundance of PCV2, albeit different to the PCV2 sequenced in this study. The PCV2 (MZ328311-PCV2-ORF1-AUS-2017 and MZ328312-PCV2-ORF2-AUS-2017) obtained from pig samples from the slaughterhouse was 99.8% identical to a PCV2 sequence (EU886638) from Australia in 2007 [126] both in the ORF1 and ORF2 region but 99% similar in the ORF1 and only 91.6% similar in the ORF2 region of the PCV2 sequence obtained from pigs in this study. It is therefore highly likely that these viruses truly came from the samples taken from the Australian pigs, and given their high identity to some viruses overseas, they possibly have a different pattern of introductions into Australia compared to the other viruses, which showed more evidence of independent evolution to related viruses sequenced elsewhere in the world. As these circular ssDNA viruses have been shown to have high levels of environmental stability [127,128,129,130], it is possible that these viruses may have been inadvertently spread internationally, including to Australia, through trade of non-animal products such as imported minerals, vitamins, other feed additives, or via transport containers. Such spread of infectious agents has been observed overseas, i.e., the spread of porcine epidemic diarrhoea virus (PEDV) in North America in 2013 [131,132].
Taken together, the association of these several bacteria and viruses with disease pathogenesis as shown by various studies fits well with the clinical signs of disease observed in the herd in this study and could potentially support the role of these identified pathogens for disease pathogenesis [26,48,94,133,134] in which the viral infection might have followed the bacterial infection or vice versa. Despite being a relatively expensive and time intensive technology, in particular for data analysis of even a single sample, NGS has become a powerful tool and has demonstrated numerous advantages over conventional targeted detection technologies by detecting co-infections, novel or unique and emerging pathogens for which detection methods do not exist. The sensitivity of the technique is very high, and it was able to detect some virus sequences at very low abundance, which indicated that some viruses were likely simply present in the environment of the pigs rather than actively replicating at the time of sampling. Nevertheless, the detection of these low abundance viruses is still useful as it provides information as to potential infections occurring in the pigs earlier in their lives. Overall, we do show here that metagenomic surveillance has the potential to be a very useful monitoring tool to characterize a wide range of pathogens associated with pig production, tease apart pathogen interactions in multifactorial disease, and has the potential to be used as a tool to look for unexpected introductions and transmission of viruses and other micro-organisms between pig farms to assist producers in identifying gaps in biosecurity to prevent the arrival of other pathogens. Further genetic sequencing of more sequences of these viruses from Australian farms will better establish their evolutionary rate and better identify the time frame of their arrival into the Australian pig herd. Moreover, further studies looking at the virus and bacterial pathogens in both healthy pigs and those with poor growth and diarrhoea are needed not only to genetically characterise these diverse potential pathogens but also to clearly understand the detailed pathological and epidemiological roles of the individual pathogens in clinical disease in pigs.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13081608/s1, Table S1. Table showing the abundance of bacterial chromosomal DNA/plasmid reads in different pig swab samples. Table S2. Ten consensus partial sequences of porcine sapelovirus (PSV) were obtained from Pig 45 and Pig 46 colon and lung swab samples. The values for each sequence length, IGV coverage, mapping quality and accession number are shown. Table S3. Twenty two consensus partial sequences of porcine enterovirus-G (PEV-G) were obtained from Pig 45 and Pig 46 colon and lung swab samples. They were compared with one of the longest obtained sequence MZ515509-EV-G-AUS-2018 of PEV-G and the positions were determined for these partial sequences. Table S4. Four and seven partial consensus sequences of porcine teschovirus (PTVs) were obtained from Pig 45 and Pig 46 colon swab samples, respectively. They were compared with one of the reference sequence AF296088-PTV3/O 2b/GER-2000 of PTV and the positions were determined for these partial sequences. Table S5. Ten partial consensus sequences of porcine astrovirus (PAstV) were obtained from Pig 46 colon swab sample. They were compared with one of the reference sequence MT470220-PoAstV-CX1-Anhui-China-2018 of PAstV and the positions were determined for these partial sequences. Table S6. Twenty two partial consensus sequences of porcine bocavirus (PBoV) were obtained from Pig 45 and Pig 46 colon and lung swab samples. They were compared with one of the reference sequence KF025384-PBoV/MN154-1/USA-2011 of PBoV and the positions were determined for these partial sequences. Table S7: Three partial consensus sequences of porcine parvovirus 2 (PPV2) were obtained from Pig 45 and Pig 46 lung swab samples. They were compared with one of the reference sequence MG345013-PPV2/S1/China-2017 of PPV2 and the positions were determined for these partial sequences. Table S8. Nine partial consensus sequences of torque teno sus virus k2a (TTSuVk2a) were obtained from Pig 45 colon and lung and Pig 46 lung swab samples. They were compared with one of the reference sequence MN272073-TTSuVk2a/TTV LV06/Brazil-2016 of TTSuVk2a and the positions were determined for these partial sequences. Table S9. One complete coding sequence (cds) torque teno sus virus k2b TTSuVk2b was obtained from Pig 45 lung swab and four partial consensus sequences of TTSuVk2b were obtained from Pig 45 colon and Pig 46 lung swab samples. The longest obtained sequence MZ544061-TTSuVk2b-AUS-2018was used to determine the given positions for the other short partial sequences. Figure S1. Phylogenetic analysis of near complete nucleotide sequence of VP4-VP3 capsid coding region of EV-G. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the General Time Reversible (GTR + G + I) [70] model with a bootstrapping of 1000 replicates. The analysis involved 22 reference sequences of the VP4- VP3 capsid region of EV-Gs genome and six EV-G sequences from this study. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. Sequence names indicated in red colour from Pig 45 colon and lung samples have been labelled with black (⬤) and white (⭘) circle, respectively while Pig 46 colon samples have been labelled with a black rhombus (♦). Figure S2. Phylogenetic analysis of nucleotide sequence of partial NS1 region of PPV2. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the (T92 + G + I) [135] model with a bootstrapping of 1000 replicates. The analysis involved 17 reference sequences of the NS1 non-structural region of PPV2 genome and one PPV2 sequence from this study. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. Sequence name indicated in red colour from Pig 45 lung sample has been labelled with white circle (⭘). Figure S3. Phylogenetic analysis of nucleotide sequence of partial VP1 region of PPV2. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the (T92 + G + I) [135] model with a bootstrapping of 1000 replicates. The analysis involved 21 reference sequences of the VP1 capsid region of PPV2 genome and one PPV2 sequence from this study. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. Sequence name indicated in red colour from Pig 45 lung sample has been labelled with white circle (⭘). Figure S4: Phylogenetic analysis of partial nucleotide sequence of ORF2 region of PCV2. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the General Time Reversible (GTR + G + I) [70] model with a bootstrapping of 1000 replicates. The analysis involved 17 reference sequences of the ORF2 capsid region of PCV2 genome and one PCV2 sequence from this study. Additionally, one other PCV2 sequence from our previous study is included. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. Sequence named indicated in red colour are from Australia. The sequence from this study (i.e., Pig 46 colon) has been labelled with a black rhombus (♦) whereas sequence from our other study (not published, but available at NCBI Accession number MZ328312) has been labelled with a black triangle (▲). Figure S5: Phylogenetic analysis of partial nucleotide sequence of ORF1 region of PCV2. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the General Time Reversible (GTR + G + I) [70] model with a bootstrapping of 1000 replicates. The analysis involved 17 reference sequences of the ORF1 replicase region of PCV2 genome and one PCV2 sequence from this study. One other PCV2 sequence from our other study was also included. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. Sequence name indicated in red colour are from Australia. The sequence from this study (i.e., Pig 46 colon) has been labelled with a black rhombus (♦) whereas a sequence from our previous study (not published, but available at NCBI Accession number MZ328311) has been labelled with a black triangle (▲).

Author Contributions

Conceptualization, S.A.; methodology, T.R.B.; software, T.R.B. and S.A.; validation, T.R.B., S.A. and A.C.; formal analysis, T.R.B. and S.A.; investigation, T.R.B., S.A. and A.C.; resources, S.A. and A.C.; data curation, T.R.B., S.A. and A.C.; writing—original draft preparation, T.R.B. and S.A.; writing—review and editing, T.R.B., S.A. and A.C.; visualization, T.R.B., S.A. and A.C.; supervision, S.A.; funding acquisition, S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Australian Government Research Training Program Scholarship through Deakin University; NHMRC Equipment Grant GNT9000413 through Deakin University; funding provided by Deakin University, Barwon Health, and CSIRO to Soren Alexandersen.

Institutional Review Board Statement

Ethical review and approval were waived for this study, due to the samples used in this investigation were surplus colon and lung tissue samples collected by the attending farm veterinarian as part of their routine necropsy investigations to determine the cause of disease in the pigs, and surplus left over after submission for diagnostic laboratory testing. The studies described here were performed in accordance with all relevant guidelines and regulations.

Data Availability Statement

With accession numbers (MZ515498–MZ515550, MZ544028–MZ544063). The sequence reads of our porcine circovirus 2 (PCV2) have been deposited in the NCBI Sequence Read Archive (SRA) under SRA accession: PRJNA745489. The fasta file used for PCV2 mapping is available in Supplementary Data 1. Additional datasets analysed in the paper can be made available from the authors upon reasonable request.

Acknowledgments

We want to acknowledge Patrick Daniel for providing the surplus pig colon and lung samples and for his helpful comments on the draft manuscript. We also acknowledge Trish Holyoake, Jason Hodge, and Jessy Vibin for their helpful comments on the draft manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Arnold, M.; Crienen, A.; Swam, H.; Berg, S.V.; Jolie, R.; Nathues, H. Correlation of Lawsonia intracellularis positivity in quantitative PCR and herd factors in European pig herds. Porc. Health Manag. 2021, 7, 1–8. [Google Scholar] [CrossRef]
  2. Cornelison, A.; Karriker, L.; Williams, N.; Haberl, B.; Stalder, K.; Schulz, L.; Patience, J. Impact of health challenges on pig growth performance, carcass characteristics, and net returns under commercial conditions. Transl. Anim. Sci. 2018, 2, 50–61. [Google Scholar] [CrossRef] [Green Version]
  3. Sjolund, M.; Zoric, M.; Wallgren, P. Financial impact on pig production III: Gastrointestinal disorders. In Proceedings of the 6th European Symposium of Porcine Health Management, Sorrento, Italy, 7–9 May 2014. [Google Scholar]
  4. Dufresne, L. Economics of Pig Health Improvements. Available online: https://porkgateway.org/wp-content/uploads/2015/07/economics-of-pig-health-improvement1.pdf (accessed on 10 June 2021).
  5. Van Breda, L.K.; Dhungyel, O.P.; Ginn, A.N.; Iredell, J.R.; Ward, M.P. Pre-and post-weaning scours in southeastern Australia: A survey of 22 commercial pig herds and characterisation of Escherichia coli isolates. PLoS ONE 2017, 12, e0172528. [Google Scholar] [CrossRef] [Green Version]
  6. Schulz, L.L.; Tonsor, G.T. Assessment of the economic impacts of porcine epidemic diarrhoea virus in the United States. J. Anim. Sci. 2015, 93, 5111–5118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Robbins, R.; Almond, G.; Byers, E. Swine Diseases and Disorders. Encycl. Agric. Food Syst. 2014, 261–276. [Google Scholar] [CrossRef]
  8. VanderWaal, K.; Deen, J. Global trends in infectious diseases of swine. Proc. Natl. Acad. Sci. USA 2018, 115, 11495–11500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Ramirez, A. Diseases affecting pigs: An overview of common bacterial, viral and parasitic pathogens of pigs. In Achieving Sustainable Production of Pig Meat: Animal Health and Welfare; Alejandro, R., Ed.; Iowa State University: Ames, IA, USA, 2018; Volume 3, pp. 3–29. [Google Scholar]
  10. Bergeland, M.E.; Henry, S.C. Infectious diarrhoeas of young pigs. Vet. Clin. N. Am. Large Anim. Pract. 1982, 4, 389–399. [Google Scholar]
  11. Matias Ferreyra, F.; Arruda, B.; Stevenson, G.; Schwartz, K.; Madson, D.; Yoon, K.-J.; Zhang, J.; Piñeyro, P.; Chen, Q.; Arruda, P. Development of polioencephalomyelitis in cesarean-derived colostrum-deprived pigs following experimental inoculation with either teschovirus a serotype 2 or serotype 11. Viruses 2017, 9, 179. [Google Scholar] [CrossRef] [Green Version]
  12. Arruda, P.; Arruda, B.; Schwartz, K.; Vannucci, F.; Resende, T.; Rovira, A.; Sundberg, P.; Nietfeld, J.; Hause, B. Detection of a novel sapelovirus in central nervous tissue of pigs with polioencephalomyelitis in the USA. Transbound. Emerg. Dis. 2017, 64, 311–315. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, S.; Wang, Y.; Shen, Q.; Zhang, W.; Hua, X. Prevalence of porcine enterovirus 9 in pigs in middle and eastern China. Virol. J. 2013, 10, 1–6. [Google Scholar] [CrossRef] [Green Version]
  14. Kim, D.-S.; Kang, M.-I.; Son, K.-Y.; Bak, G.-Y.; Park, J.-G.; Hosmillo, M.; Seo, J.-Y.; Kim, J.-Y.; Alfajaro, M.M.; Soliman, M. Pathogenesis of Korean Sapelovirus A in piglets and chicks. J. Gen. Virol. 2016, 97, 2566–2574. [Google Scholar] [CrossRef] [PubMed]
  15. Schock, A.; Gurrala, R.; Fuller, H.; Foyle, L.; Dauber, M.; Martelli, F.; Scholes, S.; Roberts, L.; Steinbach, F.; Dastjerdi, A. Investigation into an outbreak of encephalomyelitis caused by a neuroinvasive porcine sapelovirus in the United Kingdom. Vet. Microbiol. 2014, 172, 381–389. [Google Scholar] [CrossRef] [PubMed]
  16. Padmanabhan, A.; Hause, B.M. Detection and characterization of a novel genotype of porcine astrovirus 4 from nasal swabs from pigs with acute respiratory disease. Arch. Virol. 2016, 161, 2575–2579. [Google Scholar] [CrossRef]
  17. Boros, Á.; Albert, M.; Pankovics, P.; Bíró, H.; Pesavento, P.A.; Phan, T.G.; Delwart, E.; Reuter, G. Outbreaks of neuroinvasive astrovirus associated with encephalomyelitis, weakness, and paralysis among weaned pigs, Hungary. Emerg. Infect. Dis. 2017, 23, 1982–1993. [Google Scholar] [CrossRef]
  18. Rawal, G.; Ferreyra, F.M.; Macedo, N.R.; Bradner, L.K.; Harmon, K.M.; Allison, G.; Linhares, D.C.; Arruda, B.L. Ecology of Porcine Astrovirus Type 3 in a Herd with Associated Neurologic Disease. Viruses 2020, 12, 992. [Google Scholar] [CrossRef]
  19. Batista, L. Porcine Circovirus Associated Diseases (PCVAD) in Canada-Prevalence, Co-Factors, and Risk Factors. In Advances in Pork Production, Proceedings of the 2007 Banff Pork Seminar, Banff, Alberta, Canada, 16–19 January 2007; University of Alberta: Edmonton, AB, Canada, 2007; Volume 18, pp. 57–68. [Google Scholar]
  20. Segalés, J. Porcine circovirus type 2 (PCV2) infections: Clinical signs, pathology and laboratory diagnosis. Virus Res. 2012, 164, 10–19. [Google Scholar] [CrossRef] [PubMed]
  21. Phaneuf, L.R.; Ceccarelli, A.; Laing, J.R.; Moloo, B.; Turner, P.V. Porcine dermatitis and nephropathy syndrome associated with porcine circovirus 2 infection in a Yorkshire pig. J. Am. Assoc. Lab. Anim. Sci. 2007, 46, 68–72. [Google Scholar] [PubMed]
  22. Allan, G.M.; Ellis, J.A. Porcine circoviruses: A review. J. Vet. Diagn. Investig. 2000, 12, 3–14. [Google Scholar] [CrossRef] [PubMed]
  23. Qin, S.; Ruan, W.; Yue, H.; Tang, C.; Zhou, K.; Zhang, B. Viral communities associated with porcine respiratory disease complex in intensive commercial farms in Sichuan province, China. Sci. Rep. 2018, 8, 13341. [Google Scholar] [CrossRef]
  24. McMenamy, M.J.; McKillen, J.; McNair, I.; Duffy, C.; Blomström, A.-L.; Charreyre, C.; Welsh, M.; Allan, G. 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. Vet. Microbiol. 2013, 164, 293–298. [Google Scholar] [PubMed]
  25. McOrist, S.; Smith, S.; Green, L. Estimate of direct financial losses due to porcine proliferative enteropathy. Vet. Rec. 1997, 140, 579–581. [Google Scholar] [CrossRef] [PubMed]
  26. Karuppannan, A.K.; Opriessnig, T. Lawsonia intracellularis: Revisiting the disease ecology and control of this fastidious pathogen in pigs. Front. Vet. Sci. 2018, 5, 181. [Google Scholar] [CrossRef]
  27. Vannucci, F.; Gebhart, C. Recent advances in understanding the pathogenesis of Lawsonia intracellularis infections. Vet. Pathol. 2014, 51, 465–477. [Google Scholar] [CrossRef] [PubMed]
  28. Osorio, J.; Carvajal, A.; Naharro, G.; Rubio, P.; La, T.; Phillips, N.; Hampson, D. Identification of weakly haemolytic Brachyspira isolates recovered from pigs with diarrhoea in Spain and Portugal and comparison with results from other countries. Res. Vet. Sci. 2013, 95, 861–869. [Google Scholar] [CrossRef] [Green Version]
  29. Clothier, K.A.; Kinyon, J.M.; Frana, T.S.; Naberhaus, N.; Bower, L.; Strait, E.L.; Schwartz, K. Species characterization and minimum inhibitory concentration patterns of Brachyspira species isolates from swine with clinical disease. J. Vet. Diagn. Investig. 2011, 23, 1140–1145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Hampson, D.J.; La, T.; Phillips, N.D. Emergence of Brachyspira species and strains: Reinforcing the need for surveillance. Porc. Health Manag. 2015, 1, 1–6. [Google Scholar] [CrossRef] [Green Version]
  31. Burrough, E.; Terhorst, S.; Sahin, O.; Zhang, Q. Prevalence of Campylobacter spp. relative to other enteric pathogens in grow-finish pigs with diarrhoea. Anaerobe 2013, 22, 111–114. [Google Scholar] [CrossRef]
  32. Edfors-Lilja, I.; Wallgren, P.; Axford, R.; Bishop, S.; Nicholas, F.; Owen, J. Escherichia coli and Salmonella diarrhoea in pigs. In Breedings for Disease Resistance in Farm Animals, 2nd ed.; CABI Publisher: Wallingford, UK, 2000; Volume 12, pp. 253–267. [Google Scholar]
  33. Suh, D.K.; Song, J.C. Simultaneous detection of Lawsonia intracellularis, Brachyspira hyodysenteriae and Salmonella spp. in swine intestinal specimens by multiplex polymerase chain reaction. J. Vet. Sci. 2005, 6, 231–237. [Google Scholar] [CrossRef]
  34. Leblanc-Maridor, M.; Beaudeau, F.; Seegers, H.; Denis, M.; Belloc, C. Rapid identification and quantification of Campylobacter coli and Campylobacter jejuni by real-time PCR in pure cultures and in complex samples. BMC Microbiol. 2011, 11, 1–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Bhatta, T.R.; Ryt-Hansen, P.; Nielsen, J.P.; Larsen, L.E.; Larsen, I.; Chamings, A.; Goecke, N.B.; Alexandersen, S. Infection Dynamics of Swine Influenza Virus in a Danish Pig Herd Reveals Recurrent Infections with Different Variants of the H1N2 Swine Influenza A Virus Subtype. Viruses 2020, 12, 1013. [Google Scholar] [CrossRef] [PubMed]
  36. Chelli, E.; De Sabato, L.; Vaccari, G.; Ostanello, F.; Di Bartolo, I. Detection and Characterization of Porcine Sapelovirus in Italian Pig Farms. Animals 2020, 10, 966. [Google Scholar] [CrossRef] [PubMed]
  37. Palmquist, J.M.; Munir, S.; Taku, A.; Kapur, V.; Goyal, S.M. Detection of porcine teschovirus and enterovirus type II by reverse transcription–polymerase chain reaction. J. Vet. Diagn. Investig. 2002, 14, 476–480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Chen, M.; Tang, W.; Hua, X. Molecular characterization of a porcine teschovirus HuN-1 isolate proliferating in PK-15 cell. BMC Vet. Res. 2018, 14, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Mor, S.K.; Chander, Y.; Marthaler, D.; Patnayak, D.P.; Goyal, S.M. Detection and molecular characterization of Porcine astrovirus strains associated with swine diarrhoea. J. Vet. Diagn. Investig. 2012, 24, 1064–1067. [Google Scholar] [CrossRef] [Green Version]
  40. Liu, M.; Li, Y.; Sun, D.; Xia, Y.; Huang, J.; Guo, L. Detection and genetic analysis of porcine bocavirus in different swine herds in North Central China. Sci. World J. 2014, 2014. [Google Scholar] [CrossRef]
  41. Jacob, D.M.; Lee, C.Y.; Arshad, S.S.; Selvarajah, G.T.; Bande, F.; Ong, B.L.; Ooi, P.T. First molecular detection of porcine bocavirus in Malaysia. Trop. Anim. Health Prod. 2018, 50, 733–739. [Google Scholar] [CrossRef]
  42. Yang, K.; Jiao, Z.; Zhou, D.; Guo, R.; Duan, Z.; Tian, Y. Development of a multiplex PCR to detect and discriminate porcine circoviruses in clinical specimens. BMC Infect. Dis. 2019, 19, 1–10. [Google Scholar] [CrossRef] [Green Version]
  43. Blomström, A.-L. Viral metagenomics as an emerging and powerful tool in veterinary medicine. Vet. Q. 2011, 31, 107–114. [Google Scholar] [CrossRef]
  44. Bhatta, T.R.; Chamings, A.; Vibin, J.; Alexandersen, S. Detection and characterisation of canine astrovirus, canine parvovirus and canine papillomavirus in puppies using next generation sequencing. Sci. Rep. 2019, 9, 1–10. [Google Scholar]
  45. Vibin, J.; Chamings, A.; Collier, F.; Klaassen, M.; Nelson, T.M.; Alexandersen, S. Metagenomics detection and characterisation of viruses in faecal samples from Australian wild birds. Sci. Rep. 2018, 8, 1–23. [Google Scholar] [CrossRef]
  46. Vibin, J.; Chamings, A.; Klaassen, M.; Bhatta, T.R.; Alexandersen, S. Metagenomic characterisation of avian parvoviruses and picornaviruses from Australian wild ducks. Sci. Rep. 2020, 10, 1–15. [Google Scholar] [CrossRef] [PubMed]
  47. Blomström, A.-L.; Fossum, C.; Wallgren, P.; Berg, M. Viral metagenomic analysis displays the co-infection situation in healthy and PMWS affected pigs. PLoS ONE 2016, 11, e0166863. [Google Scholar] [CrossRef] [PubMed]
  48. Ito, M.; Kuroda, M.; Masuda, T.; Akagami, M.; Haga, K.; Tsuchiaka, S.; Kishimoto, M.; Naoi, Y.; Sano, K.; Omatsu, T. Whole genome analysis of porcine astroviruses detected in Japanese pigs reveals genetic diversity and possible intra-genotypic recombination. Infect. Genet. Evol. 2017, 50, 38–48. [Google Scholar] [CrossRef]
  49. Tsuchiaka, S.; Naoi, Y.; Imai, R.; Masuda, T.; Ito, M.; Akagami, M.; Ouchi, Y.; Ishii, K.; Sakaguchi, S.; Omatsu, T. Genetic diversity and recombination of enterovirus G strains in Japanese pigs: High prevalence of strains carrying a papain-like cysteine protease sequence in the enterovirus G population. PLoS ONE 2018, 13, e0190819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Harima, H.; Kajihara, M.; Simulundu, E.; Bwalya, E.; Qiu, Y.; Isono, M.; Okuya, K.; Gonzalez, G.; Yamagishi, J.; Hang’ombe, B.M. Genetic and biological diversity of porcine Sapeloviruses prevailing in Zambia. Viruses 2020, 12, 180. [Google Scholar] [CrossRef] [Green Version]
  51. Wu, Z.; Ren, X.; Yang, L.; Hu, Y.; Yang, J.; He, G.; Zhang, J.; Dong, J.; Sun, L.; Du, J. Virome analysis for identification of novel mammalian viruses in bat species from Chinese provinces. J. Virol. 2012, 86, 10999–11012. [Google Scholar] [CrossRef] [Green Version]
  52. Li, L.; Pesavento, P.A.; Shan, T.; Leutenegger, C.M.; Wang, C.; Delwart, E. Viruses in diarrhoeic dogs include novel kobuviruses and sapoviruses. J. Gen. Virol. 2011, 92, 2534–2541. [Google Scholar] [CrossRef]
  53. Moreno, P.S.; Wagner, J.; Mansfield, C.S.; Stevens, M.; Gilkerson, J.R.; Kirkwood, C.D. Characterisation of the canine faecal virome in healthy dogs and dogs with acute diarrhoea using shotgun metagenomics. PLoS ONE 2017, 12, e0178433. [Google Scholar] [CrossRef] [Green Version]
  54. Phan, T.G.; Kapusinszky, B.; Wang, C.; Rose, R.K.; Lipton, H.L.; Delwart, E.L. The fecal viral flora of wild rodents. PLoS Pathog. 2011, 7, e1002218. [Google Scholar] [CrossRef] [Green Version]
  55. Victoria, J.G.; Kapoor, A.; Li, L.; Blinkova, O.; Slikas, B.; Wang, C.; Naeem, A.; Zaidi, S.; Delwart, E. Metagenomic analyses of viruses in stool samples from children with acute flaccid paralysis. J. Virol. 2009, 83, 4642–4651. [Google Scholar] [CrossRef] [Green Version]
  56. Yang, J.; Yang, F.; Ren, L.; Xiong, Z.; Wu, Z.; Dong, J.; Sun, L.; Zhang, T.; Hu, Y.; Du, J. Unbiased parallel detection of viral pathogens in clinical samples by use of a metagenomic approach. J. Clin. Microbiol. 2011, 49, 3463–3469. [Google Scholar] [CrossRef] [Green Version]
  57. Bhatta, T.R.; Chamings, A.; Vibin, J.; Klaassen, M.; Alexandersen, S. Detection of a reassortant H9N2 avian influenza virus with intercontinental gene segments in a resident australian chestnut teal. Viruses 2020, 12, 88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  59. Mount, D.W. Using the basic local alignment search tool (BLAST). Cold Spring Harb. Protoc. 2007, 2007, 17. [Google Scholar] [CrossRef] [PubMed]
  60. Caboche, S.; Audebert, C.; Lemoine, Y.; Hot, D. Comparison of mapping algorithms used in high-throughput sequencing: Application to Ion Torrent data. BMC Genom. 2014, 15, 1–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Thorvaldsdóttir, H.; Robinson, J.T.; Mesirov, J.P. Integrative Genomics Viewer (IGV): High-performance genomics data visualization and exploration. Brief. Bioinform. 2013, 14, 178–192. [Google Scholar] [CrossRef] [Green Version]
  62. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [Green Version]
  63. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  64. Martin, D.; Posada, D.; Crandall, K.; Williamson, C. A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res. Hum. Retrovir. 2005, 21, 98–102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Oberste, M.S.; Maher, K.; Kilpatrick, D.R.; Pallansch, M.A. Molecular evolution of the human enteroviruses: Correlation of serotype with VP1 sequence and application to picornavirus classification. J. Virol. 1999, 73, 1941–1948. [Google Scholar] [CrossRef] [Green Version]
  66. Conserved Protein Domain Family DUF3724. Available online: https://www.ncbi.nlm.nih.gov/Structure/cdd/cddsrv.cgi?uid=152955 (accessed on 20 July 2021).
  67. Van Dung, N.; Anh, P.H.; Van Cuong, N.; Hoa, N.T.; Carrique-Mas, J.; Hien, V.B.; Campbell, J.; Baker, S.; Farrar, J.; Woolhouse, M.E. Prevalence, genetic diversity and recombination of species G enteroviruses infecting pigs in Vietnam. J. Gen. Virol. 2014, 95, 549–556. [Google Scholar] [CrossRef]
  68. Vilar, M.; Peralta, B.; García-Bocanegra, I.; Simon-Grifé, M.; Bensaid, A.; Casal, J.; Segalés, J.; Pina-Pedrero, S. Distribution and genetic characterization of Enterovirus G and Sapelovirus A in six Spanish swine herds. Virus Res. 2016, 215, 42–49. [Google Scholar] [CrossRef] [PubMed]
  69. Bunke, J.; Receveur, K.; Oeser, A.C.; Fickenscher, H.; Zell, R.; Krumbholz, A. High genetic diversity of porcine enterovirus G in Schleswig-Holstein, Germany. Arch. Virol. 2018, 163, 489–493. [Google Scholar] [CrossRef]
  70. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University Press: Oxford, UK, 2000. [Google Scholar]
  71. Usherwood, E.; Nash, A. Lymphocyte recognition of picornaviruses. J. Gen. Virol. 1995, 76, 499–508. [Google Scholar] [CrossRef] [PubMed]
  72. Kaku, Y.; Murakami, Y.; Sarai, A.; Wang, Y.; Ohashi, S.; Sakamoto, K. Antigenic properties of porcine teschovirus 1 (PTV-1) Talfan strain and molecular strategy for serotyping of PTVs. Arch. Virol. 2007, 152, 929–940. [Google Scholar] [CrossRef]
  73. Cano-Gómez, C.; Palero, F.; Buitrago, M.D.; García-Casado, M.A.; Fernández-Pinero, J.; Fernández-Pacheco, P.; Agüero, M.; Gómez-Tejedor, C.; Jiménez-Clavero, M.Á. Analyzing the genetic diversity of teschoviruses in Spanish pig populations using complete VP1 sequences. Infect. Genet. Evol. 2011, 11, 2144–2150. [Google Scholar] [CrossRef]
  74. De Benedictis, P.; Schultz-Cherry, S.; Burnham, A.; Cattoli, G. Astrovirus infections in humans and animals–molecular biology, genetic diversity, and interspecies transmissions. Infect. Genet. Evol. 2011, 11, 1529–1544. [Google Scholar] [CrossRef]
  75. Dong, J.; Dong, L.; Méndez, E.; Tao, Y. Crystal structure of the human astrovirus capsid spike. Proc. Natl. Acad. Sci. USA 2011, 108, 12681–12686. [Google Scholar] [CrossRef] [Green Version]
  76. Matsui, M.; Ushijima, H.; Hachiya, M.; Kakizawa, J.; Wen, L.; Oseto, M.; Morooka, K.; Kurtz, J.B. Determination of serotypes of astroviruses by reverse transcription-polymerase chain reaction and homologies of the types by the sequencing of Japanese isolates. Microbiol. Immunol. 1998, 42, 539–547. [Google Scholar] [CrossRef] [Green Version]
  77. Lau, S.K.; Woo, P.C.; Yip, C.C.; Li, K.S.; Fu, C.T.; Huang, Y.; Chan, K.-H.; Yuen, K.-Y. 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. J. Gen. Virol. 2011, 92, 2047–2059. [Google Scholar] [CrossRef] [PubMed]
  78. McKillen, J.; McNeilly, F.; Duffy, C.; McMenamy, M.; McNair, I.; Hjertner, B.; Millar, A.; McKay, K.; Lagan, P.; Adair, B. Isolation in cell cultures and initial characterisation of two novel bocavirus species from swine in Northern Ireland. Vet. Microbiol. 2011, 152, 39–45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Li, B.; Ma, J.; Xiao, S.; Fang, L.; Zeng, S.; Wen, L.; Zhang, X.; Ni, Y.; Guo, R.; Yu, Z. Complete genome sequence of a novel species of porcine bocavirus, PBoV5. Am. Soc. Microbiol. 2012. [Google Scholar] [CrossRef] [Green Version]
  80. Jiang, Y.-H.; Xiao, C.-T.; Yin, S.-H.; Gerber, P.F.; Halbur, P.G.; Opriessnig, T. High prevalence and genetic diversity of porcine bocaviruses in pigs in the USA, and identification of multiple novel porcine bocaviruses. J. Gen. Virol. 2014, 95, 453–465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Xiao, C.-T.; Halbur, P.G.; Opriessnig, T. Molecular evolutionary genetic analysis of emerging parvoviruses identified in pigs. Infect. Genet. Evol. 2013, 16, 369–376. [Google Scholar] [CrossRef] [PubMed]
  82. Yang, W.-Z.; Yu, J.-M.; Li, J.-S.; Cheng, W.-X.; Huang, C.-P.; Duan, Z.-J. Genome characterization of a novel porcine bocavirus. Arch. Virol. 2012, 157, 2125–2132. [Google Scholar] [CrossRef]
  83. Wang, E.; Liu, W.; Yang, B.; Liu, J.; Ma, X.; Lan, X. Complete sequence and phylogenetic analysis of a porcine bocavirus strain swBoV CH437. Virus Genes 2014, 48, 387–390. [Google Scholar] [CrossRef] [PubMed]
  84. Hasegawa, M.; Kishino, H.; Yano, T.-A. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 1985, 22, 160–174. [Google Scholar] [CrossRef]
  85. Yang, N.; Li, J.; Yang, Q.; Qiao, J.; Cui, D.; Liu, F.; Li, H.; Zhou, S. Reduced antigen presentation capability and modified inflammatory/immunosuppressive cytokine expression of induced monocyte-derived dendritic cells from peripheral blood of piglets infected with porcine circovirus type 2. Arch. Virol. 2018, 163, 1231–1239. [Google Scholar] [CrossRef]
  86. Olvera, A.; Cortey, M.; Segales, J. Molecular evolution of porcine circovirus type 2 genomes: Phylogeny and clonality. Virology 2007, 357, 175–185. [Google Scholar] [CrossRef] [Green Version]
  87. Ramos, N.; Mirazo, S.; Castro, G.; Arbiza, J. Molecular analysis of Porcine Circovirus Type 2 strains from Uruguay: Evidence for natural occurring recombination. Infect. Genet. Evol. 2013, 19, 23–31. [Google Scholar] [CrossRef]
  88. Cortey, M.; Pileri, E.; Segalés, J.; Kekarainen, T. Globalisation and global trade influence molecular viral population genetics of Torque Teno Sus Viruses 1 and 2 in pigs. Vet. Microbiol. 2012, 156, 81–87. [Google Scholar] [CrossRef] [PubMed]
  89. Sunaga, F.; Masuda, T.; Ito, M.; Akagami, M.; Naoi, Y.; Sano, K.; Katayama, Y.; Omatsu, T.; Oba, M.; Sakaguchi, S. Complete genomic analysis and molecular characterization of Japanese porcine sapeloviruses. Virus Genes 2019, 55, 198–208. [Google Scholar] [CrossRef] [PubMed]
  90. Prodělalová, J. The survey of porcine teschoviruses, sapeloviruses and enteroviruses B infecting domestic pigs and wild boars in the Czech Republic between 2005 and 2011. Infect. Genet. Evol. 2012, 12, 1447–1451. [Google Scholar] [CrossRef] [PubMed]
  91. Boros, Á.; Pankovics, P.; Reuter, G. Characterization of a novel porcine enterovirus in domestic pig in Hungary. Infect. Genet. Evol. 2011, 11, 1096–1102. [Google Scholar] [CrossRef]
  92. Zhang, W.; Yang, S.; Shen, Q.; Ren, L.; Shan, T.; Wei, J.; Cui, L.; Hua, X. Complete genome sequence of a novel porcine enterovirus strain in China. J. Virol. 2012, 86, 7008–7009. [Google Scholar] [CrossRef] [Green Version]
  93. Moon, H.-J.; Song, D.; Seon, B.H.; Kim, H.-K.; Park, S.-J.; An, D.-J.; Kim, J.-M.; Kang, B.-K.; Park, B.-K. Complete genome analysis of porcine enterovirus B isolated in Korea. J. Virol. 2012, 86, 10250. [Google Scholar] [CrossRef] [Green Version]
  94. Knutson, T.P.; Velayudhan, B.T.; Marthaler, D.G. A porcine enterovirus G associated with enteric disease contains a novel papain-like cysteine protease. J. Gen. Virol. 2017, 98, 1305–1310. [Google Scholar] [CrossRef] [PubMed]
  95. Lee, S.; Lee, C. First detection of novel enterovirus G recombining a torovirus papain–like protease gene associated with diarrhoea in swine in South Korea. Transbound. Emerg. Dis. 2019, 66, 1023–1028. [Google Scholar] [CrossRef] [Green Version]
  96. Malik, Y.S.; Bhat, S.; Vlasova, A.N.; Wang, F.-I.; Touil, N.; Ghosh, S.; Dhama, K.; Yadav, M.P.; Singh, R.K. Emerging and Transboundary Animal Viruses; Springer: Singapore, 2020; pp. 123–136. [Google Scholar]
  97. Forman, A.; Pass, D.; Connaughton, I. The characterisation and pathogenicity of porcine enteroviruses isolated in Victoria. Aust. Vet. J. 1982, 58, 136–142. [Google Scholar] [CrossRef]
  98. Xiao, C.-T.; Gimenez-Lirola, L.G.; Gerber, P.F.; Jiang, Y.-H.; Halbur, P.G.; Opriessnig, T. Identification and characterization of novel porcine astroviruses (PAstVs) with high prevalence and frequent co-infection of individual pigs with multiple PAstV types. J. Gen. Virol. 2013, 94, 570–582. [Google Scholar] [CrossRef]
  99. Lv, S.-L.; Zhang, H.-H.; Li, J.-Y.; Hu, W.-Q.; Song, Y.-T.; Opriessnig, T.; Xiao, C.-T. High genetic diversity and recombination events of porcine astrovirus strains identified from ill and asymptomatic pigs in 2017, Hunan Province, China. Virus Genes 2019, 55, 673–681. [Google Scholar] [CrossRef]
  100. Fang, Q.; Wang, C.; Liu, H.; Wu, Q.; Liang, S.; Cen, M.; Dong, Q.; Wei, Y.; Chen, Y.; Ouyang, K. Pathogenic characteristics of a porcine astrovirus strain isolated in China. Viruses 2019, 11, 1156. [Google Scholar] [CrossRef] [Green Version]
  101. Austermann-Busch, S.; Becher, P. RNA structural elements determine frequency and sites of nonhomologous recombination in an animal plus-strand RNA virus. J. Virol. 2012, 86, 7393–7402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Simmonds, P. Recombination and selection in the evolution of picornaviruses and other mammalian positive-stranded RNA viruses. J. Virol. 2006, 80, 11124–11140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Wohlgemuth, N.; Honce, R.; Schultz-Cherry, S. Astrovirus evolution and emergence. Infect. Genet. Evol. 2019, 69, 30–37. [Google Scholar] [CrossRef]
  104. Van Dung, N.; Anh, P.H.; Van Cuong, N.; Hoa, N.T.; Carrique-Mas, J.; Hien, V.B.; Sharp, C.; Rabaa, M.; Berto, A.; Campbell, J. Large-scale screening and characterization of enteroviruses and kobuviruses infecting pigs in Vietnam. J. Gen. Virol. 2016, 97, 378–388. [Google Scholar] [CrossRef]
  105. Qiu, Z.; Wang, Z.; Zhang, B.; Zhang, J.; Cui, S. The prevalence of porcine teschovirus in the pig population in northeast of China. J. Virol. Methods 2013, 193, 209–214. [Google Scholar] [CrossRef] [PubMed]
  106. Zeng, S.; Wang, D.; Fang, L.; Ma, J.; Song, T.; Zhang, R.; Chen, H.; Xiao, S. Complete coding sequences and phylogenetic analysis of porcine bocavirus. J. Gen. Virol. 2011, 92, 784–788. [Google Scholar] [CrossRef] [Green Version]
  107. Blomström, A.-L.; Belák, S.; Fossum, C.; Fuxler, L.; Wallgren, P.; Berg, M. 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. Virus Res. 2010, 152, 59–64. [Google Scholar] [CrossRef]
  108. Brunborg, I.M.; Moldal, T.; Jonassen, C.M. Quantitation of porcine circovirus type 2 isolated from serum/plasma and tissue samples of healthy pigs and pigs with postweaning multisystemic wasting syndrome using a TaqMan-based real-time PCR. J. Virol. Methods 2004, 122, 171–178. [Google Scholar] [CrossRef]
  109. Monini, M.; Vignolo, E.; Ianiro, G.; Ostanello, F.; Ruggeri, F.M.; Di Bartolo, I. Detection of Torque teno sus virus in pork bile and liver sausages. Food Environ. Virol. 2016, 8, 283–288. [Google Scholar] [CrossRef] [PubMed]
  110. Pfankuche, V.M.; Bodewes, R.; Hahn, K.; Puff, C.; Beineke, A.; Habierski, A.; Osterhaus, A.D.; Baumgärtner, W. Porcine bocavirus infection associated with encephalomyelitis in a pig, Germany. Emerg. Infect. Dis. 2016, 22, 1310–1312. [Google Scholar] [CrossRef] [Green Version]
  111. Correa-Fiz, F.; Franzo, G.; Llorens, A.; Huerta, E.; Sibila, M.; Kekarainen, T.; Segalés, J. Porcine circovirus 2 (PCV2) population study in experimentally infected pigs developing PCV2-systemic disease or a subclinical infection. Sci. Rep. 2020, 10, 1–9. [Google Scholar] [CrossRef] [PubMed]
  112. Ouyang, T.; Zhang, X.; Liu, X.; Ren, L. Co-infection of swine with porcine circovirus type 2 and other swine viruses. Viruses 2019, 11, 185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  113. Rammohan, L.; Xue, L.; Wang, C.; Chittick, W.; Ganesan, S.; Ramamoorthy, S. Increased prevalence of torque teno viruses in porcine respiratory disease complex affected pigs. Vet. Microbiol. 2012, 157, 61–68. [Google Scholar] [CrossRef]
  114. Kekarainen, T.; Segalés, J. Torque teno sus virus in pigs: An emerging pathogen? Transbound. Emerg. Dis. 2012, 59, 103–108. [Google Scholar] [CrossRef]
  115. Rosell, C.; Segalés, J.; Plana-Duran, J.; Balasch, M.; Rodrıguez-Arrioja, G.; Kennedy, S.; Allan, G.; McNeilly, F.; Latimer, K.; Domingo, M. Pathological, immunohistochemical, and in-situ hybridization studies of natural cases of postweaning multisystemic wasting syndrome (PMWS) in pigs. J. Comp. Pathol. 1999, 120, 59–78. [Google Scholar] [CrossRef]
  116. Kekarainen, T.; Segalés, J. Porcine circovirus 2 immunology and viral evolution. Porc. Health Manag. 2015, 1, 1–6. [Google Scholar] [CrossRef] [Green Version]
  117. Segalés, J. Best practice and future challenges for vaccination against porcine circovirus type 2. Expert Rev. Vaccines 2015, 14, 473–487. [Google Scholar] [CrossRef]
  118. Fachinger, V.; Bischoff, R.; Jedidia, S.B.; Saalmüller, A.; Elbers, K. The effect of vaccination against porcine circovirus type 2 in pigs suffering from porcine respiratory disease complex. Vaccine 2008, 26, 1488–1499. [Google Scholar] [CrossRef]
  119. Ren, X.; Tao, Y.; Cui, J.; Suo, S.; Cong, Y.; Tijssen, P. Phylogeny and evolution of porcine parvovirus. Virus Res. 2013, 178, 392–397. [Google Scholar] [CrossRef] [PubMed]
  120. Streck, A.F.; Canal, C.W.; Truyen, U. Molecular epidemiology and evolution of porcine parvoviruses. Infect. Genet. Evol. 2015, 36, 300–306. [Google Scholar] [CrossRef] [PubMed]
  121. Treacy, D. A genetic analysis of the pedigree Landrace pig breed in Australia. Aust. J. Exp. Agric. 1976, 16, 76–81. [Google Scholar] [CrossRef]
  122. Pig Breeds in Australia. Available online: https://www.dpi.nsw.gov.au/__data/assets/pdf_file/0007/872683/Pig-breeds-brochure.pdf (accessed on 21 June 2021).
  123. Firth, C.; Charleston, M.A.; Duffy, S.; Shapiro, B.; Holmes, E.C. Insights into the evolutionary history of an emerging livestock pathogen: Porcine circovirus 2. J. Virol. 2009, 83, 12813–12821. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  124. Franzo, G.; Cortey, M.; Segalés, J.; Hughes, J.; Drigo, M. Phylodynamic analysis of porcine circovirus type 2 reveals global waves of emerging genotypes and the circulation of recombinant forms. Mol. Phylogenetics Evol. 2016, 100, 269–280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  125. Naccache, S.N.; Greninger, A.L.; Lee, D.; Coffey, L.L.; Phan, T.; Rein-Weston, A.; Aronsohn, A.; Hackett, J.; Delwart, E.L.; Chiu, C.Y. The perils of pathogen discovery: Origin of a novel parvovirus-like hybrid genome traced to nucleic acid extraction spin columns. J. Virol. 2013, 87, 11966–11977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  126. O’Dea, M.; Kabay, M.; Carr, J.; Wilcox, G.; Richards, R. Porcine circovirus–associated disease in weaner pigs in Western Australia. Aust. Vet. J. 2011, 89, 122–130. [Google Scholar] [CrossRef]
  127. Verani, M.; Casini, B.; Battistini, R.; Pizzi, F.; Rovini, E.; Carducci, A. One-year monthly monitoring of Torque teno virus (TTV) in river water in Italy. Water Sci. Technol. 2006, 54, 191–195. [Google Scholar] [CrossRef]
  128. Takayama, S.; Miura, T.; Matsuo, S.; Taki, M.; Sugii, S. Prevalence and persistence of a novel DNA TT virus (TTV) infection in Japanese haemophiliacs. Br. J. Haematol. 1999, 104, 626–629. [Google Scholar] [CrossRef]
  129. López-Lorenzo, G.; Díaz-Cao, J.M.; Prieto, A.; López-Novo, C.; López, C.M.; Díaz, P.; Rodríguez-Vega, V.; Díez-Baños, P.; Fernández, G. Environmental distribution of Porcine Circovirus Type 2 (PCV2) in swine herds with natural infection. Sci. Rep. 2019, 9, 1–8. [Google Scholar] [CrossRef] [Green Version]
  130. Patterson, A.; Opriessnig, T. Epidemiology and horizontal transmission of porcine circovirus type 2 (PCV2). Anim. Health Res. Rev. 2010, 11, 217–234. [Google Scholar] [CrossRef]
  131. Scott, A.; McCluskey, B.; Brown-Reid, M.; Grear, D.; Pitcher, P.; Ramos, G.; Spencer, D.; Singrey, A. Porcine epidemic diarrhoea virus introduction into the United States: Root cause investigation. Prev. Vet. Med. 2016, 123, 192–201. [Google Scholar] [CrossRef] [Green Version]
  132. Dee, S.; Neill, C.; Singrey, A.; Clement, T.; Cochrane, R.; Jones, C.; Patterson, G.; Spronk, G.; Christopher-Hennings, J.; Nelson, E. Modeling the transboundary risk of feed ingredients contaminated with porcine epidemic diarrhoea virus. BMC Vet. Res. 2016, 12, 1–12. [Google Scholar] [CrossRef] [Green Version]
  133. Nagai, M.; Wang, Q.; Oka, T.; Saif, L.J. Porcine sapoviruses: Pathogenesis, epidemiology, genetic diversity, and diagnosis. Virus Res. 2020, 286, 198025. [Google Scholar] [CrossRef] [PubMed]
  134. Bak, G.-Y.; Kang, M.-I.; Son, K.-Y.; Park, J.-G.; Kim, D.-S.; Seo, J.-Y.; Kim, J.-Y.; Alfajaro, M.M.; Soliman, M.; Baek, Y.-B. Occurrence and molecular characterization of Sapelovirus A in diarrhoea and non-diarrhoea feces of different age group pigs in one Korean pig farm. J. Vet. Med. Sci. 2016, 78, 1911–1914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  135. Tamura, K. Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+ C-content biases. Mol. Biol. Evol. 1992, 9, 678–687. [Google Scholar] [PubMed] [Green Version]
Figure 1. Similarity plot generated in SimPlot using the sequence MZ515502-PSV-AUS-2018 (Pig 46 Colon) as the query sequence against three other sequences using a 200 nucleotide sliding window at 20 nucleotide intervals and the F84 distance [64] model with the maximum likelihood method. Percentage identities at each analysis point were plotted on a line chart. For similarity plot analysis, the y-axis shows the percentage similarity between the reference sequences and the query sequence. Different colours are indexed for different reference sequences.
Figure 1. Similarity plot generated in SimPlot using the sequence MZ515502-PSV-AUS-2018 (Pig 46 Colon) as the query sequence against three other sequences using a 200 nucleotide sliding window at 20 nucleotide intervals and the F84 distance [64] model with the maximum likelihood method. Percentage identities at each analysis point were plotted on a line chart. For similarity plot analysis, the y-axis shows the percentage similarity between the reference sequences and the query sequence. Different colours are indexed for different reference sequences.
Viruses 13 01608 g001
Figure 2. Phylogenetic analysis of complete nucleotide sequence of VP1 capsid region of EV-G. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the General Time Reversible (GTR + G) [70] model with a bootstrapping of 1000 replicates. The analysis involved 39 reference sequences of the VP1 capsid region of EV-Gs genome and four EV-G sequences from this study. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. Sequences name indicated in red colour from Pig 45 colon and lung sample have been labelled with a black (⬤) and white circle (⭘), respectively while Pig 46 colon and lung sample have been labelled with a black (♦) and white rhombus (◊), respectively.
Figure 2. Phylogenetic analysis of complete nucleotide sequence of VP1 capsid region of EV-G. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the General Time Reversible (GTR + G) [70] model with a bootstrapping of 1000 replicates. The analysis involved 39 reference sequences of the VP1 capsid region of EV-Gs genome and four EV-G sequences from this study. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. Sequences name indicated in red colour from Pig 45 colon and lung sample have been labelled with a black (⬤) and white circle (⭘), respectively while Pig 46 colon and lung sample have been labelled with a black (♦) and white rhombus (◊), respectively.
Viruses 13 01608 g002
Figure 3. Phylogenetic analysis of the partial nucleotide sequences from the NS1 region of PBoV3. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the Hasegawa–Kishino–Yano (HKY + G + I) [84] model with a bootstrapping of 1000 replicates. The analysis included 21 reference sequences of the NS1 non-structural region of PBoV3 genome and two PBoV3 sequences from this study. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. The sequence name indicated in red colour from Pig 45 lung sample has been labelled with a white circle (⭘) while the Pig 46 colon sample has been labelled with a black rhombus (♦).
Figure 3. Phylogenetic analysis of the partial nucleotide sequences from the NS1 region of PBoV3. The nucleotide sequences were aligned and analysed using the maximum likelihood method in MEGA 7.0 [63] using the Hasegawa–Kishino–Yano (HKY + G + I) [84] model with a bootstrapping of 1000 replicates. The analysis included 21 reference sequences of the NS1 non-structural region of PBoV3 genome and two PBoV3 sequences from this study. The numbers at nodes represent bootstrap values and values above 60% are shown. Branch lengths are scaled according to the numbers of nucleotide substitutions per site. The sequence name indicated in red colour from Pig 45 lung sample has been labelled with a white circle (⭘) while the Pig 46 colon sample has been labelled with a black rhombus (♦).
Viruses 13 01608 g003
Table 1. The abundance of virus reads in the pig swab samples. Nucleotide identity percentage range with the closest reference sequence in the structural and non-structural region is shown. The structural and non-structural regions of TTSuV are overlapping and hence sequence identity with the closest reference has been shown in common.
Table 1. The abundance of virus reads in the pig swab samples. Nucleotide identity percentage range with the closest reference sequence in the structural and non-structural region is shown. The structural and non-structural regions of TTSuV are overlapping and hence sequence identity with the closest reference has been shown in common.
Abundance of Virus Reads (%)Identity with Closest Nucleotide Sequence in NCBI at Structural RegionIdentity with Closest Nucleotide Sequence in NCBI at Non-Structural Region
Samples
Pig 45, ColonPig 45, LungPig 46, ColonPig 46, Lung
PC45-BC24-AUS-2018PL45-BC23-AUS-2018PC46-BC26-AUS-2018PL46-BC25-AUS-2018
Virus Name
Porcine sapelovirus (PSV)0.35650.18050.98660.0235~80–88%~89%
Porcine enterovirus G (PEV-G)0.22500.35370.10560.0238~80%~83–85%
Porcine teschovirus (PTV)0.0160-0.02440.0072~82–86%~87–89%
Porcine astrovirus (PAstV)--0.0180-~75–81%~95%
Porcine bocavirus (PBoV)0.00510.03050.22460.0018~78–90%~93–97%
Porcine parvovirus 2 (PPV2)-0.0279-0.0001~94%~95%
Porcine parvovirus 7 (PPV7)0.0002-0.00380.0041~98%~99%
Porcine bufa virus (PBuV)0.0068-0.00220.0007~94–99%~99%
Adeno associated virus (AAV)0.00130.00570.00230.0002~80–97%~84–88%
Porcine circovirus 2 (PCV2)0.00020.000070.00070.0000399.82%100%
Torque teno sus virus k2a (TTSuVk2a)0.00240.0119-0.00047~97% (overlapped)
Torque teno sus virus k2b (TTSuVk2b)0.00580.0249-0.01399.92% (overlapped)
Total0.61930.63511.36820.0749
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bhatta, T.R.; Chamings, A.; Alexandersen, S. Exploring the Cause of Diarrhoea and Poor Growth in 8–11-Week-Old Pigs from an Australian Pig Herd Using Metagenomic Sequencing. Viruses 2021, 13, 1608. https://doi.org/10.3390/v13081608

AMA Style

Bhatta TR, Chamings A, Alexandersen S. Exploring the Cause of Diarrhoea and Poor Growth in 8–11-Week-Old Pigs from an Australian Pig Herd Using Metagenomic Sequencing. Viruses. 2021; 13(8):1608. https://doi.org/10.3390/v13081608

Chicago/Turabian Style

Bhatta, Tarka Raj, Anthony Chamings, and Soren Alexandersen. 2021. "Exploring the Cause of Diarrhoea and Poor Growth in 8–11-Week-Old Pigs from an Australian Pig Herd Using Metagenomic Sequencing" Viruses 13, no. 8: 1608. https://doi.org/10.3390/v13081608

APA Style

Bhatta, T. R., Chamings, A., & Alexandersen, S. (2021). Exploring the Cause of Diarrhoea and Poor Growth in 8–11-Week-Old Pigs from an Australian Pig Herd Using Metagenomic Sequencing. Viruses, 13(8), 1608. https://doi.org/10.3390/v13081608

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop