Integrative chromosome-level genomics and metabolomics uncover regulatory networks linking monoterpenoid biosynthesis and glandular trichome formation in Mosla chinensis

Abstract

Xiangru, with Mosla chinensis (Mc, 2n = 18) and its considered cultivar M. chinensis ‘Jiangxiangru’ (McJ, 2n = 18) as original plants, is an annual herb of the Lamiaceae family, and is widely used as medicinal and edible plant due to its spleen strengthening function. However, absence of genomic resource impedes in-depth research towards Xiangru. In this study, the morphological characteristics and volatile organic compounds (VOC) contents of Mc and McJ were analyzed, showing higher trichome density and monoterpenoid accumulation obtained in Mc, whereas McJ possessed higher biomass. We assembled high-quality Mc, McJ, and their adulterant Mosla soochowensis (2n = 18) genomes of 426.1, 408.8, and 412.8 Mb, respectively, containing the repeat sequences of 57.17%, 56.33%, and 55.83%. Comparative genomics analysis indicated Mosla radiated ~13.3 Mya, supporting McJ initially as a natural naturally formed resource. Five monoterpene synthase genes were identified through comparative transcriptome and were responsible for catalyzing production of diversified monoterpene skeleton, in which TPS1 mediated formation of γ-terpinene, accompanied by CYP71D179 and SDR2, leading to the final production of carvacrol and thymol. We further explored correlation between monoterpenoids biosynthesis and trichome development, indicating MIXTA and WIN1 jointly regulate both trichome formation and VOC accumulation by directly binding promoters of TPS1 and CYP71D179, respectively. Our study fills vacancy of genus Mosla genomes, improving the biosynthetic and regulatory mechanism of volatile compounds in aromatic Traditional Chinese Medicine, also offering novel targets for quality-directed breeding in Xiangru.

Acknowledgements

This research was supported by the National Key Research and Development Program of China (2023YFC3503901), Natural Science Foundation of China (82404802, 82404785), National Administration of Traditional Chinese Medicine Science and Technology Project (GZY-KJS-2023-032), Scientific and Technological Innovation Project of China Academy of Chinese Medical Sciences (CI2024C005YN, CI2023E002), Fundamental Research Funds for the Central public welfare research institutes (ZZ15-YQ-060, ZZXT202417), earmarked fund for CARS (CARS-21), and the key project at the central government level: The ability establishment of sustainable use for valuable Chinese medicine resources (2060302).

Authors contributions

C.S.-L., J.Y., and L.H. designed the research. M.Y. and H.Z. performed the function identification, regulatory network construction, VOC determination, and data analysis. C.Y.-L. performed the genome assembly and evolutionary analysis. D.J., X.F., and S.Q. conducted transcriptional regulation validation. X.W. and C.C. collected the plant materials and analyzed the data. D.J., G.R., and X.F. helped with transcriptome analysis. Y.Z. helped with functional verification in vivo. X.L., Q.L., Y.H., and X.H. helped with sample treatment and data analysis. Z.B. and B.T. helped with the data visualization. M.Y., C.Y.-L., and X.W. wrote the manuscript. C.S.-L. and H.Z. revised the manuscript. All authors read and approved the manuscript.

Data availability

