Identification of RNA silencing components in soybean and sorghum

Function and regulation of plant ARGONAUTE proteins in response to environmental challenges: a review

Environmental stresses diversely affect multiple processes related to the growth, development, and yield of many crops worldwide. In response, plants have developed numerous sophisticated defense mechanisms at the cellular and subcellular levels to react and adapt to biotic and abiotic stressors. RNA silencing, which is an innate immune mechanism, mediates sequence-specific gene expression regulation in higher eukaryotes. ARGONAUTE (AGO) proteins are essential components of the RNA-induced silencing complex (RISC). They bind to small noncoding RNAs (sRNAs) and target complementary RNAs, causing translational repression or triggering endonucleolytic cleavage pathways. In this review, we aim to illustrate the recently published molecular functions, regulatory mechanisms, and biological roles of AGO family proteins in model plants and cash crops, especially in the defense against diverse biotic and abiotic stresses, which could be helpful in crop improvement and stress tolerance in various plants.

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.

Genome-wide analysis of key gene families in RNA silencing and their responses to biotic and drought stresses in adzuki bean

BACKGROUND

In plants, RNA silencing is an important conserved mechanism to regulate gene expression and combat against abiotic and biotic stresses. Dicer-like (DCL) and Argonaute (AGO) proteins and RNA-dependent RNA polymerase (RDR) are the core elements involved in gene silencing and their gene families have been explored in many plants. However, these genes and their responses to stresses have not yet been well characterized in adzuki bean.

RESULTS

A total of 11 AGO, 7 DCL and 6 RDR proteins were identified, and phylogenetic analyses of these proteins showed that they clustered into six, four and four clades respectively. The expression patterns of these genes in susceptible or resistant adzuki bean cultivars challenged with drought, bean common mosaic virus and Podosphaera xanthii infections were further validated by quantitative RT-PCR. The different responses of these proteins under abiotic and biotic stresses indicated their specialized regulatory mechanisms.

CONCLUSIONS

In this study, 24 genes of the DCL, AGO and RDR gene families in adzuki bean were identified, and the sequence characterization, structure of the encoded proteins, evolutionary relationship with orthologues in other legumes and gene expression patterns under drought and biotic stresses were primarily explored, which enriched our understanding of these genes in adzuki bean. Our findings provide a foundation for the comparative genomic analyses of RNA silencing elements in legume plants and further new insights into the functional complexity of RNA silencing in the response to various stresses in adzuki bean.

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.

Genome-Wide Identification of Key Components of RNA Silencing in Two Phaseolus vulgaris Genotypes of Contrasting Origin and Their Expression Analyses in Response to Fungal Infection

RNA silencing serves key roles in a multitude of cellular processes, including development, stress responses, metabolism, and maintenance of genome integrity. Dicer, Argonaute (AGO), double-stranded RNA binding (DRB) proteins, RNA-dependent RNA polymerase (RDR), and DNA-dependent RNA polymerases known as Pol IV and Pol V form core components to trigger RNA silencing. Common bean (Phaseolus vulgaris) is an important staple crop worldwide. In this study, we aimed to unravel the components of the RNA-guided silencing pathway in this non-model plant, taking advantage of the availability of two genome assemblies of Andean and Meso-American origin. We identified six PvDCLs, thirteen PvAGOs, 10 PvDRBs, 5 PvRDRs, in both genotypes, suggesting no recent gene amplification or deletion after the gene pool separation. In addition, we identified one PvNRPD1 and one PvNRPE1 encoding the largest subunits of Pol IV and Pol V, respectively. These genes were categorized into subgroups based on phylogenetic analyses. Comprehensive analyses of gene structure, genomic localization, and similarity among these genes were performed. Their expression patterns were investigated by means of expression models in different organs using online data and quantitative RT-PCR after pathogen infection. Several of the candidate genes were up-regulated after infection with the fungus Colletotrichum lindemuthianum.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Transcriptome sequencing analysis of sorghum callus with various regeneration capacities

