MicroRNAs and other small RNAs enriched in the Arabidopsis RNA-dependent RNA polymerase-2 mutant

Transposon-Associated Small RNAs Involved in Plant Defense in Poplar

Utilizing high-throughput Illumina sequencing, we examined how small RNA (sRNA) profiles vary in Chinese white poplar (Populus tomentosa) across two pivotal infection stages by the rust fungus Melampsora larici-populina: the biotrophic growth phase (T02; 48 h post infection) and the urediniospore development and dispersal phase (T03; 168 h), both essential for plant colonization and prolonged biotrophic engagement. Far exceeding random expectations, siRNA clusters predominantly arose from transposon regions, with pseudogenes also contributing significantly, and infection-stage-specific variations were notably tied to these transposon-derived siRNAs. As the infection advanced, clusters of 24 nt siRNAs in transposon and intergenic regions exhibited pronounced abundance shifts. An analysis of targets indicated that Populus sRNAs potentially regulate 95% of Melampsora larici-populina genes, with pathogen effector genes showing heightened targeting by sRNAs during the biotrophic and urediniospore phases compared to controls, pointing to selective sRNA-target interactions. In contrast to conserved miRNAs across plant species, Populus-specific miRNAs displayed a markedly greater tendency to target NB-LRR genes. These observations collectively highlight the innovative roles of sRNAs in plant defense, their evolutionary roots, and their dynamic interplay with pathogen coevolution.

Transposon-triggered epigenetic chromatin dynamics modulate EFR-related pathogen response

Infectious diseases drive wild plant evolution and impact crop yield. Plants, like animals, sense biotic threats through pattern recognition receptors (PRRs). Overly robust immune responses can harm plants; thus, understanding the tuning of defense response mechanisms is crucial for developing pathogen-resistant crops. In this study, we found that an inverted-repeat transposon (EFR-associated IR, Ea-IR) located between the loci encoding PRRs ELONGATION FACTOR-TU RECEPTOR (EFR) and myosin XI-k (XI-k) in Arabidopsis affects chromatin organization, promoting the formation of a repressive chromatin loop. Upon pathogen infection, chromatin changes around EFR and XI-k correlate with increased EFR transcription. Pathogen-induced chromatin opening causes RNA polymerase II readthrough, producing a longer, Ea-IR-containing XI-k transcript, processed by Dicer-like enzymes into small RNAs, which reset chromatin to a repressive state attenuating the immune response after infection. Arabidopsis accessions lacking Ea-IR have higher basal EFR levels and resistance to pathogens. We show a scenario in which a transposon, chromatin organization and gene expression interact to fine-tune immune responses, during both the course of infection and the course of evolution.

Regulation of plant immunity via small RNA-mediated control of NLR expression

Plants use different receptors to detect potential pathogens: membrane-anchored pattern recognition receptors (PRRs) activated upon perception of pathogen-associated molecular patterns (PAMPs) that elicit pattern-triggered immunity (PTI); and intracellular nucleotide-binding leucine-rich repeat proteins (NLRs) activated by detection of pathogen-derived effectors, activating effector-triggered immunity (ETI). The interconnections between PTI and ETI responses have been increasingly reported. Elevated NLR levels may cause autoimmunity, with symptoms ranging from fitness cost to developmental arrest, sometimes combined with run-away cell death, making accurate control of NLR dosage key for plant survival. Small RNA-mediated gene regulation has emerged as a major mechanism of control of NLR dosage. Twenty-two nucleotide miRNAs with the unique ability to trigger secondary siRNA production from target transcripts are particularly prevalent in NLR regulation. They enhance repression of the primary NLR target, but also bring about repression of NLRs only complementary to secondary siRNAs. We summarize current knowledge on miRNAs and siRNAs in the regulation of NLR expression with an emphasis on 22 nt miRNAs and propose that miRNA and siRNA regulation of NLR levels provides additional links between PTI and NLR defense pathways to increase plant responsiveness against a broad spectrum of pathogens and control an efficient deployment of defenses.

Plant polymerase IV sensitizes chromatin through histone modifications to preclude spread of silencing into protein-coding domains

Across eukaryotes, gene regulation is manifested via chromatin states roughly distinguished as heterochromatin and euchromatin. The establishment, maintenance, and modulation of the chromatin states is mediated using several factors including chromatin modifiers. However, factors that avoid the intrusion of silencing signals into protein-coding genes are poorly understood. Here we show that a plant specific paralog of RNA polymerase (Pol) II, named Pol IV, is involved in avoidance of facultative heterochromatic marks in protein-coding genes, in addition to its well-established functions in silencing repeats and transposons. In its absence, H3K27 trimethylation (me3) mark intruded the protein-coding genes, more profoundly in genes embedded with repeats. In a subset of genes, spurious transcriptional activity resulted in small(s) RNA production, leading to post-transcriptional gene silencing. We show that such effects are significantly pronounced in rice, a plant with a larger genome with distributed heterochromatin compared with Arabidopsis Our results indicate the division of labor among plant-specific polymerases, not just in establishing effective silencing via sRNAs and DNA methylation but also in influencing chromatin boundaries.

Polymorphic inverted repeats near coding genes impact chromatin topology and phenotypic traits in Arabidopsis thaliana

Transposons are mobile elements that are commonly silenced to protect eukaryotic genome integrity. In plants, transposable element (TE)-derived inverted repeats (IRs) are commonly found near genes, where they affect host gene expression. However, the molecular mechanisms of such regulation are unclear in most cases. Expression of these IRs is associated with production of 24-nt small RNAs, methylation of the IRs, and drastic changes in local 3D chromatin organization. Notably, many of these IRs differ between Arabidopsis thaliana accessions, causing variation in short-range chromatin interactions and gene expression. CRISPR-Cas9-mediated disruption of two IRs leads to a switch in genome topology and gene expression with phenotypic consequences. Our data show that insertion of an IR near a gene provides an anchor point for chromatin interactions that profoundly impact the activity of neighboring loci. This turns IRs into powerful evolutionary agents that can contribute to rapid adaptation.

Roles of microRNAs in abiotic stress response and characteristics regulation of plant

MicroRNAs (miRNAs) are a class of non-coding endogenous small RNAs (long 20-24 nucleotides) that negatively regulate eukaryotes gene expression at post-transcriptional level via cleavage or/and translational inhibition of targeting mRNA. Based on the diverse roles of miRNA in regulating eukaryotes gene expression, research on the identification of miRNA target genes has been carried out, and a growing body of research has demonstrated that miRNAs act on target genes and are involved in various biological functions of plants. It has an important influence on plant growth and development, morphogenesis, and stress response. Recent case studies indicate that miRNA-mediated regulation pattern may improve agronomic properties and confer abiotic stress resistance of plants, so as to ensure sustainable agricultural production. In this regard, we focus on the recent updates on miRNAs and their targets involved in responding to abiotic stress including low temperature, high temperature, drought, soil salinity, and heavy metals, as well as plant-growing development. In particular, this review highlights the diverse functions of miRNAs on achieving the desirable agronomic traits in important crops. Herein, the main research strategies of miRNAs involved in abiotic stress resistance and crop traits improvement were summarized. Furthermore, the miRNA-related challenges and future perspectives of plants have been discussed. miRNA-based research lays the foundation for exploring miRNA regulatory mechanism, which aims to provide insights into a potential form of crop improvement and stress resistance breeding.

Simple Webserver-Facilitated Method to Design and Synthesize Artificial miRNA Gene and Its Application in Engineering Viral Resistance

Plant viruses impose serious threats on crop production. Artificial miRNAs can mediate specific and effective gene silencing in plants and are widely used in plant gene function studies and to engineer plant viral resistance. To facilitate the design of artificial miRNA genes, we developed a webserver, AMIRdesigner, which can be used to design oligos for artificial miRNA synthesis using wild-type and permutated MIR171 and MIR164 backbones. The artificial miRNA genes designed by AMIRdesigner can be easily assembled into miRNA clusters for multiple target sites. To validate the server functionality, we designed four artificial miRNA genes targeting four conserved regions in the potato leafroll virus genome using AMIRdesigner. These genes were synthesized with the server-designed oligos and further assembled into a quadruple miRNA cluster, which was cloned into an overexpression vector and transformed into potato plants. Small RNA Northern blot and virus inoculation analyses showed that a high level of artificial miRNA expression and good viral resistance were achieved in some of the transgenic lines. These results demonstrate the utility of our webserver AMIRdesigner for engineering crop viral resistance.

miR778 mediates gene expression, histone modification, and DNA methylation during cyst nematode parasitism

Despite the known critical regulatory functions of microRNAs, histone modifications, and DNA methylation in reprograming plant epigenomes in response to pathogen infection, the molecular mechanisms underlying the tight coordination of these components remain poorly understood. Here, we show how Arabidopsis (Arabidopsis thaliana) miR778 coordinately modulates the root transcriptome, histone methylation, and DNA methylation via post-transcriptional regulation of the H3K9 methyltransferases SU(var)3-9 homolog 5 (SUVH5) and SUVH6 upon infection by the beet cyst nematode Heterodera schachtii. miR778 post-transcriptionally silences SUVH5 and SUVH6 upon nematode infection. Manipulation of the expression of miR778 and its two target genes significantly altered plant susceptibility to H. schachtii. RNA-seq analysis revealed a key role of SUVH5 and SUVH6 in reprograming the transcriptome of Arabidopsis roots upon H. schachtii infection. In addition, chromatin immunoprecipitation (ChIP)-seq analysis established SUVH5 and SUVH6 as the main enzymes mediating H3K9me2 deposition in Arabidopsis roots in response to nematode infection. ChIP-seq analysis also showed that these methyltransferases possess distinct DNA binding preferences in that they are targeting transposable elements under noninfected conditions and protein-coding genes in infected plants. Further analyses indicated that H3K9me2 deposition directed by SUVH5 and SUVH6 contributes to gene expression changes both in roots and in nematode feeding sites and preferentially associates with CG DNA methylation. Together, our results uncovered multi-layered epigenetic regulatory mechanisms coordinated by miR778 during Arabidopsis-H. schachtii interactions.

MicroRNA-mediated host defense mechanisms against pathogens and herbivores in rice: balancing gains from genetic resistance with trade-offs to productivity potential

BACKGROUND

Rice (Oryza sativa L.) is the major source of daily caloric intake for more than 30% of the human population. However, the sustained productivity of this staple food crop is continuously threatened by various pathogens and herbivores. Breeding has been successful in utilizing various mechanisms of defense by gene pyramiding in elite cultivars, but the continuous resurgence of highly resistant races of pathogens and herbivores often overcomes the inherent capacity of host plant immunity. MicroRNAs (miRNAs) are endogenous, short, single-stranded, non-coding RNA molecules that regulate gene expression by sequence-specific cleavage of target mRNA or suppressing target mRNA translation. While miRNAs function as upstream regulators of plant growth, development, and host immunity, their direct effects on growth and development in the context of balancing defenses with agronomic potential have not been extensively discussed and explored as a more viable strategy in breeding for disease and pest resistant cultivars of rice with optimal agronomic potentials.

RESULTS

Using the available knowledge in rice and other model plants, this review examines the important roles of miRNAs in regulating host responses to various fungal, bacterial, and viral pathogens, and insect pests, in the context of gains and trade-offs to crop yield. Gains from R-gene-mediated resistance deployed in modern rice cultivars are often undermined by the rapid breakdown of resistance, negative pleiotropic effects, and linkage drags with undesirable traits. In stark contrast, several classes of miRNAs are known to efficiently balance the positive gains from host immunity without significant costs in terms of losses in agronomic potentials (i.e., yield penalty) in rice. Defense-related miRNAs such as Osa-miR156, Osa-miR159, Osa-miR162, Osa-miR396, Osa-530, Osa-miR1432, Osa-miR1871, and Osa-miR1873 are critical in fine-tuning and integrating immune responses with physiological processes that are necessary to the maintenance of grain yield. Recent research has shown that many defense-related miRNAs regulate complex and agronomically important traits.

CONCLUSIONS

Identification of novel immune-responsive miRNAs that orchestrate physiological processes critical to the full expression of agronomic potential will facilitate the stacking of optimal combinations of miRNA-encoding genes to develop high-yielding cultivars with durable resistance to disease and insect pests with minimal penalties to yield.

RNA Interference: Promising Approach to Combat Plant Viruses

Plant viruses are devastating plant pathogens that severely affect crop yield and quality. Plants have developed multiple lines of defense systems to combat viral infection. Gene silencing/RNA interference is the key defense system in plants that inhibits the virulence and multiplication of pathogens. The general mechanism of RNAi involves (i) the transcription and cleavage of dsRNA into small RNA molecules, such as microRNA (miRNA), or small interfering RNA (siRNA), (ii) the loading of siRNA/miRNA into an RNA Induced Silencing Complex (RISC), (iii) complementary base pairing between siRNA/miRNA with a targeted gene, and (iv) the cleavage or repression of a target gene with an Argonaute (AGO) protein. This natural RNAi pathway could introduce transgenes targeting various viral genes to induce gene silencing. Different RNAi pathways are reported for the artificial silencing of viral genes. These include Host-Induced Gene Silencing (HIGS), Virus-Induced Gene Silencing (VIGS), and Spray-Induced Gene Silencing (SIGS). There are significant limitations in HIGS and VIGS technology, such as lengthy and time-consuming processes, off-target effects, and public concerns regarding genetically modified (GM) transgenic plants. Here, we provide in-depth knowledge regarding SIGS, which efficiently provides RNAi resistance development against targeted genes without the need for GM transgenic plants. We give an overview of the defense system of plants against viral infection, including a detailed mechanism of RNAi, small RNA molecules and their types, and various kinds of RNAi pathways. This review will describe how RNA interference provides the antiviral defense, recent improvements, and their limitations.

A unique miR775-GALT9 module regulates leaf senescence in Arabidopsis during post-submergence recovery by modulating ethylene and the abscisic acid pathway

The submergence-induced hypoxic condition negatively affects the plant growth and development, and causes early onset of senescence. Hypoxia alters the expression of a number of microRNAs (miRNAs). However, the molecular function of submergence stress-induced miRNAs in physiological or developmental changes and recovery remains poorly understood. Here, we show that miR775 is an Arabidopsis thaliana-specific young and unique miRNA that possibly evolved non-canonically. miR775 post-transcriptionally regulates GALACTOSYLTRANSFERASE 9 (GALT9) and their expression is inversely affected at 24 h of complete submergence stress. The overexpression of miR775 (miR775-Oe) confers enhanced recovery from submergence stress and reduced accumulation of RBOHD and ROS, in contrast to wild-type and MIM775 Arabidopsis shoot. A similar recovery phenotype in the galt9 mutant indicates the role of the miR775-GALT9 module in post-submergence recovery. We predicted that Golgi-localized GALT9 is potentially involved in protein glycosylation. The altered expression of senescence-associated genes (SAG12, SAG29 and ORE1), ethylene signalling (EIN2 and EIN3) and abscisic acid (ABA) biosynthesis (NCED3) pathway genes occurs in miR775-Oe, galt9 and MIM775 plants. Thus, our results indicate the role for the miR775-GALT9 module in post-submergence recovery through a crosstalk between the ethylene signalling and ABA biosynthesis pathways.

Reinforcement of CHH methylation through RNA-directed DNA methylation ensures sexual reproduction in rice

DNA methylation is an important epigenetic mark that regulates the expression of genes and transposons. RNA-directed DNA methylation (RdDM) is the main molecular pathway responsible for de novo DNA methylation in plants. Although the mechanism of RdDM has been well studied in Arabidopsis (Arabidopsis thaliana), most mutations in RdDM genes cause no remarkable developmental defects in Arabidopsis. Here, we isolated and cloned Five Elements Mountain 1 (FEM1), which encodes RNA-dependent RNA polymerase 2 (OsRDR2) in rice (Oryza sativa). Mutation in OsRDR2 abolished the accumulation of 24-nt small interfering RNAs, and consequently substantially decreased genome-wide CHH (H = A, C, or T) methylation. Moreover, male and female reproductive development was disturbed, which led to sterility in osrdr2 mutants. We discovered that OsRDR2-dependent DNA methylation may regulate the expression of multiple key genes involved in stamen development, meiosis, and pollen viability. In wild-type (WT) plants but not in osrdr2 mutants, genome-wide CHH methylation levels were greater in panicles, stamens, and pistils than in seedlings. The global increase of CHH methylation in reproductive organs of the WT was mainly explained by the enhancement of RdDM activity, which includes OsRDR2 activity. Our results, which revealed a global increase in CHH methylation through enhancement of RdDM activity in reproductive organs, suggest a crucial role for OsRDR2 in the sexual reproduction of rice.

Methodologies for Discovery and Quantitative Profiling of sRNAs in Potato

Small RNAs (sRNAs) are short noncoding RNAs involved in the regulation of a wide range of biological processes in plants. Advances in high-throughput sequencing and development of new computational tools had facilitated the discovery of different classes of sRNAs, their quantification, and elucidation of their functional role in gene expression regulation by target transcript predictions. The workflow presented here allows identification of different sRNA species: known and novel potato miRNAs, and their sequence variants (isomiRs), as well as identification of phased small interfering RNAs (phasiRNAs). Moreover, it includes steps for differential expression analysis to search for regulated sRNAs across different tested biological conditions. In addition, it describes two different methods for predicting sRNA targets, in silico prediction, and degradome sequencing data analysis. All steps of the workflow are written in a clear and user-friendly way; thus they can be followed also by the users with minimal bioinformatics knowledge. We also included several in-house scripts together with valuable notes to facilitate data (pre)processing steps and to reduce the analysis time.

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.

Combined Transcriptome and Proteome Analysis of Maize (Zea mays L.) Reveals A Complementary Profile in Response to Phosphate Deficiency

A deficiency in the macronutrient phosphate (Pi) brings about various changes in plants at the morphological, physiological and molecular levels. However, the molecular mechanism for regulating Pi homeostasis in response to low-Pi remains poorly understood, particularly in maize (Zea mays L.), which is a staple crop and requires massive amounts of Pi. Therefore, in this study, we performed expression profiling of the shoots and roots of maize seedlings with Pi-tolerant genotype at both the transcriptomic and proteomic levels using RNA sequencing and isobaric tags for relative and absolute quantitation (iTRAQ). We identified 1944 differentially expressed transcripts and 340 differentially expressed proteins under low-Pi conditions. Most of the differentially expressed genes were clustered as regulators, such as transcription factors involved in the Pi signaling pathway at the transcript level. However, the more functional and metabolism-related genes showed expression changes at the protein level. Moreover, under low-Pi conditions, Pi transporters and phosphatases were specifically induced in the roots at both the transcript and protein levels, and increased amounts of mRNA and protein of two purple acid phosphatases (PAPs) and one UDP-sulfoquinovose synthase (SQD) were specifically detected in the roots. The new insights provided by this study will help to improve the P-utilization efficiency of maize.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5′ rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Computational methods for annotation of plant regulatory non-coding RNAs using RNA-seq

Plant transcriptome encompasses numerous endogenous, regulatory non-coding RNAs (ncRNAs) that play a major biological role in regulating key physiological mechanisms. While studies have shown that ncRNAs are extremely diverse and ubiquitous, the functions of the vast majority of ncRNAs are still unknown. With ever-increasing ncRNAs under study, it is essential to identify, categorize and annotate these ncRNAs on a genome-wide scale. The use of high-throughput RNA sequencing (RNA-seq) technologies provides a broader picture of the non-coding component of transcriptome, enabling the comprehensive identification and annotation of all major ncRNAs across samples. However, the detection of known and emerging class of ncRNAs from RNA-seq data demands complex computational methods owing to their unique as well as similar characteristics. Here, we discuss major plant endogenous, regulatory ncRNAs in an RNA sample followed by computational strategies applied to discover each class of ncRNAs using RNA-seq. We also provide a collection of relevant software packages and databases to present a comprehensive bioinformatics toolbox for plant ncRNA researchers. We assume that the discussions in this review will provide a rationale for the discovery of all major categories of plant ncRNAs.

Genome-wide analysis of epigenetic and transcriptional changes associated with heterosis in pigeonpea

Hybrids are extensively used in agriculture to deliver an increase in yield, yet the molecular basis of heterosis is not well understood. Global DNA methylation analysis, transcriptome analysis and small RNA profiling were aimed to understand the epigenetic effect of the changes in gene expression level in the two hybrids and their parental lines. Increased DNA methylation was observed in both the hybrids as compared to their parents. This increased DNA methylation in hybrids showed that majority of the 24-nt siRNA clusters had higher expression in hybrids than the parents. Transcriptome analysis revealed that various phytohormones (auxin and salicylic acid) responsive hybrid-MPV DEGs were significantly altered in both the hybrids in comparison to MPV. DEGs associated with plant immunity and growth were overexpressed whereas DEGs associated with basal defence level were repressed. This antagonistic patterns of gene expression might contribute to the greater growth of the hybrids. It was also noticed that some common as well as unique changes in the regulatory pathways were associated with heterotic growth in both the hybrids. Approximately 70% and 67% of down-regulated hybrid-MPV DEGs were found to be differentially methylated in ICPH 2671 and ICPH 2740 hybrid, respectively. This reflected the association of epigenetic regulation in altered gene expressions. Our findings also revealed that miRNAs might play important roles in hybrid vigour in both the hybrids by regulating their target genes, especially in controlling plant growth and development, defence and stress response pathways. The above finding provides an insight into the molecular mechanism of pigeonpea heterosis.

Identification and expression analysis of miRNAs and elucidation of their role in salt tolerance in rice varieties susceptible and tolerant to salinity

Soil salinization is a serious problem for cultivation of rice, as among cereals rice is the most salt sensitive crop, and more than 40% of the total agricultural land amounting to approximately 80 million ha the world over is salt affected. Salinity affects a plant in a varieties of ways, including ion toxicity, osmotic stress and oxidative damage. Since miRNAs occupy the top place in biochemical events determining a trait, understanding their role in salt tolerance is highly desirable, which may allow introduction of the trait in the rice cultivars of choice through biotechnological interventions. High throughput sequencing of sRNAs in the root and shoot tissues of the seedlings of the control and NaCl treated Pokkali, a salt-tolerant rice variety, identified 75 conserved miRNAs and mapped 200 sRNAs to the rice genome as novel miRNAs. Expression of nine novel miRNAs and two conserved miRNAs were confirmed by Northern blotting. Several of both conserved and novel miRNAs that expressed differentially in root and/or shoot tissues targeted transcription factors like AP2/EREBP domain protein, ARF, NAC, MYB, NF-YA, HD-Zip III, TCP and SBP reported to be involved in salt tolerance or in abiotic stress tolerance in general. Most of the novel miRNAs expressed in the salt tolerant wild rice Oryza coarctata, suggesting conservation of miRNAs in taxonomically related species. One of the novel miRNAs, osa-miR12477, also targeted L-ascorbate oxidase (LAO), indicating build-up of oxidative stress in the plant upon salt treatment, which was confirmed by DAB staining. Thus, salt tolerance might involve miRNA-mediated regulation of 1) cellular abundance of the hormone signaling components like EREBP and ARF, 2) synthesis of abiotic stress related transcription factors, and 3) antioxidative component like LAO for mitigation of oxidative damage. The study clearly indicated importance of osa-miR12477 regulated expression of LAO in salt tolerance in the plant.

