Landscape of microRNA and target expression variation and covariation in single mouse embryonic stem cells

Abstract MicroRNAs are small RNA molecules that can repress the expression of protein coding genes post-transcriptionally. Previous studies have shown that microRNAs can also have alternative functions including target noise buffering and co-expression, but these observations have been limited to a few microRNAs. Here we systematically study microRNA alternative functions in mouse embryonic stem cells, by genetically deleting Drosha - leading to global loss of microRNAs. We apply complementary single-cell RNA-seq methods to study the variation of the targets and the microRNAs themselves, and transcriptional inhibition to measure target half-lives. We find that microRNAs form four distinct co-expression groups across single cells. In particular the mir-290 and the mir-182 clusters are abundantly, variably and inversely expressed. Intriguingly, some cells have global biases towards specific miRNAs originating from either end of the hairpin precursor, suggesting the presence of unknown regulatory cofactors. We find that miRNAs generally increase variation and covariation of their targets at the RNA level, but we also find miRNAs such as miR-182 that appear to have opposite functions. In particular, miRNAs that are themselves variable in expression, such as miR-291a, are more likely to induce covariations. In summary, we apply genetic perturbation and multi-omics to give the first global picture of microRNA dynamics at the single cell level.