Small RNAs in angiosperms: sequence characteristics, distribution and generation

Small RNAs target native and cross-kingdom transcripts on both sides of the wheat stripe rust interaction

The wheat stripe rust fungus (Puccinia striiformis f.sp. tritici) threatens global wheat production. Small RNAs (sRNAs) modulate plant defense induction, and RNA exchange between host and microbe causes cross-kingdom gene silencing, but few examples are known in rust fungi. This study combined sRNA, parallel analysis of RNA ends, and gene expression data to discover sRNA-target pairs on each side of the interaction. Specific wheat 24 nt sRNAs were suppressed, while particular 35 nt fragments were strongly induced upon infection. Wheat sRNAs cleaved fungal transcripts coding for a ribosomal protein and a glycosyl hydrolase effector. Fungal microRNA-like and phased 21 nt sRNAs originated from long inverted repeats near protein coding genes. Fungal sRNAs targeted native transcripts: transposons and kinases; and cross-kingdom transcripts: a wheat nucleotide-binding domain leucine-rich repeat receptor (NLR) and multiple defense-related transcription factor families. This work sheds light on host-microbe coevolution and delivers prospects for developing pathogen control biotechnology.

Atypical DNA methylation, sRNA-size distribution, and female gametogenesis in Utricularia gibba

The most studied DNA methylation pathway in plants is the RNA Directed DNA Methylation (RdDM), a conserved mechanism that involves the role of noncoding RNAs to control the expansion of the noncoding genome. Genome-wide DNA methylation levels have been reported to correlate with genome size. However, little is known about the catalog of noncoding RNAs and the impact on DNA methylation in small plant genomes with reduced noncoding regions. Because of the small length of intergenic regions in the compact genome of the carnivorous plant Utricularia gibba, we investigated its repertoire of noncoding RNA and DNA methylation landscape. Here, we report that, compared to other angiosperms, U. gibba has an unusual distribution of small RNAs and reduced global DNA methylation levels. DNA methylation was determined using a novel strategy based on long-read DNA sequencing with the Pacific Bioscience platform and confirmed by whole-genome bisulfite sequencing. Moreover, some key genes involved in the RdDM pathway may not represented by compensatory paralogs or comprise truncated proteins, for example, U. gibba DICER-LIKE 3 (DCL3), encoding a DICER endonuclease that produces 24-nt small-interfering RNAs, has lost key domains required for complete function. Our results unveil that a truncated DCL3 correlates with a decreased proportion of 24-nt small-interfering RNAs, low DNA methylation levels, and developmental abnormalities during female gametogenesis in U. gibba. Alterations in female gametogenesis are reminiscent of RdDM mutant phenotypes in Arabidopsis thaliana. It would be interesting to further study the biological implications of the DCL3 truncation in U. gibba, as it could represent an initial step in the evolution of RdDM pathway in compact genomes.

Comparative analysis of sRNAs, degradome and transcriptomics in sweet sorghum reveals the regulatory roles of miRNAs in Cd accumulation and tolerance

Key miRNAs including sbi-miR169p/q, sbi-miR171g/j, sbi-miR172a/c/d, sbi-miR172e, sbi-miR319a/b, sbi-miR396a/b, miR408, sbi-miR5384, sbi-miR5565e and nov_23 were identified to function in the regulation of Cd accumulation and tolerance. As an energy plant, sweet sorghum shows great potential in the phytoremediation of Cd-contaminated soils. However, few studies have focused on the regulatory roles of miRNAs and their targets under Cd stress. In this study, comparative analysis of sRNAs, degradome and transcriptomics was conducted in high-Cd accumulation (H18) and low-Cd accumulation (L69) genotypes of sweet sorghum. A total of 38 conserved and 23 novel miRNAs with differential expressions were identified under Cd stress or between H18 and L69, and 114 target genes of 41 miRNAs were validated. Furthermore, 25 miRNA-mRNA pairs exhibited negatively correlated expression profiles and sbi-miR172e together with its target might participate in the distinct Cd tolerance between H18 and L69 as well as sbi-miR172a/c/d. Additionally, two groups of them: miR169p/q-nov_23 and miR408 were focused through the co-expression analysis, which might be involved in Cd uptake and tolerance by regulating their targets associated with transmembrane transportation, cytoskeleton activity, cell wall construction and ROS (reactive oxygen species) homeostasis. Further experiments exhibited that cell wall components of H18 and L69 were different when exposed to cadmium, which might be regulated by miR169p/q, miR171g/j, miR319a/b, miR396a/b, miR5384 and miR5565e through their targets. Through this study, we aim to reveal the potential miRNAs involved in sweet sorghum in response to Cd stress and provide references for developing high-Cd accumulation or high Cd-resistant germplasm of sweet sorghum that can be used in phytoremediation.

