Spatiotemporal Developmental Trajectories in the Arabidopsis Root Revealed Using High-Throughput Single-Cell RNA Sequencing

Amino acids biosynthesis in root hair development: a mini-review

Metabolic factors are essential for developmental biology of an organism. In plants, roots fulfill important functions, in part due to the development of specific epidermal cells, called hair cells that form root hairs (RHs) responsible for water and mineral uptake. RH development consists in (a) patterning processes involved in formation of hair and non-hair cells developed from trichoblasts and atrichoblasts; (b) RH initiation; and (c) apical (tip) growth of the RH. Here we review how these processes depend on pools of different amino acids and what is known about RH phenotypes of mutants disrupted in amino acid biosynthesis. This analysis shows that some amino acids, particularly aromatic ones, are required for RH apical (tip) growth, and that not much is known about the role of amino acids at earlier stages of RH formation. We also address the role of amino acids in rhizosphere, inhibitory and stimulating effects of amino acids on RH growth, amino acids as N source in plant nutrition, and amino acid transporters and their expression in the RHs. Amino acids form conjugates with auxin, a hormone essential for RH growth, and respective genes are overviewed. Finally, we outline missing links and envision some perspectives in the field.

Streamlined spatial and environmental expression signatures characterize the minimalist duckweed Wolffia australiana

Single-cell genomics permits a new resolution in the examination of molecular and cellular dynamics, allowing global, parallel assessments of cell types and cellular behaviors through development and in response to environmental circumstances, such as interaction with water and the light-dark cycle of the Earth. Here, we leverage the smallest, and possibly most structurally reduced, plant, the semiaquatic Wolffia australiana, to understand dynamics of cell expression in these contexts at the whole-plant level. We examined single-cell-resolution RNA-sequencing data and found Wolffia cells divide into four principal clusters representing the above- and below-water-situated parenchyma and epidermis. Although these tissues share transcriptomic similarity with model plants, they display distinct adaptations that Wolffia has made for the aquatic environment. Within this broad classification, discrete subspecializations are evident, with select cells showing unique transcriptomic signatures associated with developmental maturation and specialized physiologies. Assessing this simplified biological system temporally at two key time-of-day (TOD) transitions, we identify additional TOD-responsive genes previously overlooked in whole-plant transcriptomic approaches and demonstrate that the core circadian clock machinery and its downstream responses can vary in cell-specific manners, even in this simplified system. Distinctions between cell types and their responses to submergence and/or TOD are driven by expression changes of unexpectedly few genes, characterizing Wolffia as a highly streamlined organism with the majority of genes dedicated to fundamental cellular processes. Wolffia provides a unique opportunity to apply reductionist biology to elucidate signaling functions at the organismal level, for which this work provides a powerful resource.

Single-cell transcriptome atlases of soybean root and mature nodule reveal new regulatory programs that control the nodulation process

The soybean root system is complex. In addition to being composed of various cell types, the soybean root system includes the primary root, the lateral roots, and the nodule, an organ in which mutualistic symbiosis with N-fixing rhizobia occurs. A mature soybean root nodule is characterized by a central infection zone where atmospheric nitrogen is fixed and assimilated by the symbiont, resulting from the close cooperation between the plant cell and the bacteria. To date, the transcriptome of individual cells isolated from developing soybean nodules has been established, but the transcriptomic signatures of cells from the mature soybean nodule have not yet been characterized. Using single-nucleus RNA-seq and Molecular Cartography technologies, we precisely characterized the transcriptomic signature of soybean root and mature nodule cell types and revealed the co-existence of different sub-populations of B. diazoefficiens-infected cells in the mature soybean nodule, including those actively involved in nitrogen fixation and those engaged in senescence. Mining of the single-cell-resolution nodule transcriptome atlas and the associated gene co-expression network confirmed the role of known nodulation-related genes and identified new genes that control the nodulation process. For instance, we functionally characterized the role of GmFWL3, a plasma membrane microdomain-associated protein that controls rhizobial infection. Our study reveals the unique cellular complexity of the mature soybean nodule and helps redefine the concept of cell types when considering the infection zone of the soybean nodule.

Single-cell network analysis reveals gene expression programs for Arabidopsis root development and metabolism

Single-cell RNA-sequencing datasets of Arabidopsis roots have been generated, but related comprehensive gene co-expression network analyses are lacking. We conducted a single-cell gene co-expression network analysis with publicly available scRNA-seq datasets of Arabidopsis roots using a SingleCellGGM algorithm. The analysis identified 149 gene co-expression modules, which we considered to be gene expression programs (GEPs). By examining their spatiotemporal expression, we identified GEPs specifically expressed in major root cell types along their developmental trajectories. These GEPs define gene programs regulating root cell development at different stages and are enriched with relevant developmental regulators. As examples, a GEP specific for the quiescent center (QC) contains 20 genes regulating QC and stem cell niche homeostasis, and four GEPs are expressed in sieve elements (SEs) from early to late developmental stages, with the early-stage GEP containing 17 known SE developmental regulators. We also identified GEPs for metabolic pathways with cell-type-specific expression, suggesting the existence of cell-type-specific metabolism in roots. Using the GEPs, we discovered and verified a columella-specific gene, NRL27, as a regulator of the auxin-related root gravitropism response. Our analysis thus systematically reveals GEPs that regulate Arabidopsis root development and metabolism and provides ample resources for root biology studies.

A single-cell transcriptome atlas reveals the trajectory of early cell fate transition during callus induction in Arabidopsis

The acquisition of pluripotent callus from somatic cells plays an important role in plant development studies and crop genetic improvement. This developmental process incorporates a series of cell fate transitions and reprogramming. However, our understanding of cell heterogeneity and mechanisms of cell fate transition during callus induction remains quite limited. Here, we report a time-series single-cell transcriptome experiment on Arabidopsis root explants that were induced in callus induction medium for 0, 1, and 4 days, and the construction of a detailed single-cell transcriptional atlas of the callus induction process. We identify the cell types responsible for initiating the early callus: lateral root primordium-initiating (LRPI)-like cells and quiescent center (QC)-like cells. LRPI-like cells are derived from xylem pole pericycle cells and are similar to lateral root primordia. We delineate the developmental trajectory of the dedifferentiation of LRPI-like cells into QC-like cells. QC-like cells are undifferentiated pluripotent acquired cells that appear in the early stages of callus formation and play a critical role in later callus development and organ regeneration. We also identify the transcription factors that regulate QC-like cells and the gene expression signatures that are related to cell fate decisions. Overall, our cell-lineage transcriptome atlas for callus induction provides a distinct perspective on cell fate transitions during callus formation, significantly improving our understanding of callus formation.

