Frequent occurrence of Mungbean yellow mosaic India virus in tomato leaf curl disease affected tomato in Oman

Next generation sequencing (NGS) of DNAs amplified by rolling circle amplification from 6 tomato (Solanum lycopersicum) plants with leaf curl symptoms identified a number of monopartite begomoviruses, including Tomato yellow leaf curl virus (TYLCV), and a betasatellite (Tomato leaf curl betasatellite [ToLCB]). Both TYLCV and ToLCB have previously been identified infecting tomato in Oman. Surprisingly the NGS results also suggested the presence of the bipartite, legume-adapted begomovirus Mungbean yellow mosaic Indian virus (MYMIV). The presence of MYMIV was confirmed by cloning and Sanger sequencing from four of the six plants. A wider analysis by PCR showed MYMIV infection of tomato in Oman to be widespread. Inoculation of plants with full-length clones showed the host range of MYMIV not to extend to Nicotiana benthamiana or tomato. Inoculation to N. benthamiana showed TYLCV to be capable of maintaining MYMIV in both the presence and absence of the betasatellite. In tomato MYMIV was only maintained by TYLCV in the presence of the betasatellite and then only at low titre and efficiency. This is the first identification of TYLCV with ToLCB and the legume adapted bipartite begomovirus MYMIV co-infecting tomato. This finding has far reaching implications. TYLCV has spread around the World from its origins in the Mediterranean/Middle East, in some instances, in live tomato planting material. The results here may suggest that begomoviruses which do not commonly infect tomato, such as MYMIV, could be spread as a passenger of TYLCV in tomato.

