Identification of novel and candidate miRNAs in rice by high throughput sequencing

Characterization of miRNAs responsive to exogenous ethylene in grapevine berries at whole genome level

MicroRNAs (miRNAs) are critical regulators of various biological and metabolic processes of plants. Numerous miRNAs and their functions have been identified and analyzed in many plants. However, till now, the involvement of miRNAs in the response of grapevine berries to ethylene has not been reported yet. Here, Solexa technology was employed to deeply sequence small RNA libraries constructed from grapevine berries treated with and without ethylene. A total of 124 known and 78 novel miRNAs were identified. Among these miRNAs, 162 miRNAs were clearly responsive to ethylene, with 55 downregulated, 59 upregulated, and 14 unchanged miRNAs detected only in the control. The other 35 miRNAs responsive to ethylene were induced by ethylene and detected only in the ethylene-treated grapevine materials. Expression analysis of 27 conserved and 26 novel miRNAs revealed that 13 conserved and 18 novel ones were regulated by ethylene during the whole development of grapevine berries. High-throughput sequencing and qRT-PCR assays revealed consistent results on the expression results of ethylene-responsive miRNAs. Moreover, 90 target genes for 34 novel miRNAs were predicted, most of which were involved in responses to various stresses, especially like exogenous ethylene treatment. The identified miRNAs may be mainly involved in grapevine berry development and response to various environmental conditions.

Generation of different sizes and classes of small RNAs in barley is locus, chromosome and/or cultivar-dependent

Background

Various small RNA (sRNA) sizes and varieties have been identified, but their relationship as well as relationship with their origins and allocations have not been well understood or investigated.

Results

By comparing sRNAs generated from two barley cultivars, Golden Promise (GP) and Pallas, we identified that the generation of different sizes and types of sRNAs in barley was locus-, chromosome- and/or cultivar-dependent. 20-nt sRNAs mainly comprising miRNAs and chloroplast-derived sRNAs were significantly over-expressed in Pallas vs. GP on chromosomes 3H and 6H. MiRNAs-enriched 21-nt sRNAs were significantly over-expressed in Pallas vs. GP only on chromosome 4H. On chromosome 5H this size of sRNAs was significantly under-expressed in Pallas, so were 22-nt sRNAs mainly comprising miRNAs and repeat-derived sRNAs. 24-nt sRNAs mostly derived from repeats were evenly distributed in all chromosomes and expressed similarly between GP and Pallas. Unlike other sizes of sRNAs, 24-nt sRNAs were little conserved in other plant species. Abundant sRNAs were mostly generated from 3’ terminal regions of chromosome 1H and 5’ terminal regions of chromosome 5H. Over-expressed miRNAs in GP vs. Pallas primarily function in stress responses and iron-binding.

Conclusions

Our study indicates that 23−24-nt sRNAs may be linked to repressive chromatin modifications and function in genome stability while 20−21-nt sRNAs may be important for the cultivar specificity. This study provides a novel insight into the mechanism of sRNA expression and function in barley.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3023-5) contains supplementary material, which is available to authorized users.

Time-Course Small RNA Profiling Reveals Rice miRNAs and Their Target Genes in Response to Rice Stripe Virus Infection

It has been known that many microRNAs (miRNAs) are involved in the regulation for the plant development and defense mechanism by regulating the expression of the target gene. Several previous studies has demonstrated functional roles of miRNAs in antiviral defense mechanisms. In this study, we employed high-throughput sequencing technology to identify rice miRNAs upon rice stripe virus (RSV) infection at three different time points. Six libraries from mock and RSV-infected samples were subjected for small RNA sequencing. Bioinformatic analyses revealed 374 known miRNAs and 19 novel miRNAs. Expression of most identified miRNAs was not dramatically changed at 3 days post infection (dpi) and 7 dpi by RSV infection. However, many numbers of miRNAs were up-regulated in mock and RSV-infected samples at 15 dpi by RSV infection. Moreover, expression profiles of identified miRNAs revealed that only few numbers of miRNAs were strongly regulated by RSV infection. In addition, 15 resistance genes were targets of six miRNAs suggesting that those identified miRNAs and 15 NBS-LRR resistance genes might be involved in RSV infection. Taken together, our results provide novel insight into the dynamic expression profiles of rice miRNAs upon RSV infection and clues for the understanding of the regulatory roles of miRNAs via time-course.

Identification of Ice Plant (Mesembryanthemum crystallinum L.) MicroRNAs Using RNA-Seq and Their Putative Roles in High Salinity Responses in Seedlings

The halophyte Mesembryanthemum crystallinum (common or crystalline ice plant) is a useful model for studying molecular mechanisms of salt tolerance. The morphology, physiology, metabolism, and gene expression of ice plant have been studied and large-scale analyses of gene expression profiling have drawn an outline of salt tolerance in ice plant. A rapid root growth to a sudden increase in salinity was observed in ice plant seedlings. Using a fluorescent dye to detect Na(+), we found that ice plant roots respond to an increased flux of Na(+) by either secreting or storing Na(+) in specialized cells. High-throughput sequencing was used to identify small RNA profiles in 3-day-old seedlings treated with or without 200 mM NaCl. In total, 135 conserved miRNAs belonging to 21 families were found. The hairpin precursor of 19 conserved mcr-miRNAs and 12 novel mcr-miRNAs were identified. After 6 h of salt stress, the expression of most mcr-miRNAs showed decreased relative abundance, whereas the expression of their corresponding target genes showed increased mRNA relative abundance. The cognate target genes are involved in a broad range of biological processes: transcription factors that regulate growth and development, enzymes that catalyze miRNA biogenesis for the most conserved mcr-miRNA, and proteins that are involved in ion homeostasis and drought-stress responses for some novel mcr-miRNAs. Analyses of the functions of target genes revealed that cellular processes, including growth and development, metabolism, and ion transport activity are likely to be enhanced in roots under salt stress. The expression of eleven conserved miRNAs and two novel miRNAs were correlated reciprocally with predicted targets within hours after salt stress exposure. Several conserved miRNAs have been known to regulate root elongation, root apical meristem activity, and lateral root formation. Based upon the expression pattern of miRNA and target genes in combination with the observation of Na(+) distribution, ice plant likely responds to increased salinity by using Na(+) as an osmoticum for cell expansion and guard cell opening. Excessive Na(+) could either be secreted through the root epidermis or stored in specialized leaf epidermal cells. These responses are regulated in part at the miRNA-mediated post-transcriptional level.

Identification of novel miRNAs from drought tolerant rice variety Nagina 22

MicroRNAs regulate a spectrum of developmental and biochemical processes in plants and animals. Thus, knowledge of the entire miRNome is essential to understand the complete regulatory schema of any organism. The current study attempts to unravel yet undiscovered miRNA genes in rice. Analysis of small RNA libraries from various tissues of drought-tolerant ‘aus’ rice variety Nagina 22 (N22) identified 71 novel miRNAs. These were validated based on precursor hairpin structure, small RNA mapping pattern, ‘star’ sequence, conservation and identification of targets based on degradome data. While some novel miRNAs were conserved in other monocots and dicots, most appear to be lineage-specific. They were segregated into two different classes based on the closeness to the classical miRNA definition. Interestingly, evidence of a miRNA-like cleavage was found even for miRNAs that lie beyond the classical definition. Several novel miRNAs displayed tissue-enriched and/or drought responsive expression. Generation and analysis of the degradome data from N22 along with publicly available degradome identified several high confidence targets implicated in regulation of fundamental processes such as flowering and stress response. Thus, discovery of these novel miRNAs considerably expands the dimension of the miRNA-mediated regulation in rice.

Genomics Approaches For Improving Salinity Stress Tolerance in Crop Plants

Salinity is one of the major factors which reduces crop production worldwide. Plant responses to salinity are highly complex and involve a plethora of genes. Due to its multigenicity, it has been difficult to attain a complete understanding of how plants respond to salinity. Genomics has progressed tremendously over the past decade and has played a crucial role towards providing necessary knowledge for crop improvement. Through genomics, we have been able to identify and characterize the genes involved in salinity stress response, map out signaling pathways and ultimately utilize this information for improving the salinity tolerance of existing crops. The use of new tools, such as gene pyramiding, in genetic engineering and marker assisted breeding has tremendously enhanced our ability to generate stress tolerant crops. Genome editing technologies such as Zinc finger nucleases, TALENs and CRISPR/Cas9 also provide newer and faster avenues for plant biologists to generate precisely engineered crops.

