Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

A computational-based update on microRNAs and their targets in barley (Hordeum vulgare L.)

Moreno Colaiacovo*, Annalisa Subacchi, Paolo Bagnaresi, Antonella Lamontanara, Luigi Cattivelli and Primetta Faccioli

Author Affiliations

CRA-Genomics Research Centre, via S.Protaso 302, I-29107 Fiorenzuola d'Arda (Pc), Italy

For all author emails, please log on.

BMC Genomics 2010, 11:595  doi:10.1186/1471-2164-11-595


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/11/595


Received:30 April 2010
Accepted:22 October 2010
Published:22 October 2010

© 2010 Colaiacovo et al; licensee 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.

Abstract

Background

Many plant species have been investigated in the last years for the identification and characterization of the corresponding miRNAs, nevertheless extensive studies are not yet available on barley (at the time of this writing). To extend and to update information on miRNAs and their targets in barley and to identify candidate polymorphisms at miRNA target sites, the features of previously known plant miRNAs have been used to systematically search for barley miRNA homologues and targets in the publicly available ESTs database. Matching sequences have then been related to Unigene clusters on which most of this study was based.

Results

One hundred-fifty-six microRNA mature sequences belonging to 50 miRNA families have been found to significantly match at least one EST sequence in barley. As expected on the basis of phylogenetic relations, miRNAs putatively orthologous to those of Triticum are significantly over-represented inside the set of identified barley microRNA mature sequences. Many previously known and several putatively new miRNA/target pairs have been identified. When the predicted microRNA targets were grouped into functional categories, biological processes previously known to be regulated by miRNAs, such as development and response to biotic and abiotic stress, have been highlighted and most of the target molecular functions were related to transcription regulation. Candidate microRNA coding genes have been reported and genetic variation (SNPs/indels) both in functional regions of putative miRNAs (mature sequence) and at miRNA target sites has been found.

Conclusions

This study has provided an update of the information on barley miRNAs and their targets representing a foundation for future studies. Many of previously known plant microRNAs have homologues in barley with expected important roles during development, nutrient deprivation, biotic and abiotic stress response and other important physiological processes. Putative polymorphisms at miRNA target sites have been identified and they can represent an interesting source for the identification of functional genetic variability.

Background

MicroRNAs (miRNAs) are a class of non-coding small RNAs with fundamental roles in key plant biological processes such as development, signal transduction and environmental stress response [1]. miRNAs act on gene regulation at post-transcriptional level, a phenomenon known in plants as PTGS (Post Transcriptional Gene Silencing), through sequence-based interaction with target mRNAs.