Confirmation of NGS results. RCA products from the six tomato samples were digested with various restriction endonucleases. Restriction yielded ~2.7 kb DNA fragments from all samples, which were cloned in pUC19 and Sanger end sequenced. After partial sequencing and blast analysis the sequences formed three groups, representative of three distinct begomovirus species (results not shown). For each species a single clone from each tomato sample was fully sequenced using the primer-walking strategy. Analysis of these sequences is detailed in the next section.
Specific  , obtained from RCA products digested with Xbal, showed high levels of sequence identity to isolates of TYLCV available in the databases. Four sequences (Tom 30-3, Tom 31-5, Tom 26-10 and Tom , isolated from tomato plants in which MYMIV was identified, showed between 99.1 and 100% sequence identity to each other but only 94.6 to 95.4% identity to the two clones (Tom 32-29 and Tom 34-33) obtained from the tomato plants in which MYMIV was not identified. The four sequences showed the highest levels of sequence identity (99.3 to 100%) to an isolate of TYLCV strain "Iran" (TYLCV-IR) recently obtained from common bean in Oman (MG970362) 32 whereas the two other isolates showed the highest levels of sequence identity (99.9%) to two isolates of TYLCV-IR obtained from tomato originating from Oman (DQ644565, FJ956702) 2,19 . A closer analysis of the six TYLCV-IR sequences obtained here showed the two sequences from tomato in which MYMIV was not identified to differ from the four sequences in which MYMIV was identified in the intergenic region which contain numerous sequence changes as well as several insertions and deletions (results not shown). In a phylogenetic analysis the two sequences from tomato in which MYMIV was not identified segregate apart from the other four sequences ( Fig. 2A). An analysis of the six TYLCV sequences produced here from clones obtained by PCR amplification for possible recombination using the Recombination Detection Program (RDP) 33 showed the two clones (Tom 32-29 and Tom 34-33) obtained from the tomato plants in which MYMIV was not identified to differ from the four clones from the tomato plants in which MYMIV was identified. All sequences have potential recombination events in the intergenic region and the sequences obtained from tomato plants in which MYMIV was identified additionally having a recombination event at the N-terminal end of the C4 gene ( Supplementary Fig. 1). These recombination events possibly explain the differences seen between the two groups of TYLCV sequences in the phylogenetic analysis.
The sequences of six clones (Tom 26-12, Tom 30-19, Tom 31-23, Tom 32-30, Tom 34-34 and Tom 35-26) obtained from RCA product digested with Hindlll showed the highest levels (98.8 to 100%) of sequence identity to two ChiLCV isolates available in the databases; HG969264 isolated from radish originating from Oman 34 23 . In a phylogenetic analysis these sequence segregated with previously published sequences of ToLCBrV; a begomovirus species so far only identified in Oman ( Fig. 2A).
Six potentially full-length clones, Tom26-11, Tom30-9, Tom31-14, Tom32-7, Tom 34-38 and Tom35-2, obtained from PCR amplifications with universal betasatellite primers 36 were sequenced. These sequences are between 1,375 and 1,386 nt in length and are available in the GenBank database under the accession numbers given in Table 2. The six sequences show the features typical of betasatellites consisting of a single conserved open reading frame, in the complementary-sense, with the capacity to encode a product of 118 amino acids (Table 2), a region of sequence rich in adenine (coordinates 721-967) and a sequence conserved between all betasatellites, the satellite conserved region (coordinates 1255-16). All the betasatellite sequences obtained showed high nucleotide sequence identity (>90%) to ToLCB isolates previously reported from Oman. An analysis using the Species Demarcation Tool (SDT) showed the highest levels of identity to be to an isolate of ToLCB associated with Chili leaf curl virus from watermelon (KX787940) 35 for Tom31-14 and Tom35-2 (99.3 and 99.6%, respectively), to an isolate associated with Tomato leaf curl Barka virus from tomato (KF293292) 23 for Tom26-11 (96.6%) and to an isolate from tomato (KF229727) for Tom30-9 (94.5%). A phylogenetic analysis based upon an alignment with selected other betasatellite sequences from the databases showed the four betasatellite sequences obtained here to be most closely related to other isolates obtained from Oman, with which they form a separate clade, followed by ToLCB isolates from Iran (Fig. 2B).
Efforts to identify the presence of possible alphasatellites, either by PCR with universal primers for alphsatellites 37 or by cloning from RCA products, were uniformly negative.  www.nature.com/scientificreports www.nature.com/scientificreports/ Characterization of MYMIV clones obtained from tomato with ToLCD. The PCR amplicons obtained with primers specific for MYMIV DNA A and DNA B were cloned into pTZ57R/T. Eight potentially full-length clones, four (Tom26-4, Tom30-6, Tom31-1 and Tom35-8) obtained with the DNA A primers and four (Tom26-7, Tom30-13, Tom31-15 and Tom35-18) obtained with the DNA B primers, were selected for further analysis. The eight clones were sequenced in their entirety and the sequences are available in the GenBank sequence database under the accession numbers given in Table 2.
The sequences of clones Tom26-4, Tom30-6, Tom31-1 are 2,746 nt in length and that of Tom35-8 is 2,748 nt. The sequences have all characteristics typical of the genomes/DNA A components of begomoviruses originating from the Old World, encoding four genes in the complementary-sense and two in the virion-sense orientation ( Table 2). The four sequences showed greater than 99.5% nucleotide identity with each other while and the highest nucleotide sequence identities (99 to 99.3% %) with the DNA A components of MYMIV isolates previously identified in Oman 38,39 , followed by 99 to 99.3% identity to the DNA A sequence of a MYMIV isolate from cowpea in India (AY937195) 40 . Based on the presently applicable criteria for species demarcation of begomoviruses Tom26-4, Tom30-6, Tom31-1 and Tom35-8 are the DNA A component of isolates of MYMIV 41 . In a phylogenetic analysis the four sequences from tomato segregated with the DNA A components of other MYMIV isolates available in the databases, with the closest relationship, and forming a discrete clade, with isolates previously characterized from Oman, isolates originating from Indonesia, three isolates from India and an isolate from Bangladesh (Fig. 3A).
The sequences of the four clones obtained with MYMIV DNA B specific primers (Tom26-7, Tom30-13, Tom31-15 and Tom35-18) were each of 2,652 nt in length and are available in the GenBank database under the accession numbers given in Table 2. Analysis of the four sequences with ORF Finder showed them to have an arrangement typical of the DNA B components of bipartite begomoviruses, encoding one gene in each orientation ( Table 2) and to share greater than 99.2-99.8% nucleotide sequence identity. An SDT analysis of selected MYMIV and MYMV DNA B sequences available in the nucleotide sequence databases showed the DNA B sequences obtained from Oman to have the highest nucleotide sequence identities (98.8%) with the DNA B sequences of MYMIV isolates previously reported in Oman 38,39 followed by 98-98.5% to a MYMIV DNA B sequence isolated from cowpea in India (AY937196) 40 . A phylogenetic analysis of the DNA B sequences showed the isolates from tomato to be most closely related to isolates previously identified in Oman. All the isolates from Oman were most closely related to AY937196 40 and formed a distinct clade with the DNA B sequences of MYMIV isolates from Indonesia, some isolates from India as well as some isolates of MYMV (Fig. 3B).
The MYMIV sequences obtained by NGS agreed well with the sequences obtained from Sanger sequencing of clones. For example, for plant Tom 26 the DNA A determined by NGS showed 99.7% identity to Tom26-4 whereas the DNA B sequence showed 99.2% identity to Tom26-7. The sequences of the MYMIV DNA A and DNA B determined by NGS are given in Supplementary Fig. 2. The sequence differences likely are due to natural sequence variation in the populations of viral molecules in the plant.  www.nature.com/scientificreports www.nature.com/scientificreports/ Geographical incidence of MYMiV in tomato in oman. To Table 1). The samples were collected from across the country, including the main agricultural area in Al Batinah Governorate and in the south around Salalah (Dhofar province), showing the presence of MYMIV in tomato to be widespread.

