Transcriptome-wide identification of microRNA targets in rice

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.

Characterization of small RNAs and their target genes in wheat seedlings using sequencing-based approaches

Wheat is the most highly cultivated plant species for its grain production throughout the world. Because small RNA-dependent gene regulation is critical for successful completion of plant life cycle including its productivity, identification of not only miRNAs but also confirming their targets in wheat is important. To identify small RNAs including novel miRNAs as well as miRNA targets in wheat, we constructed small RNA and degradome libraries from wheat seedlings. Small RNA analysis resulted in identification of most conserved miRNAs including novel miRNAs that can be grouped into 32 miRNA families. The sequence analysis also led to the characterization of two abundantly expressed rRNA-derived small RNAs. To identify miRNA targets, degradome library was sequenced and the bioinformatic analysis confirmed 53 genes as targets for miRNAs and Tas3-siRNAs. Degradome analysis also confirmed a conserved fine-tuning mechanism of Tas3-siRNA abundance by siRNA-mediated silencing of TAS3 transcripts in diverse plant species. These findings added additional information to the small RNA knowledge-base in wheat.

Identification of microRNA targets in tomato fruit development using high-throughput sequencing and degradome analysis

MicroRNAs (miRNAs) play important roles in plant development through regulation of gene expression by mRNA degradation or translational inhibition. Despite the fact that tomato (Solanum lycopersicum) is the model system for studying fleshy fruit development and ripening, only a few experimentally proven miRNA targets are known, and the role of miRNA action in these processes remains largely unknown. Here, by using parallel analysis of RNA ends (PARE) for global identification of miRNA targets and comparing four different stages of tomato fruit development, a total of 119 target genes of miRNAs were identified. Of these, 106 appeared to be new targets. A large part of the identified targets (56) coded for transcription factors. Auxin response factors, as well as two known ripening regulators, colorless non-ripening (CNR) and APETALA2a (SlAP2a), with developmentally regulated degradation patterns were identified. The levels of the intact messenger of both CNR and AP2a are actively modulated during ripening, by miR156/157 and miR172, respectively. Additionally, two TAS3-mRNA loci were identified as targets of miR390. Other targets such as Argonaute 1 (AGO1), shown to be involved in miRNA biogenesis in other plant species, were identified, which suggests a feedback loop regulation of this process. In this study, it is shown that miRNA-guided cleavage of mRNAs is likely to play an important role in tomato fruit development and ripening.

Boron stress responsive microRNAs and their targets in barley

Boron stress is an environmental factor affecting plant development and production. Recently, microRNAs (miRNAs) have been found to be involved in several plant processes such as growth regulation and stress responses. In this study, miRNAs associated with boron stress were identified and characterized in barley. miRNA profiles were also comparatively analyzed between root and leave samples. A total of 31 known and 3 new miRNAs were identified in barley; 25 of them were found to respond to boron treatment. Several miRNAs were expressed in a tissue specific manner; for example, miR156d, miR171a, miR397, and miR444a were only detected in leaves. Additionally, a total of 934 barley transcripts were found to be specifically targeted and degraded by miRNAs. In silico analysis of miRNA target genes demonstrated that many miRNA targets are conserved transcription factors such as Squamosa promoter-binding protein, Auxin response factor (ARF), and the MYB transcription factor family. A majority of these targets were responsible for plant growth and response to environmental changes. We also propose that some of the miRNAs in barley such as miRNA408 might play critical roles against boron exposure. In conclusion, barley may use several pathways and cellular processes targeted by miRNAs to cope with boron stress.

Evolution of variety-specific regulatory schema for expression of osa-miR408 in indica rice varieties under drought stress

Evolution of differential regulatory mechanisms can lead to quite distinct physiological attributes. In the present study, we have identified one such regulatory schema that regulates osa-miR408 and responds differentially in drought-sensitive and -tolerant indica rice varieties. A comparison of the drought stress response in drought-sensitive (Pusa Basmati 1 and IR64) and drought-tolerant (Nagina 22 and Vandana) indica rice varieties revealed that, during drought stress, levels of miR408 transcript decrease significantly in sensitive cultivars, whereas they remain elevated in the tolerant cultivars. The trend is reflected in young seedlings, as well as in flag leaf and spikelets of adult plants (heading stage). Members of the plastocyanin-like protein family targeted by miR408 also show the inverse expression profile and thus accumulate at a lower level in tolerant cultivars during drought. Interestingly, some members of this family are implicated in maintaining the cellular redox state and spikelet fertility in Arabidopsis. An investigation of miR408 loci (including promoter) in all four cultivars did not reveal any significant sequence variation indicating an involvement of the upstream regulatory schema. Indeed, a similar variety-specific stress response was found in the Oryza sativa squamosa promoter-binding-like 9 transcription factor that regulates miR408 expression. We further demonstrate that drought-mediated induction of miR408 in Nagina 22 is regulated by [Ca(2+)]cyt levels. However, [Ca(2+)]cyt does not appear to regulate miR408 levels in Pusa Basmati 1, suggesting a variety-specific evolution of regulatory schema in rice.

Identification and characterization of novel and conserved microRNAs in radish (Raphanus sativus L.) using high-throughput sequencing

MicroRNAs (miRNAs) are endogenous, non-coding, small RNAs that play significant regulatory roles in plant growth, development, and biotic and abiotic stress responses. To date, a great number of conserved and species-specific miRNAs have been identified in many important plant species such as Arabidopsis, rice and poplar. However, little is known about identification of miRNAs and their target genes in radish (Raphanus sativus L.). In the present study, a small RNA library from radish root was constructed and sequenced using the high-throughput Solexa sequencing. Through sequence alignment and secondary structure prediction, a total of 545 conserved miRNA families as well as 15 novel (with their miRNA* strand) and 64 potentially novel miRNAs were identified. Quantitative real-time PCR (qRT-PCR) analysis confirmed that both conserved and novel miRNAs were expressed in radish, and some of them were preferentially expressed in certain tissues. A total of 196 potential target genes were predicted for 42 novel radish miRNAs. Gene ontology (GO) analysis showed that most of the targets were involved in plant growth, development, metabolism and stress responses. This study represents a first large-scale identification and characterization of radish miRNAs and their potential target genes. These results could lead to the further identification of radish miRNAs and enhance our understanding of radish miRNA regulatory mechanisms in diverse biological and metabolic processes.

Transposable element-associated microRNA hairpins produce 21-nt sRNAs integrated into typical microRNA pathways in rice

microRNAs (miRNAs) are a class of small RNAs (sRNAs) of ~21 nucleotides (nt) in length processed from foldback hairpins by dicer-like1 (DCL1) or DCL4. They regulate the expression of target mRNAs by base pairing through RNA-induced silencing complex (RISC). In the RISC, Argonaute1 (AGO1) is the key protein that cleaves miRNA targets at position ten of a miRNA:target duplex. The authenticity of many annotated rice miRNA hairpins is under debate because of their homology to repeat sequences. Some of them, like miR1884b, have been removed from the current release of miRBase based on incomplete information. In this study, we investigated the association of transposable element (TE)-derived miRNAs with typical miRNA pathways (DCL1/4- and AGO1-dependent) using publicly available deep sequencing datasets. Seven miRNA hairpins with 13 unique sRNAs were specifically enriched in AGO1 immunoprecipitation samples and relatively reduced in DCL1/4 knockdown genotypes. Interestingly, these species are ~21-nt long, instead of 24-nt as annotated in miRBase and the literature. Their expression profiles meet current criteria for functional annotation of miRNAs. In addition, diagnostic cleavage tags were found in degradome datasets for predicted target mRNAs. Most of these miRNA hairpins share significant homology with miniature inverted-repeat transposable elements, one type of abundant DNA transposons in rice. Finally, the root-specific production of a 24-nt miRNA-like sRNA was confirmed by RNA blot for a novel EST that maps to the 3'-UTR of a candidate pseudogene showing extensive sequence homology to miR1884b hairpin. Our data are consistent with the hypothesis that TEs can serve as a driving force for the evolution of some MIRNAs, where co-opting of DICER-LIKE1/4 processing and integration into AGO1 could exapt transcribed TE-associated hairpins into typical miRNA pathways.

