Transcriptome and small RNAome profiling uncovers how a recombinant begomovirus evades RDRγ-mediated silencing of viral genes and outcompetes its parental virus in mixed infection.

Tomato yellow leaf curl virus (TYLCV, genus Begomovirus, family Geminiviridae) causes severe disease of cultivated tomatoes. Geminiviruses replicate circular single-stranded genomic DNA via rolling-circle and recombination-dependent mechanisms, frequently generating recombinants in mixed infections. Circular double-stranded intermediates of replication also serve as templates for Pol II bidirectional transcription. IS76, a recombinant derivative of TYLCV with a short sequence in the bidirectional promoter/origin-of-replication region acquired from a related begomovirus, outcompetes TYLCV in mixed infection and breaks disease resistance in tomato Ty-1 cultivars. Ty-1 encodes a γ-clade RNA-dependent RNA polymerase (RDRγ) implicated in Dicer-like (DCL)-mediated biogenesis of small interfering (si)RNAs directing gene silencing. Here, we profiled transcriptome and small RNAome of Ty-1 resistant and control susceptible plants infected with TYLCV, IS76 or their combination at early and late infection stages. We found that RDRγ boosts production rates of 21, 22 and 24 nt siRNAs from entire genomes of both viruses and modulates DCL activities in favour of 22 and 24 nt siRNAs. Compared to parental TYLCV, IS76 undergoes faster transition to the infection stage favouring rightward transcription of silencing suppressor and coat protein genes, thereby evading RDRγ activity and facilitating its DNA accumulation in both single and mixed infections. In coinfected Ty-1 plants, IS76 efficiently competes for host replication and transcription machineries, thereby impairing TYLCV replication and transcription and forcing its elimination associated with further increased siRNA production. RDRγ is constitutively overexpressed in Ty-1 plants, which correlates with begomovirus resistance, while siRNA-generating DCLs (DCL2b/d, DCL3, DCL4) and genes implicated in siRNA amplification (α-clade RDR1) and function (Argonaute2) are upregulated to similar levels in TYLCV- and IS76-infected susceptible plants. Collectively, IS76 recombination facilitates replication and promotes expression of silencing suppressor and coat proteins, which allows the recombinant virus to evade the negative impact of RDRγ-boosted production of viral siRNAs directing transcriptional and posttranscriptional silencing.

Refining the emergence scenario of the invasive recombinant Tomato yellow leaf curl virus -IS76.

TYLCV-IS76, a unique recombinant between tomato yellow leaf curl virus (TYLCV) and tomato yellow leaf curl Sardinia virus (TYLCSV), has replaced its parental viruses in southern Morocco. To refine its emergence scenario, its fitness was monitored experimentally in conditions aiming at reproducing natural situations, i.e. superinfection of plants already infected with parental viruses and competition with other TYLCV/TYLCSV recombinants (LSRec) automatically generated in plants coinfected with TYLCV and TYLCSV. TYLCV-IS76 accumulated significantly more than parental viruses regardless of plant age and superinfection delay. Although TYLCV-IS76 and LSRec both accumulated more than parental viruses in laboratory conditions, LSRec were displaced by TYLCV-IS76 in nature like parental viruses were. TYLCV-IS76 did not exhibit any vector transmission advantage over LSRec and TYLCV the most competitive parental virus. Thus, it is apparently only in the plant compartment that the recombination event that generated TYLCV-IS76, induced the competitiveness advantage by which the last became first.

Synergy between an emerging monopartite begomovirus and a DNA-B component.

In recent decades, a legion of monopartite begomoviruses transmitted by the whitefly Bemisia tabaci has emerged as serious threats to vegetable crops in Africa. Recent studies in Burkina Faso (West Africa) reported the predominance of pepper yellow vein Mali virus (PepYVMLV) and its frequent association with a previously unknown DNA-B component. To understand the role of this DNA-B component in the emergence of PepYVMLV, we assessed biological traits related to virulence, virus accumulation, location in the tissue and transmission. We demonstrate that the DNA-B component is not required for systemic movement and symptom development of PepYVMLV (non-strict association), but that its association produces more severe symptoms including growth arrest and plant death. The increased virulence is associated with a higher viral DNA accumulation in plant tissues, an increase in the number of contaminated nuclei of the phloem parenchyma and in the transmission rate by B. tabaci. Our results suggest that the association of a DNA-B component with the otherwise monopartite PepYVMLV is a key factor of its emergence.