The possible molecular mechanisms regulating sorghum callus regeneration were revealed by RNA-sequencing. Plant callus regeneration has been widely applied in agricultural improvement. Recently, callus regeneration has been successfully applied in the genetic transformation of sorghum by using immature sorghum embryos as explants. However, the mechanism underlying callus regeneration in sorghum is still largely unknown. Here, we describe three types of callus (Callus I-III) with different redifferentiation abilities undergoing distinct induction from immature embryos of the Hiro-1 variety. Compared with nonembryonic Callus III, Callus I produced only some identifiable roots, and embryonic Callus II was sufficient to regenerate whole plants. Genome-wide transcriptome profiles were generated to reveal the underlying mechanisms. The numbers of differentially expressed genes for the three types of callus varied from 5906 to 8029. In accordance with the diverse regeneration abilities observed for different types of callus and leaf tissues, the principal component analysis revealed that the gene expression patterns of Callus I and Callus II were different from those of Callus III and leaves regenerated from Callus II. Notably, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses, pharmacological treatment, and substance content determinations revealed that plant ribosomes, lignin metabolic processes, and metabolism of starch and sucrose were significantly enriched, suggesting that these factors are associated with callus regeneration. These results helped elucidate the molecular regulation of three types of callus with different regeneration abilities in sorghum.

Genome-wide analyses across Viridiplantae reveal the origin and diversification of small RNA pathway-related genes

Small RNAs play a major role in the post-transcriptional regulation of gene expression in eukaryotes. Despite the evolutionary importance of streptophyte algae, knowledge on small RNAs in this group of green algae is almost non-existent. We used genome and transcriptome data of 34 algal and plant species, and performed genome-wide analyses of small RNA (miRNA & siRNA) biosynthetic and degradation pathways. The results suggest that Viridiplantae started to evolve plant-like miRNA biogenesis and degradation after the divergence of the Mesostigmatophyceae in the streptophyte algae. We identified two major evolutionary transitions in small RNA metabolism in streptophyte algae; during the first transition, the origin of DCL-New, DCL1, AGO1/5/10 and AGO4/6/9 in the last common ancestor of Klebsormidiophyceae and all other streptophytes could be linked to abiotic stress responses and evolution of multicellularity in streptophytes. During the second transition, the evolution of DCL 2,3,4, and AGO 2,3,7 as well as DRB1 in the last common ancestor of Zygnematophyceae and embryophytes, suggests their possible contribution to pathogen defense and antibacterial immunity. Overall, the origin and diversification of DICER and AGO along with several other small RNA pathway-related genes among streptophyte algae suggested progressive adaptations of streptophyte algae during evolution to a subaerial environment.

Wang, Liang et al. conduct a genome-wide investigation into the origin of small RNA pathway-related genes in Viridiplantae. Their findings suggest that streptophyte algae progressively adapted to a subaerial environment through generation of these pathways.

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.

What does it take to be antiviral? An Argonaute-centered perspective on plant antiviral defense

RNA silencing is a major mechanism of constitutive antiviral defense in plants, mediated by a number of proteins, including the Dicer-like (DCL) and Argonaute (AGO) endoribonucleases. Both DCL and AGO protein families comprise multiple members. In particular, the AGO protein family has expanded considerably in different plant lineages, with different family members having specialized functions. Although the general mode of action of AGO proteins is well established, the properties that make different AGO proteins more or less efficient at targeting viruses are less well understood. In this report, we review methodologies used to study AGO antiviral activity and current knowledge about which AGO family members are involved in antiviral defense. In addition, we discuss what is known about the different properties of AGO proteins thought to be associated with this function.

LuluDB-The Database Created Based on Small RNA, Transcriptome, and Degradome Sequencing Shows the Wide Landscape of Non-coding and Coding RNA in Yellow Lupine (Lupinus luteus L.) Flowers and Pods