[Non-coding RNAs involved in plant responses to environmental constraints]

In recent years, in addition to mRNAs, the non-protein-coding RNAs (or ncRNAs) have emerged as a major part of the eukaryotic transcriptome. New genomic approaches allowed the discovery of many novel long and small ncRNAs that may be linked to the generation of evolutionary complexity in multicellular organisms. Many long ncRNAs are regulated by abiotic stresses although only very few long ncRNAs have been functionally analyzed. On the other hand, small RNAs act in the regulation of gene expression at transcriptional or post-transcriptional level and several among them have been linked to abiotic stress responses. Here we describe various ncRNAs associated with environmental stress responses such as to salt, cold or nutrient deprivation. The understanding of these RNA networks may reveal novel mechanisms involved in plant adaptation to changing environmental conditions.

Characterization and expression profiling of selected microRNAs in tomato (Solanum lycopersicon) 'Jiangshu14'

Presence of selected tomato (Solanum lycopersicon) microRNAs (sly-miRNAs) was validated and their expression profiles established in roots, stems, leaves, flowers and fruits of tomato variety Jiangshu14 by quantitative RT-PCR (qRT-PCR). In addition conservation characteristics these sly-miRNAs were analyzed and target genes predicted bioinformatically. Results indicate that some of these miRNAs are specific to tomato while most are conserved in other plant species. Predicted sly-miRNA targets genes were shown to be targeted by either by a single or more miRNAs and are involved in diverse processes in tomato plant growth and development. All the 36 miRNAs were present in the cDNA of mixed tissues and qRT-PCR revealed that some of these sly-miRNAs are ubiquitous in tomato while others have tissue-specific expression. The experimental validation and expression profiling as well target gene prediction of these miRNAs in tomato as done in this study can add to the knowledge on the important roles played by these sly-miRNAs in the growth and development, environmental stress tolerance as well as pest and disease resistance in tomatoes and related species. In addition these findings broaden the knowledge of small RNA-mediated regulation in S. lycopersicon. It is recommended that experimental validation of the target genes be done so as to give a much more comprehensive information package on these miRNAs in tomato and specifically in the selected variety.

Analyses of a Glycine max degradome library identify microRNA targets and microRNAs that trigger secondary siRNA biogenesis

Plant microRNAs (miRNAs) regulate gene expression mainly by guiding cleavage of target mRNAs. In this study, a degradome library constructed from different soybean (Glycine max (L.) Merr.) tissues was deep-sequenced. 428 potential targets of small interfering RNAs and 25 novel miRNA families were identified. A total of 211 potential miRNA targets, including 174 conserved miRNA targets and 37 soybean-specific miRNA targets, were identified. Among them, 121 targets were first discovered in soybean. The signature distribution of soybean primary miRNAs (pri-miRNAs) showed that most pri-miRNAs had the characteristic pattern of Dicer processing. The biogenesis of TAS3 small interfering RNAs (siRNAs) was conserved in soybean, and nine Auxin Response Factors were identified as TAS3 siRNA targets. Twenty-three miRNA targets produced secondary small interfering RNAs (siRNAs) in soybean. These targets were guided by five miRNAs: gma-miR393, gma-miR1508, gma-miR1510, gma-miR1514, and novel-11. Multiple targets of these secondary siRNAs were detected. These 23 miRNA targets may be the putative novel TAS genes in soybean. Global identification of miRNA targets and potential novel TAS genes will contribute to research on the functions of miRNAs in soybean.

Genome-wide identification of soybean microRNAs and their targets reveals their organ-specificity and responses to phosphate starvation

BACKGROUND

Phosphorus (P) plays important roles in plant growth and development. MicroRNAs involved in P signaling have been identified in Arabidopsis and rice, but P-responsive microRNAs and their targets in soybean leaves and roots are poorly understood.

RESULTS

Using high-throughput sequencing-by-synthesis (SBS) technology, we sequenced four small RNA libraries from leaves and roots grown under phosphate (Pi)-sufficient (+Pi) and Pi-depleted (-Pi) conditions, respectively, and one RNA degradome library from Pi-depleted roots at the genome-wide level. Each library generated ~21.45-28.63 million short sequences, resulting in ~20.56-27.08 million clean reads. From those sequences, a total of 126 miRNAs, with 154 gene targets were computationally predicted. This included 92 new miRNA candidates with 20-23 nucleotides that were perfectly matched to the Glycine max genome 1.0, 70 of which belong to 21 miRNA families and the remaining 22 miRNA unassigned into any existing miRNA family in miRBase 18.0. Under both +Pi and -Pi conditions, 112 of 126 total miRNAs (89%) were expressed in both leaves and roots. Under +Pi conditions, 12 leaf- and 2 root-specific miRNAs were detected; while under -Pi conditions, 10 leaf- and 4 root-specific miRNAs were identified. Collectively, 25 miRNAs were induced and 11 miRNAs were repressed by Pi starvation in soybean. Then, stem-loop real-time PCR confirmed expression of four selected P-responsive miRNAs, and RLM-5' RACE confirmed that a PHO2 and GmPT5, a kelch-domain containing protein, and a Myb transcription factor, respectively are targets of miR399, miR2111, and miR159e-3p. Finally, P-responsive cis-elements in the promoter regions of soybean miRNA genes were analyzed at the genome-wide scale.

CONCLUSIONS

Leaf- and root-specific miRNAs, and P-responsive miRNAs in soybean were identified genome-wide. A total of 154 target genes of miRNAs were predicted via degradome sequencing and computational analyses. The targets of miR399, miR2111, and miR159e-3p were confirmed. Taken together, our study implies the important roles of miRNAs in P signaling and provides clues for deciphering the functions for microRNA/target modules in soybean.

Characterization and expression patterns of microRNAs involved in rice grain filling

MicroRNAs (miRNAs) are upstream gene regulators of plant development and hormone homeostasis through their directed cleavage or translational repression of the target mRNAs, which may play crucial roles in rice grain filling and determining the final grain weight and yield. In this study, high-throughput sequencing was performed to survey the dynamic expressions of miRNAs and their corresponding target genes at five distinct developmental stages of grain filling. In total, 445 known miRNAs and 45 novel miRNAs were detected with most of them expressed in a developmental stage dependent manner, and the majority of known miRNAs, which increased gradually with rice grain filling, showed negatively related to the grain filling rate. Detailed expressional comparisons revealed a clear negative correlation between most miRNAs and their target genes. It was found that specific miRNA cohorts are expressed in a developmental stage dependent manner during grain filling and the known functions of these miRNAs are involved in plant hormone homeostasis and starch accumulation, indicating that the expression dynamics of these miRNAs might play key roles in regulating rice grain filling.

Identification of miRNAs and their targets using high-throughput sequencing and degradome analysis in cytoplasmic male-sterile and its maintainer fertile lines of Brassica juncea

Background

Regulatory network of cytoplasmic male sterility (CMS) occurrence is still largely unknown in plants, although numerous researches have been attempted to isolate genes involved in CMS. Here, we employed high-throughput sequencing and degradome analysis to identify microRNAs and their targets using high-throughput sequencing in CMS and its maintainer fertile (MF) lines of Brassica juncea.

Results

We identified 197 known and 78 new candidate microRNAs during reproductive development of B. juncea. A total of 47 differentially expressed microRNAs between CMS and its MF lines were discovered, according to their sequencing reads number. Different expression levels of selected microRNAs were confirmed by using real-time quantitative PCR between CMS and MF lines. Furthermore, we observed that the transcriptional patterns of these microRNAs could be mimicked by artificially inhibiting mitochondrial F₁F₀-ATPase activity and its function in MF line by using treatment with oligomycin. Targeted genes of the microRNAs were identified by high-throughput sequencing and degradome approaches, including auxin response factor, NAC (No Apical Meristem) domain transcription factor, GRAS family transcription factor, MYB transcription factor, squamosa promoter binding protein, AP2-type transcription factor, homeobox/homeobox-leucine zipper family and TCP family transcription factors, which were observed to be differentially expressed between CMS and MF.

