Global Dynamic Molecular Profiling of Stomatal Lineage Cell Development by Single-Cell RNA Sequencing

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.

Mitochondrial ATP Synthase beta-Subunit Affects Plastid Retrograde Signaling in Arabidopsis

Plastid retrograde signaling plays a key role in coordinating the expression of plastid genes and photosynthesis-associated nuclear genes (PhANGs). Although plastid retrograde signaling can be substantially compromised by mitochondrial dysfunction, it is not yet clear whether specific mitochondrial factors are required to regulate plastid retrograde signaling. Here, we show that mitochondrial ATP synthase beta-subunit mutants with decreased ATP synthase activity are impaired in plastid retrograde signaling in Arabidopsis thaliana. Transcriptome analysis revealed that the expression levels of PhANGs were significantly higher in the mutants affected in the AT5G08670 gene encoding the mitochondrial ATP synthase beta-subunit, compared to wild-type (WT) seedlings when treated with lincomycin (LIN) or norflurazon (NF). Further studies indicated that the expression of nuclear genes involved in chloroplast and mitochondrial retrograde signaling was affected in the AT5G08670 mutant seedlings treated with LIN. These changes might be linked to the modulation of some transcription factors (TFs), such as LHY (Late Elongated Hypocotyl), PIF (Phytochrome-Interacting Factors), MYB, WRKY, and AP2/ERF (Ethylene Responsive Factors). These findings suggest that the activity of mitochondrial ATP synthase significantly influences plastid retrograde signaling.

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.

Harnessing spatial transcriptomics for advancing plant regeneration research

Song et al. utilized spatial transcriptomics to study the molecular characteristics of various cells - such as shoot primordia and chlorenchyma cells - in tomato callus during shoot regeneration. This research enhances our knowledge of shoot regeneration and demonstrates the potential of spatial transcriptomics in advancing plant biology.

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.

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.

Regulation of stomatal development by epidermal, subepidermal and long-distance signals

Plant leaves consist of three layers, including epidermis, mesophyll and vascular tissues. Their development is meticulously orchestrated. Stomata are the specified structures on the epidermis for uptake of carbon dioxide (CO2) while release of water vapour and oxygen (O2), and thus play essential roles in regulation of plant photosynthesis and water use efficiency. To function efficiently, stomatal formation must coordinate with the development of other epidermal cell types, such as pavement cell and trichome, and tissues of other layers, such as mesophyll and leaf vein. This review summarizes the regulation of stomatal development in three dimensions (3D). In the epidermis, specific stomatal transcription factors determine cell fate transitions and also activate a ligand-receptor- MITOGEN-ACTIVATED PROTEIN KINASE (MAPK) signaling for ensuring proper stomatal density and patterning. This forms the core regulation network of stomatal development, which integrates various environmental cues and phytohormone signals to modulate stomatal production. Under the epidermis, mesophyll, endodermis of hypocotyl and inflorescence stem, and veins in grasses secrete mobile signals to influence stomatal formation in the epidermis. In addition, long-distance signals which may include phytohormones, RNAs, peptides and proteins originated from other plant organs modulate stomatal development, enabling plants to systematically adapt to the ever changing environment.

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.

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.

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.

Systemic stomatal responses in plants: Coordinating development, stress, and pathogen defense under a changing climate

To successfully survive, develop, grow and reproduce, multicellular organisms must coordinate their molecular, physiological, developmental and metabolic responses among their different cells and tissues. This process is mediated by cell-to-cell, vascular and/or volatile communication, and involves electric, chemical and/or hydraulic signals. Within this context, stomata serve a dual role by coordinating their responses to the environment with their neighbouring cells at the epidermis, but also with other stomata present on other parts of the plant. As stomata represent one of the most important conduits between the plant and its above-ground environment, as well as directly affect photosynthesis, respiration and the hydraulic status of the plant by controlling its gas and vapour exchange with the atmosphere, coordinating the overall response of stomata within and between different leaves and tissues plays a cardinal role in plant growth, development and reproduction. Here, we discuss different examples of local and systemic stomatal coordination, the different signalling pathways that mediate them, and the importance of systemic stomatal coordination to our food supply, ecosystems and weather patterns, under our changing climate. We further discuss the potential biotechnological implications of regulating systemic stomatal responses for enhancing agricultural productivity in a warmer and CO2 -rich environment.

