A Single-Cell RNA Sequencing Profiles the Developmental Landscape of Arabidopsis Root

MINI-EX: Integrative inference of single-cell gene regulatory networks in plants

Multicellular organisms, such as plants, are characterized by highly specialized and tightly regulated cell populations, establishing specific morphological structures and executing distinct functions. Gene regulatory networks (GRNs) describe condition-specific interactions of transcription factors (TFs) regulating the expression of target genes, underpinning these specific functions. As efficient and validated methods to identify cell-type-specific GRNs from single-cell data in plants are lacking, limiting our understanding of the organization of specific cell types in both model species and crops, we developed MINI-EX (Motif-Informed Network Inference based on single-cell EXpression data), an integrative approach to infer cell-type-specific networks in plants. MINI-EX uses single-cell transcriptomic data to define expression-based networks and integrates TF motif information to filter the inferred regulons, resulting in networks with increased accuracy. Next, regulons are assigned to different cell types, leveraging cell-specific expression, and candidate regulators are prioritized using network centrality measures, functional annotations, and expression specificity. This embedded prioritization strategy offers a unique and efficient means to unravel signaling cascades in specific cell types controlling a biological process of interest. We demonstrate the stability of MINI-EX toward input data sets with low number of cells and its robustness toward missing data, and show that it infers state-of-the-art networks with a better performance compared with other related single-cell network tools. MINI-EX successfully identifies key regulators controlling root development in Arabidopsis and rice, leaf development in Arabidopsis, and ear development in maize, enhancing our understanding of cell-type-specific regulation and unraveling the roles of different regulators controlling the development of specific cell types in plants.

Plant synthetic epigenomic engineering for crop improvement

Efforts have been directed to redesign crops with increased yield, stress adaptability, and nutritional value through synthetic biology-the application of engineering principles to biology. A recent expansion in our understanding of how epigenetic mechanisms regulate plant development and stress responses has unveiled a new set of resources that can be harnessed to develop improved crops, thus heralding the promise of "synthetic epigenetics." In this review, we summarize the latest advances in epigenetic regulation and highlight how innovative sequencing techniques, epigenetic editing, and deep learning-driven predictive tools can rapidly extend these insights. We also proposed the future directions of synthetic epigenetics for the development of engineered smart crops that can actively monitor and respond to internal and external cues throughout their life cycles.

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus

Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials - a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

AtHB40 modulates primary root length and gravitropism involving CYCLINB and auxin transporters

Gravitropism is a finely regulated tropistic response based on the plant perception of directional cues. Such perception allows them to direct shoot growth upwards, above ground, and root growth downwards, into the soil, anchoring the plant to acquire water and nutrients. Gravity sensing occurs in specialized cells and depends on auxin distribution, regulated by influx/efflux carriers. Here we report that AtHB40, encoding a transcription factor of the homeodomain-leucine zipper I family, was expressed in the columella and the root tip. Athb40 mutants exhibited longer primary roots. Enhanced primary root elongation was in agreement with a higher number of cells in the transition zone and the induction of CYCLINB transcript levels. Moreover, athb40 mutants and AtHB40 overexpressors displayed enhanced and delayed gravitropistic responses, respectively. These phenotypes were associated with altered auxin distribution and deregulated expression of the auxin transporters LAX2, LAX3, and PIN2. Accordingly, lax2 and lax3 mutants also showed an altered gravitropistic response, and LAX3 was identified as a direct target of AtHB40. Furthermore, AtHB40 is induced by AtHB53 when the latter is upregulated by auxin. Altogether, these results indicate that AtHB40 modulates cell division and auxin distribution in the root tip thus altering primary root length and gravitropism.

Single-cell transcriptome of Nepeta tenuifolia leaves reveal differentiation trajectories in glandular trichomes

The peltate glandular trichomes (PGTs) on Nepeta tenuifolia leaves can secrete and store bioactive essential oils. ScRNA-seq is a powerful tool for uncovering heterogeneous cells and exploring the development and differentiation of specific cells. Due to leaves rich in PGTs, the young leaves were used to isolated protoplasts and successfully captured 33,254 protoplasts for sequencing purposes. After cell type annotation, all the cells were partitioned into six broad populations with 19 clusters. Cells from PGTs were identified based on the expression patterns of trichome-specific genes, monoterpene biosynthetic genes, and metabolic analysis of PGT secretions. The developmental trajectories of PGTs were delineated by pseudotime analysis. Integrative analysis of scRNA-seq data from N. tenuifolia leaves and Arabidopsis thaliana shoot revealed that PGTs were specific to N. tenuifolia. Thus, our results provide a promising basis for exploring cell development and differentiation in plants, especially glandular trichome initiation and development.

Research strategies for single-cell transcriptome analysis in plant leaves

The recent and continuous improvement in single-cell RNA sequencing (scRNA-seq) technology has led to its emergence as an efficient experimental approach in plant research. However, compared with single-cell research in animals and humans, the application of scRNA-seq in plant research is limited by several challenges, including cell separation, cell type annotation, cellular function analysis, and cell-cell communication networks. In addition, the unavailability of corresponding reliable and stable analysis methods and standards has resulted in the relative decentralization of plant single-cell research. Considering these shortcomings, this review summarizes the research progress in plant leaf using scRNA-seq. In addition, it describes the corresponding feasible analytical methods and associated difficulties and problems encountered in the current research. In the end, we provide a speculative overview of the development of plant single-cell transcriptome research in the future.

Single-cell RNA-seq uncovered hemocyte functional subtypes and their differentiational characteristics and connectivity with morphological subpopulations in Litopenaeus vannamei