Identification and characterization of a novel NAC-like gene in chrysanthemum (Dendranthema lavandulifolium)

KEY MESSAGE

A NAC -like gene named DlNAC1 was identified in chrysanthemum and characterized; it may be involved in regulation of response to abiotic stressors, especially in tolerance to drought and salinity. NAC transcription factors in plants play crucial roles in tolerance to abiotic stressors, and overexpression of the NAC gene in Arabidopsis has been demonstrated to lead to improved drought tolerance. Functions of the NAC genes in chrysanthemum, however, remain poorly understood. In this study, a NAC-like gene named DlNAC1 was identified in chrysanthemum (Dendranthema lavandulifolium) and characterized. Phylogenetic analysis indicated that DlNAC1 contains a typical NAC domain and belongs to the ONAC022 subgroup. According to the subcellular localization and yeast one-hybrid assay, the DlNAC1 protein is localized to nuclei and has a transcription activation ability. Moreover, quantitative real-time PCR analyses showed that DlNAC1 was induced by low-temperature, high-salinity, and drought conditions (separately), but not by abscisic acid (ABA) and heat shock. In these experiments, the downstream genes of NAC transcription factors were found to be up-regulated, including stress-responsive genes KIN1 and AMY1. To further explore the effects of DlNAC1 in response to abiotic stressors, DlNAC1 was overexpressed in tobacco, and these transgenic plants showed significantly enhanced tolerance to drought and salinity. This study suggests that in chrysanthemum, the DlNAC1 gene is involved in regulation of the response to abiotic stressors, especially in tolerance to drought and salinity.

Identification of Known and Novel microRNAs and Their Targets in Peach (Prunus persica) Fruit by High-Throughput Sequencing

MicroRNAs (miRNAs) are a group of non-coding RNAs that have functions in post-transcriptional gene regulation in plants. Although the most important economic component of peach trees (Prunus persica) is the fruit, not much is known about miRNAs in this organ. In this study, miRNAs and their targets were identified and characterized from libraries of small RNAs of peach fruit through Solexa based-sequencing and bioinformatics approaches. A total of 557 known peach miRNAs belonging to 34 miRNA families were identified, and some of these miRNAs were found to be highly conserved in at least four other plant species. Using the most current criteria for miRNA annotation, 275 putative novel miRNAs were predicted, and the sequencing frequencies of these novel miRNAs were less than those of the conserved miRNAs. In total, 3959 and 1614 target genes for 349 known and 193 novel miRNAs, respectively, were predicted with the criteria that a single target gene can be targeted by different miRNAs and that a single miRNA can also have a large number of target genes. Three targets were even found to be targeted by 13 novel miRNAs that contained the same complete miRNA sequence at different locations and had different scaffolds. The proteins predicted to be targeted by the miRNAs identified in this study encompass a wide range of transcription factors and are involved in many biological processes and pathways, including development, metabolism, stress responses and signal transduction. A total of 115 and 101 target genes were identified to be cleaved by 60 known miRNAs and 27 novel miRNAs through degradome sequencing, respectively. These miRNAs induce cleavage of their targets precisely at the position between nucleotides 10 and 11 of the miRNA sequences from the 5’ to the 3’ end. Thirty conserved miRNAs and 19 novel miRNAs exhibited differential expression profiles in the peach, and the expression patterns of some miRNAs appeared to be tissue- or developmental stage-specific. The findings of this study provide an important basis for the analysis of miRNAs, their targets and the functions of these targets in peach fruit.

Endogenous Small-Noncoding RNAs and Potential Functions in Desiccation Tolerance in Physcomitrella Patens

Early land plants like moss Physcomitrella patens have developed remarkable drought tolerance. Phytohormone abscisic acid (ABA) protects seeds during water stress by activating genes through transcription factors such as ABSCISIC ACID INSENSITIVE (ABI3). Small noncoding RNA (sncRNA), including microRNAs (miRNAs) and endogenous small-interfering RNAs (endo-siRNAs), are key gene regulators in eukaryotes, playing critical roles in stress tolerance in plants. Combining next-generation sequencing and computational analysis, we profiled and characterized sncRNA species from two ABI3 deletion mutants and the wild type P. patens that were subject to ABA treatment in dehydration and rehydration stages. Small RNA profiling using deep sequencing helped identify 22 novel miRNAs and 6 genomic loci producing trans-acting siRNAs (ta-siRNAs) including TAS3a to TAS3e and TAS6. Data from degradome profiling showed that ABI3 genes (ABI3a/b/c) are potentially regulated by the plant-specific miR536 and that other ABA-relevant genes are regulated by miRNAs and ta-siRNAs. We also observed broad variations of miRNAs and ta-siRNAs expression across different stages, suggesting that they could potentially influence desiccation tolerance. This study provided evidence on the potential roles of sncRNA in mediating desiccation-responsive pathways in early land plants.

High-throughput sequencing reveals differential regulation of miRNAs in fenoxaprop-P-ethyl-resistant Beckmannia syzigachne

Non-target site resistance (NTSR) to herbicides is an increasing concern for weed control. The majority of previous studies have focused on metabolic resistance mechanisms of NTSR, but no research exists on gene regulation mechanisms behind herbicide resistance, such as microRNA (miRNA). Here, we identified 3 American sloughgrass (Beckmannia syzigachne Steud.) populations containing fenoxaprop-P-ethyl-resistant plants. We then constructed small RNA libraries and subjected them to deep sequencing and bioinformatics analyses. Forty known and 36 potentially novel, predicted miRNAs were successfully identified. Of these, we identified 3 conserved, predicted candidate NTSR-determinant miRNAs and their potential corresponding target genes, as well as 4 novel potential miRNAs with high count. Target gene prediction and annotation indicated that these 7 differentially expressed miRNAs potentially play a role in regulating specific stress-responsive genes, very likely related to herbicide resistance. Expression profiles were determined with quantitative real-time PCR. The present study is a novel, large-scale characterization of weed miRNAs. The results should further our understanding of miRNA expression profiles associated with herbicide resistance, allowing for the development of more effective weed management strategies.

MicroRNA Signatures of Drought Signaling in Rice Root

Background

Drought stress is one of the most important abiotic stresses and the main constraint to rice agriculture. MicroRNA-mediated post-transcriptional gene regulation is one of the ways to establish drought stress tolerance in plants. MiRNAs are 20–24-nt regulatory RNAs that play an important role in regulating plant gene expression upon exposure to biotic and abiotic stresses.

Methodology/Principal Findings

In this study, we applied a partial root drying system as well as a complete root drying system to identify miRNAs involved in conditions of drought stress, drought signaling and wet signaling using high-throughput sequencing. To this end, we produced four small RNA libraries: (1) fully-watered (WW), (2) fully-droughted (WD), and split-root systems where (3) one-half was well watered (SpWW) and (4) the other half was water-deprived (SpWD). Our analysis revealed 10,671 and 783 unique known and novel miRNA reads in all libraries, respectively. We identified, 65 (52 known + 13 novel), 72 (61 known + 11 novel) and 51 (38 known + 13 novel) miRNAs that showed differential expression under conditions of drought stress, drought signaling and wet signaling, respectively. The results of quantitative real-time PCR showed expression patterns similar to the high-throughput sequencing results. Furthermore, our target prediction led to the identification of 244, 341 and 239 unique target genes for drought-stress-, drought-signaling- and wet-signaling-responsive miRNAs, respectively.

Conclusions/Significance

Our results suggest that miRNAs that are responsive under different conditions could play different roles in the regulation of abscisic acid signaling, calcium signaling, detoxification and lateral root formation.

Characterization and primary functional analysis of a bamboo NAC gene targeted by miR164b