Secondary siRNAs in Plants: Biosynthesis, Various Functions, and Applications in Virology

The major components of RNA silencing include both transitive and systemic small RNAs, which are technically called secondary sRNAs. Double-stranded RNAs trigger systemic silencing pathways to negatively regulate gene expression. The secondary siRNAs generated as a result of transitive silencing also play a substantial role in gene silencing especially in antiviral defense. In this review, we first describe the discovery and pathways of transitivity with emphasis on RNA-dependent RNA polymerases followed by description on the short range and systemic spread of silencing. We also provide an in-depth view on the various size classes of secondary siRNAs and their different roles in RNA silencing including their categorization based on their biogenesis. The other regulatory roles of secondary siRNAs in transgene silencing, virus-induced gene silencing, transitivity, and trans-species transfer have also been detailed. The possible implications and applications of systemic silencing and the different gene silencing tools developed are also described. The details on mobility and roles of secondary siRNAs derived from viral genome in plant defense against the respective viruses are presented. This entails the description of other compatible plant-virus interactions and the corresponding small RNAs that determine recovery from disease symptoms, exclusion of viruses from shoot meristems, and natural resistance. The last section presents an overview on the usefulness of RNA silencing for management of viral infections in crop plants.

Landscape of Fluid Sets of Hairpin-Derived 21-/24-nt-Long Small RNAs at Seed Set Uncovers Special Epigenetic Features in Picea glauca

Conifers’ exceptionally large genome (20–30 Gb) is scattered with 60% retrotransposon (RT) components and we have little knowledge on their origin and evolutionary implications. RTs may impede the expression of flanking genes and provide sources of the formation of novel small RNA (sRNAs) populations to constrain events of transposon (TE) proliferation/transposition. Here we show a declining expression of 24-nt-long sRNAs and low expression levels of their key processing gene, pgRTL2 (RNASE THREE LIKE 2) at seed set in Picea glauca. The sRNAs in 24-nt size class are significantly less enriched in type and read number than 21-nt sRNAs and have not been documented in other species. The architecture of MIR loci generating highly expressed 24-/21-nt sRNAs is featured by long terminal repeat—retrotransposons (LTR-RTs) in families of Ty3/Gypsy and Ty1/Copia elements. This implies that the production of sRNAs may be predominantly originated from TE fragments on chromosomes. Furthermore, a large proportion of highly expressed 24-nt sRNAs does not have predictable targets against unique genes in Picea, suggestive of their potential pathway in DNA methylation modifications on, for instance, TEs. Additionally, the classification of computationally predicted sRNAs suggests that 24-nt sRNA targets may bear particular functions in metabolic processes while 21-nt sRNAs target genes involved in many different biological processes. This study, therefore, directs our attention to a possible extrapolation that lacking of 24-nt sRNAs at the late conifer seed developmental phase may result in less constraints in TE activities, thus contributing to the massive expansion of genome size.

Regulatory crosstalk between microRNAs and hormone signalling cascades controls the variation on seed dormancy phenotype at Arabidopsis thaliana seed set

We employed an Illumina sequencing approach to identify candidate microRNA cohorts that may greatly contribute to seed dormancy modulation and to construct a microRNA-gene regulatory network in hormone signalling cascades. MicroRNAs (miRNAs) are important signalling molecules and regulate many developmental programs of plants. Some miRNAs have been integrated into gene regulatory networks (GRNs) and coordinate developmental plasticity, but few study systematically investigated how phenotypical variations are regulated through differential expression of miRNA tags in GRNs during seed set. Using 'top-down' analyses (i.e., identify miRNAs associated with known phenotypical variations), we chose two Arabidopsis ecotypes (Cvi-0 and Col-0) with contrasting seed dormancy and sequenced miRNA reads in the first ten phases at seed set. We computationally predicted target genes of miRNAs and implemented statistical analyses for normalized relative expression of top abundant miRNA cohorts between the two ecotypes. We especially focused on miRNA cohorts targeting mRNAs encoding transcription factors in hormone signalling cascades. We report, with high confidence hits, that a cohort of 14 miRNAs (miR-156b, -159b, -160, -161*, -319a, -390a, -396, -773a, -779, -842, -852, -859, -1886*, and a novel sequence in miR8172 family) may greatly contribute to seed dormancy modulation, of which seven are involved in hormone signalling cascades. Moreover, their expression patterns indicated that 5 ± 1 days after flowering (at embryogenesis-to-maturation transition) is a critical phase at seed set. This study reinforces the notion that miRNAs that regulate seed dormancy modulation and provides a novel paradigm of studying the correlation between genotypes (miRNAs) and phenotypes.

