Transcriptomic Analysis Suggests Auxin Regulation in Dorsal-Ventral Petal Asymmetry of Wild Progenitor Sinningia speciosa

The establishment of dorsal–ventral (DV) petal asymmetry is accompanied by differential growth of DV petal size, shape, and color differences, which enhance ornamental values. Genes involved in flower symmetry in Sinningia speciosa have been identified as CYCLOIDEA (SsCYC), but which gene regulatory network (GRN) is associated with SsCYC to establish DV petal asymmetry is still unknown. To uncover the GRN of DV petal asymmetry, we identified 630 DV differentially expressed genes (DV-DEGs) from the RNA-Seq of dorsal and ventral petals in the wild progenitor, S. speciosa ‘ES’. Validated by qRT-PCR, genes in the auxin signaling transduction pathway, SsCYC, and a major regulator of anthocyanin biosynthesis were upregulated in dorsal petals. These genes correlated with a higher endogenous auxin level in dorsal petals, with longer tube length growth through cell expansion and a purple dorsal color. Over-expression of SsCYC in Nicotiana reduced petal size by regulating cell growth, suggesting that SsCYC also controls cell expansion. This suggests that auxin and SsCYC both regulate DV petal asymmetry. Transiently over-expressed SsCYC, however, could not activate most major auxin signaling genes, suggesting that SsCYC may not trigger auxin regulation. Whether auxin can activate SsCYC or whether they act independently to regulate DV petal asymmetry remains to be explored in the future.


Introduction
Floral zygomorphy (bilateral symmetry, dorsal-ventral asymmetry) has been shown to be a key innovation that is often associated with flower diversity in many angiosperms due to its role in facilitating plant-pollinator interactions. Zygomorphy has evolved independently in many angiosperm lineages [1]. Within the Lamiales, the flowers of many outcrossing zygomorphic lineages in the Gesneriaceae show diversification in size, shape, and color that are associated with pollination specializations [2][3][4]. Species in the genus Sinningia (subfamily Gesnerioideae, tribe Gesnerieae) are remarkably diverse, with flowers that have features consistent with bee, hummingbird, bat, and moth pollination syndromes [5]. Floral morphologies, composed of various tubular shapes and pigment patterns, are shown to be correlated with different types of pollination syndromes, and further affect flower diversification in Gesneriaceae [5][6][7]. The development of floral dorsalventral asymmetry therefore appears to be associated with flower pollination traits relating to size, shape, and color. Sinningia speciosa (Lodd.) Hiern 'Espírito Santo' ('ES') is a bee-pollinated wild perennial species. The zygomorphic corolla tube of S. speciosa is distinguished by an asymmetric dorsal-ventral petal contour that comprises prominent ventral swelling and dorsal gibbosity (nectar chamber) at the base of the corolla tube ( Figure 1). In addition, the corolla exhibits an asymmetrical dorsal-ventral purple pigmentation pattern, in which dark spots or stripes are present on the inner surface of the ventral corolla tube, which likely serve as a nectar guide for visiting insects. These pollination-related traits such as dorsal-ventral petal size, petal shape, and pigmentation differences in the corolla therefore appear to correlate with the development of floral zygomorphy of S. speciosa.
The key genes responsible for establishing the dorsal-ventral identity of flowers have been shown to be homologs of CYCLOIDEA (CYC) of Antirrhinum majus in Lamiales and almost all major angiosperm lineages [8,9]. CYC and its paralog DICHOTOMA (DICH) activate their downstream target MYB gene, RADIALIS (RAD), in dorsal petals by directly binding to the RAD promoter [10]. In contrast, the ventral identity is established by a MYB protein complex consisting of DIVARICATA (DIV) and DIV-and-RAD-interactingfactor (DRIF), which controls the expression of downstream genes involved in floral ventralization. In dorsal petals, the RAD protein can bind to DRIF and prevent the formation of the DIV/DRIF complex [11]. Therefore, the ventralizing effect of DIV is restricted to the ventral petals, while CYC activates RAD in the dorsal petals, thus creating the dorsal-ventral polarity. In Gesneriaceae species, diversified CYC expression patterns have been found to correlate with floral symmetry transition and flower shape diversity [12][13][14]. However, little is known about which gene regulatory network (GRN) is necessary in addition to CYC-RAD and DIV in generating the size and shape differences between dorsal and ventral petals.
CYC-like genes belong to the TCP gene family encoding plant-specific proteins sharing a basic helix-loop-helix (bHLH) motif, the so-called TCP domain that allows DNA binding and regulates plant growth and development [15]. Several functional studies have shown that CYC-like genes in the class II TCP family are involved in regulating cell proliferation or expansion to affect organ size. Transgenic plants of Arabidopsis thaliana over-expressing the CYC gene from A. majus show increased petal size as a consequence of an increase in cell expansion in petals, whereas CYC reduces both cell proliferation (by promoting the early arrest of cell division) and cell expansion in leaves [16,17]. Among Gesneriaceae species, over-expressing the SiCYC1A gene from the African violet (Streptocarpus ionanthus) in A. thaliana results in a reduction in petal size due to decreased cell proliferation [14]. The constitutive expression of CYC1C homologous genes from Primulina heterotricha and SsCYC from S. speciosa in A. thaliana, however, reduce petal size via the regulation of decreased cell growth [18,19]. Thus, there appear to be no simple rule to predict the effects of CYC homologs on petal growth across angiosperm lineages.
Although flowers of wild progenies of S. speciosa are zygomorphic, many actinomorphic flower cultivars have been favored and selected by humans, attributing to a SsCYC coding region deletion [19]. Efforts in crossing the zygomorphic wildtype (containing functional SsCYC alleles) with the actinomorphic peloria (containing mutated Sscyc alleles) resulted in the F1 generation, which are all zygomorphic [19][20][21]. The F2 generation from F1 selfing are segregated at a 3:1 ratio of dominancy of zygomorphy to actinomorphy, and this allows genotype-phenotype associations of SsCYC alleles to petal morphology [20,21]. Geometric morphometric analysis of 3D petal morphology of S. speciosa has been shown to be useful to analyze flower shape and size variation [20,22]. The associations of SsCYC alleles with 3D petal morphology reveal that the functional SsCYC allele correlates with dorsal petal outward curvature (e.g., the gibbous structure-nectar chamber) and ventral outward curvature (e.g., the ventricose chamber in the ventral tube, see our Figure 1) [21]. This means that SsCYC not only functions to establish flower dorsal-ventral identity, but may also affect petal shape, size, and color differences between dorsal and ventral petals. dorsal-ventral petal contour that comprises prominent ventral swelling and dorsal gibbosity (nectar chamber) at the base of the corolla tube ( Figure 1). In addition, the corolla exhibits an asymmetrical dorsal-ventral purple pigmentation pattern, in which dark spots or stripes are present on the inner surface of the ventral corolla tube, which likely serve as a nectar guide for visiting insects. These pollination-related traits such as dorsal-ventral petal size, petal shape, and pigmentation differences in the corolla therefore appear to correlate with the development of floral zygomorphy of S. speciosa.  (D) Corolla tube development from stages FB3 to FB16, during which the dorsal-ventral petal asymmetry is established gradually. The red and blue lines indicate a gibbous outgrowth at the base of the dorsal corolla tube (red curve), and a ventricose chamber (blue curve) at the ventral corolla tube. Refer to the main text for a detailed explanation. The dashed yellow lines demonstrate the tube length difference along the curvature of the dorsal and ventral corolla tubes. Scale bars = 5 mm. (E,F) The measurement of tube length (E) and cell number (F) between the dorsal and ventral parts of the corolla. (G) A cartoon illustration showing the gradient and spotting pigmentation patterns on the adaxial and abaxial surfaces of the corolla tube at FB16. (H) Purple pigments (anthocyanin) start to accumulate on cells of the adaxial surface of the dorsal petal at FB5 and become evident at FB8. (a,b) Purple spots also appear on the adaxial surface of the ventral petal at FB5 (arrowheads) but become obvious at FB8 (c,d). dPe, dorsal petal; lPe, lateral petal; vPe, ventral petal; dSt, dorsal staminode; lSt, lateral stamen; vSt, ventral stamen; Pi, pistil; Pe, petal; St, stamen primordia.
We therefore would like to identify what developmental modules, in addition to SsCYC, could contribute to dorsal-ventral petal shape and color asymmetry in the wild zygomorphic progeny of Sinningia speciosa 'ES'. We hypothesize that the gene regulatory networks involved in regulating cell growth and the purple color gradient between the dorsal and ventral petals are likely accompanied by the establishment of floral symmetry. There is currently little information about what GRN may be dependent or independent of SsCYC in creating dorsal-ventral petal asymmetry across angiosperms. In this study, we therefore performed RNA-Seq separately on the dorsal and ventral petals during a single floral bud developmental stage, establishing floral zygomorphy when SsCYC is dorsally highly expressed. We first characterized flower buds to ascertain the stages in which the major morphological changes occurred during the establishment of dorsal-ventral petal asymmetry using SEM and petal size measurement. These morphological observations helped to identify differentially expressed genes that relate to dorsal-ventral petal asymmetry. Auxin signaling pathway genes, together with the SsCYC-SsRAD regulatory module and the major regulator of anthocyanin biosynthesis, showed dorsal-ventral asymmetrical expressions. A model illustrating how the regulation of these genes correlates with asymmetric dorsal-ventral petal growth and pigmentation patterns in S. speciosa 'ES' was proposed.