PeSNAC1 , a stress-related NAC1 from Phyllostachys edulis , was characterized. Ectopic expression in Arabidopsis indicated that PeSNAC1 together with ped -miR164b participated in the regulation of organ boundaries and stress tolerance. NAC (NAM, ATAF1/2 and CUC2) participates in many different processes regulating plant growth, development, and stress response. A total of 125 NAC genes have been predicted in moso bamboo (Phyllostachys edulis), but their roles are poorly understood. PeSNAC1 targeted by ped-miR164b was focused for further study. The cleavage of PeSNAC1 mRNA guided by ped-miR164b was validated using RLM-5' RACE. Tissue-specific expression analysis demonstrated that ped-miR164b had a declining trend from root, sheath, leaf, to that of stem, which was opposite to that of PeSNAC1. Transgenic Arabidopsis plants overexpressing either PeSNAC1 (OX-PeSNAC1) or, ped-miR164b (OX-ped-miR164b) driven by the CaMV35S promoter were generated. OX-ped-miR164b plants showed similar phenotype of cuc2 mutants whose growth was seriously suppressed. Compared with Col-0, sense OX-PeSNAC1 plants grew rapidly and flowered earlier, whereas antisense plants grew slowly and exhibited delayed flowering. Sense OX-PeSNAC1 plants had the greatest number of lateral roots, while antisense OX-PeSNAC1 and OX-ped-miR164b plants had fewer lateral roots than Col-0. Under NaCl and PEG6000 stresses, survival rates were higher and F v/F m values declined more slowly in sense OX-PeSNAC1 plants than in Col-0, with lower survival rates and a more rapid decrease in F v/F m values conversely observed in antisense OX-PeSNAC1 and OX-ped-miR164b plants. These findings indicated that ped-miR164b-targeted PeSNAC1 may play key roles in plant development and tolerance to salinity and drought stresses.

Genome-Wide Identification and Characterization of microRNAs in Developing Grains of Zea mays L

The development and maturation of maize kernel involves meticulous and fine gene regulation at transcriptional and post-transcriptional levels, and miRNAs play important roles during this process. Although a number of miRNAs have been identified in maize seed, the ones involved in the early development of grains and in different lines of maize have not been well studied. Here, we profiled four small RNA libraries, each constructed from groups of immature grains of Zea mays inbred line Chang 7–2 collected 4–6, 7–9, 12–14, and 18–23 days after pollination (DAP). A total of 40 known (containing 111 unique miRNAs) and 162 novel (containing 196 unique miRNA candidates) miRNA families were identified. For conserved and novel miRNAs with over 100 total reads, 44% had higher accumulation before the 9^th DAP, especially miR166 family members. 42% of miRNAs had highest accumulation during 12–14 DAP (which is the transition stage from embryogenesis to nutrient storage). Only 14% of miRNAs had higher expression 18–23 DAP. Prediction of potential targets of all miRNAs showed that 165 miRNA families had 377 target genes. For miR164 and miR166, we showed that the transcriptional levels of their target genes were significantly decreased when co-expressed with their cognate miRNA precursors in vivo. Further analysis shows miR159, miR164, miR166, miR171, miR390, miR399, and miR529 families have putative roles in the embryogenesis of maize grain development by participating in transcriptional regulation and morphogenesis, while miR167 and miR528 families participate in metabolism process and stress response during nutrient storage. Our study is the first to present an integrated dynamic expression pattern of miRNAs during maize kernel formation and maturation.

Combined small RNA and degradome sequencing to identify miRNAs and their targets in response to drought in foxtail millet

Background

Foxtail millet (Setaria italica) is a diploid C₄ panicoid species. Because of its prominent drought resistance, small genome size, self-pollination, and short life cycle, foxtail millet has become an ideal model system for studying drought tolerance of crops. MicroRNAs (miRNAs) are endogenous, small RNAs that play important regulatory roles in the development and stress response in plants.

Results

In this study, we applied Illumina sequencing to systematically investigate the drought-responsive miRNAs derived from S. italica inbred An04-4783 seedlings grown under control and drought conditions. Degradome sequencing was applied to confirm the targets of these miRNAs at a global level. A total of 81 known miRNAs belonging to 28 families were identified, among which 14 miRNAs were upregulated and four were downregulated in response to drought. In addition, 72 potential novel miRNAs were identified, three of which were differentially expressed under drought conditions. Degradome sequencing analysis showed that 56 and 26 genes were identified as targets of known and novel miRNAs, respectively.

Conclusions

Our analysis revealed post-transcriptional remodeling of cell development, transcription factors, ABA signaling, and cellar homeostasis in S.italica in response to drought. This preliminary characterization provided useful information for further studies on the regulatory networks of drought-responsive miRNAs in foxtail millet.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0364-7) contains supplementary material, which is available to authorized users.

Wheat miRNA ancestors: evident by transcriptome analysis of A, B, and D genome donors

MicroRNAs are critical players of post-transcriptional gene regulation with profound effects on the fundamental processes of cellular life. Their identification and characterization, together with their targets, hold great significance in exploring and exploiting their roles on a functional context, providing valuable clues into the regulation of important biological processes, such as stress tolerance or environmental adaptation. Wheat is a hardy crop, extensively harvested in temperate regions, and is a major component of the human diet. With the advent of the next generation sequencing technologies considerably decreasing sequencing costs per base-pair, genomic, and transcriptomic data from several wheat species, including the progenitors and wild relatives have become available. In this study, we performed in silico identification and comparative analysis of microRNA repertoires of bread wheat (Triticum aestivum L.) and its diploid progenitors and relatives, Aegilops sharonensis, Aegilops speltoides, Aegilops tauschii, Triticum monococcum, and Triticum urartu through the utilization of publicly available transcriptomic data. Over 200 miRNA families were identified, majority of which have not previously been reported. Ancestral relationships expanded our understanding of wheat miRNA evolution, while T. monococcum miRNAs delivered important clues on the effects of domestication on miRNA expression. Comparative analyses on wild Ae. sharonensis accessions highlighted candidate miRNAs that can be linked to stress tolerance. The miRNA repertoires of bread wheat and its diploid progenitors and relatives provide important insight into the diversification and distribution of miRNA genes, which should contribute to the elucidation of miRNA evolution of Poaceae family. A thorough understanding of the convergent and divergent expression profiles of miRNAs in different genetic backgrounds can provide unique opportunities to modulation of gene regulation for better crop performance.

Expression Profile Analysis of miR-221 and miR-222 in Different Tissues and Head Kidney Cells of Cynoglossus semilaevis, Following Pathogen Infection

Half-smooth tongue sole (Cynoglossus semilaevis) is an important marine commercial fish species in China, which suffers from widespread disease outbreaks. Recently, in this regard, our group identified immune-related microRNAs (miRNAs) of C. semilaevis following Vibrio anguillarum infection. Furthermore, miRNA microarray was utilized to characterize the immune roles of important miRNA candidates in response to bacterial infection. Therefore, in the present study, we characterized miR-221 and miR-222 and profiled their expression after challenge. Here, miR-221 and miR-222 precursors were predicted to have a typical hairpin structure. Both miRNAs were expressed in a broad range of tissues in C. semilaevis, while miR-221 and miR-222 were significantly differentially expressed in the immune tissues of C. semilaevis among three small RNA libraries [control group (CG), bacteria-challenged fish without obvious symptoms of infection (NOSG), and bacteria-challenged fish with obvious symptoms of infection (HOSG)]. In order to further characterize and understand the immune response of miR-221 and miR-222, therefore, we profiled miR-221 and miR-222 expression in selected immune tissues after challenge with V. anguillarum. Both miR-221 and miR-222 were upregulated in the liver and spleen, while different expression patterns were observed in the head kidney. In addition, in half-smooth tongue sole head kidney cell line after challenge with lipopolysaccharide (LPS), polyinosinic:polycytidylic acid (poly I:C), peptidoglycan (PGN), and red-spotted grouper nervous necrosis virus (RGNNV), both miR-221 and miR-222 showed significant difference in expression response to pathogen. Meanwhile, the target gene of miR-221 and miR-222 was predicted, which indicated that tumor necrosis factor receptor-associated factor 6 (TRAF6) and interleukin-1 beta (IL-1β) were the target genes of miR-221 and miR-222, respectively. Collectively, these findings indicated that miR-221 and miR-222 have putative roles in innate immune response during C. semilaevis exposure to pathogens. Our findings could expand the knowledge of immune function of C. semilaevis miRNA and guide future studies on C. semilaevis immunity.

