Transcriptome-wide identification of microRNA targets in rice

miR408 is involved in abiotic stress responses in Arabidopsis

MicroRNAs (miRNAs) are small RNAs that regulate the expression of target genes post-transcriptionally; they are known to play major roles in development and responses to abiotic stress. miR408 is a highly conserved miRNA in plants that responds to the availability of copper and targets genes encoding copper-containing proteins. It was recently recognized to be an important component of the HY5-SPL7 gene network that mediates a coordinated response to light and copper, illustrating its central role in the response of plants to the environment. Expression of miR408 is significantly affected by a variety of developmental and ‏environmental conditions; however, its biological function is ‏unknown. Involvement of miR408 in the abiotic stress response was investigated in Arabidopsis. Expression of miR408, as well as its target genes, was investigated in response to salinity, cold, oxidative stress, drought and osmotic stress. Analyses of transgenic plants with modulated miR408 expression revealed that higher miR408 expression leads to improved tolerance to salinity, cold and oxidative stress, but enhanced sensitivity to drought and osmotic stress. Cellular antioxidant capacity was enhanced in plants with elevated miR408 expression, as manifested by reduced levels of reactive oxygen species and induced expression of genes associated with antioxidative functions, including Cu/Zn superoxide dismutases (CSD1 and CSD2) and glutathione-S-transferase (GST-U25), as well as auxiliary genes: the copper chaperone CCS1 and the redox stress-associated gene SAP12. Overall, the results demonstrate significant involvement of miR408 in abiotic stress responses, emphasizing the central function of miR408 in plant survival.

Identification of four functionally important microRNA families with contrasting differential expression profiles between drought-tolerant and susceptible rice leaf at vegetative stage

Background

Developing drought-tolerant rice varieties with higher yield under water stressed conditions provides a viable solution to serious yield-reduction impact of drought. Understanding the molecular regulation of this polygenic trait is crucial for the eventual success of rice molecular breeding programmes. microRNAs have received tremendous attention recently due to its importance in negative regulation. In plants, apart from regulating developmental and physiological processes, microRNAs have also been associated with different biotic and abiotic stresses. Hence here we chose to analyze the differential expression profiles of microRNAs in three drought treated rice varieties: Vandana (drought-tolerant), Aday Sel (drought-tolerant) and IR64 (drought-susceptible) in greenhouse conditions via high-throughput sequencing.

Results

Twenty-six novel microRNA candidates involved in the regulation of diverse biological processes were identified based on the detection of miRNA*. Out of their 110 predicted targets, we confirmed 16 targets from 5 novel microRNA candidates. In the differential expression analysis, mature microRNA members from 49 families of known Oryza sativa microRNA were differentially expressed in leaf and stem respectively with over 28 families having at least a similar mature microRNA member commonly found to be differentially expressed between both tissues. Via the sequence profiling data of leaf samples, we identified osa-miR397a/b, osa-miR398b, osa-miR408-5p and osa-miR528-5p as being down-regulated in two drought-tolerant rice varieties and up-regulated in the drought-susceptible variety. These microRNAs are known to be involved in regulating starch metabolism, antioxidant defence, respiration and photosynthesis. A wide range of biological processes were found to be regulated by the target genes of all the identified differentially expressed microRNAs between both tissues, namely root development (5.3–5.7 %), cell transport (13.2–18.4 %), response to stress (10.5–11.3 %), lignin catabolic process (3.8–5.3 %), metabolic processes (32.1–39.5 %), oxidation-reduction process (9.4–13.2 %) and DNA replication (5.7–7.9 %). The predicted target genes of osa-miR166e-3p, osa-miR166h-5p*, osa-miR169r-3p* and osa-miR397a/b were found to be annotated to several of the aforementioned biological processes.

Conclusions

The experimental design of this study, which features rice varieties with different drought tolerance and tissue specificity (leaf and stem), has provided new microRNA profiling information. The potentially regulatory importance of the microRNA genes mentioned above and their target genes would require further functional analyses.

Electronic supplementary material

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

Combined Small RNA and Degradome Sequencing Reveals Novel MiRNAs and Their Targets in the High-Yield Mutant Wheat Strain Yunong 3114

Wheat is one of the main food sources worldwide; large amount studies have been conducted to improve wheat production. MicroRNAs (miRNAs) with about 20-30 nucleotide are a class of regulatory small RNAs (sRNAs), which could regulate gene expression through sequence-specific base pairing with target mRNAs, playing important roles in plant growth. An ideal plant architecture (IPA) is crucial to enhance yield in bread wheat. In this study, the high-yield wheat strain Yunong 3114 was EMS-mutagenesis from the wild-type strain Yunong 201, exhibiting a preferable plant structure compared with the wild-type strain. We constructed small RNA and degradome libraries from Yunong 201 and Yunong 3114, and performed small RNA sequencing of these libraries in order identify miRNAs and their targets related to IPA in wheat. Totally, we identified 488 known and 837 novel miRNAs from Yunong 3114 and 391 known and 533 novel miRNAs from Yunong 201. The number of miRNAs in the mutant increased. A total of 37 known and 432 putative novel miRNAs were specifically expressed in the mutant strain; furthermore, 23 known and 159 putative novel miRNAs were specifically expressed in the wild-type strain. A total of 150 known and 100 novel miRNAs were differentially expressed between mutant and wild-type strains. Among these differentially expressed novel miRNAs, 4 and 8 predict novel miRNAs were evidenced by degradome sequencing and showed up-regulated and down-regulated expressions in the mutant strain Yunong 3114, respectively. Targeted gene annotation and previous results indicated that this set of miRNAs is related to plant structure. Our results further suggested that miRNAs may be necessary to obtain an optimal wheat structure.

MADS-box transcription factor OsMADS25 regulates root development through affection of nitrate accumulation in rice

MADS-box transcription factors are vital regulators participating in plant growth and development process and the functions of most of them are still unknown. ANR1 was reported to play a key role in controlling lateral root development through nitrate signal in Arabidopsis. OsMADS25 is one of five ANR1-like genes in Oryza Sativa and belongs to the ANR1 clade. Here we have investigated the role of OsMADS25 in the plant’s responses to external nitrate in Oryza Sativa. Our results showed that OsMADS25 protein was found in the nucleus as well as in the cytoplasm. Over-expression of OsMADS25 significantly promoted lateral and primary root growth as well as shoot growth in a nitrate-dependent manner in Arabidopsis. OsMADS25 overexpression in transgenic rice resulted in significantly increased primary root length, lateral root number, lateral root length and shoot fresh weight in the presence of nitrate. Down-regulation of OsMADS25 in transgenic rice exhibited significantly reduced shoot and root growth in the presence of nitrate. Furthermore, over-expression of OsMADS25 in transgenic rice promoted nitrate accumulation and significantly increased the expressions of nitrate transporter genes at high rates of nitrate supply while down-regulation of OsMADS25 produced the opposite effect. Taken together, our findings suggest that OsMADS25 is a positive regulator control lateral and primary root development in rice.

Growth-Regulating Factors (GRFs): A Small Transcription Factor Family with Important Functions in Plant Biology

Growth-regulating factors (GRFs) are plant-specific transcription factors that were originally identified for their roles in stem and leaf development, but recent studies highlight them to be similarly important for other central developmental processes including flower and seed formation, root development, and the coordination of growth processes under adverse environmental conditions. The expression of several GRFs is controlled by microRNA miR396, and the GRF-miRNA396 regulatory module appears to be central to several of these processes. In addition, transcription factors upstream of GRFs and miR396 have been discovered, and gradually downstream target genes of GRFs are being unraveled. Here, we review the current knowledge of the biological functions performed by GRFs and survey available molecular data to illustrate how they exert their roles at the cellular level.

Identification of miRNAs and Their Targets in Cotton Inoculated with Verticillium dahliae by High-Throughput Sequencing and Degradome Analysis

