Developmental landscape and asymmetric gene expression in the leaf vasculature of Brassica rapa revealed by single-cell transcriptome

Abstract

Leaf vasculature not only acts as a channel for nutrients and signaling information but also influences leaf morphology. It consists of several distinct cell types with specialized functions. Cell type-specific characterizations based on single-cell RNA sequencing technology could aid in understanding the identities of vascular tissues and their roles in leaf morphogenesis in Brassica rapa. Here, we generated a single-cell transcriptome landscape of the Chinese cabbage leaf vasculature. A total of 12 cell clusters covering seven known cell types were identified. Different vascular cell types were characterized by distinct identities. The xylem parenchyma and companion cells exhibited an active expression pattern of amino acid metabolism genes. Tracheary elements and sieve elements were enriched in many genes related to cell wall biosynthesis, and the phloem parenchyma was enriched in many sugar transporter-encoding genes. Pseudo-time analyses revealed the developmental trajectories of the xylem and phloem and the potential roles of auxin and ethylene in xylem development. Furthermore, we identified key candidate regulators along the differentiation trajectory of the sieve elements and companion cells. Most of the homoeologous genes in the syntenic triads from the three subgenomes showed asymmetric gene expression patterns in different vascular cell types. Collectively, our study revealed that Chinese cabbage leaf vasculature cells had highly heterogeneous transcriptomes, providing new insights into the complex processes of leaf vasculature development in B. rapa leafy vegetables and other Brassica crops.

Acknowledgements

This work was supported by the State Key Laboratory of Vegetable Biobreeding, the National Natural Science Foundation of China (32402585; 32102393; 32472730), the Science and Technology Program of Henan Province (242102111143).

Author contributions

X.W. and J.W. designed the research. X.G., J.Y., and Y.Z. prepared materials. X.G. and J.Y. performed data analysis. X.G., J.W., and X.W. wrote the manuscript.

Data availability

The data used in this study are included in the article.

Conflict of interest statement

The authors declare no conflicts of interest.

References

[1] Xu B, Ohtani M, Yamaguchi M. et al. Contribution of NAC tran-scription factors to plant adaptation to land. Science. 2014; 343: 1505-8

[2] Lucas WJ, Groover A, Lichtenberger R. et al. The plant vascular system: evolution, development and functions. J Integr Plant Biol. 2013; 55:294-388

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

[4] Haritatos E, Medville R, Turgeon R. Minor vein structure and sugar transport in Arabidopsis thaliana. Planta. 2000; 211:105-11

[5] Scarpella E, Barkoulas M, Tsiantis M. Control of leaf and vein development by auxin. Cold Spring Harb Perspect Biol. 2010;2: a001511

[6] Yamaguchi T, Nukazuka A, Tsukaya H. Leaf adaxial-abaxial polarity specification and lamina outgrowth: evolution and development. Plant Cell Physiol. 2012; 53:1180-94

[7] Cheng F, Sun R, Hou X. et al. Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea. Nat Genet. 2016; 48:1218-24

[8] Yue X, Su T, Xin X. et al. The adaxial/abaxial patterning of auxin and auxin gene in leaf veins functions in leafy head formation of Chinese cabbage. Front Plant Sci. 2022; 13:918112

[9] Beuchat J, Scacchi E, Tarkowska D. et al. BRX promotes Arabidop-sis shoot growth. New Phytol. 2010; 188:23-9

[10] Truernit E, Bauby H, Belcram K. et al. OCTOPUS, a polarly localised membrane-associated protein, regulates phloem dif-ferentiation entry in Arabidopsis thaliana. Development. 2012; 139: 1306-15

[11] Zhang X, Ma W, Liu M. et al. OCTOPUS regulates BIN2 to con-trol leaf curvature in Chinese cabbage. Proc Natl Acad Sci USA. 2022; 119:e2208978119

[12] Chan AC, Khan D, Girard IJ. et al. Tissue-specific laser microdis-section of the Brassica napus funiculus improves gene discov-ery and spatial identification of biological processes. JExp Bot. 2016; 67:3561-71