Hemocytes play central roles in shrimp immune system, whereas whose subclasses have not yet been completely defined. At present, the morphological classification of hemocytes is inadequate to classify the complete hemocyte repertoire and elucidate the functions and differentiation and maturation processes. Based on single-cell RNA sequencing (scRNA-seq) of hemocytes in healthy Litopenaeus vannamei, combined with RNA-FISH and flow cytometric sorting, we identified three hemocyte clusters including TGase+ cells, CTL+ cells and Crustin+ cells, and further determined their functional properties, potential differentiation trajectory and correspondence with morphological subpopulations. The TGase+ cells were mainly responsible for the coagulation, exhibiting distinguishable characteristics of hyalinocyte, and appeared to be developmentally arrested at an early stage of hemocyte differentiation. The CTL+ cells and Crustin+ cells arrested at terminal stages of differentiation mainly participated in recognizing foreign pathogens and initiating immune defense responses, owning distinctive features of granule-containing hemocytes. Furthermore, we have revealed the functional sub-clusters of three hemocyte clusters and their potential differentiation pathways according to the expression of genes involved in cell cycle, cell differentiation and immune response, and the successive differentiation and maturation of hyalinocytes to granule-containing hemocytes have also mapped. The results revealed the diversity of shrimp hemocytes and provide new theoretical rationale for hemocyte classification, which also facilitate systematic research on crustacean immunity.

Cell-Type-Specific Profiling of the Arabidopsis thaliana Membrane Protein-Encoding Genes

Membrane proteins work in large complexes to perceive and transduce external signals and to trigger a cellular response leading to the adaptation of the cells to their environment. Biochemical assays have been extensively used to reveal the interaction between membrane proteins. However, such analyses do not reveal the unique and complex composition of the membrane proteins of the different plant cell types. Here, we conducted a comprehensive analysis of the expression of Arabidopsis membrane proteins in the different cell types composing the root. Specifically, we analyzed the expression of genes encoding membrane proteins interacting in large complexes. We found that the transcriptional profiles of membrane protein-encoding genes differ between Arabidopsis root cell types. This result suggests that different cell types are characterized by specific sets of plasma membrane proteins, which are likely a reflection of their unique biological functions and interactions. To further explore the complexity of the Arabidopsis root cell membrane proteomes, we conducted a co-expression analysis of genes encoding interacting membrane proteins. This study confirmed previously reported interactions between membrane proteins, suggesting that the co-expression of genes at the single cell-type level can be used to support protein network predictions.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.

Root Cap to Soil Interface: A Driving Force Toward Plant Adaptation and Development

Land plants have developed robust roots to grow in diverse soil ecosystems. The distal end of the root tip has a specialized organ called the 'root cap'. The root cap assists the roots in penetrating the ground, absorbing water and minerals, avoiding heavy metals and regulating the rhizosphere microbiota. Furthermore, root-cap-derived auxin governs the lateral root patterning and directs root growth under varying soil conditions. The root cap formation is hypothesized as one of the key innovations during root evolution. Morphologically diversified root caps in early land plant lineage and later in angiosperms aid in improving the adaptation of roots and, thereby, plants in diverse soil environments. This review article presents a retrospective view of the root cap's important morphological and physiological characteristics for the root-soil interaction and their response toward various abiotic and biotic stimuli. Recent single-cell RNAseq data shed light on root cap cell-type-enriched genes. We compiled root cap cell-type-enriched genes from Arabidopsis, rice, maize and tomato and analyzed their transcription factor (TF) binding site enrichment. Further, the putative gene regulatory networks derived from root-cap-enriched genes and their TF regulators highlight the species-specific biological functions of root cap genes across the four plant species.

Transcriptional analysis of Ceratopteris richardii young sporophyte reveals conservation of stem cell factors in the root apical meristem

Gene expression in roots has been assessed in different plant species in studies ranging from complete organs to specific cell layers, and more recently at the single cell level. While certain genes or functional categories are expressed in the root of all or most plant species, lineage-specific genes have also been discovered. An increasing amount of transcriptomic data is available for angiosperms, while a limited amount of data is available for ferns, and few studies have focused on fern roots. Here, we present a de novo transcriptome assembly from three different parts of the Ceratopteris richardii young sporophyte. Differential gene expression analysis of the root tip transcriptional program showed an enrichment of functional categories related to histogenesis and cell division, indicating an active apical meristem. Analysis of a diverse set of orthologous genes revealed conserved expression in the root meristem, suggesting a preserved role for different developmental roles in this tissue, including stem cell maintenance. The reconstruction of evolutionary trajectories for ground tissue specification genes suggests a high degree of conservation in vascular plants, but not for genes involved in root cap development, showing that certain genes are absent in Ceratopteris or have intricate evolutionary paths difficult to track. Overall, our results suggest different processes of conservation and divergence of genes involved in root development.

A root phloem pole cell atlas reveals common transcriptional states in protophloem-adjacent cells

Single-cell sequencing has recently allowed the generation of exhaustive root cell atlases. However, some cell types are elusive and remain underrepresented. Here we use a second-generation single-cell approach, where we zoom in on the root transcriptome sorting with specific markers to profile the phloem poles at an unprecedented resolution. Our data highlight the similarities among the developmental trajectories and gene regulatory networks common to protophloem sieve element (PSE)-adjacent lineages in relation to PSE enucleation, a key event in phloem biology. As a signature for early PSE-adjacent lineages, we have identified a set of DNA-binding with one finger (DOF) transcription factors, the PINEAPPLEs (PAPL), that act downstream of PHLOEM EARLY DOF (PEAR) genes and are important to guarantee a proper root nutrition in the transition to autotrophy. Our data provide a holistic view of the phloem poles that act as a functional unit in root development.

Exploiting plant transcriptomic databases: Resources, tools, and approaches

There are now more than 300 000 RNA sequencing samples available, stemming from thousands of experiments capturing gene expression in organs, tissues, developmental stages, and experimental treatments for hundreds of plant species. The expression data have great value, as they can be re-analyzed by others to ask and answer questions that go beyond the aims of the study that generated the data. Because gene expression provides essential clues to where and when a gene is active, the data provide powerful tools for predicting gene function, and comparative analyses allow us to study plant evolution from a new perspective. This review describes how we can gain new knowledge from gene expression profiles, expression specificities, co-expression networks, differential gene expression, and experiment correlation. We also introduce and demonstrate databases that provide user-friendly access to these tools.

Transcriptional landscapes of de novo root regeneration from detached Arabidopsis leaves revealed by time-lapse and single-cell RNA sequencing analyses

