Identification of miRNAs in sorghum by using bioinformatics approach

Comprehensive Genome-Wide Identification and Expression Profiling of Eceriferum (CER) Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions

Plant surfaces are covered with cuticle wax and are the first barrier between a plant and environmental stresses. Eceriferum (CER) is an important gene family involved in wax biosynthesis and stress resistance. In this study, for the first time, 34 CER genes were identified in the passion fruit (Passiflora edulis) genome, and PeCER proteins varied in physicochemical properties. A phylogenetic tree was constructed and divided into seven clades to identify the evolutionary relationship with other plant species. Gene structure analyses revealed that conserved motifs ranged from 1 to 24, and that exons ranged from 1 to 29. The cis-element analysis provides insight into possible roles of PeCER genes in plant growth, development and stress responses. The syntenic analysis revealed that segmental (six gene pairs) and tandem (six gene pairs) gene duplication played an important role in the expansion of PeCER genes and underwent a strong purifying selection. In addition, 12 putative ped-miRNAs were identified to be targeting 16 PeCER genes, and PeCER6 was the most targeted by four miRNAs including ped-miR157a-5p, ped-miR164b-5p, ped-miR319b, and ped-miR319l. Potential transcription factors (TFs) such as ERF, AP2, MYB, and bZIP were predicted and visualized in a TF regulatory network interacting with PeCER genes. GO and KEGG annotation analysis revealed that PeCER genes were highly related to fatty acid, cutin, and wax biosynthesis, plant-pathogen interactions, and stress response pathways. The hypothesis that most PeCER proteins were predicted to localize to the plasma membrane was validated by transient expression assays of PeCER32 protein in onion epidermal cells. qRT-PCR expression results showed that most of the PeCER genes including PeCER1, PeCER11, PeCER15, PeCER17, and PeCER32 were upregulated under drought and Fusarium kyushuense stress conditions compared to controls. These findings provide a foundation for further studies on functions of PeCER genes to further facilitate the genetic modification of passion fruit wax biosynthesis and stress resistance.

Genome-Wide Identification and Expression Profiling of KCS Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions

Plant and fruit surfaces are covered with cuticle wax and provide a protective barrier against biotic and abiotic stresses. Cuticle wax consists of very-long-chain fatty acids (VLCFAs) and their derivatives. β-Ketoacyl-CoA synthase (KCS) is a key enzyme in the synthesis of VLCFAs and provides a precursor for the synthesis of cuticle wax, but the KCS gene family was yet to be reported in the passion fruit (Passiflora edulis). In this study, thirty-two KCS genes were identified in the passion fruit genome and phylogenetically grouped as KCS1-like, FAE1-like, FDH-like, and CER6-like. Furthermore, thirty-one PeKCS genes were positioned on seven chromosomes, while one PeKCS was localized to the unassembled genomic scaffold. The cis-element analysis provides insight into the possible role of PeKCS genes in phytohormones and stress responses. Syntenic analysis revealed that gene duplication played a crucial role in the expansion of the PeKCS gene family and underwent a strong purifying selection. All PeKCS proteins shared similar 3D structures, and a protein-protein interaction network was predicted with known Arabidopsis proteins. There were twenty putative ped-miRNAs which were also predicted that belong to nine families targeting thirteen PeKCS genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation results were highly associated with fatty acid synthase and elongase activity, lipid metabolism, stress responses, and plant-pathogen interaction. The highly enriched transcription factors (TFs) including ERF, MYB, Dof, C2H2, TCP, LBD, NAC, and bHLH were predicted in PeKCS genes. qRT-PCR expression analysis revealed that most PeKCS genes were highly upregulated in leaves including PeKCS2, PeKCS4, PeKCS8, PeKCS13, and PeKCS9 but not in stem and roots tissues under drought stress conditions compared with controls. Notably, most PeKCS genes were upregulated at 9th dpi under Fusarium kyushuense biotic stress condition compared to controls. This study provides a basis for further understanding the functions of KCS genes, improving wax and VLCFA biosynthesis, and improvement of passion fruit resistance.

Comparative analysis of sRNAs, degradome and transcriptomics in sweet sorghum reveals the regulatory roles of miRNAs in Cd accumulation and tolerance

