An atlas of soybean small RNAs identifies phased siRNAs from hundreds of coding genes

Genome-wide comparative analysis in Solanaceous species reveals evolution of microRNAs targeting defense genes in Capsicum spp

MicroRNAs (miRNAs) play roles in various biological processes in plants including growth, development, and disease resistance. Previous studies revealed that some plant miRNAs produce secondary small interfering RNAs (siRNAs) such as phased, secondary siRNAs (phasiRNAs), and they regulate a cascade of gene expression. We performed a genome-wide comparative analysis of miRNAs in Solanaceous species (pepper, tomato, and potato), from an evolutionary perspective. Microsynteny of miRNAs was analysed based on the genomic loci and their flanking genes and most of the well-conserved miRNA genes maintained microsynteny in Solanaceae. We identified target genes of the miRNAs via degradome analysis and found that several miRNAs target many genes encoding nucleotide-binding leucine-rich repeat (NLR) or receptor-like proteins (RLPs), which are known to be major players in defense responses. In addition, disease-resistance-associated miRNAs trigger phasiRNA production in pepper, indicating amplification of the regulation of disease-resistance gene families. Among these, miR-n033a-3p, whose target NLRs have been duplicated in pepper, targets more NLRs belonging to specific subgroup in pepper than those in potato. miRNAs targeting resistance genes might have evolved to regulate numerous targets in Solanaceae, following expansion of target resistance genes. This study provides an insight into evolutionary relationship between miRNAs and their target defense genes in plants.

Genome-wide identification and characterization of novel microRNAs in seed development of soybean

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in eukaryotes. However, the information about miRNAs population and their regulatory functions involving in soybean seed development remains incomplete. Base on the Dicer-like1-mediated cleavage signals during miRNA processing could be employed for novel miRNA discovery, a genome-wide search for miRNA candidates involved in seed development was carried out. As a result, 17 novel miRNAs, 14 isoforms of miRNA (isomiRs) and 31 previously validated miRNAs were discovered. These novel miRNAs and isomiRs represented tissue-specific expression and the isomiRs showed significantly higher abundance than that of their miRNA counterparts in different tissues. After target prediction and degradome sequencing data-based validation, 13 novel miRNA-target pairs were further identified. Besides, five targets of 22-nt iso-gma-miR393h were found to be triggered to produce secondary trans-acting siRNA (ta-siRNAs). Summarily, our results could expand the repertoire of miRNAs with potentially important functions in soybean.

Leaf rust (Puccinia triticina) mediated RNAi in wheat (Triticum aestivum L.) prompting host susceptibility

Significance of microRNAs in regulating gene expression in higher eukaryotes as well as in pathogens like fungi to suppress host defense is a well-established phenomenon. The present study focuses on leaf rust fungi Puccinia triticina (Pathotype 77-5) mediated RNAi to make wheat (Triticum aestivum L.) more susceptible. To reach such conclusions, we first confirmed the presence of argonaute (AGO) and dicer-like protein (DCL) family sequences in Puccinia. Bioinformatic tools were applied to retrieve the sequences from Puccinia genome followed by cloning and sequencing from P. triticina pathotype 77-5 cDNA to obtain the specific sequences. Their homologs were searched in other 14 Puccinia races to relate them with pathogenesis. Further, precursor sequences for three miRNA-like RNA molecules (milRs) were cloned from P. triticina cDNA. Their target genes like MAP kinase were successfully predicted and validated through degradome mapping and qRT-PCR. Gradual increase in milR2 (milR and milR*) expression over progressive time point of infection and positive expression for all the milRs within 77-5 urediniospores confirmed a complete host- independent RNAi activity by P. triticina.

Enhanced resistance to bacterial and oomycete pathogens by short tandem target mimic RNAs in tomato

Plants, like animals, have complex disease resistance systems in which receptors interact directly or indirectly with effectors of disease produced by pests and pathogens. To minimize the fitness cost of these systems to the plant, there are miRNAs that target the mRNAs of a family of receptor proteins required for disease resistance. Target site mimics of these miRNAs confer enhanced quantitative resistance in tomato against an oomycete and a bacterium. These findings are consistent with a role of the receptor proteins in quantitative disease resistance and show how blocking these miRNAs could be a useful approach in crop protection.

Nucleotide binding site leucine-rich repeat (NLR) proteins of the plant innate immune system are negatively regulated by the miR482/2118 family miRNAs that are in a distinct 22-nt class of miRNAs with a double mode of action. First, they cleave the target RNA, as with the canonical 21-nt miRNAs, and second, they trigger secondary siRNA production using the target RNA as a template. Here, we address the extent to which the miR482/2118 family affects expression of NLR mRNAs and disease resistance. We show that structural differences of miR482/2118 family members in tomato (Solanum lycopersicum) are functionally significant. The predicted target of the miR482 subfamily is a conserved motif in multiple NLR mRNAs, whereas for miR2118b, it is a noncoding RNA target formed by rearrangement of several different NLR genes. From RNA sequencing and degradome data in lines expressing short tandem target mimic (STTM) RNAs of miR482/2118, we confirm the different targets of these miRNAs. The effect on NLR mRNA accumulation is slight, but nevertheless, the tomato STTM lines display enhanced resistance to infection with the oomycete and bacterial pathogens. These data implicate an RNA cascade of miRNAs and secondary siRNAs in the regulation of NLR RNAs and show that the encoded NLR proteins have a role in quantitative disease resistance in addition to dominant gene resistance that has been well characterized elsewhere. We also illustrate the use of STTM RNA in a biotechnological approach for enhancing quantitative disease resistance in highly bred cultivars.

Evolution of Disease Defense Genes and Their Regulators in Plants

Biotic stresses do damage to the growth and development of plants, and yield losses for some crops. Confronted with microbial infections, plants have evolved multiple defense mechanisms, which play important roles in the never-ending molecular arms race of plant–pathogen interactions. The complicated defense systems include pathogen-associated molecular patterns (PAMP) triggered immunity (PTI), effector triggered immunity (ETI), and the exosome-mediated cross-kingdom RNA interference (CKRI) system. Furthermore, plants have evolved a classical regulation system mediated by miRNAs to regulate these defense genes. Most of the genes/small RNAs or their regulators that involve in the defense pathways can have very rapid evolutionary rates in the longitudinal and horizontal co-evolution with pathogens. According to these internal defense mechanisms, some strategies such as molecular switch for the disease resistance genes, host-induced gene silencing (HIGS), and the new generation of RNA-based fungicides, have been developed to control multiple plant diseases. These broadly applicable new strategies by transgene or spraying ds/sRNA may lead to reduced application of pesticides and improved crop yield.

Banana sRNAome and degradome identify microRNAs functioning in differential responses to temperature stress

Background

Temperature stress is a major environmental factor affecting not only plant growth and development, but also fruit postharvest life and quality. MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in various biological processes. Harvested banana fruit can exhibit distinct symptoms in response to different temperature stresses, but the underlying miRNA-mediated regulatory mechanisms remained unknown.

Results

