Identification of microRNAs in wild soybean (Glycine soja)

The oxidative stress 3-like protein GsOS3L, substrate of GsSnRK1, enhances salt and cadmium stresses in soybean roots

Salt and heavy metal stresses have significant impacts on crop growth and agricultural development. Wild soybean (Glycine soja) exhibits greater resistance to abiotic stresses than its cultivated counterpart (Glycine max). In this study, the oxidative stress 3-like gene GsOS3L was identified from yeast two hybridization cDNA library constructed from wild soybean RNA. 150 mM NaCl and 10 % PEG induced their expression in roots, stems and leaves, respectively. pGsOS3L:GUS activity was enhanced in various tissues with increasing NaCl and CdCl2 concentrations. Y2H, BiFC, and LCI activity assays revealed that GsOS3L interacted physically with the GsSnRK1 kinase. The GsOS3L protein, which contains potential phosphorylation and palmitoylation sites, was localized to the nucleus under normal conditions but translocated from the nucleus to the cell membrane under cadmium stress. This translocation was prevented by the palmitoylation inhibitor 2-bromopalmitate (2-BP) and by double mutation of the predicted palmitoylation sites (C4S/C21S). 2-BP treatment attenuated GsOS3L transgenic composite soybeans' cadmium stress tolerance. GsOS3L was found to be phosphorylated by GsSnRK1, which reduced the salt and cadmium stress tolerance of transgenic Arabidopsis plants. The findings of this study provide promising insights into the physiological and molecular mechanisms of GsOS3L in soybean under salt and cadmium stresses.

Analyzing the defense response mechanism of Atractylodes macrocephala to Fusarium oxysporum through small RNA and degradome sequencing

Introduction

Fusarium oxysporum is a significant soil-borne fungal pathogen that affects over 100 plant species, including crucial crops like tomatoes, bananas, cotton, cucumbers, and watermelons, leading to wilting, yellowing, growth inhibition, and ultimately plant death. The root rot disease of A. macrocephala, caused by F. oxysporum, is one of the most serious diseases in continuous cropping, which seriously affects its sustainable development.

Methods

In this study, we explored the interaction between A. macrocephala and F. oxysporum through integrated small RNA (sRNA) and degradome sequencing to uncover the microRNA (miRNA)-mediated defense mechanisms.

Results

We identified colonization of F. oxysporum in A. macrocephala roots on day 6. Nine sRNA samples were sequenced to examine the dynamic changes in miRNA expression in A. macrocephala infected by F. oxysporum at 0, 6, and 12 days after inoculation. Furthermore, we using degradome sequencing and quantitative real-time PCR (qRT-PCR), validated four miRNA/target regulatory units involved in A. macrocephala-F. oxysporum interactions.

Discussion

This study provides new insights into the molecular mechanisms underlying A. macrocephala's early defense against F. oxysporum infection, suggesting directions for enhancing resistance against this pathogen.

Using the Knowledge of Post-transcriptional Regulations to Guide Gene Selections for Molecular Breeding in Soybean

The omics approaches allow the scientific community to successfully identify genomic regions associated with traits of interest for marker-assisted breeding. Agronomic traits such as seed color, yield, growth habit, and stress tolerance have been the targets for soybean molecular breeding. Genes governing these traits often undergo post-transcriptional modifications, which should be taken into consideration when choosing elite genes for molecular breeding. Post-transcriptional regulations of genes include transcript regulations, protein modifications, and even the regulation of the translational machinery. Transcript regulations involve elements such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) for the maintenance of transcript stability or regulation of translation efficiency. Protein modifications involve molecular modifications of target proteins and the alterations of their interacting partners. Regulations of the translational machinery include those on translation factors and the ribosomal protein complex. Post-transcriptional regulations usually involve a set of genes instead of a single gene. Such a property may facilitate molecular breeding. In this review, we will discuss the post-transcriptional modifications of genes related to favorable agronomic traits such as stress tolerance, growth, and nutrient uptake, using examples from soybean as well as other crops. The examples from other crops may guide the selection of genes for marker-assisted breeding in soybean.

Differences in the metabolic profiles and antioxidant activities of wild and cultivated black soybeans evaluated by correlation analysis