Genome-scale, single-cell-type resolution of microRNA activities within a whole plant organ

Loaded into ARGONAUTE(AGO) proteins, eukaryotic micro(mi)RNAs regulate gene expression via cleavage, translational repression, and/or accelerated decay of sequence-complementary target transcripts. Despite their importance in development, cell identity maintenance and stress responses, how individual miRNAs contribute to spatial gene regulation within the complex cell mosaics formed in tissues/organs has remained inaccessible in any organism to date. We have developed a non-invasive methodology to examine, at single-cell-type resolution, the AGO-loading and activity patterns of entire miRNA cohorts in intact organs, applied here to the Arabidopsis root tip. A dual miRNAome-targetome analytical interface allowing intuitive data integration/visualization was developed as the basis for in-depth investigations via single-cell-type experimentation. These uncovered an array of so far speculative or hitherto unknown types of spatial miRNA-mediated gene regulation schemes, including via widespread cell-to-cell movement between contiguous layers of distinct identities. This study provides the proof of principle that minimally invasive, genome-scale analysis of miRNA activities within and between single-cell types of whole organs is achievable.

Transcriptome-wide identification and functional investigation of the RDR2- and DCL3-dependent small RNAs encoded by long non-coding RNAs in Arabidopsis thaliana

The involvement of the long non-coding RNAs (lncRNAs) in small RNA (sRNA)-related pathways remains elusive. Taking advantage of the public sRNA sequencing data, we searched for RNA-dependent RNA polymerase 2 (RDR2)- and Dicer-like 3 (DCL3)-dependent sRNAs generated from the lncRNAs of Arabidopsis thaliana. First, 55,162 sRNAs were identified to be RDR2- and DCL3-dependent. These sRNAs were then mapped onto the lncRNAs. As a result, a total of 26,643 sRNAs found their loci on 3,834 lncRNAs, and 29,388 sRNAs found their loci on 4,174 reverse complementary (RC) sequences of the lncRNAs. To support the formation of the double-stranded precursors for sRNA generation, double-stranded RNA sequencing (dsRNA-seq) reads were mapped onto the sense and antisense strands of the lncRNAs with RDR2- and DCL3-dependent sRNA loci. As a result, 1,075 regions longer than 100 nt were identified to be covered by dsRNA-seq reads on 390 sense strands of the lncRNAs, and 1,352 regions were identified on 544 RC strands. Besides, 2,238 out of 3,211 dsRNA-seq read-covered sRNA loci were supported by degradome sequencing data on the sense strands of the lncRNAs. Interestingly, dozens of dsRNA-seq read-covered regions with AGO4-associated sRNA loci showed site-specific chromatin modification patterns. Thus, some of the lncRNAs were integrated into the RDR2- and DCL3-dependent sRNA biogenesis pathway. Moreover, our results indicated that the site-specific chromatin modifications mediated by the AGO4-associated sRNAs might play a regulatory role on the transcription activity of the lncRNA genes.

Conserved MicroRNA Act Boldly During Sprout Development and Quality Formation in Pingyang Tezaocha (Camellia sinensis)

Tea tree [Camellia sinensis (L.) O. Kuntze] is an important leaf (sometimes tender stem)-using commercial plant with many medicinal uses. The development of newly sprouts would directly affect the yield and quality of tea product, especially significant for Pingyang Tezaocha (PYTZ) which takes up a large percent in the early spring tea market. MicroRNA (miRNA), particularly the conserved miRNAs, often position in the center of subtle and complex gene regulatory systems, precisely control the biological processes together with other factors in a spatio-temporal pattern. Here, quality-determined metabolites catechins, theanine and caffeine in PYTZ sprouts including buds (sBud), different development stages of leaves (sL1, sL2) and stems (sS1, sS2) were quantified. A total of 15 miRNA libraries of the same tissue with three repetitions for each were constructed to explore vital miRNAs during the biological processes of development and quality formation. We analyzed the whole miRNA profiles during the sprout development and defined conserved miRNA families in the tea plant. The differentially expressed miRNAs related to the expression profiles buds, leaves, and stems development stages were described. Twenty one miRNAs and eight miRNA-TF pairs that most likely to participate in regulating development, and at least two miRNA-TF-metabolite triplets that participate in both development and quality formation had been filtered. Our results indicated that conserved miRNA act boldly during important biological processes, they are (i) more likely to be linked with morphological function in primary metabolism during sprout development, and (ii) hold an important position in secondary metabolism during quality formation in tea plant, also (iii) coordinate with transcription factors in forming networks of complex multicellular organism regulation.

Directions for research and training in plant omics: Big Questions and Big Data

A key remit of the NSF-funded "Arabidopsis Research and Training for the 21st Century" (ART-21) Research Coordination Network has been to convene a series of workshops with community members to explore issues concerning research and training in plant biology, including the role that research using Arabidopsis thaliana can play in addressing those issues. A first workshop focused on training needs for bioinformatic and computational approaches in plant biology was held in 2016, and recommendations from that workshop have been published (Friesner et al., Plant Physiology, 175, 2017, 1499). In this white paper, we provide a summary of the discussions and insights arising from the second ART-21 workshop. The second workshop focused on experimental aspects of omics data acquisition and analysis and involved a broad spectrum of participants from academics and industry, ranging from graduate students through post-doctorates, early career and established investigators. Our hope is that this article will inspire beginning and established scientists, corporations, and funding agencies to pursue directions in research and training identified by this workshop, capitalizing on the reference species Arabidopsis thaliana and other valuable plant systems.

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.

Identification and characterization of known and novel microRNAs in strawberry fruits induced by Botrytis cinerea

MicroRNAs are endogenous small non-coding RNAs that negatively regulate mRNAs, mainly at the post-transcriptional level, and play an important role in resistance response of plants. To date, there are few reports on resistance response of strawberry miRNAs to pathogens. In this study, using high-throughput sequencing, 134 conserved and 35 novel miRNAs were identified in six libraries within the treatment of Botrytis cinerea. A total 497 potential target genes were predicted using Fragaria vesca genome. Most of the differential expressed miRNAs in strawberry fruits were up-regulated in early libraries and down-regulated in late libraries. PIRL, the target gene of miR5290a, showed the opposite expressed trend compared with miR5290 from T1 to T3 libraries, and functional analysis of the PIRL gene shows that it has obvious resistance to B. cinerea in the strawberry fruits with overexpressed PIRL gene. We speculate that miR5290a negatively regulates its target gene PIRL to increase resistance to pathogen infection, and further analysis of PIRL function is meaningful for studying the plant-pathogen relationship and improving strawberry fruit quality and yield.

Argonaute Proteins and Mechanisms of RNA Interference in Eukaryotes and Prokaryotes

Noncoding RNAs play essential roles in genetic regulation in all organisms. In eukaryotic cells, many small noncoding RNAs act in complex with Argonaute proteins and regulate gene expression by recognizing complementary RNA targets. The complexes of Argonaute proteins with small RNAs also play a key role in silencing of mobile genetic elements and, in some cases, viruses. These processes are collectively called RNA interference. RNA interference is a powerful tool for specific gene silencing in both basic research and therapeutic applications. Argonaute proteins are also found in prokaryotic organisms. Recent studies have shown that prokaryotic Argonautes can also cleave their target nucleic acids, in particular DNA. This activity of prokaryotic Argonautes might potentially be used to edit eukaryotic genomes. However, the molecular mechanisms of small nucleic acid biogenesis and the functions of Argonaute proteins, in particular in bacteria and archaea, remain largely unknown. Here we briefly review available data on the RNA interference processes and Argonaute proteins in eukaryotes and prokaryotes.

