Cell layer-specific expression of the homeotic MADS-box transcription factor PhDEF contributes to modular petal morphogenesis in petunia.

Floral homeotic MADS-box transcription factors ensure the correct morphogenesis of floral organs, which are organized in different cell layers deriving from distinct meristematic layers. How cells from these distinct layers acquire their respective identities and coordinate their growth to ensure normal floral organ morphogenesis is unresolved. Here, we studied petunia (Petunia × hybrida) petals that form a limb and tube through congenital fusion. We identified petunia mutants (periclinal chimeras) expressing the B-class MADS-box gene DEFICIENS in the petal epidermis or in the petal mesophyll, called wico and star, respectively. Strikingly, wico flowers form a strongly reduced tube while their limbs are almost normal, while star flowers form a normal tube but greatly reduced and unpigmented limbs, showing that petunia petal morphogenesis is highly modular. These mutants highlight the layer-specific roles of PhDEF during petal development. We explored the link between PhDEF and petal pigmentation, a well-characterized limb epidermal trait. The anthocyanin biosynthesis pathway was strongly downregulated in star petals, including its major regulator ANTHOCYANIN2 (AN2). We established that PhDEF directly binds to the AN2 terminator in vitro and in vivo, suggesting that PhDEF might regulate AN2 expression and therefore petal epidermis pigmentation. Altogether, we show that cell layer-specific homeotic activity in petunia petals differently impacts tube and limb development, revealing the relative importance of the different cell layers in the modular architecture of petunia petals.

The R2R3-MYB transcription factor EVER controls the emission of petunia floral volatiles by regulating epicuticular wax biosynthesis in the petal epidermis.

The epidermal cells of petunia (Petunia × hybrida) flowers are the main site of volatile emission. However, the mechanisms underlying the release of volatiles into the environment are still being explored. Here, using cell-layer-specific transcriptomic analysis, reverse genetics by virus-induced gene silencing and clustered regularly interspaced short palindromic repeat (CRISPR), and metabolomics, we identified EPIDERMIS VOLATILE EMISSION REGULATOR (EVER)-a petal adaxial epidermis-specific MYB activator that affects the emission of volatiles. To generate ever knockout lines, we developed a viral-based CRISPR/Cas9 system for efficient gene editing in plants. These knockout lines, together with transient-suppression assays, revealed EVER's involvement in the repression of low-vapor-pressure volatiles. Internal pools and annotated scent-related genes involved in volatile production and emission were not affected by EVER. RNA-Seq analyses of petals of ever knockout lines and EVER-overexpressing flowers revealed enrichment in wax-related biosynthesis genes. Liquid chromatography/gas chromatography-MS analyses of petal epicuticular waxes revealed substantial reductions in wax loads in ever petals, particularly of monomers of fatty acids and wax esters. These results implicate EVER in the emission of volatiles by fine-tuning the composition of petal epicuticular waxes. We reveal a petunia MYB regulator that interlinks epicuticular wax composition and volatile emission, thus unraveling a regulatory layer in the scent-emission machinery in petunia flowers.

Transcriptomic analysis implicates ABA signaling and carbon supply in the differential outgrowth of petunia axillary buds.

BACKGROUND: Shoot branching of flowering plants exhibits phenotypic plasticity and variability. This plasticity is determined by the activity of axillary meristems, which in turn is influenced by endogenous and exogenous cues such as nutrients and light. In many species, not all buds on the main shoot develop into branches despite favorable growing conditions. In petunia, basal axillary buds (buds 1-3) typically do not grow out to form branches, while more apical axillary buds (buds 6 and 7) are competent to grow. RESULTS: The genetic regulation of buds was explored using transcriptome analyses of petunia axillary buds at different positions on the main stem. To suppress or promote bud outgrowth, we grew the plants in media with differing phosphate (P) levels. Using RNA-seq, we found many (> 5000) differentially expressed genes between bud 6 or 7, and bud 2. In addition, more genes were differentially expressed when we transferred the plants from low P to high P medium, compared with shifting from high P to low P medium. Buds 6 and 7 had increased transcript abundance of cytokinin and auxin-related genes, whereas the basal non-growing buds (bud 2 and to a lesser extent bud 3) had higher expression of strigolactone, abscisic acid, and dormancy-related genes, suggesting the outgrowth of these basal buds was actively suppressed. Consistent with this, the expression of ABA associated genes decreased significantly in apical buds after stimulating growth by switching the medium from low P to high P. Furthermore, comparisons between our data and transcriptome data from other species suggest that the suppression of outgrowth of bud 2 was correlated with a limited supply of carbon to these axillary buds. Candidate genes that might repress bud outgrowth were identified by co-expression analysis. CONCLUSIONS: Plants need to balance growth of axillary buds into branches to fit with available resources while allowing some buds to remain dormant to grow after the loss of plant parts or in response to a change in environmental conditions. Here we demonstrate that different buds on the same plant with different developmental potentials have quite different transcriptome profiles.

