Skip to main content

REVIEW article

Front. Plant Sci., 11 August 2023
Sec. Technical Advances in Plant Science
This article is part of the Research Topic Omics and its Integration: A Systems Biology Approach to Understanding Plant Physiology View all 10 articles

Opportunities and challenges in the application of single-cell and spatial transcriptomics in plants

  • 1Ministry of Education Key Laboratory of Environment Remediation and Ecological Health, College of Environmental and Resource Sciences, Zhejiang University, Hangzhou, China
  • 2Key Laboratory of Agricultural Resource and Environment of Zhejiang Province, College of Environmental and Resource Sciences, Zhejiang University, Hangzhou, China

Single-cell and spatial transcriptomics have diverted researchers’ attention from the multicellular level to the single-cell level and spatial information. Single-cell transcriptomes provide insights into the transcriptome at the single-cell level, whereas spatial transcriptomes help preserve spatial information. Although these two omics technologies are helpful and mature, further research is needed to ensure their widespread applicability in plant studies. Reviewing recent research on plant single-cell or spatial transcriptomics, we compared the different experimental methods used in various plants. The limitations and challenges are clear for both single-cell and spatial transcriptomic analyses, such as the lack of applicability, spatial information, or high resolution. Subsequently, we put forth further applications, such as cross-species analysis of roots at the single-cell level and the idea that single-cell transcriptome analysis needs to be combined with other omics analyses to achieve superiority over individual omics analyses. Overall, the results of this review suggest that combining single-cell transcriptomics, spatial transcriptomics, and spatial element distribution can provide a promising research direction, particularly for plant research.

1 Introduction

A transcriptome is the collection of all transcripts in a cell under specific physiological conditions. Transcriptomics can reveal differences in gene expression under different conditions (Shojaee et al., 2021). From the beginning of transcriptome analysis of various organs to the combination of laser capture microdissection (LCM) and high-throughput sequencing for single-cell transcriptome analysis, the demand for the resolution of transcriptome analysis is increasing. Traditional simple transcriptome analysis cannot meet this demand; however, single-cell transcriptome analysis can. As its name implies, a single-cell transcriptome is the transcriptome of a single cell. Since scientists first reported single-cell transcriptome technology, it has undergone considerable development (Tang et al., 2009). In 2015, two research groups from Harvard developed the Drop-seq and InDrop technologies and applied them to study mouse cells (Klein et al., 2015; Macosko et al., 2015; Klein and Macosko, 2017). In 2016, the 10× Genomics Chromium system, which is a droplet-based technique, was developed. This single-cell sequencing system allows high throughput and can detect rare cell types (See et al., 2018; Svensson et al., 2018; Wang et al., 2021b). Since then, the rapid and low-cost gene expression analyses of many single cells have become a reality. Related single-cell sequencing technologies, such as GemCode single-cell technology (Zheng et al., 2017), MARS-seq2.0 (Keren-Shaul et al., 2019), and Paired-seq (Zhang et al., 2020), are constantly improving in many aspects, such as high cell capture capacity, high gene detection capability, low technical cell-to-cell contamination rate, and low cost, making it possible to use single-cell sequencing technology more conveniently and accurately. Currently, there are three main single-cell sequencing approaches: plate-based, combinatorial indexing-based, and bead-based. Plate-based approaches, e.g., Smart-seq2 and CEL-Seq2, sort a single cell into a well of a multi-well plate; however, these approaches have low throughput and are time consuming. Bead-based methods, such as 10× Chromium, Drop-seq, inDrop, and Seq-Well, use tiny droplets or wells to distribute cells, and are high-throughput and low-cost. Combinatorial indexing-based approaches, such as sciRNA-seq, can reverse transcribe and barcode mRNAs without physically isolating the cells (Ding et al., 2020). Single-cell transcriptomics technology opens a new avenue for molecular studies of tissues and organs from animals and plants, and avoids the lack of information on a particular cell type and heterogeneity from average data (Shulse et al., 2019; Mo and Jiao, 2022).

In plant research, single-cell transcriptomics has developed relatively late, but today it plays an essential and decisive role. Compared to traditional transcriptomics, the sequencing of single plant cells can resolve their heterogeneity. For example, when utilising traditional transcriptomics to profile drought-induced transcriptomic changes, researchers have observed the downregulation of growth-related, energy-consuming processes and the upregulation of stress defence genes at the organ level. However, single-cell transcriptome results suggested that, under mild drought stress, the mesophyll showed significant downregulation of genes, whereas most genes were upregulated in the epidermis. The response of plants to environmental stimuli is complex, and the key to understanding it depends on the study of networks at the single-cell level. The use of single-cell transcriptomics can help pinpoint tissue-specific pathways that respond to stress, and ultimately, these pathways can be specifically engineered in truly important tissues without causing unexpected side effects in other tissues or organs (Tenorio Berrio et al., 2022). There are two sequencing methods: one for sequencing protoplasts and the other for sequencing nuclei. However, the large and uncertain size of plant cells complicates single-cell sequencing. Problems persist even after removing the cell walls to obtain protoplasts, in addition to existing issues with the methods of protoplast extraction (Shulse et al., 2019). These problems can be solved using single nuclei for sequencing instead of whole cells, and the technique has a gene detection sensitivity similar to that of protoplast-based sequencing (van den Brink et al., 2017; Denisenko et al., 2020; Ding et al., 2020; Farmer et al., 2021). Both methods have the disadvantage of losing the spatial information of the cells when isolating them from tissues. Therefore, identifying ways to solve this problem is a widely discussed topic in single-cell transcriptomic research.

Spatial transcriptome technology has undergone considerable development in recent years. This provides a solution that compensates for the shortcomings of single-cell transcription analysis. Spatial transcriptome analysis differs from the method used for single-cell transcriptome analysis and faces specific difficulties. Compared to single-cell and traditional transcriptomics, spatial transcriptomics preserves spatial information while profiling the transcriptome, which provides more precise results for studying regulatory networks in plants. The data obtained from single-cell or traditional transcriptomic analyses are mapped back to the plant through cell sorting instead of using spatial data (Gurazada et al., 2021). During growth, differences in gene expression patterns at different locations in plant tissues lead to different phenotypes. Therefore, spatial transcriptomics provide a new perspective on plant growth and development.

Thus far, there has been a lack of studies to compare the existing methods for plant single-cell transcriptomics, specifically for isolating single cells and nuclei. Studies on spatial transcriptomics in plants are also rare. Here, we introduce and compare in detail the techniques used in related studies in recent years, explain the distinctions among different methods for extracting sequencing materials, and describe the differences between protoplast- and nucleus-based single-cell and spatial transcriptomics. Then, we put forth the challenges and future development directions and provide an idea for combining single-cell and spatial transcriptomics and spatial distribution analysis of elements in plant tissues to promote plant research.

2 Development of plant single-cell transcriptomics

Single-cell-related research has been conducted on animals, microorganisms, and plants. Microbiologists have used single-cell-related technologies to explore the differences in RNA accumulation between different parts of microorganisms (de Bekker et al., 2011). Zoologists have used this technology to study tumour cells (Ramsköld et al., 2012) and immune cells (Shalek et al., 2013) or to explore new methods for RNA sequencing (Grindberg et al., 2013). In recent years, plant single-cell transcriptomics has been developed. Researchers have conducted extensive plant studies using single-cell transcriptomic analyses. We reviewed the relevant literature in recent years and tabulated the results (Table 1). In 2015, Yan et al. published the first study on whole-genome sequencing of Zea mays using single plant cells (Li et al., 2015). Since then, the development of single-cell isolation and sequencing methods has made single-cell transcriptomics a popular topic in botanical research. In 2017, to explore the establishment of monoallelic gene expression in single plant cells, Jiao et al. developed a single-cell transcriptome sequencing protocol and applied it to Oryza sativa, which was the first report on the application of single-cell transcriptomics to plants (Han et al., 2017). To date, there have been four single-cell transcriptome studies of plant species: Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, and Zea mays, in the Single Cell Expression Atlas of the European Molecular Biology Laboratory (Conde et al., 2021).

TABLE 1
www.frontiersin.org

Table 1 An overview of plant single-cell transcriptomics related research in recent years.

From the published studies listed in Table 1, it is clear that most experiments were based on A. thaliana, and subsequently on O. sativa and Z. mays. Scientists typically use whole plant roots and root tips as experimental materials. A few scientists have also used other plant parts, such as leaves, stems, or germ cells, based on their research direction and purpose. Table 1 shows that most scientists have mainly constructed a single-cell atlas to explore the cell types of particular parts of plants, specifically roots, and examined developmental or differentiation trajectories combined with trajectory analysis, as transcriptomics can deliver information about gene expression.

Single-cell transcriptome technology has been used in numerous plant studies. In 2019, a group studied cell types from the primary root tissues of A. thaliana using a high-throughput Drop-seq approach. Their findings confirmed that RNA sequencing helps describe the developmental processes of plants (Shulse et al., 2019). It is worth mentioning that they also revealed that external stimuli, such as sucrose, can influence the developmental process, causing changes in sucrose response-related cells. In a similar case, in 2022, a group studied another external stimulus, salt stress, using cotton as their object and observed the influence of external stimuli on the profile and dynamic changes in gene expression (Liu et al., 2022b). For further research on the exploration of cellular heterogeneity, it is essential to compare gene expression profiles before and after external environmental stimuli by studying single-cell transcriptome sequencing results and constructing an atlas. In addition to investigating the environmental influences on plant cell gene expression, single-cell transcriptomics is a hot topic in the study of plant development and differentiation. One group has built models of cell differentiation within leaf tissues, studying stomatal lineage, and thereby determined a series of cellular programs related to tissue flexibility (Lopez-Anido et al., 2021).

With the help of single-cell transcriptomics, scientists have explored gene expression and signal transduction pathways while studying cell developmental trajectories. However, no research has been conducted on this topic. For example, researchers have observed how cytokinin signalling creates a link between the vascular perception of limited phosphate availability and epidermal responses (Wendrich et al., 2020). They drew an intersection between a single-cell atlas and target genes. They observed the distribution of transcription factor complexes in specific tissues, revealing the significance of the cytokinin signalling pathway. Comparing single-cell transcriptomic atlas and target genes in a critical signalling pathway, such as brassinosteroid signalling (Graeff et al., 2021), helps to understand the functions or effects of particular gene products or signal complexes in plants and their delivery or distribution in different tissues. Many scientists have used single-cell transcriptome techniques to create an atlas of plant species, which still have some genes that need to be explored, and have compared them with those of A. thaliana to study the functions of a few genes. This comparison reduces the difficulties and increases the hope of constructing a cell atlas for different plant species. However, studies involving target genes are limited, and even fewer studies have examined gene functions. In future studies, combining the analysis of an atlas with the target genes will be an analytical method that will attract increasing attention.