Conclusion

Taken together, from these findings we suggested microRNA might participate in the regulatory network of CMS by tuning fork in gene expressions in CMS B. juncea. The differential expression of miRNAs observed between CMS and MF lines suggested that biogenesis of miRNAs could be influenced in the CMS.

MicroRNA target prediction and validation

The accurate prediction and validation of microRNA targets is essential to understanding the function of microRNAs. Computational predictions indicate that all human genes may be regulated by microRNAs, with each microRNA possibly targeting thousands of genes. Here we discuss computational and experimental methods for identifying mammalian microRNA targets. We describe microRNA target prediction resources and procedures that are suitable for experiments where more accurate prediction of microRNA targets is more important than detecting all putative targets. We then discuss experimental methods for identifying and validating microRNA target genes, with an emphasis on the target reporter assay as the method of choice for specifically testing functional microRNA target sites.

MicroRNAs play critical roles during plant development and in response to abiotic stresses

MicroRNAs (miRNAs) have been identified as key molecules in regulatory networks. The fine-tuning role of miRNAs in addition to the regulatory role of transcription factors has shown that molecular events during development are tightly regulated. In addition, several miRNAs play crucial roles in the response to abiotic stress induced by drought, salinity, low temperatures, and metals such as aluminium. Interestingly, several miRNAs have overlapping roles with regard to development, stress responses, and nutrient homeostasis. Moreover, in response to the same abiotic stresses, different expression patterns for some conserved miRNA families among different plant species revealed different metabolic adjustments. The use of deep sequencing technologies for the characterisation of miRNA frequency and the identification of new miRNAs adds complexity to regulatory networks in plants. In this review, we consider the regulatory role of miRNAs in plant development and abiotic stresses, as well as the impact of deep sequencing technologies on the generation of miRNA data.

Identification and comparative analysis of microRNAs associated with low-N tolerance in rice genotypes

BACKGROUND

Nitrogen [N] is a critical limiting nutrient for plants and has to be exogenously supplied to many crops, to achieve high yield with significant economic and environmental costs, specifically for rice. Development of low-input nitrogen sustainable crop is necessary for sustainable agriculture. Identification of regulatory elements associated with low-N tolerance is imperative for formulating innovative approaches for developing low-N tolerant crop plants, using gene manipulation. MicroRNAs (miRNAs) are known to play crucial roles in the modulation of gene expression in plants under various environmental conditions.

METHODOLOGY/PRINCIPAL FINDINGS

MiRNAs associated with low-N tolerance have not been identified so far. In this study, we investigated microarray-based miRNA expression in low-N tolerant and low-N sensitive rice genotypes under low N condition. Expressions of 32 miRNAs differed significantly in the two genotypes. Of these 32 miRNAs, expressions of nine miRNAs were further validated experimentally in leaves as well as in roots. Of these differentially expressed miRNAs, six miRNAs (miR156, miR164, miR528, miR820, miR821 and miR1318) were reported in leaves and four (miR164, miR167, miR168 and miR528) in roots. Target genes of all the 32 miRNAs were predicted, which encode transcription factors, and proteins associated with metabolic processes or stress responses. Expression levels of some of the corresponding miRNA targets were analysed and found to be significantly higher in low N-tolerant genotype than low-N sensitive genotype. These findings suggested that miRNAs played an important role in low-N tolerance in rice.

CONCLUSIONS/SIGNIFICANCE

Genome-wide differences in expression of miRNA in low N-tolerant and low N-sensitive rice genotypes were reported. This provides a platform for selection as well as manipulation of genotypes for better N utilization efficiency.

Identification of Taxus microRNAs and their targets with high-throughput sequencing and degradome analysis

Plant microRNAs (miRNAs) have an impact in the regulation of several biological processes such as development, growth and metabolism by negatively controlling gene expression at the post-transcriptional level. However, the role of these small molecules in the medicinal gymnosperm species Taxus remained elusive. To elucidate the role of miRNAs in Taxus we used a deep sequencing approach to analyze the small RNA and degradome sequence tags of Taxus mairei leaves. For miRNAs, the sequencing library generated 14.9 million short sequences, resulting in 13.1 million clean reads. The library contains predominantly small RNAs with 21 nucleotide length, followed by 19-nt and 20-nt small RNAs. Around 29% of total small RNAs are matched to the T. mairei transcriptome. By sequence alignment, we identified 871 mature miRNAs, 15 miRNA* and 869 miRNA precursors representing known plant miRNA families. There are 547 unique small RNA matching the miRNA precursors. We predict 37 candidate novel miRNAs from the unannotated small RNAs that could be mapped to the reference transcriptome. The expression of the selected candidates was for the first time quantified by real-time reverse transcription polymerase chain reaction. The novel miRNA m0034 turns out to be from the intron sequence of the paclitaxel biosynthetic gene taxadiene synthase. The 21 potential targets of nine novel miRNAs are also predicted. Additionally, 56 targets for known miRNA families and 15 targets for novel candidate miRNA families were identified by high-throughput degradome-sequencing approach. It is found that two paclitaxel biosynthetic genes, taxane 13α hydroxylase and taxane 2α-O-benzoyltransferase, are the cleavage targets of miR164 and miR171, respectively. This study represents the first transcriptome-based analysis of miRNAs and degradome in gymnosperms.

Identification of the highly accumulated microRNA*s in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa)

Plant microRNAs (miRNAs) are crucial for the regulation of gene expression, which is involved in almost all the important biological processes. In the cytoplasm, the miRNA strand is selectively incorporated into a specific Argonaute (AGO)-associated gene silencing complex, while the miRNA* is degraded rapidly. Thus, most miRNA*s were thought to be biologically meaningless. Interestingly, several recent reports in both plants and animals have shaken this notion. Many miRNA*s were demonstrated to possess regulatory roles in gene expression. However, the low accumulation levels of most miRNA*s raise the question whether the activities of this small RNA (sRNA) species are widespread in plants. Here, by using publicly available sRNA high-throughput sequencing data, we found that the accumulation levels of several miRNA*s could be much higher than those of their miRNA partners in certain organs, mutants and/or AGO-associated silencing complexes of both Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). Based on target prediction and degradome sequencing data-based validation, some of these highly accumulated miRNA*s were indicated to possess cleavage-based potential regulatory role on certain targets. Besides, some interesting biological interpretations were obtained based on the accumulation patterns of the miRNA*s, the annotations of the target genes, and literature mining. Taken together, the expanded list of the highly accumulated miRNA*s along with their potential target genes discovered in this study further strengthened the current notion that certain members of the miRNA* species are biologically relevant, which needs further inspection.

Identification of nitrogen starvation-responsive microRNAs in Arabidopsis thaliana

microRNAs (miRNAs) are a class of negative regulators that take part in many processes such as growth and development, stress responses, and metabolism in plants. Recently, miRNAs were shown to function in plant nutrient metabolism. Moreover, several miRNAs were identified in the response to nitrogen (N) deficiency. To investigate the functions of other miRNAs in N deficiency, deep sequencing technology was used to detect the expression of small RNAs under N-sufficient and -deficient conditions. The results showed that members from the same miRNA families displayed differential expression in response to N deficiency. Upon N starvation, the expression of miR169, miR171, miR395, miR397, miR398, miR399, miR408, miR827, and miR857 was repressed, whereas those of miR160, miR780, miR826, miR842, and miR846 were induced. miR826, a newly identified N-starvation-induced miRNA, was found to target the AOP2 gene. Among these N-starvation-responsive miRNAs, several were involved in cross-talk among responses to different nutrient (N, P, S, Cu) deficiencies. miR160, miR167, and miR171 could be responsible for the development of Arabidopsis root systems under N-starvation conditions. In addition, twenty novel miRNAs were identified and nine of them were significantly responsive to N-starvation. This study represents comprehensive expression profiling of N-starvation-responsive miRNAs and advances our understanding of the regulation of N homeostasis mediated by miRNAs.

