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

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

As sessile organisms, plants constantly face challenges from the external environment. In order to meet these challenges and survive, plants have evolved a set of sophisticated adaptation strategies, including changes in leaf morphology and epidermal cell development. These developmental patterns are regulated by both light and hormonal signaling pathways. However, our mechanistic understanding of the role of these signaling pathways in regulating plant response to environmental stress is still very limited. By applying single-cell RNA-Seq, we determined the expression pattern of PHYTOCHROME INTERACTING FACTOR (PIF) 1, PIF3, PIF4, and PIF5 genes in leaf epidermal pavement cells (PCs) and guard cells (GCs). PCs and GCs are very sensitive to environmental stress, and our previous research suggests that these PIFs may be involved in regulating the development of PCs, GCs, and leaf morphology under environmental stress. Growth analysis showed that pif1/3/4/5 quadruple mutant maintained tolerance to drought and salt stress, and the length to width ratio of leaves and petiole length under normal growth conditions were similar to those of wild-type (WT) plants under drought and salt treatment. Analysis of the developmental patterns of PCs and GCs, and whole leaf morphology, further confirmed that these PIFs may be involved in mediating the development of epidermal cells under drought and salt stress, likely by regulating the expression of MUTE and TOO MANY MOUTHS (TMM) genes. These results provide new insights into the molecular mechanism of plant adaptation to adverse growth environments.

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

Epidermal cells are the main avenue for signal and material exchange between plants and the environment. Leaf epidermal cells primarily include pavement cells (PCs), guard cells, and trichomes cells (TCs), which differentiate from protodermal cells or meristemoids. The development and distribution of different epidermal cells are tightly regulated by a complex transcriptional regulatory network mediated by phytohormones, including jasmonic acid (JA), and transcription factors. Understanding 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 its 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 the bZIP25 and bZIP53, which are highly expressed in PCs and early-stage meristemoid cells. Densities of PCs and TCs 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 JA, suggesting that these transcription factors regulate the development of TCs and PCs in response to JA.IN A NUTSHELLBackgroundLeaf epidermal cells, comprised of trichome cells (TCs), guard cells (GCs), and pavement cells (PCs), are responsible for exchanging materials and information between plants and the surrounding aerial environment. Many genes have been identified in Arabidopsis thaliana and confirmed to be involved in the initiation and differentiation of TCs and PCs. The fate determination of TCs and PCs is tightly regulated by positive and negative regulators at the cellular level. The precise underlying molecular mechanisms responsible for the fate determination of TCs and PCs, however, are still unclear at this time.QuestionWhat are the transcriptomic profiles of different leaf epidermal cell types? Can we dissect the genes that are specifically expressed in certain epidermal cell types? What kinds of transcription factors are involved in regulating the fate determination of TCs and PCs?FindingsWe performed single cell RNA-seq to investigate the transcriptomic profiles of different leaf epidermal cell types and identified differentially expressed genes in each cell type. We found that genes that are involved in jasmonic acid signaling are highly expressed in early-stage meristemoid (EM) cells which can act as the precursor of PCs and perhaps of TCs. To investigate the regulatory mechanisms underlying EM development, we identified the transcription factors (TFs) in EM cells and found that two bZIP TF genes, bZIP25 and bZIP53, are highly expressed in EMs. Further analyses of these two genes using both loss-of-function and gain-of-function approaches indicated that bZIP25 and bZIP53 are functionally involved in promoting trichome formation but inhibit pavement cell development in response to jasmonic acid.Next stepsBesides of bZIP25 and bZIP53, we also identified other key genes, for example FES1B, in leaf epidermal cells. Our next step will be to explore the regulation of other key genes involved in the fate determination of different cell types in leaf epidermis.