Wen Chieh Tsai

The genome sequence of the orchid Phalaenopsis equestris

Orchidaceae, renowned for its spectacular flowers and other reproductive and ecological adaptations, is one of the most diverse plant families. Here we present the genome sequence of the tropical epiphytic orchid Phalaenopsis equestris, a frequently used parent species for orchid breeding. P. equestris is the first plant with crassulacean acid metabolism (CAM) for which the genome has been sequenced. Our assembled genome contains 29,431 predicted protein-coding genes. We find that contigs likely to be underassembled, owing to heterozygosity, are enriched for genes that might be involved in self-incompatibility pathways. We find evidence for an orchid-specific paleopolyploidy event that preceded the radiation of most orchid clades, and our results suggest that gene duplication might have contributed to the evolution of CAM photosynthesis in P. equestris. Finally, we find expanded and diversified families of MADS-box C/D-class, B-class AP3 and AGL6-class genes, which might contribute to the highly specialized morphology of orchid flowers.

Strategies for Functional Validation of Genes Involved in Reproductive Stages of Orchids

Plants in the largest family of angiosperms, Orchidaceae, are diverse in both specialized pollination and ecological strategies and provide a rich source for investigating evolutionary relationships and developmental biology. However, studies in orchids have been hindered by several challenges that include low transformation efficiency and long regeneration time. To overcome such obstacles, we selected a symptomless cymbidium mosaic virus (CymMV) isolate for constructing virus-induced gene-silencing vectors. The feasibility of the virus vectors was first assessed with use of an orchid phytoene desaturase gene. The vector was able to induce gene silencing in orchids; however, because of the slow growth of orchids, the commonly used phytoene desaturase gene was not a good visual marker in orchids. We inserted a 150-nucleotide unique region of a B-class MADS-box family gene, PeMADS6, into pCymMV-pro60. The transcription level of PeMADS6 in inoculated Phalaenopsis plants was reduced by up to 73%, but no effect was observed for other MADS-box family genes. In contrast, in Phalaenopsis plants inoculated with CymMV transcripts containing 500 nucleotides of PeMADS6, a conserved region among MADS-box genes, the transcription level of PeMADS6 and the B- and C-class MADS-box genes was reduced by up to 97.8% as compared with plants inoculated with the vector alone. Flower morphology was affected in the MADS-box family gene-silenced plants as well. This in vivo experiment demonstrates an efficient way to study genes involved in the reproductive stage of plants with a long life cycle.

OrchidBase: A Collection of Sequences of the Transcriptome Derived from Orchids

Orchids are one of the most ecological and evolutionarily significant plants, and the Orchidaceae is one of the most abundant families of the angiosperms. Genetic databases will be useful not only for gene discovery but also for future genomic annotation. For this purpose, OrchidBase was established from 37,979,342 sequence reads collected from 11 in-house Phalaenopsis orchid cDNA libraries. Among them, 41,310 expressed sequence tags (ESTs) were obtained by using Sanger sequencing, whereas 37,908,032 reads were obtained by using next-generation sequencing (NGS) including both Roche 454 and Solexa Illumina sequencers. These reads were assembled into 8,501 contigs and 76,116 singletons, resulting in 84,617 non-redundant transcribed sequences with an average length of 459 bp. The analysis pipeline of the database is an automated system written in Perl and C#, and consists of the following components: automatic pre-processing of EST reads, assembly of raw sequences, annotation of the assembled sequences and storage of the analyzed information in SQL databases. A web application was implemented with HTML and a Microsoft .NET Framework C# program for browsing and querying the database, creating dynamic web pages on the client side, analyzing gene ontology (GO) and mapping annotated enzymes to KEGG pathways. The online resources for putative annotation can be searched either by text or by using BLAST, and the results can be explored on the website and downloaded. Consequently, the establishment of OrchidBase will provide researchers with a high-quality genetic resource for data mining and facilitate efficient experimental studies on orchid biology and biotechnology. The OrchidBase database is freely available at http://lab.fhes.tn.edu.tw/est.

Comparison of transcripts in Phalaenopsis bellina and Phalaenopsis equestris (Orchidaceae) flowers to deduce monoterpene biosynthesis pathway

