Transcriptomic Analysis Reveals Transcription Factors Related to Leaf Anthocyanin Biosynthesis in Paeonia qiui

Paeonia qiui is a wild species of tree peony. P. qiui has good ornamental value owing to its leaf color change in spring. So far, the molecular mechanism of leaf color change in P. qiui is unclear. This study analyzes the anthocyanin level and transcriptome of two different color stages in P. qiui leaf. The purplish-red leaf stage is rich in anthocyanin, while the green leaf has very low levels of anthocyanin. Transcriptome analysis reveals that 6678 differentially-expressed genes (DEGs) are up-regulated, and 14,667 are down-regulated in the purplish-red leaf. Among these DEGs, 40 MYB (v-myb avian myeloblastosis viral oncogene homolog) genes, 40 bHLH (MYC-like basic helix–loop–helix) genes, and 15 WD40 (WD-repeat protein) genes were found. Based on phylogenetic and alignment analysis with the deduced amino acid sequences with known transcription factors, Unigene0024459 (MYB1) is likely the R2R3-MYB that promotes anthocyanin biosynthesis; Unigene0050761 (MYB2) is likely the R2R3-MYB that represses anthocyanin biosynthesis; Unigene0005081 (bHLH1) and Unigene0006146 (WD40-1) are likely the bHLH and WD40 that participate in regulating anthocyanin biosynthesis. Additionally, quantitative RT-PCR results confirmed the transcriptome analyses for key genes