Single-nucleus sequencing and spatial metabolomics analysis reveal the regulatory mechanism of ginkgolic acid biosynthesis in the episperm of Ginkgo biloba

Abstract

Ginkgo biloba is a singular and relict gymnosperm indigenous to China. Its distinctive fleshy episperm is rich in unique metabolites, ginkgolic acids, which protect the developing seed from biotic stresses. The unique nature of the tissue and its metabolites has made it highly challenging to elucidate the molecular and cellular mechanisms governing ginkgolic acid biosynthesis and regulation. In this study, we performed the mass spectrometry imaging of G. biloba seed, revealing that ginkgolic acids primarily accumulate in the secretory cavities of the episperm. We constructed a single-cell expression atlas of the G. biloba episperm and identified seven cellular types: meristem cells, subepidermal cells, lignified cells, trancheid cells, parenchymal cells, secretory cavity cells, and epidermis cells. Based on the analysis of upregulated gene expression in secretory cavity cells, pseudotime analysis of cell differentiation, and gene expression trajectory analysis, we precisely identified the key enzyme-encoding genes highly associated with ginkgolic acid biosynthesis. This approach elucidated the cellular and molecular mechanisms underlying secretory cell differentiation, secretory cavity formation, and ginkgolic acid biosynthesis and accumulation in response to exogenous jasmonic acid induction. By constructing a molecular interaction network, it was determined that the GbWRKY35, encoded by Gb_25334, is the core transcription factor. We further identified the signaling proteins that interact with GbWRKY35, confirming its central positive regulatory role in ginkgolic acid biosynthesis. As a core transcription factor, GbWRKY35 regulates ginkgolic acid biosynthesis through stimulating the expression of GbAAE16. This study provides the first spatially resolved investigation into the molecular and cellular regulatory mechanisms of ginkgolic acid biosynthesis in the episperm under jasmonic acid induction.

Acknowledgements

This work was supported by the Natural Science Foundation of China (32171842), National Key Research and Development Project of China (2019YFD1100403), and Hunan Graduate Research Innovation Project (Key project, Grant No. CX20240068). We are grateful to Guangzhou Genedenovo Biotechnology Co., Ltd, for assisting in sequencing and bioinformatics analysis.

Author contributions

Y.W. and M.L. in charge of this research program were the lead investigators. Z.F., Z.Y., and Q.W. performed the experiments. B.Z., H.W., Y.W., and B.A. collected and managed samples. X.Z., H.J., and Y.X. conducted analyses of transcriptomic and metabolomic. Z.F. designed the experiments. Z.F., Z.Y., and Q.W. wrote the manuscript and contributed equally to this work.

Data availability

Raw reads have been deposited as a BioProject under accession PRJNA1214489.

Conflicts of interest statement

The authors declare no conflict of interest.

Supplementary material

Supplementary material is available at Horticulture Research online.

References

[1] Lovisetto A, Baldan B, Pavanello A, et al. Characterization of an AGAMOUS gene expressed throughout development of the fleshy fruit-like structure produced by Ginkgo biloba around its seeds. BMC Evol Biol. 2015; 15: 1-13

[2] Jin B, Xie Y, Lu Y, et al. Starch granule and protein accumulation during seed development of Ginkgo biloba L. Int Sch Res Notices. 2012; 2012: 653796

[3] Yang J, Feng Z, Liu W, et al. Exogenous hormone on episperm development and ginkgolic acid accumulation in Ginkgo biloba L. Ind Crop Prod. 2021; 160: 113140

[4] Kochibe N. Allergic substances of Ginkgo biloba. In: Hori T, Ridge RW, Tulecke W, Del Tredici P, Trémouillaux-Guiller J, Tobe H, eds. Ginkgo biloba a Global Treasure: From Biology to Medicine. Springer Japan: Tokyo, 1997, 301-7

[5] Hua Z, Wu C, Fan G, et al. The antibacterial activity and mechanism of ginkgolic acid C15:1. BMC Biotechnol. 2017; 17: 5

[6] Liu J, Chen J, Ye S, et al. Inhibitory action of ginkgolic acid against pathogenic fungi and characterisation of its inhibitory activities on. Folia Hortic. 2023; 35: 49-59

