Development of a single-cell atlas for woodland strawberry (Fragaria vesca) leaves during early Botrytis cinerea infection using single-cell RNA-seq

Abstract

Pathogen invasion leads to fast, local-to-systemic signal transduction that initiates plant defense responses. Despite tremendous progress in past decades, aspects of this process remain unknown, such as which cell types respond first and how signals are transferred among cell types. Here, we used single-cell RNA-seq of >50 000 single cells to document the gene expression landscape in leaves of woodland strawberry during infection by Botrytis cinerea and identify major cell types. We constructed a single-cell atlas and characterized the distinct gene expression patterns of hydathode, epidermal, and mesophyll cells during the incubation period of B. cinerea infection. Pseudotime trajectory analysis revealed signals of the transition from normal functioning to defense response in epidermal and mesophyll cells upon B. cinerea infection. Genes related to disease resistance showed different expression patterns among cell types: disease resistance-related genes and genes encoding transcription factors were highly expressed in individual cell types and interacted to trigger plant systemic immunity to B. cinerea. This is the first report to document the single-cell transcriptional landscape of the plant pathogenic invasion process; it provides new insights into the holistic dynamics of host–pathogen interactions and can guide the identification of genes and the formulation of strategies for resistant cultivar development.

References

[1] Jones JDG, Dangl JL . The plant immune system. Nature. 2006; 444: 323-9.

[2] Nürnberger T, Brunner F, Kemmerling B et al. Innate immunity in plants and animals: striking similarities and obvious differences. Immunol Rev. 2010; 198: 249-66.

[3] Williamson B, Tudzynski B, Tudzynski P et al. Botrytis cinerea: the cause of grey mould disease. Mol Plant Pathol. 2007; 8: 561-80.

[4] Fillinger S, Elad YB . Botrytis - The Fungus, the Pathogen and its Management in Agricultural Systems. Switzerland: Springer International Publishing, 2016.

[5] Zhang H, Wang L, Dong Y et al. Postharvest biological control of gray mold decay of strawberry with Rhodotorula glutinis. Biol Control. 2007; 40: 287-92.

[6] Francisco AR, Rosario BP, Juan M-B et al. The strawberry plant defense mechanism: a molecular review. Plant Cell Physiol. 2011; 52: 1873-903.

[7] Xiong JS, Zhu HY, Bai YB et al. RNA sequencing-based transcriptome analysis of mature strawberry fruit infected by necrotrophic fungal pathogen Botrytis cinerea. Physiol Mol Plant Pathol. 2018; 104: 77-85.

[8] Hu Z, Chang X, Dai T et al. Metabolic profiling to identify the latent infection of strawberry by Botrytis cinerea. Evol Bioinform Online. 2019; 15: 1176934319838518.

[9] Xu Y, Tong Z, Zhang X et al. Unveiling the mechanisms for the plant volatile organic compound linalool to control gray mold on strawberry fruits. J Agric Food Chem. 2019; 67: 9265-76.

[10] Zhang P, Yang M, Zhang Y et al. Dissecting the single-cell transcriptome network underlying gastric premalignant lesions and early gastric cancer. Cell Rep. 2019; 27: 1934-47.

[11] Shaw R, Tian X, Xu J . Single-cell transcriptome analysis in plants: advances and challenges. Mol Plant. 2020; 14: 115-26.

[12] Seyffert HC, Renema J, Wendrich JR et al. Advances and opportunities of single-cell transcriptomics for plant research. Annu Rev Plant Biol 2021; 72: 847-66.

[13] Tom, 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-52.

[14] 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.

[15] 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-68.

[16] Jean-Baptiste K, McFaline-Figueroa JL, Alexandre CM et al. Dynamics of gene expression in single root cells of Arabidopsis thaliana. Plant Cell 2019; 31: 993-1011.

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

[18] 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.

[19] Lieckfeldt E, Simon-Rosin U, Kose F et al. Gene expression profiling of single epidermal, basal and trichome cells of Arabidopsis thaliana. J Plant Physiol 2008; 165: 1530-44.

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

[21] Zhao C, Johnson BJ, Kositsup B et al. Exploiting secondary growth in Arabidopsis. Construction of xylem and bark cDNA libraries and cloning of three xylem endopeptidases. Plant Physiol 2000; 123: 1185-96.

[22] Yagi H, Nagano AJ, Kim J et al. Fluorescent protein-based imaging and tissue-specific RNA-seq analysis of Arabidopsis hydathodes. J Exp Bot 2020; 72: 1260-70.

[23] Kim H, Choi D, Suh MC . Cuticle ultrastructure, cuticular lipid composition, and gene expression in hypoxia-stressed Arabidopsis stems and leaves. Plant Cell Rep. 2017; 36: 815-27.

