bayesReact: Expression-coupled regulatory motif analysis detects microRNA activity in cancer and at the single cell level

Motivation: Regulatory constraints are crucial in maintaining tissue and cell integrity, and play important roles during developmental processes and environmental responses. Yet many regulatory mechanisms remain unobserved at the single-cell level and statistical inference may, in some cases, help elucidate their condition-specific activity and perturbation during disease progression. Results: We introduce bayesReact (BAYESian modeling of Regular Expression ACTivity), a generative model of motif occurrence across experimentally ranked sequences to infer motif-based regulatory activities. The method is evaluated for microRNAs (miRNAs), which perform post-transcriptional regulation through target mRNA destabilization and translational repression. Inferred miRNA activities positively correlate with the observed miRNA expressions in primary tumors from The Cancer Genome Atlas (TCGA) and mouse stem cells. The top miRNA activity profiles are as informative for TCGA cancer-type cluster identification as the top miRNA or mRNA expression profiles. The activity captures tissue-specific miRNA patterns observed in the matched expression, e.g., the expression of miR-122-5p in the liver and miR-124-3p in low-grade gliomas (LGG). We observe a negative association between the activity of the two miRNAs and their target gene expressions, including between the miR-124-3p activity and the anti-neuronal REST expression in LGG. bayesReact outperforms the existing method, miReact, on sparse count data, and shows a higher correlation with the miRNA expression in single-cell data. The method recovers temporal activities of prominent miRNAs during murine stem cell differentiation, including miR-298-5p, miR-92-2-5p, and the large Sfmbt2 cluster (miR-297-669). The bayesReact model is probabilistic and quantifies the uncertainty of all provided estimates. It is unsupervised and permits screens of bulk or single-cell data to identify condition-specific regulatory motif candidates. It further improves miRNA activity inference in single-cell data. Availability and implementation: bayesReact is implemented as an R-package, uses a Hamiltonian Monte-Carlo sampler for posterior approximation, and is available at https://github.com/astamr/bayesReact.