Identification of sesquiterpene synthase genes in the genome of Aquilaria sinensis and characterization of an α-humulene synthase

doi:10.1007/s11676-022-01542-w

Abstract

Abstract:

Sesquiterpenes are the major pharmacodynamic components of agarwood, a precious traditional Chinese medicine obtained from the resinous portions of Aquilaria sinensis trees that form in response to environmental stressors. To characterize the sesquiterpene synthases responsible for sesquiterpene production in A. sinensis, a bioinformatics analysis of the genome of A. sinensis identified six new terpene synthase genes, and 16 sesquiterpene synthase genes were identified as type TPS-a in a phylogenetic analysis. The expression patterns for eight of the sesquiterpene synthase genes after treatment with various hormones or hydrogen peroxide were analyzed by real-time quantitative PCR. The results suggest that 100 μM methyl jasmonate, ethephon, ( ±)-abscisic acid or hydrogen peroxide could be effective short-term effectors to increase the expression of sesquiterpene synthase genes, while 1 mM methyl salicylate may have long-term effects on increasing the expression of specific sesquiterpene synthase genes (e.g., As-SesTPS, AsVS, AsTPS12 and AsTPS29). The expression changes in these genes under various conditions reflected their specific roles during abiotic or biotic stresses. Heterologous expression of a novel A. sinensis sesquiterpene synthase gene, AsTPS2, in Escherichia coli produced a major humulene product, so AsTPS2 is renamed AsHS1. AsHS1 is different from ASS1, AsSesTPS, and AsVS, for mainly producing α-humulene. Based on the predicted space conformation of the AsHS1 model, the small ligand molecule may bind to the free amino acid by hydrogen bonding for the catalytic function of the enzyme, while the substrate farnesyl diphosphate (FPP) probably binds to the free amino acid on one side of the RxR motif. Arg450, Asp453, Asp454, Thr457, and Glu461 from the NSE/DTE motif and D307 and D311 from the DDxxD motif were found to form a polar interaction with two Mg²⁺ clusters by docking. The Mg²⁺-bound DDxxD and NSE/DTE motifs and the free RXR motif are jointly directed into the catalytic pocket of AsHS1. Comparison of the tertiary structural models of AsHS1 with ASS1 showed that they differed in structures in several positions, such as surrounding the secondary catalytic pocket, which may lead to differences in catalytic products. Based on the results, biosynthetic pathways for specific sesquiterpenes such as α-humulene in A. sinensis are proposed. This study provides novel insights into the functions of the sesquiterpene synthases of A. sinensis and enriches knowledge on agarwood formation.

Key words: Agarwood, Aquilaria sinensis, Sesquiterpene synthase, Humulene

Jiadong Ran, Yuan Li, Xin Wen, Xin Geng, Xupeng Si, Liping Zhang, Yimian Ma, Zheng Zhang. Identification of sesquiterpene synthase genes in the genome of Aquilaria sinensis and characterization of an α-humulene synthase[J]. JOURNAL OF FORESTRY RESEARCH, 2023, 34(4): 1117-1131.

Jiadong Ran, Yuan Li, Xin Wen, Xin Geng, Xupeng Si, Liping Zhang, Yimian Ma, Zheng Zhang. [J]. 林业研究(英文版), 2023, 34(4): 1117-1131.

References 82