In addition to single-cell transcriptomics, there are many other omics and single-cell omics studies, and a few scientists have combined single-cell transcriptomic analysis with other omics analyses. For example, plant cells can be isolated for single-cell transcriptomics, which results in the loss of spatial information. Therefore, other omics techniques, such as spatial transcriptomics, are required to supplement missing data. In 2022, scientists applied scStereo-seq technology to plant research and constructed a single-cell spatial atlas (Xia et al., 2022). Owing to the combination of two omics technologies, they avoided a few disadvantages of a single omics technology and observed that the expression levels of related genes showed a gradient change trend in space. Another group combined single-cell transcriptomics technology with multi omics technology to study the molecular basis of developing Z. mays ears (Xu et al., 2021a). In addition to combining single-cell transcriptomics with spatial transcriptomics, some scientists have applied single-nucleus RNA sequencing and ATAC sequencing to plant roots, compared transcriptomic and epigenomic data at the single-cell level, and revealed how chromatin accessibility influences gene expression (Farmer et al., 2021). Using different omics technologies in combination with single-cell transcriptomics helps researchers solve problems associated with using only one technology, such as gene redundancy and spatial information.

As mentioned above, plant cells have cell walls, which make it challenging to obtain a single cell. Isolating single cells is affected by many factors, such as the parts or tissues of the plant used and the environmental conditions or genotypes of the plants (Shaw et al., 2021). Therefore, few researchers have been concerned with the methods used to perform single-cell transcriptomics. They introduced a commonly used method for isolating plant nuclei and demonstrated its universal applicability (Sunaga-Franze et al., 2021). A few researchers have introduced protoplast isolation methods for cotton (Liu et al., 2022b). New methods, not only for single-cell extraction but also for sequencing, have significance in plant research, reducing costs and obtaining enough single cells.

3 Standard methods for single-cell transcriptomics in plants

Sequencing of plant single-cell transcriptomes can be divided into sequencing of protoplasts and sequencing of nuclei. More accurately, the latter should be referred to as the single-nucleus transcriptome. Comparing the sequencing processes of the two, overall, the sequencing materials are all obtained from plant tissues, and then purified and separated for sequencing. Specifically, to protect the activity and integrity of protoplasts, the methods used are gentler and require higher freshness of plant samples during the process of obtaining protoplasts, whereas frozen or dried plant samples can be used to obtain nuclei. The following is a detailed introduction to both methods.

Based on recent protoplast-based single-cell transcriptomics studies (Table 1), it is clear that these procedures are similar (Figures 1, 2). First, the plants were grown based on the research goals. For example, to explore the single-cell transcriptome of rice roots under Cd stress, a small amount of rice was planted in both standard and Cd-stress environments. After a period of growth, the plant material is harvested from the required parts, including the root tips, whole roots, leaves, or other parts of the plant, and cut or crushed into small pieces to expand the area of the plant material exposed to the enzyme solution and shorten the time required for enzyme treatment. Plant materials are then mixed with the enzyme solution, called the digesting enzyme solution, and incubated under appropriate conditions to release protoplasts or plant cell nuclei. Plant materials must be mixed with an enzyme solution as soon as possible after being cut or crushed to maintain cell activity and improve the quality of protoplasts. This step is followed by further steps that involve filtration, centrifugation, and washing to obtain sufficient material for sequencing, relatively fewer undigested plant cells, and fewer broken cell organelles. Digestion helps remove the cell wall to isolate protoplasts, filtration helps remove substances that are not protoplasts or are larger in size, centrifugation helps concentrate protoplasts, and washing helps remove impurities and residual reagents. The quality and concentration of sequencing materials must be confirmed, as high quality and quantity are the basis for further sequencing to obtain sufficiently precise transcriptomic data (Liu et al., 2022b). A single-cell RNA sequencing (scRNA-seq) library is prepared and sequenced. Finally, the data are analysed to obtain the experimental results. This is the entire process of using single-cell transcriptomics; however, when it comes to each specific experiment, there still exist differences, particularly in the methods used to obtain ideal protoplast materials. A critical step in protoplast-based single-cell transcriptomics in plants is obtaining high-quality sequencing objects. Culture conditions, genotype, plant age, and other factors affect extraction methods to varying degrees, implying that when preparing protoplasts, the procedures must be adjusted based on the experimental conditions. In Table 2, the methods of partial studies are compared, that is, from studies listed in Table 1 or studies the authors referred to if they did not write the procedure for obtaining protoplasts.

FIGURE 1
www.frontiersin.org

Figure 1 Methods of Single-Cell Transcriptomics and Spatial Transcriptomics in Plants. The figure shows the general process for performing single-cell and spatial transcriptomics. Single-cell transcriptome technologies are divided mainly into protoplast-based and nucleus-based methods. Both approaches include collecting plant materials, separating protoplasts/nuclei through microfluidics after purification steps, i.e., filtration, centrifugation, and washing, confirming the viability and concentration of obtained protoplasts or cell nuclei, and sequencing and analysis. Spatial transcriptome methods differ from single-cell transcriptome methods and include embedding, cryosectioning, fixation, staining, and permeabilization steps. The plots a, b, c shown in the figure are visual diagrams of analysis results of protoplast-based sequencing, nucleus-based sequencing, and spatial transcriptomics sequencing respectively from published studies (Long et al., 2021; Zhang et al., 2021b; Liu et al., 2022a). The plots of single-cell transcriptomics, such as uniform manifold approximation and projection (UMAP) plots, present the cell classification results of sequencing, and those of spatial transcriptomics demonstrate the combination of transcriptome and spatial information.

FIGURE 2
www.frontiersin.org

Figure 2 Flow Chart for Protoplast Isolation. The figure shows the general process for protoplast isolation. First, obtain sequencing materials and cut or crush them into small pieces. Then, add the prepared enzyme solution. The amount and proportion of enzymes need to be appropriate. The reaction is carried out under appropriate conditions and enzymes digest plant cell walls to release protoplasts. This step is followed by further steps that involve filtration, centrifugation, and washing. Finally, the quality and quantity of the obtained protoplast solution are tested. If both meet the expectations, the preparation of the protoplast solution is completed.

TABLE 2
www.frontiersin.org

Table 2 Comparisons between methods used for protoplast-based single-cell transcriptome analysis in partially published studies.

First, the enzyme solution used during protoplasting differs among the different methods, but the types of reagents and enzymes used are similar in some methods. Enzymes are used to digest cell walls, whereas other reagents primarily serve to maintain the state of protoplasts or the efficiency of enzymes. Plant cell walls are mainly composed of cellulose, hemicellulose, and pectin, and the corresponding enzymes, including cellulase, pectolyase, and hemicellulose, are used to digest the cell walls (Somssich et al., 2016). However, most studies have used only cellulase and macerozymes, as macerozymes contain pectolyase and hemicellulase activities. A few scientists have added all four types of enzymes, that is cellulase, macerozyme, hemicellulose, and pectolyase, to the solution to obtain protoplasts from rice (Liu et al., 2021b). A few scientists have also added another enzyme, snailase, to reduce digestion time (Wang et al., 2022). Adding more enzymes is conducive to shortening the enzymatic hydrolysis time but can cause enzyme waste. The composition of cell walls differs among different plant species and organs. Therefore, determining the appropriate enzyme content and ratio is conducive for efficient enzymatic digestion. Bovine serum albumin (BSA), CaCl2, mannitol, and 2-(N-morpholino)ethanesulfonic acid (MES) are generally used to protect protoplasts. An enzyme solution with added BSA has the benefit of protecting the protoplasts, and a few researchers have used dithiothreitol (DTT) to perform the same function as BSA (Kim et al., 2021). Although BSA seems necessary to stabilise the enzyme solution, few researchers have used it (Xie et al., 2020). Salts such as CaCl2 are used to maintain the osmotic pressure of protoplasts in vivo and in vitro and prevent protoplasts from absorbing water and breaking. There do exist differences—for example, when studying A. thaliana and Z. mays, researchers are more likely to add KCl into the enzyme solution than when studying O. sativa as shown in Table 2. There are still differences in details within the same species, even when using the same parts of the materials in different studies. For example, when preparing rice protoplasts, a few studies added KCl to the digestion solution, whereas several other studies did not; however, mercaptoethanol was added (Han et al., 2017; Liu et al., 2021b; Wang et al., 2021c). NaCl has been added to enzyme solutions in a few studies (Xie et al., 2020). In addition, Table 2 shows differences among the types of mercaptoethanol used in solution: a few used 2-mercaptoethanol, others used 4-mercaptoethanol acid, and a few others used β-mercaptoethanol.

After determining the composition of the enzyme solution, reagents and enzymes were mixed. Usually, scientists dissolve cellulase, macerozyme, mannitol, MES, and certain enzyme solution reagents in water and heat the mix to activate enzymes in a warm bath at approximately 55 °C for nearly 10 min. Upon heating, the turbid enzyme solution becomes clear. After heating, it is cooled to room temperature, i.e., approximately 25 °C, and the rest of the materials, including BSA, CaCl2, and mercaptoethanol, are added before adding plant materials. Unlike the general methods of mixing plant materials and the entire digestion solution, a few techniques involve mixing plant materials and the enzyme solution and incubating the mixture in the dark for the reaction. The remaining chemicals, including NaCl, CaCl2, KCl, and MES, are added for further enzymatic reactions to release protoplasts (Yoo et al., 2007; Xie et al., 2020; Gala et al., 2021).

After mixing the enzyme solution with plant materials, it is essential to provide sufficient time and appropriate reaction conditions to ensure the release of protoplasts. Most studies listed in Table 2 did not specify the proportions of plant materials and enzyme solutions. Table 2 lists the incubation conditions. A few groups used a vacuum pump to isolate protoplasts in experiments to reduce the digestion reaction time, whereas others did not. Therefore, the reaction times of others were longer than reaction times of those that used the vacuum pump. This is because negative pressure promotes the infiltration of the enzyme solution into plant materials, such that the contact surface between the enzyme and plant cells increases, thereby reducing the reaction time. In addition, it is helpful to incubate a mixture of enzyme solution and plant materials in the dark or to shake this mixture using a machine or manually to reduce the reaction time. When shaking the enzyme solution, it is necessary to shake it gently at a speed of 40–80 rpm to avoid cell fragmentation caused by violent shaking. However, a few groups used 200 rpm shaking speed (Denyer et al., 2019). The reaction times used in different studies to release protoplasts are similar, but there are still differences in the details. For example, when digesting the leaves of A. thaliana, the digestion time used by a few groups was 30 min, whereas others chose to digest for 2 h, probably because the leaves they used were different; that is, leaves from the leaf base regions or the third leaf (Liu et al., 2022c; Tenorio Berrio et al., 2022). Although different digestion times are listed in Table 2, it is only for reference—digestion times can differ based on differences in the digestion conditions, species, and parts of plants used, and it is necessary to check whether most of the plant cell walls are digested at least every half an hour. If measures are taken to shorten the time, the enzyme content in the enzyme solution will be higher, and if the plant tissues used are young or have a short growth time, the digestion time will be shorter.