Here, we profiled temperature-responsive miRNAs in banana, using deep sequencing and computational and molecular analyses. In total 113 known miRNAs and 26 novel banana-specific miRNAs were identified. Of these miRNAs, 42 miRNAs were expressed differentially under cold and heat stresses. Degradome sequencing identified 60 target genes regulated by known miRNAs and half of these targets were regulated by 15 temperature-responsive miRNAs. The correlative expression patterns between several miRNAs and their target genes were further validated via qRT-PCR. Our data showed that miR535 and miR156 families may derive from a common ancestor during evolution and jointly play a role in fine-tuning SPL gene expression in banana. We also identified the miRNA-triggered phased secondary siRNAs in banana and found miR393-TIR1/AFB phasiRNA production displaying cold stress-specific enrichment.

Conclusions

Our results provide a foundation for understanding the miRNA-dependent temperature stress response in banana. The characterized correlations between miRNAs and their response to temperature stress could serve as markers in the breeding programs or tools for improving temperature tolerance of banana.

Electronic supplementary material

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

Molecular mechanisms that limit the costs of NLR-mediated resistance in plants

Crop diseases cause significant yield losses, and the use of resistant cultivars can effectively mitigate these losses and control many plant diseases. Most plant resistance (R) genes encode immune receptors composed of nucleotide-binding and leucine-rich repeat (NLR) domains. These proteins mediate the specific recognition of pathogen avirulence effectors to induce defence responses. However, NLR-triggered immunity can be associated with a reduction in growth and yield, so-called 'fitness costs'. Recent data have shown that plants use an elaborate interplay of different mechanisms to control NLR gene transcript levels, as well as NLR protein abundance and activity, to avoid the associated cost of resistance in the absence of a pathogen. In this review, we discuss the different levels of NLR regulation (transcriptional, post-transcriptional and at the protein level). We address the apparent need for plants to maintain diverse modes of regulation. A recent model suggesting an equilibrium 'ON/OFF state' of NLR proteins, in the absence of a pathogen, provides the context for our discussion.

Cis-directed cleavage and nonstoichiometric abundances of 21-nucleotide reproductive phased small interfering RNAs in grasses

Post-transcriptional gene silencing in plants results from independent activities of diverse small RNA types. In anthers of grasses, hundreds of loci yield noncoding RNAs that are processed into 21- and 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs); these are triggered by miR2118 and miR2275. We characterized these 'reproductive phasiRNAs' from rice (Oryza sativa) panicles and anthers across seven developmental stages. Our computational analysis identified characteristics of the 21-nt reproductive phasiRNAs that impact their biogenesis, stability, and potential functions. We demonstrate that 21-nt reproductive phasiRNAs can function in cis to target their own precursors. We observed evidence of this cis regulatory activity in both rice and maize (Zea mays). We validated this activity with evidence of cleavage and a resulting shift in the pattern of phasiRNA production. We characterize biases in phasiRNA biogenesis, demonstrating that the Pol II-derived 'top' strand phasiRNAs are consistently higher in abundance than the bottom strand. The first phasiRNA from each precursor overlaps the miR2118 target site, and this impacts phasiRNA accumulation or stability, evident in the weak accumulation of this phasiRNA position. Additional influences on this first phasiRNA duplex include the sequence composition and length, and we show that these factors impact Argonaute loading.

Cis-directed cleavage and nonstoichiometric abundances of 21-nucleotide reproductive phased small interfering RNAs in grasses

Post-transcriptional gene silencing in plants results from independent activities of diverse small RNA types. In anthers of grasses, hundreds of loci yield noncoding RNAs that are processed into 21- and 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs); these are triggered by miR2118 and miR2275. We characterized these 'reproductive phasiRNAs' from rice (Oryza sativa) panicles and anthers across seven developmental stages. Our computational analysis identified characteristics of the 21-nt reproductive phasiRNAs that impact their biogenesis, stability, and potential functions. We demonstrate that 21-nt reproductive phasiRNAs can function in cis to target their own precursors. We observed evidence of this cis regulatory activity in both rice and maize (Zea mays). We validated this activity with evidence of cleavage and a resulting shift in the pattern of phasiRNA production. We characterize biases in phasiRNA biogenesis, demonstrating that the Pol II-derived 'top' strand phasiRNAs are consistently higher in abundance than the bottom strand. The first phasiRNA from each precursor overlaps the miR2118 target site, and this impacts phasiRNA accumulation or stability, evident in the weak accumulation of this phasiRNA position. Additional influences on this first phasiRNA duplex include the sequence composition and length, and we show that these factors impact Argonaute loading.

A Short Indel-Lacking-Resistance Gene Triggers Silencing of the Photosynthetic Machinery Components Through TYLCSV-Associated Endogenous siRNAs in Tomato

Plant viruses modify gene expression in infected tissues by altering the micro (mi)RNA-mediated regulation of genes. Among conserved miRNA targets there are transcripts coding for transcription factors, RNA silencing core, and disease-resistance proteins. Paralogs in these gene families are widely present in plant genomes and are known to respond differently to miRNA-mediated regulation during plant virus infections. Using genome-wide approaches applied to Solanum lycopersicum infected by a nuclear-replicating virus, we highlighted miRNA-mediated cleavage events that could not be revealed in virus-free systems. Among them we confirmed miR6024 targeting and cleavage of RX-coiled-coil (RX-CC), nucleotide binding site (NBS), leucine-rich (LRR) mRNA. Cleavage of paralogs was associated with short indels close to the target sites, indicating a general functional significance of indels in fine-tuning gene expression in plant-virus interaction. miR6024-mediated cleavage, uniquely in virus-infected tissues, triggers the production of several 21-22 nt secondary siRNAs. These secondary siRNAs, rather than being involved in the cascade regulation of other NBS-LRR paralogs, explained cleavages of several mRNAs annotated as defence-related proteins and components of the photosynthetic machinery. Outputs of these data explain part of the phenotype plasticity in plants, including the appearance of yellowing symptoms in the viral pathosystem.

Plant 24-nt reproductive phasiRNAs from intramolecular duplex mRNAs in diverse monocots

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

A Single-Nucleotide Polymorphism in the Promoter of a Hairpin RNA Contributes to Alternaria alternata Leaf Spot Resistance in Apple (Malus × domestica)

Apple leaf spot caused by the Alternaria alternata f. sp mali (ALT1) fungus is one of the most devastating diseases of apple (Malus × domestica). We identified a hairpin RNA (hpRNA) named MdhpRNA277 that produces small RNAs and is induced by ALT1 infection in 'Golden Delicious' apple. MdhpRNA277 produces mdm-siR277-1 and mdm-siR277-2, which target five resistance (R) genes that are expressed at high levels in resistant apple variety 'Hanfu' and at low levels in susceptible variety 'Golden Delicious' following ALT1 infection. MdhpRNA277 was strongly induced in 'Golden Delicious' but not 'Hanfu' following ALT1 inoculation. MdhpRNA277 promoter activity was much stronger in inoculated 'Golden Delicious' versus 'Hanfu'. We identified a single-nucleotide polymorphism (SNP) in the MdhpRNA277 promoter region between 'Golden Delicious' (pMdhpRNA277-GD) and 'Hanfu' (pMdhpRNA277-HF). The transcription factor MdWHy binds to pMdhpRNA277-GD, but not to pMdhpRNA277-HF Transgenic 'GL-3' apple expressing pMdhpRNA277-GD:MdhpRNA277 was more susceptible to ALT1 infection than plants expressing pMdhpRNA277-HF:MdhpRNA277 due to induced mdm-siR277 accumulation and reduced expression of the five target R genes. We confirmed that the SNP in pMdhpRNA277 is associated with A. alternata leaf spot resistance by crossing. This SNP could be used as a marker to distinguish between apple varieties that are resistant or susceptible to A. alternata leaf spot.

