Integrated single-cell transcriptomics and spatial metabolomics unveil cellular differentiation and ginsenosides biosynthesis in Panax root tips

Abstract

Root tips, which represent the initial stage of taproot development, serve as an ideal model for investigating plant growth and secondary metabolism. However, studies of root tips in Panax species have been limited, restricting our understanding of cell fate transitions during early root development and the cellular heterogeneity associated with ginsenosides biosynthesis. To address this gap, we conducted single-cell RNA sequencing (scRNA-seq) and spatial metabolomics analyses on the root tips of three Panax species: Panax notoginseng, Panax ginseng, and Panax quinquefolium. Our research reconstructed the developmental trajectory of the early endodermis and revealed epidermis-specific expression patterns of key enzyme genes involved in ginsenosides biosynthesis. We identified several novel transcription factors (TFs): IAA29 (which positively regulates endodermis suberization) and MYB2/MYB78 (positive regulators of ginsenosides biosynthesis), validated by dual-LUC reporter and electrophoretic mobility shift assay (EMSA). Conserved and divergent ligand-receptor interaction patterns across the three Panax species were discovered, with the FAD gene family exhibiting tissue-and species-specific expression. Cell-specific genes expression was confirmed by RNA in situ hybridization. Mass spectrometry imaging (MSI) mapped ginsenosides spatial distribution, while LC-MS/MS verified species-specific biosynthesis. This study presents a single-cell transcriptional landscape of early differentiation and cell type-specific ginsenosides accumulation in the Panax genus.

Acknowledgments

This study was financially supported by the financial support from the Major Science and Technology Special Project of Yunnan Province (Grant No. 202202AG050021) and National Science Foundation of China (Grant No. 32360151).

Author contributions

L.Y. and Z.Y. contributed equally to this work. L.Y. and Z.Y. performed data analysis, experiment section, and manuscript writing. M.L., S.W., H.W., and Q.Y. provided help for experiment section. L.H. contributed scientific advice to manuscript discussion. Y.Y. provided sample materials. X.C. and Y.L. conceived and designed the project, and then revised the manuscript.

Data availability

All raw sequencing dataset involved in this paper have been submitted to NCBI under accession number PRJNA1046239 and PRJNA1045109.

Conflict of interest statement

The authors declare no conflict of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

[1] Hou M, Wang R, Zhao S. et al. Ginsenosides in Panax genus and their biosynthesis. Acta Pharm Sin B. 2021; 11:1813-34

[2] Liu L, Xu F, Wang Y. Traditional uses, chemical diversity and biological activities of Panax L. (Araliaceae): a review. J Ethnophar-macol. 2020; 263:112792

[3] Song Y, Zhang Y, Wang X. et al. Telomere-to-telomere reference genome for Panax ginseng highlights the evolution of saponin biosynthesis. Hortic Res. 2024;11:uhae107

[4] Xia P, Huang Y, Zhu J. et al. The complete chloroplast genome sequences of 11 Panax species: providing insights for evolution and species identification. IndCropProd. 2025; 223:120160

[5] Xu C, Ma D, Ding Q. et al. PlantPhoneDB: a manually curated pan-plant database of ligand-receptor pairs infers cell-cell com-munication. Plant Biotechnol J. 2022; 20:2123-34

[6] Shulse CN, Cole BJ, Ciobanu D. et al. High-throughput single-cell transcriptome profiling of plant cell types. Cell Rep. 2019; 27:2241-2247.e4

[7] Cohen H, Fedyuk V, Wang C. et al. SUBERMAN regulates devel-opmental suberization of the Arabidopsis root endodermis. Plant J. 2020; 102:431-47

[8] Chen JT. Advances in Ginsenosides. Biomolecules. 2020; 10:681

[9] XueL, HeZ, BiX. et al. Transcriptomic profiling reveals MEP path-way contributing to ginsenoside biosynthesis in Panax ginseng. BMC Genomics. 2019; 20:383

[10] Kim J, Yun Y, Huh J. et al. Comparative transcriptome analysis on wild-simulated ginseng of different age revealed possible mechanism of ginsenoside accumulation. Plant Physiol Biochem. 2023; 201:107870

[11] Jiang Y, Liu S, Li L. et al. Transcriptome and phenotype integrated analysis identifies genes controlling ginsenoside Rb1 biosynthe-sis and reveals their interactions in the process in Panax ginseng. Int J Mol Sci. 2022; 23:14016

[12] Huang C, Li P, Yang X. et al. Integrated transcriptome and proteome analyses reveal candidate genes for ginsenoside biosynthesis in Panax japonicus C. A. Meyer. Front Plant Sci. 2023; 13:1106145