After protoplast release, cells are filtered, centrifuged, and washed. Filtration aims to remove undigested plant cells and obtain protoplasts; therefore, the filter must be of the proper size. The size of plant protoplasts is small, e.g., A. thaliana mesophyll protoplast diameter ranges from 30 to 50 μm, and therefore, most groups used 30–80 μm nylon mesh or cell strainer to filter protoplasts (Yoo et al., 2007). A few groups filtered the protoplasts only once, whereas other groups filtered them twice or more times, and other groups filtered the protoplasts through two different sizes of filters to obtain more protoplasts and less undigested plant cells (Wang et al., 2021c; Zhang et al., 2021b; Liu et al., 2022c). Centrifugation helps concentrate protoplasts, and the centrifugation conditions differ among different plant species. Table 2 lists the centrifugation conditions used in various studies. Centrifugation speed and time depend on the characteristics of the plant materials, such as their fragility, type, and the number of cells produced during digestion (Bargmann and Birnbaum, 2010). High centrifugation speed can increase both the recovery rate of protoplasts and the risk of cell disruption (Yoo et al., 2007). As shown in Table 2, when centrifuging rice protoplasts, the centrifugation speed (e.g., 100 g, 130 g, or 200 g) is slower and the centrifugation time (e.g., 3 min or 5 min) is shorter than that used for centrifuging A. thaliana (e.g., 500 g for 5–10 min). The differences in centrifugation speeds may result from the different characteristics of the protoplasts of these two plants or from differences in the isolation steps used in published studies. During centrifugation, a few groups maintained the solution at constant low-temperature conditions, such as 4 °C or −4 °C, to protect protoplasts from damage (Denyer et al., 2019; Satterlee et al., 2020; Xu et al., 2021a). It is also possible to collect small protoplasts by centrifugation at low speeds and then increase the speed of collecting larger protoplasts, avoiding the breakage of small protoplasts and the loss of large protoplasts caused by insufficient centrifugation speeds. Washing is also necessary because it can help remove specific ions or tissue residues that affect the downstream steps, such as undigested plant tissues and Mg2+ ions (Liu et al., 2022d). During the preparation of the washing solution, the entry of impurities is avoided, and an appropriate osmotic pressure amenable to protoplasts is maintained in the washing solution. Scientists have typically used 8% mannitol or a protoplast solution without enzymes for washing; however, a few scientists have also used W5 solution, WB buffer, or other washing buffers (Table 2). It is worth mentioning that the components of W5 solutions prepared by different groups showed little difference (Chen et al., 2021b; Xu et al., 2021b). After washing, the protoplasts were centrifuged again to separate the washing solution from the protoplasts. The centrifugation speed and time can be the same as those used in the centrifugation step, or at a slower rate, shorter time, or both. The protoplasts are resuspended in mannitol, washing solution, or other solutions based on the sequencing protocol, which helps adjust the concentration of protoplasts and avoid influencing reverse transcription reactions (Chen et al., 2021b). The viability and concentration of the protoplasts are then confirmed for better sequencing results. Trypan blue, 4′,6-diamidino-2-phenylindole, fluorescein diacetate, or acridine orange-propidium iodide are used to select protoplasts, and a haemocytometer is used to measure the protoplast concentration. The viability and concentration of protoplasts need to be adjusted to reach the standard of sequencing methods, and the concentration of protoplasts can be adjusted using a resuspension solution. Next, an scRNA-seq library is prepared and sequenced according to different protocols.

The process of isolating plant nuclei is similar to that used to separate plant protoplasts. The plant material used to extract nuclei can be either fresh or frozen. Plant tissues are cut or crushed into small pieces and mixed with a nuclei extraction solution. The solutions used for nuclei extraction in the relevant studies are listed in Table 3; there are similarities, with a few differences. Compared to isolating protoplasts, plant nuclei isolation requires the lysis of plant cell membranes and protection of nuclear membranes; therefore, the solutions are relatively more destructive. The incubation is usually performed on ice. After incubation, filtration, centrifugation, and resuspension are conducted. The speed of centrifugation for isolating plant nuclei is usually faster than 1,000 g; for example, a few groups centrifuge at 2,000 g, whereas others centrifuge at 25,000 g (Giuliano et al., 1988; Dorrity et al., 2021). Plant nuclei are washed, recovered, and resuspended in a buffer solution. Some groups have conducted Percoll gradient centrifugation for impurity removal and nuclear recovery (Dorrity et al., 2021). The remaining steps are performed according to the nucleus-based sequencing protocols of the sequencing platforms.

TABLE 3
www.frontiersin.org

Table 3 Nuclei extraction solutions.

4 Development of plant spatial transcriptomics

Spatial transcriptomic technologies originated in mammalian systems and are widely used in mammals. A few published studies have applied omics technology to plant research. However, more strategies are available for mammalian spatial transcriptomics than for plants, mainly because of the differences between animal and plant cell structures such as cell walls. Currently, research on plant spatial transcriptomics is divided primarily into two categories: finding ways to use and optimise spatial transcriptomic methods, and using spatial technology to solve problems in plant research.

In 2017, Stefania and his group provided a method for generating and studying high-resolution and spatially resolved functional profiles in plants (Giacomello et al., 2017; Giacomello and Lundeberg, 2018). At the method level, research has not been limited to methodology in plants or the application of methods used in mammalian systems to plants, but has also focused on finding ways to improve spatial transcriptome technology in plant and mammalian systems. Spatial transcriptome analysis requires high spatial resolution and high spatial transcriptome data, such as maps in daily life, and provides an extensive range of accurate information on buildings located in an area. A few analytical techniques, such as single-cell transcriptomics, produce high-quality data, implying that they provide relatively high spatial resolution but lose spatial information. In addition, sorting or isolation technologies, such as fluorescence-activated cell sorting (FACS), isolation of nuclei tagged in specific cell types (INTACT), and LCM, provide spatial resolution or cell-specific expression data. However, these methods have certain limitations. Plant materials must be transgenic using FACS and INTACT; however, this is a challenge for a few plant species. In addition, many factors such as cell type, yield, and purity affect the quality of the LCM and FACS methods. Technologies such as FACS require the digestion of plant cells to isolate protoplasts, which causes loss of spatial information (Giacomello et al., 2017; Gurazada et al., 2021). Fluorescent in situ hybridisation (FISH) is an optimal option for obtaining spatial information from plant cells; however, it can only detect one gene at a time, limiting its speed (Tirichine et al., 2009; Gurazada et al., 2021). Therefore, spatial transcriptomics technology needs to be improved. Currently, there are two main types of spatial transcriptome technology: next-generation sequencing (NGS)- and imaging-based approaches. NGS-based methods involve high-throughput sequencing based on scRNA-seq technology and a spatial barcode, whereas imaging-based methods include in situ sequencing (ISS) and in situ hybridisation (ISH). Recent research on spatial transcriptomic technologies, including ISH-based MERFISH and ISS-based STARmap, has promoted genetic and genomic studies (Moffitt et al., 2018; Wang et al., 2018; Liu et al., 2022a). Spatial indexing approaches, including 10× Visium, use NGS-based methods to quantify the gene expression profiles formed by the local hybridisation of barcodes and RNA molecules, and the spatial information is stored by barcodes. These approaches can detect the entire transcriptome in an unbiased manner without knowing the target gene in advance but must ensure the integrity of mRNA in the tissues. ISS- and ISH-based methods can intuitively analyse data at organisational spatial locations, but mRNA abundance in situ is low and often degraded (Elhanani et al., 2023). Although spatial transcriptomic methods face several challenges, including resolution and sensitivity, researchers are constantly attempting to overcome these difficulties (Rao et al., 2021).

Spatial transcriptomics technologies provide new insights into plant growth, development, and molecular biology. A few researchers have mapped and visualised the related gene expression in C4 and crassulacean acid metabolism in Portulaca at the spatial level using the 10× Genomics Visium spatial transcriptomics platform (Moreno-Villena et al., 2022). One group combined spatial and temporal transcriptome information to explore the early development of tomatoes (Zhang et al., 2016). In 2022, a Chinese group studied Phalaenopsis Big Chili using 10× Visium technology (Liu et al., 2022a). A few studies have analysed transcriptome data at the spatial level; however, in these studies, researchers collected plant materials in spatial order and then sequenced them separately, which preserved the spatial information and provided help in the analysis, but the information lacked high resolution (Guo et al., 2021).

Spatial transcriptome technology is helpful, but there are still few studies using this technology in plant research owing to many reasons and challenges that will be discussed later.

5 Methods used in plant spatial transcriptomics

In plant research, technology for preserving spatial information is not widely used. This is partly because appropriate methods have not been used in plants, and the structures of plant cells differ from those of animal cells, such as the presence of cell walls and vacuoles. In a few studies, the processes involved collecting plant materials in spatial order, separating leaves from outer to inner, sorting root and aerial tissues, and sequencing them separately (Guo et al., 2021; Apelt et al., 2022). Such methods are usually not highly resolved because they involve the sequencing of many mixed cells, providing a considerably limited variety of spatial locations. One such process involves dissecting different tissues for transcriptome analysis, such that the spatial information for the entire section is maintained, but the spatial information is distinguished according to the initially selected regions (Stahl et al., 2016).

A few scientists achieved high-resolution gene expression data, even near single-cell resolution, using high-throughput methods. In general, these methods include collecting targeted sequencing materials, embedding, cryosectioning, fixation, staining, imaging, permeabilisation, and sequencing, according to the manufacturer’s instructions (Figure 1) (Giacomello et al., 2017; Giacomello and Lundeberg, 2018; Kivivirta et al., 2019; Liu et al., 2022a; Moreno-Villena et al., 2022). Embedding is usually the first step in spatial transcriptomic experiments after harvesting the targeted plant materials. The samples are embedded in cold optimal cutting temperature (OCT) compound to support the plant tissues during cryosectioning. Plant tissues have structures such as hard cell walls and large vacuoles with high water content, which increase the hardness of plant tissues after conventional freezing, making it difficult to obtain high-quality slices. Then, the plant samples are sliced at the proper temperature (−15 to −20 °C). The thickness of the sections must also be appropriate to avoid breakage and overlapping of the different tissue sections during slicing. In addition, the organisational structure of plant samples is complex, the degree of lignification varies, and the relative RNA content of the tissue is relatively low; therefore, the thickness of plant slices generally depends on the type of sample (Giacomello and Lundeberg, 2018; Liu et al., 2022a). Subsequently, the plant sections are fixed on the array at approximately 37 °C. The solutions used in the fixation step differed among different studies; a few studies used neutral formaldehyde, whereas others used methanol. The dyes used in the staining step also differ among various studies, such as toluidine blue or haematoxylin-eosin staining, and the staining time is usually 1 min, followed by imaging of the stained sections. Permeabilisation is a critical step in spatial transcriptomics following imaging. In this step, it is necessary to determine the optimum permeabilisation conditions for plant samples to release sufficient mRNA (Moreno-Villena et al., 2022). A few experiments have also included a pre-permeabilisation step before permeabilisation for partially denaturing proteins to improve the efficiency of permeabilisation (Giacomello and Lundeberg, 2018). After permeabilisation, the samples were sequenced according to the manufacturer’s protocol, including reverse transcription, tissue removal, and probe release. Finally, the water used in these experiments must be nuclease-free in specific steps to avoid degradation of the cDNA-mRNA hybrids.

