key: cord-0734625-32ienza1 authors: Mira, Francesco; Dowgier, Giulia; Purpari, Giuseppa; Vicari, Domenico; Di Bella, Santina; Macaluso, Giusi; Gucciardi, Francesca; Randazzo, Vincenzo; Decaro, Nicola; Guercio, Annalisa title: Molecular typing of a novel canine parvovirus type 2a mutant circulating in Italy date: 2018-07-31 journal: Infection, Genetics and Evolution DOI: 10.1016/j.meegid.2018.03.010 sha: 989c34b373a88768c74cff26bb57ca36c6ba415b doc_id: 734625 cord_uid: 32ienza1 Abstract Canine parvovirus (CPV) is the etiological agent of a severe viral disease of dogs. After its emergence in late 1970s, the CPV original type (CPV-2) was rapidly and totally replaced by three antigenic variants named CPV-2a, CPV-2b and CPV-2c. CPV has an evolutionary rate nearest to those of RNA viruses, with consequences on disease diagnosis and epidemiology. This paper reports the molecular characterization of eight CPV-2a strains collected from dogs in Italy in 2016–2017. Genetic analysis was conducted on a CPV genomic region encompassing both open reading frames (ORFs) encoding for nonstructural (NS1-NS2) and structural proteins (VP1-VP2). Sequence analysis indicates new and unreported sequence changes, mainly affecting the VP2 gene, which included the mutation Tyr324Leu. This study represents the first evidence of a new CPV-2a mutant (VP2 324Leu) and illustrates the importance of a continuous molecular survey in order to obtain more information on effective spread of new CPV mutants. Canine parvovirus (CPV) is a member of the Protoparvovirus genus (family Parvoviridae, subfamily Parvovirinae) and, together with feline panleukopenia virus (FPLV), it was recently included in the unique specie Carnivore protoparvovirus 1 (Cotmore et al., 2014; Tijssen et al., 2011) . CPV is a small, non-enveloped, linear single-stranded DNA virus. Its genome consists of an approximately 5200 nucleotide (nt) DNA molecule containing two open reading frames (ORFs), encoding respectively for two nonstructural proteins (NS1 and NS2) and for two structural proteins (VP1 and VP2) through alternative splicing of the same mRNAs (Decaro and Buonavoglia, 2012; Reed et al., 1988) . VP2 is the major capsid protein, representing the major determinant of host range and being subject to antibody mediated selection . For this reason, over the years most studies were focused on the evolution of the VP2 gene, with limited studies on the nonstructural genes (Hoelzer et al., 2008) . Also NS gene sequences available in GenBank are relatively few if compared to partial or total length VP2 sequences. CPV emerged in late 1970s as host variant of FPLV or FPLV-like parvovirus (Parrish et al., 1991) , causing acute gastroenteritis, leukopenia and myocarditis in dogs. Amino acids changes between CPV and FPLV affect the VP2 domain interacting with the host-cell transferrin receptor (TfR) and may have contributed to CPV ability of binding domestic dog TfR, which has been recognized as the determinant for the canine host shift Shackelton et al., 2005) . Soon after its emergence, the original type CPV-2 was rapidly and totally replaced by different CPV variants commonly named, on the basis of differences observed in the VP2 amino acid (aa) residue 426, as CPV-2a (Asn), CPV-2b (Asp) and CPV-2c (Glu) Decaro et al., 2005; Parrish et al., 1985 Parrish et al., , 1988 Parrish et al., , 1991 . Over the years many other changes have arisen in the VP2 gene and some of them (Ser297Ala, Thr440Ala, Tyr324Ile) were reported in several countries (Battilani et al., 2002; Decaro et al., 2009a; Geng et al., 2015; Ikeda et al., 2000; Jeoung et al., 2008; Mittal et al., 2014; Truyen, 1999; Truyen et al., 2000) . Although CPV-2a and CPV-2b variants showing Ala at aa 297 were designated respectively as "new CPV-2a" and "new CPV-2b" (Decaro et al., 2009a; Ohshima et al., 2008) , CPV strains with others specific changes (Thr440Ala, Tyr324Ile) did not receive any clear taxonomical or common designation, despite their presence in strains detected in large different geographical areas. For instance, the VP2 change Tyr324Ile may be considered a common change in CPV strains of Asian countries (Geng et al., 2015) . Therefore, current typing of CPV variants reported in the literature is based mainly on antigenic and genetic features at critical residues on the VP2 sequence (aa 297 and 426). The aim of this paper is the characterization of eight new CPV-2a strains collected from dogs in Italy between November 2016 and May 2017, which have shown common specific nucleotide and amino acid changes. The detection of CPV strains with unique molecular features can contribute to better understand the genetic diversity and molecular evolution of this canine virus. Moreover, the characterization of both viral ORFs provides useful data for further epidemiological and evolutionary analyses, improving the limited information on NS1 sequence of circulating CPV strains. Samples of eight dogs with clinical signs of severe haemorrhagic gastroenteritis and suspected of parvovirosis were analyzed. Samples were collected from November 2016 to May 2017 and analyzed at Istituto Zooprofilattico Sperimentale della Sicilia "A.Mirri" (Palermo, Italy) for diagnostic purposes. Dogs showed depression, fever, intense haemorrhagic vomiting, diarrhea, and severe leukopenia (from 2.2 to 0.55 × 10 3 cells/μl of blood). Seven dogs died within three days after the onset of clinical signs despite hospitalization and intensive cares, while the fate of the other dog was unknown. Samples collected (organs from necropsy, rectal swabs and blood sera) were submitted to virological or serological assays with suspect of parvovirosis. Details are summarized in Table 1 . Viral DNA was extracted from 200 μl of swab/organ homogenate using DNeasy Blood & Tissue Kit (Qiagen s.p.a., Milan, Italy), according to the manufacturer's instructions. All PCR assays were performed using the commercial kit GoTaq® G2 DNA Polymerase (Promega Italia s.r.l., Milan, Italy). Presence of CPV DNA was screened using a PCR protocol amplifying a 700-bp fragment of the VP2 gene (Touihri et al., 2009 ) in a 50-μl reaction mix, as previously described . PCR products were checked after electrophoresis on a 3% agarose gel supplemented with 0,005% ethidium bromide. Analyses were conducted amplifying a long genome sequence, encompassing both major ORFs, using primers pairs described by Pérez et al. (2014) . For this purpose, two separate reactions with the same reaction mix and thermal conditions were performed in 50 μl, using primer pairs NS-Fext/NS-Rext and 2161F/4835R, respectively. After electrophoresis on agarose gel supplemented with ethidium bromide, one positive amplicon of each dog from both screening and sequencing PCR analyses were purified with Illustra™ GFX™ PCR DNA and Gel Band Purification Kit (GE Healthcare Life Sciences, Amersham, Buckinghamshire, UK) and submitted to BMR Genomics srl (Padova, Italy) for direct Sanger sequencing with same primers used for amplification and also with additional internal primers, as previously described . According to an overlapping strategy, sequences were assembled using BioEdit ver 7.2.5 software (Hall, 1999) and two sequences for each sample (partial VP2 sequence and complete genome sequence) were obtained. Assembled nt sequences were submitted to nBLAST program (Zhang et al., 2000) to search related sequences into public domain databases (additional data files: Supplementary Dataset 1, 2 and 3). Phylogenetic relationships among the 324-VP2 mutants detected and the CPV reference sequences were evaluated in order to elucidate the relationships with previous available sequences. Due to the few available CPV complete genomes and single NS1 genes in comparison with the more abundant VP2 sequences and the higher variability of the VP2 gene, phylogeny was restricted to this gene. Phylogenetic analysis was performed using the best-fit model of nt substitution with MEGA7 software (Kumar et al., 2016) using ML method according to the Tamura 3-parameters model with discrete Gamma distribution (bootstrap 1000 replicates). Complete genome sequences data were submitted to the DDBJ/EMBL/GenBank databases under accession numbers reported in Table 1 . In order to have viral isolates of the CPV mutants for future studies, four samples (rectal swabs of dogs PA43847/2017 and PA5610/2017; intestines of dogs PA48686/2016 and PA3213/2017) ( Table 1) were processed by homogenization (10% w/v) in Minimum Essential Medium Eagle (EMEM; Lonza, BioWhittaker, VWR, UK) supplemented with antibiotics and antimycotic solution (1000 U/ml penicillin G sodium salt, 1 mg/ml streptomycin sulfate, 2,5 μg/ml amphotericin B; PAA Laboratories GmbH, Austria). Then, sample homogenates were centrifuged at low speed (1500 x g for 15 min at +4°C) and supernatants were collected, decontaminated at +37°C for 1 h and finally stored at −80°C until use. Supernatants were inoculated in 80% confluent cell monolayers (A-72, CRFK), cultivated in Minimum Essential Medium (MEM) added with 5% of bovine fetal serum, 1% of sodium pyruvate (A-72) and 0.1% of lactalbumin. Inoculated cells were monitored daily for maximum five days and viral growth was evaluated by detection of cytopathic effect (CPE) and PCR. A total of five passages were carried out before considering virus isolation as unsuccessful. In order to evaluate the immune status of the infected animals, the serum antibody titers for CPV were investigated in collected blood sera (Table 1) using the haemagglutination inhibition (HI) test. HI tests were carried out as previously described at 4°C using 1% pig erythrocytes and 10 haemagglutinating (HA) units of CPV-2b antigen. Serial two-fold serum dilutions were made in phosphate-buffered saline (pH 7.2), starting from a 1:10 dilution. Titers were expressed as the reciprocal of the highest serum dilution completely inhibiting viral haemagglutination. RNAs were extracted from samples using QIAamp Viral RNA Mini Kit (Qiagen S.p.A.), according to the manufacturer's instructions. Extracted DNAs/RNAs were amplified using a set of real-time (RT-)PCR assays for detection of canine distemper virus (CDV) , canine adenovirus (CAdV) type 1 and type 2 (Dowgier et al., 2016) , and canine coronavirus (CCoV) (Decaro et al., 2004) . All samples tested negative for CDV, CAdVs and CCoV by real-time (RT-)PCR assays. A single band of the expected size (700 bp) referred to a fragment of CPV VP2 gene was observed for all samples. One amplicon from each dog, obtained both from PCR and sequencing analyses, was sequenced. According to aa residues 297 (Ala) and 426 (Asn) in VP2 gene sequence, all strains were characterized as new CPV-2a. Comparing with sequences available in GenBank database, the highest identity rates were obtained with new CPV-2a strains B-2004 (99.56-99 Comparison of the obtained sequences revealed a high nt identity (99.97-99.95%) between the detected strains. In the common tract of the obtained sequences (4329 nt), differences were observed at five nucleotide residues: at nt residue 716, 1189 and 1608 of the NS1 gene (Table 2 ) and at nt residues 876 and 1731 of the VP2 gene (Table 3) . Comparison with 193 NS1 and 250 VP2 sequences (additional data file: Supplementary Dataset 1 and 2) retrieved from the GenBank database and with 156 VP2 sequences (additional data file: Supplementary Dataset 3) from CPV strains collected in Italy from 1994 to 2017, showed 17 nt changes resulting in 11 synonymous and 6 nonsynonymous substitutions (Tables 2 and 3) . NS1 sequence analysis evidenced four synonymous and five nonsynonymous substitutions (Table 2 ). Changes at nt residues 633, 716, 1189 and 1608 were observed only in our strains. Other three changes at nt residues 96, 1164 and 1790 had been previously observed in CPV strains collected in Canada from coyotes in 2014 (CPV/Coyote/C16/ NL_2014 and CPV/Coyote/C55/NL_2014), from dogs in Uruguay in 2010 (CPV-2c-UY235) and in USA in 1998 (strains CPV-412 and CPV-402) and in 2003 (strains CPV-435 and CPV-436) . Change at nt residue 1048 was common to 9 CPV and 2 MEV strains, while change at residue 1631 was present in 38 CPV strains collected from dogs and 2 CPV strains collected from coyotes (CPV/Coyote/C16/NL_2014 and CPV/ Coyote/C55/NL_2014). None of these changes were observed in FPLV NS1 sequences (Supplementary Dataset 1) . Four synonymous and one non-synonymous substitutions (T970C, A971T → Tyr324Leu or A970C → Ile324Leu) were observed in the VP2 sequences (Table 3) . Nucleotide changes at residues 876 and 1731 were shared only by our strains. Change T1482C was observed only in 22 new CPV-2b strains collected in Italy, France, Turkey, South Korea and mainly in USA. Change A1710G was observed only in 30 new CPV-2a strains collected mainly in Europe (Italy, France, Germany, Portugal, Hungary) as well as in Turkey, Vietnam and China. The main difference was observed at VP2 aa residue 324 where our sequences showed aa Leu in place of Tyr, as for the original type CPV-2 and for all Italian CPV strains, or of Ile, as reported for the most recent Asian CPV strains. The non-synonimous substitution Tyr324Leu resulted from two nt changes at the first and second codon positions (TAT → CTT), while change Ile324Leu was due to a first codon nt mutation (ATT → CTT) (Supplementary Dataset 2,3). The aa change at VP2 residue 324 was recently reported in CPV-2b strains from Brazil (Silva et al., 2017) , but the codon analysis showed a difference at the nt level (TTA and CTT in Brazilian and Italian sequences, respectively). Moreover, by sequence analysis of both ORFs, additional nt differences were encountered in the NS1 (twenty-four synonymous and ten non-synonymous substitutions) and in the VP2 sequence (sixteen synonymous and one non-synonymous substitutions). At key positions of the VP2 protein all our sequences showed aa Ala at residue 297 and Thr at residue 440. VP2 sequence analysis of previous collected Italian CPV type 2a/2b/ 2c strains (Supplementary Dataset 3) showed mainly five synonymous changes at nt residues 36, 303, 537, 1482 and 1710 (Table 3) . Comparison with 52 Italian CPV-2a sequences evidenced that change T1482C was unique to our sequences while other changes were already present in Italian CPV-2a strains from 1999 (accession number KF373570), 1996 (AF306445) and 1994 (AF306446) with low frequency rate (15-31%), except for residue 1710 (68%). While changes at aa positions 537 and 1710 were observed also in CPV-2b and CPV-2c Italian sequences, changes at positions 36, 303 and 1482 were not observed in these variants with exception of two CPV-2b sequences (accession numbers KF373599 and FJ005265) (Supplementary Dataset 3). The phylogenetic tree inferred from VP2 sequences (Fig. 1) shows that the eight novel sequences detected in this study cluster in a distinct clade located between two main clusters including almost all European and Asian strains, respectively. The only exceptions in these clusters are VP2 sequences of strains Ku5_08 (Thailand, 2008; VP2 324Tyr), located within the European cluster, and of strains H-Erd, H-Illatos, H-Halmai, H-41 (Hungary, 2012; VP2 324Ile), located within the Asian cluster. Sequences cluster also according to variations at aa residue 324: in the phylogenetic tree, sequences displaying Tyr (Y) or Ile (I) at residue 324 cluster separately from our sequences presenting Leu (L). However, due to the few differences among available sequences and the low bootstrap values, the tree topology could not reflect the real clustering of the two 324-VP2 mutants. The results of HI test carried out on the sera of dogs PA43847/2016 F. Mira et al. Infection, Genetics and Evolution 61 (2018) [67] [68] [69] [70] [71] [72] [73] and PA5610/2017 were respectively 1:1280 and 1:640, accounting either for active CPV infection or recent vaccination. Canine parvovirus still remains one of the most important enteric pathogens of dogs in all continents, despite the widespread use of vaccines (Parrish et al., 1991; Miranda and Thompson, 2016) . Therefore, a continuous epidemiological surveillance is needed to detect new CPV variants (Decaro and Buonavoglia, 2012) . The aim of this study was to analyze and describe unreported molecular features of CPV strains recovered from dogs in southern Italy. The genetic analysis was carried out on both CPV open reading frames, thus representing the first study on the complete genome of Italian strains. The present study led to the identification of a novel CPV-2a mutant, displaying previously unreported molecular features. Considering the established role F. Mira et al. Infection, Genetics and Evolution 61 (2018) 67-73 of key residues of the VP2 protein, the main finding was observed at the critical amino acid residue 324 (Tyr324Leu). This finding provides new insight into the evolution of the VP2 protein gene, although further studies are needed to fully understand the biological and immunological role of this amino acid change. The full-length genomes provided in the present study will help further characterization of circulating strains, with particular regards to the analysis of the NS1 gene for which only few sequences are currently available. Previous field studies indicate CPV-2a and CPV-2c as the most prevalent types in Italy, with CPV-2c having very rapidly replaced CPV-2b and displaying a higher prevalence in the last years (Decaro and Buonavoglia, 2012; Decaro et al., 2006) . Epidemiological investigations, carried out in the same geographic area of the samples of the present study, detected CPV-2c as the most prevalent type with a minimal circulation of other CPV-2 types (including CPV-2a) from 2009 to 2015 . However, those CPV-2a strains did not show a relevant sequence variability. Analogously, sequence analysis of the VP2 gene from 52 recent CPV-2a strains collected in Italy (Supplementary Dataset 3) showed a relatively low level of variability during the last years (20 aa changes) and among the detected changes only the aa residue 440 showed multiple (Thr/Ala/Ser) and frequent (84.6%/13.5%/1.9%) substitutions. The strains detected in the present study were typed as new CPV-2a according to 297Ala and 426Asn aa residues of the VP2 gene. While aa Asn at VP2 residue 426 is critical for CPV-2/CPV-2a typing, aa Ala at residue 297 is commonly used as marker of the new CPV-2a/2b variants (Geng et al., 2015) , in order to differ them from old strains. This mutation is predominant among CPV-2a/2b strains since 1990 . The detected CPV-2a strains displayed a non-synonymous change at aa residue 324 of the VP2 protein. This critical residue is located into the VP2 region that forms the GH loop located between the βG and βH strands and, along with aa residues 297 and 440, it is comprised into the VP2 region 2 (gene positions from nt 883 to 1332) that forms the 22 Å long spike of the capsid surface on the threefold axes (Horiuchi et al., 1998) . Hoelzer et al. (2008) predicted that residue 324 is under positive selection, so that mutation in this site arise independently in different geographic locations. This residue is adjacent to aa residue 323, which affects canine transferrin receptor (TfR) binding Hueffer et al., 2003; Mittal et al., 2014) , contributing to viral replication into canine cells and determining the pH dependence of haemagglutination (Chang et al., 1992) . This residue, together with aa residue 93 of VP2, has been observed to play a crucial role in controlling CPV host range and in determining a specific antigenic site on the viral capsid (Govindasamy et al., 2003) . Residue 324 is also likely to have the same effect as residue 323 on CPV host range (Mittal et al., 2014) . Although the introduction of CPV strains from other countries/ continents was previously observed , the nt differences at the codon level with CPV sequences with the same VP2 324Leu mutation described by Silva et al. (2017) seem to rule out a southern America origin of the Italian mutant strains. Instead, a convergent evolution could be speculated that resulted in the same aa substitution even with different nt mutations. Moreover, the lack of related CPV sequences from the same area prevented any further speculation about the origin of the mutant. Therefore, only extensive and retrospective sequence analyses could clarify whether these mutant emerged from Fig. 1 . Maximum Likelihood tree based on the full-length VP2 gene sequences (1755 nucleotides) of canine parvoviruses displaying the genetic relationships between analyzed strains and 90 CPV-2a strains obtained from GenBank (bootstrap 1000 replicates; bootstrap values shown greater than 70%). Black dots markings (•) indicate CPV strains analyzed in this study, white dot markings (o) indicate original CPV-2a strains (297Ser), white triangle markings (Δ) indicate CPV-2a strains with change Thr440Ala. All sequences were indicated with: country and year of collection/submission, type, isolate/ strain and accession number. F. Mira et al. Infection, Genetics and Evolution 61 (2018) 67-73 CPV strains previously circulating in the same geographic area or was imported from different countries. Another change affecting the same residue (Tyr324Ile) was observed since 2006 (Jeoung et al., 2008) and it became predominant in Asia in 2014 . This substitution was observed in 100% of recently analyzed Asian strains (Geng et al., 2015; Guo et al., 2013; Xu et al., 2015) and can result in a stronger receptor binding (Cságola et al., 2014) . Further studies are required to better address any potential biological effect of this change and to evaluate the spread of this CPV mutant. A high mutation substitution rate, similar to those of RNA viruses (about 10 −4 substitution/site/year), was evidenced for CPV (Shackelton et al., 2005) , thus accounting for its rapid evolution and occurrence of new variants. The region comprised between aa residues 267 and 498 of the VP2 protein, due to its exposure on the capsid surface, is affected by the greatest variability in parvoviruses (Decaro et al., 2009a) and this antigenic region represents the target for neutralizing antibodies (Agbandje et al., 1993; Tsao et al., 1991) . Moreover, analysis of the CPV codon usage and evolutionary process suggested that VP2 aa residue 324 is one of the most affected substitution site, potentially able to alter relevant biological characteristics of the virus and to contribute to virus immune escape via antigenic drift and consequent vaccine failure . According to this evolutionary scenario, after substitution of original Tyr by Ile at VP2 amino acid residue 324, our findings account for the occurrence of a second aa change at the same residue. Evidence for a new CPV-2a mutant increases the need of a revision of current CPV typing. Due to the presence of a single aa mutation (position 426 of the VP2 protein) among types CPV-2a, 2b and 2c types, Organtini et al. (2015) suggested to consider these types as mutants of the CPV-2a clade rather than distinct subtypes. Indeed, the presence of additional aa changes in the VP2 protein (residues 297, 324, 440) requires a clearer taxonomical designation or a common naming, even if they can be considered as markers of the worldwide distribution of circulating CPV strains. To date, a little information is available for changes in the NS1 protein and for their potential role. NS1 is a multifunctional nuclear phosphoprotein responsible for viral replication, cellular transcription and cell death inducing apoptosis (Gupta et al., 2016; Saxena et al., 2013) . Only recently some of CPV NS1 amino acids were associated to specific functions (Niskanen et al., 2010 (Niskanen et al., , 2013 . In order to contribute to future studies on this gene, we provided a complete genome NS1 sequencing for a number of CPV strains. Two dogs tested in this study (PA43847/2016 and PA5610/2017) had been vaccinated with at least four doses of CPV-2 or CPV-2b based vaccines according to the guideline of the Vaccination Guidelines Group (VGG) of the World Small Animal Veterinary Association (WSAVA) (Day et al., 2016) before displaying CPV clinical signs and shedding. Previous reports accounted for the occurrence of CPV infection in vaccinated dogs Decaro et al., 2007 Decaro et al., , 2008 Decaro et al., , 2009b Miranda and Thompson, 2016; Mittal et al., 2014; Parker et al., 2017) , whereas other studies demonstrated a high prevalence of nonresponders especially in some breeds (Houston et al., 1996; Kennedy et al., 2007) . The detection of intermediate HI antibody titers in the two CPVaffected cannot be useful to assess whether the multiple vaccine administrations resulted in the development of non-protective immunity or they were completely unsuccessful because dogs were not responders or had high maternally-derived antibody titers. Serological tests should be routinely carried out after CPV vaccination to assess the seroconversion of vaccinated dogs, especially in putative "high-risk" dog breeds that include most low-responders and non-responders. At the same time, further studies are required to evaluate any additional potential mechanisms of antigenic escape of the VP2 324Leu mutants. In conclusion, this study is the first evidence of the active circulation of a novel Carnivore protoparvovirus 1 type 2a mutant with specific molecular signatures in a canine population with high prevalence of CPV-2c. Further studies are required to better elucidate the prevalence and the potential spread of this mutant in other geographic areas. Moreover, this field study could contribute to evaluate new vaccinations strategies based on the molecular epidemiological evidences. Supplementary data to this article can be found online at https:// doi.org/10.1016/j.meegid.2018.03.010. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Structure determination of feline panleukopenia virus empty particles Genetic analysis of canine parvovirus isolates (CPV-2) from dogs in Italy Evidence for evolution of canine parvovirus type 2 in Italy Characterization of a canine parvovirus strain isolated from an adult dog Multiple amino acids in the capsid structure of canine parvovirus coordinately determine the canine host range and specific antigenic and hemagglutination properties The family Parvoviridae Analysis of the full-length VP2 protein of canine parvoviruses circulating in Hungary Vaccination Guidelines Group (VGG) of the World Small Animal Veterinary Association (WSAVA), 2016. WSAVA guidelines for the vaccination of dogs and cats Canine parvovirus-a review of epidemiological and diagnostic aspects, with emphasis on type 2c Quantitation of canine coronavirus RNA in the faeces of dogs by TaqMan RT-PCR New approaches for the molecular characterization of canine parvovirus type 2 strains Characterisation of the canine parvovirus type 2 variants using minor groove binder probe technology Occurrence of severe gastroenteritis in pups after canine parvovirus vaccine administration: a clinical and laboratory diagnostic dilemma Evidence for immunisation failure in vaccinated adult dogs infected with canine parvovirus type 2c Genetic analysis of canine parvovirus type 2c Severe parvovirus in a 12-year-old dog that had been repeatedly vaccinated A duplex real-time PCR assay based F. Mira et al. Infection on TaqMan technology for simultaneous detection and differentiation of canine adenovirus types 1 and 2 Antibody levels and protection to canine parvovirus type 2 Detection of canine distemper virus in dogs by real-time RT-PCR 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 Structures of host range-controlling regions of the capsids of canine and feline parvoviruses and mutants Identification of canine parvovirus with the Q370R point mutation in the VP2 gene from a giant panda (Ailuropoda melanoleuca) Canine parvovirus NS1 induced apoptosis involves mitochondria, accumulation of reactive oxygen species and activation of caspases Asymmetric binding of transferrin receptor to parvovirus capsids BioEdit: a user-friendly biological sequence alignment editor and analysis program fo windows 95/98/NT Phylogenetic analysis reveals the emergence, evolution and dispersal of carnivore parvoviruses Differences in the evolutionary pattern of feline panleukopenia virus and canine parvovirus Risk factors associated with parvovirus enteritis in dogs: 283 cases (1982-1991) The natural host range shift and subsequent evolution of canine parvovirus resulted from virus-specific binding to the canine transferrin receptor Predominance of canine parvovirus (CPV) in unvaccinated cat populations and emergence of new antigenic types of CPVs in cats Genetic analysis of VP2 gene of canine parvovirus isolates in Korea Factors influencing the antibody response of dogs vaccinated against rabies MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets Evolutionary and genetic analysis of the VP2 gene of canine parvovirus Introduction of Asian canine parvovirus in Europe through dog importation Canine parvovirus in vaccinated dogs: a field study Molecular typing of canine parvovirus strains circulating from 2008 to 2012 in an organized kennel in India reveals the possibility of vaccination failure Different mechanisms of antibody-mediated neutralization of parvoviruses revealed using the fab fragments of monoclonal antibodies Effect of ATP binding and hydrolysis on dynamics of canine parvovirus NS1 Mutations in DNA binding and transactivation domains affect the dynamics of parvovirus NS1 protein Chronological analysis of canine parvovirus type 2 isolates in Japan Global displacement of canine parvovirus by a host-adapted variant: structural comparison between pandemic viruses with distinct host ranges Investigation of a canine parvovirus outbreak using next generation sequencing Natural variation of canine parvovirus The global spread and replacement of canine parvovirus strains Rapid antigenic-type replacement and DNA sequence evolution of canine parvovirus Phylogenetic and genome-wide deep-sequencing analyses of canine parvovirus reveal co-infection with field variants and emergence of a recent recombinant strain Investigation on canine parvovirus circulation in dogs from Sicily (Italy) by biomolecular assay Nucleotide sequence and genome organization of canine parvovirus Apoptosis induced by NS1 gene of canine Parvovirus-2 is caspase dependent and p53 independent High rate of viral evolution associated with the emergence of carnivore parvovirus Full-length genomic and molecular characterization of Canine parvovirus in dogs from North of Brazil The family Parvoviridae Molecular characterization of canine parvovirus-2 variants circulating in Tunisia Emergence and recent evolution of canine parvovirus Distribution of antigen types of canine parvovirus in Switzerland, Austria and Germany The three-dimensional structure of canine parvovirus and its functional implications Phylogenetic analysis of canine parvovirus isolates from Sichuan and Gansu provinces of China in 2011 A greedy algorithm for aligning DNA sequences The genetic evolution of canine parvovirus -a new perspective The authors have declared that no competing interests exist. None