BackgroundFloral scent is one of the important strategies for ensuring fertilization and for determining seed or fruit set. Research on plant scents has hampered mainly by the invisibility of this character, its dynamic nature, and complex mixtures of components that are present in very small quantities. Most progress in scent research, as in other areas of plant biology, has come from the use of molecular and biochemical techniques. Although volatile components have been identified in several orchid species, the biosynthetic pathways of orchid flower fragrance are far from understood. We investigated how flower fragrance was generated in certain Phalaenopsis orchids by determining the chemical components of the floral scent, identifying floral expressed-sequence-tags (ESTs), and deducing the pathways of floral scent biosynthesis in Phalaneopsis bellina by bioinformatics analysis.ResultsThe main chemical components in the P. bellina flower were shown by gas chromatography-mass spectrometry to be monoterpenoids, benzenoids and phenylpropanoids. The set of floral scent producing enzymes in the biosynthetic pathway from glyceraldehyde-3-phosphate (G3P) to geraniol and linalool were recognized through data mining of the P. bellina floral EST database (dbEST). Transcripts preferentially expressed in P. bellina were distinguished by comparing the scent floral dbEST to that of a scentless species, P. equestris, and included those encoding lipoxygenase, epimerase, diacylglycerol kinase and geranyl diphosphate synthase. In addition, EST filtering results showed that transcripts encoding signal transduction and Myb transcription factors and methyltransferase, in addition to those for scent biosynthesis, were detected by in silico hybridization of the P. bellina unigene database against those of the scentless species, rice and Arabidopsis. Altogether, we pinpointed 66% of the biosynthetic steps from G3P to geraniol, linalool and their derivatives.ConclusionThis systems biology program combined chemical analysis, genomics and bioinformatics to elucidate the scent biosynthesis pathway and identify the relevant genes. It integrates the forward and reverse genetic approaches to knowledge discovery by which researchers can study non-model plants.

A novel homodimeric geranyl diphosphate synthase from the orchid Phalaenopsis bellina lacking a DD(X)2–4D motif

SUMMARY Geranyl diphosphate (GDP) is the precursor of monoterpenes, which are the major floral scent compounds in Phalaenopsis bellina. The cDNA of P. bellina GDP synthase (PbGDPS) was cloned, and its sequence corresponds to the second Asp-rich motif (SARM), but not to any aspartate-rich (Asp-rich) motif. The recombinant PbGDPS enzyme exhibits dual prenyltransferase activity, producing both GDP and farnesyl diphosphate (FDP), and a yeast two-hybrid assay and gel filtration revealed that PbGDPS was able to form a homodimer. Spatial and temporal expression analyses showed that the expression of PbGDPS was flower specific, and that maximal PbGDPS expression was concomitant with maximal emission of monoterpenes on day 5 post-anthesis. Homology modelling of PbGDPS indicated that the Glu-rich motif might provide a binding site for Mg(2+) and catalyze the formation of prenyl products in a similar way to SARM. Replacement of the key Glu residues with alanine totally abolished enzyme activity, whereas their mutation to Asp resulted in a mutant with two-thirds of the activity of the wild-type protein. Phylogenetic analysis indicated that plant GDPS proteins formed four clades: members of both GDPS-a and GDPS-b clades contain Asp-rich motifs, and function as homodimers. In contrast, proteins in the GDPS-c and GDPS-d clades do not contain Asp-rich motifs, but although members of the GDPS-c clade function as heterodimers, PbGDPS, which is more closely related to the GDPS-c clade proteins than to GDPS-a and GDPS-b proteins, and is currently the sole member of the GDPS-d clade, functions as a homodimer.

Research on Orchid Biology and Biotechnology

Orchidaceae constitute one of the largest families of angiosperms. They are one of the most ecological and evolutionary significant plants and have successfully colonized almost every habitat on earth. Because of the significance of plant biology, market needs and the current level of breeding technologies, basic research into orchid biology and the application of biotechnology in the orchid industry are continually endearing scientists to orchids in Taiwan. In this introductory review, we give an overview of the research activities in orchid biology and biotechnology, including the status of genomics, transformation technology, flowering regulation, molecular regulatory mechanisms of floral development, scent production and color presentation. This information will provide a broad scope for study of orchid biology and serve as a starting point for uncovering the mysteries of orchid evolution.

Interactions of B-class complex proteins involved in tepal development in Phalaenopsis orchid

