Next Article in Journal
A Dengue Virus Type 2 (DENV-2) NS4B-Interacting Host Factor, SERP1, Reduces DENV-2 Production by Suppressing Viral RNA Replication
Next Article in Special Issue
Phylodynamics of Influenza A/H1N1pdm09 in India Reveals Circulation Patterns and Increased Selection for Clade 6b Residues and Other High Mortality Mutants
Previous Article in Journal
Combinatorial Avidity Selection of Mosaic Landscape Phages Targeted at Breast Cancer Cells—An Alternative Mechanism of Directed Molecular Evolution
Previous Article in Special Issue
Proteomics Computational Analyses Suggest that the Antennavirus Glycoprotein Complex Includes a Class I Viral Fusion Protein (α-Penetrene) with an Internal Zinc-Binding Domain and a Stable Signal Peptide
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 11
Issue 9
10.3390/v11090786
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolution and Genetic Diversity of Porcine Circovirus 3 in China

by
Ye Chen
1,
Quanming Xu
2,
Hong Chen
1,
Xian Luo
1,
Qi Wu
1,
Chen Tan
1,
Qidong Pan
1 and
Ji-Long Chen
1,3,*
1
Key Laboratory of Fujian-Taiwan Animal Pathogen Biology, College of Animal Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
Fujian Agriculture and Forestry University, Fuzhou 350002, China
3
CAS Key Laboratory of Pathogenic Microbiology and Immunology, Institute of Microbiology, Chinese Academy of Sciences, Beijing 100101, China
*
Author to whom correspondence should be addressed.
Viruses 2019, 11(9), 786; https://0-doi-org.brum.beds.ac.uk/10.3390/v11090786
Submission received: 2 July 2019 / Revised: 21 August 2019 / Accepted: 21 August 2019 / Published: 27 August 2019
(This article belongs to the Special Issue Emerging Viruses: Surveillance, Prevention, Evolution and Control)
Browse Figures
Versions Notes

Abstract

:
The identification of a new circovirus (Porcine Circovirus 3, PCV3) has raised concern because its impact on swine health is not fully known. In Fujian Province in eastern China, even its circulating status and genetic characteristics are unclear. Here, we tested 127 tissue samples from swine from Fujian Province that presented respiratory symptoms. All of the PCV3 positive samples were negative for many other pathogens involved in respiratory diseases like PCV2, PRRSV, and CSFV, suggesting that PCV3 is potentially pathogenic. From phylogenetic analysis, PCV3 strains are divided into two main clades and five sub-clades; PCV3a-1, PCV3a-2, PCV3a-3, PCV3b-1, and PCV3b-2. Our identified strains belong to genotypes PCV3a-1, PCV3a-2, PCV3a-3, and PCV3b-2, indicating a high degree of genetic diversity of PCV3 in Fujian province until 2019. Interestingly, we found the time of the most recent common ancestor (tMRCA) of PCV3 was dated to the 1950s, and PCV3 has a similar evolutionary rate as PCV2 (the main epidemic genotypes PCV2b and PCV2d). In addition, positive selection sites N56D/S and S77T/N on the capsid gene are located on the PCV3 antigen epitope, indicating that PCV3 is gradually adaptive in swine. In summary, our results provide important insights into the epidemiology of PCV3.

1. Introduction

The Circovirus genus claimed veterinarians’ attention shortly after its discovery, most members infect avian species but swine also serve as natural hosts [1,2]. PCV1 and PCV2 are infectious agents in swine [3]; PCV1 was first detected as a contaminant of PK-15 cell cultures and has not been associated with clinical disease, PCV2 is a ubiquitous, economically important pathogen to the swine industry [4]. PCV2 infection leads to a diverse range of clinical diseases collectively termed PCV2-associated disease (PCVAD), which includes post-weaning multi-systemic wasting syndrome (PMWS), porcine dermatitis and nephropathy syndrome (PDNS), interstitial pneumonia, enteric disease, and reproductive failure [5]. Porcine Circovirus 3 (PCV3) is the third member of the Circovirus genus to infect swine [6], it has a circular ssDNA genome of approximatively 2000 bases containing three identified ORFs [6], although only ORF1 and ORF2 have been characterized. ORF1, located on the positive strand, apparently codes for a single replication-associated protein (Rep) of 296–297 aa [6]. ORF2 is on the negative strand and encodes the capsid protein (Cap). Despite the common genomic organization, PCV3 is distantly related to other known circoviruses, although a relatively close relationship with bat circoviruses has been suggested based on phylogenetic analysis [7]. There is also some unresolved opinion on the genotyping of PCV3. For example, Li et al. divide PCV3 into 3a and 3b, while Ouyang et al. suggest that PCV3 can be divided into PCV3a, 3b and 3c [7,8].
PCV3 was first identified in the USA in 2015; from tissues of animals presenting porcine dermatitis and nephropathy syndrome (PNDS) and reproductive disorders, but that had tested negative for PCV2; a metagenomics approach revealed an uncharacterized virus [6]. PCV3 has since been identified in Asia [9,10,11], South America, and Europe [12,13,14,15,16] in swine with various clinical syndromes [6]. Its comparable prevalence in healthy swine [17] and wild boar [18,19] however, raises questions about its role in swine pathogenesis.
In China, the presence of PCV3 has been confirmed in more than 24 provinces and regions, but its prevalence on swine farms (especially in asymptomatic herds) is not well established and specific genetic information is lacking, particularly in eastern China. The aim of this study was to determine the infection status and genetic characteristics of PCV3 in swine presenting respiratory symptoms on farms in Fujian province, eastern China. In addition, we aimed to elucidate the evolution and genetic diversity of PCV3 at the national scale.

