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

Improved chromosome-level genome assembly of Indian sandalwood (Santalum album)

Santalum album is a well-known aromatic and medicinal plant that is highly valued for the essential oil (EO) extracted from its heartwood. In this study, we present a high-quality chromosome-level genome assembly of S. album after integrating PacBio Sequel, Illumina HiSeq paired-end and high-throughput chromosome conformation capture sequencing technologies. The assembled genome size is 207.39 M with a contig N50 of 7.33 M and scaffold N50 size of 18.31 M. Compared with three previously published sandalwood genomes, the N50 length of the genome assembly was longer. In total, 94.26% of the assembly was assigned to 10 pseudo-chromosomes, and the anchor rate far exceeded that of a recently released value. BUSCO analysis yielded a completeness score of 94.91%. In addition, we predicted 23,283 protein-coding genes, 89.68% of which were functionally annotated. This high-quality genome will provide a foundation for sandalwood functional genomics studies, and also for elucidating the genetic basis of EO biosynthesis in S. album.

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.

Volatile secondary metabolome and transcriptome analysis reveals distinct regulation mechanism of aroma biosynthesis in Syringa oblata and S. vulgaris

Lilacs have high ornamental value due to their strong aroma. However, the molecular regulatory mechanisms of aroma biosynthesis and metabolism in lilac were largely unclear. In this study, two varieties with distinct aroma, Syringa oblata 'Zi Kui' (faint aroma) and Syringa vulgaris 'Li Fei' (strong aroma), were used for exploring the regulation mechanism of aroma difference. Via GC-MS analysis, a total of 43 volatile components were identified. Terpene volatiles was the most abundant volatiles constituting the aroma of two varieties. Notably, 3 volatile secondary metabolites were unique in 'Zi Kui' and 30 volatile secondary metabolites were unique in 'Li Fei'. Then, a transcriptome analysis was performed to clarify the regulation mechanism of aroma metabolism difference between these two varieties, and identified 6411 differentially expressed genes (DEGs). Interestingly, ubiquinone and other terpenoid-quinone biosynthesis genes were significantly enriched in DEGs. We further conducted a correlation analysis between the volatile metabolome and transcriptome and found that TPS, GGPPS, and HMGS genes might be the key contributors to the differences in floral fragrance composition between the two lilac varieties. Our study improves the understanding in the regulation mechanism of Lilac aroma and would help improve the aroma of ornamental crops by metabolic engineering.

Expression Analysis of CoHMGS in Cornus officinalis and Subcellular Localization of the Enzyme It Encodes

Cornus officinalis is a perennial deciduous tree or shrub. Its mature fruits are extracted and used in Traditional Chinese Medicine, called Shanzhuyu. The characteristic active components of C. officinalis include loganin and morroniside, which belong to iridoid glycosides. 3-Hydroxy-3-methylglutaryl-CoA synthase (HMGS) is a key enzyme in the cytoplasmic mevalonate pathway providing the precursor molecules isopentenyl pyrophosphate and dimethylallyl pyrophosphate for isoprenoid biosynthesis such as sterols, triterpenes, and their derivatives such as iridoid glycosides. Different concentrations of methyl jasmonate (MeJA) and ethephon (ETH) solutions were sprayed on C. officinalis seedlings, and the effect of hormones on CoHMGS gene expression was detected by real-time fluorescence quantitative PCR. The quantitative real-time PCR results showed that 750 mg/L ETH treatment had the most significant induction effect on CoHMGS gene expression. The HPLC analysis of extracts revealed that the treatment could also significantly increase the content of morroniside and loganin in the leaves of C. officinalis. By use of a CoHMGS-green fluorescent protein (GFP) fusion construct for heterologous expression in tobacco, laser scanning confocal microscopy revealed a cytoplasmic localization. This preliminary study of the CoHMGS gene could prepare the ground for more precisely elucidating the synthesis of secondary metabolite in C. officinalis.

Genome-wide transcriptome analysis reveals the regulatory network governing terpene trilactones biosynthesis in Ginkgo biloba

