The Functions of RNA-Dependent RNA Polymerases in Arabidopsis

The biogenesis, regulation and functions of transitive siRNA in plants

Small RNA (sRNA)-mediated RNA interference (RNAi) is a sequence-specific gene silencing mechanism that modulates gene expression in eukaryotes. As core molecules of RNAi, various sRNAs are encoded in the plant genome or derived from invading RNA molecules, and their biogenesis depends on distinct genetic pathways. Transitive small interfering RNAs (siRNAs), which are sRNAs produced from double-strand RNA (dsRNA) in a process that depends on RNA-dependent RNA polymerases (RDRs), can amplify and spread silencing signals to additional transcripts, thereby enabling a phenomenon termed "transitive RNAi". Members of this class of siRNAs function in various biological processes ranging from development to stress adaptation. In Arabidopsis thaliana, two RDRs participate in the generation of transitive siRNAs, acting cooperatively with various siRNA generation-related factors, such as the RNA-induced silencing complex (RISC) and aberrant RNAs. Transitive siRNAs are produced in diverse subcellular locations and structures under the control of various mechanisms, highlighting the intricacies of their biogenesis and functions. In this review, we discuss recent advances in understanding the molecular events of transitive siRNA biogenesis and its regulation, with a particular focus on factors involved in RDR recruitment. We aim to provide a comprehensive description of the generalized mechanism governing the biogenesis of transitive siRNAs. Additionally, we present an overview of the diverse biological functions of these siRNAs and raise some pressing questions in this area for further investigation.

Cy-1, a major QTL for tomato leaf curl New Delhi virus resistance, harbors a gene encoding a DFDGD-Class RNA-dependent RNA polymerase in cucumber (Cucumis sativus)

BACKGROUND

Tomato leaf curl New Delhi virus (ToLCNDV) (family Geminiviridae, genus Begomovirus) is a significant threat to cucumber (Cucumis sativus) production in many regions. Previous studies have reported the genetic mapping of loci related to ToLCNDV resistance, but no resistance genes have been identified.

RESULTS

We conducted map-based cloning of the ToLCNDV resistance gene in cucumber accession No.44. Agroinfiltration and graft-inoculation analyses confirmed the resistance of No.44 to ToLCNDV isolates from the Mediterranean and Asian countries. Initial mapping involving two rounds of phenotyping with two independent F2 populations generated by crossing the begomovirus-susceptible cultivar SHF and No.44 consistently detected major quantitative trait loci (QTLs) on chromosomes 1 and 2 that confer resistance to ToLCNDV. Fine-mapping of Cy-1, the dominant QTL on chromosome 1, using F3 populations narrowed the candidate region to a 209-kb genomic segment harboring 24 predicted genes. Among these genes, DFDGD-class RNA-dependent RNA polymerase (CsRDR3), an ortholog of Ty-1/Ty-3 of tomato and Pepy-2 of capsicum, was found to be a strong candidate conferring ToLCNDV resistance. The CsRDR3 sequence of No.44 contained multiple amino acid substitutions; the promoter region of CsRDR3 in No.44 had a large deletion; and the CsRDR3 transcript levels were greater in No.44 than in SHF. Virus-induced gene silencing (VIGS) of CsRDR3 using two chromosome segment substitution lines harboring chromosome 1 segments derived from No.44 compromised resistance to ToLCNDV.

CONCLUSIONS

Forward and reverse genetic approaches identified CsRDR3, which encodes a DFDGD-class RNA-dependent RNA polymerase, as the gene responsible for ToLCNDV resistance at the major QTL Cy-1 on chromosome 1 in cucumber. Marker-assisted breeding of ToLCNDV resistance in cucumber will be expedited by using No.44 and the DNA markers developed in this study.

Silencing of RDR1 and RDR6 genes by a single RNAi enhances lettuce's capacity to express recombinant proteins in transient assays

KEY MESSAGE

Enhanced recombinant protein expression was achieved in Salinas lettuce and commercial lettuce by designing a unique RNAi that knockdown the gene-silencing mechanism in transient assays. Improved yields of recombinant proteins (RP) are necessary for protein-production efficiency and ease of purification. Achieving high yield in non-tobacco plants will enable diverse plants to be used as hosts in transient protein-expression systems. With improved protein yield, lettuce (Lactuca sativa) could take the lead as a plant host for RP production. Therefore, this study aimed to improve RP production in lettuce var. Salinas by designing a single RNA interference (RNAi) construct targeting LsRDR1 and LsRDR6 using the Tsukuba system vector. Two RNAi constructs, RNAi-1 and RNAi-2, targeting common regions of LsRDR1 and LsRDR6 with 75% and 76% similarity, respectively, were employed to evaluate simultaneous gene silencing. Quantitative transcription analysis demonstrated that both RNAi constructs effectively knocked down LsRDR6 and LsRDR1, but not LsRDR2, at both 3 and 5 days post-infiltration (dpi), with RNAi-1 exhibited slightly higher efficiency. Based on the protein yield, co-expression of RNAi-1 with enhanced green fluorescent protein (EGFP) increased EGFP expression by approximately 4.9-fold and 3.7-fold at 3 dpi and 5 dpi, respectively, compared to control. A similar but slightly lower increase (2.4-fold and 2.33-fold) was observed in commercial lettuce at 3 and 5 dpi, respectively. To confirm these results, co-infiltration with Bet v 1, a major allergen from birch pollen, resulted in a 2.5-fold increase in expression in Salinas lettuce at 5 dpi. This study marks a significant advancement in enhancing transient protein production in lettuce, elevating its potential as a host for recombinant protein production.

Genetic Dissection of ToLCNDV Resistance in Resistant Sources of Cucumis melo

Tomato leaf curl New Delhi virus (ToLCNDV) is a begomovirus causing significant melon (Cucumis melo) crop losses globally. This study aims to map the ToLCNDV resistance in the PI 414723 melon accession, previously identified and characterized through phenotypic studies, thereby exploring shared genomic regions with the established resistant source WM-7. In the present study, WM-7 and PI 414723 were crossed with the susceptible accessions 'Rochet' and 'Blanco' respectively, to generate F1 hybrids. These hybrids were self-pollinated to generate the populations for mapping the ToLCNDV resistance region and designing markers for marker-assisted selection. Disease evaluation included visual symptom scoring, viral-load quantification and tissue printing. Genotyping-by-sequencing and SNP markers were used for quantitative trait loci (QTL) mapping. For genetic analysis, qPCR and bulked segregant RNA-seq (BSR-seq) were performed. Gene expression was assessed using RNA-seq, and qRT-PCR was used for confirmation. The research narrows the candidate region for resistance in WM-7 and identifies overlapping QTLs on chromosome 11 in PI 414723, found in the region of the DNA primase large subunit. BSR-seq and expression analyses highlight potential regulatory roles of chromosome 2 in conferring resistance. Differential expression was confirmed for six genes in the candidate region on chromosome 2. This study confirms the existence of common resistance genes in PI 414723 and WM-7.

Evidence for an RNAi-independent role of Arabidopsis DICER-LIKE2 in growth inhibition and basal antiviral resistance

Flowering plant genomes encode four or five DICER-LIKE (DCL) enzymes that produce small interfering RNAs (siRNAs) and microRNAs, which function in RNA interference (RNAi). Different RNAi pathways in plants effect transposon silencing, antiviral defense, and endogenous gene regulation. DCL2 acts genetically redundantly with DCL4 to confer basal antiviral defense. However, DCL2 may also counteract DCL4 since knockout of DCL4 causes growth defects that are suppressed by DCL2 inactivation. Current models maintain that RNAi via DCL2-dependent siRNAs is the biochemical basis of both effects. Here, we report that DCL2-mediated antiviral resistance and growth defects cannot be explained by the silencing effects of DCL2-dependent siRNAs. Both functions are defective in genetic backgrounds that maintain high levels of DCL2-dependent siRNAs, either with specific point mutations in DCL2 or with reduced DCL2 dosage because of heterozygosity for dcl2 knockout alleles. Intriguingly, all DCL2 functions require its catalytic activity, and the penetrance of DCL2-dependent growth phenotypes in dcl4 mutants correlates with DCL2 protein levels but not with levels of major DCL2-dependent siRNAs. We discuss this requirement and correlation with catalytic activity but not with resulting siRNAs, in light of other findings that reveal a DCL2 function in innate immunity activation triggered by cytoplasmic double-stranded RNA.

Genome-wide identification and in silico characterization of major RNAi gene families in date palm (Phoenix dactylifera)

BACKGROUND

