Abstract
RNA interference (RNAi) is a biotechnological tool used for gene silencing in plants, with both endogenous and exogenous applications. Endogenous approaches, such as host-induced gene silencing (HIGS), involve genetically modified (GM) plants, while exogenous methods include spray-induced gene silencing (SIGS). The RNAi mechanism hinges on the introduction of double-stranded RNA (dsRNA), which is processed into short interfering RNAs (siRNAs) that degrade specific messenger RNAs (mRNAs). However, unintended effects on non-target organisms and GM plants are a concern due to sequence homologies or siRNA-induced epigenetic changes. Regulatory bodies such as the EPA and EFSA emphasize the need for comprehensive risk assessments. Detecting unintended effects is complex, often relying on bioinformatic tools and untargeted analyses like transcriptomics and metabolomics, though these methods require extensive genomic data. This review aims to classify mechanisms of RNAi effects induced by short interfering RNA from different sources in plants and to identify technologies that can be used to detect these effects. In addition, practical case studies are summarized and discussed in which previously unintended RNAi effects in genetically modified plants have been investigated. Current literature is limited but suggests RNAi is relatively specific, with few unintended effects observed in GM crops. However, further studies are needed to fully understand and mitigate potential risks, particularly those related to transcriptional gene silencing (TGS) mechanisms, which are less predictable than post-transcriptional gene silencing (PTGS). Particularly the application of untargeted approaches such as small RNA sequencing and transcriptomics is recommended for thorough and comprehensive risk assessments.
Introduction
RNA interference (RNAi) represents a cutting-edge approach in biotechnology for gene expression silencing, applied e.g. in plant protection, leveraging molecular principles to control gene expression. This innovative strategy encompasses both endogenous and exogenous applications, each with distinct methodologies and implications. Genetically modified (GM) plants harness RNAi to target plant endogenous transcripts e.g. to regulate the gibberellin pathway (maize event MON 94804) or to alter the fatty acid profile (soy event MON 87705). Endogenous applications also involve GM plants in a process known as host-induced gene silencing (HIGS) (Nowara et al., 2010) for pesticidal applications (e.g. maize event MON 87411 containing dsSnf7 against Diabrotica), (see https://euginius.eu). Conversely, exogenous applications, such as spray-induced gene silencing (SIGS), or root soaking of RNAi involve the direct application of RNA molecules to plants (Liu et al., 2020; Werner et al., 2020).
The core mechanism of RNAi in biotechnology application, such as plant protection, lies in its ability to selectively reduce the expression of specific genes within the target organism (Koeppe et al., 2023). In the majority of cases, this is achieved through the introduction of double-stranded RNA (dsRNA), which is subsequently processed by the RNase III Dicer or related enzymes to short interfering (si)RNA, whose base pairing with the complementary sequence of the target messenger (m)RNA leads to its degradation (Guo et al., 2016; Hung and Slotkin, 2021). While this sequence-based mechanism is advantageous for targeting pests and pathogens, there is a potential for unintended effects on non-target organisms (NTOs) and the GM plant itself (Christiaens et al., 2018). These effects may arise due to sequence homologies between the dsRNA and non-target mRNAs or through mechanisms such as siRNA-induced epigenetic changes and disruption of the organism’s endogenous RNAi pathways (Kloc et al., 2008; Zaratiegui and Martienssen, 2012; Swevers et al., 2013).
Recognizing the novel challenges posed by RNAi-based plant protection, regulatory bodies such as the US Environmental Protection Agency (EPA) and the European Food Safety Authority (EFSA) have acknowledged the need for comprehensive risk assessments (Christiaens et al., 2018; Papadopoulou et al., 2020; Christiaens et al., 2022). The Chemicals Committee and the Working Party on Chemicals, Pesticides and Biotechnology of the Organisation for Economic Co-operation and Development (OECD) have compiled considerations to integrate the latest scientific understanding into the environmental risk assessment of RNAi applications (Organisation for Economic Co-operation and Development (OECD) 2020).
One significant concern is the potential for unintended effects on GM plants themselves. Detecting these effects is complex due to several factors. Current prediction methods primarily rely on bioinformatic searches for complementary sequences to the siRNA within the GM plant’s transcriptome (Good et al., 2016; Lück et al., 2019; Farooq et al., 2021). However, these analyses are often hampered by the lack of a complete and accurate reference genome for the GM plant. When available, reference genomes of closely related cultivars may be used, but these can lead to inaccuracies due to sequence polymorphisms, resulting in false positives or negatives in off-target effect predictions.
In this review, we summarize the mechanisms by which RNAi applications could induce unintended effects in plants and evaluate the technologies and approaches available to detect these effects. By assessing the relevance of RNAi-mediated cellular mechanisms to GM plants based on existing literature, we provide a comprehensive overview and aim to rank these mechanisms according to their significance. This detailed examination will contribute to a better understanding of RNAi applications and the development of more accurate risk assessment methodologies.
Mechanisms of RNAi-induced effects in plants
The principle of RNAi in plant protection relies on reducing or silencing the expression of specific essential genes in the target organism or the GM plant itself. These target genes typically belong to vital metabolic or developmental pathways, leading to a loss-of-function phenotype (Werner et al., 2020; Hernández-Soto and Chacón-Cerdas, 2021). RNAi-based pest control strategies primarily utilize two types of RNA precursors: short hairpin RNAs (shRNA), which consist of two complementary strands forming a stem-loop structure, and complementary dsRNA. The enzyme Dicer, found in nearly all eukaryotes with various isotypes (Zapletal et al., 2023), processes these precursor molecules into short, mostly 21-24 nucleotide (nt) RNA duplexes in the cytoplasm (Figure 1). In plants, Dicer-like (DCL) proteins play an important role in processing dsRNA into siRNAs of different length (Henderson et al., 2006; Mukherjee et al., 2013). The RNA duplexes include a guide strand and a passenger strand [reviewed in (Kim et al., 2009; Borges and Martienssen, 2015)]. While the passenger strand is degraded during further processing, the guide strand, which is complementary to the target gene sequence, is crucial for the silencing of the gene. In the following, we will first focus on the biogenesis of small RNAs in plants and then discuss the mechanisms of RNAi-based silencing before we discuss the implications of these mechanisms for possible off-target effects in GM plants.
Figure 1
Biogenesis of small RNAs in plants
To investigate the effects of genetic modifications on the RNAi pathway in GM plants, it is essential to consider the natural mechanisms by which RNAi can affect gene expression in plants. To this end, the cellular pathways by which siRNA molecules can be produced in plants are first described here (Figure 1; Table 1). Precursors of siRNA are almost without exception double-stranded RNA molecules, which are either synthesized by endogenous RNA polymerases (RNA Pol) or introduced exogenously (Vazquez et al., 2004; Allen et al., 2005; Axtell et al., 2006). Endogenous precursors include natural antisense transcripts (NAT) synthesized by RNA Pol II, which base-pair with the sense mRNA of the coding gene and thus form the double-stranded substrate for corresponding RNases (Figure 1A) (Borsani et al., 2005; Jen et al., 2005; Zhang et al., 2012). RNA Pol II also synthesizes shRNAs encoded in the genome, which can then be processed by Dicer into siRNA (Figure 1B) (Wesley et al., 2001; Helliwell and Waterhouse, 2003; Senthil-Kumar and Mysore, 2011) or long non-coding (lnc) RNAs (Kim and Sung, 2012; Liu et al., 2012; Wu et al., 2012), whose secondary structures can have hairpins and can thus also be converted into siRNA by corresponding RNases (Figure 1C). Endogenously encoded micro (mi)RNAs are synthesized by Dicer or DCL1 in plants (Kurihara and Watanabe, 2004) from shRNAs, the miRNA precursors (Figure 1D), and either directly regulate the expression of target genes (by miRNA) or base-pair with the precursors of so-called trans-acting (ta)siRNAs, which are then generated by DCL from a double-stranded template (Figure 1E). Exogenously introduced precursors of siRNA are molecules introduced into a cell from an external source. A natural example are viral RNAs, which are either immediately present after infection and replication (RNA viruses, in plant viruses often single-stranded (ss) RNA genome) or are generated by transcription of the viral genome (DNA viruses) and are then templates for DCLs, which produce siRNA from them (Figure 1F) (Ruiz et al., 1998; Lu et al., 2003; Burch-Smith et al., 2004). RNA Pol IV or V can also be involved in siRNA synthesis, for example in the case of the synthesis of precursors of heterochromatic (hc)siRNAs, which are then converted to siRNA by DCL3 (Figure 1G) (Law and Jacobsen, 2010; Zhang and Zhu, 2011; Matzke and Mosher, 2014). In plants, the proteins DCL2, DCL3 and DCL4 generate siRNAs of different lengths mostly with 22 nt, 24 nt and 21 nt, respectively, which in turn trigger different mechanisms of silencing (Henderson et al., 2006; Mukherjee et al., 2013). Of note, DCL2-derived 22 nt siRNAs in plants are involved in a transitive and systemic spread of siRNA especially for antiviral defense, called secondary RNAi (Bouché et al., 2006; Chen et al., 2010; Garcia-Ruiz et al., 2010; Qin et al., 2017). This spread of RNAi involves the amplification and expansion of silencing signals that are mediated by RNA-dependent RNA polymerases (RdRp) (Sanan-Mishra et al., 2021). In this process, siRNAs act on longer RNAs (such as mRNA) as primers for RdRp, whereby a new, long dsRNA is synthesized, which is then eventually processed again by the RNAi machinery into siRNA triggering secondary RNAi.
Table 1
| Category | Mechanism | Description | Source of siRNA | References |
|---|---|---|---|---|
| Post-Transcriptional Gene Silencing (PTGS) | siRNA-Mediated mRNA Degradation | siRNAs guide RISC to complementary mRNA, leading to its cleavage and degradation. | Trans-Acting siRNAs (ta-siRNAs) | (Vazquez et al., 2004; Allen et al., 2005; Axtell et al., 2006) |
| Natural Antisense Transcripts (NATs) | (Borsani et al., 2005; Jen et al., 2005; Zhang et al., 2012) | |||
| Exogenous short hairpin RNA (shRNA) | (Wesley et al., 2001; Helliwell and Waterhouse, 2003; Senthil-Kumar and Mysore, 2011) | |||
| Long Non-Coding RNAs (lncRNAs) | (Kim and Sung, 2012; Liu et al., 2012; Wu et al., 2012) | |||
| Virus-derived siRNAs | (Ruiz et al., 1998; Lu et al., 2003; Burch-Smith et al., 2004) | |||
| Heterochromatic siRNAs (hc-siRNAs) | (Law and Jacobsen, 2010; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| siRNA-Mediated Translation inhibition | siRNAs guide RISC to complementary mRNA, leading to translation inhibition | Trans-Acting siRNAs (ta-siRNAs) | (Vazquez et al., 2004; Allen et al., 2005; Axtell et al., 2006) | |
| Natural Antisense Transcripts (NATs) | (Borsani et al., 2005; Jen et al., 2005; Zhang et al., 2012) | |||
| Exogenous short hairpin RNA (shRNA) | (Wesley et al., 2001; Helliwell and Waterhouse, 2003; Senthil-Kumar and Mysore, 2011) | |||
| Long Non-Coding RNAs (lncRNAs) | (Heo and Sung, 2011; Kim and Sung, 2012; Liu et al., 2012) | |||
| Virus-Derived siRNAs (VIGS) | (Ruiz et al., 1998; Lu et al., 2003; Burch-Smith et al., 2004) | |||
| miRNA-Mediated mRNA Degradation | miRNAs guide RISC to complementary or partially complementary mRNAs, resulting in cleavage or repression. | Endogenously expressed miRNAs | (Jones-Rhoades et al., 2006; Mallory and Vaucheret, 2006; Voinnet, 2009) | |
| Transcriptional Gene Silencing (TGS) | DNA Methylation | siRNAs guide DNA methylation machinery to specific genomic regions, repressing transcription. | Long Non-Coding RNAs (lncRNAs) | (Wierzbicki et al., 2008; Zhang and Zhu, 2011; Matzke and Mosher, 2014) |
| Heterochomatic siRNAs (hc-siRNAs) | (Law and Jacobsen, 2010; Zhang and Zhu, 2011; Movahedi et al., 2015) | |||
| Pol IV/Pol V-derived siRNAs | (Law and Jacobsen, 2010; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Endogenously expressed miRNAs | (Wu et al., 2010; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Trans-Acting siRNAs (ta-siRNAs) | (Allen et al., 2005; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Natural Antisense Transcripts (NATs) | (Borsani et al., 2005; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Exogenous short hairpin RNA (shRNA) | (Waterhouse and Helliwell, 2003; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Histone Modification | siRNAs direct histone-modifying enzymes to specific loci, causing chromatin condensation and gene silencing. | Long Non-Coding RNAs (lncRNAs) | (He et al., 2011; Zhang and Zhu, 2011; Liu et al., 2012) | |
| Heterochomatic siRNAs (hc-siRNAs) | (Law and Jacobsen, 2010; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Pol IV/Pol V-derived siRNAs | (Law and Jacobsen, 2010; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Trans-Acting siRNAs (ta-siRNAs) | (Axtell et al., 2006; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Natural Antisense Transcripts (NATs) | (Borsani et al., 2005; Zhang and Zhu, 2011; Matzke and Mosher, 2014) | |||
| Exogenous short hairpin RNA (shRNA) | (Helliwell and Waterhouse, 2003; Law and Jacobsen, 2010; Senthil-Kumar and Mysore, 2011) | |||
| Virus-Derived siRNAs (VIGS) | (Baulcombe, 2004; Blevins et al., 2006; Matzke and Mosher, 2014) |
Mechanisms by which RNA interference induces gene expression changes in plants, categorized by their general mode of action, including the mechanism, the source of siRNA and corresponding references.
Mechanisms of RNAi-based silencing
Silencing mechanisms can occur in the GM plant harboring the RNAi construct, at the transcriptional level in the cell nucleus or the translational/post-transcriptional level in the cytoplasm (Figure 2). In the nucleus, siRNA can pair with the nascent mRNA of the target gene, recruiting factors to the transcription machinery that inhibit the transcription elongation by RNA polymerase (Figure 2A, left) (Guang et al., 2010). Similarly, siRNA can recruit enzymes that induce epigenetic silencing of the target gene through DNA methylation (Law and Jacobsen, 2010; Wu et al., 2010; Zhang and Zhu, 2011; Wu et al., 2012; Matzke and Mosher, 2014; Movahedi et al., 2015) or histone modification (Baulcombe, 2004; He et al., 2011; Liu et al., 2012) (Figure 2A, right) (Verdel et al., 2009). In plants, epigenetic silencing via DNA methylation is triggered by DCL3-generated ~24 nt siRNA involving a RISC complex containing the protein Argonaute (Ago)4 (Zilberman et al., 2003; Henderson et al., 2006; Qi et al., 2006; Zheng et al., 2007; Wierzbicki et al., 2009; Havecker et al., 2010; Olmedo-Monfil et al., 2010; Sarkies and Miska, 2014; Lewsey et al., 2016). The most well-studied RNAi silencing mechanism involves the degradation of the target gene’s mRNA (Figure 2B, left). In this process, the protein Ago recruits siRNA to the complementary mRNA sequence to form the RNA-induced silencing complex (RISC) (Wu et al., 2012). In plants, this is triggered by DCL4-generated ~21 nt siRNAs involving a RISC complex containing Ago1 (Xie et al., 2005; Qu et al., 2008; Chen et al., 2010; Wang et al., 2011). If there is perfect complementarity between siRNA and the target gene, the mRNA is degraded, leading to down-regulation of the target protein’s production (Valencia-Sanchez et al., 2006). With incomplete base pairing between siRNA and target mRNA, the RNA is not degraded; instead, ribosome-mediated translation is inhibited, resulting in reduced expression of the target gene (Figure 2B, right) (Brodersen et al., 2008).
Figure 2
The current literature suggests that siRNA molecules produced via different biogenesis pathways can differ in terms of their length, triggering different types of mechanisms of gene expression regulation described. While DCL4-generated 21 nt siRNA predominantly triggers PTGS via mRNA degradation, DCL3-generated 24 nt siRNA triggers TGS via epigenetic silencing and DCL2-generated 22 nt siRNA induces secondary siRNA. However, all DCL may act on long dsRNA molecules introduced into the plant. Therefore, both TGS and PTGS need to be considered when analyzing RNAi-induced effects in GM plants (Table 1).
Implications of RNAi silencing mechanisms for possible off-target effects in GM plants
With regard to the knowledge about mechanisms by which RNAi can potentially induce unintended effects in the GM plant, the existing literature shows clear bias towards PTGS. For example, a PubMed search with the search term “RNAi AND PTGS NOT TGS” in title and abstract returned 118 hits, whereas the search term “RNAi AND TGS NOT PTGS” only returned 25 hits (as of 18.10.2024). Hence, most published studies are concerned with the investigation of effects resulting from the inhibition of translation or degradation of mRNA (possibly resulting from incomplete complementarity). Relatively fewer studies deal with TGS, possibly because here effects, for example via epigenetic silencing, could also arise upstream or downstream of the gene with sequence complementarity and these cannot be clearly determined on the basis of the pure small RNA sequence by analyzing complementary sequences in the genome.
Techniques for the assessment of RNAi-induced effects in plants
RNAi can induce different types of off-target effects in the plant, which can be identified and studied using different techniques. Here we provide a brief overview of the different techniques that can be used to study the changes induced by RNAi and RNAi off-target effects. The methods employed to study RNAi effects can be divided into two main approaches: targeted and untargeted analysis (Table 2). Targeted screening of RNAi effects focuses on analyzing the intended silencing effects on specific target genes and includes, for example, validation of gene knockdown, functional assays, validation of phenotypic effects, assessment of specificity and long-term effects (e.g. stability of gene silencing). Targeted screening can also be used to analyze effects on predicted off-target genes. The corresponding techniques include molecular techniques such as RT qPCR (Chi et al., 2008; Sun and Rossi, 2009; Holmes et al., 2010; Varkonyi-Gasic and Hellens, 2011; Augustine et al., 2013; Kitzmann et al., 2013; Liu et al., 2014; Czarnecki et al., 2016; Keykha et al., 2016; Manske et al., 2017; Betti et al., 2021; Sarkar and Roy-Barman, 2021; Xu et al., 2021; Zhou et al., 2021; López-Márquez et al., 2023; Kyslík et al., 2024), northern blotting (Chi et al., 2008; Fukuhara et al., 2011; Augustine et al., 2013; Manske et al., 2017; Sarkar and Roy-Barman, 2021), western blotting (Kumar et al., 2003; Sahin et al., 2007; Sun and Rossi, 2009; Holmes et al., 2010; Liang et al., 2013; Han, 2018; Vidarsdottir et al., 2019; Kyslík et al., 2024), genetic techniques such as reporter gene assays (Kumar et al., 2003; Smart et al., 2005; Rinaldi et al., 2008; Sun and Rossi, 2009; Manske et al., 2017; López-Márquez et al., 2023) or genetic mutations (Chan et al., 2006; Czarnecki et al., 2016; Krzyszton and Kufel, 2022), phenotypic assays (Chi et al., 2008; Liu et al., 2014; Xu et al., 2021; Zhou et al., 2021; Tao et al., 2023), enzyme activity assays (Chi et al., 2008; Betti et al., 2021; Sarkar and Roy-Barman, 2021) or advanced techniques such as genome editing using CRISPR/Cas9 (Moore, 2015; Kanchiswamy et al., 2016; Peretz et al., 2018; Kleter, 2020; Mujtaba et al., 2021; Bock et al., 2022).
Table 2
| Class | Category | Technique | Description | References |
|---|---|---|---|---|
| Targeted analysis | Molecular techniques | Quantitative RT-PCR (RT-qPCR) | Quantification of specific RNA molecules to assess reduction in target or off-target genes. advantages: high specificity; disadvantages: limited to a minimum RNA size, sequence information required, single gene analysis | (Chi et al., 2008; Sun and Rossi, 2009; Holmes et al., 2010; Varkonyi-Gasic and Hellens, 2011; Augustine et al., 2013; Kitzmann et al., 2013; Liu et al., 2014; Czarnecki et al., 2016; Keykha et al., 2016; Manske et al., 2017; Betti et al., 2021; Sarkar and Roy-Barman, 2021; Xu et al., 2021; Zhou et al., 2021; López-Márquez et al., 2023; Kyslík et al., 2024) |
| Northern blotting | Quantification of specific RNA molecules to assess reduction in target or off-target genes as well as specific detection of siRNAs. advantages: high specificity, detection of fragments possible, detection of short RNA molecules; disadvantages: sequence information required, time-intensiveness, single gene analysis | (Chi et al., 2008; Fukuhara et al., 2011; Augustine et al., 2013; Manske et al., 2017; Sarkar and Roy-Barman, 2021) | ||
| Western blotting | Quantification of specific proteins to assess reduction in target or off-target genes at the protein level. advantages: quantification of gene products, integrated assessment of mRNA degradation and translation inhibition; disadvantages: requires specific antibodies, single protein analysis | (Kumar et al., 2003; Sahin et al., 2007; Sun and Rossi, 2009; Holmes et al., 2010; Liang et al., 2013; Han, 2018; Vidarsdottir et al., 2019; Kyslík et al., 2024) | ||
| Genetic techniques | Reporter gene assays | Assessment of the effect of RNAi on target gene expression in a GM reporter system. advantages: quick assessment, clear read out; disadvantages: artificial system, single gene analysis | (Kumar et al., 2003; Smart et al., 2005; Rinaldi et al., 2008; Sun and Rossi, 2009; Manske et al., 2017; López-Márquez et al., 2023) | |
| Genetic mutants | Comparison of RNAi effects with genetic mutants to validate phenotypic effects of gene knockdown. advantages: coverage of all phenotypic effects; disadvantages: mutant required, no discrimination between target and off-target effects | (Chan et al., 2006; Czarnecki et al., 2016; Krzyszton and Kufel, 2022) | ||
| Phenotypic techniques | Phenotypic assays | Measurement of physiological parameters such as photosynthetic efficiency and hormone levels comparing RNAi-based GMP with unmodified comparator. advantages: identification of physiological parameters; disadvantages: no discrimination between target and off-target effects | (Chi et al., 2008; Liu et al., 2014; Xu et al., 2021; Zhou et al., 2021; Tao et al., 2023) | |
| Biochemical techniques | Enzyme activity assays | Measurement of activity of enzymes encoded by target genes to confirm functional consequences. advantages: assessment of functional consequences; disadvantages: restriction to the target gene, no assessment of off-target effects | (Chi et al., 2008; Betti et al., 2021; Sarkar and Roy-Barman, 2021) | |
| Advanced techniques | CRISPR/Cas9 | CRISPR/Cas9 gene editing for validation of RNAi effects by knocking out target genes; Editing of potential off-target genes to assess phenotypic outcomes compared to RNAi treatments. advantages: specific analysis of phenotypic changes induced by target knockout; disadvantages: knowledge about off-target required, time consuming for a number of off-targets | (Moore, 2015; Kanchiswamy et al., 2016; Peretz et al., 2018; Kleter, 2020; Mujtaba et al., 2021; Bock et al., 2022) | |
| Untargeted analysis | Molecular techniques | RNA sequencing (RNA-Seq) | Next generation sequencing of RNA and differential gene expression analysis to assess target and off-target genes. advantages: global detection of gene expression changes; disadvantages: time consuming, ideally availability of (high quality) reference genome required (workaround: de novo transcriptome assembly) | (Chan et al., 2006; Surget-Groba and Montoya-Burgos, 2010; Haque and Nishiguchi, 2011; Narzisi and Mishra, 2011; Jiao et al., 2021; Xu et al., 2021; Gaffo et al., 2022; Krzyszton and Kufel, 2022; Nguyen et al., 2022; Tyagi et al., 2022; Dong et al., 2023; López-Márquez et al., 2023; Budnick et al., 2024; Cazares et al., 2024) |
| Small RNA sequencing (small RNA-Seq) | Next generation sequencing of small RNAs for siRNA and miRNA quantification e.g. to inform bioinformatic off-target predictions. advantages: global analysis of small RNAs; disadvantages: time consuming | (Huang et al., 2024) | ||
| Phenotypic techniques | Visual observation | Observation of plant phenotypes such as growth rate, seed weight, leaf shape, and flower development. advantages: global assessment of phenotypic effects; disadvantages: no discrimination between target and off-target effects | (Augustine et al., 2013; Manske et al., 2017; López-Márquez et al., 2023) | |
| Microscopy | Microscopic observation cellular and subcellular changes. advantages: assessment of effects at the cellular and subcellular level; disadvantages: no discrimination between target and off-target effects | (Chi et al., 2008; Kitzmann et al., 2013; Betti et al., 2021; Sarkar and Roy-Barman, 2021; Xu et al., 2021; Zhou et al., 2021; Kyslík et al., 2024) | ||
| Biochemical techniques | Metabolite profiling | Gas chromatography (GC) coupled mass spectrometry (MS) analysis of changes in metabolite levels to provide insights into affected metabolic pathways. advantages: global assessment of metabolites; disadvantages: time-consuming, no discrimination between target and off-target effects | (Chen et al., 2012; Huang et al., 2022; Baysoy et al., 2023; Bressan et al., 2023; Naik et al., 2023; Huang et al., 2024) | |
| Advanced techniques | Proteomics | Liquid chromatography coupled mass spectrometry (LC-MS) to identify changes in protein abundance and post-translational modifications. advantages: global detection of gene expression changes at the protein level; disadvantages: lower sensitivity, time consuming | (Chi et al., 2008; Lacourse et al., 2008; Asano and Nishiuchi, 2011; Chen et al., 2012; Naik et al., 2023) | |
| Chromatin immunoprecipitation sequencing (ChIP-Seq) | Assessment of changes in DNA methylation or histone modifications and transcription factor binding as a result of RNAi. advantages: global assessment of epigenetic TGS; disadvantages: time consuming | (Warnatz et al., 2011; Muhammad et al., 2020; Navarro-Mendoza et al., 2023) |
Techniques for detecting RNAi off target effects in plants categorized by class and field. Short descriptions of each technique as well as the corresponding references are given.
Untargeted screening of RNAi effects involves comprehensive analyses mainly aimed at identifying unintended consequences and potential unpredicted off-target effects of RNAi treatments. These techniques include analyzing changes in transcriptomic profiles (Chan et al., 2006; Surget-Groba and Montoya-Burgos, 2010; Haque and Nishiguchi, 2011; Narzisi and Mishra, 2011; Jiao et al., 2021; Xu et al., 2021; Gaffo et al., 2022; Krzyszton and Kufel, 2022; Nguyen et al., 2022; Tyagi et al., 2022; Dong et al., 2023; López-Márquez et al., 2023; Budnick et al., 2024; Cazares et al., 2024), changes in protein expression (Chi et al., 2008; Lacourse et al., 2008; Asano and Nishiuchi, 2011; Chen et al., 2012; Naik et al., 2023) and modifications, metabolites (Chen et al., 2012; Huang et al., 2022; Baysoy et al., 2023; Bressan et al., 2023; Naik et al., 2023; Huang et al., 2024) and epigenetic changes (Warnatz et al., 2011; Muhammad et al., 2020; Navarro-Mendoza et al., 2023) to understand the downstream effects of RNAi on cellular processes. In addition, the distribution and potential off-target interactions of RNAi (small RNAs) with unintended mRNA targets can be determined. Furthermore, there are also bioinformatic tools that utilize computational algorithms to predict potential off-target sites based on sequence complementarity and thermodynamic stability (Good et al., 2016; Lück et al., 2019). However, such bioinformatic prediction tools require extensive knowledge, for example of the plant’s genome or its RNAi machinery, in order to apply them effectively.
When studying off-target effects of RNAi, both targeted and untargeted analyses offer unique advantages and disadvantages. Targeted analysis as focuses on predefined genes or pathways, provide specific and efficient validation of RNAi-induced gene silencing. It ensures detailed understanding of intended effects but has a limited scope, potentially missing broader biological impacts and introducing bias by overlooking unexpected interactions. These techniques require fewer technical resources and their costs are reduced, making targeted analysis well suited as validation techniques. In contrast, untargeted analysis provides a comprehensive, genome/proteome/transcriptome-wide assessment, enabling the discovery of both known and unknown off-target interactions. However, this approach depends on high-quality, well-annotated genomes for precise mapping of RNAi-induced changes and understanding the broader implications of gene silencing in plants. While this unbiased method generates extensive datasets that provide deeper insights into RNAi effects, it is resource-intensive and complex, demanding substantial time, computational power, and expertise for analysis and interpretation. Additionally, the large datasets can introduce noise, probably requiring further validation to identify meaningful effects. Despite potential challenges, combining both approaches can offer a balanced perspective, profiting the specificity of targeted analysis and the breadth of untargeted analysis to achieve thorough insights into RNAi effects.
Discussion
Relevance of unintended effects of RNAi for risk assessment
Unintended effects of RNAi applications in GM plants themselves are a critical focus in the safety assessment of food and feed. Consequently, the Food and Agriculture Organization of the United Nations, for example, has issued guidelines for conducting food safety assessments of food derived from recombinant DNA plants (Food and Agriculture Organization of the United Nations, 2003). Also the OECD publishes science-based consensus documents offering information for the regulatory assessments of specific food and feed products, including those derived from transgenic organisms (Organisation for Economic Co-operation and Development, 2021). These documents gather data on the product’s nutrients, anti-nutrients and toxicants, its use as food or feed, and other factors relevant to food and feed safety. Here and in various review articles on the topic of risk assessment of RNAi-based GM crops, primarily untargeted methods for analyzing gene products and their metabolites, such as proteomics and metabolomics, are proposed to investigate RNAi-induced effects in the GM crop itself (Senthil-Kumar and Mysore, 2011; Kleter, 2020; Papadopoulou et al., 2020; Chaudhary et al., 2024).
The mechanisms by which the RNAi pathway can trigger specific gene expression changes in plants include both transcriptional and post-transcriptional regulation. These processes rely on specific base pairing, either with the nascent transcript (TGS) or with the mature target mRNA or a sequence-like mRNA (PTGS). While 21 nt siRNAs are predominantly involved in PTGS, 24 nt siRNAs often trigger TGS via epigenetic changes. In PTGS the target gene is directly known based on the sequence, whereas TGS can also affect genes located in close or distant proximity to the gene with sequence homology, making sequence-based prediction of TGS induced effects more difficult. PTGS is by far the most investigated mechanism in scientific studies to date, while the literature on RNAi-induced TGS is relatively limited. Therefore, the sheer number of scientific studies and the focus on PTGS to date does not necessarily reflect the actual relevance of the respective mechanisms in the plant, making it difficult to rank them according to their potential for causing unintended effects in plants.
Scientific literature on case studies investigating unintended effects in RNAi-based GM crops is currently scarce. However, bioinformatic tools are being dynamically developed to predict intended target genes and potential unintended effects on off-target genes in the GM crop or NTOs in case of HIGS, leveraging sequence homology to enhance the accuracy and scope of these predictions (Chen et al., 2019). While these tools often reach their limits in NTOs due to the lack or deficient annotated-genomes, high quality annotations are available for model plants or major crops, enabling such tools to predict PTGS effects on plant off-target genes with a higher probability. However, there are also mechanisms (such as TGS) that are not based on direct sequence homology to the target and whose unintended effects cannot be easily predicted bioinformatically. In most cases, it can be assumed that off-target effects manifest themselves at the transcriptome level and can be measured using sufficiently sensitive methods.
Adequate techniques to detect unintended RNAi-induced effects
To detect unintended RNAi-induced effects in GM plants for risk assessment, knowledge about the siRNAs processed in the GM plant, such as size and sequence, compared to the wild type is necessary. Since both intended and possible secondary siRNAs (such as tasiRNA) can play a role, untargeted analyses, such as small RNA sequencing, should be used to identify the sequences of all siRNAs. With this knowledge, bioinformatic tools can be used to predict both intended and unintended effects mediated by sequence homology, primarily through PTGS, and these predictions can be validated using targeted methods such as RT-qPCR. However, a comprehensive bioinformatic search for homologies requires access to the plant´s complete genome, whereas RT-qPCR analyses can also be managed with knowledge of shorter sequence segments. Unintended effects mediated by TGS, on the other hand, are not directly linked to the actual sequence of the siRNA and therefore cannot be adequately detected with targeted methods, but only with untargeted methods. RNA sequencing, for example, can be used for the direct, untargeted investigation of gene expression changes, changes in histone modifications can be detected using ChIP-Seq or altered DNA methylation patterns can be detected using bisulphite sequencing. However, all these methods require the availability of the plant´s genome for accurate analysis. Additionally, there are currently no studies that specifically address the importance of selecting appropriate plant material such as tissue type, developmental stage, and sampling time points or the sensitivities required for untargeted analyses to effectively capture RNAi-induced changes (e.g. alterations in gene expression). Most published studies have focused on using plant tissues, like leaves, without a detailed exploration on how these factors might influence the detection and interpretation of RNAi-induced effects. Likewise, unintended off-target genes may be expressed, for example, in certain tissue types and not in others. These gaps highlight the need for more comprehensive research to optimize experimental designs in RNAi studies aiming to identify unintended effects.
Case studies assessing unintended effects
Among the few studies assessing unintended effects of RNAi in GM plants, some have employed untargeted omics methods to analyze changes in gene expression and metabolite profiles. For example, Huang et al. (2022) compared the leaves of three transgenic maize RNAi lines resistant to Apolygus lucorum with those of three conventionally bred maize lines. Using untargeted omics methods at the levels of small RNAs, the transcriptome and the metabolome, the authors observed that the number of differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were greater in RNAi lines than in conventional lines. Additionally, Zörb et al. (2013) using GC-MS-based metabolite profiling showed that RNAi-mediated silencing of the sulfur-rich alpha-gliadin storage protein family in wheat grains did not induce changes in any of the 109 metabolites analyzed. Similarly, Zhang et al. (2020) investigated transcriptomic and metabolomic changes in RNAi-based GM maize resistant to Monolepta hieroglyphica compared to its unmodified variant. This study only identified a single DEG at the transcriptome level and 8 out of 5787 metabolites as DAMs, leading the authors to conclude that the RNAi variant exhibited negligible changes compared to the wild type.
Building on the insights gained from studies exploring off-target effects in RNAi-based GM plants, these findings have helped to inform regulatory approaches, including the one of the first authorization-relevant risk assessments for an RNAi-based genetically modified crop was carried out by the US Environmental Protection Agency (US EPA) for SmartStax Pro (MON 87411/Unique ID: MON-87411-9) (EPA Reg. Number: 62719-707). As part of the product characterization and human risk assessment, in 2016 the US EPA recommended a number of methods to rule out unintended side effects. These include transcriptome analyses using microarray or RNA sequencing, proteome analyses, GC-MS-based metabolomics, and the global detection of changes in DNA methylation patterns. It should be noted that certain recommended methods, such as microarray analyses for transcriptome studies or 2D gel electrophoresis coupled with MS for transcriptome analysis, are no longer state-of-the-art and should be replaced by more up-to-date methods such as RNA sequencing and LC-coupled MS, respectively. The US EPA advised that these analyses should be carried out comparatively between the GM plant containing all modification events (SmartStax Pro), the GM plant lacking the dsRNA cassette (SmartStax) as well as non-genetically modified lines across several generations. Furthermore, they recommended using a combination of different omics methods and to combine them with more sensitive methods such as RT-qPCR, to thoroughly exclude unintended effects.
Current limitations and future research
In summary, the challenges in detecting unintended RNAi effects in GM plants lie in the diversity of siRNAs that can be formed from corresponding precursor molecules and in the fact that TGS (especially via epigenetic mechanisms) can also affect the expression of nearby genes without sequence homology, indicating that targeted/biased bioinformatic methods alone are not sufficient for excluding unintended effects. The few available studies indicate that the RNAi method appears to be relatively specific with minimal unintended effects expected (Zörb et al., 2013; Zhang et al., 2020; Huang et al., 2024).
Untargeted approaches, such as RNA sequencing for transcriptome analysis, LC-MS-based proteomics or GC-MS-based metabolome profiling, offer a promising and increasingly sensitive means of investigating these effects. The current state of well-annotated plant genomes varies significantly across species, with high-quality annotations available for some model plants and major crops, while others remain underrepresented. This variability poses challenges for accurately mapping RNAi-induced changes, as comprehensive and well-annotated reference genomes are crucial for identifying both target and off-target effects, as well as for understanding the broader biological impact of RNAi in diverse plant species. One way around this problem is to perform a de novo assembly of the transcriptome of unannotated plants (Surget-Groba and Montoya-Burgos, 2010; Narzisi and Mishra, 2011). However, this depends on the quality and depth of the sequencing. In combination, bioinformatic approaches with untargeted methods, such as various omics, offer the possibility to detect specific off-target effects in GM plants.
Future research on detecting RNAi-induced effects in GM plants should focus on improving sensitivity and specificity with advanced sequencing technologies, better off-target detection through CRISPR, and more accurate quantification using methods like RT-qPCR and proteomics. Environmental impact studies, long-term monitoring, and standardizing protocols will be key for regulatory safety assessments.
Statements
Author contributions
CD: Data curation, Formal Analysis, Investigation, Writing – original draft. SA: Data curation, Investigation, Writing – review & editing. SS: Conceptualization, Writing – review & editing. LZ: Conceptualization, Writing – review & editing. AS: Data curation, Investigation, Writing – review & editing. EE: Data curation, Investigation, Writing – review & editing. SE: Data curation, Formal Analysis, Funding acquisition, Investigation, Supervision, Writing – original draft.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The project was funded by the German Federal Ministry for the Environment, Nature Conservation, Nuclear Safety and Consumer Protection and commissioned by the German Federal Agency for Nature Conservation (BfN); project number 3522 84 2100.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AllenE.XieZ.GustafsonA. M.CarringtonJ. C. (2005). microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell121, 207–221. doi: 10.1016/j.cell.2005.04.004
2
AsanoT.NishiuchiT. (2011). Comparative analysis of phosphoprotein expression using 2D-DIGE. 1940-6029744, 225–233. doi: 10.1007/978-1-61779-123-9_16
3
AugustineR.MukhopadhyayA.BishtN. C. (2013). Targeted silencing of BjMYB28 transcription factor gene directs development of low glucosinolate lines in oilseed Brassica juncea. Plant Biotechnol. J.11, 855–866. doi: 10.1111/pbi.12078
4
AxtellM. J.JanC.RajagopalanR.BartelD. P. (2006). A two-hit trigger for siRNA biogenesis in plants. Cell127, 565–577. doi: 10.1016/j.cell.2006.09.032
5
BaulcombeD. (2004). RNA silencing in plants. Nature431, 356–363. doi: 10.1038/nature02874
6
BaysoyA.BaiZ.SatijaR.FanR. (2023). The technological landscape and applications of single-cell multi-omics. Nat. Rev. Mol. Cell Biol.24, 695–713. doi: 10.1038/s41580-023-00615-w
7
BettiF.Ladera-CarmonaM. J.WeitsD. A.FerriG.IacopinoS.NoviG.et al. (2021). Exogenous miRNAs induce post-transcriptional gene silencing in plants. Nat. Plants7, 1379–1388. doi: 10.1038/s41477-021-01005-w
8
BlevinsT.RajeswaranR.ShivaprasadP. V.BeknazariantsD.Si-AmmourA.ParkH.-S.et al. (2006). Four plant Dicers mediate viral small RNA biogenesis and DNA virus induced silencing. Nucleic Acids Res.34, 6233–6246. doi: 10.1093/nar/gkl886
9
BockC.DatlingerP.ChardonF.CoelhoM. A.DongM. B.LawsonK. A.et al. (2022). High-content CRISPR screening. Nat. Rev. Methods Primers2, 1–23. doi: 10.1038/s43586-021-00093-4
10
BorgesF.MartienssenR. A. (2015). The expanding world of small RNAs in plants. Nat. Rev. Mol. Cell Biol.16, 727–741. doi: 10.1038/nrm4085
11
BorsaniO.ZhuJ.VersluesP. E.SunkarR.ZhuJ.-K. (2005). Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell123, 1279–1291. doi: 10.1016/j.cell.2005.11.035
12
BouchéN.LauresserguesD.GasciolliV.VaucheretH. (2006). An antagonistic function for Arabidopsis DCL2 in development and a new function for DCL4 in generating viral siRNAs. EMBO J.25, 3347–3356. doi: 10.1038/sj.emboj.7601217
13
BressanD.BattistoniG.HannonG. J. (2023). The dawn of spatial omics. Science381, eabq4964. doi: 10.1126/science.abq4964
14
BrodersenP.Sakvarelidze-AchardL.Bruun-RasmussenM.DunoyerP.YamamotoY. Y.SieburthL.et al. (2008). Widespread translational inhibition by plant miRNAs and siRNAs. Science320, 1185–1190. doi: 10.1126/science.1159151
15
BudnickA.FranklinM. J.UtleyD.EdwardsB.CharlesM.HornsteinE. D.et al. (2024). Long- and short-read sequencing methods discover distinct circular RNA pools in Lotus japonicus. Plant Genome17, e20429. doi: 10.1002/tpg2.20429
16
Burch-SmithT. M.AndersonJ. C.MartinG. B.Dinesh-KumarS. P. (2004). Applications and advantages of virus-induced gene silencing for gene function studies in plants. Plant J.39, 734–746. doi: 10.1111/j.1365-313X.2004.02158.x
17
CazaresT.HiggsR. E.WangJ.OzerH. G. (2024). SeedMatchR: identify off-target effects mediated by siRNA seed regions in RNA-seq experiments. Bioinformatics40, 1–4. doi: 10.1093/bioinformatics/btae011
18
ChanS. W.-L.HendersonI. R.ZhangX.ShahG.ChienJ. S.-C.JacobsenS. E. (2006). RNAi, DRD1, and histone methylation actively target developmentally important non-CG DNA methylation in arabidopsis. PloS Genet.2, e83. doi: 10.1371/journal.pgen.0020083
19
ChaudharyD.JeenaA. S.RohitGaurS.RajR.MishraS.et al. (2024). Advances in RNA interference for plant functional genomics: unveiling traits, mechanisms, and future directions. Appl. Biochem. Biotechnol. 196, 5681–5710. doi: 10.1007/s12010-023-04850-x
20
ChenH.-M.ChenL.-T.PatelK.LiY.-H.BaulcombeD. C.WuS.-H. (2010). 22-Nucleotide RNAs trigger secondary siRNA biogenesis in plants. Proc. Natl. Acad. Sci. U.S.A.107, 15269–15274. doi: 10.1073/pnas.1001738107
21
ChenL.HeikkinenL.WangC.YangY.SunH.WongG. (2019). Trends in the development of miRNA bioinformatics tools. Brief Bioinform.20, 1836–1852. doi: 10.1093/bib/bby054
22
ChenY.PangQ.-Y.HeY.ZhuN.BranstromI.YanX.-F.et al. (2012). Proteomics and metabolomics of Arabidopsis responses to perturbation of glucosinolate biosynthesis. Mol. Plant5, 1138–1150. doi: 10.1093/mp/sss034
23
ChiY. H.MoonJ. C.ParkJ. H.KimH.-S.ZulfugarovI. S.FanataW. I.et al. (2008). Abnormal chloroplast development and growth inhibition in rice thioredoxin m knock-down plants. Plant Physiol.148, 808–817. doi: 10.1104/pp.108.123547
24
ChristiaensO.DzhambazovaT.KostovK.ArpaiaS.JogaM. R.UrruI.et al. (2018). Literature review of baseline information on RNAi to support the environmental risk assessment of RNAi-based GM plants. EFS315, 1–173. doi: 10.2903/sp.efsa.2018.EN-1424
25
ChristiaensO.SweetJ.DzhambazovaT.UrruI.SmaggheG.KostovK.et al. (2022). Implementation of RNAi-based arthropod pest control: environmental risks, potential for resistance and regulatory considerations. J. Of Pest Sci.95, 1–15. doi: 10.1007/s10340-021-01439-3
26
CzarneckiO.BryanA. C.JawdyS. S.YangX.ChengZ.-M.ChenJ.-G.et al. (2016). Simultaneous knockdown of six non-family genes using a single synthetic RNAi fragment in Arabidopsis thaliana. Plant Methods12, 16. doi: 10.1186/s13007-016-0116-8
27
DongY.GaoQ.ChenY.ZhangZ.DuY.LiuY.et al. (2023). Identification of CircRNA signature associated with tumor immune infiltration to predict therapeutic efficacy of immunotherapy. Nat. Commun.14, 2540. doi: 10.1038/s41467-023-38232-y
28
FarooqR.HussainK.BashirA.RashidK.AshrafM. (2021). Databases and bioinformatics tools for genome engineering in plants using RNA interferenceNanobiotechnology for Plant Protection. 773–786. doi: 10.1016/B978-0-12-821910-2.00023-0
29
Food and Agriculture Organization of the United Nations (2003). GUIDELINE FOR THE CONDUCT OF FOOD SAFETY ASSESSMENT OF FOODS DERIVED FROM RECOMBINANT-DNA PLANTS. Available online at: https://www.fao.org/fileadmin/user_upload/gmfp/docs/CAC.GL_45_2003.pdf (Accessed December 9, 2024).
30
FukuharaT.UrayamaS.OkadaR.KiyotaE.MoriyamaH. (2011). Detection of long and short double-stranded RNAs. Springer Protocols744, 129–144. doi: 10.1007/978-1-61779-123-9_9
31
GaffoE.BuratinA.Dal MolinA.BortoluzziS. (2022). Sensitive, reliable and robust circRNA detection from RNA-seq with CirComPara2. Brief Bioinform.23, 1–12. doi: 10.1093/bib/bbab418
32
Garcia-RuizH.TakedaA.ChapmanE. J.SullivanC. M.FahlgrenN.BrempelisK. J.et al. (2010). Arabidopsis RNA-dependent RNA polymerases and dicer-like proteins in antiviral defense and small interfering RNA biogenesis during Turnip Mosaic Virus infection. Plant Cell22, 481–496. doi: 10.1105/tpc.109.073056
33
GoodR. T.VargheseT.GolzJ. F.RussellD. A.PapanicolaouA.EdwardsO.et al. (2016). OfftargetFinder: a web tool for species-specific RNAi design. Bioinformatics32, 1232–1234. doi: 10.1093/bioinformatics/btv747
34
GuangS.BochnerA. F.BurkhartK. B.BurtonN.PavelecD. M.KennedyS. (2010). Small regulatory RNAs inhibit RNA polymerase II during the elongation phase of transcription. Nature465, 1097–1101. doi: 10.1038/nature09095
35
GuoQ.LiuQ.SmithN. A.LiangG.WangM.-B. (2016). RNA silencing in plants: mechanisms, technologies and applications in horticultural crops. Curr. Genomics17, 476–489. doi: 10.2174/1389202917666160520103117
36
HanH. (2018). “RNA interference to knock down gene expression,” in Disease gene identification: methods and protocols. Ed. DiStefanoJ. K. (Springer, New York, NY), 293–302.
37
HaqueN.NishiguchiM. (2011). Bisulfite sequencing for cytosine-methylation analysis in plants. 1940-6029744, 187–197. doi: 10.1007/978-1-61779-123-9_13
38
HaveckerE. R.WallbridgeL. M.HardcastleT. J.BushM. S.KellyK. A.DunnR. M.et al. (2010). The arabidopsis RNA-directed DNA methylation argonautes functionally diverge based on their expression and interaction with target loci. Plant Cell22, 321–334. doi: 10.1105/tpc.109.072199
39
HeG.EllingA. A.DengX. W. (2011). The epigenome and plant development. Annu. Rev. Plant Biol.62, 411–435. doi: 10.1146/annurev-arplant-042110-103806
40
HelliwellC.WaterhouseP. (2003). Constructs and methods for high-throughput gene silencing in plants. Methods30, 289–295. doi: 10.1016/s1046-2023(03)00036-7
41
HendersonI. R.ZhangX.LuC.JohnsonL.MeyersB. C.GreenP. J.et al. (2006). Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning. Nat. Genet.38, 721–725. doi: 10.1038/ng1804
42
HeoJ. B.SungS. (2011). Vernalization-mediated epigenetic silencing by a long intronic noncoding RNA. Science331, 76–79. doi: 10.1126/science.1197349
43
Hernández-SotoA.Chacón-CerdasR. (2021). RNAi crop protection advances. Int. J. Mol. Sci.22, 1–15. doi: 10.3390/ijms222212148
44
HolmesK.WilliamsC. M.ChapmanE. A.CrossM. J. (2010). Detection of siRNA induced mRNA silencing by RT-qPCR: considerations for experimental design. BMC Res. Notes3, 53. doi: 10.1186/1756-0500-3-53
45
HuangC.WangZ.ZhuP.WangC.WangC.XuW.et al. (2022). RNA interference-based genetic engineering maize resistant to apolygus lucorum does not manifest unpredictable unintended effects relative to conventional breeding: short interfering RNA, transcriptome, and metabolome analysis. Front. Plant Sci.13. doi: 10.3389/fpls.2022.745708
46
HuangC.WangZ.ZhuP.XuW.DuZ.QianY.et al. (2024). RNAi-based genetically engineered rice resistant to black-streaked dwarf virus does not show adverse genetic effects: A multi-omics analysis. Plants People Planet6, 622–639. doi: 10.1002/ppp3.10477
47
HungY.-H.SlotkinR. K. (2021). The initiation of RNA interference (RNAi) in plants. Curr. Opin. Plant Biol.61, 102014. doi: 10.1016/j.pbi.2021.102014
48
JenC.-H.MichalopoulosI.WestheadD. R.MeyerP. (2005). Natural antisense transcripts with coding capacity in Arabidopsismay have a regulatory role that is not linked to double-stranded RNA degradation. Genome Biol.6, R51. doi: 10.1186/gb-2005-6-6-r51
49
JiaoS.WuS.HuangS.LiuM.GaoB. (2021). Advances in the identification of circular RNAs and research into circRNAs in human diseases. Front. In Genet.12, 665233. doi: 10.3389/fgene.2021.665233
50
Jones-RhoadesM. W.BartelD. P.BartelB. (2006). MicroRNAS and their regulatory roles in plants. Annu. Rev. Plant Biol.57, 19–53. doi: 10.1146/annurev.arplant.57.032905.105218
51
KanchiswamyC. N.MaffeiM.MalnoyM.VelascoR.KimJ.-S. (2016). Fine-tuning next-generation genome editing tools. Trends Biotechnol.34, 562–574. doi: 10.1016/j.tibtech.2016.03.007
52
KeykhaF.BagheriA.MoshtaghiN.BahramiA. R.SharifiA. (2016). RNAi-induced silencing in floral tissues of Petunia hybrida by agroinfiltration: a rapid assay for chalcone isomerase gene function analysis. Cell Mol. Biol. (Noisy-le-grand)62, 26–31. doi: 10.14715/cmb/2016.62.10.4
53
KimV. N.HanJ.SiomiM. C. (2009). Biogenesis of small RNAs in animals. Nat. Rev. Mol. Cell Biol.10, 126–139. doi: 10.1038/nrm2632
54
KimE.-D.SungS. (2012). Long noncoding RNA: unveiling hidden layer of gene regulatory networks. Trends Plant Sci.17, 16–21. doi: 10.1016/j.tplants.2011.10.008
55
KitzmannP.SchwirzJ.Schmitt-EngelC.BucherG. (2013). RNAi phenotypes are influenced by the genetic background of the injected strain. BMC Genomics14, 5. doi: 10.1186/1471-2164-14-5
56
KleterG. A. (2020). Food safety assessment of crops engineered withRNAinterference and other methods to modulate expression of endogenous and plant pest genes. Pest Manag Sci.76, 3333–3339. doi: 10.1002/ps.5957
57
KlocA.ZaratieguiM.NoraE.MartienssenR. (2008). RNA interference guides histone modification during the S phase of chromosomal replication. Curr. Biol.18, 490–495. doi: 10.1016/j.cub.2008.03.016
58
KoeppeS.KawchukL.KalischukM. (2023). RNA interference past and future applications in plants. Int. J. Mol. Sci.24, 1–18. doi: 10.3390/ijms24119755
59
KrzysztonM.KufelJ. (2022). Analysis of mRNA-derived siRNAs in mutants of mRNA maturation and surveillance pathways in Arabidopsis thaliana. Sci. Rep.12, 1474. doi: 10.1038/s41598-022-05574-4
60
KumarR.ConklinD. S.MittalV. (2003). High-throughput selection of effective RNAi probes for gene silencing. Genome Res.13, 2333–2340. doi: 10.1101/gr.1575003
61
KuriharaY.WatanabeY. (2004). Arabidopsis micro-RNA biogenesis through Dicer-like 1 protein functions. Proc. Natl. Acad. Sci. U.S.A.101, 12753–12758. doi: 10.1073/pnas.0403115101
62
KyslíkJ.Born-TorrijosA.HolzerA. S.KosakyanA. (2024). RNAi-directed knockdown in the cnidarian fish blood parasite Sphaerospora molnari. Sci. Rep.14, 3545. doi: 10.1038/s41598-024-54171-0
63
LacourseE. J.PerallyS.Hernandez-ViadelM.WrightH. A.BrophyP. M. (2008). A proteomics approach to quantify protein levels following RNA interference: case study with glutathione transferase superfamily from the model metazoan Caenorhabditis elegans. J. Proteome Res.7, 3314–3318. doi: 10.1021/pr8001035
64
LawJ. A.JacobsenS. E. (2010). Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat. Rev. Genet.11, 204–220. doi: 10.1038/nrg2719
65
LewseyM. G.HardcastleT. J.MelnykC. W.MolnarA.ValliA.UrichM. A.et al. (2016). Mobile small RNAs regulate genome-wide DNA methylation. Proc. Natl. Acad. Sci. U.S.A.113, E801–E810. doi: 10.1073/pnas.1515072113
66
LiangW.MasonA. J.LamJ. K. W. (2013). “Western Blot Evaluation of siRNA Delivery by pH-Responsive Peptides,” in Target identification and validation in drug discovery: methods and protocols. Ed. Moll (Humana Press Incorporated, Totowa, NJ), 73–87.
67
LiuJ.-X.ChiouC.-Y.ShenC.-H.ChenP.-J.LiuY.-C.JianC.-D.et al. (2014). RNA interference-based gene silencing of phytoene synthase impairs growth, carotenoids, and plastid phenotype in Oncidium hybrid orchid. Springerplus3, 478. doi: 10.1186/2193-1801-3-478
68
LiuS.JaouannetM.DempseyD. A.ImaniJ.CoustauC.KogelK.-H. (2020). RNA-based technologies for insect control in plant production. Biotechnol. Adv.39, 107463. doi: 10.1016/j.bioteChadv.2019.107463
69
LiuJ.JungC.XuJ.WangH.DengS.BernadL.et al. (2012). Genome-wide analysis uncovers regulation of long intergenic noncoding RNAs in Arabidopsis. Plant Cell24, 4333–4345. doi: 10.1105/tpc.112.102855
70
López-MárquezD.Del-EspinoÁ.Ruiz-AlbertJ.BejaranoE. R.BrodersenP.BeuzónC. R. (2023). Regulation of plant immunity via small RNA-mediated control of NLR expression. J. Exp. Bot.74, 6052–6068. doi: 10.1093/jxb/erad268
71
LuR.Martin-HernandezA. M.PeartJ. R.MalcuitI.BaulcombeD. C. (2003). Virus-induced gene silencing in plants. Methods30, 296–303. doi: 10.1016/s1046-2023(03)00037-9
72
LückS.KresziesT.StrickertM.SchweizerP.KuhlmannM.DouchkovD. (2019). siRNA-finder (si-fi) software for RNAi-target design and off-target prediction. Front. Plant Sci.10. doi: 10.3389/fpls.2019.01023
73
MalloryA. C.VaucheretH. (2006). Functions of microRNAs and related small RNAs in plants. Nat. Genet.38 Suppl, S31–S36. doi: 10.1038/ng1791
74
ManskeU.LandsmannJ.Dietz-PfeilstetterA. (2017). Comparison of different methods for the establishment of RNA silencing in plants. Plant Biotechnol. Rep.11, 115–125. doi: 10.1007/s11816-017-0436-9
75
MatzkeM. A.MosherR. A. (2014). RNA-directed DNA methylation: an epigenetic pathway of increasing complexity. Nat. Rev. Genet.15, 394–408. doi: 10.1038/nrg3683
76
MooreJ. D. (2015). The impact of CRISPR-Cas9 on target identification and validation. Drug Discovery Today20, 450–457. doi: 10.1016/j.drudis.2014.12.016
77
MovahediA.SunW.ZhangJ.WuX.MousaviM.MohammadiK.et al. (2015). RNA-directed DNA methylation in plants. Plant Cell Rep.34, 1857–1862. doi: 10.1007/s00299-015-1839-0
78
MuhammadI. I.KongS. L.Akmar AbdullahS. N.MunusamyU. (2020). RNA-seq and chIP-seq as complementary approaches for comprehension of plant transcriptional regulatory mechanism. IJMS21, 167. doi: 10.3390/ijms21010167
79
MujtabaM.WangD.CarvalhoL. B.OliveiraJ. L.Espirito Santo PereiraA.SharifR.et al. (2021). Nanocarrier-mediated delivery of miRNA, RNAi, and CRISPR-cas for plant protection: current trends and future directions. ACS Agric. Sci. Technol.1, 417–435. doi: 10.1021/acsagscitech.1c00146
80
MukherjeeK.CamposH.KolaczkowskiB. (2013). Evolution of animal and plant dicers: early parallel duplications and recurrent adaptation of antiviral RNA binding in plants. Mol. Biol. Evol.30, 627–641. doi: 10.1093/molbev/mss263
81
NaikB.KumarV.RizwanuddinS.ChauhanM.ChoudharyM.GuptaA. K.et al. (2023). Genomics, proteomics, and metabolomics approaches to improve abiotic stress tolerance in tomato plant. Int. J. Mol. Sci.24. doi: 10.3390/ijms24033025
82
NarzisiG.MishraB. (2011). Comparing de novo genome assembly: the long and short of it. PloS One6, e19175. doi: 10.1371/journal.pone.0019175
83
Navarro-MendozaM. I.Pérez-ArquesC.HeitmanJ. (2023). Heterochromatin and RNAi act independently to ensure genome stability in Mucorales human fungal pathogens. Proc. Natl. Acad. Sci. U.S.A.120, 1–12. doi: 10.1073/pnas.2220475120
84
NguyenM. H.NguyenH.-N.VuT. N. (2022). Evaluation of methods to detect circular RNAs from single-end RNA-sequencing data. BMC Genomics23, 106. doi: 10.1186/s12864-022-08329-7
85
NowaraD.GayA.LacommeC.ShawJ.RidoutC.DouchkovD.et al. (2010). HIGS: host-induced gene silencing in the obligate biotrophic fungal pathogen Blumeria graminis. Plant Cell22, 3130–3141. doi: 10.1105/tpc.110.077040
86
Olmedo-MonfilV.Durán-FigueroaN.Arteaga-VázquezM.Demesa-ArévaloE.AutranD.GrimanelliD.et al. (2010). Control of female gamete formation by a small RNA pathway in Arabidopsis. Nature464, 628–632. doi: 10.1038/nature08828
87
Organisation for Economic Co-operation and Development (2021). REVISED CONSENSUS DOCUMENT ON COMPOSITIONAL CONSIDERATIONS FOR NEW VARIETIES OF POTATO (Solanum tuberosum): key food and feed nutrients, toxicants, allergens, anti-nutrients and other plant metabolite: ENV/JM/MONO(2020). Paris: OECD.
88
Organisation for Economic Co-operation and Development (OECD). (2020). Considerations for the environmental risk assessment of the application of sprayed or externally applied ds-RNA-based pesticides (OECD).
89
PapadopoulouN.DevosY.Alvarez-AlfagemeF.LanzoniA.WaigmannE. (2020). Risk assessment considerations for genetically modified RNAi plants: EFSA’s activities and perspective. Front. Plant Sci.11. doi: 10.3389/fpls.2020.00445
90
PeretzL.BesserE.HajbiR.CasdenN.ZivD.KronenbergN.et al. (2018). Combined shRNA over CRISPR/cas9 as a methodology to detect off-target effects and a potential compensatory mechanism. Sci. Rep.8, 93. doi: 10.1038/s41598-017-18551-z
91
QiY.HeX.WangX.-J.KohanyO.JurkaJ.HannonG. J. (2006). Distinct catalytic and non-catalytic roles of ARGONAUTE4 in RNA-directed DNA methylation. Nature443, 1008–1012. doi: 10.1038/nature05198
92
QinC.LiB.FanY.ZhangX.YuZ.RyabovE.et al. (2017). Roles of dicer-like proteins 2 and 4 in intra- and intercellular antiviral silencing. Plant Physiol.174, 1067–1081. doi: 10.1104/pp.17.00475
93
QuF.YeX.MorrisT. J. (2008). Arabidopsis DRB4, AGO1, AGO7, and RDR6 participate in a DCL4-initiated antiviral RNA silencing pathway negatively regulated by DCL1. Proc. Natl. Acad. Sci. U.S.A.105, 14732–14737. doi: 10.1073/pnas.0805760105
94
RinaldiG.MoralesM. E.CancelaM.CastilloE.BrindleyP. J.TortJ. F. (2008). Development of functional genomic tools in trematodes: RNA interference and luciferase reporter gene activity in Fasciola hepatica. PloS Negl. Trop. Dis.2, e260. doi: 10.1371/journal.pntd.0000260
95
RuizM. T.VoinnetO.BaulcombeD. C. (1998). Initiation and maintenance of virus-induced gene silencing. Plant Cell10, 937–946. doi: 10.1105/tpc.10.6.937
96
SahinO.LöbkeC.KorfU.AppelhansH.SültmannH.PoustkaA.et al. (2007). Combinatorial RNAi for quantitative protein network analysis. Proc. Natl. Acad. Sci. U.S.A.104, 6579–6584. doi: 10.1073/pnas.0606827104
97
Sanan-MishraN.Abdul Kader JailaniA.MandalB.MukherjeeS. K. (2021). Secondary siRNAs in plants: biosynthesis, various functions, and applications in virology. Front. Plant Sci.12, 610283. doi: 10.1126/science.1079695
98
SarkarA.Roy-BarmanS. (2021). Spray-induced silencing of pathogenicity gene moDES1 via exogenous double-stranded RNA can confer partial resistance against fungal blast in rice. Front. Plant Sci.12. doi: 10.3389/fpls.2021.733129
99
SarkiesP.MiskaE. A. (2014). Small RNAs break out: the molecular cell biology of mobile small RNAs. Nat. Rev. Mol. Cell Biol.15, 525–535. doi: 10.1038/nrm3840
100
Senthil-KumarM.MysoreK. S. (2011). Caveat of RNAi in plants: the off-target effect. 1940-6029744, 13–25. doi: 10.1007/978-1-61779-123-9_2
101
SmartN.ScamblerP. J.RileyP. R. (2005). A rapid and sensitive assay for quantification of siRNA efficiency and specificity. Biol. Procedures Online7, 1–7. doi: 10.1251/bpo99
102
SunG.RossiJ. J. (2009). Problems associated with reporter assays in RNAi studies. RNA Biol.6, 406–411. doi: 10.4161/rna.6.4.9218
103
Surget-GrobaY.Montoya-BurgosJ. I. (2010). Optimization of de novo transcriptome assembly from next-generation sequencing data. Genome Res.20, 1432–1440. doi: 10.1101/gr.103846.109
104
SweversL.Vanden BroeckJ.SmaggheG. (2013). The possible impact of persistent virus infection on the function of the RNAi machinery in insects: a hypothesis. Front. Physiol.4. doi: 10.3389/fphys.2013.00319
105
TaoY.ChiuL.-W.HoyleJ. W.DewhirstR. A.RicheyC.RasmussenK.et al. (2023). Enhanced photosynthetic efficiency for increased carbon assimilation and woody biomass production in engineered hybrid poplar. Forests14, 827. doi: 10.3390/f14040827
106
TyagiP.SinghD.MathurS.SinghA.RanjanR. (2022). Upcoming progress of transcriptomics studies on plants: An overview. Front. Plant Sci.13. doi: 10.3389/fpls.2022.1030890
107
Valencia-SanchezM. A.LiuJ.HannonG. J.ParkerR. (2006). Control of translation and mRNA degradation by miRNAs and siRNAs. Genes Dev.20, 515–524. doi: 10.1101/gad.1399806
108
Varkonyi-GasicE.HellensR. P. (2011). Quantitative stem-loop RT-PCR for detection of microRNAs. Springer Protocols744, 145–157. doi: 10.1007/978-1-61779-123-9_10
109
VazquezF.VaucheretH.RajagopalanR.LepersC.GasciolliV.MalloryA. C.et al. (2004). Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol. Cell16, 69–79. doi: 10.1016/j.molcel.2004.09.028
110
VerdelA.VavasseurA.Le GorrecM.Touat-TodeschiniL. (2009). Common themes in siRNA-mediated epigenetic silencing pathways. Int. J. Dev. Biol.53, 245–257. doi: 10.1387/ijdb.082691av
111
VidarsdottirL.GoroshchukO.KolosenkoI.Palm-ApergiC. (2019). Designing siRNA and Evaluating Its Effect on RNA Targets Using qPCR and Western Blot. Springer Protocols2036, 53–72. doi: 10.1007/978-1-4939-9670-4_3
112
VoinnetO. (2009). Origin, biogenesis, and activity of plant microRNAs. Cell136, 669–687. doi: 10.1016/j.cell.2009.01.046
113
WangX.-B.JovelJ.UdompornP.WangY.WuQ.LiW.-X.et al. (2011). The 21-nucleotide, but not 22-nucleotide, viral secondary small interfering RNAs direct potent antiviral defense by two cooperative argonautes in Arabidopsis thaliana. Plant Cell23, 1625–1638. doi: 10.1105/tpc.110.082305
114
WarnatzH.-J.SchmidtD.MankeT.PicciniI.SultanM.BorodinaT.et al. (2011). The BTB and CNC homology 1 (BACH1) target genes are involved in the oxidative stress response and in control of the cell cycle. J. Of Biol. Chem.286, 23521–23532. doi: 10.1074/jbc.M111.220178
115
WaterhouseP. M.HelliwellC. A. (2003). Exploring plant genomes by RNA-induced gene silencing. Nat. Rev. Genet.4, 29–38. doi: 10.1038/nrg982
116
WernerB. T.GaffarF. Y.SchuemannJ.BiedenkopfD.KochA. M. (2020). RNA-spray-mediated silencing of fusarium graminearum AGO and DCL genes improve barley disease resistance. Front. Plant Sci.11. doi: 10.3389/fpls.2020.00476
117
WesleyS. V.HelliwellC. A.SmithN. A.WangM. B.RouseD. T.LiuQ.et al. (2001). Construct design for efficient, effective and high-throughput gene silencing in plants. Plant J.27, 581–590. doi: 10.1046/j.1365-313x.2001.01105.x
118
WierzbickiA. T.HaagJ. R.PikaardC. S. (2008). Noncoding transcription by RNA polymerase Pol IVb/Pol V mediates transcriptional silencing of overlapping and adjacent genes. Cell135, 635–648. doi: 10.1016/j.cell.2008.09.035
119
WierzbickiA. T.ReamT. S.HaagJ. R.PikaardC. S. (2009). RNA polymerase V transcription guides ARGONAUTE4 to chromatin. Nat. Genet.41, 630–634. doi: 10.1038/ng.365
120
WuL.MaoL.QiY. (2012). Roles of dicer-like and argonaute proteins in TAS-derived small interfering RNA-triggered DNA methylation. Plant Physiol.160, 990–999. doi: 10.1104/pp.112.200279
121
WuL.ZhouH.ZhangQ.ZhangJ.NiF.LiuC.et al. (2010). DNA methylation mediated by a microRNA pathway. Mol. Cell38, 465–475. doi: 10.1016/j.molcel.2010.03.008
122
XieZ.AllenE.WilkenA.CarringtonJ. C. (2005). DICER-LIKE 4 functions in trans-acting small interfering RNA biogenesis and vegetative phase change in Arabidopsis thaliana. Proc. Natl. Acad. Sci. U.S.A.102, 12984–12989. doi: 10.1073/pnas.0506426102
123
XuK.ZhangX.-M.ChenH.ZhangC.ZhuJ.ChengZ.et al. (2021). Fine-tuning florigen increases field yield through improving photosynthesis in soybean. Front. Plant Sci.12. doi: 10.3389/fpls.2021.710754
124
ZapletalD.KubicekK.SvobodaP.SteflR. (2023). Dicer structure and function: conserved and evolving features. EMBO Rep.24, e57215. doi: 10.15252/embr.202357215
125
ZaratieguiM.MartienssenR. A. (2012). SnapShot: small RNA-mediated epigenetic modifications. Cell151, 456–456.e1. doi: 10.1016/j.cell.2012.10.001
126
ZhangX.XiaJ.LiiY. E.Barrera-FigueroaB. E.ZhouX.GaoS.et al. (2012). Genome-wide analysis of plant nat-siRNAs reveals insights into their distribution, biogenesis and function. Genome Biol.13, R20. doi: 10.1186/gb-2012-13-3-r20
127
ZhangX. L.ZhangR. Y.LiL.YangY.DingY. J.GuanH. T.et al. (2020). Negligible transcriptome and metabolome alterations in RNAi insecticidal maize againstMonolepta hieroglyphica. Plant Cell Rep.39, 1539–1547. doi: 10.1007/s00299-020-02582-4
128
ZhangH.ZhuJ.-K. (2011). RNA-directed DNA methylation. Curr. Opin. Plant Biol.14, 142–147. doi: 10.1016/j.pbi.2011.02.003
129
ZhengX.ZhuJ.KapoorA.ZhuJ.-K. (2007). Role of Arabidopsis AGO6 in siRNA accumulation, DNA methylation and transcriptional gene silencing. EMBO J.26, 1691–1701. doi: 10.1038/sj.emboj.7601603
130
ZhouH.YangM.ZhaoL.ZhuZ.LiuF.SunH.et al. (2021). HIGH-TILLERING AND DWARF 12 modulates photosynthesis and plant architecture by affecting carotenoid biosynthesis in rice. J. Exp. Bot.72, 1212–1224. doi: 10.1093/jxb/eraa497
131
ZilbermanD.CaoX.JacobsenS. E. (2003). ARGONAUTE4 control of locus-specific siRNA accumulation and DNA and histone methylation. Science299, 716–719. doi: 10.1126/science.1079695
132
ZörbC.BeckerD.HaslerM.MühlingK.GöddeV.NiehausK.et al. (2013). Silencing of the sulfur rich α-gliadin storage protein family in wheat grains (Triticum aestivum L.) causes no unintended side-effects on other metabolites. Front. Plant Sci.4. doi: 10.3389/fpls.2013.00369
Summary
Keywords
RNAi GM plants, detection techniques, RNAi mechanism, off-target effects, RNAi pest control
Citation
Diaz C, Ayobahan SU, Simon S, Zühl L, Schiermeyer A, Eilebrecht E and Eilebrecht S (2025) Classification of and detection techniques for RNAi-induced effects in GM plants. Front. Plant Sci. 16:1535384. doi: 10.3389/fpls.2025.1535384
Received
27 November 2024
Accepted
08 February 2025
Published
07 March 2025
Volume
16 - 2025
Edited by
Ruonan Ma, Zhengzhou University, China
Reviewed by
Ming Wang, Chinese Academy of Sciences (CAS), China
Hao Wang, Northwest A&F University, China
Updates
Copyright
© 2025 Diaz, Ayobahan, Simon, Zühl, Schiermeyer, Eilebrecht and Eilebrecht.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sebastian Eilebrecht, sebastian.eilebrecht@ime.fraunhofer.de
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.