HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Free via Open Access: OA OA Abstract Full Text (PDF) Supplementary Table OA All Versions of this Article: vir.0.015537-0v1 90/12/3066    most recent Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via CrossRef Google Scholar Articles by Harkins, G. W. Articles by Varsani, A. PubMed PubMed Citation Articles by Harkins, G. W. Articles by Varsani, A. Agricola Articles by Harkins, G. W. Articles by Varsani, A.

Dating the origins of the maize-adapted strain of maize streak virus, MSV-A

¹ South African National Bioinformatics Institute, University of the Western Cape, Private Bag X17, Bellville 7535, South Africa
² Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Observatory, Cape Town, South Africa
³ Centre for High-Performance Computing, Rosebank, Cape Town, South Africa
⁴ Department of Ecology, Evolution and Natural Resources, School of Environmental and Biological Sciences, Rutgers University, New Brunswick, NJ 08901 USA
⁵ Department of Molecular and Cell Biology, University of Cape Town, Rondebosch, Cape Town 7701, South Africa
⁶ Warwick Systems Biology Centre, University of Warwick, Wellesbourne CV35 9EF, UK
⁷ South African Sugarcane Research Institute, Mount Edgecombe, KwaZulu Natal, South Africa
⁸ Crop Protection, ARC-Grain Crops Institute, Potchefstroom 2520, South Africa
⁹ CIRAD, UMR BGPI, TA A54/K, Campus International de Baillarguet, 34398 Montpellier Cedex 5, France
¹⁰ Electron Microscope Unit, University of Cape Town, Private Bag, Rondebosch 7701, South Africa
¹¹ School of Biological Sciences, University of Canterbury, Private Bag 4800, Christchurch, New Zealand

Correspondence
Darren P. Martin
darrin.martin{at}uct.ac.za

	ABSTRACT

TOP ABSTRACT INTRODUCTION METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
Maize streak virus (MSV), which causes maize streak disease (MSD), is one of the most serious biotic threats to African food security. Here, we use whole MSV genomes sampled over 30 years to estimate the dates of key evolutionary events in the 500 year association of MSV and maize. The substitution rates implied by our analyses agree closely with those estimated previously in controlled MSV evolution experiments, and we use them to infer the date when the maize-adapted strain, MSV-A, was generated by recombination between two grass-adapted MSV strains. Our results indicate that this recombination event occurred in the mid-1800s, 20 years before the first credible reports of MSD in South Africa and centuries after the introduction of maize to the continent in the early 1500s. This suggests a causal link between MSV recombination and the emergence of MSV-A as a serious pathogen of maize.

The GenBank/EMBL/DDBJ accession numbers for the genome sequences of the MSV isolates analysed in this study are FJ882089–FJ882124 and FJ882126–FJ882149.

A supplementary table, listing sampling locations, dates and full genome sequence accession numbers of MSV isolates examined in this study, is available with the online version of this paper.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
Maize streak virus (MSV; family Geminiviridae, genus Mastrevirus) is amongst the most important plant pathogens threatening food security in Africa (Martin & Shepherd, 2009). Whereas MSV is best known for causing occasionally devastating maize streak disease (MSD) epidemics throughout sub-Saharan Africa, it also infects over 80 other species, primarily wild grasses in the family Poaceae (Damsteegt, 1983; ICTVdb Management, 2006; Konate & Traore, 1992).

Of the 11 known MSV strains (Varsani et al., 2008) only one, MSV-A, is usually found infecting maize plants (Briddon et al., 1994; Hughes et al., 1992; Owor et al., 2007a; Varsani et al., 2008; Willment et al., 2001). MSV-A is also the only MSV strain that causes serious disease in maize (Martin et al., 1999, 2001). Other MSV strains (MSV-B to -K) are generally found infecting wild grass species in the genera Digitaria (e.g. MSV-B), Urochloa (e.g. MSV-F and -G) or Setaria (e.g. MSV-C, -H and -K; Varsani et al., 2008; Rybicki et al., 1998). ‘Grass-adapted’ strains such as MSV-B, -C, -D and -E cannot symptomatically infect any but the most MSV-susceptible maize genotypes (Martin et al., 2001; Schnippenkoetter et al., 2001; Willment et al., 2002).

MSD was first described in 1901 in southern Africa, where it had affected maize since at least the 1870s (Fuller, 1901), suggesting that maize-adapted MSV-A genotypes may have already evolved by the mid-1800s. Maize was only introduced to West Africa by Portuguese traders in the early 1500s and it reached southern Africa in the mid-1600s with the Dutch East India Company (McCann, 2001). It is possible that a grass-adapted progenitor of MSV-A circulating in African grasses during the 16th or 17th centuries was serendipitously capable of causing severe infections in this exotic new host. A second possibility is that grass-adapted MSV-A ancestors progressively adapted to maize over the previous four centuries. A third possibility is that the capacity to cause severe infections in maize arose instantaneously along with the genesis of MSV-A through the chance genetic recombination of two grass-adapted MSV strains. All MSV-A genomes sampled to date are apparently descended from an inter-strain recombinant that derived most of its movement (mp) and coat (cp) protein genes from an MSV-B-like parental virus, and the remainder of its genome from a parent resembling an ancestral MSV-F and -G-like virus (Varsani et al., 2008). Recombination is extremely common in geminiviruses (Lefeuvre et al., 2009) and has been directly implicated in the emergence of new crop diseases (Fondong et al., 2000; Pita et al., 2001; Monci et al., 2002; García-Andrés et al., 2007a, b). It is particularly plausible that the recombination event that generated MSV-A triggered the emergence of MSV as a serious agricultural pathogen because the event involved the transfer of a mostly self-contained movement and encapsidation ‘module’ (Martin et al., 2005; van der Walt et al., 2008a) that is known to contain the main pathogenicity determinants of MSV-A (Martin & Rybicki, 2002; van der Walt et al., 2008a).

To rationally assess these alternative MSV-A origin hypotheses, it would be useful to know, firstly, when MSV-A diverged from its grass-adapted relatives and, secondly, when the recombination event that produced the ancestral MSV-A occurred. All that is required to meaningfully estimate these dates is a reliable estimate of the MSV evolution rate.

Analyses of mutation frequencies and highly adaptive reversion mutations over very short time periods (<60 days) have indicated that MSV experiences basal mutation frequencies in excess of 2x10^–3 mutations per site per year (Shepherd et al., 2005, 2006; van der Walt et al., 2009). Longer term monitoring of small MSV populations in individual hosts over 1–6 years (Isnard et al., 1998; Harkins et al., 2009; van der Walt et al., 2008b) has indicated that the MSV substitution rate (the rate at which mutations accumulate in MSV populations over time) is in the range of 2–7x10^–4 substitutions per site per year (subs per site per year) i.e. approximately three- to tenfold lower than the basal mutation rate.

However, it is presently not known how MSV substitution rate estimates from small experimental populations compare with those from large continent-wide natural populations. Apparent co-divergence of some mastreviruses with their hosts has led to speculation that these viruses have substitution rates of only 10^–8 subs per site per year in nature (Wu et al., 2008). However, analyses of temporally structured tomato yellow leaf curl virus and East African cassava mosaic virus datasets sampled from nature have indicated global evolution rates for these geminivirus species of between 10^–3 and 10^–5 subs per site per year (Duffy & Holmes, 2008, 2009), rates that are broadly concordant with those determined experimentally for both MSV and some tomato yellow leaf curl disease-causing geminiviruses (Ge et al., 2007; Harkins et al., 2009; Urbino et al., 2008).

Here, we describe the application of a Bayesian coalescent-based approach to study temporally structured MSV-A datasets sampled from nature. The methodology used is similar to that recently applied to other important rapidly evolving plant viruses in the families Potyviridae (Gibbs et al., 2008) and Sobemoviridae (Fargette et al., 2008). Specifically, we infer both the date when the last common ancestor of all currently sampled MSV-A genomes existed and the earliest probable date when the recombination event that generated MSV-A might have occurred. With these estimates in hand, we assess three potential hypotheses of the origin of MSV-A.

	METHODS

TOP ABSTRACT INTRODUCTION METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
Assembly of temporally structured MSV-A datasets.
We obtained 41 MSV samples archived in freezers between 1979 and 1999 at the Institut de Biologie Moleculaire et Cellulaire du CNRS, Montpellier, France, and the University of Cape Town, South Africa. We cloned and sequenced full genomes from these and an additional 19 genomes sampled between 2005 and 2008 using previously described methods (Owor et al., 2007b; Shepherd et al., 2008; see Supplementary Table S1, available in JGV Online). We combined these sequences with 64 dated sequences taken from GenBank to produce a temporally structured MSV dataset comprising 124 full genome sequences. All of these sequences had been sampled directly from natural environments and none had been passaged for prolonged periods under glasshouse conditions prior to cloning and sequencing. These genomes were aligned together with representative sequences from each of the 10 other known MSV strains using POA (Grasso & Lee, 2004) with manual editing of alignments using MEGA4 (Tamura et al., 2007).