Dates contain various minerals that are essential for good health. The major RNA interference (RNAi) gene families play a vital role in plant growth and development by controlling the expression of protein-coding genes against different biotic and abiotic stresses. However, these gene families for date palm are not yet studied. Therefore, this study has explored major RNAi genes and their characteristics in date palm.

RESULTS

We have identified 4 PdDCLs, 7 PdAGOs, and 3 PdRDRs as RNAi proteins from the date palm genome by using AtRNAi genes as query sequences in BLASTp search. Domain analysis of predicted RNAi genes has revealed the Helicase_C, Dicer_dimer, PAZ, RNase III, and Piwi domains that are associated with the gene silencing mechanisms. Most PdRNAi proteins have been found in the nucleus and cytosol associated with the gene silencing actions. The gene ontology (GO) enrichment analysis has revealed some important GO terms including RNA interference, dsRNA fragmentation, and ribonuclease_III activity that are related to the protein-coding gene silencing mechanisms. Gene regulatory network (GRN) analysis has identified PAZ and SNF2 as the transcriptional regulators of PdRNAi genes. Top-ranked 10 microRNAs including Pda-miR156b, Pda-miR396a, Pda-miR166a, Pda-miR167d, and Pda-miR529a have been identified as the key post-transcriptional regulators of PdRNAi genes that are associated with different biotic/abiotic stresses. The cis-acting regulatory element analysis of PdRNAi genes has detected some vital cis-acting elements including ABRE, MBS, MYB, MYC, Box-4, G-box, I-box, and STRE that are linked with different abiotic stresses.

CONCLUSION

The results of this study might be valuable resources for the improvement of different characteristics in date palm by further studies in wet-lab.

AT Hook-Like 10 phosphorylation determines ribosomal RNA processing 6-like 1 (RRP6L1) chromatin association and growth suppression during water stress

Phosphorylation of AT Hook-Like 10 (AHL10), one of the AT-hook family of plant-specific DNA binding proteins, is critical for growth suppression during moderate severity drought (low water potential, ψw ) stress. To understand how AHL10 phosphorylation determines drought response, we identified putative AHL10 interacting proteins and further characterized interaction with RRP6L1, a protein involved in epigenetic regulation. RRP6L1 and AHL10 mutants, as well as ahl10-1rrp6l1-2, had similar phenotypes of increased growth maintenance during low ψw . Chromatin precipitation demonstrated that RRP6L1 chromatin association increased during low ψw stress and was dependent upon AHL10 phosphorylation. Transcriptome analyses showed that AHL10 and RRP6L1 have concordant effects on expression of stress- and development-related genes. Together these results indicate that stress signalling can act via AHL10 phosphorylation to control the chromatin association of the key regulatory protein RRP6L1. AHL10 and RRP6L1 interaction in meristem cells is part of a mechanism to downregulate growth during low ψw stress. Interestingly, the loss of AHL13, which is homologous to AHL10 and phosphorylated at a similar C-terminal site, blocked the enhanced growth maintenance of ahl10-1. Thus, AHL10 and AHL13, despite their close homology, are not redundant but rather have distinct roles, likely related to the formation of AHL hetero-complexes.

Phylogenetic analyses of seven protein families refine the evolution of small RNA pathways in green plants

Several protein families participate in the biogenesis and function of small RNAs (sRNAs) in plants. Those with primary roles include Dicer-like (DCL), RNA-dependent RNA polymerase (RDR), and Argonaute (AGO) proteins. Protein families such as double-stranded RNA-binding (DRB), SERRATE (SE), and SUPPRESSION OF SILENCING 3 (SGS3) act as partners of DCL or RDR proteins. Here, we present curated annotations and phylogenetic analyses of seven sRNA pathway protein families performed on 196 species in the Viridiplantae (aka green plants) lineage. Our results suggest that the RDR3 proteins emerged earlier than RDR1/2/6. RDR6 is found in filamentous green algae and all land plants, suggesting that the evolution of RDR6 proteins coincides with the evolution of phased small interfering RNAs (siRNAs). We traced the origin of the 24-nt reproductive phased siRNA-associated DCL5 protein back to the American sweet flag (Acorus americanus), the earliest diverged, extant monocot species. Our analyses of AGOs identified multiple duplication events of AGO genes that were lost, retained, or further duplicated in subgroups, indicating that the evolution of AGOs is complex in monocots. The results also refine the evolution of several clades of AGO proteins, such as AGO4, AGO6, AGO17, and AGO18. Analyses of nuclear localization signal sequences and catalytic triads of AGO proteins shed light on the regulatory roles of diverse AGOs. Collectively, this work generates a curated and evolutionarily coherent annotation for gene families involved in plant sRNA biogenesis/function and provides insights into the evolution of major sRNA pathways.

AGO2a but not AGO2b mediates antiviral defense against infection of wild-type cucumber mosaic virus in tomato

Evolutionarily conserved antiviral RNA interference (RNAi) mediates a primary antiviral innate immunity preventing infection of broad-spectrum viruses in plants. However, the detailed mechanism in plants is still largely unknown, especially in important agricultural crops, including tomato. Varieties of pathogenic viruses evolve to possess viral suppressors of RNA silencing (VSRs) to suppress antiviral RNAi in the host. Due to the prevalence of VSRs, it is still unknown whether antiviral RNAi truly functions to prevent invasion by natural wild-type viruses in plants and animals. In this research, for the first time we applied CRISPR-Cas9 to generate ago2a, ago2b, or ago2ab mutants for two differentiated Solanum lycopersicum AGO2s, key effectors in antiviral RNAi. We found that AGO2a but not AGO2b was significantly induced to inhibit the propagation of not only VSR-deficient Cucumber mosaic virus (CMV) but also wild-type CMV-Fny in tomato; however, neither AGO2a nor AGO2b regulated disease induction after infection with either virus. Our findings firstly reveal a prominent role of AGO2a in antiviral RNAi innate immunity in tomato and demonstrate that antiviral RNAi evolves to defend against infection of natural wild-type CMV-Fny in tomato. However, AGO2a-mediated antiviral RNAi does not play major roles in promoting tolerance of tomato plants to CMV infection for maintaining health.

Plant terminators: the unsung heroes of gene expression

To be properly expressed, genes need to be accompanied by a terminator, a region downstream of the coding sequence that contains the information necessary for the maturation of the mRNA 3' end. The main event in this process is the addition of a poly(A) tail at the 3' end of the new transcript, a critical step in mRNA biology that has important consequences for the expression of genes. Here, we review the mechanism leading to cleavage and polyadenylation of newly transcribed mRNAs and how this process can affect the final levels of gene expression. We give special attention to an aspect often overlooked, the effect that different terminators can have on the expression of genes. We also discuss some exciting findings connecting the choice of terminator to the biogenesis of small RNAs, which are a central part of one of the most important mechanisms of regulation of gene expression in plants.

Genome-wide identification, characterization and expression analysis of AGO, DCL, and RDR families in Chenopodium quinoa

RNA interference is a highly conserved mechanism wherein several types of non-coding small RNAs regulate gene expression at the transcriptional or post-transcriptional level, modulating plant growth, development, antiviral defence, and stress responses. Argonaute (AGO), DCL (Dicer-like), and RNA-dependent RNA polymerase (RDR) are key proteins in this process. Here, these three protein families were identified in Chenopodium quinoa. Further, their phylogenetic relationships with Arabidopsis, their domains, three-dimensional structure modelling, subcellular localization, and functional annotation and expression were analysed. Whole-genome sequence analysis predicted 21 CqAGO, eight CqDCL, and 11 CqRDR genes in quinoa. All three protein families clustered into phylogenetic clades corresponding to those of Arabidopsis, including three AGO clades, four DCL clades, and four RDR clades, suggesting evolutionary conservation. Domain and protein structure analyses of the three gene families showed almost complete homogeneity among members of the same group. Gene ontology annotation revealed that the predicted gene families might be directly involved in RNAi and other important pathways. Largely, these gene families showed significant tissue-specific expression patterns, RNA-sequencing (RNA-seq) data revealed that 20 CqAGO, seven CqDCL, and ten CqRDR genes tended to have preferential expression in inflorescences. Most of them being downregulated in response to drought, cold, salt and low phosphate stress. To our knowledge, this is the first study to elucidate these key protein families involved in the RNAi pathway in quinoa, which are significant for understanding the mechanisms underlying stress responses in this plant.

Deciphering the genomic regions governing major agronomic traits and blast resistance using genome wide association mapping in finger millet

