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

A new genus of alphasatellites associated with banana bunchy top virus in Southeast Asia

Autonomously replicating alphasatellites (family Alphasatellitidae) are frequently associated with plant single-stranded (ss)DNA viruses of the families Geminiviridae, Metaxyviridae, and Nanoviridae. Alphasatellites encode a single replication-initiator protein (Rep) similar to Rep proteins of helper viruses and depend on helper viruses for encapsidation, movement, and transmission. Costs versus benefits of alphasatellite-helper virus association are poorly understood. Our surveys in Southeast Asia (SEA) for wild and cultivated banana plants infected with banana bunchy top virus (BBTV, Nanoviridae) and Illumina sequencing reconstruction of their viromes revealed, in addition to a six-component BBTV genome, one to three distinct alphasatellites present in sixteen of twenty-four BBTV-infected plants. Comparative nucleotide and Rep protein sequence analyses classified these alphasatellites into four distinct species: two known species falling into the genus Muscarsatellite (subfamily Petromoalphasatellitinae) previously identified in SEA and two novel species falling into the tentative genus Banaphisatellite (subfamily Nanoalphasatellitinae) so far containing a single species recently identified in Africa. The banaphisatellites were found to be most related to members of the genus Fabenesatellite of subfamily Nanoalphasatellitinae and the genus Gosmusatellite of subfamily Geminialphasatellitinae, both infecting dicots. This suggests a dicot origin of banaphisatellites that got independently associated with distinct strains of monocot-infecting BBTV in Africa and SEA. Analysis of conserved sequence motifs in the common regions driving replication and gene expression of alphasatellites and BBTV strains revealed both differences and similarities, pointing at their ongoing co-evolution. An impact of alphasatellites on BBTV infection and evasion of RNA interference-based antiviral defences was evaluated by measuring relative abundance of BBTV genome components and alphasatellites and by profiling BBTV- and alphasatellite-derived small interfering RNAs. Taken together, our findings shed new light on the provenance of alphasatellites, their co-evolution with helper viruses, and potential mutual benefits of their association.

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.

Self-replicating alphasatellites are frequently associated with plant ssDNA viruses. Their origin and costs versus benefits for helper virus replication, antiviral defense evasion and transmission by insect vectors are poorly understood. Here we describe identification in Africa and in depth molecular and biological characterization of a newly emerging alphasatellite of BBTV, a multicomponent ssDNA babuvirus causing one of the most economically-important diseases of monocotyledonous bananas and plantains. Phylogenetically, this alphasatellite represents a novel genus and is more related to alphasatellites of nanoviruses infecting dicot hosts than to other BBTV alphasatellites previously identified only in Asia. Consistent with its hypothetical dicot origin, cloned alphasatellite and BBTV can establish systemic infection in a model dicot plant, followed by recovery and symptomless alphasatellite persistence. In banana plants, alphasatellite competes for the host replication and transcription machinery and accumulates at high levels, thereby reducing loads of the helper virus, modifying relative abundance of its components and interfering with its acquisition and transmission by aphids. On the other hand, plant antiviral defenses silence alphasatellite gene expression at both transcriptional and posttranscriptional levels, generating highly-abundant 21, 22 and 24 nucleotide small interfering RNAs, suggesting that alphasatellite may serve as a decoy protecting its helper virus from gene silencing.

Exploring diverse roles of micro RNAs in banana: Current status and future prospective

Micro RNAs (miRNAs) are 20-24 nucleotides long non-coding RNA sequences identified and characterized in multiple plant and animal systems. miRNAs play multifarious roles ranging from plant development to stress tolerance by synchronizing physiological processes at the level of transcription and translation. Banana is a major horticultural crop with colossal production worldwide. Despite the recent encouraging developments, the information on functions of miRNAs in banana plants is still in its infancy. The available literature pertaining to miRNAs in banana plants hints towards their contribution as master regulators in crucial physiological processes for instance abiotic stress responses, pathogenic defence response, fruit ripening and so on. This review is focused on biogenesis of miRNAs, their identification and deciphering their respective roles in banana plants with special emphasis on abiotic stress responses, plant immune responses, fruit ripening and storage. Based on the prior reports, we identified a few miRNAs with prospective roles in stress tolerance and illustrated the potential applications of miRNAs in banana crop improvement utilizing recent biotechnological tools such as CRISPR cas9, RNAi and the nano particle based foliar spray of miRNAs. The review briefly explained the future directions in banana research with a special emphasis on miRNA regulatory networks and agronomic traits improvement. Finally, future domains in miRNA research in plants and their possible applications towards crop improvement in agriculture are described briefly.

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.