Detached Arabidopsis thaliana leaves can regenerate adventitious roots, providing a platform for studying de novo root regeneration (DNRR). However, the comprehensive transcriptional framework of DNRR remains elusive. Here, we provide a high-resolution landscape of transcriptome reprogramming from wound response to root organogenesis in DNRR and show key factors involved in DNRR. Time-lapse RNA sequencing (RNA-seq) of the entire leaf within 12 h of leaf detachment revealed rapid activation of jasmonate, ethylene, and reactive oxygen species (ROS) pathways in response to wounding. Genetic analyses confirmed that ethylene and ROS may serve as wound signals to promote DNRR. Next, time-lapse RNA-seq within 5 d of leaf detachment revealed the activation of genes involved in organogenesis, wound-induced regeneration, and resource allocation in the wounded region of detached leaves during adventitious rooting. Genetic studies showed that BLADE-ON-PETIOLE1/2, which control aboveground organs, PLETHORA3/5/7, which control root organogenesis, and ETHYLENE RESPONSE FACTOR115, which controls wound-induced regeneration, are involved in DNRR. Furthermore, single-cell RNA-seq data revealed gene expression patterns in the wounded region of detached leaves during adventitious rooting. Overall, our study not only provides transcriptome tools but also reveals key factors involved in DNRR from detached Arabidopsis leaves.

PlantcircBase 7.0: Full-length transcripts and conservation of plant circRNAs

Circular RNA (circRNA) is a special type of non-coding RNA that participates in diverse biological processes in both animals and plants. Five years ago, we developed a comprehensive plant circRNA database (PlantcircBase), which has attracted much attention from the plant circRNA community. Here, we report an updated PlantcircBase (v.7.0), which contains 171,118 circRNAs from 21 plant species. Over 31,000 of the circRNAs have full-length sequences constructed based on analysis of 749 bulk RNA sequencing (RNA-seq) datasets downloaded from the public domain and Nanopore long-read sequencing results of rice RNAs newly generated in this study. A plant multiple conservation score (PMCS), based on the conservation of both sequence and expression profiles, was calculated for each circRNA to quantify and compare the conservation of all circRNAs. A new parameter, plant circRNA confidence level (PCCL), is introduced to measure the identity reliability of each circRNA based on experimental validation results and the number of references that support the circRNA. All this information and other details of circRNAs can be browsed, searched, and downloaded from PlantcircBase 7.0, which also provides online bioinformatics tools for visualization and sequence alignment. PlantcircBase 7.0 is publicly and freely accessible at http://ibi.zju.edu.cn/plantcircbase/.

A deep matrix factorization based approach for single-cell RNA-seq data clustering

The rapid development of single-cell sequencing technologies makes it possible to analyze cellular heterogeneity at the single-cell level. Cell clustering is one of the most fundamental and common steps in the heterogeneity analysis. However, due to the high noise level, high dimensionality and high sparsity, accurate cell clustering is still challengeable. Here, we present DeepCI, a new clustering approach for scRNA-seq data. Using two autoencoders to obtain cell embedding and gene embedding, DeepCI can simultaneously learn cell low-dimensional representation and clustering. In addition, the recovered gene expression matrix can be obtained by the matrix multiplication of cell and gene embedding. To evaluate the performance of DeepCI, we performed it on several real scRNA-seq datasets for clustering and visualization analysis. The experimental results show that DeepCI obtains the overall better performance than several popular single cell analysis methods. We also evaluated the imputation performance of DeepCI by a dedicated experiment. The corresponding results show that the imputed gene expression of known specific marker gene can greatly improve the accuracy of cell type classification.

REGENOMICS: A web-based application for plant REGENeration-associated transcriptOMICS analyses

In plants, differentiated somatic cells exhibit an exceptional ability to regenerate new tissues, organs, or whole plants. Recent studies have unveiled core genetic components and pathways underlying cellular reprogramming and de novo tissue regeneration in plants. Although high-throughput analyses have led to key discoveries in plant regeneration, a comprehensive organization of large-scale data is needed to further enhance our understanding of plant regeneration. Here, we collected all currently available transcriptome datasets related to wounding responses, callus formation, de novo organogenesis, somatic embryogenesis, and protoplast regeneration to construct REGENOMICS, a web-based application for plant REGENeration-associated transcriptOMICS analyses. REGENOMICS supports single- and multi-query analyses of plant regeneration-related gene-expression dynamics, co-expression networks, gene-regulatory networks, and single-cell expression profiles. Furthermore, it enables user-friendly transcriptome-level analysis of REGENOMICS-deposited and user-submitted RNA-seq datasets. Overall, we demonstrate that REGENOMICS can serve as a key hub of plant regeneration transcriptome analysis and greatly enhance our understanding on gene-expression networks, new molecular interactions, and the crosstalk between genetic pathways underlying each mode of plant regeneration. The REGENOMICS web-based application is available at http://plantregeneration.snu.ac.kr.

Advances and applications of single-cell omics technologies in plant research

Single-cell sequencing approaches reveal the intracellular dynamics of individual cells and answer biological questions with high-dimensional catalogs of millions of cells, including genomics, transcriptomics, chromatin accessibility, epigenomics, and proteomics data across species. These emerging yet thriving technologies have been fully embraced by the field of plant biology, with a constantly expanding portfolio of applications. Here, we introduce the current technical advances used for single-cell omics, especially single-cell genome and transcriptome sequencing. Firstly, we overview methods for protoplast and nucleus isolation and genome and transcriptome amplification. Subsequently, we use well-executed benchmarking studies to highlight advances made through the application of single-cell omics techniques. Looking forward, we offer a glimpse of additional hurdles and future opportunities that will introduce broad adoption of single-cell sequencing with revolutionary perspectives in plant biology.

scESI: evolutionary sparse imputation for single-cell transcriptomes from nearest neighbor cells

