Figures
Abstract
A study was conducted to investigate the serum virome of sows with and without stillbirths after farrowing. Sera from sows with at least one stillbirth or with normal litters were collected immediately after farrowing. Viral DNA was extracted from serum pools and submitted to high throughput sequencing. No differences in the proportion of virus-related reads were found in both groups (p > 0.05). A variety of viral DNA genomes were identified, mostly representative of three viral families: Anelloviridae, Circoviridae and Smacoviridae. Besides, a number of novel unclassified circular Rep-encoding single stranded DNA (CRESS DNA) viruses were also identified. These findings suggest that the presence of such viral genomes in sows’ sera bears no correlation with stillbirths’ occurrence; it seems likely that these constitute part of the normal serum microbiome of sows at farrowing.
Citation: Tochetto C, Muterle Varela AP, Alves de Lima D, Loiko MR, Mengue Scheffer C, Pinto Paim W, et al. (2020) Viral DNA genomes in sera of farrowing sows with or without stillbirths. PLoS ONE 15(3): e0230714. https://doi.org/10.1371/journal.pone.0230714
Editor: Douglas Gladue, Plum Island Animal Disease Center, UNITED STATES
Received: July 23, 2019; Accepted: March 6, 2020; Published: March 26, 2020
Copyright: © 2020 Tochetto et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The viral genomes described here were deposited in GenBank under accession numbers MH170056 to MH170073.
Funding: - CT was in receipt of a Master's' degree fellowship from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES). - PMR is a research fellow 1A from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), grant n° 303432/2015-2. - This study was also supported by CNPq grant n° 421104/2018-0.
Competing interests: The authors have declared that no competing interests exist.
Introduction
In swine reproduction, the term “stillbirth” refers to piglets with healthy appearance born dead at term. The average incidence of stillbirths in commercial farms may range between 3% to 8%, causing significant losses to swine producers [1–3]. Such condition is considered a multifactorial problem [4]; on occasions, stillbirths have been associated to a number of known pathogens, including bacteria, protozoa, and viruses. Among viruses, those which have been most often associated with stillbirths are porcine parvovirus (ungulate protoparvovirus 1, UPV1) [5], porcine circovirus type 2 (PCV2) [6], Aujeszky's disease virus (ADV), Influenza virus [7], porcine reproductive and respiratory syndrome virus (PRRSV) [8], encephalomyocarditis virus (EMCV) [9] and classical swine fever virus (CSFV) [10]. Nevertheless, frequently, no viruses or other known pathogens seem to play a role in such condition [1, 2, 11, 12].
To date, most studies focusing on viruses with potential association with stillbirths have used traditional approaches, such as serological tests on samples from sows or stillborn piglets. More recently, the polymerase chain reaction (PCR) has been employed in this kind of studies [1, 11]; however, this method is still unable to detect agents to which no previously published genome sequences are available. Currently, high throughput sequencing (HTS) has been applied to investigate the full content of nucleic acids in a sample with no need for previous knowledge of specific nucleic acid sequences. As such, HTS has become a suitable tool to investigate the participation of infectious agents, in the etiology of conditions such as stillbirths. In this study, an investigation was conducted in search for DNA viruses in the serum of sows. Serum samples were collected in commercial swine farms in Southern Brazil from sows soon after farrowing and submitted to high throughput sequencing. Associations between the detection of DNA viral genomes and the occurrence of stillbirths were investigated. This work led to the detection of a number of complete viral genomes, which are here reported.
Material and methods
Ethics statement
This study was approved by the Ethics Committee in the Use of Animals from the Veterinary Research Institute Desidério Finamor (CEUA-IPVDF)–protocol number 16/2015.
Sampling and study design
Just farrowed sows from six commercial piglet producing farms located in five municipalities in Southern Brazil were selected. The sows were derived from farms certified as free of classical swine fever, pseudorabies, brucellosis, tuberculosis and mange. All sows were vaccinated to porcine circovirus 2, porcine parvovirus (i.e. ungulate protoparvovirus), Haemophilus parasuis, Erysipelothrix rhusiopathiae, colibacilosis and leptospirosis. All animals were clinically healthy at the time of sampling. Blood samples were taken from the auricular vein using disposable syringes and transported to the laboratory on ice. Sera were separated from blood and processed according to previously described protocols [13]. Sows from each farm were divided in two groups. The first included sows which had at least one stillbirth at this farrowing (named “S”; n = 6); the second consisted of sows with no stillbirths (named “H”; n = 6). Each group consisted of a pool with different numbers of serum samples. On total, serum samples of 94 sows were collected and grouped in twelve pools (six “S” and six “H”) for sequencing (S1 Table). Sows which had abortions or mummified fetuses at farrowing were not included in this study.
Sample processing and sequencing
After centrifugation (~ 5,000 x g/5 min), serum samples were pooled and filtered through 0.45 and 0.22 μm filters (Millex-GV, PVDF, Millex®). These were then centrifuged on a sucrose cushion (25%) at an average 150,000 x g for 4 h at 4 °C. The obtained pellets were resuspended in 400 μL of milli-Q water and stored at -80 °C until further processing. Subsequently, the resuspended pellets were treated with DNase (2 U/μL, Turbo DNase Kit, Ambion) and RNase A (20 mg/mL, Invitrogen) to reduce the concentration of non-encapsidated nucleic acids [14]. The samples were then incubated for 2 h at 37 °C and each half (200 uL) used for either DNA or RNA extraction. Viral DNA was extracted using a standard phenol protocol [15]. The DNA was enriched by multiple displacement amplification (MDA) using ϕ29 DNA polymerase [16]. Viral RNA was extracted with TRIzol (Invitrogen, Carlsbad, CA) according to the manufacturer’s protocol and enriched by whole transcriptome amplification (WTA) using REPLI-g® WTA Single Cell Kit (Qiagen). Both preparations were then purified with AMPure XP magnetic beads (Agentcourt). The quality and quantity of the enriched nucleic acids was checked by microvolume spectrophotometry (Nanodrop, ThermoScientific, USA) and fluorimetry (Qubit, Invitrogen, USA), respectively. The quantity of enriched RNA obtained after purification was insufficient for HTS; therefore, only DNA was sequenced. Twelve libraries (one for each pool) were constructed using the Nextera DNA preparation kit (Illumina®) according to the manufacturers’ recommendations. The sequencing was performed using the Miseq v2 300 kit (2 x 150 paired-end).
Bioinformatics pipeline
The quality of the obtained reads was verified using FastQC software (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Low quality reads (phred quality score <20) were trimmed using PRINSEQ tool (v. 0.20.4) [17]. The trimmed reads were de novo assembled using metaSPAdes assembler (v. 3.9.0) with default settings [18]. Assembled contigs were compared to the viral protein database of GenBank using BLASTx (e-value cutoff of 10−5) with Blast2GO software. Sequences with the best BLAST score were selected and assigned into known viral families and when applicable, to the respective CRESS DNA viruses current groups. Viral hits classified as bacteriophages were further submitted to PHASTER [19] to confirm their identities. Annotation and visualization of circular genomes were performed in Geneious software (v.8.1.3). Genome coverage was obtained by mapping raw reads using low-sensitivity/fastest mode (Geneious, v.8.1.3).
Assessment of viral abundance
The number of raw reads matching selected viral sequences (contigs and concatenates previously classified by BLAST) were measured using the Geneious assembly tool, with all raw data output reads in the “low-sensitivity/fastest” mode. Virus-specific reads from each pooled sample were then normalized by the total number of viral reads generated in each pool. Differences between the data obtained from pools of sows with and without stillbirths were compared using the Wilcoxon T test. A p-value <0.05 was indicative of statistical significance. A Principal Component Analysis (PCA) was performed to compare the virome profile between the groups. All such analyses were carried out with the SPSS software version 22 (IBM, Armonk, NY, USA).
Phylogenetic analyses
Phylogenetic analyses were performed based on predicted amino acid sequences (for circular rep-encoding ssDNA viruses, CRESS DNA) or nucleotide sequences (for Torque teno sus virus–TTSuV). Multiple sequence alignments were performed with the MEGA7 software [20] using ClustalW [21] with default settings to CRESS DNA sequences and with modifications for TTSuV (gap open creation penalty of 10 and gap extension penalty of 5). A phylogenetic tree of TTSuV was inferred from p-distance matrix using neighbor-joining method with a bootstrap of 1000 interactions. For CRESS DNA, the phylogenetic tree was inferred by maximum likelihood method using the best fit substitution model (VT + G + I + F) determined by PhyML through AIC criterion with approximate likelihood branch support (aLRT) [22].
Results and discussion
Virome analysis
The twelve pools generated a total of 7,918,254 reads (S2 Table). The overall percentage of eukaryotic viral reads was 29.50% in sows with stillbirths and 29.91% in sows without stillbirths. Reads corresponding to bacteriophages represented 0.03% of obtained reads and were identified in only four pools (1S, 3S, 4S, 6H). One of the challenges of virome analysis by HTS is to obtain enough amounts of viral nucleic acid free of host or bacterial DNA contamination. This was overcome here, since the rates of viral reads were higher than previously reported in swine serum [23] and faeces [24].
Reads related to two eukaryotic viral families with DNA genomes (Anelloviridae and Circoviridae), plus a number of circular Rep-encoding single stranded DNA (CRESS DNA) viruses were identified (Fig 1). The most abundant reads were related to torque teno sus viruses (TTSuV), members of Anelloviridae family, which were detected in all examined pools. Although many studies have attempted to link TTSuVs to disease in swine, their role in pathological conditions, if any, is still under debate, as TTSuVs have been detected circulating in healthy and diseased pigs at high viral loads [25–27]. Regarding reproductive failure, fetuses infected with TTSuV have been already identified; however, no associations between TTSuV and the occurrence of abortion or other reproductive problems have, to date, been proven [28]. Therefore, the participation of such viruses in reproductive disease seems unlikely.
The number of reads corresponding to each family was normalized by the number of eukaryotic viral reads. Dark blue columns represent the average reads obtained in each of the groups. The numbers followed by a letter in parenthesis correspond to the identification of the pools analyzed in this study. Numbers correspond to each of the farms (1–6); letters S (stillbirth) or H (healthy) refer with the occurrence of stillbirths. Reads related of Circoviridae were identified in three pools: 1S, 3S and 3H.
The second most frequently detected reads were those of CRESS DNA viruses (Fig 1), which were detected in eleven out of the twelve serum pools. Currently, the International Committee on Taxonomy of Viruses (ICTV) classified CRESS DNA viruses associated with eukaryotic hosts into four families: Circoviridae, Genomoviridae, Geminiviridae, Nanoviridae, Bacilladnaviridae and Smacoviridae [29]. These have been detected in a wide variety of samples, including faeces of pigs from many countries [30–36]; however, this is the first description of CRESS DNA viruses in sows’ sera. It seems likely that such viruses constitute part of the normal swine microbiota. What would be the role for CRESS DNA viruses circulating in the blood of pregnant sows is a query that remain unanswered.
The third most frequently detected reads were representatives of the Circoviridae family (Porcine circovirus type 1 —PCV1—and Porcine circovirus type 3—PCV3). These were detected in three pools (1S, 3S, 3H), although in pools 3S and 3H these represent less than 1% of the reads (Fig 1). PCV1 reads were identified in low proportions in two pools: one of sows with stillbirths (3S) and other of sows without stillbirths (3H). On the other hand, PCV3 reads were identified in two pools from sows with stillbirths (1S and 3S). Nevertheless, no association could be drawn between detection of circovirus genomes and stillbirths. Subsequent studies revealed that PCV3 genomes were highly prevalent in sera of most sows sampled here and showed no association with the occurrence of stillbirths (Tochetto et al., 2019 submitted).
When comparing the viral serum profile on both groups (with and without stillbirths), no statistically significant differences (p > 0.05) were observed (Table 1). Likewise, no distinct patterns could be detected by PCA (Fig 2). These findings support the lack of association between serum viromes and the occurrence of stillbirths. It is important to note that due to the large diversity of virus observed, a “core” virome could not be identified in the samples, which in part may have hampered the identification of specific viral patterns according to the group. Nevertheless, some viruses were detected exclusively in the group of sows without stillbirths (porcine stool-associated circular virus, hudisavirus, PCV3 and unclassified CRESS DNA viruses). However, due to the low amount of reads recovered, no statistically significant differences were achieved between the presence of such genomes in sows with or without stillbirths (Table 1).
The first two principal components (PC1 and PC2) were maintained since they explained more than 62% of the observed variability. Black dots represent the pools of sera from sows with stillbirths (S); white diamonds represent the pools of sera from sows without stillbirths (H).
The viruses which are often associated with stillbirths, i.e. UPV1, PCV2, PRRVS, EMCV, CSFV and Influenza virus, were not detected in this study. It must be pointed that the sows were vaccinated to UPV1 and PCV2; therefore, genomes corresponding to these agents were not expected to be detected. Viral RNA genomes were not obtained since they could not be recovered from the samples. This might have been caused either by the true absence of viral RNA genomes in samples, or by its presence in such small quantities that it could have been missed by the extraction procedure. However, the protocol employed here is well stablished in the laboratory and has demonstrated efficiency at viral RNA extraction [37,38]. Considering the good health status of the sows during sampling, it is possible that RNA viruses were not circulating in the animals. Whereas PRRSV is exotic to the country [39], EMCV causes a disease of rare occurrence in swine [40] and CSFV has been eradicated in the sampled area [41]. Sows from which samples were collected were from certified CSFV-free farms. Therefore, these agents were not expected to be detected in such samples.
It must should be considered that HTS does not provide an actual quantification of the detected viromes. Biases such as those introduced by the enrichment method (MDA) must be born in mind, although the impact of MDA on HTS has been lessened and regarded as not significant enough to lead to data misinterpretation [42]. Circular ssDNA genomes are preferentially amplified by MDA, and this could have interfered with the apparent abundance of such viral genomes detected in the present study However, such abundance certainly contributed to provide additional credibility to the six new CRESS DNA genomes discovered here.
DNA viruses which are often associated with stillbirths, as UPV1 and PCV2, were not detected in this study. It must be pointed that the sows were vaccinated to UPV1 and PCV2; therefore, genomes corresponding to these agents were not expected to be detected.
It is important to highlight that, in the present study, the samples were examined in pools, which may hamper the detection of genomes present in low concentrations. It is possible that sequencing individual samples might provide some additional information. In addition, only viral genomes matching the GenBank virus database were evaluated here. Most of the reads (~ 72%) could not be related to previously published viral sequences; these may represent new genomes of viruses, as yet unknown, infectious agents. Such sequences should be subject of further investigation.
The findings reported here reveal a number of viral genomes circulating in the serum of sows at farrowing. These strongly suggest that there is no association between the presence of such viral genomes in sow’s sera at farrowing and the occurrence of stillbirths. Nevertheless, it is still possible that an as yet undetected virus (or other infectious agents), could putatively be associated to stillbirths, though just not present in sow’s sera at farrowing. Viremia is expected to precede transplacental infection [43]. However, infection might have occurred throughout the gestation period and not be detected in sows’ sera at farrowing.
Complete viral genomes
Altogether, a total of 20 complete viral genomes were recovered from sows either with or without stillbirths, including twelve TTSuV, six CRESS DNA viruses and two PCV3 genomes. PCV3 genomes were recovered from two different pools of sows with stillbirths (farms 1 and 3) and were described elsewhere [13]. However, no association could be detected between PCV3 infection and the occurrence of stillbirths.
Next, a brief report is provided on the full viral genomes detected in this study.
Anelloviruses.
All twelve genomes representative of members of the Anelloviridae corresponded to TTSuV: four TTSuV 1a, six TTSuV 1b and two TTSuV k2a. The three TTSuV species were identified in sera of sows with (S) and without (H) stillbirths (Table 2). The TTSuV genome sizes varied between 2.8 to 2.9 kb. Four open reading frames (ORF) were identified: ORF1, ORF2, ORF1/1 and ORF2/2 (Fig 3A). ORF 1/1 and ORF 2/2 share the same intron of ORF1 and ORF2, respectively. The UTR length varied between 706–819 nt among TTSuV 1a, 696–824 nt among TTSuV 1b, and 880–939 nt between TTSuV k2a. Typical domains in the UTR region as the TATA box (motif ATATAA) and the G-C rich region were also detected.
(A) Schematic representation of one out of the twelve TTSuV genome recovered in this study. (B) Phylogenetic analysis based on the entire ORF1 at nucleotide level. Neighbor-joining method with p-distance and a bootstrap of 1000 replicates. Drawings of pigs highlight the sequences recovered in this study.
Based on the criteria for TTSuV classification from ICTV (2015), multiple alignment analyses of the entire ORF1 were performed. Phylogenetic analyses were performed based on the entire ORF1. The twelve TTSuV ORF1 sequences reported here were compared to 40 previously available TTSuV sequences. In such phylogenetic reconstruction, six of the ORF1 sequences identified here clustered along with TTSuV 1b; four clustered with TTSuV 1a, whereas two clustered along with TTSuV k2a (Fig 3B).
The TTSuV recovered here share 36.6% to 84.7% nucleotide identity, revealing two distinct genera. The TTSuV 1a sequences share 64.1–82.5% nucleotide identity, whereas the TTSuV 1b sequences shared 65.1–84.7% nucleotide identity in relation to each other. The two TTSuV k2a sequences exhibited a nucleotide sequence identity of 71.5%.
CRESS DNA genomes.
Six complete CRESS DNA genomes were recovered from four different serum pools, collected in three of the farms (Table 3). Five of the full CRESS DNA genomes identified do not fall within any family previously reported CRESS DNA viruses; thus, these were named “porcine serum associated circular DNA viruses” (PoSCV-1 to -5). One genome representative of a member of a recently created family of CRESS DNA genomes, named Smacoviridae [29], genus Porprismacovirus (as in “porcine and primate smacovirus”), was identified and named Porcine associated porprismacovirus 3, 12/RS/BR (Table 3).
The CRESS DNA genomes recovered here ranged in length from 1,846 to 2,568 nt, with two major ORFs that encode the putative Rep and Cap proteins. Four types of genome architectures previously described [44] were identified: type I (n = 2), type II (n = 1), type IV (n = 2) and type VI (n = 1) (Fig 4). All genomes revealed a putative origin of replication (ori) that contains a nonanucleotide motif at the apex of a predicted stem-loop structure responsible for initiating the rolling-cycle replication [44]. A phylogenetic tree was reconstructed based on Rep amino acid sequences reported here, with the best matches of BLASTp in GenBank and those of representative CRESS DNA genomes, including members of Smacoviridae, Genomoviridae, Gemiviridae, Circoviridae and Nanoviridae (Fig 5). PoSCV-1 2A/RS/BR, PoSCV-2 2B/RS/BR, PoSCV-3 7A/RS/BR and PoSCV-5 8/RS/BR clustered within different groups of unclassified CRESS DNA. Porcine associated porprismacovirus 3 12/RS/BR (MH170061) clustered with members of Smacoviridae family, and PoSCV-4 7B/RS/BR clustered with sequences of Hudisavirus (Fig 5).
Phylogenetic tree was inferred using PhyML with the VT + G + I + F substitution model. Sequences recovered in this study are highlighted by pig drawings. The number represents the farm and the letter the group; e.g. 2S = farm 2, sows with stillbirths; 2H = farm 2, sows without stillbirths.
PoSCV-2 2B/RS/BR and PoSCV-4 7B/RS/BR exhibited a conserved nonanucleotide motif (NANTATTAC), while PoSCV-1 2A/RS/BR, PoSCV-3 7A/RS/BR and PoSCV-5 8/RS/BR displayed less common nonamers (Table 4). Rolling circle replication (RCR) motifs I, II, and III and superfamily 3 helicase (SF3) Walker A, B, and C motifs were also identified in all of them, except for PoSCV-2 2B/RS/BR and PoSCV-5 8/RS/BR, which lack the motif I (Table 4).
PoSCV-1 2A/RS/BR and PoSCV-3 7A/RS/BR were most closely to sequences of CRESS DNA detected in serum of Brazilian pigs (KU203355, KU203356) and in duck feces in New Zealand (NC030133, NC030134) (Fig 5). They shared both the RCR and SF3 helicase motifs described in the genome of faeces of ducks, named associated circular DNA virus 2 (DufCV-2), except for one variation in the second aa in the Walker-A motif in PoSCV 3 7A/RS/BR (Table 4).
The PoSCV-2 2B/RS/BR genome showed no significant nucleotide identity to any of the previously known CRESS DNA genomes available on GenBank database. Despite this, PoSCV-2 genome contains characteristics identified in others CRESS DNA viruses. This genome of 1,900 nt in length is bidirectionally organized, encoding a putative Cap on the virion-sense strand and a putative Rep on the complementary strand. The large intergenic region (181 nt) contains a predicted stem-loop located at the 5’ end of Rep, with a conserved nonanucleotide at the apex (TAGTATTAC) (Fig 4). The Rep protein have RCR motifs II and III and the SF3 helicase domains Walker A, Walker B and motif C (Table 4). The absence of the RCR motif I had already been reported in CRESS DNA viruses [32]. The function of the RCR motif I is not fully understood, but it is believed to be important in DNA binding and cleavage prior to RCR [32, 44, 45]. Additionally, the Rep of PoSCV-2 shares 63% amino acid identity with Odonata-associated circular virus-17 (KM598400), identified in dragonflies in USA.
Another unclassified CRESS DNA genome, named PoSCV-4 7B/RS/BR shares ~78% nucleotide identity with hudisavirus (MF351515). These viruses were recently discovered in samples of human diarrheal feces in Peruvian patients [46], but their pathogenic potential remains unclear. Similar to hudisavirus, PoSCV-4 contains an ambisense genome that encodes a putative Cap on the virion-sense and a Rep on the complementary sense (Fig 4). Indeed, PoSCV-4 exhibited similar characteristics with hudisavirus, such as the nanomer (NANNNTTAC) at the apex of a stem-loop, the RCR and SF3 helicase motifs (Table 4). Very similar RCR motifs were also described in an unclassified CRESS DNA recovered from bovine faeces (BofCV-1) in New Zealand (KT822230). The SF3 helicase motifs of PoSCV-4 were the same of MF351515 and BofCV-1, except for one change in the first amino acid of the Walker-A motif (G/A). Finally, PoSCV-5 8/RS/BR was most closely related to Porcine serum-associated circular virus (KU203352) identified in Brazilian pig sera, with which the Rep coding region shares 94% identity (coverage 100%; e-value 0.0).
Smacoviruses.
The recently proposed family of CRESS DNA viruses, Smacoviridae, comprises genomes that contain two major ORFs organized in opposite directions, encoding the rolling circle replication-associated protein (Rep) on the complementary strand and the capsid protein (Cap) on the virion-sense strand. Recently, smacoviruses were divided into six new genera, including one that represents smacoviruses recovered from pigs and primates, named Porprismacovirus (Porcine and primate smacovirus) [29]. Out of the six full CRESS DNA genomes identified here, one genome shares more than 82% genome-wide pairwise identity with Porcine associated porprismacovirus 3 (KC545227-30). Following the criteria for classification of smacovirus, that stablishes a threshold of 77% genome-wide pairwise identity for species demarcation, this genome was named Porcine associated porprismacovirus 3 12/RS/BR. The RCR and SF3 helicase motifs detected in this genome were the same identified in previously reported genomes representatives of the same species (Table 4). Smacovirus genomes recovered from pigs were reported in USA, South Korea and New Zealand [31,33–35]; however this is the first report of such virus genome in serum of pigs.
Conclusions
This study provides a picture of DNA viruses present in sera of just farrowed sows with or without cases of stillbirths, in which genomes of anelloviruses, circoviruses and CRESS DNA viruses were the most abundantly detected. No association could be established between any of the viral genomes here identified and the occurrence of stillbirths. Nevertheless, a number of novel viral genomes were identified; these seem to be part of the normal virome in the serum of farrowing sows.
Supporting information
S1 Table. Identification of serum pools.
Farm and city of origin of sows, occurrence/non-occurrence of stillbirths, and numbers of sows sampled in each pool are shown.
https://doi.org/10.1371/journal.pone.0230714.s001
(DOCX)
S2 Table. Number of reads obtained by high throughput sequencing.
Viral eukaryotic reads and bacteriophage reads were normalized by the number of viral reads in each pool.
https://doi.org/10.1371/journal.pone.0230714.s002
(DOCX)
Acknowledgments
The authors are grateful to Brazilian veterinary practitioners and farm owners who helped to collect the samples and farm specific information.
References
- 1. Silva GS, da Costa Lana MV, Dias GB, da Cruz RA, Lopes LL, Machado G, et al. Case-control study evaluating the sow's risk factors associated with stillbirth piglets in Midwestern in Brazil. Trop Anim Health Prod. 2015;47(2):445–9. pmid:25516004
- 2. Borges VF, Bernardi ML, Bortolozzo FP, Wentz I. Risk factors for stillbirth and foetal mummification in four Brazilian swine herds. Prev Vet Med. 2005;70(3–4):165–76. pmid:16023524
- 3. Lucia T Jr., Correa MN, Deschamps JC, Bianchi I, Donin MA, Machado AC, et al. Risk factors for stillbirths in two swine farms in the south of Brazil. Prev Vet Med. 2002;53(4):285–92. pmid:11937235
- 4. Vanderhaeghe C, Dewulf J, de Kruif A, Maes D. Non-infectious factors associated with stillbirth in pigs: a review. Anim Reprod Sci. 2013;139(1–4):76–88. pmid:23602009
- 5. Wolf VH, Menossi M, Mourao GB, Gatti MS. Molecular basis for porcine parvovirus detection in dead fetuses. Genet Mol Res. 2008;7(2):509–17. pmid:18561384
- 6. Oropeza-Moe M, Oropeza Delgado AJ, Framstad T. Porcine circovirus type 2 associated reproductive failure in a specific pathogen free (SPF) piglet producing herd in Norway: a case report. Porcine Health Manag. 2017;3:25. pmid:29085660
- 7.
Kirkwood RN, Althouse GC, Yaeger MJ, Carr J, Almond GW. Diseases of the Reproductive System. In: Zimmerman JJ, Karriker LA, Ramirez A, Schwartz KJ, Stevenson GW, editors. Diseases of Swine. 10th ed. USA: Wiley-Blackwell; 2012. p. 341–6.
- 8. Jeong J, Kang I, Park C, Kim S, Park SJ, Park KH, et al. A comparison of the severity of reproductive failure between single and dual infection with porcine reproductive and respiratory syndrome virus (PRRSV)-1 and PRRSV-2 in late-term pregnancy gilts. Transbound Emerg Dis. 2018;65(6):1641–7. pmid:29877065
- 9. Carocci M, Bakkali-Kassimi L. The encephalomyocarditis virus. Virulence. 2012;3(4):351–67. pmid:22722247
- 10. Blome S, Staubach C, Henke J, Carlson J, Beer M. Classical swine fever—an updated review. Viruses. 2017;9(4):86.
- 11. Maldonado J, Segales J, Martinez-Puig D, Calsamiglia M, Riera P, Domingo M, et al. Identification of viral pathogens in aborted fetuses and stillborn piglets from cases of swine reproductive failure in Spain. Vet J. 2005;169(3):454–6. pmid:15848788
- 12. Rocha DL, Alberton GC, Dos Santos JL. Identificação do circovírus suíno tipo 2 e do parvovírus suíno em fetos suínos natimortos e mumificados provenientes de granjas no Brasil. Ciência Animal Brasileira. 2010;11(3).
- 13. Tochetto C, Lima DA, Varela APM, Loiko MR, Paim WP, Scheffer CM, et al. Full-genome sequence of porcine circovirus type 3 recovered from serum of sows with stillbirths in Brazil. Transbound Emerg Dis. 2018;65(1):5–9. pmid:29027372
- 14. Shan T, Li L, Simmonds P, Wang C, Moeser A, Delwart E. The fecal virome of pigs on a high-density farm. J Virol. 2011;85(22):11697–708. pmid:21900163
- 15.
Sambrook J, Russell DWDW. Molecular cloning: a laboratory manual. 3rd ed. New York: Cold Spring Harbor Laboratory; 2001.
- 16. Dean FB, Nelson JR, Giesler TL, Lasken RS. Rapid amplification of plasmid and phage DNA using Phi 29 DNA polymerase and multiply-primed rolling circle amplification. Genome Res. 2001;11(6):1095–9. pmid:11381035
- 17. Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4. pmid:21278185
- 18. Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27(5):824–34. pmid:28298430
- 19. Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44(W1):W16–21. pmid:27141966
- 20. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4. pmid:27004904
- 21. Higgins DG. CLUSTAL V: multiple alignment of DNA and protein sequences. Methods Mol Biol. 1994;25:307–18. pmid:8004173
- 22. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21. pmid:20525638
- 23. Masembe C, Michuki G, Onyango M, Rumberia C, Norling M, Bishop RP, et al. Viral metagenomics demonstrates that domestic pigs are a potential reservoir for Ndumu virus. Virol J. 2012;9:218. pmid:23006778
- 24. Zhang B, Tang C, Yue H, Ren Y, Song Z. Viral metagenomics analysis demonstrates the diversity of viral flora in piglet diarrhoeic faeces in China. J Gen Virol. 2014;95(Pt 7):1603–11. pmid:24718833
- 25. Teixeira TF, Cibulski SP, dos Santos HF, Wendlant A, de Sales Lima FE, Schmidt C, et al. Torque teno sus virus 1 (TTSuV1) and 2 (TTSuV2) viral loads in serum of postweaning multisystemic wasting syndrome (PMWS)-affected and healthy pigs in Brazil. Res Vet Sci. 2015;101:38–41. pmid:26267087
- 26. Ramos N, Mirazo S, Botto G, Teixeira TF, Cibulski SP, Castro G, et al. High frequency and extensive genetic heterogeneity of TTSuV1 and TTSuVk2a in PCV2- infected and non-infected domestic pigs and wild boars from Uruguay. Vet Microbiol. 2018;224:78–87. pmid:30269794
- 27. Rogers AJ, Huang YW, Heffron CL, Opriessnig T, Patterson AR, Meng XJ. Prevalence of the novel torque teno sus virus species k2b from pigs in the United States and lack of association with post-weaning multisystemic wasting syndrome or mulberry heart disease. Transbound Emerg Dis. 2017;64(6):1877–83. pmid:27878979
- 28. Martinez-Guino L, Kekarainen T, Maldonado J, Aramouni M, Llorens A, Segales J. Torque teno sus virus (TTV) detection in aborted and slaughterhouse collected foetuses. Theriogenology. 2010;74(2):277–81. pmid:20452011
- 29. Varsani A, Krupovic M. Smacoviridae: a new family of animal-associated single-stranded DNA viruses. Arch Virol. 2018;163(7):2005–15. pmid:29572596
- 30. Cheung AK, Ng TF, Lager KM, Bayles DO, Alt DP, Delwart EL, et al. A divergent clade of circular single-stranded DNA viruses from pig feces. Arch Virol. 2013;158(10):2157–62. pmid:23612924
- 31. Cheung AK, Ng TF, Lager KM, Alt DP, Delwart E, Pogranichniy RM. Identification of several clades of novel single-stranded circular DNA viruses with conserved stem-loop structures in pig feces. Arch Virol. 2015;160(1):353–8. pmid:25248627
- 32. Cheung AK, Ng TF, Lager KM, Alt DP, Delwart EL, Pogranichniy RM. Identification of a novel single-stranded circular DNA virus in pig feces. Genome Announc. 2014;2(2)
- 33. Steel O, Kraberger S, Sikorski A, Young LM, Catchpole RJ, Stevens AJ, et al. Circular replication-associated protein encoding DNA viruses identified in the faecal matter of various animals in New Zealand. Infect Genet Evol. 2016;43:151–64. pmid:27211884
- 34. Kim AR, Chung HC, Kim HK, Kim EO, Nguyen VG, Choi MG, et al. Characterization of a complete genome of a circular single-stranded DNA virus from porcine stools in Korea. Virus genes. 2014;48(1):81–8. pmid:24170425
- 35. Sikorski A, Arguello-Astorga GR, Dayaram A, Dobson RC, Varsani A. Discovery of a novel circular single-stranded DNA virus from porcine faeces. Arch Virol. 2013;158(1):283–9. pmid:22972681
- 36. Sachsenroder J, Twardziok S, Hammerl JA, Janczyk P, Wrede P, Hertwig S, et al. Simultaneous identification of DNA and RNA viruses present in pig faeces using process-controlled deep sequencing. PloS One. 2012;7(4):e34631. pmid:22514648
- 37. Lima D, Cibulski SP, Finkler F, Teixeira T, Varela AP, Cerva C, et al. Fecal virome of healthy chickens reveals a large diversity of the eukaryote viral community, including novel circular single-stranded DNA viruses. J Gen Virol. 2017;98(4):690–703. pmid:28100302
- 38. Lima DA, Cibulski SP, Tochetto C, Varela APM, Finkler F, Teixeira TF, et al. The intestinal virome of malabsorption syndrome-affected and unaffected broilers through shotgun metagenomics. Virus Res. 2019;261:9–20. pmid:30543873
- 39.
OIE, World Organisation for Animal Health. World Animal Health Information Database (WAHIS) Interface. 2019. https://www.oie.int/wahis_2/public/wahid.php/Countryinformation/Animalsituation#BRA
- 40. Roehe PM, Rodrigues N.C., Oliveira S.C., Guizzardi I.I., Barcellos D.E.S.N. Encephalomyocarditis virus (EMCV) in swine in the state of Rio Grande do Sul, Brazil. Rev. microbiol. 1985:16(4).
- 41.
BRASIL. Ministério da Agricultura, Pecuária e Abastecimento. Situação sanitária das doenças de animais terrestres. 2017. http://www.agricultura.gov.br/assuntos/sanidade-animal-e-vegetal/saude-animal/arquivos-sisa/Situacao_sanitaria_WAHID_Brasil_2017_02.pdf
- 42. Parras-Moltó M, Rodríguez-Galet A, Suárez-Rodríguez P, López-Bueno A. Evaluation of bias induced by viral enrichment and random amplification protocols in metagenomic surveys of saliva DNA viruses. Microbiome. 2018;6(1):119. pmid:29954453
- 43. Pensaert M. Viremia and effect of fetal infection with porcine viruses with special reference to porcine circovirus 2 infection. Vet Microbiol. 2004;98(2):175–83. pmid:14741131
- 44. Rosario K, Duffy S, Breitbart M. A field guide to eukaryotic circular single-stranded DNA viruses: insights gained from metagenomics. Arch Virol. 2012;157(10):1851–71. pmid:22760663
- 45. Orozco BM, Hanley-Bowdoin L. Conserved sequence and structural motifs contribute to the DNA binding and cleavage activities of a geminivirus replication protein. J Biol Chem. 1998;273(38):24448–56.43. pmid:9733736
- 46. Altan E, Del Valle Mendoza J, Deng X, Phan TG, Sadeghi M, Delwart EL. Small circular rep-encoding single-stranded DNA genomes in peruvian diarrhea virome. Genome Announc. 2017;5(38).