Genome-wide identification, organization and phylogenetic analysis of Dicer-like, Argonaute and RNA-dependent RNA Polymerase gene families and their expression analysis during reproductive development and stress in rice

Genome-wide identification of gene families related to miRNA biogenesis in Mangifera indica L. and their possible role during heat stress

Mango is a popular tropical fruit that requires quarantine hot water treatment (QHWT) for postharvest sanitation, which can cause abiotic stress. Plants have various defense mechanisms to cope with stress; miRNAs mainly regulate the expression of these defense responses. Proteins involved in the biogenesis of miRNAs include DICER-like (DCL), ARGONAUTE (AGO), HYPONASTIC LEAVES 1 (HYL1), SERRATE (SE), HUA ENHANCER1 (HEN1), HASTY (HST), and HEAT-SHOCK PROTEIN 90 (HSP90), among others. According to our analysis, the mango genome contains five DCL, thirteen AGO, six HYL, two SE, one HEN1, one HST, and five putative HSP90 genes. Gene structure prediction and domain identification indicate that sequences contain key domains for their respective gene families, including the RNase III domain in DCL and PAZ and PIWI domains for AGOs. In addition, phylogenetic analysis indicates the formation of clades that include the mango sequences and their respective orthologs in other flowering plant species, supporting the idea these are functional orthologs. The analysis of cis-regulatory elements of these genes allowed the identification of MYB, ABRE, GARE, MYC, and MeJA-responsive elements involved in stress responses. Gene expression analysis showed that most genes are induced between 3 to 6 h after QHWT, supporting the early role of miRNAs in stress response. Interestingly, our results suggest that mango rapidly induces the production of miRNAs after heat stress. This research will enable us to investigate further the regulation of gene expression and its effects on commercially cultivated fruits, such as mango, while maintaining sanitary standards.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Barley AGO4 proteins show overlapping functionality with distinct small RNA-binding properties in heterologous complementation

KEY MESSAGE

Barley AGO4 proteins complement expressional changes of epigenetically regulated genes in Arabidopsis ago4-3 mutant and show a distinct affinity for the 5' terminal nucleotide of small RNAs, demonstrating functional conservation and divergence. The function of Argonaute 4 (AGO4) in Arabidopsis thaliana has been extensively characterized; however, its role in monocots, which have large genomes abundantly supplemented with transposable elements (TEs), remains elusive. The study of barley AGO4 proteins can provide insights into the conserved aspects of RNA-directed DNA methylation (RdDM) and could also have further applications in the field of epigenetics or crop improvement. Bioinformatic analysis of RNA sequencing data identified two active AGO4 genes in barley, HvAGO4a and HvAGO4b. These genes function similar to AtAGO4 in an Arabidopsis heterologous complementation system, primarily binding to 24-nucleotide long small RNAs (sRNAs) and triggering methylation at specific target loci. Like AtAGO4, HvAGO4B exhibits a preference for binding sRNAs with 5' adenine residue, while also accepting 5' guanine, uracil, and cytosine residues. In contrast, HvAGO4A selectively binds only sRNAs with a 5' adenine residue. The diverse binding capacity of barley AGO4 proteins is reflected in TE-derived sRNAs and in their varying abundance. Both barley AGO4 proteins effectively restore the levels of extrachromosomal DNA and transcript abundancy of the heat-activated ONSEN retrotransposon to those observed in wild-type Arabidopsis plants. Our study provides insight into the distinct binding specificities and involvement in TE regulation of barley AGO4 proteins in Arabidopsis by heterologous complementation.

Genome-Wide Identification and Characterization of Major RNAi Genes Highlighting Their Associated Factors in Cowpea (Vigna unguiculata (L.) Walp.)

In different regions of the world, cowpea (Vigna unguiculata (L.) Walp.) is an important vegetable and an excellent source of protein. It lessens the malnutrition of the underprivileged in developing nations and has some positive effects on health, such as a reduction in the prevalence of cancer and cardiovascular disease. However, occasionally, certain biotic and abiotic stresses caused a sharp fall in cowpea yield. Major RNA interference (RNAi) genes like Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are essential for the synthesis of their associated factors like domain, small RNAs (sRNAs), transcription factors, micro-RNAs, and cis-acting factors that shield plants from biotic and abiotic stresses. In this study, applying BLASTP search and phylogenetic tree analysis with reference to the Arabidopsis RNAi (AtRNAi) genes, we discovered 28 VuRNAi genes, including 7 VuDCL, 14 VuAGO, and 7 VuRDR genes in cowpea. We looked at the domains, motifs, gene structures, chromosomal locations, subcellular locations, gene ontology (GO) terms, and regulatory factors (transcription factors, micro-RNAs, and cis-acting elements (CAEs)) to characterize the VuRNAi genes and proteins in cowpea in response to stresses. Predicted VuDCL1, VuDCL2(a, b), VuAGO7, VuAGO10, and VuRDR6 genes might have an impact on cowpea growth, development of the vegetative and flowering stages, and antiviral defense. The VuRNAi gene regulatory features miR395 and miR396 might contribute to grain quality improvement, immunity boosting, and pathogen infection resistance under salinity and drought conditions. Predicted CAEs from the VuRNAi genes might play a role in plant growth and development, improving grain quality and production and protecting plants from biotic and abiotic stresses. Therefore, our study provides crucial information about the functional roles of VuRNAi genes and their associated components, which would aid in the development of future cowpeas that are more resilient to biotic and abiotic stress. The manuscript is available as a preprint at this link: doi:10.1101/2023.02.15.528631v1.

Genomic assembly, characterization, and quantification of DICER-like gene family in Okra plants under dehydration conditions

Background

Okra is a plant farmed for its pods, leaves, and stems all of which are edible. It is famous for its ability to tolerate long desiccation periods. It belongs to the Malvaceae family and is a sister species to hibiscus, cotton, and cacao plants.

Methods

In the current study, okra plants were used as a model to sequence, assemble, and analyze the evolutionary and functional characteristics of the Dicer-like protein gene family (DCL) based on DNAseq and qPCR techniques.

Results

Four Dicer-like (DCL) single-copy genes of the okra plant Abelmoschus esculentus (L.) Moench (AeDCL) were successfully assembled. The lengths of the AeDCL copies were 8,494, 5,214, 4,731, and 9,329 bp. The detected exons in these samples ranged from a single exon in AeDCL3 to 24 exons in AeDCL4. AeDCLs had five functional domains of two DEAD-like helicase superfamilies, N and C; one Dicer domain; one ribonuclease III domain (a and b); and one double-stranded RNA-binding domain. The PAZ domain was completely annotated only for AeDCL1 and AeDCL3. All AeDCLs were up-regulated under drought conditions, with leaves showing more extensive fold changes than roots. The study focused on a comprehensive genome-wide identification and analysis of the DCL gene family in naturally drought-tolerant okra plants, an orphan crop that can be used as a model for further genomic and transcriptomic studies on drought-tolerance mechanisms in plants.

