*4.8. Genes Related to the ABP and Phylogenetic Analyses*

The ABP-related structure genes and transcription factors including *CHS*, *CHI*, *F3H*, *F3 5 H*, *DFR*, *FLS*, *ANS*, *FNS*, *UFGT*, *MYB*, *bHLH*, and *WD40* were used as queries to retrieve the corresponding unigenes from the *P. limprichtii* libraries. Meanwhile, every reference gene from model plants obtained from KEGG were aligned with the libraries using BLASTx to search for more related unigenes.

For the MYB genes, to determine which of the genes belonged to which R2R3-MYB subfamily, the R2R3-MYB genes from *Arabidopsis thaliana* were used to conduct phylogenetic analyses (Amino acid sequence obtained from Gene Bank, accession number is listed in Supplementary Table S2), and MEGA7.0 software (Institute of Molecular Evolutionary Genetics, PA, USA) [52] was used to perform sequence aligning and construct a circular phylogenetic tree according to neighbor joining method with 1000 interactions.