Microfluidic platform for omics analysis on single cells with diverse morphology and size: A review

BACKGROUND

Microfluidic techniques have emerged as powerful tools in single-cell research, facilitating the exploration of omics information from individual cells. Cell morphology is crucial for gene expression and physiological processes. However, there is currently a lack of integrated analysis of morphology and single-cell omics information. A critical challenge remains: what platform technologies are the best option to decode omics data of cells that are complex in morphology and size?

RESULTS

This review highlights achievements in microfluidic-based single-cell omics and isolation of cells based on morphology, along with other cell sorting methods based on physical characteristics. Various microfluidic platforms for single-cell isolation are systematically presented, showcasing their diversity and adaptability. The discussion focuses on microfluidic devices tailored to the distinct single-cell isolation requirements in plants and animals, emphasizing the significance of considering cell morphology and cell size in optimizing single-cell omics strategies. Simultaneously, it explores the application of microfluidic single-cell sorting technologies to single-cell sequencing, aiming to effectively integrate information about cell shape and size.

SIGNIFICANCE AND NOVELTY

The novelty lies in presenting a comprehensive overview of recent accomplishments in microfluidic-based single-cell omics, emphasizing the integration of different microfluidic platforms and their implications for cell morphology-based isolation. By underscoring the pivotal role of the specialized morphology of different cells in single-cell research, this review provides robust support for delving deeper into the exploration of single-cell omics data.

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.

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.

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.

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.

Ameliorating the effects of multiple stresses on agronomic traits in crops: modern biotechnological and omics approaches

While global climate change poses a significant environmental threat to agriculture, the increasing population is another big challenge to food security. To address this, developing crop varieties with increased productivity and tolerance to biotic and abiotic stresses is crucial. Breeders must identify traits to ensure higher and consistent yields under inconsistent environmental challenges, possess resilience against emerging biotic and abiotic stresses and satisfy customer demands for safer and more nutritious meals. With the advent of omics-based technologies, molecular tools are now integrated with breeding to understand the molecular genetics of genotype-based traits and develop better climate-smart crops. The rapid development of omics technologies offers an opportunity to generate novel datasets for crop species. Identifying genes and pathways responsible for significant agronomic traits has been made possible by integrating omics data with genetic and phenotypic information. This paper discusses the importance and use of omics-based strategies, including genomics, transcriptomics, proteomics and phenomics, for agricultural and horticultural crop improvement, which aligns with developing better adaptability in these crop species to the changing climate conditions.

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 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.

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.

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.

Single-cell RNA landscape of the special fiber initiation process in Bombax ceiba

As a new source of natural fibers, the Bombax ceiba tree can provide thin, light, extremely soft and warm fiber material for the textile industry. Natural fibers are an ideal model system for studying cell growth and differentiation, but the molecular mechanisms that regulate fiber initiation are not fully understood. In B. ceiba, we found that fiber cells differentiate from the epidermis of the inner ovary wall. Each initiated cell then divides into a cluster of fiber cells that eventually develop into mature fibers, a process very different from the classical fiber initiation process of cotton. We used high-throughput single-cell RNA sequencing (scRNA-seq) to examine the special characteristics of fiber initiation in B. ceiba. A total of 15 567 high-quality cells were identified from the inner wall of the B. ceiba ovary, and 347 potential marker genes for fiber initiation cell types were identified. Two major cell types, initiated fiber cells and epidermal cells, were identified and verified by RNA in situ hybridization. A developmental trajectory analysis was used to reconstruct the process of fiber cell differentiation in B. ceiba. Comparative analysis of scRNA-seq data from B. ceiba and cotton (Gossypium hirsutum) confirmed that the additional cell division process in B. ceiba is a novel species-specific mechanism for fiber cell development. Candidate genes and key regulators that may contribute to fiber cell differentiation and division in B. ceiba were identified. This work reveals gene expression signatures during B. ceiba fiber initiation at a single-cell resolution, providing a new strategy and viewpoint for investigation of natural fiber cell differentiation and development.

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.

