Computational Identification of Novel MicroRNAs and Their Targets in Vigna unguiculata

Differential expression and global analysis of miR156/SQUAMOSA promoter binding-like proteins (SPL) module in oat

SQUAMOSA promoter binding-like proteins (SPLs) are important transcription factors that influence growth phase transition and reproduction in plants. SPLs are targeted by miR156 but the SPL/miR156 module is completely unknown in oat. We identified 28 oat SPL genes (AsSPLs) distributed across all 21 oat chromosomes except for 4C and 6D. The oat- SPL gene family represented six of eight SPL phylogenetic groups, with no AsSPLs in groups 3 and 7. A novel oat miR156 (AsmiR156) family with 21 precursors divided into 7 groups was characterized. A total of 16 AsSPLs were found to be targeted by AsmiR156. Intriguingly, AsSPL3s showed high transcript abundance during early inflorescence (GS-54), as compared to the lower abundance of AsmiR156, indicating their role in reproductive development. Unravelling the SPL/miR156 regulatory hub and alterations in expression patterns of AsSPLs could provide an essential toolbox for genetic improvement in the cultivated oat.

An Integrated Bioinformatics and Functional Approach for miRNA Validation

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

Characterization of microRNA genes from Pigeonpea (Cajanus cajan L.) and understanding their involvement in drought stress

MicroRNAs (miRNA) are non-coding 20-24 nucleotide long RNAs regulating gene expression. In this study, we have characterized and analysed expression of miRNAs in Pigeonpea by using bioinformatics and experimental tools. We identified 116 miRNAs belonging to 32 phylogenetic families. Further, transcription start sites of miRNA genes revealed abundance and unique arrangement of adenine at +1 and thymine at -1 position. Promoter analysis exhibited presence of 19 most prevalent motifs which comprises majorly of TATA box and MYC domains. In total, 252 miRNA-targets were identified and found to be involved in various developmental processes and stress responses. Moreover, genome-wide localization studies demonstrated clustering of cca-miRNA 395 and 169 genes. The tandem and segmental duplication events were observed suggesting miRNA genes have been originated parallelly with protein coding genes. The expression analysis revealed induction of cca-miR169a, 398a and 408 miRNAs under drought stress highlighting their involvement. Conversely, down-regulation of their putative targets (NFYA, SOD, and UCLA, respectively) confirmed regulatory role of miRNAs in their expression emphasising the negative relationship between these miRNAs and targets in Pigeonpea. This study reports vast repertoire of miRNA genes which further can be experimentally characterized to elucidate their functions in various biological processes and can be recommended for Pigeonpea improvement programs.

Identification, characterization, and expression analysis of cowpea (Vigna unguiculata [L.] Walp.) miRNAs in response to cowpea severe mosaic virus (CPSMV) challenge

Cowpea miRNAs and Argonaute genes showed differential expression patterns in response to CPSMV challenge Several biotic stresses affect cowpea production and yield. CPSMV stands out for causing severe negative impacts on cowpea. Plants have two main induced immune systems. In the basal system (PTI, PAMP-triggered immunity), plants recognize and respond to conserved molecular patterns associated with pathogens (PAMPs). The second type (ETI, Effector-triggered immunity) is induced after plant recognition of specific factors from pathogens. RNA silencing is another important defense mechanism in plants. Our research group has been using biochemical and proteomic approaches to learn which proteins and pathways are involved and could explain why some cowpea genotypes are resistant whereas others are susceptible to CPSMV. This current study was conducted to determine the role of cowpea miRNA in the interaction between a resistant cowpea genotype (BRS-Marataoã) and CPSMV. Previously identified and deposited plant microRNA sequences were used to find out all possible microRNAs in the cowpea genome. This search detected 617 mature microRNAs, which were distributed in 89 microRNA families. Next, 4 out of these 617 miRNAs and their possible target genes that encode the proteins Kat-p80, DEAD-Box, GST, and SPB9, all involved in the defense response of cowpea to CPSMV, had their expression compared between cowpea leaves uninoculated and inoculated with CPSMV. Additionally, the differential expression of genes that encode the Argonaute (AGO) proteins 1, 2, 4, 6, and 10 is reported. In summary, the studied miRNAs and AGO 2 and AGO4 associated genes showed differential expression patterns in response to CPSMV challenge, which indicate their role in cowpea defense.