Genomic dissection of small RNAs in wild rice (Oryza rufipogon): lessons for rice domestication

The lack of a MIRNA set and genome sequence of wild rice (Oryza rufipogon) has prevented us from determining the role of MIRNA genes in rice domestication. In this study, a genome, three small RNA populations and a degradome of O. rufipogon were sequenced by Illumina platform and the expression levels of microRNAs (miRNAs) were investigated by miRNA chips. A de novo O. rufipogon genome was assembled using c. 55× coverage of raw sequencing data and a total of 387 MIRNAs were identified in the O. rufipogon genome based on c. 5.2 million unique small RNA reads from three different tissues of O. rufipogon. Of these, O. rufipogon MIRNAs, 259 were not found in the cultivated rice, suggesting a loss of these MIRNAs in the cultivated rice. We also found that 48 MIRNAs were novel in the cultivated rice, suggesting that they were potential targets of domestication selection. Some miRNAs showed significant expression differences between wild and cultivated rice, suggesting that expression of miRNA could also be a target of domestication, as demonstrated for the miR164 family. Our results illustrated that MIRNA genes, like protein-coding genes, might have been significantly shaped during rice domestication and could be one of the driving forces that contributed to rice domestication.

Finding microRNA targets in plants: current status and perspectives

MicroRNAs (miRNAs), a class of ∼20–24 nt long non-coding RNAs, have critical roles in diverse biological processes including development, proliferation, stress response, etc. With the development and availability of experimental technologies and computational approaches, the field of miRNA biology has advanced tremendously over the last decade. By sequence complementarity, miRNAs have been estimated to regulate certain mRNA transcripts. Although it was once thought to be simple and straightforward to find plant miRNA targets, this viewpoint is being challenged by genetic and biochemical studies. In this review, we summarize recent progress in plant miRNA target recognition mechanisms, principles of target prediction, and introduce current experimental and computational tools for plant miRNA target prediction. At the end, we also present our thinking on the outlook for future directions in the development of plant miRNA target finding methods.

Identification of wild soybean miRNAs and their target genes responsive to aluminum stress

BACKGROUND

MicroRNAs (miRNAs) play important regulatory roles in development and stress response in plants. Wild soybean (Glycine soja) has undergone long-term natural selection and may have evolved special mechanisms to survive stress conditions as a result. However, little information about miRNAs especially miRNAs responsive to aluminum (Al) stress is available in wild soybean.

RESULTS

Two small RNA libraries and two degradome libraries were constructed from the roots of Al-treated and Al-free G. soja seedlings. For miRNA identification, a total of 7,287,655 and 7,035,914 clean reads in Al-treated and Al-free small RNAs libraries, respectively, were generated, and 97 known miRNAs and 31 novel miRNAs were identified. In addition, 49 p3 or p5 strands of known miRNAs were found. Among all the identified miRNAs, the expressions of 30 miRNAs were responsive to Al stress. Through degradome sequencing, 86 genes were identified as targets of the known miRNAs and five genes were found to be the targets of the novel miRNAs obtained in this study. Gene ontology (GO) annotations of target transcripts indicated that 52 target genes cleaved by conserved miRNA families might play roles in the regulation of transcription. Additionally, some genes, such as those for the auxin response factor (ARF), domain-containing disease resistance protein (NB-ARC), leucine-rich repeat and toll/interleukin-1 receptor-like protein (LRR-TIR) domain protein, cation transporting ATPase, Myb transcription factors, and the no apical meristem (NAM) protein, that are known to be responsive to stress, were found to be cleaved under Al stress conditions.

CONCLUSIONS

A number of miRNAs and their targets were detected in wild soybean. Some of them that were responsive to biotic and abiotic stresses were regulated by Al stress. These findings provide valuable information to understand the function of miRNAs in Al tolerance.

Distinctive expression patterns and roles of the miRNA393/TIR1 homolog module in regulating flag leaf inclination and primary and crown root growth in rice (Oryza sativa)

• MicroRNA (miRNA)-mediated regulation of auxin signaling components plays a critical role in plant development. miRNA expression and functional diversity contribute to the complexity of regulatory networks of miRNA/target modules. • This study functionally characterizes two members of the rice (Oryza sativa) miR393 family and their target genes, OsTIR1 and OsAFB2 (AUXIN SIGNALING F-BOX), the two closest homologs of Arabidopsis TRANSPORT INHIBITOR RESPONSE 1 (TIR1). • We found that the miR393 family members possess distinctive expression patterns, with miR393a expressed mainly in the crown and lateral root primordia, as well as the coleoptile tip, and miR393b expressed in the shoot apical meristem. Transgenic plants overexpressing miR393a/b displayed a severe phenotype with hallmarks of altered auxin signaling, mainly including enlarged flag leaf inclination and altered primary and crown root growth. Furthermore, OsAFB2- and OsTIR1-suppressed lines exhibited increased inclination of flag leaves at the booting stage, resembling miR393-overexpressing plants. Moreover, yeast two-hybrid and bimolecular fluorescence complementation assays showed that OsTIR1 and OsAFB2 interact with OsIAA1. • Expression diversification of miRNA393 implies the potential role of miRNA regulation during species evolution. The conserved mechanisms of the miR393/target module indicate the fundamental importance of the miR393-mediated regulation of auxin signal transduction in rice.

Transcriptome-wide characterization of miRNA-directed and non-miRNA-directed endonucleolytic cleavage using Degradome analysis under low ambient temperature in Phalaenopsis aphrodite subsp. formosana

Plant microRNAs (miRNAs) regulate gene expression through post-transcriptional gene silencing. Phalaenopsis aphrodite subsp. formosana is an orchid species native to Taiwan, which has high economic value and a high frequency of floral polymorphism. To date, few studies have focused on the regulatory roles of miRNAs and functional small RNAs (sRNAs) in orchids although understanding the regulation of flower development and flowering time is potentially important. Here, we combined analyses of the transcriptome, sRNAs and the degradome to identify sRNA-directed transcript cleavages in Phalaenopsis. Degradome analysis provided large-scale evidence of conserved and novel miRNA-directed cleavage of target transcripts, and 46 abundant sRNA groups and their target transcripts were identified. Low temperature-responsive sRNAs were validated with normalized reads from an sRNA library and quantitative stem-loop reverse transcription-PCR (RT-PCR) analysis. According to gene ontology (GO) categorization, target transcripts of the novel miRNAs and sRNAs are functionally involved in metabolic processes or responses to stress. One particular homologous gene, Allcontig28452, which encodes digalactosyldiacylglycerol synthase 2 (DGD2), was found to be targeted by natural antisense transcripts (NATs) unique to Phalaenopsis. In summary, comprehensive analyses of the transcriptome, sRNAs and degradome using deep sequencing technology provided a useful platform for investigating miRNA-directed and non-miRNA-directed endonucleolytic cleavage in a non-model plant, the orchid Phalaenopsis.

Exploring the evolutionary differences of SBP-box genes targeted by miR156 and miR529 in plants

The combinatorial control of one target by multiple miRNAs brings big challenges to elucidate its precise evolutionary mechanism. Squamosa promoter binding protein-like (SBP) gene family exhibits the different regulatory patterns, in which some members are only regulated by miR156 and others by miR156 and miR529. Here, we explored the different evolutionary patterns and rates between miR156 targets and miR529 ones in three species (moss, rice, and maize). Our work found that the miR529 targets were members of miR156 target dataset, indicative of cooperative control. Further phylogenetic analyses as well as gene structure features demonstrated that miR529 targets derived from a monophyletic branch of miR156 targets which evolved into two independent branches duo to the ancient gene duplication. Moreover, inspection of evolutionary rate parameters (d(N)/d(S), d(N) and d(S)) for miR156 targets and miR529 ones revealed they were under different selection strength. MiR529 targets were more constraint by strong purifying selection and evolved conservatively with a slow rate. By contrast, miR156 targets evolved more rapidly and experienced more relaxed purifying selection, which may contribute to their functional diversification. Our results will enhance the understanding of different evolutionary fates of SBP-box genes regulated by the different numbers of miRNA families before functional studies.

