Bulk and single-cell transcriptome analyses of islet tissue unravel gene signatures associated with pyroptosis and immune infiltration in type 2 diabetes

Introduction: Type 2 diabetes (T2D) is a common chronic heterogeneous metabolic disorder. However, the roles of pyroptosis and infiltrating immune cells in islet dysfunction of patients with T2D have yet to be explored. In this study, we aimed to explore potential crucial genes and pathways associated with pyroptosis and immune infiltration in T2D.

Methods: To achieve this, we performed a conjoint analysis of three bulk RNA-seq datasets of islets to identify T2D-related differentially expressed genes (DEGs). After grouping the islet samples according to their ESTIMATE immune scores, we identified immune- and T2D-related DEGs. A clinical prediction model based on pyroptosis-related genes for T2D was constructed. Weighted gene co-expression network analysis was performed to identify genes positively correlated with pyroptosis-related pathways. A protein–protein interaction network was established to identify pyroptosis-related hub genes. We constructed miRNA and transcriptional networks based on the pyroptosis-related hub genes and performed functional analyses. Single-cell RNA-seq (scRNA-seq) was conducted using the GSE153885 dataset. Dimensionality was reduced using principal component analysis and t-distributed statistical neighbor embedding, and cells were clustered using Seurat. Different cell types were subjected to differential gene expression analysis and gene set enrichment analysis (GSEA). Cell–cell communication and pseudotime trajectory analyses were conducted using the samples from patients with T2D.

Results: We identified 17 pyroptosis-related hub genes. We determined the abundance of 13 immune cell types in the merged matrix and found that these cell types were correlated with the 17 pyroptosis-related hub genes. Analysis of the scRNA-seq dataset of 1892 islet samples from patients with T2D and controls revealed 11 clusters. INS and IAPP were determined to be pyroptosis-related and candidate hub genes among the 11 clusters. GSEA of the 11 clusters demonstrated that the myc, G2M checkpoint, and E2F pathways were significantly upregulated in clusters with several differentially enriched pathways.

Discussion: This study elucidates the gene signatures associated with pyroptosis and immune infiltration in T2D and provides a critical resource for understanding of islet dysfunction and T2D pathogenesis.