MicroRNAs (miRNAs) are a group of endogenous small non-coding RNAs that play important roles in plant growth, development, and stress response processes. Verticillium wilt is a vascular disease in plants mainly caused by Verticillium dahliae Kleb., the soil-borne fungal pathogen. However, the role of miRNAs in the regulation of Verticillium defense responses is mostly unknown. This study aimed to identify new miRNAs and their potential targets that are involved in the regulation of Verticillium defense responses. Four small RNA libraries and two degradome libraries from mock-infected and infected roots of cotton (both Gossypiumhirsutum L. and Gossypiumbarbadense L.) were constructed for deep sequencing. A total of 140 known miRNAs and 58 novel miRNAs were identified. Among the identified miRNAs, many were differentially expressed between libraries. Degradome analysis showed that a total of 83 and 24 genes were the targets of 31 known and 14 novel miRNA families, respectively. Gene Ontology analysis indicated that many of the identified miRNA targets may function in controlling root development and the regulation of Verticillium defense responses in cotton. Our findings provide an overview of potential miRNAs involved in the regulation of Verticillium defense responses in cotton and the interactions between miRNAs and their corresponding targets. The profiling of these miRNAs lays the foundation for further understanding of the function of small RNAs in regulating plant response to fungal infection and Verticillium wilt in particular.

Transcriptome-Wide Identification of miRNA Targets under Nitrogen Deficiency in Populus tomentosa Using Degradome Sequencing

miRNAs are endogenous non-coding small RNAs with important regulatory roles in stress responses. Nitrogen (N) is an indispensable macronutrient required for plant growth and development. Previous studies have identified a variety of known and novel miRNAs responsive to low N stress in plants, including Populus. However, miRNAs involved in the cleavage of target genes and the corresponding regulatory networks in response to N stress in Populus remain largely unknown. Consequently, degradome sequencing was employed for global detection and validation of N-responsive miRNAs and their targets. A total of 60 unique miRNAs (39 conserved, 13 non-conserved, and eight novel) were experimentally identified to target 64 mRNA transcripts and 21 precursors. Among them, we further verified the cleavage of 11 N-responsive miRNAs identified previously and provided empirical evidence for the cleavage mode of these miRNAs on their target mRNAs. Furthermore, five miRNA stars (miRNA*s) were shown to have cleavage function. The specificity and diversity of cleavage sites on the targets and miRNA precursors in P. tomentosa were further detected. Identification and annotation of miRNA-mediated cleavage of target genes in Populus can increase our understanding of miRNA-mediated molecular mechanisms of woody plants adapted to low N environments.

Unique miRNome during anthesis in drought-tolerant indica rice var. Nagina 22

Drought-tolerant rice variety, Nagina 22 (N22), has a unique spikelet miRNome during anthesis stage drought as well as transition from heading to anthesis. Molecular characterization of genetic diversity of rice is essential to understand the evolution and molecular basis of various agronomically important traits such as drought tolerance. miRNAs play an important role in regulating plant development as well as stress response such as drought. In this study, we characterized the yet unexplored dynamics of the spikelet miRNA population during developmental transition from 'heading' to 'anthesis' as well as anthesis stage drought stress in a drought-tolerant indica rice variety, N22. A significant proportion of miRNA population (~20 %) in N22 spikelets is modulated during transition from heading to anthesis indicating a unique miRNome at anthesis, a developmental stage highly sensitive to stress (drought/heat). Based on the analysis of degradome data, majority of differentially regulated miRNAs appear to regulate transcription factors, some of which are implicated in regulation of development and fertilization. Similarly, drought during anthesis leads to a global change in miRNA expression pattern including those which regulate ROS homeostasis. It was possible to identify several miRNAs that were not reported to be drought responsive in earlier studies. Interestingly, a significant proportion of the drought-regulated miRNAs co-localize within QTLs related to drought tolerance and associated traits. Comparison of the expression profiles between N22 and Pusa Basmati 1 (drought sensitive) identified miRNAs with variety-specific expression patterns during phase transition (miR164, miR396, miR812, and miR1881) as well as drought stress (miR1881) indicating an evolution of a distinct and variety-specific regulatory mechanism. The promoters of these miRNAs contain LREs (light-responsive elements) and are induced by dark treatment. It was also possible to identify 4 novel miRNAs including an intronic miRNA that was conserved in both rice varieties.

Identification and Characterization of ABA-Responsive MicroRNAs in Rice

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that silence genes through mRNA degradation or translational inhibition. The phytohormone abscisic acid (ABA) is essential for plant development and adaptation to abiotic and biotic stresses. In Arabidopsis, miRNAs are implicated in ABA functions. However, ABA-responsive miRNAs have not been systematically studied in rice. Here high throughput sequencing of small RNAs revealed that 107 miRNAs were differentially expressed in the rice ABA deficient mutant, Osaba1. Of these, 13 were confirmed by stem-loop RT-PCR. Among them, miR1425-5P, miR169a, miR169n, miR390-5P, miR397a and miR397b were up-regulated, but miR162b reduced in expression in Osaba1. The targets of these 13 miRNAs were predicted and validated by gene expression profiling. Interestingly, the expression levels of these miRNAs and their targets were regulated by ABA. Cleavage sites were detected on 7 of the miRNA targets by 5'-Rapid Amplification of cDNA Ends (5'-RACE). Finally, miR162b and its target OsTRE1 were shown to affect rice resistance to drought stress, suggesting that miR162b increases resistance to drought by targeting OsTRE1. Our work provides important information for further characterization and functional analysis of ABA-responsive miRNAs in rice.

Construction of regulatory networks mediated by small RNAs responsive to abiotic stresses in rice (Oryza sativa)

Plants have evolved exquisite molecular mechanisms to adapt to diverse abiotic stresses. MicroRNAs play an important role in stress response in plants. However, whether the other small RNAs (sRNAs) possess stress-related roles remains elusive. In this study, thousands of sRNAs responsive to cold, drought and salt stresses were identified in rice seedlings and panicles by using high-throughput sequencing data. These sRNAs were classified into 12 categories, including "Panicle_Cold_Down", "Panicle_Cold_Up", "Panicle_Drought_Down", "Panicle_Drought_Up", "Panicle_Salt_Down", "Panicle_Salt_Up", "Seedling_Cold_Down", "Seedling_Cold_Up", "Seedling_Drought_Down", "Seedling_Drought_Up", "Seedling_Salt_Down" and "Seedling_Salt_Up". The stress-responsive sRNAs enriched in Argonaute 1 were extracted for target prediction and degradome sequencing data-based validation, which enabled network construction. Within certain subnetworks, some target genes were further supported by microarray data. Literature mining indicated that certain targets were potentially involved in stress response. These results demonstrate that the established networks are biologically meaningful. We discovered that in some cases, one sRNA sequence could be assigned to two or more categories. Moreover, within certain target-centered subnetworks, one transcript was regulated by several stress-responsive sRNAs assigned to different categories. It implies that these subnetworks are potentially implicated in stress signal crosstalk. Together, our results could advance the current understanding of the biological role of plant sRNAs in stress signaling.

MicroRNA-directed cleavage of targets: mechanism and experimental approaches

MicroRNAs (miRNAs) are a large family of post-transcriptional regulators, which are 21-24 nt in length and play a role in a wide variety of biological processes in eukaryotes. The past few years have seen rapid progress in our understanding of miRNA biogenesis and the mechanism of action, which commonly entails a combination of target degradation and translational repression. The target degradation mediated by Argonaute-catalyzed endonucleolytic cleavage exerts a significant repressive effect on target mRNA expression, particularly during rapid developmental transitions. This review outlines the current understanding of the mechanistic aspects of this important process and discusses several different experimental approaches to identify miRNA cleavage targets. [BMB Reports 2014; 47(8): 417-423]

PlantMirnaT: miRNA and mRNA integrated analysis fully utilizing characteristics of plant sequencing data

miRNA is known to regulate up to several hundreds coding genes, thus the integrated analysis of miRNA and mRNA expression data is an important problem. Unfortunately, the integrated analysis is challenging since it needs to consider expression data of two different types, miRNA and mRNA, and target relationship between miRNA and mRNA is not clear, especially when microarray data is used. Fortunately, due to the low sequencing cost, small RNA and RNA sequencing are routinely processed and we may be able to infer regulation relationships between miRNAs and mRNAs more accurately by using sequencing data. However, no method is developed specifically for sequencing data. Thus we developed PlantMirnaT, a new miRNA-mRNA integrated analysis system. To fully leverage the power of sequencing data, three major features are developed and implemented in PlantMirnaT. First, we implemented a plant-specific short read mapping tool based on recent discoveries on miRNA target relationship in plant. Second, we designed and implemented an algorithm considering miRNA targets in the full intragenic region, not just 3' UTR. Lastly but most importantly, our algorithm is designed to consider quantity of miRNA expression and its distribution on target mRNAs. The new algorithm was used to characterize rice under drought condition using our proprietary data. Our algorithm successfully discovered that two miRNAs, miRNA1425-5p, miRNA 398b, that are involved in suppression of glucose pathway in a naturally drought resistant rice, Vandana. The system can be downloaded at https://sites.google.com/site/biohealthinformaticslab/resources.