Biogenesis of a 22-nt microRNA in Phaseoleae species by precursor-programmed uridylation

Phased, secondary siRNAs (phasiRNAs) represent a class of small RNAs in plants generated via distinct biogenesis pathways, predominantly dependent on the activity of 22-nt miRNAs. Most 22-nt miRNAs are processed by DCL1 from miRNA precursors containing an asymmetric bulge, yielding a 22/21-nt miRNA/miRNA* duplex. Here we show that miR1510, a soybean miRNA capable of triggering phasiRNA production from numerous nucleotide-binding leucine-rich repeat (NB-LRRs), previously described as 21 nt in its mature form, primarily accumulates as a 22-nt isoform via monouridylation. We demonstrate that, in Arabidopsis, this uridylation is performed by HESO1. Biochemical experiments showed that the 3' terminus of miR1510 is only partially 2'-O-methylated because of the terminal mispairing in the miR1510/miR1510* duplex that inhibits HEN1 activity in soybean. miR1510 emerged in the Phaseoleae ∼41-42 million years ago with a conserved precursor structure yielding a 22-nt monouridylated form, yet a variant in mung bean is processed directly in a 22-nt mature form. This analysis of miR1510 yields two observations: (i) plants can utilize postprocessing modification to generate abundant 22-nt miRNA isoforms to more efficiently regulate target mRNA abundances; and (ii) comparative analysis demonstrates an example of selective optimization of precursor processing of a young plant miRNA.

Biogenesis of a 22-nt microRNA in Phaseoleae species by precursor-programmed uridylation

Phased, secondary siRNAs (phasiRNAs) represent a class of small RNAs in plants generated via distinct biogenesis pathways, predominantly dependent on the activity of 22-nt miRNAs. Most 22-nt miRNAs are processed by DCL1 from miRNA precursors containing an asymmetric bulge, yielding a 22/21-nt miRNA/miRNA* duplex. Here we show that miR1510, a soybean miRNA capable of triggering phasiRNA production from numerous nucleotide-binding leucine-rich repeat (NB-LRRs), previously described as 21 nt in its mature form, primarily accumulates as a 22-nt isoform via monouridylation. We demonstrate that, in Arabidopsis, this uridylation is performed by HESO1. Biochemical experiments showed that the 3' terminus of miR1510 is only partially 2'-O-methylated because of the terminal mispairing in the miR1510/miR1510* duplex that inhibits HEN1 activity in soybean. miR1510 emerged in the Phaseoleae ∼41-42 million years ago with a conserved precursor structure yielding a 22-nt monouridylated form, yet a variant in mung bean is processed directly in a 22-nt mature form. This analysis of miR1510 yields two observations: (i) plants can utilize postprocessing modification to generate abundant 22-nt miRNA isoforms to more efficiently regulate target mRNA abundances; and (ii) comparative analysis demonstrates an example of selective optimization of precursor processing of a young plant miRNA.

CRISPR/Cas9 and TALENs generate heritable mutations for genes involved in small RNA processing of Glycine max and Medicago truncatula

Processing of double-stranded RNA precursors into small RNAs is an essential regulator of gene expression in plant development and stress response. Small RNA processing requires the combined activity of a functionally diverse group of molecular components. However, in most of the plant species, there are insufficient mutant resources to functionally characterize each encoding gene. Here, mutations in loci encoding protein machinery involved in small RNA processing in soya bean and Medicago truncatula were generated using the CRISPR/Cas9 and TAL-effector nuclease (TALEN) mutagenesis platforms. An efficient CRISPR/Cas9 reagent was used to create a bi-allelic double mutant for the two soya bean paralogous Double-stranded RNA-binding2 (GmDrb2a and GmDrb2b) genes. These mutations, along with a CRISPR/Cas9-generated mutation of the M. truncatula Hua enhancer1 (MtHen1) gene, were determined to be germ-line transmissible. Furthermore, TALENs were used to generate a mutation within the soya bean Dicer-like2 gene. CRISPR/Cas9 mutagenesis of the soya bean Dicer-like3 gene and the GmHen1a gene was observed in the T0 generation, but these mutations failed to transmit to the T1 generation. The irregular transmission of induced mutations and the corresponding transgenes was investigated by whole-genome sequencing to reveal a spectrum of non-germ-line-targeted mutations and multiple transgene insertion events. Finally, a suite of combinatorial mutant plants were generated by combining the previously reported Gmdcl1a, Gmdcl1b and Gmdcl4b mutants with the Gmdrb2ab double mutant. Altogether, this study demonstrates the synergistic use of different genome engineering platforms to generate a collection of useful mutant plant lines for future study of small RNA processing in legume crops.

Biogenesis and regulatory hierarchy of phased small interfering RNAs in plants

Several varieties of small RNAs including microRNAs (miRNAs) and small interfering RNAs (siRNAs) are generated in plants to regulate development, genome stability and response to adverse environments. Phased siRNA (phasiRNA) is a type of secondary siRNA that is processed from a miRNA-mediated cleavage of RNA transcripts, increasing silencing efficiency or simultaneously suppressing multiple target genes. Trans-acting siRNAs (ta-siRNAs) are a particular class of phasiRNA produced from noncoding transcripts that silence targets in trans. It was originally thought that 'one-hit' and 'two-hit' models were essential for processing distinct TAS precursors; however, a single hit event was recently shown to be sufficient at triggering all types of ta-siRNAs. This review discusses the findings about biogenesis, targeting modes and regulatory networks of plant ta-siRNAs. We also summarize recent advances in the generation of other phasiRNAs and their possible biological benefits to plants.

Coupling of microRNA-directed phased small interfering RNA generation from long noncoding genes with alternative splicing and alternative polyadenylation in small RNA-mediated gene silencing

MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play vital regulatory roles in plant growth and development. Little is known about these small RNAs in litchi (Litchi chinensis), an economically important fruit crop widely cultivated in Southeast Asia. We profiled the litchi small RNA population with various deep-sequencing techniques and in-depth bioinformatic analyses. The genome-wide identification of miRNAs, their target genes, and phasiRNA-generating (PHAS) genes/loci showed that the function of miR482/2118 has expanded, relative to its canonical function. We also discovered that, for 29 PHAS loci, miRNA-mediated phasiRNA production was coupled with alternative splicing (AS) and alternative polyadenylation (APA). Most of these loci encoded long noncoding RNAs. An miR482/2118 targeted locus gave rise to four main transcript isoforms through AS/APA, and diverse phasiRNAs generated from these isoforms appeared to target long terminal repeat (LTR) retrotransposons and other unrelated genes. This coupling enables phasiRNA production from different exons of noncoding PHAS genes and yields diverse phasiRNA populations, both broadening and altering the range of downstream phasiRNA-regulated genes. Our results reveal the diversity of miRNA and phasiRNA in litchi, and demonstrate AS/APA as a new layer of regulation in small RNA-mediated gene silencing.