Wild soybeans are considered a potential resource for soybean domestication and an important source of genetic diversity for soybean crop improvement. Understanding metabolite-caused bioactivity differences between cultivated and wild soybeans is essential for designing a soybean with enhanced nutritional traits. In this study, the non-targeted metabolic profiling of 26 soybean varieties, 15 wild black soybeans (WBS) and 11 cultivated black soybeans (CBS), using liquid chromatography-mass spectrometry (LC-MS) in combination with multivariate analysis revealed significant differences in 25 differential metabolites. Among these, the soyasaponins Ab and Bb were found to be characteristic metabolites expressed more substantially in CBS than in WBS. Three different antioxidant assays and correlation analysis identified major and minor antioxidants that contributed to WBS having an antioxidant activity 4- to 8-fold stronger than that of CBS. Epicatechin, procyanidin B2, and cyanidin-3-O-glucoside were identified by both association analysis and the online LC-ABTS radical scavenging assay as being major antioxidants.

Extrapolative microRNA precursor based SSR mining from tea EST database in respect to agronomic traits

Tea (Camellia sinensis, (L.) Kuntze) is considered as most popular drink across the world and it is widely consumed beverage for its several health-benefit characteristics. These positive traits primarily rely on its regulatory networks of different metabolic pathways. Development of microsatellite markers from the conserved genomic regions are being worthwhile for reviewing the genetic diversity of closely related species or self-pollinated species. Although several SSR markers have been reported, in tea, the trait-specific Simple Sequence Repeat (SSR) markers, leading to be useful in marker assisted breeding technique, are yet to be identified. Micro RNAs are short, non-coding RNA molecules, involved in post transcriptional mode of gene regulation and thus effects on related phenotype. Present study deals with identification of the microsatellite motifs within the reported and predicted miRNA precursors that are effectively followed by designing of primers from SSR flanking regions in order to PCR validation. In addition to the earlier reports, two new miRNAs are predicting here from tea expressed tag sequence database. Furthermore, 18 SSR motifs are found to be in 13 of all 33 predicted miRNAs. Trinucleotide motifs are most abundant among all followed by dinucleotides. Since, miRNA based SSR markers are evidenced to have significant role on genetic fingerprinting study, these outcomes would pave the way in developing novel markers for tagging tea specific agronomic traits as well as substantiating non-conventional breeding program.

An improved method to quantitate mature plant microRNA in biological matrices using modified periodate treatment and inclusion of internal controls

MicroRNAs (miRNAs) ubiquitously exist in microorganisms, plants, and animals, and appear to modulate a wide range of critical biological processes. However, no definitive conclusion has been reached regarding the uptake of exogenous dietary small RNAs into mammalian circulation and organs and cross-kingdom regulation. One of the critical issues is our ability to assess and distinguish the origin of miRNAs. Although periodate oxidation has been used to differentiate mammalian and plant miRNAs, validation of treatment efficiency and the inclusion of proper controls for this method were lacking in previous studies. This study aimed to address: 1) the efficiency of periodate treatment in a plant or mammalian RNA matrix, and 2) the necessity of inclusion of internal controls. We designed and tested spike-in synthetic miRNAs in various plant and mammalian matrices and showed that they can be used as a control for the completion of periodate oxidation. We found that overloading the reaction system with high concentration of RNA resulted in incomplete oxidation of unmethylated miRNA. The abundant miRNAs from soy and corn were analyzed in the plasma, liver, and fecal samples of C57BL/6 mice fed a corn and soy-based chow diet using our improved methodology. The improvement resulted in the elimination of the false positive detection in the liver, and we did not detect plant miRNAs in the mouse plasma or liver samples. In summary, an improved methodology was developed for plant miRNA detection that appears to work well in different sample matrices.

Genome-wide identification and characterization of Eutrema salsugineum microRNAs for salt tolerance

