Single cell RNA sequencing can yield high-resolution cell-type-specific expression signatures that reveal new cell types and the developmental trajectories of cell lineages. Here, we apply this approach to Arabidopsis (Arabidopsis thaliana) root cells to capture gene expression in 3,121 root cells. We analyze these data with Monocle 3, which orders single cell transcriptomes in an unsupervised manner and uses machine learning to reconstruct single cell developmental trajectories along pseudotime. We identify hundreds of genes with cell-type-specific expression, with pseudotime analysis of several cell lineages revealing both known and novel genes that are expressed along a developmental trajectory. We identify transcription factor motifs that are enriched in early and late cells, together with the corresponding candidate transcription factors that likely drive the observed expression patterns. We assess and interpret changes in total RNA expression along developmental trajectories and show that trajectory branch points mark developmental decisions. Finally, by applying heat stress to whole seedlings, we address the longstanding question of possible heterogeneity among cell types in the response to an abiotic stress. Although the response of canonical heat-shock genes dominates expression across cell types, subtle but significant differences in other genes can be detected among cell types. Taken together, our results demonstrate that single cell transcriptomics holds promise for studying plant development and plant physiology with unprecedented resolution. RESULTS Single-Cell RNA-Seq of Whole A. thaliana Roots Reveals Distinct Populations of Cortex, Endodermis, Hair, Nonhair, and Stele CellsWe used whole Arabidopsis roots from 7d-old seedlings to generate protoplasts for transcriptome analysis using the 103

The scarcity of accessible sites that are dynamic or cell type-specific in plants may be due in part to tissue heterogeneity in bulk studies. To assess the effects of tissue heterogeneity, we apply single-cell ATAC-seq to Arabidopsis thaliana roots and identify thousands of differentially accessible sites, sufficient to resolve all major cell types of the root. We find that the entirety of a cell’s regulatory landscape and its transcriptome independently capture cell type identity. We leverage this shared information on cell identity to integrate accessibility and transcriptome data to characterize developmental progression, endoreduplication and cell division. We further use the combined data to characterize cell type-specific motif enrichments of transcription factor families and link the expression of family members to changing accessibility at specific loci, resolving direct and indirect effects that shape expression. Our approach provides an analytical framework to infer the gene regulatory networks that execute plant development.

Vernalization is the promotion of the competence for flowering by long periods of low temperatures such as those typically experienced during winters. In Arabidopsis, the vernalization response is, to a large extent, mediated by the repression of the floral repressor FLC, and the stable epigenetic silencing of FLC after cold treatments is essential for vernalization. In addition to FLC, other vernalization targets exist in Arabidopsis. In grasses, vernalization seems to be entirely independent of FLC. Here, the current understanding of FLC-independent branches of the vernalization pathway in Arabidopsis and vernalization without FLC in grasses is discussed. This review focuses on the role of AGL19, AGL24, and the MAF genes in Arabidopsis. Interestingly, vernalization acts through related molecular machineries on distinct targets. In particular, protein complexes similar to Drosophila Polycomb Repressive Complex 2 play a prominent role in establishing an epigenetic cellular memory for cold-regulated expression states of AGL19 and FLC. Finally, the similar network topology of the apparently independently evolved vernalization pathways of grasses and Arabidopsis is discussed.

Arabidopsis MSI1 has fundamental functions in plant development. MSI1 is a subunit of Polycomb group protein complexes and Chromatin assembly factor 1, and it interacts with the Retinoblastoma-related protein 1. Altered levels of MSI1 result in pleiotropic phenotypes, reflecting the complexity of MSI1 protein functions. In order to uncover additional functions of MSI1, we performed transcriptional profiling of wild-type and plants with highly reduced MSI1 levels (msi1-cs). Surprisingly, the known functions of MSI1 could only account for a minor part of the transcriptional changes in msi1-cs plants. One of the most striking unexpected observations was the up-regulation of a subset of ABA-responsive genes eliciting the response to drought and salt stress. We report that MSI1 can bind to the chromatin of the drought-inducible downstream target RD20 and suggest a new role for MSI1 in the negative regulation of the Arabidopsis drought-stress response.