Identification of novel drought-responsive microRNAs and trans-acting siRNAs from Sorghum bicolor (L.) Moench by high-throughput sequencing analysis

Small non-coding RNAs (sRNAs) namely microRNAs (miRNAs) and trans-acting small interfering RNAs (tasi-RNAs) play a crucial role in post-transcriptional regulation of gene expression and thus the control plant development and stress responses. In order to identify drought-responsive miRNAs and tasi-RNAs in sorghum, we constructed small RNA libraries from a drought tolerant (M35-1) and susceptible (C43) sorghum genotypes grown under control and drought stress conditions, and sequenced by Illumina Genome Analyzer IIx. Ninety seven conserved and 526 novel miRNAs representing 472 unique miRNA families were identified from sorghum. Ninety-six unique miRNAs were found to be regulated by drought stress, of which 32 were up- and 49 were down-regulated (fold change ≥ 2 or ≤ -2) at least in one genotype, while the remaining 15 miRNAs showed contrasting drought-regulated expression pattern between genotypes. A maximum of 17 and 18 miRNAs was differentially regulated under drought stress condition in the sensitive and tolerant genotypes, respectively. These results suggest that genotype dependent stress responsive regulation of miRNAs may contribute, at least in part, to the differential drought tolerance of sorghum genotypes. We also identified two miR390-directed TAS3 gene homologs and the auxin response factors as tasi-RNA targets. We predicted more than 1300 unique target genes for the novel and conserved miRNAs. These target genes were predicted to be involved in different cellular, metabolic, response to stimulus, biological regulation, and developmental processes. Genome-wide identification of stress-responsive miRNAs, tasi-RNAs and their targets identified in this study will be useful in unraveling the molecular mechanisms underlying drought stress responses and genetic improvement of biomass production and stress tolerance in sorghum.

High-throughput sequencing reveals miRNA effects on the primary and secondary production properties in long-term subcultured Taxus cells

Plant-cell culture technology is a promising alternative for production of high-value secondary metabolites but is limited by the decreased metabolite production after long-term subculture. The goal of this study was to determine the effects of miRNAs on altered gene expression profiles during long-term subculture. Two Taxus cell lines, CA (subcultured for 10 years) and NA (subcultured for 6 months), were high-throughput sequenced at the mRNA and miRNA levels. A total of 265 known (78.87% of 336) and 221 novel (79.78% of 277) miRNAs were differentially expressed. Furthermore, 67.17% of the known differentially expressed (DE) miRNAs (178) and 60.63% of the novel DE-miRNAs (134) were upregulated in NA. A total of 275 inverse-related miRNA/mRNA modules were identified by target prediction analysis. Functional annotation of the targets revealed that the high-ranking miRNA targets were those implicated in primary metabolism and abiotic or biotic signal transduction. For example, various genes for starch metabolism and oxidative phosphorylation were inversely related to the miRNA levels, thereby indicating that miRNAs have important roles in these pathways. Interestingly, only a few genes for secondary metabolism were inversely related to miRNA, thereby indicating that factors other than miRNA are present in the regulatory system. Moreover, miR8154 and miR5298b were upregulated miRNAs that targeted a mass of DE genes. The overexpression of these miRNAs in CA increased the genes of taxol, phenylpropanoid, and flavonoid biosynthesis, thereby suggesting their function as crucial factors that regulate the entire metabolic network during long-term subculture. Our current studies indicated that a positive conversion of production properties from secondary metabolism to primary metabolism occurred in long-term subcultured cells. miRNAs are important regulators in the upregulation of primary metabolism.

NPK macronutrients and microRNA homeostasis

Macronutrients are essential elements for plant growth and development. In natural, non-cultivated systems, the availability of macronutrients is not a limiting factor of growth, due to fast recycling mechanisms. However, their availability might be an issue in modern agricultural practices, since soil has been frequently over exploited. From a crop management perspective, the nitrogen (N), phosphorus (P), and potassium (K) are three important limiting factors and therefore frequently added as fertilizers. NPK are among the nutrients that have been reported to alter post-embryonic root developmental processes and consequently, impairs crop yield. To cope with nutrients scarcity, plants have evolved several mechanisms involved in metabolic, physiological, and developmental adaptations. In this scenario, microRNAs (miRNAs) have emerged as additional key regulators of nutrients uptake and assimilation. Some studies have demonstrated the intrinsic relation between miRNAs and their targets, and how they can modulate plants to deal with the NPK availability. In this review, we focus on miRNAs and their regulation of targets involved in NPK metabolism. In general, NPK starvation is related with miRNAs that are involved in root-architectural changes and uptake activity modulation. We further show that several miRNAs were discovered to be involved in plant-microbe symbiosis during N and P uptake, and in this way we present a global view of some studies that were conducted in the last years. The integration of current knowledge about miRNA-NPK signaling may help future studies to focus in good candidates genes for the development of important tools for plant nutritional breeding.

Exploration of microRNAs and their targets engaging in the resistance interaction between wheat and stripe rust

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most destructive diseases of wheat worldwide. miRNAs are important regulators, they play very central roles in plant organ development, vegetable phase change and defense responses. In this study, two miRNA libraries from wheat cultivar Xingzi 9104 (XZ) challenged with the avirulent Pst race CYR32 and sterile water were constructed, respectively. A total of 596 miRNA candidates were obtained. 420 wheat-specific candidate miRNAs were screened in adult plants challenged with Pst using microarray-based analyses. We analyzed the abundance of candidate miRNAs, and the levels of a subset of candidate miRNAs were determined by quantitative real time PCR (qRT-PCR). The qRT-PCR results indicated that some miRNAs were involved in the incompatible interaction between wheat and Pst. In addition, we identified some miRNAs differentially expressed in different leaves. Additionally, the target genes of wheat miRNAs were confirmed by using degradome sequencing technology. Most of the annotated target genes are related to signal transduction, energy metabolism, and other functions. We selected some target genes for relative expression analysis using qRT-PCR, and found that RabGAP/TBC domain-containing protein, zinc finger protein and Cysteine-rich receptor-like protein kinase 41 may play important role in the incompatible interaction between XZ and CYR32. Intriguingly, miRNAs and target gene seem to form a complicated regulation network that regulates the wheat-Pst interaction. Our data provide the foundation for evaluating the important regulatory roles of miRNAs in the wheat-Pst interaction.

Genome-wide screening and functional analysis identify a large number of long noncoding RNAs involved in the sexual reproduction of rice

Background

Long noncoding RNAs (lncRNAs) play important roles in a wide range of biological processes in mammals and plants. However, the systematic examination of lncRNAs in plants lags behind that in mammals. Recently, lncRNAs have been identified in Arabidopsis and wheat; however, no systematic screening of potential lncRNAs has been reported for the rice genome.

Results

In this study, we perform whole transcriptome strand-specific RNA sequencing (ssRNA-seq) of samples from rice anthers, pistils, and seeds 5 days after pollination and from shoots 14 days after germination. Using these data, together with 40 available rice RNA-seq datasets, we systematically analyze rice lncRNAs and definitively identify lncRNAs that are involved in the reproductive process. The results show that rice lncRNAs have some different characteristics compared to those of Arabidopsis and mammals and are expressed in a highly tissue-specific or stage-specific manner. We further verify the functions of a set of lncRNAs that are preferentially expressed in reproductive stages and identify several lncRNAs as competing endogenous RNAs (ceRNAs), which sequester miR160 or miR164 in a type of target mimicry. More importantly, one lncRNA, XLOC_057324, is demonstrated to play a role in panicle development and fertility. We also develop a source of rice lncRNA-associated insertional mutants.

Conclusions

Genome-wide screening and functional analysis enabled the identification of a set of lncRNAs that are involved in the sexual reproduction of rice. The results also provide a source of lncRNAs and associated insertional mutants in rice.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0512-1) contains supplementary material, which is available to authorized users.