[13] Jansen L, Hollunder J, Roberts I. et al. Comparative transcrip-tomics as a tool for the identification of root branching genes in maize. Plant Biotechnol J. 2013; 11:1092-102

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

[15] Tian C, Wang Y, Yu H. et al. A gene expression map of shoot domains reveals regulatory mechanisms. Nat Commun. 2019; 10:141

[16] Rich-Griffin C, Stechemesser A, Finch J. et al. Single-cell tran-scriptomics: a high-resolution avenue for plant functional genomics. Trends Plant Sci. 2020; 25:186-97

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

[18] Liu Q, Liang Z, Feng D. et al. Transcriptional landscape of rice roots at the single-cell resolution. Mol Plant. 2021; 14:384-94

[19] Xu X, Crow M, Rice BR. et al. Single-cell RNA sequencing of developing maize ears facilitates functional analysis and trait candidate gene discovery. Dev Cell. 2021; 56:557-568.e6

[20] Nelms B, Walbot V. Defining the developmental program leading to meiosis in maize. Science. 2019; 364:52-6

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

[22] Chen X, Ru Y, Takahashi H. et al. Single-cell transcriptomic anal-ysis of pea shoot development and cell-type-specific responses to boron deficiency. Plant J. 2023; 117:302-22

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

[24] Li H, Dai X, Huang X. et al. Single-cell RNA sequencing reveals a high-resolution cell atlas of xylem in Populus. J Integr Plant Biol. 2021; 63:1906-21

[25] 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-9

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

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

[28] Wang, Wang H, Wang J. et al. The genome of the mesopolyploid crop species Brassica rapa. Nat Genet. 2011; 43:1035-9

[29] Guo X, Liang J, Lin R. et al. Single-cell transcriptome reveals differentiation between adaxial and abaxial mesophyll cells in Brassica rapa. Plant Biotechnol J. 2022; 20:2233-5

[30] Zhong R, Kandasamy MK, Ye ZH. XND 1 regulates secondary wall deposition in xylem vessels through the inhibition of VND functions. Plant Cell Physiol. 2021; 62:53-65

[31] Nguyen CT, Kurenda A, Stolz S. et al. Identification of cell popula-tions necessary for leaf-to-leaf electrical signaling in a wounded plant. Proc Natl Acad Sci USA. 2018; 115:10178-83

[32] Zhao C, Craig JC, Petzold HE. et al.The xylem and phloem transcriptomes from secondary tissues of the Arabidopsis root-hypocotyl. Plant Physiol. 2005; 138:803-18

[33] Funk V, Kositsup B, Zhao C. et al. The Arabidopsis xylem peptidase XCP1 is a tracheary element vacuolar protein that may be a papain ortholog. Plant Physiol. 2002; 128:84-94

[34] Chen LQ, Qu XQ, Hou BH. et al. Sucrose efflux mediated by SWEET proteins as a key step for phloem transport. Science. 2012; 335:207-11

[35] Tian H, Baxter IR, Lahner B. et al. Arabidopsis NPCC6/NaKR1 is a phloem mobile metal binding protein necessary for phloem function and root meristem maintenance. Plant Cell. 2010; 22: 3963-79

[36] Cayla T, Batailler B, le Hir R. et al. Live imaging of compan-ion cells and sieve elements in Arabidopsis leaves. PLoS One. 2015; 10:e0118122

[37] Froelich DR, Mullendore DL, Jensen KH. et al. Phloem ultra-structure and pressure flow: sieve-element-occlusion-related agglomerations do not affect translocation. Plant Cell. 2011; 23: 4428-45

[38] Anstead JA, Froelich DR, Knoblauch M. et al. Arabidopsis P-protein filament formation requires both AtSEOR1 and AtSEOR2. Plant Cell Physiol. 2012; 53:1033-42

[39] Barratt DH, Kölling K, Graf A. et al. Callose synthase GSL7 is nec-essary for normal phloem transport and inflorescence growth in Arabidopsis. Plant Physiol. 2011; 155:328-41

[40] Shi D, Lebovka I, Lopez-Salmeron V. et al. Bifacial cambium stem cells generate xylem and phloem during radial plant growth. Development. 2019;146:dev171355

