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.

A newly emerging alphasatellite affects banana bunchy top virus replication, transcription, siRNA production and transmission by aphids.

Banana bunchy top virus (BBTV) is a six-component ssDNA virus (genus Babuvirus, family Nanoviridae) transmitted by aphids, infecting monocots (mainly species in the family Musaceae) and likely originating from South-East Asia where it is frequently associated with self-replicating alphasatellites. Illumina sequencing analysis of banana aphids and leaf samples from Africa revealed an alphasatellite that should be classified in a new genus, phylogenetically related to alphasatellites of nanoviruses infecting dicots. Alphasatellite DNA was encapsidated by BBTV coat protein and accumulated at high levels in plants and aphids, thereby reducing helper virus loads, altering relative abundance (formula) of viral genome components and interfering with virus transmission by aphids. BBTV and alphasatellite clones infected dicot Nicotiana benthamiana, followed by recovery and symptomless persistence of alphasatellite, and BBTV replication protein (Rep), but not alphasatellite Rep, induced leaf chlorosis. Transcriptome sequencing revealed 21, 22 and 24 nucleotide small interfering (si)RNAs covering both strands of the entire viral genome, monodirectional Pol II transcription units of viral mRNAs and pervasive transcription of each component and alphasatellite in both directions, likely generating double-stranded precursors of viral siRNAs. Consistent with the latter hypothesis, viral DNA formulas with and without alphasatellite resembled viral siRNA formulas but not mRNA formulas. Alphasatellite decreased transcription efficiency of DNA-N encoding a putative aphid transmission factor and increased relative siRNA production rates from Rep- and movement protein-encoding components. Alphasatellite itself spawned the most abundant siRNAs and had the lowest mRNA transcription rate. Collectively, following African invasion, BBTV got associated with an alphasatellite likely originating from a dicot plant and interfering with BBTV replication and transmission. Molecular analysis of virus-infected banana plants revealed new features of viral DNA transcription and siRNA biogenesis, both affected by alphasatellite. Costs and benefits of alphasatellite association with helper viruses are discussed.

Endogenous viral elements are targeted by RNA silencing pathways in banana.

Endogenous banana streak virus (eBSV) integrants derived from three distinct species, present in Musa balbisiana (B) but not Musa acuminata (A) banana genomes are able to reconstitute functional episomal viruses causing banana streak disease in interspecific triploid AAB banana hybrids but not in the diploid (BB) parent line, which harbours identical eBSV loci. Here, we investigated the regulation of these eBSV. In-depth characterization of siRNAs, transcripts and methylation derived from eBSV using Illumina and bisulfite sequencing were carried out on eBSV-free Musa acuminata AAA plants and BB or AAB banana plants with eBSV. eBSV loci produce low-abundance transcripts covering most of the viral sequence and generate predominantly 24-nt siRNAs. siRNA accumulation is restricted to duplicated and inverted viral sequences present in eBSV. Both siRNA-accumulating and nonaccumulating sequences of eBSV in BB plants are heavily methylated in all three CG, CHG and CHH contexts. Our data suggest that eBSVs are controlled at the epigenetic level in BB diploids. This regulation not only prevents their awakening and systemic infection of the plant but is also probably involved in the inherent resistance of the BB plants to mealybug-transmitted viral infection. These findings are thus of relevance to other plant resources hosting integrated viruses.

Extrachromosomal viral DNA produced by transcriptionally active endogenous viral elements in non-infected banana hybrids impedes quantitative PCR diagnostics of banana streak virus infections in banana hybrids.

The main edible and cultivated banana varieties are intra- and interspecific hybrids of the two main Musa species, Musa acuminata and Musa balbisiana, having diploid genomes denoted A and B, respectively. The B genome naturally hosts sequences of banana streak virus (BSV) named endogenous BSV (eBSV). Upon stress, eBSVs are identified as the origin of BSV infection for at least three BSV species, causing banana streak disease. For each of the three species, BSV and eBSV share >99.9 % sequence identity, complicating PCR-based diagnosis of viral infection in the B genome-containing bananas. Here, we designed a quantitative PCR-based method to only quantify episomal BSV particles produced, overcoming the limitation of eBSV also being detected by qPCR by using it as a 'calibrator'. However, our results revealed unexpected variation of eBSV amplification in calibrator plants composed of a clonal population of 53 replicating virus-free banana hybrids with the same AAB genotype. Our in-depth molecular analyses suggest that this calibrator variation is due to the variable abundance of non-encapsidated extrachromosomal viral DNA, likely produced via the transcription of eBSVs, followed by occasional reverse transcription. We also present evidence that accumulation of viral transcripts in AAB plants is downregulated both at post-transcriptional and transcriptional levels by an RNA interference mechanism that keeps the plants free of virus infection. Finally, we recommend that such eBSV amplification variation be taken into account to establish a quantitative viral diagnostic for banana plants with the B genome.

How endogenous plant pararetroviruses shed light on Musa evolution.

