Genome-Wide Identification and Evolution Analysis of R2R3-MYB Gene Family Reveals S6 Subfamily R2R3-MYB Transcription Factors Involved in Anthocyanin Biosynthesis in Carrot

The taproot of purple carrot accumulated rich anthocyanin, but non-purple carrot did not. MYB transcription factors (TFs) condition anthocyanin biosynthesis in many plants. Currently, genome-wide identification and evolution analysis of R2R3-MYB gene family and their roles involved in conditioning anthocyanin biosynthesis in carrot is still limited. In this study, a total of 146 carrot R2R3-MYB TFs were identified based on the carrot transcriptome and genome database and were classified into 19 subfamilies on the basis of R2R3-MYB domain. These R2R3-MYB genes were unevenly distributed among nine chromosomes, and Ka/Ks analysis suggested that they evolved under a purified selection. The anthocyanin-related S6 subfamily, which contains 7 MYB TFs, was isolated from R2R3-MYB TFs. The anthocyanin content of rhizodermis, cortex, and secondary phloem in 'Black nebula' cultivar reached the highest among the 3 solid purple carrot cultivars at 110 days after sowing, which was approximately 4.20- and 3.72-fold higher than that in the 'Deep purple' and 'Ziwei' cultivars, respectively. The expression level of 7 MYB genes in purple carrot was higher than that in non-purple carrot. Among them, DcMYB113 (DCAR_008994) was specifically expressed in rhizodermis, cortex, and secondary phloem tissues of 'Purple haze' cultivar, with the highest expression level of 10,223.77 compared with the control 'DPP' cultivar at 70 days after sowing. DcMYB7 (DCAR_010745) was detected in purple root tissue of 'DPP' cultivar and its expression level in rhizodermis, cortex, and secondary phloem was 3.23-fold higher than that of secondary xylem at 110 days after sowing. Our results should be useful for determining the precise role of S6 subfamily R2R3-MYB TFs participating in anthocyanin biosynthesis in carrot.