ARMOUR – A Rice miRNA: mRNA Interaction Resource

ARMOUR was developed as A Rice miRNA:mRNA interaction resource. This informative and interactive database includes the experimentally validated expression profiles of miRNAs under different developmental and abiotic stress conditions across seven Indian rice cultivars. This comprehensive database covers 689 known and 1664 predicted novel miRNAs and their expression profiles in more than 38 different tissues or conditions along with their predicted/known target transcripts. The understanding of miRNA:mRNA interactome in regulation of functional cellular machinery is supported by the sequence information of the mature and hairpin structures. ARMOUR provides flexibility to users in querying the database using multiple ways like known gene identifiers, gene ontology identifiers, KEGG identifiers and also allows on the fly fold change analysis and sequence search query with inbuilt BLAST algorithm. ARMOUR database provides a cohesive platform for novel and mature miRNAs and their expression in different experimental conditions and allows searching for their interacting mRNA targets, GO annotation and their involvement in various biological pathways. The ARMOUR database includes a provision for adding more experimental data from users, with an aim to develop it as a platform for sharing and comparing experimental data contributed by research groups working on rice.

Introduction

MicroRNAs are a class of non-coding RNAs that are transcribed from different loci in the genome and play a vital role in the regulation of gene expression at the transcriptional and post-transcriptional levels (Carthew and Sontheimer, 2009). The miRNAs mainly act to down-regulate cognate transcripts in a sequence specific manner. A single miRNA may regulate the expression of many target transcripts and thus acts as a probable master-switch in biological processes related to growth, development and response to environment (Djami-Tchatchou et al., 2017; Sharma et al., 2017). In plants, a miRNA gene is transcribed by RNA polymerase II into a primary miRNA transcript which is processed by DCL1 (DICER like 1) protein into a hairpin intermediate, called precursor miRNA (pre-miRNA). This is further acted upon by DCL1 to generate 20-24 nt long miRNA/miRNA^∗ duplex in the nucleus. The duplex is then transported to cytoplasm where it is incorporated into the RISC (RNA inducing silencing complex) to guide the AGO (argonaute) protein which results in mRNA cleavage or translational repression (Bartel, 2004; Beauclair et al., 2010). The biogenesis and function of miRNA has been described and reviewed in many articles (Yang and Li, 2012; Tripathi et al., 2015; Mittal et al., 2016). Increasingly, the role of miRNAs is turning out to be much more complex than predicted earlier.

The advent of next generation sequencing (NGS) technologies along with computational approaches have revolutionized the identification and prediction of miRNAs and their targets (Motameny et al., 2010). The experimental validations and related studies have detailed the criteria for miRNA classification and this has formed the basis for several computational algorithms to screen NGS data sets (Unamba et al., 2015; Axtell and Meyers, 2018). Thus NGS is actively replacing hybridization-based methods to catalog and quantify miRNAs in a comprehensive and precise manner (Tripathi et al., 2015). The NGS data sets have also been useful in unraveling the transcriptome profiles thereby providing information on miRNA targets. Our early experiments reported the identification of tissue-preferential (Mittal et al., 2013) and stress-induced (Tripathi et al., 2015, 2018; Sharma et al., 2015; Goswami et al., 2017; Khan et al., 2018) expression patterns of rice miRNAs.

