# Best practices for the execution, analysis, and data storage of plant single-cell/nucleus transcriptomics

Carolin Grones, Thomas Eekhout, Dongbo Shi, Manuel Neumann, Lea S Berg, Yuji Ke, Rachel Shahan, Kevin L Cox, Fabio Gomez-Cano, Hilde Nelissen, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Carolin Grones, Thomas Eekhout, Dongbo Shi, Manuel Neumann, Lea S Berg, et al. Best practices for the execution, analysis, and data storage of plant single-cell/nucleus transcriptomics. The Plant cell, 2024, 36 (4), pp.812-828. 10.1093/plcell/koae003. hal-04600040

HAL Id: hal-04600040

https://hal.inrae.fr/hal-04600040

Submitted on 4 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_1.png" /></figure>

# Best practices for the execution, analysis, and data storage of plant single-cell/nucleus transcriptomics

Carolin Grones,¹,² Thomas Eekhout ID,¹,²,³ Dongbo Shi,⁴,⁵ Manuel Neumann ID,⁶ Lea S. Berg ID,⁷ Yuji Ke,¹,²

Rachel Shahan ID,⁸,⁹ Kevin L. Cox Jr ID,¹⁰ Fabio Gomez-Cano,¹¹ Hilde Nelissen ID,¹,² Jan U. Lohmann ID,⁴

Stefania Giacomello,¹² Olivier C. Martin ID,¹³ Benjamin Cole ID,¹⁴ Jia-Wei Wang ID,¹⁵

Kerstin Kaufmann ID,⁶ Michael T. Raissig ID,⁷ Gergo Palfalvi,¹⁶ Thomas Greb ID,⁴ Marc Libault ID¹⁷,*

and Bert De Rybel ID¹,²,*

1 Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent 9052, Belgium

2 VIB Centre for Plant Systems Biology, Ghent 9052, Belgium

3 VIB Single Cell Core Facility, Ghent 9052, Belgium

4 Centre for Organismal Studies, Heidelberg University, 69120 Heidelberg, Germany

5 Institute of Biochemistry and Biology, University of Potsdam, 14476 Potsdam, Germany

6 Institute of Biology, Humboldt-Universität zu Berlin, 10115 Berlin, Germany

7 Institute of Plant Sciences, University of Bern, 3012 Bern, Switzerland

8 Department of Biology, Duke University, Durham, NC 27708, USA

9 Howard Hughes Medical Institute, Duke University, Durham, NC 27708, USA

10 Donald Danforth Plant Science Center, St. Louis, MO 63132, USA

11 Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI 48109, USA

12 SciLifeLab, Department of Gene Technology, KTH Royal Institute of Technology, 17165 Solna, Sweden

13 Universities of Paris-Saclay, Paris-Cité and Evry, CNRS, INRAE, Institute of Plant Sciences Paris-Saclay, Gif-sur-Yvette 91192, France

14 DOE-Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA

15 National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences (CEMPS), Institute of Plant Physiology and Ecology (SIPPE), Chinese Academy of Sciences (CAS), Shanghai 200032, China

16 Department of Comparative Development and Genetics, Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany

17 Division of Plant Science and Technology, Interdisciplinary Plant Group, College of Agriculture, Food, and Natural Resources, University of Missouri-Columbia, Columbia, MO 65201, USA

## Abstract

Single-cell and single-nucleus RNA-sequencing technologies capture the expression of plant genes at an unprecedented resolution. Therefore, these technologies are gaining traction in plant molecular and developmental biology for elucidating the transcriptional changes across cell types in a specific tissue or organ, upon treatments, in response to biotic and abiotic stresses, or between genotypes. Despite the rapidly accelerating use of these technologies, collective and standardized experimental and analytical procedures to support the acquisition of high-quality data sets are still missing. In this commentary, we discuss common challenges associated with the use of single-cell transcriptomics in plants and propose general guidelines to improve reproducibility, quality, comparability, and interpretation and to make the data readily available to the community in this fast-developing field of research.

<footnote>*Author for correspondence: libaultm@missouri.edu (M.T), bert.derybel@psb.vib-ugent.be (B.D.R.)</footnote>

<footnote>The authors responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (https://academic.oup.com/plcell/pages/General-Instructions) are: Marc Libault (libaultm@missouri.edu) and Bert De Rybel (bert.derybel@psb.vib-ugent.be)</footnote>
<figure><img src="image_3.png" /></figure>

# Best practices for the execution, analysis, and data storage of plant single-cell/nucleus transcriptomics

Carolin Grones,¹,² Thomas Eekhout,¹,²,³ Dongbo Shi,⁴,⁵ Manuel Neumann,⁶ Lea S. Berg,⁷ Yuji Ke,¹,²

Rachel Shahan,⁸,⁹ Kevin L. Cox Jr,¹⁰ Fabio Gomez-Cano,¹¹ Hilde Nelissen,¹,² Jan U. Lohmann,⁴

Stefania Giacomello,¹² Olivier C. Martin,¹³ Benjamin Cole,¹⁴ Jia-Wei Wang,¹⁵

Kerstin Kaufmann,⁶ Michael T. Raissig,⁷ Gergo Palfalvi,¹⁶ Thomas Greb,⁴ Marc Libault,¹⁷,*

and Bert De Rybel,¹,²,*

1 Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent 9052, Belgium

2 VIB Centre for Plant Systems Biology, Ghent 9052, Belgium

3 VIB Single Cell Core Facility, Ghent 9052, Belgium

4 Centre for Organismal Studies, Heidelberg University, 69120 Heidelberg, Germany

5 Institute of Biochemistry and Biology, University of Potsdam, 14476 Potsdam, Germany

6 Institute of Biology, Humboldt-Universität zu Berlin, 10115 Berlin, Germany

7 Institute of Plant Sciences, University of Bern, 3012 Bern, Switzerland

8 Department of Biology, Duke University, Durham, NC 27708, USA

9 Howard Hughes Medical Institute, Duke University, Durham, NC 27708, USA

10 Donald Danforth Plant Science Center, St. Louis, MO 63132, USA

11 Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI 48109, USA

12 SciLifeLab, Department of Gene Technology, KTH Royal Institute of Technology, 17165 Solna, Sweden

13 Universities of Paris-Saclay, Paris-Cité and Evry, CNRS, INRAE, Institute of Plant Sciences Paris-Saclay, Gif-sur-Yvette 91192, France

14 DOE-Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA

15 National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences (CEMPS), Institute of Plant Physiology and Ecology (SIPPE), Chinese Academy of Sciences (CAS), Shanghai 200032, China

16 Department of Comparative Development and Genetics, Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany

17 Division of Plant Science and Technology, Interdisciplinary Plant Group, College of Agriculture, Food, and Natural Resources, University of Missouri-Columbia, Columbia, MO 65201, USA

## Abstract

Single-cell and single-nucleus RNA-sequencing technologies capture the expression of plant genes at an unprecedented resolution. Therefore, these technologies are gaining traction in plant molecular and developmental biology for elucidating the transcriptional changes across cell types in a specific tissue or organ, upon treatments, in response to biotic and abiotic stresses, or between genotypes. Despite the rapidly accelerating use of these technologies, collective and standardized experimental and analytical procedures to support the acquisition of high-quality data sets are still missing. In this commentary, we discuss common challenges associated with the use of single-cell transcriptomics in plants and propose general guidelines to improve reproducibility, quality, comparability, and interpretation and to make the data readily available to the community in this fast-developing field of research.

<footnote>*Author for correspondence: libaultm@missouri.edu (M.T), bert.derybel@psb.vib-ugent.be (B.D.R.)</footnote>

<footnote>The authors responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (https://academic.oup.com/plcell/pages/General-Instructions) are: Marc Libault (libaultm@missouri.edu) and Bert De Rybel (bert.derybel@psb.vib-ugent.be)</footnote>
# Introduction: plant-specific challenges for single-cell approaches

Plant molecular and developmental biologists are fully embracing single-cell applications. Specifically, single-cell RNA-sequencing (scRNA-seq) and single-nucleus RNA-sequencing (snRNA-seq) are gaining a lot of traction while spatial transcriptomics is emerging as a promising complementary technology (Fig. 1). Despite an increase in the use and publication of plant single-cell experimentation (Fig. 1A), it is fair to say that the plant field has, so far, not settled on common strategies, protocols, or analysis methods. Given the high complexity of the different technologies and sample types (Fig. 1, B and C), we feel it is important to provide a best-practice workflow and guidelines that will help in establishing a collectively accepted quality cutoff. These guidelines will aid in the evaluation of experimental approaches and computational analyses of single-cell transcriptomic data, while also offering solutions to commonly observed challenges, thereby improving the reproducibility and comparability of experiments in the broader field of plant research. The present coauthors collectively accept these guidelines and commit to applying them to their research. We also highlight examples where consensus has not yet been achieved between coauthors, which will need to be resolved when both the technologies and the field develop further. As one example, single-cell multiomics and spatial transcriptomics are, in our opinion, not established enough in the plant field to propose any sort of definitive rules at this moment in time.

To date, we have identified the following 8 main challenges in the field of plant single-cell/nucleus (sc/sn) transcriptomics: (i) deciding on the best single-cell methods to answer a specific biological question; (ii) understanding experimental variability; (iii) biases in protocols and platforms; (iv) deciding on a sequencing strategy; (v) generating expression matrices and defining high-quality cells; (vi) constructing cell clusters and mapping them to cell types; (vii) trajectory inference methods and applications; and (viii) documenting and publishing data sets (Fig. 2). Each of these challenges is discussed in detail in the following sections.

## Challenge 1: selecting the best approach to answer a specific biological question

Before considering the best experimental approach to obtain single-cell transcriptomic data, it is important to evaluate the potential benefits of accessing single-cell resolution over bulk RNA-sequencing (bulk RNA-seq). This assessment depends on the biological system considered and the biological questions to answer. sc/snRNA-seq provides a snapshot of the transcriptome of each cell within an organism, offering a high spatiotemporal resolution of the dynamic gene regulation involved in plant development, cell differentiation, or responses to biotic and abiotic stresses. Single-cell transcriptomics can also offer the required resolution to study

cell-type-specific responses during cellular evolution and adaptation mechanisms among plant species (Guillotin et al. 2023). We thus advocate using sc/sn transcriptomic technologies over bulk RNA-seq when working with a celluarly complex sample or to capture dynamic transcriptomic responses to stimuli. In other cases, the question at hand might be more easily addressed using bulk RNA-seq or targeted gene expression experiments.

### scRNA-seq versus snRNA-seq

A thorough understanding of the different strategies and types of protocols (Fig. 1C) is essential before one can make an educated decision on which technology will best answer a specific biological question. A first example of an important choice is whether to profile the transcriptome from isolated nuclei or cells. When doing scRNA-seq, the most popular choice to obtain single cells from a plant organ requires the enzymatic digestion of cell walls and the generation of so-called protoplasts. There are a number of disadvantages of using protoplasts such as some tissues (e.g. sclerenchyma) and species (e.g. sorghum; Guillotin et al. 2023) are recalcitrant to cell wall digestion; enzymatic digestion affects the transcriptional status of the plant cells and could bias the outcome of experiments (Birnbaum et al. 2005); and the large size of protoplasts reduces their capture efficiency with most of the currently available commercialized single-cell platforms. Nucleus isolation followed by snRNA-seq gained traction in plant single-cell transcriptomics as well. However, the recovered data content per nucleus (e.g. Unique Molecular Identifier (UMI) or genes) is up to 10 (for UMIs) and 3 times (for genes) lower compared to scRNA-seq (Farmer et al. 2021; Guillotin et al. 2023). Furthermore, even though the transcriptome coverage is similarly efficient between scRNA-seq and snRNA-seq (e.g. 89% of all *Arabidopsis* expressed genes were captured in snRNA-seq data; Farmer et al. 2021; Guillotin et al. 2023), a nuclear transcriptome and a cellular transcriptome are not equivalent (Lee and Bailey-Serres 2019; Reynoso et al. 2019). For example, differences in abundance and composition between transcripts obtained from nuclear versus polyA RNA under hypoxia point toward nuclear transcript retention or enrichment as part of the stress response (Lee and Bailey-Serres 2019; Reynoso et al. 2019). Furthermore, the half-life of the transcripts (estimations range between 12 min and more than 24 h in *Arabidopsis* cells; Narsai et al. 2007) suggests that the cellular transcriptome is the result of the accumulation of the transcript synthesis over time, while the nuclear transcriptome is considered to accommodate faster to changes in gene activity. These differences are important to consider when selecting and later interpreting a single-cellular transcriptome versus a nuclear-based transcriptome and should be determined by the experimental system or biological question. Therefore, when studying, e.g. early stress responses of plant cells, a snRNA-seq could achieve higher resolution of rapid transcriptomic changes, while scRNA-seq might be more informative when understanding the biology
# Introduction: plant-specific challenges for single-cell approaches

Plant molecular and developmental biologists are fully embracing single-cell applications. Specifically, single-cell RNA-sequencing (scRNA-seq) and single-nucleus RNA-sequencing (snRNA-seq) are gaining a lot of traction while spatial transcriptomics is emerging as a promising complementary technology (Fig. 1). Despite an increase in the use and publication of plant single-cell experimentation (Fig. 1A), it is fair to say that the plant field has, so far, not settled on common strategies, protocols, or analysis methods. Given the high complexity of the different technologies and sample types (Fig. 1, B and C), we feel it is important to provide a best-practice workflow and guidelines that will help in establishing a collectively accepted quality cutoff. These guidelines will aid in the evaluation of experimental approaches and computational analyses of single-cell transcriptomic data, while also offering solutions to commonly observed challenges, thereby improving the reproducibility and comparability of experiments in the broader field of plant research. The present coauthors collectively accept these guidelines and commit to applying them to their research. We also highlight examples where consensus has not yet been achieved between coauthors, which will need to be resolved when both the technologies and the field develop further. As one example, single-cell multiomics and spatial transcriptomics are, in our opinion, not established enough in the plant field to propose any sort of definitive rules at this moment in time.

To date, we have identified the following 8 main challenges in the field of plant single-cell/nucleus (sc/sn) transcriptomics: (i) deciding on the best single-cell methods to answer a specific biological question; (ii) understanding experimental variability; (iii) biases in protocols and platforms; (iv) deciding on a sequencing strategy; (v) generating expression matrices and defining high-quality cells; (vi) constructing cell clusters and mapping them to cell types; (vii) trajectory inference methods and applications; and (viii) documenting and publishing data sets (Fig. 2). Each of these challenges is discussed in detail in the following sections.

## Challenge 1: selecting the best approach to answer a specific biological question

Before considering the best experimental approach to obtain single-cell transcriptomic data, it is important to evaluate the potential benefits of accessing single-cell resolution over bulk RNA-sequencing (bulk RNA-seq). This assessment depends on the biological system considered and the biological questions to answer. sc/snRNA-seq provides a snapshot of the transcriptome of each cell within an organism, offering a high spatiotemporal resolution of the dynamic gene regulation involved in plant development, cell differentiation, or responses to biotic and abiotic stresses. Single-cell transcriptomics can also offer the required resolution to study

cell-type-specific responses during cellular evolution and adaptation mechanisms among plant species (Guillotin et al. 2023). We thus advocate using sc/sn transcriptomic technologies over bulk RNA-seq when working with a celluarly complex sample or to capture dynamic transcriptomic responses to stimuli. In other cases, the question at hand might be more easily addressed using bulk RNA-seq or targeted gene expression experiments.

### scRNA-seq versus snRNA-seq

A thorough understanding of the different strategies and types of protocols (Fig. 1C) is essential before one can make an educated decision on which technology will best answer a specific biological question. A first example of an important choice is whether to profile the transcriptome from isolated nuclei or cells. When doing scRNA-seq, the most popular choice to obtain single cells from a plant organ requires the enzymatic digestion of cell walls and the generation of so-called protoplasts. There are a number of disadvantages of using protoplasts such as some tissues (e.g. sclerenchyma) and species (e.g. sorghum; Guillotin et al. 2023) are recalcitrant to cell wall digestion; enzymatic digestion affects the transcriptional status of the plant cells and could bias the outcome of experiments (Birnbaum et al. 2005); and the large size of protoplasts reduces their capture efficiency with most of the currently available commercialized single-cell platforms. Nucleus isolation followed by snRNA-seq gained traction in plant single-cell transcriptomics as well. However, the recovered data content per nucleus (e.g. Unique Molecular Identifier (UMI) or genes) is up to 10 (for UMIs) and 3 times (for genes) lower compared to scRNA-seq (Farmer et al. 2021; Guillotin et al. 2023). Furthermore, even though the transcriptome coverage is similarly efficient between scRNA-seq and snRNA-seq (e.g. 89% of all *Arabidopsis* expressed genes were captured in snRNA-seq data; Farmer et al. 2021; Guillotin et al. 2023), a nuclear transcriptome and a cellular transcriptome are not equivalent (Lee and Bailey-Serres 2019; Reynoso et al. 2019). For example, differences in abundance and composition between transcripts obtained from nuclear versus polyA RNA under hypoxia point toward nuclear transcript retention or enrichment as part of the stress response (Lee and Bailey-Serres 2019; Reynoso et al. 2019). Furthermore, the half-life of the transcripts (estimations range between 12 min and more than 24 h in *Arabidopsis* cells; Narsai et al. 2007) suggests that the cellular transcriptome is the result of the accumulation of the transcript synthesis over time, while the nuclear transcriptome is considered to accommodate faster to changes in gene activity. These differences are important to consider when selecting and later interpreting a single-cellular transcriptome versus a nuclear-based transcriptome and should be determined by the experimental system or biological question. Therefore, when studying, e.g. early stress responses of plant cells, a snRNA-seq could achieve higher resolution of rapid transcriptomic changes, while scRNA-seq might be more informative when understanding the biology
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Overview of plant single-cell and single-nucleus experiments. A) Number of publications describing sc/snRNA-seq data in the plant field per year from 2019 until 2022. B) Distribution of species used in these papers (n: 46). C) Overview of the different sc/sn technologies and their usage in the plant field. Example references used: a. (Farmer et al. 2021; Cervantes-Pérez et al. 2022; Conde et al. 2022; Neumann et al. 2022; Sun et al. 2022; Guillotin et al. 2023; Li et al. 2023; Liu et al. 2023); b. (Tian et al. 2020; Marand et al. 2021); c. (Wendrich et al. 2020; Graeff et al. 2021; Lopez-Anido et al. 2021; Ortiz-Ramírez et al. 2021; Wang, Huan, et al. 2021; Apelt et al. 2022; Otero et al. 2022; Kim et al. 2023); d. (Denyer et al. 2019; Jean-Baptiste et al. 2019; Ryu et al. 2019; Shulse et al. 2019; Turco et al. 2019; Zhang et al. 2019; Satterlee et al. 2020; Bezrutczyk et al. 2021; Chen et al. 2021; Gala et al. 2021; Kim et al. 2021; Liu et al. 2021; Ma et al. 2021; Yang et al. 2021; Zhang, Chen, et al. 2021; Zhang, Chen, and Wang 2021; Li et al. 2022; Shahan et al. 2022; Tao et al. 2022); e. (Kao et al. 2021; Picard et al. 2021; Sunaga-Franze et al. 2021; Abramson et al. 2022; Li et al. 2023); f. none; g. (Efroni et al. 2016; Lopez-Anido et al. 2021; Roszak et al. 2021; Serrano-Ron et al. 2021; Omary et al. 2022); h. (Zong et al. 2022); and i. (Nelms and Walbot 2019; Song et al. 2020; Xie et al. 2022).</figcaption></figure>

of a cell type or when studying cells that are enucleated at some stages of development (e.g. sieve element cells in the phloem cell lineage; Miyashima et al. 2019).

## Biological replicates in single-cell transcriptomics

As for all scientific observations, generating robust sc/sn data sets requires performance evaluation across multiple, independent biological replicates. We hereby note that a biological replicate relies on the independent growth, harvesting and processing of various plant samples. Any separation after protoplast or nucleus isolation cannot be classified as biological replicates and can only be reported as technical replicates. No standardized metrics are available within the community to evaluate reproducibility between replicates. We propose that a correlation coefficient of the average gene expressions among all cells would be an informative assay. Alternatively, one could compare the frequency of cell types or cell clusters across replicates. As such, we advise

analyzing cell cluster-specific differentially expressed genes and annotating each replicate separately, before merging the replicates and applying batch effect correction. Other parameters, e.g. Average Silhouette Width and Adjusted Rand Index, have been used to quantify cell-type purity assessments after batch effect correction and can also be informative to evaluate replicate robustness (Tran et al. 2020).

From a statistical point of view, independent biological replicates are unconditionally advised to increase the significance of biological data sets (Heumos et al. 2023). However, in many cases, replicates in sc/snRNA-seq experiments are currently performed to increase the total number of cells or nuclei analyzed, while the replicate information and a comparison between replicates are not necessarily incorporated in the actual statistical analysis. A statistical comparison among biological replicates is thus strongly advised to ensure high data quality and to prevent cluster formation based on replicate-specific artifacts. As such, biological
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Overview of plant single-cell and single-nucleus experiments. A) Number of publications describing sc/snRNA-seq data in the plant field per year from 2019 until 2022. B) Distribution of species used in these papers (n: 46). C) Overview of the different sc/sn technologies and their usage in the plant field. Example references used: a. (Farmer et al. 2021; Cervantes-Pérez et al. 2022; Conde et al. 2022; Neumann et al. 2022; Sun et al. 2022; Guillotin et al. 2023; Li et al. 2023; Liu et al. 2023); b. (Tian et al. 2020; Marand et al. 2021); c. (Wendrich et al. 2020; Graeff et al. 2021; Lopez-Anido et al. 2021; Ortiz-Ramírez et al. 2021; Wang, Huan, et al. 2021; Apelt et al. 2022; Otero et al. 2022; Kim et al. 2023); d. (Denyer et al. 2019; Jean-Baptiste et al. 2019; Ryu et al. 2019; Shulse et al. 2019; Turco et al. 2019; Zhang et al. 2019; Satterlee et al. 2020; Bezrutczyk et al. 2021; Chen et al. 2021; Gala et al. 2021; Kim et al. 2021; Liu et al. 2021; Ma et al. 2021; Yang et al. 2021; Zhang, Chen, et al. 2021; Zhang, Chen, and Wang 2021; Li et al. 2022; Shahan et al. 2022; Tao et al. 2022); e. (Kao et al. 2021; Picard et al. 2021; Sunaga-Franze et al. 2021; Abramson et al. 2022; Li et al. 2023); f. none; g. (Efroni et al. 2016; Lopez-Anido et al. 2021; Roszak et al. 2021; Serrano-Ron et al. 2021; Omary et al. 2022); h. (Zong et al. 2022); and i. (Nelms and Walbot 2019; Song et al. 2020; Xie et al. 2022).</figcaption></figure>

