Identification of novel and candidate miRNAs in rice by high throughput sequencing

microRNA, seeds, and Darwin?: diverse function of miRNA in seed biology and plant responses to stress

microRNAs (miRNAs) are small, single-stranded RNAs that down-regulate target genes at the post-transcriptional level. miRNAs regulate target genes by guiding mRNA cleavage or by repressing translation. miRNAs play crucial roles in a broad range of developmental processes in plants. Multiple miRNAs are present in germinating seeds and seedlings of Arabidopsis, some of which are involved in the regulation of germination and seedling growth by plant hormones such as abscisic acid (ABA) and auxin. The involvement of miRNAs in ABA responses is not limited to the early stages of plant development but seems to be important for general stress responses throughout the plant life cycle. This Darwin review summarizes recent progress in miRNA research focusing on seed and stress biology, two topics which were of interest to Charles Darwin.

Discovery and comparative profiling of microRNAs in a sweet orange red-flesh mutant and its wild type

Background

Red-flesh fruit is absent from common sweet orange varieties, but is more preferred by consumers due to its visual attraction and nutritional properties. Our previous researches on a spontaneous red-flesh mutant revealed that the trait is caused by lycopene accumulation and is regulated by both transcriptional and post-transcriptional mechanisms. However, the knowledge on post-transcriptional regulation of lycopene accumulation in fruits is rather limited so far.

Results

We used Illumina sequencing method to identify and quantitatively profile small RNAs on the red-flesh sweet orange mutant and its wild type. We identified 85 known miRNAs belonging to 48 families from sweet orange. Comparative profiling revealed that 51 known miRNAs exhibited significant expression differences between mutant (MT) and wild type (WT). We also identified 12 novel miRNAs by the presence of mature miRNAs and corresponding miRNA*s in the sRNA libraries. Comparative analysis showed that 9 novel miRNAs are differentially expressed between WT and MT. Target predictions of the 60 differential miRNAs resulted 418 target genes in sweet orange. GO and KEGG annotation revealed that high ranked miRNA-target genes are those implicated in transcription regulation, protein modification and photosynthesis. The expression profiles of target genes involved in carotenogenesis and photosynthesis were further confirmed to be complementary to the profiles of corresponding miRNAs in WT and MT.

Conclusion

This study comparatively characterized the miRNAomes between the red-flesh mutant and the wild type, the results lay a foundation for unraveling the miRNA-mediated molecular processes that regulate lycopene accumulation in the sweet orange red-flesh mutant.

DNA methylation mediated by a microRNA pathway

In plants, the known microRNAs (miRNAs) are produced as approximately 21 nucleotide (nt) duplexes from their precursors by Dicer-like 1 (DCL1). They are incorporated into Argonaute 1 (AGO1) protein to regulate target gene expression primarily through mRNA cleavage. We report here the discovery of a class of miRNAs in the model monocot rice (Oryza sativa). These are 24 nt in length and require another member of the Dicer family, DCL3, for their biogenesis. The 24 nt long miRNAs (lmiRNAs) are loaded into AGO4 clade proteins according to hierarchical rules, depending on the upstream biogenesis machinery and the 5'-terminal nucleotide. We demonstrated that lmiRNAs direct DNA methylation at loci from which they are produced as well as in trans at their target genes and play roles in gene regulation. Considered together, our findings define a miRNA pathway that mediates DNA methylation.

MicroRNA gene evolution in Arabidopsis lyrata and Arabidopsis thaliana

MicroRNAs (miRNAs) are short regulatory RNAs processed from partially self-complementary foldbacks within longer MIRNA primary transcripts. Several MIRNA families are conserved deeply through land plants, but many are present only in closely related species or are species specific. The finding of numerous evolutionarily young MIRNA, many with low expression and few if any targets, supports a rapid birth-death model for MIRNA evolution. A systematic analysis of MIRNA genes and families in the close relatives, Arabidopsis thaliana and Arabidopsis lyrata, was conducted using both whole-genome comparisons and high-throughput sequencing of small RNAs. Orthologs of 143 A. thaliana MIRNA genes were identified in A. lyrata, with nine having significant sequence or processing changes that likely alter function. In addition, at least 13% of MIRNA genes in each species are unique, despite their relatively recent speciation (approximately 10 million years ago). Alignment of MIRNA foldbacks to the Arabidopsis genomes revealed evidence for recent origins of 32 families by inverted or direct duplication of mostly protein-coding gene sequences, but less than half of these yield miRNA that are predicted to target transcripts from the originating gene family. miRNA nucleotide divergence between A. lyrata and A. thaliana orthologs was higher for young MIRNA genes, consistent with reduced purifying selection compared with deeply conserved MIRNA genes. Additionally, target sites of younger miRNA were lost more frequently than for deeply conserved families. In summary, our systematic analyses emphasize the dynamic nature of the MIRNA complement of plant genomes.

Post-transcriptional gene regulation of salinity and drought responses by plant microRNAs

In the past few years, factors involved in abscisic acid signalling have been isolated and recognized as elements related to RNA metabolism, suggesting that post-transcriptional regulation of gene expression is required for abiotic stress responses. Some of these factors can be linked to the biogenesis of microRNAs (miRNAs), small RNA molecules that are important regulators of gene expression at the posttranscriptional level by repressing mRNA expression. Here, we review the role of miRNAs in stress responses, highlighting recent advances in elucidating the role of individual miRNAs and efforts to identify stress-responsive miRNAs at a genome-wide level in different model plants. Complete understanding of miRNA action depends on the identification of its target transcripts, and recent developments in miRNA research indicate that they will be uncovered in the near future.

Cloning and validation of novel miRNA from basmati rice indicates cross talk between abiotic and biotic stresses

