Positive selection and ancient duplications in the evolution of class B floral homeotic genes of orchids and grasses

Research advances on the gene regulation of floral development and color in orchids

Orchidaceae is one of the largest monocotyledon families and contributes significantly to worldwide biodiversity, with value in the fields of landscaping, medicine, and ecology. The diverse phenotypes and vibrant colors of orchid floral organs make them excellent research objects for investigating flower development and pigmentation. In recent years, a number of orchid genomes have been published, laying the molecular foundation for revealing flower development and color presentation. In this article, we review transcription factors, the structural genes responsible for the floral pigment synthesis pathways, the molecular mechanisms of flower morphogenesis, and the potential relationship between flower type and flower color. This study provides a theoretical reference for the research on molecular mechanisms related to flower morphogenesis and color presentation, genetic improvement, and new variety creation in orchids.

Molecular mechanism underlying pseudopeloria in Habenaria radiata (Orchidaceae)

Habenaria radiata (Orchidaceae) has two whorls of perianth, comprising three greenish sepals, two white petals and one lip (labellum). By contrast, the pseudopeloric (with a decreased degree of zygomorphy) mutant cultivar of H. radiata, 'Hishou', has changes in the identities of the dorsal sepal to a petaloid organ and the two ventral sepals to lip-like organs. Here, we isolated four DEFICIENS-like and two AGL6-like genes from H. radiata, and characterized their expression. Most of these genes revealed similar expression patterns in the wild type and in the 'Hishou' cultivar, except HrDEF-C3. The HrDEF-C3 gene was expressed in petals and lip in the wild type but was ectopically expressed in sepal, petals, lip, leaf, root and bulb in 'Hishou'. Sequence analysis of the HrDEF-C3 loci revealed that the 'Hishou' genome harbored two types of HrDEF-C3 genes: one identical to wild-type HrDEF-C3 and the other carrying a retrotransposon insertion in its promoter. Genetic linkage analysis of the progeny derived from an intraspecific cross between 'Hishou' and the wild type demonstrated that the mutant pseudopeloric trait was dominantly inherited and was linked to the HrDEF-C3 gene carrying the retrotransposon. These results indicate that the pseudopeloric phenotype is caused by retrotransposon insertion in the HrDEF-C3 promoter, resulting in the ectopic expression of HrDEF-C3. As the expression of HrAGL6-C2 was limited to lateral sepals and lip, the overlapping expression of HrDEF-C3 and HrAGL6-C2 is likely to be responsible for the sepal to lip-like identity in the lateral sepals of the 'Hishou' cultivar.

A Dead Gene Walking: Convergent Degeneration of a Clade of MADS-Box Genes in Crucifers

Genes are "born," and eventually they "die." These processes shape the phenotypic evolution of organisms and are hence of great biological interest. If genes die in plants, they generally do so quite rapidly. Here, we describe the fate of GOA-like genes that evolve in a dramatically different manner. GOA-like genes belong to the subfamily of Bsister genes of MIKC-type MADS-box genes. Typical MIKC-type genes encode conserved transcription factors controlling plant development. We show that ABS-like genes, a clade of Bsister genes, are indeed highly conserved in crucifers (Brassicaceae) maintaining the ancestral function of Bsister genes in ovule and seed development. In contrast, their closest paralogs, the GOA-like genes, have been undergoing convergent gene death in Brassicaceae. Intriguingly, erosion of GOA-like genes occurred after millions of years of coexistence with ABS-like genes. We thus describe Delayed Convergent Asymmetric Degeneration, a so far neglected but possibly frequent pattern of duplicate gene evolution that does not fit classical scenarios. Delayed Convergent Asymmetric Degeneration of GOA-like genes may have been initiated by a reduction in the expression of an ancestral GOA-like gene in the stem group of Brassicaceae and driven by dosage subfunctionalization. Our findings have profound implications for gene annotations in genomics, interpreting patterns of gene evolution and using genes in phylogeny reconstructions of species.

B and E MADS-box genes determine the perianth formation in Cymbidium goeringii Rchb.f

Cymbidium goeringii Rchb.f. is an important ornamental plant with a striking well-differentiated lip. Its complex floral architecture presents an exciting opportunity to examine perianth development. In flowering plants, class A, B and E floral homeotic genes play key roles in the specification of perianth identity. In this study, we used a cDNA library of wild-type C. goeringii flower buds for transcriptome sequencing. Eighteen candidate class A, B and E genes (including AP1/FUL-, AP2-, DEF-, GLO-, SEP- and AGL6-like genes) were identified. Quantitative real time polymerase chain reaction (qRT-PCR) results showed that CgDEF1, CgSEP2 and CgAGL6-1 were strongly detected only in the sepals and petals and were significantly downregulated in the lips. CgDEF3, CgDEF4 and CgAGL6-3 were highly expressed in the lips and lip-like petals but were only minimally detected in the sepals. Yeast two-hybrid analysis indicated that CgDEF1 and CgGLO formed a heterodimer. CgAGL6-1/CgSEP2 and CgDEF1 formed higher-order protein complexes with the assistance of the CgGLO protein, and both CgAGL6-1 and CgSEP2 formed a heterodimer. CgDEF3/CgDEF4 could interact independently with CgGLO and CgAGL6-3, respectively, while CgDEF3 and CgDEF4 also formed heterodimers with the assistance of the CgGLO. Based on a comprehensive analysis relating these gene expression patterns to protein interaction profiles, the mechanism of sepal/petal/lip determination was studied in C. goeringii. Furthermore, a hypothesis explaining the sepal/petal/lip determination of C. goeringii is proposed. The lip-quartet (CgDEF3/CgDEF4/CgAGL6-3/CgGLO) promoted lip formation, whereas the sepal/petal-quartet (CgDEF1/CgAGL6-1/CgSEP2/CgGLO) promoted sepal/petal formation. These results enrich the current knowledge regarding the mechanism and pathways of perianth formation in orchids.

Transcriptome-Wide Identification and Expression Analysis of DIVARICATA- and RADIALIS-Like Genes of the Mediterranean Orchid Orchis italica