The ubiquitous dropout problem in single-cell RNA sequencing technology causes a large amount of data noise in the gene expression profile. For this reason, we propose an evolutionary sparse imputation (ESI) algorithm for single-cell transcriptomes, which constructs a sparse representation model based on gene regulation relationships between cells. To solve this model, we design an optimization framework based on nondominated sorting genetics. This framework takes into account the topological relationship between cells and the variety of gene expression to iteratively search the global optimal solution, thereby learning the Pareto optimal cell-cell affinity matrix. Finally, we use the learned sparse relationship model between cells to improve data quality and reduce data noise. In simulated datasets, scESI performed significantly better than benchmark methods with various metrics. By applying scESI to real scRNA-seq datasets, we discovered scESI can not only further classify the cell types and separate cells in visualization successfully but also improve the performance in reconstructing trajectories differentiation and identifying differentially expressed genes. In addition, scESI successfully recovered the expression trends of marker genes in stem cell differentiation and can discover new cell types and putative pathways regulating biological processes.

The molecular basis of zinc homeostasis in cereals

Plants require zinc (Zn) as an essential cofactor for diverse molecular, cellular and physiological functions. Zn is crucial for crop yield, but is one of the most limiting micronutrients in soils. Grasses like rice, wheat, maize and barley are crucial sources of food and nutrients for humans. Zn deficiency in these species therefore not only reduces annual yield but also directly results in Zn malnutrition of more than two billion people in the world. There has been good progress in understanding Zn homeostasis and Zn deficiency mechanisms in plants. However, our current knowledge of monocots, including grasses, remains insufficient. In this review, we provide a summary of our knowledge of molecular Zn homeostasis mechanisms in monocots, with a focus on important cereal crops. We additionally highlight divergences in Zn homeostasis of monocots and the dicot model Arabidopsis thaliana, as well as important gaps in our knowledge that need to be addressed in future research on Zn homeostasis in cereal monocots.

Identification of new marker genes from plant single-cell RNA-seq data using interpretable machine learning methods

An essential step in the analysis of single-cell RNA sequencing data is to classify cells into specific cell types using marker genes. In this study, we have developed a machine learning pipeline called single-cell predictive marker (SPmarker) to identify novel cell-type marker genes in the Arabidopsis root. Unlike traditional approaches, our method uses interpretable machine learning models to select marker genes. We have demonstrated that our method can: assign cell types based on cells that were labelled using published methods; project cell types identified by trajectory analysis from one data set to other data sets; and assign cell types based on internal GFP markers. Using SPmarker, we have identified hundreds of new marker genes that were not identified before. As compared to known marker genes, the new marker genes have more orthologous genes identifiable in the corresponding rice single-cell clusters. The new root hair marker genes also include 172 genes with orthologs expressed in root hair cells in five non-Arabidopsis species, which expands the number of marker genes for this cell type by 35-154%. Our results represent a new approach to identifying cell-type marker genes from scRNA-seq data and pave the way for cross-species mapping of scRNA-seq data in plants.

Single-Cell RNA Sequencing for Plant Research: Insights and Possible Benefits

In recent years, advances in single-cell RNA sequencing (scRNA-seq) technologies have continued to change our views on biological systems by increasing the spatiotemporal resolution of our analysis to single-cell resolution. Application of scRNA-seq to plants enables the comprehensive characterization of both common and rare cell types and cell states, uncovering new cell types and revealing how cell types relate to each other spatially and developmentally. This review provides an overview of scRNA-seq methodologies, highlights the application of scRNA-seq in plant science, justifies why scRNA-seq is a master player of sequencing, and explains the role of single-cell transcriptomics technologies in environmental stress adaptation, alongside the challenges and prospects of single-cell transcriptomics. Collectively, we put forward a central role of single-cell sequencing in plant research.

Applications of cell- and tissue-specific 'omics to improve plant productivity

The individual tissues and cell types of plants each have characteristic properties that contribute to the function of the plant as a whole. These are reflected by unique patterns of gene expression, protein and metabolite content, which enable cell-type-specific patterns of growth, development and physiology. Gene regulatory networks act within the cell types to govern the production and activity of these components. For the broader organism to grow and reproduce successfully, cell-type-specific activity must also function within the context of surrounding cell types, which is achieved by coordination of signalling pathways. We can investigate how gene regulatory networks are constructed and function using integrative 'omics technologies. Historically such experiments in plant biological research have been performed at the bulk tissue level, to organ resolution at best. In this review, we describe recent advances in cell- and tissue-specific 'omics technologies that allow investigation at much improved resolution. We discuss the advantages of these approaches for fundamental and translational plant biology, illustrated through the examples of specialised metabolism in medicinal plants and seed germination. We also discuss the challenges that must be overcome for such approaches to be adopted widely by the community.

A rice single cell transcriptomic atlas defines the developmental trajectories of rice floret and inflorescence meristems

Rice inflorescence development determines yield and relies on the activity of axillary meristems (AMs); however, high-resolution analysis of its early development is lacking. Here, we have used high-throughput single-cell RNA sequencing to profile 37 571 rice inflorescence cells and constructed a genome-scale gene expression resource covering the inflorescence-to-floret transition during early reproductive development. The differentiation trajectories of florets and AMs were reconstructed, and discrete cell types and groups of regulators in the highly heterogeneous young inflorescence were identified and then validated by in situ hybridization and with fluorescent marker lines. Our data demonstrate that a WOX transcription factor, DWARF TILLER1, regulates flower meristem activity, and provide evidence for the role of auxin in rice inflorescence branching by exploring the expression and biological role of the auxin importer OsAUX1. Our comprehensive transcriptomic atlas of early rice inflorescence development, supported by genetic evidence, provides single-cell-level insights into AM differentiation and floret development.

Identification of novel regulators required for early development of vein pattern in the cotyledons by single-cell RNA-sequencing

