Abstract
Plant glandular trichomes (GTs) are characterised by their ability to synthesise and store abundant secondary metabolites of significant biological importance. The medicinal plant Artemisia argyi exhibits a dense covering of GTs, yet the precise metabolic compositions and the developmental molecular dynamics of A. argyi GTs remain insufficiently characterised. Using liquid chromatography-mass spectrometry (LC-MS) and gas chromatography-mass spectrometry (GC-MS), we identified 969 differentially accumulated metabolites (DAMs) between GTs and non-glandular trichomes (NGTs) of A. argyi. GTs were significantly enriched in diverse metabolites, including terpenoids, flavonoids and fatty acyls, among which sesquiterpenoids were the most abundant terpenoid subclass in DAMs. To elucidate the molecular basis of GT secondary metabolite biosynthesis and development, we generated a single-cell transcriptomic atlas of A. argyi leaves, annotating cell populations of mesophyll cells (MCs), epidermal cells (ECs), vascular cells (VCs), stomatal guard cells (GCs) and GTs. Pseudotime trajectory analysis uncovered the continuous developmental trajectory of GTs and identified several candidate transcription factors (TFs) potentially involved in GT development. Leveraging the single-cell atlas, we constructed cell-type-specific co-expression networks for sesquiterpene biosynthesis genes. Within the GT-specific expression module, we identified β-caryophyllene synthase AarTPS77, while the EC-specific AarTPS52 was found to function as β-farnesene synthase. Furthermore, we functionally characterised AarTPS95 and AarTPS96, which catalyse the formation of germacrene A and 12 additional sesquiterpenoids. These findings provide insights into the molecular basis of trichome development and secondary metabolite accumulation in A. argyi, laying the foundation for improving the quality of medicinal materials.