Next Article in Journal
Essential Oils from Neotropical Piper Species and Their Biological Activities
Next Article in Special Issue
Hepatitis C Virus: Viral Quasispecies and Genotypes
Previous Article in Journal
In Silico Screening and In Vitro Activity Measurement of Javamide Analogues as Potential p38 MAPK Inhibitors
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 18
Issue 12
10.3390/ijms18122703
ijms-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

Phylodynamic and Genetic Diversity of Canine Parvovirus Type 2c in Taiwan

by
Yung-Cheng Lin
1,2,
Shu-Yun Chiang
3,
Hung-Yi Wu
4,
Jih-Hui Lin
5,
Ming-Tang Chiou
3,6,*,
Hsin-Fu Liu
1,7,8,* and
Chao-Nan Lin
3,6,*
1
Department of Medical Research, Mackay Memorial Hospital, Taipei 10449, Taiwan
2
Department of Nursing, Shu-Zen Junior College of Medicine and Management, Kaohsiung 82144, Taiwan
3
Department of Veterinary Medicine, College of Veterinary Medicine, National Pingtung University of Science and Technology, Pingtung 91201, Taiwan
4
Graduate Institute of Veterinary Pathobiology, College of Veterinary Medicine, National Chung-Hsing University, Taichung 40227, Taiwan
5
Center of Diagnostics and Vaccine Development, Centers for Disease Control, Taipei 11561, Taiwan
6
Animal Disease Diagnostic Center, College of Veterinary Medicine, National Pingtung University of Science and Technology, Pingtung 91201, Taiwan
7
Department of Bioscience and Biotechnology, National Taiwan Ocean University, Keelung 20224, Taiwan
8
Department of Nursing, National Taipei University of Nursing and Health Sciences, Taipei 11219, Taiwan
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2017, 18(12), 2703; https://doi.org/10.3390/ijms18122703
Submission received: 11 November 2017 / Revised: 4 December 2017 / Accepted: 9 December 2017 / Published: 13 December 2017
(This article belongs to the Special Issue Virus Comparative Genomics)
Browse Figures
Versions Notes

Abstract

:
Canine parvovirus type 2c (CPV-2c) emerged in 2000 and is known for causing a more severe disease than other CPV-2 variants in puppies. In 2015, the emerging CPV-2c variant was isolated in Taiwan and it subsequently became the predominant variant. To trace the evolution of Taiwanese CPV-2c, we compared complete VP2 genes of CPV-2c from Taiwan and sequences obtained from GenBank. The evolutionary rate of CPV-2c was estimated to be 4.586 × 10−4 substitutions per site per year (95% highest posterior density (HPD) was 3.284–6.076 × 10−4). The time to the most recent common ancestor (TMRCA) dated to 1990 (95% HPD: 1984–1996) and 2011 (95% HPD: 2010–2013) for the CPV-2c variant and Taiwanese isolates, respectively. The CPV-2c variant isolated from Taiwan was clustered with CPV-2c from China. This phylogenetic clade began to branch off in approximately 2010 (95% HPD was 3.823–6.497). Notably, two unique mutations of Taiwanese CPV-2c were found, Q383R and P410L. In summary, this is the first report on the genome evolution of CPV-2c in Taiwan, revealing that this CPV-2c variant shares a common evolutionary origin with strains from China. The demographic history inferred by the Bayesian skyline plot showed that the effective population of CPV-2c increased until 2006 and then slowly declined until 2011.

1. Introduction

Canine parvovirus type 2 (CPV-2) is a non-enveloped single-stranded linear DNA virus belonging to the genus Protoparvovirus, a member of the Parvoviridae family. The virus has two major open reading frames, one encoding the nonstructural proteins (NS1 and NS2) and the other encoding the capsid proteins (VP1 and VP2) [1]. Based on antigenical and genetic differences, CPV-2 was divided into three variable types (CPV-2a, 2b and 2c). The critical difference of variability was observed in residue 426 of the VP2 protein. The position was asparagine (Asn) in 2a, aspartic acid (Asp) in 2b and glutamic acid (Glu) in 2c [2,3]. The functions of the capsid protein VP2 include providing receptor binding, controlling the host range [4,5,6] and eliciting neutralizing antibodies [3].
The original CPV-2 was recognized in the late 1970s. Subsequently, it was replaced as the predominant form in the dog population of the United States by CPV-2a in approximately 1980. CPV-2b and CPV-2c were first detected in 1984 and 2000 in the United States and Italy, respectively [7]. Although CPV-2c has similar clinical signs to CPV-2a and CPV-2b, it has been reported to be associated with a more severe disease [8,9].
The mean emerging date estimated for the CPV-2 clade from feline panleukopenia virus was 1968. After CPV-2 diverged into CPV-2 and CPV-2a, the mean age of CPV-2a was 8 years after CPV-2 first emerged [1]. However, little information on the mean age of CPV-2c was available to develop an age estimate. Another retrospective analysis revealed that the oldest CPV-2c strain was identified in 1996 in Germany [10]. A retrospective analysis revealed that the frequency of the CPV-2 variants in Italy underwent a rapid fluctuation during 1995–2005, with CPV-2c rapidly replacing CPV-2b [11]. To date, the CPV-2c variant has been predominant worldwide except for Asian countries [8,12,13,14,15,16,17,18,19,20,21,22,23].
The first case of the CPV-2c variant in Taiwan was reported in 2015 and clinical isolates of CPV-2c significantly increased to the level of 54.6% to become the dominant variant [18]. Moreover, the novel CPV-2c variant was found to be distributed throughout Taiwan [18]. Previously, CPV-2a and CPV-2b were the major antigenic variants in Taiwan during the past two decades [18,21,22,23]. The capsid protein region on CPV-2 is highly antigenic and is the target of many neutralizing antibodies. Therefore, to trace the evolutionary history of Taiwanese CPV-2c variants, we conducted a comprehensive phylogenetic and evolutionary analysis based on full-length VP2 sequences from five Taiwanese CPV-2c strains and reference sequences retrieved from GenBank.

