Evidence for the Involvement of Fatty Acid Biosynthesis and Degradation in the Formation of Insect Sex Pheromone-Mimicking Chiloglottones in Sexually Deceptive Chiloglottis Orchids

Hundreds of orchid species secure pollination by sexually luring specific male insects as pollinators by chemical and morphological mimicry. Yet, the biochemical pathways involved in the synthesis of the insect sex pheromone-mimicking volatiles in these sexually deceptive plants remain poorly understood. Here, we explore the biochemical pathways linked to the chemical mimicry of female sex pheromones (chiloglottones) employed by the Australian sexually deceptive Chiloglottis orchids to lure their male pollinator. By strategically exploiting the transcriptomes of chiloglottone 1-producing Chiloglottis trapeziformis at distinct floral tissues and at key floral developmental stages, we identified two key transcriptional trends linked to the stage- and tissue-dependent distribution profiles of chiloglottone in the flower: (i) developmental upregulation of fatty acid biosynthesis and β-oxidation genes such as KETOACYL-ACP SYNTHASE, FATTY ACYL-ACP THIOESTERASE, and ACYL-COA OXIDASE during the transition from young to mature buds and flowers and (ii) the tissue-specific induction of fatty acid pathway genes in the callus (the insectiform odor-producing structure on the labellum of the flower) compared to the labellum remains (non-odor-producing) regardless of development stage of the flower. Enzyme inhibition experiments targeting KETOACYL-ACP SYNTHASE activity alone in three chiloglottone-producing species (C. trapeziformis, C. valida, and C. aff. valida) significantly inhibited chiloglottone biosynthesis up to 88.4% compared to the controls. These findings highlight the role of coordinated (developmental stage- and tissue-dependent) fatty acid gene expression and enzyme activities for chiloglottone production in Chiloglottis orchids.


