MYB transcription factors (TFs) are one of the largest TF families in plants to regulate numerous biological processes. However, our knowledge of the MYB family in Physcomitrella patens is limited. We identified 116 MYB genes in the P. patens genome, which were classified into the R2R3-MYB, R1R2R3-MYB, 4R-MYB, and MYB-related subfamilies. Most R2R3 genes contain 3 exons and 2 introns, whereas R1R2R3 MYB genes contain 10 exons and 9 introns. N3R-MYB (novel 3RMYB) and NR-MYBs (novel RMYBs) with complicated gene structures appear to be novel MYB proteins. In addition, we found that the diversity of the MYB domain was mainly contributed by domain shuffling and gene duplication. RNA-seq analysis suggested that MYBs exhibited differential expression to heat and might play important roles in heat stress responses, whereas CCA1-like MYB genes might confer greater flexibility to the circadian clock. Some R2R3-MYB and CCA1-like MYB genes are preferentially expressed in the archegonium and during the transition from the chloronema to caulonema stage, suggesting their roles in development. Compared with that of algae, the numbers of MYBs have significantly increased, thus our study lays the foundation for further exploring the potential roles of MYBs in the transition from aquatic to terrestrial environments.