Rhododendron pulchrum Sweet is a traditional ornamental plant cultivated in China and presents a great variation in petal coloration. However, few studies have been performed to reveal the genes involved and the regulatory mechanism of flower color formation in this plant. In this study, to explore the underlying genetic basis of flower color formation, transcriptome analysis was performed by high-throughput sequencing techniques on four petal samples of different colors: purple, pink, light pink, and white. Results show that a total of 35.55 to 40.56 million high-quality clean reads were obtained, of which 28.56 to 32.65 million reads were mapped to the reference genome. For their annotation, 28,273, 18,054, 24,301, 19,099, and 11,507 genes were allocated to Nr, Swiss-Prot, Pfam, GO, and KEGG databases, correspondingly. There were differentially expressed genes among the four different petal samples, including signal-transduction-related genes, anthocyanin biosynthesis genes, and transcription factors. We found that the higher expressed levels of genes associated with flavonol synthase (FLS) might be the key to white formation, and the formation of red color may be related to the higher expression of flavanone 4-reductase (DFR) families. Overall, our study provides some valuable information for exploring and understanding the flower color intensity variation in R. pulchrum.
