Transcriptome analysis and identification of genes associated with flower development in Rhododendron pulchrum Sweet (Ericaceae)

Flowering process is essential for plant development. However, the molecular mechanisms driving flower development of ornamental woody Rhododendron pulchrum Sweet are difficult to elucidate due to the lack of genomic data. In this research, high-throughput sequencing and comparative transcriptome analyses of R. pulchrum flowers collected at three key stages were performed: floral bud stage, early flowering stage, and full-flowering stage. Furthermore, expression of genes involved in flower development was also validated with quantitative real-time PCR (qRT-PCR). RNA-seq yielded 96,350,697 bp of clean reads, which were assembled into 98,610 unigenes with an average length of 717 bp. 58,279 (59.10%) unigenes could be annotated, including 324 major unigenes associated with floral development. In addition, ten modules (20,443 mRNAs) were dissected in the co-expression network. Especially, Flowering Locus (FLC) and Flowering Locus T (FT) were co-expressed. 9493 differentially expressed genes (DEGs) were scanned among three stages, and most DEGs existed between flower bud stage and early flowering stage. In particular, 79 DGEs associated with flowering process were enriched in 28 GO terms. Moreover, the expression levels of MYC2, EIN3, and ARR-B were all lowest at early flowering stage, while transcripts of MYC2, TIR1, CYCD3, COL-1, and EIN3 were all peaked at flower bud stage. Transcriptome profile presented here will benefit deep insights into molecular mechanism underlying R. pulchrum flowering process.