Novel insights into expansion and functional diversification of MIR169 family in tomato

MAIN CONCLUSION

Expansion of MIR169 members by duplication and new mature forms, acquisition of new promoters, differential precursor-miRNA processivity and engaging novel targets increase the functional diversification of MIR169 in tomato. MIR169 family is an evolutionarily conserved miRNA family in plants. A systematic in-depth analysis of MIR169 family in tomato is lacking. We report 18 miR169 precursors, annotating new loci for MIR169a, b and d, as well as 3 novel mature isoforms (MIR169f/g/h). The family has expanded by both tandem- and segmental-duplication events during evolution. A tandem-pair MIR169b/b-1 and MIR169b-2/h is polycistronic in nature coding for three MIR169b isoforms and a new variant miR169h, that is evidently absent in the wild relatives S. pennellii and S. pimpinellifolium. Seven novel miR169 targets including RNA-binding protein, protein-phosphatase, aminotransferase, chaperone, tetratricopeptide-repeat-protein, and transcription factors ARF-9B and SEPELLATA-3 were established by efficient target cleavage in the presence of specific precursors as well as increased target abundance upon miR169 chelation by short-tandem-target-mimic construct in transient assays. Comparative antagonistic expression profiles of MIR169:target pairs suggest MIR169 family as ubiquitous regulator of various abiotic stresses (heat, cold, dehydration and salt) and developmental pathways. This regulation is partly brought about by acquisition of new promoters as demonstrated by promoter MIR169:GUS reporter assays as well as differential processivity of different precursors and miRNA cleavage efficiencies. Thus, the current study augments the functional horizon of MIR169 family with applications for stress tolerance in crops.

In silico determination of transposon-derived miRNAs and targets in Aegilops species

Transposable elements (TEs) are found almost in all living organism, shaping organisms' genomes. miRNAs are noncoding RNA types which are especially important in gene expression regulations. Many previously determined plant miRNAs are identical/homologous to transposons (TE-MIR). The aim of this study was computational characterization of novel TE-related miRNAs and their targets in Aegilops genome by using stringent criteria. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed by BLAST2GO. Seventeen novel TE-related miRNAs in Aegilops genome were identified for the first time. GO analyses indicated that 40 targets played different roles in biological processes, cellular components and molecular functions. Moreover, these genes were involved in 10 metabolic pathways such as purine metabolism, nitrogen metabolism, oxidative phosphorylation, etc. as a result of KEGG analyses. Identification of miRNAs and their targets are significant to understand miRNA-TEs relationships and even how TEs affect plant growth and development. Obtaining results of this study are expected to provide possible new insight into Aegilops and its related species, wheat, with respect to miRNAs evolution and domestication.Communicated by Ramaswamy H. Sarma.

A study of transcriptome in leaf rust infected bread wheat involving seedling resistance gene Lr28

Leaf rust disease causes severe yield losses in wheat throughout the world. During the present study, high-throughput RNA-Seq analysis was used to gain insights into the role of Lr28 gene in imparting seedling leaf rust resistance in wheat. Differential expression analysis was conducted using a pair of near-isogenic lines (NILs) (HD 2329 and HD 2329+Lr28) at early (0h before inoculation (hbi), 24 and 48h after inoculation (hai)) and late stages (72, 96 and 168 hai) after inoculation with a virulent pathotype of pathogen Puccinia triticina. Expression of a large number of genes was found to be affected due to the presence/absence of Lr28. Gene ontology analysis of the differentially expressed transcripts suggested enrichment of transcripts involved in carbohydrate and amino acid metabolism, oxidative stress and hormone metabolism, in resistant and/or susceptible NILs. Genes encoding receptor like kinases (RLKs) (including ATP binding; serine threonine kinases) and other kinases were the most abundant class of genes, whose expression was affected. Genes involved in reactive oxygen species (ROS) homeostasis and several genes encoding transcription factors (TFs) (most abundant being WRKY TFs) were also identified along with some ncRNAs and histone variants. Quantitative real-time PCR was also used for validation of 39 representative selected genes. In the long term, the present study should prove useful in developing leaf rust resistant wheat cultivars through molecular breeding.

sRNAome and transcriptome analysis provide insight into chilling response of cowpea pods