Infectivity of an Infectious Clone of Banana Streak CA Virus in A-Genome Bananas (Musa acuminata ssp.)

We have characterized the complete genome sequence of an Australian isolate of banana streak CA virus (BSCAV). A greater-than-full-length, cloned copy of the virus genome was assembled and agroinoculated into five tissue-cultured plants of nine different Musa acuminata banana accessions. BSCAV was highly infectious in all nine accessions. All five inoculated plants from eight accessions developed symptoms by 28 weeks post-inoculation, while all five plants of M. acuminata AA subsp. zebrina remained symptomless. Symptoms were mild in six accessions but were severe in Khae Phrae (M. acuminata subsp. siamea) and the East African Highland banana accession Igisahira Gisanzwe. This is the first full-length BSCAV genome sequence reported from Australia and the first report of the infectivity of an infectious clone of banana streak virus.

Virus Detection by High-Throughput Sequencing of Small RNAs: Large-Scale Performance Testing of Sequence Analysis Strategies

Recent developments in high-throughput sequencing (HTS), also called next-generation sequencing (NGS), technologies and bioinformatics have drastically changed research on viral pathogens and spurred growing interest in the field of virus diagnostics. However, the reliability of HTS-based virus detection protocols must be evaluated before adopting them for diagnostics. Many different bioinformatics algorithms aimed at detecting viruses in HTS data have been reported but little attention has been paid thus far to their sensitivity and reliability for diagnostic purposes. Therefore, we compared the ability of 21 plant virology laboratories, each employing a different bioinformatics pipeline, to detect 12 plant viruses through a double-blind large-scale performance test using 10 datasets of 21- to 24-nucleotide small RNA (sRNA) sequences from three different infected plants. The sensitivity of virus detection ranged between 35 and 100% among participants, with a marked negative effect when sequence depth decreased. The false-positive detection rate was very low and mainly related to the identification of host genome-integrated viral sequences or misinterpretation of the results. Reproducibility was high (91.6%). This work revealed the key influence of bioinformatics strategies for the sensitive detection of viruses in HTS sRNA datasets and, more specifically (i) the difficulty in detecting viral agents when they are novel or their sRNA abundance is low, (ii) the influence of key parameters at both assembly and annotation steps, (iii) the importance of completeness of reference sequence databases, and (iv) the significant level of scientific expertise needed when interpreting pipeline results. Overall, this work underlines key parameters and proposes recommendations for reliable sRNA-based detection of known and unknown viruses.

Small RNA-Omics for Plant Virus Identification, Virome Reconstruction, and Antiviral Defense Characterization

RNA interference (RNAi)-based antiviral defense generates small interfering RNAs that represent the entire genome sequences of both RNA and DNA viruses as well as viroids and viral satellites. Therefore, deep sequencing and bioinformatics analysis of small RNA population (small RNA-ome) allows not only for universal virus detection and genome reconstruction but also for complete virome reconstruction in mixed infections. Viral infections (like other stress factors) can also perturb the RNAi and gene silencing pathways regulating endogenous gene expression and repressing transposons and host genome-integrated endogenous viral elements which can potentially be released from the genome and contribute to disease. This review describes the application of small RNA-omics for virus detection, virome reconstruction and antiviral defense characterization in cultivated and non-cultivated plants. Reviewing available evidence from a large and ever growing number of studies of naturally or experimentally infected hosts revealed that all families of land plant viruses, their satellites and viroids spawn characteristic small RNAs which can be assembled into contigs of sufficient length for virus, satellite or viroid identification and for exhaustive reconstruction of complex viromes. Moreover, the small RNA size, polarity and hotspot profiles reflect virome interactions with the plant RNAi machinery and allow to distinguish between silent endogenous viral elements and their replicating episomal counterparts. Models for the biogenesis and functions of small interfering RNAs derived from all types of RNA and DNA viruses, satellites and viroids as well as endogenous viral elements are presented and discussed.

