Identification of conserved micro-RNAs and their target transcripts in opium poppy (Papaver somniferum L.)

Promising approaches for simultaneous enhancement of medicinally significant benzylisoquinoline alkaloids in opium poppy

Benzylisoquinoline alkaloids (BIAs) produced in opium poppy have been evidenced to heal patients suffering from various diseases. They, therefore, hold an integral position in the herbal drug industry. Despite the adoption of several approaches for the large-scale production of BIAs, opium poppy remains the only platform in this purpose. The only disadvantage associated with producing BIAs in the plant is their small quantity. Thus, recruiting strategies that boost their levels is deemed necessary. All the methods which have been employed so far are just able to enhance a maximum of two BIAs. Thus, if these methods are utilized, a sizable amount of time and budget must be spent on the synthesis of all BIAs. Hence, the exploitation of strategies which increase the content of all BIAs at the same time is more commercially effective and time-saving, avoiding the laborious step of resolving the biosynthetic pathway of each compound. Exposure to biotic and abiotic elicitors, development of a synthetic auto-tetraploid, overexpression of a WRKY transcription factor, formation of an artificial metabolon, and suppression of a gene in the shikimate pathway and miRNA are strategies that turn opium poppy into a versatile bioreactor for the concurrent and massive production of BIAs. The last three strategies have never been applied for BIA biosynthetic pathways.

Integrated miRNA, target mRNA, and metabolome profiling of Tinospora cordifolia with reference to berberine biosynthesis

MicroRNAs play a central role in gene regulation and emerge as novel targets for secondary metabolites improvement in plants. The crops thus can be improved through knowledge obtained by the study of miRNAs because of their conserved nature in gene regulation. The present study has been carried out on Tinospora cordifolia (T. cordifolia), because of its illimitable application for the treatment of various diseases. This plant has tremendous medicinal properties, yet unexplored at the molecular level, and has not received much recognition in the scientific field. Thus, here computational analysis was performed to identify T. cordifolia miRNAs using EST database. Using these miRNAs, we predicted their targets which were found to be associated with the regulation of diverse gene networks including 433 berberine biosynthesis genes in T. cordifolia. Further, selected miRNAs were validated and their expression was detected in different T. cordifolia tissues followed by expression analysis of their target mRNAs. These data were then compared with the metabolic profile of T. cordifolia with an emphasis on therapeutically important compound berberine. In this study, we did simultaneous miRNA/target gene expression and metabolome analysis which opens a new way for initiating new proposition and prioritization of miRNAs/genes/metabolites for targeted follow‑up metabolic engineering experimentations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03342-9.

Identification of miRNA and their target genes in Cestrum nocturnum L. and Cestrum diurnum L. in stress responses

MicroRNAs (miRNAs) are small, highly conserved non-coding RNA molecules and products of primary miRNAs that regulate the target gene expression. Homology-based approaches were employed to identify miRNAs and their targets in Cestrum nocturnum L. and Cestrum diurnum L. A total of 32 and 12 miRNA candidates were identified in C. nocturnum and C. diurnum. These miRNAs belong to 26 and 10 miRNA families and regulate 1024 and 1007 target genes in C. nocturnum, and C. diurnum, respectively. The functional roles of these miRNAs have not been earlier elucidated in Cestrum. MiR815a, miR849, miR1089 and miR172 have a strong propensity to target genes controlling phytochrome-interacting factor 1 (PIF1), ubiquitin-specific protease 12 (UBP12), leucine-rich repeat (LRR) protein kinase and GAI, RGA, SCR (GRAS) family transcription factor in C. nocturnum. While miR5205a, miR1436 and miR530 regulate PATATIN-like protein 6 (PLP6), PHD finger transcription factor and myb domain protein 48 (MYB48) in C. diurnum. Overall, these miRNAs have regulatory responses in biotic and abiotic stresses in both plant species. Eight putative miRNAs and their target genes were selected for qRT-PCR validation. The validated results suggested the importance of miR815a, miR849, miR5205a, miR1089, miR172, miR1436, and miR530 in exerting control over stress responses in C. nocturnum and C. diurnum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01127-1.

Bm-miR172c-5p Regulates Lignin Biosynthesis and Secondary Xylem Thickness by Altering the Ferulate 5 Hydroxylase Gene in Bacopa monnieri

MicroRNAs (miRNAs) are small non-coding, endogenous RNAs containing 20-24 nucleotides that regulate the expression of target genes involved in various plant processes. A total of 1,429 conserved miRNAs belonging to 95 conserved miRNA families and 12 novel miRNAs were identified from Bacopa monnieri using small RNA sequencing. The Bm-miRNA target transcripts related to the secondary metabolism were further selected for validation. The Bm-miRNA expression in shoot and root tissues was negatively correlated with their target transcripts. The Bm-miRNA cleavage sites were mapped within the coding or untranslated region as depicted by the modified RLM-RACE. In the present study, we validate three miRNA targets, including asparagine synthetase, cycloartenol synthase and ferulate 5 hydroxylase (F5H) and elucidate the regulatory role of Bm-miR172c-5p, which cleaves the F5H gene involved in the lignin biosynthesis. Overexpression (OE) of Bm-miR172c-5p precursor in B. monnieri suppresses F5H gene, leading to reduced lignification and secondary xylem thickness under control and drought stress. By contrast, OE of endogenous target mimics (eTMs) showed enhanced lignification and secondary xylem thickness leading to better physiological response under drought stress. Taken together, we suggest that Bm-miRNA172c-5p might be a key player in maintaining the native phenotype of B. monnieri under control and different environmental conditions.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Phyto-miRNAs-based regulation of metabolites biosynthesis in medicinal plants

Medicinal plants, are known to produce a wide range of plant secondary metabolites (PSMs) applied as insecticides, drugs, dyes and toxins in agriculture, medicine, industry and bio-warfare plus bio-terrorism, respectively. However, production of PSMs is usually in small quantities, so we need to find novel ways to increase both quantity and quality of them. Fortunately, biotechnology suggests several options through which secondary metabolism in plants can be engineered in innovative ways to: 1) over-produce the useful metabolites, 2) down-produce the toxic metabolites, 3) produce the new metabolites. Among the ways, RNA interference (RNAi) technology which involves gene-specific regulation by small non-coding RNAs (sncRNAs) have been recently emerged as a promising tool for plant biotechnologist, not only to decipher the function of plant genes, but also for development of the plants with improved and novel traits through manipulation of both desirable and undesirable genes. Among sncRNAs, miRNAs have been recorded various regulatory roles in plants such as development, signal transduction, response to environmental stresses, metabolism. Certainly, the use of miRNAs in metabolic engineering requires identification of miRNAs involved in metabolites biosynthesis, understanding of the biosynthetic pathways, as well as the identification of key points of the pathways in which the miRNAs have their own effect. Thus, we firstly consider these three issues on metabolic engineering of medicinal plants. Our review shows, application of miRNAs can open a novel perspective to metabolic engineering of medicinal plants.