Cowpea is an important horticultural crop in tropical and subtropical areas of Asia, Africa, and Latin America, as well as parts of southern Europe and Central and South America. Chilling injury is a common physiological hazard of cowpea in cold chain logistics which reduce the cowpea pods nutritional quality and product value. However, the molecular mechanism involved in chilling injury remains unclear in cowpea pods. RNA-Seq and sRNA-Seq technologies were employed to decipher the miRNAs and mRNAs expression profiles and their regulatory networks in cowpea pods involved in chilling stress. Differentially expressed miRNAs and mRNA profiles were obtained based on cluster analysis, miRNAs and target genes were found to show coherent relationships in the regulatory networks of chilling injury. Furthermore, we found that numerous miRNAs and nat-siRNAs' targets were predicted to be key enzymes involved in the redox reactions such as POD, CAT, AO and LOX, energy metabolism such as ATPase, FAD and NAD related enzymes and different transcription factors such as WRKY, bHLH, MYB, ERF and NAC which play important roles in chilling injury.

Regulation of Isoflavone Biosynthesis by miRNAs in Two Contrasting Soybean Genotypes at Different Seed Developmental Stages

Owing to the presence of nutritionally important, health-promoting bioactive compounds, especially isoflavones, soybean has acquired the status of a functional food. miRNAs are tiny riboregulator of gene expression by either decreasing and/or increasing the expression of their corresponding target genes. Despite several works on identification and functional characterization of plant miRNAs, the role of miRNAs in the regulation of isoflavones metabolism is still a virgin field. In the present study, we identified a total of 31 new miRNAs along with their 245 putative target genes from soybean seed-specific ESTs using computational approach. The Kyoto Encyclopedia of Genes and Genomes pathway analyses indicated that miRNA putatively regulates metabolism and genetic information processing. Out of that, a total of 5 miRNAs (Gma-miRNA12, Gma-miRNA24, Gma-miRNA26, Gma-miRNA28, and Gma-miRNA29) were predicted and validated for their probable role during isoflavone biosynthesis. We also validated their five target genes using RA-PCR, which is as good as 5'RLM-RACE. Temporal regulation [35 days after flowering, 45, 55, and 65 DAF] of miRNAs and their targets showed differential expression schema. Differential expression of Gma-miR26 and Gma-miRNA28 along with their corresponding target genes (Glyma.10G197900 and Glyma.09G127200) showed a direct relationship with the total isoflavone content. Therefore, understanding the miRNA-based genetic regulation of isoflavone pathway would assist in selection and manipulation to get high-performing soybean genotypes with better isoflavone yield.

Functional Roles of microRNAs in Agronomically Important Plants-Potential as Targets for Crop Improvement and Protection

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have recently emerged as important regulators of gene expression, mainly through cleavage and/or translation inhibition of the target mRNAs during or after transcription. miRNAs play important roles by regulating a multitude of biological processes in plants which include maintenance of genome integrity, development, metabolism, and adaptive responses toward environmental stresses. The increasing population of the world and their food demands requires focused efforts for the improvement of crop plants to ensure sustainable food production. Manipulation of mRNA transcript abundance via miRNA control provides a unique strategy for modulating differential plant gene expression and miRNAs are thus emerging as the next generation targets for genetic engineering for improvement of the agronomic properties of crops. However, a deeper understanding of its potential and the mechanisms involved will facilitate the design of suitable strategies to obtain the desirable traits with minimum trade-offs in the modified crops. In this regard, this review highlights the diverse roles of conserved and newly identified miRNAs in various food and industrial crops and recent advances made in the uses of miRNAs to improve plants of agronomically importance so as to significantly enhance crop yields and increase tolerance to various environmental stress agents of biotic-or abiotic origin.

The Class II Trehalose 6-phosphate Synthase Gene PvTPS9 Modulates Trehalose Metabolism in Phaseolus vulgaris Nodules

Legumes form symbioses with rhizobia, producing nitrogen-fixing nodules on the roots of the plant host. The network of plant signaling pathways affecting carbon metabolism may determine the final number of nodules. The trehalose biosynthetic pathway regulates carbon metabolism and plays a fundamental role in plant growth and development, as well as in plant-microbe interactions. The expression of genes for trehalose synthesis during nodule development suggests that this metabolite may play a role in legume-rhizobia symbiosis. In this work, PvTPS9, which encodes a Class II trehalose-6-phosphate synthase (TPS) of common bean (Phaseolus vulgaris), was silenced by RNA interference in transgenic nodules. The silencing of PvTPS9 in root nodules resulted in a reduction of 85% (± 1%) of its transcript, which correlated with a 30% decrease in trehalose contents of transgenic nodules and in untransformed leaves. Composite transgenic plants with PvTPS9 silenced in the roots showed no changes in nodule number and nitrogen fixation, but a severe reduction in plant biomass and altered transcript profiles of all Class II TPS genes. Our data suggest that PvTPS9 plays a key role in modulating trehalose metabolism in the symbiotic nodule and, therefore, in the whole plant.