A Comprehensive Prescription for Plant miRNA Identification

microRNAs (miRNAs) are tiny ribo-regulatory molecules involved in various essential pathways for persistence of cellular life, such as development, environmental adaptation, and stress response. In recent years, miRNAs have become a major focus in molecular biology because of their functional and diagnostic importance. This interest in miRNA research has resulted in the development of many specific software and pipelines for the identification of miRNAs and their specific targets, which is the key for the elucidation of miRNA-modulated gene expression. While the well-recognized importance of miRNAs in clinical research pushed the emergence of many useful computational identification approaches in animals, available software and pipelines are fewer for plants. Additionally, existing approaches suffers from mis-identification and annotation of plant miRNAs since the miRNA mining process for plants is highly prone to false-positives, particularly in cereals which have a highly repetitive genome. Our group developed a homology-based in silico miRNA identification approach for plants, which utilizes two Perl scripts "SUmirFind" and "SUmirFold" and since then, this method helped identify many miRNAs particularly from crop species such as Triticum or Aegliops. Herein, we describe a comprehensive updated guideline by the implementation of two new scripts, "SUmirPredictor" and "SUmirLocator," and refinements to our previous method in order to identify genuine miRNAs with increased sensitivity in consideration of miRNA identification problems in plants. Recent updates enable our method to provide more reliable and precise results in an automated fashion in addition to solutions for elimination of most false-positive predictions, miRNA naming and miRNA mis-annotation. It also provides a comprehensive view to genome/transcriptome-wide location of miRNA precursors as well as their association with transposable elements. The "SUmirPredictor" and "SUmirLocator" scripts are freely available together with a reference high-confidence plant miRNA list.

Novel and conserved miRNAs in the halophyte Suaeda maritima identified by deep sequencing and computational predictions using the ESTs of two mangrove plants

BACKGROUND

Although miRNAs are reportedly involved in the salt stress tolerance of plants, miRNA profiling in plants has largely remained restricted to glycophytes, including certain crop species that do not exhibit any tolerance to salinity. Hence, this manuscript describes the results from the miRNA profiling of the halophyte Suaeda maritima, which is used worldwide to study salt tolerance in plants.

RESULTS

A total of 134 conserved miRNAs were identified from unique sRNA reads, with 126 identified using miRBase 21.0 and an additional eight identified using the Plant Non-coding RNA Database. The presence of the precursors of seven conserved miRNAs was validated in S. maritima. In addition, 13 novel miRNAs were predicted using the ESTs of two mangrove plants, Rhizophora mangle and Heritiera littoralis, and the precursors of seven miRNAs were found in S. maritima. Most of the miRNAs considered for characterization were responsive to NaCl application, indicating their importance in the regulation of metabolic activities in plants exposed to salinity. An expression study of the novel miRNAs in plants of diverse ecological and taxonomic groups revealed that two of the miRNAs, sma-miR6 and sma-miR7, were also expressed in Oryza sativa, whereas another two, sma-miR2 and sma-miR5, were only expressed in plants growing under the influence of seawater, similar to S. maritima.

CONCLUSION

The distribution of conserved miRNAs among only 25 families indicated the possibility of identifying a greater number of miRNAs with increase in knowledge of the genomes of more halophytes. The expression of two novel miRNAs, sma-miR2 and sma-miR5, only in plants growing under the influence of seawater suggested their metabolic regulatory roles specific to saline environments, and such behavior might be mediated by alterations in the expression of certain genes, modifications of proteins leading to changes in their activity and production of secondary metabolites as revealed by the miRNA target predictions. Moreover, the auxin responsive factor targeted by sma-miR7 could also be involved in salt tolerance because the target is conserved between species. This study also indicated that the transcriptome of one species can be successfully used to computationally predict the miRNAs in other species, especially those that have similar metabolism, even if they are taxonomically separated.

Conserved and highly expressed tRNA derived fragments in zebrafish

BACKGROUND

Small non-coding RNAs (sncRNAs) are a class of transcripts implicated in several eukaryotic regulatory mechanisms, namely gene silencing and chromatin regulation. Despite significant progress in their identification by next generation sequencing (NGS) we are still far from understanding their full diversity and functional repertoire.

RESULTS

Here we report the identification of tRNA derived fragments (tRFs) by NGS of the sncRNA fraction of zebrafish. The tRFs identified are 18-30 nt long, are derived from specific 5' and 3' processing of mature tRNAs and are differentially expressed during development and in differentiated tissues, suggesting that they are likely produced by specific processing rather than random degradation of tRNAs. We further show that a highly expressed tRF (5'tRF-Pro(CGG)) is cleaved in vitro by Dicer and has silencing ability, indicating that it can enter the RNAi pathway. A computational analysis of zebrafish tRFs shows that they are conserved among vertebrates and mining of publicly available datasets reveals that some 5'tRFs are differentially expressed in disease conditions, namely during infection and colorectal cancer.

CONCLUSIONS

tRFs constitute a class of conserved regulatory RNAs in vertebrates and may be involved in mechanisms of genome regulation and in some diseases.

Identification of microRNAs and their targets in Finger millet by high throughput sequencing

MicroRNAs are short non-coding RNAs which play an important role in regulating gene expression by mRNA cleavage or by translational repression. The majority of identified miRNAs were evolutionarily conserved; however, others expressed in a species-specific manner. Finger millet is an important cereal crop; nonetheless, no practical information is available on microRNAs to date. In this study, we have identified 95 conserved microRNAs belonging to 39 families and 3 novel microRNAs by high throughput sequencing. For the identified conserved and novel miRNAs a total of 507 targets were predicted. 11 miRNAs were validated and tissue specificity was determined by stem loop RT-qPCR, Northern blot. GO analyses revealed targets of miRNA were involved in wide range of regulatory functions. This study implies large number of known and novel miRNAs found in Finger millet which may play important role in growth and development.

Role of bioinformatics in establishing microRNAs as modulators of abiotic stress responses: the new revolution

microRNAs (miRs) are a class of 21-24 nucleotide long non-coding RNAs responsible for regulating the expression of associated genes mainly by cleavage or translational inhibition of the target transcripts. With this characteristic of silencing, miRs act as an important component in regulation of plant responses in various stress conditions. In recent years, with drastic change in environmental and soil conditions different type of stresses have emerged as a major challenge for plants growth and productivity. The identification and profiling of miRs has itself been a challenge for research workers given their small size and large number of many probable sequences in the genome. Application of computational approaches has expedited the process of identification of miRs and their expression profiling in different conditions. The development of High-Throughput Sequencing (HTS) techniques has facilitated to gain access to the global profiles of the miRs for understanding their mode of action in plants. Introduction of various bioinformatics databases and tools have revolutionized the study of miRs and other small RNAs. This review focuses the role of bioinformatics approaches in the identification and study of the regulatory roles of plant miRs in the adaptive response to stresses.

Identification of novel and conserved microRNAs in Panax notoginseng roots by high-throughput sequencing

BACKGROUND

MicroRNAs (miRNAs) are small, non-coding RNAs that are important regulators of gene expression, and play major roles in plant development and their response to the environment. Root extracts from Panax notoginseng contain triterpene saponins as their principal bioactive constituent, and demonstrate medicinal properties. To investigate the novel and conserved miRNAs in P. notoginseng, three small RNA libraries constructed from 1-, 2-, and 3-year-old roots in which root saponin levels vary underwent high-throughput sequencing.

METHODS

P. notoginseng roots, purified from 1-, 2-, and 3-year-old roots, were extracted for RNA, respectively. Three small libraries were constructed and subjected to next generation sequencing.

RESULTS

Sequencing of the three libraries generated 67,217,124 clean reads from P. notoginseng roots. A total of 316 conserved miRNAs (belonging to 67 miRNA families and one unclassified family) and 52 novel miRNAs were identified. MIR156 and MIR166 were the largest miRNA families, while miR156i and miR156g showed the highest abundance of miRNA species. Potential miRNA target genes were predicted and annotated using Cluster of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes. Comparing these miRNAs between root samples revealed 33 that were differentially expressed between 2- and 1-year-old roots (8 increased, 25 decreased), 27 differentially expressed between 3- and 1-year-old roots (7 increased, 20 decreased), and 29 differentially expressed between 3- and 2-year-old roots (8 increased, 21 decreased). Two significantly differentially expressed miRNAs and four miRNAs predicted to target genes involved in the terpenoid backbone biosynthesis pathway were selected and validated by quantitative reverse transcription PCR. Furthermore, the expression patterns of these six miRNAs were analyzed in P. notoginseng roots, stems, and leaves at different developmental stages.

CONCLUSIONS

This study identified a large number of P. notoginseng miRNAs and their target genes, functional annotations, and gene expression patterns. It provides the first known miRNA profiles of the P. notoginseng root development cycle.

Identification of microRNAs in the Toxigenic Dinoflagellate Alexandrium catenella by High-Throughput Illumina Sequencing and Bioinformatic Analysis

Micro-ribonucleic acids (miRNAs) are a large group of endogenous, tiny, non-coding RNAs consisting of 19–25 nucleotides that regulate gene expression at either the transcriptional or post-transcriptional level by mediating gene silencing in eukaryotes. They are considered to be important regulators that affect growth, development, and response to various stresses in plants. Alexandrium catenella is an important marine toxic phytoplankton species that can cause harmful algal blooms (HABs). To date, identification and function analysis of miRNAs in A. catenella remain largely unexamined. In this study, high-throughput sequencing was performed on A. catenella to identify and quantitatively profile the repertoire of small RNAs from two different growth phases. A total of 38,092,056 and 32,969,156 raw reads were obtained from the two small RNA libraries, respectively. In total, 88 mature miRNAs belonging to 32 miRNA families were identified. Significant differences were found in the member number, expression level of various families, and expression abundance of each member within a family. A total of 15 potentially novel miRNAs were identified. Comparative profiling showed that 12 known miRNAs exhibited differential expression between the lag phase and the logarithmic phase. Real-time quantitative RT-PCR (qPCR) was performed to confirm the expression of two differentially expressed miRNAs that were one up-regulated novel miRNA (aca-miR-3p-456915), and one down-regulated conserved miRNA (tae-miR159a). The expression trend of the qPCR assay was generally consistent with the deep sequencing result. Target predictions of the 12 differentially expressed miRNAs resulted in 1813target genes. Gene ontology (GO) analysis and the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG) annotations revealed that some miRNAs were associated with growth and developmental processes of the alga. These results provide insights into the roles that miRNAs play in the growth of A. catenella, and they provide the basis for further studies of the molecular mechanisms that underlie bloom growth in red tides species.

