Single-cell RNA sequencing reveals a high-resolution cell atlas of petals in Prunus mume at different flowering development stages

Abstract

Prunus mume (mei), a traditional ornamental plant in China, is renowned for its fragrant flowers, primarily emitted by its petals. However, the cell types of mei petals and where floral volatile synthesis occurs are rarely reported. The study used single-cell RNA sequencing to characterize the gene expression landscape in petals of P. mume ‘Fenhong Zhusha’ at budding stage (BS) and full-blooming stage (FS). Six major cell types of petals were identified: epidermal cells (ECs), parenchyma cells (PCs), xylem parenchyma cells, phloem parenchyma cells, xylem vessels and fibers, and sieve elements and companion cells complex. Cell-specific marker genes in each cell type were provided. Floral volatiles from mei petals were measured at four flowering development stages, and their emissions increased from BS to FS, and decreased at the withering stage. Fifty-eight differentially expressed genes (DEGs) in benzenoid/phenylpropanoid pathway were screened using bulk RNA-seq data. Twenty-eight DEGs expression increased from BS to FS, indicating that they might play roles in floral volatile synthesis in P. mume, among which PmBAHD3 would participate in benzyl acetate synthesis. ScRNA-seq data showed that 27 DEGs mentioned above were expressed variously in different cell types. In situ hybridization confirmed that PmPAL2, PmCAD1, PmBAHD3, 5, and PmEGS1 involved in floral volatile synthesis in mei petals are mainly expressed in EC, PC, and most vascular tissues, consistent with scRNA-seq data. The result indicates that benzyl acetate and eugenol, the characteristic volatiles in mei, are mostly synthesized in these cell types. The first petal single-cell atlas was constructed, offering new insights into the molecular mechanism of floral volatile synthesis.

Acknowledgements

The work supported by the National Natural Science Foundation of China (No. 32301659), the Natural Science Basic Research Program of Shaanxi (Program No. 2023-JC-QN-0188), and the Chinese Universities Scientific Fund (No. 2452021098).

Author contributions

T.Z. designed the experiments; J.L., Q.W., S.Z., N.L., T.Z., X.C., and Y.G. performed the research; T.Z., Y.G., and X.C. wrote the paper. Y.Z. revised the manuscript. All authors discussed the results, reviewed the article, and approved the final article.

Data availability

The data supporting the findings of this study are available within the paper figures and the supplementary data.

Conflict of interests

The authors declare no conflicts of interest.

Supplementary information

Supplementary data is available at Horticulture Research online.

References

[1] Dudareva N, Pichersky E. Floral scent metabolic pathways: their regulation and evolution. Biol Floral Scent. 2006; 1:55-78

[2] Dudareva N, Pichersky E. Biochemical and molecular genetic aspects of floral scents. Plant Physiol. 2000; 122:627-33

[3] Dudareva N, Cseke L, Blanc VM. et al. Evolution of floral scent in clarkia: novel patterns of S-linalool synthase gene expression in the C. breweri flower. Plant Cell. 1996; 8:1137-48

[4] Kolosova N, Sherman D, Karlson D. et al. Cellular and subcellular localization of S-adenosyl-L-methionine: benzoic acid carboxyl methyltransferase, the enzyme responsible for biosynthesis of the volatile ester methylbenzoate in snapdragon flowers. Plant Physiol. 2001; 126:956-64

[5] Scalliet G, et al.Lionnet C, le Bechec M. Role of petal-specific orcinol O-methyltransferases in the evolution of rose scent. Plant Physiol. 2006; 140:18-29

[6] Muhlemann JK, Klempien A, Dudareva N. Floral volatiles: from biosynthesis to function. Plant Cell Environ. 2014; 37:1936-49

[7] Bao F, Ding A, Zhang T. et al. Expansion of PmBEAT genes in the Prunus mume genome induces characteristic floral scent production. Hortic Res-England. 2019; 6:24

[8] Zhang T, Bao F, Yang Y. et al. A comparative analysis of floral scent compounds in intraspecific cultivars of Prunus mume with different corolla colours. Molecules. 2020; 25:145

[9] Yuan X, Ma K, Zhang M. et al. Integration of transcriptome and methylome analyses provides insight into the pathway of floral scent biosynthesis in Prunus mume. Front Genet. 2021; 12:779557

[10] Bao F, Zhang T, Ding A. et al. Metabolic, enzymatic activity, and transcriptomic analysis reveals the mechanism underlying the lack of characteristic floral scent in apricot Mei varieties. Front Plant Sci. 2020; 11:574982

[11] Xiujun W, Zhenqi S, Yujing T. et al. Comparative transcriptome analysis linked to key volatiles reveals molecular mechanisms of aroma compound biosynthesis in Prunus mume. BMC Plant Biol. 2022; 22:1-19

[12] Zhao K, Yang W, Zhou Y. et al. Comparative transcriptome reveals benzenoid biosynthesis regulation as inducer of floral scent in the woody plant Prunus mume. Front Plant Sci. 2017; 8:319

