- 1Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, China
- 2State Key Laboratory of Agrobiotechnology, Department of Biochemistry and Molecular Biology, College of Biological Sciences, China Agricultural University, Beijing, China
- 3College of Veterinary Medicine, Hebei Agricultural University, Baoding, China
Canine parvovirus type 2 (CPV-2) and feline panleukopenia virus (FPV) cause severe disease in young animals, pups, and kittens. CPV-2 evolved from FPV by altering the species-specific binding of the viral capsid to the host receptor, i.e., the transferrin receptor (TfR), and CPV-2 genetic variants have been identified by specific VP2 amino acid residues (297, 426). Early studies focused on the main capsid protein VP2; however, there have been limited studies on the non-structural protein NS1. In this study, we identified the genetic variants of clinical samples in dogs and cats in northern China during 2019–2020. The genetic characterization and phylogenetic analyses of VP2 and NS1 gene were also conducted. The results revealed that the CPV-2c was identified as the major genetic variant. One new CPV-2b and two CPV-2c strains were collected from cats. Four mutation sites (60, 630, 443, and 545 amino acid residues) were located in the functional domains of the NS1 protein. The phylogenetic analysis of VP2 and NS1 genes showed that they were clustered by geographical regions and genotypes. The gene mutation rate of CPV-2 was increasing in recent years, resulting in a complex pattern of gene evolution in terms of host preference, geographical selection, and new genetic variants. This study emphasizes that continuous molecular epidemiological surveillance is required to understand the genetic diversity of FPV and CPV-2 strains.
Introduction
Canine parvovirus type 2 (CPV-2) and feline panleukopenia virus (FPV) infections seriously threaten the health of dogs and cats, especially puppies and kittens, with high morbidity and mortality rates (1, 2). FPV was identified at the beginning of the 20th century and has maintained a certain degree of genetic stability (3). CPV-2 was first reported in the USA and Australia in 1978 (4). The emergence of CPV-2 resulted from site-specific mutations on the surface of the viral capsid, and the cross-species transmission of CPV-2 occurs by altering the species-specific binding of the viral capsid to the host receptor, i.e., the transferrin receptor TfR (5, 6). Although CPV-2 is a DNA virus, the rate of amino acid mutations almost reaches that of RNA viruses (7). After its emergence, a new antigenic variant (CPV-2a) emerged with five important amino acid mutations (M87L, I101T, A300G, D305Y, and V555I) in the VP2 gene from 1979 to 1981. These mutations affected the antigenic reactivity of monoclonal antibodies, the binding affinity to TfR, and the ability to replicate in cats (8). In 1984, another antigenic variant (CPV-2b) emerged, with two amino acid mutations (N426D and I555V). A third genetic variant (CPV-2c), with 426E, was reported in 2000. At the end of the 20th century, the new CPV-2a and new CPV-2b genetic variants were identified, carrying an S297A substitution (9). In the 2000s, the CPV-2c genetic variant was frequently found in Europe and South America but was relatively rare in Asia, where CPV-2a was more predominant (10). However, currently, the CPV-2c genetic variant is distributed worldwide including in Asia, Africa, Australia, Europe, North and South America, China, Vietnam, Mongolia, and Thailand (11). In China, the molecular surveillance carried out between 1986 and 2017 showed a clear predominance of CPV-2a. However, CPV-2c strains had become dominant in 2018 and 2021. One CPV-2c strain with typical amino acid residues (5G, 267Y, 324I, and 370R) in the VP2 protein has been recognized as the “Asian CPV-2c” genetic variant, which had never been reported prior to 2013 (12). The reasons for the emergence of different antigenic variants and their spread among different hosts have been attributed to amino acid mutation sites (13, 14).
CPV-2 and FPV are non-enveloped single-stranded linear DNA viruses of the Parvoviridae family (15), with a genome of ~5,200 nt in length (16). The genome contains two open reading frames (ORFs), which encode for two non-structural (NS1 and NS2) and two structural (VP1 and VP2) proteins (17). VP2 protein, the main component of the viral capsid, not only induces the neutralizing antibodies by acting as the main protective antigen, but also determines hemagglutination, cellular tropism, and the host ranges of the virus (18). Additionally, the NS1 protein plays an important role in virus replication, DNA packaging, cytotoxicity, and pathogenicity (18–20). In this study, we aimed to perform molecular and phylogenetic analyses of the VP2 and NS1 genes from FPV and CPV-2 variants in northern China.
Materials and methods
Clinical samples
A total of 84 samples (feces or rectal swabs) were collected from domestic dogs (n = 51) and cats (n = 33) with clinical signs, like fever, vomiting, and diarrhea. Animals were initially identified to be infected with CPV-2 or FPV using colloidal gold test strips (BioNote Rapid Test Kit; BioNote, Gyeonggi-do, Korea). The samples were collected from different areas of northern China: Xiong'an (n = 31), Jinan (n = 8), Langfang (n = 1), Haidian District in Beijing (n = 4), and Tongzhou District in Beijing (n = 40). Viral DNA was extracted by the Aidlab DNA rapid extraction kit (Aidlab Biotech, Beijing, China), according to the manufacturer's instructions, and stored at −20°C until further processing. The details are summarized in Table 1.
Table 1. Detailed description of canine parvovirus type 2 (CPV-2) and feline panleukopenia virus (FPV) characterized in this study.
PCR and DNA sequencing assays
The clinical positive samples were further confirmed by PCR, targeting the 573 bp and 286 bp VP2 fragments of CPV-2 and FPV, respectively (Table 2). PCR was performed using PrimeSTAR® Max DNA Polymerase (Takara Biotech, Dalian, China) according to the manufacturer's instructions, in a final volume of 20 μl comprising 10 μl of PrimeSTAR MAX Premix (2 × ), 1 μl of each primer (10 μM), 1 μl of DNA, and 7 μl of sterile distilled water. The PCR condition was set as an initial denaturation step at 98°C for 1 min; 35 cycles at 98°C for 10 s, annealing at 56°C for 5 s, extension at 72°C for 5 s; and a final extension at 72°C for 10 min. The positive samples were then subjected to the amplification of the complete NS1 and VP2 genes using two primer pairs NS1-F/R or VP2-F/R, respectively, in a final volume of 50 μl, comprising 25 μl of PrimeSTAR MAX Premix (2 × ), 1 μl of each primer (10 μM), 2 μl of DNA, and 21 μl of sterile distilled water. The PCR condition was set as an initial denaturation step at 98°C for 1 min; 35 cycles at 98°C for 10 s, annealing at 57°C for 15 s, extension at 72°C for 30 s; and a final extension at 72°C for 10 min. Positive amplicons were purified using a gel extraction kit (OMEGA Bio Tek, Winooski, VT, USA) and submitted to BGI biotech in Beijing for direct Sanger sequencing. The VP2 gene was sequenced using the primers pair used for amplification (VP2-F/R) and with three additional internal primers (CPV-VP2-F/R/F1). The NS1 gene was also sequenced using the primers pair used for amplification (NS1-F/R) and with three additional internal primers (CPV-NS1-658F/800R/1228F). Detailed information on the primers is shown in Table 2. Consensus sequences were assembled using the EditSeq program (DNAStar, Madison, WI, USA).
Sequence and phylogenetic analyses
The complete coding sequences of the NS1 and VP2 genes were assembled and translated into amino acid sequences (668 and 584 aa, respectively) using the EditSeq software. The translated VP2 and NS1 amino acid sequences were aligned with reference sequences (new CPV-2a, AY742935; new CPV-2a, EF011664; new CPV-2a, JN867611; new CPV-2b, EU659119; new CPV-2b, JQ268284; new CPV-2b, LC270892; CPV-2c, KM457125; CPV-2c, KU508407; CPV-2c, MK388674; CPV-2c, MN832850; FPV, KP019621; FPV, MT614366), which were obtained from the NCBI database. The sequence alignments were obtained using the ClustalW program included in the MegAlign software (DNAStar).
The phylogenetic relationships were evaluated using MEGA6.0 (21) and iTOL software (22). The phylogenetic trees were constructed using the neighbor-joining and maximum-likelihood methods. A bootstrap analysis with 1,000 replicates was performed to assess the confidence level of the branching pattern. For a more comprehensive evaluation, 41 sequences of the NS1 gene and 46 sequences of the VP2 gene were downloaded from GenBank as reference sequences (Supplementary Table 1).
Results
Antigenic typing of FPV and CPV-2 variants
A total of 84 clinical samples were confirmed to be positive for FPV or CPV-2. Sequencing of five VP2 and nine NS1 gene sequences was not successful, as shown in Table 1. Based on 79 VP2 amino acid sequences, 32 sequences from cats and 47 sequences from dogs were obtained. Among the 47 strains from dogs, 13, 1, and 33 strains were typed as new CPV-2a, new CPV-2b, and CPV-2c, respectively. Twenty-nine strains from the 32 cats were typed as FPV, and one strain (HD-cat-newCPV2b-9) was identified as new CPV-2b and two strains (HD-cat-CPV-2c-79 and TZ-cat-CPV-2c-121) were identified as CPV-2c.
Sequence analysis of VP2 gene
Sequence analysis showed that all CPV-2c strains had Arg at VP2-370 amino acid residue, except for XA-CPV-2c-0-8 which showed Gln. Meanwhile, new CPV-2a/2b strains showed Gln at the same amino acid residue. In addition, CPV-2c strains harbored 440Thr, and the new CPV-2a/2b strains contained 440Ala. Furthermore, other mutations also emerged among the CPV strains, i.e., P13S (TZ-CPV-164), E18D (XA-CPV-0-22), L218I (XA-CPV-0-14), and G441R (JN-CPV-10 and TZ-CPV-150). The VP2 gene of FPV strains was more conserved than that of CPV-2 strains. Indeed, 44.8% (13/29) of FPV strains showed the A91S mutation. Meanwhile, three other VP2 mutations were observed in FPV strains: F290Y (TZ-FPV-236), T425A (JN-FPV-12 and JN-FPV-90), and V535I (TZ-FPV-243).
Sequence analysis of NS1 gene
Based on 75 NS1 amino acid sequences, 26 sequences from cats and 49 sequences from dogs were obtained. The mutation rates of the non-synonymous amino acid mutation sites in different virus types and genetic variants (FPV/new CPV-2a/new CPV-2b/CPV-2c) are shown in Table 3.
Table 3. Statistics of mutation sites and rates in the NS1 gene from FPV and different CPV-2 genetic variants.
CPV-2 and FPV presented common mutation sites at amino acids 60, 222, 544, 545, 572, 624, and 667. Residues 247 and 545 were further analyzed due to their key role in protein function. All CPV-2 sequences showed 247Q, while 39.1% of FPV sequences showed 247R, and the remaining 60.9% of FPV sequences showed 247Q. Moreover, in this study, 8.7% of the FPV strains showed E545K, and 69.2% (36/52) of CPV-2 strains harbored E545V. Amino acid mutations 19, 187, 222, 356, 583, 596, 624, and 667 were newly discovered in this study, which awaited further investigation.
Phylogenetic analysis of VP2 gene
Figure 1 shows the phylogenetic tree inferred from the VP2 gene, which clustered into two main clades (Clades I and II; FPV and CPV-2). Fifteen FPV strains were phylogenetically closed to foreign strains and 13 FPV strains were closed to the strains of Asian origin. All the CPV-2 sequences were clustered in Clade II, and the branches were formed by genetic types. Clade II could be further divided into two parts: Clade IIA clustering new CPV-2a/2b strains and Clade IIB clustering CPV-2c strains. The phylogenetic analysis showed that the new CPV-2b strain infecting cats (HD-cat-newCPV2b-9) was close to another new CPV-2b strain from dogs (XA-newCPV-2b-150). Meanwhile, two CPV-2c strains infecting cats (HD-cat-CPV-2c-79 and TZ-cat-CPV-2c-121) also clustered with other CPV-2c strains from dogs.
Figure 1. Phylogenetic tree based on a comparison of the structural (VP2) gene of canine parvovirus type 2 (CPV-2) and feline panleukopenia virus (FPV). Phylogenetic relationships were calculated using the maximum likelihood (ML) method. Gaps were handled by pairwise deletion and bootstrap values were calculated from 1,000 replicates. Strains analyzed in this study are marked in red and reference strains are marked in black. CPV-2 strains isolated from cats in this study are marked with “.” The monophyletic clusters of new CPV-2a/2b are marked in blue (Clade IIA), the monophyletic clusters of CPV-2c are marked in yellow (Clade IIB), and the monophyletic clusters of FPV are marked in pink (Clade I).
Phylogenetic analysis of NS1 gene
Figure 2 shows the phylogenetic tree of the NS1 gene, and unlike the VP2 gene, the characteristic of clustering by genetic types was less supported (Figure 2). Seventeen FPV sequences formed the Clade I, new CPV-2a/2b, and CPV-2c strains clustered in Clades IIA and B, respectively. Meanwhile, two FPV strains and five CPV-2c strains were closed to Clade IIA, and another two FPV strains and two new CPV-2b strains were closed to Clade IIB. In addition, we found an individual clade comprising four strains with high homology (TZ-FPV-148, TZ-cat-CPV-2c-121, XA-CPV-2c-0-8, and TZ-FPV-99). The new CPV-2b and one CPV-2c strain infecting cats (HD-cat-newCPV2b-9 and HD-cat-CPV-2c-79) showed high homology with CPV-2 strains from dogs.
Figure 2. Phylogenetic tree based on a comparison of the non-structural (NS1) gene of canine parvovirus type 2 (CPV-2) and feline panleukopenia virus (FPV). Phylogenetic relationships were calculated using the maximum likelihood (ML) method. Gaps were handled by pairwise deletion and bootstrap values were calculated from 1,000 replicates. Strains analyzed in this study are marked in red and reference strains are marked in black. CPV-2 strains isolated from cats in this study are marked with “.” The monophyletic clusters of new CPV-2a/2b are marked in blue (Clade IIA), the monophyletic clusters of CPV-2c are marked in yellow (Clade IIB), and the monophyletic clusters of FPV are marked in pink (Clade I).
Discussion
The VP2 protein is the main component of the virus capsid for determining antigenicity, host range, and hemagglutination (9). It has five important domains, Loop 1 (aa 50–100), Loop 2 (aa 200–250), Loop 3 (aa 300–350), Loop 4 (aa 400–450), and Loop 5 (aa 350–400) (23). Among them, different Loops perform specific functions. Loop 1, Loop 2, and Loop 4 not only constitute the spike of the triple folding (24) but also form the antigen A site with four key amino acid sites (93, 222, 224, and 426). Residue 440 is located on the top of the 3-fold spike as a main viral genetic site, which might result in the emergence of further genetic variants (25). Loop 3 not only forms the shoulder of the triple folding but also participates in the formation of the antigen B site, with the key amino acid sites 299, 300, and 302 (26).
In the results of CPV-2 VP2 gene alignment, CPV-2c strains harbored the key mutation sites with 370R and 440T, while new CPV-2a/2b genetic variants showed 370Q and 440A. Geng and colleagues reported that the VP2 protein of CPV-2c strains exhibited a unique Q370R substitution (27). Zhuang et al. also discovered that mutation Q370R was specific for all identified CPV-2c strains, and all the new CPV-2a strains showed the T440A mutation (28). This characteristic might help to distinguish CPV-2 genetic variants in more detail. In this study, CPV-2c was identified as the primary genetic variant in Beijing and Xiong'an of northern China. CPV-2c was first reported in Jilin province in 2010, and after that, it continued to spread from the north to the south of China (29) and it was reported as the dominant strain in many regions (30). Among FPV strains, 44.8% of strains had the A91S mutation. A91S is located in Loop 1 and T425A is located in Loop 4. These two mutation sites could influence the recognition of neutralizing antibodies against the antigen A site, maybe resulting in evading host immune surveillance. Based on the 3D structure prediction, A91S mutation extended the random coil of aa residues from 92–95 to 91–95 in FPV VP2 protein. Therefore, it is possible that the A91S mutation may also affect the receptor-binding ability (31). Further studies are then necessary to analyze in vitro any biological relevance of these mutations observed in silico.
Relatively enough data are available on the VP2 gene. But unfortunately, the current research on the NS1 gene is not enough to understand the corresponding protein structure and function, only three or five amino acid mutation sites were found. In previous studies, researchers evaluated the function of NS1 protein and described the potential location of functional domains, the origin of replication (ORI) binding (16–275aa), helicase (299–486aa), and transactivation (600–667aa) functional domains (32). In this study, we highlighted the role of amino acid positions 60 and 630, which are located on the ORI binding and transactivation functional domains, respectively. Transactivation can increase the rate of gene expression (33). Residue 443 putatively lies in the β3-sheet of the Walker motif B of the helicase domain, and residues 350 and 544–545 are located between the α5- and α6-helices close to the α11-helix of the same domain (32). Helicases open the complementary double strands of DNA using energy derived from ATP hydrolysis to obtain single strands; therefore, these functional domains play an important role in early viral replication and DNA assembly (34). Nevertheless, most mutation sites tend to gather in these functional domains, which might influence the host range and virulence. Researchers considered that amino acid site 248 can be used to distinguish CPV-2 and FPV (11). In the case of a nucleotide change in the second base of the codon (c743t), FPV showed 248T (act), while CPV-2 showed 248I (att). By contrast, our research found that all CPV-2 strains showed 248I (att), while among FPV strains, 78.3% of strains showed 248T (act) and five strains (TZ-FPV-99, TZ-FPV-108, TZ-FPV-143, TZ-FPV-148, and TZ-FPV-193) showed 248I (att). This indicated that neither NS1-248 aa residue distinguished FPV from CPV-2.
At present, although various hypotheses explaining its derivation and sudden emergence have been proposed, opinions on the origin of CPV-2 were still not conclusive. The most widely accepted hypothesis for its emergence is that CPV-2 was derived from FPV in cats or FPV-like viruses in wild animals by a natural genetic mutation (35). Truyen et al. reported that ~5% of isolates from domestic cats in Germany and the USA were either CPV-2a or CPV-2b, while the original CPV-2 could not replicate in feline tissues (36). The pathogenicity of the CPV-2 genetic variants was milder and could replicate only poorly in cats, with mild clinical signs, such as leukopenia and poor growth. Subsequently, it was reported that the genetic variants of CPV-2,−2a,−2b, and−2c have acquired the feline host range, and were able to infect and replicated efficiently in cats and some wildlife (37–39). Our results revealed that CPV-2 genetic variants were circulating in dogs and, with low prevalence in cats, which demonstrated the need for continued surveillance of FPV and CPV-2 infections considering the key amino acid sites in both structural and non-structural proteins.
Conclusion
The epidemiology of CPV-2 and FPV was assessed in northern China from 2019 to 2020, and CPV-2c strains were the dominant CPV-2 genetic variant among dogs. Within the summarized mutation sites of NS1, amino acid positions 60, 630, 443, 544, and 545 maybe play an important role in genomic replication and virion packaging. Further studies are then necessary to elucidate the relevance of these mutations in the biological properties of the NS1 protein functional domains. Besides, one new CPV-2b strain and two CPV-2c strains were collected from domestic cats.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The animal study was reviewed and approved by the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences.
Author contributions
XC: resources. SL: investigation and writing—original draft. YH: formal analysis. TQ: supervision and project administration. GZ, YL, JW, and WL: writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Agricultural Science and Technology Innovation Program (Grant No: ASTIP-IAS15).
Acknowledgments
We thank Prof. Shangjin Cui for his valuable discussion and guidance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.934849/full#supplementary-material
References
1. Battilani M, Balboni A, Ustulin M, Giunti M, Scagliarini A, Prosperi S. Genetic complexity and multiple infections with more parvovirus species in naturally infected cats. Vet Res. (2011) 42:43. doi: 10.1186/1297-9716-42-43
2. Duque-Garcia Y, Echeverri-Zuluaga M, Trejos-Suarez J, Ruiz-Saenz J. Prevalence and molecular epidemiology of canine parvovirus 2 in diarrheic dogs in Colombia, South America: a possible new Cpv-2a is emerging? Vet Microbiol. (2017) 201:56–61. doi: 10.1016/j.vetmic.2016.12.039
3. Decaro N, Desario C, Miccolupo A, Campolo M, Parisi A, Martella V, et al. Genetic analysis of feline panleukopenia viruses from cats with gastroenteritis. J Gen Virol. (2008) 89:2290–8. doi: 10.1099/vir.0.2008/001503-0
4. Fagbohun OA, Omobowale TO. Sequence and phylogenetic analysis of canine parvovirus-2 isolates in dogs revealed circulation of three subtypes In Nigeria. Virusdisease. (2018) 29:411–5. doi: 10.1007/s13337-018-0475-z
5. Nguyen Manh T, Piewbang C, Rungsipipat A, Techangamsuwan S. Molecular and phylogenetic analysis of vietnamese canine parvovirus 2c originated from dogs reveals a new Asia-Iv clade. Transbound Emerg Dis. (2020) 68:1445–53. doi: 10.1111/tbed.13811
6. Cureton DK, Harbison CE, Cocucci E, Parrish CR, Kirchhausen T. Limited transferrin receptor clustering allows rapid diffusion of canine parvovirus into clathrin endocytic structures. J Virol. (2012) 86:5330–40. doi: 10.1128/JVI.07194-11
7. Voorhees IEH, Lee H, Allison AB, Lopez-Astacio R, Goodman LB, Oyesola OO, et al. Limited intrahost diversity and background evolution accompany 40 years of canine parvovirus host adaptation and spread. J Virol. (2019) 94:e01162-19. doi: 10.1128/JVI.01162-19
8. Miranda C, Thompson G. Canine parvovirus: the worldwide occurrence of antigenic variants. J Gen Virol. (2016) 97:2043–57. doi: 10.1099/jgv.0.000540
9. Li G, Ji S, Zhai X, Zhang Y, Liu J, Zhu M, et al. Evolutionary and genetic analysis of the Vp2 gene of canine parvovirus. BMC Genom. (2017) 18:534. doi: 10.1186/s12864-017-3935-8
10. Li C, Tang J, Chen Z, Niu G, Liu G. A divergent canine parvovirus type 2c (Cpv-2c) isolate circulating in China. Infect Genet Evol. (2019) 73:242–7. doi: 10.1016/j.meegid.2019.05.004
11. Mira F, Canuti M, Purpari G, Cannella V, Di Bella S, Occhiogrosso L, et al. Molecular characterization and evolutionary analyses of carnivore protoparvovirus 1 Ns1 gene. Viruses. (2019) 11:308. doi: 10.3390/v11040308
12. Liu C, Gao J, Li H, Sun F, Liang H, Liu H, et al. Phylogenetic characteristics of canine parvovirus type 2c variant endemic in Shanghai, China. Viruses. (2021) 13:2257. doi: 10.3390/v13112257
13. Spera CG, Lorenzetti E, Lavorente FLP, De Calasans Marques G, Bisca JM, Teixeira CR, et al. Canine parvovirus 2b in fecal samples of asymptomatic free-living south American coatis (Nasua Nasua, Linnaeus, 1766). Braz J Microbiol. (2020) 51:1399–403. doi: 10.1007/s42770-020-00293-2
14. He HJ, Zhang W, Liang J, Lu M, Wang R, Li G, et al. Etiology and genetic evolution of canine coronavirus circulating in five provinces of China, during 2018–2019. Microb Pathog. (2020) 145:104209. doi: 10.1016/j.micpath.2020.104209
15. Moon BY, Jang J, Kim SH, Kim YH, Lee HK, So B, et al. Genetic characterization of canine parvovirus type 2c from domestic dogs in Korea. Transbound Emerg Dis. (2020) 67:1645–53. doi: 10.1111/tbed.13501
16. Chen MR, Guo XY, Wang ZY, Jiang YT, Yuan WF, Xin T, et al. Isolation and sequence analysis of the complete Vp2 gene of canine parvovirus from chinese domestic pets and determination of the pathogenesis of these circulating strains in Beagles. Pol J Vet Sci. (2019) 22:287–96. doi: 10.24425/pjvs.2019.129219
17. Mira F, Purpari G, Di Bella S, Colaianni ML, Schiro G, Chiaramonte G, et al. Spreading of canine parvovirus type 2c mutants of asian origin in Southern Italy. Transbound Emerg Dis. (2019) 66:2297–304. doi: 10.1111/tbed.13283
18. Feng H, Hu GQ, Wang HL, Liang M, Liang H, Guo H, et al. Canine parvovirus Vp2 protein expressed in silkworm pupae self-assembles into virus-like particles with high immunogenicity. PLoS ONE. (2014) 9:E79575. doi: 10.1371/journal.pone.0079575
19. Wang D, Yuan W, Davis I, Parrish CR. Nonstructural protein-2 and the replication of canine parvovirus. Virology. (1998) 240:273–81. doi: 10.1006/viro.1997.8946
20. Wang X, Zhang J, Huo S, Zhang Y, Wu F, Cui D, et al. Development of a monoclonal antibody against canine parvovirus Ns1 protein and investigation of Ns1 dynamics and localization in Cpv-infected cells. Protein Expr Purif. (2020) 174:105682. doi: 10.1016/j.pep.2020.105682
21. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. Mega6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. (2013) 30:2725–9. doi: 10.1093/molbev/mst197
22. Letunic I, Bork P. Interactive tree of life (Itol) V5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. (2021) 49:W293–6. doi: 10.1093/nar/gkab301
23. Llamas-Saiz AL, Agbandje-Mckenna M, Parker JS, Wahid AT, Parrish CR, Rossmann MG. Structural analysis of a mutation in canine parvovirus which controls antigenicity and host range. Virology. (1996) 225:65–71. doi: 10.1006/viro.1996.0575
24. Allison AB, Organtini LJ, Zhang S, Hafenstein SL, Holmes EC, Parrish CR. Single mutations in the Vp2 300 loop region of the three-fold spike of the carnivore parvovirus capsid can determine host range. J Virol. (2016) 90:753–67. doi: 10.1128/JVI.02636-15
25. Li C, Tang J, Chen Z, Li Q, Huang Z, Wang Q, et al. Genetic characterization of the complete genome of a mutant Canine parvovirus isolated in China. Arch Virol. (2018) 163:521–5. doi: 10.1007/s00705-017-3586-8
26. Casal JI. Use of parvovirus-like particles for vaccination and induction of multiple immune responses. Biotechnol Appl Biochem. (1999) 29:141–50.
27. Geng Y, Guo D, Li C, Wang E, Wei S, Wang Z, 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. doi: 10.1371/journal.pone.0137288
28. Zhuang QY, Qiu Y, Pan ZH, Wang SC, Wang B, Wu WK, et al. Genome sequence characterization of canine parvoviruses prevalent in the sichuan province Of China. Transbound Emerg Dis. (2019) 66:897–907. doi: 10.1111/tbed.13100
29. Jiang H, Yu Y, Yang R, Zhang S, Wang D, Jiang Y, et al. Detection and molecular epidemiology of canine parvovirus type 2 (Cpv-2) circulating in Jilin Province, Northeast China. Comp Immunol Microbiol Infect Dis. (2021) 74:101602. doi: 10.1016/j.cimid.2020.101602
30. Mittal M, Chakravarti S, Mohapatra JK, Chug PK, Dubey R, Upmanuyu V, et al. Molecular typing of canine parvovirus strains circulating from 2008 To 2012 in an organized kennel in india reveals the possibility of vaccination failure. Infect Genet Evol. (2014) 23:1–6. doi: 10.1016/j.meegid.2014.01.015
31. Chen X, Wang J, Zhou Y, Yue H, Zhou N, Tang C. Circulation of heterogeneous carnivore protoparvovirus 1 in diarrheal cats and prevalence of an a91s feline panleukopenia virus variant In China. Transbound Emerg Dis. (2022) 69:e2913–25. doi: 10.1111/tbed.14641
32. Niskanen EA, Ihalainen TO, Kalliolinna O, Hakkinen MM, Vihinen-Ranta M. Effect of atp binding and hydrolysis on dynamics of canine parvovirus Ns1. J Virol. (2010) 84:5391–403. doi: 10.1128/JVI.02221-09
33. Storgaard T, Christensen J, Aasted B, Alexandersen S. Cis-acting sequences in the aleutian mink disease parvovirus late promoter important for transcription: comparison to the canine parvovirus and minute virus of mice. J Virol. (1993) 67:1887–95. doi: 10.1128/jvi.67.4.1887-1895.1993
34. Strahan R, Dabral P, Dingman K, Stadler C, Hiura K, Verma SC. Kaposi's sarcoma-associated herpesvirus deregulates host cellular replication during lytic reactivation by disrupting the mcm complex through Orf59. J Virol. (2018) 92:e00739-18. doi: 10.1128/JVI.00739-18
35. Truyen U, Gruenberg A, Chang SF, Obermaier B, Veijalainen P, Parrish CR. Evolution of the feline-subgroup parvoviruses and the control of canine host range in vivo. J Virol. (1995) 69:4702–10. doi: 10.1128/jvi.69.8.4702-4710.1995
36. Truyen U, Evermann JF, Vieler E, Parrish CR. Evolution of canine parvovirus involved loss and gain of feline host range. Virology. (1996) 215:186–9. doi: 10.1006/viro.1996.0021
37. Charoenkul K, Tangwangvivat R, Janetanakit T, Boonyapisitsopa S, Bunpapong N, Chaiyawong S, et al. Emergence of canine parvovirus type 2c in domestic dogs and cats from Thailand. Transbound Emerg Dis. (2019) 66:1518–28. doi: 10.1111/tbed.13177
38. Nakamura K, Sakamoto M, Ikeda Y, Sato E, Kawakami K, Miyazawa T, et al. Pathogenic potential of canine parvovirus types 2a and 2c in domestic cats. Clin Diagn Lab Immunol. (2001) 8:663–8. doi: 10.1128/CDLI.8.3.663-668.2001
Keywords: canine parvovirus type 2, feline panleukopenia virus, epidemiology, genetic characterization, phylogenetic analysis
Citation: Li S, Chen X, Hao Y, Zhang G, Lyu Y, Wang J, Liu W and Qin T (2022) Characterization of the VP2 and NS1 genes from canine parvovirus type 2 (CPV-2) and feline panleukopenia virus (FPV) in Northern China. Front. Vet. Sci. 9:934849. doi: 10.3389/fvets.2022.934849
Received: 03 May 2022; Accepted: 01 November 2022;
Published: 28 November 2022.
Edited by:
Moh A. Alkhamis, Kuwait University, KuwaitReviewed by:
Soumendu Chakravarti, Indian Veterinary Research Institute (IVRI), IndiaJianye Wang, Yangzhou University, China
Francesco Mira, Experimental Zooprophylactic Institute of Sicily (IZSSi), Italy
Copyright © 2022 Li, Chen, Hao, Zhang, Lyu, Wang, Liu and Qin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tong Qin, qintong@caas.cn