Regulation and Evolution of NLR Genes: A Close Interconnection for Plant Immunity

NLR (NOD-like receptor) genes belong to one of the largest gene families in plants. Their role in plants' resistance to pathogens has been clearly described for many members of this gene family, and dysregulation or overexpression of some of these genes has been shown to induce an autoimmunity state that strongly affects plant growth and yield. For this reason, these genes have to be tightly regulated in their expression and activity, and several regulatory mechanisms are described here that tune their gene expression and protein levels. This gene family is subjected to rapid evolution, and to maintain diversity at NLRs, a plethora of genetic mechanisms have been identified as sources of variation. Interestingly, regulation of gene expression and evolution of this gene family are two strictly interconnected aspects. Indeed, some examples have been reported in which mechanisms of gene expression regulation have roles in promotion of the evolution of this gene family. Moreover, co-evolution of the NLR gene family and other gene families devoted to their control has been recently demonstrated, as in the case of miRNAs.

Revisiting Criteria for Plant MicroRNA Annotation in the Era of Big Data

MicroRNAs (miRNAs) are ∼21-nucleotide-long regulatory RNAs that arise from endonucleolytic processing of hairpin precursors. Many function as essential posttranscriptional regulators of target mRNAs and long noncoding RNAs. Alongside miRNAs, plants also produce large numbers of short interfering RNAs (siRNAs), which are distinguished from miRNAs primarily by their biogenesis (typically processed from long double-stranded RNA instead of single-stranded hairpins) and functions (typically via roles in transcriptional regulation instead of posttranscriptional regulation). Next-generation DNA sequencing methods have yielded extensive data sets of plant small RNAs, resulting in many miRNA annotations. However, it has become clear that many miRNA annotations are questionable. The sheer number of endogenous siRNAs compared with miRNAs has been a major factor in the erroneous annotation of siRNAs as miRNAs. Here, we provide updated criteria for the confident annotation of plant miRNAs, suitable for the era of "big data" from DNA sequencing. The updated criteria emphasize replication and the minimization of false positives, and they require next-generation sequencing of small RNAs. We argue that improved annotation systems are needed for miRNAs and all other classes of plant small RNAs. Finally, to illustrate the complexities of miRNA and siRNA annotation, we review the evolution and functions of miRNAs and siRNAs in plants.

Small RNA differential expression and regulation in Tuxpeño maize embryogenic callus induction and establishment

Somatic embryogenesis represents an alternative developmental process used to achieve genetic transformation and to approach key questions in maize development. It is known that embryogenic callus induction and plant regeneration are accompanied by microRNA expression changes. However, small RNA (sRNA) populations have not been explored during the proliferative callus subculture establishment and their impact on maintaining the dedifferentiated status and embryogenic potential is far from being completely understood. Here we globally tested the sRNA populations in explants (immature embryos), induced and established maize embryogenic callus from the Mexican cultivar VS-535, Tuxpeño landrace. We detected readjustments in 24 nt and 21-22 nt sRNAs during the embryogenic callus (EC) establishment and maintenance. A follow up on specific microRNAs (miRNAs) indicated that miRNAs related to stress response substantially increase upon the callus proliferation establishment, correlating with a reduction in some of their target levels. On the other hand, while 24 nt-long heterochromatic small interfering RNAs (hc-siRNAs) derived from transposable retroelements transiently decreased in abundance during the EC establishment, a population of 22 nt-hc-siRNAs increased. This was accompanied by reduction in transposon expression in the established callus subcultures. We conclude that stress- and development-related miRNAs are highly expressed upon maize EC callus induction and during maintenance of the subcultures, while miRNAs involved in hormone response only transiently increase during induction. In addition, the establishment of a proliferative status in embryogenic callus is accompanied by important readjustments in hc-siRNAs mapping to long tandem repeat (LTR) retrotransposons, and their expression regulation.

Small RNA Functions Are Required for Growth and Development of Magnaporthe oryzae

RNA interference (RNAi) is conserved in eukaryotic organisms, and it has been well studied in many animal and plant species and some fungal species, yet it is not well studied in fungal plant pathogens. In the rice blast fungus Magnaporthe oryzae, we examined small RNA (sRNA) and their biogenesis in the context of growth and pathogenicity. Through genetic and genomic analyses, we demonstrate that loss of a single gene encoding Dicer, RNA-dependent RNA polymerase, or Argonaute reduces sRNA levels. These three proteins are required for the biogenesis of sRNA-matching genome-wide regions (coding regions, repeats, and intergenic regions). The loss of one Argonaute reduced both sRNA and fungal virulence on barley leaves. Transcriptome analysis of multiple mutants revealed that sRNA play an important role in transcriptional regulation of repeats and intergenic regions in M. oryzae. Together, these data support that M. oryzae sRNA regulate developmental processes including, fungal growth and virulence.

Biogenesis of phased siRNAs on membrane-bound polysomes in Arabidopsis

Small RNAs are central players in RNA silencing, yet their cytoplasmic compartmentalization and the effects it may have on their activities have not been studied at the genomic scale. Here we report that Arabidopsis microRNAs (miRNAs) and small interfering RNAs (siRNAs) are distinctly partitioned between the endoplasmic reticulum (ER) and cytosol. All miRNAs are associated with membrane-bound polysomes (MBPs) as opposed to polysomes in general. The MBP association is functionally linked to a deeply conserved and tightly regulated activity of miRNAs - production of phased siRNAs (phasiRNAs) from select target RNAs. The phasiRNA precursor RNAs, thought to be noncoding, are on MBPs and are occupied by ribosomes in a manner that supports miRNA-triggered phasiRNA production, suggesting that ribosomes on the rough ER impact siRNA biogenesis. This study reveals global patterns of cytoplasmic partitioning of small RNAs and expands the known functions of ribosomes and ER.

Epigenetic regulation and epigenomic landscape in rice

Epigenetic regulation has been implicated in the control of complex agronomic traits in rice (Oryza sativa), a staple food crop and model monocot plant. Recent advances in high-throughput sequencing and the moderately complex genome of rice have made it possible to study epigenetic regulation in rice on a genome-wide scale. This review discusses recent advances in our understanding of epigenetic regulation in rice, with an emphasis on the roles of key epigenetic regulators, the epigenomic landscape, epigenetic variation, transposon repression, and plant development.

The Diversification of Plant NBS-LRR Defense Genes Directs the Evolution of MicroRNAs That Target Them

High expression of plant nucleotide binding site leucine-rich repeat (NBS-LRR) defense genes is often lethal to plant cells, a phenotype perhaps associated with fitness costs. Plants implement several mechanisms to control the transcript level of NBS-LRR defense genes. As negative transcriptional regulators, diverse miRNAs target NBS-LRRs in eudicots and gymnosperms. To understand the evolutionary benefits of this miRNA-NBS-LRR regulatory system, we investigated the NBS-LRRs of 70 land plants, coupling this analysis with extensive small RNA data. A tight association between the diversity of NBS-LRRs and miRNAs was found. The miRNAs typically target highly duplicated NBS-LRRs In comparison, families of heterogeneous NBS-LRRs were rarely targeted by miRNAs in Poaceae and Brassicaceae genomes. We observed that duplicated NBS-LRRs from different gene families periodically gave birth to new miRNAs. Most of these newly emerged miRNAs target the same conserved, encoded protein motif of NBS-LRRs, consistent with a model of convergent evolution for these miRNAs. By assessing the interactions between miRNAs and NBS-LRRs, we found nucleotide diversity in the wobble position of the codons in the target site drives the diversification of miRNAs. Taken together, we propose a co-evolutionary model of plant NBS-LRRs and miRNAs hypothesizing how plants balance the benefits and costs of NBS-LRR defense genes.