Small RNA-based interactions between rice and the viruses which cause the tungro disease

Rice tungro disease is caused by a complex of two viruses, Rice tungro bacilliform virus (RTBV) and Rice tungro spherical virus (RTSV). To examine the RNAi-based defence response in rice during tungro disease, we characterized the virus-derived small RNAs and miRNAs by Deep Sequencing. We found that, while 21 nt/22 nt (nucleotide) siRNAs are predominantly produced in a continuous, overlapping and asymmetrical manner from RTBV, siRNA accumulation from RTSV were negligible. Additionally, 54 previously known miRNAs from rice, predicted to be regulating genes involved in plant defence, hormone signaling and developmental pathways were differentially expressed in the infected samples, compared to the healthy ones. This is the first study of sRNA profile of tungro virus complex from infected rice plants. The biased response of the host antiviral machinery against the two viruses and the differentially-expressed miRNAs are novel observations, which entail further studies.

Small RNA-Omics for Virome Reconstruction and Antiviral Defense Characterization in Mixed Infections of Cultivated Solanum Plants

In plants, RNA silencing-based antiviral defense generates viral small RNAs (sRNAs) faithfully representing the viral genomes. We employed sRNA sequencing and bioinformatics (sRNA-omics) to characterize antiviral defense and to reconstruct the full genomic sequences and their variants in the evolving viral quasispecies in cultivated solanaceous plants carrying mixed infections. In naturally infected Solanum tuberosum (potato), one case study revealed a virome comprising Potato virus Y (genus Potyvirus) and Potato virus X (genus Potexvirus), which was reconstructed by de novo-assembling separate genome-size sRNA contigs. Another case study revealed a virome comprising NTN and O strains of Potato virus Y, whose sRNAs assembled in chimeric contigs, which could be disentangled on the basis of reference genome sequences. Both viromes were stable in vegetative potato progeny. In a cross-protection trial of Solanum lycopersicum (tomato), the supposedly protective mild strain CH2 of Pepino mosaic virus (genus Potexvirus) was tested for protection against strain LP of the same virus. Reciprocal mechanical inoculations eventually resulted in co-infection of all individual plants with CH2 and LP strains, reconstructed as separate sRNA contigs. LP invasions into CH2-preinfected plants and vice versa were accompanied by alterations of consensus genome sequences in viral quasispecies, indicating a potential risk of cross-protection measures. Additionally, the study also revealed, by reconstruction from sRNAs, the presence of the mechanically nontransmissible Southern tomato virus (genus Amalgavirus) in some plants. Our in-depth analysis of sRNA sizes, 5'-nucleotide frequencies and hotspot maps revealed similarities in sRNA-generating mechanisms in potato and tomato, differential silencing responses to virome components and potential for sRNA-directed cross-targeting between viral strains which could not, however, prevent the formation of stable viromes.

Ribosome Shunting, Polycistronic Translation, and Evasion of Antiviral Defenses in Plant Pararetroviruses and Beyond

