Small RNAs Participate in Plant-Virus Interaction and Their Application in Plant Viral Defense

The Biosynthesis Process of Small RNA and Its Pivotal Roles in Plant Development

In the realm of plant biology, small RNAs (sRNAs) are imperative in the orchestration of gene expression, playing pivotal roles across a spectrum of developmental sequences and responses to environmental stressors. The biosynthetic cascade of sRNAs is characterized by an elaborate network of enzymatic pathways that meticulously process double-stranded RNA (dsRNA) precursors into sRNA molecules, typically 20 to 30 nucleotides in length. These sRNAs, chiefly microRNAs (miRNAs) and small interfering RNAs (siRNAs), are integral in guiding the RNA-induced silencing complex (RISC) to selectively target messenger RNAs (mRNAs) for post-transcriptional modulation. This regulation is achieved either through the targeted cleavage or the suppression of translational efficiency of the mRNAs. In plant development, sRNAs are integral to the modulation of key pathways that govern growth patterns, organ differentiation, and developmental timing. The biogenesis of sRNA itself is a fine-tuned process, beginning with transcription and proceeding through a series of processing steps involving Dicer-like enzymes and RNA-binding proteins. Recent advances in the field have illuminated the complex processes underlying the generation and function of small RNAs (sRNAs), including the identification of new sRNA categories and the clarification of their involvement in the intercommunication among diverse regulatory pathways. This review endeavors to evaluate the contemporary comprehension of sRNA biosynthesis and to underscore the pivotal role these molecules play in directing the intricate performance of plant developmental processes.

Asymmetric bulges within hairpin RNA transgenes influence small RNA size, secondary siRNA production and viral defence

Small RNAs (sRNAs) are essential for normal plant development and range in size classes of 21-24 nucleotides. The 22nt small interfering RNAs (siRNAs) and miRNAs are processed by Dicer-like 2 (DCL2) and DCL1 respectively and can initiate secondary siRNA production from the target transcript. 22nt siRNAs are under-represented due to competition between DCL2 and DCL4, while only a small number of 22nt miRNAs exist. Here we produce abundant 22nt siRNAs and other siRNA size classes using long hairpin RNA (hpRNA) transgenes. By introducing asymmetric bulges into the antisense strand of hpRNA, we shifted the dominant siRNA size class from 21nt of the traditional hpRNA to 22, 23 and 24nt of the asymmetric hpRNAs. The asymmetric hpRNAs effectively silenced a β-glucuronidase (GUS) reporter transgene and the endogenous ethylene insensitive-2 (EIN2) and chalcone synthase (CHS) genes. Furthermore, plants containing the asymmetric hpRNA transgenes showed increased amounts of 21nt siRNAs downstream of the hpRNA target site compared to plants with the traditional hpRNA transgenes. This indicates that these asymmetric hpRNAs are more effective at inducing secondary siRNA production to amplify silencing signals. The 22nt asymmetric hpRNA constructs enhanced virus resistance in plants compared to the traditional hpRNA constructs.

A novel ilarvirus protein CP-RT is expressed via stop codon readthrough and suppresses RDR6-dependent RNA silencing

Ilarviruses are a relatively understudied but important group of plant RNA viruses that includes a number of crop pathogens. Their genomes comprise three RNA segments encoding two replicase subunits, movement protein, coat protein (CP), and (in some ilarvirus subgroups) a protein that suppresses RNA silencing. Here we report that, in many ilarviruses, RNA3 encodes an additional protein (termed CP-RT) as a result of ribosomal readthrough of the CP stop codon into a short downstream readthrough (RT) ORF. Using asparagus virus 2 as a model, we find that CP-RT is expressed in planta where it functions as a weak suppressor of RNA silencing. CP-RT expression is essential for persistent systemic infection in leaves and shoot apical meristem. CP-RT function is dependent on a putative zinc-finger motif within RT. Replacing the asparagus virus 2 RT with the RT of an ilarvirus from a different subgroup restored the ability to establish persistent infection. These findings open up a new avenue for research on ilarvirus silencing suppression, persistent meristem invasion and vertical transmission.

