- Instituto Valenciano de Investigaciones Agrarias, Moncada, Valencia, Spain
RNA viruses have a great potential for genetic variation, rapid evolution and adaptation. Characterization of the genetic variation of viral populations provides relevant information on the processes involved in virus evolution and epidemiology and it is crucial for designing reliable diagnostic tools and developing efficient and durable disease control strategies. Here we performed an updated analysis of sequences available in Genbank and reviewed present knowledge on the genetic variability and evolutionary processes of viruses of the family Closteroviridae. Several factors have shaped the genetic structure and diversity of closteroviruses. (I) A strong negative selection seems to be responsible for the high genetic stability in space and time for some viruses. (2) Long distance migration, probably by human transport of infected propagative plant material, have caused that genetically similar virus isolates are found in distant geographical regions. (3) Recombination between divergent sequence variants have generated new genotypes and plays an important role for the evolution of some viruses of the family Closteroviridae. (4) Interaction between virus strains or between different viruses in mixed infections may alter accumulation of certain strains. (5) Host change or virus transmission by insect vectors induced changes in the viral population structure due to positive selection of sequence variants with higher fitness for host-virus or vector-virus interaction (adaptation) or by genetic drift due to random selection of sequence variants during the population bottleneck associated to the transmission process.
Introduction
There are five basic mechanisms determining the genetic structure and evolution of biological populations: mutation, recombination, natural selection, genetic drift, and migration (Moya et al., 2004). RNA viruses have a great potential for high genetic variability, rapid evolution and adaptation to new conditions and environments due to their rapid replication, generation of very large populations, and high mutation rates (at least 105 times higher than those of their hosts) as a consequence of the lack of proofreading activity of RNA polymerases (Holland et al., 1982; García-Arenal and Fraile, 2011). In many of these viruses, genome recombination and/or reassortment of genomic segments (pseudorecombination) between divergent virus strains or virus species increase genetic variability and accelerate evolution (Chare and Holmes, 2004; Nagy, 2008). The genetic variation generated by mutation and recombination is limited and structured by the other three evolutionary forces: natural selection, genetic drift, and gene flow (Roossinck, 2003; Moya et al., 2004). Natural selection is a directional process by which variants that are fittest in a certain environment will increase their frequency in the population (positive or adaptive selection), whereas variants less fit will decrease their frequency (negative or purifying selection), this process being determined by numerous specific interactions of viruses with their plant hosts (Schneider and Roossinck, 2001), vectors (Power, 2000; Chare and Holmes, 2004) and even with other viruses co-infecting the same plant. Genetic drift consists of stochastic changes in allele frequencies in a finite population due to the random sampling of genes at reproduction (Moya et al., 2004). This supposes a reduction of the genetic variability and fixation of selectively neutral variants, and it has an important effect during severe and rapid reduction in population size produced by population bottlenecks or founder events (Ali and Roossinck, 2008). Genetic drift can occur in different events of the virus life cycle such as virus movement between plant cells (Sacristan et al., 2003; Li and Roossinck, 2004), transmission between plants by vectors (Ali et al., 2006; Betancourt et al., 2008) and interaction between coinfecting viruses (Fraile et al., 1997). Finally, migration (gene flow) among distinct geographical areas, plants or different parts of the same plant is an important factor shaping the genetic structure of viral populations. High migration favors genetic uniformity between populations and thus decreases the global genetic diversity (Moya et al., 2004). While mutation and recombination are intrinsic of the virus genome and its replication and expression systems, natural selection, genetic drift and gene flow are affected by virus biology (e.g., host type and range, means and extent of dispersal), environmental conditions, and population parameters (e.g., population size and history of population bottlenecks).
The study of genetic variability and the evolutionary mechanisms related to the different aspects of the virus biology is crucial to understand virus epidemiology and emergence (Grenfell et al., 2004), designing specific diagnostic tools, and developing efficient and durable strategies of disease control (García-Arenal and McDonald, 2003; Acosta-Leal et al., 2011). Several reviews on evolution of the family Closteroviridae have been published, which were focused on macroevolution and taxonomy (Dolja et al., 1994; Agranovsky, 1996; Karasev, 2000; Dolja et al., 2006).
The family Closteroviridae is composed of viruses characterized by their long (up to 2000 nm) and flexuous, non-enveloped, polar, virions with two coat proteins, the major (CP), covering most of the genomic RNA, and the minor (CPm) located to one of the virion ends (Agranovsky et al., 1995; Febres et al., 1996; Tian et al., 1999). Its members have the largest genomes of all positive sense RNA plant viruses (up to 20 kb). Although the number and relative position of open reading frames (ORFs) vary between species, there is a common genome organization. ORFs 1a and 1b encode replication-related proteins, with protease, methyl-transferase, helicase, and RNA-dependent RNA polymerase conserved domains. Downstream ORFS include a conserved five-gene module encoding a small hydrophobic protein with affinity to cell membranes, a homolog of the plant heat shock proteins HSP70 (HSP70h), a ~60 kDa protein with a diverged coat protein motif, the CP and the CPm. The functions postulated for the HSP70h are: cell-to-cell movement, involvement in the assembly of multisubunits complexes for genome replication and/or subgenomic RNA synthesis and assembly of virus particles, whereas the ~60 kDa protein is required for incorporation of both HSP70h and CPm to the particle tail (Tian et al., 1999; Satyanarayana et al., 2000; Alzhanova et al., 2001). The genome expression strategy is based on: (I) proteolytic processing of the polyprotein encoded by ORF 1a, (II) +1 ribosomal frameshift for the expression of ORF1b, and (III) expression of the downstream ORFs via the formation of 3′ co-terminal subgenomic RNAs. Presently, there are three genera in the family (Martelli et al., 2011): Ampelovirus, Crinivirus and Closterovirus, although a new genus named Velarivirus has been proposed (Al Rwahnih et al., 2012). The characteristics of the accepted genera are:
- Ampelovirus. Mealybug-transmitted, monopartite genome, and the CPm gene is located downstream of CP gene or lacking in some species. Viruses studied here are: Grapevine leafroll-associated virus 1 (GLRaV-1), GLRaV-3, GLRaV-5, and GLRaV-11 (tentative member), and Pineapple mealybug wilt-associated virus 1 (PMWaV-1).
- Crinivirus. Whitefly-transmitted, bipartite or tripartite genome, and the CP gene is located upstream of the CPm gene. Viruses analyzed here are: Blackberry yellow vein-associated virus (BYVaV), Cucurbit yellow stunting disorder virus (CYSDV), Potato yellow vein virus (PYVV), Sweet potato chlorotic stunt virus (SPCSV), Tomato chlorosis virus (ToCV), and Tomato infectious chlorosis virus (TICV).
- Closterovirus. Mostly aphid-transmitted, monopartite genome, and the CP gene located downstream of the CPm gene. Viruses studied here are: Citrus tristeza virus (CTV) and Grapevine leafroll-associated virus 2 (GLRaV-2).
In this work we performed an updated analysis of the genetic variation of viruses in the family Closteroviridae by analysis of the coat protein genes using nucleotide sequences retrieved from Genbank and present an updated review on the genetic variability and evolutionary processes of the viral populations of members of the family Closteroviridae.
Materials and Methods
Alignement of Nucleotide Sequences
Nucleotide sequences from worldwide isolates of members of the family Closteroviridae were retrieved from GenBank (http://www.ncbi.nlm.nih.gov and http://www.dpvweb.net). The coat protein genes (CP or CPm) were selected because they are present in all viruses and because it is the genomic region for which more sequences are available. Only those viruses with sequences of five or more different isolates were analyzed (Table 1). Multiple alignment was performed with the algorithm CLUSTAL W (Larkin et al., 2007) implemented in the program MEGA 5.05 (Tamura et al., 2011).
Table 1. Genetic diversity and population genetic parameters of viruses of the family Closteroviridae.
Analysis of Nucleotide Sequences
The program MEGA 5.05 was used for: (1) inference of phylogenetic relationships between isolates of each viral species by the neighbor-joining method (Saitou and Nei, 1987), (II) estimation of nucleotide distances between sequence pairs and diversities (mean nucleotide distances) using Kimura-two-parameter as the nucleotide substitution model (Kimura, 1980), and (III) estimation of the ratio between non-synonymous and synonymous substitution (N/S) by the Pamilo-Bianchi-Li method (Pamilo and Bianchi, 1993) to study the role of natural selection at the protein level. N/S ≈ 0 indicates neutral evolution, N/S < 1 negative or purifying selection and N/S > 1 positive or adaptive selection.
The program DnaSP 5.10 (Librado and Rozas, 2009) was used to assess genetic differentiation and the gene flow level between different countries or geographic areas with the statistic Fst (Weir and Cockerham, 1984). Fst can take values from 0, no genetic differentiation and complete gene flow, to 1, complete genetic differentiation as a consequence of null gene flow. Only countries or geographical areas with more than four isolates of the virus analyzed were taken into account.
Recombination between isolates of the same virus was analyzed with the program RDP3 (Martin et al., 2010) that incorporates the recombination-detecting algorithms GENECONV (Padidam et al., 1999), BOOTSCAN (Salminen et al., 1995; Martin et al., 2005), MAXCHI (Smith, 1992; Posada and Crandall, 2001), CHIMAERA (Posada and Crandall, 2001), SISCAN (Gibbs et al., 2000), 3SEQ (Boni et al., 2007), and RDP (Martin and Rybicki, 2000), using their default parameter values. Only those events recognized by at least four different algorithms were accepted as evidence for recombination. The effect of recombination was taken into account during analysis of selection.
Results
Genetic Variation Between Virus Isolates
Figure 1 shows the phylogenetic relationships between isolates of each virus species with branch length indicating genetic distances. For each virus, isolates were classified into genetic groups considered as those clades with all isolates having nucleotide distances higher than 0.1 with respect to all isolates of the other clades. These groups are indicated in gray boxes. Table 1 shows the nucleotide diversities and other population genetic parameters.
Figure 1. Neighbor-joining phylogenetic trees of the coat protein genes of different viruses of the family Closteroviridae. Bootstrap values higher than 0.75 are indicated with asterisks. Two-letter codes indicate countries (defined in ISO 3166-1) and/or US states. Numbers preceding the codes indicate the number of isolates analyzed from each country. For each virus, genetic groups are indicated in gray boxes with Roman numerals and include virus isolates having nucleotide distances higher than 0.1 with all isolates from other clades. Subgroups or clades in CTV are indicated in darker boxes with Arabic numbers. (A) Genus Ampelovirus, (B) Genus Crinivirus, and (C) genus Closterovirus.
In the genus Ampelovirus, four out of the six viruses studied comprised three to five genetic groups with distances between isolates of up to 0.391 (Figure 1A and Table 1). GLRaV-1 sequences (Alabi et al., 2011) formed three genetic groups with isolates from: (I) Canada, (II) Iran, and (III) Brazil, Poland and South Africa. GLRaV-3 sequences (Turturo et al., 2005; Fajardo et al., 2007; Chooi et al., 2009; Fuchs et al., 2009; Jooste et al., 2010; Gouveia et al., 2011; Sharma et al., 2011; Wang et al., 2011; Bester et al., 2012; Farooq et al., 2012) also fell into three groups including isolates from: (I) South Africa, (II) California, Portugal, South Africa and China, and (III) California, Brazil, Chile, Portugal, Poland, South Africa, China and Taiwan, with group I diverging much more than groups II and III, and with a maximum nucleotide distance between isolates of 0.256 (Table 1). GLRaV-4, in spite of having CP sequences of only five isolates (Saldarelli et al., 2006) was composed of three genetic groups with isolates from: (I) Israel, Turkey and USA, (II) Japan, and (III) Japan, with group I diverging more than groups II and III and with a maximum nucleotide distance of 0.391. This was the most variable member of the family with a nucleotide diversity value of 0.254 (Table 1). GLRaV-5, with most isolates from Portugal and Argentina, showed moderate genetic variability and all isolates clustered in a single genetic group. GLRaV-11, in spite of having sequences of only seven isolates (six from Argentina and one from Greece), was composed of five genetic groups: four from Argentina and one from Greece, with a nucleotide diversity value of 0.109 (Table 1). Finally, for the pineapple-infecting PMWaV-1, the CP sequences available include four isolates from Cuba and one from Taiwan which were similar and formed a unique genetic group.
The genus Crinivirus had the lowest genetic variability of the three genera, with nucleotide diversity values lower than 0.033 (Table 1). Four viruses were composed of isolates clustered in a unique genetic group and two viruses had isolates clustered in two divergent genetic groups (Figure 1B), with isolates within each group having very low variability (diversity below 0.015). The blackberry-infecting BYVaV isolates, all from USA, presented very low genetic variation. The cucurbit-infecting CYSDV was composed of two genetic groups located in: (I) Middle East (Iran and Arabia) and (II) Mediterranean Basin (South Europe, North Africa and Near East), North America and China. Each group had a very low nucleotide diversity (Table 1; Rubio et al., 1999, 2001a; Marco and Aranda, 2005; Sweiss et al., 2007). The potato-infecting PVYV was restricted to South America and presented very low genetic variation. The sweet potato-infecting SPCSV had two very divergent genetic groups with a nucleotide distance of up to 0.379 (Table 1): (I) located predominantly in Eastern Africa, but also with isolates in Peru and China; and (II) located in Western Africa. In our analysis the Western Africa group was represented by only one isolate from the Canary Islands, but analysis of the HSP70h gene of several isolates from Nigeria showed that they were also part of the Western Africa group which diverged from the East Africa group (Fenby et al., 2002). Finally, the tomato-infecting TICV and ToCV showed very low genetic variation in spite of their wide distribution. TICV isolates were from USA, Europe and Japan; and ToCV isolates from Europe, North America, Mediterranean Basin, Africa (Mauritius) and Japan.
In the genus Closterovirus, the grape-infecting GLRaV-2 showed high genetic variation (Bertazzon et al., 2010b; Jarugula et al., 2010) like the grape-infecting ampeloviruses, with a maximum nucleotide distance between isolates of 0.375 and a nucleotide diversity of 0.128 (Table 1). GLRaV-2 sequences were classified into six genetic groups including isolates from: (I) China, USA, Brazil, Italy, Poland, and South Africa, (II) Italy and South Africa, (III) USA, Italy, and Brazil, (IV) Italy and France, (V) USA, and (VI) Italy (Figure 1C). Finally, the citrus-infecting CTV, the best studied virus in the family, with CP sequences of almost 600 worldwide isolates (Albiach-Martí et al., 2000; Rubio et al., 2001b; Alavi et al., 2005; Roy et al., 2005; Sorrentino et al., 2005; Papayiannis et al., 2007; Iglesias et al., 2008; Jiang et al., 2008; Černi et al., 2009; Herrera-Isidrón et al., 2009; Oliveros-Garay et al., 2009; Fisher et al., 2010; Harper et al., 2010; Matos et al., 2013), showed a moderate genetic variability. All isolates fell into a single genetic group, albeit this could be divided into five to seven subgroups (Figure 1C), some of them very homogeneous (isolates with almost identical sequences). There was no association between the subgroups and the geographic origin of CTV isolates. In our analysis, all biologically characterized isolates in the subgroup 1a were severe and induced stem pitting in sweet orange and/or grapefruit, whereas those characterized in the subgroup two were mild, inducing only weak symptoms in Mexican lime. Other subgroups included mild and severe isolates with no association between symptoms and genetic distance between isolates being observed. For example, about 43, 70, 40, 25, and 19% of the biologically characterized isolates in subgroups 1b, 3b, 3a, 4, and 5a, respectively, were mild, and 57, 30, 60, 75, and 81%, respectively, were severe and incited stem pitting in grapefruit or sweet orange (Alavi et al., 2005; Papayiannis et al., 2007; Černi et al., 2008; Harper et al., 2009; Nolasco et al., 2009; Biswas et al., 2012; Hančević et al., 2013).
Recombination
Most of the analyzed viruses showed no recombination in the CP (Table 2). GLRaV-3 had 14 isolates out of 191 with a recombinant CP which resulted from three different recombination events involving different parental sequences and recombination sites (approximate nucleotide positions 300, 450, and 500 of the CP gene). All recombinations involved isolates from China, except one that had a parental sequence from Chile. GLRaV-11 had four isolates out 15 with the same recombination at nucleotide position 500 involving the same parental sequences. CTV comprised 63 out of 577 isolates with a recombinant CP resulting from two recombination events with respect to the parental sequences and recombination sites (Table 2). CTV isolates from South America, North America, Africa, and Asia showed the same recombination indicating that this corresponds to an ancient event that occurred before the recombinants spread worldwide.
Table 2. Recombination events detected in the coat protein genes of viruses of the family Closteroviridae.
Natural Selection
The ratio of non-synonymous to synonymous substitutions (N/S values in Table 1), were low indicating functional or structural constraints for amino acid changes. The selection pressure was particularly strong for CYSDV, with N/S = 0.008 and composed of two genetic groups separated by a genetic distance of about 0.1. Indeed only 0.2% of the nucleotide changes between the two groups produced amino acid changes. In some cases the number of synonymous substitutions was very low (not shown) suggesting also constraints for nucleotide changes that could affect thermodynamic stability of RNA, codon usage bias for translation efficiency, secondary structure, activation of gene silencing, or RNA-RNA or RNA-protein interactions (Cuevas et al., 2012).
Migration (Gene Flow)
The presence in one geographic region of a diverse virus population with some isolates similar to those of other regions is usually indicative of a possible dispersion center for that virus (Bateson et al., 2002). Our analyses showed that GLRaV-3 had isolates from South Africa in the three genetic groups making plausible to consider South Africa or nearby areas as a dispersion center of this virus. However, this could be a secondary dispersion center since South Africa is far from the host (genus Vitis) origin region and likely initial dispersion area. GLRaV-11 had four genetic groups with isolates from Argentina and one group with one isolate from Greece, suggesting that Argentina could be a dispersion center for this virus. Finally, GLRaV-2 had six genetic groups with isolates from Italy distributed in five groups and isolates from USA in three groups which pointed at both countries as possible dispersion zones. Another explanation is that the locations with divergent isolates received them via importation of infected grapevine material from different areas. In most cases, the phylogenetic relationships among viral isolates were not geographically structured and often isolates from distant regions were genetically very close (Figure 1), indicating long distance movement probably by international traffic of virus-infected plant material (Rubio et al., 2001b; Angelini et al., 2004; Alabi et al., 2011).
To estimate the degree of virus migration or gene flow, the statistic Fst was calculated (see above). The global gene flow was high for GLRaV-3, moderate for GRLaV-5 and GRLaV-2, and low for criniviruses and CTV (Table 1). A more detailed analysis was performed by comparing virus subpopulations of different areas two by two (not shown). GRLaV-1 subpopulations from California, Washington, and New York (Alabi et al., 2011) had an infrequent genetic exchange. CYSDV has a null gene flow between the Middle East and the rest of the world, but the gene flow was very high between Spain, Near East, and North America (Fst = 0.000), albeit this could be due to a unique migration event given the genetic stability of this virus. SPCSV showed also a high geographical structure between the two genetic groups, but availability of only one isolate in one group precluded gene flow analysis between both groups. In group I, the comparison of three neighboring countries of East Africa showed a high gene flow between Uganda and Kenya (Fst = 0.027) but very low between them and Tanzania (Fst ≈ 0.450). In spite of its low genetic diversity, ToCV showed a low gene flow between Europe and North America (Fst = 0.560). GRLaV-2 with a moderate global flow showed a puzzling migration pattern which did not correspond to geographic proximity. Thus, extensive gene flow occurred between Poland and China and between Italy, Brazil and USA but this was very low between Italy and Poland. Finally, CTV, the best documented virus with isolates from 47 countries showed different degrees of gene flow which were not correlated with geographic distance. Some subpopulations had a high gene flow (Fst < 0.100), e.g., California, Mexico, Spain, Italy, and Portugal, or Brazil, Angola, China, and Portugal, whereas other subpopulations were almost isolated (Fst > 0.300) such as Cuba (although with a moderate gene flow, Fst ≈ 0.150, with Mexico and Guatemala), Trinidad and Tobago, and Argentina. Several reports showed that genetically and biologically divergent isolates of CTV have been introduced one or several times in Iran, Sicily (Italy), Cyprus, and Dominican Republic (Alavi et al., 2005; Davino et al., 2005; Papayiannis et al., 2007; Matos et al., 2013).
Discussion
The genetic variability for each virus analyzed was different, although it must be taken into account that in some cases the number of sequences available was low and/or these were from a few geographical locations. The most variable viruses were the grapevine-infecting ampeloviruses GRLaV-1, 4, and 11 and the closterovirus GLRLaV-2, whereas criniviruses showed very low genetic variation. In general, the viruses of the family Closteroviridae showed a low genetic diversity or were comprised of genetic groups with very low within-group nucleotide diversity, as found in many other plant viruses (García-Arenal et al., 2001), due to strong negative selection. The existence of a few genetically homogeneous genetic groups suggests that the sequence space of these viruses may be restricted to a few narrow peaks in the adaptive or fitness landscape (Wright, 1932). A high level of covariation at molecular level (the coordinated change of certain nucleotides in response to the change of other nucleotides to maintain biologically relevant structures and functions) could explain this discontinuous adaptive landscape (Gultyaev et al., 2000). Reduced fitness of chimeras between CTV strains from different genetic subgroups occupying distinct adaptive peaks (Satyanarayana et al., 1999) support this notion. The temporal genetic stability reported for CTV (Albiach-Martí et al., 2000; Rubio et al., 2001b; Silva et al., 2012) and CYSDV (Marco and Aranda, 2005; Rubio et al., 2001a) supports the existence of a strong negative selection. There have been many studies trying to associate genetic relationships and the pathogenicity characteristics of CTV variants (Sambade et al., 2003; Hilf et al., 2005; Nolasco et al., 2009). Although sequence of the CP gene does not appear associated with pathogenicity characteristics in many subgroups, separation between mild and severe stem pitting isolates found in other genomic regions and analysis of phylogenetic networks suggest that pathogenicity has been an important evolutionary force in CTV populations (Martin et al., 2009).
Viruses of the family of Closteroviridae are transmitted by insect vectors which favor mixed infections between different viruses or strains of the same virus. Mixed infections may have important evolutionary implications since they can affect the within-isolate population of virus variants (quasispecies) and allow interaction and/or recombination between different virus entities which can affect to pathogenicity and adaptability.
Recombination has played an important role in the evolutionary history of the family Closteroviridae. It has been postulated that the ancestor of this family was a smaller filamentous virus composed of three genes encoding replication-associated proteins, the p6-like movement protein and a single coat protein (Dolja et al., 2006). During evolution new genes were incorporated to the genome by two means: (i) recombination with cellular mRNAs, e.g., HSP70h, or with RNAs of other viruses, e.g., the leader proteinase (L-Pro) from potyviruses, and (ii) unequal recombination between two genomic RNA copies (or involving subgenomic RNAs) of the same virus which produced gene duplication, e.g., the genes encoding the ~60 kD and CPm proteins evolved after being generated as duplicates of the CP gene. Recent cases of recombination-mediated gene gain have occurred in the criniviruses SPCSV and Beet pseudoyellows virus (BPYV), both with isolates differing in the number of genes (Tzanetakis and Martin, 2004; Cuellar et al., 2008). Unequal recombination also generated the multipartite genome of criniviruses.
Our analyses indicated homologous recombination in the coat protein gene between divergent genetic variants of GLRaV-3, GLRVaV-11, and CTV. It is clear that more recombination events would be found if additional genomic regions were analyzed. Homologous recombination between divergent isolates of the same virus has been reported for CTV (Rubio et al., 2001b; Sambade et al., 2003; Hilf, 2005; Roy et al., 2005; Vives et al., 2005; Weng et al., 2007; Gomes et al., 2008; Martin et al., 2009; Roy and Brlansky, 2009, 2010; Harper et al., 2010; Melzer et al., 2010; Scott et al., 2013), GLRaV-2 (Alabi et al., 2011) and Raspberry leaf mottle virus (McGavin and MacFarlane, 2010); the crinivirus BYVaV (Poudel et al., 2012) and the ampeloviruses GLRaV-3 (Turturo et al., 2005), GLRaV-4 (Thompson et al., 2012) and GLRaV-5 (Turturo et al., 2005; Farooq et al., 2012). Moreover, phylogenetic network analysis showed that homologous recombination must be an important evolutionary force for CTV (Martin et al., 2009). Population analyses showed CTV isolates containing a heterogeneous population of diverged virus strains and recombinants at low frequency (Vives et al., 2005; Weng et al., 2007; Scott et al., 2013), suggesting that these recombinants did not have a selective advantage (more fitness) with respect to the parental sequences. In some cases artificial chimeras of two genetically and biologically divergent CTV isolates failed to infect citrus (Satyanarayana et al., 1999) suggesting that only some recombinants are viable. Also, homologous recombination seems to have occurred between CTV and another unknown closterovirus given the unusual disparity in the divergence of CTV isolates between the two halves of CTV genome: ~0.1 for the 3'-moiety and ~0.3 for the 5'-moiety (Mawassi et al., 1996; Vives et al., 1999).
Unequal recombination seems to occur frequently during replication as evidenced by the large number of defective RNAs (D RNAs) detected in the closterovirus CTV (Mawassi et al., 1995; Ayllón et al., 1999a; Mawassi et al., 2000; Che et al., 2002, 2003), the criniviruses Lettuce infectious yellows virus (LIYV) (Rubio et al., 2000b, 2002), PYVV (Livieratos et al., 2004; Eliasco et al., 2006) and SPCSV (Kreuze et al., 2002; Cuellar et al., 2008), and the ampelovirus GLRaV-3 (Ling et al., 1998). D RNAs are deletion forms of virus genomic RNAs that retain the replication signals but require the parental virus for replication. The existence of direct repeats, secondary structure or AT-rich regions at the junction site of some D RNAs suggested template-switching as a plausible mechanism for recombination (Ayllón et al., 1999a; Rubio et al., 2000b). In this model the D RNAs are generated by a translocation event in which the polymerase, together with the nascent RNA strand, falls off the template strand probably at regions of secondary structure and RNA synthesis reinitiates at a different site with identical or similar nucleotide sequence to the jumping site. Also a weak base pairing in A/U-rich regions can facilitate the release and/or the re-annealing of the incomplete nascent RNA by formation of a temporary bubble (Nagy and Bujarski, 1997). The junction site of some CTV D RNAs coincided with the transcription start site of subgenomic RNAs (Yang et al., 1997) suggesting their involvement in recombination and in the genome modular evolution of the family Closteroviridae (Bar-Joseph et al., 1997). The stochastic nature of D RNA generation as replication errors was evidenced by the great variety of D RNAs generated after protoplast inoculation with RNA transcripts of the two LIYV genomic segments (Rubio et al., 2002).
The genetic variation and structure of viruses within an infected plant (considered as a virus isolate) also provides important information to understand viral evolution. This has been studied for GRLaV-1 (Alabi et al., 2011), GRLaV-3 (Turturo et al., 2005; Esteves et al., 2012), CYSDV (Rubio et al., 1999, 2001a), ToCV (Lozano et al., 2009), GRLaV-2 (Bertazzon et al., 2010b; Jarugula et al., 2010), and CTV (Ayllón et al., 1999b, 2006; d'Urso et al., 2000; Kong et al., 2000; Rubio et al., 2001b; Davino et al., 2005; Hilf et al., 2005; Melzer et al., 2005; Silva et al., 2007; Gomes et al., 2008; Iglesias et al., 2008; Černi et al., 2009; Oliveros-Garay et al., 2009; Matos et al., 2013; Wu et al., 2013). Analysis of nucleotide sequences or molecular markers such as single-strand conformation polymorphism (Rubio et al., 1996) of a certain number of clones obtained from RT-PCR products showed that most viral isolates had populations composed of a predominant sequence variant and different one- or two-nucleotide mutants in a very low frequency. These mutant clouds are predicted by the quasispecies model as a consequence of the high error frequency in RNA replication and have been described for some animal and plant viruses (Domingo and Holland, 1997). However, some viral isolates had two or more divergent variants, some of which were genetically similar to variants predominant in other viral isolates, suggesting mixed infection by two different strains. These strain mixes that have been found in the closteroviruses CTV (Kong et al., 2000; Rubio et al., 2001b; Sambade et al., 2002; Ayllón et al., 2006; Iglesias et al., 2008) and GLRaV-2 (Bertazzon et al., 2010a), and in the ampeloviruses GLRaV-3 (Farooq et al., 2012) and GLRaV-5 (Esteves et al., 2012), must be common in vector-transmitted viruses infecting long-lived woody hosts. Infection of the same cell with diverged virus strains is a requisite for detectable recombination events to occur.
Co-infection of two viruses or virus strains in mixed infections may have additional evolutionary consequences resulting from their interactions. Sometimes interaction is synergistic, inciting more severe symptoms and increased fitness (virus accumulation) of one or both variants in comparison with single infections. This effect seems to be due to the supression of a host defense mechanism, e.g., gene silencing, by one of the viruses that inhibits accumulation of the other virus in single infections (Palukaitis, 2011). Several cases of synergism have been described between criniviruses and other plant viruses. BYVaV increases concentration of Blackberry virus Y (BVY, genus Brambyvirus, family Potyviridae) in mixed infections (Susaimuthu et al., 2008). CYSDV enhanced multiplication of Cucumber vein yellowing virus (CVYV, genus Ipomovirus, family Potyviridae) and increased symptom severity in mixed infections (Gil-Salas et al., 2012). In sweet potato, SPCSV increased multiplication of several viruses of the genera Potyvirus (Sweet potato feathery mottle virus, Sweet potato latent virus and Sweet potato mild speckling virus), Ipomovirus (Sweet potato mild mottle virus), Cucumovirus (Cucumber mosaic virus), and putative members of the genus Carlavirus (Sweet potato chlorotic fleck virus and C-6 virus) in double and triple infections, which was associated to an increase in the severity of symptoms. In some cases SPCSV titer was reduced indicating an antagonistic interaction (Karyeija et al., 2000; Mukasa et al., 2006; Untiveros et al., 2007). Co-infection of ToCV and Tomato spotted wilt virus (TSWV, genus Tospovirus) in tomato plants susceptible to both viruses resulted in rapid death of the plants, with a pronounced enhancement of ToCV accumulation, whereas TSWV accumulation was not altered. However, in tomato cultivars carrying the Sw-5 gene that confers resistance to TSWV, preinfection with ToCV resulted in TSWV resistance breakage, a phenomenon not observed when plants were simultaneously co-inoculated with both viruses. This suggested that a threshold level or a time lapse is needed for ToCV to interfere or downregulate the defense response in the TSWV-resistant plants (García-Cano et al., 2006). Finally, co-infection of the two tomato-infecting criniviruses TICV and ToCV altered accumulation of each virus in a host-specific manner. While in Nicotiana benthamiana the TICV titer increased and the ToCV titer decreased, in Physalis wrightii the titers of both TICV and ToCV decreased in comparison with the corresponding single infections (Wintermantel et al., 2008). In summary, sinergistic co-infections lead to higher accumulation of at least one of the viruses and may acelerate its adaptation to an initially difficult host.
Other times antagonist interactions may produce a fitness decrease of one or the two viruses. In some cases, previous infection by one viral isolate prevents or hamper subsequent infection by other viral isolate (superinfection exclusion). This phenomenon has been exploited as a disease control strategy named cross protection consisting of preinoculation of the plant with a mild isolate to protect it against damage caused by infection with a virulent isolate. This interaction is most common between genetically related viruses and it has been hypothesized that it might be caused by competition for host resources or because previous infection would trigger the gene silencing antiviral defence of the plant that would impair infection by the second virus (Palukaitis, 2011). CTV cross protection has been efficient for disease control only in some areas (e.g., South Africa and South America) and with some citrus varieties, whereas it has shown limited success in other areas or with other varieties, indicating that (i) cross protection probably depends on the varieties, CTV strains and environmental conditions prevalent in each region (Moreno et al., 2008), and (ii) it is unlikely that CTV cross protection is ruled only by the gene silencing reaction triggered by the pre-inoculated isolate. Indeed preinoculation of citrus plants with artificial hybrid CTV virions containing some genomic segments of the challenging isolate failed to exclude superinfection by this isolate, with only isolates of the same strain being excluded (Folimonova et al., 2010). Demostration that the CTV p33 protein is necessary for superinfection exclusion (Folimonova, 2012) further supports the hypothesis that antiviral silencing reaction triggered by mild isolate pre-inoculation may not be the main mechanism for CTV cross protection. Finally, co-inoculation with a mild and a severe isolate genetically divergent usually resulted in severe symptom expression and predominance of the severe isolate, indicating cross protection failure between divergent isolates and higher fitness of the severe isolate (Roistacher and Dodds, 1993; Sambade et al., 2007; Velazquez-Monreal et al., 2009). Therefore, whatever the cross protection mechanism will be, massive use of cross protecting CTV isolates in some citrus areas is doubtless an important determinant of the viral population structure and an evolutionary factor.
Interactions with the host are one of the most important factors in virus evolution. The tandem of leader proteases of GLRaV-2 seems to have evolved to facilitate infection of this virus in grapevine, a woody perennial host (Liu et al., 2009). Another host effect is the CTV codon usage adaptation to citrus that has also been found in other closteroviruses infecting woody plants but not in those infecting herbaceous hosts (Cheng et al., 2012). Biological and genetic variations of some CTV isolates have been also observed after host change (Ayllón et al., 2006; Scott et al., 2013). Thus, graft-transmission of the mild isolate T317 from citron to sweet orange originated the virulent isolate T318. When T318 was transmitted back to citron, it remained virulent and had the same population structure as it had in sweet orange (Rubio et al., 2000a). This suggests that a minor severe variant contained in the mild isolate T317 was established in sweet orange and became predominant in the T318 population by genetic drift during a transmission bottleneck. This severe variant was as fit in citron as the mild variant (well-adapted to citron and sweet orange). However, when 3 years latter the mild isolate T317 was transmitted again to sweet orange, the new isolate T317D was also mild, albeit the population structure had changed as detected by SSCP analysis. When T317D was transmitted back to citron, its population structure changed back to become indistinguishable from that of the original isolate T317 (Rubio et al., 2000a). This suggests that the severe variant was not sorted in this occasion but some of the mild variants were positively selected as an adaptation to sweet orange. Predominance of these mutants was reverted after back transmission to citron probably as a result of fitness trade-off by host specialization (Woolhouse et al., 2001).
Finally, interactions with the vector can also have an important effect in viral evolution. The association between phylogenetic relationship among members of the three genera of the Closteroviridae family and their type of insect vector (mealybugs for the genus Ampelovirus, whiteflies for Crinivirus and aphids for Closterovirus) likely reflects vector adaptation as a driving force in the diversification in this family (Karasev, 2000). In criniviruses, an emergent group of viruses whose expansion has been linked to the rapid spread of their whitefly vectors (Wisler et al., 1998), the specificity of the association between virus species and whitefly species is a main factor determining the geographical distribution of the different criniviruses. Thus, the displacement of the crinivirus BPYV by CYSDV in Spain has been associated with the increasing populations of Bemisia tabaci in comparison with Trialeurodes vaporariorum (Berdiales et al., 1999). Also disappearance of LIYV in Southern California has been associated with the displacement of the biotype A of B. tabaci by the biotype B, a very poor vector of LIYV (Cohen et al., 1992; McLain et al., 1998; Wisler et al., 1998). Geographical distribution of TICV and ToCV also seems to depend on distribution of the whitefly vectors (Wintermantel, 2008). Interestingly, mixed infections of ToCV and TICV allowed transmission of TICV by the non-vector T. abutilonea (Wintermantel, 2008), indicating an interaction of these viruses during transmission. In ampeloviruses, there is no evidence of vector–virus specificity in the mealybug transmission of different grapevine viruses (Tsai et al., 2010). In closteroviruses, changes in aphid transmissibility of the local CTV isolates have been observed along the years. In the early 1950s, the transmission rate of CTV by Aphis gossypii Glover was very low but this increased in 1960–1970 in Israel and California. (Bar-Joseph and Loebenstein, 1973; Roistacher et al., 1980), suggesting an adaptation and coevolution of CTV to this vector. Brazil and Dominican Republic severe CTV variants are preferentially transmitted by Toxoptera citricida in comparison with the mild components (Brlansky et al., 2011). The absence of this vector in some geographical areas could explain in part why severe CTV isolates are less common in these areas, a situation that may change after introduction of the brown citrus aphid as observed in several countries (Garnsey et al., 2000; Powell et al., 2003). Differences in transmissibility of CTV isolates by different aphid species (Raccah et al., 1976, 1980; Hermoso de Mendoza et al., 1988; Broadbent et al., 1996) indicate the vector exerts a selective pressure and it is an important factor shaping CTV populations. There are evidences that aphid transmission can induce changes in the population structure and/or biological characteristics of individual CTV isolates (Kano and Koizumi, 1991; Ayllón et al., 1999b, 2006; d'Urso et al., 2000; Van Vuuren et al., 2000; Sentandreu et al., 2006; Sambade et al., 2007; Roy and Brlansky, 2009; Velazquez-Monreal et al., 2009; Ananthakrishnan et al., 2010). Comparison of genetic diversity in plants and in aphids showed the occurrence of an important bottleneck for the CTV population during aphid transmission (Nolasco et al., 2008). All these observations indicate that vectors have an important effect in the evolution of closteroviruses by the interplay of natural selection imposed by vector-virus interactions and genetic drift by population bottlenecks during transmission between plants. Also, the movility and dispersibility of vectors determine the level of gene flow between close geographical areas. Thus, the rapid spread of some white-flies which could account in part for the low genetic variation and wide distribution areas of the criniviruses transmitted by them.
Conclusions
The increased number of nucleotide sequences and the availability of more sophisticated analytical tools allows a better understanding of the evolution, population genetics and epidemiology of the viruses in the family Closteroviridae. Analysis of the genetic variability and population structure shows a limited genetic variation as in other plant viruses. This seems to occur mainly by a strong negative selection, indicated by the low number of non-synonymous substutions with respect to the synonymous substitutions in the three genera of this family. Some viruses, had isolate groups with a great genetic divergence between groups but each group being very homogeneous. This suggests that their sequence space is restricted to a few sharp adaptive peaks and that covariation between different nucleotide positions occurs, as suggested by decreased fitness of chimeras of different CTV strains. Long distance movement or gene flow may have contributed in some cases to this low genetic variation.
Although most viral isolates were composed of one major variant and a population of genetically related variants forming a quasispecies structure, some isolates had divergent variants originated from mixed infections of different strains which can affect symptom expression. Interaction between different viruses or variants from the same virus has been found in some cases. The antagonism between close variants of the same virus is the base for cross protection to control damage from severe CTV strains (Moreno et al., 2008). Another element in the virus evolution is the interaction between different virus species, as the synergistic or antagonist interactions observed between the two tomato infecting criniviruses (ToCV and TICV) in different hosts (Wintermantel et al., 2008). Also, viruses of the genus Crinivirus produce increased symptom expression and accumulation of viruses from other families in mixed infections (Karyeija et al., 2000; Mukasa et al., 2006; Untiveros et al., 2007).
Changes in the population structure caused by vector transmission or host change, either by selection of viruses or genetic variants of the same virus or by genetic drift, could explain the appearance of virulent CTV isolates, or the emergence of some viruses in new areas, e.g., criniviruses, associated to their whitefly vector spread.
Recombination between divergent genotypes have been described in the genera Closterovirus and Ampelovirus, which have played an important role in their evolution by increasing genetic diversity and adaptability.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was supported by grants AGL2009-08226 and AGL2012-32429, co-financed by FEDER (European Fund for Regional Development) and by the Ministerio de Economia y Competitividad (MIEC), respectively.
References
Acosta-Leal, R., Duffy, S., Xiong, Z., Hammond, R., and Elena, S. F. (2011). Advances in plant virus evolution: translating evolutionary insights into better disease management. Phytopathology 101, 1136–1148. doi: 10.1094/PHYTO-01-11-0017
Agranovsky, A. (1996). Principles of molecular organization, expression, and evolution of closteroviruses: over the barriers. Adv. Virus Res. 47, 119–158. doi: 10.1016/S0065-3527(08)60735-6
Agranovsky, A. A., Lesemann, D. E., Maiss, E., Hull, R., and Atabekov, J. G. (1995). “Rattlesnake” structure of a filamentous plant RNA virus built of two capsid proteins. Proc. Natl. Acad. Sci. U.S.A. 92, 2470–2473. doi: 10.1073/pnas.92.7.2470
Alabi, O. J., Al Rwahnih, M., Karthikeyan, G., Poojari, S., Fuchs, M., Rowhani, A., et al. (2011). Grapevine leafroll-associated virus 1 occurs as genetically diverse populations. Phytopathology 101, 1446–1456. doi: 10.1094/PHYTO-04-11-0114
Alavi, V., Khatabi, B., and Salekdeh, G. (2005). Comparison of biologically distinct isolates of Citrus tristeza virus from Iran using major coat protein sequences. Australas. Plant Pathol. 34, 577–582. doi: 10.1071/AP05079
Albiach-Martí, M. R., Mawassi, M., Gowda, S., Satyanarayana, T., Hilf, M. E., Shanker, S., et al. (2000). Sequences of Citrus Tristeza Virus separated in time and space sre essentially identical. J. Virol. 74, 6856–6865. doi: 10.1128/JVI.74.15.6856-6865.2000
Ali, A., Li, H., Schneider, W. L., Sherman, D. J., Gray, S., Smith, D., et al. (2006). Analysis of genetic bottlenecks during horizontal transmission of Cucumber mosaic virus. J. Virol. 80, 8345–8350. doi: 10.1128/JVI.00568-06
Ali, A., and Roossinck, M. J. (2008). “Genetic bottlenecks,” in Plant Virus Evolution, ed M. J. Roossinck (Berlin; Heidelberg: Springer), 123–131.
Al Rwahnih, M., Dolja, V. V., Daubert, S., Koonin, E. V., and Rowhani, A. (2012). Genomic and biological analysis of Grapevine leafroll-associated virus 7 reveals a possible new genus within the family Closteroviridae. Virus Res. 163, 302–309. doi: 10.1016/j.virusres.2011.10.018
Alzhanova, D. V., Napuli, A. J., Creamer, R., and Dolja, V. V. (2001). Cell-to-cell movement and assembly of a plant closterovirus: roles for the capsid proteins and Hsp70 homolog. EMBO J. 20, 6997–7007. doi: 10.1093/emboj/20.24.6997
Ananthakrishnan, G., Venkataprasanna, T., Roy, A., and Brlansky, R. (2010). Characterization of the mixture of genotypes of a Citrus tristeza virus isolate by reverse transcription-quantitative real-time PCR. J. Virol. Methods 164, 75–82. doi: 10.1016/j.jviromet.2009.12.001
Angelini, E., Bertazzon, N., and Borgo, M. (2004). Diversity among Grapevine leafroll-associated virus 2 Isolates Detected by Heteroduplex Mobility Assay. J. Phytopathol. 152, 416–422. doi: 10.1111/j.1439-0434.2004.00865.x
Ayllón, M., Navas-Castillo, J., Mawassi, M., Dawson, W., Guerri, J., Flores, R., et al. (1999a). New defective RNAs from Citrus tristeza virus: evidence for a replicase-driven template switching mechanism in their generation. J. Gen. Virol. 80, 817–821.
Ayllón, M. A., Rubio, L., Moya, A., Guerri, J., and Moreno, P. (1999b). The haplotype distribution of two genes of Citrus tristeza virus is altered after host change or aphid transmission. Virology 255, 32–39. doi: 10.1006/viro.1998.9566
Ayllón, M. A., Rubio, L., Sentandreu, V., Moya, A., Guerri, J., and Moreno, P. (2006). Variations in two gene sequences of Citrus tristeza virus after host passage. Virus Genes 32, 119–128. doi: 10.1007/s11262-005-6866-4
Bar-Joseph, M., and Loebenstein, G. (1973). Effects of strain, source plant, and temperature on the transmissibility of Citrus tristeza virus by the melon aphid. Phytopathology 63, 716–720. doi: 10.1094/Phyto-63-716
Bar-Joseph, M., Yang, G., Gafny, R., and Mawassi, M. (1997). Subgenomic RNAs: the possible building blocks for modular recombination of Closteroviridae genomes. Semin. Virol. 8, 113–119. doi: 10.1006/smvy.1997.0116
Bateson, M. F., Lines, R. E., Revill, P., Chaleeprom, W., Ha, C. V., Gibbs, A. J., et al. (2002). On the evolution and molecular epidemiology of the potyvirus Papaya ringspot virus. J. Gen. Virol. 83, 2575–2585.
Berdiales, B., Bernal, J. J., Sáez, E., Woudt, B., Beitia, F., and Rodríguez-Cerezo, E. (1999). Occurrence of Cucurbit yellow stunting disorder virus (CYSDV) and Beet pseudo-yellows virus in cucurbit crops in Spain and transmission of CYSDV by two biotypes of Bemisia tabaci. Eur. J. Plant Pathol. 105, 211–215. doi: 10.1023/A:1008713629768
Bertazzon, N., Borgo, M., and Angelini, E. (2010a). The complete genome sequence of the BD variant of Grapevine leafroll-associated virus 2. Arch. Virol. 155, 1717–1719. doi: 10.1007/s00705-010-0769-y
Bertazzon, N., Borgo, M., Vanin, S., and Angelini, E. (2010b). Genetic variability and pathological properties of Grapevine leafroll-associated virus 2 isolates. Eur. J. Plant Pathol. 127, 185–197.
Bester, R., Maree, H., and Burger, J. (2012). Complete nucleotide sequence of a new strain of Grapevine leafroll-associated virus 3 in South Africa. Arch. Virol. 157, 1815–1819. doi: 10.1007/s00705-012-1333-8
Betancourt, M., Fereres, A., Fraile, A., and García-Arenal, F. (2008). Estimation of the effective number of founders that initiate an infection after aphid transmission of a multipartite plant virus. J. Virol. 82, 12416. doi: 10.1128/JVI.01542-08
Biswas, K., Tarafdar, A., Diwedi, S., and Lee, R. (2012). Distribution, genetic diversity and recombination analysis of Citrus tristeza virus of India. Virus Genes 45, 139–148. doi: 10.1007/s11262-012-0748-3
Boni, M. F., Posada, D., and Feldman, M. W. (2007). An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics 176, 1035–1047. doi: 10.1534/genetics.106.068874
Brlansky, R., Roy, A., and Damsteegt, V. (2011). Stem-pitting Citrus tristeza virus predominantly transmitted by the brown bitrus aphid from mixed infections containing non-stem-pitting and stem-pitting isolates. Plant Dis. 95, 913–920. doi: 10.1094/PDIS-10-10-0772
Broadbent, P., Brlansky, R., and Indsto, J. (1996). Biological characterization of Australian isolates of citrus tristeza virus and separation of subisolates by single aphid transmissions. Plant Dis. 80, 329–333. doi: 10.1094/PD-80-0329
Černi, S., Škorić, D., Ruščiæ, J., Krajaèić, M., Papic, T., Djelouah, K., et al. (2009). East Adriatic—a reservoir region of severe Citrus tristeza virus strains. Eur. J. Plant Pathol. 124, 701–706. doi: 10.1007/s10658-009-9444-0
Černi, S., Ruščić, J., Nolasco, G., Gatin, Ž, Krajačić, M., and Škorić, D. (2008). Stem pitting and seedling yellows symptoms of Citrus tristeza virus infection may be determined by minor sequence variants. Virus Genes 36, 241–249. doi: 10.1007/s11262-007-0183-z
Chare, E. R., and Holmes, E. C. (2004). Selection pressures in the capsid genes of plant RNA viruses reflect mode of transmission. J. Gen. Virol. 85, 3149–3157. doi: 10.1099/vir.0.80134-0
Che, X., Dawson, W. O., and Bar-Joseph, M. (2003). Defective RNAs of Citrus tristeza virus analogous to Crinivirus genomic RNAs. Virology 310, 298–309. doi: 10.1016/S0042-6822(03)00127-2
Che, X., Mawassi, M., and Bar-Joseph, M. (2002). A novel class of large and Infectious defective RNAs of Citrus Tristeza Virus. Virology 298, 133–145. doi: 10.1006/viro.2002.1472
Cheng, X., Wu, X., Wang, H., Sun, Y., Qian, Y., and Luo, L. (2012). High codon adaptation in citrus tristeza virus to its citrus host. Virol. J. 9, 113. doi: 10.1186/1743-422X-9-113
Chooi, K., Pearson, M., Cohen, D., and Pong, J. (2009). “Sequence variation in Grapevine leafroll-associated virus 3 (GLRaV-3) new zealand isolates,” in 16th Meeting of International Council for the Study of Virus and Virus-like Diseases of the Grapevine, Vol. 31, (Dijon: ICVG), 290–291.
Cohen, S., Duffus, J., and Liu, H. (1992). A new Bemisia tabaci biotype in the southwestern United States and its role in silverleaf of squash and transmission of Lettuce infectious yellow virus. Phytopathology 82, 86. doi: 10.1094/Phyto-82-86
Cuellar, W. J., Tairo, F., Kreuze, J. F., and Valkonen, J. P. T. (2008). Analysis of gene content in Sweet potato chlorotic stunt virus RNA1 reveals the presence of the p22 RNA silencing suppressor in only a few isolates: implications for viral evolution and synergism. J. Gen. Virol. 89, 573–582. doi: 10.1099/vir.0.83471-0
Cuevas, J. M., Domingo-Calap, P., and Sanjuán, R. (2012). The fitness effects of synonymous mutations in DNA and RNA viruses. Mol. Biol. Evol. 29, 17–20. doi: 10.1093/molbev/msr179
Davino, S., Rubio, L., and Davino, M. (2005). Molecular analysis suggests that recent Citrus tristeza virus outbreaks in Italy were originated by at least two independent introductions. Eur. J. Plant Pathol. 111, 289–293. doi: 10.1007/s10658-003-2815-z
Dolja, V. V., Karasev, A. V., and Koonin, E. V. (1994). Molecular biology and evolution of closteroviruses: sophisticated build-up of large RNA genomes. Annu. Rev. Phytopathol. 32, 261–285. doi: 10.1146/annurev.py.32.090194.001401
Dolja, V. V., Kreuze, J. F., and Valkonen, J. (2006). Comparative and functional genomics of closteroviruses. Virus Res. 117, 38–51. doi: 10.1016/j.virusres.2006.02.002
Domingo, E., and Holland, J. (1997). RNA virus mutations and fitness for survival. Annu. Rev. Microbiol. 51, 151–178. doi: 10.1146/annurev.micro.51.1.151
d'Urso, F., Ayllón, M. A., Rubio, L., Sambade, A., Hermoso de Mendoza, A., Guerri, J., et al. (2000). Contribution of uneven distribution of genomic RNA variants of citrus tristeza virus (CTV) within the plant to changes in the viral population following aphid transmission. Plant Pathol. 49, 288–294. doi: 10.1046/j.1365-3059.2000.00448.x
Eliasco, E., Livieratos, I., Müller, G., Guzman, M., Salazar, L., and Coutts, R. (2006). Sequences of defective RNAs associated with potato yellow vein virus. Arch. Virol. 151, 201–204. doi: 10.1007/s00705-005-0625-7
Esteves, F., Teixeira Santos, M., Eiras-Dias, J. E., and Fonseca, F. (2012). Occurrence of Grapevine leafroll-associated virus 5 in Portugal: genetic variability and population structure in field-grown grapevines. Arch. Virol. 157, 1747–1765. doi: 10.1007/s00705-012-1371-2
Fajardo, T. V. M., Dianese, É. C., Eiras, M., Cerqueira, D. M., Lopes, D. B., Ferreira, M. A. S. V., et al. (2007). Variability of the coat protein gene of Grapevine leafroll-associated virus 3 in Brazil. Fitopatol. Bras. 32, 335–340. doi: 10.1590/S0100-41582007000400008
Farooq, A. B. U., Ma, Y., Wang, Z., Zhuo, N., Wenxing, X., Wang, G., et al. (2012). Genetic diversity analyses reveal novel recombination events in Grapevine leafroll-associated virus 3 in China. Virus Res. 171, 15–21. doi: 10.1016/j.virusres.2012.10.014
Febres, V., Ashoulin, L., Mawassi, M., Frank, A., Bar-Joseph, M., Manjunath, K., et al. (1996). The p27 protein is present at one end of Citrus tristeza virus particles. Phytopathology 86, 1331–1335.
Fenby, N., Foster, G., Gibson, R., and Seal, S. (2002). Partial sequence of HSP70 homologue gene shows diversity between West African and East African isolates of Sweetpotato chlorotic stunt virus. Trop. Agric. 79, 26–30.
Fisher, L. C., Tennant, P. F., and McLaughlin, W. A. (2010). Detection and characterization of Citrus tristeza virus stem pitting isolates in Jamaica. Eur. J. Plant Pathol. 127, 1–6. doi: 10.1007/s10658-009-9572-6
Folimonova, S. Y. (2012). Superinfection exclusion is an active virus-controlled function that requires a specific viral protein. J. Virol. 86, 5554–5561. doi: 10.1128/JVI.00310-12
Folimonova, S. Y., Robertson, C. J., Shilts, T., Folimonov, A. S., Hilf, M. E., Garnsey, S. M., et al. (2010). Infection with strains of Citrus tristeza virus does not exclude superinfection by other strains of the virus. J. Virol. 84, 1314–1325. doi: 10.1128/JVI.02075-09
Fraile, A., Escriu, F., Aranda, M. A., Malpica, J. M., Gibbs, A. J., and Garcia-Arenal, F. (1997). A century of tobamovirus evolution in an Australian population of Nicotiana glauca. J. Virol. 71, 8316–8320.
Fuchs, M., Martinson, T., Loeb, G., and Hoch, H. (2009). Survey for the three major leafroll disease-associated viruses in Finger Lakes vineyards in New York. Plant Dis. 93, 395–401. doi: 10.1094/PDIS-93-4-0395
García-Arenal, F., and Fraile, A. (2011). “Population dynamics and genetics of plant infection by viruses,” in Recent Advances in Plant Virology, eds C. Caranta, M. A., Aranda, M. Tepfer and J. J. Lopez-Moya (Norfolk, UK: Caister Academic Press), 263–281.
García-Arenal, F., Fraile, A., and Malpica, J. M. (2001). Variability and genetic structure of plant virus populations. Annu. Rev. Phytopathol. 39, 157–186. doi: 10.1146/annurev.phyto.39.1.157
García-Arenal, F., and McDonald, B. A. (2003). An analysis of the durability of resistance to plant viruses. Phytopathology 93, 941–952. doi: 10.1094/PHYTO.2003.93.8.941
García-Cano, E., Resende, R. O., Fernández-Muñoz, R., and Moriones, E. (2006). Synergistic interaction between Tomato chlorosis virus and Tomato spotted wilt virus results in breakdown of resistance in tomato. Phytopathology 96, 1263–1269. doi: 10.1094/PHYTO-96-1263
Garnsey, S., Gottwald, T., Hilf, M., Matos, L., and Borbón, J. (2000). “Emergence and spread of severe strains of Citrus tristeza virus in the dominican republic,” in Proceedings of the 14th Conference of the International Organization of Citrus Virologists, vol. 14, eds J. V. da Graça, R. F. Lee, and R. K. Yokomi (Riverside, CA: IOCV), 57–68.
Gibbs, M. J., Armstrong, J. S., and Gibbs, A. J. (2000). Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics 16, 573–582. doi: 10.1093/bioinformatics/16.7.573
Gil-Salas, F., Peters, J., Boonham, N., Cuadrado, I., and Janssen, D. (2012). Co-infection with Cucumber vein yellowing virus and Cucurbit yellow stunting disorder virus leading to synergism in cucumber. Plant Pathol. 61, 468–478. doi: 10.1111/j.1365-3059.2011.02545.x
Gomes, C. P. C., Nagata, T., de Jesus, W. C. Jr., Neto, C. R. B., Pappas, G. J. Jr., and Martin, D. P. (2008). Genetic variation and recombination of RdRp and HSP70h genes of Citrus tristeza virus isolates from orange trees showing symptoms of citrus sudden death disease. Virol. J. 5, 9. doi: 10.1186/1743-422X-5-9
Gouveia, P., Santos, M. T., Eiras-Dias, J. E., and Nolasco, G. (2011). Five phylogenetic groups identified in the coat protein gene of Grapevine leafroll-associated virus 3 obtained from Portuguese grapevine varieties. Arch. Virol. 156, 413–420. doi: 10.1007/s00705-010-0878-7
Grenfell, B. T., Pybus, O. G., Gog, J. R., Wood, J. L. N., Daly, J. M., Mumford, J. A., et al. (2004). Unifying the epidemiological and evolutionary dynamics of pathogens. Science 303, 327–332. doi: 10.1126/science.1090727
Gultyaev, A. P., Franch, T., and Gerdes, K. (2000). Coupled nucleotide covariations reveal dynamic RNA interaction patterns. RNA 6, 1483–1491. doi: 10.1017/S1355838200990708
Hančević, K., Černi, S., Nolasco, G., Radić, T., Djelouah, K., and Škorić, D. (2013). Biological characterization of Citrus tristeza virus monophyletic isolates with respect to p25 gene. Physiol. Mol. Plant Pathol. 81, 45–53. doi: 10.1016/j.pmpp.2012.10.005
Harper, S., Dawson, T., and Pearson, M. (2009). Complete genome sequences of two distinct and diverse Citrus tristeza virus isolates from New Zealand. Arch. Virol. 154, 1505–1510. doi: 10.1007/s00705-009-0456-z
Harper, S., Dawson, T., and Pearson, M. (2010). Isolates of Citrus tristeza virus that overcome Poncirus trifoliata resistance comprise a novel strain. Arch. Virol. 155, 471–480. doi: 10.1007/s00705-010-0604-5
Hermoso de Mendoza, A., Ballester-Olmos, J., and Pina, J. (1988). “Comparative aphid transmission of a common Citrus tristeza virus isolate and a seedling yellows isolate recently introduced into Spain,” in Proceedings of the 10th Conference of the International Organization of Citrus Virologists, Vol. 10, eds L. W. Timmer, S. M. Garnsey, and L. Navarro (Riverside, CA: IOCV), 68–70.
Herrera-Isidrón, L., Ochoa-Sánchez, J. C., Rivera-Bustamante, R., and Martínez-Soriano, J. P. (2009). Sequence diversity on four ORFs of Citrus tristeza virus correlates with pathogenicity. Virol. J. 6, 1–10. doi: 10.1186/1743-422X-6-116
Hilf, M. (2005). “Partial sequence characterization of Citrus tristeza virus associated with breaking of the general resistance to CTV expressed in Poncirus trifoliata,” in Proceedings of the 16th Conference of the International Organization of Citrus Virologists, Vol. 16, eds M. E. Hilf, N. Duran-Vila, and M. A. Rocha-Peña (Riverside, CA: IOCV), 52–60.
Hilf, M. E., Mavrodieva, V. A., and Garnsey, S. M. (2005). Genetic marker analysis of a global collection of isolates of Citrus tristeza virus: characterization and distribution of CTV genotypes and association with symptoms. Phytopathology 95, 909–917. doi: 10.1094/PHYTO-95-0909
Holland, J., Spindler, K., Horodyski, F., Grabau, E., Nichol, S., and VandePol, S. (1982). Rapid evolution of RNA genomes. Science 215, 1577–1588. doi: 10.1126/science.7041255
Iglesias, N. G., Gago-Zachert, S. P., Robledo, G., Costa, N., Plata, M. I., Vera, O., et al. (2008). Population structure of Citrus tristeza virus from field Argentinean isolates. Virus Genes 36, 199–207. doi: 10.1007/s11262-007-0169-x
Jarugula, S., Alabi, O. J., Martin, R. R., and Naidu, R. A. (2010). Genetic variability of natural populations of Grapevine leafroll-associated virus 2 in Pacific Northwest vineyards. Phytopathology 100, 698–707. doi: 10.1094/PHYTO-100-7-0698
Jiang, B., Hong, N., Wang, G. P., Hu, J., Zhang, J. K., Wang, C. X., et al. (2008). Characterization of Citrus tristeza virus strains from southern China based on analysis of restriction patterns and sequences of their coat protein genes. Virus Genes 37, 185–192. doi: 10.1007/s11262-008-0254-9
Jooste, A. E. C., Maree, H., Bellstedt, D., Goszczynski, D., Pietersen, G., and Burger, J. (2010). Three genetic Grapevine leafroll-associated virus 3 variants identified from South African vineyards show high variability in their 5′ UTR. Arch. Virol. 155, 1997–2006. doi: 10.1007/s00705-010-0793-y
Kano, T., and Koizumi, M. (1991). “Separation of Citrus tristeza virus (CTV) serotypes through aphid transmission” in Proceedings of the 11th Conference of the International Organization of Citrus Virologists, Vol. 11, eds R. H. Brlansky, R. F. Lee, and L. W. Timmer (Riverside, CA: IOCV), 82–85.
Karasev, A. V. (2000). Genetic diversity and evolution of closteroviruses. Annu. Rev. Phytopathol. 38, 293–324. doi: 10.1146/annurev.phyto.38.1.293
Karyeija, R., Kreuze, J., Gibson, R., and Valkonen, J. (2000). Synergistic interactions of a potyvirus and a phloem-limited crinivirus in sweet potato plants. Virology 269, 26–36. doi: 10.1006/viro.1999.0169
Kimura, M. (1980). A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16, 111–120. doi: 10.1007/BF01731581
Kong, P., Rubio, L., Polek, M. L., and Falk, B. W. (2000). Population structure and genetic diversity within California Citrus tristeza virus (CTV) isolates. Virus Genes 21, 139–145. doi: 10.1023/A:1008198311398
Kreuze, J., Savenkov, E., and Valkonen, J. (2002). Complete genome sequence and analyses of the subgenomic RNAs of Sweet potato chlorotic stunt virus reveal several new features for the genus Crinivirus. J. Virol. 76, 9260–9270. doi: 10.1128/JVI.76.18.9260-9270.2002
Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., McGettigan, P. A., McWilliam, H., et al. (2007). Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948. doi: 10.1093/bioinformatics/btm404
Li, H., and Roossinck, M. J. (2004). Genetic bottlenecks reduce population variation in an experimental RNA virus population. J. Virol. 78, 10582–10587. doi: 10.1128/JVI.78.19.10582-10587.2004
Librado, P., and Rozas, J. (2009). DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452. doi: 10.1093/bioinformatics/btp187
Ling, K., Zhu, H., Drong, R., Slightom, J., McFerson, J., and Gonsalves, D. (1998). Nucleotide sequence of the 3′-terminal two-thirds of the Grapevine leafroll-associated virus-3 genome reveals a typical monopartite closterovirus. J. Gen. Virol. 79, 1299–1307.
Liu, Y., Peremyslov, V. V., Medina, V., and Dolja, V. V. (2009). Tandem leader proteases of Grapevine leafroll-associated virus-2: host-specific functions in the infection cycle. Virology 383, 291–299. doi: 10.1016/j.virol.2008.09.035
Livieratos, I., Eliasco, E., Müller, G., Olsthoorn, R., Salazar, L., Pleij, C., et al. (2004). Analysis of the RNA of Potato yellow vein virus: evidence for a tripartite genome and conserved 3′-terminal structures among members of the genus Crinivirus. J. Gen. Virol. 85, 2065–2075. doi: 10.1099/vir.0.79910-0
Lozano, G., Grande-Pérez, A., and Navas-Castillo, J. (2009). Populations of genomic RNAs devoted to the replication or spread of a bipartite plant virus differ in genetic structure. J. Virol. 83, 12973–12983. doi: 10.1128/JVI.00950-09
Marco, C., and Aranda, M. (2005). Genetic diversity of a natural population of Cucurbit yellow stunting disorder virus. J. Gen. Virol. 86, 815–822. doi: 10.1099/vir.0.80584-0
Martin, D., Posada, D., Crandall, K., and Williamson, C. (2005). A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. Aids Res. Hum. Retrovir. 21, 98–102. doi: 10.1089/aid.2005.21.98
Martin, D., and Rybicki, E. (2000). RDP: detection of recombination amongst aligned sequences. Bioinformatics 16, 562–563. doi: 10.1093/bioinformatics/16.6.562
Martin, D. P., Lemey, P., Lott, M., Moulton, V., Posada, D., and Lefeuvre, P. (2010). RDP3: A flexible and fast computer program for analyzing recombination. Bioinformatics 26, 2462–2463. doi: 10.1093/bioinformatics/btq467
Martin, S., Sambade, A., Rubio, L., Vives, M. C., Moya, P., Guerri, J., et al. (2009). Contribution of recombination and selection to molecular evolution of Citrus tristeza virus. J. Gen. Virol. 90, 1527. doi: 10.1099/vir.0.008193-0
Matos, L. A., Hilf, M. E., Cayetano, X. A., Feliz, A. O., Harper, S. J., and Folimonova, S. Y. (2013). Dramatic change in Citrus tristeza virus populations in the Dominican Republic. Plant Dis. 97, 339–345. doi: 10.1094/PDIS-05-12-0421-RE
Martelli, G. P., Agranovsky, A. A., Bar-Joseph, M., Boscia, D., Candresse, T., Coutts, R. H. A., et al. (2011). “Family closteroviridae,” in Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses, eds A. M. Q. King, M. J. Adams, E. B. Carstens, and E. J. Lefkowitz (London, UK: Elsevier Academic Press), 987–1001.
Mawassi, M., Mietkiewska, E., Gofman, R., Yang, G., and Bar-Joseph, M. (1996). Unusual sequence relationships between two isolates of Citrus tristeza virus. J. Gen. Virol. 77, 2359–2364. doi: 10.1099/0022-1317-77-9-2359
Mawassi, M., Mietkiewska, E., Hilf, M. E., Ashoulin, L., Karasev, A. V., Gafny, R., et al. (1995). Multiple species of defective RNAs in plants infected with Citrus tristeza virus. Virology 214, 264–268. doi: 10.1006/viro.1995.9930
Mawassi, M., Satyanarayana, T., Albiach-Martí, M. R., Gowda, S., Ayllón, M. A., Robertson, C., et al. (2000). The fitness of Citrus Tristeza Virus defective RNAs is affected by the lengths of their 5′-and 3′-termini and by the coding capacity. Virology 275, 42–56. doi: 10.1006/viro.2000.0510
McGavin, W., and MacFarlane, S. (2010). Sequence similarities between Raspberry leaf mottle virus, Raspberry leaf spot virus and the closterovirus Raspberry mottle virus. Ann. Appl. Biol. 156, 439–448. doi: 10.1111/j.1744-7348.2010.00401.x
McLain, J., Castle, S., Holmes, G., and Creamer, R. (1998). Physiochemical characterization and field assessment of Lettuce chlorosis virus. Plant Dis. 82, 1248–1252. doi: 10.1094/PDIS.1998.82.11.1248
Melzer, M. J., Borth, W. B., Sether, D. M., Ferreira, S., Gonsalves, D., and Hu, J. S. (2010). Genetic diversity and evidence for recent modular recombination in Hawaiian Citrus tristeza virus. Virus Genes 40, 111–118. doi: 10.1007/s11262-009-0409-3
Melzer, M., Borth, W., Zee, F., Hilf, M., Garnsey, S., and Hu, J. (2005). “Incidence, distribution, and diversity of Citrus tristeza virus in the Hawaiian Islands,” in Proceedings of the 16th Conference of the International Organization of Citrus Virologists, Vol. 16, eds M. E. Hilf, N. Duran-Vila, and M. A. Rocha-Peña (Riverside, CA: IOCV), 179–186.
Moreno, P., Ambrós, S., Albiach-Martí, M. R., Guerri, J., and Peña, L. (2008). Citrus tristeza virus: a pathogen that changed the course of the citrus industry. Mol. Plant Pathol. 9, 251–268. doi: 10.1111/j.1364-3703.2007.00455.x
Moya, A., Holmes, E. C., and Gonzalez-Candelas, F. (2004). The population genetics and evolutionary epidemiology of RNA viruses. Nat. Rev. Microbiol. 2, 279–288. doi: 10.1038/nrmicro863
Mukasa, S., Rubaihayo, P., and Valkonen, J. (2006). Interactions between a crinivirus, an ipomovirus and a potyvirus in coinfected sweetpotato plants. Plant Pathol. 55, 458–467. doi: 10.1111/j.1365-3059.2006.01350.x
Nagy, P. D. (2008). “Recombination in plant RNA viruses,” in Plant Virus Evolution, ed M. J. Roossinck (Berlin; Heidelberg: Springer), 133–156.
Nagy, P. D., and Bujarski, J. J. (1997). Engineering of homologous recombination hotspots with AU-rich sequences in Brome mosaic virus. J. Virol. 71, 3799–3810.
Nolasco, G., Fonseca, F., and Silva, G. (2008). Occurrence of genetic bottlenecks during Citrus tristeza virus acquisition by Toxoptera citricida under field conditions. Arch. Virol. 153, 259–271. doi: 10.1007/s00705-007-1089-8
Nolasco, G., Santos, C., Silva, G., and Fonseca, F. (2009). Development of an asymmetric PCR-ELISA typing method for Citrus tristeza virus based on the coat protein gene. J. Virol. Methods 155, 97–108. doi: 10.1016/j.jviromet.2008.09.030
Oliveros-Garay, O. A., Martinez-Salazar, N., Torres-Ruiz, Y., and Acosta, O. (2009). CPm gene diversity in field isolates of Citrus tristeza virus from Colombia. Arch. Virol. 154, 1933–1937. doi: 10.1007/s00705-009-0530-6
Padidam, M., Sawyer, S., and Fauquet, C. M. (1999). Possible emergence of new geminiviruses by frequent recombination. Virology 265, 218–225. doi: 10.1006/viro.1999.0056
Palukaitis, P. (2011). The Road to RNA Silencing is paved with plant-virus interactions. Plant Pathol. J. 27, 197–206. doi: 10.5423/PPJ.2011.27.3.197
Pamilo, P., and Bianchi, N. O. (1993). Evolution of the Zfx and Zfy genes: rates and interdependence between the genes. Mol. Biol. Evol. 10, 271–281.
Papayiannis, L., Santos, C., Kyriakou, A., Kapari, T., and Nolasco, G. (2007). Molecular characterization of Citrus tristeza virus isolates from Cyprus on the basis of the coat protein gene. J. Plant Pathol. 89, 291–295.
Posada, D., and Crandall, K. A. (2001). Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc. Natl. Acad. Sci. U.S.A. 98, 13757–11376. doi: 10.1073/pnas.241370698
Poudel, B., Sabanadzovic, S., Bujarski, J., and Tzanetakis, I. E. (2012). Population structure of Blackberry yellow vein associated virus, an emerging crinivirus. Virus Res. 169, 272–275. doi: 10.1016/j.virusres.2012.07.009
Powell, C., Pelosi, R., Rundell, P., and Cohen, M. (2003). Breakdown of cross-protection of grapefruit from decline-inducing isolates of Citrus tristeza virus following introduction of the brown citrus aphid. Plant Dis. 87, 1116–1118. doi: 10.1094/PDIS.2003.87.9.1116
Power, A. G. (2000). Insect transmission of plant viruses: a constraint on virus variability. Curr. Opin. Plant Biol. 3, 336–340. doi: 10.1016/S1369-5266(00)00090-X
Raccah, B., Loebenstein, G., and Bar-Joseph, M. (1976). Transmission of Citrus tristeza virus by the melon aphid. Phytopathology 66, 1102–1104. doi: 10.1094/Phyto-66-1102
Raccah, B., Loebenstein, G., and Singer, S. (1980). Aphid-transmissibility variants of citrus tristeza virus in infected citrus trees. Phytopathology 70, 89–93. doi: 10.1094/Phyto-70-89
Roistacher, C., Nauer, E., Kishaba, A., and Calavan, E. (1980). “Transmission of Citrus tristeza virus by Aphis gossypii reflecting changes in virus transmissibility in California,” in Proceedings of the 8th Conference of the International Organization of Citrus Virologists, Vol. 8, eds E. C. Calavan, S. M. Garnsey, and L. W. Timmer (Riverside, CA: IOCV), 76–82.
Roistacher, C. N., and Dodds, J. A. (1993). “Failure of 100 mild Citrus tristeza virus isolates from California to cross protect against a challenge by severe sweet orange stem pitting isolates,” in Proceedings of the 12th Conference of the International Organization of Citrus Virologists, Vol. 12, eds N. Duran-Vila, R. G. Milne, and J. V. da Graça (Riverside, CA: IOCV), 100–107.
Roossinck, M. J. (2003). Plant RNA virus evolution. Curr. Opin. Microbiol. 6, 406–409. doi: 10.1016/S1369-5274(03)00087-0
Roy, A., and Brlansky, R. (2009). Population dynamics of a Florida Citrus tristeza virus isolate and aphid-transmitted subisolates: identification of three genotypic groups and recombinants after aphid transmission. Phytopathology 99, 1297–1306. doi: 10.1094/PHYTO-99-11-1297
Roy, A., and Brlansky, R. (2010). Genome analysis of an orange stem pitting Citrus tristeza virus isolate reveals a novel recombinant genotype. Virus Res. 151, 118–130. doi: 10.1016/j.virusres.2010.03.017
Roy, A., Manjunath, K., and Brlansky, R. (2005). Assessment of sequence diversity in the 5′-terminal region of Citrus tristeza virus from India. Virus Res. 113, 132–142. doi: 10.1016/j.virusres.2005.04.023
Rubio, L., Abou-Jawdah, Y., Lin, H. X., and Falk, B. W. (2001a). Geographically distant isolates of the crinivirus Cucurbit yellow stunting disorder virus show very low genetic diversity in the coat protein gene. J. Gen. Virol. 82, 929–933.
Rubio, L., Ayllon, M. A., Kong, P., Fernandez, A., Polek, M. L., Guerri, J., et al. (2001b). Genetic variation of Citrus tristeza virus isolates from California and Spain: evidence for mixed infections and recombination. J. Virol. 75, 8054–8062.
Rubio, L., Ayllón, M., Guerri, J., Pappu, H., Niblett, C., and Moreno, P. (1996). Differentiation of Citrus tristeza closterovirus (CTV) isolates by single-strand conformation polymorphism analysis of the coat protein gene. Ann. Appl. Biol. 129, 479–489. doi: 10.1111/j.1744-7348.1996.tb05770.x
Rubio, L., Guerri, J., and Moreno, P. (2000a). “Characterization of Citrus tristeza virus isolates by single-strand conformation polymorphism analysis of DNA complementary to their RNA population,” in Proceedings of the 14th Conference of the International Organization of Citrus Virologists, Vol. 14, eds J. V. da Graça, R. F. Lee, and R. K. Yokomi (Riverside, CA: IOCV), 12–17.
Rubio, L., Yeh, H. H., Tian, T., and Falk, B. W. (2000b). A heterogeneous population of defective RNAs is associated with Lettuce infectious yellows virus. Virology 271, 205–212. doi: 10.1006/viro.2000.0318
Rubio, L., Soong, J., Kao, J., and Falk, B. W. (1999). Geographic distribution and molecular variation of isolates of three whitefly-borne closteroviruses of cucurbits: Lettuce infectious yellows virus, Cucurbit yellow stunting disorder virus, and Beet pseudo-yellows virus. Phytopathology 89, 707–711. doi: 10.1094/PHYTO.1999.89.8.707
Rubio, L., Tian, T., Yeh, H. H., Livieratos, Y., and Falk, B. W. (2002). De novo generation of Lettuce infectious yellows virus defective RNAs in protoplasts. Mol. Plant Pathol. 3, 321–327. doi: 10.1046/j.1364-3703.2002.00125.x
Sacristan, S., Malpica, J. M., Fraile, A., and Garcia-Arenal, F. (2003). Estimation of population bottlenecks during systemic movement of Tobacco mosaic virus in tobacco plants. J. Virol. 77, 9906–9911. doi: 10.1128/JVI.77.18.9906-9911.2003
Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425.
Saldarelli, P., Cornuet, P., Vigne, E., Talas, F., Bronnenkant, I., Dridi, A., et al. (2006). Partial characterization of two divergent variants of Grapevine leafroll-associated virus 4. J. Plant Pathol. 88, 203–214.
Salminen, M. O., Carr, J. K., Burke, D. S., and McCutchan, F. E. (1995). Identification of breakpoints in intergenotypic recombinants of HIVtype 1 by Bootscanning. AIDS Res. Hum. Retroviruses 11, 1423–1425.
Sambade, A., Ambrós, S., López, C., Ruíz-Ruíz, S., Hermoso de Mendoza, A., Flores, R., et al. (2007). Preferential accumulation of severe variants of Citrus tristeza virus in plants co-inoculated with mild and severe variants. Arch. Virol. 152, 1115–1126. doi: 10.1007/s00705-006-0932-7
Sambade, A., Lopez, C., Rubio, L., Flores, R., Guerri, J., and Moreno, P. (2003). Polymorphism of a specific region in gene p23 of Citrus tristeza virus allows discrimination between mild and severe isolates. Arch. Virol. 148, 2325–2340. doi: 10.1007/s00705-003-0191-9
Sambade, A., Rubio, L., Garnsey, S. M., Costa, N., Müller, G. W., Peyrou, M., et al. (2002). Comparison of viral RNA populations of pathogenically distinct isolates of Citrus tristeza virus: application to monitoring cross-protection. Plant Pathol. 51, 257–265. doi: 10.1046/j.1365-3059.2002.00720.x
Satyanarayana, T., Gowda, S., Boyko, V., Albiach-Marti, M., Mawassi, M., Navas-Castillo, J., et al. (1999). An engineered closterovirus RNA replicon and analysis of heterologous terminal sequences for replication. Proc. Natl. Acad. Sci. U.S.A. 96, 7433–7438. doi: 10.1073/pnas.96.13.7433
Satyanarayana, T., Gowda, S., Mawassi, M., Albiach-Martì, M. R., Ayllón, M. A., Robertson, C., et al. (2000). Closterovirus encoded HSP70 homolog and p61 in addition to both coat proteins function in efficient virion assembly. Virology 278, 253–265. doi: 10.1006/viro.2000.0638
Schneider, W. L., and Roossinck, M. J. (2001). Genetic diversity in RNA virus quasispecies is controlled by host-virus interactions. J. Virol. 75, 6566–6571. doi: 10.1128/JVI.75.14.6566-6571.2001
Scott, K. A., Hlela, Q., Zablocki, O., Read, D., van Vuuren, S., and Pietersen, G. (2013). Genotype composition of populations of grapefruit-cross-protecting Citrus tristeza virus strain GFMS12 in different host plants and aphid-transmitted sub-isolates. Arch. Virol. 158, 27–37. doi: 10.1007/s00705-012-1450-4
Sentandreu, V., Castro, J., Ayllón, M., Rubio, L., Guerri, J., Gonzalez-Candelas, F., et al. (2006). Evolutionary analysis of genetic variation observed in Citrus tristeza virus (CTV) after host passage. Arch. Virol. 151, 875–894. doi: 10.1007/s00705-005-0683-x
Sharma, A. M., Wang, J., Duffy, S., Zhang, S., Wong, M. K., Rashed, A., et al. (2011). Occurrence of grapevine leafroll-associated virus complex in Napa Valley. PLoS ONE 6:e26227. doi: 10.1371/journal.pone.0026227
Silva, G., Fonseca, F., Santos, C., and Nolasco, G. (2007). Presence of citrus tristeza virus in Angola and São Tomé e Príncipe: characterization of isolates based on coat protein gene analysis. J. Plant Pathol. 89, 149–152.
Silva, G., Marques, N., and Nolasco, G. (2012). The evolutionary rate of Citrus tristeza virus ranks among the rates of the slowest RNA viruses. J. Gen. Virol. 93, 419–429. doi: 10.1099/vir.0.036574-0
Sorrentino, G., Davino, S., Guardo, M., Davino, M., and Caruso, A. (2005). “Incidence and effect of Citrus tristeza virus on commercial orchard of “Tarocco” sweet orange OL in Italy,” in Proceedings of the 16th Conference of the International Organization of Citrus Virologists, Vol. 16, eds M. E. Hilf, N. Duran-Vila, and M. A. Rocha-Peña (Riverside, CA: IOCV), 173–178.
Susaimuthu, J., Tzanetakis, I. E., Gergerich, R. C., Kim, K. S., and Martin, R. R. (2008). Viral interactions lead to decline of blackberry plants. Plant Dis. 92, 1288–1292. doi: 10.1094/PDIS-92-9-1288
Sweiss, M., Anfoka, G., and Abou-Jawdah, Y. (2007). Molecular characterization of jordanian isolates of Cucurbit yellow stunting disorder virus. J. Phytopathol. 155, 557–562. doi: 10.1111/j.1439-0434.2007.01280.x
Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121
Thompson, J. R., Fuchs, M., and Perry, K. L. (2012). Genomic analysis of Grapevine leafroll associated virus- 5 and related viruses. Virus Res. 163, 19–27. doi: 10.1016/j.virusres.2011.08.006
Tian, T., Rubio, L., Yeh, H. H., Crawford, B., and Falk, B. W. (1999). Lettuce infectious yellows virus: in vitro acquisition analysis using partially purified virions and the whitefly Bemisia tabaci. J. Gen. Virol. 80, 1111–1117.
Tsai, C. W., Rowhani, A., Golino, D. A., Daane, K. M., and Almeida, R. P. P. (2010). Mealybug transmission of grapevine leafroll viruses: an analysis of virus-vector specificity. Phytopathology 100, 830–834. doi: 10.1094/PHYTO-100-8-0830
Turturo, C., Saldarelli, P., Yafeng, D., Digiaro, M., Minafra, A., Savino, V., et al. (2005). Genetic variability and population structure of Grapevine leafroll-associated virus 3 isolates. J. Gen. Virol. 86, 217–224. doi: 10.1099/vir.0.80395-0
Tzanetakis, I. E., and Martin, R. R. (2004). Complete nucleotide sequence of a strawberry isolate of Beet pseudoyellows virus. Virus Genes 28, 239–246. doi: 10.1023/B:VIRU.0000025771.48128.f8
Untiveros, M., Fuentes, S., and Salazar, L. F. (2007). Synergistic interaction of Sweet potato chlorotic stunt virus (Crinivirus) with carla-, cucumo-, ipomo-, and potyviruses infecting sweet potato. Plant Dis. 91, 669–676. doi: 10.1094/PDIS-91-6-0669
Van Vuuren, S., Van der Vyver, J., and Luttig, M. (2000). “Diversity among sub-isolates of cross-protecting Citrus tristeza virus isolates in South Africa,” in Proceedings of the 14th Conference of the International Organization of Citrus Virologists, Vol. 14, eds J. V. da Graça, R. F. Lee, and R. K. Yokomi (Riverside, CA: IOCV), 103–110.
Velazquez-Monreal, J., Mathews, D., and Dodds, J. (2009). Segregation of distinct variants from Citrus tristeza virus isolate SY568 using aphid transmission. Phytopathology 99, 1168–1176. doi: 10.1094/PHYTO-99-10-1168
Vives, M. C., Rubio, L., López, C., Navas-Castillo, J., Albiach-Marti, M. R., Dawson, W. O., et al. (1999). The complete genome sequence of the major component of a mild Citrus tristeza virus isolate. J. Gen. Virol. 80, 811.
Vives, M. C., Rubio, L., Sambade, A., Mirkov, T. E., Moreno, P., and Guerri, J. (2005). Evidence of multiple recombination events between two RNA sequence variants within a Citrus tristeza virus isolate. Virology 331, 232–237. doi: 10.1016/j.virol.2004.10.037
Wang, J., Sharma, A. M., Duffy, S., and Almeida, R. P. P. (2011). Genetic Diversity in the 3'Terminal 4.7-kb Region of Grapevine leafroll-associated virus 3. Phytopathology 101, 445–450. doi: 10.1094/PHYTO-07-10-0173
Weir, B. S., and Cockerham, C. C. (1984). Estimating F-statistics for the analysis of population structure. Evolution 38, 1358–1370. doi: 10.2307/2408641
Weng, Z., Barthelson, R., Gowda, S., Hilf, M. E., Dawson, W. O., Galbraith, D. W., et al. (2007). Persistent infection and promiscuous recombination of multiple genotypes of an RNA virus within a single host generate extensive diversity. PLoS ONE 2:e917. doi: 10.1371/journal.pone.0000917
Wintermantel, W. (2008). “Vector specificity of criniviruses in tomato and virus competitiveness during mixed infection,” in Journal of Insect Science 8:4, eds P. A. Stansly and C. L. McKenzie (Duck Key, FL: Fourth International Bemisia Workshop International Whitefly Genomics Workshop), 3–8.
Wintermantel, W. M., Cortez, A. A., Anchieta, A. G., Gulati-Sakhuja, A., and Hladky, L. L. (2008). Co-infection by two criniviruses alters accumulation of each virus in a host-specific manner and influences efficiency of virus transmission. Phytopathology 98, 1340–1345. doi: 10.1094/PHYTO-98-12-1340
Wisler, G., Duffus, J., Liu, H. Y., and Li, R. (1998). Ecology and epidemiology of whitefly-transmitted closteroviruses. Plant Dis. 82, 270–280. doi: 10.1094/PDIS.1998.82.3.270
Woolhouse, M. E. J., Taylor, L. H., and Haydon, D. T. (2001). Population biology of multihost pathogens. Science 292, 1109–1112. doi: 10.1126/science.1059026
Wright, S. (1932). “The roles of mutation, inbreeding, crossbreeding and selection in evolution,” in Proceedings of the Sixth International Congress of genetics, Vols. 1, 6, 356–366.
Wu, G., Pan, S., Wang, G., Tang, M., Liu, Y., Yang, F., et al. (2013). The genotypes of Citrus tristeza virus isolates from China revealed by sequence analysis of multiple molecular markers. Arch. Virol. 158, 231–235. doi: 10.1007/s00705-012-1475-8
Keywords: Closterovirus, Crinivirus, Ampelovirus, recombination, selection, phylogeny, gene flow
Citation: Rubio L, Guerri J and Moreno P (2013) Genetic variability and evolutionary dynamics of viruses of the family Closteroviridae. Front. Microbiol. 4:151. doi: 10.3389/fmicb.2013.00151
Received: 15 February 2013; Accepted: 29 May 2013;
Published online: 26 June 2013.
Edited by:
Ricardo Flores, Instituto de Biología Molecular y Celular de Plantas, SpainReviewed by:
Scott J. Harper, University of Florida, USAIoannis E. Tzanetakis, University of Arkansas, USA
Copyright © 2013 Rubio, Guerri and Moreno. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Luis Rubio, Instituto Valenciano de Investigaciones Agrarias, Ctra. Moncada-Naquera Km 5, 46113 Moncada, Valencia, Spain e-mail: lrubio@ivia.es