An integrative approach to identify hexaploid wheat miRNAome associated with development and tolerance to abiotic stress
- Zahra Agharbaoui1,
- Mickael Leclercq†2, 3,
- Mohamed Amine Remita†2,
- Mohamed A Badawi1,
- Etienne Lord2,
- Mario Houde1,
- Jean Danyluk1,
- Abdoulaye Baniré Diallo†2Email author and
- Fathey Sarhan†1Email author
© Agharbaoui et al.; licensee BioMed Central. 2015
Received: 23 January 2015
Accepted: 26 March 2015
Published: 24 April 2015
Wheat is a major staple crop with broad adaptability to a wide range of environmental conditions. This adaptability involves several stress and developmentally responsive genes, in which microRNAs (miRNAs) have emerged as important regulatory factors. However, the currently used approaches to identify miRNAs in this polyploid complex system focus on conserved and highly expressed miRNAs avoiding regularly those that are often lineage-specific, condition-specific, or appeared recently in evolution. In addition, many environmental and biological factors affecting miRNA expression were not yet considered, resulting still in an incomplete repertoire of wheat miRNAs.
We developed a conservation-independent technique based on an integrative approach that combines machine learning, bioinformatic tools, biological insights of known miRNA expression profiles and universal criteria of plant miRNAs to identify miRNAs with more confidence. The developed pipeline can potentially identify novel wheat miRNAs that share features common to several species or that are species specific or clade specific. It allowed the discovery of 199 miRNA candidates associated with different abiotic stresses and development stages. We also highlight from the raw data 267 miRNAs conserved with 43 miRBase families. The predicted miRNAs are highly associated with abiotic stress responses, tolerance and development. GO enrichment analysis showed that they may play biological and physiological roles associated with cold, salt and aluminum (Al) through auxin signaling pathways, regulation of gene expression, ubiquitination, transport, carbohydrates, gibberellins, lipid, glutathione and secondary metabolism, photosynthesis, as well as floral transition and flowering.
This approach provides a broad repertoire of hexaploid wheat miRNAs associated with abiotic stress responses, tolerance and development. These valuable resources of expressed wheat miRNAs will help in elucidating the regulatory mechanisms involved in freezing and Al responses and tolerance mechanisms as well as for development and flowering. In the long term, it may help in breeding stress tolerant plants.
Abiotic stresses such as cold, drought, salt and aluminum (Al) limit plant growth and development, causing reduction in crop yield and important economic losses for farmers. To tolerate these stresses, plants have evolved a broad spectrum of metabolic, physiological and developmentally adaptations. These adaptive changes are under the control of dynamic networks of genetic regulatory mechanisms that involve a large number of stress responsive genes. MicroRNAs (miRNAs), a major class of small non-coding RNAs, have emerged as key regulators of gene expression at the post-transcriptional level during plant growth and development [1-3]. Several studies have shown that many miRNA families are involved in response to different abiotic stresses in many species [4-7]. A large number of plant miRNAs and their targets have been identified in the plant model Arabidopsis thaliana and many other species. Recent results have shown that plant miRNA genes are dispersed throughout the genome  within protein coding genes [8,9], introns of protein coding and non-coding genes, and in intergenic regions [10,11]. Moreover, miRNAs may be produced from repetitive transposable elements [12,13]. To date, at the best of our knowledge, 2707 wheat miRNA candidates were identified by both bioinformatics and experimental approaches, using wheat expressed sequence tags (EST) database, the available genomic sequences of the hexaploid wheat genome, its individual chromosome arms and its ancestors [6,13-29]. Among the wheat miRNA published sequences, 237 are registered in miRBase, a database of experimental miRNAs , and 170 are registered in PMRD, a database of plant miRNAs identified using an in silico approach . Although the wheat genome is completely sequenced, it is not yet possible to perform a thorough genome-wide study in the hexaploid wheat T. aestivum since the genome is not completely assembled and annotated. This is caused by its large and complex genome containing a high percentage of DNA repeats (hexaploid genome AABBDD with approximately 1.7 × 1010 bp with at least 80% of DNA repeats) . In silico approaches for the prediction of miRNAs include screening genomic or EST databases for orthologous sequences of known miRNAs and analyzing their pre-miRNA hairpin structures. Although these approaches were successful in identifying conserved miRNAs in plants that have their genomes fully sequenced and annotated [10,33,34], they eliminate the potential of searching for low abundance miRNAs that are often lineage-specific  or condition-specific  or that appeared recently in evolution (young miRNAs). The challenge is bigger using polyploid species with partially sequenced and assembled genome such as the hexaploid wheat having a high content of repetitive DNA. To tackle this issue, one should develop conservation-independent techniques based on structure analyses and/or expression pattern of dicer cleavage products among pre-miRNAs .
Most computational approaches labeled as miRNA predictors are actually pre-miRNA predictors, in the sense that they identify candidate genomic regions that may form pre-miRNAs but rarely take into account the availability of candidate mature miRNA evidence within the pre-miRNA. Several tools such as miRDeep [37,38], miRanalyzer [39,40] and MiRdup  were developed to predict miRNAs from raw reads data and shown to be accurate in most cases. Furthermore, many factors that affect miRNA expression including genotypes, tissues, age, development stage, growth condition (soil, hydroponic solution, temperature, humidity and photoperiod), stress treatment, are rarely considered in previous wheat miRNA identification studies. All wheat reported miRNAs were identified in libraries produced from seedlings or plants grown under normal conditions [14,21,23,26,31], or tissue exposed to heat  or seedling  and pollen mother cells from plants  exposed to cold stress , or drought . They were identified from different genotypes of winter or spring wheat in soil, or hydroponic solution and under different photoperiod conditions, or in field conditions. Since miRNA expression is tissue specific and regulated in response to plant development and growth conditions, the miRNA repertoire of hexaploid wheat is still incomplete. Although a large number of miRNAs associated with development or some abiotic stresses in wheat were previously identified, their functional diversity in Al, freezing tolerance, and floral transition in winter wheat is still unknown. Hence, the identification of miRNAs associated with tolerance to abiotic stress and floral transition is a first step towards the elucidation of their role in wheat.
Identification of miRNA candidates and their targets in hexaploid wheat
Our miRNA discovery pipeline consists of more than twenty steps divided in three main parts: producing and sequencing small RNAs, predicting miRNAs from deep sequencing data, classifying predicted miRNAs based on their expression profiles and Gene Ontology (GO) of their target genes (Figure 1). The sequencing of ten constructed libraries from three wheat genotypes grown under different abiotic stress conditions and development stages (Additional file 1: Method S1) yielded a total of 89,105,096 redundant raw reads (66,400,401 distinct reads) with 56% of high sequence quality (Additional file 2: Table S2). Before mapping the raw reads, we collected 1.4 million wheat ESTs from several databases, clustered into 127,039 Uniref clusters yielding to the best of our knowledge, the largest well-annotated EST databank in wheat (Additional file 1: Method S2). After raw reads adapter removal, a total of 56.4 million unique raw reads were mapped to our collected EST database resulting of 5.4 million unique mapped sequences (Figure 1). We identified 168,834 small RNAs and extracted 337,668 potential pre-miRNAs. Among the extracted pre-miRNA, 17,180 and 39,144 potential hairpins satisfy the minimal structural criteria of miRNAs assessed by two well-known pre-miRNA predictors, MiPred  and HHMMiR , respectively (Figure 1). The lack of consideration of the mature miRNA localization in the hairpin, the expression profiles of reads throughout the hairpin, and/or other evidence of plant-miRNA characteristics represent major causes of overestimation of the number of candidate hairpins. We used a two-step method of miRNA identification within its hairpin from the folded pre-miRNA using MiRdup* and expression profile filtering using Meyers et al. relaxed criteria  (Figure 1). In the first step, MiRdup* trained on experimentally validated miRNAs from different datasets (all miRBase species, all plants and monocots only) localizes miRNA positions in pre-miRNAs (Additional file 1: Method S4 and Additional file 2: Table S3). The use of MiRdup* reduced the number of pre-miRNA hairpin candidates by 81% (Additional file 3: Figures S1 and S2). In the second step, the expression profile filtering was based on miRNA expression pattern using the abundance of the miRNA candidates in each library and the distribution of reads mapped to a candidate pre-miRNA according to Meyers et al. relaxed criteria . This allowed the elimination of the miRNA dicer-Like candidates and the reduction of miRNA candidates by 84%.
Overall, this method results in the identification of more candidates (Additional file 3: Figures S1 and S2) compared to Meyers et al.,  and MIRcheck . Consequently, it yields pre-miRNA candidates that have various ranges of secondary structures as shown in dotbracket notation (Additional file 4: data SD1). Taken together, our approach identifies 199 unique miRNA candidates associated with 361 pre-miRNAs (Additional file 4: data SD1). It is important to notice that the majority of reads (95%) and the predicted miRNAs (64%) have the highest quality value for their sequence (Additional file 2: Table S4). It is important to note that, MiRdup* captures 95% of the miRNAs identified using MIRcheck with the same criteria from Meyers et al.  while 151 putative candidates (containing validated candidates) are excluded by MIRcheck (Additional file 3: Figures S2). In addition, our pipeline identified miRNAs that have features that are species specific, clade specific or shared between several species. We found that among the 199 identified miRNAs, 147 were identified by MiRdup* trained on all species of miRBase; only 49 were commonly identified by MiRdup* trained on the three datasets (all miRBase, all plants, monocot only). This suggests that these miRNAs share common features with all the widely separated plant lineages recorded in the database miRBase. For instance, apMir_22246 corresponding to miR160 with perfect match in wheat and moss Physcomitrella patens (ppt-miR160) is highly expressed in our investigated conditions indicating that this miRNA may play common biological functions in plants kingdom. While 109 miRNAs were only identified by MiRdup* trained on all plant and 92 miRNAs when trained on only monocot.
The number of identified miRNAs may be an overestimation due to the redundancy created by similar but not identical ESTs in part due to the polyploid nature of wheat. In the latter scenario, two or more closely related ESTs (true homeologs or ESTs with SNP differences) could encode identical or closely related miRNAs. Furthermore expressed isomiRNAs that share the same properties with the real miRNA in one library could be the dominant functional in another library. The Additional file 3: Figure S3a highlights that the majority (about 69%) of the predicted miRNAs are associated with one pre-miRNA. Furthermore, most of the pre-miRNA candidates (93%) harbour a unique miRNA leading to an exclusive miRNA/pre-miRNA association (Additional file 3: Figure S3b). To characterize further the nature of the pre-miRNA candidates, we determined if they were associated with repetitive transposable elements and protein coding regions. The results revealed that 20% of pre-miRNAs that correspond to 6.5% of miRNA candidates overlap with transposable elements (at e-value of 5E-5 with 80% identity) from TREP database (Additional file 2: Table S5a). In addition, 15% of ESTs corresponding to less than 5% of miRNA candidates overlap partly with protein-coding regions (at e-value of 1E-20 with 75% identity) from protein plant database (Additional file 2: Table S5b).
Prediction of miRNA targets is an important step to elucidate miRNA function in regulating gene expression. Among the identified candidates, 67% (133/199) were predicted to have Uniref annotated target genes (Additional file 2: Table S6). Unlike animal target genes, it is generally accepted that plant targets adopt a perfect seed match with the corresponding miRNAs, allowing more accuracy in their prediction. We found that 37 miRNA candidates are predicted to target a unique gene identified as UniRef (Additional file 3: Figure S3c). Although the majority of miRNA candidates seem to have more than two targets, detailed analysis reveals that in many cases the targets annotated with different UniRefs have the same gene description (Additional file 2: Table S7). To better explore the functional properties of the target genes, GO analyses were performed (Additional file 3: Figure S3d).
Selected GO Slim enrichment in the different libraries and their relevant target genes
GO Slim description
Relevant associated targets
CBFIVb-B20, MIKC-type MADS-box transcription factor, Auxin-responsive protein, Homeobox-leucine zipper protein, Histones, E3 ubiquitin-protein ligase
DEAD-box ATP-dependent RNA helicase, Homeobox-leucine zipper proteins, Histones
Gibberellin 20 oxidase , Lipoxygenase, Glutathione peroxidase, Superoxide dismutase, Fructose-bisphosphate aldolase, Alcohol dehydrogenase, NADP-malic enzyme
Caffeic acid 3-O-methyltransferase, Tricetin 3',4',5'-O-trimethyltransferase, Serine/threonine-protein kinase, Glutathione-S-transferase
Non-specific lipid-transfer protein, Homeobox-leucine zipper proteins
Histones, WCOR719, DNA-directed RNA polymerases, MIKC-type MADS-box transcription factors, Serine/threonine-protein kinase, Cullin-1
Sodium/hydrogen exchanger, Triose phosphate translocator, Protein transport protein, Non-specific lipid-transfer protein, Sucrose transporters
Resp. to endogenous stimulus
Auxin-responsive protein, Homeobox-leucine zipper protein, Myo-inositol 1-phosphate synthase
Multicellular organismal dev.
Ubiquitin-like-specific protease, MIKC-type MADS-box transcription factors, Homeobox-leucine zipper protein
DL related protein, Ubiquitin-like-specific protease
DNA metabolic process
Replication factor C, Histones, BARE-1 polyprotein
Sec. metabolic process
Phenylalanine ammonia-lyase, Tricetin 3',4',5'-O-trimethyltransferase
Characteristics of the miRNA candidates
The diversity of predicted miRNA sequences is greater at 19 nt in length in all libraries (Additional file 3: Figure S7a-d) while the diversity of conserved miRNAs is greater at 21 nt (Additional file 3: Figure S7e) (unique sequences). For redundant miRNA sequences, a major peak at 21 nt was observed for both predicted and conserved miRNAs (Additional file 3: Figure S7a-e). In addition, the majority of the identified miRNA expressed in the different investigated conditions are 19 or 21 nt long depending on the tissues, stresses, growth conditions or genotypes (Figure 2d). These miRNAs were shown to regulate at least 150 targets (60 unirefs) and at most 900 targets (335 unirefs) in all the explored conditions (Figure 2e).
Confirmation of predicted miRNA candidates
Characteristics of selected miRNAs using MiRdup* and MIRcheck validated by northern blot
Relevant targets (with uniref ID)
Predicted MiRNAs in common with MiRdup* and MIRcheck
Salt, floral transition and flowering
Glutathione peroxidase (Q9SME6), putative phosphatase phospho1 (M8CZ66)
Al, floral transition and flowering
Putative membrane-associated protein (gi|22548307|gb|BU100508.1|BU100508)
Cold, Al and salt
Putative RSH disease resistance-related protein (Q8H5X7), T-complex protein 1 subunit alpha (I3RZC6)
Cold, floral transition and flowering
Predicted MiRNAs specific to MiRdup*
Not differentially expressed
Dehydrin, (CD909074, TA50415_4565), Phosphorylase (Q84P16)
Cold, floral transition and Al
Salt, Al and floral transition
Conserved miRNAs predicted in common by MiRdup* and MIRcheck
Cold, salt, Al, floral transition and flowering
Auxin response factor (M8BC98), Auxin responsive protein (R7WEP7)
Cold, salt, Al, flowering
Bifunctional 3′-phosphoadenosine 5′-phosphosulfate synthase (M7ZFX2), ATP sulfurylase (M9T1P9)
Expression of the identified miRNAs in response to different abiotic stresses and plant development in wheat
Functional classification of abiotic stress and developmentally regulated miRNAs in wheat
The miRNA group associated with cold responses (Co5) is specifically enriched for membrane (triose phosphate translocator) in the category cell component (Additional file 3: Figure S4b) and signal transduction (Auxin-responsive proteins) in the category biological process (Figure 5 and Additional file 3: Figure S6b). Consistently, enrichment is also found for nucleus (Additional file 3: Figure S4b), protein binding activity (Additional file 3: Figure S5b), nucleobase containing component metabolic process and response to endogenous stimulus (Figure 5 and Additional file 3: Figure S6b) which are all overrepresented by auxin responsive proteins. These results indicate that cold regulated miRNAs may function in carbon partitioning during photosynthesis and in auxin-activated signaling pathways. For miRNA groups associated with Al responses (Al5-Al6), an enrichment is found for hydrolase activity (protein phosphatase 2C, lipase), catalytic activity (glutathione peroxidase, phenylalanine ammonia-lyase) (Additional file 3: Figure S5b), response to endogenous stimulus (Auxin Response Factors), DNA metabolic process (histone 4) (Figure 5 and Additional file 3: Figure S6b). These results indicate that Al-regulated miRNAs may function in regulation of gene expression and signaling as well as plant defense under oxidative stress. More interestingly, many targets found for miRNA groups associated with cold (Co1, Co2, Co4) and Al (Al1, Al2 and Al3) tolerance are known for their function in stress adaptation. For miRNA groups associated with cold tolerance, the groups Co1 and Co2, showed enrichment for the GO Slim term response to stress (phosphoglycerate mutase, Defensin-like protein 1, Universal stress protein A-like protein). In addition, Co1 is enriched for response to abiotic stimulus (thaumatin-like protein, glutathione S-transferase) (Figure 5 and Additional file 3: Figure S6b) and the group Co2 is enriched for cell wall (Defensin-like protein, phospho-3-sulfolactate synthase-like), nucleus (CBFIVb-B20), hydrolase activity (Ubiquitin-like-specific protease, Serine carboxypeptidase) and catalytic activity (Gibberellin 20 oxidase, Glyceraldehyde-3-phosphate dehydrogenase) (Additional file 3: Figure S5b). This indicates that the identified cold regulated miRNAs may function in proteolysis, gibberellins biosynthesis and glucose metabolism. The group Co3 is enriched for transporter activity (Sodium/hydrogen exchanger) (Additional file 3: Figure S5b) while the group Co4 is enriched for pollen-pistil interaction (Serine/threonine-protein kinase) (Figure 5 and Additional file 3: Figure S6b) indicating that they may also have a function in ion transport and signaling.
Moreover, for groups associated with Al tolerance, the most enriched GO Slim terms are mitochondrion, transporter activity (Cytochrome b-c1 complex) and carbohydrates metabolic process (glyceraldehyde-3-phosphate dehydrogenase) indicating that miRNAs from these groups may mediate glycolysis and respiration process under Al stress conditions. Enrichment is also found for the terms catalytic activity, hydrolase activity, sequence specific transcription factor activity (Additional file 3: Figure S5b), secondary metabolic process, embryo development, anatomical structural morphogenesis and multi-organismal development (Figure 5 and Additional file 3: Figure S6b). They are overrepresented by many important metabolism enzymes involved in phenylpropanoid metabolism and ubiquitination process (phenylalanine ammonia-lyase, DEAD-box ATP-dependent RNA helicase, Ubiquitin carboxyl-terminal hydrolase) as well as transcription factors (Homeobox-leucine zipper proteins).
Interestingly, the targets in miRNA groups associated with development are involved in cell growth and flowering. The miRNA groups associated with flowering (Dev1, Dev2) show enrichment for the GO Slim terms protein binding activity, sequence specific DNA transcription binding activity (MIKC-type MADS-box transcription factors), and/or kinase activity (Serine/threonine-protein kinase) (Additional file 3: Figure S5b). They are also enriched for nucleobase containing compound metabolism and pollen pistil interaction that are specifically represented by two well characterized regulators genes (MIKC-type MADS-box and Serine/threonine-protein kinase) (Figure 5). One miRNA group associated with floral transition (Dev 4) is enriched for the GO Slim terms transport and response to endogenous stimulus, specifically represented by Auxin-responsive proteins. Furthermore, the potential function of ubiquitously expressed miRNA libraries was also investigated. They show significant enrichment for the GO Slim terms thylakoid (Additional file 3: Figure S4b), kinase activity, nucleotide binding activity (Additional file 3: Figure S5b), and cell protein modifications (Figure 5 and Additional file 3: Figure S6b). This result indicates that constitutively expressed miRNAs may modulate the basic cellular functions reflecting their vital regulatory role in other growth conditions yet to be identified in wheat.
Discussion and conclusions
The wheat miRNA pipeline
In this study, we developed a pipeline that identifies conserved as well as clade and species-specific or young miRNAs. This pipeline can be easily adapted for other plant species. To predict miRNAs from NGS and analyze their function, the steps described in Figure 1 are required. While several steps are standard in NGS analyses [4,49,53], we improved the miRNA prediction steps by integrating folded pre-miRNA candidates, expression profiling and functional analyses of differentially expressed candidates. To address the step of miRNA prediction, we decided to exploit two methods with different algorithmic schemes MiPred  and HHMMiR  to have a broad range of hairpin candidates. These methods were trained on pre-miRNAs from plants and wheat sequences available in miRBase and resulted in the identification of a large number of pre-miRNA candidates using the predictors MiPred  and HHMMiR . To address issues of latter methods for the lack of consideration of mature miRNA and their surrounding biological features, we developed a classifier that ranked the best 35 biological features of plant miRNAs that was integrated into MiRdup* (Additional file 2: Table S3). For robustness, the classifier’s models were trained separately on three datasets (all miRBase species, all plants and only monocots). This increases species-specificity and allows the discovery of features that distinguish wheat miRNAs from those of other species. The developed classifier (MiRdup*) was able to reduce the level of false prediction obtained by MiPred  and HHMMiR  by more than 81% (Figure 1, Additional file 3: Figures S1 and S2) and allowed the assessment of the position of a miRNA in a given pre-miRNA sequence. In addition, the combinatorial analysis between MiRdup* and MIRcheck  which identifies 20-nt regions of a given plant pre-miRNA using a predetermined set of rules and constraints, show that MIRcheck is too stringent and easily removed experimentally validated miRNAs (Figure 3b and Additional file 3: Figures S2).
The availability of wheat EST databases and our approach enabled us to identify with confidence 199 miRNA candidates. These candidates may include miRNA gene homeologs from the three genomes of hexaploid wheat, or ESTs with SNP differences in different wheat varieties. It is also reasonable to assume that these families represent only a fraction of the total miRNAs that may exists in hexaploid wheat since many small RNAs still remain unmapped to wheat sequences or conserved miRNAs from miRBase. The availability of the complete assembled and well-annotated hexaploid wheat genome will help to complete the discovery of the remaining miRNAs.
It is important to emphasize that among the predicted miRNAs, in spite of being derived from ESTs, less than 5% of the mature miRNAs are associated with known protein coding regions and less than 7% are related to transposable elements (Additional file 2: Table S5a and S5b). According to Dinger et al., , many transcripts are categorized as bifunctional RNAs. They can be translated into protein but also function independently as RNA. The presence of such bifunctional RNAs challenges the assumption that the RNA world can be neatly parsed between mutually exclusive protein-coding and non-coding categories.
MiRNA candidates associated with abiotic stress responses
This study represents one of the largest de novo miRNAome analyses in response to different abiotic stresses and development in hexaploid wheat. Although many cold responsive miRNAs have been identified in spring wheat using NGS , our study identified a large number of novel candidates regulated by cold, vernalization, Al and salt with dynamic and complex expression patterns (Figures 3b, 4b and Additional file 2: Table S9). Several identified miRNAs are either associated with a specific stress or common to at least two stresses (Figures 3b and 4c). Many of their targets are known to be stress-related genes (Figure 5 and Additional file 3: Figures S4b-6b) commonly regulated under abiotic stresses.
Our results show that miRNAs may mediate plant responses to Al treatment by regulating expression of stress related genes particularly those involved in auxin signaling and fatty acid metabolism. This is consistent with the fact that Al affects the relative abundance of membrane lipids and the degree of fatty acid unsaturation [55,56] and Auxin Response Factors (ARFs) that are known to inhibit root development in response to Al toxicity . In addition, the experimentally validated apMir_22246 (which corresponds to tae-miR160) is regulated by Al exposure (Figure 3c) and targets specifically ARFs. Many ARF members are known to be regulated by miR167 and miR160 and to play regulatory roles in adventitious rooting , supporting the possible role of apMir_22246 in root development under Al treatment.
MiRNA candidates associated with cold responses and freezing tolerance
Our data indicate that cold regulates the expression of several miRNAs in spring as well as in winter wheat (Figures 3b, 4b and c). Four miRNA groups associated with cold tolerance (Additional file 2: Table S10) target a set of cold regulated genes known to be involved in freezing tolerance including the transcription factors CBFs, dehydrins, DEAD-box RNA helicases, thaumatin-like protein [59-62]. Interestingly, many candidate miRNA target genes related to the ICE1–CBF major pathway that regulates freezing tolerance in cold hardy plants. This includes the targets DEAD-box ATP-dependent RNA helicase 12, CBF and dehydrin (Additional file 2: Table S7). Results from our previous studies demonstrated that genes related to the ICE1–CBF pathway play a critical role in freezing tolerance in hexaploid wheat . Here we show that the miRNA candidate apMir_16808 is regulated in response to cold (Figure 3b), and target the cold responsive genes dehydrins (Table 2) [59,62]. The candidate apMir_19532 from miRNA group associated with cold tolerance target CBFIVb-B20 gene (Additional file 2: Table S7). These results suggest that these miRNAs may contribute to freezing tolerance by regulating cold-regulated genes belonging to the CBF regulon in winter wheat.
Predicted miRNA target genes common in regulating several stresses
Plants evolved common regulatory mechanisms to adapt to environmental stresses such as oxidative stress commonly induced by both cold and Al. Our results show that many of the identified abiotic stress responsive miRNAs exhibited a common stress expression pattern (Figures 3b and c and 4c). For instance, the expression of the new member of miR395 family, miR395-21 corresponding to apMir_20968, is commonly regulated in response to cold and Al stress (Figure 3c) indicating that miR395 is not specific to sulfate starvation as previously reported in Arabidopsis and rice [49,64]. Zhao et al.,  also reported that miR395 is involved in phosphate homeostasis in wheat. This indicates that miR395 mediates not only plant response to sulfate deficiency but also may mediate responses to other nutrients that are imbalanced under abiotic stress conditions. Taken together, our results indicate that miR395 would play a common role in plant nutrient homeostasis under abiotic stress conditions. In agreement with previous suggestions, our results indicate that miRNAs coordinate crosstalk among different nutrient deficiencies. This is the first indication that crosstalk between cold, Al stress and plant nutrients could be regulated by miRNAs. Moreover, we show that the miRNA candidate apMir_20602 is also commonly regulated under cold, salt and Al (Figure 4b) and targets glutathione peroxidase (Table 2). Recent findings showed that human miRNAs regulate glutathione peroxidase expression to maintain redox homeostasis . This supports the possible role of apMir_20602 in mediating crosstalk between abiotic stress responses by regulating glutathione metabolism.
Wheat vernalization responsive miRNAs associated with floral transition and flowering
In this study, we investigate the role of miRNAs during the transition from the vegetative to the reproductive phase, and during flowering in winter wheat that requires vernalization to flower. We found that among developmentally responsive miRNAs, many candidates target cold responsive genes known for their function in flowering transition and flower development (Additional file 2: Table S7). For instance, the candidate apMir_19892 corresponding to hvu-miR444b (Additional file 4: data SD1) could target many MIKC-type MADS-box transcription factors, the homologs of TaAGL17 and OsMADS57. In wheat, MIKC-type MADS-box transcription factors control flower development and morphogenesis . In barley, this target contains both the target site for miR444b and the precursor sequence for miR444a ). In rice, OsMADS57 is involved in axillary bud development and regulation of tillering through down-regulation of miR444a . Since the miRNA variants from miR444 family are functional, and MADS-box genes are collectively regulated by the miR444 family , we suggest that apMir_19892 may mediate flowering through the regulation of MIKC-type MADS-box transcription factor gene expression. ApMir_19532 target genes encoding Ubiquitin-like-specific protease ESD4 known to regulate plant responses to cold and the time of flower initiation [71,72]. In addition, apMir_20860 corresponding to miR159 (Additional file 4: data SD1) are involved in promotion of floral transition in many species. In ornamental plants, miR159-regulated GAMYB expression is an effective pathway of flowering time control . This suggests that apMir_19532 and apMir_20860 may mediate flowering time in wheat through the regulation of Ubiquitin-like-specific protease ESD4 and GAMYB gene expression.
Plant material and small RNAs isolation
In this study, three genotypes of hexaploid wheat (Triticum aestivum L. 2n = 6x = 42, AABBDD), one spring genotype (cv Bounty, cold and Al sensitive) and two winter genotypes (cv Clair, cold tolerant and Atlas66, Al tolerant genotype) were used to construct ten different small RNA libraries from plants in vegetative and/or reproductive stages and/or exposed to different stress treatments or under normal conditions (Additional file 1: Method S1 and Additional file 2: Table S1). To identify miRNAs that are associated with different development stages, tissues from both vegetative and reproductive phases were used. Vegetative phase samples include leaves and crown from the aerial part of plants. Reproductive phase samples include leaves at flag leaf stage, developing spikes with sizes ranging from 2 to 110 mm, and spikes partially and completely-opened with and without pollen. We used also root tips to identify miRNAs associated with Al and salt stress, and aerial parts including leaves and crown to identify miRNAs associated with cold and salt. In addition, we used different genotypes of winter (tolerant) and spring (sensitive) wheat to identify miRNAs associated with Al and freezing tolerance. Small RNA extraction was initiated from 200 mg of a mixture of leaves, stem or root tip tissues from 10 to 100 seedlings for each time point. Control and treated plants were sampled at the same time of the day for each time point (except for the first day where a few samples were taken at short time points) as described in (Additional file 1: Method S1 and Additional file 2: Table S1). Small RNAs (below 200 nt) were isolated from each sample using the mirVana miRNA Isolation Kit (Ambion Inc. US). MiRNAs (small RNAs below 40 nt) from each time point were isolated using the flashPAGE fractionation kit (Ambion) (Additional file 1: Method S1 and Additional file 2: Table S1), and then purified using the flashPAGE Reaction Clean-up kit (Ambion) according to the manufacturer’s protocols. Their integrity was assessed using a DNA 1000 LabChip on an Agilent 2100 Bioanalyzer (Santa Clara, CA, USA).
MiRNA libraries construction and sequencing
Twenty five nanograms of purified miRNAs from each time point of a given condition (Additional file 1: Method S1 and Additional file 2: Table S1) were pooled and used as a template to produce the corresponding miRNA library. MiRNAs were tagged with a barcode system containing ten unique and specific amplification primers (1 barcode/library) and ten cDNA libraries were produced using the SREK kit (small RNA expression Kit, Ambion) according to the manufacturer's protocol. The libraries were sequenced on the SOLiD Analyzer according to the standard protocol (V2.1 Applied Biosystems).
Experimental validation of predicted miRNAs
For each library, identical amounts of plant tissues from each time point were ground and mixed with TRIzol Reagent (Life Technologies). The same extract volume from each time point of each library was pooled to isolate small RNAs using the mirVana miRNA Isolation Kit. Five micrograms of small RNAs from each library were analyzed by northern blot . The experiment was repeated at least twice for each selected probe. The oligonucleotide probes are presented in Additional file 2: Table S11.
Identification and extraction of potential pre-miRNA candidates from sequenced small RNAs
From the 89 million reads obtained from the ten libraries, we first removed reads that have low quality scores as recommended for SOLID sequencing [39,75] (Additional file 2: Table S2). After adapter removal using the program cutadapt v0.9 , small RNAs between 18 and 30 nt were mapped to ESTs with the MAQ v07.1 program  allowing a maximum of two mismatches (Additional file 1: Method S3). Then, sequences with low complexity or containing repeats were filtered out using RepeatMasker v3.2.9  with RepBase15.09 and Repeatmasker-Libraries-20130422 , and a slow search method against Triticum aestivum and Oryza sativa. For each mapped EST, we considered two sequences that could include a potential pre-miRNA candidate as follows: 20 nt before the start and 160 nt after the end of the mapping were extracted from both EST strands.
The secondary structures of the extracted sequences (pri-miRNA) were folded with RNAfold from ViennaRNA v1.8.4 package  to identify those having a hairpin-like shape, one of the fundamental characteristics of pre-miRNAs. Then, these sequences were submitted to two pre-miRNA predictors using different algorithmic schemes, HHMMiR  and MiPred , which were trained with pre-miRNAs of all cloned or sequenced miRNAs from miRBase . Finally, to identify conserved miRNAs, we performed a blastn of small RNAs with more than 100 reads in at least one library against miRBase (V21)  with a word_size 7, maximum e-value 0.1, percentage identity 80, gap open penalty 5, gap extension penalty 2, match score 1, mismatch score −2, filter low complexity. We restricted the blast results as follows: query coverage of 90%, subject coverage of 90%, with no gap allowed.
Filtering false positive pre-miRNAs
We adapted our previously published machine learning classifier  to the best plant features associated with the position of miRNA-miRNA* duplex in the pre-miRNA (Additional file 1: Method S4), and we named this version MiRdup*. The latter differs from the original one on the 35 retained features (Additional file 2: Table S3) that are relevant for plants and has been trained on all experimental (cloned or sequenced) miRNAs from miRBase subdivided in three datasets (all species, all plants and only monocots). This classifier computes a score of prediction scaled between 0 and 100 (more evidence). The potential pre-miRNAs from MiPred or HHMMiR that obtained a MiRdup* classification score higher than 90 and with miRNA read abundance above 100 in at least one library were selected as candidates. In addition, to help identifying the potential functional miRNAs among several candidates, we applied the relaxed expression rules derived from the update of the specific criteria for plant miRNA annotation reported by Meyers et al., . In parallel, we applied MIRcheck , a well-known tool for plant miRNA identification, on the overall predicted miRNAs to compare the differences with MiRdup*. A combinatorial analysis between the two tools is provided (Additional file 3: Figures S1 and S2). Then, the pre-miRNAs were blasted against TREP database to identify miRNAs that originate from transposable elements (Additional file 2: Table S5a). The e-value threshold used is 5.0E-05 and Hit Coverage (HC) ≥ 85 and percentage of identity ≥80. We performed also a blastx of ESTs producing the identified pre-miRNAs against proteins coding sequences from protein plant database with default parameters (Additional file 2: Table S5b). The threshold to include a blast hit of an EST into a given protein core is an e-value lower than 1.00E-20 and query coverage or hit coverage higher than 85% and percentage of identity higher than 75%.
Statistical analyses of the abundance of potential miRNAs
To quantify and compare sequence abundance across different libraries, raw read counts were normalized using rpm (reads per million). Sequences with read counts lower than 100 in all libraries are removed. Significance level of the difference of small RNA between two libraries was analyzed using a corrected Z–Score method as described in Kal et al., . An adjustment for multiple comparisons based on the false discovery rate (FDR)  was performed (FDR < 5%). The small RNAs with fold change lower than 0.5 or higher than 2.0 were retained.
MiRNA target analyses and GO enrichments
MiRNA target genes were identified using the FASTA engine of Tapir v1.0 program, with a stringent maximum score of three and minimum free energy ≥ 0.7  excluding ESTs annotated as unknown proteins. For the obtained and annotated target genes, we retrieved their classification in Gene Ontology (GO) through the GO Slim viewer on AgBase webserver . The GO Slim enrichments were performed using the standard hypergeometric test. The wheat genome GO Slim background was constructed taking into account the overall GO Slims covered by the 127,039 UniRefs id retrieved from all the collected ESTs Database. The GO Slim terms with P-value < 0.05 were considered as enriched. The same procedure was applied for the targets of the differentially expressed miRNAs. Unlike the overall analysis, the GO Slim background considered in each condition was computed from only the GO Slim of the identified target genes present in the two compared libraries from sensitive and tolerant genotypes. For functional analysis, we investigate the potential function of all identified miRNAs in the 10 investigated conditions (10 libraries) for miRNAs having at least 100 reads in at least one sequenced library. (Additional file 4: data SD3).
Availability of supporting data
All the predicted and conserved miRNAs from this study, the published miRNAs from the literature, the small RNA expression profiles are provided at the following database http://wheat.bioinfo.uqam.ca.
We would like to thank Mathieu Blanchette for his useful comments and suggestions. Thanks also to all people and institutions who gave us access to their super computers, Alix Boc for the Trex cluster, Daniel Lemire in the Licef research center for the cluster Erasme at the TeluQ, the Clumeq (Supercomputer Consortium Laval UQAM McGill and Eastern Quebec) for the clusters Colosse and Guillimin, as well as HHMMiR and Tapir authors for providing source code for this project.
This work was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grants to FS, MH and ABD; NSERC fellowship for E.L. and M.A.R., as well as the Fond de Recherche du Québec-Nature et Technologie (FRQNT) to M.L and M.A.R.
- Carrington JC, Ambros V. Role of microRNAs in plant and animal development. Science. 2003;301:336–8.View ArticlePubMedGoogle Scholar
- Jones-Rhoades MW, Bartel DP, Bartel B. MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol. 2006;57:19–53.View ArticlePubMedGoogle Scholar
- Xing S, Salinas M, Höhmann S, Berndtgen R, Huijser P. miR156-targeted and nontargeted SBP-box transcription factors act in concert to secure male fertility in Arabidopsis. Plant Cell. 2011;22:3935–50.View ArticleGoogle Scholar
- Hsieh LC, Lin SI, Shih AC, Chen JW, Lin WY, Tseng CY, et al. Uncovering small RNA-mediated responses to phosphate deficiency in Arabidopsis by deep sequencing. Plant Physiol. 2009;151:2120–32.View ArticlePubMed CentralPubMedGoogle Scholar
- Sun G, Stewart CN, Xiao P, Zhang B. MicroRNA expression analysis in the cellulosic biofuel crop Switchgrass Panicum virgatum under abiotic stress. PLoS One. 2012;7:e32017.View ArticlePubMed CentralPubMedGoogle Scholar
- Tang Z, Zhang L, Xu C, Yuan S, Zhang F, Zhen Y, et al. Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing. Plant Physiol. 2012;159:721–38.View ArticlePubMed CentralPubMedGoogle Scholar
- Zeng QY, Yang CY, Ma QB, Li XP, Dong WW, Nian H. Identification of wild soybean miRNAs and their target genes responsive to aluminum stress. BMC Plant Biol. 2012;12:182.View ArticlePubMed CentralPubMedGoogle Scholar
- Rogers K, Chen X. Biogenesis, turnover, and mode of action of plant microRNAs. Plant Cell. 2013;25:2383–99.View ArticlePubMed CentralPubMedGoogle Scholar
- Rajagopalan R, Vaucheret H, Trejo J, Bartel DP. A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006;20:3407–25.View ArticlePubMed CentralPubMedGoogle Scholar
- Sunkar R, Girke T, Jain PK, Zhu JK. Cloning and characterization of microRNAs from rice. Plant Cell. 2005;17:1397–411.View ArticlePubMed CentralPubMedGoogle Scholar
- Szarzynska B, Sobkowiak L, Pant BD, Balazadeh S, Scheible WR, Mueller-Roeber B, et al. Gene structures and processing of Arabidopsis thaliana HYL1-dependent pri-miRNAs. Nucleic Acids Res. 2009;37:3083–93.View ArticlePubMed CentralPubMedGoogle Scholar
- Piriyapongsa J, Jordan IK. Dual coding of siRNAs and miRNAs by plant transposable elements. RNA. 2008;14:814–21.View ArticlePubMed CentralPubMedGoogle Scholar
- Lucas SJ, Budak H. Sorting the wheat from the chaff: identifying miRNAs in genomic survey sequences of Triticum aestivum chromosome 1AL. PLoS One. 2012;7:e40859.View ArticlePubMed CentralPubMedGoogle Scholar
- Wei B, Cai T, Zhang R, Li A, Huo N, Li S, et al. Novel microRNAs uncovered by deep sequencing of small RNA transcriptomes in bread wheat Triticum aestivum L. and Brachypodium distachyon L Beauv. Funct Integr Genomics. 2009;9:499–511.View ArticlePubMedGoogle Scholar
- Xin M, Wang Y, Yao Y, Song N, Hu Z, Qin D, et al. Identification and characterization of wheat long non-protein coding RNAs responsive to powdery mildew infection and heat stress by using microarray analysis and SBS sequencing. BMC Plant Biol. 2011;11:61.View ArticlePubMed CentralPubMedGoogle Scholar
- Kantar M, Akpinar BA, Valárik M, Lucas SJ, Doležel J, Hernández P, et al. Subgenomic analysis of microRNAs in polyploid wheat. Funct Integr Genomics. 2012;12:465–79.View ArticlePubMedGoogle Scholar
- Ling HQ, Zhao S, Liu D, Wang J, Sun H, Zhang C, et al. Draft genome of the wheat A-genome progenitor Triticum urartu. Nature. 2013;496(7443):87–90.View ArticlePubMedGoogle Scholar
- Wang B, Sun YF, Song N, Wan XJG, Feng H, Huang LL, et al. Identification of UV-B-induced microRNAs in wheat. Genet Mol Res. 2013;12:4213–21.View ArticlePubMedGoogle Scholar
- Han J, Kong ML, Xie H, Sun QP, Nan ZJ, Zhang QZ, et al. Identification of miRNAs and their targets in wheat Triticum aestivum L. by EST analysis. Genet Mol Res. 2013;12:3793–805.View ArticlePubMedGoogle Scholar
- Kurtoglu KY, Kantar M, Lucas SJ, Budak H. Unique and conserved microRNAs in wheat chromosome 5D revealed by next-generation sequencing. PLoS One. 2013;8:1932–6203.View ArticleGoogle Scholar
- Meng F, Liu H, Wang K, Liu L, Wang S, Zha Y, et al. Development-associated microRNAs in grains of wheat Triticum aestivum L. BMC Plant Biol. 2013;13:140.View ArticlePubMed CentralPubMedGoogle Scholar
- Pandey B, Gupta OP, Pandey DM, Sharma I, Sharma P. Identification of new stress-induced microRNA and their targets in wheat using computational approach. Plant Signal Behav. 2013;8:e23932.View ArticlePubMed CentralPubMedGoogle Scholar
- Li YF, Zheng Y, Jagadeeswaran G, Sunkar R. Characterization of small RNAs and their target genes in wheat seedlings using sequencing-based approaches. Plant Sci. 2013;203–204:17–24.View ArticlePubMedGoogle Scholar
- Deng P, Nie X, Wang L, Cui L, Liu P, Tong W, et al. Computational identification and comparative analysis of miRNAs in wheat group 7 chromosomes. Plant Mol Biol Rep. 2014;32:487–500.View ArticleGoogle Scholar
- Li A, Liu D, Wu J, Zhao X, Hao M, Geng S, et al. mRNA and small RNA transcriptomes reveal Insights into dynamic homoeolog regulation of allopolyploid heterosis in nascent hexaploid Wheat. Plant Cell. 2014;26:1878–900.View ArticlePubMed CentralPubMedGoogle Scholar
- Sun F, Guo G, Du J, Guo W, Peng H, Ni Z, et al. Whole-genome discovery of miRNAs and their targets in wheat (Triticum aestivum L.). BMC Plant Bio. 2014;14:142.View ArticleGoogle Scholar
- Han R, Jian C, Lv J, Yan Y, Chi Q, Li Z, et al. Identification and characterization of microRNAs in the flag leaf and developing seed of wheat (Triticum aestivum L.). BMC Genomics. 2014;15:289.View ArticlePubMed CentralPubMedGoogle Scholar
- Pandey R, Joshi G, Bhardwaj AR, Agarwal M, Katiyar-Agarwal S. A comprehensive genome-wide study on tissue-specific and abiotic stress-specific miRNAs in Triticum aestivum. PLoS One. 2014;9(4):e95800.View ArticlePubMed CentralPubMedGoogle Scholar
- Yao Y, Guo G, Ni Z, Sunkar R, Du J, Zhu JK, et al. Cloning and characterization of microRNAs from wheat Triticum aestivum L. Genome Biol. 2007;8:R96.View ArticlePubMed CentralPubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011;39(Database issue):D152–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang Z, Yu J, Li D, Zhang Z, Liu F, Zhou X, et al. PMRD: plant microRNA database. Nucleic Acids Res. 2010;38(Database issue):D806–13.View ArticlePubMed CentralPubMedGoogle Scholar
- Brenchley R, Spannagl M, Pfeifer M, Barker GL, D’Amore R, Allen AM, et al. Analysis of the bread wheat genome using whole-genome shotgun sequencing. Nature. 2012;491:705–10.View ArticlePubMed CentralPubMedGoogle Scholar
- Adai A, Johnson C, Mlotshwa S, Archer-Evans S, Manocha V, Vance V, et al. Computational prediction of miRNAs in Arabidopsis thaliana. Genome Res. 2005;5:78–91.View ArticleGoogle Scholar
- Zhang BH, Pan XP, Wang QL, Cobb GP, Anderson TA. Identification and characterization of new plant microRNAs using EST analysis. Cell Res. 2005;15:336–60.View ArticlePubMedGoogle Scholar
- Fahlgren N, Jogdeo S, Kasschau KD, Sullivan CM, Chapman EJ, Laubinger S, et al. MicroRNA gene evolution in Arabidopsis lyrata and Arabidopsis thaliana. Plant Cell. 2010;22:1074–89.View ArticlePubMed CentralPubMedGoogle Scholar
- Breakfield NW, Corcoran DL, Petricka JJ, Shen J, Sae-Seaw J, Rubio-Somoza I, et al. High-resolution experimental and computational profiling of tissue-specific known and novel miRNAs in Arabidopsis. Genome Res. 2012;22:163–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, et al. Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol. 2008;26:407–15.View ArticlePubMedGoogle Scholar
- Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40:37–52.View ArticlePubMed CentralPubMedGoogle Scholar
- Hackenberg M, Sturm M, Langenberger D, Falcon-Perez JM, Aransay AM. miRanalyzer: a microRNA detection and analysis tool for next-generation sequencing experiments. Nucleic Acids Res. 2009;37:W68–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Hackenberg M, Rodriguez-Ezpeleta N, Aransay AM. miRanalyzer: an update on the detection and analysis of microRNAs in high-throughput sequencing experiments. Nucleic Acids Res. 2011;39:W132–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Leclercq M, Diallo AB, Blanchette M. Computational prediction of the localization of microRNAs within their pre-miRNA. Nucleic Acids Res. 2013;41:7200–11.View ArticlePubMed CentralPubMedGoogle Scholar
- Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, et al. Criteria for annotation of plant microRNAs. Plant Cell. 2008;20:3186–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Jiang P, Wu H, Wang W, Ma W, Sun X, Lu Z. MiPred: classification of real and pseudo microRNA precursors using random forest prediction model with combined features. Nucleic Acids Res. 2007;35:W339–44.View ArticlePubMed CentralPubMedGoogle Scholar
- Kadri S, Hinman V, Benos PV. HHMMiR: efficient de novo prediction of microRNAs using hierarchical hidden Markov models. BMC Bioinf. 2009;10:S35.View ArticleGoogle Scholar
- Jones-Rhoades MW, Bartel DP. Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004;14:787–99.View ArticlePubMedGoogle Scholar
- Dryanova A, Zakharov A, Gulick PJ. Data mining for miRNAs and their targets in the Triticeae. Genome. 2008;51:433–43.View ArticlePubMedGoogle Scholar
- Jin W, Li N, Zhang B, Wu F, Li W, Guo A, et al. Identification and verification of microRNA in wheat Triticum aestivum. J Plant Res. 2008;121:351–5.View ArticlePubMedGoogle Scholar
- Yin ZJ, Shen FF. Identification and characterization of conserved microRNAs and their target genes in wheat (Triticum aestivum). Genet Mol Res. 2010;9:1186–96.View ArticlePubMedGoogle Scholar
- Jeong DH, Park S, Zhai J, Gurazada SG, De Paoli E, Meyersm BC, et al. Massive analysis of rice small RNAs: Mechanistic implications of regulated miRNAs and variants for differential target RNA cleavage. Plant Cell. 2011;23:4185–207.View ArticlePubMed CentralPubMedGoogle Scholar
- Cuperus JT, Fahlgren N, Carrington JC. Evolution and functional diversification of miRNA genes. Plant Cell. 2011;23:431–42.View ArticlePubMed CentralPubMedGoogle Scholar
- Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, et al. A uniform system for microRNA annotation. RNA. 2003;9:277–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Chen HM, Chen LT, Patel K, Li YH, Baulcombe DC, Wu SH. 22-Nucleotide RNAs trigger secondary siRNA biogenesis in plants. Proc Natl Acad Sci U S A. 2010;107:15269–74.View ArticlePubMed CentralPubMedGoogle Scholar
- Ling KH, Brautigan PJ, Hahn CN, Daish T, Rayner JR, Cheah PS, et al. Deep sequencing analysis of the developing mouse brain reveals a novel microRNA. BMC Genomics. 2011;12:176.View ArticlePubMed CentralPubMedGoogle Scholar
- Dinger ME, Pang KC, Mercer TR, Mattick JS. Differentiating protein-coding and noncoding RNA: challenges and ambiguities. PLoS Comput Biol. 2008;4:e1000176.View ArticlePubMed CentralPubMedGoogle Scholar
- Wagatsuma T, Ishikawa S, Uemura M, Mitsuhashi W, Kawamura T, Khan MSH, et al. Plasma membrane lipids are the powerful components for early stage aluminum tolerance in triticale. Soil Sci Plant Nutr. 2005;51:701–4.View ArticleGoogle Scholar
- Khan MS, Tawaraya K, Sekimoto H, Koyama H, Kobayashi Y, Murayama T, et al. Relative abundance of Delta 5-sterols in plasma membrane lipids of root-tip cells correlates with aluminum tolerance of rice. Physiol Plant. 2009;135:73–83.View ArticlePubMedGoogle Scholar
- Wang JW, Wang LJ, Mao YB, Cai WJ, Xue HW, Chen XY. Control of root cap formation by microRNA-targeted auxin response factors in Arabidopsis. Plant Cell. 2005;17:2204–16.View ArticlePubMed CentralPubMedGoogle Scholar
- Gutierrez L, Mongelard G, Floková K, Pacurar DI, Novák O, Staswick P, et al. Auxin controls Arabidopsis adventitious root initiation by regulating jasmonic acid homeostasis. Plant Cell. 2012;6:2515–27.View ArticleGoogle Scholar
- Badawi M, Reddy YV, Agharbaoui Z, Tominaga Y, Danyluk J, Sarhan F, et al. Structure and functional analysis of wheat ICE Inducer of CBF expression genes. Plant Cell Physiol. 2008;49:1237–49.View ArticlePubMedGoogle Scholar
- Kim JS, Kim KA, Oh TR, Park CM, Kang H. Functional characterization of DEAD-box RNA helicases in Arabidopsis thaliana under abiotic stress conditions. Plant Cell Physiol. 2008;49:1563–71.View ArticlePubMedGoogle Scholar
- Kurepin LV, Dahal KP, Savitch LV, Singh J, Bode R, Ivanov AG, et al. Role of CBFs as integrators of chloroplast redox, phytochrome and plant hormone signaling during cold acclimation. Int J Mol Sci. 2013;14:12729–63.View ArticlePubMed CentralPubMedGoogle Scholar
- Janska A, Marsik P, Zelenkova S, Ovesna J. Cold stress and acclimation: what is important for metabolic adjustment? Plant Biol (Stuttg). 2010;12:395–405.View ArticleGoogle Scholar
- Badawi M, Danyluk J, Boucho B, Houde M, Sarhan F. The CBF gene family in hexaploid wheat and its relationship to the phylogenetic complexity of cereal CBFs. Mol Genet Genomics. 2007;277:533–54.View ArticlePubMed CentralPubMedGoogle Scholar
- Liang G, Yu D. Reciprocal regulation among miR395, APS and SULTR2,1 in Arabidopsis thaliana. Plant Signal Behav. 2010;10:1257–9.View ArticleGoogle Scholar
- Zhao X, Liu X, Guo C, Gu J, Kai X. Identification and characterization of microRNAs from wheat Triticum aestivum L. under phosphorus deprivation. J. Plant Biochem Biotechnol. 2013;22:113–23.View ArticleGoogle Scholar
- Wang L, Huang H, Fan Y, Kong B, Hu H, Hu K, et al. Effects of downregulation of microRNA-181a on H2O2-induced H9c2 cell apoptosis via the mitochondrial apoptotic pathway. Oxid Med Cell Longev. 2014;2014:960362.PubMed CentralPubMedGoogle Scholar
- Paolacci AR, Tanzarella OA, Porceddu E, Varotto S, Ciaffi M. Molecular and phylogenetic analysis of MADS-box genes of MIKC type and chromosome location of SEP-like genes in wheat (Triticum aestivum L.). Mol Genet Genomics. 2007;278:689–708.View ArticlePubMedGoogle Scholar
- Colaiacovo M, Lamontanara A, Bernardo L, Alberici R, Crosatti C, Giusti L, et al. On the complexity of miRNA-mediated regulation in plants: novel insights into the genomic organization of plant miRNAs. Biol Direct. 2012;7:15.View ArticlePubMed CentralPubMedGoogle Scholar
- Guo S, Xu Y, Liu H, Mao Z, Zhang C, Ma Y, et al. The interaction between OsMADS57 and OsTB1 modulates rice tillering via DWARF14. Nat Commun. 2013;4:1566.View ArticlePubMed CentralPubMedGoogle Scholar
- Sunkar R, Li YF, Jagadeeswaran G. Functions of microRNAs in plant stress responses. Trends Plant Sci. 2012;4:196–203.View ArticleGoogle Scholar
- Reeves PH, Murtas G, Dash S, Coupland G. Early in short days 4, a mutation in Arabidopsis that causes early flowering and reduces the mRNA abundance of the floral repressor FLC. Development. 2002;129:5349–61.View ArticlePubMedGoogle Scholar
- Murtas G, Reeves PH, Fu Y-F, Bancroft I, Dean C, Coupland G. A nuclear protease required for flowering time regulation in Arabidopsis reduces the abundance of small ubiquitin-related modifier conjugates. Plant Cell. 2003;15:2308–19.View ArticlePubMed CentralPubMedGoogle Scholar
- Li X, Hongwu B, Dafeng S, Shengyun M, Ning H, Junhui W, et al. Flowering time control in ornamental gloxinia (Sinningia speciosa). Ann Bot. 2013;111:791–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Várallyay E, Burgyán J, Havelda Z. MicroRNA detection by northern blotting using locked nucleic acid probes. Nat Protoc. 2008;3:190–6.View ArticlePubMedGoogle Scholar
- Ribeiro-dos-Santos Â, Khayat A, Silva A, Alencar D, Lobato J, Luz L, et al. Ultra-deep sequencing reveals the microRNAs expression pattern of the human stomach. PLoS One. 2010;5:e13205.View ArticlePubMed CentralPubMedGoogle Scholar
- Schulte J, Marschall T, Martin M, Rosenstie P, Mestdagh P, Schlierf S, et al. Deep sequencing reveals differential expression of microRNAs in favorable versus unfavorable neuroblastoma. Nucleic Acids Res. 2010;38:5919–28.View ArticlePubMed CentralPubMedGoogle Scholar
- Ondov BD, Varadarajan A, Passalacqua KD, Bergman NH. Efficient mapping of Applied Biosystems SOLiD sequence data to a reference genome for functional genomics applications. Bioinformatics. 2008;24:2776–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Smit AFA, Hubley R, Green P. RepeatMasker. 2010. Open-3.0.Google Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005;110:462–7.View ArticlePubMedGoogle Scholar
- Hofacker IL. Vienna RNA secondary structure server. Nucleic Acids Res. 2003;31:3429–31.View ArticlePubMed CentralPubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42:D68–73.View ArticlePubMed CentralPubMedGoogle Scholar
- Kal AJ, Van Zonneveld AJ, Benes V, Vandenberg M, Koerkamp MG, Albermann K, et al. Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources. Mol Biol Cell. 1999;10:1859–72.View ArticlePubMed CentralPubMedGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc. 1995;85:289–300.Google Scholar
- Bonnet E, He Y, Billiau K, Vandepeer Y. TAPIR, a web server for the prediction of plant microRNA targets, including target mimics. Bioinformatics. 2010;26:1566–8.View ArticlePubMedGoogle Scholar
- Binns D, Dimmer E, Huntley R, Barrell D, O’Donovan C, Apweiler R. QuickGO: a web-based tool for Gene Ontology searching. Bioinformatics. 2009;25:3045–6.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.