- 1College of Life Science, Nanjing Agricultural University, Nanjing, China
- 2School of Life Sciences, Sun Yat-Sen University, Guangzhou, China
- 3College of Agro-grassland Science, Nanjing Agricultural University, Nanjing, China
Broad diversity of flowers in Fabaceae provides a good system to investigate development and evolution of floral symmetry in higher plants. Many studies have demonstrated a conserved mechanism controlling development of zygomorphic flower during last decades. However, the molecular basis of how asymmetric flower established is largely unknown. In this study, we characterized mutants named keeled wings (kw) in mungbean (Vigna radiata L.), which is a legume species with asymmetric flowers. Compared to those in the wild type plants, the lateral petals were ventralized in the kw mutants. Map-based cloning showed that KW was VrCYC3 gene in mungbean, the ortholog of Lotus japonicus CYC3 (LjCYC3) and Pisum sativum CYC3 (PsCYC3). In addition, another two CYC-like genes named VrCYC1 and VrCYC2 were identified from mungbean genome. The three CYC-like genes displayed distinct expression patterns in dorsal, lateral and ventral petals. It was found that VrCYC3 was located in nucleus. Further analysis showed that VrCYC3 had transcription activity and could interact with VrCYC1 and VrCYC2 in yeast cell. Moreover, the deletion of two amino acid residues in the R domain of VrCYC3 protein could decrease its interaction with VrCYC1 and VrCYC2 proteins. Our results suggest that LjCYC3/VrCYC3 orthologs play conserved roles determining the lateral petal shape and identity of zygomorphic flower as well as asymmetric flower in Papilionoideae.
Introduction
Floral symmetry is a distinct character with diversity in flowering plants, with three main types including actinomorphy, zygomorphy and asymmetry (Endress, 1999; Citerne et al., 2010; Hileman, 2014). Multiple studies have revealed that evolution of floral symmetry cooperates with pollinators to improve the plant productivity (Endress, 2001; Cubas, 2004; van der Kooi et al., 2021). It has been shown that CYCLOIDEA (CYC) and its paralog DICHOTOMA (DICH), encoding TCP transcription factors, are involved in the control of zygomorphic flower development in Antirrhinum majus (Luo et al., 1996; Luo et al., 1999). Furthermore, it has been found that CYC-like TCP family proteins have been recruited to play central roles in determining zygomorphy development in different species (Howarth and Donoghue, 2006; Busch and Zachgo, 2007; Broholm et al., 2008; Kim et al., 2008; Hileman and Cubas, 2009; Jabbour et al., 2009; Preston and Hileman, 2009; Rosin and Kramer, 2009; Zhang et al., 2010; Fambrini et al., 2011; Busch et al., 2012; Yang et al., 2012; Yang X. et al.,2015; Citerne et al., 2017; Dong et al., 2018; Hsu et al., 2018; Spencer and Kim, 2018; Zhao et al., 2018; Sengupta and Hileman, 2022; Sun et al., 2022).
Fabaceae is the third largest family of plants, with more than 600 genus and 18000 species (Graham and Vance, 2003). The broad diversity of flowers in Fabaceae provides a good system to investigate the development and evolution of floral symmetry. Three CYC-like TCP genes, belonging to the CYC2 clade of ECE TCP genes (Howarth and Donoghue, 2006), with overlapping and divergent functions in the establishing dorsal-ventral patterning (three types of petals: dorsal petal, lateral petal and ventral petal) have been identified in Papilionoideae species with zygomorphic flowers including Lotus japonicus and pea (Pisum sativum) (Citerne et al., 2003; Citerne et al., 2006; Feng et al., 2006; Wang et al., 2008; Woollacott and Cronk, 2018; Zhao et al., 2019). For example, LjCYC2/PsCYC2 interacts with LjCYC3/PsCYC3 to confer dorsal petal identity in L. japonicus and pea, while LjCYC3/PsCYC3 controls the lateral petal identity (Feng et al., 2006; Wang et al., 2008). The lateral petals in flowers of ljcyc3 and pscyc3 mutants are ventralized, displaying the shape of the ventral petals in the wild type plants (Feng et al., 2006; Wang et al., 2008).
However, unlike most Papilionoideae species with zygomorphic flowers, some lineages including many species in Vigna genus such as mungbean (Vigna radiata) and V. caracalla have asymmetric flowers (Endress, 1999; Etcheverry et al., 2008; Guo et al., 2019). For example, in mungbean flower, the dorsal petal unfolds, right lateral petal is closely attached to ventral petals, the tip of left lateral petal curls inward, two ventral petals fuse together, and a conspicuous spur develops on the left ventral petal (Figure 1A). Interestingly, the three types of petals are internal asymmetric and maintain the distinct dorsal-ventral patterning in these species. But, the underlying mechanism is largely unknown in the asymmetric flower.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996239/fpls-13-996239-HTML/image_m/fpls-13-996239-g001.jpg)
Figure 1 The floral phenotype of WT and kw mutants. Petals of the flowers in the wild type of plant (A), kw-1 (B) and kw-2 (C) mutants possessed dorsal-ventral (DV) differentiation. DP, dorsal petal; LP, lateral petal; VP, ventral petal. Red arrows indicated petal spurs. White arrows indicated retarded petal spurs.
In this study, we identified two allelic mutants, keeled wings (kw) in mungbean, exhibiting the deficiency in lateral petal identity. Map-based cloning showed that KW was LjCYC3 ortholog in mungbean. Our results suggest that LjCYC3/VrCYC3 orthologs play conserved roles in determining lateral petal identity of flowers with zygomorphy as well as asymmetry.
Materials and methods
Plant materials
The mungbean kw mutants (kw-1 and kw-2) were identified from the gamma ray mutagenesis population under the accession Sulu (Jiao et al., 2019). For phenotype analysis, WT plants and kw mutants in greenhouse were grown at 30/28°C (day/night temperature) with 16 h photoperiod of 250 μmol m-2 s-1 photosynthetically active radiation. The allelic tests for two mutants were carried out by crossing the kw-1 mutant with the kw-2 mutan. Ten plants of F1 generation showed the mutated phenotype.
Molecular cloning of the KW gene
Genetic mapping of KW in mungbean was conducted as described previously (Jiao et al., 2016). We crossed the kw-1 mutants with another accession AL127 and the phenotype of F1 plants were normal. The kw-1 mapping population (F2 population) was constructed from the self-pollination of F1 flowers. The 265 individual plants with mutant flower phenotype were identified in the F2 population with a total of 1150 plants. For fine mapping of KW gene in mungbean, new markers were developed according to InDels information identified from genome re-sequencing between the accessions AL127 and Sulu (Jiao et al., 2016). Primers used for genetic mapping and gene cloning were listed in Supplementary Table 1.
Scanning electron microscopy
The samples of inflorescence and floral buds were fixed with FAA and then dissected to reveal the internal floral organs. The samples were dehydrated in an alcohol series. The samples for SEM analysis were prepared as described by Chen et al. (2002). The scanning electron microscopy (SU8010, Hitachi, Tokyo, Japan) analysis was conducted as previously described (Jiao et al., 2019).
Quantitative reverse transcription PCR
Total RNA was isolated from dorsal, lateral and ventral petals from the flowers at development stage 9 of mungbean following manuals (Omega, Shanghai, China). First strand cDNA was synthesized via Takara PrimeScript™ RT reagent Kit (TaKaRa, Dalian, China) after removing genomic DNA. qRT-PCR analysis was conducted using TB GreenTM Premix ExTM RR420A (TaKaRa), and on the ABI StepOnePlus machine (Applied Biosystems, Foster City, CA, USA). The amplification condition was set as following: initial denaturation at 95°C for 10 min, followed by 40 cycles of 95°C for 15 s, 60°C for 15 s, and 72°C for 20 s. Three biological replicates with three technical repeats were conducted. VrTUB was selected as the reference gene in mungbean (Jiao et al., 2019). Relative gene expression level was determined using the 2−ΔΔCt method (Livak and Schmittgen, 2001).
Subcellular localization of VrCYC3 protein
For analysis of the subcellular localization of VrCYC3, the coding sequence (CDS) of VrCYC3 was amplified by PCR method and inserted into the pA7-YFP vector. Mungbean protoplasts were isolated and DNA-polyethylene glycol (PEG)-calcium transformation was conducted following the protocol as described by Li et al. (2019). The transformed protoplasts were observed and photographed by the confocal laser scanning microscopy (Leica, TCS SP5, Wetzlar, Germany).
Yeast two-hybrid assay
The yeast two-hybrid (Y2H) assays were carried out following the protocol (Clontech). The coding sequences of VrCYC1, VrCYC2, VrCYC3 and the mutated form VrCYC3m with 6-bp deletion (at the position 174-178 of VrCYC3 CDS) were amplified and cloned into the pGADT7 and pGBKT7 vectors respectively. The constructs of AD and BD fused to the VrCYCs were co-transformed into AH109 cells. The clones were screened on the plates with selective media SD/-Leu-Trp and SD/-Leu-Trp-Ala-His with 80 mg/mL X-α-gal.
Phylogenetic analysis of TCP family proteins
The TCP family proteins in mungbean, L. japoniucs and Arabidopsis thaliana (Supplementary Table 2) were aligned via MEGA 7.0 (Kumar et al., 2016). The phylogenetic tree was built by the neighbor-joining method with 1,000 bootstrap replicates and performed by iTOL v6 (https://itol.embl.de).
Results
Characterization of keeled wings mutants in mungbean
In order to investigate the genetic mechanism of asymmetrical flower development in mungbean, screening for variations of flowers and petals was carried out on a large-scale gamma ray mutagenesis population about 70,000 M2 lines (Jiao et al., 2019). We identified the keeled wings (kw) locus in mungbean with two mutated alleles (kw-1 and kw-2) (Figure 1), affecting dorsal-ventral patterning during flower development. In kw-1 mutant, there were ventralized petals at the lateral position and the left lateral petal possessed a spur (Figure 1B), which developed normally on the left ventral petals in the wild type plants (Figure 1A). While, in kw-2 mutant, the ventralization was somehow partial and a retarded spur developed on the left lateral petals (Figure 1C). No other obvious phenotype was observed in dorsal and ventral petals in the kw mutants, compared with the wild type plants.
Using scanning electron microscopy, flower development process was observed and compared between the wild type plants and mutants (Figure 2). It was found that the length of the lateral petals was similar to that of the ventral petals at the developmental stage 8 (Tucker, 2003), and there was no obvious difference between the wild type of plant and kw mutant at this stage (Figure 2A). However, at stage 9, when the fused ventral petals enclosed the stamens and carpels, the length of the lateral petals in the wild type of plant was two-thirds of that of the ventral petals, but the lateral petals in kw mutant mimicked that of the ventral petals (Figure 2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996239/fpls-13-996239-HTML/image_m/fpls-13-996239-g002.jpg)
Figure 2 SEM analysis of flower development in WT and kw mutant. Floral organ development at the flower development stage 8 (A) and the stage 9 (B). DP, dorsal petal; LP, lateral petal; VP, ventral petal.
Genetic mapping of KW in mungbean
Genetic analysis of kw was carried out by backcrossing the kw-1 mutant with its wild type of plant, Sulu. The flowers of F1 plants displayed a normal phenotype. An F2 population consisting of 114 plants was constructed to investigate the segregation ratio. 88 individuals displayed the wild type flowers, and 26 individuals exhibited the mutant flowers. Thus, a segregation ratio of 3:1 (χ2= 0.29<χ20.05 = 3.84) confirmed that kw mutant was controlled by a single recessive gene in mungbean.
In previous study, we have developed a set of insertion/deletion (InDel) markers to map genetic loci controlling the traits in mungbean (Jiao et al., 2016). Using the F2 mapping population from the cross between kw-1 and AL127, KW was located on chromosome 8 (Chr.08) of VC1973A genome (Kang et al., 2014). By chromosome walking approach based on developing molecular markers (Supplementary Table 1), KW was narrowed down to a region located at 40.8 Mb and 42.17 Mb between the markers ID359 and ID371 (Figure 3A). However, it was found that there was a big gap in the KW mapping region of mungbean VC1973A genome, based on the comparative analysis with aduzki bean (Vigna angularis) genome (Figure 3A; Yang K. et al., 2015). Recently, by combining short-read, long-read, and high-throughput chromatin conformation capture (Hi-C) sequencing technologies, we generated a chromosome-level genome of mungbean M5311 (data not shown). Thus, we analyzed the corresponding region on the chromosome 3 of mungbean M5311 genome to clone the KW gene (Figure 3B). Finally, KW was narrowed down to a 240 kb interval between the markers ID303 and ID306 (Supplementary Table 1), and the marker ID301 was co-segregated with the KW gene in 265 mutated plants out of a total of 1150 individuals from the F2 mapping population (Figure 3B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996239/fpls-13-996239-HTML/image_m/fpls-13-996239-g003.jpg)
Figure 3 Genetic mapping the KW gene. (A) Genetic map of the KW gene on chromosome 8 in mungbean VC1973A genome; (B) Fine mapping of KW based on the mungbean M5311 genome. The molecular marker ID301 was co-segregated with the kw-1 mutant phenotype in the F2 mapping population. (C) Schematic diagram of the VrCYC3 gene. Exons, the rectangles; Intron, the line. The CDS in exons were highlighted in grey and the non-coding sequences in white. (D) Mutations in ORF of VrCYC3. Alignment of VrCYC3 coding sequences from WT, kw-1 and kw-2. Numbers on the top of the sequences indicated the position on the open reading frame.
KW was the ortholog of L. japonicus LjCYC3 in mungbean
In the KW mapping interval, there were 25 putative genes (Supplementary Table 3). One of the candidates (VradChr3T0011489.1) encoded a TCP family protein which was the ortholog of LjCYC3 and PsCYC3 in mungbean (Figure 3C). Combined with the similar mutant phenotype, we assumed that VradChr3T0011489.1 (VrCYC3) was the likely candidate for KW in mungbean. When the sequences of VrCYC3 were analyzed, different deletions were found in both kw-1 and kw-2 mutants (Figure 3D). A 5-bp deletion was detected at the position 174-178 of the VrCYC3 ORF (Open Reading Frame) in kw-1 mutant (Figure 3D). It maycause a frameshift and premature stop codon, thus resulting a truncated protein that lacking most parts of VrCYC3 protein (Figure 4). In addition, a 6-bp deletion at the position 668-673 of the VrCYC3 ORF was detected in kw-2 mutant (Figure 3D), which would cause the translation of a mutated protein, lacking two amino acid residues, Argnine (R) and Alanine (A), in the R domain (Figure 4; Supplementary Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996239/fpls-13-996239-HTML/image_m/fpls-13-996239-g004.jpg)
Figure 4 The predicted sequences of VrCYC3 and its mutations in mungbean. Alignment of the VrCYC3 protein sequences from the wild type plant (WT), kw-1 and kw-2 mutants. Red line, TCP domain; Blue line, R domain; Red arrow indicated the truncated site in the kw-1 mutant.
Sequence analysis indicated that VrCYC3 contained two exons and one intron and encoded a 375-AA peptide (Figure 3, Figure 4). We conducted a BLASTP search for sequences with homology to VrCYC3 in our mungbean database. It was found there were 25 TCP family proteins in mungbean M5311 genome (Supplementary Table 2), which were classified into three clades by phylogenetic analysis with 26 TCP family proteins in L. japonicus (Lj1.0v1; Li et al., 2020) and 24 TCP family proteins in Arabidopsis (Figure 5). Results showed that VrCYC3 had another two close related homologs VradChr7T0021050.1 (VrCYC1) and VradChr8T0024552.1 (VrCYC2), which also belong to the ECE CYC2 clade (Howarth and Donoghue, 2006; Zhao et al., 2019), together with VradChr1T0003910.1 and VradChr4T0012762.1, forming the CYC/TB1 clade of the TCP family proteins in mungbean (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996239/fpls-13-996239-HTML/image_m/fpls-13-996239-g005.jpg)
Figure 5 Phylogenetic analysis of TCP family proteins in mungbean. The phylogenetic tree was generated using the MEGA7.0 program and displayed by iTOL. The TCP family proteins from Vigna radiata, Lotus japonicus (Feng et al., 2006; Wang et al., 2008) and Arabidopsis thaliana were used. Different clades were indicated in a specific background colour. VrCYC3 was shown in red font.
Expression and localization of VrCYC3
The expression patterns of three CYC-like genes in three types of petals of the wild type plants were analyzed by qRT-PCR. We found that both VrCYC1 and VrCYC2 were only expressed in the dorsal petals, whereas the expression of VrCYC3 could be detected in both dorsal and lateral petals, but with a relative higher level in the lateral ones (Figure 6A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996239/fpls-13-996239-HTML/image_m/fpls-13-996239-g006.jpg)
Figure 6 Gene expression of VrCYCs and subcellular localization of VrCYC3. (A) Relative expressions of VrCYC1, VrCYC2 and VrCYC3 in the dorsal petal (DP), lateral petal (VP) and ventral petal (VP) of WT. The data were means ± SD. nd, not detected. The Student’s test was used. *p < 0.05. (B) Subcellular localization of VrCYC3-GFP in mungbean protoplasts. OsARF4-mCherry was the nuclear marker. It was analyzed for green fluorescence emission, mCherry fluorescence emission and bright-field illumination.
To investigate the subcellular localization of the VrCYC3 protein in mungbean, full length of the VrCYC3 ORF was fused with the report gene GFP and then the fused construct was transiently expressed in the mungbean protoplasts through PEG4000-meidated transformation. Results showed that VrCYC3-GFP fluorescence was associated with the nucleus in the protoplast, merging with the signal of the positive control OsARF4-mCherry (Shen et al., 2010; Figure 6B).
Protein interaction between VrCYC3 and its related homologs in mungbean
It has been previously reported that the CYC-like proteins form dimer to control petal development in L. japonicus (Xu et al., 2016; Yuan et al., 2020). To test the possibility of the interactions among three CYC-like proteins in mungbean, the GAL4 binding domain (BD) and activating domain (AD)were fused with different VrCYC proteins. The interactions between VrCYC3 and VrCYC1, VrCYC3 and VrCYC2 were analyzed in yeast two-hybrid system respectively (Figure 7). The assays showed that VrCYC3 interacted with VrCYC1 and VrCYC2 in vitro. Moreover, it was found that the loss of two amino acid residues in the R domain of VrCYC3 significantly affected its interaction with VrCYC1 and VrCYC2, suggesting that the R domain could play a role in protein-protein interaction.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996239/fpls-13-996239-HTML/image_m/fpls-13-996239-g007.jpg)
Figure 7 VrCYC3 interacted with VrCYC1 and VrCYC2 in yeast. AH109 cells were co-transformed with GAL4 BD fused VrCYC and empty GAL4 AD (-) or GAL4 AD fused VrCYC. VrCYC3m, a mutated VrCYC3 form lacking two residues, Argnine (R) and Alanine (A), in R domain as in kw-2 mutant. Transformants were spotted on the plates with SD/-Trp/-Leu and SD/-Trp/-Leu/-His/-Ade/X-α-gal in 1-, 10-, 100-, and 1000-fold dilutions. The vectors of pGBKT7-53 and pGBKT7-lam, co-transformed with pGADT7-T, were used as the positive and negative control in the assays, respectively.
In addition, the results also showed that the VrCYC3 protein exhibited high expression of the reporter gene (Figure 7), suggesting that VrCYC3 had strong transcription activating activity. Conversely, the VrCYC2 and VrCYC1 fusion proteins with GAL4 BD did not exhibit this activity in the assays (Figure 7).
Discussion
In this study, two kw mutants (kw-1 and kw-2) with centralized lateral petals were isolated by screening the large-scale gamma ray mutagenesis population about 70,000 M2 lines in mungbean, a species with asymmetric flower. Interestingly, the mutations also allowed each left ventralized lateral petal to form a ectopic spur, which developed normally on the left ventral petal in the wild type plant (Figure 1). It was found that KW gene encoded the TCP family protein VrCYC3 in mungbean, the ortholog of LjCYC3 and PsCYC3. All of these genes fall into the CYC2 clade of ECE TCP genes, and not the CYC3 clade, which has been lost in Papilionoideae (Howarth and Donoghue, 2006; Zhao et al., 2019). Our previous studies have shown that in the Papilionoideae with zygomorphic flower including L. japonicus and pea, loss of function of the orthologs of LjCYC3/VrCYC3 lead to the ventralized petals at the lateral position (Feng et al., 2006; Wang et al., 2008). Thus, our results suggested that LjCYC3/VrCYC3 orthologs have functional conservation in determining the lateral petal shape and identity of flowers with zygomorphy as well as asymmetry in Papilionoideae. Consistently, it has been proposed that the asymmetric flower in some lineages in Vigna is the derived form of zygomorphy in Papilionoideae (Endress, 2001; Tucker, 2003). In addition, in mungbean flower, the three types of petals are internal (IN) asymmetric. However, the loss of function of VrCYC3 does not show the defective IN asymmetry, suggesting that there are other factors independently determining IN asymmetry in mungbean.
The establishment of floral zygormorphy is controlled by a module that relies on regulatory interaction among CYC-like TCP proteins and MYB proteins such as RADIALIS(RAD) and DIVARICATA (DIV) in the model plant A. majus (Raimundo et al., 2013). Recently, several studies suggested that the MYB transcription factors of this module have potential conserved function in the flower symmetry in other plants such as the Orchid and Lamiales (Valoroso et al., 2019; Sengupta and Hileman, 2022). In Asteraceae, the CYC-like TCP genes are required for the flower type specification (Broholm et al., 2008; Fambrini et al., 2011; Chapman et al., 2012; Juntheikki-Palovaara et al., 2014; Garcês et al., 2016; Huang et al., 2016; Ding et al., 2020; Sun et al., 2022). Interestingly, it is found that the MADS-box protein GRCD5 activates CYC-like gene expression during flower development in Gerbera hybrida, revealing an interplay between CYC-like TCP protein and MADS-box protein in the control of the flower symmetry in Asteraceae (Zhao et al., 2020; Yu, 2020).
However, in Papilionoideae such as pea and L. japonicus, none of MYB family transcription factors has been found to control flower symmetry. So far, the upstream and downstream targets or interacted partners of CYC-like TCP proteins involved in the dorsal-ventral patterning were not identified in Papilionoideae, except that a protein kinase LjRLK1 is found to interact with and modify CYC-like proteins in L. japonicus (Xu et al., 2017). In this study, we found that in the kw-2 mutant, the ventralization of lateral petals was not complete, compared with kw-1 mutant. Sequences analysis showed that a 5-bp deletion in kw-1 would lead to a frameshift and premature stop codon, lacking most parts of VrCYC3 protein (Figure 4), while the 6-bp deletion in kw-2 mutant would result in the translation of a mutated VrCYC3 protein lacking two residues in R domain which is predicted to form a coiled coil mediating protein-protein interaction (Cubas et al., 1999). Consistently, we found that the mutation in R domain decreased the affinity of the interaction of VrCYC3 with VrCYC1 and VrCYC2. Taken together, kw-2 should not be a null mutant, which could be used for the secondary mutagenesis to screen enhancers or repressors and then to identify key components involved in the CYC pathway to regulate flower symmetry in Papilionoideae.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
XL, MS, YJ, DQ, and QP performed experiment. XL and LZ analyzed the data and prepared the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the National Natural Science Foundation of China (grant no. 31200179) and the Jiangsu Agricultural Science and Technology Innovation Fund of China (CX (20)3030).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.996239/full#supplementary-material
Supplementary Figure 1 | Alignment of CYC1/2/3 proteins in Vigna radiata, Lotus japonicus and Pisum sativum. Red line, TCP domain; Blue line, R domain.
Supplementary Table 1 | Primers used in this study.
Supplementary Table 2 | Annotated genes in the kw mapping region of mungbean M5311 genome.
Supplementary Table 3 | TCP family proteins from Vigna radiata, Lotus japonicus and Arabidopsis thaliana.
References
Broholm, S. K., Tähtiharju, S., Laitinen, R. A., Albert, V. A., Teeri, T. H., Elomaa, P. (2008). A TCP domain transcription factor controls flower type specification along the radial axis of the gerbera (Asteraceae) inflorescence. Proc. Natl. Acad. Sci. U. S. A. 105, 9117–9122. doi: 10.1073/pnas.0801359105
Busch, A., Horn, S., Mühlhausen, A., Mummenhoff, K., Zachgo, S. (2012). Corolla monosymmetry: Evolution of a morphological novelty in the brassicaceae family. Mol. Biol. Evol. 29, 1241–1254. doi: 10.1093/molbev/msr297
Busch, A., Zachgo, S. (2007). Control of corolla monosymmetry in the brassicaceae Iberis amara. Proc. Natl. Acad. Sci. U. S. A. 104, 16714–16719. doi: 10.1073/pnas.0705338104
Chapman, M. A., Tang, S., Draeger, D., Nambeesan, S., Shaffer, H., Barb, J. G., et al. (2012). Genetic analysis of floral symmetry in van gogh's sunflowers reveals independent recruitment of CYCLOIDEA genes in the asteraceae. PloS Genet. 8, e1002628. doi: 10.1371/journal.pgen.1002628
Chen, C., Wang, S., Huang, H. (2002). LEUNIG has multiple functions in gynoecium development in arabidopsis. Genesis. 26, 42–54. doi: 10.1002/(SICI)1526-968X(200001)26:1<42::AID-GENE7>3.0.CO;2-J
Citerne, H., Jabbour, F., Nadot, S., Damerval, C. (2010). The evolution of floral symmetry. Adv. Bot. Res. 54, 85–137. doi: 10.1016/S0065-2296(10)54003-5
Citerne, H. L., Luo, D., Pennington, R. T., Coen, E., Cronk, Q. C. (2003). A phylogenomic investigation of CYCLOIDEA-like TCP genes in the leguminosae. Plant Physiol. 131, 1042–1053. doi: 10.1104/pp.102.016311
Citerne, H. L., Pennington, R. T., Cronk, Q. C. (2006). An apparent reversal in floral symmetry in the legume cadia is a homeotic transformation. Proc. Natl. Acad. Sci. U. S. A. 103, 12017–12020. doi: 10.1073/pnas.0600986103
Citerne, H. L., Reyes, E., Le Guilloux, M., Delannoy, E., Simonnet, F., Sauquet, H., et al. (2017). Characterization of CYCLOIDEA-like genes in proteaceae, a basal eudicot family with multiple shifts in floral symmetry. Ann. Bot. 119, 367–378. doi: 10.1093/aob/mcw219
Cubas, P. (2004). Floral zygomorphy, the recurring evolution of a successful trait. Bioessays. 26, 1175–1184. doi: 10.1002/bies.20119
Cubas, P., Lauter, N., Doebley, J., Coen, E. (1999). The TCP domain: A motif found in proteins regulating plant growth and development. Plant J. 18, 215–222. doi: 10.1046/j.1365-313x.1999.00444.x
Ding, L., Song, A., Zhang, X., Li, S., Su, J., Xia, W., et al. (2020). The core regulatory networks and hub genes regulating flower development in Chrysanthemum morifolium. Plant Mol. Bio. 103, 669–688. doi: 10.1007/s11103-020-01017-8
Dong, Y., Liu, J., Li, P. W., Li, C. Q., Lü, T. F., Yang, X., et al. (2018). Evolution of darwin's peloric gloxinia (Sinningia speciosa) is caused by a null mutation in a pleiotropic TCP gene. Mol. Biol. Evol. 35, 1901–1915. doi: 10.1093/molbev/msy090
Endress, P. K. (1999). Symmetry in flowers: diversity and evolution. Int. J. Plant Sci. 160, 3–23. doi: 10.1086/314211
Endress, P. K. (2001). Evolution of floral symmetry. Curr. Opin. Plant Biol. 4, 86–91. doi: 10.1016/s1369-5266(00)00140-0
Etcheverry, A. V., Alemán, M. M., Fleming, T. F. (2008). Flower morphology, pollination biology and mating system of the complex flower of Vigna caracalla (Fabaceae: Papilionoideae). Ann. Bot. 102, 305–316. doi: 10.1093/aob/mcn106
Fambrini, M., Salvini, M., Pugliesi, C. (2011). A transposon-mediate inactivation of a CYCLOIDEA-like gene originates polysymmetric and androgynous ray flowers in Helianthus annuus. Genetica. 139, 1521–1529. doi: 10.1007/s10709-012-9652-y
Feng, X., Zhao, Z., Tian, Z., Xu, S., Luo, Y., Cai, Z., et al. (2006). Control of petal shape and floral zygomorphy in Lotus japonicus. Proc. Natl. Acad. Sci. U. S. A. 103, 4970–4975. doi: 10.1073/pnas.0600681103
Garcês, H. M., Spencer, V. M., Kim, M. (2016). Control of floret symmetry by RAY3, SvDIV1B, and SvRAD in the capitulum of Senecio vulgaris. Plant Physiol. 171, 2055–2068. doi: 10.1104/pp.16.00395
Graham, P. H., Vance, C. P. (2003). Legumes: importance and constraints to greater use. Plant Physiol. 131, 872–877. doi: 10.1104/pp.017004
Guo, W., Zhang, X., Peng, Q., Luo, D., Jiao, K., Su, S. (2019). Love on wings, a dof family protein regulates floral vasculature in Vigna radiata. BMC Plant Biol. 19, 495. doi: 10.1186/s12870-019-2099-x
Hileman, L. C. (2014). Trends in flower symmetry evolution revealed through phylogenetic and developmental genetic advances. Philos. Trans. R. Soc Lond. Ser. B. Biol. Sci. 369, 20130348. doi: 10.1098/rstb.2013.0348
Hileman, L. C., Cubas, P. (2009). An expanded evolutionary role for flower symmetry genes. J. Biol. 8, 90. doi: 10.1186/jbiol193
Howarth, D. G., Donoghue, M. J. (2006). Phylogenetic analysis of the "ECE" (CYC/TB1) clade reveals duplications predating the core eudicots. Proc. Natl. Acad. Sci. U. S. A. 103, 9101–9106. doi: 10.1073/pnas.0602827103
Hsu, H. J., He, C. W., Kuo, W. H., Hsin, K. T., Lu, J. Y., Pan, Z. J., et al. (2018). Genetic analysis of floral symmetry transition in African violet suggests the involvement of trans-acting factor for CYCLOIDEA expression shifts. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.01008
Huang, D., Li, X., Sun, M., Zhang, T., Pan, H., Cheng, T., et al. (2016). Identification and characterization of CYC-like genes in regulation of ray floret development in Chrysanthemum morifolium. Front. Plant Sci. 7, 1633. doi: 10.3389/fpls.2016.01633
Jabbour, F., Ronse De Craene, L. P., Nadot, S., Damerval, C. (2009). Establishment of zygomorphy on an ontogenic spiral and evolution of perianth in the tribe delphinieae (Ranunculaceae). Ann. Bot. 104, 809–822. doi: 10.1093/aob/mcp162
Jiao, K. Y., Li, X., Guo, W. X., Yuan, X. X., Cui, X. Y., Chen, X. (2016). Genome re-sequencing of two accessions and fine mapping the locus of lobed leaflet margins in mungbean. Mol. Breed. 36, 128. doi: 10.1007/s11032-016-0552-1
Jiao, K., Li, X., Su, S., Guo, W., Guo, Y., Guan, Y., et al. (2019). Genetic control of compound leaf development in the mungbean (Vigna radiata l.). Hortic. Res. 6, 23. doi: 10.1038/s41438-018-0088-0
Juntheikki-Palovaara, I., Tähtiharju, S., Lan, T., Broholm, S. K., Rijpkema, A. S., Ruonala, R., et al. (2014). Functional diversification of duplicated CYC2 clade genes in regulation of inflorescence development in Gerbera hybrida (Asteraceae). Plant J. 79, 783–796. doi: 10.1111/tpj.12583
Kang, Y. J., Kim, S. K., Kim, M. Y., Lestari, P., Kim, K. H., Ha, B. K., et al. (2014). Genome sequence of mungbean and insights into evolution within vigna species. Nat. Commun. 5, 5443. doi: 10.1038/ncomms6443
Kim, M., Cui, M. L., Cubas, P., Gillies, A., Lee, K., Chapman, M. A., et al. (2008). Regulatory genes control a key morphological and ecological trait transferred between species. Science. 322, 1116–1119. doi: 10.1126/science.1164371
Kumar, S., Stecher, G., Tamura, K. (2016). MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Li, H., Jiang, F., Wu, P., Wang, K., Cao, Y. (2020). A high-quality genome sequence of model legume Lotus japonicus (MG-20) provides insights into the evolution of root nodule symbiosis. Genes. 11, 483. doi: 10.3390/genes11050483
Li, X., Liu, W., Zhuang, L., Zhu, Y., Wang, F., Chen, T., et al. (2019). BIGGER ORGANS and ELEPHANT EAR-LIKE LEAF1 control organ size and floral organ internal asymmetry in pea. J. Exp. Bot. 70, 179–191. doi: 10.1093/jxb/ery352
Livak, K. J., Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and 2–ΔΔCt method. Methods 25, 402–428. doi: 10.1006/meth.2001.1262
Luo, D., Carpenter, R., Copsey, L., Vincent, C., Clark, J., Coen, E. (1999). Control of organ asymmetry in flowers of antirrhinum. Cell. 99, 367–376. doi: 10.1016/s0092-8674(00)81523-8
Luo, D., Carpenter, R., Vincent, C., Copsey, L., Coen, E. (1996). Origin of floral asymmetry in antirrhinum. Nature. 383, 794–799. doi: 10.1038/383794a0
Preston, J. C., Hileman, L. C. (2009). Developmental genetics of floral symmetry evolution. Trends Plant Sci. 14, 147–154. doi: 10.1016/j.tplants.2008.12.005
Raimundo, J., Sobral, R., Bailey, P., Azevedo, H., Galego, L., Almeida, J., et al. (2013). A subcellular tug of war involving three MYB-like proteins underlies a molecular antagonism in antirrhinum flower asymmetry. Plant J. 75, 527–538. doi: 10.1111/tpj.12225
Rosin, F. M., Kramer, E. M. (2009). Old dogs, new tricks: Regulatory evolution in conserved genetic modules leads to novel morphologies in plants. Dev. Biol. 332, 25–35. doi: 10.1016/j.ydbio.2009.05.542
Sengupta, A., Hileman, L. C. (2022). A CYC-RAD-DIV-DRIF interaction likely predates the origin of floral monosymmetry in lamiales. Evodevo. 13, 3. doi: 10.1186/s13227-021-00187-w
Shen, C., Wang, S., Bai, Y., Wu, Y., Zhang, S., Chen, M., et al. (2010). Functional analysis of the structural domain of ARF proteins in rice (Oryza sativa l.). J. Exp. Bot. 61, 3971–3981. doi: 10.1093/jxb/erq208
Spencer, V., Kim, M. (2018). Re"CYC"ling molecular regulators in the evolution and development of flower symmetry. Semin. Cell Dev. Biol. 79, 16–26. doi: 10.1016/j.semcdb.2017.08.052
Sun, P., Bao, Y., Zhu, Y., Huang, N., Wang, X., Wu, Z. (2022). Possible role of the CYC2c gene in the cornflower-like ray floret phenotype of gaillardia cultivars. J. Plant Res. 135, 465–472. doi: 10.1007/s10265-022-01379-8
Tucker, S. C. (2003). Floral development in legumes. Plant Physiol. 131, 911–926. doi: 10.1104/pp.102.017459
Valoroso, M. C., Sobral, R., Saccone, G., Salvemini, M., Costa, M. M. R., Aceto, S. (2019). Evolutionary conservation of the orchid MYB transcription factors DIV, RAD, and DRIF. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.01359
van der Kooi, C. J., Vallejo-Marín, M., Leonhardt, S. D. (2021). Mutualisms and (A)symmetry in plant-pollinator interactions. Curr. Biol. 31, R91–R99. doi: 10.1016/j.cub.2020.11.020
Wang, Z., Luo, Y., Li, X., Wang, L., Xu, S., Yang, J., et al. (2008). Genetic control of floral zygomorphy in pea (Pisum sativum l.). Proc. Natl. Acad. Sci. U. S. A. 105, 10414–10419. doi: 10.1073/pnas.0803291105
Woollacott, C., Cronk, Q. C. B. (2018). The hooded mutant of Lathyrus odoratus (Fabaceae) is associated with a cycloidea gene mutation. Botany. 96, 47–55. doi: 10.1139/CJB-2017-0097
Xu, Z., Cheng, K., Li, X., Yang, J., Xu, S., Cao, X., et al. (2016). Transcriptional and post-transcriptional modulation of SQU and KEW activities in the control of dorsal-ventral asymmetric flower development in Lotus japonicus. Mol. Plant. 9, 722–736. doi: 10.1016/j.molp.2016.01.013
Yang, X., Pang, H. B., Liu, B. L., Qiu, Z. J., Gao, Q., Wei, L., et al. (2012). Evolution of double positive autoregulatory feedback loops in CYCLOIDEA2 clade genes is associated with the origin of floral zygomorphy. Plant Cell. 24, 1834–1847. doi: 10.1105/tpc.112.099457
Yang, K., Tian, Z., Chen, C., Luo, L., Zhao, B., Wang, Z., et al. (2015). Genome sequencing of adzuki bean (Vigna angularis) provides insight into high starch and low fat accumulation and domestication. Proc. Natl. Acad. Sci. U. S. A. 112, 13213–13218. doi: 10.1073/pnas.1420949112
Yang, X., Zhao, X. G., Li, C. Q., Liu, J., Qiu, Z. J., Dong, Y., et al. (2015). Distinct regulatory changes underlying differential expression of TEOSINTE BRANCHED1-CYCLOIDEA-PROLIFERATING CELL FACTOR genes associated with petal variations in zygomorphic flowers of petrocosmea spp. of the family gesneriaceae. Plant Physiol. 169, 2138–2151. doi: 10.1104/pp.15.01181
Yu, Y. (2020). CYCLOIDEA3 is targeted by disparate transcription factors in patterning flowers in gerbera. Plant Physiol. 184, 1214–1216. doi: 10.1104/pp.20.01315
Yuan, C., Huang, D., Yang, Y., Sun, M., Cheng, T., Wang, J., et al. (2020). CmCYC2-like transcription factors may interact with each other or bind to the promoter to regulate floral symmetry development in chrysanthemum morifolium. Plant Mol. Biol. 103, 159–171. doi: 10.1007/s11103-020-00981-5
Zhang, W., Kramer, E. M., Davis, C. C. (2010). Floral symmetry genes and the origin and maintenance of zygomorphy in a plant-pollinator mutualism. Proc. Natl. Acad. Sci. U. S. A. 107, 6388–6393. doi: 10.1073/pnas.0910155107
Zhao, Y., Broholm, S. K., Wang, F., Rijpkema, A. S., Lan, T., Albert, V. A., et al. (2020). TCP And MADS-box transcription factor networks regulate heteromorphic flower type identity in Gerbera hybrida. Plant Physiol. 184, 1455–1468. doi: 10.1104/pp.20.00702
Zhao, Z., Hu, J., Chen, S., Luo, Z., Luo, D., Wen, J., et al. (2019). Evolution of CYCLOIDEA-like genes in fabales: Insights into duplication patterns and the control of floral symmetry. Mol. Phylogenet. Evol. 132, 81–89. doi: 10.1016/j.ympev.2018.11.007
Zhao, Y., Pfannebecker, K., Dommes, A. B., Hidalgo, O., Becker, A., Elomaa, P. (2018). Evolutionary diversification of CYC/TB1-like TCP homologs and their recruitment for the control of branching and floral morphology in papaveraceae (basal eudicots). New Phytol. 220, 317–331. doi: 10.1111/nph.15289
Keywords: Fabaceae, Vigna radiata L, flower asymmetry, lateral petal identity, VrCYC3
Citation: Li X, Sun M, Jia Y, Qiu D, Peng Q and Zhuang L (2022) Genetic control of the lateral petal shape and identity of asymmetric flowers in mungbean (Vigna radiata L.). Front. Plant Sci. 13:996239. doi: 10.3389/fpls.2022.996239
Received: 17 July 2022; Accepted: 14 September 2022;
Published: 29 September 2022.
Edited by:
Elena M Kramer, Harvard University, United StatesReviewed by:
Wenheng Zhang, Virginia Commonwealth University, United StatesDianella G Howarth, St. John’s University, United States
Copyright © 2022 Li, Sun, Jia, Qiu, Peng and Zhuang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lili Zhuang, nauzll@njau.edu.cn