Molecular regulation of santalol biosynthesis in Santalum album L

Identification of Transcription Factors of Santalene Synthase Gene Promoters and SaSSY Cis-Elements through Yeast One-Hybrid Screening in Santalum album L

The main components of sandalwood heartwood essential oil are terpenoids, approximately 80% of which are α-santalol and β-santalol. In the synthesis of the main secondary metabolites of sandalwood heartwood, the key gene, santalene synthase (SaSSY), can produce α-santalene and β-santalene by catalyzed (E, E)-FPP. Furthermore, santalene is catalyzed by the cytochrome monooxygenase SaCYP736A167 to form sandalwood essential oil, which then produces a fragrance. However, the upstream regulatory mechanism of the key gene santalene synthase remains unclear. In this study, SaSSY (Sal3G10690) promoter transcription factors and SaSSY cis-elements were screened. The results showed that the titer of the sandalwood cDNA library was 1.75 × 107 CFU/mL, 80% of the inserted fragments identified by PCR were over 750 bp in length, and the positivity rate of the library was greater than 90%. The promoter region of the SaSSY gene was shown to have the structural basis for potential regulatory factor binding. After sequencing and bioinformatics analysis, we successfully obtained 51 positive clones and identified four potential SaSSY transcriptional regulators. Sal6G03620 was annotated as the transcription factor MYB36-like, and Sal8G07920 was annotated as the small heat shock protein HSP20 in sandalwood. Sal1G00910 was annotated as a hypothetical protein of sandalwood. Sal4G10880 was annotated as a homeobox-leucine zipper protein (ATHB-15) in sandalwood. In this study, a cDNA library of sandalwood was successfully constructed using a yeast one-hybrid technique, and the transcription factors that might interact with SaSSY gene promoters were screened. This study provides a foundation for exploring the molecular regulatory mechanism involved in the formation of sandalwood heartwood.

Integrated mRNA and small RNA sequencing reveals post-transcriptional regulation of the sesquiterpene pathway in Santalum album L. (Indian sandalwood)

Key message

In sandalwood, negative pattern of regulation by miRNAs was documented in key genes from the sesquiterpene pathway, with cytochrome P450 reductase showing maximum miRNA targets, followed by sesquisabianene synthase 1.

Abstract

A comprehensive knowledge of the molecular regulation of sesquiterpene biosynthetic pathway through transcriptomic studies is well established in Santalum album (Indian Sandalwood). However, the post-transcriptional regulation of the genes regulating the pathway is still elusive in this genus. In the present study, an integrated analysis of wood transcriptome and small RNA datasets was conducted to investigate the role of miRNAs in regulating the expression of transcripts involved in santalol production mediated by the sesquiterpene biosynthesis pathway. A total of 24,237 transcripts were annotated from the wood transcriptome, and 45 transcripts were mapped to the sesquiterpenoid pathway. Small RNA data analysis identified 257 conserved miRNAs belonging to 50 families and 7 novel putative miRNAs. Sa-miR156, Sa-miR396, Sa-miR166, and Sa-miR319 had the most number of members among the miRNA families. An integrated analysis predicted 69 miRNA members belonging to 12 families that targeted 12 transcripts from the sesquiterpene pathway, with a maximum of 24 miRNAs regulating cytochrome P450 reductase, followed by sesquisabianene synthase 1, which was targeted by 23 miRNAs. Validation of miRNA-mRNA interaction by qRT-PCR revealed a negative pattern of regulation in six miRNA-mRNA target pairs across wood tissues sourced from four genotypes. The present study provides the first crucial insight into the post-transcriptional regulation of the sesquiterpene pathway genes in the genus Santalum and opens up a new perspective in metabolite engineering for enhanced essential oil production in sandalwood.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03816-4.

Sesquiterpene Synthase Engineering and Targeted Engineering of α-Santalene Overproduction in Escherichia coli