[41] Schuetz M, Smith R, Ellis B. Xylem tissue specification, pattern-ing, and differentiation mechanisms. JExp Bot. 2013; 64:11-31

[42] Wallner ES, López-Salmerón V, Belevich I. et al. Strigolactone-and karrikin-independent SMXL proteins are central regulators of phloem formation. Curr Biol. 2017; 27:1241-7

[43] Smit ME, McGregor SR, Sun H. et al. A PXY-mediated tran-scriptional network integrates signaling mechanisms to con-trol vascular development in Arabidopsis. Plant Cell. 2020; 32: 319-35

[44] Sivitz AB, Reinders A, Ward JM. Arabidopsis sucrose trans-porter AtSUC1 is important for pollen germination and sucrose-induced anthocyanin accumulation. Plant Physiol. 2008; 147: 92-100

[45] Cao MJ, Wang Z, Wirtz M. et al. SULTR3;1 is a chloroplast-localized sulfate transporter in Arabidopsis thaliana. Plant J. 2013; 73:607-16

[46] Madsen SR, Olsen CE, Nour-Eldin HH. et al. Elucidating the role of transport processes in leaf glucosinolate distribution. Plant Physiol. 2014; 166:1450-62

[47] Nour-Eldin HH, Andersen TG, Burow M. et al. NRT/PTR trans-porters are essential for translocation of glucosinolate defence compounds to seeds. Nature. 2012; 488:531-4

[48] Lee YH, Foster J, Chen J. et al. AAP1 transports uncharged amino acids into roots of Arabidopsis. Plant J. 2007; 50:305-19

[49] Kuang B, Zhao X, Zhou C. et al. Role of UDP-glucuronic acid decarboxylase in xylan biosynthesis in Arabidopsis. Mol Plant. 2016; 9:1119-31

[50] Takenaka Y, Watanabe Y, Schuetz M. et al. Patterned deposition of xylan and lignin is independent from that of the secondary wall cellulose of Arabidopsis xylem vessels. Plant Cell. 2018; 30: 2663-76

[51] Brown DM, Zeef LA, Ellis J. et al. Identification of novel genes in Arabidopsis involved in secondary cell wall formation using expression profiling and reverse genetics. Plant Cell. 2005; 17: 2281-95

[52] BauerS, VasuP, PerssonS. et al. Development and application of a suite of polysaccharide-degrading enzymes for analyzing plant cell walls. Proc Natl Acad Sci USA. 2006; 103:11417-22

[53] Vilaine F, Kerchev P, Clément G. et al. Increased expression of a phloem membrane protein encoded by NHL26 alters phloem export and sugar partitioning in Arabidopsis. Plant Cell. 2013; 25: 1689-708

[54] Truernit E. Sieve elements and their cell neighbours in the Arabidopsis root - roles and relationships. J Plant Physiol. 2022; 268:153569

[55] Chisholm ST, Parra MA, Anderberg RJ. et al. Arabidopsis RTM1 and RTM2 genes function in phloem to restrict long-distance movement of tobacco etch virus. Plant Physiol. 2001; 127:1667-75

[56] Sanchez P, Nehlin L, Greb T. From thin to thick: major transitions during stem development. Trends Plant Sci. 2012; 17:113-21

[57] Roszak P, Heo JO, Blob B. et al. Cell-by-cell dissection of phloem development links a maturation gradient to cell specialization. Science. 2021;374:eaba5531

[58] Gardiner J, Sherr I, Scarpella E. Expression of DOF genes identi-fies early stages of vascular development in Arabidopsis leaves. Int J Dev Biol. 2010; 54:1389-96

[59] Marhava P, Bassukas AEL, Zourelidou M. et al. A molecular rheostat adjusts auxin flux to promote root protophloem dif-ferentiation. Nature. 2018; 558:297-300

[60] Cai X, Chang L, Zhang T. et al. Impacts of allopolyploidization and structural variation on intraspecific diversification in Bras-sica rapa. Genome Biol. 2021; 22:166

