A de novo floral transcriptome reveals clues into Phalaenopsis orchid flower development

Breeding of ornamental orchids with focus on Phalaenopsis: current approaches, tools, and challenges for this century

Ornamental orchid breeding programs have been conducted to develop commercially valuable cultivars with improved characteristics of commercial interest, such as size, flower color, pattern, shape, and resistance to pathogens. Conventional breeding, including sexual hybridization followed by selection of desirable characteristics in plants, has so far been the main method for ornamental breeding, but other techniques, including mutation induction by polyploidization and gamma irradiation, and biotechnological techniques, such as genetic transformation, have also been studied and used in ornamental breeding programs. Orchids are one of the most commercially important families in floriculture industry, having very particular reproductive biology characteristics and being a well-studied group of ornamentals in terms of genetic improvement. The present review focuses on the conventional and biotechnological techniques and approaches specially employed in breeding Phalaenopsis orchids, the genus with highest worldwide importance as an ornamental orchid, highlighting the main limitations and strengths of the approaches. Furthermore, new opportunities and future prospects for ornamental breeding in the CRISPR/Cas9 genome editing era are also discussed. We conclude that conventional hybridization remains the most used method to obtain new cultivars in orchids. However, the emergence of the first biotechnology-derived cultivars, as well as the new biotechnological tools available, such as CRISPR-Cas9, rekindled the full potential of biotechnology approaches and their importance for improve ornamental orchid breeding programs.

Transcriptomic Analysis Suggests Auxin Regulation in Dorsal-Ventral Petal Asymmetry of Wild Progenitor Sinningia speciosa

The establishment of dorsal-ventral (DV) petal asymmetry is accompanied by differential growth of DV petal size, shape, and color differences, which enhance ornamental values. Genes involved in flower symmetry in Sinningia speciosa have been identified as CYCLOIDEA (SsCYC), but which gene regulatory network (GRN) is associated with SsCYC to establish DV petal asymmetry is still unknown. To uncover the GRN of DV petal asymmetry, we identified 630 DV differentially expressed genes (DV-DEGs) from the RNA-Seq of dorsal and ventral petals in the wild progenitor, S. speciosa 'ES'. Validated by qRT-PCR, genes in the auxin signaling transduction pathway, SsCYC, and a major regulator of anthocyanin biosynthesis were upregulated in dorsal petals. These genes correlated with a higher endogenous auxin level in dorsal petals, with longer tube length growth through cell expansion and a purple dorsal color. Over-expression of SsCYC in Nicotiana reduced petal size by regulating cell growth, suggesting that SsCYC also controls cell expansion. This suggests that auxin and SsCYC both regulate DV petal asymmetry. Transiently over-expressed SsCYC, however, could not activate most major auxin signaling genes, suggesting that SsCYC may not trigger auxin regulation. Whether auxin can activate SsCYC or whether they act independently to regulate DV petal asymmetry remains to be explored in the future.

Whole-transcriptome analysis of differentially expressed genes between ray and disc florets and identification of flowering regulatory genes in Chrysanthemum morifolium

Chrysanthemum morifolium has ornamental and economic values. However, there has been minimal research on the morphology of the chrysanthemum florets and related genes. In this study, we used the leaves as a control to screen for differentially expressed genes between ray and disc florets in chrysanthemum flowers. A total of 8,359 genes were differentially expressed between the ray and disc florets, of which 3,005 were upregulated and 5,354 were downregulated in the disc florets. Important regulatory genes that control flower development and flowering determination were identified. Among them, we identified a TM6 gene (CmTM6-mu) that belongs to the Class B floral homeotic MADS-box transcription factor family, which was specifically expressed in disc florets. We isolated this gene and found it was highly similar to other typical TM6 lineage genes, but a single-base deletion at the 3' end of the open reading frame caused a frame shift that generated a protein in which the TM6-specific paleoAP3 motif was missing at the C terminus. The CmTM6-mu gene was ectopically expressed in Arabidopsis thaliana. Petal and stamen developmental processes were unaffected in transgenic A. thaliana lines; however, the flowering time was earlier than in the wild-type control. Thus, the C-terminal of paleoAP3 appears to be necessary for the functional performance in regulating the development of petals or stamens and CmTM6-mu may be involved in the regulation of flowering time in chrysanthemum. The results of this study will be useful for future research on flowering molecular mechanisms and for the breeding of novel flower types.