Computational Identification, Target Prediction, and Validation of Conserved miRNAs in Insulin Plant (Costus pictus D. Don)

Insulin plant (Costus pictus D. Don) is an economically important medicinal plant for the content of its high value secondary metabolites, bioactive compounds, and remarkable flowering features. MicroRNAs are a class of short (∼21 nucleotides), endogenous, noncoding RNA molecules that play a vital role in regulating gene expression. Here, we used a computer-based homology approach to identify conserved miRNAs in Transcribed Sequence Assemblies (TSA) of C. pictus. It led us to identify 42 miRNAs of 13 different families in C. pictus for the first time. Using quantitative polymerase chain reaction (qPCR) assays, we further confirmed the expression of 8 miRNAs (miR394, miR159b, miR166k, miR172, miR159f, miR166, miR144, and miR858) in young and mature leaf tissues. A total of 109 potential target genes of the identified miRNAs were subsequently predicted in rice (Oryza sativa L.) genome. The target genes encode transcription factors, enzymes, and various functional proteins involved in the regulation of several metabolic pathways. The findings in the present study lay the foundation for further research on miRNAs and miRNA-mediated gene regulation in this important medicinal plant.

An approach to identify the novel miRNA encoded from H. Annuus EST sequences

MicroRNAs are a newly discovered class of non-protein small RNAs with 22-24 nucleotides. They play multiple roles in biological processes including development, cell proliferation, apoptosis, stress responses and many other cell functions. In this research, several approaches were combined to make a computational prediction of potential miRNAs and their targets in Helianthus annuus (H. annuus). The already available information of the plant miRNAs present in miRBase v21 was used against expressed sequence tags (ESTs). A total of three miRNAs were detected from which one potential novel miRNA was identified following a range of strict filtering criteria. The target prediction was carried out for these three miRNAs having various targets. These targets were functionally annotated and GO terms were assigned. To study the conserved nature of the miRNAs, predicted phylogenetic analysis was carried out. These findings will significantly provide the broader picture for understanding the functions in H. annuus.

Computational prediction of miRNAs and their targets in Phaseolus vulgaris using simple sequence repeat signatures

BACKGROUND

MicroRNAs (miRNAs) are endogenous, noncoding, short RNAs directly involved in regulating gene expression at the post-transcriptional level. In spite of immense importance, limited information of P. vulgaris miRNAs and their expression patterns prompted us to identify new miRNAs in P. vulgaris by computational methods. Besides conventional approaches, we have used the simple sequence repeat (SSR) signatures as one of the prediction parameter. Moreover, for all other parameters including normalized Shannon entropy, normalized base pairing index and normalized base-pair distance, instead of taking a fixed cut-off value, we have used 99% probability range derived from the available data.

RESULTS

We have identified 208 mature miRNAs in P. vulgaris belonging to 118 families, of which 201 are novel. 97 of the predicted miRNAs in P. vulgaris were validated with the sequencing data obtained from the small RNA sequencing of P. vulgaris. Randomly selected predicted miRNAs were also validated using qRT-PCR. A total of 1305 target sequences were identified for 130 predicted miRNAs. Using 80% sequence identity cut-off, proteins coded by 563 targets were identified. The computational method developed in this study was also validated by predicting 229 miRNAs of A. thaliana and 462 miRNAs of G. max, of which 213 for A. thaliana and 397 for G. max are existing in miRBase 20.

CONCLUSIONS

There is no universal SSR that is conserved among all precursors of Viridiplantae, but conserved SSR exists within a miRNA family and is used as a signature in our prediction method. Prediction of known miRNAs of A. thaliana and G. max validates the accuracy of our method. Our findings will contribute to the present knowledge of miRNAs and their targets in P. vulgaris. This computational method can be applied to any species of Viridiplantae for the successful prediction of miRNAs and their targets.