Yellow lupine (Lupinus luteus L.) belongs to a legume family that benefits from symbiosis with nitrogen-fixing bacteria. Its seeds are rich in protein, which makes it a valuable food source for animals and humans. Yellow lupine is also the model plant for basic research on nodulation or abscission of organs. Nevertheless, the knowledge about the molecular regulatory mechanisms of its generative development is still incomplete. The RNA-Seq technique is becoming more prominent in high-throughput identification and expression profiling of both coding and non-coding RNA sequences. However, the huge amount of data generated with this method may discourage other scientific groups from making full use of them. To overcome this inconvenience, we have created a database containing analysis-ready information about non-coding and coding L. luteus RNA sequences (LuluDB). LuluDB was created on the basis of RNA-Seq analysis of small RNA, transcriptome, and degradome libraries obtained from yellow lupine cv. Taper flowers, pod walls, and seeds in various stages of development, flower pedicels, and pods undergoing abscission or maintained on the plant. It contains sequences of miRNAs and phased siRNAs identified in L. luteus, information about their expression in individual samples, and their target sequences. LuluDB also contains identified lncRNAs and protein-coding RNA sequences with their organ expression and annotations to widely used databases like GO, KEGG, NCBI, Rfam, Pfam, etc. The database also provides sequence homology search by BLAST using, e.g., an unknown sequence as a query. To present the full capabilities offered by our database, we performed a case study concerning transcripts annotated as DCL 1–4 (DICER LIKE 1–4) homologs involved in small non-coding RNA biogenesis and identified miRNAs that most likely regulate DCL1 and DCL2 expression in yellow lupine. LuluDB is available at http://luluseqdb.umk.pl/basic/web/index.php.

Dynamic distribution of ARGONAUTE1 (AGO1) and ARGONAUTE4 (AGO4) in Hyacinthus orientalis L. pollen grains and pollen tubes growing in vitro

The transcriptional and posttranscriptional AGO-mediated control of gene expression may play important roles during male monocot gametophyte development. In this report, we demonstrated dynamic changes in the spatiotemporal distribution of AGO1 and AGO4, which are key proteins of the RNA-induced silencing complex (RISC) in Hyacinthus orientalis male gametophyte development. During maturation of the bicellular pollen grains and in vitro pollen tube growth, the pattern of AGO1 localization was correlated with previously observed transcriptional activity of the cells. During the period of high transcriptional activity, AGO1 is associated with chromatin while the clustered distribution of AGO1 in the interchromatin areas is accompanied by condensation of chromatin and the gradual transcriptional silencing of both cells in mature, dehydrated pollen. During pollen tube growth and the restarting of RNA synthesis in the vegetative nucleus, AGO1 is dispersed in the chromatin. Additionally, the gradual increase in the cytoplasmic pool of AGO1 in the elongating pollen tube indicates the activation of the posttranscriptional gene silencing (PTGS) pathway. During pollen tube growth in the generative cell and in the sperm cells, AGO1 is present mainly in the areas between highly condensed chromatin clusters. Changes in the distribution of AGO4 that indicated the possibility of spatiotemporal organization in the RNA-directed DNA methylation (RdDM) process (cytoplasmic and nuclear steps) were also observed during hyacinth male gametophyte development. Based on our findings, we propose that in the germinating pollen tube, the cytoplasmic assembly of AGO4/siRNA takes place and that the mature complexes could be transported to the nucleus to carry out their function during the next steps of pollen tube growth.

A putative AGO protein, OsAGO17, positively regulates grain size and grain weight through OsmiR397b in rice

