Transcriptomic and hormone analyses reveal mechanisms underlying petal elongation in Chrysanthemum morifolium 'Jinba'

Heterologous Expression of Chrysanthemum TCP Transcription Factor CmTCP13 Enhances Salinity Tolerance in Arabidopsis

Plant-specific TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) proteins play critical roles in plant development and stress responses; however, their functions in chrysanthemum (Chrysanthemum morifolium) have not been well-studied. In this study, we isolated and characterized the chrysanthemum TCP transcription factor family gene CmTCP13, a homolog of AtTCP13. This gene encoded a protein harboring a conserved basic helix-loop-helix motif, and its expression was induced by salinity stress in chrysanthemum plants. Subcellular localization experiments showed that CmTCP13 localized in the nucleus. Sequence analysis revealed the presence of multiple stress- and hormone-responsive cis-elements in the promoter region of CmTCP13. The heterologous expression of CmTCP13 in Arabidopsis plants enhanced their tolerance to salinity stress. Under salinity stress, CmTCP13 transgenic plants exhibited enhanced germination, root length, seedling growth, and chlorophyll content and reduced relative electrical conductivity compared with those exhibited by wild-type (WT) plants. Moreover, the expression levels of stress-related genes, including AtSOS3, AtP5CS2, AtRD22, AtRD29A, and AtDREB2A, were upregulated in CmTCP13 transgenic plants than in WT plants under salt stress. Taken together, our results demonstrate that CmTCP13 is a critical regulator of salt stress tolerance in plants.

FUL homologous gene CmFL1 is involved in regulating flowering time and floret numbers in Chrysanthemum morifolium

Flowering time and floret numbers are important ornamental characteristics of chrysanthemums that control their adaptability and inflorescence morphology, respectively. The FRUITFULL (FUL) gene plays a key role in inducing flowering and inflorescence meristem development. In this study, we isolated a homolog of the MADS-box gene FUL, CmFUL-Like 1 (CmFL1), from chrysanthemum inflorescence buds. Quantitative RT-PCR and in situ analyses showed that CmFL1 was strongly expressed in young inflorescence buds. Overexpression of CmFL1 caused early flowering while co-suppression expression of CmFL1 increased the number of florets. Furthermore, the floral promoting factors CmSOC1, CmFDL1, and CmLFY were up-regulated in the shoot tips of transgenic plants. In addition, RNA-seq analysis of the transgenic plants suggested that certain differentially expressed genes (DEGs)-such as MADS-box, homeobox family, and ethylene pathway genes-may be involved in the inflorescence meristem development. GO pathway enrichment analysis showed that the differentially transcribed genes enriched the representation of the carbohydrate metabolic pathway, which is critical for flower development. Overall, our findings revealed the conserved function of CmFL1 in controlling flowering time along with a novel function in regulating the number of florets.

The ARF2-MYB6 module mediates auxin-regulated petal expansion in rose

In cut rose (Rosa hybrida), the flower-opening process is closely associated with vase life. Auxin induces the expression of transcription factor genes that function in petal growth via cell expansion. However, the molecular mechanisms underlying the auxin effect during flower opening are not well understood. Here, we identified the auxin-inducible transcription factor gene RhMYB6, whose expression level is high during the early stages of flower opening. Silencing of RhMYB6 delayed flower opening by controlling petal cell expansion through down-regulation of cell expansion-related genes. Furthermore, we demonstrated that the auxin response factor RhARF2 directly interacts with the promoter of RhMYB6 and represses its transcription. Silencing of RhARF2 resulted in larger petal size and delayed petal movement. We also showed that the expression of genes related to ethylene and petal movement showed substantial differences in RhARF2-silenced petals. Our results indicate that auxin-regulated RhARF2 is a critical player that controls flower opening by governing RhMYB6 expression and mediating the crosstalk between auxin and ethylene signaling.

Combined Transcriptome and Proteome Analysis Provides Insights into Petaloidy in Pomegranate

Petaloidy leads to a plump floral pattern and increases the landscape value of ornamental pomegranates; however, research on the mechanism of petaloidy in ornamental pomegranates is limited. In this study, we aimed to screen candidate genes related to petaloidy. We performed transcriptomic and proteomic sequencing of the stamens and petals of single-petal and double-petal flowers of ornamental pomegranates. Briefly, 24,567 genes and 5865 proteins were identified, of which 5721 genes were quantified at both transcriptional and translational levels. In the petal and stamen comparison groups, the association between differentially abundant proteins (DAPs) and differentially expressed genes (DEGs) was higher than that between all genes and all proteins, indicating that petaloidy impacts the correlation between genes and proteins. The enrichment results of transcriptome, proteome, and correlation analyses showed that cell wall metabolism, jasmonic acid signal transduction, redox balance, and transmembrane transport affected petaloidy. Nine hormone-related DEGs/DAPs were selected, among which ARF, ILR1, LAX2, and JAR1 may promote petal doubling. Sixteen transcription factor DEGs/DAPs were selected, among which EREBP, LOB, MEF2, MYB, C3H, and trihelix may promote petal doubling. Our results provide transcriptomic and proteomic data on the formation mechanism of petaloidy and a theoretical basis for breeding new ornamental pomegranate varieties.

Analyses of a chromosome-scale genome assembly reveal the origin and evolution of cultivated chrysanthemum

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a globally important ornamental plant with great economic, cultural, and symbolic value. However, research on chrysanthemum is challenging due to its complex genetic background. Here, we report a near-complete assembly and annotation for C. morifolium comprising 27 pseudochromosomes (8.15 Gb; scaffold N50 of 303.69 Mb). Comparative and evolutionary analyses reveal a whole-genome triplication (WGT) event shared by Chrysanthemum species approximately 6 million years ago (Mya) and the possible lineage-specific polyploidization of C. morifolium approximately 3 Mya. Multilevel evidence suggests that C. morifolium is likely a segmental allopolyploid. Furthermore, a combination of genomics and transcriptomics approaches demonstrate the C. morifolium genome can be used to identify genes underlying key ornamental traits. Phylogenetic analysis of CmCCD4a traces the flower colour breeding history of cultivated chrysanthemum. Genomic resources generated from this study could help to accelerate chrysanthemum genetic improvement.

Floral Development Stage-Specific Transcriptomic Analysis Reveals the Formation Mechanism of Different Shapes of Ray Florets in Chrysanthemum