[7] Wang GX, Jiang D, Zhou Z, et al. In vivo assessment of anthelmintic efficacy of ginkgolic acids (C13:0, C15:1) on removal of Pseudodactylogyrus in European eel. Aquaculture. 2009; 297: 38-43

[8] Li Q, Sun W, Jiang Y, et al. Extraction and biodegradation of ginkgolic acids from Ginkgo biloba sarcotestae. Front Agric Sci Eng. 2017; 4: 465-72

[9] Hamad FB, Mubofu EB . Potential biological applications of bio-based anacardic acids and their derivatives. Int J Mol Sci. 2015; 16: 8569-90

[10] Cartayrade A, Bourgeois G, Balz JP, et al. The secretory apparatus of Ginkgo biloba: structure, differentiation and analysis of the secretory product. Trees. 1990; 4: 171-8

[11] Feng Z, Fan K, Yao Z, et al. Integrative metabolome and transcriptome analysis reveals GbKCS and GbMYB involved in the biosynthesis of ginkgolic acids. Ind Crop Prod. 2024; 220: 119225

[12] He B, Kun Q, Xin H, et al. Novel mechanisms for the synthesis of important secondary metabolites in Ginkgo biloba seed revealed by multi-omics data. Front Plant Sci. 2023; 14: 1196609

[13] van Beek TA . Chemical analysis of Ginkgo biloba leaves and extracts. J Chromatogr A. 2002; 967: 21-55

[14] Goldberg RB . Plants: novel developmental processes. Science. 1988; 240: 1460-7

[15] Ma S, Leng Y, Li X, et al. High spatial resolution mass spectrometry imaging for spatial metabolomics: advances, challenges, and future perspectives. TrAC Trends Anal Chem. 2023; 159: 116902

[16] Wang J, Yang E, Chaurand P, et al. Visualizing the distribution of strawberry plant metabolites at different maturity stages by MALDI-TOF imaging mass spectrometry. Food Chem. 2021; 345: 128838

[17] Zhao WH, Zhang YD, Shi YP . Visualizing the spatial distribution of endogenous molecules in wolfberry fruit at different development stages by matrix-assisted laser desorption/ionization mass spectrometry imaging. Talanta. 2021; 234: 122687

[18] Sun C, Cui L, Zhou B, et al. Visualizing the spatial distribution and alteration of metabolites in continuously cropped Salvia miltiorrhiza Bge using MALDI-MSI. J Pharm Anal. 2022; 12: 719-24

[19] Dai W, Hu Z, Xie D, et al. A novel spatial-resolution targeted metabolomics method in a single leaf of the tea plant (Camellia sinensis). Food Chem. 2020; 311: 126007

[20] 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

[21] Bawa G, Liu Z, Yu X, et al. Single-cell RNA sequencing for plant research: insights and possible benefits. Int J Mol Sci. 2022; 23: 4497

[22] Cervantes-Pérez SA, Thibivillliers S, Tennant S, et al. Challenges and perspectives in applying single nuclei RNA-seq technology in plant biology. Plant Sci. 2022; 325: 111486

[23] Song Q, Ando A, Jiang N, et al. Single-cell RNA-seq analysis reveals ploidy-dependent and cell-specific transcriptome changes in Arabidopsis female gametophytes. Genome Biol. 2020; 21: 1-18

[24] Li X, Zhang X, Gao S, et al. Single-cell RNA sequencing reveals the landscape of maize root tips and assists in identification of cell type-specific nitrate-response genes. Crop J. 2022; 10: 1589-600

[25] Yue H, Chen G, Zhang Z, et al. Single-cell transcriptome landscape elucidates the cellular and developmental responses to tomato chlorosis virus infection in tomato leaf. Plant Cell Environ. 2024; 47: 2658-72

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

[27] Song S, Wang J, Zhou J, et al. Single-cell RNA-sequencing of soybean reveals transcriptional changes and antiviral functions of GmGSTU23 and GmGSTU24 in response to soybean mosaic virus. Plant Cell Environ. 2024;

[28] Chen H, Li X, Cheng Q, et al. Single-cell landscape of long and short glandular trichomes in Nicotiana tabacum leaves. iScience. 2024; 27: 110650

[29] Li G, Zhao Q, Shi X, et al. Full-length RNA sequencing and single-nucleus sequencing deciphers programmed cell death and developmental trajectories in laticiferous canals of Decaisnea insignis fruits. Front Plant Sci. 2024; 15: 1446561

