Abstract
In plants, microRNAs (miRNAs) participate in complex gene regulatory networks together with the transcription factors (TFs) in response to biotic and abiotic stresses. To date, analyses of miRNAs-induced transcriptome remodeling are at the whole plant or tissue levels. Here, Arabidopsis's ABA-induced single-cell RNA-seq (scRNA-seq) is performed at different stages of time points-early, middle, and late. Single-cell level primary miRNAs (pri-miRNAs) atlas supported the rapid, dynamic, and cell-type specific miRNA responses under ABA treatment. MiRNAs respond rapidly and prior to target gene expression dynamics, and these rapid response miRNAs are highly cell-type specific, especially in mesophyll and vascular cells. MiRNA-TF-mRNA regulation modules are identified by identifying miRNA-contained feed-forward loops (M-FFLs) in the regulatory network, and regulatory networks with M-FFLs have higher co-expression and clustering coefficient (CC) values than those without M-FFLs, suggesting the hub role of miRNAs in regulatory networks. The cell-type-specific M-FFLs are regulated by these hub miRNAs rather than TFs through sc-RNA-seq network analysis. MiR858a-FBH3-MYB module inhibited the expression of MYB63 and MYB20, which related to the formation of plant secondary wall and the production of lignin, through M-FFL specifically in vascular. These results can provide prominent insights into miRNAs' dynamic and cell-type-specific roles in plant development and stress responses.