The formation mechanism of different ray floret shapes of chrysanthemum (Chrysanthemum × morifolium) remains elusive due to its complex genetic background. C. vestitum, with the basic ray floret shapes of the flat, spoon, and tubular types, is considered a model material for studying ray floret morphogenesis. In this study, the flat and tubular type lines of C. vestitum at specific stages were used to investigate the key genes that regulate morphological differences in ray florets. We found that the expression levels of genes related to auxin synthesis, transport, and response were generally higher in the tubular type than in the flat type. CvARF3 was highly expressed in the flat type, while CvARF5 and CvARF6 were highly expressed in the tubular type. Additionally, the transcription levels of Class B and E genes closely related to petal development, including CvPI, CvAP3, Cvdefh21, CvSEP3, and CvCDM77, were expressed at higher levels in the tubular type than the flat type. Based on the results, it is proposed that auxin plays a key role in the development of ray florets, and auxin-related genes, especially CvARFs, may be key genes to control the morphological difference of ray florets. Simultaneously, MADS-box genes are involved in the co-regulation of ray floret morphogenesis. The results provide novel insights into the molecular mechanism of different petal type formation and lay a theoretical foundation for the directional breeding of petal type in chrysanthemums.

Transcriptomic analysis reveals the formation mechanism of anemone-type flower in chrysanthemum

BACKGROUND

The ray and disc florets on the chrysanthemum capitulum are morphologically diverse and have remarkably abundant variant types, resulting in a rich variety of flower types. An anemone shape with pigmented and elongated disk florets is an important trait in flower shape breeding of chrysanthemums. The regulatory mechanism of their anemone-type disc floret formation was not clear, thus limiting the directional breeding of chrysanthemum flower types. In this study, we used morphological observation, transcriptomic analysis, and gene expression to investigate the morphogenetic processes and regulatory mechanisms of anemone-type chrysanthemum.

RESULT

Scanning electron microscopy (SEM) observation showed that morphological differences between non-anemone-type disc florets and anemone-type disc florets occurred mainly during the petal elongation period. The anemone-type disc florets elongated rapidly in the later stages of development. Longitudinal paraffin section analysis revealed that the anemone-type disc florets were formed by a great number of cells in the middle layer of the petals with vigorous division. We investigated the differentially expressed genes (DEGs) using ray and disc florets of two chrysanthemum cultivars, 082 and 068, for RNA-Seq and their expression patterns of non-anemone-type and anemone-type disc florets. The result suggested that the CYCLOIDEA2 (CYC2s), MADS-box genes, and phytohormone signal-related genes appeared significantly different in both types of disc florets and might have important effects on the formation of anemone-type disc florets. In addition, it is noteworthy that the auxin and jasmonate signaling pathways might play a vital role in developing anemone-type disc florets.

CONCLUSIONS

Based on our findings, we propose a regulatory network for forming non-anemone-type and anemone-type disc florets. The results of this study lead the way to further clarify the mechanism of the anemone-type chrysanthemum formation and lay the foundation for the directive breeding of chrysanthemum petal types.

The TIFY family protein CmJAZ1-like negatively regulates petal size via interaction with the bHLH transcription factor CmBPE2 in Chrysanthemum morifolium

Petals are the second floral whorl of angiosperms, exhibiting astonishing diversity in their size between and within species. This variation is essential for protecting their inner reproductive organs and attracting pollinators for fertilization. However, currently, the genetic and developmental control of petal size remains unexplored. Chrysanthemum (Chrysanthemum morifolium) belongs to the Asteraceae family, the largest group of angiosperms, and the extraordinary diversity of petal size in chrysanthemums makes it an ideal model for exploring the regulation mechanism of petal size. Here, we reveal that overexpression of a JAZ repressor CmJAZ1-like exhibits decreased petal size compared to that of the wild-type as a result of repressed cell expansion. Through further in-depth exploration, we confirm an interaction pair between CmJAZ1-like and the bHLH transcription factor CmBPE2. The inhibition of CmBPE2 expression negatively regulates petal size by downregulating the expression of genes involved in cell expansion. Furthermore, CmJAZ1-like significantly reduced the activation ability of CmBPE2 on its target gene CmEXPA7 by directly interacting with it, thus participating in the regulation of petal size development in chrysanthemum. Our results will provide insights into the molecular mechanisms of petal size regulation in flowering plants.

Combination of long-read and short-read sequencing provides comprehensive transcriptome and new insight for Chrysanthemum morifolium ray-floret colorization

Chrysanthemum morifolium is one of the most popular ornamental plants globally. Owing to its large and complex genome (around 10 Gb, segmental hexaploid), it has been difficult to obtain comprehensive transcriptome, which will promote to perform new breeding technique, such as genome editing, in C. morifolium. In this study, we used single-molecule real-time (SMRT) sequencing and RNA-seq technologies, combined them with an error-correcting process, and obtained high-coverage ray-floret transcriptome. The SMRT-seq data increased the ratio of long mRNAs containing complete open-reading frames, and the combined dataset provided a more complete transcriptomic data than those produced from either SMRT-seq or RNA-seq-derived transcripts. We finally obtained 'Sei Arabella' transcripts containing 928,645 non-redundant mRNA, which showed 96.6% Benchmarking Universal Single-Copy Orthologs (BUSCO) score. We also validated the reliability of the dataset by analyzing a mapping rate, annotation and transcript expression. Using the dataset, we searched anthocyanin biosynthesis gene orthologs and performed a qRT-PCR experiment to assess the usability of the dataset. The assessment of the dataset and the following analysis indicated that our dataset is reliable and useful for molecular biology. The combination of sequencing methods provided genetic information and a way to analyze the complicated C. morifolium transcriptome.

Cytological, physiological and transcriptomic analysis of variegated Leaves in Primulina pungentisepala offspring

BACKGROUND

Primulina pungentisepala is suitable for use as a potted plant because of its beautiful leaf variegation, which is significantly different in its selfed offspring. However, the mechanism of P. pungentisepala leaf variegation is unclear. In this study, two types of offspring showing the greatest differences were compared in terms of leaf structure, chlorophyll contents, chlorophyll fluorescence parameters and transcriptomes to provide a reference for studying the molecular mechanism of structural leaf variegation.

RESULTS

