Water deficit has a significant impact on growth, development and yield of fava bean (Vicia fava L.) in arid and semi-arid climates. The aim of this study was to identify differentially expressed genes in the Qinghai 13 genotype under soil drought through leaf transcriptome analysis. A total of 256.95 M clean reads were obtained and assembled into 176334 unigenes, with an average length of 766 bp. A total of 9126 (4439 upregulated and 4687 downregulated) differentially expressed genes (DEGs) were identified in faba bean leaves under soil drought. In total, 324 putative transcription factors were identified and classified as belonging to different transcription factor families. According to GO and KEGG analysis, the soil drought stress-inducible DEGs encoded proteins mainly involved in regulating photosynthesis, osmotic adjustment, detoxification, autophagy and other functions. In addition, a large portion of DEGs appeared to be novel because they could not be annotated in any functional databases, therefore, suggesting a specific response to soil drought in faba bean. Finally, RNA-seq analysis was validated by quantitative reverse-transcription PCR analysis. This work provides comprehensive and valuable information for understanding the molecular mechanisms which faba bean uses to respond to soil drought.

Keywords: Differentially expressed genes; Faba bean; Function; Soil drought; Transcriptome.