Natural recombination between Tomato yellow leaf curl virus-Is and Tomato leaf curl virus

J. Navas-Castillo1, S. Sánchez-Campos1, E. Noris2, D. Louro3, G. P. Accotto2 and E. Moriones1

Estación Experimental ‘La Mayora’, Consejo Superior de Investigaciones Científicas, 29750 Algarrobo-Costa, Málaga, Spain1
Istituto di Fitovirologia Applicata, CNR, Strada delle Cacce 73, 10135 Torino, Italy2
Direcçâo Geral de Protecçâo das Culturas, Quinta do Marquês, 2780 Oeiras, Portugal3

Author for correspondence: Enrique Moriones. Fax +34 952552677. e-mail moriones{at}eelm.csic.es


   Abstract
Top
Abstract
Main text
References
 
The complete genome sequences (2791 and 2793 nt) of isolates of Tomato yellow leaf curl virus-Is (TYLCV-Is) from Spain (SP72/97) and Portugal (Port2/95) were determined. These isolates are closely related to TYLCV-Is isolates reported in Japan (Japan-A and Japan-S) and Israel (Israel/Mild). Comparison of all sequenced isolates of TYLCV-Is showed that part of the genome comprising the intergenic region and the 5'-end of the rep gene of the Iran and Israel isolates was not closely related to that of other isolates. Phylogenetic analyses suggest that the Israel and Iran isolates may have chimeric genomes that have arisen by recombination between TYLCV-Is-like and tomato leaf curl virus (ToLCV)-like ancestors. The TYLCV-Is donors of the Iran and the Israel genomes were closely related to each other and to other known TYLCV-Is isolates. However, the ToLCV donors differed from each other, although both were related to ToLCV isolates from India (Bangalore-2 and Bangalore-4).


   Main text
Top
Abstract
Main text
References
 
Tomato yellow leaf curl is one of the most devastating virus diseases of tomato crops in tropical and subtropical regions (Czosnek & Laterrot, 1997 ). The causal agents are a group of virus species of the genus Begomovirus, family Geminiviridae, all of them named Tomato yellow leaf curl virus (TYLCV), that are transmitted by the whitefly Bemisia tabaci Gen. (Briddon & Markham, 1995 ). The TYLCV species have a circular, single-stranded DNA genome that, except in one case in which two components were isolated (Rochester et al., 1994 ), is composed of a single molecule. Monopartite TYLCVs have a genome (about 2·8 kb) that contains six partially overlapping open reading frames (ORFs), two (V1 and V2) on the virion-sense strand and four (C1 to C4) on the complementary-sense strand, separated by an intergenic region (IR) of about 300 nt. Biological, molecular, epidemiological and control aspects of the TYLCV complex have been reviewed in Picó et al. (1996 ) and, more recently, in Moriones & Navas-Castillo (2000 ).

TYLCV has been known in the Iberian Peninsula since 1992, when the species TYLCV-Sar was reported as the causal agent of yellow leaf curl epidemics in tomatoes in southern Spain (Noris et al., 1994 ). The TYLCV-Is species was also found associated with epidemics in tomato in 1995 in southern Portugal and 1997 in southern Spain, either alone or in mixed infections with TYLCV-Sar (Louro et al., 1996 ; Navas-Castillo et al., 1997 ). In this paper, we report the genetic characterization of infectious clones from TYLCV-Is isolates collected in Spain and Portugal and we describe evidence of natural recombination between TYLCV-Is and the begomovirus Tomato leaf curl virus (ToLCV).