[30] Pang X, Sun J, Jia J, et al. Integration of single-cell and spatial RNA sequencing uncovers spatiotemporal transition of fruit senescence trajectory from exocarp to mesocarp in pitaya (Hylocereus undatus). Postharvest Biol Technol. 2024; 213: 112954

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

[32] Jia J, Chen E, Tian Y, et al. Fine localization and functional differentiation of exocarp cell clusters during fruit senescence revealed by single-cell and spatial transcriptomics. Postharvest Biol Technol. 2025; 219: 113223

[33] Wang L, Qi Y, Xie Y, et al. Observation on the formation and structural development of secretory cavities in integument of Ginkgo biloba L. J Yangzhou Univ (Agric Life Sci Ed). 2010; 31: 86-90

[34] Hou Z, Liu Y, Zhang M, et al. High-throughput single-cell transcriptomics reveals the female germline differentiation trajectory in Arabidopsis thaliana. Commun Biol. 2021; 4: 1149

[35] Serrano-Ron L, Perez-Garcia P, Sanchez-Corrionero A, et al. Reconstruction of lateral root formation through single-cell RNA sequencing reveals order of tissue initiation. Mol Plant. 2021; 14: 1362-78

[36] Shim Y, Lim C, Seong G, et al. The AP2/ERF transcription factor LATE FLOWERING SEMI-DWARF suppresses long-day-dependent repression of flowering. Plant Cell Environ. 2022; 45: 2446-59

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

[38] Douglas AW, Stevenson DW, Little DP . Ovule development in Ginkgo biloba L., with emphasis on the collar and nucellus. Int J Plant Sci. 2007; 168: 1207-36

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

[40] Shahan R, Hsu C, Nolan T, 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

[41] Sun G, Xia M, Li J, et al. The maize single-nucleus transcriptome comprehensively describes signaling networks governing movement and development of grass stomata. Plant Cell. 2022; 34: 1890-911

[42] Kubo H, Peeters AJM, Aarts MGM, et al. ANTHOCYANINLESS2, a Homeobox gene affecting anthocyanin distribution and root development in Arabidopsis. Plant Cell. 1999; 11: 1217-26

[43] Bai Y, Liu H, Lyu H, et al. Development of a single-cell atlas for woodland strawberry (Fragaria vesca) leaves during early Botrytis cinerea infection using single-cell RNA-seq. Hortic Res. 2022; 9: uhab055

[44] Farmer A, Thibivilliers S, Ryu KH, et al. Single-nucleus RNA and ATAC sequencing reveals the impact of chromatin accessibility on gene expression in Arabidopsis roots at the single-cell level. Mol Plant. 2021; 14: 372-83

[45] Procko C, Lee T, Borsuk A, et al. Leaf cell-specific and single-cell transcriptional profiling reveals a role for the palisade layer in UV light protection. Plant Cell. 2022; 34: 3261-79

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

[47] Torii K, Inoue K, Bekki K, et al. A guiding role of the Arabidopsis circadian clock in cell differentiation revealed by time-series single-cell RNA sequencing. Cell Rep. 2022; 40: 111059

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

[49] Liu Z, Zhou Y, Guo J, et al. Global dynamic molecular profiling of stomatal lineage cell development by single-cell RNA sequencing. Mol Plant. 2020; 13: 1178-93

[50] Kim JY, Symeonidi E, Pang T, et al. Distinct identities of leaf phloem cells revealed by single cell transcriptomics. Plant Cell. 2021; 33: 511-30

[51] Wendrich JR, Yang B, Vandamme N, et al. Vascular transcription factors guide plant epidermal responses to limiting phosphate conditions. Science. 2020; 370: eaay4970

[52] Wang Y, Mostafa S, Zeng W, et al. Function and mechanism of jasmonic acid in plant responses to abiotic and biotic stresses. Int J Mol Sci. 2021; 22: 8568

[53] Saha B, Nayak J, Srivastava R, et al. Unraveling the involvement of WRKY TFs in regulating plant disease defense signaling. Planta. 2024; 259: 7

[54] Jiang J, Ma S, Ye N, et al. WRKY transcription factors in plant responses to stresses. J Integr Plant Biol. 2017; 59: 86-101

