Calcium-dependent transcriptional changes in human pancreatic islet cells reveal functional diversity in islet cell subtypes.

AIMS/HYPOTHESIS: Pancreatic islets depend on cytosolic calcium (Ca2+) to trigger the secretion of glucoregulatory hormones and trigger transcriptional regulation of genes important for islet response to stimuli. To date, there has not been an attempt to profile Ca2+-regulated gene expression in all islet cell types. Our aim was to construct a large single-cell transcriptomic dataset from human islets exposed to conditions that would acutely induce or inhibit intracellular Ca2+ signalling, while preserving biological heterogeneity. METHODS: We exposed intact human islets from three donors to the following conditions: (1) 2.8 mmol/l glucose; (2) 16 mmol/l glucose and 40 mmol/l KCl to maximally stimulate Ca2+ signalling; and (3) 16 mmol/l glucose, 40 mmol/l KCl and 5 mmol/l EGTA (Ca2+ chelator) to inhibit Ca2+ signalling, for 1 h. We sequenced 68,650 cells from all islet cell types, and further subsetted the cells to form an endocrine cell-specific dataset of 59,373 cells expressing INS, GCG, SST or PPY. We compared transcriptomes across conditions to determine the differentially expressed Ca2+-regulated genes in each endocrine cell type, and in each endocrine cell subcluster of alpha and beta cells. RESULTS: Based on the number of Ca2+-regulated genes, we found that each alpha and beta cell cluster had a different magnitude of Ca2+ response. We also showed that polyhormonal clusters expressing both INS and GCG, or both INS and SST, are defined by Ca2+-regulated genes specific to each cluster. Finally, we identified the gene PCDH7 from the beta cell clusters that had the highest number of Ca2+-regulated genes, and showed that cells expressing cell surface PCDH7 protein have enhanced glucose-stimulated insulin secretory function. CONCLUSIONS/INTERPRETATION: Here we use our large-scale, multi-condition, single-cell dataset to show that human islets have cell-type-specific Ca2+-regulated gene expression profiles, some of them specific to subpopulations. In our dataset, we identify PCDH7 as a novel marker of beta cells having an increased number of Ca2+-regulated genes and enhanced insulin secretory function. DATA AVAILABILITY: A searchable and user-friendly format of the data in this study, specifically designed for rapid mining of single-cell RNA sequencing data, is available at https://lynnlab.shinyapps.io/Human_Islet_Atlas/ . The raw data files are available at NCBI Gene Expression Omnibus (GSE196715).

Spatial and transcriptional heterogeneity of pancreatic beta cell neogenesis revealed by a time-resolved reporter system.

AIMS/HYPOTHESIS: While pancreatic beta cells have been shown to originate from endocrine progenitors in ductal regions, it remains unclear precisely where beta cells emerge from and which transcripts define newborn beta cells. We therefore investigated characteristics of newborn beta cells extracted by a time-resolved reporter system. METHODS: We established a mouse model, 'Ins1-GFP; Timer', which provides spatial information during beta cell neogenesis with high temporal resolution. Single-cell RNA-sequencing (scRNA-seq) was performed on mouse beta cells sorted by fluorescent reporter to uncover transcriptomic profiles of newborn beta cells. scRNA-seq of human embryonic stem cell (hESC)-derived beta-like cells was also performed to compare newborn beta cell features between mouse and human. RESULTS: Fluorescence imaging of Ins1-GFP; Timer mouse pancreas successfully dissected newly generated beta cells as green fluorescence-dominant cells. This reporter system revealed that, as expected, some newborn beta cells arise close to the ducts (ßduct); unexpectedly, the others arise away from the ducts and adjacent to blood vessels (ßvessel). Single-cell transcriptomic analyses demonstrated five distinct populations among newborn beta cells, confirming spatial heterogeneity of beta cell neogenesis such as high probability of glucagon-positive ßduct, musculoaponeurotic fibrosarcoma oncogene family B (MafB)-positive ßduct and musculoaponeurotic fibrosarcoma oncogene family A (MafA)-positive ßvessel cells. Comparative analysis with scRNA-seq data of mouse newborn beta cells and hESC-derived beta-like cells uncovered transcriptional similarity between mouse and human beta cell neogenesis including microsomal glutathione S-transferase 1 (MGST1)- and synaptotagmin 13 (SYT13)-highly-expressing state. CONCLUSIONS/INTERPRETATION: The combination of time-resolved histological imaging with single-cell transcriptional mapping demonstrated novel features of spatial and transcriptional heterogeneity in beta cell neogenesis, which will lead to a better understanding of beta cell differentiation for future cell therapy. DATA AVAILABILITY: Raw and processed single-cell RNA-sequencing data for this study has been deposited in the Gene Expression Omnibus under accession number GSE155742.

Type 2 diabetes susceptibility gene GRK5 regulates physiological pancreatic ß-cell proliferation via phosphorylation of HDAC5.

Restoring functional ß cell mass is a potential therapy for those with diabetes. However, the pathways regulating ß cell mass are not fully understood. Previously, we demonstrated that Sox4 is required for ß cell proliferation during prediabetes. Here, we report that Sox4 regulates ß cell mass through modulating expression of the type 2 diabetes (T2D) susceptibility gene GRK5. ß cell-specific Grk5 knockout mice showed impaired glucose tolerance with reduced ß cell mass, which was accompanied by upregulation of cell cycle inhibitor gene Cdkn1a. Furthermore, we found that Grk5 may drive ß cell proliferation through a pathway that includes phosphorylation of HDAC5 and subsequent transcription of immediate-early genes (IEGs) such as Nr4a1, Fosb, Junb, Arc, Egr1, and Srf. Together, these studies suggest GRK5 is linked to T2D through regulation of ß cell growth and that it may be a target to preserve ß cells during the development of T2D.