Arbuscular mycorrhizal fungi (AMF) establish a dynamic and asynchronous symbiosis with a wide range of land plants, which involves distinct stages of root colonization and associated cellular responses that co-occur within the same root. While decades of research have significantly advanced our understanding of the plant's symbiotic gene repertoire, this spatial and temporal complexity has hindered a detailed dissection of the molecular mechanisms underlying fungal accommodation. Here, we present the first single-nucleus RNA-sequencing (snRNA-seq) dataset of Solanum lycopersicum roots colonized by Rhizophagus irregularis. Unsupervised subclustering of an arbuscular mycorrhiza (AM)-specific cell population resolves AM-responsive root epidermal cells as well as a developmental gradient of cortical cells across distinct stages of arbuscule formation, thus unveiling stage-specific transcriptional signatures during AMF colonization. Moreover, using motif-informed network inference based on single-cell expression data (MINI-EX), we put forward candidate transcription factors orchestrating these stage-specific transcriptional programs. Together, our data support novel hypotheses on how diverse plant developmental and physiological processes-including localized cell-cycle reactivation and the integration of multiple nutritional cues-are coordinated to facilitate the establishment of a functional symbiosis. As such, this high-resolution dataset serves as a valuable resource for candidate gene prioritization and future reverse genetic studies.

The isolation of single plant cells from complex tissues is prone to selective enrichment and sampling biases, which complicates accurate profiling of the large diversity in cell types. Optimizing methodologies for cell enrichment and single-cell transcriptomics is therefore critical for single-cell studies addressing plant cell heterogeneity. Here, we systematically compared protoplast enrichment technologies (including conventional and image-based flow cytometry, as well as magnetic cell sorting) and single-cell RNA sequencing (scRNA-seq) platforms (10X Genomics Chromium, BD Rhapsody) using Arabidopsis roots. Image-based flow cytometry offered increased precision due to customizable gating strategies, while magnetic sorting provided faster processing and enhanced representation of cell size heterogeneity. Both scRNA-seq platforms captured root cell heterogeneity and yielded reproducible gene expression profiles, but showed platform-associated differences in cell type composition. Notably, single-nucleotide polymorphism analysis of a mixed ecotype sample revealed that, among cells identified as doublets by computational algorithms, two-thirds were likely to have been misclassified. These insights identify key biases in plant cell purification and scRNA-seq workflows and provide practical guidance for improving data quality across plant species.

The endosperm is a transient nutritive tissue in plant seeds. During maize (Zea mays) grain development, 2 distinct endosperm cell death processes occur: in 1 process, the endosperm adjacent to the embryo scutellum (EAS) is completely dismantled; in the other, the starchy endosperm (SE) retains nutrient-packed cell corpses after grain filling. Here, we show that SE cell death degrades some organelles including the mitochondria and the endoplasmic reticulum, while preserving protein bodies, starch granules, and chromatin. In contrast, EAS cells undergo lytic cell death to remobilize stored nutrients through a complex corpse clearance process. Using single-cell transcriptome analysis, we identified 2 NAC transcription factors, KIRA-LIKE 1 (KIL1) and 2 (KIL2), as specifically upregulated in the EAS. Analyses using dominant and recessive loss-of-function kil mutants demonstrate that these genes redundantly promote cell death and corpse clearance in the EAS, but are not required for SE cell death. Reduced EAS cell death in kil loss-of-function mutants strongly impeded embryo growth, indicating that EAS elimination is crucial for optimal embryo development. Notably, kil1 and kil2 expression is regulated by DOSAGE-EFFECT DEFECTIVE1, an imprinted paternally expressed endosperm transcription factor. Our findings suggest paternal control over EAS cell death and the embryo-endosperm size ratio in maize, providing new leads to modulate this agronomically important trait.