As a natural sesquiterpene compound with numerous biological activities, α-santalene has extensive applications in the cosmetic and pharmaceutical industries. Although several α-santalene-producing microbial strains have been constructed, low productivity still hampers large-scale fermentation. Herein, we present a case of engineered sesquiterpene biosynthesis where the insufficient downstream pathway capacity limited high-level α-santalene production in Escherichia coli. The initial strain was constructed, and it produced 6.4 mg/L α-santalene. To increase α-santalene biosynthesis, we amplified the flux toward farnesyl diphosphate (FPP) precursor by screening and choosing the right FPP synthase and reprogrammed the rate-limiting downstream pathway by generating mutations in santalene synthase (Clausena lansium; ClSS). Santalene synthase was engineered by site-directed mutagenesis, resulting in the improved soluble expression of ClSS and an α-santalene titer of 887.5 mg/L; the α-santalene titer reached 1078.8 mg/L after adding a fusion tag to ClSS. The most productive pathway, which included combining precursor flux amplification and mutant synthases, conferred an approximate 169-fold increase in α-santalene levels. Maximum titers of 1272 and 2916 mg/L were achieved under shake flask and fed-batch fermentation, respectively, and were among the highest levels reported using E. coli as the host.

Comparative transcriptome profiling of high and low oil yielding Santalum album L

East Indian Sandalwood (Santalum album L.) is highly valued for its heartwood and its oil. There have been no efforts to comparative study of high and low oil yielding genetically identical sandalwood trees grown in similar climatic condition. Thus we intend to study a genome wide transcriptome analysis to identify the corresponding genes involved in high oil biosynthesis in S. album. In this study, 15 years old S. album (SaSHc and SaSLc) genotypes were targeted for analysis to understand the contribution of genetic background on high oil biosynthesis in S. album. A total of 28,959187 and 25,598869 raw PE reads were generated by the Illumina sequencing. 2.12 million and 1.811 million coding sequences were obtained in respective accessions. Based on the GO terms, functional classification of the CDS 21262, & 18113 were assigned into 26 functional groups of three GO categories; (4,168; 3,641) for biological process (5,758;4,971) cellular component and (5,108;4,441) for molecular functions. Total 41,900 and 36,571 genes were functionally annotated and KEGG pathways of the DEGs resulted 213 metabolic pathways. In this, 14 pathways were involved in secondary metabolites biosynthesis pathway in S. album. Among 237 cytochrome families, nine groups of cytochromes were participated in high oil biosynthesis. 16,665 differentially expressed genes were commonly detected in both the accessions (SaHc and SaSLc). The results showed that 784 genes were upregulated and 339 genes were downregulated in SaHc whilst 635 upregulated 299 downregulated in SaSLc S. album. RNA-Seq results were further validated by quantitative RT-PCR. Maximum Blast hits were found to be against Vitis vinifera. From this study, we have identified additional number of cytochrome family in high oil yielding sandalwood accessions (SaHc). The accessibility of a RNA-Seq for high oil yielding sandalwood accessions will have broader associations for the conservation and selection of superior elite samples/populations for further genetic improvement program.

Optimized biosynthesis of santalenes and santalols in Saccharomyces cerevisiae

Santalenes and santalols from Santalum album are the main components of the valuable spice sandalwood essential oil, which also has excellent pharmacological activities such as antibacterial, anti-inflammatory, and antitumor. Firstly, we constructed biosynthesis pathways of santalenes by synthetic biology strategy. The assembled biosynthetic cassettes were integrated into the multiple copy loci of δ gene in S. cerevisiae BY4742 with assistance of pDi-CRISPR, and 94.6 mg/L santalenes was obtained by shake flask fermentation of engineered yeast. Secondly, a selected optimized P450-CPR redox system was integrated into the chromosome of the santalenes-producing strain with a single copy, and 24.6 mg/L santalols were obtained. Finally, the yields of santalenes and santalols were increased to 164.7 and 68.8 mg/L, respectively, by downregulating ERG9 gene. This is the first report on the de novo synthesis of santalols by P450-CPR chimera in S. cerevisiae. Meanwhile, the optimized chimeric CYP736A167^opt-46tATR1^opt exhibits higher activity to oxidize santalenes into santalols. It would provide a feasible solution for the optimal biosynthesis of santalols. KEY POINTS: • First-time de novo synthesis of santalols by P450-CPR chimera in S. cerevisiae. • Truncated 46tATR1 has higher activity than that of CPR2. • Yields of santalenes and santalols were increased by downregulating ERG9 gene.