Eutrema salsugineum, a close relative of Arabidopsis thaliana, is a valuable halophytic model plant that has extreme tolerance to salinity. As posttranscriptional gene regulators, microRNAs (miRNAs) control gene expression and a variety of biological processes, including plant-stress responses. To identify salt-stress responsive miRNAs in E. salsugineum and reveal their possible roles in the adaptive response to salt stress, we chose the Solexa sequencing platform to screen the miRNAs in 4-week-old E. salsugineum seedlings under salt treatment. A total of 82 conserved miRNAs belonging to 27 miRNA families and 17 novel miRNAs were identified and 11 conserved miRNA families and 4 novel miRNAs showed a significant response to salt stress. To investigate the possible biological roles of miRNAs, 1060 potential targets were predicted. Moreover, 35 gene ontology (GO) categories and 1 pathway, including a few terms that were directly and indirectly related to salt stress, were significantly enriched in the salt-stress-responsive miRNAs targets. The relative expression analysis of six target genes was analyzed using quantitative real-time polymerase chain reaction (PCR) and showed a negative correlation with their corresponding miRNAs. Many stress regulatory and phytohormone regulatory cis-regulatory elements were widely present in the promoter region of the salt-responsive miRNA precursors. This study describes the large-scale characterization of E. salsugineum miRNAs and provides a useful resource for further understanding of miRNA functions in the regulation of the E. salsugineum salt-stress response.

Comparative Profiling of microRNA Expression in Soybean Seeds from Genetically Modified Plants and their Near-Isogenic Parental Lines

MicroRNAs (miRNAs) have been widely demonstrated to play fundamental roles in gene regulation in most eukaryotes. To date, there has been no study describing the miRNA composition in genetically modified organisms (GMOs). In this study, small RNAs from dry seeds of two GM soybean lines and their parental cultivars were investigated using deep sequencing technology and bioinformatic approaches. As a result, several differentially expressed gma-miRNAs were found between the GM and non-GM soybeans. Meanwhile, more differentially expressed gma-miRNAs were identified between distantly relatednon-GM soybeans, indicating that the miRNA components of soybean seeds varied among different soybean lines, including the GM and non-GM soybeans, and the extent of difference might be related to their genetic relationship. Additionally, fourteen novel gma-miRNA candidates were predicted in soybean seeds including a potential bidirectionally transcribed miRNA family with two genomic loci (gma-miR-N1). Our findings firstly provided useful data for miRNA composition in edible GM crops and also provided valuable information for soybean miRNA research.

Identification of Novel and Conserved miRNAs in Leaves of In vitro Grown Citrus reticulata "Lugan" Plantlets by Solexa Sequencing

MicroRNAs (miRNAs) play essential roles in plant development, but the roles in the in vitro plant development are unknown. Leaves of ponkan plantlets derived from mature embryos at in vitro culture conditions were used to sequence small RNA fraction via Solexa sequencing, and the miRNAs expression was analyzed. The results showed that there were 3,065,625 unique sequences in ponkan, of which 0.79% were miRNAs. The RNA sequences with lengths of 18-25 nt derived from the library were analyzed, leading to the identification of 224 known miRNAs, of which the most abundant were miR157, miR156, and miR166. Three hundred and fifty-eight novel miRNA candidates were also identified, and the number of reads of ponkan novel miRNAs varied from 5 to 168,273. The expression of the most known miRNAs obtained was at low levels, which varied from 5 to 4,946,356. To better understand the role of miRNAs during the preservation of ponkan in vitro plantlet, the expression patterns of cre-miR156a/159b/160a/166a/167a/168a/171/398b were validated by quantitative real-time PCR (qPCR). The results showed that not only the development-associated miRNAs, e.g., cre-miR156/159/166/396, expressed highly at the early preservation period in the in vitro ponkan plantlet leaves but also the stress-related miRNAs, e.g., cre-miR171 and cre-miR398b, expressed highly at the same time. The expression levels of most tested miRNAs were found to decrease after 6 months and the amounts of these miRNAs were kept at low levels at 18 months. After analyzing the expression level of their targets during the reservation of the ponkan in vitro plantlet, development-associated cre-ARF6 and stress-related cre-CSD modules exhibited negative correlation with miR167 and miR398, respectively, indicating an involvement of the miRNAs in the in vitro development of ponkan and function in the conservation of ponkan germplasm.

Identification of miRNAs and their targets by high-throughput sequencing and degradome analysis in cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B of soybean

BACKGROUND

Cytoplasmic male sterility (CMS) provides crucial breeding materials that facilitate hybrid seed production in various crops, and thus plays an important role in the study of hybrid vigor (heterosis), in plants. However, the CMS regulatory network in soybean remains unclear. MicroRNAs (miRNAs) play crucial roles in flower and pollen development by targeting genes that regulate their expression in plants. To identify the miRNAs and their targets that exist in the soybean CMS line NJCMS1A and its maintainer NJCMS1B, high-throughput sequencing and degradome analysis were conducted in this study.