Bilateral symmetry of flowers is a relevant novelty that has occurred many times throughout the evolution of flowering plants. In Antirrhinum majus, establishment of flower dorso-ventral asymmetry is mainly due to interaction of TCP (CYC and DICH) and MYB (DIV, RAD, and DRIF) transcription factors. In the present study, we characterized 8 DIV-, 4 RAD-, and 2 DRIF-like genes from the transcriptome of Orchis italica, an orchid species with bilaterally symmetric and resupinate flowers. We found a similar number of DIV- and RAD-like genes within the genomes of Phalaenopsis equestris and Dendrobium catenatum orchids. Orchid DIV- and RAD-like proteins share conserved motifs whose distribution reflects their phylogeny and analysis of the genomic organization revealed a single intron containing many traces of transposable elements. Evolutionary analysis has shown that purifying selection acts on the DIV- and RAD-like coding regions in orchids, with relaxation of selective constraints in a branch of the DIV-like genes. Analysis of the expression patterns of DIV- and RAD-like genes in O. italica revealed possible redundant functions for some of them. In the perianth of O. italica, the ortholog of DIV and DRIF of A. majus are expressed in all tissues, whereas RAD is mainly expressed in the outer tepals and lip. These data allow for proposal of an evolutionary conserved model in which the expression of the orthologs of the DIV, RAD, and DRIF genes might be related to establishment of flower bilateral symmetry in the nonmodel orchid species O. italica.

Exploring the evolutionary origin of floral organs of Erycina pusilla, an emerging orchid model system

BACKGROUND

Thousands of flowering plant species attract pollinators without offering rewards, but the evolution of this deceit is poorly understood. Rewardless flowers of the orchid Erycina pusilla have an enlarged median sepal and incised median petal ('lip') to attract oil-collecting bees. These bees also forage on similar looking but rewarding Malpighiaceae flowers that have five unequally sized petals and gland-carrying sepals. The lip of E. pusilla has a 'callus' that, together with winged 'stelidia', mimics these glands. Different hypotheses exist about the evolutionary origin of the median sepal, callus and stelidia of orchid flowers.

RESULTS

The evolutionary origin of these organs was investigated using a combination of morphological, molecular and phylogenetic techniques to a developmental series of floral buds of E. pusilla. The vascular bundle of the median sepal indicates it is a first whorl organ but its convex epidermal cells reflect convergence of petaloid features. Expression of AGL6 EpMADS4 and APETALA3 EpMADS14 is low in the median sepal, possibly correlating with its petaloid appearance. A vascular bundle indicating second whorl derivation leads to the lip. AGL6 EpMADS5 and APETALA3 EpMADS13 are most highly expressed in lip and callus, consistent with current models for lip identity. Six vascular bundles, indicating a stamen-derived origin, lead to the callus, stelidia and stamen. AGAMOUS is not expressed in the callus, consistent with its sterilization. Out of three copies of AGAMOUS and four copies of SEPALLATA, EpMADS22 and EpMADS6 are most highly expressed in the stamen. Another copy of AGAMOUS, EpMADS20, and the single copy of SEEDSTICK, EpMADS23, are most highly expressed in the stelidia, suggesting EpMADS22 may be required for fertile stamens.

CONCLUSIONS

The median sepal, callus and stelidia of E. pusilla appear to be derived from a sepal, a stamen that gained petal identity, and stamens, respectively. Duplications, diversifying selection and changes in spatial expression of different MADS-box genes shaped these organs, enabling the rewardless flowers of E. pusilla to mimic an unrelated rewarding flower for pollinator attraction. These genetic changes are not incorporated in current models and urge for a rethinking of the evolution of deceptive flowers.

The Evolution of Haploid Chromosome Numbers in the Sunflower Family

Chromosome number changes during the evolution of angiosperms are likely to have played a major role in speciation. Their study is of utmost importance, especially now, as a probabilistic model is available to study chromosome evolution within a phylogenetic framework. In the present study, likelihood models of chromosome number evolution were fitted to the largest family of flowering plants, the Asteraceae. Specifically, a phylogenetic supertree of this family was used to reconstruct the ancestral chromosome number and infer genomic events. Our approach inferred that the ancestral chromosome number of the family is n = 9. Also, according to the model that best explained our data, the evolution of haploid chromosome numbers in Asteraceae was a very dynamic process, with genome duplications and descending dysploidy being the most frequent genomic events in the evolution of this family. This model inferred more than one hundred whole genome duplication events; however, it did not find evidence for a paleopolyploidization at the base of this family, which has previously been hypothesized on the basis of sequence data from a limited number of species. The obtained results and potential causes of these discrepancies are discussed.

Evolutionary Dynamics of Floral Homeotic Transcription Factor Protein-Protein Interactions

Protein–protein interactions (PPIs) have widely acknowledged roles in the regulation of development, but few studies have addressed the timing and mechanism of shifting PPIs over evolutionary history. The B-class MADS-box transcription factors, PISTILLATA (PI) and APETALA3 (AP3) are key regulators of floral development. PI-like (PI^L) and AP3-like (AP3^L) proteins from a number of plants, including Arabidopsis thaliana (Arabidopsis) and the grass Zea mays (maize), bind DNA as obligate heterodimers. However, a PI^L protein from the grass relative Joinvillea can bind DNA as a homodimer. To ascertain whether Joinvillea PI^L homodimerization is an anomaly or indicative of broader trends, we characterized PI^L dimerization across the Poales and uncovered unexpected evolutionary lability. Both obligate B-class heterodimerization and PI^L homodimerization have evolved multiple times in the order, by distinct molecular mechanisms. For example, obligate B-class heterodimerization in maize evolved very recently from PI^L homodimerization. A single amino acid change, fixed during domestication, is sufficient to toggle one maize PI^L protein between homodimerization and obligate heterodimerization. We detected a signature of positive selection acting on residues preferentially clustered in predicted sites of contact between MADS-box monomers and dimers, and in motifs that mediate MADS PPI specificity in Arabidopsis. Changing one positively selected residue can alter PI^L dimerization activity. Furthermore, ectopic expression of a Joinvillea PI^L homodimer in Arabidopsis can homeotically transform sepals into petals. Our results provide a window into the evolutionary remodeling of PPIs, and show that novel interactions have the potential to alter plant form in a context-dependent manner.

Non-canonical structure, function and phylogeny of the Bsister MADS-box gene OsMADS30 of rice (Oryza sativa)

