*3.1. Identification of MYB Gene Family Members in H. hamabo*

In order to identify MYB family genes in the genome of *H. hamabo*, we used the hidden Markov model profile (HMM) of the MYB domain (PF00249) to conduct an HMM-search based against the predicted gene coding sequences. A total of 204 genes were obtained with E-values < 0.001 and labeled as *HhMYB1*-*HhMYB204* (Table S2). The length range of the HhMYB sequences was from 300 aa to 1731 aa, of which HhMYB186 had the longest amino acid sequence. In addition, the molecular weights of HhMYBs ranged from 32.81 kDa (HhMYB4) to 196.98 kDa (HhMYB186), and the range of theoretical isoelectric points was 4.51 (HhMYB195) to 10.01 (HhMYB30) (Table S2).