Single-nucleus RNA and ATAC sequencing analyses provide molecular insights into early pod development of peanut fruit

Peanut (Arachis hypogaea L.) is an important leguminous oil and economic crop that produces flowers aboveground and fruits belowground. Subterranean fruit-pod development, which significantly affects peanut production, involves complex molecular mechanisms that likely require the coordinated regulation of multiple genes in different tissues. To investigate the molecular mechanisms that underlie peanut fruit-pod development, we characterized the anatomical features of early fruit-pod development and integrated single-nucleus RNA-sequencing (snRNA-seq) and single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) data at the single-cell level. We identified distinct cell types, such as meristem, embryo, vascular tissue, cuticular layer, and stele cells within the shell wall. These specific cell types were used to examine potential molecular changes unique to each cell type during pivotal stages of fruit-pod development. snRNA-seq analyses of differentially expressed genes revealed cell-type-specific insights that were not previously obtainable from transcriptome analyses of bulk RNA. For instance, we identified MADS-box genes that contributes to the formation of parenchyma cells and gravity-related genes that are present in the vascular cells, indicating an essential role for the vascular cells in peg gravitropism. Overall, our single-nucleus analysis provides comprehensive and novel information on specific cell types, gene expression, and chromatin accessibility during the early stages of fruit-pod development. This information will enhance our understanding of the mechanisms that underlie fruit-pod development in peanut and contribute to efforts aimed at improving peanut production.

Artificial intelligence in plant breeding

Harnessing cutting-edge technologies to enhance crop productivity is a pivotal goal in modern plant breeding. Artificial intelligence (AI) is renowned for its prowess in big data analysis and pattern recognition, and is revolutionizing numerous scientific domains including plant breeding. We explore the wider potential of AI tools in various facets of breeding, including data collection, unlocking genetic diversity within genebanks, and bridging the genotype-phenotype gap to facilitate crop breeding. This will enable the development of crop cultivars tailored to the projected future environments. Moreover, AI tools also hold promise for refining crop traits by improving the precision of gene-editing systems and predicting the potential effects of gene variants on plant phenotypes. Leveraging AI-enabled precision breeding can augment the efficiency of breeding programs and holds promise for optimizing cropping systems at the grassroots level. This entails identifying optimal inter-cropping and crop-rotation models to enhance agricultural sustainability and productivity in the field.

Unraveling plant-microbe symbioses using single-cell and spatial transcriptomics

Plant-microbe symbioses require intense interaction and genetic coordination to successfully establish in specific cell types of the host and symbiont. Traditional RNA-seq methodologies lack the cellular resolution to fully capture these complexities, but single-cell and spatial transcriptomics (ST) are now allowing scientists to probe symbiotic interactions at an unprecedented level of detail. Here, we discuss the advantages that novel spatial and single-cell transcriptomic technologies provide in studying plant-microbe endosymbioses and highlight key recent studies. Finally, we consider the remaining limitations of applying these approaches to symbiosis research, which are mainly related to the simultaneous capture of both plant and microbial transcripts within the same cells.

Are We There Yet? Assessing the Readiness of Single-Cell Proteomics to Answer Biological Hypotheses

Single-cell analysis is an active area of research in many fields of biology. Measurements at single-cell resolution allow researchers to study diverse populations without losing biologically meaningful information to sample averages. Many technologies have been used to study single cells, including mass spectrometry-based single-cell proteomics (SCP). SCP has seen a lot of growth over the past couple of years through improvements in data acquisition and analysis, leading to greater proteomic depth. Because method development has been the main focus in SCP, biological applications have been sprinkled in only as proof-of-concept. However, SCP methods now provide significant coverage of the proteome and have been implemented in many laboratories. Thus, a primary question to address in our community is whether the current state of technology is ready for widespread adoption for biological inquiry. In this Perspective, we examine the potential for SCP in three thematic areas of biological investigation: cell annotation, developmental trajectories, and spatial mapping. We identify that the primary limitation of SCP is sample throughput. As proteome depth has been the primary target for method development to date, we advocate for a change in focus to facilitate measuring tens of thousands of single-cell proteomes to enable biological applications beyond proof-of-concept.

Multiplexed in situ hybridization reveals distinct lineage identities for major and minor vein initiation during maize leaf development

Leaves of flowering plants are characterized by diverse venation patterns. Patterning begins with the selection of vein-forming procambial initial cells from within the ground meristem of a developing leaf, a process which is considered to be auxin-dependent, and continues until veins are anatomically differentiated with functional xylem and phloem. At present, the mechanisms responsible for leaf venation patterning are primarily characterized in the model eudicot Arabidopsis thaliana which displays a reticulate venation network. However, evidence suggests that vein development may proceed via a different mechanism in monocot leaves where venation patterning is parallel. Here, we employed Molecular Cartography, a multiplexed in situ hybridization technique, to analyze the spatiotemporal localization of a subset of auxin-related genes and candidate regulators of vein patterning in maize leaves. We show how different combinations of auxin influx and efflux transporters are recruited during leaf and vein specification and how major and minor vein ranks develop with distinct identities. The localization of the procambial marker PIN1a and the spatial arrangement of procambial initial cells that give rise to major and minor vein ranks further suggests that vein spacing is prepatterned across the medio-lateral leaf axis prior to accumulation of the PIN1a auxin transporter. In contrast, patterning in the adaxial-abaxial axis occurs progressively, with markers of xylem and phloem gradually becoming polarized as differentiation proceeds. Collectively, our data suggest that both lineage- and position-based mechanisms may underpin vein patterning in maize leaves.

Opportunities and Challenges in Advancing Plant Research with Single-cell Omics

Plants possess diverse cell types and intricate regulatory mechanisms to adapt to the ever-changing environment of nature. Various strategies have been employed to study cell types and their developmental progressions, including single-cell sequencing methods which provide high-dimensional catalogs to address biological concerns. In recent years, single-cell sequencing technologies in transcriptomics, epigenomics, proteomics, metabolomics, and spatial transcriptomics have been increasingly used in plant science to reveal intricate biological relationships at the single-cell level. However, the application of single-cell technologies to plants is more limited due to the challenges posed by cell structure. This review outlines the advancements in single-cell omics technologies, their implications in plant systems, future research applications, and the challenges of single-cell omics in plant systems.