Carica papaya microRNAs are responsive to Papaya meleira virus infection

MicroRNAs are implicated in the response to biotic stresses. Papaya meleira virus (PMeV) is the causal agent of sticky disease, a commercially important pathology in papaya for which there are currently no resistant varieties. PMeV has a number of unusual features, such as residence in the laticifers of infected plants, and the response of the papaya to PMeV infection is not well understood. The protein levels of 20S proteasome subunits increase during PMeV infection, suggesting that proteolysis could be an important aspect of the plant defense response mechanism. To date, 10,598 plant microRNAs have been identified in the Plant miRNAs Database, but only two, miR162 and miR403, are from papaya. In this study, known plant microRNA sequences were used to search for potential microRNAs in the papaya genome. A total of 462 microRNAs, representing 72 microRNA families, were identified. The expression of 11 microRNAs, whose targets are involved in 20S and 26S proteasomal degradation and in other stress response pathways, was compared by real-time PCR in healthy and infected papaya leaf tissue. We found that the expression of miRNAs involved in proteasomal degradation increased in response to very low levels of PMeV titre and decreased as the viral titre increased. In contrast, miRNAs implicated in the plant response to biotic stress decreased their expression at very low level of PMeV and increased at high PMeV levels. Corroborating with this results, analysed target genes for this miRNAs had their expression modulated in a dependent manner. This study represents a comprehensive identification of conserved miRNAs inpapaya. The data presented here might help to complement the available molecular and genomic tools for the study of papaya. The differential expression of some miRNAs and identifying their target genes will be helpful for understanding the regulation and interaction of PMeV and papaya.

Identification and characterization of microRNAs in the crab-eating macaque (Macaca fascicularis) using transcriptome analysis

MicroRNAs (miRNAs), with an average length between 16 nt and 26 nt, are small non-coding RNAs that can repress gene expression on the post-transcriptional level. Macaca fascicularis (M. fascicularis), one of the most important nonhuman primate animal models, is widely used in basic and applied preclinical research, especially in studies that involve neuroscience and disease. However, due to the lack of a complete genome sequence, the miRNAs in M. fascicularis have not been completely characterized. In this study, 86 putative M. fascicularis miRNAs were identified using a strategy of our design. The expression of some of these miRNAs in the tissue was confirmed by qRT-PCR. The function and pathway of their targeted genes were analyzed to reveal the potential relevance of miRNA regulation on diseases and physiological processes. The current study provides insight into potential miRNAs and forms a useful knowledge base for the future understanding of the function of miRNAs in M. fascicularis.

Computational Identification of MicroRNAs from the Expressed Sequence Tags of Toxic Dinoflagellate Alexandrium Tamarense

Micro ribonucleic acids (miRNAs) represent a class of small noncoding RNAs that play important roles in multiple biological processes by degrading targeted mRNAs or by repressing mRNA translation. In the case of algal lineages, especially dinoflagellates, knowledge regarding the miRNA system is still limited and its regulatory role remains unclear. In the present study, a computational approach was employed to screen miRNAs from the expressed sequence tags (ESTs) of Alexandrium tamarense. A total of 18 potential miRNAs were identified according to a range of filtering criteria. In addition, unique evolutionary features, such as miRNA gene duplication and sequence similarity to metazoan miRNAs, implied that the miRNA system in dinoflagellates is complex. Moreover, based on these 18 miRNA sequences, 42 potential target genes showing diverse functions in regulating growth and development were predicted in Thalassiosira pseudonana and Phaeodactylum tricornutum. Taken together, our data suggest the existence of miRNAs in dinoflagellates and provide clues for further functional studies on these predicted miRNAs.

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.

MicroRNA prediction and its function in regulating drought-related genes in cowpea

Cowpea has indigenous drought-tolerant characteristics, but the molecular mechanisms underlying the drought-tolerance are largely unknown. Drought sensitive and resistant cowpea have different responses regarding to drought stress. We applied homology search to predict miRNAs and their corresponding targets. The newly identified cowpea miRNAs were validated by real-time quantitative PCR in the leaves and roots of cowpea plants under drought treatment. Target gene prediction shows that a set of miRNA target genes are involved in the metabolic pathways regarding the physiological changes that are highly related to drought stress. We analyzed the expression levels of some important genes that participate in the physiological responses to drought stress and found that variations in their expression levels correspond well to the different responses of drought sensitive and resistant cowpea to drought stress. The expression levels of the target genes were negatively correlated to those of miRNAs. The same miRNA in different tissues responds differently to drought stress. Our results indicate that miRNAs play important roles in response to drought stress by regulating the expression levels of drought-related genes in cowpea.