Most of the physiological processes are controlled by the small RNAs in several organisms including plants. A huge database exists on one type of small RNA, i.e., microRNAs (miRs) identified from diverse species. However, the processes of data-mining of miRs in most of the species are still incomplete. Rice feeds the hungry trillions and hence understanding its developmental processes as well as its stress biology, which might be largely controlled by the small RNA pathways, is certainly a worthwhile task. Here, we report the cloning and identification of approximately 40 new putative miRs from local basmati rice variety in accordance to the annotation suggested by Meyers et al. (Plant Cell 20:3186-3190, 2008). About 23 sequences were derived from rice exposed to salt stress while 18 were derived from rice infected with tungro virus. A few of these putative miRs were common to both. Our data showed that at least two of these miRs were up-regulated in response to both abiotic and biotic stresses. The miR target predictions indicate that most of the putative miRs target specific metabolic processes. The up-regulation of similar miRs in response to two entirely different types of stresses suggests a converging functional role of miRs in managing various stresses. Our findings suggest that more rice miRs need to be identified and a thorough understanding of the function of such miRs will help unravel the mysteries of rice stress biology.

Correction: High throughput approaches reveal splicing of primary microRNA transcripts and tissue specific expression of mature microRNAs in Vitis vinifera

The version of this article published in BMC Genomics 2009, 10:558, contains data in Table 1 which are now known to be unreliable, and an illustration, in Figure 1, of unusual miRNA processing events predicted by these unreliable data. In this full-length correction, new data replace those found to be unreliable, leading to a more straightforward interpretation without altering the principle conclusions of the study. Table 1 and associated methods have been corrected, Figure 1 deleted, supplementary file 1 added, and modifications made to the sections "Deep sequencing of small RNAs from grapevine leaf tissue" and "Microarray analysis of miRNA expression". The editors and authors regret the inconvenience caused to readers by premature publication of the original paper.

Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development

BACKGROUND

In eukaryotes, microRNAs (miRNAs) have emerged as critical regulators of gene expression. The Silkworm (Bombyx mori L.) is one of the most suitable lepidopteran insects for studying the molecular aspects of metamorphosis because of its large size, availability of mutants and genome sequence. Besides, this insect also has been amply studied from a physiological and biochemical perspective. Deep sequencing of small RNAs isolated from different stages of silkworm is a powerful tool not only for measuring the changes in miRNA profile but also for discovering novel miRNAs.

RESULTS

We generated small RNA libraries from feeding larvae, spinning larvae, pupae and adults of B. mori and obtained approximately 2.5 million reads of 18-30 nt. Sequence analysis identified 14 novel and 101 conserved miRNAs. Most novel miRNAs are preferentially expressed in pupae, whereas more than 95% of the conserved miRNAs are dynamically regulated during different developmental stages. Remarkably, the miRNA-star (miR*) of four miRNAs are expressed at much higher levels than their corresponding miRNAs, and their expression profiles are distinct from their corresponding miRNA profiles during different developmental stages. Additionally, we detected two antisense miRNA loci (miR-263-S and miR-263-AS; miR-306-S and miR-306-AS) that are expressed in sense and antisense directions. Interestingly, miR-263 and miR-306 are preferentially and abundantly expressed in pupae and adults, respectively.

CONCLUSIONS

We identified 101 homologs of conserved miRNAs, 14 species-specific and two antisense miRNAs in the silkworm. Our results provided deeper insights into changes in conserved and novel miRNA and miRNA* accumulation during development.

Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development

BACKGROUND

In eukaryotes, microRNAs (miRNAs) have emerged as critical regulators of gene expression. The Silkworm (Bombyx mori L.) is one of the most suitable lepidopteran insects for studying the molecular aspects of metamorphosis because of its large size, availability of mutants and genome sequence. Besides, this insect also has been amply studied from a physiological and biochemical perspective. Deep sequencing of small RNAs isolated from different stages of silkworm is a powerful tool not only for measuring the changes in miRNA profile but also for discovering novel miRNAs.

RESULTS

We generated small RNA libraries from feeding larvae, spinning larvae, pupae and adults of B. mori and obtained approximately 2.5 million reads of 18-30 nt. Sequence analysis identified 14 novel and 101 conserved miRNAs. Most novel miRNAs are preferentially expressed in pupae, whereas more than 95% of the conserved miRNAs are dynamically regulated during different developmental stages. Remarkably, the miRNA-star (miR*) of four miRNAs are expressed at much higher levels than their corresponding miRNAs, and their expression profiles are distinct from their corresponding miRNA profiles during different developmental stages. Additionally, we detected two antisense miRNA loci (miR-263-S and miR-263-AS; miR-306-S and miR-306-AS) that are expressed in sense and antisense directions. Interestingly, miR-263 and miR-306 are preferentially and abundantly expressed in pupae and adults, respectively.

CONCLUSIONS

We identified 101 homologs of conserved miRNAs, 14 species-specific and two antisense miRNAs in the silkworm. Our results provided deeper insights into changes in conserved and novel miRNA and miRNA* accumulation during development.

Over-expression of miR172 causes loss of spikelet determinacy and floral organ abnormalities in rice (Oryza sativa)

BACKGROUND

Regulation of gene expression by microRNAs (miRNAs) plays a crucial role in many developmental and physiological processes in plants. miRNAs act to repress expression of their target genes via mRNA cleavage or translational repression. Dozens of miRNA families have been identified in rice, 21 of which are conserved between rice and Arabidopsis. miR172 is a conserved miRNA family which has been shown to regulate expression of APETALA2 (AP2)-like transcription factors in Arabidopsis and maize. The rice genome encodes five AP2-like genes predicted to be targets of miR172. To determine whether these rice AP2-like genes are regulated by miR172 and investigate the function of the target genes, we studied the effect of over-expressing two members of the miR172 family on rice plant development.

RESULTS

Analysis of miR172 expression showed that it is most highly expressed in late vegetative stages and developing panicles. Analyses of expression of three miR172 targets showed that SUPERNUMERARY BRACT (SNB) and Os03g60430 have high expression in developing panicles. Expression of miR172 was not inversely correlated with expression of its targets although miR172-mediated cleavage of SNB was detected by 5' rapid amplification of cDNA ends (RACE). Over-expression of miR172b in rice delayed the transition from spikelet meristem to floral meristem, and resulted in floral and seed developmental defects, including changes to the number and identity of floral organs, lower fertility and reduced seed weight. Plants over-expressing miR172b not only phenocopied the T-DNA insertion mutant of SNB but showed additional defects in floret development not seen in the snb mutant. However SNB expression was not reduced in the miR172b over-expression plants.