Diel dynamics of multi-omics in elkhorn fern provide new insights into weak CAM photosynthesis

Crassulacean acid metabolism (CAM) has high water-use efficiency (WUE) and is widely recognized to have evolved from C3 photosynthesis. Different plant lineages have convergently evolved CAM, but the molecular mechanism that underlies C3-to-CAM evolution remains to be clarified. Platycerium bifurcatum (elkhorn fern) provides an opportunity to study the molecular changes underlying the transition from C3 to CAM photosynthesis because both modes of photosynthesis occur in this species, with sporotrophophyll leaves (SLs) and cover leaves (CLs) performing C3 and weak CAM photosynthesis, respectively. Here, we report that the physiological and biochemical attributes of CAM in weak CAM-performing CLs differed from those in strong CAM species. We investigated the diel dynamics of the metabolome, proteome, and transcriptome in these dimorphic leaves within the same genetic background and under identical environmental conditions. We found that multi-omic diel dynamics in P. bifurcatum exhibit both tissue and diel effects. Our analysis revealed temporal rewiring of biochemistry relevant to the energy-producing pathway (TCA cycle), CAM pathway, and stomatal movement in CLs compared with SLs. We also confirmed that PHOSPHOENOLPYRUVATE CARBOXYLASE KINASE (PPCK) exhibits convergence in gene expression among highly divergent CAM lineages. Gene regulatory network analysis identified candidate transcription factors regulating the CAM pathway and stomatal movement. Taken together, our results provide new insights into weak CAM photosynthesis and new avenues for CAM bioengineering.

Single-cell RNA sequencing profiles reveal cell type-specific transcriptional regulation networks conditioning fungal invasion in maize roots

Stalk rot caused by Fusarium verticillioides (Fv) is one of the most destructive diseases in maize production. The defence response of root system to Fv invasion is important for plant growth and development. Dissection of root cell type-specific response to Fv infection and its underlying transcription regulatory networks will aid in understanding the defence mechanism of maize roots to Fv invasion. Here, we reported the transcriptomes of 29 217 single cells derived from root tips of two maize inbred lines inoculated with Fv and mock condition, and identified seven major cell types with 21 transcriptionally distinct cell clusters. Through the weighted gene co-expression network analysis, we identified 12 Fv-responsive regulatory modules from 4049 differentially expressed genes (DEGs) that were activated or repressed by Fv infection in these seven cell types. Using a machining-learning approach, we constructed six cell type-specific immune regulatory networks by integrating Fv-induced DEGs from the cell type-specific transcriptomes, 16 known maize disease-resistant genes, five experimentally validated genes (ZmWOX5b, ZmPIN1a, ZmPAL6, ZmCCoAOMT2, and ZmCOMT), and 42 QTL or QTN predicted genes that are associated with Fv resistance. Taken together, this study provides not only a global view of maize cell fate determination during root development but also insights into the immune regulatory networks in major cell types of maize root tips at single-cell resolution, thus laying the foundation for dissecting molecular mechanisms underlying disease resistance in maize.

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.

The maturation and aging trajectory of Marchantia polymorpha at single-cell resolution

Bryophytes represent a sister to the rest of land plants. Despite their evolutionary importance and relatively simple body plan, a comprehensive understanding of the cell types and transcriptional states that underpin the temporal development of bryophytes has not been achieved. Using time-resolved single-cell RNA sequencing, we define the cellular taxonomy of Marchantia polymorpha across asexual reproduction phases. We identify two maturation and aging trajectories of the main plant body of M. polymorpha at single-cell resolution: the gradual maturation of tissues and organs along the tip-to-base axis of the midvein and the progressive decline of meristem activities in the tip along the chronological axis. Specifically, we observe that the latter aging axis is temporally correlated with the formation of clonal propagules, suggesting an ancient strategy to optimize allocation of resources to producing offspring. Our work thus provides insights into the cellular heterogeneity that underpins the temporal development and aging of bryophytes.

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.