Identification and expression analysis of a microRNA cluster derived from pre-ribosomal RNA in Papaver somniferum L. and Papaver bracteatum L

Opium poppy (Papaver somniferum L.) is one of the ancient medical crops, which produces several important alkaloids such as morphine, noscapine, sanguinarine and codeine. MicroRNAs are endogenous non-coding RNAs that play important regulatory roles in plant diverse biological processes. Many plant miRNAs are encoded as single transcriptional units, in contrast to animal miRNAs, which are often clustered. Herein, using computational approaches, a total of 22 miRNA precursors were identified, which five of them were located as a clustered in pre-ribosomal RNA. Afterward, the transcript level of the precursor and the mature of clustered miRNAs in two species of the Papaveraceae family, i.e. P. somniferum L. and P. bracteatum L, were quantified by RT-PCR. With respect to obtained results, these clustered miRNAs were expressed differentially in different tissues of these species. Moreover, using target prediction and Gene Ontology (GO)-based on functional classification indicated that these miRNAs might play crucial roles in various biological processes as well as metabolic pathways. In this study, we discovered the clustered miRNA derived from pre-rRNA, which may shed some light on the importance of miRNAs in the plant kingdom.

Identification of microRNAs and long non-coding RNAs involved in fatty acid biosynthesis in tree peony seeds

MicroRNAs (miRNAs) and long noncoding RNAs (lncRNAs) act as important molecular regulators in a wide range of biological processes during plant development and seed formation, including oil production. Tree peony seeds contain >90% unsaturated fatty acids (UFAs) and high proportions of α-linolenic acid (ALA, > 40%). To dissect the non-coding RNAs (ncRNAs) pathway involved in fatty acids synthesis in tree peony seeds, we construct six small RNA libraries and six transcriptome libraries from developing seeds of two cultivars (J and S) containing different content of fatty acid compositions. After deep sequencing the RNA libraries, the ncRNA expression profiles of tree peony seeds in two cultivars were systematically and comparatively analyzed. A total of 318 known and 153 new miRNAs and 22,430 lncRNAs were identified, among which 106 conserved and 9 novel miRNAs and 2785 lncRNAs were differentially expressed between the two cultivars. In addition, potential target genes of the microRNA and lncRNAs were also predicted and annotated. Among them, 9 miRNAs and 39 lncRNAs were predicted to target lipid related genes. Results showed that all of miR414, miR156b, miR2673b, miR7826, novel-m0027-5p, TR24651|c0_g1, TR24544|c0_g15, and TR27305|c0_g1 were up-regulated and expressed at a higher level in high-ALA cultivar J when compared to low-ALA cultivar S, suggesting that these ncRNAs and target genes are possibly involved in different fatty acid synthesis and lipid metabolism through post-transcriptional regulation. These results provide a better understanding of the roles of ncRNAs during fatty acid biosynthesis and metabolism in tree peony seeds.

Development of a Method to Extract Opium Poppy (Papaver somniferum L.) DNA from Heroin

This study is the first to report the successful development of a method to extract opium poppy (Papaver somniferum L.) DNA from heroin samples. Determining of the source of an unknown heroin sample (forensic geosourcing) is vital to informing domestic and foreign policy related to counter-narcoterrorism. Current profiling methods focus on identifying process-related chemical impurities found in heroin samples. Changes to the geographically distinct processing methods may lead to difficulties in classifying and attributing heroin samples to a region/country. This study focuses on methods to optimize the DNA extraction and amplification of samples with low levels of degraded DNA and inhibiting compounds such as heroin. We compared modified commercial-off-the-shelf extraction methods such as the Qiagen Plant, Stool and the Promega Maxwell-16 RNA-LEV tissue kits for the ability to extract opium poppy DNA from latex, raw and cooked opium, white and brown powder heroin and black tar heroin. Opium poppy DNA was successfully detected in all poppy-derived samples, including heroin. The modified Qiagen stool method with post-extraction purification and a two-stage, dual DNA polymerase amplification procedure resulted in the highest DNA yield and minimized inhibition. This paper describes the initial phase in establishing a DNA-based signature method to characterize heroin.

Turmeric (Curcuma longa): miRNAs and their regulating targets are involved in development and secondary metabolite pathways

Turmeric has been used as a therapeutic herb over centuries in traditional medicinal systems due to the presence of several secondary metabolite compounds. microRNAs are known to regulate gene expression at the post-transcriptional level by transcriptional cleavage or translation repression. miRNAs have been demonstrated to play an active role in secondary metabolism regulation. The present work was focused on the identification of the miRNAs involved in the regulation of secondary metabolite and development process of turmeric. Eighteen miRNA families were identified for turmeric. Sixteen miRNA families were observed to regulate 238 target transcripts. LncRNAs targets of the putative miRNA candidates were also predicted. Our results indicated their role in binding, reproduction, stress, and other developmental processes. Gene annotation and pathway analysis illustrated the biological function of the targets regulated by the putative miRNAs. The miRNA-mediated gene regulatory network also revealed co-regulated targets that were regulated by two or more miRNA families. miR156 and miR5015 were observed to be involved in rhizome development. miR5021 showed regulation for terpenoid backbone biosynthesis and isoquinoline alkaloid biosynthesis pathways. The flavonoid biosynthesis pathway was observed to be regulated by miR2919. The analysis revealed the probable involvement of three miRNAs (miR1168.2, miR156b and miR1858) in curcumin biosynthesis. Other miRNAs were found to be involved in the growth and developmental process of turmeric. Phylogenetic analysis of selective miRNAs was also performed.

Computational Identification of MicroRNAs and Their Transcript Target(s) in Field Mustard (Brassica rapa L.)

BACKGROUND

Micro RNAs (miRNAs) are a pivotal part of non-protein-coding endogenous small RNA molecules that regulate the genes involved in plant growth and development, and respond to biotic and abiotic environmental stresses posttranscriptionally.