Argonaute (AGO) proteins and small RNAs (sRNAs) are core components of the RNA-induced silencing complex (RISC). It has been reported that miRNAs regulate plant height and grain size in rice, but which AGO is involved in grain size regulation remains unclear. Here, we report that enhanced expression of OsAGO17, a putative AGO protein, could improve grain size and weight and promote stem development in rice. Cytological evidence showed that these effects are mainly caused by alteration of cell elongation. Expression analyses showed that OsAGO17 was highly expressed in young panicles and nodes, which was consistent with the expression pattern of OsmiR397b. SRNA sequencing, stem-loop RT-PCR and sRNA blotting showed that the expression of OsmiR397b was reduced in ago17 and enhanced in the OsAGO17 OE lines. Four OsmiR397b target laccase (LAC) genes showed complementary expression patterns with OsAGO17 and OsmiR397b. Combined with the results of immunoprecipitation (IP) analysis, we suggested that OsAGO17 formed an RISC with OsmiR397b and affected rice development by suppression of LAC expression. In conclusion, OsAGO17 might be a critical protein in the sRNA pathway and positively regulates grain size and weight in rice.

Argonaute Proteins: Why Are They So Important for the Legume-Rhizobia Symbiosis?

Unlike most other land plants, legumes can fulfill their nitrogen needs through the establishment of symbioses with nitrogen-fixing soil bacteria (rhizobia). Through this symbiosis, fixed nitrogen is incorporated into the food chain. Because of this ecological relevance, the genetic mechanisms underlying the establishment of the legume-rhizobia symbiosis (LRS) have been extensively studied over the past decades. During this time, different types of regulators of this symbiosis have been discovered and characterized. A growing number of studies have demonstrated the participation of different types of small RNAs, including microRNAs, in the different stages of this symbiosis. The involvement of small RNAs also indicates that Argonaute (AGO) proteins participate in the regulation of the LRS. However, despite this obvious role, the relevance of AGO proteins in the LRS has been overlooked and understudied. Here, we discuss and hypothesize the likely participation of AGO proteins in the regulation of the different steps that enable the establishment of the LRS. We also briefly review and discuss whether rhizobial symbiosis induces DNA damages in the legume host. Understanding the different levels of LRS regulation could lead to the development of improved nitrogen fixation efficiency to enhance sustainable agriculture, thereby reducing dependence on inorganic fertilizers.

Rhizobial tRNA-derived small RNAs are signal molecules regulating plant nodulation

Rhizobial infection and root nodule formation in legumes require recognition of signal molecules produced by the bacteria and their hosts. Here, we show that rhizobial transfer RNA (tRNA)-derived small RNA fragments (tRFs) are signal molecules that modulate host nodulation. Three families of rhizobial tRFs were confirmed to regulate host genes associated with nodule initiation and development through hijacking the host RNA-interference machinery that involves ARGONAUTE 1. Silencing individual tRFs with the use of short tandem target mimics or by overexpressing their targets represses root hair curling and nodule formation, whereas repressing these targets with artificial microRNAs identical to the respective tRFs or mutating these targets with CRISPR-Cas9 promotes nodulation. Our findings thus uncover a bacterial small RNA-mediated mechanism for prokaryote-eukaryote interaction and may pave the way for enhancing nodulation efficiency in legumes.

Genome-wide identification of AGO18b-bound miRNAs and phasiRNAs in maize by cRIP-seq

Background

Argonaute proteins (AGOs) are important players in the regulation of plant development by directing sRNAs to target mRNAs. In maize (Zea mays), AGO18b is a tassel-enriched and grass-specific AGO. Previous studies have shown that AGO18b is highly expressed in tassels during meiosis and negatively regulates determinacy of spikelet meristems. However, binding profile on RNAs and acting mechanisms of AGO18b remain unknown.

Results

In this study, we explored the binding profile of AGO18b in maize tassel by UV cross-linking RNA immunoprecipitation, followed by deep sequencing of these cDNA libraries (cRIP-seq), and systematically studied AGO18b-associated small RNAs and mRNAs by bioinformatics analysis. By globally analyzing the phased small-interfering RNA (phasiRNA) and miRNA abundance bound by AGO18b, we found AGO18b primarily binds to 21-nt phasiRNAs/miRNAs with a 5′-uridine and binds less strongly to 24-nt phasiRNAs with a 5′-adenosine in the premeiotic tassels. The abundance profile of AGO18b-associated miRNAs was different from their expression profile. Moreover, AGO18b strongly binds to miR166a-3p. We then obtained the AGO18b-bound mRNA targets of miR166a-3p by cRIP-seq, and confirmed the molecular function of AGO18b in regulating spikelet meristems.

