Computational identification of microRNAs and their targets in wheat (Triticum aestivum L.)

Expression Patterns of miR398, miR167, and miR159 in the Interaction between Bread Wheat (Triticum aestivum L.) and Pathogenic Fusarium culmorum and Beneficial Trichoderma Fungi

Bread wheat (Triticum aestivum L.) is an agronomically significant cereal cultivated worldwide. Wheat breeding is limited by numerous abiotic and biotic stresses. One of the most deleterious factors is biotic stress provoked by the Fusarium culmorum fungus. This pathogen is a causative agent of Fusarium root rot and Fusarium head blight. Beneficial fungi Trichoderma atroviride and T. cremeum are strong antagonists of mycotoxigenic Fusarium spp. These fungi promote plant growth and enhance their tolerance of negative environmental conditions. The aim of the study was to determine and compare the spatial (in above- and underground organs) and temporal (early: 6 and 22 hpi; and late: 5 and 7 dpi reactions) expression profiles of three mature miRNAs (miR398, miR167, and miR159) in wheat plants inoculated with two strains of F. culmorum (KF846 and EW49). Moreover, the spatial expression patterns in wheat response between plants inoculated with beneficial T. atroviride (AN35) and T. cremeum (AN392) were assessed. Understanding the sophisticated role of miRNAs in wheat-fungal interactions may initiate a discussion concerning the use of this knowledge to protect wheat plants from the harmful effects of fungal pathogens. With the use of droplet digital PCR (ddPCR), the absolute quantification of the selected miRNAs in the tested material was carried out. The differential accumulation of miR398, miR167, and miR159 in the studied groups was observed. The abundance of all analyzed miRNAs in the roots demonstrated an increase in the early and reduction in late wheat response to F. culmorum inoculation, suggesting the role of these particles in the initial wheat reaction to the studied fungal pathogen. The diverse expression patterns of the studied miRNAs between Trichoderma-inoculated or F. culmorum-inoculated plants and control wheat, as well as between Trichoderma-inoculated and F. culmorum-inoculated plants, were noticed, indicating the need for further analysis.

Elucidating the Functional Role of Predicted miRNAs in Post- Transcriptional Gene Regulation Along with Symbiosis in Medicago truncatula

Background:

microRNAs are small non-coding RNAs which inhibit translational and post-transcriptional processes whereas long non-coding RNAs are found to regulate both transcriptional and post-transcriptional gene expression. Medicago truncatula is a well-known model plant for studying legume biology and is also used as a forage crop. In spite of its importance in nitrogen fixation and soil fertility improvement, little information is available about Medicago non-coding RNAs that play important role in symbiosis.

Objective:

In this study we have tried to understand the role of Medicago ncRNAs in symbiosis and regulation of transcription factors.

Methods:

We have identified novel miRNAs by computational methods considering various parameters like length, MFEI, AU content, SSR signatures and tried to establish an interaction model with their targets obtained through psRNATarget server.

Results:

149 novel miRNAs are predicted along with their 770 target proteins. We have also shown that 51 of these novel miRNAs are targeting 282 lncRNAs.

Conclusion:

In this study role of Medicago miRNAs in the regulation of various transcription factors are elucidated. Knowledge gained from this study will have a positive impact on the nitrogen fixing ability of this important model plant, which in turn will improve the soil fertility.

Genome-wide identification and characterization of gene family for RWP-RK transcription factors in wheat (Triticum aestivum L.)