Evolutionary walks through flower colour space driven by gene expression in Petunia and allies (Petunieae).

The structure and function of biochemical and developmental pathways determine the range of accessible phenotypes, which are the substrate for evolutionary change. Accordingly, we expect that observed phenotypic variation across species is strongly influenced by pathway structure, with different phenotypes arising due to changes in activity along pathway branches. Here, we use flower colour as a model to investigate how the structure of pigment pathways shapes the evolution of phenotypic diversity. We focus on the phenotypically diverse Petunieae clade in the nightshade family, which contains ca 180 species of Petunia and related genera, as a model to understand how flavonoid pathway gene expression maps onto pigment production. We use multivariate comparative methods to estimate co-expression relationships between pathway enzymes and transcriptional regulators, and then assess how expression of these genes relates to the major axes of variation in floral pigmentation. Our results indicate that coordinated shifts in gene expression predict transitions in both total anthocyanin levels and pigment type, which, in turn, incur trade-offs with the production of UV-absorbing flavonol compounds. These findings demonstrate that the intrinsic structure of the flavonoid pathway and its regulatory architecture underlies the accessibility of pigment phenotypes and shapes evolutionary outcomes for floral pigment production.

Hormonal Signaling during dPCD: Cytokinin as the Determinant of RNase-Based Self-Incompatibility in Solanaceae.

Research into molecular mechanisms of self-incompatibility (SI) in plants can be observed in representatives of various families, including Solanaceae. Earlier studies of the mechanisms of S-RNase-based SI in petunia (Petunia hybrida E. Vilm.) demonstrate that programmed cell death (PCD) is an SI factor. These studies suggest that the phytohormon cytokinin (CK) is putative activator of caspase-like proteases (CLPs). In this work, data confirming this hypothesis were obtained in two model objects-petunia and tomato (six Solanaceae representatives). The exogenous zeatin treatment of tomato and petunia stigmas before a compatible pollination activates CLPs in the pollen tubes in vivo, as shown via the intravital imaging of CLP activities. CK at any concentration slows down the germination and growth of petunia and tomato male gametophytes both in vitro and in vivo; shifts the pH of the cytoplasm (PHc) to the acid region, thereby creating the optimal conditions for CLP to function and inhibiting the F-actin formation and/or destructing the cytoskeleton in pollen tubes to point foci during SI-induced PCD; and accumulates in style tissues during SI response. The activity of the ISOPENTENYLTRANSFERASE 5 (IPT5) gene at this moment exceeds its activity in a cross-compatible pollination, and the levels of expression of the CKX1 and CKX2 genes (CK OXIDASE/DEHYDROGENASE) are significantly lower in self-incompatible pollination. All this suggests that CK plays a decisive role in the mechanism underlying SI-induced PCD.

Functional identification of PhMATE1 in flower color formation in petunia.

Multidrug and toxic compound extrusion (MATE) transporter proteins are a class of secondary transporter proteins that can transport flavonoids. Anthocyanins, a kind of flavonoid, are important secondary metabolites widely found in higher plants; they determine the flower color of most angiosperms. TT12 in Arabidopsis was the first MATE protein identified to be involved in flavonoid transport. Petunia (Petunia hybrida) is an important ornamental plant and is one of the ideal plants for studying plant flower color. However, there are few reports on anthocyanin transport in petunia. In this study, we characterized a homolog of Arabidopsis TT12 in the petunia genome, PhMATE1, that shares the highest amino acid sequence identity with Arabidopsis TT12. PhMATE1 protein contained 11 transmembrane helices. PhMATE1 showed a high transcription level in corollas. The silencing of PhMATE1 mediated by both virus-induced gene silence and RNA interference changed flower color and reduced anthocyanin content in petunia, suggesting that PhMATE1 is involved in anthocyanin transport in petunia. Furthermore, PhMATE1 silencing downregulated the expression of the structural genes of the anthocyanin synthesis pathway. The results of this study supported the hypothesis that MATEs are involved in the sequestration of anthocyanins during flower color formation.

Altered profile of floral volatiles and lignin content by down-regulation of Caffeoyl Shikimate Esterase in Petunia.

BACKGROUND: The floral volatile profile of Petunia x hybrida 'Mitchell diploid' (MD) is dominated by phenylpropanoids, many of which are derived from p-coumaric acid. However, the downstream processes involved in the production of caffeoyl-CoA and feruloyl-CoA from p-coumaric acid are complex, as the genes and biosynthesis steps are associated with flavonoids and lignin synthesis as well as floral volatiles benzenoid/phenylpropanoid (FVBP). Caffeoyl shikimate esterase (CSE) converts caffeoyl shikimate to caffeic acid and is considered one of the essential regulators in lignin production. Moreover, CSE in involved in phenylpropanoid production. To investigate the roles of CSE in FVBP biosynthesis, we used RNAi-mediated CSE down-regulated (ir-PhCSE) petunias. RESULTS: Lowered CSE transcript accumulation in ir-PhCSE plants resulted in reduced lignin layers in the stems and stunted growth, suggesting a positive correlation between lignin layers and lignin content. The altered CSE level influenced the expression of many FVBP genes, including elevated transcripts of p-coumarate-3-hydroxylase (C3H), hydroxycinnamoyl transferase (HCT), and 4-coumaric acid: CoA ligase (4CL). In particular, the expression of C4H in ir-PhCSE plants was more than twice the expression in MD plants. Moreover, the production of volatile compounds was alterend in ir-PhCSE plants. Most floral volatiles decreased, and the amounts of phenylalanine and caffeic acid were significantly lower. CONCLUSIONS: Reduced lignin layers in the stems and stunted growth in ir-PhCSE plants suggest that PhCSE is essential for lignin production and plant growth in petunia. The decreased CSE level influenced the expression of many FVBP genes, and interference of shikimate derivates altered volatile compound production. Significantly decreased caffeic acid, but not ferulic acid, in ir-PhCSE plants suggest that CSE is primarily involved in the reaction of caffeoyl shikimate. Higher C3H and C4H transcripts seem to alleviate accumulated p-coumaric acid resulting from altered CSE. Finally, alteration in C3H, HCT, and 4CL in CSE down-regulated plants suggests an interaction of the FVBP genes, leading to the regulation of floral volatiles of petunia.

Genetic architecture of a pollinator shift and its fate in secondary hybrid zones of two Petunia species.

BACKGROUND: Theory suggests that the genetic architecture of traits under divergent natural selection influences how easily reproductive barriers evolve and are maintained between species. Divergently selected traits with a simple genetic architecture (few loci with major phenotypic effects) should facilitate the establishment and maintenance of reproductive isolation between species that are still connected by some gene flow. While empirical support for this idea appears to be mixed, most studies test the influence of trait architectures on reproductive isolation only indirectly. Petunia plant species are, in part, reproductively isolated by their different pollinators. To investigate the genetic causes and consequences of this ecological isolation, we deciphered the genetic architecture of three floral pollination syndrome traits in naturally occurring hybrids between the widespread Petunia axillaris and the highly endemic and endangered P. exserta. RESULTS: Using population genetics, Bayesian linear mixed modelling and genome-wide association studies, we found that the three pollination syndrome traits vary in genetic architecture. Few genome regions explain a majority of the variation in flavonol content (defining UV floral colour) and strongly predict the trait value in hybrids irrespective of interspecific admixture in the rest of their genomes. In contrast, variation in pistil exsertion and anthocyanin content (defining visible floral colour) is controlled by many genome-wide loci. Opposite to flavonol content, the genome-wide proportion of admixture between the two species predicts trait values in their hybrids. Finally, the genome regions strongly associated with the traits do not show extreme divergence between individuals representing the two species, suggesting that divergent selection on these genome regions is relatively weak within their contact zones. CONCLUSIONS: Among the traits analysed, those with a more complex genetic architecture are best maintained in association with the species upon their secondary contact. We propose that this maintained genotype-phenotype association is a coincidental consequence of the complex genetic architectures of these traits: some of their many underlying small-effect loci are likely to be coincidentally linked with the actual barrier loci keeping these species partially isolated upon secondary contact. Hence, the genetic architecture of a trait seems to matter for the outcome of hybridization not only then when the trait itself is under selection.

