Leaf cell-specific and single-cell transcriptional profiling reveals a role for the palisade layer in UV light protection

Beauty is more than epidermis deep: How cell division and expansion sculpt the leaf spongy mesophyll

As the main location of photosynthesis, leaf mesophyll cells are one of the most abundant and essential cell types on earth. Forming the bulk of the internal tissues of the leaf, their size, shape, and patterns of interconnectivity define the internal structure and surface area of the leaf, which in turn determines the efficiency of light capture and carbon fixation. Understanding how these cellular traits are controlled and translated into tissue- and organ-scale traits, and how they influence photosynthetic performance will be key to our ability to improve crop plants in the face of a changing climate. In contrast to the extensive literature on the anatomical and physiological aspects of mesophyll function, our understanding of the cell-level morphogenetic processes underpinning mesophyll cell growth and differentiation is scant. In this review, we focus on how cell division, expansion, and separation are coordinated to create the intricate architecture of the spongy mesophyll.

Large language models in plant biology

Large language models (LLMs), such as ChatGPT, have taken the world by storm. However, LLMs are not limited to human language and can be used to analyze sequential data, such as DNA, protein, and gene expression. The resulting foundation models can be repurposed to identify the complex patterns within the data, resulting in powerful, multipurpose prediction tools able to predict the state of cellular systems. This review outlines the different types of LLMs and showcases their recent uses in biology. Since LLMs have not yet been embraced by the plant community, we also cover how these models can be deployed for the plant kingdom.

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.

Multiple factors and features dictate the selective production of ct-siRNA in Arabidopsis

Coding transcript-derived siRNAs (ct-siRNAs) produced from specific endogenous loci can suppress the translation of their source genes to balance plant growth and stress response. In this study, we generated Arabidopsis mutants with deficiencies in RNA decay and/or post-transcriptional gene silencing (PTGS) pathways and performed comparative sRNA-seq analysis, revealing that multiple RNA decay and PTGS factors impede the ct-siRNA selective production. Genes that produce ct-siRNAs often show increased or unchanged expression and typically have higher GC content in sequence composition. The growth and development of plants can perturb the dynamic accumulation of ct-siRNAs from different gene loci. Two nitrate reductase genes, NIA1 and NIA2, produce massive amounts of 22-nt ct-siRNAs and are highly expressed in a subtype of mesophyll cells where DCL2 exhibits higher expression relative to DCL4, suggesting a potential role of cell-specific expression of ct-siRNAs. Overall, our findings unveil the multifaceted factors and features involved in the selective production and regulation of ct-siRNAs and enrich our understanding of gene silencing process in plants.

Graphene oxide-based aerogel stimulates growth, mercury accumulation, photosynthesis-related gene expression, antioxidant efficiency and redox status in wheat under mercury exposure

Mercury (Hg) pollution is a global concern in cropland systems. Hg contamination causes a disruption in the growth, energy metabolism, redox balance, and photosynthetic activity of plants. In the removal of Hg toxicity, a recent critical strategy is the use of aerogels with biodegradability and biocompatibility. However, it is unknown how graphene oxide-based aerogels stimulate the defense systems in wheat plants exposed to Hg toxicity. Therefore, in this study, the photosynthetic, genetic, and biochemical effects of reduced graphene oxide aerogel treatments (gA; 50-100-250 mg L-1) were examined in wheat (Triticum aestivum) under Hg stress (50 μM HgCl2). The relative growth rate (RGR) significantly decreased (84%) in response to Hg stress. However, the reduced RGR and water relations (RWC) of wheat were improved by gA treatments. The impaired gas exchange levels (stomatal conductance, carbon assimilation rate, intercellular CO2 concentrations, and transpiration rate) caused by stress were reversed under Hg plus gAs. Additionally, stress hampered chlorophyll fluorescence (Fv/Fo, Fv/Fm), and under Hg toxicity the expression of psaA genes was reduced (>0.4-fold), but psaB gene was significantly up-regulated (>3-fold) which are the genes involved in PSI. By increasing expression patterns of both genes relating to PSI, gAs reversed the adverse consequences on Fv/Fo and Fv/Fm in the presence of excessive Hg concentration. The activities of glutathione S-transferase (GST), glutathione reductase (GR), catalase (CAT), ascorbate peroxidase (APX) and dehydroascorbate reductase (DHAR) decreased under Hg toxicity. On the other hand, the activities of superoxide dismutase (SOD), APX, GST, and glutathione peroxidase (GPX) increased following gA treatments against stress, leading to the successful elimination of toxic levels of H2O2 and lipid peroxidation (TBARS content) by decreasing the levels by about 30%, and 40%, respectively. By modulating enzyme/non-enzyme activity/contents including the AsA-GSH cycle, gAs contributed to the protection of the cellular redox state. Most important of all, gA applications were able to reduce Hg intake by approximately 66%. Therefore, these results showed that gAs were effective in highly inhibiting Hg uptake and could significantly increase wheat tolerance to toxicity by eliminating Hg-induced oxidative damage and inhibiting metabolic processes involved in photosynthesis. The findings obtained from the study provide a new perspective on the alleviation roles of reduced graphene oxide aerogels as an effective adsorbent for decreasing damages of mercury toxicity in wheat plants.

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.

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

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