Single-cell transcriptomic analyses reveal cellular and molecular patterns of rubber tree response to early powdery mildew infection

As a perennial woody plant, the rubber tree (Hevea brasiliensis) must adapt to various environmental challenges through gene expression in multiple cell types. It is still unclear how genes in this species are expressed at the cellular level and the precise mechanisms by which cells respond transcriptionally to environmental stimuli, especially in the case of pathogen infection. Here, we characterized the transcriptomes in Hevea leaves during early powdery mildew infection using single-cell RNA sequencing. We identified 10 cell types and constructed the first single-cell atlas of Hevea leaves. Distinct gene expression patterns of the cell clusters were observed under powdery mildew infection, which was especially significant in the epidermal cells. Most of the genes involved in host-pathogen interactions in epidermal cells exhibited a pattern of dramatically increased expression with increasing pseudotime. Interestingly, we found that the HbCNL2 gene, encoding a nucleotide-binding leucine-rich repeat protein, positively modulated the defence of rubber leaves against powdery mildew. Overexpression of the HbCNL2 gene triggered a typical cell death phenotype in tobacco leaves and a higher level of reactive oxygen species in the protoplasts of Hevea leaves. The HbCNL2 protein was located in the cytomembrane and nucleus, and its leucine-rich repeat domain interacted with the histidine kinase-like ATPase domain of the molecular chaperone HbHSP90 in the nucleus. Collectively, our results provide the first observation of the cellular and molecular responses of Hevea leaves to biotrophic pathogen infection and can guide the identification of disease-resistance genes in this important tree species.

CsBPC2 is a key regulator of root growth and development

BASIC PENTACYSTEINE (BPCs) transcription factors are important regulators of plant growth and development. However, the regulatory mechanism of BPC2 in roots remains unclear. In our previous study, we created Csbpc2 cucumber mutants by the CRISPR/Cas9 system, and our studies on the phenotype of Csbpc2 mutants showed that the root growth was inhibited compared with wide-type (WT). Moreover, the surface area, volume and number of roots decreased significantly, with root system architecture changing from dichotomous branching to herringbone branching. Compared with WT, the leaf growth of the Csbpc2 mutants was not affected. However, the palisade and spongy tissue were significantly thinner, which was not beneficial for photosynthesis. The metabolome of root exudates showed that compared with WT, amino acids and their derivatives were significantly decreased, and the enriched pathways were mainly regulated by amino acids and their derivatives, indicating that knockout of CsBPC2 mainly affected the amino acid content in root exudates. Importantly, transcriptome analysis showed that knockout of CsBPC2 mainly affected root gene expression. Knockout of CsBPC2 significantly reduced the gene expression of gibberellins synthesis. However, the expression of genes related to amino acid synthesis, nitrogen fixation and PSII-related photosynthesis increased significantly, which may be due to the effect of knocking out CsBPC2 on gibberellins synthesis, resulting in the inhibition of seedling growth, thus forming negative feedback regulation. Generally, we showed for the first time that BPC2 is a key regulator gene of root growth and development, laying the foundation for future mechanisms of BPC2 regulation in roots.

Plant BBR/BPC transcription factors: unlocking multilayered regulation in development, stress and immunity

MAIN CONCLUSION