INTRODUCTION
Flowers have evolved diverse strategies to attract potential pollinators by conveying visual and olfactory signals (Raguso, 2004;Borghi et al., 2017). Some plants even use specific "semiochemical" signals that mimic the sex pheromone of female insects to sexually lure their highly specific male insect pollinators. This fascinating pollination strategy, commonly known as sexual deception (SD), has evolved multiple times on several continents (Gaskett, 2011). Amongst the Orchidaceace, hundreds of species spanning more than 20 genera across different subtribes are involved. Other cases are known in two other plant families (Bohman et al., 2016). In orchids, pollination by SD is often achieved when the male insect pollinator attempts to copulate with the flower, with volatiles holding the key to this interaction (Bohman et al., 2016;Wong et al., 2017b).
The Australian SD orchids of the genus Chiloglottis employ chiloglottones, a class of 2,5-dialkylcyclohexan-1,3-dione natural products (Peakall et al., 2010), to lure their specific male thynnine wasp pollinator (See Supplementary Figures 1A,B). For example, chiloglottone 1 is employed by Chiloglottis trapeziformis to attract its primary male pollinator, Neozeleboria cryptoides while the mostly allopatric C. valida employs chiloglottone 1 but attracts N. monticola for pollination. Conversely, a morphologically cryptic species with taxonomic affinity to C. valida (here after C. aff. valida secures pollination using a specific blend of chiloglottones 1 and 2 (Peakall et al., 2010;Peakall and Whitehead, 2014). In C. trapeziformis, chiloglottone 1 production has been shown to be specific to the densely clustered "insectiform" calli structure (callus) on the labellum, and its production is unexpectedly UV-B dependent (Falara et al., 2013). Flowers that open in the field become depleted of chiloglottone within 3-5 days when they are grown under light lacking in the UV-B range. Re-exposure to UV-B light, or natural sunlight, rapidly initiates chiloglottone production. Chiloglottone 1 is also dependent on the flowers' developmental stage (Amarasinghe et al., 2015). In all but the very young buds, chiloglottone 1 can be detected when flowers are exposed to UV-B or sunlight, with levels increasing from trace amounts from manually opened young buds to appreciable amounts in very mature buds and in fully open flowers.
The formation of chiloglottones has been predicted to involve the condensation of intermediates of varying chain lengths in the fatty acid (FA) biosynthetic and degradative (β-oxidation) pathways, followed by decarbonylation (Franke et al., 2009;Bohman et al., 2016). For example, the condensation of activated 3-oxohexanoic acid (CoA or ACP) and activated 2-hexenoic acid precursors and subsequent decarboxylation may give rise to chiloglottone 1 (See Supplementary Figure 1C). As such, formation of chiloglottone 2 may involve activated 3-oxohexanoic acid (CoA or ACP) and activated 2-octenoic acid precursors. In plants, such intermediates are ubiquitous, existing in the plastids as ACP derivatives or the peroxisomes as CoA derivatives. However, just how these precursors are made available for chiloglottone formation in an otherwise iterative pathway is unknown. One hypothesis for this may involve the mid-cycle termination (i.e., after a specified number of cycles) of FA biosynthesis and/or degradation (β-oxidation) pathway enzymes (Bohman et al., 2016).
Previously, we reported that tissue-specific floral transcriptome analysis, comparing the chiloglottone-emitting callus with the non-active labellum of C. trapeziformis, revealed that FA biosynthesis and β-oxidation pathways were highly coordinated in a tissue-specific manner whereby pathway transcripts are often highly or exclusively expressed in the callus (Wong et al., 2017a). This finding matched with the known tissue-specific distribution of chiloglottone 1 (Falara et al., 2013;Amarasinghe et al., 2015) and supported the current hypothesis of chiloglottones formation (Franke et al., 2009;Bohman et al., 2016). Interestingly, no induction of FA pathways was observed under UV-B treatment in the callus of mature buds and flowers despite the activation of UVR8-mediated signaling pathways, suggesting an unknown mechanism operating on chiloglottone 1 production at the transcriptional and metabolic level.
Building upon these previous biochemical, chemical, and transcriptional observations linked to chiloglottone formation, here we address three outstanding questions using targeted gene expression analysis and molecular inhibitor experiments. We ask the following three specific questions: (1) Are FA metabolic pathway genes developmentally regulated during the transition from young buds to mature flowers in C. trapeziformis? (2) Are FA metabolic pathway genes regulated in a tissue-specific manner in C. trapeziformis? (3) Does the inhibition of FA biosynthesis impair chiloglottone production in a cross section of Chiloglottis orchids?

Study Species
Whole Chiloglottis trapeziformis plants with single flowers at two developmental stages, namely very young buds (vyb) and naturally opened flowers (flw) (Amarasinghe et al., 2015), were sampled from a colony growing naturally within the Australian National Botanic Gardens (Canberra, ACT, Australia) in September 2014. Whole C. valida and C. aff. valida plants with single naturally opened flowers were also sourced from wild populations within the Kosciuszko National Park in NSW, Australia, in November 2014. See supplementary materials for plant growth conditions.

RNA Extraction, Library Construction, RNA Sequencing, and Transcriptome Analysis
Fifteen plants (three biological replicates with each replicate containing five individual plants) were used for each treatment. Floral tissues of vyb were carefully dissected to separate the stalked callus from the remaining labellum then immediately snap-frozen in liquid nitrogen. RNA extractions, library construction, and RNA sequencing was performed as previously described (Wong et al., 2017a;Xu et al., 2017). All raw sequence reads obtained in this study have been added to the existing BioProject accession PRJNA390683 and SRA study accession SRP109328 1 . Paired-end reads from (Wong et al., 2017a), and from non UV-B treated callus and labellum at very mature buds (vmb) and flw stages, were used in conjunction with sequenced tissue-specific reads at vyb stage obtained in this study and analyzed according to previously pipelines (Wong et al., 2017a).
The classifications of bud and flower stages used in this study follows Amarasinghe et al. (2015). The vyb stage is characterized by a small and very tightly closed buds with both green stalk and bud. The vmb stage is characterized by larger buds that are about to open with sepals and petals beginning to separate. Fully open flowers define the flw stage. See supplementary materials for detailed methods.

Pharmacological Inhibition Experiments and Chemical Analysis
All inhibition experiments were conducted on chiloglottonedepleted cut flowers from growth-chamber acclimatized whole plants (C. trapeziformis, C. valida, and C. aff. valida). Stock solutions (40 mM) of Cerulenin (Cayman Chemical, United States) were dissolved in Dimethyl Sulfoximine (Sigma, United States) and diluted to a final concentration of 100 µM with water (assay buffer). An assay buffer without the inhibitor served as the solvent control. The stalks of cut flowers were immersed into either 100 µM of Cerulenin (assay buffer) or solvent control. Next, the tip of each test tube was sealed with parafilm to ensure no direct contact between the solution and the flowers (3 flowers/tube) and then held in a test tube for 3 days. To induce chiloglottone production, UV-B treatments were conducted over a 2-h period on inhibitor-treated and control plants using a custom made light-box following Amarasinghe et al. (2015). The calli of all three species were immediately excised and assayed for chiloglottones as previously described (Falara et al., 2013). See supplementary materials for detailed chemical analysis methods.

Statistical Analysis
Differential expression analysis was performed using DESeq2 (Love et al., 2014). Transcripts differentially expressed between any given contrasts are defined as having an absolute log2 fold change (log2FC) >0.5 with a false discovery rate (FDR) <0.05. For all inhibition experiments, at least six flowers were used in each treatment or control group. The outcomes of inhibition treatments were analyzed with Student's t-test in R 2 .

NEW INSIGHTS INTO THE FORMATION OF INSECT SEX PHEROMONE-MIMICKING CHILOGLOTTONES
Spatio-Temporal Gene Expression of Fatty Acid Pathways Coincides With Chiloglottone 1 Production in C. trapeziformis In this study, we specifically addressed the first two objectives by interrogating the callus and labellum transcriptomes of vyb (obtained in this study; Supplementary Table 1) in conjunction with vmb and flw transcriptomes collected from 2 https://www.r-project.org/ the same colony in the same year obtained from our previous study (Wong et al., 2017a). Principal component analysis revealed that major differences among treatments was mainly driven by developmental stage followed by tissue specificity. (Supplementary Figure 2A). Enrichment (FDR <0.05) of gene ontology categories such as carbohydrate and lipid metabolic process were often observed in both developmental stage and tissue-specific contrasts ( Supplementary Figures 2A-C). Our findings revealed that many genes of the fatty acid (FA) biosynthesis and β-oxidation pathways were only upregulated in the callus during the transition from vyb to vmb. This coincides with a strong enrichment of lipid metabolism processes (FDR < 6.43 × 10 −11 ). Meanwhile, tissue-specific pairwise contrasts (callus vs. labellum) revealed that lipid metabolism process were also enriched in the callus regardless of developmental stage (FDR vmb < 3.50 × 10 −6 , FDR vyb < 1.81 × 10 −3 , FDR flw < 2.02 × 10 −2 ).

Three Fatty Acid Biosynthesis and Two β-Oxidation Pathway Genes Exhibit Callus-Specific Developmental Upregulation
Most of the FA biosynthesis such as one KETOACYL-ACP SYNTHASE III (CtKASIII), four KETOACYL-ACP REDUCTASE (CtKAR-L1 -3, mitochondrial CtKAR), two KETOACYL-ACP SYNTHASE I (CtKASI-1 and CtKASI-2), and one FATTY ACYL-ACP THIOESTERASE (CtFATB2) and FA β-oxidation pathway genes such as three ACYL-COA OXIDASE (CtACX2/3, CtACX4, and CtACX1/5) and one MULTIFUNCTIONAL PROTEIN (CtMFP3) were consistently upregulated in the callus compared to the labellum (Figures 1A,B). In addition, the transition from vyb to vmb stages generally involved coordinated upregulation of pathway genes, with transcripts remaining high in flw (no significant change from vmb to flw). Interestingly, four genes exhibited striking tissue (callus)-specific developmental upregulation (e.g., CtKASI-2, CtFATB, CtACX2/3, and CtACX4). These genes showed consistently (i) higher expression in the callus compared to the labellum and (ii) exhibited callus-specific developmental upregulation during the transition from vyb to vmb. We hypothesized that they may have direct involvement in chiloglottone 1 formation compared to those that are expressed at higher levels in the callus but are non-specifically upregulated in both callus and labellum tissues during the transitions from vyb to vmb (i.e., CtKASI-1, CtACX1/5, CtMFP1, CtMFP2, and CtMFP3) and vmb to flw (i.e., CtKAR-L1, CtMFP3, and CtKAT2L).

Inhibition of Fatty Acid Biosynthesis Significantly Block Chiloglottone Production in Chiloglottis Flowers
In light of the striking developmental stage-and tissues-specific upregulation of several key FA pathway genes (Figures 1A,B), we hypothesized that inhibition of their activities may affect chiloglottone 1 biosynthesis in C. trapeziformis. Therefore, we tested the effect of KAS inhibition using Cerulenin. Studies in several plant species have shown that Cerulenin specifically aff. valida. The number of flowers used in each treatment group are the following: CTR CER = 10, CTR Control = 9, CVA CER = 12, CVA Control = 8, CAV CER = 6, CAV Control = 7. Bars represent ± SE. Asterisks indicate significant differences between treatments at * P < 0.05, * * P < 0.01, * * * P < 0.001 based on Student's t-test. Light and dark blue colors indicate chiloglottone 1 and 2, respectively. All images have been reproduced with permission from the respective copyright holders. Please refer to the Section "Acknowledgments" for image credits.

PERSPECTIVES ON THE BIOSYNTHESIS OF INSECT SEX PHEROMONE-MIMICKING CHILOGLOTTONES
Parallel temporal changes in enzyme activities, protein content, and their corresponding structural gene expression are often pivotal for developmentally regulated and/or tissue-specific volatile production in flowers (Dudareva et al., 2013). Confirmed cases where the biosynthesis of the semiochemicals involved in the sexual mimicry is known in SD orchids is presently limited to the biosynthesis of 7-, 9-, and 12-alkenes in Ophrys orchids (Schlüter et al., 2011;Xu et al., 2012;Sedeek et al., 2016) and (S)-β-citronellol in Caladenia plicata (Xu et al., 2017). In this study, a strategic developmental stage and tissue differential expression analysis of C. trapeziformis floral transcriptomes identified two key transcriptional trends linked to the distribution of chiloglottone in the flowers: (i) During the transition from vyb to vmb, large suites of FA biosynthesis and β-oxidation genes are upregulated. (ii) In all development stages tested (especially vmb), FA pathway genes were consistently induced in the callus compared to the labellum (Figure 1). Despite this strong coordinated regulation, no genes from other major biosynthetic pathway (e.g., glycolysis and tricarboxylic acid) showed similar preferential expression (Supplementary Tables 2, 3). Our findings lend further support to the existing hypothesis of chiloglottone formation in planta via a FA biosynthetic route and demonstrate a key role of floral developmental transitions (i.e., vyb to vmb) in the priming of FA pathway gene expression to initiate chiloglottone biosynthesis.
While in plants there is no evidence that Cerulenin directly/indirectly inhibits FA β-oxidation (Shimakata and Stumpf, 1982;Dehesh et al., 1998;Yasuno et al., 2004), we cannot yet rule out the possibility that FA β-oxidation as an alternative route for chiloglottone biosynthesis (Bohman et al., 2016). This is because any inhibition of KAS enzymes by Cerulenin may have an indirect effect -reducing the levels of de novo synthesized FAs (e.g., C16/C18) -and thus restricting further catabolism (via β-oxidation) to yield alternative activated (CoA) FA precursors of appropriate chain lengths for chiloglottone biosynthesis. Transcriptome observations showing coordinated (developmental stage-and tissue-dependent) expression of several ACX in C. trapeziformis also support the possibility that one or both putative chiloglottone precursors could be obtained by FA β-oxidation. Nonetheless, our results reinforce the current hypothesis that FA biosynthesis serves as the crucial "starting point" for chiloglottones formation. Future studies will require targeted knockdown of pathway candidates and metabolomics analysis for activated (ACP/CoA) FA precursors to ascertain whether chiloglottone production is largely determined by precursor availability (via the FA biosynthesis and/or β-oxidation pathway) or during the condensation of activated precursors that may involve novel KAS activities. Together, these new findings highlight the role of coordinated (developmental stageand tissue-dependent) FA gene expression and enzyme activities for chiloglottone production in our study species and may have widespread implications for Chiloglottis and other orchid genera employing chiloglottones for SD pollination (Peakall et al., 2010).

AUTHOR CONTRIBUTIONS
DW and RA performed the experiments and analyzed the data. RP and EP secured funding, designed the study, and coordinated the experiments and data analysis. DW wrote the article with assistance from EP and RP. All authors have read and approved the paper.

FUNDING
This work was supported by Australian Research Council projects DP1094453 and DP150102762 to RP and EP.