*2.1. Identification of MYB Gene Members in the H. hamabo Genome*

HMMER software (v. 3.1) was used with the hidden Markov model (HMM) of the MYB domain (PF00249), obtained from the Pfam database (http://pfam.xfam.org/, accessed on 18 March 2021), to identify MYB genes in the *H. hamabo* genome (Bioproject accession number: PRJNA759075) with E-values < 0.001 [2,4]. ClustalX software (v2.0) was used for alignment of protein sequences and eliminating redundant sequences. The Pfam website (http://pfam.xfam.org/search#tabview=tab1, accessed on 18 March 2021) was used to check the MYB domain structure, and sequences with a typical conserved domain were identified as *H. hamabo* MYB members (*HhMYB*), which were used in further analyses.
