Single-cell resolution of lineage trajectories in the Arabidopsis stomatal lineage and developing leaf

Single-cell RNA-sequencing provides new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass leaf blades

As an important warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) flourishes in warm areas around the world due to the existence of the C4 photosynthetic pathway. However, how C4 photosynthesis operates in bermudagrass leaves is still poorly understood. In this study, we performed single-cell RNA-sequencing on 5296 cells from bermudagrass leaf blades. Eight cell clusters corresponding to mesophyll, bundle sheath, epidermis and vascular bundle cells were successfully identified using known cell marker genes. Expression profiling indicated that genes encoding NADP-dependent malic enzymes (NADP-MEs) were highly expressed in bundle sheath cells, whereas NAD-ME genes were weakly expressed in all cell types, suggesting C4 photosynthesis of bermudagrass leaf blades might be NADP-ME type rather than NAD-ME type. The results also indicated that starch synthesis-related genes showed preferential expression in bundle sheath cells, whereas starch degradation-related genes were highly expressed in mesophyll cells, which agrees with the observed accumulation of starch-filled chloroplasts in bundle sheath cells. Gene co-expression analysis further revealed that different families of transcription factors were co-expressed with multiple C4 photosynthesis-related genes, suggesting a complex transcription regulatory network of C4 photosynthesis might exist in bermudagrass leaf blades. These findings collectively provided new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass.

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.

ScRNA-seq reveals dark- and light-induced differentially expressed gene atlases of seedling leaves in Arachis hypogaea L

Although the regulatory mechanisms of dark and light-induced plant morphogenesis have been broadly investigated, the biological process in peanuts has not been systematically explored on single-cell resolution. Herein, 10 cell clusters were characterized using scRNA-seq-identified marker genes, based on 13 409 and 11 296 single cells from 1-week-old peanut seedling leaves grown under dark and light conditions. 6104 genes and 50 transcription factors (TFs) displayed significant expression patterns in distinct cell clusters, which provided gene resources for profiling dark/light-induced candidate genes. Further pseudo-time trajectory and cell cycle evidence supported that dark repressed the cell division and perturbed normal cell cycle, especially the PORA abundances correlated with 11 TFs highly enriched in mesophyll to restrict the chlorophyllide synthesis. Additionally, light repressed the epidermis cell developmental trajectory extending by inhibiting the growth hormone pathway, and 21 TFs probably contributed to the different genes transcriptional dynamic. Eventually, peanut AHL17 was identified from the profile of differentially expressed TFs, which encoded protein located in the nucleus promoted leaf epidermal cell enlargement when ectopically overexpressed in Arabidopsis through the regulatory phytohormone pathway. Overall, our study presents the different gene atlases in peanut etiolated and green seedlings, providing novel biological insights to elucidate light-induced leaf cell development at the single-cell level.

AraLeTA: An Arabidopsis leaf expression atlas across diurnal and developmental scales

Mature plant leaves are a composite of distinct cell types, including epidermal, mesophyll, and vascular cells. Notably, the proportion of these cells and the relative transcript concentrations within different cell types may change over time. While gene expression data at a single-cell level can provide cell-type-specific expression values, it is often too expensive to obtain these data for high-resolution time series. Although bulk RNA-seq can be performed in a high-resolution time series, RNA-seq using whole leaves measures average gene expression values across all cell types in each sample. In this study, we combined single-cell RNA-seq data with time-series data from whole leaves to assemble an atlas of cell-type-specific changes in gene expression over time for Arabidopsis (Arabidopsis thaliana). We inferred how the relative transcript concentrations of different cell types vary across diurnal and developmental timescales. Importantly, this analysis revealed 3 subgroups of mesophyll cells with distinct temporal profiles of expression. Finally, we developed tissue-specific gene networks that form a community resource: an Arabidopsis Leaf Time-dependent Atlas (AraLeTa). This allows users to extract gene networks that are confirmed by transcription factor-binding data and specific to certain cell types at certain times of day and at certain developmental stages. AraLeTa is available at https://regulatorynet.shinyapps.io/araleta/.

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.

Young guard cells function dynamically despite low mechanical anisotropy but gain efficiency during stomatal maturation in Arabidopsis thaliana

Stomata are pores at the leaf surface that enable gas exchange and transpiration. The signaling pathways that regulate the differentiation of stomatal guard cells and the mechanisms of stomatal pore formation have been characterized in Arabidopsis thaliana. However, the process by which stomatal complexes develop after pore formation into fully mature complexes is poorly understood. We tracked the morphogenesis of young stomatal complexes over time to establish characteristic geometric milestones along the path of stomatal maturation. Using 3D-nanoindentation coupled with finite element modeling of young and mature stomata, we found that despite having thicker cell walls than young guard cells, mature guard cells are more energy efficient with respect to stomatal opening, potentially attributable to the increased mechanical anisotropy of their cell walls and smaller changes in turgor pressure between the closed and open states. Comparing geometric changes in young and mature guard cells of wild-type and cellulose-deficient plants revealed that although cellulose is required for normal stomatal maturation, mechanical anisotropy appears to be achieved by the collective influence of cellulose and additional wall components. Together, these data elucidate the dynamic geometric and biomechanical mechanisms underlying the development process of stomatal maturation.

