Multiome in the Same Cell Reveals the Impact of Osmotic Stress on Arabidopsis Root Tip Development at Single‐Cell Level

Abstract

Cell‐specific transcriptional regulatory networks (TRNs) play vital roles in plant development and response to environmental stresses. However, traditional single‐cell mono‐omics techniques are unable to directly capture the relationships and dynamics between different layers of molecular information within the same cells. While advanced algorithm facilitates merging scRNA‐seq and scATAC‐seq datasets, accurate data integration remains a challenge, particularly when investigating cell‐type‐specific TRNs. By examining gene expression and chromatin accessibility simultaneously in 16,670 Arabidopsis root tip nuclei, the TRNs are reconstructed that govern root tip development under osmotic stress. In contrast to commonly used computational integration at cell‐type level, 12,968 peak‐to‐gene linkage is captured at the bona fide single‐cell level and construct TRNs at an unprecedented resolution. Furthermore, the unprecedented datasets allow to more accurately reconstruct the coordinated changes of gene expression and chromatin states during cellular state transition. During root tip development, chromatin accessibility of initial cells precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for subsequent differentiation steps. Pseudo‐time trajectory analysis reveal that osmotic stress can shift the functional differentiation of trichoblast. Candidate stress‐related gene‐linked cis‐regulatory elements (gl‐cCREs) as well as potential target genes are also identified, and uncovered large cellular heterogeneity under osmotic stress.