Cell‐Type Specific miRNA Regulatory Network Responses to ABA Stress Revealed by Time Series Transcriptional Atlases in Arabidopsis

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.