Long non-coding RNAs are major contributors to transcriptome changes in sunflower meiocytes with different recombination rates

BACKGROUND

Meiosis is a form of specialized cell division that marks the transition from diploid meiocyte to haploid gamete, and provides an opportunity for genetic reassortment through recombination. Experimental data indicates that, relative to their wild ancestors, cultivated sunflower varieties show a higher recombination rate during meiosis. To better understand the molecular basis for this difference, we compared gene expression in male sunflower meiocytes in prophase I isolated from a domesticated line, a wild relative, and a F1 hybrid of the two.

RESULTS

Of the genes that showed differential expression between the wild and domesticated genotypes, 63.62 % could not be identified as protein-coding genes, and of these genes, 70.98 % passed stringent filters to be classified as long non-coding RNAs (lncRNAs). Compared to the sunflower somatic transcriptome, meiocytes express a higher proportion of lncRNAs, and the majority of genes with exclusive expression in meiocytes were lncRNAs. Around 40 % of the lncRNAs showed sequence similarity with small RNAs (sRNA), while 1.53 % were predicted to be sunflower natural antisense transcripts (NATs), and 9.18 % contained transposable elements (TE). We identified 6895 lncRNAs that are exclusively expressed in meiocytes, these lncRNAs appear to have higher conservation, a greater degree of differential expression, a higher proportion of sRNA similarity, and higher TE content relative to lncRNAs that are also expressed in the somatic transcriptome.

CONCLUSIONS

lncRNAs play important roles in plant meiosis and may participate in chromatin modification processes, although other regulatory functions cannot be excluded. lncRNAs could also be related to the different recombination rates seen for domesticated and wild sunflowers.

NATpipe: an integrative pipeline for systematical discovery of natural antisense transcripts (NATs) and phase-distributed nat-siRNAs from de novo assembled transcriptomes

Nat-siRNAs (small interfering RNAs originated from natural antisense transcripts) are a class of functional small RNA (sRNA) species discovered in both plants and animals. These siRNAs are highly enriched within the annealed regions of the NAT (natural antisense transcript) pairs. To date, great research efforts have been taken for systematical identification of the NATs in various organisms. However, developing a freely available and easy-to-use program for NAT prediction is strongly demanded by researchers. Here, we proposed an integrative pipeline named NATpipe for systematical discovery of NATs from de novo assembled transcriptomes. By utilizing sRNA sequencing data, the pipeline also allowed users to search for phase-distributed nat-siRNAs within the perfectly annealed regions of the NAT pairs. Additionally, more reliable nat-siRNA loci could be identified based on degradome sequencing data. A case study on the non-model plant Dendrobium officinale was performed to illustrate the utility of NATpipe. Finally, we hope that NATpipe would be a useful tool for NAT prediction, nat-siRNA discovery, and related functional studies. NATpipe is available at www.bioinfolab.cn/NATpipe/NATpipe.zip.

Genome-wide discovery and validation of Eucalyptus small RNAs reveals variable patterns of conservation and diversity across species of Myrtaceae

Background

Micro RNAs are a class of small non coding RNAs of 20–24 nucleotides transcribed as single stranded precursors from MIR gene loci. Initially described as post-transcriptional regulators involved in development, two decades ago, miRNAs have been proven to regulate a wide range of processes in plants such as germination, morphology and responses to biotic and abiotic stress. Despite wide conservation in plants, a number of miRNAs are lineage specific. We describe the first genome wide survey of Eucalyptus miRNAs based on high throughput sequencing.

Results

