The MYB transcription factor family has numerous members, and is involved in biological activities, such as ABA signaling, which plays an important role in a plant's resistance to abiotic stresses such as drought. However, the diversity of MYB members that respond to drought stress and their regulatory mechanisms in different flax varieties were unclear. In this study, we obtained 855.69 Gb of clean data from 120 flax root samples from 20 flax (Linum usitatissimum L.) varieties, assembled 92,861 transcripts, and identified 434 MYB family members in each variety. The expression profiles of the MYB transcription factor family from 20 flax varieties under drought stress were analyzed. The results indicated that there are four strategies by which the MYB family responds to drought stress in these 20 flax varieties, each of which has its own specific processes, such as development, reproduction, and localization processes. The four strategies also include common biological processes, such as stimulus responses, metabolic processes, and biological regulation. The WGCNA method was subsequently employed to identify key members of the MYB family involved in response strategies to drought stress. The results demonstrated that a 1R-MYB subfamily gene co-expression network is significantly related to the gibberellin response and cytokinin-activated signaling pathway processes in the 'Strategy 4' for MYB family response to drought, identifying core genes such as Lus.scaffold70.240. Our results showed a diversity of MYB family responses to drought stress within flax varieties, and these results contribute to deciphering the mechanisms of the MYB family regulation of drought resistance. This will promote the more accurate breeding development of flax to adapt to agricultural production under drought conditions.
Keywords: WGCNA; expression profile; functional annotation; heat map; transcription factor.