MicroRNA profiles and their control of male gametophyte development in rice

Plant microRNAs (miRNAs) act as negative regulators of gene expression by slicing target transcripts or inhibiting translation. A number of miRNAs play important roles in development. In order to investigate the potential function of miRNAs during male gametogenesis in rice, we obtained both gene and small RNA expression profiles by combining microarray and high-throughput sequencing technologies. From the microarray datasets, 2,925 male gametophyte-specific genes were identified, including 107 transcription factors and three significant Argonaute genes (AGO12, AGO13, and AGO17). From the sRNA-Seq datasets, 104 unique miRNAs (miRus) were identified, including 47 known miRus and 57 novel miRus; interestingly, most of the new miRus are pollen-specific and not conserved among species. Furthermore, an interactive network of miRNA-target was constructed based on the two datasets. By employing enrichment analysis, the miRNA-regulated targets were found to be involved in both the up and down pathways, but predominantly in the down pathways, including 37 GO biological processes and 32 KEGG pathways. These findings indicate that miRNAs play a broad regulatory role during male gametophyte development in rice.

Conservation and divergence in plant microRNAs

MicroRNAs (miRNAs) are a class of small, non-coding RNAs that regulate gene expression in eukaryotic cells. The past decade has seen an explosion in our understanding of the sets of miRNA genes encoded in the genomes in different species of plants and the mechanisms by which miRNAs interact with target RNAs. A subset of miRNA families (and their binding sites in target RNAs) are conserved between angiosperms and basal plants, suggesting they predate the divergence of existing lineages of plants. However, the majority of miRNA families expressed by any given plant species have a narrow phylogenetic distribution. As a group, these "young" miRNAs genes appear to be evolutionarily fluid and lack clearly understood biological function. The goal of this review is to summarize our understanding of the sets of miRNA genes and miRNA targets that exist in various plant species and to discuss hypotheses that explain the patterns of conservation and divergence observed among microRNAs in plants.

Identification of miRNAs and their targets from Brassica napus by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are endogenous regulators of a broad range of physiological processes and act by either degrading mRNA or blocking its translation. Oilseed rape (Brassica napus) is one of the most important crops in China, Europe and other Asian countries with publicly available expressed sequence tags (ESTs) and genomic survey sequence (GSS) databases, but little is known about its miRNAs and their targets. To date, only 46 miRNAs have been identified in B. napus.

RESULTS

Forty-one conserved and 62 brassica-specific candidate B. napus miRNAs, including 20 miRNA* sequences, were identified using Solexa sequencing technology. Furthermore, 33 non-redundant mRNA targets of conserved brassica miRNAs and 19 new non-redundant mRNA targets of novel brassica-specific miRNAs were identified by genome-scale sequencing of mRNA degradome.

CONCLUSIONS

This study describes large scale cloning and characterization of B. napus miRNAs and their potential targets, providing the foundation for further characterization of miRNA function in the regulation of diverse physiological processes in B. napus.

High throughput sequencing reveals novel and abiotic stress-regulated microRNAs in the inflorescences of rice

BACKGROUND

MicroRNAs (miRNAs) are small RNA molecules that play important regulatory roles in plant development and stress responses. Identification of stress-regulated miRNAs is crucial for understanding how plants respond to environmental stimuli. Abiotic stresses are one of the major factors that limit crop growth and yield. Whereas abiotic stress-regulated miRNAs have been identified in vegetative tissues in several plants, they are not well studied in reproductive tissues such as inflorescences.

RESULTS

We used Illumina deep sequencing technology to sequence four small RNA libraries that were constructed from the inflorescences of rice plants that were grown under control condition and drought, cold, or salt stress. We identified 227 miRNAs that belong to 127 families, including 70 miRNAs that are not present in the miRBase. We validated 62 miRNAs (including 10 novel miRNAs) using published small RNA expression data in DCL1, DCL3, and RDR2 RNAi lines and confirmed 210 targets from 86 miRNAs using published degradome data. By comparing the expression levels of miRNAs, we identified 18, 15, and 10 miRNAs that were regulated by drought, cold and salt stress conditions, respectively. In addition, we identified 80 candidate miRNAs that originated from transposable elements or repeats, especially miniature inverted-repeat elements (MITEs).

CONCLUSION

We discovered novel miRNAs and stress-regulated miRNAs that may play critical roles in stress response in rice inflorescences. Transposable elements or repeats, especially MITEs, are rich sources for miRNA origination.

miRNA regulation in the early development of barley seed

BACKGROUND

During the early stages of seed development many genes are under dynamic regulation to ensure the proper differentiation and establishment of the tissue that will constitute the mature grain. To investigate how miRNA regulation contributes to this process in barley, a combination of small RNA and mRNA degradome analyses were used to identify miRNAs and their targets.

RESULTS

Our analysis identified 84 known miRNAs and 7 new miRNAs together with 96 putative miRNA target genes regulated through a slicing mechanism in grain tissues during the first 15 days post anthesis. We also identified many potential miRNAs including several belonging to known miRNA families. Our data gave us evidence for an increase in miRNA-mediated regulation during the transition between pre-storage and storage phases. Potential miRNA targets were found in various signalling pathways including components of four phytohormone pathways (ABA, GA, auxin, ethylene) and the defence response to powdery mildew infection. Among the putative miRNA targets we identified were two essential genes controlling the GA response, a GA3oxidase1 and a homolog of the receptor GID1, and a homolog of the ACC oxidase which catalyses the last step of ethylene biosynthesis. We found that two MLA genes are potentially miRNA regulated, establishing a direct link between miRNAs and the R gene response.

CONCLUSION

Our dataset provides a useful source of information on miRNA regulation during the early development of cereal grains and our analysis suggests that miRNAs contribute to the control of development of the cereal grain, notably through the regulation of phytohormone response pathways.

A reversed framework for the identification of microRNA-target pairs in plants

Most plant microRNAs (miRNAs) perform their repressive regulation through target cleavages. The resulting slicing sites on the target transcripts could be mapped by sequencing of the 3'-cleavage remnants, called degradome sequencing. The high sequence complementarity between miRNAs and their targets has greatly facilitated the development of the target prediction tools for plant miRNAs. The prediction results were then subjected to degradome sequencing data-based validation, through which numerous miRNA-target interactions have been extracted. However, some drawbacks are unavoidable when using this forward approach. Essentially, a known list of plant miRNAs should be obtained in advance of target prediction and validation. This becomes an obstacle to discover novel miRNAs and their targets. Here, after reviewing the current available algorithms for reverse identification of miRNA-target pairs in plants, a case study was performed by using a newly established framework with adjustable parameters. In this workflow, integration of degradome and ARGONAUTE 1-enriched small RNA sequencing data was recommended to do a relatively comprehensive and reliable search. Besides, several computational algorithms such as BLAST, target plots and RNA secondary structure prediction were used. The results demonstrated the prevalent utility of the reversed approach for uncovering miRNA-target interactions in plants.

Identification of soybean seed developmental stage-specific and tissue-specific miRNA targets by degradome sequencing

BACKGROUND

MicroRNAs (miRNAs) regulate the expression of target genes by mediating gene silencing in both plants and animals. The miRNA targets have been extensively investigated in Arabidopsis and rice using computational prediction, experimental validation by overexpression in transgenic plants, and by degradome or PARE (parallel analysis of RNA ends) sequencing. However, miRNA targets mostly remain unknown in soybean (Glycine max). More specifically miRNA mediated gene regulation at different seed developmental stages in soybean is largely unexplored. In order to dissect miRNA guided gene regulation in soybean developing seeds, we performed a transcriptome-wide experimental method using degradome sequencing to directly detect cleaved miRNA targets.

