Sandra Thibivilliers
Marc Libault- Department of Agronomy and Horticulture, Center for Plant Science Innovation, University of Nebraska-Lincoln, Lincoln, NE, United States
Plants are composed of cells that physically interact and constantly adapt to their environment. To reveal the contribution of each plant cells to the biology of the entire organism, their molecular, morphological, and physiological attributes must be quantified and analyzed in the context of the morphology of the plant organs. The emergence of single-cell/nucleus omics technologies now allows plant biologists to access different modalities of individual cells including their epigenome and transcriptome to reveal the unique molecular properties of each cell composing the plant and their dynamic regulation during cell differentiation and in response to their environment. In this manuscript, we provide a perspective regarding the challenges and strategies to collect plant single-cell biological datasets and their analysis in the context of cellular interactions. As an example, we provide an analysis of the transcriptional regulation of the Arabidopsis genes controlling the differentiation of the root hair cells at the single-cell level. We also discuss the perspective of the use of spatial profiling to complement existing plant single-cell omics.
Introduction: How Single-Cell Approach Can Help to Enhance Our Understanding of Plants as Biological Systems?
Plants are complex and very dynamic biological systems composed of various cell types that communicate together and constantly respond to their environment. Therefore, to better understand plants as biological systems, there is a need to understand the contribution of each cell to the biology of the organism, to reveal the unique and dynamic response of each plant cell to environmental stimuli, and to characterize how cell-to-cell communication plays a role in controlling these responses. To reach such knowledge, plant scientists must gain molecular information from each cell composing the tissue/organ/plant and analyze this information in the context of the spatial organization of the organ and interaction of the cells (Figure 1).
Figure 1. Schematic representation of the analysis of plant cell-to-cell interaction using single-cell omics technologies. Plant organs are composed by physically interacting somatic cells that constantly adapt to their biological environment (A). The analysis of the epigenomic and transcriptomic profiles of isolated plant cells (i.e., epigenome and transcriptome) requires the isolation of individual cells and nuclei (B) before the construction of DNA sequencing libraries (C). Therefore, the relative position of the cell in the organ and its interaction with somatic, microbial, and pathogenic cells is lost. A systems view of plant cellular communication would require the capture of cellular modalities in the context of the morphology of the plant sample before their integration using computational tools (D).
The emergence of single-cell omics technologies and their recent application on plant organs now enable the characterization of the molecular attributes of thousands of cells in one single experiment (Rich-Griffin et al., 2020; Shaw et al., 2021). For instance, recent studies reported the establishment of Arabidopsis roots, stomatal cells, and maize anther cell transcriptomes at the single-cell level using isolated protoplasts or nuclei as input (Denyer et al., 2019; Jean-Baptiste et al., 2019; Nelms and Walbot, 2019; Ryu et al., 2019; Shulse et al., 2019; Zhang et al., 2019; Liu et al., 2020b; Farmer et al., 2021). To complement this first set of information, single-cell ATAC-seq was also applied on nuclei isolated from Arabidopsis roots and various maize organs to reveal the differential chromatin accessibility between plant cell types (Dorrity et al., 2020; Marand et al., 2020; Farmer et al., 2021). Recently, Farmer et al. (2021) integrated single-nuclei RNA-seq and ATAC-seq datasets to reveal the impact of chromatin accessibility in controlling gene expression and the differential organization of the Arabidopsis genome between cell types. While informative, these analyses require the abolition of the interactions between plant cells, a pre-requirement to the study of the molecular attributes of each plant cell. In this manuscript, we provide an update about the current limitations and advantages associated with the use of single-cell omics technologies to study plant cell biology notably in the context of cell-to-cell communication. As an example, we provide a comprehensive analysis of the differential transcriptional regulation between H and N cells, epidermal root cells differentiating into trichoblasts and atrichoblasts, respectively.
Accessing Isolated Plant Biological Entities to Conduct Omics Analyses
Plant biologists face unique challenges when characterizing the molecular attributes of each plant cell. First, the presence of the cell wall prevents the isolation of plant cells. To overcome this first difficulty, several groups used enzymatic cocktails to digest the plant cell wall and release plant protoplasts. These protoplasts were used as input into microfluidic systems to generate single-cell barcoded cDNA libraries and establish the transcriptomes of the Arabidopsis root cells (Denyer et al., 2019; Jean-Baptiste et al., 2019; Ryu et al., 2019; Shulse et al., 2019; Turco et al., 2019; Zhang et al., 2019; Shahan et al., 2020; Song et al., 2020; Wendrich et al., 2020), stomatal cells (Liu et al., 2020b; Lopez-Anido et al., 2020), leaf phloem cells (Kim et al., 2020), and sperm cells (Misra et al., 2019), and the transcriptomes of the cells composing the rice root tip (Liu et al., 2021) and the maize shoot apical meristem (Satterlee et al., 2020), developing ears (Xu et al., 2021) and anther (Nelms and Walbot, 2019). However, the release of a representative and viable population of protoplasts would require constant optimization by taking into consideration the unique biochemical composition of the cell wall without compromising protoplast viability. Therefore, the establishment of a protoplast-based single-cell transcriptome might be restricted to a limited number of plant species, organs, and cell types. For instance, Shulse et al. (2019) noticed that the transcriptome of differentiated endodermal cells that are characterized by the suberization of their cell wall was missing in Arabidopsis root protoplast single-cell transcriptomic datasets (Shulse et al., 2019). It is also important to acknowledge that plant protoplasts are prone to bursting and that protoplastization leads to the activation of the expression of protoplast-induced genes (Birnbaum et al., 2003).
As an alternative to the use of plant protoplasts, plant nuclei were recently utilized to establish biologically meaningful transcriptomic information from the Arabidopsis root (Farmer et al., 2021; Long et al., 2021), inflorescences (Sunaga-Franze et al., 2020), and seeds (Picard et al., 2020), and the tomato shoot apex (Tian et al., 2020). Besides, the use of isolated nuclei revealed the transcriptome of Arabidopsis roots cells that were not captured by protoplast-based single-cell transcriptomes likely due to the limited digestibility of their cell wall (Farmer et al., 2021). Therefore, this method represents an alternative to the use of plant protoplasts to broadly access the transcriptome of plant cells from different species and organs (Sunaga-Franze et al., 2020; Thibivilliers et al., 2020; Thibivilliers and Libault, 2021). However, accessing a more limited pool of polyadenylated transcripts from an isolated nucleus compared to an entire cell would necessarily lead to the detection of a lower number of expressed genes per nucleus vs. per protoplast. Considering that isolated nuclei are also prone to RNA leakage when not properly manipulated, the depth of the nuclear transcriptome might be low when using damaged nuclei. On the other hand, the nuclear transcriptome could be considered as a snapshot of the dynamic transcriptional activity of the genes while the cellular transcriptome may represent an integration of gene activity over time [i.e., the half-life of the cellular mRNA is estimated at 9 h in human cells (Schwanhausser et al., 2011)]. Hence, the nuclear transcriptome should allow the characterization of the early and more subtle responses of the cell in response to stress. Nevertheless, despite its challenges, it has been demonstrated that the nuclear transcriptome is sufficient to decode the tissue heterogeneity to a similar level to the cellular transcriptome.
Multi-Omics Approaches to Reveal Plant Cell Dynamics
Accessing one molecular modality of plant cells, such as their transcriptome, is a major milestone. However, to capture the entire diversity and subtle differences existing between cells and to reveal cell-type-specific regulatory networks and biological processes, there is a need to characterize and integrate different modalities at the single-cell level (Iacono et al., 2019; Hu et al., 2020; Jackson et al., 2020). Previous studies revealed changes in the patterns of histone modifications and gene expression of the Arabidopsis stomatal cells (Lee et al., 2019), and in the profile of methylation of the Arabidopsis root cell types (Kawakatsu et al., 2016). More recently, Dorrity et al. (2020) and Farmer et al. (2021) characterized the profiles of chromatin accessibility of Arabidopsis root cells at the single-cell level using microfluidic technology on isolated plant nuclei (Dorrity et al., 2020; Farmer et al., 2021). A similar study was also conducted on maize axillary buds, inflorescences, whole seedling, embryonic root tips, and post-embryonic crown roots cells (Marand et al., 2020). As a first effort in integrating various molecular markers of plant cells, Horvath et al. (2019) revealed that CG-methylated Arabidopsis genes are constitutively expressed (Horvath et al., 2019). More recently, Farmer et al. (2021) revealed the impact of the profiles of chromatin accessibility in regulating the transcriptional activity of the Arabidopsis root cells (Farmer et al., 2021). Such approaches should be expanded to other modalities, potentially gained at the same time from the same cell, to maximize dataset integration and to highlight the relationships existing between structural and chemical changes on the genomic DNA, somatic mutations, their impact on controlling gene expression and protein abundance. Considering the recent emergence of real multi-omics technology (e.g., analysis of gene expression and profile of chromatin accessibility from the same cell/nucleus using 10x Genomics technology), such technology needs to be expanded to other biological modalities to reveal the diverse and dynamic use of genomic information, proteome and metabolome of thousands of individual plant cells.
Characterize Molecular Modalities at the Single-Cell Level in the Context of Plant Cell-to-Cell Interactions
Proximal interactions and distal communication between cells and organs play critical roles in plant biology. Local communication between two cells depends on the formation of plasmodesmata that connect the cytoplasm of neighboring cells to allow the exchange of proteins, metabolites, and nucleotidic sequences. For instance, the transportation of auxin and cytokinin between root cells via plasmodesmata plays a major role in controlling plant organ differentiation, such as the initiation of the formation of lateral root (Mellor et al., 2020) and legume nodule, plant organ resulting from the symbiotic interaction between legumes and Rhizobia (Fisher et al., 2018). Distal communication also plays a critical role in regulating biological processes. For instance, legume nodulation is controlled by the autoregulation of nodulation mechanism, a distal communication system between the canopy and the root of legume plants that regulates the formation of nodules and, as a consequence, nitrogen fixation efficiency (Kassaw et al., 2015; Wang et al., 2018; Suzaki and Nishida, 2019; de Bruijn, 2020).
Besides their proximal and distal interactions with other somatic cells, plant cells are also subject to interactions with a diverse population of microbes. Therefore, plant cells must constantly adapt their response upon recognition of symbiotic and pathogenic microbes notably by regulating cell-to-cell trafficking (Aung et al., 2020). Decades of work on plant microbes interactions revealed that the response of plant cells to microbial infection is complex and sequential. Upon recognition of the pathogen, a first response, the MAMP-triggered immunity response, is initiated. Later, the effector-triggered immunity response will allow the infected plant cell to initiate a more specific response to the microbe (Jones and Dangl, 2006; Naveed et al., 2020). These two types of immunity form the plant cell-autonomous immunity. Increasing the complexity of these interactions, the neighboring plant cells to an event of infection will trigger a non-cell-autonomous immunity response to minimize new events of infection by the same microbe. This immunity is activated upon communication between plant cells (Yan et al., 2019; Aung et al., 2020; Li et al., 2020; Zeng et al., 2020). Considering that different events of interaction and infection co-occur in a complex organ, the molecular characterization of the cell-autonomous and non-cell-autonomous responses of plant cells to a pathogenic infection remains challenging when conducted at the level of complex tissues and organs. Similar challenges are also faced by plant scientists when considering other complex cellular interactions including mutualistic symbiotic interactions between plants and microorganisms (e.g., legume nodulation and arbuscular mycorrhization) or when considering the interactions of plants with multicellular organisms, such as nematodes and insects. Single-cell approaches represent an attractive solution to reveal the cell-autonomous and non-cell-autonomous regulatory programs activated and repressed by the plant in response to microbial infection and pathogenic organisms. Ultimately, considering that the transcriptome of plant cells at different stages of infection will be captured and considering the development of performant computational tools to create transcriptomic trajectories (Qiu et al., 2017; Hao et al., 2020), the use of single-cell omics technology will clarify the sequential transcriptomic response of the plant cell to pathogenic infections.
To enhance our understanding of the biology of the plants as a complex and organized cellular system, the molecular attributes of each cell should be characterized in the context of the morphology of the tissue/organ. However, as mentioned above, a pre-requirement to the use of single-cell omics technologies is the dissociation of the tissue to access independent biological entities (i.e., cells or nuclei) and, consequently, the loss of the spatial organization of the cells in the tissue. This limitation is partially recovered by the use of performant dimensionality reduction methods (e.g., UMAP and t-SNE) that allow the annotation of the plant cells based on their molecular attribute. Therefore, the dividing/differentiating cells have been reported to be located in the center of dimensionality reduction maps while differentiated cells are located at their periphery (Ryu et al., 2019; Farmer et al., 2021).
The Arabidopsis Root Single-Cell Transcriptome Highlights the Role of Cell-to-Cell Interactions in Controlling Root Hair Differentiation
To evaluate the usefulness of single-cell transcriptomes in the context of cell-to-cell communication, we looked at the transcriptional activity of the Arabidopsis genes involved in the differentiation process and patterning of the root epidermal cells, a biological process that depends on intercellular communication between cortical and epidermal cells. The genes controlling the differentiation of the epidermal cells into H and N cells (i.e., trichoblasts and atrichoblasts, respectively) have been well characterized through a series of functional genomic studies (Salazar-Henao et al., 2016). Upon detection of a signal generated in large quantities by the two cortical cells underlying an H cell, the root epidermal leucine-rich repeat receptor SCRAMBLED (SCM) repressed the expression of the MYB transcription factor WEREWOLF exclusively in the H cells (Kwak et al., 2005; Kwak and Schiefelbein, 2007; Wang et al., 2019). Mining the recently published Arabidopsis single-cell/nucleus RNA-seq UMAP projections (Farmer et al., 2021), we found WER mostly transcriptionally active only in a subpopulation of atrichoblasts, and in the lower branch of the cortical cells (Figure 2). As expected, WER was not significantly expressed in the H cells. In the N cells, WER interacts with GLABRA3 (GL3), ENHANCER OF GLABRA3 (EGL3), and TRANSPARENT TESTA GLABRA (TTG) to induce the expression of CAPRICE (CPC), TRIPTYCHON (TRY), and GLABRA2 (GL2) (Koshino-Kimura et al., 2005). Mining the single-cell/nucleus transcriptome, we observed the co-expression of TTG, CPC, TRY, and GL2 genes, and, at a lower level, the expression of GL3 and EGL3 in the atrichoblast cluster (Figure 2). Besides, GL3, EGL3, TTG, and CPC are also expressed in the trichoblast cluster. This observation is supported by the role of these genes in controlling the differentiation of H cells notably by repressing the expression of GL2 (Kurata et al., 2005). Indeed, we did not detect any GL2 transcripts, nor TRY, into the cells and nuclei composing the trichoblast cluster as supported by previously published works (Rerie et al., 1994; Di Cristina et al., 1996; Masucci et al., 1996). This transcriptomic analysis at the single-cell level supports functional genomic studies showing the co-expression of major regulatory genes controlling the differentiation process and patterning of the root epidermal cells. However, the isolations of protoplasts or nuclei before conducting single-cell omics analyses necessarily lead to the loss of the physical interactions between cells. Therefore, the molecular information collected on isolated cells cannot be analyzed in the context of the relative position of the cells in the organ and their interactions with their neighboring cells.
Figure 2. Transcriptional activity of the Arabidopsis genes playing a major role during root epidermal cell differentiation, a biological process that depends on intercellular communication (top left panel). The relative expression levels of the genes controlling Arabidopsis epidermal cell differentiation are highlighted in yellow/red color. The cortical, N, and H cell clusters are, respectively, highlighted in blue, green, and red on the top right panel. Shortly, a biological signal produced in high quantities by two cortical cells underlying a single H root epidermal cell is detected by the epidermal cell-localized leucine-rich repeat receptor SCM. This recognition notably leads to the repression of the expression of the WER gene and the initiation of epidermal cell differentiation. Based on the previous studies, WER and GL2 are expected to be specifically expressed in the N cells (green cluster). Besides, we found that WER is also expressed in cortical cells (blue cells). As expected, the expression of GL3 and EGL3 is almost exclusively restricted to the H cells (red cells). Sc/sNucRNA-seq datasets also confirmed the transcriptional activity of TTG1 in both N and H cells (green and red cells). However, CPC seems to have a ubiquitous transcriptional activity.
Perspectives
To gain a systems view of plant cellular communication, there is a need to quantify molecular modalities of individual cells in the context of the morphology of the organ analyzed, or at least to bridge information gained from single-cell approaches with the relative position of the cells.
Spatially resolved transcriptomics technologies offer an opportunity to access the transcriptome of cells or groups of cells in the context of the morphology of the organ analyzed. While spatial methods are now routinely applied in animal science, the plant scientific community starts to discuss their use to better understand gene regulation (Giacomello and Lundeberg, 2018; Gurazada et al., 2021). Spatial transcriptomics technologies will also facilitate the characterization of cell-type-specific molecular markers. While marker genes are well characterized in model plant species, such as Arabidopsis thaliana (e.g., Denyer et al., 2019; Jean-Baptiste et al., 2019; Ryu et al., 2019; Shulse et al., 2019; Zhang et al., 2019; Farmer et al., 2021), non-model species suffer from limited access to the cell-type-specific marker genes needed to properly annotate the various cell types composing an organ based in their transcriptomic information. The use of spatial transcriptomics technology on plant organ cross-sections will allow the characterization of cell-type-specific marker genes in the morphological context of the organ analyzed. This knowledge can be used to enhance the functional annotation of plant single-cell clusters; especially from plant species or organs where the number of functionally validated single-cell marker genes is limited.
Two strategies have been used to enable spatial transcriptomics analyses. On the one hand, Slide-seq (Rodriques et al., 2019) and Visium technology from 10x Genomics® are based on the use of nucleotide spatial barcodes arrayed on a slide. On the other hand, High-Definition Spatial Transcriptomics (Vickovic et al., 2019) and Spatial Molecular Imaging technology from Nanostring® offer very high-level resolution transcriptomes of complex organs. Applied to plant samples, these technologies will enable the accurate analysis of the differential use of the genomic information between plant cells and the impact of cell-to-cell interactions in controlling biological processes. However, it is important to acknowledge that spatial omics technologies suffer from several limitations. First, the resolution of the information gained might requires the use of computational methods for deconvolution of the Visium pots (e.g., the 10x Genomics Visium Gene Expression spots are 55 μm in diameter leading to the analysis of the transcriptome of several plant cells per spots; Andersson et al., 2020; Elosua-Bayes et al., 2021). Second, spatial omics technologies are currently almost exclusively restricted to the analysis of the transcriptome. To overcome this limitation, Liu et al. (2020a) recently developed Deterministic Barcoding in Tissue (DBiT-seq) technology allowing the quantification of transcripts abundance and the detection of proteins of interest in the context of the morphology of tissue (Liu et al., 2020a). To date, this method has been applied on mouse embryos and will likely require substantial optimization before implementing its use on plant samples. Such an approach will need to be expanded to cover additional molecular modalities to gain a deeper understanding of plant cell biology and to reveal the impact of these modalities on cell biology, physiology, and morphology.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/geo/, GSE155304.
Author Contributions
ST drafted the manuscript. ML edited the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the US National Science Foundation (IOS# 1339194 and 1854326), the Center for Plant Science Innovation, and the Department of Agronomy and Horticulture at the University of Nebraska-Lincoln.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Andersson, A., Bergenstråhle, J., Asp, M., Bergenstråhle, L., Jurek, A., Fernández Navarro, J., et al. (2020). Single-cell and spatial transcriptomics enables probabilistic inference of cell type topography. Commun. Biol. 3, 1–8. doi: 10.1038/s42003-020-01247-y
Aung, K., Kim, P., Li, Z., Joe, A., Kvitko, B., Alfano, J. R., et al. (2020). pathogenic bacteria target plant plasmodesmata to colonize and invade surrounding tissues. Plant Cell 32, 595–611. doi: 10.1105/tpc.19.00707
Birnbaum, K., Shasha, D. E., Wang, J. Y., Jung, J. W., Lambert, G. M., Galbraith, D. W., et al. (2003). A gene expression map of the Arabidopsis root. Science 302, 1956–1960. doi: 10.1126/science.1090022
de Bruijn, F. J. (ed.) (2020). “Autoregulation of nodule numbers (AON) in Medicago truncatula: introduction,” in The Model Legume Medicago truncatula (Wiley-Blackwell), 809–810.
Denyer, T., Ma, X., Klesen, S., Scacchi, E., Nieselt, K., and Timmermans, M. C. P. (2019). Spatiotemporal developmental trajectories in the Arabidopsis root revealed using high-throughput single-cell RNA sequencing. Dev. Cell 48, 840–852. doi: 10.1016/j.devcel.2019.02.022
Di Cristina, M., Sessa, G., Dolan, L., Linstead, P., Baima, S., Ruberti, I., et al. (1996). The Arabidopsis Athb-10 (GLABRA2) is an HD-Zip protein required for regulation of root hair development. Plant J. 10, 393–402. doi: 10.1046/j.1365-313X.1996.10030393.x
Dorrity, M. W., Alexandre, C., Hamm, M., Vigil, A.-L., Fields, S., Queitsch, C., et al. (2020). The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution. bioRxiv [Preprint]. doi: 10.1101/2020.07.17.204792
Elosua-Bayes, M., Nieto, P., Mereu, E., Gut, I., and Heyn, H. (2021). SPOTlight: seeded NMF regression to deconvolute spatial transcriptomics spots with single-cell transcriptomes. Nucleic Acids Res. 49:e50. doi: 10.1093/nar/gkab043
Farmer, A., Thibivilliers, S., Ryu, K. H., Schiefelbein, J., and Libault, M. (2021). 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 14, 372–383. doi: 10.1016/j.molp.2021.01.001
Fisher, J., Gaillard, P., Fellbaum, C. R., Subramanian, S., and Smith, S. (2018). Quantitative 3D imaging of cell level auxin and cytokinin response ratios in soybean roots and nodules. Plant Cell Environ. 41, 2080–2092. doi: 10.1111/pce.13169
Giacomello, S., and Lundeberg, J. (2018). Preparation of plant tissue to enable spatial transcriptomics profiling using barcoded microarrays. Nat. Protoc. 13, 2425–2446. doi: 10.1038/s41596-018-0046-1
Gurazada, S. G. R., Cox, K. L., Czymmek, K. J., and Meyers, B. C. (2021). Space: the final frontier — achieving single-cell, spatially resolved transcriptomics in plants. Emerg. Top. Life Sci. 5, 179–188. doi: 10.1042/ETLS20200274
Hao, Y., Hao, S., Andersen-Nissen, E., Mauck, W. M., Zheng, S., Butler, A., et al. (2020). Integrated analysis of multimodal single-cell data. bioRxiv [Preprint]. doi: 10.1101/2020.08.27.271130
Horvath, R., Laenen, B., Takuno, S., and Slotte, T. (2019). Single-cell expression noise and gene-body methylation in Arabidopsis thaliana. Heredity 123, 81–91. doi: 10.1038/s41437-018-0181-z
Hu, X., Hu, Y., Wu, F., Leung, R. W. T., and Qin, J. (2020). Integration of single-cell multi-omics for gene regulatory network inference. Comput. Struct. Biotechnol. J. 18, 1925–1938. doi: 10.1016/j.csbj.2020.06.033
Iacono, G., Massoni-Badosa, R., and Heyn, H. (2019). Single-cell transcriptomics unveils gene regulatory network plasticity. Genome Biol. 20, 1–20. doi: 10.1186/s13059-019-1713-4
Jackson, C. A., Castro, D. M., Saldi, G.-A., Bonneau, R., and Gresham, D. (2020). Gene regulatory network reconstruction using single-cell RNA sequencing of barcoded genotypes in diverse environments. elife 9:e51254. doi: 10.7554/eLife.51254
Jean-Baptiste, K., McFaline-Figueroa, J. L., Alexandre, C. M., Dorrity, M. W., Saunders, L., Bubb, K. L., et al. (2019). Dynamics of gene expression in single root cells of A. thaliana. Plant Cell 31, 993–1011. doi: 10.1105/tpc.18.00785
Jones, J. D. G., and Dangl, J. L. (2006). The plant immune system. Nature 444, 323–329. doi: 10.1038/nature05286
Kassaw, T., Bridges, W. Jr., and Frugoli, J. (2015). Multiple autoregulation of nodulation (AON) signals identified through split root analysis of Medicago truncatula sunn and rdn1 Mutants. Plants 4, 209–224. doi: 10.3390/plants4020209
Kawakatsu, T., Stuart, T., Valdes, M., Breakfield, N., Schmitz, R. J., Nery, J. R., et al. (2016). Unique cell-type-specific patterns of DNA methylation in the root meristem. Nat. Plants 2:16058. doi: 10.1038/nplants.2016.58
Kim, J.-Y., Symeonidi, E., Pang, T. Y., Denyer, T., Weidauer, D., Bezrutczyk, M., et al. (2020). Unique and distinct identities and functions of leaf phloem cells revealed by single cell transcriptomics. bioRxiv [Preprint]. doi: 10.1101/2020.07.24.205583
Koshino-Kimura, Y., Wada, T., Tachibana, T., Tsugeki, R., Ishiguro, S., and Okada, K. (2005). Regulation of CAPRICE transcription by MYB proteins for root epidermis differentiation in Arabidopsis. Plant Cell Physiol. 46, 817–826. doi: 10.1093/pcp/pci096
Kurata, T., Ishida, T., Kawabata-Awai, C., Noguchi, M., Hattori, S., Sano, R., et al. (2005). Cell-to-cell movement of the CAPRICE protein in Arabidopsis root epidermal cell differentiation. Development 132, 5387–5398. doi: 10.1242/dev.02139
Kwak, S. H., and Schiefelbein, J. (2007). The role of the SCRAMBLED receptor-like kinase in patterning the Arabidopsis root epidermis. Dev. Biol. 302, 118–131. doi: 10.1016/j.ydbio.2006.09.009
Kwak, S. H., Shen, R., and Schiefelbein, J. (2005). Positional signaling mediated by a receptor-like kinase in Arabidopsis. Science 307, 1111–1113. doi: 10.1126/science.1105373
Lee, L. R., Wengier, D. L., and Bergmann, D. C. (2019). Cell-type–specific transcriptome and histone modification dynamics during cellular reprogramming in the Arabidopsis stomatal lineage. Proc. Natl. Acad. Sci. 116, 21914–21924. doi: 10.1073/pnas.1911400116
Li, Z., Variz, H., Chen, Y., Liu, S.-L., and Aung, K. (2020). Plasmodesmata-dependent intercellular movement of bacterial effectors. bioRxiv 12:464 [Preprint]. doi: 10.1101/2020.12.24.424262
Liu, Q., Liang, Z., Feng, D., Jiang, S., Wang, Y., Du, Z., et al. (2021). Transcriptional landscape of rice roots at the single-cell resolution. Mol. Plant 14, 384–394. doi: 10.1016/j.molp.2020.12.014
Liu, Y., Yang, M., Deng, Y., Su, G., Enninful, A., Guo, C. C., et al. (2020a). High-spatial-resolution multi-omics sequencing via deterministic barcoding in tissue. Cell 183, 1665–1681. doi: 10.1016/j.cell.2020.10.026
Liu, Z., Zhou, Y., Guo, J., Li, J., Tian, Z., Zhu, Z., et al. (2020b). Global dynamic molecular profiling of stomatal lineage cell development by single-cell RNA sequencing. Mol. Plant 13, 1178–1193. doi: 10.1016/j.molp.2020.06.010
Long, Y., Liu, Z., Jia, J., Mo, W., Fang, L., Lu, D., et al. (2021). FlsnRNA-seq: protoplasting-free full-length single-nucleus RNA profiling in plants. Genome Biol. 22, 1–14. doi: 10.1186/s13059-021-02288-0
Lopez-Anido, C. B., Vatén, A., Smoot, N. K., Sharma, N., Guo, V., Gong, Y., et al. (2020). Single-cell resolution of lineage trajectories in the Arabidopsis stomatal lineage and developing leaf. bioRxiv 56 [Preprint]. doi: 10.1101/2020.09.08.288498
Marand, A. P., Chen, Z., Gallavotti, A., and Schmitz, R. J. (2020). A cis−regulatory atlas in maize at single-cell resolution. bioRxiv [Preprint]. doi: 10.1101/2020.09.27.315499
Masucci, J. D., Rerie, W. G., Foreman, D. R., Zhang, M., Galway, M. E., Marks, M. D., et al. (1996). The homeobox gene GLABRA2 is required for position-dependent cell differentiation in the root epidermis of Arabidopsis thaliana. Development 122, 1253–1260. doi: 10.1242/dev.122.4.1253
Mellor, N. L., Voß, U., Janes, G., Bennett, M. J., Wells, D. M., and Band, L. R. (2020). Auxin fluxes through plasmodesmata modify root-tip auxin distribution. Development 147:dev181669. doi: 10.1242/dev.181669
Misra, C. S., Santos, M. R., Rafael-Fernandes, M., Martins, N. P., Monteiro, M., and Becker, J. D. (2019). Transcriptomics of Arabidopsis sperm cells at single-cell resolution. Plant Reprod. 32, 29–38. doi: 10.1007/s00497-018-00355-4
Naveed, Z. A., Wei, X., Chen, J., Mubeen, H., and Ali, G. S. (2020). The PTI to ETI continuum in phytophthora-plant interactions. Front. Plant Sci. 11:230. doi: 10.3389/fpls.2020.593905
Nelms, B., and Walbot, V. (2019). Defining the developmental program leading to meiosis in maize. Science 364, 52–56. doi: 10.1126/science.aav6428
Picard, C. L., Povilus, R. A., Williams, B. P., and Gehring, M. (2020). Single nucleus analysis of arabidopsis seeds reveals new cell types and imprinting dynamics. bioRxiv [Preprint]. doi: 10.1101/2020.08.25.267476
Qiu, X., Mao, Q., Tang, Y., Wang, L., Chawla, R., Pliner, H. A., et al. (2017). Reversed graph embedding resolves complex single-cell trajectories. Nat. Methods 14, 979–982. doi: 10.1038/nmeth.4402
Rerie, W. G., Feldmann, K. A., and Marks, M. D. (1994). The GLABRA2 gene encodes a homeo domain protein required for normal trichome development in arabidopsis. Genes Dev. 8, 1388–1399. doi: 10.1101/gad.8.12.1388
Rich-Griffin, C., Stechemesser, A., Finch, J., Lucas, E., Ott, S., and Schafer, P. (2020). Single-cell transcriptomics: a high-resolution avenue for plant functional genomics. Trends Plant Sci. 25, 186–197. doi: 10.1016/j.tplants.2019.10.008
Rodriques, S. G., Stickels, R. R., Goeva, A., Martin, C. A., Murray, E., Vanderburg, C. R., et al. (2019). Slide-seq: a scalable technology for measuring genome-wide expression at high spatial resolution. Science 363, 1463–1467. doi: 10.1126/science.aaw1219
Ryu, K. H., Huang, L., Kang, H. M., and Schiefelbein, J. (2019). Single-cell RNA sequencing resolves molecular relationships among individual plant cells. Plant Physiol. 179, 1444–1456. doi: 10.1104/pp.18.01482
Salazar-Henao, J. E., Vélez-Bermúdez, I. C., and Schmidt, W. (2016). The regulation and plasticity of root hair patterning and morphogenesis. Development 143, 1848–1858. doi: 10.1242/dev.132845
Satterlee, J. W., Strable, J., and Scanlon, M. J. (2020). Plant stem cell organization and differentiation at single-cell resolution. bioRxiv 117 [Preprint], 33689–33699. doi: 10.1101/2020.08.25.267427
Schwanhausser, B., Busse, D., Li, N., Dittmar, G., Schuchhardt, J., Wolf, J., et al. (2011). Global quantification of mammalian gene expression control. Nature 473, 337–342. doi: 10.1038/nature10098
Shahan, R., Hsu, C.-W., Nolan, T. M., Cole, B. J., Taylor, I. W., Vlot, A. H. C., et al. (2020). A single cell Arabidopsis root atlas reveals developmental trajectories in wild type and cell identity mutants. bioRxiv [Preprint]. doi: 10.1101/2020.06.29.178863
Shaw, R., Tian, X., and Xu, J. (2021). Single-cell transcriptome analysis in plants: advances and challenges. Mol. Plant 14, 115–126. doi: 10.1016/j.molp.2020.10.012
Shulse, C. N., Cole, B. J., Ciobanu, D., Lin, J., Yoshinaga, Y., Gouran, M., et al. (2019). High-throughput single-cell transcriptome profiling of plant cell types. Cell Rep. 27, 2241–2247. doi: 10.1016/j.celrep.2019.04.054
Song, Q., Ando, A., Jiang, N., Ikeda, Y., and Chen, Z. J. (2020). Single-cell RNA-seq analysis reveals ploidy-dependent and cell-specific transcriptome changes in arabidopsis female gametophytes. Genome Biol. 21, 1–8. doi: 10.1186/s13059-020-02094-0
Sunaga-Franze, D. Y., Muino, J. M., Braeuning, C., Xu, X., Zong, M., Smaczniak, C., et al. (2020). Single-nuclei RNA-sequencing of plants. bioRxiv [Preprint]. doi: 10.1101/2020.11.14.382812
Suzaki, T., and Nishida, H. (2019). Autoregulation of legume nodulation by sophisticated transcriptional regulatory networks. Mol. Plant 12, 1179–1181. doi: 10.1016/j.molp.2019.07.008
Thibivilliers, S., Anderson, D., and Libault, M. (2020). Isolation of plant root nuclei for single cell RNA sequencing. Curr. Protoc. Plant. Biol. 5:e20120. doi: 10.1002/cppb.20120
Thibivilliers, S., and Libault, M. (2021). Plant single-cell multiomics: cracking the molecular profiles of plant cells. Trends Plant. Sci. 26, 662–663. doi: 10.1016/j.tplants.2021.03.001
Tian, C., Du, Q., Xu, M., Du, F., and Jiao, Y. (2020). Single-nucleus RNA-seq resolves spatiotemporal developmental trajectories in the tomato shoot apex. bioRxiv [Preprint]. doi: 10.1101/2020.10.04.325423
Turco, G. M., Rodriguez-Medina, J., Siebert, S., Han, D., Valderrama-Gómez, M. Á., Vahldick, H., et al. (2019). Molecular mechanisms driving switch behavior in xylem cell differentiation. Cell Rep. 28, 342–351. doi: 10.1016/j.celrep.2019.06.041
Vickovic, S., Eraslan, G., Salmén, F., Klughammer, J., Stenbeck, L., Schapiro, D., et al. (2019). High-definition spatial transcriptomics for in situ tissue profiling. Nat. Methods 16, 987–990. doi: 10.1038/s41592-019-0548-y
Wang, C., Reid, J. B., and Foo, E. (2018). The art of self-control – autoregulation of plant–microbe symbioses. Front. Plant Sci. 9:988. doi: 10.3389/fpls.2018.00988
Wang, W., Ryu, K. H., Barron, C., and Schiefelbein, J. (2019). Root epidermal cell patterning is modulated by a critical residue in the WEREWOLF transcription factor. Plant Physiol. 181, 1239–1256. doi: 10.1104/pp.19.00458
Wendrich, J. R., Yang, B., Vandamme, N., Verstaen, K., Smet, W., Van de Velde, C., et al. (2020). Vascular transcription factors guide plant epidermal responses to limiting phosphate conditions. Science 370:eaay4970. doi: 10.1126/science.aay4970
Xu, X., Crow, M., Rice, B. R., Li, F., Harris, B., Liu, L., et al. (2021). Single-cell RNA sequencing of developing maize ears facilitates functional analysis and trait candidate gene discovery. Dev. Cell 56, 557–568. doi: 10.1016/j.devcel.2020.12.015
Yan, J., Yu, H., Li, B., Fan, A., Melkonian, J., Wang, X., et al. (2019). Cell autonomous and non-autonomous functions of plant intracellular immune receptors in stomatal defense and apoplastic defense. PLoS Pathog. 15:e1008094. doi: 10.1371/journal.ppat.1008094
Zeng, S. M., Lo, E. K. W., Hazelton, B. J., Morales, M. F., and Torii, K. U. (2020). Effective range of non-cell autonomous activator and inhibitor peptides specifying plant stomatal patterning. Development 147:dev192237. doi: 10.1242/dev.192237
Keywords: transcriptomics, single-cell omics, multi-omics analyses, spatial transcriptomics, cell-to-cell interactions
Citation: Thibivilliers S and Libault M (2021) Enhancing Our Understanding of Plant Cell-to-Cell Interactions Using Single-Cell Omics. Front. Plant Sci. 12:696811. doi: 10.3389/fpls.2021.696811
Edited by:
Dario Di Silvestre, Italian National Research Council, ItalyReviewed by:
Chuang Ma, Northwest A & F University, ChinaKoh Aoki, Osaka Prefecture University, Japan
Copyright © 2021 Thibivilliers and Libault. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marc Libault, Marc.libault@unl.edu