1	Alquezar B, Rodriguez A, de la Pena M, Pena L. Genomic analysis of terpene synthase family and functional characterization of seven sesquiterpene synthases from Citrus sinensis. Front Plant Sci, 2017, 8: 1481, DOI
2	Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, Studholme DJ, Yeats C, Eddy SR. The Pfam protein families database. Nucleic Acids Res, 2004, 32: 138D-141, DOI
3	Bleeker PM, Spyropoulou EA, Diergaarde PJ, Volpin H, De Both MT, Zerbe P, Bohlmann J, Falara V, Matsuba Y, Pichersky E, Haring MA, Schuurink RC. RNA-seq discovery, functional characterization, and comparison of sesquiterpene synthases from Solanum lycopersicum and Solanum habrochaites trichomes. Plant Mol Biol, 2011, 77(4–5): 323-336, DOI
4	Bohlmann J, Meyer-Gauen G, Croteau R. Plant terpenoid synthases: molecular biology and phylogenetic analysis. P Natl Acad Sci USA, 1998, 95(8): 4126-4133, DOI
5	Booth JK, Yuen MMS, Jancsik S, Madilao LL, Page JE, Bohlmann J. Terpene synthases and terpene variation in Cannabis sativa. Plant Physiol, 2020, 184(1): 130-147, DOI
6	Chen F, Tholl D, Bohlmann J, Pichersky E. The family of terpene synthases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. Plant J, 2011, 66(1): 212-229, DOI
7	Chen HQ, Wei JH, Yang JS, Zhang Z, Yang Y, Gao ZH, Sui C, Gong B. Chemical constituents of agarwood originating from the endemic genus Aquilaria plants. Chem Biodivers, 2012, 9(2): 236-250, DOI
8	Chen X, Sui C, Liu YY, Yang Y, Liu PW, Zhang Z, Wei JH. Agarwood formation induced by fermentation liquid of Lasiodiplodia theobromae, the dominating fungus in wounded wood of Aquilaria sinensis. Curr Microbiol, 2017, 74(4): 460-468, DOI
9	Chen X, Yang W, Zhang LQ, Wu XM, Cheng T, Li GL. Genome-wide identification, functional and evolutionary analysis of terpene synthases in pineapple. Comput Biol Chem, 2017, 70: 40-48, DOI
10	Chen CJ, Chen H, Zhang Y, Thomas HR, Frank MH, He YH, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant, 2020, 13(8): 1194-1202, DOI
11	Chhipa H, Chowdhary K, Kaushik N. Artificial production of agarwood oil in Aquilaria sp. by fungi: a review. Phytochem Rev, 2017, 16: 835-860, DOI
12	Chia JSM, Omar Farouk AA, Mohamad AS, Sulaiman MR, Perimal EK. Zerumbone alleviates chronic constriction injury-induced allodynia and hyperalgesia through serotonin 5-HT receptors. Biomed Pharmacother, 2016, 83: 1303-1310, DOI
13	The Pharmacopoeia Commission and others (2020) The Pharmacopoeia of People’s Republic of China (I). China pharmaceutical science and technology press, Beijing
14	Cites (2004) Amendments to Appendix I and II of CITES. Proceedings of thirteenth meeting of the conference of the parties. Bangkok, Thailand
15	Das A, Begum K, Akhtar S, Ahmed R, Kulkarni R, Banu S. Genome-wide detection and classification of terpene synthase genes in Aquilaria agallochum. Physiol Mol Biol Plants, 2021, 27(8): 1711-1729, DOI
16	Ding XP, Mei JJ, Peng SQ, Mei WL, Dai HF. Cloning and expression analysis of sesquiterpene synthase gene AsVS in Aquilaria sinensis. Mol Plant Breeding, 2019, 17(7): 2152-2158, DOI
17	Ding XP, Mei WL, Lin Q, Wang H, Wang J, Peng SQ, Li HL, Zhu JH, Li W, Wang P, Chen HQ, Dong WH, Guo D, Cai CH, Huang SZ, Cui P, Dai HF. Genome sequence of the agarwood tree Aquilaria sinensis (Lour.) Spreng: the first chromosome-level draft genome in the Thymelaeceae family. Gigascience, 2020, 9(3): giaa013, DOI
18	Dudareva N, Cseke L, Blanc VM, Pichersky E. Evolution of floral scent in Clarkia: novel patterns of S-linalool synthase gene expression in the C. breweri flower. Plant Cell, 1996, 8(7): 1137-1148, DOI
19	Durairaj J, Di Girolamo A, Bouwmeester HJ, de Ridder D, Beekwilder J, van Dijk AD. An analysis of characterized plant sesquiterpene synthases. Phytochemistry, 2019, 158: 157-165, DOI
20	El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, Qureshi M, Richardson LJ, Salazar GA, Smart A, Sonnhammer ELL, Hirsh L, Paladin L, Piovesan D, Tosatto SCE, Finn RD. The Pfam protein families database in 2019. Nucleic Acids Res, 2019, 47(D1): D427-D432, DOI
21	Emanuelsson O, Nielsen H, von Heijne G. ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites. Protein Sci, 1999, 8(5): 978-984, DOI
22	Emanuelsson O, Brunak S, von Heijne G, Nielsen H. Locating proteins in the cell using TargetP. SignalP and Related Tools Nat Protoc, 2007, 2(4): 953-971, DOI
23	Falara V, Alba JM, Kant MR, Schuurink RC, Pichersky E. Geranyllinalool synthases in solanaceae and other angiosperms constitute an ancient branch of diterpene synthases involved in the synthesis of defensive compounds. Plant Physiol, 2014, 166(1): 428-441, DOI
24	Fernandes ES, Passos GF, Medeiros R, da Cunha FM, Ferreira J, Campos MM, Pianowski LF, Calixto JB. Anti-inflammatory effects of compounds alpha-humulene and (-)-trans-caryophyllene isolated from the essential oil of Cordia verbenacea. Eur J Pharmacol, 2007, 569(3): 228-236, DOI
25	Gao Y, Honzatko RB, Peters RJ. Terpenoid synthase structures: a so far incomplete view of complex catalysis. Nat Prod Rep, 2012, 29(10): 1153-1175, DOI
26	Gao ZH, Wei JH, Yang Y, Zhang Z, Zhao WT. Selection and validation of reference genes for studying stress-related agarwood formation of Aquilaria sinensis. Plant Cell Rep, 2012, 31(9): 1759-1768, DOI
27	Gao M, Han XM, Sun Y, Chen HJ, Yang Y, Liu YY, Meng H, Gao ZH, Xu YH, Zhang Z, Han JP. Overview of sesquiterpenes and chromones of agarwood originating from four main species of the genus Aquilaria. RSC Adv, 2019, 9(8): 4113-4130, DOI
28	Halim VA, Vess A, Scheel D, Rosahl S. The role of salicylic acid and jasmonic acid in pathogen defence. Plant Biol, 2008, 8: 307-313, DOI
29	Herde M, Gartner K, Kollner TG, Fode B, Boland W, Gershenzon J, Gatz C, Tholl D. Identification and regulation of TPS04/GES, an Arabidopsis geranyllinalool synthase catalyzing the first step in the formation of the insect-induced volatile C16-homoterpene TMTT. Plant Cell, 2008, 20(4): 1152-1168, DOI
30	Huang H, Kuo YW, Chuang YC, Yang YP, Huang LM, Jeng MF, Chen WH, Chen HH. Terpene synthase-b and terpene synthase-e/f genes produce monoterpenes for Phalaenopsis bellina floral scent. Front Plant Sci, 2021, 12: 700958, DOI
31	Jones CG, Moniodis J, Zulak KG, Scaffidi A, Plummer JA, Ghisalberti EL, Barbour EL, Bohlmann J. Sandalwood fragrance biosynthesis involves sesquiterpene synthases of both the terpene synthase (TPS)-a and TPS-b subfamilies, including santalene synthases. J Biol Chem, 2011, 286(20): 17445-17454, DOI
32	Kapustin Y, Souvorov A, Tatusova T, Lipman D. Splign: algorithms for computing spliced alignments with identification of paralogs. Biol Direct, 2008, 3: 20, DOI
33	Keeling CI, Weisshaar S, Ralph SG, Jancsik S, Hamberger B, Dullat HK, Bohlmann J. Transcriptome mining, functional characterization, and phylogeny of a large terpene synthase gene family in spruce (Picea spp.). BMC Plant Biol, 2011, 11: 43, DOI
34	Koo HJ, Vickery CR, Xu Y, Louie GV, O'Maille PE, Bowman M, Nartey CM, Burkart MD, Noel JP. Biosynthetic potential of sesquiterpene synthases: product profiles of Egyptian Henbane premnaspirodiene synthase and related mutants. J Antibiot (tokyo), 2016, 69(7): 524-533, DOI
35	Kristanti AN, Tanjung M, Aminah NS. Review: secondary metabolites of Aquilaria, a Thymelaeaceae genus. Mini Rev Org Chem, 2018, 15(1): 36-55, DOI
36	Kumeta Y, Ito M. Characterization of δ-guaiene synthases from cultured cells of Aquilaria, responsible for the formation of the sesquiterpenes in agarwood. Plant Physiol, 2010, 154(4): 1998-2007, DOI
37	Kumeta Y, Ito M. Characterization of α-humulene synthases responsible for the production of sesquiterpenes induced by methyl jasmonate in Aquilaria cell culture. J Nat Med, 2016, 70(3): 452-459, DOI
38	Kurosaki F, Hirohashi S, Katoh T, Taura F, Lee JB. Cloning and characterization of δ-guaiene synthase genes encoding a sesquiterpene cyclase from Aquilaria microcarpa cell cultures. Am J Plant Sci, 2015, 6(16): 2603-2611, DOI
39	Lee S, Chappell J. Biochemical and genomic characterization of terpene synthases in Magnolia grandiflora. Plant Physiol, 2008, 147(3): 1017-1033, DOI
40	Lee JB, Hirohashi S, Yamamura Y, Taura F, Kurosaki F. Induction, cloning and functional expression of a sesquiterpene biosynthetic enzyme, δ-guaiene synthase, of Aquilaria microcarpa cell cultures. Nat Prod Commun, 2014, 9(9): 1231-1235, DOI
41	Li ZH, Peng JY, Wen X, Guo HW. Ethylene-Insensitive3 is a senescence-associated gene that accelerates age-dependent leaf senescence by directly repressing miR164 transcription in Arabidopsis. Plant Cell, 2013, 25: 3311-3128, DOI
42	Li RS, Zhu JH, Guo D, Li HL, Wang Y, Ding XP, Mei WL, Chen ZB, Dai HF, Peng SQ. Genome-wide identification and expression analysis of terpene synthase gene family in Aquilaria sinensis. Plant Physiol Bioch, 2021, 164: 185-194, DOI
43	Li W, Chen HQ, Wang H, Mei WL, Dai HF. Natural products in agarwood and Aquilaria plants: chemistry, biological activities and biosynthesis. Nat Prod Rep, 2021, 38(3): 528-565, DOI
44	Liang L, Guo QM, Zhang Z, Xu YH, Han XM, Liu J. Cloning, prokaryotic expression, and functional identification of a sesquiterpene synthase gene (AsSS4) from Aquilaria sinensis. Yao Xue Xue Bao, 2014, 49(12): 1724-1729, DOI
45	Liu YY, Chen HQ, Yang Y, Zhang Z, Wei JH, Meng H, Chen WP, Feng JD, Gan BC, Chen XY, Gao ZH, Huang JQ, Chen B, Chen HJ. Whole-tree agarwood-inducing technique: an efficient novel technique for producing high-quality agarwood in cultivated Aquilaria sinensis trees. Molecules, 2013, 18(3): 3086-3106, DOI
46	Liu JY, Huang F, Wang X, Zhang M, Zheng R, Wang J, Yu DY. Genome-wide analysis of terpene synthases in soybean: functional characterization of GmTPS3. Gene, 2014, 544(1): 83-92, DOI
47	Liu J, Xu YH, Zhang Z, Wei JH. Hydrogen peroxide promotes programmed cell death and salicylic acid accumulation during the induced production of sesquiterpenes in cultured cell suspensions of Aquilaria sinensis. Funct Plant Biol, 2015, 42(4): 337-346, DOI
48	Lv FF, Li SS, Feng J, Liu PW, Gao ZH, Yang Y, Xu YH, Wei JH. Hydrogen peroxide burst triggers accumulation of jasmonates and salicylic acid inducing sesquiterpene biosynthesis in wounded Aquilaria sinesis. J Plant Physiol, 2019, 234–235: 167-175, DOI
49	Martin DM, Aubourg S, Schouwey MB, Daviet L, Schalk M, Toub O, Lund ST, Bohlmann J. Functional annotation, genome organization and phylogeny of the grapevine (Vitis vinifera) terpene synthase gene family based on genome assembly, FLcDNA cloning, and enzyme assays. BMC Plant Biol, 2010, 10: 226, DOI
50	Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R. Corrigendum to: IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol, 2020, 37(8): 2461, DOI
51	Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol, 2020, 37(5): 1530-1534, DOI
52	Monggoot S, Popluechai S, Gentekaki E, Pripdeevech P. Fungal endophytes: an alternative source for production of volatile compounds from agarwood oil of Aquilaria subintegra. Microb Ecol, 2017, 74(1): 54-61, DOI
53	Nieuwenhuizen NJ, Wang MY, Matich AJ, Green SA, Chen X, Yauk YK, Beuning LL, Nagegowda DA, Dudareva N, Atkinson RG. Two terpene synthases are responsible for the major sesquiterpenes emitted from the flowers of kiwifruit (Actinidia deliciosa). J Exp Bot, 2009, 60(11): 3203-3219, DOI
54	Nong W, Law STS, Wong AYP, Baril T, Swale T, Chu LM, Hayward A, Lau DTW, Hui JHL. Chromosomal-level reference genome of the incense tree Aquilaria sinensis. Mol Ecol Resour, 2020, 20(4): 971-979, DOI
55	Ogita S, Lee JB, Kurosaki F, Kato Y. The biosynthetic activities of primary and secondary metabolites in suspension cultures of Aquilaria microcarpa. Nat Prod Commun, 2015, 10(5): 779-782, DOI
56	Rahman HS, Rasedee A, Abdul AB, Zeenathul NA, Othman HH, Yeap SK, How CW, Hafiza WA. Zerumbone-loaded nanostructured lipid carrier induces G2/M cell cycle arrest and apoptosis via mitochondrial pathway in a human lymphoblastic leukemia cell line. Int J Nanomed, 2014, 9: 527-538, DOI
57	Rasool S, Mohamed R. Mohamed R. Understanding agarwood formation and its challenge. Agarwood: science behind the fragrance, 2016 Berlin Springer, DOI
58	Rizvi SM, Shakil S, Haneef M. A simple click by click protocol to perform docking: AutoDock 4.2 made easy for non-bioinformaticians. Excli J, 2013, 12: 831-857, DOI
59	Sen S, Dehingia M, Talukdar NC, Khan M. Chemometric analysis reveals links in the formation of fragrant bio-molecules during agarwood (Aquilaria malaccensis) and fungal interactions. Sci Rep, 2017, 7: 44406, DOI
60	Siah CH, Namasivayam P, Mohamed R. Transcriptome reveals senescing callus tissue of Aquilaria malaccensis, an endangered tropical tree, triggers similar response as wounding with respect to terpenoid biosynthesis. Tree Genet Genomes, 2016, 12(2): 1-10, DOI
61	Singh S, Thulasiram HV, Sengupta D, Kulkarni K. Dynamic coupling analysis on plant sesquiterpene synthases provides leads for the identification of product specificity determinants. Biochem Bioph Res Co, 2021, 536: 107-114, DOI
62	Starks CM, Back K, Chappell J, Noel JP. Structural basis for cyclic terpene biosynthesis by tobacco 5-epi-aristolochene synthase. Science, 1997, 277(5333): 1815-1820, DOI
63	Subasinghe S, Hitihamu H, Fernando K. Use of two fungal species to induce agarwood resin formation in Gyrinops walla. J Forestry Res, 2019, 30: 345-350, DOI
64	Sun PW, Xu YH, Yu CC, Lv FF, Tang XL, Gao ZH, Zhang Z, Wang H, Liu Y, Wei JH. WRKY44 represses expression of the wound-induced sesquiterpene biosynthetic gene ASS1 in Aquilaria sinensis. J Exp Bot, 2020, 71(3): 1128-1138, DOI
65	Tan CS, Isa NM, Ismail I, Zainal Z. Agarwood induction: current developments and future perspectives. Front Plant Sci, 2019, 10: 122, DOI
66	Wang GD, Tian L, Aziz N, Broun P, Dai XB, He J, King A, Zhao PX, Dixon RA. Terpene biosynthesis in glandular trichomes of hop. Plant Physiol, 2008, 148(3): 1254-1266, DOI
67	Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, Heer FT, de Beer TAP, Rempfer C, Bordoli L, Lepore R, Schwede T. Swiss-Model: homology modelling of protein structures and complexes. Nucleic Acids Res, 2018, 46(W1): W296-W303, DOI
68	Wendt KU, Schulz GE. Isoprenoid biosynthesis: manifold chemistry catalyzed by similar enzymes. Structure, 1998, 6(2): 127-133, DOI
69	Wu ZQ, Liu S, Li JF, Li MC, Du HF, Qi LK, Lin L. Analysis of gene expression and quality of agarwood using Agar-bit in Aquilaria sinensis. J Trop for Sci, 2017, 29: 380-388, DOI
70	Xu YH, Zhang Z, Wang MX, Wei JH, Chen HJ, Gao ZH, Sui C, Luo HM, Zhang XL, Yang Y, Meng H, Li WL. Identification of genes related to agarwood formation: transcriptome analysis of healthy and wounded tissues of Aquilaria sinensis. BMC Genomics, 2013, 14: 227, DOI
71	Xu YH, Liao YC, Zhang Z, Liu J, Sun PW, Gao ZH, Sui C, Wei JH. Jasmonic acid is a crucial signal transducer in heat shock induced sesquiterpene formation in Aquilaria sinensis. Sci Rep, 2016, 6: 21843, DOI
72	Yan H, Ren MY, Wang ZX, Feng SJ, Li S, Cheng Y, Hu CX, Gao SQ, Zhang GQ. Zerumbone inhibits melanoma cell proliferation and migration by altering mitochondrial functions. Oncol Lett, 2017, 13(4): 2397-2402, DOI
73	Yan TT, Yang S, Chen Y, Wang Q, Li GY. Chemical profiles of cultivated agarwood induced by different techniques. Molecules, 2019, 24(10): 1990, DOI
74	Yang CQ, Wu XM, Ruan JX, Hu WL, Mao YB, Chen XY, Wang LJ. Isolation and characterization of terpene synthases in cotton (Gossypium hirsutum). Phytochemistry, 2013, 96: 46-56, DOI
75	Ye W, He X, Wu HQ, Wang L, Zhang WM, Fan YF, Li HH, Liu TM, Gao XX. Identification and characterization of a novel sesquiterpene synthase from Aquilaria sinensis: an important gene for agarwood formation. Int J Biol Macromol, 2018, 108: 884-892, DOI
76	Yin YF, Jiao LC, Dong MY, Jiang XM, Zhang SJ. Mohamed R. Wood resources, identification, and utilization of Agarwood in China. Agarwood, 2016 Singapore Springer 21-38, DOI
77	Yu F, Okamto S, Nakasone K, Adachi K, Matsuda S, Harada H, Misawa N, Utsumi R. Molecular cloning and functional characterization of α-humulene synthase, a possible key enzyme of zerumbone biosynthesis in shampoo ginger (Zingiber zerumbet Smith). Planta, 2008, 227(6): 1291-1299, DOI
78	Zhang Z, Yang Y, Wei JH, Meng H, Sui C, Chen HQ. Advances in studies on mechanism of agarwood formation in Aquilaria sinensis and its hypothesis of agarwood formation induced by defense response. Chin Tradit Herb Drugs, 2010, 41: 156-160 (in Chinese)
79	Zhang Z, Zhang XL, Yang Y, Wei JH, Meng H, Gao ZH, Xu YH. Hydrogen peroxide induces vessel occlusions and stimulates sesquiterpenes accumulation in stems of Aquilaria sinensis. Plant Growth Regul, 2014, 72(1): 81-87, DOI
80	Zhou F, Pichersky E. The complete functional characterisation of the terpene synthase family in tomato. New Phytol, 2020, 226(5): 1341-1360, DOI
81	Zhou HC, Shamala LF, Yi XK, Yan Z, Wei S. Analysis of terpene synthase family genes in Camellia sinensis with an emphasis on abiotic stress conditions. Sci Rep, 2020, 10(1): 933, DOI
82	Zhou ZY, Xian JC, Wei WK, Xu C, Yang JF, Zhan RT, Ma DM. Volatile metabolic profiling and functional characterization of four terpene synthases reveal terpenoid diversity in different tissues of Chrysanthemum indicum L. Phytochemistry, 2021, 185: 112687, DOI