Key miRNAs including sbi-miR169p/q, sbi-miR171g/j, sbi-miR172a/c/d, sbi-miR172e, sbi-miR319a/b, sbi-miR396a/b, miR408, sbi-miR5384, sbi-miR5565e and nov_23 were identified to function in the regulation of Cd accumulation and tolerance. As an energy plant, sweet sorghum shows great potential in the phytoremediation of Cd-contaminated soils. However, few studies have focused on the regulatory roles of miRNAs and their targets under Cd stress. In this study, comparative analysis of sRNAs, degradome and transcriptomics was conducted in high-Cd accumulation (H18) and low-Cd accumulation (L69) genotypes of sweet sorghum. A total of 38 conserved and 23 novel miRNAs with differential expressions were identified under Cd stress or between H18 and L69, and 114 target genes of 41 miRNAs were validated. Furthermore, 25 miRNA-mRNA pairs exhibited negatively correlated expression profiles and sbi-miR172e together with its target might participate in the distinct Cd tolerance between H18 and L69 as well as sbi-miR172a/c/d. Additionally, two groups of them: miR169p/q-nov_23 and miR408 were focused through the co-expression analysis, which might be involved in Cd uptake and tolerance by regulating their targets associated with transmembrane transportation, cytoskeleton activity, cell wall construction and ROS (reactive oxygen species) homeostasis. Further experiments exhibited that cell wall components of H18 and L69 were different when exposed to cadmium, which might be regulated by miR169p/q, miR171g/j, miR319a/b, miR396a/b, miR5384 and miR5565e through their targets. Through this study, we aim to reveal the potential miRNAs involved in sweet sorghum in response to Cd stress and provide references for developing high-Cd accumulation or high Cd-resistant germplasm of sweet sorghum that can be used in phytoremediation.

Regulatory protein genes and microRNAs in response to selenium stimuli in Pueraria lobata (Willd.) Ohwi

Regulatory protein genes and microRNAs (miRNAs) play important roles in response to abiotic and biotic stress, and the biosynthesis of secondary metabolites in plants. However, their responses to selenium (Se) stimuli have not been comprehensively studied in Pueraria lobata (Willd.) Ohwi, a selenocompound-rich medicinal and edible plant. In this study, we identified a total of 436/556/1161/624 transcription factors, 134/157/308/172 transcriptional regulators, and 341/456/250/518 protein kinases, which were co-expressed with at least one selenocompound-related structural gene/sulfate transporter or phosphate transporter/reactive oxygen species (ROS) scavenging structural gene/isoflavone-related structural gene, respectively. Then, we identified a total of 87 expressed miRNAs by Se disposure, in which 11 miRNAs, including miR171f-3p, miR390b-3P, miR-N111b, miR-N118, miR-N30, miR-N38-3P, miR-N61a, miR-N61b, miR-N80-3p, miR-N84-3P, and miR-N90.2-3P, were significantly upregulated. We also identified a total of 1172 target genes for the 87 expressed miRNAs. Gene Ontology enrichment analysis of these target genes showed that regulation of transcription, DNA-templated, integral component of membrane, nucleus, ATP binding, and plasma membrane are the top five subclassifications. Finally, we revealed that 5 miRNAs targeted 10 regulatory protein genes, which are highly correlated with at least one selenocompound-related structural gene or transporter gene; 5 miRNAs targeted 10 regulatory protein genes, which are highly correlated with at least one ROS scavenging structural gene; and 5 miRNAs targeted 9 regulatory protein genes, which are potentially involved in the isoflavone biosynthesis. Overall, the study provides us the comprehensive insight into the roles of regulatory proteins and miRNAs in response to Se stimuli in P. lobata.

Identification of genuine and novel miRNAs in Amaranthus hypochondriacus from high-throughput sequencing data

Amaranth has been proposed as an exceptional alternative for food security and climate change mitigation. Information about the distribution, abundance, or specificity of miRNAs in amaranth species is scare. Here, small RNAs from seedlings under control, drought, heat, and cold stress conditions of the Amaranthus hypocondriacus variety "Gabriela" were sequenced and miRNA loci identified in the amaranth genome using the ShortStack software. Fifty-three genuine miRNA clustersthirty-nine belonging to conserved families, and fourteen novel, were identified. Identification of their target genes suggests that conserved amaranth miRNAs are involved in growth and developmental processes, as well as stress responses. MiR0005, an amaranth-specific miRNA, exhibited an unusual high level of expression, akin to that of conserved miRNAs. Overall, our results broaden our knowledge regarding the distribution, abundance and expression of miRNAs in amaranth, providing the basis for future research on miRNAs and their functions in this important species.

Potential Prophylactic Effect of Berberine against Rat Colon Carcinoma Induce by 1,2-Dimethyl Hydrazine