Bsister MADS-box genes play key roles in female reproductive organ and seed development throughout seed plants. This view is supported by their high conservation in terms of sequence, expression and function. In grasses, there are three subclades of Bsister genes: the OsMADS29-, the OsMADS30- and the OsMADS31-like genes. Here, we report on the evolution of the OsMADS30-like genes. Our analyses indicate that these genes evolved under relaxed purifying selection and are rather weakly expressed. OsMADS30, the representative of the OsMADS30-like genes from rice (Oryza sativa), shows strong sequence deviations in its 3' region when compared to orthologues from other grass species. We show that this is due to a 2.4-kbp insertion, possibly of a hitherto unknown helitron, which confers a heterologous C-terminal domain to OsMADS30. This putative helitron is not present in the OsMADS30 orthologues from closely related wild rice species, pointing to a relatively recent insertion event. Unlike other Bsister mutants O. sativa plants carrying a T-DNA insertion in the OsMADS30 gene do not show aberrant seed phenotypes, indicating that OsMADS30 likely does not have a canonical 'Bsister function'. However, imaging-based phenotyping of the T-DNA carrying plants revealed alterations in shoot size and architecture. We hypothesize that sequence deviations that accumulated during a period of relaxed selection in the gene lineage that led to OsMADS30 and the alteration of the C-terminal domain might have been a precondition for a potential neo-functionalization of OsMADS30 in O. sativa.

Digital Gene Expression Analysis Based on De Novo Transcriptome Assembly Reveals New Genes Associated with Floral Organ Differentiation of the Orchid Plant Cymbidium ensifolium

Cymbidium ensifolium belongs to the genus Cymbidium of the orchid family. Owing to its spectacular flower morphology, C. ensifolium has considerable ecological and cultural value. However, limited genetic data is available for this non-model plant, and the molecular mechanism underlying floral organ identity is still poorly understood. In this study, we characterize the floral transcriptome of C. ensifolium and present, for the first time, extensive sequence and transcript abundance data of individual floral organs. After sequencing, over 10 Gb clean sequence data were generated and assembled into 111,892 unigenes with an average length of 932.03 base pairs, including 1,227 clusters and 110,665 singletons. Assembled sequences were annotated with gene descriptions, gene ontology, clusters of orthologous group terms, the Kyoto Encyclopedia of Genes and Genomes, and the plant transcription factor database. From these annotations, 131 flowering-associated unigenes, 61 CONSTANS-LIKE (COL) unigenes and 90 floral homeotic genes were identified. In addition, four digital gene expression libraries were constructed for the sepal, petal, labellum and gynostemium, and 1,058 genes corresponding to individual floral organ development were identified. Among them, eight MADS-box genes were further investigated by full-length cDNA sequence analysis and expression validation, which revealed two APETALA1/AGL9-like MADS-box genes preferentially expressed in the sepal and petal, two AGAMOUS-like genes particularly restricted to the gynostemium, and four DEF-like genes distinctively expressed in different floral organs. The spatial expression of these genes varied distinctly in different floral mutant corresponding to different floral morphogenesis, which validated the specialized roles of them in floral patterning and further supported the effectiveness of our in silico analysis. This dataset generated in our study provides new insights into the molecular mechanisms underlying floral patterning of Cymbidium and supports a valuable resource for molecular breeding of the orchid plant.

The C-Terminal Sequence and PI motif of the Orchid (Oncidium Gower Ramsey) PISTILLATA (PI) Ortholog Determine its Ability to Bind AP3 Orthologs and Enter the Nucleus to Regulate Downstream Genes Controlling Petal and Stamen Formation

This study focused on the investigation of the effects of the PI motif and C-terminus of the Oncidium Gower Ramsey MADS box gene 8 (OMADS8), a PISTILLATA (PI) ortholog, on floral organ formation. 35S::OMADS8 completely rescued and 35S::OMADS8-PI (with the PI motif deleted) partially rescued petal/stamen formation, whereas these deficiencies were not rescued by 35S::OMADS8-C (C-terminal 29 amino acids deleted) in pi-1 mutants. OMADS8 could interact with Arabidopsis APETALA3 (AP3) and enter the nucleus. The nuclear entry efficiency was reduced for OMADS8-PI/AP3 and OMADS8-C/AP3. OMADS8 could also interact with OMADS5/OMADS9 (the Oncidium AP3 ortholog) and enter the nucleus with an efficiency only slightly affected by the deletion of the C-terminal sequence or PI motif. However, the stability of the OMADS8/OMADS5 and OMADS8/OMADS9 complexes was significantly reduced by deletion of the C-terminal sequence or PI motif. Further analysis indicated that the expression of genes downstream of AP3/PI (BNQ1/BNQ2/GNC/At4g30270) was compensated by 35S::OMADS8 and 35S::OMADS8-PI to a level similar to wild-type plants but was not affected by 35S::OMADS8-C in the pi-1 mutants. A similar FRET (fluorescence resonance energy transfer) efficiency was observed for Arabidopsis AGAMOUS (AG) and the Oncidium AG ortholog OMADS4 for OMADS8, OMADS8-PI and OMADS8-C. These results indicated that the OMADS8 PI motif and C-terminus were valuable for the interaction of OMADS8 with the AP3 orthologs to form higher order heterotetrameric complexes that regulated petal/stamen formation in both Oncidium orchids and transgenic Arabidopsis. However, the C-terminal sequence and PI motif were dispensable for the interaction of OMADS8 with the AG orthologs.

Role of transcriptional regulation in the evolution of plant phenotype: A dynamic systems approach

A growing body of evidence suggests that alterations in transcriptional regulation of genes involved in modulating development are an important part of phenotypic evolution, and this can be documented among species and within populations. While the effects of differential transcriptional regulation in organismal development have been preferentially studied in animal systems, this phenomenon has also been addressed in plants. In this review, we summarize evidence for cis-regulatory mutations, trans-regulatory changes and epigenetic modifications as molecular events underlying important phenotypic alterations, and thus shaping the evolution of plant development. We postulate that a mechanistic understanding of why such molecular alterations have a key role in development, morphology and evolution will have to rely on dynamic models of complex regulatory networks that consider the concerted action of genetic and nongenetic components, and that also incorporate the restrictions underlying the genotype to phenotype mapping process. Developmental Dynamics 244:1074-1095, 2015. © 2015 Wiley Periodicals, Inc.

Histone acetylation accompanied with promoter sequences displaying differential expression profiles of B-class MADS-box genes for phalaenopsis floral morphogenesis

