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

Identification and expression profiling of selected MADS-box family genes in Dendrobium officinale

Dendrobium officinale, a herb with highly medicinal and ornamental value, is widely distributed in China. MADS-box genes encode transcription factors that regulate various growth and developmental processes in plants, particular in flowering. However, the MADS-box genes in D. officinale are largely unknown. In our study, expression profiling analyses of selected MADS-box genes in D. officinale were performed. In total, 16 DnMADS-box genes with full-length ORF were identified and named according to their phylogenetic relationships with model plants. The transient expression of eight selected MADS-box genes in the epidermal cells of tobacco leaves showed that these DnMADS-box proteins localized to the nucleus. Tissue-specific expression analysis pointed out eight flower-specific expressed MADS-box genes in D. officinale. Furthermore, expression patterns of DnMADS-box genes were investigated during the floral transition process. DnMADS3, DnMADS8 and DnMADS22 were significantly up-regulated in the reproductive phase compared with the vegetative phase, suggesting putative roles of these DnMADS-box genes in flowering. Our data showed that the expressions of MADS-box genes in D. officinale were controlled by diverse exogenous phytohormones. Together, these findings will facilitate further studies of MADS-box genes in Orchids and broaden our understanding of the genetics of flowering.

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

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 MIKC^C 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 MIKC^C 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.

Evolutionary Conservation of the Orchid MYB Transcription Factors DIV, RAD, and DRIF

The MYB transcription factors DIVARICATA (DIV), DIV-and-RAD-Interacting-Factor (DRIF), and the small interfering peptide RADIALIS (RAD) can interact, forming a regulatory module that controls different plant developmental processes. In the snapdragon Antirrhinum majus, this module, together with the TCP transcription factor CYCLOIDEA (CYC), is responsible for the establishment of floral dorsoventral asymmetry. The spatial gene expression pattern of the OitDIV, OitDRIF, and OitRAD homologs of Orchis italica, an orchid with zygomorphic flowers, has suggested a possible conserved role of these genes in bilateral symmetry of the orchid flower. Here, we have identified four DRIF genes of orchids and have reconstructed their genomic organization and evolution. In addition, we found snapdragon transcriptional cis-regulatory elements of DIV and RAD loci generally conserved within the corresponding orchid orthologues. We have tested the biochemical interactions among OitDIV, OitDRIF1, and OitRAD of O. italica, showing that OitDRIF1 can interact both with OitDIV and OitRAD, whereas OitDIV and OitRAD do not directly interact, as in A. majus. The analysis of the quantitative expression profile of these MYB genes revealed that in zygomorphic orchid flowers, the DIV, DRIF1, and RAD transcripts are present at higher levels in the lip than in lateral inner tepals, whereas in peloric orchid flowers they show similar expression levels. These results indicate that MYB transcription factors could have a role in shaping zygomorphy of the orchid flower, potentially enriching the underlying orchid developmental code.

Floral Induction and Flower Development of Orchids

Orchids comprise one of the largest, most highly evolved angiosperm families, and form an extremely peculiar group of plants. Various orchids are available through traditional breeding and micro-propagation since they are valuable as potted plants and/or cut flowers in horticultural markets. The flowering of orchids is generally influenced by environmental signals such as temperature and endogenous developmental programs controlled by genetic factors as is usual in many flowering plant species. The process of floral transition is connected to the flower developmental programs that include floral meristem maintenance and floral organ specification. Thanks to advances in molecular and genetic technologies, the understanding of the molecular mechanisms underlying orchid floral transition and flower developmental processes have been widened, especially in several commercially important orchids such as Phalaenopsis, Dendrobium and Oncidium. In this review, we consolidate recent progress in research on the floral transition and flower development of orchids emphasizing representative genes and genetic networks, and also introduce a few successful cases of manipulation of orchid flowering/flower development through the application of molecular breeding or biotechnology tools.

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.

Analysis of the TCP genes expressed in the inflorescence of the orchid Orchis italica