Introduction: Colon Cancer remains one of the major worldwide causes of cancer related morbidity and mortality in both genders. Berberine (BBR), a major component of alkaloids that possess a variety of pharmacological properties. Objective: This study shows the ameliorating roles of berberine on 1,2 Di methyl hydrazine (DMH) induced colon cancer in male Swiss albino rats. Methods: The rats were segregated into four groups: group 1, control rats; group 2, rats were orally received berberine (75 mg/kg b.wt./day) daily for ten weeks; group 3,rats were subcutaneously injected with DMH (20 mg/kg b.wt) once a week for 8 weeks ,group 4, rats were treated firstly with berberine for two weeks before DMH intoxication and concurrently with DMH over 8 weeks. Result: DMH injection decreased the antioxidants levels (GSH and SOD) and increased inflammatory markers (MPO, MAPK and COX-2). Moreover, it downregulated apoptotic markers (Caspase-3 and P53) expression that confirmed by colon cell proliferation. The prophylactic effect of berberine was noticed as its pre-and co-administration increased antioxidants status and apoptotic markers expression that associated with inflammatory markers down-regulation with absence of proliferated colon cells. Conclusion: Therefore, the overall findings proved that the anti-proliferative effect of berberine return to its antioxidants and anti-inflammatory properties that activated the programmed cell death process.

Abiotic Stress Responsive miRNA-Target Network and Related Markers (SNP, SSR) in Brassica juncea

Abiotic stress is one of the major factors responsible for huge yield loss in crop plants. MicroRNAs play a key role in adaptive responses of plants under abiotic stress conditions through post-transcriptional gene regulations. In present study, 95 potential miRNAs were predicted in Brassica juncea using comparative genomics approach. It was noted that these miRNAs, target several transcription factors (TFs), transporter family proteins, signaling related genes, and protease encoding genes. Nineteen distinct miRNA-target regulatory networks were observed with significant involvement in regulation of transcription, response to stimulus, hormone and auxin mediated signaling pathway related gene ontology (GO) term. The sucrose-starch metabolism and pentose-gluconate interconversion pathways were found significantly enriched for these target genes. Molecular markers such as Simple Sequence Repeats (SSR) and Single Nucleotide Polymorphism (SNPs) were identified on miRNAs (miR-SSRs and miR-SNPs) and their target genes in B. juncea. Notably, one of the miR-SNP (C/T) was found at the 5^th position on mature region of miR2926. This C/T transition led to the distorted and unstable hairpin structure of miR2926, consequently complete loss of target function. Hence, findings from this study will lay a foundation for marker assisted breeding for abiotic stress tolerant varieties of B. juncea.

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.

Large-scale bioinformatic analysis of the regulation of the disease resistance NBS gene family by microRNAs in Poaceae

In the present study, we have screened 71, 713, 525, 119 and 241 mature miRNA variants from Hordeum vulgare, Oryza sativa, Brachypodium distachyon, Triticum aestivum, and Sorghum bicolor, respectively, and classified them with respect to their conservation status and expression levels. These Poaceae non-redundant miRNA species (1,669) were distributed over a total of 625 MIR families, among which only 54 were conserved across two or more plant species, confirming the relatively recent evolutionary differentiation of miRNAs in grasses. On the other hand, we have used 257 H. vulgare, 286T. aestivum, 119 B. distachyon, 269 O. sativa, and 139 S. bicolor NBS domains, which were either mined directly from the annotated proteomes, or predicted from whole genome sequence assemblies. The hybridization potential between miRNAs and their putative NBS genes targets was analyzed, revealing that at least 454 NBS genes from all five Poaceae were potentially regulated by 265 distinct miRNA species, most of them expressed in leaves and predominantly co-expressed in additional tissues. Based on gene ontology, we could assign these probable miRNA target genes to 16 functional groups, among which three conferring resistance to bacteria (Rpm1, Xa1 and Rps2), and 13 groups of resistance to fungi (Rpp8,13, Rp3, Tsn1, Lr10, Rps1-k-1, Pm3, Rpg5, and MLA1,6,10,12,13). The results of the present analysis provide a large-scale platform for a better understanding of biological control strategies of disease resistance genes in Poaceae, and will serve as an important starting point for enhancing crop disease resistance improvement by means of transgenic lines with artificial miRNAs.

MicroRNA expression profiles in response to drought stress in Sorghum bicolor