2. Results

2.1. Sequence Comparison and Amino Acid Sequence Analysis

A total of 156 complete CPV-2c VP2 gene sequences including Taiwanese isolates and reference strains from the National Center for Biotechnology Information (NCBI) database (Supplementary Materials Table S1) were analyzed with MEGA software version 7 [24]. The result showed the amino acid distance was significantly low (1.50–0.23%) between our isolates and reference strains. Furthermore, several substitution sites were found at A5G, F267Y, Y324I, Q370R Q383R and P410L. These substitution positions have only been observed in Taiwanese and Chinese isolates. Notably, Q383R and P410L were determined to be two unique substitution positions in Taiwanese isolates (KU244254, KX421787) (Table 1).

2.2. Evolutionary Rates and the Most Recent Common Ancestor of CPV-2 and CPV-2c

The dataset of CPV-2 and CPV-2c showed a positive correlation between genetic divergence and sampling time. Therefore, the dataset is suitable for phylogenetic molecular clock analysis in BEAST software. The uncorrelated exponential relaxed clock model and expansion growth population model were determined as the best fit models. The estimated nucleotide substitution rates of the VP2 region were 4.586 × 10−4 substitutions per site per year (95% highest posterior density (HPD) was 3.284–6.076 × 10−4) and 6.071 × 10−4 (95% HPD was 4.277–8.132) substitutions per site per year on CPV-2 and CPV-2c, respectively. The time to the most recent common ancestor (TMRCA) of CPV-2 was estimated to 1973 (95% HPD: 1963–1978). For the CPV-2c group, the overall TMRCA was dated to 1990 (95% HPD: 1984–1996) and to 2011 in Taiwanese isolates (95% HPD: 2010–2013). The evolutionary rate of the synonymous positions (third codon position) was 2.199, significantly higher than that of the nonsynonymous positions (first and second codon positions) (0.4) in the VP2 gene (Table 2). The strain C104-014/2015 (KU244254) was the first CPV-2c isolated in Taiwan and strains were clustered together with Chinese CPV-2c variants. The divergence of the Taiwanese lineage occurred in approximately 2010 (95% HPD was 2009–2012) (Figure 1).

2.3. Phylodynamics of CPV-2c

The phylodynamics of CPV-2c were estimated by a Bayesian skyline plot (BSP) based on the VP2 protein gene. The effective population size of CPV-2c was slowly rising until 2003, subsequently followed by a dramatic increase until 2005. After 2006, the effective population size decreased before remaining steady (Figure 2).

2.4. Selection Pressures in the CPV-2 VP2 Protein

The selection pressures of the VP2 protein in CPV-2 were estimated by the ratio of nonsynonymous substitutions (dN) to synonymous substitutions (dS). A dN/dS ratio of <1 represents negative selection; dN/dS ratio of 1 represents neutrality and dN/dS ratio of >1 represents positive selection. The overall mean dN/dS ratio in the CPV-2 VP2 gene was lower than 1 (0.108), implying that the gene was under negative selection. One positive selection was detected at codon 426 (Glu to Asn or Asp) using the FUBAR method (posterior probability was 0.97) (Table 3).

3. Discussion