Many plant species have been investigated during recent years for miRNAs identification and characterization. The current information available on barley refers to two papers [2,3]. In particular, the paper of Dryanova et al. reports detailed information on both targets and miRNA coding sequences from Hordeum vulgare and for other members of Triticeae tribe, to which barley belongs [2]. However, extensive studies describing the organization of miRNA families, specifically in barley, are not yet available (at the time of this writing) and no miRNAs have been deposited in the publicly available miRNA database (miRBase, http://www.mirbase.org webcite), this despite the economic importance of barley and its role as model species for Triticeae [4].

The conservation of miRNA sequences across species provides a powerful tool for the identification of novel miRNA genes based on homology with miRNAs previously described in other species. Search based on evolutionary conservation has allowed the identification of miRNA families in many plant species, including those where the complete genome sequence is not available, as it is currently the case of barley. Without genome sequence information a powerful alternative data source comes from ESTs (Expressed Sequence Tags): currently 501,616 ESTs are available in barley http://www.ncbi.nlm.nih.gov/dbEST/dbEST_summary.html webcite[5].

The identification of target genes is a fundamental step for the determination of the biological function of microRNAs, besides being an indirect evidence for their existence. Evolutionary conserved targets have proven very helpful to test the effectiveness of miRNA target detection. The perfect or near perfect complementarity between a miRNA and its target mRNA, that is a peculiar feature of plant miRNAs, gives a powerful tool for the identification of target genes through BLAST analysis of miRNA mature sequences vs EST/genomic sequences. A large part of the "in silico" predicted targets have then been confirmed as bona fide targets by experimental approaches including Northern, 5'-RACE and, more recently, degradome analysis via NGS (Next Generation Sequencing) [6,7].

The correct binding of miRNA to its cognate mRNA is critical for regulating the mRNA level and protein expression. This binding can be affected by single-nucleotide polymorphisms or indels in the miRNA target site leading to the suppression of existing binding sites or the generation of illegitimate ones. Therefore, small polymorphisms in miRNA targets can have a relevant effect on gene and protein expression and represent a type of genetic variability that can influence agronomical traits. As an example, overexpression of miR156b and miR156h in rice results in severe dwarfism, strongly reduced panicle size and delayed heading date [8].

To extend and to update information about miRNAs and their targets in barley and to identify candidate polymorphisms at microRNA target sites, barley EST sequences have been screened and related to Unigene clusters. UniGene is an experimental system for partitioning transcript sequences into a non-redundant set of gene-oriented clusters. Thus each UniGene cluster contains sequences that appear to come from the same transcription locus (gene or expressed pseudogene) http://www.ncbi.nlm.nih.gov/UniGene/index.html webcite. Mining SNPs from ESTs allows the exploitation of genetic variability based on published sequences and the analysis of Unigene clusters can be very helpful for this purpose [9].

Results and Discussion

Barley miRNAs

Since only mature miRNA sequences rather than precursor sequences are conserved among plant species, mature miRNA sequences have been used as queries for BLAST search against Hordeum vulgare ESTs [10]. One hundred-fifty-six microRNA mature sequences belonging to 50 miRNA families have been found to significantly match at least one EST sequence in barley (the total number of matching ESTs was 855 - as reported in additional files 1 and 2) and could actually be related both to target or miRNA sequences, even if the probability is lower for the latter. Indeed the estimated frequency of pri-miRNAs in T. aestivum EST collection is as low as 0.003% [2].

Additional file 1. BLAST results. Alignments identified between plant miRNAs and barley ESTs.

Format: TXT Size: 613KB Download fileOpen Data

Additional file 2. Barley mature miRNA sequences. Aligned barley mature miRNA sequences grouped on the basis of the miRNA family. Families for which only one barley EST has been found to match are not reported in this file.

Format: DOC Size: 198KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

The results illustrated above have been compared with those reported by Dryanova et al. where miRNAs and their targets have been searched in the Triticeae tribe [2]. Among the 33 miRNA families identified by Dryanova et al. in at least one species of the Triticeae tribe, 22 families were found in barley and 17 of them overlap with the present findings. Regarding barley, some miRNA families were found in just one of the two papers. Dryanova et al. found evidences for 5 additional miRNA families while the present work has found evidences in barley for miR390 and miR396 previously reported only in T. aestivum, and for additional 31 families not found by Dryanova et al. in anyone of the investigated species (i.e. miR442, miR529). The reasons for these discrepancies can be ascribed to the different miRBase release used (miRBase Release 8.0 for Dryanova et al., 2008 and miRBase Release 13.0 in the present work) and partially to differences in the BLAST settings adopted. Monocot-specific miRNAs (i.e. miR444) have also been found in both works [11].

Statistical analysis was employed to identify over and under-represented plant species from which the corresponding barley miRNA comes from. As reported in table 1 and 2, barley miRNA sequences putative orthologous to those of Triticum are significantly over-represented in our data also when very stringent p-value, e.g. 0.001, was used. Hordeum and Triticum genera are both members of the Poaceae family, Pooideae subfamily, Triticeae tribe. H. vulgare is often used as a model species for Triticeae, thanks to its diploid genome that could facilitate genome-wide searches of miRNAs.

Table 1. Statistical analysis for the identification of over and under-represented plant species.

Table 2. Over and under-represented plant species within barley miRNAs identified with respect to the stringency chosen for the p-value.

Zea mays is also closely related to barley being part of monocot group and Poaceae family. Oryza sativa although is part of Poaceae family is under-represented, when a low stringent p-value (0.05) was used.

Some ESTs have matched to more than one miRNAs belonging either to the same family or to different families (additional file 3). The first case can be due to the high level of similarity among mature sequences from different members of the same family, while ESTs matching to different miRNA families could represent examples of multi-microRNA based control.

Additional file 3. microRNA matching ESTs. List of ESTs matching the mature miRNA sequences. ESTs that present more than 1 target site for the same miRNA family or for different miRNA families have been highlighted.

Format: XLS Size: 238KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Transcripts targeted by more than one miRNA have also been found also in other plant species such as rice [12]. These findings are common in animals where many different miRNAs recognize the same target mRNA, usually at the 3'UTR [13].

To identify and annotate potential microRNA-regulated genes in barley, the 855 matching ESTs were related to Unigene clusters. Clusters annotated as protein-coding sequences were then selected for subsequent analysis and listed in tables 3 and 4. A total of 121 different Unigene clusters putatively representing the targets for 37 miRNA families has been found. Similar results (e.g. on average more than 1 putative target/miRNA family) were reported by Zhang et al. in maize (115 target for 26 miRNA families) [14]. Sometimes different targets for a specific miRNA are members of the same gene family (e.g. miR156-SBP family), while in other cases there is no evident relationship among the putative targets of a given miRNA (e.g. miR1121). Previous studies report six targets or fewer for most Arabidopsis miRNAs, a number significantly lower than in animals, for example, in Drosophila each miRNA has on average over 50 predicted targets [13,15].

Table 3. miRNA target genes identified in barley and confirmed by previous studies

Table 4. Novel miRNA target genes identified

Although several of the candidate miRNA/target pairs here identified have the same functional annotation reported in previously studied species (table 3) and specifically in barley some putative novel microRNA/target pairs have been discovered (table 4) [2]. Actually, some of these novel targets were reported by literature as regulated by a different microRNA. Most of the novel miRNA/target pairs refer to miRNAs recently discovered and thus probably less studied (i.e. miR1120, miR1122, miR1134). The Argonaute-like protein found as a novel target for miR408 in H.vulgare by Dryanova et al. has been confirmed also in the present work.

Transcription factor families comprise most of the highly conserved miRNA targets (see table 3) such as SBP family for miRNA 156, AP2 family for miR172, GRAS family for miR171, myb family for miR159, GRF family for miR396 and ARF family for miR160. These results confirmed what previously observed in Triticeae and in other species [2]. In rice about 70% of conserved miRNA targets are transcription factors, while in wheat one-third of the predicted targets was found to encode for transcription factors [11,12]. Conserved miRNAs also target genes involved in their own biogenesis and function: as an example miR168 targets AGO1 which is part of the RISC complex responsible for the miRNA-mediated mRNA cleavage [15]. miRNA regulate gene expression also by targeting enzymes of the ubiquitination pathway (ubiquitin conjugating enzyme E2 and TIR1/ubiquitin ligase): barley miR393, miR399, miR1128, miR1133, miR1135 can be considered putative regulators of gene expression at protein level.

The number of target genes identified as different Unigene clusters (tables 3-4) is very different among the miRNA families. In rice Zhou et al. have found a high number of targets for miR156 and miR396 and a low number for miR162, miR167, miR395, miR398 and miR399 [12]. This finding could indicate that the former miRNAs are nodes in gene regulation networks, while the latter could act on specialized pathways.

The predicted targets have been grouped into functional categories and reported in figures 1 and 2 where the target annotations based on GO terms are shown. Biological processes known to be regulated by miRNAs, such as development and response to biotic and abiotic stress, have been highlighted both in known (figure 1a) and in novel targets (figure 2a). Moreover, most of the molecular functions are related to transcriptional regulation and DNA/nucleotide binding in both groups (figures 1b-2b). These findings suggest that the predicted target genes can be considered a reliable dataset to be used in subsequent analysis.

thumbnailFigure 1. Functional enrichment for the miRNA targets identified. For each GO term it is shown the number of targets annotated with that term with respect to the total number of targets (%). Figure 1a refers to the biological process, while figure 1b refers to the molecular function.

thumbnailFigure 2. Functional enrichment for the novel miRNA targets identified. For each GO term it is shown the number of targets annotated with that term with respect to the total number of novel targets (%). Figure 2a refers to the biological process, while figure 2b refers to the molecular function.

For some Unigene clusters the annotation was related to transcribed genes rather than protein coding sequences. These Unigenes could represent miRNA-coding genes as shown by other authors [16,17]. Table 5 reports the Unigene clusters candidate to encode miRNA coding genes on the basis of the precursor sequence secondary structure (MFEI >0.85, see Materials and Methods) and of the presence of the miRNA* (miRNA passenger sequence). It cannot be excluded that the clusters unable to fold with a miRNA-like structure (e.g. Hv.8579, Hv.11623) are false negatives for several reasons, such as truncated precursor sequences in EST database. Putative microRNA sequences have also been BLASTed against previously known precursors available from mirBASE: the analysis found similarities with 6 different miRNA families. The secondary structures of the putative microRNA precursors are reported in the additional file 4. Linking together sequences containing miRNA precursors from Dryanova et al. and from the present work, information on several microRNA putative secondary structures, belonging to 10 miRNA families are now available [2]. The mature miRNAs predicted from these data are 18 to 24 nt long, with a higher frequency for 20 and 21 nt.

Table 5. Unigene clusters candidate to encode for miRNAs

Additional file 4. Secondary structures of the putative miRNA precursors. Predicted structures of the identified putative miRNA precursors. The prediction was performed with mfold.

Format: DOC Size: 407KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Genetic variation at miRNA target sites

A single nucleotide change in the sequence of a target site can affect miRNA regulation: as a consequence naturally occurring SNPs in target sites are candidates for relevant functional variations. Nair et al. established a perfect association between a SNP at the miR172 targeting site and cleistogamy in barley [18]. Overall few papers have been published to date describing variations among plant genotypes at miRNAs and their target sites, while plenty of information is available for humans [19-23]. Genome-wide studies in humans have shown that the levels of polymorphism at miRNA and miRNA target sites are lower than at coding or neutral regions, however beneficial miRNA-target site polymorphisms also exist [19].

In this study, publicly available SNP data have been analyzed in context with miRNAs and their target sites. EST-derived SNPs can provide a rich source of biologically useful genetic variation due to the redundancy of gene sequence, the diversity of genotypes present in the databases and the fact that each putative polymorphism is associated with an expressed gene. Variations both in functional regions of putative miRNAs (mature sequence) and at miRNA target sites have been found. Previous works in human have highlighted a relatively low level of variation in functional microRNA regions and an appreciable level of variation at target sites [21].

Hv.5064, the candidate for miR1137 coding sequence, has been tested for modifications of pre-miRNA structure due to a base substitution in position 13 (C/G, table 6, figure 3). To evaluate the possible impact of this SNP on pre-miRNA secondary structure, Gibbs free energy (ΔG) and MFEI from each version of pre-miRNA were calculated using mfold program. Data in figure 3 show the structural variation obtained when moving from "C variant" to "G variant" with a higher MFEI for the second one and thus a greater stability of the molecule (miRNA-miRNA* pairing enhanced in the G variant). Difference in ΔG moving from C to G and vice versa were calculated according to Ehrenreich and Purugganan [19]. ΔΔG was +1.3 for the former change and -1.3 for the latter suggesting that some SNPs can stabilize/destabilize pre-miRNA structure. No target gene has been reported in literature for miR1137.

Table 6. Putative polymorphisms identified at miRNA target sites and inside miRNA mature sequences

thumbnailFigure 3. Predicted secondary structures of the two variants of the miR1137 precursor identified in the Unigene cluster Hv.5064. The variants with a C and a G in the 13th position are respectively reported in the left and in the right side of the figure. The table shows for each variant: the free energy ΔG, the length of the precursor, the GC content, the MFEI (Minimal Folding Energy Index), the number of mismatches between the mature sequence and the paired miRNA* passenger and the arm of the hairpin where the mature sequence is located.

In plants most of the miRNA-based regulation relies on the cleavage of target mRNAs that normally occurs at the tenth nucleotide of the complementary region and numerous studies on miRNA-target interaction have highlighted the importance of positions 2 to 12, more frequently 10 and 11 [24]. Although most of the putative polymorphisms highlighted in this work are outside those critical positions, several examples of putative functionally relevant polymorphisms have been detected. Table 6 reports the putative polymorphisms detected after comparison among EST sequences inside Unigene clusters, without any selection against false positives. Some of these nucleotide variation could be due to sequencing errors or related to very similar genes belonging to a specific family, nevertheless when the SNPs/indels rely on two or more copies of independent sequences it can be considered a good candidate for a true positive polymorphic target site [25]. For example, a polymorphism in miRNA 408 target site detected by AutoSNP in contig 2094 (coding for a plastocyanin) is based on sequences from two different cultivars reporting the same allelic variant as part of a haplotype where a SSR (Simple Sequence Repeat) polymorphism is located upstream the target sequence (figure 4). Some polymorphisms also showed an evolutionary conserved position, the nucleotide variation identified in Hv.2498 (targeted by miR393) has also been found in the orthologous gene of Arabidopsis in the same position by Ehrenreich and Purugganan [19].

thumbnailFigure 4. SNP identified in contig2094 within the target site for miRNA 408. In this multiple alignment performed with AutoSNP, two cultivars (Optic and Morex) report the same allelic variant as part of a haplotype including a SSR polymorphism located upstream the target sequence.

The Squamosa-promoter Binding Protein (SBP) is a known target family for miR156. Many plant transcription factors involved in the regulation of the transition from the vegetative to the reproductive phase belong to this family and it has been shown that overexpressing SBP genes can lead to increased leaf initiation, decreased apical dominance and delayed flowering time [15]. The increase of the activity of some miRNAs (among which miR156) is part of the infection strategy performed by the Turnip mosaic virus in Arabidopsis [26,27]. miR156 performs a critical function in mediating developmental processes and it is also related to the response to biotic stress. The screening of barley databases has identified two SBP genes targeted by miR156 for which two nucleotide variations occur in critical positions (11-12). If these SNPs will be experimentally confirmed, they could have the effect of destabilizing the interaction between the miRNA and the mRNA, which could consequently avoids cleavage and lead to phenotypical variations in developmental features or in the resistance to viral infection.

A SNP also occurs in a crucial point of the experimentally confirmed NAC1 target for miR164. NAC1 is a transcription factor involved in shoot apical meristem formation and auxin-mediated lateral root formation. Guo et al. showed that the overexpression of miR164 leads to reduced lateral rooting; conversely the disruption of the regulation mediated by this miRNA increases the number of lateral roots [28]. The authors have reported that miR164 directs cleavage in vivo at a position complementary to the 10th nucleotide from the 5' end of the mature sequence [28]. The SNP found in barley is in the 11th position, therefore it is likely to prevent the cleavage and produce phenotypic effects on root development.

SNPs have been identified also in other two conserved miRNA targets, TIR1 and AGO4, targeted respectively by miR393 and miR408. TIR1 is an auxin-receptor negatively regulated by miRNAs in response to bacterial flagellin, as a defence mechanism against Pseudomonas syringae [29]. AGO4 is a protein involved in the siRNA mediated gene silencing, and it is required for the resistance to the same pathogen [30]. Therefore, miR393 and miR408 are likely to work in a coupled manner during P. syringae infection. The two SNPs identified are in the 12th position and could potentially alter the levels of pathogens resistance.

SNPs were also found in previously not reported miRNA targets, such as the AWPM-19-like protein matching to the miRNA 1134. AWPM-19 accumulates in wheat plasma membrane during cold acclimation in response to abscisic acid [31]. If this miRNA really controls the synthesis of this protein, a deleterious SNP in the 11th position could then change resistance to cold stress.

Conclusions

This study has thus provided an update of the information on barley miRNAs and their targets representing a foundation for future studies.

Novel putative target genes have been identified and most of them are involved in stress and hormone response. Indeed, the role of plant miRNAs in abiotic and biotic stress response as well as in auxin signalling is well known [32,33]. In particular, protein kinases such as protein kinase C and serine/threonine kinase, known to be important regulator on abiotic stress resistance, are largely present in novel microRNA/target pairs identified.

The results have also shown that microRNA target sites can be an interesting source for the identification of functional genetic variability, representing an interesting source of candidate molecular markers for application in barley breeding. Putative polymorphisms have now to be verified by amplification and sequencing of the target sequences on a larger set of genotypes.

Sequence analysis based on known miRNAs can obviously give insights only on conserved mRNAs and related targets. Future work will thus be based on the construction of a degradome library for parallel analysis of RNA end (PARE), a powerful approach for high-throughput identification/validation of conserved and non conserved targets.

Methods

miRNA reference dataset

The initial miRNA dataset has been obtained by extracting the mature sequences (1929 entries) of the Viridiplantae group from the miRBase release 13 http://www.mirbase.org webcite[34]. By removing identical mature sequences, the size of this dataset has been subsequently reduced to 1014 non-redundant sequences related to 468 miRNA families.

Searching for mature miRNAs matching sequences in barley

The full collection of non-redundant mature miRNA sequences was used in a BLASTn search against dbEST http://www.ncbi.nlm.nih.gov webcite, accepting a number of mismatch lower than 4.

The set of miRNA mature sequences (including the identical sequences removed at the first step of the work) with at least one matching EST have been classified on the basis of the species of origin. The binomial distribution was used to assess the statistical significance for the represented plant species; this allowed identifying those species chosen from the initial dataset more or less frequently than random. Four different thresholds for the p-values were applied (0.05, 0.01, 0.005, 0.001).

Matching ESTs have then been related to Unigene clusters and the corresponding annotations were recorded (if available). The GO slimmer tool available on the Gene Ontology website http://www.geneontology.org webcite has been used to identify the GO slim terms more represented in the set of potential targets on the basis of the Unigene cluster annotations. For this analysis the Plant GO Slim subset has been used.

Identification of putative miRNA precursors

True miRNA precursors should have both a mature sequence on one arm of the hairpin and a paired passenger sequence (called miRNA*) on the opposite arm. To assess these features the precursor sequences were extracted from the consensus sequences, obtained by the Sequencer Software (Gene Codes) on Unigene cluster assemblies, by cutting 13 nt before the 5' hit and 13 nt after the 3' hit, since this region (called the pri-extension region of the hairpin) was recently shown to have this average length in plants [35]. In order to predict the secondary structure of the precursors, the software mfold 3.2, free available at http://mfold.bioinfo.rpi.edu/cgi-bin/rna-form1.cgi webcite, was used [36]. The minimal folding free energy index (MFEI) and the GC content were calculated for each sequence.

All the sequences with a MFEI greater than 0.85 were considered potential miRNA precursors [37]; besides, only 4 mismatches were allowed between the mature sequence and the passenger sequence, and only few and small asymmetric bulges were accepted [38].

Identification of SNPs/indels at miRNA target sites

Polymorphisms in target genes have been searched through a comparison of the ESTs belonging to the same Unigene cluster. Each cluster has been assembled by Sequencer Software (Gene Codes) and polymorphisms have been searched on miRNA complementary sequence sites.

AutoSNP database http://autosnpdb.qfab.org.au webcite was also screened using target gene annotations as contig-searching keywords.

Authors' contributions

MC identified barley microRNAs and targets and highlighted the putative novel ones. He was also involved in polymorphism searching. PF conceived, designed and coordinated the work and wrote the manuscript. AS gave a contribution for polymorphisms searching and alignment of barley mature microRNA sequences. PB and AL prepared the additional material and gave a contribution to bioinformatics analysis. LC contributed in the design and discussion of the work and assisted in drafting the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work has been supported by "MIRNA" project and "Mappa 5A" project (Italian Ministry for Agriculture). The authors thank Mr. Renzo Alberici and Miss Alice Martini for the technical assistance provided.

References

  1. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism and function.

    Cell 2004, 116:281-297. PubMed Abstract | Publisher Full Text OpenURL

  2. Dryanova A, Zakharov A, Gulick PJ: Data mining for miRNAs and their targets in the Triticeae.

    Genome 2008, 51:433-443. PubMed Abstract | Publisher Full Text OpenURL

  3. Sunkar R, Jagadeeswaran G: In silico identification of conserved microRNAs in large number of diverse plant species.

    BMC Plant Biol 2008, 8:37. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  4. Schreiber AW, Sutton T, Caldo RA, Kalashyan E, Lovell B, Mayo G, Muehlbauer GJ, Druka A, Waugh R, Wise RP, Langridge P, Baumann U: Comparative transcriptomics in the Triticeae.

    BMC Genomics 2009, 10:285. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  5. Zhang BH, Pan XP, Wang QL, Cobb GP, Anderson TA: Identification and characterization of new plant microRNAs using EST analysis.

    Cell Res 2005, 15(5):336-360. PubMed Abstract | Publisher Full Text OpenURL

  6. German MA, Pillay M, Jeong D, Hetawal A, Luo S, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R, De Paoli E, Lu C, Schroth G, Meyers BC, Green PJ: Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends.

    Nat Biotechnol 2008, 26(8):941-946. PubMed Abstract | Publisher Full Text OpenURL

  7. Zeng C, Wang W, Zheng Y, Chen X, Bo W, Song S, Zhang W, Peng M: Conservation and divergence of microRNAs and their functions in Euphorbiaceous plants.

    Nucleic Acids Res 2009, 38(3):981-995. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Xie K, Wu C, Xiong L: Genomic organization, differential expression and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice.

    Plant Physiol 2006, 142:280-293. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Picoult-Newberg L, Ideker TE, Pohl MG, Taylor SL, Donaldson MA, Nickerson DA, Boyce-Jacino M: Mining SNPs from EST databases.

    Genome Res 1999, 9:167-174. PubMed Abstract | PubMed Central Full Text OpenURL

  10. Zhang B, Pan X, Cannon CH, Cobb GP, Anderson TA: Conservation and divergence of plant microRNA genes.

    Plant J 2006, 46:243-259. PubMed Abstract | Publisher Full Text OpenURL

  11. 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:R96. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  12. Zhou M, Gu L, Li P, Song X, Wei L, Chen Z, Cao X: Degradome sequencing reveals endogenous small RNA targets in rice (Oryza sativa L., ssp. Indica).

    Front Biol 2010, 5(1):67-90. OpenURL

  13. Grun D, Wang YL, Langenberger D, Gunsalus KC, Rajewsky N: microRNA target predictions across seven Drosophila species and comparison to mammalian targets.

    PLoS Comput Biol 2005, 1:e13. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Zhang B, Pan X, Anderson TA: Identification of 188 conserved maize microRNAs and their targets.

    FEBS Lett 2006, 580:3753-3762. PubMed Abstract | Publisher Full Text OpenURL

  15. Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAs and their regulatory roles in plants.

    Annu Rev Plant Biol 2006, 57:19-53. PubMed Abstract | Publisher Full Text OpenURL

  16. Guo Q, Xiang A, Yang Q, Yang Z: Bioinformatic identification of microRNAs and their target genes from Solanum tuberosum expressed sequence tags.

    Chinese Sci Bull 2007, 52(17):2380-2389. Publisher Full Text OpenURL

  17. Xie FL, Huang SQ, Guo K, Xiang AL, Zhu YY, Nie L, Yang ZM: Computational identification of novel microRNAs and targets in Brassica napus.

    FEBS Lett 2007, 581:1464-1474. PubMed Abstract | Publisher Full Text OpenURL

  18. Nair SK, Wang N, Turuspekov Y, Pourkheirandish M, Sinsuwongwat S, Chen G, Sameri M, Tagiri A, Honda I, Watanabe Y, Kanamori H, Wicker T, Stein N, Nagamura Y, Matsumoto T, Komatsuda T: Cleistogamous flowering in barley arises from the suppression of microRNA-guided HvAP2 mRNA cleavage.

    PNAS 2010, 107:480-496. OpenURL

  19. Ehrenreich IM, Purugganan MD: Sequence variation of microRNAs and their binding sites in Arabidopsis.

    Plant Physiol 2008, 146:1974-1982. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Guo X, Gui Y, Wang Y, Zhu Q, Helliwell C, Fan L: Selection and mutation on microRNA target sequences during rice evolution.

    BMC Genomics 2008, 9:454. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  21. Saunders MA, Liang H, Li WH: Human polymorphisms at microRNAs and microRNA target sites.

    PNAS 2007, 104:3300-3305. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Sethupathy P, Collins FS: MicroRNA target site polymorphisms and human disease.

    Trends Genet 2008, 24:489-497. PubMed Abstract | Publisher Full Text OpenURL

  23. Gardner PP, Vinther J: Mutation of miRNA target sequences during human evolution.

    Trends Genet 2008, 24:262-265. PubMed Abstract | Publisher Full Text OpenURL

  24. Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant microRNA targets.

    Cell 2002, 110:513-520. PubMed Abstract | Publisher Full Text OpenURL

  25. Batley J, Barker G, O'Sullivan H, Edwards KJ, Edwards D: Mining for single nucleotide polymorphisms and insertions/deletions in maize expressed sequence tag data.

    Plant Physiol 2003, 132:84-91. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Kasschau KD, Xie Z, Allen E, Llave C, Chapman EJ, Krizan KA, Carrington JC: P1/HC-Pro, a viral suppressor of RNA silencing, interferes with Arabidposis development and miRNA function.

    Dev Cell 2003, 4:205-217. PubMed Abstract | Publisher Full Text OpenURL

  27. Chen J, Li WX, Xie D, Peng JR, Ding SW: Viral virulence protein suppresses RNA silencing-mediated defense but upregulates the role of microRNA in host gene expression.

    Plant Cell 2004, 16:1302-1313. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Guo HS, Xie Q, Fei JF, Chua NH: MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development.

    Plant Cell 2005, 17:1376-1386. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, Voinnet O, Jones JDG: A plant miRNA contributes to antibacterial resistance by repressing auxin signalling.

    Science 2006, 312:436-439. PubMed Abstract | Publisher Full Text OpenURL

  30. Agorio A, Vera P: ARGONAUTE4 is required for resistance to Pseudomonas syringae in Arabidopsis.

    Plant Cell 2007, 19:3778-3790. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Koike M, Takezawa D, Arakawa K, Yoshida S: Accumulation of 19-kDa plasma membrane polypeptide during induction of freezing tolerance in wheat suspension-cultured cells by abscisic acid.

    Plant Cell Physiol 1997, 38:707-716. PubMed Abstract | Publisher Full Text OpenURL

  32. Martin RC, Liu PP, Goloviznina NA, Nonogaki H: microRNA, seeds, and Darwin?: diverse function of microRNA in seed biology and plant response to stress.

    J Exp Bot 2010, 69:2229-2234. Publisher Full Text OpenURL

  33. Teotia PS, Mukherjee SK, Mishra NS: Fine tuning of auxin signalling by miRNAs.

    Physiol Mol Biol Plants 2008, 14:81-90. Publisher Full Text OpenURL

  34. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics.

    NAR 2008, (36 Database):D154-D158. PubMed Abstract | PubMed Central Full Text OpenURL

  35. Kadri S, Hinman V, Benos PV: HHMMiR: efficient de novo prediction of microRNAs using hierarchical hidden Markov models.

    Bioinformatics 2009, 10(Suppl I):S35. PubMed Abstract | PubMed Central Full Text OpenURL

  36. Zuker M: mfold web server for nucleic acid folding and hybridization prediction.

    Nucleic Acids Res 2003, 31(13):3406-15. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Zhang BH, Pan XP, Cox SB, Cobb GP, Anderson TA: Evidence that miRNAs are different from other RNAs.

    Cell Mol Life Sci 2006, 63:246-254. PubMed Abstract | Publisher Full Text OpenURL

  38. Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Vonnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant microRNAs.

    Plant Cell 2008, 20:3186-3190. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA.

    Molecular Cell 2004, 14:787-799. PubMed Abstract | Publisher Full Text OpenURL

  40. Puzey JR, Kramer EM: Identification of conserved Aquilegia coerulea microRNAs and their targets.

    Gene 2009, 448(1):46-56. PubMed Abstract | Publisher Full Text OpenURL

  41. Unver T, Budak H: Conserved microRNAs and their targets in model grass species Brachypodium distachyon.

    Planta 2009, 230(4):659-69. PubMed Abstract | Publisher Full Text OpenURL

  42. 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.

    Current Biology 2002, 12(17):1484-1495. PubMed Abstract | Publisher Full Text OpenURL

  43. Trindade I, Capitão C, Dalmay T, Fevereiro MP, Dos Santos DM: miR398 and miR408 are up-regulated in response to water deficit in Medicago truncatula.

    Planta 2010, 231(3):705-716. PubMed Abstract | Publisher Full Text OpenURL

  44. Lu S, Sun YH, Shi R, Clark C, Li LG, Chiang VL: Novel and mechanical stress-responsive microRNAs in Populus trichocarpa that are absent from Arabidopsis.

    Plant Cell 2005, 17:2186-2203. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Li YF, Zheng Y, Addo-Quaye C, Zhang L, Saini A, Jagadeeswaran G, Axtell MJ, Zhang WX, Sunkar R: Transcriptome-wide identification of microRNA targets in rice.

    Plant J 2010, 62(5):742-759. PubMed Abstract | Publisher Full Text OpenURL

  46. Zhang LF, Chia JM, Kumari S, Stein JC, Liu ZJ, Narechania A, Maher CA, Guill K, McMullen MD, Ware D: A genome-wide characterization of microRNA genes in maize.

    PLoS Genet 2009, 5(11):e1000716. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana.

    Genes and Development 2006, 20(24):3407-3425. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Lacombe S, Nagasaki H, Santi C, Duval D, Piégu B, Bangratz M, Breitler JC, Guiderdoni E, Brugidou C, Hirsch J, Cao XF, 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL