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.
Keywords: Glycine max; Glycine soja; PHAS loci; PHAS-phasiRNA pathway; miRNA; phasiRNA; sliding window analysis.