Recombination and pseudorecombination driving the evolution of the begomoviruses Tomato severe rugose virus (ToSRV) and Tomato rugose mosaic virus (ToRMV): two recombinant DNA-A components sharing the same DNA-B
Virology Journal volume 11, Article number: 66 (2014)
Begomoviruses are dicot-infecting, whitefly-transmitted viruses with a genome comprised of one or two molecules of circular, single-stranded DNA. In Brazil, tomato-infecting begomoviruses have emerged as serious pathogens since the introduction of a new biotype of the insect vector in the mid-1990’s. Tomato rugose mosaic virus (ToRMV) and Tomato severe rugose virus (ToSRV) are often found in tomato fields. The complete sequence of the DNA-B components of ToSRV and ToRMV show an identity of 98.2%. Additionally, the high nucleotide identity (96.2%) between their common regions indicates that these two viruses may share the same DNA-B.
Tomato seedlings were biolistically inoculated with ToSRV (DNA-A and DNA-B) and ToRMV (DNA-A and DNA-B) infectious clones in every possible combination of single or mixed infection. Symptom expression was evaluated for up to 35 days post-inoculation (dpi). DNA was extracted at 28 dpi and the presence of each viral genomic component was examined by rolling circle amplification (RCA) followed by digestion, as well as by quantitative, real-time PCR. Sequence comparisons, recombination and phylogenetic analyzes were performed using EMBOSS needle, RDP program and maximum likelihood inference, respectively.
Symptoms in tomato plants inoculated with the different combinations of ToRMV and ToSRV DNA-A and DNA-B components consisted of a typical mosaic in all combinations. Pseudorecombinants were formed in all possible combinations. When two DNA-A or two DNA-B components were inoculated simultaneously, the ToRMV components were detected preferentially in relation to the ToSRV components. The combination of minor changes in both the Rep protein and the CR may be involved in the preferential replication of ToRMV components. Recombination and phylogenetic analyzes support the exchange of genetic material between ToRMV and ToSRV.
ToRMV and ToSRV form viable pseudorecombinants in their natural host (Solanum lycopersicum) and share the same DNA-B. ToRMV DNA components are preferentially replicated over ToSRV components. These results indicate that the emergence of ToRMV involved both recombination and pseudorecombination, further highlighting the importance of these mechanisms in the emergence and adaptation of begomoviruses.
The genus Begomovirus (family Geminiviridae) contains viruses with a genome comprised of circular, single-stranded DNA, transmitted by the whitefly Bemisia tabaci (Homoptera: Aleyrodidae) to dicot plant species . Most New World begomoviruses have bipartite genomes consisting of two ssDNA components named DNA-A and DNA-B, each with approximately 2600 nucleotides (nt). These components share a common region (CR) of approximately 200 nt with high sequence identity (>85%, with a few exceptions) . The DNA-A encodes viral proteins associated with replication, control of gene expression, suppression of host defense responses and encapsidation, whereas the DNA-B encodes proteins associated with intra- and intercellular movement .
The CR includes the origin of replication as well as conserved sequences that are recognized by the replication-associated protein (Rep) [4, 5]. Recognition by Rep is considered to be virus-specific , so that Rep only initiates replication of cognate DNA components. However, in some cases the Rep protein encoded by a given DNA-A will recognize the DNA-B of a different virus, a phenomenon known as pseudorecombination . Pseudorecombinants are usually formed with DNA components from strains of the same viral species [7, 8], although, more rarely, they can also be formed with components from different species [9–11].
The viability of a pseudorecombinat is, at least partially, a function of the conservation of the high affinity binding sites for the Rep protein. For most begomoviruses these binding sites include two direct repeats and one inverted repeat known as iterons . DNA components with identical iterons will usually form viable pseudorecombinants, whereas those with different iterons will not [12, 13]. However, in both cases there are exceptions [9, 14], indicating that other, less understood factors also play a role in pseudorecombinant formation. The Rep protein of begomoviruses has specificity determinants (SPDs), which confer specificity to the recognition of iterons . Mutation or deletion of SPDs located adjacent to motif 1, known as the iteron-related domain (IRD) , eliminates the ability of the Rep protein to bind DNA in a sequence-specific manner [16, 17].
Genetic exchanges through recombination allow for a rapid evolution of plant viruses, in many cases promoting changes in virulence or host range . The role of recombination in the emergence, establishment and evolution of novel, better adapted begomoviruses is well established [19–27]. In Spain, a recombinant between Tomato yellow leaf curl virus (TYLCV) and Tomato yellow leaf curl Sardinia virus (TYLCSV) displayed a wider host range compared to both parental viruses, and became established in the field . In Uganda, a recombinant between African cassava mosaic virus (ACMV) and East African cassava mosaic virus (EACMV) displayed increased virulence and aggressiveness, and was responsible for devastating epidemics of cassava mosaic disease during the 1990’s [25, 26]. The combined effect of pseudorecombination and recombination in the emergence of novel begomoviruses was elegantly demonstrated by the work of Hou & Gilbertson . Working with Bean dwarf mosaic virus (BDMV) and Tomato mottle virus (ToMoV), the authors showed that, after five consecutive passages in Nicotiana benthamiana, a pseudorecombinant with ToMoV DNA-A and BDMV DNA-B underwent a recombination event in which part of the ToMoV DNA-A CR was transferred to the BDMV DNA-B. This resulted in an marked increase in the titer of the DNA-B and in the induction of more severe symptoms by the pseudorecombinant.
The incidence and severity of begomovirus epidemics in tomatoes has greatly increased in Brazil since the mid-1990’s, following the introduction and dissemination of the B biotype of B. tabaci . The initial characterization of begomoviruses associated with these epidemics indicated a high degree of genetic diversity [30, 31]. Several new species have been described, including Tomato rugose mosaic virus (ToRMV) , Tomato chlorotic mottle virus (ToCMoV)  and Tomato yellow spot virus (ToYSV) . Mixed infections, a pre-requisite for the occurrence of recombination and pseudorecombination events, are common in the field [27, 35, 36]. Our comprehensive analysis of Brazilian begomovirus populations provided strong support for the hypothesis that tomato-infecting begomoviruses evolved from indigenous viral populations present in non-cultivated hosts, and indicated that recombination played a major role in the emergence of these viruses as tomato pathogens .
ToRMV and Tomato severe rugose virus (ToSRV) are often found in tomato fields in Southeastern Brazil [27, 32, 35]. The complete DNA-B sequences of ToRMV (GenBank: AF291705) and ToSRV (DQ207749) show an identity of 98.2%, suggesting that these two viruses may share the same DNA-B. Additionally, recombination analysis indicated that part of the CR and most of the Rep gene were transferred from ToSRV to ToRMV .
The results presented here indicate that ToRMV and ToSRV isolates are capable of forming viable pseudorecombinants in their natural host, tomato (Solanum lycopersicum). Symptom severity was equivalent in single or mixed infections, indicating that synergism does not occur between these two viruses. However, ToRMV DNA-A and -B components were preferentially replicated over ToSRV components. These results further highlight the importance of recombination and pseudorecombination in the emergence and adaptation of begomoviruses, and pose a challenging taxonomical question: To what extent should a recombinant be considered a new species?
Latent period, symptoms and infectivity of the different combinations between ToRMV and ToSRV DNA-A and DNA-B
Symptoms in tomato plants inoculated with the different combinations of ToRMV and ToSRV DNA-A and DNA-B first appeared at 14 dpi (Table 1), and consisted of a typical mosaic in all combinations (Figure 1). In plants inoculated with ToYSV (used as a positive control) the first symptoms appeared at 10 dpi (Table 1), and included yellow mosaic, yellow spots and leaf distortion (Figure 1). The presence of each inoculated DNA component was verified by RCA followed by digestion with component-specific enzymes (Table 1 and Additional file 1: Figure S1). In single infections, infectivity of the two homologous combinations (ToRMV and ToSRV) was 92% and 100%, respectively. Pseudorecombinants were formed in both possible combinations with high infectivity (100% and 92%). Interestingly, when two DNA-A or two DNA-B components were inoculated simultaneously, the ToRMV component was detected preferentially in relation to the ToSRV component (Table 1, underlined values; Additional file 1: Figure S1). This was true even in the two cases in which ToSRV DNA-A and -B were inoculated with either ToRMV DNA-A or -B, suggesting that the Rep protein encoded by ToSRV may have a higher affinity for the ToRMV binding sites. Moreover, when all four components were inoculated, the two ToRMV components were detected in 100% of the plants, but the ToSRV components were detected in only 15% and 8% of the plants (DNA-A and DNA-B, respectively), again suggesting a preferential replication of ToRMV DNA components over ToSRV components.
To confirm the negative interference of ToRMV over ToSRV, primers for quantitative, real-time PCR (qPCR) were designed based on the most divergent region of the DNA-A of each virus. The primers designed for the detection of ToRMV DNA-A were not specific, yielding multiple, non-specific fragments from ToRMV-infected samples (data not shown). Therefore, only ToSRV DNA-A accumulation was quantified in single or dual infected plants (in which both DNA-A’s were co-inoculated). The qPCR assay revealed a strong reduction of ToSRV DNA-A accumulation in plants co-inoculated with ToRMV DNA-A (Additional file 2: Figure S2). Due to the high nucleotide sequence identity between the DNA-B’s of both viruses (98.2%), it was not possible to design DNA-B-specific primers. However, the RCA-RFLP assay for the DNA-B strongly suggests a similar pattern to that of DNA-A accumulation.
Together, these results indicate that ToSRV and ToRMV form viable pseudorecombinats in tomato plants and that synergism does not occur between these two viruses. Quite the contrary, probably ToRMV has a negative interference over ToSRV, although ToSRV can replicate and induce symptoms as efficiently as ToRMV in a single infection.
Nucleotide and deduced amino acid sequence comparisons
In an attempt to identify genomic differences that may explain the negative interference of ToRMV over ToSRV, comparisons of their coding and non-coding regions were carried out. The CRs of the two viruses show sequence identities ranging from 95% (ToSRV-A vs. ToRMV-B, and ToSRV-B vs. ToRMV-A) to 97.5% (ToSRV-A vs. -B), with identical iterons (Figure 2A). Alignment of the four CRs indicated nine positions with nucleotide sequence divergence (Figure 2A). At divergent positions 1, 5 and 6, the sequence is conserved between the cognate DNA components of ToRMV and ToSRV. At other positions one component differs from the other three: ToRMV-A at positions 4 and 7, ToRMV-B at positions 3 and 8, ToSRV-A at position 2, and ToSRV-B at position 9 (Figure 2A).
The comparison of nucleotide and amino acid sequences of the Rep gene indicates identities of 95.3% and 93.5%, respectively (Table 2). This was expected since the region involved in the recombination event comprises practically the entire Rep gene. The domains associated with sequence-specific recognition of the origin of replication, named “iteron-related domain” (IRD) , “Motif 1”, “Motif 2” and specificity determinants (SPDs) , conserved in rolling-circle replication-initiator proteins including geminivirus Rep proteins, are shown in Figure 2B. ToRMV has four amino acid changes compared with ToSRV. Two changes are in the IRD (Pro → Thr and Lys → Arg). These differences are also present in 85 aditional isolates which have been sequence in our laboratory over the years (data not shown).
To further investigate the occurrence of recombination events involving ToSRV and ToRMV, recombination analysis was carried out using the RDP3 package. Besides ToRMV and ToSRV, the analysis included seventeen closely related begomoviruses that occur in cultivated and non-cultivated plants (Table 3).
A recombination event involving ToSRV and ToRMV previously proposed by Ribeiro et al.  was confirmed by six detection methods contained in the RDP3 package (P-values: rdp = 1.667 × 10−14, Geneconv = 1.405 × 10−25, Bootscan = 3.538 × 10−28, maximum χ2 = 1.024 × 10−20, Chimaera = 2.975.45 × 10−12 and Sister scan = 1.171 × 10−38). The transferred region comprises practically the entire Rep gene as well as part of the CR including the iterons (Figure 3). Also as previously reported , a recombinant event was detected in which ToCMoV received a fragment of ToRMV DNA-A comprising the entire CP and Ren genes, as well as most of the Trap gene. This event was detected by all detection methods (P-values: rdp = 3.039 × 10−39, Geneconv = 3.734 × 10−37, Bootscan = 1.157 × 10−40, maximum χ2 = 2.085 × 10−21, Chimaera = 3.101 × 10−25, Sister scan = 4.771 × 10−43 and 3Seq = 5.053 × 10−64).
Phylogenetic analysis based on non-recombinant and recombinant nucleotide sequences between ToRMV and ToSRV in the DNA-A demonstrates topological incongruence, supporting the exchange of genetic material between the two viruses. ToRMV was placed in the same branch with ToCMoV in the tree based on the non-recombinant region (Figure 4A) and with ToSRV in the tree based on the recombinant region (Figure 4B), in both cases with high bootstrap values, indicating a close phylogenetic relationship. In the phylogenetic tree based on the DNA-B, ToRMV and ToSRV formed a well-supported clade with branches of the same size and a bootstrap value of 100% (Figure 4C). In addition, the high identity between their DNA-B sequences (nt and deduced aa; Table 2) suggests that the capture of the ToSRV DNA-B by ToRMV was a recent event.
Here, we report the formation of pseudorecombinants between two distinct, albeit closely related, begomoviruses. The reciprocal pseudorecombinants ToRMV-A + ToSRV-B and ToRMV-B + ToSRV-A were both viable. Their viability is most likely due to the high sequence identity of their CRs and Rep proteins, as well as the identical sequences of their iterons (GGTAG). The latent period and symptoms were the same in all combinations compared to the parental viruses, which is not surprising since the severity of symptoms is associated with the DNA-B, and these are almost identical between the two viruses.
When DNA-A and DNA-B components of both ToSRV and ToRMV were inoculated, mosaic symptoms identical to those induced by each virus alone appeared at 14 dpi (Figure 1 and Table 1). This suggests that synergism does not occur between these two viruses. On the contrary, ToRMV had a negative interference over ToSRV. In most cases in which three or four DNA components were inoculated, the ToRMV components were always detected, whereas the ToSRV components were rarely detected. Negative interference exerted by ToRMV in mixed infection with ToYSV has already been reported . Early in infection ToRMV negatively affected the accumulation of ToYSV, but when the infection was established the negative interference ceased . Interestingly, although ToRMV negatively affected ToSRV in our experiments, ToSRV behaved similarly to ToRMV in a single infection, both in terms of latent period and symptoms.
Considering that their DNA-B components are almost identical, and that their iterons are identical, we hypothesized that the preferential replication of ToRMV over ToSRV would be due to differences in their Rep proteins, most likely in their virus-specific recognition domains. Begomovirus iterons are composed of an invariable GG sequence followed by three nt (identified as N1, N2 and N3) that vary for each virus. The virus-specific recognition domain of the Rep protein maps to its N-terminal region [40, 41] and includes the conserved motif 1 of rolling-circle replication-initiator proteins . The IRD sequence is conserved among begomoviruses with identical iterons, but varies among species with different iteron sequences. Predicted nt-aa pairing would occur between N1 and the eighth aa of the IRD, N2 with the sixth IRD aa, and N3 with the first or third IRD aa, depending on the iteron sequence .
Although ToRMV and ToSRV have identical iterons, their IRDs are different at positions 2 and 3. The change of lysine to arginine at IRD position 3, which pairs with N3, should not promote major changes in the protein, since the two amino acids have similar properties. Therefore, it is reasonable to assume that this difference is not responsible for the preferential replication of ToRMV over ToSRV. In fact, this change does not affect iteron recognition, since the two Rep proteins have the ability to replicate viral DNA. This is observed when any combination of a DNA-A and a DNA-B component (cognate or pseudorecombinant) is inoculated: both components are always detected (Table 1).
Further analysis revealed small differences between ToRMV and ToSRV in the CR and in the Rep protein. In the CR, nine nucleotide changes occur among the four components, with ToRMV DNA components having the largest number of modifications compared to ToSRV. It is possible that these changes in the CR sequence, as well as changes in IRD aa 2 and 3, have enabled the ToRMV Rep protein to recognize the iteron sequence GGTAG with higher affinity compared to ToSRV. These same differences could cause the ToSRV Rep protein to replicate the DNA-A and -B of ToRMV preferentially. More likely is that the combination of minor changes in both Rep and CR may be involved in the preferential replication of ToRMV components. A detailed mutagenesis analysis of the IRD region and of the iterons of ToRMV and ToSRV will be required to confirm this hypothesis.
Interestingly, ToSRV is more prevalent in the field compared to ToRMV [27, 35]. This could be due to differences in transmission by the insect vector, since the coat protein, which is considerably distinct between the two viruses, is the major viral determinant of insect transmission.
The detection of a recombination event involving the DNA-A of the two viruses was not surprising, since it is well established that begomoviruses are recombination-prone. The region transferred from ToSRV to ToRMV comprises practically the entire Rep ORF as well as part of the CR including the TATA box, G-box, the conserved nonanucleotide (TAATATTAC) and the iterons. This recombination event between ToSRV and ToRMV was detected also by Ribeiro et al. , who proposed that ToRMV would be a recombinant between ToCMoV and ToSRV.
The specific genomic region involved in a recombination event is an important aspect in the maintenance and adaptation of the recombinant virus in the population. It is known that the begomovirus genome has a number of recombination hot spots, including the CR, the CP/Ren interface and the 3′ region of the Rep gene, while the remaining coding regions are recombination cold spots [28, 43–45]. The breakpoints of the recombination event involving ToSRV and ToRMV map to two hotspots, the CR and to the 3′ region of the Rep gene, near the interface with the Trap gene. Recombination in coding regions probably occurs due to the conflict between the machineries of transcription and replication [45, 46]. Lefeuvre et al.  suggest that selection apparently disfavors recombinants with breakpoints in coding regions, and that the probability of such a recombinant to remain in the population is dependent on the specific position of the breakpoint: the peripheries of coding regions tend to be more favorable to recombinantion than central regions. However, recent experimental evidence demonstrated that the coding regions of the begomoviruses Tomato yellow leaf curl virus-Mld (TYLCV) and Tomato leaf curl Mayotte virus (ToLCYTV) are recombination prone, and that all recombinants between TYLCV and ToLCYTV generated using L-DNA-shuffling technology displayed infectivity and viral accumulation similar to the parental begomoviruses .
Recombination events involving DNA-A components are considered a major source of molecular variation for begomoviruses and may result in a gain of virulence [49, 50]. However, predicting the epidemiological outcome of a recombination event is not a straightforward task. For example, van der Walt et al.  demonstrated that infection by a recombinant between different MSV isolates led to increased symptom severity in maize, indicating greater adaptation of the recombinant. On the other hand, Davino et al.  demonstrated that a recombinant between TYLCV and TYLCSV coexists with the parental viruses and has a lower replicative capacity, a fact that might help to explain why plants infected only by the recombinant have not been found. In the case of the recombinant ToRMV, there was an adaptive advantage over the parental ToSRV reflected in the preferential replication of ToRMV over ToSRV, nevertheless ToSRV is more prevalent in the field compared to ToRMV. We suggest that ToSRV is maintained in nature by infecting plants in the absence of ToRMV. If the occurrence of mixed infection between ToRMV and ToSRV was frequent, ToRMV would be the prevalent virus.
Another interesting aspect is that the two DNA-B components are almost identical (98.2% nt identity). This supports the hypothesis that, following the transfer of practically the entire Rep ORF as well as part of the CR from ToSRV to ToRMV, ToRMV captured the DNA-B of ToSRV. According to the current ICTV-supported criteria for species demarcation, a begomovirus is considered to represent a new species when its DNA-A sequence is less than 89% identical to that of other begomoviruses . ToSRV and ToRMV DNA-A components have nucleotide sequence identity of 86.3%, which clearly separates them as members of different species based on this criterion. However, these two viruses are very similar in terms of biological and molecular properties, including: (i) high sequence identity of their CRs and their Rep coding regions; (ii) nearly identical Rep binding sites in their CRs, with identical iterons; (iii) the same DNA-B; and (iv) induction of nearly identical symptoms in tomato, with the same latent period. We thus propose that ToRMV and ToSRV should be considered to be strains of the same viral species, named ToRMV since this virus was described first.
The fact that begomovirus taxonomy ignores the DNA-B was addressed by Briddon et al. , who performed a hypothetical taxonomy based on the DNA-B sequence alone. The authors show that 85% of the sequences analyzed remain unaltered in their classification based on the DNA-A, providing good support for a taxonomy based only on the DNA-A. However, in specific cases such as the one presented in this work, the criterion based on the DNA-A could be less restrictive and/or the DNA-B could be considered in the analysis for species assignment.
The recombination analysis together with the formation of viable pseudorecombinants between ToSRV and ToRMV support the hypothesis that ToRMV captured the DNA-B of ToSRV after acquiring, by recombination, the DNA-A region that includes the Rep ORF and the CR. Therefore, the emergence of ToRMV involves both recombination and pseudorecombination events with ToSRV. We propose that ToRMV and ToSRV should be considered to be members of the same viral species, retaining the name ToRMV since it was first described.
Pseudorecombination between ToRMV and ToSRV
Full-length DNA-A and DNA-B clones of the viral isolates ToRMV-[BR:Ub1:96]  and ToSRV-[BR:PG1:Pep:03]  were used in the experiments. Tomato (S. lycopersicum cv. ‘Santa Clara’) seedlings were biolistically inoculated  using two micrograms of each genomic component (DNA-A and DNA-B) in every possible combination of single or mixed infection. Inoculated plants were kept in a greenhouse with average daily temperatures of 26 ± 2°C. Two independent experiments were carried out, with six plants inoculated for each treatment in the first experiment, and seven in the second experiment.
Inoculated plants were evaluated for symptom expression for up to 35 days post-inoculation (dpi). DNA from all plants was extracted at 28 dpi as described  and the presence of each viral genomic component was examined by rolling circle amplification (RCA) , followed by digestion with component specific restriction enzymes: XhoI and BglII to detect ToRMV-A and -B, respectively; HindIII and BglII + SacII to detect ToSRV-A and -B, respectively.
Quantification of ToRMV and ToSRV by quantitative, real time PCR (qPCR)
Based on the results obtained from RCA-RFLP assays, control plants (mock-inoculated, inoculated with ToRMV DNA-A + DNA-B and inoculated with ToSRV DNA-A + DNA-B) and plants in which the DNA-A’s of both viruses were co-inoculated in different combinations (ToRMV DNA-A + DNA-B and ToSRV DNA-A + DNA-B; ToRMV DNA-A + DNA-B and ToSRV DNA-A; ToRMV DNA-A and ToSRV DNA-A + DNA-B) were selected for quantitative, real-time PCR (qPCR) assays. Primers were designed based on the sequence of the common region of the DNA-A (in the most divergent region between the two viruses) using the qPCR tool implemented in http://www.idtdna.com/scitools/Applications/RealTimePCR: ToRMV-A(For), 5′CAT CGG GCC TCT GTT GG3′ and ToRMV-A (Rev), 5′GTT ATG CAA CTT GGG CGT TAA G3′; ToSRV-A(For), 5′AAA GTA AAG TGA TTG TCT GTG G3′ and ToSRV-A(Rev), 5′GCC GTT CAA CAA ATT GGG3′. Primer specificity was tested by conventional PCR using positive samples for each virus, followed by electrophoresis in 2.0% agarose gel. Reactions were prepared in a final volume of 10 μL, using Fast SYBR Green Master Mix (Applied Biosystems) following the manufacturer’s instructions, and were analyzed in a CFX96 Real-Time PCR System (Bio-Rad). The PCR protocol included an initial denaturing step at 95°C for 20s, followed by 39 cycles of 95°C for 3 s and 60°C for 30s, followed by a dissociation stage.
Standard curves were prepared using serial dilutions of plasmid DNA containing an insert corresponding to the complete ToSRV DNA-A (10° to 106 copies of the viral genome per reaction). Standard curves were obtained by regression analysis of cycle threshold (Ct) values of each one of the three replications of a given dilution in relation to the log of the amount of DNA in each dilution. For absolute quantification of the number of viral DNA molecules in the different treatments, 100 ng of total DNA was extracted as described previously and used in reactions containing virus-specific primers. Each sample was analyzed in triplicates, and three different plants of each treatment described above were assayed. DNA-B accumulation was not evaluated due to the high sequence similarity between the two viruses, which prevented the design of virus-specific primers for the qPCR assay.
Sequence comparisons and recombination analysis
Complete genomes, individual coding regions and the CRs of ToRMV and ToSRV were compared using EMBOSS pairwise alignment algorithms (http://www.ebi.ac.uk/Tools/emboss/align/index.html), using the EMBOSS needle (Global) with default settings. Nineteen New World begomovirus DNA components (Table 3) that occur in weeds and tomato plants in Brazil were aligned for analysis of recombination events involving ToRMV and ToSRV, using RDP3 . Recombination events detected by at least four of the analysis methods available in the program were considered reliable. Alignments were scanned using default settings for each analysis method using a Bonferroni-corrected P value cutoff of 0.05.
Multiple sequence alignments were performed for the DNA-A sequences of Brazilian begomoviruses (Table 3) based on the recombinant and non-recombinant regions of the ToRMV genome (determined in RDP3) and for the cognate full-length DNA-B sequences using Muscle . Maximum likelihood (ML) trees were inferred using PAUP v. 4.0  with the best fit model of nucleotide substitution determined using Modeltest  by the Akaike Information Criterion (AIC). A heuristic search was initiated with a neighbor-joining tree using the tree-bisection-reconnection (TBR) algorithm to optimize the ML tree. The robustness of each internal branch was estimated from 3,000 bootstrap replications. The Nearest Neighbor Interchange (NNI) algorithm was used to optimize the bootstrap replications. Trees were visualized and edited using the FigTree (http://www.tree.bio.ed.ac.uk/software/figtree). The begomovirus Macroptilium mosaic Puerto Rico virus (MaMPRV) [GenBank: AY044133] was used as an outgroup.
Brown JK, Fauquet CM, Briddon RW, Zerbini FM, Moriones E, Navas-Castillo J: Family Geminiviridae. Virus Taxonomy 9th Report of the International committee on taxonomy of viruses. Edited by: King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. 2012, London, UK: Elsevier Academic Press, 351-373.
Briddon RW, Patil BL, Bagewadi B, Nawaz-ul-Rehman MS, Fauquet CM: Distinct evolutionary histories of the DNA-A and DNA-B components of bipartite begomoviruses. BMC Evol Biol. 2010, 10: 97-10.1186/1471-2148-10-97.
Rojas MR, Hagen C, Lucas WJ, Gilbertson RL: Exploiting chinks in the plant’s armor: Evolution and emergence of geminiviruses. Annu Rev Phytopathol. 2005, 43: 361-394. 10.1146/annurev.phyto.43.040204.135939.
Fontes EPB, Eagle PA, Sipe PS, Luckow VA, Hanley-Bowdoin L: Interaction between a geminivirus replication protein and origin DNA is essential for viral replication. J Biol Chem. 1994, 269: 8459-8465.
Chatterji A, Padidam M, Beachy RN, Fauquet CM: Identification of replication specificity determinants in two strains of tomato leaf curl virus from New Delhi. J Virol. 1999, 73: 5481-5489.
Arguello-Astorga G, Herrera-Estrella L, Rivera-Bustamante R: Experimental and theoretical definition of geminivirus origin of replication. Plant Mole Biol. 1994, 26: 553-556. 10.1007/BF00013742.
Pita JS, Fondong VN, Sangare 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.
Faria JC, Gilbertson RL, Hanson SF, Morales FJ, Ahlquist PG, Loniello AO, Maxwell DP: Bean golden mosaic geminivirus type II isolates from the Dominican Republic and Guatemala: Nucleotide sequences, infectious pseudorecombinants, and phylogenetic relationships. Phytopathol. 1994, 84: 321-329. 10.1094/Phyto-84-321.
Andrade EC, Manhani GG, Alfenas PF, Calegario RF, Fontes EPB, Zerbini FM: Tomato yellow spot virus, a tomato-infecting begomovirus from Brazil with a closer relationship to viruses from Sida sp., forms pseudorecombinants with begomoviruses from tomato but not from Sida. J Gen Virol. 2006, 87: 3687-3696. 10.1099/vir.0.82279-0.
Gilbertson RL, Hidayat SH, Paplomatas EJ, Rojas MR, Hou Y-H, Maxwell DP: Pseudorecombination between infectious cloned DNA components of tomato mottle and bean dwarf mosaic geminiviruses. J Gen Virol. 1993, 74: 23-31. 10.1099/0022-1317-74-1-23.
Garrido-Ramirez ER, Sudarshana M, Gilbertson RL: Bean golden yellow mosaic virus from Chiapas, Mexico: characterization, pseudorecombination with other bean-infecting geminiviruses and germ plasm screening. Phytopathol. 2000, 90: 1224-1232. 10.1094/PHYTO.2000.90.11.1224.
Bull SE, Briddon RW, Sserubombwe WS, Ngugi K, Markham PG, Stanley J: Infectivity, pseudorecombination and mutagenesis of Kenyan cassava mosaic begomoviruses. J Gen Virol. 2007, 88: 1624-1633. 10.1099/vir.0.82662-0.
Ramos PL, Guevara-Gonzalez RG, Peral R, Ascencio-Ibanez JT, Polston JE, Arguello-Astorga GR, Vega-Arreguin JC, Rivera-Bustamante RF: Tomato mottle Taino virus pseudorecombines with PYMV but not with ToMoV: Implications for the delimitation of cis- and trans-acting replication specificity determinants. Arch Virol. 2003, 148: 1697-1712. 10.1007/s00705-003-0136-3.
Chakraborty S, Vanitharani R, Chattopadhyay B, Fauquet CM: Supervirulent pseudorecombination and asymmetric synergism between genomic components of two distinct species of begomovirus associated with severe tomato leaf curl disease in India. J Gen Virol. 2008, 89: 818-828. 10.1099/vir.0.82873-0.
Arguello-Astorga GR, Ruiz-Medrano R: An iteron-related domain is associated to Motif 1 in the replication proteins of geminiviruses: Identification of potential interacting amino acid-base pairs by a comparative approach. Arch Virol. 2001, 146: 1465-1485. 10.1007/s007050170072.
Orozco BM, Hanley-Bowdoin L: Conserved sequence and structural motifs contribute to the DNA binding and cleavage activities of a geminivirus replication protein. J Biol Chem. 1998, 273: 24448-24456. 10.1074/jbc.273.38.24448.
Orozco BM, Miller AB, Settlage SB, Hanley-Bowdoin L: Functional domains of a geminivirus replication protein. J Biol Chem. 1997, 272: 9840-9846. 10.1074/jbc.272.15.9840.
García-Arenal F, Fraile A, Malpica JM: Variability and genetic structure of plant virus populations. Annu Rev Phytopathol. 2001, 39: 157-186. 10.1146/annurev.phyto.39.1.157.
Lima ATM, Sobrinho RR, Gonzalez-Aguilera J, Rocha CS, Silva SJC, Xavier CAD, Silva FN, Duffy S, Zerbini FM: Synonymous site variation due to recombination explains higher genetic variability in begomovirus populations infecting non-cultivated hosts. J Gen Virol. 2013, 94: 418-431. 10.1099/vir.0.047241-0.
Tiendrebeogo F, Lefeuvre P, Hoareau M, Harimalala MA, De Bruyn A, Villemot J, Traore VS, Konate G, Traore AS, Barro N, Reynaud B, Traoré O, Lett J-M: Evolution of African cassava mosaic virus by recombination between bipartite and monopartite begomoviruses. Virol J. 2012, 9: 67-10.1186/1743-422X-9-67.
Davino S, Napoli C, Dellacroce C, Miozzi L, Noris E, Davino M, Accotto GP: Two new natural begomovirus recombinants associated with the tomato yellow leaf curl disease co-exist with parental viruses in tomato epidemics in Italy. Virus Res. 2009, 143: 15-23. 10.1016/j.virusres.2009.03.001.
Monci F, Sanchez-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.
Patil BL, Fauquet CM: Cassava mosaic geminiviruses: actual knowledge and perspectives. Mole Plant Pathol. 2009, 10: 685-701. 10.1111/j.1364-3703.2009.00559.x.
van der Walt E, Rybicki EP, Varsani A, Polston JE, Billharz R, Donaldson L, Monjane AL, Martin DP: Rapid host adaptation by extensive recombination. J Gen Virol. 2009, 90: 734-746. 10.1099/vir.0.007724-0.
Harrison BD, Zhou X, Otim Nape GW, Liu Y, Robinson DJ: Role of a novel type of double infection in the geminivirus-induced epidemic of severe cassava mosaic in Uganda. Ann Appl Biol. 1997, 131: 437-448. 10.1111/j.1744-7348.1997.tb05171.x.
Zhou X, Liu Y, Calvert L, Munoz C, Otim-Nape GW, Robinson DJ, Harrison BD: Evidence that DNA-A of a geminivirus associated with severe cassava mosaic disease in Uganda has arisen by interspecific recombination. J Gen Virol. 1997, 78: 2101-2111.
Rocha CS, Castillo-Urquiza GP, Lima ATM, Silva FN, Xavier CAD, Hora-Junior BT, Beserra-Junior JEA, Malta AWO, Martin DP, Varsani A, Alfenas-Zerbini P, Mizubuti ESG, Zerbini FM: Brazilian begomovirus populations are highly recombinant, rapidly evolving, and segregated based on geographical location. J Virol. 2013, 87: 5784-5799. 10.1128/JVI.00155-13.
Hou YM, Gilbertson RL: Increased pathogenicity in a pseudorecombinant bipartite geminivirus correlates with intermolecular recombination. J Virol. 1996, 70: 5430-5436.
Zerbini FM, Andrade EC, Barros DR, Ferreira SS, Lima ATM, Alfenas PF, Mello RN: Traditional and novel strategies for geminivirus management in Brazil. Australas Plant Pathol. 2005, 34: 475-480. 10.1071/AP05088.
Ambrozevicius LP, Calegario RF, Fontes EPB, Carvalho MG, Zerbini FM: Genetic diversity of begomoviruses infecting tomato and associated weeds in Southeastern Brazil. Braz Phytopathol. 2002, 27: 372-377.
Ribeiro SG, Ambrozevicius LP, Ávila AC, Bezerra IC, Calegario RF, Fernandes JJ, Lima MF, Mello RN, Rocha H, Zerbini FM: Distribution and genetic diversity of tomato-infecting begomoviruses in Brazil. Arch Virol. 2003, 148: 281-295. 10.1007/s00705-002-0917-0.
Fernandes JJ, Carvalho MG, Andrade EC, Brommonschenkel SH, Fontes EPB, Zerbini FM: Biological and molecular properties of Tomato rugose mosaic virus (ToRMV), a new tomato-infecting begomovirus from Brazil. Plant Pathol. 2006, 55: 513-522. 10.1111/j.1365-3059.2006.01395.x.
Ribeiro SG, Martin DP, Lacorte C, Simões IC, Orlandini DRS, Inoue-Nagata AK: Molecular and biological characterization of Tomato chlorotic mottle virus suggests that recombination underlies the evolution and diversity of Brazilian tomato begomoviruses. Phytopathology. 2007, 97: 702-711. 10.1094/PHYTO-97-6-0702.
Calegario RF, Ferreira SS, Andrade EC, Zerbini FM: Characterization of Tomato yellow spot virus, (ToYSV), a novel tomato-infecting begomovirus from Brazil. Braz J Agric Res. 2007, 42: 1335-1343.
Fernandes FR, Albuquerque LC, Giordano LB, Boiteux LS, Ávila AC, Inoue-Nagata AK: Diversity and prevalence of Brazilian bipartite begomovirus species associated to tomatoes. Virus Genes. 2008, 36: 251-258. 10.1007/s11262-007-0184-y.
Albuquerque LC, Varsani A, Fernandes FR, Pinheiro B, Martin DP, Oliveira Ferreira PT, Lemos TO, Inoue-Nagata AK: Further characterization of tomato-infecting begomoviruses in Brazil. Arch Virol. 2012, 157: 747-752. 10.1007/s00705-011-1213-7.
Londono A, Riego-Ruiz L, Arguello-Astorga GR: DNA-binding specificity determinants of replication proteins encoded by eukaryotic ssDNA viruses are adjacent to widely separated RCR conserved motifs. Arch Virol. 2010, 155: 1033-1046. 10.1007/s00705-010-0674-4.
Galvão RM, Mariano AC, Luz DF, Alfenas PF, Andrade EC, Zerbini FM, Almeida MR, Fontes EPB: A naturally occurring recombinant DNA-A of a typical bipartite begomovirus does not require the cognate DNA-B to infect Nicotiana benthamiana systemically. J Gen Virol. 2003, 84: 715-726. 10.1099/vir.0.18783-0.
Alves-Junior M, Alfenas-Zerbini P, Andrade EC, Esposito DA, Silva FN, Cruz ACF, Ventrella MC, Otoni WC, Zerbini FM: Synergism and negative interference during co-infection of tomato and Nicotiana benthamiana with two bipartite begomoviruses. Virology. 2009, 387: 257-266. 10.1016/j.virol.2009.01.046.
Choi IR, Stenger DC: Strain-specific determinants of beet curly top geminivirus DNA replication. Virology. 1995, 206: 904-912. 10.1006/viro.1995.1013.
Gladfelter HJ, Eagle PA, Fontes EPB, Batts L, Hanley-Bowdoin L: Two domains of the AL 1 protein mediate geminivirus origin recognition. Virology. 1997, 239: 186-197. 10.1006/viro.1997.8869.
Ilyina TV, Koonin EV: Conserved sequence motifs in the initiator proteins for rolling circle DNA replication encoded by diverse replicons from eubacteria, eucaryotes and archaebacteria. Nucleic Acids Res. 1992, 20: 3279-3285. 10.1093/nar/20.13.3279.
Stanley J: Analysis of African cassava mosaic virus recombinants suggest strand nicking occurs within the conserved nonanucleotide motif during the initiation of rolling circle DNA replication. Virology. 1995, 206: 707-712. 10.1016/S0042-6822(95)80093-X.
Sanz AI, Fraile A, Gallego JM, Malpica JM, García-Arenal F: Genetic variability of natural populations of cotton leaf curl geminivirus, a single-stranded DNA virus. J Mole Evol. 1999, 49: 672-681. 10.1007/PL00006588.
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.
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 J-M, 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.
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.
Vuillaume F, Thebaud G, Urbino C, Forfert N, Granier M, Froissart R, Blanc S, Peterschmitt M: Distribution of the phenotypic effects of random homologous recombination between two virus species. PLoS Pathog. 2011, 7: e1002028-10.1371/journal.ppat.1002028.
Padidam M, Sawyer S, Fauquet CM: Possible emergence of new geminiviruses by frequent recombination. Virology. 1999, 265: 218-224. 10.1006/viro.1999.0056.
Nawaz-Ul-Rehman MS, Fauquet CM: Evolution of geminiviruses and their satellites. FEBS Lett. 2009, 583: 1825-1832. 10.1016/j.febslet.2009.05.045.
Fauquet CM, Briddon RW, Brown JK, Moriones E, Stanley J, Zerbini FM, Zhou X: Geminivirus strain demarcation and nomenclature. Arch Virol. 2008, 153: 783-821. 10.1007/s00705-008-0037-6.
Bezerra-Agasie IC, Ferreira GB, Ávila AC, Inoue-Nagata AK: First report of Tomato severe rugose virus in chili pepper in Brazil. Plant Dis. 2006, 90:
Aragão FJL, Barros LMG, Brasileiro ACM, Ribeiro SG, Smith FD, Sanford JC, Faria JC, Rech EL: Inheritance of foreign genes in transgenic bean (Phaseolus vulgaris L.) co-transformed via particle bombardment. Theor Appl Genet. 1996, 93: 142-150. 10.1007/BF00225739.
Dellaporta SL, Wood J, Hicks JB: A plant DNA minipreparation: Version II. Plant Mol Biol Rep. 1983, 1: 19-21. 10.1007/BF02712670.
Inoue-Nagata AK, Albuquerque LC, Rocha WB, Nagata T: A simple method for cloning the complete begomovirus genome using the bacteriophage phi 29 DNA polymerase. J Virol Methods. 2004, 116: 209-211. 10.1016/j.jviromet.2003.11.015.
Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P: RDP3: A flexible and fast computer program for analyzing recombination. Bioinformatics. 2010, 26: 2462-2463. 10.1093/bioinformatics/btq467.
Edgar RC: MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5: 1-19. 10.1186/1471-2105-5-1.
Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. 2003, Sunderland, Massachusetts: Sinauer Associates
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
We thank Alice K. Inoue-Nagata for providing the ToSRV clones. This work was funded by CNPq grant 484447/2011-4 and FAPEMIG grants CAG-838-04 and APQ-949-09 to FMZ. GPCU was a CAPES-PNPD postdoctoral fellow. ATML is a FAPEMIG postdoctoral fellow.
The authors declare that they have no competing interests.
FMZ, FNS and ATML designed the study. FNS, ATML, GPCU and MAJ performed the experiments. FNS, ATML and CSR analyzed the data. FNS and FMZ wrote the manuscript. All authors have read and approved the final version of the manuscript.
Fábio N Silva, Alison TM Lima contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1: Detection by rolling circle amplification followed by restriction digests (RCA-RFLP) of each DNA component from ToRMV-[BR:Ub1:96] and ToSRV-[BR:PG1:Pep:03] after biolistic inoculation of tomato plants. A. Digestion with Hind III, generating a fragment of approximately 2600 nucleotides (nt) corresponding to ToSRV DNA-A (the other three components are not cleaved by this enzyme). B. Digestion with Xho I generating a fragment of approximately 2600 nt corresponding to ToRMV DNA-A (the other three components are not cleaved by this enzyme). C. Digestion with Bgl II + Sac II, generating fragments of 1800 and 800 nt corresponding to ToSRV DNA-B, and approximately 2600 nt corresponding to ToRMV DNA-B (the two DNA-A components are not cleaved by these enzymes). Underlined means DNA-A + DNA-B. (–) negative control. M, 1 kb plus DNA ladder (Invitrogen), in bp. (TIFF 474 KB)
Additional file 2: Figure S2: ToSRV-[BR:PG1:Pep:03] DNA-A accumulation in single or dual infections with ToRMV-[BR:Ub1:96] in tomato plants. Total DNA was extracted from systemically infected leaves at 28 days post-inoculation and used as a template for quantitative, real-time PCR (qPCR) using ToSRV DNA-A primers. 1. Mock-inoculated plants; 2. ToRMV DNA-A + DNA-B; 3. ToSRV DNA-A + DNA-B; 4. ToRMV DNA-A + DNA-B and ToSRV DNA-A + DNA-B; 5. ToRMV DNA-A + DNA-B and ToSRV DNA-A; 6. ToRMV DNA-A and ToSRV DNA-A + DNA-B. In treatment 2, the ToSRV primers weakly amplified ToRMV DNA-A. Values are given as the mean ± confidence interval of three independent tomato plants. Asterisks indicate statistically significant differences by the t test at P ≤ 0.05. (JPEG 483 KB)
About this article
Cite this article
Silva, F.N., Lima, A.T., Rocha, C.S. et al. Recombination and pseudorecombination driving the evolution of the begomoviruses Tomato severe rugose virus (ToSRV) and Tomato rugose mosaic virus (ToRMV): two recombinant DNA-A components sharing the same DNA-B. Virol J 11, 66 (2014). https://doi.org/10.1186/1743-422X-11-66