Skip to main page content

• HOME
• Current Issue
• Archive
• Alerts
• About ASM
• Contact us
• Tech Support
• Journals.ASM.org

This Article Abstract Full Text (PDF) Supplemental material Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Fargette, D. Articles by Konaté, G. Search for Related Content PubMed PubMed Citation Articles by Fargette, D. Articles by Konaté, G.

 Previous Article  |  Next Article 

Journal of Virology, April 2008, p. 3584-3589, Vol. 82, No. 7
0022-538X/08/$08.00+0     doi:10.1128/JVI.02506-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Rice Yellow Mottle Virus, an RNA Plant Virus, Evolves as Rapidly as Most RNA Animal Viruses ^,

D. Fargette,^1* A. Pinel,¹ M. Rakotomalala,² E. Sangu,³ O. Traoré,⁴ D. Sérémé,⁴ F. Sorho,¹ S. Issaka,⁵ E. Hébrard,¹ Y. Séré,⁶ Z. Kanyeka,³ and G. Konaté⁴

Institut de Recherche pour le Développement (IRD), UMR RPB, BP 64501, 34394 Montpellier cedex 5, France,¹ Centre National de la Recherche Appliquée au Développement Rural (FOFIFA), BP 289, Mahajanga 401, Madagascar,² Botany Department, Dar es Salaam University, P.O. Box 35091, Dar es Salaam, Tanzania,³ Institut de l'Environnement et de Recherches Agricoles (INERA), Laboratoire de Biotechnologie et de Virologie Végétale, 01 BP 476, Kamboinsé, Ouagadougou, Burkina Faso,⁴ Institut de Recherche Agronomique du Niger, BP 60, Kollo, Niger,⁵ Africa Rice Center (WARDA), 01 BP 2031, Cotonou, Bénin⁶

Received 22 November 2007/ Accepted 4 January 2008

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 
The rate of evolution of an RNA plant virus has never been estimated using temporally spaced sequence data, by contrast to the information available on an increasing range of animal viruses. Accordingly, the evolution rate of Rice yellow mottle virus (RYMV) was calculated from sequences of the coat protein gene of isolates collected from rice over a 40-year period in different parts of Africa. The evolution rate of RYMV was estimated by pairwise distance linear regression on five phylogeographically defined groups comprising a total of 135 isolates. It was further assessed from 253 isolates collected all over Africa by Bayesian coalescent methods under strict and relaxed molecular clock models and under constant size and skyline population genetic models. Consistent estimates of the evolution rate between 4 x 10^–4 and 8 x 10^–4 nucleotides (nt)/site/year were obtained whatever method and model were applied. The synonymous evolution rate was between 8 x 10^–4 and 11 x 10^–4 nt/site/year. The overall and synonymous evolution rates of RYMV were within the range of the rates of 50 RNA animal viruses, below the average but above the distribution median. Experimentally, in host change studies, substitutions accumulated at an even higher rate. The results show that an RNA plant virus such as RYMV evolves as rapidly as most RNA animal viruses. Knowledge of the molecular clock of plant viruses provides methods for testing a wide range of biological hypotheses.

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 
The mutation rates of RNA viruses (i.e., the number of nucleotide misincorporations per site and per round of replication) are 10⁴ to 10⁵ times higher than those of their DNA hosts (8, 9). Such a high mutation rate is attributed to the lack of repair function of the RNA polymerase of these viruses, the short replication times, and the large populations in infected hosts (7). A high mutation rate often results in rapid evolution of RNA animal viruses. This allowed the measure of the evolution rates of an extensive range of animal viruses through the analysis of heterochronous sequences, i.e., sequences of viral genes isolated at different times (11, 36). A large variation in the evolution rates of RNA animal viruses was subsequently found and was attributed mostly to differences in replication rates (24, 26). Estimates of evolution rates are used increasingly to date the emergence and to reconstruct the population dynamics of major viral epidemics (6).

Interestingly, some RNA viruses change little or not at all over time. The best-documented example is an RNA plant virus, Tobacco mild green mosaic virus, which showed no increase in genetic diversity over the 90 years considered, in the longest series of isolates with known isolation times for any virus (20). Indeed, many studies have shown the remarkable genetic stability of RNA plant virus populations from different geographical regions, hosts, and collection times (21). It was claimed that most tobamovirus populations are very stable and do not evolve at a measurable rate (22). It was even observed that populations of Turnip yellow mosaic virus from Europe and Australia that probably separated more than 12,000 years ago differed by less than 1% (4).

Actually, the lack of estimates of evolution rates of RNA plant viruses over time may merely reflect the absence of a large enough number of isolates collected over a sufficiently long period. Heterochronous data for plant virus isolates are particularly difficult to gather compared to data for animal viruses, for which isolates are readily recovered from blood samples that have been collected over many decades and stored for medical purposes. Even then, the temporal component of variation can be blurred or biased by other sources of diversity, such as long range dispersal, recombinant events, and subpopulation division, which are common features of plant viruses.

