Deciphering the diversity of small RNAs in plants: the long and short of it

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

fIdentification of B. napus small RNAs responsive to infection by a necrotrophic pathogen

BACKGROUND

Small RNAs are short non-coding RNAs that are key gene regulators controlling various biological processes in eukaryotes. Plants may regulate discrete sets of sRNAs in response to pathogen attack. Sclerotinia sclerotiorum is an economically important pathogen affecting hundreds of plant species, including the economically important oilseed B. napus. However, there are limited studies on how regulation of sRNAs occurs in the S. sclerotiorum and B. napus pathosystem.

RESULTS

We identified different classes of sRNAs from B. napus using high throughput sequencing of replicated mock and infected samples at 24 h post-inoculation (HPI). Overall, 3999 sRNA loci were highly expressed, of which 730 were significantly upregulated during infection. These 730 up-regulated sRNAs targeted 64 genes, including disease resistance proteins and transcriptional regulators. A total of 73 conserved miRNA families were identified in our dataset. Degradome sequencing identified 2124 cleaved mRNA products from these miRNAs from combined mock and infected samples. Among these, 50 genes were specific to infection. Altogether, 20 conserved miRNAs were differentially expressed and 8 transcripts were cleaved by the differentially expressed miRNAs miR159, miR5139, and miR390, suggesting they may have a role in the S. sclerotiorum response. A miR1885-triggered disease resistance gene-derived secondary sRNA locus was also identified and verified with degradome sequencing. We also found further evidence for silencing of a plant immunity related ethylene response factor gene by a novel sRNA using 5'-RACE and RT-qPCR.

CONCLUSIONS

The findings in this study expand the framework for understanding the molecular mechanisms of the S. sclerotiorum and B. napus pathosystem at the sRNA level.

Transient expression of intron-containing transgenes generates non-spliced aberrant pre-mRNAs that are processed into siRNAs

MAIN CONCLUSION

In this study, we show that aberrant pre-mRNAs from non-spliced and non-polyadenylated intron-containing transgenes are channelled to the RNA silencing pathway. In plants, improperly processed transcripts are called aberrant RNAs (ab-RNAs) and are eliminated by either RNA silencing or RNA decay mechanisms. Ab-RNAs transcribed from intronless genes are copied by RNA-directed RNA polymerases (RDRs) into double-stranded RNAs which are subsequently cleaved by DICER-LIKE endonucleases into small RNAs (sRNAs). In contrast, ab-RNAs from intron-containing genes are suggested to be channelled post-splicing to exonucleolytic degradation. Yet, it is not clear how non-spliced aberrant pre-mRNAs are eliminated. We reasoned that transient expression of agroinfiltrated intron-containing transgenes in Nicotiana benthamiana would allow us to study the steady-state levels of non-spliced pre-mRNAs. SRNA deep sequencing of the agroinfiltrated transgenes revealed the presence of sRNAs mapping to the entire non-spliced pre-mRNA suggesting that RDRs (most likely RDR6) processed aberrant non-spliced pre-mRNAs. Primary and secondary sRNAs with lengths of 18-25 nucleotides (nt) were detected, with the most prominent sRNA size class of 22 nt. SRNAs also mapped to the terminator sequence, indicating that RDR substrates also comprised read-through transcripts devoid of polyadenylation tail. Importantly, the occurring sRNAs efficiently targeted cognate mRNA for degradation but failed to cleave the non-spliced pre-mRNA, corroborating the notion that sRNAs are not triggering RNA cleavage in the nucleus.

Transposon-associated epigenetic silencing during Pleurotus ostreatus life cycle