[55] Cambra I, Martinez M, Dader B, et al. A cathepsin F-like peptidase involved in barley grain protein mobilization, HvPap-1, is modulated by its own propeptide and by cystatins. J Exp Bot. 2012; 63: 4615-29

[56] Bernoux M, Timmers T, Janueau A, et al. RD19, an Arabidopsis cysteine protease required for RRS1-R-mediated resistance, is relocalized to the nucleus by the Ralstonia solanacearum PopP2 effector. Plant Cell. 2008; 20: 2252-64

[57] Cerri M, Reale L . Anatomical traits of the principal fruits: an overview. Sci Hortic. 2020; 270: 109390

[58] Burns JK, Achor DS, Echeverria E . Ultrastructural studies on the ontogeny of grapefruit juice vesicles (Citrus paradisi Macf. cv. Star Ruby). Int J Plant Sci. 1992; 153: 14-25

[59] Izquierdo L, Sendra J . Citrus fruits| composition and characterization. Encyclopedia of Food Sciences and Nutrition. 2003; 1335-1341

[60] Beck CB . An Introduction to Plant Structure and Development: Plant Anatomy for the Twenty-First Century. Cambridge, England, UK: Cambridge University Press; 2010:

[61] Fahn A . Secretory tissues in vascular plants. New Phytol. 1988; 108: 229-57

[62] Alain C, Bourgeois G, Balz J . The secretory apparatus of Ginkgo biloba: structure, differentiation and analysis of the secretions. Trees. 1990; 4: 171-8

[63] Ameele RJ . Developmental anatomy of secretory cavities in the microsporophylls of Ginkgo biloba L. Am J Bot. 1980; 67: 912-7

[64] Li B, Neumann E, Ge J, et al. Interrogation of spatial metabolome of Ginkgo biloba with high-resolution matrix-assisted laser desorption/ionization and laser desorption/ionization mass spectrometry imaging. Plant Cell Environ. 2018; 41: 2693-703

[65] Yao X, Zhou G, Tang Y, et al. Comparison of Total Ginkgolic Acids Content in Leaf of Fruit-Utilization Type of Ginkgo biloba with Different Tree Ages and from Different Locations. J Plant Resour Environ. 2012; 21: 108-10

[66] Liang SJ, Wu H, Lun X, et al. Secretory cavity development and its relationship with the accumulation of essential oil in fruits of Citrus medica L. var. sarcodactylis (Noot.) Swingle. J Integr Plant Biol. 2006; 48: 573-83

[67] Bosabalidis A, Tsekos I . Ultrastructural studies on the secretory cavities of Citrus deliciosa ten. II. Development of the essential oil-accumulating central space of the gland and process of active secretion. Protoplasma. 1982; 112: 63-70

[68] Liu X, Liu B, Jiang D, et al. The accumulation and composition of essential oil in kumquat peel. Sci Hortic. 2019; 252: 121-9

[69] Horisberger JD . Recent insights into the structure and mechanism of the sodium pump. Physiology. 2004; 19: 377-87

[70] Malcolm SB . Cardenolide-mediated interactions between plants and herbivores. In: Rosenthal GA, Berenbaum MR, eds. Herbivores: Their Interactions with Secondary Plant Metabolites. San Diego, California, USA: Elsevier, 1991, 251-96

[71] Brower L, Seiber J, Nelson C, et al. Plant-determined variation in cardenolide content and thin-layer chromatography profiles of monarch butterflies, Danaus plexippus reared on milkweed plants in California: 3. Asclepias californica. J Chem Ecol. 1984; 10: 1823-57

[72] van Loon JJ, Schoonhoven LM . Specialist deterrent chemoreceptors enable Pieris caterpillars to discriminate between chemically different deterrents. In: Simpson SJ, Mordue (Luntz) AJ, Hardie J (eds.), Proceedings of the 10th International Symposium on Insect-Plant Relationships. 1999, pp. 29-35, Dordrecht, Springer Netherlands

[73] Zumajo-Cardona C, Frangos S, Stevenson DW . Seed anatomy and development in Cycads and Ginkgo, keys for understanding the evolution of seeds. Flora. 2021; 285: 151951

[74] Jin J, Lu P, Xu Y, et al. PCMDB: a curated and comprehensive resource of plant cell markers. Nucleic Acids Res. 2022; 50: D1448-55