Analysis of the infectivity of the MYMIV, TYLCV and ToLCB clones isolated from tomato.
Agrobacterium-mediated inoculation of just the DNA A component of the MYMIV isolate from tomato to N. benthamiana plants resulted no apparent symptoms of infection at 24 days post inoculation (dpi) ( Fig. 4; Table 3). However, for a single plant, out of seventeen inoculated plants, the presence of the DNA A component was shown by PCR diagnostics. Similarly there were no apparent symptoms following inoculation with both the DNA A and DNA B components of MYMIV, although more plants (n = 3) were found to harbour the DNA A at 24 dpi ( Fig. 5; Table 3). Of these only two plants were also found to harbour the DNA B component.
Agroinoculation  (Fig. 5). However, the symptoms were slightly milder, with less pronounced upward leaf curling, and appeared one day earlier (Table 3) Table 3). For plants inoculated with TYLCV and either the DNA A or DNA B of MYMIV, the symptoms of infection were indistinguishable from plants inoculated with only TYLCV and in none of the plants were the MYMIV components detected by diagnostic PCR. However, the latent period was possibly one day longer than for plants inoculated with only TYLCV (10 days in comparison to 9 days, respectively). Co-inoculation of TYLCV with both the DNA A and DNA B components of MYMIV induced symptoms that were qualitatively similar to N. benthamiana plants inoculated with only TYLCV. However, the symptoms were a little less severe and were delayed in comparison to plant inoculated with only TYLCV (12 days compared to 9 days, respectively; Fig. 5; Table 3). When ToLCB was coinoculated with TYLCV and MYMIV DNA A and DNA B, N. benthamiana plants developed symptoms indistinguishable from the symptoms induced by a TYLCV with ToLCB infection, although the symptoms were delayed by approx. 2 days (Fig. 5; Table 3).
Agroinoculation of tomato plants with either DNA A alone or with MYMIV DNA A and DNA B did not lead to symptoms and neither DNA A nor DNA B were detected by diagnostic PCR or by PCR reactions using RCA products, obtained from amplifications with DNA extracted from inoculated plants, as template (RCA-PCR; Fig. 5; Table 3). In contrast, agroinoculation of tomato with either MYMIV DNA A and ToLCB or with MYMIV DNA A, DNA B and ToLCB led to very late developing symptoms (approx. 28 dpi) consisting of mild downward leaf curling and mild yellowing for a few plants ( Fig. 5; Table 3). Neither of the virus components nor the betasatellite were detected in diagnostic PCR. However RCA-PCR showed the presence of the virus components in the plants with symptoms. The TYLCV isolate obtained here was highly infectious to tomato and induced severe downward leaf curling and yellowing symptoms by 12 dpi in all inoculated plants. Co-inoculation of TYLCV with ToLCB reduced the latent period (10 days) and altered the symptoms to mild downward leaf curling with slightly more yellowing. Maintenance of the betasatellite by TYLCV was shown in 14 out of 16 tomato plants inoculated by diagnostic PCR (Fig. 5; Table 3  www.nature.com/scientificreports www.nature.com/scientificreports/ A, MYMIV DNA B and ToLCB) induced symptoms that were qualitatively equivalent to plants inoculated with TYLCV with ToLCB but slightly milder and with a slightly longer latent period (12 days) than in the absence of the betasatellite. Again MYMIV DNA A (in four plants) and MYMIV DNA B (in two plants) were detected in these plants although both TYLCV and ToLCB were detected in the majority of plants.