Single-cell transcriptome landscape elucidates the cellular and developmental responses to tomato chlorosis virus infection in tomato leaf

Plant viral diseases compromise the growth and yield of the crop globally, and they tend to be more serious under extreme temperatures and drought climate changes. Currently, regulatory dynamics during plant development and in response to virus infection at the plant cell level remain largely unknown. In this study, single-cell RNA sequencing on 23 226 individual cells from healthy and tomato chlorosis virus-infected leaves was established. The specific expression and epigenetic landscape of each cell type during the viral infection stage were depicted. Notably, the mesophyll cells showed a rapid function transition in virus-infected leaves, which is consistent with the pathological changes such as thinner leaves and decreased chloroplast lamella in virus-infected samples. Interestingly, the F-box protein SKIP2 was identified to play a pivotal role in chlorophyll maintenance during virus infection in tomato plants. Knockout of the SlSKIP2 showed a greener leaf state before and after virus infection. Moreover, we further demonstrated that SlSKIP2 was located in the cytomembrane and nucleus and directly regulated by ERF4. In conclusion, with detailed insights into the plant responses to viral infections at the cellular level, our study provides a genetic framework and gene reference in plant-virus interaction and breeding in the future research.

Coexpression enhances cross-species integration of single-cell RNA sequencing across diverse plant species

Single-cell RNA sequencing is increasingly used to investigate cross-species differences driven by gene expression and cell-type composition in plants. However, the frequent expansion of plant gene families due to whole-genome duplications makes identification of one-to-one orthologues difficult, complicating integration. Here we demonstrate that coexpression can be used to trim many-to-many orthology families down to identify one-to-one gene pairs with proxy expression profiles, improving the performance of traditional integration methods and reducing barriers to integration across a diverse array of plant species.

Global dynamics and cytokinin participation of salt gland development trajectory in recretohalophyte Limonium bicolor

Salt gland is an epidermal Na+ secretory structure that enhances salt resistance in the recretohalophyte sea lavender (Limonium bicolor). To elucidate the salt gland development trajectory and related molecular mechanisms, we performed single-cell RNA sequencing of L. bicolor protoplasts from young leaves at salt gland initiation and differentiation stages. Dimensionality reduction analyses defined 19 transcriptionally distinct cell clusters, which were assigned into 4 broad populations-promeristem, epidermis, mesophyll, and vascular tissue-verified by in situ hybridization. Cytokinin was further proposed to participate in salt gland development by the expression patterns of related genes and cytological evidence. By comparison analyses of Single-cell RNA sequencing with exogenous application of 6-benzylaminopurine, we delineated 5 salt gland development-associated subclusters and defined salt gland-specific differentiation trajectories from Subclusters 8, 4, and 6 to Subcluster 3 and 1. Additionally, we validated the participation of TRIPTYCHON and the interacting protein Lb7G34824 in salt gland development, which regulated the expression of cytokinin metabolism and signaling-related genes such as GLABROUS INFLORESCENCE STEMS 2 to maintain cytokinin homeostasis during salt gland development. Our results generated a gene expression map of young leaves at single-cell resolution for the comprehensive investigation of salt gland determinants and cytokinin participation that helps elucidate cell fate determination during epidermis formation and evolution in recretohalophytes.

Advances in Single-Cell Transcriptome Sequencing and Spatial Transcriptome Sequencing in Plants

"Omics" typically involves exploration of the structure and function of the entire composition of a biological system at a specific level using high-throughput analytical methods to probe and analyze large amounts of data, including genomics, transcriptomics, proteomics, and metabolomics, among other types. Genomics characterizes and quantifies all genes of an organism collectively, studying their interrelationships and their impacts on the organism. However, conventional transcriptomic sequencing techniques target population cells, and their results only reflect the average expression levels of genes in population cells, as they are unable to reveal the gene expression heterogeneity and spatial heterogeneity among individual cells, thus masking the expression specificity between different cells. Single-cell transcriptomic sequencing and spatial transcriptomic sequencing techniques analyze the transcriptome of individual cells in plant or animal tissues, enabling the understanding of each cell's metabolites and expressed genes. Consequently, statistical analysis of the corresponding tissues can be performed, with the purpose of achieving cell classification, evolutionary growth, and physiological and pathological analyses. This article provides an overview of the research progress in plant single-cell and spatial transcriptomics, as well as their applications and challenges in plants. Furthermore, prospects for the development of single-cell and spatial transcriptomics are proposed.

Spatiotemporal transcriptomic atlas of rhizome formation in Oryza longistaminata

Rhizomes are modified stems that grow underground and produce new individuals genetically identical to the mother plant. Recently, a breakthrough has been made in efforts to convert annual grains into perennial ones by utilizing wild rhizomatous species as donors, yet the developmental biology of this organ is rarely studied. Oryza longistaminata, a wild rice species featuring strong rhizomes, provides a valuable model for exploration of rhizome development. Here, we first assembled a double-haplotype genome of O. longistaminata, which displays a 48-fold improvement in contiguity compared to the previously published assembly. Furthermore, spatiotemporal transcriptomics was performed to obtain the expression profiles of different tissues in O. longistaminata rhizomes and tillers. Two spatially reciprocal cell clusters, the vascular bundle 2 cluster and the parenchyma 2 cluster, were determined to be the primary distinctions between the rhizomes and tillers. We also captured meristem initiation cells in the sunken area of parenchyma located at the base of internodes, which is the starting point for rhizome initiation. Trajectory analysis further indicated that the rhizome is regenerated through de novo generation. Collectively, these analyses revealed a spatiotemporal transcriptional transition underlying the rhizome initiation, providing a valuable resource for future perennial crop breeding.

Multiome in the Same Cell Reveals the Impact of Osmotic Stress on Arabidopsis Root Tip Development at Single-Cell Level

Cell-specific transcriptional regulatory networks (TRNs) play vital roles in plant development and response to environmental stresses. However, traditional single-cell mono-omics techniques are unable to directly capture the relationships and dynamics between different layers of molecular information within the same cells. While advanced algorithm facilitates merging scRNA-seq and scATAC-seq datasets, accurate data integration remains a challenge, particularly when investigating cell-type-specific TRNs. By examining gene expression and chromatin accessibility simultaneously in 16,670 Arabidopsis root tip nuclei, the TRNs are reconstructed that govern root tip development under osmotic stress. In contrast to commonly used computational integration at cell-type level, 12,968 peak-to-gene linkage is captured at the bona fide single-cell level and construct TRNs at an unprecedented resolution. Furthermore, the unprecedented datasets allow to more accurately reconstruct the coordinated changes of gene expression and chromatin states during cellular state transition. During root tip development, chromatin accessibility of initial cells precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for subsequent differentiation steps. Pseudo-time trajectory analysis reveal that osmotic stress can shift the functional differentiation of trichoblast. Candidate stress-related gene-linked cis-regulatory elements (gl-cCREs) as well as potential target genes are also identified, and uncovered large cellular heterogeneity under osmotic stress.