SCARECROW-like GRAS protein PES positively regulates petunia floral scent production.

Emission of scent volatiles by flowers is important for successful pollination and consequently, reproduction. Petunia (Petunia hybrida) floral scent is formed mainly by volatile products of the phenylpropanoid pathway. We identified and characterized a regulator of petunia scent production: the GRAS protein PHENYLPROPANOID EMISSION-REGULATING SCARECROW-LIKE (PES). Its expression increased in petals during bud development and was highest in open flowers. Overexpression of PES increased the production of floral volatiles, while its suppression resulted in scent reduction. We showed that PES upregulates the expression of genes encoding enzymes of the phenylpropanoid and shikimate pathways in petals, and of the core regulator of volatile biosynthesis ODORANT1 by activating its promoter. PES is an ortholog of Arabidopsis (Arabidopsis thaliana) PHYTOCHROME A SIGNAL TRANSDUCTION 1, involved in physiological responses to far-red (FR) light. Analyses of the effect of nonphotosynthetic irradiation (low-intensity FR light) on petunia floral volatiles revealed FR light as a scent-activating factor. While PHYTOCHROME A regulated scent-related gene expression and floral scent production under FR light, the influence of PES on volatile production was not limited by FR light conditions.

PhAAT1, encoding an anthocyanin acyltransferase, is transcriptionally regulated by PhAN2 in petunia.

Anthocyanins widely exist in plants and they are important pigments for color of petals and fruits. They are produced through a multi-step pathway controlled by transcription factor complexes. The anthocyanin skeleton modification is the last reaction in the anthocyanin synthesis pathway, which improves the stability of anthocyanins. Acylation modification is an important modification of anthocyanins. However, the identification and function of anthocyanin acyltransferase genes and their expression regulation are rarely reported. In this study, we identified the petunia anthocyanin acyltransferase gene, PhAAT1. PhAAT1 is located in the cytoplasm and PhAAT1 silencing changed flower color and reduced the stability of anthocyanin. Metabolomics analysis showed that PhAAT1 silencing led to the reduction of p-coumaroylated and caffeoylated anthocyanins. In addition, PhAAT1 was positively regulated by the MYB transcription factor, PhAN2, which directly interacts with the promoter of PhAAT1.

Physiological changes besides the enhancement of pigmentation in Petunia hybrida caused by overexpression of PhAN2, an R2R3-MYB transcription factor.

KEY MESSAGE: Ectopic expression of PhAN2 in vegetative tissue can improve regeneration and adventitious rooting but inhibit axillary bud outgrowth of petunia, while overexpression specifically in flowers could shorten longevity. Anthocyanin 2 has been only treated as a critical positive regulation factor of anthocyanin biosynthesis in petunia flowers. To determine if this gene had other functions in plant growth, we overexpressed this gene in an an2 mutant petunia cultivar driven by promoters with different strengths or tissue specificity. Various physiological processes of transformants in different growth stages and environments were analyzed. Besides the expected pigmentation improvement in different tissues, the results also showed that ectopic expression of AN2 could improve the regeneration skill but inhibit the axillary bud germination of in vitro plants. Moreover, the rooting ability of shoot tips of transformants was significantly improved, while some transgenic lines' flower longevity was shortened. Gene expression analysis showed that the transcripts level of AN2, partner genes anthocyanin 1 (AN1), anthocyanin 11 (AN11), and target gene dihydroflavonol 4-reductase (DFR) was altered in the different transgenic lines. In addition, ethylene biosynthesis-related genes 1-aminocyclopropane-1-carboxylic acid synthase (ACS1) and ACC oxidase (ACO1) were upregulated in rooting and flower senescence processes but at different time points. Overall, our data demonstrate that the critical role of this AN2 gene in plant growth physiology may extend beyond that of a single activator of anthocyanin biosynthesis.

Protoplast isolation and transient gene expression in different petunia cultivars.