TCP proteins are plant-specific transcription factors involved in many different processes. Because of their involvement in a large number of developmental pathways, their roles have been investigated in various plant species. However, there are almost no studies of this transcription factor family in orchids. Based on the available transcriptome of the inflorescence of the orchid Orchis italica, in the present study we identified 12 transcripts encoding TCP proteins. The phylogenetic analysis showed that they belong to different TCP classes (I and II) and groups (PCF, CIN and CYC/TB1), and that they display a number of conserved motifs when compared with the TCPs of Arabidopsis and Oryza. The presence of a specific cleavage site for the microRNA miRNA319, an important post-transcriptional regulator of several TCP genes in other species, was demonstrated for one transcript of O. italica, and the analysis of the expression pattern of the TCP transcripts in different inflorescence organs and in leaf tissue suggests that some TCP transcripts of O. italica exert their role only in specific tissues, while others may play multiple roles in different tissues. In addition, the evolutionary analysis showed a general purifying selection acting on the coding region of these transcripts.

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.

The analysis of the inflorescence miRNome of the orchid Orchis italica reveals a DEF-like MADS-box gene as a new miRNA target

Plant microRNAs (miRNAs) are small, regulatory non-coding RNAs involved in a wide range of biological processes, from organ development to response to stimuli. In recent years, an increasing number of studies on model plant species have highlighted the evolutionary conservation of a high number of miRNA families and the existence of taxon-specific ones. However, few studies have examined miRNAs in non-model species such as orchids, which are characterized by highly diversified floral structures and pollination strategies. Therefore, we analysed a small RNA library of inflorescence tissue of the Mediterranean orchid Orchis italica to increase the knowledge on miRNAs in a non-model plant species. The high-throughput sequencing and analysis of a small RNA library of inflorescence of O. italica revealed 23 conserved and 161 putative novel miRNA families. Among the putative miRNA targets, experimental validation demonstrated that a DEF-like MADS-box transcript is cleaved by the homolog of miR5179 of O. italica. The presence of conserved miRNA families in the inflorescence of O. italica indicates that the basic developmental flower regulatory mechanisms mediated by miRNAs are maintained through evolution. Because, according to the "orchid code" theory, DEF-like genes exert a key function in the diversification of tepals and lip, the cleavage-mediated inhibitory activity of miR5179 on a OitaDEF-like transcript suggests that, in orchids, miRNAs play an important role in the diversification of the perianth organs.

A modified ABCDE model of flowering in orchids based on gene expression profiling studies of the moth orchid Phalaenopsis aphrodite

Previously we developed genomic resources for orchids, including transcriptomic analyses using next-generation sequencing techniques and construction of a web-based orchid genomic database. Here, we report a modified molecular model of flower development in the Orchidaceae based on functional analysis of gene expression profiles in Phalaenopsis aphrodite (a moth orchid) that revealed novel roles for the transcription factors involved in floral organ pattern formation. Phalaenopsis orchid floral organ-specific genes were identified by microarray analysis. Several critical transcription factors including AP3, PI, AP1 and AGL6, displayed distinct spatial distribution patterns. Phylogenetic analysis of orchid MADS box genes was conducted to infer the evolutionary relationship among floral organ-specific genes. The results suggest that gene duplication MADS box genes in orchid may have resulted in their gaining novel functions during evolution. Based on these analyses, a modified model of orchid flowering was proposed. Comparison of the expression profiles of flowers of a peloric mutant and wild-type Phalaenopsis orchid further identified genes associated with lip morphology and peloric effects. Large scale investigation of gene expression profiles revealed that homeotic genes from the ABCDE model of flower development classes A and B in the Phalaenopsis orchid have novel functions due to evolutionary diversification, and display differential expression patterns.

The AP2-like gene OitaAP2 is alternatively spliced and differentially expressed in inflorescence and vegetative tissues of the orchid Orchis italica