UV-B responses in the spotlight: Dynamic photoreceptor interplay and cell-type specificity

Plants are constantly exposed to a multitude of external signals, including light. The information contained within the full spectrum of light is perceived by a battery of photoreceptors, each with specific and shared signalling outputs. Recently, it has become clear that UV-B radiation is a vital component of the electromagnetic spectrum, guiding growth and being crucial for plant fitness. However, given the large overlap between UV-B specific signalling pathways and other photoreceptors, understanding how plants can distinguish UV-B specific signals from other light components deserves more scrutiny. With recent evidence, we propose that UV-B signalling and other light signalling pathways occur within distinct tissues and cell-types and that the contribution of each pathway depends on the type of response and the developmental stage of the plant. Elucidating the precise site(s) of action of each molecular player within these signalling pathways is key to fully understand how plants are able to orchestrate coordinated responses to light within the whole plant body. Focusing our efforts on the molecular study of light signal interactions to understand plant growth in natural environments in a cell-type specific manner will be a next step in the field of photobiology.

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.

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.

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

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

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.

Leaf hydraulic maze: Abscisic acid effects on bundle sheath, palisade, and spongy mesophyll conductance

Leaf hydraulic conductance (Kleaf) facilitates the supply of water, enabling continual CO2 uptake while maintaining plant water status. We hypothesized that bundle sheath and mesophyll cells play key roles in regulating the radial flow of water out of the xylem by responding to abscisic acid (ABA). Thus, we generated transgenic Arabidopsis thaliana plants that are insensitive to ABA in their bundle sheath (BSabi) and mesophyll (MCabi) cells. We also introduced tissue-specific fluorescent markers to distinguish between cells of the palisade mesophyll, spongy mesophyll, and bundle sheath. Both BSabi and MCabi plants showed greater Kleaf and transpiration under optimal conditions. MCabi plants had larger stomatal apertures, higher stomatal index, and greater vascular diameter and biomass relative to the wild-type (WT) and BSabi plants. In response to xylem-fed ABA, both transgenic and WT plants reduced their Kleaf and transpiration. The membrane osmotic water permeability (Pf) of the WT's spongy mesophyll was higher than that of the WT's palisade mesophyll. While the palisade mesophyll maintained a low Pf in response to high ABA, the spongy mesophyll Pf was reduced. Compared to the WT, BSabi bundle sheath cells had a higher Pf, but MCabi spongy mesophyll had an unexpected lower Pf. These results suggest that tissue-specific regulation of Pf by ABA may be confounded by whole-leaf hydraulics and transpiration. ABA increased the symplastic permeability, but its contribution to Kleaf was negligible. We suggest that the bundle sheath spongy mesophyll pathway dynamically responds to the fluctuations in water availability, while the palisade mesophyll serves as a hydraulic buffer.

Integrated mass spectrometry imaging and single-cell transcriptome atlas strategies provide novel insights into taxoid biosynthesis and transport in Taxus mairei stems