[75] Wagner G, Wang E, Shepherd R . New approaches for studying and exploiting an old protuberance, the plant trichome. Ann Bot. 2004; 93: 3-11

[76] Ishida T, Kurata T, Okada K, et al. A genetic regulatory network in the development of trichomes and root hairs. Annu Rev Plant Biol. 2008; 59: 365-86

[77] Zhou YF, Mao SL, Li SF, et al. Programmed cell death: a mechanism for the lysigenous formation of secretory cavities in leaves of Dictamnus dasycarpus. Plant Sci. 2014; 225: 147-60

[78] Ciccarelli D, Garbari F, Pagni AM . The flower of Myrtus communis (Myrtaceae): secretory structures, unicellular papillae, and their ecological role. Flora. 2008; 203: 85-93

[79] Seo SY, Kim YJ, Kim J, et al. Phytosphingosine induces systemic acquired resistance through activation of sphingosine kinase. Plant Direct. 2021; 5: e351

[80] Glenz R, Schmalhaus D, Krischke M, et al. Elevated levels of phosphorylated sphingobases do not antagonize sphingobase-or Fumonisin B1-induced plant cell death. Plant Cell Physiol. 2019; 60: 1109-19

[81] Magnin-Robert M, Bourse D, Markham J, et al. Modifications of sphingolipid content affect tolerance to hemibiotrophic and necrotrophic pathogens by modulating plant defense responses in Arabidopsis. Plant Physiol. 2015; 169: 2255-74

[82] Werker E . Seed Anatomy. Berlin and Stuttgart, Germany: Gebrüder Borntraeger, 1997

[83] Fangren P, Hongyan G, Mingzhuo H, et al. The ultrastructure characteristics of secretory cavities associated with the secretory products of Ginkgo biloba. Am J Plant Sci. 2012; 3: 102-9

[84] Qi J, Zhou G, Yang L, et al. The chloroplast-localized phospholipases D α4 and α5 regulate herbivore-induced direct and indirect defenses in rice. Plant Physiol. 2011; 157: 1987-99

[85] Nakata M, Mitsuda N, Herde M, et al. A bHLH-type transcription factor, ABA-INDUCIBLE BHLH-TYPE TRANSCRIPTION FACTOR/JA-ASSOCIATED MYC2-LIKE1, acts as a repressor to negatively regulate jasmonate signaling in Arabidopsis. Plant Cell. 2013; 25: 1641-56

[86] Hao BZ, Wu JL . Laticifer differentiation in Hevea brasiliensis: induction by exogenous jasmonic acid and linolenic acid. Ann Bot. 2000; 85: 37-43

[87] Shockey J, Browse J . Genome-level and biochemical diversity of the acyl-activating enzyme superfamily in plants. Plant J. 2011; 66: 143-60

[88] Carlberg C, Velleuer E, Molnár F. Transcription factors and signal transduction. In: Carlberg C, Velleuer E, Molnár F (eds.), Molecular Medicine. Springer International Publishing: Cham, 2020, 35-57

[89] Kendirgi F, Rexer DJ, Alcázar-Román AR, et al. Interaction between the shuttling mRNA export factor Gle1 and the nucleoporin hCG1: a conserved mechanism in the export of Hsp70 mRNA. Mol Biol Cell. 2005; 16: 4304-15

[90] Sung P, Guzder SN, Prakash L, et al. Reconstitution of TFIIH and requirement of its DNA helicase subunits, Rad3 and Rad25, in the incision step of nucleotide excision repair (*). J Biol Chem. 1996; 271: 10821-6

[91] Franklin KA, Lee S, Patel D, et al. Phytochrome-interacting factor 4 (PIF4) regulates auxin biosynthesis at high temperature. Proc Natl Acad Sci. 2011; 108: 20231-5

[92] Hornitschek P, Kohnen M, Lorrain S, et al. Phytochrome interacting factors 4 and 5 control seedling growth in changing light conditions by directly controlling auxin signaling. Plant J. 2012; 71: 699-711

[93] Javed T, Gao SJ . WRKY transcription factors in plant defense. Trends Genet. 2023; 39: 787-801

[94] Jiang W, Chen J, Duan X, et al. Comparative transcriptome profiling reveals two WRKY transcription factors positively regulating polysaccharide biosynthesis in Polygonatum cyrtonema. Int J Mol Sci. 2023; 24: 12943