Cloning and expression analysis of mevalonate kinase and phosphomevalonate kinase genes associated with the MVA pathway in Santalum album

Sandalwood (Santalum album L.) is highly valued for its fragrant heartwood and extracted oil. Santalols, which are the main components of that oil, are terpenoids, and these are biosynthesized via the mevalonic acid (MVA) pathway. Mevalonate kinase (MK) and phosphomevalonate kinase (PMK) are key enzymes in the MVA pathway. Little is known about the genes that encode MK and PMK in S. album or the mechanism that regulates their expression. To isolate and identify the functional genes involved in santalol biosynthesis in S. album, an MK gene designated as SaMK, and a PMK gene designated as SaPMK, were cloned from S. album. The sequences of these genes were analyzed. A bioinformatics analysis was conducted to assess the homology of SaMK and SaPMK with MK and PMK genes from other plants. The subcellular localization of SaMK and SaPMK proteins was also investigated, as was the functional complementation of SaMK and SaPMK in yeast. Our results show that the full-length cDNA sequences of SaMK and SaPMK were 1409 bp and 1679 bp long, respectively. SaMK contained a 1381 bp open reading frame (ORF) encoding a polypeptide of 460 amino acids and SaPMK contained a 1527 bp ORF encoding a polypeptide of 508 amino acids. SaMK and SaPMK showed high homology with MK and PMK genes of other plant species. Functional complementation of SaMK in a MK-deficient mutant yeast strain YMR208W and SaPMK in a PMK-deficient mutant yeast strain YMR220W confirmed that cloned SaMK and SaPMK cDNA encode a functional MK and PMK, respectively, mediating MVA biosynthesis in yeast. An analysis of tissue expression patterns revealed that SaMK and SaPMK were constitutively expressed in all the tested tissues. SaMK was highly expressed in young leaves but weakly expressed in sapwood. SaPMK was highly expressed in roots and mature leaves, but weakly expressed in young leaves. Induction experiments with several elicitors showed that SaMK and SaPMK expression was upregulated by methyl jasmonate. These results will help to further study the role of MK and PMK genes during santalol biosynthesis in S. album.

Molecular Cloning and Functional Analysis of 1-Deoxy-D-Xylulose 5-Phosphate Reductoisomerase from Santalum album

Sandalwood (Santalum album L.) heartwood-derived essential oil contains a high content of sesquiterpenoids that are economically highly valued and widely used in the fragrance industry. Sesquiterpenoids are biosynthesized via the mevalonate acid and methylerythritol phosphate (MEP) pathways, which are also the sources of precursors for photosynthetic pigments. 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) is a secondary rate-limiting enzyme in the MEP pathway. In this paper, the 1416-bp open reading frame of SaDXR and its 897-bp promoter region, which contains putative conserved cis-elements involved in stress responsiveness (HSE and TC-rich repeats), hormone signaling (abscisic acid, gibberellin and salicylic acid) and light responsiveness, were cloned from 7-year-old S. album trees. A bioinformatics analysis suggested that SaDXR encodes a functional and conserved DXR protein. SaDXR was widely expressed in multiple tissues, including roots, twigs, stem sapwood, leaves, flowers, fruit and stem heartwood, displaying significantly higher levels in tissues with photosynthetic pigments, like twigs, leaves and flowers. SaDXR mRNA expression increased in etiolated seedlings exposed to light, and the content of chlorophylls and carotenoids was enhanced in all 35S::SaDXR transgenic Arabidopsis thaliana lines, consistent with the SaDXR expression level. SaDXR was also stimulated by MeJA and H₂O₂ in seedling roots. α-Santalol content decreased in response to fosmidomycin, a DXR inhibitor. These results suggest that SaDXR plays an important role in the biosynthesis of photosynthetic pigments, shifting the flux to sandalwood-specific sesquiterpenoids.

Cloning, characterization, and functional analysis of acetyl-CoA C-acetyltransferase and 3-hydroxy-3-methylglutaryl-CoA synthase genes in Santalum album