RESULTS

In this study, degradome libraries were separately prepared from immature soybean cotyledons representing three stages of development and from seed coats of two stages. Sequencing and analysis of 10 to 40 million reads from each library resulted in identification of 183 different targets for 53 known soybean miRNAs. Among these, some were found only in the cotyledons representing cleavage by 25 miRNAs and others were found only in the seed coats reflecting cleavage by 12 miRNAs. A large number of targets for 16 miRNAs families were identified in both tissues irrespective of the stage. Interestingly, we identified more miRNA targets in the desiccating cotyledons of late seed maturation than in immature seed. We validated four different auxin response factor genes as targets for gma-miR160 via RNA ligase mediated 5' rapid amplification of cDNA ends (RLM-5'RACE). Gene Ontology (GO) analysis indicated the involvement of miRNA target genes in various cellular processes during seed development.

CONCLUSIONS

The miRNA targets in both the cotyledons and seed coats of several stages of soybean seed development have been elucidated by experimental evidence from comprehensive, high throughput sequencing of the enriched fragments resulting from miRNA-guided cleavage of messenger RNAs. Nearly 50% of the miRNA targets were transcription factors in pathways that are likely important in setting or maintaining the developmental program leading to high quality soybean seeds that are one of the dominant sources of protein and oil in world markets.

Genome-wide identification of Brassica napus microRNAs and their targets in response to cadmium

MicroRNAs (miRNAs) are a distinct class of small RNAs in plants that not only regulate biological processes but also regulate response to environmental stresses. The toxic heavy metal cadmium (Cd) induces expression of several miRNAs in rapeseed (Brassica napus), but it is not known on a genome-wide scale how the expression of miRNAs and their target genes, is regulated by Cd. In this study, four small RNA libraries and four degradome libraries were constructed from Cd-treated and non-Cd-treated roots and shoots of B. napus seedlings. Using high-throughput sequencing, the study identified 84 conserved and non-conserved miRNAs (belonging to 37 miRNA families) from Cd-treated and non-treated B. napus, including 19 miRNA members that were not identified before. Some of the miRNAs were validated by RNA gel blotting. Most of the identified miRNAs were found to be differentially expressed in roots/shoots or regulated by Cd exposure. The study simultaneously identified 802 targets for the 37 (24 conserved and 13 non-conserved) miRNA families, from which there are 200, 537, and 65 targets, belonging to categories I, II, and III, respectively. In category I alone, many novel targets for miRNAs were identified and shown to be involved in plant response to Cd.

Apple miRNAs and tasiRNAs with novel regulatory networks

BACKGROUND

MicroRNAs (miRNAs) and their regulatory functions have been extensively characterized in model species but whether apple has evolved similar or unique regulatory features remains unknown.

RESULTS

We performed deep small RNA-seq and identified 23 conserved, 10 less-conserved and 42 apple-specific miRNAs or families with distinct expression patterns. The identified miRNAs target 118 genes representing a wide range of enzymatic and regulatory activities. Apple also conserves two TAS gene families with similar but unique trans-acting small interfering RNA (tasiRNA) biogenesis profiles and target specificities. Importantly, we found that miR159, miR828 and miR858 can collectively target up to 81 MYB genes potentially involved in diverse aspects of plant growth and development. These miRNA target sites are differentially conserved among MYBs, which is largely influenced by the location and conservation of the encoded amino acid residues in MYB factors. Finally, we found that 10 of the 19 miR828-targeted MYBs undergo small interfering RNA (siRNA) biogenesis at the 3' cleaved, highly divergent transcript regions, generating over 100 sequence-distinct siRNAs that potentially target over 70 diverse genes as confirmed by degradome analysis.

CONCLUSIONS

Our work identified and characterized apple miRNAs, their expression patterns, targets and regulatory functions. We also discovered that three miRNAs and the ensuing siRNAs exploit both conserved and divergent sequence features of MYB genes to initiate distinct regulatory networks targeting a multitude of genes inside and outside the MYB family.

Expression and tissue-specific localization of nitrate-responsive miRNAs in roots of maize seedlings

Nitrogen availability seriously affects crop productivity and environment. The knowledge of post-transcriptional regulation of plant response to nutrients is important to improve nitrogen use efficiency of crop. This research was aimed at understanding the role of miRNAs in the molecular control of plant response to nitrate. The expression profiles of six mature miRNAs were deeply studied by quantitative real time polymerase chain reaction and in situ hybridization (ISH). To this aim, a novel optimized protocol was set up for the use of digoxygenin-labelled Zip Nucleic Acid-modified oligonucleotides as probes for ISH. Significant differences in miRNAs' transcripts accumulation were evidenced between nitrate-supplied and nitrate-depleted roots. Real-time PCR analyses and in situ detection of miRNA confirmed the array data and allowed us to evidence distinct miRNAs spatio-temporal expression patterns in maize roots. Our results suggest that a prolonged nitrate depletion may induce post-transcriptionally the expression of target genes by repressing the transcription of specific miRNAs. In particular, the repression of the transcription of miR528a/b, miR528a*/b*, miR169i/j/k, miR169i*/j*/k*, miR166j/k/n and miR408/b upon nitrate shortage could represent a crucial step integrating nitrate signals into developmental changes in maize roots.

Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing

The male sterility of thermosensitive genic male sterile (TGMS) lines of wheat (Triticum aestivum) is strictly controlled by temperature. The early phase of anther development is especially susceptible to cold stress. MicroRNAs (miRNAs) play an important role in plant development and in responses to environmental stress. In this study, deep sequencing of small RNA (smRNA) libraries obtained from spike tissues of the TGMS line under cold and control conditions identified a total of 78 unique miRNA sequences from 30 families and trans-acting small interfering RNAs (tasiRNAs) derived from two TAS3 genes. To identify smRNA targets in the wheat TGMS line, we applied the degradome sequencing method, which globally and directly identifies the remnants of smRNA-directed target cleavage. We identified 26 targets of 16 miRNA families and three targets of tasiRNAs. Comparing smRNA sequencing data sets and TaqMan quantitative polymerase chain reaction results, we identified six miRNAs and one tasiRNA (tasiRNA-ARF [for Auxin-Responsive Factor]) as cold stress-responsive smRNAs in spike tissues of the TGMS line. We also determined the expression profiles of target genes that encode transcription factors in response to cold stress. Interestingly, the expression of cold stress-responsive smRNAs integrated in the auxin-signaling pathway and their target genes was largely noncorrelated. We investigated the tissue-specific expression of smRNAs using a tissue microarray approach. Our data indicated that miR167 and tasiRNA-ARF play roles in regulating the auxin-signaling pathway and possibly in the developmental response to cold stress. These data provide evidence that smRNA regulatory pathways are linked with male sterility in the TGMS line during cold stress.

Cucumis melo microRNA expression profile during aphid herbivory in a resistant and susceptible interaction

Aphis gossypii resistance in melon (Cucumis melo) is due to the presence of a single dominant virus aphid transmission (Vat) gene belonging to the nucleotide-binding site leucine-rich repeat family of resistance genes. Significant transcriptional reprogramming occurs in Vat(+) plants during aphid infestation as metabolism shifts to respond to this biotic stress. MicroRNAs (miRNAs) are involved in the regulation of many biotic stress responses. The role of miRNAs was investigated in response to aphid herbivory during both resistant and susceptible interactions. Small RNA (smRNA) libraries were constructed from bulked leaf tissues of a Vat(+) melon line following early and late aphid infestations. Sequence analysis indicated that the expression profiles of conserved and newly identified miRNAs were altered during different stages of aphid herbivory. These results were verified by quantitative polymerase chain reaction experiments in both resistant Vat(+) and susceptible Vat(-) interactions. The comparative analyses revealed that most of the conserved miRNA families were differentially regulated during the early stages of aphid infestation in the resistant and susceptible interactions. Along with the conserved miRNA families, 18 cucurbit-specific miRNAs were expressed during the different stages of aphid herbivory. The comparison of the miRNA profiles in the resistant and susceptible interactions provides insight into the miRNA-dependent post-transcriptional gene regulation in Vat-mediated resistance.

