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

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.

Cellular export of sugars and amino acids: role in feeding other cells and organisms

Sucrose, hexoses, and raffinose play key roles in the plant metabolism. Sucrose and raffinose, produced by photosynthesis, are translocated from leaves to flowers, developing seeds and roots. Translocation occurs in the sieve elements or sieve tubes of angiosperms. But how is sucrose loaded into and unloaded from the sieve elements? There seem to be two principal routes: one through plasmodesmata and one via the apoplasm. The best-studied transporters are the H+/SUCROSE TRANSPORTERs (SUTs) in the sieve element-companion cell complex. Sucrose is delivered to SUTs by SWEET sugar uniporters that release these key metabolites into the apoplasmic space. The H+/amino acid permeases and the UmamiT amino acid transporters are hypothesized to play analogous roles as the SUT-SWEET pair to transport amino acids. SWEETs and UmamiTs also act in many other important processes-for example, seed filling, nectar secretion, and pollen nutrition. We present information on cell type-specific enrichment of SWEET and UmamiT family members and propose several members to play redundant roles in the efflux of sucrose and amino acids across different cell types in the leaf. Pathogens hijack SWEETs and thus represent a major susceptibility of the plant. Here, we provide an update on the status of research on intercellular and long-distance translocation of key metabolites such as sucrose and amino acids, communication of the plants with the root microbiota via root exudates, discuss the existence of transporters for other important metabolites and provide potential perspectives that may direct future research activities.

Arabidopsis bioinformatics: tools and strategies

The sequencing of the Arabidopsis thaliana genome 21 years ago ushered in the genomics era for plant research. Since then, an incredible variety of bioinformatic tools permit easy access to large repositories of genomic, transcriptomic, proteomic, epigenomic and other '-omic' data. In this review, we cover some more recent tools (and highlight the 'classics') for exploring such data in order to help formulate quality, testable hypotheses, often without having to generate new experimental data. We cover tools for examining gene expression and co-expression patterns, undertaking promoter analyses and gene set enrichment analyses, and exploring protein-protein and protein-DNA interactions. We will touch on tools that integrate different data sets at the end of the article.

Transcriptional landscape of highly lignified poplar stems at single-cell resolution

BACKGROUND

Plant secondary growth depends on the activity of the vascular cambium, which produces xylem and phloem. Wood derived from xylem is the most abundant form of biomass globally and has played key socio-economic and subsistence roles throughout human history. However, despite intensive study of vascular development, the full diversity of cell types and the gene networks engaged are still poorly understood.

RESULTS

Here, we have applied an optimized protoplast isolation protocol and RNA sequencing to characterize the high-resolution single-cell transcriptional landscape of highly lignified poplar stems. We identify 20 putative cell clusters with a series of novel cluster-specific marker genes and find that these cells are highly heterogeneous based on the transcriptome. Analysis of these marker genes' expression dynamics enables reconstruction of the cell differentiation trajectories involved in phloem and xylem development. We find that different cell clusters exhibit distinct patterns of phytohormone responses and emphasize the use of our data to predict potential gene redundancy and identify candidate genes related to vascular development in trees.

CONCLUSIONS

These findings establish the transcriptional landscape of major cell types of poplar stems at single-cell resolution and provide a valuable resource for investigating basic principles of vascular cell specification and differentiation in trees.

Methods of In Situ Quantitative Root Biology

When dealing with plant roots, a multiscale description of the functional root structure is needed. Since the beginning of 21st century, new devices such as laser confocal microscopes have been accessible for coarse root structure measurements, including three-dimensional (3D) reconstruction. Most researchers are familiar with using simple 2D geometry visualization that does not allow quantitative determination of key morphological features from an organ-like perspective. We provide here a detailed description of the quantitative methods available for 3D analysis of root features at single-cell resolution, including root asymmetry, lateral root analysis, cell size and nuclear organization, cell-cycle kinetics, and chromatin structure analysis. Quantitative maps of the root apical meristem (RAM) are shown for different species, including Arabidopsis thaliana (L.), Heynh, Nicotiana tabacum L., Medicago sativa L., and Setaria italica (L.) P. Beauv. The 3D analysis of the RAM in these species showed divergence in chromatin organization and cell volume distribution that might be used to study root zonation for each root tissue. Detailed protocols and possible pitfalls in the usage of the marker lines are discussed. Therefore, researchers who need to improve their quantitative root biology portfolio can use them as a reference.

Plant PhysioSpace: a robust tool to compare stress response across plant species