Transposable elements constitute an important fraction of eukaryotic genomes. Given their mutagenic potential, host-genomes have evolved epigenetic defense mechanisms to limit their expansion. In fungi, epigenetic modifications have been widely studied in ascomycetes, although we lack a global picture of the epigenetic landscape in basidiomycetes. In this study, we analysed the genome-wide epigenetic and transcriptional patterns of the white-rot basidiomycete Pleurotus ostreatus throughout its life cycle. Our results performed by using high-throughput sequencing analyses revealed that strain-specific DNA methylation profiles are primarily involved in the repression of transposon activity and suggest that 21 nt small RNAs play a key role in transposon silencing. Furthermore, we provide evidence that transposon-associated DNA methylation, but not sRNA production, is directly involved in the silencing of genes surrounded by transposons. Remarkably, we found that nucleus-specific methylation levels varied in dikaryotic strains sharing identical genetic complement but different subculture conditions. Finally, we identified key genes activated in the fruiting process through the comparative analysis of transcriptomes. This study provides an integrated picture of epigenetic defense mechanisms leading to the transcriptional silencing of transposons and surrounding genes in basidiomycetes. Moreover, our findings suggest that transcriptional but not methylation reprogramming triggers fruitbody development in P. ostreatus.

MicroRNA Regulation of Abiotic Stress Response in 7B-1 Male-Sterile Tomato Mutant

The 7B-1 tomato (Solanum lycopersicum L. 'Rutgers') is a male-sterile mutant with enhanced tolerance to abiotic stress in a blue-light (BL) specific manner compared with its wild-type (WT). This makes the 7B-1 a potential candidate for hybrid seed breeding and stress engineering. To identify small RNAs (sRNAs) linked to stress tolerance of 7B-1, two sRNA libraries from BL-grown 7B-1 and WT seedlings treated simultaneously with abscisic acid (ABA) and mannitol were sequenced, and sRNA profiles were compared. Twenty nine families of known microRNAs (miRNAs) and 27 putative novel miRNAs were identified from the two libraries. MiR5300, miR5301, miR2916, and a novel miRNA denoted miR#C were upregulated, while miR159, miR166, miR472, miR482, and two novel miRNAs, miR#A and miR#D, were downregulated in stress-treated 7B-1 seedlings. MiRNA targets with potential roles in stress regulation were validated by rapid amplification of 5' complementary DNA ends (5'-RACE) analysis. Expression of miR159, miR166, miR472, miR482, miR#A, and miR#D together with their targets were further investigated in response to ABA, mannitol, NaCl, and cold treatments and a strong negative correlation was observed between the levels of these miRNAs and expression of their targets. Only miR159 and miR166 responded to cold treatment. MiR#A and its target were regulated by ABA and mannitol as early as 0.5 h after the treatments, while other miRNAs and targets were regulated only after 2 h. This suggests a role in early response to stress for miR#A. Our data suggests that miR159, miR166, miR472, miR482, miR#A, and miR#D are likely to facilitate the BL-specific enhanced tolerance of 7B-1 to abiotic stress.

Identification of miRNAs with potential roles in regulation of anther development and male-sterility in 7B-1 male-sterile tomato mutant

BACKGROUND

The 7B-1 tomato line (Solanum lycopersicum cv. Rutgers) is a photoperiod-sensitive male-sterile mutant, with potential application in hybrid seed production. Small RNAs (sRNAs) in tomato have been mainly characterized in fruit development and ripening, but none have been studied with respect to flower development and regulation of male-sterility. Using sRNA sequencing, we identified miRNAs that are potentially involved in anther development and regulation of male-sterility in 7B-1 mutant.

RESULTS

Two sRNA libraries from 7B-1 and wild type (WT) anthers were sequenced and thirty two families of known miRNAs and 23 new miRNAs were identified in both libraries. MiR390, miR166, miR159 were up-regulated and miR530, miR167, miR164, miR396, miR168, miR393, miR8006 and two new miRNAs, miR#W and miR#M were down-regulated in 7B-1 anthers. Ta-siRNAs were not differentially expressed and likely not associated with 7B-1 male-sterility. miRNA targets with potential roles in anther development were validated using 5'-RACE. QPCR analysis showed differential expression of miRNA/target pairs of interest in anthers and stem of 7B-1, suggesting that they may regulate different biological processes in these tissues. Expression level of most miRNA/target pairs showed negative correlation, except for few. In situ hybridization showed predominant expression of miR159, GAMYBL1, PMEI and cystatin in tapetum, tetrads and microspores.