BACKGROUND AND AIMS: Banana genomes harbour numerous copies of viral sequences derived from banana streak viruses (BSVs) - dsDNA viruses belonging to the family Caulimoviridae.These viral integrants (eBSVs) are mostly defective, probably as a result of 'pseudogenization' driven by host genome evolution. However, some can give rise to infection by releasing a functional viral genome following abiotic stresses. These distinct infective eBSVs correspond to the three main widespread BSV species (BSOLV, BSGFV and BSIMV), fully described within the Musa balbisiana B genomes of the seedy diploid 'Pisang Klutuk Wulung' (PKW). METHODS: We characterize eBSV distribution among a Musa sampling including seedy BB diploids and interspecific hybrids with Musa acuminate exhibiting different levels of ploidy for the B genome (ABB, AAB, AB). We used representative samples of the two areas of sympatry between M. acuminate and M. balbisiana species representing the native area of the most widely cultivated AAB cultivars (in India and in East Asia, ranging from the Philippines to New Guinea). Seventy-seven accessions were characterized using eBSV-related PCR markers and Southern hybridization approaches. We coded both sets of results to create a common dissimilarity matrix with which to interpret eBSV distribution. KEY RESULTS: We propose a Musa phylogeny driven by the M. balbisiana genome based on a dendrogram resulting from a joint neighbour-joining analysis of the three BSV species, showing for the first time lineages between BB and ABB/AAB hybrids. eBSVs appear to be relevant phylogenetic markers that can illustrate theM. balbisiana phylogeography story. CONCLUSION: The theoretical implications of this study for further elucidation of the historical and geographical process of Musa domestication are numerous. Discovery of banana plants with B genome non-infective for eBSV opens the way to the introduction of new genitors in programmes of genetic banana improvement.

Evasion of short interfering RNA-directed antiviral silencing in Musa acuminata persistently infected with six distinct banana streak pararetroviruses.

UNLABELLED: Vegetatively propagated crop plants often suffer from infections with persistent RNA and DNA viruses. Such viruses appear to evade the plant defenses that normally restrict viral replication and spread. The major antiviral defense mechanism is based on RNA silencing generating viral short interfering RNAs (siRNAs) that can potentially repress viral genes posttranscriptionally through RNA cleavage and transcriptionally through DNA cytosine methylation. Here we examined the RNA silencing machinery of banana plants persistently infected with six pararetroviruses after many years of vegetative propagation. Using deep sequencing, we reconstructed consensus master genomes of the viruses and characterized virus-derived and endogenous small RNAs. Consistent with the presence of endogenous siRNAs that can potentially establish and maintain DNA methylation, the banana genomic DNA was extensively methylated in both healthy and virus-infected plants. A novel class of abundant 20-nucleotide (nt) endogenous small RNAs with 5'-terminal guanosine was identified. In all virus-infected plants, 21- to 24-nt viral siRNAs accumulated at relatively high levels (up to 22% of the total small RNA population) and covered the entire circular viral DNA genomes in both orientations. The hotspots of 21-nt and 22-nt siRNAs occurred within open reading frame (ORF) I and II and the 5' portion of ORF III, while 24-nt siRNAs were more evenly distributed along the viral genome. Despite the presence of abundant viral siRNAs of different size classes, the viral DNA was largely free of cytosine methylation. Thus, the virus is able to evade siRNA-directed DNA methylation and thereby avoid transcriptional silencing. This evasion of silencing likely contributes to the persistence of pararetroviruses in banana plants. IMPORTANCE: We report that DNA pararetroviruses in Musa acuminata banana plants are able to evade DNA cytosine methylation and transcriptional gene silencing, despite being targeted by the host silencing machinery generating abundant 21- to 24-nucleotide short interfering RNAs. At the same time, the banana genomic DNA is extensively methylated in both healthy and virus-infected plants. Our findings shed light on the siRNA-generating gene silencing machinery of banana and provide a possible explanation why episomal pararetroviruses can persist in plants whereas true retroviruses with an obligatory genome-integration step in their replication cycle do not exist in plants.

Cannabis Virome Reconstruction and Antiviral RNAi Characterization through Small RNA Sequencing.

Viral infections pose an emerging threat to hemp (Cannabis sativa) cultivation. We used Illumina small (s)RNA sequencing for virome reconstruction and characterization of antiviral RNA interference (RNAi) in monoecious and dioecious hemp varieties, which exhibited different virus-like symptoms. Through de novo and reference-based sRNA assembly, we identified and reconstructed Cannabis cryptic virus (family Partitiviridae), Cannabis sativa mitovirus 1 (Mitoviridae) and Grapevine line pattern virus (Bromoviridae) as well as a novel virus tentatively classified into Partitiviridae. Members of both Partitiviridae and Bromoviridae were targeted by antiviral RNAi, generating 21 nt and, less abundant, 22 nt sRNAs from both strands of the entire virus genome, suggesting the involvement of Dicer-like (DCL) 4 and DCL2 in viral sRNA biogenesis, respectively. Mitovirus sRNAs represented predominantly the positive-sense strand and had a wider size range, with the 21 nt class being most abundant on both strands. For all viruses, 21 and 22 nt sRNAs had predominantly 5'-terminal uridine or cytosine, suggesting their binding to antiviral Argonaute (AGO) 1 and AGO5, respectively. As no clear association of any virus with symptoms was observed, further studies should clarify if these viruses individually or in combination can cause hemp diseases.