The Brassicaceae Family Displays Divergent, Shoot-Skewed NLR Resistance Gene Expression

Nucleotide-binding site leucine-rich repeat resistance genes (NLRs) allow plants to detect microbial effectors. We hypothesized that NLR expression patterns could reflect organ-specific differences in effector challenge and tested this by carrying out a meta-analysis of expression data for 1,235 NLRs from nine plant species. We found stable NLR root/shoot expression ratios within species, suggesting organ-specific hardwiring of NLR expression patterns in anticipation of distinct challenges. Most monocot and dicot plant species preferentially expressed NLRs in roots. In contrast, Brassicaceae species, including oilseed rape (Brassica napus) and the model plant Arabidopsis (Arabidopsis thaliana), were unique in showing NLR expression skewed toward the shoot across multiple phylogenetically distinct groups of NLRs. The Brassicaceae are also outliers in the sense that they have lost the common symbiosis signaling pathway, which enables intracellular infection by root symbionts. While it is unclear if these two events are related, the NLR expression shift identified here suggests that the Brassicaceae may have evolved unique pattern-recognition receptors and antimicrobial root metabolites to substitute for NLR protection. Such innovations in root protection could potentially be exploited in crop rotation schemes or for enhancing root defense systems of non-Brassicaceae crops.

Revisiting Criteria for Plant MicroRNA Annotation in the Era of Big Data

MicroRNAs (miRNAs) are ∼21-nucleotide-long regulatory RNAs that arise from endonucleolytic processing of hairpin precursors. Many function as essential posttranscriptional regulators of target mRNAs and long noncoding RNAs. Alongside miRNAs, plants also produce large numbers of short interfering RNAs (siRNAs), which are distinguished from miRNAs primarily by their biogenesis (typically processed from long double-stranded RNA instead of single-stranded hairpins) and functions (typically via roles in transcriptional regulation instead of posttranscriptional regulation). Next-generation DNA sequencing methods have yielded extensive data sets of plant small RNAs, resulting in many miRNA annotations. However, it has become clear that many miRNA annotations are questionable. The sheer number of endogenous siRNAs compared with miRNAs has been a major factor in the erroneous annotation of siRNAs as miRNAs. Here, we provide updated criteria for the confident annotation of plant miRNAs, suitable for the era of "big data" from DNA sequencing. The updated criteria emphasize replication and the minimization of false positives, and they require next-generation sequencing of small RNAs. We argue that improved annotation systems are needed for miRNAs and all other classes of plant small RNAs. Finally, to illustrate the complexities of miRNA and siRNA annotation, we review the evolution and functions of miRNAs and siRNAs in plants.

Phased secondary small interfering RNAs in Panaxnotoginseng

Background

Recent results demonstrated that either non-coding or coding genes generate phased secondary small interfering RNAs (phasiRNAs) guided by specific miRNAs. Till now, there is no studies for phasiRNAs in Panax notoginseng (Burk.) F.H. Chen (P. notoginseng), an important traditional Chinese herbal medicinal plant species.

Methods

Here we performed a genome-wide discovery of phasiRNAs and its host PHAS loci in P. notoginseng by analyzing small RNA sequencing profiles. Degradome sequencing profile was used to identify the trigger miRNAs of these phasiRNAs and potential targets of phasiRNAs. We also used RLM 5’-RACE to validate some of the identified phasiRNA targets.

Results

After analyzing 24 small RNA sequencing profiles of P. notoginseng, 204 and 90 PHAS loci that encoded 21 and 24 nucleotide (nt) phasiRNAs, respectively, were identified. Furthermore, we found that phasiRNAs produced from some pentatricopeptide repeat-contain (PPR) genes target another layer of PPR genes as validated by both the degradome sequencing profile and RLM 5’-RACE analysis. We also found that miR171 with 21 nt triggers the generations of 21 nt phasiRNAs from its conserved targets.

Conclusions

We validated that some phasiRNAs generated from PPRs and TASL genes are functional by targeting other PPRs in trans. These results provide the first set of PHAS loci and phasiRNAs in P. notoginseng, and enhance our understanding of PHAS in plants.

Electronic supplementary material

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

Small RNAs, emerging regulators critical for the development of horticultural traits

Small RNAs (sRNAs) have been recently recognized as key genetic and epigenetic regulators in various organisms, ranging from the modification of DNA and histone methylations to the modulation of the abundance of coding or non-coding RNAs. In plants, major regulatory sRNAs are classified as respective microRNA (miRNA) and small interfering RNA (siRNA) species, with the former primarily engaging in posttranscriptional regulation while the latter in transcriptional one. Many of these characterized sRNAs are involved in regulation of diverse biological programs, processes, and pathways in response to developmental cues, environmental signals/stresses, pathogen infection, and pest attacks. Recently, sRNAs-mediated regulations have also been extensively investigated in horticultural plants, with many novel mechanisms unveiled, which display far more mechanistic complexity and unique regulatory features compared to those studied in model species. Here, we review the recent progress of sRNA research in horticultural plants, with emphasis on mechanistic aspects as well as their relevance to trait regulation. Given that major and pioneered sRNA research has been carried out in the model and other plants, we also discuss ongoing sRNA research on these plants. Because miRNAs and phased siRNAs (phasiRNAs) are the most studied sRNA regulators, this review focuses on their biogenesis, conservation, function, and targeted genes and traits as well as the mechanistic relation between them, aiming at providing readers comprehensive information instrumental for future sRNA research in horticulture crops.

Salicylic Acid Perturbs sRNA-Gibberellin Regulatory Network in Immune Response of Potato to Potato virus Y Infection

Potato virus Y is the most economically important potato viral pathogen. We aimed at unraveling the roles of small RNAs (sRNAs) in the complex immune signaling network controlling the establishment of tolerant response of potato cv. Désirée to the virus. We constructed a sRNA regulatory network connecting sRNAs and their targets to link sRNA level responses to physiological processes. We discovered an interesting novel sRNAs-gibberellin regulatory circuit being activated as early as 3 days post inoculation (dpi) before viral multiplication can be detected. Two endogenous sRNAs, miR167 and phasiRNA931 were predicted to regulate gibberellin biosynthesis genes GA20-oxidase and GA3-oxidase. The increased expression of phasiRNA931 was also reflected in decreased levels of GA3-oxidase transcripts. Moreover, decreased concentration of gibberellin confirmed this regulation. The functional relation between lower activity of gibberellin signaling and reduced disease severity was previously confirmed in Arabidopsis-virus interaction using knockout mutants. We further showed that this regulation is salicylic acid-dependent as the response of sRNA network was attenuated in salicylic acid-depleted transgenic counterpart NahG-Désirée expressing severe disease symptoms. Besides downregulation of gibberellin signaling, regulation of immune receptor transcripts by miR6022 as well as upregulation of miR164, miR167, miR169, miR171, miR319, miR390, and miR393 in tolerant Désirée, revealed striking similarities to responses observed in mutualistic symbiotic interactions. The intertwining of different regulatory networks revealed, shows how developmental signaling, disease symptom development, and stress signaling can be balanced.

