SciELO - Scientific Electronic Library Online

vol.34 número4Expressed sequence tag analysis of khat (Catha edulis) provides a putative molecular biochemical basis for the biosynthesis of phenylpropylamino alkaloidsOptimizing expression and purification of an ATP-binding gene gsiA from Escherichia coli k-12 by using GFP fusion índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados



Links relacionados


Genetics and Molecular Biology

versão impressa ISSN 1415-4757

Genet. Mol. Biol. vol.34 no.4 São Paulo  2011 Epub 16-Set-2011 

Positive selection, molecular recombination structure and phylogenetic reconstruction of members of the family Tombusviridae: implication in virus taxonomy



Moncef Boulila

Institut de l'Olivier, Sousse Ibn-khaldoun, Tunisia

Send Correspondence to




A detailed study of putative recombination events and their evolution frequency in the whole genome of the currently known members of the family Tombusviridae, comprising 79 accessions retrieved from the international databases, was carried out by using the RECCO and RDP version 3.31β algorithms. The first program allowed the detection of potential recombination sites in seven out of eight virus genera (Aureusvirus, Avenavirus, Carmovirus, Dianthovirus, Necrovirus, Panicovirus, and Tombusvirus), the second program provided the same results except for genus Dianthovirus. On the other hand, both methods failed to detect recombination breakpoints in the genome of members of genus Machlomovirus. Furthermore, based on Fisher's Exact Test of Neutrality, positive selection exerted on protein-coding genes was detected in 17 accession pairs involving 15 different lineages. Except genera Machlomovirus, and Panicovirus along with unclassified Tombusviridae, all the other taxonomical genera and the unassigned Tombusviridae encompassed representatives under positive selection. The evolutionary history of all members of the Tombusviridae family showed that they segregated into eight distinct groups corresponding to the eight genera which constitute this family. The inferred phylogeny reshuffled the classification currently adopted by the International Committee on Taxonomy of Viruses. A reclassification was proposed.

Key words: bioinformatics, phylogeny, recombination, positive selection, Tombusviridae, sequence, taxonomy.




RNA recombination is one of the major factors responsible for the generation of new RNA viruses and retroviruses. The biological mechanisms of recombination differ across organisms, but in broad terms recombination results in the creation of mosaic sequences where the evolutionary history at each site may be different. Recombination, defined as the exchange of genetic information between two nucleotide sequences, is an important process that influences biological evolution at many different levels. Recombination explains a considerable amount of genetic diversity in natural populations and, in general, genes located in regions of the genome with low levels of recombination have low levels of polymorphism (Posada and Crandall, 2001). Recombination reshuffles existing variation and even creates new variants. It has been shown that RNA recombination enables the exchange of genetic material, not only between the same or similar viruses but also between distinctly different viruses (Worobey and Holmes, 1999). Sometimes, it also permits crossovers between viral and host RNA (Greene and Allison, 1994; Aaziz and Tepfer, 1999; Baroth et al., 2000; Nagai et al., 2003). Taking into account the structure of viral genomic molecules and the location of crossover sites, three basic types of RNA recombination were distinguished: homologous, aberrant homologous and non-homologous (Lai, 1992; Alejska et al., 2001). The former two occur between two identical or similar RNAs (or between molecules displaying local homology), while the latter involves two different molecules. Most of the collected data suggest that RNA recombinants are formed according to a copy choice model (Alejska et al., 2001). A viral replication complex starts nascent RNA strand synthesis on one template, called RNA donor, and then switches to another template, called RNA acceptor. Accordingly, two main factors are thought to affect RNA recombination: the structure of recombining molecules and the ability of the viral replicase to switch templates. Through generations, viral populations evolve under various selective forces at different regions and sites that display different functional constraints. A stringent and robust criterion for detecting adaptive evolution in a proteincoding gene is an accelerated nonsynonymous (dN, amino acid replacing) rate relative to the synonymous (dS, silent) rate of substitutions, with the rate ratio 0 = dN/dS >1.Assilent mutations do not change the amino acid whereas replacement mutations do, the difference in their fixation rates provides a measure of selective pressure on the protein.

Amongst positive-strand plant RNA viruses, the family Tombusviridae encompasses several viruses with an important economical impact. According to the 8th ICTV (International Committee on Taxonomy of Viruses) report (Fauquet et al., 2005), the family Tombusviridae includes the following genera: Tombusvirus, Carmovirus, Necrovirus, Dianthovirus, Machlomovirus, Avenavirus, Aureusvirus and Panicovirus. According to the Baltimore classification, the viruses in this family are classified as Type IV viruses, and are part of the luteovirus supergroup (Habili and Symons, 1989). The RNA is contained in an icosahedral (T = 3) capsid, composed of 180 units of a single coat protein 27-42 kDa in size; the virion measures 28-35 nm in diameter, and is not enveloped. All Tombusviridae have a positive-sense, single-stranded linear genome, with the exception of dianthoviruses, whose genome is bipartite. The genome is approximately 4-5.4 kb in length, depending on the genus. The 3' terminus is not polyadenylated. The 5' terminus is capped only in Carnation mottle carmovirus, Red clover necrotic mosaic dianthovirus and Maize chlorotic mottle machlomovirus. The genome encodes 4-6 ORFs. The polymerase ORF encodes an amber stop codon that is the site of a readthrough event within ORF 1 (except in dianthoviruses, where readthrough occurs via a frameshift), producing two products necessary for replication. There is no helicase encoded by the virus. The replication process of members of family Tombusviridae comprises the following steps: (i) the virus penetrates into the host cell, (ii) the viral genomic RNA is uncoated and released into the cytoplasm, (iii) the viral RNA is translated to produce the two proteins necessary for RNA synthesis (replication and transcription), (iv) a negative-sense complementary ssRNA is synthesized using the genome RNA as a template, (v) a new genomic RNA is synthesized using the negative-sense RNA as a template, (vi) the RNA-dependant RNA polymerase (RdRp) recognizes internal subgenomic promoters on the negative-sense RNA, to transcribe the 3' co-terminal subgenomic RNAs that will generate the capsid and movement protein, (vii) new virus particles are formed (White and Nagy, 2004).

The main objective of this work was to determine and characterize virus evolution mechanisms of the Tombusviridae based on the occurrence of putative recombination events and positive selection in their full-length genome. This was achieved by the analysis of 79 accessions obtained from GenBank. As a result, we propose a reclassification according to their predicted evolutionary history.


Materials and Methods

The sequences of the entire genome of 79 accessions cataloged in GenBank were used in this study (Table 1).

The nucleotide sequences were aligned using programs CLUSTALW 2.0.9 and CLUSTALX 2.0.9 (Larkin et al., 2007) with default configuration. Their phylogenetic relationships were determined with the Maximum-likelihood (ML) algorithm incorporated in the MEGA version 5 program (Tamura et al., 2011) under assumption of the substitution models proposed by Jukes and Cantor (1969) (JC), Hasegawa et al. (1985) (HKY85), and Tamura and Nei (1993) (TN93). Bootstrap analyses with 500 replicates were performed to assess the robustness of the branches.

Using the MEGA4.1β program (Kumar et al., 2008), positive selection was inferred by the counting method described by Nei and Gojobori (1986) and, later on, by Suzuki and Gojobori (1999). According to this method, the phylogenetic tree of sequences analyzed was used. For the parsimony method, the total numbers of synonymous (cS) and nonsynonymous (cN) substitutions as well as the average numbers of synonymous (sS) and nonsynonymous (sN) sites per codon over the phylogenetic tree for each codon site were computed according to the maximum parsimony principle (Fitch, 1971; Hartigan, 1973). The null hypothesis of selective neutrality (rS= rN or 0 = 1) was tested for each site by computing the probability (p) of obtaining the observed or more biased values for cS and cN, which were assumed to follow a binomial distribution with the probabilities of occurrence of synonymous and nonsynonymous substitutions given by sS/(sS+ sN) and sN/(sS+ sN), respectively. Positive selection is inferred when p < 0.05 and cN/sN> cS/sS (Suzuki, 2006).