CONCLUSION

Overall, we identified miRNAs with potential roles in anther development and regulation of male-sterility in 7B-1. A number of new miRNAs were also identified from tomato for the first time. Our data could be used as a benchmark for future studies of the molecular mechanisms of male-sterility in other crops.

Drawing siRNAs of viral origin out from plant siRNAs libraries

Viruses are obligate intracellular entities that infect all forms of life. In plants, invading viral nucleic acids trigger RNA silencing machinery and it results in the accumulation of viral short interfering RNAs (v-siRNAs). The study of v-siRNAs population in biological samples has become a major part of many research projects aiming to identify viruses infecting them, including unknown viruses, even at extremely low titer. Currently, siRNA populations are investigated by high-throughput sequencing approaches, which generate very large data sets. The major difficulty in these studies is to properly analyze such huge amount of data. In this regard, easy-to-use bioinformatics tools to groom and decipher siRNA libraries and to draw out v-siRNAs are needed. Here we describe a workflow, which permit users with little experience in bioinformatics to draw out v-siRNAs from raw data sequences obtained by Illumina technology. Such pipeline has been released in the context of Galaxy, an open source Web-based platform for bioinformatics analyses.

Novel miRNA and phasiRNA biogenesis networks in soybean roots from two sister lines that are resistant and susceptible to SCN race 4

The soybean cyst nematode (SCN), Heterodera glycines, is the most devastating pathogen of soybean worldwide. SiRNAs (small interfere RNAs) have been proven to induce the silencing of cyst nematode genes. However, whether small RNAs from soybean root have evolved a similar mechanism against SCN is unknown. Two genetically related soybean sister lines (ZP03-5373 and ZP03-5413), which are resistant and susceptible, respectively, to SCN race 4 infection were selected for small RNA deep sequencing to identify small RNAs targeted to SCN. We identified 71 less-conserved miRNAs-miRNAs* counterparts belonging to 32 families derived from 91 loci, and 88 novel soybean-specific miRNAs with distinct expression patterns. The identified miRNAs targeted 42 genes representing a wide range of enzymatic and regulatory activities. Roots of soybean conserved one TAS (Trans-acting siRNA) gene family with a similar but unique trans-acting small interfering RNA (tasiRNA) biogenesis profile. In addition, we found that six miRNAs (gma-miR393, 1507, 1510, 1515, 171, 2118) guide targets to produce secondary phasiRNAs (phased, secondary, small interfering RNAs) in soybean root. Multiple targets of these phasiRNAs were predicted and detected. Importantly, we also found that the expression of 34 miRNAs differed significantly between the two lines. Seven ZP03-5373-specific miRNAs were differentially expressed after SCN infection. Forty-four transcripts from SCN were predicted to be potential targets of ZP03-5373-specific differential miRNAs. These findings suggest that miRNAs play an important role in the soybean response to SCN.

The diversity of small non-coding RNAs in the diatom Phaeodactylum tricornutum

Background

Marine diatoms constitute a major component of eukaryotic phytoplankton and stand at the crossroads of several evolutionary lineages. These microalgae possess peculiar genomic features and novel combinations of genes acquired from bacterial, animal and plant ancestors. Furthermore, they display both DNA methylation and gene silencing activities. Yet, the biogenesis and regulatory function of small RNAs (sRNAs) remain ill defined in diatoms.

Results

Here we report the first comprehensive characterization of the sRNA landscape and its correlation with genomic and epigenomic information in Phaeodactylum tricornutum. The majority of sRNAs is 25 to 30 nt-long and maps to repetitive and silenced Transposable Elements marked by DNA methylation. A subset of this population also targets DNA methylated protein-coding genes, suggesting that gene body methylation might be sRNA-driven in diatoms. Remarkably, 25-30 nt sRNAs display a well-defined and unprecedented 180 nt-long periodic distribution at several highly methylated regions that awaits characterization. While canonical miRNAs are not detectable, other 21-25 nt sRNAs of unknown origin are highly expressed. Besides, non-coding RNAs with well-described function, namely tRNAs and U2 snRNA, constitute a major source of 21-25 nt sRNAs and likely play important roles under stressful environmental conditions.