In silico identification of papaya genome-encoded microRNAs to target begomovirus genes in papaya leaf curl disease

Papaya leaf curl disease (PaLCuD) is widespread and classified in the genus begomovirus (Geminiviridae), disseminated by the vector whitefly Bemisia tabaci. RNA interference (RNAi)-based antiviral innate immunity stands as a pivotal defense mechanism and biological process in limiting viral genomes to manage plant diseases. The current study aims to identify and analyze Carica Papaya locus-derived capa-microRNAs with predicted potential for targeting divergent begomovirus species-encoded mRNAs using a 'four integrative in silico algorithms' approach. This research aims to experimentally activate the RNAi catalytic pathway using in silico-predicted endogenous capa-miRNAs and create papaya varieties capable of assessing potential resistance against begomovirus species and monitoring antiviral capabilities. This study identified 48 predicted papaya locus-derived candidates from 23 miRNA families, which were further investigated for targeting begomovirus genes. Premised all the four algorithms combined, capa-miR5021 was the most anticipated miRNA followed by capa-miR482, capa-miR5658, capa-miR530b, capa-miR3441.2, and capa-miR414 'effective' papaya locus-derived candidate capa-miRNA and respected putative binding sites for targets at the consensus nucleotide position. It was predicted to bind and target mostly to AC1 gene of the complementary strand and the AV1 gene of the virion strand of different begomovirus isolates, which were associated with replication-associated protein and encapsidation, respectively, during PaLCuD. These miRNAs were also found targeting betaC1 gene of betasatellite which were associated with retardation in leaf growth and developmental abnormalities with severe symptoms during begomovirus infection. To validate target prediction accuracy, we created an integrated Circos plot for comprehensive visualization of host-virus interaction. In silico-predicted papaya genome-wide miRNA-mediated begomovirus target gene regulatory network corroborated interactions that permit in vivo analysis, which could provide biological material and valuable evidence, leading to the development of begomovirus-resistant papaya plants. The integrative nature of our research positions it at the forefront of efforts to ensure the sustainable cultivation of papaya, particularly in the face of evolving pathogenic threats. As we move forward, the knowledge gained from this study provides a solid foundation for continued exploration and innovation in the field of papaya virology, and to the best of our knowledge, this study represents a groundbreaking endeavor, undertaken for the first time in the context of PaLCuD research.

Differential regulation of miRNAs involved in the susceptible and resistance responses of wheat cultivars to wheat streak mosaic virus and Triticum mosaic virus

BACKGROUND

Wheat streak mosaic virus (WSMV) and Triticum mosaic virus (TriMV) are components of the wheat streak mosaic virus disease complex in the Great Plains region of the U.S.A. and elsewhere. Co-infection of wheat with WSMV and TriMV causes synergistic interaction with more severe disease symptoms compared to single infections. Plants are equipped with multiple antiviral mechanisms, of which regulation of microRNAs (miRNAs) is a potentially effective constituent. In this investigation, we have analyzed the total and relative expression of miRNA transcriptome in two wheat cultivars, Arapahoe (susceptible) and Mace (temperature-sensitive-resistant), that were mock-inoculated or inoculated with WSMV, TriMV, or both at 18 °C and 27 °C.

RESULTS

