Developmentally regulated expression and complex processing of barley pri-microRNAs
- Katarzyna Kruszka†1,
- Andrzej Pacak†1,
- Aleksandra Swida-Barteczka†1,
- Agnieszka K Stefaniak1,
- Elzbieta Kaja1,
- Izabela Sierocka1,
- Wojciech Karlowski2,
- Artur Jarmolowski1 and
- Zofia Szweykowska-Kulinska1, 2Email author
© Kruszka et al.; licensee BioMed Central Ltd. 2013
Received: 12 October 2012
Accepted: 4 January 2013
Published: 16 January 2013
MicroRNAs (miRNAs) regulate gene expression via mRNA cleavage or translation inhibition. In spite of barley being a cereal of great economic importance, very little data is available concerning its miRNA biogenesis. There are 69 barley miRNA and 67 pre-miRNA sequences available in the miRBase (release 19). However, no barley pri-miRNA and MIR gene structures have been shown experimentally. In the present paper, we examine the biogenesis of selected barley miRNAs and the developmental regulation of their pri-miRNA processing to learn more about miRNA maturation in barely.
To investigate the organization of barley microRNA genes, nine microRNAs - 156g, 159b, 166n, 168a-5p/168a-3p, 171e, 397b-3p, 1120, and 1126 - were selected. Two of the studied miRNAs originate from one MIR168a-5p/168a-3p gene. The presence of all miRNAs was confirmed using a Northern blot approach. The miRNAs are encoded by genes with diverse organizations, representing mostly independent transcription units with or without introns. The intron-containing miRNA transcripts undergo complex splicing events to generate various spliced isoforms. We identified miRNAs that were encoded within introns of the noncoding genes MIR156g and MIR1126. Interestingly, the intron that encodes miR156g is spliced less efficiently than the intron encoding miR1126 from their specific precursors. miR397b-3p was detected in barley as a most probable functional miRNA, in contrast to rice where it has been identified as a complementary partner miRNA*. In the case of miR168a-5p/168a-3p, we found the generation of stable, mature molecules from both pre-miRNA arms, confirming evolutionary conservation of the stability of both species, as shown in rice and maize. We suggest that miR1120, located within the 3′ UTR of a protein-coding gene and described as a functional miRNA in wheat, may represent a siRNA generated from a mariner-like transposable element.
Seven of the eight barley miRNA genes characterized in this study contain introns with their respective transcripts undergoing developmentally specific processing events prior to the dicing out of pre-miRNA species from their pri-miRNA precursors. The observed tendency to maintain the intron encoding miR156g within the transcript, and preferences in splicing the miR1126-harboring intron, may suggest the existence of specific regulation of the levels of intron-derived miRNAs in barley.
KeywordsMicroRNA Pri-microRNA processing MicroRNA genes Splicing Alternative splicing Introns Barley
MicroRNAs are small, single-stranded RNAs, usually 21 nucleotides in length, for the first time identified in Ceanorhabditis elegans, and then in various other eukaryotic species which play a key regulatory role in gene expression at the posttranscriptional level . Arabidopsis thaliana was the first plant species in which the existence of miRNAs was demonstrated [2–4]. Further studies have confirmed the existence of miRNAs in all plant species studied . The regulatory roles of miRNAs have been demonstrated in plant development, signal transduction, protein degradation, response to environmental stress, and pathogen invasion . Additionally, miRNAs can regulate their own biogenesis, as shown in the case of miR838. The miR838 precursor is localized in the DICER-LIKE 1 (DCL1) intron 14. Dicing out of the pre-miR838 leads to DCL1 mRNA degradation, which decreases the level of DCL1, a key protein in miRNA biogenesis .
miRNAs, together with their almost perfect complementary partners, called miRNAs*, form a duplex located in the stem of a hairpin structure (pre-miRNA). Pre-miRNAs are embedded within primary-miRNAs (pri-miRNAs), which are long products of RNA polymerase II activity that possess their characteristic 5′ cap and 3′ polyA tail . miRNAs can be located in either arm of a pre-miRNA stem. In plants, the enzyme engaged in trimming the pri-miRNA hairpins, as well as dicing out miRNA-miRNA* duplexes, is called DCL1 [9, 10]. The DCL1 together with SERRATE (SE) and HYPONASTIC LEAVES 1 (HYL1), forms microprocessing complex [11, 12] that ensures efficiency and accuracy of pri-miRNA to miRNA processing . The efficiency of pri-miRNA recruitment to DCL1-HYL1-SE complex is stimulated by a RNA binding protein TOUGH (TGH) , whereas the accuracy of pri-miRNA processing requires HYL1 dephosphorylation triggered by a C-TERMINAL DOMAIN PHOSPHATASE-LIKE 1 (CPL1) . SE cooperates with a cap-binding complex (CBC) to ensure the proper processing of pri-miRNAs [16, 17]. Another protein involved in a proper plant miRNA biogenesis is called SICLE (SIC). However, its exact function is elusive . The 3′ termini of miRNA/miRNA* duplexes are 2′-O-methylated by HUA ENHANCER 1 (HEN1) methyltransferase to prevent 3′–5′ degradation or 3′ uridylation [19–22]. After HASTY (HST)-driven export of the duplex to the cytoplasm, the miRNA is loaded into the RNA-induced silencing complex (RISC) [23, 24], while the miRNA* is usually degraded . The miRNA-loaded RISC directs posttranscriptional silencing of the target mRNA, or triggers microRNA-directed phasing during trans-acting siRNA biogenesis . Due to the almost perfect complementarity of miRNA to its target mRNA, it is widely assumed the target mRNA is predominantly cleaved [27, 28]. However, there are examples of translation suppression without mRNA cleavage, as has been shown for the ath-miR172-triggered downregulation of APETALA2 expression level . Similar to animal models, studies on A. thaliana miRNA-action mutants have revealed the existence of translational repression by miRNAs .
As previously mentioned, miRNA biogenesis produces obligatory side products from DCL1-triggered cleavage of pre-miRNAs that results in miRNAs* originating from the strand opposite to the mature miRNA. While the expression of miRNAs* is rarely detected due to their rapid degradation , there are some examples of animal miRNA* that remain stable. It has been shown that these stable miRNA* are incorporated into RISC complexes to posttranscriptionally downregulate mRNA translation . It has been demonstrated in the case of human miRNA-155 and its partner miRNA-155*, that both of them regulate type I interferon production . An increasing amount of deep sequencing data for plant small RNAs has provided a basis for describing some examples of substantial miRNA* representation. The physiological roles of these molecules have yet to be established. Supposing miRNAs* participate in the posttranscriptional silencing of targeted mRNAs, it is still not known whether they are involved in regulation of the same biological pathway as their cognate miRNAs. Nevertheless, it has been demonstrated for the miR393/miR393* pair that they both regulate plant immune responses through different cellular pathways .
Although hundreds of plant miRNAs have been identified, their pri-miRNAs and genes are mostly unknown. As was demonstrated for Arabidopsis, plant miRNA precursors can be hundreds or thousands of nucleotides long, with many containing one or more introns that undergo constitutive or alternative splicing . While plant MIR genes are predominately located in intergenic loci, there have been an increasing number of examples of intragenic MIR loci located within introns of protein coding genes [7, 35, 36]. The mirtrons represent a class of intron-encoded miRNAs which are processed from spliced-out introns and constitute hairpin substrates for the dicing machinery [37–39]. There are 24 mirtrons identified in different plant species: five in A. thaliana, 18 in rice, and one in cassava (Manihot esculenta) [38, 40, 41]. The existence of plant polycistronic miRNA genes that carry multiple miRNAs has been shown in A. thaliana, Oryza sativa and Physcomitrella patens[42–44]. Bioinformatic analyses suggest that MIR genes may also overlap with protein coding genes. In A. thaliana, the MIR777 gene partly covers the 5′ UTR of the protein-coding gene At1g70650. Another surprising finding was the discovery of osa-miR3981, which is located in the last exon of the putative glyoxylase mRNA in rice . These reports suggest there might be additional examples of MIR genes overlapping with protein-coding genes. Regulated co-expression of protein/miRNA-coding genes remains unknown.
Barley is an economically important monocotyledonous crop plant. However, little is known about barley miRNA precursors and there is no data on MIR gene structures. High-throughput sequencing of conserved and novel small RNAs is widely used for computational scanning of genomic sequences in search of miRNA genes. A bioinformatic search of EST sequences resulted in a set of predicted potential pre-miRNAs, but until now there has been no experimental evidence validating these data [46, 47]. Currently, deep sequencing of barley small RNAs gained a set of detailed information concerning mature barley miRNAs [48, 49]. Lv et al.  has identified 259 mature miRNAs (126 conserved and 133 novel miRNAs) from which 46 miRNAs were deposited into miRBase (realease 19, http://www.mirbase.org/index.shtml) [50, 51]. Currently, there are 69 barley miRNA sequences and 67 pre-miRNA structures deposited in miRBase. However, none of the precursor structures have been validated using experimental approaches. Moreover, primary miRNA transcripts, as well as the genes for barley miRNAs, remain unknown.
In this study, we show experimental evidence for nine conserved barley miRNAs and their pre- and pri-miRNA precursors, as well as their gene structures. Our findings reveal a diverse organization of barley MIRNA genes, as well as expression levels of pre-miRNA and miRNA primary transcripts in five barley developmental stages studied. Furthermore, we show complex processing events of pri-miRNAs, with different pathways leading to various levels of mature miRNAs being observed in particular developmental stages. A regulatory role of posttranscriptional processing of pri-miRNAs is postulated.
Results and discussion
To determine the structures of miRNA genes and their transcripts, we selected eight barley cDNA nucleotide sequences for microRNAs precursors - 156g, 159b, 166n, 168a-5p/168a-3p, 171e, 397b-3p, 1120 and 1126 - deposited in GeneBank, http://www.ncbi.nlm.nih.gov/ (see Methods) for which computationally predicted hairpin structures, carrying conserved or newly estimated miRNA homologues, were described [46, 47]. Moreover, we followed their expression profiles at the level of pri-miRNAs, pre-miRNAs, and mature miRNAs in the following developmental stages: 1-, 2-, 3-, and 6-week-old, and 68-day-old plants. We designed 5′ and 3′ RACE primers to determine the full-length sequences of the pri-miRNA transcripts for all pri-miRNAs analyzed. On the basis of the nucleotide sequences and hairpin structural similarities, we have identified the barley MIRNA genes in this study as the orthologues of corresponding genes in rice or wheat. We named the barley MIRNA genes according to their matching rice or wheat orthologues. The total length of the pri-miRNA precursors was calculated on the basis of the longest 5′ and 3′ RACE products. For all barley pri-miRNAs, RT-PCR was carried out using primers designed against the 5′ and 3′ ends of the longest pri-miRNA RACE products to confirm that the longest pri-miRNA 5′ and 3′ ends belong to the same precursor molecule. In addition, qRT-PCR was performed to compare the level of all pri-miRNAs in the developmental stages studied.
The lengths and structures of eight characterized H. vulgare MIRNA genes
H. vulgare MIRNAgene
Position of mature miRNA and miRNA* sequences within gene#
Number of exons [length in bp]
Number of introns [length in bp]
Rice (osa)/wheat (tae) orthologues^
255-275 (intron 1) 349-370*
6 [126, 116, 68, 83, 145, 513]
5 [851, 1024, 4666, 76, 2091]
392-412 (exon 1) 241-261*
2 [596, 249]
359-379 (exon 1) 304-324*
2 [510, 469]
133-153 (exon 1) 180-200*
2 [261, 502]
150-170 (exon 1) 95-115*
2 [296, 486]
588-608 (exon 1) 520-540*
3288-3311 (exon 7) 3238-3260*
7 [97, 163, 304, 201, 144, 179, 697]
6 [111, 284, 275, 151, 94, 906]
1242-1264 (intron 3) 1309-1331*
7 [65, 74, 40, 190, 218, 76, 120]
6 [150, 203, 1305, 91, 605, 160]
miR397b-3p: miR* or functional miRNA?
Detailed analysis showed that the lowest level of miR397b-3p was observed in 1-week-old plants, while the highest level was seen in 3-week-old plants (Figure 1E, left middle panel). Using a hybridization probe complementary to the 5′ arm of the stem and loop of the hairpin structure, we were able to detect the precursor of miR397b-3p, which is approximately 110 nucleotides (nt) long. Amplification of the whole pri-miR397b-3p (Figure 1C) as well as real-time PCR experiments (Figure 1D) revealed the presence of a single transcript expressed almost equally in all developmental stages examined. However, there was no correlation between the level of precursor and mature miR397b-3p, most likely due to detection of putative miR397 molecules belonging to the same microRNA397 family encoded by other loci.
pre-miR159b belongs to the longest known barley precursors’ hairpins
MIR166ngenerates two transcripts with heterogeneous, developmentally specific 5′ ends
The qRT-PCR analysis confirmed the presence of the miR166n precursor in all developmental stages tested (Figure 3D, upper graph). The lowest expression level of the pri-miR166n was detected in 6-week-old plants as it was also detected in RT-PCR experiments (Figure 3C). In addition, we analyzed the existence of spliced (ΔIVS) and unspliced (+IVS) isoforms of the miR166n precursor (Figure 3D, lower graph). During barley development, both the spliced and unspliced precursor variants were present on an almost equal level in each growth stage tested.
Using Northern hybridization, we identified two pre-miR166n forms varying in length by approximately 30 nt (Figure 3E). The shorter form (pre-miR166n_S) was around 75 nt long, while the other (pre-miR166n_L) was about 100 nt long. The shorter form of the miR166n precursor may represent the hairpin structure, with the stem having miRNA and miRNA* at the base of this hairpin, while the longer form may correspond to the same hairpin form with an extended stem. The shorter form of pre-mR166n is more abundant compared to the longer one. The highest expression level of pre-miR166n detected in 1-week-old plants corresponded to the lowest level of mature miR166 observed in the same growth stage. The level of the mature miR166 reached its highest level in 2-week-old plants. In general, our observations show that the miR166 level fluctuates when various developmental stages are compared (Figure 3E).
miR168a-5p and miR168a-3p are generated from the same precursor and differ in their expression levels
Interestingly, in addition to the expression of mature miR168-5p, we also observed the expression of miR168-3p (Figure 4E). There are two possible explanations for this observation: (i) as postulated for rice and maize, both miR168-5p and miR168-3p molecules could be functional in barley [57, 63] or (ii) miR168-3p represents a relatively stable molecule of miRNA*. It should be noted that a target mRNA sequence for the miR168-3p molecule has not yet been identified. Using psRNATarget software (expectation=3.0), we identified a potential target barley mRNA [GeneBank: AK364646.1] encoding a protein similar to a ubiquitin-like protein from Triticum aestivum. The expression level of the mature miR168-5p is notably stronger than the miR168-3p. The expression profiles of miRNA168-5p/168-3p in particular developmental stages also differed, which may suggest a functional role of both miRNA species, as was reported for miR319b and miR319b.2 in Arabidopsis.
Intron-containing pri-miR171e undergoes complex splicing events during development
Using Northern hybridization, we detected the mature miR171 in all growth stages tested, with a notably elevated level in 3-week-old plants (Figure 5E). Northern analysis also revealed two pre-miR171e variants - a shorter one of approximately 75 nt (pre-miR171e_S), which might correspond to the precursor of the hairpin structure with miR171e/miR171e* at the base of the stem, and a longer variant, approximately 110 nt long (pre-miR171e_L), which might represent the same precursor with the extended stem. Interestingly, the level of shortened pre-miRNA171e was highest in the 3-week-old plants, which corresponds to the highest level of the mature miRNA observed in the same developmental stage. Taken together, the developmental analysis of the processing of the miR171e transcript shows a complex landscape of various splicing events that suggests potential regulatory role of splicing in mature microRNA biogenesis.
miR156g and miR1126 are encoded within introns of noncoding genes
The MIR156g gene consists of six exons and five introns with both miR156g and miR156g* sequences localized within the first intron (Table 1, Figure 6A). All introns carry U2-type signatures. Based on nucleotide sequence and structural similarities, we classify barley MIR156 as an orthologue of rice MIR156g (Figure 6B). Our analysis revealed that miR156g from rice is also intron-encoded, within the P0701F11.20 gene, described as encoding a hypothetical 132 aa protein. This finding suggests an evolutionary conservation of intron-encoded miR156g organization.
Using 5′ and 3′ RACE results and full transcript analyses, we detected eight splice isoforms (I–VIII, Figure 6C). Isoform I is polyadenylated within the second intron, while isoforms II to VIII are polyadenylated within the last exon (Figure 6A and C). Isoforms I–V maintain the first intron in which miR156g is embedded, but the remaining bodies of the precursors represent various alternatively spliced variants. Isoforms VI–VIII do not contain the microRNA-encoding intron, and while still representing splice variants of the same gene transcript, they cannot be named as genuine miRNA precursors. This high number of splice variants of the MIR156g gene indicates a very complex processing of its transcript.
The level of pri-miR156g is the lowest in the 6-week-old plants as revealed by qRT-PCR analysis (Figure 6G, upper graph). Using primers anchored in exons 1 and 2, we confirmed the presence of the unspliced isoform containing pre-miR156g in all developmental stages (Figure 6D). Real-time PCR analyses revealed that the expression level of the transcript with spliced intron 1 was gradually decreasing, while the level of the intron 1 containing product fluctuates during development (Figure 6G, lower graph).
Using primers anchored in exon 1 and the second intron, we detected RT-PCR products containing intron 1 in 1- and 2-week-old plants, while there was no product in the other three developmental stages (Figure 6E). These results suggest that intron 2 is spliced from the precursor more efficiently than intron 1 (see Figure 6D). Using primers specific for the first and last exons, we amplified various RT-PCR products of partially spliced MIR156g transcripts (Figure 6F). The obtained pattern of RT-PCR splice isoform products is very complex because of dynamic processing events. Therefore the observed complex pattern of multiple splice isoforms can slightly vary between biological replicates. The introns numbered 2–5 were spliced in all identified products. The only intron maintained was the first one containing pre-miR156g. These results confirm and extend our previous suggestion that intron 1 is spliced less efficiently, not only in comparison to the second intron, but also compared to the other introns of pri-miR156g.
Northern hybridization confirmed expression of the approximately 120nt long pre-miR156g in all growth stages tested. Hybridization also revealed the presence of two mature miR156, 20 and 21 nt long (Figure 6H). The 20 nt long mature miR156 was previously identified in barley using deep sequencing . A 21 nt long mature miR156 with an additional adenosine residue at the 3′ end is annotated in the databases of many eukaryotic species [50, 51]. Both the 20 and 21 nt miR156 species were equally represented in 6-week- and 68-day-old plants; however, in 1- and 2-week-old plants, primarily the 20 nt long miR156 was detectable. Both 20 and 21 nt long miR156 were expressed at the highest level in 68-day-old plants.
The MIR1126 gene consists of seven exons and six introns, with miR1126 and miR1126* localized within the third intron (Table 1, Figure 7A and B). All introns are U2 type. A search for the presence of ORF regions within the gene sequence resulted in the identification of three putative sequences: 80 aa, 87 aa, and the longest 91 aa (see Additional file 1). No conserved domains were detected in these coding sequences (CDS) and no significant similarities to known proteins were found. Thus, we concluded that the identified ORFs most probably do not represent CDS and do not encode proteins.
In the case of the intron-encoded miR1126 precursor, similar to pri-miR156g, a plethora of splicing isoforms was observed when 5′ and 3' RACE experiments, as well as RT-PCR amplifications of the full transcripts, were carried out. Among the many RT-PCR products, we identified five splice isoforms (I–V, Figure 7C and D). We found the fully spliced transcript (splice isoform III), as well as transcripts retaining the last intron (splice isoform II), and an isoform retaining both the first and the last introns (splice isoform I) (Figure 7C and D, upper panel). In contrast to MIR156g, none of the identified spliced isoforms contained miR1126 or miR1126*, thus they cannot be named pri-miRNAs. Using a 5′ primer anchored in the third intron upstream of miR1126/miR1126*, and the most peripheral 3′ primer anchored in the last exon, we were able to amplify precursors containing the intron sequence with miR1126 and miR1126* in all developmental stages studied (Figure 7D, middle panel). Real-time PCR analyses confirmed that splice isoforms I–III lacking intron 3 were present in higher amount in comparison to the precursor isoforms IV and V containing intron 3 (Figure 7E, lower graph). The observed tendency to maintain the intron-containing miR156g within the transcript, and preferences in splicing the miR1126-harboring intron, may suggest the existence of special regulation of the levels of intron-derived miRNAs in barley.
qRT-PCR of pri-miR1126 shows the highest expression level in 3-week-old plants (Figure 7E, upper graph), which is in agreement with the highest level of pri-miR1126 in 3-week-old plants detected by Northern hybridization (Figure 7F). The mature barley miR1126 molecule is 23 nt long, the same as reported for wheat miR1126 (Figure 7F). A sequence comparison between barley and wheat miR1126 revealed differences in two nucleotide positions. We confirmed the presence of the mature miR1126 and its corresponding precursor molecule in all developmental stages analyzed (Figure 7F).
Is miRNA1120, located in the 3′ UTR of a putative protein encoding gene, a functional microRNA?
The sequence and structure of the barley pre-miR1120 precursor show high similarity to its only known wheat orthologue pre-miR1120 (Figure 8B). For mRNA/pri-miR1120 transcripts, we were able to detect only the fully spliced RNA, which is probably due to rapid and efficient splicing of all introns from the primary transcript (Figure 8C). The expression level of the mRNA/pri-miR1120 was almost equal in all developmental stages tested (Figure 8D). However, the level of the mature miR1120 varies during development with the lowest amount in 2-week- and 68-day-old plants (Figure 8E). This suggests the presence of the posttranscriptional mechanisms regulating the miR1120 biogenesis. The pre-miR1120 detected by Northern blot is about 80 nt long and might correspond to the stem-loop structure predicted for pre-miR1120.
Unexpectedly, we found an 85 nt long region which included miR1120/miR1120* and displayed almost 80% sequence similarity to the short transposon element DNA/TcMar-Stowaway . Bioinformatic analysis revealed that this DNA transposable element is overrepresented in the barley genome and exists in around 400 copies. The DNA/TcMar-Stowaway transposon is widespread among both monocot and dicot plants [70, 71]. Our finding raises the question whether miR1120 is a true miRNA molecule or it represents a small noncoding RNA such as a siRNA, especially considering its size of 24 nucleotides.
In this study, we provide experimental evidence for selected mature miRNAs and their pre- and pri-miRNA structures in barley. Seven of the eight analyzed miRNA genes contain one or more introns, and their transcripts are the subject of complex processing events before hairpin pre-miRNA species are diced out from their pri-miRNA precursors. Complex alternative splicing of intron-containing transcripts generated various splice isoforms, which show variations in expression level when studied across five stages of barley growth. Two interesting examples of miRNAs encoded within introns of noncoding genes were identified. The observed tendency to maintain the intron encoding miR156g within the transcript, and preferences in splicing the miR1126-harboring intron, may suggest the existence of special regulation of the levels of intron-derived miRNAs in barley.
The discovery of developmental regulation at the level of expression of mature miRNA species as well as their precursors could help explain the regulatory role of miRNAs in economically important traits of the barley plant.
Plant material and growth conditions
Spring barley seeds, cultivar Rolap  were obtained from the Institute of Plant Genetics of the Polish Academy of Sciences (Poznan, Poland). Plants were grown in a greenhouse between August and October of 2009 with seasonal photoperiod and light conditions. Plants were grown in 5 liter pots in medium composed of 2/3 Klasmann TS1 substrate (Klasmann-Deilmann GmbH, Geeste, Germany) and 1/3 sand, and were watered to maintain optimal growth conditions. Whole plants from five growth stages and three biological replicates for each growth stadium were used in experiments. Zadoks decimal code was used to identify the developmental stages . Plants were collected when the first leaf developed, code 11 of Zadoks system (1-week-old plants, Additional file 2: Figure S1 A); after the third leaf developed, code 13 (2-week-old plants, Additional file 2: Figure S1 B); at the beginning of tillering, code 20–21 (3-week-old plants, Additional file 2: Figure S1 C); during stem elongation, code 32–36 (6-week-old plants, Additional file 2: Figure S1 D); when kernels reached milk ripeness, code 75–77 (68-day-old plants, Additional file 2: Figure S1 E). Ten plants from every growth stadium were pooled together and treated as one biological replicate.
DNA and RNA isolation
Genomic DNA (gDNA) was isolated from 1 g of 6-week-old barley plant tissue using DNeasy Plant Maxi Kit (Qiagen, Hilden, Germany); the concentration and quality of the gDNA were evaluated using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) and confirmed on a 0.6% agarose/EtBr gel. 100 mg of tissue from whole plants, collected 1, 2, 3, and 6 weeks, and 68 days after sowing, was used for total RNA isolation, using a modified method that allows for enrichment of small RNAs . The RNA for Northern blot analyses was extracted twice with 38% phenol solution saturated with 0.1 M sodium acetate (Roti Aqua Phenol, Roth, Karlsruhe, Germany), supplemented with 0.8 M guanidine thiocyanate, 0.4 M ammonium thiocyanate, 0.1 M sodium acetate, 5% glycerol, 0.5% sodium lauroylsarcosine, and 5 mM EDTA. To remove polysaccharides, the Ambion Plant RNA Isolation Aid (Life Technologies, Carlsbad, CA, USA) was used during phenol extraction. Next, three phenol/chloroform and two chloroform extractions were performed. RNA was precipitated in the presence of glycogen using 1.25 vol. of ethanol and 0.5 vol. of 0.8 M sodium citrate in a 1.2 M sodium chloride solution. The quality and quantity of RNA were measured with a NanoDrop ND-1000 spectrophotometer and an Infinite M200 PRO multimode reader (Tecan), RNA integrity was estimated on agarose gels. RNA for RT-PCRs was isolated as described above except for the additional phenol/chloroform and chloroform extractions, which were omitted, and precipitation was achieved with one vol. of isopropanol. DNA contaminants from these samples were removed with RQ1 RNase-free DNase (Promega, Madison, WI, USA). To prove the purity of RNA samples depleted of DNA traces, PCR reactions (thermal profile detailed in “Full-length cDNA of pri-miRNAs amplification”) with 1 μg of DNase-treated RNA as templates and primers amplifying the MIR171 gene fragment were performed for all biological replicates. In a positive control reaction, 1 ng of gDNA was used (Additional file 3: Figure S2 A).
Northern blot analysis of pre-miRNAs and mature miRNAs
Denaturing 8 M urea PAGE (15%) was used to separate 30 μg of RNA; the electrophoresis was run in 1x TBE buffer at a temperature of about 55°C. Both the Low Range GeneRuler DNA Ladder (Thermo Scientific, Lithuania) and 10bp DNA ladder (Invitrogen, Carlsbad, CA, USA) were loaded as length markers. RNA was transferred with the aid of a Trans-Blot Electrophoretic Transfer Cell (Bio-Rad) onto Amersham Hybond-NX nitrocellulose (GE Healthcare, Little Chalfont, Buckinghamshire, UK) and fixed using CL-1000 Ultraviolet Crosslinker (UVP). A 1-h pre-hybridization and a 16-h hybridization were performed in hybridization buffer (3.5% SDS, 0.375 M sodium phosphate dibasic, 0.125 M sodium phosphate monobasic) at 50°C for pre-miRNA analysis and at 42°C for miRNAs, with γ32P ATP-labeled (6000Ci/mmol; NEN-PerkinElmer, Boston, MA, USA) DNA oligo probes (Sigma). Pre-miRNAs and their respective mature miRNAs were detected on the same blot; a DNA probe complementary to U6 snRNA was used, and the U6 hybridization signal was taken as a loading control. Excess radioactive probe was washed out with 2x SSC, 0.1% SDS buffer, and the blots were exposed for one week to phosphorimaging screen (Fujifilm) and scanned with Fujifilm FLA5100 reader (Fujifilm Co., Ltd., Tokyo, Japan). Blots were quantified with Multi Gauge V2.2 software.
pri-miRNA 3′ RACE and 5′ RACE experiments
The 5′ and 3′ RACE cDNA template synthesis and two-step RACE experiments were conducted with the SMARTer RACE cDNA Amplification Kit (Clontech, Mountain View, CA, USA) according to the manufacturer’s protocol. PCR reactions were carried out using the Advantage 2 PCR Enzyme System (Clontech, Mountain View, CA, USA) in a Veriti thermal cycler (Applied Biosystems). Primers were designed for chosen ESTs [GenBank: BG415888.2, AJ475696.1, BQ760548.1, CA003609.1, CA009309.1, BG300360.1, BU974512.1, BJ486588.1] carrying computationally predicted hairpin structure sequences with conserved miRNAs (397b-3p, 159b, 166n, 168a-5p/168a-3p, 171e, 156g, 1126 and 1120, respectively) [46, 47]. The primer sequences are listed in Additional file 4. PCR products were cloned into the pGEM T-Easy vector (Promega, Madison, WI, USA) and sequenced (Genomed S.A., Warsaw, Poland).
5′ and 3′ genome walking
Four genome walking gDNA libraries were prepared according to the Genome Walker Universal Kit protocol (Clontech, Mountain View, CA, USA), and the PCR reactions were performed with the Advantage 2 PCR Enzyme System (Clontech, Mountain View, CA, USA). Products cloned into the pGEM T-Easy vector (Promega, Madison, WI, USA) were sequenced and compared to 5′ and 3′ RACE-obtained cDNA fragments using MAFFT software version 6 online . Intron positions were predicted by comparisons of genomic and cDNA sequences using FSPLICE software, http://linux1.softberry.com.
Full-length cDNA of pri-miRNAs amplification
cDNA templates were synthesized with oligo(dT)15 (Novazym, Poland) primer and SuperScript III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) using 1 μg of DNase-treated RNA as template. 1-, 2-, 3-, and 6-week and 68-day plant cDNAs were diluted 10–15 times, depending on the reverse transcription reaction efficiency, which was estimated by PCR amplification of a ubiquitin cDNA fragment [GenBank: X04133.1]. The purity of cDNA samples containing no gDNA was controlled by PCR amplification of a barley phosphate transporter 1 (HvPht1-1) [GenBank: AF543197.1] promoter fragment of 977 bp with primers anchored upstream of the first P1BS-like motif  (Additional file 3: Figure S2 B). Control of gDNA contamination was carried out for all biological replicates. The pri-miRNA amplifications and cDNA purity control reactions were performed with Taq DNA polymerase (Thermo Fisher Scientific, formerly Fermentas, Lithuania) or Expand High Fidelity PCR system (Roche, Mannheim, Germany) and two pri-miRNA specific primers (500 nM each) using the following thermal profile - 1 cycle: denaturation at 94°C/1 min, annealing at 65°C/30 s, elongation at 72°C/2 min; 29 cycles: denaturation at 94°C/30 s, annealing at 63°C/30 s (Δ -0.5°C/cycle), elongation at 72°C/2 min; 10 to 13 cycles, depending on the expression level of the pri-miRNA: denaturation at 94°C/30 s, annealing at 53°C/30 s, elongation 72°C/2 min. To improve amplification, Q-Solution (Qiagen, Hilden, Germany) was added to the RT-PCR mix. Genomic DNA template was used as a positive PCR control. Products of the PCR reactions were visualized with ethidium bromide on 1.2% agarose gels with GeneRuler 100 bp Plus or 1kb Plus DNA Ladders (Thermo Fisher Scientific, formerly Fermentas, Lithuania) as length markers. Primer sequences can be found in Additional file 4. Two additional biological replicates were performed for each pri-miRNA amplification presented in the Results and Discussion section (Additional file 5: Figure S3). RT-PCRs were only performed for the qualitative visualization of pri-miRNA processing products.
Quantitative real-time PCR
Three μg of DNA-free RNA was reverse-transcribed with SuperScript III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) and oligo(dT)15 (Novazym, Poland) primer. cDNA samples were diluted 2-times and 1μl was used as a template. qPCR was performed with Power SYBR® Green PCR Master MIX (Applied Biosystems, Warrington, UK) and two pri-miRNA-specific primers (final concentration 200 nM each) on 7900HT Fast Real-Time PCR System (Applied Biosystems) in 10 μl reaction volumes in 384-well plates. To amplify the pri-miR168a-5p/miR168a-3p precursor, a reverse primer complementary to exon-exon junction was used due to lack of specific product amplification using the other primers. The following thermal profile parameters were used: 10 min at 90°C, 45 cycles (or 40 cycles for pri-miRNA159b, pri-miRNA166n, pri-miRNA1126 and pri-miRNA1120) of 15 s at 95°C, and 1 min at 60°C. Each real-time PCR reaction was performed independently for three biological replicates, and for every biological replicate three (splicing isoforms analysis) or two (pri-miRNA abundance analysis) technical replicates were performed. The barley ADP-ribosylation factor 1-like [GenBank: AJ508228.2] gene fragment of 61 nt was simultaneously amplified and detected as an internal reference . Expression levels were calculated with the relative quantification method (2-ΔCt) as fold-change value and presented in a form of log102-ΔCt. The R2 values of analyzed data (≥0.997) were calculated with LinRegPCR software . Since the pri-miRNAs expression levels are lower than the reference gene, we have shown the expression profiles in a positive data range without changing the actual values by shifting the zero value of the graph’s y-axis to the basal expression level of the whole experiment . qRT-PCR primers were designed and used for pri-miRNA expression level validation. Primers designed and used for the validation of splice isoform levels were complementary to the exon-intron or exon-exon junctions. Primers are listed in Additional file 4.
The sequences of barley miRNA genes - 156g, 159b, 166n, 168a-5p/168a-3p, 171e, 397b-3p, 1120 and 1126 - were deposited in GeneBank  [GeneBank: JX121292, JX195499, JX195500, JX195501, JX195502, JX195498, JX195503 and JX195504, respectively]. Sequence analysis was performed with MAFFT version 6, http://mafft.cbrc.jp/alignment/server/index.html and NCBI Blast software, http://blast.ncbi.nlm.nih.gov/Blast.cgi. Secondary structures of pre-miRNAs were predicted using Folder Version 1.11 BETA software with RNAfold, Version 1.6.3 algorithm, http://www.tbi.univie.ac.at/RNA/[81, 82]. A pri-miRNA fragment covering at least 120 nt downstream and upstream of miRNA was used to determine miRNA/miRNA* pairing stability. Structures with the lowest minimal folding free energy (ΔG kcal/mol) were shown in this paper. Target prediction was performed with online tools available at http://plantgrn.noble.org/psRNATarget/, and Hordeum vulgare DFCI gene index HVGI release 11 database http://compbio.dfci.harvard.edu/ was used as a software input. Default parameter settings were used.
The work was supported by the European Regional Development Fund through the Innovative Economy for Poland 2007–2013, project WND-POIG.01.03.01-00-101/08 POLAPGEN-BD “Biotechnological tools for breeding cereals with increased resistance to drought” and by The National Science Centre grant No. UMO-2011/01/M/NZ2/01435.
- Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75: 843-854. 10.1016/0092-8674(93)90529-Y.View ArticlePubMed
- Llave C, Kasschau KD, Rector MA, Carringon JC: Endogenous and silencing-associated small RNAs in plants. Plant Cell. 2002, 14: 1605-1619. 10.1105/tpc.003210.PubMed CentralView ArticlePubMed
- Reinhart BJ, Weinstein EG, Rhoades MW, Bartel B, Bartel DP: MicroRNAs in plants. Genes Dev. 2002, 16: 1616-1626. 10.1101/gad.1004402.PubMed CentralView ArticlePubMed
- Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14: 787-799. 10.1016/j.molcel.2004.05.027.View ArticlePubMed
- Nozawa M, Miura S, Nei M: Origins and evolution of microRNA genes in plant species. Genome Biol Evol. 2012, 4: 230-239. 10.1093/gbe/evs002.PubMed CentralView ArticlePubMed
- Kruszka K, Pieczynski M, Windels D, Bielewicz D, Jarmolowski A, Szweykowska-Kulinska Z, Vazquez F: Role of microRNAs and other sRNAs of plants in their changing environments. J Plant Physiol. 2012, 169: 1664-1672. 10.1016/j.jplph.2012.03.009.View ArticlePubMed
- Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006, 20: 3407-3425. 10.1101/gad.1476406.PubMed CentralView ArticlePubMed
- Xie Z, Allen E, Fahlgren N, Calamar A, Givan SA, Carrington JC: Expression of Arabidopsis MIRNA genes. Plant Physiol. 2005, 138: 2145-2154. 10.1104/pp.105.062943.PubMed CentralView ArticlePubMed
- Papp J, Mette F, Aufsatz W, Daxinger L, Schauer S, Ray A, van der Winden J, Matzke M, Matzke AJM: Evidence for nuclear processing of plant micro RNA and short interfering RNA precursors. Plant Physiol. 2003, 132: 1382-1390. 10.1104/pp.103.021980.PubMed CentralView ArticlePubMed
- Kurihara Y, Watanabe Y: Arabidopsis micro-RNA biogenesis through Dicer-like 1 protein functions. Proc Natl Acad Sci USA. 2004, 101: 12753-12758. 10.1073/pnas.0403115101.PubMed CentralView ArticlePubMed
- Yang L, Liu Z, Lu F, Dong A, Huang H: SERRATE is a novel nuclear regulator in primary microRNA processing in Arabidopsis. Plant J. 2006, 47: 41-850.
- Kurihara Y, Takashi Y, Watanabe Y: The interaction between DCL1 and HYL1 is important for efficient and precise processing of pri-miRNA in plant microRNA biogenesis. RNA. 2006, 12: 206-212.PubMed CentralView ArticlePubMed
- Dong Z, Han M-H, Fedoroff N: The RNA-binding proteins HYL1 and SE promote accurate in vitro processing of pri-miRNA by DCL1. Proc Natl Acad Sci USA. 2008, 105: 9970-9975. 10.1073/pnas.0803356105.PubMed CentralView ArticlePubMed
- Ren G, Xie M, Dou Y, Zhang S, Zhang C, Yu B: Regulation of miRNA abundance by RNA binding protein TOUGH in Arabidopsis. Proc Natl Acad Sci USA. 2012, 109: 12817-12821. 10.1073/pnas.1204915109.PubMed CentralView ArticlePubMed
- Manavella PA, Hagmann J, Ott F, Laubinger S, Franz M, Macek B, Weigel D: Fast-forward genetics identifies plant CPL phosphatases as regulators of miRNA processing factor HYL1. Cell. 2012, 151: 859-870. 10.1016/j.cell.2012.09.039.View ArticlePubMed
- Laubinger S, Sachsenberg T, Zeller G, Busch W, Lohmann JU, Ratsch G, Weigel D: Dual roles of the nuclear cap-binding complex and SERRATE in pre-mRNA splicing and micro RNA processing in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2008, 105: 8795-8800. 10.1073/pnas.0802493105.PubMed CentralView ArticlePubMed
- Kim S, Yang J-Y, Xu J, Jang I-C, Prigge MJ, Chua N-H: Two Cap-Binding Proteins CBP20 and CBP80 are involved in processing primary MicroRNAs. Plant Cell Physiol. 2008, 49: 1634-1644. 10.1093/pcp/pcn146.PubMed CentralView ArticlePubMed
- Zhan X, Wang B, Li H, Lin R, Kalia RK, Zhu J-K, Chinnusamy V: Arabidopsis proline–rich protein important for development and abiotic stress tolerance is involved in micro RNA biogenesis. Proc Natl Acad Sci USA. 2012, 109: 18198-18203. 10.1073/pnas.1216199109.PubMed CentralView ArticlePubMed
- Park W, Li J, Song R, Messing J, Chen X: CARPEL FACTORY, a dicer homolog, and HEN1, a novel protein, act in microRNA metabolism in Arabidopsis thaliana. Curr Biol. 2002, 12: 1484-1495. 10.1016/S0960-9822(02)01017-5.View ArticlePubMed
- Yu B, Yang Z, Li J, Minakhina S, Yang M, Padgett RW, Steward R, Chen X: Methylation as a crucial step in plant microRNA biogenesis. Science. 2005, 307: 932-935. 10.1126/science.1107130.View ArticlePubMed
- Yang Z, Ebright YW, Yu B, Chen X: HEN1 recognizes 21–24 nt small RNA duplexes and deposits a methyl group onto the 2′ OH of the 3′ terminal nucleotide. Nucleic Acids Res. 2006, 34: 667-675. 10.1093/nar/gkj474.PubMed CentralView ArticlePubMed
- Ji L, Chen X: Regulation of small RNA stability: methylation and beyond. Cell Res. 2012, 22: 624-636. 10.1038/cr.2012.36.PubMed CentralView ArticlePubMed
- Park MY, Wu G, Gonzalez-Sulser A, Vaucheret H, Poethig RS: Nuclear processing and export of microRNAs in Arabidopsis. Proc Natl Acad Sci USA. 2005, 102: 3691-3696. 10.1073/pnas.0405570102.PubMed CentralView ArticlePubMed
- Vaucheret H, Vazquez F, Crété P, Bartel DP: The action of ARGONAUTE1 in the miRNA pathway and its regulation by the miRNA pathway are crucial for plant development. Genes Dev. 2004, 18: 1187-1197. 10.1101/gad.1201404.PubMed CentralView ArticlePubMed
- Eamens AL, Smith NA, Curtin SJ, Wang MB, Waterhouse PM: The Arabidopsis thaliana double-stranded RNA binding protein DRB1 directs guide strand selection from microRNA duplexes. RNA. 2009, 15: 2219-2235. 10.1261/rna.1646909.PubMed CentralView ArticlePubMed
- Allen E, Xie Z, Gustafson AM, Carrington JC: MicroRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005, 121: 207-221. 10.1016/j.cell.2005.04.004.View ArticlePubMed
- Llave C, Xie Z, Kasschau KD, Carrington JC: Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science. 2002, 297: 2053-2056. 10.1126/science.1076311.View ArticlePubMed
- Tang G, Reinhart BJ, Bartel DP, Zamore PD: A biochemical framework for RNA silencing in plants. Genes Dev. 2003, 17: 49-63. 10.1101/gad.1048103.PubMed CentralView ArticlePubMed
- Aukerman MJ, Sakai H: Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell. 2003, 15: 2730-2741. 10.1105/tpc.016238.PubMed CentralView ArticlePubMed
- 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: 1185-1190. 10.1126/science.1159151.View ArticlePubMed
- Okamura K, Philips MD, Duan H, Chou Y, Lai EC: The regulatory activity of microRNA* species has substantial influence on microRNA and 3′ UTR evolution. Nature Struct Mol Biol. 2008, 15: 354-363. 10.1038/nsmb.1409.View Article
- Zhou H, Huang X, Cui H, Luo X, Tang Y, Chen S: Wu L, an Shen: miR-155 and its star-form partner miR-155* cooperatively regulate type I interferon production by human plasmacytoid dendritic cells. Blood. 2010, 116: 5885-5894. 10.1182/blood-2010-04-280156.View ArticlePubMed
- Zhang X, Zhao H, Gao S, Wang WC, Katiyar-Agarwal S, Huang HD, Raikhel N, Jin H: Arabidopsis Argonaute 2 regulates innate immunity via miRNA393(*)-mediated silencing of a Golgi-localized SNARE gene, MEMB12. Mol Cell. 2011, 42: 356-366. 10.1016/j.molcel.2011.04.010.PubMed CentralView ArticlePubMed
- Szarzynska B, Sobkowiak L, Pant BD, Balazadeh S, Scheible WR, Mueller-Roeber B, Jarmolowski A, Szweykowska-Kulinska Z: Gene structures and processing of Arabidopsis thaliana HYL1-dependent pri-miRNAs. Nucleic Acids Res. 2009, 37: 3083-3093. 10.1093/nar/gkp189.PubMed CentralView ArticlePubMed
- Brown JWS, Marshall DF, Echeverria M: Intronic noncoding RNAs and splicing. Trends Plant Sci. 2008, 13: 335-342. 10.1016/j.tplants.2008.04.010.View ArticlePubMed
- Szarzynska B, Sobkowiak L, Jarmolowski A, Szweykowska-Kulinska Z: Gene structures and processing of plant pri-miRNAs. Res Adv in Nucleic Acids Res. 2011, 1: 1-12.
- Ruby JG, Jan CH, Bartel DP: Intronic microRNA precursors that bypass Drosha processing. Nature. 2007, 448: 83-86. 10.1038/nature05983.PubMed CentralView ArticlePubMed
- Zhu Q-H, Sproggs 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: 1456-1465. 10.1101/gr.075572.107.PubMed CentralView ArticlePubMed
- Jan CH, Friedman RC, Ruby JG, Bartel DP: Formation, regulation and evolution ofCaenorhabditis elegans3′ UTRs. Nature. 2011, 469: 97-101. 10.1038/nature09616.PubMed CentralView ArticlePubMed
- Meng Y, Shao C: Large-scale identification of mirtrons in Arabidopsis and rice. PLoS One. 2012, 7: e31163-10.1371/journal.pone.0031163.PubMed CentralView ArticlePubMed
- Patancun O, Lertpanyasampatha M, Sojikul P, Viboonjun U, Narangajavaha J: Computational identification of microRNAs and their targets in cassava (Manihot esculenta Crantz). Mol Biotechnol. 2012, 10.1007/s12033-012-9521-z.
- Zhang W, Gao S, Zhou X, Xia J, Chellappan P, Zhou X, Zhang X, Jin H: Multiple distinct small RNAs originate from the same microRNA precursors. Genome Biol. 2010, 11: R81-10.1186/gb-2010-11-8-r81.PubMed CentralView ArticlePubMed
- Lacombe S, Nagasaki H, Santi C, Duval D, Piegu B, Bangratz M, Breitler JC, Guiderdoni E, Brugidou C, Hirsch J, Cao X, Brice C, Panaud O, Karlowski WM, Sato Y, Echeverria M: Identification of precursor transcripts for 6 novel miRNAs expands the diversity on the genomic organisation and expression of miRNA genes in rice. BMC Plant Biol. 2008, 8: 123-10.1186/1471-2229-8-123.PubMed CentralView ArticlePubMed
- Talmor-Neiman M, Stav R, Frank W, Voss B, Arazi T: Novel micro-RNAs and intermediates of micro-RNA biogenesis from moss. Plant J. 2006, 47: 25-37. 10.1111/j.1365-313X.2006.02768.x.View ArticlePubMed
- Li T, Li H, Zhang YX, Liu JY: Identification and analysis of seven H2O2-responsive miRNAs and 32 new miRNAs in the seedlings of rice (Oryza sativaL. ssp. indica). Nucleic Acids Res. 2011, 39: 2821-2833. 10.1093/nar/gkq1047.PubMed CentralView ArticlePubMed
- Dryanova A, Zakharov A, Gulick PJ: Data mining for miRNAs and their targets in the Triticeae. Genome. 2008, 51: 433-443. 10.1139/G08-025.View ArticlePubMed
- Kantar M, Unver T, Budak H: Regulation of barley miRNAs upon dehydration stress correlated with target gene expression. Funct Integr Genomics. 2010, 10: 493-507. 10.1007/s10142-010-0181-4.View ArticlePubMed
- Schreiber AW, Shi BJ, Huang CY, Langridge P, Baumann U: Discovery of barley miRNAs through deep sequencing of short reads. BMC Genomics. 2011, 25: 129-View Article
- Lv S, Nie X, Wang L, Du X, Biradar SS, Jia X, Weining S: Identification and characterization of MicroRNAs from barley (Hordeum vulgare L.) by high-throughput sequencing. Int J Mol Sci. 2012, 13: 2973-2984. 10.3390/ijms13032973.PubMed CentralView ArticlePubMed
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34: D140-D144. 10.1093/nar/gkj112.PubMed CentralView ArticlePubMed
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (Database Issue): D152-D157.PubMed CentralView ArticlePubMed
- Benson DA, Karsch-Mizrachi I, Clark K, Lipman DJ, Ostell J, Sayers EW: GenBank. Nucleic Acids Res. 2012, 40 (Database issue): D48-D53.PubMed CentralView ArticlePubMed
- FSPLICE software: [http://linux1.softberry.com]
- Jeong DH, Park S, Zhai J, Gurazada SGR, De Paoli E, Meyers BC, Green PJ: Massive analysis of rice small RNAs: mechanistic implications of regulated microRNAs and variants for differential target RNA cleavage. Plant Cell. 2011, 23: 4185-4207. 10.1105/tpc.111.089045.PubMed CentralView ArticlePubMed
- Mette MF, van der Winden J, Matzke M, Matzke AJ: Short RNAs can identify new candidate transposable element families in Arabidopsis. Plant Physiol. 2002, 130: 6-9. 10.1104/pp.007047.PubMed CentralView ArticlePubMed
- Dezulian T, Palatnik JF, Huson DH, Weigel D: Conservation and divergence of microRNA families in plants. Genome Biol. 2005, 6: P13-10.1186/gb-2005-6-11-p13.View Article
- Zhang L, Chia JM, Kumari S, Stein JC, Liu Z, Narechania A, Maher CA, Guill K, McMullen MD, Ware D: A genome-wide characterization of microRNA genes in maize. PLoS Genet. 2009, 5: e1000716-10.1371/journal.pgen.1000716.PubMed CentralView ArticlePubMed
- Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant microRNA targets. Cell. 2002, 110: 513-520. 10.1016/S0092-8674(02)00863-2.View ArticlePubMed
- Jung JH, Park CM: MIR166/165 genes exhibit dynamic expression patterns in regulating shoot apical meristem and floral development in Arabidopsis. Planta. 2007, 225: 1327-1338. 10.1007/s00425-006-0439-1.View ArticlePubMed
- Li Y, Li W, Jin YX: Computational identification of novel family members of microRNA genes in Arabidopsis thaliana and Oryza sativa. Acta Biochim Biophys Sin. 2005, 37: 75-87. 10.1093/abbs/37.2.75.View ArticlePubMed
- Maher C, Timmermans M, Stein L, Ware D: Identifying microRNAs in plant genomes. CSB '04 Proceedings of the 2004 IEEE Computational Systems Bioinformatics Conference. 2004, DC, USA: IEEE Computer Society Washington, 718-723.
- Vaucheret H, Mallory AC, Bartel DP: AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. Mol Cell. 2006, 22: 129-136. 10.1016/j.molcel.2006.03.011.PubMed CentralView ArticlePubMed
- Wu L, Zhang Q, Zhou H, Ni FWX, Qi Y: Rice microRNA effector complexes and targets. Plant Cell. 2009, 21: 3421-3435. 10.1105/tpc.109.070938.PubMed CentralView ArticlePubMed
- Sobkowiak L, Bielewicz D, Malecka EM, Jakobsen I, Albrechtsen M, Szweykowska-Kulinska Z, Pacak A: The role of the P1BS element containing promoter-driven genes in Pi transport and homeostasis in plants. Front Plant Sci. 2012, 3: 58-PubMed CentralPubMed
- Rohrig H, Schmidt J, Miklashevichs E, Schell J, John M: Soybean ENOD40 encodes two peptides that bind to sucrose synthase. Proc Natl Acad Sci USA. 2002, 99: 1915-1920. 10.1073/pnas.022664799.PubMed CentralView ArticlePubMed
- Sharp PA, Burge CB: Classification of introns: U2-type or U12-type. Cell. 1997, 91: 875-879. 10.1016/S0092-8674(00)80479-1.View ArticlePubMed
- Lin CF, Mount SM, Jarmołowski A, Makałowski W: Evolutionary dynamics of U12-type spliceosomal introns. BMC Evol Biol. 2010, 10: 47-10.1186/1471-2148-10-47.PubMed CentralView ArticlePubMed
- Simpson CG, Brown JWS: U12-dependent intron splicing in plants. Nuclear pre-mRNA Processing in Plants: Current Topics in Microbiology and Immunology 326. Edited by: Reddy ASN, Golovkin M. 2008, Berlin Heidelberg: Springer-Verlag, 61-82.View Article
- Feschotte C, Zhang X, Wessler S: Miniature inverted-repeat transposable elements (MITEs) and their relationship with established DNA transposons. Mobile DNA II. Edited by: Craig NL, Craigie R, Gellert M, Lambowitz AM. 2002, Washington, DC: American Society for Microbiology Press, 1147-1158.View Article
- Bureau TE, Wessler SR: Stowaway: a new family of inverted-repeat elements associated with genes of both monocotyledonous and dicotyledonous plants. Plant Cell. 1994, 6: 907-916.PubMed CentralView ArticlePubMed
- Feschotte C, Jiang N, Wessler SR: Plant transposable elements: where genetics meets genomics. Nature Rev Genet. 2002, 3: 329-341.View ArticlePubMed
- Devaux P, Adamski T, Surma M: Inheritance of seed set in crosses of spring barley and Hordeum bulbosum L. Crop Sci. 1992, 32: 269-271. 10.2135/cropsci1992.0011183X003200010054x.View Article
- Zadoks JC, Chang TT, Konzak CF: A decimal code for the growth stages of cereals. Weed Res. 1974, 14: 415-421. 10.1111/j.1365-3180.1974.tb01084.x.View Article
- Pant BD, Musialak-Lange M, Nuc P, May P, Buhtz A, Kehr J, Walther D, Scheible WR: Identification of nutrient-responsive Arabidopsis and rapeseed microRNAs by comprehensive real-time polymerase chain reaction profiling and small RNA sequencing. Plant Physiol. 2009, 150: 1541-1555. 10.1104/pp.109.139139.PubMed CentralView ArticlePubMed
- Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9: 286-298. 10.1093/bib/bbn013.View ArticlePubMed
- Schumann PHD, Richardson AE, Smith FW, Delhaize E: Characterization of promoter expression patterns derived from Pht1 transporter genes of barley (Hordeum vulgare L.). J Exp Bot. 2004, 55: 855-865. 10.1093/jxb/erh103.View Article
- Rapacz M, Stepien A, Skorupa K: Internal standards for quantitative RT-PCR studies of gene expression under drought treatment in barley (Hordeum vulgare L.): the effects of developmental stage and leaf age. Acta Physiol Plant. 2012, 34: 1723-1733. 10.1007/s11738-012-0967-1.View Article
- Ramakers C, Ruijter JM, Lekanne Deprez RH, Moorman AFM: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003, 339: 62-66. 10.1016/S0304-3940(02)01423-4.View ArticlePubMed
- Bielewicz D, Dolata J, Zielezinski A, Alaba S, Szarzynska B, Szczesniak MW, Jarmolowski J, Szweykowska-Kulinska Z, Karlowski W: mirEX: a platform for comparative exploration of plant pri-miRNA expression data. Nucleic Acid Res. 2012, 40: 191-197. 10.1093/nar/gkr878.View Article
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.View ArticlePubMed
- Hansen T: Folder version 1.11 beta. [http://www.ncrnalab.dk/rnafolder/]
- Hofacker IL, Fontana W, Stadler PF, Bonhoeffer LS, Tacker M, Schuster P: Fast folding and comparison of RNA secondary structures (The Vienna RNA Package). Monatsh Chem. 1994, 125: 167-188. 10.1007/BF00818163.View Article
- Dai X, Zhao PX: psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011, 39: 155-159. 10.1093/nar/gkq766.View Article
- Quackenbush J, Cho J, Lee D, Liang F, Holt I, Karamycheva S, Parvizi B, Pertea G, Sultana R, White J: The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species. Nucleic Acids Res. 2001, 29: 159-164. 10.1093/nar/29.1.159.PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.