Potential recombination events between diverged nucleotide sequences were explored using two programs: RDP v3.31β (Martin et al., 2005b) and RECCO (Maydt and Lengauer, 2006). RDP incorporates several published recombination detection methods into a single suite of tools: RDP (Martin and Rybicki, 2000), GENECONV (Padidam et al., 1999), BOOTSCAN (Martin et al., 2005a), MAXCHI (Smith, 1992), CHIMAERA (Posada and Crandall, 2001), SISCAN (Gibbs et al., 2000), and 3SEQ (Boni et al., 2007). In all cases, default parameters were used. Only events predicted by more than half of the methods are considered as significant. The algorithm developed and described by Maydt and Lengauer (2006) as being a fast, simple and sensitive method for detecting recombination in a set of sequences and locating putative recombination breakpoints is based on cost minimization. This method has only two tunable parameters, recombination and mutation cost. In practice the only parameter considered is α, representing the cost of mutation relative to recombination. When a changes from 0 to 1, the cost of mutation weighted by a increases, and the cost for recombination weighted by 1-α decreases. In other words, parameter α controls the ambiguity between mutation and recombination.



Recombination events during Tombusviridae evolution

Examination of the RECCO program output regarding the occurrence of recombination events in the complete genome of the Tombusviridae family, revealed that three out of five aureusviruses were putative recombinants (PoLV.Pigeonpea, JCSMV.Iran, MaWLMV.USA). In contrast, CLSV (unknown isolate) and CLSV.Canada did not show any recombinant signal (Table 2). Within the genus Aureusvirus, the most frequently recombining virus was PoLV.Pigeonpea (33 putative recombination sites), whereas only 28 possible recombination signals were detected in the genome of viruses JCSMV.Iran and MaWLMV.USA. Similarly, the only representative of the genus Avenavirus (OCSV) was a potential recombinant with 175 putative sites. The RDP package confirmed these results for both genera. Among the carmoviruses, 14 out of 30 members were possible recombinants. According to RECCO, the most frequently recombining virus was JINRSV with 134 putative events, while MeNSV.Nagasaki and MeNSV.NK had the lowest number of recombination sites (two putative sites). The RDP v3.31β algorithm confirmed the occurrence of possible recombination events only for accessions MeNSV.NK and MeNSV.NH. Recombination investigations of genus Dianthovirus based on RECCO analysis showed that only 80% of them were possible recombinants (CarRSV-RNA 1 and 2, RCNMV-RNA 1 and 2, RCNMV.Can. RNA 1 and 2, SCNMV.59.RNA 1, and RVX.RNA 1). in contrast, SCNMV.59.RNA 2 and SCNMV.38.RNA did not show any putative recombination signals. However, the RDP package did not predict recombination in the dianthoviruses (Table 3). Although the most frequently recombining necrovirus was RVX (166 putative sites), RCNMV.RNA 2 had only two putative sites. Based on RECCO analysis, 50% of the necroviruses (BBSV.Val25.Iran, LWSV, TNV-A.FMB, TNV-A.C, TNV-D.Hungarian, and OMMV-GP.Portugal,) were possible recombinants. Conversely, BBSV, BBSV.CO.USA, BBSV.Xinjiang, OLV-1.Citrus, OLV-1.GM6.Portugal, and TNV-D.Rhotamsted were not possible recombinants. These results were congruent with those obtained with the RDP package. While the most frequently recombining virus was LWSV (39 sites), TNV.A.C recombined into two sites. Regarding the sole representative of genus Panicovirus (PMV), the results obtained by the two methods (RECCO and RDP v3.31β) were incongruent. Indeed, with RECCO, 108 possible sites were detected, whereas no recombination signals were found with the RDP package. A similar situation was observed with regard to the newly proposed carmoviruses (NLVCV.Alaska, PCRPV.GR 57, PLPV.PV-0193). According to RECCO analysis, although NLVCV.Alaska was the most frequently recombining virus (65 sites), PLPV.PV-0193 recombined only into 38 sites (Table 4). Regarding the members of genus Tombusvirus, there was an agreement between the two methods indicating that 80% of the analyzed accessions were putative recombinants. While CBLV had the highest number of putative recombination signals (67 sites), TBSV.Cherry had only two recombination sites. Furthermore, it is noteworthy that the two representatives of genus Machlomovirus (MCMoV, and MCMoV.Nebraska) were not recombinants as assessed by the two methods of analysis used in this study. Seeking for the recombination frequency in the genome of the Tombusviridae, two-thirds of the aureusviruses (JCSMV.Iran, and MaWLMV.USA) showed that in most cases, their breakpoint length was a single residue. In contrast, the breakpoint length of most putative recombination sites of PoLV.Pigeonpea was between three and 37 nucleotides (Table 2). Also, the breakpoint length of the major recombination sites of the single representative of genus Avenavirus (OCSV) consisted of a single residue. In about 50% of the members of the genus Carmovirus, the length of their most detected recombination sites was a single residue. As opposed to that, the breakpoint interval of the remaining members exceeded three residues reaching a size as long as 82 residues (MeNSV.NH). In 62% of the investigated dianthoviruses, the breakpoint length exceeded three nucleotides reaching 100 residues (CarRSV.RNA 1) (Table 3). In the necroviruses, the breakpoint interval distribution was similar i.e., 50% of the breakpoints consisted of a single residue, while the remaining breakpoints were between three and 77 nucleotides. For the sole member of the genus Panicovirus (PMV), most of the recombination sites had a breakpoint length of a single residue (45) (Table 3). As for the tombusviruses, 75% showed a breakpoint length exceeding three residues up to 161 nucleotides (AMoCV.Bari) (Table 4).