RWP-RKs represent a small family of transcription factors (TFs) that are unique to plants and function particularly under conditions of nitrogen starvation. These RWP-RKs have been classified in two sub-families, NLPs (NIN-like proteins) and RKDs (RWP-RK domain proteins). NLPs regulate tissue-specific expression of genes involved in nitrogen use efficiency (NUE) and RKDs regulate expression of genes involved in gametogenesis/embryogenesis. During the present study, using in silico approach, 37 wheat RWP-RK genes were identified, which included 18 TaNLPs (2865 to 7340 bp with 4/5 exons), distributed on 15 chromosomes from 5 homoeologous groups (with two genes each on 4B,4D and 5A) and 19 TaRKDs (1064 to 5768 bp with 1 to 6 exons) distributed on 12 chromosomes from 4 homoeologous groups (except groups 1, 4 and 5); 2–3 splice variants were also available in 9 of the 37 genes. Sixteen (16) of these genes also carried 24 SSRs (simple sequence repeats), while 11 genes had targets for 13 different miRNAs. At the protein level, MD simulation analysis suggested their interaction with nitrate-ions. Significant differences were observed in the expression of only two (TaNLP1 and TaNLP2) of the nine representative genes that were used for in silico expression analysis under varying levels of N at post-anthesis stage (data for other genes was not available for in silico expression analysis). Differences in expression were also observed during qRT-PCR, when expression of four representative genes (TaNLP2, TaNLP7, TaRKD6 and TaRKD9) was examined in roots and shoots of seedlings (under different conditions of N supply) in two contrasting genotypes which differed in NUE (C306 with low NUE and HUW468 with high NUE). These four genes for qRT-PCR were selected on the basis of previous literature, level of homology and the level of expression (in silico study). In particular, the TaNLP7 gene showed significant up-regulation in the roots and shoots of HUW468 (with higher NUE) during N-starvation; this gene has already been characterized in Arabidopsis and tobacco, and is known to be involved in nitrate-signal transduction pathway.

Computational Identification of MicroRNAs and Their Targets from Finger Millet (Eleusine coracana)

MicroRNAs are endogenous small RNAs regulating intrinsic normal growth and development of plant. Discovering miRNAs, their targets and further inferring their functions had become routine process to comprehend the normal biological processes of miRNAs and their roles in plant development. In this study, we used homology-based analysis with available expressed sequence tag of finger millet (Eleusine coracana) to predict conserved miRNAs. Three potent miRNAs targeting 88 genes were identified. The newly identified miRNAs were found to be homologous with miR166 and miR1310. The targets recognized were transcription factors and enzymes, and GO analysis showed these miRNAs played varied roles in gene regulation. The identification of miRNAs and their targets is anticipated to hasten the pace of key epigenetic regulators in plant development.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

MicroRNA-based biotechnology for plant improvement

MicroRNAs (miRNAs) are an extensive class of newly discovered endogenous small RNAs, which negatively regulate gene expression at the post-transcription levels. As the application of next-generation deep sequencing and advanced bioinformatics, the miRNA-related study has been expended to non-model plant species and the number of identified miRNAs has dramatically increased in the past years. miRNAs play a critical role in almost all biological and metabolic processes, and provide a unique strategy for plant improvement. Here, we first briefly review the discovery, history, and biogenesis of miRNAs, then focus more on the application of miRNAs on plant breeding and the future directions. Increased plant biomass through controlling plant development and phase change has been one achievement for miRNA-based biotechnology; plant tolerance to abiotic and biotic stress was also significantly enhanced by regulating the expression of an individual miRNA. Both endogenous and artificial miRNAs may serve as important tools for plant improvement.

Genome-wide discovery of novel and conserved microRNAs in white shrimp (Litopenaeus vannamei)

Of late years, a large amount of conserved and species-specific microRNAs (miRNAs) have been performed on identification from species which are economically important but lack a full genome sequence. In this study, Solexa deep sequencing and cross-species miRNA microarray were used to detect miRNAs in white shrimp. We identified 239 conserved miRNAs, 14 miRNA* sequences and 20 novel miRNAs by bioinformatics analysis from 7,561,406 high-quality reads representing 325,370 distinct sequences. The all 20 novel miRNAs were species-specific in white shrimp and not homologous in other species. Using the conserved miRNAs from the miRBase database as a query set to search for homologs from shrimp expressed sequence tags (ESTs), 32 conserved computationally predicted miRNAs were discovered in shrimp. In addition, using microarray analysis in the shrimp fed with Panax ginseng polysaccharide complex, 151 conserved miRNAs were identified, 18 of which were significant up-expression, while 49 miRNAs were significant down-expression. In particular, qRT-PCR analysis was also performed for nine miRNAs in three shrimp tissues such as muscle, gill and hepatopancreas. Results showed that these miRNAs expression are tissue specific. Combining results of the three methods, we detected 20 novel and 394 conserved miRNAs. Verification with quantitative reverse transcription (qRT-PCR) and Northern blot showed a high confidentiality of data. The study provides the first comprehensive specific miRNA profile of white shrimp, which includes useful information for future investigations into the function of miRNAs in regulation of shrimp development and immunology.

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

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

