Skip to main content
  • Research article
  • Open access
  • Published:

Identification and characterization of microRNAs in the flag leaf and developing seed of wheat (Triticum aestivum L.)

Abstract

Background

MicroRNAs (miRNAs) regulate various biological processes in plants. Considerable data are available on miRNAs involved in the development of rice, maize and barley. In contrast, little is known about miRNAs and their functions in the development of wheat. In this study, five small RNA (sRNA) libraries from wheat seedlings, flag leaves, and developing seeds were developed and sequenced to identify miRNAs and understand their functions in wheat development.

Results

Twenty-four known miRNAs belonging to 15 miRNA families were identified from 18 MIRNA loci in wheat in the present study, including 15 miRNAs (9 MIRNA loci) first identified in wheat, 13 miRNA families (16 MIRNA loci) being highly conserved and 2 (2 MIRNA loci) moderately conserved. In addition, fifty-five novel miRNAs were also identified. The potential target genes for 15 known miRNAs and 37 novel miRNAs were predicted using strict criteria, and these target genes are involved in a wide range of biological functions. Four of the 15 known miRNA families and 22 of the 55 novel miRNAs were preferentially expressed in the developing seeds with logarithm (log2) of the fold change of 1.0 ~ 7.6, and half of them were seed-specific, suggesting that they participate in regulating wheat seed development and metabolism. From 5 days post-anthesis to 20 days post-anthesis, miR164 and miR160 increased in abundance in the developing seeds, whereas miR169 decreased, suggesting their coordinating functions in the different developmental stages of wheat seed. Moreover, 8 known miRNA families and 28 novel miRNAs exhibited tissue-biased expression in wheat flag leaves, with the logarithm of the fold changes of 0.1 ~ 5.2. The putative targets of these tissue-preferential miRNAs were involved in various metabolism and biological processes, suggesting complexity of the regulatory networks in different tissues. Our data also suggested that wheat flag leaves have more complicated regulatory networks of miRNAs than developing seeds.

Conclusions

Our work identified and characterised wheat miRNAs, their targets and expression patterns. This study is the first to elucidate the regulatory networks of miRNAs involved in wheat flag leaves and developing seeds, and provided a foundation for future studies on specific functions of these miRNAs.

Background

Wheat (Triticum aestivum L.) is a major staple crop for human diet. With the increase in global population, the shortage of food has become increasingly serious. Therefore, improvement of wheat yield is vital for relieving food shortages. Conventional breeding approaches via manipulation of genetic variation have been successful in improving the important agronomic traits of cereals; however, further wheat improvement requires deep understanding of the molecular mechanisms that control wheat growth and development [1].

