Inter- and intracontinental migrations and local differentiation have shaped the contemporary epidemiological landscape of canine parvovirus in South America

Abstract Canine parvovirus (CPV) is a fast-evolving single-stranded DNA virus that causes one of the most significant infectious diseases of dogs. Although the virus dispersed over long distances in the past, current populations are considered to be spatially confined and with only a few instances of migration between specific localities. It is unclear whether these dynamics occur in South America where global studies have not been performed. The aim of this study is to analyze the patterns of genetic variability in South American CPV populations and explore their evolutionary relationships with global strains. Genomic sequences of sixty-three strains from South America and Europe were generated and analyzed using a phylodynamic approach. All the obtained strains belong to the CPV-2a lineage and associate with global strains in four monophyletic groups or clades. European and South American strains from all the countries here analyzed are representative of a widely distributed clade (Eur-I) that emerged in Southern Europe during 1990–98 to later spread to South America in the early 2000s. The emergence and spread of the Eur-I clade were correlated with a significant rise in the CPV effective population size in Europe and South America. The Asia-I clade includes strains from Asia and Uruguay. This clade originated in Asia during the late 1980s and evolved locally before spreading to South America during 2009–10. The third clade (Eur-II) comprises strains from Italy, Brazil, and Ecuador. This clade appears in South America as a consequence of an early introduction from Italy to Ecuador in the middle 1980s and has experienced extensive local genetic differentiation. Some strains from Argentina, Uruguay, and Brazil constitute an exclusive South American clade (SA-I) that emerged in Argentina in the 1990s. These results indicate that the current epidemiological scenario is a consequence of inter- and intracontinental migrations of strains with different geographic and temporal origins that set the conditions for competition and local differentiation of CPV populations. The coexistence and interaction of highly divergent strains are the main responsible for the drastic epidemiological changes observed in South America in the last two decades. This highlights the threat of invasion from external sources and the importance of whole-genome resolution to robustly infer the origin and spread of new CPV variants. From a taxonomic standpoint, the findings herein show that the classification system that uses a single amino acid to identify variants (2a, 2b, and 2c) within the CPV-2a lineage does not reflect phylogenetic relationships and is not suitable to analyze CPV evolution. In this regard, the identification of clades or sublineages within circulating CPV strains is the first step towards a genetic and evolutionary classification of the virus.


