Abstract
MicroRNAs (miRNAs) are post-transcriptional, non-coding regulatory RNAs that function coordinately with transcription factors (TFs) in gene regulatory networks. TFs and their targets are often co-regulated by miRNAs, forming composite feedforward circuits (cFFCs) with varying degrees of redundancy, primarily mediated by miRNAs. However, the maintenance of miRNA-mediated regulatory redundancy and its impact on gene expression evolution remain elusive. By integrating ChIP-seq data from ENCODE and miRNA targeting from TargetScanFly, we quantified miRNA-mediated cFFC redundancy in Drosophila melanogaster embryos and larvae, revealing more than three quarters of miRNA targets are involved in redundant cFFCs. Higher cFFC redundancy, where more miRNAs target the same gene within a cFFC, is correlated with stronger purifying selection, reduced expression divergence between species, and increased expression stability under heat shock stress. Redundant cFFCs primarily regulate older or broadly expressed young genes. These findings highlight the role of miRNA-mediated cFFC redundancy in enhancing gene expression robustness through natural selection.