Validation and characterization of Citrus sinensis microRNAs and their target genes

Background

MicroRNAs play vital role in plant growth and development by changeable expression of their target genes with most plant microRNAs having perfect or near-perfect complementarities with their target genes but miRNAs in Citrus sinensis (csi-miRNAs) and their function have not been widely studied.

Findings

In this study, 15 potential microRNAs in Citrus sinensis (csi-miRNAs) were identified and bioinformatically validated using miR-RACE, a newly developed method for determination of miRNAs prediction computationally. The expression of these fifteen C. sinensis miRNAs can be detected in leaves, stems, flowers and fruits of C. sinensis by QRT-PCR with some of them showed tissue-specific expression. Six potential target genes were identified for six csi-miRNAs and also experimentally verified by Poly (A) polymerase -mediated 3′ rapid amplification of cDNA ends (PPM-RACE) and RNA ligase-mediated 5′ rapid amplification of cDNA ends (RLM-RACE) which mapped the cleavage site of target mRNAs and detected expression patterns of cleaved fragments that indicate the regulatory function of the miRNAs on their target genes.

Conclusions

Our results confirm that small RNA-mediated regulation whereby all csi-miRNAs regulate their target genes by degradation.

Functions of microRNAs in plant stress responses

The discovery of microRNAs (miRNAs) as gene regulators has led to a paradigm shift in the understanding of post-transcriptional gene regulation in plants and animals. miRNAs have emerged as master regulators of plant growth and development. Evidence suggesting that miRNAs play a role in plant stress responses arises from the discovery that miR398 targets genes with known roles in stress tolerance. In addition, the expression profiles of most miRNAs that are implicated in plant growth and development are significantly altered during stress. These later findings imply that attenuated plant growth and development under stress may be under the control of stress-responsive miRNAs. Here we review recent progress in the understanding of miRNA-mediated plant stress tolerance.

A combined approach of high-throughput sequencing and degradome analysis reveals tissue specific expression of microRNAs and their targets in cucumber

MicroRNAs (miRNAs) are endogenous small RNAs playing an important regulatory function in plant development and stress responses. Among them, some are evolutionally conserved in plant and others are only expressed in certain species, tissue or developmental stages. Cucumber is among the most important greenhouse species in the world, but only a limited number of miRNAs from cucumber have been identified and the experimental validation of the related miRNA targets is still lacking. In this study, two independent small RNA libraries from cucumber leaves and roots were constructed, respectively, and sequenced with the high-throughput Illumina Solexa system. Based on sequence similarity and hairpin structure prediction, a total of 29 known miRNA families and 2 novel miRNA families containing a total of 64 miRNA were identified. QRT-PCR analysis revealed that some of the cucumber miRNAs were preferentially expressed in certain tissues. With the recently developed 'high throughput degradome sequencing' approach, 21 target mRNAs of known miRNAs were identified for the first time in cucumber. These targets were associated with development, reactive oxygen species scavenging, signaling transduction and transcriptional regulation. Our study provides an overview of miRNA expression profile and interaction between miRNA and target, which will help further understanding of the important roles of miRNAs in cucumber plants.

Identification of microRNAs from Amur grape (Vitis amurensis Rupr.) by deep sequencing and analysis of microRNA variations with bioinformatics

BACKGROUND

MicroRNA (miRNA) is a class of functional non-coding small RNA with 19-25 nucleotides in length while Amur grape (Vitis amurensis Rupr.) is an important wild fruit crop with the strongest cold resistance among the Vitis species, is used as an excellent breeding parent for grapevine, and has elicited growing interest in wine production. To date, there is a relatively large number of grapevine miRNAs (vv-miRNAs) from cultivated grapevine varieties such as Vitis vinifera L. and hybrids of V. vinifera and V. labrusca, but there is no report on miRNAs from Vitis amurensis Rupr, a wild grapevine species.

RESULTS

A small RNA library from Amur grape was constructed and Solexa technology used to perform deep sequencing of the library followed by subsequent bioinformatics analysis to identify new miRNAs. In total, 126 conserved miRNAs belonging to 27 miRNA families were identified, and 34 known but non-conserved miRNAs were also found. Significantly, 72 new potential Amur grape-specific miRNAs were discovered. The sequences of these new potential va-miRNAs were further validated through miR-RACE, and accumulation of 18 new va-miRNAs in seven tissues of grapevines confirmed by real time RT-PCR (qRT-PCR) analysis. The expression levels of va-miRNAs in flowers and berries were found to be basically consistent in identity to those from deep sequenced sRNAs libraries of combined corresponding tissues. We also describe the conservation and variation of va-miRNAs using miR-SNPs and miR-LDs during plant evolution based on comparison of orthologous sequences, and further reveal that the number and sites of miR-SNP in diverse miRNA families exhibit distinct divergence. Finally, 346 target genes for the new miRNAs were predicted and they include a number of Amur grape stress tolerance genes and many genes regulating anthocyanin synthesis and sugar metabolism.

CONCLUSIONS

Deep sequencing of short RNAs from Amur grape flowers and berries identified 72 new potential miRNAs and 34 known but non-conserved miRNAs, indicating that specific miRNAs exist in Amur grape. These results show that a number of regulatory miRNAs exist in Amur grape and play an important role in Amur grape growth, development, and response to abiotic or biotic stress.

Identification of microRNAs from Amur grape (Vitis amurensis Rupr.) by deep sequencing and analysis of microRNA variations with bioinformatics

Background

MicroRNA (miRNA) is a class of functional non-coding small RNA with 19-25 nucleotides in length while Amur grape (Vitis amurensis Rupr.) is an important wild fruit crop with the strongest cold resistance among the Vitis species, is used as an excellent breeding parent for grapevine, and has elicited growing interest in wine production. To date, there is a relatively large number of grapevine miRNAs (vv-miRNAs) from cultivated grapevine varieties such as Vitis vinifera L. and hybrids of V. vinifera and V. labrusca, but there is no report on miRNAs from Vitis amurensis Rupr, a wild grapevine species.

Results

A small RNA library from Amur grape was constructed and Solexa technology used to perform deep sequencing of the library followed by subsequent bioinformatics analysis to identify new miRNAs. In total, 126 conserved miRNAs belonging to 27 miRNA families were identified, and 34 known but non-conserved miRNAs were also found. Significantly, 72 new potential Amur grape-specific miRNAs were discovered. The sequences of these new potential va-miRNAs were further validated through miR-RACE, and accumulation of 18 new va-miRNAs in seven tissues of grapevines confirmed by real time RT-PCR (qRT-PCR) analysis. The expression levels of va-miRNAs in flowers and berries were found to be basically consistent in identity to those from deep sequenced sRNAs libraries of combined corresponding tissues. We also describe the conservation and variation of va-miRNAs using miR-SNPs and miR-LDs during plant evolution based on comparison of orthologous sequences, and further reveal that the number and sites of miR-SNP in diverse miRNA families exhibit distinct divergence. Finally, 346 target genes for the new miRNAs were predicted and they include a number of Amur grape stress tolerance genes and many genes regulating anthocyanin synthesis and sugar metabolism.

Conclusions

Deep sequencing of short RNAs from Amur grape flowers and berries identified 72 new potential miRNAs and 34 known but non-conserved miRNAs, indicating that specific miRNAs exist in Amur grape. These results show that a number of regulatory miRNAs exist in Amur grape and play an important role in Amur grape growth, development, and response to abiotic or biotic stress.

PAREsnip: a tool for rapid genome-wide discovery of small RNA/target interactions evidenced through degradome sequencing

