Transcriptome landscape of Rafflesia cantleyi floral buds reveals insights into the roles of transcription factors and phytohormones in flower development

Bud development, flower phenology and life history of holoparasitic Rafflesia cantleyi

Despite being the world's largest single-flower, Rafflesia's biology and life history are still poorly understood due to its cryptic growth strategy on Tetrastigma vines. Previous studies have been mostly short-term, contrary to Rafflesia's long development period before blooming. Bud development and flower phenology of R. cantleyi was studied in a dipterocarp forest in Lata Jarum, Peninsular Malaysia. Seven populations, consisting of 247 buds, were monitored fortnightly for 65 months in two discrete studies between 2009 and 2018. The bud size distribution of R. cantleyi is dynamic, progressively changing from small flower buds to larger buds before flowering. Buds < 5.0 cm across had the slowest growth and highest mortality rates, while those > 15.0 cm across demonstrated accelerated growth. The bud growth profiles of the same site clustered distinctively regardless of sex with successful blooming rate that varied greatly between sites, prompting speculation about their relatedness to the sites' physical attributes. We reported the first female-dominated population in Rafflesia's life history. Rafflesia cantleyi's blooming rate at Lata Jarum is moderate to high, with non-seasonal flowering phenology as evident by the lack of synchronisation and consistency between flowering and local rainfall patterns. Based on the field data of the present study and the published information of other Rafflesia species, R. cantleyi's life cycle was estimated between 4.0 and 5.3 years. Our findings further explain Rafflesia's biology and life history and highlight the gap in knowledge of the natural habitats on the endoparasite's growth and fate potentially for future conservation and study.

OfBFT genes play an essential role in the proliferate flower formation of Osmanthus fragrans

Floral organ development is one of the most vital events in flowering plants and is closely related to ornamental properties. The proliferate flower (a new branch or flower occurring in the centre of a flower) in plants is an interesting type, while the specific molecular mechanism remains largely unknown. Osmanthus fragrans 'Tianxiang Taige' has two different flower morphologies: normal flower and proliferate flower. Phenotypic observation suggested that a normal flower was composed of calyx, petal, stamen and pistil (reduced to leaf-like carpel). While in proliferate flower, the leaf-like carpel continued to grow and was replaced by a new branch. Paraffin section indicated that the re-growth of leaf carpels might be the main reason for proliferate flower formation. Transcriptome sequencing of normal and proliferate flower was performed, and the expression levels of related genes were analysed. Among the differentially expressed genes, OfBFT-a and OfBFT-b had differential expression during the proliferate flower formation process. The expression patterns revealed that both OfBFT-a and OfBFT-b were highly accumulated in carpels, and were significantly downregulated during the proliferate flower development process. Subcellular localization indicated that OfBFT-a and OfBFT-b proteins were located in the nucleus. Functional studies in 'Tianxiang Taige' and Arabidopsis showed that OfBFT-a and OfBFT-b had important roles in floral organ development, especially the proliferate flower formation process by downregulating the accumulation of AG and SEP3 homologous genes. These results may shed new light on the study of proliferate flower formation and flower morphology breeding in flowering plants.

Cross-Talk between Transcriptome Analysis and Dynamic Changes of Carbohydrates Identifies Stage-Specific Genes during the Flower Bud Differentiation Process of Chinese Cherry (Prunus pseudocerasus L.)

Flower bud differentiation is crucial to reproductive success in plants. In the present study, RNA-Seq and nutrients quantification were used to identify the stage-specific genes for flower bud differentiation with buds which characterize the marked change during flower bud formation from a widely grown Chinese cherry (Prunus pseudocerasus L.) cultivar 'Manaohong'. A KEGG enrichment analysis revealed that the sugar metabolism pathways dynamically changed. The gradually decreasing trend in the contents of total sugar, soluble sugar and protein implies that the differentiation was an energy-consuming process. Changes in the contents of D-glucose and sorbitol were conformed with the gene expression trends of bglX and SORD, respectively, which at least partially reflects a key role of the two substances in the transition from physiological to morphological differentiation. Further, the WRKY and SBP families were also significantly differentially expressed during the vegetative-to-reproductive transition. In addition, floral meristem identity genes, e.g., AP1, AP3, PI, AGL6, SEP1, LFY, and UFO demonstrate involvement in the specification of the petal and stamen primordia, and FPF1 might promote the onset of morphological differentiation. Conclusively, the available evidence justifies the involvement of sugar metabolism in the flower bud differentiation of Chinese cherry, and the uncovered candidate genes are beneficial to further elucidate flower bud differentiation in cherries.

Transcriptomic Analysis of Sex-Associated DEGs in Female and Male Flowers of Kiwifruit (Actinidia deliciosa [A. Chev] C. F. Liang & A. R. Ferguson)

