Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

The MADS-box genes expressed in the inflorescence of Orchis italica (Orchidaceae)

Abstract

The Orchidaceae family, which is one of the most species-rich flowering plant families, includes species with highly diversified and specialized flower shapes. The aim of this study was to analyze the MADS-box genes expressed in the inflorescence of Orchis italica, a wild Mediterranean orchid species. MADS-box proteins are transcription factors involved in various plant biological processes, including flower development. In the floral tissues of O. italica, 29 MADS-box genes are expressed that are classified as both class I and II. Class I MADS-box genes include one Mβ-type gene, thereby confirming the presence of this type of MADS-box genes in orchids. The class II MIKC* gene is highly expressed in the column, which is consistent with the conserved function of the MIKC* genes in gametophyte development. In addition, homologs of the SOC, SVP, ANR1, AGL12 and OsMADS32 genes are expressed. Compared with previous knowledge on class II MIKCC genes of O. italica involved in the ABCDE model of flower development, the number of class B and D genes has been confirmed. In addition, 4 class A (AP1/FUL) transcripts, 2 class E (SEP) transcripts, 2 new class C (AG) transcripts and 1 new AGL6 transcript have been identified. Within the AP1/FUL genes, the sequence divergence, relaxation of purifying selection and expression profiles suggest a possible functional diversification within these orchid genes. The detection of only two SEP transcripts in O. italica, in contrast with the 4 genes found in other orchids, suggests that only two SEP genes could be present in the subfamily Orchidoideae. The expression pattern of the MIKCC genes of O. italica indicates that low levels at the boundary of the domain of a given MADS-box gene can overlap with the expression of genes belonging to a different functional A-E class in the adjacent domain, thereby following a “fading borders” model.

Citation: Valoroso MC, Censullo MC, Aceto S (2019) The MADS-box genes expressed in the inflorescence of Orchis italica (Orchidaceae). PLoS ONE 14(3): e0213185. https://doi.org/10.1371/journal.pone.0213185

Editor: Keqiang Wu, National Taiwan University, TAIWAN

Received: January 2, 2019; Accepted: February 15, 2019; Published: March 1, 2019

Copyright: © 2019 Valoroso et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the manuscript and its Supporting Information files.

Funding: The authors received no specific funding for this work.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Among the flowering plants, the monocot family Orchidaceae is one of the most species-rich and widespread; this family has adapted to different habitats and exhibits highly specialized reproductive strategies [1]. The Orchidaceae family includes five subfamilies (Apostasioideae, Cypripedioideae, Vanilloideae, Epidendroideae and Orchidoideae) and numerous tribes and subtribes [2]. One of the most attractive orchid structures is the flower that assumes an enormous variety of shapes and colors among the species although it has a generally conserved structural organization. Three outer tepals, two lateral inner tepals and an inner median tepal (labellum or lip) define a zygomorphic perianth. Male and female reproductive tissues are fused to form the gynostemium or column, and pollen grains are located at the top of this structure. The ovary is located at the base of the column, and its maturation is triggered by pollination [3].

MADS-box genes play a crucial role in the evolution of flower architecture. This family of transcription factors is divided into two lineages, types I and II, that differ in genomic organization, developmental roles, evolutionary rate and level of functional redundancy [4].

Type I MADS-box proteins contain the MADS domain and are divided in Mα, Mβ and Mγ based on sequence divergence at the C-terminus [4]. These proteins are involved in seed, embryo and female gametophyte development [5]. Type II MADS-box proteins are the most studied lineage of MADS-box genes given their involvement in different plant developmental processes, including flower formation. These transcription factors are characterized by four domains, including three variably conserved domains and one variable domain, that form the MIKC structure. The highly conserved MADS domain has DNA-binding activity and recognizes the CArG-box motifs in the target genes [6]. The less conserved I and K domains are involved in protein-protein interactions and the formation of protein complexes [7, 8]. The variable C domain has a role in the formation of protein complexes and confers specific functional activity [8, 9].

A duplication event involving the 5’ region of the exon encoding for the K domain followed by neofunctionalization gave rise to two classes of MIKC-type genes: MIKCC and MIKC* [10]. The MIKC* genes are involved in male gametophyte development [11], whereas the MIKCC, the most studied MIKC genes, play different roles in various processes of plant growth and the establishment and maintenance of floral organs [12].

The MIKCC genes involved in flower organ formation are divided into five functional classes (from A to E), and their activity is described by the ABCDE model of flower development. The MADS-box genes involved in this model form homo- and heterodimers (floral quartets) that regulate specific expression programs in different floral whorls [13].

The ABCDE floral quartet model is applicable to different plant species, where it is well conserved [14–17]; however, analyses of non-model species have revealed differences related to the different structures of the flower [18]. For example, in orchids, the class B MADS-box genes exhibit an expression profile expanded to the first floral whorl, explaining the presence of petaloid sepals [19–21]. In addition, the orchid AP3/DEF lineage of the class B genes underwent two subsequent duplication events that have played an important role in the evolutionary origin of the current structure of the orchid flower described by the “orchid code” and the “homeotic orchid tepal” (HOT) models [22–24]. The more recent “P-code” model integrates the function of class B and AGL6 MADS-box genes to explain the formation of the orchid perianth, particularly the lip [25].

Orchis italica is an orchid species belonging to the subfamily Orchidoideae. It is one of the most widespread Mediterranean orchids with a white-purple cluster inflorescence and a lip with flaps that assume an anthropomorphic form (Fig 1). Previous studies have analyzed the structure, expression and evolution of some genes of O. italica involved in flower development. In particular, the floral meristem identity gene LFY [26], the AP2/ERF gene AP2 (class A) [27] and the MADS-box genes PI/GLO, AP3/DEF (class B), AG (class C), STK (class D) and AGL6 [20, 21, 27–33] have been assessed. However, the other MADS-box genes of O. italica have not yet been studied in contrast to other orchid species belonging to different subfamilies, mainly Epidendroideae [34–39]. Thus, the aim of this work was to identify the MADS-box genes expressed in the inflorescence of O. italica, analyze their expression in floral tissues, evaluate their evolutionary rate and compare the results with those reported in other orchids.