The protocol optimized for Petunia hybrida cv. Mirage Rose produced high protoplast yields in 3 out of other 11 cultivars (Damask White, Dreams White, and Opera Supreme White). Factors optimized in the protoplast transfection process showed that the best transfection efficiency (80%) was obtained using 2.5 × 105 protoplast density, 40% polyethylene glycol (PEG) concentration, 10 µg plasmid DNA, and 15 min of transfection time. Assessing the usability of the protocol for other cultivars (Damask White, Dreams White, and Opera Supreme White), a reasonable protoplast transfection efficiency (⁓50%) was observed in the cultivars Dreams White and Opera Supreme White, with lower efficiency (⁓50%) observed in the cv. Damask White. The transient expression of enhanced green fluorescent protein (eGFP) in the nucleus of the transfected protoplasts of all cultivars was confirmed using PCR. This system could be valuable for genome editing of unwanted genes in petunias using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated 9 (Cas9) technology. Furthermore, it could contribute to other studies on protein subcellular localization, protein-protein interactions, and functional gene expression in the petunias.

Shikimate Kinase Plays Important Roles in Anthocyanin Synthesis in Petunia.

In plants, the shikimate pathway is responsible for the production of aromatic amino acids L-tryptophan, L-phenylalanine, and L-tyrosine. L-Phenylalanine is the upstream substrate of flavonoid and anthocyanin synthesis. Shikimate kinase (SK) catalyzes the phosphorylation of the C3 hydroxyl group of shikimate to produce 3-phosphate shikimate (S3P), the fifth step of the shikimate pathway. However, whether SK participates in flavonoid and anthocyanin synthesis is unknown. This study characterized the single-copy PhSK gene in the petunia (Petunia hybrida) genome. PhSK was localized in chloroplasts. PhSK showed a high transcription level in corollas, especially in the coloring stage of flower buds. Suppression of PhSK changed flower color and shape, reduced the content of anthocyanins, and changed the flavonoid metabolome profile in petunia. Surprisingly, PhSK silencing caused a reduction in the shikimate, a substrate of PhSK. Further qPCR analysis showed that PhSK silencing resulted in a reduction in the mRNA level of PhDHQ/SDH, which encodes the protein catalyzing the third and fourth steps of the shikimate pathway, showing a feedback regulation mechanism of gene expression in the shikimate pathway.

PhMYB37 Promotes Shoot Branching in Petunia.

Petunia is one of the world's most important flowers, and its branch development has long been a source of discussion. MYB transcription factors have been identified as important plant branching regulators. In this study, 113 R2R3-MYB genes were identified from the petunia genome. PhMYB genes, closely related to RAXs, were expressed at greater levels in axillary buds and roots. Decapitation and 6-BA did not regulate the expression of PhMYB37. PhMYB37 was localized in the nucleus. Heterologous overexpression of PhMYB37 promoted shoot branching in transgenic Arabidopsis while silencing of PhMYB37 inhibited shoot branching. These results suggest that PhMYB37 plays a critical and positive role in petunia shoot branching.

Crotonylation versus acetylation in petunia corollas with reduced acetyl-CoA due to PaACL silencing.

Protein acetylation and crotonylation are important posttranslational modifications of lysine. In animal cells, the correlation of acetylation and crotonylation has been well characterized and the lysines of some proteins are acetylated or crotonylated depending on the relative concentrations of acetyl-CoA and crotonyl-CoA. However, in plants, the correlation of acetylation and crotonylation and the effects of the relative intracellular concentrations of crotonyl-CoA and acetyl-CoA on protein crotonylation and acetylation are not well known. In our previous study, PaACL silencing changed the content of acetyl-CoA in petunia (Petunia hybrida) corollas, and the effect of PaACL silencing on the global acetylation proteome in petunia was analyzed. In the present study, we found that PaACL silencing did not significantly alter the content of crotonyl-CoA. We performed a global crotonylation proteome analysis of the corollas of PaACL-silenced and control petunia plants; we found that protein crotonylation was closely related to protein acetylation and that proteins with more crotonylation sites often had more acetylation sites. Crotonylated proteins and acetylated proteins were enriched in many common Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. However, PaACL silencing resulted in different KEGG pathway enrichments of proteins with different levels of crotonylation sites and acetylation sites. PaACLB1-B2 silencing did not led to changes in the opposite direction in crotonylation and acetylation levels at the same lysine site in cytoplasmic proteins, which indicated that cytoplasmic lysine acetylation and crotonylation might not depend on the relative concentrations of acetyl-CoA and crotonyl-CoA. Moreover, the global crotonylome and acetylome were weakly positively correlated in the corollas of PaACL-silenced and control plants.

