Soybean miRNAs Research Articles

Genome-Wide Investigation Using sRNA-Seq, Degradome-Seq and Transcriptome-Seq Reveals Regulatory Networks of microRNAs and Their Target Genes in Soybean during Soybean mosaic virus Infection.

MicroRNAs (miRNAs) play key roles in a variety of cellular processes through regulation of their target gene expression. Accumulated experimental evidence has demonstrated that infections by viruses are associated with the altered expression profile of miRNAs and their mRNA targets in the host. However, the regulatory network of miRNA-mRNA interactions during viral infection remains largely unknown. In this study, we performed small RNA (sRNA)-seq, degradome-seq and as well as a genome-wide transcriptome analysis to profile the global gene and miRNA expression in soybean following infections by three different Soybean mosaic virus (SMV) isolates, L (G2 strain), LRB (G2 strain) and G7 (G7 strain). sRNA-seq analyses revealed a total of 253 soybean miRNAs with a two-fold or greater change in abundance compared with the mock-inoculated control. 125 transcripts were identified as the potential cleavage targets of 105 miRNAs and validated by degradome-seq analyses. Genome-wide transcriptome analysis showed that total 2679 genes are differentially expressed in response to SMV infection including 71 genes predicted as involved in defense response. Finally, complex miRNA-mRNA regulatory networks were derived using the RNAseq, small RNAseq and degradome data. This work represents a comprehensive, global approach to examining virus-host interactions. Genes responsive to SMV infection are identified as are their potential miRNA regulators. Additionally, regulatory changes of the miRNAs themselves are described and the regulatory relationships were supported with degradome data. Taken together these data provide new insights into molecular SMV-soybean interactions and offer candidate miRNAs and their targets for further elucidation of the SMV infection process.

MicroRNA Maturation and MicroRNA Target Gene Expression Regulation Are Severely Disrupted in Soybean dicer-like1 Double Mutants.

Small nonprotein-coding microRNAs (miRNAs) are present in most eukaryotes and are central effectors of RNA silencing-mediated mechanisms for gene expression regulation. In plants, DICER-LIKE1 (DCL1) is the founding member of a highly conserved family of RNase III-like endonucleases that function as core machinery proteins to process hairpin-like precursor transcripts into mature miRNAs, small regulatory RNAs, 21–22 nucleotides in length. Zinc finger nucleases (ZFNs) were used to generate single and double-mutants of putative soybean DCL1 homologs, DCL1a and DCL1b, to confirm their functional role(s) in the soybean miRNA pathway. Neither DCL1 single mutant, dcl1a or dcl1b plants, exhibited a pronounced morphological or molecular phenotype. However, the dcl1a/dcl1b double mutant expressed a strong morphological phenotype, characterized by reduced seed size and aborted seedling development, in addition to defective miRNA precursor transcript processing efficiency and deregulated miRNA target gene expression. Together, these findings indicate that the two soybean DCL1 paralogs, DCL1a and DCL1b, largely play functionally redundant roles in the miRNA pathway and are essential for normal plant development.

Identification of miRNAs and their targets by high-throughput sequencing and degradome analysis in cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B of soybean.