Air spaces were found between water storage tissue, and the palisade tissue cells were spherical in the white type. The content of chlorophyll a and total chlorophyll (chlorophyll a + b) was significantly lower in the white type, but there were no significant differences in the content of chlorophyll b, chlorophyll a/b or chlorophyll fluorescence parameters between the white and green types. We performed transcriptomic sequencing to identify differentially expressed genes (DEGs) involved in cell division and differentiation, chlorophyll metabolism and photosynthesis. Among these genes, the expression of the cell division- and differentiation-related leucine-rich repeat receptor-like kinases (LRR-RLKs), xyloglucan endotransglycosylase/hydrolase (XET/H), pectinesterase (PE), expansin (EXP), cellulose synthase-like (CSL), VARIEGATED 3 (VAR3), and ZAT10 genes were downregulated in the white type, which might have promoted the development air spaces and variant palisade cells. Chlorophyll biosynthesis-related hydroxymethylbilane synthase (HEMC) and the H subunit of magnesium chelatase (CHLH) were downregulated, while chlorophyll degradation-related chlorophyllase-2 (CHL2) was upregulated in the white type, which might have led to lower chlorophyll accumulation.

CONCLUSION

Leaf variegation in P. pungentisepala was caused by a combination of mechanisms involving structural variegation and low chlorophyll levels. Our research provides significant insights into the molecular mechanisms of structural leaf variegation.

Integrated Transcriptome and Targeted Metabolite Analysis Reveal miRNA-mRNA Networks in Low-Light-Induced Lotus Flower Bud Abortion

Most Nelumbo nucifera (lotus) flower buds were aborted during the growing season, notably in low-light environments. How lotus produces so many aborted flower buds is largely unknown. An integrated transcriptome and targeted metabolite analysis was performed to reveal the genetic regulatory networks underlying lotus flower bud abortion. A total of 233 miRNAs and 25,351 genes were identified in lotus flower buds, including 68 novel miRNAs and 1108 novel genes. Further enrichment analysis indicated that sugar signaling plays a potential central role in regulating lotus flower bud abortion. Targeted metabolite analysis showed that trehalose levels declined the most in the aborting flower buds. A potential regulatory network centered on miR156 governs lotus flower bud abortion, involving multiple miRNA-mRNA pairs related to cell integrity, cell proliferation and expansion, and DNA repair. Genetic analysis showed that miRNA156-5p-overexpressing lotus showed aggravated flower bud abortion phenotypes. Trehalose-6-P synthase 1 (TPS1), which is required for trehalose synthase, had a negative regulatory effect on miR156 expression. TPS1-overexpression lotus showed significantly decreased flower bud abortion rates both in normal-light and low-light environments. Our study establishes a possible genetic basis for how lotus produces so many aborted flower buds, facilitating genetic improvement of lotus’ shade tolerance.

CRISPR/Cas9 System: A Potential Tool for Genetic Improvement in Floricultural Crops

Demand of flowers is increasing with time worldwide. Floriculture has become one of the most important commercial trades in agriculture. Although traditional breeding methods like hybridization and mutation breeding have contributed significantly to the development of important flower varieties, flower production and quality of flowers can be significantly improved by employing modern breeding approaches. Novel traits of significance have interest to consumers and producers, such as fragrance, new floral color, change in floral architecture and morphology, vase life, aroma, and resistance to biotic and abiotic stresses, have been introduced by genetic manipulation. The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) system has recently emerged as a powerful genome-editing tool for accurately changing DNA sequences at specific locations. It provides excellent means of genetically improving floricultural crops. CRISPR/Cas system has been utilized in gene editing in horticultural cops. There are few reports on the utilization of the CRISPR/Cas9 system in flowers. The current review summarizes the research work done by employing the CRISPR/Cas9 system in floricultural crops including improvement in flowering traits such as color modification, prolonging the shelf life of flowers, flower initiation, and development, changes in color of ornamental foliage by genome editing. CRISPR/Cas9 gene editing could be useful in developing novel cultivars with higher fragrance and enhanced essential oil and many other useful traits. The present review also highlights the basic mechanism and key components involved in the CRISPR/Cas9 system.

An Aux/IAA Family Member, RhIAA14, Involved in Ethylene-Inhibited Petal Expansion in Rose (Rosa hybrida)

Flower size, a primary agronomic trait in breeding of ornamental plants, is largely determined by petal expansion. Generally, ethylene acts as an inhibitor of petal expansion, but its effect is restricted by unknown developmental cues. In this study, we found that the critical node of ethylene-inhibited petal expansion is between stages 1 and 2 of rose flower opening. To uncover the underlying regulatory mechanism, we carried out a comparative RNA-seq analysis. Differentially expressed genes (DEGs) involved in auxin-signaling pathways were enriched. Therefore, we identified an auxin/indole-3-acetic acid (Aux/IAA) family gene, RhIAA14, whose expression was development-specifically repressed by ethylene. The silencing of RhIAA14 reduced cell expansion, resulting in diminished petal expansion and flower size. In addition, the expressions of cell-expansion-related genes, including RhXTH6, RhCesA2, RhPIP2;1, and RhEXPA8, were significantly downregulated following RhIAA14 silencing. Our results reveal an Aux/IAA that serves as a key player in orchestrating petal expansion and ultimately contributes to flower size, which provides new insights into ethylene-modulated flower opening and the function of the Aux/IAA transcription regulator.

Characterization of WRKY Gene Family in Whole-Genome and Exploration of Flowering Improvement Genes in Chrysanthemum lavandulifolium

Chrysanthemum is a well-known ornamental plant with numerous uses. WRKY is a large family of transcription factors known for a variety of functions ranging from stress resistance to plant growth and development. Due to the limited research on the WRKY family in chrysanthemums, we examined them for the first time in Chrysanthemum lavandulifolium. A total of 138 ClWRKY genes were identified, which were classified into three groups. Group III in C. lavandulifolium contains 53 members, which is larger than group III of Arabidopsis. The number of introns varied from one to nine in the ClWRKY gene family. The "WRKYGQK" motif is conserved in 118 members, while other members showed slight variations. AuR and GRE responsive cis-acting elements were located in the promoter region of WRKY members, which are important for plant development and flowering induction. In addition, the W box was present in most genes; the recognition site for the WRKY gene may play a role in autoregulation and cross-regulation. The expression of the most variable 19 genes in terms of different parameters was observed at different stages. Among them, 10 genes were selected due to the presence of CpG islands, while nine genes were selected based on their close association with important Arabidopsis genes related to floral traits. ClWRKY36 and ClWRKY45 exhibit differential expression at flowering stages in the capitulum, while methylation is detected in three genes, including ClWRKY31, ClWRKY100, and ClWRKY129. Our results provide a basis for further exploration of WRKY members to find their functions in plant growth and development, especially in flowering traits.

A regulatory GhBPE-GhPRGL module maintains ray petal length in Gerbera hybrida

