Genome-Wide Identification of microRNAs in Response to Salt/Alkali Stress in Medicago truncatula through High-Throughput Sequencing

Saline-alkaline stress is a universal abiotic stress that adversely affects plant growth and productivity. Saline-alkaline conditions results in plant abnormal transcriptome expression finally manifesting as defective phenotypes. Considerable research has revealed the active role of microRNA in various stress conditions. This study was aimed to identify novel miRNAs and the miRNA expression patterns in the leguminous model plant R108 (Medicago truncatula). The miRNA contained in the total RNA extracted from Medicago truncatula seedlings (72 h) that had been treated with solutions mimicking saline and alkaline soils was subjected to miRNA deep sequencing. The Illumina HiSeq sequencing platform was used to analyze nine small RNA libraries of three treatment groups: distilled water, 20 mM NaCl + Na2SO4 and 5 mM Na2CO3 + NaHCO3. Sequencing revealed that 876 miRNAs including 664 known miRNAs and 212 potential novel miRNAs were present in all the libraries. The miR159 family, miR156 family, miR2086-3p, miR396, miR166, miR319, miR167, miR5213-5p, miR1510 and miR2643 were among the most expressed miRNAs in all libraries. The results of miRNAs expression under treatments were validated by reverse-transcription quantitative PCR (RT-qPCR). Target gene prediction through computational analysis and pathway annotation analysis revealed that the primary pathways affected by stress were related to plant development, including metabolic processes, single-organism processes and response to the stimulus. Our results provide valuable information towards elucidating the molecular mechanisms of salt/alkali tolerance in Medicago truncatula and provide insight into the putative role of miRNAs in plant stress resistance.