Identification of four functionally important microRNA families with contrasting differential expression profiles between drought-tolerant and susceptible rice leaf at vegetative stage

Background

Developing drought-tolerant rice varieties with higher yield under water stressed conditions provides a viable solution to serious yield-reduction impact of drought. Understanding the molecular regulation of this polygenic trait is crucial for the eventual success of rice molecular breeding programmes. microRNAs have received tremendous attention recently due to its importance in negative regulation. In plants, apart from regulating developmental and physiological processes, microRNAs have also been associated with different biotic and abiotic stresses. Hence here we chose to analyze the differential expression profiles of microRNAs in three drought treated rice varieties: Vandana (drought-tolerant), Aday Sel (drought-tolerant) and IR64 (drought-susceptible) in greenhouse conditions via high-throughput sequencing.

Results

Twenty-six novel microRNA candidates involved in the regulation of diverse biological processes were identified based on the detection of miRNA*. Out of their 110 predicted targets, we confirmed 16 targets from 5 novel microRNA candidates. In the differential expression analysis, mature microRNA members from 49 families of known Oryza sativa microRNA were differentially expressed in leaf and stem respectively with over 28 families having at least a similar mature microRNA member commonly found to be differentially expressed between both tissues. Via the sequence profiling data of leaf samples, we identified osa-miR397a/b, osa-miR398b, osa-miR408-5p and osa-miR528-5p as being down-regulated in two drought-tolerant rice varieties and up-regulated in the drought-susceptible variety. These microRNAs are known to be involved in regulating starch metabolism, antioxidant defence, respiration and photosynthesis. A wide range of biological processes were found to be regulated by the target genes of all the identified differentially expressed microRNAs between both tissues, namely root development (5.3–5.7 %), cell transport (13.2–18.4 %), response to stress (10.5–11.3 %), lignin catabolic process (3.8–5.3 %), metabolic processes (32.1–39.5 %), oxidation-reduction process (9.4–13.2 %) and DNA replication (5.7–7.9 %). The predicted target genes of osa-miR166e-3p, osa-miR166h-5p*, osa-miR169r-3p* and osa-miR397a/b were found to be annotated to several of the aforementioned biological processes.

Conclusions

The experimental design of this study, which features rice varieties with different drought tolerance and tissue specificity (leaf and stem), has provided new microRNA profiling information. The potentially regulatory importance of the microRNA genes mentioned above and their target genes would require further functional analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1851-3) contains supplementary material, which is available to authorized users.

Promoter methylation regulates the abundance of osa-miR393a in contrasting rice genotypes under salinity stress

MicroRNAs (miRNAs) are important molecules that regulate gene expression under salinity stress. Despite their evolutionary conservation, these regulatory elements have been shown to behave differently in different plant species under a particular environmental stress. In this study, we investigated the behavior of salt responsive osa-miR393a and its target gene (TIR1, LOC_Os05g05800) in salt-tolerant (FL478) and salt-sensitive (IR29) rice genotypes. It was found that the mature and precursor sequences of osa-miR393a as well as its cleavage site in TIR1 were conserved among salt tolerant and sensitive genotypes. Promoters of different salt-responsive miRNAs were also found to be less variable between salt-tolerant and salt-susceptible genotypes. Analysis of gene expression, promoter methylation, and cis-element abundance showed that osa-miR393a behaves differently in FL478 and IR29. Salt stress altered the expression pattern of osa-miR393a-TIR1 module in a time-dependent manner in the roots and shoots of two genotypes. Promoter methylation of this regulatory module was also altered at different time points under salt stress. Expression analysis in two genotypes indicated the overall down-regulation of osa-miR393a and up-regulation of TIR1 in FL478 and their reciprocal regulation in IR29. The expression results were complemented by the differential promoter methylation and cis-element abundance of this regulatory module. Together, the results of transcript abundance and promoter methylation of osa-miR393a-TIR1 module signified the association between these two processes which is reported for the first time in plants to the best of our knowledge.

Root precursors of microRNAs in wild emmer and modern wheats show major differences in response to drought stress

MicroRNAs, small regulatory molecules with significant impacts on the transcriptional network of all living organisms, have been the focus of several studies conducted mostly on modern wheat cultivars. In this study, we investigated miRNA repertoires of modern durum wheat and its wild relatives, with differing degrees of drought tolerance, to identify miRNA candidates and their targets involved in drought stress response. Root transcriptomes of Triticum turgidum ssp. durum variety Kızıltan and two Triticum turgidum ssp. dicoccoides genotypes TR39477 and TTD-22 under control and drought conditions were assembled from individual RNA-Seq reads and used for in silico identification of miRNAs. A total of 66 miRNAs were identified from all species, across all conditions, of which 46 and 38 of the miRNAs identified from modern durum wheat and wild genotypes, respectively, had not been previously reported. Genotype- and/or stress-specific miRNAs provide insights into our understanding of the complex drought response. Particularly, miR1435, miR5024, and miR7714, identified only from drought-stress roots of drought-tolerant genotype TR39477, can be candidates for future studies to explore and exploit the drought response to develop tolerant varieties.

miRNA-based drought regulation in wheat