The leaf veins of higher plants contain a highly specialized vascular system comprised of xylem and phloem cells that transport water, organic compounds and mineral nutrients. The development of the vascular system is controlled by phytohormones that interact with complex transcriptional regulatory networks. Before the emergence of true leaves, the cotyledons of young seedlings perform photosynthesis that provides energy for the sustainable growth and survival of seedlings. However, the mechanisms underlying the early development of leaf veins in cotyledons are still not fully understood, in part due to the complex cellular composition of this tissue. To better understand the development of leaf veins, we analyzed 14 117 single cells from 3-day-old cotyledons using single-cell RNA sequencing. Based on gene expression patterns, we identified 10 clusters of cells and traced their developmental trajectories. We discovered multiple new marker genes and developmental features of leaf veins. The transcription factor networks of some cell types indicated potential roles of CYCLING DOF FACTOR 5 (CDF5) and REPRESSOR OF GA (RGA) in the early development and function of the leaf veins in cotyledons. These new findings lay a foundation for understanding the early developmental dynamics of cotyledon veins. The mechanisms underlying the early development of leaf veins in cotyledons are still not fully understood. In this study, we comprehensively characterized the early differentiation and development of leaf veins in 3-day-old cotyledons based on single-cell transcriptome analysis. We identified the cell types and novel marker genes of leaf veins and characterized the novel regulators of leaf vein.

Single-cell genomics revolutionizes plant development studies across scales

Understanding the development of tissues, organs and entire organisms through the lens of single-cell genomics has revolutionized developmental biology. Although single-cell transcriptomics has been pioneered in animal systems, from an experimental perspective, plant development holds some distinct advantages: cells do not migrate in relation to one another, and new organ formation (of leaves, roots, flowers, etc.) continues post-embryonically from persistent stem cell populations known as meristems. For a time, plant studies lagged behind animal or cell culture-based, single-cell approaches, largely owing to the difficulty in dissociating plant cells from their rigid cell walls. Recent intensive development of single-cell and single-nucleus isolation techniques across plant species has opened up a wide range of experimental approaches. This has produced a rapidly expanding diversity of information across tissue types and species, concomitant with the creative development of methods. In this brief Spotlight, we highlight some of the technical developments and how they have led to profiling single-cell genomics in various plant organs. We also emphasize the contribution of single-cell genomics in revealing developmental trajectories among different cell types within plant organs. Furthermore, we present efforts toward comparative analysis of tissues and organs at a single-cell level. Single-cell genomics is beginning to generate comprehensive information relating to how plant organs emerge from stem cell populations.

Protoplast Dissociation and Transcriptome Analysis Provides Insights to Salt Stress Response in Cotton

As one of the pioneer crops widely planted in saline-alkaline areas, Gossypium provides daily necessities, including natural fiber, vegetable proteins, and edible oils. However, cotton fiber yield and quality are highly influenced by salt stress. Therefore, elucidating the molecular mechanisms of cotton in response to salinity stress is importance to breed new cultivars with high tolerance. In this study, we first developed a method for single-cell RNA-seq based on isolating protoplast from cotton root tips; then, we studied the impact of salinity stress on gene expression profiling and their dynamic changes using the developed high-efficiency method for protoplast dissociation suitable for single-cell RNA-seq. A total of 3391 and 2826 differentially expressed genes (DEGs) were identified in salt-treated samples before and after protoplast dissociation, respectively, which were enriched into several molecular components, including response to stimulus, response to stress, and cellular macromolecule metabolic process by gene ontology (GO) analysis. Plant hormone signal transduction, phenylpropanoid biosynthesis, and MAPK signaling pathway were found to be enriched via Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Twenty-two and nine salinity-responsive DEGs participated in plant hormone signaling and MAPK signaling in roots, before and after protoplast dissociation, respectively; six upregulated DEGs were involved in ABA signaling transduction, namely, Ga04G2111, Ga07G0142, Ga09G2061, Ga10G0262, Ga01G0063, and Ga08G1915 which indicates their potential functions on plants adapting to salt stress. Additionally, 384 and 257 transcription factors (TFs) were differentially expressed in salt-stress roots before and after protoplast dissociation, respectively, of which significantly up-regulated TFs mainly belonged to the AP2/ERF-ERF family, which implied their potential roles responding to salt stress. These results not only provide novel insights to reveal the regulatory networks in plant’s root response to salt stress, but also lay the solid foundation for further exploration on cellular heterogeneity by single-cell transcriptome sequencing.

Single-cell RNA sequencing technologies and applications: A brief overview

Single-cell RNA sequencing (scRNA-seq) technology has become the state-of-the-art approach for unravelling the heterogeneity and complexity of RNA transcripts within individual cells, as well as revealing the composition of different cell types and functions within highly organized tissues/organs/organisms. Since its first discovery in 2009, studies based on scRNA-seq provide massive information across different fields making exciting new discoveries in better understanding the composition and interaction of cells within humans, model animals and plants. In this review, we provide a concise overview about the scRNA-seq technology, experimental and computational procedures for transforming the biological and molecular processes into computational and statistical data. We also provide an explanation of the key technological steps in implementing the technology. We highlight a few examples on how scRNA-seq can provide unique information for better understanding health and diseases. One important application of the scRNA-seq technology is to build a better and high-resolution catalogue of cells in all living organism, commonly known as atlas, which is key resource to better understand and provide a solution in treating diseases. While great promises have been demonstrated with the technology in all areas, we further highlight a few remaining challenges to be overcome and its great potentials in transforming current protocols in disease diagnosis and treatment.

Shaping the root system architecture in plants for adaptation to drought stress

Root system architecture plays an important role in plant adaptation to drought stress. The root system architecture (RSA) consists of several structural features, which includes number and length of main and lateral roots along with the density and length of root hairs. These features exhibit plasticity under water-limited environments and could be critical to developing crops with efficient root systems for adaptation under drought. Recent advances in the omics approaches have significantly improved our understanding of the regulatory mechanisms of RSA remodeling under drought and the identification of genes and other regulatory elements. Plant response to drought stress at physiological, morphological, biochemical, and molecular levels in root cells is regulated by various phytohormones and their crosstalk. Stress-induced reactive oxygen species play a significant role in regulating root growth and development under drought stress. Several transcription factors responsible for the regulation of RSA under drought have proven to be beneficial for developing drought tolerant crops. Molecular breeding programs for developing drought-tolerant crops have been greatly benefitted by the availability of quantitative trait loci (QTLs) associated with the RSA regulation. In the present review, we have discussed the role of various QTLs, signaling components, transcription factors, microRNAs and crosstalk among various phytohormones in shaping RSA and present future research directions to better understand various factors involved in RSA remodeling for adaptation to drought stress. We believe that the information provided herein may be helpful in devising strategies to develop crops with better RSA for efficient uptake and utilization of water and nutrients under drought conditions.