Finger millet (Eleusine coracana L.) is climate resilient minor millet of Asia and Africa with wide adaptation and unparallel nutritional profile. To date, genomic resources available in finger millet are scanty and genetic control of agronomic traits remains elusive. Here, a collection of eco-geographically diverse 186 genotypes was quantified for variation in 13 agronomic traits and reaction to blast to identify marker-trait associations (MTAs) using genotyping-by-sequencing (GBS) and genome-wide association study (GWAS). GBS generated 2977 high quality single nucleotide polymorphism (SNPs) markers and identified three subpopulations with varying admixture levels. General linear and mixed model approaches of GWAS to correct for population structure and genetic relatedness identified 132 common MTAs for agronomic traits across the years. The phenotypic variance explained by the makers varied from 4.8% (TP692389-flag leaf width) to 20% (TP714446-green fodder weight). Of these, 26 MTAs showed homology with candidate genes having role in plant growth, development and photosynthesis in the genomes of foxtail millet, rice, maize, wheat and barley. We also found 4 common MTAs for neck blast resistance, which explained 5.9-15.1% phenotypic variance. Three MTAs for neck blast resistance showed orthologues in related genera having putative functions in pathogen defense in plants. The results of this work lay a foundation for understanding the genetic architecture of agronomic traits and blast resistance in finger millet and provide a framework for genomics assisted breeding.

Genome-Wide Identification and Characterization of Argonaute, Dicer-like and RNA-Dependent RNA Polymerase Gene Families and Their Expression Analyses in Fragaria spp

In the growth and development of plants, some non-coding small RNAs (sRNAs) not only mediate RNA interference at the post-transcriptional level, but also play an important regulatory role in chromatin modification at the transcriptional level. In these processes, the protein factors Argonaute (AGO), Dicer-like (DCL), and RNA-dependent RNA polymerase (RDR) play very important roles in the synthesis of sRNAs respectively. Though they have been identified in many plants, the information about these gene families in strawberry was poorly understood. In this study, using a genome-wide analysis and a phylogenetic approach, 13 AGO, six DCL, and nine RDR genes were identified in diploid strawberry Fragaria vesca. We also identified 33 AGO, 18 DCL, and 28 RDR genes in octoploid strawberry Fragaria × ananassa, studied the expression patterns of these genes in various tissues and developmental stages of strawberry, and researched the response of these genes to some hormones, finding that almost all genes respond to the five hormone stresses. This study is the first report of a genome-wide analysis of AGO, DCL, and RDR gene families in Fragaria spp., in which we provide basic genomic information and expression patterns for these genes. Additionally, this study provides a basis for further research on the functions of these genes and some evidence for the evolution between diploid and octoploid strawberries.

Transcriptome analyses unveiled differential regulation of AGO and DCL genes by pepino mosaic virus strains

Pepino mosaic virus (PepMV) is a single-stranded (ss), positive-sense (+) RNA potexvirus that affects tomato crops worldwide. We have described an in planta antagonistic interaction between PepMV isolates of two strains in which the EU isolate represses the accumulation of the CH2 isolate during mixed infections. Reports describing transcriptomic responses to mixed infections are scant. We carried out transcriptomic analyses of tomato plants singly and mixed-infected with two PepMV isolates of both strains. Comparison of the transcriptomes of singly infected plants showed that deeper transcriptomic alterations occurred at early infection times, and also that each of the viral strains modulated the host transcriptome differentially. Mixed infections caused transcriptomic alterations similar to those for the sum of single infections at early infection times, but clearly differing at later times postinfection. We next tested the hypothesis that PepMV-EU, in either single or mixed infections, deregulates host gene expression differentially so that virus accumulation of both strains gets repressed. That seemed to be the case for the genes AGO1a, DCL2d, AGO2a, and DCL2b, which are involved in the antiviral silencing pathway and were upregulated by PepMV-EU but not by PepMV-CH2 at early times postinfection. The pattern of AGO2a expression was validated by reverse transcription-quantitative PCR in tomato and Nicotiana benthamiana plants. Using an N. benthamiana ago2 mutant line, we showed that AGO2 indeed plays an important role in the antiviral defence against PepMV, but it is not the primary determinant of the outcome of the antagonistic interaction between the two PepMV strains.

Expansion and contraction of small RNA and methylation machinery throughout plant evolution

The revolution in sequencing has created a wealth of plant genomes that can be mined to understand the evolution of biological complexity. Complexity is often driven by gene duplication, which allows paralogs to specialize in an activity of the ancestral gene or acquire novel functions. Angiosperms encode a variety of gene silencing pathways that share related machinery for small RNA biosynthesis and function. Recent phylogenetic analysis of these gene families plots the expansion, specialization, and occasional contraction of this core machinery. This analysis reveals the ancient origin of RNA-directed DNA Methylation in early land plants, or possibly their algal ancestors, as well as ongoing duplications that evolve novel small RNA pathways.

Concepts and considerations for enhancing RNAi efficiency in phytopathogenic fungi for RNAi-based crop protection using nanocarrier-mediated dsRNA delivery systems

Existing, emerging, and reemerging strains of phytopathogenic fungi pose a significant threat to agricultural productivity globally. This risk is further exacerbated by the lack of resistance source(s) in plants or a breakdown of resistance by pathogens through co-evolution. In recent years, attenuation of essential pathogen gene(s) via double-stranded (ds) RNA-mediated RNA interference (RNAi) in host plants, a phenomenon known as host-induced gene silencing, has gained significant attention as a way to combat pathogen attack. Yet, due to biosafety concerns regarding transgenics, country-specific GMO legislation has limited the practical application of desirable attributes in plants. The topical application of dsRNA/siRNA targeting essential fungal gene(s) through spray-induced gene silencing (SIGS) on host plants has opened up a transgene-free avenue for crop protection. However, several factors influence the outcome of RNAi, including but not limited to RNAi mechanism in plant/fungi, dsRNA/siRNA uptake efficiency, dsRNA/siRNA design parameters, dsRNA stability and delivery strategy, off-target effects, etc. This review emphasizes the significance of these factors and suggests appropriate measures to consider while designing in silico and in vitro experiments for successful RNAi in open-field conditions. We also highlight prospective nanoparticles as smart delivery vehicles for deploying RNAi molecules in plant systems for long-term crop protection and ecosystem compatibility. Lastly, we provide specific directions for future investigations that focus on blending nanotechnology and RNAi-based fungal control for practical applications.

PepYLCIV and PepYLCAV resistance gene Pepy-2 encodes DFDGD-Class RNA-dependent RNA polymerase in Capsicum

A begomovirus resistance gene Pepy-2 encoding the DFDGD-Class RNA-dependent RNA polymerase 3a was identified in pepper (C. annuum) through the forward and reverse genetic analyses. In several countries throughout the world, the whitefly-transmitted begomovirus causes massive yield losses in pepper (Capsicum spp.) production. Although introgression of the genetic resistance against begomovirus to commercial cultivars is strongly required, the recently discovered recessive resistance gene pepy-1, which encodes the messenger RNA surveillance factor Pelota, is the only begomovirus resistance gene identified in Capsicum so far. In this study, we fine-mapped another begomovirus resistance gene from PG1-1 (C. annuum), which is resistant to pepper yellow leaf curl Indonesia virus (PepYLCIV) and pepper yellow leaf curl Aceh virus (PepYLCAV), to further speed up the marker-assisted breeding of begomovirus resistance in peppers. A single dominant locus, Pepy-2, conferring resistance against PepYLCIV in PG1-1 was identified on chromosome 7 by screening recombinants from the F₂ and F₃ segregating populations derived from a cross between PG1-1 and begomovirus susceptible SCM334. In the target region spanning 722 kb, a strong candidate gene, the RNA-dependent RNA polymerase 3a (CaRDR3a), was identified. The whole-genome and transcriptome sequences of PG1-1 and SCM334 revealed a single Guanine (G) deletion in CaRDR3a first exon, causing a frameshift resulting in loss-of-function in SCM334. In addition, multiple loss-of-function alleles of CaRDR3a were identified in the reference sequences of C. annuum, C. chinense, and C. baccatum in the public database. Furthermore, virus-induced gene silencing of CaRDR3a in PG1-1 resulted in the loss of resistance against PepYLCIV. PG1-1 and the DNA marker developed in this study will be useful to breeders using Pepy-2 in their breeding programs.

Roles of RNA silencing in viral and non-viral plant immunity and in the crosstalk between disease resistance systems