Conclusions

Our results indicated that AGO18b binds to phasiRNAs with obvious 5 prime end bias under different sRNA length. MiRNAs and their target mRNAs associated with AGO18b indicated the molecular mechanisms of AGO18b as a negative regulator of inflorescence meristem and tassel development through integrating both phasiRNAs and miRNA pathways, which extended our view of sRNA regulation in flower development and provided potential methods to control pollination in the future.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6028-z) contains supplementary material, which is available to authorized users.

RNA Interference: A Natural Immune System of Plants to Counteract Biotic Stressors

During plant-pathogen interactions, plants have to defend the living transposable elements from pathogens. In response to such elements, plants activate a variety of defense mechanisms to counteract the aggressiveness of biotic stressors. RNA interference (RNAi) is a key biological process in plants to inhibit gene expression both transcriptionally and post-transcriptionally, using three different groups of proteins to resist the virulence of pathogens. However, pathogens trigger an anti-silencing mechanism through the expression of suppressors to block host RNAi. The disruption of the silencing mechanism is a virulence strategy of pathogens to promote infection in the invaded hosts. In this review, we summarize the RNA silencing pathway, anti-silencing suppressors, and counter-defenses of plants to viral, fungal, and bacterial pathogens.

Evolution and Diversification of Small RNA Pathways in Flowering Plants

Small regulatory RNAs guide gene silencing at the DNA or RNA level through repression of complementary sequences. The two main forms of small RNAs are microRNA (miRNA) and small interfering RNA (siRNAs), which are generated from the processing of different forms of double-stranded RNA (dsRNA) precursors. These two forms of small regulatory RNAs function in distinct but overlapping gene silencing pathways in plants. Gene silencing pathways in eukaryotes evolved from an ancient prokaryotic mechanism involved in genome defense against invasive genetic elements, but has since diversified to also play a crucial role in regulation of endogenous gene expression. Here, we review the biogenesis of the different forms of small RNAs in plants, including miRNAs, phased, secondary siRNAs (phasiRNAs) and heterochromatic siRNAs (hetsiRNAs), with a focus on their functions in genome defense, transcriptional and post-transcriptional gene silencing, RNA-directed DNA methylation, trans-chromosomal methylation and paramutation. We also discuss the important role that gene duplication has played in the functional diversification of gene silencing pathways in plants, and we highlight recently discovered components of gene silencing pathways in plants.

Elements Involved in the Rsv3-Mediated Extreme Resistance against an Avirulent Strain of Soybean Mosaic Virus

Extreme resistance (ER) is a type of R-gene-mediated resistance that rapidly induces a symptomless resistance phenotype, which is different from the phenotypical R-resistance manifested by the programmed cell death, accumulation of reactive oxygen species, and hypersensitive response. The Rsv3 gene in soybean cultivar L29 is responsible for ER against the avirulent strain G5H of soybean mosaic virus (SMV), but is ineffective against the virulent strain G7H. Rsv3-mediated ER is achieved through the rapid accumulation of callose, which arrests SMV-G5H at the point of infection. Callose accumulation, however, may not be the lone mechanism of this ER. Analyses of RNA-seq data obtained from infected soybean plants revealed a rapid induction of the abscisic acid pathway at 8 h post infection (hpi) in response to G5H but not to G7H, which resulted in the down-regulation of transcripts encoding β-1,3 glucanases that degrade callose in G5H-infected but not G7H-infected plants. In addition, parts of the autophagy and the small interfering (si) RNA pathways were temporally up-regulated at 24 hpi in response to G5H but not in response to G7H. The jasmonic acid (JA) pathway and many WRKY factors were clearly up-regulated only in G7H-infected plants. These results suggest that ER against SMV-G5H is achieved through the quick and temporary induction of ABA, autophagy, and the siRNA pathways, which rapidly eliminate G5H. The results also suggest that suppression of the JA pathway in the case of G5H is important for the Rsv3-mediated ER.