Arabidopsis RNASE THREE LIKE2 Modulates the Expression of Protein-Coding Genes via 24-Nucleotide Small Interfering RNA-Directed DNA Methylation

RNaseIII enzymes catalyze the cleavage of double-stranded RNA (dsRNA) and have diverse functions in RNA maturation. Arabidopsis thaliana RNASE THREE LIKE2 (RTL2), which carries one RNaseIII and two dsRNA binding (DRB) domains, is a unique Arabidopsis RNaseIII enzyme resembling the budding yeast small interfering RNA (siRNA)-producing Dcr1 enzyme. Here, we show that RTL2 modulates the production of a subset of small RNAs and that this activity depends on both its RNaseIII and DRB domains. However, the mode of action of RTL2 differs from that of Dcr1. Whereas Dcr1 directly cleaves dsRNAs into 23-nucleotide siRNAs, RTL2 likely cleaves dsRNAs into longer molecules, which are subsequently processed into small RNAs by the DICER-LIKE enzymes. Depending on the dsRNA considered, RTL2-mediated maturation either improves (RTL2-dependent loci) or reduces (RTL2-sensitive loci) the production of small RNAs. Because the vast majority of RTL2-regulated loci correspond to transposons and intergenic regions producing 24-nucleotide siRNAs that guide DNA methylation, RTL2 depletion modifies DNA methylation in these regions. Nevertheless, 13% of RTL2-regulated loci correspond to protein-coding genes. We show that changes in 24-nucleotide siRNA levels also affect DNA methylation levels at such loci and inversely correlate with mRNA steady state levels, thus implicating RTL2 in the regulation of protein-coding gene expression.

Comparative Profiling of miRNAs and Target Gene Identification in Distant-Grafting between Tomato and Lycium (Goji Berry)

Local translocation of small RNAs between cells is proved. Long distance translocation between rootstock and scion is also well documented in the homo-grafting system, but the process in distant-grafting is widely unexplored where rootstock and scion belonging to different genera. Micro RNAs are a class of small, endogenous, noncoding, gene silencing RNAs that regulate target genes of a wide range of important biological pathways in plants. In this study, tomato was grafted onto goji (Lycium chinense Mill.) to reveal the insight of miRNAs regulation and expression patterns within a distant-grafting system. Goji is an important traditional Chinese medicinal plant with enriched phytochemicals. Illumina sequencing technology has identified 68 evolutionary known miRNAs of 37 miRNA families. Moreover, 168 putative novel miRNAs were also identified. Compared with control tomato, 43 (11 known and 32 novels) and 163 (33 known and 130 novels) miRNAs were expressed significantly different in shoot and fruit of grafted tomato, respectively. The fruiting stage was identified as the most responsive in the distant-grafting approach and 123 miRNAs were found as up-regulating in the grafted fruit which is remarkably higher compare to the grafted shoot tip (28). Potential targets of differentially expressed miRNAs were found to be involved in diverse metabolic and regulatory pathways. ADP binding activities, molybdopterin synthase complex and RNA helicase activity were found as enriched terms in GO (Gene Ontology) analysis. Additionally, "metabolic pathways" was revealed as the most significant pathway in KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis. The information of the small RNA transcriptomes that are obtained from this study might be the first miRNAs elucidation for a distant-grafting system, particularly between goji and tomato. The results from this study will provide the insights into the molecular aspects of miRNA-mediated regulation in the medicinal plant goji, and in grafted tomato. Noteworthy, it would provide a basis how miRNA signals could exchange between rootstock and scion, and the relevance to diverse biological processes.

Plants Encode a General siRNA Suppressor That Is Induced and Suppressed by Viruses

Small RNAs play essential regulatory roles in genome stability, development, and responses to biotic and abiotic stresses in most eukaryotes. In plants, the RNaseIII enzyme DICER-LIKE1 (DCL1) produces miRNAs, whereas DCL2, DCL3, and DCL4 produce various size classes of siRNAs. Plants also encode RNASE THREE-LIKE (RTL) enzymes that lack DCL-specific domains and whose function is largely unknown. We found that virus infection induces RTL1 expression, suggesting that this enzyme could play a role in plant–virus interaction. To first investigate the biochemical activity of RTL1 independent of virus infection, small RNAs were sequenced from transgenic plants constitutively expressing RTL1. These plants lacked almost all DCL2-, DCL3-, and DCL4-dependent small RNAs, indicating that RTL1 is a general suppressor of plant siRNA pathways. In vivo and in vitro assays revealed that RTL1 prevents siRNA production by cleaving dsRNA prior to DCL2-, DCL3-, and DCL4-processing. The substrate of RTL1 cleavage is likely long-perfect (or near-perfect) dsRNA, consistent with the RTL1-insensitivity of miRNAs, which derive from DCL1-processing of short-imperfect dsRNA. Virus infection induces RTL1 mRNA accumulation, but viral proteins that suppress RNA silencing inhibit RTL1 activity, suggesting that RTL1 has evolved as an inducible antiviral defense that could target dsRNA intermediates of viral replication, but that a broad range of viruses counteract RTL1 using the same protein toolbox used to inhibit antiviral RNA silencing. Together, these results reveal yet another level of complexity in the evolutionary battle between viruses and plant defenses.

Viral infection of plants triggers the expression of an RNaseIII enzyme that represses the production of siRNAs by cleaving their long dsRNA precursors. Read the accompanying Synopsis.

Most eukaryotes produce essential regulatory molecules called small RNAs. These molecules are produced primarily by a class of RNaseIII enzymes called DICER, which excises small RNA duplexes from long double-stranded (ds)RNA precursor molecules. Plants also encode several RNaseIII enzymes called RNASE THREE-LIKE (RTL), but the function of these proteins is largely unknown. Here, we show that RTL1 represses small RNA production by cleaving dsRNA before DICER can process them. RTL1 appears to specifically act on the templates of a class of small RNAs called siRNAs, but not on miRNA precursors, suggesting that it cleaves long-perfect (or near-perfect) dsRNA, but not short-imperfect dsRNA. We also found that RTL1 expression is induced after virus infection, suggesting that RTL1 could act as an inducible antiviral defense by destroying dsRNA intermediates of viral replication. Our findings suggest that viruses have evolved to inhibit RTL1 activity, ultimately resulting in successful viral infection.

Promoter-based identification of novel non-coding RNAs reveals the presence of dicistronic snoRNA-miRNA genes in Arabidopsis thaliana

Background

In the past few decades, non-coding RNAs (ncRNAs) have emerged as important regulators of gene expression in eukaryotes. Most studies of ncRNAs in plants have focused on the identification of silencing microRNAs (miRNAs) and small interfering RNAs (siRNAs). Another important family of ncRNAs that has been well characterized in plants is the small nucleolar RNAs (snoRNAs) and the related small Cajal body-specific RNAs (scaRNAs). Both target chemical modifications of ribosomal RNAs (rRNAs) and small nuclear RNAs (snRNAs). In plants, the snoRNA genes are organized in clusters, transcribed by RNA Pol II from a common promoter and subsequently processed into mature molecules. The promoter regions of snoRNA polycistronic genes in plants are highly enriched in two conserved cis-regulatory elements (CREs), Telo-box and Site II, which coordinate the expression of snoRNAs and ribosomal protein coding genes throughout the cell cycle.

Results

In order to identify novel ncRNA genes, we have used the snoRNA Telo-box/Site II motifs combination as a functional promoter indicator to screen the Arabidopsis genome. The predictions generated by this process were tested by detailed exploration of available RNA-Seq and expression data sets and experimental validation. As a result, we have identified several snoRNAs, scaRNAs and 'orphan' snoRNAs. We also show evidence for 16 novel ncRNAs that lack similarity to any reported RNA family. Finally, we have identified two dicistronic genes encoding precursors that are processed to mature snoRNA and miRNA molecules. We discuss the evolutionary consequences of this result in the context of a tight link between snoRNAs and miRNAs in eukaryotes.