Construction of mostly recombination-free datasets.
As recombination can confound substitution rate estimation (Schierup & Hein, 2000), the alignment was analysed, as described by Heath et al. (2006), for the presence of recombination using a battery of recombination analysis methods implemented in the program RDP3. Two mostly recombination-free MSV datasets were assembled: (i) comprised 104 sequences each containing 900 mp and cp gene nucleotides (dataset MPCP) and (ii) comprised 125 MSV-A sequences and six MSV-B sequences each containing 450 cp gene nucleotides (dataset CP). These datasets were only ‘mostly recombination-free’ because despite removing all detectable evidence of recombination, it is very probable that, due to limitations of the methods applied, other very difficult to detect trace recombination signals remained (van der Walt et al., 2009). For the full genome dataset, we removed all nucleotides from predominantly MSV-A sequences that were apparently derived through recombination from other strains. This involved the removal of two short sequence tracts: one of 85 nt from 48 sequences (corresponding to genome coordinates 1346–1431 in MSV-A [Zm-MatA-1994]) and the other of 112 nt from six sequences (corresponding to genome coordinates 155–267 in MSV-A [Zm-MatA-1994]). Except for six MSV-B sequences left in the CP dataset, non-MSV-A sequences were removed from all the datasets following recombination analysis. MSV-B sequences were retained in the CP dataset so that we could later use it to date the recombination event that yielded MSV-A (the tract of MSV-A sequence retained in this dataset was derived from its MSV-B parent). All analysed alignments are available on request.

Estimation of evolution rates and dating of ancestral sequences.
For all three datasets (whole genome, MPCP and CP) a co-estimate of the nucleotide substitution model [Hasegawa, Kishino and Yano with a calculated proportion of invariable sites with a gamma distribution (HKY+ I+G) or general time reversible (GTR)+I+G] parameters, phylogeny and time to the most recent common ancestor (tMRCA) was obtained for the various major MSV-A lineages using the Bayesian Markov chain Monte Carlo (MCMC) method implemented in BEAST v1.4.8 (Drummond & Rambaut, 2007). The HKY+I+G₄ and GTR+I+G₄ nucleotide substitution models both yielded similar results and we only report here the analyses using HKY+I+G_4. For each dataset, six different coalescent models were tested including both parametric (constant population size, exponential growth) and non-parametric [Bayesian skyline plot (BSP)] models with both a strict and a relaxed (uncorrelated LogNormal prior) molecular clock.

For each evolutionary model, between 5 and 15 independent runs, each with between 2.0x10⁷ and 5.0x10⁷ steps in the Markov chain were performed using BEAST and checked for convergence using Tracer v1.4 (Drummond & Rambaut 2007). The effective sample sizes (ESSs) were always >200, indicating sufficient mixing of the Markov chain and parameter sampling. When similar results were produced from independent runs of the Markov chain, the log files were combined with the program LogCombiner v1.4.7 available in the BEAST package (Drummond & Rambaut 2007).

Identification of the best-fit clock and demographic models was achieved by calculating the Bayes factor (BF), which is the ratio of the marginal likelihoods of the two models being compared (Kass & Raftery 1995; Suchard et al., 2001). This method permits the comparison of non-nested models (such as the non-parametric BSP versus the parametric constant or exponential growth demographic models) that cannot be validly compared using the mean log posterior probabilities.

	RESULTS AND DISCUSSION

TOP ABSTRACT INTRODUCTION METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
Determining the MSV-A substitution rate
The mean substitution rates estimated for the full genome and mostly recombination-free MPCP datasets with the various clock (strict or relaxed) and demographic (constant population size, exponential growth or BSP) models tested were broadly similar and ranged from 1.8x10^–4 subs per site per year (for the full genome dataset analysed with the strict clock+BSP model) to 3.9x10^–4 subs per site per year (for the full genome dataset analysed with the relaxed clock+exponential growth model). For the mostly recombination-free CP dataset, which contained six MSV-B sequences in addition to MSV-A sequences, estimated mean substitution rates were a little higher, ranging from 4.9x10^–4 subs per site per year (strict clock+constant population size model) to 7.0x10^–4 subs per site per year (relaxed clock+BSP model). Also, for each individual dataset, 95 % highest probability density (HPD) intervals were largely overlapping irrespective of the demographic and clock model combinations tested. For all three datasets, the BF tests significantly supported models enforcing relaxed over strict molecular clock models (Table 1; note that 2lnBF scores >3 reflect a significantly better fit of model H₁ versus H₀). Consistent with these results, the estimated values of the coefficient of variation parameter from the relaxed clock runs (n=9 models) ranged from 0.54 to 0.92, suggesting the presence of significant substitution rate heterogeneity among lineages (Drummond et al., 2006). The BF tests also generally supported exponential growth and/or constant population size models over BSP models. For the whole genome and CP datasets, the exponential growth and constant population size relaxed clock models fitted the data equally well, while for the MPCP dataset, both the exponential growth and BSP demographic models fitted the data equally well. We noticed, however, that in the maximum clade credibility trees that we constructed for the CP dataset using the relaxed clock+exponential growth/BSP models, the six MSV-B isolates included in the CP dataset were not basal to the MSV-A sequences. Despite the CP exponential growth and constant population size models being approximately equally supported by BF tests, failure to recover a credible tree under the exponential growth model disqualified this model from further consideration in analyses of the CP dataset.