RESULTS

Two small RNA libraries were constructed from the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. A total of 105 new miRNAs present on the other arm of known pre-miRNAs, 23 new miRNA members, 158 novel miRNAs and 160 high-confidence soybean miRNAs were identified using high-throughput sequencing. Among the identified miRNAs, 101 differentially expressed miRNAs with greater than two-fold changes between NJCMS1A and NJCMS1B were discovered. The different expression levels of selected miRNAs were confirmed by stem-loop quantitative real-time PCR. A degradome analysis showed that 856 targets were predicted to be targeted by 296 miRNAs, including a squamosa promoter-binding protein-like transcription factor family protein, a pentatricopeptide repeat-containing protein, and an auxin response factor, which were previously shown to be involved in floral organ or anther development in plants. Additionally, some targets, including a MADS-box transcription factor, NADP-dependent isocitrate dehydrogenase and NADH-ubiquinone oxidoreductase 24 kDa subunit, were identified, and they may have some relationship with the programmed cell death, reactive oxygen species accumulation and energy deficiencies, which might lead to soybean male sterility.

CONCLUSIONS

The present study is the first to use deep sequencing technology to identify miRNAs and their targets in the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. The results revealed that the miRNAs might participate in flower and pollen development, which could facilitate our understanding of the molecular mechanisms behind CMS in soybean.

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.

Identification of MicroRNAs in Response to Different Day Lengths in Soybean Using High-Throughput Sequencing and qRT-PCR

MicroRNAs (miRNAs) are short, non-coding single-strand RNA molecules that play important roles in plant growth, development and stress responses. Flowering time affects the seed yield and quality of soybean. However, the miRNAs involved in the regulation of flowering time in soybean have not been reported until recently. Here, high-throughput sequencing and qRT-PCR were used to identify miRNAs involved in soybean photoperiodic pathways. The first trifoliate leaves of soybean that receive the signal of light treatment were used to construct six libraries (0, 8, and 16 h under short-day (SD) treatment and 0, 8, and 16 h under long-day (LD) treatment). The libraries were sequenced using Illumina Solexa. A total of 318 known plant miRNAs belonging to 163 miRNA families and 81 novel predicted miRNAs were identified. Among these, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h, showed differences in abundance between LD and SD treatments. Furthermore, the results of GO and KEGG analyses indicated that most of the miRNA targets were transcription factors. Seven miRNAs at 0 h, 23 miRNAs (including four novel predicted miRNAs) at 8 h, 16 miRNAs (including one novel predicted miRNA) at 16 h and miRNA targets were selected for qRT-PCR analysis to assess the accuracy of the sequencing and target prediction. The results indicated that the expression patterns of the selected miRNAs and miRNA targets showed no differences between the qRT-PCR and sequencing results. In addition, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h responded to day length changes in soybean, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h. These results provided an important molecular basis to understand the regulation of flowering time through photoperiodic pathways in soybean.

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.

Exploring microRNA-like small RNAs in the filamentous fungus Fusarium oxysporum

RNA silencing such as quelling and meiotic silencing by unpaired DNA (MSUD) and several other classes of special small RNAs have been discovered in filamentous fungi recently. More than four different mechanisms of microRNA-like RNAs (milRNAs) production have been illustrated in the model fungus Neurospora crassa including a dicer-independent pathway. To date, very little work focusing on small RNAs in fungi has been reported and no universal or particular characteristic of milRNAs were defined clearly. In this study, small RNA and degradome libraries were constructed and subsequently deep sequenced for investigating milRNAs and their potential cleavage targets on the genome level in the filamentous fungus F. oxysporum f. sp. lycopersici. As a result, there is no intersection of conserved miRNAs found by BLASTing against the miRBase. Further analysis showed that the small RNA population of F. oxysporum shared many common features with the small RNAs from N. crassa and other fungi. According to the known standards of miRNA prediction in plants and animals, milRNA candidates from 8 families (comprising 19 members) were screened out and identified. However, none of them could trigger target cleavage based on the degradome data. Moreover, most major signals of cleavage in transcripts could not match appropriate complementary small RNAs, suggesting that other predominant modes for milRNA-mediated gene regulation could exist in F. oxysporum. In addition, the PAREsnip program was utilized for comprehensive analysis and 3 families of small RNAs leading to transcript cleavage were experimentally validated. Altogether, our findings provided valuable information and important hints for better understanding the functions of the small RNAs and milRNAs in the fungal kingdom.

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.