Stomatal maturomics: hunting genes regulating guard cell maturation and function formation from single-cell transcriptomes

Stomata play critical roles in gas exchange and immunity to pathogens. While many genes regulating early stomatal development up to the production of young guard cells (GCs) have been identified in Arabidopsis, much less is known about how young GCs develop into mature functional stomata. Here we perform a maturomics study on stomata, with "maturomics" defined as omics analysis of the maturation process of a tissue or organ. We develop an integrative scheme to analyze three public stomata-related single-cell RNA-seq datasets and identify a list of 586 genes that are specifically up-regulated in all three datasets during stomatal maturation and function formation. The list, termed sc_586, is enriched with known regulators of stomatal maturation and functions. To validate the reliability of the dataset, we selected two candidate G2-like transcription factor genes, MYS1 and MYS2, to investigate their roles in stomata. These two genes redundantly regulate the size and hoop rigidity of mature GCs, and the mys1 mys2 double mutants cause mature GCs with severe defects in regulating their stomatal apertures. Taken together, our results provide a valuable list of genes for studying GC maturation and function formation.

Stomatal evolution and plant adaptation to future climate

Global climate change is affecting plant photosynthesis and transpiration processes, as well as increasing weather extremes impacting socio-political and environmental events and decisions for decades to come. One major research challenge in plant biology and ecology is the interaction of photosynthesis with the environment. Stomata control plant gas exchange and their evolution was a crucial innovation that facilitated the earliest land plants to colonize terrestrial environments. Stomata couple homoiohydry, together with cuticles, intercellular gas space, with the endohydric water-conducting system, enabling plants to adapt and diversify across the planet. Plants control stomatal movement in response to environmental change through regulating guard cell turgor mediated by membrane transporters and signaling transduction. However, the origin, evolution, and active control of stomata remain controversial topics. We first review stomatal evolution and diversity, providing fossil and phylogenetic evidence of their origins. We summarize functional evolution of guard cell membrane transporters in the context of climate changes and environmental stresses. Our analyses show that the core signaling elements of stomatal movement are more ancient than stomata, while genes involved in stomatal development co-evolved de novo with the earliest stomata. These results suggest that novel stomatal development-specific genes were acquired during plant evolution, whereas genes regulating stomatal movement, especially cell signaling pathways, were inherited ancestrally and co-opted by dynamic functional differentiation. These two processes reflect the different adaptation strategies during land plant evolution.

Multifaceted role and regulation of aquaporins for efficient stomatal movements

Stomata are micropores on the leaf epidermis that allow carbon dioxide (CO2) uptake for photosynthesis at the expense of water loss through transpiration. Stomata coordinate the plant gas exchange of carbon and water with the atmosphere through their opening and closing dynamics. In the context of global climate change, it is essential to better understand the mechanism of stomatal movements under different environmental stimuli. Aquaporins (AQPs) are considered important regulators of stomatal movements by contributing to membrane diffusion of water, CO2 and hydrogen peroxide. This review compiles the most recent findings and discusses future directions to update our knowledge of the role of AQPs in stomatal movements. After highlighting the role of subsidiary cells (SCs), which contribute to the high water use efficiency of grass stomata, we explore the expression of AQP genes in guard cells and SCs. We then focus on the cellular regulation of AQP activity at the protein level in stomata. After introducing their post-translational modifications, we detail their trafficking as well as their physical interaction with various partners that regulate AQP subcellular dynamics towards and within specific regions of the cell membranes, such as microdomains and membrane contact sites.

Single-cell transcriptome profiling reveals the spatiotemporal distribution of triterpenoid saponin biosynthesis and transposable element activity in Gynostemma pentaphyllum shoot apexes and leaves

Gynostemma pentaphyllum (Thunb.) Makino is an important producer of dammarene-type triterpenoid saponins. These saponins (gypenosides) exhibit diverse pharmacological benefits such as anticancer, antidiabetic, and immunomodulatory effects, and have major potential in the pharmaceutical and health care industries. Here, we employed single-cell RNA sequencing (scRNA-seq) to profile the transcriptomes of more than 50,000 cells derived from G. pentaphyllum shoot apexes and leaves. Following cell clustering and annotation, we identified five major cell types in shoot apexes and four in leaves. Each cell type displayed substantial transcriptomic heterogeneity both within and between tissues. Examining gene expression patterns across various cell types revealed that gypenoside biosynthesis predominantly occurred in mesophyll cells, with heightened activity observed in shoot apexes compared to leaves. Furthermore, we explored the impact of transposable elements (TEs) on G. pentaphyllum transcriptomic landscapes. Our findings the highlighted the unbalanced expression of certain TE families across different cell types in shoot apexes and leaves, marking the first investigation of TE expression at the single-cell level in plants. Additionally, we observed dynamic expression of genes involved in gypenoside biosynthesis and specific TE families during epidermal and vascular cell development. The involvement of TE expression in regulating cell differentiation and gypenoside biosynthesis warrant further exploration. Overall, this study not only provides new insights into the spatiotemporal organization of gypenoside biosynthesis and TE activity in G. pentaphyllum shoot apexes and leaves but also offers valuable cellular and genetic resources for a deeper understanding of developmental and physiological processes at single-cell resolution in this species.

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.

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.

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.

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.