Our results showed that the most abundant miRNA family among all the treatments was miRNA166, followed by 159a and 168a, although the order of the latter two changed depending on the infections. When comparing infected and control groups, twenty miRNAs showed significant upregulation, while eight miRNAs were significantly downregulated. Among them, miRNAs 9670-3p, 397-5p, and 5384-3p exhibited the most significant upregulation, whereas miRNAs 319, 9773, and 9774 were the most downregulated. The comparison of infection versus the control group for the cultivar Mace showed temperature-dependent regulation of these miRNAs. The principal component analysis confirmed that less abundant miRNAs among differentially expressed miRNAs were strongly correlated with the inoculated symptomatic wheat cultivars. Notably, miRNAs 397-5p, 398, and 9670-3p were upregulated in response to WSMV and TriMV infections, an observation not yet reported in this context. The significant upregulation of these three miRNAs was further confirmed with RT-qPCR analysis; in general, the RT-qPCR results were in agreement with our computational analysis. Target prediction analysis showed that the miRNAs standing out in our analysis targeted genes involved in defense response and regulation of transcription.

CONCLUSION

Investigation into the roles of these miRNAs and their corresponding targets holds promise for advancing our understanding of the mechanisms of virus infection and possible manipulation of these factors for developing durable virus resistance in crop plants.

Combined Omics Approaches Reveal Distinct Mechanisms of Resistance and/or Susceptibility in Sugar Beet Double Haploid Genotypes at Early Stages of Beet Curly Top Virus Infection

Sugar beet is susceptible to Beet curly top virus (BCTV), which significantly reduces yield and sugar production in the semi-arid growing regions worldwide. Sources of genetic resistance to BCTV is limited and control depends upon insecticide seed treatments with neonicotinoids. Through double haploid production and genetic selection, BCTV resistant breeding lines have been developed. Using BCTV resistant (R) [KDH13; Line 13 and KDH4-9; Line 4] and susceptible (S) [KDH19-17; Line 19] lines, beet leafhopper mediated natural infection, mRNA/sRNA sequencing, and metabolite analyses, potential mechanisms of resistance against the virus and vector were identified. At early infection stages (2- and 6-days post inoculation), examples of differentially expressed genes highly up-regulated in the 'R' lines (vs. 'S') included EL10Ac5g10437 (inhibitor of trypsin and hageman factor), EL10Ac6g14635 (jasmonate-induced protein), EL10Ac3g06016 (ribosome related), EL10Ac2g02812 (probable prolyl 4-hydroxylase 10), etc. Pathway enrichment analysis showed differentially expressed genes were predominantly involved with peroxisome, amino acids metabolism, fatty acid degradation, amino/nucleotide sugar metabolism, etc. Metabolite analysis revealed significantly higher amounts of specific isoflavonoid O-glycosides, flavonoid 8-C glycosides, triterpenoid, and iridoid-O-glycosides in the leaves of the 'R' lines (vs. 'S'). These data suggest that a combination of transcriptional regulation and production of putative antiviral metabolites might contribute to BCTV resistance. In addition, genome divergence among BCTV strains differentially affects the production of small non-coding RNAs (sncRNAs) and small peptides which may potentially affect pathogenicity and disease symptom development.

Argonaute 5-mediated antiviral defense and viral counter-defense in Nicotiana benthamiana

The argonaute (AGO) family proteins play a crucial role in preventing viral invasions through the plant antiviral RNA silencing pathway, with distinct AGO proteins recruited for specific antiviral mechanisms. Our previous study revealed that Nicotiana benthamiana AGO5 (NbAGO5) expression was significantly upregulated in response to bamboo mosaic virus (BaMV) infection. However, the roles of NbAGO5 in antiviral mechanisms remained to be explored. In this research, we examined the antiviral functions of NbAGO5 in the infections of different viruses. It was found that the accumulation of NbAGO5 was induced not only at the RNA but also at the protein level following the infections of BaMV, potato virus X (PVX), tobacco mosaic virus (TMV), and cucumber mosaic virus (CMV) in N. benthamiana. To explore the antiviral mechanism and regulatory function of NbAGO5, we generated NbAGO5 overexpression (OE-NbAGO5) and knockout (nbago5) transgenic N. benthamiana lines. Our findings reveal that NbAGO5 provides defense against BaMV, PVX, TMV, and a mutant CMV deficient in 2b gene, but not against the wild-type CMV and turnip mosaic virus (TuMV). Through affinity purification and small RNA northern blotting, we demonstrated that NbAGO5 exerts its antiviral function by binding to viral small interfering RNAs (vsiRNAs). Moreover, we observed that CMV 2b and TuMV HC-Pro interact with NbAGO5, triggering its degradation via the 26S proteasome and autophagy pathways, thereby allowing these viruses to overcome NbAGO5-mediated defense. In addition, TuMV HC-Pro provides another line of counter-defense by interfering with vsiRNA binding by NbAGO5. Our study provides further insights into the antiviral RNA interference mechanism and the complex interplay between NbAGO5 and plant viruses.