In addition to discovering small RNA sequences, MIR loci were mapped onto the reference genome and interspecific variability investigated. Sequencing was carried out for the two most world widely planted species, E. grandis and E. globulus. To maximize discovery, E. grandis samples were from BRASUZ1, the same tree whose genome provided the reference sequence. Interspecific analysis reinforces the variability in small RNA repertoire even between closely related species. Characterization of Eucalyptus small RNA sequences showed 95 orthologous to conserved miRNAs and 193 novel miRNAs. In silico target prediction confirmed 163 novel miRNAs and degradome sequencing experimentally confirmed several hundred targets. Experimental evidence based on the exclusive expression of a set of small RNAs across 16 species within Myrtaceae further highlighted variable patterns of conservation and diversity of these regulatory elements.

Conclusions

The description of miRNAs in Eucalyptus contributes to scientific knowledge of this vast genre, which is the most widely planted hardwood crop in the tropical and subtropical world, adding another important element to the annotation of Eucalyptus grandis reference genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2322-6) contains supplementary material, which is available to authorized users.

Identification of mRNA-like non-coding RNAs and validation of a mighty one named MAR in Panax ginseng

Increasing evidence suggests that long non-coding RNAs (lncRNAs) play significant roles in plants. However, little is known about lncRNAs in Panax ginseng C. A. Meyer, an economically significant medicinal plant species. A total of 3,688 mRNA-like non-coding RNAs (mlncRNAs), a class of lncRNAs, were identified in P. ginseng. Approximately 40% of the identified mlncRNAs were processed into small RNAs, implying their regulatory roles via small RNA-mediated mechanisms. Eleven miRNA-generating mlncRNAs also produced siRNAs, suggesting the coordinated production of miRNAs and siRNAs in P. ginseng. The mlncRNA-derived small RNAs might be 21-, 22-, or 24-nt phased and could be generated from both or only one strand of mlncRNAs, or from super long hairpin structures. A full-length mlncRNA, termed MAR (multiple-function-associated mlncRNA), was cloned. It generated the most abundant siRNAs. The MAR siRNAs were predominantly 24-nt and some of them were distributed in a phased pattern. A total of 228 targets were predicted for 71 MAR siRNAs. Degradome sequencing validated 68 predicted targets involved in diverse metabolic pathways, suggesting the significance of MAR in P. ginseng. Consistently, MAR was detected in all tissues analyzed and responded to methyl jasmonate (MeJA) treatment. It sheds light on the function of mlncRNAs in plants.

GC content fluctuation around plant small RNA-generating sites

GC content of small RNA-generating sites and their flanking sequences in Arabidopsis thaliana and rice was systematically analyzed in silico. High GC content fluctuation (GCF) is observed at the borders of sRNA sites, while the GCF within sRNA sites is low. Furthermore, the GC content along sequences of some miniature inverted-repeat transposable element (MITE) families coincides with the abundance of MITE-derived small RNAs. The GCF within tasiRNA clusters is negatively correlated with its phasing score. We conclude that high GC content and low GCF could increase the expression of small RNA. Our results provide further insights on small RNA expression, which may be applied to improve the silencing efficiency of RNAi and virus-induced gene silencing.

Mechanisms of small RNA generation from cis-NATs in response to environmental and developmental cues

A large proportion of eukaryotic genomes is transcribed from both positive and negative strands of DNA and thus may generate overlapping sense and antisense transcripts. Some of these so-called natural antisense transcripts (NATs) are possibly co-regulated. When the overlapping sense and antisense transcripts are expressed at the same time in the same cell in response to various developmental and environmental cues; they may form double-stranded RNAs, which could be recognized by the small RNA biogenesis machinery and processed into small interfering RNAs (siRNAs). cis-NAT-derived siRNAs (nat-siRNAs) are present in plants, animals, and fungi. In plants, the presence of nat-siRNAs is supported not only by Northern blot and genetic analyses, but also by the fact that there is an overall sixfold enrichment of siRNAs in the overlapping regions of cis-NATs and 19%-29% of the siRNA-generating cis-NATs in plants give rise to siRNAs only in their overlapping regions. Silencing mediated by nat-siRNAs is one of the mechanisms for regulating the expression of the cis-NATs. This review focuses on challenging issues related to the biogenesis mechanisms as well as regulation and detection of nat-siRNAs. The advantages and limitations of new technologies for detecting cis-NATs, including direct RNA sequencing and strand-specific RNA sequencing, are also discussed.

Prediction and identification of natural antisense transcripts and their small RNAs in soybean (Glycine max)