Ginkgo biloba L. is currently the only remaining gymnosperm of the Ginkgoaceae Ginkgo genus, and its history can be traced back to the Carboniferous 200 million years ago. Terpene trilactones (TTLs) are one of the main active ingredients in G. biloba, including ginkgolides and bilobalide. They have a good curative effect on cardiovascular and cerebrovascular diseases because of their special antagonistic effect on platelet-activating factors. Therefore, it is necessary to deeply mine genes related to TTLs and to analyze their transcriptional regulation mechanism, which will hold vitally important scientific and practical significance for quality improvement and regulation of G. biloba. In this study, we performed RNA-Seq on the root, stem, immature leaf, mature leaf, microstrobilus, ovulate strobilus, immature fruit and mature fruit of G. biloba. The TTL regulatory network of G. biloba in different organs was revealed by different transcriptomic analysis strategies. Weighted gene co-expression network analysis (WGCNA) revealed that the five modules were closely correlated with organs. The 12 transcription factors, 5 structural genes and 24 Cytochrome P450 (CYP450) were identified as candidate regulators for TTL accumulation by WGCNA and cytoscape visualization. Finally, 6 APETALA2/ethylene response factors, 2 CYP450s and bHLH were inferred to regulate the metabolism of TTLs by correlation analysis. This study is the comprehensive in authenticating transcription factors, structural genes and CYP450 involved in TTL biosynthesis, thereby providing molecular evidence for revealing the comprehensive regulatory network involved in TTL metabolism in G. biloba.

Cloning and Characterization of the Gene Encoding HMGS Synthase in Polygonatum sibiricum

The saponins of Polygonatum sibiricum had many pharmacological activities such as antitumor, antioxidation, and blood sugar lowering, which were synthesized by two pathways: mevalonate (MVA) and methylerythritol phosphate (MEP). 3-Hydroxy-3-methylglutaryl coenzyme A synthase (HMGS) was the key enzyme in the MVA synthesis pathway, and its expression level may affect the accumulation of saponins which were the main active ingredients of P. sibiricum. In this study, we successfully cloned HMGS1 and HMGS2 from P. sibiricum and their sequence similarity was 93.71% with 89 different sites. The multiple sequence alignment results indicated that the N-terminal sequences of HMGS were conserved. Phylogenetic analysis showed that P. sibiricum, A. officinalis, N. tazetta, D. nobile, and other relatives had a common evolutionary ancestor. The expression levels of both HMGSs and the total saponin content in different tissues revealed that HMGS expression in rhizomes was positively correlated with total saponin content. Further study of the abiotic stress effect of Methyl Jasmonate (MeJA) demonstrated that the expression of HMGS1 and HMGS2 genes was induced by MeJA, peaked at 24 h, and fell by 48 h. Our present findings would provide a blueprint for future studies of HMGS and its role in triterpenoid biosynthesis in P. sibiricum.

Cloning and functional analysis of 1-deoxy-D-xylulose-5-phosphate synthase (DXS) in Santalum album L

The commercial value of Santalum album L. lies in its aromatic heartwood and essential oil. Sesquiterpenes are the main components of sandal essential oil, and these are synthesized through the plant's mevalonate (MVA) and methylerythritol phosphate (MEP) pathways. In this study, the first key rate-limiting enzyme, 1-deoxy-D-xylulose-5-phosphate synthase (SaDXS), was investigated to provide a theoretical molecular basis for the sandalwood MEP sesquiterpene biosynthetic pathway. The biofunctions of SaDXS were also analyzed. SaDXS promoters were successfully cloned from a seven-year-old S. album tree. SaDXS1A/1B promoter activity was verified by a β-glucuronidase (GUS) assay and by analyzing cis-acting elements of the promoters, which carried light- and methyl jasmonate (MeJA)-responsive signals. In an experiment involving yellow S. album seedlings, exposure to light upregulated SaDXS1A/1B expression and increased chlorophyll and carotenoid contents when overexpressed in Arabidopsis thaliana. Analysis of the expression of SaDXS1A/1B and SaSSy, key genes of santalol biosynthesis, revealed SaDXS1A expression in all tissues whereas SaDXS1B was expressed in tissues that contained photosynthetic pigments, such as stems, leaves and flowers. Sandal seedlings exogenously treated with two hormones, MeJA and ethylene, revealed similar expression patterns for SaDXS1A/1B and SaSSy. Sandal seedlings were treated with an inhibitor of DXS, clomazone, but showed no significant changes in the contents of α-santalene, β-santalene and α-santalol between treatment and control groups. These results suggest that SaDXS1A/1B play a role in the synthesis of sandalwood sesquiterpenes, providing carbon for downstream secondary metabolites. SaDXS1A/1B also play a role in the biosynthesis of chlorophyll, carotenoids, and primary metabolites.