RNA silencing is a well-established antiviral immunity system in plants, in which small RNAs guide Argonaute proteins to targets in viral RNA or DNA, resulting in virus repression. Virus-encoded suppressors of silencing counteract this defence system. In this Review, we discuss recent findings about antiviral RNA silencing, including the movement of RNA through plasmodesmata and the differentiation between plant self and viral RNAs. We also discuss the emerging role of RNA silencing in plant immunity against non-viral pathogens. This immunity is mediated by transkingdom movement of RNA into and out of the infected plant cells in vesicles or as extracellular nucleoproteins and, like antiviral immunity, is influenced by the silencing suppressors encoded in the pathogens' genomes. Another effect of RNA silencing on general immunity involves host-encoded small RNAs, including microRNAs, that regulate NOD-like receptors and defence signalling pathways in the innate immunity system of plants. These RNA silencing pathways form a network of processes with both positive and negative effects on the immune systems of plants.

RNAi-Based Antiviral Innate Immunity in Plants

Multiple antiviral immunities were developed to defend against viral infection in hosts. RNA interference (RNAi)-based antiviral innate immunity is evolutionarily conserved in eukaryotes and plays a vital role against all types of viruses. During the arms race between the host and virus, many viruses evolve viral suppressors of RNA silencing (VSRs) to inhibit antiviral innate immunity. Here, we reviewed the mechanism at different stages in RNAi-based antiviral innate immunity in plants and the counteractions of various VSRs, mainly upon infection of RNA viruses in model plant Arabidopsis. Some critical challenges in the field were also proposed, and we think that further elucidating conserved antiviral innate immunity may convey a broad spectrum of antiviral strategies to prevent viral diseases in the future.

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

Small RNAs are significant regulators of gene expression, which play multiple roles in plant development, growth, reproductive and stress response. It is generally believed that the regulation of plants’ endogenous genes by small RNAs has evolved from a cellular defense mechanism for RNA viruses and transposons. Most small RNAs have well-established roles in the defense response, such as viral response. During viral infection, plant endogenous small RNAs can direct virus resistance by regulating the gene expression in the host defense pathway, while the small RNAs derived from viruses are the core of the conserved and effective RNAi resistance mechanism. As a counter strategy, viruses evolve suppressors of the RNAi pathway to disrupt host plant silencing against viruses. Currently, several studies have been published elucidating the mechanisms by which small RNAs regulate viral defense in different crops. This paper reviews the distinct pathways of small RNAs biogenesis and the molecular mechanisms of small RNAs mediating antiviral immunity in plants, as well as summarizes the coping strategies used by viruses to override this immune response. Finally, we discuss the current development state of the new applications in virus defense based on small RNA silencing.

Small RNA Profiling of Susceptible and Resistant Ty-1 Encoding Tomato Plants Upon Tomato Yellow Leaf Curl Virus Infection

Ty-1 presents an atypical dominant resistance gene that codes for an RNA-dependent RNA polymerase (RDR) of the gamma class and confers resistance to tomato yellow leaf curl virus (TYLCV) and other geminiviruses. Tomato lines bearing Ty-1 not only produce relatively higher amounts of viral small interfering (vsi)RNAs, but viral DNA also exhibits a higher amount of cytosine methylation. Whether Ty-1 specifically enhances posttranscriptional gene silencing (PTGS), leading to a degradation of RNA target molecules and primarily relying on 21-22 nucleotides (nts) siRNAs, and/or transcriptional gene silencing (TGS), leading to the methylation of cytosines within DNA target sequences and relying on 24-nts siRNAs, was unknown. In this study, small RNAs were isolated from systemically TYLCV-infected leaves of Ty-1 encoding tomato plants and susceptible tomato Moneymaker (MM) and sequence analyzed. While in susceptible tomato plants vsiRNAs of the 21-nt size class were predominant, their amount was drastically reduced in tomato containing Ty-1. The latter, instead, revealed elevated levels of vsiRNAs of the 22- and 24-nt size classes. In addition, the genomic distribution profiles of the vsiRNAs were changed in Ty-1 plants compared with those from susceptible MM. In MM three clear hotspots were seen, but these were less pronounced in Ty-1 plants, likely due to enhanced transitive silencing to neighboring viral genomic sequences. The largest increase in the amount of vsiRNAs was observed in the intergenic region and the V1 viral gene. The results suggest that Ty-1 enhances an antiviral TGS response. Whether the elevated levels of 22 nts vsiRNAs contribute to an enhanced PTGS response or an additional TGS response involving a noncanonical pathway of RNA dependent DNA methylation remains to be investigated.

N. S, Pachamuthu K, Chenna S, Nair A, Shivaprasad PV. Essential role of γ-clade RNA-dependent RNA polymerases in rice development and yield-related traits is linked to their atypical polymerase activities regulating specific genomic regions

RNA-dependent RNA polymerases (RDR) generate double-stranded (ds)RNA triggers for RNA silencing across eukaryotes. Among the three clades, α-clade and β-clade members are key components of RNA silencing and mediators of stress responses across eukaryotes. However, γ-clade members are unusual in that they are represented in phylogenetically distant plants and fungi, and their functions are unknown. Using genetic, bioinformatic and biochemical methods, we show that γ-clade RDRs from Oryza sativa L. are involved in plant development as well as regulation of expression of coding and noncoding RNAs. Overexpression of γ-clade RDRs in transgenic rice and tobacco plants resulted in robust growth phenotype, whereas their silencing in rice displayed strong inhibition of growth. Small (s)RNA and RNA-seq analysis of OsRDR3 mis-expression lines suggested that it is specifically involved in the regulation of repeat-rich regions in the genome. Biochemical analysis confirmed that OsRDR3 has robust polymerase activities on both single stranded (ss)RNA and ssDNA templates similar to the activities reported for α-clade RDRs such as AtRDR6. Our results provide the first evidence of the importance of γ-clade RDRs in plant development, their atypical biochemical activities and their contribution to the regulation of gene expression.

Cucumber RDR1s and cucumber mosaic virus suppressor protein 2b association directs host defence in cucumber plants

RNA-dependent RNA polymerases (RDRs) regulate important aspects of plant development and resistance to pathogens. The role of RDRs in virus resistance has been demonstrated using siRNA signal amplification and through the methylation of viral genomes. Cucumber (Cucumis sativus) has four RDR1 genes that are differentially induced during virus infection: CsRDR1a, CsRDR1b, and duplicated CsRDR1c1/c2. The mode of action of CsRDR1s during viral infection is unknown. Transient expression of the cucumber mosaic virus (CMV)-2b protein (the viral suppressor of RNA silencing) in cucumber protoplasts induced the expression of CsRDR1c, but not of CsRDR1a/1b. Results from the yeast two-hybrid system showed that CsRDR1 proteins interacted with CMV-2b and this was confirmed by bimolecular fluorescence complementation assays. In protoplasts, CsRDR1s localized in the cytoplasm as punctate spots. Colocalization experiments revealed that CsRDR1s and CMV-2b were uniformly dispersed throughout the cytoplasm, suggesting that CsRDR1s are redistributed as a result of interactions. Transient overexpression of individual CsRDR1a/1b genes in protoplasts reduced CMV accumulation, indicating their antiviral role. However, overexpression of CsRDR1c in protoplasts resulted in relatively higher accumulation of CMV and CMVΔ2b. In single cells, CsRDR1c enhances viral replication, leading to CMV accumulation and blocking secondary siRNA amplification of CsRDR1c by CMV-2b protein. This suggests that CMV-2b acts as both a transcription factor that induces CsRDR1c (controlling virus accumulation) and a suppressor of CsRDR1c activity.

RNA-dependent RNA polymerase 1 delays the accumulation of viroids in infected plants

RNA-dependent RNA polymerase 1 (RDR1) is essential for plant antiviral defence, but its role in plant defence against viroid infection remains unknown. The present study aimed to identify the function and mechanism of RDR1 in plant resistance to viroid infection. Overexpression of Nicotiana tabacum RDR1 (NtRDR1) delayed the accumulation of potato spindle tuber viroid (PSTVd) genomic RNA and PSTVd-derived small RNA (sRNA) in Nicotiana benthamiana plants at the early invasion stage, but not in the late stage of infection. Conversely, virus-induced gene silencing of tomato RDR1 (SlRDR1a) increased the susceptibility to PSTVd infection (increased viroid accumulation). Salicylic acid (SA) pretreatment induced SlRDR1a expression and enhanced the defence against PSTVd infection in tomato plants. Our study demonstrated that RDR1 is involved in SA-mediated defence and restricts the early systemic invasion by PSTVd in plants. The decreased PSTVd accumulation in N. benthamiana was not caused by efficient accumulation of PSTVd sRNAs. These results deepen our understanding of the mechanism of RDR1 in plant defence responses to viroid attack.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory elements analyses in banana (Musa acuminata)