Taxol, which is a widely used important chemotherapeutic agent, was originally isolated from Taxus stem barks. However, little is known about the precise distribution of taxoids and the transcriptional regulation of taxoid biosynthesis across Taxus stems. Here, we used MALDI-IMS analysis to visualize the taxoid distribution across Taxus mairei stems and single-cell RNA sequencing to generate expression profiles. A single-cell T. mairei stem atlas was created, providing a spatial distribution pattern of Taxus stem cells. Cells were reordered using a main developmental pseudotime trajectory which provided temporal distribution patterns in Taxus stem cells. Most known taxol biosynthesis-related genes were primarily expressed in epidermal, endodermal, and xylem parenchyma cells, which caused an uneven taxoid distribution across T. mairei stems. We developed a single-cell strategy to screen novel transcription factors (TFs) involved in taxol biosynthesis regulation. Several TF genes, such as endodermal cell-specific MYB47 and xylem parenchyma cell-specific NAC2 and bHLH68, were implicated as potential regulators of taxol biosynthesis. Furthermore, an ATP-binding cassette family transporter gene, ABCG2, was proposed as a potential taxoid transporter candidate. In summary, we generated a single-cell Taxus stem metabolic atlas and identified molecular mechanisms underpinning the cell-specific transcriptional regulation of the taxol biosynthesis pathway.

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.

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.

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 RNA Sequencing of Arabidopsis Leaf Tissues Identifies Multiple Specialized Cell Types: Idioblast Myrosin Cells and Potential Glucosinolate-Producing Cells

The glucosinolate-myrosinase defense system (GMDS), characteristic of Brassicales, is involved in plant defense. Previous single-cell transcriptomic analyses have reported the expression profiles of multiple GMDS-related cell types (i.e. myrosinase-rich myrosin idioblasts and multiple types of potential glucosinolate synthetic cells as well as a candidate S-cell for glucosinolate accumulation). However, differences in plant stages and cell-type annotation methods have hindered comparisons among studies. Here, we used the single-cell transcriptome profiles of extended Arabidopsis leaves and verified the distribution of previously used markers to refine the expression profiles of GMDS-associated cell types. Moreover, we performed beta-glucuronidase promoter assays to confirm the histological expression patterns of newly obtained markers for GMDS-associated candidates. As a result, we found a set of new specific reporters for myrosin cells and potential glucosinolate-producing cells.

Water deficit aggravated the inhibition of photosynthetic performance of maize under mercury stress but is alleviated by brassinosteroids

Mercury (Hg) significantly inhibits maize (Zea mays L.) production, which could be aggravated by water deficit (WD) due to climate change. However, there is no report on the maize in response to combined their stresses. This work was conducted for assessing the response and adaptive mechanism of maize to combined Hg and WD stress using two maize cultivars, Xianyu (XY) 335 and Yudan (YD) 132. The analysis was based on plant growth, physiological function, and transcriptomic data. Compared with the single Hg stress, Hg accumulation in whole plant and translocation factor (TF) under Hg+WD were increased by 64.51 % (1.44 mg kg^-1) and 260.00 %, respectively, for XY 335; and 50.32 % (0.62 mg kg^-1) and 220.02 %, respectively, for YD 132. Combined Hg and WD stress further increased the reactive oxygen species accumulation, aggravated the damage of the thylakoid membrane, and decreased chlorophyll content compared with single stress. For example, Chl a and Chl b contents of XY 335 were significantly decreased by 48.67 % and 28.08 %, respectively at 48 h after Hg+WD treatment compared with Hg stress. Furthermore, transcriptome analysis revealed that most of down-regulated genes were enriched in photosynthetic-antenna proteins, photosynthesis, chlorophyll and porphyrin metabolism pathways (PsbS1, PSBQ1 and FDX1 etc.) under combined stress, reducing light energy capture and electron transport. However, most genes related to the brassinosteroids (BRs) signaling pathway were up-regulated under Hg+WD stress. Correspondingly, exogenous BRs significantly enhanced the maize tolerance to stress by decreasing Hg accumulation and TF, and raising activities of antioxidant enzyme, the content of chlorophyll and photosynthetic performance. The PI, Fv/Fm and Fv/Fo of Hg+WD+BR treatment were increased by 29.88 %, 32.06 %, and 14.56 %, respectively, for XY 335 compared to Hg+WD. Overall, combined Hg and WD stress decreased photosynthetic efficiency by adversely affecting light absorption and electron transport, especially in stress-sensitive variety, but BRs could alleviate the inhibition of photosynthesis, providing a novel strategy for enhancing crop Hg and WD tolerance and food safety.