The regulatory role of small non-coding RNAs that are 20-24 nucleotides in length has become the foremost area of research for biologists. A major class of small RNAs represented by the microRNAs (miRNAs), has been implicated in various aspects of plant development including leaf pattering, meristem function, root patterning etc. Recent findings support that miRNAs are regulated by drought and other abiotic stresses in various plant species. In this study, were report the expression profiling of 8 known abiotic stress deregulated miRNAs in 11 elite sorghum genotypes, under watered and drought conditions. Significant deregulation was observed with miR396, miR393, miR397-5p, miR166, miR167 and miR168. Among these, the expression levels of sbi-miR396 and sbi-miR398 were the highest in all the genotypes. The expression of sbi-miR396 was maximum in the grain sorghum HSD3226 under well-watered conditions and the profile shifted towards HSD3221 under drought stress. Forage accessions, N98 and Atlas, showed an opposite behavior in expression patterns of miR397-5p in drought physiologies. Such dynamic expression patterns could be indicative of prevailing drought tolerant mechanisms present in these sorghum accessions. This data provides insights into sorghum miRNAs which may have potential use in improving drought tolerance in sorghum and other cereal crops.

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

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

Manipulating microRNAs for improved biomass and biofuels from plant feedstocks

Petroleum-based fuels are nonrenewable and unsustainable. Renewable sources of energy, such as lignocellulosic biofuels and plant metabolite-based drop-in fuels, can offset fossil fuel use and reverse environmental degradation through carbon sequestration. Despite these benefits, the lignocellulosic biofuels industry still faces many challenges, including the availability of economically viable crop plants. Cell wall recalcitrance is a major economic barrier for lignocellulosic biofuels production from biomass crops. Sustainability and biomass yield are two additional, yet interrelated, foci for biomass crop improvement. Many scientists are searching for solutions to these problems within biomass crop genomes. MicroRNAs (miRNAs) are involved in almost all biological and metabolic process in plants including plant development, cell wall biosynthesis and plant stress responses. Because of the broad functions of their targets (e.g. auxin response factors), the alteration of plant miRNA expression often results in pleiotropic effects. A specific miRNA usually regulates a biologically relevant bioenergy trait. For example, relatively low miR156 overexpression leads to a transgenic feedstock with enhanced biomass and decreased recalcitrance. miRNAs have been overexpressed in dedicated bioenergy feedstocks such as poplar and switchgrass yielding promising results for lignin reduction, increased plant biomass, the timing of flowering and response to harsh environments. In this review, we present the status of miRNA-related research in several major biofuel crops and relevant model plants. We critically assess published research and suggest next steps for miRNA manipulation in feedstocks for increased biomass and sustainability for biofuels and bioproducts.

Profiling microRNA expression during multi-staged date palm (Phoenix dactylifera L.) fruit development

MicroRNAs (miRNAs) play crucial roles in multiple stages of plant development and regulate gene expression at posttranscriptional and translational levels. In this study, we first identified 238 conserved miRNAs in date palm (Phoenix dactylifera) based on a high-quality genome assembly and defined 78 fruit-development-associated (FDA) miRNAs, whose expression profiles are variable at different fruit development stages. Using experimental data, we subsequently detected 276 novel P. dactylifera-specific FDA miRNAs and predicted their targets. We also revealed that FDA miRNAs function mainly in regulating genes involved in starch/sucrose metabolisms and other carbon metabolic pathways; among them, 221 FDA miRNAs exhibit negative correlation with their corresponding targets, which suggests their direct regulatory roles on mRNA targets. Our data define a comprehensive set of conserved and novel FDA miRNAs along with their expression profiles, which provide a basis for further experimentation in assigning discrete functions of these miRNAs in P. dactylifera fruit development.

Identification of novel drought-responsive microRNAs and trans-acting siRNAs from Sorghum bicolor (L.) Moench by high-throughput sequencing analysis

Small non-coding RNAs (sRNAs) namely microRNAs (miRNAs) and trans-acting small interfering RNAs (tasi-RNAs) play a crucial role in post-transcriptional regulation of gene expression and thus the control plant development and stress responses. In order to identify drought-responsive miRNAs and tasi-RNAs in sorghum, we constructed small RNA libraries from a drought tolerant (M35-1) and susceptible (C43) sorghum genotypes grown under control and drought stress conditions, and sequenced by Illumina Genome Analyzer IIx. Ninety seven conserved and 526 novel miRNAs representing 472 unique miRNA families were identified from sorghum. Ninety-six unique miRNAs were found to be regulated by drought stress, of which 32 were up- and 49 were down-regulated (fold change ≥ 2 or ≤ -2) at least in one genotype, while the remaining 15 miRNAs showed contrasting drought-regulated expression pattern between genotypes. A maximum of 17 and 18 miRNAs was differentially regulated under drought stress condition in the sensitive and tolerant genotypes, respectively. These results suggest that genotype dependent stress responsive regulation of miRNAs may contribute, at least in part, to the differential drought tolerance of sorghum genotypes. We also identified two miR390-directed TAS3 gene homologs and the auxin response factors as tasi-RNA targets. We predicted more than 1300 unique target genes for the novel and conserved miRNAs. These target genes were predicted to be involved in different cellular, metabolic, response to stimulus, biological regulation, and developmental processes. Genome-wide identification of stress-responsive miRNAs, tasi-RNAs and their targets identified in this study will be useful in unraveling the molecular mechanisms underlying drought stress responses and genetic improvement of biomass production and stress tolerance in sorghum.

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

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