Role of Plant Virus Movement Proteins in Suppression of Host RNAi Defense

One of the systems of plant defense against viral infection is RNA silencing, or RNA interference (RNAi), in which small RNAs derived from viral genomic RNAs and/or mRNAs serve as guides to target an Argonaute nuclease (AGO) to virus-specific RNAs. Complementary base pairing between the small interfering RNA incorporated into the AGO-based protein complex and viral RNA results in the target cleavage or translational repression. As a counter-defensive strategy, viruses have evolved to acquire viral silencing suppressors (VSRs) to inhibit the host plant RNAi pathway. Plant virus VSR proteins use multiple mechanisms to inhibit silencing. VSRs are often multifunctional proteins that perform additional functions in the virus infection cycle, particularly, cell-to-cell movement, genome encapsidation, or replication. This paper summarizes the available data on the proteins with dual VSR/movement protein activity used by plant viruses of nine orders to override the protective silencing response and reviews the different molecular mechanisms employed by these proteins to suppress RNAi.

Evaluation of reference genes for miRNA and mRNA normalization in tobacco infected with PVYNTN, PVYN-Wi and PVYZ-NTN strains

This is the first report on identification of the most suitable reference genes for RT-qPCR quantification of miRNA and mRNA in tobacco response to the prevalent recombinant potato virus Y (PVY) strains PVY^NTN, PVY^N-Wi and the newly identified PVY^Z-NTN. Of 10 tested genes, the expression levels of neIF5C, nU2af and nPP2A were the most stable in samples taken from non-inoculated, mock-inoculated, and infected plants at three days post-inoculation (dpi) and 14 dpi. While the homologues of eIF5 were most stably expressed in tobacco in this study and in potato in our previous study (Yin et al., 2021) following inoculation with the same three PVY strains, the homologues of other two genes PP2A and U2af were stably expressed only in tobacco but unstable in potato. The tobacco homologue of PP2A, which was the most stably expressed one in tobacco interaction with PVY^NTN, PVY^N-Wi and PVY^Z-NTN strains in this study, was the least stable one in tobacco interaction with the non-recombinant PVY^O strain in a previous study (Baek et al., 2017). This study provides evidence on the influence of host species on expression of housekeeping genes and points out virus strain as a new factor influencing expression stability of reference gene. Caution should be taken when choosing reference genes in gene expression study in Solanaceae hosts response to different PVY strains.

Small Talk: On the Possible Role of Trans-Kingdom Small RNAs during Plant-Virus-Vector Tritrophic Communication

Small RNAs (sRNAs) are the hallmark and main effectors of RNA silencing and therefore are involved in major biological processes in plants, such as regulation of gene expression, antiviral defense, and plant genome integrity. The mechanisms of sRNA amplification as well as their mobile nature and rapid generation suggest sRNAs as potential key modulators of intercellular and interspecies communication in plant-pathogen-pest interactions. Plant endogenous sRNAs can act in cis to regulate plant innate immunity against pathogens, or in trans to silence pathogens' messenger RNAs (mRNAs) and impair virulence. Likewise, pathogen-derived sRNAs can act in cis to regulate expression of their own genes and increase virulence towards a plant host, or in trans to silence plant mRNAs and interfere with host defense. In plant viral diseases, virus infection alters the composition and abundance of sRNAs in plant cells, not only by triggering and interfering with the plant RNA silencing antiviral response, which accumulates virus-derived small interfering RNAs (vsiRNAs), but also by modulating plant endogenous sRNAs. Here, we review the current knowledge on the nature and activity of virus-responsive sRNAs during virus-plant interactions and discuss their role in trans-kingdom modulation of virus vectors for the benefit of virus dissemination.