CONCLUSIONS

The phenotypes resulting from over-expression of miR172b suggests it represses SNB and at least one of the other miR172 targets, most likely Os03g60430, indicating roles for other AP2-like genes in rice floret development. miR172 and the AP2-like genes had overlapping expression patterns in rice and their expression did not show an obvious negative correlation. There was not a uniform decrease in the expression of the AP2-like miR172 target mRNAs in the miR172b over-expression plants. These observations are consistent with miR172 functioning via translational repression or with expression of the AP2-like genes being regulated by a negative feedback loop.

Carrot alternative oxidase gene AOX2a demonstrates allelic and genotypic polymorphisms in intron 3

Single nucleotide polymorphisms (SNPs) and insertion-deletions (InDels) are becoming important genetic markers for major crop species. In this study, we focus on variations at genomic level of the Daucus carota L. AOX2a gene. The use of gene-specific primers designed in exon regions on the boundaries of introns permitted to recognize intron length polymorphism (ILP) in intron 3 AOX2a by simple polymerase chain reaction (PCR) assays. The length of intron 3 can vary in individual carrot plants. Thus, allelic variation can be used as a tool to discriminate between single plant genotypes. Using this approach, individual plants from cv. Rotin and from diverse breeding lines and cultivars were identified that showed genetic variability by AOX2a ILPs. Repetitive patterns of intron length variation have been observed which allows grouping of genotypes. Polymorphic and identical PCR fragments revealed underlying high levels of sequence polymorphism. Variability was due to InDel events and intron single nucleotide polymorphisms (ISNPs), with a repetitive deletion in intron 3 affecting a putative pre-miRNA site. The results suggest that high AOX2a gene diversity in D. carota can be explored for the development of functional markers related to agronomic traits.

High throughput approaches reveal splicing of primary microRNA transcripts and tissue specific expression of mature microRNAs in Vitis vinifera

Background

MicroRNAs are short (~21 base) single stranded RNAs that, in plants, are generally coded by specific genes and cleaved specifically from hairpin precursors. MicroRNAs are critical for the regulation of multiple developmental, stress related and other physiological processes in plants. The recent annotation of the genome of the grapevine (Vitis vinifera L.) allowed the identification of many putative conserved microRNA precursors, grouped into multiple gene families.

Results

Here we use oligonucleotide arrays to provide the first indication that many of these microRNAs show differential expression patterns between tissues and during the maturation of fruit in the grapevine. Furthermore we demonstrate that whole transcriptome sequencing and deep-sequencing of small RNA fractions can be used both to identify which microRNA precursors are expressed in different tissues and to estimate genomic coordinates and patterns of splicing and alternative splicing for many primary miRNA transcripts.

Conclusion

Our results show that many microRNAs are differentially expressed in different tissues and during fruit maturation in the grapevine. Furthermore, the demonstration that whole transcriptome sequencing can be used to identify candidate splicing events and approximate primary microRNA transcript coordinates represents a significant step towards the large-scale elucidation of mechanisms regulating the expression of microRNAs at the transcriptional and post-transcriptional levels.

A genome-wide characterization of microRNA genes in maize

MicroRNAs (miRNAs) are small, non-coding RNAs that play essential roles in plant growth, development, and stress response. We conducted a genome-wide survey of maize miRNA genes, characterizing their structure, expression, and evolution. Computational approaches based on homology and secondary structure modeling identified 150 high-confidence genes within 26 miRNA families. For 25 families, expression was verified by deep-sequencing of small RNA libraries that were prepared from an assortment of maize tissues. PCR-RACE amplification of 68 miRNA transcript precursors, representing 18 families conserved across several plant species, showed that splice variation and the use of alternative transcriptional start and stop sites is common within this class of genes. Comparison of sequence variation data from diverse maize inbred lines versus teosinte accessions suggest that the mature miRNAs are under strong purifying selection while the flanking sequences evolve equivalently to other genes. Since maize is derived from an ancient tetraploid, the effect of whole-genome duplication on miRNA evolution was examined. We found that, like protein-coding genes, duplicated miRNA genes underwent extensive gene-loss, with approximately 35% of ancestral sites retained as duplicate homoeologous miRNA genes. This number is higher than that observed with protein-coding genes. A search for putative miRNA targets indicated bias towards genes in regulatory and metabolic pathways. As maize is one of the principal models for plant growth and development, this study will serve as a foundation for future research into the functional roles of miRNA genes.

Rice MicroRNA effector complexes and targets

MicroRNAs (miRNAs) are small silencing RNAs with regulatory roles in gene expression. miRNAs interact with Argonaute (AGO) proteins to form effector complexes that cleave target mRNAs or repress translation. Rice (Oryza sativa) encodes four AGO1 homologs (AGO1a, AGO1b, AGO1c, and AGO1d). We used RNA interference (RNAi) to knock down the four AGO1s. The RNAi lines displayed pleiotropic developmental phenotypes and had increased accumulation of miRNA targets. AGO1a, AGO1b, and AGO1c complexes were purified and further characterized. The three AGO1s all have a strong preference for binding small RNAs (sRNAs) with 5' U and have Slicer activity. We cataloged the sRNAs in each AGO1 complex by deep sequencing and found that all three AGO1s predominantly bound known miRNAs. Most of the miRNAs were evenly distributed in the three AGO1 complexes, suggesting a redundant role for the AGO1s. Intriguingly, a subset of miRNAs were specifically incorporated into or excluded from one of the AGO1s, suggesting functional specialization among the AGO1s. Furthermore, we identified rice miRNA targets at a global level. The validated targets include transcription factors that control major stages of development and also genes involved in a variety of physiological processes, indicating a broad regulatory role for miRNAs in rice.

miRExpress: analyzing high-throughput sequencing data for profiling microRNA expression

Background

