Previous Article | Next Article ![]()
Journal of Virology, November 2008, p. 11096-11105, Vol. 82, No. 22
0022-538X/08/$08.00+0 doi:10.1128/JVI.01003-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Baker Institute for Animal Health, Department of Microbiology and Immunology, College of Veterinary Medicine, Cornell University, Ithaca, New York 14853,1 Center for Infectious Disease Dynamics, Department of Biology, The Pennsylvania State University, Mueller Laboratory, University Park, Pennsylvania 16802,2 Fogarty International Center, National Institutes of Health, Bethesda, Maryland 208923
Received 13 May 2008/ Accepted 2 September 2008
|
|
|---|
|
|
|---|
The extent of genetic diversity present within viral populations is determined largely by a balance between erroneous replication (measured as the mutation rate) and purifying selection (reflected in the population substitution rate, defined as the number of fixed mutations/nucleotide site/year), with particularly high error rates in RNA viruses that replicate using RNA-dependent RNA polymerases (15), or, in the case of retroviruses, reverse transcriptases. Although large double-stranded DNA viruses possess mutation rates far lower than those seen in RNA viruses, some small single-stranded DNA (ssDNA) viruses appear to both mutate (14, 39) and have substitution rates closer to those of RNA viruses than to those of double-stranded DNA viruses (reviewed in reference 16). The parvoviruses represent a particularly well-studied case, with the mean substitution rate for the canine parvovirus (CPV) capsid protein gene at
1 x 10–4 substitutions/nucleotide/year (38, 42). However, the extent and structure of within-host variability in ssDNA viruses, including coinfection as a potential source of genetic diversity, have rarely been analyzed.
Various phenomena such as recombination, nonrandom codon usage, and host effects such as hypermutation mediated by APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3B) might also affect standing genetic diversity, although their role in the evolution of ssDNA viruses is largely unclear. In addition, some rapidly evolving viruses such as human immunodeficiency virus undergo tissue-specific variation (24, 47). Tissue-specific effects may also occur during infections with parvoviruses, but this has not been analyzed previously.
CPV emerged in the 1970s as a new virus of dogs, derived from either feline panleukopenia virus (FPV) or a very closely related virus of another host (Fig. 1). After circulating undetected in dogs in Europe or Eurasia for a few years, the virus spread globally in 1978 (35). This virus is referred to as CPV type 2 (CPV-2) to distinguish it from the distantly related minute virus of canines. In 1979, a variant strain of CPV-2 (CPV-2a) emerged and replaced CPV-2 globally by the end of 1980 (37). The emerging strain differed antigenically and readily infected cats, while CPV-2 did not replicate in felines (44). In 1984, an antigenic variant of CPV-2a with apparently identical host range, which is referred to as CPV-2b, arose. CPV-2a and CPV-2b are currently cocirculating in the global dog population, but their relative frequencies appear to vary among geographic regions and are potentially also subject to temporal fluctuation (reviewed in reference 20). Natural infections of cats and wild felines with CPV have been reported (4, 18, 31), but FPV has remained the more prevalent parvovirus causing disease in cats.
![]() View larger version (35K): [in a new window] |
FIG. 1. Overview of natural virus samples analyzed. The host from which the sample was collected, the virus type, and year of collection are indicated. The likely natural host range and year of global spread are indicated. Viruses used for inoculation of kittens are marked.
|
![]() View larger version (26K): [in a new window] |
FIG. 4. Analysis of the sequence recovered from experimental CPV-13.us.81 infection of cats. (A) Summary of the expected pathogenesis of FPV in cats, including the viral location at various days postinfection (p.i.) and organs infected. (B) Characterization of viral sequences recovered from different organs, including the type of nucleotide substitution and the attained protein (for nonsynonymous changes) of each infected cats.
|
![]() View larger version (41K): [in a new window] |
FIG. 2. Mutations in the consensus sequences. The isolate consensus sequences from clinical samples of FPV (A) and CPV (B) infections were compared by using the oldest FPV and CPV samples as reference sequences in each case. The nucleotide position is indicated at the top, and the character state of individual nucleotides is indicated below each sequence. Changes in amino acid sequence are characterized above the sequence, where appropriate. Synonymous changes are indicated by black boxes, and nonsynonymous changes are indicated by red boxes, while the location in the viral genome can be inferred by referring to the genome map shown below the sequences.
|
![]() View larger version (31K): [in a new window] |
FIG. 3. Intrahost diversity in FPV and CPV samples. (A) The gene regions covered by PCR amplification and a corresponding translation map of the parvovirus genome are indicated. (B and C) Divergent viral sequences detected in animals naturally infected with FPV (B) or CPV (C) are shown. The location of mutations in the parvovirus genome and the type of nucleotide substitution are indicated for each divergent sequence, the number of individual clones analyzed for each sample and the gene product affected by nonsynonymous mutations are indicated, and the year of isolation is identified. For comparison, in the case of the CPV-infected cat described previously by Battilani et al. (4), 10 out of 14 clones analyzed for the VP2 gene harbored one or more mutations each, which distinguished the clone from the consensus sequence.
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Characterization of FPV and CPV isolates from virus samples analyzed from natural infectionsa
|
qPCR analysis. A TaqMan real-time PCR assay for the quantification of parvovirus genomes was established, targeting a conserved 60-nt region located between nt 1039 and 1114 of the CPV genome (quantitative PCR [qPCR] forward primer AAATGAAACCAGAAACCGTTGAA and qPCR reverse primer TCCCGCGCTTTGTTTCC). A TaqMan minor-groove binding probe (ABI) was designed to bind to the amplified region (TaqMan minor-groove binding probe ACAGTGACGACAGCAC). The linearized sequence of a CPV-2b genome, cloned into a plasmid, was used as an external standard for the quantification of copy numbers.
Sequence analysis. Isolate consensus sequences were published previously (19) and are available under GenBank accession numbers EU659111 to EU659121. Mutations detected in the course of this study were compared to all publicly available FPV and CPV sequences covering the respective genome regions (alignments were identical to those published previously) (19).
Some nucleotide sequences contained deletions or insertions, which were further analyzed using BLAST (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi). These sequences were aligned separately to the consensus sequence. Gene rearrangements were confirmed only if the exact boundaries of the deletion or duplication event could be inferred from the sequence information. To confirm that those rearranged sequences were present in the original samples before PCR amplification, we chose one of the gene rearrangements, detected in the FPV-kai.us.106 sample, and designed PCR primer pairs for which one primer each binds within the duplicated gene region (the sequence is available from the authors upon request).
Experimental cat infections.
All animal infections were approved by the Cornell University Animal Care and Use Committee. The original virus-containing specimens (not tissue culture passaged) were resuspended in phosphate-buffered saline (pH 7.2) and sterile filtered through a 0.22-µm filter prior to inoculation. The viral titers in CPV isolates CPV-13.us.81 and CPV-410.us.100 were determined as 50% tissue culture infective doses (TCID50) in Nordon Laboratories feline kidney (NLFK) cells as described previously (44), and the viral genome copy numbers were determined by qPCR. Groups of two
5-week-old parvovirus-seronegative kittens (Liberty Research, Waverly, NY) were inoculated through the oronasal route with 5 x 105 TCID50 of CPV-13.us.81 or with 3 x 105 TCID50 of CPV-410.us.100. Those inocula contained 5.5 x 1014 or 3.4 x 1014 viral copy numbers, respectively, as determined by qPCR. The kittens were monitored daily for clinical symptoms, and viral shedding was assessed by PCR analysis of rectal swabs. On days 0 and 6 postinoculation, serum samples were collected to confirm the absence of antiviral antibodies in hemagglutination inhibition tests as described elsewhere previously (34). On day 6 (for the CPV13 virus-challenged kittens) or 8 (for the CPV410 virus-challenged kittens) after inoculation, the kittens were euthanized, and thymus, bone marrow, and fecal samples were then collected and examined for viral DNA by PCR. PCR products (where these were generated) were cloned and sequenced as described above.
|
|
|---|
Analysis of natural FPV infections. We analyzed viral genetic diversity in clinical FPV samples collected over a 43-year period (Fig. 1), representing a variety of host tissues (Table 1). Most mutations that distinguished isolate consensus sequences were synonymous (Fig. 2A). No mutations were detected in the NS1 carboxy terminus where NS2 is encoded in a different reading frame, and only 1 nonsynonymous mutation was observed among the 12 mutations in the VP2 protein gene (8.3%).
The population structures within individual FPV samples were very homogenous and comparable across samples (Fig. 3A). Eleven of 90 sequences (12.2%) that covered the capsid protein gene harbored a mutation, with 3 of those being nonsynonymous. The VP2 Ile 101-to-Thr change was present in three of seven sequences from the same sample. Six of 93 (6.5%) NS1-covering sequences harbored one or more mutations each, 64% of which were nonsynonymous. One mutation at nt 2166 resulted in nonsynonymous mutations in both the NS1 and C-terminal NS2 reading frames (Fig. 3A).
In most FPV samples, one mutation was detected every 4 x 10–5 to 6 x 10–5 nt, and the mutations were somewhat more frequently located in the NS1 than in the VP2 coding region (see Table 3). About half of the mutations detected in viral genomes from infected feline hosts can also be found in FPV GenBank sequences, such as 10 of 20 mutations in regions with sufficient sequence data available in GenBank. Of the 12 nonsynonymous mutations, 4 (33%) were also present in GenBank FPV sequences (Table 2), indicating a likely relatively high associated fitness of these mutations. Moreover, in the cases of FPV-8.us.89 and CPV-411.us.98, the possibility of coinfections is underlined by the presence of these respective mutations in the global FPV population.
|
View this table: [in a new window] |
TABLE 3. Analysis of heterogeneity in the viral populationa
|
|
View this table: [in a new window] |
TABLE 2. Characterization of within-host mutations detected in individual FPV or CPV sequencesa
|
Analysis of natural CPV infections. Clinical CPV samples were isolated over a 22-year period, starting the year after CPV spread worldwide (Fig. 1). Most samples were in feces, except for CPV-5.us.79, which was a sample from the spleen (Table 1). The sequences were assigned to the previously defined CPV-2 (CPV-5.us.79 and CPV-6.us.80), CPV-2a (CPV-13.us.81), and CPV-2b (CPV-410.us. 00 and CPV-411.is.98) antigenic types (Fig. 1) (35, 37). Most mutations that distinguished the isolate consensus sequences clustered in the capsid protein gene (Fig. 2B). Eleven of 20 (55%) mutations in this gene were nonsynonymous, while genome-wide, 14 nonsynonymous mutations were found among the 30 mutations (46.7%). Three mutations were located in the carboxy-terminal region of NS1 where NS2 is encoded in a different reading frame, changing residue 641 of NS1 or residues 94 and 152 of NS2.
The CPV-411.us.98 sample, collected from a puppy infected during a shelter outbreak, harbored two clearly distinct viral populations (Fig. 3B). Both strains were of the CPV-2b type, and all but one of the distinguishing mutations were also detected in GenBank sequences (Table 2), indicating that this likely represented a coinfection by two different viruses. The sample also contained a single capsid protein gene clone that was of the CPV-2 type, possibly a remnant of vaccine virus. However, since the possibility of contamination in this case could not be ruled out, the sequence was excluded from further analyses.
In most CPV samples, one mutation was detected every 2 x 10–5 to 4 x 10–5 nt, and mutations again appeared to be somewhat more frequent in the NS1 than VP2 ORF (Table 3). Four of the 97 NS1 clones (4.1%) and two of the 83 VP2 clones (2.4%) harbored mutations from the consensus sequences, with individual clones harboring one or two mutations. Only one of the seven mutations detected within CPV-infected dogs (14.3%) was also detected in sequences from GenBank (Table 2). We found evidence of one gene duplication and deletion in the VP2 gene fragment of CPV-13.us.81 (Table 1), and 3 of 21 CPV-13.us.81 sequences differed in the frequencies of a repeated 62-nt sequence located in the 3' end of the viral genome (Fig. 3B), a phenomenon described previously for CPVs and FPVs (36).
Analysis of experimental cat infections. Only low levels of viral replication and little evidence of clinical disease were seen in the two kittens inoculated with CPV-13.us.81. Kitten 1 developed mild clinical signs on day 4 after inoculation and shed small amounts of parvovirus in the feces the following day, while kitten 2 showed mild clinical signs and viral shedding on day 6 postinoculation. Neither of the kittens inoculated with CPV-410.us.100 showed clinical signs, and no clear evidence of viral fecal shedding was detected by conventional PCR or qPCR (data not shown).
We analyzed a total of 29 NS1- and 22 VP2-spanning clones isolated from different tissues of the two CPV-13.us.81 inoculated animals (Fig. 4B). In general, all viral populations were highly homogeneous. Mutations were detected in the bone marrow (VP2 Asp528 to Val) and feces of kitten 1, with the later affecting the start codon of VP2. Two distinct sequences were recovered from the feces of kitten 1, with four of the eight clones being identical to the challenge virus (Fig. 3B). The other four clones harbored the same set of four mutations in the NS1 region, two synonymous changes in the NS1 carboxy-terminal sequence and nonsynonymous in the NS2 ORF (Thr94 to Ala and Met152 to Val) and two synonymous changes located in the NS1 amino terminus. All mutations except the one in the VP2 gene of kitten 1 were also present in GenBank sequences (Table 4). Gene rearrangements were detected in sequences from the thymus of kitten 1 (nt 2134 to 2875, 2474 to 3342, and 2676 to 3770), but such rearrangements were not detected in the challenge virus. Both forms of 62-nt repeat arrangements were present in the challenge virus and were also detected after experimental passage in cats.
|
View this table: [in a new window] |
TABLE 4. Mutations and rearrangements detected in individual CPV sequences after experimental passage in catsa
|
|
|
|---|
Intrahost sequence diversity. In general, we observed low levels of sequence variation during natural infections by either FPV or CPV and also during the experimental cross-species transfer of CPV to cats. Only 6.3% of viral sequences cloned from natural infections harbored any mutations. This contrasts with most previous studies of intrahost population structure during parvovirus infections but has also been suggested by one other recent study of CPV in dogs (3). Battilani et al. analyzed individual virus genomes isolated from a single CPV-infected cat and detected high levels of sequence diversity within a 1,745-nt fragment of the VP2 gene (4), with 10 distinct sequences observed among 14 analyzed viral clones (71%). Two antigenically distinct CPV variants (CPV-2a and a variant with the VP2 D426E substitution, referred to as CPV-2c) were isolated from this single animal, indicating that at least part of the observed variation was likely due to superinfection rather than newly arising mutations.
Evidence for coinfections. Notably, our study provided strong evidence for multiple infections. In particular, three of seven clones analyzed for FPV-8.us.89 contained the same change (VP2 residue Ile101 to Thr) in the capsid protein gene while being otherwise identical to the consensus sequence. This particular mutation is circulating in the FPV (Table 2) and CPV (19) populations, thereby making a coinfection event a plausible explanation. The CPV-411.us.98 isolate presented particularly strong evidence of coinfection. This sample contained two genetically distinct CPV-2b viruses (CPV-411a.us.98 and CPV-411b.us.98), which are distinguished by 6 and 12 mutations in the NS1 and VP1 genes, respectively. All but one of the mutations differentiating those two viral strains are present in GenBank sequences, again supporting the idea of coinfection since the mutations are clearly circulating in the global CPV population. If we assumed that coinfection was the source of some heterogeneity, the average number of mutations per nucleotide was in the range of 2 x 10–5 to 6 x 10–6 in all but one of the analyzed samples (the estimate for kitten 1 was excluded here, as the artificial host switching might impact the developing population structure, and the estimates might therefore not be comparable). Conversely, the estimates were considerably higher for the questionable samples if it was assumed that all variation arose by de novo mutation. These higher estimates were in the range estimated based on a previous report by Battilani et al. (3), so similar reasoning might explain the high diversity reported by those authors. Coinfections have been described for vaccine and field strains of CPV (13) and for several related parvovirus family members, including human parvovirus B19 (5). Coinfections with multiple parvovirus strains may thus occur frequently, potentially facilitating recombination (41).
Gene rearrangements. We detected several gene duplication and deletion events among the FPV and CPV genomes examined. Such rearranged sequences have long been recognized among the parvoviruses (1, 10, 21) and likely result from template switching of the polymerase during replication. This process, which is also thought to be responsible for nonhomologous and homologous recombination among RNA viruses (23, 28), may be facilitated by the complex secondary structures of the viral ssDNA genome. Recombination among CPV genomes was suggested to explain the origins of genomes containing various combinations of mutations after extended tissue culture passage (2).
Effect of host species, organ, or sampling time after CPV emergence. The viral population structures in all isolates were generally characterized by low levels of sequence diversity regardless of their host species, organ, time of isolation, or specific viral strain. On average, one mutation was detected every 104 to 105 nt, which is comparable to the annual substitution rate estimates of these viruses obtained on a population level (19, 43). Since variant CPV sequences arise very readily during tissue culture passage (2), strong purifying selection appears to quickly purge most arising mutants during natural infections, resulting in the low degree of sequence diversity within infected hosts. We detected no observable differences in intrahost population structure between CPV isolates collected during the first wave of spread (from 1979 and 1980) and those collected after the virus had been circulating in dogs for longer times, suggesting that the dynamics of intrahost mutation and selection did not change markedly during this time period.
To better determine the effects of a single cross-species transmission event on the virus, we examined experimental CPV infections of cats, which revealed only limited viral replication in the susceptible cats. The CPV sequences recovered from the cats showed low levels of heterogeneity, but parts of the viral capsid protein gene were deleted in some cases, which is indicative of defective genomes. In addition, one virus isolated from the feces of cat 1 harbored a likely lethal mutation that altered the start codon of VP2. Mutations were detected in the NS1 and NS2 regions of viruses isolated from bone marrow and feces, with approximately 50% of viruses in the feces of cat 1 harboring the same set of four mutations. Two of the changes (at residues 94 and 152 of NS2) mapped to the gene region where NS2 is derived from a different reading frame overlapping the NS1 reading frame and were identified as being potentially positively selected at a population level (19). Since these changes represent two of four apparently linked mutations observed after artificial host switching, a role of some or all of these residues in host adaptation is likely. Little is known about the function of NS2 in CPV, and NS2 knockout mutants showed no obvious differences in replication in cell culture or dogs (46). In rodent parvovirus strains MVM and LuIII, NS2 appears to be required for efficient translation, capsid assembly, and nuclear transport. In these viruses, host-specific effects modulate the functions of NS2 (9, 11, 12, 17, 25, 30). Interestingly, upon infection of severe combined immunodeficient (SCID) mice with MVM in the presence of polyclonal anticapsid antibodies, the viral population harbored nonsynonymous changes in the NS2 C terminus, which likely affected CRM1 binding (27). In this case, nonsynonymous mutations were absent from the capsid protein genes (27), although such mutations arose readily during infections in the absence of polyclonal antibodies (26). The role of NS2 in the host adaptation of CPV therefore merits further study.
Distribution of mutations in the genomes. The majority of mutations involved in the emergence of the CPV-2 cluster are in the capsid protein ORF, and this region has continued to evolve more rapidly than the nonstructural ORF, at least in CPV (19). Among the FPVs, on the contrary, marked differences in mutation accumulation levels among the genome regions have not been observed (19). Interestingly, the mutations that we detected within individual infected animals showed no clear clustering in any specific genome region. The scarcity of mutations in the capsid protein region after experimental cross-species transmission to cats is also surprising and highlights the weak effect of positive selection within individual hosts.
This work was supported by National Institutes of Health grants GM080533 to E.C.H. and AI028385 to C.R.P. K.H. is supported by a graduate assistantship from the College of Veterinary Medicine at Cornell.
Published ahead of print on 3 September 2008. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»