RNA silencing is mediated through RNA interference (RNAi) pathway gene families, i.e., Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) and their cis-acting regulatory elements. The RNAi pathway is also directly connected with the post-transcriptional gene silencing (PTGS) mechanism, and the pathway controls eukaryotic gene regulation during growth, development, and stress response. Nevertheless, genome-wide identification of RNAi pathway gene families such as DCL, AGO, and RDR and their regulatory network analyses related to transcription factors have not been studied in many fruit crop species, including banana (Musa acuminata). In this study, we studied in silico genome-wide identification and characterization of DCL, AGO, and RDR genes in bananas thoroughly via integrated bioinformatics approaches. A genome-wide analysis identified 3 MaDCL, 13 MaAGO, and 5 MaRDR candidate genes based on multiple sequence alignment and phylogenetic tree related to the RNAi pathway in banana genomes. These genes correspond to the Arabidopsis thaliana RNAi silencing genes. The analysis of the conserved domain, motif, and gene structure (exon-intron numbers) for MaDCL, MaAGO, and MaRDR genes showed higher homogeneity within the same gene family. The Gene Ontology (GO) enrichment analysis exhibited that the identified RNAi genes could be involved in RNA silencing and associated metabolic pathways. A number of important transcription factors (TFs), e.g., ERF, Dof, C2H2, TCP, GATA and MIKC_MADS families, were identified by network and sub-network analyses between TFs and candidate RNAi gene families. Furthermore, the cis-acting regulatory elements related to light-responsive (LR), stress-responsive (SR), hormone-responsive (HR), and other activities (OT) functions were identified in candidate MaDCL, MaAGO, and MaRDR genes. These genome-wide analyses of these RNAi gene families provide valuable information related to RNA silencing, which would shed light on further characterization of RNAi genes, their regulatory elements, and functional roles, which might be helpful for banana improvement in the breeding program.

Atypical DNA methylation, sRNA-size distribution, and female gametogenesis in Utricularia gibba

The most studied DNA methylation pathway in plants is the RNA Directed DNA Methylation (RdDM), a conserved mechanism that involves the role of noncoding RNAs to control the expansion of the noncoding genome. Genome-wide DNA methylation levels have been reported to correlate with genome size. However, little is known about the catalog of noncoding RNAs and the impact on DNA methylation in small plant genomes with reduced noncoding regions. Because of the small length of intergenic regions in the compact genome of the carnivorous plant Utricularia gibba, we investigated its repertoire of noncoding RNA and DNA methylation landscape. Here, we report that, compared to other angiosperms, U. gibba has an unusual distribution of small RNAs and reduced global DNA methylation levels. DNA methylation was determined using a novel strategy based on long-read DNA sequencing with the Pacific Bioscience platform and confirmed by whole-genome bisulfite sequencing. Moreover, some key genes involved in the RdDM pathway may not represented by compensatory paralogs or comprise truncated proteins, for example, U. gibba DICER-LIKE 3 (DCL3), encoding a DICER endonuclease that produces 24-nt small-interfering RNAs, has lost key domains required for complete function. Our results unveil that a truncated DCL3 correlates with a decreased proportion of 24-nt small-interfering RNAs, low DNA methylation levels, and developmental abnormalities during female gametogenesis in U. gibba. Alterations in female gametogenesis are reminiscent of RdDM mutant phenotypes in Arabidopsis thaliana. It would be interesting to further study the biological implications of the DCL3 truncation in U. gibba, as it could represent an initial step in the evolution of RdDM pathway in compact genomes.

Transgenerational effect of mutants in the RNA-directed DNA methylation pathway on the triploid block in Arabidopsis

BACKGROUND

Hybridization of plants that differ in number of chromosome sets (ploidy) frequently causes endosperm failure and seed arrest, a phenomenon referred to as triploid block. In Arabidopsis, loss of function of NRPD1, encoding the largest subunit of the plant-specific RNA polymerase IV (Pol IV), can suppress the triploid block. Pol IV generates short RNAs required to guide de novo methylation in the RNA-directed DNA methylation (RdDM) pathway. Recent work suggests that suppression of the triploid block by mutants in RdDM components differs, depending on whether the diploid pollen is derived from tetraploid plants or from the omission in second division 1 (osd1) mutant. This study aims to understand this difference.

RESULTS

In this study, we find that the ability of mutants in the RdDM pathway to suppress the triploid block depends on their degree of inbreeding. While first homozygous generation mutants in RdDM components NRPD1, RDR2, NRPE1, and DRM2 have weak or no ability to rescue the triploid block, they are able to suppress the triploid block with successive generations of inbreeding. Inbreeding of nrpd1 was connected with a transgenerational loss of non-CG DNA methylation on sites jointly regulated by CHROMOMETHYLASES 2 and 3.

CONCLUSIONS

Our data reveal that loss of RdDM function differs in its effect in early and late generations, which has important implications when interpreting the effect of RdDM mutants.

Global Analysis of RNA-Dependent RNA Polymerase-Dependent Small RNAs Reveals New Substrates and Functions for These Proteins and SGS3 in Arabidopsis

RNA silencing pathways control eukaryotic gene expression transcriptionally or posttranscriptionally in a sequence-specific manner. In RNA silencing, the production of double-stranded RNA (dsRNA) gives rise to various classes of 20-24 nucleotide (nt) small RNAs (smRNAs). In Arabidopsis thaliana, smRNAs are often derived from long dsRNA molecules synthesized by one of the six genomically encoded RNA-dependent RNA Polymerase (RDR) proteins. However, the full complement of the RDR-dependent smRNAs and functions that these proteins and their RNA-binding cofactors play in plant RNA silencing has not been fully uncovered. To address this gap, we performed a global genomic analysis of all six RDRs and two of their cofactors to find new substrates for RDRs and targets of the resulting RDR-derived siRNAs to uncover new functions for these proteins in plants. Based on these analyses, we identified substrates for the three RDRγ clade proteins (RDR3-5), which had not been well-characterized previously. We also identified new substrates for the other three RDRs (RDR1, RDR2, and RDR6) as well as the RDR2 cofactor RNA-directed DNA methylation 12 (RDM12) and the RDR6 cofactor suppressor of gene silencing 3 (SGS3). These findings revealed that the target substrates of SGS3 are not limited to those solely utilized by RDR6, but that this protein seems to be a more general cofactor for the RDR family of proteins. Additionally, we found that RDR6 and SGS3 are involved in the production of smRNAs that target transcripts related to abiotic stresses, including water deprivation, salt stress, and ABA response, and as expected the levels of these mRNAs are increased in rdr6 and sgs3 mutant plants. Correspondingly, plants that lack these proteins (rdr6 and sgs3 mutants) are hypersensitive to ABA treatment, tolerant to high levels of PEG8000, and have a higher survival rate under salt treatment in comparison to wild-type plants. In total, our analyses have provided an extremely data-rich resource for uncovering new functions of RDR-dependent RNA silencing in plants, while also revealing a previously unexplored link between the RDR6/SGS3-dependent pathway and plant abiotic stress responses.

Response of bitter and sweet Chenopodium quinoa varieties to cucumber mosaic virus: Transcriptome and small RNASeq perspective

Saponins are secondary metabolites with antiviral properties. Low saponin (sweet) varieties of quinoa (Chenopodium quinoa) have been developed because seeds high in saponins taste bitter. The aim of this study was to elucidate the role of saponin in resistance of quinoa to Cucumber mosaic virus (CMV). Differential gene expression was studied in time-series study of CMV infection. High-throughput transcriptome sequence data were obtained from 36 samples (3 varieties × +/- CMV × 1 or 4 days after inoculation × 3 replicates). Translation, lipid, nitrogen, amino acid metabolism, and mono- and sesquiterpenoid biosynthesis genes were upregulated in CMV infections. In 'Red Head' (bitter), CMV-induced systemic symptoms were concurrent with downregulation of a key saponin biosynthesis gene, TSARL1, four days after inoculation. In local lesion responses (sweet and semi-sweet), TSARL1 levels remained up-regulated. Known microRNAs (miRNA) (81) from 11 miR families and 876 predicted novel miRNAs were identified. Differentially expressed miRNA and short interfering RNA clusters (24nt) induced by CMV infection are predicted to target genomic and intergenic regions enriched in repetitive elements. This is the first report of integrated RNASeq and sRNASeq data in quinoa-virus interactions and provides comprehensive understanding of involved genes, non-coding regions, and biological pathways in virus resistance.

