Cell-type-specific characterization of miRNA gene dynamics in immune cell subpopulations during aging and atherosclerosis disease development at single-cell resolution

MicroRNAs (miRNAs) are a class of regulatory non-coding RNAs that finetune cellular functions by modulating the stability and abundance of their target mRNAs, thereby contributing to regulation of tissue homeostasis. MiRNA genes are transcribed similarly to protein-coding genes and recent studies have enabled their annotation and quantification genome-wide from bulk nascent transcriptomes. Here, we developed an approach to quantify and integrate miRNA gene signatures into single-cell studies. To characterize miRNA gene expression dynamics, we first compared the suitability of droplet and plate-based single-cell RNA-sequencing (scRNA-seq) platforms using the matched datasets provided by the Tabula Muris Senis and Tabula Sapiens consortiums. We found high concordance between the platforms and with cell type-specific bulk expression data. Based on the comprehensive aging profiles, our analysis comparing spleen immune cells between young and old mice revealed a concordant regulation of miRNAs involved in senescence and inflammatory pathways in multiple immune cell types, including up-regulation of mmu-mir-146a, mmu-mir-101a and mmu-mir-30 family genes. To study the aberrant regulation of immune cell homeostasis and tissue inflammation that pre-dispose to aging-related disease development, we collected transcriptome profiles from atherosclerosis development in LDLR-/-ApoB100/100 mice. We found an elevated myeloid cell proportion in the adipose tissue and further characterized the cell subtypes based on reproducible transcriptome clusters. We then compared miRNA gene expression in early versus late disease and upon inflammatory challenge to monitor different stages during disease progression. At atherosclerotic stage, pro-inflammatory mmu-mir-511 expression increased in several macrophage subtypes, while immunosuppressive mmu-mir-23b[~]mir-24-2[~]mir-27b up-regulation was specific to Trem2+ lipid-associated macrophages. The infiltrating monocytes up-regulated mmu-mir-1938 and mmu-mir-22 expression and in classical monocytes maturation further increased mmu-mir-221[~]222, mmu-mir-511 and mmu-mir-155 expression. To validate that these changes detected from single cell profiles represent miRNA gene transcriptional regulation, we used nascent transcriptomics data from ex vivo macrophage cultures with pro-inflammatory stimulation, confirming both rapid and long-lasting transcriptional activation of the miRNA loci studied. Collectively, our work enables integrating miRNA gene analysis to current single cell genomics pipelines and facilitates characterization of miRNA regulatory networks during aging and disease development.