Identification of bZIP Transcription Factors That Regulate the Development of Leaf Epidermal Cells in Arabidopsis thaliana by Single-Cell RNA Sequencing

Epidermal cells are the main avenue for signal and material exchange between plants and the environment. Leaf epidermal cells primarily include pavement cells, guard cells, and trichome cells. The development and distribution of different epidermal cells are tightly regulated by a complex transcriptional regulatory network mediated by phytohormones, including jasmonic acid, and transcription factors. How the fate of leaf epidermal cells is determined, however, is still largely unknown due to the diversity of cell types and the complexity of their regulation. Here, we characterized the transcriptional profiles of epidermal cells in 3-day-old true leaves of Arabidopsis thaliana using single-cell RNA sequencing. We identified two genes encoding BASIC LEUCINE-ZIPPER (bZIP) transcription factors, namely bZIP25 and bZIP53, which are highly expressed in pavement cells and early-stage meristemoid cells. Densities of pavement cells and trichome cells were found to increase and decrease, respectively, in bzip25 and bzip53 mutants, compared with wild-type plants. This trend was more pronounced in the presence of jasmonic acid, suggesting that these transcription factors regulate the development of trichome cells and pavement cells in response to jasmonic acid.

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.

Single-cell transcriptome analysis dissects lncRNA-associated gene networks in Arabidopsis