Virus and Viroid-Derived Small RNAs as Modulators of Host Gene Expression: Molecular Insights Into Pathogenesis

Virus-derived siRNAs (vsiRNAs) generated by the host RNA silencing mechanism are effectors of plant’s defense response and act by targeting the viral RNA and DNA in post-transcriptional gene silencing (PTGS) and transcriptional gene silencing (TGS) pathways, respectively. Contrarily, viral suppressors of RNA silencing (VSRs) compromise the host RNA silencing pathways and also cause disease-associated symptoms. In this backdrop, reports describing the modulation of plant gene(s) expression by vsiRNAs via sequence complementarity between viral small RNAs (sRNAs) and host mRNAs have emerged. In some cases, silencing of host mRNAs by vsiRNAs has been implicated to cause characteristic symptoms of the viral diseases. Similarly, viroid infection results in generation of sRNAs, originating from viroid genomic RNAs, that potentially target host mRNAs causing typical disease-associated symptoms. Pathogen-derived sRNAs have been demonstrated to have the propensity to target wide range of genes including host defense-related genes, genes involved in flowering and reproductive pathways. Recent evidence indicates that vsiRNAs inhibit host RNA silencing to promote viral infection by acting as decoy sRNAs. Nevertheless, it remains unclear if the silencing of host transcripts by viral genome-derived sRNAs are inadvertent effects due to fortuitous pairing between vsiRNA and host mRNA or the result of genuine counter-defense strategy employed by viruses to enhance its survival inside the plant cell. In this review, we analyze the instances of such cross reaction between pathogen-derived vsiRNAs and host mRNAs and discuss the molecular insights regarding the process of pathogenesis.

RNA-Seq Transcriptome Analysis Provides Candidate Genes for Resistance to Tomato Leaf Curl New Delhi Virus in Melon

Tomato leaf curl New Delhi virus (ToLCNDV) emerged in the Mediterranean Basin in 2012 as the first DNA bipartite begomovirus (Geminiviridae family), causing severe yield and economic losses in cucurbit crops. A major resistance locus was identified in the wild melon accession WM-7 (Cucumis melo kachri group), but the mechanisms involved in the resistant response remained unknown. In this work, we used RNA-sequencing to identify disease-associated genes that are differentially expressed in the course of ToLCNDV infection and could contribute to resistance. Transcriptomes of the resistant WM-7 genotype and the susceptible cultivar Piñonet Piel de Sapo (PS) (C. melo ibericus group) in ToLCNDV and mock inoculated plants were compared at four time points during infection (0, 3, 6, and 12 days post inoculation). Different gene expression patterns were observed over time in the resistant and susceptible genotypes in comparison to their respective controls. Differentially expressed genes (DEGs) in ToLCNDV-infected plants were classified using gene ontology (GO) terms, and genes of the categories transcription, DNA replication, and helicase activity were downregulated in WM-7 but upregulated in PS, suggesting that reduced activity of these functions reduces ToLCNDV replication and intercellular spread and thereby contributes to resistance. DEGs involved in the jasmonic acid signaling pathway, photosynthesis, RNA silencing, transmembrane, and sugar transporters entail adverse consequences for systemic infection in the resistant genotype, and lead to susceptibility in PS. The expression levels of selected candidate genes were validated by qRT-PCR to corroborate their differential expression upon ToLCNDV infection in resistant and susceptible melon. Furthermore, single nucleotide polymorphism (SNPs) with an effect on structural functionality of DEGs linked to the main QTLs for ToLCNDV resistance have been identified. The obtained results pinpoint cellular functions and candidate genes that are differentially expressed in a resistant and susceptible melon line in response to ToLCNDV, an information of great relevance for breeding ToLCNDV-resistant melon cultivars.

Secondary siRNAs in Plants: Biosynthesis, Various Functions, and Applications in Virology

The major components of RNA silencing include both transitive and systemic small RNAs, which are technically called secondary sRNAs. Double-stranded RNAs trigger systemic silencing pathways to negatively regulate gene expression. The secondary siRNAs generated as a result of transitive silencing also play a substantial role in gene silencing especially in antiviral defense. In this review, we first describe the discovery and pathways of transitivity with emphasis on RNA-dependent RNA polymerases followed by description on the short range and systemic spread of silencing. We also provide an in-depth view on the various size classes of secondary siRNAs and their different roles in RNA silencing including their categorization based on their biogenesis. The other regulatory roles of secondary siRNAs in transgene silencing, virus-induced gene silencing, transitivity, and trans-species transfer have also been detailed. The possible implications and applications of systemic silencing and the different gene silencing tools developed are also described. The details on mobility and roles of secondary siRNAs derived from viral genome in plant defense against the respective viruses are presented. This entails the description of other compatible plant-virus interactions and the corresponding small RNAs that determine recovery from disease symptoms, exclusion of viruses from shoot meristems, and natural resistance. The last section presents an overview on the usefulness of RNA silencing for management of viral infections in crop plants.

In silico identification and characterization of AGO, DCL and RDR gene families and their associated regulatory elements in sweet orange (Citrus sinensis L.)

RNA interference (RNAi) plays key roles in post-transcriptional and chromatin modification levels as well as regulates various eukaryotic gene expressions which are involved in stress responses, development and maintenance of genome integrity during developmental stages. The whole mechanism of RNAi pathway is directly involved with the gene-silencing process by the interaction of Dicer-Like (DCL), Argonaute (AGO) and RNA-dependent RNA polymerase (RDR) gene families and their regulatory elements. However, these RNAi gene families and their sub-cellular locations, functional pathways and regulatory components were not extensively investigated in the case of economically and nutritionally important fruit plant sweet orange (Citrus sinensis L.). Therefore, in silico characterization, gene diversity and regulatory factor analysis of RNA silencing genes in C. sinensis were conducted by using the integrated bioinformatics approaches. Genome-wide comparison analysis based on phylogenetic tree approach detected 4 CsDCL, 8 CsAGO and 4 CsRDR as RNAi candidate genes in C. sinensis corresponding to the RNAi genes of model plant Arabidopsis thaliana. The domain and motif composition and gene structure analyses for all three gene families exhibited almost homogeneity within the same group members. The Gene Ontology enrichment analysis clearly indicated that the predicted genes have direct involvement into the gene-silencing and other important pathways. The key regulatory transcription factors (TFs) MYB, Dof, ERF, NAC, MIKC_MADS, WRKY and bZIP were identified by their interaction network analysis with the predicted genes. The cis-acting regulatory elements associated with the predicted genes were detected as responsive to light, stress and hormone functions. Furthermore, the expressed sequence tag (EST) analysis showed that these RNAi candidate genes were highly expressed in fruit and leaves indicating their organ specific functions. Our genome-wide comparison and integrated bioinformatics analyses provided some necessary information about sweet orange RNA silencing components that would pave a ground for further investigation of functional mechanism of the predicted genes and their regulatory factors.

Messenger RNA 5' NAD+ Capping Is a Dynamic Regulatory Epitranscriptome Mark That Is Required for Proper Response to Abscisic Acid in Arabidopsis

Although eukaryotic messenger RNAs (mRNAs) normally possess a 5' end N⁷-methyl guanosine (m⁷G) cap, a non-canonical 5' nicotinamide adenine dinucleotide (NAD⁺) cap can tag certain transcripts for degradation mediated by the NAD⁺ decapping enzyme DXO1. Despite this importance, whether NAD⁺ capping dynamically responds to specific stimuli to regulate eukaryotic transcriptomes remains unknown. Here, we reveal a link between NAD⁺ capping and tissue- and hormone response-specific mRNA stability. In the absence of DXO1 function, transcripts displaying a high proportion of NAD⁺ capping are instead processed into RNA-dependent RNA polymerase 6-dependent small RNAs, resulting in their continued turnover likely to free the NAD⁺ molecules. Additionally, the NAD⁺-capped transcriptome is significantly remodeled in response to the essential plant hormone abscisic acid in a mechanism that is primarily independent of DXO1. Overall, our findings reveal a previously uncharacterized and essential role of NAD⁺ capping in dynamically regulating transcript stability during specific physiological responses.

RdDM pathway components differentially modulate Tobamovirus symptom development