[13] Di P, Yan Y, Wang P. et al. Integrative SMRT sequencing and ginsenoside profiling analysis provide insights into the biosyn-thesis of ginsenoside in Panax quinquefolium. Chin J Nat Medicines. 2022; 20:614-26

[14] Wei G, Dong L, Yang J. et al. Integrated metabolomic and tran-scriptomic analyses revealed the distribution of saponins in Panax notoginseng. Acta Pharm Sin B. 2018; 8:458-65

[15] Wei G, Yang F, Wei F. et al. Metabolomes and transcriptomes revealed the saponin distribution in root tissues of Panax quin-quefolius and Panax notoginseng. J Ginseng Res. 2020; 44:757-69

[16] Yin R Chen,R, Xia K. et al. A single-cell transcriptome atlas reveals the trajectory of early cell fate transition during callus induction in Arabidopsis. Plant Commun 2024; 5:100941.

[17] Lopez-Anido CB, Vatén A, Smoot NK. et al. Single-cell resolution of lineage trajectories in the Arabidopsis stomatal lineage and developing leaf. Dev Cell. 2021; 56:1043-1055.e4

[18] Han E, Geng Z, Qin Y. et al. Single-cell network analysis reveals gene expression programs for Arabidopsis root development and metabolism. Plant Commun. 2024; 5:100978

[19] Zhang L, He C, Lai Y. et al. Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis. Genome Biol. 2023; 24:65

[20] Guo Y, Chen X, Li J. et al. Single-cell RNA sequencing reveals a high-resolution cell atlas of petals in Prunus mume at different flowering development stages. Hortic Res. 2024;11:uhae189

[21] Li X, Li B, Gu S. et al. Single-cell and spatial RNA sequencing reveal the spatiotemporal trajectories of fruit senescence. Nat Commun. 2024; 15:3108

[22] Qin Y, Sun M, Li W. et al. Single-cell RNA-seq reveals fate deter-mination control of an individual fibre cell initiation in cotton (Gossypium hirsutum). Plant Biotechnol. 2022; 20:2372-88

[23] Ding Y, Gao W, Qin Y. et al. Single-cell RNA landscape of the special fiber initiation process in Bombax ceiba. Plant Commun. 2023; 4:100554

[24] Xie J, Li M, Zeng J. et al. Single-cell RNA sequencing profiles of stem-differentiating xylem in poplar. Plant Biotechnol. 2022; 20: 417-9

[25] Tung CC, Kuo SC, Yang CL. et al. Single-cell transcriptomics unveils xylem cell development and evolution. Genome Biol. 2023; 24:3

[26] Cervantes-Pérez SA, Zogli P, Amini S. et al. Single-cell transcrip-tome atlases of soybean root and mature nodule reveal new regulatory programs that control the nodulation process. Plant Commun. 2024; 5:100984

[27] Frank M, Fechete LI, Tedeschi F. et al. Single-cell analysis identi-fies genes facilitating rhizobium infection in Lotus japonicus. Nat Commun. 2023; 14:7171

[28] Sun X, Feng D, Liu M. et al. Single-cell transcriptome reveals dominant subgenome expression and transcriptional response to heat stress in Chinese cabbage. Genome Biol. 2022; 23:262

[29] Gao S, Li F, Zeng Z. et al. A single-cell transcriptomic atlas reveals the cell differentiation trajectory and the response to virus invasion in swelling clove of garlic. Hortic Res. 2025;12: uhae365

[30] Zhan X, Qiu T, Zhang H. et al. Mass spectrometry imaging and single-cell transcriptional profiling reveal the tissue-specific regulation of bioactive ingredient biosynthesis in Taxus leaves. Plant commun. 2023; 4:100630

[31] Li C, Wood JC, Vu AH. et al. Single-cell multi-omics in the medic-inal plant Catharanthus roseus. Nat Chem Biol. 2023; 19:1031-41

[32] Wang Q, Wu Y, Peng A. et al. Single-cell transcriptome atlas reveals developmental trajectories and a novel metabolic path-way of catechin esters in tea leaves. Plant Biotechnol J. 2022; 20: 2089-106

[33] Zhang T, Chen Y, Wang J. A single-cell analysis of the Arabidopsis vegetative shoot apex. Dev Cell. 2021; 56:1056-1074.e8

[34] Zhang T, Chen Y, Liu Y. et al. Single-cell transcriptome atlas and chromatin accessibility landscape reveal differentiation trajec-tories in the rice root. Nat Commun. 2021; 12:2053

[35] Dai Y, Zhang S, Guan J. et al. Single-cell transcriptomic analysis of flowering regulation and vernalization in Chinese cabbage shoot apex. Hortic Res. 2024;11:uhae214

[36] Sun Y, Han Y, Sheng K. et al. Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regu-latory networks of pigment glands in Gossypium bickii. Mol Plant. 2023; 16:694-708