[95] Chen Z, Xu M, Quan C, et al. Genome-wide identification and expression analysis of the WRKY gene family reveal essential roles in abiotic stress responses and polysaccharides and flavonoids biosynthesis in Platostoma palustre (Blume) AJ Paton. BMC Plant Biol. 2024; 24: 1122

[96] Yan Y, Wang P, He C, et al. MeWRKY20 and its interacting and activating autophagy-related protein 8 (MeATG8) regulate plant disease resistance in cassava. Biochem Biophys Res Commun. 2017; 494: 20-6

[97] Song I, Hong S, Huh SU . Identification and expression analysis of the Solanum tuberosum StATG8 family associated with the WRKY transcription factor. Plants. 2022; 11: 2858

[98] Huh SU . Interaction of potato autophagy-related StATG8 family proteins with pathogen effector and WRKY transcription factor in the nucleus. Microorganisms. 2025; 13: 1589

[99] Wani SH, Anand S, Singh B, et al. WRKY transcription factors and plant defense responses: latest discoveries and future prospects. Plant Cell Rep. 2021; 40: 1071-85

[100] Wang D, Wei L, Liu T, et al. Suppression of ETI by PTI priming to balance plant growth and defense through an MPK3/MPK6-WRKYs-PP2Cs module. Mol Plant. 2023; 16: 903-18

[101] Deslandes L, Oliver J, Theulieres F, et al. Resistance to Ralstonia solanacearum in Arabidopsis thaliana is conferred by the recessive RRS1-R gene, a member of a novel family of resistance genes. Proc Natl Acad Sci. 2002; 99: 2404-9

[102] Chen B, OuYang C, Tian Z, et al. A high resolution atmospheric pressure matrix-assisted laser desorption/ionization-quadrupole-orbitrap MS platform enables in situ analysis of biomolecules by multi-mode ionization and acquisition. Anal Chim Acta. 2018; 1007: 16-25

[103] Keller C, Maeda J, Jayaraman D, et al. Comparison of vacuum MALDI and AP-MALDI platforms for the mass spectrometry imaging of metabolites involved in salt stress in Medicago truncatula. Front Plant Sci. 2018; 9: 1238

[104] Sunaga-Franze DY, Muino J, Breauning C, et al. Single-nucleus RNA sequencing of plant tissues using a nanowell-based system. Plant J. 2021; 108: 859-69

[105] Lun AT, Riesenfeld S, Andrews T, et al. EmptyDrops: distinguishing cells from empty droplets in droplet-based single-cell RNA sequencing data. Genome Biol. 2019; 20: 1-9

[106] 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

[107] 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

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

[109] Chung NC, Storey JD . Statistical significance of variables driving systematic variation in high-dimensional data. Bioinformatics. 2014; 31: 545-54

[110] McInnes L, Healy J, Melville J . Umap: uniform manifold approximation and projection for dimension reduction. J Open Source Softw. 2018; 3: 861

[111] Rotta R, Noack A . Multilevel local search algorithms for modularity clustering. J Exp Algorithmics. 2011; 16: 2-1

[112] Camp JG, Sekine K, Gerber T, et al. Multilineage communication regulates human liver bud development from pluripotency. Nature. 2017; 546: 533-8

[113] Trapnell C, Cacchiarelli D, Grimsby J, et al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol. 2014; 32: 381-6

[114] Cui Y, Su Y, Bian J, et al. Single-nucleus RNA and ATAC sequencing analyses provide molecular insights into early pod development of peanut fruit. Plant Commun. 2024; 5: 100979

[115] Want EJ, Masson P, Michopoulos F, et al. Global metabolic profiling of animal and human tissues via UPLC-MS. Nat Protoc. 2013; 8: 17-32

[116] Luo P, Dai W, Yin P, et al. MRM-ion pair finder: a systematic approach to transform nontargeted mode to pseudo-targeted mode for metabolomics study based on liquid chromatography-mass spectrometry. Anal Chem. 2015; 87: 5050-5

[117] Feng Z, Yang T, Li M, et al. Identification and cloning of GbMADS6, a SOC1 homolog gene involved in floral development in Ginkgo biloba. J Plant Biochem Biotechnol. 2021; 30: 554-63

[118] Schmittgen TD, Livak KJ . Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008; 3: 1101-8