[61] Ramirez-Gonzalez RH, Borrill P, Lang D. et al. The transcriptional landscape of polyploid wheat. Science. 2018;361:eaar6089

[62] Burian A, Paszkiewicz G, Nguyen KT. et al. Specification of leaf dorsiventrality via a prepatterned binary readout of a uniform auxin input. Nat Plants. 2022; 8:269-80

[63] Yu C, Yan C, Liu Y. et al. Upregulation of a KN1 homolog by trans-poson insertion promotes leafy head development in lettuce. Proc Natl Acad Sci USA. 2020; 117:33668-78

[64] RenW, WuF, Bai J. et al. BcpLH organizes a specific subset of microRNAs to form a leafy head in Chinese cabbage (Brassica rapa ssp. pekinensis). Hortic Res. 2020; 7:1

[65] An G, Yu C, Yan C. et al. Loss-of-function of SAWTOOTH 1 affects leaf dorsiventrality genes to promote leafy heads in lettuce. Plant Cell. 2022; 34:4329-47

[66] Sun W, Zhang Z, Bonnema G. et al. Gene expression, tran-scription factor binding and histone modification predict leaf adaxial-abaxial polarity related genes. Hortic Plant J. 2024; 10: 971-82

[67] Alemán-Báez J, Acevedo-Zamora JF, Bucher J. et al. Expression changes of miRNA-regulated genes associated with the forma-tion of the leafy head in cabbage. Hortic Plant J. 2024; 10:1007-19

[68] Zhang Y, Liang J, Cai X. et al. Divergence of three BRX homoeologs in Brassica rapa and its effect on leaf morphology. Hortic Res. 2021; 8:68

[69] Rupp HM. Effects of ethylene on secondary xylem formation in Arabidopsis thaliana. Michigan Technological University; 2013

[70] Etchells JP, Provost CM, Turner SR. Plant vascular cell division is maintained by an interaction between PXY and ethylene signalling. PLoS Genet. 2012; 8:e1002997

[71] Miller AR, Pengelly WL, Roberts LW. Introduction of xylem differentiation in lactuca by ethylene. Plant Physiol. 1984; 75: 1165-6

[72] Seyfferth C, Wessels BA, Vahala J. et al. Populus PtERF85 balances xylem cell expansion and secondary cell wall formation in hybrid aspen. Cells. 2021; 10:1971

[73] Furuta KM, Yadav SR, Lehesranta S. et al. Plant development. Arabidopsis NAC45/86 direct sieve element morphogenesis cul-minating in enucleation. Science. 2014; 345:933-7

[74] Kim H, Zhou J, Kumar D. et al. SHORTROOT-mediated intercellu-lar signals coordinate phloem development in Arabidopsis roots. Plant Cell. 2020; 32:1519-35

[75] Hughes TE, Langdale JA, Kelly S. The impact of widespread reg-ulatory neofunctionalization on homeolog gene evolution fol-lowing whole-genome duplication in maize. Genome Res. 2014; 24: 1348-55

[76] Pophaly SD, Tellier A. Population level purifying selection and gene expression shape subgenome evolution in maize. Mol Biol Evol. 2015; 32:3226-35

[77] Salman-Minkov A, Sabath N, Mayrose I. Whole-genome duplica-tion as a key factor in crop domestication. Nature Plants. 2016; 2: 16115

[78] Cheng F, Wu J, Fang L. et al. Biased gene fractionation and dominant gene expression among the subgenomes of Brassica rapa. PLoS One. 2012; 7:e36442

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

[80] Satija R, Farrell JA, Gennert D. et al. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015; 33:495-502

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

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

[83] Bu D, Luo H, Huo P. et al. KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. 2021;49:W317-25

[84] Chen H, Wang T, He X. et al.BRAD V3.0:an upgraded Brassi-caceae database. Nucleic Acids Res. 2021;50:D1432-41

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

[86] Xu HX, Meng D, Yang Q. et al. Sorbitol induces flower bud for-mation via the MADS-box transcription factor EjCAL in loquat. J Integr Plant Biol. 2023; 65:1241-61