Introduction
The parvoviruses (Parvoviridae family) are small (25 nm), nonenveloped viruses with icosahedral symmetry that infect a wide range of animal species. The subfamily Densovirinae infects invertebrates, while the subfamily Parvovirinae has diverse vertebrate hosts, including humans, pigs, rodents, cattle, and several domestic and wild carnivore species (Allison et al. 2013). Canine parvovirus (CPV) belongs to the species Carnivore protoparvovirus 1 (genus Protoparvovirus), a highly genetically related cluster of viruses that also comprises raccoon parvovirus, mink enteritis virus (MEV), and feline panleukopenia virus (FPV) (Cotmore et al. 2013).
CPV is a well-recognized cause of acute and often fatal hemorrhagic gastroenteritis in puppies and needs to be controlled by early vaccination (Decaro and Buonavoglia 2012). The virus has a linear single-stranded DNA genome (5.2 kb) with two open reading frames (ORFs) (Reed et al. 1988). The left ORF encodes nonstructural proteins 1 and 2 (NS1 and NS2), which are involved in DNA replication, capsid assembly, and viral cellular transport (Wang et al. 1998). The N-terminal region of NS1 and NS2 are identical in sequence, whereas the C-terminal region of NS2 is generated from differential splicing and is translated in a different reading frame. The right ORF encodes the viral capsid proteins 1 and 2 (VP1 and VP2), which are the main antigens that induce protective antibodies (Tsao et al. 1991;Agbandje et al.1995;Nelson et al. 2007). VP1 and VP2 are splice variants and are identical in sequence, except for a 143-amino-acid Nterminal region that is unique to VP1 and essential for viral cell infection (Tattersall and Ward 1976). In infectious virions, VP2 may be proteolytically cleaved near the amino terminus to form VP3 (Weichert et al. 1998). At both ends of the genome, there are non-translated regions with hairpin structures that prime viral DNA replication (Tattersall and Ward 1976).
CPV emerged in the late 1970s as a consequence of an interspecies jump from other carnivores to dogs, a process that involved a small number of amino acid changes in the VPs (Parrish 1999). The newly emerged CPV, referred as CPV type-2 (CPV-2), rapidly spread worldwide causing a pandemic in the new canine host (Decaro and Buonavoglia 2012). In 1979 appeared CPV-2a, a new genetic lineage that differs at five residues in VP2 and regained the ability to infect cats and other carnivores (Parrish et al. 1988;Truyen et al. 1996;Decaro et al. 2010Decaro et al. , 2011. This CPV-2a lineage replaced the original CPV-2 strains within two years and became the ancestor of all the variants that currently circulate worldwide. This early period of CPV emergence was characterized by global spreading and strain replacement (Hoelzer et al. 2008).
The CPV-2a lineage underwent further evolution, showing an intrinsically high substitution rate that resembled fastevolving RNA viruses (Shackelton et al. 2005). Changes in the residue 426 of VP2 have received considerable interest because they alter the antigenic properties of the capsid and have a geographic and temporal pattern. According to the amino acid in the 426 position, variants are denoted 2a (Asn), 2b (Asp), and 2c (Glu). The CPV-2a lineage was comprised of only 2a variants until the 2b variant was detected in the United States in 1984 (Parrish et al. 1991). The newly emerged 2b variant reached high prevalence in some countries and is currently circulating worldwide, but with different frequencies (Decaro and Buonavoglia 2012;Miranda and Thompson 2016). The 2c variant was identified around 2000 in Italy and is now detected mainly in Europe and America (Buonavoglia et al. 2001;Miranda and Thompson 2016).
Although the virus dispersed over long distances in the past, several current CPV populations have acquired particular nucleotide substitutions, which is evidence of local differentiation. This fact leads to the hypothesis that current populations are spatially confined and show only a few instances of migration between specific localities (Doki et al. 2006;Hoelzer et al. 2008;Clegg et al. 2011). However, it is unclear whether these dynamics persist in other scenarios due to the difficulty of comparing strains using only the 2a/2b/2c classification and partial VP2 sequences. To infer the emergence and spread of the genetic variants, it is critical to extend the analysis to other genomic regions and focus on virus populations from different places and periods.
In South America (Argentina, Brazil, Colombia, Ecuador, and Uruguay), CPV has shown singular characteristics that supports a significant temporal and spatial genetic variability (Pereira et al. 2007;Streck et al. 2009;Gallo Calderó n et al. 2012Pé rez et al. 2012;Aldaz et al. 2013;Duque-García et al. 2017). Events of local diversification were detected in Brazil during 1980-2000 by analyzing positively selected residues in VP1 (Pereira et al. 2007). In the early 2000s, the emergence and spread of 2c strains seemingly changed the epidemiological scenario of CPV in South America (Pé rez et al. 2007). More recently, a divergent 2a strain emerged in the Uruguayan 2c homogenous population and reached a frequency of 85% within 2 years (Maya et al. 2013).
Questions remain about the origin and spread of South American variants as most strains have been only partially sequenced; complete genomes available in South America are exclusively from Uruguay (Pé rez et al. 2014).
In this study, the complete coding genome of global CPV populations was analyzed to infer the origin and patterns of genetic variability in South America.

Genome dataset
To generate a comprehensive genome dataset, most full-length sequences available in the GenBank database were retrieved (Supplementary Table S2). The only genomes available from South America corresponded to forty-one Uruguayan strains that were previously reported by our research group (Pé rez et al. 2014).

DNA extraction, PCR, and sequencing
Viral DNA was extracted from dog feces using the QIAamp DNA Mini kit (Qiagen, Valencia, CA, USA). Full-length genome amplification was performed with a Long PCR Enzyme Mix (Thermo Scientific, Waltham, MA, USA), and PCR products were directly sequenced using internal primers according to Pé rez et al.

Sequence assembly and alignment
Sequence assembly and consensus sequences were obtained using the Lasergene Genomics Suite (DNASTAR). DNA sequences were aligned using the MAFFT algorithm (Katoh et al. 2002).

Recombination analysis
Identification of potential recombinant, parental sequences and localization of possible recombinant breakpoints were performed using the recombination detection program software, which implements seven distinct algorithms for characterization of recombinant sequences (Martin et al. 2015). Only potential recombination events detected by two or more algorithms, coupled with phylogenetic evidence of recombination, were considered significant using the highest acceptable P-value cutoff of 0.05.

Evolutionary rates and population dynamics
The evolutionary rates (measured as nucleotide substitutions per site per year, s/s/y), the ages of the most recent common ancestor (tMRCA) and internal nodes, and the geographic movement patterns of the CPV population were jointly estimated using the Bayesian Markov Chain Monte Carlo approach implemented in BEAST v1.7.5 package (Drummond and Rambaut 2007;Drummond et al. 2012). In order to minimize sampling bias, some countries were pooled together: China, India, and Japan were pooled as Asia, Albania, and Italy as Southern Europe and France and Germany as Western Europe. A single sample from Paraguay was pooled with Argentina. The Bayesian skyline plot method was used as the coalescent tree prior, under a GTRþIþG4 substitution rate model and a relaxed uncorrelated lognormal molecular clock model (Drummond et al. 2006). Migratory events were inferred from a Continuous Time Markov Chain model for the reconstruction of discrete ancestral characters. For comparison of reversible and nonreversible matrices of geographic movement, path, and stepping-stone samplings were used to determine Bayes factors (BFs). The Bayesian stochastic search variable selection approach was used to identify significant migratory events over non-zero rates. Markov Chain Monte Carlo was ran for 100 million generations and convergence of all parameters was checked by observing Effective Sample Size > 200 with TRACER v1.5 (Drummond and Rambaut 2007). The maximum clade credibility tree was obtained from the posterior distribution of trees using TreeAnnotator v1.8 (part of the BEAST package) after discarding 10% of steps as burn-in. For each node, posterior probability (PP) and posterior state probability (PSP) was collected. Phylogenetic uncertainty for each parameter was informed as 95% highest probability density (HPD) values. Maps, migration routes, trees and population effective number through time were plotted using in-house R scripts.

Informative amino acid residues and selective pressures
Identification of variable residues in the protein alignments was performed using DIVEIN web server (Deng et al. 2010). To evaluate selective pressures on individuals sites of the codon alignments, the synonymous (dS) and non-synonymous (dN) substitution rates per site were estimated using different methods available in DATAMONKEY web server: Single-Likelihood Ancestor Counting (SLAC), Fixed Effects Likelihood (FEL), Fast Unconstrained Bayesian AppRoximation (FUBAR), Internal Fixed Effects Likelihood (IFEL) and Mixed Effects Model of Evolution (MEME) (Pond and Frost 2005).

Strain sequencing
Sequences of sixty-three CPV genomes (4629 nt) were obtained from South American and European strains (Supplementary Table S1). The sixty-three newly sequenced genomes were added to the genomes retrieved from the GenBank (n ¼ 110) (Supplementary Table S2), to compile a dataset of 173 full-length sequences of Carnivore protoparvovirus 1. Sequences came from strains isolated mainly from dogs but also from raccoons (n ¼ 9), felines (n ¼ 11), and minks (n ¼ 7).

Recombination analysis
As recombination may lead to confounding the signal of vertical evolution when reconstructing phylogenetic trees, the first approach was to analyze the strains to detect recombination events. Thirteen recombinant strains were detected and characterized (Supplementary Tables S2 and S3).

Phylogenetic analysis using full-length genomes
An initial phylogenetic analysis was performed with the nonrecombinant dataset that included 160 genomes of Carnivore protoparvovirus 1. This analysis shows two distinct and highly supported genetic groups, which are referred as the FPV-and the CPV-like viruses (Hoelzer and Parrish 2010) (Fig. 1). The FPVlike group contains the viruses from felines and minks going back to the first isolates that were collected in the early 1960s. The CPV-like group splits into two branches: the original pandemic CPV-2 strains that were collected in the USA during 1978-79 and the CPV-2a lineage that includes the variants collected from 1980 onward (Fig. 1).

Phylodynamic analysis and identification of amino acid signatures
The dataset for the phylodynamic analysis (n ¼ 144) included only the viruses from the CPV-like group (CPV-2 and CPV-2a lineage). Most of the strains of the CPV-2a lineage come from domestic dogs but there are also some raccoon strains that appear grouped in the phylogeny. The tMRCA of the 2a lineage was estimated to be 36 years, calculated from the most recent analyzed strain (2014) with a 95% HPD ranging from 34 to 38 years. South American CPVs associate with global strains in four monophyletic groups or clades (Europe I, Asia I, Europe II, and South America I) with high PP ( Fig. 2A, Supplementary Figs S1 and S2).
The Europe I clade (Eur-I) includes European (Italy, France, and Albania) and South American (Argentina, Brazil, Ecuador, Paraguay, and Uruguay) strains ( Fig. 2A). This clade has Southern Europe as its most probable ancestral geographic location (PSP ¼ 0.57) ( Supplementary Fig. S3); and has strains from the United States and Western Europe (Germany) as sister branches ( Fig. 2A). The tMRCA of this clade was estimated to be 20 years with a 95% HPD interval of 16-24 years. The strains from this clade are relatively homogenous (p-distance: 2.7 Â 10 À3 and 1.6 Â 10 À3 for nucleotides and amino acids, respectively). This clade has a characteristic combination of residues that includes the 426E, which identifies strains originally denoted as 2c, but there is not a single amino acid in the NS and VP sequences that is unique (i.e. exclusively found) for the Eur-I clade ( Table 1). The strains of the Eur-I clade group with no geographic pattern, with the exception of Ecuadorian strains that cluster in a single branch and associate with Southern European strains ( Fig. 2A).
The Europe II clade (Eur-II) comprises strains from Europe (Italy) and South American strains from Ecuador and Northeast Brazil. The Eur-II clade has Southern Europe as the geographic origin (PSP ¼ 0.98) ( Supplementary Fig. S3) and a tMRCA of 30 years with a 95% HPD interval of 25-31 years ( Fig. 2A). This clade has high variability (5.6 Â 10 À3 and 4.4 Â 10 À3 for nucleotide and amino acids, respectively). The NS1 residue 544F is shared by all the Italian and South American strains (Table 1). Brazilian strains have characteristic residues in NS1 (369K, 583K, and 616N) and NS2 (92K, 125K); they also show a one nucleotide (G) deletion at position 2343 of the intronic region of VP. Ecuadorian strains have a 668E substitution in NS1, and a T440S change in VP2; some Ecuadorian strains have 139I and 426E residues in VP2 (Supplementary Table S4).

Effective population size through time
The global full-length sequences were also analyzed using the Bayesian skyline plot method to determine the population size fluctuations that occurred during the 1978-2014 period of the CPV epidemic. The Bayesian skyline plot reveals that CPV has undergone two relatively recent events of population size expansion. The first event occurred in the 1990s and the second one in the 2000s (Fig. 2B).

Significant dispersal pathways involving South American strains
The analysis herein dated the emergence of the CPV-2a lineage in the USA during 1978 (95% HPD: 1976-80). From then, the virus began expanding across continents, including Asia (BF ¼ 62) and Southern Europe (BF ¼ 56) (Fig. 3A and B). Southern Europe was a secondary hub from where the virus disseminated to Western Europe and South America (Fig. 3A). There are significant inter-and intracontinental migration events that involved South American strains (Fig. 3C). Highly significant intercontinental movements include the migration from Southern Europe to Ecuador (BF ¼ 139) and Argentina (BF ¼ 12), and the migration of the Asia-I clade from Asia to Uruguay (BF ¼ 91). The most significant regional movements (i.e. across countries borders) were from Argentina to Uruguay (BF ¼ 91) and Brazil (BF ¼ 49) (Fig. 3B and C).

Discussion
The original pandemic CPV (type 2) was described from dog samples collected simultaneously in 1978 from the USA, Europe and Australia (Appel et al. 1979). By 1980, the disease had spread worldwide because of the transportation of dogs or contaminated fomites (Hoelzer and Parrish 2010). CPV-2 was replaced 2 years later by the CPV-2a lineage (Parrish et al. 1988). The CPV-2a lineage has been incredibly successful and originated all the strains that are actually circulating in the field (Hoelzer and Parrish 2010). Data here obtained indicate that the CPV-2a lineage emerged in the USA (PSP ¼ 1) (Supplementary Fig. S3) during the late 1970 s and spread to Eurasia; this constitutes the first report of the geographic origin and the dispersal routes that this lineage underwent during its intercontinental spreading (Fig. 3A). Europe served as a second dispersion hub from where certain CPV strains arrived in South America (Figs 2A and 3C). Early records of CPV in South America dated back to 1980 in Brazil and Chile (Abalos et al. 1982;Pereira et al. 2007). Brazilian CPV strains belonged to the 2a lineage and were extremely similar to strains collected in other continents during the same period, representing a clear example of the global dispersion occurred during the early period of the CPV-2a lineage emergence (Pereira et al. 2007;Hoelzer et al. 2008). The present pattern of genetic variability in the continent is remarkably different; the heterogeneity in the viral population and the phylogenetic relationship with foreign strains suggests that long migrations and local differentiation have been driving forces for CPV evolution in South America in the recent past. Viral migrations seem to be acquiring a more predominant role also in other continents, as evidenced by the increase in the number of reports of strains that share high homology with variants from distant geographic regions (Zhao et al. 2015;Zhou et al. 2017).

Eur-I clade: origin and spread in South America
The findings here reported indicate that some European and South American strains from all the countries here analyzed are representative of a relevant clade (Eur-I) within the 2a lineage. The Eur-I clade was circulating in Europe during 1990-98 to later spread to South American territories. Southern Europe appears as the ancestral location of this clade. The first representative strains of this clade were reported in Italy in 2001 (Buonavoglia et al. 2001) and were denoted as '2c strains' based on a single amino acid position (426E). Notable, the Eur-I clade is not related with Italian strains collected during the late 1990s (1997-2000) ( Fig. 2A), suggesting that Eur-I clade would have emerged in another European country and later introduced in Italy around 2001. Strains with the 426E residue were circulating in Germany in 1996 and in Poland in 1995 but genome sequences are not available to confirm their precise origin (Decaro et al. 2007;Majer-Dziedzic et al. 2011). Soon after emerging, strains of the Eur-I clade successfully spread in several European dog populations in few years, replacing previous variants (Martella et al. 2005;Decaro et al. 2007Decaro et al. , 2009Decaro and Buonavoglia 2012). Interestingly, there was a significant rise in the CPV effective population size just before the emergence of this clade in Europe (Fig. 2B). It is possible that the ancestor of the Eur-I clade had a better fitness than previously circulating variants and infect dogs more effectively, leading to virus immune escape, vaccine failure and the subsequent population increase (Decaro et al. 2008, Decaro andBuonavoglia 2012). To our knowledge, there has been no previous phylodynamic evidence of a rapid expansion growth in the global CPV population associated with the emergence of a CPV strain with particular amino acid combinations (Fig. 2B, Table 1). Similar expansion events have been observed in other viruses with worldwide distribution, such as HIV and influenza, associated with the emergence of genetic variants (Ekong et al. 2017;Li et al. 2017).
The Eur-I strains invaded and spread in all South American countries analyzed so far (Argentina, Brazil, Ecuador, Paraguay, and Uruguay) (Pé rez et al. 2007(Pé rez et al. , 2012Streck et al. 2009;Calderó n et al. 2011;Pinto et al. 2012;Aldaz et al. 2013). Some regions of Uruguay, Argentina and Brazil have seen Eur-I strains becoming the only variant circulating in some period during the last decade. This spreading and replacement behavior was even more prominent than in Europe, where this clade has been circulating in some countries for almost 20 years but remains scarce or absent in others (Decaro and Buonavoglia 2012). The introduction in South America of the Eur-I clade was soon after its emergence in Europe and occurred independently in Argentina and Ecuador (Fig. 3C). From Argentina, the strains migrated to Uruguay and Brazil, revealing their successful spread and invasion capability. The Argentinean origin is coincident with epidemiological data, as strains from this clade were reported in Argentina during 2003(Calderó n et al. 2009), before the first description in Uruguay (2006) andBrazil (2008) (Pé rez et al. 2007;Streck et al. 2009). The second migration wave to Ecuador of the Eur-I clade is evidenced by an Ecuadorian subclade that is closely related with some strains of the ancestral South European CPV population.
The invasion and replacement dynamics of the Eur-I clade in South America are unique and had been only observed for the CPV-2a lineage that primary replaced CPV-2 and, in a lesser extent, for the 2b antigenic type that reached high frequencies in the USA during the 1980s (Parrish et al. 1991). Notable, the South American spread of the Eur-I clade correlates with a second CPV population expansion during the early 2000s (Fig. 2B). The occurrence of the population expansion inferred from genomic data is supported by an apparent increase in the number of cases in South America (Gallo Calderó n et al. 2012).
The Glu 426 residue, which was used to describe the variant 2c, is a retained ancestral character for the Eur-I clade but also occurs in unrelated strains of the Eur-II clade from Ecuador (Table 1) and therefore cannot be used as a diagnostic character for the group. The group identity is based on non-coding changes and particular residue combinations (molecular signature) in the NS and VP viral proteins (Table 1).

The introduction of the Asia-I clade in South America
Some Uruguayan strains are representative of a large Asiatic clade (Asia-I) within the CPV-2a lineage (Fig. 2A). This clade originated in Asia during the late 1980s and evolved in this continent by the acquisition of particular nucleotide and amino acid changes. The amino acid signature of the Asia-I clade is very characteristic and includes unique residues in the NS and VP (Table 1).
After emerging and evolving in Asia, the Asia-I clade was introduced to South America around 2009 (Fig. 3C). This scenario is completely concordant with previous studies; Asiatic strains from this clade were already reported in 2004 (Ying et al. 2009) but in Uruguay, they were first detected in 2010 (Pé rez et al. 2012). Interestingly, the Uruguayan strains reached a frequency of 85% in 2 years, replacing the previous variant of the Eur-I clade that was prevalent in the territory (Pé rez et al. 2012). Although the spreading was fast, this strain has not been described in the neighboring country of Argentina (Gallo Calderó n et al. 2015). This could be indicating that migration and replacement events are restricted to certain strains with particular characteristics and that not all strains have the same fitness in different regions. CPV dynamics seem to follow the changing connection between the countries as the introduction of the Asiatic clade in South America was associated with an increase in the cultural and economic connections between China and Uruguay.
The Uruguayan strains of the Asia-I clade are genetically homogenous, most likely because its presence in South America is the consequence of a single and recent introduction event. During the co-circulation of the Eur-I and Asia-I clades in South America, one recombinant strain was detected in Uruguay (Pé rez et al. 2014). This strain did not increase its frequency in the population to acquire a more significant role in CPV evolution. Similar sporadic events of recombination have been identified within CPV strains or between CPV and FPV or MEV (Mochizuki et al. 2008;Ohshima and Mochizuki 2009;Wang et al. 2012) (Supplementary Table S3). Remarkably, nine of the recombinant strains here detected are closely related and should be considered as circulating recombinant forms (Supplementary Table S3). These strains were collected in China from 2009 to 2014 and originated from a single recombinant event between an Asiatic strain of the Asia-I clade and another undetermined strain of the CPV-2a lineage (Wang et al. 2016). This circulating recombinant form is a clear evidence that recombination could strongly influence CPV evolution by creating new lineages.

Early introduction of European strains to South America
Some strains from Ecuador, Italy, and Brazil form a third clade (Eur-II) within the CPV-2a lineage. This clade appears in South America as a consequence of an early introduction from Europe to Ecuador in the middle 1980s ( Fig. 3B and C). The strains of the Eur-II clade have significant levels of variability (5.6 Â 10 À3 for nucleotide and 4.4 Â 10 À3 for amino acids) and signs of local differentiation (Table 1). The existence of relevant amino acid changes indicates that some viruses were undergoing diversification that may be related to selection pressures, a fact that is supported by the presence of positively selected residues in the molecular signature (Table 1) (Truyen et al. 1995;Decaro et al. 2009;Yoon et al. 2009). One remarkable change is the occurrence of the 426E residue in some Ecuadorian strains of the Eur-II clade. This residue emerged de novo in Ecuador and does not have the same ancestral origin that the 426E change of the Eur-I clade ( Fig. 2A, Table 1). Notably, Brazilian strains also show a conserved deletion in the VP1 intron that occurs in most isolates after 1990 but was not present in strains from other regions of the world (Pereira et al. 2007).

An exclusive South American genetic clade
Some South American strains from Argentina, Brazil and Uruguay constitute a clade (SA-I) that has evolved exclusively in South America and harbors a particular amino acid signature ( Table 1). The strains of the SA-I clade can be considered as excellent representatives of the evolution of CPV strains in the continent. The strains of the SA-I clade acquired changes by local differentiation and may have been circulating before the Eur-I and Asia-I introductions to South America. A similar situation occur in other countries that show spatial clustering and a limited international spread in different continents (Doki et al. 2006;Hoelzer et al. 2008;Clegg et al. 2011).

CPV classification and phylodynamics
The classification system that uses a single amino acid to identify variants (2a, 2b, and 2c) within the CPV-2a lineage does not reflect phylogenetic relationships of the strains and it is not suitable to analyze CPV evolution. In this regard, the identification of clades or sublineages within circulating CPV strains is the first step towards a genetic and evolutionary classification of the virus. The new classification system should take into account the phylogenetic clustering using the coding genome, the presence of unique amino acid signatures, and the ancestral geographic location of the genetic clade. The clade nomenclature can be based on the continent or territory of origin to refer to them in a consistent and standardized way. The wellcharacterized complete genome sequences of CPV here analyzed will serve as the reference or representative sequences for future studies that deal with newly identified CPV strains.

Conclusion
CPV evolution in South America is the consequence of the coexistence of different strains from multiple origins that sets the condition for competition and local genetic differentiation. Evolutionary dynamics involve locally differentiated populations that spread within or between continents and began to interact with CPV strains of the receptive population. One potential advantage of this type of evolution is the coexistence of highly divergent strains in the same territory that can compete and even recombine to allow the virus to evolve into more adapted forms. This could eventually produce impactful evolutionary events that can drastically change the epidemiological scenario in a short period. From a taxonomic standpoint, the results here obtained support that the amino acid changes that have been traditionally used for CPV typing (i.e. 426 for 2a/2b/2c strains) occurred independently in the populations and are therefore not suitable markers to study CPV evolution. Viral epidemiology requires a new classification system that takes into account phylodynamic data and the inclusion of recombinant circulating forms.

Supplementary data
Supplementary data are available at Virus Evolution online.
Conflict of interest: None declared.