The MYB gene family constitutes one of the largest transcription factors (TFs) modulating various biological processes in plants. Although genome-wide analysis of this gene family has been carried out in some species, only three MYB members have been functionally characterized heretofore in sesame (Sesamum indicum L.). Here, we identified a relatively high number (287) of sesame MYB genes (SIMYBs) with an uncommon overrepresentation of the 1R-subfamily. A total of 95% of SIMYBs was mapped unevenly onto the 16 linkage groups of the sesame genome with 55 SIMYBs tandemly duplicated. In addition, molecular characterization, gene structure, and evolutionary relationships of SIMYBs were established. Based on the close relationship between sesame and Arabidopsis thaliana, we uncovered that the functions of SIMYBs are highly diverse. A total of 65% of SIMYBs were commonly detected in five tissues, suggesting that they represent key TFs modulating sesame growth and development. Moreover, we found that SIMYBs regulate sesame responses to drought and waterlogging, which highlights the potential of SIMYBs towards improving stress tolerance in sesame. This work presents a comprehensive picture of the MYB gene family in sesame and paves the way for further functional validation of the members of this versatile gene family.

Sesame is a source of a healthy vegetable oil, attracting a growing interest worldwide. Abiotic stresses have devastating effects on sesame yield; hence, studies have been performed to understand sesame molecular responses to abiotic stresses, but the core abiotic stress-responsive genes (CARG) that the plant reuses in response to an array of environmental stresses are unknown. We performed a meta-analysis of 72 RNA-Seq datasets from drought, waterlogging, salt and osmotic stresses and identified 543 genes constantly and differentially expressed in response to all stresses, representing the sesame CARG. Weighted gene co-expression network analysis of the CARG revealed three functional modules controlled by key transcription factors. Except for salt stress, the modules were positively correlated with the abiotic stresses. Network topology of the modules showed several hub genes predicted to play prominent functions. As proof of concept, we generated over-expressing Arabidopsis lines with hub and non-hub genes. Transgenic plants performed better under drought, waterlogging, and osmotic stresses than the wild-type plants but did not tolerate the salt treatment. As expected, the hub gene was significantly more potent than the non-hub gene. Overall, we discovered several novel candidate genes, which will fuel investigations on plant responses to multiple abiotic stresses.