Sandalwood (Santalum album L.) is famous for its unique fragrance derived from the essential oil of heartwood, whose major components are santalols. To understand the mechanism underlying the biosynthesis of santalols, in this study, we cloned two related genes involved in the mevalonate pathway in S. album coding for acetyl-CoA C-acetyl transferase (AACT) and 3-hydroxy-3-methyglutary-CoA synthase (HMGS). These genes were characterized and functionally analyzed, and their expression profiles were also assessed. An AACT gene designated as SaAACT (GenBank accession No. MH018694) and a HMGS gene designated as SaHMGS (GenBank accession No. MH018695) were successfully cloned from S. album. The deduced SaAACT and SaHMGS proteins contain 415 and 470 amino acids, and the corresponding size of their open-reading frames is 1538 bp and 1807 bp, respectively. Phylogenetic trees showed that the SaAACT protein had the closest relationship with AACT from Hevea brasiliensis and the SaHMGS proteins had the highest homology with HMGS from Siraitia grosvenorii. Functional complementation of SaAACT and SaHMGS in a mutant yeast strain deficient in these proteins confirmed that SaAACT and SaHMGS cDNA encodes functional SaAACT and SaHMGS that mediate mevalonate biosynthesis in yeast. Tissue-specific expression analysis revealed that both genes were constitutively expressed in all examined tissues (roots, sapwood, heartwood, young leaves, mature leaves and shoots) of S. album, both genes showing highest expression in roots. After S. album seedlings were treated with 100 μM methyl jasmonate, the expression levels of SaAACT and SaHMGS genes increased, suggesting that these genes were responsive to this elicitor. These studies provide insight that would allow further analysis of the role of genes related to the sandalwood mevalonate pathway in the regulation of biosynthesis of sandalwood terpenoids and a deeper understanding of the molecular mechanism of santalol biosynthesis.

Advances in biotechnological production of santalenes and santalols

Sandalwood essential oil has been widely used not only as natural medicines but also in perfumery and food industries, with sesquiterpenoids as its major components including (Z)- α-santalol and (Z)-β-santalol and so on. The mature heartwoods of Santalum album, Santalum austrocaledonicum and Santalum spicatum are the major plant resources for extracting sandalwood essential oil, which have been overexploited. Synthetic biology approaches have been successfully applied to produce natural products on large scale. In this review, we summarize biosynthetic enzymes of santalenes and santalols, including various santalene synthases (STSs) and cytochrome P450 monooxygenases (CYPs), and then highlight the advances of biotechnological production of santalenes and santalols in heterologous hosts, especially metabolic engineering strategies for constructing santalene- and santalol-producing Saccharomyces cerevisiae.

Assessment of morphological and genetic variability through genic microsatellite markers for essential oil in Sandalwood (Santalum album L.)

Sandalwood (Santalum album L; family Santalaceae) is a highly significant aromatic oil yielding tree. It is valued for two important traits, heartwood and essential oil obtained from the heartwood. This study was proposed to assess the morphological and genetic variability of sandalwood accessions. For this, genotypes were randomly selected (n = 177) from the 14 populations from three states in southern India. The total heartwood oil content and quality was estimated by UV method and GC-MS. Total 14 oil-specific genic SSR markers were procured to evaluate the genetic diversity among the sandalwood accessions. Total core size, heartwood content, and oil of S. album ranged from 4.4 to 19.1 cm; 0.0 to 17.3 cm; and 0.0 to 5.96% with covariance 27.61, 85.25, and 73.12% followed by mean 9.74, 3.77, and 2.71, respectively. Genetic diversity estimates were highly polymorphic in terms of Na 7.28, Ne 5.89, He 8.0 PIC 0.891, with little Ho, and F-0.922. AMOVA revealed that minimal genetic variation among populations and highest variation was found among individuals with Nm (58.4). The UPGMA reveals the cluster favored the grouping pattern by the PCA analysis. Structure and PCA analysis clustered the entire populations into two major groups with F ST 0.046 in which population of Kerala and Karnataka were pure and Telangana accessions were found admixtures. No significant correlation (r 2 = 0.23, P = 0.00) was observed between heartwood oil and genetic structures. A high degree of transferability of genic markers would facilitate the assessment of novel genotypes for future tree improvement and conservation of Sandalwood populations.

Identification and functional characterization of three new terpene synthase genes involved in chemical defense and abiotic stresses in Santalum album

BACKGROUND