The molecular mechanism regulating petal length in flowers is not well understood. Here we used transient transformation assays to confirm that GhPRGL (proline-rich and GASA-like)-a GASA (gibberellic acid [GA] stimulated in Arabidopsis) family gene-promotes the elongation of ray petals in gerbera (Gerbera hybrida). Yeast one-hybrid screening assay identified a bHLH transcription factor of the jasmonic acid (JA) signaling pathway, here named GhBPE (BIGPETAL), which binds to the GhPRGL promoter and represses its expression, resulting in a phenotype of shortened ray petal length when GhBPE is overexpressed. Further, the joint response to JA and GA of GhBPE and GhPRGL, together with their complementary expression profiles in the early stage of petal growth, suggests a novel GhBPE-GhPRGL module that controls the size of ray petals. GhPRGL promotes ray petal elongation in its early stage especially, while GhBPE inhibits ray petal elongation particularly in the late stage by inhibiting the expression of GhPRGL. JA and GA operate in concert to regulate the expression of GhBPE and GhPRGL genes, providing a regulatory mechanism by which ray petals could grow to a fixed length in gerbera species.

Comparison of chrysanthemum flowers grown under hydroponic and soil-based systems: yield and transcriptome analysis

BACKGROUND

Flowers of Chrysanthemum × morifolium Ramat. are used as tea in traditional Chinese cuisine. However, with increasing population and urbanization, water and land availability have become limiting for chrysanthemum tea production. Hydroponic culture enables effective, rapid nutrient exchange, while requiring no soil and less water than soil cultivation. Hydroponic culture can reduce pesticide residues in food and improve the quantity or size of fruits, flowers, and leaves, and the levels of active compounds important for nutrition and health. To date, studies to improve the yield and active compounds of chrysanthemum have focused on soil culture. Moreover, the molecular effects of hydroponic and soil culture on chrysanthemum tea development remain understudied.

RESULTS

Here, we studied the effects of soil and hydroponic culture on yield and total flavonoid and chlorogenic acid contents in chrysanthemum flowers (C. morifolium 'wuyuanhuang'). Yield and the total flavonoids and chlorogenic acid contents of chrysanthemum flowers were higher in the hydroponic culture system than in the soil system. Transcriptome profiling using RNA-seq revealed 3858 differentially expressed genes (DEGs) between chrysanthemum flowers grown in soil and hydroponic conditions. Gene Ontology (GO) enrichment annotation revealed that these differentially transcribed genes are mainly involved in "cytoplasmic part", "biosynthetic process", "organic substance biosynthetic process", "cell wall organization or biogenesis" and other processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed enrichment in "metabolic pathways", "biosynthesis of secondary metabolites", "ribosome", "carbon metabolism", "plant hormone signal transduction" and other metabolic processes. In functional annotations, pathways related to yield and formation of the main active compounds included phytohormone signaling, secondary metabolism, and cell wall metabolism. Enrichment analysis of transcription factors also showed that under the hydroponic system, bHLH, MYB, NAC, and ERF protein families were involved in metabolic pathways, biosynthesis of secondary metabolites, and plant hormone signal transduction.

CONCLUSIONS

Hydroponic culture is a simple and effective way to cultivate chrysanthemum for tea production. A transcriptome analysis of chrysanthemum flowers grown in soil and hydroponic conditions. The large number of DEGs identified confirmed the difference of the regulatory machinery under two culture system.

Molecular understanding of postharvest flower opening and senescence

Flowers are key organs in many ornamental plants, and various phases of flower development impact their economic value. The final stage of petal development is associated with flower senescence, which is an irreversible process involving programmed cell death, and premature senescence of cut flowers often results in major losses in quality during postharvest handling. Flower opening and senescence are two sequential processes. As flowers open, the stamens are exposed to attract pollinators. Once pollination occurs, flower senescence is initiated. Both the opening and senescence processes are regulated by a range of endogenous phytohormones and environmental factors. Ethylene acts as a central regulator for the ethylene-sensitive flowers. Other phytohormones, including auxin, gibberellin, cytokinin, jasmonic acid and abscisic acid, are also involved in the control of petal expansion and senescence. Water status also directly influences postharvest flower opening, while pollination is a key event in initiating the onset flower senescence. Here, we review the current understanding of flower opening and senescence, and propose future research directions, such as the study of interactions between hormonal and environmental signals, the application of new technology, and interdisciplinary research.

Overexpression of the CmJAZ1-like gene delays flowering in Chrysanthemum morifolium

Chrysanthemum (Chrysanthemum morifolium) is one of the four major cut-flower plants worldwide and possesses both high ornamental value and cultural connotation. As most chrysanthemum varieties flower in autumn, it is costly to achieve annual production. JAZ genes in the TIFY family are core components of the jasmonic acid (JA) signaling pathway; in addition to playing a pivotal role in plant responses to defense, they are also widely implicated in regulating plant development processes. Here, we characterized the TIFY family gene CmJAZ1-like from the chrysanthemum cultivar 'Jinba'. CmJAZ1-like localizes in the nucleus and has no transcriptional activity in yeast. Tissue expression pattern analysis indicated that CmJAZ1-like was most active in the root and shoot apex. Overexpressing CmJAZ1-like with Jas domain deletion in chrysanthemum resulted in late flowering. RNA-Seq analysis of the overexpression lines revealed some differentially expressed genes (DEGs) involved in flowering, such as the homologs of the flowering integrators FT and SOC1, an FUL homolog involved in flower meristem identity, AP2 domain-containing transcription factors, MADS box genes, and autonomous pathway-related genes. Based on KEGG pathway enrichment analysis, the differentially transcribed genes were enriched in carbohydrate metabolic and fatty acid-related pathways, which are notable for their role in flowering in plants. This study preliminarily verified the function of CmJAZ1-like in chrysanthemum flowering, and the results can be used in molecular breeding programs aimed at flowering time regulation of chrysanthemum.

Comparative transcriptome analysis of flower bud transition and functional characterization of EjAGL17 involved in regulating floral initiation in loquat

Floral initiation plays a critical role for reproductive success in plants, especially fruit trees. However, little information is known on the mechanism of the initiation in loquat (Eriobotrya japonica Lindl.). Here, we used transcriptomic, expression and functional analysis to investigate the candidate genes in floral initiation in loquat. Comparative transcriptome analysis showed differentially expressed genes (DEGs) were mainly enriched in the metabolic pathways of plant hormone signal transduction. The DEGs were mainly involved in the gibberellin, auxin, cytokinin, abscisic acid, salicylic acid and ethylene signaling pathways. Meanwhile, some transcription factors, including MADS-box (MCM1, AGAMOUS, DEFICIENS and SRF), MYB (Myeloblastosis), TCP (TEOSINTE BRANCHED 1, CYCLOIDEA and PCF1), WOX (WUSCHEL-related homeobox) and WRKY (WRKY DNA-binding protein), were significantly differentially expressed. Among these key DEGs, we confirmed that an AGL17 ortholog EjAGL17 was significantly upregulated at the flower bud transition stage. Phylogenetic tree analysis revealed that EjAGL17 was grouped into an AGL17 clade of MADS-box transcription factors. Protein sequence alignment showed that EjAGL17 included a distinctive C-terminal domain. Subcellular localization of EjAGL17 was found only in the nucleus. Expression levels of EjAGL17 reached the highest at the development stage of flower bud transition. Moreover, ectopic expression of EjAGL17 in Arabidopsis significantly exhibited early flowering. Our study provides abundant resources of candidate genes for studying the mechanisms underlying the floral initiation in loquat and other Rosaceae species.