MicroRNAs (miRNAs), small non-coding RNAs of 19 to 25 nt, play important roles in gene regulation in both animals and plants. In the last few years, the oligonucleotide microarray is one high-throughput and robust method for detecting miRNA expression. However, the approach is restricted to detecting the expression of known miRNAs. Second-generation sequencing is an inexpensive and high-throughput sequencing method. This new method is a promising tool with high sensitivity and specificity and can be used to measure the abundance of small-RNA sequences in a sample. Hence, the expression profiling of miRNAs can involve use of sequencing rather than an oligonucleotide array. Additionally, this method can be adopted to discover novel miRNAs.

Results

This work presents a systematic approach, miRExpress, for extracting miRNA expression profiles from sequencing reads obtained by second-generation sequencing technology. A stand-alone software package is implemented for generating miRNA expression profiles from high-throughput sequencing of RNA without the need for sequenced genomes. The software is also a database-supported, efficient and flexible tool for investigating miRNA regulation. Moreover, we demonstrate the utility of miRExpress in extracting miRNA expression profiles from two Illumina data sets constructed for the human and a plant species.

Conclusion

We develop miRExpress, which is a database-supported, efficient and flexible tool for detecting miRNA expression profile. The analysis of two Illumina data sets constructed from human and plant demonstrate the effectiveness of miRExpress to obtain miRNA expression profiles and show the usability in finding novel miRNAs.

PMRD: plant microRNA database

MicroRNAs (miRNA) are approximately 21 nucleotide-long non-coding small RNAs, which function as post-transcriptional regulators in eukaryotes. miRNAs play essential roles in regulating plant growth and development. In recent years, research into the mechanism and consequences of miRNA action has made great progress. With whole genome sequence available in such plants as Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, Glycine max, etc., it is desirable to develop a plant miRNA database through the integration of large amounts of information about publicly deposited miRNA data. The plant miRNA database (PMRD) integrates available plant miRNA data deposited in public databases, gleaned from the recent literature, and data generated in-house. This database contains sequence information, secondary structure, target genes, expression profiles and a genome browser. In total, there are 8433 miRNAs collected from 121 plant species in PMRD, including model plants and major crops such as Arabidopsis, rice, wheat, soybean, maize, sorghum, barley, etc. For Arabidopsis, rice, poplar, soybean, cotton, medicago and maize, we included the possible target genes for each miRNA with a predicted interaction site in the database. Furthermore, we provided miRNA expression profiles in the PMRD, including our local rice oxidative stress related microarray data (LC Sciences miRPlants_10.1) and the recently published microarray data for poplar, Arabidopsis, tomato, maize and rice. The PMRD database was constructed by open source technology utilizing a user-friendly web interface, and multiple search tools. The PMRD is freely available at http://bioinformatics.cau.edu.cn/PMRD. We expect PMRD to be a useful tool for scientists in the miRNA field in order to study the function of miRNAs and their target genes, especially in model plants and major crops.

Deep sequencing of Brachypodium small RNAs at the global genome level identifies microRNAs involved in cold stress response

Background

MicroRNAs (miRNAs) are endogenous small RNAs having large-scale regulatory effects on plant development and stress responses. Extensive studies of miRNAs have only been performed in a few model plants. Although miRNAs are proved to be involved in plant cold stress responses, little is known for winter-habit monocots. Brachypodium distachyon, with close evolutionary relationship to cool-season cereals, has recently emerged as a novel model plant. There are few reports of Brachypodium miRNAs.

Results

High-throughput sequencing and whole-genome-wide data mining led to the identification of 27 conserved miRNAs, as well as 129 predicted miRNAs in Brachypodium. For multiple-member conserved miRNA families, their sizes in Brachypodium were much smaller than those in rice and Populus. The genome organization of miR395 family in Brachypodium was quite different from that in rice. The expression of 3 conserved miRNAs and 25 predicted miRNAs showed significant changes in response to cold stress. Among these miRNAs, some were cold-induced and some were cold-suppressed, but all the conserved miRNAs were up-regulated under cold stress condition.

Conclusion

Our results suggest that Brachypodium miRNAs are composed of a set of conserved miRNAs and a large proportion of non-conserved miRNAs with low expression levels. Both kinds of miRNAs were involved in cold stress response, but all the conserved miRNAs were up-regulated, implying an important role for cold-induced miRNAs. The different size and genome organization of miRNA families in Brachypodium and rice suggest that the frequency of duplication events or the selection pressure on duplicated miRNAs are different between these two closely related plant species.

Deep sequencing of Brachypodium small RNAs at the global genome level identifies microRNAs involved in cold stress response

BACKGROUND

MicroRNAs (miRNAs) are endogenous small RNAs having large-scale regulatory effects on plant development and stress responses. Extensive studies of miRNAs have only been performed in a few model plants. Although miRNAs are proved to be involved in plant cold stress responses, little is known for winter-habit monocots. Brachypodium distachyon, with close evolutionary relationship to cool-season cereals, has recently emerged as a novel model plant. There are few reports of Brachypodium miRNAs.

RESULTS

High-throughput sequencing and whole-genome-wide data mining led to the identification of 27 conserved miRNAs, as well as 129 predicted miRNAs in Brachypodium. For multiple-member conserved miRNA families, their sizes in Brachypodium were much smaller than those in rice and Populus. The genome organization of miR395 family in Brachypodium was quite different from that in rice. The expression of 3 conserved miRNAs and 25 predicted miRNAs showed significant changes in response to cold stress. Among these miRNAs, some were cold-induced and some were cold-suppressed, but all the conserved miRNAs were up-regulated under cold stress condition.

CONCLUSION

Our results suggest that Brachypodium miRNAs are composed of a set of conserved miRNAs and a large proportion of non-conserved miRNAs with low expression levels. Both kinds of miRNAs were involved in cold stress response, but all the conserved miRNAs were up-regulated, implying an important role for cold-induced miRNAs. The different size and genome organization of miRNA families in Brachypodium and rice suggest that the frequency of duplication events or the selection pressure on duplicated miRNAs are different between these two closely related plant species.

Pyrosequencing of small non-coding RNAs in HIV-1 infected cells: evidence for the processing of a viral-cellular double-stranded RNA hybrid