Analysis of promoters of microRNAs from a Glycine max degradome library

OBJECTIVE

MicroRNAs (miRNAs) are genome-encoded, small non-coding RNAs that play important functions in development, biotic and abiotic stress responses, and other processes. Our aim was to explore the regulation of miRNA expression.

METHODS

We used bioinformatics methods to predict the core promoters of 440 miRNAs identified from a soybean (Glycine max) degradome library and to analyze cis-acting elements for 369 miRNAs.

RESULTS

The prediction results showed that 83.86% of the 440 miRNAs contained promoters in their upstream sequences, and 8.64% (38 loci) in their downstream sequences. The distributions of two core promoter elements, TATA-boxes and transcription start sites (TSSs), were similar. The cis-acting elements were examined to provide clues to the function and regulation of spatiotemporal expression of the miRNAs. Analyses of miRNA cis-elements and targets indicated a potential auxin response factor (ARF)- and gibberellin response factor (GARF)-mediated negative feedback loop for miRNA expression.

CONCLUSIONS

The features of miRNAs from a Glycine max degradome library obtained here provide insights into the transcription regulation and functions of miRNAs in soybean.

Deep sequencing discovery of novel and conserved microRNAs in wild type and a white-flesh mutant strawberry

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression by base pairing to mRNA target sequences, and play crucial roles in plant development and stress responses. The knowledge on post-transcriptional regulation by miRNAs in strawberry is rather limited so far. In order to understand the role of miRNA in the molecular control during strawberry fruit development, small RNA libraries were constructed from fruits at the turning stage of strawberry cultivar 'Sachinoka' and its white-flesh mutant by using the Solexa platform. One hundred and twenty conserved miRNAs belonging to 27 miRNA families and 33 putative novel strawberry miRNAs were identified in both libraries. Their target genes were predicted using the Fragaria vesca genome. Nine of all miRNAs showed significant expression differences between two types of samples. Four miRNAs were up-regulated and five were down-regulated in white-flesh mutant. The sequencing results were partially validated by quantitative RT-PCR. Among them, the expression of miR399a shows the biggest change between the two samples. The prediction of its target gene showed that miR399 may play an important role in phosphate homeostasis of strawberry fruits. Furthermore, we deduce that the expression of miR399 has negative correlation with the content of sugars.

Identification and comparative analysis of low phosphate tolerance-associated microRNAs in two maize genotypes

Low phosphate (Pi) availability is a major constraint on maize growth and productivity. MicroRNAs (miRNAs) are known to play critical roles in plant responses to various environmental conditions. The identification of low Pi tolerance-associated miRNAs will accelerate the development of Pi starvation tolerant maize plants. However, miRNAs associated with low Pi tolerance have not been identified. In this study, we compared deep sequencing small RNA reads from two maize genotypes, the wild type, Qi319, and the low Pi tolerant mutant, 99038, under normal and low Pi conditions. Six known miRNA families and seven novel miRNAs were found differently expressed by the two genotypes. All these miRNAs were confirmed by sequencing a second batch of small RNA libraries constructed in the same way as those used in the first sequencing. The expression profiles of some of these miRNAs were further confirmed by real-time PCR. The predicted target genes of the low Pi tolerance-associated miRNAs were involved in root development or stress responses. Expression levels of some of target genes were significantly different between Qi319 and 99038. These findings suggested that miRNAs may play important roles in low Pi tolerance in maize and may be a key factor in determining the level of low Pi tolerance in different maize genotypes. This study provides an approach for identifying low Pi tolerance-associated miRNAs and can help in the selection and manipulation of high performing maize genotypes under low Pi conditions.

Elucidation of miRNAs-mediated responses to low nitrogen stress by deep sequencing of two soybean genotypes

