Comprehensive sequencing profile and functional analysis of IsomiRs in human pancreatic islets and beta cells

Abstract

Aims/hypothesis

MiRNAs regulate gene expression, influencing beta cell function and pathways. Isoforms of miRNA (isomiRs), sequence variants of miRNAs with post-transcriptional modifications, exhibit cell-type-specific expression and functions. Despite their biological significance, a comprehensive isomiR profile in human pancreatic islets and beta cells remains unexplored. This study aims to profile isomiR expression in four beta cell sources: (1) laser capture microdissected human islets (LCM-HI); (2) collagenase-isolated human islets (CI-HI); (3) sorted beta cells; and (4) the EndoC-βH1 beta cell line, and to investigate their potential role in beta cell function.

Methods

Small RNA-seq and/or small RNA dataset analysis was conducted on human pancreatic islets and beta cells. Data were processed using the sRNAbench bioinformatics pipeline to classify isomiRs based on sequence variations. A beta cell-specific isomiR signature was identified via cross-validation across datasets. Correlations between LCM-HI isomiR expression and in vivo clinical parameters were analysed using regression models. Functional validation of isomiR-411-5p-Ext5p(+1) was performed via overexpression in EndoC-βH1 cells and CI-HI, followed by glucose-stimulated insulin secretion (GSIS) assays and/or transcriptomic analysis.

Results

IsomiRs constituted 59.2 ± 1.9% (LCM-HI), 59.6 ± 2.4% (CI-HI), 42.3 ± 7.2% (sorted beta cells) and 43.8 ± 1.2% (EndoC-βH1) of total miRNA reads (data represented as mean ± SD), with 3′ end trimming (Trim3p) being the predominant modification. A beta cell-specific isomiR signature of 30 sequences was identified, with isomiR-411-5p-Ext5p(+1) showing a significant inverse correlation with basal insulin secretion (p=0.0009, partial R²=0.68) and total insulin secretion (p=0.005, partial R²=0.54). Overexpression of isomiR-411-5p-Ext5p(+1), but not of its canonical counterpart, importantly reduced GSIS by 51% ( ± 15.2%; mean ± SD) (p=0.01) in EndoC-βH1 cells. Transcriptomic analysis performed in EndoC-βH1 cells and CI-HI identified 47 genes significantly downregulated by isomiR-411-5p-Ext5p(+1) (false discovery rate [FDR]<0.05) but not by the canonical miRNA, with enriched pathways related to Golgi vesicle biogenesis (FDR=0.017) and trans-Golgi vesicle budding (FDR=0.018). TargetScan analysis confirmed seed sequence-dependent target specificity for 81 genes uniquely regulated by the isomiR (p=1.1 × 10⁻⁹).

Conclusions/interpretation

This study provides the first comprehensive isomiR profiling in human islets and beta cells, revealing their substantial contribution to miRNA regulation. IsomiR-411-5p-Ext5p(+1) emerges as a distinct key modulator of insulin secretion and granule dynamics in beta cells. These findings highlight isomiRs as potential biomarkers and therapeutic targets for diabetes, warranting further exploration of their roles in beta cell biology.

Graphical Abstract