A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants

In all multicellular organisms, transcriptional networks orchestrate organ development. The Arabidopsis root, with its simple structure and indeterminate growth, is an ideal model for investigating the spatiotemporal transcriptional signatures underlying developmental trajectories. To map gene expression dynamics across root cell types and developmental time, we built a comprehensive, organ-scale atlas at single-cell resolution. In addition to estimating developmental progressions in pseudotime, we employed the mathematical concept of optimal transport to infer developmental trajectories and identify their underlying regulators. To demonstrate the utility of the atlas to interpret new datasets, we profiled mutants for two key transcriptional regulators at single-cell resolution, shortroot and scarecrow. We report transcriptomic and in vivo evidence for tissue trans-differentiation underlying a mixed cell identity phenotype in scarecrow. Our results support the atlas as a rich community resource for unraveling the transcriptional programs that specify and maintain cell identity to regulate spatiotemporal organ development.

A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants

In all multicellular organisms, transcriptional networks orchestrate organ development. The Arabidopsis root, with its simple structure and indeterminate growth, is an ideal model for investigating the spatiotemporal transcriptional signatures underlying developmental trajectories. To map gene expression dynamics across root cell types and developmental time, we built a comprehensive, organ-scale atlas at single-cell resolution. In addition to estimating developmental progressions in pseudotime, we employed the mathematical concept of optimal transport to infer developmental trajectories and identify their underlying regulators. To demonstrate the utility of the atlas to interpret new datasets, we profiled mutants for two key transcriptional regulators at single-cell resolution, shortroot and scarecrow. We report transcriptomic and in vivo evidence for tissue trans-differentiation underlying a mixed cell identity phenotype in scarecrow. Our results support the atlas as a rich community resource for unraveling the transcriptional programs that specify and maintain cell identity to regulate spatiotemporal organ development.

Single-Cell Transcriptome and Network Analyses Unveil Key Transcription Factors Regulating Mesophyll Cell Development in Maize

Background: Maize mesophyll (M) cells play important roles in various biological processes such as photosynthesis II and secondary metabolism. Functional differentiation occurs during M-cell development, but the underlying mechanisms for regulating M-cell development are largely unknown. Results: We conducted single-cell RNA sequencing (scRNA-seq) to profile transcripts in maize leaves. We then identified coregulated modules by analyzing the resulting pseudo-time-series data through gene regulatory network analyses. WRKY, ERF, NAC, MYB and Heat stress transcription factor (HSF) families were highly expressed in the early stage, whereas CONSTANS (CO)-like (COL) and ERF families were highly expressed in the late stage of M-cell development. Construction of regulatory networks revealed that these transcript factor (TF) families, especially HSF and COL, were the major players in the early and later stages of M-cell development, respectively. Integration of scRNA expression matrix with TF ChIP-seq and Hi-C further revealed regulatory interactions between these TFs and their targets. HSF1 and COL8 were primarily expressed in the leaf bases and tips, respectively, and their targets were validated with protoplast-based ChIP-qPCR, with the binding sites of HSF1 being experimentally confirmed. Conclusions: Our study provides evidence that several TF families, with the involvement of epigenetic regulation, play vital roles in the regulation of M-cell development in maize.

Single-cell genomics in plants: current state, future directions, and hurdles to overcome

Single-cell genomics has the potential to revolutionize the study of plant development and tissue-specific responses to environmental stimuli by revealing heretofore unknown players and gene regulatory processes. Here, I focus on the current state of single-cell genomics in plants, emerging technologies and applications, in addition to outlining possible future directions for experiments. I describe approaches to enable cheaper and larger experiments and technologies to measure multiple types of molecules to better model and understand cell types and their different states and trajectories throughout development. Lastly, I discuss the inherent limitations of single-cell studies and the technological hurdles that need to be overcome to widely apply single-cell genomics in crops to generate the greatest possible knowledge gain.

To the proteome and beyond: advances in single-cell omics profiling for plant systems

Recent advances in single-cell proteomics for animal systems could be adapted for plants to increase our understanding of plant development, response to stimuli, and cell-to-cell signaling.

Single-cell transcriptomics sheds light on the identity and metabolism of developing leaf cells

As the main photosynthetic instruments of vascular plants, leaves are crucial and complex plant organs. A strict organization of leaf mesophyll and epidermal cell layers orchestrates photosynthesis and gas exchange. In addition, water and nutrients for leaf growth are transported through the vascular tissue. To establish the single-cell transcriptomic landscape of these different leaf tissues, we performed high-throughput transcriptome sequencing of individual cells isolated from young leaves of Arabidopsis (Arabidopsis thaliana) seedlings grown in two different environmental conditions. The detection of approximately 19,000 different transcripts in over 1,800 high-quality leaf cells revealed 14 cell populations composing the young, differentiating leaf. Besides the cell populations comprising the core leaf tissues, we identified subpopulations with a distinct identity or metabolic activity. In addition, we proposed cell-type-specific markers for each of these populations. Finally, an intuitive web tool allows for browsing the presented dataset. Our data present insights on how the different cell populations constituting a developing leaf are connected via developmental, metabolic, or stress-related trajectories.

Shoot and root single cell sequencing reveals tissue- and daytime-specific transcriptome profiles

Although several large-scale single-cell RNA sequencing (scRNAseq) studies addressing the root of Arabidopsis (Arabidopsis thaliana) have been published, there is still need for a de novo reference map for both root and especially above-ground cell types. As the plants' transcriptome substantially changes throughout the day, shaped by the circadian clock, we performed scRNAseq on both Arabidopsis root and above-ground tissues at defined times of the day. For the root scRNAseq analysis, we used tissue-specific reporter lines grown on plates and harvested at the end of the day (ED). In addition, we submitted above-ground tissues from plants grown on soil at ED and end of the night to scRNAseq, which allowed us to identify common cell types/markers between root and shoot and uncover transcriptome changes to above-ground tissues depending on the time of the day. The dataset was also exploited beyond the traditional scRNAseq analysis to investigate non-annotated and di-cistronic transcripts. We experimentally confirmed the predicted presence of some of these transcripts and also addressed the potential function of a previously unidentified marker gene for dividing cells. In summary, this work provides insights into the spatial control of gene expression from nearly 70,000 cells of Arabidopsis for below- and whole above-ground tissue at single-cell resolution at defined time points.