[37] Sun Z, Jiang S, Wang D. et al. Single-cell RNA-seq of Lotus japonicus provide insights into identification and function of root cell types of legume. J Integr Plant Biol. 2023; 65:1147-52

[38] Du J, Wang Y, Chen W. et al. High-resolution anatomical and spatial transcriptome analyses reveal two types of meristematic cell pools within the secondary vascular tissue of poplar stem. Mol Plant. 2023; 16:809-28

[39] Yuan Y, Huo Q, Zhang Z. et al. Decoding the gene regulatory network of endosperm differentiation in maize. Nat Commun. 2024; 15:34

[40] Zhou P, Chen H, Dang J. et al. Single-cell transcriptome of Nepeta tenuifolia leaves reveal differentiation trajectories in glandular trichomes. Front Plant Sci. 2022; 2022:13

[41] Zhao B, Gao Y, Ma Q. et al. Global dynamics and cytokinin partic-ipation of salt gland development trajectory in recretohalophyte Limonium bicolor. Plant Physiol. 2024; 195:2094-110

[42] Horn PJ, Chapman KD. Imaging plant metabolism in situ. JExp Bot. 2024; 75:1654-70

[43] Korsunsky I, Millard N, Fan J. et al. Fast, sensitive and accu-rate integration of single-cell data with harmony. Nat Methods. 2019; 16:1289-96

[44] Chen Y, Tong S, Jiang Y. et al. Transcriptional landscape of highly lignified poplar stems at single-cell resolution. Genome Biol. 2021; 22:319

[45] Lee J, Kim H, Park SG. et al. Brassinosteroid-BZR1/2-WAT1 mod-ule determines the high level of auxin signalling in vascu-lar cambium during wood formation. New Phytol. 2021; 230: 1503-16

[46] Conde D, Triozzi PM, Pereira WJ. et al. Single-nuclei transcrip-tome analysis of the shoot apex vascular system differentiation in Populus. Development. 2022;149:dev200632

[47] Otero S, Gildea I, Roszak P. et al. A root phloem pole cell atlas reveals common transcriptional states in protophloem-adjacent cells. Nat Plants. 2022; 8:954-70

[48] Denyer T, Ma X, Klesen S. et al. Spatiotemporal developmen-tal trajectories in the Arabidopsis root revealed using high-throughput single-cell RNA sequencing. Dev Cell. 2019; 48:840-852. e5

[49] Ryu KH, Huang L, Kang HM. et al. Single-cell RNA sequencing resolves molecular relationships among individual plant cells. Plant Physiol. 2019; 179:1444-56

[50] Shahan R, Hsu CW, Nolan TM. et al. A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants. Dev Cell. 2022; 57:543-560.e9

[51] Aichinger E, Kornet N, Friedrich T. et al. Plant stem cell niches. Annu Rev Plant Biol. 2012; 63:615-36

[52] Lu KJ, Rybel BD, Mourik HV. et al. Regulation of intercellular TARGET OF MONOPTEROS 7 protein transport in the Arabidopsis root. Development. 2018;145:dev152892

[53] Goldy C, Pedroza-Garcia JA, Breakfield N. et al. The Arabidop-sis GRAS-type SCL28 transcription factor controls the mitotic cell cycle and division plane orientation. Proc Natl Acad Sci. 2021; 118:e2005256118

[54] Kamiya T, Borghi M, Wang P. et al. The MYB36 transcription factor orchestrates casparian strip formation. Proc Natl Acad Sci USA. 2015; 112:10533-8

[55] Zhang TQ, Xu ZG, Shang GD. et al. A single-cell RNA sequencing profiles the developmental landscape of Arabidopsis root. Mol Plant. 2019; 12:648-60

[56] Fu J, Zhang X, Liu J. et al. A mechanism coordinating root elonga-tion, endodermal differentiation, redox homeostasis and stress response. Plant J. 2021; 107:1029-39

[57] Winter CM, Szekely P, Popov V. et al. SHR and SCR coordinate root patterning and growth early in the cell cycle. Nature. 2024; 626: 611-6

[58] Liberman LM, Sparks EE, Moreno-Risueno MA. et al. MYB36 reg-ulates the transition from proliferation to differentiation in the Arabidopsis root. Proc Natl Acad Sci USA. 2015; 112:12099-104

[59] Piao X, Zhang H, Kang JP. et al. Advances in saponin diversity of Panax ginseng. Molecules. 2020; 25:3452

[60] Gonneau M, Desprez T, Martin M. et al. Receptor kinase THE-SEUS1 is a rapid alkalinization factor 34 receptor in Arabidopsis. Curr Biol. 2018; 28:2452-2458.e4