Infectious clones that produced Bemisia tabaci-transmissible progeny were obtained as described in Navas-Castillo et al. (1999 ) from the TYLCV-Is-infected tomato samples SP72/97, collected in Spain in June 1997, and Port2/95, collected in Portugal in September 1995. The complete genome sequences of these clones were determined from both strands of the full-length inserts cloned in pSP72/97 and pPort2 [GenBank accession numbers AF071228 (2791 nt) and AF105975 (2793 nt)] by automated sequencing with pUC/M13 universal primers and specific primers based on partial sequences. The DNA sequences of SP72/97 and Port2/95 were compared with those of TYLCV-Sar isolate Sardinia (GenBank accession no. X61153) and TYLCV-Is isolates for which complete genome information was available: Japan Shizuoka (Japan-S), GenBank accession no. AB014346; Japan Aichi (Japan-A), AB014347; Israel, X15656; Israel/Mild, X76319; Cuba, AJ223505; Dominican Republic, AF024715; and Iran, AJ132711. Isolates SP72/97 and Port2/95 appeared to be closely related variants of TYLCV-Is and their genomes strongly resembled those of the Japan-A and Japan-S isolates (Kato et al., 1998 ) (Table 1). In addition, the isolates from the Iberian Peninsula and Japan were closely related to the Israel/Mild isolate and more distantly related to the Israel, Cuba, Dominican Republic and Iran isolates. Sequence identities among the IRs or complete genomes followed a similar pattern. Nucleotide sequence comparisons of the different ORFs provided further support for the observed relationships between the TYLCV-Is isolates; interestingly, ORFs C1 and C4 were less conserved than V1, V2, C2 and C3 (Table 1). The presence of such closely related variants in geographical regions so distantly separated can only be explained by movement of infected plant material or viruliferous B. tabaci individuals via plant trade or human movement. In fact, this was argued as the cause of the primary spread of Middle East TYLCV-Is to the Caribbean and then to USA (Polston & Anderson, 1997 ; Polston et al., 1999 ).


View this table:
[in this window]
[in a new window]
 