Although differences exist between different methods and there are only a few relevant studies in plant research, we can learn from the methods used in mammalian systems or follow the manufacturer’s instructions to optimise these methods in plant systems (Giacomello and Lundeberg, 2018).

6 Comparisons among protoplast-based and nucleus-based single-cell transcriptomics and spatial transcriptomics

Although protoplasts and plant nuclei can be used for single-cell sequencing, there are differences, and both have advantages and disadvantages. Comparing these two sequencing methods and the scope of their application is beneficial for attaining a deeper understanding of single-cell transcriptomics in plant systems and for exploring future research directions.

As shown in Tables 2, 3, enzyme solutions, centrifugation speeds, and other parameters differ between the two methods. For example, because of differences in sequencing materials, these two methods use different solutions to digest or isolate plant cells. When separating protoplasts, researchers usually use enzymes to digest the cell walls of plant cells and add various reagents to protect the cell membrane and maintain protoplast shape. In contrast, when isolating nuclei, researchers usually do not use enzymes, but instead use other reagents, such as Triton X-100, to separate plant nuclei, and the cell membranes need to be destroyed (Sikorskaite et al., 2013).

When separating protoplasts, artefacts affecting the transcriptome are introduced by the addition of enzymes. Scientists have confirmed that protoplasts affect transcriptomics (Birnbaum et al., 2003; Denisenko et al., 2020; Conde et al., 2021; Sunaga-Franze et al., 2021). The addition of enzymes necessitates an appropriate incubation temperature, such that the enzymes are in the most suitable state for digestion. However, cellular machinery is also active at the appropriate temperature for enzymatic action, leading to alterations in gene expression. A few scientists have found that when obtaining single cells from animals, the addition of enzymes at the incubation temperature can lead to a dramatic increase in gene expression after incubation at 37 °C (Adam et al., 2017; van den Brink et al., 2017; Potter, 2018; Denisenko et al., 2020). Similar results were obtained when plant cells were digested to obtain protoplasts. Isolating plant nuclei for single-cell transcriptomics would not face this challenge, as no enzymes are used and incubation is performed on ice, avoiding the activation of gene expression. Therefore, from the perspective of introducing artefacts and anomalous gene expression, plant nucleus-based single-cell transcriptomics sequencing is more precise, convenient, and applicable than protoplast-based single-cell transcriptomics. Researchers have also developed methods to overcome the limitations of protoplast release, such as performing independent bulk RNA sequencing experiments to eliminate the effect of enzyme addition or using transcription inhibitors (Potter, 2018; Denyer et al., 2019; Sunaga-Franze et al., 2021). Specifically, α-amanitin inhibits transcription via RNA polymerase II, thereby preserving gene expression patterns (Potter, 2018). Some researchers have extracted RNA from protoplast and unprotoplast plant tissues, identified genes induced by protoplast, and removed them from analyses (Denyer et al., 2019; Jean-Baptiste et al., 2019; Shulse et al., 2019). Protoplast-based single-cell transcriptome methods are not applicable, whereas nucleus-based techniques can. Protoplasts are used for plant single-cell transcriptomics; however, the methods for isolating protoplasts differ for different plant materials, even for the same plant and their exact position in the plant. As mentioned above, this is mainly because many factors affect the enzyme solution and incubation conditions adopted in the experiments. The main components of cell walls are similar, but differences in the functions of cells result in the production of distinct cell wall polymers surrounding each cell. Therefore, cell walls are considered heterogeneous when comparing plant cells of different types and developmental stages (Somssich et al., 2016). Thus, the methods for obtaining protoplasts lack applicability, leading to limitations in protoplast-based single-cell transcriptomics. In addition, these protoplast-based methods rely heavily on the quality and quantity of protoplasts and their compatibility with sequencing systems (Thibivilliers et al., 2020). Protoplast suspensions need to contain more intact protoplasts and less or no cell debris and damaged cells to avoid low cell numbers and quality, mRNA cell leakage, and obstruction of microfluidic systems, thereby ensuring that the microfluidic systems run correctly (Conde et al., 2021). Nucleus-based single-cell transcriptomics can overcome these problems owing to the introduction of a general method for plant nucleus isolation, which is an alternative to protoplast isolation methods and has been tested for its applicability to several different types of plants (Sunaga-Franze et al., 2021).

The applicability of these two procedures is not only reflected in isolation techniques, but also in the requirements for cell sequencing, such as the size of sequencing materials or freshness of sequencing samples (Hwang et al., 2018; Denisenko et al., 2020; Nadelmann et al., 2021). Generally, to avoid obstruction of sequencing systems, the sequencing platform restricts the cell size. The protoplast was larger than the nucleus, which increased the possibility of obstruction. The freshness of protoplasts is also critical for sequencing, and it is necessary to maintain their shape to avoid membrane breakage. However, nuclei can be obtained from any plant cell, regardless of its freshness and size, and even from frozen plant tissues (Krishnaswami et al., 2016; Lake et al., 2016; Bakken et al., 2018; Wu et al., 2019; Nadelmann et al., 2021). Protoplast-based sequencing is a stringent requirement for sequencing materials.

Nuclei seem to be more appropriate for single-cell transcriptomics than protoplasts; however, isolating nuclei from plant cells results in the loss of cell cytoplasm, chloroplasts, and mitochondria, implying that single-nucleus RNA-seq (snRNA-seq) usually contains less transcription information. Therefore, snRNA-seq has limitations in obtaining information from particular cell types because it only uses a part of the cell, that is, the nucleus, instead of an entire cell for transcriptional profiling (Hu et al., 2018; Maitra et al., 2021). The cell nucleus contains only 10–20% of all cellular transcripts; therefore, nuclear RNA does not entirely reflect the RNA from the whole cell, which may increase the difficulties in the analysis (Bakken et al., 2018; Potter, 2018; Nadelmann et al., 2021). Moreover, snRNA-seq excludes critical gene information from outside the nucleus and does not completely represent the single-cell transcriptome. From this point of view, protoplasts are more accurate than plant nuclei for single-cell transcriptome analysis. However, in research on frozen samples, precious frozen samples, complex plant tissues that are difficult to digest, and those that are easily influenced by gene expression changes during protoplast digestion, plant nuclei are still more appropriate for sequencing.

The most significant difference between single-cell and spatial transcriptomics is that the former is highly resolved but loses spatial information, whereas the latter preserves spatial information but fails to achieve precision. Therefore, if more detailed research is preferred, such as further exploration of cell types and their molecular characteristics in response to different environmental stimuli, single-cell transcriptomics is often used. If the characteristics of different spatial regions or the molecular characteristics of cells in a certain region of plant organs are preferred, for example, by constructing a spatiotemporal atlas of the development of floral organs, spatial transcriptomics is often used. Microfluidic technology is typically used to separate single cells in single-cell transcriptomics. Microarray technology is typically used to store spatial information in transcriptomics. Each spot on the microarray has millions of oligonucleotide probes. An oligonucleotide probe sequence used in spatial transcriptomics usually has a spatial barcode to store spatial information, whereas the barcode in single-cell transcriptomics is used to distinguish different cells (Gurazada et al., 2021).

7 Challenges and perspectives

Single-cell transcriptomics is a mature analytical technology in animal research, but a rising and naive front field in plant research, which still needs further exploration and development. Earlier sections revealed problems in both protoplast-based and plant nucleus-based single-cell transcriptomics, and these difficulties need to be urgently solved. In contrast, this is the era of multi omics, and finding ways to combine single-cell transcriptomics with other omics is noteworthy. Integrating single-cell transcriptome analysis with other single-cell technologies or omics technologies may provide a list of cell types, and more importantly, new insights into the regulatory logic and spatial organisation among cells (Stuart et al., 2019).

Plant materials contain different cell types, a few of which are resistant to digestion owing to their secondary cell walls, implying that, based on currently used methods, such cells pose challenges in protoplast isolation, causing a biased analysis of cell types (Bezrutczyk et al., 2021; Conde et al., 2021; Farmer et al., 2021). The different positions of cells in plant materials also introduce differences in their area exposed to digestive enzymes; for example, a few cells located inside the plant materials and surrounded by other cells are challenging to digest. These cells are readily missed during the analysis because they are removed during filtration. A few groups have discovered that central stele cells are more difficult to capture than epidermal cells; thus, the preparation of plant tissue protoplasts is biased for cells located on the outer surface of plant materials and sufficiently exposed to digestion enzymes (Denyer et al., 2019; Jean-Baptiste et al., 2019; Shulse et al., 2019; Farmer et al., 2021). In addition, a few groups have observed that not only the position of cells, but also their developmental stage influences digestion, as younger cells are more readily digested than mature cells, leading to increased representation in the analysis (Jean-Baptiste et al., 2019; Ryu et al., 2019). Therefore, a few general and effective methods are needed to isolate protoplasts that can avoid errors caused by cell types, relative positions of cells, or developmental stages.

In addition, single-cell sequencing technology requires improvement to achieve more precise sequencing results. For example, droplet-based methods such as the 10× Chromium method help generate an RNA-seq library. In this method, a single cell or nucleus is encapsulated inside a gel bead that can supply a barcoded oligonucleotide for reverse transcription. The droplet-based method can be used to analyse thousands of cells in a single experiment. However, in a few cases, there is more than one cell or nucleus inside a droplet or a suboptimal number of beads inside a droplet, thereby causing problems in the analysis. One group tested the probability of this event and observed that approximately 20% of the droplets contained either no beads or more than one bead in their experiments (Lareau et al., 2020; Sunaga-Franze et al., 2021). Nanowell-based systems, such as Takara iCELL8, can also generate an RNA-seq library and do not suffer from this problem because they capture single cells or nuclei in a nanowell. Researchers can check the number of cells or nuclei by using microscopy. However, the chips used in this method have only a fixed number of nanowells, resulting in limited scalability (Sunaga-Franze et al., 2021). Therefore, it is important to improve the relevant technology to ensure that one droplet contains one bead and one cell or nucleus and provide better scalability.