BackgroundCytoplasmic male sterility (CMS) provides crucial breeding materials that facilitate hybrid seed production in various crops, and thus plays an important role in the study of hybrid vigor (heterosis), in plants. However, the CMS regulatory network in soybean remains unclear. MicroRNAs (miRNAs) play crucial roles in flower and pollen development by targeting genes that regulate their expression in plants. To identify the miRNAs and their targets that exist in the soybean CMS line NJCMS1A and its maintainer NJCMS1B, high-throughput sequencing and degradome analysis were conducted in this study.ResultsTwo small RNA libraries were constructed from the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. A total of 105 new miRNAs present on the other arm of known pre-miRNAs, 23 new miRNA members, 158 novel miRNAs and 160 high-confidence soybean miRNAs were identified using high-throughput sequencing. Among the identified miRNAs, 101 differentially expressed miRNAs with greater than two-fold changes between NJCMS1A and NJCMS1B were discovered. The different expression levels of selected miRNAs were confirmed by stem-loop quantitative real-time PCR. A degradome analysis showed that 856 targets were predicted to be targeted by 296 miRNAs, including a squamosa promoter-binding protein-like transcription factor family protein, a pentatricopeptide repeat-containing protein, and an auxin response factor, which were previously shown to be involved in floral organ or anther development in plants. Additionally, some targets, including a MADS-box transcription factor, NADP-dependent isocitrate dehydrogenase and NADH-ubiquinone oxidoreductase 24 kDa subunit, were identified, and they may have some relationship with the programmed cell death, reactive oxygen species accumulation and energy deficiencies, which might lead to soybean male sterility.ConclusionsThe present study is the first to use deep sequencing technology to identify miRNAs and their targets in the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. The results revealed that the miRNAs might participate in flower and pollen development, which could facilitate our understanding of the molecular mechanisms behind CMS in soybean.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2352-0) contains supplementary material, which is available to authorized users.

Novel and conserved microRNAs in soybean floral whorls

MicroRNAs (miRNAs) correspond to a class of endogenous small non-coding RNAs (19–24nt) that regulates the gene expression, through mRNA target cleavage or translation inhibition. In plants, miRNAs have been shown to play pivotal roles in a wide variety of metabolic and biological processes like plant growth, development, and response to biotic and abiotic stress. Soybean is one of the most important crops worldwide, due to the production of oil and its high protein content. The reproductive phase is considered the most important for soybean yield, which is mainly intended to produce the grains. The identification of miRNAs is not yet saturated in soybean, and there are no studies linking them to the different floral organs. In this study, three different mature soybean floral whorls were used in the construction of sRNA libraries. The sequencing of petal, carpel and stamen libraries generated a total of 10,165,661 sequences. Subsequent analyses identified 200 miRNAs sequences, among which, 41 were novel miRNAs, 80 were conserved soybean miRNAs, 31 were new antisense conserved soybean miRNAs and 46 were soybean miRNAs isoforms. We also found a new miRNA conserved in other plant species, and finally one miRNA-sibling of a soybean conserved miRNA. Conserved and novel miRNAs were evaluated by RT-qPCR. We observed a differential expression across the three whorls for six miRNAs. Computational predicted targets for miRNAs analyzed by RT-qPCR were identified and present functions related to reproductive process in plants. In summary, the increased accumulation of specific and novel miRNAs in different whorls indicates that miRNAs are an important part of the regulatory network in soybean flower.

Functional marker development of miR1511-InDel and allelic diversity within the genus Glycine.

BackgroundSingle-stranded non-protein coding small RNAs, 18–25 nucleotides in length, are ubiquitous throughout plants genomes and are involved in post-transcriptional gene regulation. Several types of DNA markers have been reported for the detection of genetic diversity or sequence variation in soybean, one of the most important legume crops in worldwide for seed protein and oil content. Recently, with the available of public genomic databases, there has been a shift from the labor-intensive development of PCR-based markers to sequence-based genotyping and the development of functional markers within genes, often coupled with the use of RNA information. But thus far miRNA-based markers have been only developed in rice and tobacco. Here we report the first functional molecular miRNA marker, miR1511-InDel, in soybean for a specific single copy locus used to assess genetic variation in domesticated soybean (Glycine max [L.] Merr) and its wild progenitor (Glycine soja Sieb. & Zucc.).ResultsWe genotyped a total of 1,669 accessions of domesticated soybean (G. max) and its wild progenitor G. soja which are native throughout the China and parts of Korea, Japan and Russia. The results indicate that the miR1511 locus is distributed in cultivated soybean and has three alleles in annual wild soybean. Based on this result, we proposed that miR-InDel marker technology can be used to assess genetic variation. The inclusion of geo-reference data with miR1511-InDel marker data corroborated that accessions from the Yellow River basin (Huanghuai) exhibited high genetic diversity which provides more molecular evidence for gene diversity in annual wild soybean and domestication of soybean.ConclusionsThese results provide evidence for the use of RNA marker, miRNA1511-InDel, as a soybean-specific functional maker for the study of genetic diversity, genotyping of germplasm and evolution studies. This is also the first report of functional marker developed from soybean miRNA located within the functional region of pre-miRNA1511.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1665-3) contains supplementary material, which is available to authorized users.