Nucleotide sequence analysis

Maximum composite likelihood estimate of the nucleotide substitution pattern were made using the MEGA4.1β program. The results for Tombusviridae showed that the rates of different transitional substitutions varied from 3.18 to 14.61, and those of transversional substitutions varied from 6.6 to 8.57. The nucleotide frequencies were: 0.269 (A), 0.258 (T/U), 0.207 (C), and 0.266 (G). The transition/transversion rate ratios were k1 = 1.705 (purines) and k2 = 0.482 (pyrimidines). The overall transition/transversion bias was R = 0.547, where R = [AGk1+ TCk2]/[(A+G)(T+C)]. There were a total of 1218 positions in the final dataset. In all these analyses, the codon positions included were first + second + third + noncoding. All positions containing gaps and missing data were excluded from the dataset (complete deletion option).

The MEGA4.1β program also incorporates the Tajima's Neutrality Test. The purpose of this test is to indentify sequences which do not fit the neutral theory model at equilibrium between mutation and genetic drift. Tajima's test compares a standardized measure of the total number of segregating sites (the polymorphic DNA sites) in the sampled DNA and the average number of mutations between pairs in the sample. Tajima's D was determined (D = 5.280926).

Positive selection

The high genetic stability of viruses can be attributed to negative or purifying selection to maintain the functional integrity of the viral genome. The degree of negative selection in genes, or the degree of functional constraint for the maintenance of the encoded protein sequence, can be estimated, as mentioned above, by the ratio between the nucleotide diversities in nonsynonymous and synonymous positions (dN/dS). For most coding genes the dN/dS ratio is < 1 which is consistent with negative selection against protein change. In contrast, a dN/dS ratio > 1 may be an indication that adaptive or positive selection is driving gene divergence. In this study, pairwise comparisons of all screened accessions showed that, none of the members of the genera Machlomovirus and Panicovirus, and unclassified Tombusviridae was under positive selection. On the contrary, the genera Aureusvirus (JCMSV.Iran), Avenavirus (OCSV), Carmovirus (CarMoV.China, CarMoV.Indian), Dianthovirus (CarRSV-RNA 2, RCNMV-RNA 2, SCNMV-RNA 2.59, SCNMV-RNA 2.38), Necrovirus (BBSV, BBSV.Val25.Iran), Tombuvirus (GALV.nipplefruit, PeLV, TBSV.Statice, PNSV) along with the unassigned Tombusviridae (PLPV.PV.0193) were under positive selection (Table 5). It is worth pointing out that, in the viruses with a segmented genome, positive selection was detected only in RNA 2, suggesting that probably reassortment events occurred. All these results were obtained by testing neutrality in sequence pairs with Fisher's Exact Test. The probability of rejecting the null hypothesis of strict-neutrality (dN= dS) in favor of positive selection for each sequence pair was determined. Values of p less than 0.05 were considered significant at the 5% level. The variance of the difference (dN-dS) was computed using the bootstrap method (500 replicates). All analyses were made using the Nei-Gojobori method incorporated in the MEGA program. All positions containing gaps and missing data were excluded from the dataset (complete deletion option). The final dataset comprised a total of 234 positions.