Dicer-like and RNA-dependent RNA polymerase gene family identification and annotation in the cultivated Solanum tuberosum and its wild relative S. commersonii

We provide advances in DCL and RDR gene diversity in Solanaceae. We also shed light on DCL and RDR gene expression in response to cold stress. DICER-like (DCL) and RNA-dependent RNA polymerase (RDR) genes form the core components to trigger small non-coding RNA (ncRNA) production. In spite of this, little is known about the two gene families in non-model plant species. As their genome sequences are now available, the cultivated potato (Solanum tuberosum) and its cold-tolerant wild relative Solanum commersonii offer a valuable opportunity to advance our understanding of the above genes. To determine the extent of diversification and evolution of DCLs and RDRs in these species, we performed a comparative analysis. Seven DCLs were identified in the two species, whereas seven and six RDR genes were found in S. tuberosum and S. commersonii, respectively. Based on phylogenetic analysis with DCLs and RDRs from several species, we provide evidence for an increase in their number in both potato species. We also disclosed that tandem duplications played a major role in the evolution of these gene families in Solanaceae. DCL and RDR expression was investigated in different tissues and under cold and virus stresses, with divergent profiles of the tandem duplicated genes being found in different tissues. DCL paralogs showed a contrasting expression in S. tuberosum and S. commersonii following cold stress and virus infection. By contrast, no change in RDR transcript activity was detected following both stresses. Overall, this study provides the first comparative genomic analysis of the core components of the RNAi machinery in Solanaceae and offers a scaffold for future functional analysis of these gene families.

Expression analysis of argonaute, Dicer-like, and RNA-dependent RNA polymerase genes in cucumber (Cucumis sativus L.) in response to abiotic stress

Posttranscriptional control of gene expression can be achieved through RNA interference when the activities of Dicer-like (DCL), argonaute (AGO) and RNA-dependent RNA polymerase (RDR) proteins are significant. In this study, we analysed the expression of seven AGO, five DCL and eight RDR genes in cucumber under cold, heat, hormone, salinity and dehydration treatments using quantitative reverse-transcription PCR (qRT-PCR). All CsAGO, CsDCL and CsRDR genes were differentially expressed under abiotic stress treatment. In response to abiotic stress treatment, most genes were expressed at higher levels in flowers or stems than in other organs, whereas some CsAGOs (CsAGO1c, CsAGO6 and CsAGO7) and CsRDRs (CsRDR1d and CsRDR2) were highly expressed in roots during dehydration treatment. The expression patterns indicate that most CsDCLs, CsAGOs and CsRDRs respond to abiotic stress, and stems or flowers are the most sensitive organs, followed by roots. This is the first report of expression analysis of all CsDCL, CsAGO and CsRDR family genes in cucumber under abiotic stress, which provides basic information and insights into the putative roles of these genes in abiotic stress. The results of this study should serve as a basis for further functional characterization of these gene families in cucumber and related Cucurbitaceae species.

A genome-wide analysis of the RNA-guided silencing pathway in coffee reveals insights into its regulatory mechanisms