DNA-dependent RNA polymerases in plants

DNA-dependent RNA polymerases (Pols) transfer the genetic information stored in genomic DNA to RNA in all organisms. In eukaryotes, the typical products of nuclear Pol I, Pol II, and Pol III are ribosomal RNAs, mRNAs, and transfer RNAs, respectively. Intriguingly, plants possess two additional Pols, Pol IV and Pol V, which produce small RNAs and long noncoding RNAs, respectively, mainly for silencing transposable elements. The five plant Pols share some subunits, but their distinct functions stem from unique subunits that interact with specific regulatory factors in their transcription cycles. Here, we summarize recent advances in our understanding of plant nucleus-localized Pols, including their evolution, function, structures, and transcription cycles.

ZmAGO18b negatively regulates maize resistance against southern leaf blight

KEY MESSAGE

Here, we report that ZmAGO18b encoding an argonaute protein is a negative regulator of maize resistance against southern leaf blight. Southern leaf blight caused by fungal pathogen Cochliobolus heterostrophus is a destructive disease on maize throughout the world. Argonaute (AGO) proteins, key regulators in small RNA pathway, play important roles in plant defense. But whether they have function in maize resistance against C. heterostrophus is unknown. Association analysis between the nucleic variation of 18 ZmAGO loci with disease phenotype against C. heterostrophus was performed, and the ZmAGO18b locus was identified to be associated with resistance against C. heterostrophus. Overexpression of ZmAGO18b gene suppresses maize resistance against C. heterostrophus, and mutation of ZmAGO18b enhances maize resistance against C. heterostrophus. Further, we identified the resistant haplotype of ZmAGO18b by association analysis of natural variation in ZmAGO18b genomic DNA sequences with seedling resistance phenotypes against C. heterostrophus and confirmed the resistant haplotype is co-segregated with resistance phenotypes against C. heterostrophus in two F₂ populations. In sum, this study reports that ZmAGO18b negatively regulates maize resistance against C. heterostrophus.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory element analyses in sunflower (Helianthus annuus)

RNA interference (RNAi) regulates a variety of eukaryotic gene expressions that are engaged in response to stress, growth, and the conservation of genomic stability during developmental phases. It is also intimately connected to the post-transcriptional gene silencing (PTGS) process and chromatin modification levels. The entire process of RNA interference (RNAi) pathway gene families mediates RNA silencing. The main factors of RNA silencing are the Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) gene families. To the best of our knowledge, genome-wide identification of RNAi gene families like DCL, AGO, and RDR in sunflower (Helianthus annuus) has not yet been studied despite being discovered in some species. So, the goal of this study is to find the RNAi gene families like DCL, AGO, and RDR in sunflower based on bioinformatics approaches. Therefore, we accomplished an inclusive in silico investigation for genome-wide identification of RNAi pathway gene families DCL, AGO, and RDR through bioinformatics approaches such as (sequence homogeneity, phylogenetic relationship, gene structure, chromosomal localization, PPIs, GO, sub-cellular localization). In this study, we have identified five DCL (HaDCLs), fifteen AGO (HaAGOs), and ten RDR (HaRDRs) in the sunflower genome database corresponding to the RNAi genes of model plant Arabidopsis thaliana based on genome-wide analysis and a phylogenetic method. The analysis of the gene structure that contains exon-intron numbers, conserved domain, and motif composition analyses for all HaDCL, HaAGO, and HaRDR gene families indicated almost homogeneity among the same gene family. The protein-protein interaction (PPI) network analysis illustrated that there exists interconnection among identified three gene families. The analysis of the Gene Ontology (GO) enrichment showed that the detected genes directly contribute to the RNA gene-silencing and were involved in crucial pathways. It was observed that the cis-acting regulatory components connected to the identified genes were shown to be responsive to hormone, light, stress, and other functions. That was found in HaDCL, HaAGO, and HaRDR genes associated with the development and growth of plants. Finally, we are able to provide some essential information about the components of sunflower RNA silencing through our genome-wide comparison and integrated bioinformatics analysis, which open the door for further research into the functional mechanisms of the identified genes and their regulatory elements.

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.

Mobile ARGONAUTE 1d binds 22-nt miRNAs to generate phasiRNAs important for low-temperature male fertility in rice

Phased small interfering RNAs (phasiRNAs) are abundantly expressed in anthers and linked to environment-related male fertility in grasses, yet how they function under different environmental conditions remains unclear. Here, we identified a rice (Oryza sativa) low temperature-induced Argonaute (AGO) protein, OsAGO1d, that is responsible for generating phasiRNAs and preserving male fertility at low temperature. Loss of OsAGO1d function causes low-temperature male sterility associated with delayed programmed cell death of tapetal cells during anther development. OsAGO1d binds miR2118 and miR2275 family members and triggers phasiRNA biogenesis; it also binds 21-nt phasiRNAs with a 5' terminal U. In total, phasiRNAs from 972 loci are OsAGO1d-dependent. OsAGO1d protein moves from anther wall cells into meiocytes, where it loads miR2275 to produce 24-nt phasiRNAs. Together, our results show that OsAGO1d acts as a mobile signal to fine-tune phasiRNA production and this function is important for male fertility at low temperature.

Genome-Wide Identification and Posttranscriptional Regulation Analyses Elucidate Roles of Key Argonautes and Their miRNA Triggers in Regulating Complex Yield Traits in Rapeseed

Argonautes (AGOs) interact with microRNAs (miRNAs) to form the RNA-induced silencing complex (RISC), which can posttranscriptionally regulate the expression of targeted genes. To date, however, the AGOs and their miRNA triggers remain elusive in rapeseed (Brassica napus). Here, we systematically performed a phylogenetic analysis and examined the collinear relationships of the AGOs among four Brassicaceae species. Their physicochemical properties, gene structures, and expression patterns among 81 tissues from multiple materials and developmental stages were further analyzed. Additionally, their posttranscriptional regulation was analyzed using psRNATarget prediction, miRNA-/mRNA-Seq analyses, and a qRT-PCR verification. We finally identified 10 AtAGOs, 13 BolAGOs, 11 BraAGOs, and 24 BnaAGOs. An expression analysis of the BnaAGOs in the B. napus cultivar ZS11, as well as genotypes with extreme phenotypes in various yield-related traits, revealed the conservation and diversity of these genes. Furthermore, we speculated the posttranscriptional regulation of the B. napus miR168a-AGO1s and miR403-AGO2s modules. Combining miRNA-Seq and mRNA-Seq analyses, we found that the B. napus miR168a-AGO1s module may play an essential role in negatively regulating yield traits, whereas the miR403-AGO2s module positively impacts yield. This is the first attempt to comprehensively analyze the AGOs and their miRNA triggers in B. napus and provides a theoretical basis for breeding high-yielding varieties through the manipulation of the miRNA-AGOs modules.

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

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