The Development of Zygomorphy Is Due to Shape and Size Differences in the Dorsal and Ventral Parts of the Corolla Tube
To characterize when the floral zygomorphy in S. speciosa 'ES' was initially established, we examined the developing floral primordia throughout floral bud initiation stages using SEM. Floral zygomorphy, however, was not established in the early bud stage until FB1-7 (ca. 0.5 mm in diameter) (Figure 1(Ca-Cc)). At stage FB1-4, each of the five petal primordia and five stamen primordia were equal in size (Figure 1(Ca)). Floral symmetry thus appears to be actinomorphic at the beginning of petal and stamen organogenesis. Later during development, at stages of FB1-5 and FB1-6, all five petal primordia continued to grow at the same rate and started to differentiate into the tube at the base and the petal lobes towards the tip (Figure 1(Cb,Cc). At stage FB1-7, while the tube further elongated, the corolla lobes enlarged to enclose the stamens and pistils. The arrangement of five lobes developed into a specified aestivation pattern showing two lateral lobes enfolding the dorsal (middle) and ventral (innermost) petal lobes (Figure 1(Cd)). The floral zygomorphy along the dorsal-ventral axis therefore initiates at FB1-7 (the dashed line in Figure 1(Cd)). The development of the aborted dorsal staminode probably also follows a dorsal-ventral asymmetry pattern as early as at FB1-8 ( Figure S1A and Table S1).
In the later floral stages, FB3-FB16, the flower buds gradually assumed the differences in shape between the dorsal and ventral corolla parts ( Figure 1D). The development of flower dorsal-ventral asymmetry was characterized by a differential degree of curvature between the dorsal and ventral parts of the corolla tube (yellow dashed curved lines). Additionally, this asymmetry was accompanied by a convex, gibbous outgrowth (nectar chamber) at the base of the dorsal corolla tube (red curve), and a ventricose chamber (blue curve) in the ventral part of the corolla tube ( Figure 1D). The difference in the tube lengths between the dorsal and ventral parts of the corolla became evident starting at developmental stage FB10 ( Figure 1E). The dorsal tube length gradually became longer than the ventral part of the tube during FB10 and FB12. However, the length difference between the dorsal tube and ventral tube narrowed again at FB16. Moreover, the number of tube cells in the dorsal and ventral parts remained almost equal across all flower developmental stages ( Figure 1F). Thus, we can conclude that the dorsal-ventral corolla size difference was mainly due to the larger degree of cell expansion in length in the dorsal part compared to the ventral part of the corolla. We further examined surface cellular micromorphology of the petal epidermal cells with a focus on cell size and shape differences using cryo-SEM at FB5 and FB8 stages (Figure 2A,B). The dorsal part of the corolla had larger, flattened, and rougher epidermal cells ( Figure 2C,D vs. Figure 2K,L), while the ventral part had relatively smaller, more protruding conical cells and a smooth epidermis, particularly in the proximal and middle regions ( Figure 2G,H vs. Figure 2O,P). Moreover, there are more granule-like epicuticular wax crystals deposited on the epidermal surface of the dorsal part of the corolla tube than on the ventral part ( Figure 2K-R and Figure S1B). Thus, the development of petal epidermal micromorphology in zygomorphic S. speciosa 'ES' corolla petals followed a dorsal-ventral asymmetry pattern. cryo-SEM at FB5 and FB8 stages (Figure 2A,B). The dorsal part of the corolla had larger, flattened, and rougher epidermal cells ( Figure 2C,D vs. Figure 2K,L), while the ventral part had relatively smaller, more protruding conical cells and a smooth epidermis, particularly in the proximal and middle regions ( Figure 2G,H vs. Figure 2O,P). Moreover, there are more granule-like epicuticular wax crystals deposited on the epidermal surface of the dorsal part of the corolla tube than on the ventral part (Figures 2K-R and S1B). Thus, the development of petal epidermal micromorphology in zygomorphic S. speciosa 'ES' corolla petals followed a dorsal-ventral asymmetry pattern.

Flower Coloration in Zygomorphic S. speciosa Displays a Distinct Dorsal-Ventral Pigmentation Pattern
The development of flower dorsal-ventral asymmetry in S. speciosa 'ES' is not only apparent in the differences in corolla shape and size, but is also evident in pigmentation pattern. The corolla begins to show its purple gradient coloration along the dorsal-ventral and proximal-distal axes at stages FB8, especially on the dorsal side, while at FB5 and earlier stages, the entire corolla is green ( Figure 1D). The color becomes more intense and spreads from the dorsal tube to the dorsal lobes beginning in stage FB10 ( Figure 1G). The color is a sign that purple anthocyanin has accumulated. This can be observed by the accumulation of purple pigments in corolla epidermal cells ( Figure 1H). On the dorsal tube of stage FB8 flowers, a number of epidermal cells, particularly those at the tip of the trichomes, show the purple color ( Figure 1(Hb)). At anthesis (stage FB16), almost all epidermal and trichome cells on the dorsal surface of the corolla highly accumulate purple pigments ( Figure S1C). On the inner (adaxial) surface of the ventral tube, certain epidermal cells may also produce purple anthocyanin starting at FB5 (Figure 1(Hc)). These purple-colored epidermal cells later accumulate in clusters from stage FB8 to anthesis at FB16 (Figure 1(Hd) and Figure S1C), thus exhibiting a spotting pattern on the inner surface of the ventral tube. Therefore, the display of the ventral petal spotting pattern also follows a position-dependent pattern along the dorsal-ventral axis of the flower in S. speciosa 'ES'.

Transcriptome Analysis of DEGs between the Dorsal and Ventral Petals
RNA-Seq data between the dorsal and ventral petals obtained from FB5 were chosen because dorsal-ventral petal asymmetry was established and high SsCYC expression in the dorsal part was observed ( Figure S2A,B). A total of 262,599,569 clean reads from four libraries (ZD-1, ZD-2, ZV-1, and ZV-2) were obtained (Table S2) and were individually mapped to the predicted transcripts from draft genome sequences of S. speciosa 'Avenida Niemeyer' ('AN'), which is closely related to S. speciosa 'ES'. Of the clean reads, 84.5-85.3% aligned (Table S2), indicating that the draft genome sequence was an acceptable reference for our transcriptomic analysis. A total of 630 dorsal-ventral differentially expressed genes (DV-DEGs) were obtained from the transcriptomic comparison between ZDs and ZVs. The heatmap shows the relative expression levels of the 294 dorsally high and 336 ventrally high DEGs ( Figure S3).

Functional Annotation and GO Enrichment Analysis of the Dorsal-Ventral DEGs
The DV-DEGs were then functionally annotated against the GeneBank NR, GO, and KEGG public databases (Table S3). GO functional classification assigned the DV-DEGs to 796 GO terms in the three main GO categories and 39 sub-categories. GO terms in the 'biological process' category, including 'biological regulation', 'localization', and 'cellular component organization or biogenesis', were highly abundant ( Figure S4). We next performed GO enrichment analysis to further investigate the over-represented GO terms of the DV-DEGs. DEGs annotated as GO terms related to 'regulation of size', 'cell wall modification', 'hormone-related', and 'growth and development' were the most significantly enriched (Figure 3, blue, green, and cyan blocks, respectively). Hormone-related GO terms and auxin-related terms such as 'auxin mediated signaling pathway', 'response to auxin stimulus', and 'auxin transport' were significantly enriched ( Figure 3, pink block). In addition, the DEGs associated with the regulation of secondary metabolic processes such as 'lipid metabolic', 'lignin metabolic processes', and 'flavonoid biosynthesis' were also significant. There were also some enriched GO terms correlated with the response to light (UV and blue light) and stress (biotic and abiotic). These results imply that auxin signaling, cell wall expansion, and developmental processes may be associated with the differential growth of dorsal-ventral petal size and shape, while secondary metabolic pathways probably contribute to asymmetric pigmentation patterns in the dorsal and ventral petals.

Auxin Signaling Transduction Genes Identified as Major DV-DEGs
To further understand the biological functions of the DV-DEGs, pathway-based analysis was conducted. A total of 301 DV-DEGs (47.8%) were successfully annotated to 177 KEGG pathways. Notably, 18 DV-DEGs were assigned to hormone signaling transduction pathways (Table S4). A total of 9 of these 18 DV-DEGs could be mapped to auxin signaling . Node size is proportional to the number of genes, and the color scale represents the significance levels: white, no significant difference; yellow, p-value < 5 × 10 −2 ; orange, p-value < 5 × 10 −7 . Color shadings refer to the GO terms that relate to six broad biological functions. Arrows indicate connections between GO annotations.
The qRT-PCR results of selected DEGs were quite consistent with the gene expression profiles of RNA-Seq (Supplementary Table S8). Among the 27 DV-DEGs we validated, 26 of them shared the same DV expression patterns between RNA-Seq and qRT-PCR, supporting the accuracy of RNA-Seq. Since only two biological replicates were used in RNA-Seq, genes that did not show a large differential expression between dorsal and ventral petals may not have been identified. Most of these signaling transduction pathway DEGs including SsAUX1, SsAUX/IAAs, SsARF2, and SsGH3 were found to be dorsally highly expressed, whereas only SsARF16 had a high ventral expression, especially at stage FB5 ( Figure 5A,D). Although SsARF1 and SsARF3 were not identified as DV-DEGs in the comparative transcriptomes of FB5, SsARF1 and SsARF3 also had a high dorsal transcript at the FB3 stage and FB8 stage, respectively ( Figure 5D). Figure 3. The Gene Ontology (GO) enrichment analysis for the significantly enriched GO terms belonging to the 'biological process' GO category for the 630 DV-DEGs in S. speciosa 'ES'. Node size is proportional to the number of genes, and the color scale represents the significance levels: white, no significant difference; yellow, p-value < 5 × 10 −2 ; orange, p-value < 5 × 10 −7 . Color shadings refer to the GO terms that relate to six broad biological functions. Arrows indicate connections between GO annotations.

Auxin Signaling Transduction Genes Identified as Major DV-DEGs
To further understand the biological functions of the DV-DEGs, pathway-based analysis was conducted. A total of 301 DV-DEGs (47.8%) were successfully annotated to 177 KEGG pathways. Notably, 18 DV-DEGs were assigned to hormone signaling transduction pathways (Table S4). A total of 9 of these 18 DV-DEGs could be mapped to auxin signaling Based on our reconstructed phylogeny of the ARF family, both SsARF1 and SsARF2 in clade II were close to SsARF3, while SsARF3 and SsARF16 were grouped to clade IIIc ( Figure S6A). MicroRNA-target tool analysis (see Supplementary Methods) further predicted that both SsARF2 and SsARF16 were putative miRNA targets. SsARF16 formed a near complementary sequence pair with miR160 with a 1 nt mismatch at the 3 end, while SsARF3 base-paired with TRANS-ACTING SIRNA3 (TAS3)-derived trans-acting short-interfering RNA (ta-siRNA) ( Figure S6B). This implied that the small RNA regulatory mechanism on auxin response may also be involved in dorsal-ventral petal asymmetric growth.
Apart from auxin signaling transduction, auxin conjugation and deconjugation genes were also among the dorsally highly expressed DEGs. SsILL6-1 and SsILL6-2 were classified in the amidohydrolase ILR1-like 6 (ILL6) group for deconjugating auxin amides, whereas SsGH3 was in the group II GH3 family that catalyzes the synthesis of IAA amide conjugates in the regulation of free cellular auxin levels. The homeostasis of auxin storage and release was thus also associated with dorsal-ventral petal asymmetric growth.
While many auxin pathway genes were found to be DV-DEGs, very few other hormone pathways genes were identified in DV-DEGs ( Figure S7). These include SsGID1 (GA receptor), SsABA2 (ABA biosynthesis), and SsOPR3 (JA biosynthesis). They also showed dorsally high transcriptions ( Figure S8).

Cell Wall Growth Genes Identified in DV-DEGs
The cell wall modification pathway is a putative response to auxin for regulating cell growth. We also identified several major cell wall modification pathway genes among the DV-DEGs. These included SsEXPs (EXPANSIN), SsXTH (xyloglucan endo transglucosylase/hydrolases), SsPEX1 (pollen-specific LRR/EXTENSIN1), and SsPME (pectin methylesterase). They encode proteins homologous to cell wall modulators such as structural proteins and proteins that modify interactions related to the regulation of auxininduced cell wall loosening during cell growth [23]. Specifically, there were six SsEXP family members identified in the DV-DEGs (Table S5). The SsEXP genes were in clade I (EXPA1 and EXPA2), clade IV (EXPA4, EXPA5, and EXPA6), and clade III (EXPA3) of the EXPA subfamily based on our reconstructed phylogeny ( Figure S9). Int. J. Mol. Sci. 2022, 23, x FOR PEER REVIEW 9 of 21  qRT-PCR validation of these DV-DEGs revealed that most of these cell wall modification genes were dorsally highly expressed. It was found that two SsEXP family genes, SsEXPA1 and SsEXPA2, and two cell wall modification genes, SsPME and SsPEX1, had a high dorsal expression pattern, while two SsEXP genes, SsEXPA4 and SsEXPA5, were ventrally high ( Figure 5B).

Transcription Factors and CYC-RAD Modules Identified in DV-DEGs
To identify the possible master regulators responsible for asymmetric dorsal-ventral petal growth in S. speciose 'ES', we performed transcription factor (TF) predictions using the iTAK v1.6 and the Plant Transcription Factor Database v4.0 (see Supplementary Methods) of these DV-DEGs. A total of 42 putative differentially expressed TF genes (DETFs) were predicted. This included 20 and 22 DETF genes that were expressed at high levels in the dorsal and ventral petals, respectively. As expected, homologs of floral symmetry genes CYC (SsCYC) and RAD (SsRAD1 and SsRAD2) were all identified in the DETF, with higher expression in the dorsal petals (Table S6; Figure S2A,C and Figure 5C). Phylogenetic analysis revealed that SsRAD2 was clustered in the RAD2 clade, homologous to RAD from A. majus, while SsRAD1 was in the distantly related RAD1 clade ( Figure S10A). SsRAD2 has two TCP binding sites in the 5 regulatory region ( Figure S10B), which suggests that  Table S6), all of which are known to function in regulating cell growth for plant development.