Phylogenetic relationships

The phylogenetic relationships among members of the family Tombusviridae, based on the sequences of their complete genome, were inferred using a Maximum Likelihood algorithm under the assumption of three models of substitution (JC, HKY85, TN93). The topologies of the constructed trees were identical. The inferred phylogeny showed that each taxonomical genus in the family Tombusviridae constituted a homogenous group clearly distinct from the others. However, the results obtained in this study evidenced a few differences in terms of virus species composition within each taxonomical genus compared to the current classification adopted by the ICTV. In fact, three viruses considered by the ICTV as unassigned (PLPV.PV-0193) and unclassified Tombusviridae (NLVCV.Alaska, PCRPV.GR 57) showed a close phylogenetic relationship to known members of the genus Carmovirus. Moreover, the viruses belonging to this genus were divided into two distinct subgroups. The first subgroup comprised viruses: TCV, CCFV, JINRV, HCRSV, PLPV, PCRPV, NLVCV, SCV, AFBV, PFBV, CPMoV, SYMoMV and CarMoV, and the second subgroup encompassed viruses: MeNSV, and PSNV. Furthermore, it was proposed that genus Necrovirus should be constituted by two distinct subgroups named tentative Subgroup I (BBSV, TNV.D, LWSV) and tentative subgroup II (OMMV, TNV.A, OLV-1) (Figure 1). It should be noted that here OMMV is an integral part of subgroup I rather than an unclassified Necrovirus. In contrast, genus Aureusvirus encompassed members that evolved in a homogenous manner: CLSV, PoLV, MaWLMV, and JCSMV. Similarly, the following members of genus Tombusvirus also formed a coherent ensemble: MaNSV, CBLV, LNV.L, LNV.Zantedeschia, PeLV, CNV, CymRSV, AMoCV, TBSV.Statice, TBSV.Nipplefruit, TBSV.Pepper, TBSV.Cherry, GALV, PNSV, and CarIRSV. Their evolutionary history reshuffled the existing classification adopted by the ICTV since 2009. In fact, according to this classification, MaNSV was considered as an unassigned Tombusviridae, whereas LNV and PNSV were included in the unclassified Tombusvirus group. Concerning genus Dianthovirus which clearly was not monophyletic, the clustering pattern showed two distinct clades representing their RNAs 1 and 2, as illustrated in Figure 1. Originally, RVX was considered as an unclassified virus within genus Dianthovirus.



