# miRglmm: a generalized linear mixed model of isomiR-level counts improves estimation of miRNA-level differential expression and uncovers variable differential expression between isomiRs

**Authors:** Andrea M. Baran, Arun H. Patil, Ernesto Aparicio-Puerta, Seong-Hwan Jun, Marc K. Halushka, Matthew N. McCall

PMC · DOI: 10.1186/s13059-025-03549-y · Genome Biology · 2025-04-22

## TL;DR

This paper introduces miRglmm, a new method for analyzing microRNA sequencing data that improves detection of differential expression and isomiR usage.

## Contribution

miRglmm introduces a generalized linear mixed model for isomiR-level counts, improving miRNA differential expression estimation.

## Key findings

- miRglmm outperforms existing methods in estimating miRNA differential expression.
- miRglmm provides estimates of isomiR-level differential expression.
- The method handles differential isomiR usage effectively.

## Abstract

MicroRNA-seq data is produced by aligning small RNA sequencing reads of different microRNA transcript isoforms, called isomiRs, to known microRNAs. Aggregation to microRNA-level counts discards information and violates core assumptions of differential expression methods developed for mRNA-seq data. We establish miRglmm, a differential expression method for microRNA-seq data, that uses a generalized linear mixed model of isomiR-level counts, facilitating detection of miRNA with differential expression or differential isomiR usage. We demonstrate that miRglmm outperforms current differential expression methods in estimating differential expression for miRNA, whether or not there is differential isomiR usage, and simultaneously provides estimates of isomiR-level differential expression.

The online version contains supplementary material available at 10.1186/s13059-025-03549-y.

## Full-text entities

- **Genes:** Mir26a-1 (microRNA 26a-1) [NCBI Gene 387218] {aka Mirn26a, Mirn26a-1, miR-26a, mir-26a-1}, Cd19 (CD19 antigen) [NCBI Gene 12478], Mir191 (microRNA 191) [NCBI Gene 387186] {aka Mirn191, mir-191, mmu-mir-191}, Mir223 (microRNA 223) [NCBI Gene 723814] {aka Mirn223, miR-223, mmu-mir-223}, Mir378a (microRNA 378a) [NCBI Gene 723889] {aka Mir378, Mirn378, mmu-mir-378, mmu-mir-378a}, MIRLET7B (microRNA let-7b) [NCBI Gene 406884] {aka LET7B, MIRNLET7B, hsa-let-7b, let-7b}, MIRLET7C (microRNA let-7c) [NCBI Gene 406885] {aka LET7C, MIRNLET7C, hsa-let-7c, let-7c}, Mir181a-2 (microRNA 181a-2) [NCBI Gene 387176] {aka Mirn181, Mirn181a, Mirn181a-2, miR-181, mir-181a, mir-181a-2}, Mir150 (microRNA 150) [NCBI Gene 387168] {aka Mirn150, mir-150, mmu-mir-150}, Cd4 (CD4 antigen) [NCBI Gene 12504] {aka L3T4, Ly-4}
- **Diseases:** cancer (MESH:D009369), hypoxia (MESH:D000860)
- **Chemicals:** adenosine (MESH:D000241), uracil (MESH:D014498)
- **Species:** Mus musculus (house mouse, species) [taxon 10090], Homo sapiens (human, species) [taxon 9606]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12016310/full.md

## Figures

3 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12016310/full.md

## References

10 references — full list in the complete paper: https://tomesphere.com/paper/PMC12016310/full.md

---
Source: https://tomesphere.com/paper/PMC12016310