Genome-Wide Identification of bZIP Transcription Factors in Faba Bean Based on Transcriptome Analysis and Investigation of Their Function in Drought Response

Faba bean is an important cool-season edible legume crop that is constantly threatened by abiotic stresses such as drought. The basic leucine zipper (bZIP) gene family is one of the most abundant and diverse families of transcription factors in plants. It regulates plant growth and development and plays an important role in the response to biotic and abiotic stresses. In this study, we identified 18 members of the faba bean bZIP transcription factor family at the genome-wide level based on previous faba bean drought stress transcriptome sequencing data. A phylogenetic tree was constructed to group the 18 VfbZIP proteins into eight clades. Analysis of cis-acting elements in the promoter region suggested that these 18 VfbZIPs may be involved in regulating abiotic stress responses such as drought. Transcriptome data showed high expression of seven genes (VfbZIP1, VfbZIP2, VfbZIP5, VfbZIP7, VfbZIP15, VfbZIP17, and VfbZIP18) in the drought-tolerant cultivar under drought stress, in which VfbZIP1, VfbZIP2, and VfbZIP5 were consistently expressed as detected by quantitative real-time polymerase chain reaction (qRT-PCR) compared to the transcriptome data. Ectopic overexpression of the three VfbZIPs in tobacco, based on the potato Virus X (PVX) vector, revealed that VfbZIP5 enhanced the drought tolerance. Overexpressed VfbZIP5 in plants showed lower levels of proline (PRO), malondialdehyde (MDA), and peroxidase (POD) compared to those overexpressing an empty vector under 10 days of drought stress. Protein-protein interaction (PPI) analysis showed that VfbZIP5 interacted with seven proteins in faba bean, including VfbZIP7 and VfbZIP10. The results depict the importance of VfbZIPs in response to drought stress, and they would be useful for the improvement of drought tolerance.