It is well known that aromatic essential oils extracted from the heartwood of Santalum album L. have wide economic value. However, little is known about the role of terpenoids in response to various adverse environmental stresses as other plants do in the form of signals during plant-environment interactions.

RESULTS

In this study, trace amounts of volatiles consisting of α-santalene, epi-β-santalene, β-santalene, α-santalol, β-santalol, (E)-α-bergamotene, (E)-β-farnesene and β-bisabolene were found in the leaves of mature S. album trees. We identified more than 40 candidate terpene synthase (TPS) unigenes by mining publicly-available RNA-seq data and characterized the enzymes encoded by three cDNAs: one mono-TPS catalyzes the formation of mostly α-terpineol, and two multifunctional sesqui-TPSs, one of which produces (E)-α-bergamotene and sesquisabinene as major products and another which catalyzes the formation of (E)-β-farnesene, (E)-nerolidol and (E,E)-farnesol as main products. Metabolite signatures and gene expression studies confirmed that santalol content is closely related with santalene synthase (SaSSY) transcripts in heartwood, which is key enzyme responsible for santalol biosynthesis. However, the expression of three new SaTPS genes differed significantly from SaSSY in the essential oil-producing heartwood. Increased activities of antioxidant enzymes, superoxide dismutase, catalase, peroxidase and ascorbate peroxidase, were detected in different tissues of S. album plants after applying 1 mM methyl jasmonate (MeJA) and 1 mM salicylic acid (SA), or exposure to 4°C, 38°C and high light intensity. MeJA and SA dramatically induced the expression of SaTPS1 and SaTPS2 in leaves. SaTPS1 to 3 transcripts were differentially activated among different tissues under adverse temperature and light stresses. In contrast, almost all SaSSY transcripts decreased in response to these environmental stresses, unlike SaTPS1 to 3.

CONCLUSIONS

Multifunctional enzymes were biochemically characterized, including one chloroplastic mono-TPS and two cytosolic sesqui-TPSs in sandalwood. Our results suggest the ecological importance of these three new SaTPS genes in defensive response to biotic attack and abiotic stresses in S. album.

Selection and Validation of Novel RT-qPCR Reference Genes under Hormonal Stimuli and in Different Tissues of Santalum album

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a widely used technique to investigate gene expression levels due to its high throughput, specificity, and sensitivity. An appropriate reference gene is essential for RT-qPCR analysis to obtain accurate and reliable results. To date, no reliable reference gene has been validated for the economically tropical tree, sandalwood (Santalum album L.). In this study, 13 candidate reference genes, including 12 novel putative reference genes selected from a large set of S. album transcriptome data, as well as the currently used β-actin gene (ACT), were validated in different tissues (stem, leaf, root and callus), as well as callus tissue under salicylic acid (SA), jasmonic acid methyl ester (MeJA), and gibberellin (GA) treatments using geNorm, NormFinder, BestKeeper, Delta Ct and comprehensive RefFinder algorithms. Several novel candidate reference genes were much more stable than the currently used traditional gene ACT. ODD paired with Fbp1 for SA treatment, CSA and Fbp3 for MeJA treatment, PP2C and Fbp2 for GA treatment, as well as Fbp1 combined with Fbp2 for the total of three hormone treatments were the most accurate reference genes, respectively. FAB1A, when combined with PP2C, was identified as the most suitable reference gene combination for the four tissues tested, while the combination of HLMt, PPR and FAB1A were the most optimal reference genes for all of the experimental samples. In addition, to verify our results, the relative expression level of the SaSSy gene was evaluated by the validated reference genes and their combinations in the three S. album tissues and under MeJA treatment. The evaluated reference genes in this study will improve the accuracy of RT-qPCR analysis and will benefit S. album functional genomics studies in different tissues and under hormone stimuli in the future.

Physiological and transcriptomic analyses reveal a response mechanism to cold stress in Santalum album L. leaves