microRNAs (miRNAs) are derived from self-complementary hairpin structures, while small-interfering RNAs (siRNAs) are derived from double-stranded RNA (dsRNA) or hairpin precursors. The core mechanism of sRNA production involves DICER-like (DCL) in processing the smallRNAs (sRNAs) and ARGONAUTE (AGO) as effectors of silencing, and siRNA biogenesis also involves action of RNA-Dependent RNA Polymerase (RDR), Pol IV and Pol V in biogenesis. Several other proteins interact with the core proteins to guide sRNA biogenesis, action, and turnover. We aimed to unravel the components and functions of the RNA-guided silencing pathway in a non-model plant species of worldwide economic relevance. The sRNA-guided silencing complex members have been identified in the Coffea canephora genome, and they have been characterized at the structural, functional, and evolutionary levels by computational analyses. Eleven AGO proteins, nine DCL proteins (which include a DCL1-like protein that was not previously annotated), and eight RDR proteins were identified. Another 48 proteins implicated in smallRNA (sRNA) pathways were also identified. Furthermore, we identified 235 miRNA precursors and 317 mature miRNAs from 113 MIR families, and we characterized ccp-MIR156, ccp-MIR172, and ccp-MIR390. Target prediction and gene ontology analyses of 2239 putative targets showed that significant pathways in coffee are targeted by miRNAs. We provide evidence of the expansion of the loci related to sRNA pathways, insights into the activities of these proteins by domain and catalytic site analyses, and gene expression analysis. The number of MIR loci and their targeted pathways highlight the importance of miRNAs in coffee. We identified several roles of sRNAs in C. canephora, which offers substantial insight into better understanding the transcriptional and post-transcriptional regulation of this major crop.

Mutations in Argonaute5 Illuminate Epistatic Interactions of the K1 and I Loci Leading to Saddle Seed Color Patterns in Glycine max

The soybean (Glycine max) seed coat has distinctive, genetically programmed patterns of pigmentation, and the recessive k1 mutation can epistatically overcome the dominant I and ii alleles, which inhibit seed color by producing small interfering RNAs (siRNAs) targeting chalcone synthase (CHS) mRNAs. Small RNA sequencing of dissected regions of immature seed coats demonstrated that CHS siRNA levels cause the patterns produced by the ii and ik alleles of the I locus, which restrict pigment to the hilum or saddle region of the seed coat, respectively. To identify the K1 locus, we compared RNA-seq data from dissected regions of two Clark isolines having similar saddle phenotypes mediated by CHS siRNAs but different genotypes (homozygous ik K1 versus homozygous ii k1). By examining differentially expressed genes, mapping information, and genome resequencing, we identified a 129-bp deletion in Glyma.11G190900 encoding Argonaute5 (AGO5), a member of the Argonaute family. Amplicon sequencing of several independent saddle pattern mutants from different genetic backgrounds revealed independent lesions affecting AGO5, thus establishing Glyma.11G190900 as the K1 locus. Nonfunctional AGO5 from k1 alleles leads to altered distributions of CHS siRNAs, thus explaining how the k1 mutation reverses the phenotype of the seed coat regions from yellow to pigmented, even in the presence of the normally dominant I or ii alleles.

Six Medicago truncatula Dicer-like protein genes are expressed in plant cells and upregulated in nodules

Here we report the existence of six putative Dicer-like genes in the Medicago truncatula genome. They are ubiquitously expressed throughout the plant and significantly induced in root nodules. Over the past decade, small noncoding RNAs (sncRNA) have emerged as widespread and important regulatory molecules influencing both the structure and expression of plant genomes. One of the key factors involved in sncRNA biogenesis in plants is a group of RNase III-type nucleases known as Dicer-like (DCL) proteins. Based on functional analysis of DCL proteins identified in Arabidopsis thaliana, four types of DCLs were distinguished (DCL1-4). DCL1 mainly produces 21 nt miRNAs. The products generated by DCL2, DCL3, and DCL4 belong to various classes of siRNAs that are 22, 24 and 21 nt in length, respectively. M. truncatula is a model legume plant closely related to many economically important cultivable species. By screening the recent M. truncatula genome assembly, we were able to identify three new DCL genes in addition to the MtDCL1-3 genes that had been earlier characterized. The newly found genes include MtDCL4 and two MtDCL2 homologs. We showed that all six M. truncatula DCL genes are expressed in plant cells. The first of the identified MtDCL2 paralogs encodes a truncated version of the DCL2 protein, while the second undergoes substantial and specific upregulation in the root nodules. Additionally, we identified an alternative splicing variant of MtDCL1 mRNA, similar to the one found in Arabidopsis. Our results indicate that DCL genes are differently activated during Medicago symbiosis with nitrogen fixing bacteria and upon pathogen infection. In addition, we hypothesize that the alternative splicing variant of MtDCL1 mRNA may be involved in tissue-specific regulation of the DCL1 level.