During secondary growth, the vascular cambium produces conductive xylem and phloem cells, while the phellogen (cork cambium) deposits phellem (cork) as the outermost protective barrier. Although most of the secondary tissues are made up of parenchyma cells, which are also produced by both cambia, their diversity and function are poorly understood. Here we combined single-cell RNA sequencing analysis with lineage tracing to recreate developmental trajectories of the cell types in the Arabidopsis root undergoing secondary growth. By analysing 93 reporter lines, we were able to identify 20 different cell types or cell states, many of which have not been described before. We additionally observed distinct transcriptome signatures of parenchyma cells depending on their maturation state and proximity to the conductive cell types. Our data show that both xylem and phloem parenchyma tissues are required for normal formation of conductive tissue cell types. Furthermore, we show that mature phloem parenchyma gradually obtains periderm identity, and this transformation can be accelerated by jasmonate treatment or wounding. Our study thus reveals the diversity of parenchyma cells and their capacity to undergo considerable identity changes during secondary growth.

Somatic embryogenesis allows a somatic plant cell to develop into an embryo, and potentially into a fertile plant. Transcription factors such as BABY BOOM can induce somatic embryogenesis when ectopically expressed and are widely used for aiding regeneration in tissue culture for transformation and gene editing of crops. Nevertheless, regeneration remains a bottleneck and alternative morphogenic genes are highly desired. Here, we co-expressed BABY BOOM and WUSCHEL2 in zygotic maize ( Zea may s L.) embryos and studied gene regulatory networks in induced somatic embryos at the single-cell level. By inferring cell-type-specific regulons, we prioritized candidate regulators and confirmed functionality of four transcription factors, bHLH48, EREB152, GRF4, and HB77, for enhanced maize transformation frequency, leading to fertile, transgenic plants. Interestingly, the basic helix-loop-helix and homeodomain-leucine zipper families had previously not been associated with induced somatic embryogenesis. Our work will contribute to more efficient transformation, much needed to deliver on the promise of gene editing for agriculture. ### Competing Interest Statement The authors (J.R., S.L., K.V., and L.P.) declare that a patent application related to the content of this manuscript has been filed. Industrieel Onderzoeksfonds, F2020/IOF-StarTT/151; IOF.PRO.2021.0017