Conclusions

P. tricornutum has evolved diversified sRNA pathways, likely implicated in the regulation of largely still uncharacterized genetic and epigenetic processes. These results uncover an unexpected complexity of diatom sRNA population and previously unappreciated features, providing new insights into the diversification of sRNA-based processes in eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-698) contains supplementary material, which is available to authorized users.

Deciphering small noncoding RNAs during the transition from dormant embryo to germinated embryo in Larches (Larix leptolepis)

Small RNAs (sRNAs), as a key component of molecular biology, play essential roles in plant development, hormone signaling, and stress response. However, little is known about the relationships among sRNAs, hormone signaling, and dormancy regulation in gymnosperm embryos. To investigate the roles of sRNAs in embryo dormancy maintenance and release in Larix leptolepis, we deciphered the endogenous “sRNAome” in dormant and germinated embryos. High-throughput sequencing of sRNA libraries showed that dormant embryos exhibited a length bias toward 24-nt while germinated embryos showed a bias toward 21-nt lengths. This might be associated with distinct levels of RNA-dependent RNA polymerase2 (RDR2) and/or RDR6, which is regulated by hormones. Proportions of miRNAs to nonredundant and redundant sRNAs were higher in germinated embryos than in dormant embryos, while the ratio of unknown sRNAs was higher in dormant embryos than in germinated embryos. We identified a total of 160 conserved miRNAs from 38 families, 3 novel miRNAs, and 16 plausible miRNA candidates, of which many were upregulated in germinated embryos relative to dormant embryos. These findings indicate that larches and possibly other gymnosperms have complex mechanisms of gene regulation involving miRNAs and other sRNAs operating transcriptionally and posttranscriptionally during embryo dormancy and germination. We propose that abscisic acid modulates embryo dormancy and germination at least in part through regulation of the expression level of sRNA-biogenesis genes, thus changing the sRNA components.

CoLIde: a bioinformatics tool for CO-expression-based small RNA Loci Identification using high-throughput sequencing data

Small RNAs (sRNAs) are 20-25 nt non-coding RNAs that act as guides for the highly sequence-specific regulatory mechanism known as RNA silencing. Due to the recent increase in sequencing depth, a highly complex and diverse population of sRNAs in both plants and animals has been revealed. However, the exponential increase in sequencing data has also made the identification of individual sRNA transcripts corresponding to biological units (sRNA loci) more challenging when based exclusively on the genomic location of the constituent sRNAs, hindering existing approaches to identify sRNA loci. To infer the location of significant biological units, we propose an approach for sRNA loci detection called CoLIde (Co-expression based sRNA Loci Identification) that combines genomic location with the analysis of other information such as variation in expression levels (expression pattern) and size class distribution. For CoLIde, we define a locus as a union of regions sharing the same pattern and located in close proximity on the genome. Biological relevance, detected through the analysis of size class distribution, is also calculated for each locus. CoLIde can be applied on ordered (e.g., time-dependent) or un-ordered (e.g., organ, mutant) series of samples both with or without biological/technical replicates. The method reliably identifies known types of loci and shows improved performance on sequencing data from both plants (e.g., A. thaliana, S. lycopersicum) and animals (e.g., D. melanogaster) when compared with existing locus detection techniques. CoLIde is available for use within the UEA Small RNA Workbench which can be downloaded from: http://srna-workbench.cmp.uea.ac.uk.

An endogene-resembling transgene delays the onset of silencing and limits siRNA accumulation

In plants, transgenes are generally more sensitive against RNA silencing than endogenes are. In this study, we generated a transgene that structurally mimicks an endogene. It is composed of endogenous promoter, 5'-UTR, introns, 3'-UTR and terminator elements. Our data revealed that, in contrast to a conventional transgene, an endogene-resembling transgene was more stably expressed and poorly processed into small RNAs. In addition, although both constructs triggered methylation of homologous DNA sequences at similar levels, the endogene-resembling transgene exhibited significantly delayed onset of local and systemic silencing.