Kiwifruit (Actinidia deliciosa [A. Chev.], C.V. Liang & A. R. Ferguson, 1984) is a perennial plant, with morphologically hermaphroditic and functionally dioecious flowers. Fruits of this species are berries of great commercial and nutritional importance. Nevertheless, few studies have analyzed the molecular mechanisms involved in sexual differentiation in this species. To determine these mechanisms, we performed RNA-seq in floral tissue at stage 60 on the BBCH scale in cultivar ‘Hayward’ (H, female) and a seedling from ‘Green Light’ × ‘Tomuri’ (G × T, male). From these analyses, we obtained expression profiles of 24,888 (H) and 27,027 (G × T) genes, of which 6413 showed differential transcript abundance. Genetic ontology (GO) and KEGG analysis revealed activation of pathways associated with the translation of hormonal signals, plant-pathogen interaction, metabolism of hormones, sugars, and nucleotides. The analysis of the protein-protein interaction network showed that the genes ERL1, AG, AGL8, LFY, WUS, AP2, WRKY, and CO, are crucial elements in the regulation of the hormonal response for the formation and development of anatomical reproductive structures and gametophytes. On the other hand, genes encoding four Putative S-adenosyl-L-methionine-dependent methyltransferases (Achn201401, Achn281971, Achn047771 and Achn231981) were identified, which were up-regulated mainly in the male flowers. Moreover, the expression profiles of 15 selected genes through RT-qPCR were consistent with the results of RNA-seq. Finally, this work provides gene expression-based interactions between transcription factors and effector genes from hormonal signaling pathways, development of floral organs, biological and metabolic processes or even epigenetic mechanisms which could be involved in the kiwi sex-determination. Thus, in order to decode the nature of these interactions, it could be helpful to propose new models of flower development and sex determination in the Actinidia genus.

Transcriptome analysis of Rafflesia cantleyi flower stages reveals insights into the regulation of senescence

Rafflesia is a unique plant species existing as a single flower and produces the largest flower in the world. While Rafflesia buds take up to 21 months to develop, its flowers bloom and wither within about a week. In this study, transcriptome analysis was carried out to shed light on the molecular mechanism of senescence in Rafflesia. A total of 53.3 million high quality reads were obtained from two Rafflesia cantleyi flower developmental stages and assembled to generate 64,152 unigenes. Analysis of this dataset showed that 5,166 unigenes were differentially expressed, in which 1,073 unigenes were identified as genes involved in flower senescence. Results revealed that as the flowers progress to senescence, more genes related to flower senescence were significantly over-represented compared to those related to plant growth and development. Senescence of the R. cantleyi flower activates senescence-associated genes in the transcription activity (members of the transcription factor families MYB, bHLH, NAC, and WRKY), nutrient remobilization (autophagy-related protein and transporter genes), and redox regulation (CATALASE). Most of the senescence-related genes were found to be differentially regulated, perhaps for the fine-tuning of various responses in the senescing R. cantleyi flower. Additionally, pathway analysis showed the activation of genes such as ETHYLENE RECEPTOR, ETHYLENE-INSENSITIVE 2, ETHYLENE-INSENSITIVE 3, and ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR, indicating the possible involvement of the ethylene hormone response pathway in the regulation of R. cantleyi senescence. Our results provide a model of the molecular mechanism underlying R. cantleyi flower senescence, and contribute essential information towards further understanding the biology of the Rafflesiaceae family.

Functional Characterization of the EMBRYONIC FLOWER 2 Gene Involved in Flowering in Ginkgo biloba

Ginkgo biloba has edible, medicinal, and ornamental value. However, the long juvenile phase prevents the development of the G. biloba industry, and there are few reports on the identification and functional analysis of genes regulating the flowering time of G. biloba. EMBRYONIC FLOWER 2 (EMF), an important protein in flower development, functions to promote vegetative growth and repress flowering. In this study, a novel EMF gene (GbEMF2) was cloned and characterized from G. biloba. GbEMF2 contains a 2,193 bp open reading frame (ORF) encoding 730 amino acids. GbEMF2 harbors conserved VEFS-Box domain by the plant EMF protein. The phylogenic analysis showed that GbEMF2 originated from a polycomb-group (Pc-G) protein ancestor and was a member of the EMF2 protein. The quantitative real-time PCR (qRT-PCR) analysis revealed that GbEMF2 was expressed in all detected organs, and it showed a significantly higher level in ovulating strobilus and microstrobilus than in other organs. Compared with emf2 mutant plants, overexpression of GbEMF2 driven by the CaMV 35S promoter in emf2 mutant Arabidopsis plants delayed flowering but earlier than wild-type (WT) plants. This result indicated that GbEMF2 repressed flowering in G. biloba. Moreover, the RNA-seq analysis of GbEMF2 transgenic Arabidopsis plants (GbEMF2-OE/emf2), WT plants, and emf2 mutants screened out 227 differentially expressed genes (DEGs). Among these DEGs, FLC, MAF5, and MAF5-1 genes were related to flower organ development and regulated by GbEMF2. In addition, some genes participating in sugar metabolism, such as Alpha-amylase 1 (AMY1), BAM1, and Sucrose synthase 3 (SUS3) genes, were also controlled by GbEMF2. Overall, our results suggested that GbEMF2 negatively regulates flowering development in G. biloba. This finding provided a foundation and target gene for shortening the Ginkgo juvenile period by genetic engineering technology.