[13] Kang M, Choi Y, Kim H. et al. Single-cell RNA-sequencing of Nicotiana attenuata corolla cells reveals the biosynthetic pathway of a floral scent. New Phytol. 2022; 234:527-44

[14] Ke M, Elshenawy B, Sheldon H. et al. Single cell RNA-sequencing: a powerful yet still challenging technology to study cellular heterogeneity. BioEssays. 2022; 44:2200084

[15] Wang Y, Huan Q, Li K. et al. Single-cell transcriptome atlas of the leaf and root of rice seedlings. J Genet Genomics. 2021; 48:881-98

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

[17] 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-England. 2022; 9:uhab055

[18] Li-Beisson Y, Pollard M, Sauveplane V. et al. Nanoridges that characterize the surface morphology of flowers require the synthesis of cutin polyester. Proc Natl Acad Sci. 2009; 106:22008-13

[19] DeBono A, Yeats TH, JKC R. et al. Arabidopsis LTPG is a glycosylphosphatidylinositol-anchored lipid transfer protein required for export of lipids to the plant surface. Plant Cell. 2009; 21:1230-8

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

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

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

[23] Liu Z, Wang J, Zhou Y. et al. Identification of novel regulators required for early development of vein pattern in the cotyledons by single-cell RNA-sequencing. Plant J. 2022; 110:7-22

[24] Ranocha P, Denancé N, Vanholme R. et al. Walls are thin 1 (WAT1), an Arabidopsis homolog of Medicago truncatula NODULIN21, is a tonoplast-localized protein required for secondary wall formation in fibers. Plant J. 2010; 63:469-83

[25] Dinant S, Clark AM, Zhu Y. et al. Diversity of the superfamily of phloem lectins (phloem protein 2) in angiosperms. Plant Physiol. 2003; 131:114-28

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

[27] Yu D, Liu X, Cui Y. et al. Comparative transcriptome combined with morpho-physiological analyses revealed candidate genes potentially for differential cold tolerance in two contrasting apricot (Prunus armeniaca L.) cultivars. Trees. 2020; 34:1205-17

[28] Han S, Liu H, Yan M. et al. Differential gene expression in leaf tissues between mutant and wild-type genotypes response to late leaf spot in peanut (Arachis hypogaea L.). PLoS One. 2017; 12:e0183428

[29] Qiu L, Zheng T, Liu W. et al. Integration of transcriptome and metabolome reveals the formation mechanism of red stem in Prunus mume. Front Plant Sci. 2022; 13:884883

[30] Stefanovic A, Ribot C, Rouached H. et al. Members of the PHO1 gene family show limited functional redundancy in phosphate transfer to the shoot, and are regulated by phosphate deficiency via distinct pathways. Plant J. 2007; 50:982-94

[31] Zhang T, Huo T, Ding A. et al. Genome-wide identification, characterization, expression and enzyme activity analysis of coniferyl alcohol acetyltransferase genes involved in eugenol biosynthesis in Prunus mume. PLoS One. 2019; 14:e0223974

[32] Pyke KA, Page AM. Plastid ontogeny during petal development in Arabidopsis. Plant Physiol. 1998; 116:797-803

[33] Zhou Q, Shi M, Zhang H. et al. Comparative study of the petal structure and fragrance components of the Nymphaea hybrid, a precious water lily. Molecules. 2022; 27:408

[34] Effmert U, Große J, Röse US. et al. Volatile composition, emission pattern, and localization of floral scent emission in Mirabilis jalapa (Nyctaginaceae). Am J Bot. 2005; 92:2-12

[35] Sulborska-Ró˙zycka A, Weryszko-Chmielewska E, Polak B. et al. Secretory products in petals of Centaurea cyanus L. flowers: a histochemistry, ultrastructure, and phytochemical study of volatile compounds. Molecules. 2022; 27:1371

[36] Hao R, du D, Wang T. et al. A comparative analysis of characteristic floral scent compounds in Prunus mume and related species. Biosci Biotechnol Biochem. 2014; 78:1640-7

[37] Koeduka T, Louie GV, Orlova I. et al. The multiple phenylpropene synthases in both Clarkia breweri and Petunia hybrida represent two distinct protein lineages. Plant J. 2008; 54:362-74

[38] Hao R, Yang S, Zhang Z. et al. Identification and specific expression patterns in flower organs of ABCG genes related to floral scent from Prunus mume. Sci Hortic. 2021; 288:110218

[39] Kim D, Paggi JM, Park C. et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37:907-15

[40] Pertea M, Kim D, Pertea GM. et al. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and ballgown. Nat Protoc. 2016; 11:1650-67

[41] Zang Y, Pei Y, Cong X. et al. Single-cell RNA-sequencing profiles reveal the developmental landscape of the Manihot esculenta Crantz leaves. Plant Physiol. 2023; 194:456-74