Small non-coding RNAs of 18-25 nt in length can regulate gene expression through the RNA interference (RNAi) pathway. To characterize small RNAs in HIV-1-infected cells, we performed linker-ligated cloning followed by high-throughput pyrosequencing. Here, we report the composition of small RNAs in HIV-1 productively infected MT4 T-cells. We identified several HIV-1 small RNA clones and a highly abundant small 18-nt RNA that is antisense to the HIV-1 primer-binding site (PBS). This 18-nt RNA apparently originated from the dsRNA hybrid formed by the HIV-1 PBS and the 3' end of the human cellular tRNAlys3. It was found to associate with the Ago2 protein, suggesting its possible function in the cellular RNAi machinery for targeting HIV-1.

Identification of conserved Aquilegia coerulea microRNAs and their targets

Aquilegia is an emerging model organism that is phylogenetically intermediate between the core eudicot and monocot models, Arabidopsis and Oryza. In this study, we have used a comparative genomics approach to identify 45 Aquilegia microRNAs that comprise 20 separate plant microRNA families. We have predicted 84 targets of these newly identified Aquilegia microRNAs including transcription factors and loci involved in metabolism, stress responses, transport, and auxin signaling. microRNA families from 16 plant species and the newly identified microRNAs from Aquilegia were analyzed in a phylogenetic context revealing 40 distantly conserved microRNA families. In addition to these highly conserved plant microRNA families, several families with disjointed phylogenetic distribution were identified. This study provides a phylogenetically important dataset for plant microRNA evolution studies. The current study is the first to identify miRNAs in a lower eudicot in which comprehensive genomic resources are becoming available.

Clusters and superclusters of phased small RNAs in the developing inflorescence of rice

To address the role of small regulatory RNAs in rice development, we generated a large data set of small RNAs from mature leaves and developing roots, shoots, and inflorescences. Using a spatial clustering algorithm, we identified 36,780 genomic groups of small RNAs. Most consisted of 24-nt RNAs that are expressed in all four tissues and enriched in repeat regions of the genome; 1029 clusters were composed primarily of 21-nt small RNAs and, strikingly, 831 of these contained phased RNAs and were preferentially expressed in developing inflorescences. Thirty-eight of the 24-mer clusters were also phased and preferentially expressed in inflorescences. The phased 21-mer clusters derive from nonprotein coding, nonrepeat regions of the genome and are grouped together into superclusters containing 10-46 clusters. The majority of these 21-mer clusters (705/831) are flanked by a degenerate 22-nt motif that is offset by 12 nt from the main phase of the cluster. Small RNAs complementary to these flanking 22-nt motifs define a new miRNA family, which is conserved in maize and expressed in developing reproductive tissues in both plants. These results suggest that the biogenesis of phased inflorescence RNAs resembles that of tasiRNAs and raise the possibility that these novel small RNAs function in early reproductive development in rice and other monocots.

Conserved and novel miRNAs in the legume Phaseolus vulgaris in response to stress

MicroRNAs (miRNAs) are small RNA molecules recognized as important regulators of gene expression. Although plant miRNAs have been extensively studied in model systems, less is known in other plants with limited genome sequence data. We are interested in the identification of miRNAs in Phaseolus vulgaris (common bean) to uncover different plant strategies to cope with adverse conditions and because of its relevance as a crop in developing countries. Here we present the identification of conserved and candidate novel miRNAs in P. vulgaris present in different organs and growth conditions, including drought, abscisic acid treatment, and Rhizobium infection. We also identified cDNA sequences in public databases that represent the corresponding miRNA precursors. In addition, we predicted and validated target mRNAs amongst reported EST and cDNAs for P. vulgaris. We propose that the novel miRNAs present in common bean and other legumes, are involved in regulation of legume-specific processes including adaptation to diverse external cues.

Cloning and characterization of small RNAs from Medicago truncatula reveals four novel legume-specific microRNA families

MicroRNAs (miRNAs) and small-interfering RNAs (siRNAs) have emerged as important regulators of gene expression in higher eukaryotes. Recent studies indicate that genomes in higher plants encode lineage-specific and species-specific miRNAs in addition to the well-conserved miRNAs. Leguminous plants are grown throughout the world for food and forage production. To date the lack of genomic sequence data has prevented systematic examination of small RNAs in leguminous plants. Medicago truncatula, a diploid plant with a near-completely sequenced genome has recently emerged as an important model legume. We sequenced a small RNA library generated from M. truncatula to identify not only conserved miRNAs but also novel small RNAs, if any. Eight novel small RNAs were identified, of which four (miR1507, miR2118, miR2119 and miR2199) are annotated as legume-specific miRNAs because these are conserved in related legumes. Three novel transcripts encoding TIR-NBS-LRR proteins are validated as targets for one of the novel miRNA, miR2118. Small RNA sequence analysis coupled with the small RNA blot analysis, confirmed the expression of around 20 conserved miRNA families in M. truncatula. Fifteen transcripts have been validated as targets for conserved miRNAs. We also characterized Tas3-siRNA biogenesis in M. truncatula and validated three auxin response factor (ARF) transcripts that are targeted by tasiRNAs. These findings indicate that M. truncatula and possibly other related legumes have complex mechanisms of gene regulation involving specific and common small RNAs operating post-transcriptionally.

Methodologies for in vitro cloning of small RNAs and application for plant genome(s)

The "RNA revolution" that started at the end of the 20th century with the discovery of post-transcriptional gene silencing and its mechanism via RNA interference (RNAi) placed tiny 21-24 nucleotide long noncoding RNAs (ncRNAs) in the forefront of biology as one of the most important regulatory elements in a host of physiologic processes. The discovery of new classes of ncRNAs including endogenous small interfering RNAs, microRNAs, and PIWI-interacting RNAs is a hallmark in the understanding of RNA-dependent gene regulation. New generation high-throughput sequencing technologies further accelerated the studies of this "tiny world" and provided their global characterization and validation in many biological systems with sequenced genomes. Nevertheless, for the many "yet-unsequenced" plant genomes, the discovery of small RNA world requires in vitro cloning from purified cellular RNAs. Thus, reproducible methods for in vitro small RNA cloning are of paramount importance and will remain so into the foreseeable future. In this paper, we present a description of existing small RNA cloning methods as well as next-generation sequencing methods that have accelerated this research along with a description of the application of one in vitro cloning method in an initial small RNA survey in the "still unsequenced" allotetraploid cotton genome.

