Abstract
Horsegram (Macrotyloma uniflorum (Lam.) Verdc.) is a drought hardy food and fodder legume of Indo-African continents with diverse germplasm sources demonstrating alternating mechanisms depicting contrasting adaptations to different climatic zones. Tissue specific expression of genes contributes substantially to location specific adaptations. Regulatory networks of such adaptive genes are elucidated for downstream translational research. MicroRNAs are small endogenous regulatory RNAs which alters the gene expression profiles at a particular time and type of tissue. Identification of such small regulatory RNAs in low moisture stress hardy crops can help in cross species transfer and validation confirming stress tolerance ability. This study outlined prediction of conserved miRNAs from transcriptome shotgun assembled sequences and EST sequences of horsegram. We could validate eight out of 15 of the identified miRNAs to demonstrate their role in deficit moisture stress tolerance mechanism of horsegram variety Paiyur1 with their target networks. The putative mumiRs were related to other food legumes indicating the presence of gene regulatory networks. Differential miRNA expression among drought specific tissues indicted the probable energy conservation mechanism. Targets were identified for functional characterization and regulatory network was constructed to find out the probable pathways of post-transcriptional regulation. The functional network revealed mechanism of biotic and abiotic stress tolerance, energy conservation and photoperiod responsiveness.
Similar content being viewed by others
Introduction
MicroRNAs (miRNAs/ miRs) are 19–24 nucleotides long endogenous molecular bigwigs in post-transcriptional gene regulatory networks1. These miR genes are capped, polyadenylated and spliced like other RNA polymerase II transcripts. The mature miRs are located in a hairpin structure within the primary transcript (pri-transcript) and are preprocessed by at least two RNase mediated steps to mature miRNA. Pri-miRNAs range from 50–100 nucleotides that coil into hairpin loop structures encompassing paired stems and unpaired loops2. Intensive investigation over decades has enriched the miRNA’s knowhow in biogenesis and explicit regulatory machinery. Plant miRNA precursors are less conserved, whereas mature miRNAs are conserved at higher magnitude in comparison to animal miRNAs3. Hence, exploring conserved miRNAs makes sense to identify the potential miRNAs from flora with targeted traits.
Macrotyloma uniflorum (Lam.) Verdc. (Horsegram) originated from South-Western India is a multipurpose pulse crop mainly cultivated for providing nutritional security in the form of food, livestock supplement and green manure in India and Africa4,5,6. Being a diploid (2n = 20, 22 or 24], short duration (120–180 days to maturity) plant species and adapted to grow on wide range of agro-climatic conditions, horsegram can be weighed as appropriate model for moisture stress tolerance genes/QTLs investigation challenging undernourishment in drought prone regions6,7. Further, it can be envisaged as nutraceuticals, forage crop8 and anti-calcifying inhibitors9. Recently, substantial improvement in accumulation of EST5,10 and transcriptome data11 has compelled its position as a future crop with versatile utility. The present study details the comprehensive computational approach12,13,14 to predict miRNAs fromavailable ESTs and Transcriptome Shotgun Assembly (TSA) sequences of horsegram based on miRNA homolog search. Potential pathways contributing to drought tolerance is studied as inherent trait of horsegram, which could have common or divergent gene regulatory networks6,10,15.
Flora without whole genome sequences are having alternate resource sequence sources in public databases such as Genome Survey Sequences (GSS), Expressed Sequence Tags (EST) and Bacterial Artificial Chromosome (BAC) sequences rendering plentiful resources to mine conserved microRNAs16,17,18. The database miRBase hosts 8746 reported miRNAs belonging to four phyla of plant species enlisted both mature and precursor miRNAs (https://www.miRbase.org/) (Release 21, June 2014). Plant miRs share functional similarity with small interfering RNAs (siRNAs) in guided target cleavage as microRNA targets sites of coding mRNA sequences3,19. Currently, the data generated from next generation sequencing (NGS) studies have been employed in miRNAs prediction and their impact on multiple traits in related species. The in-silico homology based prediction of miRNA is advantageous and orthologs of previously reported miRNAs could be deciphered with their evolutionary significance among species20,21,22. In recent years, the multi-faceted functionalities of miRs in plants are being understood effectively despite the fact that miRNAs are less studied in plants.
This is the first report in horsegram with comprehensive analysis of conservation and phylogeny of miRs. Previously, a comprehensive study was performed in horsegram from ESTs in predicting eight novel miRs11. Here we report differential expression of identified miRs through a stringent in-silico schema elucidating conserved miRNAs, their characterization, validation with their target annotation and networking.
Materials and methods
Query and reference datasets
Totalling 27,997 shotgun assembled contigs of transcriptome were collected from TSA sequence set: GANR01000001-GANR01027997 from European Nucleotide Archive (www.ebi.ac.uk/ena), SSH-Mu library of moisture stressed cDNA of Macrotyloma (https://www.ncbi.nlm.nih.gov/nucest NCBI dbEST ID 75866463 from) and 1008 ESTs were downloaded (www.ncbi.nlm.nih.gov/dbEST/) from publically available NCBI EST database to represent query sequences for miRNA homolog search. About 8496 miRBase mature Viridiplantae miRNAs were used (https://www.miRbase.org) (Released 21: June, 2014) database23 in the present investigation and clustered by CD- HIT-EST, with threshold value of 10024. Out of 3777 clustered sequences, only non-redundant miRNAs were used as reference miRNAs for finding the homologs in M. uniflorum (Lam.) Verdc. candidate sequences to create a local nucleotide sequence database. To elucidate likely role of miRNAs involved in drought stress tolerance trait, the predicted microRNAs were screened for target genes listed in DroughtdB database (https://pgsb.helmholtz-muenchen.de/droughtdb/).
Bioinformatics tools employed
For conserved miRNA prediction of horsegram candidate sequences from Viridiplantae miRNAs reported in miRBase database, NCBI BLAST version 2.2.2725 an alignment tool was used (https://ftp.ncbi.nih.gov/blast/executables/blast+). The representative sequences of all plant miRNAs were obtained after clustering with CD-HIT-EST with threshold value of 10024. MFOLD online tool26 was (https://unafold.rna.albany.edu/?q=mfold) used for pre-miRNAs secondary structures prediction. The miRNA targets were deduced using psRNAtarget server27. Online circoletto28 tool was used to illustrate the miRNAs and its target genes by circos plot29. The phylogenetic analysis of predicted miRNAs with their closely related miRNA families was performed with MEGA730.
Computational prediction of potential miRNA homologs
The prediction of miRNA homologs was performed using TSA and EST sequences of M. uniflorum available at NCBI and EBI databases. The plant miRNAs were obtained from NCBI was clustered and the representative sequences generated were employed as reference miRNA database using BLAST-2.2.27+. The query sequence consisting of TSA and EST of horsegram were subjected to nucleotide blast (blastn) against reference sequences of plant miRNAs with the following set parameters: (i) word match size-7; (ii) length of mature miRNA sequence ≥ 18 nt without gap; (iii) mismatch range- 0–2 (iv) e-value-0.1. The filtered sequences were utilized to find out the coding and non-coding candidate sequences by performing blastx against non-redundant protein database. The protein coding regions were removed whereas the non-coding regions were exploited for secondary structure prediction and validation.
Secondary structure prediction and validation
With maximum stringency only non-coding sequences were subjected to Zuker’s folding algorithm based secondary structure prediction in MFold26. A sliding window of 100nt from ~ 80nt upstream and ~ 80nt nucleotide downstream of the mature miRNA were set to find the precursors33 with parameters set as: (a) Folding temperature- 37 °C, (b) Ionic conditions of 1 M NaCl without divalent ions, (c) linear RNA sequence, (d) percent sub-optimality number of 5 and (e) maximum interior/bulge loop size-30. As reported elsewhere, to validate the structures of pre-miRs adjusted minimal folding free energy (AMFE) and minimal folding free energy index (MFEI) were calculated34 (Table 1).
The secondary structures were screened for35: (i) Less than three nucleotide substitutions within predicted mature miRNA to reference miR, (ii) The candidate sequence must fold into an appropriate and proper stem-loop hairpin secondary structure, (iii) The localization of mature miRNA must be in one arm of the stem-loop structure, (iv) Less than six mismatches between mature miRNA and its corresponding star sequence (miRNA*) and (v) The secondary structure of putative pre-miR must contain high negative MFE and MFEI levels. Putative horsegram miRNAs were designated based on the standard nomenclature system23,36.
Phylogeny of horsegram miRNA
To enumerate the conserved nature and its phylogenetic relationship of miRs and its precursors were aligned with reported plant miRNAs using BLASTn with e-value, maximum mismatch and hits number set as 10, 3 and 5 respectively. The homologous precursor miRNAs were aligned to predicted miRNAs with Phylogeny.fr web tool (https://www.phylogeny.fr/index.cgi)37 by integrating multiple sequence alignment tools, MUSCLE38 and GBlocks39 to refine the alignment; phylogenetic tree construction by PhyML40 and Tree rendering by TreeDyn41. The alignment output is used in MEGA730 for interpretation of molecular clock to estimate interfamily evolution.
Phylogenetic analysis
The analysis involved 146 nucleotide sequences. All positions containing gaps and missing data were eliminated. Evolutionary analyses were also conducted in MEGA730.
Functional annotation of miRNA targets
As annotations are not available for horsegram, the closely related soybean was used as a reference to annotate miR targets of the predicted miRNAs. Putative identified query miRs were hit against horsegram mRNA sequences using psRNATarget tool27. The parameters set target prediction were (i) 9–11 nucleotide mismatch range for translational inhibition, (ii) maximum 4 mismatches without gap at complementary site, (iii) multiplicity of target sites -2, (iv) maximum expectation value -3 and (v) number of hits -10. The homologs of putative miRNAs and their targets were represented in Glycine max genome as circos plot to demonstrate multiple targets of identified putative miRNAs in horsegram28. In addition, to elucidate the role of identified miRs in drought stress response, specific miRNA target genes among drought tolerance contributing genes reported in Droughtdb (https://pgsb.helmholtz-muenchen.de/droughtdb/) were confirmed. The genes reported for horsegram are insufficient, thus, A. thaliana was used as reference organism to determine gene targets for the new miRNAs predicted. The identified putative miRNAs were used as query against the A. thaliana DFCI gene index (AGI) release 15 and A. thaliana TAIR10, cDNA, removed miRNA gene (release date 14th December 2010) using psRNATarget tool. The parameters set for the target gene prediction were as follows: (i) range of central mismatch for translational inhibition 9–11 nucleotide; (ii) a maximum of 64 mismatches without gaps at complementary site; (iii) multiplicity of target sites 2; (iv) maximum expectation value 3 and (v) Number of hits -5. Similarly, to understand the specific possibilities of putative horsegram miRNAs involved in drought stress, the drought genes were downloaded from Droughtdb in order to search for drought miRNA target genes. The visual representation of miRNA and its targets was represented by drawing circus plot to show the multiple targets of horsegram microRNAs predicted.
Functional annotation and metabolic pathway analysis
Mercator42 AgriGO43 and B2g44 were used to determine the functional roles and to classify Gene Ontology (GO) terms into molecular functions, biological processes and cellular components. Further, corresponding pathways were mined through KASS server, which assigns KEGG Orthology (KO) terms and employs KEGG pathways for miRNA target genes. The simplified schema of workflow adopted in present investigation for miRNA prediction is given in Fig. 1.
Quantitative expression analysis of identified miRNAs
Stress responsive horsegram germplasm were identified from the germplasm core developed as described earlier (6). Selected stress responsive variety Paiyur 1 was raised in glass house under control and moisture stress conditions. Plants were maintained at field condition for control and at temporary wilting point for stress condition. Leaf tissues of control and stress plants were collected and fixed in RNALater and stored at − 80 °C till RNA isolation.
High quality RNA samples were extracted from 100 mg (wet weight) leaf tissue using Plant RNA Easy mini spin coloumn kit following the manufacturer’s guidelines (Qiagen). Quantity and quality of the RNA was quantified with the advanced Qiaexpert.
RNA samples of one microgram each were 3′ polyadenylated using Poly A RNA polymerase (Sigma Aldrich). Forward primers sequences of all 15 miRs and a common Poly T reverse primers were designed (Supplementary material 1).
Quantitative rtPCR using One Step PrimeScript™ RT-PCR Kit with SYBR green (Takara) was performed for all 15 identified miRs. qPCR capturing was conducted using the Illumina Eco RT PCR machine.
For qPCR cycling conditions set were 50 °C for 10 min for cDNA sysnthesis; 95 °C for 3 min for polymerase activation followed by 45° cycles of 95 °C 15 s, 60 °C for 30 s capturing with SYBR green filter at end of each cycle.
Results
miRNAs in M.uniflorum
Significant miRNA homologs within reported 8496 miRNAs were identified by executing nucleotide blast (BLASTn) with 27,997 TSA contigs, SSH-Mu library sequences of moisture stressed horsegram cDNA (NCBI dbEST ID 75866463) and 1008 EST sequences as query resulted in 16 ESTs and 6303 TSA hits (E value 0.1) for further analysis. Stringent filtering (mismatch < 3 and length > 17) was incorporated to narrow down the hits to 5807 (8 ESTs and 5799 TSA) sequences for putative candidate sequences identification with non-coding regions. Based on coding potential, 214 noncoding sequences (7 ESTs and 207 TSA) were filtered out to predict ten distinct pre-miRs coding for 15 conserved mature miRs (Tables 1 and 2) clustering into nine different miR-families (Figs. 2 and 3).
Secondary structure prediction and validation
The predicted miRNAs were analyzed for various structural features to distinguish from other small RNAs such as tRNAs, rRNAs and mRNAs (Table 1 and Fig. 2). Most crucial characteristic feature of stable secondary structure is its minimal folding energy (MFE) ranging from 37.3 to 59.6 (−kcal/mol) for ten in-silico predicted premiRNAs. The MFE of the precursor miRNA was retained low to achieve thermodynamic stability45. Owing to their sequence length polymorphism, pre-miRNAs characterization was based only on MFEI (Minimum free energy index) following previous reports22,46. The MFEI ranges from 0.57 to 1.03 for ten in-silico predicted pre-miRNAs (Fig. 2). The (A + U) % of precursor horsegram miRNAs ranges from 39.24–60.38 satisfies the criteria included by Zhang et al.18. Nucleotide distribution (A = 31.6%, U = 32.31%, G = 26% and C = 23.8%) of the predicted pre-miR are heterogeneous as given Tables 1 and 3. On an average pre-miRs were of 111.4 bp length and mature miRs ranged from 18 to 24 (Fig. 3a) with mean of 20.73 (Table 3). Of predicted miRs, miR482 family had more members (six) representing its presence as one of the supreme molecule in the consortia of horsegram regulatory miRs. Similarly, miR156 and miR390 were found to have two members each whereas, other miRs have singlets from each family (Fig. 3b). Although, no miRs were identified from ESTs and SSH library validating them as coding sequences, the TSA sequences showed miRNA frequency of 1 in 1937 contigs (Table 4). The results of the predicted miRs and its targets in horsegram were statistically analyzed and summarized in Tables 3 and 4.
Evolutionary relationship with soybean
Soybean (G. max L.), the closely related model legume crop with complete genome information cracked has been used to establish a comparative syntenic map of predicted horsegram miRNAs to determine the putative regions of homologous miR genes. From mapping results, except mun-miR-482a-3p and miR2673a, rest of the putative horsegram miRNAs have orthologs widespread in G. max genome. miR156r and miR156k had more number of orthologs. The mapped miRs are depicted in soybean genome (Fig. 4). Comparative genomics enabled us to infer miR gene function, which could enumerate future research focus contributed by horsegram miRNAs in related species and in distant crops as well.
Conservation and phylogenetic analysis
The pre-miR homologs were identified by performing blast of horsegram pre-miRNA against the miRBase database. The hits were filtered to retrieve only miRNAs of same family. In this study, high degree of conservation was exemplified by comparison of mun-miR168 to other plant precursor miRNAs (Figs. 5 and 6). It is evident that horsegram miRNAs share similarity with related legume species (Fig. 7). It is conclusively apparent that, horsegram miRNAs seem to have evolved at different rates in different time period similar to other plants.
Putative target genes, their functions and networks
Inferring the function of miRNA targets is crucial to significantly substantiate the functional role of miRNAs in gene expression and regulation. For predicted 15 horsegram miRs, 39 target genes were identified and classified into different groups based on their functional annotation. As revealed by Mercator annotation results (Fig. 8) of predicted miR targets, diverse processes ranging from RNA transcription regulation, protein posttranscriptional modifications, development, signalling, biotic and abiotic stress tolerance and glycolysis were being regulated by the identified 15 miRs. In addition, target genes also regulate cell wall degradation, hormone biosynthesis and redox components like ascorbate and glutathione synthesis (Table 5). The annotation of target genes using the Mercator tool categorized them into12 broad biochemical processes (Fig. 9). Differential expression of drought responsive miRNAs and gene networks Transcriptome analysis of Illumina sequence data from eight samples representing shoot and root tissues of contrasting horsegram genotypes M-191 (drought sensitive) and M- 249 (drought tolerant) (NCBI Bioproject PRJNA216977) was used to decipher the horsegram miRNA expression in root and shoot samples each under normal and drought stress condition respectively11. The miRNA abundance estimated from blast output indicates the behavioral miRNA expression (Fig. 10). The clustering simplifies the differential gene expression levels of horsegram miRNA and samples based on the similar expression. Network predicted and depicted as target- target interaction (Fig. 11) designates machineries of energy conservation confirming earlier hypothesis of structural compaction and energy conservation to survive under stress conditions14,15. This hypothesis may best fit and appropriate for other crops as well. To extend this hypothesis we identified target homologs in other plants also (Fig. 12).
Discussion
Comprehensive studies on plant miRNAs endorse stage specificity and multitude of targets47,48,49. The non-coding miRNAs play a regulatory role in its target protein coding gene expression50. The miRNA multiplexes with RNA induced silencing complex (RISC) guiding the repression or cleavage of its target messenger RNA by seed nuclei base-pairing3. Complementarity between miRNAs and their target genes are high to regulate developmental processes, metabolism and stress responses51,52. Hence, there is enormous necessity to identify and validate miRNAs for further downstream applications in plants22. In horsegram like deficit moisture stress tolerant crop, major genes interaction network influence in expression of a tolerant phenotype. From contrasting stress responsive genotype data, and qPCR of identified miRNAs from a single stress responsive genotype under irrigated and deficit stress (Fig. 13), we could identify differential expression of miR genes. These predicted 15 miRs were clustered in to 11 different families and are conserved. To validated identified miRNAs qPCR was performed as reported earlier (53). Of the total eight miRs validated, two clusters were formed. Mun-miR 482 was found in both the clusters; whereas, mun-miR1507 was found in one cluster which distinguish the tested variety Paiyur1 for its ability to react for stress (Fig. 14). Of the 15 identified mun-miRs, miR 156, miR 171 and miR 390 were found to be differentially expressed in germinating seeds of halophyte Reaumuria soongorica under salt stress conditions (54). Legume genomes are rich in SNPs. Saturated map of SNPs were reported in legumes like vigna (55) and In pigeonpea (56). Of which highest haplotype desnsity of 0.7380 was reported for serine threonine kinase coding disease resistance gene (56) which was identified as an important miR target in the present ivestigation. Variations in these miR sequences are expected among germplasm to change their expression. This could be a reason for non expression of other genes. The conserved nature of maximum plant miRNAs bolstered the miRNA search in different plant species utilizing available genomic resources like ESTs, GSS and mRNA sequences34,57,58,59 to a great extent and we followed the trend in identifying them in a drought tolerant crop. Conserved miRNAs in ginger, garlic, coffee and tea were identified60,61,62,63. For predicted pre-miRNAs, MFEI resolution for length variation extends from 0.57 to 1.03 in horsegram. Pre-miRs illustrate the absence of large internal loops/bulges and at least 20 nucleotides for Wobble base pairing (G/U base pairings) or Watson–Crick base pairing between the miRNA and the star sequence62,64. Estimated (A + U) % of pre-miRs range satisfies the predetermined criteria18. The putative genomic region of these non-coding sequences were determined and intergenic sequences were utilized for secondary structure prediction with the stringent filtering criteria to attain potential precursor miRNA26. Resultant secondary structures were inspected for precursor miRNA and the positioning of the mature miRNA within its stem-loop structure manually. The sequences with suitable secondary structure66 were characterized for structure attributes and its possible existence as potential precursor of the predicted miRs. The transcription of miRNAs from sense and antisense strands of genes were already reported67. Our results stand by the possibility of the miRNA in both sense and antisense strands (Table 2) similar results were already reported in potato, tobacco and B. rapa16,68,69. The frequency of miRNA using expressed sequence tags is 1 in 1000 confirming previous reports62. The interaction of miRNAs with its target genes can enumerate evolutionary role of microRNAs48,70,71. Utilization of miRNAs in RNA interference (RNAi) mediated gene regulation has emerged as an important tool in novel traits engineering either by over expression of a miRNA or target genes by synthetic miRNAs. Gene silencing/knockout may help in understanding miRNA in plant responses to stress resulting increased productivity with improved nutritional value72. Cellular functions were not known for eight of the target proteins identified in the present investigation. Biotic/Abiotic stresses are critically affecting growth and development of plants. There are reports on contrasting mechanisms confirming grades of deficit moisture stress condition to control among different horsegram germplasm sources6,15. Among the identified miRNAs mun-miR156r, mun-miR156k, mun-miR482a, mun-miR390b-5p and mun-miR482a-5p were noticed to be involved in RNA processing, protein synthesis and modifications as well as plant development by altering the expression of their respective targets, whereas the genes encoding abiotic stress, biotic stress (NBS-LLR resistant class) and signalling associated proteins were targeted by mun-miR482a-5p, mun-miR482d-3p and mun-miR1507a respectively. Thus, there is concurrent substantiation that microRNAs can play pivotal role in crop improvement73 and the miRNAs predicted from this investigation can be useful for future research. Closely related homologs of predicted putative miRs exhibited high degree of conservation as in mum-miR482 with related plant mature miRNAs (Fig. 6). The precursor sequences from miR482, miR390 and miR150 represented a strong candidate promulgating its necessary importance at post-transcriptional gene regulation in horsegram. The conserved nature of precursor and mature miRNAs has been reported in various plant groups from earlier studies16,34,73,74. The potential of miRNAs to bind corresponding target mRNA are comparable with their complementarity to degrade target mRNA. Therefore, to infer the contribution of microRNAs in cellular functions and regulatory gene networks, the miRNA target gene prediction is a crucial step64. The majority of the predicted miRNA targets in horsegram depict energy conservation mechanism which play important role in survival during adverse conditions. Target gene prediction reveals the regulation of multiple genes by single miRNA with different levels of regulation22,47,75. Also, the genes targeted belong to more than one gene family, which shows the multitude of miRNA functions in various metabolic processes (Table 4). We could spot 39 target genes of predicted miRs in horsegram and its homologs in other plant species. This indicates that, drought tolerance phenotype can be manipulated by adjusting the expression of identified miRs and their targets.
Conclusions
Great concern is being publicized recently to unravel various mechanisms involved in drought tolerance with emphasis to changing climatic conditions. Even though, there are extensive inquiries on miRNAs discovery and their functional prediction were completed, some of the non-model plants with considerable traits were not subjected to systematic investigation to elucidate contrasting mechanisms and their key players. Eight novel miRNAs were76 predicted from horsegram. However, these non-validated miRNAs were predicted with low stringency and standard nomenclature was not followed. Inappropriately, the drought hardy horsegram miRNAs still remain unknown and lacks extensive validation. In this study, 15 conserved miRNAs belonging to nine different families were identified from EST and TSA sequences and validated in this first report with its tissue specific differential expression (Fig. 10). Mun-miR482a-5p and mun-miR482d-3p are miR482 family miRNAs reported to regulate stress response in soybean and apple77,78. The same trend is observed in our results confirming the earlier report. Mun-miR390b-5p is a miR390b family miRNA earlier found for protein degradation and post transcriptional modifications in Arabidopsis, maize and soybean79,80,81. Upregulation of miR390b and downregulation of 1507 family miRNA was observed (Fig. 14) in qPCR of Paiyur1 variety which is a stress responsive accession. The present investigation indirectly links stress tolerance to energy conservation as indicated by the network rather than direct response to stress conditions (Fig. 15). Hence, we confirm the interplay of miRNA-stress response-NBS-LLR class R-protein response77 in energy conservation to survive the stress conditions. Additionally, the study reveals that identified miRNA regulated target genes have differential biological functions including cell wall degradation, hormone synthesis and synthesis of redox component like ascorbate and glutathione (Table 5). These findings can further help in translational research of related legumes and may result in selection of better germplasm for higher productivity and nutritional security. The network predicted from this investigation confirms the earlier hypotheses: structural compaction to overcome stress tolerance14,15. Altogether, the outcomes of the present investigation deliver clarity that, horsegram is a drought adapted crop and can be considered as a model crop for drought tolerance research. Thus, the predicted horsegram miRNAs may unravel the unique tolerance capability associated to its metabolic pathways and the present workflow represents simple and straightforward approach for the prediction and characterization of miRNAs in those plants for which genomes are yet to be splintered.
References
Jones-Rhoades, M. W., Bartel, D. P. & Bartel, B. MicroRNAs and their regulatory roles in plants. Annu. Rev. Plant Biol. 57, 19–53 (2006).
Pritchard, C. C., Cheng, H. H. & Tewari, M. MicroRNA profiling: approaches and considerations. Nat. Rev. Genet. 13, 358 (2012).
Bartel, D. P. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116, 281–297 (2004).
Arora, R. K. & Chandel, K. P. S. Botanical source areas of wild herbage legumes in India. Trop. Grasslands. 6, 213–221 (1972).
Yasin, J. K. et al. Identification and validation of genes responsible for moisture stress tolerance in Horsegram (Macrotyloma uniflorum(Lam.) Verdc.), LIBEST_027585.2011. https://www.ncbi.nlm.nih.gov/biosample/SAMN00750282 (2011).
Yasin, J. K. et al. Alternate antioxidant defence system in moisture stress responsive accessions of horse gram. Legume Res. 37, 145–154 (2014).
Cook, B.G., Pengelly, B.C., Brown, S.D., Donnelly, J.L., Eagles, D.A., Franco, M.A., Hanson, J., Mullen, B.F., Partridge, I.J. & Peters, M. Tropical Forages: An interactive selection tool (2005).
Morris, J. B. Macrotyloma axillare and M. uniflorum: descriptor analysis, anthocyanin indexes, and potential uses. Genet. Resour. Crop Evol. 55, 5–8 (2008).
Peshin, A. & Singla, S. K. Anticalcifying properties of Dolichos biflorus(horse gram) Seeds. Indian J. Exp. Biol. 32, 889–891 (1994).
Reddy, P. C. O. et al. Identification of stress-induced genes from the drought tolerant semi-arid legume crop horsegram (Macrotyloma uniflorum (Lam.) Verdc.) through analysis of subtracted expressed sequence tags. Plant Sci. 175, 372–384 (2008).
Bhardwaj, J. et al. Comprehensive transcriptomic study on horse gram (Macrotyloma uniflorum): De novo assembly, functional characterization and comparative analysis in relation to drought stress. BMC Genomics 14, 1 (2013).
Yasin, J.K., Mishra, B. K., Chaudhary, S., Magadum, S., Chinnusamy, V. & Singh, N.K. Transcriptome analyses for genome wide identification of ncRNAs and miRNA from Pigeonpea (Cajanus cajanL.). In “Plant Biology 2016” 13–18th July, 2016 by ASPB, Austin, Texas, USA (2016).
Yasin, J.K., Sreevathsa, R., Vivek, T., Nager, R., Lal, S.K., Pillai, M.A. & Chinnusamy, V. lncRNA plays a major role in Susceptibility to resistance: Insights into pigeonpea RGA. In “Plant Biology 2016” 13–18th July, 2016 by ASPB, Austin, Texas, USA (2016).
Yasin, J.K. & Magadum, S. Structural compaction to conserve energy: ncRNA expression directs pH flux of floral parts and yield loss in pigeonpea (Cajanus cajanL.). In Royal Society Theo murphy meeting on “Evolution brings Ca2+ and ATP together to control life and death. March 16–17th 2016 at Royal Society of UK, London (2016).
Yasin, J.K., Bhat, K.V., Rajkumar, S., Subalakshmi, Ramya, K.T. & Fiyaz, A.R. Structural compaction, Mechanism of acid tolerance in moisture stress responsive accessions of horse gram. In The 8TH International symposium on “Plant soil interactions at low pH”, October 18-22, 2012, Bangalore, India (2012); 170
Dhandapani, V. et al. Identification of potential microRNAs and their targets in Brassica rapa L. Mol. Cells. 32, 21–37 (2011).
Xie, F. L. et al. Computational identification of novel microRNAs and targets in Brassica napus. FEBS Lett. 581, 1464–1474 (2007).
Zhang, B. H., Pan, X., Cobb, G. P. & Anderson, T. A. Evidence that miRNAs are different from other RNAs. Cell Mol. Life Sci. 63, 246–254 (2006).
Yang, T., Xue, L. & An, L. Functional diversity of miRNA in plants. Plant Sci. 172, 423–432 (2007).
Gupta, P. K. MicroRNAs and target mimics for crop improvement. Curr. Sci. 108, 1624–1633 (2015).
Zhang, B. H., Pan, X. P., Wang, Q. L., Cobb, G. P. & Anderson, T. A. Identification and characterization of new plant microRNAs using EST analysis. Cell Res. 15, 336–360 (2005).
Zhang, B. H., Pan, X., Cobb, G. P. & Anderson, T. A. Plant microRNA: a small regulatory molecule with big impact. Dev. Biol. 289, 3–16 (2006).
Jones, G. S., Grocock, R. J., Van Dongen, S., Bateman, A. & Enright, A. J. MiRBase: microRNA sequences, targets and gene nomenclature. Nucl. Acids Res. 34, D140–D144 (2006).
Li, W. & Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659 (2006).
Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410 (1990).
Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucl. Acids Res. 31, 3406–3415 (2003).
Dai, X. & Zhao, P. X. psRNATarget: a plant small RNA target analysis server. Nucl. Acids Res. 39, W155–W159 (2011).
Darzentas, N. Circoletto visualizing sequence similarity with Circos. Bioinformatics 26, 20 (2010).
Krzywinski, M. et al. Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639–1645 (2009).
Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874 (2016).
Kong, L. et al. CPC, assess the protein coding potential of transcripts using sequence features and support vector machine. Nucl. Acids Res. 35, 345-W349 (2007).
Yi, X., Zhang, Z., Ling, Y., Xu, W. & Su, Z. PNRD: a plant non-coding RNA database. Nucl. Acids Res. 43, D982–D989 (2015).
Singh, J. & Nagaraju, J. In silico prediction and characterization of microRNAs from red flour beetle (Tribolium castaneum). Insect Mol. Biol. 17, 427–436 (2008).
Zhang, B. H., Pan, X. P., Cox, S. B., Cobb, G. P. & Anderson, T. A. Evidence that miRNAs are different from other RNAs. Cell. Mol. Life Sci. 63, 246–254 (2006).
Wang, J., Hou, X. & Yang, X. Identification of conserved microRNAs and their targets in Chinese cabbage (Brassica rapa subsp. pekinensis). Genome. 52, 1029–1040 (2011).
Ambros, V., Bartel, B. & Bartel, D. P. A uniform system for microRNA annotation. RNA 9, 277–279 (2003).
Dereeper, V. et al. Phylogeny. Fr: robust phylogenetic analysis for the non-specialist. Nucl. Acids Res. 36, W465–W469 (2008).
Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl. Acids Res. 32, 1792–1797 (2004).
Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17, 540–552 (2000).
Guindon, S. & Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52, 696–704 (2003).
Chevenet, F., Brun, C., Banuls, A. L., Jacq, B. & Chisten, R. TreeDyn: towards dynamic graphics and annotations for analyses of trees. BMC Bioinform. 7, 439 (2006).
Lohse, M. et al. Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data. Plant Cell Environ. 37, 1250–1258 (2014).
Tian, T. et al. agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucl. Acids Res. 45, W122–W129 (2017).
Conesa, A. & Götz, S. Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int. J. Plant Genomics. https://doi.org/10.1155/2008/619832 (2008).
Bonnet, E., Wuyts, J., Rouze, P. & Van de Peer, Y. Detection of 91 potential conserved plant microRNAs in Arabidopsis thaliana and Oryza sativaidentifies important target genes. Proc. Natl. Acad. Sci. USA 101, 11511–11516 (2004).
Bonnet, E., Wuyts, J., Rouze, P. & Van de Peer, Y. Evidence that microRNA precursors, unlike other non-coding RNAs, have lower folding free energies than random sequences. Bioinformatics 20, 2911–2917 (2004).
Yin, V. P. et al. Fgf dependent depletion of microRNA-133 promotes appendage regeneration in zebrafish. Genes Dev. 22, 728–733 (2008).
Sunkar, R., Li, Y. F. & Jagadeeswaran, G. Functions of microRNAs in plant stress responses. Trends Plant Sci. 17, 196–203 (2012).
Zhang, B. MicroRNA: a new target for improving plant tolerance to abiotic stress. J. Exp. Bot. 66, 1749–1761 (2015).
Franco, N. B. et al. The micro-RNA172c-APETALA2-1 node as a key regulator of the common bean-rhizobium etli Nitrogen fixation symbiosis. Plant Physiol. 168, 273–291 (2015).
Kim, V. N. MicroRNA biogenesis: coordinated cropping and dicing. Nat. Rev. Mol. Cell Biol. 6, 376–385 (2005).
Chen, X. MicroRNA biogenesis and function in plants. FEBS Lett. 579, 5923–5931 (2005).
Kang, K. et al. A novel real-time PCR assay of microRNAs using S-Poly (T), a specific oligo (dT) reverse transcription primer with excellent sensitivity and specificity. PLoS ONE 7(11), e48536 (2012).
Zhang, H. et al. miRNA–mRNA integrated analysis reveals roles for miRNAs in a typical halophyte, Reaumuria soongorica, during seed germination under salt stress. Plants 9(3), 351 (2020).
Yasin, J.K. High density SSR and SNP saturated physical maps of Vigna radiata. https://legumeinfo.org/genomes/gbrowse/Vr1.0 (2018).
Singh, S. et al. A 62K genic-SNP chip array for genetic studies and breeding applications in pigeonpea (Cajanus cajan L. Millsp.). Sci. Rep. 10, 4960. https://doi.org/10.1038/s41598-020-61889-0 (2020).
Rogers, K. & Chen, X. Biogenesis, turnover, and mode of action of plant microRNAs. Plant Cell. 25, 2383–2399 (2013).
Dezulian, T., Palatnik, J. F., Huson, D. & Weigel, D. Conservation and divergence of microRNA families in plants. Genome Biol. 6, 1 (2005).
Dezulian, T., Schaefer, M., Wiese, R., Weigel, D. & Huson, D. H. CrossLink: visualization and exploration of sequence relationships between (micro) RNAs. Nucl. Acids Res. 34, W400–W404 (2006).
Weber, M. J. New human and mouse microRNA genes found by homology search. FEBS J. 272, 59–73 (2005).
Singh, N., Srivastava, S. & Sharma, A. Identification and analysis of miRNAs and their targets in ginger using bioinformatics approach. Gene 575, 570–576 (2016).
Panda, D. et al. Computational identification and characterization of conserved miRNAs and their target genes in garlic (Allium sativum L.) expressed sequence tags. Gene. 537, 333–342 (2014).
Akter, M. M. et al. Computational identification of miRNA and targets from expressed sequence tags of coffee (Coffea arabica). Saudi J. Biol. Sci. 21, 3–12 (2014).
Das, A. & Mondal, T. K. Computational identification of conserved microRNAs and their targets in tea (Camellia sinensis). Am. J. Plant Sci. 1, 77 (2010).
Zhang, B. H., Pan, X. P., Cox, S. B., Cobb, G. P. & Anderson, T. A. Conservation and divergence of plant microRNA genes. Plant J. 46, 243–259 (2006).
Wang, L., Liu, H., Li, D. & Chen, H. Identification and characterization of maize microRNAs involved in the very early stage of seed germination. BMC Genomics. 12, 154 (2011).
Stark, K. L. et al. Altered brain microRNA biogenesis contributes to phenotypic deficits in a 22q11-deletion mouse model. Nat. Genet. 40, 751–760 (2008).
Xie, Y. F., Shu, R., Jiang, S. Y., Liu, D. L. & Zhang, X. L. Comparison of microRNA profiles of human periodontal diseased and healthy gingival tissues. Int. J. Oral Sci. 3, 125–134 (2011).
Frazier, T. P., Xie, F., Freistaedter, A., Burklew, C. E. & Zhang, B. Identification and characterization of microRNAs and their target genes in tobacco (Nicotiana tabacum). Planta 232, 1289–1308 (2010).
Jha, A. & Shankar, R. Employing machine learning for reliable miRNA target identification in plants. BMC Genomics 12, 1 (2011).
Ding, J., Li, D., Ohler, U., Guan, J. & Zhou, S. Genome-wide search for miRNA-target interactions in Arabidopsis thaliana with an integrated approach. BMC Genomics 13, 1 (2012).
Kamthan, A., Chaudhuri, A., Kamthan, M. & Datta, A. Small RNAs in plants: recent development and application for crop improvement. Front. Plant Sci. 6, 208 (2015).
Li, C. & Zhang, B. MicroRNAs in control of plant development. J. Cell. Physiol. 231, 303–313 (2016).
Catalano, D., Pignone, D., Sonnante, G. & Finetti-Sialer, M. M. In-silico and in-vivo analyses of EST databases unveil conserved miRNAs from Carthamus tinctorius and Cynara cardunculus. BMC Bioinform. 13(S4), S12 (2012).
Rhoades, M. W. J. & Bartel, D. P. Computational identification of plant microRNAs and their targets, including a stress- induced miRNA. Mol. Cell. 14, 787–799 (2004).
Bhardwaj, S., Singh, A. & Singh, P. MicroRNA-based cancer therapeutics: big hope from small RNAs. Mol. Cell. Pharmacol. 2, 213–219 (2010).
Zhai, J. et al. MicroRNAs as master regulators of the plant NB-LRR defense gene family via the production of phased, trans-acting siRNAs. Genes Dev. 25, 2540–2553 (2011).
Xia, R., Zhu, H., An, Y. Q., Beers, E. P. & Liu, Z. Apple miRNAs and tasiRNAs with novel regulatory networks. Genome Biol. 13, 1 (2012).
Adai, A. et al. Computational prediction of miRNAs in Arabidopsis thaliana. Genome Res. 15, 78–91 (2005).
Zhang, L. et al. A genome-wide characterization of microRNA genes in maize. PLoS Genet. 5, e1000716 (2009).
Joshi, T. et al. Prediction of novel miRNAs and associated target genes in Glycine max. BMC Bioinform. 11, S14 (2010).
Szklarczyk, D. et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucl. Acids Res. 47(D1), D607–D613 (2019).
Xu, L. et al. OrthoVenn2: a web server for whole-genome comparison and annotation of orthologous clusters across multiple species. Nucleic Acids Res. 47(W1), W52–W58 (2019).
Babicki, S. et al. Heatmapper: web-enabled heat mapping for all. Nucl. Acids Res. 44(W1), W147–W153 (2016).
Acknowledgements
We thank Director, Indian Council of Agricultural Research -National Bureau of Plant Genetic Resources and Officer-In-Charge, Division of Genomic Resources for providing funds (PGR/DFP-BUR-DEL-01.01) and facilities to conduct this research. The authors extend their appreciation to King Saud University, Researchers Supportig Project (RSP-2020/118).
Funding
This work was supported by the Indian Council of Agricultural Research –National Bureau of Plant Genetic Resources [PGR/DFP-BUR-DEL-01.01] and funded by the King Saud University through Researchers Supporting Project number (RSP-2020/118).
Author information
Authors and Affiliations
Contributions
Y.J.K. the corresponding author; Y.J.K., A.P. and V.C. developed the concept, planned, coordinated and executed the research. Y.J.K., B.K.M. carried out research. S.H.W. and Y.J.K. and V.C. drafted the manuscript; Y.J.K. conducted pot experiments and wet lab experiments; M.A.P. performed quantitative expression of miRs; N.V., S.H.W., E.H.O., E.D.O. and P.S.P. edited the manuscript. Y.J.K. modified the manuscript. All authors read, approved and accepted the manuscript submission.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yasin, J.K., Mishra, B.K., Pillai, M.A. et al. Genome wide in-silico miRNA and target network prediction from stress responsive Horsegram (Macrotyloma uniflorum) accessions. Sci Rep 10, 17203 (2020). https://doi.org/10.1038/s41598-020-73140-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-73140-x
This article is cited by
-
Navigating Through Harsh Conditions: Coordinated Networks of Plant Adaptation to Abiotic Stress
Journal of Plant Growth Regulation (2024)
-
Comprehensive in silico analysis of the underutilized crop tef (Eragrostis tef (Zucc.) Trotter) genome reveals drought tolerance signatures
BMC Plant Biology (2023)
-
Orphan legumes: harnessing their potential for food, nutritional and health security through genetic approaches
Planta (2022)
-
Integrative small RNA and transcriptome analysis provides insight into key role of miR408 towards drought tolerance response in cowpea
Plant Cell Reports (2022)
-
Physical map of lncRNAs and lincRNAs linked with stress responsive miRs and genes network of pigeonpea (Cajanus cajan L.)
Journal of Plant Biochemistry and Biotechnology (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.