key: cord-0002398-mmd8ru3r authors: Chiang, Shu-Yun; Wu, Hung-Yi; Chiou, Ming-Tang; Chang, Min-Chen; Lin, Chao-Nan title: Identification of a novel canine parvovirus type 2c in Taiwan date: 2016-09-23 journal: Virol J DOI: 10.1186/s12985-016-0620-5 sha: b3d7fd02de63c9e801c668f9784a1f31be59cc97 doc_id: 2398 cord_uid: mmd8ru3r BACKGROUND: Taiwan has been considered free from canine parvovirus type 2c (CPV-2c) based on the last report of canine parvovirus type 2 (CPV-2) surveillance. However, since January 2015, the first report of CPV-2c in a puppy has occurred in Taiwan. There is currently limited information about the CPV-2c variant in Taiwan. In the present study, we characterized the previously unidentified CPV-2c variant and investigated the distribution of CPV-2 variants in Taiwan. METHODS: During January 2014 to April 2016, fecal or rectal swab samples from 99 dogs with suspected CPV-2 infection in Taiwan were collected. Eighty-eight were identified as being either CPV-2a, −2b or -2c variants positive by real-time PCR and sequence analysis. RESULTS: Sequence analysis of the 88 isolates confirmed CPV-2c as the dominant variant (54.6 %), followed by CPV-2b (26.1 %) and CPV-2a (19.3 %). Phylogenetic analysis demonstrated that the recent CPV-2c variants are similar to the Chinese CPV-2c strain but can be considered as novel Asian CPV-2c isolates. CONCLUSION: The present study provides evidence for the existence of a novel CPV-2c variant in Taiwan. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12985-016-0620-5) contains supplementary material, which is available to authorized users. puppy in Taiwan occurred, and there is to date limited information about the CPV-2c variant in Taiwan. In the present study, we examined this CPV-2c variant and investigated the distribution of CPV-2 variants in Taiwan. Clinical specimens (feces and/or rectal swab) were collected from 99 dogs with suspected CPV-2 infection from northern, central, southern, and eastern Taiwan from January 2014 to April 2016. These samples were mainly acquired from dogs with diarrhea and/or bloody diarrhea. The year of sampling and the age, clinical history, and CPV-2 types of the sampled dogs are summarized in Additional file 1. Viral DNA was extracted from the clinical samples (either feces or rectal swab) and screened for CPV-2 by real-time PCR, as described by Lin et al. [43] . Samples showing positive results for either type of specimen were included in this study. The partial VP2 gene of CPV-2 was amplified by PCR, as described by Buonavoglia et al. [1] , and the DNA fragments were purified and sequenced as described by Lin et al. [41] . The VP2 DNA sequences of our samples were compared to those of reference FPV (M38246), CPV-2 (M38245), CPV-2a (M24003), CPV-2b (M74849), new CPV-2a (JX 048605), new CPV-2b (JX048607), and CPV-2c (JF41 4818, JF414820, JF414822, FJ005247, GU380303, GU38 0305, KR611522, KT074339, KT162005, KF149962, FJ005196, GQ865518, FJ222821, FJ005213, FJ005238, FJ005214, KC196099, KM457119, AB120727, KP071956, KR559893, FJ005235). Multiple alignments of the nucleic acid and amino acid sequences were performed using the Clustal W method and the MegAlign program (DNASTAR, Madison, WI, USA). Phylogenetic analyses were performed with the maximum likelihood method using MEGA 6, version 6.06. Polymerase chain reaction (PCR) amplification and genotype analysis A total of 88 samples from 99 dogs were positive for CPV-2. Of the 88 CPV-2 isolates, 17, 23, and 48 were identified as CPV-2a (19.3 %), CPV-2b (26.1 %), and CPV-2c (54.6 %), respectively. First, we observed the CPV-2c variant collected in January 2015 (Additional file 1, Fig. 1 ). This variant rapidly spread throughout the Taiwanese dog population, and detection rates were 53.1 % (26/49) and 68.8 % (22/32) in 2015 (January to December) and 2016 (January to April), respectively ( Fig. 1 ). Taken together, our data revealed co-circulation of CPV-2a, −2b, and -2c on the island, and CPV-2c was identified as the current dominant variant. Partial VP2 nucleotide sequences were analyzed using DNASTAR software, revealing 96.6-100 %, 98.3-100 %, and 98.5-100 % homology within local CPV-2a isolates, CPV-2b isolates, and CPV-2c isolates, respectively ( Table 1 ). The same nucleotide sequences for local CPV-2c isolates exhibited 96.9-99 % and 97.5-99.2 % homology with those of the CPV-2a and -2b isolates, respectively (Table 1 ). In contrast to the low nucleotide sequence similarity between prototype CPV-2c (FJ222821, from Italy) and CPV-2c from Taiwan (97.9-99.0 %), the homology levels between our analyzed CPV-2c isolates (97.7-100 %) and Chinese CPV-2c (GU380303, GU380305, KR611522, KT074339, KT162005) appeared to be much higher (Table 1) . Comparisons among the 88 isolates and 12 reference strains at amino acids 267-440 are presented in Table 2 . For the partial VP2 analysis, 17 sequences with an Asn at position 426 were classified as CPV-2a. All of the CPV-2a strains show substitution at position 324 (Tyr to Ile) caused by mutation of TAT to ATT at nucleotide positions 970-972 of the VP2 gene (Table 2) . Interestingly, five CPV-2a strains from eastern Taiwan are identical in nucleotide sequence but have distinctive residues (Phe267Tyr, Tyr324Ile, and Thr440Ala) ( Table 2) . Four CPV-2b strains are identical to the prototype of CPV-2b ( Table 2 ). Fifteen of the 19 CPV-2b strains show substitution at position 267 (Phe to Tyr) and 324 (Tyr to Ile), and four CPV-2b strains only show substitution at position 324 (Tyr to Ile). Glu426, which is unique to strain CPV-2c, was first observed in 48 samples in this study. The phylogenetic tree of the partial VP2 gene for the 88 isolates and 28 reference strains was generated using the maximum likelihood method with MEGA 6, version 6.06. Two clusters based on the phylogenetic relationship of the partial VP2 gene of CPV-2c were detected (Fig. 3 ). The first cluster comprises prototype isolates, European isolates, American isolates, and Asian isolates of CPV-2c. The second cluster consists of all 48 Taiwanese isolates from this study and 4 Chinese isolates of CPV-2c (Fig. 3) . To date, these novel CPV-2c isolates have only been detected in China and Taiwan. Our findings demonstrated that the recent CPV-2c isolates in Taiwan are more genetically similar within the VP2 gene to Chinese CPV-2c strains than to prototype strains of CPV-2c. This is the first study to demonstrate the CPV-2c variant in Taiwan. Previous studies have shown that both the CPV-2a and -2b genotypes constitute the prevalent CPV-2 field strains circulating in Taiwan, and no CPV-2c cases were reported in the last two decades (Table 3 ). In our continuous surveillance and sequence analysis, Taiwan was considered to be free of CPV-2c. However, since January 2015, one case of CPV-2c occurred in a puppy in Taiwan. The CPV-2c variant appears thus far to have the highest detection rate in the dog population of Taiwan. Similar results have shown that CPV-2c replaced the previous circulation of CPV-2 strains in Italy [8] , Uruguay [19] , Argentina [21] , Brazil [20] , and the United States [44] . Interestingly, comparative VP2 genome analysis of the isolated CPV-2c reference revealed that the partial VP2 genome sequence of the Taiwanese strain is similar (97.7-100 %) to that of Chinese CPV-2c strains. Our results indicate that the recent CPV-2c isolate from Taiwan shares a common evolutionary origin with the Chinese CPV-2c strains and should be classified as novel Asian CPV-2c isolates. According to genotype Fig. 3 Phylogenetic analysis of nucleotide sequences of the partial VP2 genes of CPV-2 strains (522 bp). The phylogenetic tree was constructed using the maximum likelihood method with bootstrap analysis (n = 1,000) to determine the best fitting tree. CPV variants are indicated by □, •, ○, and ▲ for CPV-2, 2a, 2b, and 2c, respectively. The gray diagram represents the Taiwanese CPV-2 variants in this study surveillance between Taiwan [30, 31, 39, 41] and China [26-28, 37, 45-48] , CPV-2c was first detected in China in 2009 [28] . In contrast, no CPV-2c variant was observed in studies conducted over the last two decades in Taiwan (Table 3 ). In addition, the Taiwanese CPV-2c variant is more closely related to the Chinese CPV-2c strains than the recent Taiwanese CPV-2a (JX048605) and -2b (JX048607) isolates (Fig. 3) . Taken together, our results suggest that the Taiwanese CPV-2c may be present due to import from China at some time between the end of 2014 to early 2015 rather than to evolution of existing CPV-2a or -2b genotypes. An ongoing investigation and complete VP2 genome sequence analysis is needed to trace the genetic evolution of this novel CPV-2c variant. The recent Taiwanese CPV-2a is composed of two divergent lineages that have different ancestors. Most Taiwanese CPV-2a strains belong to the recent Taiwanese lineage of CPV-2a, sharing a common amino acid substitution (Tyr324Ile). The second lineage of the Taiwanese CPV-2a variant is more closely related to recent Uruguayan and Chinese CPV-2a strains and has distinctive amino acid substitutions of Phe267Tyr, Tyr324Ile, and Thr440Ala. This new CPV-2a variant was discovered in China and Uruguay between 2006 and 2009 and in 2010 [49] , respectively. This is the first detection of this lineage of CPV-2a in eastern Taiwan in 2015. However, this new CPV-2a variant recently emerged in Uruguay and underwent clonal expansion [49] . An ongoing investigation is aimed at determining whether this new CPV-2a variant will replaced the CPV-2c variant in the Taiwanese dog population. Amino acid substitution of Tyr324Ile has been observed in Korea [32, 34] , China [26, 27, 37, [45] [46] [47] [48] , Thailand [36] , Uruguay [49, 50] , Japan [38] , Taiwan [39, 41] , and India [40, 42] . Interestingly, the frequency of the Ile324 variant has reached a high prevalence among Taiwanese CPV-2 isolates (94.5 %), and our results revealed that this variant is not only present among Taiwanese CPV-2a strains but also CPV-2b and -2c strains. In addition, this study reports for the first time the amino acid substitution of Phe267Tyr in Taiwan. Surprisingly, all of the Tyr267 variants among Taiwanese CPV-2a, −2b, and −2c strains contain the amino acid substitution of Tyr324Ile. Our review of sequence analysis in the literature indicated that this phenomenon is also found in Uruguay [51] . The functions of CPV-2 residues 267 and 324 are still unknown and remain to be elucidated. The substitution of Gln370Arg is unique to the Taiwanese CPV-2c strains, and this mutation is also observed in Chinese panda parvovirus [52] and Chinese CPV-2c strains [26, 27] . Residue 359 and 375 constitute a flexible surface loop of the capsid protein that is adjacent to a double Ca 2+ -binding site; this region is essential for virus infectivity, and changes are correlated with the ability of the virus to cause erythrocyte hemagglutination [53] . Therefore, it remains to be investigated whether Gln370Arg substitution causes antigenic alterations. Mutation of residue 420 had been reported in Brazilian reference CPV-2c strains (Phe420Leu) [54] . In the present study, we detected 10 CPV-2c strains with a unique change at the same position, yet Phe420Ser was unique in these 10 Taiwanese CPV-2c strains. Therefore, further studies focusing on potential variants of the CPV-2c strains should be conducted to elucidate the relationship between Phe420Ser substitution and viral pathogenicity. Although several studies have demonstrated the efficacy of the current CPV-2 vaccine against CPV-2c infection [55, 56] , some evidence suggests that dogs with the complete vaccination program still suffer from CPV-2c [6] . In the present study, four of 22 CPV-2c-diseased dogs died despite vaccination (C104-030, C104-031, C104-042, C104-216) (Additional file 1). Among those that died, three were under 6 months of age. Surprisingly, despite having undergone the complete vaccination program, one adult dog (strain no. C104-216) was infected by this novel CPV-2c variant. Therefore, coinfection with other diseases needs to examined, and the efficacy of the current vaccine against this novel CPV-2c variant remains to be evaluated, especially in regard to the amino acid substitutions observed in this novel CPV-2c variant compared to the CPV-2c prototype. This is the first report to identify a novel CPV-2c variant in Taiwan. The novel CPV-2c variant was found to be Evidence for evolution of canine parvovirus type 2 in Italy Evolution of CPV-2 and implication for antigenic/genetic characterization Fine mapping of canine parvovirus B cell epitopes Multiple amino acids in the capsid structure of canine parvovirus coordinately determine the canine host range and specific antigenic and hemagglutination properties High local genetic diversity of canine parvovirus from Ecuador Evidence for immunisation failure in vaccinated adult dogs infected with canine parvovirus type 2c The study molecular epidemiology of canine parvovirus Characterisation of the canine parvovirus type 2 variants using minor groove binder probe technology Canine parvovirus: a review of epidemiological and diagnostic aspects, with emphasis on type 2c Western European epidemiological survey for parvovirus and coronavirus infections in dogs Epidemiological evolution of canine parvovirus in the Portuguese domestic dog population First detection of canine parvovirus type 2c in pups with haemorrhagic enteritis in Spain Characterization of Canine parvovirus 2 variants circulating in Greece A Ten-Year Molecular Survey on Parvoviruses Infecting Carnivores in Bulgaria Canine parvovirus type 2c identified from an outbreak of severe gastroenteritis in a litter in Sweden Characterization of the partial VP2 gene region of canine parvoviruses in domestic cats from Turkey Molecular characterization of canine parvovirus-2 variants circulating in Tunisia Occurrence of canine parvovirus type 2c in the United States First detection of canine parvovirus type 2c in South America Typing of canine parvovirus strains circulating in Brazil between Evolution of canine parvovirus in Argentina between years 2003 and 2010: CPV2c has become the predominant variant affecting the domestic dog population Genotyping of Canine parvovirus in western Mexico Molecular epidemiology of canine parvovirus in Morocco A novel antigenic variant of Canine parvovirus from a Vietnamese dog Occurrence of canine parvovirus type 2c in the dogs with haemorrhagic enteritis in India 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 Typing of Canine Parvovirus Strains Circulating in North-East China The first detection of canine parvovirus type 2c in China Chronological analysis of canine parvovirus type 2 isolates in Japan Antigenic types of canine parvoviruses prevailing in Taiwan Phylogenetic analysis of canine parvovirus VP2 gene in Taiwan Genetic analysis of VP2 gene of canine parvovirus isolates in Korea Prevalence and genetic characterization of canine parvoviruses in Korea Molecular insights into the phylogeny of canine parvovirus 2 (CPV-2) with emphasis on Korean isolates: a Bayesian approach Isolation, molecular characterization and phylogenetic analysis of canine parvovirus Sequence analysis of VP2 gene of canine parvovirus isolates in Thailand Phylogenetic analysis of the VP2 gene of canine parvoviruses circulating in China Analysis of the VP2 protein gene of canine parvovirus strains from affected dogs in Japan Genotyping of canine parvovirus type 2 VP2 gene in southern Taiwan in 2011 Phylogenetic analysis of canine parvovirus partial VP2 gene in India Genetic characterization of type 2a canine parvoviruses from Taiwan reveals the emergence of an Ile324 mutation in VP2 Molecular typing of canine parvovirus strains circulating from 2008-2012 in an organized kennel in India reveals the possibility of vaccination failure Development of SYBR Green-Based Real-Time PCR for the Detection of Canine, Feline and Porcine Parvoviruses Canine parvovirus types 2c and 2b circulating in North American dogs in 2006 and Phylogenetic Analysis of Canine Parvovirus VP2 Gene in China 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 Full-length genomic characterizations of two canine parvoviruses prevalent in Northwest China Phylogenetic analysis of canine parvovirus isolates from Sichuan and Gansu provinces of China in 2011 Recent spreading of a divergent canine parvovirus type 2a (CPV-2a) strain in a CPV-2c homogenous population Phylogenetic and genome-wide deep-sequencing analyses of canine parvovirus reveal co-infection with field variants and emergence of a recent recombinant strain Phylodynamics analysis of canine parvovirus in Uruguay: evidence of two successive invasions by different variants Identification of canine parvovirus with the Q370R point mutation in the VP2 gene from a giant panda Host range and variability of calcium binding by surface loops in the capsids of canine and feline parvoviruses Genetic analysis of canine parvovirus type 2c Do two current canine parvovirus type 2 and 2b vaccines provide protection against the new type 2c variant? Canine parvovirus type 2 vaccine protects against virulent challenge with type 2c virus The authors would like to thank Prof. Ling-Ling Chueh and Prof. Hsin-Fu Liu for comments that greatly improved the manuscript. No funding was obtained for this study.Availability of data and materials Not applicable. Authors' contributions SYC analyzed the experimental data and wrote the manuscript. HYW and MTC managed the study, provided materials and reagents, contributed to the interpretation of the data, and co-wrote the manuscript. MCC contributed to the DNA extraction and PCR. CNL designed and analyzed the experimental data and wrote the manuscript. All of the authors read and approved the final manuscript. The authors declare that they have no competing interests. Not applicable. The study did not involve any animal experiment. The Institutional Animal Care and Use Committee (IACUC) of National Pingtung University of Science and Technology did not deem it necessary for this research group to obtain formal approval to conduct this study. Submit your next manuscript to BioMed Central and we will help you at every step: