A non-canonical plant microRNA target site

MicroRNA164e suppresses NAC100 transcription factor-mediated synthesis of seed storage proteins in chickpea

Development of protein-enriched chickpea varieties necessitates an understanding of specific genes and key regulatory circuits that govern the synthesis of seed storage proteins (SSPs). Here, we demonstrated the novel involvement of Ca-miR164e-CaNAC100 in regulating SSP synthesis in chickpea. Ca-miRNA164e was significantly decreased during seed maturation, especially in high-protein accessions. The miRNA was found to directly target the transactivation conferring C-terminal region of a nuclear-localized transcription factor, CaNAC100 as revealed using RNA ligase-mediated-rapid amplification of cDNA ends and target mimic assays. The functional role of CaNAC100 was demonstrated through seed-specific overexpression (NACOE) resulting in significantly augmented seed protein content (SPC) consequential to increased SSP transcription. Further, NACOE lines displayed conspicuously enhanced seed weight but reduced numbers and yield. Conversely, a downregulation of CaNAC100 and SSP transcripts was evident in seed-specific overexpression lines of Ca-miR164e that culminated in significantly lowered SPC. CaNAC100 was additionally demonstrated to transactivate the SSP-encoding genes by directly binding to their promoters as demonstrated using electrophoretic mobility shift and dual-luciferase reporter assays. Taken together, our study for the first time established a distinct role of CaNAC100 in positively influencing SSP synthesis and its critical regulation by CamiR164e, thereby serving as an understanding that can be utilized for developing SPC-rich chickpea varieties.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Target Landscape of Conserved Plant MicroRNAs and the Complexities of Their Ancient MicroRNA-Binding Sites

In plants, microRNA (miRNA)-target interactions (MTIs) require high complementarity, a feature from which bioinformatic programs have predicted numerous and diverse targets for any given miRNA, promoting the idea of complex miRNA networks. Opposing this is a hypothesis of constrained miRNA specificity, in which functional MTIs are restricted to the few targets whose required expression output is compatible with the expression of the miRNA. To explore these opposing views, the bioinformatic pipeline Targets Ranked Using Experimental Evidence was applied to strongly conserved miRNAs to identity their high-evidence (HE) targets across species. For each miRNA family, HE targets predominantly consisted of homologs from one conserved target gene family (primary family). These primary families corresponded to the known canonical miRNA-target families, validating the approach. Very few additional HE target families were identified (secondary family), and if so, they were likely functionally related to the primary family. Many primary target families contained highly conserved nucleotide sequences flanking their miRNA-binding sites that were enriched in HE homologs across species. A number of these flanking sequences are predicted to form conserved RNA secondary structures that preferentially base pair with the miRNA-binding site, implying that these sites are highly structured. Our findings support a target landscape view that is dominated by the conserved primary target families, with a minority of either secondary target families or non-conserved targets. This is consistent with the constrained hypothesis of functional miRNA specificity, which potentially in part is being facilitated by features beyond complementarity.

Small RNAs target native and cross-kingdom transcripts on both sides of the wheat stripe rust interaction

The wheat stripe rust fungus (Puccinia striiformis f.sp. tritici) threatens global wheat production. Small RNAs (sRNAs) modulate plant defense induction, and RNA exchange between host and microbe causes cross-kingdom gene silencing, but few examples are known in rust fungi. This study combined sRNA, parallel analysis of RNA ends, and gene expression data to discover sRNA-target pairs on each side of the interaction. Specific wheat 24 nt sRNAs were suppressed, while particular 35 nt fragments were strongly induced upon infection. Wheat sRNAs cleaved fungal transcripts coding for a ribosomal protein and a glycosyl hydrolase effector. Fungal microRNA-like and phased 21 nt sRNAs originated from long inverted repeats near protein coding genes. Fungal sRNAs targeted native transcripts: transposons and kinases; and cross-kingdom transcripts: a wheat nucleotide-binding domain leucine-rich repeat receptor (NLR) and multiple defense-related transcription factor families. This work sheds light on host-microbe coevolution and delivers prospects for developing pathogen control biotechnology.

Degradome Assisted Plant MicroRNA Prediction Under Alternative Annotation Criteria

Current microRNA (miRNA) prediction methods are generally based on annotation criteria that tend to miss potential functional miRNAs. Recently, new miRNA annotation criteria have been proposed that could lead to improvements in miRNA prediction methods in plants. Here, we investigate the effect of the new criteria on miRNA prediction in Arabidopsis thaliana and present a new degradome assisted functional miRNA prediction approach. We investigated the effect by applying the new criteria, and a more permissive criteria on miRNA prediction using existing miRNA prediction tools. We also developed an approach to miRNA prediction that is assisted by the functional information extracted from the analysis of degradome sequencing. We demonstrate the improved performance of degradome assisted miRNA prediction compared to unassisted prediction and evaluate the approach using miRNA differential expression analysis. We observe how the miRNA predictions fit under the different criteria and show a potential novel miRNA that has been missed within Arabidopsis thaliana. Additionally, we introduce a freely available software 'PAREfirst' that employs the degradome assisted approach. The study shows that some miRNAs could be missed due to the stringency of the former annotation criteria, and combining a degradome assisted approach with more permissive miRNA criteria can expand confident miRNA predictions.

Parallel analysis of RNA ends reveals global microRNA-mediated target RNA cleavage in maize

MicroRNAs (miRNAs) are endogenous 20-24-nucleotide non-coding RNAs that play important regulatory roles in many biological processes in eukaryotes. miRNAs modulate the expression of target genes at the post-transcriptional level by transcript cleavage or translational inhibition. The identification of miRNA target genes has been extensively investigated in Arabidopsis and rice, but an in-depth global analysis of miRNA-mediated target regulation is still lacking in maize. Here, we report a transcriptome-wide identification of miRNA targets by analyzing parallel analysis of RNA ends (PARE) datasets derived from nine different tissues at five developmental stages of the maize (Zea mays L.) B73 cultivar. In total, 246 targets corresponding to 60 miRNAs from 25 families were identified, including transcription factors and other genes. In addition, PARE analysis revealed that miRNAs guide specific target transcript cleavage in a tissue-preferential manner. Primary transcripts of MIR159c and MIR169e were found to be cleaved by mature miR159 and miR169, respectively, indicating a negative-feedback regulatory mechanism in miRNA biogenesis. Moreover, several miRNA-target gene pairs involved in seed germination were identified and experimentally validated. Our PARE analyses generated a wide and detailed miRNA-target interaction atlas, which provides a valuable resource for investigating the roles of miRNAs and their targets in maize.

Copper microRNAs modulate the formation of giant feeding cells induced by the root knot nematode Meloidogyne incognita in Arabidopsis thaliana

Root-knot nematodes (RKNs) are root endoparasites that induce the dedifferentiation of a few root cells and the reprogramming of their gene expression to generate giant hypermetabolic feeding cells. We identified two microRNA families, miR408 and miR398, as upregulated in Arabidopsis thaliana and Solanum lycopersicum roots infected by RKNs. In plants, the expression of these two conserved microRNA families is known to be activated by the SPL7 transcription factor in response to copper starvation. By combining functional approaches, we deciphered the network involving these microRNAs, their regulator and their targets. MIR408 expression was located within nematode-induced feeding cells like its regulator SPL7 and was regulated by copper. Moreover, infection assays with mir408 and spl7 knockout mutants or lines expressing targets rendered resistant to cleavage by miR398 demonstrated the essential role of the SPL7/MIR408/MIR398 module in the formation of giant feeding cells. Our findings reveal how perturbation of plant copper homeostasis, via the SPL7/MIR408/MIR398 module, modulates the development of nematode-induced feeding cells.

MicroRNA398: A Master Regulator of Plant Development and Stress Responses

MicroRNAs (miRNAs) play crucial roles in plant development and stress responses, and a growing number of studies suggest that miRNAs are promising targets for crop improvement because they participate in the regulation of diverse, important agronomic traits. MicroRNA398 (miR398) is a conserved miRNA in plants and has been shown to control multiple stress responses and plant growth in a variety of species. There are many studies on the stress response and developmental regulation of miR398. To systematically understand its function, it is necessary to summarize the evolution and functional roles of miR398 and its target genes. In this review, we analyze the evolution of miR398 in plants and outline its involvement in abiotic and biotic stress responses, in growth and development and in model and non-model plants. We summarize recent functional analyses, highlighting the role of miR398 as a master regulator that coordinates growth and diverse responses to environmental factors. We also discuss the potential for fine-tuning miR398 to achieve the goal of simultaneously improving plant growth and stress tolerance.

TRUEE; a bioinformatic pipeline to define the functional microRNA targetome of Arabidopsis

Central to plant microRNA (miRNA) biology is the identification of functional miRNA-target interactions (MTIs). However, the complementarity basis of bioinformatic target prediction results in mostly false positives, and the degree of complementarity does not equate with regulation. Here, we develop a bioinformatic workflow named TRUEE (Targets Ranked Using Experimental Evidence) that ranks MTIs on the extent to which they are subjected to miRNA-mediated cleavage. It sorts predicted targets into high (HE) and low evidence (LE) groupings based on the frequency and strength of miRNA-guided cleavage degradome signals across multiple degradome experiments. From this, each target is assigned a numerical value, termed a Category Score, ranking the extent to which it is subjected to miRNA-mediated cleavage. As a proof-of-concept, the 428 Arabidopsis miRNAs annotated in miRBase were processed through the TRUEE pipeline to determine the miRNA 'targetome'. The majority of high-ranking Category Score targets corresponded to highly conserved MTIs, validating the workflow. Very few Arabidopsis-specific, Brassicaceae-specific, or Conserved-passenger miRNAs had HE targets with high Category Scores. In total, only several hundred MTIs were found to have Category Scores characteristic of currently known physiologically significance MTIs. Although non-exhaustive, clearly the number of functional MTIs is much narrower than many studies claim. Therefore, using TRUEE to numerically rank targets  directly on experimental evidence has given insights into the scope of the functional miRNA targetome of Arabidopsis.

The Multiverse of Plant Small RNAs: How Can We Explore It?

Plant small RNAs (sRNAs) are a heterogeneous group of noncoding RNAs with a length of 20–24 nucleotides that are widely studied due to their importance as major regulators in various biological processes. sRNAs are divided into two main classes—microRNAs (miRNAs) and small interfering RNAs (siRNAs)—which differ in their biogenesis and functional pathways. Their identification and enrichment with new structural variants would not be possible without the use of various high-throughput sequencing (NGS) techniques, allowing for the detection of the total population of sRNAs in plants. Classifying sRNAs and predicting their functional role based on such high-performance datasets is a nontrivial bioinformatics task, as plants can generate millions of sRNAs from a variety of biosynthetic pathways. Over the years, many computing tools have been developed to meet this challenge. Here, we review more than 35 tools developed specifically for plant sRNAs over the past few years and explore some of their basic algorithms for performing tasks related to predicting, identifying, categorizing, and quantifying individual sRNAs in plant samples, as well as visualizing the results of these analyzes. We believe that this review will be practical for biologists who want to analyze their plant sRNA datasets but are overwhelmed by the number of tools available, thus answering the basic question of how to choose the right one for a particular study.

MiR398-regulated antioxidants contribute to Bamboo mosaic virus accumulation and symptom manifestation

Virus infections that cause mosaic or mottling in leaves commonly also induce increased levels of reactive oxygen species (ROS). However, how ROS contributes to symptoms is less well documented. Bamboo mosaic virus (BaMV) causes chlorotic mosaic symptoms in both Brachypodium distachyon and Nicotiana benthamiana. The BaMV △CPN35 mutant with an N-terminal deletion of its coat protein gene exhibits asymptomatic infection independently of virus titer. Histochemical staining of ROS in mock-, BaMV-, and BaMV△CPN35-infected leaves revealed that hydrogen peroxide (H2O2) accumulated solely in BaMV-induced chlorotic spots. Moreover, exogenous H2O2 treatment enhanced yellowish chlorosis in BaMV-infected leaves. Both BaMV and BaMV△CPN35 infection could induce the expression of Cu/Zu superoxide dismutase (CSD) antioxidants at messenger RNA and protein level. However, BaMV triggered the abundant accumulation of full-length NbCSD2 preprotein (prNbCSD2, without transit peptide cleavage), whereas BaMV△CPN35 induced a truncated prNbCSD2. Confocal microscopy showed that majority of NbCSD2-green fluorescent protein (GFP) predominantly localized in the cytosol upon BaMV infection, but BaMV△CPN35 infection tended to cause NbCSD2-GFP to remain in chloroplasts. By 5'-RNA ligase-mediated rapid amplification of cDNA ends, we validated CSDs are the targets of miR398 in vivo. Furthermore, BaMV infection increased the level of miR398, while the level of BaMV titer was regulated positively by miR398 but negatively by CSD2. In contrast, overexpression of cytosolic form NbCSD2, impairing the transport into chloroplasts, greatly enhanced BaMV accumulation. Taken together, our results indicate that induction of miR398 by BaMV infection may facilitate viral titer accumulation, and cytosolic prNbCSD2 induction may contribute to H2O2 accumulation, resulting in the development of BaMV chlorotic symptoms in plants.

An Integrated Bioinformatics and Functional Approach for miRNA Validation

MicroRNAs (miRNAs) are small (20-24 nucleotides) non-coding ribo-regulatory molecules with significant roles in regulating target mRNA and long non-coding RNAs at transcriptional and post-transcriptional levels. Rapid advancement in the small RNA sequencing methods with integration of degradome sequencing has accelerated the understanding of miRNA-mediated regulatory hubs in plants and yielded extensive annotation of miRNAs and corresponding targets. However, it is becoming clear that large numbers of such annotations are questionable. Therefore, it is imperative to adopt reliable and strict bioinformatics pipelines for miRNA identification. Furthermore, sensitive methods are needed for validation and functional characterization of miRNA and its target(s). In this chapter, we have provided a comprehensive and streamlined methodology for miRNA identification and its functional validation in plants. This includes a combination of various in silico and experimental methodologies. To identify miRNA compendium from large-scale Next-Generation Sequencing (NGS) small RNA datasets, the miR-PREFeR (miRNA PREdiction From small RNA-Seq data) bioinformatics tool has been described. Also, a homology-based search protocol for finding members of a specific miRNA family has been discussed. The chapter also includes techniques to ascertain miRNA:target pair specificity using in silico target prediction from degradome NGS libraries using CleaveLand pipeline, miRNA:target validation by in planta transient assays, 5' RLM-RACE and expression analysis as well as functional techniques like miRNA overexpression, short tandem target mimic and resistant target approaches. The proposed strategy offers a reliable and sensitive way for miRNA:target identification and validation. Additionally, we strongly promulgate the use of multiple methodologies to validate a miRNA as well as its target.

Computational Detection of MicroRNA Targets

MicroRNAs (miRNAs) are small noncoding RNAs that are recognized as posttranscriptional regulators of gene expression. These molecules have been shown to play important roles in several cellular processes. MiRNAs act on their target by guiding the RISC complex and binding to the mRNA molecule. Thus, it is recognized that the function of a miRNA is determined by the function of its target (s). By using high-throughput methodologies, novel miRNAs are being identified, but their functions remain uncharted. Target validation is crucial to properly understand the specific role of a miRNA in a cellular pathway. However, molecular techniques for experimental validation of miRNA-target interaction are expensive, time-consuming, laborious, and can be not accurate in inferring true interactions. Thus, accurate miRNA target predictions are helpful to understand the functions of miRNAs. There are several algorithms proposed for target prediction and databases containing miRNA-target information. However, these available computational tools for prediction still generate a large number of false positives and fail to detect a considerable number of true targets, which indicates the necessity of highly confident approaches to identify bona fide miRNA-target interactions. This chapter focuses on tools and strategies used for miRNA target prediction, by providing practical insights and outlooks.

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.