Drought stress imposes severe challenges on agriculture by impacting crop performance. Understanding drought responses in plants at a cellular level is a crucial first step toward engineering improved drought resilience. However, the molecular responses to drought are complex as they depend on multiple factors, including the severity of drought, the profiled organ, its developmental stage or even the cell types therein. Thus, deciphering the transcriptional responses to drought is especially challenging. In this study, we investigated tissue-specific responses to mild drought (MD) in young Arabidopsis thaliana (Arabidopsis) leaves using single-cell RNA sequencing (scRNA-seq). To preserve transcriptional integrity during cell isolation, we inhibited RNA synthesis using the transcription inhibitor actinomycin D, and demonstrated the benefits of transcriptome fixation for studying mild stress responses at a single-cell level. We present a curated and validated single-cell atlas, comprising 50 797 high-quality cells from almost all known cell types present in the leaf. All cell type annotations were validated with a new library of reporter lines. The curated data are available to the broad community in an intuitive tool and a browsable single-cell atlas (http://www.single-cell.be/plant/leaf-drought). We show that the mesophyll contains two spatially separated cell populations with distinct responses to drought: one enriched in canonical abscisic acid-related drought-responsive genes, and another one enriched in genes involved in iron starvation responses. Our study thus reveals a dual adaptive mechanism of the leaf mesophyll in response to MD and provides a valuable resource for future research on stress responses.

Plant-parasitic nematodes are responsible for important annual losses in crop productivity worldwide[1][1],[2][2]. Although the formation of feeding organs within the roots is essential for successful sedentary parasitism[3][3], the molecular mechanisms underlying their development are poorly understood. This is partly because these organs originate from a limited number of root cells[4][4]–[7][5], making difficult to capture the transcriptional reprogramming that occurs during the early stages of the infection. Here, we first developed a comparative host-pathogen framework to study the nematode infection process in Arabidopsis and rice. Using a cross-species single-cell transcriptomics approach, we identified a unique molecular signature in infected root cells and show that the cellular reprogramming during these early stages is highly conserved across both host-pathogen interactions. This transcriptional cell reprogramming is associated with stemness acquisition related to de novo organogenesis process. By cell-type specific gene regulatory network analysis, we identified AtATHB2/OsHOX28 as an evolutionary conserved and key regulator of the nematode infection process. Loss-of-function of this regulator in both species results in nematode resistance without affecting root growth. This discovery opens up new avenues for the development of sustainable nematode control strategies that could be translated across crop species. ### Competing Interest Statement The authors have declared no competing interest. [1]: #ref-1 [2]: #ref-2 [3]: #ref-3 [4]: #ref-4 [5]: #ref-7

Despite the broad use of single-cell/nucleus RNA sequencing in plant research, accurate cluster annotation in less-studied plant species remains a major challenge due to the lack of validated marker genes. Here, we generated a single-cell RNA sequencing atlas of soil-grown wheat roots and annotated cluster identities by transferring annotations from publicly available datasets in wheat, rice, maize, and Arabidopsis. The predictions from our orthology-based annotation approach were next validated using untargeted spatial transcriptomics. These results allowed us to predict evolutionarily conserved tissue-specific markers and generate cell type-specific gene regulatory networks for root tissues of wheat and the other species used in our analysis. In summary, we generated a single-cell and spatial transcriptomics resource for wheat root apical meristems, including numerous known and uncharacterized cell type-specific marker genes and developmental regulators. These data and analyses will facilitate future cell type annotation in non-model plant species.

MOTIVATION:Existing nanopore single-cell data analysis tools showed severe limitations in handling current data sizes. RESULTS:We introduce scywalker, an innovative and scalable package developed to comprehensively analyze long-read sequencing data of full-length single-cell or single-nuclei cDNA. We developed novel scalable methods for cell barcode demultiplexing and single-cell isoform calling and quantification and incorporated these in an easily deployable package. Scywalker streamlines the entire analysis process, from sequenced fragments in FASTQ format to demultiplexed pseudobulk isoform counts, into a single command suitable for execution on either server or cluster. Scywalker includes data quality control, cell type identification, and an interactive report. Assessment of datasets from the human brain, Arabidopsis leaves, and previously benchmarked data from mixed cell lines demonstrate excellent correlation with short-read analyses at both the cell-barcoding and gene quantification levels. At the isoform level, we show that scywalker facilitates the direct identification of cell-type-specific expression of novel isoforms. AVAILABILITY AND IMPLEMENTATION:Scywalker is available on github.com/derijkp/scywalker under the GNU General Public License (GPL) and at https://zenodo.org/records/13359438/files/scywalker-0.108.0-Linux-x86_64.tar.gz.

Aluminum-dependent stoppage of root growth requires the DNA damage response (DDR) pathway including the p53-like transcription factor SUPPRESSOR OF GAMMA RADIATION 1 (SOG1), which promotes terminal differentiation of the root tip in response to Al dependent cell death. Transcriptomic analyses identified Al-induced SOG1-regulated targets as candidate mediators of this growth arrest. Analysis of these factors either as loss-of-function mutants or by overexpression in the als3-1 background shows ERF115, which is a key transcription factor that in other scenarios is rate-limiting for damaged stem cell replenishment, instead participates in transition from an actively growing root to one that has terminally differentiated in response to Al toxicity. This is supported by a loss-of-function erf115 mutant raising the threshold of Al required to promote terminal differentiation of Al hypersensitive als3-1. Consistent with its key role in stoppage of root growth, a putative ERF115 barley ortholog is also upregulated following Al exposure, suggesting a conserved role for this ATR-dependent pathway in Al response. In contrast to other DNA damage agents, these results show that ERF115 and likely related family members are important determinants of terminal differentiation of the root tip following Al exposure and central outputs of the SOG1-mediated pathway in Al response.

Non-small cell lung cancer (NSCLC) is known for high relapse rates despite resection in early stages. present the results of a phase I clinical trial in which a dendritic cell (DC) vaccine targeting patient neoantigens is evaluated in patients with resected NSCLC. Vaccine manufacturing is feasible enrolled patients. Toxicity is limited to grade 1-2 adverse events. Systemic T cell responses in five out of six vaccinated patients, with T cell responses remaining detectable up to 19 months nation. Single-cell analysis indicates that the responsive T cell population is polyclonal and exhibits entire spectrum of T cell differentiation states, including a naive-like state, but excluding exhausted states. Three of six vaccinated patients experience disease recurrence during the follow-up period Collectively, these data support the feasibility, safety, and immunogenicity of this treatment NSCLC.

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

Droplet-based single-cell RNA-sequencing (scRNA-seq) empowers transcriptomic profiling with an unprecedented resolution, facilitating insights into the cellular heterogeneity of tissues, developmental progressions, stress-response dynamics, and more at single-cell level. In this chapter, we describe the experimental workflow of processing Arabidopsis root tissue into protoplasts and generating single-cell transcriptomes. We also describe the general computational workflow of visualizing and utilizing scRNA-seq data. This protocol can be used as a starting point for establishing a scRNA-seq workflow.

Survival of living organisms is fully dependent on their maintenance of genome integrity, being permanently threatened by replication stress in proliferating cells. Although the plant DNA damage response (DDR) regulator SOG1 has been demonstrated to cope with replication defects, accumulating evidence points to other pathways functioning independent of SOG1. Here, we report the roles of the Arabidopsis E2FA and EF2B transcription factors, two well-characterized regulators of DNA replication, in plant response to replication stress. Through a combination of reverse genetics and chromatin immunoprecipitation approaches, we show that E2FA and E2FB share many target genes with SOG1, providing evidence for their involvement in the DDR. Analysis of double-and triple-mutant combinations revealed that E2FB, rather than E2FA, plays the most prominent role in sustaining plant growth in the presence of replication defects, either operating antagonistically or synergistically with SOG1. Conversely, SOG1 aids in overcoming the replication defects of E2FA/E2FB-deficient plants. Collectively, our data reveal a complex transcriptional network controlling the replication stress response in which E2Fs and SOG1 act as key regulatory factors.

Plant roots possess remarkable regenerative potential owing to their ability to replenish damaged or lost stem cells. ETHYLENE RESPONSE FACTOR 115 (ERF115), one of the key molecular elements linked to this potential, plays a predominant role in the activation of regenerative cell divisions. However, the downstream operating molecular machinery driving wound-activated cell division is largely unknown. Here, we biochemically and genetically identified the GRAS-domain transcription factor SCARECROW-LIKE 5 (SCL5) as an interaction partner of ERF115 in Arabidopsis thaliana. Although non-essential under control growth conditions, SCL5 acts redundantly with the related PHYTOCHROME A SIGNAL TRANSDUCTION 1 (PAT1) and SCL21 transcription factors to activate the expression of the DNA-BINDING ONE FINGER 3.4 (DOF3.4) transcription factor gene. DOF3.4 expression is wound-inducible in an ERF115-dependent manner, and in turn activates D3-type cyclin expression. Accordingly, ectopic DOF3.4 expression drives periclinal cell division, while its downstream D3-type cyclins are essential for the regeneration of a damaged root. Our data highlight the importance and redundant roles of the SCL5, SCL21 and PAT1 transcription factors in wound-activated regeneration processes and pinpoint DOF3.4 as a key downstream element driving regenerative cell division.

Plant roots possess a remarkable regenerative potential owing to their ability to replenish damaged or lost stem cells. The ETHYLENE RESPONSE FACTOR 115 (ERF115), one of the key molecular elements that have been linked to this potential, plays a predominant role in the activation of regenerative cell divisions. However, the downstream operating molecular machinery driving wound-activated cell division is largely unknown. Here, we biochemically and genetically identified the GRAS-domain transcription factor SCARECROW-LIKE 5 (SCL5) as an interaction partner of ERF115. Although being non-essential under control growth conditions, SCL5 controls redundantly with the related PHYTOCHROME A SIGNAL TRANSDUCTION 1 (PAT1) and SCL21 transcription factors the expression of the DNA-BINDING ONE FINGER 3.4 ( DOF3.4 ) transcription factor gene. DOF3.4 expression is wound-inducible in an ERF115-dependent manner, and in turn activates D3-type cyclin expression. Accordingly, ectopic DOF3.4 expression drives periclinal cell division, while its downstream D3-type cyclins are essential for regeneration of a damaged root. Our data highlight the importance and redundant role of the SCL5, SCL21 and PAT1 transcription factors in the process of wound-activated regeneration processes and pinpoint the DOF3.4 as a key downstream element driving regenerative cell division.

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

The anaphase-promoting complex/cyclosome (APC/C) marks key cell cycle proteins for proteasomal breakdown, thereby ensuring unidirectional progression through the cell cycle. Its target recognition is temporally regulated by activating subunits, one of which is called CELL CYCLE SWITCH 52 A2 (CCS52A2). We sought to expand the knowledge on the APC/C by using the severe growth phenotypes of CCS52A2 -deficient Arabidopsis thaliana plants as a readout in a suppressor mutagenesis screen, resulting in the identification of the previously undescribed gene called PIKMIN1 ( PKN1 ). PKN1 deficiency rescues the disorganized root stem cell phenotype of the ccs52a2-1 mutant, whereas an excess of PKN1 inhibits growth of ccs52a2-1 plants, indicating the need for control of PKN1 abundance for proper development. Accordingly, the lack of PKN1 in a wild-type background negatively impacts cell division, while its systemic overexpression promotes proliferation. PKN1 shows a cell cycle phase-dependent accumulation pattern, localizing to microtubular structures, including the preprophase band, the mitotic spindle, and phragmoplast. PKN1 is conserved throughout the plant kingdom, with its function in cell division being evolutionary conserved in the liverwort Marchantia polymorpha . Our data thus demonstrate that PKN1 represents a novel, plant-specific gene with a role in cell division, which is likely proteolytically controlled by the CCS52A2-activated APC/C.

The anaphase-promoting complex/cyclosome (APC/C) marks key cell cycle proteins for proteasomal breakdown, thereby ensuring unidirectional progression through the cell cycle. Its target recognition is temporally regulated by activating subunits, one of which is called CELL CYCLE SWITCH 52 A2 (CCS52A2). We sought to expand the knowledge on the APC/C by using the severe growth phenotypes of CCS52A2-deficient Arabidopsis (Arabidopsis thaliana) plants as a readout in a suppressor mutagenesis screen, resulting in the identification of the previously undescribed gene called PIKMIN1 (PKN1). PKN1 deficiency rescues the disorganized root stem cell phenotype of the ccs52a2-1 mutant, whereas an excess of PKN1 inhibits the growth of ccs52a2-1 plants, indicating the need for control of PKN1 abundance for proper development. Accordingly, the lack of PKN1 in a wild-type background negatively impacts cell division, while its systemic overexpression promotes proliferation. PKN1 shows a cell cycle phase-dependent accumulation pattern, localizing to microtubular structures, including the preprophase band, the mitotic spindle, and the phragmoplast. PKN1 is conserved throughout the plant kingdom, with its function in cell division being evolutionarily conserved in the liverwort Marchantia polymorpha. Our data thus demonstrate that PKN1 represents a novel, plant-specific protein with a role in cell division that is likely proteolytically controlled by the CCS52A2-activated APC/C.

In Arabidopsis thaliana , the negative brassinosteroid (BR) signalling regulator, BR INSENSITIVE2 (BIN2) promotes and restricts stomatal asymmetric cell division (ACD) depending on its subcellular localization, which is regulated by the stomatal lineage-specific scaffolding protein POLAR. BRs inactivate BIN2, but how they govern stomatal development remains unclear. Mapping the single-cell transcriptome of stomatal lineages with either exogenous BRs or the specific BIN2 inhibitor revealed that the two modes of BR activation triggered spatiotemporally distinct transcriptional responses. We established that when in a complex with POLAR and its closest homolog POLAR-LIKE1, BIN2 is insulated from BR-mediated inactivation, nevertheless, it remains sensitive to the inhibitor. Subsequently, BR signalling is attenuated in ACD precursors, whereas it remains active in epidermal cells that would undergo differentiation. Our study demonstrates how scaffold proteins contribute to cellular signal specificity of hormonal responses in plants.