of a cell type or when studying cells that are enucleated at some stages of development (e.g. sieve element cells in the phloem cell lineage; Miyashima et al. 2019).

## Biological replicates in single-cell transcriptomics

As for all scientific observations, generating robust sc/sn data sets requires performance evaluation across multiple, independent biological replicates. We hereby note that a biological replicate relies on the independent growth, harvesting and processing of various plant samples. Any separation after protoplast or nucleus isolation cannot be classified as biological replicates and can only be reported as technical replicates. No standardized metrics are available within the community to evaluate reproducibility between replicates. We propose that a correlation coefficient of the average gene expressions among all cells would be an informative assay. Alternatively, one could compare the frequency of cell types or cell clusters across replicates. As such, we advise

analyzing cell cluster-specific differentially expressed genes and annotating each replicate separately, before merging the replicates and applying batch effect correction. Other parameters, e.g. Average Silhouette Width and Adjusted Rand Index, have been used to quantify cell-type purity assessments after batch effect correction and can also be informative to evaluate replicate robustness (Tran et al. 2020).

From a statistical point of view, independent biological replicates are unconditionally advised to increase the significance of biological data sets (Heumos et al. 2023). However, in many cases, replicates in sc/snRNA-seq experiments are currently performed to increase the total number of cells or nuclei analyzed, while the replicate information and a comparison between replicates are not necessarily incorporated in the actual statistical analysis. A statistical comparison among biological replicates is thus strongly advised to ensure high data quality and to prevent cluster formation based on replicate-specific artifacts. As such, biological
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Challenges in plant single-cell and single-nucleus transcriptomics. Overview of the different steps of performing sc/snRNA-seq in plant samples and summary of how the most important challenges discussed here influence this flowchart.</figcaption></figure>

replicates are imperative to add certainty on the reproducibility of the experiment. However, merely adding biological replicates does not remove transcriptional artifacts introduced during sample preparation in each of the replicates. One example is the effect of the enzymatic digestion needed to generate protoplasts or the procedures to extract nuclei on the transcriptome. Therefore, performing replicates by themselves does not provide sufficient confidence in the data to draw biological conclusions. To achieve confident biological interpretation, extensive downstream experimental validation is always required in the form of e.g. reporter lines, in situ hybridization, or spatial transcriptomics.

## Challenge 2: experimental variability during sample preparation

While the potential of sc/snRNA-seq for plant research is evident, its applicability depends largely on establishing reliable cell and nucleus isolation protocols. These protocols must support the generation of high-quality, high-yield nuclear and visible cellular suspensions within a short amount of time and must be compatible with downstream procedures (e.g. limited usage of PCR inhibitors like CaCl₂). The efficiency of protoplast generation from tomato roots for example was increased by optimizing the pH of the enzyme-containing buffer and, in part, also by using hand sections instead of intact tissues (Omary

et al. 2022). Preincubation in L-cysteine and sorbitol for roots of maize, sorghum, and Setaria improves enzymatic cell wall digestion and protoplast generation (Ortiz-Ramírez et al. 2018, 2021), while L-arginine positively influenced the survival rate of maize meristem protoplasts (Satterlee et al. 2020). In contrast, nuclei can be isolated from fresh (Farmer et al. 2021; Picard et al. 2021; Cervantes-Pérez et al. 2022; Conde et al. 2022; Sun et al. 2022; Liu et al. 2023), frozen (Sunaga-Franze et al. 2021; Abramson et al. 2022; Neumann et al. 2022; Li et al. 2023), or fixed (Kao et al. 2021) starting material, offering flexibility in terms of sample handling and preparation, while simultaneously securing dynamic transcriptional changes upon their rapid fixation. While nucleus isolation seems more straightforward to conduct than protoplast isolation, the assessment of nucleus quality prior to snRNA-seq library construction remains a difficult task. The leaking and clumping of isolated nuclei should be seen as a sign of breakage of the nuclear membrane leading to RNA leakage and the generation of low-quality libraries.

Overall, careful workflow optimization should include the following:

1) Visual assessment of tissue digestion or nucleus release through e.g. the observation of protoplasts/nuclei produced from all desired cell types, via cell wall digestion of fluorescently tagged cells of a particular cell type (if available) as a proxy or via gene expression quantification of cell-type markers in a pilot experiment.

2) Rapid and nondestructive sample cleanup strategies including washing steps (e.g. centrifugation and filtration), fluorescence (and image-based) activated cell/nucleus sorting (FACS/FANS), or microfluidic cell enrichment devices can increase the population of viable cells and the purity of cellular/nuclear suspensions. A careful analysis of nucleus shape will help to identify problems with RNA leakage.

3) Careful analysis of cell sizes for protoplasts since most commercial platforms have a cell size restriction that might introduce a bias in cell capturing and a preference for incorporating smaller over larger cells (see Challenge 3).

4) Quantification of cell viability by manual cell counting (upon staining with trypan blue or fluorescein diacetate) or with the help of automated cell counters.

The procedure of cell wall digestion itself (Birnbaum et al. 2005)—as well as external factors introduced during sample collection and generation (e.g. growth and harvesting conditions, enzyme concentration and activity, and temperature and timing)—affects cell viability, cell wall digestion efficiency, cell-type representation, and the transcriptional profiles of cells. One of the most promising developments for reducing experimental biases is the inclusion of a fixation step. Until now, scRNA-seq-compatible cell fixation protocols have mainly been described in mammalian research (Attar et al. 2018; Wohnhaas et al. 2019; Phan et al. 2021; Wang, Yu, and Wu 2021), but its application could drastically boost
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Challenges in plant single-cell and single-nucleus transcriptomics. Overview of the different steps of performing sc/snRNA-seq in plant samples and summary of how the most important challenges discussed here influence this flowchart.</figcaption></figure>

replicates are imperative to add certainty on the reproducibility of the experiment. However, merely adding biological replicates does not remove transcriptional artifacts introduced during sample preparation in each of the replicates. One example is the effect of the enzymatic digestion needed to generate protoplasts or the procedures to extract nuclei on the transcriptome. Therefore, performing replicates by themselves does not provide sufficient confidence in the data to draw biological conclusions. To achieve confident biological interpretation, extensive downstream experimental validation is always required in the form of e.g. reporter lines, in situ hybridization, or spatial transcriptomics.

## Challenge 2: experimental variability during sample preparation

While the potential of sc/snRNA-seq for plant research is evident, its applicability depends largely on establishing reliable cell and nucleus isolation protocols. These protocols must support the generation of high-quality, high-yield nuclear and visible cellular suspensions within a short amount of time and must be compatible with downstream procedures (e.g. limited usage of PCR inhibitors like CaCl₂). The efficiency of protoplast generation from tomato roots for example was increased by optimizing the pH of the enzyme-containing buffer and, in part, also by using hand sections instead of intact tissues (Omary

et al. 2022). Preincubation in L-cysteine and sorbitol for roots of maize, sorghum, and Setaria improves enzymatic cell wall digestion and protoplast generation (Ortiz-Ramírez et al. 2018, 2021), while L-arginine positively influenced the survival rate of maize meristem protoplasts (Satterlee et al. 2020). In contrast, nuclei can be isolated from fresh (Farmer et al. 2021; Picard et al. 2021; Cervantes-Pérez et al. 2022; Conde et al. 2022; Sun et al. 2022; Liu et al. 2023), frozen (Sunaga-Franze et al. 2021; Abramson et al. 2022; Neumann et al. 2022; Li et al. 2023), or fixed (Kao et al. 2021) starting material, offering flexibility in terms of sample handling and preparation, while simultaneously securing dynamic transcriptional changes upon their rapid fixation. While nucleus isolation seems more straightforward to conduct than protoplast isolation, the assessment of nucleus quality prior to snRNA-seq library construction remains a difficult task. The leaking and clumping of isolated nuclei should be seen as a sign of breakage of the nuclear membrane leading to RNA leakage and the generation of low-quality libraries.

Overall, careful workflow optimization should include the following:

1) Visual assessment of tissue digestion or nucleus release through e.g. the observation of protoplasts/nuclei produced from all desired cell types, via cell wall digestion of fluorescently tagged cells of a particular cell type (if available) as a proxy or via gene expression quantification of cell-type markers in a pilot experiment.

2) Rapid and nondestructive sample cleanup strategies including washing steps (e.g. centrifugation and filtration), fluorescence (and image-based) activated cell/nucleus sorting (FACS/FANS), or microfluidic cell enrichment devices can increase the population of viable cells and the purity of cellular/nuclear suspensions. A careful analysis of nucleus shape will help to identify problems with RNA leakage.

3) Careful analysis of cell sizes for protoplasts since most commercial platforms have a cell size restriction that might introduce a bias in cell capturing and a preference for incorporating smaller over larger cells (see Challenge 3).

4) Quantification of cell viability by manual cell counting (upon staining with trypan blue or fluorescein diacetate) or with the help of automated cell counters.

The procedure of cell wall digestion itself (Birnbaum et al. 2005)—as well as external factors introduced during sample collection and generation (e.g. growth and harvesting conditions, enzyme concentration and activity, and temperature and timing)—affects cell viability, cell wall digestion efficiency, cell-type representation, and the transcriptional profiles of cells. One of the most promising developments for reducing experimental biases is the inclusion of a fixation step. Until now, scRNA-seq-compatible cell fixation protocols have mainly been described in mammalian research (Attar et al. 2018; Wohnhaas et al. 2019; Phan et al. 2021; Wang, Yu, and Wu 2021), but its application could drastically boost
the plant single-cell field by massively reducing the effects of external factors during sample processing, including the generation of protoplasts. Indeed, protoplast isolation efficiency increased when plant tissues were fixed and digested at optimal enzyme activity temperature (Marchant et al. 2022). However, concerns about tissue fixation on protoplast shape (Marchant et al. 2022) and the sequencing results have been reported, motivated in particular by the reduction of cDNA yields and biases toward 3'-end enrichment (Wang, Yu, and Wu 2021). Despite these limitations, the potential gains for the field could be major, warranting dedicated investments in tissue fixation approaches.

## Challenge 3: biases and specificities of commercial platforms for plant single-cell transcriptomic samples

The most popular commercial platforms and scRNA-seq protocols used for plant samples rely on microfluidic droplet-based cell compartmentalization or nanowell-based cell separation (Fig. 1C). Techniques that can be performed by manual handling, such as combinatorial barcoding (Cao et al. 2017; Rosenberg et al. 2018), are rapidly expanding in the animal field but have yet to be shown in use for plant samples. The choice of the sample processing method or platform must be taken carefully to allow uniform cell size capture rate, resolvability, and, if necessary, a sample multiplexing option or flexibility toward cell capturing and lysis steps. Droplet-based platforms allow fast cell/nucleus processing but offer limited flexibility regarding the cell preparation workflow. Also, the level of pressure imposed on the sample when creating the emulsion could cause the bursting of cells into the droplet-based platform. Well-based methods like SMART-Seq2 (Lopez-Anido et al. 2021) and platforms used with plant samples, such as BD Rhapsody (Zong et al. 2022) or iCELL8 (Sunaga-Franze et al. 2021), require longer cell processing protocols but offer more flexibility during the sample processing. However, the compatibility of commercial platforms to handle the size and fragility of plant protoplasts is not necessarily evaluated. Plant cell sizes typically lie in the range of 10 to 80 µm, with even larger values observed for endoreduplicated cells, which is far above the recommended cell size maxima from current technology providers (~40 µm). Furthermore, cell size heterogeneity can create cell capture biases, because droplet-based techniques favor smaller cells, while the well sizes in nanowell techniques must be fine-tuned to reduce the possibilities of doublets from smaller cells while still allowing capture and processing of larger cells. The consequences could be high multiplet rates and/or imbalanced cell-type/stage representation. Careful optimization of the maximum cell loading capacity, loading speed, cell compartmentalization time, and the number of washing steps is necessary depending on the platform of choice. Identification of cell type or stage capture rates, however, requires in vivo experimental

validation by quantification of cell types or tracing of developmental cell stages. This validation has been done by comparing cell numbers per cell type between scRNA-seq data and cell counting via imaging (Wendrich et al. 2020), but it could also be achieved by spiking in a fixed ratio of cell types using transgenic marker lines.

Furthermore, a detailed plant-specific benchmark study comparing the commercially available platforms and kits is urgently needed to evaluate the benefits and pitfalls when applied to plant samples. Similar benchmark studies using human and mouse cell lines allowed practicality and financial comparisons of common methods, while also comparing cell capture rates and technical bias across cells with distinct cell properties (Mereu et al. 2020). Ideally, a benchmark study should cover numerous species and tissue combinations to allow the establishment of quality standards independent of the species or tissue used. A plant-specific benchmark study might in addition also focus on comparing the ability to recreate developmental or spatial cell states, because most plant single-cell samples offer a chance to also capture and compare developmental cell trajectories. Given the sensitive nature of protoplasts, it would also be important to address the effect of total sample processing durations, as well as the abundance of ambient RNA due to for example protoplast bursting.

## Challenge 4: establishing an efficient sequencing strategy

### Full-length versus 3'- or 5'-end transcript sequencing
Two approaches for library preparation are currently used in single-cell methodology, namely full-length transcript coverage and 3'- or 5'-end transcript coverage. While most reported plant single-cell transcriptomic data sets today (Fig. 1C) were generated using 3'end transcript coverage, full-length transcript sequencing (used for example in Lopez-Anido et al. 2021) bares huge application potential in plant research, as it can help in improving transcriptome and epitranscriptome resources (reviewed in Shen et al. 2023) and in providing cell-type-selective isoform usage (shown for example during neuronal maturation in mouse embryos; Lebrigand et al. 2020).

### Cell number and sequencing depth

Although cDNA amount and profile after library preparation are used as a proxy for the overall quality of the sc/snRNA-seq library, sequencing followed by analysis currently remains the only way to fully estimate its quality and biological value. Two major issues that require careful consideration during the experimental planning and sample optimization are the number of cells/nuclei needed for optimal coverage of the cell type(s) of interest and the aimed sequencing depth per cell/nucleus. The number of cells/nuclei in published cell atlases is moving from thousands to hundreds of thousands. Increasing numbers is indeed beneficial for predicting novel marker genes, because it allows better
the plant single-cell field by massively reducing the effects of external factors during sample processing, including the generation of protoplasts. Indeed, protoplast isolation efficiency increased when plant tissues were fixed and digested at optimal enzyme activity temperature (Marchant et al. 2022). However, concerns about tissue fixation on protoplast shape (Marchant et al. 2022) and the sequencing results have been reported, motivated in particular by the reduction of cDNA yields and biases toward 3'-end enrichment (Wang, Yu, and Wu 2021). Despite these limitations, the potential gains for the field could be major, warranting dedicated investments in tissue fixation approaches.

## Challenge 3: biases and specificities of commercial platforms for plant single-cell transcriptomic samples

The most popular commercial platforms and scRNA-seq protocols used for plant samples rely on microfluidic droplet-based cell compartmentalization or nanowell-based cell separation (Fig. 1C). Techniques that can be performed by manual handling, such as combinatorial barcoding (Cao et al. 2017; Rosenberg et al. 2018), are rapidly expanding in the animal field but have yet to be shown in use for plant samples. The choice of the sample processing method or platform must be taken carefully to allow uniform cell size capture rate, resolvability, and, if necessary, a sample multiplexing option or flexibility toward cell capturing and lysis steps. Droplet-based platforms allow fast cell/nucleus processing but offer limited flexibility regarding the cell preparation workflow. Also, the level of pressure imposed on the sample when creating the emulsion could cause the bursting of cells into the droplet-based platform. Well-based methods like SMART-Seq2 (Lopez-Anido et al. 2021) and platforms used with plant samples, such as BD Rhapsody (Zong et al. 2022) or iCELL8 (Sunaga-Franze et al. 2021), require longer cell processing protocols but offer more flexibility during the sample processing. However, the compatibility of commercial platforms to handle the size and fragility of plant protoplasts is not necessarily evaluated. Plant cell sizes typically lie in the range of 10 to 80 µm, with even larger values observed for endoreduplicated cells, which is far above the recommended cell size maxima from current technology providers (~40 µm). Furthermore, cell size heterogeneity can create cell capture biases, because droplet-based techniques favor smaller cells, while the well sizes in nanowell techniques must be fine-tuned to reduce the possibilities of doublets from smaller cells while still allowing capture and processing of larger cells. The consequences could be high multiplet rates and/or imbalanced cell-type/stage representation. Careful optimization of the maximum cell loading capacity, loading speed, cell compartmentalization time, and the number of washing steps is necessary depending on the platform of choice. Identification of cell type or stage capture rates, however, requires in vivo experimental

validation by quantification of cell types or tracing of developmental cell stages. This validation has been done by comparing cell numbers per cell type between scRNA-seq data and cell counting via imaging (Wendrich et al. 2020), but it could also be achieved by spiking in a fixed ratio of cell types using transgenic marker lines.

Furthermore, a detailed plant-specific benchmark study comparing the commercially available platforms and kits is urgently needed to evaluate the benefits and pitfalls when applied to plant samples. Similar benchmark studies using human and mouse cell lines allowed practicality and financial comparisons of common methods, while also comparing cell capture rates and technical bias across cells with distinct cell properties (Mereu et al. 2020). Ideally, a benchmark study should cover numerous species and tissue combinations to allow the establishment of quality standards independent of the species or tissue used. A plant-specific benchmark study might in addition also focus on comparing the ability to recreate developmental or spatial cell states, because most plant single-cell samples offer a chance to also capture and compare developmental cell trajectories. Given the sensitive nature of protoplasts, it would also be important to address the effect of total sample processing durations, as well as the abundance of ambient RNA due to for example protoplast bursting.

## Challenge 4: establishing an efficient sequencing strategy

### Full-length versus 3'- or 5'-end transcript sequencing

Two approaches for library preparation are currently used in single-cell methodology, namely full-length transcript coverage and 3'- or 5'-end transcript coverage. While most reported plant single-cell transcriptomic data sets today (Fig. 1C) were generated using 3'end transcript coverage, full-length transcript sequencing (used for example in Lopez-Anido et al. 2021) bares huge application potential in plant research, as it can help in improving transcriptome and epitranscriptome resources (reviewed in Shen et al. 2023) and in providing cell-type-selective isoform usage (shown for example during neuronal maturation in mouse embryos; Lebrigand et al. 2020).

### Cell number and sequencing depth

Although cDNA amount and profile after library preparation are used as a proxy for the overall quality of the sc/snRNA-seq library, sequencing followed by analysis currently remains the only way to fully estimate its quality and biological value. Two major issues that require careful consideration during the experimental planning and sample optimization are the number of cells/nuclei needed for optimal coverage of the cell type(s) of interest and the aimed sequencing depth per cell/nucleus. The number of cells/nuclei in published cell atlases is moving from thousands to hundreds of thousands. Increasing numbers is indeed beneficial for predicting novel marker genes, because it allows better
<table><caption>**Table 1.** Necessary reported information to allow evaluation and repetition of a plant sc/sn experiment</caption><thead><tr><th></th><th>Details</th><th>Recommendations</th></tr></thead><tbody><tr><td><strong>Biological material</strong></td><td>Species<br>Accession<br>Genotype<br>Tissue type<br>Detailed growth conditions<br>Harvest conditions</td><td>e.g. <em>Arabidopsis thaliana</em>, <em>Zea mays</em><br>e.g. Col-0<br>e.g. WT or mutant background<br>e.g. root, leaf, stem, and seed<br>e.g. temperature, light conditions, and medium<br>e.g. age of plants, time of day, and amount harvested</td></tr><tr><td><strong>Sample preparation</strong></td><td>Isolation protocol<br>Tissue dissection<br>Fixation<br>Cell/nucleus enrichment<br>Total sample preparation time</td><td>Short description of the way the sample was isolated<br>e.g. razor blades, needles, and tissue homogenizer<br>Short description of the way the sample was fixed if this was done<br>e.g. sucrose gradient, FACS (incl. model, nozzle size, and temperature)<br>For cells: &lt;90 min for <em>Arabidopsis</em> roots (from material harvest to cell loading)*<br>*Duration may increase depending on starting material and time needed for optimal tissue digestion<br>For nuclei: 30 to 60 min (depending whether a nucleus enrichment step is included)</td></tr><tr><td></td><td>Estimated cell/nucleus number loaded</td><td>An estimation of the amount of cells or nuclei loaded, based on the cell/nucleus concentration and volume that was loaded</td></tr><tr><td></td><td>Instrument/method/kit<br>Cell viability test</td><td>e.g. 10× Genomics 3' v3.1, BD Rhapsody WTA<br>For cells: trypan blue, fluorescein diacetate, calcein, propidium iodide, 4',6-diamidino-2-phenylindol<br>For nuclei: not applicable</td></tr><tr><td><strong>Libraries</strong></td><td>Library construction<br>Amplification method<br>End bias</td><td>Protocol and revision/version that was followed, e.g. CG000204 Rev D for 3' v3.1<br>e.g. number of PCR cycles used for cDNA amplification<br>e.g. 3'-end mainly; excess of rRNA or TSO sequences</td></tr><tr><td><strong>Sequence results</strong></td><td>Instrument/method<br>Library layout/paired-end<br>N° sequenced reads</td><td>e.g. NovaSeq, NextSeq, ONT, and DNBSEQ<br>Consider to use standardized library structures (Booeshaghi et al. 2023)<br>20,000 to 50,000/cell for RNA or more<br>20,000 to 40,000/nucleus for RNA or more</td></tr><tr><td><strong>Raw data</strong></td><td>Reference genome<br>Annotation version<br>Mapping method (incl. software, customized settings)<br>Mapping efficiency</td><td>Link to Ensembl Plant fasta file, JGI, NCBI, PLAZA<br>If custom annotation, also include.gtf/.gff/.gff3 files<br>e.g. STAR (cellranger)<br>&gt;85% for <em>Arabidopsis</em><br>*Value may be lower in other species</td></tr><tr><td></td><td>Sequencing saturation<br>Estimation of ambient RNA</td><td>&gt;50%<br>Fraction of reads in cells &gt; 60% for scRNA-seq<br>Fraction of reads in cells &gt; 50% for snRNA-seq</td></tr><tr><td><strong>Processed data</strong></td><td>Imputation method and settings<br>N° captured cells/nuclei<br>N° high quality cells/nuclei<br>Filter criteria: % mitochondrial reads/cell or nucleus<br>Filter criteria: % chloroplast reads/cell or nucleus<br>Filter criteria:<br>Minimum N° UMI/cell or nucleus<br>N° total detected transcripts</td><td>If relevant<br>e.g. 60% of estimated number<br>e.g. 20% of estimated number<br>&lt;10% for scRNA-seq*<br>close to 0% for snRNA-seq*<br>*value may deviate depending on biological context<br>&lt;5% to 10% for scRNA-seq*<br>close to 0% for snRNA-seq*<br>*values may be higher depending on biological context<br>&gt;1000 for scRNA-seq&gt;400 for snRNA-seq<br>Dependent on the species but e.g. 60% of total number of transcripts in the annotated genome</td></tr><tr><td></td><td>Doublet rate<br>Replicate comparisons<br>Batch correction method for merging (incl. reasoning for batch correction)<br>Additional processing</td><td>Estimates according to 10× Genomics user guide based on number of loaded cells/ nuclei<br>Provide coefficient correlation of the most variable genes or cluster-specific genes between independent replicates or compare pseudo-bulk from sc/snRNA-seq versus bulk RNA-seq<br>e.g. Seurat, Harmony</td></tr><tr><td><strong>Validation</strong></td><td>Method of automatic annotation of clusters<br>Method of manual annotation (markers, gene function info)</td><td>e.g. removal of protoplast-induced genes, cell cycle regression, noise (ambient RNA) removal, cluster membership bias between replicates, removal of low quality clusters<br>e.g. Label transfer, spatial transcriptomics<br>e.g. Marker genes, orthologous, correlation with bulk RNA-seq and microarray data; Index of Cell Identity calculation</td></tr></tbody></table>