Small RNAs (sRNAs) have been found to control cellular metabolism, growth and differentiation, to maintain genome integrity, and to combat viruses and mobile genetics elements in eukaryotes [2]. According to the present knowledge, these regulatory sRNAs are classified into microRNAs (miRNAs), short interfering RNAs (siRNAs), and Piwi-interacting RNAs in metazoans. miRNAs are produced from non-coding single-stranded RNA precursors (pri-miRNAs) that are transcribed by RNA polymerase II to generate an imperfect self-complementary, stem-loop secondary structure [3, 4]. In higher plants, a pri-miRNA is trimmed and spliced into an sRNA duplex consisting of an miRNA and its complementary sequence (miRNA*) by Dicer-like1 (DCL1) with the help of the dsRNA binding protein HYL1 and the dsRNA methylase HEN1 [5–7]. A mature miRNA, which is typically 20 nt to 22 nt long, is incorporated in the RNA-induced silencing complex to down-regulate the expression of its target genes in trans through a base pairing mechanism, whereas miRNA* is assumed to be degraded [5, 8]. miRNAs display near-perfect complementarity to their target mRNAs and interfere with target gene expression by mRNA cleavage, which occurs at the 10th and 11th positions of miRNAs, or by inhibition of translation in plants [9, 10]. In plants, miRNAs are the second most abundant sRNAs [11], acting as powerful endogenous regulators. For example, many distinct miRNAs target transcripts encoding for various transcription factors that control plant development and phase transition in Arabidopsis[12], rice [13] and maize [14], whereas others are involved in biotic and abiotic stress responses [15, 16]. The number of annotated miRNAs in miRBase database has exponentially increased in recent years. To date, 7,390 miRNAs have been identified in the plant kingdom, and the information has been deposited in miRBase (Release 20.0, June 2013; http://www.mirbase.org). However, most of these miRNAs have been identified in plants whose whole genome sequences are available; some of these plants include Oryza sativa (713), Populus trichocarpa (401), Arabidopsis thaliana (337), Brachypodium distachyon (464) and Zea mays (321). Similar to many gene regulatory systems, miRNAs are both conserved and diverse among plant lineages. Some miRNAs are conserved in angiosperms or even in embryophyta [16], whereas others are species specific, reflecting their fast-evolving and functionally diverging nature [16, 17]. Recent studies have investigated the functions of miRNAs in the seed development of various plants, such as Arabidopsis[18, 19], rice [13, 20, 21], maize [22], barley [23] and wheat [24]. Plants expressing the miR160-resistant auxin response factor 17 (ARF17) may cause abnormal embryo symmetry [18]. miR164 regulates NAC-domain target genes in Arabidopsis, and perturbation of miR164-directed regulation causes developmental abnormalities in embryonic, vegetative and floral organs [12]. miRNAs are involved in many regulatory pathways that control seed development in Arabidopsis; therefore, loss-of-function miRNA mutations may lead to developmental defects or death [19]. miR156 targets squamosa promoter-binding protein-like 10 (SPL10) and SPL11, and the regulation of these targets prevents premature gene expression during early embryogenesis [19]. In barley, it was found that the regulatory functions of miRNAs peak during the transition between the pre-storage and storage phases in seed development; miRNAs regulate the development of cereal grain through phytohormone response pathways [23]. A recent study has found that miR397 overexpression improves rice yield by increasing grain size and promoting panicle branching [21]. Therefore, understanding the involvement of miRNAs in plant development is necessary to elucidate their functions in seed development and to promote yield improvement not only in rice but also in other cereal crops by manipulating related miRNA genes.

To date, less than 220 miRNAs have been identified in wheat, 42 of which have been registered in the miRBase/Triticum aestivum (Release 20.0, June 2013; http://www.mirbase.org) [25–29]. Considerable data are available on the miRNAs involved in the development of rice, maize and barley. In contrast, little is known about the miRNAs and their regulatory functions in the development of wheat. Therefore, identification and functional analysis of miRNAs in wheat are urgently needed.

Wheat development is a complex event wherein the expansion and specialisation of different tissues are controlled by complicated interactions of signalling and gene regulation networks. The complex developmental events crucial to the development of a mature grain composed of the embryo, endosperm (starchy endosperm and aleurone) and outside layer (seed coat and pericarp) take approximately 40 days-post-anthesis (DPA). Grain development can be divided into three stages based on morphological changes and metabolite accumulation: pre-storage, storage (or maturation) and desiccation [30, 31]. The pre-storage phase corresponding to the first 5 DPA is featured by extensive mitotic activity in both embryo and endosperm. The transition to storage phase, characterised by dramatic transcriptional changes, occurs at 5 DPA to 10 DPA; this transition initiates the differentiation of tissues that will constitute the mature grain [30]. Throughout the storage phase, which lasts up to approximately 25 DPA, aleurone and embryonic tissues acquire desiccation tolerance, whereas endosperm cells undergo endoreduplication and accumulate storage metabolites (mainly starch and proteins) [31]. Flag leaf is essential in wheat reproduction development because it contributes approximately 45% to 58% of the total photosynthetic activity [32] and approximately 41% to 43% of the carbohydrates for grain filling [33].

Our long-term goal is to reveal the roles of miRNAs in wheat development. The objective of this study is to systematically identify miRNA species in wheat tissues during different developmental stages. To achieve this, we sequenced five sRNA libraries from wheat seedlings, flag leaves and immature seeds at 5 DPA (5-d seed), 10 DPA (10-d seed) and 20 DPA (20-d seed), respectively, to identify the miRNAs and understand their functions in wheat flag leaf and seed development. We identified known and novel miRNAs based on the presence of their precursors in wheat EST databases or wheat genome shotgun-sequence assemblies and their abundance as well as based on the detection of miRNA*. We also predicated the potential targets for the identified miRNAs and then analysed the expression of their profiles. This study is the first to provide useful information for uncovering the regulatory networks of miRNAs in wheat flag leaves and developing seeds.

Results and discussion

Diverse sRNA population in different wheat tissues or developmental stages

All sequencing data were first processed by filtering adaptor sequences and removing low-quality reads, and reads with sequences larger than 30 nt and smaller than 18 nt using the SOAPnuke software (http://soap.genomics.org.cn/) developed by BGIA, and clean reads were generated for each sRNA library. Identical reads were subsequently pooled to create a list of non-redundant sequences (unique sequences). A total of 74,590,133 clean reads representing 19,872,325 unique reads 18 nt to 30 nt long were obtained from deep sequencing of the five sRNA libraries (Table 1). The number of unique reads increased in flag leaves and developing seeds of 5-d, 10-d and 20-d. This observation is consistent with a previous study, which found the increasing unique signatures in developing barley grain samples at 1 DPA to 5 DPA, 6 DPA to 10 DPA, and 11 DPA to 15 DPA [23]. These data suggest that wheat developing seeds have a more complex sRNA population than vegetative tissues.

Table 1 Read statistics in five sRNA libraries

Only 41,352,236 of the 74,590,133 redundant reads (55.39%) and 11,028,735 of the 19,872,325 unique reads (55.45%) were matched perfectly to wheat genome shotgun-sequence assemblies (http://mips.helmholtz-muenchen.de/plant/wheat/uk454survey/index.jsp) (Table 1). The summary of sRNA sequencing in the five libraries or tissues is listed in Additional file 1.

We focused on the distribution of 18 nt to 26 nt reads in the five libraries or tissues because most of the sRNAs with known functions are 20 nt to 24 nt long. The size distribution patterns of the sRNA reads in the five wheat tissues suggest a distinct sRNA population in a particular wheat tissue or developmental stage (Figure 1). In the seedlings, the most abundant sRNA reads were 24 nt long (about 30%), followed by reads of 21 nt long (15%), the typical length of canonical miRNAs (Figure 1). These results are consistent with previous findings in Arabidopsis[34], rice [13], soybean [35] and apple [36]. However, sRNA reads from flag leaves were characterised with 24, 21 and 20 nt (27%, 30% and 22%, respectively) as the major size classes (Figures 1). The distribution of sRNA reads from flag leaves is similar to what has been observed in a previous study [37]. sRNAs with different sizes perform different functions. For example, post-transcriptional gene silencing is usually mediated by 21 nt sRNAs, whereas gene silencing mediated by RNA-dependent DNA methylation and heterochromatin maintenance are usually performed by 24 nt sRNAs [38, 39].

Figure 1
figure 1

Size distribution of redundant and unique short RNA sequences. Number of redundant and unique small RNA (sRNA) sequences from seedlings, flag leaves and immature seeds at 5 (5-d seeds), 10 (10-d seeds), and 20 (20-d seeds) days post-anthesis are shown separately. The frequency of redundant and unique 20 to 24 nt sRNAs in the seedlings, flag leaves and 5-d, 10-d and 20-d seeds are also shown, respectively. The frequency is expressed as a percentage of the total number of clean reads for each tissue.

The sRNA reads from the 5-d, 10-d and 20-d seeds shared similar size distribution patterns, with 24 nt as the dominant class of sRNAs (Figure 1). The frequencies of 24 reads decreased from the 5-d seeds to the 20-d seeds, while 21 nt increased. This result is in agreement with previous data on the early development of barley seed samples at 1 DPA to 5 DPA, 6 DPA to 10 DPA, and 11 DPA to 15 DPA [23]. However, if unique sequences are considered, 24 nt reads composed the predominant class of sRNAs (approximately 55% to 70%) in all the five libraries or tissues, as compared with the 20 to 23 nt reads (3% to 15%) (Figure 1). The five sRNA libraries shared less than 1.2% of the 21 nt sRNAs and less than 0.4% of the 24 nt sRNAs, indicating diverse sRNA populations in the tested tissues. Overall, these results indicate that distinct pools of sRNA work in different wheat tissues or developmental stages.

Characterization of highly to moderately conserved known miRNAs

In this study, sRNA reads with sequence similarity above 90% to rRNAs, tRNAs, snRNAs and snoRNAs were removed by BLASTN alignment against the data deposited in Rfam 10.0 (http://rfam.janelia.org/), and unique sRNA sequences perfectly matched to wheat genome shotgun-sequence assemblies were aligned to known plant miRNAs in miRBase (Release 20.0, http://www.mirbase.org). As a result, 24 known miRNAs representing 15 miRNA families were identified from 18 MIRNA loci in wheat (Table 2 and Additional file 2). Of these 24 known miRNAs, 9 (9 MIRNA loci) have been registered in miRBase/Triticum aestivum (Release 20.0, http://www.mirbase.org), and the remaining 15 (9 MIRNA loci) were first identified in wheat in the present study (Additional files 2, 3 and 4). Among the 15 known miRNAs (from 9 MIRNA loci in wheat), 6 pairs are miRNA/miRNA*. In the present study, these known miRNAs/families were sorted into two categories: highly conserved miRNAs and moderately conserved miRNAs. The former miRNAs are conserved in all three dicots (Arabidopsis, soybean and Populus) and three monocots (rice, maize and Brachypodium) whose genome sequences are available; the latter miRNAs are only conserved in some of these plant species, but not in the all six plant species described above (miRBase 20.0). Of the 15 known miRNA families (18 MIRNA loci), 13 (16 MIRNA loci) are highly conserved and 2 (2 MIRNA loci) are moderately conserved (Table 2 and Additional file 2).

Table 2 Known miRNAs identified in the five wheat sRNA libraries or tissues

The abundance of these known miRNAs, as reflected in normalised reads (reads per million of total miRNA reads, RPM), was further compared in this study. The highly conserved miRNAs or families showed significant variation among the families in all the five tissues (Table 2 and Additional file 2). The highest read abundance (approximately 238,000 RPM) was detected in the miR168 family and was 3.8 to 78 times more abundant than the other miRNA families, including miR156, miR166, miR167 and miR172, whose abundance ranged from about 2,900 RPM to 62,000 RPM (Table 2). Although low expression (976 RPM and 921 RPM, respectively) was observed for both miR164 and miR396 families, their expression level was still about 4 to 200 times greater than any of the 6 remaining highly conserved miRNA families (Table 2 and Additional file 2). Similarly, the two moderately conserved miRNA families exhibited great variation in abundance between the families (Table 2). Variation in expression levels of different miRNA species was also found in other plant species, such as grapevine [40] and apple [36].

Novel miRNA species in wheat

sRNA reads with no sequence similarity to known miRNAs were further analyzed for potential new miRNA species. These unique sequences with 20 to 23 nt were analyzed against wheat EST database (http://www.ncbi.nlm.nih.gov/nucest/?term=wheat) or wheat genome shotgun-sequence assemblies (http://mips.helmholtz-muenchen.de/plant/wheat/uk454survey/index.jsp). The EST sequences or the contigs from the wheat genome shotgun-sequence assemblies that matched to the sRNA reads were extracted from the databases and analyzed for the existence of pre-miRNA sequences that can form a hairpin secondary structure using miRNA prediction software MIREAP (http://sourceforge.net/projects/mireap/), the potential new miRNAs predicted meeting the common criteria (see details in Methods). Then, novel miRNAs were further identified based on both their abundance (at least 5 RPM in at least one of the five tissues examined) and the detection of miRNA*s, because accumulations of miRNA*s are strong supporting evidence for the cleavage by DCLs to produce functional mature miRNAs during miRNA biogenesis in plants [5]. A total of 55 novel miRNAs corresponding to 54 precursors which can form qualified secondary structures were identified in the present study (Table 3, Additional file 5), and of which, 6 mature miRNA sequences were previously described in wheat [24, 26]. The detailed information of all the 55 novel miRNAs, including the pre-miRNA sequences and the structures, the loci of the mature miRNA sequences and their miRNA*s, and the alignments of all novel miRNAs to their precursor sequences, were shown in Additional files 6 and 7. While, 53 miRNAs with abundance more than 5 RPM in at least one of the five tissues tested but without miRNA*s detected were identified as candidate miRNAs (Additional file 8), their loci, pre-miRNA sequences and structures, and reads in deep sequencing were also list in Additional files 6 and 7. A recent review [41] has reported that 65.5% of the known miRNAs in wheat begin with a 5′ uridine and that 56.5% are 21 nt long; these properties are typical in miRNAs from other species [42, 43]. In the present study, 56.4% (31/55) of the novel miRNAs started with a U at their 5′-end, and that 65.5% (36/55) are 21 nt long (Additional file 5). As been found in Arabidopsis and apple, new species-specific miRNAs are young miRNAs that have recently evolved and are often expressed at a lower level than conserved miRNAs [34, 36, 44]. This observation is also true for most of the 55 novel miRNAs identified in the present study (Table 3, Additional file 5).

Table 3 Summary of newly identified 55 novel miRNAs in the five wheat libraries or tissues

Target gene prediction for novel miRNAs

To understand the functions of the novel miRNAs identified, the putative targets of these miRNAs were predicted by using a web-based psRNA Target program (http://plantgrn.noble.org/psRNATarget/) with default parameters given in the Methods section. Thirty-seven of the 55 novel miRNAs (67.3%) had predicted targets that met the criteria (Additional file 9). According to the data in miRBase, 34 of the 42 (80.9%) wheat miRNAs have predicted targets, which is a little more than our result. As shown in Additional file 9, the predicted targets included transcript factors, protein kinases, enzymes, cellular components, receptors and transporters involved in multiple cellular processes, indicating the extensive functions of miRNAs in gene regulation networks. For example, seed-specific tae-miR1127b which were only present in developing seed (Additional file 5) targeted riboflavin biosynthesis protein Rib gene/amino acid permease gene, which function in seed development [45, 46]. This result is consistent with its functions reported earlier. Seed-specific tae-miR9653 and tae-miR9661 targeted Zinc finger transcription factor-like protein and F-box domain containing protein, respectively (Additional file 9). These genes have important functions throughout plant development [47, 48].

Preferential expression of different miRNA species in specific wheat tissues

Knowledge about the temporal and spatial expression of miRNAs might provide clues on where these miRNAs function. In the present study, the availability of the five sRNA datasets from the seedlings, flag leaves, 5-d seeds, 10-d seeds and 20-d seeds of wheat provided an opportunity to compare the expression profiles of these miRNAs in different tissues and developmental stages. The datasets from Solexa sequencing indicated that the majority of the miRNAs including the known and the novel miRNAs showed different levels of tissue-biased expression with logarithm (log2) of the fold changes of −4.6 ~ 5.2 for the known miRNAs and −8.6 ~ 7.6 for the novel miRNAs (Table 2 and Additional file 5), whereas the fold change was the ratio of abundance of a miRNA or family in flag leaves, 5-d seeds, 10-d seeds or 20-d seeds, to abundance of the same miRNA or family in seedlings.

The expression levels of randomly selected representative miRNA species in the five tissues were further determined through quantitative real-time RT-PCR (qPCR) using UBQ as the internal reference gene. The UBQ gene was confirmed to be relatively stable in the tested wheat tissues (data not shown). Results showed that the relative abundance of almost all the miRNAs determined by qPCR followed similar trends as the read numbers in the libraries (Figure 2). However, we also observed a discrepancy between the qPCR and the sequencing data for novel miRNA tae-miR1127b (Figure 2b), which was expressed at extremely low level in the tissues tested (Additional file 5). This difference might be ascribed to biases introduced to the very low abundant miRNA in some RNA samples when polyadenylation and reverse-transcription were performed with poly (T) adapters into cDNAs using an miRNA cDNA synthesis kit (Takara, Inc., Dalian, China). Thus, the miRNAs’ cDNAs used for qPCR contained biases. Biases could have also been introduced in the sequences of some samples during library generation or sequencing. A contradiction between in vivo RNA levels and sequencing results for miRNAs was also reported in grapevine [40] and apple [36]. Further analysis showed that a very significant correlation exists between the sequencing and the qPCR data (Pearson coefficient R2 = 0.892, P < 0.01) (Additional file 10), suggesting that the sequencing results are reliable.

Figure 2
figure 2

Comparison of the miRNA expression profiles determined by quantitative real-time RT-PCR (qPCR) and deep sequencing. a, known miRNAs; b, novel miRNAs and candidate miRNAs. In qPCR, UBQ was used as the internal reference gene, and the relative expression of each miRNA was calculated using a comparative CT (ΔΔCT) method. The miRNA sample with the lowest CT value that corresponds to the highest expression level was selected as the calibrator, in which the expression level was set as 1.0. The relative expression levels of the same miRNA in the other four samples were then normalised by comparing with the highest one in the tested tissues. Three independent biological replicates were performed in this experiment. For each sample, qPCR was performed in triplicate. Each column represents the mean of three samples, and error bars represent the standard deviation. In deep sequencing technology, read counts for each miRNA in one sample were normalised to reads per million of total miRNA reads (RPM). The relative expression of each miRNA was calculated by setting the highest RPM of each miRNA across the five samples as 1.0, and the relative expression of the same miRNA in the other four samples was its RPM divided by the highest RPM.

Therefore, the expression patterns of all the known and the novel miRNAs identified were profiled based on their sequencing results. Most of the known and the novel miRNAs showed various degrees of tissue-biased expression in the flag leaves or developing seeds, with the logarithm of the fold changes between 0.1 and 7.6 (Figure 3, Table 2 and Additional file 5). Of the 15 known miRNA families, 4 (miR169, miR166, miR164 and miR160) were preferentially expressed in the developing seeds (with the logarithm of the fold changes of 0.3 ~ 3.0 in the developing seeds, more than those in the flag leaves) (Figure 3a, Table 2). From 5 DPA to 20 DPA, miR164 and miR160 increased in abundance, whereas miR169 decreased (Figure 3a, Table 2). This result suggests that these miRNAs might have coordinating functions in different developmental stages of wheat seed. The potential targets of these known miRNAs were predicted (Additional file 11). The expression patterns of these miRNAs were closely related to their functions. miR160 and miR164 targeted Auxin response factor (ARF) and NAC transcription factor (NAC), respectively (Additional file 11), which control key steps in plant development. Auxin is crucial in seed development processes, including pattern formation, cell division and cell expansion [49]. Plants expressing the miR160-resistant ARF17 may cause abnormal embryo symmetry [18]. NAC genes have important functions in developmental processes, such as auxin signalling and defence [50, 51]. A recent research has revealed that NAC genes that regulate senescence improve protein, zinc, and iron contents in wheat grain [52]. In the present study, the expression level of miR164 increased with wheat grain development, from 135 RPM in the 5-d seeds to more than 240 RPM in the 10-d and 20-d seeds (Table 2). This result is consistent with the previously reported functions of miR164. miR169 targets a CCAAT-box transcription factor, which is involved in diverse processes, such as embryo development, flowering time control and root development [53]. The decreased abundance of miR169 from the 5-d seeds to the 20-d seeds was coordinated with its functions in seed development. Of the 55 novel miRNAs, 22 showed preferential expression in the different developmental stages of wheat seed (Figure 3b), with the logarithm of the fold change of 1.0 ~ 7.6, and half of these miRNAs (tae-miR1122b, tae-miR9653, tae-miR9654a, tae-miR9656, tae-miR9657a, tae-miR9659, tae-miR2275, tae-miR9665, tae-miR1127b, tae-miR9660, tae-miR9657b and tae-miR9667) were seed specific (Figure 3b, Additional file 5). This result suggests that these novel miRNAs might participate in regulating wheat seed development and metabolism. The predicted targets for these seed-specific miRNAs included transcript factors, enzymes, nucleosome/chromatin assembly factor, ribosome recycling factor and cellular components (Additional file 9), indicating the extensive functions of miRNAs in wheat seed development. For examples, tae-miR9657a, with an increasing abundance from the 5-d seeds to the 20-d seeds, targets a nucleosome/chromatin assembly factor which is essential for DNA replication during cell proliferation. tae-miR1127b which increased in abundance from the 5-d seeds to the 20-d seeds targets an amino acid permease gene, which is important for early seed development and plays a crucial for the uptake of amino acids into the endosperm and supplying amino acids for the developing embryo during early embryogenesis [45]. While, tae-miR9661, with a decreasing expression pattern from the 5-d seeds to the 20-d seeds, targets a F-box domain containing protein in wheat (Additional files 3 and 5). This result is in accordance with a previous report that F-box protein-encoding genes showed different transcript levels during seed development, suggesting the involvement of F-box proteins in rice seed development [54].

Figure 3
figure 3

Expression patterns of the known and the novel miRNAs based on deep-sequencing datasets. a, Known miRNAs; b, Novel miRNAs. The bars represent the scale of the relative expression levels of miRNAs (MEAN centred).

Of the 15 known miRNA families, 8 (miR396, miR168, miR156, miR172, miR159, miR398, miR1318 and miR167) showed different levels of preferential expression in wheat flag leaves, with the logarithm of the fold changes ranged from 0.5 to 5.2 as well as more than those in the developing seeds (Figure 3a, Table 2). Putative targets of these miRNAs encode a wide range of proteins, including various transcription factors, enzymes and argonaute (AGO) protein (Additional file 11), which are involved in diverse developmental and metabolic processes. Of the 55 novel miRNAs, 28 were characterised to have different degrees of flag leaf-biased expression, with the logarithm of the fold changes ranged from 0.1 to 5.2, whereas 4 (tae-miR1120c, tae-miR1130b, tae-miR5384, and tae-miR9675) were detected only in flag leaves (Figure 3b, Additional file 5), suggesting that these novel miRNAs might participate in regulating the development and metabolism in wheat flag leaves. Fourteen of the 28 novel miRNAs with flag leaf-biased expression had putative targets, including protein kinases, glycosyltransferase, functional proteins (Histone H2B.1), and enzymes (Additional file 9), all of which are important in plant signalling pathway and metabolism. The above results imply that the regulation networks of miRNAs in wheat flag leaves are very complex. Overall, all our data suggested that wheat flag leaves have more complicated regulatory networks of miRNAs than developing seeds.

To the best of our knowledge, this study is the first to identify miRNAs with tissue-biased expression in wheat flag leaves or seeds of different developmental stage on a large scale. Sufficient wheat genome sequence data are necessary to elucidate the functions of all miRNAs preferentially expressed in wheat flag leaves or developing seeds and to reveal the involved regulatory networks.

miRNA*s accumulation in specific wheat tissues

miRNA:miRNA* is excised by DCL1 to produce functional mature miRNA, whereas miRNA* is assumed to be gradually degraded during miRNA biogenesis [5]. However, previous studies have also suggested that miRNA*s can accumulate to a considerable level and down-regulate their target genes in both plants and animals [34, 55, 56]. Evidence supporting miRNA*-mediated mRNA cleavage has been recently found in Medicagotruncatula and maize [57, 58]. In the present study, 17 miRNA*s for known mature miRNAs registered in the miRbase/Tricum aestivum (Release 20.0) were detected in at least one of the five libraries or tissues tested (Additional file 12). About half of these miRNA*s had relatively high tissue-biased accumulation in the developing seeds, flag leaves, or seedling. For examples, the reads of miR171b* in the 5-d, 10-d and 20-d seeds were 88, 46 and 46, respectively, which are considerably higher than those in the seedlings (1) and flag leaves (3) (Additional file 12). The highest accumulation of miR1127* was observed in the flag leaves (106 reads), followed by the 20-d seeds (53 reads). Only 1 and 2 reads were detected in the seedlings and 5-d seeds, respectively. The highest level of miR169* was observed in the seedlings (28 reads), followed by the flag leaves (5 reads). However, mature miRNA sequences for miR171b*, miR1127* and miR169* were not found in any of the five tissues tested. These results suggested that miR171b*, miR1127* and miR169* might be de facto miRNAs with important regulatory functions in specific tissues and developmental stages. Among the 17 miRNA*s detected, 12 had no mature miRNA sequences detected in the five libraries (Additional file 12).

The differential accumulation of miRNAs and miRNA*s in specific tissues may be attributed to the different activities of proteins involved in sRNA biogenesis. Different AGO proteins in Arabidopsis harbour sRNAs with different 5′-terminal nucleotides; thus, miRNAs and miRNA*s are directed in different AGO complexes [59]. A recent study revealed that different AGO-coding genes in rice have different levels of tissue-biased expression [20]. In the present study, miR168, which targets an AGO-coding gene (Ta.71657) and is the most abundant miRNA or family in the five sRNA libraries, exhibited a tissue-preferential expression, with the highest read number in the flag leaves (Table 2 and Additional file 11).

Conclusion

Five sRNA libraries from wheat seedlings, flag leaves, 5-d, 10-d and 20-d seeds were sequenced in this study. Twenty-four known miRNAs belonging to 15 miRNA families were identified from 18 MIRNA loci in wheat, including 15 (9 MIRNA loci) first identified in wheat in the present study, 13 miRNA families (16 MIRNA loci) being highly conserved and 2 (2 MIRNA loci) moderately conserved. In addition, fifty-five novel miRNAs were also identified. The Potential target genes for 15 known miRNAs and 37 novel miRNAs were predicted using strict criteria, and these target genes are involved in a wide range of biological functions. Four of the 15 known miRNA families, including miR169, miR166, miR164 and miR160 were preferentially expressed in the developing seeds with the logarithm of the fold changes of 0.3 ~ 3.0. From 5 days post-anthesis to 20 days post-anthesis, miR164 and miR160 increased in abundance, whereas miR169 decreased, suggesting that these miRNAs have coordinating functions in the different developmental stages of wheat seed. Twenty-two of the 55 novel miRNAs also preferentially expressed in the different developmental stages of wheat seed with the logarithm of the fold change of 1.0 ~ 7.6, and half of them were seed-specific, suggesting that they participate in regulating wheat seed development and metabolism. Eight known miRNA families and 28 novel miRNAs exhibited different levels of tissue-biased expression in wheat flag leaves, with the logarithm of the fold changes ranged from 0.1 to 5.2. The potential targets of these miRNAs were involved in a wide range of biological functions, suggesting the complexity of the regulatory networks in this tissue. Our data also suggested that wheat flag leaves have more complicated regulatory networks of miRNA than developing seeds.

Sufficient wheat genome sequence annotation data are necessary to identify all miRNAs and their regulatory functions. Our dataset provides a useful source of information on miRNA regulation in wheat flag leaves and developing seeds.

Methods

Small RNA library construction and RNA sequencing

Winter wheat (T. aestivum L.) cultivar Xiaoyan 6 was used in this experiment. Wheat plants were grown under natural conditions in the experimental field of Northwest A & F University, Yangling, China (longitude 108°E, latitude 34°15′N) in 2009 to 2010. Wheat seedlings (at five-leaf stage post-vernalisation), flag leaves of heading plants and immature seeds at 5 DPA (5-d seeds), 10 DPA (10-d seeds) and 20 DPA (20-d seeds) were collected. Three independent biological replicates were included in the experiments. All tissues mentioned above were frozen in liquid nitrogen immediately after collection and then stored at −80°C until further use. Total RNA was isolated from each of these samples using TRIzol reagent (Invitrogen, Grand Island, NY, USA) according to the manufacturer’s instructions. Total RNA was extracted separately from the leaves, stems and roots of the seedlings and then mixed with equal amounts of individual parts. The integrity of the RNA samples was checked by 1% agarose gel electrophoresis. RNA samples that passed the quality check were sent to BGI (Shenzhen, China) for sRNA library construction and Solexa sequencing using standard protocols on the Illumina Hiseq 2000 platform.

Small RNA data analysis

Small RNA libraries were constructed and sequenced for the five wheat tissues. The sequencing data were deposited in NCBI Gene Expression Omnibus (GEO, http://www.ncbi.nlm.gov.geo/) under the accession number GSE50524.

All sequencing data were first processed by filtering 5' adaptor contaminants from the 50 nt tags, removing the low-quality reads and getting rid of the sequences larger than 30 nt and sequences smaller than 18 nt using the SOAPnuke software (http://soap.genomics.org.cn/) developed by BGI, sequences with over 50% homopolymer or dinucleotide repeats were filtered out, and clean reads were generated for each sRNA library. Identical reads were subsequently pooled to create a list of non-redundant sequences (unique sequences). The unique sequences were aligned against wheat genome shotgun-sequence assemblies (http://mips.helmholtz-muenchen.de/plant/wheat/uk454survey/index.jsp) [60], and those matched to the wheat genome shotgun-sequence assemblies were kept for further analysis. Unique sequences homologous (with sequence similarity above 90%) to non-coding RNAs, including rRNAs, tRNAs, siRNAs, snRNAs and snoRNAs, were removed by BLASTN alignment against the data deposited in Rfam 10.0 (http://rfam.janelia.org/). The remaining sRNA sequences were retained for miRNA identification. To identify previously known miRNAs, these sRNA sequences were checked for an exact match (in terms of sequence and length) to a known miRNA present in miRBase (Release 20.0, http://www.mirbase.org). For those known non-wheat reference miRNAs in the miRBase, their potential miRNA precursors in wheat genome and their processing patterns (miRNA/miRNA* etc.) were further detected in the same way as novel miRNAs which was described in following section. The rest sRNA sequences were retained to predicate novel miRNAs. The potential miRNA precursors were searched by extracting the wheat EST sequences (http://www.ncbi.nlm.nih.gov/nucest/?term=wheat) or the contig sequences from wheat genome shotgun-sequence assemblies surrounding the aligned sRNA sequences, and testing their potential to form a hairpin secondary structure using RNAfold (http://rna.tbi.univie.ac.at/cgi-bin/RNAfold.cgi) [61]. The characteristic hairpin structures of miRNA precursors were used to predict novel miRNA candidates using miRNA prediction software MIREAP (http://sourceforge.net/projects/mireap/) by exploring the secondary structures, the DCL1 cleavage sites and the minimum free energy of the unannotated small RNA tags. The program was run with following parameters: Minimal miRNA sequence length (20); Maximal miRNA sequence length (23); Minimal miRNA reference sequence length (20); Maximal miRNA reference sequence length (23); Maximal copy number of miRNAs on reference (20); Maximal free energy allowed for a miRNA precursor (−18 kcal/mol); Maximal space between miRNA and miRNA* (300); Minimal base pairs of miRNA and miRNA* (16); Maximal bulge of miRNA and miRNA* (2); Maximal asymmetry of miRNA/miRNA* duplex (2); Flank sequence length of miRNA precursor (20). After this analysis, novel miRNAs have observed the following rules: (1) The set of reads from the novel miRNA locus should account for more than 95% of all the precursor mapped small RNA reads, and reliable novel miRNA reads should account for more than 75% of the corresponding set of reads; (2) the miRNA* reads should have two-nucleotide 3’ overhangs; or (3) base-pairing between the miRNA and the other arm of the hairpin, which includes the miRNA*, is extensive such that there are typically four or fewer mismatched miRNA bases, five mismatched bases being allowed if the miRNA* was detected, mean-while, no asymmetric bulges larger than two nucleotides and no more than two asymmetric bulges should be present within the miRNA/miRNA* duplex [62–64]. In addition, novel miRNAs were further identified depending on both the abundance of each sequence, which was normalized as reads per million of total miRNA reads (RPM) [65], and detection of miRNA*s. The sRNA sequences with abundance at least 5 RPM in at least one of the five tissues examined and with miRNA* detected were considered as novel miRNAs, while those had abundance at least 5 RPM in at least one of the five tissues tested but without miRNA*s detected were considered as candidate miRNAs.

Quantification of miRNAs by qPCR

Total RNA from each of the five tissues was extracted as described above. miRNA abundance was detected according to previous reports [66, 67]. Briefly, total RNA (3 μg), including miRNAs, was first polyadenylated and then reverse transcribed with poly (T) adapters into cDNA using an miRNA cDNA synthesis kit (Takara, Inc., Dalian, China) according to the manufacturer’s instructions. The cDNA products of each tissue were normalised using UBQ [AF517839] as the internal reference gene [68], which was confirmed to be relatively stable in the tested wheat tissues in our laboratory. These products were used as templates for qPCR. The qPCR was performed on a CFX96 Real-time System (BIO-RAD, USA) using SYBR® Premix Ex TaqTM II (TaKaRa, Dalian, China). The miRNA-specific forward primer for each miRNA was designed based on the entire miRNA sequence, and the universal reverse primer was designed based on the adapter sequence, which was provided by miRNA cDNA synthesis kit (Takara, Inc., Dalian, China) (Additional file 13). The following qPCR program was used: denaturation at 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. Melting curve analysis with a programmed temperature ramp from 45°C to 95°C in 5 min was also performed to produce a dissociation curve for verification of amplification specificity. All reactions were run in triplicate. All primers used in this experiment are shown in Additional file 9.

The relative abundance of each miRNA was calculated by a comparative CT method (ΔΔCT) using the formula 2-ΔΔCT, where ΔΔCT= (CT miRNA − CT reference RNA) − (CT calibrator − CT reference RNA) [69]. The miRNA sample with the lowest CT value that corresponds to the highest expression level was selected as the calibrator, in which the expression level was set as 1.0. The relative expression levels of the same miRNA in the other four samples were then normalised by comparing with the highest one in the tested tissues.

miRNA target prediction

We predicted the potential targets of the newly identified miRNAs using the web-based psRNA Target program (http://plantgrn.noble.org/psRNATarget/). The following default parameters were used: Maximum expectation, 3; Length for complementarity (between the miRNA and its target) scoring, 20 [70]; Target accessibility-allowed maximum energy to unpair the target site (UPE), 25 [71]; Flanking length around target site for target accessibility analysis, 17 bp in upstream/13 bp in downstream [72]; Range of central mismatch leading to translational inhibition, 9 nt to 11 nt [73]. The custom plant transcript databases include EST database (http://www.ncbi.nlm.nih.gov/nucest/?term=wheat) and UniGene database (http://www.ncbi.nlm.nih.gov/unigene/?term=wheat) of T.aestivum in NCBI and contigs of T.aestivum [wheat gene index, release 12, http://compbio.dfci.harvard.edu/cgi-bin/tgi/gimain.pl?gudb=wheat].

References

  1. Cakir C, Gillespie ME, Scofield SR: Rapid determination of gene function by virus-induced gene silencing in wheat and barley. Crop Sci. 2010, 50: S77-S84.

    Article  CAS  Google Scholar 

  2. Jinek M, Doudna JA: A three-dimensional view of the molecular machinery of RNA interference. Nature. 2008, 457 (7228): 405-412.

    Article  Google Scholar 

  3. Xie Z, Allen E, Fahlgren N, Calamar A, Givan SA, Carrington JC: Expression of Arabidopsis MIRNA genes. Plant Physiol. 2005, 138 (4): 2145-2154. 10.1104/pp.105.062943.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  4. Lee Y, Kim M, Han J, Yeom KH, Lee S, Baek SH, Kim VN: MicroRNA genes are transcribed by RNA polymerase II. EMBO J. 2004, 23 (20): 4051-4060. 10.1038/sj.emboj.7600385.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  5. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.

    Article  CAS  PubMed  Google Scholar 

  6. Kurihara Y, Watanabe Y: Arabidopsis micro-RNA biogenesis through Dicer-like 1 protein functions. Proc Natl Acad Sci U S A. 2004, 101 (34): 12753-10.1073/pnas.0403115101.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Vazquez F, Gasciolli V, Crété P, Vaucheret H: The nuclear dsRNA binding protein HYL1 is required for microRNA accumulation and plant development, but not posttranscriptional transgene silencing. Curr Biol. 2004, 14 (4): 346-351. 10.1016/j.cub.2004.01.035.

    Article  CAS  PubMed  Google Scholar 

  8. Park MY, Wu G, Gonzalez-Sulser A, Vaucheret H, Poethig RS: Nuclear processing and export of microRNAs in Arabidopsis. Proc Natl Acad Sci U S A. 2005, 102 (10): 3691-10.1073/pnas.0405570102.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  9. Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136 (4): 669-687. 10.1016/j.cell.2009.01.046.

    Article  CAS  PubMed  Google Scholar 

  10. Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (6): 787-799. 10.1016/j.molcel.2004.05.027.

    Article  CAS  PubMed  Google Scholar 

  11. Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ: Elucidation of the small RNA component of the transcriptome. Science. 2005, 309 (5740): 1567-1569. 10.1126/science.1114112.

    Article  CAS  PubMed  Google Scholar 

  12. Mallory AC, Dugas DV, Bartel DP, Bartel B: MicroRNA regulation of NAC-domain targets is required for proper formation and separation of adjacent embryonic, vegetative, and floral organs. Current Biol: CB. 2004, 14 (12): 1035-10.1016/j.cub.2004.06.022.

    Article  CAS  Google Scholar 

  13. Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C: A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 2008, 18 (9): 1456-1465. 10.1101/gr.075572.107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  14. Lauter N, Kampani A, Carlson S, Goebel M, Moose SP: microRNA172 down-regulates glossy15 to promote vegetative phase change in maize. Proc Natl Acad Sci U S A. 2005, 102 (26): 9412-9417. 10.1073/pnas.0503927102.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. Sunkar R, Chinnusamy V, Zhu J, Zhu J-K: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci. 2007, 12 (7): 301-309. 10.1016/j.tplants.2007.05.001.

    Article  CAS  PubMed  Google Scholar 

  16. Cuperus JT, Fahlgren N, Carrington JC: Evolution and functional diversification of MIRNA genes. Plant Cell Online. 2011, 23 (2): 431-442. 10.1105/tpc.110.082784.

    Article  CAS  Google Scholar 

  17. Axtell MJ, Snyder JA, Bartel DP: Common functions for diverse small RNAs of land plants. Plant Cell. 2007, 19 (6): 1750-1769. 10.1105/tpc.107.051706.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  18. Mallory AC, Bartel DP, Bartel B: MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. Plant Cell Online. 2005, 17 (5): 1360-1375. 10.1105/tpc.105.031716.

    Article  CAS  Google Scholar 

  19. Nodine MD, Bartel DP: MicroRNAs prevent precocious gene expression and enable pattern formation during plant embryogenesis. Genes Dev. 2010, 24 (23): 2678-2692. 10.1101/gad.1986710.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. Xue LJ, Zhang JJ, Xue HW: Characterization and expression profiles of miRNAs in rice seeds. Nucleic Acids Res. 2009, 37 (3): 916-930. 10.1093/nar/gkn998.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Zhang YC, Yu Y, Wang CY, Li ZY, Liu Q, Xu J, Liao JY, Wang XJ, Qu LH, Chen F, Xin P, Yan C, Chu J, Li HQ, Chen YQ: Overexpression of microRNA OsmiR397 improves rice yield by increasing grain size and promoting panicle branching. Nat Biotechnol. 2013, 31 (9): 848-852. 10.1038/nbt.2646.

    Article  CAS  PubMed  Google Scholar 

  22. Kang M, Zhao Q, Zhu D, Yu J: Characterization of microRNAs expression during maize seed development. BMC Genomics. 2012, 13 (1): 360-10.1186/1471-2164-13-360.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  23. Curaba J, Spriggs A, Taylor J, Li Z, Helliwell C: miRNA regulation in the early development of barley seed. BMC Plant Biol. 2012, 12 (1): 120-10.1186/1471-2229-12-120.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Meng F, Liu H, Wang K, Liu L, Wang S, Zhao Y, Yin J, Li Y: Development-associated microRNAs in grains of wheat (Triticum aestivum L.). BMC Plant Biol. 2013, 13 (1): 140-10.1186/1471-2229-13-140.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Yao Y, Guo G, Ni Z, Sunkar R, Du J, Zhu JK, Sun Q: Cloning and characterization of microRNAs from wheat (Triticum aestivum L.). Genome Biol. 2007, 8 (6): R96-10.1186/gb-2007-8-6-r96.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Wei B, Cai T, Zhang R, Li A, Huo N, Li S, Gu YQ, Vogel J, Jia J, Qi Y: Novel microRNAs uncovered by deep sequencing of small RNA transcriptomes in bread wheat (Triticum aestivum L.) and Brachypodium distachyon (L.) Beauv. Funct Integr Genom. 2009, 9 (4): 499-511. 10.1007/s10142-009-0128-9.

    Article  CAS  Google Scholar 

  27. Xin M, Wang Y, Yao Y, Xie C, Peng H, Ni Z, Sun Q: Diverse set of microRNAs are responsive to powdery mildew infection and heat stress in wheat (Triticum aestivum L.). BMC Plant Biol. 2010, 10 (1): 123-10.1186/1471-2229-10-123.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Tang Z, Zhang L, Xu C, Yuan S, Zhang F, Zheng Y, Zhao C: Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing. Plant Physiol. 2012, 159 (2): 721-738. 10.1104/pp.112.196048.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. 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: 17-24.

    Article  PubMed  Google Scholar 

  30. Bennett M, Smith J, Barclay I: Early seed development in the Triticeae. Phil Trans Roy Soc Lond B Biol Sci. 1975, 272 (916): 199-227. 10.1098/rstb.1975.0083.

    Article  Google Scholar 

  31. Sabelli PA, Larkins BA: The development of endosperm in grasses. Plant Physiol. 2009, 149 (1): 14-26. 10.1104/pp.108.129437.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Xu JS, Zhao HY: Canopy photosynthesis capacity and the contribution from different organs in high-yielding winter wheat. Acta Agronomica Sinica. 1995, 2: 012-

    Google Scholar 

  33. Sharma S, Sain R, Sharma R: The genetic control of flag leaf length in normal and late sown durum wheat. J Agric Sci. 2003, 141: 323-331. 10.1017/S0021859603003642.

    Article  CAS  Google Scholar 

  34. Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006, 20 (24): 3407-3425. 10.1101/gad.1476406.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Song QX, Liu YF, Hu XY, Zhang WK, Ma B, Chen SY, Zhang JS: Identification of miRNAs and their target genes in developing soybean seeds by deep sequencing. BMC Plant Biol. 2011, 11 (1): 5-10.1186/1471-2229-11-5.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  36. Xia R, Zhu H, An Y, Beers EP, Liu Z: Apple miRNAs and tasiRNAs with novel regulatory networks. Genome Biol. 2012, 13 (6): R47-10.1186/gb-2012-13-6-r47.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Cantu D, Vanzetti LS, Sumner A, Dubcovsky M, Matvienko M, Distelfeld A, Michelmore RW, Dubcovsky J: Small RNAs, DNA methylation and transposable elements in wheat. BMC Genomics. 2010, 11 (1): 408-10.1186/1471-2164-11-408.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Slotkin RK, Martienssen R: Transposable elements and the epigenetic regulation of the genome. Nat Rev Genet. 2007, 8 (4): 272-285.

    Article  CAS  PubMed  Google Scholar 

  39. Schreiber AW, Shi BJ, Huang CY, Langridge P, Baumann U: Discovery of barley miRNAs through deep sequencing of short reads. BMC Genomics. 2011, 12 (1): 129-10.1186/1471-2164-12-129.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Pantaleo V, Szittya G, Moxon S, Miozzi L, Moulton V, Dalmay T, Burgyan J: Identification of grapevine microRNAs and their targets using high‒throughput sequencing and degradome analysis. Plant J. 2010, 62 (6): 960-976.

    CAS  PubMed  Google Scholar 

  41. Yao YY, Sun QX: Exploration of small non coding RNAs in wheat (Triticum aestivum L.). Plant Mol Biol. 2012, 80 (1): 67-73. 10.1007/s11103-011-9835-4.

    Article  CAS  PubMed  Google Scholar 

  42. Krol J, Sobczak K, Wilczynska U, Drath M, Jasinska A, Kaczynska D, Krzyzosiak WJ: Structural features of microRNA (miRNA) precursors and their relevance to miRNA biogenesis and small interfering RNA/short hairpin RNA design. J Biol Chem. 2004, 279 (40): 42230-42239. 10.1074/jbc.M404931200.

    Article  CAS  PubMed  Google Scholar 

  43. Adai A, Johnson C, Mlotshwa S, Archer-Evans S, Manocha V, Vance V, Sundaresan V: Computational prediction of miRNAs in Arabidopsis thaliana. Genome Res. 2005, 15 (1): 78-91. 10.1101/gr.2908205.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  44. Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Law TF, Grant SR, Dangl JL: High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes. PLoS One. 2007, 2 (2): e219-10.1371/journal.pone.0000219.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Schmidt R, Stransky H, Koch W: The amino acid permease AAP8 is important for early seed development in Arabidopsis thaliana. Planta. 2007, 226 (4): 805-813. 10.1007/s00425-007-0527-x.

    Article  CAS  PubMed  Google Scholar 

  46. D'Ovidio R, Masci S: The low-molecular-weight glutenin subunits of wheat gluten. J Cereal Sci. 2004, 39 (3): 321-339. 10.1016/j.jcs.2003.12.002.

    Article  Google Scholar 

  47. Moura DS, Bergey DR, Ryan CA: Characterization and localization of a wound-inducible type I serine-carboxypeptidase from leaves of tomato plants (Lycopersicon esculentum Mill.). Planta. 2001, 212 (2): 222-230. 10.1007/s004250000380.

    Article  CAS  PubMed  Google Scholar 

  48. Imaizumi T, Schultz TF, Harmon FG, Ho LA, Kay SA: FKF1 F-box protein mediates cyclic degradation of a repressor of CONSTANS in Arabidopsis. Science. 2005, 309 (5732): 293-297. 10.1126/science.1110586.

    Article  CAS  PubMed  Google Scholar 

  49. Schruff MC, Spielman M, Tiwari S, Adams S, Fenby N, Scott RJ: The AUXIN RESPONSE FACTOR 2 gene of Arabidopsis links auxin signalling, cell division, and the size of seeds and other organs. Development. 2006, 133 (2): 251-261.

    Article  CAS  PubMed  Google Scholar 

  50. Olsen AN, Ernst HA, Leggio LL, Skriver K: NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci. 2005, 10 (2): 79-87. 10.1016/j.tplants.2004.12.010.

    Article  CAS  PubMed  Google Scholar 

  51. Guo Y, Gan S: AtNAP, a NAC family transcription factor, has an important role in leaf senescence. Plant J. 2006, 46 (4): 601-612. 10.1111/j.1365-313X.2006.02723.x.

    Article  CAS  PubMed  Google Scholar 

  52. Uauy C, Distelfeld A, Fahima T, Blechl A, Dubcovsky J: A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science. 2006, 314 (5803): 1298-1301. 10.1126/science.1133649.

    Article  CAS  PubMed  Google Scholar 

  53. Siefers N, Dang KK, Kumimoto RW, Bynum WE, Tayrose G, Holt BF: Tissue-specific expression patterns of Arabidopsis NF-Y transcription factors suggest potential for extensive combinatorial complexity. Plant Physiol. 2009, 149 (2): 625-641.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  54. Jain M, Nijhawan A, Arora R, Agarwal P, Ray S, Sharma P, Kapoor S, Tyagi AK, Khurana JP: F-box proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress. Plant Physiol. 2007, 143 (4): 1467-1483. 10.1104/pp.106.091900.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Okamura K, Phillips MD, Tyler DM, Duan H, Chou YT, Lai EC: The regulatory activity of microRNA* species has substantial influence on microRNA and 3′ UTR evolution. Nat Struct Mol Biol. 2008, 15 (4): 354-363. 10.1038/nsmb.1409.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  56. Yang JS, Phillips MD, Betel D, Mu P, Ventura A, Siepel AC, Chen KC, Lai EC: Widespread regulatory activity of vertebrate microRNA* species. RNA. 2011, 17 (2): 312-326. 10.1261/rna.2537911.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Zhao M, Tai H, Sun S, Zhang F, Xu Y, Li W: Cloning and characterization of maize miRNAs involved in responses to nitrogen deficiency. PLoS One. 2012, 7 (1): e29669-10.1371/journal.pone.0029669.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  58. Devers EA, Branscheid A, May P, Krajinski F: Stars and symbiosis: microRNA-and microRNA*-mediated transcript cleavage involved in arbuscular mycorrhizal symbiosis. Plant Physiol. 2011, 156 (4): 1990-2010. 10.1104/pp.111.172627.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, Wu L, Li S, Zhou H, Long C: Sorting of Small RNAs into Arabidopsis Argonaute Complexes Is Directed by the 5'Terminal Nucleotide. Cell. 2008, 133 (1): 116-127. 10.1016/j.cell.2008.02.034.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  60. Brenchley R, Spannagl M, Pfeifer M, Barker GL, D’Amore R, Allen AM, McKenzie N, Kramer M, Kerhornou A, Bolser D, Kay S, Waite D, Trick M, Bancroft I, Gu Y, Huo N, Luo MC, Sehgal S, Gill B, Kianian S, Anderson O, Kersey P, Dvorak J, McCombie WR, Hall A, Mayer KF, Edwards KJ, Bevan MW, Hall N: Analysis of the bread wheat genome using whole-genome shotgun sequencing. Nature. 2012, 491 (7426): 705-710. 10.1038/nature11650.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  61. Hofacker IL, Fontana W, Stadler PF, Bonhoeffer LS, Tacker M, Schuster P: Fast folding and comparison of RNA secondary structures. Monatshefte für Chemie/Chemical Monthly. 1994, 125 (2): 167-188. 10.1007/BF00818163.

    Article  CAS  Google Scholar 

  62. Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths JS, Jacobsen SE, Mallory AC, Martienssen RA, Poething RS, Qi YJ, Vaucheret H, Voinnet O, Watanabe Y, Weigel D: Criteria for annotation of plant MicroRNAs. Plant Cell. 2008, 20 (12): 3186-3190. 10.1105/tpc.108.064311.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  63. Kozomara A, Griffiths-Jones S: miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42 (D1): D68-D73. 10.1093/nar/gkt1181.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  64. Taylor RS, Tarver JE, Hiscock SJ, Donoghue PC: Evolutionary history of plant microRNAs. Trends Plant Sci. 2014, 19 (3): 175-182. 10.1016/j.tplants.2013.11.008.

    Article  CAS  PubMed  Google Scholar 

  65. AC’t Hoen P, Friedländer MR, Almlöf J, Sammeth M, Pulyakhina I, Anvar SY, Laros JF, Buermans HP, Karlberg O, Brännvall M: Reproducibility of high-throughput mRNA and small RNA sequencing across laboratories. Nat Biotechnol. 2013, 31 (11): 1015-1022. 10.1038/nbt.2702.

    Article  Google Scholar 

  66. Shi R, Chiang VL: Facile means for quantifying microRNA expression by real-time PCR. Biotechniques. 2005, 39 (4): 519-525. 10.2144/000112010.

    Article  CAS  PubMed  Google Scholar 

  67. Han R, Yan Y, Jian C, Zhao HX: Comparison of two miRNA quantification methods for assaying expression profiles of miRNAs in wheat (Triticum aestivum L.). J Integrative Agr. 2014, 13 (4): 60345-60354.

    Article  Google Scholar 

  68. Beccari G, Covarelli L, Nicholson P: Infection processes and soft wheat response to root rot and crown rot caused by Fusarium culmorum. Plant Pathol. 2011, 60 (4): 671-684. 10.1111/j.1365-3059.2011.02425.x.

    Article  CAS  Google Scholar 

  69. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

  70. Zhang Y: miRU: an automated plant miRNA target prediction server. Nucleic Acids Res. 2005, 33: W701-W704. 10.1093/nar/gki383.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  71. Muckstein U, Tafer H, Hackermuller J, Bernhart SH, Stadler PF, Hofacker IL: Thermodynamics of RNA-RNA binding. Bioinformatics. 2006, 22 (10): 1177-1182. 10.1093/bioinformatics/btl024.

    Article  PubMed  Google Scholar 

  72. Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nat Genet. 2007, 39 (10): 1278-1284. 10.1038/ng2135.

    Article  CAS  PubMed  Google Scholar 

  73. Brodersen P, Sakvarelidze-Achard L, Bruun-Rasmussen M, Dunoyer P, Yamamoto YY, Sieburth L, Voinnet O: Widespread translational inhibition by plant miRNAs and siRNAs. Science. 2008, 320 (5880): 1185-1190. 10.1126/science.1159151.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was financially supported by grants from the National Nature Science Foundation of China (30871578) and Key Projects of National Plant Transgenic Genes of China (2011ZX08002-004).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Huixian Zhao.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HZ conceived and designed the experiments. RH, YY and CJ performed the experiments. RH, JL, YY, QC, ZL, JZ and XL analysed the data. RH, QW and HZ wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Summary of small RNA sequencing in the five wheat sRNA libraries or tissues.(XLS 26 KB)

Additional file 2: Known miRNAs identified in five wheat sRNA libraries or tissues.(XLSX 14 KB)

12864_2014_5982_MOESM3_ESM.txt

Additional file 3: The detail information of the loci, the precursor sequences and the structures for 15 known miRNAs first identified in wheat.(TXT 14 KB)

Additional file 4: The predicted pre-miRNA structures for 15 known miRNAs first identified in wheat.(DOCX 22 KB)

Additional file 5: The details of the 55 novel miRNAs identified in five wheat sRNA libraries or tissues.(XLS 48 KB)

12864_2014_5982_MOESM6_ESM.txt

Additional file 6: The detail information of the loci, the corresponding precursor sequences and the structures for 55 novel miRNAs and 53 candidate miRNAs in wheat.(TXT 273 KB)

12864_2014_5982_MOESM7_ESM.doc

Additional file 7: The predicted pre-miRNA structures for 55 novel miRNAs (tae-miR1120b~tae-miR9679) and 53 candidate miRNAs (Tae-1~Tae-53).(DOC 218 KB)

Additional file 8: The details of 53 candidate miRNAs identified in five wheat sRNA libraries or tissues.(XLSX 23 KB)

Additional file 9: The Potential wheat targets predicted for novel miRNAs identified in five wheat tissues.(XLSX 18 KB)

12864_2014_5982_MOESM10_ESM.docx

Additional file 10: Correlation between the deep sequencing data and the quantitative real time RT-PCR (qPCR) data.(DOCX 151 KB)

Additional file 11: The potential targets predicted for the known miRNAs in wheat.(XLSX 13 KB)

12864_2014_5982_MOESM12_ESM.xlsx

Additional file 12: The abundance of the miRNA star sequences for the known wheat miRNAs detected in five wheat small RNA libraries or tissues.(XLSX 13 KB)

Additional file 13: The primers used in this study.(DOCX 16 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.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 ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Han, R., Jian, C., Lv, J. et al. Identification and characterization of microRNAs in the flag leaf and developing seed of wheat (Triticum aestivum L.). BMC Genomics 15, 289 (2014). https://doi.org/10.1186/1471-2164-15-289

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-289

Keywords