Abstract
Rafflesia possesses unique biological features and known primarily for producing the world's largest and existing as a single flower. However, to date, little is known about key regulators participating in Rafflesia flower development. In order to further understand the molecular mechanism that regulates Rafflesia cantleyi flower development, RNA-seq data from three developmental stages of floral bud, representing the floral organ primordia initiation, floral organ differentiation, and floral bud outgrowth, were analysed. A total of 89,890 transcripts were assembled of which up to 35% could be annotated based on homology search. Advanced transcriptome analysis using K-mean clustering on the differentially expressed genes (DEGs) was able to identify 12 expression clusters that reflect major trends and key transitional states, which correlate to specific developmental stages. Through this, comparative gene expression analysis of different floral bud stages identified various transcription factors related to flower development. The members of WRKY, NAC, bHLH, and MYB families are the most represented among the DEGs, suggesting their important function in flower development. Furthermore, pathway enrichment analysis also revealed DEGs that are involved in various phytohormone signal transduction events such as auxin and auxin transport, cytokinin and gibberellin biosynthesis. Results of this study imply that transcription factors and phytohormone signalling pathways play major role in Rafflesia floral bud development. This study provides an invaluable resource for molecular studies of the flower development process in Rafflesia and other plant species.