Plant MYB transcription factors (TFs) play crucial roles in regulating the biosynthesis of flavonoids but current analysis on their role in Zanthoxylum bungeanum Maxim. (ZBM) is far from comprehensive. In this study, we identified 270 MYB genes in ZBM and divided them into four subfamilies. The R2R3-MYB (ZbMYB) category contained 251 genes and was classified into 33 subfamilies according to their phylogenetic results and sequence similarity. These subfamilies included 24 subgroups containing both MYBs of ZBM plants and AtMYBs, and nine subgroups containing only ZBM MYBs or AtMYBs. ZbMYBs with similar functions clustered into the same subgroup, indicating functional conservation. The subcellular localization analysis predicted that most ZbMYB genes were found in the nucleus. The transposed duplications appeared to play a major role in the expansion of the MYB gene family in ZBM. Through phylogenetic analysis and transcriptome profiling, it was found that 28 ZbMYB genes may regulate the biosynthesis of flavonoids in ZBM, and these genes expression presented distinct temporal and spatial expression patterns. In different fruit development stages of ZBM, the expression patterns of EVM0042160 and EVM0033809 genes obtained by qRT-PCR analysis are very similar to the flavonoid and anthocyanin content curves in ZBM. Further correlation analysis showed that the content of flavonoids in different fruit development stages and the transcript abundance levels of 28 ZbMYB genes have different degrees of correlation relationship. These results indicated that the ZbMYB genes might be involved in the flavonoid metabolic pathway. This comprehensive and systematic analysis of MYB family genes provided a solid foundation for further functional analysis of MYB TFs in ZBM.