Skip to main content

Comparative analysis of Panicum streak virus and Maize streak virus diversity, recombination patterns and phylogeography

Abstract

Background

Panicum streak virus (PanSV; Family Geminiviridae; Genus Mastrevirus) is a close relative of Maize streak virus (MSV), the most serious viral threat to maize production in Africa. PanSV and MSV have the same leafhopper vector species, largely overlapping natural host ranges and similar geographical distributions across Africa and its associated Indian Ocean Islands. Unlike MSV, however, PanSV has no known economic relevance.

Results

Here we report on 16 new PanSV full genome sequences sampled throughout Africa and use these together with others in public databases to reveal that PanSV and MSV populations in general share very similar patterns of genetic exchange and geographically structured diversity. A potentially important difference between the species, however, is that the movement of MSV strains throughout Africa is apparently less constrained than that of PanSV strains. Interestingly the MSV-A strain which causes maize streak disease is apparently the most mobile of all the PanSV and MSV strains investigated.

Conclusion

We therefore hypothesize that the generally increased mobility of MSV relative to other closely related species such as PanSV, may have been an important evolutionary step in the eventual emergence of MSV-A as a serious agricultural pathogen.

The GenBank accession numbers for the sequences reported in this paper are GQ415386-GQ415401

Background

Panicum streak virus (PanSV) is one of seven known African streak virus species within the Mastrevirus genus of the Geminiviridae. The best studied and most economically relevant species amongst the African streak viruses is Maize streak virus (MSV) which seriously constrains maize production throughout most of sub-Saharan Africa [1]. Like MSV, African streak virus species such as Panicum streak virus (PanSV), Sugarcane streak virus (SSV), Sugarcane streak Reunion virus (SSRV) and Sugarcane streak Egypt virus (SSEV) are transmitted by various leafhopper species in the genus Cicadulina and have geographical ranges that are apparently restricted to Africa and its neighboring islands [1–7].

Whereas African streak virus species such as Eragrostis streak virus (ESV), Saccharum streak virus (SacSV), Urochloa streak virus (USV) and SSEV have been relatively poorly sampled and have therefore only ever been found in individual African countries [2, 8–10], better sampling of MSV and PanSV has indicated that these species occur throughout sub-Saharan Africa [11, 12]. PanSV and MSV display similar degrees of genetic diversity characterized by the existence of multiple discrete strains, many of which have distinctive geographical ranges [11, 12]. Both species also have what appear to be largely overlapping host ranges. Unlike MSV, however, PanSV has no known economic relevance in that it has only ever been found in nature infecting wild grass species in the genera Urochloa, Ehrharta and Panicum [3, 11, 13].

Despite it not having any direct impact on African agriculture, the diversity and phylogeography of PanSV could still provide potentially useful information on other more economically important African streak viruses such as those that cause maize and sugarcane diseases. For example a recent comparative phylogeographic analysis of different MSV strains has indicated that the economically relevant maize adapted MSV-A strain is probably moving around Africa more freely than the closely related but Digitaria adapted MSV-B strain [12]. Comparative analyses of the diversity and phylogeography of different African streak virus species could therefore help identify the characteristics of MSV that facilitated its emergence as an important agricultural pathogen.

It has also been determined that African streak virus species such as PanSV have contributed indirectly to the evolution of MSV through genetic recombination [11, 14]. Recombination is a major force in geminivirus evolution [15, 16] and it appears to have played at least some role in the emergence of a number of serious geminiviral crop diseases [17–22]. At least seven of the eleven currently described MSV strains (including the important MSV-A strain) have apparently come into existence through recombination between two or more other strains [12]. It would be of great interest to determine whether such inter-strain recombination has featured as prominently in the diversification of other African streak viruses such as PanSV.

Here we use 23 full PanSV genome sequences sampled throughout Africa and one of its neighboring islands to show that there generally exist very similar patterns of diversity, recombination and geographical structure within PanSV and MSV populations. Our results indicate, however, that the maize adapted MSV-A strain is possibly unique amongst PanSV and MSV strains in both its total geographical range and the rates at which individual virus variants within the strain are moving across Africa.

Results and discussion

Discovery of five new PanSV strains

Sixteen full mastrevirus genome sequences were cloned and sequenced from Brachiaria deflexa, Panicum maximum, Panicum trichocladium, Urochloa maxima and Ehrharta calycina plants sampled from South Africa, Mozambique, Kenya, Nigeria, the Central African Republic and the Indian Ocean island of Mayotte (Table 1). All shared greater than 80% genome-wide identity with PanSV genomes currently deposited in public databases and were therefore all classified as being PanSV isolates. After confirming that plots of pairwise genetic similarity between all fully sequenced PanSV genomes closely matched those previously determined for MSV (Additional file 1), we used the 93% identity rule that has been used as a MSV strain demarcation criterion [14] to tentatively classify the PanSV isolates. This 93% identity threshold represents a logical, if not natural cutoff for classifying MSV and PanSV strains (Additional file 1) and it indicated that amongst the new sequences there potentially existed five new PanSV-strains (named PanSV-E to -I; Figure 1). It should be noted, however, that this classification scheme relied on the use of similarity measurements that exclude alignment gaps as missing data. Many other geminivirus classification schemes, such as the 75% and 89% thresholds endorsed by the ICTV for respectively demarcating mastrevirus and begomovirus species [23], do not specify how alignment gaps should be handled during similarity measurements. If we had included alignment gaps as a fifth character state - as is often done either accidentally or by design when arguments are made for or against new isolates being considered as new species - the 93% MSV/PanSV strain demarcation threshold would drop to between 90 and 91%.