In single-cell transcriptomics, cells are digested and isolated from their original position in the tissue, and their spatial information is lost. Scientists have successfully mapped scRNA-seq data by inferring cellular sources and locations from scRNA-seq data and in situ RNA patterns (Achim et al., 2015; Satija et al., 2015). Scientists have also constructed a 3D atlas in plant research by combining scRNA-seq and microscopy-based 3D spatial reconstruction (Neumann et al., 2022). Recent studies have used spatial transcriptomics to preserve spatial information. Single-cell transcriptome analysis can be strengthened by combining it with a spatial transcriptome analysis. The cell types of whole plant tissues can be classified using single-cell transcriptome analysis, and the spatial information of whole tissues can be obtained using a spatial transcriptome. Researchers have established an in situ single-cell spatial transcriptome method in plants based on stereo-seq and applied it to study A. thaliana leaves (Chen et al., 2022; Xia et al., 2022). Their method differed from the single-cell transcriptome methods described, but was similar to spatial transcriptome methods. They sectioned frozen plant tissues and adhered these plant sections to a stereo-seq chip for downstream operations. There are a few other integration strategies, and Longo et al. identified two main types of algorithms: deconvolution and mapping (Longo et al., 2021). Methods based on deconvolution algorithms separate discrete cellular subpopulations from the mRNA transcript mix at each capture point based on single-cell data. SPOTlight is such a method, integrating single-cell and spatial transcriptome, centering on a seeded regression using non-negative matrix factorisation, initialising with cell type marker gene, and non-negative least squares, and deconvoluting capture locations (Elosua-Bayes et al., 2021). There are other algorithms for deconvoluting spatial information, such as dampened weighted least-squares, negative binomial distribution, and Poisson distribution models (Andersson et al., 2020; Dong and Yuan, 2021; Cable et al., 2022). Some scientists believe that single-cell transcriptome data and spatial transcriptome data follow a particular probability distribution, and based on these hypotheses, a few probabilistic models have been proposed to help integrate these two types of data. Mapping methods combine these two types of transcriptome data, mapping the designated scRNA-based cell subtypes to each cell on a high-plex RNA imaging map, and mapping each scRNA-seq cell to a specific niche or region of tissue. Similar to deconvolution methods, a few mapping methods use probabilistic models, such as the variation Bayesian mean-field approximation (Qian et al., 2020). In addition, other methods, such as CellTrek, combine single-cell and spatial transcriptome data using co-embedding and metric learning approaches (Wei et al., 2022).

Similar to single-cell transcriptomics, research on spatial transcriptomics is more prevalent in mammals than in plants, partially because of differences between mammalian and plant cells. Different plant cell structures and organelles, including cell walls, vacuoles, chloroplasts, and the presence of secondary metabolites, make it challenging to obtain spatial transcriptomic data from plant cells or lead to cryosectioning problems (Bourgaud et al., 2001; Cosgrove, 2005; Giacomello et al., 2017; Gurazada et al., 2021). Specifically, owing to the cell walls, plant cell sections can neither be extremely thick nor extremely thin to avoid reducing mRNA release or capture and causing sectioning artefacts or low transcript levels (Gurazada et al., 2021). In addition, because a few plant tissues, such as leaves and roots, have curvatures, many sections are required to obtain complete information. Currently, there are few solutions to the challenges caused by these differences, and it is often necessary to optimise the method for specific plant species or tissues (Giacomello et al., 2017; Giacomello and Lundeberg, 2018; Gurazada et al., 2021). Furthermore, challenges have been introduced by spatial transcriptome technology, including resolution, sensitivity, throughput, and accessibility (Rao et al., 2021).

In addition, the related applications of single-cell transcriptomics and spatial transcriptomics can be explored further, whether using one of these two techniques or in combination with other analytical techniques. For example, although there has been comprehensive research on the root cell types of plants such as A. thaliana and O. sativa, there is almost no relevant research on the cross-species analysis of roots at the single-cell level. This cross-species analysis may provide a basis for analysing differences in elemental absorption among different plants. Relevant animal studies have identified a new cell subtype, defined signature genes, and revealed differences in drug absorption among species through cross-species analysis (Li et al., 2022). Consideration should be given to the spatial distribution of elements in plants to capture thorough spatial information. Combining single-cell and spatial transcriptome analyses with the distribution patterns and speciation of one or more plant tissue elements offers several opportunities, particularly for studying the effects of stress caused by these elements. In recent years, techniques for analysing elemental distribution and speciation in plant materials have been developed, making it possible to draw an element-based map with high resolution (Lu et al., 2013; Lu et al., 2017). The combination of these techniques will provide new insights into plant research. For example, a combination of these three techniques could be used to explore the effects of different elemental stresses on plant roots. Analysis of the elemental distribution would reveal the characteristics of the absorption of this element and its storage and transportation mechanisms in plant roots. Spatial transcriptome and single-cell transcriptome analyses help analyse the effect of this element on gene expression in plant roots, or the difference in absorption and storage caused by the difference in gene expression, to provide a gene reference for cultivating low-accumulation plant varieties.

8 Conclusion

Both single-cell and spatial transcriptome analyses are helpful but still developing technologies in plant research, although they are relatively mature in mammalian systems. Single-cell transcriptomics provides insights into each cell type in plant tissues, such as classifying the cell type and analysing the cell developmental trajectory, and there is likely to be more research investigating functional genes in the future. Some problems in plant single-cell transcriptomics are common to single-cell transcriptomics in all organisms, such as those in sequencing techniques, whereas other difficulties are unique to plant research, including the lack of standard methods for isolating protoplasts. Owing to differences in cellular features and culture conditions, digestion solutions, and incubation conditions for releasing protoplasts differ among studies, leading to a lack of applicability. The snRNA-seq method, which has several advantages over scRNA-seq using protoplasts, can overcome this problem; however, it lacks information outside the nuclei. Spatial transcriptome technology provides spatial information on targeted tissues; however, it fails to reach single-cell resolution. Therefore, in the future, the single-cell transcriptome technique requires more effort to increase the applicability of methods for releasing protoplasts, and the spatial transcriptome technique requires more effort to achieve high resolution. More importantly, the application of single-cell transcriptomics and spatial transcriptomics needs to be explored, such as cross-species analysis using single-cell transcriptomics in the roots of different plants, and the combination of single-cell analysis and other analyses, such as spatial transcriptome analysis and spatial element distribution analysis, also needs to be emphasised to compensate for their limitations and create ample possibilities for research.

Author contributions

CC and LL conceived the manuscript. CC wrote the manuscript. YG and LL reviewed and revised the manuscript. All the authors contributed to the manuscript and agreed to the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (22276164), the Key Research and Development Project of Science and Technology Department of Zhejiang Province (2023C02020), and the National Natural Science Foundation of Zhejiang Province (LZ22C150004).

Acknowledgments

We thank Editage for their assistance with English language editing.

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

Achim, K., Pettit, J. B., Saraiva, L. R., Gavriouchkina, D., Larsson, T., Arendt, D., et al. (2015). High-throughput spatial mapping of single-cell RNA-seq data to tissue of origin. Nat. Biotechnol. 33 (5), 503–509. doi: 10.1038/nbt.3209

PubMed Abstract | CrossRef Full Text | Google Scholar

Adam, M., Potter, A. S., Potter, S. S. (2017). Psychrophilic proteases dramatically reduce single-cell RNA-seq artifacts: a molecular atlas of kidney development. Development 144 (19), 3625–3632. doi: 10.1242/dev.151142

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, A., Bergenstrahle, J., Asp, M., Bergenstrahle, L., Jurek, A., Fernandez Navarro, J., et al. (2020). Single-cell and spatial transcriptomics enables probabilistic inference of cell type topography. Commun. Biol. 3 (1), 565. doi: 10.1038/s42003-020-01247-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Apelt, F., Mavrothalassiti, E., Gupta, S., Machin, F., Olas, J. J., Annunziata, M. G., et al. (2022). Shoot and root single cell sequencing reveals tissue- and daytime-specific transcriptome profiles. Plant Physiol. 188 (2), 861–878. doi: 10.1093/plphys/kiab537

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakken, T. E., Hodge, R. D., Miller, J. A., Yao, Z., Nguyen, T. N., Aevermann, B., et al. (2018). Single-nucleus and single-cell transcriptomes compared in matched cortical cell types. PloS One 13 (12), e0209648. doi: 10.1371/journal.pone.0209648

PubMed Abstract | CrossRef Full Text | Google Scholar

Bargmann, B. O., Birnbaum, K. D. (2010). Fluorescence activated cell sorting of plant protoplasts. J. Vis. Exp. (36), e1673. doi: 10.3791/1673

CrossRef Full Text | Google Scholar

Bezrutczyk, M., Zollner, N. R., Kruse, C. P. S., Hartwig, T., Lautwein, T., Kohrer, K., et al. (2021). Evidence for phloem loading via the abaxial bundle sheath cells in maize leaves. Plant Cell 33 (3), 531–547. doi: 10.1093/plcell/koaa055

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (5652), 1956–1960. doi: 10.1126/science.1090022

PubMed Abstract | CrossRef Full Text | Google Scholar

Bourgaud, F., Gravot, A., Milesi, S., Gontier, E. (2001). Production of plant secondary metabolites: a historical perspective. Plant Sci. 161 (5), 839–851. doi: 10.1016/s0168-9452(01)00490-3

CrossRef Full Text | Google Scholar

Cable, D. M., Murray, E., Zou, L. S., Goeva, A., Macosko, E. Z., Chen, F., et al. (2022). Robust decomposition of cell type mixtures in spatial transcriptomics. Nat. Biotechnol. 40 (4), 517–526. doi: 10.1038/s41587-021-00830-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, A., Liao, S., Cheng, M., Ma, K., Wu, L., Lai, Y., et al. (2022). Spatiotemporal transcriptomic atlas of mouse organogenesis using DNA nanoball-patterned arrays. Cell 185 (10), 1777–1792.e1721. doi: 10.1016/j.cell.2022.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Lv, Y., Yin, X., Chen, X., Chu, Q., Zhu, Q. H., et al. (2021a). Effects of Sample Size on Plant Single-Cell RNA Profiling. Curr. Issues Mol. Biol. 43 (3), 1685–1697. doi: 10.3390/cimb43030119

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Tong, S., Jiang, Y., Ai, F., Feng, Y., Zhang, J., et al. (2021b). Transcriptional landscape of highly lignified poplar stems at single-cell resolution. Genome Biol. 22 (1), 319. doi: 10.1186/s13059-021-02537-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Conde, D., Triozzi, P. M., Balmant, K. M., Doty, A. L., Miranda, M., Boullosa, A., et al. (2021). A robust method of nuclei isolation for single-cell RNA sequencing of solid tissues from the plant genus Populus. PloS One 16 (5), e0251149. doi: 10.1371/journal.pone.0251149

PubMed Abstract | CrossRef Full Text | Google Scholar

Cosgrove, D. J. (2005). Growth of the plant cell wall. Nat. Rev. Mol. Cell Biol. 6 (11), 850–861. doi: 10.1038/nrm1746

PubMed Abstract | CrossRef Full Text | Google Scholar