Potassium transporter TRH1/KUP4 contributes to distinct auxin-mediated root system architecture responses

In plants, auxin transport and development are tightly coupled, just as hormone and growth responses are intimately linked in multicellular systems. Here we provide insights into uncoupling this tight control by specifically targeting the expression of TINY ROOT HAIR 1 (TRH1), a member of plant high-affinity potassium (K+)/K+ uptake/K+ transporter (HAK/KUP/KT) transporters that facilitate K+ uptake by co-transporting protons, in Arabidopsis root cell files. Use of this system pinpointed specific root developmental responses to acropetal versus basipetal auxin transport. Loss of TRH1 function shows TRHs and defective root gravitropism, associated with auxin imbalance in the root apex. Cell file-specific expression of TRH1 in the central cylinder rescued trh1 root agravitropism, whereas positional TRH1 expression in peripheral cell layers, including epidermis and cortex, restored trh1 defects. Applying a system-level approach, the role of RAP2.11 and ROOT HAIR DEFECTIVE-LIKE 5 transcription factors (TFs) in root hair development was verified. Furthermore, ERF53 and WRKY51 TFs were overrepresented upon restoration of root gravitropism supporting involvement in gravitropic control. Auxin has a central role in shaping root system architecture by regulating multiple developmental processes. We reveal that TRH1 jointly modulates intracellular ionic gradients and cell-to-cell polar auxin transport to drive root epidermal cell differentiation and gravitropic response. Our results indicate the developmental importance of HAK/KUP/KT proton-coupled K+ transporters.

The biology of time: dynamic responses of cell types to developmental, circadian and environmental cues

As sessile organisms, plants are finely tuned to respond dynamically to developmental, circadian and environmental cues. Genome-wide studies investigating these types of cues have uncovered the intrinsically different ways they can impact gene expression over time. Recent advances in single-cell sequencing and time-based bioinformatic algorithms are now beginning to reveal the dynamics of these time-based responses within individual cells and plant tissues. Here, we review what these techniques have revealed about the spatiotemporal nature of gene regulation, paying particular attention to the three distinct ways in which plant tissues are time sensitive. (i) First, we discuss how studying plant cell identity can reveal developmental trajectories hidden in pseudotime. (ii) Next, we present evidence that indicates that plant cell types keep their own local time through tissue-specific regulation of the circadian clock. (iii) Finally, we review what determines the speed of environmental signaling responses, and how they can be contingent on developmental and circadian time. By these means, this review sheds light on how these different scales of time-based responses can act with tissue and cell-type specificity to elicit changes in whole plant systems.

Integration of nutrient and water availabilities via auxin into the root developmental program

In most soils, the spatial distribution of nutrients and water in the rooting zone of plants is heterogeneous and changes over time. To access localized resources more efficiently, plants induce foraging responses by modulating individual morphological root traits, such as the length of the primary root or the number and length of lateral roots. These adaptive responses require the integration of exogenous and endogenous nutrient- or water-related signals into the root developmental program. Recent studies corroborated a central role of auxin in shaping root architectural traits in response to fluctuating nutrient and water availabilities. In this review, we highlight current knowledge on nutrient- and water-related developmental processes that impact root foraging and involve auxin as a central player. A deeper understanding and exploitation of these auxin-related processes and mechanisms promises advances in crop breeding for higher resource efficiency.

Advancing root developmental research through single-cell technologies

Single-cell RNA-sequencing has greatly increased the spatiotemporal resolution of root transcriptomics data, but we are still only scratching the surface of its full potential. Despite the challenges that remain in the field, the orderly aligned structure of the Arabidopsis root meristem makes it specifically suitable for lineage tracing and trajectory analysis. These methods will become even more potent by increasing resolution and specificity using tissue-specific single-cell RNA-sequencing and spatial transcriptomics. Feeding multiple single-cell omics data sets into single-cell gene regulatory networks will accelerate the discovery of regulators of root development in multiple species. By providing transcriptome atlases for virtually any species, single-cell technologies could tempt many root developmental biologists to move beyond the comfort of the well-known Arabidopsis root meristem.

Quality control and evaluation of plant epigenomics data

Epigenomics is the study of molecular signatures associated with discrete regions within genomes, many of which are important for a wide range of nuclear processes. The ability to profile the epigenomic landscape associated with genes, repetitive regions, transposons, transcription, differential expression, cis-regulatory elements, and 3D chromatin interactions has vastly improved our understanding of plant genomes. However, many epigenomic and single-cell genomic assays are challenging to perform in plants, leading to a wide range of data quality issues; thus, the data require rigorous evaluation prior to downstream analyses and interpretation. In this commentary, we provide considerations for the evaluation of plant epigenomics and single-cell genomics data quality with the aim of improving the quality and utility of studies using those data across diverse plant species.

Crafting a blueprint for single-cell RNA sequencing

Droplet-based single-cell RNA sequencing (scRNA-Seq) has moved rapidly from a technology with great potential to a method applied to ever-broader questions. The detailed information that scRNA-Seq offers has proven incredibly powerful in resolving cell responses to developmental and environmental cues. However, to maximize the potential of this technology, a panoply of upstream, practical points require consideration. Principal among these are the optimization of cell-isolation procedures, accommodating biotic/abiotic stress responses, and discerning the number of cells and sequencing reads needed. To complement excellent reviews outlining applications and data analysis tools for scRNA-Seq, we here discuss these considerations and provide practical tips to tailor experimental design and ensure the best possible outcome.

Identification of new proteins in mature sieve elements