In our previous studies, we identified four DEFICIENS (DEF)-like genes and one GLOBOSA (GLO)-like gene involved in floral organ development in Phalaenopsis equestris. Revealing the DNA binding properties and protein-protein interactions of these floral homeotic MADS-box protein complexes (PeMADS) in orchids is crucial for the elucidation of the unique orchid floral morphogenesis. In this study, the interactome of B-class PeMADS proteins was assayed by the yeast two-hybrid system (Y2H) and glutathione S-transferase (GST) pull-down assays. Furthermore, the DNA binding activities of these proteins were assessed by using electrophoretic mobility shift assay (EMSA). All four DEF-like PeMADS proteins interacted individually with the GLO-like PeMADS6 in Y2H assay, yet with different strengths of interaction. Generally, the PeMADS3/PeMADS4 lineage interacted more strongly with PeMADS6 than the PeMADS2/PeMADS5 lineage did. In addition, independent homodimer formation for both PeMADS4 (DEF-like) and PeMADS6 (GLO-like) was detected. The protein-protein interactions between pairs of PeMADS proteins were further confirmed by using a GST pull-down assay. Furthermore, both the PeMADS4 homodimer and the PeMADS6 homodimer/homomultimer per se were able to bind to the MADS-box protein-binding motif CArG. The heterodimeric complexes PeMADS2-PeMADS6, PeMADS4-PeMADS6 and PeMADS5-PeMADS6 showed CArG binding activity. Taken together, these results suggest that various complexes formed among different combinations of the five B-class PeMADS proteins may increase the complexity of their regulatory functions and thus specify the molecular basis of whorl morphogenesis and combinatorial interactions of floral organ identity genes in orchids.

An overview of the Phalaenopsis orchid genome through BAC end sequence analysis

BackgroundPhalaenopsis orchids are popular floral crops, and development of new cultivars is economically important to floricultural industries worldwide. Analysis of orchid genes could facilitate orchid improvement. Bacterial artificial chromosome (BAC) end sequences (BESs) can provide the first glimpses into the sequence composition of a novel genome and can yield molecular markers for use in genetic mapping and breeding.ResultsWe used two BAC libraries (constructed using the BamHI and HindIII restriction enzymes) of Phalaenopsis equestris to generate pair-end sequences from 2,920 BAC clones (71.4% and 28.6% from the BamHI and HindIII libraries, respectively), at a success rate of 95.7%. A total of 5,535 BESs were generated, representing 4.5 Mb, or about 0.3% of the Phalaenopsis genome. The trimmed sequences ranged from 123 to 1,397 base pairs (bp) in size, with an average edited read length of 821 bp. When these BESs were subjected to sequence homology searches, it was found that 641 (11.6%) were predicted to represent protein-encoding regions, whereas 1,272 (23.0%) contained repetitive DNA. Most of the repetitive DNA sequences were gypsy- and copia-like retrotransposons (41.9% and 12.8%, respectively), whereas only 10.8% were DNA transposons. Further, 950 potential simple sequence repeats (SSRs) were discovered. Dinucleotides were the most abundant repeat motifs; AT/TA dimer repeats were the most frequent SSRs, representing 253 (26.6%) of all identified SSRs. Microsynteny analysis revealed that more BESs mapped to the whole-genome sequences of poplar than to those of grape or Arabidopsis, and even fewer mapped to the rice genome. This work will facilitate analysis of the Phalaenopsis genome, and will help clarify similarities and differences in genome composition between orchids and other plant species.ConclusionUsing BES analysis, we obtained an overview of the Phalaenopsis genome in terms of gene abundance, the presence of repetitive DNA and SSR markers, and the extent of microsynteny with other plant species. This work provides a basis for future physical mapping of the Phalaenopsis genome and advances our knowledge thereof.

Gene discovery using next-generation pyrosequencing to develop ESTs for Phalaenopsis orchids

BackgroundOrchids are one of the most diversified angiosperms, but few genomic resources are available for these non-model plants. In addition to the ecological significance, Phalaenopsis has been considered as an economically important floriculture industry worldwide. We aimed to use massively parallel 454 pyrosequencing for a global characterization of the Phalaenopsis transcriptome.ResultsTo maximize sequence diversity, we pooled RNA from 10 samples of different tissues, various developmental stages, and biotic- or abiotic-stressed plants. We obtained 206,960 expressed sequence tags (ESTs) with an average read length of 228 bp. These reads were assembled into 8,233 contigs and 34,630 singletons. The unigenes were searched against the NCBI non-redundant (NR) protein database. Based on sequence similarity with known proteins, these analyses identified 22,234 different genes (E-value cutoff, e-7). Assembled sequences were annotated with Gene Ontology, Gene Family and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Among these annotations, over 780 unigenes encoding putative transcription factors were identified.ConclusionPyrosequencing was effective in identifying a large set of unigenes from Phalaenopsis. The informative EST dataset we developed constitutes a much-needed resource for discovery of genes involved in various biological processes in Phalaenopsis and other orchid species. These transcribed sequences will narrow the gap between study of model organisms with many genomic resources and species that are important for ecological and evolutionary studies.

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.