Insights into multilevel spatial regulation within the root stem cell niche

All differentiated root cells derive from stem cells spatially organized within the stem cell niche (SCN), a microenvironment located within the root tip. Here, we compiled recent advances in the understanding of how the SCN drives the establishment and maintenance of cell types. The quiescent center (QC) is widely recognized as the primary driver of cell fate determination, but it is recently considered a convergence center of multiple signals. Cell identity of the cortex endodermis initials is mainly driven by the regulatory feedback loops between transcription factors (TFs), acting as mobile signals between neighboring cells, including the QC. As exemplified in the vascular initials, the precise spatial expression of these regulatory TFs is connected with a dynamic hormonal interplay. Thus, stem cell maintenance and cell differentiation are regulated by a plethora of signals forming a complex, multilevel regulatory network. Integrating the transcriptional and post-translational regulations, protein-protein interactions, and mobile signals into models will be fundamental for the comprehensive understanding of SCN maintenance and differentiation.

Comprehensive integration of single-cell transcriptomic data illuminates the regulatory network architecture of plant cell fate specification

Plant morphogenesis relies on precise gene expression programs at the proper time and position which is orchestrated by transcription factors (TFs) in intricate regulatory networks in a cell-type specific manner. Here we introduced a comprehensive single-cell transcriptomic atlas of Arabidopsis seedlings. This atlas is the result of meticulous integration of 63 previously published scRNA-seq datasets, addressing batch effects and conserving biological variance. This integration spans a broad spectrum of tissues, including both below- and above-ground parts. Utilizing a rigorous approach for cell type annotation, we identified 47 distinct cell types or states, largely expanding our current view of plant cell compositions. We systematically constructed cell-type specific gene regulatory networks and uncovered key regulators that act in a coordinated manner to control cell-type specific gene expression. Taken together, our study not only offers extensive plant cell atlas exploration that serves as a valuable resource, but also provides molecular insights into gene-regulatory programs that varies from different cell types.

Nitrogen sensing and regulatory networks: it's about time and space

A plant's response to external and internal nitrogen signals/status relies on sensing and signaling mechanisms that operate across spatial and temporal dimensions. From a comprehensive systems biology perspective, this involves integrating nitrogen responses in different cell types and over long distances to ensure organ coordination in real time and yield practical applications. In this prospective review, we focus on novel aspects of nitrogen (N) sensing/signaling uncovered using temporal and spatial systems biology approaches, largely in the model Arabidopsis. The temporal aspects span: transcriptional responses to N-dose mediated by Michaelis-Menten kinetics, the role of the master NLP7 transcription factor as a nitrate sensor, its nitrate-dependent TF nuclear retention, its "hit-and-run" mode of target gene regulation, and temporal transcriptional cascade identified by "network walking." Spatial aspects of N-sensing/signaling have been uncovered in cell type-specific studies in roots and in root-to-shoot communication. We explore new approaches using single-cell sequencing data, trajectory inference, and pseudotime analysis as well as machine learning and artificial intelligence approaches. Finally, unveiling the mechanisms underlying the spatial dynamics of nitrogen sensing/signaling networks across species from model to crop could pave the way for translational studies to improve nitrogen-use efficiency in crops. Such outcomes could potentially reduce the detrimental effects of excessive fertilizer usage on groundwater pollution and greenhouse gas emissions.

Sugar transporters spatially organize microbiota colonization along the longitudinal root axis of Arabidopsis

Plant roots are functionally heterogeneous in cellular architecture, transcriptome profile, metabolic state, and microbial immunity. We hypothesized that axial differentiation may also impact spatial colonization by root microbiota along the root axis. We developed two growth systems, ArtSoil and CD-Rhizotron, to grow and then dissect Arabidopsis thaliana roots into three segments. We demonstrate that distinct endospheric and rhizosphere bacterial communities colonize the segments, supporting the hypothesis of microbiota differentiation along the axis. Root metabolite profiling of each segment reveals differential metabolite enrichment and specificity. Bioinformatic analyses and GUS histochemistry indicate microbe-induced accumulation of SWEET2, 4, and 12 sugar uniporters. Profiling of root segments from sweet mutants shows altered spatial metabolic profiles and reorganization of endospheric root microbiota. This work reveals the interdependency between root metabolites and microbial colonization and the contribution of SWEETs to spatial diversity and stability of microbial ecosystem.

Merit of integrating in situ transcriptomics and anatomical information for cell annotation and lineage construction in single-cell analyses of Populus

Cell type annotation and lineage construction are two of the most critical tasks conducted in the analyses of single-cell RNA sequencing (scRNA-seq). Four recent scRNA-seq studies of differentiating xylem propose four models on differentiating xylem development in Populus. The differences are mostly caused by the use of different strategies for cell type annotation and subsequent lineage interpretation. Here, we emphasize the necessity of using in situ transcriptomes and anatomical information to construct the most plausible xylem development model.

Lost in space: what single-cell RNA sequencing cannot tell you

Plant scientists are rapidly integrating single-cell RNA sequencing (scRNA-seq) into their workflows. Maximizing the potential of scRNA-seq requires a proper understanding of the spatiotemporal context of cells. However, positional information is inherently lost during scRNA-seq, limiting its potential to characterize complex biological systems. In this review we highlight how current single-cell analysis pipelines cannot completely recover spatial information, which confounds biological interpretation. Various strategies exist to identify the location of RNA, from classical RNA in situ hybridization to spatial transcriptomics. Herein we discuss the possibility of utilizing this spatial information to supervise single-cell analyses. An integrative approach will maximize the potential of each technology, and lead to insights which go beyond the capability of each individual technology.

Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis

The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

Studying plant vascular development using single-cell approaches