This review provides a detailed structural and functional understanding of BBR/BPC TF, their conservation across the plant lineage, and their comparative study with animal GAFs. Plant-specific Barley B Recombinant/Basic PentaCysteine (BBR/BPC) transcription factor (TF) family binds to "GA" repeats similar to animal GAGA Factors (GAFs). These GAGA binding proteins are among the few TFs that regulate the genes at multiple steps by modulating the chromatin structure. The hallmark of the BBR/BPC TF family is the presence of a conserved C-terminal region with five cysteine residues. In this review, we present: first, the structural distinct yet functional similar relation of plant BBR/BPC TF with animal GAFs, second, the conservation of BBR/BPC across the plant lineage, third, their role in planta, fourth, their potential interacting partners and structural insights. We conclude that BBR/BPC TFs have multifaceted roles in plants. Besides the earliest identified function in homeotic gene regulation and developmental processes, presently BBR/BPC TFs were identified in hormone signaling, stress, circadian oscillation, and sex determination processes. Understanding how plants' development and stress processes are coordinated is central to divulging the growth-immunity trade-off regulation. The BBR/BPC TFs may hold keys to divulge the interactions between development and immunity. Moreover, the conservation of BBR/BPC across plant lineage makes it an evolutionary vital gene family. Consequently, BBR/BPCs are prospective to attract the increasing attention of the scientific communities as they are probably at the crossroads of diverse fundamental processes.

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 and spatial multi-omics in the plant sciences: Technical advances, applications, and perspectives

Plants contain a large number of cell types and exhibit complex regulatory mechanisms. Studies at the single-cell level have gradually become more common in plant science. Single-cell transcriptomics, spatial transcriptomics, and spatial metabolomics techniques have been combined to analyze plant development. These techniques have been used to study the transcriptomes and metabolomes of plant tissues at the single-cell level, enabling the systematic investigation of gene expression and metabolism in specific tissues and cell types during defined developmental stages. In this review, we present an overview of significant breakthroughs in spatial multi-omics in plants, and we discuss how these approaches may soon play essential roles in plant research.

An integrase toolbox to record gene-expression during plant development

There are many open questions about the mechanisms that coordinate the dynamic, multicellular behaviors required for organogenesis. Synthetic circuits that can record in vivo signaling networks have been critical in elucidating animal development. Here, we report on the transfer of this technology to plants using orthogonal serine integrases to mediate site-specific and irreversible DNA recombination visualized by switching between fluorescent reporters. When combined with promoters expressed during lateral root initiation, integrases amplify reporter signal and permanently mark all descendants. In addition, we present a suite of methods to tune the threshold for integrase switching, including: RNA/protein degradation tags, a nuclear localization signal, and a split-intein system. These tools improve the robustness of integrase-mediated switching with different promoters and the stability of switching behavior over multiple generations. Although each promoter requires tuning for optimal performance, this integrase toolbox can be used to build history-dependent circuits to decode the order of expression during organogenesis in many contexts.

Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of pigment glands in Gossypium bickii

Comprehensive utilization of cottonseeds is limited by the presence of pigment glands and its inclusion gossypol. The ideal cotton has glandless seeds but a glanded plant, a trait found in only a few Australian wild cotton species, including Gossypium bickii. Introgression of this trait into cultivated species has proved to be difficult. Understanding the biological processes toward pigment gland morphogenesis and the associated underlying molecular mechanisms will facilitate breeding of cultivated cotton varieties with the trait of glandless seeds and glanded plant. In this study, single-cell RNA sequencing (scRNA-seq) was performed on 12 222 protoplasts isolated from cotyledons of germinating G. bickii seeds 48 h after imbibition. Clustered into 14 distinct clusters unsupervisedly, these cells could be grouped into eight cell populations with the assistance of known cell marker genes. The pigment gland cells were well separated from others and could be separated into pigment gland parenchyma cells, secretory cells, and apoptotic cells. By integrating the pigment gland cell developmental trajectory, transcription factor regulatory networks, and core transcription factor functional validation, we established a model for pigment gland formation. In this model, light and gibberellin were verified to promote the formation of pigment glands. In addition, three novel genes, GbiERF114 (ETHYLENE RESPONSE FACTOR 114), GbiZAT11 (ZINC FINGER OF ARABIDOPSIS THALIANA 11), and GbiNTL9 (NAC TRANSCRIPTION FACTOR-LIKE 9), were found to affect pigment gland formation. Collectively, these findings provide new insights into pigment gland morphogenesis and lay the cornerstone for future cotton scRNA-seq investigations.