Santalum album L. (Indian sandalwood) is an economically important plant species because of its ability to produce highly valued perfume oils. Little is known about the mechanisms by which S. album adapts to low temperatures. In this study, we obtained 100,445,724 raw reads by paired-end sequencing from S. album leaves. Physiological and transcriptomic changes in sandalwood seedlings exposed to 4 °C for 0-48 h were characterized. Cold stress induced the accumulation of malondialdehyde, proline and soluble carbohydrates, and increased the levels of antioxidants. A total of 4,424 differentially expressed genes were responsive to cold, including 3,075 cold-induced and 1,349 cold-repressed genes. When cold stress was prolonged, there was an increase in the expression of cold-responsive genes coding for transporters, responses to stimuli and stress, regulation of defense response, as well as genes related to signal transduction of all phytohormones. Candidate genes in the terpenoid biosynthetic pathway were identified, eight of which were significantly involved in the cold stress response. Gene expression analyses using qRT-PCR showed a peak in the accumulation of SaCBF2 to 4, 50-fold more than control leaves and roots following 12 h and 24 h of cold stress, respectively. The CBF-dependent pathway may play a crucial role in increasing cold tolerance.

A Chemometrics Approach to the Investigation of the Intraspecific Variability of the Volatile Oil of Eupatorium tremulum from Southern Brazil

Eupatorium tremulum is a South American shrub reported to cause cattle digestive intoxication, of which the volatile oil, mainly composed by bisabolane- and amorphane-type sesquiterpenoids, exhibits high quantitative variability. This report describes the application of chemometric tools for the identification of volatile compounds that characterize phenophasical changes in the plant. Preblooming, blooming, and postblooming specimens were paired-sampled and submitted to hydrodistillation and GC-MS analysis. Differential results were analyzed by orthogonal projection to latent structures-discriminant analysis (OPLS-DA), and the substances with different distribution in each phase were highlighted. Mean results between phases were submitted to factor analysis (FA), and correlations between the variables were demonstrated. Preblooming to blooming phase change was characterized by decreased levels of amorpha-4-en-7-ol (13) and epi-α-bisabolol (19) and increased amounts of amorpha-4,7(11)-diene (1). Blooming to postblooming change was characterized by decreases in 1, germacrene D (2), and β-bisabolene (4) and increases in 13 and 19. Finally, enhanced levels of 1, 4, and 2 reflected the change from the postblooming to the preblooming phase. FA revealed a strong correlation in the variability between the bisabolane hydrocarbons, possibly related to its common enzymatic origin. Another strong source of negative correlation showed bisabolane- and amorphane-type alcohols, on one side, and amorphane-type furans, on the other side, to occur in two alternative oxidation routes. Finally, 1 was strongly negatively correlated to its oxidized furan and ketofuran derivatives [verboccidentafuran (16) and 3-oxo-verboccidentafuran (23)] and additionally to a third compound, putatively identified as a biosynthetic intermediate between this hydrocarbon and the furans, amorpha-4,7(11)-dien-8-one (20).

Functional Characterization of Novel Sesquiterpene Synthases from Indian Sandalwood, Santalum album

Indian Sandalwood, Santalum album L. is highly valued for its fragrant heartwood oil and is dominated by a blend of sesquiterpenes. Sesquiterpenes are formed through cyclization of farnesyl diphosphate (FPP), catalyzed by metal dependent terpene cyclases. This report describes the cloning and functional characterization of five genes, which encode two sesquisabinene synthases (SaSQS1, SaSQS2), bisabolene synthase (SaBS), santalene synthase (SaSS) and farnesyl diphosphate synthase (SaFDS) using the transcriptome sequencing of S. album. Using Illumina next generation sequencing, 33.32 million high quality raw reads were generated, which were assembled into 84,094 unigenes with an average length of 494.17 bp. Based on the transcriptome sequencing, five sesquiterpene synthases SaFDS, SaSQS1, SaSQS2, SaBS and SaSS involved in the biosynthesis of FPP, sesquisabinene, β-bisabolene and santalenes, respectively, were cloned and functionally characterized. Novel sesquiterpene synthases (SaSQS1 and SaSQS2) were characterized as isoforms of sesquisabinene synthase with varying kinetic parameters and expression levels. Furthermore, the feasibility of microbial production of sesquisabinene from both the unigenes, SaSQS1 and SaSQS2 in non-optimized bacterial cell for the preparative scale production of sesquisabinene has been demonstrated. These results may pave the way for in vivo production of sandalwood sesquiterpenes in genetically tractable heterologous systems.