ABSTRACT
Pummelo (Citrus maxima) is one of important fruit trees, which belongs to Citrus species. The fruits of different pummelo cultivars have different colors and differ in the contents of carotenoid. Our results clearly showed that 'Huangjinmiyou' (HJMY) has the highest content of ß-carotene, followed by 'Hongroumiyou' (HRMY) and 'Guanximiyou' (GXMY). Lycopene is dominantly accumulated in HRMY. However, the molecular mechanism underlying the carotenoid accumulation in pummelo flesh is not fully understood. In this study, we used the RNA-Seq technique to investigate the candidate genes of carotenoid metabolism in the flesh of pummelo cv. GXMY and its mutants HRMY and HJMY in three development periods of fruit. After data assembly and bioinformatic analysis, a total of 357 genes involved in biosynthesis of secondary metabolites were isolated, of which 12 differentially expressed genes (DEGs) are involved in carotenoid biosynthesis. Among these 12 DEGs, phytoene synthase (PSY2), lycopene ß-cyclase (LYCB2), lycopene Æ-cyclase (LYCE), carotenoid cleavage dioxygenases (CCD4), 9-cis-epoxycarotenoid dioxygenase (NCED2), aldehyde oxidase 3 (AAO3), and ABA 8'-hydroxylases (CYP707A1) are the most distinct DEGs in three pummelo cultivars. The co-expression analysis revealed that the expression patterns of several transcription factors such as bHLH, MYB, ERF, NAC and WRKY are highly correlated with DEGs, which are involved in carotenoid biosynthesis. In addition, the expression patterns of 22 DEGs were validated by real-time quantitative PCR (RT-qPCR) and the results are highly concordant with the RNA-Seq results. Our results provide a global vision of transcriptomic profile among three pummelo cultivars with different pulp colors. These results would be beneficial to further study the molecular mechanism of carotenoid accumulation in pummelo flesh and help the breeding of citrus with high carotenoid content.