An expression atlas of miRNAs in Arabidopsis thaliana

MicroRNAs (miRNAs) are small non-coding RNAs that regulate a variety of biological processes. MiRNA expression often exhibits spatial and temporal specificity. However, genome-wide miRNA expression patterns in different organs during development of Arabidopsis thaliana have not yet been systemically investigated. In this study, we sequenced small RNA libraries generated from 27 different organ/tissue types, which cover the entire life cycle of Arabidopsis. Analysis of the sequencing data revealed that most miRNAs are ubiquitously expressed, whereas a small set of miRNAs display highly specific expression patterns. In addition, different miRNA members within the same family have distinct spatial and temporal expression patterns. Moreover, we found that some miRNAs are produced from different arms of their hairpin precursors at different developmental stages. This work provides new insights into the regulation of miRNA biogenesis and a rich resource for future investigation of miRNA functions in Arabidopsis.

Identification and molecular characterization of a trans-acting small interfering RNA producing locus regulating leaf rust responsive gene expression in wheat (Triticum aestivum L.)

A novel leaf rust responsive ta-siRNA-producing locus was identified in wheat showing similarity to 28S rRNA and generated four differentially expressing ta-siRNAs by phasing which targeted stress responsive genes. Trans-acting-small interfering RNAs (Ta-siRNAs) are plant specific molecules generally involved in development and are also stress responsive. Ta-siRNAs identified in wheat till date are all responsive to abiotic stress only. Wheat cultivation is severely affected by rusts and leaf rust particularly affects grain filling. This study reports a novel ta-siRNA producing locus (TAS) in wheat which is a segment of 28S ribosomal RNA but shows differential expression during leaf rust infestation. Four small RNA libraries prepared from wheat Near Isogenic Lines were treated with leaf rust pathogen and compared with untreated controls. A TAS with the ability to generate four ta-siRNAs by phasing events was identified along with the microRNA TamiR16 as the phase initiator. The targets of the ta-siRNAs included α-gliadin, leucine rich repeat, trans-membrane proteins, glutathione-S-transferase, and fatty acid desaturase among others, which are either stress responsive genes or are essential for normal growth and development of plants. Expression of the TAS, its generated ta-siRNAs, and their target genes were profiled at five different time points after pathogen inoculation of susceptible and resistant wheat isolines and compared with mock-inoculated controls. Comparative analysis of expression unveiled differential and reciprocal relationship as well as discrete patterns between susceptible and resistant isolines. The expression profiles of the target genes of the identified ta-siRNAs advocate more towards effector triggered susceptibility favouring pathogenesis. The study helps in discerning the functions of wheat genes regulated by ta-siRNAs in response to leaf rust.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Conservation and divergence of small RNA pathways and microRNAs in land plants

Background

As key regulators of gene expression in eukaryotes, small RNAs have been characterized in many seed plants, and pathways for their biogenesis, degradation, and action have been defined in model angiosperms. However, both small RNAs themselves and small RNA pathways are not well characterized in other land plants such as lycophytes and ferns, preventing a comprehensive evolutionary perspective on small RNAs in land plants.

Results

Using 25 representatives from major lineages of lycophytes and ferns, most of which lack sequenced genomes, we characterized small RNAs and small RNA pathways in these plants. We identified homologs of DICER-LIKE (DCL), ARGONAUTE (AGO), and other genes involved in small RNA pathways, predicted over 2600 conserved microRNA (miRNA) candidates, and performed phylogenetic analyses on small RNA pathways as well as miRNAs. Pathways underlying miRNA biogenesis, degradation, and activity were established in the common ancestor of land plants, but the 24-nucleotide siRNA pathway that guides DNA methylation is incomplete in sister species of seed plants, especially lycophytes. We show that the functional diversification of key gene families such as DCL and AGO as observed in angiosperms occurred early in land plants followed by parallel expansion of the AGO family in ferns and angiosperms. We uncovered a conserved AGO subfamily absent in angiosperms.

Conclusions

Our phylogenetic analyses of miRNAs in bryophytes, lycophytes, ferns, and angiosperms refine the time-of-origin for conserved miRNA families as well as small RNA machinery in land plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1291-2) contains supplementary material, which is available to authorized users.

The pivotal role of small non-coding RNAs in the regulation of seed development

Seeds represent a crucial stage of the seed plants life cycle. It is during seed development that the foundations of the future plant body, and the ability to give rise to a new plant capable of growing under sometimes adverse environmental conditions, are established. Small non-coding RNAs are major regulators of gene expression both at the post-transcriptional and transcriptional levels and, not surprisingly, these elements play major roles in seed development and germination. We review here the current knowledge about small RNA expression and functions in seed development, going from the morphogenesis phase comprehending embryo development and patterning, to the several steps of the maturation phase, ending in the transition to the germination. A special focus is given to the small RNAs for which functional studies have been conducted and their participation in regulatory networks operating in seeds. Many challenges remain ahead for dissecting the complex small RNA landscape in seeds, but this is a highly relevant issue in plant biology and advances in this area will most certainly impact plant breeding.

The functions of plant small RNAs in development and in stress responses

Like metazoans, plants use small regulatory RNAs (sRNAs) to direct gene expression. Several classes of sRNAs, which are distinguished by their origin and biogenesis, exist in plants. Among them, microRNAs (miRNAs) and trans-acting small interfering RNAs (ta-siRNAs) mainly inhibit gene expression at post-transcriptional levels. In the past decades, plant miRNAs and ta-siRNAs have been shown to be essential for numerous developmental processes, including growth and development of shoots, leaves, flowers, roots and seeds, among others. In addition, miRNAs and ta-siRNAs are also involved in the plant responses to abiotic and biotic stresses, such as drought, temperature, salinity, nutrient deprivation, bacteria, virus and others. This review summarizes the roles of miRNAs and ta-siRNAs in plant physiology and development.

Host-Pathogen interactions modulated by small RNAs

Biological processes such as defense mechanisms and microbial offense strategies are regulated through RNA induced interference in eukaryotes. Genetic mutations are modulated through biogenesis of small RNAs which directly impacts upon host development. Plant defense mechanisms are regulated and supported by a diversified group of small RNAs which are involved in streamlining several RNA interference pathways leading toward the initiation of pathogen gene silencing mechanisms. In the similar context, pathogens also utilize the support of small RNAs to launch their offensive attacks. Also there are strong evidences about the active involvement of these RNAs in symbiotic associations. Interestingly, small RNAs are not limited to the individuals in whom they are produced; they also show cross kingdom influences through variable interactions with other species thus leading toward the inter-organismic gene silencing. The phenomenon is understandable in the microbes which utilize these mechanisms to overcome host defense line. Understanding the mechanism of triggering host defense strategies can be a valuable step toward the generation of disease resistant host plants. We think that the cross kingdom trafficking of small RNA is an interesting insight that is needed to be explored for its vitality.

Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce

Lettuce (Lactuca sativa) is a major crop and a member of the large, highly successful Compositae family of flowering plants. Here we present a reference assembly for the species and family. This was generated using whole-genome shotgun Illumina reads plus in vitro proximity ligation data to create large superscaffolds; it was validated genetically and superscaffolds were oriented in genetic bins ordered along nine chromosomal pseudomolecules. We identify several genomic features that may have contributed to the success of the family, including genes encoding Cycloidea-like transcription factors, kinases, enzymes involved in rubber biosynthesis and disease resistance proteins that are expanded in the genome. We characterize 21 novel microRNAs, one of which may trigger phasiRNAs from numerous kinase transcripts. We provide evidence for a whole-genome triplication event specific but basal to the Compositae. We detect 26% of the genome in triplicated regions containing 30% of all genes that are enriched for regulatory sequences and depleted for genes involved in defence.

Mutations in Argonaute5 Illuminate Epistatic Interactions of the K1 and I Loci Leading to Saddle Seed Color Patterns in Glycine max

The soybean (Glycine max) seed coat has distinctive, genetically programmed patterns of pigmentation, and the recessive k1 mutation can epistatically overcome the dominant I and ii alleles, which inhibit seed color by producing small interfering RNAs (siRNAs) targeting chalcone synthase (CHS) mRNAs. Small RNA sequencing of dissected regions of immature seed coats demonstrated that CHS siRNA levels cause the patterns produced by the ii and ik alleles of the I locus, which restrict pigment to the hilum or saddle region of the seed coat, respectively. To identify the K1 locus, we compared RNA-seq data from dissected regions of two Clark isolines having similar saddle phenotypes mediated by CHS siRNAs but different genotypes (homozygous ik K1 versus homozygous ii k1). By examining differentially expressed genes, mapping information, and genome resequencing, we identified a 129-bp deletion in Glyma.11G190900 encoding Argonaute5 (AGO5), a member of the Argonaute family. Amplicon sequencing of several independent saddle pattern mutants from different genetic backgrounds revealed independent lesions affecting AGO5, thus establishing Glyma.11G190900 as the K1 locus. Nonfunctional AGO5 from k1 alleles leads to altered distributions of CHS siRNAs, thus explaining how the k1 mutation reverses the phenotype of the seed coat regions from yellow to pigmented, even in the presence of the normally dominant I or ii alleles.

Using genomic information to improve soybean adaptability to climate change

Climate change has brought severe challenges to agriculture. It is anticipated that there will be a drop in crop yield - including that of soybean - due to climatic stress factors that include drastic fluctuations in temperature, drought, flooding and high salinity. Genomic information on soybean has been accumulating rapidly since initial publication of its reference genome, providing a valuable tool for the improvement of cultivated soybean. Not only are many molecular markers that are associated with important quantitative trait loci now identified, but we also have a more detailed picture of the genomic variations among soybean germplasms, enabling us to utilize these as tools to assist crop breeding. In this review, we will summarize and discuss the currently available soybean genomic approaches, including whole-genome sequencing, sequencing-based genotyping, functional genomics, proteomics, and epigenomics. The information uncovered through these techniques will help further pinpoint important gene candidates and genomic loci associated with adaptive traits, as well as achieving a better understanding of how soybeans cope with the changing climate.

The legume miR1514a modulates a NAC transcription factor transcript to trigger phasiRNA formation in response to drought

Recent studies have identified microRNAs as post-transcriptional regulators involved in stress responses. miR1514a is a legume microRNA that is induced in response to drought stress in Phaseolus vulgaris (common bean) and shows differential accumulation levels in roots during water deficit in two cultivars with different drought tolerance phenotypes. A recent degradome analysis revealed that miR1514a targets the transcripts of two NAC transcription factors (TFs), Phvul.010g121000 and Phvul.010g120700. Furthermore, expression studies and small RNA-seq data indicate that only Phvul.010g120700 generates phasiRNAs, which also accumulate under water deficit conditions. To confirm these results, we over-expressed miR1514a in transgenic hairy roots, and observed a reduced accumulation of Phvul.010g120700 and an increase in NAC-derived phasiRNAs; inhibition of miR1514a activity resulted in the opposite effect. Moreover, we determined that a NAC-derived phasiRNA associates with ARGONAUTE 1 (AGO1), suggesting that it is functional. In addition, a transcriptome analysis of transgenic hairy roots with reduced miR1514a levels revealed several differentially expressed transcripts, mainly involved in metabolism and stress responses, suggesting they are regulated by the NAC TF and/or by phasiRNAs. This work therefore demonstrates the participation of miR1514 in the regulation of a NAC transcription factor transcript through phasiRNA production during the plant response to water deficit.

Insights into the function of the phasiRNA-triggering miR1514 in response to stress in legumes

We recently described the activity of miR1514a in response to water deficit in Phaseolus vulgaris. Pvu-miR1514a targets a NAC transcription factor mRNA for cleavage and subsequently triggers NAC-derived phasiRNA formation. Here we show that accumulation and activity of miR1514a are also conserved in the model legume Medicago truncatula. Consistently, we identified Mtr-miR1514a and detected its increased accumulation in response to stress conditions, targeting a NAC TF mRNA for cleavage and triggering phasiRNA production. In P. vulgaris, miR1514a inhibition in transgenic hairy roots was reported to increase NAC 700 mRNA levels and to affect expression patterns of several genes, including that of a Sec 14 homolog. We report here that in adult plant roots exposed to dehydration conditions, where miR1514a levels increased and NAC 700 mRNA decreased, there was a reduction of Sec 14 homolog mRNA levels, suggesting a direct transcriptional effect. The functions of miR1514a, NAC 700 and derived phasiRNAs have just begun to be elucidated in common bean; future understanding of their activities in this and other legumes species will advance our knowledge of microRNA functions in plants.

Biogenesis of diverse plant phasiRNAs involves an miRNA-trigger and Dicer-processing

It has been almost 30 years since RNA interference (RNAi) was shown to silence genes via double-stranded RNAs (dsRNAs) in Caenorhabditis elegans (Fire et al. 1998). 20-30-nucleotide (nt) small non-coding RNAs are a key element of the RNAi machinery. Recently, phased small interfering RNAs (phasiRNAs), small RNAs that are generated from a long RNA precursor at intervals of 21 to 26-nt, have been identified in plants and animals. In Drosophila, phasiRNAs are generated by the endonuclease, Zucchini (Zuc), in germlines. These phasiRNAs, known as one of PIWI-interacting RNAs (piRNAs), mainly repress transposable elements. Similarly, reproduction-specific phasiRNAs have been identified in the family Poaceae, although DICER LIKE (DCL) protein-dependent phasiRNA biogenesis in rice is distinct from piRNA biogenesis in animals. In plants, phasiRNA biogenesis is initiated when 22-nt microRNAs (miRNAs) cleave single-stranded target RNAs. Subsequently, RNA-dependent RNA polymerase (RDR) forms dsRNAs from the cleaved RNAs, and dsRNAs are further processed by DCLs into 21 to 24-nt phasiRNAs. Finally, the phasiRNAs are loaded to ARGONAUTE (AGO) proteins to induce RNA-silencing. There are diverse types of phasiRNA precursors and the miRNAs that trigger the biogenesis. Their expression patterns also differ among plant species, suggesting that species-specific combinations of these triggers dictate the spatio-temporal pattern of phasiRNA biogenesis during development, or in response to environmental stimuli.