Five B-class MADS-box genes, including four APETALA3 (AP3)-like PeMADS2∼5 and one PISTILLATA (PI)-like PeMADS6, specify the spectacular flower morphology in orchids. The PI-like PeMADS6 ubiquitously expresses in all floral organs. The four AP3-like genes, resulted from two duplication events, express ubiquitously at floral primordia and early floral organ stages, but show distinct expression profiles at late floral organ primordia and floral bud stages. Here, we isolated the upstream sequences of PeMADS2∼6 and studied the regulatory mechanism for their distinct gene expression. Phylogenetic footprinting analysis of the 1.3-kb upstream sequences of AP3-like PeMADS2∼5 showed that their promoter regions have sufficiently diverged and contributed to their subfunctionalization. The amplified promoter sequences of PeMADS2∼6 could drive beta-glucuronidase (GUS) gene expression in all floral organs, similar to their expression at the floral primordia stage. The promoter sequence of PeMADS4, exclusively expressed in lip and column, showed a 1.6∼3-fold higher expression in lip/column than in sepal/petal. Furthermore, we noted a 4.9-fold increase in histone acetylation (H3K9K14ac) in the translation start region of PeMADS4 in lip as compared in petal. All these results suggest that the regulation via the upstream sequences and increased H3K9K14ac level may act synergistically to display distinct expression profiles of the AP3-like genes at late floral organ primordia stage for Phalaenopsis floral morphogenesis.

DEF- and GLO-like proteins may have lost most of their interaction partners during angiosperm evolution

BACKGROUND AND AIMS

DEFICIENS (DEF)- and GLOBOSA (GLO)-like proteins constitute two sister clades of floral homeotic transcription factors that were already present in the most recent common ancestor (MRCA) of extant angiosperms. Together they specify the identity of petals and stamens in flowering plants. In core eudicots, DEF- and GLO-like proteins are functional in the cell only as heterodimers with each other. There is evidence that this obligate heterodimerization contributed to the canalization of the flower structure of core eudicots during evolution. It remains unknown as to whether this strict heterodimerization is an ancient feature that can be traced back to the MRCA of extant flowering plants or if it evolved later during the evolution of the crown group angiosperms.

METHODS

The interactions of DEF- and GLO-like proteins of the early-diverging angiosperms Amborella trichopoda and Nuphar advena and of the magnoliid Liriodendron tulipifera were analysed by employing yeast two-hybrid analysis and electrophoretic mobility shift assay (EMSA). Character-state reconstruction, including data from other species as well, was used to infer the ancestral interaction patterns of DEF- and GLO-like proteins.

KEY RESULTS

The yeast two-hybrid and EMSA data suggest that DEF- and GLO-like proteins from early-diverging angiosperms both homo- and heterodimerize. Character-state reconstruction suggests that the ability to form heterodimeric complexes already existed in the MRCA of extant angiosperms and that this property remained highly conserved throughout angiosperm evolution. Homodimerization of DEF- and GLO-like proteins also existed in the MRCA of all extant angiosperms. DEF-like protein homodimerization was probably lost very early in angiosperm evolution and was not present in the MRCA of eudicots and monocots. GLO-like protein homodimerization might have been lost later during evolution, but very probably was not present in the MRCA of eudicots.

CONCLUSIONS

The flexibility of DEF- and GLO-like protein interactions in early-diverging angiosperms may be one reason for the highly diverse flower morphologies observed in these species. The results strengthen the hypothesis that a reduction in the number of interaction partners of DEF- and GLO-like proteins, with DEF-GLO heterodimers remaining the only DNA-binding dimers in core eudicots, contributed to developmental robustness, canalization of flower development and the diversification of angiosperms.

Co-modification of class B genes TfDEF and TfGLO in Torenia fournieri Lind. alters both flower morphology and inflorescence architecture

The class B genes DEFICIENS (DEF)/APETALA3 (AP3) and GLOBOSA (GLO)/PISTILLATA (PI), encoding MADS-box transcription factors, and their functions in petal and stamen development have been intensely studied in Arabidopsis and Antirrhinum. However, the functions of class B genes in other plants, including ornamental species exhibiting floral morphology different from these model plants, have not received nearly as much attention. Here, we examine the cooperative functions of TfDEF and TfGLO on floral organ development in the ornamental plant torenia (Torenia fournieri Lind.). Torenia plants co-overexpressing TfDEF and TfGLO showed a morphological alteration of sepals to petaloid organs. Phenotypically, these petaloid sepals were nearly identical to petals but had no stamens or yellow patches like those of wild-type petals. Furthermore, the inflorescence architecture in the co-overexpressing torenias showed a characteristic change in which, unlike the wild-types, their flowers developed without peduncles. Evaluation of the petaloid sepals showed that these attained a petal-like nature in terms of floral organ phenotype, cell shape, pigment composition, and the expression patterns of anthocyanin biosynthesis-related genes. In contrast, torenias in which TfDEF and TfGLO were co-suppressed exhibited sepaloid petals in the second whorl. The sepaloid petals also attained a sepal-like nature, in the same way as the petaloid sepals. The results clearly demonstrate that TfDEF and TfGLO play important cooperative roles in petal development in torenia. Furthermore, the unique transgenic phenotypes produced create a valuable new way through which characteristics of petal development and inflorescence architecture can be investigated in torenia.

De novo transcriptome assembly from inflorescence of Orchis italica: analysis of coding and non-coding transcripts

The floral transcriptome of Orchis italica, a wild orchid species, was obtained using Illumina RNA-seq technology and specific de novo assembly and analysis tools. More than 100 million raw reads were processed resulting in 132,565 assembled transcripts and 86,079 unigenes with an average length of 606 bp and N50 of 956 bp. Functional annotation assigned 38,984 of the unigenes to records present in the NCBI non-redundant protein database, 32,161 of them to Gene Ontology terms, 15,775 of them to Eukaryotic Orthologous Groups (KOG) and 7,143 of them to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The in silico expression analysis based on the Fragments Per Kilobase of transcript per Million mapped reads (FPKM) was confirmed by real-time RT-PCR experiments on 10 selected unigenes, which showed high and statistically significant positive correlation with the RNA-seq based expression data. The prediction of putative long non-coding RNAs was assessed using two different software packages, CPC and Portrait, resulting in 7,779 unannotated unigenes that matched the threshold values for both of the analyses. Among the predicted long non-coding RNAs, one is the homologue of TAS3, a long non-coding RNA precursor of trans-acting small interfering RNAs (ta-siRNAs). The differential expression pattern observed for the selected putative long non-coding RNAs suggests their possible functional role in different floral tissues.