Transcription Factors of Anthocyanin Biosynthesis Identified in DV-DEGs
We also identified DV-DEGs that encode transcription factors that regulate the flavonoid biosynthesis pathway (Figure 3). In particular, an MYB transcription factor containing the MYBL2-like domain, SsMYBL2, was found to be mainly expressed in dorsal petals at the FB5 and FB11 stages ( Figure 5C) when the dorsal corolla turned from green to purple. MYBL2 has been reported to be the major regulator for anthocyanin synthesis/flavonoid biosynthesis [24]. In addition, several structural genes in the flavonoid and anthocyanin biosynthetic pathways, such as SsFLS (flavonol synthase), SsF3 H (flavonoid 3 -hydroxylase), and SsGT (UDP-glucose: flavonoid 3-O-glucosyltransferase), were also identified in the DV-DEGs ( Figure S11).

SsCYC Localized to Nucleus and Regulated Petal Size in Nicotiana
As a major regulator for floral symmetry, we asked whether the SsCYC protein acts as a TF to enter the nuclear region. Using the petal protoplast transient expression system, we found that more than 95% of the fluorescence signal was clearly localized to the nucleus in the 35S:SsCYC-GFP transfected cells ( Figure 6, upper row). Notably, 4.87% of the fusion protein was also located in the nucleus and cytoplasm ( Figure 6, middle row), similar to a previous report in strawberries which found that the TCP protein can also localize to the cytoplasm [25].
Given the fact that most of the DV-DEGs identified in SsCYC that were highly expressed in the floral bud stage are those genes relating to auxin signaling and cell growth, we hypothesized that SsCYC and auxin play a combined role in determining dorsal-ventral petal size and/or shape differences. To investigate SsCYC's function in petal growth, we ectopically expressed it in transgenic plants of N. benthamiana. Genomic insertions were validated by PCR, and elevated mRNA levels of SsCYC were confirmed in transgenic lines ( Figure S12). 35S:SsCYC transgenic plants generally had various degrees of dwarf phenotypes (OX-0~OX-3) ( Figure 7A). Noticeably, all petals of the transgenic plants were greatly reduced in size with narrower lobes towards their distal ends ( Figure 7B). The measurement of petal epidermal cells in transgenic flowers revealed that cell areas were smaller (969.14 ± 49.19 µm 2 and 721.03 ± 37.75 µm 2 ) in the middle and distal regions compared to wildtype flowers (1483.13 ± 78.60 µm 2 and 1524.65 ± 59.49 µm 2 ), but not in the proximal region ( Figure 7C,D). This suggests that SsCYC may function to repress cell expansion to affect petal size and shape.