The crop yield losses induced by phytoviruses are mainly associated with the symptoms of the disease. DNA modifications as methylation can modulate the information coded by the sequence, process named epigenetics. Viral infection can change the expression patterns of different genes linked to defenses and symptoms. This work represents the initial step to expose the role of epigenetic process, in the production of symptoms associated with plants-virus interactions. Small RNAs (sRNAs) are important molecules for gene regulation in plants and play an essential role in plant-pathogen interactions. Researchers have evaluated the relationship between viral infections as well as the endogenous accumulation of sRNAs and the transcriptional changes associated with the production of symptoms, but little is known about a possible direct role of epigenetics, mediated by 24-nt sRNAs, in the induction of these symptoms. Using different RNA directed DNA methylation (RdDM) pathway mutants and a triple demethylase mutant; here we demonstrate that the disruption of RdDM pathway during viral infection produce alterations in the plant transcriptome and in consequence changes in plant symptoms. This study represents the initial step in exposing that DNA methylation directed by endogenous sRNAs has an important role, uncoupled to defense, in the production of symptoms associated with plant-virus interactions.

Genome-Wide Identification and Coexpression Network Analysis of DNA Methylation Pathway Genes and Their Differentiated Functions in Ginkgo biloba L

DNA methylation plays a vital role in diverse biological processes. DNA methyltransferases (DNMTs) genes and RNA-directed DNA methylation (RdDM)-related genes are key genes responsible for establishing and maintaining genome DNA methylation in plants. In the present study, we systematically identified nine GbDNMTs in Ginkgo biloba, including the three common families of GbMET1a/1b, GbCMT2, and GbDRMa/b/2a/2b/2c, and a fourth family—GbDNMT3—which is absent in most angiosperms. We also identified twenty RdDM-related genes, including four GbDCLs, six GbAGOs, and ten GbRDRs. Expression analysis of the genes showed the different patterns of individual genes, and 15 of 29 genes displayed expression change under five types of abiotic stress. Gene coexpression analysis and weighted gene co-expression network analysis (WGCNA) using 126 public transcriptomic datasets revealed that these genes were clustered into two groups. In group I, genes covered members from all six families which were preferentially expressed in the ovulate strobile and fruit. A gene ontology (GO) enrichment analysis of WGCNA modules indicated that group I genes were most correlated with the biological process of cell proliferation. Group II only consisted of RdDM-related genes, including GbDRMs, GbAGOs, and GbRDRs, but no GbDCLs, and these genes were specifically expressed in the cambium, suggesting that they may function in a dicer-like (DCL)-independent RdDM pathway in specific tissues. The gene module related to group II was most enriched in signal transduction, cell communication, and the response to the stimulus. These results demonstrate that gene family members could be conserved or diverged across species, and multi-member families in the same pathway may cluster into different modules to function differentially. The study provides insight into the DNA methylation genes and their possible functions in G. biloba, laying a foundation for the further study of DNA methylation in gymnosperms.

Versatile in vitro assay to recognize Cas9-induced mutations

The discovery of CRISPR/Cas9 has revolutionized molecular biology, and its impact on plant biotechnology and plant breeding cannot be over-estimated. In many plant species, its application for mutagenesis is now a routine procedure--if suitable target sites, sufficient expression of the Cas9 protein, and functioning sgRNAs are combined. sgRNAs differ in their efficiency, depending on parameters that are only poorly understood. Several software tools and experience from growing databases are supporting the design of sgRNAs, but some seemingly perfect sgRNAs turn out to be inefficient or fail entirely, and most data bases stem from work with mammalian cells. Different in vitro assays testing sgRNAs in reconstituted Cas9 complexes are available and useful to reduce the risk of failure, especially in plants when CRISPR/Cas9 application requires modifications within the germ line and laborious transformation protocols. Low sgRNA efficiency and long generation times in plants can also contribute to the workload and costs of screening for the wanted genome edits. Here, we present a protocol in which a simple, initial in vitro test for suitable sgRNAs is modified to accelerate genotyping of Cas9-induced mutations. We demonstrate applicability of our protocol for mutagenesis and mutation screen for specific genes in Arabidopsis, but the principle should be universally suitable to provide a simple, low-cost, and rapid method to identify edited genes also in other plants and other organisms.

Genome-wide identification and characterization of DCL, AGO and RDR gene families in Saccharum spontaneum

RNA silencing is a conserved mechanism in eukaryotic organisms to regulate gene expression. Argonaute (AGO), Dicer-like (DCL) and RNA-dependent RNA polymerase (RDR) proteins are critical components of RNA silencing, but how these gene families’ functions in sugarcane were largely unknown. Most stress-resistance genes in modern sugarcane cultivars (Saccharum spp.) were originated from wild species of Saccharum, for example S. spontaneum. Here, we used genome-wide analysis and a phylogenetic approach to identify four DCL, 21 AGO and 11 RDR genes in the S. spontaneum genome (termed SsDCL, SsAGO and SsRDR, respectively). Several genes, particularly some of the SsAGOs, appeared to have undergone tandem or segmental duplications events. RNA-sequencing data revealed that four SsAGO genes (SsAGO18c, SsAGO18b, SsAGO10e and SsAGO6b) and three SsRDR genes (SsRDR2b, SsRDR2d and SsRDR3) tended to have preferential expression in stem tissue, while SsRDR5 was preferentially expressed in leaves. qRT-PCR analysis showed that SsAGO10c, SsDCL2 and SsRDR6b expressions were strongly upregulated, whereas that of SsAGO18b, SsRDR1a, SsRDR2b/2d and SsRDR5 was significantly depressed in S. spontaneum plants exposed to PEG-induced dehydration stress or infected with Xanthomonas albilineans, causal agent of leaf scald disease of sugarcane, suggesting that these genes play important roles in responses of S. spontaneum to biotic and abiotic stresses.

Genome-Wide Identification of RNA Silencing-Related Genes and Their Expressional Analysis in Response to Heat Stress in Barley (Hordeum vulgare L.)

Barley (Hordeum vulgare L.) is an economically important crop cultivated in temperate climates all over the world. Adverse environmental factors negatively affect its survival and productivity. RNA silencing is a conserved pathway involved in the regulation of growth, development and stress responses. The key components of RNA silencing are the Dicer-like proteins (DCLs), Argonautes (AGOs) and RNA-dependent RNA polymerases (RDRs). Despite its economic importance, there is no available comprehensive report on barley RNA silencing machinery and its regulation. In this study, we in silico identified five DCL (HvDCL), eleven AGO (HvAGO) and seven RDR (HvRDR) genes in the barley genome. Genomic localization, phylogenetic analysis, domain organization and functional/catalytic motif identification were also performed. To understand the regulation of RNA silencing, we experimentally analysed the transcriptional changes in response to moderate, persistent or gradient heat stress treatments: transcriptional accumulation of siRNA- but not miRNA-based silencing factor was consistently detected. These results suggest that RNA silencing is dynamically regulated and may be involved in the coordination of development and environmental adaptation in barley. In summary, our work provides information about barley RNA silencing components and will be a ground for the selection of candidate factors and in-depth functional/mechanistic analyses.

Synergy between the anthocyanin and RDR6/SGS3/DCL4 siRNA pathways expose hidden features of Arabidopsis carbon metabolism

Anthocyanin pigments furnish a powerful visual output of the stress and metabolic status of Arabidopsis thaliana plants. Essential for pigment accumulation is TRANSPARENT TESTA19 (TT19), a glutathione S-transferase proposed to bind and stabilize anthocyanins, participating in their vacuolar sequestration, a function conserved across the flowering plants. Here, we report the identification of genetic suppressors that result in anthocyanin accumulation in the absence of TT19. We show that mutations in RDR6, SGS3, or DCL4 suppress the anthocyanin defect of tt19 by pushing carbon towards flavonoid biosynthesis. This effect is not unique to tt19 and extends to at least one other anthocyanin pathway gene mutant. This synergy between mutations in components of the RDR6-SGS3-DCL4 siRNA system and the flavonoid pathway reveals genetic/epigenetic mechanisms regulating metabolic fluxes.

TRANSPARENT TESTA19 (TT19) encodes a glutathione S-transferase which functions in anthocyanin stabilization and vacuolar transport. Here, by tt19 suppressor screening, the authors show that RDR6/SGS3/DCL4 siRNA pathway constituents synergistically interact with components of the flavonoid pathway to control carbon metabolism.

DIRS retrotransposons amplify via linear, single-stranded cDNA intermediates