Fitness advantage of inter-species TYLCV recombinants induced by beneficial intra-genomic interactions rather than by specific mutations.

Tomato yellow leaf curl virus (TYLCV) and its related viruses are prone to recombination. It was reported that random homologous recombination between 20% diverging TYLCV related species is rarely deleterious and may be associated with a fitness advantage. Indeed, TYLCV-IS76, a recombinant between the 20% divergent TYLCV and tomato yellow leaf curl Sardinia virus (TYLCSV), exhibited a higher fitness than that of parental viruses. As this typical fitness advantage was observed with TYLCV-IS76 representatives of different pedigrees, it was thought that it is induced by beneficial intra-genomic interactions rather than by specific mutations. This hypothesis was further supported with TYLCV-IS141, a TYLCV recombinant with a short TYLCSV inherited fragment of around 141 nts, slightly longer than that of TYLCV-IS76. Indeed, the typical fitness advantage was detected irrespective of the position of the recombination breakpoint (loci 76 or 141) and the sequences of the TYLCV and TYLCSV inherited fragments.

A multicellular way of life for a multipartite virus.

A founding paradigm in virology is that the spatial unit of the viral replication cycle is an individual cell. Multipartite viruses have a segmented genome where each segment is encapsidated separately. In this situation the viral genome is not recapitulated in a single virus particle but in the viral population. How multipartite viruses manage to efficiently infect individual cells with all segments, thus with the whole genome information, is a long-standing but perhaps deceptive mystery. By localizing and quantifying the genome segments of a nanovirus in host plant tissues we show that they rarely co-occur within individual cells. We further demonstrate that distinct segments accumulate independently in different cells and that the viral system is functional through complementation across cells. Our observation deviates from the classical conceptual framework in virology and opens an alternative possibility (at least for nanoviruses) where the infection can operate at a level above the individual cell level, defining a viral multicellular way of life.

Characterization of Begomoviruses Sampled during Severe Epidemics in Tomato Cultivars Carrying the Ty-1 Gene.

Tomato yellow leaf curl virus (TYLCV, genus Begomovirus, family Geminiviridae) is a major species that causes a tomato disease for which resistant tomato hybrids (mainly carriers of the Ty-1/Ty-3 gene) are being used widely. We have characterized begomoviruses severely affecting resistant tomato crops in Southeast Spain. Circular DNA was prepared from samples by rolling circle amplification, and sequenced by massive sequencing (2015) or cloning and Sanger sequencing (2016). Thus, 23 complete sequences were determined, all belonging to the TYLCV Israel strain (TYLCV-IL). Massive sequencing also revealed the absence of other geminiviral and beta-satellite sequences. A phylogenetic analysis showed that the Spanish isolates belonged to two groups, one related to early TYLCV-IL isolates in the area (Group 1), and another (Group 2) closely related to El Jadida (Morocco) isolates, suggesting a recent introduction. The most parsimonious evolutionary scenario suggested that the TYLCV isolates of Group 2 are back recombinant isolates derived from TYLCV-IS76, a recombinant virus currently predominating in Moroccan epidemics. Thus, an infectious Group 2 clone (TYLCV-Mu15) was constructed and used in in planta competition assays against TYLCV-IS76. TYLCV-Mu15 predominated in single infections, whereas TYLCV-IS76 did so in mixed infections, providing credibility to a scenario of co-occurrence of both types of isolates.

Accumulation and transmission of alphasatellite, betasatellite and tomato yellow leaf curl virus in susceptible and Ty-1-resistant tomato plants.