PMS1T, producing phased small-interfering RNAs, regulates photoperiod-sensitive male sterility in rice

Phased small-interfering RNAs (phasiRNAs) are a special class of small RNAs, which are generated in 21- or 24-nt intervals from transcripts of precursor RNAs. Although phasiRNAs have been found in a range of organisms, their biological functions in plants have yet to be uncovered. Here we show that phasiRNAs generated by the photopheriod-sensetive genic male sterility 1 (Pms1) locus were associated with photoperiod-sensitive male sterility (PSMS) in rice, a germplasm that started the two-line hybrid rice breeding. The Pms1 locus encodes a long-noncoding RNA PMS1T that was preferentially expressed in young panicles. PMS1T was targeted by miR2118 to produce 21-nt phasiRNAs that preferentially accumulated in the PSMS line under long-day conditions. A single nucleotide polymorphism in PMS1T nearby the miR2118 recognition site was critical for fertility change, likely leading to differential accumulation of the phasiRNAs. This result suggested possible roles of phasiRNAs in reproductive development of rice, demonstrating the potential importance of this RNA class as regulators in biological processes.

Comparison of Small RNA Profiles of Glycine max and Glycine soja at Early Developmental Stages

Small RNAs, including microRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs; from PHAS loci), play key roles in plant development. Cultivated soybean, Glycine max, contributes a great deal to food production, but, compared to its wild kin, Glycine soja, it may lose some genetic information during domestication. In this work, we analyzed the sRNA profiles of different tissues in both cultivated (C08) and wild soybeans (W05) at three stages of development. A total of 443 known miRNAs and 15 novel miRNAs showed varying abundances between different samples, but the miRNA profiles were generally similar in both accessions. Based on a sliding window analysis workflow that we developed, 50 PHAS loci generating 55 21-nucleotide phasiRNAs were identified in C08, and 46 phasiRNAs from 41 PHAS loci were identified in W05. In germinated seedlings, phasiRNAs were more abundant in C08 than in W05. Disease resistant TIR-NB-LRR genes constitute a very large family of PHAS loci. PhasiRNAs were also generated from several loci that encode for NAC transcription factors, Dicer-like 2 (DCL2), Pentatricopeptide Repeat (PPR), and Auxin Signaling F-box 3 (AFB3) proteins. To investigate the possible involvement of miRNAs in initiating the PHAS-phasiRNA pathway, miRNA target predictions were performed and 17 C08 miRNAs and 15 W05 miRNAs were predicted to trigger phasiRNAs biogenesis. In summary, we provide a comprehensive description of the sRNA profiles of wild versus cultivated soybeans, and discuss the possible roles of sRNAs during soybean germination.

Small RNA profiles in soybean primary root tips under water deficit

BACKGROUND

Soybean (Glycine max) production is significantly hampered by frequent droughts in many regions of the world including the United States. Identifying microRNA (miRNA)-controlled posttranscriptional gene regulation under drought will enhance our understanding of molecular basis of drought tolerance in this important cash crop. Indeed, miRNA profiles in soybean exposed to drought were studied but not from the primary root tips, which is not only a main zone of water uptake but also critical for water stress sensing and signaling.

METHODS

Here we report miRNA profiles specifically from well-watered and water-stressed primary root tips (0 to 8 mm from the root apex) of soybean. Small RNA sequencing confirmed the expression of vastly diverse miRNA (303 individual miRNAs) population, and, importantly several conserved miRNAs were abundantly expressed in primary root tips.

RESULTS

Notably, 12 highly conserved miRNA families were differentially regulated in response to water-deficit; six were upregulated while six others were downregulated at least by one fold (log2) change. Differentially regulated soybean miRNAs are targeting genes include auxin response factors, Cu/Zn Superoxide dismutases, laccases and plantacyanin and several others.

CONCLUSIONS

These results highlighted the importance of miRNAs in primary root tips both under control and water-deficit conditions; under control conditions, miRNAs could be important for cell division, cell elongation and maintenance of the root apical meristem activity including quiescent centre whereas under water stress differentially regulated miRNAs could decrease auxin signaling and oxidative stress as well as other metabolic processes that save energy and water.

The Diversification of Plant NBS-LRR Defense Genes Directs the Evolution of MicroRNAs That Target Them

High expression of plant nucleotide binding site leucine-rich repeat (NBS-LRR) defense genes is often lethal to plant cells, a phenotype perhaps associated with fitness costs. Plants implement several mechanisms to control the transcript level of NBS-LRR defense genes. As negative transcriptional regulators, diverse miRNAs target NBS-LRRs in eudicots and gymnosperms. To understand the evolutionary benefits of this miRNA-NBS-LRR regulatory system, we investigated the NBS-LRRs of 70 land plants, coupling this analysis with extensive small RNA data. A tight association between the diversity of NBS-LRRs and miRNAs was found. The miRNAs typically target highly duplicated NBS-LRRs In comparison, families of heterogeneous NBS-LRRs were rarely targeted by miRNAs in Poaceae and Brassicaceae genomes. We observed that duplicated NBS-LRRs from different gene families periodically gave birth to new miRNAs. Most of these newly emerged miRNAs target the same conserved, encoded protein motif of NBS-LRRs, consistent with a model of convergent evolution for these miRNAs. By assessing the interactions between miRNAs and NBS-LRRs, we found nucleotide diversity in the wobble position of the codons in the target site drives the diversification of miRNAs. Taken together, we propose a co-evolutionary model of plant NBS-LRRs and miRNAs hypothesizing how plants balance the benefits and costs of NBS-LRR defense genes.

Identification of novel miRNAs from drought tolerant rice variety Nagina 22

MicroRNAs regulate a spectrum of developmental and biochemical processes in plants and animals. Thus, knowledge of the entire miRNome is essential to understand the complete regulatory schema of any organism. The current study attempts to unravel yet undiscovered miRNA genes in rice. Analysis of small RNA libraries from various tissues of drought-tolerant ‘aus’ rice variety Nagina 22 (N22) identified 71 novel miRNAs. These were validated based on precursor hairpin structure, small RNA mapping pattern, ‘star’ sequence, conservation and identification of targets based on degradome data. While some novel miRNAs were conserved in other monocots and dicots, most appear to be lineage-specific. They were segregated into two different classes based on the closeness to the classical miRNA definition. Interestingly, evidence of a miRNA-like cleavage was found even for miRNAs that lie beyond the classical definition. Several novel miRNAs displayed tissue-enriched and/or drought responsive expression. Generation and analysis of the degradome data from N22 along with publicly available degradome identified several high confidence targets implicated in regulation of fundamental processes such as flowering and stress response. Thus, discovery of these novel miRNAs considerably expands the dimension of the miRNA-mediated regulation in rice.

Regulation of Small RNAs and Corresponding Targets in Nod Factor-Induced Phaseolus vulgaris Root Hair Cells