Expression of paralogous SEP-, FUL-, AG- and STK-like MADS-box genes in wild-type and peloric Phalaenopsis flowers

The diverse flowers of Orchidaceae are the result of several major morphological transitions, among them the most studied is the differentiation of the inner median tepal into the labellum, a perianth organ key in pollinator attraction. Type A peloria lacking stamens and with ectopic labella in place of inner lateral tepals are useful for testing models on the genes specifying these organs by comparing their patterns of expression between wild-type and peloric flowers. Previous studies focused on DEFICIENS- and GLOBOSA-like MADS-box genes because of their conserved role in perianth and stamen development. The "orchid code" model summarizes this work and shows in Orchidaceae there are four paralogous lineages of DEFICIENS/AP3-like genes differentially expressed in each floral whorl. Experimental tests of this model showed the conserved, higher expression of genes from two specific DEF-like gene lineages is associated with labellum development. The present study tests whether eight MADS-box candidate SEP-, FUL-, AG-, and STK-like genes have been specifically duplicated in the Orchidaceae and are also differentially expressed in association with the distinct flower organs of Phalaenopsis hyb. "Athens." The gene trees indicate orchid-specific duplications. In a way analogous to what is observed in labellum-specific DEF-like genes, a two-fold increase in the expression of SEP3-like gene PhaMADS7 was measured in the labellum-like inner lateral tepals of peloric flowers. The overlap between SEP3-like and DEF-like genes suggests both are associated with labellum specification and similar positional cues determine their domains of expression. In contrast, the uniform messenger levels of FUL-like genes suggest they are involved in the development of all organs and their expression in the ovary suggests cell differentiation starts before pollination. As previously reported AG-like and STK-like genes are exclusively expressed in gynostemium and ovary, however no evidence for transcriptional divergence was found in the stage investigated. Gene expression suggests a developmental regulatory system based on the combined activity of duplicate MADS-box genes. We discuss its feasibility based on documented protein interactions and patterns of expression.

Perspectives on MADS-box expression during orchid flower evolution and development

The diverse morphology of orchid flowers and their complex, often deceptive strategies to become pollinated have fascinated researchers for a long time. However, it was not until the 20th century that the ontogeny of orchid flowers, the genetic basis of their morphology and the complex phylogeny of Orchidaceae were investigated. In parallel, the improvement of techniques for in vitro seed germination and tissue culture, together with studies on biochemistry, physiology, and cytology supported the progress of what is now a highly productive industry of orchid breeding and propagation. In the present century both basic research in orchid flower evo-devo and the interest for generating novel horticultural varieties have driven the characterization of many members of the MADS-box family encoding key regulators of flower development. This perspective summarizes the picture emerging from these studies and discusses the advantages and limitations of the comparative strategy employed so far. I address the growing role of natural and horticultural mutants in these studies and the emergence of several model species in orchid evo-devo and genomics. In this context, I make a plea for an increasingly integrative approach.

bZIP transcription factors in the oomycete phytophthora infestans with novel DNA-binding domains are involved in defense against oxidative stress

Transcription factors of the basic leucine zipper (bZIP) family control development and stress responses in eukaryotes. To date, only one bZIP has been described in any oomycete; oomycetes are members of the stramenopile kingdom. In this study, we describe the identification of 38 bZIPs from the Phytophthora infestans genome. Half contain novel substitutions in the DNA-binding domain at a site that in other eukaryotes is reported to always be Asn. Interspecific comparisons indicated that the novel substitutions (usually Cys, but also Val and Tyr) arose after oomycetes diverged from other stramenopiles. About two-thirds of P. infestans bZIPs show dynamic changes in mRNA levels during the life cycle, with many of the genes being upregulated in sporangia, zoospores, or germinated zoospore cysts. One bZIP with the novel Cys substitution was shown to reside in the nucleus throughout growth and development. Using stable gene silencing, the functions of eight bZIPs with the Cys substitution were tested. All but one were found to play roles in protecting P. infestans from hydrogen peroxide-induced injury, and it is proposed that the novel Cys substitution serves as a redox sensor. A ninth bZIP lacking the novel Asn-to-Cys substitution, but having Cys nearby, was also shown through silencing to contribute to defense against peroxide. Little effect on asexual development, plant pathogenesis, or resistance to osmotic stress was observed in transformants silenced for any of the nine bZIPs.

Organ homologies in orchid flowers re-interpreted using the Musk Orchid as a model

Background and Aims. The presence of novel structures in orchid flowers, including auricles, rostellum and bursicles on the gynostemium and a lobed labellum, has prompted long-standing homology disputes, fuelled by conflicting evidence from a wide range of sources. Re-assessment of this debate using an improved model is timely, following recent phylogenetic insights and on the cusp of a revolution in developmental genetics.

Methods. We use new data from floral development and anatomy in the small-flowered terrestrial orchid Herminium monorchis as a model to explore organ homologies in orchid flowers within the context of a review of recent literature on developmental genetics.

Key Results. The apex of the median carpel of Herminium is trilobed, and the bursicles develop from its lateral lobes, relatively late in flower ontogeny. The bursicles enclose the viscidia, which adhere to the tapetal remnants to form a caudicle linking the viscidium with the pollinium. The auricles are initiated earlier than the bursicles, but they also remain unvascularized. The deeply trilobed labellum possesses three vascular traces, in contrast with the lateral petals, each of which contains a single vascular trace. The two lateral labellum traces diverge from the traces supplying the two adjacent lateral sepals. Data from flower ontogeny and anatomy conflict with respect to organ homologies.

Conclusions. Much progress has recently been made in understanding the exceptional differentiation shown by orchids among perianth segments, focusing on multiple copies of the DEF/AP3 subclass of B-class MADS-box genes. In contrast, untangling homologies of profound congenital union of multiple floral organs forming the orchid gynostemium is hampered by their profound congenital union, which we ascribe to overlap in gene expression between organs. Thus, the functional morphology of the orchid flower could ultimately reflect extreme synorganization and associated genetic integration. Analogizing the deeply lobed orchid labellum with a compound leaf, we speculate that KNOX genes could be implicated not only in their demonstrated role in spur development but also in the development of both the characteristic lobed morphology of the orchid labellum and the lobing of the median carpel that differentiates the bursicles and rostellum.