View larger version (27K): [in this window] [in a new window]   Fig. 1. Estimated substitution rates used to date key events in the evolution of MSV-A. (a) Mean substitution rates estimated with temporally structured MSV-A datasets sampled from nature using the best-fit demographic and clock models (highlighted in black) (constant population size+relaxed LogNormal molecular clock for the full genome and CP datasets and exponential population growth+relaxed LogNormal molecular clock for the MPCP dataset) compared with those estimated either experimentally for MSV (van der Walt et al., 2008b; Harkins et al., 2009), sugarcane streak Réunion virus (SSRV) (Harkins et al., 2009), tomato yellow leaf curl China virus (TYLCCV) (Ge et al., 2007) and tomato yellow leaf curl virus (TYLCV) (Urbino et al., 2008) or from temporally structured genome sequence datasets sampled from nature [East African cassava mosaic virus (EACMV) and TYLCV; Duffy & Holmes, 2008, 2009]. Squares indicate most probable substitution rates, vertical bars indicate the 95 % highest probability densities of the substitution rate estimates. (b) Most probable estimates (squares) of dates when the last common MSV-A ancestor existed (estimated with the full genome and MPCP datasets) and when the MSV-A/MSV-B lineages split (indicative of the earliest date when MSV-A could have arisen as estimated with the CP dataset). Vertical lines indicate 95 % HPD intervals, the date estimates determined with the best-fit models (see Table 1) are highlighted in black. The three estimates with the best-fit models that we focus on the most in the text are indicated by asterisks (*). For clock models: S, strict clock; R, relaxed clock. For demographic models: C, constant population size; E, exponential population growth; B, BSP (allowing for varying phases of population growth/decline). Whereas some recombination is detectable in the full genome dataset none was detectable in either the MPCP or CP datasets. The dates of key historical events are indicated for comparative purposes. Dotted lines indicate that a portion of some of the branches and the timeline have been removed for display purposes.

Despite these broad confidence intervals, this analysis implied that the last common ancestor of all currently sampled MSV-A viruses only arose at least 200 years after maize was first introduced into West Africa by the Portuguese and at least 100 years after it was introduced to southern Africa by the Dutch (McCann, 2001). As there are narrower credibility intervals attained with the recombination-free MPCP dataset (Fig. 1), it is very probable that the last common ancestor of all the MSV-A viruses occurred more recently than 1791 [the lowest 95 % HPD limit for any of the clock and demographic model combinations used to analyse this dataset (constant population size+relaxed clock model)]. Considering the whole genome and MPCP analyses together, this date was most probably between 1839 (full genome+constant population size+relaxed clock) and 1905 (MPCP+exponential growth+relaxed clock); both of these are within 62 years of when MSV was first unambiguously described by Fuller (1901).

Determining the earliest date when MSV-A could have emerged
While the most recent common ancestor of all currently sampled MSV-A lineages (i.e. the root of the MSV-A phylogenetic tree) probably existed at some time between the mid-1800s and early 1900s, these estimated dates do not necessarily indicate when MSV-A emerged as a serious maize pathogen. The date of this emergence is expected to be older, and to map somewhere on the branch of the phylogenetic tree that separates the MSV-A clade from the other MSV strains. In fact, Fuller (1901) cited reports indicating that something resembling MSD was endemic in South Africa as early as the 1870s, indicating that the origin of MSV-A probably precedes this date.

The earliest possible date for the emergence of MSV-A was when it diverged from its nearest grass-adapted relatives. Estimating this date is paradoxically helped by the fact that all currently sampled MSV-As are apparently recombinants of MSV-B- and MSV-F-like viruses (Varsani et al., 2008). The coat protein gene (cp) of MSV-A very closely resembles that of MSV-B, whereas the rest of its genome is derived from a currently unsampled strain most closely resembling MSV-G- or MSV-F-like viruses. The problem of inferring the earliest possible date of emergence of MSV-A is therefore equivalent to inferring the earliest date when the recombination event that yielded MSV-A could have occurred. This is, in turn, equivalent to using the cp region of MSV-A and MSV-B viruses to infer the time when the last common ancestor of these strains existed.