Overexpression of an Osa-miR162a Derivative in Rice Confers Cross-Kingdom RNA Interference-Mediated Brown Planthopper Resistance without Perturbing Host Development

Rice is a main food crop for more than half of the global population. The brown planthopper (BPH, Nilaparvata lugens) is one of the most destructive insect pests of rice. Currently, repeated overuse of chemical insecticides represents a common practice in agriculture for BPH control, which can induce insect tolerance and provoke environmental concerns. This situation calls for innovative and widely applicable strategies for rice protection against BPH. Here we report that the rice osa-miR162a can mediate cross-kingdom RNA interference (RNAi) by targeting the NlTOR (Target of rapamycin) gene of BPH that regulates the reproduction process. Through artificial diet or injection, osa-miR162a mimics repressed the NlTOR expression and impaired the oviposition of BPH adults. Consistently, overproduced osa-miR162a in transgenic rice plants compromised the fecundity of BPH adults fed with these plants, but meanwhile perturbed root and grain development. To circumvent this issue, we generated osa-miR162a-m1, a sequence-optimized osa-miR162a, by decreasing base complementarity to rice endogenous target genes while increasing base complementarity to NlTOR. Transgenic overexpression of osa-miR162a-m1 conferred rice resistance to BPH without detectable developmental penalty. This work reveals the first cross-kingdom RNAi mechanism in rice-BPH interactions and inspires a potentially useful approach for improving rice resistance to BPH. We also introduce an effective strategy to uncouple unwanted host developmental perturbation from desirable cross-kingdom RNAi benefits for overexpressed plant miRNAs.

The GRAS gene family and its roles in seed development in litchi (Litchi chinensis Sonn)

BACKGROUND

The GRAS gene family plays crucial roles in multiple biological processes of plant growth, including seed development, which is related to seedless traits of litchi (Litchi chinensis Sonn.). However, it hasn't been fully identified and analyzed in litchi, an economic fruit tree cultivated in subtropical regions.

RESULTS

In this study, 48 LcGRAS proteins were identified and termed according to their chromosomal location. LcGRAS proteins can be categorized into 14 subfamilies through phylogenetic analysis. Gene structure and conserved domain analysis revealed that different subfamilies harbored various motif patterns, suggesting their functional diversity. Synteny analysis revealed that the expansion of the GRAS family in litchi may be driven by their tandem and segmental duplication. After comprehensively analysing degradome data, we found that four LcGRAS genes belong to HAM subfamily were regulated via miR171-mediated degradation. The various expression patterns of LcGRAS genes in different tissues uncovered they were involved in different biological processes. Moreover, the different temporal expression profiles of LcGRAS genes between abortive and bold seed indicated some of them were involved in maintaining the normal development of the seed.

CONCLUSION

Our study provides comprehensive analyses on GRAS family members in litchi, insight into a better understanding of the roles of GRAS in litchi development, and lays the foundation for further investigations on litchi seed development.

Eating microRNAs: pharmacological opportunities for cross-kingdom regulation and implications in host gene and gut microbiota modulation

Cross-kingdom communication via non-coding RNAs is a recent discovery. Exogenous microRNAs (exog-miRNAs) mainly enter the host via the diet. Generally considered unstable in the gastrointestinal tract, some exogenous RNAs may resist these conditions, especially if transported in extracellular vesicles. They could then reach the intestines and more probably exert a regulatory effect. We give an overview of recent discoveries concerning dietary miRNAs, possible ways of enhancing their resistance to food processing and gut conditions, their transport in extracellular vesicles (animal- and plant-origin) and possible biological effects on recipient cells after ingestion. We critically focus on what we believe are the most relevant data for future pharmacological development of dietary miRNAs as therapeutic agents. Finally, we discuss the miRNA-mediated cross-kingdom regulation between diet, host and the gut microbiota. We conclude that, despite many obstacles and challenges, extracellular miRNAs are serious candidates to be targeted pharmacologically for development of new therapeutic agents.

ASCORBATE PEROXIDASE6 delays the onset of age-dependent leaf senescence

Age-dependent changes in reactive oxygen species (ROS) levels are critical in leaf senescence. While H2O2-reducing enzymes such as catalases and cytosolic ASCORBATE PEROXIDASE1 (APX1) tightly control the oxidative load during senescence, their regulation and function are not specific to senescence. Previously, we identified the role of ASCORBATE PEROXIDASE6 (APX6) during seed maturation in Arabidopsis (Arabidopsis thaliana). Here, we show that APX6 is a bona fide senescence-associated gene. APX6 expression is specifically induced in aging leaves and in response to senescence-promoting stimuli such as abscisic acid (ABA), extended darkness, and osmotic stress. apx6 mutants showed early developmental senescence and increased sensitivity to dark stress. Reduced APX activity, increased H2O2 level, and altered redox state of the ascorbate pool in mature pre-senescing green leaves of the apx6 mutants correlated with the early onset of senescence. Using transient expression assays in Nicotiana benthamiana leaves, we unraveled the age-dependent post-transcriptional regulation of APX6. We then identified the coding sequence of APX6 as a potential target of miR398, which is a key regulator of copper redistribution. Furthermore, we showed that mutants of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE7 (SPL7), the master regulator of copper homeostasis and miR398 expression, have a higher APX6 level compared with the wild type, which further increased under copper deficiency. Our study suggests that APX6 is a modulator of ROS/redox homeostasis and signaling in aging leaves that plays an important role in developmental- and stress-induced senescence programs.

Small RNA and degradome deep sequencing reveals important roles of microRNAs in cotton (Gossypium hirsutum L.) response to root-knot nematode Meloidogyne incognita infection

Investigation of cotton response to nematode infection will allow us to better understand the cotton immune defense mechanism and design a better biotechnological approach for efficiently managing pest nematodes in cotton. In this study, we firstly treated cotton by root knot nematode (RKN, Meloidogyne incognita) infections, then we employed the high throughput deep sequencing technology to sequence and genome-widely identify all miRNAs in cotton; finally, we analyzed the functions of these miRNAs in cotton response to RKN infections. A total of 266 miRNAs, including 193 known and 73 novel miRNAs, were identified by deep sequencing technology, which belong to 67 conserved and 66 novel miRNA families, respectively. A majority of identified miRNA families only contain one miRNA; however, miR482 family contains 14 members and some others contain 2-13 members. Certain miRNAs were specifically expressed in RKN-infected cotton roots and others were completely inhibited by RKN infection. A total of 50 miRNAs were differentially expressed after RKN infection, in which 28 miRNAs were up-regulated and 22 were inhibited by RKN treatment. Based on degradome sequencing, 87 gene targets were identified to be targeted by 57 miRNAs. These miRNA-targeted genes are involved in the interaction of cotton plants and nematode infection. Based on GO (gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, 466 genes from all 636 miRNA targets were mapped to 6340 GO terms, 181 genes from 228 targets of differentially expressed miRNAs were mapped to 1588 GO terms. The GO terms were then categorized into the three main GO classes: biological processes, cellular components, and molecular functions. The targets of differentially expressed miRNAs were enriched in 43 GO terms, including 22 biological processes, 10 cellular components, and 11 molecular functions (p < 0.05). Many identified processes were associated with organism responses to the environmental stresses, including regulation of nematode larval development, response to nematode, and response to flooding. Our results will enhance the study and application of developing new cotton cultivars for nematode resistance.

Identification of Small RNAs During High Light Acclimation in Arabidopsis thaliana

The biological significance of non-coding RNAs (ncRNAs) has been firmly established to be important for the regulation of genes involved in stress acclimation. Light plays an important role for the growth of plants providing the energy for photosynthesis; however, excessive light conditions can also cause substantial defects. Small RNAs (sRNAs) are a class of non-coding RNAs that regulate transcript levels of protein-coding genes and mediate epigenetic silencing. Next generation sequencing facilitates the identification of small non-coding RNA classes such as miRNAs (microRNAs) and small-interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), but changes in the ncRNA transcriptome in response to high light are poorly understood. We subjected Arabidopsis plants to high light conditions and performed a temporal in-depth study of the transcriptome data after 3 h, 6 h, and 2 days of high light treatment. We identified a large number of high light responsive miRNAs and sRNAs derived from NAT gene pairs, lncRNAs and TAS transcripts. We performed target predictions for differentially expressed miRNAs and correlated their expression levels through mRNA sequencing data. GO analysis of the targets revealed an overrepresentation of genes involved in transcriptional regulation. In A. thaliana, sRNA-mediated regulation of gene expression in response to high light treatment is mainly carried out by miRNAs and sRNAs derived from NAT gene pairs, and from lncRNAs. This study provides a deeper understanding of sRNA-dependent regulatory networks in high light acclimation.

Cca-miR398 increases copper sulfate stress sensitivity via the regulation of CSD mRNA transcription levels in transgenic Arabidopsis thaliana

MicroRNAs play crucial roles during the process of plant development under stress conditions. Copper is an essential micronutrient for most organisms and serves as an important redox-active cofactor for various functional proteins. In the present study, we investigated the effects of copper sulfate stress on hickory (Carya cathayensis) root development. We identified that hickory cca-miR398 was related to copper sulfate stress response, targeting Copper/Zinc superoxide dismutases (cytosolic (CSD1) and chloroplastic (CSD2)) and a 5b subunit of mitochondrial cytochrome C oxidase (COX5b.1) that are linked directly to stress regulatory networks. The sequence of hickory cca-miR398 is highly similar to that of Arabidopsis miR398b and miR398c, regardless of one nucleotide variation. Therefore, target genes of cca-miR398 were investigated by using 5′-Rapid-amplification of cDNA ends. An overexpression of cca-miR398 in Arabidopsis caused a reduction not only in root length and cotyledon greening, but also in the CSD1, CSD2, and CSD3 transcription levels. These reductions had greater significance in transgenic Arabidopsis than in wild-type Arabidopsis under copper sulfate stress. The level of physiological indicators also changed in transgenic Arabidopsis. In addition, the expressions of copper-responsive microRNAs, such as miR397 and miR408, were affected by the copper sulfate stress. These results showed that CSD possesses the ability to enhance copper sulfate stress response in both transgenic Arabidopsis and hickory roots by increasing the production of superoxide dismutase. Our results also demonstrated that cca-miR398 weakens hickory tolerance to copper sulfate by regulating CSD targets.

Computational tools for plant small RNA detection and categorization

Small RNAs (sRNAs) are important short-length molecules with regulatory functions essential for plant development and plasticity. High-throughput sequencing of total sRNA populations has revealed that the largest share of sRNA remains uncategorized. To better understand the role of sRNA-mediated cellular regulation, it is necessary to create accurate and comprehensive catalogues of sRNA and their sequence features, a task that currently relies on nontrivial bioinformatic approaches. Although a large number of computational tools have been developed to predict features of sRNA sequences, these tools are mostly dedicated to microRNAs and none integrates the functionalities necessary to describe units from all sRNA pathways thus far discovered in plants. Here, we review the different classes of sRNA found in plants and describe available bioinformatics tools that can help in their detection and categorization.

PAREameters: a tool for computational inference of plant miRNA-mRNA targeting rules using small RNA and degradome sequencing data

MicroRNAs (miRNAs) are short, non-coding RNAs that modulate the translation-rate of messenger RNAs (mRNAs) by directing the RNA-induced silencing complex to sequence-specific targets. In plants, this typically results in cleavage and subsequent degradation of the mRNA. Degradome sequencing is a high-throughput technique developed to capture cleaved mRNA fragments and thus can be used to support miRNA target prediction. The current criteria used for miRNA target prediction were inferred on a limited number of experimentally validated A. thaliana interactions and were adapted to fit these specific interactions; thus, these fixed criteria may not be optimal across all datasets (organisms, tissues or treatments). We present a new tool, PAREameters, for inferring targeting criteria from small RNA and degradome sequencing datasets. We evaluate its performance using a more extensive set of experimentally validated interactions in multiple A. thaliana datasets. We also perform comprehensive analyses to highlight and quantify the differences between subsets of miRNA-mRNA interactions in model and non-model organisms. Our results show increased sensitivity in A. thaliana when using the PAREameters inferred criteria and that using data-driven criteria enables the identification of additional interactions that further our understanding of the RNA silencing pathway in both model and non-model organisms.

Regulation of micoRNA2111 and its target IbFBK in sweet potato on wounding

Plant microRNAs (miRNAs) are non-coding RNAs, which are composed of 20-24 nucleotides. MiRNAs play important roles in plant growth and responses to biotic and abiotic stresses. Wounding is one of the most serious stresses for plants; however, the regulation of miRNAs in plants upon wounding is not well studied. In this study, miR2111, a wound-repressed miRNA, identified previously in sweet potato (Ipomoea batatas cv Tainung 57) by small RNA deep sequencing was chosen for further analysis. Based on sweet potato transcriptome database, F-box/kelch repeat protein (IbFBK), a target gene of miR2111, was identified. IbFBK is a wound-inducible gene, and the miR2111-induced cleavage site in IbFBK mRNA is between the 10th and 11th nucleotides of miR2111. IbFBK is a component of the E3 ligase SCF (SKP1-Cullin-F-box) complex participating in protein ubiquitination and degradation. The results of yeast two-hybrid and bimolecular fluorescence complementation assays demonstrate that IbFBK was conjugated with IbSKP1 through the F-box domain in IbFBK N-terminus to form SCF complex, and interacted with IbCNR8 through the kelch-repeat domain in IbFBK C-terminus. The interaction of IbFBK and IbCNR8 may lead to the ubiquitination and degradation of IbCNR8. In conclusion, the suppression of miR2111 resulted in the increase of IbFBK, and may regulate protein degradation of IbCNR8 in sweet potato responding to wounding.

Small RNA Bidirectional Crosstalk During the Interaction Between Wheat and Zymoseptoria tritici

Cross-kingdom RNA interference (RNAi) has been shown to play important roles during plant-pathogen interactions, and both plants and pathogens can use small RNAs (sRNAs) to silence genes in each other. This bidirectional cross-kingdom RNAi was still unexplored in the wheat-Zymoseptoria tritici pathosystem. Here, we performed a detailed analysis of the sRNA bidirectional crosstalk between wheat and Z. tritici. Using a combination of small RNA sequencing (sRNA-seq) and microRNA sequencing (mRNA-seq), we were able to identify known and novel sRNAs and study their expression and their action on putative targets in both wheat and Z. tritici. We predicted the target genes of all the sRNAs in either wheat or Z. tritici transcriptome and used degradome analysis to validate the cleavage of these gene transcripts. We could not find any clear evidence of a cross-kingdom RNAi acting by mRNA cleavage in this pathosystem. We also found that the fungal sRNA enrichment was lower in planta than during in vitro growth, probably due to the lower expression of the only Dicer gene of the fungus during plant infection. Our results support the recent finding that Z. tritici sRNAs cannot play important roles during wheat infection. However, we also found that the fungal infection induced wheat sRNAs regulating the expression of specific wheat genes, including auxin-related genes, as an immune response. These results indicate a role of sRNAs in the regulation of wheat defenses during Z. tritici infection. Our findings contribute to improve our understanding of the interactions between wheat and Z. tritici.

Identification and expression analysis of conserved microRNAs during short and prolonged chromium stress in rice (Oryza sativa)

MicroRNAs (miRNAs) are one of the most critical epigenetic regulators of gene expression which modulate a spectrum of development and defence response processes in plants. Chromium (Cr) contamination in rice imposes a serious concern to human health as rice is used as staple food throughout the world. Although several studies have established the differential response of miRNAs in rice during heavy metal (arsenic, cadmium) and heat or cold stress, no report is available about the response of miRNAs during Cr stress. In the present study, we identified 512 and 568 known miRNAs from Cr treated and untreated samples, respectively. Expression analysis revealed that 13 conserved miRNAs (miR156, miR159, miR160, miR166, miR169, miR171, miR396, miR397, miR408, miR444, miR1883, miR2877, miR5072) depicted preferential up- or down-regulation (> 4-fold change; P value < 0.05). Target gene prediction of differentially expressed miRNAs and their functional annotation suggested the important role of miRNAs in defence and detoxification of Cr though ATP-binding cassette transporters (ABC transporters), transcription factors, heat shock proteins, auxin response, and metal ion transport. Real-time PCR analysis validated the differential expression of selected miRNAs and their putative target genes. In conclusion, our study identifies and predicts miRNA-mediated regulation of signalling pathway in rice during Cr stress.

Up-regulation of microRNA targets correlates with symptom severity in Citrus sinensis plants infected with two different isolates of citrus psorosis virus

MAIN CONCLUSION: miRNA targets from Citrus sinensis are predicted and validated using degradome data. They show an up-regulation upon infection with CPsV, with a positive correlation between target expression and symptom severity. Sweet orange (Citrus sinensis) may suffer from disease symptoms induced by virus infections, thus resulting in drastic economic losses. Infection of sweet orange plants with two isolates of citrus psorosis virus (CPsV), expressing different symptomatologies, alters the accumulation of a set of endogenous microRNAs (miRNAs). Here, we predicted ten putative targets from four down-regulated miRNAs: three belonging to the CCAAT-binding transcription factor family (CBFAs); an Ethylene-responsive transcription factor (RAP2-7); an Integrase-type DNA-binding superfamily protein (AP2B); Transport inhibitor response 1 (TIR1); GRR1-like protein 1-related (GRR1); Argonaute 2-related (AGO2), Argonaute 7 (AGO7), and a long non-coding RNA (ncRNA). We validated six of them through analysis of leaf degradome data. Expressions of the validated targets increase in infected samples compared to healthy tissue, showing a more striking up-regulation those samples with higher symptom severity. This study contributes to the understanding of the miRNA-mediated regulation of important transcripts in Citrus sinensis through target validation and shed light in the manner a virus can alter host regulatory mechanisms leading to symptom expression.

WPMIAS: Whole-degradome-based Plant MicroRNA-Target Interaction Analysis Server

A critical aspect for exploring the biological function of a microRNA (miRNA) lies on exact detection and validation of its target mRNAs. However, no convenient and efficient web-based server is available for plant biologists to identify the experimentally verified target mRNAs of miRNAs. In this work, we built a comprehensive web-based platform for miRNA-target analysis, named as Whole-degradome-based Plant MiRNA-target Interaction Analysis Server (WPMIAS), for validation of predicted interactions of miRNAs and their target mRNAs (MTIs) by user-submitted data or all available pre-loaded degradome data. Besides, the server can construct degradome-based miRNA regulatory networks (MRNs) based on the validated MTIs to help study the functions and relations among miRNAs and target mRNAs. WPMIAS is also suitable for other small RNAs (sRNAs), such as 21-nt phased siRNAs (phasiRNAs) and natural antisense siRNAs (nat-siRNAs), which direct cleavage of target mRNAs. Currently, WPMIAS supports 64 plant species with ∼200 cDNA libraries and 274 pre-loaded plant degradome datasets. The user can identify all validated MTIs by analyzing all degradome data at a time and understand when and where MTIs take place and their cleavage levels. With the data obtained from WPMIAS, the user can build a plant miRNA-target map, where it is convenient to find interesting research ideas on miRNAs. In summary, WPMIAS is able to support a comprehensive web-based plant miRNA-target analysis and expected to greatly promote future research on plant miRNAs.

AVAILABILITY

It can be freely accessed at https://cbi.njau.edu.cn/WPMIAS/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

psRNATarget: a plant small RNA target analysis server (2017 release)

Plant regulatory small RNAs (sRNAs), which include most microRNAs (miRNAs) and a subset of small interfering RNAs (siRNAs), such as the phased siRNAs (phasiRNAs), play important roles in regulating gene expression. Although generated from genetically distinct biogenesis pathways, these regulatory sRNAs share the same mechanisms for post-translational gene silencing and translational inhibition. psRNATarget was developed to identify plant sRNA targets by (i) analyzing complementary matching between the sRNA sequence and target mRNA sequence using a predefined scoring schema and (ii) by evaluating target site accessibility. This update enhances its analytical performance by developing a new scoring schema that is capable of discovering miRNA–mRNA interactions at higher ‘recall rates’ without significantly increasing total prediction output. The scoring procedure is customizable for the users to search both canonical and non-canonical targets. This update also enables transmitting and analyzing ‘big’ data empowered by (a) the implementation of multi-threading chunked file uploading, which can be paused and resumed, using HTML5 APIs and (b) the allocation of significantly more computing nodes to its back-end Linux cluster. The updated psRNATarget server has clear, compelling and user-friendly interfaces that enhance user experiences and present data clearly and concisely. The psRNATarget is freely available at http://plantgrn.noble.org/psRNATarget/.

PAREsnip2: a tool for high-throughput prediction of small RNA targets from degradome sequencing data using configurable targeting rules

Small RNAs (sRNAs) are short, non-coding RNAs that play critical roles in many important biological pathways. They suppress the translation of messenger RNAs (mRNAs) by directing the RNA-induced silencing complex to their sequence-specific mRNA target(s). In plants, this typically results in mRNA cleavage and subsequent degradation of the mRNA. The resulting mRNA fragments, or degradome, provide evidence for these interactions, and thus degradome analysis has become an important tool for sRNA target prediction. Even so, with the continuing advances in sequencing technologies, not only are larger and more complex genomes being sequenced, but also degradome and associated datasets are growing both in number and read count. As a result, existing degradome analysis tools are unable to process the volume of data being produced without imposing huge resource and time requirements. Moreover, these tools use stringent, non-configurable targeting rules, which reduces their flexibility. Here, we present a new and user configurable software tool for degradome analysis, which employs a novel search algorithm and sequence encoding technique to reduce the search space during analysis. The tool significantly reduces the time and resources required to perform degradome analysis, in some cases providing more than two orders of magnitude speed-up over current methods.

Genome-wide comparative analysis in Solanaceous species reveals evolution of microRNAs targeting defense genes in Capsicum spp

MicroRNAs (miRNAs) play roles in various biological processes in plants including growth, development, and disease resistance. Previous studies revealed that some plant miRNAs produce secondary small interfering RNAs (siRNAs) such as phased, secondary siRNAs (phasiRNAs), and they regulate a cascade of gene expression. We performed a genome-wide comparative analysis of miRNAs in Solanaceous species (pepper, tomato, and potato), from an evolutionary perspective. Microsynteny of miRNAs was analysed based on the genomic loci and their flanking genes and most of the well-conserved miRNA genes maintained microsynteny in Solanaceae. We identified target genes of the miRNAs via degradome analysis and found that several miRNAs target many genes encoding nucleotide-binding leucine-rich repeat (NLR) or receptor-like proteins (RLPs), which are known to be major players in defense responses. In addition, disease-resistance-associated miRNAs trigger phasiRNA production in pepper, indicating amplification of the regulation of disease-resistance gene families. Among these, miR-n033a-3p, whose target NLRs have been duplicated in pepper, targets more NLRs belonging to specific subgroup in pepper than those in potato. miRNAs targeting resistance genes might have evolved to regulate numerous targets in Solanaceae, following expansion of target resistance genes. This study provides an insight into evolutionary relationship between miRNAs and their target defense genes in plants.

Fast-Forward Identification of Highly Effective Artificial Small RNAs Against Different Tomato spotted wilt virus Isolates

Artificial small RNAs (sRNAs), including artificial microRNAs (amiRNAs) and synthetic trans-acting small interfering RNAs (syn-tasiRNAs), are used to silence viral RNAs and confer antiviral resistance in plants. Here, the combined use of recent high-throughput methods for generating artificial sRNA constructs and the Tomato spotted wilt virus (TSWV)-Nicotiana benthamiana pathosystem allowed for the simple and rapid identification of amiRNAs with high anti-TSWV activity. A comparative analysis between the most effective amiRNA construct and a syn-tasiRNA construct including the four most effective amiRNA sequences showed that both were highly effective against two different TSWV isolates. These results highlight the usefulness of this high-throughput methodology for the fast-forward identification of artificial sRNAs with high antiviral activity prior to time-consuming generation of stably transformed plants.

Evaluation of mature miR398 family, expression analysis and the post-transcriptional regulation evidence in gamma-irradiated and nitrogen-stressed Medicago sativa seedlings

BACKGROUND

The miR398 (microRNA398) posttranscriptionally regulates the superoxide dismutase (Cu/Zn-SOD).

MATERIALS AND METHODS

The miR398 level was quantitated in gamma-irradiated (20Gy and 200Gy) and nitrogen-stressed seedlings by stem-loop RT-PCR (Real Time-Polymerase chain reaction). The positional preference of nucleotide (nt) for miR398 families and new targets was done. The SOD enzyme was assayed in native PAGE (Polyacrylamide gel electrophoresis).

RESULTS

A relative increase in miR398-3p expression in the roots, and reduction in shoots for gamma-irradiated tissue and downregulation in miR398-5p in shoots is noted. The nitrogen stress shows upregulation of miR398-3p in roots and shoots, whereas the expression of miR398-5p is upregulated in roots and downregulated in shoots. Positional preference in miR398-3p for 1-14 nt is 90% conserved unlike miR398-5p where no nucleotide positional preference is seen. Targets obtained were functionally characterized.

CONCLUSIONS

The mature miR398-5p and miR398-3p levels estimated in gamma-irradiated (20 and 200Gy) and nitrogen-stressed Medicago sativa seedlings show differential levels in roots and shoots. Native PAGE of Cu/Zn-SOD suggests its negative correlation with miR398 in shoots of irradiated and nitrogen-stressed samples. The nucleotide preferences for the nucleotide for a given position and functional characterization of targets are reported herein.

Identification and characterisation of microRNAs and their target genes in phosphate-starved Nicotiana benthamiana by small RNA deep sequencing and 5'RACE analysis

BACKGROUND

Phosphorus is an important macronutrient that is severely lacking in soils. In plants, specific microRNAs (miRNAs) essential for nutrient management and the regulation of stress responses are responsible for the control of many phosphate starvation responses. Further understanding of conserved and species-specific microRNA species has potential implications for the development of crops tolerant to soils with low phosphate.

RESULTS

This study identified and characterised phosphate starvation-responsive miRNAs in the native Australian tobacco Nicotiana benthamiana. Small RNA libraries were constructed and sequenced from phosphate-starved plant leaves, stems and roots. Twenty-four conserved miRNA families and 36 species-specific miRNAs were identified. The majority of highly phosphate starvation-responsive miRNAs were highly conserved, comprising of members from the miR399, miR827, and miR2111 families. In addition, two miRNA-star species were identified to be phosphate starvation-responsive. A total of seven miRNA targets were confirmed using RLM-5'RACE to be cleaved by five miRNA families, including two confirmed cleavage targets for Nbe-miR399 species, one for Nbe-miR2111, and two for Nbe-miR398. A number of N. benthamiana-specific features for conserved miRNAs were identified, including species-specific miRNA targets predicted or confirmed for miR399, miR827, and miR398.

CONCLUSIONS

Our results give an insight into the phosphate starvation-responsive miRNAs of Nicotiana benthamiana, and indicate that the phosphate starvation response pathways in N. benthamiana contain both highly conserved and species-specific components.

In Silico Identification of Conserved MiRNAs from Physcomitrella patens ESTs and their Target Characterization

Background:

MicroRNAs (miRNAs) are groups of small non-protein-coding endogenous single stranded RNAs with approximately 18-24 nucleotides in length. High evolutionary sequence conservation of miRNAs among plant species and availability of powerful computational tools allow identification of new orthologs and paralogs.

Methods:

New conserved miRNAs in P. patens were found by EST-based homology search approaches. All candidates were screened according to a series of miRNA filtering criteria. Unigene, DFCI Gene Index (PpspGI) databases and psRNATarget algorithm were applied to identify target transcripts using P. patens putative conserved miRNA sequences.

Results:

Nineteen conserved P. patens miRNAs were identified. The sequences were homologous to known reference plant mature miRNA from 10 miRNA families. They could be folded into the typical miRNA secondary structures. RepeatMasker algorithm demonstrated that ppt-miR2919e and pptmiR1533 had simple sequence repeats in their sequences. Target sites (49 genes) were identified for 7 out of 19 miRNAs. GO and KEGG analysis of targets indicated the involvement of some in important multiple biological and metabolic processes.

Conclusion:

The majority of the registered miRNAs in databases were predicted by computational approaches while many more have remained unknown. Due to the conserved nature of miRNAs in plant species from closely to distantly related, homology search-based approaches between plants species could lead to the identification of novel miRNAs in other plant species providing baseline information for further search about the biological functions and evolution of miRNAs.

Ambient temperature regulates the expression of a small set of sRNAs influencing plant development through NF-YA2 and YUC2

Plants substantially alter their developmental programme upon changes in the ambient temperature. The 21-24 nt small RNAs (sRNAs) are important gene expression regulators, which play a major role in development and adaptation. However, little is known about how the different sRNA classes respond to changes in the ambient temperature. We profiled the sRNA populations in four different tissues of Arabidopsis thaliana plants grown at 15°C, 21°C, and 27°C. We found that only a small fraction (0.6%) of the sRNA loci are ambient temperature-controlled. We identified thermoresponsive microRNAs and identified their target genes using degradome libraries. We verified that the target of the thermoregulated miR169, NF-YA2, is also ambient temperature-regulated. NF-YA2, as the component of the conserved transcriptional regulator NF-Y complex, binds the promoter of the flowering time regulator FT and the auxin biosynthesis gene YUC2. Other differentially expressed loci include thermoresponsive phased siRNA loci that target various auxin pathway genes and tRNA fragments. Furthermore, a temperature-dependent 24-nt heterochromatic siRNA locus in the promoter of YUC2 may contribute to the epigenetic regulation of auxin homeostasis. This holistic approach facilitated a better understanding of the role of different sRNA classes in ambient temperature adaptation of plants.

MicroRNAs from plants to animals, do they define a new messenger for communication?

