Contribution of recombination and selection to molecular evolution of Citrus tristeza virus

The genetic variation of Citrus tristeza virus (CTV) was analysed by comparing the predominant sequence variants in seven genomic regions (p33, p65, p61, p18, p13, p20 and p23) of 18 pathogenically distinct isolates from seven different countries. Analyses of the selective constraints acting on each codon suggest that most regions were under purifying selection. Phylogenetic analysis shows diverse patterns of molecular evolution for different genomic regions. A first clade composed of isolates that are genetically close to the reference mild isolates T385 or T30 was inferred from all genomic regions. A second clade, mostly comprising virulent isolates, was defined from regions p33, p65, p13 and p23. For regions p65, p61, p18, p13 and p23, a third clade that mostly included South American isolates could not be related to any reference genotype. Phylogenetic relationships among isolates did not reflect their geographical origin, suggesting significant gene flow between geographically distant areas. Incongruent phylogenetic trees for different genomic regions suggested recombination events, an extreme that was supported by several recombination-detecting methods. A phylogenetic network incorporating the effect of recombination showed an explosive radiation pattern for the evolution of some isolates and also grouped isolates by virulence. Taken together, the above results suggest that negative selection, gene flow, sequence recombination and virulence may be important factors driving CTV evolution.