Single-nucleus sequencing deciphers developmental trajectories in rice pistils

Angiosperms possess a life cycle with an alternation of sporophyte and gametophyte generations, which happens in plant organs like pistils. Rice pistils contain ovules and receive pollen for successful fertilization to produce grains. The cellular expression profile in rice pistils is largely unknown. Here, we show a cell census of rice pistils before fertilization through the use of droplet-based single-nucleus RNA sequencing. The ab initio marker identification validated by in situ hybridization assists with cell-type annotation, revealing cell heterogeneity between ovule- and carpel-originated cells. A comparison of 1N (gametophyte) and 2N (sporophyte) nuclei identifies the developmental path of germ cells in ovules with typical resetting of pluripotency before the sporophyte-gametophyte transition, while trajectory analysis of carpel-originated cells suggests previously neglected features of epidermis specification and style function. These findings gain a systems-level view of cellular differentiation and development of rice pistils before flowering and lay a foundation for understanding female reproductive development in plants.

Improved cotton yield: Can we achieve this goal by regulating the coordination of source and sink?

Cotton is one of the major cash crops globally. It is characterized by determinate growth and multiple fruiting, which makes the source-sink contradiction more obvious. Coordination between source and sink is crucial for normal growth, yield, and quality of cotton. Numerous studies reported how the assimilate transport and distribution under varying environmental cues affected crop yields. However, less is known about the functional mechanism underlying the assimilate transport between source and sink, and how their distribution impacts cotton growth. Here, we provided an overview of the assimilate transport and distribution mechanisms , and discussed the regulatory mechanisms involved in source-sink balance in relation to cotton yield. Therefore, this review enriched our knowledge of the regulatory mechanism involved in source-sink relationship for improved cotton yield.

Cell-specific clock-controlled gene expression program regulates rhythmic fiber cell growth in cotton

BACKGROUND

The epidermis of cotton ovule produces fibers, the most important natural cellulose source for the global textile industry. However, the molecular mechanism of fiber cell growth is still poorly understood.

RESULTS

Here, we develop an optimized protoplasting method, and integrate single-cell RNA sequencing (scRNA-seq) and single-cell ATAC sequencing (scATAC-seq) to systematically characterize the cells of the outer integument of ovules from wild type and fuzzless/lintless (fl) cotton (Gossypium hirsutum). By jointly analyzing the scRNA-seq data from wildtype and fl, we identify five cell populations including the fiber cell type and construct the development trajectory for fiber lineage cells. Interestingly, by time-course diurnal transcriptomic analysis, we demonstrate that the primary growth of fiber cells is a highly regulated circadian rhythmic process. Moreover, we identify a small peptide GhRALF1 that circadian rhythmically controls fiber growth possibly through oscillating auxin signaling and proton pump activity in the plasma membrane. Combining with scATAC-seq, we further identify two cardinal cis-regulatory elements (CREs, TCP motif, and TCP-like motif) which are bound by the trans factors GhTCP14s to modulate the circadian rhythmic metabolism of mitochondria and protein translation through regulating approximately one third of genes that are highly expressed in fiber cells.

CONCLUSIONS

We uncover a fiber-specific circadian clock-controlled gene expression program in regulating fiber growth. This study unprecedentedly reveals a new route to improve fiber traits by engineering the circadian clock of fiber cells.

Single-cell analysis opens a goldmine for plant functional studies

Functional studies in biology require the identification of genes, regulatory elements, and networks, followed by a deep understanding of how they orchestrate to specify cell types, mediate signaling, and respond to internal and external cues over evolutionary timescales. Advances in single-cell analysis have enabled biologists to tackle these questions at the resolution of the individual cell. Here, we highlight recent studies in plants that have embraced single-cell analyses to facilitate functional studies. This review will provide guidance and perspectives for incorporating these advanced approaches in plant research for the coming decades.