BACKGROUND

Natural antisense transcripts (NATs) are a class of RNAs that contain a sequence complementary to other transcripts. NATs occur widely in eukaryotes and play critical roles in post-transcriptional regulation. Soybean NAT sequences are predicted in the PlantNATsDB, but detailed analyses of these NATs remain to be performed.

RESULTS

A total of 26,216 NATs, including 994 cis-NATs and 25,222 trans-NATs, were predicted in soybean. Each sense transcript had 1-177 antisense transcripts. We identified 21 trans-NATs using RT-PCR amplification. Additionally, we identified 179 cis-NATs and 6,629 trans-NATs that gave rise to small RNAs; these were enriched in the NAT overlapping region. The most abundant small RNAs were 21, 22, and 24 nt in length. The generation of small RNAs was biased to one stand of the NATs, and the degradation of NATs was biased. High-throughput sequencing of the degradome allowed for the global identification of NAT small interfering RNAs (nat-siRNAs) targets. 446 target genes for 165 of these nat-siRNAs were identified. The nat-siRNA target could be one transcript of a given NAT, or from other gene transcripts. We identified five NAT transcripts containing a hairpin structure that is characteristic of pre-miRNA. We identified a total of 86 microRNA (miRNA) targets that had antisense transcripts in soybean.

CONCLUSIONS

We globally identified nat-siRNAs, and the targets of nat-siRNAs in soybean. It is likely that the cis-NATs, trans-NATs, nat-siRNAs, miRNAs, and miRNA targets form complex regulatory networks.

Genome-wide identification and analysis of microRNA responding to long-term waterlogging in crown roots of maize seedlings

MicroRNAs (miRNAs) are critical post-transcriptional modulators of gene expression involving in plant responses to abiotic stress. However, the regulation of miRNA in the morphological response to waterlogging is poorly understood in maize. In this study, we detected miRNAs and their targets that expressed in waterlogged crown roots of maize seedlings in two inbred lines (Hz32 and Mo17) by RNA sequencing. A total of 61 mature miRNAs were found including 36 known maize (zma) miRNAs and 25 potential novel miRNA candidates. Comparison of miRNA expression in both waterlogged and control crown roots revealed 32 waterlogging-responsive miRNAs, most were consistently downregulated under waterlogging in the two inbred lines. We identified the miRNA targets through degradome sequencing. Many known miRNA targets involving in transcription regulation and reactive oxygen species elimination were found in the degradome libraries, and 17 targets of 10 newly detected miRNAs were identified as well. Moreover, the miRNA-mediated pathways that respond to waterlogging and regulate the induction of crown roots were discussed. This study is a comprehensive survey of responsive miRNAs in waterlogged maize crown roots. The results will help to understand the miRNA expression in response to waterlogging and miRNA-mediated regulation of morphological adaptation to waterlogging in maize.

Construction of gene regulatory networks mediated by vegetative and reproductive stage-specific small RNAs in rice (Oryza sativa)

Although huge amounts of high-throughput sequencing (HTS) data are available, limited systematic analyses have been performed by integrating these valuable resources. Based on small RNA (sRNA), RNA and degradome HTS data, the sRNAs specifically expressed at vegetative and reproductive stages were identified separately in rice. Two distinct groups of sRNA HTS data, which were prepared during the vegetative and the reproductive stages, were utilized to extract stage-specific sRNAs. Degradome sequencing data were employed for sRNA target validation. RNA sequencing data were used to construct expression-based, sRNA-mediated networks. As a result, 26 microRNAs and 413 sRNAs were specifically expressed at the vegetative stage, and 79 microRNAs and 539 sRNAs were specifically expressed at the reproductive stage. In addition to the microRNAs, numerous stage-specific sRNAs enriched in ARGONAUTE1 showed great potential to perform cleavage-based repression on the targets. Several stage-specific sRNAs were indicated to result from the wobble effect of Dicer-like 1-mediated processing of microRNA precursors. The expression patterns of the sRNA targets, and the stage-specific cleavage signals strongly indicated the reliability of the constructed networks. A set of rice stage-specific sRNAs along with the regulatory cascades, which have great potential in regulating specific developmental stages, were provided for further investigation.

PlantNATsDB: a comprehensive database of plant natural antisense transcripts