The MADS and the Beauty: Genes Involved in the Development of Orchid Flowers

Since the time of Darwin, biologists have studied the origin and evolution of the Orchidaceae, one of the largest families of flowering plants. In the last two decades, the extreme diversity and specialization of floral morphology and the uncoupled rate of morphological and molecular evolution that have been observed in some orchid species have spurred interest in the study of the genes involved in flower development in this plant family. As part of the complex network of regulatory genes driving the formation of flower organs, the MADS-box represents the most studied gene family, both from functional and evolutionary perspectives. Despite the absence of a published genome for orchids, comparative genetic analyses are clarifying the functional role and the evolutionary pattern of the MADS-box genes in orchids. Various evolutionary forces act on the MADS-box genes in orchids, such as diffuse purifying selection and the relaxation of selective constraints, which sometimes reveals a heterogeneous selective pattern of the coding and non-coding regions. The emerging theory regarding the evolution of floral diversity in orchids proposes that the diversification of the orchid perianth was a consequence of duplication events and changes in the regulatory regions of the MADS-box genes, followed by sub- and neo-functionalization. This specific developmental-genetic code is termed the "orchid code."

When ABC becomes ACB

Understanding how the information contained in genes is mapped onto the phenotypes, and deriving formal frameworks to search for generic aspects of developmental constraints and evolution remains one of the main challenges of contemporary biological research. The Mexican endemic triurid Lacandonia schismatica (Lacandoniaceae), a mycoheterotrophic monocotyledonous plant with hermaphroditic reproductive axes is alone among 250,000 species of angiosperms, as it has central stamens surrounded by a peripheral gynoecium, representing a natural instance of a homeotic mutant. Based on the classical ABC model of flower development, it has recently been shown that the B-function gene APETALA3 (AP3), essential for stamen identity, was displaced toward the flower centre in L. schismatica (ABC to ACB) from the early stages of flower development. A functional conservation of B-function genes from L. schismatica through the rescue of B-gene mutants in Arabidopsis thaliana, as well as conserved protein interactions, has also been demonstrated. Thus, it has been shown that relatively simple genetic alterations may underlie large morphological shifts fixed in extant natural populations. Nevertheless, critical questions remain in order to have a full and sufficient explanation of the molecular genetic mechanisms underlying L. schismatica's unique floral arrangement. Evolutionary approaches to developmental mechanisms and systems biology, including high-throughput functional genomic studies and models of complex developmental gene regulatory networks, constitute two main approaches to meet such a challenge. In this review, the aim is to address some of the pending questions with the ultimate goal of investigating further the mechanisms of L. schismatica's unique homeotic flower arrangement and its evolution.

Parallel evolution of TCP and B-class genes in Commelinaceae flower bilateral symmetry

BACKGROUND

Flower bilateral symmetry (zygomorphy) has evolved multiple times independently across angiosperms and is correlated with increased pollinator specialization and speciation rates. Functional and expression analyses in distantly related core eudicots and monocots implicate independent recruitment of class II TCP genes in the evolution of flower bilateral symmetry. Furthermore, available evidence suggests that monocot flower bilateral symmetry might also have evolved through changes in B-class homeotic MADS-box gene function.

METHODS

In order to test the non-exclusive hypotheses that changes in TCP and B-class gene developmental function underlie flower symmetry evolution in the monocot family Commelinaceae, we compared expression patterns of teosinte branched1 (TB1)-like, DEFICIENS (DEF)-like, and GLOBOSA (GLO)-like genes in morphologically distinct bilaterally symmetrical flowers of Commelina communis and Commelina dianthifolia, and radially symmetrical flowers of Tradescantia pallida.

RESULTS

Expression data demonstrate that TB1-like genes are asymmetrically expressed in tepals of bilaterally symmetrical Commelina, but not radially symmetrical Tradescantia, flowers. Furthermore, DEF-like genes are expressed in showy inner tepals, staminodes and stamens of all three species, but not in the distinct outer tepal-like ventral inner tepals of C. communis.

CONCLUSIONS

Together with other studies, these data suggest parallel recruitment of TB1-like genes in the independent evolution of flower bilateral symmetry at early stages of Commelina flower development, and the later stage homeotic transformation of C. communis inner tepals into outer tepals through the loss of DEF-like gene expression.

Functional analysis reveals the possible role of the C-terminal sequences and PI motif in the function of lily (Lilium longiflorum) PISTILLATA (PI) orthologues

Two lily (Lilium longiflorum) PISTILLATA (PI) genes, Lily MADS Box Gene 8 and 9 (LMADS8/9), were characterized. LMADS9 lacked 29 C-terminal amino acids including the PI motif that was present in LMADS8. Both LMADS8/9 mRNAs were prevalent in the first and second whorl tepals during all stages of development and were expressed in the stamen only in young flower buds. LMADS8/9 could both form homodimers, but the ability of LMADS8 homodimers to bind to CArG1 was relatively stronger than that of LMADS9 homodimers. 35S:LMADS8 completely, and 35S:LMADS9 only partially, rescued the second whorl petal formation and partially converted the first whorl sepal into a petal-like structure in Arabidopsis pi-1 mutants. Ectopic expression of LMADS8-C (with deletion of the 29 amino acids of the C-terminal sequence) or LMADS8-PI (with only the PI motif deleted) only partially rescued petal formation in pi mutants, which was similar to what was observed in 35S:LMADS9/pi plants. In contrast, 35:LMADS9+L8C (with the addition of the 29 amino acids of the LMADS8 C-terminal sequence) or 35S:LMADS9+L8PI (with the addition of the LMADS8 PI motif) demonstrated an increased ability to rescue petal formation in pi mutants, which was similar to what was observed in 35S:LMADS8/pi plants. Furthermore, ectopic expression of LMADS8-M (with the MADS domain truncated) generated more severe dominant negative phenotypes than those seen in 35S:LMADS9-M flowers. These results revealed that the 29 amino acids including the PI motif in the C-terminal region of the lily PI orthologue are valuable for its function in regulating perianth organ formation.

Agave tequilana MADS genes show novel expression patterns in meristems, developing bulbils and floral organs

