Background {#Sec1}
==========

Flowers emit a large group of volatile organic compounds (VOCs) that play crucial roles in interactions with other organisms. In general, floral VOCs are considered not only to serve to attract pollinators to ensure successful reproduction, but also to act as defense agents against microbes and herbivores \[[@CR1]\]. Most of floral VOCs that are synthesized on the petals of the flowers fall into three classes of VOCs namely terpenoids, phenylpropanoids/benzenoids and volatile esters \[[@CR2]\]. The quality and quantity of floral VOCs emitted determine floral scents and are also vital in determining the economic value of flowering plants as well as its usage in the flavors and fragrances industry.

Terpenoids, the largest class of VOCs, are produced and emitted by a number of floral species such as snapdragon and ylang ylang \[[@CR3], [@CR4]\]. In plants, biosynthesis of terpenoids involves three steps: (1) production of C~5~ isoprenoid precursors isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) by two compartmentally separated pathways: the methylerythritol 4-phosphate (MEP) pathway in plastids and the mevalonate (MVA) pathway in the cytosol, (2) condensation of IPP and DMAPP into geranyl diphosphate (GPP, C~10~), farnesyl diphosphate (FPP, C~15~), and geranylgeranyl diphosphate (GGPP, C~20~), the precursors for the production of mono-, sesqui- and diterpenes, respectively and (3) generation of diverse terpene structures by the final cyclization and oxidation steps carried out by the terpene synthases (TPS) and cytochrome P450s (CYP450) \[[@CR5]\]. A large number of TPSs involved in floral scent production have been characterized in a variety of plants \[[@CR1]\].

Volatile esters are also known to contribute to the aroma of many fruits and flowers. For instance, volatile esters such as methyl hexanoate, ethyl hexanoate, methyl jasmonate and methyl tiglate contribute profoundly to the total floral scent composition of many plant species \[[@CR2]\] whereas hexyl acetate, ethyl 2-methylbutanoate and methyl 3-methylvalerate have been reported to be key flavor constituents of many fruits including apple \[[@CR6]\] and strawberry \[[@CR7]\]. They can be derived from fatty acids or from branched-chain and aromatic amino acids \[[@CR8], [@CR9]\]. Conversion of fatty acids to volatile esters occur via three processes: α-oxidation, β-oxidation and lipoxygenase (LOX) pathway \[[@CR10]\]. On the other hand, the first step in the conversion of amino acids to volatile esters is a deamination step catalyzed by aminotransferases resulting in α-ketoacids. α-ketoacids can be converted into (1) carboxylic acids via oxidative decarboxylation, (2) aldehydes via decarboxylation and (3) α-hydroxyacids via reduction by a multi-subunit complex of enzymes named as α-ketoacid dehydrogenase complex \[[@CR11]\]. The last step in the formation of volatile esters is the esterification of the aldehydes, acids and alcohols formed by the above-mentioned processes by a class of enzymes called alcohol acyl transferases (AAT) \[[@CR12]\].

The initial substrate for the biosynthesis of phenylpropanoids and benzenoids is phenylalanine, which is provided by the shikimate pathway \[[@CR13]\]. Production and emission of phenylpropanoids/benzenoids are highly regulated both temporally and spatially \[[@CR14]\]. Many phenylpropanoids/benzenoids biosynthetic enzymes have been identified and functionally characterized \[[@CR14], [@CR15]\]. However, the regulatory mechanisms involved in their production are relatively unknown as compared to terpenoids biosynthesis.

*Magnolia champaca*, commonly known as champak is a tall evergreen aromatic tree of the Magnoliaceae family. Since various tissues of champak have been reported to possess a broad range of medicinal properties, they are traditionally used in the treatment of various diseases such as cephalalgia, ophthalmia, gout, rheumatism and microbial infections \[[@CR16]\]. In addition, the tree has commercial value for its strongly diffusive fragrant flowers which are widely used in cosmetic industries for the production of perfumes and essential oils due to the high amounts of VOCs \[[@CR17]\]. The composition of champak essential oil has been investigated previously showing that it varied widely due to various factors including extraction method and flower maturity \[[@CR18]--[@CR20]\]. Although the reported percentages of the compounds were different, champak flowers mainly produced linalool and linalool oxides, 2-phenyl ethanol, methyl anthranilate, indole and methyl linoleate \[[@CR18]--[@CR20]\].

Despite the overall popularity of the champak, to our knowledge very little DNA sequences are available on GenBank for this species. The availability of DNA sequences may provide an opportunity to identify genes involved in biosynthesis of floral VOCs. Identifying and investigating the regulatory mechanisms behind floral VOCs biosynthesis is essential to understand their roles in plants and to improve VOC production through metabolic engineering. Here, we present the chemical composition and the transcriptome data from fully opened champak flowers, which were comprehensive enough to discover and analyze major secondary metabolite pathways associated with floral scent production. In addition, we identified the genes encoding the first and the last enzymes of the MEP pathway, DXS and HDR along with the (R)-linalool synthase, which makes substantial contribution to the fragrance of champak essential oil.