Transient Over-Expression of the SsCYC in Petal Protoplasts
Since we identified DV-DEGs involved in the auxin pathway at the key bud stage when SsCYC was also highly expressed, we also tested whether transient over-expression of SsCYC can regulate these DV-DEGs in petal protoplasts. The expression levels of SsRADs, auxin homeostasis, auxin signaling, and cell growth genes were examined when 35S:SsCYC was transiently over-expressed in petal protoplasts of S. speciose 'ES'. Including SsRAD2, the expression of one auxin signaling gene (SsAUX/IAA1), one auxin homeostasis gene (SsILL6-2), and two cell growth genes (SsEXPA1 and SsEXPA4) was upregulated with SsCYC over-expression ( Figure S13A-F). However, the expression of other auxin pathway genes such as SsILL6-1(auxin homeostasis), SsARF2 (auxin response factor), and one cell growth gene (SsEXPA2) was not induced by SsCYC over-expression ( Figure S13G-I). Therefore, it remains inconclusive as to whether SsCYC regulates the auxin pathway and cell growth genes.

Endogenous Auxin Level Exhibited a Dorsally High Pattern
To estimate the difference in endogenous auxin levels between the dorsal and ventral petals during floral zygomorphy development, IAA was extracted from FB3, FB5, and FB8 ( Figure 8). IAA levels showed no difference between dorsal and ventral petals at FB3 when dorsal-ventral petal asymmetry was just about to initiate. However, at FB5 when zygomorphy was establishing ( Figure 1D), the IAA level was significantly higher in dorsal petals than in ventral petals (twice as much as that in the ventral petals). This correlated with the fact that SsCYC also showed the highest dorsal expression at stage FB5 ( Figure 5C). The difference in IAA levels between dorsal and ventral petals was no longer obvious at FB8 when flower zygomorphy was fully developed. Similarly, the DV differential expression pattern of SsCYC had almost disappeared by stage FB8. This indicates that the dorsalventral auxin level difference correlated with asymmetric dorsal-ventral petal growth, and also the DV differential expression of SsCYC during floral zygomorphy development.