Nitrogen (N) is a major limiting factor in crop production, and plant adaptive responses to low N are involved in many post-transcriptional regulation. Recent studies indicate that miRNAs play important roles in adaptive responses. However, miRNAs in soybean adaptive responses to N limitation have been not reported. We constructed sixteen libraries to identify low N-responsive miRNAs on a genome-wide scale using samples from 2 different genotypes (low N sensitive and low N tolerant) subjected to various periods of low nitrogen stress. Using high-throughput sequencing technology (Illumina-Solexa), we identified 362 known miRNAs variants belonging to 158 families and 90 new miRNAs belonging to 55 families. Among these known miRNAs variants, almost 50% were not different from annotated miRNAs in miRBase. Analyses of their expression patterns showed 150 known miRNAs variants as well as 2 novel miRNAs with differential expressions. These differentially expressed miRNAs between the two soybean genotypes were compared and classified into three groups based on their expression patterns. Predicted targets of these miRNAs were involved in various metabolic and regulatory pathways such as protein degradation, carbohydrate metabolism, hormone signaling pathway, and cellular transport. These findings suggest that miRNAs play important roles in soybean response to low N and contribute to the understanding of the genetic basis of differences in adaptive responses to N limitation between the two soybean genotypes. Our study provides basis for expounding the complex gene regulatory network of these miRNAs.

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 of wild soybean miRNAs and their target genes responsive to aluminum stress

BACKGROUND

MicroRNAs (miRNAs) play important regulatory roles in development and stress response in plants. Wild soybean (Glycine soja) has undergone long-term natural selection and may have evolved special mechanisms to survive stress conditions as a result. However, little information about miRNAs especially miRNAs responsive to aluminum (Al) stress is available in wild soybean.

RESULTS

Two small RNA libraries and two degradome libraries were constructed from the roots of Al-treated and Al-free G. soja seedlings. For miRNA identification, a total of 7,287,655 and 7,035,914 clean reads in Al-treated and Al-free small RNAs libraries, respectively, were generated, and 97 known miRNAs and 31 novel miRNAs were identified. In addition, 49 p3 or p5 strands of known miRNAs were found. Among all the identified miRNAs, the expressions of 30 miRNAs were responsive to Al stress. Through degradome sequencing, 86 genes were identified as targets of the known miRNAs and five genes were found to be the targets of the novel miRNAs obtained in this study. Gene ontology (GO) annotations of target transcripts indicated that 52 target genes cleaved by conserved miRNA families might play roles in the regulation of transcription. Additionally, some genes, such as those for the auxin response factor (ARF), domain-containing disease resistance protein (NB-ARC), leucine-rich repeat and toll/interleukin-1 receptor-like protein (LRR-TIR) domain protein, cation transporting ATPase, Myb transcription factors, and the no apical meristem (NAM) protein, that are known to be responsive to stress, were found to be cleaved under Al stress conditions.

CONCLUSIONS

A number of miRNAs and their targets were detected in wild soybean. Some of them that were responsive to biotic and abiotic stresses were regulated by Al stress. These findings provide valuable information to understand the function of miRNAs in Al tolerance.

High-throughput sequencing and characterization of the small RNA transcriptome reveal features of novel and conserved microRNAs in Panax ginseng

microRNAs (miRNAs) play vital regulatory roles in many organisms through direct cleavage of transcripts, translational repression, or chromatin modification. Identification of miRNAs has been carried out in various plant species. However, no information is available for miRNAs from Panax ginseng, an economically significant medicinal plant species. Using the next generation high-throughput sequencing technology, we obtained 13,326,328 small RNA reads from the roots, stems, leaves and flowers of P. ginseng. Analysis of these small RNAs revealed the existence of a large, diverse and highly complicated small RNA population in P. ginseng. We identified 73 conserved miRNAs, which could be grouped into 33 families, and 28 non-conserved ones belonging to 9 families. Characterization of P. ginseng miRNA precursors revealed many features, such as production of two miRNAs from distinct regions of a precursor, clusters of two precursors in a transcript, and generation of miRNAs from both sense and antisense transcripts. It suggests the complexity of miRNA production in P. gingseng. Using a computational approach, we predicted for the conserved and non-conserved miRNA families 99 and 31 target genes, respectively, of which eight were experimentally validated. Among all predicted targets, only about 20% are conserved among various plant species, whereas the others appear to be non-conserved, indicating the diversity of miRNA functions. Consistently, many miRNAs exhibited tissue-specific expression patterns. Moreover, we identified five dehydration- and ten heat-responsive miRNAs and found the existence of a crosstalk among some of the stress-responsive miRNAs. Our results provide the first clue to the elucidation of miRNA functions in P. ginseng.

