An Atlas of Soybean Small RNAs Identifies Phased siRNAs from Hundreds of Coding Genes

Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.     

Comprehensive Annotation of Physcomitrella patens Small RNA Loci Reveals That the Heterochromatic Short Interfering RNA Pathway Is Largely Conserved in Land Plants

Many plant small RNAs are sequence-specific negative regulators of target mRNAs and/or chromatin. In angiosperms, the two most abundant endogenous small RNA populations are usually 21-nucleotide microRNAs (miRNAs) and 24-nucleotide heterochromatic short interfering RNAs (siRNAs). Heterochromatic siRNAs are derived from repetitive regions and reinforce DNA methylation at targeted loci. The existence and extent of heterochromatic siRNAs in other land plant lineages has been unclear. Using small RNA-sequencing (RNA-seq) of the moss Physcomitrella patens, we identified 1090 loci that produce mostly 23- to 24-nucleotide siRNAs. These loci are mostly in intergenic regions with dense DNA methylation. Accumulation of siRNAs from these loci depends upon P. patens homologs of DICER-LIKE3 (DCL3), RNA-DEPENDENT RNA POLYMERASE2, and the largest subunit of DNA-DEPENDENT RNA POLYMERASE IV, with the largest subunit of a Pol V homolog contributing to expression at a smaller subset of the loci. A MINIMAL DICER-LIKE (mDCL) gene, which lacks the N-terminal helicase domain typical of DCL proteins, is specifically required for 23-nucleotide siRNA accumulation. We conclude that heterochromatic siRNAs, and their biogenesis pathways, are largely identical between angiosperms and P. patens, with the notable exception of the P. patens-specific use of mDCL to produce 23-nucleotide siRNAs.     

Arabidopsis RNASE THREE LIKE2 Modulates the Expression of Protein-Coding Genes via 24-Nucleotide Small Interfering RNA-Directed DNA Methylation

RNaseIII enzymes catalyze the cleavage of double-stranded RNA (dsRNA) and have diverse functions in RNA maturation. Arabidopsis thaliana RNASE THREE LIKE2 (RTL2), which carries one RNaseIII and two dsRNA binding (DRB) domains, is a unique Arabidopsis RNaseIII enzyme resembling the budding yeast small interfering RNA (siRNA)-producing Dcr1 enzyme. Here, we show that RTL2 modulates the production of a subset of small RNAs and that this activity depends on both its RNaseIII and DRB domains. However, the mode of action of RTL2 differs from that of Dcr1. Whereas Dcr1 directly cleaves dsRNAs into 23-nucleotide siRNAs, RTL2 likely cleaves dsRNAs into longer molecules, which are subsequently processed into small RNAs by the DICER-LIKE enzymes. Depending on the dsRNA considered, RTL2-mediated maturation either improves (RTL2-dependent loci) or reduces (RTL2-sensitive loci) the production of small RNAs. Because the vast majority of RTL2-regulated loci correspond to transposons and intergenic regions producing 24-nucleotide siRNAs that guide DNA methylation, RTL2 depletion modifies DNA methylation in these regions. Nevertheless, 13% of RTL2-regulated loci correspond to protein-coding genes. We show that changes in 24-nucleotide siRNA levels also affect DNA methylation levels at such loci and inversely correlate with mRNA steady state levels, thus implicating RTL2 in the regulation of protein-coding gene expression.     

Failure of the Tomato Trans-Acting Short Interfering RNA Program to Regulate AUXIN RESPONSE FACTOR3 and ARF4 Underlies the Wiry Leaf Syndrome

Interfering with small RNA production is a common strategy of plant viruses. A unique class of small RNAs that require microRNA and short interfering (siRNA) biogenesis for their production is termed trans-acting short interfering RNAs (ta-siRNAs). Tomato (Solanum lycopersicum) wiry mutants represent a class of phenotype that mimics viral infection symptoms, including shoestring leaves that lack leaf blade expansion. Here, we show that four WIRY genes are involved in siRNA biogenesis, and in their corresponding mutants, levels of ta-siRNAs that regulate AUXIN RESPONSE FACTOR3 (ARF3) and ARF4 are reduced, while levels of their target ARFs are elevated. Reducing activity of both ARF3 and ARF4 can rescue the wiry leaf lamina, and increased activity of either can phenocopy wiry leaves. Thus, a failure to negatively regulate these ARFs underlies tomato shoestring leaves. Overexpression of these ARFs in Arabidopsis thaliana, tobacco (Nicotiana tabacum), and potato (Solanum tuberosum) failed to produce wiry leaves, suggesting that the dramatic response in tomato is exceptional. As negative regulation of orthologs of these ARFs by ta-siRNA is common to land plants, we propose that ta-siRNA levels serve as universal sensors for interference with small RNA biogenesis, and changes in their levels direct species-specific responses.