A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome

An ordered draft sequence of the 17-gigabase hexaploid bread wheat (Triticum aestivum) genome has been produced by sequencing isolated chromosome arms. We have annotated 124,201 gene loci distributed nearly evenly across the homeologous chromosomes and subgenomes. Comparative gene analysis of wheat subgenomes and extant diploid and tetraploid wheat relatives showed that high sequence similarity and structural conservation are retained, with limited gene loss, after polyploidization. However, across the genomes there was evidence of dynamic gene gain, loss, and duplication since the divergence of the wheat lineages. A high degree of transcriptional autonomy and no global dominance was found for the subgenomes. These insights into the genome biology of a polyploid crop provide a springboard for faster gene isolation, rapid genetic marker development, and precise breeding to meet the needs of increasing food demand worldwide.

Computational identification of novel microRNAs and targets in Glycine max

Plant miRNAs, the endogenous non-coding small RNAs of about 20-24 nucleotides, play important roles in multiple biological processes by acting as negative regulators of their targeted mRNAs. Soybean (Glycine max (L.) Merr.) is one of the important oil crops of the world, in which many miRNAs have been obtained through the computational prediction or experiments. However, the miRNA genes identified for soybean are still far from saturation, and their biological functions are largely unknown. Here, a total of 48 candidates of miRNAs were identified following a range of strict filtering criteria. Detailed sequence analysis showed that G. max pre-miRNAs vary in length from 47 to 380 nt, embody mature miRNAs that differ in their physical location within the pre-miRNAs. In this study, twenty miRNAs were confirmed by microarray and three miRNAs were further validated by poly(A)-tailed RT-PCR. Comparative sequence analysis of soybean miRNA sequences showed that uracil is the dominant base in the first position at the 5' end of the mature miRNAs, and the base may have an important functional role in miRNA biogenesis and/or miRNA-mediated gene regulation. Finally, we predicted potential targets of these miRNAs. These target genes were predicted to encode transcription factors, resistance protein, heat shock protein, protein kinase, transporter, zinc finger protein and others, which might play important roles in soybean development and stress response.

A comprehensive genome-wide study on tissue-specific and abiotic stress-specific miRNAs in Triticum aestivum

Productivity of wheat crop is largely dependent on its growth and development that, in turn, is mainly regulated by environmental conditions, including abiotic stress factors. miRNAs are key regulators of gene expression networks involved in diverse aspects of development and stress responses in plants. Using high-throughput sequencing of eight small RNA libraries prepared from diverse abiotic stresses and tissues, we identified 47 known miRNAs belonging to 20 families, 49 true novel and 1030 candidate novel miRNAs. Digital gene expression analysis revealed that 257 miRNAs exhibited tissue-specific expression and 74 were associated with abiotic stresses. Putative target genes were predicted for miRNAs identified in this study and their grouping into functional categories indicated that the putative targets were involved in diverse biological processes. RLM-RACE of predicted targets of three known miRNAs (miR156, miR160 and miR164) confirmed their mRNA cleavage, thus indicating their regulation at post-transcriptional level by the corresponding miRNAs. Mapping of the sequenced data onto the wheat progenitors and closely related monocots revealed a large number of evolutionary conserved miRNAs. Additional expression profiling of some of these miRNAs in other abiotic stresses underline their involvement in multiple stresses. Our findings provide valuable resource for an improved understanding of the role of miRNAs in stress tolerance as well as plant development.

Computational identification of microRNAs and their targets in perennial Ryegrass (Lolium perenne)

MicroRNAs (miRNAs) are small non-coding RNA molecules of 22 nucleotides in length that have been characterized as regulators of messenger RNA (mRNA) regulating a number of developmental processes in plants and animals by silencing genes using multiple mechanisms. miRNAs have been extensively studied in various plant species; however, few information are available about miRNAs in perennial ryegrass, animal feed, and industrial raw materials. In this study, the 12 potential perennial ryegrass miRNAs were identified for the first time by computational approach. Using the newly identified miRNA sequences, the perennial ryegrass mRNA database was further used for BLAST search and detected 33 potential targets of miRNAs. Prediction of potential miRNA target genes revealed their functions involved in various important plant biological processes. Our result should be useful for further investigation into the biological functions of miRNAs in perennial ryegrass. The selected miRNAs representing four families were verified by RT-PCR experiment, indicating that the prediction method that we used to identify the miRNAs was effective.