Vascular cells form a highly complex and heterogeneous tissue. Its composition, function, shape, and arrangement vary with the developmental stage and between organs and species. Understanding the transcriptional regulation underpinning this complexity thus requires a high-resolution technique that is capable of capturing rapid events during vascular cell formation. Single-cell and single-nucleus RNA sequencing (sc/snRNA-seq) approaches provide powerful tools to extract transcriptional information from these lowly abundant and dynamically changing cell types, which allows the reconstruction of developmental trajectories. Here, we summarize and reflect on recent studies using single-cell transcriptomics to study vascular cell types and discuss current and future implementations of sc/snRNA-seq approaches in the field of vascular development.

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

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

Unraveling the molecular mechanisms governing axillary meristem initiation in plants

MAIN CONCLUSION

Axillary meristems (AMs) located in the leaf axils determine the number of shoots or tillers eventually formed, thus contributing significantly to the plant architecture and crop yields. The study of AM initiation is unavoidable and beneficial for crop productivity. Shoot branching is an undoubted determinant of plant architecture and thus greatly impacts crop yield due to the panicle-bearing traits of tillers. The emergence of the AM is essential for the incipient bud formation, and then the bud is dormant or outgrowth immediately to form a branch or tiller. While numerous reviews have focused on plant branching and tillering development networks, fewer specifically address AM initiation and its regulatory mechanisms. This review synthesizes the significant advancements in the genetic and hormonal factors governing AM initiation, with a primary focus on studies conducted in Arabidopsis (Arabidopsis thaliana L.) and rice (Oryza sativa L.). In particular, by elaborating on critical genes like LATERAL SUPPRESSOR (LAS), which specifically regulates AM initiation and the networks in which they are involved, we attempt to unify the cascades through which they are positioned. We concentrate on clarifying the precise mutual regulation between shoot apical meristem (SAM) and AM-related factors. Additionally, we examine challenges in elucidating AM formation mechanisms alongside opportunities provided by emerging omics approaches to identify AM-specific genes. By expanding our comprehension of the genetic and hormonal regulation of AM development, we can develop strategies to optimize crop production and address global food challenges effectively.

Single-cell transcriptomics is revolutionizing the improvement of plant biotechnology research: recent advances and future opportunities

Single-cell approaches are a promising way to obtain high-resolution transcriptomics data and have the potential to revolutionize the study of plant growth and development. Recent years have seen the advent of unprecedented technological advances in the field of plant biology to study the transcriptional information of individual cells by single-cell RNA sequencing (scRNA-seq). This review focuses on the modern advancements of single-cell transcriptomics in plants over the past few years. In addition, it also offers a new insight of how these emerging methods will expedite advance research in plant biotechnology in the near future. Lastly, the various technological hurdles and inherent limitations of single-cell technology that need to be conquered to develop such outstanding possible knowledge gain is critically analyzed and discussed.

Rapid alkalinization factor 22 has a structural and signalling role in root hair cell wall assembly

Pressurized cells with strong walls make up the hydrostatic skeleton of plants. Assembly and expansion of such stressed walls depend on a family of secreted RAPID ALKALINIZATION FACTOR (RALF) peptides, which bind both a membrane receptor complex and wall-localized LEUCINE-RICH REPEAT EXTENSIN (LRXs) in a mutually exclusive way. Here we show that, in root hairs, the RALF22 peptide has a dual structural and signalling role in cell expansion. Together with LRX1, it directs the compaction of charged pectin polymers at the root hair tip into periodic circumferential rings. Free RALF22 induces the formation of a complex with LORELEI-LIKE-GPI-ANCHORED PROTEIN 1 and FERONIA, triggering adaptive cellular responses. These findings show how a peptide simultaneously functions as a structural component organizing cell wall architecture and as a feedback signalling molecule that regulates this process depending on its interaction partners. This mechanism may also underlie wall assembly and expansion in other plant cell types.

Advances in the Application of Single-Cell Transcriptomics in Plant Systems and Synthetic Biology

Plants are complex systems hierarchically organized and composed of various cell types. To understand the molecular underpinnings of complex plant systems, single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for revealing high resolution of gene expression patterns at the cellular level and investigating the cell-type heterogeneity. Furthermore, scRNA-seq analysis of plant biosystems has great potential for generating new knowledge to inform plant biosystems design and synthetic biology, which aims to modify plants genetically/epigenetically through genome editing, engineering, or re-writing based on rational design for increasing crop yield and quality, promoting the bioeconomy and enhancing environmental sustainability. In particular, data from scRNA-seq studies can be utilized to facilitate the development of high-precision Build-Design-Test-Learn capabilities for maximizing the targeted performance of engineered plant biosystems while minimizing unintended side effects. To date, scRNA-seq has been demonstrated in a limited number of plant species, including model plants (e.g., Arabidopsis thaliana), agricultural crops (e.g., Oryza sativa), and bioenergy crops (e.g., Populus spp.). It is expected that future technical advancements will reduce the cost of scRNA-seq and consequently accelerate the application of this emerging technology in plants. In this review, we summarize current technical advancements in plant scRNA-seq, including sample preparation, sequencing, and data analysis, to provide guidance on how to choose the appropriate scRNA-seq methods for different types of plant samples. We then highlight various applications of scRNA-seq in both plant systems biology and plant synthetic biology research. Finally, we discuss the challenges and opportunities for the application of scRNA-seq in plants.

The single-cell transcriptome program of nodule development cellular lineages in Medicago truncatula

Legumes establish a symbiotic relationship with nitrogen-fixing rhizobia by developing nodules. Nodules are modified lateral roots that undergo changes in their cellular development in response to bacteria, but the transcriptional reprogramming that occurs in these root cells remains largely uncharacterized. Here, we describe the cell-type-specific transcriptome response of Medicago truncatula roots to rhizobia during early nodule development in the wild-type genotype Jemalong A17, complemented with a hypernodulating mutant (sunn-4) to expand the cell population responding to infection and subsequent biological inferences. The analysis identifies epidermal root hair and stele sub-cell types associated with a symbiotic response to infection and regulation of nodule proliferation. Trajectory inference shows cortex-derived cell lineages differentiating to form the nodule primordia and, posteriorly, its meristem, while modulating the regulation of phytohormone-related genes. Gene regulatory analysis of the cell transcriptomes identifies new regulators of nodulation, including STYLISH 4, for which the function is validated.

Plant biotechnology research with single-cell transcriptome: recent advancements and prospects

KEY MESSAGE