2. Material and Methods

2.1. Tissue Samples, DNA Extraction, PCR Amplification, and Sequencing

One hundred and twenty-seven lung samples were collected from one hundred and twenty-seven swine on fifty-five farms in Fujian province, China. Samples were homogenated then subjected to three cycles of freeze–thaw followed by centrifugation at 5000 rpm for 10 min at 4 °C, after which supernatants were collected for DNA extraction. Viral DNA was extracted using the Virus Genomics DNA Isolation Kit (Tianlong Biotech, Suzhou, China) following the manufacturer’s instructions. PCR with pair specific primers (Table 1) was used for detecting PCV3, PCV2, CSFV, and PRRSV [20,21]. PCR conditions were as follows: for each sample 12.5 μL of 2× Taq Master mix (Vazyme Biotech, Shanghai, China), 9.5 μL of double distilled water (ddH2O), 1 μL of template DNA, 0.4 nmol/L forward and 0.4 nmol/L reverse primer were combined. Thermocycler conditions were 95 °C for 5 min, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 59 °C for 30 s, and extension at 72 °C for 30 s, with a final extension at 72 °C for 10 min. Samples were held at 4 °C until further use. Complete genomes were amplified using the Phanta Max Super-Fidelity DNA polymerase (Vazyme Biotech, Shanghai, China) with the same amplification conditions, subjected to 1% agarose gel electrophoresis, then sent to Tsingke (Nanjing, China) for sequencing. Sequences were assembled using the BioEdit software [22].

2.2. Sequence Collection, Alignment, and Phylogenetic Analysis