Novel microRNAs uncovered by deep sequencing of small RNA transcriptomes in bread wheat (Triticum aestivum L.) and Brachypodium distachyon (L.) Beauv

The small RNA transcriptomes of bread wheat and its emerging model Brachypodium distachyon were obtained by using deep sequencing technology. Small RNA compositions were analyzed in these two species. In addition to 70 conserved microRNAs (miRNAs) from 25 families, 23 novel wheat miRNAs were identified. For Brachypodium, 12 putative miRNAs were predicted from a limited number of expressed sequence tags, of which one was a potential novel miRNA. Also, 94 conserved miRNAs from 28 families were identified in this species. Expression validation was performed for several novel wheat miRNAs. RNA ligase-mediated 5' rapid amplification of complementary DNA ends experiments demonstrated their capability to cleave predicted target genes including three disease-resistant gene analogs. Differential expression of miRNAs was observed between Brachypodium vegetative and reproductive tissues, suggesting their different roles at the two growth stages. Our work significantly increases the novel miRNA numbers in wheat and provides the first set of small RNAs in B. distachyon.

MicroRNA expression profiles in conventional and micropropagated strawberry (Fragaria x ananassa Duch.) plants

MicroRNAs (miRNAs) are a class of small non-coding RNAs which play a critical role in plant growth and development. To detect strawberry miRNAs and discover the expression difference between conventional and micropropagated strawberry plants, we carried out the detection and quantification of strawberry miRNAs by microarray. The main findings were that 74 miRNAs were checked in strawberry plants and four miRNA genes displayed clear expression difference between conventional and micropropagated strawberry plants, including two up-regulated genes (miR535 and miR390) and two down-regulated genes (miR169a and miR169d). The ratios of conventionally propagated strawberry plant/micropropagated strawberry plant for miR535, miR390, miR169a and miR169d were 2.6884, 2.2673, 0.2496 and 0.3814, respectively. Quantitative reverse transcription polymerase chain reaction applied to the two up-regulated genes (miR535 and miR390) validated the microarray result. This is the first report on differential expression of miRNAs in conventional and micropropagated plants.

Parallel DNA pyrosequencing unveils new zebrafish microRNAs

Background

MicroRNAs (miRNAs) are a new class of small RNAs of approximately 22 nucleotides in length that control eukaryotic gene expression by fine tuning mRNA translation. They regulate a wide variety of biological processes, namely developmental timing, cell differentiation, cell proliferation, immune response and infection. For this reason, their identification is essential to understand eukaryotic biology. Their small size, low abundance and high instability complicated early identification, however cloning/Sanger sequencing and new generation genome sequencing approaches overcame most technical hurdles and are being used for rapid miRNA identification in many eukaryotes.

Results

We have applied 454 DNA pyrosequencing technology to miRNA discovery in zebrafish (Danio rerio). For this, a series of cDNA libraries were prepared from miRNAs isolated at different embryonic time points and from fully developed organs. Each cDNA library was tagged with specific sequences and was sequenced using the Roche FLX genome sequencer. This approach retrieved 90% of the 192 miRNAs previously identified by cloning/Sanger sequencing and bioinformatics. Twenty five novel miRNAs were predicted, 107 miRNA star sequences and also 41 candidate miRNA targets were identified. A miRNA expression profile built on the basis of pyrosequencing read numbers showed high expression of most miRNAs throughout zebrafish development and identified tissue specific miRNAs.

Conclusion

This study increases the number of zebrafish miRNAs from 192 to 217 and demonstrates that a single DNA mini-chip pyrosequencing run is effective in miRNA identification in zebrafish. This methodology also produced sufficient information to elucidate miRNA expression patterns during development and in differentiated organs. Moreover, some zebrafish miRNA star sequences were more abundant than their corresponding miRNAs, suggesting a functional role for the former in gene expression control in this vertebrate model organism.

Next is now: new technologies for sequencing of genomes, transcriptomes, and beyond

The sudden availability of DNA sequencing technologies that rapidly produce vast amounts of sequence information has triggered a paradigm shift in genomics, enabling massively parallel surveying of complex nucleic acid populations. The diversity of applications to which these technologies have already been applied demonstrates the immense range of cellular processes and properties that can now be studied at the single-base resolution. These include genome resequencing and polymorphism discovery, mutation mapping, DNA methylation, histone modifications, transcriptome sequencing, gene discovery, alternative splicing identification, small RNA profiling, DNA-protein, and possibly even protein-protein interactions. Thus, these deep sequencing technologies offer plant biologists unprecedented opportunities to increase the understanding of the functions and dynamics of plant cells and populations.

Current tools for the identification of miRNA genes and their targets

The discovery of microRNAs (miRNAs), almost 10 years ago, changed dramatically our perspective on eukaryotic gene expression regulation. However, the broad and important functions of these regulators are only now becoming apparent. The expansion of our catalogue of miRNA genes and the identification of the genes they regulate owe much to the development of sophisticated computational tools that have helped either to focus or interpret experimental assays. In this article, we review the methods for miRNA gene finding and target identification that have been proposed in the last few years. We identify some problems that current approaches have not yet been able to overcome and we offer some perspectives on the next generation of computational methods.

Cloning, characterization and expression analysis of porcine microRNAs

Background

MicroRNAs (miRNAs) are small ~22-nt regulatory RNAs that can silence target genes, by blocking their protein production or degrading the mRNAs. Pig is an important animal in the agriculture industry because of its utility in the meat production. Besides, pig has tremendous biomedical importance as a model organism because of its closer proximity to humans than the mouse model. Several hundreds of miRNAs have been identified from mammals, humans, mice and rats, but little is known about the miRNA component in the pig genome. Here, we adopted an experimental approach to identify conserved and unique miRNAs and characterize their expression patterns in diverse tissues of pig.