By applying an approach to the CP dataset that is similar to that recently used to date the origin of a human immunodeficiency virus circulating recombinant form (Tee et al., 2009), we determined that the last common MSV-A/MSV-B ancestor in the cp region most probably existed in 1853 (95 % HPD=17451941; relaxed clock+constant population size model). This indicates that MSV-A and MSV-B share a common ancestor that existed at approximately the same time as the last common ancestor of the MSV-A viruses. In fact, the last common MSV-A/MSV-B ancestor probably did not exist more that 50 years prior to the last common MSV-A ancestor. Considering the lower bounds of 95 % HPD intervals obtained using the best-fit model (relaxed clock+constant population size), this prototypical recombinant MSV-A virus almost certainly existed at some time after 1744 – at least 200 years after maize was first introduced to West Africa and at most 130 years before MSV-A-like infection symptoms were first noted in 1870 by maize farmers in South Africa. It is very interesting that our analysis indicates that 1853 is the most probable date of the MSV-A/MSV-B split, since this implies that the recombination event that produced the first MSV-A genome most probably occurred within 20 years of the first reports of MSD.

The spread of MSV-A throughout Africa
Given all the available data, we formulated a most-probable timeline for the emergence and spread of MSV-A. Although the recombination event that yielded MSV-A almost certainly occurred between the 1740s and 1870s, it most probably occurred in the 1850s. In the 1870s, symptoms resembling those produced by MSV-A were reported in South Africa, marking this as the latest date at which MSV-A-like viruses might have emerged. Given the high degree of overlap between the estimated dates of the last common MSV-A ancestor and the last common MSV-A/MSV-B ancestor, it is probable that this earliest MSV-A sequence was either reasonably well adapted to maize (possibly due to the initial recombination event that produced it having been selectively favoured in this host) or that its adaptation to maize probably took less than 130 years (the time between the earliest credible date of the MSV-A/MSV-B split and 1870).

At some time between the 1870s and 1951 (the upper 95 % HPD date considering all best-fit models) and most probably around 1904, an MSV-A variant existed that was to become the last common ancestor of all the MSV-A viruses that have been sampled since 1979 (node 1 in Fig. 2). Whereas the MSV-A₆ lineage of the descendents of this ancestral virus colonized the Indian Ocean islands off the coast of Africa, another lineage diversified after 1904 (95 % HPD=1861–1938; node 2 in Fig. 2) to yield all the MSV-A variants that are currently found throughout mainland Africa.

View larger version (21K): [in this window] [in a new window]   Fig. 2. The maximum clade credibility MSV-A phylogenetic tree (for the full genome dataset) recovered under a best-fit demographic (constant population size) and clock (LogNormal relaxed) model. Viruses sampled from southern Africa, East Africa, West Africa and the island of La Réunion are indicated by orange, blue, green and pink branches, respectively. Nodes labelled 1 to 6 indicate the dated most recent common ancestors of major MSV lineages referred to in the text. Branches recovered in >99 % of sampled trees are indicated by filled circles, whereas those recovered in 90–99 % of sampled trees are indicated by empty circles. Branches recovered in less than 80 % of sampled trees were collapsed. An approximate timeline is given beneath the tree and major MSV-A lineages (A1, A2, A3, A4 and A6, as defined by Owor et al., 2007a) are labelled on the right.

On mainland Africa, MSV-A spread and differentiated into a southern African lineage (MSV-A₄) sometime around 1924 (95 % HPD=1893–1949; node 3 in Fig. 2), into an East African lineage (MSV-A₃) sometime around 1967 (95 % HPD=1939–1980; node 4 in Fig. 2) and a West African lineage (MSV-A₂) sometime around 1969 (95 % HPD=1940–1981; node 5 in Fig. 2). The last common ancestor of the economically most relevant pan-African lineage, MSV-A₁, probably existed sometime around 1922 (95 % HPD=1882–1953; node 6 in Fig. 2). Unlike the other MSV-A lineages, the mixing of MSV-A₁ isolates sampled from East, West and southern Africa within several well-supported clades in the MSV-A phylogenetic tree (Fig. 2), suggests that in the past 80 years, there were multiple waves of continent-wide MSV-A₁ dispersal.

	CONCLUSIONS

TOP ABSTRACT INTRODUCTION METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
While our results are inconsistent with the first origin hypothesis that MSV-A already existed in Africa at the time when maize was first brought to the continent in the early 1500s, it remains possible that, after MSV-A arose by recombination during the 18th or 19th centuries, it required up to 130 years to adapt to maize before being noticed by farmers in the 1870s (the second origin hypothesis). We have, however, also shown that the recombination event that yielded the MSV-A prototype most probably occurred within 20 years of MSD being first reported. Our data are therefore clearly also consistent with, and strongly supportive of, the third origin hypothesis that the adaptation of MSV-A to maize was very rapid and that the emergence of MSV-A as a serious agricultural pathogen was temporally, and therefore possibly causally, linked with the recombination event that produced this strain.

	ACKNOWLEDGEMENTS

 
D. P. M. is supported by the Wellcome Trust. A. V. was supported by the Carnegie Corporation of New York. The research was funded by the National Research Foundation of South Africa.

	REFERENCES

TOP ABSTRACT INTRODUCTION METHODS RESULTS AND DISCUSSION CONCLUSIONS REFERENCES

 
Briddon, R. W., Lunness, P., Chamberlin, L. C. & Markham, P. G. (1994). Analysis of the genetic variability of maize streak virus. Virus Genes 9, 93–100.[CrossRef][Medline]

Damsteegt, V. D. (1983). Maize streak virus: I. Host range and vulnerability of maize germ plasm. Plant Dis 67, 734–737.[CrossRef]

Drummond, A. J. & Rambaut, A. (2007). BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7, 214[CrossRef][Medline]

Drummond, A. J., Ho, S. Y., Phillips, M. J. & Rambaut, A. (2006). Relaxed phylogenetics and dating with confidence. PLoS Biol 4, e88[CrossRef][Medline]

Duffy, S. & Holmes, E. C. (2008). Phylogenetic evidence for rapid rates of molecular evolution in the single-stranded DNA begomovirus Tomato yellow leaf curl virus (TYLCV). J Virol 82, 957–965.[Abstract/Free Full Text]

Duffy, S. & Holmes, E. C. (2009). Validation of high rates of nucleotide substitution in geminiviruses: phylogenetic evidence from East African cassava mosaic viruses. J Gen Virol 90, 1539–1547.[Abstract/Free Full Text]

Fargette, D., Pinel-Galzi, A., Sérémé, D., Lacombe, S., Hébrard, E., Traoré, O. & Konaté, G. (2008). Diversification of rice yellow mottle virus and related viruses spans the history of agriculture from the neolithic to the present. PLoS Pathog 4, e1000125[CrossRef][Medline]

Fondong, V. N., Pita, J. S., Rey, M. E., de Kochko, A., Beachy, R. N. & Fauquet, C. M. (2000). Evidence of synergism between African cassava mosaic virus and a new double-recombinant geminivirus infecting cassava in Cameroon. J Gen Virol 81, 287–297.[Abstract/Free Full Text]

Fuller, C. (1901). Mealie Variegation. First Report of the Government Entomologist, 1899–1900. Pietermaritzburg, Natal: P. Davis & Sons, Government Printers.

García-Andrés, S., Tomas, D. M., Sanchez-Campos, S., Navas-Castillo, J. & Moriones, E. (2007a). Frequent occurrence of recombinants in mixed infections of tomato yellow leaf curl disease-associated begomoviruses. Virology 365, 210–219.[CrossRef][Medline]

García-Andrés, S., Accotto, G. P., Navas-Castillo, J. & Moriones, E. (2007b). Founder effect, plant host, and recombination shape the emergent population of begomoviruses that cause the tomato yellow leaf curl disease in the Mediterranean basin. Virology 359, 302–312.[CrossRef][Medline]

Ge, L., Zhang, J., Zhou, X. & Li, H. (2007). Genetic structure and population variability of tomato yellow leaf curl China virus. J Virol 81, 5902–5907.[Abstract/Free Full Text]

Gibbs, A. J., Ohshima, K., Phillips, M. J. & Gibbs, M. J. (2008). The prehistory of potyviruses: their initial radiation was during the dawn of agriculture. PLoS One 3, e2523[CrossRef][Medline]

Grasso, C. & Lee, C. (2004). Combining partial order alignment and progressive multiple sequence alignment increases alignment speed and scalability to very large alignment problems. Bioinformatics 20, 1546–1556.[Abstract/Free Full Text]

Harkins, G. W., Delport, W., Duffy, S., Wood, N., Monjane, A. L., Owor, B. E., Donaldson, L., Saumtally, S., Triton, G. & other authors (2009). Experimental evidence indicating that mastreviruses probably did not co-diverge with their hosts. Virol J 6, 104[CrossRef][Medline]

Heath, L., van der Walt, E., Varsani, A. & Martin, D. P. (2006). Recombination patterns in aphthoviruses mirror those found in other picornaviruses. J Virol 80, 11827–11832.[Abstract/Free Full Text]

Hughes, F. L., Rybicki, E. P. & von Wechmar, M. B. (1992). Genome typing of southern African subgroup 1 geminiviruses. J Gen Virol 73, 1031–1040.[Abstract/Free Full Text]

ICTVdB Management (2006). 00.029.0.01.001. Maize streak virus. In ICTVdB – The Universal Virus Database, version 4. Edited by C. Büchen-Osmond. http://www.ncbi.nlm.nih.gov/ICTVdb/ICTVdB/