Identification and characterization of cold-responsive microRNAs in tea plant (Camellia sinensis) and their targets using high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are approximately 19 ~ 21 nucleotide noncoding RNAs produced by Dicer-catalyzed excision from stem-loop precursors. Many plant miRNAs have critical functions in development, nutrient homeostasis, abiotic stress responses, and pathogen responses via interaction with specific target mRNAs. Camellia sinensis is one of the most important commercial beverage crops in the world. However, miRNAs associated with cold stress tolerance in C. sinensis remains unexplored. The use of high-throughput sequencing can provide a much deeper understanding of miRNAs. To obtain more insight into the function of miRNAs in cold stress tolerance, Illumina sequencing of C. sinensis sRNA was conducted.

RESULT

Solexa sequencing technology was used for high-throughput sequencing of the small RNA library from the cold treatment of tea leaves. To align the sequencing data with known plant miRNAs, we characterized 106 conserved C. sinensis miRNAs. In addition, 215 potential candidate miRNAs were found, among, which 98 candidates with star sequences were chosen as novel miRNAs. Both congruously and differentially regulated miRNAs were obtained, and cultivar-specific miRNAs were identified by microarray-based hybridization in response to cold stress. The results were also confirmed by quantitative real-time polymerase chain reaction. To confirm the targets of miRNAs, two degradome libraries from two treatments were constructed. According to degradome sequencing, 455 and 591 genes were identified as cleavage targets of miRNAs from cold treatments and control libraries, respectively, and 283 targets were present in both libraries. Functional analysis of these miRNA targets indicated their involvement in important activities, such as development, regulation of transcription, and stress response.

CONCLUSIONS

We discovered 31 up-regulated miRNAs and 43 down-regulated miRNAs in 'Yingshuang', and 46 up-regulated miRNA and 45 down-regulated miRNAs in 'Baiye 1' in response to cold stress, respectively. A total of 763 related target genes were detected by degradome sequencing. The RLM-5'RACE procedure was successfully used to map the cleavage sites in six target genes of C. sinensis. These findings reveal important information about the regulatory mechanism of miRNAs in C. sinensis, and promote the understanding of miRNA functions during the cold response. The miRNA genotype-specific expression model might explain the distinct cold sensitivities between tea lines.

Small RNA and degradome profiling reveals a role for miRNAs and their targets in the developing fibers of Gossypium barbadense

microRNAs (miRNAs) are 20-24 nucleotide non-coding small RNAs that play important roles in plant development. The stages of cotton fiber development include initiation, elongation, secondary wall thickening (SWT) and maturation. We constructed seven fiber RNA libraries representing the initiation, elongation and SWT stages. In total, 47 conserved miRNA families and seven candidate miRNAs were profiled using small RNA sequencing. Northern blotting and real-time polymerase chain reaction (PCR) analyses revealed the dynamic expression of miRNAs during fiber development. In addition, 140 targets of 30 conserved miRNAs and 38 targets of five candidate miRNAs were identified through degradome sequencing. Analysis of correlated expression between miRNAs and their targets demonstrated that specific miRNAs suppressed the expression of transcription factors, SBP and MYB, a leucine-rich receptor-like protein kinase, a pectate lyase, α-tubulin, a UDP-glucuronic acid decarboxylase and cytochrome C oxidase subunit 1 to affect fiber development. Histochemical analyses detected the biological activity of miRNA156/157 in ovule and fiber development. Suppressing miRNA156/157 function resulted in the reduction of mature fiber length, illustrating that miRNA156/157 plays an essential role in fiber elongation.

Advances in identification and validation of plant microRNAs and their target genes

Developments in the field of molecular biology and genetics, such as microarray, gene transfer and discovery of small regulatory RNAs, have led to significant advances in plant biotechnology. Among the small RNAs, microRNAs (miRNAs) have elicited much interest as key post-transcriptional regulators in eukaryotic gene expression. Advances in genome and transcriptome sequencing of plants have facilitated the generation of a huge wealth of sequence information that can find much use in the discovery of novel miRNAs and their target genes. In this review, we present an overview of the developments in the strategies and methods used to identify and study miRNAs, their target genes and the mechanisms by which these miRNAs interact with their target genes since the discovery of the first miRNA. The approaches discussed include both reverse and forward genetics. We observed that despite the availability of advanced methods, certain limitations ranging from the cost of materials, equipment and personnel to the availability of genome sequences for many plant species present a number of challenges for the development and utilization of modern scientific methods for the elucidation and development of miRNAs in many important plant species.

MiR397b regulates both lignin content and seed number in Arabidopsis via modulating a laccase involved in lignin biosynthesis

Plant laccase (LAC) enzymes belong to the blue copper oxidase family and polymerize monolignols into lignin. Recent studies have established the involvement of microRNAs in this process; however, physiological functions and regulation of plant laccases remain poorly understood. Here, we show that a laccase gene, LAC4, regulated by a microRNA, miR397b, controls both lignin biosynthesis and seed yield in Arabidopsis. In transgenic plants, overexpression of miR397b (OXmiR397b) reduced lignin deposition. The secondary wall thickness of vessels and the fibres was reduced in the OXmiR397b line, and both syringyl and guaiacyl subunits are decreased, leading to weakening of vascular tissues. In contrast, overexpression of miR397b-resistant laccase mRNA results in an opposite phenotype. Plants overexpressing miR397b develop more than two inflorescence shoots and have an increased silique number and silique length, resulting in higher seed numbers. In addition, enlarged seeds and more seeds are formed in these miR397b overexpression plants. The study suggests that miR397-mediated development via regulating laccase genes might be a common mechanism in flowering plants and that the modulation of laccase by miR397 may be potential for engineering plant biomass production with less lignin.

Nitrogen signalling pathways shaping root system architecture: an update

Root system architecture is a fundamentally important trait for resource acquisition in both ecological and agronomic contexts. Because of the plasticity of root development and the almost infinite complexity of the soil, root system architecture is shaped by environmental factors to a much greater degree than shoot architecture. In attempting to understand how roots sense and respond to environmental cues, the striking effects of nitrate and other forms of nitrogen on root growth and branching have received particular attention. This minireview focuses on the latest advances in our understanding of the diverse nitrogen signalling pathways that are now known to act at multiple stages in the process of lateral root development, as well as on primary root growth.

sPARTA: a parallelized pipeline for integrated analysis of plant miRNA and cleaved mRNA data sets, including new miRNA target-identification software

Parallel analysis of RNA ends (PARE) is a technique utilizing high-throughput sequencing to profile uncapped, mRNA cleavage or decay products on a genome-wide basis. Tools currently available to validate miRNA targets using PARE data employ only annotated genes, whereas important targets may be found in unannotated genomic regions. To handle such cases and to scale to the growing availability of PARE data and genomes, we developed a new tool, 'sPARTA' (small RNA-PARE target analyzer) that utilizes a built-in, plant-focused target prediction module (aka 'miRferno'). sPARTA not only exhibits an unprecedented gain in speed but also it shows greater predictive power by validating more targets, compared to a popular alternative. In addition, the novel 'seed-free' mode, optimized to find targets irrespective of complementarity in the seed-region, identifies novel intergenic targets. To fully capitalize on the novelty and strengths of sPARTA, we developed a web resource, 'comPARE', for plant miRNA target analysis; this facilitates the systematic identification and analysis of miRNA-target interactions across multiple species, integrated with visualization tools. This collation of high-throughput small RNA and PARE datasets from different genomes further facilitates re-evaluation of existing miRNA annotations, resulting in a 'cleaner' set of microRNAs.

Differentially expressed microRNA cohorts in seed development may contribute to poor grain filling of inferior spikelets in rice

BACKGROUND

The inferior spikelets are defined to be those at portions where the grains receive less photosynthetic products during the seed development. The typical inferior spikelets are physically located on the proximal secondary branches in a rice panicle and traditionally characterized by a later flowering time and a slower grain-filling rate, compared to those so-called superior spikelets. Grains produced on the inferior spikelets are consequently under-developed and lighter in weight than those formed on the superior spikelets. MicroRNAs (miRNAs) are recognized as key players in regulating plant development through post-transcriptional gene regulations. We previously presented the evidence that miRNAs may influence grain-filling rate and played a role in determining the grain weight and yield in rice.

RESULTS