Transient Over-Expression of the SsCYC in Petal Protoplasts
Since we identified DV-DEGs involved in the auxin pathway at the key bud stage when SsCYC was also highly expressed, we also tested whether transient over-expression of SsCYC can regulate these DV-DEGs in petal protoplasts. The expression levels of SsRADs, auxin homeostasis, auxin signaling, and cell growth genes were examined when 35S:SsCYC was transiently over-expressed in petal protoplasts of S. speciose 'ES'. Including SsRAD2, the expression of one auxin signaling gene (SsAUX/IAA1), one auxin homeostasis gene (SsILL6-2), and two cell growth genes (SsEXPA1 and SsEXPA4) was upregulated with SsCYC over-expression ( Figure S13A-F). However, the expression of other auxin pathway genes such as SsILL6-1(auxin homeostasis), SsARF2 (auxin response factor), and one cell growth gene (SsEXPA2) was not induced by SsCYC over-expression ( Figure  S13G-I). Therefore, it remains inconclusive as to whether SsCYC regulates the auxin pathway and cell growth genes.

Endogenous Auxin Level Exhibited a Dorsally High Pattern
To estimate the difference in endogenous auxin levels between the dorsal and ventral petals during floral zygomorphy development, IAA was extracted from FB3, FB5, and FB8 ( Figure 8). IAA levels showed no difference between dorsal and ventral petals at FB3 when dorsal-ventral petal asymmetry was just about to initiate. However, at FB5 when zygomorphy was establishing ( Figure 1D), the IAA level was significantly higher in dorsal petals than in ventral petals (twice as much as that in the ventral petals). This correlated with the fact that SsCYC also showed the highest dorsal expression at stage FB5 ( Figure  5C). The difference in IAA levels between dorsal and ventral petals was no longer obvious at FB8 when flower zygomorphy was fully developed. Similarly, the DV differential expression pattern of SsCYC had almost disappeared by stage FB8. This indicates that the dorsal-ventral auxin level difference correlated with asymmetric dorsal-ventral petal growth, and also the DV differential expression of SsCYC during floral zygomorphy development. The mean values ± SD are from three biological repeats. Comparisons between the two groups were performed using a Student's t test, *** p < 0.001.

Expression of the Auxin Pathway and Cell Growth Genes Correlate to DV Petal Asymmetry
Our transcriptomic analysis identified major auxin signaling pathway genes which are dorsally highly expressed, and this was associated with a longer dorsal tube length The mean values ± SD are from three biological repeats. Comparisons between the two groups were performed using a Student's t test, *** p < 0.001.

Expression of the Auxin Pathway and Cell Growth Genes Correlate to DV Petal Asymmetry
Our transcriptomic analysis identified major auxin signaling pathway genes which are dorsally highly expressed, and this was associated with a longer dorsal tube length ( Figure 1D,E and Figure 4). Larger petal size has been shown to be enhanced by auxin via the triggering of auxin signaling pathway genes, both in Arabidopsis and orchids [26]. Our qRT-PCR validation revealed that major auxin homeostasis (two SsILL6s and one SsGH3), transport (SsPIN and SsAUX1), signaling (SsAUX/IAAs and SsARF2), and early auxin response (SsSAUR) genes all had high dorsal expression levels in bud stage FB3 (zygomorphy initiating), but their expression increased even further at FB5 (zygomorphy establishing) ( Figure 5A). However, the difference in their expression in ventral petals became less evident towards FB11 (when zygomorphy was already fully developed). Similarly, the expression of SsCYC also exhibited high dorsal expression from FB3 to FB5, but this trend was no longer evident at FB11 ( Figure 5C). The expression of auxin signaling pathway genes thus correlates with the expression of SsCYC, particularly during the zygomorphy establishing stage (FB5) ( Figure 1D). These results suggest that auxin pathway regulation, together with SsCYC, is involved in DV petal asymmetry.
EXP family genes have been reported to be regulated by auxin signaling to control cell growth (cell wall modification) by promoting cell expansion [23]. Our DV-DEG and qRT-PCR validation also found that two of these cell growth genes (SsEXP1 and SsEXP2) had high dorsal expression ( Figure 5B). Although this suggests that the activation of auxin signaling in dorsal petals may mediate larger cell growth, this is arguable because there were also two SsEXPs genes (SsEXP4 and SsEXP5) that exhibited high ventral expression ( Figure 5B). One possible explanation could be that these SsEXPs genes each regulate the local outgrowth of corolla tubes such as the dorsal gibbous structure and ventral ventricose chamber. Previously, a phenotype-genotype association from geometric morphometric analysis suggested that the functional SsCYC allele correlates with the dorsal gibbous structure and ventral ventricose chamber [21]. These SsEXPs may act similarly to SsCYC in controlling these local outgrowths of the corolla tube. It would be a good idea in the future to check whether the expression of these SsEXPs is concentrated in these petal outgrowth regions.