Single-cell transcriptomic techniques have emerged as powerful tools in plant biology, offering high-resolution insights into gene expression at the individual cell level. This review highlights the rapid expansion of single-cell technologies in plants, their potential in understanding plant development, and their role in advancing plant biotechnology research. Single-cell techniques have emerged as powerful tools to enhance our understanding of biological systems, providing high-resolution transcriptomic analysis at the single-cell level. In plant biology, the adoption of single-cell transcriptomics has seen rapid expansion of available technologies and applications. This review article focuses on the latest advancements in the field of single-cell transcriptomic in plants and discusses the potential role of these approaches in plant development and expediting plant biotechnology research in the near future. Furthermore, inherent challenges and limitations of single-cell technology are critically examined to overcome them and enhance our knowledge and understanding.

Evolution of plant cell-type-specific cis -regulatory elements

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species ( Zea mays, Sorghum bicolor, Panicum miliaceum , and Urochloa fusca ). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1 epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Light controls mesophyll-specific post-transcriptional splicing of photoregulatory genes by AtPRMT5

Light plays a central role in plant growth and development, providing an energy source and governing various aspects of plant morphology. Previous study showed that many polyadenylated full-length RNA molecules within the nucleus contain unspliced introns (post-transcriptionally spliced introns, PTS introns), which may play a role in rapidly responding to changes in environmental signals. However, the mechanism underlying post-transcriptional regulation during initial light exposure of young, etiolated seedlings remains elusive. In this study, we used FLEP-seq2, a Nanopore-based sequencing technique, to analyze nuclear RNAs in Arabidopsis (Arabidopsis thaliana) seedlings under different light conditions and found numerous light-responsive PTS introns. We also used single-nucleus RNA sequencing (snRNA-seq) to profile transcripts in single nucleus and investigate the distribution of light-responsive PTS introns across distinct cell types. We established that light-induced PTS introns are predominant in mesophyll cells during seedling de-etiolation following exposure of etiolated seedlings to light. We further demonstrated the involvement of the splicing-related factor A. thaliana PROTEIN ARGININE METHYLTRANSFERASE 5 (AtPRMT5), working in concert with the E3 ubiquitin ligase CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1), a critical repressor of light signaling pathways. We showed that these two proteins orchestrate light-induced PTS events in mesophyll cells and facilitate chloroplast development, photosynthesis, and morphogenesis in response to ever-changing light conditions. These findings provide crucial insights into the intricate mechanisms underlying plant acclimation to light at the cell-type level.

Single-cell RNA-seq reveals a link of ovule abortion and sugar transport in Camellia oleifera

Camellia oleifera is the most important woody oil crop in China. Seed number per fruit is an important yield trait in C. oleifera. Ovule abortion is generally observed in C. oleifera and significantly decreases the seed number per fruit. However, the mechanisms of ovule abortion remain poorly understood at present. Single-cell RNA sequencing (scRNA-seq) was performed using mature ovaries of two C. oleifera varieties with different ovule abortion rates (OARs). In total, 20,526 high-quality cells were obtained, and 18 putative cell clusters were identified. Six cell types including female gametophyte, protoxylem, protophloem, procambium, epidermis, and parenchyma cells were identified from three main tissue types of ovule, placenta, and pericarp inner layer. A comparative analysis on scRNA-seq data between high- and low-OAR varieties demonstrated that the overall expression of CoSWEET and CoCWINV in procambium cells, and CoSTP in the integument was significantly upregulated in the low-OAR variety. Both the infertile ovule before pollination and the abortion ovule producing after compatible pollination might be attributed to selective abortion caused by low sugar levels in the apoplast around procambium cells and a low capability of hexose uptake in the integument. Here, the first single-cell transcriptional landscape is reported in woody crop ovaries. Our investigation demonstrates that ovule abortion may be related to sugar transport in placenta and ovules and sheds light on further deciphering the mechanism of regulating sugar transport and the improvement of seed yield in C. oleifera.

Introducing single cell stereo-sequencing technology to transform the plant transcriptome landscape

Single cell RNA-sequencing (scRNA-seq) advancements have helped detect transcriptional heterogeneities in biological samples. However, scRNA-seq cannot currently provide high-resolution spatial transcriptome information or identify subcellular organs in biological samples. These limitations have led to the development of spatially enhanced-resolution omics-sequencing (Stereo-seq), which combines spatial information with single cell transcriptomics to address the challenges of scRNA-seq alone. In this review, we discuss the advantages of Stereo-seq technology. We anticipate that the application of such an integrated approach in plant research will advance our understanding of biological process in the plant transcriptomics era. We conclude with an outlook of how such integration will enhance crop improvement.

Single-cell transcriptomic analysis of pea shoot development and cell-type-specific responses to boron deficiency

Understanding how nutrient stress impacts plant growth is fundamentally important to the development of approaches to improve crop production under nutrient limitation. Here we applied single-cell RNA sequencing to shoot apices of Pisum sativum grown under boron (B) deficiency. We identified up to 15 cell clusters based on the clustering of gene expression profiles and verified cell identity with cell-type-specific marker gene expression. Different cell types responded differently to B deficiency. Specifically, the expression of photosynthetic genes in mesophyll cells (MCs) was down-regulated by B deficiency, consistent with impaired photosynthetic rate. Furthermore, the down-regulation of stomatal development genes in guard cells, including homologs of MUTE and TOO MANY MOUTHS, correlated with a decrease in stomatal density under B deficiency. We also constructed the developmental trajectory of the shoot apical meristem (SAM) cells and a transcription factor interaction network. The developmental progression of SAM to MC was characterized by up-regulation of genes encoding histones and chromatin assembly and remodeling proteins including homologs of FASCIATA1 (FAS1) and SWITCH DEFECTIVE/SUCROSE NON-FERMENTABLE (SWI/SNF) complex. However, B deficiency suppressed their expression, which helps to explain impaired SAM development under B deficiency. These results represent a major advance over bulk-tissue RNA-seq analysis in which cell-type-specific responses are lost and hence important physiological responses to B deficiency are missed. The reported findings reveal strategies by which plants adapt to B deficiency thus offering breeders a set of specific targets for genetic improvement. The reported approach and resources have potential applications well beyond P. sativum species and could be applied to various legumes to improve their adaptability to multiple nutrient or abiotic stresses.

A Single-Nucleus Resolution Atlas of Transcriptome and Chromatin Accessibility for Peanut (Arachis Hypogaea L.) Leaves