Three hundred and six complete PCV3 genomes, collected between 1996 and 2019 from different countries, were downloaded from GenBank at the National Center for Biotechnology Information (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/). All sequences were separated by Rep and Cap according coding region, then the coding area was spliced together (Supplementary Table S1). MUSCLE [23] was used to align sequences, and manual adjustments were made to the alignment using MEGA7.0 [24]. After poor quality sequences were removed, 272 reference sequences remained. The best fit nucleotide substitution model was selected using ModelFind in IQ-tree [25]. The maximum likelihood (ML) tree was reconstructed using RAxML (version 2.8.10) [26]. For nucleotide substitution, we used the General Time Reversible parameter plus the GAMMA distribution model (GTR + G) pre-estimated, with 1000 bootstrap replicates. To estimate the time elapsed since the most recent common ancestor and the rate of PCV3 evolution, we constructed a maximum clade credibility (MCC) tree using Markov chain Monte Carlo (MCMC) methods, using BEAST (v1.8.4) [27] and the Cap protein sequences. The substitution model was set to GTR + G with a relaxed molecular clock (lognormal), and the tree prior was coalescent (Bayesian skyline) according to model comparison built in Tracer (http://tree.bio.ed.ac.uk/software/tracer/). The total chain length was 1 × 108, with sampling every 10,000 steps. Two independent runs were performed and then combined using LogCombiner. The final MCC tree was annotated using TreeAnnotator and visualized using Figtree.

2.3. Selection Model Analysis

The ML tree based on full-length sequences (spliced ORF1 and ORF2 coding region) was uploaded to DataMonkey (www.datamonkey.org) to estimate branching and sites under selection. Selected sites were identified using four algorithms: single-likelihood ancestor counting (SLAC), mixed effects model of evolution (MEME), fixed effects likelihood (FEL), and fast unconstrained Bayesian approximation (FUBAR). Positively selected branches were identified using the adaptive branch-site REL test for episodic diversification (aBSREL) [28,29,30,31,32]. A site was considered to be under positive selection only if it satisfied at least two algorithms (p < 0.1 in SLAC, p < 0.05 in FEL and MEME, p > 0.9 in FUBAR). Cap sequences from strains PCV3a and PCV3b were used to predict protein structures in I-TASSER (https://zhanglab.ccmb.med.umich.edu/I-TASSER/) [33]. Structures were visualized using PyMoL.

3. Results

3.1. Characterization of PCV3 Strains Circulating in Fujian

Seven of the 127 samples were positive for PCV3 (positive rate of 5.51%), these samples are henceforth referred to as PCV3-CN-FJ1-2018, PCV3-CN-FJ2-2018, PCV3-CN-FJ22-2018, PCV3-CN-FJ27-2018, PCV3-CN-FJ33-2018, PCV3-CN-FJ37-2018, and PCV3-CN-FJ65-2018 (Supplementary Table S2). It should be noted that these seven samples were negative for PCV2, CSFV and PRRSV. All the PCV3 genomes we sequenced here were 2000 bp. Sequencing analysis of the full-length genomes, and Rep and Cap genes revealed maximum nucleotide sequence similarities amongst each other of 98.85–99.80%, 99.00–100% and 98.60–99.80% respectively, and to reference strains from GenBank of 91.24–99.90%, 98.2–100% and 96.7–100% respectively. Maximum amino acid sequence similarities for Rep and Cap were 99.3–100% and 98.1–100% within the strains, and compared to reference strains were 97.2–100% and 96.7–100% (Table 2), respectively. The 24th amino acid of the Cap protein of PCV3 is a key amino acid that distinguishes PCV3a from PCV3b; this amino acid is valine on the capsids of PCV3-CN-FJ2-2018, PCV3-CN-FJ27-2018 and PCV3-CN-FJ37-2018, and alanine on the capsids of PCV3-CN-FJ1-2018, PCV3-CN-FJ22-2018, PCV3-CN-FJ33-2018, and PCV3-CN-FJ65-2018.

3.2. Phylogenetic and Evolution Analysis of PCV3 in China and Worldwide

After removal of poor-quality data from an initial set of 306 PCV3 sequences obtained from GenBank, 272 reference sequences from the 14 countries remained (Figure 1). We used these, along with PCV3 sequences generated in our study, to reconstruct phylogenetic trees. Based on the ML tree (Figure 2), PCV3 can be classified into two main clades, PCV3a and PCV3b. At a higher level of resolution, PCV3a can be divided into PCV3a-1, PCV3a-2, and PCV3a-3, while PCV3b is comprised of PCV3b-1 and PCV3b-2. Globally, PCV3a strains outnumber PCV3b strains. Similarly, of the sequences from the PCV3 strains we isolated, five of seven were classified as PCV3a and two as PCV3b. From the tree, we found that the strains we isolated in Fujian were PCV3a-1, PCV3a-2, PCV3a-3 and PCV3b-2. In addition, these strains were close to the strains found in Jilin, Guangxi, Jiangxi, and Jiangsu provinces. The tMRCA was estimated to be 1954.25 with 95% highest probability density (HPD) between 1940.70 and 1984.90, while the evolutionary rate was estimated to be 2.41 × 10−4 with 95% HPD between 1.53 × 10−4 and 3.40 × 10−4 (Figure 3A).
The scaled effective population size (Neτ) for PCV3 was also estimated. Figure 3B shows that the population size increased from 1985 to 2012, then stabilized until 2016. From 2016 to 2019, the population size declined slightly but was still at a high level (about 5 × 102).

3.3. Selection and Amino Acid Function Analysis

Selective pressure analysis demonstrated a clear dominance of sites under negative pressure both in the Rep and Cap proteins [34]. We detected four positively selected amino acid sites, one at residue 122 in the Rep protein, and three at residues 24, 56, and 77 in Cap (Table 3). Residues 122 in Rep and 24 in Cap were detected by all methods (SLAC, FEL, FUBAR and MEME). In contrast, residues 56 and 77 were detected only by FEL and FUBAR. However, there was no positively selected branch in the PCV3 ML tree. Interestingly, several positively selected sites were associated with functional activities. Cap residues 24, 56, and 77 were in the epitope region and were classified as positively selected sites by Li et al. and Sun et al. [7,9] (Figure 4). Residue 122 in Rep was also described as positively selected by Li et al.

4. Discussion

The identification of PCV3, a new porcine circovirus resembling the significantly pathogenic PCV2, has raised a great interest in its evolution and epidemiology. In China, the number of emerging swine viruses has increased in the last years [7,35] raising two pressing questions. (i) Are swine novel hosts to these viruses? (ii) What is the epidemic character of PCV3 in China? Although PCV3 strains from most parts of China have been well studied, there is no systematic analysis of the prevalence and evolutionary characteristics of strains in eastern China, especially in Fujian Province. Fujian is a large aquaculture province with a subtropical monsoon climate that provides favorable conditions for the prevalence and evolution of disease.
Although some evidence supports an association between PCV3 infection and clinical disease [6], contradictory reports have been published [36], and more extensive studies are needed. Fu et al., 2018 found that in southern China 22.3% of samples were co-infected with PCV2. Ku et al., 2017 found that 45.4% of PCV3 positive samples were also positive for PCV2, while Zhang et al., 2018 found co-infection of PCV2 and PCV3 in 6.8% of the 265 clinical samples screened in that study. In a survey of samples from Shandong province we found 59.4% were infected with PCV3 alone and 39.3% were co-infected with PCV2 [37]. In this study, we found that all the samples from swine with respiratory symptoms were positive for PCV3 only. Although there are some people think PCV3 is suspicious of the pathogenicity of the herd [38], research on PCV3 pathogenicity should be further studied. Additionally, there are differing views of the phylogenetic analysis of PCV3; we found that PCV3 can be divided into two clades: PCV3a, PCV3b, and five sub-clades PCV3a-1, PCV3a-2, PCV3a-3, PCV3b-1, and PCV3b-2. This agrees with the report by Li et al. [7], but not fully with Ouyang et al., who are using ORF2 gene for analysis [8]. PCV3a is the most wide spread genotype in Fujian province consisting of PCV3a-1, PC3a-2 and PCV3a-3. Meanwhile, PCV3b-2but not PCV3b-1 was also found in Fujian province. Thus, in Fujian province, PCV3 has a high genetic diversity, in addition, the tMRCA of PCV3 is estimated later than in Saraiva et al. [39]. From previous important studies [7], it is known that an alanine at amino acid number 24 distinguishes PCV3a from PCV3b, which has a valine at that site. In our study the strain PCV3-CN-FJ27-2018 had a valine at position 24 but it belongs to the PCV3a group. This prompts us to propose that PCV3 subtypes might recombine with each other during infections, similar to PCV2. After MCMC analysis we found that the evolutionary rate of PCV3 exceeds that of canine parvovirus or porcine pseudorabies virus but is similar to PCV2b and PCV2d [40,41], this result is similar to that reported by Li et al. [7]. However, based on more reference sequences than used in the studies of Li et al., the tMRCA of PCV3 is earlier than reported in those studies. The estimated Neτ value of PCV3 is greater than that of Middle East respiratory syndrome coronavirus. We expect scaled effective population size Neτ to follow I/2β, where β is the equilibrium rate of transmission and I is the equilibrium number of infections. Although the duration of infection time of PCV3 in swine is unknown, given the large population size, its true value should also be similar to or greater than the estimated value.
Additionally, we found some amino acids on the antigen epitope that have been subject to positive selection during evolution. Similar to previous reports, amino acid 24 and 77 were found to be subjected to positive selection [7,8], however, amino acid 56 was first found in the Cap protein. Our results indicate the action of the host immune response in shaping PCV3 evolution. Swine with respiratory symptoms were positive only for PCV3, which may indicate that PCV3 has gradually adapted to swine and even become pathogenic to swine.
Although the number of sequences in PCV3 is increasing, we found that the analysis results are very similar to those reported previously after removal of poor quality sequences [7]. Pcv3 is still one of the DNA viruses with the highest evolutionary rate. In addition, it is of great concern that the inaccuracy of early samples’ epidemiological results may be caused by many laboratory pollutions. This situation will greatly misleading the evolution rate analysis of certain viruses, especially for new emerging viruses, so the evolutionary dynamics of PCV3 deserve further study. Our study has increased the information on the genetics and molecular characteristics of PCV3. These results will contribute to the evaluation of the relevance of PCV3 to the swine industry, and to the planning of effective control strategies in China.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/11/9/786/s1, Table S1: Reference sequences information used in this study, Table S2: Sequences information sequenced in this study.

Author Contributions

Data curation, C.T.; Formal analysis, Y.C. and X.L.; Funding acquisition, Q.W.; Methodology, Q.X. and Q.P.; Project administration, H.C.; Writing—original draft, J.-L.C.

Funding

This work was supported by National Key Research and Development Program of China (2017YFD0500203), Natural Science Foundation of Fujian Province of China (2018J01595).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tischer, I.; Gelderblom, H.; Vettermann, W.; Koch, M.A. A very small porcine virus with circular single-stranded DNA. Nature 1982, 295, 64–66. [Google Scholar] [CrossRef] [PubMed]
  2. Todd, D. Avian circovirus diseases: Lessons for the study of PMWS. Vet. Microbiol. 2004, 98, 169–174. [Google Scholar] [CrossRef] [PubMed]
  3. 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]
  4. Wang, H.J.; Gu, J.Y.; Xing, G.; Qiu, X.H.; An, S.T.; Wang, Y.X.; Zhang, C.; Liu, C.M.; Gong, W.J.; Tu, C.C.; et al. Genetic diversity of porcine circovirus type 2 in China between 1999–2017. Transbound. Emerg. Dis. 2019, 66, 599–605. [Google Scholar] [CrossRef] [PubMed]
  5. Ssemadaali, M.A.; Ilha, M.; Ramamoorthy, S. Genetic diversity of porcine circovirus type 2 and implications for detection and control. Res. Vet. Sci. 2015, 103, 179–186. [Google Scholar] [CrossRef] [PubMed]
  6. Palinski, R.; Piñeyro, P.; Shang, P.; Yuan, F.; Guo, R.; Fang, Y.; Byers, E.; Hause, B.M. A novel porcine circovirus distantly related to known circoviruses is associated with porcine dermatitis and nephropathy syndrome and reproductive failure. J. Virol. 2016, 91, e01879-16. [Google Scholar] [CrossRef] [PubMed]
  7. Li, G.; He, W.; Zhu, H.; Bi, Y.; Wang, R.; Xing, G.; Zhang, C.; Zhou, J.; Yuen, K.Y.; Gao, G.F.; et al. Origin, Genetic Diversity, and Evolutionary Dynamics of Novel Porcine Circovirus 3. Adv. Sci. 2018, 5, 1800275. [Google Scholar] [CrossRef] [PubMed]
  8. Ouyang, T.; Niu, G.; Liu, X.; Zhang, X.; Zhang, Y.; Ren, L. Recent progress on porcine circovirus type 3. Infect. Genet. Evol. 2019, 73, 227–233. [Google Scholar] [CrossRef]
  9. Sun, J.; Wei, L.; Lu, Z.; Mi, S.; Bao, F.; Guo, H.; Tu, C.; Zhu, Y.; Gong, W. Retrospective study of porcine circovirus 3 infection in China. Transbound. Emerg. Dis. 2018, 65, 607–613. [Google Scholar] [CrossRef]
  10. Ku, X.; Chen, F.; Li, P.; Wang, Y.; Yu, X.; Fan, S.; Qian, P.; Wu, M.; He, Q. Identification and genetic characterization of porcine circovirus type 3 in China. Transbound. Emerg. Dis. 2017, 64, 703–708. [Google Scholar] [CrossRef]
  11. Li, G.R.; Wang, H.J.; Wang, S.L.; Xing, G.; Zhang, C.; Zhang, W.Y.; Liu, J.; Zhang, J.Y.; Su, S.; Zhou, J.Y. Insights into the genetic and host adaptability of emerging porcine circovirus 3. Virulence 2018, 9, 1301–1313. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Franzo, G.; Legnardi, M.; Hjulsager, C.K.; Klaumann, F.; Larsen, L.E.; Segales, J.; Drigo, M. Full-genome sequencing of porcine circovirus 3 field strains from Denmark, Italy and Spain demonstrates a high within-Europe genetic heterogeneity. Transbound. Emerg. Dis. 2018, 65, 602–606. [Google Scholar] [CrossRef] [PubMed]
  13. Ye, X.; Berg, M.; Fossum, C.; Wallgren, P.; Blomstrom, A.L. Detection and genetic characterisation of porcine circovirus 3 from pigs in Sweden. Virus Genes 2018, 54, 466–469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Stadejek, T.; Wozniak, A.; Milek, D.; Biernacka, K. First detection of porcine circovirus type 3 on commercial pig farms in Poland. Transbound. Emerg. Dis. 2017, 64, 1350–1353. [Google Scholar] [CrossRef] [PubMed]
  15. Collins, P.J.; McKillen, J.; Allan, G. Porcine circovirus type 3 in the UK. Vet. Rec. 2017, 181, 599. [Google Scholar] [CrossRef]
  16. Tochetto, C.; Lima, D.A.; Varela, A.P.M.; Loiko, M.R.; Paim, W.P.; Scheffer, C.M.; Herpich, J.I.; Cerva, C.; Schmitd, C.; Cibulski, S.P.; et al. Full-Genome Sequence of Porcine Circovirus type 3 recovered from serum of sows with stillbirths in Brazil. Transbound. Emerg. Dis. 2018, 65, 5–9. [Google Scholar] [CrossRef]
  17. Zheng, S.; Shi, J.; Wu, X.; Peng, Z.; Xin, C.; Zhang, L.; Liu, Y.; Gao, M.; Xu, S.; Han, H.; et al. Presence of Torque teno sus virus 1 and 2 in porcine circovirus 3-positive pigs. Transbound. Emerg. Dis. 2018, 65, 327–330. [Google Scholar] [CrossRef]
  18. Franzo, G.; Tucciarone, C.M.; Drigo, M.; Cecchinato, M.; Martini, M.; Mondin, A.; Menandro, M.L. First report of wild boar susceptibility to Porcine circovirus type 3: High prevalence in the Colli Euganei Regional Park (Italy) in the absence of clinical signs. Transbound. Emerg. Dis. 2018, 65, 957–962. [Google Scholar] [CrossRef]
  19. Klaumann, F.; Dias-Alves, A.; Cabezon, O.; Mentaberre, G.; Castillo-Contreras, R.; Lopez-Bejar, M.; Casas-Diaz, E.; Sibila, M.; Correa-Fiz, F.; Segales, J. Porcine circovirus 3 is highly prevalent in serum and tissues and may persistently infect wild boar (Sus scrofa scrofa). Transbound. Emerg. Dis. 2019, 66, 91–101. [Google Scholar] [CrossRef]
  20. Sun, W.; Wang, W.; Xin, J.; Cao, L.; Zhuang, X.; Zhang, C.; Zhu, Y.; Zhang, H.; Qin, Y.; Du, Q.; et al. An epidemiological investigation of porcine circovirus 3 infection in dogs in the Guangxi Province from 2015 to 2017, China. Virus Res. 2019, 270, 197663. [Google Scholar] [CrossRef]
  21. Yue, F.X. Establishment and preliminary application of multiplex PCR detection methods for PCV2, PPV, PRV and PRRSV. Chin. Vet. Sci. 2008, 8, 691–696. [Google Scholar]
  22. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  23. Edgar, R.C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004, 5, 113. [Google Scholar] [CrossRef]
  24. 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]
  25. Lam-Tung, N.; Schmidt, H.A.; Arndt, V.H.; Bui Quang, M. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar]
  26. Stamatakis, A. RAxML Version 8: A tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  27. Drummond, A.J.; Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 2007, 7, 214. [Google Scholar] [CrossRef]
  28. Delport, W.; Poon, A.F.Y.; Frost, S.D.W.; Pond, S.L.K. Datamonkey 2010: A suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics 2010, 26, 2455–2457. [Google Scholar] [CrossRef]
  29. Kosakovsky Pond, S.L.; Frost, S.D. Not so different after all: A comparison of methods for detecting amino acid sites under selection. Mol. Biol. Evol. 2005, 22, 1208–1222. [Google Scholar] [CrossRef]
  30. Murrell, B.; Moola, S.; Mabona, A.; Weighill, T.; Sheward, D.; Pond, S.L.K.; Scheffler, K. FUBAR: A Fast, Unconstrained Bayesian AppRoximation for Inferring Selection. Mol. Biol. Evol. 2013, 30, 1196–1205. [Google Scholar] [CrossRef] [Green Version]
  31. Murrell, B.; Wertheim, J.O.; Moola, S.; Weighill, T.; Scheffler, K.; Pond, S.L.K. Detecting Individual Sites Subject to Episodic Diversifying Selection. PLoS Genet. 2012, 8, e1002764. [Google Scholar] [CrossRef]
  32. Smith, M.D.; Wertheim, J.O.; Weaver, S.; Murrell, B.; Scheffler, K.; Pond, S.L.K. Less Is More: An Adaptive Branch-Site Random Effects Model for Efficient Detection of Episodic Diversifying Selection. Mol. Biol. Evol. 2015, 32, 1342–1353. [Google Scholar] [CrossRef] [Green Version]
  33. Ambrish, R.; Alper, K.; Yang, Z. I-TASSER: A unified platform for automated protein structure and function prediction. Nat. Protoc. 2010, 5, 725–738. [Google Scholar]
  34. Franzo, G.; He, W.T. A Shift in Porcine Circovirus 3 (PCV-3) History Paradigm: Phylodynamic Analyses Reveal an Ancient Origin and Prolonged Undetected Circulation in the Worldwide Swine Population. Adv. Sci. in press.
  35. He, W.T.; Auclert, L.Z.; Zhai, X.; Wong, G.; Zhang, C.; Zhu, H.; Xing, G.; Wang, S.; He, W.; Li, K.; et al. Interspecies transmission, genetic diversity, and evolutionary dynamics of pseudorabies virus. J. Infect. Dis. 2018, 219, 1705–1715. [Google Scholar] [CrossRef]
  36. Klaumann, F.; Correa-Fiz, F.; Franzo, G.; Sibila, M.; Nunez, J.I.; Segales, J. Current Knowledge on Porcine circovirus 3 (PCV-3): A Novel Virus with a Yet Unknown Impact on the Swine Industry. Front. Vet. Sci. 2018, 5, 315. [Google Scholar] [CrossRef]
  37. Zheng, S.; Wu, X.; Zhang, L.; Xin, C.; Liu, Y.; Shi, J.; Peng, Z.; Xu, S.; Fu, F.; Yu, J.; et al. The occurrence of porcine circovirus 3 without clinical infection signs in Shandong Province. Transbound. Emerg. Dis. 2017, 64, 1337–1341. [Google Scholar] [CrossRef]
  38. Franzo, G.; Legnardi, M.; Tucciarone, C.M.; Drigo, M.; Klaumann, F.; Sohrmann, M.; SegalEs, J. Porcine circovirus type 3: A threat to the pig industry? Vet. Rec. 2018, 182, 83. [Google Scholar] [CrossRef]
  39. Saraiva, G.L.; Vidigal, P.M.P.; Fietto, J.L.R.; Bressan, G.C.; Silva Junior, A.; de Almeida, M.R. Evolutionary analysis of Porcine circovirus 3 (PCV3) indicates an ancient origin for its current strains and a worldwide dispersion. Virus Genes 2018, 54, 376–384. [Google Scholar] [CrossRef]
  40. Franzo, G.; Cortey, M.; Segales, 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. Phylogenet. Evol. 2016, 100, 269–280. [Google Scholar] [CrossRef] [Green Version]
  41. He, W.; Zhao, J.; Xing, G.; Li, G.; Wang, R.; Wang, Z.; Zhang, C.; Franzo, G.; Su, S.; Zhou, J. Genetic analysis and evolutionary changes of Porcine circovirus 2. Mol. Phylogenet. Evol. 2019, 139, 106520. [Google Scholar] [CrossRef]