The raw sequences and genome assemblies are deposited in China National Center for Bioinformation (CNCB, https://www.cncb.ac.cn/) under BioProject ID PRJCA032919, while the annotations are available at figshare (https://doi.org/10.6084/m9.figshare.29635061.v1).

Conflicts of interest statement

The authors declare no competing interests.

Supplementary material

Supplementary material is available at Horticulture Research online.

References

[1] Hu S, Cui F, Yuan C. Materia medica textual research of herba Jiangxiangru and studies on its resources. J Jiangxi Univ Tradit Chin Med. 1994; 6:31-4

[2] Wang Y, Wu X, Chen Y. et al. Phylogenomic analyses revealed widely occurring hybridization events across Elsholtzieae (Lamiaceae). Mol Phylogenet Evol. 2024; 198:108112

[3] Crocoll C, Asbach J, Novak J. et al. Terpene synthases of oregano ( and their roles in the path-way and regulation of terpene biosynthesis. Plant Mol Biol. 2010; 73:587-603

[4] Lima AS, Schimmel J, Lukas B. et al. Genomic characteriza-tion, molecular cloning and expression analysis of two ter-pene synthases from Thymus caespititius (Lamiaceae). Planta. 2013; 238:191-204

[5] Krause ST, Liao P, Crocoll C. et al. The biosynthesis of thymol, carvacrol, and thymohydroquinone in Lamiaceae proceeds via cytochrome P450s and a short-chain dehydrogenase. Proc Natl Acad SciUSA. 2021; 118:e2110092118

[6] Sun M, Zhang Y, Zhu L. et al. Chromosome-level assembly and analysis of the Thymus genome provide insights into glandular secretory trichome formation and monoterpenoid biosynthesis in thyme. Plant Commun. 2022; 3:100413

[7] Liu C, Srividya N, Parrish AN. et al. Morphology of glandular trichomes of Japanese catnip (Schizonepeta tenuifolia Bri-quet) and developmental dynamics of their secretory activity. Phytochemistry. 2018; 150:23-30

[8] Amrehn E, Aschenbrenner AK, Heller A. et al. Localization of sesquiterpene lactone biosynthesis in cells of capitate glan-dular trichomes of Helianthus annuus (Asteraceae). Proto-plasma. 2016; 253:447-55

[9] Sagawa T, Ikeda H, Hiraoka T. et al. Study of rosemary peltate glandular trichomes using combined morphologi-cal and chemical approach. Food Sci Technol Res. 2013; 19: 491-5

[10] Chang J, Yu T, Yang Q. et al. Hair, encoding a single C2H2 zinc-finger protein, regulates multicellular trichome formation in tomato. Plant J. 2018; 96:90-102

[11] Jae-In C, Seong-Min K, Heejin K. et al. SlHair2 regulates the initiation and elongation of type I trichomes on tomato leaves and stems. Plant Cell Physiol. 2021; 62:1446-59

[12] Yang C, Li H, Zhang J. et al. A regulatory gene induces trichome formation and embryo lethality in tomato. Proc Natl Acad Sci USA. 2011; 108:11836-41

[13] Ewas M, Gao Y, Ali F. et al. RNA-seq reveals mechanisms of SlMX1 for enhanced carotenoids and terpenoids accu-mulation along with stress resistance in tomato. Sci Bull. 2017; 62:476-85

[14] Yu X, Chen G, Tang B. et al. The jasmonate ZIM-domain pro-tein gene SlJAZ2 regulates plant morphology and accelerates flower initiation in Solanum lycopersicum plants. Plant Sci. 2018; 267:65-73

[15] Liu H, He W, Yao X. et al. The light- and jasmonic acid-induced AaMYB108-like positive regulates the initiation of glandu-lar secretory trichome in Artemisia annua L. Int J Mol Sci. 2023; 24:12929

[16] Yan T, Li L, Xie L. et al. A novel HD-ZIP IV/MIXTA complex pro-motes glandular trichome initiation and cuticle development in Artemisia annua. New Phytol. 2018; 218:567-78

[17] He Y, Fu X, Li L. et al. AaSPL9 affects glandular trichomes initia-tion by positively regulating expression of AaHD1 in Artemisia annua L. Plant Sci. 2022; 317:111172

[18] Chen T, Liu H, Li Y. et al. AaSEPALLATA1 integrates jasmonate and light-regulated glandular secretory trichome initiation in Artemisia annua. Plant Physiol. 2023; 192:1483-97

[19] XieL YanT,LiL. et al. The WRKY transcription factor AaGSW2 promotes glandular trichome initiation in Artemisia annua. J Exp Bot. 2020; 72:1691-701

[20] Matías-Hernández L, Jiang W, Yang K. et al. AaMYB1 and its orthologue AtMYB61 affect terpene metabolism and trichome development in Artemisia annua and Arabidopsis thaliana. Plant J. 2017; 90:520-34

[21] Lv Z, Li JX, Qiu S. et al. The transcription factors TLR1 and TLR2 negatively regulate trichome density and artemisinin levels in Artemisia annua. J Integr Plant Biol. 2022; 64:1212-28

[22] Yu L, Zhang Y, Ding Q. et al. The SlMYC1-TOR module regulates trichome formation and terpene biosynthesis in tomatoes (Solanum lycopersicum L.). J Plant Growth Regul. 2024; 43:3282-94

[23] Wang Y, Fu X, Xie L. et al. Stress associated protein 1 regulates the development of glandular trichomes in Artemisia annua. Plant Cell Tissue Org. 2019; 139:249-59

[24] Zheng F, Cui L, Li C. et al. Hair interacts with SlZFP8-like to regulate the initiation and elongation of trichomes by modu-lating SlZFP6 expression in tomato. J Exp Bot. 2021; 73:228-44

[25] Qi X, Chen Z, Yu X. et al. Characterisation of the Mentha canadensis R2R3-MYB transcription factor gene McMIXTA and its involvement in peltate glandular trichome development. BMC Plant Biol. 2022; 22:219

[26] Zhou P, Dang J, Shao Y. et al. Cloning of StHD1 and StHD8 from Schizonepeta tenuifolia and function of regulating glan-dular trichome development. Chin J Chin Mater Med. 2022; 47: 5838-48

[27] Liu Z, Zhang C, He J. et al. plantGIR: a genomic database of plants. Hortic Res. 2024;11:uhae342

[28] Gao S, Gao Y, Xiong C. et al. The tomato B-type cyclin gene, SlCycB2, plays key roles in reproductive organ development, trichome initiation, terpenoids biosynthesis and Prodenia litura defense. Plant Sci. 2017; 262:103-14

[29] Berman P, Haro LAD, Jozwiak A. et al. Parallel evolution of cannabinoid biosynthesis. Nat Plants. 2023; 9:817-831

[30] Yang S, Chu H, Wang Y. et al. Analysis of monoterpene biosyn-thesis and functional TPSs of Perilla frutescens based on tran-scriptome and metabolome. Med Plant Biol. 2024; 3:e017

[31] Bae SH, Lee MH, Lee JH. et al. The genome of the Korean island-originated Perilla citriodora ’Jeju17’ sheds light on its environmental adaptation and fatty acid and lipid production pathways. Genes. 2023; 14:17

[32] Zhou Y, Tang N, Huang L. et al. Effects of salt stress on plant growth, antioxidant capacity, glandular trichome density, and volatile exudates of Schizonepeta tenuifolia Briq. Int J Mol Sci. 2018; 19:252

[33] Schilmiller AL, Last RL, Pichersky E. Harnessing plant tri-chome biochemistry for the production of useful compounds. Plant J. 2010; 54:702-11

[34] Shen Q, Zhang L, Liao Z. et al. The genome of Artemisia annua provides insight into the evolution of Asteraceae family and artemisinin biosynthesis. Mol Plant. 2018; 11:776-88

[35] Tiwari P. Recent advances and challenges in trichome research and essential oil biosynthesis in Mentha arvensis L. Ind Crop Prod. 2016; 82:141-8

[36] Zhang Y, Wang D, Li H. et al. Formation mechanism of glan-dular trichomes involved in the synthesis and storage of ter-penoids in lavender. BMC Plant Biol. 2023; 23:307

[37] Feng S, Liu Z, Chen H. et al. PHGD: an integrative and user-friendly database for plant hormone-related genes. iMeta. 2024; 3:e164

[38] Huang H, Qiao H, Ma X. et al. Roles of jasmonates in tomato growth, development and defense. Veg Res. 2023; 3:14

[39] Yuan S, Li Q, Shen H. et al. Advances in the regulatory mech-anisms of multicellular trichome formation and its secondary metabolite synthesis in vegetable crops. Veg Res. 2023; 3:24

[40] Jiao H, Hua Z, Zhou J. et al. Genome-wide analysis of Panax MADS-box genes reveals role of PgMADS41 and PgMADS44 in modulation of root development and ginsenoside synthesis. Int J Biol Macromol. 2023; 233:123648

[41] Yu M, Ma C, Tai B. et al. Unveiling the regulatory mech- anisms of nodules development and quality formation in Panax notoginseng using multi-omics and MALDI-MSI. J Adv Sci. 2025; 69:463-75

[42] Zhang H, Chen H, Hou Z. et al. Overexpression of Ath-MIR160b increased the biomass while reduced the content of tan-shinones in Salvia miltiorrhiza hairy roots by targeting ARFs genes. Plant Cell Tiss Org. 2020; 142:327-38

[43] Hua B, Chang J, Wu M. et al. Mediation of JA signalling in glandular trichomes by the woolly/SlMYC1 regulatory mod-ule improves pest resistance in tomato. Plant Biotechnol J. 2021; 19:375-93

[44] Guo Z, Hao K, Lv Z. et al. Profiling of phytohormone-specific microRNAs and characterization of the miR160-ARF1 module involved in glandular trichome development and artemisinin biosynthesis in Artemisia annua. Plant Biotechnol J. 2023; 21:591-605

[45] Kind T, Wohlgemuth G, Lee DY. et al. FiehnLib: mass spec-tral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal Chem. 2009; 81:10038-48

[46] Chen F, Tholl D, Tholl D. et al. The family of terpene syn-thases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the king-dom. Plant J. 2011; 66:212-29

[47] Fu XQ, Peng BW, Hassani D. et al. AaWRKY9 contributes to light- and jasmonate-mediated to regulate the biosynthesis of artemisinin in Artemisia annua. New Phytol. 2021; 231:1858-74