Extending the Toolkit for Beauty: Differential Co-Expression of DROOPING LEAF-Like and Class B MADS-Box Genes during Phalaenopsis Flower Development

The molecular basis of orchid flower development is accomplished through a specific regulatory program in which the class B MADS-box AP3/DEF genes play a central role. In particular, the differential expression of four class B AP3/DEF genes is responsible for specification of organ identities in the orchid perianth. Other MADS-box genes (AGL6 and SEP-like) enrich the molecular program underpinning the orchid perianth development, resulting in the expansion of the original “orchid code” in an even more complex gene regulatory network. To identify candidates that could interact with the AP3/DEF genes in orchids, we conducted an in silico differential expression analysis in wild-type and peloric Phalaenopsis. The results suggest that a YABBY DL-like gene could be involved in the molecular program leading to the development of the orchid perianth, particularly the labellum. Two YABBY DL/CRC homologs are present in the genome of Phalaenopsis equestris, PeDL1 and PeDL2, and both express two alternative isoforms. Quantitative real-time PCR analyses revealed that both genes are expressed in column and ovary. In addition, PeDL2 is more strongly expressed the labellum than in the other tepals of wild-type flowers. This pattern is similar to that of the AP3/DEF genes PeMADS3/4 and opposite to that of PeMADS2/5. In peloric mutant Phalaenopsis, where labellum-like structures substitute the lateral inner tepals, PeDL2 is expressed at similar levels of the PeMADS2-5 genes, suggesting the involvement of PeDL2 in the development of the labellum, together with the PeMADS2-PeMADS5 genes. Although the yeast two-hybrid analysis did not reveal the ability of PeDL2 to bind the PeMADS2-PeMADS5 proteins directly, the existence of regulatory interactions is suggested by the presence of CArG-boxes and other MADS-box transcription factor binding sites within the putative promoter of the orchid DL2 gene.

Comparative Transcriptome Reveals the Genes' Adaption to Herkogamy of Lumnitzera littorea (Jack) Voigt

Lumnitzera littorea (Jack) Voigt is among the most endangered mangrove species in China. The morphology and evolution of L. littorea flowers have received substantial attention for their crucial reproductive functions. However, little is known about the genomic regulation of flower development in L. littorea. In this study, we characterized the morphology of two kinds of L. littorea flowers and performed comparative analyses of transcriptome profiles of the two different flowers. Morphological observation showed that some flowers have a column embedded in the petals while others produce a stretched flower style during petal unfolding in flowering. By using RNA-seq, we obtained 138,857 transcripts that were assembled into 82,833 unigenes with a mean length of 1055.48 bp. 82,834 and 34,997 unigenes were assigned to 52 gene ontology (GO) functional groups and 364 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. A total of 4,267 differentially expressed genes (DEGs), including 1,794 transcription factors (TFs), were identified between two types of flowers. These TFs are mainly involved in bHLH, B3, bZIP, MYB-related, and NAC family members. We further validated that 12 MADS-box genes, including 4 MIKC-type and 8 M-type TFs, were associated with the pollinate of L. littorea by herkogamy. Our current results provide valuable information for genetic analysis of L. littorea flowering and may be useful for illuminating its adaptive evolutionary mechanisms.

Proteomic Analysis of the Early Development of the Phalaenopsis amabilis Flower Bud under Low Temperature Induction Using the iTRAQ/MRM Approach

Phalaenopsis amabilis, one of the most important plants in the international flower market due to its graceful shape and colorful flowers, is an orchid that undergoes vernalization and requires low-temperature treatment for flowering. There have been few reports on the proteomics of the development of flower buds. In this study, isobaric tags for relative and absolute quantification (iTRAQ) were used to identify 5064 differentially expressed proteins in P. amabilis under low-temperature treatment; of these, 42 were associated with early floral induction, and 18 were verified by mass spectrometry multi-reaction monitoring (MRM). The data are available via ProteomeXchange under identifier PXD013908. Among the proteins associated with the vernalization pathway, PEQU_11434 (glycine-rich RNA-binding protein GRP1A-like) and PEQU_19304 (FT, VRN3 homolog) were verified by MRM, and some other important proteins related to vernalization and photoperiod pathway that were detected by iTRAQ but not successfully verified by MRM, such as PEQU_11045 (UDP-N-acetylglucosamine diphosphorylase), phytochromes A (PEQU_13449, PEQU_35378), B (PEQU_09249), and C (PEQU_41401). Our data revealed a regulation network of the early development of flower buds in P. amabilis under low temperature induction.