de Bekker, C., Bruning, O., Jonker, M. J., Breit, T. M., Wosten, H. A. (2011). Single cell transcriptomics of neighboring hyphae of Aspergillus niger. Genome Biol. 12 (8), R71. doi: 10.1186/gb-2011-12-8-r71

PubMed Abstract | CrossRef Full Text | Google Scholar

Denisenko, E., Guo, B. B., Jones, M., Hou, R., de Kock, L., Lassmann, T., et al. (2020). Systematic assessment of tissue dissociation and storage biases in single-cell and single-nucleus RNA-seq workflows. Genome Biol. 21 (1), 130. doi: 10.1186/s13059-020-02048-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Denyer, T., Ma, X., Klesen, S., Scacchi, E., Nieselt, K., Timmermans, M. C. P. (2019). Spatiotemporal Developmental Trajectories in the Arabidopsis Root Revealed Using High-Throughput Single-Cell RNA Sequencing. Dev. Cell 48840–52 (6), e845. doi: 10.1016/j.devcel.2019.02.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, J., Adiconis, X., Simmons, S. K., Kowalczyk, M. S., Hession, C. C., Marjanovic, N. D., et al. (2020). Systematic comparison of single-cell and single-nucleus RNA-sequencing methods. Nat. Biotechnol. 38 (6), 737–746. doi: 10.1038/s41587-020-0465-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, R., Yuan, G. C. (2021). SpatialDWLS: accurate deconvolution of spatial transcriptomic data. Genome Biol. 22 (1), 145. doi: 10.1186/s13059-021-02362-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorrity, M. W., Alexandre, C. M., Hamm, M. O., Vigil, A. L., Fields, S., Queitsch, C., et al. (2021). The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution. Nat. Commun. 12 (1), 3334. doi: 10.1038/s41467-021-23675-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Elhanani, O., Ben-Uri, R., Keren, L. (2023). Spatial profiling technologies illuminate the tumor microenvironment. Cancer Cell 41 (3), 404–420. doi: 10.1016/j.ccell.2023.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Elosua-Bayes, M., Nieto, P., Mereu, E., Gut, I., Heyn, H. (2021). SPOTlight: seeded NMF regression to deconvolute spatial transcriptomics spots with single-cell transcriptomes. Nucleic Acids Res. 49 (9), e50. doi: 10.1093/nar/gkab043

PubMed Abstract | CrossRef Full Text | Google Scholar

Farmer, A., Thibivilliers, S., Ryu, K. H., Schiefelbein, J., 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 (3), 372–383. doi: 10.1016/j.molp.2021.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Gala, H. P., Lanctot, A., Jean-Baptiste, K., Guiziou, S., Chu, J. C., Zemke, J. E., et al. (2021). A single-cell view of the transcriptome during lateral root initiation in Arabidopsis thaliana. Plant Cell 33 (7), 2197–2220. doi: 10.1093/plcell/koab101

PubMed Abstract | CrossRef Full Text | Google Scholar

Giacomello, S., Lundeberg, J. (2018). Preparation of plant tissue to enable Spatial Transcriptomics profiling using barcoded microarrays. Nat. Protoc. 13 (11), 2425–2446. doi: 10.1038/s41596-018-0046-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Giacomello, S., Salmen, F., Terebieniec, B. K., Vickovic, S., Navarro, J. F., Alexeyenko, A., et al. (2017). Spatially resolved transcriptome profiling in model plant species. Nat. Plants 3, 17061. doi: 10.1038/nplants.2017.61

PubMed Abstract | CrossRef Full Text | Google Scholar

Giuliano, G., Pichersky, E., Malik, V. S., Timko, M. P., Scolnik, P. A., Cashmore, A. R. (1988). An evolutionarily conserved protein binding sequence upstream of a plant light-regulated gene. Proc. Natl. Acad. Sci. U.S.A. 85 (19), 7089–7093. doi: 10.1073/pnas.85.19.7089

PubMed Abstract | CrossRef Full Text | Google Scholar

Graeff, M., Rana, S., Wendrich, J. R., Dorier, J., Eekhout, T., Aliaga Fandino, A. C., et al. (2021). A single-cell morpho-transcriptomic map of brassinosteroid action in the Arabidopsis root. Mol. Plant 14 (12), 1985–1999. doi: 10.1016/j.molp.2021.07.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Grindberg, R. V., Yee-Greenbaum, J. L., McConnell, M. J., Novotny, M., O'Shaughnessy, A. L., Lambert, G. M., et al. (2013). RNA-sequencing from single nuclei. Proc. Natl. Acad. Sci. U.S.A. 110 (49), 19802–19807. doi: 10.1073/pnas.1319700110

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Liang, J., Lin, R., Zhang, L., Wu, J., Wang, X. (2021). Series-Spatial Transcriptome Profiling of Leafy Head Reveals the Key Transition Leaves for Head Formation in Chinese Cabbage. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.787826

CrossRef Full Text | Google Scholar

Gurazada, S. G. R., Cox, K. L., Czymmek, K. J., Meyers, B. C. (2021). Space: the final frontier - achieving single-cell, spatially resolved transcriptomics in plants. Emerg. Top. Life Sci. 5 (2), 179–188. doi: 10.1042/ETLS20200274

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., Chu, X., Yu, H., Ma, Y.-K., Wang, X.-J., Qian, W., et al. (2017). Single-cell transcriptome analysis reveals widespread monoallelic gene expression in individual rice mesophyll cells. Sci. Bull. 62 (19), 1304–1314. doi: 10.1016/j.scib.2017.09.011

CrossRef Full Text | Google Scholar

Hey, S., Baldauf, J., Opitz, N., Lithio, A., Pasha, A., Provart, N., et al. (2017). Complexity and specificity of the maize (Zea mays L.) root hair transcriptome. J. Exp. Bot. 68 (9), 2175–2185. doi: 10.1093/jxb/erx104

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, Z., Liu, Y., Zhang, M., Zhao, L., Jin, X., Liu, L., et al. (2021). High-throughput single-cell transcriptomics reveals the female germline differentiation trajectory in Arabidopsis thaliana. Commun. Biol. 4 (1), 1149. doi: 10.1038/s42003-021-02676-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, P., Liu, J., Zhao, J., Wilkins, B. J., Lupino, K., Wu, H., et al. (2018). Single-nucleus transcriptomic survey of cell diversity and functional maturation in postnatal mammalian hearts. Genes Dev. 32 (19-20), 1344–1357. doi: 10.1101/gad.316802.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, B., Lee, J. H., Bang, D. (2018). Single-cell RNA sequencing technologies and bioinformatics pipelines. Exp. Mol. Med. 50 (8), 1–14. doi: 10.1038/s12276-018-0071-8

PubMed Abstract | CrossRef Full Text | Google Scholar

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 Arabidopsis thaliana. Plant Cell 31 (5), 993–1011. doi: 10.1105/tpc.18.00785

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C., Sun, J., Li, R., Yan, S., Chen, W., Guo, L., et al. (2022). A reactive oxygen species burst causes haploid induction in maize. Mol. Plant 15 (6), 943–955. doi: 10.1016/j.molp.2022.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Keren-Shaul, H., Kenigsberg, E., Jaitin, D. A., David, E., Paul, F., Tanay, A., et al. (2019). MARS-seq2.0: an experimental and analytical pipeline for indexed sorting combined with single-cell RNA sequencing. Nat. Protoc. 14 (6), 1841–1862. doi: 10.1038/s41596-019-0164-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. Y., Symeonidi, E., Pang, T. Y., Denyer, T., Weidauer, D., Bezrutczyk, M., et al. (2021). Distinct identities of leaf phloem cells revealed by single cell transcriptomics. Plant Cell 33 (3), 511–530. doi: 10.1093/plcell/koaa060

PubMed Abstract | CrossRef Full Text | Google Scholar

Kivivirta, K., Herbert, D., Lange, M., Beuerlein, K., Altmuller, J., Becker, A. (2019). A protocol for laser microdissection (LMD) followed by transcriptome analysis of plant reproductive tissue in phylogenetically distant angiosperms. Plant Methods 15, 151. doi: 10.1186/s13007-019-0536-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Klein, A. M., Macosko, E. (2017). InDrops and Drop-seq technologies for single-cell sequencing. Lab. Chip 17 (15), 2540–2541. doi: 10.1039/c7lc90070h

PubMed Abstract | CrossRef Full Text | Google Scholar

Klein, A. M., Mazutis, L., Akartuna, I., Tallapragada, N., Veres, A., Li, V., et al. (2015). Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell 161 (5), 1187–1201. doi: 10.1016/j.cell.2015.04.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Krishnaswami, S. R., Grindberg, R. V., Novotny, M., Venepally, P., Lacar, B., Bhutani, K., et al. (2016). Using single nuclei for RNA-seq to capture the transcriptome of postmortem neurons. Nat. Protoc. 11 (3), 499–524. doi: 10.1038/nprot.2016.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Lake, B. B., Ai, R., Kaeser, G. E., Salathia, N. S., Yung, Y. C., Liu, R., et al. (2016). Neuronal subtypes and diversity revealed by single-nucleus RNA sequencing of the human brain. Science 352 (6293), 1586–1590. doi: 10.1126/science.aaf1204

PubMed Abstract | CrossRef Full Text | Google Scholar

Lareau, C. A., Ma, S., Duarte, F. M., Buenrostro, J. D. (2020). Inference and effects of barcode multiplets in droplet-based single-cell assays. Nat. Commun. 11 (1), 866. doi: 10.1038/s41467-020-14667-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lhamo, D., Luan, S. (2021). Potential Networks of Nitrogen-Phosphorus-Potassium Channels and Transporters in Arabidopsis Roots at a Single Cell Resolution. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.689545

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Wang, X., Wang, Y., Zhang, M., Hong, F., Wang, H., et al. (2022). Cross-species single-cell transcriptomic analysis reveals divergence of cell composition and functions in mammalian ileum epithelium. Cell Regener. 11 (1), 19. doi: 10.1186/s13619-022-00118-7

CrossRef Full Text | Google Scholar