Results {#Sec2}
=======

Chemical composition profile of champak flowers {#Sec3}
-----------------------------------------------

A previous report on the floral VOCs of *Michelia alba*, one of the species in *Magnoliaceae* family, showed the characteristic *M. alba* fragrance to be intense from half bloom flowers to fully opened flowers \[[@CR19]\]. To investigate the VOCs profile of fully opened champak flowers grown in Singapore, we analyzed the headspace volatiles emitted from champak flowers by GC-MS. It was observed that the VOCs emitted from champak flowers consisted of approximately 43 chemical compounds, which can contribute towards the floral scent (Fig. [1a](#Fig1){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}). These diverse VOCs emitted from champak flowers can be broadly classified into three classes: terpenoids (46.9%), volatile esters (38.9%) and phenylpropanoids/benzenoids (5.2%) (Fig. [1b](#Fig1){ref-type="fig"}). Out of 46.9% of terpenoids identified, 25.8% were monoterpenes and their derivatives, whereas 21.1% were sesquiterpenes. Five monoterpenoids, β-myrcene, β-ocimene, β-linalool, neo-allo-ocimene and α-terpineol were found in champak VOCs profile. Of these, β-linalool and its oxides such as *trans*-furanoid linalool oxide, *cis*-furanoid linalool oxide and *trans*-pyranoid linalool oxide were the most abundant components comprising about 17.5% of the total volatile composition and were followed by *cis*- and *trans*-forms of β-ocimene (6.8%). These dominant monoterpenes may contribute significantly to floral scent of champak. Although we found 15 sesquiterpenes accounting for \~21% of total volatiles, many of them were emitted in minor quantities (\< 1%) except β-cubebene (8.6%), δ-elemene (3.2%), isogermacrene-D (2.3%), β-copaene (1.5%) and β-bourbonene (1.4%).Fig. 1VOCs emitted from champak flowers. **a** Gas chromatogram of champak floral VOCs. The peak numbers counted in the gas chromatogram are identical to the compounds listed in Table [1](#Tab1){ref-type="table"}. IS, internal standard, camphor. **b** Relative abundance of different classes of champak floral VOCs. Data used in the percentage calculation were obtained from the average of three replicates. Error bars indicate standard deviation (SD; *n* = 3) Table 1Volatile organic compounds composition emitted from champak flowersNo.^a^CompoundRT (min)^b^RI^c^Formulang.g^−1^ flower.h^-1\ d^RA (%)1Ethyl propanoate6.492703C~5~H~10~O~2~1.96 ± 0.130.5622-Methylbutanoate6.593713C~5~H~10~O~2~4.14 ± 0.401.193Methyl 2-methylbutanoate7.188744C~6~H~12~O~2~80.69 ± 4.3423.164Ethyl butanoate7.383792C~6~H~12~O~2~1.93 ± 0.400.555Ethyl 2-methylbutanoate8.048841C~7~H~14~O~2~32.66 ± 3.649.3863-Methylbutanal oxime8.189851C~5~H~11~NO0.57 ± 0.010.167Methyl tiglate8.297859C~6~H~10~O~2~0.43 ± 0.010.128Methyl hexanoate9.122906C~7~H~14~O~2~5.17 ± 0.571.489Propyl 2-methylbutanoate9.477937C~8~H~16~O~2~1.44 ± 0.020.4110β-Myrcene10.261982C~10~H~16~0.84 ± 0.010.2411Ethyl hexanoate10.338989C~8~H~16~O~2~2.16 ± 0.530.6212(Z)-β-Ocimene11.0451031C~10~H~16~7.71 ± 0.782.2113(E)-β-Ocimene11.251043C~10~H~16~16.10 ± 1.424.6214(Z)-Furanoid linalool oxide11.7631063C~10~H~18~O~2~2.91 ± 0.10.8415(E)-Furanoid linalool oxide12.0471079C~10~H~18~O~2~27.07 ± 0.927.7716β-Linalool12.1681092C~10~H~18~O24.71 ± 1.547.0917Phenylethyl alcohol12.4971103C~8~H~10~O17.12 ± 1.834.9118Neo-allo-ocimene12.6561125C~10~H~16~1.91 ± 0.260.5519Benzyl nitrile12.9121140C~8~H~7~N0.68 ± 0.140.1920(E)-Pyranoid linalool oxide13.4611173C~10~H~18~O~2~6.34 ± 0.291.8221α-Terpineol13.8091185C~10~H~18~O0.61 ± 0.040.1722Geranyl methyl ether14.2251219C~11~H~20~O0.99 ± 0.030.2823Geranyl ethyl ether 115.2251264C~12~H~22~O0.73 ± 0.020.2124Indole15.5191285C~8~H~7~N30.03 ± 2.818.6225δ-Elemene16.1941342C~15~H~24~11.13 ± 0.713.1926Methyl anthranilate16.251346C~8~H~9~NO~2~0.54 ± 0.010.1527α-Copaene16.3821355C~15~H~24~0.41 ± 0.010.1228α-Ylangene16.7711380C~15~H~24~2.50 ± 0.230.7229β-Bourbonene17.0211396C~15~H~24~4.89 ± 0.301.4030Isogermacrene-D17.5181430C~15~H~24~7.99 ± 0.322.2931β-Copaene17.671441C~15~H~24~5.38 ± 0.261.5532γ-Amorphene17.8981456C~15~H~24~1.72 ± 0.140.4933γ-Muurolene18.0481466C~15~H~24~2.28 ± 0.180.6534(Z)-Muurola-4(14),5-diene18.1721475C~15~H~24~1.73 ± 0.070.5035β-Cubebene18.4431478C~15~H~24~29.87 ± 3.428.5736δ-Guaiene18.5421500C~15~H~24~0.84 ± 0.050.2437α-Muurolene18.6221506C~15~H~24~1.48 ± 0.040.4338δ-Cadinene18.7171513C~15~H~24~1.3 ± 0.060.3739γ-Cadinene18.8621523C~15~H~24~0.99 ± 0.030.2840β-Cadinene18.9391529C~15~H~24~1.22 ± 0.020.3541Methyl palmitate23.8361908C~17~H~34~O~2~0.58 ± 0.060.1742Methyl linoleate25.7922080C~19~H~34~O~2~2.83 ± 0.300.8143Methyl linolenate25.8722087C~19~H~32~O~2~1.83 ± 0.170.52^a^Compounds are listed in the same order as their elution from HP-5MS column^b^Retention Time in minutes^c^Retention Indices calculated using C~7~-C~30~ n-alkanes on HP-5MS column^d^Calculated from three independent readings

The pleasant fragrance of flowers such as *Wisteria*, orchids and lilies is often due to the presence of a particular volatile ester or a mixture of several volatile esters \[[@CR2]\]. We found volatile esters as the second major class of champak floral volatiles, which constituted 38.9% of total VOCs (Fig. [1](#Fig1){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}). Particularly, different derivatives of butanoic acid contributed significantly (35.9%) to the total VOCs of champak flowers. Of these, methyl 2-methylbutanoate solely constituted 23.2% of total floral volatiles. In addition, ester forms of hexanoic acid, propanoic acid, linoleic acid and linolenic acid were also among the champak floral VOCs. Phenylethyl alcohol (4.9%) was the only phenylpropanoid identified and minor amounts of two benzenoids, benzyl nitrile and methyl anthranilate were also part of the champak floral volatile composition.

VOCs composition of fully opened champak flowers was analyzed at several time points of the year. However, we could not see significant variation in volatile profile of these flowers, as the weather conditions in Singapore do not fluctuate considerably throughout the year due to the geographical location of Singapore.

RNA-seq, de novo assembly, and annotation of transcriptome {#Sec4}
----------------------------------------------------------

To identify genes involved in the biosynthesis of floral VOCs from champak flowers, we sequenced RNA libraries derived from flowers using an Illumina HiSeq2000. Illumina sequencing from RNA-seq libraries yielded approximately 169.6 million reads of 101 bp length. After eliminating the adaptor sequences, 122.8 million reads with 47% GC content were obtained. FastQC \[[@CR21]\] analysis showed that 73.2% of the total sequences were of quality above Q30 (Additional file [1](#MOESM1){ref-type="media"}: Figure S1a). As the genomic sequence of champak is unavailable, the reads were assembled de novo into a total of 47,688 unigenes with N50 of 1814 bp using Trinity method (\[[@CR22]\]; Additional file [2](#MOESM2){ref-type="media"}; Table [2](#Tab2){ref-type="table"}). Total number of assembled unigenes might be overestimated due to the absence of reference genome. It should be noted that more than 5000 unigenes were observed in the length range of 200--300 bp, most of which were not annotated by public databases. A large number of unannotated unigenes may indicate (1) novel genes specifically expressed in champak flowers, (2) incomplete sequencing of some very low expressed genes yielding two/more unigenes, (3) presence of non-coding RNAs \[[@CR23], [@CR24]\].Table 2Overview of RNA-seq de novo assembly resultsNo. ReadsNo. UnigenesN50 (bp)No. Annotated% Annotation169,602,59447,688181427,04356.7

Of the 47,688 unigenes, 27,043 (56.7%) were annotated using BLASTX (E-value cut-off of 10^−5^) against four protein databases: National Centre for Biotechnology Information (NCBI) non-redundant (nr), *Arabidopsis* (*Arabidopsis thaliana*), rice (*Oryza sativa*) and grape (*Vitis vinifera*). BLASTX analysis showed that 41% of unigenes had extremely high homology (E-value \<10^−100^), 20% had high homology (10^−100^ \< E-value \<10^−50^) and 39% had moderate homology (10^−50^ \< E-value \<10^−5^) to nr database (Fig. [2a](#Fig2){ref-type="fig"}). The species distribution of the best-match result showed that the top hits for 58% of unigenes were from *Vitis vinifera*, followed by *Ricinus communis* (9%), *Populus trichocarpa* (9%), *Glycine max* (5%), *Oryza sativa* (5%) and *Arabidopsis* (1.1%) (Fig. [2b](#Fig2){ref-type="fig"}).Fig. 2Analysis of champak unigenes. **a** E-value distribution of top BLASTX hits. **b** Species distribution based on the top BLASTX hits. **c** Top 10 GO terms from all three categories

Functional classification of genes {#Sec5}
----------------------------------

Champak unigenes were functionally classified into different Gene Ontology (GO) terms using Trinotate \[[@CR25]\]. Classification showed that 38% of the annotated genes were involved in biological process, 29.5% in cellular component and 32.5% in molecular function (Additional file [2](#MOESM2){ref-type="media"}). We further looked at the top ten GO-terms in each of the three GO categories (Fig. [2c](#Fig2){ref-type="fig"}). Within biological process, transcription (19.6%) and regulation of transcription (14.2%) were the two dominating GO terms and they were followed by defense response, which suggests that champak flowers are probably an active tissue for secondary metabolism. Among the category cellular component, 51.6%, 50%, 26.9%, 25.6% of the annotated genes were classified into the GO-terms nucleus, integral component of membrane, plasma membrane and cytoplasm respectively. In the group of molecular function, ATP binding, metal ion binding, DNA binding and zinc binding were the principal GO-terms of molecular function comprising of 50.5%, 23.7%, 20.2% and 17.8% annotated genes respectively.

Analysis of highly expressed genes in champak flowers {#Sec6}
-----------------------------------------------------

We investigated the top 20 transcripts that were highly expressed in champak flowers and found many of them to be involved in the biosynthesis of floral VOCs (Fig. [3a](#Fig3){ref-type="fig"}). For instance, transcript encoding for geranyl diphosphate synthase small subunit (GPS.SSU), S-adenosyl-L-methionine-dependent methyltransferases (SAMT), methionine synthase, S-adenosylmethionine synthase, 13S--lipoxygenase, medium chain fatty-acid CoA ligase and alcohol acyl transferase were all highly expressed and are known to be involved in the biosynthesis of terpenoids, benzenoids and volatile esters. This expression pattern correlates well with the result showing VOCs profile emitted from champak flowers in Fig. [1](#Fig1){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}. Many fatty acid-derived methyl esters (FAMEs) are believed to be synthesized from polyunsaturated fatty acids via the lipoxygenase pathway \[[@CR26]\]. Several FAMEs such as methyl palmitate, methyl linoleate and methyl linolenate were found in champak VOCs (Fig. [1](#Fig1){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}). The highly expressed unigene encoding 13S--lipoxygenase might be responsible for their synthesis.Fig. 3Expression and pathway analysis of champak unigenes. **a** Top 20 genes most highly expressed in champak flowers. **b** Top 15 mapped pathways annotated by Plant metabolic pathway database (PlantCyc, \[[@CR30]\])

The second set of abundant transcripts encoded for proteins that function in plant defense. These were the metallothionein (MT), methyl salicylate esterase, stress-related protein and serine hydroxymethyltransferase (SHMT), which play a role in detoxification of heavy metals and reactive oxygen species \[[@CR27]\], to confer immunity to plants \[[@CR28]\] and in controlling cell damage from stress conditions \[[@CR29]\], respectively.

Pathway analysis of unigenes expressed in champak flowers {#Sec7}
---------------------------------------------------------

The top fifteen pathways with highest average Fragments Per Kilobase of transcript per Million mapped reads (FPKM) of tagged enzymes were identified using the plant metabolic pathway database (PlantCyc, \[[@CR30]\]; Additional file [2](#MOESM2){ref-type="media"}). Figure [3b](#Fig3){ref-type="fig"} shows that majority of the top fifteen pathways were involved in secondary metabolites biosynthesis. It was apparent that pathways involved in the production of substrates (GPP, FPP, and GGPP) essential for terpenoids biosynthesis were highly expressed. In addition, genes involved in pathways leading to the biosynthesis of S-adenosyl-L-methionine, a common substrate for methyl group transfers, and its degradation were also highly expressed in the champak flowers. These results are not surprising as champak flowers produced diverse mono- and sesqui-terpenes as well as methylated volatile esters (Fig. [3b](#Fig3){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}). Interestingly, pathways for the indole-3-acetate (IAA) and indole-3-acetyl-ester conjugate biosynthesis were among the top 20. This suggests that auxin signaling probably impedes floral VOCs biosynthesis as IAA conjugates are known to be endogenous auxin inhibitors \[[@CR31]\].

The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database was used to identify the biological pathways active in champak flowers \[[@CR32]\]. When the unigenes were assigned KEGG Orthology (KO) numbers using the bi-directional best-hit method in KEGG Automatic Annotation Server (KAAS) \[[@CR33]\], a total of 6597 unigenes were mapped onto 373 KEGG pathways (Additional file [2](#MOESM2){ref-type="media"}). Carbon metabolism (ko01200), biosynthesis of amino acids (ko01230) and protein processing in endoplasmic reticulum (ko04141) were identified as the top three KEGG pathways (Additional file [1](#MOESM1){ref-type="media"}: Figure S1b). Among the top twenty KEGG pathways, terpenoid backbone biosynthesis, pentose phosphate pathway and glycolysis were found, supporting the results obtained by the plant metabolic pathway database (Fig. [3b](#Fig3){ref-type="fig"} and Additional file [1](#MOESM1){ref-type="media"}: Figure S1b). Additionally, pathway analysis using KEGG database identified 2-oxocarboxylic acid metabolism (ko01210) as one of the top pathways. This pathway includes the conversion of L-isoleucine to (S)-3-methyl-2-oxopentanoic acid, which is the first step in the proposed biosynthetic pathway of 2-methylbutanoate, methyl 2-methylbutanoate, ethyl 2-methylbutanoate and propyl 2-methylbutanoate (Additional file [1](#MOESM1){ref-type="media"}: Figure S1b and Figure S2).

Next, we investigated the transcripts encoding enzymes for the precursor pathways related to champak floral VOCs biosynthesis. Figure [4a](#Fig4){ref-type="fig"} shows the expression profile of genes involved in MEP, MVA, shikimate, phenylpropanoids/benzenoids, lipoxygenase and amino acid catabolic pathways, which are main pathways for champak floral VOCs such as terpenoids, phenylpropanoids/benzenoids and volatile esters. Most of the genes involved in these pathways showed high abundance in the RNA-seq of champak flowers and included full length open reading frames (ORFs) with high homology to their orthologous genes from other plants (Additional file [1](#MOESM1){ref-type="media"}: Figure S3).Fig. 4Analysis of expression levels of genes involved in floral VOCs biosynthesis. **a** Number of transcripts for individual genes involved in different pathways. Transcript counts are presented in a logarithmic scale of 10. **b** qRT-PCR analysis of representative genes from (i) MEP (ii) MVA (iii) Shikimate (iv) Phenylpropanoid/Benzenoid (v) Lipoxygenase and (vi) Amino acid catabolic pathways in leaves (L) and flowers (F). DXR, 1-deoxy-D-xylulose 5-phosphate reductoisomerase; CMS, 2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase; CMK, 4-(cytidine 5′-diphospho)-2-C-methyl-D-erythritol kinase; MCS, 2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase; HDS, 4- hydroxy-3-methylbut-2-en-1-yl diphosphate synthase; AACT, acetyl-CoA acetyltransferase; HMGS, hydroxymethylglutaryl-CoA synthase; HMGR, hydroxymethylglutaryl-CoA reductase; MVK, mevalonate kinase; PMK, phosphomevalonate kinase; MPDC, mevalonate diphosphate decarboxylase; DAHPS, 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase; DHQS, 3-dehydroquinate synthase; SDH, shikimate dehydrogenase; SK, shikimate kinase; CM, chorismate mutase; CS, chorismate synthase; PAL, phenylalanine ammonia lyase; C4H, cinnamate-4-hydroxylase; C3H, p-coumarate-3-hydroxylase; COMT, caffeic acid/5-hydroxyferulic acid O-methyltransferase; 4CL, 4-coumaroyl-CoA ligase; CCoAOMT, caffeoyl-CoA 3-O-methyltransferase; CCR, cinnamoyl-CoA reductase; CAD, cinnamyl alcohol dehydrogenase; 9-LOX, linoleate 9S--lipoxygenase 5; 13-LOX, linoleate 13S--lipoxygenase; HPL, hydroperoxide lyase; AOS, allene oxide synthase; AOC, allene oxide cyclase; OPR, 12-oxophytodienoate reductase; BCAT, branched-chain amino-acid aminotransferase; BCKDC_E1, pyruvate dehydrogenase E1 alpha; BCKDC_E2, dihydrolipoamide acyltransferase; BCKDC_E3, dihydrolipoamide dehydrogenase; AAT, alcohol acyl transferase

The expression levels of these transcripts were checked in champak flowers and leaves by quantitative real-time PCR (qRT-PCR) in order to see if they are specifically up-regulated in flowers, thus contributing to floral scent (Fig. [4b](#Fig4){ref-type="fig"}). Overall, approximately 70% of these transcripts analyzed showed at least 2-fold higher expression in flowers than leaves. Interestingly, transcripts encoding branched-chain amino acid aminotransferase (BCAT) and alcohol acyl transferase (AAT), which are involved in the biosynthesis of several volatile esters including methyl 2-methylbutanoate (Additional file [1](#MOESM1){ref-type="media"}: Figure S2) showed preferential expression in flowers than in leaves.

Characterization of MEP pathway genes {#Sec8}
-------------------------------------

Since monoterpenes are one of most abundant VOCs in champak flowers, we looked at genes encoding MEP pathway enzymes that are necessary to provide IPP and DMAPP for monoterpene biosynthesis in the plastids \[[@CR34]\]. Based on the search of the champak transcriptome against orthologous sequences from other plant species, we found single copy genes of all MEP pathway enzymes except *DXS* (1-deoxy-D-xylulose 5-phosphate synthase) that encodes four isoforms in champak flower (Additional file [1](#MOESM1){ref-type="media"}: Table S1).

We further investigated the champak *DXS* family genes and the *HDR* gene, as they are the first and last enzyme in MEP pathway, respectively. Four isoforms of DXS were found in champak flowers, which were designated *McDXS1--4* based on transcript abundance. Phylogenetic analysis based on deduced amino acid sequences of four *McDXS* cDNAs showed that McDXS1 and McDXS2 belonged to clade II and their transcripts were predominantly expressed in champak flowers than in leaves (Fig. [5a](#Fig5){ref-type="fig"} and Additional file [1](#MOESM1){ref-type="media"}: Figure S4a). Clade II DXS is involved in secondary metabolism. On the other hand, McDXS4 was the only member of clade I and McDXS3, the most divergent among the four McDXSs, was part of the most distant group - clade III. The amino acid residues essential for binding of glyceraldehyde-3-phosphate (G3P) and thiamine pyrophosphate (TPP) were conserved in McDXS1, McDXS2 and McDXS4, but not in McDXS3 (Additional file [1](#MOESM1){ref-type="media"}: Figure S4b).Fig. 5Characterization of *McDXSs* and *McHDR*. **a** Expression of champak 1-deoxy-D-xylulose 5- phosphate synthase (*McDXSs*) and 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase (*McHDR*) genes in leaves (L) and flowers (F). **b** Subcellular localization of McDXSs and McHDR. Auto, chlorophyll autofluorescence; YFP, YFP channel image; Merged, merged image of Auto and YFP; Light, light microscopy image. Bars, 20 μm. **c** and **d** Complementation assays of champak DXS isoforms (**c**) and HDR (**d**) using *E.coli dxs* ^*−*^ and MG1655 *ara\<\>ispH* mutants, respectively. For DXS complementation assay, *E.coli* cells containing *Arabidopsis thaliana* DXS (AtDXS1) and empty vector (EV) were used as a positive and negative control, respectively. For the HDR complementation assay, *Arabidopsis thaliana* HDR (AtHDR) was used as a positive control while empty vector (EV) and McDXS1 were used as negative controls

Analysis of deduced amino acid sequence of *McHDR* showed that four-cysteine residues, which might be essential for the coordination of iron-sulfur bridge were well conserved (\[[@CR35]\]; Additional file [1](#MOESM1){ref-type="media"}: Figure S5).

All four DXSs and HDR contained a transit peptide (TP) sequence at their N-terminus for plastidic targeting (Additional file [1](#MOESM1){ref-type="media"}: Table S1). In order to validate this result, we observed the subcellular localization of McDXS1--4 and McHDR. Yellow fluorescent protein (YFP) fused DXS or HDR construct was transiently expressed in *N. benthamiana* leaves using an *Agrobacterium*-mediated infiltration. Figure [5b](#Fig5){ref-type="fig"} shows that all McDXSs as well as McHDR were clearly localized in chloroplasts of *N. benthamiana* leaves.

To find out if all the four McDXSs are functionally active, we carried out complementation assay using a *dxs*-deficient *Escherichia coli* strain. The *dxs* ^*−*^ *E. coli* cells can grow normally in the presence of mevalonate but on mevalonate-free media, the cells require functional DXS for their viability \[[@CR36]\]. Figure [5c](#Fig5){ref-type="fig"} shows that transformants harbouring McDXS1 and McDXS2 were able to complement the *dxs* ^*−*^ lines on mevalonate-free media similar to the positive control *Arabidopsis* DXS1. However, transformants containing McDXS3, MCDXS4 and empty vector as a negative control were unable to grow on mevalonate-free media even after 2d of incubation.

The enzymatic activity of McHDR was analyzed by complementation assay using a *hdr*-deficient *E. coli* strain, *MG1655 ara\<\>ispH* that grows normally in the presence of arabinose, but require a functional exogenous HDR to survive in the presence of glucose \[[@CR37]\]. Figure [5d](#Fig5){ref-type="fig"} shows that *MG1655 ara\<\>ispH* cells transformed with *McHDR* were capable of growing well on media containing glucose, similar to the cells transformed with the positive control *Arabidopsis HDR.* But, the cells transformed with either the empty vector or *McDXS1* did not grow on media containing glucose. This result demonstrated that the *McHDR* encodes a functionally active HDR enzyme.

Terpene synthases from champak flowers {#Sec9}
--------------------------------------

Champak flowers emit at least twenty six different types of mono- and sesqui-terpenes and their derivatives (Fig. [1a](#Fig1){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}). We found approximately nine candidate *TPS* unigenes from the transcirptome data of champak flowers (Additional file [1](#MOESM1){ref-type="media"}: Table S2). While most candidate *TPS* transcripts were partial mRNA sequences, only one transcript contained a full-length ORF. Deduced amino acid sequence of full-length TPS showed 67% similarity with *trans*-ocimene synthase from mountain pepper (*Litsea cubeba*)*.* Interestingly, we found four partial mRNA sequences that showed high homology up to 98% to the single gene encoding β-cubebene synthase (Mg25) from *Magnolia grandiflora* \[[@CR38]\]. Since β-cubebene was the most abundant sesqui-terpene found in champak flowers accounting for 8.57% of the total VOCs, 18.2% of terpenes, and 40.5% of sesqui-terpenes (Fig. [1a](#Fig1){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}), they were predicted to encode β-cubebene synthase.

Next, we further investigated the *TPS* gene (designated here as *McTPS1*) containing the complete ORF. *McTPS1* was predominantly expressed in flowers than in leaves (Fig. [6a](#Fig6){ref-type="fig"}). According to the phylogenetic analysis, McTPS1 was grouped into TPS-b subfamily, which commonly represents mono-TPSs (Fig. [6b](#Fig6){ref-type="fig"}; \[[@CR39]\]). In addition to the aspartate-rich domains, DDXXD and NSE/DTE motif, which is highly conserved in plant TPSs, McTPS1 contained two distinct structural domains of the TPS-b group, a TP sequence for plastid targeting and the R(R)X~8~W motif for monoterpene cyclization at N-terminal region (Additional file [1](#MOESM1){ref-type="media"}: Figure S6). As expected, the McTPS1-YFP was localized in chloroplast of *N. benthamiana* cells (Fig. [6c](#Fig6){ref-type="fig"}), suggesting McTPS1 is a mono-TPS.Fig. 6Characterization of *McTPS1*. **a** Expression of *McTPS1* in leaves (L) and flowers (F). **b** Phylogenetic analysis of McTPS1. A maximum likelihood tree was constructed using amino acid sequences of McTPS1 with other plant TPSs. Accession numbers of proteins used in phylogenetic analysis are listed in Additional file [1](#MOESM1){ref-type="media"}: Table S3. **c** Subcellular localization of McTPS1. Auto, chlorophyll autofluorescence; YFP, YFP channel image; Merged, merged image of Auto and YFP; Light, light microscopy image. Bars, 20 μm

To figure out the exact function of McTPS1, 6His-tagged McTPS1 recombinant protein was used for in vitro assays in the presence of geranyl diphosphate (GPP) and farnesyl diphosphate (FPP) as the common substrate for mono- and sesqui-TPS, respectively. Figure [7a](#Fig7){ref-type="fig"} shows that McTPS1 reacted with GPP to produce (R)-linalool, a monoterpene alcohol, whereas it did not react with FPP. As a negative control, heat denatured McTPS1 protein was added to GPP and FPP reaction mixtures and they failed to produce any terpene. Therefore, our in vitro assays identified McTPS1 as the (R)-linalool synthase, responsible for the production of the most abundant monoterpene of champak flowers. In vivo functional studies were carried out by infiltrating the leaves of *N. benthamiana* with *Agrobacterium* harboring *McTPS1*. Plants infiltrated with green fluorescent protein (GFP) served as a negative control. At 3 days post-infiltration (dpi), the volatiles were collected from the plants by push-pull headspace method and analyzed by GC-MS. *N. benthamiana* plants expressing McTPS1 emitted (R)-linalool which is consistent with those obtained in vitro whereas no linalool was detected in *N. benthamiana* plants expressing GFP (Fig. [7b](#Fig7){ref-type="fig"}). When chiral column was used to analyze the volatile composition of champak flowers, it was found that champak flowers emitted (R)-linalool almost exclusively (Additional file [1](#MOESM1){ref-type="media"}: Figure S7).Fig. 7Functional characterization of McTPS1. **a** In vitro enzymatic assay of recombinant 6His-McTPS1 protein using either GPP or FPP as substrate. Control, TPS assay of the heat-inactivated 6His-McTPS1 recombinant protein with GPP. The reaction products were analyzed by chiral phase GC-MS. The peak obtained from in vitro assay was identified with the authentic standard, (±)-linalool. Mass spectra of the peak from the assay and (R)-linalool are shown at the bottom of chromatograms. **b** In vivo characterization of McTPS1 by transient expression of McTPS1 in *N. benthamiana* leaves. VOCs were collected from McTPS1 or GFP infiltrated *N. benthamiana* leaves and analyzed by chiral phase GC-MS. **c** (R)-linalool emission levels in transgenic *N. tabacum* lines overexpressing *McTPS1* (\#1 and \#2, McTPS1-OX) or wild-type (WT). Volatiles were collected by headspace method and analyzed by chiral phase GC-MS method. The quantity of (R)-linalool was determined by calculating the peak area of internal standard camphor (100 ng/μl). Error bars indicate SD (*n* = 2). nd, not detectable

Additionally, we generated more than 10 lines of transgenic *N. tabacum* plants expressing *McTPS1* (Additional file [1](#MOESM1){ref-type="media"}: Figure S8). One-month-old T1 progenies from two independent transgenic lines expressing *McTPS1* were used for headspace collection of VOCs using a push-pull headspace collection system. Figure [7c](#Fig7){ref-type="fig"} show that (R)-linalool emission in the transgenic lines were positively correlated with expression levels of *McTPS1* while no linalool was detected in wild type (WT), which was the same result as that from transient expression in the leaves of *N. benthamiana*.

Discussion {#Sec10}
==========

Our analysis on VOCs of champak flowers showed that over 85% of VOCs is constituted by terpenoids and volatile esters. Among the terpenoids, β-linalool and its derivatives (17.5%), β-ocimene (6.8%) and β-cubebene (8.6%) were the major compounds. However, we could not detect any volatile diterpenoids emitted from champak flowers suggesting that diterpenoids may not contribute to floral scent of champak. Among five monoterpenoids identified, β-linalool, β-ocimene and their derivatives were dominant accounting for more than 94% of total monoterpenoids. Both monoterpenes are not only known to be major components of floral fragrance \[[@CR4], [@CR40]\] but are also very common plant VOCs which are released after herbivore damage \[[@CR41]\]. Furanoid and pyranoid linalool oxides are found as prominent constituents of the floral scent along with β-linalool in *Clarkia breweri* \[[@CR42]\] and *Magnolia kobus* \[[@CR43]\]. Further studies in *Arabidopsis* CYP76C1 mutants, which is a major linalool metabolizing oxygenase demonstrated linalool oxides as a repellent for a number of insect taxa \[[@CR44]\]. Linalool metabolism might also serve as a detoxification mechanism of linalool, which is shown to be cytotoxic \[[@CR45]\].

2-Methylbutanoate and its derivatives such as methyl 2-methylbutanoate, ethyl 2-methyl butanoate and propyl 2-methylbutanoate were some of the major components of champak floral VOCs. These are known to possess a fruity odor and contribute to flavors of strawberry and pineapple \[[@CR46], [@CR47]\]. Methyl 2-methylbutanoate was identified as the key volatile of *Magnolia ovata* in attracting cyclocephaline scarab beetles \[[@CR48]\], suggesting that champak might be a cyclocephaline-pollinated flower. It is known that cyclocephaline and moth-pollinated flowers produce the attractant volatile in much larger quantities than the amount of attractant chemical produced by bee-pollinated flowers \[[@CR49]\]. Our VOCs analysis suggests that champak flowers produce large quantities of 2-methylbutanoate derivatives to attract pollinators, at the same time, they convert linalool to linalool oxides to protect their reproductive tissues from insects.

Our metabolite-guided RNA-seq approach provided a framework to identify pathways and genes for the biosynthesis of floral VOCs from champak flowers. GPP is the universal precursor for all monoterpenes whose synthesis is catalyzed by gernayl diphosphate synthase. This enzyme can exist as homodimers or heterodimers depending on the plant species \[[@CR50]\]. Our RNA-seq data indicates the GPS of champak is a heterodimer. GPS.SSU was one of the most abundantly expressed genes in champak flowers. The non-catalytic GPS.SSU is necessary to interact with the catalytic large subunit of GPS (GPS.LSU) to form the functional enzyme \[[@CR50]\]. The unigene encoding the GPS.LSU was found to be expressed at lower levels when compared to the small subunit. Two transcripts encoding SAMT were among the top 20 highly expressed champak unigenes. SAMTs catalyze the transfer of methyl group from S-adenosyl-L-methionine (SAM) to various substrates and are essential for the formation of many volatile methyl esters such as methyl benzoate, methyl salicylate, and methyl jasmonate \[[@CR51]\]. Although we could not identify methyl benzoate, methyl salicylate, or methyl jasmonate from VOCs of champak flowers, we were able to detect other volatile methyl esters such as 2-methylbutanoate, methyl 2-methylbutanoate, methyl tiglate, methyl hexanoate, methyl anthranilate, methyl palmitate, methyl linoleate and methyl linolenate (Fig. [1](#Fig1){ref-type="fig"} and Table [1](#Tab1){ref-type="table"}), all or some of which might be produced by the action of SAMTs. Not surprisingly, the genes encoding methionine synthase and SAM synthase that are crucial for biosynthesis of the methyl donor SAM were also among the most expressed unigenes. *Alcohol acyltransferase* (*AAT*) gene was also among the abundant transcripts, which catalyze the final step in the production of volatile esters by accepting a large range of alcohols and acyl-CoAs as substrates \[[@CR52]\]. Therefore, AAT could be involved in the production of ethyl 2-methylbutanoate and propyl 2-methylbutanoate in champak flowers, respectively (Additional file [1](#MOESM1){ref-type="media"}: Figure S2).

Among the four DXS subfamily unigenes analyzed, McDXS1 and McDXS2 but not McDXS3 and McDXS4 encoded functional DXS enzymes when tested in a *dxs*-deficient *E. coli* strain. Members of DXS clade II are shown to play a significant role in secondary metabolites synthesis \[[@CR53]\]. The two active DXS in champak, McDXS1 and McDXS2, belong to clade II suggesting that these 2 DXSs contribute towards the production of monoterpenes in champak flowers. McDXS3 and McDXS4 might catalyze a reaction distinct from that of McDXS1 and McDXS2. Since the amino acids essential for binding of G3P and TPP were conserved in McDXS4, it could catalyze a distinct reaction involving DXS substrates similar to CmDXS3/DXL from melon \[[@CR54]\]. Unlike most plants that possess two isoforms of *HDR* \[[@CR54]\], only single *HDR* (*McHDR*) was found from our champak RNA-seq (Additional file [1](#MOESM1){ref-type="media"}: Table S1). *McHDR* transcript level was more than 4 times higher in champak flowers than in leaves (Fig. [5a](#Fig5){ref-type="fig"}). This observation could either mean that champak has a single copy of *HDR* or the two isoforms may display spatial and temporal differential expression.

Linalool is one of the most common components of floral scent in a number of plant species \[[@CR2]\]. Linalool synthases have been characterized from many plants \[[@CR55]--[@CR57]\]. Unlike many TPSs that form multiple products from a single substrate \[[@CR4], [@CR58]\], linalool synthases are likely to catalyze the formation of linalool exclusively \[[@CR59]\]. The exception is the linalool synthase (PaLinS) from gymnosperm *Picea abies* (Norway spruce), which also produces very low quantities of nine monoterpenes including trans-β-ocimene, myrcene, α-terpinolene and 3-carene in vitro \[[@CR58]\]. However, these minor compounds were not detected in transgenic tobacco plants expressing PaLinS \[[@CR60]\].

Linalool synthases from multiple species do not share high amino acid sequence similarity. Owing to their substantially different sequences, linalool synthases from angiosperms were reported to belong to TPS-b, TPS-g and TPS-f subgroups of the TPS phylogenetic tree \[[@CR61]\]. Even though McTPS1 belongs to TPS-b subgroup with the R(R)X~8~W motif for monoterpene cyclization and shares highest sequence similarity to trans-β-ocimene synthase from *Litsea cubeba*, it only produced an acyclic monoterpene alcohol, (R)-linalool in in vitro, in vivo assays and in transgenic plants. Linalool occurs as two enantiomers in nature: (R)-linalool in lavender and bay laurel \[[@CR62]\]*,* whereas (S)-linalool is found in coriander \[[@CR63]\] and *Clarkia breweri* \[[@CR42]\]. (R)- and (S)-forms differ considerably in their olfactory qualities: (R)-form has a woody and lavender-like scent whereas the (S)-form as petitgrain-like and floral \[[@CR64]\].

In addition to (R)-linalool, we found trans-furanoid and trans-pyranoid linalool oxides to contribute majorly to champak floral scent. Trans-furanoid linalool oxide was identified as the major product of biotransformation of (R)-linalool in *Aspergillus niger* \[[@CR65]\]*.* Moreover, in fungus, biotransformation of linalool to furanoid and pyranoid linalool oxides was postulated to have epoxylinalool as an intermediate \[[@CR66]\]. These results indicate that (R)-linalool could be converted into trans-furanoid and trans-pyranoid linalool oxides in champak flowers. The oxidation of (S)-linalool to epoxylinalool involved a single cytochrome P450 enzyme (CYP71B31) in *Arabidopsis* \[[@CR67]\]. Additionally, cytochrome P450 76 family genes were also shown to be involved in linalool metabolism in *Arabidopsis* \[[@CR44]\]. From our champak transcriptome data, we found few homologues of *Arabidopsis* cytochrome P450s involved in linalool metabolism. These would be the ideal candidates to analyze and identify the genes involved in the formation of furanoid and pyranoid linalool oxides.

Conclusions {#Sec11}
===========

Availability of extensive genome resources for the commercially significant non-model plant *Magnolia champaca* is lacking. In this study, we performed de novo transcriptome assembly of high quality reads generated through Illumina paired end sequencing. Transcripts for the enzymes involved in biosynthesis of terpenes, phenylpropanoids/benzenoids and volatile esters, predominant VOCs in champak flowers were identified. Expression levels analysis of these unigenes showed that most of them were expressed higher in flowers than in leaves. Moreover, functions of the enzymes involved in the first and the last step of MEP pathway were validated by *E. coli* complementation assay. Finally, we also characterized the function of one of the most highly expressed terpene synthase as (R)-linalool synthase using in vitro, in vivo and transgenic studies. Our work will facilitate new gene discovery and provide crucial information for future genetic studies in champak. It will also serve as platform for metabolic engineering of champak floral-scent related genes.

Methods {#Sec12}
=======

Plant materials {#Sec13}
---------------

Fully opened champak flowers were collected from the trees (10--12 m tall and 25--30 cm in diameter) grown in the National University of Singapore campus during the month of November for VOCs analysis and RNA-seq. Young leaves of about 10--12 cm length were collected for qRT-PCR analysis. All sample collection was carried out at 7:00--8:30 am. Fresh samples unattacked by herbivores/aphids were chosen for the study. The samples were used immediately for VOCs analysis using a push-pull headspace collection system and total RNA extraction. *Nicotiana benthamiana* plants were grown in greenhouse under long day condition (16 h L/8 h D) for 4 weeks before using them for in vivo assays and subcellular localization experiments. *Nicotiana tabacum* plants were used for *Agrobacterium*-mediated transformation \[[@CR68]\].

VOCs collection and analysis {#Sec14}
----------------------------

For analysis of headspace chemical composition of champak flowers, floral VOCs were collected from 100 g of freshly collected flowers for 6 h in a tissue culture room at 25 °C using a push-pull headspace collection system \[[@CR69]\]. Champak flowers were placed in a glass jar (15 cm diameter, 20 cm height; UFO Labglass, Singapore). A compressed air pump was used to pull headspace air through a sorbent trap filled with HayeSep-Q trap (80/100 mesh size; Restek, USA). One μl (100 μg/ml) of camphor was added as an internal standard to each sorbent trap. The floral VOCs along with camphor were extracted twice from sorbent traps with 200 μl of hexane and analyzed by Agilent GC 7890A with 5975C inert mass selective detector, equipped with a HP-5MS column (30 m × 0.25 mm, 0.25 μm film thickness; Agilent Technologies, USA). 5 μl sample was injected into the column heated to 250 °C and the temperature was increased from 50 °C (1 min hold) to 300 °C (1 min hold) at the rate of 8 °C min^−1^. Retention indices (RI) were calculated by using C~7~-C~30~ n-alkanes standard. The compounds were identified by comparison with mass spectra reference library NIST MS 2014 and by using RI match. The data were processed by MSD ChemStation Data Analysis (Agilent Technologies). The internal standard camphor was used in calculating the quantity of other compounds.

In order to determine the linalool enantiomer emitted by champak, gas chromatography (GC) was carried out with a CP-Chirasil Dex CB column (25 m × 0.25 mm, 0.25 μm film thickness; Agilent Technologies, USA) and mass spectrometry (MS) detector at 220 °C. The headspace extracts from champak flowers, in vitro and in vivo assay products were analyzed using split injection into the column heated to 250 °C and the temperature was increased from 40 °C (1 min hold) to 200 °C at the rate of 4 °C min^−1^. The enantiomeric identity of linalool was confirmed by comparison of GC data to authentic (±)-linalool and (R)-linalool standards.

RNA isolation for RNA sequencing (RNA-seq) {#Sec15}
------------------------------------------

Total RNA was isolated from homogenized flower sample using the Spectrum™ Plant Total RNA Kit (Sigma-Aldrich). The quantity and quality of RNA were measured by a Nanodrop spectrophotometer (ND-1000, Thermo Fisher Scientific) and Agilent 2100 Bioanalyzer and RNA 6000 Nano Labchip Kit (Agilent Technologies), respectively. RNA sample with RNA Integrity Number (RIN) of \>7 was sent to the Rockefeller University Genomics Resource Center (New York, USA) to carry out next generation sequencing using Illumina HiSeq 2000. De novo assembly of the transcripts were performed as described in \[[@CR70]\].

cDNA synthesis and quantitative real time PCR (qRT-PCR) {#Sec16}
-------------------------------------------------------

1 μg of total RNA was used for cDNA synthesis with M-MLV reverse transcriptase (Promega), dNTP and oligo dT. The resulting 25 μl of RT reaction product was diluted to a total volume of 100 μl with RNase-free water and stored at −20 °C until use.

The expression levels of selected genes were analyzed by qRT-PCR. The cDNA sequences obtained from the RNA-seq data were exploited for designing primers using the Primer3 program \[[@CR71]\]. The primers used in this study are listed in Additional file [1](#MOESM1){ref-type="media"}: Table S4. The PCR reaction mixtures were subjected to the following conditions: (50 °C for 2 min, 95 °C for 10 min; 45 cycles of 95 °C for 15 s and 60 °C for 1 min) in the Applied Biosystems 7900HT fast real-time PCR system. *Actin* amplification was used as an internal normalization. Non-template control and non-RTase treated templates were included for each gene in order to eliminate the possibility of primer dimer formation and random genomic DNA contamination. All qRT-PCR experiments were carried out in triplicates with biological replicates. The results obtained were analyzed using SDS 2.4 software (Applied Biosystems).

Phylogenetic analysis {#Sec17}
---------------------

The deduced amino acid sequence of McDXSs and McTPS1 were aligned with DXSs and TPSs from other plant species using Clustal W under the following parameters: gap open - 3; gap extension - 1.8; Gonnet; penalties - on; gap separation - 4; cut off - 30%. Best-fit maximum likelihood models were identified for the alignment files using MEGA 6.0 software. The maximum likelihood trees for McDXS and McTPS1 were constructed with JTT + G5 and LG + G5 models respectively. 1000 bootstrap replications were made for each of the above trees.

Vector construction for *Agrobacterium*-mediated gene expression {#Sec18}
----------------------------------------------------------------

Full-length ORFs of *McDXS1--4*, *McHDR* and *McTPS1* were amplified from champak flower cDNAs by PCR. Gene specific primers were based on cDNA sequences of *McDXS1--4*, *McHDR* and *McTPS1* from champak RNA-seq data. The purified PCR products were cloned into pDONR221 vector using BP clonase (Invitrogen). All constructs were verified by DNA sequencing.

In order to generate YFP fusion constructs, pDONR221 containing *McDXS1*, *McDXS2*, *McDXS3*, *McDXS4*, *McHDR* and *McTPS1* were integrated into pBA-DC-YFP expression vector, which contains the cauliflower mosaic virus (CaMV) 35S promoter using LR clonase (Invitrogen).

Subcellular localization of McDXSs, McHDR and McTPS1 {#Sec19}
----------------------------------------------------

For subcellular localization, plasmids harbouring *McDXS1-YFP*, *McDXS2-YFP, McDXS3-YFP, McDXS4-YFP, McHDR-YFP* and *McTPS1-YFP* constructs were transformed into *Agrobacterium tumefaciens* GV3101 strain. The transformed GV3101 cells were then spread on LB plates with spectinomycin (100 μg/ml) and gentamycin (20 μg/ml) and incubated at 28 °C for two days. Cultures obtained from single colony of the above plates were infiltrated into the leaves of 4 weeks old *N. benthamiana* plants using a 1 ml syringe. The infiltrated leaves were excised three days post infiltration (dpi) and observed under a confocal laser-scanning microscope (Carl Zeiss LSM 5 Exciter) with a standard filter set. Images were analyzed by the Carl Zeiss's LSM image browser.

Functional complementation assays using *Escherichia coli* mutants {#Sec20}
------------------------------------------------------------------

Full-length ORFs of *McDXS1*, *McDXS2*, *McDXS3*, *McDXS4, MsHDR* along with *Arabidopsis DXS1* and *HDR* were inserted into gateway destination vector pDEST17 (Thermo Fisher Scientific) for expression in *E. coli*. Constructs containing *DXS* and *HDR* genes were transformed into *E. coli dxs* ^−^ strain, defective in DXS activity \[[@CR36]\] and *E. coli* HDR mutant strain MD1655 *ara*\<\>*ispH* \[[@CR37]\], respectively.

For DXS complementation assay, the transformed cells were grown overnight at 37 °C on LB agar plates containing 1 mM mevalonate and ampicillin (100 μg.ml^−1^). Colonies were then transferred to LB agar-ampicillin plates lacking mevalonate and incubated at 37 °C overnight. *E. coli dxs* ^−^ strain transformed with *AtDXS1* and empty pDEST17 vector served as positive and negative controls, respectively.

HDR complementation assay was done with same procedure as above except media composition. The transformants were selected on LB agar plates containing 0.2% (*w*/*v*) arabinose, ampicillin (100 μg.ml^−1^) and kanamycin (50 μg.ml^−1^), then shifted to LB agar-ampicillin/kanamycin plates containing 0.2% (*w*/*v*) glucose instead of arabinose to suppress the expression of endogenous *HDR* and incubated at 37 °C overnight. *AtHDR* and *McDXS1* were used as positive and negative controls, respectively.

Preparation of recombinant proteins {#Sec21}
-----------------------------------

To construct vector for the recombinant N-terminal 6His-tagged protein, pDONR221 clone possessing the ORF of *McTPS1* was integrated into the destination vector, pDEST17 to generate pDEST-McTPS1. The final construct was transformed into *E. coli* BL21 pLysS strain. *E. coli* extract after isopropyl β-D-1-thiogalactopyranoside (IPTG) induction was incubated with Ni-NTA Sepharose resin (Qiagen). The bound proteins were then eluted using 250 mM imidazole.

In vitro and in vivo TPS assay {#Sec22}
------------------------------

In vitro TPS assay was carried out by mixing 250 μl of 2X reaction buffer (50 mM HEPES pH 7.4, 200 mM KCl, 15 mM MgCl~2~, 10% glycerol, 10 mM DTT) with 20 μg purified recombinant protein and 10 μg of substrate (GPP and FPP) in an inert glass bottle. The reaction was mixed well and overlaid slowly with 250 μl of hexane. The reaction bottle was tightly closed and sealed with parafilm before incubating at 30 °C for 2 h. After 2 h incubation, the reaction mixture was vortexed for 1 min and centrifuged at 1200 rpm for 30 min. The hexane layer was then transferred to a fresh glass GC bottle (Agilent Technologies) and subjected to GC-MS (see [VOCs collection and analysis](#Sec14){ref-type="sec"}).

In vivo characterization of TPS using an *Agrobacterium*-mediated transient assay in *N. benthamiana* was performed according to Jin et al., (2015) with exception to the VOCs collection procedure. VOCs were collected and analyzed from four *N. benthamiana* plants infiltrated either with McTPS1-GFP or GFP using headspace method as described in [VOCs collection and analysis](#Sec14){ref-type="sec"}. The compound obtained was identified with the authentic (±)-linalool standard (Sigma-Aldrich) and by mass spectra reference library.

Additional file 1 {#Sec23}
=================

Additional file 1: Figure S1.Quality and KEGG analysis of champak RNA-seq. **Figure S2.** Proposed biosynthesis pathways of volatile ester 2-methylbutanoate and its derivatives via the catabolism of branched-chain amino acid L-isoleucine. **Figure S3.** Comparison of deduced amino acid sequence of representative genes from pathways responsible for the production of VOCs. **Figure S4.** Phylogenetic analysis and amino acid alignment of champak DXSs. **Figure S5.** Amino acid sequence alignment of McHDR. **Figure S6.** Amino acid sequence alignment of McTPS1. **Figure S7.** GC-MS chiral analysis of β-linalool emitted from champak flowers. **Figure S8.** Analysis of transgenic *N. tabacum* overexpressing *McTPS1.* **Table S1.** MEP pathway genes from champak RNA-seq. **Table S2.** TPS genes from champak RNA-seq.**Table S3.** Accession numbers of proteins used in the TPS phylogenetic analysis. **Table S4.** List of primers used in this study. **Table S5.** Accession numbers of proteins used in the amino acid sequence alignments. **Table S6.** Accession numbers of proteins used in the DXS phylogenetic analysis. (PDF 1708 kb) Additional file 2:Functional annotation of champak unigenes using BLASTX against the multiple protein databases. (XLS 13237 kb)

CaMV

:   Cauliflower mosaic virus

CYP450

:   Cytochrome P450s

FAME

:   Fatty acid-derived methyl esters

FPKM

:   Fragments per kilobase of transcript per million mapped reads

GC-MS

:   Gas chromatography-mass spectrometry

GFP

:   Green fluorescent protein

GO

:   Gene ontology

IPTG

:   Isopropyl β-D-1-thiogalactopyranoside

KAAS

:   KEGG automatic annotation server

KEGG

:   Kyoto encyclopedia of genes and genomes

KO

:   KEGG orthology

MEP

:   Methylerythritol 4-phosphate

MVA

:   Mevalonate

NCBI

:   National centre for biotechnology information

NR

:   Non-redundant

ORF

:   Open reading frame

qRT-PCR

:   Quantitative real-time polymerase chain reaction

RNA-seq

:   RNA sequencing

TP

:   Transit peptide

TPS

:   Terpene synthase

VOC

:   Volatile organic compound

WT

:   Wild type

YFP

:   Yellow fluorescent protein

**Electronic supplementary material**

The online version of this article (doi:10.1186/s12864-017-3846-8) contains supplementary material, which is available to authorized users.

We thank Dr. Micheal A. Phillips (Centre for Research in Agricultural Genomics, Barcelona, Spain) for providing us the *E.coli* strains, *dxs-* and *MG1655 ara\<\>ispH*. We are also grateful to Deepa Panicker for champak RNA isolation, Penny Yeo Liu Qing for the technical support on chiral column chromatography, and the Temasek Life Sciences Laboratory central facility for confocal microscopy.

Funding {#FPar1}
=======

This work was supported by the Singapore National Research Foundation (Competitive Research Programme Award No: NRF-CRP8--2011-02) and Economic Development Board (COY-15-RISC/120014 S12--1146-RISC). These funding sources had no role in study design, data collection and interpretation and in writing the manuscript.

Availability of data and materials {#FPar2}
==================================

The raw RNA-seq data supporting the result of this work is available in the DDBJ: DNA Bank of Japan, under the accession number: DRA005500 (<http://trace.ddbj.nig.ac.jp/DRASearch/submission?acc=DRA005500>). Nucleotide sequences from this study are available in the NCBI (<https://www.ncbi.nlm.nih.gov>) with accession numbers McDXS1 (KY502270), McDXS2 (KY502271), McDXS3 (KY502272), McDXS4 (KY502273), McHDR1 (KY502274) and McTPS1 (KY502275).

Authors' contributions {#FPar3}
======================

RS, NHC and ICJ conceived the project; JJ performed the de novo assembly of RNA-seq; ICJ supervised the experiments; VS and SD analyzed the transcriptome data; SD, RS, NHC and ICJ designed the experiments and interpreted data; SD performed the molecular and biochemical experiments; SD and ICJ wrote the manuscript. All authors read and approved the final manuscript.

Competing interests {#FPar4}
===================

The authors declare that they have no competing interests.

Consent for publication {#FPar5}
=======================

Not applicable.

Ethics approval and consent to participate {#FPar6}
==========================================

*Magnolia champaca,* previously classified as *Michelia champaca,* is widely cultivated in Singapore. *Magnolia champaca* is not listed in the appendices I, II and III of the Convention on the Trade in Endangered Species of Wild Fauna and Flora, that has been valid from 4 April 2017 (<https://cites.org/eng/app/appendices.php>). Champak samples used in this study were collected from the trees grown at Temasek Life Sciences Laboratory in the National University of Singapore. Collection of plant materials complied with the institutional, national and international guidelines. No specific permits were required. Vouchers of this species have been placed in the Singapore herbarium.

Publisher's Note {#FPar7}
================

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