Viruses have compact genomes and usually translate more than one protein from polycistronic RNAs using leaky scanning, frameshifting, stop codon suppression or reinitiation mechanisms. Viral (pre-)genomic RNAs often contain long 5′-leader sequences with short upstream open reading frames (uORFs) and secondary structure elements, which control both translation initiation and replication. In plants, viral RNA and DNA are targeted by RNA interference (RNAi) generating small RNAs that silence viral gene expression, while viral proteins are recognized by innate immunity and autophagy that restrict viral infection. In this review we focus on plant pararetroviruses of the family Caulimoviridae and describe the mechanisms of uORF- and secondary structure-driven ribosome shunting, leaky scanning and reinitiation after translation of short and long uORFs. We discuss conservation of these mechanisms in different genera of Caulimoviridae, including host genome-integrated endogenous viral elements, as well as in other viral families, and highlight a multipurpose use of the highly-structured leader sequence of plant pararetroviruses in regulation of translation, splicing, packaging, and reverse transcription of pregenomic RNA (pgRNA), and in evasion of RNAi. Furthermore, we illustrate how targeting of several host factors by a pararetroviral effector protein can lead to transactivation of viral polycistronic translation and concomitant suppression of antiviral defenses. Thus, activation of the plant protein kinase target of rapamycin (TOR) by the Cauliflower mosaic virus transactivator/viroplasmin (TAV) promotes reinitiation of translation after long ORFs on viral pgRNA and blocks antiviral autophagy and innate immunity responses, while interaction of TAV with the plant RNAi machinery interferes with antiviral silencing.

Viral small RNAs reveal the genomic variations of three grapevine vein clearing virus quasispecies populations

Viral small RNAs (vsRNAs) include viral small interfering RNAs (vsiRNAs) that are initiators and products of RNA silencing, and small RNAs that are derived from viral RNAs with function still unknown. Sequencing of vsRNAs allows assembling of viral genomes and revelation of viral population variations at genomic levels. Grapevine vein clearing virus (GVCV) is a new member of the family Caulimoviridae whose DNA genome is replicated by reverse transcription of pre-genomic RNA molecules. In this short report, three genomic sequences of GVCV were assembled from vsRNAs that were isolated and sequenced from three individual grapevines in commercial vineyards and compared to the GVCV-CHA reference genome. Profiles of single nucleotide polymorphism among three viral populations indicated a closer relatedness between two populations in different grape cultivars at the same location than those in the same grape cultivar at different locations, suggesting the spread of GVCV populations among vineyards of close proximity. Classic types of vsiRNAs (21-nt, 22-nt, and 24-nt) were found in the three GVCV vsiRNA populations, but these did not produce alignment hotspots on the GVCV-CHA reference genome. The number of 36-nt reads is the highest among vsRNAs, the role of these vsRNAs remains unclear. The analysis of vsRNAs provides a first holistic picture of genomic variations among GVCV viral quasispecies populations that help monitor epidemics and evolution of GVCV populations, an emerging virus that is becoming a threat to grape production in the Midwest region of the USA.

The Evolution of Antiviral Defense Systems

Self-replicating genetic material presumably provided the architecture necessary for generating the last universal ancestor of all nucleic-acid-based life. As biological complexity increased in the billions of years that followed, the same genetic material also morphed into a wide spectrum of viruses and other parasitic genetic elements. The resulting struggle for existence drove the evolution of host defenses, giving rise to a perpetual arms race. This Perspective summarizes the antiviral mechanisms evident across the tree of life, discussing each in their evolutionary context to postulate how the coevolution of host and pathogen shaped the cellular antiviral defenses we know today.

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.

MISIS-2: A bioinformatics tool for in-depth analysis of small RNAs and representation of consensus master genome in viral quasispecies

In most eukaryotes, small RNA (sRNA) molecules such as miRNAs, siRNAs and piRNAs regulate gene expression and repress transposons and viruses. AGO/PIWI family proteins sort functional sRNAs based on size, 5'-nucleotide and other sequence features. In plants and some animals, viral sRNAs are extremely diverse and cover the entire viral genome sequences, which allows for de novo reconstruction of a complete viral genome by deep sequencing and bioinformatics analysis of viral sRNAs. Previously, we have developed a tool MISIS to view and analyze sRNA maps of viruses and cellular genome regions which spawn multiple sRNAs. Here we describe a new release of MISIS, MISIS-2, which enables to determine and visualize a consensus sequence and count sRNAs of any chosen sizes and 5'-terminal nucleotide identities. Furthermore we demonstrate the utility of MISIS-2 for identification of single nucleotide polymorphisms (SNPs) at each position of a reference sequence and reconstruction of a consensus master genome in evolving viral quasispecies. MISIS-2 is a Java standalone program. It is freely available along with the source code at the website http://www.fasteris.com/apps.

