Comparative transcriptomics and weighted gene co-expression correlation network analysis (WGCNA) reveal potential regulation mechanism of carotenoid accumulation in Chrysanthemum × morifolium

The variation of flower color of chrysanthemum (Chrysanthemum×morifolium) is extremely rich, and carotenoids, which are mainly stored in the plastid, are important pigments that determine the color of chrysanthemum. However, the genetic regulation of the carotenoid metabolism pathway in this species still remains unclear. In this study, a pink chrysanthemum cultivar, 'Jianliuxiang Pink', and its three bud sport mutants (including white, yellow and red color mutants, 'Jianliuxiang White', 'Jianliuxiang Yellow' and 'Jianliuxiang Red', respectively) were used as experimental materials to analyze the dynamic changes of carotenoid components and plastid ultrastructure at different developmental stages of ray florets. We found that the carotenoid components and plastid ultrastructure of the four color cultivars in the early developmental stage of the chrysanthemum capitulum (S1) were almost identical, and the carotenoids mainly included violaxanthin, lutein and β-carotene, which exist in proplastids and immature chloroplasts. With the development of capitulum, the chloroplasts in 'Jianliuxiang White' and 'Jianliuxiang Pink' were degraded, and the protoplasts did not transform but rather formed vesicles that accumulated trace amounts of carotenoids. The proplastids and chloroplasts in 'Jianliuxiang Yellow' and 'Jianliuxiang Red' were all transformed into chromoplasts and consist of lutein as well as lutein's isomer and derivatives. Using comparative transcriptomics combined with gene expression analysis, we found that CmPg-1, CmPAP10, and CmPAP13, which were involved in chromoplast transformation, CmLCYE, which was involved in carotenoid biosynthesis, and CmCCD4a-2, which was involved in carotenoid degradation, were differentially expressed between four cultivars, and these key genes therefore should affect the accumulation of carotenoids in chrysanthemum. In addition, six transcription factors, CmMYB305, CmMYB29, CmRAD3, CmbZIP61, CmAGL24, CmNAC1, were screened using weighted gene co-expression correlation network analysis (WGCNA) combined with correlative analysis to determine whether they play an important role in carotenoid accumulation by regulating structural genes related to the carotenoid metabolism pathway and plastid development. This study analyzed dynamic changes of carotenoid components and plastid ultrastructure of the four bud mutation cultivars of chrysanthemum and identified structural genes and transcription factors that may be involved in carotenoid accumulation. The above results laid a solid foundation for further analysis of the regulatory mechanism of the carotenoid biosynthesis pathway in chrysanthemum.