Differential Gene Expression with an Emphasis on Floral Organ Size Differences in Natural and Synthetic Polyploids of Nicotiana tabacum (Solanaceae)

Floral organ size, especially the size of the corolla, plays an important role in plant reproduction by facilitating pollination efficiency. Previous studies have outlined a hypothesized organ size pathway. However, the expression and function of many of the genes in the pathway have only been investigated in model diploid species; therefore, it is unknown how these genes interact in polyploid species. Although correlations between ploidy and cell size have been shown in many systems, it is unclear whether there is a difference in cell size between naturally occurring and synthetic polyploids. To address these questions comparing floral organ size and cell size across ploidy, we use natural and synthetic polyploids of Nicotiana tabacum (Solanaceae) as well as their known diploid progenitors. We employ a comparative transcriptomics approach to perform analyses of differential gene expression, focusing on candidate genes that may be involved in floral organ size, both across developmental stages and across accessions. We see differential expression of several known floral organ candidate genes including ARF2, BIG BROTHER, and GASA/GAST1. Results from linear models show that ploidy, cell width, and cell number positively influence corolla tube circumference; however, the effect of cell width varies by ploidy, and diploids have a significantly steeper slope than both natural and synthetic polyploids. These results demonstrate that polyploids have wider cells and that polyploidy significantly increases corolla tube circumference.

The core regulatory networks and hub genes regulating flower development in Chrysanthemum morifolium

The study has facilitated important insights into the regulatory networks involved in flower development in chrysanthemum (Asteraceae), and is informative with respect to the mechanism of flower shape determination. Chrysanthemum morifolium, valued as an ornamental species given the diversity of its inflorescence form, is viewed as a model for understanding flower development in the Asteraceae. Yet, the underlying regulatory networks remain largely unexplored. Here, a transcriptomic survey of the Chrysanthemum morifolium variety 'Jinba' was undertaken to uncover the global gene expression profiles and identify the modules of co-transcribed genes associated with flower development. The weighted gene coexpression network analysis revealed important networks and hub genes including ray floret petals-specific coexpression network, disc floret petals-specific network, B and E class genes involved network and CYC2 genes network. Three ray floret petal-specific hub genes were also strongly transcribed in the ray florets of a selection of six diverse varieties and especially so in those which form ligulate ray floret petals. CmCYC2c was strongly transcribed in the distal and lateral regions of the ray floret petals, and also, along with CmCYC2d, in the tubular ray florets. Furthermore, CmOFP, belonging to the family of ovate proteins, was identified in the CYC2 genes network. CmOFP can interact with CmCYC2d that physically interact with CmCYC2c. This work provides important insights into the regulatory networks involved in flower development in chrysanthemum, and is informative with respect to the mechanistic basis of the regulation of flower shape.

An Integrative Analysis of Transcriptome, Proteome and Hormones Reveals Key Differentially Expressed Genes and Metabolic Pathways Involved in Flower Development in Loquat

Flower development is a vital developmental process in the life cycle of woody perennials, especially fruit trees. Herein, we used transcriptomic, proteomic, and hormone analyses to investigate the key candidate genes/proteins in loquat (Eriobotrya japonica) at the stages of flower bud differentiation (FBD), floral bud elongation (FBE), and floral anthesis (FA). Comparative transcriptome analysis showed that differentially expressed genes (DEGs) were mainly enriched in metabolic pathways of hormone signal transduction and starch and sucrose metabolism. Importantly, the DEGs of hormone signal transduction were significantly involved in the signaling pathways of auxin, gibberellins (GAs), cytokinin, ethylene, abscisic acid (ABA), jasmonic acid, and salicylic acid. Meanwhile, key floral integrator genes FLOWERING LOCUS T (FT) and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) and floral meristem identity genes SQUAMOSA PROMOTER BINDING LIKE (SPL), LEAFY (LFY), APETALA1 (AP1), and AP2 were significantly upregulated at the FBD stage. However, key floral organ identity genes AGAMOUS (AG), AP3, and PISTILLATA (PI) were significantly upregulated at the stages of FBE and FA. Furthermore, transcription factors (TFs) such as bHLH (basic helix-loop-helix), NAC (no apical meristem (NAM), Arabidopsis transcription activation factor (ATAF1/2) and cup-shaped cotyledon (CUC2)), MYB_related (myeloblastosis_related), ERF (ethylene response factor), and C2H2 (cysteine-2/histidine-2) were also significantly differentially expressed. Accordingly, comparative proteomic analysis of differentially accumulated proteins (DAPs) and combined enrichment of DEGs and DAPs showed that starch and sucrose metabolism was also significantly enriched. Concentrations of GA3 and zeatin were high before the FA stage, but ABA concentration remained high at the FA stage. Our results provide abundant sequence resources for clarifying the underlying mechanisms of the flower development in loquat.

Investigation of genes associated with petal variations between diploid and autotetraploid in Chinese cabbage (Brassica rapa L. ssp. pekinensis) by RNA-seq and sRNA-seq

Polyploidy promotes morphological, physiological, and reproductive diversity in plants. The imminent effect of chromosome doubling in plants is the enlargement of organs such as flowers and fruits, which increases the commercial value of crops. Flowering plays a vital role in the growth and development of angiosperms. Here, we prepared an isolated microspore culture of 'FT', a doubled haploid (DH) line of Chinese cabbage (Brassica rapa L. ssp. pekinensis), and obtained diploid and autotetraploid plants with the same genetic background. Compared with diploids, the autotetraploids were characterized by large floral organs, dark petals, delayed flowering, and reduced fertility. The indole-3-acetic acid (IAA) and jasmonic acid (JA) levels in autotetraploid petals were significantly higher and the abscisic acid (ABA) level was significantly lower than those in the diploid petals. The lutein level in autotetraploid petals was nearly two times higher than that in the diploid petals. A comparative transcriptome analysis revealed 14,412 differentially expressed genes (DEGs) between the diploids and autotetraploids, and they were enriched in 117 Gene Ontology terms and 110 Kyoto Encyclopedia of Genes and Genomes pathways. We detected 231 DEGs related to phytohormone signal transduction and 29 DEGs involved in carotenoid biosynthesis. An miRNA-target mRNA analysis showed that 32 DEGs regulated by 16 DEMs were associated with flowering timing (BraA03000336, BraA09004319, and BraA09000515), petal development (BraA05002408, BraA01004006, BraA09004069, and BraA04000966), flower opening (BraA07000350), and pollen development (BraA01000720, BraA09005727, and BraA01000253). This study provides information to help elucidate the molecular mechanisms underlying phenotypic variations induced by autopolyploidy in Chinese cabbage.