Characterization by Small RNA Sequencing of Taro Bacilliform CH Virus (TaBCHV), a Novel Badnavirus

RNA silencing is an antiviral immunity that regulates gene expression through the production of small RNAs (sRNAs). In this study, deep sequencing of small RNAs was used to identify viruses infecting two taro plants. Blast searching identified five and nine contigs assembled from small RNAs of samples T1 and T2 matched onto the genome sequences of badnaviruses in the family Caulimoviridae. Complete genome sequences of two isolates of the badnavirus determined by sequence specific amplification comprised of 7,641 nucleotides and shared overall nucleotide similarities of 44.1%‒55.8% with other badnaviruses. Six open reading frames (ORFs) were identified on the plus strand, showed amino acid similarities ranging from 59.8% (ORF3) to 10.2% (ORF6) to the corresponding proteins encoded by other badnaviruses. Phylogenetic analysis also supports that the virus is a new member in the genus Badnavirus. The virus is tentatively named as Taro bacilliform CH virus (TaBCHV), and it is the second badnavirus infecting taro plants, following Taro bacilliform virus (TaBV). In addition, analyzes of viral derived small RNAs (vsRNAs) from TaBCHV showed that almost equivalent number of vsRNAs were generated from both strands and the most abundant vsRNAs were 21 nt, with uracil bias at 5' terminal. Furthermore, TaBCHV vsRNAs were asymmetrically distributed on its entire circular genome at both orientations with the hotspots mainly generated in the ORF5 region.

RNA based viral silencing suppression in plant pararetroviruses

The 35S promoter of cauliflower mosaic virus and that of other plant pararetroviruses gives rise to an RNA, which is both a pre-genome and a polycistronic mRNA. The 600 nucleotide long very structured leader of this RNA is also transcribed separately. The resulting 8S RNA is then converted to a double strand giving rise to a huge set of siRNAs, which suppress silencing. In this Mini-Review I discuss how this versatile stretch of 600 nts constitutes a masterpiece of evolution.

Genome Sequence of Banana Streak MY Virus from the Pacific Ocean Island of Tonga

Banana streak disease is caused by a variety of banana-infecting badnaviruses. A genome of the episomal form of a banana streak MY virus was recovered from an infected banana plant sampled on Vava'u Island, Tonga, and shares >98% pairwise identity with the six other genomes available in public databases.

Interactions of Rice tungro bacilliform pararetrovirus and its protein P4 with plant RNA-silencing machinery

Small interfering RNA (siRNA)-directed gene silencing plays a major role in antiviral defense. Virus-derived siRNAs inhibit viral replication in infected cells and potentially move to neighboring cells, immunizing them from incoming virus. Viruses have evolved various ways to evade and suppress siRNA production or action. Here, we show that 21-, 22-, and 24-nucleotide (nt) viral siRNAs together constitute up to 19% of total small RNA population of Oryza sativa plants infected with Rice tungro bacilliform virus (RTBV) and cover both strands of the RTBV DNA genome. However, viral siRNA hotspots are restricted to a short noncoding region between transcription and reverse-transcription start sites. This region generates double-stranded RNA (dsRNA) precursors of siRNAs and, in pregenomic RNA, forms a stable secondary structure likely inaccessible to siRNA-directed cleavage. In transient assays, RTBV protein P4 suppressed cell-to-cell spread of silencing but enhanced cell-autonomous silencing, which correlated with reduced 21-nt siRNA levels and increased 22-nt siRNA levels. Our findings imply that RTBV generates decoy dsRNA that restricts siRNA production to the structured noncoding region and thereby protects other regions of the viral genome from repressive action of siRNAs, while the viral protein P4 interferes with cell-to-cell spread of antiviral silencing.