Identification and functional characterization of soybean root hair microRNAs expressed in response to Bradyrhizobium japonicum infection.

Three soybean [Glycine max (L) Merr.] small RNA libraries were generated and sequenced using the Illumina platform to examine the role of miRNAs during soybean nodulation. The small RNA libraries were derived from root hairs inoculated with Bradyrhizobium japonicum (In_RH) or mock-inoculated with water (Un_RH), as well as from the comparable inoculated stripped root samples (i.e. inoculated roots with the root hairs removed). Sequencing of these libraries identified a total of 114 miRNAs, including 22 novel miRNAs. A comparison of miRNA abundance among the 114 miRNAs identified 66 miRNAs that were differentially expressed between root hairs and stripped roots, and 48 miRNAs that were differentially regulated in infected root hairs in response to B. japonicum when compared to uninfected root hairs (P ≤ 0.05). A parallel analysis of RNA ends (PARE) library was constructed and sequenced to reveal a total of 405 soybean miRNA targets, with most predicted to encode transcription factors or proteins involved in protein modification, protein degradation and hormone pathways. The roles of gma-miR4416 and gma-miR2606b during nodulation were further analysed. Ectopic expression of these two miRNAs in soybean roots resulted in significant changes in nodule numbers. miRNA target information suggested that gma-miR2606b regulates a Mannosyl-oligosaccharide 1, 2-alpha-mannosidase gene, while gma-miR4416 regulates the expression of a rhizobium-induced peroxidase 1 (RIP1)-like peroxidase gene, GmRIP1, during nodulation.

Identification of microRNAs and their mRNA targets during soybean nodule development: functional analysis of the role of miR393j-3p in soybean nodulation.

Plant microRNAs (miRNAs) play important regulatory roles in a number of developmental processes. The present work investigated the roles of miRNAs during nodule development in the crop legume soybean (Glycine max). Fifteen soybean small RNA libraries were sequenced from different stages of nodule development, including young nodules, mature nodules and senescent nodules. In order to identify the regulatory targets of the miRNAs, five parallel analysis of RNA ends (PARE) libraries were also sequenced from the same stages of nodule development. Sequencing identified 284 miRNAs, including 178 novel soybean miRNAs. Analysis of miRNA abundance identified 139 miRNAs whose expression was significantly regulated during nodule development, including 12 miRNAs whose expression changed > 10-fold. Analysis of the PARE libraries identified 533 miRNA targets, including three nodulation-related genes and eight nodule-specific genes. miR393j-3p was selected for detailed analysis as its expression was significantly regulated during nodule formation, and it targeted a nodulin gene, Early Nodulin 93 (ENOD93). Strong, ectopic expression of miR393j-3p, as well as RNAi silencing of ENOD93 expression, significantly reduced nodule formation. The data indicate that miR393j-3p regulation of ENOD93 mRNA abundance is a key control point for soybean nodule formation.

Evolution Analysis About Soybean MIR166 Family

MicroRNA (miRNA) is a class of important regulating non-coding small molecular RNA. The gma-MIR166 gene family consists of 21 members and their expression patterns diversify widely. It is important to analyze the evolution of gma-MIR166 gene family in order to understand the evolutionary mechanisms of miRNAs in soybean. In this study, we implemented soybean wide genome block analysis, the molecular phylogeny of gma-MIR166 and block analysis of gma-MIR166 family. The results showed that both chromosome big segmental duplications and tandem duplications were main reasons contributed to the expanding of gma-MIR166 gene family. These findings suggested that gma-MIR166 gene family might originate from one or two ancient miRNA genes. The results of research provided a support for evolutionary study of miRNAs in soybean and related species in Fabaceae.

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.