The phloem enables vascular plants to transport photoassimilates from source tissues to heterotrophic sink tissues. In the phloem, unbroken strings of enucleated sieve elements, which lose the majority of their cellular contents upon maturation, provide a low resistance path for mass flow. The protein machinery in mature sieve elements performs vital functions to maintain the flow, transmit systemic signals and defend the sugar stream against pests. However, our knowledge of this particular protein population is very limited since mature sieve elements are difficult to isolate and not amenable to transcriptomic analysis due to their enucleate nature. Here, we used co-expression analysis and published gene clusters from transcriptomic studies to generate a list of sieve element proteins that potentially survive the enucleation process to reside in mature sieve elements. We selected seven candidates and show that they all localize in sieve elements in Arabidopsis roots and six of them in bolting stems. Our results support the idea that nascent sieve elements prior to enucleation translate part of the protein machinery found in mature sieve elements. Our co-expression list and the publicly available gene clusters expressed in late proto- and meta-phloem sieve elements are valuable resources for uncharacterized genes that may function in mature sieve elements.

The flowering transition pathways converge into a complex gene regulatory network that underlies the phase changes of the shoot apical meristem in Arabidopsis thaliana

Post-embryonic plant development is characterized by a period of vegetative growth during which a combination of intrinsic and extrinsic signals triggers the transition to the reproductive phase. To understand how different flowering inducing and repressing signals are associated with phase transitions of the Shoot Apical Meristem (SAM), we incorporated available data into a dynamic gene regulatory network model for Arabidopsis thaliana. This Flowering Transition Gene Regulatory Network (FT-GRN) formally constitutes a dynamic system-level mechanism based on more than three decades of experimental data on flowering. We provide novel experimental data on the regulatory interactions of one of its twenty-three components: a MADS-box transcription factor XAANTAL2 (XAL2). These data complement the information regarding flowering transition under short days and provides an example of the type of questions that can be addressed by the FT-GRN. The resulting FT-GRN is highly connected and integrates developmental, hormonal, and environmental signals that affect developmental transitions at the SAM. The FT-GRN is a dynamic multi-stable Boolean system, with 2²³ possible initial states, yet it converges into only 32 attractors. The latter are coherent with the expression profiles of the FT-GRN components that have been experimentally described for the developmental stages of the SAM. Furthermore, the attractors are also highly robust to initial states and to simulated perturbations of the interaction functions. The model recovered the meristem phenotypes of previously described single mutants. We also analyzed the attractors landscape that emerges from the postulated FT-GRN, uncovering which set of signals or components are critical for reproductive competence and the time-order transitions observed in the SAM. Finally, in the context of such GRN, the role of XAL2 under short-day conditions could be understood. Therefore, this model constitutes a robust biological module and the first multi-stable, dynamical systems biology mechanism that integrates the genetic flowering pathways to explain SAM phase transitions.

Cell-by-cell dissection of phloem development links a maturation gradient to cell specialization

In the plant meristem, tissue-wide maturation gradients are coordinated with specialized cell networks to establish various developmental phases required for indeterminate growth. Here, we used single-cell transcriptomics to reconstruct the protophloem developmental trajectory from the birth of cell progenitors to terminal differentiation in the Arabidopsis thaliana root. PHLOEM EARLY DNA-BINDING-WITH-ONE-FINGER (PEAR) transcription factors mediate lineage bifurcation by activating guanosine triphosphatase signaling and prime a transcriptional differentiation program. This program is initially repressed by a meristem-wide gradient of PLETHORA transcription factors. Only the dissipation of PLETHORA gradient permits activation of the differentiation program that involves mutual inhibition of early versus late meristem regulators. Thus, for phloem development, broad maturation gradients interface with cell-type-specific transcriptional regulators to stage cellular differentiation.

A Transcriptomic Atlas of the Ectomycorrhizal Fungus Laccaria bicolor

Trees are able to colonize, establish and survive in a wide range of soils through associations with ectomycorrhizal (EcM) fungi. Proper functioning of EcM fungi implies the differentiation of structures within the fungal colony. A symbiotic structure is dedicated to nutrient exchange and the extramatricular mycelium explores soil for nutrients. Eventually, basidiocarps develop to assure last stages of sexual reproduction. The aim of this study is to understand how an EcM fungus uses its gene set to support functional differentiation and development of specialized morphological structures. We examined the transcriptomes of Laccaria bicolor under a series of experimental setups, including the growth with Populus tremula x alba at different developmental stages, basidiocarps and free-living mycelium, under various conditions of N, P and C supply. In particular, N supply induced global transcriptional changes, whereas responses to P supply seemed to be independent from it. Symbiosis development with poplar is characterized by transcriptional waves. Basidiocarp development shares transcriptional signatures with other basidiomycetes. Overlaps in transcriptional responses of L. bicolor hyphae to a host plant and N/C supply next to co-regulation of genes in basidiocarps and mature mycorrhiza were detected. Few genes are induced in a single condition only, but functional and morphological differentiation rather involves fine tuning of larger gene sets. Overall, this transcriptomic atlas builds a reference to study the function and stability of EcM symbiosis in distinct conditions using L. bicolor as a model and indicates both similarities and differences with other ectomycorrhizal fungi, allowing researchers to distinguish conserved processes such as basidiocarp development from nutrient homeostasis.

Isolation of Plant Nuclei Compatible with Microfluidic Single-nucleus ATAC-sequencing

Gene expression depends on the binding of transcription factors with DNA regulatory sequences. The level of accessibility for these sequences varies between cells and cell types. Until recently, using the Tn5 assay for transposase-accessible chromatin for sequencing (ATAC-seq) technology allowed assessing the profiles of chromatin from an entire organ or, when coupled with the isolation of nuclei tagged in specific cell types (INTACT) method, from a cell-type. Recently, ATAC-seq experiments were conducted at the level of individual plant nuclei. Applying single nuclei ATAC-seq (sNucATAC-seq) technology to thousands of individual cells revealed more finely tuned profiles of chromatin accessibility. In this manuscript, we describe a method to isolate nuclei fom plant roots and green tissues, permeabilize the nuclear membrane using detergent to allow the penetration of the Tn5 transposase, and re-suspend them in a nuclei resuspension buffer compatible with the construction of sNucATAC-seq libraries using the 10× Genomic's Chromium technology. This protocol was successfully applied on Arabidopsis thaliana and Glycine max root nuclei.