Identification of miRNAs and their targets from Brassica napus by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are endogenous regulators of a broad range of physiological processes and act by either degrading mRNA or blocking its translation. Oilseed rape (Brassica napus) is one of the most important crops in China, Europe and other Asian countries with publicly available expressed sequence tags (ESTs) and genomic survey sequence (GSS) databases, but little is known about its miRNAs and their targets. To date, only 46 miRNAs have been identified in B. napus.

RESULTS

Forty-one conserved and 62 brassica-specific candidate B. napus miRNAs, including 20 miRNA* sequences, were identified using Solexa sequencing technology. Furthermore, 33 non-redundant mRNA targets of conserved brassica miRNAs and 19 new non-redundant mRNA targets of novel brassica-specific miRNAs were identified by genome-scale sequencing of mRNA degradome.

CONCLUSIONS

This study describes large scale cloning and characterization of B. napus miRNAs and their potential targets, providing the foundation for further characterization of miRNA function in the regulation of diverse physiological processes in B. napus.

Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing

The male sterility of thermosensitive genic male sterile (TGMS) lines of wheat (Triticum aestivum) is strictly controlled by temperature. The early phase of anther development is especially susceptible to cold stress. MicroRNAs (miRNAs) play an important role in plant development and in responses to environmental stress. In this study, deep sequencing of small RNA (smRNA) libraries obtained from spike tissues of the TGMS line under cold and control conditions identified a total of 78 unique miRNA sequences from 30 families and trans-acting small interfering RNAs (tasiRNAs) derived from two TAS3 genes. To identify smRNA targets in the wheat TGMS line, we applied the degradome sequencing method, which globally and directly identifies the remnants of smRNA-directed target cleavage. We identified 26 targets of 16 miRNA families and three targets of tasiRNAs. Comparing smRNA sequencing data sets and TaqMan quantitative polymerase chain reaction results, we identified six miRNAs and one tasiRNA (tasiRNA-ARF [for Auxin-Responsive Factor]) as cold stress-responsive smRNAs in spike tissues of the TGMS line. We also determined the expression profiles of target genes that encode transcription factors in response to cold stress. Interestingly, the expression of cold stress-responsive smRNAs integrated in the auxin-signaling pathway and their target genes was largely noncorrelated. We investigated the tissue-specific expression of smRNAs using a tissue microarray approach. Our data indicated that miR167 and tasiRNA-ARF play roles in regulating the auxin-signaling pathway and possibly in the developmental response to cold stress. These data provide evidence that smRNA regulatory pathways are linked with male sterility in the TGMS line during cold stress.

Sculpting the maturation, softening and ethylene pathway: the influences of microRNAs on tomato fruits

Background

MicroRNAs (miRNAs), a ubiquitous class of short RNAs, play vital roles in physiological and biochemical processes in plants by mediating gene silencing at post-transcriptional (PTGS) level. Tomato is a model system to study molecular basis of fleshy fruit ripening and senescence, ethylene biosynthesis and signal transduction owing to its genetic and molecular tractability. To study the functions of miRNAs in tomato fruit ripening and senescence, and their possible roles in ethylene response, the next generation sequencing method was employed to identify miRNAs in tomato fruit. Bioinformatics and molecular biology approaches were combined to profile the miRNAs expression patterns at three different fruit ripening stages and by exogenous ethylene treatment.

Results

In addition to 7 novel miRNA families, 103 conserved miRNAs belonging to 24 families and 10 non-conserved miRNAs matching 9 families were identified in our libraries. The targets of many these miRNAs were predicted to be transcriptional factors. Other targets are known to play roles in the regulation of metabolic processes. Interestingly, some targets were predicted to be involved in fruit ripening and softening, such as Pectate Lyase, beta-galactosidase, while a few others were predicted to be involved in ethylene biosynthesis and signaling pathway, such as ACS, EIN2 and CTR1. The expression patterns of a number of such miRNAs at three ripening stages were confirmed by stem-loop RT-PCR, which showed a strong negative correlation with that of their targets. The regulation of exogenous ethylene on miRNAs expression profiles were analyzed simultaneously, and 3 down-regulated, 5 up-regulated miRNAs were found in this study.