The Dictyostelium Intermediate Repeat Sequence 1 (DIRS-1) is the name-giving member of the DIRS order of tyrosine recombinase retrotransposons. In Dictyostelium discoideum, DIRS-1 is highly amplified and enriched in heterochromatic centromers of the D. discoideum genome. We show here that DIRS-1 it tightly controlled by the D. discoideum RNA interference machinery and is only mobilized in mutants lacking either the RNA dependent RNA polymerase RrpC or the Argonaute protein AgnA. DIRS retrotransposons contain an internal complementary region (ICR) that is thought to be required to reconstitute a full-length element from incomplete RNA transcripts. Using different versions of D. discoideum DIRS-1 equipped with retrotransposition marker genes, we show experimentally that the ICR is in fact essential to complete retrotransposition. We further show that DIRS-1 produces a mixture of single-stranded, mostly linear extrachromosomal cDNA intermediates. If this cDNA is isolated and transformed into D. discoideum cells, it can be used by DIRS-1 proteins to complete productive retrotransposition. This work provides the first experimental evidence to propose a general retrotransposition mechanism of the class of DIRS like tyrosine recombinase retrotransposons.

Ty-1, a universal resistance gene against geminiviruses that is compromised by co-replication of a betasatellite

Tomato yellow leaf curl virus (TYLCV), a begomovirus, causes large yield losses and breeding for resistance is an effective way to combat this viral disease. The resistance gene Ty-1 codes for an RNA-dependent RNA polymerase and has recently been shown to enhance transcriptional gene silencing of TYLCV. Whereas Ty-1 was earlier shown to also confer resistance to a bipartite begomovirus, here it is shown that Ty-1 is probably generic to all geminiviruses. A tomato Ty-1 introgression line, but also stable transformants of susceptible tomato cv. Moneymaker and Nicotiana benthamiana (N. benthamiana) expressing the Ty-1 gene, exhibited resistance to begomoviruses as well as to the distinct, leafhopper-transmitted beet curly top virus, a curtovirus. Stable Ty-1 transformants of N. benthamiana and tomato showed fewer symptoms and reduced viral titres on infection compared to wild-type plants. TYLCV infections in wild-type N. benthamiana plants in the additional presence of a betasatellite led to increased symptom severity and a consistent, slightly lowered virus titre relative to the high averaged levels seen in the absence of the betasatellite. On the contrary, in Ty-1 transformed N. benthamiana viral titres increased in the presence of the betasatellite. The same was observed when these Ty-1-encoding plants were challenged with TYLCV and a potato virus X construct expressing the RNA interference suppressor protein βC1 encoded by the betasatellite. The resistance spectrum of Ty-1 and the durability of the resistance are discussed in light of antiviral RNA interference and viral counter defence strategies.

The Whys and Wherefores of Transitivity in Plants

Transitivity in plants is a mechanism that produces secondary small interfering RNAs (siRNAs) from a transcript targeted by primary small RNAs (sRNAs). It expands the silencing signal to additional sequences of the transcript. The process requires RNA-dependent RNA polymerases (RDRs), which convert single-stranded RNA targets into a double-stranded (ds) RNA, the precursor of siRNAs and is critical for effective and amplified responses to virus infection. It is also important for the production of endogenous secondary siRNAs, such as phased siRNAs (phasiRNAs), which regulate several genes involved in development and adaptation. Transitivity on endogenous transcripts is very specific, utilizing special primary sRNAs, such as miRNAs with unique features, and particular ARGONAUTEs. In contrast, transitivity on transgene and virus (exogenous) transcripts is more generic. This dichotomy of responses implies the existence of a mechanism that differentiates self from non-self targets. In this work, we examine the possible mechanistic process behind the dichotomy and the intriguing counter-intuitive directionality of transitive sequence-spread in plants.

Current overview on viroid-host interactions

Viroids are one of the most enigmatic highly structured, circular, single-stranded RNA phytopathogens. Although they are not known to code for any peptide, viroids induce visible symptoms in susceptible host plants that resemble those associated with many plant viruses. It is known that viroids induce disease symptoms by direct interaction with host factors; however, the precise mechanism by which this occurs remains poorly understood. Studies on the host's responses to viroid infection, host susceptibility and nonhost resistance have been underway for several years, but much remains to be done in order to fully understand the complex nature of viroid-host interactions. Recent progress using molecular biology techniques combined with computational algorithms, in particular evidence of the role of viroid-derived small RNAs in the RNA silencing pathways of a disease network, has widened the knowledge of viroid pathogenicity. The complexity of viroid-host interactions has been revealed in the past decades to include, but not be limited to, the involvement of host factors, viroid structural complexity, and viroid-induced ribosomal stress, which is further boosted by the discovery of long noncoding RNAs (lncRNAs). In this review, the current understanding of the viroid-host interaction has been summarized with the goal of simplifying the complexity of viroid biology for future research. This article is categorized under: RNA in Disease and Development > RNA in Disease.

CRISPR/Cas9-mediated knockout of the RDR6 gene in Nicotiana benthamiana for efficient transient expression of recombinant proteins

MAIN CONCLUSION

RDR6 gene knockout Nicotiana benthamiana plant was successfully produced using CRISPR/Cas9 technology. The production of recombinant proteins in plants has many advantages, such as safety and reduced costs. However, there are several problems with this technology, especially low levels of protein production. The dysfunction of the RNA silencing mechanism in plant cells would be effective to improve recombinant protein production because the RNA silencing mechanism efficiently degrades transgene-derived mRNAs. Therefore, to overcome this problem, clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 technology was used to develop RNA silencing-related gene knockout transgenic Nicotiana benthamiana. We successfully produced RNA-dependent RNA polymerase 6 (RDR6), one of the most important components of the RNA silencing mechanism-knockout N. benthamiana (ΔRDR6 plants). The ΔRDR6 plants had abnormal flowers and were sterile, as with the Arabidopsis RDR6 mutants. However, a transient gene expression assay showed that the ΔRDR6 plants accumulated larger amounts of green fluorescent protein (GFP) and GFP mRNA than the wild-type (WT) plants. Small RNA sequencing analysis revealed that levels of small interfering RNA against the GFP gene were greatly reduced in the ΔRDR6 plants, as compared to that of the WT plants. These findings demonstrate that the ΔRDR6 plants can express larger amounts of recombinant proteins than WT plants and, therefore, would be useful for recombinant protein production and understanding the contributions of RDR6 to genetic and physiological events in plants.

Identification of transcription factor binding sites on promoter of RNA dependent RNA polymerases (RDRs) and interacting partners of RDR proteins through in silico analysis

RNA silencing phenomenon in plants provides resistance to various pathogens and also, it maintains genome integrity. The process of RNA silencing is regulated by diverse proteins, among which RNA dependent RNA polymerases (RDRs) are very crucial for the amplification of small RNAs (sRNAs). Out of various RDR proteins present in plants, role of RDR1, RDR2 and RDR6 for providing resistance against various biotic stresses have been well documented. In contrast, very few information is available regarding the role of RDR3, RDR4 and RDR5 proteins in plant biology and stress response. Furthermore, the regulation of RDRs is not yet known. Here, we have carried out in silico studies for identification of the transcription factor (TF) binding sites on the promoter of RDR1-6 genes of various plant species. Among the TFs predicted to bind on the promoter of RDRs, MYB44, AS1/AS2, WRKY1 are the major one. Furthermore, putative interacting protein partners of RDRs proteins of tomato and rice were also predicted by STRING database which suggests that DCL (Dicer-like) proteins are strong candidate proteins as the interacting partners of RDRs. The knowledge of regulation of RDRs and its interacting protein partners might help in developing resistant plants to biotic stresses.

Suppression of RNA-Dependent RNA Polymerase6 Favors the Accumulation of Potato Spindle Tuber Viroid in Nicotiana Benthamiana

To date, two plant genes encoding RNA-dependent RNA polymerases (RdRs) that play major roles in the defense against RNA viruses have been identified: (i) RdR1, which is responsible for the viral small RNAs (vsRNAs) found in virus-infected plants, and, (ii) RdR6, which acts as a surrogate in the absence of RdR1. In this study, the role of RdR6 in the defense against viroid infection was examined by knock-down of RdR6 followed by potato spindle tuber viroid (PSTVd) infection. The suppression of RdR6 expression increased the plant's growth, as was illustrated by the plant's increased height. PSTVd infection of RdR6 compromised plants resulted in an approximately three-fold increase in the accumulation of viroid RNA as compared to that seen in control plants. Additionally, RNA gel blot assay revealed an increase in the number of viroids derived small RNAs in RdR6 suppressed plants as compared to control plants. These data provide a direct correlation between RdR6 and viroid accumulation and indicate the role of RDR6 in the plant's susceptibility to viroid infection.