MicroRNAs (miRNAs), a class of single-stranded non-coding RNA of about 22 nucleotides, are potent regulators of gene expression existing in both plants and animals. Recent studies showed that plant miRNAs could enter mammalian bloodstream via gastrointestinal tract, through which access a variety of tissues and cells of recipients to exert therapeutic effects. This intriguing phenomenon indicates that miRNAs of diet/plant origin may act as a new class of bioactive ingredients communicating with mammalian systems. In this review, in order to pinpoint the reason underlying discrepancies of miRNAs transmission from diet/plant to animals, the pathways that generate miRNAs and machineries involved in the functions of miRNAs in both kingdoms were outlined and compared. Then, the current controversies concerning cross-kingdom regulations and the potential mechanisms responsible for absorption and transfer of diet/plant-derived miRNAs were interpreted. Furthermore, the hormone-like action of miRNAs and the intricate interplay between miRNAs and hormones were implicated. Finally, how these findings may impact nutrition and medicine were briefly discussed.

Transcriptomic analyses of rice (Oryza sativa) genes and non-coding RNAs under nitrogen starvation using multiple omics technologies

BACKGROUND

Nitrogen (N) is a key macronutrient essential for plant growth, and its availability has a strong influence on crop development. The application of synthetic N fertilizers on crops has increased substantially in recent decades; however, the applied N is not fully utilized due to the low N use efficiency of crops. To overcome this limitation, it is important to understand the genome-wide responses and functions of key genes and potential regulatory factors in N metabolism.

RESULTS

Here, we characterized changes in the rice (Oryza sativa) transcriptome, including genes, newly identified putative long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) and their target mRNAs in response to N starvation using four different transcriptome approaches. Analysis of rice genes involved in N metabolism and/or transport using strand-specific RNA-Seq identified 2588 novel putative lncRNA encoding loci. Analysis of previously published RNA-Seq datasets revealed a group of N starvation-responsive lncRNAs showing differential expression under other abiotic stress conditions. Poly A-primed sequencing (2P-Seq) revealed alternatively polyadenylated isoforms of N starvation-responsive lncRNAs and provided precise 3' end information on the transcript models of these lncRNAs. Analysis of small RNA-Seq data identified N starvation-responsive miRNAs and down-regulation of miR169 family members, causing de-repression of NF-YA, as confirmed by strand-specific RNA-Seq and qRT-PCR. Moreover, we profiled the N starvation-responsive down-regulation of root-specific miRNA, osa-miR444a.4-3p, and Degradome sequencing confirmed MADS25 as a novel target gene.

CONCLUSIONS

In this study, we used a combination of multiple RNA-Seq analyses to extensively profile the expression of genes, newly identified lncRNAs, and microRNAs in N-starved rice roots and shoots. Data generated in this study provide an in-depth understanding of the regulatory pathways modulated by N starvation-responsive miRNAs. The results of comprehensive, large-scale data analysis provide valuable information on multiple aspects of the rice transcriptome, which may be useful in understanding the responses of rice plants to changes in the N supply status of soil.

Integration of proteomic and transcriptomic profiles reveals multiple levels of genetic regulation of salt tolerance in cotton

BACKGROUND

Salinity is a major abiotic stress that limits upland cotton growth and reduces fibre production worldwide. To reveal genetic regulation via transcript and protein levels after salt stress, we comprehensively analysed the global changes in mRNA, miRNA, and protein profiles in response to salt stress in two contrasting salt-tolerant cotton genotypes.

RESULTS

In the current study, proteomic and mRNA-seq data were combined to reveal that some genes are differentially expressed at both the proteomic and mRNA levels. However, we observed no significant change in mRNA corresponding to most of the strongly differentially abundant proteins. This finding may have resulted from global changes in alternative splicing events and miRNA levels under salt stress conditions. Evidence was provided indicating that several salt stress-responsive proteins can alter miRNAs and modulate alternative splicing events in upland cotton. The results of the stringent screening of the mRNA-seq and proteomic data between the salt-tolerant and salt-sensitive genotypes identified 63 and 85 candidate genes/proteins related to salt tolerance after 4 and 24 h of salt stress, respectively, between the tolerant and sensitive genotype. Finally, we predicted an interaction network comprising 158 genes/proteins and then discovered that two main clusters in the network were composed of ATP synthase (CotAD_74681) and cytochrome oxidase (CotAD_46197) in mitochondria. The results revealed that mitochondria, as important organelles involved in energy metabolism, play an essential role in the synthesis of resistance proteins during the process of salt exposure.

CONCLUSION

We provided a plausible schematic for the systematic salt tolerance model; this schematic reveals multiple levels of gene regulation in response to salt stress in cotton and provides a list of salt tolerance-related genes/proteins. The information here will facilitate candidate gene discovery and molecular marker development for salt tolerance breeding in cotton.

Tamarix microRNA Profiling Reveals New Insight into Salt Tolerance

The halophyte tamarisk (Tamarix) is extremely salt tolerant, making it an ideal material for salt tolerance-related studies. Although many salt-responsive genes of Tamarix were identified in previous studies, there are no reports on the role of post-transcriptional regulation in its salt tolerance. We constructed six small RNA libraries of Tamarix chinensis roots with NaCl treatments. High-throughput sequencing of the six libraries was performed and microRNA expression profiles were constructed. We investigated salt-responsive microRNAs to uncover the microRNA-mediated genes regulation. From these analyses, 251 conserved and 18 novel microRNA were identified from all small RNAs. From 191 differentially expressed microRNAs, 74 co-expressed microRNAs were identified as salt-responsive candidate microRNAs. The most enriched GO (gene ontology) terms for the 157 genes targeted by differentially expressed microRNAs suggested that transcriptions factors were highly active. Two hub microRNAs (miR414, miR5658), which connected by several target genes into an organic microRNA regulatory network, appeared to be the key regulators of post-transcriptional salt-stress responses. As the first survey on the tamarisk small RNAome, this study improves the understanding of tamarisk salt-tolerance mechanisms and will contribute to the molecular-assisted resistance breeding.

Wheat miR9678 Affects Seed Germination by Generating Phased siRNAs and Modulating Abscisic Acid/Gibberellin Signaling

Seed germination is important for grain yield and quality and rapid, near-simultaneous germination helps in cultivation; however, cultivars that germinate too readily can undergo preharvest sprouting (PHS), which causes substantial losses in areas that tend to get rain around harvest time. Moreover, our knowledge of mechanisms regulating seed germination in wheat (Triticum aestivum) remains limited. In this study, we analyzed function of a wheat-specific microRNA 9678 (miR9678), which is specifically expressed in the scutellum of developing and germinating seeds. Overexpression of miR9678 delayed germination and improved resistance to PHS in wheat through reducing bioactive gibberellin (GA) levels; miR9678 silencing enhanced germination rates. We provide evidence that miR9678 targets a long noncoding RNA (WSGAR) and triggers the generation of phased small interfering RNAs that play a role in the delay of seed germination. Finally, we found that abscisic acid (ABA) signaling proteins bind the promoter of miR9678 precursor and activate its expression, indicating that miR9678 affects germination by modulating the GA/ABA signaling.

IIKmTA: Inter and Intra Kingdom miRNA-Target Analyzer

BACKGROUND

Growing evidences suggest that microRNAs (miRNAs) can efficiently regulate gene expression at intracellular and extracellular levels. It has been previously reported that plant/food-derived miRNAs are highly enriched in human serum or serum from phytophagous animals, and they are responsible for regulating mammalian gene expression. Thus, miRNAs could function as active signaling molecules, which carry information across distinct species or even kingdoms. However, the mode of miRNA shuttling among various organisms is still a mystery to unravel. The intra and inter kingdom miRNA transfer has boosted up the hypothesis about the potential impact of plant or animal miRNAs on each other. To our knowledge, the software for analyzing cross-kingdom miRNA-targets is lacking.

RESULTS