Transcriptomic analysis of flower opening response to relatively low temperatures in Osmanthus fragrans

BACKGROUND

Sweet osmanthus (Osmanthus fragrans Lour.) is one of the top ten traditional ornamental flowers in China. The flowering time of once-flowering cultivars in O. fragrans is greatly affected by the relatively low temperature, but there are few reports on its molecular mechanism to date. A hypothesis had been raised that genes related with flower opening might be up-regulated in response to relatively low temperature in O. fragrans. Thus, our work was aimed to explore the underlying molecular mechanism of flower opening regulated by relatively low temperature in O. fragrans.

RESULTS

The cell size of adaxial and abaxial petal epidermal cells and ultrastructural morphology of petal cells at different developmental stages were observed. The cell size of adaxial and abaxial petal epidermal cells increased gradually with the process of flower opening. Then the transcriptomic sequencing was employed to analyze the differentially expressed genes (DEGs) under different number of days' treatments with relatively low temperatures (19 °C) or 23 °C. Analysis of DEGs in Gene Ontology analysis showed that "metabolic process", "cellular process", "binding", "catalytic activity", "cell", "cell part", "membrane", "membrane part", "single-organism process", and "organelle" were highly enriched. In KEGG analysis, "metabolic pathways", "biosynthesis of secondary metabolites", "plant-pathogen interaction", "starch and sucrose metabolism", and "plant hormone signal transduction" were the top five pathways containing the greatest number of DEGs. The DEGs involved in cell wall metabolism, phytohormone signal transduction pathways, and eight kinds of transcription factors were analyzed in depth.

CONCLUSIONS

Several unigenes involved in cell wall metabolism, phytohormone signal transduction pathway, and transcription factors with highly variable expression levels between different temperature treatments may be involved in petal cell expansion during flower opening process in response to the relatively low temperature. These results could improve our understanding of the molecular mechanism of relatively-low-temperature-regulated flower opening of O. fragrans, provide practical information for the prediction and regulation of flowering time in O. fragrans, and ultimately pave the way for genetic modification in O. fragrans.

Homeotic transformation from stamen to petal in Eriobotrya japonica is associated with hormone signal transduction and reduction of the transcriptional activity of EjAG

Double-flower loquat (Eriobotrya japonica) is a new germplasm with homeotic transformation of stamen into petal in whorl 3. However, little information is available on the molecular mechanism of this transformation. Herein, we analyzed the transcriptome, candidate genes and endogenous hormones to investigate the mechanisms underlying this homeotic transformation. Some transcription factors, such as MADS-box, TCP and MYB, were significantly differentially expressed. Importantly, we confirmed that one of these (DN39625_c0_g1), which encoded a C-class floral homeotic protein referred to as AGAMOUS ortholog (EjAG), was significantly downregulated. Subcellular localization of EjAG was found to be in the nucleus. Ectopic expression of EjAG rescued the development of stamens and carpels from the double-flower phenotype in an Arabidopsis ag mutant, suggesting that EjAG expression is associated with double-flower formation. Meanwhile, enrichment analyses showed that the differentially expressed genes (DEGs) were mainly involved in the metabolic pathways of hormone signal transduction. The DEGs of auxin, gibberellin A (GA) and cytokinin signaling pathways were mainly upregulated. However, the DEGs of abscisic acid (ABA) and the ethylene signaling pathway were mainly downregulated. Accordingly, the concentrations of indoleacetic acid, kinetin and GA₃ were high at the petaloid stamen stage, but the ABA concentration remained low. The identified genes and pathways provide abundant sequence resources for studying the mechanisms underlying the homeotic transformation in loquat and other Rosaceae species.

Comprehensive characterization of a floral mutant reveals the mechanism of hooked petal morphogenesis in Chrysanthemum morifolium

The diversity of form of the chrysanthemum flower makes this species an ideal model for studying petal morphogenesis, but as yet, the molecular mechanisms underlying petal shape development remain largely unexplored. Here, a floral mutant, which arose as a bud sport in a plant of the variety 'Anastasia Dark Green', and formed straight, rather than hooked petals, was subjected to both comparative morphological analysis and transcriptome profiling. The hooked petals only became discernible during a late stage of flower development. At the late stage of 'Anastasia Dark Green', genes related to chloroplast, hormone metabolism, cell wall and microtubules were active, as were cell division-promoting factors. Auxin concentration was significantly reduced, and a positive regulator of cell expansion was down-regulated. Two types of critical candidates, boundary genes and adaxial-abaxial regulators, were identified from 7937 differentially expressed genes in pairwise comparisons, which were up-regulated at the late stage in 'Anastasia Dark Green' and another two hooked varieties. Ectopic expression of a candidate abaxial gene, CmYAB1, in chrysanthemum led to changes in petal curvature and inflorescence morphology. Our findings provide new insights into the regulatory networks underlying chrysanthemum petal morphogenesis.

Current achievements and future prospects in the genetic breeding of chrysanthemum: a review

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a leading flower with applied value worldwide. Developing new chrysanthemum cultivars with novel characteristics such as new flower colors and shapes, plant architectures, flowering times, postharvest quality, and biotic and abiotic stress tolerance in a time- and cost-efficient manner is the ultimate goal for breeders. Various breeding strategies have been employed to improve the aforementioned traits, ranging from conventional techniques, including crossbreeding and mutation breeding, to a series of molecular breeding methods, including transgenic technology, genome editing, and marker-assisted selection (MAS). In addition, the recent extensive advances in high-throughput technologies, especially genomics, transcriptomics, proteomics, metabolomics, and microbiomics, which are collectively referred to as omics platforms, have led to the collection of substantial amounts of data. Integration of these omics data with phenotypic information will enable the identification of genes/pathways responsible for important traits. Several attempts have been made to use emerging molecular and omics methods with the aim of accelerating the breeding of chrysanthemum. However, applying the findings of such studies to practical chrysanthemum breeding remains a considerable challenge, primarily due to the high heterozygosity and polyploidy of the species. This review summarizes the recent achievements in conventional and modern molecular breeding methods and emerging omics technologies and discusses their future applications for improving the agronomic and horticultural characteristics of chrysanthemum.