The plant genome produces an extremely large collection of long noncoding RNAs (lncRNAs) that are generally expressed in a context-specific manner and have pivotal roles in regulation of diverse biological processes. Here, we mapped the transcriptional heterogeneity of lncRNAs and their associated gene regulatory networks at single-cell resolution. We generated a comprehensive cell atlas at the whole-organism level by integrative analysis of 28 published single-cell RNA sequencing (scRNA-seq) datasets from juvenile Arabidopsis seedlings. We then provided an in-depth analysis of cell-type-related lncRNA signatures that show expression patterns consistent with canonical protein-coding gene markers. We further demonstrated that the cell-type-specific expression of lncRNAs largely explains their tissue specificity. In addition, we predicted gene regulatory networks on the basis of motif enrichment and co-expression analysis of lncRNAs and mRNAs, and we identified putative transcription factors orchestrating cell-type-specific expression of lncRNAs. The analysis results are available at the single-cell-based plant lncRNA atlas database (scPLAD; https://biobigdata.nju.edu.cn/scPLAD/). Overall, this work demonstrates the power of integrative single-cell data analysis applied to plant lncRNA biology and provides fundamental insights into lncRNA expression specificity and associated gene regulation.

Leaf growth - complex regulation of a seemingly simple process

Understanding the underlying mechanisms of plant development is crucial to successfully steer or manipulate plant growth in a targeted manner. Leaves, the primary sites of photosynthesis, are vital organs for many plant species, and leaf growth is controlled by a tight temporal and spatial regulatory network. In this review, we focus on the genetic networks governing leaf cell proliferation, one major contributor to final leaf size. First, we provide an overview of six regulator families of leaf growth in Arabidopsis: DA1, PEAPODs, KLU, GRFs, the SWI/SNF complexes, and DELLAs, together with their surrounding genetic networks. Next, we discuss their evolutionary conservation to highlight similarities and differences among species, because knowledge transfer between species remains a big challenge. Finally, we focus on the increase in knowledge of the interconnectedness between these genetic pathways, the function of the cell cycle machinery as their central convergence point, and other internal and environmental cues.

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.

Genome-wide analysis of the TIFY family in Lycium and the negative regulation of stomatal development by LrJAZ2 gene

Stomata are ports that facilitate gas and water vapor exchange during plant photosynthesis and transpiration. Stomatal development is strictly regulated by endogenous hormone. Jasmonate, an important signal that modulates multiple physiological processes in plants, has been found to negatively regulate stomatal development in Arabidopsis thaliana, yet the molecular mechanisms underlying stomata development signaling remain to be understood. Jasmonate ZIM-domain (JAZ) proteins are the members of TIFY family and the key component of JA signaling pathway. Its function in stomatal development is unclear to data. Here, we screened out 24 TIFY family members against the genome of Lycium, and identified a JAZ member by combination analyses of evolutionary tree, cis-elements in promoter and gene expression patterns. Overexpression of this gene (LrJAZ2) in Lycium ruthenicum and Arabidopsis thaliana indicated LrJAZ2 negatively regulates stomatal development. Microscopic observations revealed that overexpression of LrJAZ2 negatively regulated stomatal development by decreasing stomatal density and index, which may lead to lower leaf transpiration rates. Transcriptome data indicated the overexpression of LrJAZ2 up-regulated the stomatal related genes such as LrERL2, LrPYL4, and down-regulated the LrSPCH. Collectively, our study found that LrJAZ2 is a key gene in stomatal development regulation in L. ruthenicum and provided new insights into the regulation of stomatal development.

Distinct guard cell-specific remodeling of chromatin accessibility during abscisic acid- and CO2-dependent stomatal regulation

In plants, epidermal guard cells integrate and respond to numerous environmental signals to control stomatal pore apertures, thereby regulating gas exchange. Chromatin structure controls transcription factor (TF) access to the genome, but whether large-scale chromatin remodeling occurs in guard cells during stomatal movements, and in response to the hormone abscisic acid (ABA) in general, remains unknown. Here, we isolate guard cell nuclei from Arabidopsis thaliana plants to examine whether the physiological signals, ABA and CO2 (carbon dioxide), regulate guard cell chromatin during stomatal movements. Our cell type-specific analyses uncover patterns of chromatin accessibility specific to guard cells and define cis-regulatory sequences supporting guard cell-specific gene expression. We find that ABA triggers extensive and dynamic chromatin remodeling in guard cells, roots, and mesophyll cells with clear patterns of cell type specificity. DNA motif analyses uncover binding sites for distinct TFs enriched in ABA-induced and ABA-repressed chromatin. We identify the Abscisic Acid Response Element (ABRE) Binding Factor (ABF) bZIP-type TFs that are required for ABA-triggered chromatin opening in guard cells and roots and implicate the inhibition of a clade of bHLH-type TFs in controlling ABA-repressed chromatin. Moreover, we demonstrate that ABA and CO2 induce distinct programs of chromatin remodeling, whereby elevated atmospheric CO2 had only minimal impact on chromatin dynamics. We provide insight into the control of guard cell chromatin dynamics and propose that ABA-induced chromatin remodeling primes the genome for abiotic stress resistance.

Single-cell transcriptome sequencing of plant leaf expressing anti-HER2 VHH-FcK cancer therapeutic protein

The transgenic plant is a promising strategy for the production of highly valuable biotherapeutic proteins such as recombinant vaccines and antibodies. To achieve an efficient level of protein production, codon sequences and expression cassette elements need to be optimized. However, the systematical expression of recombinant proteins in plant biomass can generally be controlled for the production of therapeutic proteins after the generation of transgenic plants. Without understanding the transgene expression patterns in plant tissue, it is difficult to enhance further production levels. In this study, single-cell RNA-sequencing (scRNA-seq) analysis of transgenic tobacco (Nicotiana tabacum) leaf, expressing an immunotherapeutic llama antibody against breast cancer, anti-HER2 VHH-Fc, was conducted to obtain data on the expression pattern of tissue-specific cells. These high-quality scRNA-seq data enabled the identification of gene expression patterns by cell types, which can be applied to select the best cell types or tissues for the high production of these recombinant antibodies. These data provide a foundation to elucidate the mechanisms that regulate the biosynthesis of recombinant proteins in N. tabacum.

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.

ABA guides stomatal proliferation and patterning through the EPF-SPCH signaling pathway in Arabidopsis thaliana

Adaptation to dehydration stress requires plants to coordinate environmental and endogenous signals to inhibit stomatal proliferation and modulate their patterning. The stress hormone abscisic acid (ABA) induces stomatal closure and restricts stomatal lineage to promote stress tolerance. Here, we report that mutants with reduced ABA levels, xer-1, xer-2 and aba2-2, developed stomatal clusters. Similarly, the ABA signaling mutant snrk2.2/2.3/2.6, which lacks core ABA signaling kinases, also displayed stomatal clusters. Exposure to ABA or inhibition of ABA catabolism rescued the increased stomatal density and spacing defects observed in xer and aba2-2, suggesting that basal ABA is required for correct stomatal density and spacing. xer-1 and aba2-2 displayed reduced expression of EPF1 and EPF2, and enhanced expression of SPCH and MUTE. Furthermore, ABA suppressed elevated SPCH and MUTE expression in epf2-1 and epf1-1, and partially rescued epf2-1 stomatal index and epf1-1 clustering defects. Genetic analysis demonstrated that XER acts upstream of the EPF2-SPCH pathway to suppress stomatal proliferation, and in parallel with EPF1 to ensure correct stomatal spacing. These results show that basal ABA and functional ABA signaling are required to fine-tune stomatal density and patterning.

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.

The stomatal fates: Understanding initiation and enforcement of stomatal cell fate transitions

In the stomatal lineage, repeated arcs of initiation, stem-cell proliferation, and terminal cell fate commitment are displayed on the surface of aerial organs. Over the past two decades, the core transcription and signaling elements that guide cell divisions, patterning, and fate transitions were defined. Here we highlight recent work that extends the core using a variety of cutting-edge techniques in different plant species. New work has discovered transcriptional circuits that initiate and reinforce stomatal fate transitions, while also enabling the lineage to interpret and respond to environmental inputs. Recent developments show that some key stomatal factors are more flexible or potentially even interchangeable, opening up avenues to explore stomatal fates and regulatory networks.

Stomatal cell fate commitment via transcriptional and epigenetic control: Timing is crucial

The formation of stomata presents a compelling model system for comprehending the initiation, proliferation, commitment and differentiation of de novo lineage-specific stem cells. Precise, timely and robust cell fate and identity decisions are crucial for the proper progression and differentiation of functional stomata. Deviations from this precise specification result in developmental abnormalities and nonfunctional stomata. However, the molecular underpinnings of timely cell fate commitment have just begun to be unravelled. In this review, we explore the key regulatory strategies governing cell fate commitment, emphasizing the distinctions between embryonic and postembryonic stomatal development. Furthermore, the interplay of transcription factors and cell cycle machineries is pivotal in specifying the transition into differentiation. We aim to synthesize recent studies utilizing single-cell as well as cell-type-specific transcriptomics, epigenomics and chromatin accessibility profiling to shed light on how master-regulatory transcription factors and epigenetic machineries mutually influence each other to drive fate commitment and maintenance.

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.

Distinct guard cell specific remodeling of chromatin accessibility during abscisic acid and CO2 dependent stomatal regulation

In plants, epidermal guard cells integrate and respond to numerous environmental signals to control stomatal pore apertures thereby regulating gas exchange. Chromatin structure controls transcription factor access to the genome, but whether large-scale chromatin remodeling occurs in guard cells during stomatal movements, and in response to the hormone abscisic acid (ABA) in general, remain unknown. Here we isolate guard cell nuclei from Arabidopsis thaliana plants to examine whether the physiological signals, ABA and CO2, regulate guard cell chromatin during stomatal movements. Our cell type specific analyses uncover patterns of chromatin accessibility specific to guard cells and define novel cis-regulatory sequences supporting guard cell specific gene expression. We find that ABA triggers extensive and dynamic chromatin remodeling in guard cells, roots, and mesophyll cells with clear patterns of cell-type specificity. DNA motif analyses uncover binding sites for distinct transcription factors enriched in ABA-induced and ABA-repressed chromatin. We identify the ABF/AREB bZIP-type transcription factors that are required for ABA-triggered chromatin opening in guard cells and implicate the inhibition of a set of bHLH-type transcription factors in controlling ABA-repressed chromatin. Moreover, we demonstrate that ABA and CO2 induce distinct programs of chromatin remodeling. We provide insight into the control of guard cell chromatin dynamics and propose that ABA-induced chromatin remodeling primes the genome for abiotic stress resistance.

Cell-type-specific responses to fungal infection in plants revealed by single-cell transcriptomics

Pathogen infection is a dynamic process. Here, we employ single-cell transcriptomics to investigate plant response heterogeneity. By generating an Arabidopsis thaliana leaf atlas encompassing 95,040 cells during infection by a fungal pathogen, Colletotrichum higginsianum, we unveil cell-type-specific gene expression, notably an enrichment of intracellular immune receptors in vasculature cells. Trajectory inference identifies cells that had different interactions with the invading fungus. This analysis divulges transcriptional reprogramming of abscisic acid signaling specifically occurring in guard cells, which is consistent with a stomatal closure dependent on direct contact with the fungus. Furthermore, we investigate the transcriptional plasticity of genes involved in glucosinolate biosynthesis in cells at the fungal infection sites, emphasizing the contribution of the epidermis-expressed MYB122 to disease resistance. This work underscores spatially dynamic, cell-type-specific plant responses to a fungal pathogen and provides a valuable resource that supports in-depth investigations of plant-pathogen interactions.

Improving photosynthetic efficiency toward food security: Strategies, advances, and perspectives

Photosynthesis in crops and natural vegetation allows light energy to be converted into chemical energy and thus forms the foundation for almost all terrestrial trophic networks on Earth. The efficiency of photosynthetic energy conversion plays a crucial role in determining the portion of incident solar radiation that can be used to generate plant biomass throughout a growth season. Consequently, alongside the factors such as resource availability, crop management, crop selection, maintenance costs, and intrinsic yield potential, photosynthetic energy use efficiency significantly influences crop yield. Photosynthetic efficiency is relevant to sustainability and food security because it affects water use efficiency, nutrient use efficiency, and land use efficiency. This review focuses specifically on the potential for improvements in photosynthetic efficiency to drive a sustainable increase in crop yields. We discuss bypassing photorespiration, enhancing light use efficiency, harnessing natural variation in photosynthetic parameters for breeding purposes, and adopting new-to-nature approaches that show promise for achieving unprecedented gains in photosynthetic efficiency.

Exploring the identity of individual plant cells in space and time

In recent years, single-cell genomics, coupled to imaging techniques, have become the state-of-the-art approach for characterising biological systems. In plant sciences, a variety of tissues and species have been profiled, providing an enormous quantity of data on cell identity at an unprecedented resolution, but what biological insights can be gained from such data sets? Using recently published studies in plant sciences, we will highlight how single-cell technologies have enabled a better comprehension of tissue organisation, cell fate dynamics in development or in response to various stimuli, as well as identifying key transcriptional regulators of cell identity. We discuss the limitations and technical hurdles to overcome, as well as future directions, and the promising use of single-cell omics to understand, predict, and manipulate plant development and physiology.

Experimental validation of the mechanism of stomatal development diversification

Stomata are the structures responsible for gas exchange in plants. The established framework for stomatal development is based on the model plant Arabidopsis, but diverse patterns of stomatal development have been observed in other plant lineages and species. The molecular mechanisms behind these diversified patterns are still poorly understood. We recently proposed a model for the molecular mechanisms of the diversification of stomatal development based on the genus Callitriche (Plantaginaceae), according to which a temporal shift in the expression of key stomatal transcription factors SPEECHLESS and MUTE leads to changes in the behavior of meristemoids (stomatal precursor cells). In the present study, we genetically manipulated Arabidopsis to test this model. By altering the timing of MUTE expression, we successfully generated Arabidopsis plants with early differentiation or prolonged divisions of meristemoids, as predicted by the model. The epidermal morphology of the generated lines resembled that of species with prolonged or no meristemoid divisions. Thus, the evolutionary process can be reproduced by varying the SPEECHLESS to MUTE transition. We also observed unexpected phenotypes, which indicated the participation of additional factors in the evolution of the patterns observed in nature. This study provides novel experimental insights into the diversification of meristemoid behaviors.

Application of single-cell multi-omics approaches in horticulture research

Cell heterogeneity shapes the morphology and function of various tissues and organs in multicellular organisms. Elucidation of the differences among cells and the mechanism of intercellular regulation is essential for an in-depth understanding of the developmental process. In recent years, the rapid development of high-throughput single-cell transcriptome sequencing technologies has influenced the study of plant developmental biology. Additionally, the accuracy and sensitivity of tools used to study the epigenome and metabolome have significantly increased, thus enabling multi-omics analysis at single-cell resolution. Here, we summarize the currently available single-cell multi-omics approaches and their recent applications in plant research, review the single-cell based studies in fruit, vegetable, and ornamental crops, and discuss the potential of such approaches in future horticulture research.

Arabidopsis stomatal lineage cells establish bipolarity and segregate differential signaling capacity to regulate stem cell potential

Cell polarity combined with asymmetric cell divisions (ACDs) generates cellular diversity. In the Arabidopsis stomatal lineage, a single cortical polarity domain marked by BASL orients ACDs and is segregated to the larger daughter to enforce cell fate. We discovered a second, oppositely positioned polarity domain defined by OCTOPUS-LIKE (OPL) proteins, which forms prior to ACD and is segregated to the smaller (meristemoid) daughter. Genetic and misexpression analyses show that OPLs promote meristemoid-amplifying divisions and delay stomatal fate progression. Polarity mediates OPL segregation into meristemoids but is not required for OPL function. OPL localization and activity are largely independent of other stomatal polarity genes and of the brassinosteroid signaling components associated with OPLs in other contexts. While OPLs are unique to seed plants, ectopic expression in the liverwort Marchantia suppressed epidermal fate progression, suggesting that OPLs engage ancient and broadly conserved pathways to regulate cell division and cell fate.

A cell size threshold triggers commitment to stomatal fate in Arabidopsis

How flexible developmental programs integrate information from internal and external factors to modulate stem cell behavior is a fundamental question in developmental biology. Cells of the Arabidopsis stomatal lineage modify the balance of stem cell proliferation and differentiation to adjust the size and cell type composition of mature leaves. Here, we report that meristemoids, one type of stomatal lineage stem cell, trigger the transition from asymmetric self-renewing divisions to commitment and terminal differentiation by crossing a critical cell size threshold. Through computational simulation, we demonstrate that this cell size-mediated transition allows robust, yet flexible termination of stem cell proliferation, and we observe adjustments in the number of divisions before the differentiation threshold under several genetic manipulations. We experimentally evaluate several mechanisms for cell size sensing, and our data suggest that this stomatal lineage transition is dependent on a nuclear factor that is sensitive to DNA content.

scRNA-seq Reveals the Mechanism of Fatty Acid Desaturase 2 Mutation to Repress Leaf Growth in Peanut (Arachis hypogaea L.)

Fatty Acid Desaturase 2 (FAD2) controls the conversion of oleic acids into linoleic acids. Mutations in FAD2 not only increase the high-oleic content, but also repress the leaf growth. However, the mechanism by which FAD2 regulates the growth pathway has not been elucidated in peanut leaves with single-cell resolution. In this study, we isolated fad2 mutant leaf protoplast cells to perform single-cell RNA sequencing. Approximately 24,988 individual cells with 10,249 expressed genes were classified into five major cell types. A comparative analysis of 3495 differentially expressed genes (DEGs) in distinct cell types demonstrated that fad2 inhibited the expression of the cytokinin synthesis gene LOG in vascular cells, thereby repressing leaf growth. Further, pseudo-time trajectory analysis indicated that fad2 repressed leaf cell differentiation, and cell-cycle evidence displayed that fad2 perturbed the normal cell cycle to induce the majority of cells to drop into the S phase. Additionally, important transcription factors were filtered from the DEG profiles that connected the network involved in high-oleic acid accumulation (WRKY6), activated the hormone pathway (WRKY23, ERF109), and potentially regulated leaf growth (ERF6, MYB102, WRKY30). Collectively, our study describes different gene atlases in high-oleic and normal peanut seedling leaves, providing novel biological insights to elucidate the molecular mechanism of the high-oleic peanut-associated agronomic trait at the single-cell level.

scPlant: A versatile framework for single-cell transcriptomic data analysis in plants

Single-cell transcriptomics has been fully embraced in plant biological research and is revolutionizing our understanding of plant growth, development, and responses to external stimuli. However, single-cell transcriptomic data analysis in plants is not trivial, given that there is currently no end-to-end solution and that integration of various bioinformatics tools involves a large number of required dependencies. Here, we present scPlant, a versatile framework for exploring plant single-cell atlases with minimum input data provided by users. The scPlant pipeline is implemented with numerous functions for diverse analytical tasks, ranging from basic data processing to advanced demands such as cell-type annotation and deconvolution, trajectory inference, cross-species data integration, and cell-type-specific gene regulatory network construction. In addition, a variety of visualization tools are bundled in a built-in Shiny application, enabling exploration of single-cell transcriptomic data on the fly.

Mass spectrometry imaging and single-cell transcriptional profiling reveal the tissue-specific regulation of bioactive ingredient biosynthesis in Taxus leaves

Taxus leaves provide the raw industrial materials for taxol, a natural antineoplastic drug widely used in the treatment of various cancers. However, the precise distribution, biosynthesis, and transcriptional regulation of taxoids and other active components in Taxus leaves remain unknown. Matrix-assisted laser desorption/ionization-mass spectrometry imaging analysis was used to visualize various secondary metabolites in leaf sections of Taxus mairei, confirming the tissue-specific accumulation of different active metabolites. Single-cell sequencing was used to produce expression profiles of 8846 cells, with a median of 2352 genes per cell. Based on a series of cluster-specific markers, cells were grouped into 15 clusters, suggesting a high degree of cell heterogeneity in T. mairei leaves. Our data were used to create the first Taxus leaf metabolic single-cell atlas and to reveal spatial and temporal expression patterns of several secondary metabolic pathways. According to the cell-type annotation, most taxol biosynthesis genes are expressed mainly in leaf mesophyll cells; phenolic acid and flavonoid biosynthesis genes are highly expressed in leaf epidermal cells (including the stomatal complex and guard cells); and terpenoid and steroid biosynthesis genes are expressed specifically in leaf mesophyll cells. A number of novel and cell-specific transcription factors involved in secondary metabolite biosynthesis were identified, including MYB17, WRKY12, WRKY31, ERF13, GT_2, and bHLH46. Our research establishes the transcriptional landscape of major cell types in T. mairei leaves at a single-cell resolution and provides valuable resources for studying the basic principles of cell-type-specific regulation of secondary metabolism.

Cell type-specific attenuation of brassinosteroid signaling precedes stomatal asymmetric cell division

In Arabidopsis thaliana, brassinosteroid (BR) signaling and stomatal development are connected through the SHAGGY/GSK3-like kinase BR INSENSITIVE2 (BIN2). BIN2 is a key negative regulator of BR signaling but it plays a dual role in stomatal development. BIN2 promotes or restricts stomatal asymmetric cell division (ACD) depending on its subcellular localization, which is regulated by the stomatal lineage-specific scaffold protein POLAR. BRs inactivate BIN2, but how they govern stomatal development remains unclear. Mapping the single-cell transcriptome of stomatal lineages after triggering BR signaling with either exogenous BRs or the specific BIN2 inhibitor, bikinin, revealed that the two modes of BR signaling activation generate spatiotemporally distinct transcriptional responses. We established that BIN2 is always sensitive to the inhibitor but, when in a complex with POLAR and its closest homolog POLAR-LIKE1, it becomes protected from BR-mediated inactivation. Subsequently, BR signaling in ACD precursors is attenuated, while it remains active in epidermal cells devoid of scaffolds and undergoing differentiation. Our study demonstrates how scaffold proteins contribute to cellular signal specificity of hormonal responses in plants.

Cell Fate Programming by Transcription Factors and Epigenetic Machinery in Stomatal Development

The development of multi-cellular organisms requires coordinated changes in gene expression that are often mediated by the interaction between transcription factors (TFs) and their corresponding cis-regulatory elements (CREs). During development and differentiation, the accessibility of CREs is dynamically modulated by the epigenome. How the epigenome, CREs and TFs together exert control over cell fate commitment remains to be fully understood. In the Arabidopsis leaf epidermis, meristemoids undergo a series of stereotyped cell divisions, then switch fate to commit to stomatal differentiation. Newly created or reanalyzed scRNA-seq and ChIP-seq data confirm that stomatal development involves distinctive phases of transcriptional regulation and that differentially regulated genes are bound by the stomatal basic-helix-loop-helix (bHLH) TFs. Targets of the bHLHs often reside in repressive chromatin before activation. MNase-seq evidence further suggests that the repressive state can be overcome and remodeled upon activation by specific stomatal bHLHs. We propose that chromatin remodeling is mediated through the recruitment of a set of physical interactors that we identified through proximity labeling - the ATPase-dependent chromatin remodeling SWI/SNF complex and the histone acetyltransferase HAC1. The bHLHs and chromatin remodelers localize to overlapping genomic regions in a hierarchical order. Furthermore, plants with stage-specific knock-down of the SWI/SNF components or HAC1 fail to activate specific bHLH targets and display stomatal development defects. Together these data converge on a model for how stomatal TFs and epigenetic machinery cooperatively regulate transcription and chromatin remodeling during progressive fate specification.

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

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

Analysis of meristems and plant regeneration at single-cell resolution

Rapid development of high-throughput single-cell RNA sequencing (scRNA-seq) technologies offers exciting opportunities to reveal new and rare cell types, previously hidden cell states, and continuous developmental trajectories. In this review, we first illustrate the ways in which scRNA-seq enables researchers to distinguish between distinct plant cell populations, delineate cell cycle continuums, and infer continuous differentiation trajectories of diverse cell types in shoots, roots, and floral and vascular meristems with unprecedented resolution. We then highlight the emerging power of scRNA-seq to dissect cell heterogeneity in regenerating tissues and uncover the cellular basis of cell reprogramming and stem cell commitment during plant regeneration. We conclude by discussing related outstanding questions in the field.

Single-cell profiling of Arabidopsis leaves to Pseudomonas syringae infection

Plant response to pathogen infection varies within a leaf, yet this heterogeneity is not well resolved. We expose Arabidopsis to Pseudomonas syringae or mock treatment and profile >11,000 individual cells using single-cell RNA sequencing. Integrative analysis of cell populations from both treatments identifies distinct pathogen-responsive cell clusters exhibiting transcriptional responses ranging from immunity to susceptibility. Pseudotime analyses through pathogen infection reveals a continuum of disease progression from an immune to a susceptible state. Confocal imaging of promoter-reporter lines for transcripts enriched in immune cell clusters shows expression surrounding substomatal cavities colonized or in close proximity to bacterial colonies, suggesting that cells within immune clusters represent sites of early pathogen invasion. Susceptibility clusters exhibit more general localization and are highly induced at later stages of infection. Overall, our work shows cellular heterogeneity within an infected leaf and provides insight into plant differential response to infection at a single-cell level.

Signaling in plant development and immunity through the lens of the stomata

The proper development and function of stomata - turgor-driven valves for efficient gas-exchange and water control - impact plant survival and productivity. It has become apparent that various receptor kinases regulate stomatal development and immunity. Although stomatal development and immunity occur over different cellular time scales, their signaling components and regulatory modules are strikingly similar, and often shared. In this review, we survey the current knowledge of stomatal development and immunity signaling components, and provide a synthesis and perspectives on the key concepts to further understand the conservation and specificity of these two signaling pathways.

Single-Cell RNA sequencing of leaf sheath cells reveals the mechanism of rice resistance to brown planthopper (Nilaparvata lugens)

The brown planthopper (BPH) (Nilaparvata lugens) sucks rice sap causing leaves to turn yellow and wither, often leading to reduced or zero yields. Rice co-evolved to resist damage by BPH. However, the molecular mechanisms, including the cells and tissues, involved in the resistance are still rarely reported. Single-cell sequencing technology allows us to analyze different cell types involved in BPH resistance. Here, using single-cell sequencing technology, we compared the response offered by the leaf sheaths of the susceptible (TN1) and resistant (YHY15) rice varieties to BPH (48 hours after infestation). We found that the 14,699 and 16,237 cells (identified via transcriptomics) in TN1 and YHY15 could be annotated using cell-specific marker genes into nine cell-type clusters. The two rice varieties showed significant differences in cell types (such as mestome sheath cells, guard cells, mesophyll cells, xylem cells, bulliform cells, and phloem cells) in the rice resistance mechanism to BPH. Further analysis revealed that although mesophyll, xylem, and phloem cells are involved in the BPH resistance response, the molecular mechanism used by each cell type is different. Mesophyll cell may regulate the expression of genes related to vanillin, capsaicin, and ROS production, phloem cell may regulate the cell wall extension related genes, and xylem cell may be involved in BPH resistance response by controlling the expression of chitin and pectin related genes. Thus, rice resistance to BPH is a complicated process involving multiple insect resistance factors. The results presented here will significantly promote the investigation of the molecular mechanisms underlying the resistance of rice to insects and accelerate the breeding of insect-resistant rice varieties.

Deciphering plant cell-cell communications using single-cell omics data

Plants have various cell types that respond to different environmental factors, and cell-cell communication is the fundamental process that controls these plant responses. The emergence of single-cell techniques provides opportunities to explore features unique to each cell type and construct a comprehensive cell-cell communication (CCC) network. Although the most current successes of CCC inference were achieved in animal research, computational methods can also be directly applied to plants. This review describes the current major models for cell-cell communication inference and summarizes the computational tools based on single-cell omics datasets. In addition, we discuss the limitations of plant cell-cell communication research and propose new directions to expand the field in meaningful ways.

Single-cell multi-omics in the medicinal plant Catharanthus roseus

Advances in omics technologies now permit the generation of highly contiguous genome assemblies, detection of transcripts and metabolites at the level of single cells and high-resolution determination of gene regulatory features. Here, using a complementary, multi-omics approach, we interrogated the monoterpene indole alkaloid (MIA) biosynthetic pathway in Catharanthus roseus, a source of leading anticancer drugs. We identified clusters of genes involved in MIA biosynthesis on the eight C. roseus chromosomes and extensive gene duplication of MIA pathway genes. Clustering was not limited to the linear genome, and through chromatin interaction data, MIA pathway genes were present within the same topologically associated domain, permitting the identification of a secologanin transporter. Single-cell RNA-sequencing revealed sequential cell-type-specific partitioning of the leaf MIA biosynthetic pathway that, when coupled with a single-cell metabolomics approach, permitted the identification of a reductase that yields the bis-indole alkaloid anhydrovinblastine. We also revealed cell-type-specific expression in the root MIA pathway.