Knockdown of the 7S globulin subunits shifts distribution of nitrogen sources to the residual protein fraction in transgenic soybean seeds.

A platform of gene silencing by amiRNA had been established in fertile transgenic soybean. We demonstrated that knockdown of storage protein shifted the distribution of nitrogen sources in soybean seeds. Artificial microRNAs (amiRNAs) were designed using the precursor sequence of the endogenous soybean (Glycine max L. Merrill) miRNA gma-miR159a and expressed in transgenic soybean plants to suppress the biosynthesis of 7S globulin, which is one of the major storage proteins. Seed-specific expression of these amiRNAs (amiR-7S) resulted in a strong suppression of 7S globulin subunit genes and decreased accumulation of the 7S globulin subunits in seeds. Thus, the results demonstrate that a platform for gene silencing by amiRNA was first developed in fertile transgenic soybean plants. There was no difference in nitrogen, carbon, and lipid contents between amiR-7S and control seeds. Four protein fractions were collected from defatted mature seeds on the basis of solubility at different pH to examine the distribution of nitrogen sources and compensatory effects. In the whey and lipophilic fractions, nitrogen content was similar in amiR-7S and control seeds. Nitrogen content was significantly decreased in the major soluble protein fraction and increased in the residual fraction (okara) of the amiR-7S seeds. Amino acid analysis revealed that increased nitrogen compounds in okara were proteins or peptides rather than free amino acids. Our study indicates that the decrease in 7S globulin subunits shifts the distribution of nitrogen sources to okara in transgenic soybean seeds.

Computational identification of novel microRNAs and targets in Glycine max.

Plant miRNAs, the endogenous non-coding small RNAs of about 20-24 nucleotides, play important roles in multiple biological processes by acting as negative regulators of their targeted mRNAs. Soybean (Glycine max (L.) Merr.) is one of the important oil crops of the world, in which many miRNAs have been obtained through the computational prediction or experiments. However, the miRNA genes identified for soybean are still far from saturation, and their biological functions are largely unknown. Here, a total of 48 candidates of miRNAs were identified following a range of strict filtering criteria. Detailed sequence analysis showed that G. max pre-miRNAs vary in length from 47 to 380 nt, embody mature miRNAs that differ in their physical location within the pre-miRNAs. In this study, twenty miRNAs were confirmed by microarray and three miRNAs were further validated by poly(A)-tailed RT-PCR. Comparative sequence analysis of soybean miRNA sequences showed that uracil is the dominant base in the first position at the 5' end of the mature miRNAs, and the base may have an important functional role in miRNA biogenesis and/or miRNA-mediated gene regulation. Finally, we predicted potential targets of these miRNAs. These target genes were predicted to encode transcription factors, resistance protein, heat shock protein, protein kinase, transporter, zinc finger protein and others, which might play important roles in soybean development and stress response.

Analysis of promoters of microRNAs from a Glycine max degradome library

MicroRNAs (miRNAs) are genome-encoded, small non-coding RNAs that play important functions in development, biotic and abiotic stress responses, and other processes. Our aim was to explore the regulation of miRNA expression. We used bioinformatics methods to predict the core promoters of 440 miRNAs identified from a soybean (Glycine max) degradome library and to analyze cis-acting elements for 369 miRNAs. The prediction results showed that 83.86% of the 440 miRNAs contained promoters in their upstream sequences, and 8.64% (38 loci) in their downstream sequences. The distributions of two core promoter elements, TATA-boxes and transcription start sites (TSSs), were similar. The cis-acting elements were examined to provide clues to the function and regulation of spatiotemporal expression of the miRNAs. Analyses of miRNA cis-elements and targets indicated a potential auxin response factor (ARF)- and gibberellin response factor (GARF)-mediated negative feedback loop for miRNA expression. The features of miRNAs from a Glycine max degradome library obtained here provide insights into the transcription regulation and functions of miRNAs in soybean.