Natural antisense transcripts (NATs), as one type of regulatory RNAs, occur prevalently in plant genomes and play significant roles in physiological and pathological processes. Although their important biological functions have been reported widely, a comprehensive database is lacking up to now. Consequently, we constructed a plant NAT database (PlantNATsDB) involving approximately 2 million NAT pairs in 69 plant species. GO annotation and high-throughput small RNA sequencing data currently available were integrated to investigate the biological function of NATs. PlantNATsDB provides various user-friendly web interfaces to facilitate the presentation of NATs and an integrated, graphical network browser to display the complex networks formed by different NATs. Moreover, a ‘Gene Set Analysis’ module based on GO annotation was designed to dig out the statistical significantly overrepresented GO categories from the specific NAT network. PlantNATsDB is currently the most comprehensive resource of NATs in the plant kingdom, which can serve as a reference database to investigate the regulatory function of NATs. The PlantNATsDB is freely available at http://bis.zju.edu.cn/pnatdb/.

Plant siRNAs from introns mediate DNA methylation of host genes

Small RNAs (sRNAs), largely known as microRNAs (miRNAs) and short interfering RNAs (siRNAs), emerged as the critical components of genetic and epigenetic regulation in eukaryotic genomes. In animals, a sizable portion of miRNAs reside within the introns of protein-coding genes, designated as mirtron genes. Recently, high-throughput sequencing (HTS) revealed a huge amount of sRNAs that derived from introns in plants, such as the monocot rice (Oryza sativa). However, the biogenesis and the biological functions of this kind of sRNAs remain elusive. Here, we performed a genome-scale survey of intron-derived sRNAs in rice based on HTS data. Several introns were found to have great potential to form internal hairpin structures, and the short hairpins could generate miRNAs while the larger ones could produce siRNAs. Furthermore, 22 introns, termed "sirtrons," were identified from the rice protein-coding genes. The single-stranded sirtrons produced a diverse set of siRNAs from long hairpin structures. These sirtron-derived siRNAs are dominantly 21 nt, 22 nt, and 24 nt in length, whose production relied on DCL4, DCL2, and DCL3, respectively. We also observed a strong tendency for the sirtron-derived siRNAs to be coexpressed with their host genes. Finally, the 24-nt siRNAs incorporated with Argonaute 4 (AGO4) could direct DNA methylation on their host genes. In this regard, homeostatic self-regulation between intron-derived siRNAs and their host genes was proposed.

The Selaginella genome identifies genetic changes associated with the evolution of vascular plants

Vascular plants appeared ~410 million years ago, then diverged into several lineages of which only two survive: the euphyllophytes (ferns and seed plants) and the lycophytes. We report here the genome sequence of the lycophyte Selaginella moellendorffii (Selaginella), the first nonseed vascular plant genome reported. By comparing gene content in evolutionarily diverse taxa, we found that the transition from a gametophyte- to a sporophyte-dominated life cycle required far fewer new genes than the transition from a nonseed vascular to a flowering plant, whereas secondary metabolic genes expanded extensively and in parallel in the lycophyte and angiosperm lineages. Selaginella differs in posttranscriptional gene regulation, including small RNA regulation of repetitive elements, an absence of the trans-acting small interfering RNA pathway, and extensive RNA editing of organellar genes.

A simple and efficient method for isolating small RNAs from different plant species

BACKGROUND

Small RNAs emerged over the last decade as key regulators in diverse biological processes in eukaryotic organisms. To identify and study small RNAs, good and efficient protocols are necessary to isolate them, which sometimes may be challenging due to the composition of specific tissues of certain plant species. Here we describe a simple and efficient method to isolate small RNAs from different plant species.

RESULTS

We developed a simple and efficient method to isolate small RNAs from different plant species by first comparing different total RNA extraction protocols, followed by streamlining the best one, finally resulting in a small RNA extraction method that has no need of first total RNA extraction and is not based on the commercially available TRIzol® Reagent or columns. This small RNA extraction method not only works well for plant tissues with high polysaccharide content, like cactus, agave, banana, and tomato, but also for plant species like Arabidopsis or tobacco. Furthermore, the obtained small RNA samples were successfully used in northern blot assays.

CONCLUSION

Here we provide a simple and efficient method to isolate small RNAs from different plant species, such as cactus, agave, banana, tomato, Arabidopsis, and tobacco, and the small RNAs from this simplified and low cost method is suitable for downstream handling like northern blot assays.