Identification and Functional Analysis of SabHLHs in Santalum album L

Santalum album L., a semi-parasitic evergreen tree, contains economically important essential oil, rich in sesquiterpenoids, such as (Z) α- and (Z) β-santalol. However, their transcriptional regulations are not clear. Several studies of other plants have shown that basic-helix-loop-helix (bHLH) transcription factors (TFs) were involved in participating in the biosynthesis of sesquiterpene synthase genes. Herein, bHLH TF genes with similar expression patterns and high expression levels were screened by co-expression analysis, and their full-length ORFs were obtained. These bHLH TFs were named SaMYC1, SaMYC3, SaMYC4, SaMYC5, SabHLH1, SabHLH2, SabHLH3, and SabHLH4. All eight TFs had highly conserved bHLH domains and SaMYC1, SaMYC3, SaMYC4, and SaMYC5, also had highly conserved MYC domains. It was indicated that the eight genes belonged to six subfamilies of the bHLH TF family. Among them, SaMYC1 was found in both the nucleus and the cytoplasm, while SaMYC4 was only localized in the cytoplasm and the remaining six TFs were localized in nucleus. In a yeast one-hybrid experiment, we constructed decoy vectors pAbAi-SSy1G-box, pAbAi-CYP2G-box, pAbAi-CYP3G-box, and pAbAi-CYP4G-box, which had been transformed into yeast. We also constructed pGADT7-SaMYC1 and pGADT7-SabHLH1 capture vectors and transformed them into bait strains. Our results showed that SaMYC1 could bind to the G-box of SaSSy, and the SaCYP736A167 promoter, which SaSSy proved has acted as a key enzyme in the synthesis of santalol sesquiterpenes and SaCYP450 catalyzed the ligation of santalol sesquiterpenes into terpene. We have also constructed pGreenII 62-SK-SaMYC1, pGreenII 0800-LUC-SaSSy and pGreenII 0800-LUC-SaCYP736A167 via dual-luciferase fusion expression vectors and transformed them into Nicotiana benthamiana using an Agrobacterium-mediated method. The results showed that SaMYC1 was successfully combined with SaSSy or SaCYP736A167 promoter and the LUC/REN value was 1.85- or 1.55-fold higher, respectively, than that of the control group. Therefore, we inferred that SaMYC1 could activate both SaSSy and SaCYP736A167 promoters.

Capturing acyltransferase(s) transforming final step in the biosynthesis of a major Iridoid Glycoside, (Picroside-II) in a Himalayan Medicinal Herb, Picrorhiza kurroa

BACKGROUND

Picrorhiza kurroa has been reported as an age-old ayurvedic hepato-protection to treat hepatic disorders due to the presence of iridoids such as picroside-II (P-II), picroside-I, and kutkoside. The acylation of catalpol and vanilloyl coenzyme A by acyltransferases (ATs) is critical step in P-II biosynthesis. Since accumulation of P-II occurs only in roots, rhizomes and stolons in comparison to leaves uprooting of this critically endangered herb has been the only source of this compound. Recently, we reported that P-II acylation likely happen in roots, while stolons serve as the vital P-II storage compartment. Therefore, developing an alternate engineered platform for P-II biosynthesis require identification of P-II specific AT/s.

METHODS AND RESULTS

In that direction, egg-NOG function annotated 815 ATs from de novo RNA sequencing of tissue culture based 'shoots-only' system and nursery grown shoots, roots, and stolons varying in P-II content, were cross-compared in silico to arrive at ATs sequences unique and/or common to stolons and roots. Verification for organ and accession-wise upregulation in gene expression of these ATs by qRT-PCR has shortlisted six putative 'P-II-forming' ATs. Further, six-frame translation, ab initio protein structure modelling and protein-ligand molecular docking of these ATs signified one MBOAT domain containing AT with preferential binding to the vanillic acid CoA thiol ester as well as with P-II, implying that this could be potential AT decorating final structure of P-II.