Figure 1
figure 1

Maximum likelihood phylogenetic trees (both constructed with the GTR+G 4 nucleotide substitution model) indicating possible evolutionary relationships between 23 PanSV isolates. (a) Tree constructed using complete genome sequences. Virus names take the form "Strain [country-region-year of isolation]" (b) A tree constructed using recombination-free portions of the genome indicated beneath the genome map and recombination mosaic cartoons in c. (c) Linearised genome cartoons depicting unique recombinant mosaics detected amongst the PanSV sequences. Colours represent as best as possible the origins of different genome regions. Letters below the depicted recombination events correspond to detailed descriptions of each of the events given in Additional file 7. For labels on the genome map: mp = movement protein gene, cp = coat protein gene, rep = replication associated protein gene, repA = RepA gene. Whereas branches marked with filled and open circles were supported in >90% and 70-89% of bootstrap replicates, respectively, branches with <50% bootstrap support have been collapsed. The tree was rooted on the Sugarcane streak Reunion virus isolate SSRV-Bas (not shown).

Table 1 Accession numbers and sampling coordinates of various mastrevirus isolates used to detect recombination in PanSV.

Although the newly described PanSV-I strain represented the most divergent group of PanSV isolates yet discovered, we found no major genomic features that could distinguish this or any of the other newly described PanSV strains from those already represented in public sequence databases (see additional files 2, 3, 4, 5 and 6 for annotated genome and probable expressed protein maps).

Recombination between PanSV strains

It has been previously determined that recombination has featured prominently in the evolution of MSV strains [12] and that it may have also contributed substantially to the diversification of PanSV [11]. We therefore analysed the PanSV sequences for evidence of inter-species and inter-strain recombination events using a battery of recombination detection and analysis methods implemented in the program RDP3 [24]. We identified clear evidence of three inter-species (labeled a, b, d, e, f, and h in Figure 1 and Additional File 7) and six inter-strain recombination events (labeled a, b, d, e, f, and h in Figure 1) within the PanSV sequences.

The pattern of recombination we observed in PanSV is very similar to that which has been described for MSV [12]. For both species most detectable recombination events have involved intra-species sequence exchanges. The few inter-species recombination events that have been detected in both species have also all involved the exchange of small (<200 nt) tracts of sequence.

Another similarity between the two species is that many of the described strains have apparently arisen through inter-strain recombination events. For MSV all currently sampled isolates of seven of the eleven described strains (MSV-A, -F, -H, -J, K, C and D) share evidence of ancestral inter-strain recombination events that involved exchanges of genome fragments >30% of the full genome [12]. Likewise, exchanges of >30% genome size fragments are evident in all sampled representatives of five of the nine PanSV strains (PanSV-B, -C, -H, F and G).

The patterns of recombination seen in PanSV and MSV, where inter-species recombination events generally involve exchanges of only small genomic fragments (<10% of the full genome length), is quite different to that seen amongst related whitefly transmitted geminiviruses in the genus Begomovirus [15, 16]. In these viruses inter-species recombination is very common and often involves exchanges of large (>30% of the full genome length) genome fragments. This difference is due, at least in part, to differences between the species classification criteria used for mastreviruses and begomoviruses. Whereas the main begomovirus species demarcation criterion is that DNA-A or DNA-A-like sequences (begomoviruses often have two component genomes where the DNA-A component of such genomes is largely homologous to mastrevirus genomes) sharing <89% identity belong to different species, the analogous mastrevirus species demarcation threshold is 75%. If the begomovirus classification scheme were applied to PanSV and MSV then, many of the inter-strain recombination events detectable in these species would be "upgraded" to inter-species recombination events.

It is still noteworthy, however, that detectable recombination events between more distantly related PanSV and MSV genomes have been less frequent and have tended to involve smaller sequence exchanges than recombination events between more closely related genomes. It is possible that the observed ratios of intra:inter species recombination events in PanSV and MSV might be partially attributable to mixed infections involving different mastrevirus species being rarer than mixed infections involving different strains of the same species. Although many of the different African streak virus species share hosts such as Urochloa and Eragrostis species, there are probably greater host-range differences between viruses in different species than there are between viruses within the same species. Such differences should surely influence the relative frequencies of mixed species and mixed strain infections and should therefore also influence the relative rates of inter-species and intra-species recombination events.

The most striking difference between the inter- and intra-species recombination events in these viruses is, however, not their relative frequencies, but rather the relative amounts of sequence that have been exchanged in these different recombination event categories. This pattern of recombination in fact conforms very well with the hypothesis that a major determinant of recombinant fitness is how well foreign DNA fragments function within the context of genomic backgrounds that they did not co-evolve within [25–29]. Functional nucleotide sequences tend to work best within genomes that are similar to the ones in which they evolved [25, 27, 30]. The probable reason for this is simply that the interaction networks that define the functionality of a particular nucleotide sequence within any given genomic context could potentially be disrupted if that sequence were placed into a genome where it was forced to interact with nucleotide sequences different from those it co-evolved with. As the relatedness between prospective parental sequences drops so too should the proportion of their genomes that could be exchanged without disrupting the delicate intra-genomic interactions required for optimal fitness [27, 31]. The net effect of this process should be that amongst (presumably high fitness) genomes sampled from nature, one should tend to observe larger sequence exchanges between more closely related genomes than are detectable between less closely related ones. This is the exact pattern of recombination seen in both PanSV and MSV, suggesting that rather than inter-species recombination events being uncommon due entirely to different mastrevirus species only rarely infecting the same hosts, they are uncommon because of genetic constraints on the relative viability of inter-species recombinants.

