Identification of R2R3-MYB gene family reveal candidate genes for anthocyanin biosynthesis in Lonicera caerulea fruit based on RNA-seq data

BACKGROUND: R2R3-MYB transcription factor (TF) family plays important roles in various biological processes in many plants, especially in the regulation of plant flavonoid accumulation. The fruit of Lonicera caerulea contains abundant anthocyanin. OBJECTIVE: The R2R3-MYB TF family was systematically analyzed according to the RNA-seq data, and the R2R3-MYB candidate genes that were involved in anthocyanin biosynthesis in the fruit of Lonicera caerulea were screened. METHODS: The R2R3-MYB TFs in Lonicera caerulea were identified, and the physical and chemical properties, protein conserved sequence alignment and motifs of each R2R3-MYB TFs were analyzed using bioinformatics methods. The expression levels of these genes and anthocyanin levels in different tissues and different developmental stages of fruit were determined by RT-qPCR and pH shift method. RESULTS: A total of 59 genes encoding R2R3-MYB TFs in Lonicera caerulea were identified and clustered into 20 subgroups (C1 to C20) based on the relationship to AtR2R3-MYBs. Expression profiles showed that the expression of CL6086 and CL552 in fruit were higher than other tissues, and upregulated in the veraison fruit compared to the green ripe fruit. As the expression of the two genes was concurrent with the anthocyanin content, and showed high correlation with anthocyanin biosynthetic structural genes, they were considered as closely related to anthocyanin biosynthesis in the fruit. CONCLUSION: The results provide a systematic analysis of LcR2R3-MYBs, and the foundation for further molecular mechanisms research of anthocyanin biosynthesis regulated by R2R3-MYB in the fruit of Lonicera caerulea.