Genome-Wide Analysis of Long Non-Coding RNAs (lncRNAs) in Two Contrasting Soybean Genotypes Subjected to Phosphate Starvation
Abstract Background: Phosphorus (P) is essential for plant growth and development, and low-phosphorus (LP) stress is a major factor limiting the growth and yield of soybean. Recently, long noncoding RNAs (lncRNAs) have been reported to be key regulators in response to stress conditions in plants. In soybean, however, how LP stress mediates biogenesis of lncRNAs remains unclear. Results: In this study, to explore the response mechanisms of lncRNAs to LP stress, we used the roots of two representative soybean genotypes with opposite P deficiency responsiveness, a P-sensitive genotype (Bogao) and a P-tolerant genotype (NN94156), to construct RNA sequencing (RNA-seq) libraries. In total, 4,166 novel lncRNAs including 525 differently expressed (DE) lncRNAs were identified across two genotypes at different P levels. GO and KEGG analyses indicated that numerous DE lncRNAs might be involved in diverse biological processes of phosphate, such as lipid metabolic process, catalytic activity, cell membrane formation, signal transduction, nitrogen fixation. Moreover, lncRNA-mRNA-miRNA and lncRNA-mRNA networks were constructed and several promising lncRNAs were identified, which may have highly valuable for further analysis the mechanism in response to LP stress in soybean.Conclusions: These results revealed that LP stress can significantly alter the genome-wide profiles of lncRNAs, especially for P sensitive genotype Bogao. Our findings increase the understanding and provides new insights into the function of lncNRAs responses to P stress in soybean.