PanSV and MSV phylogenies display similar patterns of geographical structure

It has been previously demonstrated that there are strong signals of geographical structure within the phylogenetic trees of both the maize adapted MSV strain, MSV-A [14, 22], and the grass adapted MSV strain, MSV-B [12]. These two strains differ, however, in the degree to which viruses have been moving across Africa [12]. Whereas no MSV-B isolates have ever been detected in West Africa, there have apparently been no movements of MSV-B isolates between East Africa, southern Africa and the Indian Ocean island of La Reunion since the initial spread of this strain to these three locations. Conversely, in the time since MSV-A first spread throughout the continent there have apparently been multiple instances where these viruses have moved between the major regions of Africa [12, 22].

We sought to determine whether similar phylogeographic patterns exist amongst the currently sampled PanSV sequences. Taking note of the locations from which sequences were sampled, we compared the PanSV and MSV phylogenies (Figure 2) and noted some striking similarities between them with respect to the geographical ranges of the various distinct strain groupings represented.

Figure 2
figure 2

Maximum likelihood phylogenetic trees (best fit nucleotide substitution models = GTR+G 4 for PanSV tree and GTR+I+G 4 for the MSV tree) depicting the sampling locations of 23 PanSV (tree a) and 181 MSV isolates (trees b and c). Branches are colored according to sampling locations (blue = southern African lineages, orange = East African lineages, green = West African lineages, purple = Indian Ocean island lineages). Wherever under a maximum parsimony criterion, it is <60% certain that ancestral sequences represented by tree nodes are from one of the four regions, the branches basal to that node have been left uncoloured. Whereas branches marked with filled and open circles were supported in >90% and 70-89% of bootstrap replicates, respectively, branches with <50% bootstrap support have been collapsed. The arrow on tree c indicates a clade of MSV-A sequences from West Africa nested within a clade of MSV-A sequences from East Africa indicating an instance of recent east to west movement of MSV-A isolates. All the trees were rooted on the Sugarcane streak Reunion virus isolate SSRV-Bas (not shown).

Besides both MSV and PanSV strains clearly grouping according to their sampling locations, it is evident that PanSV (Figure 2a) and grass adapted MSV strains (Figure 2b) from East Africa, Southern Africa and the Indian Ocean islands are generally more closely related to one another than they are to viruses from West Africa. Taken together with the MSV data, the PanSV sample therefore provides additional evidence that, in general, African streak viruses may move more freely between East Africa, southern Africa and the Indian Ocean islands than they do between these regions and West Africa [12]. This pattern is notably different from that seen for both the maize adapted MSV-A strain and whitefly transmitted geminiviruses in the genus Begomovirus. Whereas island begomovirus populations display strong evidence of extensive isolation from mainland lineages [32, 33], there is good phylogeographic evidence, particularly for cassava infecting begomoviruses, of lineages moving between the major regions of the continent [34]. Cassava infecting begomoviruses might, however, represent a special case in that cassava is propagated from cuttings and its viruses might therefore be moved more extensively by humans than viruses infecting seed propagated hosts.

Unlike with the wild-grass infecting MSV and PanSV strains, the maize adapted MSV-A strain (Figure 2c) has apparently moved quite extensively throughout the continent with the Indian Ocean islands being relatively more isolated than West Africa [12, 35]. There is in fact clear evidence of at least one fairly recent movement of a MSV-A lineage from East Africa to West Africa (see arrow indicating the green clade nested within the orange clade in Figure 2c). Similar east to west movements across Africa have been detected in various other vector-born plant viruses including whitefly transmitted cassava infecting geminivirus species [34, 36] and the beetle transmitted sobemovirus species, Rice yellow mottle virus [37, 38]. It remains to be determined, however, whether movement of MSV-A and perhaps these other viruses too is natural or whether it is facilitated by human trafficking of infected plant material/viruliferous vectors [12, 39]. It is also currently unknown whether PanSV, the grass adapted MSV strains and MSV-A are adapted to transmission by either different Cicadulina species or different biotypes within these species. MSV-A is transmitted with varying efficiencies by different Cicadulina species [40], and it remains a strong possibility that differences in the geographical distribution and migration routes of different preferred vector species might also account for differences in the movement patterns of these virus groups across East and West Africa.

The final subtle difference between the grass adapted MSV-strains and the PanSV dataset are the genetic distances between viruses found in different regions. Since the same demarcation threshold was used in both the MSV and PanSV strain classifications it is perhaps interesting that in no case was any PanSV strain detected in more than one of the surveyed regions. Members of each of the MSV-B -C, -E,-F and -K strains have been isolated across multiple geographical regions (three for MSV-B and two each for the rest). Assuming that the PanSV and grass adapted MSV strains are evolving at approximately the same rate, this observation indicates that over the time-scales represented by these phylogenies, MSV moves more frequently than PanSV between the regions examined. Without comparative analysis of PanSV and MSV substitution rates, however, it cannot be discounted that, rather than moving at different rates, the two species are simply evolving at different rates. It is also possible that with better sampling the isolates of different PanSV strains will, as is the case for MSV, be found in multiple different regions of Africa.