Evolution, functions, and mysteries of plant ARGONAUTE proteins

ARGONAUTE (AGO) proteins bind small RNAs (sRNAs) to form RNA-induced silencing complexes for transcriptional and post-transcriptional gene silencing. Genomes of primitive plants encode only a few AGO proteins. The Arabidopsis thaliana genome encodes ten AGO proteins, designated AGO1 to AGO10. Most early studies focused on these ten proteins and their interacting sRNAs. AGOs in other flowering plant species have duplicated and diverged from this set, presumably corresponding to new, diverged or specific functions. Among these, the grass-specific AGO18 family has been discovered and implicated as playing important roles during plant reproduction and viral defense. This review covers our current knowledge about functions and features of AGO proteins in both eudicots and monocots and compares their similarities and differences. On the basis of these features, we propose a new nomenclature for some plant AGOs.

New insights into pri-miRNA processing and accumulation in plants

MicroRNAs (miRNAs) regulate many biological processes such as development, metabolism, and others. They are processed from their primary transcripts called primary miRNA transcripts (pri-miRNAs) by the processor complex containing the RNAse III enzyme, DICER-LIKE1 (DCL1), in plants. Consequently, miRNA biogenesis is controlled through altering pri-miRNA accumulation and processing, which is crucial for plant development and adaptation to environmental changes. Plant pri-miRNAs are transcribed by DNA-dependent RNA polymerase II (Pol II) and their levels are determined through transcription and degradation, whereas pri-miRNA processing is affected by its structure, splicing, alternative splicing, loading to the processor and the processor activity, which involve in many accessory proteins. Here, we summarize recent progresses related to pri-miRNA transcription, stability, and processing in plants.

microRNA biogenesis, degradation and activity in plants

microRNAs (miRNAs) are important regulators of gene expression. After excised from primary miRNA transcript by dicer-like1 (DCL1, an RNAse III enzyme), miRNAs bind and guide their effector protein named argonaute 1 (AGO1) to silence the expression of target RNAs containing their complementary sequences in plants. miRNA levels and activities are tightly controlled to ensure their functions in various biological processes such as development, metabolism and responses to abiotic and biotic stresses. Studies have identified many factors that involve in miRNA accumulation and activities. Characterization of these factors in turn greatly improves our understanding of the processes related to miRNAs. Here, we review recent progress of mechanisms underlying miRNA expression and functions in plants.

De novo transcriptome sequence assembly from coconut leaves and seeds with a focus on factors involved in RNA-directed DNA methylation

Coconut palm (Cocos nucifera) is a symbol of the tropics and a source of numerous edible and nonedible products of economic value. Despite its nutritional and industrial significance, coconut remains under-represented in public repositories for genomic and transcriptomic data. We report de novo transcript assembly from RNA-seq data and analysis of gene expression in seed tissues (embryo and endosperm) and leaves of a dwarf coconut variety. Assembly of 10 GB sequencing data for each tissue resulted in 58,211 total unigenes in embryo, 61,152 in endosperm, and 33,446 in leaf. Within each unigene pool, 24,857 could be annotated in embryo, 29,731 could be annotated in endosperm, and 26,064 could be annotated in leaf. A KEGG analysis identified 138, 138, and 139 pathways, respectively, in transcriptomes of embryo, endosperm, and leaf tissues. Given the extraordinarily large size of coconut seeds and the importance of small RNA-mediated epigenetic regulation during seed development in model plants, we used homology searches to identify putative homologs of factors required for RNA-directed DNA methylation in coconut. The findings suggest that RNA-directed DNA methylation is important during coconut seed development, particularly in maturing endosperm. This dataset will expand the genomics resources available for coconut and provide a foundation for more detailed analyses that may assist molecular breeding strategies aimed at improving this major tropical crop.