Identification of soybean seed developmental stage-specific and tissue-specific miRNA targets by degradome sequencing

BACKGROUND

MicroRNAs (miRNAs) regulate the expression of target genes by mediating gene silencing in both plants and animals. The miRNA targets have been extensively investigated in Arabidopsis and rice using computational prediction, experimental validation by overexpression in transgenic plants, and by degradome or PARE (parallel analysis of RNA ends) sequencing. However, miRNA targets mostly remain unknown in soybean (Glycine max). More specifically miRNA mediated gene regulation at different seed developmental stages in soybean is largely unexplored. In order to dissect miRNA guided gene regulation in soybean developing seeds, we performed a transcriptome-wide experimental method using degradome sequencing to directly detect cleaved miRNA targets.

RESULTS

In this study, degradome libraries were separately prepared from immature soybean cotyledons representing three stages of development and from seed coats of two stages. Sequencing and analysis of 10 to 40 million reads from each library resulted in identification of 183 different targets for 53 known soybean miRNAs. Among these, some were found only in the cotyledons representing cleavage by 25 miRNAs and others were found only in the seed coats reflecting cleavage by 12 miRNAs. A large number of targets for 16 miRNAs families were identified in both tissues irrespective of the stage. Interestingly, we identified more miRNA targets in the desiccating cotyledons of late seed maturation than in immature seed. We validated four different auxin response factor genes as targets for gma-miR160 via RNA ligase mediated 5' rapid amplification of cDNA ends (RLM-5'RACE). Gene Ontology (GO) analysis indicated the involvement of miRNA target genes in various cellular processes during seed development.

CONCLUSIONS

The miRNA targets in both the cotyledons and seed coats of several stages of soybean seed development have been elucidated by experimental evidence from comprehensive, high throughput sequencing of the enriched fragments resulting from miRNA-guided cleavage of messenger RNAs. Nearly 50% of the miRNA targets were transcription factors in pathways that are likely important in setting or maintaining the developmental program leading to high quality soybean seeds that are one of the dominant sources of protein and oil in world markets.

Identification of soybean microRNAs involved in soybean cyst nematode infection by deep sequencing

Soybean cyst nematode (SCN), Heterodera glycines, is the most devastating pathogen of soybean worldwide. MicroRNAs (miRNAs) are a class of small, non-coding RNAs that are known to play important role in plant stress response. However, there are few reports profiling the miRNA expression patterns during pathogen stress. We sequenced four small RNA libraries from two soybean cultivar (Hairbin xiaoheidou, SCN race 3 resistant, Liaodou 10, SCN race 3 susceptible) that grown under un-inoculated and SCN-inoculated soil. Small RNAs were mapped to soybean genome sequence, 364 known soybean miRNA genes were identified in total. In addition, 21 potential miRNA candidates were identified. Comparative analysis of miRNA profiling indicated 101 miRNAs belong to 40 families were SCN-responsive. We also found 20 miRNAs with different express pattern even between two cultivars of the same species. These findings suggest that miRNA paly important role in soybean response to SCN and have important implications for further identification of miRNAs under pathogen stress.

Genome-wide profiling of novel and conserved Populus microRNAs involved in pathogen stress response by deep sequencing

MicroRNAs (miRNAs) are small RNAs, generally of 20-23 nt, that down-regulate target gene expression during development, differentiation, growth, and metabolism. In Populus, extensive studies of miRNAs involved in cold, heat, dehydration, salinity, and mechanical stresses have been performed; however, there are few reports profiling the miRNA expression patterns during pathogen stress. We obtained almost 38 million raw reads through Solexa sequencing of two libraries from Populus inoculated and uninoculated with canker disease pathogen. Sequence analyses identified 74 conserved miRNA sequences belonging to 37 miRNA families from 154 loci in the Populus genome and 27 novel miRNA sequences from 35 loci, including their complementary miRNA* strands. Intriguingly, the miRNA* of three conserved miRNAs were more abundant than their corresponding miRNAs. The overall expression levels of conserved miRNAs increased when subjected to pathogen stress, and expression levels of 33 miRNA sequences markedly changed. The expression trends determined by sequencing and by qRT-PCR were similar. Finally, nine target genes for three conserved miRNAs and 63 target genes for novel miRNAs were predicted using computational analysis, and their functions were annotated. Deep sequencing provides an opportunity to identify pathogen-regulated miRNAs in trees, which will help in understanding the regulatory mechanisms of plant defense responses during pathogen infection.

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.