Research progress on plant noncoding RNAs in response to low-temperature stress

Low temperature (LT) is an important factor limiting plant growth and distribution. Plants have evolved sophisticated adaptive mechanisms to cope with hypothermia. RNA silencing is the orchestrator of these cellular responses. RNA silencing, which modifies gene expression through noncoding RNAs (ncRNAs), is a strategy used by plants to combat environmental stress. ncRNAs, which have very little protein-coding capacity, work by binding reverse complementary endogenous transcripts. In plants, ncRNAs include small non-coding RNAs (sncRNAs), medium-sized non-coding RNAs (mncRNAs), and long non-coding RNAs (lncRNAs). Apart from describing the biogenesis of different ncRNAs (miRNAs, siRNAs, and lncRNAs), we thoroughly discuss the functions of these ncRNAs during cold acclimation. Two major classes of sncRNAs, microRNAs and siRNAs, play essential regulatory roles in cold response processes through the posttranscriptional gene silencing (PTGS) pathway or transcriptional gene silencing (TGS) pathway. Microarray or transcriptome sequencing analysis can reveal a large number of cold-responsive miRNAs in plants. In this review, the cold-response patterns of miRNAs verified by Northern blotting or quantitative PCR in Arabidopsis thaliana, rice, and many other important crops are discussed. The detailed molecular mechanisms of several miRNAs in Arabidopsis (miR397, miR408, miR402, and miR394) and rice (Osa-miR156, Osa-miR319, and Osa-miR528) that regulate plant cold resistance are elucidated. In addition, the regulatory mechanism of the lncRNA SVALKA in the cold signaling pathway is explained in detail. Finally, we present the challenges for understanding the roles of small ncRNAs in cold signal transduction.

Argonaute5 and its associated small RNAs modulate the transcriptional response during the rhizobia-Phaseolus vulgaris symbiosis

Both plant- and rhizobia-derived small RNAs play an essential role in regulating the root nodule symbiosis in legumes. Small RNAs, in association with Argonaute proteins, tune the expression of genes participating in nodule development and rhizobial infection. However, the role of Argonaute proteins in this symbiosis has been overlooked. In this study, we provide transcriptional evidence showing that Argonaute5 (AGO5) is a determinant genetic component in the root nodule symbiosis in Phaseolus vulgaris. A spatio-temporal transcriptional analysis revealed that the promoter of PvAGO5 is active in lateral root primordia, root hairs from rhizobia-inoculated roots, nodule primordia, and mature nodules. Transcriptional analysis by RNA sequencing revealed that gene silencing of PvAGO5 affected the expression of genes involved in the biosynthesis of the cell wall and phytohormones participating in the rhizobial infection process and nodule development. PvAGO5 immunoprecipitation coupled to small RNA sequencing revealed the small RNAs bound to PvAGO5 during the root nodule symbiosis. Identification of small RNAs associated to PvAGO5 revealed miRNAs previously known to participate in this symbiotic process, further supporting a role for AGO5 in this process. Overall, the data presented shed light on the roles that PvAGO5 plays during the root nodule symbiosis in P. vulgaris.

Comprehensive In Silico Analysis of RNA Silencing-Related Genes and Their Regulatory Elements in Wheat (Triticum aestivum L.)

Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are known as the three major gene families that act as the critical components of RNA interference or silencing mechanisms through the noncoding small RNA molecules (miRNA and siRNA) to regulate the expressions of protein-coding genes in eukaryotic organisms. However, most of their characteristics including structures, chromosomal location, subcellular locations, regulatory elements, and gene networking were not rigorously studied. Our analysis identified 7 TaDCL, 39 TaAGO, and 16 TaRDR genes as RNA interference (RNAi) genes from the wheat genome. Phylogenetic analysis of predicted RNAi proteins with the RNAi proteins of Arabidopsis and rice showed that the predicted proteins of TaDCL, TaAGO, and TaRDR groups are clustered into four, eight, and four subgroups, respectively. Domain, 3D protein structure, motif, and exon-intron structure analyses showed that these proteins conserve identical characteristics within groups and maintain differences between groups. The nonsynonymous/synonymous mutation ratio (Ka/Ks) < 1 suggested that these protein sequences conserve some purifying functions. RNAi genes networking with TFs revealed that ERF, MIKC-MADS, C2H2, BBR-BPC, MYB, and Dof are the key transcriptional regulators of the predicted RNAi-related genes. The cis-regulatory element (CREs) analysis detected some important CREs of RNAi genes that are significantly associated with light, stress, and hormone responses. Expression analysis based on an online database exhibited that almost all of the predicted RNAi genes are expressed in different tissues and organs. A case-control study from the gene expression level showed that some RNAi genes significantly responded to the drought and heat stresses. Overall results would therefore provide an excellent basis for in-depth molecular investigation of these genes and their regulatory elements for wheat crop improvement against different stressors.

Genome-Wide Identification and Expression Analyses of the Cotton AGO Genes and Their Potential Roles in Fiber Development and Stress Response

Argonaute proteins (AGOs) are indispensable components of RNA silencing. However, systematic characterization of the AGO genes have not been completed in cotton until now. In this study, cotton AGO genes were identified and analyzed with respect to their evolution and expression profile during biotic and abiotic stresses. We identified 14 GaAGO, 14 GrAGO, and 28 GhAGO genes in the genomes of Gossypium arboreum, Gossypium raimondii, and Gossypium hirsutum. Cotton AGO proteins were classified into four subgroups. Structural and functional conservation were observed in the same subgroups based on the analysis of the gene structure and conserved domains. Twenty-four duplicated gene pairs were identified in GhAGO genes, and all of them exhibited strong purifying selection during evolution. Moreover, RNA-seq analysis showed that most of the GhAGO genes exhibit high expression levels in the fiber initiation and elongation processes. Furthermore, the expression profiles of GhAGO genes tested by quantitative real-time polymerase chain reaction (qPCR) demonstrated that they were sensitive to Verticillium wilt infection and salt and drought stresses. Overall, our results will pave the way for further functional investigation of the cotton AGO gene family, which may be involved in fiber development and stress response.

Genome-wide H3K9 acetylation level increases with age-dependent senescence of flag leaf in rice

