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 ) . 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 . 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 protein - coding gene is an accelerated nonsynonymous ( dn , amino acid replacing ) rate relative to the synonymous ( ds , silent ) rate of substitutions , with the rate ratio  = dn / ds > 1 . as silent 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 8 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 2742 kda in size ; the virion measures 2835 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 5 terminus is capped only in carnation mottle carmovirus , red clover necrotic mosaic dianthovirus and maize chlorotic mottle machlomovirus . 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 . 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 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  = 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 , 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 . the only parameter considered is  , representing the cost of mutation relative to recombination . when  changes from 0 to 1 , the cost of mutation weighted by  increases , and the cost for recombination weighted by 1 -  decreases . 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 . 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 , 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 ) . 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 ) ] . 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 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 . 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 . 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 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 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 . 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 . 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 . 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 . 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 , 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 ) . 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 ) ] . 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 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 . 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 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 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 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 . 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 . 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.31 ) , 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 ) the most frequently recombining viruses were : ocsv , rvx.rna 1 , jinrsv , and pmv with 175 , 166 , 134 , and 108 putative recombination sites , respectively . 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 , bbsv-val25.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 . 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 ) . these rna silencing suppressors have been identified in the following genera : aureusvirus , carmovirus , tombusvirus , and dianthovirus ( voinnet e al . , 1999 ; 2003 ; mrai et al . , 2005 ; takeda et al . , 2005 ) . 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 . 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 .