PeERF1, a SHINE-Like Transcription Factor, Is Involved in Nanoridge Development on Lip Epidermis of Phalaenopsis Flowers

Phalaenopsis orchids have a spectacular floral morphology with a highly evolved lip that offers a landing platform for pollinators. The typical morphological orchid lip features are essential for the special pollination mechanism of Phalaenopsis flowers. Previously, we found that in the lip, a member of the AP2/EREBP protein family was highly expressed. Here, we further confirmed its high expression and characterized its function during lip development. Phylogenetic analysis showed that AP2/EREBP belongs to the Va2 subgroup of ERF transcription factors. We named it PeERF1. We found that PeERF1 was only expressed at stage 5, as flowers opened. This coincided with both thickening of the cuticle and development of nanoridges. We performed knockdown expression of PeERF1 using CymMV-based virus-induced gene silencing in either the AP2 conserved domain, producing PeERF1_AP2-silenced plants, or the SHN specific domain, producing PeERF1_SHN-silenced plants. Using cryo-SEM, we found that the number of nanoridges was reduced only in the PeERF1_AP2-silenced group. This change was found on both the abaxial and adaxial surfaces of the central lip lobe. Expression of PeERF1 was reduced significantly in PeERF1_AP2-silenced plants. In cutin biosynthesis genes, expression of both PeCYP86A2 and PeDCR was significantly decreased in both groups. The expression of PeCYP77A4 was reduced significantly only in the PeERF1_AP2-silenced plants. Although PeGPAT expression was reduced in both silenced plants, but to a lesser degree. The expression of PeERF1 was significantly reduced in the petal-like lip of a big-lip variant. PeCYP77A4 and PeGPAT in the lip were also reduced, but PeDCR was not. Furthermore, heterologous overexpression of PeERF1 in the genus Arabidopsis produced leaves that were shiny on the adaxial surface. Taken together, our results show that in Phalaenopsis orchids PeERF1 plays an important role in formation of nanoridges during lip epidermis development.

Comparative transcriptomics provides insight into the molecular basis of species diversification of section Trigonopedia (Cypripedium) on the Qinghai-Tibetan Plateau

Deceptive pollination is key to the species richness of Orchidaceae. However, the genetic basis of species diversification is still under study. Section Trigonopedia is a monophyletic clade of genus Cypripedium distributed in the southwest of China. The species of this section are pollinated by different flies. Pollinator differentiation makes section Trigonopedia an ideal group for studying the genetic basis underlying species diversification. Here, we sequenced the transcriptomes of eight species of the genus Cypripedium, including six co-flowering species of section Trigonopedia and two species outside this section as an outgroup. We reconstructed the phylogeny of the section with the combined 1572 single-copy genes extracted from the eight species and produced a highly resolved tree of the section. Furthermore, we combined substitution rate estimation and differential expression analysis to identify candidate genes, including genes related to floral scent synthesis and environmental adaptation, involved in species differentiation. Field investigations showed that these species have adapted to different habitats. We propose that the species diversification in this section is initiated by floral scent differentiation, followed by habitat differentiation, finally leading to speciation. This study sheds novel light on the diversification of closely related orchid species in the Qinghai-Tibetan region.

Post genomics era for orchid research

Among 300,000 species in angiosperms, Orchidaceae containing 30,000 species is one of the largest families. Almost every habitats on earth have orchid plants successfully colonized, and it indicates that orchids are among the plants with significant ecological and evolutionary importance. So far, four orchid genomes have been sequenced, including Phalaenopsis equestris, Dendrobium catenatum, Dendrobium officinale, and Apostaceae shengen. Here, we review the current progress and the direction of orchid research in the post genomics era. These include the orchid genome evolution, genome mapping (genome-wide association analysis, genetic map, physical map), comparative genomics (especially receptor-like kinase and terpene synthase), secondary metabolomics, and genome editing.

Tissue-Specific Floral Transcriptome Analysis of the Sexually Deceptive Orchid Chiloglottis trapeziformis Provides Insights into the Biosynthesis and Regulation of Its Unique UV-B Dependent Floral Volatile, Chiloglottone 1