Identification of microRNA genes and their mRNA targets in Festuca arundinacea

MicroRNAs (miRNAs) have emerged as a novel class of endogenous, small, non-coding RNAs of 22 nucleotides (nts) in length, which plays important roles in post-transcriptional degradation of target mRNA or inhibition of protein synthesis through binding the specific sites of target mRNA. Growing evidences have shown that miRNAs play an important role in various biological processes, including growth and development, signal transduction, apoptosis, proliferation, stress responses, maintenance of genome stability, and so on. In our study, we used bioinformatic tools to predict miRNA and the corresponding target genes of Festuca arundinacea. We used known miRNAs of other plants from miRBase to search against expressed sequence tags (EST) databases and genome survey sequences (GSS) of F. arundinacea. A total of 8 potential miRNAs were predicted. Phylogenetic analysis of the predicted miRNAs revealed that miRNA398c of F. arundinacea species was evolutionary highly conserved with Populus trichocarpa. The 8 potential miRNAs corresponding to 20 target genes were found. Most of the miRNA target genes were predicted to encode transcription factors that regulate cell growth and development, signaling, metabolism, and other biology processes. By bioinformatics methods, we can effectively predict novel miRNAs and its target genes and add information to F. arundinacea miRNA database. Moreover, it shows a path for the prediction and analysis of miRNAs to those species whose genomes are not available through bioinformatics tools.

Computational identification of miRNA and targets from expressed sequence tags of coffee (Coffea arabica)

MicroRNAs (miRNAs) are the group of ∼22 nucleotides long noncoding small endogenous and evolutionary conserved post-transcriptional regulatory RNAs, which show an enormous role in various biological and metabolic processes in both animals and plants. To date not a single miRNA has been identified in coffee (Coffea arabica), which is an economically important plant of Rubiaceae family. In this study a well-developed, powerful and comparative computational approach, EST-based homology search is applied to find potential miRNA of coffee. We blasted publicly available EST sequences obtained from NCBI GenBank against previously known plant miRNAs. For the first time, one potential miRNA from a large miRNA family with appropriate fold back structures was identified through a series of filtration criteria. A total of six potential target genes in Arabidopsis were identified based on their sequence complementarities. The target genes mainly encode transport inhibitor like protein, transcription factor, DNA-binding protein, and GRR1-like protein, and these genes play an important role in various biological processes like response to chitin, cold, salt stress, water deprivation etc. Overall, findings from this study will accelerate the way for further researches of miRNAs and their functions in coffee.

Identification of new stress-induced microRNA and their targets in wheat using computational approach

MicroRNAs (miRNAs) are a class of short endogenous non-coding small RNA molecules of about 18-22 nucleotides in length. Their main function is to downregulate gene expression in different manners like translational repression, mRNA cleavage and epigenetic modification. Computational predictions have raised the number of miRNAs in wheat significantly using an EST based approach. Hence, a combinatorial approach which is amalgamation of bioinformatics software and perl script was used to identify new miRNA to add to the growing database of wheat miRNA. Identification of miRNAs was initiated by mining the EST (Expressed Sequence Tags) database available at National Center for Biotechnology Information. In this investigation, 4677 mature microRNA sequences belonging to 50 miRNA families from different plant species were used to predict miRNA in wheat. A total of five abiotic stress-responsive new miRNAs were predicted and named Ta-miR5653, Ta-miR855, Ta-miR819k, Ta-miR3708 and Ta-miR5156. In addition, four previously identified miRNA, i.e., Ta-miR1122, miR1117, Ta-miR1134 and Ta-miR1133 were predicted in newly identified EST sequence and 14 potential target genes were subsequently predicted, most of which seems to encode ubiquitin carrier protein, serine/threonine protein kinase, 40S ribosomal protein, F-box/kelch-repeat protein, BTB/POZ domain-containing protein, transcription factors which are involved in growth, development, metabolism and stress response. Our result has increased the number of miRNAs in wheat, which should be useful for further investigation into the biological functions and evolution of miRNAs in wheat and other plant species.

Identification and characterization of microRNAs in Asiatic cotton (Gossypium arboreum L.)