Viruses 11 00786 g001 550
Figure 1. Map of PCV3 positive farms by country, and by province in China. Colored rectangles correspond to indicated PCV3 positive country and Chinese provinces. The yellow star indicates Fujian province.
Figure 1. Map of PCV3 positive farms by country, and by province in China. Colored rectangles correspond to indicated PCV3 positive country and Chinese provinces. The yellow star indicates Fujian province.
Viruses 11 00786 g001
Viruses 11 00786 g002 550
Figure 2. Maximum likelihood trees based on full-length PCV3 sequences. Trees were reconstructed using RAxML (Version 8.4.10) with the general time reversible plus GAMMA (GTR + G) distribution substitution model (1000 bootstraps). Blue dots represent sequenced strains from this study. Colored lines indicate country or region of origin.
Figure 2. Maximum likelihood trees based on full-length PCV3 sequences. Trees were reconstructed using RAxML (Version 8.4.10) with the general time reversible plus GAMMA (GTR + G) distribution substitution model (1000 bootstraps). Blue dots represent sequenced strains from this study. Colored lines indicate country or region of origin.
Viruses 11 00786 g002
Viruses 11 00786 g003 550
Figure 3. (A) Evolutionary rate for the PCV3 Cap gene. The evolutionary rate was estimated using BEAST (v1.8.4) with the GTR + G model and a relaxed lognormal molecular clock. The tree prior was set at coalescent: Bayesian skyline, and total chain length was 1 × 108 with sampling every 10,000 steps. (B) Bayesian skyline plot for the PCV3 Cap gene. The mean genetic diversity through time (Net) for the PCV3 ORF2 gene is shown by the black boldface line, while the 95% highest probability density (HPD) of Net is plotted on the y-axis (grey region).
Figure 3. (A) Evolutionary rate for the PCV3 Cap gene. The evolutionary rate was estimated using BEAST (v1.8.4) with the GTR + G model and a relaxed lognormal molecular clock. The tree prior was set at coalescent: Bayesian skyline, and total chain length was 1 × 108 with sampling every 10,000 steps. (B) Bayesian skyline plot for the PCV3 Cap gene. The mean genetic diversity through time (Net) for the PCV3 ORF2 gene is shown by the black boldface line, while the 95% highest probability density (HPD) of Net is plotted on the y-axis (grey region).
Viruses 11 00786 g003
Viruses 11 00786 g004 550
Figure 4. Predicted structure of PCV3 Cap protein. Cap protein structure was predicted using I-TASSER (https://zhanglab.ccmb.med.umich.edu/I-TASSER/). The yellow filled circles indicate the epitope of PCV3 cap protein, the red filled circles indicate amino acids subjected to positive selection, and green lines indicate α-helixes and β-folds. (A) PCV3a (MF318452:PCV3-BJ-2), (B) PCV3b (KX458235:2164).
Figure 4. Predicted structure of PCV3 Cap protein. Cap protein structure was predicted using I-TASSER (https://zhanglab.ccmb.med.umich.edu/I-TASSER/). The yellow filled circles indicate the epitope of PCV3 cap protein, the red filled circles indicate amino acids subjected to positive selection, and green lines indicate α-helixes and β-folds. (A) PCV3a (MF318452:PCV3-BJ-2), (B) PCV3b (KX458235:2164).
Viruses 11 00786 g004
Table
Table 1. Primers used in this study.
Table 1. Primers used in this study.
PrimersAmplicon Length (bp)
PCV3-D-FACTTAGAGAACGGACTTGTAACGAA649
PCV3-D-RAAATGAGACACAGAGCTATATTCAG
PCV3-1-FATTATGGATGCTCCTCATCGTG553
PCV3-1-RCATCTTCTCCGCAACTTCAGTC
PCV3-2-FGACTGAAGTTGCGGAGAAGATG789
PCV3-2-RCGGCACGAAAGAAGTTTGGATT
PCV3-3-FCCCACATGCGAGGGCGTTTACC895
PCV3-4-RCGAGGCCGCTTCATCATCCACT
PCV2-D-FAGAAGCTCTCTATCGGAG569
PCV2-D-RAAGGTTGAATTCTGGCCC
CSFV-D-FTAGGGTGGACGGGTGTCATAGAGT566
CSFV-D-RAAGCATATATTGCTGGAAGTAGCT
PRRSV-D-FGCCTCGTGTTGGGTGGCAGAA532
PRRSV-D-RCGCCCTAATTGAATAGGTGACTT
Classical swine fever virus (CSFV), porcine reproductive and respiratory syndrome virus (PRRSV).
Table
Table 2. Maximum nucleotide and amino acid sequence divergence of Porcine Circovirus 3 (PCV3).
Table 2. Maximum nucleotide and amino acid sequence divergence of Porcine Circovirus 3 (PCV3).
Isolated StrainsSequence Similarity among Isolated StrainsSequence Similarity between Isolated Strains and Reference Strain
Rep gene (nt)99.0–100%98.2–100%
Rep gene (aa)99.3–100%97.2–100%
Cap gene (nt)98.6–99.8%97.8–100%
Cap gene (aa)98.1–100%96.7–100%
Full-length (nt)98.85–99.80%97.26–99.90%
Replication (Rep), Capsid (Cap), amino acid (aa), nucleotide (nt).
Table
Table 3. Selection analysis of full-length PCV3.
Table 3. Selection analysis of full-length PCV3.
AAFELSLACFUBARMEME
dN-dSp-ValuedN-dSp-ValuedN-dSPost.Prow+p-Value
1226.5150.00213.80.0021711.7310.999>1000
320(24)2.7840.0276.290.05854.930.993>1000.02
373(77)3.440.0395.020.2854.1640.982>1000.1
352(56)3.7390.0555.660.2655.5780.987>1000.12

Share and Cite

MDPI and ACS Style

Chen, Y.; Xu, Q.; Chen, H.; Luo, X.; Wu, Q.; Tan, C.; Pan, Q.; Chen, J.-L. Evolution and Genetic Diversity of Porcine Circovirus 3 in China. Viruses 2019, 11, 786. https://0-doi-org.brum.beds.ac.uk/10.3390/v11090786

AMA Style

Chen Y, Xu Q, Chen H, Luo X, Wu Q, Tan C, Pan Q, Chen J-L. Evolution and Genetic Diversity of Porcine Circovirus 3 in China. Viruses. 2019; 11(9):786. https://0-doi-org.brum.beds.ac.uk/10.3390/v11090786

Chicago/Turabian Style

Chen, Ye, Quanming Xu, Hong Chen, Xian Luo, Qi Wu, Chen Tan, Qidong Pan, and Ji-Long Chen. 2019. "Evolution and Genetic Diversity of Porcine Circovirus 3 in China" Viruses 11, no. 9: 786. https://0-doi-org.brum.beds.ac.uk/10.3390/v11090786

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