Auxin Dorsal-Ventral Asymmetric Distribution Correlate to DV Petal Asymmetry
Local micro-applications of auxin to floral primordia have been shown to convert typical zygomorphic flowers to ventralized actinomorphic peloric flowers in A. majus, probably because of the disturbance of normal endogenous auxin distribution, resulting in transient local auxin concentrations that may vary in floral primordia [27]. In Asteraceae species Matricaria inodora, the spraying of auxin on the capitulum inflorescence can convert the actinomorphic disc florets into zygomorphic ray florets by enhancing the expression of the CYC homolog (RAY2) [28]. Recently, Liu et al. [29] identified one putative ARE (auxin responsive element) site in a CpCYC promoter of Chirita pumila. The EMSA result showed that this ARE site can be bound by floral nucleoproteins, implying that auxin may regulate CpCYC expression. These findings suggest that auxin may regulate dorsal-ventral petal asymmetry through the activation of CYC.
Our finding of a higher endogenous auxin level in dorsal petals thus provides new insight that the asymmetric distribution of auxin may regulate asymmetric dorsal-ventral petal growth. At bud stage FB3, when petal zygomorphy is just about to develop, there is no difference in endogenous auxin levels between dorsal and ventral petals (Figure 8). However, at bud stage FB5, the auxin in dorsal petals increased to a significantly higher level, which coincides with the key stage of establishing petal zygomorphy ( Figures 1D and 8). This also correlates with the highest SsCYC expression level at bud stage FB5 ( Figure 5C). We further tested the hypothesis that SsCYC regulates auxin pathway genes. However, when we transiently over-expressed SsCYC in petal protoplasts, not all selected auxin pathway and cell growth genes were activated by SsCYC ( Figure S13). This suggests that over-expressing SsCYC was not able to activate the expression of most of the auxin pathway genes. Thus, it is less likely that SsCYC regulates the auxin pathway. Rather, our results, together with the prediction from previous studies, suggest that auxin may regulate dorsal-ventral petal asymmetry, perhaps by promoting SsCYC expression.
In the future, we will treat auxin or auxin inhibitor on flower buds of S. speciosa 'ES' and subsequently examine whether SsCYC can be upregulated or downregulated. Moreover, we will examine whether the auxin level is increased in SsCYC-over-expressing Nicotiana flowers. These studies will help to clarify whether auxin and SsCYC work dependently or independently in determining zygomorphy.

SsCYC May Regulate Petal Size in a Context-Dependent Manner
Our results show that ectopically expressing SsCYC in Nicotiana resulted in a reduction in petal size by repressing cell expansion (smaller cell area, see Figure 7B,C). This is similar to the result from another S. speciosa cultivar ('Pink Flower', 'PF') that shows that the overexpression of SsCYC in A. thaliana also repressed petal cell expansion [19]. However, this is different from the observation in S. speciosa 'ES' itself, where the SsCYC-expressed dorsal petals have longer tube lengths due to the expanded epidermal cell size in the proximal part of the dorsal tube ( Figures 1E and 2C,D,K,L). This can be explained by the fact that the regulation of CYC has been suggested as being context-dependent [30]. In different species and tissues, CYC might function in protein complexes by incorporating different co-factors to activate or repress genes, thus controlling cell expansion or cell proliferation. Indeed, in Antirrhinum majus, the expression of AmCYC promotes dorsal petal growth by increasing levels of cell expansion [8], while in Iberis amara, IaCYC reduces petal cell growth by repressing cell proliferation [16]. CYC may also recruit different downstream targets in different lineages. For example, CYC has been documented to activate RAD for A. majus, but RAD was not found downstream of TCP1 (CYC homolog) for A. thaliana [17,31]. Therefore, it is possible that SsCYC still plays a role in regulating dorsal-ventral asymmetry in S. speciosa 'ES', since a strong association between SsCYC expression pattern ( Figure 5D) and floral asymmetry development ( Figure 1E) was observed, and more importantly, some of the DV-DEGs were induced by SsCYC over-expression (Figure S13A-F). However, SsCYC might have activated different downstream targets between S. speciosa 'ES' and Nicotiana because of context-dependent regulation.

SsARF Genes May Function to Modulate Petal Size and Shape
Auxin response factors (ARFs) are TFs that transduce the auxin signal, and each ARF has its own distinct function and expression pattern in organ development. In the Phalaenopsis orchid, several PeARFs were found to be differently expressed between flowers of zygomorphic wildtype and peloric actinomorphy during the stage of floral symmetry establishment [32]. Among the DV-DEGs of the auxin pathway, we identified SsARF2, which showed high dorsal expression at FB5 (zygomorphy establishing) to FB8 ( Figure 5D). As a higher auxin level was detected in the dorsal petals at FB5 (Figure 8), the upregulation of SsARFs may act in response to auxin for establishing zygomorphy.
ARF3 has been well-characterized by its roles in organ polarity as a leaf abaxialdetermining gene [33]. TAS3 ta-siRNA-mediated regulation of ARF3/ARF4 plays a critical role in establishing leaf adaxial-abaxial polarity [34]. Preliminary data shows that ta-siARF is expressed in ventral petals, and thus perhaps has a role in eliminating ARF3 expression in the ventral side. It remains to be explored whether the development of dorsal-ventral petal asymmetry shares a similar mechanism of leaf adaxial-abaxial polarity in S. speciosa 'ES'.
The role of ARF3 in flower development has been reported to function in retarding petal growth and disrupting petal fusion into corolla [35,36]. Over-expression of tomato ARF3 siRNA-resistant form (35S: mSlARF3) in tobacco, Nicotiana benthamiana, in wiry tomato (a siRNA mutant), and in Mimulus lewisii caused the whole flower to become smaller in size, and the petals to become narrowly unfused and separated [35,36]. Thus, ARF3 in different angiosperm lineages has common activities in regulating petal growth.
Our results also reveal that SsARF16 is a ventrally expressed DV-DEG ( Figure 5D), which is phylogenetically grouped with A. thaliana ARF10, ARF16, and ARF17 ( Figure S6A). These three ARF genes have been reported to be regulated by miR160 to control petal size and shape in Arabidopsis [37]. When the miR160-resistant ARF genes are driven by their own promoters (e.g., Pro ARF16 :mARF16), the flowers of the transgenic plants have defects with inwards-curled petals, suggesting that ARF16 negatively regulates cell division and cell expansion.