The Australian sexually deceptive orchid, Chiloglottis trapeziformis, employs a unique UV-B-dependent floral volatile, chiloglottone 1, for specific male wasp pollinator attraction. Chiloglottone 1 and related variants (2,5-dialkylcyclohexane-1,3-diones), represent a unique class of specialized metabolites presumed to be the product of cyclization between two fatty acid (FA) precursors. However, the genes involved in the biosynthesis of precursors, intermediates, and transcriptional regulation remains to be discovered. Chiloglottone 1 production occurs in the aggregation of calli (callus) on the labellum under continuous UV-B light. Therefore, deep sequencing, transcriptome assembly, and differential expression (DE) analysis were performed across different tissue types and UV-B treatments. Transcripts expressed in the callus and labellum (∼23,000 transcripts) were highly specialized and enriched for a diversity of known and novel metabolic pathways. DE analysis between chiloglottone-emitting callus versus the remainder of the labellum showed strong coordinated induction of entire FA biosynthesis and β-oxidation pathways including genes encoding Ketoacyl-ACP Synthase, Acyl-CoA Oxidase, and Multifunctional Protein. Phylogenetic analysis revealed potential gene duplicates with tissue-specific differential regulation including two Acyl-ACP Thioesterase B and a Ketoacyl-ACP Synthase genes. UV-B treatment induced the activation of UVR8-mediated signaling and large-scale transcriptome changes in both tissues, however, neither FA biosynthesis/β-oxidation nor other lipid metabolic pathways showed clear indications of concerted DE. Gene co-expression network analysis identified three callus-specific modules enriched with various lipid metabolism categories. These networks also highlight promising candidates involved in the cyclization of chiloglottone 1 intermediates (e.g., Bet v I and dimeric α,β barrel proteins) and orchestrating regulation of precursor pathways (e.g., AP2/ERF) given a strong co-regulation with FA biosynthesis/β-oxidation genes. Possible alternative biosynthetic routes for precursors (e.g., aldehyde dehydrogenases) were also indicated. Our comprehensive study constitutes the first step toward understanding the biosynthetic pathways involved in chiloglottone 1 production in Chiloglottis trapeziformis - supporting the roles of FA metabolism in planta, gene duplication as a potential source of new genes, and co-regulation of novel pathway genes in a tissue-specific manner. This study also provides a new and valuable resource for future discovery and comparative studies in plant specialized metabolism of other orchids and non-model plants.

De novo transcriptome assembly databases for the butterfly orchid Phalaenopsis equestris

Orchids are renowned for their spectacular flowers and ecological adaptations. After the sequencing of the genome of the tropical epiphytic orchid Phalaenopsis equestris, we combined Illumina HiSeq2000 for RNA-Seq and Trinity for de novo assembly to characterize the transcriptomes for 11 diverse P. equestris tissues representing the root, stem, leaf, flower buds, column, lip, petal, sepal and three developmental stages of seeds. Our aims were to contribute to a better understanding of the molecular mechanisms driving the analysed tissue characteristics and to enrich the available data for P. equestris. Here, we present three databases. The first dataset is the RNA-Seq raw reads, which can be used to execute new experiments with different analysis approaches. The other two datasets allow different types of searches for candidate homologues. The second dataset includes the sets of assembled unigenes and predicted coding sequences and proteins, enabling a sequence-based search. The third dataset consists of the annotation results of the aligned unigenes versus the Nonredundant (Nr) protein database, Kyoto Encyclopaedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG) databases with low e-values, enabling a name-based search.

The genome and transcriptome of Phalaenopsis yield insights into floral organ development and flowering regulation

The Phalaenopsis orchid is an important potted flower of high economic value around the world. We report the 3.1 Gb draft genome assembly of an important winter flowering Phalaenopsis ‘KHM190’ cultivar. We generated 89.5 Gb RNA-seq and 113 million sRNA-seq reads to use these data to identify 41,153 protein-coding genes and 188 miRNA families. We also generated a draft genome for Phalaenopsis pulcherrima ‘B8802,’ a summer flowering species, via resequencing. Comparison of genome data between the two Phalaenopsis cultivars allowed the identification of 691,532 single-nucleotide polymorphisms. In this study, we reveal that the key role of PhAGL6b in the regulation of labellum organ development involves alternative splicing in the big lip mutant. Petal or sepal overexpressing PhAGL6b leads to the conversion into a lip-like structure. We also discovered that the gibberellin pathway that regulates the expression of flowering time genes during the reproductive phase change is induced by cool temperature. Our work thus depicted a valuable resource for the flowering control, flower architecture development, and breeding of the Phalaenopsis orchids.