The Heterologous Expression of a Chrysanthemum nankingense TCP Transcription Factor Blocks Cell Division in Yeast and Arabidopsis thaliana

Both the presence of, and the important contribution to growth and development made by TCP transcription factors, have been established in various plant species. Here, a TCP4 homolog isolated from Chrysanthemum nankingense was shown to be more strongly transcribed in the diploid than in the autotetraploid form of the species. CnTCP4 was shown to encode a member of the class II TCP family and to be transcribed most strongly in the leaf and ligulate flowers. Its transcription was found to be substantially inhibited by spraying the plant with the synthetic cytokinin 6-benzylaminopurine. The transient expression of CnTCP4 in onion epidermal cells showed that its product localized to the nucleus, and a yeast one hybrid assay suggested that its product had transcriptional activation ability. The constitutive expression of CnTCP4 in fission yeast suppressed cell proliferation, inducing the formation of longer and a higher frequency of multinuclated cells. Its constitutive expression in Arabidopsis thaliana reduced the size of the leaves. The presence of the transgene altered the transcription of a number of cell division-related genes. A yeast one hybrid assay identified a second TCP gene (CnTCP2) able to interact with the CnTCP4 promoter. A transient expression experiment in Nicotiana benthamiana leaves showed that CnTCP2 was able to activate the CnTCP4 promoter. Like CnTCP4, CnTCP2 was shown to encode a member of the class II TCP family, to be transcribed most strongly in the leaf and ligulate flowers, and to be suppressed by exogenous 6-benzylaminopurine treatment. The CnTCP2 protein also localized to the nucleus, but had no transcriptional activation ability. Its constitutive expression in A. thaliana had similar phenotypic consequences to those induced by CnTCP4.

Transcriptome analysis during floral organ development provides insights into stamen petaloidy in Lagerstroemia speciosa

As one of the most popular woody species that blooms in summer, Lagerstroemia speciosa has been used abundantly in urban landscape for its excellent floral beauty. For the first time, we discovered a double-flower variant with all petaloid stamens. To understand the molecular basis of this variation, we contrasted the transcriptomes of single- and double-flower buds at three stamen development stages. In total, 73,536 unigenes were mapped and 30,714 differently expressed genes (DEGs) were identified in the tissues. We focused on the DEGs expressing in both phenotypes and investigated the association of their expression profiles with their functions in transcription pathways. Furthermore, we performed WGCNA and identified co-expressed genes with four floral homeotic genes as hubs (MADS16, Unigene0026169; AP2, Unigene0042732; SOC1, Unigene0046314; AG, Unigene0056437). The expression of these hub genes has been conserved across the three developmental stages but significantly different between the two floral phenotypes. As a result, the robust transcriptional regulation of stamen petaloidy in double flowers was deduced. These findings will help to unravel the regulatory mechanisms of several specific genes, thereby providing a basis to study double-flower molecular breeding in L. speciosa.

Comparative transcriptome and metabolite profiling of four tissues from Alisma orientale (Sam.) Juzep reveals its inflorescence developmental and medicinal characteristics

Alisma orientale (Sam.) Juzep (A. orientale) is an important medicinal plant in traditional Chinese medicine. In this study, de novo RNA-seq of A. orientale was performed based on the cDNA libraries from four different tissues, roots, leaves, scapes and inflorescences. A total of 41,685 unigenes were assembled, 25,024 unigene functional annotations were obtained by searching against the five public sequence databases, and 3,411 simple sequence repeats in A. orientale were reported for the first time. 15,402 differentially expressed genes were analysed. The morphological characteristics showed that compared to the other tissues, the leaves had more chlorophyll, the scapes had more vascular bundles, and the inflorescences contained more starch granules and protein. In addition, the metabolic profiles of eight kinds of alisols metabolite profiling, which were measured by ultra-Performance liquid chromatography-triple quadrupole-mass spectrometry showed that alisol B 23-acetate and alisol B were the major components of the four tissues at amounts of 0.068~0.350 mg/g and 0.046~0.587 mg/g, respectively. In addition, qRT-PCR validated that farnesyl pyrophosphate synthase and 3-hydroxy-3-methylglutaryl-CoA reductase should be considered the critical candidate genes involved in alisol biosynthesis. These transcriptome and metabolic profiles of A. orientale may help clarify the molecular mechanisms underlying the medicinal characteristics of A. orientale.

The CmTCP20 gene regulates petal elongation growth in Chrysanthemum morifolium

Chrysanthemum morifolium is one of the most popular ornamental species worldwide, with high ornamental and economic value. Petal size is an important factor that influences the ornamental value. CmTCP20 is a member of TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTORs (TCPs) gene family, which is closely associated with the growth and development of plants. Our previous study found that the expression of CmTCP20 was obviously down-regulated during chrysanthemum petal elongation, but its function in petal elongation has not yet been revealed. We show here that the overexpression CmTCP20 in Arabidopsis and chrysanthemum leads to similar phenotypes, including larger flower buds (or inflorescences) and longer petals. Interestingly, ectopic expression in Schizosaccharomyces pombe yeast cells showed that CmTCP20 could repress cell division and promote cell elongation. Moreover, the yeast two-hybrid, BiFC and pull-down experimental results indicated that CmTCP20 may regulate petal size via interacting with CmJAZ1-like and inducing down-regulation of CmBPE2 gene expression. This study preliminarily clarifies the function of CmTCP20 on chrysanthemum petal elongation, providing the basic theory for improving the ornamental characteristic of chrysanthemum.

TCP Transcription Factors in Moso Bamboo (Phyllostachys edulis): Genome-Wide Identification and Expression Analysis

TEOSINTE BRANCHED 1, CYCLOIDEA, and PROLIFERATING CELL FACTORS (T), members of a plant-specific gene family, play significant roles during plant growth and development, as well as in response to environmental stress. However, knowledge about this family in moso bamboo (Phyllostachys edulis) is limited. Therefore, in this study, the first genome-wide identification, classification, characterization, and expression pattern analysis of the TCP transcription factor family in moso bamboo was performed. Sixteen TCP members were identified from the moso bamboo genome using a BLASTP algorithm-based method and verified using the Pfam database. Based on a multiple-sequence alignment, the members were divided into two subfamilies, and members of the same family shared highly conserved motif structures. Subcellular localization and transactivation activity analyses of four selected genes revealed that they were nuclear localized and had self-activation activities. Additionally, the expression levels of several PeTCP members were significantly upregulated under abscisic acid, methyl jasmonate, and salicylic acid treatments, indicating that they play crucial plant hormone transduction roles in the processes of plant growth and development, as well as in responses to environmental stresses. Thus, the current study provides previously lacking information on the TCP family in moso bamboo and reveals the potential functions of this gene family in growth and development.