Musa balbisiana genome reveals subgenome evolution and functional divergence.

Banana cultivars (Musa ssp.) are diploid, triploid and tetraploid hybrids derived from Musa acuminata and Musa balbisiana. We presented a high-quality draft genome assembly of M. balbisiana with 430 Mb (87%) assembled into 11 chromosomes. We identified that the recent divergence of M. acuminata (A-genome) and M. balbisiana (B-genome) occurred after lineage-specific whole-genome duplication, and that the B-genome may be more sensitive to the fractionation process compared to the A-genome. Homoeologous exchanges occurred frequently between A- and B-subgenomes in allopolyploids. Genomic variation within progenitors resulted in functional divergence of subgenomes. Global homoeologue expression dominance occurred between subgenomes of the allotriploid. Gene families related to ethylene biosynthesis and starch metabolism exhibited significant expansion at the pathway level and wide homoeologue expression dominance in the B-subgenome of the allotriploid. The independent origin of 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) homoeologue gene pairs and tandem duplication-driven expansion of ACO genes in the B-subgenome contributed to rapid and major ethylene production post-harvest in allotriploid banana fruits. The findings of this study provide greater context for understanding fruit biology, and aid the development of tools for breeding optimal banana cultivars.

A reliable IC One-step RT-PCR method for the detection of BBrMV to ensure safe exchange of Musa germplasm.

An immunocapture (IC) One-step RT-PCR assay was developed to improve the detection of Banana bract mosaic virus (BBrMV) in single and bulked samples of banana plants. In this paper, an atypical strain of BBrMV was described, the BBrMV "Ref" strain, and we showed that detection with available BBrMV tools using ELISA and RT-PCR approaches was not reliable. Primer sets Bract N1/NR and N2/NR specific to BBrMV were designed and used in RT-PCR and IC-RT-PCR assays with two commercial kits that allow the RT and the PCR reactions to take place simultaneously in the same tube. The new assay enabled detection of BBrMV in leaf extract diluted up to 1 x 10(-10) and in bulked samples of 10 plants, and was proposed as a new international standard to index BBrMV.

The genome of African yam (Dioscorea cayenensis-rotundata complex) hosts endogenous sequences from four distinct Badnavirus species.

Several endogenous viral elements (EVEs) have been identified in plant genomes, including endogenous pararetroviruses (EPRVs). Here, we report the first characterization of EPRV sequences in the genome of African yam of the Dioscorea cayenensis-rotundata complex. We propose that these sequences should be termed 'endogenous Dioscorea bacilliform viruses' (eDBVs). Molecular characterization of eDBVs shows that they constitute sequences originating from various parts of badnavirus genomes, resulting in a mosaic structure that is typical of most EPRVs characterized to date. Using complementary molecular approaches, we show that eDBVs belong to at least four distinct Badnavirus species, indicating multiple, independent, endogenization events. Phylogenetic analyses of eDBVs support and enrich the current taxonomy of yam badnaviruses and lead to the characterization of a new Badnavirus species in yam. The impact of eDBVs on diagnosis, yam germplasm conservation and movement, and breeding is discussed.

High genetic variability and evidence for plant-to-plant transfer of Banana mild mosaic virus.

A total of 154 partial nucleotide sequences within the Banana mild mosaic virus (BanMMV) ORF1, which encodes the viral RNA-dependent RNA polymerase (RdRp), was obtained from 68 distinct infected banana accessions originating from various locations worldwide. The 310 nt sequences displayed a high level of variability with a mean pairwise nucleotide sequence divergence level of 20.4 %. This situation resulted essentially from a high rate of synonymous mutations. A similar analysis was performed for a limited selection of 10 banana accessions (30 sequences) on the region comprising approximately the last 310 nt of the BanMMV genome. This region corresponds to the 3' end of ORF5, which encodes the coat protein (234 nt), and to the 3' non-coding region. This analysis confirmed the high level of diversity observed in the RdRp dataset, characterized by a high level of synonymous mutations. Analysis of intra-host diversity indicated the existence of two distinct situations, with some plants containing only closely related sequence variants, whereas others contained widely divergent isolates. Analyses indicated that BanMMV genetic diversity is not structured by the geographical origin of the infected Musa accessions or by their genotype. This situation may be, in part, explained by the exchange of banana germplasm between different parts of the world and also by plant-to-plant transfer of virus isolates, the evidence for which is, for the first time, provided by this study.