De novo genome assembly of faba bean (Vicia faba L.) leaf transcriptome provide novel insights into the drought stress tolerance mechanisms

Drought is a major limiting constraint to faba bean production worldwide, including Tunisia. However, molecular mechanisms underlying faba bean responses to drought stress are not well understood. In this context, transcriptome analysis by RNA-seq was performed to investigate drought-related genes and construct a network of faba bean drought stress response and tolerance. De novo assembly of the transcriptome generated 26,728 differentially expressed genes (DEGs). Of these, 13,920 were up-regulated and 12,808 down-regulated in faba bean drought-stressed leaves. Moreover, a total of 10,800 simple sequence repeats (SSRs) and 2130 transcription factors involved in major metabolic pathways including abscisic acid (ABA)-dependent and -independent signaling pathway were identified. GO, KOG and KEGG enrichment analyses revealed that these DEGs were involved in several important processes including photosynthesis, flavonoid biosynthesis, response to stimulus and abiotic stress, reactive oxygen species (ROS) scavengers, signal transduction, biosynthesis of secondary metabolites and transporters, suggesting the involvement of these important pathways in faba bean response to water deficit. Various stress proteins such as late embryogenesis abundant proteins (LEA), dehydrins (DHNs) and heat shock proteins (HSPs) have been identified and their expression was robustly upregulated in drought-stressed leaves, indicating their key contribution to drought response and adaptation by conferring protection and providing stability to faba bean plant cellular processes under water deficit. The reliability of the RNA-seq results was confirmed by the analysis of 10 randomly selected genes using qRT-PCR. Taken together, these findings help advancing our knowledge and can guide breeding programs aimed at improving the tolerance of faba bean to drought stress. The goal of the study was to investigate the leaf transcriptomic changes of faba bean (Vicia faba) in response to drought stress and to perform de novo assembly. The collected leaf samples (stressed and control) were ground in liquid nitrogen and the resulting powder was then used for RNA extraction using NucleoSpin RNA Plant kit (Macherey-Nagel, Germany). Library construction and RNA sequencing were carried out at Beijing Genomics Institute (BGI), China, at a sequencing depth of 40 million paired end reads.