Conclusions

We present an alternative computational approach for non-coding RNA detection. Instead of depending on sequence or structure similarity in the whole genome screenings, we have explored the properties of promoter regions of well-characterized ncRNAs. Interestingly, besides expected ncRNAs predictions we were also able to recover single precursor arrangement for snoRNA-miRNA. Accompanied by analyses performed on rice sequences, we conclude that such arrangement might have interesting functional and evolutionary consequences and discuss this result in the context of a tight link between snoRNAs and miRNAs in eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2221-x) contains supplementary material, which is available to authorized users.

Identification of miRNAs with potential roles in regulation of anther development and male-sterility in 7B-1 male-sterile tomato mutant

BACKGROUND

The 7B-1 tomato line (Solanum lycopersicum cv. Rutgers) is a photoperiod-sensitive male-sterile mutant, with potential application in hybrid seed production. Small RNAs (sRNAs) in tomato have been mainly characterized in fruit development and ripening, but none have been studied with respect to flower development and regulation of male-sterility. Using sRNA sequencing, we identified miRNAs that are potentially involved in anther development and regulation of male-sterility in 7B-1 mutant.

RESULTS

Two sRNA libraries from 7B-1 and wild type (WT) anthers were sequenced and thirty two families of known miRNAs and 23 new miRNAs were identified in both libraries. MiR390, miR166, miR159 were up-regulated and miR530, miR167, miR164, miR396, miR168, miR393, miR8006 and two new miRNAs, miR#W and miR#M were down-regulated in 7B-1 anthers. Ta-siRNAs were not differentially expressed and likely not associated with 7B-1 male-sterility. miRNA targets with potential roles in anther development were validated using 5'-RACE. QPCR analysis showed differential expression of miRNA/target pairs of interest in anthers and stem of 7B-1, suggesting that they may regulate different biological processes in these tissues. Expression level of most miRNA/target pairs showed negative correlation, except for few. In situ hybridization showed predominant expression of miR159, GAMYBL1, PMEI and cystatin in tapetum, tetrads and microspores.

CONCLUSION

Overall, we identified miRNAs with potential roles in anther development and regulation of male-sterility in 7B-1. A number of new miRNAs were also identified from tomato for the first time. Our data could be used as a benchmark for future studies of the molecular mechanisms of male-sterility in other crops.

miRNA778 and SUVH6 are involved in phosphate homeostasis in Arabidopsis

microRNAs (miRNAs) play an important role in plant adaptation to phosphate (Pi) starvation. Histone methylation can remodel chromatin structure and mediate gene expression. This study identified Arabidopsis miR778, a Pi-responsive miRNA, and its target gene Su(var) 3-9 homologs 6 (SUVH6) encoding a histone H3 lysine 9 (H3K9) methyltransferase. Overexpression of miR778 moderately enhanced primary and lateral root growth, free phosphate accumulation in shoots, and accumulation of anthocyanin under Pi deficient conditions. miR778 overexpression relieved the arrest of columella cell development under Pi starvation. Conversely, transgenic plants overexpressing a miR778-target mimic (35S::MIM778), that act as a sponge and sequesters miR778, showed opposite phenotypes of 35S::miR778 plants under Pi deficiency. Expression of several Pi deficiency-responsive genes such as miR399, Phosphate Transporter (PHT1;4), Low Phosphate-Resistant1 (LPR1) and Production of Anthocyanin Pigment 1 (PAP1) were elevated in the miR778 overexpressing plants, suggesting that both miR778 and SUVH6 are involved in phosphate homeostasis in plants. This study has provided a basis for further investigation on how SUVH6 regulates its downstream genes through chromatin remodeling and DNA methylation in plants stressed by Pi deficiency.

JMJ24 binds to RDR2 and is required for the basal level transcription of silenced loci in Arabidopsis

Transposable elements (TEs) and repetitive sequences are ubiquitously present in eukaryotic genomes which are in general epigenetically silenced by DNA methylation and/or histone 3 lysine 9 methylation (H3K9me). RNA-directed DNA methylation (RdDM) is the major pathway that initiates de novo DNA methylation in Arabidopsis and sets up a self-reinforcing silencing loop between DNA methylation and H3K9me. However, a key issue is the requirement of a basal level transcript from the target loci to initiate the RNA-based silencing. How the heterochromatic silenced loci are transcribed remains largely unknown. Here, we show that JMJ24, a JmjC domain-containing protein counteracts H3K9me to promote basal level transcription of endogenous silenced loci in Arabidopsis. JMJ24 functionally resembles the fission yeast JmjC protein Epe1. The transcript promoted by JMJ24 is, at least in part, processed to small RNA to initiate the RdDM. Genome-wide transcriptome profiling indicates that transcript levels of TEs are more likely regulated by JMJ24, compared with protein-coding genes. Our data suggest that JMJ24 plays a conserved role in promoting basal level transcription of endogenous silenced loci to reinforce the silencing. We also provide evidence of a physical association between JMJ24 and RNA-dependent RNA polymerase 2 (RDR2), which represents an evolved property of the RNA silencing pathway.

Identification of MicroRNAs in Response to Different Day Lengths in Soybean Using High-Throughput Sequencing and qRT-PCR

MicroRNAs (miRNAs) are short, non-coding single-strand RNA molecules that play important roles in plant growth, development and stress responses. Flowering time affects the seed yield and quality of soybean. However, the miRNAs involved in the regulation of flowering time in soybean have not been reported until recently. Here, high-throughput sequencing and qRT-PCR were used to identify miRNAs involved in soybean photoperiodic pathways. The first trifoliate leaves of soybean that receive the signal of light treatment were used to construct six libraries (0, 8, and 16 h under short-day (SD) treatment and 0, 8, and 16 h under long-day (LD) treatment). The libraries were sequenced using Illumina Solexa. A total of 318 known plant miRNAs belonging to 163 miRNA families and 81 novel predicted miRNAs were identified. Among these, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h, showed differences in abundance between LD and SD treatments. Furthermore, the results of GO and KEGG analyses indicated that most of the miRNA targets were transcription factors. Seven miRNAs at 0 h, 23 miRNAs (including four novel predicted miRNAs) at 8 h, 16 miRNAs (including one novel predicted miRNA) at 16 h and miRNA targets were selected for qRT-PCR analysis to assess the accuracy of the sequencing and target prediction. The results indicated that the expression patterns of the selected miRNAs and miRNA targets showed no differences between the qRT-PCR and sequencing results. In addition, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h responded to day length changes in soybean, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h. These results provided an important molecular basis to understand the regulation of flowering time through photoperiodic pathways in soybean.

Are ta-siRNAs only originated from the cleavage site of miRNA on its target RNAs and phased in 21-nt increments?

Trans-acting siRNAs (ta-siRNAs) are a class of small RNAs playing crucial roles in the regulation of plant gene expression. According to the canonical model, specific miRNA-guided cleavage of a TAS transcript triggers and sets the registry for the subsequent production of ta-siRNAs at 21-nt increments from the cleavage site. However, a previously validated 22-nt ta-siR2140 indicated that ta-siRNAs might be initiated from other phase increments and registers, which resulted in massive ta-siRNAs missing in the canonical model. To test this hypothesis, we employed high-throughput sequencing data to thoroughly identify the miR173-triggered ta-siRNAs from TAS1/TAS2 transcripts. As a result, thousands of phased siRNAs not generated through the canonical pathway were identified and 110 novel siRNA-target interactions were further validated based on degradome sequencing data. Based on these results, we propose that the canonical biogenesis model of ta-siRNAs should be modified in order to recruit the previously unidentified ta-siRNA candidates.

siRNA-directed DNA Methylation in Plants