[61] Tang Y, Wang L, Qu Z. et al. BSISTER transcription factors directly binds to the promoter of IAA19 and IAA29 genes to up-regulate gene expression and promote the root development. Plant Sci. 2022; 321:111324

[62] Gao H, Song W, Severing E. et al. PIF4 enhances DNA bind-ing of CDF2 to co-regulate target gene expression and pro-mote Arabidopsis hypocotyl cell elongation. Nat Plants. 2022; 8: 1082-93

[63] Xu H, Liu P, Wang C. et al. Transcriptional networks regulating suberin and lignin in endodermis link development and ABA response. Plant Physiol. 2022b; 190:1165-81

[64] Al-Khayri JM, Rashmi R, Toppo V. et al. Plant secondary metabo-lites: the weapons for biotic stress management. Metabolites. 2023; 13:716

[65] Isah T. Stress and defense responses in plant secondary metabo-lites production. Biol Res. 2019; 52:39

[66] Mahajan M, Kuiry R, Pal PK. Understanding the consequence of environmental stress for accumulation of secondary metabo-lites in medicinal and aromatic plants. J Appl Res Med Aroma. 2020; 18:100255

[67] Ma L, Liu X, Guo L. et al. Discovery of plant chemical defence mediated by a two-component system involving β-glucosidase in Panax species. Nat Commun. 2024; 15:602

[68] Chu Y, Xiao S, Su H. et al. Genome-wide characterization and analysis of bHLH transcription factors in Panax ginseng. Acta Pharm Sin B. 2018; 8:666-77

[69] Jiang Z, Tu L, Yang W. et al. The chromosome-level reference genome assembly for Panax notoginseng and insights into gin-senoside biosynthesis. Plant Commun. 2021; 2:100113

[70] Lim GH, Kim SW, Ryu J. et al. Upregulation of the MYB2 tran-scription factor is associated with increased accumulation of anthocyanin in the leaves of Dendrobium bigibbum. Int J Mol Sci. 2020; 21:5653

[71] Zhang H, Li G, Jiang M. et al. Longitudinal expression pro-files of flavonoid synthesis-related genes and their correla-tions with contents of flavonoids including anthocyanins in purplish-green aerial stems of Panax notoginseng. Acta Physiol Plant. 2023; 45:3

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

[73] Bien T, Koerfer K, Schwenzfeier J. et al. Mass spectrometry imag-ing to explore molecular heterogeneity in cell culture. Proc Natl Acad Sci USA. 2022; 119:e2114365119

[74] Kim NH, Jayakodi M, Lee SC. et al. Genome and evolution of the shade-requiring medicinal herb Panax ginseng. Plant Biotechnol. 2018; 16:1904-17

[75] Zhang H, Zhang L, Ji Y. et al. Arabidopsis SIGMA FACTOR BINDING PROTEIN1 (SIB1) and SIB2 inhibit WRKY75 function in abscisic acid-mediated leaf senescence and seed germination. JExp Bot. 2021; 73:182-96

[76] Butler A, Hoffman P, Smibert P. et al. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018; 36:411-20

[77] Yadav RK, Girke T, Pasala S. et al. Gene expression map of the Arabidopsis shoot apical meristem stem cell niche. Proc Natl Acad Sci. 2009; 106:4941-6

[78] McGinnis CS, Murrow LM, Gartner ZJ. DoubletFinder: doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst. 2019; 8:329-337.e4

[79] Qiu X, Mao Q, Tang Y. et al. Reversed graph embedding resolves complex single-cell trajectories. Nat Methods. 2017; 14:979-82

[80] Manno GL, Soldatov R, Zeisel A. et al. RNA velocity of single cells. Nature. 2018; 560:494-8

[81] Bergen V, Lange M, Peidli S. et al. Generalizing RNA velocity to transient cell states through dynamical modeling. Nat Biotechnol. 2020; 38:1408-14

[82] Yu G, Wang L, Han Y. et al. ClusterProfiler: an R package for com-paring biological themes among gene clusters. OMICS. 2012; 16: 284-7

[83] Kumar L, Futschik ME. Mfuzz: a software package for soft clus-tering of microarray data. Bioinformation. 2007; 2:5-7

[84] Zheng Y, Jiao C, Sun H. et al. iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Mol Plant. 2016; 9: 1667-70

[85] Tian F, Yang D, Meng Y. et al. PlantRegMap: charting functional regulatory maps in plants. Nucleic Acids Res. 2019;48:D1104-13

[86] Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30:2114-20

[87] Dobin A, Davis CA, Schlesinger F. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2012; 29:15-21

[88] Liao Y, Smyth GK, Shi W. featureCounts: an efficient general pur-pose program for assigning sequence reads to genomic features. Bioinformatics. 2013; 30:923-30