Aquilegia oxysepala is a characteristic wildflower of Northeast China, which has great application potential owing to its unique floral morphology and strong adaptability. The R2R3-MYB is one of the largest transcription factor families and has important roles in anthocyanin biosynthesis. A total of 90 AoMYBs were identified from A. oxysepala and divided into 29 subfamilies via phylogenetic analysis. Using bioinformatics methods, a systematic analysis was conducted on the physicochemical properties, gene structures, conserved motifs, and other characteristics of these genes. Based on the transcriptomic data of A. oxysepala and A. oxysepala f. pallidiflora across different floral developmental stages, 32 differentially expressed genes (DEGs) belonging to the R2R3-MYB family were screened, and their expression patterns were analyzed. The results indicated that AoMYB32 and AoMYB34 were significantly upregulated only at the pre-anthesis stage (PrA) in A. oxysepala, whereas their expression levels remained relatively low in A. oxysepala f. pallidiflora during the same stage. Correlation network analysis indicated that AoMYB34 was positively correlated with six functional genes associated with anthocyanin synthesis, namely PAL-2, C4H, 4CL, CHS, F3′5′H, and ANS. Reverse transcription quantitative PCR (RT-qPCR) validation revealed that AoMYB32 and AoMYB34 exhibited significant upregulation at the PrA stage in A. oxysepala, while their expression showed no obvious fluctuations across the entire flowering period in A. oxysepala f. pallidiflora. In summary, it is postulated that AoMYB32 and AoMYB34 are candidate genes associated with flower color variation formation of A. oxysepala. This study provides a theoretical basis and genetic resources for deciphering the mechanism of flower color formation and facilitating flower color improvement breeding in A. oxysepala.
Ma et al. (Sun,) studied this question.