Discussion
During a small survey of tomato crops conducted during 2015 in the Al-Batinah region of Oman using NGS, the presence in some plants of monopartite begomoviruses (TYLCV, ChiLCV, ToLCBrV), a betasatellite (ToLCB) and a bipartite begomovirus (MYMIV) was shown. The NGS results were confirmed by cloning the viruses/virus components and betasatellite, either from RCA product or by PCR-mediated amplification. TYLCV (specifically the "Iran" and "Oman" strains), ChiLCV and ToLCBrV are common monopartite begomoviruses infecting tomato, causing ToLCD, in Oman and are frequently associated with ToLCB, the only betasatellite so far shown to infect tomato in Oman 19,20 . The finding of these viruses and betasatellite in tomato is thus not surprising. However, the identification of co-infection of these viruses and betasatellite with MYMIV in tomato is surprising.
The incidence of ToLCD in the survey fields in Al-Batinah during 2015 was between 65 and 70% and of the six tomato plants examined by NGS four were found to harbour MYMIV. A more extensive survey for MYMIV in tomato, for samples collected in 2015 and 2016, indicated that approx. 25% of symptomatic plants harboured MYMIV and showed that the phenomenon was occurring not only in the major agricultural area of Al-Batinah (north Oman) but also in the agricultural area in southern Oman (Dhofar region).
In the phylogenetic analyses the MYMIV DNA A and DNA B components originating from Oman, including those identified earlier 38,39 , form a distinct clade with MYMIV isolates from Indonesia and some isolates from India and Bangladesh. All these isolates are unusual in that, as first noted by John et al. 40 and later confirmed for the isolates from Indonesia by Tsai et al. 30 , they are MYMIV but with a DNA B component derived from the closely related MYMV, likely by a process known as pseudo-recombination (component exchange) 26 . This was confirmed for the Oman isolates by alignments, which show the majority of the DNA B sequence to have the highest levels of identity to the sequences of MYMV DNA B components but the intergenic region to have higher levels of sequence identity to MYMIV isolates (results not shown).
Based on their results Tsai et al. 30 concluded that MYMIV was introduced into Indonesia only very recently. The collective results from Oman similarly suggest that MYMIV has only been introduced quite recently. The very close relationship of the isolates in Oman and Indonesia also suggests that they have a common origin, likely having been introduced from South Asia. The likely mechanism of introduction is unclear, although the demonstration for an increasing number of geminiviruses that seed transmission can occur, including for the LYMVs DoYMV 42 and possibly MYMV 43 , which are closely related to MYMIV, may suggest that the mechanism of geographic dispersal was via seed.
In Oman both TYLCV-IR and MYMIV have previously been identified separately in legumes; common bean and kidney bean, respectively 32,39 . Although generally considered a tomato virus, TYLCV can infect and causes problems in legumes; for example the type strain of TYLCV is one of two viruses shown to cause bean leaf crumple disease in Spain 44 . The unusual aspect of the findings of Shahid et al. 39 was the identification of MYMIV with a betasatellite (ToLCB) in kidney bean. This virus is rarely encountered in the presence of a betasatellite.   Betasatellites have previously been shown to complement DNA B functions suggesting that the only protein encoded by betasatellites, βC1, has movement functions [50][51][52] . However, the demonstration that βC1has suppressor of both post-transcriptional and transcriptional gene silencing 17 may indicate that rather than providing movement functions, betasatellites instead overcome a host resistance to movement based on gene silencing. This is supported by the finding here that inoculation with the betasatellite reduces the latent period (faster spread of the virus following inoculation) and alters the phenotype of TYLCV infection (indicating differing plant tissues affected and/or timing of the infection of those tissues), as noted previously for a number of monopartite begomovirus-betasatellite interactions including TYLCV/ToLCB 2,53-55 .
In www.nature.com/scientificreports www.nature.com/scientificreports/ and TYLCV/ToLCB which may suggest that the MYMIV components are affecting the infection. However, this aspect will require further investigation.
The earlier identification of both TYLCV-IR/ToLCB 32 and MYMIV/ToLCB 39 in Phaseolus vulgaris (kidney bean and common bean, respectively) in the same geographic area provides a likely route (insect transmission from legumes) for the co-infection of tomato occurring. This is supported by the finding here that the TYLCV-IR isolates from tomato in which MYMIV was not identified show closer relatedness to TYLCV-IR isolates earlier identified in tomato, whereas TYLCV-IR isolates from tomato plants in which MYMIV was identified show closer relatedness to TYLCV-IR isolates very recently identified in a legume. Whether this indicates that maintenance of MYMIV in tomato relies on a specific variant of TYLCV will require further investigation.
The host range of MYMIV does not extend to tomato 24 . In co-infected plants it is likely that proteins which suppress host defenses, including gene silencing (also known as RNA interference), encoded by TYLCV and/ or the betasatellite allow MYMIV to infect tomato; a phenomenon known as trans-complementation 17,[56][57][58][59][60] . However, the likely effects and stability of the TYLCV/ToLCB-MYMIV interaction are less clear. Previous studies of the interaction of bipartite begomoviruses with betasatellites has suggested that these might not be stable, either because of a conflict between the betasatellite and the DNA B 51 or due to a conflict between the virus encoded coat protein and the betasatellite 52 . Whether the TYLCV/ToLCB-MYMIV co-infection will prove to be stable in the field is unclear, although the low levels of MYMIV components detected in tomato in the field and by agroinoculation would suggest not.
The findings here with MYMIV in tomato in Oman very much mirrors the situation with Tomato leaf curl New Delhi virus (ToLCNDV) in cotton in Pakistan reported recently 61 . Cotton leaf curl disease in South Asia is caused by various monopartite begomoviruses with a specific betasatellite, Cotton leaf curl Multan betasatellite (CLCuMuB; reviewed by Sattar et al. 62 ). Zaidi et al. 61 showed that recently, at high incidence and occurring across a wide area, CLCuD affected cotton also harbours the bipartite begomovirus ToLCNDV. The study additionally showed that, in co-infected cotton, the levels of CLCuMuB were raised. Since the betasatellite encodes a dominant symptom determinant, enhanced betsatellite titres could result in greater losses to cotton production. Whether the enhanced levels of betasatellite are also occurring in the TYLCV/ToLCB-MYMIV interaction will require further investigation.
Although the potential losses to legume and tomato production from MYMIV infection in Oman are clearly of concern, the results presented here have more far reaching implications. TYLCV has been spread around the world due to the spread of the vector whitefly, the international trade in tomato seedlings and possibly even tomato fruit 4,63,64 . The fear that a disease enhancing, TYLCV-adapted betasatellite such as ToLCB could spread with the virus, has been noted previously 18,65 . The results here raise the additional possibility that other viruses, including viruses not adapted to plants of the Solanaceae such as MYMIV, could be moved around the world as a passenger to TYLCV in tomato. In that regard it is pertinent to mention Watermelon chlorotic stunt virus (WmCSV); a bipartite OW begomovirus that has recently been shown to be present, and causing problems, in Mexico 66 . WmCSV occurs across the Arabian Peninsula, North Africa and the Middle East and has been identified co-infecting tomato with Squash leaf curl virus, a bipartite begomovirus originating from the NW, in Jordan 67 . The study, however, did not show the presence of TYLCV. Some evidence was obtained in the study here to indicate that one tomato plant (Tom31) was infected with WmCSV -partial sequences of the DNA A identified in the NGS (Supplementary Fig. 3), although this was not confirmed by cloning. The mechanism of introduction of WmCSV into the NW is unclear and it is intriguing to speculate that it might have been by co-infection with TYLCV in tomato.
An important question that needs to be answered with respect to tomato acting as a "vector" for the spread of viruses other than TYLCV is whether coinfected plants can act as a source of, in this case, MYMIV for B. tabaci insects to transmit to other plant species when the "passenger" virus is apparently at such low titre. This will be the focus of future studies.