MicroRNAs (miRNAs) are a class of small non-coding regulatory RNAs that regulate gene expression by guiding target mRNA cleavage or translational inhibition. Drought is a common environmental stress influencing crop growth and development. To date, it has been reported that a number of plant miRNA are involved in drought stress response. In this study, we comparatively investigated drought stress-responsive miRNAs in the root and leaf of bread wheat (Triticum aestivum cv. Sivas 111/33) by miRNA microarray screening. miRNA microarray analysis showed that 285 miRNAs (207 upregulated and 78 downregulated) and 244 miRNAs (115 upregulated and 129 downregulated) were differentially expressed in leaf and root tissues, respectively. Among the differentially expressed miRNAs, 23 miRNAs were only expressed in the leaf and 26 miRNAs were only expressed in the root of wheat growth under drought stress. Upon drought treatment, expression of miR159, miR160, miR166, miR169, miR172, miR395, miR396, miR408, miR472, miR477, miR482, miR1858, miR2118, and miR5049 were found to be significantly differentiated in bread wheat. The regulatory network analysis showed that miR395 has connections with a number of target transcripts, and miR159 and miR319 share a number of target genes. Drought-tolerant and drought-sensitive wheat cultivars showed altered expression pattern upon drought stress in terms of investigated miRNA and their target transcript expression level.

Genome-wide identification of turnip mosaic virus-responsive microRNAs in non-heading Chinese cabbage by high-throughput sequencing

Turnip mosaic virus (TuMV) is the most prevalent viral pathogen infecting most cruciferous plants. MicroRNAs (miRNAs) are around 22 nucleotides long non-protein-coding RNAs that play key regulatory roles in plants. Recent research findings show that miRNAs are involved in plant-virus interaction. However we know little about plant defense and viral offense system networks throughout microRNA regulation pathway. In this study, two small RNA libraries were constructed based on non-heading Chinese cabbage (Brassica campestris ssp. chinensis L. Makino, NHCC) leaves infected by TuMV and healthy leaves, and sequenced using the Illumina-Solexa high-throughput sequencing technology. A total of 86 conserved miRNAs belonging to 25 known miRNA families and 45 novel ones were identified. Among them, twelve conserved and ten new miRNAs were validated by real-time fluorescence quantitative PCR (qPCR). Differential expression analysis showed that 42 miRNAs were down-regulated and 27 miRNAs were up-regulated in response to TuMV stress. A total of 271 target genes were predicted using a bioinformatics approach, these genes are mainly involved in growth and resistance to various stresses. We further selected 13 miRNAs and their corresponding target genes to explore their expression pattern under TuMV and/or cold (4°C) stresses, and the results indicated that some of the identified miRNAs could link TuMV response with cold response of NHCC. The characterization of these miRNAs could contribute to a better understanding of plant-virus interaction throughout microRNA regulation pathway. This can lead to finding new approach to defend virus infection using miRNA in Chinese cabbage.

High throughput sequencing of small RNA component of leaves and inflorescence revealed conserved and novel miRNAs as well as phasiRNA loci in chickpea

Among legumes, chickpea (Cicer arietinum L.) is the second most important crop after soybean. MicroRNAs (miRNAs) play important roles by regulating target gene expression important for plant development and tolerance to stress conditions. Additionally, recently discovered phased siRNAs (phasiRNAs), a new class of small RNAs, are abundantly produced in legumes. Nevertheless, little is known about these regulatory molecules in chickpea. The small RNA population was sequenced from leaves and flowers of chickpea to identify conserved and novel miRNAs as well as phasiRNAs/phasiRNA loci. Bioinformatics analysis revealed 157 miRNA loci for the 96 highly conserved and known miRNA homologs belonging to 38 miRNA families in chickpea. Furthermore, 20 novel miRNAs belonging to 17 miRNA families were identified. Sequence analysis revealed approximately 60 phasiRNA loci. Potential target genes likely to be regulated by these miRNAs were predicted and some were confirmed by modified 5' RACE assay. Predicted targets are mostly transcription factors that might be important for developmental processes, and others include superoxide dismutases, plantacyanin, laccases and F-box proteins that could participate in stress responses and protein degradation. Overall, this study provides an inventory of miRNA-target gene interactions for chickpea, useful for the comparative analysis of small RNAs among legumes.

Identification and Characterization of ABA-Responsive MicroRNAs in Rice

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that silence genes through mRNA degradation or translational inhibition. The phytohormone abscisic acid (ABA) is essential for plant development and adaptation to abiotic and biotic stresses. In Arabidopsis, miRNAs are implicated in ABA functions. However, ABA-responsive miRNAs have not been systematically studied in rice. Here high throughput sequencing of small RNAs revealed that 107 miRNAs were differentially expressed in the rice ABA deficient mutant, Osaba1. Of these, 13 were confirmed by stem-loop RT-PCR. Among them, miR1425-5P, miR169a, miR169n, miR390-5P, miR397a and miR397b were up-regulated, but miR162b reduced in expression in Osaba1. The targets of these 13 miRNAs were predicted and validated by gene expression profiling. Interestingly, the expression levels of these miRNAs and their targets were regulated by ABA. Cleavage sites were detected on 7 of the miRNA targets by 5'-Rapid Amplification of cDNA Ends (5'-RACE). Finally, miR162b and its target OsTRE1 were shown to affect rice resistance to drought stress, suggesting that miR162b increases resistance to drought by targeting OsTRE1. Our work provides important information for further characterization and functional analysis of ABA-responsive miRNAs in rice.

Identification and Characterization of MicroRNAs in Ginkgo biloba var. epiphylla Mak

Ginkgo biloba, a dioecious plant known as a living fossil, is an ancient gymnosperm that stands distinct from other gymnosperms and angiosperms. Ginkgo biloba var. epiphylla (G. biloba var. epiphylla), with ovules borne on the leaf blade, is an unusual germplasm derived from G. biloba. MicroRNAs (miRNAs) are post-transcriptional gene regulators that play critical roles in diverse biological and metabolic processes. Currently, little is known about the miRNAs involved in the key stage of partly epiphyllous ovule germination in G. biloba var. epiphylla. Two small RNA libraries constructed from epiphyllous ovule leaves and normal leaves of G. biloba var. epiphylla were sequenced on an Illumina/Solexa platform. A total of 82 miRNA sequences belonging to 23 families and 53 putative novel miRNAs were identified in the two libraries. Differential expression analysis showed that 25 conserved and 21 novel miRNAs were differentially expressed between epiphyllous ovule leaves and normal leaves. The expression patterns of partially differentially expressed miRNAs and the transcript levels of their predicted target genes were validated by quantitative real time RT-PCR. All the expression profiles of the 21 selected miRNAs were similar to those detected by Solexa deep sequencing. Additionally, the transcript levels of almost all the putative target genes of 9 selected miRNAs were opposite to those of the corresponding miRNAs. The putative target genes of the differentially expressed miRNAs were annotated with Gene Ontology terms related to reproductive process, metabolic process and responding to stimulus. This work presents a broad range of small RNA transcriptome data obtained from epiphyllous ovule and normal leaves of G. biloba var. epiphylla, which may provide insights into the miRNA-mediated regulation in the epiphyllous ovule germination process.

Identification and characterisation of tobacco microRNA transcriptome using high-throughput sequencing

MicroRNAs (miRNAs) are post-transcriptional regulators that are involved in numerous biological processes in plants. In this study, we investigate miRNAs in Honghua Dajinyuan, an agronomically important species of tobacco in China. Here, we report a comprehensive analysis of miRNA expression profiles in the leaf, stem and root using a high-throughput sequencing approach. A total of 165 miRNAs, representing 55 conserved families, and 50 novel miRNAs, representing 19 families, were identified in three libraries. In addition, 12 miRNAs were randomly selected from a differentially expressed conserved miRNA family in three libraries with expression alterations and subjected to qRT-PCR validation. Of these, the expression level of nta-miR167d is highly enriched in the leaf tissue. In addition, the expression level of nta-miR319a is prominently enriched in the stem, while nta-miR160c is highly enriched in the root. Moreover, the target prediction showed that most of the targets coded for transcription factors that are involved in cellular and metabolic processes. GO analysis showed that most of the targets were involved in organelle function, served binding functions, and take part in cellular and metabolic processes. This study helps shed new light on understanding the role of miRNAs in different parts of the tobacco plant and adds a significant number of novel miRNAs to the tobacco miRNA transcriptome.

Small RNA deep sequencing identifies novel and salt-stress-regulated microRNAs from roots of Medicago sativa and Medicago truncatula