Conclusions

A combination of high throughput sequencing and molecular biology approaches was used to explore the involvement of miRNAs during fruit ripening. Several miRNAs showed differential expression profiles during fruit ripening, and a number of miRNAs were influenced by ethylene treatment. The results suggest the importance of miRNAs in fruit ripening and ethylene response.

Deep sequencing of small RNAs in plants: applied bioinformatics

Small RNAs, including microRNA and short-interfering RNAs, play important roles in plants. In recent years, developments in sequencing technology have enabled the large-scale discovery of sRNAs in various cells, tissues and developmental stages and in response to various stresses. This review describes the bioinformatics challenges to analysing these large datasets of short-RNA sequences and some of the solutions to those challenges.

Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns

Plants feature a particularly diverse population of short (s)RNAs, the central component of all RNA silencing pathways. Next generation sequencing techniques enable deeper insights into this complex and highly conserved mechanism and allow identification and quantification of sRNAs. We employed deep sequencing to monitor the sRNAome of developing tomato fruits covering the period between closed flowers and ripened fruits by profiling sRNAs at 10 time-points. It is known that microRNAs (miRNAs) play an important role in development but very little information is available about the majority of sRNAs that are not miRNAs. Here we show distinctive patterns of sRNA expression that often coincide with stages of the developmental process such as flowering, early and late fruit maturation. Moreover, thousands of non-miRNA sRNAs are differentially expressed during fruit development and ripening. Some of these differentially expressed sRNAs derived from transposons but many derive from protein coding genes or regions that show homology to protein coding genes, several of which are known to play a role in flower and fruit development. These findings raise the possibility of a regulative role of these sRNAs during fruit onset and maturation in a crop species. We also identified six new miRNAs and experimentally validated two target mRNAs. These two mRNAs are targeted by the same miRNA but do not belong to the same gene family, which is rare for plant miRNAs. Expression pattern and putative function of these targets indicate a possible role in glutamate accumulation, which contributes to establishing the taste of the fruit.

Deep sequencing on genome-wide scale reveals the unique composition and expression patterns of microRNAs in developing pollen of Oryza sativa

Background

Pollen development in flowering plants requires strict control of the gene expression program and genetic information stability by mechanisms possibly including the miRNA pathway. However, our understanding of the miRNA pathway in pollen development remains limited, and the dynamic profile of miRNAs in developing pollen is unknown.

Results

Using next-generation sequencing technology, we pyrosequenced small RNA populations from rice uninucleate microspores to tricellular pollen and control sporophytic tissues at the genome-wide level. We identified 292 known miRNAs, including members of all 20 families conserved in plants, and 75 novel miRNAs. Of the 292 known miRNAs, 202 were expressed, with 103 enriched, in developing pollen. More than half of these novel miRNAs displayed pollen-or stage-specific expression. Furthermore, analyzing the 367 miRNAs and their predicted targets indicated that correlation in expression profiles of pollen-enriched known miRNAs and their targets significantly differs from that of sporophyte-enriched known miRNAs and their targets in some functional terms, while novel miRNAs appeared to negatively regulate their targets. Importantly, gene ontology abundance analysis demonstrated chromatin assembly and disassembly was important in the targets of bicellular pollen-expressed novel miRNAs. Principal component analysis revealed pollen of all three stages was discriminated from sporophytes, largely because of the novel and non-conserved known miRNAs.

Conclusions

Our study, for the first time, revealed the differences in composition and expression profiles of miRNAs between developing pollen and sporophytes, with novel and non-conserved known miRNAs the main contributors. Our results suggest the important roles of the miRNA pathway in pollen development.

Discovery of barley miRNAs through deep sequencing of short reads

Background

MicroRNAs are important components of the regulatory network of biological systems and thousands have been discovered in both animals and plants. Systematic investigations performed in species with sequenced genomes such as Arabidopsis, rice, poplar and Brachypodium have provided insights into the evolutionary relationships of this class of small RNAs among plants. However, miRNAs from barley, one of the most important cereal crops, remain unknown.