Results

By sequencing a small RNA library generated using pooled RNA from the pig heart, liver and thymus; we identified a total of 120 conserved miRNA homologs in pig. Expression analysis of conserved miRNAs in 14 different tissue types revealed heart-specific expression of miR-499 and miR-208 and liver-specific expression of miR-122. Additionally, miR-1 and miR-133 in the heart, miR-181a and miR-142-3p in the thymus, miR-194 in the liver, and miR-143 in the stomach showed the highest levels of expression. miR-22, miR-26b, miR-29c and miR-30c showed ubiquitous expression in diverse tissues. The expression patterns of pig-specific miRNAs also varied among the tissues examined.

Conclusion

Identification of 120 miRNAs and determination of the spatial expression patterns of a sub-set of these in the pig is a valuable resource for molecular biologists, breeders, and biomedical investigators interested in post-transcriptional gene regulation in pig and in related mammals, including humans.

Short-read sequencing technologies for transcriptional analyses

The technological advances in DNA sequencing over the past five years have changed our approaches to gene expression analysis, fundamentally altering the basic methods used and in most cases driving a shift from hybridization-based approaches to sequencing-based approaches. Quantitative, tag-based studies of gene expression were one of the earliest applications of these next-generation technologies, but the tremendous depth of sequencing facilitates de novo transcript discovery, which replaces traditional expressed sequence tag (EST) sequencing. In addition, these technologies have created new opportunities for understanding the generation, stability, and decay of RNA and the impacts of chromatin differences on gene expression. As we review the impact of these methods on plant biology, we also mention published studies from animal systems when the methods are broadly applicable. We can anticipate that the published work over the past few years is a harbinger of much broader studies that are yet to be published and are sure to further advance our understanding of plant genomes in a field changing at a dizzying pace.

Differential expression of miRNAs in response to salt stress in maize roots

BACKGROUND AND AIMS

Corn (Zea mays) responds to salt stress via changes in gene expression, metabolism and physiology. This adaptation is achieved through the regulation of gene expression at the transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) have been found to act as key regulating factors of post-transcriptional gene expression. However, little is known about the role of miRNAs in plants' responses to abiotic stresses.

METHODS