Isnard, M., Granier, M., Frutos, R., Reynaud, B. & Peterschmitt, M. (1998). Quasispecies nature of three Maize streak virus isolates obtained through different modes of selection from a population used to assess response to infection of maize cultivars. J Gen Virol 79, 3091–3099.[Abstract]

Kass, R. E. & Raftery, A. E. (1995). Bayes Factors. J Am Stat Assoc 90, 773–795.[CrossRef]

Konate, G. & Traore, O. (1992). Reservoir hosts of maize streak virus (MSV) in the Sudan Sahel zone: identification and spatio temporal distribution. Phytoprotection 73, 111–117.

Kopp, A. (1930). Les maladies des plantes à La Réunion. Rev Bot Appl Agric Trop 105, 1–7 (in French).

Kopp, A. & d'Emmerez de Charmoy, D. (1931). Nouvelles constatations sur les maladies à virus de la canne à sucre et du maïs. C R Acad Sciences 193, 876–878 (in French).

Kopp, A. & d'Emmerez de Charmoy, D. (1932). Observations nouvelles concernant la mosaïque de la canne à sucre et le streak du maïs à la Réunion. Station Agronomique Réunion. Travaux Techniques 3, 1–10 (in French).

Lefeuvre, P., Lett, J. M., Varsani, A. & Martin, D. P. (2009). Widely conserved recombination patterns among single-stranded DNA viruses. J Virol 83, 2697–2707.[Abstract/Free Full Text]

Martin, D. P. & Rybicki, E. P. (2002). Investigation of maize streak virus pathogenicity determinants using chimaeric genomes. Virology 300, 180–188.[CrossRef][Medline]

Martin, D. P. & Shepherd, D. N. (2009). The epidemiology, economic impact and control of maize streak disease. Food Security 1, 305–315.[CrossRef]

Martin, D. P., Willment, J. A. & Rybicki, E. P. (1999). Evaluation of maize streak virus pathogenicity in differentially resistant Zea mays genotypes. Phytopathology 89, 695–700.[CrossRef][Medline]

Martin, D. P., Willment, J. A., Billharz, R., Velders, R., Odhiambo, B., Njuguna, J., James, D. & Rybicki, E. P. (2001). Sequence diversity and virulence in Zea mays of Maize streak virus isolates. Virology 288, 247–255.[CrossRef][Medline]

Martin, D. P., van der Walt, E., Posada, D. & Rybicki, E. P. (2005). The evolutionary value of recombination is constrained by genome modularity. PLoS Genet 1, e51[CrossRef][Medline]

McCann, J. (2001). Maize and grace: history, corn, and Africa's new landscapes, 1500–1999. Comp Stud Soc Hist 43, 246–272.[CrossRef]

Monci, F., Sánchez-Campos, S., Navas-Castillo, J. & Moriones, E. (2002). A natural recombinant between the geminiviruses Tomato yellow leaf curl Sardinia virus and Tomato yellow leaf curl virus exhibits a novel pathogenic phenotype and is becoming prevalent in Spanish populations. Virology 303, 317–326.[CrossRef][Medline]

Owor, B. E., Martin, D. P., Shepherd, D. N., Edema, R., Monjane, A. L., Rybicki, E. P., Thomson, J. A. & Varsani, A. (2007a). Genetic analysis of maize streak virus isolates from Uganda reveals widespread distribution of a recombinant variant. J Gen Virol 88, 3154–3165.[Abstract/Free Full Text]

Owor, B. E., Shepherd, D. N., Taylor, N. J., Edema, R., Monjane, A. L., Thomson, J. A., Martin, D. P. & Varsani, A. (2007b). Successful application of FTA Classic Card technology and use of bacteriophage 29 DNA polymerase for large-scale field sampling and cloning of complete maize streak virus genomes. J Virol Methods 140, 100–105.[CrossRef][Medline]

Pita, J. S., Fondong, V. N., Sangaré, A., Otim-Nape, G. W., Ogwal, S. & Fauquet, C. M. (2001). Recombination, pseudorecombination and synergism of geminiviruses are determinant keys to the epidemic of severe cassava mosaic disease in Uganda. J Gen Virol 82, 655–665.[Abstract/Free Full Text]

Rybicki, E. P., Dennis, S., Napier, G. & Hughes, F. L. (1998). Novel grass and wheat strains of maize streak virus detected by DNA amplification and sequencing. Mol Plant Pathol On-Line http://www.bspp.org.uk/mppol/1998/0929rybicki/

Schierup, M. H. & Hein, J. (2000). Recombination and the molecular clock. Mol Biol Evol 17, 1578–1579.[Free Full Text]