(continued)
<table><caption>**Table 1.** Necessary reported information to allow evaluation and repetition of a plant sc/sn experiment</caption><thead><tr><th></th><th>Details</th><th>Recommendations</th></tr></thead><tbody><tr><td><strong>Biological material</strong></td><td>Species<br>Accession<br>Genotype<br>Tissue type<br>Detailed growth conditions<br>Harvest conditions</td><td>e.g. <em>Arabidopsis thaliana</em>, <em>Zea mays</em><br>e.g. Col-0<br>e.g. WT or mutant background<br>e.g. root, leaf, stem, and seed<br>e.g. temperature, light conditions, and medium<br>e.g. age of plants, time of day, and amount harvested</td></tr><tr><td><strong>Sample preparation</strong></td><td>Isolation protocol<br>Tissue dissection<br>Fixation<br>Cell/nucleus enrichment<br>Total sample preparation time</td><td>Short description of the way the sample was isolated<br>e.g. razor blades, needles, and tissue homogenizer<br>Short description of the way the sample was fixed if this was done<br>e.g. sucrose gradient, FACS (incl. model, nozzle size, and temperature)<br>For cells: &lt;90 min for <em>Arabidopsis</em> roots (from material harvest to cell loading)*<br>*Duration may increase depending on starting material and time needed for optimal tissue digestion<br>For nuclei: 30 to 60 min (depending whether a nucleus enrichment step is included)</td></tr><tr><td></td><td>Estimated cell/nucleus number loaded</td><td>An estimation of the amount of cells or nuclei loaded, based on the cell/nucleus concentration and volume that was loaded</td></tr><tr><td></td><td>Instrument/method/kit<br>Cell viability test</td><td>e.g. 10× Genomics 3' v3.1, BD Rhapsody WTA<br>For cells: trypan blue, fluorescein diacetate, calcein, propidium iodide, 4',6-diamidino-2-phenylindol<br>For nuclei: not applicable</td></tr><tr><td><strong>Libraries</strong></td><td>Library construction<br>Amplification method<br>End bias</td><td>Protocol and revision/version that was followed, e.g. CG000204 Rev D for 3' v3.1<br>e.g. number of PCR cycles used for cDNA amplification<br>e.g. 3'-end mainly; excess of rRNA or TSO sequences</td></tr><tr><td><strong>Sequence results</strong></td><td>Instrument/method<br>Library layout/paired-end<br>N° sequenced reads</td><td>e.g. NovaSeq, NextSeq, ONT, and DNBSEQ<br>Consider to use standardized library structures (Booeshaghi et al. 2023)<br>20,000 to 50,000/cell for RNA or more<br>20,000 to 40,000/nucleus for RNA or more</td></tr><tr><td><strong>Raw data</strong></td><td>Reference genome<br>Annotation version<br>Mapping method (incl. software, customized settings)<br>Mapping efficiency</td><td>Link to Ensembl Plant fasta file, JGI, NCBI, PLAZA<br>If custom annotation, also include.gtf/.gff/.gff3 files<br>e.g. STAR (cellranger)<br>&gt;85% for <em>Arabidopsis</em><br>*Value may be lower in other species</td></tr><tr><td></td><td>Sequencing saturation<br>Estimation of ambient RNA</td><td>&gt;50%<br>Fraction of reads in cells &gt; 60% for scRNA-seq<br>Fraction of reads in cells &gt; 50% for snRNA-seq</td></tr><tr><td><strong>Processed data</strong></td><td>Imputation method and settings<br>N° captured cells/nuclei<br>N° high quality cells/nuclei<br>Filter criteria: % mitochondrial reads/cell or nucleus<br>Filter criteria: % chloroplast reads/cell or nucleus<br>Filter criteria:<br>Minimum N° UMI/cell or nucleus<br>N° total detected transcripts</td><td>If relevant<br>e.g. 60% of estimated number<br>e.g. 20% of estimated number<br>&lt;10% for scRNA-seq*<br>close to 0% for snRNA-seq*<br>*value may deviate depending on biological context<br>&lt;5% to 10% for scRNA-seq*<br>close to 0% for snRNA-seq*<br>*values may be higher depending on biological context<br>&gt;1000 for scRNA-seq&gt;400 for snRNA-seq<br>Dependent on the species but e.g. 60% of total number of transcripts in the annotated genome</td></tr><tr><td></td><td>Doublet rate<br>Replicate comparisons<br>Batch correction method for merging (incl. reasoning for batch correction)<br>Additional processing</td><td>Estimates according to 10× Genomics user guide based on number of loaded cells/ nuclei<br>Provide coefficient correlation of the most variable genes or cluster-specific genes between independent replicates or compare pseudo-bulk from sc/snRNA-seq versus bulk RNA-seq<br>e.g. Seurat, Harmony</td></tr><tr><td><strong>Validation</strong></td><td>Method of automatic annotation of clusters<br>Method of manual annotation (markers, gene function info)</td><td>e.g. removal of protoplast-induced genes, cell cycle regression, noise (ambient RNA) removal, cluster membership bias between replicates, removal of low quality clusters<br>e.g. Label transfer, spatial transcriptomics<br>e.g. Marker genes, orthologous, correlation with bulk RNA-seq and microarray data; Index of Cell Identity calculation</td></tr></tbody></table>

(continued)
<table><caption>**Table 1.** (continued)</caption><thead><tr><th></th><th>Details</th><th>Recommendations</th></tr></thead><tbody><tr><td rowspan="6"><strong>Data availability</strong></td><td>Verification in planta (e.g. number of markers used for validation)</td><td>e.g. Spatial transcriptomics; RNA in situ hybridization; promoter fusions</td></tr><tr><td>Analysis scripts & codes (GitHub)</td><td>If relevant</td></tr><tr><td>Excel Tables DEG for each cluster</td><td>e.g. Lists for each cell type/developmental stage from FindMarkers (Seurat)</td></tr><tr><td>Objects/count matrix in repository (which one, where?)</td><td>e.g. use NCBI GEO to store count matrices and Seurat object</td></tr><tr><td>On-line tool/browser URL</td><td>List the URL if relevant</td></tr><tr><td>Cell-level metadata table</td><td>Include cell-type annotations for each cell barcode</td></tr><tr><td><strong>Additional</strong></td><td>additional comments from the authors</td><td>e.g. annotation/counting of rRNA, allow for rRNA estimation, consideration of intronic reads</td></tr></tbody></table>

A downloadable empty version for use in publications can be found as Supplementary Table S1. If numerical values deviate from the recommended numbers provided below, an explanation should be provided. The numerical values in this table are derived from available studies in plants, most of which originated from *Arabidopsis*.

coverage of rare populations and lowly expressed genes. This higher coverage in turn helps to outweigh the bias of differential expression analysis tools toward highly expressed genes (Squair et al. 2021) and data set-specific noise (Fischer and Gillis 2021). However, at what point does a cell atlas fully capture the cellular diversity of its samples? The meristematic region of an *Arabidopsis* root has about 3,000 to 4,000 cells, meaning that a data set of 100,000 cells sufficiently covers each cell about 20 to 30 times, assuming that all cells are equally represented in the data set. Alternatively, when a specific cell type is isolated from a tissue by upstream cell sorting, an atlas of 2,000 cells could already be saturating. As a rule of thumb, when the cell-type composition of the tissue of interest is known, the minimum number of cells that need to be analyzed can be estimated from the probability to robustly capture the rarest cell type(s) (e.g. https://satijalab.org/howmanycells/). For example, to obtain at least 10 quiescent center cells (estimated to represent 0.1% of all cells within the *Arabidopsis* root meristem according to Cartwright et al. 2009; Shahan et al. 2022) with 95% confidence, one would need to profile 15K to 20K cells.

Additionally, the required sequencing depth must be adapted depending on the biological question, the tissue complexity, and the sample quality. However, it is recommended that optimal coverage is given with 1 read/cell/gene (Zhang et al. 2020). Alternatively, sequencing can be staged by first initiating a shallow sequencing of the library (e.g. up to 10,000 reads/cell) before deeper sequencing (e.g. aiming for 50,000 reads/cell). Such shallow sequencing allows evaluating the performance of cell cluster analysis and annotation and is sufficient to capture the entire cell-type heterogeneity of the sample (Zhang et al. 2020). Another common suggestion for a preliminary sample quality control is to analyze the expression of a gene subset related to a biological question (Zhang et al. 2020). Sequencing even less for testing the quality of the library is possible but will affect the retrieval of cell types with a lower number of transcripts, which might be lost within the background of empty droplets if the sequencing is too shallow. The desired final read depth will depend on the goal of the experiment but should ensure sequencing enough cells at a sequencing depth that captures individual events robustly (e.g. 50% sequencing saturation; Table 1).

When assessing the most cost-effective sc/snRNA-seq technologies, should one profile many cells/nuclei but have shallow sequencing or should one profile fewer cells/nuclei but with deeper sequencing? In many cases, the ideal scenario will be something in-between assuming that the researchers are working on high-quality cells/nuclei. Nevertheless, a choice toward either a higher number of cells/nuclei or higher sequencing depth can be made depending on the biological question, the quality of the biological entities used for the analysis, and the relative abundance of each cell type composing the organ. If the aim is to generate an atlas potentially uncovering rare cell types, a better strategy would be to profile many cells/nuclei with a lower sequencing saturation. However, a minimal depth in sequencing (Table 1) upon maximizing cell/nucleus quality is still necessary to ensure that low-abundance transcripts that define rare cell types are captured and to saturate the transcriptome of the sample. Validation of high-throughput technologies in the plant field that enable to access the transcriptomes of hundreds of thousands and even up to 1 million cells or nuclei, combined with the ongoing expansion of sequencing capabilities and the decrease in sequencing costs, could help to overcome this dilemma. If the goal is to do functional gene discovery and generate gene regulatory networks for example (Ferrari et al. 2022), a high sequencing saturation per cell/nucleus favors the discovery of low-abundance gene transcripts. To achieve this goal, one can select a subset of genes related to the biological question and adjust the sequencing depth until at least 1 read/cell for each of those genes is reached (Zhang et al. 2020). In all cases, we recommend optimizing cell and nucleus isolation methods to ensure the capture of the largest number of transcripts from each biological entity.

### Challenge 5: generating an expression matrix from high-quality cells/nuclei

Once mRNA sequence reads are obtained, reads are mapped to the genome and ultimately to genes and cells of origin using a reference genome and UMI and cellular barcode information. Standard data analysis workflows further include quality control filtering, quantification of gene expression in each cell,
<table><caption>**Table 1.** (continued)</caption><thead><tr><th></th><th>Details</th><th>Recommendations</th></tr></thead><tbody><tr><td rowspan="6"><strong>Data availability</strong></td><td>Verification in planta (e.g. number of markers used for validation)</td><td>e.g. Spatial transcriptomics; RNA in situ hybridization; promoter fusions</td></tr><tr><td>Analysis scripts & codes (GitHub)</td><td>If relevant</td></tr><tr><td>Excel Tables DEG for each cluster</td><td>e.g. Lists for each cell type/developmental stage from FindMarkers (Seurat)</td></tr><tr><td>Objects/count matrix in repository (which one, where?)</td><td>e.g. use NCBI GEO to store count matrices and Seurat object</td></tr><tr><td>On-line tool/browser URL</td><td>List the URL if relevant</td></tr><tr><td>Cell-level metadata table</td><td>Include cell-type annotations for each cell barcode</td></tr><tr><td><strong>Additional</strong></td><td>additional comments from the authors</td><td>e.g. annotation/counting of rRNA, allow for rRNA estimation, consideration of intronic reads</td></tr></tbody></table>

A downloadable empty version for use in publications can be found as Supplementary Table S1. If numerical values deviate from the recommended numbers provided below, an explanation should be provided. The numerical values in this table are derived from available studies in plants, most of which originated from *Arabidopsis*.

coverage of rare populations and lowly expressed genes. This higher coverage in turn helps to outweigh the bias of differential expression analysis tools toward highly expressed genes (Squair et al. 2021) and data set-specific noise (Fischer and Gillis 2021). However, at what point does a cell atlas fully capture the cellular diversity of its samples? The meristematic region of an *Arabidopsis* root has about 3,000 to 4,000 cells, meaning that a data set of 100,000 cells sufficiently covers each cell about 20 to 30 times, assuming that all cells are equally represented in the data set. Alternatively, when a specific cell type is isolated from a tissue by upstream cell sorting, an atlas of 2,000 cells could already be saturating. As a rule of thumb, when the cell-type composition of the tissue of interest is known, the minimum number of cells that need to be analyzed can be estimated from the probability to robustly capture the rarest cell type(s) (e.g. https://satijalab.org/howmanycells/). For example, to obtain at least 10 quiescent center cells (estimated to represent 0.1% of all cells within the *Arabidopsis* root meristem according to Cartwright et al. 2009; Shahan et al. 2022) with 95% confidence, one would need to profile 15K to 20K cells.

Additionally, the required sequencing depth must be adapted depending on the biological question, the tissue complexity, and the sample quality. However, it is recommended that optimal coverage is given with 1 read/cell/gene (Zhang et al. 2020). Alternatively, sequencing can be staged by first initiating a shallow sequencing of the library (e.g. up to 10,000 reads/cell) before deeper sequencing (e.g. aiming for 50,000 reads/cell). Such shallow sequencing allows evaluating the performance of cell cluster analysis and annotation and is sufficient to capture the entire cell-type heterogeneity of the sample (Zhang et al. 2020). Another common suggestion for a preliminary sample quality control is to analyze the expression of a gene subset related to a biological question (Zhang et al. 2020). Sequencing even less for testing the quality of the library is possible but will affect the retrieval of cell types with a lower number of transcripts, which might be lost within the background of empty droplets if the sequencing is too shallow. The desired final read depth will depend on the goal of the experiment but should ensure sequencing enough cells at a sequencing depth that captures individual events robustly (e.g. 50% sequencing saturation; Table 1).

When assessing the most cost-effective sc/snRNA-seq technologies, should one profile many cells/nuclei but have shallow sequencing or should one profile fewer cells/nuclei but with deeper sequencing? In many cases, the ideal scenario will be something in-between assuming that the researchers are working on high-quality cells/nuclei. Nevertheless, a choice toward either a higher number of cells/nuclei or higher sequencing depth can be made depending on the biological question, the quality of the biological entities used for the analysis, and the relative abundance of each cell type composing the organ. If the aim is to generate an atlas potentially uncovering rare cell types, a better strategy would be to profile many cells/nuclei with a lower sequencing saturation. However, a minimal depth in sequencing (Table 1) upon maximizing cell/nucleus quality is still necessary to ensure that low-abundance transcripts that define rare cell types are captured and to saturate the transcriptome of the sample. Validation of high-throughput technologies in the plant field that enable to access the transcriptomes of hundreds of thousands and even up to 1 million cells or nuclei, combined with the ongoing expansion of sequencing capabilities and the decrease in sequencing costs, could help to overcome this dilemma. If the goal is to do functional gene discovery and generate gene regulatory networks for example (Ferrari et al. 2022), a high sequencing saturation per cell/nucleus favors the discovery of low-abundance gene transcripts. To achieve this goal, one can select a subset of genes related to the biological question and adjust the sequencing depth until at least 1 read/cell for each of those genes is reached (Zhang et al. 2020). In all cases, we recommend optimizing cell and nucleus isolation methods to ensure the capture of the largest number of transcripts from each biological entity.

### Challenge 5: generating an expression matrix from high-quality cells/nuclei

Once mRNA sequence reads are obtained, reads are mapped to the genome and ultimately to genes and cells of origin using a reference genome and UMI and cellular barcode information. Standard data analysis workflows further include quality control filtering, quantification of gene expression in each cell,
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Workflow for sc/snRNA-seq data analysis.</figcaption></figure>

clustering, and visualization of cells based on transcriptomic similarity (Fig. 3). While the recommendations below are given with the intention to standardize sample quality parameters within the field (see recommendations in Table 1), we want to highlight that a sn/scRNA-seq experiment (with a read depth of 30K reads/nucleus or cell) should typically result in the capture of 1,000 expressed genes/nucleus or 3,000 expressed genes/cell, respectively. Furthermore, utilizing universal preprocessing pipelines for single-cell genomic data (e.g. Booshaghi et al. 2023) can help in streamlining cell quality filtering and enhance data reproducibility in the future.

## Read mapping

Plant genomes in general and crop genomes in particular are poorly annotated compared to e.g. human or mouse genomes. Moreover, due to frequent whole-genome duplications (Fox et al. 2020), many plants are polyploid and thus contain multiple similar copies of each gene. Structural annotations of genes are especially important for most droplet-based technologies using a 3' capturing strategy to properly map the sequencing reads and to quantify transcript abundance. The distribution of the mapped reads on the genome can give an indication about the quality of the annotations. Poor mapping efficiency consequently causes gene loss, which can be dramatic especially for popular 3'-based single-cell technologies. Even for *Arabidopsis* scRNA-seq data sets, mapping rates vary but should be e.g. >85% (Table 1). A high percentage of reads mapping to intergenic regions (e.g. >20%) can be an indication either that not all genes are annotated or that the annotated 3' UTR regions should be longer. As a note, it is important to mention that the sequencing reads generated upon conducting a snRNA-seq experiment should be mapped against the exonic and intronic sequences of the annotated transcripts, a reflection of the capture of spliced and unspliced transcripts.

## Removal of low-quality cells/nuclei

After read mapping, low-quality cells/nuclei (e.g. cells/nuclei with low number of UMs and genes) need to be filtered out (see numerical recommendations in Table 1). Impairments in applying low-quality cell filtering may increase the noise in the data set and reduce the accuracy in downstream analysis, including cell clustering or erroneous identification of cell types (Fig. 4, A and B). General filter parameters will depend on the application and sample type (i.e. cells versus nuclei). For plants, it is good practice to exclude cells with high mitochondrial (e.g. >10%) and chloroplast (e.g. >5% to 10%) reads (Table 1). Such cases may indicate cells under stress because of perturbations during the sample preparation. These values might need to be adapted when studying highly or lowly metabolically active cells or cells undergoing e.g. programmed cell death. If no mitochondria or chloroplast gene is available, plotting the number of genes versus the total UMI count can be used instead to show cells with low data content. The barcode-rank-plot ("knee plot") is a commonly used tool to determine sample quality by ranking all barcodes according to their UMI content (Fig. 4A). A sample of low quality can be identified if there is no clearly defined boundary between barcodes with high UMI content and barcodes with significantly lower UMI content (Fig. 4C). Similarly, poor cluster separation in a sample that fails to differentiate cell types indicates an insufficient amount of transcript content per cell/nuclei.

Ambient RNA and the presence of empty droplets and doublets—both associated with droplet and nanowell-based technologies—can also lead to noise in the expression matrices, inaccurate cell clusters, and falsely differentially expressed genes. Therefore, it is important to ensure low amount of ambient RNA introduced during sample preparation (e.g. by mixing cells/nuclei from multiple species) and to optimize cell loading concentration (Fig. 4C).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Workflow for sc/snRNA-seq data analysis.</figcaption></figure>

clustering, and visualization of cells based on transcriptomic similarity (Fig. 3). While the recommendations below are given with the intention to standardize sample quality parameters within the field (see recommendations in Table 1), we want to highlight that a sn/scRNA-seq experiment (with a read depth of 30K reads/nucleus or cell) should typically result in the capture of 1,000 expressed genes/nucleus or 3,000 expressed genes/cell, respectively. Furthermore, utilizing universal preprocessing pipelines for single-cell genomic data (e.g. Booshaghi et al. 2023) can help in streamlining cell quality filtering and enhance data reproducibility in the future.

## Read mapping

Plant genomes in general and crop genomes in particular are poorly annotated compared to e.g. human or mouse genomes. Moreover, due to frequent whole-genome duplications (Fox et al. 2020), many plants are polyploid and thus contain multiple similar copies of each gene. Structural annotations of genes are especially important for most droplet-based technologies using a 3' capturing strategy to properly map the sequencing reads and to quantify transcript abundance. The distribution of the mapped reads on the genome can give an indication about the quality of the annotations. Poor mapping efficiency consequently causes gene loss, which can be dramatic especially for popular 3'-based single-cell technologies. Even for *Arabidopsis* scRNA-seq data sets, mapping rates vary but should be e.g. >85% (Table 1). A high percentage of reads mapping to intergenic regions (e.g. >20%) can be an indication either that not all genes are annotated or that the annotated 3' UTR regions should be longer. As a note, it is important to mention that the sequencing reads generated upon conducting a snRNA-seq experiment should be mapped against the exonic and intronic sequences of the annotated transcripts, a reflection of the capture of spliced and unspliced transcripts.

## Removal of low-quality cells/nuclei

After read mapping, low-quality cells/nuclei (e.g. cells/nuclei with low number of UMs and genes) need to be filtered out (see numerical recommendations in Table 1). Impairments in applying low-quality cell filtering may increase the noise in the data set and reduce the accuracy in downstream analysis, including cell clustering or erroneous identification of cell types (Fig. 4, A and B). General filter parameters will depend on the application and sample type (i.e. cells versus nuclei). For plants, it is good practice to exclude cells with high mitochondrial (e.g. >10%) and chloroplast (e.g. >5% to 10%) reads (Table 1). Such cases may indicate cells under stress because of perturbations during the sample preparation. These values might need to be adapted when studying highly or lowly metabolically active cells or cells undergoing e.g. programmed cell death. If no mitochondria or chloroplast gene is available, plotting the number of genes versus the total UMI count can be used instead to show cells with low data content. The barcode-rank-plot ("knee plot") is a commonly used tool to determine sample quality by ranking all barcodes according to their UMI content (Fig. 4A). A sample of low quality can be identified if there is no clearly defined boundary between barcodes with high UMI content and barcodes with significantly lower UMI content (Fig. 4C). Similarly, poor cluster separation in a sample that fails to differentiate cell types indicates an insufficient amount of transcript content per cell/nuclei.

Ambient RNA and the presence of empty droplets and doublets—both associated with droplet and nanowell-based technologies—can also lead to noise in the expression matrices, inaccurate cell clusters, and falsely differentially expressed genes. Therefore, it is important to ensure low amount of ambient RNA introduced during sample preparation (e.g. by mixing cells/nuclei from multiple species) and to optimize cell loading concentration (Fig. 4C).
A Ideal scenario

<figure><img src="image_5.png" /><figcaption>UMAP representation</figcaption></figure>

<figure><img src="image_6.png" /></figure>

B Poor filter parameters

<figure><img src="image_7.png" /><figcaption>technical variation replicates</figcaption></figure>

<figure><img src="image_8.png" /></figure>

* Apply stringent filter criteria

* Batch effect correction

* Remove repeats from analysis

C Low sample quality

<figure><img src="image_9.png" /><figcaption>low transcript content</figcaption></figure>

<figure><img src="image_10.png" /></figure>

* Optimize sample preparation

* Reduce protoplast or nuclei breakage and debris

* Optimize cell loading concentration

<figure><img src="image_11.png" /><figcaption>high doublet rate</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 4.** Overview of common problems in sc/snRNA-seq analysis and possible solutions.</figcaption></figure>

Bioinformatic tools (e.g. SoupX, Young and Behjati 2020; or CellBender, Fleming et al. 2023) can be used to computationally remove transcriptional noise introduced by ambient RNA. Experimentally, even though not yet shown suitable for plant single-cell transcriptomics, doublets can be detected using antibody (cell hashing; Stoeckius et al. 2018) or lipid-tagged indices (MULTI-Seq; McGinnis et al. 2019) sample multiplexing strategies. Here, doublets are identified if the cell-specific barcode is connected to multiple antibody or lipid-tagged indices, respectively. Mixing cells/nuclei from multiple species in 1 sample offers another experimental setup to identify doublets (Shulse et al. 2019; Fig. 4C). In this case, the doublet is identified due to the cell-specific barcode being linked to multiple species. While an experimental setup that allows doublet identification would be the best practice for identifying doublets, it will not help in removing doublets from already existing data sets. Even though significant advances have been made in identifying doublets computationally (see benchmark study; Xi and Li 2021), it remains a major challenge in general and even more in plant sc/sn transcriptomic analysis due to the presence of endoreduplication. Such polyploid cells will appear as outliers when plotting the gene content or UMI/cell, but additional expression quantification of ploidy marker genes (if available) allows to distinguish endoreduplicating cells (Wendrich et al. 2020).

## Identification of protoplast-induced genes

In order to exclude that the invasive enzymatic treatments needed to generate protoplasts might influence the observed transcriptional status of certain cells or cell populations, the overall transcriptomic responses induced during these procedures should be determined using bulk RNA-seq in each experimental setup (Birnbaum et al. 2005; Brady et al. 2007). Indeed, protoplast isolation adds a definite stressing factor to each cell type (Denyer et al. 2019; Xu et al. 2021; Wang, Huan, et al. 2021). Although the absence or presence of protoplast-induced genes did not alter cell clustering or annotation in *Arabidopsis* root scRNA-seq data (Denyer et al. 2019), these genes should at least be flagged in the data set to avoid misinterpretation.

## Challenge 6: cell cluster identification and annotation

### Data normalization, dimensionality reduction, and cell cluster visualization

Considering the limited number of transcripts per cell, single-cell transcriptomics fail to detect transcripts for most genes in a given cell. This sparsity is further enhanced by intrinsic noise from stochastic transcript fluctuations, cell-cycle state, and cell heterogeneity among other biological factors. As a
<figure><img src="image_7.png" /><figcaption>**B** Poor filter parameters</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**C** Low sample quality</figcaption></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

* Apply stringent filter criteria

* Batch effect correction

* Remove repeats from analysis

* Optimize sample preparation

* Reduce protoplast or nuclei breakage and debris

* Optimize cell loading concentration

**Figure 4.** Overview of common problems in sc/snRNA-seq analysis and possible solutions.

Bioinformatic tools (e.g. SoupX, Young and Behjati 2020; or CellBender, Fleming et al. 2023) can be used to computationally remove transcriptional noise introduced by ambient RNA. Experimentally, even though not yet shown suitable for plant single-cell transcriptomics, doublets can be detected using antibody (cell hashing; Stoeckius et al. 2018) or lipid-tagged indices (MULTI-Seq; McGinnis et al. 2019) sample multiplexing strategies. Here, doublets are identified if the cell-specific barcode is connected to multiple antibody or lipid-tagged indices, respectively. Mixing cells/nuclei from multiple species in 1 sample offers another experimental setup to identify doublets (Shulse et al. 2019; Fig. 4C). In this case, the doublet is identified due to the cell-specific barcode being linked to multiple species. While an experimental setup that allows doublet identification would be the best practice for identifying doublets, it will not help in removing doublets from already existing data sets. Even though significant advances have been made in identifying doublets computationally (see benchmark study; Xi and Li 2021), it remains a major challenge in general and even more in plant sc/sn transcriptomic analysis due to the presence of endoreduplication. Such polyploid cells will appear as outliers when plotting the gene content or UMI/cell, but additional expression quantification of ploidy marker genes (if available) allows to distinguish endoreduplicating cells (Wendrich et al. 2020).

## Identification of protoplast-induced genes

In order to exclude that the invasive enzymatic treatments needed to generate protoplasts might influence the observed transcriptional status of certain cells or cell populations, the overall transcriptomic responses induced during these procedures should be determined using bulk RNA-seq in each experimental setup (Birnbaum et al. 2005; Brady et al. 2007). Indeed, protoplast isolation adds a definite stressing factor to each cell type (Denyer et al. 2019; Xu et al. 2021; Wang, Huan, et al. 2021). Although the absence or presence of protoplast-induced genes did not alter cell clustering or annotation in *Arabidopsis* root scRNA-seq data (Denyer et al. 2019), these genes should at least be flagged in the data set to avoid misinterpretation.

## Challenge 6: cell cluster identification and annotation

### Data normalization, dimensionality reduction, and cell cluster visualization

Considering the limited number of transcripts per cell, single-cell transcriptomics fail to detect transcripts for most genes in a given cell. This sparsity is further enhanced by intrinsic noise from stochastic transcript fluctuations, cell-cycle state, and cell heterogeneity among other biological factors. As a
<figure><img src="image_13.png" /><figcaption>**Figure 5.** Cluster annotation and appropriate labels. A common technique to represent single-cell RNA-seq data requires mapping of each cell's transcriptome to a low (typically 2-) dimensional domain (e.g. tSNE or UMAP), after which highly similar cells group together into initially unannotated clusters. Annotating cells to a cell type or developmental stage is important for further interpretation of transcriptomic signatures. Traditional, manual annotation methods screen differentially expressed genes within each cluster for the presence of individual marker genes or by transferring knowledge on cluster annotation from a reference data set to an unannotated data set (reviewed in Clarke et al. 2021). Without prior knowledge of marker genes and reference data sets, automated annotation methods (e.g. eager, lazy, and marker learning methods; Xie et al. 2021) can aid in assigning labels based on comparing cell cluster-specific genes and their biological functions. Cells and clusters that cannot be annotated with high confidence to (only) one cell type must be analyzed carefully to determine if they represent a mix of subcell types and/or cell states. Although there is currently no standardized definition of a "cell type" or a "cell state" (discussed in more detail in Amini et al. 2023), it has been proposed that a cluster with homogenous marker gene expression among all cells likely represents a cell type, while expression gradients among the cells within a cluster represent cell states (Clarke et al. 2021). However, as these definitions are still evolving, caution should be taken to not mix different anatomical levels (cell types, tissues, or organs) and cellular processes (e.g. cell division and cell cycle) within the same visual representation. To avoid confusion, we recommend using multiple figures with different levels of labels (e.g. separate cell types, cell division states, tissues, and so on).</figcaption></figure>

result, it is necessary to implement scRNA-seq-specific normalization and batch correction protocols (see Luecken and Theis 2019 for a review on this specifically). Clusters of cells (i.e. cells that share similar expression profiles) are constructed using community detection algorithms, which control the degree to which similar cells should be grouped together or stay separate based on preset parameters. To visualize the data, "dimensional reduction" algorithms are applied to the data, typically via Principle Component Analysis (PCA), t-distributed Stochastic Neighbor Embedding (t-SNE), or Uniform Manifold Approximation and Projection (UMAP). However, a word of caution is necessary here: whereas PCA involves linear projections, t-SNE and UMAP are nonlinear transformations introducing significant distortions (Chari and Pachter 2023). Indeed, visualizing a synthetic data set with tSNE or UMAP revealed that cluster distances and locations in discrete and trajectory simulations are inaccurately represented compared to the defined distances in the original data (Wu et al. 2018). These results indicate that cluster distances or locations cannot be used alone to draw biological conclusions. If different clustering outcomes (produced with customized parameters) seem equally plausible, it is appropriate to apply the next step—mapping clusters to cell types—to each possibility and then use that extra information to decide which clustering makes most sense from a biological point of view.

## Cell cluster annotation

The annotation of clusters (meaning mapping each cluster to a cell type or state) is facilitated using manual or automated cell annotation methods (Fig. 5). Manual annotation requires

previous transcriptome knowledge gained from as many different cell types as possible. Cell-type marker genes for plant tissues in well-studied species can come from manually curated lists that are based on bulk RNA-seq data from purified cell populations (Jin et al. 2022) or from already annotated single-cell data sets. However, to date, databases containing this information are rare and restricted to a few plant species and tissues (e.g. Jin et al. 2022). One can appeal to interspecies correspondences to support the functional annotations of the clusters, assuming that a substantial fraction of orthologous genes share similar cell-type transcriptional specificities. For instance, if an atlas for another species is available, one can attempt to extend the orthology between genes to orthology between cell types. Such an approach has been very successful, leading to atlases for instance in rice (Wang, Huan, et al. 2021; Li et al. 2023), Medicago truncatula (Cervantes-Pérez et al. 2022), maize, sorghum, and Setaria (Guillotin et al. 2023). However, this approach also has its limitations: genes often exist in multigene families rendering the orthology mapping ambiguous, and the conservation of marker gene/cell type pairs is far from perfect (Movahedi et al. 2012). An additional limitation of these approaches is that they require multiple marker genes per cluster to ensure proper annotation.

Another option is so-called label transfer: transferring cluster labels between existing single-cell expression atlases (e.g. for roots: https://rootcellatlas.org/, or other organs: http://neomorph.salk.edu:9000/) to unlabeled data sets. Common tools used are scmap (Kiselev et al. 2018), SingleCellNet (Tan and Cahan 2019), SingleR (Aran et al. 2019), and Seurat (Stuart et al. 2019), which differ in their accuracy
<figure><img src="image_15.png" /><figcaption>**Figure 5.** Cluster annotation and appropriate labels. A common technique to represent single-cell RNA-seq data requires mapping of each cell's transcriptome to a low (typically 2-) dimensional domain (e.g. tSNE or UMAP), after which highly similar cells group together into initially unannotated clusters. Annotating cells to a cell type or developmental stage is important for further interpretation of transcriptomic signatures. Traditional, manual annotation methods screen differentially expressed genes within each cluster for the presence of individual marker genes or by transferring knowledge on cluster annotation from a reference data set to an unannotated data set (reviewed in Clarke et al. 2021). Without prior knowledge of marker genes and reference data sets, automated annotation methods (e.g. eager, lazy, and marker learning methods; Xie et al. 2021) can aid in assigning labels based on comparing cell cluster-specific genes and their biological functions. Cells and clusters that cannot be annotated with high confidence to (only) one cell type must be analyzed carefully to determine if they represent a mix of subcell types and/or cell states. Although there is currently no standardized definition of a "cell type" or a "cell state" (discussed in more detail in Amini et al. 2023), it has been proposed that a cluster with homogenous marker gene expression among all cells likely represents a cell type, while expression gradients among the cells within a cluster represent cell states (Clarke et al. 2021). However, as these definitions are still evolving, caution should be taken to not mix different anatomical levels (cell types, tissues, or organs) and cellular processes (e.g. cell division and cell cycle) within the same visual representation. To avoid confusion, we recommend using multiple figures with different levels of labels (e.g. separate cell types, cell division states, tissues, and so on).</figcaption></figure>

result, it is necessary to implement scRNA-seq-specific normalization and batch correction protocols (see Luecken and Theis 2019 for a review on this specifically). Clusters of cells (i.e. cells that share similar expression profiles) are constructed using community detection algorithms, which control the degree to which similar cells should be grouped together or stay separate based on preset parameters. To visualize the data, "dimensional reduction" algorithms are applied to the data, typically via Principle Component Analysis (PCA), t-distributed Stochastic Neighbor Embedding (t-SNE), or Uniform Manifold Approximation and Projection (UMAP). However, a word of caution is necessary here: whereas PCA involves linear projections, t-SNE and UMAP are nonlinear transformations introducing significant distortions (Chari and Pachter 2023). Indeed, visualizing a synthetic data set with tSNE or UMAP revealed that cluster distances and locations in discrete and trajectory simulations are inaccurately represented compared to the defined distances in the original data (Wu et al. 2018). These results indicate that cluster distances or locations cannot be used alone to draw biological conclusions. If different clustering outcomes (produced with customized parameters) seem equally plausible, it is appropriate to apply the next step—mapping clusters to cell types—to each possibility and then use that extra information to decide which clustering makes most sense from a biological point of view.

## Cell cluster annotation

The annotation of clusters (meaning mapping each cluster to a cell type or state) is facilitated using manual or automated cell annotation methods (Fig. 5). Manual annotation requires

previous transcriptome knowledge gained from as many different cell types as possible. Cell-type marker genes for plant tissues in well-studied species can come from manually curated lists that are based on bulk RNA-seq data from purified cell populations (Jin et al. 2022) or from already annotated single-cell data sets. However, to date, databases containing this information are rare and restricted to a few plant species and tissues (e.g. Jin et al. 2022). One can appeal to interspecies correspondences to support the functional annotations of the clusters, assuming that a substantial fraction of orthologous genes share similar cell-type transcriptional specificities. For instance, if an atlas for another species is available, one can attempt to extend the orthology between genes to orthology between cell types. Such an approach has been very successful, leading to atlases for instance in rice (Wang, Huan, et al. 2021; Li et al. 2023), Medicago truncatula (Cervantes-Pérez et al. 2022), maize, sorghum, and Setaria (Guillotin et al. 2023). However, this approach also has its limitations: genes often exist in multigene families rendering the orthology mapping ambiguous, and the conservation of marker gene/cell type pairs is far from perfect (Movahedi et al. 2012). An additional limitation of these approaches is that they require multiple marker genes per cluster to ensure proper annotation.

Another option is so-called label transfer: transferring cluster labels between existing single-cell expression atlases (e.g. for roots: https://rootcellatlas.org/, or other organs: http://neomorph.salk.edu:9000/) to unlabeled data sets. Common tools used are scmap (Kiselev et al. 2018), SingleCellNet (Tan and Cahan 2019), SingleR (Aran et al. 2019), and Seurat (Stuart et al. 2019), which differ in their accuracy
and their ability to handle sample- and protocol-related nested batch effect removal and the presence or absence of cell types in the reference atlas or target data set (Luecken et al. 2022). The advantage of this approach is that already published data can be directly reused. With both automated annotation approaches, it is important to consult several tools and select one final annotation, for instance using a majority rule.

As an alternative or complementary approach, clusters can be annotated manually. To achieve this, marker gene expression of cells can be visualized in a UMAP plot or via a dot plot showing cluster-specific expression of marker genes. Cell clusters with conflicting annotations or no annotation, due to e.g. low quality of the transcriptomic information or capture of uncharacterized cell (sub)type or cell transition state, should be marked as “unknown.”

## Refining cell cluster annotation

Ideally, each cluster will be clearly associated with one cell type, using any of the approaches described above. However, the current literature often mixes different levels of anatomical annotations. The ambiguity between classical anatomical descriptions and the new molecular characterizations of these cell types, tissues, and structures make it even more challenging to navigate among these definitions. Hence, to resolve these ambiguities, multiple hierarchically structured annotations can be used in which cells/clusters can be annotated according to e.g. broad expression domains, the tissue, or the cell type level (Michielsen et al. 2021). Until there is a consensus in the community, we recommend to separate the different levels of annotations (cell type, tissue, and cell cycle) into distinct plots. This will avoid ambiguity and misinterpretation of the data (Fig. 5).

## Validation of cell cluster annotation

Independent of the annotation approach used, we strongly recommend experimentally validating the main annotations by e.g. generating corresponding reporters, performing RNA in situ hybridizations, or complementing the data with spatial transcriptomics. This independent experimental approach is the only way to assess if differential gene expression observed in a high-throughput single-cell experiment is relevant in vivo or an artifact introduced by one of the many steps, e.g. protoplasting.

Beyond these more classical approaches, spatial transcriptomics can be used to support the annotation of clusters identified in plant sc/snRNA-seq data sets (Guillotin et al. 2023; Lee et al. 2023; Nobori et al. 2023). The technologies that allow probing tissue gene expression and simultaneously retaining its spatial location can be divided into two main categories: targeted and untargeted. The division into those two categories is based on the type of approach applied to analyze the tissue gene expression information. Specifically, targeted methods, i.e. in situ sequencing (Ke et al. 2013; Laureyns et al. 2021), MERFISH (Moffitt et al. 2016; Lee et al. 2023; Nobori et al. 2023), Xenium (Ke et al. 2013; Lee et al. 2015;

Janesick et al. 2022; Liu et al. 2022), NanoString CosMx (He et al. 2022), Molecular Cartography (Groiss et al. 2021; Guillotin et al. 2023; Yang et al. 2023) to list a few, require a priori knowledge on which genes to study, since these approaches use gene-specific probes to fluorescently visualize and count the gene transcripts of interest in the tissue. In contrast, untargeted methods, e.g. Visium (Giacomello et al. 2017; Liu et al. 2022, Liu et al. 2023; Peirats-Llobet et al. 2023), DBiT-Seq (Liu et al. 2020), Slide-seq v2/Curio (Stickels et al. 2021; Lee et al. 2023), and Stereo-seq (Xia et al. 2022), leverage a localized capture of polyadenylated transcripts, thus allowing to obtain 2D whole transcriptomic maps. Both targeted and untargeted approaches can aid in determining the spatial location of specific cell types or stages of interest (e.g. indeterminate and determinate SAM cells in maize, defined by their expression levels of PLASTOCHRON1; Laureyns et al. 2021), annotate cluster identities obtained in sc/snRNA-seq experiments, and/or validate marker genes (Guillotin et al. 2023; Lee et al. 2023; Nobori et al. 2023). As such, spatial transcriptomics and sc/snRNA-seq are complementary technologies for obtaining high-resolution spatiotemporal expression data.

## Challenge 7: application of trajectory inference in plant single-cell transcriptomics

Trajectory inference can predict developmental or stress-response trajectories in sc/snRNA-seq data sets, allowing one to pinpoint e.g. cell cycle transitions or bifurcations when new specific cell identities branch off from another lineage. While trajectory interference is a very promising tool to identify novel biological phenomena and questions, it is important to understand that the pseudotemporal ordering of cells along a trajectory is purely based on transcriptomic similarities, meaning that a sufficient, unbiased sampling of cells is required, as well as prior knowledge to assign a developmental direction (see Tritschler et al. 2019 for general recommendations). Furthermore, no conclusions about the spatial organization can be drawn, and it is advised to confirm and complement the trajectory output with other methods. Multiple trajectory interference methodologies have been developed, benchmarked, and used in the animal field (Saelens et al. 2019), but they show special promise in plants because plants have many continuous developmental programs and plant cells have remarkable capabilities for differentiation and adaptation. Trajectory analysis gave significant insight into recreating cell (type)-specific transcriptional events during plant developmental processes, such as lateral root development (Serrano-Ron et al. 2021), stomata development (Kim et al. 2023), root hair (Denyer et al. 2019; Shulse et al. 2019), pistil (Li et al. 2023), or phloem development (Roszak et al. 2021; Otero et al. 2022). For example, when combined with live imaging, trajectory-predicted expression gradients and cell-type-specific transcriptional networks allowed a complete reconstruction of the developmental process and a precise, cell-by-cell lineage tracing during protophloem development (Roszak et al.
and their ability to handle sample- and protocol-related nested batch effect removal and the presence or absence of cell types in the reference atlas or target data set (Luecken et al. 2022). The advantage of this approach is that already published data can be directly reused. With both automated annotation approaches, it is important to consult several tools and select one final annotation, for instance using a majority rule.

As an alternative or complementary approach, clusters can be annotated manually. To achieve this, marker gene expression of cells can be visualized in a UMAP plot or via a dot plot showing cluster-specific expression of marker genes. Cell clusters with conflicting annotations or no annotation, due to e.g. low quality of the transcriptomic information or capture of uncharacterized cell (sub)type or cell transition state, should be marked as “unknown.”

## Refining cell cluster annotation

Ideally, each cluster will be clearly associated with one cell type, using any of the approaches described above. However, the current literature often mixes different levels of anatomical annotations. The ambiguity between classical anatomical descriptions and the new molecular characterizations of these cell types, tissues, and structures make it even more challenging to navigate among these definitions. Hence, to resolve these ambiguities, multiple hierarchically structured annotations can be used in which cells/clusters can be annotated according to e.g. broad expression domains, the tissue, or the cell type level (Michielsen et al. 2021). Until there is a consensus in the community, we recommend to separate the different levels of annotations (cell type, tissue, and cell cycle) into distinct plots. This will avoid ambiguity and misinterpretation of the data (Fig. 5).

## Validation of cell cluster annotation

Independent of the annotation approach used, we strongly recommend experimentally validating the main annotations by e.g. generating corresponding reporters, performing RNA in situ hybridizations, or complementing the data with spatial transcriptomics. This independent experimental approach is the only way to assess if differential gene expression observed in a high-throughput single-cell experiment is relevant in vivo or an artifact introduced by one of the many steps, e.g. protoplasting.

Beyond these more classical approaches, spatial transcriptomics can be used to support the annotation of clusters identified in plant sc/snRNA-seq data sets (Guillotin et al. 2023; Lee et al. 2023; Nobori et al. 2023). The technologies that allow probing tissue gene expression and simultaneously retaining its spatial location can be divided into two main categories: targeted and untargeted. The division into those two categories is based on the type of approach applied to analyze the tissue gene expression information. Specifically, targeted methods, i.e. in situ sequencing (Ke et al. 2013; Laureyns et al. 2021), MERFISH (Moffitt et al. 2016; Lee et al. 2023; Nobori et al. 2023), Xenium (Ke et al. 2013; Lee et al. 2015;

Janesick et al. 2022; Liu et al. 2022), NanoString CosMx (He et al. 2022), Molecular Cartography (Groiss et al. 2021; Guillotin et al. 2023; Yang et al. 2023) to list a few, require a priori knowledge on which genes to study, since these approaches use gene-specific probes to fluorescently visualize and count the gene transcripts of interest in the tissue. In contrast, untargeted methods, e.g. Visium (Giacomello et al. 2017; Liu et al. 2022, Liu et al. 2023; Peirats-Llobet et al. 2023), DBiT-Seq (Liu et al. 2020), Slide-seq v2/Curio (Stickels et al. 2021; Lee et al. 2023), and Stereo-seq (Xia et al. 2022), leverage a localized capture of polyadenylated transcripts, thus allowing to obtain 2D whole transcriptomic maps. Both targeted and untargeted approaches can aid in determining the spatial location of specific cell types or stages of interest (e.g. indeterminate and determinate SAM cells in maize, defined by their expression levels of PLASTOCHRON1; Laureyns et al. 2021), annotate cluster identities obtained in sc/snRNA-seq experiments, and/or validate marker genes (Guillotin et al. 2023; Lee et al. 2023; Nobori et al. 2023). As such, spatial transcriptomics and sc/snRNA-seq are complementary technologies for obtaining high-resolution spatiotemporal expression data.

## Challenge 7: application of trajectory inference in plant single-cell transcriptomics

Trajectory inference can predict developmental or stress-response trajectories in sc/snRNA-seq data sets, allowing one to pinpoint e.g. cell cycle transitions or bifurcations when new specific cell identities branch off from another lineage. While trajectory interference is a very promising tool to identify novel biological phenomena and questions, it is important to understand that the pseudotemporal ordering of cells along a trajectory is purely based on transcriptomic similarities, meaning that a sufficient, unbiased sampling of cells is required, as well as prior knowledge to assign a developmental direction (see Tritschler et al. 2019 for general recommendations). Furthermore, no conclusions about the spatial organization can be drawn, and it is advised to confirm and complement the trajectory output with other methods. Multiple trajectory interference methodologies have been developed, benchmarked, and used in the animal field (Saelens et al. 2019), but they show special promise in plants because plants have many continuous developmental programs and plant cells have remarkable capabilities for differentiation and adaptation. Trajectory analysis gave significant insight into recreating cell (type)-specific transcriptional events during plant developmental processes, such as lateral root development (Serrano-Ron et al. 2021), stomata development (Kim et al. 2023), root hair (Denyer et al. 2019; Shulse et al. 2019), pistil (Li et al. 2023), or phloem development (Roszak et al. 2021; Otero et al. 2022). For example, when combined with live imaging, trajectory-predicted expression gradients and cell-type-specific transcriptional networks allowed a complete reconstruction of the developmental process and a precise, cell-by-cell lineage tracing during protophloem development (Roszak et al.
2021). An interesting computational analysis that can be applied to verify predicted trajectories in plant tissues is making use of the ploidy increase in plant cells as they mature from meristematic cells with 2C/4C content into a more differentiated stage, marked by 8C/16C or higher content. In *Arabidopsis* root tissues, these ploidy states have been linked to specific markers, leading to predictions of the ploidy status of each individual cell. This increase in ploidy level allows to pinpoint or validate the more meristematic cells with low ploidy levels and the more differentiated cells in a sc/snRNA-seq data set with higher ploidy levels (Bhosale et al. 2018; Wendrich et al. 2020; Shahan et al. 2022).

## Challenge 8: documentation and publication of plant single-cell data sets

Recent advances in single-cell omic technologies in plants have enabled insights into diverse aspects of physiology and development and have been the centerpiece of a growing number of elegant studies. However, there is additional potential for single-cell resources through their reuse in integrative meta-analyses (showcased e.g. in Leote et al. 2022). The multiple applications would increase their power and depth through greater numbers of cells, a more comprehensive assessment of biological variation, and enhanced enrichment for different cell types or states that are targeted by individual studies. Popular methods enabled within software packages (e.g. Seurat, Monocle, Scanpy, and Harmony) (Wolf et al. 2018; Korsunsky et al. 2019; Van den Berge et al. 2020; Hao et al. 2021) have streamlined the process of stitching data sets together across samples, studies, and experimental platforms, even across tissues and species, allowing for a more expansive use of single-cell data for describing and understanding biological organization across scales.

To make optimal use of the generated data in the plant field, there is an urgent need for published single-cell data and their associated metadata to be more easily accessible and usable. First steps toward establishing a suitable infrastructure that allows data storage and comparison have been taken within the framework of the Plant Cell Atlas (Fahlgren et al. 2023). To a large extent, this FAIR (Findable, Accessible, Interoperable, Reproducible) principle has not been an issue for raw single-cell data. Unprocessed FASTQ files are deposited routinely through major data portals, such as the National Center for Biotechnology Information Short Read Archive (NCBI SRA; http://www.ncbi.nlm.nih.gov/sra), and accessibility to raw data is usually mandated by journals, funding agencies, and institutions. In most cases, well-indexed raw data are available and cited in journals. However, care must be taken for single-cell data generated through popular platforms in that only reads from the paired-end sequencing strategy, as well as indexed and/or UMI FASTQ, must be deposited. These genomic reads alone are insufficient for reconstructing the single-cell count matrix necessary for nearly all analysis steps. Processed data should also be stored in publicly accessible repositories such as the

NCBI Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo). Processed data not only include the cell/gene counts matrix at minimum but can also include more complex data objects such as those generated by the Seurat (R) (Hao et al. 2021) or AnnData (Python) (Wolf et al. 2018) packages. The EBI Single-cell Expression Atlas is another repository that will accept processed count matrices and allows for the deposition of cell metadata (see below).

Less prevalent is the public accessibility of metadata associated with single-cell studies. This metadata takes several forms: (i) experimental metadata describing how samples were treated, what tissues were harvested, how the samples were processed, and what platforms/versions single-cell sequencing was performed on them; and (ii) imputed metadata describing attributes of cells defined by downstream analysis of the raw data, including the cell's identifier (usually a sequence barcode), the number of distinct genes detected, total transcripts detected, assigned cell type or cluster number, or other features that are described in a study but not immediately available from the raw data. Experimental metadata, including protocols used to generate samples, should be well documented in the manuscript that presents them, and we encourage depositing the protocol in a public repository. Protocols.io is emerging as a standard for this (see e.g. the repository for the Plant Cell Atlas; https://www.protocols.io/workspaces/plant-cell-atlas), and detailed sample information should go along with the raw data when submitting to a repository like GEO/SRA. Imputed metadata can be provided along with the raw data in the NCBI GEO as a simple machine-readable table (.csv or .tsv). All forms of metadata are vitally important for integrative analysis, as it is difficult, if not impossible, to recreate exactly the analysis steps performed from a published study to replicate results. In Table 1, we have listed important metadata parameters that should be accompanying all publications using sc/snRNA-seq data to ensure transparency in data quality and will allow for a better interpretation of the results and their use in larger meta-analyses. Although it is challenging to give an exact number for these parameters due to the vast differences in experimental systems, species, tissues, and technologies, we do attempt to provide a range based on the collective author's experience in processing and analyzing multiple species and organs to guide the less experienced user. An empty version of this table can be downloaded as Supplementary Table S1 and freely used for publication. All coauthors collectively commit to start using this table in all publications.

Finally, analysis scripts, software environments, and (ideally) visualization portals should be made publicly available on established repositories/portals. Where possible, a well-documented code can be stored in dedicated repositories. Analysis environments should at minimum be well documented (with versions specified for all software packages used for analysis).

Many of these tools are only valid for single-cell transcriptomic or chromatin accessibility methods that use DNA sequence as an output and are thus incompatible with other data modalities such as proteomics or metabolomics.
2021). An interesting computational analysis that can be ap-
plied to verify predicted trajectories in plant tissues is making
use of the ploidy increase in plant cells as they mature from
meristematic cells with 2C/4C content into a more differen-
tiated stage, marked by 8C/16C or higher content. In
Arabidopsis root tissues, these ploidy states have been linked
to specific markers, leading to predictions of the ploidy status
of each individual cell. This increase in ploidy level allows to
pinpoint or validate the more meristematic cells with low
ploidy levels and the more differentiated cells in a sc/
snRNA-seq data set with higher ploidy levels (Bhosale et al.
2018; Wendrich et al. 2020; Shahan et al. 2022).

Challenge 8: documentation and publication of plant single-cell data sets

Recent advances in single-cell omic technologies in plants have enabled insights into diverse aspects of physiology and development and have been the centerpiece of a growing number of elegant studies. However, there is additional potential for single-cell resources through their reuse in integrative meta-analyses (showcased e.g. in Leote et al. 2022). The multiple applications would increase their power and depth through greater numbers of cells, a more comprehensive assessment of biological variation, and enhanced enrichment for different cell types or states that are targeted by individual studies. Popular methods enabled within software packages (e.g. Seurat, Monocle, Scanpy, and Harmony) (Wolf et al. 2018; Korsunsky et al. 2019; Van den Berge et al. 2020; Hao et al. 2021) have streamlined the process of stitching data sets together across samples, studies, and experimental platforms, even across tissues and species, allowing for a more expansive use of single-cell data for describing and understanding biological organization across scales.

To make optimal use of the generated data in the plant
field, there is an urgent need for published single-cell data
and their associated metadata to be more easily accessible
and usable. First steps toward establishing a suitable infra-
structure that allows data storage and comparison have
been taken within the framework of the Plant Cell Atlas
(Fahlgren et al. 2023). To a large extent, this FAIR (Findable,
Accessible, Interoperable, Reproducible) principle has not
been an issue for raw single-cell data. Unprocessed FASTQ
files are deposited routinely through major data portals,
such as the National Center for Biotechnology Information
Short Read Archive (NCBI SRA; http://www.ncbi.nlm.nih.
gov/sra), and accessibility to raw data is usually mandated
by journals, funding agencies, and institutions. In most cases,
well-indexed raw data are available and cited in journals.
However, care must be taken for single-cell data generated
through popular platforms in that only reads from the
paired-end sequencing strategy, as well as indexed and/or
UMI FASTQ, must be deposited. These genomic reads alone
are insufficient for reconstructing the single-cell count matrix
necessary for nearly all analysis steps. Processed data should
also be stored in publicly accessible repositories such as the

NCBI Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo). Processed data not only include the cell/gene counts matrix at minimum but can also include more complex data objects such as those generated by the Seurat (R) (Hao et al. 2021) or AnnData (Python) (Wolf et al. 2018) packages. The EBI Single-cell Expression Atlas is another repository that will accept processed count matrices and allows for the deposition of cell metadata (see below).

Less prevalent is the public accessibility of metadata asso-
ciated with single-cell studies. This metadata takes several
forms: (i) experimental metadata describing how samples
were treated, what tissues were harvested, how the samples
were processed, and what platforms/versions single-cell sequen-
cing was performed on them; and (ii) imputed metadata de-
scribing attributes of cells defined by downstream analysis of
the raw data, including the cell's identifier (usually a sequence
barcode), the number of distinct genes detected, total tran-
scripts detected, assigned cell type or cluster number, or other
features that are described in a study but not immediately avail-
able from the raw data. Experimental metadata, including pro-
tocols used to generate samples, should be well documented in
the manuscript that presents them, and we encourage deposit-
ing the protocol in a public repository. Protocols.io is emerging
as a standard for this (see e.g. the repository for the Plant Cell
Atlas; https://www.protocols.io/workspaces/plant-cell-atlas),
and detailed sample information should go along with the
raw data when submitting to a repository like GEO/SRA.
Imputed metadata can be provided along with the raw data
in the NCBI GEO as a simple machine-readable table (.csv
or .tsv). All forms of metadata are vitally important for integra-
tive analysis, as it is difficult, if not impossible, to recreate exactly
the analysis steps performed from a published study to replicate
results. In Table 1, we have listed important metadata para-
meters that should be accompanying all publications using
sc/snRNA-seq data to ensure transparency in data quality and
will allow for a better interpretation of the results and their
use in larger meta-analyses. Although it is challenging to give
an exact number for these parameters due to the vast differ-
ences in experimental systems, species, tissues, and technologies,
we do attempt to provide a range based on the collective
author's experience in processing and analyzing multiple species
and organs to guide the less experienced user. An empty version
of this table can be downloaded as Supplementary Table S1 and
freely used for publication. All coauthors collectively commit to
start using this table in all publications.

Finally, analysis scripts, software environments, and (ideally) visualization portals should be made publicly available on established repositories/portals. Where possible, a well-documented code can be stored in dedicated repositories. Analysis environments should at minimum be well documented (with versions specified for all software packages used for analysis).

Many of these tools are only valid for single-cell transcrip-
tomic or chromatin accessibility methods that use DNA se-
quence as an output and are thus incompatible with other
data modalities such as proteomics or metabolomics.
While these data types are outside the scope of this review, they too should adhere to FAIR principles (Wilkinson et al. 2016) to ensure that they can be integrated with other data sets when the computational infrastructure develops. These data management tools are also not guaranteed to be future proof. New technologies may arise that will render existing data sets/tools obsolete. Thus, it is imperative that, as the field evolves, we do not forget about legacy data and ensure that it is preserved in a way that will be useful before this becomes impossible.

## Conclusion and outlook

Even for more advanced users, it is challenging to keep up with the rapidly evolving field of sc/snRNA-seq. Continuous advancements are not limited to the actual technology and also include the choice of which technology to use for a specific problem, which method to use to isolate high-quality cells or nuclei, and how to analyze, compare, and store these vast amounts of data. In this commentary, we discuss recommendations regarding data generation, analysis, storage, and documentation to ensure transparency in publication and optimal use of the generated data across experiments, tissues, species, and laboratories. We collectively commit to following the guidelines and recommendations (Table 1). Future research and method maturation will allow us to fine-tune and expand these guidelines and recommendations.

## Acknowledgments

We thank the Plant Cell Atlas Research Coordination Network funded by the United States National Science Foundation grant MCB-2052590 for fostering environments for scientists to discuss, network, and collaborate on the emerging field of plant single-cell omics. We also thank and acknowledge helpful input and discussion from the following endorsers of this manuscript (in random order): Dominique Bergmann, David Jackson, Xiaosa Xu, Ronan O'Malley, Idan Efroni, Robert Schmitz, Victoria Mironova, Klaas Vandepoele, Alexis Maizel, Jian Xu, Marja Timmermans, Thomas Denyer, Sue Rhee, Zhe Yan, Ykä Helariutta, Joshua Cuperus, Siobhan Brady, Joe Ecker, and Ken Birnbaum.

## Author contributions

C.S., T.E., J.L., H.N., K.K., M.T.R., T.G., M.L., and B.D.R. conceived this initiative. C.S., T.E., D.S., M.N., L.S.B., Y.K., R.S., H.N., J.L., S.G., O.C.M., B.C., K.K., M.T.R., G.P., T.G., M.L., and B.D.R. wrote a first draft version that was commented on and finished by all coauthors.

## Supplemental data

**Supplementary Table S1.** Necessary reported information to allow evaluation and repetition of a plant single cell/nucleus experiment.

## Funding

This work was supported by the European Research Council (ERC StG TORPEDO: 714055 and ERC CoG PIPELINES: 101043257 to B.D.R. and T.E.; ERC CoG PLANTSTEMS: 647148 to T.G. and ERC SyG: 810296 DECODE to J.U.L.); the Flemish Government VLAIO fund (HBC.2019.2917) to Y.K.; VIB TechWatch Funding to C.S.; EMBO YIP to B.D.R.; the Nebraska Soybean Board (award #708); the National Science Foundation (awards #2127485 and #1854326); the USDA-NIFA (award #2022-67012-36144); the Nebraska Research Initiative to M.L.; the DFG FOR5235 "Cereal Stem Cell Systems" to M.T.R.; and the US National Institutes of Health (NRSA postdoctoral fellowship 1F32GM136030) to R.S.; the work conducted by the U.S. Department of Energy Joint Genome Institute (https://ror.org/04xm1d337), a DOE Office of Science User Facility, is supported by the Office of Science of the U.S. Department of Energy operated under Contract No. DE-AC02-05CH11231. The work of O.C.M. has benefited from a French State grant (Saclay Plant Sciences, reference no ANR-17-EUR-0007, EUR SPS-GSR). J.W.W. is supported by grants from the National Natural Science Foundation of China (31788103; 31721001), Strategic Priority Research Program of the Chinese Academy of Sciences (XDB27030101). D.S. was supported by the Japan Science and Technology Agency (JPMJPR2046) and Deutsche Forschungsgemeinschaft (504847536). K.K. is supported by DFG projects 438774542, 458750707, and 270050988.

*Conflict of interest statement.* None declared.

## References

* Abramson BW, Novotny M, Hartwick NT, Colt K, Aevermann BD, Scheuermann RH, Michael TP. The genome and preliminary single-nuclei transcriptome of *Lemna minuta* reveals mechanisms of invasiveness. Plant Physiol. 2022;188(2):879-897. https://doi.org/10.1093/plphys/kiab564

* Amini S, Doyle JJ, Libault M. The evolving definition of plant cell type. Front Plant Sci. 2023;14:1271070. https://doi.org/10.3389/fpls.2023.1271070

* Apelt F, Mavrothalassiti E, Gupta S, Machin F, Olas JJ, Annunziata MG, Schindelasch D, Kragler F. Shoot and root single cell sequencing reveals tissue- and daytime-specific transcriptome profiles. Plant Physiol. 2022;188(2):861-878. https://doi.org/10.1093/plphys/kiab537

* Aran D, Looney AP, Liu L, Wu E, Fong V, Hsu A, Chak S, Naikawadi RP, Wolters PJ, Abate AR, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. 2019;20(2):163-172. https://doi.org/10.1038/s41590-018-0276-y

* Attar M, Sharma E, Li S, Bryer C, Cubitt L, Broxholme J, Lockstone H, Kinchen J, Simmons A, et al. A practical solution for preserving single cells for RNA sequencing. Sci Rep. 2018;8(1):2151-2161. https://doi.org/10.1038/s41598-018-20372-7

* Bezrutczyk M, Zöllner NR, Kruse CPS, Hartwig T, Lautwein T, Köhrer K, Frommer WB, Kim J-Y. Evidence for phloem loading via the abaxial bundle sheath cells in maize leaves. Plant Cell. 2021;33(3):531-547. https://doi.org/10.1093/plcell/koaa055

* Bhosale R, Boudolf V, Cuevas F, Lu R, Eekhout T, Hu Z, Van Isterdael G, Lambert GM, Xu F, Nowack MK, et al. A spatiotemporal DNA endoploidy map of the *Arabidopsis* root reveals roles for the
While these data types are outside the scope of this review, they too should adhere to FAIR principles (Wilkinson et al. 2016) to ensure that they can be integrated with other data sets when the computational infrastructure develops. These data management tools are also not guaranteed to be future proof. New technologies may arise that will render existing data sets/tools obsolete. Thus, it is imperative that, as the field evolves, we do not forget about legacy data and ensure that it is preserved in a way that will be useful before this becomes impossible.

## Conclusion and outlook

Even for more advanced users, it is challenging to keep up with the rapidly evolving field of sc/snRNA-seq. Continuous advancements are not limited to the actual technology and also include the choice of which technology to use for a specific problem, which method to use to isolate high-quality cells or nuclei, and how to analyze, compare, and store these vast amounts of data. In this commentary, we discuss recommendations regarding data generation, analysis, storage, and documentation to ensure transparency in publication and optimal use of the generated data across experiments, tissues, species, and laboratories. We collectively commit to following the guidelines and recommendations (Table 1). Future research and method maturation will allow us to fine-tune and expand these guidelines and recommendations.

## Acknowledgments

We thank the Plant Cell Atlas Research Coordination Network funded by the United States National Science Foundation grant MCB-2052590 for fostering environments for scientists to discuss, network, and collaborate on the emerging field of plant single-cell omics. We also thank and acknowledge helpful input and discussion from the following endorsers of this manuscript (in random order): Dominique Bergmann, David Jackson, Xiaosa Xu, Ronan O'Malley, Idan Efroni, Robert Schmitz, Victoria Mironova, Klaas Vandepoele, Alexis Maizel, Jian Xu, Marja Timmermans, Thomas Denyer, Sue Rhee, Zhe Yan, Ykä Helariutta, Joshua Cuperus, Siobhan Brady, Joe Ecker, and Ken Birnbaum.

## Author contributions

C.S., T.E., J.L., H.N., K.K., M.T.R., T.G., M.L., and B.D.R. conceived this initiative. C.S., T.E., D.S., M.N., L.S.B., Y.K., R.S., H.N., J.L., S.G., O.C.M., B.C., K.K., M.T.R., G.P., T.G., M.L., and B.D.R. wrote a first draft version that was commented on and finished by all coauthors.

## Supplemental data

**Supplementary Table S1.** Necessary reported information to allow evaluation and repetition of a plant single cell/nucleus experiment.

## Funding

This work was supported by the European Research Council (ERC StG TORPEDO: 714055 and ERC CoG PIPELINES: 101043257 to B.D.R. and T.E.; ERC CoG PLANTSTEMS: 647148 to T.G. and ERC SyG: 810296 DECODE to J.U.L.); the Flemish Government VLAIO fund (HBC.2019.2917) to Y.K.; VIB TechWatch Funding to C.S.; EMBO YIP to B.D.R.; the Nebraska Soybean Board (award #708); the National Science Foundation (awards #2127485 and #1854326); the USDA-NIFA (award #2022-67012-36144); the Nebraska Research Initiative to M.L.; the DFG FOR5235 "Cereal Stem Cell Systems" to M.T.R.; and the US National Institutes of Health (NRSA postdoctoral fellowship 1F32GM136030) to R.S.; the work conducted by the U.S. Department of Energy Joint Genome Institute (https://ror.org/04xm1d337), a DOE Office of Science User Facility, is supported by the Office of Science of the U.S. Department of Energy operated under Contract No. DE-AC02-05CH11231. The work of O.C.M. has benefited from a French State grant (Saclay Plant Sciences, reference no ANR-17-EUR-0007, EUR SPS-GSR). J.W.W. is supported by grants from the National Natural Science Foundation of China (31788103; 31721001), Strategic Priority Research Program of the Chinese Academy of Sciences (XDB27030101). D.S. was supported by the Japan Science and Technology Agency (JPMJPR2046) and Deutsche Forschungsgemeinschaft (504847536). K.K. is supported by DFG projects 438774542, 458750707, and 270050988.

*Conflict of interest statement.* None declared.

## References

* Abramson BW, Novotny M, Hartwick NT, Colt K, Aevermann BD, Scheuermann RH, Michael TP. The genome and preliminary single-nuclei transcriptome of *Lemna minuta* reveals mechanisms of invasiveness. Plant Physiol. 2022;188(2):879-897. https://doi.org/10.1093/plphys/kiab564

* Amini S, Doyle JJ, Libault M. The evolving definition of plant cell type. Front Plant Sci. 2023;14:1271070. https://doi.org/10.3389/fpls.2023.1271070

* Apelt F, Mavrothalassiti E, Gupta S, Machin F, Olas JJ, Annunziata MG, Schindelasch D, Kragler F. Shoot and root single cell sequencing reveals tissue- and daytime-specific transcriptome profiles. Plant Physiol. 2022;188(2):861-878. https://doi.org/10.1093/plphys/kiab537

* Aran D, Looney AP, Liu L, Wu E, Fong V, Hsu A, Chak S, Naikawadi RP, Wolters PJ, Abate AR, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. 2019;20(2):163-172. https://doi.org/10.1038/s41590-018-0276-y

* Attar M, Sharma E, Li S, Bryer C, Cubitt L, Broxholme J, Lockstone H, Kinchen J, Simmons A, et al. A practical solution for preserving single cells for RNA sequencing. Sci Rep. 2018;8(1):2151-2161. https://doi.org/10.1038/s41598-018-20372-7

* Bezrutczyk M, Zöllner NR, Kruse CPS, Hartwig T, Lautwein T, Köhrer K, Frommer WB, Kim J-Y. Evidence for phloem loading via the abaxial bundle sheath cells in maize leaves. Plant Cell. 2021;33(3):531-547. https://doi.org/10.1093/plcell/koaa055

* Bhosale R, Boudolf V, Cuevas F, Lu R, Eekhout T, Hu Z, Van Isterdael G, Lambert GM, Xu F, Nowack MK, et al. A spatiotemporal DNA endoploidy map of the *Arabidopsis* root reveals roles for the
end-cycle in root development and stress adaptation. Plant Cell.
2018:30(10):2330-2351. https://doi.org/10.1105/tpc.17.00983

Birnbaum K, Jung JW, Wang JY, Lambert GM, Hirst JA, Galbraith DW, Benfey PN. Cell type-specific expression profiling in plants via cell sorting of protoplasts from fluorescent reporter lines. Nat Methods. 2005;2(8):615–619. https://doi.org/10.1038/ nmeth0805-615

Booeshaghi AS, Chen X, Pachter L. A machine-readable specification for genomics assays. bioRxiv 533215. https://doi.org/10.1101/2023.03. 17.533215, 21 March 2023, preprint: not peer reviewed.

Brady SM, Orlando DA, Lee J-Y, Wang JY, Koch J, Dinneny JR, Mace D, Ohler U, Benfey PN. A high-resolution root spatiotemporal map reveals dominant expression patterns. Science. 2007:318(5851): 801–806. https://doi.org/10.1126/science.1146265

Cao J, Packer JS, Ramani V, Cusanovich DA, Huynh C, Daza R, Qiu X, Lee C, Furlan SN, Steemers FJ, et al. Comprehensive single-cell tran- scriptional profiling of a multicellular organism. Science. 2017:357(6352):661–667. https://doi.org/10.1126/science.aam8940

Cartwright DA, Brady SM, Orlando DA, Sturmfels B, Benfey PN. Reconstructing spatiotemporal gene expression data from partial ob- servations. Bioinformatics. 2009:25(19):2581–2587. https://doi.org/ 10.1093/bioinformatics/btp437

Cervantes-Pérez SA, Thibivilliers S, Laffont C, Farmer AD, Frugier F, Libault M. Cell-specific pathways recruited for symbiotic nodulation in the *Medicago truncatula* legume. Mol Plant. 2022:15(12):1868–1888. https://doi.org/10.1016/j.molp.2022.10.021

Chari T, Pachter L. The specious art of single-cell genomics. PLoS Comput Biol. 2023:19(8):e1011288. https://doi.org/10.1371/journal. pcbi.1011288

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

Clarke ZA, Andrews TS, Atif J, Pouyabahar D, Innes BT, MacParland SA, Bader GD. Tutorial: guidelines for annotating single-cell tran- scriptomic maps using automated and manual methods. Nat Protoc. 2021:16(6):2749–2764. https://doi.org/10.1038/s41596-021- 00534-0

Conde D, Triozzi PM, Pereira WJ, Schmidt HW, Balmant KM, Knaack SA, Redondo-López A, Roy S, Dervinis C, Kirst M. Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in *Populus*. Development. 2022:149(21):dev200632. https://doi.org/10.1242/dev.200632

Denyer T, Ma X, Klesen S, Scacchi E, Nieselt K, Timmermans MCP. Spatiotemporal developmental trajectories in the *Arabidopsis* root revealed using high-throughput single-cell RNA sequencing. Dev Cell. 2019:48(6):840–852.e5. https://doi.org/10.1016/j.devcel. 2019.02.022

Efroni I, Mello A, Nawy T, Ip P-L, Rahni R, DelRose N, Powers A, Satija R, Birnbaum KD. Root regeneration triggers an embryo-like sequence guided by hormonal interactions. Cell. 2016:165(7): 1721–1733. https://doi.org/10.1016/j.cell.2016.04.046

Fahlgren N, Kapoor M, Yordanova G, Papatheodorou I, Waese J, Cole B, Harrison P, Ware D, Tickle T, Paten B, et al. Toward a data infrastructure for the Plant Cell Atlas. Plant Physiol. 2023:191(1):35–46. https://doi.org/10.1093/plphys/kiac468

Farmer A, Thibivilliers S, Ryu KH, Schiefelbein J, Libault M. 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. 2021:14(3):372–383. https://doi.org/ 10.1016/j.molp.2021.01.001

Ferrari C, Pérez NM, Vandepoele K. MINI-EX: integrative inference of single-cell gene regulatory networks in plants. Mol Plant. 2022:15(11): 1807–1824. https://doi.org/10.1016/j.molp.2022.10.016

Fischer S, Gillis J. How many markers are needed to robustly determine a cell's type? iScience. 2021:24(11):103292. https://doi.org/10.1016/j.isci.2021.103292

Fleming SJ, Chaffin MD, Arduini A, Akkad A-D, Banks E, Marioni JC, Philippakis AA, Ellinor PT, Babadi M. Unsupervised removal of systematic background noise from droplet-based single-cell experiments using CellBender. Nat Methods. 2023:20(9):1323–1335. https://doi.org/10.1038/s41592-023-01943-7

Fox DT, Soltis DE, Soltis PS, Ashman T-L, Van de Peer Y. Polyploidy: a biological force from cells to ecosystems. Trends Cell Biol. 2020:30(9): 688–694. https://doi.org/10.1016/j.tcb.2020.06.006

Gala HP, Lancot A, Jean-Baptiste K, Guiziou S, Chu JC, Zemke JE, George W, Queitsch C, Cuperus JT, Nemhauser JL. A single-cell view of the transcriptome during lateral root initiation in *Arabidopsis thaliana*. Plant Cell. 2021:33(7):2197–2220. https://doi.org/10.1093/plcell/koab101

Giacomello S, Salmén F, Terebieniec BK, Vickovic S, Navarro JF, Alexeyenko A, Reimegård J, McKee LS, Mannapperuma C, Bulone V, et al. Spatially resolved transcriptome profiling in model plant species. Nat Plants. 2017:3(6):1–11. https://doi.org/10.1038/ nplants.2017.61

Graeff M, Rana S, Wendrich JR, Dorier J, Eekhout T, Aliaga Fandino AC, Guex N, Bassel GW, De Rybel B, Hardtke CS. A single-cell morpho-transcriptomic map of brassinosteroid action in the *Arabidopsis* root. Mol Plant. 2021:14(12):1985–1999. https://doi.org/10.1016/j.molp.2021.07.021

Groiss S, Pabst D, Faber C, Meier A, Bogdoll A, Unger C, Nilges B, Strauss S, Föderl-Höbenreich E, Hardt M, et al. Highly resolved spa- tial transcriptomics for detection of rare events in cell. bioRxiv 463936. https://doi.org/10.1101/2021.10.11.463936, 12 October 2021, preprint: not peer reviewed.

Guillotin B, Rahni R, Passalacqua M, Mohammed MA, Xu X, Raju SK, Ramírez CO, Jackson D, Groen SC, Gillis J, et al. A pan-grass transcriptome reveals patterns of cellular divergence in crops. Nature. 2023:617(7962):785–791. https://doi.org/10.1038/s41586- 023-06053-0

Hao Y, Hao S, Andersen-Nissen E, Mauck WM III, Zheng S, Butler A, Lee MJ, Wilk AJ, Darby C, Zager M, et al. Integrated analysis of multimodal single-cell data. Cell. 2021:184(13):3573–3587.e29. https://doi.org/10.1016/j.cell.2021.04.048

He S, Bhatt R, Brown C, Brown EA, Buhr DL, Chantranuvatana K, Danaher P, Dunaway D, Garrison RG, Geiss G, et al. High-plex im- aging of RNA and proteins at subcellular resolution in fixed tissue by spatial molecular imaging. Nat Biotechnol. 2022:40(12):1794–1806. https://doi.org/10.1038/s41587-022-01483-z

Heumos L, Schaar AC, Lance C, Litinetskaya A, Drost F, Zappia L, Lücken MD, Strobl DC, Henao J, Curion F, et al. Best practices for single-cell analysis across modalities. Nat Rev Genet. 2023:24(8): 550–572. https://doi.org/10.1038/s41576-023-00586-w

Janesick A, Shelansky R, Gottscho AD, Wagner F, Rouault M, Beliakoff G, de Oliveira MF, Kohlway A, Abousoud J, Morrison CA, et al. High resolution mapping of the breast cancer tumor microenvironment using integrated single cell, spatial and in situ analysis of FFPE tissue. bioRxiv 510405. https://doi.org/10.1101/ 2022.10.06.510405, 07 October 2022, preprint: not peer reviewed.

Jean-Baptiste K, McFaline-Figueroa JL, Alexandre CM, Dorrity MW, Saunders L, Bubb KL, Trapnell C, Fields S, Queitsch C, Cuperus JT. Dynamics of gene expression in single root cells of *Arabidopsis thaliana*. Plant Cell. 2019:31(5):993–1011. https://doi.org/10.1105/ tpc.18.00785

Jin J, Lu P, Xu Y, Tao J, Li Z, Wang S, Yu S, Wang C, Xie X, Gao J, et al. PCMDB: a curated and comprehensive resource of plant cell markers. Nucleic Acids Res. 2022:50(D1):D1448–D1455. https://doi.org/10. 1093/nar/gkab949

Kao P, Schon MA, Mosiolek M, Enugutti B, Nodine MD. Gene expres- sion variation in *Arabidopsis* embryos at single-nucleus resolution. Development. 2021:148(13):dev199589. https://doi.org/10.1242/dev. 199589

Ke R, Mignardi M, Pacureanu A, Svedlund J, Botling J, Wählby C, Nilsson M. In situ sequencing for RNA analysis in preserved tissue
end-cycle in root development and stress adaptation. Plant Cell.
2018;30(10):2330-2351. https://doi.org/10.1105/tpc.17.00983

Birnbaum K, Jung JW, Wang JY, Lambert GM, Hirst JA, Galbraith DW, Benfey PN. Cell type-specific expression profiling in plants via cell sorting of protoplasts from fluorescent reporter lines. Nat Methods. 2005;2(8):615–619. https://doi.org/10.1038/ nmeth0805-615

Booeshaghi AS, Chen X, Pachter L. A machine-readable specification for genomics assays. bioRxiv 533215. https://doi.org/10.1101/2023.03. 17.533215, 21 March 2023, preprint: not peer reviewed.

Brady SM, Orlando DA, Lee J-Y, Wang JY, Koch J, Dinneny JR, Mace D, Ohler U, Benfey PN. A high-resolution root spatiotemporal map reveals dominant expression patterns. Science. 2007;318(5851): 801–806. https://doi.org/10.1126/science.1146265

Cao J, Packer JS, Ramani V, Cusanovich DA, Huynh C, Daza R, Qiu X, Lee C, Furlan SN, Steemers FJ, et al. Comprehensive single-cell tran- scriptional profiling of a multicellular organism. Science. 2017;357(6352):661–667. https://doi.org/10.1126/science.aam8940

Cartwright DA, Brady SM, Orlando DA, Sturmfels B, Benfey PN. Reconstructing spatiotemporal gene expression data from partial ob- servations. Bioinformatics. 2009;25(19):2581–2587. https://doi.org/ 10.1093/bioinformatics/btp437

Cervantes-Pérez SA, Thibivilliers S, Laffont C, Farmer AD, Frugier F, Libault M. Cell-specific pathways recruited for symbiotic nodulation in the *Medicago truncatula* legume. Mol Plant. 2022;15(12):1868–1888. https://doi.org/10.1016/j.molp.2022.10.021

Chari T, Pachter L. The specious art of single-cell genomics. PLoS Comput Biol. 2023;19(8):e1011288. https://doi.org/10.1371/journal. pcbi.1011288

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

Clarke ZA, Andrews TS, Atif J, Pouyabahar D, Innes BT, MacParland SA, Bader GD. Tutorial: guidelines for annotating single-cell tran- scriptomic maps using automated and manual methods. Nat Protoc. 2021;16(6):2749–2764. https://doi.org/10.1038/s41596-021- 00534-0

Conde D, Triozzi PM, Pereira WJ, Schmidt HW, Balmant KM, Knaack SA, Redondo-López A, Roy S, Dervinis C, Kirst M. Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in *Populus*. Development. 2022;149(21):dev200632. https://doi.org/10.1242/dev.200632

Denyer T, Ma X, Klesen S, Scacchi E, Nieselt K, Timmermans MCP. Spatiotemporal developmental trajectories in the *Arabidopsis* root revealed using high-throughput single-cell RNA sequencing. Dev Cell. 2019;48(6):840–852.e5. https://doi.org/10.1016/j.devcel. 2019.02.022

Efroni I, Mello A, Nawy T, Ip P-L, Rahni R, DelRose N, Powers A, Satija R, Birnbaum KD. Root regeneration triggers an embryo-like sequence guided by hormonal interactions. Cell. 2016;165(7): 1721–1733. https://doi.org/10.1016/j.cell.2016.04.046

Fahlgren N, Kapoor M, Yordanova G, Papatheodorou I, Waese J, Cole B, Harrison P, Ware D, Tickle T, Paten B, et al. Toward a data infrastructure for the Plant Cell Atlas. Plant Physiol. 2023;191(1):35–46. https://doi.org/10.1093/plphys/kiac468

Farmer A, Thibivilliers S, Ryu KH, Schiefelbein J, Libault M. 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. 2021;14(3):372–383. https://doi.org/ 10.1016/j.molp.2021.01.001

Ferrari C, Pérez NM, Vandepoele K. MINI-EX: integrative inference of single-cell gene regulatory networks in plants. Mol Plant. 2022;15(11): 1807–1824. https://doi.org/10.1016/j.molp.2022.10.016

Fischer S, Gillis J. How many markers are needed to robustly determine a cell's type? iScience. 2021;24(11):103292. https://doi.org/10.1016/j.isci.2021.103292

Fleming SJ, Chaffin MD, Arduini A, Akkad A-D, Banks E, Marioni JC, Philippakis AA, Ellinor PT, Babadi M. Unsupervised removal of systematic background noise from droplet-based single-cell experiments using CellBender. Nat Methods. 2023;20(9):1323–1335. https://doi.org/10.1038/s41592-023-01943-7

Fox DT, Soltis DE, Soltis PS, Ashman T-L, Van de Peer Y. Polyploidy: a biological force from cells to ecosystems. Trends Cell Biol. 2020;30(9): 688–694. https://doi.org/10.1016/j.tcb.2020.06.006

Gala HP, Lancot A, Jean-Baptiste K, Guiziou S, Chu JC, Zemke JE, George W, Queitsch C, Cuperus JT, Nemhauser JL. A single-cell view of the transcriptome during lateral root initiation in *Arabidopsis thaliana*. Plant Cell. 2021;33(7):2197–2220. https://doi.org/10.1093/plcell/koab101

Giacomello S, Salmén F, Terebieniec BK, Vickovic S, Navarro JF, Alexeyenko A, Reimegård J, McKee LS, Mannapperuma C, Bulone V, et al. Spatially resolved transcriptome profiling in model plant species. Nat Plants. 2017;3(6):1–11. https://doi.org/10.1038/ nplants.2017.61

Graeff M, Rana S, Wendrich JR, Dorier J, Eekhout T, Aliaga Fandino AC, Guex N, Bassel GW, De Rybel B, Hardtke CS. A single-cell morpho-transcriptomic map of brassinosteroid action in the *Arabidopsis* root. Mol Plant. 2021;14(12):1985–1999. https://doi.org/10.1016/j.molp.2021.07.021

Groiss S, Pabst D, Faber C, Meier A, Bogdoll A, Unger C, Nilges B, Strauss S, Föderl-Höbenreich E, Hardt M, et al. Highly resolved spa- tial transcriptomics for detection of rare events in cell. bioRxiv 463936. https://doi.org/10.1101/2021.10.11.463936, 12 October 2021, preprint: not peer reviewed.

Guillotin B, Rahni R, Passalacqua M, Mohammed MA, Xu X, Raju SK, Ramírez CO, Jackson D, Groen SC, Gillis J, et al. A pan-grass transcriptome reveals patterns of cellular divergence in crops. Nature. 2023;617(7962):785–791. https://doi.org/10.1038/s41586- 023-06053-0

Hao Y, Hao S, Andersen-Nissen E, Mauck WM III, Zheng S, Butler A, Lee MJ, Wilk AJ, Darby C, Zager M, et al. Integrated analysis of multimodal single-cell data. Cell. 2021;184(13):3573–3587.e29. https://doi.org/10.1016/j.cell.2021.04.048

He S, Bhatt R, Brown C, Brown EA, Buhr DL, Chantranuvatana K, Danaher P, Dunaway D, Garrison RG, Geiss G, et al. High-plex im- aging of RNA and proteins at subcellular resolution in fixed tissue by spatial molecular imaging. Nat Biotechnol. 2022;40(12):1794–1806. https://doi.org/10.1038/s41587-022-01483-z

Heumos L, Schaar AC, Lance C, Litinetskaya A, Drost F, Zappia L, Lücken MD, Strobl DC, Henao J, Curion F, et al. Best practices for single-cell analysis across modalities. Nat Rev Genet. 2023;24(8): 550–572. https://doi.org/10.1038/s41576-023-00586-w

Janesick A, Shelansky R, Gottscho AD, Wagner F, Rouault M, Beliakoff G, de Oliveira MF, Kohlway A, Abousoud J, Morrison CA, et al. High resolution mapping of the breast cancer tumor microenvironment using integrated single cell, spatial and in situ analysis of FFPE tissue. bioRxiv 510405. https://doi.org/10.1101/ 2022.10.06.510405, 07 October 2022, preprint: not peer reviewed.

Jean-Baptiste K, McFaline-Figueroa JL, Alexandre CM, Dorrity MW, Saunders L, Bubb KL, Trapnell C, Fields S, Queitsch C, Cuperus JT. Dynamics of gene expression in single root cells of *Arabidopsis thaliana*. Plant Cell. 2019;31(5):993–1011. https://doi.org/10.1105/ tpc.18.00785

Jin J, Lu P, Xu Y, Tao J, Li Z, Wang S, Yu S, Wang C, Xie X, Gao J, et al. PCMDB: a curated and comprehensive resource of plant cell markers. Nucleic Acids Res. 2022;50(D1):D1448–D1455. https://doi.org/10. 1093/nar/gkab949

Kao P, Schon MA, Mosiolek M, Enugutti B, Nodine MD. Gene expres- sion variation in *Arabidopsis* embryos at single-nucleus resolution. Development. 2021;148(13):dev199589. https://doi.org/10.1242/dev. 199589

Ke R, Mignardi M, Pacureanu A, Svedlund J, Botling J, Wählby C, Nilsson M. In situ sequencing for RNA analysis in preserved tissue
and cells. Nat Methods. 2013:10(9):857–860. https://doi.org/10.1038/nmeth.2563

Kim J-Y, Symeonidi E, Pang TY, Denyer T, Weidauer D, Bezrutczyk M, Miras M, Zöllner N, Hartwig T, Wudick MM, et al. Distinct identities of leaf phloem cells revealed by single cell transcriptomics. Plant Cell. 2021:33(3):511–530. https://doi.org/10.1093/plcell/koaa060

Kim E-J, Zhang C, Guo B, Eekhout T, Houbaert A, Wendrich JR, Vandamme N, Tiwari M, Simon-Vezo C, Vanhoutte I, et al. Cell type-specific attenuation of brassinosteroid signaling precedes stomatal asymmetric cell division. Proc Natl Acad Sci U S A. 2023:120(36):e2303758120. https://doi.org/10.1073/pnas.2303758120

Kiselev VY, Yiu A, Hemberg M. Smap: projection of single-cell RNA-seq data across data sets. Nat Methods. 2018:15(5):359–362. https://doi.org/10.1038/nmeth.4644

Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, Baglaenko Y, Brenner M, Loh PR, Raychaudhuri S. Fast, sensitive and accurate integration of single-cell data with harmony. Nat Methods. 2019:16(12):1289–1296. https://doi.org/10.1038/s41592-019-0619-0

Laureyns R, Joossens J, Herwegh D, Pevernagie J, Pavie B, Demuynck K, Debray K, Coussens G, Pauwels L, Van Hautegem T, et al. An in situ sequencing approach maps PLASTOCHRON1 at the boundary between indeterminate and determinate cells. Plant Physiol. 2021:188(2):782–794. https://doi.org/10.1093/plphys/kiab533

Lebrigand K, Magnone V, Barbry P, Waldmann R. High throughput error corrected Nanopore single cell transcriptome sequencing. Nat Commun. 2020:11(1):4025. https://doi.org/10.1038/s41467-020-17800-6

Lee TA, Bailey-Serres J. Integrative analysis from the epigenome to translatome uncovers patterns of dominant nuclear regulation during transient stress. Plant Cell. 2019:31(11):2573–2595. https://doi.org/10.1105/tpc.19.00463

Lee JH, Daugharthy ER, Scheiman J, Kalhor R, Ferrante TC, Terry R, Turczyk BM, Yang JL, Lee HS, Aach J, et al. Fluorescent in situ sequencing (FISSEQ) of RNA for gene expression profiling in intact cells and tissues. Nat Protoc. 2015:10(3):442–458. https://doi.org/10.1038/nprot.2014.191

Lee TA, Nobori T, Illouz-Eliaz N, Xu J, Jow B, Nery JR, Ecker JR. A single-nucleus Atlas of seed-to-seed development in *Arabidopsis*. bioRxiv 533992. https://doi.org/10.1101/2023.03.23.533992, 21 March 2023, preprint: not peer reviewed.

Leote AC, Wu X, Beyer A. Regulatory network-based imputation of dropouts in single-cell RNA sequencing data. PLoS Comput Biol. 2022:18(2):e1009849. https://doi.org/10.1371/journal.pcbi.1009849

Li X, Zhang X, Gao S, Cui F, Chen W, Fan L, Qi Y. Single-cell RNA sequencing reveals the landscape of maize root tips and assists in identification of cell type-specific nitrate-response genes. Crop J. 2022:10(6):1589–1600. https://doi.org/10.1016/j.cj.2022.02.004

Li C, Zhang S, Yan X, Cheng P, Yu H. Single-nucleus sequencing deciphers developmental trajectories in rice pistils. Dev Cell. 2023:58(8):694–708.e4. https://doi.org/10.1016/j.devcel.2023.03.004

Liu Z, Kong X, Long Y, Zhang H, Jia J, Qiu L, Zhai J, Yan Z. Integrated single-nucleus and spatial transcriptomics captures transitional states in soybean nodule symbiosis establishment. Nat Plants. 2023:9(4):515–524. https://doi.org/10.1038/s41477-023-01387-z, 02 July 2022b.

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

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

Liu Y, Yang M, Deng Y, Su G, Enninful A, Guo CC, Tebaldi T, Zhang D, Kim D, Bai Z, et al. High-spatial-resolution multi-omics

sequencing via deterministic barcoding in tissue. Cell. 2020:183(6):1665–1681.e18. https://doi.org/10.1016/j.cell.2020.10.026

Lopez-Anido CB, Vatén A, Smoot NK, Sharma N, Guo V, Gong Y, Anleu Gil MX, Weimer AK, Bergmann DC. Single-cell resolution of lineage trajectories in the *Arabidopsis* stomatal lineage and developing leaf. Dev Cell. 2021:56(7):1043–1055.e4. https://doi.org/10.1016/j.devcel.2021.03.014

Luecken MD, Büttner M, Chaichoompu K, Danese A, Interlandi M, Mueller MF, Strobl DC, Zappia L, Dugas M, Colomé-Tatché M, et al. Benchmarking atlas-level data integration in single-cell genomics. Nat Methods. 2022:19(1):41–50. https://doi.org/10.1038/s41592-021-01336-8

Luecken MD, Theis FJ. Current best practices in single-cell RNA-seq analysis: a tutorial. Mol Syst Biol. 2019:15(6):e8746. https://doi.org/10.15252/msb.20188746

Ma F, Salomé PA, Merchant SS, Pellegrini M. Single-cell RNA sequencing of batch Chlamydomonas cultures reveals heterogeneity in their diurnal cycle phase. Plant Cell. 2021:33(4):1042–1057. https://doi.org/10.1093/plcell/koab025

Marand AP, Chen Z, Gallavotti A, Schmitz RJ. A cis-regulatory atlas in maize at single-cell resolution. Cell. 2021:184(11):3041–3055.e21. https://doi.org/10.1016/j.cell.2021.04.014

Marchant DB, Nelms B, Walbot V. FX-Cell: quantitative cell release from fixed plant tissues for single-cell genomics. bioRxiv 463960. https://doi.org/10.1101/2021.10.11.463960, 14 September 2022, preprint: not peer reviewed.

McGinnis CS, Patterson DM, Winkler J, Conrad DN, Hein MY, Srivastava V, Hu JL, Murrow LM, Weissman JS, Werb Z, et al. MULTI-seq: sample multiplexing for single-cell RNA sequencing using lipid-tagged indices. Nat Methods. 2019:16(7):619–626. https://doi.org/10.1038/s41592-019-0433-8

Mereu E, Lafzi A, Moutinho C, Ziegenhain C, McCarthy DJ, Álvarez-Varela A, Batlle E, Sagar, Grün D, Lau JK, et al. Benchmarking single-cell RNA-sequencing protocols for cell atlas projects. Nat Biotechnol. 2020:38(6):747–755. https://doi.org/10.1038/s41587-020-0469-4

Michielsen L, Reinders MJT, Mahfouz A. Hierarchical progressive learning of cell identities in single-cell data. Nat Commun. 2021:12(1):2799. https://doi.org/10.1038/s41467-021-23196-8

Miyashima S, Roszak P, Sevilem I, Toyokura K, Blob B, Heo JO, Mellor N, Help-Rinta-Rahko H, Otero S, Smet W, et al. Mobile PEAR transcription factors integrate positional cues to prime cambial growth. Nature. 2019:565(7740):490–494. https://doi.org/10.1038/s41586-018-0839-y

Moffitt JR, Hao J, Wang G, Chen KH, Babcock HP, Zhuang X. High-throughput single-cell gene-expression profiling with multiplexed error-robust fluorescence in situ hybridization. Proc Natl Acad Sci U S A. 2016:113(39):11046–11051. https://doi.org/10.1073/pnas.1612826113

Movahedi S, Van Bel M, Heyndrickx KS, Vandepoele K. Comparative co-expression analysis in plant biology. Plant Cell Environ. 2012:35(10):1787–1798. https://doi.org/10.1111/j.1365-3040.2012.02517.x

Narsai R, Howell KA, Millar AH, O'Toole N, Small I, Whelan J. Genome-wide analysis of mRNA decay rates and their determinants in *Arabidopsis thaliana*. Plant Cell. 2007:19(11):3418–3436. https://doi.org/10.1105/tpc.107.055046

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

Neumann M, Xu X, Smacznik C, Schumacher J, Yan W, Blüthgen N. Greb T, Jönsson H, Traas J, Kaufmann K, et al. A 3D gene expression atlas of the floral meristem based on spatial reconstruction of single nucleus RNA sequencing data. Nat Commun. 2022:13(1):2838. https://doi.org/10.1038/s41467-022-30177-y

Nobori T, Monell A, Lee TA, Zhou J, Nery JR, Ecker JR. Time-resolved single-cell and spatial gene regulatory atlas of plants under pathogen
and cells. Nat Methods. 2013:10(9):857–860. https://doi.org/10.1038/nmeth.2563

Kim J-Y, Symeonidi E, Pang TY, Denyer T, Weidauer D, Bezrutczyk M, Miras M, Zöllner N, Hartwig T, Wudick MM, et al. Distinct identities of leaf phloem cells revealed by single cell transcriptomics. Plant Cell. 2021:33(3):511–530. https://doi.org/10.1093/plcell/koaa060

Kim E-J, Zhang C, Guo B, Eekhout T, Houbaert A, Wendrich JR, Vandamme N, Tiwari M, Simon-Vezo C, Vanhoutte I, et al. Cell type-specific attenuation of brassinosteroid signaling precedes stomatal asymmetric cell division. Proc Natl Acad Sci U S A. 2023:120(36):e2303758120. https://doi.org/10.1073/pnas.2303758120

Kiselev VY, Yiu A, Hemberg M. Smap: projection of single-cell RNA-seq data across data sets. Nat Methods. 2018:15(5):359–362. https://doi.org/10.1038/nmeth.4644

Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, Baglaenko Y, Brenner M, Loh PR, Raychaudhuri S. Fast, sensitive and accurate integration of single-cell data with harmony. Nat Methods. 2019:16(12):1289–1296. https://doi.org/10.1038/s41592-019-0619-0

Laureyns R, Joossens J, Herwegh D, Pevernagie J, Pavie B, Demuynck K, Debray K, Coussens G, Pauwels L, Van Hautegem T, et al. An in situ sequencing approach maps PLASTOCHRON1 at the boundary between indeterminate and determinate cells. Plant Physiol. 2021:188(2):782–794. https://doi.org/10.1093/plphys/kiab533

Lebrigand K, Magnone V, Barbry P, Waldmann R. High throughput error corrected Nanopore single cell transcriptome sequencing. Nat Commun. 2020:11(1):4025. https://doi.org/10.1038/s41467-020-17800-6

Lee TA, Bailey-Serres J. Integrative analysis from the epigenome to translatome uncovers patterns of dominant nuclear regulation during transient stress. Plant Cell. 2019:31(11):2573–2595. https://doi.org/10.1105/tpc.19.00463

Lee JH, Daugharthy ER, Scheiman J, Kalhor R, Ferrante TC, Terry R, Turczyk BM, Yang JL, Lee HS, Aach J, et al. Fluorescent in situ sequencing (FISSEQ) of RNA for gene expression profiling in intact cells and tissues. Nat Protoc. 2015:10(3):442–458. https://doi.org/10.1038/nprot.2014.191

Lee TA, Nobori T, Illouz-Eliaz N, Xu J, Jow B, Nery JR, Ecker JR. A single-nucleus Atlas of seed-to-seed development in *Arabidopsis*. bioRxiv 533992. https://doi.org/10.1101/2023.03.23.533992, 21 March 2023, preprint: not peer reviewed.

Leote AC, Wu X, Beyer A. Regulatory network-based imputation of dropouts in single-cell RNA sequencing data. PLoS Comput Biol. 2022:18(2):e1009849. https://doi.org/10.1371/journal.pcbi.1009849

Li X, Zhang X, Gao S, Cui F, Chen W, Fan L, Qi Y. Single-cell RNA sequencing reveals the landscape of maize root tips and assists in identification of cell type-specific nitrate-response genes. Crop J. 2022:10(6):1589–1600. https://doi.org/10.1016/j.cj.2022.02.004

Li C, Zhang S, Yan X, Cheng P, Yu H. Single-nucleus sequencing deciphers developmental trajectories in rice pistils. Dev Cell. 2023:58(8):694–708.e4. https://doi.org/10.1016/j.devcel.2023.03.004

Liu Z, Kong X, Long Y, Zhang H, Jia J, Qiu L, Zhai J, Yan Z. Integrated single-nucleus and spatial transcriptomics captures transitional states in soybean nodule symbiosis establishment. Nat Plants. 2023:9(4):515–524. https://doi.org/10.1038/s41477-023-01387-z, 02 July 2022b.

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

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

Liu Y, Yang M, Deng Y, Su G, Enninful A, Guo CC, Tebaldi T, Zhang D, Kim D, Bai Z, et al. High-spatial-resolution multi-omics

sequencing via deterministic barcoding in tissue. Cell. 2020:183(6):1665–1681.e18. https://doi.org/10.1016/j.cell.2020.10.026

Lopez-Anido CB, Vatén A, Smoot NK, Sharma N, Guo V, Gong Y, Anleu Gil MX, Weimer AK, Bergmann DC. Single-cell resolution of lineage trajectories in the *Arabidopsis* stomatal lineage and developing leaf. Dev Cell. 2021:56(7):1043–1055.e4. https://doi.org/10.1016/j.devcel.2021.03.014

Luecken MD, Büttner M, Chaichoompu K, Danese A, Interlandi M, Mueller MF, Strobl DC, Zappia L, Dugas M, Colomé-Tatché M, et al. Benchmarking atlas-level data integration in single-cell genomics. Nat Methods. 2022:19(1):41–50. https://doi.org/10.1038/s41592-021-01336-8

Luecken MD, Theis FJ. Current best practices in single-cell RNA-seq analysis: a tutorial. Mol Syst Biol. 2019:15(6):e8746. https://doi.org/10.15252/msb.20188746

Ma F, Salomé PA, Merchant SS, Pellegrini M. Single-cell RNA sequencing of batch Chlamydomonas cultures reveals heterogeneity in their diurnal cycle phase. Plant Cell. 2021:33(4):1042–1057. https://doi.org/10.1093/plcell/koab025

Marand AP, Chen Z, Gallavotti A, Schmitz RJ. A cis-regulatory atlas in maize at single-cell resolution. Cell. 2021:184(11):3041–3055.e21. https://doi.org/10.1016/j.cell.2021.04.014

Marchant DB, Nelms B, Walbot V. FX-Cell: quantitative cell release from fixed plant tissues for single-cell genomics. bioRxiv 463960. https://doi.org/10.1101/2021.10.11.463960, 14 September 2022, preprint: not peer reviewed.

McGinnis CS, Patterson DM, Winkler J, Conrad DN, Hein MY, Srivastava V, Hu JL, Murrow LM, Weissman JS, Werb Z, et al. MULTI-seq: sample multiplexing for single-cell RNA sequencing using lipid-tagged indices. Nat Methods. 2019:16(7):619–626. https://doi.org/10.1038/s41592-019-0433-8

Mereu E, Lafzi A, Moutinho C, Ziegenhain C, McCarthy DJ, Álvarez-Varela A, Batlle E, Sagar, Grün D, Lau JK, et al. Benchmarking single-cell RNA-sequencing protocols for cell atlas projects. Nat Biotechnol. 2020:38(6):747–755. https://doi.org/10.1038/s41587-020-0469-4

Michielsen L, Reinders MJT, Mahfouz A. Hierarchical progressive learning of cell identities in single-cell data. Nat Commun. 2021:12(1):2799. https://doi.org/10.1038/s41467-021-23196-8

Miyashima S, Roszak P, Sevilem I, Toyokura K, Blob B, Heo JO, Mellor N, Help-Rinta-Rahko H, Otero S, Smet W, et al. Mobile PEAR transcription factors integrate positional cues to prime cambial growth. Nature. 2019:565(7740):490–494. https://doi.org/10.1038/s41586-018-0839-y

Moffitt JR, Hao J, Wang G, Chen KH, Babcock HP, Zhuang X. High-throughput single-cell gene-expression profiling with multiplexed error-robust fluorescence in situ hybridization. Proc Natl Acad Sci U S A. 2016:113(39):11046–11051. https://doi.org/10.1073/pnas.1612826113

Movahedi S, Van Bel M, Heyndrickx KS, Vandepoele K. Comparative co-expression analysis in plant biology. Plant Cell Environ. 2012:35(10):1787–1798. https://doi.org/10.1111/j.1365-3040.2012.02517.x

Narsai R, Howell KA, Millar AH, O'Toole N, Small I, Whelan J. Genome-wide analysis of mRNA decay rates and their determinants in *Arabidopsis thaliana*. Plant Cell. 2007:19(11):3418–3436. https://doi.org/10.1105/tpc.107.055046

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

Neumann M, Xu X, Smacznik C, Schumacher J, Yan W, Blüthgen N. Greb T, Jönsson H, Traas J, Kaufmann K, et al. A 3D gene expression atlas of the floral meristem based on spatial reconstruction of single nucleus RNA sequencing data. Nat Commun. 2022:13(1):2838. https://doi.org/10.1038/s41467-022-30177-y

Nobori T, Monell A, Lee TA, Zhou J, Nery JR, Ecker JR. Time-resolved single-cell and spatial gene regulatory atlas of plants under pathogen
attack. bioRxiv 536170. https://doi.org/10.1101/2023.04.10.536170, 10 April 2023, preprint: not peer reviewed.

**Omary M, Gil-Yarom N, Yahav C, Steiner E, Hendelman A, Efroni I.** A conserved superlocus regulates above- and belowground root initiation. Science. 2022;375(6584):993. https://doi.org/10.1126/science.abf4368

**Ortiz-Ramírez C, Arevalo ED, Xu X, Jackson DP, Birnbaum KD.** An efficient cell sorting protocol for maize protoplasts. Curr Protoc Plant Biol. 2018;3(3):e20072. https://doi.org/10.1002/cppb.20072

**Ortiz-Ramírez C, Guillotin B, Xu X, Rahni R, Zhang S, Yan Z, Coqueiro Dias Araujo P, Demesa-Arevalo E, Lee L, Van Eck J, et al.** Ground tissue circuitry regulates organ complexity in maize and Setaria. Science. 2021;374(6572):1247–1252. https://doi.org/10.1126/science.abj2327

**Otero S, Gildea I, Roszak P, Lu Y, Di Vittori V, Bourdon M, Kalmbach L, Blob B, Heo JO, Peruzzo F, et al.** A root phloem pole cell atlas reveals common transcriptional states in protophloem-adjacent cells. Nat Plants. 2022;8(8):954–970. https://doi.org/10.1038/s41477-022-01178-y

**Peirats-Llobet M, Yi C, Liew LC, Berkowitz O, Narsai R, Lewsey MG, Whelan J.** Spatially resolved transcriptomic analysis of the germinating barley grain. Nucleic Acids Res. 2023;51(15):7798–7819. https://doi.org/10.1093/nar/gkad521

**Phan HV, van Gent M, Drayman N, Basu A, Gack MU, Tay S.** High-throughput RNA sequencing of paraformaldehyde-fixed single cells. Nat Commun. 2021;12(1):5636. https://doi.org/10.1038/s41467-021-25871-2

**Picard CL, Povilus RA, Williams BP, Gehring M.** Transcriptional and imprinting complexity in *Arabidopsis* seeds at single-nucleus resolution. Nat Plants. 2021;7(6):730–738. https://doi.org/10.1038/s41477-021-00922-0

**Reynoso MA, Kajala K, Bajic M, West DA, Pauluzzi G, Yao AI, Hatch K, Zumstein K, Woodhouse M, Rodriguez-Medina J, et al.** Evolutionary flexibility in flooding response circuitry in angiosperms. Science. 2019;365(6459):1291–1295. https://doi.org/10.1126/science.aax8862

**Rosenberg AB, Roco CM, Muscat RA, Kuchina A, Sample P, Yao Z, Graybuck LT, Peeler DJ, Mukherjee S, Chen W, et al.** Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding. Science. 2018;360(6385):176–182. https://doi.org/10.1126/science.aam8999

**Roszak P, Heo JO, Blob B, Toyokura K, Sugiyama Y, de Luis Balaguer MA, Lau WWY, Hamey F, Cirrone J, Madej E, et al.** Cell-by-cell dissection of phloem development links a maturation gradient to cell specialization. Science. 2021;374(6575):eaba5531. https://doi.org/10.1126/science.aba5531

**Ryu KH, Huang L, Kang HM, Schiefelbein J.** Single-cell RNA sequencing resolves molecular relationships among individual plant cells. Plant Physiol. 2019;179(4):1444–1456. https://doi.org/10.1104/pp.18.01482

**Saelens W, Cannoodt R, Todorov H, Saeys Y.** A comparison of single-cell trajectory inference methods. Nat Biotechnol. 2019;37(5):547–554. https://doi.org/10.1038/s41587-019-0071-9

**Satterlee JW, Strable J, Scanlon MJ.** Plant stem-cell organization and differentiation at single-cell resolution. Proc Natl Acad Sci U S A. 2020;117(52):33689–33699. https://doi.org/10.1073/pnas.201878117

**Serrano-Ron L, Perez-Garcia P, Sanchez-Corrieronero A, Gude I, Cabrera J, Ip P-L, Birnbaum KD, Moreno-Risueno MA.** Reconstruction of lateral root formation through single-cell RNA sequencing reveals order of tissue initiation. Mol Plant. 2021;14(8):1362–1378. https://doi.org/10.1016/j.molp.2021.05.028

**Shahan R, Hsu C-W, Nolan TM, Cole BJ, Taylor IW, Greenstreet L, Zhang S, Afanassiev A, Vlot AHC, Schiebinger G, et al.** A single-cell *Arabidopsis* root atlas reveals developmental trajectories in wild-type and cell identity mutants. Dev Cell. 2022;57(4):543–560.e9. https://doi.org/10.1016/j.devcel.2022.01.008

**Shen L, Ma J, Li P, Wu Y, Yu H.** Recent advances in the plant epitranscriptome. Genome Biol. 2023;24(1):43. https://doi.org/10.1186/s13059-023-02872-6

**Shulse CN, Cole BJ, Ciobanu D, Lin J, Yoshinaga Y, Gouran M, Turco GM, Zhu Y, O'Malley RC, Brady SM, et al.** High-throughput single-cell transcriptome profiling of plant cell types. Cell Rep. 2019;27(7):2241–2247.e4. https://doi.org/10.1016/j.celrep.2019.04.054

**Song Q, Ando A, Jiang N, Ikeda Y, Chen ZJ.** Single-cell RNA-seq analysis revealsddy-dependent and cell-specific transcriptome changes in *Arabidopsis* female gametophytes. Genome Biol. 2020;21(1):178. https://doi.org/10.1186/s13059-020-02094-0

**Squair JW, Gautier M, Kathe C, Anderson MA, James ND, Hutson TH, Hudelle R, Qaiser T, Matson KJE, Barraud Q, et al.** Confronting false discoveries in single-cell differential expression. Nat Commun. 2021;12(1):S692. https://doi.org/10.1038/s41467-021-25960-2

**Stickels RR, Murray E, Kumar P, Li J, Marshall JL, Di Bella DJ, Arlotta P, Macosko EZ, Chen F.** Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2. Nat Biotechnol. 2021;39(3):313–319. https://doi.org/10.1038/s41587-020-0739-1

**Stoeckius M, Zheng S, Houck-Loomis B, Hao S, Yeung BZ, Mauck WM III, Smibert P, Satija R.** Cell hashing with barcoded antibodies enables multiplexing and doublet detection for single cell genomics. Genome Biol. 2018;19(1):224. https://doi.org/10.1186/s13059-018-1603-1

**Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM III, Hao Y, Stoeckius M, Smibert P, Satija R.** Comprehensive integration of single-cell data. Cell. 2019;177(7):1888–1902.e21. https://doi.org/10.1016/j.cell.2019.05.031

**Sun G, Xia M, Li J, Ma W, Li Q, Xie J, Bai S, Fang S, Sun T, Feng X, et al.** The maize single-nucleus transcriptome comprehensively describes signaling networks governing movement and development of grass stomata. Plant Cell. 2022;34(5):1890–1911. https://doi.org/10.1093/plcell/koac047

**Sunaga-Franze DY, Muino JM, Braeuning C, Xu X, Zong M, Smaczniak C, Yan W, Fischer C, Vidal R, Kliem M, et al.** Single-nucleus RNA sequencing of plant tissues using a nanowell-based system. Plant J. 2021;108(3):859–869. https://doi.org/10.1111/tpj.15458

**Tan Y, Cahan P.** SingleCellNet: a computational tool to classify single cell RNA-seq data across platforms and across species. Cell Syst. 2019;9:207–213.e3. https://doi.org/10.1016/j.cels.2019.06.004

**Tao S, Liu P, Shi Y, Feng Y, Gao J, Chen L, Zhang A, Cheng X, Wei H, Zhang T, et al.** Single-cell transcriptome and network analyses unveil key transcription factors regulating mesophyll cell development in maize. Genes (Basel). 2022;13(2):374. https://doi.org/10.3390/genes13020374

**Tian C, Du Q, Xu M, Du F, Jiao Y.** Single-nucleus RNA-seq resolves spatiotemporal developmental trajectories in the tomato shoot apex. bioRxiv 305029. https://doi.org/10.1101/2020.09.20.305029, 20 September 2020, preprint: not peer reviewed.

**Tran HTN, Ang KS, Chevrier M, Zhang X, Lee NYS, Goh M, Chen J.** A benchmark of batch-effect correction methods for single-cell RNA sequencing data. Genome Biol. 2020;21(1):12. https://doi.org/10.1186/s13059-019-1850-9

**Tritschler S, Büttner M, Fischer DS, Lange M, Bergen V, Lickert H, Theis FJ.** Concepts and limitations for learning developmental trajectories from single cell genomics. Development. 2019;146(12):dev170506. https://doi.org/10.1242/dev.170506

**Turco GM, Rodriguez-Medina J, Siebert S, Han D, Valderrama-Gómez MÁ, Vahldick H, Shulse CN, Cole BJ, Juliano CE, Dickel DE, et al.** Molecular mechanisms driving switch behavior in xylem cell differentiation. Cell Rep. 2019;28(2):342–351.e4. https://doi.org/10.1016/j.celrep.2019.06.041

**Van den Berge K, Roux de Bézieux H, Street K, Saelens W, Cannoodt R, Saeys Y, Dudoit S, Clement L.** Trajectory-based differential expression analysis for single-cell sequencing data. Nat Commun. 2020;11(1):1201. https://doi.org/10.1038/s41467-020-14766-3
attack. bioRxiv 536170. https://doi.org/10.1101/2023.04.10.536170, 10 April 2023, preprint: not peer reviewed.

**Omary M, Gil-Yarom N, Yahav C, Steiner E, Hendelman A, Efroni I.** A conserved superlocus regulates above- and belowground root initiation. Science. 2022;375(6584):993. https://doi.org/10.1126/science.abf4368

**Ortiz-Ramírez C, Arevalo ED, Xu X, Jackson DP, Birnbaum KD.** An efficient cell sorting protocol for maize protoplasts. Curr Protoc Plant Biol. 2018;3(3):e20072. https://doi.org/10.1002/cppb.20072

**Ortiz-Ramírez C, Guillotin B, Xu X, Rahni R, Zhang S, Yan Z, Coqueiro Dias Araujo P, Demesa-Arevalo E, Lee L, Van Eck J, et al.** Ground tissue circuitry regulates organ complexity in maize and Setaria. Science. 2021;374(6572):1247–1252. https://doi.org/10.1126/science.abj2327

**Otero S, Gildea I, Roszak P, Lu Y, Di Vittori V, Bourdon M, Kalmbach L, Blob B, Heo JO, Peruzzo F, et al.** A root phloem pole cell atlas reveals common transcriptional states in protophloem-adjacent cells. Nat Plants. 2022;8(8):954–970. https://doi.org/10.1038/s41477-022-01178-y

**Peirats-Llobet M, Yi C, Liew LC, Berkowitz O, Narsai R, Lewsey MG, Whelan J.** Spatially resolved transcriptomic analysis of the germinating barley grain. Nucleic Acids Res. 2023;51(15):7798–7819. https://doi.org/10.1093/nar/gkad521

**Phan HV, van Gent M, Drayman N, Basu A, Gack MU, Tay S.** High-throughput RNA sequencing of paraformaldehyde-fixed single cells. Nat Commun. 2021;12(1):5636. https://doi.org/10.1038/s41467-021-25871-2

**Picard CL, Povilus RA, Williams BP, Gehring M.** Transcriptional and imprinting complexity in *Arabidopsis* seeds at single-nucleus resolution. Nat Plants. 2021;7(6):730–738. https://doi.org/10.1038/s41477-021-00922-0

**Reynoso MA, Kajala K, Bajic M, West DA, Pauluzzi G, Yao AI, Hatch K, Zumstein K, Woodhouse M, Rodriguez-Medina J, et al.** Evolutionary flexibility in flooding response circuitry in angiosperms. Science. 2019;365(6459):1291–1295. https://doi.org/10.1126/science.aax8862

**Rosenberg AB, Roco CM, Muscat RA, Kuchina A, Sample P, Yao Z, Graybuck LT, Peeler DJ, Mukherjee S, Chen W, et al.** Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding. Science. 2018;360(6385):176–182. https://doi.org/10.1126/science.aam8999

**Roszak P, Heo JO, Blob B, Toyokura K, Sugiyama Y, de Luis Balaguer MA, Lau WWY, Hamey F, Cirrone J, Madej E, et al.** Cell-by-cell dissection of phloem development links a maturation gradient to cell specialization. Science. 2021;374(6575):eaba5531. https://doi.org/10.1126/science.aba5531

**Ryu KH, Huang L, Kang HM, Schiefelbein J.** Single-cell RNA sequencing resolves molecular relationships among individual plant cells. Plant Physiol. 2019;179(4):1444–1456. https://doi.org/10.1104/pp.18.01482

**Saelens W, Cannoodt R, Todorov H, Saeys Y.** A comparison of single-cell trajectory inference methods. Nat Biotechnol. 2019;37(5):547–554. https://doi.org/10.1038/s41587-019-0071-9

**Satterlee JW, Strable J, Scanlon MJ.** Plant stem-cell organization and differentiation at single-cell resolution. Proc Natl Acad Sci U S A. 2020;117(52):33689–33699. https://doi.org/10.1073/pnas.201878117

**Serrano-Ron L, Perez-Garcia P, Sanchez-Corrieronero A, Gude I, Cabrera J, Ip P-L, Birnbaum KD, Moreno-Risueno MA.** Reconstruction of lateral root formation through single-cell RNA sequencing reveals order of tissue initiation. Mol Plant. 2021;14(8):1362–1378. https://doi.org/10.1016/j.molp.2021.05.028

**Shahan R, Hsu C-W, Nolan TM, Cole BJ, Taylor IW, Greenstreet L, Zhang S, Afanassiev A, Vlot AHC, Schiebinger G, et al.** A single-cell *Arabidopsis* root atlas reveals developmental trajectories in wild-type and cell identity mutants. Dev Cell. 2022;57(4):543–560.e9. https://doi.org/10.1016/j.devcel.2022.01.008

**Shen L, Ma J, Li P, Wu Y, Yu H.** Recent advances in the plant epitranscriptome. Genome Biol. 2023;24(1):43. https://doi.org/10.1186/s13059-023-02872-6

**Shulse CN, Cole BJ, Ciobanu D, Lin J, Yoshinaga Y, Gouran M, Turco GM, Zhu Y, O'Malley RC, Brady SM, et al.** High-throughput single-cell transcriptome profiling of plant cell types. Cell Rep. 2019;27(7):2241–2247.e4. https://doi.org/10.1016/j.celrep.2019.04.054

**Song Q, Ando A, Jiang N, Ikeda Y, Chen ZJ.** Single-cell RNA-seq analysis revealsddy-dependent and cell-specific transcriptome changes in *Arabidopsis* female gametophytes. Genome Biol. 2020;21(1):178. https://doi.org/10.1186/s13059-020-02094-0

**Squair JW, Gautier M, Kathe C, Anderson MA, James ND, Hutson TH, Hudelle R, Qaiser T, Matson KJE, Barraud Q, et al.** Confronting false discoveries in single-cell differential expression. Nat Commun. 2021;12(1):5692. https://doi.org/10.1038/s41467-021-25960-2

**Stickels RR, Murray E, Kumar P, Li J, Marshall JL, Di Bella DJ, Arlotta P, Macosko EZ, Chen F.** Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2. Nat Biotechnol. 2021;39(3):313–319. https://doi.org/10.1038/s41587-020-0739-1

**Stoeckius M, Zheng S, Houck-Loomis B, Hao S, Yeung BZ, Mauck WM III, Smibert P, Satija R.** Cell hashing with barcoded antibodies enables multiplexing and doublet detection for single cell genomics. Genome Biol. 2018;19(1):224. https://doi.org/10.1186/s13059-018-1603-1

**Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM III, Hao Y, Stoeckius M, Smibert P, Satija R.** Comprehensive integration of single-cell data. Cell. 2019;177(7):1888–1902.e21. https://doi.org/10.1016/j.cell.2019.05.031

**Sun G, Xia M, Li J, Ma W, Li Q, Xie J, Bai S, Fang S, Sun T, Feng X, et al.** The maize single-nucleus transcriptome comprehensively describes signaling networks governing movement and development of grass stomata. Plant Cell. 2022;34(5):1890–1911. https://doi.org/10.1093/plcell/koac047

**Sunaga-Franze DY, Muino JM, Braeuning C, Xu X, Zong M, Smaczniak C, Yan W, Fischer C, Vidal R, Kliem M, et al.** Single-nucleus RNA sequencing of plant tissues using a nanowell-based system. Plant J. 2021;108(3):859–869. https://doi.org/10.1111/tpj.15458

**Tan Y, Cahan P.** SingleCellNet: a computational tool to classify single cell RNA-seq data across platforms and across species. Cell Syst. 2019;9:207–213.e3. https://doi.org/10.1016/j.cels.2019.06.004

**Tao S, Liu P, Shi Y, Feng Y, Gao J, Chen L, Zhang A, Cheng X, Wei H, Zhang T, et al.** Single-cell transcriptome and network analyses unveil key transcription factors regulating mesophyll cell development in maize. Genes (Basel). 2022;13(2):374. https://doi.org/10.3390/genes13020374

**Tian C, Du Q, Xu M, Du F, Jiao Y.** Single-nucleus RNA-seq resolves spatiotemporal developmental trajectories in the tomato shoot apex. bioRxiv 305029. https://doi.org/10.1101/2020.09.20.305029, 20 September 2020, preprint: not peer reviewed.

**Tran HTN, Ang KS, Chevrier M, Zhang X, Lee NYS, Goh M, Chen J.** A benchmark of batch-effect correction methods for single-cell RNA sequencing data. Genome Biol. 2020;21(1):12. https://doi.org/10.1186/s13059-019-1850-9

**Tritschler S, Büttner M, Fischer DS, Lange M, Bergen V, Lickert H, Theis FJ.** Concepts and limitations for learning developmental trajectories from single cell genomics. Development. 2019;146(12):dev170506. https://doi.org/10.1242/dev.170506

**Turco GM, Rodriguez-Medina J, Siebert S, Han D, Valderrama-Gómez MÁ, Vahldick H, Shulse CN, Cole BJ, Juliano CE, Dickel DE, et al.** Molecular mechanisms driving switch behavior in xylem cell differentiation. Cell Rep. 2019;28(2):342–351.e4. https://doi.org/10.1016/j.celrep.2019.06.041

**Van den Berge K, Roux de Bézieux H, Street K, Saelens W, Cannoodt R, Saeys Y, Dudoit S, Clement L.** Trajectory-based differential expression analysis for single-cell sequencing data. Nat Commun. 2020;11(1):1201. https://doi.org/10.1038/s41467-020-14766-3
Wang Y, Huan Q, Li K, Qian W. Single-cell transcriptome atlas of the leaf and root of rice seedlings. J Genet Genomics. 2021:48(10): 881–898. https://doi.org/10.1016/j.jgg.2021.06.001

Wang X, Yu L, Wu AR. The effect of methanol fixation on single-cell RNA sequencing data. BMC Genomics. 2021:22(1):420. https://doi.org/10.1186/s12864-021-07744-6

Wendrich JR, Yang B, Vandamme N, Verstaen K, Smet W, Van de Velde C, Minne M, Wybouw B, Mor E, Arents HE, et al. Vascular transcription factors guide plant epidermal responses to limiting phosphate conditions. Science. 2020:370(6518):eaay4970. https://doi.org/10.1126/science.aay4970

Wilkinson MD, Dumontier M, Aalbersberg IJ, Appleton G, Axton M, Baak A, Blomberg N, Boiten J-W, da Silva Santos LB, Bourne PE, et al. The FAIR Guiding Principles for scientific data management and stewardship. Sci Data. 2016:3(1):160018. https://doi.org/10.1038/sdata.2016.18

Wohnhaas CT, Leparc GG, Fernandez-Albert F, Kind D, Gantner F, Viollet C, Hildebrandt T, Baum P. DMSO cryopreservation is the method of choice to preserve cells for droplet-based single-cell RNA sequencing. Sci Rep. 2019:9(1):10699. https://doi.org/10.1038/s41598-019-46932-z

Wolf FA, Angerer P, Theis FJ. SCANPY: large-scale single-cell gene expression data analysis. Genome Biol. 2018:19(1):15. https://doi.org/10.1186/s13059-017-1382-0

Wu Y, Tamayo P, Zhang K. Visualizing and interpreting single-cell gene expression datasets with similarity weighted nonnegative embedding. Cell Syst. 2018:7(6):656–666.e4. https://doi.org/10.1016/j.cels.2018.10.015

Xi NM, Li JJ. Benchmarking computational doublet-detection methods for single-cell RNA sequencing data. Cell Syst. 2021:12(2):176–194.e6. https://doi.org/10.1016/j.cels.2020.11.008

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

Xie B, Jiang Q, Mora A, Li X. Automatic cell type identification methods for single-cell RNA sequencing. Comput Struct Biotechnol J. 2021:19:5874–5887. https://doi.org/10.1016/j.csbj.2021.10.027

Xie J, Li M, Zeng J, Li X, Zhang D. Single-cell RNA sequencing profiles of stem-differentiating xylem in poplar. Plant Biotechnol J. 2022:20(3): 417–419. https://doi.org/10.1111/pbi.13763

Xu X, Crow M, Rice BR, Li F, Harris B, Liu L, Demesa-Arevalo E, Lu Z, Wang L, Fox N, et al. Single-cell RNA sequencing of developing maize ears facilitates functional analysis and trait candidate gene discovery. Dev Cell. 2021:56(4):557–568.e6. https://doi.org/10.1016/j.devcel. 2020.12.015

Yang B, Minne M, Brunoni F, Plačková L, Petřík I, Sun Y, Nolf J, Smet W, Verstaen K, Wendrich JR, et al. Non-cell autonomous and spatiotemporal signalling from a tissue organizer orchestrates root vascular development. Nat Plants. 2021:7(11):1485–1494. https:// doi.org/10.1038/s41477-021-01017-6

Yang X, Poelmans W, Grones C, Lakehal A, Pevernagie J, Van Bel M, Njo M, Xu L, Nelissen H, De et al. Spatial transcriptomics of a lycophyte root sheds light on root evolution. Curr Biol. 2023:33(19): 4069–4084.e8. https://doi.org/10.1016/j.cub.2023.08.030

Young MD, Behjati S. Soupx removes ambient RNA contamination from droplet-based single-cell RNA sequencing data. Gigascience. 2020:9(12):giaa151. https://doi.org/10.1093/gigascience/giaa151

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

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

Zhang MJ, Ntranos V, Tse D. Determining sequencing depth in a single-cell RNA-seq experiment. Nat Commun. 2020:11(1):774. https://doi.org/10.1038/s41467-020-14482-y

Zhang T-Q, Xu Z-G, Shang G-D, Wang J-W. A single-cell RNA sequencing profiles the developmental landscape of *Arabidopsis* root. Mol Plant. 2019:12(5):648–660. https://doi.org/10.1016/j.molp.2019.04.004

Zong J, Wang L, Zhu L, Bian L, Zhang B, Chen X, Huang G, Zhang X, Fan J, Cao L, et al. A rice single cell transcriptomic atlas defines the developmental trajectories of rice floret and inflorescence meristems. New Phytol. 2022:234(2):494–512. https://doi.org/10.1111/nph.18008
Wang Y, Huan Q, Li K, Qian W. Single-cell transcriptome atlas of the leaf and root of rice seedlings. J Genet Genomics. 2021:48(10): 881–898. https://doi.org/10.1016/j.jgg.2021.06.001

Wang X, Yu L, Wu AR. The effect of methanol fixation on single-cell RNA sequencing data. BMC Genomics. 2021:22(1):420. https://doi.org/10.1186/s12864-021-07744-6

Wendrich JR, Yang B, Vandamme N, Verstaen K, Smet W, Van de Velde C, Minne M, Wybouw B, Mor E, Arents HE, et al. Vascular transcription factors guide plant epidermal responses to limiting phosphate conditions. Science. 2020:370(6518):eaay4970. https://doi.org/10.1126/science.aay4970

Wilkinson MD, Dumontier M, Aalbersberg IJ, Appleton G, Axton M, Baak A, Blomberg N, Boiten J-W, da Silva Santos LB, Bourne PE, et al. The FAIR Guiding Principles for scientific data management and stewardship. Sci Data. 2016:3(1):160018. https://doi.org/10.1038/sdata.2016.18

Wohnhaas CT, Leparc GG, Fernandez-Albert F, Kind D, Gantner F, Viollet C, Hildebrandt T, Baum P. DMSO cryopreservation is the method of choice to preserve cells for droplet-based single-cell RNA sequencing. Sci Rep. 2019:9(1):10699. https://doi.org/10.1038/s41598-019-46932-z

Wolf FA, Angerer P, Theis FJ. SCANPY: large-scale single-cell gene expression data analysis. Genome Biol. 2018:19(1):15. https://doi.org/10.1186/s13059-017-1382-0

Wu Y, Tamayo P, Zhang K. Visualizing and interpreting single-cell gene expression datasets with similarity weighted nonnegative embedding. Cell Syst. 2018:7(6):656–666.e4. https://doi.org/10.1016/j.cels.2018.10.015

Xi NM, Li JJ. Benchmarking computational doublet-detection methods for single-cell RNA sequencing data. Cell Syst. 2021:12(2):176–194.e6. https://doi.org/10.1016/j.cels.2020.11.008

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

Xie B, Jiang Q, Mora A, Li X. Automatic cell type identification methods for single-cell RNA sequencing. Comput Struct Biotechnol J. 2021:19:5874–5887. https://doi.org/10.1016/j.csbj.2021.10.027

Xie J, Li M, Zeng J, Li X, Zhang D. Single-cell RNA sequencing profiles of stem-differentiating xylem in poplar. Plant Biotechnol J. 2022:20(3): 417–419. https://doi.org/10.1111/pbi.13763

Xu X, Crow M, Rice BR, Li F, Harris B, Liu L, Demesa-Arevalo E, Lu Z, Wang L, Fox N, et al. Single-cell RNA sequencing of developing maize ears facilitates functional analysis and trait candidate gene discovery. Dev Cell. 2021:56(4):557–568.e6. https://doi.org/10.1016/j.devcel. 2020.12.015

Yang B, Minne M, Brunoni F, Plačková L, Petřík I, Sun Y, Nolf J, Smet W, Verstaen K, Wendrich JR, et al. Non-cell autonomous and spatiotemporal signalling from a tissue organizer orchestrates root vascular development. Nat Plants. 2021:7(11):1485–1494. https:// doi.org/10.1038/s41477-021-01017-6

Yang X, Poelmans W, Grones C, Lakehal A, Pevernagie J, Van Bel M, Njo M, Xu L, Nelissen H, De et al. Spatial transcriptomics of a lycophyte root sheds light on root evolution. Curr Biol. 2023:33(19): 4069–4084.e8. https://doi.org/10.1016/j.cub.2023.08.030

Young MD, Behjati S. Soupx removes ambient RNA contamination from droplet-based single-cell RNA sequencing data. Gigascience. 2020:9(12):giaa151. https://doi.org/10.1093/gigascience/giaa151

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

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

Zhang MJ, Ntranos V, Tse D. Determining sequencing depth in a single-cell RNA-seq experiment. Nat Commun. 2020:11(1):774. https://doi.org/10.1038/s41467-020-14482-y

Zhang T-Q, Xu Z-G, Shang G-D, Wang J-W. A single-cell RNA sequencing profiles the developmental landscape of *Arabidopsis* root. Mol Plant. 2019:12(5):648–660. https://doi.org/10.1016/j.molp.2019.04.004

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