SsMYBL2, a Potential Regulator of Anthocyanin Biosynthesis, May Contribute to Asymmetric Dorsal-Ventral Petal Pigmentation Pattern
The flowers of S. speciosa 'ES' display a distinct dorsal-ventral petal pigmentation pattern after zygomorphy is established from FB10 to FB16 ( Figure 1D). The violet to purple corolla coloration has been reported to be mainly due to the presence of anthocyanins in Sinningia species [6,38]. This relates to our findings that SsMYBL2, a potential regulator of the anthocyanin biosynthesis pathway gene, was expressed more in dorsal petals in FB3 and FB5, but that the highest DV expression difference was at FB11 ( Figure 5C). FB11 is a stage at which zygomorphy is fully developed. Thus, the development of asymmetric DV petal pigmentation is perhaps a downstream regulator of flower zygomorphy.
SsMYBL2 is a MYB-bHLH-WD40 (MBW) transcriptional activation complex that is highly conserved across angiosperms and includes members of the R2R3-MYB or R3-MYB, bHLH, and WDR transcription factor gene families [39][40][41]. In Arabidopsis, R3and R2R3-MYBs are reported to interact with AtTCP3 to form various combinations of TCP-MYB heterodimers that activate anthocyanin accumulation [42]. A recent finding in Torenia fournieri suggests that dorsally expressed TfCYC2 can directly regulate TfMYB1, which promotes anthocyanin pigmentation in the epidermal cells of petals, leading to the asymmetric dorsal-ventral petal pigmentation pattern [43]. In the future, it would be interesting to know whether a similar regulatory relationship between SsCYC and SsMYBL2 is also established in S. speciosa 'ES'.
The major finding of this paper suggests that auxin regulation is important for establishing dorsal-ventral petal asymmetry in S. speciosa 'ES', based on the results of transcriptomic analysis, the validation of DV-DEG, and the high endogenous auxin level in dorsal petals. In dorsal petals, due to the higher endogenous auxin levels detected, auxin-signaling and homeostasis genes may be upregulated in response to elevated auxin (brown box in Figure 9). This may subsequently modulate cell wall loosening genes to regulate cell size (green boxes in Figure 9). As over-expressing SsCYC in Nicotiana also affects cell expansion, SsCYC may act directly or be mediated by auxin to regulate petal cell size (dashed lines in Figure 9). We also found a major regulator of anthocyanin biosynthesis, SsMYBL2, to be highly expressed in dorsal petals, perhaps functioning to generate dorsal-ventral petal purple pigmentation differences.

Plant Materials
Sinningia speciosa 'ES' was originally from M. Peixoto of Mogi das Cruzes, SP, Brazil [44]. Sinningia speciosa 'ES' is a wild accession with a zygomorphic flower. A total of 120 individual plants of S. speciosa 'ES' F2 population were grown in a walk-in greenhouse at the National Taiwan University (Taipei City, Taiwan) under an LED light intensity of 200 μmol·m -2 ·s -1 , and a 16 h light/8 h dark cycle with 70-80% relative humidity at 22-25 °C.

Morphological Observation
We divided floral buds into FB1 to FB16 stages according to the length of the dorsal corolla tube (Figure 1A,B). To examine petal and stamen primordia initiation, we further

Plant Materials
Sinningia speciosa 'ES' was originally from M. Peixoto of Mogi das Cruzes, SP, Brazil [44]. Sinningia speciosa 'ES' is a wild accession with a zygomorphic flower. A total of 120 individual plants of S. speciosa 'ES' F2 population were grown in a walk-in greenhouse at the National Taiwan University (Taipei City, Taiwan) under an LED light intensity of 200 µmol·m -2 ·s -1 , and a 16 h light/8 h dark cycle with 70-80% relative humidity at 22-25 • C.

Morphological Observation
We divided floral buds into FB1 to FB16 stages according to the length of the dorsal corolla tube (Figure 1A,B). To examine petal and stamen primordia initiation, we further divided stage 1 floral buds (when the length of the dorsal tube was smaller than 2 mm) into seven sub-stages, namely FB1-1-FB1-7, using scanning electron microscopy (SEM). FB5 (when the dorsal tube length = 8-10 mm) is the key stage during which dorsal-ventral petal asymmetry is established-that is, the length and shape of the dorsal tube starts to differentiate from those of the ventral tube. FB8 is the stage at which purple color starts to accumulate in dorsal petals (tubes and lobes) but not in ventral petals. Flower buds at stages FB5 and FB8 (when the dorsal tube length = 15-18 mm) were therefore harvested to observe the epidermal cell morphology of dorsal and ventral petals by means of cryo-scanning electron microscopy (Cryo-SEM). Additional information on the SEM and Cryo-SEM methods are described in the Supplementary Materials. For tube size and cell number determination, we measured tube length following the outline curvature of the corolla tube on fresh floral samples, and counted epidermal cell numbers along the main vein of the corolla tubes after the tissue was cleared.

RNA Extraction, cDNA Library Construction, Sequencing, and Read Mapping
To identify GRN associated with floral symmetry, we separated RNA-Seq between the dorsal and ventral petals (comprising both the tube and lobe) at FB5 when dorsalventral petal asymmetry was established and SsCYC was highly dorsally expressed. Two dorsal petals and one ventral petal, excluding the two lateral petals of each flower, were dissected and collected. For each biological replicate, petals were sampled from 10 to 20 buds. The samples were harvested immediately and frozen in liquid nitrogen and stored at −80 • C before they were used for RNA extraction. Total RNA was extracted using Trizol Reagent according to the manufacturer's instructions (Invitrogen, Carlsbad, CA). RNA integrity was checked using an Agilent 2100 Bioanalyzer with an RNA 6000 Chip and 2100 Expert software (Agilent Technologies). rRNA depleted mRNAs were prepared using the Ribo-Zero Gold rRNA Removal Kit. All RNA sequencing libraries were constructed using the NEBNext Ultra Directional RNA Library Prep Kit for Illumina according to the manufacturer's instructions. Four cDNA libraries (with two biological replicates) of the dorsal (ZD-1, ZD-2) and ventral (ZV-1, ZV-2) petals were constructed for RNA-Seq. The double-strand cDNA libraries were sequenced on an Illumina Hi-Seq 2000 sequencing platform. More than 60 million paired-end reads of 150 nucleotides in length were obtained for each sample. The raw RNA-Seq data are available at the National Center for Biotechnology Information Short Reads Archive (NCBI SRA) accession number SRR16636567-SRR16636570 at the following URL: https://www.ncbi.nlm.nih.gov/Traces/ study/?acc=PRJNA776263, accessed on 12 November 2021). Adaptor sequences and lowquality reads were removed from the raw reads using Trimmomatic version 0.35 [45]. Clean reads were then mapped to the predicted transcripts (31,581 predicted transcripts from 8078 assembled scaffolds) and the draft genome assembly of S. speciosa 'Avenida Niemeyer' from Rio de Janeiro city, which is closely related to S. speciosa 'ES', kindly provided by Dr. Aureliano Bombarely at Virginia Tech University (Blacksburg, VA, USA) (Bombarely and Zaitlin, unpublished data), using Bowtie2 version 2.0.6 [46] and BLAT version 34 [47], respectively. Genomic regions not annotated with the aforementioned transcripts with a maximum mapping depth greater than 10 and a length longer than 300 bps were identified as novel gene loci and were included in the following differential gene discovery (Table S9).

Differential Gene Expression Analysis and Functional Annotation
Transcript abundance was computed based on read counts normalized by the trimmed mean of M-values (TMM) method using the Bioconductor edgeR package [48]. Differentially expressed genes (DEGs) between dorsal and ventral petal samples (ZDs vs. ZVs) were identified if the t-test p-value was less than 0.05 and |log2 fold-change| > 1 (see Table  S10 for the comparison table). The dorsal-ventral DEGs (DV-DEGs) were annotated on the basis of their sequence similarities with known protein annotations in public databases. Heatmaps of DEG expression were generated with the ClustVis Web tool [49]. Functional annotations were performed using Gene Ontology (GO) classification, KEGG pathways, and transcription factor prediction. The details of these experiments are described in the Supplementary Material.  Figure 1D). RNA was treated with RQ1 RNase-Free DNase (Promega, Madison, WI, USA) to remove contaminating genomic DNA. First-strand cDNA was synthesized using MMLV reverse transcriptase (Invitrogen, Carlsbad, CA). Gene-specific primers for the DEGs (Supplementary Table S7) were designed using Primer3Plus [50]. qRT-PCR assays were performed using the CFX Connect Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA) with the KAPA SYBR ® FAST qPCR Master Mix (2X) Kit (Sigma-Aldrich, St. Louis, MO, USA). Two biological replicates were performed in the qRT-PCR, except for SsARFs, for which three biological repeats were performed. The data were analyzed using CFX Maestro Software (Bio-Rad, Hercules, CA). The expression level of S. speciosa 18S ribosomal RNA (Ss18S rRNA) was used as the internal control [51,52]. Normalized mRNA expression levels were calculated as 2 −∆Ct .

S Phylogenetic Analysis of the ARF, EXP, and RAD Families
The dorsiventral DEGs annotated as encoding ARFs, EXPs, and RADs were newly identified from S. speciosa 'ES' using primers that were designed based on the sequences of cultivar of 'AN'. The translated coding sequences were aligned by MUSCLE in MEGA5.2 with the default settings [53]. Phylogenetic relationships were reconstructed using the maximum likelihood (ML) criteria [54] (Method S1). The best-fit substitution model of each dataset was implemented in PhyML 3.0 [55]. The best-fit substitution model of each dataset was evaluated by smart model selection (SMS), which was implemented in PhyML 3.0. For inferring ARF relationships, 22 AtARFs from A. thaliana and 4 SsARFs from S. speciosa 'ES' were included, and the best-fit model was JTT + G + I + F. For the EXPs dataset in S. speciosa 'ES', 34 AtEXPs (A. thaliana), 11 LeEXPs (Solanum lycopersicum), 5 NtEXPs (Nicotiana tabacum), and 6 SsEXPs were included, and the best-fit model was the LG model. Details of EXP sequence information are listed in Supplementary Table S8. For the RAD family, 11 SsRADs and 60 RADs homologs from monocots and dicots were included, and the best-fit model was GTR + G + I model.

Subcellular Localization, Functional Characterization, and Transient Activation of SsCYC
As the SsCYC protein is a transcription factor, we transiently transfected a DNA construct containing 35S:SsCYC-GFP into a petal protoplast to confirm its subcellular localization into the nuclear region. 35S:GFP construct was also transfected as a control. The details of these procedures are described in the Supplementary Materials. To further investigate the biological function of SsCYC, we also created 35S:SsCYC transgenic plants of Nicotiana benthamiana to characterize the effect of SsCYC on flower development, particularly petal growth. Confirmation of T-DNA insertion of SsCYC by Ti plasmid and mRNA expressions in the T 1 transgenic plants was performed using PCR and RT-PCR. The details of the transfection procedures are provided in the Supplementary Materials (section "petal protoplast transfection system"). To explore whether certain DV-DEGs can be regulated by SsCYC, we transfected the 35S:SsCYC-35S:GFP effector construct into petal protoplasts, then quantified whether the transcriptional levels of candidate genes increased or decreased (putatively induced or repressed) using qRT-PCR.

Endogenous IAA Determination by Liquid Chromatography Mass Spectrometry (LC/MS)
To determine the endogenous level of indole-3-acetic acid (IAA), we collected 0.1 g fresh weight each for dorsal and ventral petals. The petals were ground in liquid nitrogen then dissolved in 1 mL of cold sodium phosphate buffer (50 mM, pH 7.0) with 5 ng/mL 13 C-IAA isotope as standard. The detailed steps followed [56]. The LC/MS (Xevo TM TQ-S, Waters) analysis was conducted by the Metabolomic Core Facility at the Agricultural Biotechnology Research Center, Academia Sinica, Taiwan.

Data Availability Statement:
The data presented in this study are openly available in Supplementary Materials. The RNA-Seq data are available at the National Center for Biotechnology Information Short Reads Archive (NCBI SRA) accession number SRR16636567-SRR16636570 at the following link: https://www.ncbi.nlm.nih.gov/Traces/study/?acc=PRJNA776263, accessed on 12 November 2021.