Systems and Evolutionary Characterization of MicroRNAs and Their Underlying Regulatory Networks in Soybean Cotyledons

MicroRNAs (miRNAs) are an emerging class of small RNAs regulating a wide range of biological processes. Soybean cotyledons evolved as sink tissues to synthesize and store seed reserves which directly affect soybean seed yield and quality. However, little is known about miRNAs and their regulatory networks in soybean cotyledons. We sequenced 292 million small RNA reads expressed in soybean cotyledons, and discovered 130 novel miRNA genes and 72 novel miRNA families. The cotyledon miRNAs arose at various stages of land plant evolution. Evolutionary analysis of the miRNA genes in duplicated genome segments from the recent Glycine whole genome duplication revealed that the majority of novel soybean cotyledon miRNAs were young, and likely arose after the duplication event 13 million years ago. We revealed the evolutionary pathway of a soybean cotyledon miRNA family (soy-miR15/49) that evolved from a neutral invertase gene through an inverted duplication and a series of DNA amplification and deletion events. A total of 304 miRNA genes were expressed in soybean cotyledons. The miRNAs were predicted to target 1910 genes, and form complex miRNA networks regulating a wide range of biological pathways in cotyledons. The comprehensive characterization of the miRNAs and their underlying regulatory networks at gene, pathway and system levels provides a foundation for further studies of miRNAs in cotyledons.

SoyFN: a knowledge database of soybean functional networks

Many databases for soybean genomic analysis have been built and made publicly available, but few of them contain knowledge specifically targeting the omics-level gene–gene, gene–microRNA (miRNA) and miRNA–miRNA interactions. Here, we present SoyFN, a knowledge database of soybean functional gene networks and miRNA functional networks. SoyFN provides user-friendly interfaces to retrieve, visualize, analyze and download the functional networks of soybean genes and miRNAs. In addition, it incorporates much information about KEGG pathways, gene ontology annotations and 3′-UTR sequences as well as many useful tools including SoySearch, ID mapping, Genome Browser, eFP Browser and promoter motif scan. SoyFN is a schema-free database that can be accessed as a Web service from any modern programming language using a simple Hypertext Transfer Protocol call. The Web site is implemented in Java, JavaScript, PHP, HTML and Apache, with all major browsers supported. We anticipate that this database will be useful for members of research communities both in soybean experimental science and bioinformatics.Database URL: http://nclab.hit.edu.cn/SoyFN

Identification and Comparative Analysis of Cadmium Tolerance-Associated miRNAs and Their Targets in Two Soybean Genotypes

MicroRNAs (miRNAs) play crucial roles in regulating the expression of various stress responses genes in plants. To investigate soybean (Glycine max) miRNAs involved in the response to cadmium (Cd), microarrays containing 953 unique miRNA probes were employed to identify differences in the expression patterns of the miRNAs between different genotypes, Huaxia3 (HX3, Cd-tolerant) and Zhonghuang24 (ZH24, Cd-sensitive). Twenty six Cd-responsive miRNAs were identified in total. Among them, nine were detected in both cultivars, while five were expressed only in HX3 and 12 were only in ZH24. The expression of 16 miRNAs was tested by qRT-PCR and most of the identified miRNAs were found to have similar expression patterns with microarray. Three hundred and seventy six target genes were identified for 204 miRNAs from a mixture degradome library, which was constructed from the root of HX3 and ZH24 with or without Cd treatment. Fifty five genes were identified to be cleaved by 14 Cd-responsive miRNAs. Gene ontology (GO) annotations showed that these target transcripts are implicated in a broad range of biological processes. In addition, the expression patterns of ten target genes were validated by qRT-PCR. The characterization of the miRNAs and the associated target genes in response to Cd exposure provides a framework for understanding the molecular mechanism of heavy metal tolerance in plants.

Inferring the soybean (Glycine max) microRNA functional network based on target gene network