CPV-2c is a new variant of CPV-2 that is now spreading worldwide. It causes gastroenteritis in puppy. The main clinical symptom among puppies in Taiwan is diarrhea or bloody diarrhea that also affects vaccinated dogs [18]. Studies have shown that both the CPV-2a and 2b variants constitute the prevalent CPV-2 field strains circulating in Taiwan in the past two decades [19,21,22,23]. However, this situation has since changed as the first case of CPV-2c was detected in 2015 and has since become the predominant strain (54.6%) in Taiwan [18].
Notably, a maximum likelihood (ML) tree revealed that CPV-2c isolated from Taiwan was clustered with Chinese strains, suggesting a closer phylogenetic relationship [18]. Results from a maximum clade credibility (MCC) tree (Figure 1) also support this observation, confirming that the recent CPV-2c isolate from Taiwan shares a common evolutionary origin with the Chinese CPV-2c and that the branching of that clade occurred in approximately 2010 (95% HPD was 3.823–6.497). Our results suggest that the Taiwanese isolate might have been imported from China and subsequently localized in Taiwan.
Selection pressure estimation results reveal a dN/dS ratio of 0.108. This means that the VP2 protein gene was under purifying selection. A positive selection was found in position 426 on the capsid VP2 protein gene. The residue 426 was located in the top of loop 4 and has been defined as the major mutation site for CPV-2 variants. Such results are consistent with previous reports [25,26]. Our BSP result shows a steady state for the effective population size of CPV-2c in recent years (Figure 2). This implies that the genetic evolution of the VP2 gene for CPV-2 has stabilized. By dating the emergent time of CPV-2, we estimated the TMRCA to be in 1973 (95% HPD: 1968–1978). This result is in agreement with reports that CPV had been in the canine population for several years before it was first recognized in 1978 [1,26,27]. The CPV-2c clade has a mean age of 26 years, suggesting that it originated in 1990 (95% HPD: 1984–1996), approximately 10 years before the first discovery of CPV-2c in 2000 [7]. Indeed, a retrospective analysis revealed that the oldest CPV-2c strain was identified in 1996 in Germany [10]. Comparison of the substitution rates of CPV-2 and CPV-2c revealed that our results (4.586 × 10−4 and 6.071 × 10−4 substitutions per site per year) are slightly higher than those observed in Pereira’s (1.2 × 10−4 substitutions per site per year) [27], Shackelton’s (1.7 × 10−4 substitutions per site per year) [1] and Hoelzer’s (2.2 × 10−4 substitutions per site per year) [26] studies. This may be due to the differing amounts of CPV-2 strains analyzed (163 CPV-2 strains and 156 CPV-2c strains in the present study vs. 78, 56 and 90 CPV-2a strains in Pereira’s, Shackelton’s and Hoelzer’s studies, respectively). Nevertheless, our data agree that CPV-2 is a rapidly evolved DNA virus, particularly the CPV-2c variant (6.071 × 10−4 substitutions per site per year) and that its evolution rate has approached that of RNA viruses at approximately 10−4 substitutions per site per year [1,28].
VP2 encodes a viral capsid protein that is the major structural protein of CPV-2 and is involved in providing receptor binding, controlling the host range [4,5,6] and eliciting neutralizing antibodies [3]. Therefore, a few mutations may result in increased pathogenicity [29]. The recent Taiwanese CPV-2c strains showed several amino acid changes compared with prototype CPV-2c strains (FJ222821) such as A5G, F267Y, Y324I and Q370R in the VP2 protein. The mutation of F267Y has also been observed in Vietnam [30], China [31,32,33,34], India [35], Uruguay [36,37] and Portugal [16]. Residue 267 is not exposed on the capsid surface [38,39] and substitutions in this position may not affect the antigenicity of the virus. However, a previous study showed that the binding of DNA to the internal surface of the parvovirus protein shell inflicts specific conformational changes on the protein [38]. Therefore, the function of residue 267 remains to be elucidated. Y324I with two nucleotide changes (TAT→ATT) has been reported in China [31,32,33,34,40,41,42], Korea [43,44], Thailand [45], Japan [46], Taiwan [19,23], India [35,47], Uruguay [36,37] and Hungary [48]. Residue 324 of VP2 is subject to positive selection [26] and is adjacent to a residue (amino acid 323) known to be involved in host range and tropism via canine transferrin receptor binding [49]. Therefore, the function of residue 324 remains to be elucidated. The substitution of Q370R in the Taiwanese CPV-2c strain is within the VP2 protein and the mutation of Q370R was also observed in the giant panda in China [50] and Chinese CPV-2c strains [32,42]. Residue 370 is located between residues 359 and 375, which is a flexible surface loop of the capsid protein that is adjacent to a double Ca2+-binding site; these were found to be essential for virus infectivity. Changes in these residues are correlated with the ability of the virus to hemagglutinate erythrocytes [51]. Therefore, whether substitution of Q370R causes antigenic change remains to be investigated. In addition to the aforementioned mutations, A5G substitutions were identified in our CPV-2c and Chinese strain (KR611522, KT162005), which had not been detected before 2014. Therefore, further studies focusing on the potential variant CPV-2c strains should be conducted to elucidate the relationship between the A5G substitution and viral pathogenicity. Q383R and P410L are two unique substitutions sites in Taiwanese isolates. Residue 383 is not located in any loop. Therefore, that substitution may not affect the antigenicity of the virus. The amino acid 410 residue is located in loop 4, whose top protrusion site consists of residues 421–428 and 433–443, indicating that 410 is not in the top protrusion [39]. Whether these two unique substitution sites benefit the localization of the virus in Taiwan remains unclear, although Q383R and P410L would not change the polarity of amino acids and is not located in the top protrusion. Therefore, ongoing research on the functional effect of these two unique sites and continued monitoring of the gene evolution of CPV-2c are necessary.

4. Materials and Methods

4.1. Ethics Statement and Study Design

Ethical approval was obtained from the Institutional Animal Care and Use Committee of National Pingtung University of Science and Technology. This was a retrospective study without intervention or obtaining extra clinical specimens from dogs. A total of 163 complete CPV-2 VP2 sequences including Taiwanese isolates and reference strains from the NCBI database were used in this study, which contained 156 CPV-2c strains (Supplementary Materials Table S1).

4.2. Complete VP2 Gene Amplification and Sequencing

The antigenic type of the CPV-2 isolate was determined by partial VP2 gene PCR, as described by Buonavoglia et al. [7]. Five CPV-2c positive samples were amplified in two overlapping fragments of 1436 and 1031 bp to obtain the full length of the VP2 gene sequence by using two sets of newly designed and reference primers [7]: first set: forward, 5′-CGGTGCAGGACAAGTAAAA-3′ and reverse (Hrev), 5′-CATTTGGATAAACTGGTGGT-3′; second set: forward (Hfor), 5′-CAGGTGATG AATTTGCTACA-3′ and reverse, 5′-AATCTTAAAATAATRTGTAATAAAC-3′. The purified PCR products were cloned into the T & ATM cloning Vector (Yeastern Biotech Co., Ltd., Taipei, Taiwan) and the presence of the desired insert in the recombinant plasmid DNA was confirmed by colony PCR with M13F and M13R universal primers. The target nucleotide sequences were determined in both orientations using an ABI automated sequencer (ABI 3730XL, Foster City, CA, USA). The viral nucleotide sequences determined in this study were deposited in GenBank with accession numbers KU244254 and KX421786-KX421789.

4.3. Phylodynamic Analysis