A genome-wide analysis identified the set of small RNAs (sRNAs) from the agronomical important legume Phaseolus vulgaris (common bean), including novel P. vulgaris-specific microRNAs (miRNAs) potentially important for the regulation of the rhizobia-symbiotic process. Generally, novel miRNAs are difficult to identify and study because they are very lowly expressed in a tissue- or cell-specific manner. In this work, we aimed to analyze sRNAs from common bean root hairs (RH), a single-cell model, induced with pure Rhizobium etli nodulation factors (NF), a unique type of signal molecule. The sequence analysis of samples from NF-induced and control libraries led to the identity of 132 mature miRNAs, including 63 novel miRNAs and 1984 phasiRNAs. From these, six miRNAs were significantly differentially expressed during NF induction, including one novel miRNA: miR-RH82. A parallel degradome analysis of the same samples revealed 29 targets potentially cleaved by novel miRNAs specifically in NF-induced RH samples; however, these novel miRNAs were not differentially accumulated in this tissue. This study reveals Phaseolus vulgaris-specific novel miRNA candidates and their corresponding targets that meet all criteria to be involved in the regulation of the early nodulation events, thus setting the basis for exploring miRNA-mediated improvement of the common bean–rhizobia symbiosis.

A naïve Bayesian classifier for identifying plant microRNAs

MicroRNAs (miRNAs) are important regulatory molecules in eukaryotic organisms. Existing methods for the identification of mature miRNA sequences in plants rely extensively on the search for stem-loop structures, leading to high false negative rates. Here, we describe a probabilistic method for ranking putative plant miRNAs using a naïve Bayes classifier and its publicly available implementation. We use a number of properties to construct the classifier, including sequence length, number of observations, existence of detectable predicted miRNA* sequences, the distribution of nearby reads and mapping multiplicity. We apply the method to small RNA sequence data from soybean, peach, Arabidopsis and rice and provide experimental validation of several predictions in soybean. The approach performs well overall and strongly enriches for known miRNAs over other types of sequences. By utilizing a Bayesian approach to rank putative miRNAs, our method is able to score miRNAs that would be eliminated by other methods, such as those that have low counts or lack detectable miRNA* sequences. As a result, we are able to detect several soybean miRNA candidates, including some that are 24 nucleotides long, a class that is almost universally eliminated by other methods.

[Letter to the Editor] Comparison of small RNA next-generation sequencing with and without isolation of small RNA fraction

Address correspondence to Ferenc Marincs, Agricultural Biotechnology Institute, NARIC, Szent-Györgyi Albert u. 4., Gödöllö, Hungary. E-mail: marincs.ferenc@abc.naik.hu.

Small RNA and PARE sequencing in flower bud reveal the involvement of sRNAs in endodormancy release of Japanese pear (Pyrus pyrifolia 'Kosui')

Background

In woody perennial plants, including deciduous fruit trees, such as pear, endodormancy is a strategy for surviving the cold winter. A better understanding of the mechanism underlying the endodormancy phase transition is necessary for developing countermeasures against the effects of global warming. In this study, we analyzed the sRNAome of Japanese pear flower buds in endodormant and ecodormant stages over two seasons by implementing of RNA-seq and degradome-sequencing.

Results

We identified 137 conserved or less conserved miRNAs and 50 pear-specific miRNAs. However, none of the conserved microRNAs or pear-specific miRNAs was differentially expressed between endodormancy and ecodormancy stages. On the contrast, 1540 of 218,050 loci that produced sRNAs were differentially expressed between endodormancy and ecodormancy, suggesting their potential roles on the phase transition from endodormancy to ecodomancy. We also characterized a multifunctional miRNA precursor MIR168, which produces two functional miR168 transcripts, namely miR168.1 and miR168.2; cleavage events were predominantly mediated by the non-conserved variant miR168.2 rather than the conserved variant miR168.1. Finally, we showed that a TAS3 trans-acting siRNA triggered phased siRNA within the ORF of one of its target genes, AUXIN RESPONSE FACTOR 4, via the analysis of phased siRNA loci, indicating that siRNAs are able to trigger phased siRNAs in pear.

Conclusion

We analyzed the sRNAome of pear flower bud during dormant phase transition. Our work described the sRNA profiles of pear winter buds during dormant phase transition, showing that dormancy release is a highly coordinated physiological process involving the regulation of sRNAs.

Electronic supplementary material

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

Small RNAs Add Zing to the Zig-Zag-Zig Model of Plant Defenses

Plant small RNAs play important roles in transcriptional and posttranscriptional regulation, with ongoing work demonstrating their functions in diverse pathways. Their roles in defense responses are a topic of active investigation, particularly the rich set of micro (mi)RNAs that target disease resistance genes such as nucleotide binding/leucine-rich repeat (NB-LRR) genes. The miRNA-NB-LRR interactions result in the production of phased, secondary small interfering (phasi)RNAs, and phasiRNAs function in both cis and trans to propagate negative regulatory effects across additional members of the target gene family. Yet, while phasiRNAs have the capacity to trigger targeted decay of specific targets, both in cis and trans, their functional relevance in NB-LRR regulation remains largely a matter of speculation.

Noncoding RNAs, Emerging Regulators in Root Endosymbioses

Endosymbiosis interactions allow plants to grow in nutrient-deficient soil environments. The arbuscular mycorrhizal (AM) symbiosis is an ancestral interaction between land plants and fungi, whereas nitrogen-fixing symbioses are highly specific for certain plants, notably major crop legumes. The signaling pathways triggered by specific lipochitooligosaccharide molecules involved in these interactions have common components that also overlap with plant root development. These pathways include receptor-like kinases, transcription factors (TFs), and various intermediate signaling effectors, including noncoding (nc)RNAs. These latter molecules have emerged as major regulators of gene expression and small ncRNAs, composed of micro (mi)RNAs and small interfering (si)RNAs, are known to control gene expression at transcriptional (chromatin) or posttranscriptional levels. In this review, we describe exciting recent data connecting variants of conserved si/miRNAs with the regulation of TFs, such as NSP2, NFY-A1, auxin-response factors, and AP2-like proteins, known to be involved in symbiosis. The link between hormonal regulations and these si- and miRNA-TF nodes is proposed in a model in which different feedback loops or regulations controlling endosymbiosis signaling are integrated. The diversity and emerging regulatory networks of young legume miRNAs are also highlighted.

Evolution of plant genome architecture

We have witnessed an explosion in our understanding of the evolution and structure of plant genomes in recent years. Here, we highlight three important emergent realizations: (1) that the evolutionary history of all plant genomes contains multiple, cyclical episodes of whole-genome doubling that were followed by myriad fractionation processes; (2) that the vast majority of the variation in genome size reflects the dynamics of proliferation and loss of lineage-specific transposable elements; and (3) that various classes of small RNAs help shape genomic architecture and function. We illustrate ways in which understanding these organism-level and molecular genetic processes can be used for crop plant improvement.

Global investigation of the co-evolution of MIRNA genes and microRNA targets during soybean domestication

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome-wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA-target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA-target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA-target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA-target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co-evolution of MIRs and miRNA targets during soybean domestication.