thumbnail
Download:
Fig 1. Orchis italica.

A, inflorescence of O. italica; B, dissected outer and lateral inner tepals (Te_out and Te_inn, respectively) and lip of O. italica.

https://doi.org/10.1371/journal.pone.0213185.g001

Materials and methods

Plant material

The Orchis italica plants used in this study (Fig 1A) are grown in the greenhouse at the Department of Biology, University of Naples Federico II. After anthesis, single florets were collected and outer tepals, lateral inner tepals, lips (Fig 1B), column and ovary (before pollination) were dissected and stored in RNALaterTM (Ambion) at −80°C for subsequent RNA extraction. Tissues from 10 flowers of two different plants were separately pooled and used for expression analysis.

Isolation of the MADS-box transcripts of O. italica

To identify the transcripts encoding for the MADS-box proteins expressed in the floral tissues, a TBLASTN search was conducted against the inflorescence transcriptome of O. italica [40] using the MADS domain of the proteins PI/GLO, AP3/DEF, AG and STK of O. italica as queries [20, 21, 30–32]. The transcripts with significant hits were virtually translated to exclude those missing the domains downstream the MADS and/or containing premature stop codons. The MADS-box coding sequences (CDSs) identified in the genome of the orchids Phalaenopsis equestris (Epidendroideae) [36], P. aphrodite [37] and Apostasia shenzenica (Apostasioideae) [35] were downloaded from GenBank and Orchidstra2 (http://orchidstra2.abrc.sinica.edu.tw/orchidstra2/pagenome.php) along with MADS-box protein and nucleotide sequences of other plant species (S1 Table).

To annotate the selected transcripts of O. italica, BLASTN and BLASTX searches were conducted using the transcripts as queries against the Viridiplantae nucleotide and protein database, respectively.

The nucleotide sequence of all the transcripts of O. italica analyzed in this study was verified by PCR amplification and sequencing using the cDNA of inflorescence of O. italica as template [33]. The sequences of the primers used are reported in S2 Table.

Phylogenetic and evolutionary analysis

The amino acid alignment of the MADS-box proteins of O. italica, Phalaenopsis and A. shenzenica was performed using MAFFT [41] and manually adjusted. The maximum likelihood (ML) tree was constructed using PhyML [42], computing the branch support values with approximate likelihood ratio tests.

The amino acid sequences of the MADS-box proteins of O. italica belonging to the AP1/FUL, SEP and AGL6 subfamilies were separately aligned to the homolog proteins of different orchids and other plant species (S1 Table), and phylogenetic trees were generated as described above. Protein-based nucleotide alignment of the orchid AP1/FUL, SEP and AGL6 subfamilies was produced using PAL2NAL [43].

The evolutionary rates of the AP1/FUL, SEP and AGL6 CDSs of O. italica and other orchids were analyzed using CODEML program from PAML v.4.9 [44]. Branch and branch-sites evolutionary models were tested, and the ratio between nonsynonymous and synonymous substitution rates (ω) was measured. The branch models assume one (one-ratio) or different (two- and three-ratios) ω values among the tree branches. The branch-site models test for positive selection on individual codons in specific branches of the tree and are compared with the respective null models that assume purifying or nearly neutral selection in the same branches. To establish which model best fits the data, a likelihood ratio test was applied comparing the alternative branch and branch-site models (allowing different ω or positive selection, respectively) with the corresponding null models.

miRNA putative target analysis

The psRNATarget web server [45] was used to predict the presence of miRNA putative target sites on the MADS-box transcripts of O. italica. The analysis (maximum expectation value 0.0) was conducted using the miRNAs of O. italica [30] as small RNA library and all the transcripts of O. italica analyzed in this study as target library.

Expression analysis

Total RNA was extracted from outer tepals (Te_out), inner lateral tepals (Te_inn), labellum (Lip), column (Co), and ovary (Ov) of O. italica after anthesis using the PureLink RNA Mini Kit (Invitrogen). Total RNA was also extracted from leaf tissue (Le). To evaluate the expression profile of the MADS-box transcripts identified in O. italica, reverse transcription and real time PCR reactions were conducted as previously reported [33, 46] using the specific primer pairs listed in S2 Table. For each transcript and tissue, the mean relative expression and standard error (SE) were calculated among three technical and two biological replicates. ANOVA followed by Tukey’s post hoc test was performed to assess the statistical significance of the mean relative expression differences among the tissues.

Results and discussion

The inflorescence transcriptome of O. italica [40] includes twenty-nine transcripts encoding for MADS-box proteins. This number is lower than that of the MADS-box genes present in the genome of P. equestris (51), D. catenatum (63) and A. shenzenica (36) [35, 36, 38] because it does not include the transcripts specifically expressed in not-floral tissues (e.g., leaf, root, stem, etc.). BLAST analysis and phylogenetic reconstruction (Fig 2) demonstrated that the MADS-box genes expressed in the inflorescence of O. italica belong to both class I and II. Excluding three class I and three class II genes (SOC, ANR and one class A), they show floral specific expression or are expressed in floral tissues at levels higher than in leaves (Fig 3).

thumbnail
Download:
Fig 2. ML tree of the orchid MADS-box proteins.

The different classes and groups of MADS-box proteins of O. italica, P. aphrodite, P. equestris and A. shenzenica are highlighted with different colors. The accessions of O. italica are noted in pink.

https://doi.org/10.1371/journal.pone.0213185.g002

thumbnail
Download:
Fig 3. Heatmap of the relative gene expression of the MADS-box genes in the floral tissues of O. italica.

The scale for each gene is independent of each other and set to 1 as the highest value. Te_out, outer tepals; Te_inn, lateral inner tepals; Co, column; Ov, ovary; Le, leaf.

https://doi.org/10.1371/journal.pone.0213185.g003

Phylogenetic and expression analysis of Class I MADS-box in O. italica

The four class I proteins of O. italica contain a conserved MADS domain (S1A Fig). Among them, two belong to the Mα type, one to Mγ and one to Mβ (Fig 2), confirming the presence of Mβ genes in orchids recently reported in P. aphrodite [37] and previously proposed to be absent from the orchid genomes [35]. The class I MADS-box transcripts of O. italica are poorly expressed in all the floral tissues (S2 Fig), as expected based on their role in embryo and endosperm maturation [5] and consistent with the low expression in the floral organs of most of the class I MADS-box genes in the orchid A. shenzenica [35]. Despite at low levels, all the class I genes identified in O. italica (especially the two Mα) are expressed in the column (at whose top are located the pollen grains) (Fig 3, S3 Fig). This pattern supports their involvement in the development of the orchid pollinia, which is consistent with the results obtained in P. aphrodite, where two Mα and three Mγ genes are specifically expressed in the pollen [37]. Their different expression in tepals, lip, ovary and leaf suggests a possible role also in other orchid floral and vegetative organs.

Phylogenetic and expression analysis of Class II MADS-box in O. italica MIKC*

The only MIKC* transcript found in the transcriptome of O. italica belongs to the P-subclade (Fig 2). This transcript encodes for a MIKC protein exhibiting a K domain longer than that of the MIKCC transcription factors (S1B Fig) as reported for the MIKC* proteins [10]. This MIKC* gene is expressed at high levels in the column (Fig 3, S3 Fig), which is consistent with the conserved function of the MIKC* genes in gametophyte development [10, 47], as hypothesized for the orchid species Erycina pusilla [34] and P. aphrodite [37].

SOC, SVP, ANR1, AGL12 and OsMADS32

Among the twenty-four MIKCC transcripts identified in the inflorescence transcriptome of O. italica, two are homologs of MADS-box genes involved in the regulation of flowering (SOC1 and SVP). Their expression in the floral tissues of O. italica is weak (S2 and S4 Figs), which is consistent with the expression profile detected in the flowers of the orchids Dendrobium [48, 49] and Erycina pusilla [34] and other plant species [50–52]. Two other transcripts only weakly expressed in the inflorescence of O. italica are homologs of MADS-box genes involved in root development (AGL12 and ANR1) [53, 54]. The AGL12 gene is present in the genome of terrestrial orchids (e.g., Apostasia) and seems to be absent in epiphytic orchids [35]. Its presence in O. italica, a terrestrial orchid species, supports its role in the formation of true roots typical of terrestrial plant species. Another MIKCC transcript expressed at low levels in the floral tissues of O. italica is a homolog of OsMADS32 (also named CFO1), a monocot-specific MADS-box gene that regulates flower organ identity [55, 56]. In rice and wheat, this gene is expressed in the early stages of the inflorescence and late seed development [50, 57], which could explain the low expression detected in the mature floral tissues of O. italica.

The MIKCC transcripts most expressed in the floral tissues of O. italica belong to the A-E classes of MADS-box genes involved in the ABCDE model of flower development (S2 Fig).

PI/GLO and AP3/DEF

Six transcript homologs of the class B genes, including two PI/GLO and four AP3/DEF, have been previously analyzed in O. italica [20, 21, 28–31]. Their expression (in particular of the AP3/DEF genes) is not restricted to the second and third floral whorls. Indeed, their expression is extended to the outer tepals (first whorl) (Fig 3), as described for orchids and other species exhibiting petaloid tepals in the first floral whorl [58]. In addition, the orchid AP3/DEF genes have been extensively studied from the evolutionary perspective and are considered the main genes responsible of the diversification of the orchid perianth, as described in the “orchid code” and HOT models [22–24].

AG and STK

Also class C (AG) and D (STK) MADS box genes of O. italica have been previously analyzed [32]. However, in the present work, two additional transcripts belonging to class C have been identified (Fig 2). Both these transcripts encode for proteins that contain AG-motifs I and II at the C-terminus (S1C Fig) and share 69% of amino acid residues. The protein encoded by the transcript OIcomp1784 exhibits 85% identity with the class C OitaAG (and 62% with the class D OitaSTK), whereas the protein encoded by the transcript OIcomp16674 exhibits 71% identity with OitaAG (and 52% with OitaSTK). A transcript similar to OIcomp16674 (84% identity) is also present in the transcriptome of Ophrys sphegodes, an orchid belonging to the Orchidoideae subfamily, as O. italica. The expression of class C and D transcripts in O. italica is high in the column and ovary (Fig 3, S5 Fig), and this pattern is consistent with the expression of the C and D genes in other orchid species, as they are involved in development of the female reproductive structure [59–61]. The newly identified transcript OIcomp16674 is also detectable in the other floral organs albeit at low levels (S5 Fig). The presence of three AG genes and one STK gene in O. italica is consistent with the copy number of class C and D MADS-box genes of Erycina pusilla [34, 39], as well as their expression pattern [39].

Excluding two AGL6 transcripts [33], class A and E genes of O. italica have never been studied. In the inflorescence of O. italica, four class A, two class E, and three AGL6 transcripts are expressed, increasing the number of AGL6 transcripts previously reported in O. italica.

AP1/FUL

The four class A transcripts of O. italica belong to the two monocot FUL-like lineages previously identified [62] (Fig 4). Within each FUL-like lineage, orchid-specific duplications expanded the copy number of the AP1/FUL genes, defining four well supported orchid FUL clades (1a, 1b, 2 and 3). Three AP1/FUL transcripts of O. italica encode for proteins exhibiting the complete FUL-like motif LPPWML at the C-terminus (S1D Fig). The protein encoded by the transcript OIcomp3679 lacks the FUL-like motif as observed in other AP1/FUL orchid proteins of species belonging to Epidendroideae [59] and Apostasioideae (S1D Fig). This finding suggests that the divergence of the C-terminus of some orchid AP1/FUL proteins is ancient given that Apostasioideae are the most basal orchid subfamily [2]. The expression of the AP1/FUL transcripts OIcomp3679, 9283 and 2508 is high in column and ovary, and the two latter transcripts are also detectable in the outer and inner tepals (Fig 3, S5 Fig). This expression profile is very similar to that of EpMADS10, 11 and 12 of E. pusilla [34, 39]. In contrast, expression of the OIcomp11046 transcript is increased in outer and inner tepals compared with column and ovary (Fig 3, S5 Fig). Although evolutionary analysis did not reveal positive selection signals within the AP1/FUL coding sequences of orchids, relaxation of purifying selection is strongly supported (Table 1, S3 Table). In fact, the ω value of the orchid AP1/FUL genes (0.23) is significantly increased compared with orchid SEP (0.11) and AGL6 (0.14) genes. This result together with a previous report of diversifying selection detected in the orchid AP1/FUL genes [39], the divergence of the C-terminus sequence of some members of this clade in orchids and the peculiar expression profile of OIcomp11046 suggest a possible functional diversification after duplication within these orchid genes.

thumbnail
Download:
Fig 4. ML tree of the Class A MADS-box proteins.

The orchid branches are highlighted in pink, and orchid FUL clades 1–3 are indicated. The accessions of O. italica are noted in pink. The numbers indicate the statistical support of the branches.

https://doi.org/10.1371/journal.pone.0213185.g004

thumbnail
Download:
Table 1. Models used and parameters estimated under different conditions for the evolutionary analysis of the orchid class A, AGL6 and class E MADS-box genes.

The parameter estimates of the branch and branch-site models were obtained from the CODEML analysis. Fore, foreground branch; back, background branch; param, parameters. The results of the LRT are reported in S3 Table.

https://doi.org/10.1371/journal.pone.0213185.t001

SEP

Two transcripts expressed in the inflorescence of O. italica encode for SEP proteins that harbor both the SEP I and II motifs (S1E Fig) and share 60% identity. In orchids, four clades of SEP genes have been reported that were generated from orchid-specific duplications within monocot SEP clades [34, 39, 59]. One SEP sequence of O. italica belongs to the clade that includes EpMADS9 and PeSEP3 (SEP clade 3), whereas the other SEP sequence is grouped with EpMADS7 and PeSEP4 (SEP clade 4) (Fig 5). BLAST search for SEP transcripts in the transcriptome of Ophrys sphegodes (Orchidoideae) [63] produced two significant hits, both exhibiting high sequence identity with the SEP transcripts identified in O. italica. In addition, a recent study identified two SEP transcripts expressed in the floral buds of Habenaria radiata (Orchidoideae) [64]. This evidence suggests that only two SEP genes are present or that only two SEP genes are expressed in the floral tissues in the subfamily Orchidoideae. Unfortunately, genome data are not currently available for orchid species belonging to Orchidoideae, rendering it difficult to discriminate between the two hypotheses. In both cases, the expression of SEP genes in O. italica is detectable in all floral organs (Fig 3, S6 Fig), and OIcomp1006 exhibits significantly increased expression in lip. This pattern is similar to that reported in wild type H. radiata [64] and other orchid species, such as the Epidendroideae P. equestris [59] and E. pusilla [34, 39], and this finding is consistent with the expression pattern in all the floral whorls of class E genes involved in the formation of all the organs of the flower.

thumbnail
Download:
Fig 5. ML tree of the class E MADS-box proteins.

The orchid branches are highlighted in pink, and the orchid SEP clades 1–4 are indicated. The accessions of O. italica are noted in pink. The numbers indicate the statistical support of the branches.

https://doi.org/10.1371/journal.pone.0213185.g005

AGL6

The three AGL6 transcripts expressed in the inflorescence of O. italica encode for proteins harboring both conserved AGL6-motifs I and II at the C-terminus (S1F Fig). In orchids, two main AGL6 clades have been reported: one includes other monocot species and one includes only orchid species [25, 37]. The three AGL6 transcripts found in O. italica belong to these two clades, and the tree topology clearly show a duplication within the orchid-specific AGL6 clade that gave rise to two different, well supported paralog groups: AGL6 clades 2 and 3 (Fig 6). Both the clades include basal orchid species (Apostasioideae, Cypripedioideae or Vanilloideae) in addition to Epidendroideae and Orchidoideae, demonstrating the ancient origin of these paralog groups. The newly identified AGL6 transcript OIcomp8204 is the paralog of OIcomp1386 [33]. These genes have an overlapping expression profile (Fig 3, S6 Fig), suggesting a redundant functional role. They are mostly expressed in lip and outer tepals, whereas OIcomp4335 exhibits high expression in the ovary and outer tepals and lower expression in lip. However, comparing the expression of these three AGL6 genes, similar levels are detected in lips, whereas OIcomp4335 expression is significantly increased in tepals (Fig 7). Orchid AGL6 genes together with AP3/DEF genes drive the formation of the orchid perianth. In particular, the “P-code” model proposes the formation of SP and L protein complexes that include combinations of different AGL6 and AP3/DEF proteins. The SP complex determines the formation of the tepals, and the L complex determines the formation of lip [25]. Given that the AGL6 transcripts of O. italica are expressed at similar levels in lip, whereas the four AP3/DEF genes exhibit different expression levels (Fig 7) [30], it is possible to hypothesize that the modulation of the expression level of the genes encoding for the different AP3/DEF proteins is responsible for the formation of the L or SP complex within the specific floral whorl.

thumbnail
Download:
Fig 6. ML tree of the AGL6 MADS-box proteins.

The orchid branches are highlighted in pink, and the orchid AGL6 clades 1–3 are indicated. The accessions of O. italica are noted in pink. The numbers indicate the statistical support of the branches.

https://doi.org/10.1371/journal.pone.0213185.g006

thumbnail
Download:
Fig 7. Expression pattern of the AGL6 and AP3/DEF MADS-box genes in the perianth of O. italica.

Each column of the AGL6 genes indicates the relative expression of 10 floral organs in two cDNA pools (10 floral organs from two different plants), both of which are amplified in triplicate. The error bars represent the standard error of the mean. Te_out, outer tepals; Te_inn, inner tepals; Le, leaf.

https://doi.org/10.1371/journal.pone.0213185.g007

The the presence of putative miRNA target sites on the twenty-nine MADS-box transcripts of O. italica identified in this study was predicted using the psRNATarget web server with a stringent parameter setting (maximum expectation 0.0). The analysis gave positive results only for the transcripts DEF2 and DEF4, known targets of the O. italica miRNAs homologs of miR5179 [30]. The other MADS-box transcripts do not have any putative target site, showing that they are not regulated by miRNAs.

Conclusions

Current advances in transcriptome and genome sequencing are highlighting the molecular programs that underlie floral evolution and development of non-model species. For many species belonging to basal angiosperms, magnoliids and basal eudicots, the “fading borders model” proposes a gradient of the expression levels of floral homeotic genes to explain flower development [65]. In particular, low expression levels at the boundary of the domain of a given MADS-box gene overlap with the expression of a different homeotic gene in the adjacent domain. The analysis of all the MADS-box genes expressed in the inflorescence of O. italica demonstrate that they follow a “fading borders” scheme and that their expression is generally conserved among orchids of different subfamilies. Given that class I and MIKC* genes are understudied compared with MIKCC genes, it will be interesting to focus forthcoming studies on these classes of MADS-box genes in orchids to understand their developmental role and evolution.

Supporting information

S1 Fig. Amino acid motifs of the orchid MADS-box proteins.

The sequences of selected orchid species are grouped according to the MADS-box type. A, class I MADS domain; B, I domain and part of the K domain of the MIKC* proteins; C-F, C-terminus of AG/STK, AP1/FUL, SEP and AGL6 proteins, respectively. The sequences of A and B are also aligned with the corresponding region of one DEF protein of O. italica. The conserved motifs of each clade of MADS-box proteins are indicated with a black line. The accessions of O. italica are noted in bold.

https://doi.org/10.1371/journal.pone.0213185.s001

(TIF)

S2 Fig. Relative expression of the MADS-box genes expressed in the inflorescence of O. italica.

Each column of the MADS-box genes indicates the relative expression of 10 floral organs in two cDNA pools (10 floral organs from two different plants), both of which are amplified in triplicate. The error bars represent the standard error of the mean. Te_out, outer tepals; Te_inn, inner tepals; Co, column; Ov, ovary.

https://doi.org/10.1371/journal.pone.0213185.s002

(TIF)

S3 Fig. Relative expression of the class I and MIKC* MADS-box genes in floral tissues of O. italica.

Each column of the class I and MIKC* genes represents the relative expression of 10 floral organs in two cDNA pools (10 floral organs from two different plants), both of which are amplified in triplicate. The error bars represent the standard error of the mean. Te_out, outer tepals; Te_inn, inner tepals; Co, column; Ov, ovary.

https://doi.org/10.1371/journal.pone.0213185.s003

(TIF)

S4 Fig. Relative expression of the SOC, SVP, ANR1, AGL12 and OsMADS32 genes in floral tissues of O. italica.

Each column of the SOC, SVP, ANR1, AGL12 and OsMADS32 genes represents the relative expression of 10 floral organs in two cDNA pools (10 floral organs from two different plants), both of which are amplified in triplicate. The error bars represent the standard error of the mean. Te_out, outer tepals; Te_inn, inner tepals; Co, column; Ov, ovary.

https://doi.org/10.1371/journal.pone.0213185.s004

(TIF)

S5 Fig. Relative expression of the AG and class A MADS-box genes in floral tissues of O. italica.

Each column of the class C AG and class A AP1/FUL genes indicates the relative expression of 10 floral organs in two cDNA pools (10 floral organs from two different plants), both of which are amplified in triplicate. The error bars represent the standard error of the mean. Te_out, outer tepals; Te_inn, inner tepals; Co, column; Ov, ovary.

https://doi.org/10.1371/journal.pone.0213185.s005

(TIF)

S6 Fig. Relative expression of the AGL6 and class E MADS-box genes in floral tissues of O. italica.

Each column of the AGL6 and the class E SEP genes shows the relative expression of 10 floral organs in two cDNA pools (10 floral organs from two different plants), both of which are amplified in triplicate. The error bars represent the standard error of the mean. Te_out, outer tepals; Te_inn, inner tepals; Co, column; Ov, ovary.

https://doi.org/10.1371/journal.pone.0213185.s006

(TIF)

S1 Table. List of sequences used in the alignments, phylogenetic and evolutionary analyses.

https://doi.org/10.1371/journal.pone.0213185.s007

(XLSX)

S2 Table. List of the primers used in the present work.

https://doi.org/10.1371/journal.pone.0213185.s008

(XLSX)

S3 Table. Likelihood ratio statistics for the comparison of the evolutionary models of the orchid class A, AGL6 and class E coding regions.

df, degrees of freedom; ns, not significant.

https://doi.org/10.1371/journal.pone.0213185.s009

(XLSX)

Acknowledgments

The authors are grateful to Prof. Giovanni Scopece for plant material and help in flower dissection and to Prof. Luciano Gaudio for critical reading of the manuscript.

References

  1. 1. Cozzolino S, Widmer A. Orchid diversity: an evolutionary consequence of deception? Trends Ecol Evol. 2005;20(9):487–94. WOS:000232128300008. pmid:16701425
  2. 2. Gustafsson ALS, Verola CF, Antonelli A. Reassessing the temporal evolution of orchids with new fossils and a Bayesian relaxed clock, with implications for the diversification of the rare South American genus Hoffmannseggella (Orchidaceae: Epidendroideae). Bmc Evol Biol. 2010;10. Artn 177 10.1186/1471-2148-10-177. WOS:000280368000001.
  3. 3. Rudall PJ, Bateman RM. Roles of synorganisation, zygomorphy and heterotopy in floral evolution: the gynostemium and labellum of orchids and other lilioid monocots. Biol Rev. 2002;77(3):403–41. WOS:000178058500004. pmid:12227521
  4. 4. Parenicova L, de Folter S, Kieffer M, Horner DS, Favalli C, Busscher J, et al. Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis: new openings to the MADS world. Plant Cell. 2003;15(7):1538–51. pmid:12837945; PubMed Central PMCID: PMC165399.
  5. 5. Masiero S, Colombo L, Grini PE, Schnittger A, Kater MM. The emerging importance of type I MADS box transcription factors for plant reproduction. The Plant cell. 2011;23(3):865–72. pmid:21378131; PubMed Central PMCID: PMC3082269.
  6. 6. Tilly JJ, Allen DW, Jack T. The CArG boxes in the promoter of the Arabidopsis floral organ identity gene APETALA3 mediate diverse regulatory effects. Development. 1998;125(9):1647–57. pmid:9521903.
  7. 7. Riechmann JL, Krizek BA, Meyerowitz EM. Dimerization specificity of Arabidopsis MADS domain homeotic proteins APETALA1, APETALA3, PISTILLATA, and AGAMOUS. Proceedings of the National Academy of Sciences of the United States of America. 1996;93(10):4793–8. pmid:8643482; PubMed Central PMCID: PMC39358.
  8. 8. Kaufmann K, Melzer R, Theissen G. MIKC-type MADS-domain proteins: structural modularity, protein interactions and network evolution in land plants. Gene. 2005;347(2):183–98. pmid:15777618.
  9. 9. Riechmann JL, Meyerowitz EM. MADS domain proteins in plant development. Biological chemistry. 1997;378(10):1079–101. pmid:9372178.
  10. 10. Kwantes M, Liebsch D, Verelst W. How MIKC* MADS-box genes originated and evidence for their conserved function throughout the evolution of vascular plant gametophytes. Molecular biology and evolution. 2012;29(1):293–302. pmid:21813465.
  11. 11. Zobell O, Faigl W, Saedler H, Munster T. MIKC* MADS-box proteins: conserved regulators of the gametophytic generation of land plants. Molecular biology and evolution. 2010;27(5):1201–11. pmid:20080864.
  12. 12. Aceto S, Gaudio L. The MADS and the Beauty: Genes Involved in the Development of Orchid Flowers. Curr Genomics. 2011;12(5):342–56. pmid:22294877; PubMed Central PMCID: PMC3145264.
  13. 13. Theissen G, Melzer R, Rumpler F. MADS-domain transcription factors and the floral quartet model of flower development: linking plant development and evolution. Development. 2016;143(18):3259–71. WOS:000393450700007. pmid:27624831
  14. 14. Ambrose BA, Lerner DR, Ciceri P, Padilla CM, Yanofsky MF, Schmidt RJ. Molecular and genetic analyses of the silky1 gene reveal conservation in floral organ specification between eudicots and monocots. Molecular cell. 2000;5(3):569–79. pmid:10882141.
  15. 15. Whipple CJ, Ciceri P, Padilla CM, Ambrose BA, Bandong SL, Schmidt RJ. Conservation of B-class floral homeotic gene function between maize and Arabidopsis. Development. 2004;131(24):6083–91. pmid:15537689.
  16. 16. Whipple CJ, Zanis MJ, Kellogg EA, Schmidt RJ. Conservation of B class gene expression in the second whorl of a basal grass and outgroups links the origin of lodicules and petals. Proceedings of the National Academy of Sciences of the United States of America. 2007;104(3):1081–6. pmid:17210918; PubMed Central PMCID: PMC1783367.
  17. 17. Ferrario S, Immink RG, Angenent GC. Conservation and diversity in flower land. Current opinion in plant biology. 2004;7(1):84–91. pmid:14732446.
  18. 18. Irish VF, Litt A. Flower development and evolution: gene duplication, diversification and redeployment. Current opinion in genetics & development. 2005;15(4):454–60. pmid:15964755.
  19. 19. Kanno A, Nakada M, Akita Y, Hirai M. Class B gene expression and the modified ABC model in nongrass monocots. TheScientificWorldJournal. 2007;7:268–79. pmid:17334618; PubMed Central PMCID: PMC5901336.
  20. 20. Cantone C, Sica M, Gaudio L, Aceto S. The OrcPI locus: Genomic organization, expression pattern, and noncoding regions variability in Orchis italica (Orchidaceae) and related species. Gene. 2009;434(1–2):9–15. WOS:000264589800002. pmid:19162144
  21. 21. Salemme M, Sica M, Gaudio L, Aceto S. Expression pattern of two paralogs of the PI/GLO-like locus during Orchis italica (Orchidaceae, Orchidinae) flower development. Development genes and evolution. 2011;221(4):241–6. pmid:21892737.
  22. 22. Mondragon-Palomino M, Theissen G. Why are orchid flowers so diverse? Reduction of evolutionary constraints by paralogues of class B floral homeotic genes. Ann Bot-London. 2009;104(3):583–94. WOS:000268807600020. pmid:19141602
  23. 23. Mondragon-Palomino M, Theissen G. Conserved differential expression of paralogous DEFICIENS- and GLOBOSA-like MADS-box genes in the flowers of Orchidaceae: refining the 'orchid code'. Plant J. 2011;66(6):1008–19. WOS:000292103900008. pmid:21435045
  24. 24. Pan ZJ, Cheng CC, Tsai WC, Chung MC, Chen WH, Hu JM, et al. The duplicated B-class MADS-box genes display dualistic characters in orchid floral organ identity and growth. Plant Cell Physiol. 2011;52(9):1515–31. pmid:21757456.
  25. 25. Hsu HF, Hsu WH, Lee YI, Mao WT, Yang JY, Li JY, et al. Model for perianth formation in orchids. Nat Plants. 2015;1(5). Artn 15046 WOS:000364398500001.
  26. 26. Montieri S, Gaudio L, Aceto S. Isolation of the LFY/FLO homologue in Orchis italica and evolutionary analysis in some European orchids. Gene. 2004;333:101–9. WOS:000222074800012. pmid:15177685
  27. 27. Salemme M, Sica M, Iazzetti G, Gaudio L, Aceto S. The AP2-like gene OitaAP2 is alternatively spliced and differentially expressed in inflorescence and vegetative tissues of the orchid Orchis italica. PLoS One. 2013;8(10):e77454. pmid:24204832; PubMed Central PMCID: PMC3804621.
  28. 28. Aceto S, Cantone C, Chiaiese P, Ruotolo G, Sica M, Gaudio L. Isolation and Phylogenetic Footprinting Analysis of the 5'-Regulatory Region of the Floral Homeotic Gene OrcPI from Orchis italica (Orchidaceae). J Hered. 2010;101(1):124–31. WOS:000272687700016. pmid:19861638
  29. 29. Aceto S, Montieri S, Sica M, Gaudio L. Molecular evolution of the OrcPI locus in natural populations of Mediterranean orchids. Gene. 2007;392(1–2):299–305. WOS:000246078300034. pmid:17321699
  30. 30. Aceto S, Sica M, De Paolo S, D'Argenio V, Cantiello P, Salvatore F, et al. The analysis of the inflorescence miRNome of the orchid Orchis italica reveals a DEF-like MADS-box gene as a new miRNA target. PloS one. 2014;9(5):e97839. pmid:24832004; PubMed Central PMCID: PMC4022656.
  31. 31. Cantone C, Gaudio L, Aceto S. The PI/GLO-like locus in orchids: duplication and purifying selection at synonymous sites within Orchidinae (Orchidaceae). Gene. 2011;481(1):48–55. pmid:21549812.
  32. 32. Salemme M, Sica M, Gaudio L, Aceto S. The OitaAG and OitaSTK genes of the orchid Orchis italica: a comparative analysis with other C- and D-class MADS-box genes. Molecular biology reports. 2013;40(5):3523–35. pmid:23277396.
  33. 33. Valoroso MC, De Paolo S, Iazzetti G, Aceto S. Transcriptome-Wide Identification and Expression Analysis of DIVARICATA- and RADIALIS-Like Genes of the Mediterranean Orchid Orchis italica. Genome biology and evolution. 2017;9(6). pmid:28541415; PubMed Central PMCID: PMC5499889.
  34. 34. Lin CS, Hsu CT, Liao DC, Chang WJ, Chou ML, Huang YT, et al. Transcriptome-wide analysis of the MADS-box gene family in the orchid Erycina pusilla. Plant biotechnology journal. 2016;14(1):284–98. pmid:25917508.
  35. 35. Zhang GQ, Liu KW, Li Z, Lohaus R, Hsiao YY, Niu SC, et al. The Apostasia genome and the evolution of orchids. Nature. 2017;549(7672):379–83. pmid:28902843.
  36. 36. Cai J, Liu X, Vanneste K, Proost S, Tsai WC, Liu KW, et al. The genome sequence of the orchid Phalaenopsis equestris. Nature genetics. 2015;47(1):65–72. pmid:25420146.
  37. 37. Chao YT, Chen WC, Chen CY, Ho HY, Yeh CH, Kuo YT, et al. Chromosome-level assembly, genetic and physical mapping of Phalaenopsis aphrodite genome provides new insights into species adaptation and resources for orchid breeding. Plant biotechnology journal. 2018. pmid:29704444.
  38. 38. Zhang GQ, Xu Q, Bian C, Tsai WC, Yeh CM, Liu KW, et al. The Dendrobium catenatum Lindl. genome sequence provides insights into polysaccharide synthase, floral development and adaptive evolution. Scientific reports. 2016;6:19029. pmid:26754549; PubMed Central PMCID: PMC4709516.
  39. 39. Dirks-Mulder A, Butot R, van Schaik P, Wijnands JWPM, van den Berg R, Krol L, et al. Exploring the evolutionary origin of floral organs of Erycina pusilla, an emerging orchid model system. Bmc Evol Biol. 2017;17. ARTN 89 WOS:000397338900005.
  40. 40. De Paolo S, Salvemini M, Gaudio L, Aceto S. De novo transcriptome assembly from inflorescence of Orchis italica: analysis of coding and non-coding transcripts. PloS one. 2014;9(7):e102155. pmid:25025767; PubMed Central PMCID: PMC4099010.
  41. 41. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular biology and evolution. 2013;30(4):772–80. pmid:23329690; PubMed Central PMCID: PMC3603318.
  42. 42. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic biology. 2010;59(3):307–21. pmid:20525638.
  43. 43. Suyama M, Torrents D, Bork P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic acids research. 2006;34(Web Server issue):W609–12. pmid:16845082; PubMed Central PMCID: PMC1538804.
  44. 44. Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Molecular biology and evolution. 2007;24(8):1586–91. pmid:17483113.
  45. 45. Dai XB, Zhuang ZH, Zhao PXC. psRNATarget: a plant small RNA target analysis server (2017 release). Nucleic acids research. 2018;46(W1):W49–W54. WOS:000438374100009. pmid:29718424
  46. 46. De Paolo S, Gaudio L, Aceto S. Analysis of the TCP genes expressed in the inflorescence of the orchid Orchis italica. Scientific reports. 2015;5:16265. pmid:26531864; PubMed Central PMCID: PMC4632031.
  47. 47. Liu Y, Cui S, Wu F, Yan S, Lin X, Du X, et al. Functional conservation of MIKC*-Type MADS box genes in Arabidopsis and rice pollen maturation. The Plant cell. 2013;25(4):1288–303. pmid:23613199; PubMed Central PMCID: PMC3663268.
  48. 48. Ding L, Wang Y, Yu H. Overexpression of DOSOC1, an ortholog of Arabidopsis SOC1, promotes flowering in the orchid Dendrobium Chao Parya Smile. Plant & cell physiology. 2013;54(4):595–608. pmid:23396600.
  49. 49. Liu XR, Pan T, Liang WQ, Gao L, Wang XJ, Li HQ, et al. Overexpression of an Orchid (Dendrobium nobile) SOC1/TM3-Like Ortholog, DnAGL19, in Arabidopsis Regulates HOS1-FT Expression. Frontiers in plant science. 2016;7:99. pmid:26904066; PubMed Central PMCID: PMC4746357.
  50. 50. Arora R, Agarwal P, Ray S, Singh AK, Singh VP, Tyagi AK, et al. MADS-box gene family in rice: genome-wide identification, organization and expression profiling during reproductive development and stress. Bmc Genomics. 2007;8. Artn 242 WOS:000248702600001. pmid:17640358
  51. 51. Diaz-Riquelme J, Lijavetzky D, Martinez-Zapater JM, Carmona MJ. Genome-Wide Analysis of MIKCC-Type MADS Box Genes in Grapevine. Plant Physiol. 2009;149(1):354–69. WOS:000262261500040. pmid:18997115
  52. 52. Hu LF, Liu SQ. Genome-wide analysis of the MADS-box gene family in cucumber. Genome. 2012;55(3):245–56. WOS:000301043300007. pmid:22376137
  53. 53. Tapia-Lopez R, Garcia-Ponce B, Dubrovsky JG, Garay-Arroyo A, Perez-Ruiz RV, Kim SH, et al. An AGAMOUS-related MADS-box gene, XAL1 (AGL12), regulates root meristem cell proliferation and flowering transition in Arabidopsis. Plant Physiol. 2008;146(3):1182–92. WOS:000256412200042. pmid:18203871
  54. 54. Zhang HM, Forde BG. An Arabidopsis MADS box gene that controls nutrient-induced changes in root architecture. Science. 1998;279(5349):407–9. WOS:000071570800053. pmid:9430595
  55. 55. Becker A, Theissen G. The major clades of MADS-box genes and their role in the development and evolution of flowering plants. Mol Phylogenet Evol. 2003;29(3):464–89. WOS:000186738000008. pmid:14615187
  56. 56. Sang XC, Li YF, Luo ZK, Ren DY, Fang LK, Wang N, et al. CHIMERIC FLORAL ORGANS1, Encoding a Monocot-Specific MADS Box Protein, Regulates Floral Organ Identity in Rice. Plant Physiol. 2012;160(2):788–807. WOS:000309523800018. pmid:22891238
  57. 57. Zhao T, Ni ZF, Dai Y, Yao YY, Nie XL, Sun QX. Characterization and expression of 42 MADS-box genes in wheat (Triticum aestivum L.). Mol Genet Genomics. 2006;276(4):334–50. WOS:000241553200003. pmid:16858583
  58. 58. Dodsworth S., Petal Sepal, or Tepal? B-Genes and Monocot Flowers. Trends Plant Sci. 2017;22(1):8–10. WOS:000392892400005. pmid:27894712
  59. 59. Acri-Nunes-Miranda R, Mondragon-Palomino M. Expression of paralogous SEP-, FUL-, AG- and STK-like MADS-box genes in wild-type and peloric Phalaenopsis flowers. Frontiers in plant science. 2014;5. Artn 76 WOS:000332670000001. pmid:24659990
  60. 60. Hsu HF, Hsieh WP, Chen MK, Chang YY, Yang CH. C/D Class MADS Box Genes from Two Monocots, Orchid (Oncidium Gower Ramsey) and Lily (Lilium longiflorum), Exhibit Different Effects on Floral Transition and Formation in Arabidopsis thaliana. Plant and Cell Physiology. 2010;51(6):1029–45. WOS:000279534200016. pmid:20395287
  61. 61. Skipper M, Johansen LB, Pedersen KB, Frederiksen S, Johansen BB. Cloning and transcription analysis of an AGAMOUS- and SEEDSTICK ortholog in the orchid Dendrobium thyrsiflorum (Reichb. f.). Gene. 2006;366(2):266–74. WOS:000236312600008. pmid:16236468
  62. 62. Litt A, Irish VF. Duplication and diversification in the APETALA1/FRUITFULL floral homeotic gene lineage: Implications for the evolution of floral development. Genetics. 2003;165(2):821–33. WOS:000186224900032. pmid:14573491
  63. 63. Sedeek KEM, Qi WH, Schauer MA, Gupta AK, Poveda L, Xu SQ, et al. Transcriptome and Proteome Data Reveal Candidate Genes for Pollinator Attraction in Sexually Deceptive Orchids. PloS one. 2013;8(5). ARTN e64621 WOS:000319725500058. pmid:23734209
  64. 64. Mitoma M, Kanno A. The Greenish Flower Phenotype of Habenaria radiata (Orchidaceae) Is Caused by a Mutation in the SEPALLATA-Like MADS-Box Gene HrSEP-1. Frontiers in plant science. 2018;9. Artn 831 WOS:000435670800004. pmid:29971084
  65. 65. Chanderbali AS, Berger BA, Howarth DG, Soltis PS, Soltis DE. Evolving Ideas on the Origin and Evolution of Flowers: New Perspectives in the Genomic Era. Genetics. 2016;202(4):1255–65. WOS:000373959100004. pmid:27053123