To date, no miRNAs have been identified in the important diploid cotton species although there are several reports on miRNAs in upland cotton. In this study, we identified 73 miRNAs, belonging to 49 families, from Asiatic cotton using a well-developed comparative genome-based homologue search. Several of the predicted miRNAs were validated using quantitative real time PCR (qRT-PCR). The length of miRNAs varied from 18 to 22 nt with an average of 20 nt. The length of miRNA precursors also varied from 46 to 684 nt with an average of 138 ±120 nt. For a majority of Asiatic cotton miRNAs, there is only one member per family; however, multiple members were identified for miRNA 156, 414, 837, 838, 1044, 1533, 2902, 2868, 5021 and 5142 families. Nucleotides A and U were dominant, accounted for 62.95%, in the Asiatic cotton pre-miRNAs. The Asiatic cotton pre-miRNAs had high negative minimal folding free energy (MFE) and adjusted MFE (AMFE) and high MFE index (MFEI). Many miRNAs identified in Asiatic cotton suggest that miRNAs also play a similar regulatory mechanism in diploid cotton.

First survey of the wheat chromosome 5A composition through a next generation sequencing approach

Wheat is one of the world's most important crops and is characterized by a large polyploid genome. One way to reduce genome complexity is to isolate single chromosomes using flow cytometry. Low coverage DNA sequencing can provide a snapshot of individual chromosomes, allowing a fast characterization of their main features and comparison with other genomes. We used massively parallel 454 pyrosequencing to obtain a 2x coverage of wheat chromosome 5A. The resulting sequence assembly was used to identify TEs, genes and miRNAs, as well as to infer a virtual gene order based on the synteny with other grass genomes. Repetitive elements account for more than 75% of the genome. Gene content was estimated considering non-redundant reads showing at least one match to ESTs or proteins. The results indicate that the coding fraction represents 1.08% and 1.3% of the short and long arm respectively, projecting the number of genes of the whole chromosome to approximately 5,000. 195 candidate miRNA precursors belonging to 16 miRNA families were identified. The 5A genes were used to search for syntenic relationships between grass genomes. The short arm is closely related to Brachypodium chromosome 4, sorghum chromosome 8 and rice chromosome 12; the long arm to regions of Brachypodium chromosomes 4 and 1, sorghum chromosomes 1 and 2 and rice chromosomes 9 and 3. From these similarities it was possible to infer the virtual gene order of 392 (5AS) and 1,480 (5AL) genes of chromosome 5A, which was compared to, and found to be largely congruent with the available physical map of this chromosome.

Computational Identification of Novel MicroRNAs and Their Targets in Vigna unguiculata

MicroRNAs (miRNAs) are a class of endogenous, noncoding, short RNAs directly involved in regulating gene expression at the posttranscriptional level. High conservation of miRNAs in plant provides the foundation for identification of new miRNAs in other plant species through homology alignment. Here, previous known plant miRNAs were BLASTed against the Expressed Sequence Tag (EST) and Genomic Survey Sequence (GSS) databases of Vigna unguiculata, and according to a series of filtering criteria, a total of 47 miRNAs belonging to 13 miRNA families were identified, and 30 potential target genes of them were subsequently predicted, most of which seemed to encode transcription factors or enzymes participating in regulation of development, growth, metabolism, and other physiological processes. Overall, our findings lay the foundation for further researches of miRNAs function in Vigna unguiculata.

Identification and characterization of conserved microRNAs and their target genes in wheat (Triticum aestivum)

MicroRNAs (miRNAs) are non-coding small RNAs that regulate gene expression by translational repression or transcript degradation. A large number of miRNAs have been identified from model plant species; however, the character of conserved miRNAs is poorly understood. We studied 42 miRNA families that are conserved within the plant kingdom, using the miRBase database. Some conserved miRNA families were found to be preferentially expressed in dicots relative to monocots, especially miR403, miR472 and miR479. Using an improved homology search-based approach and the conserved miRNAs as the query set, 34 conserved miRNAs and the miR482 family were identified in wheat. Forty-six wheat mRNAs were predicted as their putative target genes. Most conserved wheat miRNAs were found to retain homologous target interactions and have analogous molecular functions. The miR172 displayed a wheat-specific function and was found to have an additional target interaction with succinyl-CoA ligase. We concluded that although miRNAs are conserved, the expression and function of some have drifted during long periods of plant evolution.