Li, X., Li, L., Yan, J. (2015). Dissecting meiotic recombination based on tetrad analysis by single-microspore sequencing in maize. Nat. Commun. 6, 6648. doi: 10.1038/ncomms7648

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y. C., Li, W., Chen, H., Li, Q., Sun, Y. H., Shi, R., et al. (2014). A simple improved-throughput xylem protoplast system for studying wood formation. Nat. Protoc. 9 (9), 2194–2205. doi: 10.1038/nprot.2014.147

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Leng, J., Li, Y., Ge, T., Li, J., Chen, Y., et al. (2022a). A spatiotemporal atlas of organogenesis in the development of orchid flowers. Nucleic Acids Res. 50 (17), 9724–9737. doi: 10.1093/nar/gkac773

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Hu, D., Du, P., Wang, L., Liang, X., Li, H., et al. (2021a). Single-cell RNA-seq describes the transcriptome landscape and identifies critical transcription factors in the leaf blade of the allotetraploid peanut (Arachis hypogaea L.). Plant Biotechnol. J. 19 (11), 2261–2276. doi: 10.1111/pbi.13656

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Li, P., Cheng, S., Zhao, Z., Liu, Y., Wei, Y., et al. (2022b). Protoplast dissociation and transcriptome analysis provides insights to salt stress response in cotton. Int. J. Mol. Sci. 23 (5), 2845. doi: 10.3390/ijms23052845

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Liang, Z., Feng, D., Jiang, S., Wang, Y., Du, Z., et al. (2021b). Transcriptional landscape of rice roots at the single-cell resolution. Mol. Plant 14 (3), 384–394. doi: 10.1016/j.molp.2020.12.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., Zhang, Y., Fang, X., Tran, S., Zhai, N., Yang, Z., et al. (2022c). Transcriptional landscapes of de novo root regeneration from detached Arabidopsis leaves revealed by time-lapse and single-cell RNA sequencing analyses. Plant Commun. 3 (4), 100306. doi: 10.1016/j.xplc.2022.100306

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Wang, J., Zhou, Y., Zhang, Y., Qin, A., Yu, X., et al. (2022d). Identification of novel regulators required for early development of vein pattern in the cotyledons by single-cell RNA-sequencing. Plant J. 110 (1), 7–22. doi: 10.1111/tpj.15719

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Zhou, Y., Guo, J., Li, J., Tian, Z., Zhu, Z., et al. (2020). Global Dynamic Molecular Profiling of Stomatal Lineage Cell Development by Single-Cell RNA Sequencing. Mol. Plant 13 (8), 1178–1193. doi: 10.1016/j.molp.2020.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

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), 66. doi: 10.1186/s13059-021-02288-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Longo, S. K., Guo, M. G., Ji, A. L., Khavari, P. A. (2021). Integrating single-cell and spatial transcriptomics to elucidate intercellular tissue dynamics. Nat. Rev. Genet. 22 (10), 627–644. doi: 10.1038/s41576-021-00370-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez-Anido, C. B., Vaten, A., Smoot, N. K., Sharma, N., Guo, V., Gong, Y., et al. (2021). Single-cell resolution of lineage trajectories in the Arabidopsis stomatal lineage and developing leaf. Dev. Cell 56 (7), 1043–1055.e1044. doi: 10.1016/j.devcel.2021.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, L., Tian, S., Liao, H., Zhang, J., Yang, X., Labavitch, J. M., et al. (2013). Analysis of metal element distributions in rice (Oryza sativa L.) seeds and relocation during germination based on X-ray fluorescence imaging of Zn, Fe, K, Ca, and Mn. PloS One 8 (2), e57360. doi: 10.1371/journal.pone.0057360

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, L., Xie, R., Liu, T., Wang, H., Hou, D., Du, Y., et al. (2017). Spatial imaging and speciation of Cu in rice (Oryza sativa L.) roots using synchrotron-based X-ray microfluorescence and X-ray absorption spectroscopy. Chemosphere 175, 356–364. doi: 10.1016/j.chemosphere.2017.02.082

PubMed Abstract | CrossRef Full Text | Google Scholar

Macosko, E. Z., Basu, A., Satija, R., Nemesh, J., Shekhar, K., Goldman, M., et al. (2015). Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell 161 (5), 1202–1214. doi: 10.1016/j.cell.2015.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Maitra, M., Nagy, C., Chawla, A., Wang, Y. C., Nascimento, C., Suderman, M., et al. (2021). Extraction of nuclei from archived postmortem tissues for single-nucleus sequencing applications. Nat. Protoc. 16 (6), 2788–2801. doi: 10.1038/s41596-021-00514-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Marand, A. P., Chen, Z., Gallavotti, A., Schmitz, R. J. (2021). A cis-regulatory atlas in maize at single-cell resolution. Cell 184 (11), 3041–3055.e3021. doi: 10.1016/j.cell.2021.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Misra, C. S., Santos, M. R., Rafael-Fernandes, M., Martins, N. P., Monteiro, M., Becker, J. D. (2019). Transcriptomics of Arabidopsis sperm cells at single-cell resolution. Plant Reprod. 32 (1), 29–38. doi: 10.1007/s00497-018-00355-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Mo, Y., Jiao, Y. (2022). Advances and applications of single-cell omics technologies in plant research. Plant J. 110 (6), 1551–1563. doi: 10.1111/tpj.15772

PubMed Abstract | CrossRef Full Text | Google Scholar

Moffitt, J. R., Bambah-Mukku, D., Eichhorn, S. W., Vaughn, E., Shekhar, K., Perez, J. D., et al. (2018). Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region. Science 362 (6416), eaau5324. doi: 10.1126/science.aau5324

PubMed Abstract | CrossRef Full Text | Google Scholar

Moreno-Romero, J., Santos-Gonzalez, J., Hennig, L., Kohler, C. (2017). Applying the INTACT method to purify endosperm nuclei and to generate parental-specific epigenome profiles. Nat. Protoc. 12 (2), 238–254. doi: 10.1038/nprot.2016.167

PubMed Abstract | CrossRef Full Text | Google Scholar

Moreno-Villena, J. J., Zhou, H., Gilman, I. S., Tausta, S. L., Cheung, C. Y. M., Edwards, E. J. (2022). Spatial resolution of an integrated C4+CAM photosynthetic metabolism. Sci. Adv. 8 (31), eabn2349. doi: 10.1126/sciadv.abn2349

PubMed Abstract | CrossRef Full Text | Google Scholar

Nadelmann, E. R., Gorham, J. M., Reichart, D., Delaughter, D. M., Wakimoto, H., Lindberg, E. L., et al. (2021). Isolation of nuclei from mammalian cells and tissues for single-nucleus molecular profiling. Curr. Protoc. 1 (5), e132. doi: 10.1002/cpz1.132

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelms, B., Walbot, V. (2019). Defining the developmental program leading to meiosis in maize. Science 364 (6435), 52–56. doi: 10.1126/science.aav6428

PubMed Abstract | CrossRef Full Text | Google Scholar

Neumann, M., Xu, X., Smaczniak, C., Schumacher, J., Yan, W., Bluthgen, N., et al. (2022). A 3D gene expression atlas of the floral meristem based on spatial reconstruction of single nucleus RNA sequencing data. Nat. Commun. 13 (1), 2838. doi: 10.1038/s41467-022-30177-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Ortiz-Ramirez, C., Arevalo, E. D., Xu, X., Jackson, D. P., Birnbaum, K. D. (2018). An Efficient Cell Sorting Protocol for Maize Protoplasts. Curr. Protoc. Plant Biol. 3 (3), e20072. doi: 10.1002/cppb.20072

PubMed Abstract | CrossRef Full Text | Google Scholar

Potter, S. S. (2018). Single-cell RNA sequencing for the study of development, physiology and disease. Nat. Rev. Nephrol. 14 (8), 479–492. doi: 10.1038/s41581-018-0021-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Qian, X., Harris, K. D., Hauling, T., Nicoloutsopoulos, D., Munoz-Manchado, A. B., Skene, N., et al. (2020). Probabilistic cell typing enables fine mapping of closely related cell types in situ. Nat. Methods 17 (1), 101–106. doi: 10.1038/s41592-019-0631-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramsköld, D., Luo, S., Wang, Y.-C., Li, R., Deng, Q., Faridani, O. R., et al. (2012). Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells. Nat. Biotechnol. 30 (8), 777–782. doi: 10.1038/nbt.2282

PubMed Abstract | CrossRef Full Text | Google Scholar

Rao, A., Barkley, D., Franca, G. S., Yanai, I. (2021). Exploring tissue architecture using spatial transcriptomics. Nature 596 (7871), 211–220. doi: 10.1038/s41586-021-03634-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Roszak, P., Heo, J. O., Blob, B., Toyokura, K., Sugiyama, Y., de Luis Balaguer, M. A., et al. (2021). Cell-by-cell dissection of phloem development links a maturation gradient to cell specialization. Science 374 (6575), eaba5531. doi: 10.1126/science.aba5531

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryu, K. H., Huang, L., Kang, H. M., Schiefelbein, J. (2019). Single-Cell RNA Sequencing Resolves Molecular Relationships Among Individual Plant Cells. Plant Physiol. 179 (4), 1444–1456. doi: 10.1104/pp.18.01482

PubMed Abstract | CrossRef Full Text | Google Scholar

Salguero-Linares, J., Serrano, I., Ruiz-Solani, N., Salas-Gomez, M., Phukan, U. J., Gonzalez, V. M., et al. (2022). Robust transcriptional indicators of immune cell death revealed by spatiotemporal transcriptome analyses. Mol. Plant 15 (6), 1059–1075. doi: 10.1016/j.molp.2022.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Satija, R., Farrell, J. A., Gennert, D., Schier, A. F., Regev, A. (2015). Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 33 (5), 495–502. doi: 10.1038/nbt.3192

PubMed Abstract | CrossRef Full Text | Google Scholar

Satterlee, J. W., Strable, J., Scanlon, M. J. (2020). Plant stem-cell organization and differentiation at single-cell resolution. Proc. Natl. Acad. Sci. U.S.A. 117 (52), 33689–33699. doi: 10.1073/pnas.2018788117

PubMed Abstract | CrossRef Full Text | Google Scholar

See, P., Lum, J., Chen, J., Ginhoux, F. (2018). A Single-Cell Sequencing Guide for Immunologists. Front. Immunol. 9. doi: 10.3389/fimmu.2018.02425

PubMed Abstract | CrossRef Full Text | Google Scholar

Shahan, R., Hsu, C. W., Nolan, T. M., Cole, B. J., Taylor, I. W., Greenstreet, L., et al. (2022). A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants. Dev. Cell 57 (4), 543–560.e549. doi: 10.1016/j.devcel.2022.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Shalek, A. K., Satija, R., Adiconis, X., Gertner, R. S., Gaublomme, J. T., Raychowdhury, R., et al. (2013). Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells. Nature 498 (7453), 236–240. doi: 10.1038/nature12172

PubMed Abstract | CrossRef Full Text | Google Scholar

Shaw, R., Tian, X., Xu, J. (2021). Single-cell transcriptome analysis in plants: Advances and challenges. Mol. Plant 14 (1), 115–126. doi: 10.1016/j.molp.2020.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Shojaee, A., Saavedra, M., Huang, S. C. (2021). Potentials of single-cell genomics in deciphering cellular phenotypes. Curr. Opin. Plant Biol. 63, 102059. doi: 10.1016/j.pbi.2021.102059

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (7), 2241–2247.e2244. doi: 10.1016/j.celrep.2019.04.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Sikorskaite, S., Rajamaki, M. L., Baniulis, D., Stanys, V., Valkonen, J. P. (2013). Protocol: Optimised methodology for isolation of nuclei from leaves of species in the Solanaceae and Rosaceae families. Plant Methods 9, 31. doi: 10.1186/1746-4811-9-31

