Previous Article | Next Article ![]()
Journal of Virology, February 2007, p. 1961-1971, Vol. 81, No. 4
0022-538X/07/$08.00+0 doi:10.1128/JVI.01981-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Departments of Veterinary Clinical Science,1 Veterinary Pathology, University of Liverpool, Leahurst, Chester High Road, Neston, South Wirral CH64 7TE, United Kingdom2
Received 12 September 2006/ Accepted 24 November 2006
|
|
|---|
|
|
|---|
Feline calicivirus (FCV) belongs to the genus Vesivirus and is a highly infectious oral and respiratory pathogen of domestic cats (18). Following acute infection, some cats remain persistently infected with the virus (11, 40, 48, 59). Such carriers appear to be relatively common in the general population, with prevalences ranging from 15 to 91%, and they play an important role in the epidemiology of the disease (1, 3, 7, 8, 21, 22, 42). Persistent infections are also recognized in other caliciviruses, although in general these are not as well characterized as in feline calicivirus, and their role in the epidemiology of the disease is not as clear (15, 17, 35, 63).
The mechanisms of viral persistence in the Caliciviridae, both within an individual and at the population level, are not fully elucidated. Persistence within an individual is thought to be largely via immune-mediated positive selection (29, 35, 48), although other mechanisms may also play a role. Key questions that need to be resolved are the following: what contribution do such carriers make to the long-term survival of the virus at the population level, and to what extent do they drive virus evolution and diversification.
In order to address these issues, we have been exploiting endemic infection of feline calicivirus within stable household groups of cats. We have previously reported that such households tend to be infected with a small number of distinct strains of feline calicivirus (8) and provide a useful model system to monitor calicivirus evolution and persistence over long periods of time, in a defined geographic space within a natural population. In contrast, in other caliciviruses, most sampling is limited to investigation of disease outbreaks over relatively short periods, where opportunities for virus evolution are likely to be limited (15, 16, 32, 57, 64). The shape of resulting phylogenies therefore tends to be a reflection of the sampling strategy employed and, thus, intermediate evolutionary steps can only be conjectured (20). Here we use structured long-term sampling of endemically infected groups of cats to identify the relative contribution of the different evolutionary processes by which caliciviruses generate viral diversity and ensure their long-term survival in the host population.
|
|
|---|
RNA extraction, reverse transcription-PCR, and consensus sequencing. RNA was extracted from positive samples (second passage or less; QIAmp viral RNA mini kit; QIAGEN) and transcribed into cDNA (Superscript III; Invitrogen) using 500 ng of specific reverse primer (Table 1) as the first-strand primer according to the manufacturers' instructions. A 529-nucleotide region of the capsid gene, equivalent to residues 6406 to 6934 of the FCV strain F9 (5) and incorporating immunodominant region E (34, 49, 53), was amplified using Pfu DNA polymerase (Stratagene) according to the manufacturer's instructions. Each 50-µl reaction mixture contained 2 µl of cDNA and 100 ng of each PCR primer (Table 1). Thermal cycling conditions consisted of DNA denaturation at 95°C for 2 min, followed by 40 cycles of denaturation at 95°C for 30 seconds, primer annealing at 50°C for 30 seconds, and primer extension at 72°C for 90 seconds. A final extension was performed at 72°C for 5 min.
|
View this table: [in a new window] |
TABLE 1. Primers used for amplification and sequencing of variable regions C and E of the FCV capsid gene
|
Cloning of PCR products. In order to characterize the diversity of FCV within individual cats and to attempt to identify transmission events, multiple clones of 20 FCV isolates were obtained from three representative cats for which a number of serial samples were available. Cloning of PCR products was performed using the Zero Blunt TOPO PCR cloning kit (Invitrogen life technologies) essentially according to the manufacturer's instructions. Plasmids were purified using a Hamilton Microlab Star robot running the Millipore plasmid extraction kit (Montage Plasmid MiniprepB96B kit) as per the manufacturers' recommendations. Plasmids containing an insert of the appropriate size were sequenced using standard protocols.
Nucleotide sequence and phylogenetic analysis. Nucleotide distance calculations and sequence alignments were performed using the programs Distance and Pileup from version 8 of the Wisconsin package, Genetics Computer Group (13), using the Jukes-Cantor model. This model accounts for multiple substitutions and ambiguities in the sequence data, which for FCV are attributable to the quasispecies nature of the virus (48). All distance values are corrected unless otherwise specified. Where appropriate, uncorrected analyses were also carried out to allow comparison with previous epidemiological studies (44, 45). Phylogenetic analysis was performed using the program GROWTREE (Jukes-Cantor distance analysis and neighbor joining) also available in the Genetics Computer Group package and viewed using TREEVIEW (38). Support for individual nodes was sought by bootstrap analysis using 1,000 repetitions (MEGA version 3.0) (30).
Clonal sequence analysis. For those isolates for which clonal sequences were available, phylogenetic analyses and average pair-wise genetic diversity were performed using the Jukes-Cantor and neighbor-joining models (MEGA version 3.0) (30). No ambiguous sites were present in the cloned sequence data. Support for individual nodes was sought by bootstrap analysis using 1,000 repetitions (MEGA version 3.0) (30).
Simpson's equitability index.
Simpson's equitability index (E) was used in addition to the average pair-wise genetic diversity to give a further measure of intraisolate quasispecies diversity (2, 54).
![]() |
Statistical analysis.
Statistical analysis was carried out using Minitab for Windows 14.1 (Minitab Inc.). Comparisons between nucleotide distance values were analyzed using the Mann-Whitney test. The critical probability was taken as a P value of
0.05 for a two-sided alternative hypothesis.
Estimating evolution rates. Where appropriately sized subdata sets of sequences were available, evolution rates were calculated using Tipdate v. 1.01 (50). These data sets were based on consensus sequences from those cats that were persistently infected with their own strain of virus (progressive evolvers; group D, cat 4 and cat 2) (8) and from those virus strains that arrived in the group after the start of the study and for which the group ancestor sequence was available (strains D2 and E2).
Evidence for positive selection. In order to identify the mechanisms underlying the diversification of strains within individual groups, evidence of positive selection was sought using a codon-based approach as implemented in Datamonkey (41). These offer several advantages over previously described methods, including not needing to assume equal synonymous substitution rates throughout the sequence and allowing the user to choose the most appropriate model for nucleotide substitution from the original data set. Because of the different sizes of the data sets, a combination of three approaches was used. For larger data sets, a single likelihood ancestor counting (SLAC) analysis was performed suing a P value of 0.1. For smaller data sets, where type I errors may be more common, an integrative approach was taken using SLAC, fixed effects likelihood (FEL), and random effects likelihood (REL) (P = 0.25 for SLAC and FEL; Bayes Factor of 10 for REL). In all cases, ambiguities in the consensus sequence were averaged in the analysis.
Predicted structure studies. The regions of the FCV capsid predicted to be evolving under positive selection were mapped onto a predicted structure for the FCV F9 capsid protein (GenBank accession number M86379), which was modeled against the known structures for San Miguel sea lion virus (SMSV), which also belongs to the Vesivirus genus (PDB entry 2gh8.pdb) (6) and human norovirus (PDB entry1ihm.pdb) (43) using SWISS-MODEL via the ExPASy web server (52). The predicted FCV capsid was viewed and manipulated using PyMol 0.99 (12).
Nucleotide sequence accession numbers. All consensus sequences have been submitted to GenBank (accession numbers DQ397674 to DQ397814). Each sequence is identified by the following strain identification key: group and strain number, cat number, and time in months, e.g., b1c7t09 indicates group B, strain 1, from cat 7, 9 months after the start of the study.
|
|
|---|
Strain identification and variation. Phylogenetic analysis based on consensus sequences of variable immunodominant regions of the FCV capsid gene identified a total of 10 distinct viral strains among the five households. Each of these strains was separated by distance values of 26 to 51% (Fig. 1A). Their uncorrected distance values (which allow direct comparison with previous studies) ranged from 21 to 39% (data not shown), as previously reported for distinct epidemiologically unrelated strains of FCV (44). Within households, one to three distinct strains of FCV were identified. The maximum variability within each strain ranged from 5 to 24% (Fig. 1A) (with uncorrected distances of 4 to 20% [data not shown]). This is in contrast to previous studies where variability within a strain, during either disease outbreaks or persistent infection, was generally relatively low (up to 6% uncorrected distance) (44, 46). Viral strains present within households generally sustained considerable levels of diversity over time (Fig. 1B), though an occasional reduction in diversity (i.e., a putative transmission bottleneck), followed by an increase, also occurred (e.g., strain A1 at 19 months [Fig. 1B]). In each of two households (D and E), new strains were identified during the course of the study. Virus variability in these two strains increased from 0 to 5% (strain E2) (Fig. 1A) and from 0 to 14% (strain D2) (Fig. 1A and B) over 6 and 27 months, respectively, showing that in this situation viral diversity tends to increase over time.
![]() View larger version (36K): [in a new window] |
FIG. 1. (A) Unrooted neighbor-joining tree of 143 FCV consensus sequences (including the published sequence FCV F9, GenBank accession number M86379), for a 420-nucleotide region of the FCV capsid gene. One consensus sequence was excluded, as it contained two strains. Shaded areas represent the 10 distinct virus strains identified in this study. Strain names, percent diversity, and the number of variants within each strain are indicated next to the corresponding cluster. Evolutionary distances were calculated using the Jukes-Cantor model, with the scale bar indicating percent divergence. Numbers at major nodes are the bootstrap values of 80% out of 1,000 replicates (MEGA 3.0). (B) Changes in viral diversity over time for the three strains (A1, D1, and D2) for which the most data were available. Shaded areas represent the ranges of viral strain diversity (Jukes-Cantor distance analysis) present in the cat groups at the time of sampling, with putative bottlenecks at 19 months (strain A) and 6 months (strain D2).
|
![]() View larger version (38K): [in a new window] |
FIG. 2. (A) Unrooted neighbor-joining tree of 91 FCV consensus sequences from samples isolated from colony D. Phylogenetic relationships are highlighted for three representative cats, numbers 1, 4, and 6. Colored circles represent viruses isolated from each cat, with the time at which each isolate was collected shown in months. Evolutionary distances were calculated using the Jukes-Cantor model, with the scale bar indicating percent divergence. (B) All possible pair-wise comparisons of genetic distances against time for all consensus sequences for cats 1, 4, and 6. Symbols below the line represent comparisons made within a strain and were used to calculate lines of best fit. Symbols above the line for cat 1 represent comparisons made between strains. (C) Unrooted neighbor-joining trees of cloned sequences derived for viruses isolated from cats 1, 4, and 6. Colored circles represent times at which each FCV sample was obtained. The mean interisolate variability is shown in bold, and the pair-wise ranking score is in italics. Mean intraisolate variability and E values are shown in brackets. Numbers at major nodes are the bootstrap values of 80% out of 1,000 replicates (MEGA 3.0). The scale bar indicates percent divergence.
|
Interestingly, consensus and cloned sequence analyses showed that at least one cat was infected with two distinct strains of virus simultaneously (e.g., cat 1 at 6 months [Fig. 2C]).
Further evidence was sought for these two mechanisms of persistence by ranking interisolate distances between consecutive pairs of isolates from individual cats. Pairs were ranked as 1 if they were more closely related to each other than to any other viruses cocirculating in the household, suggestive of progressive evolution of a given variant within an individual. Pairs were ranked as 2 if they were more closely related to other viruses of the same strain cocirculating in the household than to each other, suggestive of transmission of variants within a strain. Pairs were ranked as 3 if they were distinct strains, suggesting between-strain transmission (Fig. 3). Significant differences were found between all three categories (P < 0.001). For the four cats that we hypothesized to be persistently infected by progressive evolution, interisolate distances were consistently ranked one (e.g., cat 4 [Fig. 2c and 3]). In contrast, for the remaining 27 cats which we hypothesized to be persistently infected by the processes of progressive evolution as well as sequential reinfection, their interisolate ranks were a mixture of 1, 2, and 3 (e.g., cats 1 and 6 [Fig. 2c and 3]).
![]() View larger version (17K): [in a new window] |
FIG. 3. (A) Bar chart showing nucleotide distance comparisons for 65 pairs of FCV isolates that were isolated on consecutive samplings from 24 cats (refer to the text for an explanation of the ranking system). (B) Box and whisker plot comparing the distributions of nucleotide distances for sequential sequences ranked 1 (progressive evolution), 2 (within-strain transmissions), and 3 (between-strain transmissions). The number of paired sequences that were placed within each ranking system is indicated. Median upper quartile (Q3, where 75% of data are less than or equal to this value) and lower quartile (Q1, where 25% of data are less than or equal to this value) distance values are shown for each rank. *, distance values outside the whiskers are classified as outliers.
|
Analysis of positive selection. Patterns of evolution were clearly associated with positive selection at known antigenic sites, suggesting an immune-mediated mechanism for viral evolution and persistence (Fig. 4). This was clear not only in the small number of individuals progressively evolving virus (e.g., cats 2 and 4), but also in those cats undergoing variant reinfection (e.g., cat 6) and for all variants of a strain at the population level (strains D2, A1, and A2) (Fig. 4). In contrast, when published sequences from epidemiologically unrelated FCV strains were compared, there was little evidence for positive selection (Fig. 4). Codons predicted to be under positive selection mapped to the 5' and 3' hypervariable regions (HVRs) of region E of the FCV capsid (Fig. 4), which have been shown to contain linear and conformational B-cell epitopes, respectively (49, 55). Positive selection was also evident in variable region C, providing for the first time evidence that this region of the capsid may also be antigenic.
![]() View larger version (43K): [in a new window] |
FIG. 4. Codon-based analysis of selection pressures on regions C and E of the FCV capsid gene (Datamonkey). The main figure shows a plot of dN-dS values for each predicted codon of the consensus sequences. The red line represents all variants of strain D1. The blue line represents previously published sequences from epidemiologically unrelated strains of FCV. For each data set, the codons predicted to be evolving under positive selection by SLAC (P = 0.1) are indicated by red and blue asterisks above the plot. Codon mutations known to disrupt neutralizing B-cell epitopes are indicated by downward-pointing arrows above the plot. The blue shaded areas correspond to FCV capsid variable regions (C, 5' and 3' HVRs). Codon 1 is equivalent to codon 383 of the FCV capsid protein (53). Bars above and below this main plot represent the results of integrated analyses (SLAC, FEL, and REL) for three other strains (D2, A1, and A2) and for three individual cats (progressive evolvers, cat 4 and cat 2 [8], and sequential reinfection with a variant of a strain [cat 6]), respectively. Gray areas correspond to those regions predicted to be evolving under positive selection in each data set (P = 0.25 for SLAC and FEL; Bayes factor of 10 for REL).
|
![]() View larger version (39K): [in a new window] |
FIG. 5. Predicted structure of the FCV F9 capsid P2 surface domain (6), highlighting regions identified in this study that are predicted to be under positive selection. The yellow region corresponds to the 5' HVR, known to contain linear epitopes. The magenta and green regions correspond to the left and right boundaries of the 3' HVR, known to contain conformational epitopes. The cyan region contains variable region C and maps close to the 3' HVR. The FCV P2 domain as modeled contains six ß-strands, labeled A, B, C, D-E, E, and F.
|
|
|
|---|
The first evolutionary mechanism identified was progressive evolution of a given variant of a strain through mutation accumulation within an individual. This was consistently present in only a small number of cats, suggesting that although a large proportion of the cat population appear to be viral carriers, only a minority of animals appear to be truly persistently infected, progressively evolving their own variant of a strain over prolonged periods of time. This helps explain the anomaly that has been observed in experimental studies, whereby in small groups of cats, persistent infections cannot always be induced (11, 28, 40, 60, 61).
The second evolutionary mechanism, sequential reinfection, which in a number of cases was interspersed with periods of progressive evolution, accounted for the majority of animals that appeared to be persistently infected. In this situation individual animals were undergoing reinfection from within the household either with a diverse variant of the same strain or with a distinct strain of virus. Reinfection with distinct strains of FCV is consistent with previous experimental studies which have shown that cats do appear to eliminate virus and that they remain susceptible to reinfection with a different strain (28, 60). However, this is the first report where sequential reinfection has been shown to occur with closely related variants of the same strain. This has important implications for both viral evolution and the epidemiology of the disease, particularly within groups of cats. Clearly, however, we cannot exclude the possibility that the original virus may still be present in such cats and is being out-competed at intervals by a new incoming virus.
Mixed infection with two distinct strains within an individual cat was also identified in one household, and although putative mixed infections have previously been reported (47), this is the first time this has been confirmed for caliciviruses by analysis of a cloned sequence. Such mixed infections are a prerequisite for the generation of recombinant viral strains, which together with the putative parent viruses have also been identified in this household and have been reported in detail elsewhere (10). Recombination in other calicivirus genera has also been described, but parental strains have not been identified (4, 25, 26, 37, 58). Recombination between variants of the same strain may also occur, but because of the close relationship between variants within a strain, such events would be more difficult to detect.
This study also demonstrated the high levels of viral diversity that can occur within a strain circulating within a population over extended periods of time. Using a corrected substitution model, levels of up to 24% (20% uncorrected) were observed within a strain within a household and up to 21% (19% uncorrected) was observed within individual animals that were persistently infected over extended periods of time. This suggests that within the overlapping population of viral variants present, new strains of FCV appear to be emerging, presumably through a combination of the evolutionary mechanisms described above. In cross-sectional studies, levels of diversity of up to 16% (uncorrected distance values) have been reported previously in some of these households (45). It is clear that such closed or semiclosed communities of endemically infected animals provide ideal conditions for the generation of viral biodiversity, and possible increased virulence, and may have parallels in other caliciviruses.
How long these strains had been circulating in each household and how this contributed to viral diversity are not known. However, in two colonies it appeared that a new strain was introduced during the course of the study, thus providing a natural opportunity to monitor viral evolution and diversification from a presumed single ancestor. In these two cases, virus variability increased from 0 to 5% and from 0 to 14% over 6 and 27 months, respectively. Thus, it appears that transmission of a new strain of virus into an individual within a household may represent a bottleneck to virus variability. Such reductions in variability may be caused either by low-dose virus transmission or by preferential selection of the fittest variant from the viral population at the time of transmission. Reductions in viral diversity by transmission bottlenecks have not previously been reported for caliciviruses but have been shown during natural transmission events for other RNA viruses (62, 65).
This study also allowed the calculation of viral evolution rates for the partial capsid region analyzed; the rate appeared to be higher (3.84 x 102 to 4.56 x 102 substitutions per nucleotide per year) for a strain circulating within a household (i.e., including transmission) compared to that for an individual cat (1.32 x 102 to 2.64 x 102 substitutions per nucleotide per year, i.e., progressive evolver, no transmission). The rate at which feline calicivirus evolved within an individual cat was lower than has been previously reported (48), although this latter study was based on the 5' hypervariable region alone. There are no other estimates available for the rate at which caliciviruses evolve within populations, although estimates have been made for other RNA viruses (24). Evolution rates depend on a number of factors, including the region of the genome selected and the various selection pressures that may be operating. In this study, the evolution rates for feline calicivirus appear to be relatively high compared to other RNA viruses, which may be in part a reflection of the impact of positive selection acting on these regions of the genome.
Positive selection in the immunodominant 5' and 3' hypervariable regions of the capsid gene was found within strains at both the individual and population level. Positive selection in the 5' HVR has been shown previously in an individual infected cat, where it was also demonstrated that such changes alter the neutralization profile (48). In other caliciviruses, such as human norovirus, immune-mediated positive selection in the P2 domain (analogous to the hypervariable regions of the FCV capsid gene) has also been suggested as a mechanism for facilitating persistent infection in an individual (35).
Positive selection within a strain at the population level has not previously been reported but can probably occur because of the closely matched immune response within the endemically infected host population. In contrast, there was little evidence of positive selection between distinct FCV strains at the population level, which is similar to the situation in other RNA viruses (31). It has been suggested that viral evolution at the population level may be a reflection of the demographic and spatial history of transmission within such populations rather than positive immune selection (23, 51). However, such conclusions may to some extent be a consequence of the sampling strategy used (20). Where in-depth sampling of host populations is carried out, as in this study, the intermediate evolutionary steps associated with positive immune selection can be clearly identified.
Interestingly, evidence of positive selection was also found in region C of the capsid, which has not previously been reported, and may suggest a possible antigenic role for this region. Previous studies have failed to identify any linear epitopes in region C (49), suggesting that if this region is antigenic it may contain either conformational B-cell or T-cell epitopes. Further support for the conformational nature of such epitopes was found in the predicted structural studies, where the codons predicted to be under positive selection within this region mapped structurally close to other known conformational epitopes in the 3' HVR (49, 55).
This study has generated insights into the evolutionary mechanisms of FCV persistence and diversification, which are likely to be generalizable to other caliciviruses. It appears that true long-term persistent infection in an individual, characterized by continuous progressive evolution, is relatively rare and that the majority of apparent carriers in the population are undergoing a combination of progressive evolution and cycles of reinfection. In human norovirus infections, persistently infected people also appear to be relatively rare, and population persistence is largely thought to be through environmental contamination and epidemic spread (14, 56). It may be, however, that in some situations the two viruses employ a similar range of evolutionary strategies and that reinfection with variants of the same strain or with a different strain of norovirus might be more common than was previously thought. Nevertheless, the minority of individuals that are true carriers may provide an important fail-safe mechanism to ensure the virus's long-term survival within a population. It will be critical, therefore, to identify these individuals to control disease and further understand the evolution of the Caliciviridae.
We thank the cat owners for allowing us to sample their cats, Chris McCracken and Ruth Ryvar for technical assistance, and Sergei Kosakovsky Pond (University of California), Dan Haydon (University of Glasgow), and Bryan Grenfell (Pennsylvania State University) for helpful discussions.
Published ahead of print on 6 December 2006. ![]()
|
|
|---|
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»