Begomoviruses (family Geminiviridae) are frequently associated with alphasatellites and betasatellites in the Old World. Tomato yellow leaf curl virus, one of the most damaging begomovirus species worldwide, was recently found associated with betasatellites in the eastern coast of the Mediterranean Sea, and in the Middle East region. Tomato yellow leaf curl virus (TYLCV)/betasatellite associations were shown to increase TYLCV virulence in experimental conditions. The sustainability of TYLCV/satellite associations in tomato was assessed here by estimating accumulation levels of satellites in comparison to TYLCV, vector transmission efficiency, and by testing how far the popular Ty-1 resistance gene used in most TYLCV-resistant tomato cultivars in the Mediterranean Basin is effective against betasatellites. Three satellites previously isolated from okra in Burkina Faso-of the species Cotton leaf curl Gezira betasatellite, Cotton leaf curl Gezira alphasatellite and Okra leaf curl Burkina Faso alphasatellite-were shown to accumulate at levels similar to, or higher than, the helper virus TYLCV-Mld in tomato plants from 32 to 150 days post inoculation (dpi). Cotton leaf curl Gezira betasatellite (CLCuGB) reduced TYLCV-Mld accumulation whereas alphasatellites did not. Transmission tests were performed with B. tabaci from plants infected with TYLCV-Mld/CLCuGB- or TYLCV-Mld/Okra leaf curl Burkina Faso alphasatellite. At 32 dpi, both satellites were transmitted to more than 50% of TYLCV-infected test plants. Betasatellite transmission, tested further with 150 dpi source plants was successful in more than 30% of TYLCV-infected test plants. Ty-1 resistant tomato plants co-infected with TYLCV (-Mld or -IL) and CLCuGB exhibited mild leaf curling and mosaic symptoms at the early stage of infection associated with a positive effect on TYLCV-IL accumulation, while resistant plants infected with TYLCV only, were asymptomatic. Together with previous experimental studies, these results further emphasize the potential risk of betasatellites to tomato cultivation, including with Ty-1 resistant cultivars.

The non-canonical tomato yellow leaf curl virus recombinant that displaced its parental viruses in southern Morocco exhibits a high selective advantage in experimental conditions.

Recombination events are frequently inferred from the increasing number of sequenced viral genomes, but their impact on natural viral populations has rarely been evidenced. TYLCV-IS76 is a recombinant (Begomovirus,Geminiviridae) between the Israel strain of tomato yellow leaf curl virus (TYLCV-IL) and the Spanish strain of tomato yellow leaf curl Sardinia virus (TYLCSV-ES) that was generated most probably in the late 1990s in southern Morocco (Souss). Its emergence in the 2000s coincided with the increasing use of resistant tomato cultivars bearing the Ty-1 gene, and led eventually to the entire displacement of both parental viruses in the Souss. Here, we provide compelling evidence that this viral population shift was associated with selection of TYLCV-IS76 viruses in tomato plants and particularly in Ty-1-bearing cultivars. Real-time quantitative PCR (qPCR) monitoring revealed that TYLCV-IS76 DNA accumulation in Ty-1-bearing plants was significantly higher than that of representatives of the parental virus species in single infection or competition assays. This advantage of the recombinant in Ty-1-bearing plants was not associated with a fitness cost in a susceptible, nearly isogenic, cultivar. In competition assays in the resistant cultivar, the DNA accumulation of the TYLCV-IL clone - the parent less affected by the Ty-1 gene in single infection - dropped below the qPCR detection level at 120 days post-infection (p.i.) and below the whitefly vector (Bemisia tabaci) transmissibility level at 60 days p.i. The molecular basis of the selective advantage of TYLCV-IS76 is discussed in relation to its non-canonical recombination pattern, and the RNA-dependent RNA polymerase encoded by the Ty-1 gene.

Mapping genetic determinants of viral traits with FST and quantitative trait locus (QTL) approaches.

The genetic determinism of viral traits can generally be dissected using either forward or reverse genetics because the clonal reproduction of viruses does not require the use of approaches based on laboratory crosses. Nevertheless, we hypothesized that recombinant viruses could be analyzed as sexually reproducing organisms, using either a quantitative trait loci (QTL) approach or a locus-by-locus fixation index (FST). Locus-by-locus FST analysis, and four different regressions and interval mapping algorithms of QTL analysis were applied to a phenotypic and genotypic dataset previously obtained from 47 artificial recombinant genomes generated between two begomovirus species. Both approaches assigned the determinant of within-host accumulation-previously identified using standard virology approaches-to a region including the 5׳ end of the replication-associated protein (Rep) gene and the upstream intergenic region. This study provides a proof of principle that QTL and population genetics tools can be extended to characterize the genetic determinants of viral traits.

Within-host dynamics of the emergence of Tomato yellow leaf curl virus recombinants.