Flag leaf senescence is an important biological process that drives the remobilization of nutrients to the growing organs of rice. Leaf senescence is controlled by genetic information via gene expression and histone modification, but the precise mechanism is as yet unclear. Here, we analysed genome-wide acetylated lysine residue 9 of histone H3 (H3K9ac) enrichment by chromatin immunoprecipitation-sequencing (ChIP-seq), and examined its association with transcriptomes by RNA-seq during flag leaf aging in rice (Oryza sativa). We found that genome-wide H3K9 acetylation levels increased with age-dependent senescence in rice flag leaf, and there was a positive correlation between the density and breadth of H3K9ac with gene expression and transcript elongation. During flag leaf aging, we observed 1249 up-regulated differentially expressed genes (DEGs) and 996 down-regulated DEGs, showing a strong relationship between temporal changes in gene expression and gain/loss of H3K9ac. We produced a landscape of H3K9 acetylation-modified gene expression targets that include known senescence-associated genes, metabolism-related genes, as well as miRNA biosynthesis-related genes. Our findings reveal a complex regulatory network of metabolism- and senescence-related pathways mediated by H3K9ac, and elucidate patterns of H3K9ac-mediated regulation of gene expression during flag leaf aging in rice.

Integrative Analysis of the DICER-like (DCL) Genes From Peach (Prunus persica): A Critical Role in Response to Drought Stress

DICER-likes (DCLs) proteins are the core component for non-coding RNA (ncRNA) biogenesis, playing essential roles in some biological processes. The DCL family has been characterized in model plants, such as Arabidopsis, rice, and poplar. However, the evolutionary aspect and the expression mechanism under drought stress were scarce and have never been reported and characterized in one of the most important worldwide cultivated fruit trees, peach (Prunus persica). Eight DCLs genes in the Prunus persica genome were detected, in addition to 51 DCLs in the other seven Rosaceae genomes. The phylogenetic analysis with Arabidopsis thaliana and RTL1 gene as outgroups suggested that DCL members are divided into four clades: DCL1, DCL2, DCL3, and DCL4 with several gene gain/loss events of DCL gene copies through the evolutionary tract of the Rosacea family. The number of homologous DCL copies within each clade, along with the chromosomal location indicated gene duplication event of the DCL2 gene occurred once for the subfamily Amygdaloideae and twice for Pyrus communis and Prunus dulics and trice for the P. persica on Chromosome number 7 genes. Another duplication event was found for the DCL3 gene that occurred once for all the eight Rosaceae species with no match in A. thaliana. The DCL genetic similarity and activity was evaluated using BLASTp and previously published RNA-seq data among different tissues and over different time points of peach trees exposed to drought conditions. Finally, the expression pattern of PrupeDCLs in response to drought stress was identified, and two of these members, Prupe.7G047900 and Prupe.6G363600, were found as main candidate genes for response to drought stress. Our data presented here provide useful information for a better understanding of the molecular evolution of DCL genes in Rosaceae genomes, and the function of DCLs in P. persica.

Proteomic and Phosphoproteomic Analyses Reveal a Complex Network Regulating Pollen Abortion and Potential Candidate Proteins in TCMS Wheat

Thermosensitive sterile lines are natural materials for exploring the effects of anther development on male fertility. To study the possible molecular mechanisms regulating protein activity during the induction of male sterility, proteomic and phosphoproteomic analyses with tandem mass tags (TMTs) were used to study the binucleate anther of the thermosensitive sterile wheat line YS3038. A total of 9072 proteins, including 5019 phosphoproteins, were identified. Enrichment analyses of differentially abundant proteins (DAPs) and phosphoproteins (DAPPs) in metabolic pathways showed that both were mainly related to energy metabolism. Soluble sugar and ATP content were significantly decreased, free fatty acid content was significantly increased, and ROS was abnormally accumulated in male sterile YS3038-A. In addition, 233 kinase–substrate pairs involved in potential phosphorylation control networks were predicted to regulate fertility. Candidate proteins were identified, and a quantitative real-time polymerase chain reaction (qRT-PCR) analysis was used to validate the TMT results. TaPDCD5 is likely to be involved in fertility conversion of YS3038 by barley stripe mosaic virus-induced gene silencing (BSMV-VIGS). Our data provide new insights into the mechanism of TCMS, which has value for identifying potential candidate proteins associated with the formation or abortion of pollen and promotion of wheat heterosis utilization.

KNOX II transcription factor HOS59 functions in regulating rice grain size

Plant Knotted1-like homeobox (KNOX) genes encode homeodomain-containing transcription factors. In rice (Oryza sativa L.), little is known about the downstream target genes of KNOX Class II subfamily proteins. Here we generated chromatin immunoprecipitation (ChIP)-sequencing datasets for HOS59, a member of the rice KNOX Class II subfamily, and characterized the genome-wide binding sites of HOS59. We conducted trait ontology (TO) analysis of 9705 identified downstream target genes, and found that multiple TO terms are related to plant structure morphology and stress traits. ChIP-quantitative PCR (qPCR) was conducted to validate some key target genes. Meanwhile, our IP-MS datasets showed that HOS59 was closely associated with BELL family proteins, some grain size regulators (OsSPL13, OsSPL16, OsSPL18, SLG, etc.), and some epigenetic modification factors such as OsAGO4α and OsAGO4β, proteins involved in small interfering RNA-mediated gene silencing. Furthermore, we employed CRISPR/Cas9 editing and transgenic approaches to generate hos59 mutants and overexpression lines, respectively. Compared with wild-type plants, the hos59 mutants have longer grains and increased glume cell length, a loose plant architecture, and drooping leaves, while the overexpression lines showed smaller grain size, erect leaves, and lower plant height. The qRT-PCR results showed that mutation of the HOS59 gene led to upregulation of some grain size-related genes such as OsSPL13, OsSPL18, and PGL2. In summary, our results indicate that HOS59 may be a repressor of the downstream target genes, negatively regulating glume cell length, rice grain size, plant architecture, etc. The identified downstream target genes and possible interaction proteins of HOS59 improve our understanding of the KNOX regulatory networks.

Silencing of Dicer-like protein 2a restores the resistance phenotype in the rice mutant, sxi4 (suppressor of Xa21-mediated immunity 4)

The rice immune receptor XA21 confers resistance to Xanthomonas oryzae pv. oryzae (Xoo), and upon recognition of the RaxX21-sY peptide produced by Xoo, XA21 activates the plant immune response. Here we screened 21 000 mutant plants expressing XA21 to identify components involved in this response, and reported here the identification of a rice mutant, sxi4, which is susceptible to Xoo. The sxi4 mutant carries a 32-kb translocation from chromosome 3 onto chromosome 7 and displays an elevated level of DCL2a transcript, encoding a Dicer-like protein. Silencing of DCL2a in the sxi4 genetic background restores resistance to Xoo. RaxX21-sY peptide-treated leaves of sxi4 retain the hallmarks of XA21-mediated immune response. However, WRKY45-1, a known negative regulator of rice resistance to Xoo, is induced in the sxi4 mutant in response to RaxX21-sY peptide treatment. A CRISPR knockout of a short interfering RNA (TE-siRNA815) in the intron of WRKY45-1 restores the resistance phenotype in sxi4. These results suggest a model where DCL2a accumulation negatively regulates XA21-mediated immunity by altering the processing of TE-siRNA815.

