Abstract
The growth, yield, and seed quality of rapeseed are negatively affected by drought stress. Therefore, it is of great value to understand the molecular mechanism behind this phenomenon. In a previous study, long non-coding RNAs (lncRNAs) were found to play a key role in the response of rapeseed seedlings to drought stress. However, many questions remained unanswered. This study was the first to investigate the expression profile of lncRNAs not only under control and drought treatment, but also under the rehydration treatment. A total of 381 differentially expressed lncRNA and 10,253 differentially expressed mRNAs were identified in the comparison between drought stress and control condition. In the transition from drought stress to rehydration, 477 differentially expressed lncRNAs and 12,543 differentially expressed mRNAs were detected. After identifying the differentially expressed (DE) lncRNAs, the comprehensive lncRNAs-engaged network with the co-expressed mRNAs in leaves under control, drought and rehydration was investigated. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of co-expressed mRNAs identified the most significant pathways related with plant hormones (expecially abscisic acid, auxin, cytokinins, and gibberellins) in the signal transduction. The genes, co-expressed with the most-enriched DE-lncRNAs, were considered as the most effective candidates in the water-loss and water-recovery processes, including protein phosphatase 2 C (PP2C), ABRE-binding factors (ABFs), and SMALL AUXIN UP-REGULATED RNAs (SAURs). In summary, these analyses clearly demonstrated that DE-lncRNAs can act as a regulatory hub in plant-water interaction by controlling phytohormone signaling pathways and provided an alternative way to explore the complex mechanisms of drought tolerance in rapeseed.