This study evidenced the prediction of putative recombination events in the genome of several members of the family Tombusviridae and demonstrated that tombusviruses and carmoviruses are highly recombinant compared to viruses of the other genera. For this purpose, two methods were chosen (RECCO and RDP v.3.311), based on the fact that they are appropriate for the mosaic structure of viruses as reported in previous works (Boulila, 2009; 2010). In this study, using the RECCO algorithm, it was demonstrated that the viruses belonging to the following genera contained putative recombination signals in their genome: Aureusvirus, Avenavirus, Carmovirus, Dianthovirus, Necrovirus, Panicovirus, and Tombusvirus. These results were in good agreement with those obtained by the RDP package except for members of genus Diantovirus. By both methods, the two representatives of genus Machlomovirus (MCMoV, MCMoV.Nebraska) were found to be nonrecombinant. As revealed by RECCO, the most frequently recombining viruses were: OCSV, RVX.RNA 1, JINRSV, and PMV with 175, 166, 134, and 108 putative recombination sites, respectively. All of these recombination signals were constituted by a single residue. MeNSV.Nagasaki, MeNSV.NK, RCNMV.RNA 2, TNV-A.C, and TBSV.Cherry (2 sites), MeNSV.MNSV-Al, CarRSV-RNA 2, and TNV-A.FM1B (3 sites), MeNSV.Kochi, BBSVVal25.Iran, and TNV-D.Hungarian (4 sites), CarMoV, MeNSV.MNSV-264, and MeNSV.NH (5 sites), RCNMV.RNA 1.Can, TBSV.nipplefruit, and TBSV.pepper (7 sites), and RCNMV.RNA 1 (8 sites) showed the lowest frequency of recombination breakpoints. In contrast, most of these breakpoints had an interval exceeding three nucleotides. Furthermore, this study showed that recombination may occur between viruses belonging to different genera. For example: Oat chlorotic stunt avenavirus (OCSV) and Melon necrotic spot carmovirus (MeNSV) may give rise to Pothos latent aureusvirus (PoLV). Similarly, OCSV itself may result from a recombination between Turnip crinkle carmovirus (TCV) and Maize white line mosaic aureusvirus (MaWLMV) (Table 2). Seemingly, these viruses could contain part of their sequences particularly in the coat protein-encoding gene of each other. Such an event was largely studied for Cucumber necrosis tombusvirus (CNV) and Melon necrotic spot carmovirus (MeNSV) (Riviere and Rochon, 1990).

On the other hand, investigations of selective pressure acting on protein expression of virus genes led to the identification of positive selection in 17 accession pairs involving 15 different lineages. It is worth mentioning that numerous viruses: JCSMV.Iran, OCSV, CarRSV-RNA 2, RCNMV-RNA 2.Can, BBSV.Val25.Iran, GALV.nipplefruit, TBSV.Statice, PNSV, and PLPV.PV.0193 evolved under both mechanisms: recombination and positive selection between which synergism might be occurring. Such a synergism between recombination and natural selection may have played a major role in Darwinian molecular evolution.

The evolutionary history of the Tombusviridae has shown that the 79 accessions split into eight clearly separated clusters representing the eight genera of the Tombusviridae family. From the present phylogenetic study, at least two taxonomic implications can be drawn: (i) three viruses (NLVCV.Alaska, PCRPV.GR 57, PLPV.PV-0193) currently considered by the ICTV as: one unassigned Tombusviridae (PLPV.PV-0193), and two unclassified Tombusviridae (NLVCV.Alaska and PCRPV.GR 57). All of them should be included in genus Carmovirus; (ii) In addition to the viruses belonging to genus Carmovirus which have formed two separated subgroups, the members of genera Necrovirus, and Dianthovirus evolved separately and divided into two distinct subgroups as shown in Figure 1. In contrast, members of genera Aureusvirus, and Tombusvirus formed separately a single ensemble. The evolutionary relationships among viruses are a reliable approach for classification. As stated by Stuart et al. (2004) (who reported similar results regarding the genetic divergence of components of genus Necrovirus), the comparison of complete genomes is a more balanced approach that should provide a more precise scheme of relatedness. On the other hand, it should be pointed out that, in genus Dianthovirus, the genetic divergence between RNAs 1 and 2 is correlated to the final products synthesized and their use by the virus to survive. For example: RNA silencing is a small RNA-guided sequence-specific gene activation mechanism in eukaryotes that is involved in different biological phenomena (e.g. development, heterochromatin formation and defense against molecular parasites such as viruses). Many viruses express suppressors to counteract RNA-silencingmediated antiviral defenses. These RNA silencing suppressors have been identified in the following genera: Aureusvirus, Carmovirus, Tombusvirus, and Dianthovirus (Voinnet eal., 1999; Qu et al., 2003; Mérai et al., 2005; Takeda et al., 2005). Dianthovirus uses a unique strategy to suppress RNA silencing. The dianthoviral suppressor consists of multiple components including P27, P88 (encoded by two ORFs in RNA 1) and viral RNA (Takeda et al., 2005). Moreover, sequence variability of the coat protein-coding gene (RNA 1) may be linked to the interaction between this structural protein and the host and vector which themselves show a major diversity among diathoviruses. In contrast, the ORF in RNA 2 encodes the movement protein. All these factors can influence the divergence between the two RNAs.