Origin, evolution and diversification of plant ARGONAUTE proteins

Argonaute (AGO) proteins are central players in RNA interference in eukaryotes. They associate with small RNAs (sRNA) and lead to transcriptional or posttranscriptional silencing of targets, thereby regulating diverse biological processes. The molecular and biological functions of AGO proteins have been extensively characterized, particularly in a few angiosperm species, leading to the recognition that the AGO family has expanded to accommodate diverse sRNAs thereby performing diverse biological functions. However, understanding of the expansion of AGO proteins in plants is still limited, due to a dearth of knowledge of AGO proteins in green algal groups. Here, we identified more than 2900 AGO proteins from 244 plant species, including green algae, and performed a large-scale phylogenetic analysis. The phylogeny shows that the plant AGO family gave rise to four clades after the emergence of hydrobiontic algae and prior to the emergence of land plants. Subsequent parallel expansion in ferns and angiosperms resulted in eight main clades in angiosperms: AGO2, AGO7, AGO6, AGO4, AGO1, AGO10a, AGO10b and AGO5. On the basis of this phylogeny, we identified two novel AGO4 orthologs that Arabidopsis does not have, and redefined AGO10, which is composed of AGO10a and AGO10b. Finally, we propose a hypothetical evolutionary model of AGO proteins in plants. Our studies provide a deeper understanding of the phylogenetic relationships of AGO family members in the green lineage, which would help to further reveal their roles as RNAi effectors.

RNAi-Based Antiviral Innate Immunity in Plants

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

Big Role of Small RNAs in Female Gametophyte Development

In living organisms, sexual reproduction relies on the successful development of the gametes. Flowering plants produce gametes in the specialized organs of the flower, the gametophytes. The female gametophyte (FG), a multicellular structure containing female gametes (egg cell and central cell), is often referred to as an embryo sac. Intriguingly, several protein complexes, molecular and genetic mechanisms participate and tightly regulate the female gametophyte development. Recent evidence indicates that small RNA (sRNA) mediated pathways play vital roles in female gametophyte development and specification. Here, we present an insight into our understanding and the recent updates on the molecular mechanism of different players of small RNA-directed regulatory pathways during ovule formation and growth.

CRISPR-Cas9 mediated OsMIR168a knockout reveals its pleiotropy in rice

MicroRNA168 (MIR168) is a key miRNA that targets the main RNA-induced silencing complex component Argonaute 1 (AGO1) to regulate plant growth and environmental stress responses. However, the regulatory functions of MIR168 need to be further elucidated in rice. In this paper, we generated clean OsMIR168a deletion mutants by CRISPR-Cas9 strategy. We then phenotypically and molecularly characterized these mutants. The rice OsMIR168a mutants grew rapidly at the seedling stage, produced more tillers and matured early. Compared to the wild-type plants, the mutants were shorter at maturity and produced smaller spikelets and seeds. Analysis of gene expression showed that the transcription levels of OsMIR168a's target genes such as OsAGO1a, OsAGO1b and OsAGO1d were elevated significantly in the OsMIR168a mutants. Intriguingly, OsAGO18, a member of a new AGO clade that is conserved in monocots, was confirmed to be a target of OsMIR168a not only by informatic prediction but also by expression analysis and a cell-based cleavage assay in the OsMIR168a mutants. Many protein-coding genes and miRNAs showed differential expression in the OsMIR168a mutants, suggesting OsMIR168a exerts a major transcriptional regulatory role, likely through its potential target genes such as OsAGO1s and OsAGO18. KEGG enrichment analysis of these differentially expressed genes pointed to OsMIR168a's involvement in important processes such as plant hormone signalling transduction and plant-pathogen interaction. These data collectively support that the complex regulation module of OsMIR168a-OsAGO1/OsAGO18-miRNAs-target genes contributes to agronomically important traits, which sheds light on miRNA-mediated crop breeding.

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.

An atypical class of non-coding small RNAs is produced in rice leaves upon bacterial infection

Non-coding small RNAs (sRNA) act as mediators of gene silencing and regulate plant growth, development and stress responses. Early insights into plant sRNAs established a role in antiviral defense and they are now extensively studied across plant-microbe interactions. Here, sRNA sequencing discovered a class of sRNA in rice (Oryza sativa) specifically associated with foliar diseases caused by Xanthomonas oryzae bacteria. Xanthomonas-induced small RNAs (xisRNAs) loci were distinctively upregulated in response to diverse virulent strains at an early stage of infection producing a single duplex of 20-22 nt sRNAs. xisRNAs production was dependent on the Type III secretion system, a major bacterial virulence factor for host colonization. xisRNA loci overlap with annotated transcripts sequences, with about half of them encoding protein kinase domain proteins. A number of the corresponding rice cis-genes have documented functions in immune signaling and xisRNA loci predominantly coincide with the coding sequence of a conserved kinase motif. xisRNAs exhibit features of small interfering RNAs and their biosynthesis depend on canonical components OsDCL1 and OsHEN1. xisRNA induction possibly mediates post-transcriptional gene silencing but they do not broadly suppress cis-genes expression on the basis of mRNA-seq data. Overall, our results identify a group of unusual sRNAs with a potential role in plant-microbe interactions.

Use of deep sequencing to profile small RNAs derived from tomato spotted wilt orthotospovirus and hippeastrum chlorotic ringspot orthotospovirus in infected Capsicum annuum

Virus-derived small RNAs are one of the key factors of RNA silencing in plant defence against viruses. We obtained virus-derived small interfering RNA profiles from Tomato spotted wilt orthotospovirus and Hippeastrum chlorotic ringspot orthotospovirus infected Capsicum annuum XX19 and XY11 by deep sequencing one day after inoculation. The vsiRNAs data were mapped to the TSWV and HCRV genomes, and the results showed that the vsiRNAs measured 19-24 nucleotides in length. Most of the vsiRNAs were mapped to the S segment of the viral genome. For XX19 and XY11 infected with HCRV, the distribution range of vsiRNAs in S RNA was 52.06-55.20%, while for XX19 and XY11 infected with TSWV, the distribution range of vsiRNAs in S RNA was 87.76-89.07%. The first base at the 5' end of the siRNA from TSWV and HCRV was primarily biased towards A, U, or C. Compared with mock-inoculated XX19 and XY11, the expression level of CaRDR1 was upregulated in TSWV- and HCRV-inoculated XX19 and XY11. CaAGO2 and CaAGO5 were upregulated in XY11 against HCRV infection, and CaRDR2 was downregulated in TSWV-infected XY11 and XX19. The profile of HCRV and TSWV vsiRNA verified in this study could be useful for selecting key vsiRNA such as those in disease-resistant varieties by artificially synthesizing amiRNA.

