Genome-wide analysis of the MYB gene family and functional analysis of BhMYB79 in wax gourd

Wax gourd ( Benincasa hispida ) is an important cucurbit crop with economic and medicinal value. The myeloblastosis (MYB) gene family is one of the largest gene families in plants and regulates various biological processes, whereas the MYB gene family has not been systematically studied in wax gourd. In this study, we performed genome-wide identification of the MYB gene family in wax gourd and analyzed their phylogenetic relationship, MYB DNA-binding domain (MYB DBD), gene structure, protein motif, synteny, duplication mode and expression pattern. As a result, a total of 215 BhMYB genes (BhMYBs) were identified, belonging to four subfamilies: 1R-, 2R-, 3R- and 4R-MYB subfamilies. Genes of 1R-MYB subfamily and 2R-MYB subfamily were subdivided into different subgroups respectively. The analysis of MYB DBD, gene structure and protein motif showed that the most genes in the same subgroup had similar characteristics and the 2R-MYB genes were more conserved than the 1R-MYB genes. Interestingly, the long terminal retrotransposons (LTR-RTs) were found in the long introns of several BhMYBs. The results of synteny analysis showed that there were more syntenic gene pairs between wax gourd and other cucurbit crops, while the least number of syntenic gene pairs existed between wax gourd and rice. Gene duplication was the main reason for the expansion of the MYB gene family in wax gourd, with the transposed duplication (TRD) mode contributing more. All duplication BhMYB genes were under purifying selection pressure. Further expression analysis showed that many BhMYBs exhibited obvious tissue-specific expression and several BhMYBs were significantly induced by one or more abiotic stresses. BhMYB79 was particularly expressed in roots and significantly induced by salt, drought, cold and heat stresses, overexpression of which led to reduced tolerance to salt stress in Arabidopsis . In conclusion, our results provide a systematic analysis of wax gourd MYB gene family and facilitate the biological role study of BhMYB79 during wax gourd salt stress response process.