Computational prediction of candidate miRNAs and their targets from the completed Linum ussitatissimum genome and EST database

Flax is an important agronomic crop grown for its fiber (linen) and oil (linseed oil). In spite of many thousands of years of breeding some fiber varieties have been shown to rapidly respond to environmental stress with heritable changes to its genome. Many miRNAs appear to be induced by abiotic or biotic conditions experienced through the plant life cycle. Computational miRNA analysis of the flax genome provides a foundation for subsequent research on miRNA function in Linum usitatissimum and may also provide novel insight into any regulatory role the RNAi pathway may play in generating adaptive structural variation in response to environmental stress. Here a bioinformatics approach is used to screen for miRNAs previously identified in other plant species, as well as to predict putative miRNAs unique to a particular species which may not have been identified as they are less abundant or dependent upon a specific set of environmental conditions. Twelve miRNA genes were identified in flax on the basis of unique pre-miRNA positions with structural homology to plant pre-miRNAs and complete sequence homology to published plant miRNAs. These miRNAs were found to belong to 7 miRNA families, with an additional 2 matches corresponding to as yet unnamed poplar miRNAs and a parologous miRNA with partial sequence homology to mtr-miR4414b. An additional 649 novel and distinct flax miRNA genes were identified to form from canonical hairpin structures and to have putative targets among the ~30,000 flax Unigenes.