With these reservations noted, it is nevertheless interesting that just as MSV-A seems to be moving across Africa with less restraint than grass adapted MSV strains, the grass adapted MSV strains are in turn apparently moving more freely across the continent than PanSV strains. It is therefore possible that the evolution of epidemiological traits enabling MSV to move more rapidly than PanSV across Africa was important for the eventual evolution of still faster rates of MSV-A movement.

Conclusion

Among 16 new PanSV isolates sampled across Africa and the Indian Ocean island of Mayotte we have potentially discovered five new PanSV strains. Together with other currently sampled PanSV genome sequences these new PanSV isolates indicate that there exist striking similarities between PanSV and MSV with respect to both detectable recombination patterns and degrees of geography-associated population structure. Although similarities between PanSV and MSV are perhaps unsurprising considering that these viruses share common leafhopper vector species and partially overlapping host-ranges, it remains interesting that both MSV strains in general, and the MSV-A strain in particular, seem to be less constrained in their movements across Africa than PanSV strains.

Methods

Virus Isolates

Sixteen grasses presenting with mild streak symptoms characteristic of African streak virus infections were sampled from various locations in Africa and the Indian ocean island of Mayotte: one Brachiaria deflexa sample from Central African Republic; four Panicum maximum samples from Mayotte; one Urochloa maxima and one Panicum maximum sample from Nigeria; four Panicum maximum, two Brachiaria deflexa and one Panicum trichocladium samples from Kenya; one Ehrharta calycina sample from South Africa and one Panicum maximum sample from Mozambique (locations and names are provided in Table 1). The Nigerian isolate, PanSV-H [NG-Jic15-PNP-1997], and the Kenyan isolates, PanSV-E [KE-Jic10-PKPM-1997] and PanSV-I [KE-Jic13-PKPB-1997] (respectively referred to as P(N)P, P(K)P-M and P(K)P-B in [5, 41]), were sampled in ~1987 but were maintained for approximately ten years within Panicum maximum under glasshouse conditions at the John Innes Centre in Norwich prior to leaf tissues being harvested and frozen. Full length PanSV genomes were amplified from leaf tissues using rolling circle amplification, cloned and sequenced using methods described previously [42–44]. Briefly, total DNA was either extracted from leaf tissues using Extract-n-Amp™ Plant PCR Kit (Sigma-Aldrich Corporation, USA) or using a Qiagen Plant miniprep DNA kit (Qiagen, Germany) and circular DNA molecules were amplified using φ 29 DNA polymerase (TempliPhi™, GE Healthcare). The amplified concatamers were digested with Bam HI, Sal I or Xho I restriction enzymes to release ~2.7 kb PanSV genomes which were subsequently ligated to similarly linearised pGEM3 Zf(+) (from Promega Biotech). The cloned PanSV genomes were sequenced by Macrogen Inc (Korea) using primer walking. Sequences were assembled and edited using DNAMAN (version 5.2.9; Lynnon Biosoft) and MEGA (version 4)[45].

Diversity Analysis

34 African streak virus full genome sequences, including all those available in GenBank for PanSV [3, 11], and representative selections of MSV [12], USV [10], ESV [8], Sugarcane streak virus [4, 8], Sugarcane streak Egypt virus [2], and Sugarcane streak Reunion virus [2, 8], were obtained from GenBank. These were aligned together with the 16 new PanSV sequences using POA (vesion 2) [46] and edited by eye using MEGA. For purposes of assigning PanSV sequences to different strain groupings using the 93% rule of Martin et al [14], MEGA was also used to calculate the pair-wise differences between aligned PanSV genomes using p-distances with pair-wise deletion of gaps (as opposed to scoring gaps as a fifth nucleotide state). Alignments used in earlier phylogeographic analyses described in [12] and [22] were merged (with duplicate sequences being discarded) and realigned with MEGA.

Recombination and phylogenetic analysis

Maximum likelihood phylogenetic trees were constructed using PHYML (version 1)[47] with automated best-fit model selection under the Akaike information criterion as described in [48].

Discreet recombination events were detected using the RDP [49], GENECONV [16], BOOTSCAN [50], MAXCHI [51], CHIMAERA [52], SISCAN [53], and 3SEQ [54] methods implemented in the program RDP3 (version 3.32; available from http://darwin.uvigo.es/rdp/rdp.html)[24]. Only potential recombination signals detected by at least three of the seven applied recombination detection methods, coupled with phylogenetic evidence of recombination were considered significant evidence of the signals representing genuine recombination events. Parental and recombinant sequences were identified from the sets of sequences used to detect recombination events as outlined in [55] and [56]. Recombination breakpoint positions and recombinant/parental designations were manually checked and adjusted where necessary using the extensive phylogenetic and recombination signal analysis features implemented in RDP3.

Abbreviations

ESV:

Eragrostis streak virus

cp :

coat protein gene

LIR:

long intergenic region

mp :

movement protein gene

MSV:

Maize streak virus

PanSV:

Panicum streak virus

rep :

replication associate protein gene

repA :

RepA gene

SIR:

short intergenic region

SSEV:

Sugarcane streak Egypt virus

SSRV:

Sugarcane streak Reunion virus

SSV:

Sugarcane streak virus

USV:

Urochloa streak virus.sacSV.

References

  1. Shepherd DN, Martin DP, Walt E, Dent K, Varsani A, Rybicki EP: Maize streak virus: an old and complex 'emerging' pathogen. Mol Plant Path 2009, in press.

    Google Scholar 

  2. Bigarré L, Salah M, Granier M, Frutos R, Thouvenel JC, Peterschmitt M: Nucleotide sequence evidence for three distinct sugarcane streak mastreviruses. Arch Virol 1999, 144: 2331-2344. 10.1007/s007050050647

    Article  PubMed  Google Scholar 

  3. Briddon RW, Lunness P, Chamberlin LC, Pinner MS, Brundish H, Markham PG: The nucleotide sequence of an infectious insect-transmissible clone of the geminivirus Panicum streak virus. J Gen Virol 1992, 73: 1041-1047. 10.1099/0022-1317-73-5-1041

    Article  CAS  PubMed  Google Scholar 

  4. Hughes FL, Rybicki EP, Kirby R: Complete nucleotide sequence of sugarcane streak Monogeminivirus. Arch Virol 1993, 132: 171-182. 10.1007/BF01309851

    Article  CAS  PubMed  Google Scholar 

  5. Pinner MS, Markham PG, Markham RH, Dekker EL: Characterization of maize streak virus:description of strains; symptoms. Plant Path 1988, 37: 74-87. 10.1111/j.1365-3059.1988.tb02198.x

    Article  Google Scholar 

  6. Peterschmitt M, Reynaud B, Sommermeyer G, Baudin P: Characterization of maize streak virus isolates using monoclonal and polyclonal antibodies and by transmission to a few hosts. Plant Dis 1991, 75: 27-32.

    Article  Google Scholar 

  7. Storey HH, McClean APD: The transmission of streak disease between maize, sugarcane and wild grasses. Ann Appl Biol 1930, 17: 691-719. 10.1111/j.1744-7348.1930.tb07240.x

    Article  Google Scholar 

  8. Shepherd DN, Varsani A, Windram O, Lefeuvre P, Monjane AL, Owor B, Martin DP: Novel Sugarcane streak virus and Sugarcane streak Reunion virus mastrevirus isolates from Southern Africa and La Reunion. Arch Virol 2008, 153: 605-609. 10.1007/s00705-007-0016-3

    Article  CAS  PubMed  Google Scholar 

  9. Lawry R, Martin DP, Shepherd DN, van Antwerpen T, Varsani A: A novel sugarcane-infecting mastrevirus species from South Africa. Arch Virol 2009, 154: 1699-1703. 10.1007/s00705-009-0490-x

    Article  CAS  PubMed  Google Scholar 

  10. Oluwafemi S, Varsani A, Monjane AL, Shepherd DN, Owor BE, Rybicki EP, Martin DP: A new African streak virus species from Nigeria. Arch Virol 2008, 153: 1407-1410. 10.1007/s00705-008-0123-9

    Article  CAS  PubMed  Google Scholar 

  11. Varsani A, Oluwafemi S, Shepherd DN, Monjane AL, Owor B, Windram O, Rybicki EP, Lefeuvre P, Martin DP: Panicum streak virus diversity is similar to that observed for Maize streak virus . Arch Virol 2008, 153: 601-604. 10.1007/s00705-007-0020-7

    Article  CAS  PubMed  Google Scholar 

  12. 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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Schnippenkoetter WH, Martin DP, Hughes F, Fyvie M, Willment JA, James D, von Wechmar B, Rybicki EP: The biological and genomic characterisation of three mastreviruses. Arch Virol 2001, 146: 1075-1088. 10.1007/s007050170107

    Article  CAS  PubMed  Google Scholar 

  14. Martin DP, Willment JA, Billharz R, Velders R, Odhiambo B, Njuguna J, James D, Rybicki EP: Sequence diversity and virulence in Zea mays of Maize streak virus isolates. Virology 2001, 288: 247-255. 10.1006/viro.2001.1075

    Article  CAS  PubMed  Google Scholar 

  15. Lefeuvre P, Lett JM, Varsani A, Martin DP: Widely conserved recombination patterns among single-stranded DNA viruses. J Virol 2009, 83: 2697-2707. 10.1128/JVI.02152-08

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Padidam M, Sawyer S, Fauquet CM: Possible emergence of new geminiviruses by frequent recombination. Virology 1999, 265: 218-225. 10.1006/viro.1999.0056

    Article  CAS  PubMed  Google Scholar 

  17. Fondong VN, Pita JS, Rey ME, de Kochko A, Beachy RN, Fauquet CM: Evidence of synergism between African cassava mosaic virus and a new double-recombinant geminivirus infecting cassava in Cameroon. J Gen Virol 2000, 81: 287-297.

    Article  CAS  PubMed  Google Scholar 

  18. Pita JS, Fondong VN, Sangaré A, Otim-Nape GW, Ogwal S, Fauquet CM: Recombination, pseudorecombination and synergism of geminiviruses are determinant keys to the epidemic of severe cassava mosaic disease in Uganda. J Gen Virol 2001, 82: 655-665.

    Article  CAS  PubMed  Google Scholar 

  19. Monci F, Sánchez-Campos S, Navas-Castillo J, Moriones E: 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 2002, 303: 317-326. 10.1006/viro.2002.1633

    Article  CAS  PubMed  Google Scholar 

  20. García-Andrés S, Tomas DM, Sanchez-Campos S, Navas-Castillo J, Moriones E: Frequent occurrence of recombinants in mixed infections of tomato yellow leaf curl disease-associated begomoviruses. Virology 2007, 365: 210-219. 10.1016/j.virol.2007.03.045

    Article  PubMed  Google Scholar 

  21. García-Andrés S, Accotto GP, Navas-Castillo J, Moriones E: 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 2007, 359: 302-312. 10.1016/j.virol.2006.09.030

    Article  PubMed  Google Scholar 

  22. Harkins GW, Martin DP, Duffy S, Monjane AL, Shepherd DN, Windram OP, Owor BE, Donaldson L, van Antwerpen T, Sayed RA, Flett B, Ramusi M, Rybicki EP, Peterschmitt M, Varsani A: Dating the origins of the maize-adapted strain of maize streak virus, MSV-A. J Gen Virol 2009, 90: 3066-3074. 10.1099/vir.0.015537-0

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Stanley J, Bisaro DM, Briddon RW, Brown JK, Fauquet CM, Harrison BD, Rybicki EP, Stenger DC: Geminiviridae. Edited by: Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA. Virus Taxonomy (VIIIth Report of the ICTV). Elsevier/Academic Press, London; 2005:301-306.

    Google Scholar 

  24. Martin DP, Williamson C, Posada D: RDP2: recombination detection and analysis from sequence alignments. Bioinformatics 2005, 21: 260-262. 10.1093/bioinformatics/bth490

    Article  CAS  PubMed  Google Scholar 

  25. Escriu F, Fraile A, García-Arenal F: Constraints to genetic exchange support gene coadaptation in a tripartite RNA virus. PLoS Pathog 2007, 3: e8. 10.1371/journal.ppat.0030008

    Article  PubMed Central  PubMed  Google Scholar 

  26. Lefeuvre P, Lett JM, Reynaud B, Martin DP: Avoidance of protein fold disruption in natural virus recombinants. PLoS Pathog 2007, 3: e181. 10.1371/journal.ppat.0030181

    Article  PubMed Central  PubMed  Google Scholar 

  27. Martin DP, Walt E, Posada D, Rybicki EP: The evolutionary value of recombination is constrained by genome modularity. PLoS Genet 2005, 1: e51. 10.1371/journal.pgen.0010051

    Article  PubMed Central  PubMed  Google Scholar 

  28. Moreno IM, Malpica JM, Diaz-Pendon JA, Moriones E, Fraile A, García-Arenal F: Variability and genetic structure of the population of watermelon mosaic virus infecting melon in Spain. Virology 2004, 318: 451-460. 10.1016/j.virol.2003.10.002

    Article  CAS  PubMed  Google Scholar 

  29. Martin DP, Rybicki EP: Investigation of Maize streak virus pathogenicity determinants using chimaeric genomes. Virology 2002, 300: 180-188. 10.1006/viro.2002.1458

    Article  CAS  PubMed  Google Scholar 

  30. Walt E, Palmer K, Martin DP, Rybicki EP: Viable chimaeric viruses confirm the biological importance of sequence specific maize streak virus movement protein and coat protein interactions. Virol J 2008, 5: 61. 10.1186/1743-422X-5-61

    Article  PubMed Central  PubMed  Google Scholar 

  31. Meyer MM, Silberg JJ, Voigt CA, Endelman JB, Mayo SL, Wang ZG, Arnold FH: Library analysis of SCHEMA-guided protein recombination. Protein Sci 2003, 12: 1686-1693. 10.1110/ps.0306603

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Lefeuvre P, Martin DP, Hoareau M, Naze F, Delatte H, Thierry M, Varsani A, Becker N, Reynaud B, Lett JM: Begomovirus 'melting pot' in the south-west Indian Ocean islands: molecular diversity and evolution through recombination. J Gen Virol 2007, 88: 3458-3468. 10.1099/vir.0.83252-0

    Article  CAS  PubMed  Google Scholar 

  33. Delatte H, Martin DP, Naze F, Goldbach R, Reynaud B, Peterschmitt M, Lett JM: South West Indian Ocean islands tomato begomovirus populations represent a new major monopartite begomovirus group. J Gen Virol 2005, 86: 1533-1542. 10.1099/vir.0.80805-0

    Article  CAS  PubMed  Google Scholar 

  34. Ndunguru J, Legg JP, Aveling TA, Thompson G, Fauquet CM: Molecular biodiversity of cassava begomoviruses in Tanzania: evolution of cassava geminiviruses in Africa and evidence for East Africa being a center of diversity of cassava geminiviruses. Virol J 2005, 2: 21. 10.1186/1743-422X-2-21

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Peterschmitt M, Granier M, Frutos R, Reynaud B: Infectivity and complete nucleotide sequence of the genome of a genetically distinct strain of maize streak virus from Réunion Island. Arch Virol 1996, 141: 1637-1650. 10.1007/BF01718288

    Article  CAS  PubMed  Google Scholar 

  36. Legg JP, Fauquet CM: Cassava mosaic geminiviruses in Africa. Plant Mol Biol 2004, 56: 585-599. 10.1007/s11103-004-1651-7

    Article  CAS  PubMed  Google Scholar 

  37. Traoré O, Pinel-Galzi A, Sorho F, Sarra S, Rakotomalala M, Sangu E, Kanyeka Z, Séré Y, Konaté G, Fargette D: A reassessment of the epidemiology of Rice yellow mottle virus following recent advances in field and molecular studies. Virus Res 2009, 141: 258-267. 10.1016/j.virusres.2009.01.011

    Article  PubMed  Google Scholar 

  38. Traore O, Sorho F, Pinel A, Abubakar Z, Banwo O, Maley J, Hebrard E, Winter S, Sere Y, Konate G, Fargette D: Processes of diversification and dispersion of rice yellow mottle virus inferred from large-scale and high-resolution phylogeographical studies. Mol Ecol 2005, 14: 2097-2110. 10.1111/j.1365-294X.2005.02578.x

    Article  CAS  PubMed  Google Scholar 

  39. 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-260. 10.1146/annurev.phyto.44.120705.104644

    Article  CAS  PubMed  Google Scholar 

  40. Oluwafemi S, Jackai LE, Alegbejo MD: Comparison of transmission abilities of four Cicadulina species vectors of maize streak virus from Nigeria. Ent Exp Et App 2007, 124: 235-239. 10.1111/j.1570-7458.2007.00577.x

    Article  Google Scholar 

  41. Pinner MS, Markham PG: Serotyping and strain identification of maize streak virus isolates. J Gen Virol 1990, 71: 1635-1640. 10.1099/0022-1317-71-8-1635

    Article  PubMed  Google Scholar 

  42. Inoue-Nagata AK, Albuquerque LC, Rocha WB, Nagata T: A simple method for cloning the complete begomovirus genome using the bacteriophage phi29 DNA polymerase. J Virol Methods 2004, 116: 209-211. 10.1016/j.jviromet.2003.11.015

    Article  CAS  PubMed  Google Scholar 

  43. Owor BE, Shepherd DN, Taylor NJ, Edema R, Monjane AL, Thomson JA, Martin DP, Varsani A: Successful application of FTA Classic Card technology and use of bacteriophage phi29 DNA polymerase for large-scale field sampling and cloning of complete maize streak virus genomes. J Virol Methods 2007, 140: 100-105. 10.1016/j.jviromet.2006.11.004

    Article  CAS  PubMed  Google Scholar 

  44. Shepherd DN, Martin DP, Lefeuvre P, Monjane AL, Owor BE, Rybicki EP, Varsani A: A protocol for the rapid isolation of full geminivirus genomes from dried plant tissue. J Virol Methods 2008, 149: 97-102. 10.1016/j.jviromet.2007.12.014

    Article  CAS  PubMed  Google Scholar 

  45. 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

    Article  CAS  PubMed  Google Scholar 

  46. Grasso C, Lee C: Combining partial order alignment and progressive multiple sequence alignment increases alignment speed and scalability to very large alignment problems. Bioinformatics 2004, 20: 1546-1556. 10.1093/bioinformatics/bth126

    Article  CAS  PubMed  Google Scholar 

  47. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 2003, 52: 696-704. 10.1080/10635150390235520

    Article  PubMed  Google Scholar 

  48. Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol 2008, 25: 1253-1256. 10.1093/molbev/msn083

    Article  CAS  PubMed  Google Scholar 

  49. Martin D, Rybicki E: RDP: detection of recombination amongst aligned sequences. Bioinformatics 2000, 16: 562-563. 10.1093/bioinformatics/16.6.562

    Article  CAS  PubMed  Google Scholar 

  50. Martin DP, Posada D, Crandall KA, Williamson C: A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS and Human Retroviruses 2005, 21: 98-102. 10.1089/aid.2005.21.98

    Article  CAS  Google Scholar 

  51. Maynard Smith J: Analyzing the mosaic structure of genes. J Mol Evol 1992, 34: 126-129.

    Google Scholar 

  52. Posada D, Crandall KA: The effect of recombination on the accuracy of phylogeny estimation. J Mol Evol 2002, 54: 396-402.

    Article  CAS  PubMed  Google Scholar 

  53. 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

    Article  CAS  PubMed  Google Scholar 

  54. Boni MF, Posada D, Feldman MW: An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics 2007, 176: 1035-1047. 10.1534/genetics.106.068874

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Heath L, Walt E, Varsani A, Martin DP: Recombination patterns in aphthoviruses mirror those found in other picornaviruses. J Virol 2006, 80: 11827-11832. 10.1128/JVI.01100-06

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Lemey P, Lott M, Martin DP, Moulton V: Identifying recombinants in human and primate immunodeficiency virus sequence alignments using quartet scanning. BMC Bioinformatics 2009, 10: 126. 10.1186/1471-2105-10-126

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

This work was funded by the South African National Research Foundation. AV was supported by the Carnegie Corporation of New York. DPM was supported by the Wellcome Trust. PL and JML were supported by the Conseil Régional de La Réunion, CIRAD, European Union (FEDER) and GIS Centre de recherche et de veille sanitaire sur les maladies émergentes dans l'océan Indien (N°PRAO/AIRD/CRVOI/08/03). RWB is supported by the Higher Education Commission, Government of Pakistan, under the "Foreign Faculty Hiring Program".

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Darren P Martin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AV, DPM, ALM, JML, SO, JZ, EKK, DPL, NM. JM, SS PGM and RWB collected isolates. AV, ALM, and LD cloned and sequenced the viruses. AV, PL and DPM conceived the project, AV and DPM analysed the data. AV and DPM prepared the manuscript and AV, PL, JML and DPM secured funding for the project's execution. EPR, PL and RWB provided ideas and comments during manuscript preparation. All authors other than PGM read and approved the final manuscript.

Electronic supplementary material

12985_2009_708_MOESM1_ESM.ppt

Additional file 1: Graph rationalizing the use with PanSV of the same 93% sequence identity strain demarcation threshold used for MSV. Graph rationalizing the use with PanSV of the same 93% sequence identity strain demarcation threshold used for MSV. The red and blue splines respectively plots the frequencies of pairwise sequence identities shared amongst 23 PanSV (253 pairwise distances) and 99 MSV isolates (corresponding to the MSV dataset used in Varsani et al [2008] and accounting for 4851 pairwise distances) at a resolution of 1% identity. Identity values were calculated with pairwise exclusion of alignment gaps (as opposed to counting gaps as a fifth state as is often currently done, either by accident or design, by many geminivirologists). (PPT 68 KB)

12985_2009_708_MOESM2_ESM.doc

Additional file 2: Full genome sequence alignments of 23 PanSV isolates. Annotated full genome sequence alignments of 23 PanSV isolates. Sequences either known or believed to have some role in mastrevirus replication and transcription are marked together with a corresponding label on the nucleotide sequence alignments. To highlight differences between the sequences, wherever nucleotides in a particular alignment column are identical to that of PanSV-A [ZM-Nya-g180-2007] they are replaced with a "-" character. In columns where they differ from PanSV-A [ZM-Nya-g180-2007] they are shown in lower case. "." characters indicate where gaps were inserted to align the sequences. [1] Stenger, et al., 1991. Proc. Natl. Acad. Sci. USA 88:8029; [2] Sunter, et al. 1985. Nucl. Acids Res. 13:4645; [3] Argüello-Astorga et al. 1994. Virology 203:90; [4] Suárez-López et al. 1995. Virology. 208:303; [5] Morris-Krsinich et al. 1984. Nucleic Acids Res. 13:7237; [6] Boulton et al. 1989. J. Gen. Virol. 70:2309; [7] Wright et al.1997 Plant J. 12:1285; [8] Donson et al. 1984. EMBO J. 3:3069; [9] Dekker et al.1991. Nucl. Acids Res. 19:4075; [10] Fenoll et al. 1990. Plant. Mol. Biol. 15:865. (DOC 238 KB)

12985_2009_708_MOESM3_ESM.doc

Additional file 3: Annotated predicted movement protein amino acid sequence alignments. Annotated predicted movement protein amino acid sequence alignments of 23 PanSV isolates. The hydrophobic, potentially membrane spanning internal domain of the sequences is highlighted. [1] Wright et al.1997. Plant J. 12:1285. (DOC 24 KB)

12985_2009_708_MOESM4_ESM.doc

Additional file 4: Annotated predicted coat protein amino acid sequence alignments. Annotated predicted coat protein amino acid sequence alignments of 23 PanSV isolates. The potential nuclear localization signal and DNA binding domains (inferred by analogy with those determined for MSV) are highlighted on the sequence. [1] Liu et al. 1999. Mol. Plant Microbe Interact. 12:894; [2] Liu et al. 1997. J. Gen. Virol. 78:1265. (DOC 29 KB)

12985_2009_708_MOESM5_ESM.doc

Additional file 5: Annotated predicted replication-associated protein amino acid sequence alignments. Annotated predicted replication-associated protein amino acid sequence alignments of 23 PanSV isolates. Potential rolling-circle replication motifs and interaction domains inferred by analogy with MSV and Wheat dwarf virus are highlighted. [1] Koonin & Ilyina. 1992. J Gen Virol, 73:2763; [2] Horvath et al. 1998. Plant Mol. Biol. 38:699; [3] Xie et al. 1995. EMBO J. 14:4073; [4] Gorbalenya & Koonin. 1989. Nucl. Acids Res. 17:8413. (DOC 34 KB)

12985_2009_708_MOESM6_ESM.doc

Additional file 6: Annotated predicted RepA amino acid sequence alignments. Annotated predicted RepA amino acid sequence alignments of 23 PanSV isolates. Potential rolling circle replication motifs and interaction domains inferred by analogy with MSV and Wheat dwarf virus are highlighted. [1] Koonin & Ilyina. 1992. J Gen Virol, 73:2763; [2] Horvath et al. 1998. Plant Mol. Biol. 38:699; [3] Xie et al. 1995. EMBO J. 14:4073; [4] Xie et al. 1999. Plant Mol. Biol. 39:647. (DOC 32 KB)

12985_2009_708_MOESM7_ESM.xls

Additional file 7: Recombination in PanSV. Details of the PanSV recombination events detected in this study including approximate breakpoint positions, parental-like sequences, and p-values for various recombination detection tests. (XLS 20 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Varsani, A., Monjane, A.L., Donaldson, L. et al. Comparative analysis of Panicum streak virus and Maize streak virus diversity, recombination patterns and phylogeography. Virol J 6, 194 (2009). https://doi.org/10.1186/1743-422X-6-194

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1743-422X-6-194

Keywords