Schnippenkoetter, W. H., Martin, D. P., Hughes, F., Fyvie, M., Willment, J. A., James, D., von Wechmar, B. & Rybicki, E. P. (2001). The biological and genomic characterisation of three mastreviruses. Arch Virol 146, 1075–1088.[CrossRef][Medline]

Shepherd, D. N., Martin, D. P., McGivern, D. R., Boulton, M. I., Thomson, J. A. & Rybicki, E. P. (2005). A three-nucleotide mutation altering the maize streak virus Rep pRBR-interaction motif reduces symptom severity in maize and partially reverts at high frequency without restoring pRBR-Rep binding. J Gen Virol 86, 803–813.[Abstract/Free Full Text]

Shepherd, D. N., Martin, D. P., Varsani, A., Thomson, J. A., Rybicki, E. P. & Klump, H. H. (2006). Restoration of native folding of single-stranded DNA sequences through reverse mutations: an indication of a new epigenetic mechanism. Arch Biochem Biophys 453, 108–122.[CrossRef][Medline]

Shepherd, D. N., Martin, D. P., Lefeuvre, P., Monjane, A. L., Owor, B. E., Rybicki, E. P. & Varsani, A. (2008). A protocol for the rapid isolation of full geminivirus genomes from dried plant tissue. J Virol Methods 149, 97–102.[CrossRef][Medline]

Suchard, M. A., Weiss, R. E. & Sinsheimer, J. S. (2001). Bayesian selection of continuous-time Markov chain evolutionary models. Mol Biol Evol 18, 1001–1013.[Abstract/Free Full Text]

Tamura, K., Dudley, J., Nei, M. & Kumar, S. (2007). MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24, 1596–1599.[Abstract/Free Full Text]

Tee, K. K., Pybus, O. G., Parker, J., Ng, K. P., Kamarulzaman, A. & Takebe, Y. (2009). Estimating the date of origin of an HIV-1 circulating recombinant form. Virology 387, 229–234.[CrossRef][Medline]

Urbino, C., Thébaud, G., Granier, M., Blanc, S. & Peterschmitt, M. (2008). A novel cloning strategy for isolating, genotyping and phenotyping genetic variants of geminiviruses. Virol J 5, 135[CrossRef][Medline]

van der Walt, E., Palmer, K. E., Martin, D. P. & Rybicki, E. P. (2008a). Viable chimaeric viruses confirm the biological importance of sequence specific maize streak virus movement protein and coat protein interactions. Virol J 5, 61[CrossRef][Medline]

van der Walt, E., Martin, D. P., Varsani, A., Polston, J. E. & Rybicki, E. P. (2008b). Experimental observations of rapid maize streak virus evolution reveal a strand-specific nucleotide substitution bias. Virol J 5, 104[CrossRef][Medline]

van der Walt, E., Rybicki, E. P., Varsani, A., Polston, J. E., Billharz, R., Donaldson, L., Monjane, A. L. & Martin, D. P. (2009). Rapid host adaptation by extensive recombination. J Gen Virol 90, 734–746.[Abstract/Free Full Text]

Varsani, A., Shepherd, D. N., Monjane, A. L., Owor, B. E., Erdmann, J. B., Rybicki, E. P., Peterschmitt, M., Briddon, R. W., Markham, P. G. & other authors (2008). Recombination, decreased host specificity and increased mobility may have driven the emergence of maize streak virus as an agricultural pathogen. J Gen Virol 89, 2063–2074.[Abstract/Free Full Text]

Willment, J. A., Martin, D. P. & Rybicki, E. P. (2001). Analysis of the diversity of African streak mastreviruses using PCR-generated RFLPs and partial sequence data. J Virol Methods 93, 75–87.[CrossRef][Medline]

Willment, J. A., Martin, D. P., Van der Walt, E. & Rybicki, E. P. (2002). Biological and genomic sequence characterization of maize streak virus isolates from wheat. Phytopathology 92, 81–86.[CrossRef][Medline]

Wu, B., Melcher, U., Guo, X., Wang, X., Fan, L. & Zhou, G. (2008). Assessment of codivergence of Mastreviruses with their plant hosts. BMC Evol Biol 8, 335[CrossRef][Medline]

Received 30 July 2009; accepted 14 August 2009.

This Article Free via Open Access: OA OA Abstract Full Text (PDF) Supplementary Table OA All Versions of this Article: vir.0.015537-0v1 90/12/3066    most recent Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via CrossRef Google Scholar Articles by Harkins, G. W. Articles by Varsani, A. PubMed PubMed Citation Articles by Harkins, G. W. Articles by Varsani, A. Agricola Articles by Harkins, G. W. Articles by Varsani, A.