OsbZIP62/OsFD7, a functional ortholog of FLOWERING LOCUS D, regulates floral transition and panicle development in rice

We have characterized a rice bZIP protein-coding gene OsbZIP62/OsFD7 that is expressed preferentially in the shoot apical meristem and during early panicle developmental stages in comparison with other OsFD genes characterized to date. Surprisingly, unlike OsFD1, OsFD7 interacts directly and more efficiently with OsFTLs; the interaction is strongest with OsFTL1 followed by Hd3a and RFT1, as confirmed by fluorescence lifetime imaging-Förster resonant energy transfer (FLIM-FRET) analysis. In addition, OsFD7 is phosphorylated at its C-terminal end by OsCDPK41 and OsCDPK49 in vitro, and this phosphorylated moiety is recognized by OsGF14 proteins. OsFD7 RNAi transgenics were late flowering; the transcript levels of some floral meristem identity genes (e.g. OsMADS14, OsMADS15, and OsMADS18) were also down-regulated. RNAi lines also exhibited dense panicle morphology with an increase in the number of primary and secondary branches resulting in longer panicles and more seeds, probably due to down-regulation of SEPALLATA family genes. In comparison with other FD-like proteins previously characterized in rice, it appears that OsFD7 may have undergone diversification during evolution, resulting in the acquisition of newer functions and thus playing a dual role in floral transition and panicle development in rice.

Spatiotemporal regulation and roles of reproductive phasiRNAs in plants

Since co-suppression was discovered as a pioneer silencing phenomenon of RNA interference (RNAi) in petunia in 1990, many types of small RNAs have been identified in the RNAi pathway among various eukaryotes. In plants, a large number of 21- or 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs) are produced via processing of long RNA precursors by Dicer-like proteins. However, the roles of phasiRNAs remain largely unknown. The development of imaging technology and RNA profiling has clarified the spatiotemporal regulation of phasiRNAs, and subsequently the different functions of 21-nt trans-acting phasiRNAs and 24-nt cis-regulatory phasiRNAs during male organ development. This review focuses on the biogenesis, diversification, spatiotemporal expression pattern and function of phasiRNAs in plants.

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

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

Beyond identity: Understanding the contribution of the 5' nucleotide of the antisense strand to RNAi activity

In both the pharmaceutical and agricultural fields, RNA-based products have capitalized upon the mechanism of RNA interference for targeted reduction of gene expression to improve phenotypes and traits. Reduction in gene expression by RNAi is the result of a small interfering RNA (siRNA) molecule binding to an ARGONAUTE (AGO) protein and directing the effector complex to a homologous region of a target gene's mRNA. siRNAs properties that govern RNA-AGO association have been studied in detail. The siRNA 5' nucleotide (nt) identity has been demonstrated in plants to be an important property responsible for directing association of endogenous small RNAs with different AGO effector proteins. However, it has not been investigated whether the 5' nt identity is an efficacious determinant for topically-applied chemically synthesized siRNAs. In this study, we employed a sandpaper abrasion method to study the silencing efficacies of topically-applied 21 base-pair siRNA duplexes. The MAGNESIUM CHELATASE and GREEN FLUORESCENT PROTEIN genes were selected as endogenous and transgenic gene targets, respectively, to assess the molecular and phenotypic effects of gene silencing. Collections of siRNA variants with different 5' nt identities and different pairing states between the 5' antisense nt and its match in the sense strand of the siRNA duplex were tested for their silencing efficacy. Our results suggest a flexibility in the 5' nt requirement for topically applied siRNA duplexes in planta and highlight the similarity of 5' thermodynamic rules governing topical siRNA efficacy across plants and animals.

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

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

Ubiquitin-dependent Argonauteprotein MEL1 degradation is essential for rice sporogenesis and phasiRNA target regulation

MEIOSIS ARRESTED AT LEPTOTENE1 (MEL1), a rice (Oryza sativa) Argonaute (AGO) protein, has been reported to function specifically at premeiotic and meiotic stages of germ cell development and is associated with a novel class of germ cell-specific small noncoding RNAs called phased small RNAs (phasiRNAs). MEL1 accumulation is temporally and spatially regulated and is eliminated after meiosis. However, the metabolism and turnover (i.e. the homeostasis) of MEL1 during germ cell development remains unknown. Here, we show that MEL1 is ubiquitinated and subsequently degraded via the proteasome pathway in vivo during late sporogenesis. Abnormal accumulation of MEL1 after meiosis leads to a semi-sterile phenotype. We identified a monocot-specific E3 ligase, XBOS36, a CULLIN RING-box protein, that is responsible for the degradation of MEL1. Ubiquitination at four K residues at the N terminus of MEL1 by XBOS36 induces its degradation. Importantly, inhibition of MEL1 degradation either by XBOS36 knockdown or by MEL1 overexpression prevents the formation of pollen at the microspore stage. Further mechanistic analysis showed that disrupting MEL1 homeostasis in germ cells leads to off-target cleavage of phasiRNA target genes. Our findings thus provide insight into the communication between a monocot-specific E3 ligase and an AGO protein during plant reproductive development.

Genome-Wide Identification and Evolutionary Analysis of Argonaute Genes in Hexaploid Bread Wheat

Argonaute (AGO) proteins play a pivotal role in plant growth and development as the core components of RNA-induced silencing complex (RISC). However, no systematic characterization of AGO genes in wheat has been reported to date. In this study, a total number of 69 TaAGO genes in the hexaploid bread wheat (Triticum aestivum cv. Chinese Spring) genome, divided into 10 subfamilies, were identified. Compared to all wheat genes, TaAGOs showed a significantly lower evolutionary rate, which is consistent with their high conservation in eukaryotes. However, the homoeolog retention was remarkably higher than the average, implying the nonredundant biological importance of TaAGO genes in bread wheat. Further homoeologous gene expression bias analyses revealed that TaAGOs may have undergone neofunctionalization after polyploidization and duplication through the divergent expression of homoeologous gene copies, to provide new opportunities for the generation of adaptive traits. Moreover, quantitative real-time polymerase chain reaction (qRT-PCR) analyses indicated that TaAGO gene expression was involved in response to heat, drought, and salt stresses. Our results would provide a theoretical basis for future studies on the biological functions of TaAGO genes in wheat and other gramineous species.