OBJECTIVE

In the present study, we report the results of a systemic search for identification of new miRNAs in B. rapa using homology-based ESTs (Expressed Sequence Tags) analysis and considering a series of fi ltration criteria.

MATERIALS AND METHODS

Plant mature miRNA sequences were searched in non-protein coding ESTs registered in NCBI EST database. Zuker RNA folding algorithm was used to generate the secondary structures of the ESTs. Potential sequences were candidate as miRNA genes and characterized evolutionarily only and if only they fi t some described criteria. Also, the web tool psRNATarget was applied to predict candidate B. rapa miRNA targets.

RESULTS

In this study, 10 novel miRNAs from B. rapa belonging to 6 miRNA families were identified using EST-based homology analysis by considering a series of fi ltration criteria. All potent miRNAs appropriate fold back structure. Several potential targets with known/unknown functions for these novel miRNAs were identified. The target genes mainly encode transcription factors, enzymes, DNA binding proteins, disease resistance proteins, carrier proteins and other biological processes.

CONCLUSIONS

MicroRNA having diverse functions in plant species growth, development and evolution by posttranscriptionally regulating the levels of specific transcriptome so by effecting on their translation products. Research in miRNA led to the identification of many miRNAs and their regulating genes from diverse plant species.

Integrated analysis of mRNA and miRNA expression profiling in rice backcrossed progenies (BC2F12) with different plant height

Inter-specific hybridization and backcrossing commonly occur in plants. The use of progeny generated from inter-specific hybridization and backcrossing has been developed as a novel model system to explore gene expression divergence. The present study investigated the analysis of gene expression and miRNA regulation in backcrossed introgression lines constructed from cultivated and wild rice. High-throughput sequencing was used to compare gene and miRNA expression profiles in three progeny lines (L1710, L1817 and L1730), with different plant heights resulting from the backcrossing of introgression lines (BC₂F₁₂) and their parents (O. sativa and O. longistaminata). A total of 25,387 to 26,139 mRNAs and 379 to 419 miRNAs were obtained in these rice lines. More differentially expressed genes and miRNAs were detected in progeny/O. longistaminata comparison groups than in progeny/O. sativa comparison groups. Approximately 80% of the genes and miRNAs showed expression level dominance to O. sativa, indicating that three progeny lines were closer to the recurrent parent, which might be influenced by their parental genome dosage. Approximately 16% to 64% of the differentially expressed miRNAs possessing coherent target genes were predicted, and many of these miRNAs regulated multiple target genes. Most genes were up-regulated in progeny lines compared with their parents, but down-regulated in the higher plant height line in the comparison groups among the three progeny lines. Moreover, certain genes related to cell walls and plant hormones might play crucial roles in the plant height variations of the three progeny lines. Taken together, these results provided valuable information on the molecular mechanisms of hybrid backcrossing and plant height variations based on the gene and miRNA expression levels in the three progeny lines.

Genome-Wide Identification of MicroRNAs and Their Targets in the Leaves and Fruits of Eucommia ulmoides Using High-Throughput Sequencing

MicroRNAs (miRNAs), a group of endogenous small non-coding RNAs, play important roles in plant growth, development, and stress response processes. Eucommia ulmoides Oliver (hardy rubber tree) is one of the few woody plants capable of producing trans-1, 4-polyisoprene (TPI), also known as Eu-rubber, which has been utilized as an industrial raw material and is extensively cultivated in China. However, the mechanism of TPI biosynthesis has not been identified in E. ulmoides. To characterize small RNAs and their targets with potential biological roles involved in the TPI biosynthesis in E. ulmoides, in the present study, eight small RNA libraries were constructed and sequenced from young and mature leaves and fruits of E. ulmoides. Further analysis identified 34 conserved miRNAs belonging to 20 families (two unclassified families), and 115 novel miRNAs seemed to be specific to E. ulmoides. Among these miRNAs, fourteen conserved miRNAs and 49 novel miRNAs were significantly differentially expressed and identified as Eu-rubber accumulation related miRNAs. Based on the E. ulmoides genomic data, 202 and 306 potential target genes were predicted for 33 conserved and 92 novel miRNAs, respectively; the predicted targets are mostly transcription factors and functional genes, which were enriched in metabolic pathways and biosynthesis of secondary metabolites. Noticeably, based on the expression patterns of miRNAs and their target genes in combination with the Eu-rubber accumulation, the negative correlation of expression of six miRNAs (Eu-miR14, Eu-miR91, miR162a, miR166a, miR172c, and miR396a) and their predicted targets serving as potential regulators in Eu-rubber accumulation. This study is the first to detect conserved and novel miRNAs and their potential targets in E. ulmoides and identify several candidate genes potentially controlling rubber accumulation, and thus provide molecular evidence for understanding the roles of miRNAs in regulating the TPI biosynthesis in E. ulmoides.

Functional Characterization of 4'OMT and 7OMT Genes in BIA Biosynthesis