Agave tequilana is a monocarpic perennial species that flowers after 5-8 years of vegetative growth signaling the end of the plant's life cycle. When fertilization is unsuccessful, vegetative bulbils are induced on the umbels of the inflorescence near the bracteoles from newly formed meristems. Although the regulation of inflorescence and flower development has been described in detail for monocarpic annuals and polycarpic species, little is known at the molecular level for these processes in monocarpic perennials, and few studies have been carried out on bulbils. Histological samples revealed the early induction of umbel meristems soon after the initiation of the vegetative to inflorescence transition in A. tequilana. To identify candidate genes involved in the regulation of floral induction, a search for MADS-box transcription factor ESTs was conducted using an A. tequilana transcriptome database. Seven different MIKC MADS genes classified into 6 different types were identified based on previously characterized A. thaliana and O. sativa MADS genes and sequences from non-grass monocotyledons. Quantitative real-time PCR analysis of the seven candidate MADS genes in vegetative, inflorescence, bulbil and floral tissues uncovered novel patterns of expression for some of the genes in comparison with orthologous genes characterized in other species. In situ hybridization studies using two different genes showed expression in specific tissues of vegetative meristems and floral buds. Distinct MADS gene regulatory patterns in A. tequilana may be related to the specific reproductive strategies employed by this species.

Expression pattern of two paralogs of the PI/GLO-like locus during Orchis italica (Orchidaceae, Orchidinae) flower development

The class B MADS-box genes belong to two distinct functional groups: the AP3/DEF-like and the PI/GLO-like sub-families. In orchids, AP3/DEF-like genes are present in four copies, each with a different role in floral organ formation, which is described in the "orchid code" model. Interestingly, the orchid PI/GLO-like genes are present in two copies in Orchidinae, whereas they are described as single copy in the other orchid lineages. The two PI/GLO-like paralogs have site-specific different selective constraints; in addition, they show relaxation of purifying selection when compared to the single-copy lineages. In this study, we present a comparative analysis of the expression patterns of the two PI/GLO-like paralogs, OrcPI and OrcPI2, in floral tissues of Orchis italica in different developmental stages using real-time PCR. The two genes show similar expression profiles in the tissue examined, with differences detectable between immature and mature inflorescence. In all cases, OrcPI2 is expressed at a higher level than OrcPI. Real-time PCR results reveal that the co-expression of the two duplicated loci could have a fully or partially redundant function. The possible evolutionary fate of OrcPI and OrcPI2 is discussed as well as their involvement in ovary development.

The duplicated B-class MADS-box genes display dualistic characters in orchid floral organ identity and growth

Orchidaceae are an excellent model to examine perianth development because of their sophisticated floral architecture. In this study, we identified 24 APETALA3 (AP3)-like and 13 PISTILLA (PI)-like genes from 11 species of orchids and characterized them into four AP3- and two PI-duplicated homologs. The first duplication event in AP3 homologs occurring in the early evolutionary history of the Orchidaceae gave rise to AP3A and AP3B clades. Further duplication events resulted in four subclades, namely AP3A1, AP3A2, AP3B1 and AP3B2, during the evolution of Orchidaceae. The AP3 paralogous genes were expressed throughout inflorescence and floral bud development. From the in situ hybridization results, we noticed that the transition timings from ubiquitous to constrained expression in floral organs for both clades are different. The transition point of expression of the AP3A clade (clades 3 and 4) was at the late floral organ primordia stage. In contrast, that for the AP3B clade (clades 1 and 2) was not observed until the late inflorescence and floral bud stages. In addition, the AP3 orthologous genes revealed diverse expression patterns in various species of orchids, whereas the PI homologs were uniformly expressed in all floral whorls. AP3A2 orthologs play a noticeable role in lip formation because of their exclusive expression in the lip. Further evidence comes from the ectopic expression of AP3A2 detected in the lip-like petals extending from the lip in four sets of peloric mutants. Finally, a Homeotic Orchid Tepal (HOT) model is proposed, in which dualistic characters of duplicated B-class MADS-box genes are involved in orchid perianth development and growth.

Conserved differential expression of paralogous DEFICIENS- and GLOBOSA-like MADS-box genes in the flowers of Orchidaceae: refining the 'orchid code'

In flowering plants, class-B floral homeotic genes encode MADS-domain transcription factors, which are key in the specification of petal and stamen identity, and have two ancient clades: DEF-like and GLO-like genes. Many species have one gene of each clade, but orchids have typically four DEF-like genes, representing ancient gene clades 1, 2, 3 and 4. We tested the 'orchid code', a combinatorial genetic model suggesting that differences between the organs of the orchid perianth (outer tepals, inner lateral tepals and labellum) are generated by the combinatorial differential expression of four DEF-like genes. Our experimental test involves highly sensitive and specific measurements, with qRT-PCR of the expression of DEF- and GLO-like genes from the distantly related Vanilla planifolia and Phragmipedium longifolium, as well as from wild-type and peloric Phalaenopsis hybrid flowers. Our findings support the first 'orchid code' hypothesis, in that absence of clade-3 and -4 gene expression distinguishes the outer tepals from the inner tepals. In contrast to the original hypothesis, however, mRNA of both clade-3 and -4 genes accumulates in wild-type inner lateral tepals and the labellum, and in labellum-like inner lateral tepals of peloric flowers, albeit in different quantities. Our data suggest a revised hypothesis where high levels of clade-1 and -2, and low levels of clade-3 and -4, gene expression specify inner lateral tepals, whereas labellum development requires low levels of clade-1 and -2 expression and high levels of clade-3 and -4 expression.

Robustness and evolvability in the B-system of flower development

BACKGROUND

Gene duplication has often been invoked as a key mechanism responsible for evolution of new morphologies. The floral homeotic B-group gene family has undergone a number of gene duplication events, and yet the functions of these genes appear to be largely conserved. However, detailed comparative analysis has indicated that such duplicate genes have considerable cryptic variability in their functions. In the Solanaceae, two duplicate B-group gene lineages have been retained in three subfamilies. Comparisons of orthologous genes across members of the Solanaceae have demonstrated that the combined function of all four B-gene members is to establish petal and stamen identity, but that this function was partitioned differently in each species. These observations emphasize both the robustness and the evolvability of the B-system.

SCOPE

We provide an overview of how the B-function genes can robustly specify petal and stamen identity and at the same time evolve through changes in protein-protein interaction, gene expression patterns, copy number variation or alterations in the downstream genes they control. By using mathematical models we explore regulatory differences between species and how these impose constraints on downstream gene regulation.