The AP2/ERF proteins are plant-specific transcription factors involved in multiple regulatory pathways, from plant organ development to response to various environmental stresses. One of the mechanisms that regulates the AP2-like genes involves the microRNA miR172, which controls their activity at the post-transcriptional level. Extensive studies on AP2-like genes are available in many different species; however, in orchids, one of the largest plant families, studies are restricted to a few species, all belonging to the Epidendroideae subfamily. In the present study, we report the isolation of an AP2-like gene in the Mediterranean orchid Orchis italica (Orchidoideae). The OitaAP2 locus includes 10 exons and 9 introns, and its transcript is alternatively spliced, resulting in the long OitaAP2 and the short OitaAP2_ISO isoforms, with the latter skipping exon 9. Both isoforms contain the conserved target site for miR172, whose action is demonstrated by the presence of cleaved OitaAP2 mRNA. The OitaAP2 and OitaAP2_ISO mRNAs are present in the tepals and lip before and after anthesis at different expression levels. In addition, the OitaAP2_ISO isoform is expressed in the ovary before pollination and in the root and stem. The isoform-specific expression pattern suggests a functional differentiation of the OitaAP2 alternatively spliced transcripts. The expression profile of miR172 is complementary to that of the OitaAP2 isoforms in inflorescence tissues before anthesis, whereas after anthesis and in ovary tissue before and after pollination, this relationship disappears, suggesting the existence of OitaAP2 inhibitory mechanisms in these tissues that differ from that involving miR172.

The OitaAG and OitaSTK genes of the orchid Orchis italica: a comparative analysis with other C- and D-class MADS-box genes

According to the ABCDE model of flower development, the C- and D- class MADS box genes are involved in the formation of male and female reproductive organs (fused to form the column in orchids) and in ovule maturation (triggered by fertilization in orchids). In the present study, we report the isolation of the Orchis italica genes OitaAG and OitaSTK, homologs of the C-class AGAMOUS and the D-class SEEDSTICK genes of Arabidopsis, respectively. Analysis of their expression profiles reveals high levels of mRNA in columns and ovaries, particularly after pollination. However, weak expression is also detectable in the inner tepals (OitaAG) and the lip and root (OitaSTK). This expression profile is only partially overlapping with those reported in other orchid species and may be the consequence of a different evolutionary history of these functional gene classes in orchids. The genomic characterization of the OitaAG and OitaSTK genes shows that a high number of traces of mobile elements are present in introns and could have contributed to the size expansion of some of them (e.g., intron 2 and 3 of OitaAG and intron 3, 4 and 5 of OitaSTK). Nucleotide sequences of intron 1 of the OitaSTK gene and other STK-like genes do not share regulatory motifs, whereas sequence comparison of intron 2 of the OitaAG gene with that of intron 2 of other AG-like genes reveals, for the first time in an orchid species, the presence of conserved cis-regulatory boxes and binding sites for transcription factors that positively (e.g., LEAFY and WUSCHEL) or negatively (e.g., BELLRINGER) regulate the expression of the AG homologs in dicots and monocots.

C- and D-class MADS-box genes from Phalaenopsis equestris (Orchidaceae) display functions in gynostemium and ovule development

Gynostemium and ovule development in orchid are unique developmental processes in the plant kingdom. Characterization of C- and D-class MADS-box genes could help reveal the molecular mechanisms underlying gynostemium and ovule development in orchids. In this study, we isolated and characterized a C- and a D-class gene, PeMADS1 and PeMADS7, respectively, from Phalaenopsis equestris. These two genes showed parallel spatial and temporal expression profiles, which suggests their cooperation in gynostemium and ovule development. Furthermore, only PeMADS1 was ectopically expressed in the petals of the gylp (gynostemium-like petal) mutant, whose petals were transformed into gynostemium-like structures. Protein-protein interaction analyses revealed that neither PeMADS1 and PeMADS7 could form a homodimer or a heterodimer. An E-class protein was needed to bridge the interaction between these two proteins. A complementation test revealed that PeMADS1 could rescue the phenotype of the AG mutant. Overexpression of PeMADS7 in Arabidopsis caused typical phenotypes of the D-class gene family. Together, these results indicated that both C-class PeMADS1 and D-class PeMADS7 play important roles in orchid gynostemium and ovule development.