Abstract
The MYB transcription factors DIVARICATA (DIV), DIV-and-RAD-Interacting-Factor (DRIF), and the small interfering peptide RADIALIS (RAD) can interact, forming a regulatory module that controls different plant developmental processes. In the snapdragon Antirrhinum majus, this module, together with the TCP transcription factor CYCLOIDEA (CYC), is responsible for the establishment of floral dorsoventral asymmetry. The spatial gene expression pattern of the OitDIV, OitDRIF, and OitRAD homologs of Orchis italica, an orchid with zygomorphic flowers, has suggested a possible conserved role of these genes in bilateral symmetry of the orchid flower. Here, we have identified four DRIF genes of orchids and have reconstructed their genomic organization and evolution. In addition, we found snapdragon transcriptional cis-regulatory elements of DIV and RAD loci generally conserved within the corresponding orchid orthologues. We have tested the biochemical interactions among OitDIV, OitDRIF1, and OitRAD of O. italica, showing that OitDRIF1 can interact both with OitDIV and OitRAD, whereas OitDIV and OitRAD do not directly interact, as in A. majus. The analysis of the quantitative expression profile of these MYB genes revealed that in zygomorphic orchid flowers, the DIV, DRIF1, and RAD transcripts are present at higher levels in the lip than in lateral inner tepals, whereas in peloric orchid flowers they show similar expression levels. These results indicate that MYB transcription factors could have a role in shaping zygomorphy of the orchid flower, potentially enriching the underlying orchid developmental code.
Introduction
The MYB proteins DIVARICATA (DIV), RADIALIS (RAD), and DIV-and-RAD-Interacting-Factor (DRIF) are part of a regulatory module involved in distinct developmental processes of plants (; ). DIV and DRIF belong to ancient gene families that emerged in the green algae lineage, whereas the RAD genes are more recent as their origin can be dated back to gymnosperms (). Canonical DIV transcription factors have two MYB domains (MYBI and II) (), in contrast to RAD and DRIF, both containing a single MYB domain (; ). In addition to the N-terminal MYB domain, DRIF proteins share the conserved DUF3755 domain at the C-terminus, found only in this protein family and whose ability to bind WUSCHEL-RELATED HOMEOBOX (WOX) and KNOTTED1-LIKE HOMEOBOX (KNOX) proteins has been recently described in Populus trichocarpa (). During evolution, the MYB domain has undergone successive rearrangements resulting in the acquisition of specific interaction abilities: the MYB domain of DRIF can interact with the MYBI domain of DIV or with the MYB domain of RAD (; ; ). In such interaction module, the small RAD proteins (less than 100 amino acids in size) have an antagonistic effect on the formation of the DIV/DRIF complex and thus have been classified as small-interfering peptides (siPEP) or microproteins (Seo et al., 2011; Staudt and Wenkel, 2011; ).
The function of the DIV, DRIF, and RAD proteins has been described in different plant species, where they control distinct developmental processes. For example, RAD-like proteins regulate photomorphogenesis and floral transition of Arabidopsis thaliana (; ), DIV-like proteins are involved in sugar and hormone regulation of Oryza sativa (), and the protein complexes DIV/DRIF and RAD/DRIF control cell expansion of the fruit pericarp of Solanum lycopersicum (). However, the majority of studies regarding the DIV, DRIF, and RAD genes focused on their role in the establishment of flower zygomorphy, an evolutionary novelty that emerged several times in flowering plants from the ancestral condition of radial symmetry (; ). The first comprehensive analysis of the molecular pathway underlying floral symmetry was conducted in the snapdragon Antirrhinum majus (Figure 1), showing that mutations of the genes CYCLOIDEA (CYC), DIV, and RAD have an effect on symmetry of the flower. The TCP transcription factor CYC is expressed in the dorsal part of the flower and activates the expression of RAD (; ; ; ; ) through the interaction with 5′-GGNCCC-3′ binding sites in the RAD promoter and intron (). The DIV and DRIF genes are expressed both in the dorsal and ventral domains of the flower of A. majus (; ; ). In the ventral domain, the protein complex DIV/DRIF controls downstream genes involved in the ventralization of the flower. In the dorsal domain, the siPEP RAD binds to DRIF and prevents its interaction with DIV, thus inhibiting ventralization (). In addition to its ability to activate ventralization genes, the DIV/DRIF protein dimers can bind the sequence 5′-GATAA-3′ () within the DIV promoter, possibly autoregulating its transcriptional activity (Sengupta and Hileman, 2018).
Figure 1
The role of the CYC, DIV, and RAD genes in controlling flower bilateral symmetry outside A. majus has been reported in other Lamiales (Zhou et al., 2008; ; ; ; ; Su et al., 2017) and in Dipsacales (; ), whereas only limited knowledge about these genes is available in monocots and basal angiosperms. However, a recent study suggested the recruitment of this molecular network for the establishment of floral zygomorphy before the diversification between monocots and dicots ().
Among monocots, Orchidaceae is one of the most species-rich families, adapted to many different habitats (; ). Most of the orchid species have zygomorphic flowers sharing a common organization of the perianth into three outer and three inner tepals. Zygomorphy of the orchid flower is evident in the diversified and complex morphology of the median inner tepal (labellum or lip; Figure 2) (). The vast majority of studies concerning orchid flower development focused on MADS-box transcription factors (Salemme et al., 2011; ; ; ; ; Zhang et al., 2016; Zhang et al., 2017; Valoroso et al., 2019), with particular attention to the DEFICIENS (DEF) and AGAMOUS-LIKE 6 (AGL6) genes. The pivotal role of these genes in the evolution and formation of the orchid perianth is well explained by the “orchid code” theory and its successive modifications (; ; ; ; ). On the contrary, the MYB transcription factors are largely understudied in orchids and their potential involvement in the establishment of orchid flower symmetry has been only recently proposed in Orchis italica and Cattleya trianae (Valoroso et al., 2017; ).
Figure 2
To date, eight DIV, four RAD, and two DRIF genes have been reported in O. italica and, among them, the corresponding homologs responsible for floral symmetry of A. majus have been identified (Valoroso et al., 2017). Phylogeny and genomic organization of the orchid DIV and RAD genes has also been studied (Valoroso et al., 2017; ), whereas a description of the DRIF gene family is still missing.
The aim of the present study was to expand knowledge on DIV, RAD, and DRIF genes of orchids and to obtain more evidence supporting their involvement in the establishment of flower zygomorphy. We firstly focused on the orchid DRIF genes, searching for homologs within the orchid genomes and reconstructing their phylogeny. Then, we scanned the putative promoter and intron of the DIV and RAD genes to identify known cis-regulatory elements conserved between orchids and snapdragon. Finally, we analyzed the interaction ability of the OitDIV, OitRAD, and OitDRIF1 proteins of O. italica and examined their transcript abundance in the perianth tissues of zygomorphic and peloric orchid flowers.
Materials and Methods
Plant Material
The orchids used in this study were grown under natural light and temperature in the greenhouse of the Department of Biology of the University of Naples Federico II (Napoli, Italy). O. italica Poir. plants are part of the Orchidaceae collection of the Department of Biology. Phalaenopsis equestris (Schauer) Rchb.f. and Phalaenopsis Joy Fairy Tale (Phal. Ho’s Princess Arai × Phal. Coral Isles) are commercially available orchids (Giulio Celandroni Orchidee, San Giuliano Terme, Pisa, Italy). O. italica and P. equestris display flower zygomorphy as the second floral whorl is clearly distinguished into two lateral inner tepals and one median inner tepal (lip) (Figures 2A–E). The peloric perianth of Phalaenopsis Joy Fairy Tale shows two lip-like structures in substitution of the lateral inner tepals, conferring radial symmetry to the flower (Figures 2F, G).
Single flowers from three different plants of each orchid were collected before (single floret length, ∼1 cm) and soon after anthesis (Figure 2). The perianth tissues (outer tepals, inner lateral tepals, and lip) were dissected and stored in RNA-later (Ambion) until RNA extraction.
Sequence Retrieval and Phylogenetic Analysis
In order to identify DRIF transcripts expressed in flower tissues of O. italica, the amino acid sequences corresponding to the DUF3755 domain of the known O. italica OitDRIF1 and two proteins (GenBank accession numbers MK834277 and MK834278, respectively) (Valoroso et al., 2017) were used as queries to scan the inflorescence transcriptome of O. italica () through tBLASTn searches. Using the same approach, DRIF homologs were searched within the transcriptomes of other orchids present in the Orchidstra database () and within the genome of Apostasia shenzenica, P. equestris, and Dendrobium catenatum (; Zhang et al., 2016; Zhang et al., 2017). Recently identified DRIF nucleotide and amino acid sequences () were downloaded from SustainPine (http://www.scbi.uma.es/sustainpinedb/home_page), Monocots PLAZA (https://bioinformatics.psb.ugent.be/plaza/versions/plaza_v4_monocots/), NCBI (https://www.ncbi.nlm.nih.gov/), and Araport (https://www.araport.org/). The species and the corresponding accession numbers of the DRIF sequences used in the present work are listed in Table S1.
The amino acid sequences of the DRIF homologs identified were aligned with MAFFT () and the resulting alignment was manually adjusted. Poorly aligned positions were removed using GBLOCKS (Talavera and Castresana, 2007) and maximum likelihood (ML) phylogenetic tree was constructed with RAxML v 8.2.10 (Stamatakis, 2014) using the default settings, with 1,000 bootstrap replicates.
Analysis of Conserved Transcription Factor Binding Sites
Approximately 3 kb noncoding sequences upstream of the translation start site of the DIV and RAD genes of the orchids A. shenzenica, D. catenatum, and P. equestris were downloaded from the corresponding genomes deposited at NCBI (Table S1). These putative promoter sequences were scanned for the presence of conserved transcription factor binding sites (TFBSs) using the PlantPAN 3.0 database (). The Multiple Promoter Analysis search mode was applied to identify known conserved plant TFBSs shared within each gene group of putative promoters. Specific nucleotide motifs known as TFBSs of CYC (5′-GGNCCC-3′) (; ; Yang et al., 2012; ; Sengupta and Hileman, 2018) and DIV (5′-VGATAMSV-3′) (; Sengupta and Hileman, 2018) of A. majus were searched within the orchid RAD and DIV putative promoters, respectively. In addition, the intron sequences of the RAD genes of O. italica, P. equestris, D. catenatum (Valoroso et al., 2017), and A. shenzenica were scanned for the presence of the CYC TFBS, as described above.
Expression Analysis
Total RNA was extracted from the perianth tissues (outer tepals, inner tepals, and lip, before and after anthesis) of O. italica, P. equestris, and Phalaenopsis Joy Fairy Tale using Trizol (Ambion) followed by DNase treatment. After RNA extraction and quantification, 500 ng of total RNA from each tissue were reverse-transcribed using the Advantage RT-PCR kit (Clontech) and a mix of oligo dT and random hexamer primers.
In order to validate the nucleotide sequence of the four OitDRIF transcripts identified in the inflorescence transcriptome of O. italica, specific primer pairs were designed (Table S2) and used to amplify the cDNA of O. italica inflorescence. The amplification products obtained were cloned into pSC-A-amp/kan vector (Agilent), sequenced using the T3 and T7 primers, and analyzed using an ABI 310 Automated Sequencer (Applied Biosystems). Their sequence was compared with that of the transcripts identified in the transcriptome of O. italica.
Relative expression of the orchid DIV, RAD, and DRIF1 genes was evaluated by real-time PCR experiments, using 18S as reference gene, as previously described (; Valoroso et al., 2017). Primer pairs are listed in the Table S2. Reactions were conducted in biological triplicates and technical duplicates. Mean ± SEM was calculated for each duplicate and biological triplicate. Gene relative expression level in inner lateral tepals and lip was normalized relative to outer tepals. Two-tailed t test was conducted to assess the statistical significance of the relative expression differences between lateral inner tepals and lip of each species, before and after anthesis.
Yeast Two-Hybrid Analysis
The coding sequences (CDSs) of the OitDIV (KY089088), OitRAD (KY089097), and OitDRIF1 (MK834277) homologs of O. italica were PCR amplified using the primer pairs listed in Table S2 and 500 ng of cDNA of O. italica inflorescence. To analyse protein–protein interactions between OitDIV, OitRAD, and OitDRIF1, the GAL4-based yeast two-hybrid (Y2H) system (Matchmaker two-hybrid system; Clontech) was used. The amplified CDSs of OitDIV, OitRAD, and OitDRIF1 were cloned into bait (pGBT9) and prey (pGAD424) vectors (Clontech). All the prey and bait recombinant vector combinations were used to transform Saccharomyces cerevisiae strain AH109 (), conducting each experiment in triplicate. Plasmid presence after double yeast transformations was checked by growing cells in Synthetic Defined (SD) medium lacking tryptophan and leucine. Protein ability to interact with each other was evaluated in SD medium lacking tryptophan, leucine, and histidine. Possible transcriptional activation activity of OitDIV, OitRAD, and OitDRIF1 proteins fused to the binding domain of GAL4 (pGBT9 vector) was verified by monitoring growth of yeast transformed cells in SD medium without histidine, in the presence of 10 mM 3-aminotriazole. Empty vectors pGBT9 or pGAD424 were transformed in combination with the recombinant vectors as negative controls.
Results and Discussion
Identification and Phylogenetic Analysis of the Orchid DRIF Genes
To date, the DRIF genes of orchids have been identified only in O. italica, where the expression pattern of OitDRIF1 and OitDRIF2 was analysed in floral tissues (Valoroso et al., 2017). Evolutionary analysis has demonstrated the ancient origin of the DRIF genes: they have been found (together with the DIV genes) from green algae to angiosperms. In angiosperms, the DRIF homolog number in the examined species is generally five (). These findings led us to search for other DRIF genes expressed in O. italica, to identify their homologs in other orchids and to verify the number of DRIF genes within the available genome of orchid species, currently restricted to the only subfamilies Epidendroideae and Apostasioideae.
Among plants, the DUF3755 domain is unique to DRIF proteins (), and when we used it as query to identify other DRIFs expressed in the inflorescence transcriptome of O. italica we found two different transcripts, in addition to those previously identified, named OitDRIF3 and OitDRIF4 (accession numbers MK834279 and MK834280, respectively). All the four OitDRIF transcripts of O. italica encode for proteins containing the N-terminal MYB domain and the C-terminal DUF3755 domain. Within the genome of P. equestris and D. catenatum (both belonging to the subfamily Epidendroideae), we found four and three DRIF genes, respectively. In the genome of A. shenzenica, belonging to the basal subfamily Apostasioideae, we found three DRIF genes. All these orchid DRIF genes encode for proteins containing the MYB and the DUF3755 domain.
To cover all the five subfamilies of Orchidaceae, we scanned the transcriptomes of Ophrys sphegodes, belonging to the same subfamily of O. italica (Orchidoideae), Cypripedium formosanum (Cypripedioideae), and Vanilla planifolia (Vanilloideae) present in the database Orchidstra (), a transcriptomics collection for Orchidaceae. In C. formosanum we found four DRIF transcripts, three in O. sphegodes and two in V. planifolia. Although some of them are not full-length transcripts, missing part of the N-terminus, all contain both the MYB and DUF3755 domains. The graphical view of the amino acid alignment of the orchid DRIF proteins is reported in Figure S1, where the consensus sequences of the MYB and DUF3755 domains are shown.
We reconstructed the genomic organization of the orchid DRIF genes from the assembled genomes of P. equestris, D. catenatum, and A. shenzenica and compared it to that of the DRIF genes of A. majus, whose genome assembly has been recently released (). Based on the exon/intron number, it is possible to divide the DRIF genes into two structural types: seven exons–six introns, shared by the DRIF1-4 genes of A. majus and one DRIF gene of A. shenzenica, D. catenatum, and P. equestris; eight exons–seven introns, displayed by the DRIF5-6 genes of A. majus, two DRIF genes of A. shenzenica and D. catenatum, and three of P. equestris (Figure 3 and Table S3). Exon size is quite well conserved both among orchids and between orchids and snapdragon, whereas intron size is variable, with very large introns in orchids, reflecting a common feature of the orchid genomes due to the high number of transposable elements (Salemme et al., 2013a; Salemme et al., 2013b; ; Zhang et al., 2016; Zhang et al., 2017).
Figure 3
To understand the evolutionary relationships among the orchid DRIF proteins and the DRIFs of other plant species, we constructed the phylogenetic tree shown in Figure 4. The DRIF proteins of Tracheophyta are included in the clade I, whereas the three DRIF sequences of the moss Physcomitrella patens (Bryophyta) belong to the ancestral clade II, in agreement with the DRIF phylogeny recently described (). Within the clade I, the orchid DRIFs (in the phylogenetic tree highlighted in different shades of pink) form three distinct groups, two belonging to the subclade Ia and one to the subclade Ib (Figure 4). Based on the tree topology, the orchid DRIFs seem to have originated by duplication events predating the diversification between monocots and dicots. In fact, within both subclades Ia and Ib, the orchid branches are grouped with other monocots (O. sativa, highlighted in green) and dicots (e.g., A. majus and A. thaliana, highlighted in different shades of blue). However, lineage-specific duplications have occurred during the evolution of DRIFs in orchids, leading to the formation of two paralog groups within the subclade Ia. In particular, an orchid-specific duplication generated the paralog groups that include, among others, OitDRIF2 and 3. The orchid homologs of DRIF1 of A. majus belong to the subclade Ib, where orchid-specific duplication events seem to be absent. In fact, within the subclade Ib, a single DRIF is present for each orchid species, including OitDRIF1 of O. italica, and the orchid DRIF phylogeny reflects the evolutionary relationships existing within Orchidaceae: Apostasioideae, Cypripedioideae, and Vanilloideae are the most basal subfamilies, Epidendroideae and Orchidoideae the most derivates (). The phylogenetic relationships among DRIF proteins coincide with the division proposed based on the genomic organization of the DRIF loci examined, showing that the genes belonging to the subclade Ia have eight exons, whereas those belonging to the subclade Ib have seven exons. The genomic organization of the DRIF genes of the moss P. patens, belonging to the basal clade II, shows the presence of seven exons and six introns (Figure 3 and Table S3). The exon/intron structure shared by genes belonging to clade II and subclade Ib indicates that this genomic organization might be the ancestral condition and that the subclade Ia might have originated through a split of exon 7, followed by lineage-specific evolution of exon size.
Figure 4
Compared to the DRIF genomic organization, DIV and RAD genes have a significantly different structure with two exons and one intron (Valoroso et al., 2017). This difference supports the hypothesis previously proposed on the evolutionary origin of the DIV, RAD, and DRIF genes based on the comparison of their MYB domain (). The RAD genes might have originated through the loss of the region of the DIV genes encoding for the MYBII domain. The DIV and DRIF genes might have evolved from a common ancestral gene through lineage-specific rearrangements (e.g., duplications or gene fusion events), resulting in different exon/intron organizations and acquisition of the region encoding for the MYBII (DIV lineage) or DUF3755 domain (DRIF lineage). Phylogenetic trees of the DIV and RAD proteins of orchids and other plant species are reported in Figures S2 and S3, respectively. Both the DIV and RAD trees have been produced by the same approach used for the construction of the DRIF tree (see Materials and Methods).
Analysis of the Conserved TFBS
Some aspects of transcriptional regulation of the genes involved in floral bilateral symmetry are known for DIV and RAD of A. majus. In snapdragon, the expression of the RAD gene is activated by direct interaction of the TCP transcription factor CYC possibly through the binding to three conserved TFBSs, two located within the promoter and one within the intron of RAD (). The presence of conserved CYC TFBSs within the promoter of RAD has been recently reported also in other Lamiales species (Sengupta and Hileman, 2018).
Within the genomic sequence upstream of the translation start site of the RAD gene of P. equestris, D. catenatum, and A. shenzenica, we found distinct conserved TFBSs (Table S4), among which TCP binding sites. The sequence 5′-GGNCCN-3′, very similar to the A. majus CYC consensus binding site 5′-GGNCCC-3′, is present in the putative RAD promoter of A. shenzenica (three sites) and D. catenatum (two sites) (Table 1). Its absence in P. equestris is possibly due to the lack of a complete sequence information of the upstream region (only 958 bp are currently available) of the RAD gene in the corresponding genomic scaffold. Within the RAD intron, the sequence 5′-GGNCCN-3′ is present in D. catenatum (four sites), P. equestris (one site), and O. italica (two sites), whereas in A. shenzenica it is not present. In D. catenatum and P. equestris one of the sites exactly matches the canonical CYC TFBS of A. majus. The presence of putative CYC target sequences within the promoter and intron of the RAD gene (Table 1) may suggest a conserved direct transcriptional regulation of RAD by CYC in zygomprphic orchid flowers. In P. equestris, D. catenatum, and possibly O. italica both promoter and intron cis-regulatory motifs might be necessary to activate the transcription of RAD in the specific spatial domain linked to zygomorphy, as in A. majus. The absence of these regulatory sequences within the RAD intron of A. shenzenica is in agreement with this hypothesis, A. shenzenica being a basal orchid species with radially symmetric perianth (Zhang et al., 2017). The putative CYC binding sequence of orchids diverged to some extent from that of A. majus and other Lamiales and possibly the orchid CYC protein co-evolved to recognize slightly different sequences.
Table 1
| Gene | Species | Sequence | Position | Strand | Feature |
|---|---|---|---|---|---|
| RAD | Dendrobium catenatum | GGTCCA | −1331 | + | Putative promoter |
| GGTCCA | −2924 | + | |||
| Apostasia shenzhenica | GGTCCT | −799 | + | ||
| GGACCA | −1026 | + | |||
| GGGCCG | −2376 | + | |||
| Phalaenopsis equestris | GGACCC | 1019 | + | Intron | |
| Dendrobium catenatum | GGACCT | 123 | + | ||
| GGACCT | 218 | + | |||
| GGTCCC | 789 | + | |||
| GGACCG | 809 | + | |||
| Orchis_italica | GGCCCG | 503 | + | ||
| GGCCCG | 527 | + | |||
| DIV | Phalaenopsis equestris | AGATAAAG | −573 | − | Putative promoter |
| AGATAATA | −1551 | + | |||
| AGATAAAA | −1685 | − | |||
| Dendrobium catenatum | CGATAACC | −2210 | − | ||
| Apostasia shenzhenica | AGATAAGA | −739 | + | ||
| CGATAAGA | −1008 | + | |||
| GGATAAGA | −2898 | + |
Predicted binding sites of CYC (5′-GGNCCC-3′) and DIV (5′-VGATAMSV-3′) of Antirrhinum majus in the putative promoter and intron of RAD and in the putative promoter of DIV, respectively, in the orchids Phalaenopsis equestris, Dendrobium catenatum, and Apostasia shenzenica.
The analysis of the RAD intron was conducted also in Orchis italica. The nucleotide positions of the putative promoters are indicated with negative numbers, considering as +1 the first nucleotide of the translation start site codon ATG. The nucleotide positions of the orchid RAD introns are numbered considering as +1 the first nucleotide of the intron sequence. The binding sites exactly conserved among orchids and snapdragon are in bold.
In A. majus, two putative DIV binding sites have been identified in silico within the DIV promoter, suggesting the existence of an autoregulatory loop that maintains the transcription of DIV (Sengupta and Hileman, 2018). Within the genomic sequence upstream the translation start site of the DIV gene of P. equestris, D. catenatum, and A. shenzenica many conserved TFBSs are present, among which MYB binding sites (Table S4). The target binding sequence of DIV of A. majus 5′-VGATAMSV-3′ is present in A. shenzenica (three sites) and D. catenatum (one site), whereas the three sequences found in P. equestris have A or T instead of C or G in the seventh position (Table 1). Although the DIV binding site of orchids is only partially conserved, these results suggest that also in orchids the transcriptional activity of DIV might be regulated by a positive feedback.
Protein Interactions and Expression Pattern
In flowering plants, the involvement of the DRIF/DIV and DRIF/RAD complexes in floral zygomorphy has been demonstrated in A. majus () and inferred in a few other species (; ).
To date, the interaction among the DRIF, DIV, and RAD proteins in orchids has never been tested. We used the Y2H assay and found that in yeast the OitDRIF1 protein of O. italica can interact both with OitDIV and OitRAD, whereas OitDIV and OitRAD do not directly interact (Figure 5). This result is in agreement with the ancient evolutionary origin of this interaction module. In fact, the ability of the DIV and DRIF proteins to interact has evolved early, coincident with their origin in the green algae lineage. Later, with the emergence of the RAD genes in gymnosperms, the DRIF–RAD interaction has evolved ().
Figure 5
To obtain evidences about the possible involvement of the orchid DIV, RAD, and DRIF genes in zygomorphy of orchid flowers, we examined their expression profile in the perianth tissues of O. italica and P. equestris, both with zygomorphic flower, and in a peloric Phalaenopsis missing the bilateral symmetry due to the presence of three lips in the second floral whorl (Figure 2G). Figure 6 shows the expression levels, before and after anthesis, of the orchid DIV, RAD, and DRIF1 genes in inner tepals and lips normalized with respect to outer tepals. In the lip of O. italica and P. equestris, the expression level of the orchid DIV, RAD, and DRIF1 is significantly higher than in inner tepals after anthesis. This same pattern is observed before anthesis for the three genes in P. equestris and only for RAD in O. italica, where OitDIV and OitDRIF1 are expressed at similar levels in lip and inner tepals. This expression profile is in agreement with the previous reports about the expression of these genes in O. italica and C. trianae that suggested a possible involvement of the orchid DIV, RAD, and DRIF1 in zygomorphy of the orchid perianth (Valoroso et al., 2017; ). Interestingly, the expression levels of DIV, RAD, and DRIF1 in the peloric Phalaenopsis are similar in the lip and in the lip-like structures present in substitution of inner tepals, both before and after anthesis. The absence of a clear morphological differentiation among the structures of the second whorl in this peloric Phalaenopsis is associated with similar expression levels of the three genes, in particular of RAD. These results support the model in which in the zygomorphic flower of O. italica and P. equestris, DIV and DRIF, expressed in all the perianth organs as in A. majus, can interact in inner tepals, in which RAD expression is lower than in lip, and control ventralization. Regardless of DIV and DRIF1 levels, in the lip, RAD competes with DIV for the binding to DRIF1 and prevents the formation of the DIV–DRIF complex, thus inhibiting ventralization (Valoroso et al., 2017). Both in O. italica and P. equestris, the levels of RAD are lower in inner tepals than in lip. Consequently, in inner tepals, DIV could interact with DRIF1 and activate ventralization. On the contrary, the presence of higher levels of RAD in lip could allow the formation of the RAD/DRIF1 complex and inhibit the interaction between DIV and DRIF1, thus preventing ventralization. In the peloric Phalaenopsis, there is competition between DIV and RAD also in the lip-like structures of the second whorl and consequently ventralization is suppressed. The apparent rotation of the orchid model with respect to that of A. majus is due to resupination: the 180° rotation of the pedicel shifts the lip (a dorsal structure) to a ventral position during orchid flower development.
Figure 6
The results here obtained suggest the involvement of the orchid DIV, RAD, and DRIF genes in the zygomorphy of the orchid perianth and their conserved function, in addition to their conserved interaction ability, in species displaying flower zygomorphy.
Previous studies have demonstrated the role of the MADS-box genes in orchid flower development and various models have been proposed to explain the evolution and formation of the orchid perianth (; Salemme et al., 2013a; ; ; ; Valoroso et al., 2019). All these models propose interaction among different MADS-box transcription factors and attribute a crucial role to the expression levels of four DEF and three AGL6 genes that permit or prevent the formation of specific protein quartets that drive the formation of specific parts of the perianth (outer tepals, inner tepals, lip) (; ; ; ; ). Our findings suggest the existence of a second, MYB-based pathway underlying flower organ arrangement in orchids.
The very recent advances in functional genetic studies of orchids () open the road to investigate the potential link between these two different molecular pathways (MADS- and MYB-based) in the formation of the orchid flower and its zygomorphy and to clarify if and what kind of interaction exists, to obtain an integrated view of this complex developmental process.
Funding
This study was financially supported by grant Ricerca dipartimentale 2018 from the University of Naples Federico II and by Fundação para a Ciência e Tecnologia/Ministério da Ciência, Tecnologia e Ensino Superior through national funds (Programa de Investimento e Despesas de Desenvolvimento da Administração Central) with a project grant PTDC/BIA-PLA/1402/2014 and by FCT/MCTES/PIDDAC (Portugal) under the project PEst-OE/BIA/UI4046/2014; UID/MULTI/04046/2013.
Statements
Data availability statement
The datasets generated for this study can be found in the NCBI nucleotide, DRIF1: MK834277, DRIF2: MK834278, DRIF3: MK834279, DRIF4: MK834280.
Author contributions
MV performed the research, analyzed the data, and participated to write the paper. RS performed the research and analyzed the data. GS, MS, and MMRC analyzed the data and participated in the paper writing process. SA designed the research, analyzed the data, and wrote the paper.
Acknowledgments
We thank Prof. Giovanni Scopece for plant material and Prof. Luciano Gaudio for critical reading of the manuscript.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.01359/full#supplementary-material
References
1
AcetoS.GaudioL. (2011). The MADS and the Beauty: Genes Involved in the Development of Orchid Flowers. Curr Genomics12, 342–356. doi: 10.2174/138920211796429754
2
AcetoS.SicaM.De PaoloS.D’argenioV.CantielloP.SalvatoreF.et al. (2014). The analysis of the inflorescence miRNome of the orchid Orchis italica reveals a DEF-like MADS-box gene as a new miRNA target. PLoS One9, e97839. doi: 10.1371/journal.pone.0097839
3
Acri-Nunes-MirandaR.Mondragon-PalominoM. (2014). Expression of paralogous SEP-, FUL-, AG- and STK-like MADS-box genes in wild-type and peloric Phalaenopsis flowers. Front. Plant Sci.5:76. doi: 10.3389/fpls.2014.00076
4
AlmeidaJ.RochetaM.GalegoL. (1997). Genetic control of flower shape in Antirrhinum majus. Development124, 1387–1392.
5
BoydenG. S.DonoghueM. J.HowarthD. G. (2012). Duplications and expression of radialis-like genes in dipsacales. Int. J. Plant Sci.173, 971–983. doi: 10.1086/667626
6
CaiJ.LiuX.VannesteK.ProostS.TsaiW. C.LiuK. W.et al. (2015). The genome sequence of the orchid Phalaenopsis equestris. Nat. Genet47, 65–72. doi: 10.1038/ng.3149
7
ChaoY. T.YenS. H.YehJ. H.ChenW. C.ShihM. C. (2017). Orchidstra 2.0-a transcriptomics resource for the orchid family. Plant Cell Physiol.58 (1):e9. doi: 10.1093/pcp/pcw220
8
ChowC. N.LeeT. Y.HungY. C.LiG. Z.TsengK. C.LiuY. H.et al. (2019). PlantPAN3.0: a new and updated resource for reconstructing transcriptional regulatory networks from ChIP-seq experiments in plants. Nucleic Acids Res.47, D1155–D1163. doi: 10.1093/nar/gky1081
9
CiterneH.JabbourF.NadotS.DamervalC. (2010). The evolution of floral symmetry. Adv. Bot. Res.54 54, 85–137. doi: 10.1016/S0065-2296(10)54003-5
10
CorleyS. B.CarpenterR.CopseyL.CoenE. (2005). Floral asymmetry involves an interplay between TO and MYB transcription factors in Antirrhinum. Proc. Natl. Acad. Sci. U. S. A.102, 5068–5073. doi: 10.1073/pnas.0501340102
11
CostaM. M. R.FoxS.HannaA. I.BaxterC.CoenE. (2005). Evolution of regulatory interactions controlling floral asymmetry. Development132, 5093–5101. doi: 10.1242/dev.02085
12
CozzolinoS.WidmerA. (2005). Orchid diversity: an evolutionary consequence of deception? Trends Ecol. Evol.20, 487–494. doi: 10.1016/j.tree.2005.06.004
13
CubasP.LauterN.DoebleyJ.CoenE. (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
14
De PaoloS.GaudioL.AcetoS. (2015). Analysis of the TCP genes expressed in the inflorescence of the orchid Orchis italica. Sci. Rep.5, 16265. doi: 10.1038/srep16265
15
De PaoloS.SalveminiM.GaudioL.AcetoS. (2014). De novo transcriptome assembly from inflorescence of Orchis italica: analysis of coding and non-coding transcripts. PLoS One9, e102155. doi: 10.1371/journal.pone.0102155
16
Dirks-MulderA.ButotR.Van SchaikP.WijnandsJ. W. P. M.Van Den BergR.KrolL.et al. (2017). Exploring the evolutionary origin of floral organs of Erycina pusilla, an emerging orchid model system. BMC Evol. Biol.17:89. doi: 10.1186/s12862-017-0938-7
17
EguenT.StraubD.GraeffM.WenkelS. (2015). MicroProteins: small size - big impact. Trends Plant Sci.20, 477–482. doi: 10.1016/j.tplants.2015.05.011
18
EndressP. K. (2012). The immense diversity of floral monosymmetry and asymmetry across angiosperms. Bot. Rev.78, 345–397. doi: 10.1007/s12229-012-9106-3
19
GalegoL.AlmeidaJ. (2002). Role of DIVARICATA in the control of dorsoventral asymmetry in Antirrhinum flowers. Genes Dev.16, 880–891. doi: 10.1101/gad.221002
20
GaoY. H.ZhangD. Z.LiJ. (2015). TCP1 Modulates DWF4 Expression via Directly Interacting with the GGNCCC Motifs in the Promoter Region of DWF4 in Arabidopsis thaliana. J. Genet. Genomics42, 383–392. doi: 10.1016/j.jgg.2015.04.009
21
GarcesH. M.SpencerV. M.KimM. (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
22
GietzR. D.SchiestlR. H.WillemsA. R.WoodsR. A. (1995). Studies on the Transformation of Intact Yeast-Cells by the Liac/S-DNA/Peg Procedure. Yeast11, 355–360. doi: 10.1002/yea.320110408
23
GivnishT. J.SpalinkD.AmesM.LyonS. P.HunterS. J.ZuluagaA.et al. (2015). Orchid phylogenomics and multiple drivers of their extraordinary diversification. Proc. Biol. Sci.282:20151553. doi: 10.1098/rspb.2015.1553
24
HamaguchiA.YamashinoT.KoizumiN.KibaT.KojimaM.SakakibaraH.et al. (2008). A Small Subfamily of Arabidopsis RADIALIS-LIKE SANT/MYB Genes: A Link to HOOKLESS1-Mediated Signal Transduction during Early Morphogenesis. Biosci. Biotechnol. Biochem.72, 2687–2696. doi: 10.1271/bbb.80348
25
HowarthD. G.DonoghueM. J. (2009). Duplications and Expression of DIVARICATA-Like Genes in Dipsacales. Mol. Biol. Evol.26, 1245–1258. doi: 10.1093/molbev/msp051
26
HsuH. F.HsuW. H.LeeY. I.MaoW. T.YangJ. Y.LiJ. Y.et al. (2015). Model for perianth formation in orchids. Nat. Plants1:15046. doi: 10.1038/nplants.2015.46
27
KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol.30, 772–780. doi: 10.1093/molbev/mst010
28
KosugiS.OhashiY. (2002). DNA binding and dimerization specificity and potential targets for the TCP protein family. Plant J.30, 337–348. doi: 10.1046/j.1365-313X.2002.01294.x
29
KuiL.ChenH. T.ZhangW. X.HeS. M.XiongZ. J.ZhangY. S.et al. (2017). Building a genetic manipulation tool box for orchid biology: identification of constitutive promoters and application of CRISPR/Cas9 in the orchid, dendrobium officinale. Front. Plant Sci.7:2036. doi: 10.3389/fpls.2016.02036
30
LiC. N.ZhouY. Y.FanL. M. (2015). A novel repressor of floral transition, MEE3, an abiotic stress regulated protein, functions as an activator of FLC by binding to its promoter in Arabidopsis. Environ. Exp. Bot.113, 1–10. doi: 10.1016/j.envexpbot.2014.12.003
31
LiM.ZhangD.GaoQ.LuoY.ZhangH.MaB.et al. (2019). Genome structure and evolution of Antirrhinum majus L. Nat. Plants5, 174–183. doi: 10.1038/s41477-018-0349-9
32
LinC. S.HsuC. T.LiaoD. C.ChangW. J.ChouM. L.HuangY. T.et al. (2016). Transcriptome-wide analysis of the MADS-box gene family in the orchid Erycina pusilla. Plant Biotechnol J14, 284–298. doi: 10.1111/pbi.12383
33
LuC. A.HoT. H. D.HoS. L.YuS. M. (2002). Three novel MYB proteins with one DNA binding repeat mediate sugar and hormone regulation of alpha-amylase gene expression. Plant Cell14, 1963–1980. doi: 10.1105/tpc.001735
34
LuoD.CarpenterR.CopseyL.VincentC.ClarkJ.CoenE. (1999). Control of organ asymmetry in flowers of Antirrhinum. Cell99, 367–376. doi: 10.1016/S0092-8674(00)81523-8
35
LuoD.CarpenterR.VincentC.CopseyL.CoenE. (1996). Origin of floral asymmetry in Antirrhinum. Nature383, 794–799. doi: 10.1038/383794a0
36
MachemerK.ShaimanO.SaltsY.ShabtaiS.SobolevI.BelausovE.et al. (2011). Interplay of MYB factors in differential cell expansion, and consequences for tomato fruit development. Plant J.68, 337–350. doi: 10.1111/j.1365-313X.2011.04690.x
37
MadrigalY.AlzateJ. F.GonzalezF.Pabon-MoraN. (2019). Evolution of RADIALIS and DIVARICATA gene lineages in flowering plants with an expanded sampling in non-core eudicots. Am J. Bot. 106 (3): 1–18. doi: 10.1002/ajb2.1243
38
Mondragon-PalominoM.TheissenG. (2009). Why are orchid flowers so diverse? Reduction of evolutionary constraints by paralogues of class B floral homeotic genes. Ann. Bot.104, 583–594. doi: 10.1093/aob/mcn258
39
Mondragon-PalominoM.TheissenG. (2011). Conserved differential expression of paralogous DEFICIENS- and GLOBOSA-like MADS-box genes in the flowers of Orchidaceae: refining the ‘orchid code’. Plant J.66, 1008–1019. doi: 10.1111/j.1365-313X.2011.04560.x
40
PanZ. J.ChengC. C.TsaiW. C.ChungM. C.ChenW. H.HuJ. M.et al. (2011). The duplicated B-class MADS-box genes display dualistic characters in orchid floral organ identity and growth. Plant Cell Physiol.52, 1515–1531. doi: 10.1093/pcp/pcr092
41
PetzoldH. E.ChandaB.ZhaoC. S.RigoulotS. B.BeersE. P.BrunnerA. M. (2018). DIVARICATA AND RADIALIS INTERACTING FACTOR (DRIF) also interacts with WOX and KNOX proteins associated with wood formation in Populus trichocarpa. Plant J.93, 1076–1087. doi: 10.1111/tpj.13831
42
PrestonJ. C.KostM. A.HilemanL. C. (2009). Conservation and diversification of the symmetry developmental program among close relatives of snapdragon with divergent floral morphologies. New Phytol.182, 751–762. doi: 10.1111/j.1469-8137.2009.02794.x
43
PrestonJ. C.MartinezC. C.HilemanL. C. (2011). Gradual disintegration of the floral symmetry gene network is implicated in the evolution of a wind-pollination syndrome. Proc. Natl. Acad. Sci. U. S. A.108, 2343–2348. doi: 10.1073/pnas.1011361108
44
RaimundoJ.SobralR.BaileyP.AzevedoH.GalegoL.AlmeidaJ.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
45
RaimundoJ.SobralR.LaranjeiraS.CostaM. M. R. (2018). Successive domain rearrangements underlie the evolution of a regulatory module controlled by a small interfering peptide. Mol. Biol. Evol.35, 2873–2885. doi: 10.1093/molbev/msy178
46
ReardonW.FitzpatrickD. A.FaresM. A.NugentJ. M. (2009). Evolution of flower shape in Plantago lanceolata. Plant Mol. Biol.71, 241–250. doi: 10.1007/s11103-009-9520-z
47
ReardonW.GallagherP.NolanK. M.WrightH.Cardenosa-RubioM. C.BragaliniC.et al. (2014). Different outcomes for the MYB floral symmetry genes DIVARICATA and RADIALIS during the evolution of derived actinomorphy in Plantago. New Phytol.202, 716–725. doi: 10.1111/nph.12682
48
RudallP. J.BatemanR. M. (2002). Roles of synorganisation, zygomorphy and heterotopy in floral evolution: the gynostemium and labellum of orchids and other lilioid monocots. Biol. Rev.77, 403–441. doi: 10.1017/S1464793102005936
49
SalemmeM.SicaM.GaudioL.AcetoS. (2011). Expression pattern of two paralogs of the PI/GLO-like locus during Orchis italica (Orchidaceae, Orchidinae) flower development. Dev. Genes Evol.221, 241–246. doi: 10.1007/s00427-011-0372-6
50
SalemmeM.SicaM.GaudioL.AcetoS. (2013a). The OitaAG and OitaSTK genes of the orchid Orchis italica: a comparative analysis with other C- and D-class MADS-box genes. Mol. Biol. Rep.40, 3523–3535. doi: 10.1007/s11033-012-2426-x
51
SalemmeM.SicaM.IazzettiG.GaudioL.AcetoS. (2013b). The AP2-like gene OitaAP2 is alternatively spliced and differentially expressed in inflorescence and vegetative tissues of the orchid Orchis italica. PLoS One8, e77454. doi: 10.1371/journal.pone.0077454
52
SenguptaA.HilemanL. C. (2018). Novel traits, flower symmetry, and transcriptional autoregulation: new hypotheses from bioinformatic and experimental data. Front. Plant Sci.9:1561. doi: 10.3389/fpls.2018.01561
53
SeoP. J.HongS. Y.KimS. G.ParkC. M. (2011). Competitive inhibition of transcription factors by small interfering peptides. Trends Plant Sci.16, 541–549. doi: 10.1016/j.tplants.2011.06.001
54
StamatakisA. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics30, 1312–1313. doi: 10.1093/bioinformatics/btu033
55
StaudtA. C.WenkelS. (2011). Regulation of protein function by ‘microProteins’. Embo Rep.12, 35–42. doi: 10.1038/embor.2010.196
56
SuS. H.XiaoW.GuoW. X.YaoX. R.XiaoJ. Q.YeZ. Q.et al. (2017). The CYCLOIDEA-RADIALIS module regulates petal shape and pigmentation, leading to bilateral corolla symmetry in Torenia fournieri (Linderniaceae). New Phytol.215, 1582–1593. doi: 10.1111/nph.14673
57
TalaveraG.CastresanaJ. (2007). Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol.56, 564–577. doi: 10.1080/10635150701472164
58
ValorosoM. C.CensulloM. C.AcetoS. (2019). The MADS-box genes expressed in the inflorescence of Orchis italica (Orchidaceae). PLoS One14, e0213185. doi: 10.1371/journal.pone.0213185
59
ValorosoM. C.De PaoloS.IazzettiG.AcetoS. (2017). Transcriptome-Wide Identification and Expression Analysis of DIVARICATA- and RADIALIS-Like Genes of the Mediterranean Orchid Orchis italica. Genome Biol. Evol.9 (6): 1418–1431. doi: 10.1093/gbe/evx101
60
YangX.PangH. B.LiuB. L.QiuZ. J.GaoQ.WeiL.et al. (2012). Evolution of double positive autoregulatory feedback loops in CYCLOIDEA2 clade genes is associated with the origin of floral zygomorphy. Plant Cell24, 1834–1847. doi: 10.1105/tpc.112.099457
61
ZhangG. Q.LiuK. W.LiZ.LohausR.HsiaoY. Y.NiuS. C.et al. (2017). The Apostasia genome and the evolution of orchids. Nature549, 379–383. doi: 10.1038/nature23897
62
ZhangG. Q.XuQ.BianC.TsaiW. C.YehC. M.LiuK. W.et al. (2016). The Dendrobium catenatum Lindl. genome sequence provides insights into polysaccharide synthase, floral development and adaptive evolution. Sci. Rep.6, 19029. doi: 10.1038/srep19029
63
ZhouX. R.WangY. Z.SmithJ. F.ChenR. J. (2008). Altered expression patterns of TCP and MYB genes relating to the floral developmental transition from initial zygomorphy to actinomorphy in Bournea (Gesneriaceae). New Phytol.178, 532–543. doi: 10.1111/j.1469-8137.2008.02384.x
Summary
Keywords
DIVARICATA, RADIALIS, DRIF, MYB, Orchidaceae
Citation
Valoroso MC, Sobral R, Saccone G, Salvemini M, Costa MMR and Aceto S (2019) Evolutionary Conservation of the Orchid MYB Transcription Factors DIV, RAD, and DRIF. Front. Plant Sci. 10:1359. doi: 10.3389/fpls.2019.01359
Received
04 June 2019
Accepted
02 October 2019
Published
01 November 2019
Volume
10 - 2019
Edited by
Jen-Tsung Chen, National University of Kaohsiung, Taiwan
Reviewed by
Minsung Kim, University of Manchester, United Kingdom; Chao Bian, Beijing Genomics Institute (BGI), China
Updates
Copyright
© 2019 Valoroso, Sobral, Saccone, Salvemini, Costa and Aceto.
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: Serena Aceto, serena.aceto@unina.it
This article was submitted to Plant Development and EvoDevo, a section of the journal Frontiers in Plant Science
Disclaimer
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.