A custom microparaflo microfluidic array containing release version 10.1 plant miRNA probes (http://microrna.sanger.ac.uk/) was used to discover salt stress-responsive miRNAs using the differences in miRNA expression between the salt-tolerant maize inbred line 'NC286' and the salt-sensitive maize line 'Huangzao4'. Key Results miRNA microarray hybridization revealed that a total of 98 miRNAs, from 27 plant miRNA families, had significantly altered expression after salt treatment. These miRNAs displayed different activities in the salt response, and miRNAs belonging to the same miRNA family showed the same behaviour. Interestingly, 18 miRNAs were found which were only expressed in the salt-tolerant maize line, and 25 miRNAs that showed a delayed regulation pattern in the salt-sensitive line. A gene model was proposed that showed how miRNAs could regulate the abiotic stress-associated process and the gene networks coping with the stress.

CONCLUSIONS

Salt-responsive miRNAs are involved in the regulation of metabolic, morphological and physiological adaptations of maize seedlings at the post-transcriptional level. The miRNA genotype-specific expression model might explain the distinct salt sensitivities between maize lines.

Characterization and expression profiles of miRNAs in rice seeds

Small RNAs (sRNAs) are common and effective modulators of gene expression in eukaryotic organisms. To characterize the sRNAs expressed during rice seed development, massively parallel signature sequencing (MPSS) was performed, resulting in the obtainment of 797 399 22-nt sequence signatures, of which 111 161 are distinct ones. Analysis on the distributions of sRNAs on chromosomes showed that most sRNAs originate from interspersed repeats that mainly consist of transposable elements, suggesting the major function of sRNAs in rice seeds is transposon silencing. Through integrative analysis, 26 novel miRNAs and 12 miRNA candidates were identified. Further analysis on the expression profiles of the known and novel miRNAs through hybridizing the generated chips revealed that most miRNAs were expressed preferentially in one or two rice tissues. Detailed comparison of the expression patterns of miRNAs and corresponding target genes revealed the negative correlation between them, while few of them are positively correlated. In addition, differential accumulations of miRNAs and corresponding miRNA*s suggest the functions of miRNA*s other than being passenger strands of mature miRNAs, and in regulating the miRNA functions.

Identification of precursor transcripts for 6 novel miRNAs expands the diversity on the genomic organisation and expression of miRNA genes in rice

BACKGROUND

The plant miRNAs represent an important class of endogenous small RNAs that guide cleavage of an mRNA target or repress its translation to control development and adaptation to stresses. MiRNAs are nuclear-encoded genes transcribed by RNA polymerase II, producing a primary precursor that is subsequently processed by DCL1 an RNase III Dicer-like protein. In rice hundreds of miRNAs have been described or predicted, but little is known on their genes and precursors which are important criteria to distinguish them from siRNAs. Here we develop a combination of experimental approaches to detect novel miRNAs in rice, identify their precursor transcripts and genes and predict or validate their mRNA targets.

RESULTS

We produced four cDNA libraries from small RNA fractions extracted from distinct rice tissues. By in silico analysis we selected 6 potential novel miRNAs, and confirmed that their expression requires OsDCL1. We predicted their targets and used 5'RACE to validate cleavage for three of them, targeting a PPR, an SPX domain protein and a GT-like transcription factor respectively. In addition, we identified precursor transcripts for the 6 miRNAs expressed in rice, showing that these precursors can be efficiently processed using a transient expression assay in transfected Nicotiana benthamiana leaves. Most interestingly, we describe two precursors producing tandem miRNAs, but in distinct arrays. We focus on one of them encoding osa-miR159a.2, a novel miRNA produced from the same stem-loop structure encoding the conserved osa-miR159a.1. We show that this dual osa-miR159a.2-osa-miR159a.1 structure is conserved in distant rice species and maize. Finally we show that the predicted mRNA target of osa-miR159a.2 encoding a GT-like transcription factor is cleaved in vivo at the expected site.

CONCLUSION

The combination of approaches developed here identified six novel miRNAs expressed in rice which can be clearly distinguished from siRNAs. Importantly, we show that two miRNAs can be produced from a single precursor, either from tandem stem-loops or tandemly arrayed in a single stem-loop. This suggests that processing of these precursors could be an important regulatory step to produce one or more functional miRNAs in plants and perhaps coordinate cleavage of distinct targets in the same plant tissue.

Genome-wide identification and analysis of small RNAs originated from natural antisense transcripts in Oryza sativa

Natural antisense transcripts (NATs) have been shown to play important roles in post-transcriptional regulation through the RNA interference pathway. We have combined pyrophosphate-based high-throughput sequencing and computational analysis to identify and analyze, in genome scale, cis-NAT and trans-NAT small RNAs that are derived under normal conditions and in response to drought and salt stresses in the staple plant Oryza sativa. Computationally, we identified 344 cis-NATs and 7,142 trans-NATs that are formed by protein-coding genes. From the deep sequencing data, we found 108 cis-NATs and 7,141 trans-NATs that gave rise to small RNAs from their overlapping regions. Consistent with early findings, the majority of these 108 cis-NATs seem to be associated with specific conditions or developmental stages. Our analyses also revealed several interesting results. The overlapping regions of the cis-NATs and trans-NATs appear to be more enriched with small RNA loci than non-overlapping regions. The small RNAs generated from cis-NATs and trans-NATs have a length bias of 21 nt, even though their lengths spread over a large range. Furthermore, >40% of the small RNAs from cis-NATs and trans-NATs carry an A as their 5'-terminal nucleotides. A substantial portion of the transcripts are involved in both cis-NATs and trans-NATs, and many trans-NATs can form many-to-many relationships, indicating that NATs may form complex regulatory networks in O. sativa. This study is the first genome-wide investigation of NAT-derived small RNAs in O. sativa. It reveals the importance of NATs in biogenesis of small RNAs and broadens our understanding of the roles of NAT-derived small RNAs in gene regulation, particularly in response to environmental stimuli.

Functional genetic variation of human miRNAs and phenotypic consequences

A large number of human protein-coding genes are finely regulated by one or more microRNAs. Members of this small noncoding RNA family have emerged as important post-transcriptional regulators of gene expression and are involved in a number of disease phenotypes. Variability in the human genome is extensive and includes the common and rare single nucleotide polymorphisms (SNPs) and copy number variations (CNVs). The functional significance of the genome's variability is under intense investigation. In this article we review the emerging literature on how human genomic variation influences the outcome of microRNA targeting and the associated phenotypic effects. Illustrative examples are discussed that demonstrate the biological importance of functional polymorphisms affecting miRNA-mediated gene regulation.

A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains

Endogenous small RNAs, including microRNAs (miRNAs) and short-interfering RNAs (siRNAs), function as post-transcriptional or transcriptional regulators in plants. miRNA function is essential for normal plant development and therefore is likely to be important in the growth of the rice grain. To investigate the roles of miRNAs in rice grain development, we carried out deep sequencing of the small RNA populations of rice grains at two developmental stages. In a data set of approximately 5.5 million sequences, we found representatives of all 20 conserved plant miRNA families. We used an approach based on the presence of miRNA and miRNA* sequences to identify 39 novel, nonconserved rice miRNA families expressed in grains. Cleavage of predicted target mRNAs was confirmed for a number of the new miRNAs. We identified a putative mirtron, indicating that plants may also use spliced introns as a source of miRNAs. We also identified a miRNA-like long hairpin that generates phased 21 nt small RNAs, strongly expressed in developing grains, and show that these small RNAs act in trans to cleave target mRNAs. Comparison of the population of miRNAs and miRNA-like siRNAs in grains to those in other parts of the rice plant reveals that many are expressed in an organ-specific manner.

In silico identification of conserved microRNAs in large number of diverse plant species

BACKGROUND

MicroRNAs (miRNAs) are recently discovered small non-coding RNAs that play pivotal roles in gene expression, specifically at the post-transcriptional level in plants and animals. Identification of miRNAs in large number of diverse plant species is important to understand the evolution of miRNAs and miRNA-targeted gene regulations. Now-a-days, publicly available databases play a central role in the in-silico biology. Because, at least ~21 miRNA families are conserved in higher plants, a homology based search using these databases can help identify orthologs or paralogs in plants.

RESULTS

We searched all publicly available nucleotide databases of genome survey sequences (GSS), high-throughput genomics sequences (HTGS), expressed sequenced tags (ESTs) and nonredundant (NR) nucleotides and identified 682 miRNAs in 155 diverse plant species. We found more than 15 conserved miRNA families in 11 plant species, 10 to14 families in 10 plant species and 5 to 9 families in 29 plant species. Nineteen conserved miRNA families were identified in important model legumes such as Medicago, Lotus and soybean. Five miRNA families - miR319, miR156/157, miR169, miR165/166 and miR394 - were found in 51, 45, 41, 40 and 40 diverse plant species, respectively. miR403 homologs were found in 16 dicots, whereas miR437 and miR444 homologs, as well as the miR396d/e variant of the miR396 family, were found only in monocots, thus providing large-scale authenticity for the dicot- and monocot-specific miRNAs. Furthermore, we provide computational and/or experimental evidence for the conservation of 6 newly found Arabidopsis miRNA homologs (miR158, miR391, miR824, miR825, miR827 and miR840) and 2 small RNAs (small-85 and small-87) in Brassica spp.

CONCLUSION

Using all publicly available nucleotide databases, 682 miRNAs were identified in 155 diverse plant species. By combining the expression analysis with the computational approach, we found that 6 miRNAs and 2 small RNAs that have been identified only in Arabidopsis thus far, are also conserved in Brassica spp. These findings will be useful for tracing the evolution of small RNAs by examining their expression in common ancestors of the Arabidopsis-Brassica lineage.