Generalization of transcriptomics results can be achieved by comparison across experiments. This generalization is based on integration of interrelated transcriptomics studies into a compendium. Such a focus on the bigger picture enables both characterizations of the fate of an organism and distinction between generic and specific responses. Numerous methods for analyzing transcriptomics datasets exist. Yet, most of these methods focus on gene-wise dimension reduction to obtain marker genes and gene sets for, for example, pathway analysis. Relying only on isolated biological modules might result in missing important confounders and relevant contexts. We developed a method called Plant PhysioSpace, which enables researchers to compute experimental conditions across species and platforms without a priori reducing the reference information to specific gene sets. Plant PhysioSpace extracts physiologically relevant signatures from a reference dataset (i.e. a collection of public datasets) by integrating and transforming heterogeneous reference gene expression data into a set of physiology-specific patterns. New experimental data can be mapped to these patterns, resulting in similarity scores between the acquired data and the extracted compendium. Because of its robustness against platform bias and noise, Plant PhysioSpace can function as an inter-species or cross-platform similarity measure. We have demonstrated its success in translating stress responses between different species and platforms, including single-cell technologies. We have also implemented two R packages, one software and one data package, and a Shiny web application to facilitate access to our method and precomputed models.

Plant Metabolic Network 15: A resource of genome-wide metabolism databases for 126 plants and algae

