Genome-Wide Comparative Analysis of the R2R3-MYB Gene Family in Five Solanaceae Species and Identification of Members Regulating Carotenoid Biosynthesis in Wolfberry

The R2R3-MYB is a large gene family involved in various plant functions, including carotenoid biosynthesis. However, this gene family lacks a comprehensive analysis in wolfberry (Lycium barbarum L.) and other Solanaceae species. The recent sequencing of the wolfberry genome provides an opportunity for investigating the organization and evolutionary characteristics of R2R3-MYB genes in wolfberry and other Solanaceae species. A total of 610 R2R3-MYB genes were identified in five Solanaceae species, including 137 in wolfberry. The LbaR2R3-MYB genes were grouped into 31 subgroups based on phylogenetic analysis, conserved gene structures, and motif composition. Five groups only of Solanaceae R2R3-MYB genes were functionally divergent during evolution. Dispersed and whole duplication events are critical for expanding the R2R3-MYB gene family. There were 287 orthologous gene pairs between wolfberry and the other four selected Solanaceae species. RNA-seq analysis identified the expression level of LbaR2R3-MYB differential gene expression (DEGs) and carotenoid biosynthesis genes (CBGs) in fruit development stages. The highly expressed LbaR2R3-MYB genes are co-expressed with CBGs during fruit development. A quantitative Real-Time (qRT)-PCR verified seven selected candidate genes. Thus, Lba11g0183 and Lba02g01219 are candidate genes regulating carotenoid biosynthesis in wolfberry. This study elucidates the evolution and function of R2R3-MYB genes in wolfberry and the four Solanaceae species.