We have developed a web-tool "IIKmTA: Inter and Intra Kingdom miRNA-Target Analyzer" utilizing a database; the data of which have been collected from another web server. Here, user can analyze the targeting potential of (i) plant miRNAs on animal UTRs (Untranslated regions), and vice versa (i.e., inter kingdom), (ii) plant miRNAs on plant UTRs and animal miRNAs on animal UTRs (i.e., intra kingdom). Further, user can analyze (i) miRNAs to targets, (ii) targets to miRNAs, and (iii) miRNA sets targeting sets of targets. For a wide variety of animal and plant species, IIKmTA can identify the miRNA binding sites in the probable target UTRs. Moreover, GC% and AU% of miRNAs will be calculated. All the results can be saved as .csv file.

CONCLUSIONS

Recent researches identified miRNAs in plants and human secretions and their role in regulating the human genes. Such findings indicate the therapeutic role of secretory miRNAs of such plants which exhibits medicinal value and in near future many diseases may be treated by consumption of these plant miRNAs through food. Using our newly developed database and analyzing tool, one can easily determine the different relationships between miRNAs and their targets across kingdoms. IIKmTA is freely available at http://www.bioinformatics.org/iikmta/ .

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.

Comparative Analysis of Fruit Ripening-Related miRNAs and Their Targets in Blueberry Using Small RNA and Degradome Sequencing

MicroRNAs (miRNAs) play vital roles in the regulation of fruit development and ripening. Blueberry is an important small berry fruit crop with economical and nutritional value. However, nothing is known about the miRNAs and their targets involved in blueberry fruit ripening. In this study, using high-throughput sequencing of small RNAs, 84 known miRNAs belonging to 28 families and 16 novel miRNAs were identified in white fruit (WF) and blue fruit (BF) libraries, which represent fruit ripening onset and in progress, respectively. Among them, 41 miRNAs were shown to be differentially expressed during fruit maturation, and 16 miRNAs representing 16 families were further chosen to validate the sRNA sequencing data by stem-loop qRT-PCR. Meanwhile, 178 targets were identified for 41 known and 7 novel miRNAs in WF and BF libraries using degradome sequencing, and targets of miR160 were validated using RLM-RACE (RNA Ligase-Mediated (RLM)-Rapid Amplification of cDNA Ends) approach. Moreover, the expression patterns of 6 miRNAs and their targets were examined during fruit development and ripening. Finally, integrative analysis of miRNAs and their targets revealed a complex miRNA-mRNA regulatory network involving a wide variety of biological processes. The findings will facilitate future investigations of the miRNA-mediated mechanisms that regulate fruit development and ripening in blueberry.

Identification and expression profiling of microRNAs involved in the stigma exsertion under high-temperature stress in tomato

Background

Autogamy in cultivated tomato varieties is a derived trait from wild type tomato plants, which are mostly allogamous. However, environmental stresses can cause morphological defects in tomato flowers and hinder autogamy. Under elevated temperatures, tomato plants usually exhibit the phenotype of stigma exsertion, with severely hindered self-pollination and fruit setting, whereas the inherent mechanism of stigma exsertion have been hitherto unknown. Numerous small RNAs (sRNAs) have been shown to play significant roles in plant development and stress responses, however, none of them have been studied with respect to stamen and pistil development under high-temperature conditions. We investigated the associations between stigma exsertion and small RNAs using high-throughput sequencing technology and molecular biology approaches.

Results

Sixteen sRNA libraries of Micro-Tom were constructed from plants stamen and pistil samples and sequenced after 2 d and 12 d of exposure to heat stress, respectively, from which a total of 110 known and 84 novel miRNAs were identified. Under heat stress conditions, 34 known and 35 novel miRNAs were differentially expressed in stamens, and 20 known and 10 novel miRNAs were differentially expressed in pistils. GO and KEGG pathway analysis showed that the predicted target genes of differentially expressed miRNAs were significantly enriched in metabolic pathways in both stamen and pistil libraries. Potential miRNA-target cleavage cascades that correlated with the regulation of stigma exsertion under heat stress conditions were found and validated through qRT-PCR and RLM-5′ RACE.

Conclusion

Overall, a global spectrum of known and novel miRNAs involved in tomato stigma exsertion and induced by high temperatures were identified using high-throughput sequencing and molecular biology approaches, laying a foundation for revealing the miRNA-mediated regulatory network involved in the development of tomato stamens and pistils under high-temperature conditions.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4238-9) contains supplementary material, which is available to authorized users.

Global Analysis of Small RNA Dynamics during Seed Development of Picea glauca and Arabidopsis thaliana Populations Reveals Insights on their Evolutionary Trajectories

While DNA methylation carries genetic signals and is instrumental in the evolution of organismal complexity, small RNAs (sRNAs), ~18-24 ribonucleotide (nt) sequences, are crucial mediators of methylation as well as gene silencing. However, scant study deals with sRNA evolution via featuring their expression dynamics coupled with species of different evolutionary time. Here we report an atlas of sRNAs and microRNAs (miRNAs, single-stranded sRNAs) produced over time at seed-set of two major spermatophytes represented by populations of Picea glauca and Arabidopsis thaliana with different seed-set duration. We applied diverse profiling methods to examine sRNA and miRNA features, including size distribution, sequence conservation and reproduction-specific regulation, as well as to predict their putative targets. The top 27 most abundant miRNAs were highly overlapped between the two species (e.g., miR166,-319 and-396), but in P. glauca, they were less abundant and significantly less correlated with seed-set phases. The most abundant sRNAs in libraries were deeply conserved miRNAs in the plant kingdom for Arabidopsis but long sRNAs (24-nt) for P. glauca. We also found significant difference in normalized expression between populations for population-specific sRNAs but not for lineage-specific ones. Moreover, lineage-specific sRNAs were enriched in the 21-nt size class. This pattern is consistent in both species and alludes to a specific type of sRNAs (e.g., miRNA, tasiRNA) being selected for. In addition, we deemed 24 and 9 sRNAs in P. glauca and Arabidopsis, respectively, as sRNA candidates targeting known adaptive genes. Temperature had significant influence on selected gene and miRNA expression at seed development in both species. This study increases our integrated understanding of sRNA evolution and its potential link to genomic architecture (e.g., sRNA derivation from genome and sRNA-mediated genomic events) and organismal complexity (e.g., association between different sRNA expression and their functionality).

Unraveling multifaceted contributions of small regulatory RNAs to photomorphogenic development in Arabidopsis

Background

Post-transcriptional control of gene expression mediated by small regulatory RNAs (sRNAs) is vital for growth and development of diverse organisms. The biogenesis of sRNAs is regulated by both positive and negative regulators known to regulate photomorphogenic development. Two microRNAs (miRNAs), miR157 and miR319, also regulate photomorphogenesis. However, genome-wide profiling of sRNAs and their regulation of target genes during photomorphogenesis has been missing. We provide a comprehensive view of sRNA-controlled gene expression in this developmental process.

Results

By profiling sRNAs and the 5′ ends of degraded mRNAs during the first 24 h of photomorphogenic development in Arabidopsis, we identified 335 sRNA-mediated mRNA cleavage events in de-etiolating seedlings. These cleavage events are primarily resulted from actions of highly expressed miRNAs and irrelevant to the abundance of target mRNAs. In the light, the expression of the slicer protein gene ARGONAUTE1 in the miRNA functioning pathway could be fine-tuned by miRNA168a/b. We also found that miR396a/b positively regulates de-etiolation by suppressing GROWTH REGULATING FACTORs. Our results suggest that the miRNAs are required to tune down the target mRNAs and regulate photomorphogenesis.

Conclusion

sRNAs may have a broad impact on gene expression regulation for optimized photomorphogenic development. With both positive and negative regulators under the control of sRNAs, young Arabidopsis seedlings can have a timely but not exaggerated developmental adaptation to light.

Electronic supplementary material

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