The rapid accumulation of microRNAs (miRNAs) and experimental evidence for miRNA interactions has ushered in a new area of miRNA research that focuses on network more than individual miRNA interaction, which provides a systematic view of the whole microRNome. So it is a challenge to infer miRNA functional interactions on a system-wide level and further draw a miRNA functional network (miRFN). A few studies have focused on the well-studied human species; however, these methods can neither be extended to other non-model organisms nor take fully into account the information embedded in miRNA-target and target-target interactions. Thus, it is important to develop appropriate methods for inferring the miRNA network of non-model species, such as soybean (Glycine max), without such extensive miRNA-phenotype associated data as miRNA-disease associations in human. Here we propose a new method to measure the functional similarity of miRNAs considering both the site accessibility and the interactive context of target genes in functional gene networks. We further construct the miRFNs of soybean, which is the first study on soybean miRNAs on the network level and the core methods can be easily extended to other species. We found that miRFNs of soybean exhibit a scale-free, small world and modular architecture, with their degrees fit best to power-law and exponential distribution. We also showed that miRNA with high degree tends to interact with those of low degree, which reveals the disassortativity and modularity of miRFNs. Our efforts in this study will be useful to further reveal the soybean miRNA-miRNA and miRNA-gene interactive mechanism on a systematic level. A web tool for information retrieval and analysis of soybean miRFNs and the relevant target functional gene networks can be accessed at SoymiRNet: http://nclab.hit.edu.cn/SoymiRNet.

Analyses of a Glycine max Degradome Library Identify microRNA Targets and MicroRNAs that Trigger Secondary SiRNA Biogenesis

Plant microRNAs (miRNAs) regulate gene expression mainly by guiding cleavage of target mRNAs. In this study, a degradome library constructed from different soybean (Glycine max (L.) Merr.) tissues was deep-sequenced. 428 potential targets of small interfering RNAs and 25 novel miRNA families were identified. A total of 211 potential miRNA targets, including 174 conserved miRNA targets and 37 soybean-specific miRNA targets, were identified. Among them, 121 targets were first discovered in soybean. The signature distribution of soybean primary miRNAs (pri-miRNAs) showed that most pri-miRNAs had the characteristic pattern of Dicer processing. The biogenesis of TAS3 small interfering RNAs (siRNAs) was conserved in soybean, and nine Auxin Response Factors were identified as TAS3 siRNA targets. Twenty-three miRNA targets produced secondary small interfering RNAs (siRNAs) in soybean. These targets were guided by five miRNAs: gma-miR393, gma-miR1508, gma-miR1510, gma-miR1514, and novel-11. Multiple targets of these secondary siRNAs were detected. These 23 miRNA targets may be the putative novel TAS genes in soybean. Global identification of miRNA targets and potential novel TAS genes will contribute to research on the functions of miRNAs in soybean.

Genome-wide identification of soybean microRNAs and their targets reveals their organ-specificity and responses to phosphate starvation