Since the 1920s, experimental evidence has established that RNA plant viruses can evolve rapidly, especially under selection pressures such as a change of host (21, 22). The evolution rate of Wheat streak mosaic virus was extrapolated from serial passage experiments (38). This method assumed that the mean rate of change measured in the laboratory reflected that of the natural viral populations, although constraints on evolution in nature and in experiments may differ. Recently, the evolution rate of Barley yellow dwarf virus was calculated by comparing an isolate preserved in old herbarium specimens to present-day specimens (29). This method postulates that the genetic diversity of the population at the time of sampling is negligible, so that sequences isolated at different times differ only by substitutions accumulated during the time interval. If this condition is not met, this method has an upward bias and overestimates the evolution rate (11). In addition, both attempts assumed that the molecular clock remained constant during the evolution of the viruses. The estimates of the evolution rates of these two RNA plant viruses therefore rely on critical but nontestable assumptions. Estimates of the evolution rates of plant viruses based on historical evidence such as outbreak records can be tentative only. Altogether, the evolution rate of an RNA plant virus has never been estimated by applying statistical methods developed to analyze temporally spaced sequences. This contrasts with the recent advances made for RNA animal viruses by using this statistical approach (11).

Rice yellow mottle virus (RYMV), of the Sobemovirus genus, was used to estimate the evolution rate of an RNA plant virus. RYMV has a high natural molecular diversity (18) and reaches a high concentration in rice (19). Experimentally, RYMV adapts rapidly to alternative hosts through accumulation of point mutations (25, 34). These features make RYMV an appropriate model with which to estimate the evolution rate of an RNA plant virus species. RYMV is an emergent virus, indigenous to Africa. It was first noticed in Kenya in East Africa in 1966 (3) and since then in almost all African countries where rice is grown, including Madagascar (39). RYMV is a major threat to rice cultivation (28). It has a narrow host range, restricted to wild and cultivated rice species and a few related grasses (2). RYMV is transmitted primarily by coleopterous beetles of the family Chrysomelidae and is disseminated by contact during cultural practices (28). Its genome contains four open reading frames (ORFs) (18). ORF1, located at the 5' end of the genome, encodes a protein involved in virus movement and gene silencing suppression. ORF2, which encodes the central polyprotein, comprises two overlapping ORFs. ORF2a encodes a serine protease and a viral-genome-linked protein (VPg). ORF2b, which is translated through a –1 ribosomal frameshift mechanism as a fusion protein, encodes the RNA-dependent RNA polymerase. The coat protein (CP) gene (ORF4) is expressed by a subgenomic RNA at the 3' end of the genome.

RYMV isolates were collected in 16 African countries between 1966 and 2006. The CP genes of 253 isolates were sequenced in this or in earlier studies (1, 18, 33, 39). Such a large collection of sequences over a 40-year period is unique for a plant virus and is used here to estimate the rate of change of RYMV. RYMV diversity is geographically structured, with different strains in East, West, and Central Africa (1, 33, 39). Substitution rates were first estimated by pairwise distance linear regressions from five phylogeographically based groups of isolates. Each group comprised isolates collected over the longest possible period of epidemiological record while other factors influencing diversity, which might adversely affect the analysis of temporally spaced viral sequences, were minimized (11). The five groups comprised 135 isolates in total. Rates were further assessed by Bayesian coalescent methods using 253 isolates originating from all regions of Africa. The rates were calculated under strict and relaxed molecular clock hypotheses (10) and under constant-size and skyline population models (14). The synonymous evolution rate was estimated from the evolution rate of the third codon position. The overall and synonymous evolution rates of the CP gene of RYMV were compared to the evolution rates of 50 RNA animal viruses (24, 26). Experimentally, the number of changes was calculated by comparing the full sequence of each isolate at inoculation with that 1 to 6 months later. Then the number of changes was compared to that estimated from the evolution rate. Altogether, we found that the overall and synonymous evolution rates of RYMV were within the ranges of those of RNA animal viruses. This shows that an RNA plant virus such as RYMV evolves as rapidly as most RNA animal viruses.

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 
Data sets. The CP genes (720 nucleotides [nt]) of 253 isolates from 16 countries in Africa collected over a 40-year period (70 isolates collected in this study and 183 isolates in earlier studies) were sequenced. The isolates were inoculated once into the susceptible cultivar IR64 to increase the virus content. Extraction, transcription, and amplification sequencing were carried out as described previously (33). For convenience, the strain nomenclature was that used previously based on serotyping and partial sequencing (32, 33). The designation, sampling year, and sequence accession number of each of the 253 isolates are given in Table S1 in the supplemental material.