In Silico Identification of Cassava Genome-Encoded MicroRNAs with Predicted Potential for Targeting the ICMV-Kerala Begomoviral Pathogen of Cassava

Cassava mosaic disease (CMD) is caused by several divergent species belonging to the genus Begomovirus (Geminiviridae) transmitted by the whitefly Bemisia tabaci cryptic species group. In India and other parts of Asia, the Indian cassava mosaic virus-Kerala (ICMV-Ker) is an emergent begomovirus of cassava causing damage that results in reduced yield loss and tuber quality. Double-stranded RNA-mediated interference (RNAi) is an evolutionary conserved mechanism in eukaryotes and highly effective, innate defense system to inhibit plant viral replication and/or translation. The objective of this study was to identify and characterize cassava genome-encoded microRNAs (mes-miRNA) that are predicted to target ICMV-Ker ssDNA-encoded mRNAs, based on four in silico algorithms: miRanda, RNA22, Tapirhybrid, and psRNA. The goal is to deploy the predicted miRNAs to trigger RNAi and develop cassava plants with resistance to ICMV-Ker. Experimentally validated mature cassava miRNA sequences (n = 175) were downloaded from the miRBase biological database and aligned with the ICMV-Ker genome. The miRNAs were evaluated for base-pairing with the cassava miRNA seed regions and to complementary binding sites within target viral mRNAs. Among the 175 locus-derived mes-miRNAs evaluated, one cassava miRNA homolog, mes-miR1446a, was identified to have a predicted miRNA target binding site, at position 2053 of the ICMV-Ker genome. To predict whether the cassava miRNA might bind predicted ICMV-Ker mRNA target(s) that could disrupt viral infection of cassava plants, a cassava locus-derived miRNA-mRNA regulatory network was constructed using Circos software. The in silico-predicted cassava locus-derived mes-miRNA-mRNA network corroborated interactions between cassava mature miRNAs and the ICMV-Ker genome that warrant in vivo analysis, which could lead to the development of ICMV-Ker resistant cassava plants.

Implementation of GA-VirReport, a Web-Based Bioinformatics Toolkit for Post-Entry Quarantine Screening of Virus and Viroids in Plants

High-throughput sequencing (HTS) of host plant small RNA (sRNA) is a popular approach for plant virus and viroid detection. The major bottlenecks for implementing this approach in routine virus screening of plants in quarantine include lack of computational resources and/or expertise in command-line environments and limited availability of curated plant virus and viroid databases. We developed: (1) virus and viroid report web-based bioinformatics workflows on Galaxy Australia called GA-VirReport and GA-VirReport-Stats for detecting viruses and viroids from host plant sRNA extracts and (2) a curated higher plant virus and viroid database (PVirDB). We implemented sRNA sequencing with unique dual indexing on a set of plants with known viruses. Sequencing data were analyzed using GA-VirReport and PVirDB to validate these resources. We detected all known viruses in this pilot study with no cross-sample contamination. We then conducted a large-scale diagnosis of 105 imported plants processed at the post-entry quarantine facility (PEQ), Australia. We detected various pathogens in 14 imported plants and discovered that de novo assembly using 21–22 nt sRNA fraction and the megablast algorithm yielded better sensitivity and specificity. This study reports the successful, large-scale implementation of HTS and a user-friendly bioinformatics workflow for virus and viroid screening of imported plants at the PEQ.