Finally, to the author's best knowledge, this is the largest study in the literature so far on recombination potentially occurring in the entire genome of all currently known members of the family Tombusviridae as well as positive selection operating on protein expression and their phylogenetic reconstruction. In addition, a reclassification based on their predicted evolutionary history, is proposed.



Aaziz R and Tepfer M (1999) Recombination in RNA viruses and in virus-resistant transgenic plants. J Gen Virol 80:1339-1346.         [ Links ]

Alejska M, Kurzyniska-Kokorniak A, Broda M, Kierzek R and Figlerowicz M (2001) How RNA viruses exchange their genetic material. Acta Bioch Polon 48:391-407.         [ Links ]

Baroth M, Orlich M, Thiel HJ and Becher P (2000) Insertion of cellular NEDD8 coding sequences in a pestivirus. Virology 278:456-466.         [ Links ]

Boni MF, Posada D and Feldman MW (2007) An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics 176:1035-1047.         [ Links ]

Boulila M (2009) Recombination structure and genetic relatedness among members of the family Bromoviridae based on their RNAs 1 and 2 sequence analyses. Virus Genes 38:435-444.         [ Links ]

Boulila M (2010) Putative recombination events and evolutionary history of five economically important viruses of fruit trees based on the coat protein-encoding gene sequence analysis. Biochem Genet 48:357-375.         [ Links ]

Fauquet CM, Mayo MA, Maniloff J, Desselberger U and Ball LA (2005) Virus taxonomy: Classification and nomenclature of viruses. In: Eighth Report of the International Committee on Taxonomy of Viruses. Elsevier/Academic Press, London, 1259 pp.         [ Links ]

Fitch WN (1971) Toward defining the course of evolution minimum change for a specific tree topology. Syst Zool 20:406-416.         [ Links ]

Gibbs MJ, Armstrong JS and Gibbs AJ (2000) Sister-scanning: A Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics 16:573-582.         [ Links ]

Greene AE and Allison RF (1994) Recombination between viral RNA and transgenic plants transcripts. Science 263:1423-1425.         [ Links ]

Habili N and Symons RH (1989) Evolutionary relationship between luteoviruses and other RNA plant viruses based on sequence motifs in their putative RNA polymerase and nucleic acid helicases. Nucleic Acids Res 17:9543-9555.         [ Links ]

Hartigan JA (1973) Minimum mutation fits to a given tree. Biometrics 29:53-63.         [ Links ]

Hasegawa M, Kishino H and Yano T (1985) Dating of human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol 22:160-174.         [ Links ]

Jukes T and Cantor C (1969) Evolution of protein molecules. In: Munro HN (ed) Mammalian Protein Metabolism. Academic Press, New York, pp 21-132.         [ Links ]

Kumar S, Nei M, Dudley J and Tamura K (2008) MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinfor 9:299-306.         [ Links ]

Lai MMC (1992) RNA recombination in animal and plant viruses. Microbiol Rev 56:61-79.         [ Links ]

Larkin MA, Blackshileds G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R et al. (2007) Clustal W and Clustal X v. 2.0. Bioinformatics 23:2947-2948.         [ Links ]