In this study, further analyses of the expressed small RNAs in superior and inferior spikelets were conducted at five distinct developmental stages of grain development. Totally, 457 known miRNAs and 13 novel miRNAs were analyzed, showing a differential expression of 141 known miRNAs between superior and inferior spikelets with higher expression levels of most miRNAs associated with the superior than the inferior spikelets during the early stage of grain filling. Genes targeted by those differentially expressed miRNAs (i.e. miR156, miR164, miR167, miR397, miR1861, and miR1867) were recognized to play roles in multiple developmental and signaling pathways related to plant hormone homeostasis and starch accumulation.

CONCLUSIONS

Our data established a complicated link between miRNA dynamics and the traditional role of hormones in grain filling and development, providing new insights into the widely accepted concepts of the so-called superior and inferior spikelets in rice production.

Bioinformatic identification and experimental validation of miRNAs from foxtail millet (Setaria italica)

MiRNAs are a novel group of non-coding small RNAs that negatively regulate gene expression. Many miRNAs have been identified and investigated extensively in plant species with sequenced genomes. However, few miRNAs have been identified in foxtail millet (Setaria italica), which is an ancient cereal crop of great importance for dry land agriculture. In this study, 271 foxtail millet miRNAs belonging to 44 families were identified using a bioinformatics approach. Twenty-three pairs of sense/antisense miRNAs belonging to 13 families, and 18 miRNA clusters containing members of 8 families were discovered in foxtail millet. We identified 432 potential targets for 38 miRNA families, most of which were predicted to be involved in plant development, signal transduction, metabolic pathways, disease resistance, and environmental stress responses. Gene ontology (GO) analysis revealed that 101, 56, and 23 target genes were involved in molecular functions, biological processes, and cellular components, respectively. We investigated the expression patterns of 43 selected miRNAs using qRT-PCR analysis. All of the miRNAs were expressed ubiquitously with many exhibiting different expression levels in different tissues. We validated five predicted targets of four miRNAs using the RNA ligase mediated rapid amplification of cDNA end (5'-RLM-RACE) method.

Whole-genome discovery of miRNAs and their targets in wheat (Triticum aestivum L.)

BACKGROUND

MicroRNAs (miRNAs) are small, non-coding RNAs playing essential roles in plant growth, development, and stress responses. Sequencing of small RNAs is a starting point for understanding their number, diversity, expression and possible roles in plants.

RESULTS

In this study, we conducted a genome-wide survey of wheat miRNAs from 11 tissues, characterizing a total of 323 novel miRNAs belonging to 276 families in wheat. A miRNA conservation analysis identified 191 wheat-specific miRNAs, 2 monocot-specific miRNAs, and 30 wheat-specific variants from 9 highly conserved miRNA families. To understand possible roles of wheat miRNAs, we determined 524 potential targets for 124 miRNA families through degradome sequencing, and cleavage of a subset of them was validated via 5' RACE. Based on the genome-wide identification and characterization of miRNAs and their associated target genes, we further identified 64 miRNAs preferentially expressing in developing or germinating grains, which could play important roles in grain development.

CONCLUSION

We discovered 323 wheat novel miRNAs and 524 target genes for 124 miRNA families in a genome-wide level, and our data will serve as a foundation for future research into the functional roles of miRNAs in wheat.

A comparison of performance of plant miRNA target prediction tools and the characterization of features for genome-wide target prediction

BACKGROUND

Deep-sequencing has enabled the identification of large numbers of miRNAs and siRNAs, making the high-throughput target identification a main limiting factor in defining their function. In plants, several tools have been developed to predict targets, majority of them being trained on Arabidopsis datasets. An extensive and systematic evaluation has not been made for their suitability for predicting targets in species other than Arabidopsis. Nor, these have not been evaluated for their suitability for high-throughput target prediction at genome level.

RESULTS

We evaluated the performance of 11 computational tools in identifying genome-wide targets in Arabidopsis and other plants with procedures that optimized score-cutoffs for estimating targets. Targetfinder was most efficient [89% 'precision' (accuracy of prediction), 97% 'recall' (sensitivity)] in predicting 'true-positive' targets in Arabidopsis miRNA-mRNA interactions. In contrast, only 46% of true positive interactions from non-Arabidopsis species were detected, indicating low 'recall' values. Score optimizations increased the 'recall' to only 70% (corresponding 'precision': 65%) for datasets of true miRNA-mRNA interactions in species other than Arabidopsis. Combining the results of Targetfinder and psRNATarget delivers high true positive coverage, whereas the intersection of psRNATarget and Tapirhybrid outputs deliver highly 'precise' predictions. The large number of 'false negative' predictions delivered from non-Arabidopsis datasets by all the available tools indicate the diversity in miRNAs-mRNA interaction features between Arabidopsis and other species. A subset of miRNA-mRNA interactions differed significantly for features in seed regions as well as the total number of matches/mismatches.

CONCLUSION

Although, many plant miRNA target prediction tools may be optimized to predict targets with high specificity in Arabidopsis, such optimized thresholds may not be suitable for many targets in non-Arabidopsis species. More importantly, non-conventional features of miRNA-mRNA interaction may exist in plants indicating alternate mode of miRNA target recognition. Incorporation of these divergent features would enable next-generation of algorithms to better identify target interactions.

Conserved miR164-targeted NAC genes negatively regulate drought resistance in rice

MicroRNAs constitute a large group of endogenous small RNAs of ~22 nt that emerge as vital regulators, mainly by targeting mRNAs for post-transcriptional repression. Previous studies have revealed that the miR164 family in Arabidopsis is comprised of three members which guide the cleavage of the mRNAs of five NAC genes to modulate developmental processes. However, the functions of the miR164-targeted NAC genes in crops are poorly deciphered. In this study, the conserved features of six miR164-targeted NAC genes (OMTN1-OMTN6) in rice are described, and evidence is provided that four of them confer a negative regulatory role in drought resistance. OMTN proteins have the characteristics of typical NAC transcriptional factors. The miR164 recognition sites of the OMTN genes are highly conserved in rice germplasms. Deletion of the recognition sites impaired the transactivation activity, indicating that the conserved recognition sites play a crucial role in maintaining the function of the OMTN proteins. The OMTN genes were responsive to abiotic stresses, and showed diverse spatio-temporal expression patterns in rice. Overexpression of OMTN2, OMTN3, OMTN4, and OMTN6 in rice led to negative effects on drought resistance at the reproductive stage. The expression of numerous genes related to stress response, development, and metabolism was altered in OMTN2-, OMTN3-, OMTN4-, and OMTN6-overexpressing plants. Most of the up-regulated genes in the OMTN-overexpressing plants were down-regulated by drought stress. The results suggest that the conserved miR164-targeted NAC genes may be negative regulators of drought tolerance in rice, in addition to their reported roles in development.

Long non-coding RNAs: a novel endogenous source for the generation of Dicer-like 1-dependent small RNAs in Arabidopsis thaliana

The biological relevance of long non-coding RNAs (lncRNAs) is emerging. Whether the lncRNAs could form structured precursors for small RNAs (sRNAs) production remains elusive. Here, 172 713 DCL1 (Dicer-like 1)-dependent sRNAs were identified in Arabidopsis. Except for the sRNAs mapped onto the microRNA precursors, the remaining ones led us to investigate their originations. Intriguingly, 65 006 sRNAs found their loci on 5891 lncRNAs. These sRNAs were sent to AGO (Argonaute) enrichment analysis. As a result, 1264 sRNAs were enriched in AGO1, which were then subjected to target prediction. Based on degradome sequencing data, 109 transcripts were validated to be targeted by 96 sRNAs. Besides, 44 lncRNAs were targeted by 23 sRNAs. To further support the origination of the DCL1-dependent sRNAs from lncRNAs, we searched for the degradome-based cleavage signals at either ends of the sRNA loci, which were supposed to be produced during DCL1-mediated processing of the long-stem structures. As a result, 63 612 loci were supported by degradome signatures. Among these loci, 6606 reside within the dsRNA-seq (double-stranded RNA sequencing) read-covered regions of 100 nt or longer. These regions were subjected to secondary structure prediction. And, 43 regions were identified to be capable of forming highly complementary long-stem structures. We proposed that these local long-stem structures could be recognized by DCL1 for cropping, thus serving as the sRNA precursors. We hope that our study could inspire more research efforts to study on the biological roles of the lncRNAs in plants.

miR444a has multiple functions in the rice nitrate-signaling pathway

Nitrate (NO3-) is a key signaling molecule in plant metabolism and development, in addition to its role as a nutrient. It has been shown previously in Arabidopsis that ANR1, a MADS-box transcription factor, is a major component in the NO3--signaling pathway that triggers lateral root growth and that miR444, which is specific to monocots, targets four genes that are homologous to ANR1 in rice. Here, we show that miR444a plays multiple roles in the rice NO3--signaling pathway - not only in root development, but also involving nitrate accumulation and even Pi -starvation responses. miR444a overexpression resulted in reduced rice lateral root elongation, but promoted rice primary and adventitious root growth, in a nitrate-dependent manner. In addition, overexpression of miR444a improved nitrate accumulation and expression of nitrate transporter genes under high nitrate concentration conditions, but reduced the remobilization of nitrate from old leaves to young leaves thus affecting the plant's ability to adapt to nitrogen-limitating conditions. Intriguingly, we found that Pi starvation strongly induced miR444 accumulation in rice roots and that overexpression of miR444a altered Pi -starvation-induced root architecture and enhanced Pi accumulation and expression of three Pi transporter genes. We further provide evidence that miR444a is involved in the interaction between the NO3--signaling and Pi -signaling pathways in rice. Taken together, our observations demonstrated that miR444a plays multiple roles in the rice NO3--signaling pathway in nitrate-dependent root growth, nitrate accumulation and phosphate-starvation responses.

Mla- and Rom1-mediated control of microRNA398 and chloroplast copper/zinc superoxide dismutase regulates cell death in response to the barley powdery mildew fungus

• Barley (Hordeum vulgare L.) Mildew resistance locus a (Mla) confers allele-specific interactions with natural variants of the ascomycete fungus Blumeria graminis f. sp. hordei (Bgh), the causal agent of powdery mildew disease. Significant reprogramming of Mla-mediated gene expression occurs upon infection by this obligate biotrophic pathogen. • We utilized a proteomics-based approach, combined with barley mla, required for Mla12 resistance1 (rar1), and restoration of Mla resistance1 (rom1) mutants, to identify components of Mla-directed signaling. • Loss-of-function mutations in Mla and Rar1 both resulted in the reduced accumulation of chloroplast copper/zinc superoxide dismutase 1 (HvSOD1), whereas loss of function in Rom1 re-established HvSOD1 levels. In addition, both Mla and Rom1 negatively regulated hvu-microRNA398 (hvu-miR398), and up-regulation of miR398 was coupled to reduced HvSOD1 expression. Barley stripe mosaic virus (BSMV)-mediated over-expression of both barley and Arabidopsis miR398 repressed accumulation of HvSOD1, and BSMV-induced gene silencing of HvSod1 impeded Mla-triggered H₂O₂ and hypersensitive reaction (HR) at barley-Bgh interaction sites. • These data indicate that Mla- and Rom1-regulated hvu-miR398 represses HvSOD1 accumulation, influencing effector-induced HR in response to the powdery mildew fungus.

Characterization of regulatory mechanism of Poncirus trifoliata microRNAs on their target genes with an integrated strategy of newly developed PPM-RACE and RLM-RACE

MicroRNAs (miRNAs) play an important role in post-transcriptional gene regulation that involved various biological and metabolic processes. Many extensive studies have been done in model plant species, to discover miRNAs' regulating expression of their target genes and analyze their functions. But, the function of Poncirus trifoliata miRNAs has not been properly investigated. In this study, we employed the RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-RACE) and the newly developed method called poly (A) polymerase-mediated 3' rapid amplification of cDNA ends (PPM-RACE), which mapped the cleavage site of target mRNAs and detected expression patterns of cleaved fragments that could in turn indicate the regulatory functions of the miRNAs on their target genes. Furthermore, the spatiotemporal expression levels of target genes were analyzed by qRT-PCR, with exhibiting different expression trends from their corresponding miRNAs, thus indicating the cleavage mode of miRNAs on their target genes. The expression patterns of miRNAs, their target mRNAs and cleaved target mRNAs in different organs of juvenile and adult trifoliate orange were studied. The results showed that the expression of miRNAs and their target mRNAs was in a trade-off trend. When the miRNA expression was high, its corresponding target mRNA expression was low, while the cleaved target mRNA expression was high; when the miRNA expression was low, its target mRNA expression was high, while the expression of cleaved target mRNAs follows that of the miRNA. The validation of the cleavage site of target mRNAs and the detection of expression patterns of cleaved fragments can further broaden the knowledge of small RNA-mediated regulation in P. trifoliate.

GC content fluctuation around plant small RNA-generating sites

GC content of small RNA-generating sites and their flanking sequences in Arabidopsis thaliana and rice was systematically analyzed in silico. High GC content fluctuation (GCF) is observed at the borders of sRNA sites, while the GCF within sRNA sites is low. Furthermore, the GC content along sequences of some miniature inverted-repeat transposable element (MITE) families coincides with the abundance of MITE-derived small RNAs. The GCF within tasiRNA clusters is negatively correlated with its phasing score. We conclude that high GC content and low GCF could increase the expression of small RNA. Our results provide further insights on small RNA expression, which may be applied to improve the silencing efficiency of RNAi and virus-induced gene silencing.

Beyond cleaved small RNA targets: unraveling the complexity of plant RNA degradome data

Background

Degradation is essential for RNA maturation, turnover, and quality control. RNA degradome sequencing that integrates a modified 5′-rapid amplification of cDNA ends protocol with next-generation sequencing technologies is a high-throughput approach for profiling the 5′-end of uncapped RNA fragments on a genome-wide scale. The primary application of degradome sequencing has been to identify the truncated transcripts that result from endonucleolytic cleavage guided by microRNAs or small interfering RNAs. As many pathways are involved in RNA degradation, degradome data should contain other RNA species besides the cleavage remnants of small RNA targets. Nevertheless, no systematic approaches have been established to explore the hidden complexity of plant degradome.

Results

Through analyzing Arabidopsis and rice RNA degradome data, we recovered 11 short motifs adjacent to predominant and abundant uncapped 5′-ends. Uncapped ends associated with several of these short motifs were more prevalent than those targeted by most miRNA families especially in the 3′ untranslated region of transcripts. Through genome-wide analysis, five motifs showed preferential accumulation of uncapped 5′-ends at the same position in Arabidopsis and rice. Moreover, the association of uncapped 5′-ends with a CA-repeat motif and a motif recognized by Pumilio/Fem-3 mRNA binding factor (PUF) proteins was also found in non-plant species, suggesting that common mechanisms are present across species. Based on these motifs, potential sources of RNA ends that constitute degradome data were proposed and further examined. The 5′-end of small nucleolar RNAs could be precisely captured by degradome sequencing. Position-specific enrichment of uncapped 5′-ends was seen upstream of motifs recognized by several RNA binding proteins especially for the binding site of PUF proteins. False uncapped 5′-ends produced from capped transcripts through non-specific PCR amplification were common artifacts among degradome datasets.

Conclusions

The complexity of plant RNA degradome data revealed in this study may contribute to the alternative applications of degradome in RNA research.

A computational model for non-conserved mature miRNAs from the rice genome

Several computational approaches employ the high complementarity of plant miRNAs to target mRNAs as a filter to recognize miRNA. Numerous non-conserved miRNAs are known with more recent evolutionary origin as a result of target gene duplication events. We present here a computational model with knowledge inputs from reported non-conserved mature miRNAs of Oryza sativa (rice). Sequence- and structure-based approaches were used to retrieve miRNA features based on rice Argonaute protein and develop a multiple linear regression (MLR) model (r(2) = 0.996, q(2)cv = 0.989) which scored mature miRNAs as predicted by the MaturePred program. The model was validated by scoring test set (q(2) = 0.990) and computationally predicted mature miRNAs as external test set (q(2)test = 0.895). This strategy successfully enhanced the confidence of retrieving most probable non-conserved miRNAs from the rice genome. We anticipate that this computational model would recognize unknown non-conserved miRNA candidates and nurture the current mechanistic understanding of miRNA sorting to unveil the role of non-conserved miRNAs in gene silencing.

Redox signaling mediates the expression of a sulfate-deprivation-inducible microRNA395 in Arabidopsis

MicroRNA395 (miR395) is a conserved miRNA that targets a low-affinity sulfate transporter (AST68) and three ATP sulfurylases (APS1, APS3 and APS4) in higher plants. In this study, At2g28780 was confirmed as another target of miR395 in Arabidopsis. Interestingly, several dicots contained genes homologous to At2g28780 and a cognate miR395 complementary site but possess a gradient of mismatches at the target site. It is well established that miR395 is induced during S deprivation in Arabidopsis; however, the signaling pathways that mediate this regulation are unknown. Several findings in the present study demonstrate that redox signaling plays an important role in induction of miR395 during S deprivation. These include the following results: (i) glutathione (GSH) supplementation suppressed miR395 induction in S-deprived plants (ii) miR395 is induced in Arabidopsis seedlings exposed to Arsenate or Cu(2+) , which induces oxidative stress (iii), S deprivation-induced oxidative stress, and (iv) compromised induction of miR395 during S deprivation in cad2 mutant (deficient in GSH biosynthesis) that is defective in glutaredoxin-dependent redox signaling and ntra/ntrb (defective in thioredoxin reductases a and b) double mutants that are defective in thioredoxin-dependent redox signaling. Collectively, these findings strongly support the involvement of redox signaling in inducing the expression of miR395 during S deprivation in Arabidopsis.

Novel insights from non-conserved microRNAs in plants

Plant microRNAs (miRNAs), a class of small non-coding regulatory RNAs, are canonically 20-24 nucleotides in length and bind to complementary target RNA sequences, guiding target attenuation via mRNA degradation or translation inhibition. Of the annotated miRNA families, evolutionarily conserved families have been well known to extensively regulate analogous targets and play critical roles in plant development and adaptation to adverse environments. By contrast, majority of these families that are merely present in a specific lineage or in a few closely related species have not been well functionally explored until recently. The fast-growing progresses being made in the actions of non-conserved miRNAs nowadays in diverse plant species may represent a highly promising research field in future. This review thereby summarizes the emerging advances in our understanding of the biogenesis, associated effectors, modes to targets, and biological functions of plant non-conserved miRNAs. In addition, it outlines the regulatory units recently discovered between conserved miRNAs and their alternative targets.

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon

BACKGROUND

The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon.

RESULTS

B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner.

CONCLUSIONS

B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.

Identification and comparative analysis of cadmium tolerance-associated miRNAs and their targets in two soybean genotypes

MicroRNAs (miRNAs) play crucial roles in regulating the expression of various stress responses genes in plants. To investigate soybean (Glycine max) miRNAs involved in the response to cadmium (Cd), microarrays containing 953 unique miRNA probes were employed to identify differences in the expression patterns of the miRNAs between different genotypes, Huaxia3 (HX3, Cd-tolerant) and Zhonghuang24 (ZH24, Cd-sensitive). Twenty six Cd-responsive miRNAs were identified in total. Among them, nine were detected in both cultivars, while five were expressed only in HX3 and 12 were only in ZH24. The expression of 16 miRNAs was tested by qRT-PCR and most of the identified miRNAs were found to have similar expression patterns with microarray. Three hundred and seventy six target genes were identified for 204 miRNAs from a mixture degradome library, which was constructed from the root of HX3 and ZH24 with or without Cd treatment. Fifty five genes were identified to be cleaved by 14 Cd-responsive miRNAs. Gene ontology (GO) annotations showed that these target transcripts are implicated in a broad range of biological processes. In addition, the expression patterns of ten target genes were validated by qRT-PCR. The characterization of the miRNAs and the associated target genes in response to Cd exposure provides a framework for understanding the molecular mechanism of heavy metal tolerance in plants.

miRTarBase update 2014: an information resource for experimentally validated miRNA-target interactions

MicroRNAs (miRNAs) are small non-coding RNA molecules capable of negatively regulating gene expression to control many cellular mechanisms. The miRTarBase database (http://mirtarbase.mbc.nctu.edu.tw/) provides the most current and comprehensive information of experimentally validated miRNA-target interactions. The database was launched in 2010 with data sources for >100 published studies in the identification of miRNA targets, molecular networks of miRNA targets and systems biology, and the current release (2013, version 4) includes significant expansions and enhancements over the initial release (2010, version 1). This article reports the current status of and recent improvements to the database, including (i) a 14-fold increase to miRNA-target interaction entries, (ii) a miRNA-target network, (iii) expression profile of miRNA and its target gene, (iv) miRNA target-associated diseases and (v) additional utilities including an upgrade reminder and an error reporting/user feedback system.

SRNAome parsing yields insights into tomato fruit ripening control

Small RNAs have emerged as critical regulators in the expression and function of eukaryotic genomes at the post-transcriptional level. To elucidate the functions of microRNA (miRNAs) and endogenous small-interfering RNAs (siRNAs) in tomato fruit ripening process, the deep sequencing and bioinformatics methods were combined to parse the small RNAs landscape in three fruit-ripening stages (mature green, breaker and red-ripe) on a whole genome. Two species-specific miRNAs and two members of TAS3 family were identified, 590 putative phased small RNAs and 125 cis-natural antisense (nat-siRNAs) were also found in our results which enriched the tomato small RNAs repository and all of them showed differential expression patterns during fruit ripening. A large amount of the targets of the small RNAs were predicted to be involved in fruit ripening and ethylene pathway. Furthermore, the promoters of the conserved and novel miRNAs were found to contain the conserved motifs of TATA-box and CT microsatellites which were also found in Arabidopsis and rice, and several species-specific motifs were found in parallel.

Genome-wide identification and characterization of cadmium-responsive microRNAs and their target genes in radish (Raphanus sativus L.) roots

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that play vital regulatory roles in plant growth, development, and environmental stress responses. Cadmium (Cd) is a non-essential heavy metal that is highly toxic to living organisms. To date, a number of conserved and non-conserved miRNAs have been identified to be involved in response to Cd stress in some plant species. However, the miRNA-mediated gene regulatory networks responsive to Cd stress in radish (Raphanus sativus L.) remain largely unexplored. To dissect Cd-responsive miRNAs and their targets systematically at the global level, two small RNA libraries were constructed from Cd-treated and Cd-free roots of radish seedlings. Using Solexa sequencing technology, 93 conserved and 16 non-conserved miRNAs (representing 26 miRNA families) and 28 novel miRNAs (representing 22 miRNA families) were identified. In all, 15 known and eight novel miRNA families were significantly differently regulated under Cd stress. The expression patterns of a set of Cd-responsive miRNAs were validated by quantitative real-time PCR. Based on the radish mRNA transcriptome, 18 and 71 targets for novel and known miRNA families, respectively, were identified by the degradome sequencing approach. Furthermore, a few target transcripts including phytochelatin synthase 1 (PCS1), iron transporter protein, and ABC transporter protein were involved in plant response to Cd stress. This study represents the first transcriptome-based analysis of miRNAs and their targets responsive to Cd stress in radish roots. These findings could provide valuable information for functional characterization of miRNAs and their targets in regulatory networks responsive to Cd stress in radish.

Microarray and degradome sequencing reveal microRNA differential expression profiles and their targets in Pinellia pedatisecta

MicroRNAs (miRNAs) are endogenous small non-coding RNAs which play a critical role in gene regulation in plants. Pinelliapedatisecta is one of the most important herbs in traditional Chinese medicine, but there are no microRNAs of Pinelliapedatisecta were deposited in miRBase and the research of the related miRNA biological functions is still insufficient. To detect Pinelliapedatisecta miRNAs and discover their expression difference with Pinelliaternata, we carried out a microarray profiling. A total of 101 miRNAs belonging to 22 miRNA families were detected both in Pinelliapedatisecta and Pinelliaternata respectively, among them 21 miRNAs showed their differentially expression. GO (gene ontology) term enrichment analysis of the target genes of differential expression miRNAs reveal that these miRNAs mainly affect the reproduction, transcription factor activity and plant developmental process. To elucidate the target function of miRNAs, we constructed a degradome library from Pinellia pedatisecta leaf. The result showed that a total of 18 transcript were identified as targets of miRNAs and further analysis indicated that miR156 and miR529 may function together to repress SPL14.

A global profiling of uncapped mRNAs under cold stress reveals specific decay patterns and endonucleolytic cleavages in Brachypodium distachyon

Background

mRNA degradation is a critical factor in determining mRNA abundance and enables rapid adjustment of gene expression in response to environmental stress. The involvement of processing bodies in stress response suggests a role for decapping-mediated mRNA degradation. However, little is known about the role of mRNA degradation under stressful environmental conditions.

Results

Here, we perform a global study of uncapped mRNAs, via parallel analysis of RNA ends (PARE), under cold stress in Brachypodium distachyon. Enrichment analysis indicates that degradation products detected by PARE are mainly generated by the decapping pathway. Endonucleolytic cleavages are detected, uncovering another way of modulating gene expression. PARE and RNA-Seq analyses identify four types of mRNA decay patterns. Type II genes, for which light-harvesting processes are over-represented in gene ontology analyses, show unchanged transcript abundance and altered uncapped transcript abundance. Uncapping-mediated transcript stability of light harvesting-related genes changes significantly in response to cold stress, which may allow rapid adjustments in photosynthetic activity in response to cold stress. Transcript abundance and uncapped transcript abundance for type III genes changes in opposite directions in response to cold stress, indicating that uncapping-mediated mRNA degradation plays a role in regulating gene expression.

Conclusion

To our knowledge, this is the first global analysis of mRNA degradation under environmental stress conditions in Brachypodium distachyon. We uncover specific degradation and endonucleolytic cleavage patterns under cold stress, which will deepen our understanding of mRNA degradation under stressful environmental conditions, as well as the cold stress response mechanism in monocots.

Construction of small RNA-mediated gene regulatory networks in the roots of rice (Oryza sativa)

Background

The root systems play essential roles for plants to anchorage to the soil, and to exploit the mineral and water resources. The molecular mechanisms underlying root development have been extensively studied to improve root system architecture, especially for the crops. Several microRNA (miRNA) families have been demonstrated to be involved in plant root development. However, whether the other small RNA (sRNA) species, which occupy a dominant portion of the plant endogenous sRNA population, possess potential roles in root development remains unclear.

Results

In this study, by using sRNA high-throughput sequencing data, we made a comparison of the sRNA accumulation levels between the rice root tips and the whole roots. The sRNAs highly accumulated in the root tips and in the whole roots were extracted respectively. After Argonaute 1 (AGO1) enrichment analysis, the sRNAs with great potential of performing target cleavages were included for target prediction and degradome sequencing data-based validation. As a result, lists of the targets regulated by the AGO1-enriched sRNAs were obtained for both the root tips and the whole roots. Further evidences were identified from microarray data of the target genes to support some of the sRNA—target interactions. Specifically, the expression patterns of certain target genes in the root tips and the whole roots were contrary to those of the regulating sRNAs. Besides, several targets were indicated to play important roles in root development based on literature mining.

Conclusions

Taken together, the regulatory networks mediated by the sRNAs highly accumulated in the root tips or in the whole roots could advance our current understanding of the sRNA-involved molecular mechanisms underlying rice root development. And, the sRNA—target lists could serve as the basis for further functional investigations.

Identification of a novel microRNA (miRNA) from rice that targets an alternatively spliced transcript of the Nramp6 (Natural resistance-associated macrophage protein 6) gene involved in pathogen resistance

Plants have evolved efficient defence mechanisms to defend themselves from pathogen attack. Although many studies have focused on the transcriptional regulation of defence responses, less is known about the involvement of microRNAs (miRNAs) as post-transcriptional regulators of gene expression in plant immunity. This work investigates miRNAs that are regulated by elicitors from the blast fungus Magnaporthe oryzae in rice (Oryza sativa). Small RNA libraries were constructed from rice tissues and subjected to high-throughput sequencing for the identification of elicitor-responsive miRNAs. Target gene expression was examined by microarray analysis. Transgenic lines were used for the analysis of miRNA functioning in disease resistance. Elicitor treatment is accompanied by dynamic alterations in the expression of a significant number of miRNAs, including new members of annotated miRNAs. Novel miRNAs from rice are proposed. We report a new rice miRNA, osa-miR7695, which negatively regulates an alternatively spliced transcript of OsNramp6 (Natural resistance-associated macrophage protein 6). This novel miRNA experienced natural and domestication selection events during evolution, and its overexpression in rice confers pathogen resistance. This study highlights an miRNA-mediated regulation of OsNramp6 in disease resistance, whilst illustrating the existence of a novel regulatory network that integrates miRNA function and mRNA processing in plant immunity.

Evolutionary conservation of microRNA regulatory programs in plant flower development

MicroRNAs (miRNAs) are post-transcriptional regulators of growth and development in both plants and animals. Flowering is critical for the reproduction of angiosperms. Flower development entails the transition from vegetative growth to reproductive growth, floral organ initiation, and the development of floral organs. These developmental processes are genetically regulated by miRNAs, which participate in complex genetic networks of flower development. A survey of the literature shows that miRNAs, their specific targets, and the regulatory programs in which they participate are conserved throughout the plant kingdom. This review summarizes the role of miRNAs and their targets in the regulation of gene expression during the floral developmental phase, which includes the floral transition stage, followed by floral patterning, and then the development of floral organs. The conservation patterns observed in each component of the miRNA regulatory system suggest that these miRNAs play important roles in the evolution of flower development.

Introns targeted by plant microRNAs: a possible novel mechanism of gene regulation

BACKGROUND

In plant cells, most microRNAs (miRNAs) perform cleavages of target mature mRNAs in the cytoplasm. A recent report of a miRNA pathway involved in DNA methylation in the rice nucleus raises the possibility that plant miRNAs could cleave intron-containing pre-mRNAs (the precursor of messenger RNAs) located in the nucleus.

RESULTS

In this study, we searched for the miRNA binding sites present within the introns of Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) genes. All miRNA-intron interactions predicted to result in cleavages were validated by using the public degradome sequencing data. As a result, 40 miRNA-intron pairs involving 25 miRNAs in Arabidopsis and 1912 pairs involving 91 miRNAs in rice were identified. For several rice genes, not all transcription forms (alternative splicing variants) were under similar regulation by specific miRNAs. Certain transcripts could escape cleavages due to the absence of intronic miRNA binding sites within these sequences. In some instances, specific cleaved intron remnants could be converted to double-stranded RNAs (dsRNAs) by RNA-dependent RNA polymerase 2. These dsRNAs could then be processed into 21- and 24-nt phased sRNAs by the activity of Dicer-like 1 and 3, respectively. The resultant siRNAs have the potential to be incorporated into Argonaute (AGO)-associated silencing complexes and result in cleavages of target pre-mRNA sequences.

CONCLUSIONS

A regulatory model, miRNA-targeting of intron-containing pre-mRNAs-phased sRNAs-targeting of mature mRNAs is proposed, which further expands the potential modes of action of plant miRNAs.

Comparative analysis reveals dynamic changes in miRNAs and their targets and expression during somatic embryogenesis in longan (Dimocarpus longan Lour.)

Somatic embryogenesis (SE), which resembles zygotic embryogenesis, is an essential component of the process of plant cell differentiation and embryo development. Although microRNAs (miRNAs) are important regulators of many plant develop- mental processes, their roles in SE have not been thoroughly investigated. In this study, we used deep-sequencing, computational, and qPCR methods to identify, profile, and describe conserved and novel miRNAs involved in longan (Dimocarpus longan) SE. A total of 643 conserved and 29 novel miRNAs (including star strands) from more than 169 miRNA families were identified in longan embryogenic tissue using Solexa sequencing. By combining computational and degradome sequencing approaches, we were able to predict 2063 targets of 272 miRNAs and verify 862 targets of 181 miRNAs. Target annotation revealed that the putative targets were involved in a broad variety of biological processes, including plant metabolism, signal transduction, and stimulus response. Analysis of stage- and tissue-specific expressions of 20 conserved and 4 novel miRNAs indicated their possible roles in longan SE. These miRNAs were dlo-miR156 family members and dlo-miR166c* associated with early embryonic culture developmental stages; dlo-miR26, dlo-miR160a, and families dlo-miR159, dlo-miR390, and dlo-miR398b related to heart-shaped and torpedo- shaped embryo formation; dlo-miR4a, dlo-miR24, dlo-miR167a, dlo-miR168a*, dlo-miR397a, dlo-miR398b.1, dlo-miR398b.2, dlo-miR808 and dlo-miR5077 involved in cotyledonary embryonic development; and dlo-miR17 and dlo-miR2089*-1 that have regulatory roles during longan SE. In addition, dlo-miR167a, dlo-miR808, and dlo-miR5077 may be required for mature embryo formation. This study is the first reported investigation of longan SE involving large-scale cloning, characterization, and expression profiling of miRNAs and their targets. The reported results contribute to our knowledge of somatic embryo miRNAs and provide insights into miRNA biogenesis and expression in plant somatic embryo development.