StCaM2, a calcium binding protein, alleviates negative effects of salinity and drought stress in tobacco

KEY MESSAGE

Overexpression of StCaM2 in tobacco promotes plant growth and confers increased salinity and drought tolerance by enhancing the photosynthetic efficiency, ROS scavenging, and recovery from membrane injury. Calmodulins (CaMs) are important Ca2+ sensors that interact with effector proteins and drive a network of signal transduction pathways involved in regulating the growth and developmental pattern of plants under stress. Herein, using in silico analysis, we identified 17 CaM isoforms (StCaM) in potato. Expression profiling revealed different temporal and spatial expression patterns of these genes, which were modulated under abiotic stress. Among the identified StCaM genes, StCaM2 was found to have the largest number of abiotic stress responsive promoter elements. In addition, StCaM2 was upregulated in response to some of the selected abiotic stress in potato tissues. Overexpression of StCaM2 in transgenic tobacco plants enhanced their tolerance to salinity and drought stress. Accumulation of reactive oxygen species was remarkably decreased in transgenic lines compared to that in wild type plants. Chlorophyll a fluorescence analysis suggested better performance of photosystem II in transgenic plants under stress compared to that in wild type plants. The increase in salinity stress tolerance in StCaM2-overexpressing plants was also associated with a favorable K+/Na+ ratio. The enhanced tolerance to abiotic stresses correlated with the increase in the activities of anti-oxidative enzymes in transgenic tobacco plants. Overall, our results suggest that StCaM2 can be a novel candidate for conferring salt and drought tolerance in plants.

Genome-wide identification, evolutionary relationship and expression analysis of AGO, DCL and RDR family genes in tea

Three gene families in plants viz. Argonaute (AGOs), Dicer-like (DCLs) and RNA dependent RNA polymerase (RDRs) constitute the core components of small RNA mediated gene silencing machinery. The present study endeavours to identify members of these gene families in tea and to investigate their expression patterns in different tissues and various stress regimes. Using genome-wide analysis, we have identified 18 AGOs, 5 DCLs and 9 RDRs in tea, and analyzed their phylogenetic relationship with orthologs of Arabidopsis thaliana. Gene expression analysis revealed constitutive expression of CsAGO1 in all the studied tissues and stress conditions, whereas CsAGO10c showed most variable expression among all the genes. CsAGO10c gene was found to be upregulated in tissues undergoing high meristematic activity such as buds and roots, as well as in Exobasidium vexans infected samples. CsRDR2 and two paralogs of CsAGO4, which are known to participate in biogenesis of hc-siRNAs, showed similarities in their expression levels in most of the tea plant tissues. This report provides first ever insight into the important gene families involved in biogenesis of small RNAs in tea. The comprehensive knowledge of these small RNA biogenesis purveyors can be utilized for tea crop improvement aimed at stress tolerance and quality enhancement.

Core RNA Interference Genes Involved in miRNA and Ta-siRNA Biogenesis in Hops and Their Expression Analysis after Challenging with Verticillium nonalfalfae

RNA interference is an evolutionary conserved mechanism by which organisms regulate the expression of genes in a sequence-specific manner to modulate defense responses against various abiotic or biotic stresses. Hops are grown for their use in brewing and, in recent years, for the pharmaceutical industry. Hop production is threatened by many phytopathogens, of which Verticillium, the causal agent of Verticillium wilt, is a major contributor to yield losses. In the present study, we performed identification, characterization, phylogenetic, and expression analyses of three Argonaute, two Dicer-like, and two RNA-dependent RNA polymerase genes in the susceptible hop cultivar Celeia and the resistant cultivar Wye Target after infection with Verticillium nonalfalfae. Phylogeny results showed clustering of hop RNAi proteins with their orthologues from the closely related species Cannabis sativa, Morus notabilis and Ziziphus jujuba which form a common cluster with species of the Rosaceae family. Expression analysis revealed downregulation of argonaute 2 in both cultivars on the third day post-inoculation, which may result in reduced AGO2-siRNA-mediated posttranscriptional gene silencing. Both cultivars may also repress ta-siRNA biogenesis at different dpi, as we observed downregulation of argonaute 7 in the susceptible cultivar on day 1 and downregulation of RDR6 in the resistant cultivar on day 3 after inoculation.

Integrated microRNA and transcriptome profiling reveal key miRNA-mRNA interaction pairs associated with seed development in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat seed development is an extremely complex process involving many gene regulatory pathways. MicroRNAs (miRNAs) have been identified as the important negative regulators of gene expression and performed crucial regulatory roles in various plant biological processes. However, whether miRNAs participate in Tartary buckwheat seed development remains unexplored.

RESULTS

In this study, we first identified 26 miRNA biosynthesis genes in the Tartary buckwheat genome and described their phylogeny and expression profiling. Then we performed small RNA (sRNA) sequencing for Tartary buckwheat seeds at three developmental stages to identify the miRNAs associated with seed development. In total, 230 miRNAs, including 101 conserved and 129 novel miRNAs, were first identified in Tartary buckwheat, and 3268 target genes were successfully predicted. Among these miRNAs, 76 exhibited differential expression during seed development, and 1534 target genes which correspond to 74 differentially expressed miRNAs (DEMs) were identified. Based on integrated analysis of DEMs and their targets expression, 65 miRNA-mRNA interaction pairs (25 DEMs corresponding to 65 target genes) were identified that exhibited significantly opposite expression during Tartary buckwheat seed development, and 6 of the miRNA-mRNA pairs were further verified by quantitative real-time polymerase chain reaction (qRT-PCR) and ligase-mediated rapid amplification of 5' cDNA ends (5'-RLM-RACE). Functional annotation of the 65 target mRNAs showed that 56 miRNA-mRNA interaction pairs major involved in cell differentiation and proliferation, cell elongation, hormones response, organogenesis, embryo and endosperm development, seed size, mineral elements transport, and flavonoid biosynthesis, which indicated that they are the key miRNA-mRNA pairs for Tartary buckwheat seed development.

CONCLUSIONS

Our findings provided insights for the first time into miRNA-mediated regulatory pathways in Tartary buckwheat seed development and suggested that miRNAs play important role in Tartary buckwheat seed development. These findings will be help to study the roles and regulatory mechanism of miRNAs in Tartary buckwheat seed development.

Contribution of Small RNA Pathway to Interactions of Rice with Pathogens and Insect Pests