The peanut is an important worldwide cash-crop for edible oil and protein. However, the kinetic mechanisms that determine gene expression and chromatin accessibility during leaf development in peanut represented allotetraploid leguminous crops are poorly understood at single-cell resolution. Here, a single-nucleus atlas of peanut leaves is developed by simultaneously profiling the transcriptome and chromatin accessibility in the same individual-cell using fluorescence-activated sorted single-nuclei. In total, 5930 cells with 50 890 expressed genes are classified into 18 cell-clusters, and 5315 chromatin fragments are enriched with 26 083 target genes in the chromatin accessible landscape. The developmental trajectory analysis reveals the involvement of the ethylene-AP2 module in leaf cell differentiation, and cell-cycle analysis demonstrated that genome replication featured in distinct cell-types with circadian rhythms transcription factors (TFs). Furthermore, dual-omics illustrates that the fatty acid pathway modulates epidermal-guard cells differentiation and providescritical TFs interaction networks for understanding mesophyll development, and the cytokinin module (LHY/LOG) that regulates vascular growth. Additionally, an AT-hook protein AhAHL11 is identified that promotes leaf area expansion by modulating the auxin content increase. In summary, the simultaneous profiling of transcription and chromatin accessibility landscapes using snRNA/ATAC-seq provides novel biological insights into the dynamic processes of peanut leaf cell development at the cellular level.

Single-cell RNA sequencing reveals a hierarchical transcriptional regulatory network of terpenoid biosynthesis in cotton secretory glandular cells

Plants can synthesize a wide range of terpenoids in response to various environmental cues. However, the specific regulatory mechanisms governing terpenoid biosynthesis at the cellular level remain largely elusive. In this study, we employed single-cell RNA sequencing to comprehensively characterize the transcriptome profile of cotton leaves and established a hierarchical transcriptional network regulating cell-specific terpenoid production. We observed substantial expression levels of genes associated with the biosynthesis of both volatile terpenes (such as β-caryophyllene and β-myrcene) and non-volatile gossypol-type terpenoids in secretory glandular cells. Moreover, two novel transcription factors, namely GoHSFA4a and GoNAC42, are identified to function downstream of the Gossypium PIGMENT GLAND FORMATION genes. Both transcription factors could directly regulate the expression of terpenoid biosynthetic genes in secretory glandular cells in response to developmental and environmental stimuli. For convenient retrieval of the single-cell RNA sequencing data generated in this study, we developed a user-friendly web server . Our findings not only offer valuable insights into the precise regulation of terpenoid biosynthesis genes in cotton leaves but also provide potential targets for cotton breeding endeavors.

Time series single-cell transcriptional atlases reveal cell fate differentiation driven by light in Arabidopsis seedlings

Light serves as the energy source for plants as well as a signal for growth and development during their whole life cycle. Seedling de-etiolation is the most dramatic manifestation of light-regulated plant development processes, as massive reprogramming of the plant transcriptome occurs at this time. Although several studies have reported about organ-specific development and expression induced by light, a systematic analysis of cell-type-specific differentiation and the associated transcriptional regulation is still lacking. Here we obtained single-cell transcriptional atlases for etiolated, de-etiolating and light-grown Arabidopsis thaliana seedlings. Informative cells from shoot and root tissues were grouped into 48 different cell clusters and finely annotated using multiple markers. With the determination of comprehensive developmental trajectories, we demonstrate light modulation of cell fate determination during guard cell specialization and vasculature development. Comparison of expression atlases between wild type and the pifq mutant indicates that phytochrome-interacting factors (PIFs) are involved in distinct developmental processes in endodermal and stomatal lineage cells via controlling cell-type-specific expression of target genes. These results provide information concerning the light signalling networks at the cell-type resolution, improving our understanding of how light regulates plant development at the cell-type and genome-wide levels. The obtained information could serve as a valuable resource for comprehensively investigating the molecular mechanism of cell development and differentiation in response to light.

Resolving plant development in space and time with single-cell genomics

Single-cell genomics technologies are ushering in a new research era. In this review, we summarize the benefits and current challenges of using these technologies to probe the transcriptional regulation of plant development. In addition to profiling cells at a single snapshot in time, researchers have recently produced time-resolved datasets to map cell responses to stimuli. Live-imaging and spatial transcriptomic techniques are rapidly being adopted to link a cell's transcriptional profile with its spatial location within a tissue. Combining these technologies is a powerful spatiotemporal approach to investigate cell plasticity and developmental responses that contribute to plant resilience. Although there are hurdles to overcome, we conclude by discussing how single-cell genomics is poised to address developmental questions in the coming years.

Coexpression enhances cross-species integration of scRNA-seq across diverse plant species

Single-cell RNA sequencing is increasingly used to investigate cross-species differences driven by gene expression and cell-type composition in plants. However, the frequent expansion of plant gene families due to whole genome duplications makes identification of one-to-one orthologs difficult, complicating integration. Here, we demonstrate that coexpression can be used to identify non-orthologous gene pairs with proxy expression profiles, improving the performance of traditional integration methods and reducing barriers to integration across a diverse array of plant species.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

How to explore what is hidden? A review of techniques for vascular tissue expression profile analysis

The evolution of plants to efficiently transport water and assimilates over long distances is a major evolutionary success that facilitated their growth and colonization of land. Vascular tissues, namely xylem and phloem, are characterized by high specialization, cell heterogeneity, and diverse cell components. During differentiation and maturation, these tissues undergo an irreversible sequence of events, leading to complete protoplast degradation in xylem or partial degradation in phloem, enabling their undisturbed conductive function. Due to the unique nature of vascular tissue, and the poorly understood processes involved in xylem and phloem development, studying the molecular basis of tissue differentiation is challenging. In this review, we focus on methods crucial for gene expression research in conductive tissues, emphasizing the importance of initial anatomical analysis and appropriate material selection. We trace the expansion of molecular techniques in vascular gene expression studies and discuss the application of single-cell RNA sequencing, a high-throughput technique that has revolutionized transcriptomic analysis. We explore how single-cell RNA sequencing will enhance our knowledge of gene expression in conductive tissues.

NRTPredictor: identifying rice root cell state in single-cell RNA-seq via ensemble learning

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) measurements of gene expression show great promise for studying the cellular heterogeneity of rice roots. How precisely annotating cell identity is a major unresolved problem in plant scRNA-seq analysis due to the inherent high dimensionality and sparsity.

RESULTS