Results

We performed a large scale study of barley miRNAs through deep sequencing of small RNAs extracted from leaves of two barley cultivars. By using the presence of miRNA precursor sequences in related genomes as one of a number of supporting criteria, we identified up to 100 miRNAs in barley. Of these only 56 have orthologs in wheat, rice or Brachypodium that are known to be expressed, while up to 44 appear to be specifically expressed in barley.

Conclusions

Our study, the first large scale investigation of small RNAs in barley, has identified up to 100 miRNAs. We demonstrate that reliable identification of miRNAs via deep sequencing in a species whose genome has not been sequenced requires a more careful analysis of sequencing errors than is commonly performed. We devised a read filtering procedure for dealing with errors. In addition, we found that the use of a large dataset of almost 35 million reads permits the use of read abundance distributions along putative precursor sequences as a practical tool for isolating miRNAs in a large background of reads originating from other non-coding and coding RNAs. This study therefore provides a generic approach for discovering novel miRNAs where no genome sequence is available.

Small RNA discovery and characterisation in eukaryotes using high-throughput approaches

RNA silencing is a mechanism of genetic regulation that is mediated by short noncoding RNAs, or small RNAs (sRNAs). Regulatory interactions are established based on nucleotide sequence complementarity between the sRNAs and their targets. The development of new high-throughput sequencing technologies has accelerated the discovery of sRNAs in a variety of plants and animals. The use of these and other high-throughput technologies, such as microarrays, to measure RNA and protein concentrations of gene products potentially regulated by sRNAs has also been important for their functional characterisation. mRNAs targeted by sRNAs can produce new sRNAs or the protein encoded by the target mRNA can regulate other mRNAs. In either case the targeting sRNAs are parts of complex RNA networks therefore identifying and characterising sRNAs contribute to better understanding of RNA networks. In this chapter we will review RNA silencing, the different types of sRNAs that mediate it and the computational methods that have been developed to use high-throughput technologies in the study of sRNAs and their targets.

Small RNAs, DNA methylation and transposable elements in wheat

Background

More than 80% of the wheat genome is composed of transposable elements (TEs). Since active TEs can move to different locations and potentially impose a significant mutational load, their expression is suppressed in the genome via small non-coding RNAs (sRNAs). sRNAs guide silencing of TEs at the transcriptional (mainly 24-nt sRNAs) and post-transcriptional (mainly 21-nt sRNAs) levels. In this study, we report the distribution of these two types of sRNAs among the different classes of wheat TEs, the regions targeted within the TEs, and their impact on the methylation patterns of the targeted regions.

Results

We constructed an sRNA library from hexaploid wheat and developed a database that included our library and three other publicly available sRNA libraries from wheat. For five completely-sequenced wheat BAC contigs, most perfectly matching sRNAs represented TE sequences, suggesting that a large fraction of the wheat sRNAs originated from TEs. An analysis of all wheat TEs present in the Triticeae Repeat Sequence database showed that sRNA abundance was correlated with the estimated number of TEs within each class. Most of the sRNAs perfectly matching miniature inverted repeat transposable elements (MITEs) belonged to the 21-nt class and were mainly targeted to the terminal inverted repeats (TIRs). In contrast, most of the sRNAs matching class I and class II TEs belonged to the 24-nt class and were mainly targeted to the long terminal repeats (LTRs) in the class I TEs and to the terminal repeats in CACTA transposons. An analysis of the mutation frequency in potentially methylated sites revealed a three-fold increase in TE mutation frequency relative to intron and untranslated genic regions. This increase is consistent with wheat TEs being preferentially methylated, likely by sRNA targeting.

Conclusions

Our study examines the wheat epigenome in relation to known TEs. sRNA-directed transcriptional and post-transcriptional silencing plays important roles in the short-term suppression of TEs in the wheat genome, whereas DNA methylation and increased mutation rates may provide a long-term mechanism to inactivate TEs.