Identification of miRNAs in sorghum by using bioinformatics approach

MicroRNAs (miRNAs) regulate gene expression mainly by post-transcriptional gene silencing (PTGS) and in some cases by transcriptional genes silencing (TGS). miRNAs play critical roles in developmental processes, nutrient homeostasis, abiotic stress and pathogen responses of plants. In contrast to the large number of miRNAs predicted in cereal model plant rice, only 148 miRNAs were predicted in sorghum till date (miRBase release 17). This suggested that miRNAs identified in sorghum is far from saturation. Hence, we developed a bioinformatics pipeline using an in-house PERL script and publicly available structure prediction tools to identify miRNAs and their target genes from publically available Expressed Sequence Tags (EST) and Genomic Survey Sequence (GSS). About 1379 known and unique plant miRNAs from 33 different crops were used to predict new miRNAs in sorghum. We identified 31 new miRNAs belonging to 10 different miRNA families. We predicted 72 potential target genes for 31 miRNAs, and most of these target genes are predicted to be involved in plant growth and development.These newly identified miRNAs add to the growing database of miRNA and lay the foundation for further understanding of miRNA function in sorghum plant development.

Identification of MicroRNAs and their targets in Helianthus

MicroRNAs (miRNAs) are small non-coding RNAs about 20-24 nucleotides long, playing regulatory role. The conserved nature among the various organisms makes them a good source of new miRNAs discovery by comparative genomics approach using bioinformatic tools. A systematic search approach was used for inter-species homologs of miRNA precursors (pre-miRNAs), from known Helianthus expressed sequence tags (ESTs). The study resulted in 61 novel miRNAs belonging to 34 families from Helianthus ESTs. The 28 miRNA families; mir 159,160, 164, 170, 390, 393, 413, 415, 419, 426, 446, 530, 822, 842, 846, 1310, 1888, 2086, 2657, 2667, 2678, 2659, 2911, 2938, 3440, 3521, 3623, and 3630 are reporting for the first time in Helianthus. In the 61 new miRNAs, 20 are from H. tuberosus, 17 miRNAs belong to H. annus, 8 are from H. ciliaris, 5 to H. exilis, 4 is from H. argophyllous, H. petiolaris each and 3 are from H. paradoxus. All the pre-miRNAs form stable minimum free energy (mfe) stem-loop structure as their orthologues form and the mature miRNAs reside in the stem portion of the stem-loop structures. Their targets consist of growth and development related, transcription factors, signalling pathway kinases, stress resistant proteins and transport related proteins.

Identification of novel soybean microRNAs involved in abiotic and biotic stresses

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Identification of miRNAs and their target genes in developing soybean seeds by deep sequencing

BACKGROUND

MicroRNAs (miRNAs) regulate gene expression by mediating gene silencing at transcriptional and post-transcriptional levels in higher plants. miRNAs and related target genes have been widely studied in model plants such as Arabidopsis and rice; however, the number of identified miRNAs in soybean (Glycine max) is limited, and global identification of the related miRNA targets has not been reported in previous research.

RESULTS

In our study, a small RNA library and a degradome library were constructed from developing soybean seeds for deep sequencing. We identified 26 new miRNAs in soybean by bioinformatic analysis and further confirmed their expression by stem-loop RT-PCR. The miRNA star sequences of 38 known miRNAs and 8 new miRNAs were also discovered, providing additional evidence for the existence of miRNAs. Through degradome sequencing, 145 and 25 genes were identified as targets of annotated miRNAs and new miRNAs, respectively. GO analysis indicated that many of the identified miRNA targets may function in soybean seed development. Additionally, a soybean homolog of Arabidopsis suppressor of gene silencing 3 (AtSGS3) was detected as a target of the newly identified miRNA Soy_25, suggesting the presence of feedback control of miRNA biogenesis.

CONCLUSIONS

We have identified large numbers of miRNAs and their related target genes through deep sequencing of a small RNA library and a degradome library. Our study provides more information about the regulatory network of miRNAs in soybean and advances our understanding of miRNA functions during seed development.