A full-length VP2 gene was aligned by T-coffee software [52]. The phylogenetic tree was constructed using the ML method in PhyML 3.0 (Available online: http://www.atgc-montpellier.fr/phyml/) [53]. Translation of the nucleotide sequences and estimation of the genetic distance were performed by MEGA software version 7 (Available online: http://www.megasoftware.net/) [24]. TempEst (TEMPoral Exploration of Sequences and Trees) was used to investigate the temporal signal [54]. The evolution rates of CPV-2 and CPV-2c were determined using the Bayesian Markov chain Monte Carlo (MCMC) method offered in BEAST v.1.8.4 (Available online: http://beast.bio.ed.ac.uk/) [55]. SRD06 was used as the best fit nucleotide substitution model because of its better resolution for coding regions for Bayesian analysis [56]. The demographic model, including constant size, expansion growth, Bayesian skyline and logistic growth, was used to estimate the evolutionary and population dynamics, under each as strict, lognormal relaxed and exponential relaxed molecular clock models [55]. The MCMC chains were run for a sufficient time to achieve convergence (Effective Sample Size > 200). The best fit demographic and clock model was estimated by Akaike’s information criterion (AICM) in the Tracer program v.1.6 [57]. An MCC tree was constructed by Tree Annotator v.1.8.4, with 10% burn-in. The final phylogenetic trees were edited by Figtree v.1.4.2. The accession numbers of the sequences used for evolutionary analysis in this study are listed in Supplementary Table S1.

4.4. Selection Pressure of CPV2 VP2 Protein Genes

To determine the selection pressures on the VP2 protein, we estimated the dN/dS ratio per site based on the maximum likelihood (ML) trees under the appropriate substitution model, using single-likelihood ancestor counting and fixed-effects likelihood methods, with the significance level set at 0.05. The Bayesian test for selection acting on individual sites was conducted using FUBAR, with the posterior probabilities set at 0.95 [58]. All the methods were implemented in the Datamonkey web server interface (http://www.datamonkey.org) [59,60].

5. Conclusions

This is the first report of CPV-2c genome evolution in Taiwan. Our results indicate that the recent CPV-2c isolate from Taiwan shares a common evolutionary origin with Chinese CPV-2c strains of and that a split of the Taiwanese lineage occurred in approximately 2010. The demographic history inferred from the BSP showed that the effective population of CPV-2c increased until 2006 and then slowly declined until 2011.

Supplementary Materials

Supplementary materials can be found at www.mdpi.com/1422-0067/18/12/2703/s1.

Acknowledgments

We thank Ling-Ling Chueh for comments that greatly improved the manuscript.

Author Contributions

Ming-Tang Chiou, Hsin-Fu Liu and Chao-Nan Lin conceived and designed the experiments; Yung-Cheng Lin, Shu-Yun Chiang, Jih-Hui Lin and Hsin-Fu Liu analyzed the experimental data and wrote the manuscript. Hung-Yi Wu, Ming-Tang Chiou, and Chao-Nan Lin managed the study, provided materials and reagents, contributed to the interpretation of the data, and co-wrote the manuscript. All of the authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Shackelton, L.A.; Parrish, C.R.; Truyen, U.; Holmes, E.C. High rate of viral evolution associated with the emergence of carnivore parvovirus. Proc. Natl. Acad. Sci. USA 2005, 102, 379–384. [Google Scholar] [CrossRef] [PubMed]
  2. Martella, V.; Decaro, N.; Buonavoglia, C. Evolution of CPV-2 and implication for antigenic/genetic characterization. Virus Genes 2006, 33, 11–13. [Google Scholar] [CrossRef] [PubMed]
  3. Lopez de Turiso, J.A.; Cortes, E.; Ranz, A.; Garcia, J.; Sanz, A.; Vela, C.; Casal, J.I. Fine mapping of canine parvovirus b cell epitopes. J. Gen. Virol. 1991, 72 Pt 10, 2445–2456. [Google Scholar] [CrossRef] [PubMed]
  4. Chang, S.F.; Sgro, J.Y.; Parrish, C.R. Multiple amino acids in the capsid structure of canine parvovirus coordinately determine the canine host range and specific antigenic and hemagglutination properties. J. Virol. 1992, 66, 6858–6867. [Google Scholar] [PubMed]
  5. Palermo, L.M.; Hueffer, K.; Parrish, C.R. Residues in the apical domain of the feline and canine transferrin receptors control host-specific binding and cell infection of canine and feline parvoviruses. J. Virol. 2003, 77, 8915–8923. [Google Scholar] [CrossRef] [PubMed]
  6. Parker, J.S.; Parrish, C.R. Canine parvovirus host range is determined by the specific conformation of an additional region of the capsid. J. Virol. 1997, 71, 9214–9222. [Google Scholar] [PubMed]
  7. Buonavoglia, C.; Martella, V.; Pratelli, A.; Tempesta, M.; Cavalli, A.; Buonavoglia, D.; Bozzo, G.; Elia, G.; Decaro, N.; Carmichael, L. Evidence for evolution of canine parvovirus type 2 in Italy. J. Gen. Virol. 2001, 82, 3021–3025. [Google Scholar] [CrossRef] [PubMed]
  8. Decaro, N.; Buonavoglia, C. Canine parvovirus: A review of epidemiological and diagnostic aspects, with emphasis on type 2c. Vet. Microbiol. 2012, 155, 1–12. [Google Scholar] [CrossRef] [PubMed]
  9. Decaro, N.; Desario, C.; Lucente, M.S.; Amorisco, F.; Campolo, M.; Elia, G.; Cavalli, A.; Martella, V.; Buonavoglia, C. Specific identification of feline panleukopenia virus and its rapid differentiation from canine parvoviruses using minor groove binder probes. J. Virol. Methods 2008, 147, 67–71. [Google Scholar] [CrossRef] [PubMed]
  10. Decaro, N.; Desario, C.; Addie, D.D.; Martella, V.; Vieira, M.J.; Elia, G.; Zicola, A.; Davis, C.; Thompson, G.; Thiry, E.; et al. The study molecular epidemiology of canine parvovirus, europe. Emerg. Infect. Dis. 2007, 13, 1222–1224. [Google Scholar] [CrossRef] [PubMed]
  11. Decaro, N.; Elia, G.; Martella, V.; Campolo, M.; Desario, C.; Camero, M.; Cirone, F.; Lorusso, E.; Lucente, M.S.; Narcisi, D.; et al. Characterisation of the canine parvovirus type 2 variants using minor groove binder probe technology. J. Virol. Methods 2006, 133, 92–99. [Google Scholar] [CrossRef] [PubMed]
  12. Kapil, S.; Cooper, E.; Lamm, C.; Murray, B.; Rezabek, G.; Johnston, L., 3rd; Campbell, G.; Johnson, B. Canine parvovirus types 2c and 2b circulating in north american dogs in 2006 and 2007. J. Clin. Microbiol. 2007, 45, 4044–4047. [Google Scholar] [CrossRef] [PubMed]
  13. Touihri, L.; Bouzid, I.; Daoud, R.; Desario, C.; El Goulli, A.F.; Decaro, N.; Ghorbel, A.; Buonavoglia, C.; Bahloul, C. Molecular characterization of canine parvovirus-2 variants circulating in tunisia. Virus Genes 2009, 38, 249–258. [Google Scholar] [CrossRef] [PubMed]
  14. Gallo Calderon, M.; Wilda, M.; Boado, L.; Keller, L.; Malirat, V.; Iglesias, M.; Mattion, N.; La Torre, J. Study of canine parvovirus evolution: Comparative analysis of full-length VP2 gene sequences from argentina and international field strains. Virus Genes 2012, 44, 32–39. [Google Scholar] [CrossRef] [PubMed]
  15. Muz, D.; Oguzoglu, T.C.; Timurkan, M.O.; Akin, H. Characterization of the partial VP2 gene region of canine parvoviruses in domestic cats from turkey. Virus Genes 2012, 44, 301–308. [Google Scholar] [CrossRef] [PubMed]
  16. Miranda, C.; Parrish, C.R.; Thompson, G. Epidemiological evolution of canine parvovirus in the portuguese domestic dog population. Vet. Microbiol. 2016, 183, 37–42. [Google Scholar] [CrossRef] [PubMed]
  17. Miranda, C.; Thompson, G. Canine parvovirus: The worldwide occurrence of antigenic variants. J. Gen. Virol. 2016, 97, 2043–2057. [Google Scholar] [CrossRef] [PubMed]
  18. Chiang, S.Y.; Wu, H.Y.; Chiou, M.T.; Chang, M.C.; Lin, C.N. Identification of a novel canine parvovirus type 2c in Taiwan. Virol. J. 2016, 13, 160. [Google Scholar] [CrossRef] [PubMed]
  19. Lin, C.N.; Chien, C.H.; Chiou, M.T.; Chueh, L.L.; Hung, M.Y.; Hsu, H.S. Genetic characterization of type 2a canine parvoviruses from Taiwan reveals the emergence of an ILE324 mutation in VP2. Virol. J. 2014, 11, 39. [Google Scholar] [CrossRef] [PubMed]
  20. Calderon, M.G.; Romanutti, C.; A, D.A.; Keller, L.; Mattion, N.; La Torre, J. Evolution of canine parvovirus in argentina between years 2003 and 2010: CPV2C has become the predominant variant affecting the domestic dog population. Virus Res. 2011, 157, 106–110. [Google Scholar] [CrossRef] [PubMed]
  21. Chang, W.L.; Chang, A.C.; Pan, M.J. Antigenic types of canine parvoviruses prevailing in Taiwan. Vet. Rec. 1996, 138, 447. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, H.C.; Chen, W.D.; Lin, S.L.; Chan, J.P.; Wong, M.L. Phylogenetic analysis of canine parvovirus VP2 gene in Taiwan. Virus Genes 2005, 31, 171–174. [Google Scholar] [CrossRef] [PubMed]
  23. Chou, S.J.; Lin, H.T.; Wu, J.T.; Yang, W.C.; Chan, K.W. Genotyping of canine parvovirus type 2 VP2 gene in southern Taiwan in 2011. Taiwan Vet. J. 2013, 39, 81–92. [Google Scholar]
  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] [PubMed]
  25. Silva, S.P.; Silva, L.; Rodrigues, E.D.L.; Cardoso, J.F.; Tavares, F.N.; Souza, W.M.; Santos, C.M.P.; Martins, F.M.S.; Jesus, I.S.; Brito, T.C.; et al. Full-length genomic and molecular characterization of canine parvovirus in dogs from north of brazil. Genet. Mol. Res. 2017, 16. [Google Scholar] [CrossRef] [PubMed]
  26. Hoelzer, K.; Shackelton, L.A.; Parrish, C.R.; Holmes, E.C. Phylogenetic analysis reveals the emergence, evolution and dispersal of carnivore parvoviruses. J. Gen. Virol. 2008, 89, 2280–2289. [Google Scholar] [CrossRef] [PubMed]
  27. Pereira, C.A.; Leal, E.S.; Durigon, E.L. Selective regimen shift and demographic growth increase associated with the emergence of high-fitness variants of canine parvovirus. Infect. Genet. Evol. 2007, 7, 399–409. [Google Scholar] [CrossRef] [PubMed]
  28. Duffy, S.; Shackelton, L.A.; Holmes, E.C. Rates of evolutionary change in viruses: Patterns and determinants. Nat. Rev. Genet. 2008, 9, 267–276. [Google Scholar] [CrossRef] [PubMed]
  29. Aldaz, J.; Garcia-Diaz, J.; Calleros, L.; Sosa, K.; Iraola, G.; Marandino, A.; Hernandez, M.; Panzera, Y.; Perez, R. High local genetic diversity of canine parvovirus from ecuador. Vet. Microbiol. 2013, 166, 214–219. [Google Scholar] [CrossRef] [PubMed]
  30. Nakamura, M.; Tohya, Y.; Miyazawa, T.; Mochizuki, M.; Phung, H.T.; Nguyen, N.H.; Huynh, L.M.; Nguyen, L.T.; Nguyen, P.N.; Nguyen, P.V.; et al. A novel antigenic variant of canine parvovirus from a vietnamese dog. Arch. Virol. 2004, 149, 2261–2269. [Google Scholar] [CrossRef] [PubMed]
  31. Zhang, R.; Yang, S.; Zhang, W.; Zhang, T.; Xie, Z.; Feng, H.; Wang, S.; Xia, X. Phylogenetic analysis of the VP2 gene of canine parvoviruses circulating in china. Virus Genes 2010, 40, 397–402. [Google Scholar] [CrossRef] [PubMed]
  32. Zhao, H.; Wang, J.; Jiang, Y.; Cheng, Y.; Lin, P.; Zhu, H.; Han, G.; Yi, L.; Zhang, S.; Guo, L.; et al. Typing of canine parvovirus strains circulating in North-East China. Transbound. Emerg. Dis. 2015, 64, 495–503. [Google Scholar] [CrossRef] [PubMed]
  33. Han, S.C.; Guo, H.C.; Sun, S.Q.; Shu, L.; Wei, Y.Q.; Sun, D.H.; Cao, S.Z.; Peng, G.N.; Liu, X.T. Full-length genomic characterizations of two canine parvoviruses prevalent in northwest china. Arch. Microbiol. 2015, 197, 621–626. [Google Scholar] [CrossRef] [PubMed]
  34. Xu, J.; Guo, H.C.; Wei, Y.Q.; Shu, L.; Wang, J.; Li, J.S.; Cao, S.Z.; Sun, S.Q. Phylogenetic analysis of canine parvovirus isolates from Sichuan and Gansu provinces of china in 2011. Transbound. Emerg. Dis. 2015, 62, 91–95. [Google Scholar] [CrossRef] [PubMed]
  35. Mukhopadhyay, H.K.; Matta, S.L.; Amsaveni, S.; Antony, P.X.; Thanislass, J.; Pillai, R.M. Phylogenetic analysis of canine parvovirus partial VP2 gene in India. Virus Genes 2013, 48, 89–95. [Google Scholar] [CrossRef] [PubMed]
  36. Perez, R.; Bianchi, P.; Calleros, L.; Francia, L.; Hernandez, M.; Maya, L.; Panzera, Y.; Sosa, K.; Zoller, S. Recent spreading of a divergent canine parvovirus type 2a (CPV-2a) strain in a CPV-2c homogenous population. Vet. Microbiol. 2012, 155, 214–219. [Google Scholar] [CrossRef] [PubMed]
  37. Perez, R.; Calleros, L.; Marandino, A.; Sarute, N.; Iraola, G.; Grecco, S.; Blanc, H.; Vignuzzi, M.; Isakov, O.; Shomron, N.; et al. Phylogenetic and genome-wide deep-sequencing analyses of canine parvovirus reveal co-infection with field variants and emergence of a recent recombinant strain. PLoS ONE 2014, 9, e111779. [Google Scholar] [CrossRef] [PubMed]
  38. Agbandje, M.; McKenna, R.; Rossmann, M.G.; Strassheim, M.L.; Parrish, C.R. Structure determination of feline Panleukopenia virus empty particles. Proteins 1993, 16, 155–171. [Google Scholar] [CrossRef] [PubMed]
  39. Tsao, J.; Chapman, M.S.; Agbandje, M.; Keller, W.; Smith, K.; Wu, H.; Luo, M.; Smith, T.J.; Rossmann, M.G.; Compans, R.W.; et al. The three-dimensional structure of canine parvovirus and its functional implications. Science 1991, 251, 1456–1464. [Google Scholar] [CrossRef] [PubMed]
  40. Zhong, Z.; Liang, L.; Zhao, J.; Xu, X.; Cao, X.; Liu, X.; Zhou, Z.; Ren, Z.; Shen, L.; Geng, Y.; et al. First isolation of new canine parvovirus 2a from Tibetan mastiff and global analysis of the full-length VP2 gene of canine parvoviruses 2 in china. Int. J. Mol. Sci. 2014, 15, 12166–12187. [Google Scholar] [CrossRef] [PubMed]
  41. Yi, L.; Tong, M.; Cheng, Y.; Song, W.; Cheng, S. Phylogenetic analysis of canine parvovirus VP2 gene in china. Transbound. Emerg. Dis. 2014, 63, e262–e269. [Google Scholar] [CrossRef] [PubMed]
  42. Geng, Y.; Guo, D.; Li, C.; Wang, E.; Wei, S.; Wang, Z.; Yao, S.; Zhao, X.; Su, M.; Wang, X.; et al. Co-circulation of the rare CPV-2c with unique Gln370Arg substitution, new CPV-2b with unique Thr440Ala substitution and new CPV-2a with high prevalence and variation in Heilongjiang province, northeast china. PLoS ONE 2015, 10, e0137288. [Google Scholar] [CrossRef] [PubMed]
  43. Jeoung, S.Y.; Ahn, S.J.; Kim, D. Genetic analysis of VP2 gene of canine parvovirus isolates in Korea. J. Vet. Med. Sci. 2008, 70, 719–722. [Google Scholar] [CrossRef] [PubMed]
  44. Yoon, S.H.; Jeong, W.; Kim, H.J.; An, D.J. Molecular insights into the phylogeny of canine parvovirus 2 (CPV-2) with emphasis on Korean isolates: A bayesian approach. Arch. Virol. 2009, 154, 1353–1360. [Google Scholar] [CrossRef] [PubMed]
  45. Phromnoi, S.; Sirinarumitr, K.; Sirinarumitr, T. Sequence analysis of VP2 gene of canine parvovirus isolates in thailand. Virus Genes 2010, 41, 23–29. [Google Scholar] [CrossRef] [PubMed]
  46. Soma, T.; Taharaguchi, S.; Ohinata, T.; Ishii, H.; Hara, M. Analysis of the VP2 protein gene of canine parvovirus strains from affected dogs in japan. Res. Vet. Sci. 2013, 94, 368–371. [Google Scholar] [CrossRef] [PubMed]
  47. Mittal, M.; Chakravarti, S.; Mohapatra, J.K.; Chug, P.K.; Dubey, R.; Narwal, P.S.; Kumar, A.; Churamani, C.P.; Kanwar, N.S. Molecular typing of canine parvovirus strains circulating from 2008–2012 in an organized kennel in India reveals the possibility of vaccination failure. Infect. Genet. Evol. 2014. [Google Scholar] [CrossRef] [PubMed]
  48. Csagola, A.; Varga, S.; Lorincz, M.; Tuboly, T. Analysis of the full-length VP2 protein of canine parvoviruses circulating in Hungary. Arch. Virol. 2014, 159, 2441–2444. [Google Scholar] [CrossRef] [PubMed]
  49. Hueffer, K.; Parrish, C.R. Parvovirus host range, cell tropism and evolution. Curr. Opin. Microbiol. 2003, 6, 392–398. [Google Scholar] [CrossRef]
  50. Guo, L.; Yang, S.L.; Chen, S.J.; Zhang, Z.; Wang, C.; Hou, R.; Ren, Y.; Wen, X.; Cao, S.; Guo, W.; et al. Identification of canine parvovirus with the Q370R point mutation in the VP2 gene from a giant panda (Ailuropoda melanoleuca). Virol. J. 2013, 10, 163. [Google Scholar] [CrossRef] [PubMed]
  51. Simpson, A.A.; Chandrasekar, V.; Hebert, B.; Sullivan, G.M.; Rossmann, M.G.; Parrish, C.R. Host range and variability of calcium binding by surface loops in the capsids of canine and feline parvoviruses. J. Mol. Biol. 2000, 300, 597–610. [Google Scholar] [CrossRef] [PubMed]
  52. Notredame, C.; Higgins, D.G.; Heringa, J. T-coffee: A novel method for fast and accurate multiple sequence alignment. J. Mol. Biol. 2000, 302, 205–217. [Google Scholar] [CrossRef] [PubMed]
  53. Guindon, S.; Dufayard, J.F.; 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, 307–321. [Google Scholar] [CrossRef] [PubMed]
  54. Rambaut, A.; Lam, T.T.; Max Carvalho, L.; Pybus, O.G. Exploring the temporal structure of heterochronous sequences using tempest (formerly Path-O-Gen). Virus Evol. 2016, 2, vew007. [Google Scholar] [CrossRef] [PubMed]
  55. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian phylogenetics with beauti and the beast 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [PubMed]
  56. Shapiro, B.; Rambaut, A.; Drummond, A.J. Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences. Mol. Biol. Evol. 2006, 23, 7–9. [Google Scholar] [CrossRef] [PubMed]
  57. Baele, G.; Lemey, P.; Bedford, T.; Rambaut, A.; Suchard, M.A.; Alekseyenko, A.V. Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol. Biol. Evol. 2012, 29, 2157–2167. [Google Scholar] [CrossRef] [PubMed]
  58. Murrell, B.; Moola, S.; Mabona, A.; Weighill, T.; Sheward, D.; Kosakovsky Pond, S.L.; Scheffler, K. Fubar: A fast, unconstrained bayesian approximation for inferring selection. Mol. Biol. Evol. 2013, 30, 1196–1205. [Google Scholar] [CrossRef] [PubMed]
  59. Poon, A.F.; Frost, S.D.; Pond, S.L. Detecting signatures of selection from DNA sequences using datamonkey. Methods Mol. Biol. 2009, 537, 163–183. [Google Scholar] [PubMed]
  60. Delport, W.; Poon, A.F.; Frost, S.D.; Kosakovsky Pond, S.L. Datamonkey 2010: A suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics 2010, 26, 2455–2457. [Google Scholar] [CrossRef] [PubMed]
Ijms 18 02703 g001 550
Figure 1. Maximum clade credibility (MCC) tree of CPV-2 inferred from 163 complete VP2 sequences. The MCC tree was constructed with 10% burn-in by Tree Annotator v 1.8 implemented in the BEAST software package. Red text represents the Taiwanese CPV-2c isolate in the present study. Numbers beside the branches are posterior probability values and branch time. Only posterior probability values above 0.95 are shown.
Figure 1. Maximum clade credibility (MCC) tree of CPV-2 inferred from 163 complete VP2 sequences. The MCC tree was constructed with 10% burn-in by Tree Annotator v 1.8 implemented in the BEAST software package. Red text represents the Taiwanese CPV-2c isolate in the present study. Numbers beside the branches are posterior probability values and branch time. Only posterior probability values above 0.95 are shown.
Ijms 18 02703 g001
Ijms 18 02703 g002 550
Figure 2. Bayesian skyline plot for the complete VP2 gene of the CPV-2. The x-axis is in units of year before 2016 and the y-axis represents the virus effective population size. The thicker bold line represents the median estimate of the effective number of infections over time and the thinner blue lines indicate the upper and lower bounds of the 95% HPD.
Figure 2. Bayesian skyline plot for the complete VP2 gene of the CPV-2. The x-axis is in units of year before 2016 and the y-axis represents the virus effective population size. The thicker bold line represents the median estimate of the effective number of infections over time and the thinner blue lines indicate the upper and lower bounds of the 95% HPD.
Ijms 18 02703 g002
Table
Table 1. Characteristics of the Amino Acid Substitution of CPV-2c Worldwide.
Table 1. Characteristics of the Amino Acid Substitution of CPV-2c Worldwide.
CountryYearAmino Acid at PositionAccession Numbers
5267324370383410440
Taiwan2015AY aI aR aR aL aTKU244254
2015AY aI aR aR aPTKX421787
2015G aY aI aR aQPTKX421786
2016G aY aI aR aQPTKX421788, KX421789
China2009AYIQQPTGU380303
2009AFYQQPTGU380305
2014GYIRQPTKR611522, KT162005
Argentina2008–2009AFYQQPTJF414818–JF414820
2009–2010AFYQQPAJF414821–JF414825
2010AFYQQPTJF414826
Belgium2008AFYQQPTFJ005247
Ecuador2012AFYQQPTKF149962–KF149965,
KF149967–KF149971,
KF149984,
2012AFYQQPSKF149966
Germany1997–1999AFYQQPTFJ005196–FJ005204
Greece2008–2009AFYQQPTGQ865518, GQ865519
Italy2000–2010AFYQQPTFJ222821, FJ005195,
FJ005205–FJ005213,
FJ005215–FJ005234,
FJ005237–FJ005245,
FJ005248–FJ005251,
HQ025913
India2014AFYQQPTKP071956
Spain2006–2008AFYQQPTFJ005214, FJ005246
Uruguay2006–2011AFYQQPTKC196079–KC196093,
KC196080–KC196093,
KC196095–KC196109,
KC196111–KC196114,
KM457103–KM457110,
KM457115–KM457131,
KM457142
2010AFYQQPAKC196094
USA2007AFYQQPTFJ005235
2007AFYQQPAFJ005236
Vietnam2004AFYQQPTAB120727
a Bold represents the amino acid differing from the prototype CPV-2c strain (FJ222821).
Table
Table 2. Mean Relative Evolutionary Rates for Codon Positions and TMRCA in VP2 gene.
Table 2. Mean Relative Evolutionary Rates for Codon Positions and TMRCA in VP2 gene.
TMRCASubstitution Rates Sub/Site/Year (10−4)Mean Relative Substitution RateSE of Mean
VP2 gene 4.586 (3.284~6.076)
CPV-21973 (1963–1978) a
CPV-2c1990 (1984–1996)6.071 (4.277~8.132)
CPV-2c Taiwanese strains2011 (2010–2013)
1st + 2nd codon position 0.4 (0.302~0.495)5.215 × 10−4
3rd codon position 2.199 (2.011~2.396)1.043 × 10−3
a ( ) Lower and upper 95% of highest posterior density (HPD).
Table
Table 3. Selection Pressure Detected in CPV-2 VP2 Protein Gene.
Table 3. Selection Pressure Detected in CPV-2 VP2 Protein Gene.
Positively Selected SitesNo. of Negatively Selected SitesMean dN/dS
SLAC aFEL aFUBA bSLAC aFEL aFUBAR b
CPV-2 VP2 geneNonNon426846310.108
a p value of <0.05. b Posterior probability of ≥0.95.

Share and Cite

MDPI and ACS Style

Lin, Y.-C.; Chiang, S.-Y.; Wu, H.-Y.; Lin, J.-H.; Chiou, M.-T.; Liu, H.-F.; Lin, C.-N. Phylodynamic and Genetic Diversity of Canine Parvovirus Type 2c in Taiwan. Int. J. Mol. Sci. 2017, 18, 2703. https://doi.org/10.3390/ijms18122703

AMA Style

Lin Y-C, Chiang S-Y, Wu H-Y, Lin J-H, Chiou M-T, Liu H-F, Lin C-N. Phylodynamic and Genetic Diversity of Canine Parvovirus Type 2c in Taiwan. International Journal of Molecular Sciences. 2017; 18(12):2703. https://doi.org/10.3390/ijms18122703

Chicago/Turabian Style

Lin, Yung-Cheng, Shu-Yun Chiang, Hung-Yi Wu, Jih-Hui Lin, Ming-Tang Chiou, Hsin-Fu Liu, and Chao-Nan Lin. 2017. "Phylodynamic and Genetic Diversity of Canine Parvovirus Type 2c in Taiwan" International Journal of Molecular Sciences 18, no. 12: 2703. https://doi.org/10.3390/ijms18122703

APA Style

Lin, Y. -C., Chiang, S. -Y., Wu, H. -Y., Lin, J. -H., Chiou, M. -T., Liu, H. -F., & Lin, C. -N. (2017). Phylodynamic and Genetic Diversity of Canine Parvovirus Type 2c in Taiwan. International Journal of Molecular Sciences, 18(12), 2703. https://doi.org/10.3390/ijms18122703

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