Genome-wide profiling of populus small RNAs

Exploring the drought-responsive miRNAs and their corresponding target genes in chickpea root tissue

Chickpea is an important pulse crop globally, with major production in Southeast Asia. However, the production of chickpea is hampered due to various biotic and abiotic stressors. In response to such stressors, microRNAs which are small non-coding regulatory RNA molecules have been observed as key players. The present study evaluates the role of drought-responsive microRNAs in the root tissues of chickpea genotypes contrasting for drought tolerance. This study led to the generation of 146.7 million short-read sequences from small RNA libraries constructed from the root tissues of the two genotypes. Upon analysis, 224 conserved and 155 novel miRNA sequences were identified. The miR156 family was found to be the most abundant among the 51 families identified for the conserved miRNAs. Quantitative real-time PCR (qRT-PCR) was used to validate eleven conserved and six novel miRNAs. The identification of drought-induced expression of specific miRNAs and their related target genes suggests miRNA-mediated response mechanisms in chickpea. Furthermore, this research investigated the role of drought-responsive miRNAs, specifically miR171 and miR166 and their target genes, SCL27 (scarecrow-like protein 27) and ATHB15 (Homeobox-leucine zipper family protein), respectively. The study validated the miR171 and miR166 directed cleavage of SCL27 and ATHB15, respectively, in drought-stressed root tissues using 5´RLM-RACE (5' RNA Ligase-Mediated Rapid Amplification of cDNA Ends) analysis. The study highlights the role of diverse miRNAs in chickpea for mitigating drought.

Genome-Wide Identification and Functional Analysis of Salvia miltiorrhiza MicroRNAs Reveal the Negative Regulatory Role of Smi-miR159a in Phenolic Acid Biosynthesis

MicroRNAs (miRNAs) are a group of endogenous small non-coding RNAs in plants. They play critical functions in various biological processes during plant growth and development. Salvia miltiorrhiza is a well-known traditional Chinese medicinal plant with significant medicinal, economic, and academic values. In order to elucidate the role of miRNAs in S. miltiorrhiza, six small RNA libraries from mature roots, young roots, stems, mature leaves, young leaves and flowers of S. miltiorrhiza and one degradome library from mixed tissues were constructed. A total of 184 miRNA precursors, generating 137 known and 49 novel miRNAs, were genome-widely identified. The identified miRNAs were predicted to play diversified regulatory roles in plants through regulating 891 genes. qRT-PCR and 5' RLM-RACE assays validated the negative regulatory role of smi-miR159a in SmMYB62, SmMYB78, and SmMYB80. To elucidate the function of smi-miR159a in bioactive compound biosynthesis, smi-miR159a transgenic hairy roots were generated and analyzed. The results showed that overexpression of smi-miR159a caused a significant decrease in rosmarinic acid and salvianolic acid B contents. qRT-PCR analysis showed that the targets of smi-miR159a, including SmMYB62, SmMYB78, and SmMYB80, were significantly down-regulated, accompanied by the down-regulation of SmPAL1, SmC4H1, Sm4CL1, SmTAT1, SmTAT3, SmHPPR1, SmRAS, and SmCYP98A14 genes involved in phenolic acid biosynthesis. It suggests that smi-miR159a is a significant negative regulator of phenolic acid biosynthesis in S. miltiorrhiza.

Identification of microRNA and analysis of target genes in Panax ginseng

Objective

Ginsenosides, polysaccharides and phenols, the main active ingredients in Panax ginseng, are not different significantly in content between 3 and 5 years old of ginsengs called Yuan ginseng and more than ten years old ones called Shizhu ginseng. The responsible chemical compounds cannot fully explain difference in efficacy between them. According to reports in Lonicerae Japonicae Flos (Jinyinhua in Chinese) and Glycyrrhizae Radix et Rhizoma (Gancao in Chinese), microRNA may play a role in efficacy, so we identified microRNAs in P. ginseng at the different growth years and analyzed their target genes.

Methods

Using high-throughput sequencing, the RNA-Seq, small RNA-Seq and degradome databases of P. ginseng were constructed. The differentially expressed microRNAs was identified by qRT-PCR.

Results

A total of 63,875 unigenes and 24,154,579 small RNA clean reads were obtained from the roots of P. ginseng. From these small RNAs, 71 miRNA families were identified by bioinformatics target prediction software, including 34 conserved miRNAs, 37 non-conserved miRNA families, as well as 179 target genes of 17 known miRNAs. Through degradome sequencing and computation, we finally verified 13 targets of eight miRNAs involved in transcription, energy metabolism, biological stress and disease resistance, suggesting the significance of miRNAs in the development of P. ginseng. Consistently, major miRNA targets exhibited tissue specificity and complexity in expression patterns.

Conclusion

Differential expression microRNAs were found in different growth years of ginsengs (Shizhu ginseng and Yuan ginseng), and the regulatory roles and functional annotations of miRNA targets in P. ginseng need further investigation.

Phaseolus vulgaris MIR1511 genotypic variations differentially regulate plant tolerance to aluminum toxicity

The common-bean (Phaseolus vulgaris), a widely consumed legume, originated in Mesoamerica and expanded to South America, resulting in the development of two geographically distinct gene pools. Poor soil condition, including metal toxicity, are often constraints to common-bean crop production. Several P. vulgaris miRNAs, including miR1511, respond to metal toxicity. The MIR1511 gene sequence from the two P. vulgaris model sequenced genotypes revealed that, as opposed to BAT93 (Mesoamerican), the G19833 (Andean) accession displays a 58-bp deletion, comprising the mature and star miR1511 sequences. Genotyping-By-Sequencing data analysis from 87 non-admixed Phaseolus genotypes, comprising different Phaseolus species and P. vulgaris populations, revealed that all the P. vulgaris Andean genotypes and part of the Mesoamerican (MW1) genotypes analyzed displayed a truncated MIR1511 gene. The geographic origin of genotypes with a complete versus truncated MIR1511 showed a distinct distribution. The P. vulgaris ALS3 (Aluminum Sensitive Protein 3) gene, known to be important for aluminum detoxification in several plants, was experimentally validated as the miR1511 target. Roots from BAT93 plants showed decreased miR1511 and increased ALS3 transcript levels at early stages under aluminum toxicity (AlT), while G19833 plants, lacking mature miR1511, showed higher and earlier ALS3 response. Root architecture analyses evidenced higher tolerance of G19833 plants to AlT. However, G19833 plants engineered for miR1511 overexpression showed lower ALS3 transcript level and increased sensitivity to AlT. Absence of miR1511 in Andean genotypes, resulting in a diminished ALS3 transcript degradation, appears to be an evolutionary advantage to high Al levels in soils with increased drought conditions.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Newly identified miRNAs may contribute to aerenchyma formation in sugarcane roots

Small RNAs comprise three families of noncoding regulatory RNAs that control gene expression by blocking mRNA translation or leading to mRNA cleavage. Such post-transcriptional negative regulation is relevant for both plant development and environmental adaptations. An important biotechnological application of miRNA identification is the discovery of regulators and effectors of cell wall degradation, which can improve/facilitate hydrolysis of cell wall polymers for second-generation bioethanol production. The recent characterization of plant innate cell wall modifications occurring during root aerenchyma development triggered by ethylene led to the possibility of prospection for mechanisms of cell wall disassembly in sugarcane. By using next-generation sequencing, 39 miRNAs were identified in root segments along the process of aerenchyma development. Among them, 31 miRNAs were unknown to the sugarcane miRBase repository but previously identified as produced by its relative Sorghum bicolor. Key putative targets related to signal transduction, carbohydrate metabolic process, and cell wall organization or biogenesis were among the most representative gene categories targeted by miRNA. They belong to the subclasses of genes associated with the four modules of cell wall modification in sugarcane roots: cell expansion, cell separation, hemicellulose, and cellulose hydrolysis. Thirteen miRNAs possibly related to ethylene perception and signaling were also identified. Our findings suggest that miRNAs may be involved in the regulation of cell wall degradation during aerenchyma formation. This work also points out to potential molecular tools for sugarcane improvement in the context of second-generation biofuels.

Transcriptome-wide identification and profiling of miRNAs in a stress-tolerant conifer Sabina chinensis

miRNAs are important regulatory components involving in many biological processes, including plant development, vegetative and reproductive growth, and stress response. However, identification and characterization of miRNAs still remain limited for conifer species. In this study, with deep sequencing, we obtained 1,314,450 unique reads with 18-30 nt length from a stress-tolerant conifer, Sabina chinensis. We identified 37 conserved and 103 novel miRNAs, their unique characteristics were further analyzed, and 10 randomly selected were validated by qRT-PCR. Through miRNA target predictions and annotations, we found miRNA may have several targets as well a target could be regulated by several miRNAs, and a total of 2,397 mRNAs were predicted to be targets of the 140 miRNAs. These targets included not only important transcription factors such as auxin response factors, but also indispensable non-transcriptional factor proteins. Pathway-based analysis showed that S. chinensis miRNAs are involved in 172 metabolic pathways, of which 3 were discovered in adaptation-related pathways, indicating their possible relevance to the species' stress-tolerance characteristics. This study is expected to lay the foundation for exploring the regulative roles of miRNAs in development, growth, and response to environmental stresses of S. chinensis.

RNAi-Mediated Resistance Against Viruses in Perennial Fruit Plants

Small RNAs (sRNAs) are 20–30-nucleotide-long, regulatory, noncoding RNAs that induce silencing of target genes at the transcriptional and posttranscriptional levels. They are key components for cellular functions during plant development, hormone signaling, and stress responses. Generated from the cleavage of double-stranded RNAs (dsRNAs) or RNAs with hairpin structures by Dicer-like proteins (DCLs), they are loaded onto Argonaute (AGO) protein complexes to induce gene silencing of their complementary targets by promoting messenger RNA (mRNA) cleavage or degradation, translation inhibition, DNA methylation, and/or histone modifications. This mechanism of regulating RNA activity, collectively referred to as RNA interference (RNAi), which is an evolutionarily conserved process in eukaryotes. Plant RNAi pathways play a fundamental role in plant immunity against viruses and have been exploited via genetic engineering to control disease. Plant viruses of RNA origin that contain double-stranded RNA are targeted by the RNA-silencing machinery to produce virus-derived small RNAs (vsRNAs). Some vsRNAs serve as an effector to repress host immunity by capturing host RNAi pathways. High-throughput sequencing (HTS) strategies have been used to identify endogenous sRNA profiles, the “sRNAome”, and analyze expression in various perennial plants. Therefore, the review examines the current knowledge of sRNAs in perennial plants and fruits, describes the development and implementation of RNA interference (RNAi) in providing resistance against economically important viruses, and explores sRNA targets that are important in regulating a variety of biological processes.

Identification of novel miRNAs and their target genes from Populus szechuanica infected with Melampsora larici-populina

Two novel miRNAs were selected from a pre-constructed RNA library of Populus szechuanica infected with the foliar rust fungus Melampsora larici-populina in order to detect the genes regulated as targets of the miRNAs novel_mir_11 and novel_mir_357. The novel miRNAs were identified from P. szechuanica using stem-loop methods and their precursors were able to fold into a complete stem loop structure. The predicted target genes of the novel miRNAs were verified with RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-5'RACE). The full-length sequences of target genes, RPM1 and RPS2/5, in P. szechuanica were obtained through rapid amplification of cDNA ends (RACE) and officially named PsRPM1 and PsRPS2/5. These genes contain nucleotide binding site-leucine-rich repeats (NBS-LRR) domains typical of resistance genes. The expression levels of miRNAs and their target genes in different periods post infection were analysed with quantitative real-time PCR (qRT-PCR). After infection with the foliar rust fungus, the expression levels of the novel miRNAs and their target genes were dynamic. Both novel_mir_11 and novel_mir_357 negatively regulated the expression of their target genes. In this study, the regulatory effects of two novel miRNAs through their target genes were characterized to provide further mechanistic information regarding the interaction between Populus and a foliar rust fungus. Results of this study improve our understanding of the defence response mechanisms of Populus and will stimulate future work to characterize strategies to prevent and control Populus diseases.

Transcript-level expression control of plant NLR genes

Plant NLR genes encode sensitive immune receptors that can mediate the specific recognition of pathogen avirulence effectors and activate a strong defence response, termed effector-triggered immunity. The expression of NLRs requires strict regulation, as their ability to trigger immunity is dependent on their dose, and overexpression of NLRs results in autoimmunity and massive fitness costs. An elaborate interplay of different mechanisms controlling NLR transcript levels allows plants to maximize their defence capacity, whilst limiting negative impact on their fitness. Global suppression of NLR transcripts may be a prerequisite for the fast evolution of new NLR variants and the expansion of this gene family. Here, we summarize recent progress made towards a comprehensive understanding of NLR transcript-level expression control. Multiple mechanistic steps, including transcription as well as co-/post-transcriptional processing and transcript turn-over, contribute to balanced base levels of NLR transcripts and allow for dynamic adjustments to defence situations.

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.

MicroRNA profiling provides insights into post-transcriptional regulation of gene expression in chickpea root apex under salinity and water deficiency

Activity of root apical meristem (RAM) at the root apex is critical for stress-mediated modulation of root-architecture. Chickpea, like other legumes, possesses a basic open root meristem. Deep sequencing was used to perform microRNA expression profiling in root apex of chickpea (Cicer arietinum L.) in order to investigate post-transcriptional regulation of gene expression in this tissue in response to salinity and water deficit. Five small RNA libraries prepared from chickpea root apices at different stages of stress treatments were sequenced to obtain 284 unique miRNA sequences including 60 novel miRNAs belonging to total 255 families. Two hundred and fiftynine miRNAs were differentially expressed in stress. Six hundred and nine mRNA targets involved in diverse cellular processes were predicted for 244 miRNAs. Stress-responsive expression patterns of selected miRNAs, inverse expression patterns of their target genes and the target-cleavage sites were validated. Three candidate miRNA-target gene relationships were validated in transient expression system in chickpea. The miRNA expression profiling under salinity and water deficiency in a legume root apex and the reported function of their target genes suggested important roles of miRNA-mediated post-transcriptional regulation of gene expression involved in re-patterning of root hair cells, lateral root formation and high-affinity K+-uptake under these stresses.

Expression Variations of miRNAs and mRNAs in Rice (Oryza sativa)

Differences in expression levels are an important source of phenotypic variation within and between populations. MicroRNAs (miRNAs) are key players in post-transcriptional gene regulation that are important for plant development and stress responses. We surveyed expression variation of miRNAs and mRNAs of six accessions from two rice subspecies Oryza sativa L. ssp. indica and Oryza sativa L. ssp. japonica using deep sequencing. While more than half (53.7%) of the mature miRNAs exhibit differential expression between grains and seedlings of rice, only 11.0% show expression differences between subspecies, with an additional 2.2% differentiated for the development-by-subspecies interaction. Expression variation is greater for lowly conserved miRNAs than highly conserved miRNAs, whereas the latter show stronger negative correlation with their targets in expression changes between subspecies. Using a permutation test, we identified 51 miRNA–mRNA pairs that correlate negatively or positively in expression level among cultivated rice. Genes involved in various metabolic processes and stress responses are enriched in the differentially expressed genes between rice indica and japonica subspecies. Our results indicate that stabilizing selection is the major force governing miRNA expression in cultivated rice, albeit positive selection may be responsible for much of the between-subspecies expression divergence.

Causes and consequences of large clonal assemblies in a poplar hybrid zone

Asexual reproduction is a common and fundamental mode of reproduction in plants. Although persistence in adverse conditions underlies most known cases of clonal dominance, proximal genetic drivers remain unclear, in particular for populations dominated by a few large clones. In this study, we studied a clonal population of the riparian tree Populus alba in the Douro river basin (northwestern Iberian Peninsula) where it hybridizes with Populus tremula, a species that grows in highly contrasted ecological conditions. We used 73 nuclear microsatellites to test whether genomic background (species ancestry) is a relevant cause of clonal success, and to assess the evolutionary consequences of clonal dominance by a few genets. Additional genotyping-by-sequencing data were produced to estimate the age of the largest clones. We found that a few ancient (over a few thousand years old) and widespread genets dominate the population, both in terms of clone size and number of sexual offspring produced. Interestingly, large clones possessed two genomic regions introgressed from P. tremula, which may have favoured their spread under stressful environmental conditions. At the population level, the spread of large genets was accompanied by an overall ancient (>0.1 Myr) but soft decline of effective population size. Despite this decrease, and the high clonality and dominance of sexual reproduction by large clones, the Douro hybrid zone still displays considerable genetic diversity and low inbreeding. This suggests that even in extreme cases as in the Douro, asexual and sexual dominance of a few large, geographically extended individuals does not threaten population survival.

Identification of phasiRNAs and their drought- responsiveness in Populus trichocarpa

Phased, secondary, small interfering RNA (phasiRNA) perform essential biological functions in plants. However, limited information is available on the role of phasiRNA genes in Populus (poplar), especially during drought stress. In this study, we identified 20 PHAS loci generating 91 phasiRNA in the genome of the model forest tree Populus trichocarpa (P. trichocarpa; western balsam-poplar) using the control and drought libraries. Our analysis indicated that six PHAS (PtPHA14-20) initiated by two Populus-specific miRNAs (miR6445 and miR6427) were specific to Populus. In addition, a total of 47 phasiRNA were found to be drought responsive, and five of them were confirmed by RT-qPCR. The phase cleavage of three PHAS loci by miRNA, and degradation of nine transcript targets by phasiRNA were experimentally confirmed based on degradome data. The identification of these Populus phasiRNA will contribute to a better understanding of their function and regulation during drought stress.

Large-scale bioinformatic analysis of the regulation of the disease resistance NBS gene family by microRNAs in Poaceae

In the present study, we have screened 71, 713, 525, 119 and 241 mature miRNA variants from Hordeum vulgare, Oryza sativa, Brachypodium distachyon, Triticum aestivum, and Sorghum bicolor, respectively, and classified them with respect to their conservation status and expression levels. These Poaceae non-redundant miRNA species (1,669) were distributed over a total of 625 MIR families, among which only 54 were conserved across two or more plant species, confirming the relatively recent evolutionary differentiation of miRNAs in grasses. On the other hand, we have used 257 H. vulgare, 286T. aestivum, 119 B. distachyon, 269 O. sativa, and 139 S. bicolor NBS domains, which were either mined directly from the annotated proteomes, or predicted from whole genome sequence assemblies. The hybridization potential between miRNAs and their putative NBS genes targets was analyzed, revealing that at least 454 NBS genes from all five Poaceae were potentially regulated by 265 distinct miRNA species, most of them expressed in leaves and predominantly co-expressed in additional tissues. Based on gene ontology, we could assign these probable miRNA target genes to 16 functional groups, among which three conferring resistance to bacteria (Rpm1, Xa1 and Rps2), and 13 groups of resistance to fungi (Rpp8,13, Rp3, Tsn1, Lr10, Rps1-k-1, Pm3, Rpg5, and MLA1,6,10,12,13). The results of the present analysis provide a large-scale platform for a better understanding of biological control strategies of disease resistance genes in Poaceae, and will serve as an important starting point for enhancing crop disease resistance improvement by means of transgenic lines with artificial miRNAs.

Simple gene silencing using the trans-acting siRNA pathway

In plants, particular micro-RNAs (miRNAs) induce the production of a class of small interfering RNAs (siRNA) called trans-acting siRNA (ta-siRNA) that lead to gene silencing. A single miRNA target is sufficient for the production of ta-siRNAs, which target can be incorporated into a vector to induce the production of siRNAs, and ultimately gene silencing. The term miRNA-induced gene silencing (MIGS) has been used to describe such vector systems in Arabidopsis. Several ta-siRNA loci have been identified in soybean, but, prior to this work, few of the inducing miRNAs have been experimentally validated, much less used to silence genes. Nine ta-siRNA loci and their respective miRNA targets were identified, and the abundance of the inducing miRNAs varies dramatically in different tissues. The miRNA targets were experimentally verified by silencing a transgenic GFP gene and two endogenous genes in hairy roots and transgenic plants. Small RNAs were produced in patterns consistent with the utilization of the ta-siRNA pathway. A side-by-side experiment demonstrated that MIGS is as effective at inducing gene silencing as traditional hairpin vectors in soybean hairy roots. Soybean plants transformed with MIGS vectors produced siRNAs and silencing was observed in the T1 generation. These results complement previous reports in Arabidopsis by demonstrating that MIGS is an efficient way to produce siRNAs and induce gene silencing in other species, as shown with soybean. The miRNA targets identified here are simple to incorporate into silencing vectors and offer an effective and efficient alternative to other gene silencing strategies.

Identification of microRNAs and their targets in Finger millet by high throughput sequencing

MicroRNAs are short non-coding RNAs which play an important role in regulating gene expression by mRNA cleavage or by translational repression. The majority of identified miRNAs were evolutionarily conserved; however, others expressed in a species-specific manner. Finger millet is an important cereal crop; nonetheless, no practical information is available on microRNAs to date. In this study, we have identified 95 conserved microRNAs belonging to 39 families and 3 novel microRNAs by high throughput sequencing. For the identified conserved and novel miRNAs a total of 507 targets were predicted. 11 miRNAs were validated and tissue specificity was determined by stem loop RT-qPCR, Northern blot. GO analyses revealed targets of miRNA were involved in wide range of regulatory functions. This study implies large number of known and novel miRNAs found in Finger millet which may play important role in growth and development.

Genome-wide expression profiling of microRNAs in poplar upon infection with the foliar rust fungus Melampsora larici-populina

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that regulate the gene expression of target mRNAs involved in plant growth, development, and abiotic stress and pathogen responses. Previous studies have reported miRNAs in Populus that respond to abiotic stresses, such as cold, heat, drought, flooding, high salt and mechanical stress. However, little is known about the regulatory roles of these molecules in the Populus response to the stress of foliar rust fungal infection. Here, we identified the miRNA profiles of Populus after inoculation with Melampsora larici-populina using high-throughput sequencing and bioinformatics analysis. Quantitative real-time PCR (qRT-PCR) was used to validate the expression levels of 10 miRNAs.

RESULTS

A total of 90 known miRNAs belonging to 42 families and 378 novel miRNAs were identified from three small-RNA libraries of Populus szechuanica infected with M. larici-populina isolates Sb052 and Th053 and a control. Comparative analysis revealed that the expression of 38 known miRNAs and 92 novel miRNAs in P. szechuanica after infection with different rust fungus isolates showed significant differences, and more miRNAs were suppressed during rust infection. Among the differentially expressed miRNAs, 7 known and 20 novel miRNAs were relevant to the rust fungus infection, and according to KEGG (Kyoto Encyclopaedia of Genes and Genomes) pathway analysis, these miRNAs primarily regulate genes encoding disease-resistance proteins, serine/threonine protein kinases, transcription factors, and related proteins. QRT-PCR analysis indicated that most miRNAs were up-regulated in the Sb052 library and down-regulated in the Th053 library at 48 h post-inoculation (hpi).

CONCLUSIONS

These results demonstrate that the expression of miRNAs was altered in poplar under stress associated with M. larici-populina infection, and different temporal dynamics were observed in incompatible and compatible libraries. These findings suggest important roles for miRNA regulation in Populus upon infection with foliar rust fungus.

Characterization and comparative profiling of the small RNA transcriptomes in two phases of flowering in Cymbidium ensifolium

Background

Cymbidium ensifolium is one of the most important ornamental flowers in China, with an elegant shape, beautiful appearance, and a fragrant aroma. Its unique flower shape has long attracted scientists. MicroRNAs (miRNAs) are critical regulators in plant development and physiology, including floral development. However, to date, few studies have examined miRNAs in C. ensifolium.

Results

In this study, we employed Solexa technology to sequence four small RNA libraries from two flowering phases to identify miRNAs related to floral development. We identified 48 mature conserved miRNA and 71 precursors. These conserved miRNA belonged to 20 families. We also identified 45 novel miRNA which includes 21 putative novel miRNAs*, and 28 hairpin forming precursors. Two trans-acting small interfering RNAs (ta-siRNAs) were identified, one of which was homologous to TAS3a1. TAS3a1 belongs to the TAS3 family, which has been previously reported to target auxin response factors (ARF) and be involved in plant growth and floral development. Moreover, we built a C. ensifolium transctriptome database to identify genes targeted by miRNA, which resulted in 790 transcriptomic target unigenes. The target unigenes were annotated with information from the non-redundant (Nr), gene ontology database (GO), eukaryotic orthologous groups (KOGs) and Kyoto encyclopedia of genes and genomes (KEGG) database. The unigenes included MADS-box transcription factors targeted by miR156, miR172 and miR5179, and various hormone responding factors targeted by miR159. The MADS-box transcription factors are well known to determine the identity of flower organs and hormone responding factors involved in floral development. In expression analysis, three novel and four conserved miRNA were differentially expressed between two phases of flowering. The results were confirmed by RNA-seq and qRT-PCR. The differential expression of two miRNA, miR160 and miR396, targeted ARFs and growth regulating factor (GRF), respectively. However, most of these small RNA were clustered in the uncharacterized group, which suggests there may be many novel small non-coding RNAs yet to be discovered.

Conclusion

Our study provides a diverse set of miRNAs related to cymbidium floral development and serves as a useful resource for investigating miRNA-mediated regulatory mechanisms of floral development.

Electronic supplementary material

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

High-Throughput Sequencing Reveals Diverse Sets of Conserved, Nonconserved, and Species-Specific miRNAs in Jute

MicroRNAs play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. A group of plant microRNAs are evolutionarily conserved; however, others are expressed in a species-specific manner. Jute is an agroeconomically important fibre crop; nonetheless, no practical information is available for microRNAs in jute to date. In this study, Illumina sequencing revealed a total of 227 known microRNAs and 17 potential novel microRNA candidates in jute, of which 164 belong to 23 conserved families and the remaining 63 belong to 58 nonconserved families. Among a total of 81 identified microRNA families, 116 potential target genes were predicted for 39 families and 11 targets were predicted for 4 among the 17 identified novel microRNAs. For understanding better the functions of microRNAs, target genes were analyzed by Gene Ontology and their pathways illustrated by KEGG pathway analyses. The presence of microRNAs identified in jute was validated by stem-loop RT-PCR followed by end point PCR and qPCR for randomly selected 20 known and novel microRNAs. This study exhaustively identifies microRNAs and their target genes in jute which will ultimately pave the way for understanding their role in this crop and other crops.

Profiling of small RNAs involved in plant-pathogen interactions

Small RNA (sRNA)-mediated gene silencing is an important gene expression regulatory mechanism conserved in eukaryotes. Such sRNAs, first discovered in plants, are involved in diverse biological processes. In plants, sRNAs participate in many growth and developmental processes, such as embryo development, seed germination, flowering, hormone synthesis and distribution, and nutrient assimilation. However, the significance of sRNA in shaping the relationship between plants and their symbiotic microbes or pathogens has been underestimated. Recent progress in profiling sRNA, especially advances in next-generation sequencing technology, has revealed its extensive and complicated involvement in interactions between plants and viruses, bacteria, and fungi. In this review, we will summarize recent findings regarding sRNA in plant-pathogen interactions.

Novel and conserved heat-responsive microRNAs in wheat (Triticum aestivum L.)

MicroRNAs (miRNAs) are small endogenous RNAs of ~22 nucleotides that have been shown to play regulatory role by negatively affecting the expression of genes at the post-transcriptional level. Information of miRNAs on some important crops like soybean, Arabidopsis, and rice, etc. are available, but no study on heat-responsive novel miRNAs has yet been reported in wheat (Triticum aestivum L.). In the present investigation, a popular wheat cultivar HD2985 was used in small RNA library construction and Illumina HiSeq 2000 was used to perform high-throughput sequencing of the library after cluster generation; 110,896,604 and 87,743,861 reads were generated in the control (22 °C) and heat-treated (42 °C for 2 h) samples, respectively. Forty-four precursor and mature miRNAs were found in T. aestivum from miRBase v 19. The frequencies of the miRNA families varied from 2 (tae-miR1117) to 60,672 (tae-miR159b). We identify 1052 and 902 mature miRNA sequences in HD2985 control and HS-treated samples by mapping on reference draft genome of T. aestivum. Maximum identified miRNAs were located on IWGSC_CSS_3B_scaff (chromosome 3B). We could identify 53 and 46 mature miRNA in the control and HS samples and more than 516 target genes by mapping on the reference genome of Oryza sativa, Zea mays, and Sorghum bicolor. Using different pipelines and plant-specific criteria, 37 novel miRNAs were identified in the control and treated samples. Six novel miRNA were validated using qRT-PCR to be heat-responsive. A negative correlation was, however, observed between the expression of novel miRNAs and their targets. Target prediction and pathway analysis revealed their involvement in the heat stress tolerance. These novel miRNAs are new additions to miRNA database of wheat, and the regulatory network will be made use of in deciphering the mechanism of thermotolerance in wheat.

A Medicago truncatula rdr6 allele impairs transgene silencing and endogenous phased siRNA production but not development

RNA-dependent RNA polymerase 6 (RDR6) and suppressor of gene silencing 3 (SGS3) act together in post-transcriptional transgene silencing mediated by small interfering RNAs (siRNAs) and in biogenesis of various endogenous siRNAs including the tasiARFs, known regulators of auxin responses and plant development. Legumes, the third major crop family worldwide, has been widely improved through transgenic approaches. Here, we isolated rdr6 and sgs3 mutants in the model legume Medicago truncatula. Two sgs3 and one rdr6 alleles led to strong developmental defects and impaired biogenesis of tasiARFs. In contrast, the rdr6.1 homozygous plants produced sufficient amounts of tasiARFs to ensure proper development. High throughput sequencing of small RNAs from this specific mutant identified 354 potential MtRDR6 substrates, for which siRNA production was significantly reduced in the mutant. Among them, we found a large variety of novel phased loci corresponding to protein-encoding genes or transposable elements. Interestingly, measurement of GFP expression revealed that post-transcriptional transgene silencing was reduced in rdr6.1 roots. Hence, this novel mis-sense mutation, affecting a highly conserved amino acid residue in plant RDR6s, may be an interesting tool both to analyse endogenous pha-siRNA functions and to improve transgene expression, at least in legume species.

Identification of novel and conserved microRNAs in Coffea canephora and Coffea arabica

As microRNAs (miRNAs) are important regulators of many biological processes, a series of small RNAomes from plants have been produced in the last decade. However, miRNA data from several groups of plants are still lacking, including some economically important crops. Here microRNAs from Coffea canephora leaves were profiled and 58 unique sequences belonging to 33 families were found, including two novel microRNAs that have never been described before in plants. Some of the microRNA sequences were also identified in Coffea arabica that, together with C. canephora, correspond to the two major sources of coffee production in the world. The targets of almost all miRNAs were also predicted on coffee expressed sequences. This is the first report of novel miRNAs in the genus Coffea, and also the first in the plant order Gentianales. The data obtained establishes the basis for the understanding of the complex miRNA-target network on those two important crops.

Advances in identification and validation of plant microRNAs and their target genes

Developments in the field of molecular biology and genetics, such as microarray, gene transfer and discovery of small regulatory RNAs, have led to significant advances in plant biotechnology. Among the small RNAs, microRNAs (miRNAs) have elicited much interest as key post-transcriptional regulators in eukaryotic gene expression. Advances in genome and transcriptome sequencing of plants have facilitated the generation of a huge wealth of sequence information that can find much use in the discovery of novel miRNAs and their target genes. In this review, we present an overview of the developments in the strategies and methods used to identify and study miRNAs, their target genes and the mechanisms by which these miRNAs interact with their target genes since the discovery of the first miRNA. The approaches discussed include both reverse and forward genetics. We observed that despite the availability of advanced methods, certain limitations ranging from the cost of materials, equipment and personnel to the availability of genome sequences for many plant species present a number of challenges for the development and utilization of modern scientific methods for the elucidation and development of miRNAs in many important plant species.

Deep sequencing on a genome-wide scale reveals diverse stage-specific microRNAs in cambium during dormancy-release induced by chilling in poplar

BACKGROUND

Trees in temperate zones show periodicity by alternating active and dormant states to adapt to environmental conditions. Although phytohormones and transcriptional regulation were found to be involved in growth cessation and dormancy transition, little is known about the mechanisms of the dormancy-active growth transition, especially dormancy maintenance and release. Small RNAs are a group of short non-coding RNAs regulating gene expressions at the post-transcriptional level during plant development and the responses to environmental stress. No report on the expression profiling of small RNAs in the cambial meristem during the dormancy-active growth transition has been reported to date.

RESULTS

Three small RNA libraries from the cambium of poplar, representing endodormancy induced by short day conditions, ecodormancy induced by chilling and active growth induced by long day conditions, respectively, were generated and sequenced by Illumina high-throughput sequencing technology. This yielded 123 known microRNAs (miRNAs) with significant expression changes, which included developmental-, phytohormone- and stress-related miRNAs. Interestingly, miR156 and miR172 showed opposite expression patterns in the cambial dormancy-active growth transition. Additionally, miR160, which is involved in the auxin signaling pathway, was expressed specifically during endodormancy release by chilling. Furthermore, 275 novel miRNAs expressed in the cambial zone were identified, and 34 of them had high detection frequencies and unique expression patterns. Finally, the target genes of these novel miRNAs were predicted and some were validated experimentally by 5'RACE.

CONCLUSIONS

Our results provided a comprehensive analysis of small RNAs in the cambial meristem during dormancy-release at the genome-wide level and novel evidence of miRNAs involved in the regulation of this biological process.

Discovery of Novel Leaf Rust Responsive microRNAs in Wheat and Prediction of Their Target Genes

MicroRNAs are endogenous small noncoding RNAs which play critical roles in gene regulation. Few wheat (Triticum aestivum L.) miRNA sequences are available in miRBase repertoire and knowledge of their biological functions related to biotic stress is limited. We identified 52 miRNAs, belonging to 19 families, from next-generation transcriptome sequence data based on homology search. One wheat specific novel miRNA was identified but could not be ascribed or assigned to any known miRNA family. Differentially expressed 22 miRNAs were found between susceptible and resistant wheat near-isogenic lines inoculated with leaf rust pathogen Puccinia triticina and compared with mock inoculated controls. Most miRNAs were more upregulated in susceptible NIL compared to resistant NIL. We identified 1306 potential target genes for these 52 miRNAs with vital roles in response to stimuli, signaling, and diverse metabolic and cellular processes. Gene ontology analysis showed 66, 20, and 35 target genes to be categorized into biological process, molecular function, and cellular component, respectively. A miRNA-mediated regulatory network revealed relationships among the components of the targetome. The present study provides insight into potential miRNAs with probable roles in leaf rust pathogenesis and their target genes in wheat which establish a foundation for future studies.

Genome-wide analysis of salt-responsive and novel microRNAs in Populus euphratica by deep sequencing

BACKGROUND

Populus euphratica is a representative model woody plant species for studying resistance to abiotic stresses such as drought and salt. Salt stress is one of the most common environmental factors that affect plant growth and development. MicroRNAs (miRNAs) are small, noncoding RNAs that have important regulatory functions in plant growth, development, and response to abiotic stress.

RESULTS

To investigate the miRNAs involved in the salt-stress response, we constructed four small cDNA libraries from P. euphratica plantlets treated with or without salt (300 mM NaCl, 3 days) in either the root or leaf. Using high-throughput sequencing to identify miRNAs, we found 164 conserved miRNAs belonging to 44 families. Of these, 136 novel miRNAs were from the leaf, and 128 novel miRNAs were from the root. In response to salt stress, 95 miRNAs belonging to 46 conserved miRNAs families changed significantly, with 56 miRNAs upregulated and 39 miRNAs downregulated in the leaf. A comparison of the leaf and root tissues revealed 155 miRNAs belonging to 63 families with significantly altered expression, including 84 upregulated and 71 downregulated miRNAs. Furthermore, 479 target genes in the root and 541 targets of novel miRNAs in the leaf were predicted, and functional information was annotated using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases.

CONCLUSIONS

This study provides a novel visual field for understanding the regulatory roles of miRNAs in response to salt stress in Populus.

Evasion of short interfering RNA-directed antiviral silencing in Musa acuminata persistently infected with six distinct banana streak pararetroviruses

Vegetatively propagated crop plants often suffer from infections with persistent RNA and DNA viruses. Such viruses appear to evade the plant defenses that normally restrict viral replication and spread. The major antiviral defense mechanism is based on RNA silencing generating viral short interfering RNAs (siRNAs) that can potentially repress viral genes posttranscriptionally through RNA cleavage and transcriptionally through DNA cytosine methylation. Here we examined the RNA silencing machinery of banana plants persistently infected with six pararetroviruses after many years of vegetative propagation. Using deep sequencing, we reconstructed consensus master genomes of the viruses and characterized virus-derived and endogenous small RNAs. Consistent with the presence of endogenous siRNAs that can potentially establish and maintain DNA methylation, the banana genomic DNA was extensively methylated in both healthy and virus-infected plants. A novel class of abundant 20-nucleotide (nt) endogenous small RNAs with 5'-terminal guanosine was identified. In all virus-infected plants, 21- to 24-nt viral siRNAs accumulated at relatively high levels (up to 22% of the total small RNA population) and covered the entire circular viral DNA genomes in both orientations. The hotspots of 21-nt and 22-nt siRNAs occurred within open reading frame (ORF) I and II and the 5' portion of ORF III, while 24-nt siRNAs were more evenly distributed along the viral genome. Despite the presence of abundant viral siRNAs of different size classes, the viral DNA was largely free of cytosine methylation. Thus, the virus is able to evade siRNA-directed DNA methylation and thereby avoid transcriptional silencing. This evasion of silencing likely contributes to the persistence of pararetroviruses in banana plants.

IMPORTANCE

We report that DNA pararetroviruses in Musa acuminata banana plants are able to evade DNA cytosine methylation and transcriptional gene silencing, despite being targeted by the host silencing machinery generating abundant 21- to 24-nucleotide short interfering RNAs. At the same time, the banana genomic DNA is extensively methylated in both healthy and virus-infected plants. Our findings shed light on the siRNA-generating gene silencing machinery of banana and provide a possible explanation why episomal pararetroviruses can persist in plants whereas true retroviruses with an obligatory genome-integration step in their replication cycle do not exist in plants.

Genome-wide analysis of salt-responsive and novel microRNAs in Populus euphratica by deep sequencing

Background

Populus euphratica is a representative model woody plant species for studying resistance to abiotic stresses such as drought and salt. Salt stress is one of the most common environmental factors that affect plant growth and development. MicroRNAs (miRNAs) are small, noncoding RNAs that have important regulatory functions in plant growth, development, and response to abiotic stress.

Results

To investigate the miRNAs involved in the salt-stress response, we constructed four small cDNA libraries from P. euphratica plantlets treated with or without salt (300 mM NaCl, 3 days) in either the root or leaf. Using high-throughput sequencing to identify miRNAs, we found 164 conserved miRNAs belonging to 44 families. Of these, 136 novel miRNAs were from the leaf, and 128 novel miRNAs were from the root. In response to salt stress, 95 miRNAs belonging to 46 conserved miRNAs families changed significantly, with 56 miRNAs upregulated and 39 miRNAs downregulated in the leaf. A comparison of the leaf and root tissues revealed 155 miRNAs belonging to 63 families with significantly altered expression, including 84 upregulated and 71 downregulated miRNAs. Furthermore, 479 target genes in the root and 541 targets of novel miRNAs in the leaf were predicted, and functional information was annotated using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases.

Conclusions

This study provides a novel visual field for understanding the regulatory roles of miRNAs in response to salt stress in Populus.

Extending the sRNAome of apple by next-generation sequencing

The global importance of apple as a fruit crop necessitates investigations into molecular aspects of the processes that influence fruit quality and yield, including plant development, fruit ripening and disease resistance. In order to study and understand biological processes it is essential to recognise the range of molecules, which influence these processes. Small non-coding RNAs are regulatory agents involved in diverse plant activities, ranging from development to stress response. The occurrence of these molecules in apple leaves was studied by means of next-generation sequencing. 85 novel microRNA (miRNA) gene loci were predicted and characterized along with known miRNA loci. Both cis- and trans-natural antisense transcript pairs were identified. Although the trans-overlapping regions were enriched in small RNA (sRNA) production, cis-overlaps did not seem to agree. More than 150 phased regions were also identified, and for a small subset of these, potential miRNAs that could initiate phasing, were revealed. Repeat-associated siRNAs, which are generated from repetitive genomic regions such as transposons, were also analysed. For this group almost all available repeat sequences, associated with the apple genome and present in Repbase, were found to produce siRNAs. Results from this study extend our current knowledge on apple sRNAs and their precursors significantly. A rich molecular resource has been created and is available to the research community to serve as a baseline for future studies.

Cloning and characterization of miRNAs and their targets, including a novel miRNA-targeted NBS-LRR protein class gene in apple (Golden Delicious)

MicroRNA (miRNA) has emerged as an important regulator of gene expression in plants. 146 miRNAs were identified from apple (Malus domestica cv. Golden Delicious) by bioinformatic analysis and RNA library sequencing. From these, 135 were conserved and 11 were novel miRNAs. Target analysis predicted one of the novel miRNAs, Md-miRLn11 (Malus domestica microRNA Ln11), targeted an apple nucleotide-binding site (NBS)-leucine-rich repeat (LRR) class protein coding gene (Md-NBS). 5' RACE assay confirmed the ability of Md-miRLn11 to cleave Md-NBS at the 11-12-nt position. Analysis of the expression of Md-miRLn11 and Md-NBS during the optimum invasion period in 40 apple varieties showed that the expression of Md-NBS gene in resistant varieties is higher than in susceptible varieties, with an inverse pattern for Md-miRLn11. Seedlings from the resistant apple variety 'JiGuan' were used to carry out an Agrobacterium infiltration assay, and then inoculated with the apple leaf spot disease. The result showed a clear decline of disease resistance in JiGuan apples. In contrast, the susceptible variety 'FuJi' infiltrated with the Md-NBS gene showed a significant increase in disease resistance. Based on the above results, we propose that Md-miRLn11 regulates Md-NBS gene expression in particular under the condition of pathogen infection, and that the Md-miRLn11 targeting P-loop site may regulate many NBS-LRR protein class genes in woody plants.

Trans-acting small interfering RNA4: key to nutraceutical synthesis in grape development?

The facility and versatility of microRNAs (miRNAs) to evolve and change likely underlies how they have become dominant constituents of eukaryotic genomes. In this opinion article I propose that trans-acting small interfering RNA gene 4 (TAS4) evolution may be important for biosynthesis of polyphenolics, arbuscular symbiosis, and bacterial pathogen etiologies. Expression-based and phylogenetic evidence shows that TAS4 targets two novel grape (Vitis vinifera L.) MYB transcription factors (VvMYBA6, VvMYBA7) that spawn phased small interfering RNAs (siRNAs) which probably function in nutraceutical bioflavonoid biosynthesis and fruit development. Characterization of the molecular mechanisms of TAS4 control of plant development and integration into biotic and abiotic stress- and nutrient-signaling regulatory networks has applicability to molecular breeding and the development of strategies for engineering healthier foods.

Identification of microRNAs in Caragana intermedia by high-throughput sequencing and expression analysis of 12 microRNAs and their targets under salt stress

142 miRNAs were identified and 38 miRNA targets were predicted, 4 of which were validated, in C. intermedia . The expression of 12 miRNAs in salt-stressed leaves was assessed by qRT-PCR. MicroRNAs (miRNAs) are endogenous small RNAs that play important roles in various biological and metabolic processes in plants. Caragana intermedia is an important ecological and economic tree species prominent in the desert environment of west and northwest China. To date, no investigation into C. intermedia miRNAs has been reported. In this study, high-throughput sequencing of small RNAs and analysis of transcriptome data were performed to identify both conserved and novel miRNAs, and also their target mRNA genes in C. intermedia. Based on sequence similarity and hairpin structure prediction, 132 putative conserved miRNAs (12 of which were confirmed to form hairpin precursors) belonging to 31 known miRNA families were identified. Ten novel miRNAs (including the miRNA* sequences of three novel miRNAs) were also discovered. Furthermore, 36 potential target genes of 17 known miRNA families and 2 potential target genes of 1 novel miRNA were predicted; 4 of these were validated by 5' RACE. The expression of 12 miRNAs was validated in different tissues, and these and five target mRNAs were assessed by qRT-PCR after salt treatment. The expression levels of seven miRNAs (cin-miR157a, cin-miR159a, cin-miR165a, cin-miR167b, cin-miR172b, cin-miR390a and cin-miR396a) were upregulated, while cin-miR398a expression was downregulated after salt treatment. The targets of cin-miR157a, cin-miR165a, cin-miR172b and cin-miR396a were downregulated and showed an approximately negative correlation with their corresponding miRNAs under salt treatment. These results would help further understanding of miRNA regulation in response to abiotic stress in C. intermedia.

isomiRex: web-based identification of microRNAs, isomiR variations and differential expression using next-generation sequencing datasets

We present an open-access web platform isomiRex, to identify isomiRs and on the fly graphical visualization of the differentially expressed miRNAs in control as well as treated library. The open-access web-platform is not restricted only to NGS sequence dataset from animals and potentially analyzes a wider dataset for plants, animals and viral NGS dataset supporting miRBase (version 19 supporting 193 species). The platform can handle the bloated amount of the read counts and reports the annotated microRNAs from plant, animal and viral NGS datasets. isomiRex also provides an estimation of the the isomiRs, of miRNAs with higher copy number relative to their mature reference sequences indexed in miRBase (version 19 supporting 193 species). Visually enhanced graphs potentially display differentially expressed isomiRs, which will help the user to demonstrate and correlate the abundance of the isomiR as a signature event to the specific condition. An additional module for estimating the differential expression has been implemented allowing the users to postulate the differential expression across the user input samples. The developed web-platform can be accessed at http://bioinfo1.uni-plovdiv.bg/isomiRex/.

Deep sequencing discovery of novel and conserved microRNAs in strawberry (Fragaria×ananassa)

In plants, microRNAs (miRNAs) play a critical role in post-transcriptional gene regulation and have been shown to control many genes involved in various biological and metabolic processes. There have been extensive studies to discover miRNAs and analyze their functions in model plant species, such as Arabidopsis and rice. Deep sequencing technologies have facilitated identification of species-specific or lowly expressed as well as conserved or highly expressed miRNAs in plants. In this research, we used Solexa sequencing to discover new miRNAs in cultivated strawberry (Fragaria×ananassa). A total of 23,282 ,309 reads representing 22,500 ,402 distinct sequences were obtained from a short RNA library generated from small RNAs extracted from strawberry fruit tissues. On the basis of sequence similarity and hairpin structure prediction, we found that 156,639 reads representing 153 sequences have good matches to known miRNAs. We also identified 37 novel miRNA candidates. These sequences had not been previously described in other plant species. Potential target genes were predicted for the majority of and novel miRNAs. These results show that regulatory miRNAs exist in the agriculturally important cultivated strawberry and may play an important role in its growth, development and response to disease.

A significant fraction of 21-nucleotide small RNA originates from phased degradation of resistance genes in several perennial species

Small RNAs (sRNAs), including microRNA (miRNA) and short-interfering RNA (siRNA), are important in the regulation of diverse biological processes. Comparative studies of sRNAs from plants have mainly focused on miRNA, even though they constitute a mere fraction of the total sRNA diversity. In this study, we report results from an in-depth analysis of the sRNA population from the conifer spruce (Picea abies) and compared the results with those of a range of plant species. The vast majority of sRNA sequences in spruce can be assigned to 21-nucleotide-long siRNA sequences, of which a large fraction originate from the degradation of transcribed sequences related to nucleotide-binding site-leucine-rich repeat-type resistance genes. Over 90% of all genes predicted to contain either a Toll/interleukin-1 receptor or nucleotide-binding site domain showed evidence of siRNA degradation. The data further suggest that this phased degradation of resistance-related genes is initiated from miRNA-guided cleavage, often by an abundant 22-nucleotide miRNA. Comparative analysis over a range of plant species revealed a huge variation in the abundance of this phenomenon. The process seemed to be virtually absent in several species, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and nonvascular plants, while particularly high frequencies were observed in spruce, grape (Vitis vinifera), and poplar (Populus trichocarpa). This divergent pattern might reflect a mechanism to limit runaway transcription of these genes in species with rapidly expanding nucleotide-binding site-leucine-rich repeat gene families. Alternatively, it might reflect variation in a counter-counter defense mechanism between plant species.

Sexual dimorphism floral microRNA profiling and target gene expression in andromonoecious poplar (Populus tomentosa)

Although the molecular basis of poplar sex-specific flower development remains largely unknown, increasing evidence indicates an essential role for microRNAs (miRNAs). The specific miRNA types and precise miRNA expression patterns in dioecious plant flower development remain unclear. Here, we used andromonoecious poplar, an exceptional model system, to eliminate the confounding effects of genetic background of dioecious plants. This system, combined with high-throughput sequencing and computational analysis, allowed us to characterize sex-specific miRNAomes from female and male flowers. Comparative miRNAome analysis combined with quantitative real-time PCR revealed the expression patterns of 27 miRNAs in poplar flower and showed that the targets of these miRNAs are involved in flower organogenesis, Ca(2+) transport, phytohormone synthesis and metabolism, and DNA methylation. This paper describes a complex regulatory network consisting of these miRNAs expressed in sex-specific flower development in a dioecious plant. The conserved and novel miRNA locations were annotated in the Populus trichocarpa genome. Among these, miRNA Pto-F70 and 4 targets are located in the sex-determination regions of chromosome XIX. Furthermore, two novel miRNAs, Pto-F47 and Pto-F68, were shown for the first time to be regulatory factors in phytohormone interactions. To our knowledge, this report is the first systematic investigation of sex-specific flower-related miRNAs and their targets in poplar, and it deepens our understanding of the important regulatory functions of miRNAs in female and male flower development in this dioecious plant.

MicroRNA superfamilies descended from miR390 and their roles in secondary small interfering RNA Biogenesis in Eudicots

Trans-acting small interfering RNAs (tasiRNAs) are a major class of small RNAs performing essential biological functions in plants. The first reported tasiRNA pathway, that of miR173-TAS1/2, produces tasiRNAs regulating a set of pentatricopeptide repeat (PPR) genes and has been characterized only in Arabidopsis thaliana to date. Here, we demonstrate that the microRNA (miRNA)-trans-acting small interfering RNA gene (TAS)-pentatricopeptide repeat-containing gene (PPR)-small interfering RNA pathway is a highly dynamic and widespread feature of eudicots. Nine eudicot plants, representing six different plant families, have evolved similar tasiRNA pathways to initiate phased small interfering RNA (phasiRNA) production from PPR genes. The PPR phasiRNA production is triggered by different 22-nucleotide miRNAs, including miR7122, miR1509, and fve-PPRtri1/2, and through distinct mechanistic strategies exploiting miRNA direct targeting or indirect targeting through TAS-like genes (TASL), one-hit or two-hit, or even two layers of tasiRNA-TASL interactions. Intriguingly, although those miRNA triggers display high sequence divergence caused by the occurrence of frequent point mutations and splicing shifts, their corresponding MIRNA genes show pronounced identity to the Arabidopsis MIR173, implying a common origin of this group of miRNAs (super-miR7122). Further analyses reveal that super-miR7122 may have evolved from a newly defined miR4376 superfamily, which probably originated from the widely conserved miR390. The elucidation of this evolutionary path expands our understanding of the course of miRNA evolution, especially for relatively conserved miRNA families.

Identification of drought-responsive and novel Populus trichocarpa microRNAs by high-throughput sequencing and their targets using degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are endogenous small RNAs (sRNAs) with a wide range of regulatory functions in plant development and stress responses. Although miRNAs associated with plant drought stress tolerance have been studied, the use of high-throughput sequencing can provide a much deeper understanding of miRNAs. Drought is a common stress that limits the growth of plants. To obtain more insight into the role of miRNAs in drought stress, Illumina sequencing of Populus trichocarpa sRNAs was implemented.

RESULTS

Two sRNA libraries were constructed by sequencing data of control and drought stress treatments of poplar leaves. In total, 207 P. trichocarpa conserved miRNAs were detected from the two sRNA libraries. In addition, 274 potential candidate miRNAs were found; among them, 65 candidates with star sequences were chosen as novel miRNAs. The expression of nine conserved miRNA and three novel miRNAs showed notable changes in response to drought stress. This was also confirmed by quantitative real time polymerase chain reaction experiments. To confirm the targets of miRNAs experimentally, two degradome libraries from the two treatments were constructed. According to degradome sequencing results, 53 and 19 genes were identified as targets of conserved and new miRNAs, respectively. Functional analysis of these miRNA targets indicated that they are involved in important activities such as the regulation of transcription factors, the stress response, and lipid metabolism.

CONCLUSIONS

We discovered five upregulated miRNAs and seven downregulated miRNAs in response to drought stress. A total of 72 related target genes were detected by degradome sequencing. These findings reveal important information about the regulation mechanism of miRNAs in P. trichocarpa and promote the understanding of miRNA functions during the drought response.

Emerging role for RNA-based regulation in plant immunity

Infection by phytopathogenic bacteria triggers massive changes in plant gene expression, which are thought to be mostly a result of transcriptional reprogramming. However, evidence is accumulating that plants additionally use post-transcriptional regulation of immune-responsive mRNAs as a strategic weapon to shape the defense-related transcriptome. Cellular RNA-binding proteins regulate RNA stability, splicing or mRNA export of immune-response transcripts. In particular, mutants defective in alternative splicing of resistance genes exhibit compromised disease resistance. Furthermore, detection of bacterial pathogens induces the differential expression of small non-coding RNAs including microRNAs that impact the host defense transcriptome. Phytopathogenic bacteria in turn have evolved effector proteins to inhibit biogenesis and/or activity of cellular microRNAs. Whereas RNA silencing has long been known as an antiviral defense response, recent findings also reveal a major role of this process in antibacterial defense. Here we review the function of RNA-binding proteins and small RNA-directed post-transcriptional regulation in antibacterial defense. We mainly focus on studies that used the model system Arabidopsis thaliana and also discuss selected examples from other plants.

Methylome of DNase I sensitive chromatin in Populus trichocarpa shoot apical meristematic cells: a simplified approach revealing characteristics of gene-body DNA methylation in open chromatin state

DNA methylation is involved in the control of plant development and adaptation to the environment through modifications of chromatin compaction and gene expression. In poplar (Populus trichocarpa), a perennial plant, variations in DNA methylation have been reported between genotypes and tissues or in response to drought. Nevertheless, the relationships between gene-body DNA methylation, gene expression and chromatin compaction still need clarification. Here, DNA methylation was mapped in the noncondensed chromatin fraction from P. trichocarpa shoot apical meristematic cells, the center of plant morphogenesis, where DNA methylation variations could influence the developmental trajectory. DNase I was used to isolate the noncondensed chromatin fraction. Methylated sequences were immunoprecipitated, sequenced using Illumina/Solexa technology and mapped on the v2.0 poplar genome. Bisulfite sequencing of candidate sequences was used to confirm mapping data and to assess cytosine contexts and methylation levels. While the methylated DNase I hypersensitive site fraction covered 1.9% of the poplar genome, it contained sequences corresponding to 74% of poplar gene models, mostly exons. The level and cytosine context of gene-body DNA methylation varied with the structural characteristics of the genes. Taken together, our data show that DNA methylation is widespread and variable among genes in open chromatin of meristematic cells, in agreement with a role in their developmental trajectory.

Evidence for small RNAs homologous to effector-encoding genes and transposable elements in the oomycete Phytophthora infestans

Phytophthora infestans is the oomycete pathogen responsible for the devastating late blight disease on potato and tomato. There is presently an intense research focus on the role(s) of effectors in promoting late blight disease development. However, little is known about how they are regulated, or how diversity in their expression may be generated among different isolates. Here we present data from investigation of RNA silencing processes, characterized by non-coding small RNA molecules (sRNA) of 19-40 nt. From deep sequencing of sRNAs we have identified sRNAs matching numerous RxLR and Crinkler (CRN) effector protein genes in two isolates differing in pathogenicity. Effector gene-derived sRNAs were present in both isolates, but exhibited marked differences in abundance, especially for CRN effectors. Small RNAs in P. infestans grouped into three clear size classes of 21, 25/26 and 32 nt. Small RNAs from all size classes mapped to RxLR effector genes, but notably 21 nt sRNAs were the predominant size class mapping to CRN effector genes. Some effector genes, such as PiAvr3a, to which sRNAs were found, also exhibited differences in transcript accumulation between the two isolates. The P. infestans genome is rich in transposable elements, and the majority of sRNAs of all size classes mapped to these sequences, predominantly to long terminal repeat (LTR) retrotransposons. RNA silencing of Dicer and Argonaute genes provided evidence that generation of 21 nt sRNAs is Dicer-dependent, while accumulation of longer sRNAs was impacted by silencing of Argonaute genes. Additionally, we identified six microRNA (miRNA) candidates from our sequencing data, their precursor sequences from the genome sequence, and target mRNAs. These miRNA candidates have features characteristic of both plant and metazoan miRNAs.

Identification and characterization of salt-responsive microRNAs in Populus tomentosa by high-throughput sequencing

Salt is one of the main environmental factors limiting plant growth and a better understanding of mechanisms of salt stress would aid efforts to bolster plant salt tolerance. MicroRNAs are well known for their important regulatory roles in response to abiotic stress in plants. In this study, high-throughput sequencing was employed to identify miRNAs in Populus tomentosa plantlets treated or not with salt (200 mM for 10 h). We found 141 conserved miRNAs belonging to 31 families, 29 non-conserved but previously-known miRNAs belonging to 26 families, and 17 novel miRNAs. Under salt stress, 19 miRNAs belonging to seven conserved miRNA families were significantly downregulated, and two miRNAs belonging to two conserved miRNA families were upregulated. Of seven non-conserved miRNAs with significantly altered expression, five were downregulated and two were upregulated. Furthermore, eight miRNAs were validated by qRT-PCR and their dynamic differential expressions were analyzed. In addition, 269 target genes of identified miRNAs were predicted and categorized by function. These results provide new insights into salt-responsive miRNAs in Populus.

smRNAome profiling to identify conserved and novel microRNAs in Stevia rebaudiana Bertoni

BACKGROUND

MicroRNAs (miRNAs) constitute a family of small RNA (sRNA) population that regulates the gene expression and plays an important role in plant development, metabolism, signal transduction and stress response. Extensive studies on miRNAs have been performed in different plants such as Arabidopsis thaliana, Oryza sativa etc. and volume of the miRNA database, mirBASE, has been increasing on day to day basis. Stevia rebaudiana Bertoni is an important perennial herb which accumulates high concentrations of diterpene steviol glycosides which contributes to its high indexed sweetening property with no calorific value. Several studies have been carried out for understanding molecular mechanism involved in biosynthesis of these glycosides, however, information about miRNAs has been lacking in S. rebaudiana. Deep sequencing of small RNAs combined with transcriptomic data is a powerful tool for identifying conserved and novel miRNAs irrespective of availability of genome sequence data.

RESULTS

To identify miRNAs in S. rebaudiana, sRNA library was constructed and sequenced using Illumina genome analyzer II. A total of 30,472,534 reads representing 2,509,190 distinct sequences were obtained from sRNA library. Based on sequence similarity, we identified 100 miRNAs belonging to 34 highly conserved families. Also, we identified 12 novel miRNAs whose precursors were potentially generated from stevia EST and nucleotide sequences. All novel sequences have not been earlier described in other plant species. Putative target genes were predicted for most conserved and novel miRNAs. The predicted targets are mainly mRNA encoding enzymes regulating essential plant metabolic and signaling pathways.

CONCLUSIONS

This study led to the identification of 34 highly conserved miRNA families and 12 novel potential miRNAs indicating that specific miRNAs exist in stevia species. Our results provided information on stevia miRNAs and their targets building a foundation for future studies to understand their roles in key stevia traits.

Global alteration of microRNAs and transposon-derived small RNAs in cotton (Gossypium hirsutum) during Cotton leafroll dwarf polerovirus (CLRDV) infection

Small RNAs (sRNAs) are a class of non-coding RNAs ranging from 20- to 40-nucleotides (nts) that are present in most eukaryotic organisms. In plants, sRNAs are involved in the regulation of development, the maintenance of genome stability and the antiviral response. Viruses, however, can interfere with and exploit the silencing-based regulatory networks, causing the deregulation of sRNAs, including small interfering RNAs (siRNAs) and microRNAs (miRNAs). To understand the impact of viral infection on the plant sRNA pathway, we deep sequenced the sRNAs in cotton leaves infected with Cotton leafroll dwarf virus (CLRDV), which is a member of the economically important virus family Luteoviridae. A total of 60 putative conserved cotton miRNAs were identified, including 19 new miRNA families that had not been previously described in cotton. Some of these miRNAs were clearly misregulated during viral infection, and their possible role in symptom development and disease progression is discussed. Furthermore, we found that the 24-nt heterochromatin-associated siRNAs were quantitatively and qualitatively altered in the infected plant, leading to the reactivation of at least one cotton transposable element. This is the first study to explore the global alterations of sRNAs in virus-infected cotton plants. Our results indicate that some CLRDV-induced symptoms may be correlated with the deregulation of miRNA and/or epigenetic networks.

MicroR828 regulates lignin and H2O2 accumulation in sweet potato on wounding

MicroRNAs (miRNAs) are small noncoding RNAs which post-transcriptionally regulate gene expression by directing mRNA cleavage or translational inhibition. miRNAs play multiple roles in the growth, development and stress responses in plants. However, little is known of the wounding-responsive miRNAs and their regulation. Here, we investigated the expression patterns of microR828 (miR828) on wounding in sweet potato (Ipomoea batatas cv Tainung 57). The expression of miR828 was only detected in leaves, and was induced by wounding rather than by ethylene, hydrogen peroxide (H2O2), methyl jasmonate or nitric oxide (NO). Moreover, cyclic guanosine monophosphate (cGMP) was necessary for miR828 accumulation in leaves on wounding. Two miR828 target candidates, named IbMYB and IbTLD, were obtained by cDNA cloning, and their mRNA cleavage caused by miR828 was confirmed by cleavage site mapping, agro-infiltration and transgenics studies. The reduction in IbMYB and IbTLD expression coincided with the induction of miR828, demonstrating that IbMYB and IbTLD might be miR828 targets. Furthermore, transgenic sweet potato overexpressing miR828 precursor affected lignin and H2O2 contents. These results showed that cGMP could regulate wounding-responsive miR828, which repressed the expression of IbMYB and IbTLD. Subsequently, lignin and H2O2 were accumulated to participate in defense mechanisms.

High-throughput sequencing and characterization of the small RNA transcriptome reveal features of novel and conserved microRNAs in Panax ginseng

microRNAs (miRNAs) play vital regulatory roles in many organisms through direct cleavage of transcripts, translational repression, or chromatin modification. Identification of miRNAs has been carried out in various plant species. However, no information is available for miRNAs from Panax ginseng, an economically significant medicinal plant species. Using the next generation high-throughput sequencing technology, we obtained 13,326,328 small RNA reads from the roots, stems, leaves and flowers of P. ginseng. Analysis of these small RNAs revealed the existence of a large, diverse and highly complicated small RNA population in P. ginseng. We identified 73 conserved miRNAs, which could be grouped into 33 families, and 28 non-conserved ones belonging to 9 families. Characterization of P. ginseng miRNA precursors revealed many features, such as production of two miRNAs from distinct regions of a precursor, clusters of two precursors in a transcript, and generation of miRNAs from both sense and antisense transcripts. It suggests the complexity of miRNA production in P. gingseng. Using a computational approach, we predicted for the conserved and non-conserved miRNA families 99 and 31 target genes, respectively, of which eight were experimentally validated. Among all predicted targets, only about 20% are conserved among various plant species, whereas the others appear to be non-conserved, indicating the diversity of miRNA functions. Consistently, many miRNAs exhibited tissue-specific expression patterns. Moreover, we identified five dehydration- and ten heat-responsive miRNAs and found the existence of a crosstalk among some of the stress-responsive miRNAs. Our results provide the first clue to the elucidation of miRNA functions in P. ginseng.