PubMed Abstract | CrossRef Full Text | Google Scholar

Somssich, M., Khan, G. A., Persson, S. (2016). Cell wall heterogeneity in root development of arabidopsis. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.01242

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Q., Ando, A., Jiang, N., Ikeda, Y., 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), 178. doi: 10.1186/s13059-020-02094-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Stahl, P. L., Salmen, F., Vickovic, S., Lundmark, A., Navarro, J. F., Magnusson, J., et al. (2016). Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science 353 (6294), 78–82. doi: 10.1126/science.aaf2403

PubMed Abstract | CrossRef Full Text | Google Scholar

Stuart, T., Butler, A., Hoffman, P., Hafemeister, C., Papalexi, E., Mauck, W. M., 3rd, et al. (2019). Comprehensive integration of single-cell data. Cell 177 (7), 1888–1902.e1821. doi: 10.1016/j.cell.2019.05.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Sunaga-Franze, D. Y., Muino, J. M., Braeuning, C., Xu, X., Zong, M., Smaczniak, C., et al. (2021). Single-nucleus RNA sequencing of plant tissues using a nanowell-based system. Plant J. 108 (3), 859–869. doi: 10.1111/tpj.15458

PubMed Abstract | CrossRef Full Text | Google Scholar

Svensson, V., Vento-Tormo, R., Teichmann, S. A. (2018). Exponential scaling of single-cell RNA-seq in the past decade. Nat. Protoc. 13 (4), 599–604. doi: 10.1038/nprot.2017.149

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, F., Barbacioru, C., Wang, Y., Nordman, E., Lee, C., Xu, N., et al. (2009). mRNA-Seq whole-transcriptome analysis of a single cell. Nat. Methods 6 (5), 377–382. doi: 10.1038/nmeth.1315

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Z., Fan, Y., Zhang, L., Zheng, C., Chen, A., Sun, Y., et al. (2021). Quantitative metabolome and transcriptome analysis reveals complex regulatory pathway underlying photoinduced fiber color formation in cotton. Gene 767, 145180. doi: 10.1016/j.gene.2020.145180

PubMed Abstract | CrossRef Full Text | Google Scholar

Tenorio Berrio, R., Verstaen, K., Vandamme, N., Pevernagie, J., Achon, I., Van Duyse, J., et al. (2022). Single-cell transcriptomics sheds light on the identity and metabolism of developing leaf cells. Plant Physiol. 188 (2), 898–918. doi: 10.1093/plphys/kiab489

PubMed Abstract | CrossRef Full Text | Google Scholar

Thibivilliers, S., Anderson, D., Libault, M. (2020). Isolation of Plant Root Nuclei for Single Cell RNA Sequencing. Curr. Protoc. Plant Biol. 5 (4), e20120. doi: 10.1002/cppb.20120

PubMed Abstract | CrossRef Full Text | Google Scholar

Tirichine, L., Andrey, P., Biot, E., Maurin, Y., Gaudin, V. (2009). 3D fluorescent in situ hybridization using Arabidopsis leaf cryosections and isolated nuclei. Plant Methods 5, 11. doi: 10.1186/1746-4811-5-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Torii, K., Kubota, A., Araki, T., Endo, M. (2020). Time-Series Single-Cell RNA-Seq Data Reveal Auxin Fluctuation during Endocycle. Plant Cell Physiol. 61 (2), 243–254. doi: 10.1093/pcp/pcz228

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Brink, S. C., Sage, F., Vertesy, A., Spanjaard, B., Peterson-Maduro, J., Baron, C. S., et al. (2017). Single-cell sequencing reveals dissociation-induced gene expression in tissue subpopulations. Nat. Methods 14 (10), 935–936. doi: 10.1038/nmeth.4437

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Li, Y., Wu, T., Miao, C., Xie, M., Ding, B., et al. (2021a). Single-cell-type transcriptomic analysis reveals distinct gene expression profiles in wheat guard cells in response to abscisic acid. Funct. Plant Biol. 48 (11), 1087–1099. doi: 10.1071/FP20368

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Wu, Y., Peng, A., Cui, J., Zhao, M., Pan, Y., et al. (2022). Single-cell transcriptome atlas reveals developmental trajectories and a novel metabolic pathway of catechin esters in tea leaves. Plant Biotechnol. J. 20, 2089–2106. doi: 10.1111/pbi.13891

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Allen, W. E., Wright, M. A., Sylwestrak, E. L., Samusik, N., Vesuna, S., et al. (2018). Three-dimensional intact-tissue sequencing of single-cell transcriptional states. Science 361 (6400), eaat5691. doi: 10.1126/science.aat5691

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., He, Y., Zhang, Q., Ren, X., Zhang, Z. (2021b). Direct Comparative Analyses of 10X Genomics Chromium and Smart-seq2. Genomics Proteomics Bioinf. 19 (2), 253–266. doi: 10.1016/j.gpb.2020.02.005

CrossRef Full Text | Google Scholar

Wang, Y., Huan, Q., Li, K., Qian, W. (2021c). Single-cell transcriptome atlas of the leaf and root of rice seedlings. J. Genet. Genomics 48 (10), 881–898. doi: 10.1016/j.jgg.2021.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, R., He, S., Bai, S., Sei, E., Hu, M., Thompson, A., et al. (2022). Spatial charting of single-cell transcriptomes in tissues. Nat. Biotechnol. 40 (8), 1190–1199. doi: 10.1038/s41587-022-01233-1

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (6518), eaay4970. doi: 10.1126/science.aay4970

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, H., Kirita, Y., Donnelly, E. L., Humphreys, B. D. (2019). Advantages of Single-Nucleus over Single-Cell RNA Sequencing of Adult Kidney: Rare Cell Types and Novel Cell States Revealed in Fibrosis. J. Am. Soc. Nephrol. 30 (1), 23–32. doi: 10.1681/ASN.2018090912

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, K., Sun, H. X., Li, J., Li, J., Zhao, Y., Chen, L., et al. (2022). The single-cell stereo-seq reveals region-specific cell subtypes and transcriptome profiling in Arabidopsis leaves. Dev. Cell 57 (10), 1299–1310.e1294. doi: 10.1016/j.devcel.2022.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Y., Jiang, S., Li, L., Yu, X., Wang, Y., Luo, C., et al. (2020). Single-Cell RNA Sequencing Efficiently Predicts Transcription Factor Targets in Plants. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.603302

CrossRef Full Text | Google Scholar

Xu, X., Crow, M., Rice, B. R., Li, F., Harris, B., Liu, L., et al. (2021a). Single-cell RNA sequencing of developing maize ears facilitates functional analysis and trait candidate gene discovery. Dev. Cell 56 (4), 557–568.e556. doi: 10.1016/j.devcel.2020.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X. F., Zhu, H. Y., Ren, Y. F., Feng, C., Ye, Z. H., Cai, H. M., et al. (2021b). Efficient isolation and purification of tissue-specific protoplasts from tea plants (Camellia sinensis (L.) O. Kuntze). Plant Methods 17 (1), 84. doi: 10.1186/s13007-021-00783-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, B., Minne, M., Brunoni, F., Plackova, L., Petrik, I., Sun, Y., et al. (2021). Non-cell autonomous and spatiotemporal signalling from a tissue organizer orchestrates root vascular development. Nat. Plants 7 (11), 1485–1494. doi: 10.1038/s41477-021-01017-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoo, S. D., Cho, Y. H., Sheen, J. (2007). Arabidopsis mesophyll protoplasts: a versatile cell system for transient gene expression analysis. Nat. Protoc. 2 (7), 1565–1572. doi: 10.1038/nprot.2007.199

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Yang, H., Hu, D., Li, B., Lin, Y., Yao, W., et al. (2021a). Single-cell RNA sequencing of meiocytes and microspores reveals the involvement of the Rf4 gene in redox homeostasis of CMS-C maize. Crop J. 9 (6), 1237–1247. doi: 10.1016/j.cj.2021.06.012

CrossRef Full Text | Google Scholar

Zhang, M., Zou, Y., Xu, X., Zhang, X., Gao, M., Song, J., et al. (2020). Highly parallel and efficient single cell mRNA sequencing with paired picoliter chambers. Nat. Commun. 11 (1), 2118. doi: 10.1038/s41467-020-15765-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Xu, M., Qiu, Z., Wang, K., Du, Y., Gu, L., et al. (2016). Spatiotemporal transcriptome provides insights into early fruit development of tomato (Solanum lycopersicum). Sci. Rep. 6, 23173. doi: 10.1038/srep23173

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, T. Q., Chen, Y., Liu, Y., Lin, W. H., Wang, J. W. (2021b). Single-cell transcriptome atlas and chromatin accessibility landscape reveal differentiation trajectories in the rice root. Nat. Commun. 12 (1), 2053. doi: 10.1038/s41467-021-22352-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, T. Q., Chen, Y., Wang, J. W. (2021c). A single-cell analysis of the Arabidopsis vegetative shoot apex. Dev. Cell 56 (7), 1056–1074.e1058. doi: 10.1016/j.devcel.2021.02.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, T. Q., Xu, Z. G., Shang, G. D., Wang, J. W. (2019). A Single-Cell RNA Sequencing Profiles the Developmental Landscape of Arabidopsis Root. Mol. Plant 12 (5), 648–660. doi: 10.1016/j.molp.2019.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, G. X., Terry, J. M., Belgrader, P., Ryvkin, P., Bent, Z. W., Wilson, R., et al. (2017). Massively parallel digital transcriptional profiling of single cells. Nat. Commun. 8, 14049. doi: 10.1038/ncomms14049

PubMed Abstract | CrossRef Full Text | Google Scholar

Zong, J., Wang, L., Zhu, L., Bian, L., Zhang, B., Chen, X., et al. (2022). A rice single cell transcriptomic atlas defines the developmental trajectories of rice floret and inflorescence meristems. New Phytol. 234 (2), 494–512. doi: 10.1111/nph.18008

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: single-cell transcriptomics, spatial transcriptomics, single-nucleus RNA-seq, single-cell RNA-seq, protoplast

Citation: Chen C, Ge Y and Lu L (2023) Opportunities and challenges in the application of single-cell and spatial transcriptomics in plants. Front. Plant Sci. 14:1185377. doi: 10.3389/fpls.2023.1185377

Received: 13 March 2023; Accepted: 26 July 2023;
Published: 11 August 2023.

Edited by:

Kashif Ali, Shaheed Zulfiqar Ali Bhutto Institute of Science and Technology, Pakistan

Reviewed by:

Congting Ye, Xiamen University, China
Tri Rini Nuringtyas, Gadjah Mada University, Indonesia

Copyright © 2023 Chen, Ge and Lu. 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: Lingli Lu, bHVsaW5nbGlAemp1LmVkdS5jbg==

Disclaimer: 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.