Small 21- to 24-nucleotide (nt) ribonucleic acids (RNAs), notably the microRNA (miRNA), are emerging as a posttranscriptional regulation mechanism. Salt stress is one of the primary abiotic stresses that cause the crop losses worldwide. In saline lands, root growth and function of plant are determined by the action of environmental salt stress through specific genes that adapt root development to the restrictive condition. To elucidate the role of miRNAs in salt stress regulation in Medicago, we used a high-throughput sequencing approach to analyze four small RNA libraries from roots of Zhongmu-1 (Medicago sativa) and Jemalong A17 (Medicago truncatula), which were treated with 300 mM NaCl for 0 and 8 h. Each library generated about 20 million short sequences and contained predominantly small RNAs of 24-nt length, followed by 21-nt and 22-nt small RNAs. Using sequence analysis, we identified 385 conserved miRNAs from 96 families, along with 68 novel candidate miRNAs. Of all the 68 predicted novel miRNAs, 15 miRNAs were identified to have miRNA*. Statistical analysis on abundance of sequencing read revealed specific miRNA showing contrasting expression patterns between M. sativa and M. truncatula roots, as well as between roots treated for 0 and 8 h. The expression of 10 conserved and novel miRNAs was also quantified by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). The miRNA precursor and target genes were predicted by bioinformatics analysis. We concluded that the salt stress related conserved and novel miRNAs may have a large variety of target mRNAs, some of which might play key roles in salt stress regulation of Medicago.

Identification of MicroRNAs and their Targets Associated with Embryo Abortion during Chrysanthemum Cross Breeding via High-Throughput Sequencing

Background

MicroRNAs (miRNAs) are important regulators in plant development. They post-transcriptionally regulate gene expression during various biological and metabolic processes by binding to the 3’-untranslated region of target mRNAs to facilitate mRNA degradation or inhibit translation. Chrysanthemum (Chrysanthemum morifolium) is one of the most important ornamental flowers with increasing demand each year. However, embryo abortion is the main reason for chrysanthemum cross breeding failure. To date, there have been no experiments examining the expression of miRNAs associated with chrysanthemum embryo development. Therefore, we sequenced three small RNA libraries to identify miRNAs and their functions. Our results will provide molecular insights into chrysanthemum embryo abortion.

Results

Three small RNA libraries were built from normal chrysanthemum ovules at 12 days after pollination (DAP), and normal and abnormal chrysanthemum ovules at 18 DAP. We validated 228 miRNAs with significant changes in expression frequency during embryonic development. Comparative profiling revealed that 69 miRNAs exhibited significant differential expression between normal and abnormal embryos at 18 DAP. In addition, a total of 1037 miRNA target genes were predicted, and their annotations were defined by transcriptome data. Target genes associated with metabolic pathways were most highly represented according to the annotation. Moreover, 52 predicted target genes were identified to be associated with embryonic development, including 31 transcription factors and 21 additional genes. Gene ontology (GO) annotation also revealed that high-ranking miRNA target genes related to cellular processes and metabolic processes were involved in transcription regulation and the embryo developmental process.

Conclusions

The present study generated three miRNA libraries and gained information on miRNAs and their targets in the chrysanthemum embryo. These results enrich the growing database of new miRNAs and lay the foundation for the further understanding of miRNA biological function in the regulation of chrysanthemum embryo abortion.

Signal transduction during wheat grain development

This review examines the signaling pathways from the developmental and environmental point of view and the interactions among external conditions, hormonal regulations, and sugarsensing in wheat. Grain development is the key phase of reproductive growth that is closely associated with vegetative organ senescence, initiation of grain filling, pre-stored assimilates remobilization, and maturation. Senescence is characterized by loss of chlorophyll and the degradation of proteins, nucleic acids, lipids as well as nutrient exports to the sink. The initiation and progression of vegetative organ senescence are under the control of an array of environmental signals (such as biotic and abiotic stresses, darkness, and nutrient availability) and endogenous factors (including aging, multiple hormones, and sugar availability). This review will discuss the major breakthroughs in signal transduction for the wheat (Triticum aestivum) grain development achieved in the past several years, with focuses on the regulation of senescence, reserves remobilization and biosynthesis of main components of the grain. Different mechanisms of diverse signals in controlling different phrases of wheat grain development, and cross talks between different signaling pathways will also be discussed. For perspectives, key signaling networks for grain development remain to be elucidated, including cross talks and the interactions between various environmental factors and internal signals.

Small RNA sequencing identifies miRNA roles in ovule and fibre development

MicroRNAs (miRNAs) have been found to be differentially expressed during cotton fibre development. However, which specific miRNAs and how they are involved in fibre development is unclear. Here, using deep sequencing, 65 conserved miRNA families were identified and 32 families were differentially expressed between leaf and ovule. At least 40 miRNAs were either leaf or ovule specific, whereas 62 miRNAs were shared in both leaf and ovule. qRT-PCR confirmed these miRNAs were differentially expressed during fibre early development. A total of 820 genes were potentially targeted by the identified miRNAs, whose functions are involved in a series of biological processes including fibre development, metabolism and signal transduction. Many predicted miRNA-target pairs were subsequently validated by degradome sequencing analysis. GO and KEGG analyses showed that the identified miRNAs and their targets were classified to 1027 GO terms including 568 biological processes, 324 molecular functions and 135 cellular components and were enriched to 78 KEGG pathways. At least seven unique miRNAs participate in trichome regulatory interaction network. Eleven trans-acting siRNA (tasiRNA) candidate genes were also identified in cotton. One has never been found in other plant species and two of them were derived from MYB and ARF, both of which play important roles in cotton fibre development. Sixteen genes were predicted to be tasiRNA targets, including sucrose synthase and MYB2. Together, this study discovered new miRNAs in cotton and offered evidences that miRNAs play important roles in cotton ovule/fibre development. The identification of tasiRNA genes and their targets broadens our understanding of the complicated regulatory mechanism of miRNAs in cotton.

Identification of reference genes for quantitative expression analysis of microRNAs and mRNAs in barley under various stress conditions

For accurate and reliable gene expression analysis using quantitative real-time reverse transcription PCR (qPCR), the selection of appropriate reference genes as an internal control for normalization is crucial. We hypothesized that non-coding, small nucleolar RNAs (snoRNAs) would be stably expressed in different barley varieties and under different experimental treatments, in different tissues and at different developmental stages of plant growth and therefore might prove to be suitable reference genes for expression analysis of both microRNAs (miRNAs) and mRNAs. In this study, we examined the expression stability of ten candidate reference genes in six barley genotypes under five experimental stresses, drought, fungal infection, boron toxicity, nutrient deficiency and salinity. We compared four commonly used housekeeping genes; Actin (ACT), alpha-Tubulin (α-TUB), Glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ADP-ribosylation factor 1-like protein (ADP), four snoRNAs; (U18, U61, snoR14 and snoR23) and two microRNAs (miR168, miR159) as candidate reference genes. We found that ADP, snoR14 and snoR23 were ranked as the best of these candidates across diverse samples. For accurate and reliable gene expression analysis using quantitative real-time reverse transcription PCR (qPCR), the selection of appropriate reference genes as an internal control for normalization is crucial. We hypothesized that non-coding, small nucleolar RNAs (snoRNAs) would be stably expressed in different barley varieties and under different experimental treatments, in different tissues and at different developmental stages of plant growth and therefore might prove to be suitable reference genes for expression analysis of both microRNAs (miRNAs) and mRNAs. In this study, we examined the expression stability of ten candidate reference genes in six barley genotypes under five experimental stresses, drought, fungal infection, boron toxicity, nutrient deficiency and salinity. We compared four commonly used housekeeping genes; Actin (ACT), alpha-Tubulin (α-TUB), Glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ADP-ribosylation factor 1-like protein (ADP), four snoRNAs; (U18, U61, snoR14 and snoR23) and two microRNAs (miR168, miR159) as candidate reference genes. We found that ADP, snoR14 and snoR23 were ranked as the best of these candidates across diverse samples. Additionally, we found that miR168 was a suitable reference gene for expression analysis in barley. Finally, we validated the performance of our stable and unstable candidate reference genes for both mRNA and miRNA qPCR data normalization under different stress conditions and demonstrated the superiority of the stable candidates. Our data demonstrate the suitability of barley snoRNAs and miRNAs as potential reference genes for miRNA and mRNA qPCR data normalization under different stress treatments.