Martin D and Rybicki E (2000) RDP: Detection of recombination amongst aligned sequences. Bioinformatics 16:562-563.         [ Links ]

Martin DP, Posada D, Crandall KA and Williamson C (2005a) A modified bootscan algorithm for automated identification of recombination sequences and recombination breakpoints. AIDS Res H Retrovir 21:98-102.         [ Links ]

Martin DP, Williamson C and Posada D (2005b) RDP2: Recombination detection and analysis from sequence alignments. Bioinformatics 21:260-262.         [ Links ]

Maydt J and Lengauer T (2006) Recco: Recombination analysis using cost optimization. Bioinformatics 22:1064-1071.         [ Links ]

Mérai Z, Kerényi Z, Molnar A, Barta E, Valoczi A, Bisztray G, Havelda Z, Burgyan J and Silhavy D (2005) Aureuvirus P14 is an efficient RNA silencing suppressor that binds doublestranded RNAs without size specificity. J Virol 79:7217-7226.         [ Links ]

Nagai M, Sakoda Y, Mori M, Hayashi M, Kida H and Akashi H (2003) Insertion of a cellular sequence and RNA recombination in the structural protein coding region of cytopathogenic bovine viral diarrhea virus. J Gen Virol 84Pt2:447-452.         [ Links ]

Nei M and Gojobori T (1986) Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 3:418-426.         [ Links ]

Padidam M, Sawyer S and Fauquet CM (1999) Possible emergence of new geminiviruses by frequent recombination. Virology 265:218-225.         [ Links ]

Posada D and Crandall K (2001) Evaluation of methods for detecting recombination from DNA sequences. Computer simulation. Proc Natl Acad Sci USA 98:13757-13762.         [ Links ]

Qu F, Ren T and Morris TJ (2003) The coat protein of Turnip crinkle virus suppresses post transcriptional gene silencing at an early initiation step. J Virol 77:511-522.         [ Links ]

Riviere CJ and Rochon DM (1990) Nucleotide sequence and genomic organization of melon necrotic spot virus. J Gen Virol 71:1887-1896.         [ Links ]

Smith JM (1992) Analyzing the mosaic structure of genes. J Mol Evol 34:126-129.         [ Links ]

Stuart G, Moffet K and Bozarth RF (2004) A whole genome perspective on the phylogeny of the plant virus family Tombusviridae. Arch Virol 149:1595-1610.         [ Links ]

Suzuki Y (2006) Statistical properties of the methods for detecting positively selected amino acid sites. Gene 365:125-129.         [ Links ]

Suzuki Y and Gojobori T (1999) A method for detecting positive selection at single amino acid sites. Mol Biol Evol 16:1315-1328.         [ Links ]

Tajima F (1989) Statistical-method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585-595.         [ Links ]

Takeda A, Tsukuda M, Mizumoto H, Okamoto K, Kaido M, Mise K and Okuno T (2005) A plant RNA virus suppresses RNA silencing through viral RNA replication EMBO J 24:3147-3157.         [ Links ]

Tamura K and Nei M (1993) Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10:512-526.         [ Links ]

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. Molecular Biology and Evolution. Doi 10.1093/molbev/msr121.         [ Links ]

Voinnet O, Pinto YM and Baulcombe DC (1999) Suppression of gene silencing: A general strategy used by diverse DNA and RNA viruses of plants. Proc Natl Acad Sci USA 96:14147-14152.         [ Links ]

White KA and Nagy PD (2004) Advances in the molecular biology of tombusviruses: Gene expression, genome replication and recombination. Prog Nucleic Acid Res Mol Biol 78:187-226.         [ Links ]

Worobey M and Holmes EC (1999) Evolutionary aspects of recombination in RNA viruses. J Gen Virol 80:2535-2543.         [ Links ]



Send Correspondence to:
Moncef Boulila
Institut de l'Olivier, B.P. 14
4061 Sousse Ibn-Khaldoun, Tunisia

Received: November 22, 2010; Accepted: May 13, 2011.



Associate Editor: Carlos F.M. Menck
License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.