CONCLUSIONS

Evolvability of the B-genes can be understood through the multiple ways in which the B-system can be robust. Quantitative approaches should allow for the incorporation of more biological realism in the representations of these regulatory systems and this should contribute to understanding the constraints under which different B-systems can function and evolve. This, in turn, can provide a better understanding of the ways in which B-genes have contributed to flower diversity.

The PI/GLO-like locus in orchids: duplication and purifying selection at synonymous sites within Orchidinae (Orchidaceae)

Positive selection and relaxation of purifying constraints after duplication events have driven the functional diversification of gene families involved in development. One example of this occurred within the plant MADS-box genes. The evolution of the orchid flower was driven by duplication events followed by sub- and neo-functionalization of class B DEF-like MADS-box genes, which are present at three to four copies in the orchid genome. In contrast, the orchid PI/GLO-like class B MADS-box genes have been reported thus far as single-copy loci, with the only exception of Habenaria radiata. We isolated a novel PI/GLO-like gene (OrcPI2) in Orchis italica, which is different than the previously characterized OrcPI locus. The presence of two functional paralogs of PI/GLO-like genes in orchids is detectable only within the tribe Orchidinae. Evolutionary analyses revealed an apparent relaxation of purifying selection acting on the two PI/GLO-like paralogs of the Orchidinae when compared to the single-copy PI/GLO-like genes found in other orchid species. Furthermore, by measuring dN/dS (ω) ratios, we show that a high percentage of sites between the two PI/GLO-like paralogs have different evolutionary pressures. Interestingly, the apparent relaxation of selective constraints on the two PI/GLO-like paralogs is due to strong purifying selection at synonymous sites rather than to a high value of nonsynonymous substitution rate. This peculiar evolutionary pattern might be related to molecular processes such as mRNA folding and/or translational efficiency control. These processes could potentially be involved in or predate the functional diversification of the two PI/GLO-like paralogs within Orchidinae.

Evidence for the involvement of Globosa-like gene duplications and expression divergence in the evolution of floral morphology in the Zingiberales

*The MADS box transcription factor family has long been identified as an important contributor to the control of floral development. It is often hypothesized that the evolution of floral development across angiosperms and within specific lineages may occur as a result of duplication, functional diversification, and changes in regulation of MADS box genes. Here we examine the role of Globosa (GLO)-like genes, members of the B-class MADS box gene lineage, in the evolution of floral development within the monocot order Zingiberales. *We assessed changes in perianth and stamen whorl morphology in a phylogenetic framework. We identified GLO homologs (ZinGLO1-4) from 50 Zingiberales species and investigated the evolution of this gene lineage. Expression of two GLO homologs was assessed in Costus spicatus and Musa basjoo. *Based on the phylogenetic data and expression results, we propose several family-specific losses and gains of GLO homologs that appear to be associated with key morphological changes. The GLO-like gene lineage has diversified concomitant with the evolution of the dimorphic perianth and the staminodial labellum. *Duplications and expression divergence within the GLO-like gene lineage may have played a role in floral diversification in the Zingiberales.

Isolation and phylogenetic footprinting analysis of the 5'-regulatory region of the floral homeotic gene OrcPI from Orchis italica (Orchidaceae)

The nucleotide sequences of regulatory elements from homologous genes can be strongly divergent. Phylogenetic footprinting, a comparative analysis of noncoding regions, can detect putative transcription factor binding sites (TFBSs) shared among the regulatory regions of 2 or more homologous genes. These conserved motifs have the potential to serve the same regulatory function in distantly related taxa. We isolated the 5'-noncoding region of the OrcPI gene, a MADS-box transcription factor involved in flower development in Orchis italica, using the thermal asymmetric interlaced polymerase chain reaction technique. This region (comprising 1352 bp) induced transient beta-glucuronidase expression in the petal tissue of white Rosa hybrida flowers and represents the 5'-regulatory sequence of the OrcPI gene. Phylogenetic footprinting analysis detected conserved regions within the 5'-regulatory sequence of OrcPI and the homologous regions of Oryza sativa, Lilium regale, and Arabidopsis thaliana. Some of these sequences are known TFBSs described in databases of plant regulatory elements. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the following accession numbers: AF198055 promoter region of the PISTILLATA (PI) gene of A. thaliana; AB094985 cDNA of OrcPI (PI/GLOBOSA [PI/GLO] homologue) of O. italica; AB378089 5'-regulatory region of the OrcPI gene of O. italica; AP008211 putative promoter region of OSMADS2 (PI/GLO homologue) of O. sativa; AP008207 putative promoter region of OSMADS4 (PI/GLO homologue) of O. sativa; and AB158292 putative promoter region of the PI/GLO homologue of L. regale.

Characterization of the possible roles for B class MADS box genes in regulation of perianth formation in orchid

To investigate sepal/petal/lip formation in Oncidium Gower Ramsey, three paleoAPETALA3 genes, O. Gower Ramsey MADS box gene5 (OMADS5; clade 1), OMADS3 (clade 2), and OMADS9 (clade 3), and one PISTILLATA gene, OMADS8, were characterized. The OMADS8 and OMADS3 mRNAs were expressed in all four floral organs as well as in vegetative leaves. The OMADS9 mRNA was only strongly detected in petals and lips. The mRNA for OMADS5 was only strongly detected in sepals and petals and was significantly down-regulated in lip-like petals and lip-like sepals of peloric mutant flowers. This result revealed a possible negative role for OMADS5 in regulating lip formation. Yeast two-hybrid analysis indicated that OMADS5 formed homodimers and heterodimers with OMADS3 and OMADS9. OMADS8 only formed heterodimers with OMADS3, whereas OMADS3 and OMADS9 formed homodimers and heterodimers with each other. We proposed that sepal/petal/lip formation needs the presence of OMADS3/8 and/or OMADS9. The determination of the final organ identity for the sepal/petal/lip likely depended on the presence or absence of OMADS5. The presence of OMADS5 caused short sepal/petal formation. When OMADS5 was absent, cells could proliferate, resulting in the possible formation of large lips and the conversion of the sepal/petal into lips in peloric mutants. Further analysis indicated that only ectopic expression of OMADS8 but not OMADS5/9 caused the conversion of the sepal into an expanded petal-like structure in transgenic Arabidopsis (Arabidopsis thaliana) plants.