High-Throughput Sequencing Reveals Diverse Sets of Conserved, Nonconserved, and Species-Specific miRNAs in Jute

MicroRNAs play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. A group of plant microRNAs are evolutionarily conserved; however, others are expressed in a species-specific manner. Jute is an agroeconomically important fibre crop; nonetheless, no practical information is available for microRNAs in jute to date. In this study, Illumina sequencing revealed a total of 227 known microRNAs and 17 potential novel microRNA candidates in jute, of which 164 belong to 23 conserved families and the remaining 63 belong to 58 nonconserved families. Among a total of 81 identified microRNA families, 116 potential target genes were predicted for 39 families and 11 targets were predicted for 4 among the 17 identified novel microRNAs. For understanding better the functions of microRNAs, target genes were analyzed by Gene Ontology and their pathways illustrated by KEGG pathway analyses. The presence of microRNAs identified in jute was validated by stem-loop RT-PCR followed by end point PCR and qPCR for randomly selected 20 known and novel microRNAs. This study exhaustively identifies microRNAs and their target genes in jute which will ultimately pave the way for understanding their role in this crop and other crops.

To bloom or not to bloom: role of microRNAs in plant flowering

During the course of their life cycles, plants undergo various morphological and physiological changes underlying juvenile-to-adult and adult-to-flowering phase transitions. To flower or not to flower is a key step of plasticity of a plant toward the start of its new life cycle. In addition to the previously revealed intrinsic genetic programs, exogenous cues, and endogenous cues, a class of small non-coding RNAs, microRNAs (miRNAs), plays a key role in plants making the decision to flower by integrating into the known flowering pathways. This review highlights the age-dependent flowering pathway with a focus on a number of timing miRNAs in determining such a key process. The contributions of other miRNAs which exist mainly outside the age pathway are also discussed. Approaches to study the flowering-determining miRNAs, their interactions, and applications are presented.

MicroRNAs and drought responses in sugarcane

There is a growing demand for renewable energy, and sugarcane is a promising bioenergy crop. In Brazil, the largest sugarcane producer in the world, sugarcane plantations are expanding into areas where severe droughts are common. Recent evidence has highlighted the role of miRNAs in regulating drought responses in several species, including sugarcane. This review summarizes the data from miRNA expression profiles observed in a wide array of experimental conditions using different sugarcane cultivars that differ in their tolerance to drought. We uncovered a complex regulation of sugarcane miRNAs in response to drought and discussed these data with the miRNA profiles observed in other plant species. The predicted miRNA targets revealed different transcription factors, proteins involved in tolerance to oxidative stress, cell modification, as well as hormone signaling. Some of these proteins might regulate sugarcane responses to drought, such as reduction of internode growth and shoot branching and increased leaf senescence. A better understanding on the regulatory network from miRNAs and their targets under drought stress has a great potential to contribute to sugarcane improvement, either as molecular markers as well as by using biotechnological approaches.

Genome-Wide Characterization of miRNAs Involved in N Gene-Mediated Immunity in Response to Tobacco Mosaic Virus in Nicotiana benthamiana

microRNAs (miRNAs) are a class of endogenous small RNAs (sRNAs) that play pivotal roles in plant development, abiotic stress response, and pathogen response. miRNAs have been extensively studied in plants, but rarely in Nicotiana benthamiana, despite its wide use in plant virology studies, particularly for studying N protein-tobacco mosaic virus (TMV) interactions. We report an efficient method using high-throughput sequencing and bioinformatics to identify genome-wide miRNAs in N. benthamiana. A total of 30 conserved miRNA families and 113 novel miRNAs belonging to 93 families were identified. Some miRNAs were clustered on chromosomes, and some were embedded in host gene introns. The predicted miRNA targets were involved in diverse biological processes, such as metabolism, signaling, and responses to stimuli. miRNA expression profiling revealed that most of them were differentially expressed during N-mediated immunity to TMV. This study provides a framework for further analysis of miRNA functions in plant immunity.

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 cold-inducible microRNAs in grapevine

Low temperature is one of the most important environmental factors that limits the geographical distribution and productivity of grapevine. However, the molecular mechanisms on how grapevine responds to cold stress remains to be elucidated. MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play an essential role during plant development and stress responses. Although miRNAs and their targets have been identified in several Vitis species, their participation during cold accumulation in grapevine remains unknown. In this study, two small RNA libraries were generated from micropropagated 'Muscat Hamburg' (V. vinifera) plantlets under normal and low temperatures (4°C). A total of 163 known miRNAs and 67 putative novel miRNAs were detected from two small RNA libraries by Solexa sequencing. Forty-four cold-inducible miRNAs were identified through differentially expressed miRNAs (DEMs) analysis; among which, 13 belonged to upregulated DEMs while 31 belonged downregulated DEMs. The expression patterns of the 13 DEMs were verified by real-time RT-PCR analysis. The prediction of the target genes for DEMs indicated that miRNA may regulate transcription factors, including AP2, SBP, MYB, bHLH, GRAS, and bZIP under cold stress. The 5'-RLM RACE were conducted to verify the cleavage site of predicted targets. Seven predicted target genes for four known and three novel vvi-miRNAs showed specific cleavage sites corresponding to their miRNA complementary sequences. The expression pattern of these seven target genes revealed negative correlation with the expression level of the corresponding vvi-miRNAs. Our results indicated that a diverse set of miRNAs in V. vinifera are cold-inducible and may play an important role in cold stress response.

Orthologous plant microRNAs: microregulators with great potential for improving stress tolerance in plants

Small RNAs that are highly conserved across many plant species are involved in stress responses. Plants are exposed to many types of unfavorable conditions during their life cycle that result in some degree of stress. Recent studies on microRNAs (miRNAs) have highlighted their great potential as regulators of stress tolerance in plants. One of the possible ways in which plants counter environmental stresses is by altering their gene expression by the action of miRNAs. miRNAs regulate the expression of target genes by hybridizing to their nascent reverse complementary sequences marking them for cleavage in the nucleus or translational repression in the cytoplasm. Some miRNAs have been reported to be key regulators in biotic as well as abiotic stress responses across many species. The present review highlights some of the regulatory roles of orthologous plant miRNAs in response to various types of stress conditions.

The high-throughput sequencing of small RNAs profiling in wide hybridisation and allopolyploidisation between Brassica rapa and Brassica nigra

Small RNAs play an important role in maintaining the genome reconstruction and stability in the plant. However, little is known regarding the role of small RNAs during the process of wide hybridisation and chromosome doubling. Therefore, the changes in the small RNAs were assessed during the formation of an allodiploid (genome: AB) and its allotetraploid (genome: AABB) between Brassica rapa (♀) and Brassica nigra (♂) in the present study. Here, the experimental methods described in details, RNA-seq data (available at Gene Expression Omnibus database under GSE61872) and analysis published by Ghani et al. [1]. The study showed that small RNAs play an important role in maintaining the genome stability, and regulate gene expression which induces the phenotype variation in the formation of an allotetraploid. This may play an important role in the occurrence of heterosis in the allotetraploid.

Discovery and analysis of microRNAs in Leymus chinensis under saline-alkali and drought stress using high-throughput sequencing

Leymus chinensis (Trin.) Tzvel. is a perennial rhizome grass of the Poaceae (also called Gramineae) family, which adapts well to drought, saline and alkaline conditions. However, little is known about the stress tolerance of L. chinensis at the molecular level. microRNAs (miRNAs) are known to play critical roles in nutrient homeostasis, developmental processes, pathogen responses, and abiotic stress in plants. In this study, we used Solexa sequencing technology to generate high-quality small RNA data from three L. chinensis groups: a control group, a saline-alkaline stress group (100 mM NaCl and 200 mM NaHCO3), and a drought stress group (20% polyethylene glycol 2000). From these data we identified 132 known miRNAs and 16 novel miRNAs candidates. For these miRNAs we also identified target genes that encode a broad range of proteins that may be correlated with abiotic stress regulation. This is the first study to demonstrate differentially expressed miRNAs in L. chinensis under saline-alkali and drought stress. These findings may help explain the saline-alkaline and drought stress responses in L. chinensis.