Materials and Methods
Detection of viruses and betasatellite. During a survey for begomoviruses in January 2015, tomato plants showing typical ToLCD infections were collected from the Barka Governorate, Al-Batinah region of Oman. Total nucleic acid was isolated from tomato leaf samples using a CTAB-based method 68 . DNA extracts were resuspended in sterile distilled water and stored at −20 °C. The nucleic extracts were used as a template in polymerase chain reaction (PCR) in initial screening with primer pair TYLCD-356/TYLCD-1044 69 , designed to conserved sequences of tomato-infecting begomoviruses in Oman.
next-generation sequencing. Circular DNA molecules in nucleic acid samples were amplified by rolling circle amplification (RCA) using a TempliPhi 100 Amplification System kit (GE Healthcare) as described by the manufacturer. Resulting high molecular weight products were sent to the Beijing Genomics Institute (Hong Kong) for sample preparation and paired-end 101 sequencing on an Illumina HighSeq. 4000 platform. NGS data received in fastq format was processed by using Trimgalore (https://www.bioinformatics.babraham. ac.uk/projects/trim_galore/) on a High Performance Computing cluster and adapter sequences and low quality (Phred score <20) sequence reads were removed. Sequence reads were mapped to the tomato genome sequence (Solanum lycopersicum iTAG2.4) using Bowtie2 70 to remove host related sequences and enrich viral sequences. Sequence reads not mapping to the tomato genome sequence were isolated and de-novo assembled using SPAdes version 3.12.0 71 .Contigs obtained from this assembly were identified by using stand-alone Blast search tool. These primers for the amplification of full-length components were designed to the sequences obtained from the NGS analysis. Betasatellites were amplified with primer pair beta01/beta02 36 . The amplification products were cloned in pTZ257R/T (Fermentas). Concatameric RCA products were digested with restriction endonucleases XbaI or PstI. Resulting ~2.7 kb fragments were purified from agarose gels using a GeneJet Gel Extraction Kit (Thermo Fisher Scientific) and cloned in pUC19 72 .
Sanger sequencing, sequence assembly and sequence analysis. Selected clones harbouring potentially full-length begomovirus or satellite inserts were sequenced commercially (Macrogen Inc., South Korea). Sequence reads were assembled using SeqMan, part of the Lasergene package of sequence analysis software (DNA Star Inc., Madison, WI, USA). The resultant sequences were initially analyzed using BLASTn (http://blast.ncbi. nlm.nih.gov/Blast.cgi) 73 to identify closely related sequences in the nucleotide sequence databases. Percentage sequence identity values quoted were determined using SDT with the MUSCLE alignment option 74,75 . Pairwise multiple sequence alignments using the MUSCLE algorithm were produced using MEGA6 76 . Evolutionary relationships were determined by constructing phylogenetic trees using Clustal X (neighbor-joining method 77 ) and displayed using Treeview 78 . Possible recombination events in sequences were determined using RDP4 33 . Open reading frames (ORFs) in sequences were identified using ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder).

Production of constructs for Agrobacterium-mediated inoculation. A full-length clone of MYMIV
DNA A (Tom26-4) was digested with HindIII and XhoI to release a fragment of ~1,600 bp that was ligated to into the binary vector pGreen0029 79 to yield pGMYA-0.4. Then the full-length insert of Tom26-4 was released using HindIII and ligated into pGMYA-0.4, linearised with XhoI, to yield the partial direct repeat construct pGMYA-1.4. Similarly, MYMIV DNA B (Tom26-7) was digested with BamHI and ClaI to release a fragment of ~1, 700 bp and ligated to into pGreen0029 to yield pGMYB-0.4 into which the full-length insert of Tom26-7, released with BamHI, was ligated to yield pGMYB-1.4. A partial direct repeat construct of TYLCV (Tom26-10) in pGreen0029 was produced using an ~1.249 bp BamHI -XbaI fragment and the full-length insert of Tom26-10, released using XbaI, to yield pGTY-1.4. A partial direct repeat construct of ToLCB (Tom26-11) in pGreen0029 was produced using an ~600 bp BamHI-XbaI fragment and the full-length insert released using KpnI to yield pGTB-1.4. The pGreen0029 constructs were finally electroporated into Agrobacterium tumefaciens strain LBA4444.

Agrobacterium-mediated inoculation and maintenance of plants. Agrobacterium cultures har-
bouring the pGreen0029 constructs were grown for 48 h in 50 ml LB liquid medium supplemented with which antibiotics (kanamycin 25 µg/ul, tetracycline 10 µg/ul, rifampicin 50 µg/ul). The bacteria were pelleted by centrifugation at 4800 g for 10 min at 4 °C and resuspended in infiltration solution (10 mM MgCl 2 , 200 µM acetosyringone) and adjusted to a final optical density at 600 nm of ~1.0. Agroinoculation was performed on 4-5 week old N. benthamiana and tomato (Solanum lycopersicum L. cv. Pusa Ruby) plants as described previously 2 . For each combination of begomovirus and/or betasatellite, 5-6 plants were used for inoculation. Inoculated plants were maintained in a secure, insect-free growth room with a day length of 18 h, a day/night temperature of 25/26 °C and a relative humidity of 65%. The plants were observed daily for the appearance of symptoms of virus infection.