Unique and conserved microRNAs in wheat chromosome 5D revealed by next-generation sequencing

MicroRNAs are a class of short, non-coding, single-stranded RNAs that act as post-transcriptional regulators in gene expression. miRNA analysis of Triticum aestivum chromosome 5D was performed on 454 GS FLX Titanium sequences of flow-sorted chromosome 5D with a total of 3,208,630 good quality reads representing 1.34x and 1.61x coverage of the short (5DS) and long (5DL) arms of the chromosome respectively. In silico and structural analyses revealed a total of 55 miRNAs; 48 and 42 miRNAs were found to be present on 5DL and 5DS respectively, of which 35 were common to both chromosome arms, while 13 miRNAs were specific to 5DL and 7 miRNAs were specific to 5DS. In total, 14 of the predicted miRNAs were identified in wheat for the first time. Representation (the copy number of each miRNA) was also found to be higher in 5DL (1,949) compared to 5DS (1,191). Targets were predicted for each miRNA, while expression analysis gave evidence of expression for 6 out of 55 miRNAs. Occurrences of the same miRNAs were also found in Brachypodium distachyon and Oryza sativa genome sequences to identify syntenic miRNA coding sequences. Based on this analysis, two other miRNAs: miR1133 and miR167 were detected in B. distachyon syntenic region of wheat 5DS. Five of the predicted miRNA coding regions (miR6220, miR5070, miR169, miR5085, miR2118) were experimentally verified to be located to the 5D chromosome and three of them : miR2118, miR169 and miR5085, were shown to be 5D specific. Furthermore miR2118 was shown to be expressed in Chinese Spring adult leaves. miRNA genes identified in this study will expand our understanding of gene regulation in bread wheat.

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.

Identification of miRNA encoded by Jatropha curcas from EST and GSS

miRNAs are endogenous approx 22 nucleotide RNA which mediates transcriptional or Post-transcriptional gene regulation and play a critical role in diverse aspects of plant development. miRNA identification in wet lab have various constraints, it is time consuming and expensive. It also faces the limitation of identifying miRNAs expressed at specific time and/or at special conditions. Due to the nature of strong conservation of miRNA in plant species, the use of comparative genomics approach for expressed sequence tags (ESTs), Genome Survey Sequence (GSS) and structural feature criteria filter has paved the way toward the identification of conserved miRNAs from the plant species whose genomes are not yet available in public domain. To identify the novel miRNA from Jatropha curcas, a total of 46862 EST sequences and 1569 GSS were searched for homology to previously known viridiplantae 2502 mature miRNA. After predicting the RNA secondary structure, 24 new potential miRNA were identified in J. curcas. Using the newly identified miRNA sequences, a total of 78 potential target genes were identified for 3 miRNA families. Most of the miRNA targeted genes were predicted to encode transcription factors that regulate cell growth and development, signaling, and metabolism. These findings considerably broaden the scope of understanding the functions of miRNA in J. curcas.

Identification and Characterization of MicroRNAs in Macaca fascicularis by EST Analysis

MicroRNAs (miRNAs) are small noncoding RNAs which repress gene expression at the posttranscriptional level. In this study, an expressed sequence tag (EST)-based combined method was applied for the detection of miRNAs in Macaca fascicularis which is used as a model animal extensively in medical experiments, particularly those involved with neuroscience and disease. Initially, previously known miRNA sequences from metazoans were used to blast with the EST databases of Macaca fascicularis, and then a range of filtering criteria was conducted to remove some pseudo ones. At last a total of 8 novel conserved miRNAs were identified; their functions were further predicted and analyzed. Together, our study provides insight into miRNAs and their functions in Macaca fascicularis, indicating that the EST analysis is an efficient and affordable alternative approach for identifying novel miRNA candidates.

Computational identification of miRNAs in medicinal plant Senecio vulgaris (Groundsel)