CONCLUSIONS

Organ-wise comparative transcriptome mining coupled with reverse transcription real time qRT-PCR and protein-ligand docking led to the identification of an acyltransferases, contributing to the final structure of P-II.

Full-Length Transcriptome and RNA-Seq Analyses Reveal the Mechanisms Underlying Waterlogging Tolerance in Kiwifruit (Actinidia valvata)

Actinidia valvata possesses waterlogging tolerance; however, the mechanisms underlying this trait are poorly characterized. Here, we performed a transcriptome analysis by combining single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing and investigated the physiological responses of the roots of KR5 (A. valvata, a tolerant genotype) after 0, 12, 24 and 72 h of waterlogging stress. KR5 roots responded to waterlogging stress mainly via carbohydrate and free amino acids metabolism and reactive oxygen species (ROS) scavenging pathways. Trehalose-6-phosphate synthase (TPS) activity, alcohol dehydrogenase (ADH) activity and the total free amino acid content increased significantly under waterlogging stress. The nicotinamide adenine dinucleotide-dependent glutamate synthase/alanine aminotransferase (NADH-GOGAT/AlaAT) cycle was correlated with alanine accumulation. Levels of genes encoding peroxidase (POD) and catalase (CAT) decreased and enzyme activity increased under waterlogging stress. Members of the LATERAL ORGAN BOUNDARIES (LOB), AP2/ERF-ERF, Trihelix and C3H transcription factor families were identified as potential regulators of the transcriptional response. Several hub genes were identified as key factors in the response to waterlogging stress by a weighted gene co-expression network analysis (WGCNA). Our results provide insights into the factors contributing to waterlogging tolerance in kiwifruit, providing a basis for further studies of interspecific differences in an important plant trait and for molecular breeding.

Cloning, prokaryotic expression, and purification of acetyl-CoA C-acetyltransferase from Atractylodes lancea

BACKGROUND

Cangzhu (Atractylodes lancea), a valuable and common traditional Chinese medicinal herb, is primarily used as an effective medicine with various health-promoting effects. The main pharmacological bioactive ingredients in the rhizome of A. lancea are terpenoids. Acetyl-CoA C-acetyltransferase (AACT) is the first enzyme in the terpenoid synthesis pathway and catalyzes two units of acetyl-CoA into acetoacetyl-CoA.

OBJECTIVE

The objective of the present work was to clone and identify function of AlAACT from Atractylodes lancea.

METHOD

A full-length cDNA clone of AlAACT was isolated using PCR and expressed in Escherichia coli. The expressed protein was purified using Ni-NTA agarose column using standard protocols. AlAACT was transiently expressed in N. benthamiana leaves to determine their subcellular location. The difference in growth between recombinant bacteria and control bacteria under different stresses was observed using the droplet plate experiment. Results：In this study, a full-length cDNA of AACT (AlAACT) was cloned from A. lancea, which contains a 1,227 bp open reading frame and encodes a protein with 409 amino acids. Bioinformatic and phylogenetic analysis clearly suggested that AlAACT shared high similarity with AACTs from other plants. The recombinant protein pET32a(+)/AlAACT was successfully expressed in Escherichia coli BL21(DE3) cells induced with 0.4 mM IPTG at 30°C as the optimized condition. The recombinant enzyme pET-32a-AlAACT was purified using the Ni-NTA column based on the His-tag, and the molecular weight was determined to be 62 kDa through SDS-PAGE and Western Blot analysis. The recombinant protein was eluted with 100, 300, and 500 mM imidazole; most of the protein was eluted with 300 mM imidazole. Under mannitol stress, the recombinant pET-32a-AlAACT protein showed a substantial advantage in terms of growth rates compared to the control. However, this phenomenon was directly opposite under NaCl abiotic stress. Subcellular localization showed that AlAACT localizes to the nucleus and cytoplasm. Conclusion：The expression and purification of recombinant enzyme pET-32a-AlAACT were successful, and the recombinant strain pET-32a-AlAACT in showed better growth in a drought stress. The expression of AlAACT-EGFP fusion protein revealed its localization in both nuclear and cytoplasm compartments. This study provides an important foundation for further research into the effects of terpenoid biosynthesis in A. lancea.

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.