Alkaloids are diverse group of secondary metabolites generally found in plants. Opium poppy (Papaver somniferum L.), the only commercial source of morphinan alkaloids, has been used as a medicinal plant since ancient times. It produces benzylisoquinoline alkaloids (BIA) including the narcotic analgesic morphine, the muscle relaxant papaverine, and the anti-cancer agent noscapine. Though BIAs play crucial roles in many biological mechanisms their steps in biosynthesis and the responsible genes remain to be revealed. In this study, expressions of 3-hydroxy-N-methylcoclaurine 4'-methyltransferase (4'OMT) and reticuline 7-O-methyltransferase (7OMT) genes were subjected to manipulation to functionally characterize their roles in BIA biosynthesis. Measurements of alkaloid accumulation were performed in leaf, stem, and capsule tissues accordingly. Suppression of 4'OMT expression caused reduction in the total alkaloid content in stem tissue whereas total alkaloid content was significantly induced in the capsule. Silencing of the 7OMT gene also caused repression in total alkaloid content in the stem. On the other hand, over-expression of 4'OMT and 7OMT resulted in higher morphine accumulation in the stem but suppressed amount in the capsule. Moreover, differential expression in several BIA synthesis genes (CNMT, TYDC, 6OMT, SAT, COR, 4'OMT, and 7OMT) were observed upon manipulation of 4'OMT and 7OMT expression. Upon silencing and overexpression applications, tissue specific effects of these genes were identified. Manipulation of 4'OMT and 7OMT genes caused differentiated accumulation of BIAs including morphine and noscapine in capsule and stem tissues.

Identification and analysis of miRNAs and their targets in ginger using bioinformatics approach

MicroRNAs (miRNAs) are a large family of endogenous small RNAs derived from the non-protein coding genes. miRNA regulates the gene expression at the post-transcriptional level and plays an important role in plant development. Zingiber officinale is an important medicinal plant having numerous therapeutic properties. Its bioactive compound gingerol and essential oil posses important pharmacological and physiological activities. In this study, we used a homology search based computational approach for identifying miRNAs in Z. officinale. A total of 16 potential miRNA families (miR167, miR407, miR414, miR5015, miR5021, miR5644, miR5645, miR5656, miR5658, miR5664, miR827, miR838, miR847, miR854, miR862 and miR864) were predicted in ginger. Phylogenetic and conserved analyses were performed for predicted miRNAs. Thirteen miRNA families were found to regulate 300 target transcripts and play an important role in cell signaling, reproduction, metabolic process and stress. To understand the miRNA mediated gene regulatory control and to validate miRNA target predictions, a biological network was also constructed. Gene ontology and pathway analyses were also done. miR5015 was observed to regulate the biosynthesis of gingerol by inhibiting phenyl ammonia lyase (PAL), a precursor enzyme in the biosynthesis of gingerol. Our results revealed that most of the predicted miRNAs were involved in the regulation of rhizome development. miR5021, miR854 and miR838 were identified to regulate the rhizome development and the essential oil biosynthesis in ginger.

Transcriptome-wide identification of Rauvolfia serpentina microRNAs and prediction of their potential targets

MicroRNAs (miRNAs) are small non-coding RNAs of ∼ 19-24 nucleotides (nt) in length and considered as potent regulators of gene expression at transcriptional and post-transcriptional levels. Here we report the identification and characterization of 15 conserved miRNAs belonging to 13 families from Rauvolfia serpentina through in silico analysis of available nucleotide dataset. The identified mature R. serpentina miRNAs (rse-miRNAs) ranged between 20 and 22nt in length, and the average minimal folding free energy index (MFEI) value of rse-miRNA precursor sequences was found to be -0.815 kcal/mol. Using the identified rse-miRNAs as query, their potential targets were predicted in R. serpentina and other plant species. Gene Ontology (GO) annotation showed that predicted targets of rse-miRNAs include transcription factors as well as genes involved in diverse biological processes such as primary and secondary metabolism, stress response, disease resistance, growth, and development. Few rse-miRNAs were predicted to target genes of pharmaceutically important secondary metabolic pathways such as alkaloids and anthocyanin biosynthesis. Phylogenetic analysis showed the evolutionary relationship of rse-miRNAs and their precursor sequences to homologous pre-miRNA sequences from other plant species. The findings under present study besides giving first hand information about R. serpentina miRNAs and their targets, also contributes towards the better understanding of miRNA-mediated gene regulatory processes in plants.

Identification of novel and conserved microRNAs in Panax notoginseng roots by high-throughput sequencing

BACKGROUND

MicroRNAs (miRNAs) are small, non-coding RNAs that are important regulators of gene expression, and play major roles in plant development and their response to the environment. Root extracts from Panax notoginseng contain triterpene saponins as their principal bioactive constituent, and demonstrate medicinal properties. To investigate the novel and conserved miRNAs in P. notoginseng, three small RNA libraries constructed from 1-, 2-, and 3-year-old roots in which root saponin levels vary underwent high-throughput sequencing.

METHODS

P. notoginseng roots, purified from 1-, 2-, and 3-year-old roots, were extracted for RNA, respectively. Three small libraries were constructed and subjected to next generation sequencing.

RESULTS

Sequencing of the three libraries generated 67,217,124 clean reads from P. notoginseng roots. A total of 316 conserved miRNAs (belonging to 67 miRNA families and one unclassified family) and 52 novel miRNAs were identified. MIR156 and MIR166 were the largest miRNA families, while miR156i and miR156g showed the highest abundance of miRNA species. Potential miRNA target genes were predicted and annotated using Cluster of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes. Comparing these miRNAs between root samples revealed 33 that were differentially expressed between 2- and 1-year-old roots (8 increased, 25 decreased), 27 differentially expressed between 3- and 1-year-old roots (7 increased, 20 decreased), and 29 differentially expressed between 3- and 2-year-old roots (8 increased, 21 decreased). Two significantly differentially expressed miRNAs and four miRNAs predicted to target genes involved in the terpenoid backbone biosynthesis pathway were selected and validated by quantitative reverse transcription PCR. Furthermore, the expression patterns of these six miRNAs were analyzed in P. notoginseng roots, stems, and leaves at different developmental stages.

CONCLUSIONS

This study identified a large number of P. notoginseng miRNAs and their target genes, functional annotations, and gene expression patterns. It provides the first known miRNA profiles of the P. notoginseng root development cycle.

Bioinformatic identification and expression analysis of banana microRNAs and their targets

MicroRNAs (miRNAs) represent a class of endogenous non-coding small RNAs that play important roles in multiple biological processes by degrading targeted mRNAs or repressing mRNA translation. Thousands of miRNAs have been identified in many plant species, whereas only a limited number of miRNAs have been predicted in M. acuminata (A genome) and M. balbisiana (B genome). Here, previously known plant miRNAs were BLASTed against the Expressed Sequence Tag (EST) and Genomic Survey Sequence (GSS), a database of banana genes. A total of 32 potential miRNAs belonging to 13 miRNAs families were detected using a range of filtering criteria. 244 miRNA:target pairs were subsequently predicted, most of which encode transcription factors or enzymes that participate in the regulation of development, growth, metabolism, and other physiological processes. In order to validate the predicted miRNAs and the mutual relationship between miRNAs and their target genes, qRT-PCR was applied to detect the tissue-specific expression levels of 12 putative miRNAs and 6 target genes in roots, leaves, flowers, and fruits. This study provides some important information about banana pre-miRNAs, mature miRNAs, and miRNA target genes and these findings can be applied to future research of miRNA functions.

Regulation of the alkaloid biosynthesis by miRNA in opium poppy

Opium poppy (Papaver somniferum) is an important medicinal plant producing benzylisoquinoline alkaloids (BIA). MicroRNAs (miRNAs) are endogenous small RNAs (sRNAs) of approximately 21 nucleotides. They are noncoding, but regulate gene expression in eukaryotes. Although many studies have been conducted on the identification and functions of plant miRNA, scarce researches on miRNA regulation of alkaloid biosynthesis have been reported. In this study, a total of 316 conserved and 11 novel miRNAs were identified in opium poppy using second-generation sequencing and direct cloning. Tissue-specific regulation of miRNA expression was comparatively analysed by miRNA microarray assays. A total of 232 miRNAs were found to be differentially expressed among four tissues. Likewise, 1469 target transcripts were detected using in silico and experimental approaches. The Kyoto Encyclopedia of Genes and Genomes pathway analyses indicated that miRNA putatively regulates carbohydrate metabolism and genetic-information processing. Additionally, miRNA target transcripts were mostly involved in response to stress against various factors and secondary-metabolite biosynthesis processes. Target transcript identification analyses revealed that some of the miRNAs might be involved in BIA biosynthesis, such as pso-miR13, pso-miR2161 and pso-miR408. Additionally, three putatively mature miRNA sequences were predicted to be targeting BIA-biosynthesis genes.

Genome-wide fungal stress responsive miRNA expression in wheat

MicroRNAs (miRNAs) are small non-coding class of RNAs. They were identified in many plants with their diverse regulatory roles in several cellular and metabolic processes. A number of miRNAs were involved in biotic and abiotic stress responses. Here, fungal stress responsive wheat miRNAs were analyzed by using miRNA-microarray strategy. Two different fungi (Fusarium culmorum and Bipolaris sorokiniana) were inoculated on resistant and sensitive wheat cultivars. A total of 87 differentially regulated miRNAs were detected in the 8 × 15 K array including all of the available plant miRNAs. Using bioinformatics tools, the target transcripts of responsive miRNAs were predicted, and related biological processes and mechanisms were assessed. A number of the miRNAs such as miR2592s, miR869.1, miR169b were highly differentially regulated showing more than 200-fold change upon fungal-inoculation. Some of the miRNAs were identified as fungal-inoculation responsive for the first time. The analyses showed that some of the differentially regulated miRNAs targeted resistance-related genes such as LRR, glucuronosyl transferase, peroxidase and Pto kinase. The comparison of the two miRNA-microarray analyses indicated that fungal-responsive wheat miRNAs were differentially regulated in pathogen- and cultivar-specific manners.

Deep sequencing identifies tissue-specific microRNAs and their target genes involving in the biosynthesis of tanshinones in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most popular traditional medicinal herbs in Asian nations. Its dried root contains a number of tanshinones, protocatechuic aldehyde, salvianolic acid B and rosmarinic, and is used for the treatment of various diseases. The finding of microRNAs (miRNAs) and their target genes will help understand their biological role on the biosynthesis of tanshinones in S. miltiorrhiza. In the present study, a total of 452 known miRNAs corresponding to 589 precursor miRNAs (pre-miRNAs), and 40 novel miRNAs corresponding to 24 pre-miRNAs were identified in different tissues of S. miltiorrhiza by high-throughput sequencing, respectively. Among them, 62 miRNAs express only in root, 95 miRNAs express only in stem, 19 miRNAs express only in leaf, and 71 miRNAs express only in flower, respectively. By the degradome analysis, 69 targets potentially cleaved by 25 miRNAs were identified. Among them, acetyl-CoA C-acetyltransferase was cleaved by miR5072, and involved in the biosynthesis of tanshinones. This study provided valuable information for understanding the tissue-specific expression patterns of miRNAs in S. miltiorrhiza, and offered a foundation for future studies of the miRNA-mediated biosynthesis of tanshinones.

Identification of conserved and novel microRNAs in Catharanthus roseus by deep sequencing and computational prediction of their potential targets

MicroRNAs are small endogenous non-coding RNAs of ~19-24 nucleotides and perform regulatory roles in many plant processes. To identify miRNAs involved in regulatory networks controlling diverse biological processes including secondary metabolism in Catharanthus roseus, an important medicinal plant, we employed deep sequencing of small RNA from leaf tissue. A total of 88 potential miRNAs comprising of 81 conserved miRNAs belonging to 35 families and seven novel miRNAs were identified. Precursors for 16 conserved and seven novel cro-miRNAs were identified, and their stem-loop hairpin structures were predicted. Selected cro-miRNAs were analyzed by stem-loop qRT-PCR and differential expression patterns were observed in different vegetative tissues of C. roseus. Targets were predicted for conserved and novel cro-miRNAs, which were found to be involved in diverse biological role(s) including secondary metabolism. Our study enriches available resources and information regarding miRNAs and their potential targets for better understanding of miRNA-mediated gene regulation in plants.

Discovery of Novel Leaf Rust Responsive microRNAs in Wheat and Prediction of Their Target Genes

MicroRNAs are endogenous small noncoding RNAs which play critical roles in gene regulation. Few wheat (Triticum aestivum L.) miRNA sequences are available in miRBase repertoire and knowledge of their biological functions related to biotic stress is limited. We identified 52 miRNAs, belonging to 19 families, from next-generation transcriptome sequence data based on homology search. One wheat specific novel miRNA was identified but could not be ascribed or assigned to any known miRNA family. Differentially expressed 22 miRNAs were found between susceptible and resistant wheat near-isogenic lines inoculated with leaf rust pathogen Puccinia triticina and compared with mock inoculated controls. Most miRNAs were more upregulated in susceptible NIL compared to resistant NIL. We identified 1306 potential target genes for these 52 miRNAs with vital roles in response to stimuli, signaling, and diverse metabolic and cellular processes. Gene ontology analysis showed 66, 20, and 35 target genes to be categorized into biological process, molecular function, and cellular component, respectively. A miRNA-mediated regulatory network revealed relationships among the components of the targetome. The present study provides insight into potential miRNAs with probable roles in leaf rust pathogenesis and their target genes in wheat which establish a foundation for future studies.

Differentiated expression of microRNAs may regulate genotype-dependent traits in cotton

miRNA is an exogenous non-coding RNA with 21-24nt in length, which plays a crucial role in almost all biological processes. In plants, miRNAs regulate organ development, phase change, signal transduction and response to different biotic and abiotic stresses at the post-transcriptional levels. Although there are many studies on plant miRNAs, no studies have been focused on the genotype dependence. Genotype-dependent traits may be controlled by the differential expression of certain miRNAs. To test this hypothesis, we investigated the expression profile patterns of 11 selected miRNAs in 5 different organs in 5 different cotton cultivars and their implication on plant development. Our results demonstrate that miRNAs have different expression patterns in different plant organs in different genotypes, which implicate their different traits, including early flowering. miR172 is a miRNA controlling floral development and phase change; our results show that miR172 has a higher expression level in the flower bud than in any other organ, our results also show that Baimian cultivars have a higher expression of miR172 than TM-1. This suggests that Baimian cultivars have an earlier transition from vegetable growth to reproductive growth, which is confirmed by our development data on floral branch development. Our result also shows that several miRNAs, including miR159 and miR162, were highly expressed in Baimian cultivars. The results obtained in this study would provide new insight for improving cotton using miRNA-based biotechnology.

Transcriptome and metabolome analysis of plant sulfate starvation and resupply provides novel information on transcriptional regulation of metabolism associated with sulfur, nitrogen and phosphorus nutritional responses in Arabidopsis

Sulfur is an essential macronutrient for plant growth and development. Reaching a thorough understanding of the molecular basis for changes in plant metabolism depending on the sulfur-nutritional status at the systems level will advance our basic knowledge and help target future crop improvement. Although the transcriptional responses induced by sulfate starvation have been studied in the past, knowledge of the regulation of sulfur metabolism is still fragmentary. This work focuses on the discovery of candidates for regulatory genes such as transcription factors (TFs) using 'omics technologies. For this purpose a short term sulfate-starvation/re-supply approach was used. ATH1 microarray studies and metabolite determinations yielded 21 TFs which responded more than 2-fold at the transcriptional level to sulfate starvation. Categorization by response behaviors under sulfate-starvation/re-supply and other nutrient starvations such as nitrate and phosphate allowed determination of whether the TF genes are specific for or common between distinct mineral nutrient depletions. Extending this co-behavior analysis to the whole transcriptome data set enabled prediction of putative downstream genes. Additionally, combinations of transcriptome and metabolome data allowed identification of relationships between TFs and downstream responses, namely, expression changes in biosynthetic genes and subsequent metabolic responses. Effect chains on glucosinolate and polyamine biosynthesis are discussed in detail. The knowledge gained from this study provides a blueprint for an integrated analysis of transcriptomics and metabolomics and application for the identification of uncharacterized genes.

Candidate loci for phenology and fruitfulness contributing to the phenotypic variability observed in grapevine

In this study, we identified several genes, which potentially contribute to phenological variation in the grapevine. This may help to maintain consistent yield and suitability of particular varieties in future climatic conditions. The timing of major developmental events in fruit crops differs with cultivar, weather conditions and ecological site. This plasticity results also in diverse levels of fruitfulness. Identifying the genetic factors responsible for phenology and fertility variation may help to improve these traits to better match future climates. Two Vitis vinifera populations, an F1 progeny of Syrah × Pinot Noir and a phenological core collection composed of 163 cultivars, were evaluated for phenology and fertility subtraits during three to six growing seasons in the same geographical location. The phenotypic variability in the core collection mostly overlapped with that observed in the F1 progeny and several accessions had exceeding values of phenological response. The progeny population was used together with SSR and SNP markers to map quantitative trait loci (QTLs). This allowed us to detect nine QTLs related to budburst, flowering beginning, the onset of ripening (véraison) and total fertility, explaining from 8 to 44 % of phenotypic variation. A genomic region on chromosome 15 was associated with budburst and véraison and two QTLs for fruitfulness were located on chromosomes 3 and 18. Several genes potentially affecting fertility and the timing of fruit development were proposed, based on their position and putative function. Allelic variation at these candidate loci may be explored by sampling accessions from the core collection.

Computational identification of microRNAs and their targets in cassava (Manihot esculenta Crantz.)

MicroRNAs (miRNAs) are a newly discovered class of noncoding endogenous small RNAs involved in plant growth and development as well as response to environmental stresses. miRNAs have been extensively studied in various plant species, however, only few information are available in cassava, which serves as one of the staple food crops, a biofuel crop, animal feed and industrial raw materials. In this study, the 169 potential cassava miRNAs belonging to 34 miRNA families were identified by computational approach. Interestingly, mes-miR319b was represented as the first putative mirtron demonstrated in cassava. A total of 15 miRNA clusters involving 7 miRNA families, and 12 pairs of sense and antisense strand cassava miRNAs belonging to six different miRNA families were discovered. Prediction of potential miRNA target genes revealed their functions involved in various important plant biological processes. The cis-regulatory elements relevant to drought stress and plant hormone response were identified in the promoter regions of those miRNA genes. The results provided a foundation for further investigation of the functional role of known transcription factors in the regulation of cassava miRNAs. The better understandings of the complexity of miRNA-mediated genes network in cassava would unravel cassava complex biology in storage root development and in coping with environmental stresses, thus providing more insights for future exploitation in cassava improvement.

Identification of microRNAs in Caragana intermedia by high-throughput sequencing and expression analysis of 12 microRNAs and their targets under salt stress

142 miRNAs were identified and 38 miRNA targets were predicted, 4 of which were validated, in C. intermedia . The expression of 12 miRNAs in salt-stressed leaves was assessed by qRT-PCR. MicroRNAs (miRNAs) are endogenous small RNAs that play important roles in various biological and metabolic processes in plants. Caragana intermedia is an important ecological and economic tree species prominent in the desert environment of west and northwest China. To date, no investigation into C. intermedia miRNAs has been reported. In this study, high-throughput sequencing of small RNAs and analysis of transcriptome data were performed to identify both conserved and novel miRNAs, and also their target mRNA genes in C. intermedia. Based on sequence similarity and hairpin structure prediction, 132 putative conserved miRNAs (12 of which were confirmed to form hairpin precursors) belonging to 31 known miRNA families were identified. Ten novel miRNAs (including the miRNA* sequences of three novel miRNAs) were also discovered. Furthermore, 36 potential target genes of 17 known miRNA families and 2 potential target genes of 1 novel miRNA were predicted; 4 of these were validated by 5' RACE. The expression of 12 miRNAs was validated in different tissues, and these and five target mRNAs were assessed by qRT-PCR after salt treatment. The expression levels of seven miRNAs (cin-miR157a, cin-miR159a, cin-miR165a, cin-miR167b, cin-miR172b, cin-miR390a and cin-miR396a) were upregulated, while cin-miR398a expression was downregulated after salt treatment. The targets of cin-miR157a, cin-miR165a, cin-miR172b and cin-miR396a were downregulated and showed an approximately negative correlation with their corresponding miRNAs under salt treatment. These results would help further understanding of miRNA regulation in response to abiotic stress in C. intermedia.

Identification of potential miRNAs and their targets in Vriesea carinata (Poales, Bromeliaceae)

The miRNAs play important roles in regulation of gene expression at the post-transcriptional level. A small RNA and RNA-seq of libraries were constructed to identify miRNAs in Vriesea carinata, a native bromeliad species from Brazilian Atlantic Rainforest. Illumina technology was used to perform high throughput sequencing and data was analyzed using bioinformatics tools. We obtained 2,191,509 mature miRNAs sequences representing 54 conserved families in plant species. Further analysis allowed the prediction of secondary structures for 19 conserved and 16 novel miRNAs. Potential targets were predicted from pre-miRNAs by sequence homology and validated using RTqPCR approach. This study provides the first identification of miRNAs and their potential targets of a bromeliad species.

Insilico profiling of microRNAs in Korean ginseng (Panax ginseng Meyer)

MicroRNAs (miRNAs) are a class of recently discovered non-coding small RNA molecules, on average approximately 21 nucleotides in length, which underlie numerous important biological roles in gene regulation in various organisms. The miRNA database (release 18) has 18,226 miRNAs, which have been deposited from different species. Although miRNAs have been identified and validated in many plant species, no studies have been reported on discovering miRNAs in Panax ginseng Meyer, which is a traditionally known medicinal plant in oriental medicine, also known as Korean ginseng. It has triterpene ginseng saponins called ginsenosides, which are responsible for its various pharmacological activities. Predicting conserved miRNAs by homology-based analysis with available expressed sequence tag (EST) sequences can be powerful, if the species lacks whole genome sequence information. In this study by using the EST based computational approach, 69 conserved miRNAs belonging to 44 miRNA families were identified in Korean ginseng. The digital gene expression patterns of predicted conserved miRNAs were analyzed by deep sequencing using small RNA sequences of flower buds, leaves, and lateral roots. We have found that many of the identified miRNAs showed tissue specific expressions. Using the insilico method, 346 potential targets were identified for the predicted 69 conserved miRNAs by searching the ginseng EST database, and the predicted targets were mainly involved in secondary metabolic processes, responses to biotic and abiotic stress, and transcription regulator activities, as well as a variety of other metabolic processes.

Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng C. A. Meyer

BACKGROUND

Panax ginseng C. A. Meyer is one of the most widely used medicinal plants. Complete genome information for this species remains unavailable due to its large genome size. At present, analysis of expressed sequence tags is still the most powerful tool for large-scale gene discovery. The global expressed sequence tags from P. ginseng tissues, especially those isolated from stems, leaves and flowers, are still limited, hindering in-depth study of P. ginseng.

RESULTS

Two 454 pyrosequencing runs generated a total of 2,423,076 reads from P. ginseng roots, stems, leaves and flowers. The high-quality reads from each of the tissues were independently assembled into separate and shared contigs. In the separately assembled database, 45,849, 6,172, 4,041 and 3,273 unigenes were only found in the roots, stems, leaves and flowers database, respectively. In the jointly assembled database, 178,145 unigenes were observed, including 86,609 contigs and 91,536 singletons. Among the 178,145 unigenes, 105,522 were identified for the first time, of which 65.6% were identified in the stem, leaf or flower cDNA libraries of P. ginseng. After annotation, we discovered 223 unigenes involved in ginsenoside backbone biosynthesis. Additionally, a total of 326 potential cytochrome P450 and 129 potential UDP-glycosyltransferase sequences were predicted based on the annotation results, some of which may encode enzymes responsible for ginsenoside backbone modification. A BLAST search of the obtained high-quality reads identified 14 potential microRNAs in P. ginseng, which were estimated to target 100 protein-coding genes, including transcription factors, transporters and DNA binding proteins, among others. In addition, a total of 13,044 simple sequence repeats were identified from the 178,145 unigenes.

CONCLUSIONS

This study provides global expressed sequence tags for P. ginseng, which will contribute significantly to further genome-wide research and analyses in this species. The novel unigenes identified here enlarge the available P. ginseng gene pool and will facilitate gene discovery. In addition, the identification of microRNAs and the prediction of targets from this study will provide information on gene transcriptional regulation in P. ginseng. Finally, the analysis of simple sequence repeats will provide genetic makers for molecular breeding and genetic applications in this species.

Nutrition metabolism plays an important role in the alternate bearing of the olive tree (Olea europaea L.)

The olive tree (Olea europaea L.) is widely known for its strong tendency for alternate bearing, which severely affects the fruit yield from year to year. Microarray based gene expression analysis using RNA from olive samples (on-off years leaves and ripe-unripe fruits) are particularly useful to understand the molecular mechanisms influencing the periodicity in the olive tree. Thus, we carried out genome wide transcriptome analyses involving different organs and temporal stages of the olive tree using the NimbleGen Array containing 136,628 oligonucleotide probe sets. Cluster analyses of the genes showed that cDNAs originated from different organs could be sorted into separate groups. The nutritional control had a particularly remarkable impact on the alternate bearing of olive, as shown by the differential expression of transcripts under different temporal phases and organs. Additionally, hormonal control and flowering processes also played important roles in this phenomenon. Our analyses provide further insights into the transcript changes between "on year" and "off year" leaves along with the changes from unrpipe to ripe fruits, which shed light on the molecular mechanisms underlying the olive tree alternate bearing. These findings have important implications for the breeding and agriculture of the olive tree and other crops showing periodicity. To our knowledge, this is the first study reporting the development and use of an olive array to document the gene expression profiling associated with the alternate bearing in olive tree.

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.

Genome-wide identification of alternate bearing-associated microRNAs (miRNAs) in olive (Olea europaea L.)

BACKGROUND

Alternate bearing is a widespread phenomenon among crop plants, defined as the tendency of certain fruit trees to produce a high-yield crop one year ("on-year"), followed by a low-yield or even no crop the following year ("off-year"). Several factors may affect the balance between such developmental phase-transition processes. Among them are the microRNA (miRNA), being gene-expression regulators that have been found to be involved as key determinants in several physiological processes.

RESULTS

Six olive (Olea europaea L. cv. Ayvalik variety) small RNA libraries were constructed from fruits (ripe and unripe) and leaves ("on year" and "off year" leaves in July and in November, respectively) and sequenced by high-throughput Illumina sequencing. The RNA was retrotranscribed and sequenced using the high-throughput Illumina platform. Bioinformatics analyses of 93,526,915 reads identified 135 conserved miRNA, belonging to 22 miRNA families in the olive. In addition, 38 putative novel miRNAs were discovered in the datasets. Expression of olive tree miRNAs varied greatly among the six libraries, indicating the contribution of diverse miRNA in balancing between reproductive and vegetative phases. Predicted targets of miRNA were categorized into 108 process ontology groups with significance abundance. Among those, potential alternate bearing-associated processes were found, such as development, hormone-mediated signaling and organ morphogenesis. The KEGG analyses revealed that the miRNA-targeted genes are involved in seven main pathways, belonging to carbohydrate metabolism and hormone signal-transduction pathways.

CONCLUSION

A comprehensive study on olive miRNA related to alternate bearing was performed. Regulation of miRNA under different developmental phases and tissues indicated that control of nutrition and hormone, together with flowering processes had a noteworthy impact on the olive tree alternate bearing. Our results also provide significant data on the miRNA-fruit development interaction and advance perspectives in the miRNA profile of the olive tree.

Genome-wide identification of miRNAs responsive to drought in peach (Prunus persica) by high-throughput deep sequencing

Peach (Prunus persica L.) is one of the most important worldwide fresh fruits. Since fruit growth largely depends on adequate water supply, drought stress is considered as the most important abiotic stress limiting fleshy fruit production and quality in peach. Plant responses to drought stress are regulated both at transcriptional and post-transcriptional level. As post-transcriptional gene regulators, miRNAs (miRNAs) are small (19–25 nucleotides in length), endogenous, non-coding RNAs. Recent studies indicate that miRNAs are involved in plant responses to drought. Therefore, Illumina deep sequencing technology was used for genome-wide identification of miRNAs and their expression profile in response to drought in peach. In this study, four sRNA libraries were constructed from leaf control (LC), leaf stress (LS), root control (RC) and root stress (RS) samples. We identified a total of 531, 471, 535 and 487 known mature miRNAs in LC, LS, RC and RS libraries, respectively. The expression level of 262 (104 up-regulated, 158 down-regulated) of the 453 miRNAs changed significantly in leaf tissue, whereas 368 (221 up-regulated, 147 down-regulated) of the 465 miRNAs had expression levels that changed significantly in root tissue upon drought stress. Additionally, a total of 197, 221, 238 and 265 novel miRNA precursor candidates were identified from LC, LS, RC and RS libraries, respectively. Target transcripts (137 for LC, 133 for LS, 148 for RC and 153 for RS) generated significant Gene Ontology (GO) terms related to DNA binding and catalytic activites. Genome-wide miRNA expression analysis of peach by deep sequencing approach helped to expand our understanding of miRNA function in response to drought stress in peach and Rosaceae. A set of differentially expressed miRNAs could pave the way for developing new strategies to alleviate the adverse effects of drought stress on plant growth and development.

Identification of miR414 and expression analysis of conserved miRNAs from Stevia rebaudiana

MicroRNAs (miRNAs) usually contain 19-24 nucleotides and have been identified as important eukaryotic gene regulators. Applications of various computational approaches have simplified the task by predicting miRNAs from available sequence data sources. In this study, we identified a conserved miR414 from a computational analysis of EST sequence data available from Stevia rebaudiana. In addition, we also identified six conserved miRNAs namely miR169, miR319, miR414, miR164, miR167 and miR398 using stem-loop RT-PCR analysis. Hence, miR414 was commonly identified using both methods. The expression analysis of these miRNAs documented their roles in growth and development of Stevia. Furthermore, the detected miRNAs were found to target genes involved in plant growth, development, metabolism and signal transduction. This is the first study reporting these conserved miRNAs and their expression in Stevia.

Identification and characterization of microRNAs and their target genes in Brassica oleracea

The microRNAs are a new class of small non-coding endogenous RNAs with lengths of approximately ~21 nt. MicroRNAs perform their biological function via the degradation of the target mRNAs or by inhibiting protein translation. Until recently, only limited numbers of miRNAs were identified in Brassica oleracea, a vegetable widely cultivated around the world. In present study, 193 potential miRNA candidates were identified from 17 expressed sequence tag (ESTs) and 152 genome survey sequences (GSSs) in B. oleracea. These miRNA candidates were classified into 70 families using a well-defined comparative genome-based computational analysis. Most miRNAs belong to the miRNA169, miR5021, miR156 and miR158 families. Of these, 36 miRNA families are firstly found in Brassica species. Around 1393 B. oleracea genes were predicted as candidate targets of 175 miRNAs. The mutual relationship between miRNAs and the candidate target genes was verified by checking differentially expression levels using quantitative real-time polymerase chain reaction (qRT-PCR) and 5' RLM-RACE analyses. These target genes participate in multiple biological and metabolic processes, including signal transduction, stress response, and plant development. Gene Ontology analysis shows that the 818, 514, and 265 target genes are involved in molecular functions, biological processes, and cellular component respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis suggests that these miRNAs might regulate 186 metabolic pathways, including those of lipid, energy, starch and sucrose, fatty acid and nitrogen.

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.

Discovering numerical laws of plant microRNA by evolution

MicroRNAs (miRNAs) are small non-coding RNA molecules with important post-transcriptional regulatory roles. They are found in most eukaryotic cell types, including a diverse range of plant species. Previous research has demonstrated that numerical features of miRNA genes, based on sequence and structure, can be indicative of a particular miRNA's specialty. We carried out a large-scale analysis to discover the laws of the relationship between numerical features and evolution. Basic statistical methods were applied to 132 numerical features of 3343 miRNA genes from phylum levels to family levels. Results demonstrated that there were three significant classes among candidate numerical features, the conserved, the increased and the decreased with evolutionary time based on phylum levels. Furthermore, the conserved numerical features among families are usually also conserved among phyla and the features increasing with evolutionary time tended to be related to A or U nucleotides, while the features decreasing with evolutionary time tended to be related to G or C nucleotides. Our results confirm the strong relationship between numerical features of miRNA and species evolution. This research will greatly assist in elucidating the numerical laws of plant miRNA based on evolution.