To address this challenge, we present NRTPredictor, an ensemble-learning system, to predict rice root cell stage and mine biomarkers through complete model interpretability. The performance of NRTPredictor was evaluated using a test dataset, with 98.01% accuracy and 95.45% recall. With the power of interpretability provided by NRTPredictor, our model recognizes 110 marker genes partially involved in phenylpropanoid biosynthesis. Expression patterns of rice root could be mapped by the above-mentioned candidate genes, showing the superiority of NRTPredictor. Integrated analysis of scRNA and bulk RNA-seq data revealed aberrant expression of Epidermis cell subpopulations in flooding, Pi, and salt stresses.

CONCLUSION

Taken together, our results demonstrate that NRTPredictor is a useful tool for automated prediction of rice root cell stage and provides a valuable resource for deciphering the rice root cellular heterogeneity and the molecular mechanisms of flooding, Pi, and salt stresses. Based on the proposed model, a free webserver has been established, which is available at https://www.cgris.net/nrtp .

Single cell RNA-seq in phytohormone signaling: a promising future

Phytohormone signaling regulates plant growth and development. Single cell RNA sequencing (scRNA-seq) provides unprecedented opportunities to decipher hormone-mediated spatiotemporal gene regulatory networks. In a recent study, Nolan et al. used time-series scRNA-seq to identify the cortex as a key site for brassinosteroid (BR)-mediated gene expression and revealed a signaling network during cell phase transition.

Cella: 3D data visualization for plant single-cell transcriptomics in Blender

AIMS

Recent advancements in single-cell transcriptomics have facilitated the possibility of acquiring vast amounts of data at single-cell resolution. This development has provided a broader and more comprehensive understanding of complex biological processes. The growing datasets require a visualization tool that transforms complex data into an intuitive representation. To address this challenge, we have utilized an open-source 3D software Blender to design Cella, a cell atlas visualization tool, which transforms data into 3D heatmaps that can be rendered into image libraries. Our tool is designed to support especially research on plant development.

DATA RESOURCES GENERATED

To validate our method, we have created a 3D model representing the Arabidopsis thaliana root meristem and mapped an existing single-cell RNA-seq dataset into the 3D model. This provided a user-friendly visual representation of the expression profiles of 21,489 genes from two perspectives (42,978 images).

UTILITY OF THE RESOURCE

This approach is not limited to single-cell RNA-seq data of the Arabidopsis root meristem. We provide detailed step-by-step instructions to generate 3D models and a script that can be customized to project data onto different tissues.

KEY RESULTS

Our tool provides a proof-of-concept method for how increasingly complex single-cell RNA-seq datasets can be visualized in a simple and cohesive manner.

Genome and transcriptome of Selaginella kraussiana reveal evolution of root apical meristems in vascular plants

The evolution of roots allowed vascular plants to adapt to land environments. Fossil evidence indicates that roots evolved independently in euphyllophytes (ferns and seed plants) and lycophytes, the two lineages of extant vascular plants. Based on a high-quality genome assembly, mRNA sequencing (mRNA-seq) data, and single-cell RNA-seq data for the lycophyte Selaginella kraussiana, we show that the two root origin events in lycophytes and euphyllophytes adopted partially similar molecular modules in the regulation of root apical meristem (RAM) development. In S. kraussiana, the RAM initiates from the rhizophore primordium guided by auxin and duplicates itself by dichotomous branching. The auxin signaling pathway directly upregulates euAINTEGUMENTAb (SkeuANTb), and then SkeuANTb directly promotes the expression of SkeuANTa and the WUSCHEL-RELATED HOMEOBOX13b (SkWOX13b) for RAM maintenance, partially similar to the molecular pathway involving the euANT-branch PLETHORA (AtPLT) genes and AtWOX5 in root initiation in the seed plant Arabidopsis thaliana. Other molecular modules, e.g., SHORT-ROOT and SCARECROW, also have partially similar expression patterns in the RAMs of S. kraussiana and A. thaliana. Overall, our study not only provides genome and transcriptome tools of S. kraussiana but also indicates the employment of some common molecular modules in RAMs during root origins in lycophytes and euphyllophytes.

Spatially resolved gene regulatory networks in Asian rice (Oryza sativa cv. Nipponbare) leaves

Transcriptome profiles in plants are heterogenous at every level of morphological organization. Even within organs, cells of the same type can have different patterns of gene expression depending on where they are positioned within tissues. This heterogeneity is associated with non-uniform distribution of biological processes within organs. The regulatory mechanisms that establish and sustain the spatial heterogeneity are unknown. Here, we identify regulatory modules that support functional specialization of different parts of Oryza sativa cv. Nipponbare leaves by leveraging transcriptome data, transcription factor binding motifs and global gene regulatory network prediction algorithms. We generated a global gene regulatory network in which we identified six regulatory modules that were active in different parts of the leaf. The regulatory modules were enriched for genes involved in spatially relevant biological processes, such as cell wall deposition, environmental sensing and photosynthesis. Strikingly, more than 86.9% of genes in the network were regulated by members of only five transcription factor families. We also generated targeted regulatory networks for the large MYB and bZIP/bHLH families to identify interactions that were masked in the global prediction. This analysis will provide a baseline for future single cell and array-based spatial transcriptome studies and for studying responses to environmental stress and demonstrates the extent to which seven coarse spatial transcriptome analysis can provide insight into the regulatory mechanisms supporting functional specialization within leaves.

Single-nucleus transcriptomes reveal spatiotemporal symbiotic perception and early response in Medicago

Establishing legume-rhizobial symbiosis requires precise coordination of complex responses in a time- and cell type-specific manner. Encountering Rhizobium, rapid changes of gene expression levels in host plants occur in the first few hours, which prepare the plants to turn off defence and form a symbiotic relationship with the microbes. Here, we applied single-nucleus RNA sequencing to characterize the roots of Medicago truncatula at 30 min, 6 h and 24 h after nod factor treatment. We found drastic global gene expression reprogramming at 30 min in the epidermis and cortex and most of these changes were restored at 6 h. Moreover, plant defence response genes are activated at 30 min and subsequently suppressed at 6 h in non-meristem cells. Only in the cortical cells but not in other cell types, we found the flavonoid synthase genes required to recruit rhizobia are highly expressed 30 min after inoculation with nod factors. A gene module enriched for symbiotic nitrogen fixation genes showed that MtFER (MtFERONIA) and LYK3 (LysM domain receptor-like kinase 3) share similar responses to symbiotic signals. We further found that MtFER can be phosphorylated by LYK3 and it participates in rhizobial symbiosis. Our results expand our understanding of dynamic spatiotemporal symbiotic responses at the single-cell level.