Five phylogeographically based groups of isolates were defined for pairwise distance linear regression estimates (see Table S1 in the supplemental material). Within each group, isolates belonged to the same phylogenetic cluster, originated from the same region, had the least population subdivision, and showed no evidence of long-distance movement or of recent recombination events. Group I, also referred to as S2/S3-CI, comprised 40 isolates from Côte d'Ivoire collected over 27 years. Strains S2 and S3, considered earlier to be distinct, were combined as a more intensive survey revealed a continuum between the two strains and are now referred to as strain S2/S3 (34). Group I had the densest sampling survey, the least residual variability, and minimal population subdivision. Group II (S2/S3-WA), with 68 isolates collected over 31 years, was an extension of group I that included 28 additional S2/S3 isolates from Guinea, Mali, and Sierra Leone, three countries that neighbor Côte d'Ivoire in West Africa. Group III (S1-WA) comprised 23 S1 isolates from Burkina Faso, Côte d'Ivoire, and Mali in West Africa, collected over 11 years, a shorter period, and displaying a higher residual diversity. Group IV (S1'-CA) included 31 isolates from Benin, Cameroon, Chad, Niger, Nigeria, and Togo collected over 23 years. These isolates formed a monophyletic cluster referred to earlier as the S1'-Central African strain to designate isolates from Cameroon and Chad in Central Africa; then, by extension, this designation was applied to isolates of the same lineage from countries to the west (Benin, Niger, Nigeria, Togo), although geographically these countries belong to West Africa. The isolates originated from a large geographic area and had a relatively high residual diversity. Group V (S4-LV) comprised 13 isolates of strain S4 collected over 34 years that originated from western Kenya, northern Tanzania, and southern Uganda around Lake Victoria. The sampling time was the longest but the number of isolates the smallest. The five groups comprised a total of 135 isolates. An additional set of 118 isolates from all parts of Africa was included to make a group of 253 isolates for the Bayesian coalescent methods.