Small RNAs (sRNAs) are a class of short (20–25 nt) non-coding RNAs that play important regulatory roles in gene expression. An essential first step in understanding their function is to confidently identify sRNA targets. In plants, several classes of sRNAs such as microRNAs (miRNAs) and trans-acting small interfering RNAs have been shown to bind with near-perfect complementarity to their messenger RNA (mRNA) targets, generally leading to cleavage of the mRNA. Recently, a high-throughput technique known as Parallel Analysis of RNA Ends (PARE) has made it possible to sequence mRNA cleavage products on a large-scale. Computational methods now exist to use these data to find targets of conserved and newly identified miRNAs. Due to speed limitations such methods rely on the user knowing which sRNA sequences are likely to target a transcript. By limiting the search to a tiny subset of sRNAs it is likely that many other sRNA/mRNA interactions will be missed. Here, we describe a new software tool called PAREsnip that allows users to search for potential targets of all sRNAs obtained from high-throughput sequencing experiments. By searching for targets of a complete ‘sRNAome’ we can facilitate large-scale identification of sRNA targets, allowing us to discover regulatory interaction networks.

Genome wide expression profiling of two accession of G. herbaceum L. in response to drought

BACKGROUND

Genome-wide gene expression profiling and detailed physiological investigation were used for understanding the molecular mechanism and physiological response of Gossypium herbaceum, which governs the adaptability of plants in drought conditions. Recently, microarray-based gene expression analysis is commonly used to decipher genes and genetic networks controlling the traits of interest. However, the results of such an analysis are often plagued due to a limited number of genes (probe sets) on microarrays. On the other hand, pyrosequencing of a transcriptome has the potential to detect rare as well as a large number of transcripts in the samples quantitatively. We used Affymetrix microarray as well as Roche's GS-FLX transcriptome sequencing for a comparative analysis of cotton transcriptome in leaf tissues under drought conditions.

RESULTS

Fourteen accessions of Gossypium herbaceum were subjected to mannitol stress for preliminary screening; two accessions, namely Vagad and RAHS-14, were selected as being the most tolerant and most sensitive to osmotic stress, respectively. Affymetrix cotton arrays containing 24,045 probe sets and Roche's GS-FLX transcriptome sequencing of leaf tissue were used to analyze the gene expression profiling of Vagad and RAHS-14 under drought conditions. The analysis of physiological measurements and gene expression profiling showed that Vagad has the inherent ability to sense drought at a much earlier stage and to respond to it in a much more efficient manner than does RAHS-14. Gene Ontology (GO) studies showed that the phenyl propanoid pathway, pigment biosynthesis, polyketide biosynthesis, and other secondary metabolite pathways were enriched in Vagad under control and drought conditions as compared with RAHS-14. Similarly, GO analysis of transcriptome sequencing showed that the GO terms responses to various abiotic stresses were significantly higher in Vagad. Among the classes of transcription factors (TFs) uniquely expressed in both accessions, RAHS-14 showed the expression of ERF and WRKY families. The unique expression of ERFs in response to drought conditions reveals that RAHS-14 responds to drought by inducing senescence. This was further supported by transcriptome analysis which revealed that RAHS-14 responds to drought by inducing many transcripts related to senescence and cell death.

CONCLUSION

The comparative genome-wide gene expression profiling study of two accessions of G.herbaceum under drought stress deciphers the differential patterns of gene expression, including TFs and physiologically relevant processes. Our results indicate that drought tolerance observed in Vagad is not because of a single molecular reason but is rather due to several unique mechanisms which Vagad has developed as an adaptation strategy.

Genome-wide identification of Medicago truncatula microRNAs and their targets reveals their differential regulation by heavy metal

We adopted a deep sequencing approach developed by Solexa (Illumina Inc., San Diego, CA, USA) to investigate global expression and complexity of microRNAs (miRNAs) and their targets from Medicago truncatula. Two small RNA libraries and two degradome libraries were constructed from mercury (Hg)-treated and Hg-free M. truncatula seedlings. For miRNAs, each library generated 18.5-18.6 million short sequences, resulting in 10.2-10.8 million clean reads. At least 52 new miRNA candidates with ≈ 21 nucleotides are perfectly matched to the M. truncatula genome. Statistical analysis on transcript abundance of the new candidate miRNAs revealed that most of them were differentially regulated by the heavy metal mercury Hg(II), with 12 miRNAs being specifically induced by Hg exposure. Additionally, we identified 201 individual miRNAs representing 63 known M. truncatula miRNA families, including 12 new conserved and one non-conserved miRNAs that have not been described before. Finally, 130 targets for 58 known (37 conserved and 21 non-conserved) miRNA families and 37 targets for 18 new M. truncatula-specific candidate miRNA families were identified by high-throughput degradome sequencing approach.

Identification of conserved and novel microRNAs that are responsive to heat stress in Brassica rapa

The species Brassica rapa includes various vegetable crops. Production of these vegetable crops is usually impaired by heat stress. Some microRNAs (miRNAs) in Arabidopsis have been considered to mediate gene silencing in plant response to abiotic stress. However, it remains unknown whether or what miRNAs play a role in heat resistance of B. rapa. To identify genomewide conserved and novel miRNAs that are responsive to heat stress in B. rapa, we defined temperature thresholds of non-heading Chinese cabbage (B. rapa ssp. chinensis) and constructed small RNA libraries from the seedlings that had been exposed to high temperature (46 °C) for 1 h. By deep sequencing and data analysis, we selected a series of conserved and novel miRNAs that responded to heat stress. In total, Chinese cabbage shares at least 35 conserved miRNA families with Arabidopsis thaliana. Among them, five miRNA families were responsive to heat stress. Northern hybridization and real-time PCR showed that the conserved miRNAs bra-miR398a and bra-miR398b were heat-inhibitive and guided heat response of their target gene, BracCSD1; and bra-miR156h and bra-miR156g were heat-induced and its putative target BracSPL2 was down-regulated. According to the criteria of miRNA and miRNA* that form a duplex, 21 novel miRNAs belonging to 19 miRNA families were predicted. Of these, four were identified to be heat-responsive by Northern blotting and/or expression analysis of the putative targets. The two novel miRNAs bra-miR1885b.3 and bra-miR5718 negatively regulated their putative target genes. 5'-Rapid amplification of cDNA ends PCR indicated that three novel miRNAs cleaved the transcripts of their target genes where their precursors may have evolved from. These results broaden our perspective on the important role of miRNA in plant responses to heat.

Employing machine learning for reliable miRNA target identification in plants

BACKGROUND

miRNAs are ~21 nucleotide long small noncoding RNA molecules, formed endogenously in most of the eukaryotes, which mainly control their target genes post transcriptionally by interacting and silencing them. While a lot of tools has been developed for animal miRNA target system, plant miRNA target identification system has witnessed limited development. Most of them have been centered around exact complementarity match. Very few of them considered other factors like multiple target sites and role of flanking regions.

RESULT

In the present work, a Support Vector Regression (SVR) approach has been implemented for plant miRNA target identification, utilizing position specific dinucleotide density variation information around the target sites, to yield highly reliable result. It has been named as p-TAREF (plant-Target Refiner). Performance comparison for p-TAREF was done with other prediction tools for plants with utmost rigor and where p-TAREF was found better performing in several aspects. Further, p-TAREF was run over the experimentally validated miRNA targets from species like Arabidopsis, Medicago, Rice and Tomato, and detected them accurately, suggesting gross usability of p-TAREF for plant species. Using p-TAREF, target identification was done for the complete Rice transcriptome, supported by expression and degradome based data. miR156 was found as an important component of the Rice regulatory system, where control of genes associated with growth and transcription looked predominant. The entire methodology has been implemented in a multi-threaded parallel architecture in Java, to enable fast processing for web-server version as well as standalone version. This also makes it to run even on a simple desktop computer in concurrent mode. It also provides a facility to gather experimental support for predictions made, through on the spot expression data analysis, in its web-server version.

CONCLUSION

A machine learning multivariate feature tool has been implemented in parallel and locally installable form, for plant miRNA target identification. The performance was assessed and compared through comprehensive testing and benchmarking, suggesting a reliable performance and gross usability for transcriptome wide plant miRNA target identification.