[24] Pollard M, Beisson F, Li Y et al. Building lipid barriers: biosynthesis of cutin and suberin. Trends Plant Sci. 2008; 13: 236-46.

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

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

[27] Huang D, Lin W, Deng B et al. Dual-located WHIRLY1 interacting with LHCA1 alters photochemical activities of photosystem I and is involved in light adaptation in Arabidopsis. Int J Mol Sci. 2017; 18: 2352.

[28] Kondo Y, Fujita T, Sugiyama M et al. A novel system for xylem cell differentiation in Arabidopsis thaliana. Mol Plant. 2015; 8: 612-21.

[29] Susek RE, Ausubel FM, Chory J . Signal transduction mutants of Arabidopsis uncouple nuclear CAB and RBCS gene expression from chloroplast development. Cell. 1993; 74: 787-99.

[30] Yang C, Ye Z . Trichomes as models for studying plant cell differentiation. Cell Mol Life Sci. 2013; 70: 1937-48.

[31] 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.

[32] Kubo M, Udagawa M, Nishikubo N et al. Transcription switches for protoxylem and metaxylem vessel formation. Genes Dev 2005; 19: 1855-60.

[33] Petrasch S, Knapp SJ, Kan JALV et al. Grey mould of strawberry, a devastating disease caused by the ubiquitous necrotrophic fungal pathogen Botrytis cinerea. Mol Plant Pathol 2019; 20: 877-92.

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

[35] Guimil S, Dunand C . Patterning of Arabidopsis epidermal cells: epigenetic factors regulate the complex epidermal cell fate pathway. Trends Plant Sci 2006; 11: 601-9.

[36] Gandotra N, Coughlan SJ, Nelson T . The Arabidopsis leaf provascular cell transcriptome is enriched in genes with roles in vein patterning. Plant J 2013; 74: 48-58.

[37] Vanneste S, Coppens F, Lee E et al. Developmental regulation of CYCA2s contributes to tissue-specific proliferation in Arabidopsis. EMBO J 2011; 30: 3430-41.

[38] Wang Y, Huan Q, Chu X et al. Single-cell transcriptome analyses recapitulate the cellular and developmental responses to abiotic stresses in rice. BioRxiv 2020; http://doi.org/10.1101/2020.01.30.926329.

[39] Carmit Z, Zhao Z, Gao YG et al. Multifunctional roles of plant cuticle during plant-pathogen interactions. Front Plant Sci 2018; 9: 1088.

[40] Yuan M, Jiang Z, Bi G et al. Pattern-recognition receptors are required for NLR-mediated plant immunity. Nature 2021; 592: 105-9.

[41] Chisholm ST, Coaker G, Day B et al. Host-microbe interactions: shaping the evolution of the plant immune response. Cell. 2006; 124: 803-14.

[42] Vadassery J, Reichelt M, Hause B et al. CML42-mediated calcium signaling coordinates responses to Spodoptera herbivory and abiotic stresses in Arabidopsis. Plant Physiol. 2012; 159: 1159-75.

[43] Chen L, Zhang L, Xiang S et al. The transcription factor WRKY75 positively regulates jasmonate-mediated plant defense to necrotrophic fungal pathogens. J Exp Bot. 2020; 72: 1473-89.

[44] Dorrity MW, Alexandre CM, Hamm MO et al. The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution. Nat Commun. 2021; 12: 3334.

[45] Liu X-M, An J, Han HJ et al. ZAT11, a zinc finger transcription factor, is a negative regulator of nickel ion tolerance in Arabidopsis. Plant Cell Rep 2014; 33: 2015-21.

[46] Smit ME, McGregor SR, Sun H et al. A PXY-mediated transcriptional network integrates signaling mechanisms to control vascular development in Arabidopsis. Plant Cell 2019; 32: 319-35.

[47] Luo C, Fernie AR, Yan J . Single-cell genomics and epigenomics: technologies and applications in plants. Trends Plant Sci 2020; 25: 1030-40.

[48] Zhu P, Li Q, Zhang C et al. Bcpks12 gene inactivation substantiates biological functions of sclerotium melanization in Botrytis cinerea. Physiol Mol Plant Pathol 2017; 98: 80-4.

[49] Edger PP, VanBuren R, Colle M et al. Single-molecule sequencing and optical mapping yields an improved genome of woodland strawberry (Fragaria vesca) with chromosome-scale contiguity. GigaScience. 2018; 7: 1-7.

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

[51] Anders S, Huber W . Differential expression analysis for sequence count data. Genome Biol 2010; 11: R106.

[52] Conway JR, Lex A, Gehlenborg N . UpSetR: an R package for the visualization of intersecting sets and their properties. Bioinformatics. 2017; 33: 2938-40.

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

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

[55] 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.