RNAs Interference plays a very important role in gene silencing. In vitro identification of miRNAs is a slow process as it is difficult to isolate them. Nucleotide sequences of miRNAs are highly conserved among the plants and, this form the key feature behind the identification of miRNAs in plant species by homology alignment. In silico identification of miRNAs from EST database is emerging as a novel, faster and reliable approach. Here EST sequences of Senecio vulgaris (Groundsel) were searched against known miRNA sequences by using BLASTN tool. A total of 10 miRNAs were identified from 1956 EST sequences and 115 GSS sequences. The most stable miRNA identified is svu-mir-1. This approach will accelerate advance research in regulation of gene expression in Groundsel by interfering RNAs.

Identification and comparative analysis of drought-associated microRNAs in two cowpea genotypes

BACKGROUND

Cowpea (Vigna unguiculata) is an important crop in arid and semi-arid regions and is a good model for studying drought tolerance. MicroRNAs (miRNAs) are known to play critical roles in plant stress responses, but drought-associated miRNAs have not been identified in cowpea. In addition, it is not understood how miRNAs might contribute to different capacities of drought tolerance in different cowpea genotypes.

RESULTS

We generated deep sequencing small RNA reads from two cowpea genotypes (CB46, drought-sensitive, and IT93K503-1, drought-tolerant) that grew under well-watered and drought stress conditions. We mapped small RNA reads to cowpea genomic sequences and identified 157 miRNA genes that belong to 89 families. Among 44 drought-associated miRNAs, 30 were upregulated in drought condition and 14 were downregulated. Although miRNA expression was in general consistent in two genotypes, we found that nine miRNAs were predominantly or exclusively expressed in one of the two genotypes and that 11 miRNAs were drought-regulated in only one genotype, but not the other.

CONCLUSIONS

These results suggest that miRNAs may play important roles in drought tolerance in cowpea and may be a key factor in determining the level of drought tolerance in different cowpea genotypes.

Identification of novel soybean microRNAs involved in abiotic and biotic stresses

BACKGROUND

Small RNAs (19-24 nt) are key regulators of gene expression that guide both transcriptional and post-transcriptional silencing mechanisms in eukaryotes. Current studies have demonstrated that microRNAs (miRNAs) act in several plant pathways associated with tissue proliferation, differentiation, and development and in response to abiotic and biotic stresses. In order to identify new miRNAs in soybean and to verify those that are possibly water deficit and rust-stress regulated, eight libraries of small RNAs were constructed and submitted to Solexa sequencing.

RESULTS

The libraries were developed from drought-sensitive and tolerant seedlings and rust-susceptible and resistant soybeans with or without stressors. Sequencing the library and subsequent analyses detected 256 miRNAs. From this total, we identified 24 families of novel miRNAs that had not been reported before, six families of conserved miRNAs that exist in other plants species, and 22 families previously reported in soybean. We also observed the presence of several isomiRNAs during our analyses. To validate novel miRNAs, we performed RT-qPCR across the eight different libraries. Among the 11 miRNAs analyzed, all showed different expression profiles during biotic and abiotic stresses to soybean. The majority of miRNAs were up-regulated during water deficit stress in the sensitive plants. However, for the tolerant genotype, most of the miRNAs were down regulated. The pattern of miRNAs expression was also different for the distinct genotypes submitted to the pathogen stress. Most miRNAs were down regulated during the fungus infection in the susceptible genotype; however, in the resistant genotype, most miRNAs did not vary during rust attack. A prediction of the putative targets was carried out for conserved and novel miRNAs families.

CONCLUSIONS

Validation of our results with quantitative RT-qPCR revealed that Solexa sequencing is a powerful tool for miRNA discovery. The identification of differentially expressed plant miRNAs provides molecular evidence for the possible involvement of miRNAs in the process of water deficit- and rust-stress responses.

Computational identification of miRNAs in medicinal plant Senecio vulgaris (Groundsel)

RNAs Interference plays a very important role in gene silencing. In vitro identification of miRNAs is a slow process as it is difficult to isolate them. Nucleotide sequences of miRNAs are highly conserved among the plants and, this form the key feature behind the identification of miRNAs in plant species by homology alignment. In silico identification of miRNAs from EST database is emerging as a novel, faster and reliable approach. Here EST sequences of Senecio vulgaris (Groundsel) were searched against known miRNA sequences by using BLASTN tool. A total of 10 miRNAs were identified from 1956 EST sequences and 115 GSS sequences. The most stable miRNA identified is svu-mir-1. This approach will accelerate advance research in regulation of gene expression in Groundsel by interfering RNAs.