Characterizing pancreatic β-cell heterogeneity in the streptozotocin model by single-cell transcriptomic analysis

We identified the heterogeneity of β-cells in both physiological and pathological conditions. However, we did not observe conversion of Glut2low to Glut2high β-cells, transcriptomic changes in the non-β endocrine cells, or direct trans-differentiation from the α-cell lineage to the β-cell lineage in the STZ model. Our results clearly define the states of islet cells treated with STZ and allow us to re-evaluate the STZ model widely used in diabetes studies.

We performed single-cell RNA sequencing of pancreatic endocrine cells from STZ-treated mice. High-quality sequencing data from 2,999 cells were used to identify clusters via Louvain clustering analysis. Principal component analysis (PCA), t-distributed stochastic neighbor embedding (t-SNE), uniform manifold approximation and projection (UMAP), force-directed layout (FDL), and differential expression analysis were used to define the heterogeneity and transcriptomic changes in islet cells. In addition, qPCR and immunofluorescence staining were used to confirm findings from the sequencing data.

Untreated β-cells were divided into two populations at the transcriptomic level, a large high-Glut2 expression (Glut2high) population and a small low-Glut2 expression (Glut2low) population. At the transcriptomic level, Glut2low β-cells in adult mice did not represent a developmentally immature state, although a fraction of genes associated with β-cell maturation and function were downregulated in Glut2low cells. After a single high-dose STZ treatment, most Glut2high cells were killed, but Glut2low cells survived and over time changed to a distinct cell state. We did not observe conversion of Glut2low to Glut2high β-cells up to 9 months after STZ treatment. In addition, we did not detect transcriptomic changes in the non-β endocrine cells or a direct trans-differentiation pathway from the α-cell lineage to the β-cell lineage in the STZ model. Conclusions: We identified the heterogeneity of β-cells in both physiological and pathological conditions. However, we did not observe conversion of Glut2low to Glut2high β-cells, transcriptomic changes in the non-β endocrine cells, or direct trans-differentiation from the α-cell lineage to the β-cell lineage in the STZ model. Our results clearly define the states of islet cells treated with STZ and allow us to re-evaluate the STZ model widely used in diabetes studies.