Ph5GT silencing alters flower color and flavonoids metabolome profile in petunia.

Anthocyanins are important pigments in plants and glycosylation plays an important role in the stability of anthocyanins. Anthocyanin 5-O-glucosyltransferase (5GT) can glycosylate anthocyanin at the 5-O-position. Till now, the enzymatic activity characteristics of 5GT had been studied in vitro in a variety of plants. However, the subcellular localization of 5GT protein still remained unclear, and little genetic evidence on the roles of 5GT in plants has been reported. The full-length Ph5GT gene from petunia (Petunia hybrida) was isolated in this study. Green fluorescent fusion protein assays revealed that Ph5GT protein was localized to the cytoplasm. Ph5GT was found to be highly expressed in flowers, with highest levels of expression occurring during the coloring stage of flower development. Furthermore, Ph5GT silencing led to the change in flower color from purple to light purple and a significant reduction in total anthocyanin content. The metabolome analysis revealed that the content of malvidins and petunidins modified by glycosylation at the 5-O-position was significantly reduced, while the content of their precursor without glycosylation was significantly increased, implying that Ph5GT could glycosylate malvidin and petunidin derivatives and that the substrate types of Ph5GT were expanded in comparison to previous studies.

Overexpression of acdS in Petunia hybrida Improved Flower Longevity and Cadmium-Stress Tolerance by Reducing Ethylene Production in Floral and Vegetative Tissues.

The role of acdS, which encodes the 1-aminocyclopropane-1-carboxylic acid (ACC) deaminase enzyme, in extending flower longevity and improving tolerance to cadmium (Cd) stress was assessed using transgenic Petunia hybrida cv. 'Mirage Rose' overexpressing acdS and wild-type (WT) plants. The overexpression of acdS reduced ethylene production in floral tissue via suppression of ethylene-related genes and improved flower longevity, approximately 2 to 4 days longer than WT flowers. Under Cd stress, acdS significantly reduced Cd-induced ethylene production in vegetable tissues of transgenic plants through suppression of ethylene-related genes. This resulted in a lower accumulation of ethylene-induced reactive oxygen species (ROS) in the transgenic plants than in WT plants. In addition, expression of the genes involved in the activities of antioxidant and proline synthesis as well as the metal chelation process was also higher in the former than in the latter. Moreover, Cd accumulation was significantly higher in WT plants than in the transgenic plants. These results are linked to the greater tolerance of transgenic plants to Cd stress than the WT plants, which was determined based on plant growth and physiological performance. These results highlight the potential applicability of using acdS to extend flower longevity of ornamental bedding plants and also reveal the mechanism by which acdS improves Cd-stress tolerance. We suggest that acdS overexpression in plants can extend flower longevity and also help reduce the negative impact of Cd-induced ethylene on plant growth when the plants are unavoidably cultivated in Cd-contaminated soil.

PhRHMs play important roles in leaf and flower development and anthocyanin synthesis in petunia.

Anthocyanins, vital metabolites in plants, are formed by anthocyanidins combined with various monosaccharides, including glucose, rhamnose, and arabinose. Rhamnose contributes greatly to the glycosylation of anthocyanidins. There are two kinds of rhamnose synthase (RS): rhamnose biosynthesis (RHM), and nucleotide-RS/epimerase-reductase (UER1). Nevertheless, no RS isoform was reported to be involved in anthocyanin synthesis. Here, three homologous PhRHM genes, namely PhRHM1, PhRHM2, and PhRHM3, and one PhUER1 gene from petunia were cloned and characterized. Green fluorescent protein fusion protein assays revealed that PhRHMs and PhUER1 are localized in the cytoplasm. We obtained PhRHM1 or/and PhRHM2 or PhUER1 silenced petunia plants and did not attempt to obtain PhRHM3 silenced plants since PhRHM3 mRNA was not detected in petunia organs examined. PhRHM1 and PhRHM2 (PhRHM1-2) silencing induced abnormal plant growth and decreased the contents of l-rhamnose, photosynthetic pigments and total anthocyanins, while PhUER1 silencing did not cause any visible phenotypic changes. Flavonoid metabolome analysis further revealed that PhRHM1-2 silencing reduced the contents of anthocyanins with rhamnose residue. These results revealed that PhRHMs contribute to the biosynthesis of rhamnose and that PhRHMs participate in the anthocyanin rhamnosylation in petunia, while PhUER1 does not.