Assessment of codivergence of Mastreviruses with their plant hosts
© Wu et al; licensee BioMed Central Ltd. 2008
Received: 15 May 2008
Accepted: 18 December 2008
Published: 18 December 2008
Viruses that have spent most of their evolutionary time associated with a single host lineage should have sequences that reflect codivergence of virus and host. Several examples for RNA viruses of host-virus tree congruence are being challenged. DNA viruses, such as mastreviruses, are more likely than RNA viruses to have maintained a record of host lineage association.
The full genomes of 28 isolates of Wheat dwarf virus (WDV), a member of the Mastrevirus genus, from different regions of China were sequenced. The analysis of these 28 entire genomes and 18 entire genome sequences of cereal mastreviruses from other countries support the designation of wheat, barley and oat mastrevirus isolates as separate species. They revealed that relative divergence times for the viruses WDV, Barley dwarf virus (BDV), Oat dwarf virus (ODV) and Maize streak virus (MSV) are proportional to divergence times of their hosts, suggesting codivergence. Considerable diversity among Chinese isolates was found and was concentrated in hot spots in the Rep A, SIR, LIR, and intron regions in WDV genomes. Two probable recombination events were detected in Chinese WDV isolates. Analysis including further Mastrevirus genomes concentrated on coding regions to avoid difficulties due to recombination and hyperdiversity. The analysis demonstrated congruence of trees in two branches of the genus, but not in the third. Assuming codivergence, an evolutionary rate of 10-8 substitutions per site per year was calculated. The low rate implies stronger constraints against change than are obtained by other methods of estimating the rate.
We report tests of the hypothesis that mastreviruses have codiverged with their monocotyledonous hosts over 50 million years of evolution. The tests support the hypothesis for WDV, BDV and ODV, but not for MSV and other African streak viruses.
Viruses are a class of genetic elements dependent on suitable host cells for their propagation. Viruses belonging to diverse viral groups have been proposed to have codiverged with their hosts based on congruence of phylogenetic trees for the viruses with those for their hosts. In codivergence, congruence results from long association of the viral and host lineages. The term codivergence is preferred to describe this situation since, unlike the term "coevolution", it does not imply that the association necessarily provides mutual benefits to the partners .
The best studied examples of codivergence  include the hantaviruses  and arenaviruses [4, 5] in their murid hosts, and potyviruses  and tobamoviruses [7, 8] in plant hosts. Three of these examples have recently been challenged. A reanalysis of hantaviruses, including data on shrew hantaviruses, has called the codivergence of hantaviruses with hosts into question . The role of recombination in generating arenavirus phylogenetic trees is under dispute . A recent reanalysis of Potyvirus divergence suggested that this genus emerged shortly after the beginning of agriculture , much later than was earlier proposed. The recent analysis of the evolution of rymoviruses  supports the view that evolution of sobemoviruses is more rapid than that of their hosts. In contrast, the addition of further tobamoviral sequences to the Tobamovirus tree has supported the congruence of host and virus trees . Further, there is evidence that, in very long-term analyses, such as through studies of viruses in herbarium specimens  and Greenland ice cores , viruses in the genus are evolving extremely slowly, such that codivergence is a possibility, while having nucleotide substitution frequencies of the order of 10-5 substitutions per site in the shorter term .
The viruses in these major examples of putative codivergence of viruses and hosts have RNA genomes that replicate using error prone RNA-dependent RNA polymerases encoded by the viral genomes. Such viruses are expected to evolve more rapidly than viruses with DNA genomes which use host DNA-dependent DNA polymerases with proof-reading ability for replication of their genomes and could be subject to the action of DNA repair systems on replication errors or spontaneous mutations . Thus, for DNA-containing viruses, mutation frequencies similar to those of host genomes are expected, making the observation of codivergence more likely for these viruses than for viruses with RNA genomes.
Members of the Geminiviridae replicate their DNA using a host DNA polymerase and encapsidate circular single-stranded DNAs . This plant virus family is one of the largest, represented by four genera: Mastrevirus, Curtovirus, Topocuvirus and Begomovirus, classified depending on their vectors, host range and genomic characteristics [18–20]. During the last two decades these viruses have emerged as devastating pathogens, threatening crop production and causing huge economic losses . Today, geminivirus-induced diseases are among the most economically important in vegetable and field crops, including beans, cassava, cotton, maize, pepper, tomato and wheat [20–24].
In the process of studying populations from China of Wheat dwarf virus (WDV) from the Mastrevirus genus of Geminiviridae, we observed patterns that suggested that viruses in this genus have substitution frequencies consistent with their replication by host DNA polymerases. The genus Mastrevirus consists of viruses with circular single-stranded (ss) DNA genomes in geminate (twinned) virions , and has 11 recognized species including WDV. WDV is transmitted in a persistent circulative manner by the leafhopper Psammotetix striatus L. to barley, wheat, oats, rye and many wild grasses [25, 26]. It was first described by Vacke  in the western parts of the former Czechoslovak Socialist Republic (CSR) and then found in many parts of the world [27, 28]. Its distribution areas are increasing and it has recently been detected in Germany , Tunisia , Turkey , Finland , Zambia  and China . The complete genome sequences of 18 isolates, 10, 7 and 1 from wheat, barley and oats respectively, have been determined from the CSR, Sweden, Hungary, France, Germany, Turkey and China . Comparisons of these sequences showed that the isolates which infected wheat, barley and oats respectively, formed three distinct clades [25, 28]. Schubert et al.  suggested reclassifying WDV into three species according to sequence differences and host range studies: WDV, Barley dwarf virus (BDV), and Oat dwarf virus (ODV), designations used in this paper.
Phylogenetic analysis of viruses
In the spring of 2004, 2005 and 2006, several diseased wheat plants showing extreme dwarfing, various types of yellowing, and reduced or no heading were found during field surveys in many wheat fields of China . Wheat samples collected from northern, central, northwestern and southwestern areas tested positive by PCR for WDV, suggesting that WDV was widely distributed throughout China. The full genomes of 28 isolates from different regions of China were sequenced in this study. Details of these, together with those of the 18 complete WDV, BDV and ODV genomes already published, are provided in supplemental material [see Additional file 1 and 2]. Phylogenetic trees were constructed by neighbor-joining (NJ) (Figure 1) and maximum-parsimony as described in Methods using Maize streak virus as an outgroup. The topologies of the two types of trees were identical at all branch points that were well supported by bootstrap analysis (> 70%) but differed at some branch points with low statistical support (Figure 1).
Unscaled evolutionary distances from the common ancestors to their sequence progeny were deduced from the NJ tree. The distances indicated that the barley isolates diverged from wheat isolates much more recently than the oat isolate did. The most basal barley isolate was TR2 from Turkey, followed by those from Germany and the CSR. Wheat isolates were separated into two groups, the first one from Europe and the second one from France, the CSR, and China. The grouping of all of the Chinese isolates in one clade with the CSR and French isolates suggested that they evolved from the same ancestor. The Chinese wheat isolates were divided into two clades, but were not clustered by geographical source or collection time (Figure 1). For example, 10 isolates from Shijiazhuang, Hebei province, were classified into four different sub-clades in the phylogenetic trees (Figure 1).
Correlation of host lineage and WDV divergence estimates.
Mastrevirus sequence alignment
To explore further the possibility of the codivergence of virus and host, we included additional members of the Mastrevirus family in the phylogenetic analysis. To align these additional sequences with the previously analyzed sequences in an unambiguous and reliable fashion, we considered elimination of regions of high diversity and regions prone to be different due to recombination. Regions with high nucleotide diversity were difficult to align reliably. To identify regions that could be reliably aligned, the level of genetic diversity (θw and π) [38, 39] of WDV was examined in genomic regions (Table 2) or along the entire genomes of Chinese isolates (Figure 2). Examination of nucleotide diversity values (Table 2) revealed few differences among the coding regions. For all WDV isolates, θw was significantly higher for the Rep region than for the other three regions, but this difference was not significant for the π nucleotide diversity. The higher Rep diversity was most apparent among non-Chinese isolates, for which both measures showed significance. Overall, as expected from Figure 1, diversity values were lower for Chinese isolates than for non-Chinese isolates. For the Chinese isolates, the π diversity values were slightly higher for Rep and RepA regions than for CP and MP regions. Theta values were not significantly different from one another. As is generally expected for mastreviruses, the most diverse regions over the 37 genomes were LIR and SIR [25, 28]. Considering only the Chinese isolates, however, the diversity of the SIR region was not distinguishable from those of the coding regions. In the non-Chinese isolates, MP diversity could not be distinguished from diversity for the non-coding regions. Among the non-Chinese isolates, the SIR region was more diverse than LIR, while the opposite was true for the Chinese isolates. Overall, SIR was more diverse by the θw measure, but not significantly so using the π diversity value.
WDV nucleotide diversities.
Nucleotide substitution for coding regions of the WDV genome.
Phylogenetic comparison of viruses and hosts
The ICTV proposed a reduced similarity value for the demarcation of species in the case of mastreviruses (75%), compared to that for the majority of geminiviruses (89%) . In our analysis, we found 68% nucleotide sequence identity between wheat and barley isolates, and 59.7%–69.3% between oat and barley isolates. The nucleotide sequence identities between oat and wheat isolates were 68.3%–69.4%, similar to the results of others . The phylogenetic trees also supported strongly the proposal that WDV should be divided into three mastrevirus species , a proposal with which we agree.
Several observations favor the hypothesis that mastreviruses codiverged with their hosts. First, despite the known propensity for mastreviruses to recombine during evolution , exemplified by the evidence of recombination in the evolution of two WDV isolates from China (Figure 3), the similarity of Rep-RepA and MP-CP Mastrevirus trees to one another (Figure 4) suggests that recombination had been eliminated as a major factor shaping trees by concentrating on the coding regions. Second, the topologies of maize, oat, barley and wheat lineages and the lineages of their mastreviruses (MSV, ODV, BDV, and WDV) were identical (Figure 1) as were the topologies of dicotyledonous plants and their viruses (Figure 4). Third, not only did the topologies agree, but there was an excellent correlation between estimated divergence times of the plant lineages and the relative divergence times of the mastreviruses, as calculated from the sequences (Table 1).
Nevertheless, support for the hypothesis of codivergence of host and virus is not conclusive. That the molecular clock has not been ticking uniformly in all lineages and in both halves of the viral genome removes the molecular clock as a tool for testing the codivergence hypothesis. Evidence for uneven ticking of the clock, such as in MisSV and ChStrMV branches, may be attributable to recombination. Yet, it is interesting that there was statistical support for a clock operating in one part of the overall tree. That part specifically included two of the branch points for which there are dates associated with divergences of plant lineages (of Triticum from Hordeum and of Triticeae from Aveneae). It is also noteworthy that adding the third datable branch (PACCAD clade from BEP clade) led to a near perfect correlation of host divergence times with relative viral divergence times (Table 1).
Non-congruence of PACCAD clade host branching patterns with branching patterns of viruses named for those hosts seems to argue against the codivergence hypothesis. However, it must be remembered in this regard that virus names are based on the host plant from which the virus was first isolated. Thus, the names do not necessarily reflect the plant lineage in which the virus spent most of its time evolving. Indeed, non-congruence is expected for viruses that can efficiently infect many species of plants. MSV isolates have been found in many grass genera including Zea, Panicum, Setaria, Urochloa and even Triticum . Further, their presence in maize is of recent origin since the streak viruses are indigenous to Africa and Indian Ocean islands and maize was only introduced to Africa after the European discovery of America . Similarly, sugar cane is not native to Africa but is infected by a complex of related virus species (Sugarcane streak virus, Sugarcane streak Reunion virus, Sugarcane streak Egypt virus and Eragrostis streak virus) which have been isolated from native grasses of the genera Setaria, Cenchrus, Paspalum and Eragrostis, indicating a wide host range . In contrast, the BEP clade viruses, WDV, BDV and ODV specialize in infecting their respective host plants and thus likely have evolved entirely in the lineage for which they are named .
Estimates of short term and longer term evolution rates are available for viruses in the sister genus, Begomovirus. Inoculation of plants with infectious cloned DNA of Tomato yellow leaf curl China virus (TYLCCV) resulted in the subsequent recovery of viral sequences with substitutions at a frequency of about 10-4 substitutions per site during a 60 day growth period in plants . Consistent with the experimental result, phylogenetic analysis suggested a substitution frequency of 5 × 10-4 per site per year for the related Tomato yellow leaf curl virus . Thus on these time scales, the genomes change as frequently as most RNA virus genomes [48, 49]. These time scales sample sites that evolve rapidly. The distribution of nucleotide diversity along the WDV genome (Figure 2) shows that areas of the genome with high diversity represent only a small percentage of the total genome. In the TYLCCV study , four nucleotide positions accounted for close to half (18 of 41) of the observed substitutions. Thus, it is likely that highly mutable positions gave rise to the substitutions in the Begomovirus investigations, while the deeper phylogenetic tree construction employed in our Mastrevirus work focuses on the areas of the genome with low nucleotide diversity. These are also the areas subjected to strong purifying selection. The apparent discrepancy between very long-term and short or long-term evolution rates in the Geminiviridae is reminiscent of similar findings in the Tobamovirus analysis. Understanding this apparent difference awaits further analysis.
WDV was collected throughout China during field surveys in the growing seasons 2004 to 2006. The 28 isolates described here originated from wheat planted in different agro-ecological areas in China, including the northwestern (Shaanxi and Gansu provinces), northern (Shanxi and Hebei provinces), central (Henan province), and southwestern (Yunan province) areas. All the field isolates were inoculated to the susceptible wheat (Triticum aestivum L.) cultivar Fengkang No. 8 by vector leafhoppers (Psammotetix alienus L.) to increase virus concentration and to allow serological typing or sequencing of polymerase chain reaction (PCR) products. The wheat plants were later tested for WDV with ELISA using an antiserum (Bio-Rad, Marnes la Coquette, France). Leaves were collected from WDV-positive plants displaying typical symptoms of WDV infection, and stored at -80°C. Details of the isolates, their names, provinces of collection, original host plant, and years of collection are shown in additional file [see Additional file 1].
Cloning of entire genomes
Total DNA was extracted from systemically WDV-infected wheat leaves . DNA extracts were used as template for PCR amplification, performed in a 50 μL reaction solution containing 10×Taq Buffer, 2.5 mM dNTP (each), 0.4 mM of the viral sense and complementary sense primers designed according to the conserved sequences of WDV genomes [see Additional file 3], and Ampli Taq DNA polymerase (Applied Biosystems, Foster City, CA, USA). PCR reactions were carried out for 35 cycles, each consisting of denaturation at 94°C for 1 min, annealing at 55°C for 1 min, and extension at 72°C for 1 min, with 95°C for 2 min at the beginning and 72°C for 10 min at the final step. The expected PCR products were 767 bp, 1152 bp and 1041 bp, using the primer pairs of 40F/806R, 735F/1886R and 1828F/118R, respectively, and together covered the entire length of the viral genome. The PCR product segments were electrophoresed in 1.0% agarose gels, bands were excised using a razor blade and purified using the BioTeq PCR quick Gel Extraction Kit (BioTeq, Inc, USA).
Nucleotide sequences of the entire genome of each isolate were determined using the above PCR fragments. The purified fragments were cloned into the pMD18-T vector (Takara, Dalian, China). The plasmids were transformed into Escherichia coli strain JM110 and plasmid DNA was isolated from overnight cultures by alkaline lysis. Insert sequences were determined on at least three clones for each PCR fragment using the dideoxynucleotide chain termination method by an automated sequencer (ABI BigDye 3.1, Applied Biosystems, Foster City, CA). Sequence data were assembled using DNASIS version 3.5 (Hitachi) or BIOEDIT version 5.0.9 . The nucleotide sequence data have been submitted to GenBank databases and assigned accession numbers EF536859 through EF536886.
Phylogenetic and molecular diversity analysis
Complete genomes of the 28 WDV isolates sequenced in this study and 18 entire sequences of other WDV, BDV and ODV isolates obtained from the NCBI database (National Center for Biotechnology Information, Bethesda, MD, USA) were analyzed. The coding and intergenic regions were annotated by reading frame or following NCBI's annotations. The complete genome sequences of the WDV, BDV and ODV genomes were aligned with CLUSTAL W V.1.8. MEGA V.4.0  determined the number of nucleotide substitutions per site (evolutionary distance) between the strains. Phylogenetic trees were constructed by neighbor-joining (NJ), and maximum parsimony (MP) as implemented by MEGA version 4.0  and DNAPARS of PHYLIP package version 3.5 , respectively, based on the Kimura 2-parameter distance matrix model. Bootstrap confidence values were obtained for 1000 replicates (Figure 1). The homologous regions of the genome of an isolate of Maize streak virus (MSV) (NC_001346)  were used as the outgroup for these analyses, as BLAST searches had shown them to be the sequences in the international sequence databases most closely related to those of MSV. Treemap 4.1.1  was used to test and display the correspondence between plant and virus trees. The Watterson's estimator of θ (θw)  and the average pairwise nucleotide diversity Pi(π) , were estimated using DnaSP version 4.10.2 . Also, the program was used to estimate the proportions of synonymous and nonsynonymous substitutions by the Jukes-Cantor one-parameter model.
To evaluate the sequence relationships among mastrevirus genomes, a selection of mastrevirus sequences available at the time was obtained [see Additional file 2]. A manually adjusted multiple sequence alignment based on encoded amino acid sequences was generated using Se-Al such that MP-CP and Rep-RepA regions were satisfactorily aligned. These regions were separately excised from the alignment for further manipulation. In cases where multiple sequences were available for the same named virus, a consensus sequence was generated by Se-Al. For examination of host phylogenetic relationships, rbcL sequences were retrieved from GenBank/DDBJ/EMBL. They are identified in the legend of Figure 4. These were also aligned. Aligned sequences were examined using PAUP  by testing models supplied by Modeltest . The parameters for the best model were used to construct neighbor joining trees as implemented in PAUP. To test for consistency, bootstrapped neighbor joining was also performed using Phylip package programs, Seqboot, DNAdist, Fitch, and Consense . Resulting trees were manually manipulated to minimize tangles between host and virus trees. The validity of a molecular clock for several assemblages of sequences was tested using the log ratio test as described by Posada . To determine selection models acting on the WDV, BDV and ODV genes, nucleotide substitutions at synonymous (Ks) and non-synonymous (Ka) positions of genes were calculated by DnaSP version 4.10.2 .
Nucleotide Substitution Frequency and Divergence Times
The average frequency with which mastrevirus sequence sites mutate in evolution is unknown. For highly mutable sites the substitution frequency has been estimated  at 3 × 10-4/site during a 60 day growing period for Tomato yellow leaf curl China virus (TYLCCNV). That value is clearly an overestimate of the average frequency.
Nevertheless, since our interest was in determining the relative ratios of divergence times of WDV, BDV and ODV from MSV, of ODV from WDV and BDV, and of WDV from BDV, the number was used to obtain divergence time estimates with MEGA software . Host divergence times were obtained from literature. Resulting calculated virus divergence times were normalized to 100 for the (WDV-BDV-ODV)-MSV split and values were plotted against corresponding host divergences. Linear regression was used to evaluate the correspondences and to determine an appropriate conversion factor that could be applied, assuming uniformity of the molecular clock, to the relative divergence times of the viruses.
Detection of recombination and mutation bias
To investigate the extent of recombination within the data set, the aligned sequences were examined using the Recombination Detection Program (RDP3) , GENECONV , BOOTSCAN , MAXIMUM CHISQUARE , CHIMERA , SISTER SCAN  and and phylpro  recombination detection methods as implemented in RDP3 , (details of program settings available from http://darwin.uvigo.es/rdp/heath2006.zip). The transversion and transition differences of all pairs of sequences were calculated using the discalc program (kindly supplied by G. F. Weiller, Australian National University) and these were compared in diplomo scatter plots .
This work was supported by the National Key Basic Research Program of China (973 program No.2006CB101903), Special Funding of State Key Laboratory, China (SKL2007SR05) and the Oklahoma Agricultural Experiment Station. We thank Dr. Andrew Doust for supplying the aligned nucleotide sequences of rbcL genes.
- Charleston MA, Perkins SL: Traversing the tangle: algorithms and applications for cophylogenetic studies. J Biomed Inform. 2006, 39: 62-71. 10.1016/j.jbi.2005.08.006.View ArticlePubMedGoogle Scholar
- Jackson AP, Charleston MA: A cophylogenetic perspective of RNA-virus evolution. Mol Biol Evol. 2004, 21: 45-57. 10.1093/molbev/msg232.View ArticlePubMedGoogle Scholar
- Plyusnin A, Morzunov SP: Virus evolution and genetic diversity of hantaviruses and their rodent hosts. Curr Top Microbiol Immunol. 2001, 256: 47-75.PubMedGoogle Scholar
- Mills JN, Bowen MD, Nichol ST: African arenaviruses – Coevolution between virus and murid host?. Belg J Zool. 1997, 127 (Suppl 1): 19-28.Google Scholar
- Bowen MD, Peters CJ, Nichol ST: Phylogenetic analysis of the Arenaviridae: patterns of virus evolution and evidence for cospeciation between Arenaviruses and their rodent hosts. Molec Phylogen Evol. 1997, 8: 301-316. 10.1006/mpev.1997.0436.View ArticleGoogle Scholar
- Ward CW, Weiller GF, Shukla DD, Gibbs A: Molecular systematics of the Potyviridae, the largest plant virus family. Molecular Basis of Virus Evolution. Edited by: Gibbs A, Calisher CH, García-Arenal F. 1995, New York: Cambridge University Press, 477-500.View ArticleGoogle Scholar
- Gibbs A: Evolution and origins of tobamoviruses. Philos Trans R Soc Lond Ser B-Biol Sci. 1999, 354: 593-602. 10.1098/rstb.1999.0411.View ArticleGoogle Scholar
- Lartey RT, Voss TC, Melcher U: Tobamovirus evolution: Gene overlaps, recombination, and taxonomic implications. Mol Biol Evol. 1996, 13: 1327-1338.View ArticlePubMedGoogle Scholar
- Ramsden C, Holmes EC, Charleston MA: Hantavirus evolution in relation to its rodent and insectivore hosts: no evidence for co-divergence. Mol Biol Evol. 2009, 26 (1): 143-153. 10.1093/molbev/msn234.View ArticlePubMedGoogle Scholar
- Cajimat MNB, Milazzo ML, Hess BD, Rood MP, Fulhorst CF: Principal host relationships and evolutionary history of the North American arenaviruses. Virology. 2007, 367: 235-243. 10.1016/j.virol.2007.05.031.PubMed CentralView ArticlePubMedGoogle Scholar
- Gibbs AJ, Ohshima K, Phillips MJ, Gibbs MJ: The prehistory of potyviruses: their initial radiation was during the dawn of agriculture. PLoS ONE. 2008, 3: e2523-10.1371/journal.pone.0002523.PubMed CentralView ArticlePubMedGoogle Scholar
- Fargette D, Pinel-Galzi A, Sérémé D, Lacombe S, Hébrard E, Traoré O, Konaté G: Diversification of rice yellow mottle virus and related viruses spans the history of agriculture from the neolithic to the present. PLoS Pathog. 2008, 4: e1000125-10.1371/journal.ppat.1000125.PubMed CentralView ArticlePubMedGoogle Scholar
- Srinivasan KG, Min BE, Ryu KH, Adkins S: Determination of complete nucleotide sequence of Hibiscus latent Singapore virus: Evidence for the presence of an internal poly(A) tract. Arch Virol. 2005, 150: 153-166. 10.1007/s00705-004-0404-x.View ArticlePubMedGoogle Scholar
- Fraile A, Escriu F, Aranda MA, Malpica JM, Gibbs AJ, García-Arenal F: A century of tobamovirus evolution in an Australian population of Nicotiana glauca. J Virol. 1997, 71: 8316-8320.PubMed CentralPubMedGoogle Scholar
- Castello JD, Rogers SO, Starmer WT, Catranis CM, Ma LJ, Bachand GD, Zhao YH, Smith JE: Detection of tomato mosaic tobamovirus RNA in ancient glacial ice. Polar Biol. 1999, 22: 207-212. 10.1007/s003000050411.View ArticleGoogle Scholar
- Malpica JM, Fraile A, Moreno I, Obies CI, Drake JW, García-Arenal F: The rate and character of spontaneous mutation in an RNA virus. Genetics. 2002, 162: 1505-1511.PubMed CentralPubMedGoogle Scholar
- Duffy S, Holmes EC: Phylogenetic evidence for rapid rates of molecular evolution in the single-stranded DNA Begomovirus Tomato Yellow Leaf Curl Virus. J Virol. 2008, 82: 957-965. 10.1128/JVI.01929-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Fauquet CM, Bisaro DM, Briddon RW, Brown JK, Harrison BD, Rybicki EP, Stenger DC, Stanley J: Revision of taxonomic criteria for species demarcation in the family Geminiviridae, and an updated list of begomovirus species. Arch Virol. 2003, 148: 405-421. 10.1007/s00705-002-0957-5.View ArticlePubMedGoogle Scholar
- Rojas MR, Hagen C, Lucas WJ, Gilbertson RL: Exploiting chinks in the plant's armor: evolution and emergence of geminivirus. Annu Rev Phytopathol. 2005, 43: 361-394. 10.1146/annurev.phyto.43.040204.135939.View ArticlePubMedGoogle Scholar
- Seal SE, Jeger MJ, Bosch Van den F: Begomovirus evolution and disease management. Adv Virus Res. 2006, 67: 297-316. 10.1016/S0065-3527(06)67008-5.View ArticlePubMedGoogle Scholar
- Varma A, Malathi VG: Emerging geminivirus problems: A serious threat to crop production. Ann appl Biol. 2003, 142: 145-164. 10.1111/j.1744-7348.2003.tb00240.x.View ArticleGoogle Scholar
- Frischmuth T, Stanley J: Strategies for the control of geminivirus diseases. Semin Virol. 1993, 4: 329-337. 10.1006/smvy.1993.1031.View ArticleGoogle Scholar
- Polston JE, Anderson PK: The emergence of whitefly-transmitted geminiviruses in tomato in the western hemisphere. Plant Dis. 1997, 81: 1358-1369. 10.1094/PDIS.19184.108.40.2068.View ArticleGoogle Scholar
- Kvarnheden A, Lindblad M, Lindsten K, Valkonen JPT: Genetic diversity of Wheat dwarf virus. Arch Virol. 2002, 147: 205-261. 10.1007/s705-002-8313-x.View ArticlePubMedGoogle Scholar
- Schubert J, Habekuß A, Kazmaier K, Jeske H: Surveying cereal-infecting geminiviruses in Germany–Diagnostics and direct sequencing using rolling circle amplification. Virus Res. 2007, 127: 61-70. 10.1016/j.virusres.2007.03.018.View ArticlePubMedGoogle Scholar
- Ramsell JNE, Lemmetty A, Jonasson J, Andersson A, Sigvald R, Kvarnheden A: Sequence analyses of Wheat dwarf virus isolates from different hosts reveal low genetic diversity within the wheat strain. Plant Pathol. 2008, 57: 834-841. 10.1111/j.1365-3059.2008.01862.x.View ArticleGoogle Scholar
- Vacke J: Wheat dwarf virus disease. Biol Plant Praha. 1961, 3: 228-233. 10.1007/BF02933566.View ArticleGoogle Scholar
- Kõklû G, Ramsell JNE, Kvarnheden A: The complete genome sequence for a Turkish isolate of Wheat dwarf virus (WDV) from barley confirms the presence of two distinct WDV strains. Virus Genes. 2007, 34: 359-366. 10.1007/s11262-006-0029-0.View ArticlePubMedGoogle Scholar
- Huth W: Viruses of Graminae in Germany – a short overview. J Plant Dis Protection. 2000, 107: 406-14.Google Scholar
- Najar AKM, Makkouk H, Boudhir SG Kumari, Zarouk R, Bessai R, Othman FB: Viral diseases of cultivated legume and cereal crops in Tunisia. Phytopathol Mediterr. 2000, 39: 423-432.Google Scholar
- Lemmetty A, Huusela-Veistola E: First report of WDV in winter wheat in Finland. Plant Dis. 2005, 89: 912-10.1094/PD-89-0912B.View ArticleGoogle Scholar
- Xie J, Wang X, Liu Y, Peng Y, Zhou G: First report of the occurrence of WDV in wheat in China. Plant Dis. 2007, 91: 111-10.1094/PD-91-0111B.View ArticleGoogle Scholar
- Grass Phylogeny Working Group, Barker NP, Clark LG, Davis JI, Duvall MR, Guala GF, Hsiao C, Kellogg EA, Linder HP: Phylogeny and Subfamilial Classification of the Grasses (Poaceae). Ann Missouri Bot Garden. 2001, 88: 373-457. 10.2307/3298585.View ArticleGoogle Scholar
- Wolfe KH, Gouy MY, Yang W, Sharp PM, Li WH: Date of the monocot-dicot divergence estimated from chloroplast DNA sequence data. Proc Natl Acad Sci USA. 1989, 86: 6201-6205. 10.1073/pnas.86.16.6201.PubMed CentralView ArticlePubMedGoogle Scholar
- Vicentini A, Barber JC, Aliscioni SS, Giussani LM, Kellogg EA: The age of the grasses and clusters of origins of C4 photosynthesis. Global Change Biol. 2008, 14: 2963-2977. 10.1111/j.1365-2486.2008.01688.x.View ArticleGoogle Scholar
- Huang S, Sirikhachornkit A, Su X, Faris J, Gill B, Haselkorn R, Gornicki P: Genes encoding plastid acetyl-CoA carboxylase and 3-phosphoglycerate kinase of the Triticum_Aegilops complex and the evolutionary history of polyploid wheat. Proc Natl Acad Sci USA. 2002, 99: 8133-8138. 10.1073/pnas.072223799.PubMed CentralView ArticlePubMedGoogle Scholar
- Gaut BS: Evolutionary dynamics of grass genomes. New Phytologist. 2002, 154: 15-28.Google Scholar
- Watterson GA: On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975, 7: 256-276. 10.1016/0040-5809(75)90020-9.View ArticlePubMedGoogle Scholar
- Tajima F: Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983, 105: 437-460.PubMed CentralPubMedGoogle Scholar
- Varsani A, Shepherd DN, Monjane AL, Owor BE, Erdmann JB, Rybicki EP, Peterschmitt M, Briddon RW, Markham PG, Oluwafemi S, Windram OP, Lefeuvre P, Lett JM, Martin DP: Recombination, decreased host specificity and increased mobility may have driven the emergence of maize streak virus as an agricultural pathogen. J Gen Virol. 2008, 89: 2063-2074. 10.1099/vir.0.2008/003590-0.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin DP, Williamson C, Posada D: RDP2: recombination detection and analysis from sequence alignments. Bioinformatics. 2005, 21: 260-262. 10.1093/bioinformatics/bth490.View ArticlePubMedGoogle Scholar
- Gibbs MJ, Armstrong JS, Gibbs AJ: Sister-Scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000, 16: 573-582. 10.1093/bioinformatics/16.7.573.View ArticlePubMedGoogle Scholar
- Page RDM: Parallel phylogenies: reconstructing the history of host-parasite assemblages. Cladistics. 1995, 10: 155-173. 10.1111/j.1096-0031.1994.tb00170.x.View ArticleGoogle Scholar
- Charleston MA: Jungles: A new solution to the host/parasite phylogeny reconciliation problem. Math Biosci. 1998, 149: 191-223. 10.1016/S0025-5564(97)10012-8.View ArticlePubMedGoogle Scholar
- Fargette D, Konaté G, Fauquet C, Muller E, Peterschmitt M, Thresh JM: Molecular ecology and emergence of tropical plant viruses. Annu Rev Phytopathol. 2006, 44: 235-60. 10.1146/annurev.phyto.44.120705.104644.View ArticlePubMedGoogle Scholar
- Shepherd DN, Varsani A, Windram OP, Lefeuvre P, Monjane AL, Owor BE, Martin DP: Novel sugarcane streak and sugarcane streak Reunion mastreviruses from southern Africa and La Reunion. Arch Virol. 2008, 153: 605-609. 10.1007/s00705-007-0016-3.View ArticlePubMedGoogle Scholar
- Ge L, Zhang J, Zhou X, Li H: Genetic structure and population variability of Tomato Yellow Leaf Curl China Virus. J Virol. 2007, 81: 5902-5907. 10.1128/JVI.02431-06.PubMed CentralView ArticlePubMedGoogle Scholar
- French R, Stenger DC: Population structure within lineages of Wheat streak mosaic virus derived from a common founding event exhibits stochastic variation inconsistent with the deterministic quasi-species model. Virology. 2005, 343: 179-189. 10.1016/j.virol.2005.08.037.View ArticlePubMedGoogle Scholar
- Duffy S, Shackelton LA, Holmes EC: Rates of evolutionary change in viruses: patterns and determinants. Nat Rev Genet. 2008, 9 (4): 267-276. 10.1038/nrg2323.View ArticlePubMedGoogle Scholar
- Woolston CJ, Barker R, Gunn H, Boulton MI, Mullineaux PM: Agroinfection and nucleotide sequence of cloned wheat dwarf virus DNA. Plant Mole Biol. 1988, 11: 35-43. 10.1007/BF00016012.View ArticleGoogle Scholar
- Hall TA: BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.5c. Distributed by the author. 1993, Department of Genetics, University of Washington, Seattle, WA, USAGoogle Scholar
- Lazarowitz SG: Infectivity and complete nucleotide sequence of the genome of a South African isolate of maize streak virus. Nucleic Acids Res. 1988, 16: 229-249. 10.1093/nar/16.1.229.PubMed CentralView ArticlePubMedGoogle Scholar
- Rozas J, Juan C, Sánchez-Dei MX, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Se-Al. [http://tree.bio.ed.ac.uk/software/seal/]
- Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. 1998, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Posada D: Selecting models of evolution. The Phylogenetic Handbook. Edited by: Salemi M, Vandamme A-M. 2003, Cambridge: Cambridge University Press, 256-282.Google Scholar
- Padidam M, Sawyer S, Fauquet CM: Possible emergence of new geminiviruses by frequent recombination. Virology. 1999, 265: 218-225. 10.1006/viro.1999.0056.View ArticlePubMedGoogle Scholar
- Martin DP, Posada D, Crandall KA, Williamson C: A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retrovir. 2005, 21: 98-102. 10.1089/aid.2005.21.98.View ArticlePubMedGoogle Scholar
- Maynard Smith J: Analysing the mosaic structure of genes. J Mol Evol. 1992, 34: 126-129.Google Scholar
- Weiller GF: Phylogenetic profiles: A graphical method for detecting genetic recombinations in homologous sequences. Mol Biol Evol. 1998, 15: 326-335.View ArticlePubMedGoogle Scholar
- Weiller GF, Gibbs A: Diplomo: The tools for a new type of evolutionary analysis. Comp Appl Biosci. 1995, 11: 535-540.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.