Tomato yellow leaf curl virus (TYLCV) is a highly damaging begomovirus native to the Middle East. TYLCV has recently spread worldwide, recombining with other begomoviruses. Recent analysis of mixed infections between TYLCV and Tomato leaf curl Comoros begomovirus (ToLCKMV) has shown that, although natural selection preserves certain co-evolved intra-genomic interactions, numerous and diverse recombinants are produced at 120 days post-inoculation (dpi), and recombinant populations from different tomato plants are very divergent. Here, we investigate the population dynamics that lead to such patterns in tomato plants co-infected with TYLCV and ToLCKMV either by agro-inoculation or using the natural whitefly vector Bemisia tabaci. We monitored the frequency of parental and recombinant genotypes independently in 35 plants between 18 and 330 dpi and identified 177 recombinants isolated at different times. Recombinants were detected from 18 dpi and their frequency increased over time to reach about 50% at 150 dpi regardless of the inoculation method. The distribution of breakpoints detected on 96 fully sequenced recombinants was consistent with a continuous generation of new recombinants as well as random and deterministic effects in their maintenance. A severe population bottleneck of around 10 genomes was estimated during early systemic infection-a phenomenon that could account partially for the heterogeneity in recombinant patterns observed among plants. The detection of the same recombinant genome in six of the thirteen plants analysed beyond 30 dpi supported the influence of selection on observed recombination patterns. Moreover, a highly virulent recombinant genotype dominating virus populations within one plant has, apparently, the potential to be maintained in the natural population according to its infectivity, within-host accumulation, and transmission efficiency - all of which were similar or intermediate to those of the parent genotypes. Our results anticipate the outcomes of natural encounters between TYLCV and ToLCKMV.

Distribution of the phenotypic effects of random homologous recombination between two virus species.

Recombination has an evident impact on virus evolution and emergence of new pathotypes, and has generated an immense literature. However, the distribution of phenotypic effects caused by genome-wide random homologous recombination has never been formally investigated. Previous data on the subject have promoted the implicit view that most viral recombinant genomes are likely to be deleterious or lethal if the nucleotide identity of parental sequences is below 90%. We decided to challenge this view by creating a bank of near-random recombinants between two viral species of the genus Begomovirus (Family Geminiviridae) exhibiting 82% nucleotide identity, and by testing infectivity and in planta accumulation of recombinant clones randomly extracted from this bank. The bank was created by DNA-shuffling-a technology initially applied to the random shuffling of individual genes, and here implemented for the first time to shuffle full-length viral genomes. Together with our previously described system allowing the direct cloning of full-length infectious geminivirus genomes, it provided a unique opportunity to generate hundreds of "mosaic" virus genomes, directly testable for infectivity. A subset of 47 randomly chosen recombinants was sequenced, individually inoculated into tomato plants, and compared with the parental viruses. Surprisingly, our results showed that all recombinants were infectious and accumulated at levels comparable or intermediate to that of the parental clones. This indicates that, in our experimental system, despite the fact that the parental genomes differ by nearly 20%, lethal and/or large deleterious effects of recombination are very rare, in striking contrast to the common view that has emerged from previous studies published on other viruses.

A novel cloning strategy for isolating, genotyping and phenotyping genetic variants of geminiviruses.

BACKGROUND: Viruses of the genus Begomovirus (Geminiviridae) are emerging economically important plant viruses with a circular, single-stranded DNA genome. Previous studies have shown that geminiviruses and RNA viruses exhibit similar mutation frequencies, although geminiviruses are replicated by host DNA polymerases and RNA viruses by their own virus-encoded error-prone RNA-dependent RNA-polymerase. However, the phenotypic effects of naturally occurring mutations have never been extensively investigated in geminiviruses, particularly because, to be infectious, cloned viral genomes usually require sub-cloning as complete or partial tandem repeats into a binary vector from Agrobacterium tumefaciens. RESULTS: Using Tomato yellow leaf curl virus (TYLCV), we show here that infectivity can be obtained when only a 41-nucleotide region containing a highly conserved stem-loop is repeated. A binary vector containing this 41-nt region and a unique restriction site was created, allowing direct cloning of infectious monomeric viral genomes provided that they harbour the same restriction site at the corresponding nucleotide position. This experimental system, which can be transferable to other geminiviruses, was validated by analysis of the phenotypic effect of mutations appearing in TYLCV genomes in a single tomato host plant originally inoculated with a unique viral sequence. Fourteen full-length infectious genomes extracted from this plant were directly cloned and sequenced. The mutation frequency was 1.38 x 10-4 mutation per nucleotide sequenced, similar to that found previously for another begomovirus by sequencing PCR-amplified partial sequences. Interestingly, even in this minimal pool of analysed genomes, mutants with altered properties were readily identified, one of them being fitter and reducing plant biomass more drastically than the parental clone. CONCLUSION: The cloning strategy presented here is useful for any extensive phenotyping of geminivirus variants and particularly of artificially generated mutants or recombinants.