Pairwise distance linear regressions. Pairwise distance linear regressions are fast and useful methods to estimate the substitution rate from large data sets of heterochronous sequences. These methods are also used to determine if the sequence diversity exhibits an adequate temporal structure and to verify that no outlier substantially alters the substitution rate (26). With RYMV, the genetic distance (percentage of nucleotide differences calculated over the 720-nt CP gene) and the separation time (sampling interval in number of years in absolute value) of each pair of isolates were determined. The genetic distance was plotted against the separation time. The significance of the correlation coefficient between genetic distance and separation time (30) was tested after 100,000 permutations by using GENETIX (version 4.05; Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université Montpellier II [http://www.genetix.univ-montp2.fr/genetix/intro.htm]). The slope of the linear regression of the pairwise genetic distance against the separation time is an estimate of the evolution rate (in nucleotide substitutions per site per year) (11, 12). A parametric bootstrap method for obtaining the variance of the rate was implemented in STATISTICA (StatSoft). The confidence interval was calculated with the percentile method: the 26th and 975th estimates after 1,000 replicates (when ranked) were, respectively, the upper and lower 95% confidence limits of the original estimate. However, this bootstrap procedure underestimates the true confidence interval, since it does not consider the correlation of bootstrap replicates due to shared ancestry (11). Thus, pairwise distance linear regression is used here as an exploratory analysis of the evolution rate of the CP gene of RYMV.

Bayesian estimates. The best-fitting nucleotide substitution model was evaluated by hierarchical likelihood ratio tests, using ModelTest on the basis of hierarchical likelihood ratio testing (35), as implemented in HyPhy (27). The best-fitting model was the Hasegawa Kino Yano model with gamma rate heterogeneity. The evolution rate was estimated from the full group of 253 heterochronous sequences within a Bayesian coalescent framework by a Markov chain Monte Carlo (MCMC) method using the BEAST program (13). The Bayesian MCMC method estimates the parameter as the mean of its marginal posterior distribution while simultaneously incorporating uncertainty in the underlying genealogy and other parameters. Evolutionary rates were estimated using both strict and relaxed molecular clocks (uncorrelated lognormal model) as implemented in BEAST (10). Two population genetic models were applied to analyze the data. The first model assumed that the population size was constant over time. The second, the skyline population model, determined the population growth model using the data supplied (14). The unweighted-pair group method using average linkages was used to construct the starting tree. Uncertainty in the estimated parameter values is summarized by the highest posterior density interval that contains 95% of the marginal posterior distribution. The length and number of MCMC chains were chosen so that the effective sample size for each parameter was >100, indicating that parameter space was sufficiently explored.

Currently implemented coalescent methods assume panmictic populations, free of recombination and selection. Substantial selection, recombination, or population subdivision may adversely affect the analysis of temporally spaced viral sequences. The first two assumptions are realistic, since earlier studies showed that RYMV evolves under a marked purifying selection, with few sites under diversifying selection (18). Consequently, positive selection makes a minor contribution to the molecular evolution of the genome. There is no evidence of recombination, either (18). By contrast, the third assumption is not verified, because the diversity of RYMV has a strong spatial basis (1, 39). Therefore, it was checked whether or not population subdivision affected the analysis of temporally spaced viral sequences by comparing the estimates from the group of 253 isolates to estimates from the geographically based groups with minimal population subdivision. The synonymous evolution rate of RYMV was computed using the nucleotides at the third codon position, which are fourfold degenerate at the amino acid level and where most synonymous substitutions occur. This approach was chosen in order to compare the synonymous evolution rate of RYMV to the rates of 50 RNA animal viruses, calculated similarly (26). The evolution rate of the CP gene was extrapolated to the other ORFs and to the full genome, assuming proportionality between the evolution rate and nucleotide diversity. Then the nucleotide diversities of the ORFs and of the full genome were calculated from a set of 22 isolates representative of the geographic and molecular diversity by using DnaSP (37).

Experimental studies. The highly susceptible cultivar Oryza sativa indica cv. IR64, the partially resistant cultivar Oryza sativa japonica cv. Azucena, and the highly resistant cultivars Oryza sativa indica cv. Gigante and Oryza sativa indica cv. Bekarosaka were inoculated with RYMV isolates of different strains. Infection after host change was monitored over time. Isolates were fully sequenced before inoculation and 1 to 6 months later. In all, six experiments were conducted. For each experiment, the number of nucleotide changes between the two stages was counted, and the numbers of synonymous and nonsynonymous substitutions were distinguished. The total number of changes observed experimentally for the six experiments pooled was subsequently compared to that estimated from the evolution rate of the full genome for the same time period.

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 
Overall and synonymous rates of evolution. Geographical origins, strains, numbers of isolates, sampling times, and genetic distances differed widely among the five groups (see Table S1 in the supplemental material). Nevertheless, for each group of isolates, the genetic distance increased with separation time (Fig. 1). No outliers that could substantially alter the substitution rates were detected in any of the groups. The correlation between genetic distance and sampling time was highly significant within each group, except for group III (P = 0.07), which had the shortest sampling interval, and ranged from 0.42 to 0.62 (Fig. 1), indicating that 18% to 38% of the data variance was explained by temporal variation. Whatever the group, the substitution rates were within the range of 3.9 x 10^–4 to 7.7 x 10^–4 nt/site/year.

View larger version (29K): [in this window] [in a new window]

FIG. 1. Genetic distance (percent nucleotide difference) and sampling interval (number of years in absolute value) of each pair of isolates for each of the five groups and the regression line between the two variables. r, estimate of the evolution rate; R, correlation coefficient.

For Bayesian estimates, all 253 sequences were used. Under the strict molecular clock model, the evolution rates were 5.2 x 10^–4 and 5.1 x 10^–4 nt/site/year under the constant-size and Bayesian skyline population genetic models, respectively (Table 1). Under the lognormal relaxed models, the rates were 7.3 x 10^–4 and 7.7 x 10^–4 nt/site/year under the constant-size and Bayesian skyline population genetic models, respectively. The rates were therefore more dependent on the molecular clock model than on the population genetic model. Comparison of the posterior probabilities indicated that the relaxed molecular clock was a better fit than the strict molecular clock, whatever the population model. However, the estimates are close and the confidence intervals overlap. Overall, the range of evolution rates was between 3.9 x 10^–4 and 7.7 x 10^–4 nt/site/year when estimated by linear regression and between 5.1 x 10^–4 and 7.7 x 10^–4 nt/site/year with the Bayesian approach. These ranges are quite similar. A conservative range of the evolution rate of the CP gene of RYMV would be between 4 x 10^–4 and 8 x 10^–4 nt/site/year, a narrow range altogether, i.e., within a 1:2 ratio. Incidentally, our results show that population subdivision does not substantially affect the estimates of the substitution rate, an assumption hardly ever checked in these kinds of analyses.

View this table: [in this window] [in a new window]

TABLE 1. Estimates of the evolution rate of Rice yellow mottle virus by Bayesian analysisa with different molecular clock and population genetic models

Estimates of evolution rates derived from pairs of isolates differing by their sampling dates gave inconsistent results. For instance, isolates CIa and CI12 of the S2/S3 strain in Côte d'Ivoire, collected, respectively, in 1975 and 1991 (see Table S1 in the supplemental material), differed by 4.3% in the 720-nt CP gene. This leads to an estimated evolution rate of 2.7 x 10^–3 nt/site/year. Isolates CI17 and CI104, of the same strain and country, collected in 1975 and 2001, respectively, differed by 3.4%, giving an evolution rate of 1.3 x 10^–3 nt/site/year. Such estimates are highly dependent on the pair of isolates selected and are two to four times greater than estimates derived from the statistical analyses of temporally spaced sequences (4 x 10^–4 to 8 x 10^–4 nt/site/year). This inconsistency reflects the uncertainty about the ancestral diversity of these isolates. It emphasizes that this approach is appropriate for neither plant viruses nor animal viruses or other organisms (11).

Comparison to RNA animal viruses. The RYMV evolution rate of 5.2 x 10^–4 nt/site/year, estimated under the strict molecular clock and constant population models, was compared to the rates of 50 RNA animal viruses calculated with similar assumptions and using broadly similar analytical techniques (26). The evolution rate of RYMV was below the average (7.0 x 10^–4 nt/site/year) but above the distribution median (3.6 x 10^–4 nt/site/year) (Fig. 2). A similar conclusion was reached when the evolution rate of the CP gene (ORF4) was extrapolated to the other ORFs. The nucleotide diversities calculated from 22 isolates representative of the genetic and geographic diversities of RYMV were 0.098 for ORF4, 0.117 for ORF1, 0.057 for ORF2a, and 0.062 for ORF2b. Assuming proportionality between the nucleotide diversity and the substitution rate gave the following estimates: ORF1, 6.1 x 10^–4 nt/site/year; ORF2a, 3.0 x 10^–4 nt/site/year; ORF2b, 3.3 x 10^–4 nt/site/year.

View larger version (36K): [in this window] [in a new window]

FIG. 2. Overall (top) and synonymous (bottom) evolution rates of a range of 50 RNA animal viruses (Jenkins et al., 2002 [26]) and of the RYMV CP gene (solid bars and arrows). Horizontal dashed lines indicate the median and average evolution rates of the 50 RNA animal viruses.

The synonymous evolution rate of RYMV, estimated from the third codon position, was 8.0 x 10^–4 and 11.1 x 10^–4 nt/site/year under the strict and lognormal relaxed molecular clock models, respectively. The evolution rate of 8.0 x 10^–4 nt/site/year is within the range of the evolution rates of 50 RNA animal viruses estimated under the strict molecular clock model (26) and of 49 other species from another study that included more slowly or rapidly evolving viruses (24). The synonymous evolution rate of RYMV was below the average (1.6 x 10^–3 nt/site/year) but above the distribution median (7.1 x 10^–4 nt/site/year) for the 50 RNA animal viruses (Fig. 2). The synonymous-to-nonsynonymous substitution (dS/dN) ratio, estimated from the third codon position versus the first two positions, was 5.9. This value is close to the dS/dN ratio of 5.8 estimated from the diversity of all 253 sequences.

The overall and synonymous rates of evolution of RYMV fell within the range of RNA animal viruses, with a rate higher than those of many animal viruses. It shows that an RNA plant virus can evolve as rapidly as most RNA animal viruses. The substitution rate of RYMV was similar to that of Tomato yellow curl virus (TYLCV), a fast-evolving plant DNA virus, recently calculated by the same Bayesian coalescent method (15). Altogether, the evolution rates of RYMV and TYLCV show that plant viruses can evolve as rapidly as most animal viruses but not as fast as the fastest-evolving animal viruses, such as human immunodeficiency virus type 1 or influenza A virus (24, 26). Evolution rates of RNA plant viruses, like those of RNA animal viruses (24, 26), may encompass a wide range of values. Similar studies should be conducted to determine the range of evolution rates of plant viruses.

Nucleotide changes in experimental studies. The number of changes observed experimentally was compared to that expected from the estimated evolution rate. The number of nucleotide changes observed experimentally on the full genome is quite variable between experiments (Table 2). A large variation in the number of changes was also observed when ORFs 2a and 4 of several other isolates were sequenced 1 to 6 months after inoculation (data not shown). Thus, inferring the rate of evolution from such low and variable numbers of experimental changes is speculative. Furthermore, the total number of changes observed experimentally greatly exceeded that estimated from the evolution rate. The ratio of the nucleotide diversity of the full genome to that of the CP gene for 22 representative isolates of RYMV is 0.8 (see above). Based on this ratio, the evolution rate of the full-length genome would be between 3.2 x 10^–4 and 6.5 x 10^–4 nt/site/year by extrapolation from the evolution rate of the CP gene (4 x 10^–4 to 8 x 10^–4 nt/site/year). Three to 6 nucleotide changes are expected when these rates are applied to the 4,450-nt genome for the length of time of the experiments, a number three to six times lower than the 19 changes observed experimentally (Table 2).

View this table: [in this window] [in a new window]

TABLE 2. Numbers of nucleotide changes observed in experimental studies with several RYMV isolates, rice cultivars, and durations of infection

Several biological features account for these differences. First, the rice seedlings were grown under optimal experimental conditions (light, temperature, nutrients). This is likely to allow a higher rate of virus replication—known to be critical in the accumulation of substitutions (24)—than that in field-grown rice plants, where virus particles survive part of the year in senescent plants or in stubble fields, i.e., with low or no replication, respectively. Second, some substitutions assayed experimentally may represent transient deleterious polymorphisms that will be removed subsequently by purifying selection. Third, in these experiments, RYMV was under a selection pressure (host change), and several nonsynonymous changes emerged (8/19), some of them known to be adaptive (34). Thus, the dS/dN ratio found experimentally was lower than that estimated from field isolates (1.4 versus 5.8, respectively). In nature, the frequency of host changes and, more generally, the extent of adaptation are unknown. However, the high rate of synonymous substitutions, the low dN/dS ratio, and the small number of sites under diversifying selection (18) suggest that RYMV evolution in nature is due to rapid mutational dynamics rather than to the frequent adaptive evolution that is observed experimentally. Taken together, the experimental studies confirm that RYMV can change rapidly over time, support the view that RNA plant viruses—like RNA animal viruses—are measurably evolving populations, and show the need for caution in extrapolating experimental results to field situations.

Plant viruses as measurably evolving populations. Molecular studies have enlarged the scope of plant virus epidemiology (17, 21, 31). Recently, renewed attention has been paid to ecological studies dealing with topics such as plant viruses in wild ecosystems, virus spread from wild to cultivated plants, virus emergence, and host plant domestication or movement (5, 16, 23, 29, 40). The resulting ecological and epidemiological scenarios are related to plant virus history or even prehistory. They raise the recurrent question of the role of agriculture in the emergence and diversification of plant viruses. Indeed, a sound time scale of virus evolution is necessary in order to validate and to calibrate these scenarios. Our results show that the evolution rates of plant viruses can be calculated and consequently that these time scales can be estimated. For instance, the substitution rate of TYLCV was used recently to date the diversification of TYLCV species (15).

 
We are grateful to W. Bakker, L. Bos, I. Buddenhagen, C. Fauquet, S. Fomba, J. L. Notteghem, Y. Pinto, D. Taylor, G. Thottappilly, and R. van der Vlugt for providing information on the dates of collection of some of the early RYMV isolates. We thank R. Froissart, N. Galtier, B. Lafay, B. Moury, O. Pybus, and A. Rodrigo for helpful discussions and M. Thresh, André Fargette, and three anonymous referees for constructive criticisms of the manuscript. We acknowledge the support of C. Tranchant, O. Plechakoya, and the computing department of IRD/Montpellier for providing us access to adequate computational capacities.

This work was supported by the ECOGER action of the French National Program ANR "Ecosphère continentale: processus, modélisation et risques environnementaux," supervised by INRA.

 
* Corresponding author. Mailing address: IRD, BP 64501, 34394 Montpellier cedex 5, France. Phone: 33-467416227. Fax: 33-467416330. E-mail: Denis.Fargette{at}mpl.ird.fr

Published ahead of print on 16 January 2008.

Supplemental material for this article may be found at http://jvi.asm.org/.

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 

1
1. Abubakar, Z., F. Ali, A. Pinel, O. Traoré, P. N'Guessan, J. Notteghem, F. Kimmins, G. Konaté, and D. Fargette. 2003. Phylogeography of Rice yellow mottle virus in Africa. J. Gen. Virol. 84:733-743.[Abstract/Free Full Text]
2
2. Allarangaye, M. D., O. Traoré, E. V. S. Traoré, R. J. Millogo, S. Guinko, and G. Konaté. 2007. Host range of Rice yellow mottle virus in Sudano-Sahelian savannahs. Pakistan J. Biol. Sci. 10:1414-1421.
3
3. Bakker, W. 1974. Characterization and ecological aspects of rice yellow mottle virus in Kenya. Agricultural research report no. 829. Wageningen Agricultural University, Wageningen, The Netherlands.
4
4. Blok, J., A. Mackenzie, P. Guy, and A. Gibbs. 1987. Nucleotide sequence comparisons of turnip yellow mosaic virus isolates from Australia and Europe. Arch. Virol. 97:283-295.[CrossRef][Medline]
5
5. Borer, E., P. Hosseini, E. Seablom, and A. Dobson. 2007. Pathogen-induced reversal of native dominance in a grassland community. Proc. Natl. Acad. Sci. USA 104:5473-5478.[Abstract/Free Full Text]
6
6. Childs, J. E., J. S. Mackenzie, and J. A. Richt (ed.). 2007. Current topics in microbiology and immunology, vol. 315. Wildlife and emerging zoonotic diseases: the biology, circumstances and consequences of cross-species transmission. Springer, London, United Kingdom.
7
7. Domingo, E., C. Escarmis, N. Sevilla, A. Moya, S. Elena, J. Quer, I. Novella, and J. Holland. 1996. Basic concepts in RNA evolution. FASEB J. 10:859-864.[Abstract]
8
8. Domingo, E., and J. Holland. 1997. RNA virus mutations and fitness for survival. Annu. Rev. Microbiol. 51:151-178.[CrossRef][Medline]
9
9. Drake, J., and J. Holland. 1999. Mutation rates among RNA viruses. Proc. Natl. Acad. Sci. USA 96:13910-13913.[Abstract/Free Full Text]
10
10. Drummond, A., S. Ho, M. Phillips, and A. Rambaut. 2006. Relaxed phylogenetics and dating with confidence. PLoS Biol. 4:e88.[CrossRef][Medline]
11
11. Drummond, A., O. G. Pybus, and A. Rambaut. 2003. Inference of viral evolutionary rates from molecular sequences. Adv. Parasitol. 54:331-358.[CrossRef][Medline]
12
12. Drummond, A., and A. Rodrigo. 2000. Reconstructing genealogies of serial samples under the assumption of a molecular clock using serial-sample UPGMA. Mol. Biol. Evol. 17:1807-1815.[Abstract/Free Full Text]
13
13. Drummond, A., and A. Rambaut. 2007. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 7:214.[CrossRef][Medline]
14
14. Drummond, A., A. Rambaut, B. Shapiro, and O. Pybus. 2005. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol. Biol. Evol. 22:1185-1192.[Abstract/Free Full Text]
15
15. Duffy, S., and E. C. Holmes. 31 October 2007. Phylogenetic evidence for rapid rates of molecular evolution in the single-stranded DNA begomovirus Tomato yellow leaf curl virus. J. Virol. doi:10.1128/JVI.01929-07. (Subsequently published, J. Virol. 82:957-965, 2008.)[Abstract/Free Full Text]
16
16. Dwyer, G. I., M. J. Gibbs, A. J. Gibbs, and R. A. C. Jones. 2007. Wheat streak mosaic virus in Australia: relationships to isolates from the Pacific Northwest of the USA and its dispersion via seed transmission. Plant Dis. 91:164-170.[CrossRef]
17
17. Fargette, D., G. Konaté, C. Fauquet, E. Muller, M. Peterschmitt, and J. M. Thresh. 2006. Molecular ecology and emergence of tropical plant viruses. Annu. Rev. Phytopathol. 44:235-260.[CrossRef][Medline]
18
18. Fargette, D., A. Pinel, Z. Abubakar, O. Traoré, C. Brugidou, S. Fatogoma, E. Hébrard, M. Choisy, S. Yacouba, C. Fauquet, and G. Konaté. 2004. Inferring the evolutionary history of Rice yellow mottle virus from genomic, phylogenetic and phylogeographic studies. J. Virol. 78:3252-3261.[Abstract/Free Full Text]
19
19. Fauquet, C. M., and J. C. Thouvenel. 1977. Isolation of the rice yellow mottle virus in Ivory Coast. Plant Dis. Reporter 61:443-446.
20
20. Fraile, A., F. Escriu, M. Aranda, J. Malpica, A. Gibbs, and F. García-Arenal. 1997. A century of tobamovirus evolution in an Australian population of Nicotiana glauca. J. Virol. 71:8316-8320.[Abstract]
21
21. García-Arenal, F., A. Fraile, and J. Malpica. 2001. Variability and genetic structure of plant virus populations. Annu. Rev. Phytopathol. 39:157-186.[CrossRef][Medline]
22
22. Gibbs, A., P. Keese, M. Gibbs, and F. García-Arenal. 1999. Plant virus evolution: past, present and future, p 263-285. In E. Domingo, R. Webster, and J. Holland (ed.), Origin and evolution of viruses. Academic Press, San Diego, CA.
23
23. Godfree, R., P. Thrall, and A. Young. 2007. Enemy release after introduction of disease-resistant genotypes into plant-pathogen systems. Proc. Natl. Acad. Sci. USA 104:2756-2760.[Abstract/Free Full Text]
24
24. Hanada, K., Y. Suzuki, and T. Gojobori. 2004. A large variation in the rates of synonymous substitution for RNA viruses and its relationships to a diversity of viral infection and transmission modes. Mol. Biol. Evol. 21:1074-1080.[Abstract/Free Full Text]
25
25. Hébrard, E., A. Pinel-Galzi, A. Bersoult, C. Siré, and D. Fargette. 2006. Emergence of a resistance-breaking isolate of Rice yellow mottle virus during serial inoculations is due to a single substitution in the genome-linked viral protein VPg. J. Gen. Virol. 87:1369-1373.[Abstract/Free Full Text]
26
26. Jenkins, G. M., A. Rambaut, O. G. Pybus, and E. C. Holmes. 2002. Rates of molecular evolution in RNA viruses: a quantitative phylogenetic analysis. J. Mol. Evol. 54:156-165.[CrossRef][Medline]
27
27. Kosakovsky Pond, S. L., S. D. Frost, and S. V. Muse. 2005. HyPhy: hypothesis testing using phylogenies. Bioinformatics 21:676-679.[Abstract/Free Full Text]
28
28. Kouassi, N. K., P. N'Guessan, L. Albar, C. M. Fauquet, and C. Brugidou. 2005. Distribution and characterization of Rice yellow mottle virus: a threat to African farmers. Plant Dis. 89:124-133.[CrossRef]
29
29. Malmstrom, C. M., R. Shu, E. W. Linton, L. A. Newton, and M. A. Cook. 2007. Barley yellow dwarf viruses (BYDVs) preserved in herbarium specimens illuminate historical disease ecology of invasive and native grasses. J. Ecol. 95:1153-1166.[CrossRef]
30
30. Mantel, N. 1967. The detection of disease clustering and a generalized regression approach. Cancer Res. 27:209-220.[Abstract/Free Full Text]
31
31. Moury, B., C. Desbiez, M. Jacquemond, and H. Lecoq. 2006. Genetic diversity of plant virus populations: towards hypothesis testing in molecular epidemiology. Adv. Virus Res. 67:49-87.[CrossRef][Medline]
32
32. N'Guessan, P., A. Pinel, M. Caruana, R. Frutos, A. Sy, A. Ghesquière, and D. Fargette. 2000. Evidence of the presence of two serotypes of rice yellow mottle sobemovirus in Côte d'Ivoire. Eur. J. Plant Pathol. 106:167-178.[CrossRef]
33
33. Pinel, A., P. N'Guessan, M. Bousalem, and D. Fargette. 2000. Molecular variability of geographically distinct isolates of Rice yellow mottle virus in Africa. Arch. Virol. 145:1621-1638.[CrossRef][Medline]
34
34. Pinel-Galzi, A., M. Rakotomalala, E. Sangu, F. Sorho, Z. Kanyeka, O. Traoré, D. Sérémé, N. Poulicard, Y. Rabenantoandro, Y. Séré, G. Konaté, A. Ghesquière, E. Hébrard, and D. Fargette. 2007. Theme and variations in the evolutionary pathways to virulence of an RNA plant virus species. PLoS Pathog. 3:e180.[CrossRef][Medline]
35
35. Posada, D., and K. Crandall. 1998. MODELTEST: testing the model of DNA substitution. Bioinformatics 14:817-818.[Abstract/Free Full Text]
36
36. Rodrigo, A., G. Ewing, and A. Drummond. 2007. The evolutionary analysis of measurably evolving populations using serially sampled gene sequences, p 30-62. In O. Gascuel and M. Steel (ed.), Reconstructing evolution: new mathematical and computational advances. Oxford University Press, Oxford, United Kingdom.
37
37. Rozas, J., J. Sánchez-DelBarrio, X. Messeguer, and R. Rozas. 2003. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19:2496-2497.[Abstract/Free Full Text]
38
38. Stenger, D., D. Seifers, and R. French. 2002. Patterns of polymorphism in Wheat streak mosaic virus: sequence space explored by a clade of closely related viral genotypes rivals that between the most divergent strains. Virology 302:58-70.[CrossRef][Medline]
39
39. Traoré, O., F. Sorho, A. Pinel, Z. Abubakar, O. Banwo, J. Maley, E. Hébrard, S. Winter, Y. Séré, G. Konaté, and D. Fargette. 2005. Processes of diversification and dispersion of Rice yellow mottle virus inferred from large-scale and high-resolution phylogeographic studies. Mol. Ecol. 14:2097-2110.[CrossRef][Medline]
40
40. Wren, J. D., M. J. Roossinck, R. S. Nelson, K. Scheets, M. W. Palmer, and U. Melcher. 2006. Plant virus biodiversity and ecology. PLoS Biol. 4:e80.[CrossRef][Medline]

---

Journal of Virology, April 2008, p. 3584-3589, Vol. 82, No. 7
0022-538X/08/$08.00+0     doi:10.1128/JVI.02506-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

This article has been cited by other articles:

• Sanjuan, R., Agudelo-Romero, P., Elena, S. F. (2009). Upper-limit mutation rate estimation for a plant RNA virus. Biol Lett 5: 394-396 [Abstract] [Full Text]  
• Duffy, S., Holmes, E. C. (2009). Validation of high rates of nucleotide substitution in geminiviruses: phylogenetic evidence from East African cassava mosaic viruses. J. Gen. Virol. 90: 1539-1547 [Abstract] [Full Text]  

This Article Abstract Full Text (PDF) Supplemental material Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Fargette, D. Articles by Konaté, G. Search for Related Content PubMed PubMed Citation Articles by Fargette, D. Articles by Konaté, G.