In this study, we conducted a comprehensive meta-analysis of differentially expressed microRNAs (miRNAs) in Arabidopsis thaliana under various abiotic stresses, including heat, drought, and salinity. Categorization of miRNA targets enabled us to investigate regulatory mechanisms and transcriptional pathways associated with each stress type.
Under heat stress, we observed a significant upregulation of miRNAs, including miR156, miR159, and miR398, while miR164 and miR169 were downregulated. These miRNAs have been previously associated with stress adaptation, including the modulation of ROS scavenging (miR398) (Chowdhury et al. 2025) and the regulation of flowering time (Wang, 2014) (miR156). Transcription factor (TF) family enrichment analysis using the TDT hub revealed an overrepresentation of AP2 and BBR/BPC TFs, implying that these families may act as key regulators of heat-responsive miRNA targets (Fig. 1a). STRING-based protein-protein interaction analysis of the miRNA target genes under heat stress identified several core proteins, including BCAT3 (BRANCHED-CHAIN-AMINO-TRANSFERASE-3) and BCAT5 (BRANCHED-CHAIN-AMINO-TRANSFERASE-5), which are involved in branched-chain amino acid metabolism and potentially aid in osmo-protection and energy balance during heat exposure (Diebold et al. 2002). Additionally, SERK (SOMATIC EMBRYOGENESIS RECEPTOR KINASE) and NAC054, both known for their roles in developmental reprogramming and stress signal transduction (Santos et al. 2009, Xu et al. 2024), were found as core components emphasizing their functional relevance in heat-stress adaptation (Fig. 1b).
miRNAs such as miR851, miR771, and miR157b were upregulated under drought stress, potentially contributing to protective gene regulation (Camargo-Ramirez et al. 2018, Shikata et al. 2012), while miR4228, miR172, and miR5020 were downregulated, possibly removing repression on stress-adaptive targets. TF family enrichment of the targets was distinctively lower compared to other stress conditions. However, BBR/BPC remained the most enriched, suggesting a selective regulatory response (Fig. 2a). The STRING PPI network for drought-responsive miRNA targets highlighted key proteins, including LBD37 (LOB-DOMAIN-CONTAINING-37) (Chai et al. 2022), a negative regulator of nitrogen-responsive genes; GPAT9 (GLYCEROL-3-PHOSPHATE-ACYL-TRANSFERASE-9), involved in membrane lipid remodelling (Gong et al. 2023); and PGD1 (PLASTID-GALACTOLIPID-DEGRADATION-1), which contributes to redox balance through the oxidative pentose phosphate pathway (Du et al. 2018). These targets perhaps indicate metabolic adjustments and signalling events required during drought tolerance.
Under salt stress, miR156a (Wang et al. 2023) and miR169e showed marked upregulation, aligning with their known roles in developmental plasticity (Sorin et al. 2014), whereas miR857 and miR778 were downregulated, suggesting a fine-tuned modulation of lignin biosynthesis and transcriptional repression pathways under salinity (Zhao et al. 2015, Bennet et al. 2022). The transcriptional landscape of the targets was visibly wider, with extensive enrichment across multiple TF families (Fig. 3a). Families such as WRKY, YABBY, and zf-HD displayed the most substantial enrichment (with the highest –log10[FDR] values), indicating widespread transcriptional rewiring under salinity conditions. The PPI network highlighted several well-characterized stress-responsive proteins (Fig. 3b), including NADP-ME2 (NADP-MALIC-ENZYME-2), which regulates malate metabolism and cellular pH homeostasis (Badia et al. 2015); ACO1 (1-AMINOCYCLOPROPANE-1-CARBOXYLATE-OXIDASE), a key enzyme in ethylene biosynthesis (Houben et al. 2019); and MYB101 and MYB33 (Wang et al. 2025), both transcription factors involved in ABA signalling and the osmotic stress response. Additionally, LOX2 (LIPOXYGENASE-2) (Yang et al. 2020) and TIFY family proteins, which are repressors of jasmonate signalling (Chung et al. 2009), were observed as network hubs, suggesting that both ethylene and jasmonate pathways are highly active during salt stress responses. Our study also indicates that these pathways also interact with miRNA regulation to shape a plant’s response towards the stress.
Gene ontology enrichment analyses of the miRNA targets revealed a striking and consistent overrepresentation of developmental processes across all three abiotic stresses: heat (Fig. 4), drought (Fig. 5) and salt (Fig. 6). In the Biological Process category (sub part a ), terms such as embryonic meristem development, floral organ development, were prominently enriched, indicating that stress-responsive miRNAs may intersect with developmental reprogramming pathways (Yang et al. 2023) The Cellular Component (subpart b) showed a concentration of targets in compartments mostly associated with transcriptional activity and intercellular communication, such as the secretory vesicles, and plasmodesmata. In the Molecular function category (subpart c), key terms like DNA-binding transcription factor activity, mRNA binding, and ligase activity suggest regulatory flexibility at both transcriptional and post-transcriptional levels. KEGG (Kyoto Encyclopedia of genes and genomes ) analyses (subpart d) of the targets specifically from heat stress highlighted enriched terms including biosynthesis of secondary metabolites and amino acid metabolism, pointing toward active metabolic adjustments in response to stress. Taken together, these results suggest that stress-induced miRNAs not only modulate canonical stress pathways but may also regulate quintessential developmental and metabolic functions, likely to optimize growth and survival under adverse stressors.