To understand and engineer plant metabolism, we need a comprehensive and accurate annotation of all metabolic information across plant species. As a step towards this goal, we generated genome-scale metabolic pathway databases of 126 algal and plant genomes, ranging from model organisms to crops to medicinal plants (https://plantcyc.org). Of these, 104 have not been reported before. We systematically evaluated the quality of the databases, which revealed that our semi-automated validation pipeline dramatically improves the quality. We then compared the metabolic content across the 126 organisms using multiple correspondence analysis and found that Brassicaceae, Poaceae, and Chlorophyta appeared as metabolically distinct groups. To demonstrate the utility of this resource, we used recently published sorghum transcriptomics data to discover previously unreported trends of metabolism underlying drought tolerance. We also used single-cell transcriptomics data from the Arabidopsis root to infer cell type-specific metabolic pathways. This work shows the quality and quantity of our resource and demonstrates its wide-ranging utility in integrating metabolism with other areas of plant biology.

Stationary sieve element proteins

Vascular plants use the phloem to move sugars and other molecules from source leaves to sink organs such as roots and fruits. Within the phloem, enucleate sieve elements provide the low-resistance pipe system that enable bulk flow of sap. In this review, we provide an overview of the highly specific protein machinery that localize to mature sieve elements without entering the phloem translocation stream. Generally, the proteins either maintain the flow, protect the sieve element against pathogens or transmit system wide signals. A notable exception is found in poppy, where part of the opium biosynthesis is compartmentalized in sieve elements. Biosynthesis of sieve element proteins happens either continuously in companion cell or transiently in immature sieve elements before nuclear disintegration. The latter population is translated during differentiation and stays functional without turnover during the entire lifespan of sieve elements. We discuss how protein longevity imposes some interesting restrictions on plants, especially in arborescent monocots with long living sieve elements.

Single-cell RNA-seq describes the transcriptome landscape and identifies critical transcription factors in the leaf blade of the allotetraploid peanut (Arachis hypogaea L.)

Single-cell RNA-seq (scRNA-seq) has been highlighted as a powerful tool for the description of human cell transcriptome, but the technology has not been broadly applied in plant cells. Herein, we describe the successful development of a robust protoplast cell isolation system in the peanut leaf. A total of 6,815 single cells were divided into eight cell clusters based on reported marker genes by applying scRNA-seq. Further, a pseudo-time analysis was used to describe the developmental trajectory and interaction network of transcription factors (TFs) of distinct cell types during leaf growth. The trajectory enabled re-investigation of the primordium-driven development processes of the mesophyll and epidermis. These results suggest that palisade cells likely differentiate into spongy cells, while the epidermal cells originated earlier than the primordium. Subsequently, the developed method integrated multiple technologies to efficiently validate the scRNA-seq result in a homogenous cell population. The expression levels of several TFs were strongly correlated with epidermal ontogeny in accordance with obtained scRNA-seq values. Additionally, peanut AHL23 (AT-HOOK MOTIF NUCLEAR LOCALIZED PROTEIN 23), which is localized in nucleus, promoted leaf growth when ectopically expressed in Arabidopsis by modulating the phytohormone pathway. Together, our study displays that application of scRNA-seq can provide new hypotheses regarding cell differentiation in the leaf blade of Arachis hypogaea. We believe that this approach will enable significant advances in the functional study of leaf blade cells in the allotetraploid peanut and other plant species.

Single-nucleus RNA sequencing of plant tissues using a nanowell-based system

Single-cell genomics provides unprecedented potential for research on plant development and environmental responses. Here, we introduce a generic procedure for plant nucleus isolation combined with nanowell-based library preparation. Our method enables the transcriptome analysis of thousands of individual plant nuclei. It serves as an alternative to the use of protoplast isolation, which is currently a standard methodology for plant single-cell genomics, although it can be challenging for some plant tissues. We show the applicability of our nucleus isolation method by using different plant materials from different species. The potential of our single-nucleus RNA sequencing method is shown through the characterization of transcriptomes of seedlings and developing flowers from Arabidopsis thaliana. We evaluated the transcriptome dynamics during the early stages of anther development, identified stage-specific activities of transcription factors regulating this process, and predicted potential target genes of these transcription factors. Our nucleus isolation procedure can be applied in different plant species and tissues, thus expanding the toolkit for plant single-cell genomics experiments.

Single-cell RNA sequencing reveals a high-resolution cell atlas of xylem in Populus

High-throughput single-cell RNA sequencing (scRNA-seq) has advantages over traditional RNA-seq to explore spatiotemporal information on gene dynamic expressions in heterogenous tissues. We performed Drop-seq, a method for the dropwise sequestration of single cells for sequencing, on protoplasts from the differentiating xylem of Populus alba × Populus glandulosa. The scRNA-seq profiled 9,798 cells, which were grouped into 12 clusters. Through characterization of differentially expressed genes in each cluster and RNA in situ hybridizations, we identified vessel cells, fiber cells, ray parenchyma cells and xylem precursor cells. Diffusion pseudotime analyses revealed the differentiating trajectory of vessels, fiber cells and ray parenchyma cells and indicated a different differentiation process between vessels and fiber cells, and a similar differentiation process between fiber cells and ray parenchyma cells. We identified marker genes for each cell type (cluster) and key candidate regulators during developmental stages of xylem cell differentiation. Our study generates a high-resolution expression atlas of wood formation at the single cell level and provides valuable information on wood formation.

Single-nucleus RNA sequencing of plant tissues using a nanowell-based system

Single-cell genomics provides unprecedented potential for research on plant development and environmental responses. Here, we introduce a generic procedure for plant nucleus isolation combined with nanowell-based library preparation. Our method enables the transcriptome analysis of thousands of individual plant nuclei. It serves as an alternative to the use of protoplast isolation, which is currently a standard methodology for plant single-cell genomics, although it can be challenging for some plant tissues. We show the applicability of our nucleus isolation method by using different plant materials from different species. The potential of our single-nucleus RNA sequencing method is shown through the characterization of transcriptomes of seedlings and developing flowers from Arabidopsis thaliana. We evaluated the transcriptome dynamics during the early stages of anther development, identified stage-specific activities of transcription factors regulating this process, and predicted potential target genes of these transcription factors. Our nucleus isolation procedure can be applied in different plant species and tissues, thus expanding the toolkit for plant single-cell genomics experiments.

Maximizing the Utility of Transcriptomics Data in Inflammatory Skin Diseases

Inflammatory skin diseases are induced by disorders of the host defense system of the skin, which is composed of a barrier, innate and acquired immunity, as well as the cutaneous microbiome. These disorders are characterized by recurrent cutaneous lesions and intense itch, which seriously affecting life quality of people across all ages and ethnicities. To elucidate molecular factors for typical inflammatory skin diseases (such as psoriasis and atopic dermatitis), transcriptomic profiling assays have been largely performed. Additionally, single-cell RNA sequencing (scRNA-seq) as well as spatial transcriptomic profiling have revealed multiple potential translational targets and offered guides to improve diagnosis and treatment strategies for inflammatory skin diseases. High-throughput transcriptomics data has shown unprecedented power to disclose the complex pathophysiology of inflammatory skin diseases. Here, we will summarize discoveries from transcriptomics data and discuss how to maximize the transcriptomics data to propel the development of diagnostic biomarkers and therapeutic targets in inflammatory skin diseases.

PCMDB: a curated and comprehensive resource of plant cell markers

The advent of single-cell sequencing opened a new era in transcriptomic and genomic research. To understand cell composition using single-cell studies, a variety of cell markers have been widely used to label individual cell types. However, the specific database of cell markers for use by the plant research community remains very limited. To overcome this problem, we developed the Plant Cell Marker DataBase (PCMDB, http://www.tobaccodb.org/pcmdb/), which is based on a uniform annotation pipeline. By manually curating over 130 000 research publications, we collected a total of 81 117 cell marker genes of 263 cell types in 22 tissues across six plant species. Tissue- and cell-specific expression patterns can be visualized using multiple tools: eFP Browser, Bar, and UMAP/TSNE graph. The PCMDB also supports several analysis tools, including SCSA and SingleR, which allows for user annotation of cell types. To provide information about plant species currently unsupported in PCMDB, potential marker genes for other plant species can be searched based on homology with the supported species. PCMDB is a user-friendly hierarchical platform that contains five built-in search engines. We believe PCMDB will constitute a useful resource for researchers working on cell type annotation and the prediction of the biological function of individual cells.

Single-cell transcriptome atlas of the leaf and root of rice seedlings

As a multicellular organism, rice flourishes relying on gene expression diversity among cells of various functions. However, cellular-resolution transcriptome features are yet to be fully recognized, let alone cell-specific transcriptional responses to environmental stimuli. In this study, we apply single-cell RNA sequencing to both shoot and root of rice seedlings growing in Kimura B nutrient solution or exposed to various abiotic stresses and characterize transcriptomes for a total of 237,431 individual cells. We identify 15 and 9 cell types in the leaf and root, respectively, and observe that common transcriptome features are often shared between leaves and roots in the same tissue layer, except for endodermis or epidermis. Abiotic stress stimuli alter gene expression largely in a cell type-specific manner, but for a given cell type, different stresses often trigger transcriptional regulation of roughly the same set of genes. Besides, we detect proportional changes in cell populations in response to abiotic stress and investigate the underlying molecular mechanisms through single-cell reconstruction of the developmental trajectory. Collectively, our study represents a benchmark-setting data resource of single-cell transcriptome atlas for rice seedlings and an illustration of exploiting such resources to drive discoveries in plant biology.

Effects of Sample Size on Plant Single-Cell RNA Profiling

Single-cell RNA (scRNA) profiling or scRNA-sequencing (scRNA-seq) makes it possible to parallelly investigate diverse molecular features of multiple types of cells in a given plant tissue and discover cell developmental processes. In this study, we evaluated the effects of sample size (i.e., cell number) on the outcome of single-cell transcriptome analysis by sampling different numbers of cells from a pool of ~57,000 Arabidopsis thaliana root cells integrated from five published studies. Our results indicated that the most significant principal components could be achieved when 20,000–30,000 cells were sampled, a relatively high reliability of cell clustering could be achieved by using ~20,000 cells with little further improvement by using more cells, 96% of the differentially expressed genes could be successfully identified with no more than 20,000 cells, and a relatively stable pseudotime could be estimated in the subsample with 5000 cells. Finally, our results provide a general guide for optimizing sample size to be used in plant scRNA-seq studies.

A perspective on cell proliferation kinetics in the root apical meristem

Organogenesis in plants is primarily postembryonic and relies on a strict balance between cell division and cell expansion. The root is a particularly well-suited model to study cell proliferation in detail since the two processes are spatially and temporally separated for all the different tissues. In addition, the root is amenable to detailed microscopic analysis to identify cells progressing through the cell cycle. While it is clear that cell proliferation activity is restricted to the root apical meristem (RAM), understanding cell proliferation kinetics and identifying its parameters have required much effort over many years. Here, we review the main concepts, experimental settings, and findings aimed at obtaining a detailed knowledge of how cells proliferate within the RAM. The combination of novel tools, experimental strategies, and mathematical models has contributed to our current view of cell proliferation in the RAM. We also discuss several lines of research that need to be explored in the future.

Single-cell analysis of cell identity in the Arabidopsis root apical meristem: insights and opportunities

A fundamental question in developmental biology is how the progeny of stem cells become differentiated tissues. The Arabidopsis root is a tractable model to address this question due to its simple organization and defined cell lineages. In particular, the zone of dividing cells at the root tip-the root apical meristem-presents an opportunity to map the gene regulatory networks underlying stem cell niche maintenance, tissue patterning, and cell identity acquisition. To identify molecular regulators of these processes, studies over the last 20 years employed global profiling of gene expression patterns. However, these technologies are prone to information loss due to averaging gene expression signatures over multiple cell types and/or developmental stages. Recently developed high-throughput methods to profile gene expression at single-cell resolution have been successfully applied to plants. Here, we review insights from the first published single-cell mRNA sequencing and chromatin accessibility datasets generated from Arabidopsis roots. These studies successfully reconstruct developmental trajectories, phenotype cell identity mutants at unprecedented resolution, and reveal cell type-specific responses to environmental stimuli. The experimental insight gained from Arabidopsis paves the way to profile roots from additional species.

Potentials of single-cell genomics in deciphering cellular phenotypes

Single-cell genomics, particularly single-cell transcriptome profiling by RNA sequencing have transformed the possibilities to relate genes to functions, structures, and eventually phenotypes. We can now observe changes in each cell's transcriptome and among its neighborhoods, interrogate the sequence of transcriptional events, and assess their influence on subsequent events. This paradigm shift in biology enables us to infer causal relationships in these events with high accuracy. Here we review the latest single-cell studies in plants that uncover how cellular phenotypes emerge as a result of the transcriptome process such as waves of expression, trajectories of development and responses to the environment, and spatial information. With an eye on the advances made in animal and human studies, we further highlight some of the needed areas for future research and development, including computational methods.

High-throughput single-cell transcriptomics reveals the female germline differentiation trajectory in Arabidopsis thaliana

Female germline cells in flowering plants differentiate from somatic cells to produce specialized reproductive organs, called ovules, embedded deep inside the flowers. We investigated the molecular basis of this distinctive developmental program by performing single-cell RNA sequencing (scRNA-seq) of 16,872 single cells of Arabidopsis thaliana ovule primordia at three developmental time points during female germline differentiation. This allowed us to identify the characteristic expression patterns of the main cell types, including the female germline and its surrounding nucellus. We then reconstructed the continuous trajectory of female germline differentiation and observed dynamic waves of gene expression along the developmental trajectory. A focused analysis revealed transcriptional cascades and identified key transcriptional factors that showed distinct expression patterns along the germline differentiation trajectory. Our study provides a valuable reference dataset of the transcriptional process during female germline differentiation at single-cell resolution, shedding light on the mechanisms underlying germline cell fate determination.

Zhimin Hou, Yanhui Liu et al. used single cell RNA-seq to analyze the model organism, Arabidopsis thaliana, at three stages during female germline differentiation. They reconstructed the continuous trajectory of female germline differentiation, providing a valuable reference for future investigation of germline cell fate determination in A. thaliana.

Plant Cell Identity in the Era of Single-Cell Transcriptomics

High-throughput single-cell transcriptomic approaches have revolutionized our view of gene expression at the level of individual cells, providing new insights into their heterogeneity, identities, and functions. Recently, technical challenges to the application of single-cell transcriptomics to plants have been overcome, and many plant organs and tissues have now been subjected to analyses at single-cell resolution. In this review, we describe these studies and their impact on our understanding of the diversity, differentiation, and activities of plant cells. We particularly highlight their impact on plant cell identity, including unprecedented views of cell transitions and definitions of rare and novel cell types. We also point out current challenges and future opportunities for the application and analyses of single-cell transcriptomics in plants. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Epigenome plasticity in plants

Plant intra-individual and inter-individual variation can be determined by the epigenome, a set of covalent modifications of DNA and chromatin that can alter genome structure and activity without changes to the genome sequence. The epigenome of plant cells is plastic, that is, it can change in response to internal or external cues, such as during development or due to environmental changes, to create a memory of such events. Ongoing advances in technologies to read and write epigenomic patterns with increasing resolution, scale and precision are enabling the extent of plant epigenome variation to be more extensively characterized and functionally interrogated. In this Review, we discuss epigenome dynamics and variation within plants during development and in response to environmental changes, including stress, as well as between plants. We review known or potential functions of such plasticity and emphasize the importance of investigating the causality of epigenomic changes. Finally, we discuss emerging technologies that may underpin future research into plant epigenome plasticity.

Vision, challenges and opportunities for a Plant Cell Atlas

With growing populations and pressing environmental problems, future economies will be increasingly plant-based. Now is the time to reimagine plant science as a critical component of fundamental science, agriculture, environmental stewardship, energy, technology and healthcare. This effort requires a conceptual and technological framework to identify and map all cell types, and to comprehensively annotate the localization and organization of molecules at cellular and tissue levels. This framework, called the Plant Cell Atlas (PCA), will be critical for understanding and engineering plant development, physiology and environmental responses. A workshop was convened to discuss the purpose and utility of such an initiative, resulting in a roadmap that acknowledges the current knowledge gaps and technical challenges, and underscores how the PCA initiative can help to overcome them.

Plant single-cell solutions for energy and the environment

Progress in sequencing, microfluidics, and analysis strategies has revolutionized the granularity at which multicellular organisms can be studied. In particular, single-cell transcriptomics has led to fundamental new insights into animal biology, such as the discovery of new cell types and cell type-specific disease processes. However, the application of single-cell approaches to plants, fungi, algae, or bacteria (environmental organisms) has been far more limited, largely due to the challenges posed by polysaccharide walls surrounding these species' cells. In this perspective, we discuss opportunities afforded by single-cell technologies for energy and environmental science and grand challenges that must be tackled to apply these approaches to plants, fungi and algae. We highlight the need to develop better and more comprehensive single-cell technologies, analysis and visualization tools, and tissue preparation methods. We advocate for the creation of a centralized, open-access database to house plant single-cell data. Finally, we consider how such efforts should balance the need for deep characterization of select model species while still capturing the diversity in the plant kingdom. Investments into the development of methods, their application to relevant species, and the creation of resources to support data dissemination will enable groundbreaking insights to propel energy and environmental science forward.

Enhancing Our Understanding of Plant Cell-to-Cell Interactions Using Single-Cell Omics

Plants are composed of cells that physically interact and constantly adapt to their environment. To reveal the contribution of each plant cells to the biology of the entire organism, their molecular, morphological, and physiological attributes must be quantified and analyzed in the context of the morphology of the plant organs. The emergence of single-cell/nucleus omics technologies now allows plant biologists to access different modalities of individual cells including their epigenome and transcriptome to reveal the unique molecular properties of each cell composing the plant and their dynamic regulation during cell differentiation and in response to their environment. In this manuscript, we provide a perspective regarding the challenges and strategies to collect plant single-cell biological datasets and their analysis in the context of cellular interactions. As an example, we provide an analysis of the transcriptional regulation of the Arabidopsis genes controlling the differentiation of the root hair cells at the single-cell level. We also discuss the perspective of the use of spatial profiling to complement existing plant single-cell omics.

Reconstruction of lateral root formation through single-cell RNA sequencing reveals order of tissue initiation

Postembryonic organogenesis is critical for plant development. Underground, lateral roots (LRs) form the bulk of mature root systems, yet the ontogeny of the LR primordium (LRP) is not clear. In this study, we performed the single-cell RNA sequencing through the first four stages of LR formation in Arabidopsis. Our analysis led to a model in which a single group of precursor cells, with a cell identity different from their pericycle origins, rapidly reprograms and splits into a mixed ground tissue/stem cell niche fate and a vascular precursor fate. The ground tissue and stem cell niche fates soon separate and a subset of more specialized vascular cells form sucrose transporting phloem cells that appear to connect to the primary root. We did not detect cells resembling epidermis or root cap, suggesting that outer tissues may form later, preceding LR emergence. At this stage, some remaining initial precursor cells form the primordium flanks, while the rest create a reservoir of pluripotent cells that are able to replace the LR if damaged. Laser ablation of the central and lateral LRP regions showed that remaining cells restart the sequence of tissue initiation to form a LR. Collectively, our study reveals an ontological hierarchy for LR formation with an early and sequential split of main root tissues and stem cells.

Network biology to uncover functional and structural properties of the plant immune system

In the last two decades, advances in network science have facilitated the discovery of important systems' entities in diverse biological networks. This graph-based technique has revealed numerous emergent properties of a system that enable us to understand several complex biological processes including plant immune systems. With the accumulation of multiomics data sets, the comprehensive understanding of plant-pathogen interactions can be achieved through the analyses and efficacious integration of multidimensional qualitative and quantitative relationships among the components of hosts and their microbes. This review highlights comparative network topology analyses in plant-pathogen co-expression networks and interactomes, outlines dynamic network modeling for cell-specific immune regulatory networks, and discusses the new frontiers of single-cell sequencing as well as multiomics data integration that are necessary for unraveling the intricacies of plant immune systems.

RNA-seq analysis of laser microdissected Arabidopsis thaliana leaf epidermis, mesophyll and vasculature defines tissue-specific transcriptional responses to multiple stress treatments

Acclimation of plants to adverse conditions requires the coordination of gene expression and signalling pathways between tissues and cell types. As the energy and carbon capturing organs, leaves are significantly affected by abiotic and biotic stresses. However, tissue- or cell type-specific analyses of stress responses have focussed on the Arabidopsis root. Here, we comparatively explore the transcriptomes of three leaf tissues (epidermis, mesophyll, vasculature) after induction of diverse stress pathways by chemical stimuli (antimycin A, 3-amino-1,2,4-triazole, methyl viologen, salicylic acid) and ultraviolet light in Arabidopsis using laser capture microdissection followed by RNA sequencing. Stimulation of stress pathways caused an overall reduction in the number of genes expressed in a tissue-specific manner, though a small subset gained or changed their tissue specificity. We find no evidence of a common stress response, with only a few genes consistently responsive to two or more treatments in the analysed tissues. However, differentially expressed genes overlap between tissues for individual treatments. A focussed analysis provided evidence for an interaction of auxin and ethylene that mediates retrograde signalling during mitochondrial dysfunction specifically in the epidermis, and a gene regulatory network defined the hierarchy of interactions. Taken together, we have generated an extensive reference dataset that will be valuable for future experiments analysing transcriptional responses on a tissue or single-cell level. Our results will enable the tailoring of the tissue-specific engineering of stress-tolerant plants.

Gene expression variation in Arabidopsis embryos at single-nucleus resolution

Soon after fertilization of egg and sperm, plant genomes become transcriptionally activated and drive a series of coordinated cell divisions to form the basic body plan during embryogenesis. Early embryonic cells rapidly diversify from each other, and investigation of the corresponding gene expression dynamics can help elucidate underlying cellular differentiation programs. However, current plant embryonic transcriptome datasets either lack cell-specific information or have RNA contamination from surrounding non-embryonic tissues. We have coupled fluorescence-activated nuclei sorting together with single-nucleus mRNA-sequencing to construct a gene expression atlas of Arabidopsis thaliana early embryos at single-cell resolution. In addition to characterizing cell-specific transcriptomes, we found evidence that distinct epigenetic and transcriptional regulatory mechanisms operate across emerging embryonic cell types. These datasets and analyses, as well as the approach we devised, are expected to facilitate the discovery of molecular mechanisms underlying pattern formation in plant embryos. This article has an associated 'The people behind the papers' interview.

Glandular trichomes: new focus on horticultural crops

Plant glandular trichomes (GTs) are epidermal outgrowths with the capacity to biosynthesize and secrete specialized metabolites, that are of great scientific and practical significance. Our understanding of the developmental process of GTs is limited, and no single plant species serves as a unique model. Here, we review the genetic mechanisms of GT initiation and development and provide a summary of the biosynthetic pathways of GT-specialized metabolites in nonmodel plant species, especially horticultural crops. We discuss the morphology and classification of GT types. Moreover, we highlight technological advancements in methods employed for investigating GTs. Understanding the molecular basis of GT development and specialized metabolites not only offers useful avenues for research in plant breeding that will lead to the improved production of desirable metabolites, but also provides insights for plant epidermal development research.

Single cell gene regulatory networks in plants: Opportunities for enhancing climate change stress resilience

Global warming poses major challenges for plant survival and agricultural productivity. Thus, efforts to enhance stress resilience in plants are key strategies for protecting food security. Gene regulatory networks (GRNs) are a critical mechanism conferring stress resilience. Until recently, predicting GRNs of the individual cells that make up plants and other multicellular organisms was impeded by aggregate population scale measurements of transcriptome and other genome-scale features. With the advancement of high-throughput single cell RNA-seq and other single cell assays, learning GRNs for individual cells is now possible, in principle. In this article, we report on recent advances in experimental and analytical methodologies for single cell sequencing assays especially as they have been applied to the study of plants. We highlight recent advances and ongoing challenges for scGRN prediction, and finally, we highlight the opportunity to use scGRN discovery for studying and ultimately enhancing abiotic stress resilience in plants.

Genome-wide expression and network analyses of mutants in key brassinosteroid signaling genes

BACKGROUND

Brassinosteroid (BR) signaling regulates plant growth and development in concert with other signaling pathways. Although many genes have been identified that play a role in BR signaling, the biological and functional consequences of disrupting those key BR genes still require detailed investigation.

RESULTS

Here we performed phenotypic and transcriptomic comparisons of A. thaliana lines carrying a loss-of-function mutation in BRI1 gene, bri1-5, that exhibits a dwarf phenotype and its three activation-tag suppressor lines that were able to partially revert the bri1-5 mutant phenotype to a WS2 phenotype, namely bri1-5/bri1-1D, bri1-5/brs1-1D, and bri1-5/bak1-1D. From the three investigated bri1-5 suppressors, bri1-5/bak1-1D was the most effective suppressor at the transcriptional level. All three bri1-5 suppressors showed altered expression of the genes in the abscisic acid (ABA signaling) pathway, indicating that ABA likely contributes to the partial recovery of the wild-type phenotype in these bri1-5 suppressors. Network analysis revealed crosstalk between BR and other phytohormone signaling pathways, suggesting that interference with one hormone signaling pathway affects other hormone signaling pathways. In addition, differential expression analysis suggested the existence of a strong negative feedback from BR signaling on BR biosynthesis and also predicted that BRS1, rather than being directly involved in signaling, might be responsible for providing an optimal environment for the interaction between BRI1 and its ligand.

CONCLUSIONS

Our study provides insights into the molecular mechanisms and functions of key brassinosteroid (BR) signaling genes, especially BRS1.

Time-Based Systems Biology Approaches to Capture and Model Dynamic Gene Regulatory Networks

All aspects of transcription and its regulation involve dynamic events. However, capturing these dynamic events in gene regulatory networks (GRNs) offers both a promise and a challenge. The promise is that capturing and modeling the dynamic changes in GRNs will allow us to understand how organisms adapt to a changing environment. The ability to mount a rapid transcriptional response to environmental changes is especially important in nonmotile organisms such as plants. The challenge is to capture these dynamic, genome-wide events and model them in GRNs. In this review, we cover recent progress in capturing dynamic interactions of transcription factors with their targets-at both the local and genome-wide levels-and how they are used to learn how GRNs operate as a function of time. We also discuss recent advances that employ time-based machine learning approaches to forecast gene expression at future time points, a key goal of systems biology.

Advances and Opportunities in Single-Cell Transcriptomics for Plant Research

Single-cell approaches are quickly changing our view on biological systems by increasing the spatiotemporal resolution of our analyses to the level of the individual cell. The field of plant biology has fully embraced single-cell transcriptomics and is rapidly expanding the portfolio of available technologies and applications. In this review, we give an overview of the main advances in plant single-cell transcriptomics over the past few years and provide the reader with an accessible guideline covering all steps, from sample preparation to data analysis. We end by offering a glimpse of how these technologies will shape and accelerate plant-specific research in the near future.

Potential Networks of Nitrogen-Phosphorus-Potassium Channels and Transporters in Arabidopsis Roots at a Single Cell Resolution

Nitrogen (N), phosphorus (P), and potassium (K) are three major macronutrients essential for plant life. These nutrients are acquired and transported by several large families of transporters expressed in plant roots. However, it remains largely unknown how these transporters are distributed in different cell-types that work together to transfer the nutrients from the soil to different layers of root cells and eventually reach vasculature for massive flow. Using the single cell transcriptomics data from Arabidopsis roots, we profiled the transcriptional patterns of putative nutrient transporters in different root cell-types. Such analyses identified a number of uncharacterized NPK transporters expressed in the root epidermis to mediate NPK uptake and distribution to the adjacent cells. Some transport genes showed cortex- and endodermis-specific expression to direct the nutrient flow toward the vasculature. For long-distance transport, a variety of transporters were shown to express and potentially function in the xylem and phloem. In the context of subcellular distribution of mineral nutrients, the NPK transporters at subcellular compartments were often found to show ubiquitous expression patterns, which suggests function in house-keeping processes. Overall, these single cell transcriptomic analyses provide working models of nutrient transport from the epidermis across the cortex to the vasculature, which can be further tested experimentally in the future.

The lncRNA APOLO interacts with the transcription factor WRKY42 to trigger root hair cell expansion in response to cold

Plant long noncoding RNAs (lncRNAs) have emerged as important regulators of chromatin dynamics, impacting on transcriptional programs leading to different developmental outputs. The lncRNA AUXIN-REGULATED PROMOTER LOOP (APOLO) directly recognizes multiple independent loci across the Arabidopsis genome and modulates their three-dimensional chromatin conformation, leading to transcriptional shifts. Here, we show that APOLO recognizes the locus encoding the root hair (RH) master regulator ROOT HAIR DEFECTIVE 6 (RHD6) and controls RHD6 transcriptional activity, leading to cold-enhanced RH elongation through the consequent activation of the transcription factor gene RHD6-like RSL4. Furthermore, we demonstrate that APOLO interacts with the transcription factor WRKY42 and modulates its binding to the RHD6 promoter. WRKY42 is required for the activation of RHD6 by low temperatures and WRKY42 deregulation impairs cold-induced RH expansion. Collectively, our results indicate that a novel ribonucleoprotein complex with APOLO and WRKY42 forms a regulatory hub to activate RHD6 by shaping its epigenetic environment and integrate signals governing RH growth and development.

The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution

The scarcity of accessible sites that are dynamic or cell type-specific in plants may be due in part to tissue heterogeneity in bulk studies. To assess the effects of tissue heterogeneity, we apply single-cell ATAC-seq to Arabidopsis thaliana roots and identify thousands of differentially accessible sites, sufficient to resolve all major cell types of the root. We find that the entirety of a cell's regulatory landscape and its transcriptome independently capture cell type identity. We leverage this shared information on cell identity to integrate accessibility and transcriptome data to characterize developmental progression, endoreduplication and cell division. We further use the combined data to characterize cell type-specific motif enrichments of transcription factor families and link the expression of family members to changing accessibility at specific loci, resolving direct and indirect effects that shape expression. Our approach provides an analytical framework to infer the gene regulatory networks that execute plant development.

The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution

The scarcity of accessible sites that are dynamic or cell type-specific in plants may be due in part to tissue heterogeneity in bulk studies. To assess the effects of tissue heterogeneity, we apply single-cell ATAC-seq to Arabidopsis thaliana roots and identify thousands of differentially accessible sites, sufficient to resolve all major cell types of the root. We find that the entirety of a cell's regulatory landscape and its transcriptome independently capture cell type identity. We leverage this shared information on cell identity to integrate accessibility and transcriptome data to characterize developmental progression, endoreduplication and cell division. We further use the combined data to characterize cell type-specific motif enrichments of transcription factor families and link the expression of family members to changing accessibility at specific loci, resolving direct and indirect effects that shape expression. Our approach provides an analytical framework to infer the gene regulatory networks that execute plant development.

Anno genominis XX: 20 years of Arabidopsis genomics

Twenty years ago, the Arabidopsis thaliana genome sequence was published. This was an important moment as it was the first sequenced plant genome and explicitly brought plant science into the genomics era. At the time, this was not only an outstanding technological achievement, but it was characterized by a superb global collaboration. The Arabidopsis genome was the seed for plant genomic research. Here, we review the development of numerous resources based on the genome that have enabled discoveries across plant species, which has enhanced our understanding of how plants function and interact with their environments.

CRISPR screens in plants: approaches, guidelines, and future prospects

Clustered regularly interspaced short palindromic repeat (CRISPR)-associated systems have revolutionized genome engineering by facilitating a wide range of targeted DNA perturbations. These systems have resulted in the development of powerful new screens to test gene functions at the genomic scale. While there is tremendous potential to map and interrogate gene regulatory networks at unprecedented speed and scale using CRISPR screens, their implementation in plants remains in its infancy. Here we discuss the general concepts, tools, and workflows for establishing CRISPR screens in plants and analyze the handful of recent reports describing the use of this strategy to generate mutant knockout collections or to diversify DNA sequences. In addition, we provide insight into how to design CRISPR knockout screens in plants given the current challenges and limitations and examine multiple design options. Finally, we discuss the unique multiplexing capabilities of CRISPR screens to investigate redundant gene functions in highly duplicated plant genomes. Combinatorial mutant screens have the potential to routinely generate higher-order mutant collections and facilitate the characterization of gene networks. By integrating this approach with the numerous genomic profiles that have been generated over the past two decades, the implementation of CRISPR screens offers new opportunities to analyze plant genomes at deeper resolution and will lead to great advances in functional and synthetic biology.

Single-cell RNA sequencing of batch Chlamydomonas cultures reveals heterogeneity in their diurnal cycle phase

The photosynthetic unicellular alga Chlamydomonas (Chlamydomonas reinhardtii) is a versatile reference for algal biology because of its ease of culture in the laboratory. Genomic and systems biology approaches have previously described transcriptome responses to environmental changes using bulk data, thus representing the average behavior from pools of cells. Here, we apply single-cell RNA sequencing (scRNA-seq) to probe the heterogeneity of Chlamydomonas cell populations under three environments and in two genotypes differing by the presence of a cell wall. First, we determined that RNA can be extracted from single algal cells with or without a cell wall, offering the possibility to sample natural algal communities. Second, scRNA-seq successfully separated single cells into nonoverlapping cell clusters according to their growth conditions. Cells exposed to iron or nitrogen deficiency were easily distinguished despite a shared tendency to arrest photosynthesis and cell division to economize resources. Notably, these groups of cells not only recapitulated known patterns observed with bulk RNA-seq but also revealed their inherent heterogeneity. A substantial source of variation between cells originated from their endogenous diurnal phase, although cultures were grown in constant light. We exploited this result to show that circadian iron responses may be conserved from algae to land plants. We document experimentally that bulk RNA-seq data represent an average of typically hidden heterogeneity in the population.