History and current status of wheat miRNAs using next-generation sequencing and their roles in development and stress

As small molecules that aid in posttranscriptional silencing, microRNA (miRNA) discovery and characterization have vastly benefited from the recent development and widespread application of next-generation sequencing (NGS) technologies. Several miRNAs were identified through sequencing of constructed small RNA libraries, whereas others were predicted by in silico methods using the recently accumulating sequence data. NGS was a major breakthrough in efforts to sequence and dissect the genomes of plants, including bread wheat and its progenitors, which have large, repetitive and complex genomes. Availability of survey sequences of wheat whole genome and its individual chromosomes enabled researchers to predict and assess wheat miRNAs both in the subgenomic and whole genome levels. Moreover, small RNA construction and sequencing-based studies identified several putative development- and stress-related wheat miRNAs, revealing their differential expression patterns in specific developmental stages and/or in response to stress conditions. With the vast amount of wheat miRNAs identified in recent years, we are approaching to an overall knowledge on the wheat miRNA repertoire. In the following years, more comprehensive research in relation to miRNA conservation or divergence across wheat and its close relatives or progenitors should be performed. Results may serve valuable in understanding both the significant roles of species-specific miRNAs and also provide us information in relation to the dynamics between miRNAs and evolution in wheat. Furthermore, putative development- or stress-related miRNAs identified should be subjected to further functional analysis, which may be valuable in efforts to develop wheat with better resistance and/or yield.

Bioinformatic identification and experimental validation of miRNAs from foxtail millet (Setaria italica)

MiRNAs are a novel group of non-coding small RNAs that negatively regulate gene expression. Many miRNAs have been identified and investigated extensively in plant species with sequenced genomes. However, few miRNAs have been identified in foxtail millet (Setaria italica), which is an ancient cereal crop of great importance for dry land agriculture. In this study, 271 foxtail millet miRNAs belonging to 44 families were identified using a bioinformatics approach. Twenty-three pairs of sense/antisense miRNAs belonging to 13 families, and 18 miRNA clusters containing members of 8 families were discovered in foxtail millet. We identified 432 potential targets for 38 miRNA families, most of which were predicted to be involved in plant development, signal transduction, metabolic pathways, disease resistance, and environmental stress responses. Gene ontology (GO) analysis revealed that 101, 56, and 23 target genes were involved in molecular functions, biological processes, and cellular components, respectively. We investigated the expression patterns of 43 selected miRNAs using qRT-PCR analysis. All of the miRNAs were expressed ubiquitously with many exhibiting different expression levels in different tissues. We validated five predicted targets of four miRNAs using the RNA ligase mediated rapid amplification of cDNA end (5'-RLM-RACE) method.

Identification of conserved microRNAs and their targets in chickpea (Cicer arietinum L.)

The microRNAs (miRNAs) are a new class of non-protein coding small RNAs that regulate gene expression at the post-transcriptional level in plants. Although thousands of miRNAs have been identified in many plant species, little studies have been reported about chickpea microRNAs. In this study, 28 potential miRNA candidates belonging to 20 families were identified from 16 ESTs and 12 GSSs in chickpea using a comparative genome-based computational analysis. A total of 664 miRNA targets were predicted and some of them encoded transcription factors as well as genes that function in stress response, signal transduction, methylation and a variety of other metabolic processes. These findings lay the foundation for further understanding of miRNA function in the development of chickpea.

MicroRNAs and their cross-talks in plant development

Plant development is a complex process influenced by exogenous and endogenous elements. A series of postembryonic developmental events is involved to form the final architecture and contend with the changing environment. MicroRNA (miRNA) is one of endogenous non-coding RNAs, which plays an important role in plant developmental regulation. In this review, we summarized 34 miRNA families that are closely associated with plant development. Among these families, nine are expressed only in specific organs, whereas 20 families are expressed in at least two different organs. It is known that some miRNAs are expressed across most processes of plant growth, while some appear only at particular developmental stages or under special environmental conditions such as drought, waterlogging and short-day time. These miRNAs execute their diverse functions by regulating developmental gene expression levels, interacting with phytohormone signaling response, participating in the biogenesis of ta-siRNAs and affecting the production of miRNAs.