Genome engineering in ornamental plants: Current status and future prospects

Ornamental plants, like roses, carnations, and chrysanthemums, are economically important and are sold all over the world. In addition, numerous cut and garden flowers add colors to homes and gardens. Various strategies of plant breeding have been employed to improve traits of many ornamental plants. These approaches span from conventional techniques, such as crossbreeding and mutation breeding, to genetically modified plants. Recently, genome editing has become available as an efficient means for modifying traits in plant species. Genome editing technology is useful for genetic analysis and is poised to become a common breeding method for ornamental plants. In this review, we summarize the benefits and limitations of conventional breeding techniques and genome editing methods and discuss their future potential to accelerate the rate breeding programs in ornamental plants.

Transcriptome analysis of differentially expressed unigenes involved in flavonoid biosynthesis during flower development of Chrysanthemum morifolium 'Chuju'

Chrysanthemum morifolium is an ornamentally and medicinally important plant species. Up to date, molecular and genetic investigations have largely focused on determination of flowering time in the ornamental species. However, little is known about gene regulatory networks for the biosynthesis of flavonoids in the medicinal species. In the current study, we employed the high-throughput sequencing technology to profile the genome-wide transcriptome of C. morifolium 'Chuju', a famous medicinal species in traditional Chinese medicine. A total of 63,854 unigenes with an average length of 741 bp were obtained. Bioinformatic analysis has identified a great number of structural and regulatory unigenes potentially participating in the flavonoid biosynthetic pathway. According to the comparison of digital gene expression, 8,370 (3,026 up-regulated and 5,344 down-regulated), 1,348 (717 up-regulated and 631 down-regulated) and 944 (206 up-regulated and 738 down-regulated) differentially expressed unigenes (DEUs) were detected in the early, middle and mature growth phases, respectively. Among them, many DEUs were implicated in controlling the biosynthesis and composition of flavonoids from the budding to full blooming stages during flower development. Furthermore, the expression patterns of 12 unigenes involved in flavonoid biosynthesis were generally validated by using quantitative real time PCR. These findings could shed light on the molecular basis of flavonoid biosynthesis in C. morifolium 'Chuju' and provide a genetic resource for breeding varieties with improved nutritional quality.

Transcriptomic analysis identifies the key genes involved in stamen petaloid in lotus (Nelumbo nucifera)

Background

Flower morphology, a phenomenon regulated by a complex network, is one of the vital ornamental features in Nelumbo nucifera. Stamen petaloid is very prevalent in lotus flowers. However, the mechanism underlying this phenomenon is still obscure.

Results

Here, the comparative transcriptomic analysis was performed among petal, stamen petaloid and stamen through RNA-seq. Using pairwise comparison analysis, a large number of genes involved in hormonal signal transduction pathways and transcription factors, especially the MADS-box genes, were identified as candidate genes for stamen petaloid in lotus.

Conclusions

Taken together, these results provide an insight into the molecular networks underlying lotus floral organ development and stamen petaloid.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4950-0) contains supplementary material, which is available to authorized users.

De novo transcriptome analysis of Rhododendron molle G. Don flowers by Illumina sequencing

Rhododendron molle G. Don occupies an important phylogenetic node in the genus rhododendron with unique yellow flower and medicinal functions. However, only limited genetic resources and their genome information are available for the generation of rhododendron flowers. The next generation sequencing technologies enables generation of genomic resources in a short time and at a minimal cost, and therefore provide a turning point for rhododendron research. Our goal is to use the genetic information to facilitate the relevant research on flowering and flower color formation in R. molle. In total, 66,026 unigenes were identified, among which 31,298 were annotated in the NCBI non-redundant protein database and 22,410 were annotated in the Swiss-Prot database. Of these annotated unigenes, 9490 and 18,680 unigenes were assigned to clusters of orthologous groups and gene ontology categories, respectively. A total of 7177 genes were mapped to 118 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. In addition, 8266 simple sequence repeats (SSRs) were detected, and these SSRs will undoubtedly benefit rhododendron breeding work. Metabolic pathway analysis revealed that 32 unigenes were predicted to be involved in carotenoid biosynthesis. Our transcriptome revealed 32 engines that encode key enzymes in the carotenoid biosynthesis pathway, including PSY, PDS, LCYB, LCYE, etc. The content of β-carotene was much higher than the other carotenoids throughout the flower development. It was consistent with the key genes expression level in the carotenoid biosynthesis pathway by the Illumina expression profile analysis and the qRT-PCR analysis. Our study identified genes associated with carotenoid biosynthesis in R. molle and provides a valuable resource for understanding the flowering and flower color formation mechanisms in R. molle.

A Mini Zinc-Finger Protein (MIF) from Gerbera hybrida Activates the GASA Protein Family Gene, GEG, to Inhibit Ray Petal Elongation

Petal appearance is an important horticultural trail that is generally used to evaluate the ornamental value of plants. However, knowledge of the molecular regulation of petal growth is mostly derived from analyses of Arabidopsis thaliana, and relatively little is known about this process in ornamental plants. Previously, GEG (Gerbera hybrida homolog of the gibberellin [GA]-stimulated transcript 1 [GAST1] from tomato), a gene from the GA stimulated Arabidopsis (GASA) family, was reported to be an inhibitor of ray petal growth in the ornamental species, G. hybrida. To explore the molecular regulatory mechanism of GEG in petal growth inhibition, a mini zinc-finger protein (MIF) was identified using yeast one-hybrid (Y1H) screen. The direct binding of GhMIF to the GEG promoter was verified by using an electrophoretic mobility shift assay and a dual-luciferase assay. A yeast two-hybrid (Y2H) revealed that GhMIF acts as a transcriptional activator. Transient transformation assay indicated that GhMIF is involved in inhibiting ray petal elongation by activating the expression of GEG. Spatiotemporal expression analyses and hormone treatment assay showed that the expression of GhMIF and GEG is coordinated during petal development. Taken together, these results suggest that GhMIF acts as a direct transcriptional activator of GEG, a gene from the GASA protein family to regulate the petal elongation.