DNA cytosine methylationis an important epigenetic process that is correlated with transgene silencing, transposon suppression, and gene imprinting. In plants, small interfering RNAs (siRNAs) can trigger DNA methylation at loci containing their homolog sequences through a process called RNA-directed DNA methylation (RdDM). In canonical RdDM, 24 nucleotide (nt) siRNAs (ra-siRNAs) will be loaded into their effector protein called ARGONAUTE 4 (AGO4) and subsequently targeted to RdDM loci through base-pairing with the non-coding transcripts produced by DNA-directed RNA Polymerase V. Then, the AGO4-ra-siRNA will recruit the DNA methyltransferase to catalyze de novo DNA methylation. Recent studies also identified non-canonical RdDM pathways that involve microRNAs or 21 nt siRNAs. These RdDM pathways are biologically important since they control responses biotic and abiotic stresses, maintain genome stability and regulate development. Here, we summarize recent pro-gresses of mechanisms governing canonical and non-canonical RdDM pathways.

Intronic regions of plant genes potentially encode RDR (RNA-dependent RNA polymerase)-dependent small RNAs

Recent research has linked the non-coding intronic regions of plant genes to the production of small RNAs (sRNAs). Certain introns, called 'mirtrons' and 'sirtrons', could serve as the single-stranded RNA precursors for the generation of microRNA and small interfering RNA, respectively. However, whether the intronic regions could serve as the template for double-stranded RNA synthesis and then for sRNA biogenesis through an RDR (RNA-dependent RNA polymerase)-dependent pathway remains unclear. In this study, a genome-wide search was made for the RDR-dependent sRNA loci within the intronic regions of the Arabidopsis genes. Hundreds of intronic regions encoding three or more RDR-dependent sRNAs were found to be covered by dsRNA-seq (double-stranded RNA sequencing) reads, indicating that the intron-derived sRNAs were indeed generated from long double-stranded RNA precursors. More interestingly, phase-distributed sRNAs were discovered on some of the dsRNA-seq read-covered intronic regions, and those sRNAs were largely 24 nt in length. Based on these results, the opinion is put forward that the intronic regions might serve as the genomic origins for the RDR-dependent sRNAs. This opinion might add a novel layer to the current biogenesis model of the intron-derived sRNAs.

Primer-dependent and primer-independent initiation of double stranded RNA synthesis by purified Arabidopsis RNA-dependent RNA polymerases RDR2 and RDR6

Cellular RNA-dependent RNA polymerases (RDRs) are fundamental components of RNA silencing in plants and many other eukaryotes. In Arabidopsis thaliana genetic studies have demonstrated that RDR2 and RDR6 are involved in the synthesis of double stranded RNA (dsRNA) from single stranded RNA (ssRNA) targeted by RNA silencing. The dsRNA is subsequently cleaved by the ribonuclease DICER-like into secondary small interfering RNAs (siRNAs) that reinforce and/or maintain the silenced state of the target RNA. Models of RNA silencing propose that RDRs could use primer-independent and primer-dependent initiation to generate dsRNA from a transcript targeted by primary siRNA or microRNA (miRNA). However, the biochemical activities of RDR proteins are still partly understood. Here, we obtained active recombinant RDR2 and RDR6 in a purified form. We demonstrate that RDR2 and RDR6 have primer-independent and primer-dependent RNA polymerase activities with different efficiencies. We further show that RDR2 and RDR6 can initiate dsRNA synthesis either by elongation of 21- to 24- nucleotides RNAs hybridized to complementary RNA template or by elongation of self-primed RNA template. These findings provide new insights into our understanding of the molecular mechanisms of RNA silencing in plants.

Evolutionary patterns and coevolutionary consequences of MIRNA genes and microRNA targets triggered by multiple mechanisms of genomic duplications in soybean

The evolutionary dynamics of duplicated protein-encoding genes (PEGs) is well documented. However, the evolutionary patterns and consequences of duplicated MIRNAs and the potential influence on the evolution of their PEG targets are poorly understood. Here, we demonstrate the evolution of plant MIRNAs subsequent to a recent whole-genome duplication. Overall, the retention of MIRNA duplicates was correlated to the retention of adjacent PEG duplicates, and the retained MIRNA duplicates exhibited a higher level of interspecific preservation of orthologs than singletons, suggesting that the retention of MIRNA duplicates is related to their functional constraints and local genomic stability. Nevertheless, duplication status, rather than local genic collinearity, was the primary determinant of levels of nucleotide divergence of MIRNAs. In addition, the retention of duplicated MIRNAs appears to be associated with the retention of their corresponding duplicated PEG targets. Furthermore, we characterized the evolutionary novelty of a legume-specific microRNA (miRNA) family, which resulted from rounds of genomic duplication, and consequent dynamic evolution of its NB-LRR targets, an important gene family with primary roles in plant-pathogen interactions. Together, these observations depict evolutionary patterns and novelty of MIRNAs in the context of genomic duplication and evolutionary interplay between MIRNAs and their PEG targets mediated by miRNAs.

Identification of microRNAs and target genes in the fruit and shoot tip of Lycium chinense: a traditional Chinese medicinal plant

Although Lycium chinense (goji berry) is an important traditional Chinese medicinal plant, little genome information is available for this plant, particularly at the small-RNA level. Recent findings indicate that the evolutionary role of miRNAs is very important for a better understanding of gene regulation in different plant species. To elucidate small RNAs and their potential target genes in fruit and shoot tissues, high-throughput RNA sequencing technology was used followed by qRT-PCR and RLM 5’-RACE experiments. A total of 60 conserved miRNAs belonging to 31 families and 30 putative novel miRNAs were identified. A total of 62 significantly differentially expressed miRNAs were identified, of which 15 (14 known and 1 novel) were shoot-specific, and 12 (7 known and 5 novel) were fruit-specific. Additionally, 28 differentially expressed miRNAs were recorded as up-regulated in fruit tissues. The predicted potential targets were involved in a wide range of metabolic and regulatory pathways. GO (Gene Ontology) enrichment analysis and the KEGG (Kyoto Encyclopedia of Genes and Genomes) database revealed that “metabolic pathways” is the most significant pathway with respect to the rich factor and gene numbers. Moreover, five miRNAs were related to fruit maturation, lycopene biosynthesis and signaling pathways, which might be important for the further study of fruit molecular biology. This study is the first, to detect known and novel miRNAs, and their potential targets, of L. chinense. The data and findings that are presented here might be a good source for the functional genomic study of medicinal plants and for understanding the links among diversified biological pathways.

Distinctive profiles of small RNA couple inverted repeat-induced post-transcriptional gene silencing with endogenous RNA silencing pathways in Arabidopsis

The experimental induction of RNA silencing in plants often involves expression of transgenes encoding inverted repeat (IR) sequences to produce abundant dsRNAs that are processed into small RNAs (sRNAs). These sRNAs are key mediators of post-transcriptional gene silencing (PTGS) and determine its specificity. Despite its application in agriculture and broad utility in plant research, the mechanism of IR-PTGS is incompletely understood. We generated four sets of 60 Arabidopsis plants, each containing IR transgenes expressing different configurations of uidA and CHALCONE Synthase (At-CHS) gene fragments. Levels of PTGS were found to depend on the orientation and position of the fragment in the IR construct. Deep sequencing and mapping of sRNAs to corresponding transgene-derived and endogenous transcripts identified distinctive patterns of differential sRNA accumulation that revealed similarities among sRNAs associated with IR-PTGS and endogenous sRNAs linked to uncapped mRNA decay. Detailed analyses of poly-A cleavage products from At-CHS mRNA confirmed this hypothesis. We also found unexpected associations between sRNA accumulation and the presence of predicted open reading frames in the trigger sequence. In addition, strong IR-PTGS affected the prevalence of endogenous sRNAs, which has implications for the use of PTGS for experimental or applied purposes.