As on date around 10 rice miRNA database are available (Supplementary Table S1). The primary resources for rice miRNAs include miRBase¹ (Griffiths-Jones, 2004; Griffiths-Jones et al., 2006, 2008) and Plant microRNA Database (PMRD²) (Zhang et al., 2009). In addition NCBI GEO/SRA has more than 30 experimental series comprising of hundreds of samples profiled for miRNA using high throughput sequencing approach. The miRBase contains 592 precursors and 713 mature rice miRNAs in the latest release (Version 21) (Kozomara and Griffiths-Jones, 2014) while PMRD contains 2773 miRNAs containing both experimentally validated and computationally predicted sequences (Zhang et al., 2010). PMRD also contains the promoter and miRNA target information for rice miRNAs. Studies based on miRNA sequences and that of their validated targets have facilitated identification of specific binding principles that led to the development of specific target prediction algorithms including psRNATarget³, PicTar (Krek et al., 2005), TargetScan (Lewis et al., 2003) and their combinations. These programs, however, are prone to report a large number of false positive targets. A number of related tools are also available for predicting the biological function of miRNAs and their targets genes like rice gene expression database (ROAD) (Cao et al., 2012), experimentally validated target gene expression database (PMTED⁴) (Sun et al., 2013), plant miRNA expression atlas database and web applications (PmiRExAT⁵) (Gurjar et al., 2016), rice expression profile database (RiceXPro⁶) (Sato et al., 2010), plant miRNA target identification tool (Target-align⁷ (Xie and Zhang, 2010), web server for the prediction of plant miRNA targets (TAPIR⁸) (Bonnet et al., 2010). These and similar databases provide the information relating to specific aspects only, for one or more plants. ARMOUR was designed to provide a cohesive database for complete analysis related to miRNAs in different rice varieties, that is not available so far.

The study of miRNA guided regulatory networks can spread a new light on genetic enhancements of stress tolerance in plants. Undoubtedly NGS has proven to be vital in detecting known and novel miRNAs and measuring their expression changes in rice (Fahlgren et al., 2007; Yang et al., 2011). The accumulation of this data has generated the need for a comprehensive and integrated database of miRNA:mRNA expression profile information and target information, implemented with biologist friendly user interface (UI) or user experience (UX). We present ARMOUR database that consolidates extensive datasets of miRNA deep sequencing studies in different rice varieties under various experimental conditions. It’s UX is designed with four different ways to interact with the database to examine the miRNA:mRNA interaction, with readily accessible information on expression levels. It provides a single platform for integrating information relating to genomic location, expression profiles, biological features, gene ontology and KEGG association for rice miRNA and their targets.

Materials and Methods

Database Design

The UI was designed using HTML5 (Hyper Text Markup Language) and CSS (Cascading Style Sheets). All data for ARMOUR was stored in MySQL v5.6 database, which is a Relational Database Management System (RDBMS) and the most preferred choice for biological databases. In ARMOUR the UX is powered by Hypertext Preprocessor (PHP) and jQuery, while SQL queries were optimized for memory efficient data retrieval and specialized scripts were used for hassle free database updating. The application programming interface (API) is unique feature of the database design and schema (Figure 1). ER (Entity Relationship) tool was used to create an ER diagram that places a total of 13 tables representing comprehensive and heterogeneous information related to miRNA and mRNA that are seamlessly connected in ARMOUR. The UI or UX was designed in four different ways to interact with the database to examine the miRNA:mRNA interaction with expression level information readily accessible (Figure 2). ARMOUR database is designed to be both potable (multiple applications can be developed on the database) and scalable (additional data can be updated time to time) without changing the design.

FIGURE 1

FIGURE 1. Overall schema of database development and functionality representation.

FIGURE 2

FIGURE 2. Database design and user interaction.

Database Content

The database includes miRNAs identified from seven varieties of rice under 38 different developmental and stressed conditions (Table 1) representing leaf, root, flag-leaf, and panicle tissues of seven varieties of rice grown under normal, salt stress and heat stress conditions. The cultivars used are traditional Pusa Basmati 1, dwarf Annapurna, high yielding Lalat, dry season Satabdi, heat-susceptible BPT 5206, salt-tolerant Pokkali and drought/heat-tolerant Nagina 22. The known or identified rice miRNAs were retrieved from miRBase Rel 21 and searched in the NGS libraries while the novel putative miRNAs were predicted using miRcat tool in sRNA Workbench from rice genes or transcripts (Stocks et al., 2012). The database includes 689 known and 1664 predicted novel miRNAs (Table 2). All miRNAs may not be necessarily present in each of the libraries, however, their large numbers indicate the abundance and variability of miRNAs.

TABLE 1

TABLE 1. Libraries details.

TABLE 2

TABLE 2. Key database statistics.

Nomenclature and Gene Annotation

Annotation data type was restricted to VARCHAR to ensure simplicity, fast query and retrieval while all images were generated on the fly through server side scripting. The nomenclature for known miRNA and their target transcripts, in ARMOUR, complies to the standards of miRBase and Rice Genome Annotation Project nomenclature (Kawahara et al., 2013). miRNA identifiers include miRBase ID, miRBase Accession and miRNA Family. All the predicted miRNAs are represented as “Novel-n” (where n is an integer ranging from 1 to 1664) designated as miRNA ID. The predicted novel miRNA ID will be updated periodically as they are validated and entered into miRBase. Using miRNA ID, the miRNA and the precursor sequences can be fetched out, information on the target transcripts can be obtained and the expression profiles of the miRNAs can be obtained. The gene or transcript identifiers that are covered include Entrez gene ID, MSU7 transcript ID, RAPDB gene ID, Uniprot ID, and RefSEQ mRNA ID. The gene or transcript annotation includes gene description, RNA type (coding or non-coding), chromosome number, gene start, gene end, orientation, transcript length, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway information (Kanehisa and Goto, 2000).

Expression Data and Fold Change Analysis

To study the accurate expression profiles of miRNA, the read counts were normalized and provided as log2 of reads per million (RPM; fold expression) in the ARMOUR database. Data normalization was done to ensure minimal redundancy levels since miRNA:mRNA relationship is one is to many and vice versa. Fold Change (FC) analysis can be performed by selecting a reference sample or condition followed by comparison against one or more samples or conditions. The formula used for FC is same as used earlier (Campbell et al., 2015).

$$\begin{array}{c}\text{Fold}\text{expression}\text{=}\text{log2}\left(\text{RPM}\right)\\ \text{FC}\text{=}\text{log2}\left(\text{RPM}\right)\text{test}(-)\text{log2}\left(\text{RPM}\right)\text{control}\end{array}$$

miRNA Targets Identification and Scoring

Transcripts targeted by both known and predicted miRNAs were identified using psRNATarget: A Plant Small RNA Target Analysis Server (Dai and Zhao, 2011) using the rice transcript dataset from Rice Annotation Project Version 7. Targets were qualified based on satisfying all of the following criteria (a) Expectation scores between 0 to 3 (Zhang, 2005), (b) Unpaired Energy Score (UPE) ranging from 0.0 to 25.0, (c) High scoring Segment Pair (HSP) size of 15 to 20 nucleotides, (d) Maximum number of transcripts targeted less than or equal to 200 hits, (e) Flanking length around target for accessibility analysis ranging from 17 nucleotides upstream and 13 nucleotides downstream around target size (Kertesz et al., 2007), and (f) central mismatch leading to translational inhibition centered between 9 to 11 nucleotide from the 5′ end of the miRNA (Brodersen et al., 2008). The information pertaining to transcripts targeted by the miRNA includes type of inhibition (cleavage or translation repression) and multiplicity (number of sites in the transcript likely to be targeted by the miRNA). Later the targets were also analyzed using PsRNA Target ver 2017. However, the 2017 version algorithm gave almost five times more target per miRNAs (∼30–40 mRNA per miRNA) in comparison to 2011 version due to the modified algorithm, indicating the possibility of including false positives (comparison results for 10 miRNAs is provided in Supplementary Table S2). So it was taken care to incorporate the targets that were common to both predictions.

The identified target transcripts were further annotated by identifying their GO categories, that include well planned terminologies, which describe the biological process (P), molecular function (F), and cellular components (C) of gene products, obtained from Rice Genome Annotation Project⁹. To find the associated pathways which are being targeted and affected by these miRNAs through their target transcripts, the information on associated KEGG pathway annotations to the selected target transcripts are also mentioned.

Results and Discussion

Expression Data and Analysis

The behavior of the known and predicted miRNAs listed in the database can be inferred from their expression status in the various libraries representing response to different environmental conditions in different rice varieties. miRNA expression analysis can be performed by filtering the results based on the cut off for fold expression using conditional filter, either in one sample or condition or in multiple samples or conditions. Fold expression can be easily calculated and the results are provided as a flexible Fold Change (FC) cut off to identify up and down regulated miRNA. Fold expression and FC analysis results are provided in an interactive table with expression values and annotation (Figure 3). Within the database there is a provision of using the miRNAs (provided as examples) for retrieving and comparing their expression patterns. The results can be downloaded and used for presentation in many ways (Sharma et al., 2015; Goswami et al., 2017; Khan et al., 2018; Tripathi et al., 2018).

FIGURE 3

FIGURE 3. Accessing of DATA: (A) Searching of data using Ids, (B) searching of data using sequence, (C) Expression based search, (D) Advance search.

ARMOUR will not just act as a knowledge reservoir but also an analytic tool for the researchers to analyze and compare their data with the existing data. The users can send a request to add their own experimental data and compare with the available experimental datasets, thereby allowing sharing and comparision of data between research groups working on rice.

Access to the Database

ARMOUR home page provides an interactive view to the user to become familiar with the database. Key aspects of UI include responsive front end (Device independent), floating frames, fluid design for tables and reader friendly color coding. Various aspects like pre-loaded sample identifiers for each query interface, client side validation of queries, simplified BLAST search, keyword filter on the query results, sorting functionality on table headers, interactive matrix table representation in advanced search, summary charts for GO and pathway queries, sequence level highlighting of miRNA structure and floating buttons ensure ease of navigation on all screens. The user can access any of these search options to retrieve all information related to the query. The results are reported in a unique correlation matrix format, which provides all the information related to the query. The matrix is populated with the number of gene hits that can be re-queried into the database.

Sequence Based Search

Local installation of NCBI-BLAST v 2.2.3 (Johnson et al., 2008) is implemented in ARMOUR database. Backend databases available for search comprise of miRNA precursor sequence database and miRNA mature sequence database from miRBase as well as transcript database obtained from MSU rice database release 7 (Ouyang et al., 2007). BLAST variants that are available for users include blastn, tblastn and tblastx with user defined expect threshold (E-Value) cut-off option for searching (Figure 3). The search can be performed directly from the home page by using the miRNA IDs and by other given keys to obtain information on the miRNA, precursor sequences, expression profiles and target transcripts.

Query Builder Interface

ARMOUR has a unique query builder interface for the users that allows the database to be queried using unlimited keywords or in combination with a list of miRNA or transcript identifiers. Keywords can be searched on Gene description, GO and KEGG pathways and results are reported in a unique correlation matrix format. The matrix is populated with the number of gene hits that can be re-queried into the database. The matrix also shows how many genes are unique and being shared within and across multiple key words (Figure 3).

Using an ID of any miRNA the user can analyze all related data to the query in the database. The GO output includes information on the details of GO category to which selected targets belong and the number of the selected targets belonging to the particular GO. To find the associated pathways which are being targeted and affected by these miRNAs through their target transcripts, the information on associated KEGG pathway annotations to the selected target transcripts are also mentioned (Figure 4).

FIGURE 4

FIGURE 4. Overall representation of a sample miRNA with an example of “novel-1640” miRNA. (A) Result of miRNA search for the selected miRNA IDs. (B) Precursor sequence information. (C) Predicted target gene information. (D) Gene ontology of the predicted target. (E) Pathway analysis of the predicted target. (F) Graphical representation of the target’s Gene ontology.

Conclusion and Perspectives

Rice (Oryza sativa) is one of the most important crops in the world, with a well mapped and well characterized genome. There are several publications relating to role of rice miRNAs in plant development and response to stress (among others Sanan-Mishra et al., 2009; Lv et al., 2010; Ding et al., 2011; Atwell et al., 2014; Mittal et al., 2016). miRNAs are small, non-protein-coding RNA molecules that in association with specific protein complexes, regulate the expression of specific gene products in a sequence-specific manner (Wahid et al., 2010). Number of computational algorithms or tools are available that allow fast and confident prediction of miRNAs and their targets (Jones-Rhoades et al., 2006). The miRNAs have been implicated in numerous developmental and disease states, but their function in distinct biological pathways and phenotypes remains largely unknown (Bazzini et al., 2007; Sharma et al., 2017). Hence, a comprehensive integrome resource of rice miRNA and mRNA would fasten the meta-analysis of miRNA mediated gene regulation in rice at various conditions. This entailed the development of, a rice miRNA and mRNA integrated analysis and user-friendly resource, ARMOUR. It serves as a huge reservoir of known and new rice miRNAs and their expression atlas.

The database is unique in its specific design as it integrates miRNA expression data, from different tissues and varieties of rice, with the predicted target information. This enables analyzing phenotypes and biological functions by associating gene level information from conventional canonical pathways and GO. It provides users with multiple interfaces to interact with the backend database through gene expression, sequence based search and custom query builder. Therefore, it provides a useful resource for researchers investigating the miRNAs and their functional impacts in rice or related cereal crops.

Current version of ARMOUR is focused only on sRNA sequencing for miRNA expression analysis. The existing data in the current version of ARMOUR will be further enriched by adding more experimental data. As future extensions we intend to include and use RNA-seq data and the PARE data or degradome analysis from same tissues and rice varieties for identification and expression analysis of targets. This will give a complete picture and an easy way to compare the expression of miRNAs and their corresponding targets. We envision that the users will also share and integrate their experimental data into ARMOUR to formulate it into a collaborative platform to access and share compiled experimentally validated data resource exclusively for studying rice miRNA and its interaction with its target mRNA. ARMOUR will be constantly updated as soon as new data is available.

Availability and Requirements

Project name: ARMOUR – a rice miRNA: mRNA interaction resource

Project home pages: http://armour.icgeb.trieste.it/login and https://www.icgeb.org/armour.html

Operating system(s): Platform independent

Programming languages: HTML, CSS, MySQL, jQuery.SQL

License: Not required. Each new user can register with his/her details using a valid e-mail ID

Any restrictions to use by non-academics: None.

Author Contributions

NS-M: conceived and designed the experiments. AT, KG, and RS: analyzed the data. RS, MV, HG, and NS-M: contributed in designing of database and analysis tools. NS-M, AT, KG, and MV: wrote the paper.

Funding

This work was supported by grants from ICGEB and Department of Biotechnology, Government of India.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. Deepti Mittal, Dr. Mohammed Aslam, and Dr. Neha Sharma for help with library preparation. We are grateful to Ms. Rashmi Renu Sahoo and Mr. Yusuf Khan for assistance with analyzing the sequencing data. We also thank Prof. Mauro Giacca and Prof. S. K. Sopory for their advice and encouragement on developing the database. We acknowledge the help of Mr. Dario Palmisano in hosting the database on website.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00602/full#supplementary-material

TABLE S1 | List of rice miRNA database.

TABLE S2 | Target predictions results for 10 miRNAs as determined using psRNA Target versions 2011 and 2017, respectively.

Abbreviations

ARMOUR, A rice miRNA: mRNA interaction resource; miRNA, microRNA.

Footnotes

1. ^ http://www.mirbase.org/
2. ^ http://bioinformatics.cau.edu.cn/PMRD/
3. ^ http://plantgrn.noble.org/psRNATarget
4. ^ http://pmted.agrinome.org/index.jsp
5. ^ http://pmirexat.nabi.res.in
6. ^ http://ricexpro.dna.affrc.go.jp/index.html
7. ^ http://www.leonxie.com/targetAlign.php
8. ^ http://bioinformatics.psb.ugent.be/webtools/tapir/
9. ^ http://rice.plantbiology.msu.edu/annotation_pseudo_goslim.shtml

References

Atwell, B. J., Wang, H., and Scafaro, A. P. (2014). Could abiotic stress tolerance in wild relatives of rice be used to improve Oryza sativa? Plant Sci. 21, 48–58. doi: 10.1016/j.plantsci.2013.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Axtell, M. J., and Meyers, B. C. (2018). Guidelines for plant miRNA annotation. Plant Cell 30, 272–284. doi: 10.1105/tpc.17.00851

PubMed Abstract | CrossRef Full Text

Bartel, D. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116, 281–297. doi: 10.1016/S0092-8674(04)00045-5

CrossRef Full Text | Google Scholar

Bazzini, A., Hopp, H., Beachy, R., and Asurmendi, S. (2007). Infection and coaccumulation of tobacco mosaic virus proteins alter microRNA levels, correlating with symptom and plant development. Proc. Natl. Acad. Sci. U.S.A. 104, 12157–12162. doi: 10.1073/pnas.0705114104

PubMed Abstract | CrossRef Full Text | Google Scholar

Beauclair, L., Yu, A., and Bouché, N. (2010). microRNA-directed cleavage and translational repression of the copper chaperone for superoxide dismutase mRNA in Arabidopsis. Plant J. 62, 454–462. doi: 10.1111/j.1365-313X.2010.04162.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonnet, E., He, Y., Billiau, K., and Van De Peer, Y. (2010). TAPIR, a web server for the prediction of plant microRNA targets, including target mimics. Bioinformatics 26, 1566–1568. doi: 10.1093/bioinformatics/btq233

PubMed Abstract | CrossRef Full Text | Google Scholar

Brodersen, P., Sakvarelidze-Achard, L., Bruun-Rasmussen, M., Dunoyer, P., Yamamoto, Y. Y., Sieburth, L., et al. (2008). Widespread translational inhibition by plant miRNAs and siRNAs. Science 320, 1185–1190. doi: 10.1126/science.1159151

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, J. D., Liu, G., Luo, L., Xiao, J., Gerrein, J., Juan-Guardela, B., et al. (2015). Assessment of microRNA differential expression and detection in multiplexed small RNA sequencing data. RNA 21, 164–171. doi: 10.1261/rna.046060.114

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, P., Jung, K.-H., Choi, D., Hwang, D., Zhu, J., and Ronald, P. C. (2012). The rice oligonucleotide array database: an atlas of rice gene expression. Rice 5:17. doi: 10.1186/1939-8433-5-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Carthew, R. W., and Sontheimer, E. J. (2009). Origins and mechanisms of miRNAs and siRNAs. Cell 136, 642–655. doi: 10.1016/j.cell.2009.01.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, X., and Zhao, P. X. (2011). psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 39, W155–W159. doi: 10.1093/nar/gkr319

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, Y., Chen, Z., and Zhu, C. (2011). Microarray-based analysis of cadmium-responsive microRNAs in rice (Oryza sativa). J. Exp. Bot. 62, 3563–3573. doi: 10.1093/jxb/err046

PubMed Abstract | CrossRef Full Text | Google Scholar

Djami-Tchatchou, A. T., Sanan-Mishra, N., Ntushelo, K., and Dubery, I. A. (2017). Functional roles of micrornas in agronomically important plants—potential as targets for crop improvement and protection. Front. Plant Sci. 8:378. doi: 10.3389/fpls.2017.00378

PubMed Abstract | CrossRef Full Text | Google Scholar

Fahlgren, N., Howell, M. D., Kasschau, K. D., Chapman, E. J., Sullivan, C. M., Cumbie, J. S., et al. (2007). High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes. PLoS One 2:e219. doi: 10.1371/journal.pone.0000219

PubMed Abstract | CrossRef Full Text | Google Scholar

Goswami, K., Tripathi, A., and Sanan-Mishra, N. (2017). Comparative miRomics of salt-tolerant and salt-sensitive rice. J. Integr. Bioinform. 14:j/jib.2017. doi: 10.1515/jib-2017-0002

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffiths-Jones, S. (2004). The microRNA Registry. Nucleic Acids Res. 32, D109–D111. doi: 10.1093/nar/gkh023

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffiths-Jones, S., Grocock, R. J., Van Dongen, S., Bateman, A., and Enright, A. J. (2006). miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 34, D140–D144. doi: 10.1093/nar/gkj112

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffiths-Jones, S., Saini, H. K., Van Dongen, S., and Enright, A. J. (2008). miRBase: tools for microRNA genomics. Nucleic Acids Res. 36, D154–D158.

Google Scholar

Gurjar, A. K. S., Panwar, A. S., Gupta, R., and Mantri, S. S. (2016). PmiRExAt: plant miRNA expression atlas database and web applications. Database 2016:baw060. doi: 10.1093/database/baw060

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, M., Zaretskaya, I., Raytselis, Y., Merezhuk, Y., Mcginnis, S., and Madden, T. L. (2008). NCBI BLAST: a better web interface. Nucleic Acids Res. 36, W5–W9. doi: 10.1093/nar/gkn201

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones-Rhoades, M. W., Bartel, D. P., and Bartel, B. (2006). MicroRNAS and their regulatory roles in plants. Annu. Rev. Plant Biol. 57, 19–53. doi: 10.1146/annurev.arplant.57.032905.105218

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., and Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27

CrossRef Full Text | Google Scholar

Kawahara, Y., De La Bastide, M., Hamilton, J. P., Kanamori, H., Mccombie, W. R., Ouyang, S., et al. (2013). Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice 6:4. doi: 10.1186/1939-8433-6-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Kertesz, M., Iovino, N., Unnerstall, U., Gaul, U., and Segal, E. (2007). The role of site accessibility in microRNA target recognition. Nat. Genet. 39, 1278–1284. doi: 10.1038/ng2135

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, A., Goswami, K., and Sanan-Mishra, N. (2018). Mirador” on the potential role of miRNAs in synergy of light and heat networks. Ind. J. Plant Physiol. 22, 1–21.

Google Scholar

Kozomara, A., and Griffiths-Jones, S. (2014). miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 42, D68–D73. doi: 10.1093/nar/gkt1181

PubMed Abstract | CrossRef Full Text | Google Scholar

Krek, A., Grun, D., Poy, M. N., Wolf, R., Rosenberg, L., Epstein, E. J., et al. (2005). Combinatorial microRNA target predictions. Nat. Genet. 37, 495–500. doi: 10.1038/ng1536

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, B. P., Shih, I. H., Jones-Rhoades, M. W., Bartel, D. P., and Burge, C. B. (2003). Prediction of mammalian microRNA targets. Cell 115, 787–798. doi: 10.1016/S0092-8674(03)01018-3

CrossRef Full Text | Google Scholar

Lv, D. K., Bai, X., Li, Y., Ding, X. D., Ge, Y., Cai, H., et al. (2010). Profiling of cold-stress-responsive miRNAs in rice by microarrays. Gene 459, 39–47. doi: 10.1016/j.gene.2010.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittal, D., Mukherjee, S. K., Vasudevan, M., and Sanan-Mishra, N. (2013). Identification of tissue-preferential expression patterns of rice miRNAs. J. Cell. Biochem. 114, 2071–2081. doi: 10.1002/jcb.24552

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittal, D., Sharma, N., Sharma, V., Sopory, S. K., and Sanan-Mishra, N. (2016). Role of microRNAs in rice plant under salt-stress. Ann. Appl. Biol. 168, 2–18. doi: 10.1111/aab.12241

CrossRef Full Text | Google Scholar

Motameny, S., Wolters, S., Nurnberg, P., and Schumacher, B. (2010). Next generation sequencing of miRNAs - strategies. Resource methods. Genes 1, 70–84. doi: 10.3390/genes1010070

PubMed Abstract | CrossRef Full Text | Google Scholar

Ouyang, S., Zhu, W., Hamilton, J., Lin, H., Campbell, M., Childs, K., et al. (2007). The TIGR rice genome annotation resource: improvements and new features. Nucleic Acids Res. 35, D883–D887. doi: 10.1093/nar/gkl976

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanan-Mishra, N., Kumar, V., Sopory, S. K., and Mukherjee, S. K. (2009). Cloning and validation of novel miRNA from basmati rice indicates cross talk between abiotic and biotic stresses. Mol. Genet. Genomics 282, 463–474. doi: 10.1007/s00438-009-0478-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, Y., Antonio, B. A., Namiki, N., Takehisa, H., Minami, H., Kamatsuki, K., et al. (2010). RiceXPro: a platform for monitoring gene expression in japonica rice grown under natural field conditions. Nucleic Acids Res. 39, D1141–D1148. doi: 10.1093/nar/gkq1085

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, N., Mittal, D., and Sanan-Mishra, N. (2017). Micro-regulators of hormones and stress. Mech. Plant Horm. Signal. Under Stress 2, 319–351. doi: 10.1002/9781118889022.ch29

CrossRef Full Text | Google Scholar

Sharma, N., Tripathi, A., and Sanan-Mishra, N. (2015). Profiling the expression domains of a rice-specific microRNA under stress. Front. Plant Sci. 6:333. doi: 10.3389/fpls.2015.00333

PubMed Abstract | CrossRef Full Text | Google Scholar

Stocks, M. B., Moxon, S., Mapleson, D., Woolfenden, H. C., Mohorianu, I., Folkes, L., et al. (2012). The UEA sRNA workbench: a suite of tools for analysing and visualizing next generation sequencing microRNA and small RNA datasets. Bioinformatics 28, 2059–2061. doi: 10.1093/bioinformatics/bts311

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X., Dong, B., Yin, L., Zhang, R., Du, W., Liu, D., et al. (2013). PMTED: a plant microRNA target expression database. BMC Bioinformatics 14:174. doi: 10.1186/1471-2105-14-174

PubMed Abstract | CrossRef Full Text | Google Scholar

Tripathi, A., Goswami, K., and Sanan-Mishra, N. (2015). Role of bioinformatics in establishing microRNAs as modulators of abiotic stress responses: the new revolution. Front. Physiol. 6:286. doi: 10.3389/fphys.2015.00286

PubMed Abstract | CrossRef Full Text | Google Scholar

Tripathi, A., Goswami, K., Tiwari, M., Mukherjee, S. K., and Sanan-Mishra, N. (2018). Identification and comparative analysis of novel microRNAs from tomato varieties showing contrasting response to ToLCV infections. Physiol. Mol. Biol. Plants 24, 185–202. doi: 10.1007/s12298-017-0482-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Unamba, C. I., Nag, A., and Sharma, R. K. (2015). Next generation sequencing technologies: the doorway to the unexplored genomics of non-model plants. Front. Plant Sci. 6:1074. doi: 10.3389/fpls.2015.01074

PubMed Abstract | CrossRef Full Text | Google Scholar

Wahid, F., Shehzad, A., Khan, T., and Kim, Y. Y. (2010). MicroRNAs: synthesis, mechanism, function, and recent clinical trials. Biochim. Biophys. Acta 1803, 1231–1243. doi: 10.1016/j.bbamcr.2010.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, F., and Zhang, B. (2010). Target-align: a tool for plant microRNA target identification. Bioinformatics 26, 3002–3003. doi: 10.1093/bioinformatics/btq568

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., and Li, L. (2012). Analyzing the microRNA transcriptome in plants using deep sequencing data. Biology 1, 297–310. doi: 10.3390/biology1020297

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Zhang, H., and Li, L. (2011). Global analysis of gene-level microRNA expression in Arabidopsis using deep sequencing data. Genomics 98, 40–46. doi: 10.1016/j.ygeno.2011.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y. (2005). miRU: an automated plant miRNA target prediction server. Nucleic Acids Res. 33, W701–W704. doi: 10.1093/nar/gki383

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Yu, J., Li, D., Liu, F., Zhou, X., Wang, T., et al. (2010). PMRD: plant microRNA database. Nucleic Acids Res. 38, D806–D813. doi: 10.1093/nar/gkp818

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Yu, J., Li, D., Zhang, Z., Liu, F., Zhou, X., et al. (2009). PMRD: plant microRNA database. Nucleic Acids Res. 38, D806–D813. doi: 10.1093/nar/gkp818

PubMed Abstract | CrossRef Full Text | Google Scholar