Table 1. Percentage identity between the nucleotide sequence of TYLCV-Is SP72/97 and those of other TYLCV-Is isolates and TYLCV-Sar isolate Sardinia

 
The IRs of SP72/97 and Port2/95 show features characteristic of begomoviruses, including the stem–loop structure with the conserved nonanucleotide sequence 5'-TAATATTAC-3' in the loop, where the breaking and joining site (ori) for rolling-circle replication occurs (for a review see Gutiérrez, 1999 ). In addition, a perfect conservation of the stem–loop structure with high thermodynamic stability ({Delta}G=-16·1 kcal/mol or -3·8 kJ/mol) (analyses performed using the DNA mfold server: M. Zuker, Washington University, St Louis, MO, USA; http://www.ibc.wustl.edu/~zuker/) was observed in the IRs of SP72/97, Port2/95 and all TYLCV-Is isolates listed above, except for Israel/Mild. In the latter, a less stable stem–loop structure was predicted ({Delta}G=-6·8 kcal/mol or -1·6 kJ/mol), owing to incomplete base-pairing in the stem. As point mutations that destroy base-pairing in the stem have been shown to be deleterious for virus replication (Orozco & Hanley-Bowdoin, 1966 ), the mildness reported for the Israel/Mild isolate in certain tomato cultivars (Antignus & Cohen, 1994 ) might be due to an altered stem–loop structure.

Phylogenetic relationships of SP72/97 and Port2/95 to other TYLCV-Is isolates (listed above) were studied on the basis of the nucleotide sequences of the IRs and ORFs V1, V2 and C1 to C4. As changes in the topological position of certain isolates occurred depending on which part of the genome was compared (not shown), a more detailed comparison throughout the genome was done with the aid of PLOTSIMILARITY diagrams (Wisconsin GCG software package; Devereux et al., 1984 ). In these comparisons, the nucleotide sequences of the SP72/97, Port2/95, Japan-A and Japan-S isolates on the one hand and the Israel, Cuba and Dominican Republic isolates on the other proved to be almost identical throughout the genome (not shown). Therefore, in Fig. 1, the first group of isolates is represented by SP72/97 while the second is represented by Israel. Comparisons revealed that a large portion of the genome, hereafter referred to as region I (see below), was quite well conserved among all TYLCV-Is isolates (Fig. 1). The nucleotide sequence of the Israel isolate differed from that of SP72/97 in the 5'-proximal two-thirds of the C1 ORF and the 5'-half of the IR (regions II and III; Fig. 1A). The Israel and Iran sequences also differed in the same part of the genome (Fig. 1B). However, when the sequences of SP72/97 and Iran were compared, differences were only observed in region III (Fig. 1C). The Israel/Mild sequence was highly similar to that of SP72/97 except for the central part of the IR (region IV; Fig. 1D). Nucleotide sequences in regions II to IV of TYLCV-Is isolates were then compared with other available geminivirus sequences, the limits of each region established by visual inspection of sequence alignments (indicated in Fig. 1). The sequence in region III of the Iran isolate was more closely related (91% identity) to the sequence of the equivalent region of the Bangalore-4 isolate of ToLCV (GenBank accession no. AF165098) than to that of any other TYLCV-Is isolate (below 75% identity). Sequences in regions II and III of the Israel, Cuba and Dominican Republic isolates were 97–100% identical, their closest relationships (between 73 and 83% identity) being to sequences in the equivalent regions of the Bangalore-2 (GenBank accession no. U38239) and Bangalore-4 isolates of ToLCV.



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 1. PLOTSIMILARITY diagrams (scanning window=50) comparing the nucleotide sequences of TYLCV-Is isolates SP72/97 and Israel (A), Israel and Iran (B), SP72/97 and Iran (C) and SP72/97 and Israel/Mild (D). Separation between regions I to IV, for which differential distribution of nucleotide identity is observed, is indicated by vertical dotted lines. The first nucleotide of the region in parentheses (numbers refer to nucleotide positions in the sequence of SP72/97 TYLCV-Is reported in this work). Positions of the ORFs and of the IR are indicated at the top of the figure. Horizontal broken lines indicate the mean similarities between the sequences compared.

 
Regions I to IV were then considered separately to infer the phylogenetic relationships among the TYLCV-Is isolates and the two ToLCV isolates indicated above. TYLCV-Sar and the A component of a New World begomovirus species, Tomato golden mosaic virus (TGMV) (GenBank accession no. K02029), were also used in the analyses, the latter as outgroup. The software packages PHYLIP 3.5 (J. Felsenstein, University of Washington, Seattle, WA, USA) and/or CLUSTAL W (Thompson et al., 1994 ) were used for these analyses. Both the neighbour-joining (Saitou & Nei, 1987 ) and unrooted parsimony (Fitch, 1971 ) methods yielded similar tree topologies in all cases. Representative phylogenetic trees obtained are shown in Fig. 2. Comparisons of sequences in regions I and IV produced similar trees except that a greater genetic divergence was observed for region IV, contained in the IR, a highly variable area between geminiviruses (Padidam et al., 1995 ). In both cases, all TYLCV-Is isolates formed a tight cluster related to TYLCV-Sar isolate Sardinia, all of them well separated from the Bangalore-2 and Bangalore-4 isolates of ToLCV [shown in Fig. 2(A) for region I], supporting the admitted species separation (Briddon & Markham, 1995 ). However, when comparing sequences in regions II and III, significant topological changes occurred in the position of certain TYLCV-Is isolates. For region II, most TYLCV isolates could be grouped together, but the isolates from Israel, Cuba and the Dominican Republic grouped best with Bangalore-2 ToLCV (not shown). When region III was considered, a similar tree topology was obtained, except that, in this case, the Iran isolate also grouped with ToLCV isolates, closely related to Bangalore-4 (Fig. 2B).



View larger version (43K):
[in this window]
[in a new window]
 
Fig. 2. Neighbour-joining phylogenetic trees depicting topology and distance relationships of sequences in regions I (A) and III (B) defined in Fig. 1. Trees were inferred using programs in the CLUSTAL W software package. The significance of the nodes in a bootstrap analysis with 1000 replicates is shown. Vertical lengths are arbitrary, while horizontal lengths reflect genetic distance between branch nodes. Scale bars indicate the horizontal distance equivalent to 0·1 replacements per position. TGMV was used as an outgroup. Isolates reported as members of the TYLCV-Is species are shaded.

 
The analysis described above indicated that sequences in regions II and III of the Israel isolate and region III of the Iran isolates were most closely related to the corresponding sequences of the Bangalore-2 and Bangalore-4 isolates of ToLCV (Fig. 2B). Harrison & Robinson (1999 ) have already suggested that ORF C1 and the IR of the Israel and Israel/Mild isolates may have been acquired from different (unidentified) parents, based on the low sequence similarity between them. In addition, Padidam et al. (1999 ) claimed to have statistical support for recombination events having occurred between TYLCV-Is and ToLCV. Our results suggest that the Israel and Iran isolates are hybrids between TYLCV-Is-like and ToLCV-like ancestors, the exchange having occurred between the 5' half of ORF C1 and the 5' half of the IR. The TYLCV-Is donors of the Iran and Israel chimeric genomes were closely related to each other and to other known TYLCV-Is isolates. However, the ToLCV donors differed from each other, although both were related to ToLCV isolates from India (Bangalore-2 and Bangalore-4). In nature, TYLCV-Is and ToLCV have common hosts such as tomato, so mixed infections can occur and they may replicate simultaneously in the same cell, which is a prerequisite for recombination.

The retention of the 5'-proximal part of the IR with the cognate 5' region of ORF C1 might be due to functional constraints. In fact, ORF C1 (rep gene) encodes the replication-associated protein (Rep) and it is known that the N-terminal region of Rep interacts specifically with the cognate 5' region of the IR for DNA replication and autoregulation of the rep gene (Fontes et al., 1994 ; Eagle et al., 1994 ; Orozco et al., 1997 ). The 3' limit of the putative ToLCV recombined fragments occurred in a highly variable region of the IR in both the Iran and Israel isolate genomes (not shown), so it was difficult to define the precise putative 3' cross-over site. However, in both cases, it occurred close to the conserved stem–loop structure of the IR, in which ori occurs (Fig. 1), and ori has been suggested as a hot-spot for recombination in the geminiviruses (Stanley, 1995 ; Sanz et al., 1999 ). The precise putative 5' cross-over site could not be defined by comparison of the Iran and Israel genomes with their closest related ancestors because stretches of identical nucleotides occurred in the region involved.

The existence of additional recombinations cannot be excluded. For example, in the genome of the Israel isolate, a short stretch of 113 nt was found within the central part of ToLCV-related region III that was almost identical (>96% identity) to the equivalent sequence of related TYLCV-Is isolates (see Fig. 1A for comparison with SP72/97).

Recombination is not a rare phenomenon among begomoviruses (Zhou et al., 1997 ; Harrison & Robinson, 1999 ; Padidam et al., 1999 ; Sanz et al., 1999 ), providing additional sources of variation with unpredictable effects on virus pathogenicity. Knowledge of sequence variants of TYLCV-Is present in a given geographical region is essential for the correct establishment of effective control measures, particularly because breeding is the best way to control this virus.


   Acknowledgments
 
We wish to thank E. R. Bejarano for providing some of the primers used for sequencing and R. G. Milne for critical reading of the manuscript. This work was financed by the Spain–Italy Integrated Action Programme (HI1997-0159), the CICYT, Spain (AGF98-0439-C05-02), the European Union (1FD97-0271) and Caja Rural de Málaga. S. Sánchez-Campos is the recipient of a predoctoral fellowship from the Ministerio de Educación y Cultura, Spain.


   Footnotes
 
The GenBank accession numbers of the sequences reported in this paper are AF071228 and AF105975.


   References
Top
Abstract
Main text
References
 
Antignus, Y. & Cohen, S.(1994). Complete nucleotide sequence of an infectious clone of a mild isolate of tomato yellow leaf curl virus (TYLCV).Phytopathology84, 707-712.

Briddon, R. W. & Markham, P. G.(1995). Family Geminiviridae. In Virus Taxonomy. Sixth Report of the International Committee on Taxonomy of Viruses, pp. 158-165. Edited by F. A. Murphy, C. M. Fauquet, D. H. L. Bishop, S. A. Ghabrial, A. W. Jarvis, G. P. Martelli, M. A. Mayo & M. D. Summers. Vienna & New York:Springer-Verlag.

Czosnek, H. & Laterrot, H.(1997). A worldwide survey of tomato yellow leaf curl viruses.Archives of Virology142, 1391-1406.[Medline]

Devereux, J., Haeberli, P. & Smithies, O.(1984). A comprehensive set of sequence analysis programs for the VAX.Nucleic Acids Research12, 387-395.[Abstract]

Eagle, P. A., Orozco, B. M. & Hanley-Bowdoin, L.(1994). A DNA sequence required for geminivirus replication also mediates transcriptional regulation.Plant Cell6, 1157-1170.[Abstract/Free Full Text]

Fitch, W. M.(1971). Toward defining the course of evolution: minimum change for a specified tree topology.Systematic Zoology20, 406-416.

Fontes, E. P. B., Eagle, P. A., Sipe, P. S., Luckow, V. A. & Hanley-Bowdoin, L.(1994). Interaction between a geminivirus replication protein and origin DNA is essential for viral replication.Journal of Biological Chemistry269, 8459-8465.[Abstract/Free Full Text]

Gutiérrez, C.(1999). Geminivirus DNA replication.Cellular and Molecular Life Sciences56, 313-329.[Medline]

Harrison, B. D. & Robinson, D. J.(1999). Natural genomic and antigenic variation in whitefly-transmitted geminiviruses (begomoviruses).Annual Review of Phytopathology37, 369-398.[Medline]

Kato, K., Onuki, M., Fuji, S. & Hanada, K.(1998). The first occurrence of tomato yellow leaf curl virus in tomato (Lycopersicon esculentum Mill.) in Japan.Annals of the Phytopathological Society of Japan64, 552-559.

Louro, D., Noris, E., Veratti, F. & Accotto, G. P.(1996). First report of tomato yellow leaf curl virus in Portugal.Plant Disease80, 1079.

Moriones, E. & Navas-Castillo, J. (2000). Tomato yellow leaf curl virus, an emerging virus complex causing epidemics worldwide. Virus Research (in press).

Navas-Castillo, J., Sánchez-Campos, S., Díaz, J. A., Sáez-Alonso, E. & Moriones, E.(1997). First report of tomato yellow leaf curl virus-Is in Spain: coexistence of two different geminiviruses in the same epidemic outbreak.Plant Disease81, 1461.

Navas-Castillo, J., Sánchez-Campos, S., Díaz, J. A., Sáez-Alonso, E. & Moriones, E.(1999). Tomato yellow leaf curl virus-Is causes a novel disease of common bean and severe epidemics in tomato in Spain.Plant Disease83, 29-32.

Noris, E., Hidalgo, E., Accotto, G. P. & Moriones, E.(1994). High similarity among the tomato yellow leaf curl virus isolates from the West Mediterranean Basin: the nucleotide sequence of an infectious clone from Spain.Archives of Virology135, 165-170.[Medline]

Orozco, B. M. & Hanley-Bowdoin, L.(1996). A DNA structure is required for geminivirus replication origin function.Journal of Virology70, 148-158.[Abstract]

Orozco, B. M., Miller, A. B., Settlage, S. B. & Hanley-Bowdoin, L.(1997). Functional domains of a geminivirus replication protein.Journal of Biological Chemistry272, 9840-9846.[Abstract/Free Full Text]

Padidam, M., Beachy, R. N. & Fauquet, C. M.(1995). Classification and identification of geminiviruses using sequence comparisons.Journal of General Virology76, 249-263.[Abstract]

Padidam, M., Sawyer, S. & Fauquet, C. M.(1999). Possible emergence of new geminiviruses by frequent recombination.Virology265, 218-225.[Medline]

Picó, B., Díez, M. J. & Nuez, F.(1996). Viral diseases causing the greatest economic losses to the tomato crop. II. The Tomato yellow leaf curl virus – a review.Scientia Horticulturae67, 151-196.

Polston, J. E. & Anderson, P. K.(1997). The emergence of whitefly-transmitted geminiviruses in tomato in the Western Hemisphere.Plant Disease81, 1358-1369.

Polston, J. E., McGovern, R. J. & Brown, L. G.(1999). Introduction of tomato yellow leaf curl virus in Florida and implications for the spread of this and other geminiviruses of tomato.Plant Disease83, 984-988.

Rochester, D. E., DePaulo, J. J., Fauquet, C. M. & Beachy, R. N.(1994). Complete nucleotide sequence of the geminivirus tomato yellow leaf curl virus, Thailand isolate.Journal of General Virology75, 477-485.[Abstract]

Saitou, N. & Nei, M.(1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees.Molecular Biology and Evolution4, 406-425.[Abstract]

Sanz, A. I., Fraile, A., Gallego, J. M., Malpica, J. M. & García-Arenal, F.(1999). Genetic variability of natural populations of cotton leaf curl geminivirus, a single-stranded DNA virus.Journal of Molecular Evolution49, 672-681.[Medline]

Stanley, J.(1995). Analysis of African cassava mosaic virus recombinants suggests strand nicking occurs within the conserved nonanucleotide motif during the initiation of rolling circle DNA replication.Virology206, 707-712.[Medline]

Thompson, J. D., Higgins, D. G. & Gibson, T. J.(1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.Nucleic Acids Research22, 4673-4680.[Abstract]

Zhou, X., Liu, Y., Calvert, L., Munoz, C., Otim-Nape, G. W., Robinson, D. J. & Harrison, B. D.(1997). Evidence that DNA-A of a geminivirus associated with severe cassava mosaic disease in Uganda has arisen by interspecific recombination.Journal of General Virology78, 2101-2111.[Abstract]

Received 12 May 2000; accepted 28 July 2000.