BackgroundPhosphorus (P) plays important roles in plant growth and development. MicroRNAs involved in P signaling have been identified in Arabidopsis and rice, but P-responsive microRNAs and their targets in soybean leaves and roots are poorly understood.ResultsUsing high-throughput sequencing-by-synthesis (SBS) technology, we sequenced four small RNA libraries from leaves and roots grown under phosphate (Pi)-sufficient (+Pi) and Pi-depleted (-Pi) conditions, respectively, and one RNA degradome library from Pi-depleted roots at the genome-wide level. Each library generated ∼21.45−28.63 million short sequences, resulting in ∼20.56−27.08 million clean reads. From those sequences, a total of 126 miRNAs, with 154 gene targets were computationally predicted. This included 92 new miRNA candidates with 20-23 nucleotides that were perfectly matched to the Glycine max genome 1.0, 70 of which belong to 21 miRNA families and the remaining 22 miRNA unassigned into any existing miRNA family in miRBase 18.0. Under both +Pi and -Pi conditions, 112 of 126 total miRNAs (89%) were expressed in both leaves and roots. Under +Pi conditions, 12 leaf- and 2 root-specific miRNAs were detected; while under -Pi conditions, 10 leaf- and 4 root-specific miRNAs were identified. Collectively, 25 miRNAs were induced and 11 miRNAs were repressed by Pi starvation in soybean. Then, stem-loop real-time PCR confirmed expression of four selected P-responsive miRNAs, and RLM-5’ RACE confirmed that a PHO2 and GmPT5, a kelch-domain containing protein, and a Myb transcription factor, respectively are targets of miR399, miR2111, and miR159e-3p. Finally, P-responsive cis-elements in the promoter regions of soybean miRNA genes were analyzed at the genome-wide scale.ConclusionsLeaf- and root-specific miRNAs, and P-responsive miRNAs in soybean were identified genome-wide. A total of 154 target genes of miRNAs were predicted via degradome sequencing and computational analyses. The targets of miR399, miR2111, and miR159e-3p were confirmed. Taken together, our study implies the important roles of miRNAs in P signaling and provides clues for deciphering the functions for microRNA/target modules in soybean.

Detection and evolutionary analysis of soybean miRNAs responsive to soybean mosaic virus

MicroRNAs (miRNA) are a class of non-coding RNAs that have important gene regulatory roles in various organisms. However, the miRNAs involved in soybean's response to soybean mosaic virus (SMV) are unknown. To identify novel miRNAs and biotic-stress regulated small RNAs that are involved in soybean's response to SMV, two small RNA libraries were constructed from mock-inoculated and SMV-infected soybean leaves and sequenced. This led to the discovery of 179 miRNAs, representing 52 families, among which five miRNAs belonging to three families were novel miRNAs in soybean. A large proportion (71.5%) of miRNAs arose from segmental duplication, similar to the process that drives the evolution of protein-coding genes. In addition, we predicted 346 potential targets of these identified miRNAs, and verified 12 targets by modified 5'-RACE analysis. Finally, three miRNAs (miR160, miR393 and miR1510) that are involved in plant resistance were observed to respond to SMV infection. The interaction between miRNAs and resistance-related genes provides a novel mechanism for pathogens to evade host recognition.

Identification of wild soybean miRNAs and their target genes responsive to aluminum stress

BackgroundMicroRNAs (miRNAs) play important regulatory roles in development and stress response in plants. Wild soybean (Glycine soja) has undergone long-term natural selection and may have evolved special mechanisms to survive stress conditions as a result. However, little information about miRNAs especially miRNAs responsive to aluminum (Al) stress is available in wild soybean.ResultsTwo small RNA libraries and two degradome libraries were constructed from the roots of Al-treated and Al-free G. soja seedlings. For miRNA identification, a total of 7,287,655 and 7,035,914 clean reads in Al-treated and Al-free small RNAs libraries, respectively, were generated, and 97 known miRNAs and 31 novel miRNAs were identified. In addition, 49 p3 or p5 strands of known miRNAs were found. Among all the identified miRNAs, the expressions of 30 miRNAs were responsive to Al stress. Through degradome sequencing, 86 genes were identified as targets of the known miRNAs and five genes were found to be the targets of the novel miRNAs obtained in this study. Gene ontology (GO) annotations of target transcripts indicated that 52 target genes cleaved by conserved miRNA families might play roles in the regulation of transcription. Additionally, some genes, such as those for the auxin response factor (ARF), domain-containing disease resistance protein (NB-ARC), leucine-rich repeat and toll/interleukin-1 receptor-like protein (LRR-TIR) domain protein, cation transporting ATPase, Myb transcription factors, and the no apical meristem (NAM) protein, that are known to be responsive to stress, were found to be cleaved under Al stress conditions.ConclusionsA number of miRNAs and their targets were detected in wild soybean. Some of them that were responsive to biotic and abiotic stresses were regulated by Al stress. These findings provide valuable information to understand the function of miRNAs in Al tolerance.