Small RNAs (sRNAs) are mainly classified into microRNAs (miRNAs) and small interfering RNAs (siRNAs) according to their origin. miRNAs originate from single-stranded RNA precursors, whereas siRNAs originate from double-stranded RNA precursors that are synthesized by RNA-dependent RNA polymerases. Both of single-stranded and double-stranded RNA precursors are processed into sRNAs by Dicer-like proteins. Then, the sRNAs are loaded into ARGONAUTE proteins, forming RNA-induced silencing complexes (RISCs). The RISCs repress the expression of target genes with sequences complementary to the sRNAs through the cleavage of transcripts, the inhibition of translation or DNA methylation. Here, we summarize the recent progress of sRNA pathway in the interactions of rice with various parasitic organisms, including fungi, viruses, bacteria, as well as insects. Besides, we also discuss the hormone signal in sRNA pathway, and the emerging roles of circular RNAs and long non-coding RNAs in rice immunity. Obviously, small RNA pathway may act as a part of rice innate immunity to coordinate with growth and development.

Rice-specific Argonaute 17 controls reproductive growth and yield-associated phenotypes

This manuscript describes the functions of an Argonaute protein named AGO17 in rice. AGO17 is required for the development of rice reproductive tissues. Argonaute (AGO) proteins are a well-conserved multigene family of regulators mediating gene silencing across eukaryotes. Monocot plants have additional members of AGO, the functions of which are poorly understood. Among the non-dicot AGO1 clade members in monocots, AGO17 expresses highly in reproductive tissues. Here we show that overexpression of Oryza sativa indica AGO17 in rice resulted in robust growth and increased yield, whereas its silencing resulted in reduced panicle length, less fertility, and poor growth. Small (s)RNA transcriptome analysis revealed misregulation of several miRNAs and other categories of sRNAs in silenced and overexpression lines, in agreement with its likely competition with other AGO1 clade members. Targets of differentially expressed miRNAs included previously unreported target RNAs coding for proteins involved in development, phase transition, and transport. Our results indicate a distinctive role for OsAGO17 in rice reproductive development that could be harnessed to improve yield.

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.

Genome-wide identification of Argonautes in Solanaceae with emphasis on potato

Regulatory small RNAs (sRNAs) play important roles in many fundamental processes in plant biology such as development, fertilization and stress responses. The AGO protein family has here a central importance in gene regulation based on their capacity to associate with sRNAs followed by mRNA targeting in a sequence-complementary manner. The present study explored Argonautes (AGOs) in the Solanaceae family, with emphasis on potato, Solanum tuberosum (St). A genome-wide monitoring was performed to provide a deeper insight into gene families, genomic localization, gene structure and expression profile against the potato late blight pathogen Phytophthora infestans. Among 15 species in the Solanaceae family we found a variation from ten AGOs in Nicotiana obtusifolia to 17 in N. tabacum. Comprehensive analyses of AGO phylogeny revealed duplication of AGO1, AGO10 and AGO4 paralogs during early radiation of Solanaceae. Fourteen AGOs were identified in potato. Orthologs of AGO8 and AGO9 were missing in the potato genome. However, AGO15 earlier annotated in tomato was identified. StAGO15 differs from the other paralogs having residues of different physico-chemical properties at functionally important amino acid positions. Upon pathogen challenge StAGO15 was significantly activated and hence may play a prominent role in sRNA-based regulation of potato defense.

Premeiotic, 24-Nucleotide Reproductive PhasiRNAs Are Abundant in Anthers of Wheat and Barley But Not Rice and Maize

Two classes of premeiotic (21-nucleotides [nt]) and meiotic (24-nt) phased small interfering RNAs (phasiRNAs) and their patterns of accumulation have been described in maize (Zea mays) and rice (Oryza sativa) anthers. Their precise function remains unclear, but studies have shown that they support male fertility. The important role of phasiRNAs in anthers underpins our present study to characterize these small RNAs in wheat (Triticum aestivum) and barley (Hordeum vulgare) anthers. We staged anthers at every 0.2 mm of development for one wheat and two barley varieties. We isolated premeiotic (0.2, 0.4, and 0.6 mm), meiotic (0.8, 1.0, and 1.4 mm), and postmeiotic (1.8 mm) anthers, for which we then investigated accumulation patterns of RNAs, including reproductive phasiRNAs. We annotated a total of 12,821 and 2,897 PHAS loci in the wheat and barley genomes, respectively. By comparing the total number of PHAS loci in genomes of maize, rice, barley, and wheat, we identified an expansion of reproductive PHAS loci in the genomes of Poaceae subfamilies from Panicoideae to Oryzoideae and to Poideae. In addition to the two classes of premeiotic (21-nt) and meiotic (24-nt) phasiRNAs, previously described in maize and rice anthers, we characterized a group of 24-nt phasiRNAs that accumulate in premeiotic anthers. The absence of premeiotic 24-nt phasiRNAs in maize and rice suggests a divergence in grass species of the Poideae subfamily. Additionally, we performed a gene coexpression analysis describing the regulation of phasiRNA biogenesis in wheat and barley anthers. We highlight Argonaute 9 (AGO9) and Argonaute 6 (AGO6) as candidate binding partners of premeiotic and meiotic 24-nt phasiRNAs, respectively.

Premeiotic, 24-Nucleotide Reproductive PhasiRNAs Are Abundant in Anthers of Wheat and Barley But Not Rice and Maize

Two classes of premeiotic (21-nucleotides [nt]) and meiotic (24-nt) phased small interfering RNAs (phasiRNAs) and their patterns of accumulation have been described in maize (Zea mays) and rice (Oryza sativa) anthers. Their precise function remains unclear, but studies have shown that they support male fertility. The important role of phasiRNAs in anthers underpins our present study to characterize these small RNAs in wheat (Triticum aestivum) and barley (Hordeum vulgare) anthers. We staged anthers at every 0.2 mm of development for one wheat and two barley varieties. We isolated premeiotic (0.2, 0.4, and 0.6 mm), meiotic (0.8, 1.0, and 1.4 mm), and postmeiotic (1.8 mm) anthers, for which we then investigated accumulation patterns of RNAs, including reproductive phasiRNAs. We annotated a total of 12,821 and 2,897 PHAS loci in the wheat and barley genomes, respectively. By comparing the total number of PHAS loci in genomes of maize, rice, barley, and wheat, we identified an expansion of reproductive PHAS loci in the genomes of Poaceae subfamilies from Panicoideae to Oryzoideae and to Poideae. In addition to the two classes of premeiotic (21-nt) and meiotic (24-nt) phasiRNAs, previously described in maize and rice anthers, we characterized a group of 24-nt phasiRNAs that accumulate in premeiotic anthers. The absence of premeiotic 24-nt phasiRNAs in maize and rice suggests a divergence in grass species of the Poideae subfamily. Additionally, we performed a gene coexpression analysis describing the regulation of phasiRNA biogenesis in wheat and barley anthers. We highlight Argonaute 9 (AGO9) and Argonaute 6 (AGO6) as candidate binding partners of premeiotic and meiotic 24-nt phasiRNAs, respectively.

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.

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

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