Patch-Seq Links Single-Cell Transcriptomes to Human Islet Dysfunction in Diabetes

Loss of ATP-Sensitive Potassium Channel Expression and Function in the Nervous System Decreases Opioid Sensitivity in a High-Fat Diet-Fed Mouse Model of Diet-Induced Obesity

During diabetes progression, β-cell dysfunction due to loss of potassium channels sensitive to ATP, known as KATP channels, occurs, contributing to hyperglycemia. The aim of this study was to investigate if KATP channel expression or activity in the nervous system was altered in a high-fat diet (HFD)-fed mouse model of diet-induced obesity. Expression of two KATP channel subunits, Kcnj11 (Kir6.2) and Abcc8 (SUR1), were decreased in the peripheral and central nervous system of mice fed HFD, which was significantly correlated with mechanical paw-withdrawal thresholds. HFD mice had decreased antinociception to systemic morphine compared with control diet (CON) mice, which was expected because KATP channels are downstream targets of opioid receptors. Mechanical hypersensitivity in HFD mice was exacerbated after systemic treatment with glyburide or nateglinide, KATP channel antagonists clinically used to control blood glucose levels. Upregulation of SUR1 and Kir6.2, through an adenovirus delivered intrathecally, increased morphine antinociception in HFD mice. These data present a potential link between KATP channel function and neuropathy during early stages of diabetes. There is a need for increased knowledge of how diabetes affects structural and molecular changes in the nervous system, including ion channels, to lead to the progression of chronic pain and sensory issues.

ARTICLE HIGHLIGHTS

Longitudinal integrative cell-free DNA analysis in gestational diabetes mellitus

Gestational diabetes mellitus (GDM) presents varied manifestations throughout pregnancy and poses a complex clinical challenge. High-depth cell-free DNA (cfDNA) sequencing analysis holds promise in advancing our understanding of GDM pathogenesis and prediction. In 299 women with GDM and 299 matched healthy pregnant women, distinct cfDNA fragment characteristics associated with GDM are identified throughout pregnancy. Integrating cfDNA profiles with lipidomic and single-cell transcriptomic data elucidates functional changes linked to altered lipid metabolism processes in GDM. Transcription start site (TSS) scores in 50 feature genes are used as the cfDNA signature to distinguish GDM cases from controls effectively. Notably, differential coverage of the islet acinar marker gene PRSS1 emerges as a valuable biomarker for GDM. A specialized neural network model is developed, predicting GDM occurrence and validated across two independent cohorts. This research underscores the high-depth cfDNA early prediction and characterization of GDM, offering insights into its molecular underpinnings and potential clinical applications.

Single-cell transcriptomic profiling of human pancreatic islets reveals genes responsive to glucose exposure over 24 h

AIMS/HYPOTHESIS

Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycaemia, beta cell glucotoxicity and subsequently type 2 diabetes. In this study, we explored the effects of in vitro hyperglycaemic conditions on human pancreatic islet gene expression across 24 h in six pancreatic cell types: alpha; beta; gamma; delta; ductal; and acinar. We hypothesised that genes associated with hyperglycaemic conditions may be relevant to the onset and progression of diabetes.

METHODS

We exposed human pancreatic islets from two donors to low (2.8 mmol/l) and high (15.0 mmol/l) glucose concentrations over 24 h in vitro. To assess the transcriptome, we performed single-cell RNA-seq (scRNA-seq) at seven time points. We modelled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Additionally, we integrated genomic features and genetic summary statistics to nominate candidate effector genes. For three of these genes, we functionally characterised the effect on insulin production and secretion using CRISPR interference to knock down gene expression in EndoC-βH1 cells, followed by a glucose-stimulated insulin secretion assay.

RESULTS

In the discrete time models, we identified 1344 genes associated with time and 668 genes associated with glucose exposure across all cell types and time points. In the continuous time models, we identified 1311 genes associated with time, 345 genes associated with glucose exposure and 418 genes associated with interaction effects between time and glucose across all cell types. By integrating these expression profiles with summary statistics from genetic association studies, we identified 2449 candidate effector genes for type 2 diabetes, HbA1c, random blood glucose and fasting blood glucose. Of these candidate effector genes, we showed that three (ERO1B, HNRNPA2B1 and RHOBTB3) exhibited an effect on glucose-stimulated insulin production and secretion in EndoC-βH1 cells.

CONCLUSIONS/INTERPRETATION

The findings of our study provide an in-depth characterisation of the 24 h transcriptomic response of human pancreatic islets to glucose exposure at a single-cell resolution. By integrating differentially expressed genes with genetic signals for type 2 diabetes and glucose-related traits, we provide insights into the molecular mechanisms underlying glucose homeostasis. Finally, we provide functional evidence to support the role of three candidate effector genes in insulin secretion and production.

DATA AVAILABILITY

The scRNA-seq data from the 24 h glucose exposure experiment performed in this study are available in the database of Genotypes and Phenotypes (dbGap; https://www.ncbi.nlm.nih.gov/gap/ ) with accession no. phs001188.v3.p1. Study metadata and summary statistics for the differential expression, gene set enrichment and candidate effector gene prediction analyses are available in the Zenodo data repository ( https://zenodo.org/ ) under accession number 11123248. The code used in this study is publicly available at https://github.com/CollinsLabBioComp/publication-islet_glucose_timecourse .

Identification of unique cell type responses in pancreatic islets to stress

Diabetes involves the death or dysfunction of pancreatic β-cells. Analysis of bulk sequencing from human samples and studies using in vitro and in vivo models suggest that endoplasmic reticulum and inflammatory signaling play an important role in diabetes progression. To better characterize cell type-specific stress response, we perform multiplexed single-cell RNA sequencing to define the transcriptional signature of primary human islet cells exposed to endoplasmic reticulum and inflammatory stress. Through comprehensive pair-wise analysis of stress responses across pancreatic endocrine and exocrine cell types, we define changes in gene expression for each cell type under different diabetes-associated stressors. We find that β-, α-, and ductal cells have the greatest transcriptional response. We utilize stem cell-derived islets to study islet health through the candidate gene CIB1, which was upregulated under stress in primary human islets. Our findings provide insights into cell type-specific responses to diabetes-associated stress and establish a resource to identify targets for diabetes therapeutics.

Optogenetic β cell interrogation in vivo reveals a functional hierarchy directing the Ca2+ response to glucose supported by vitamin B6

Coordination of cellular activity through Ca2+ enables β cells to secrete precise quantities of insulin. To explore how the Ca2+ response is orchestrated in space and time, we implement optogenetic systems to probe the role of individual β cells in the glucose response. By targeted β cell activation/inactivation in zebrafish, we reveal a hierarchy of cells, each with a different level of influence over islet-wide Ca2+ dynamics. First-responder β cells lie at the top of the hierarchy, essential for initiating the first-phase Ca2+ response. Silencing first responders impairs the Ca2+ response to glucose. Conversely, selective activation of first responders demonstrates their increased capability to raise pan-islet Ca2+ levels compared to followers. By photolabeling and transcriptionally profiling β cells that differ in their thresholds to a glucose-stimulated Ca2+ response, we highlight vitamin B6 production as a signature pathway of first responders. We further define an evolutionarily conserved requirement for vitamin B6 in enabling the Ca2+ response to glucose in mammalian systems.

HumanIslets: An integrated platform for human islet data access and analysis

Comprehensive molecular and cellular phenotyping of human islets can enable deep mechanistic insights for diabetes research. We established the Human Islet Data Analysis and Sharing (HI-DAS) consortium to advance goals in accessibility, usability, and integration of data from human islets isolated from donors with and without diabetes at the Alberta Diabetes Institute (ADI) IsletCore. Here we introduce HumanIslets.com, an open resource for the research community. This platform, which presently includes data on 547 human islet donors, allows users to access linked datasets describing molecular profiles, islet function and donor phenotypes, and to perform various statistical and functional analyses at the donor, islet and single-cell levels. As an example of the analytic capacity of this resource we show a dissociation between cell culture effects on transcript and protein expression, and an approach to correct for exocrine contamination found in hand-picked islets. Finally, we provide an example workflow and visualization that highlights links between type 2 diabetes status, SERCA3b Ca2+-ATPase levels at the transcript and protein level, insulin secretion and islet cell phenotypes. HumanIslets.com provides a growing and adaptable set of resources and tools to support the metabolism and diabetes research community.

Expression of FAM159B in Humans, Rats, and Mice: A Cross-species Examination

Little is known about the adaptor protein FAM159B. To determine whether FAM159B expression findings in rats or mice can be extrapolated to humans, we compared FAM159B expression in healthy tissue samples from all three species using immunohistochemistry. Despite variations in expression intensity, similar FAM159B expression patterns were observed in most organs across species. The most prominent species difference was noted in pancreatic islets; while FAM159B expression was limited to single cells on the outer edges in mice and rats, it was detectable across entire islets in humans. Double-labeling immunohistochemistry revealed partial overlap of FAM159B expression with that of insulin, glucagon, and somatostatin in human islets. By contrast, FAM159B showed complete colocalization with only somatostatin in rats and mice. An additional analysis of FAM159B expression in lean and obese Zucker rats revealed larger islet areas due to increased β-cell mass in obese rats, which was accompanied by a smaller percentage of FAM159B-positive δ-cells per islet area. Beyond the known differences in islet architecture across species, our results point to larger dissimilarities in blood glucose regulation between rodents and humans than generally assumed. Moreover, findings regarding FAM159B expression (and function) cannot be directly transferred between rodents and humans.

Differential CpG methylation at Nnat in the early establishment of beta cell heterogeneity

AIMS/HYPOTHESIS

Beta cells within the pancreatic islet represent a heterogenous population wherein individual sub-groups of cells make distinct contributions to the overall control of insulin secretion. These include a subpopulation of highly connected 'hub' cells, important for the propagation of intercellular Ca2+ waves. Functional subpopulations have also been demonstrated in human beta cells, with an altered subtype distribution apparent in type 2 diabetes. At present, the molecular mechanisms through which beta cell hierarchy is established are poorly understood. Changes at the level of the epigenome provide one such possibility, which we explore here by focusing on the imprinted gene Nnat (encoding neuronatin [NNAT]), which is required for normal insulin synthesis and secretion.

METHODS

Single-cell RNA-seq datasets were examined using Seurat 4.0 and ClusterProfiler running under R. Transgenic mice expressing enhanced GFP under the control of the Nnat enhancer/promoter regions were generated for FACS of beta cells and downstream analysis of CpG methylation by bisulphite sequencing and RNA-seq, respectively. Animals deleted for the de novo methyltransferase DNA methyltransferase 3 alpha (DNMT3A) from the pancreatic progenitor stage were used to explore control of promoter methylation. Proteomics was performed using affinity purification mass spectrometry and Ca2+ dynamics explored by rapid confocal imaging of Cal-520 AM and Cal-590 AM. Insulin secretion was measured using homogeneous time-resolved fluorescence imaging.

RESULTS

Nnat mRNA was differentially expressed in a discrete beta cell population in a developmental stage- and DNA methylation (DNMT3A)-dependent manner. Thus, pseudo-time analysis of embryonic datasets demonstrated the early establishment of Nnat-positive and -negative subpopulations during embryogenesis. NNAT expression is also restricted to a subset of beta cells across the human islet that is maintained throughout adult life. NNAT+ beta cells also displayed a discrete transcriptome at adult stages, representing a subpopulation specialised for insulin production, and were diminished in db/db mice. 'Hub' cells were less abundant in the NNAT+ population, consistent with epigenetic control of this functional specialisation.

CONCLUSIONS/INTERPRETATION

These findings demonstrate that differential DNA methylation at Nnat represents a novel means through which beta cell heterogeneity is established during development. We therefore hypothesise that changes in methylation at this locus may contribute to a loss of beta cell hierarchy and connectivity, potentially contributing to defective insulin secretion in some forms of diabetes.

DATA AVAILABILITY

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD048465.

Advancing Diabetes Research: A Novel Islet Isolation Method from Living Donors

Pancreatic islet isolation is critical for type 2 diabetes research. Although -omics approaches have shed light on islet molecular profiles, inconsistencies persist; on the other hand, functional studies are essential, but they require reliable and standardized isolation methods. Here, we propose a simplified protocol applied to very small-sized samples collected from partially pancreatectomized living donors. Islet isolation was performed by digesting tissue specimens collected during surgery within a collagenase P solution, followed by a Lympholyte density gradient separation; finally, functional assays and staining with dithizone were carried out. Isolated pancreatic islets exhibited functional responses to glucose and arginine stimulation mirroring donors' metabolic profiles, with insulin secretion significantly decreasing in diabetic islets compared to non-diabetic islets; conversely, proinsulin secretion showed an increasing trend from non-diabetic to diabetic islets. This novel islet isolation method from living patients undergoing partial pancreatectomy offers a valuable opportunity for targeted study of islet physiology, with the primary advantage of being time-effective and successfully preserving islet viability and functionality. It enables the generation of islet preparations that closely reflect donors' clinical profiles, simplifying the isolation process and eliminating the need for a Ricordi chamber. Thus, this method holds promises for advancing our understanding of diabetes and for new personalized pharmacological approaches.

Causal network perturbation analysis identifies known and novel type-2 diabetes driver genes

The molecular pathogenesis of diabetes is multifactorial, involving genetic predisposition and environmental factors that are not yet fully understood. However, pancreatic β-cell failure remains among the primary reasons underlying the progression of type-2 diabetes (T2D) making targeting β-cell dysfunction an attractive pathway for diabetes treatment. To identify genetic contributors to β-cell dysfunction, we investigated single-cell gene expression changes in β-cells from healthy (C57BL/6J) and diabetic (NZO/HlLtJ) mice fed with normal or high-fat, high-sugar diet (HFHS). Our study presents an innovative integration of the causal network perturbation assessment (ssNPA) framework with meta-cell transcriptome analysis to explore the genetic underpinnings of type-2 diabetes (T2D). By generating a reference causal network and in silico perturbation, we identified novel genes implicated in T2D and validated our candidates using the Knockout Mouse Phenotyping (KOMP) Project database.

Autoimmune CD8+ T cells in type 1 diabetes: from single-cell RNA sequencing to T-cell receptor redirection

Type 1 diabetes (T1D) is an organ-specific autoimmune disease caused by pancreatic β cell destruction and mediated primarily by autoreactive CD8+ T cells. It has been shown that only a small number of stem cell-like β cell-specific CD8+ T cells are needed to convert normal mice into T1D mice; thus, it is likely that T1D can be cured or significantly improved by modulating or altering self-reactive CD8+ T cells. However, stem cell-type, effector and exhausted CD8+ T cells play intricate and important roles in T1D. The highly diverse T-cell receptors (TCRs) also make precise and stable targeted therapy more difficult. Therefore, this review will investigate the mechanisms of autoimmune CD8+ T cells and TCRs in T1D, as well as the related single-cell RNA sequencing (ScRNA-Seq), CRISPR/Cas9, chimeric antigen receptor T-cell (CAR-T) and T-cell receptor-gene engineered T cells (TCR-T), for a detailed and clear overview. This review highlights that targeting CD8+ T cells and their TCRs may be a potential strategy for predicting or treating T1D.

Acinar to β-like cell conversion through inhibition of focal adhesion kinase

Insufficient functional β-cell mass causes diabetes; however, an effective cell replacement therapy for curing diabetes is currently not available. Reprogramming of acinar cells toward functional insulin-producing cells would offer an abundant and autologous source of insulin-producing cells. Our lineage tracing studies along with transcriptomic characterization demonstrate that treatment of adult mice with a small molecule that specifically inhibits kinase activity of focal adhesion kinase results in trans-differentiation of a subset of peri-islet acinar cells into insulin producing β-like cells. The acinar-derived insulin-producing cells infiltrate the pre-existing endocrine islets, partially restore β-cell mass, and significantly improve glucose homeostasis in diabetic mice. These findings provide evidence that inhibition of the kinase activity of focal adhesion kinase can convert acinar cells into insulin-producing cells and could offer a promising strategy for treating diabetes.

Checkpoint kinase 2 controls insulin secretion and glucose homeostasis

After the discovery of insulin, a century ago, extensive work has been done to unravel the molecular network regulating insulin secretion. Here we performed a chemical screen and identified AZD7762, a compound that potentiates glucose-stimulated insulin secretion (GSIS) of a human β cell line, healthy and type 2 diabetic (T2D) human islets and primary cynomolgus macaque islets. In vivo studies in diabetic mouse models and cynomolgus macaques demonstrated that AZD7762 enhances GSIS and improves glucose tolerance. Furthermore, genetic manipulation confirmed that ablation of CHEK2 in human β cells results in increased insulin secretion. Consistently, high-fat-diet-fed Chk2-/- mice show elevated insulin secretion and improved glucose clearance. Finally, untargeted metabolic profiling demonstrated the key role of the CHEK2-PP2A-PLK1-G6PD-PPP pathway in insulin secretion. This study successfully identifies a previously unknown insulin secretion regulating pathway that is conserved across rodents, cynomolgus macaques and human β cells in both healthy and T2D conditions.

Pancreatic β-Cell Identity Change through the Lens of Single-Cell Omics Research

The main hallmark in the development of both type 1 and type 2 diabetes is a decline in functional β-cell mass. This decline is predominantly attributed to β-cell death, although recent findings suggest that the loss of β-cell identity may also contribute to β-cell dysfunction. This phenomenon is characterized by a reduced expression of key markers associated with β-cell identity. This review delves into the insights gained from single-cell omics research specifically focused on β-cell identity. It highlights how single-cell omics based studies have uncovered an unexpected level of heterogeneity among β-cells and have facilitated the identification of distinct β-cell subpopulations through the discovery of cell surface markers, transcriptional regulators, the upregulation of stress-related genes, and alterations in chromatin activity. Furthermore, specific subsets of β-cells have been identified in diabetes, such as displaying an immature, dedifferentiated gene signature, expressing significantly lower insulin mRNA levels, and expressing increased β-cell precursor markers. Additionally, single-cell omics has increased insight into the detrimental effects of diabetes-associated conditions, including endoplasmic reticulum stress, oxidative stress, and inflammation, on β-cell identity. Lastly, this review outlines the factors that may influence the identification of β-cell subpopulations when designing and performing a single-cell omics experiment.

The Importance of Intra-Islet Communication in the Function and Plasticity of the Islets of Langerhans during Health and Diabetes

Islets of Langerhans are anatomically dispersed within the pancreas and exhibit regulatory coordination between islets in response to nutritional and inflammatory stimuli. However, within individual islets, there is also multi-faceted coordination of function between individual beta-cells, and between beta-cells and other endocrine and vascular cell types. This is mediated partly through circulatory feedback of the major secreted hormones, insulin and glucagon, but also by autocrine and paracrine actions within the islet by a range of other secreted products, including somatostatin, urocortin 3, serotonin, glucagon-like peptide-1, acetylcholine, and ghrelin. Their availability can be modulated within the islet by pericyte-mediated regulation of microvascular blood flow. Within the islet, both endocrine progenitor cells and the ability of endocrine cells to trans-differentiate between phenotypes can alter endocrine cell mass to adapt to changed metabolic circumstances, regulated by the within-islet trophic environment. Optimal islet function is precariously balanced due to the high metabolic rate required by beta-cells to synthesize and secrete insulin, and they are susceptible to oxidative and endoplasmic reticular stress in the face of high metabolic demand. Resulting changes in paracrine dynamics within the islets can contribute to the emergence of Types 1, 2 and gestational diabetes.

q-Diffusion leverages the full dimensionality of gene coexpression in single-cell transcriptomics

Unlocking the full dimensionality of single-cell RNA sequencing data (scRNAseq) is the next frontier to a richer, fuller understanding of cell biology. We introduce q-diffusion, a framework for capturing the coexpression structure of an entire library of genes, improving on state-of-the-art analysis tools. The method is demonstrated via three case studies. In the first, q-diffusion helps gain statistical significance for differential effects on patient outcomes when analyzing the CALGB/SWOG 80405 randomized phase III clinical trial, suggesting precision guidance for the treatment of metastatic colorectal cancer. Secondly, q-diffusion is benchmarked against existing scRNAseq classification methods using an in vitro PBMC dataset, in which the proposed method discriminates IFN-γ stimulation more accurately. The same case study demonstrates improvements in unsupervised cell clustering with the recent Tabula Sapiens human atlas. Finally, a local distributional segmentation approach for spatial scRNAseq, driven by q-diffusion, yields interpretable structures of human cortical tissue.

Single-cell sequencing: A promising approach for uncovering the characteristic of pancreatic islet cells in type 2 diabetes

Single-cell sequencing is a novel and rapidly advancing high-throughput technique that can be used to investigating genomics, transcriptomics, and epigenetics at a single-cell level. Currently, single-cell sequencing can not only be used to draw the pancreatic islet cells map and uncover the characteristics of cellular heterogeneity in type 2 diabetes, but can also be used to label and purify functional beta cells in pancreatic stem cells, improving stem cells and islet organoids therapies. In addition, this technology helps to analyze islet cell dedifferentiation and can be applied to the treatment of type 2 diabetes. In this review, we summarize the development and process of single-cell sequencing, describe the potential applications of single-cell sequencing in the field of type 2 diabetes, and discuss the prospects and limitations of single-cell sequencing to provide a new direction for exploring the pathogenesis of type 2 diabetes and finding therapeutic targets.

Endocrine islet β-cell subtypes with differential function are derived from biochemically distinct embryonic endocrine islet progenitors that are regulated by maternal nutrients

Endocrine islet b cells comprise heterogenous cell subsets. Yet when/how these subsets are produced and how stable they are remain unknown. Addressing these questions is important for preventing/curing diabetes, because lower numbers of b cells with better secretory function is a high risk of this disease. Using combinatorial cell lineage tracing, scRNA-seq, and DNA methylation analysis, we show here that embryonic islet progenitors with distinct gene expression and DNA methylation produce b-cell subtypes of different function and viability in adult mice. The subtype with better function is enriched for genes involved in vesicular production/trafficking, stress response, and Ca2+-secretion coupling, which further correspond to differential DNA methylation in putative enhancers of these genes. Maternal overnutrition, a major diabetes risk factor, reduces the proportion of endocrine progenitors of the b-cell subtype with better-function via deregulating DNA methyl transferase 3a. Intriguingly, the gene signature that defines mouse b-cell subtypes can reliably divide human cells into two sub-populations while the proportion of b cells with better-function is reduced in diabetic donors. The implication of these results is that modulating DNA methylation in islet progenitors using maternal food supplements can be explored to improve b-cell function in the prevention and therapy of diabetes.

RFX6 Maintains Gene Expression and Function of Adult Human Islet α-Cells

Mutations in the gene encoding the transcription factor regulatory factor X-box binding 6 (RFX6) are associated with human diabetes. Within pancreatic islets, RFX6 expression is most abundant in islet α-cells, and α-cell RFX6 expression is altered in diabetes. However, the roles of RFX6 in regulating gene expression, glucagon output, and other crucial human adult α-cell functions are not yet understood. We developed a method for selective genetic targeting of human α-cells and assessed RFX6-dependent α-cell function. RFX6 suppression with RNA interference led to impaired α-cell exocytosis and dysregulated glucagon secretion in vitro and in vivo. By contrast, these phenotypes were not observed with RFX6 suppression across all islet cells. Transcriptomics in α-cells revealed RFX6-dependent expression of genes governing nutrient sensing, hormone processing, and secretion, with some of these exclusively expressed in human α-cells. Mapping of RFX6 DNA-binding sites in primary human islet cells identified a subset of direct RFX6 target genes. Together, these data unveil RFX6-dependent genetic targets and mechanisms crucial for regulating adult human α-cell function.

ARTICLE HIGHLIGHTS

Gene expression analysis reveals diabetes-related gene signatures

BACKGROUND

Diabetes is a spectrum of metabolic diseases affecting millions of people worldwide. The loss of pancreatic β-cell mass by either autoimmune destruction or apoptosis, in type 1-diabetes (T1D) and type 2-diabetes (T2D), respectively, represents a pathophysiological process leading to insulin deficiency. Therefore, therapeutic strategies focusing on restoring β-cell mass and β-cell insulin secretory capacity may impact disease management. This study took advantage of powerful integrative bioinformatic tools to scrutinize publicly available diabetes-associated gene expression data to unveil novel potential molecular targets associated with β-cell dysfunction.

METHODS

A comprehensive literature search for human studies on gene expression alterations in the pancreas associated with T1D and T2D was performed. A total of 6 studies were selected for data extraction and for bioinformatic analysis. Pathway enrichment analyses of differentially expressed genes (DEGs) were conducted, together with protein-protein interaction networks and the identification of potential transcription factors (TFs). For noncoding differentially expressed RNAs, microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), which exert regulatory activities associated with diabetes, identifying target genes and pathways regulated by these RNAs is fundamental for establishing a robust regulatory network.

RESULTS

Comparisons of DEGs among the 6 studies showed 59 genes in common among 4 or more studies. Besides alterations in mRNA, it was possible to identify differentially expressed miRNA and lncRNA. Among the top transcription factors (TFs), HIPK2, KLF5, STAT1 and STAT3 emerged as potential regulators of the altered gene expression. Integrated analysis of protein-coding genes, miRNAs, and lncRNAs pointed out several pathways involved in metabolism, cell signaling, the immune system, cell adhesion, and interactions. Interestingly, the GABAergic synapse pathway emerged as the only common pathway to all datasets.

CONCLUSIONS

This study demonstrated the power of bioinformatics tools in scrutinizing publicly available gene expression data, thereby revealing potential therapeutic targets like the GABAergic synapse pathway, which holds promise in modulating α-cells transdifferentiation into β-cells.

Integrating single-cell transcriptomics with cellular phenotypes: cell morphology, Ca2+ imaging and electrophysiology

I review recent technological advancements in coupling single-cell transcriptomics with cellular phenotypes including morphology, calcium signaling, and electrophysiology. Single-cell RNA sequencing (scRNAseq) has revolutionized cell type classifications by capturing the transcriptional diversity of cells. A new wave of methods to integrate scRNAseq and biophysical measurements is facilitating the linkage of transcriptomic data to cellular function, which provides physiological insight into cellular states. I briefly discuss critical factors of these phenotypical characterizations such as timescales, information content, and analytical tools. Dedicated sections focus on the integration with cell morphology, calcium imaging, and electrophysiology (patch-seq), emphasizing their complementary roles. I discuss their application in elucidating cellular states, refining cell type classifications, and uncovering functional differences in cell subtypes. To illustrate the practical applications and benefits of these methods, I highlight their use in tissues with excitable cell-types such as the brain, pancreatic islets, and the retina. The potential of combining functional phenotyping with spatial transcriptomics for a detailed mapping of cell phenotypes in situ is explored. Finally, I discuss open questions and future perspectives, emphasizing the need for a shift towards broader accessibility through increased throughput.

Identification of type 2 diabetes- and obesity-associated human β-cells using deep transfer learning

Diabetes affects >10% of adults worldwide and is caused by impaired production or response to insulin, resulting in chronic hyperglycemia. Pancreatic islet β-cells are the sole source of endogenous insulin and our understanding of β-cell dysfunction and death in type 2 diabetes (T2D) is incomplete. Single-cell RNA-seq data supports heterogeneity as an important factor in β-cell function and survival. However, it is difficult to identify which β-cell phenotypes are critical for T2D etiology and progression. Our goal was to prioritize specific disease-related β-cell subpopulations to better understand T2D pathogenesis and identify relevant genes for targeted therapeutics. To address this, we applied a deep transfer learning tool, DEGAS, which maps disease associations onto single-cell RNA-seq data from bulk expression data. Independent runs of DEGAS using T2D or obesity status identified distinct β-cell subpopulations. A singular cluster of T2D-associated β-cells was identified; however, β-cells with high obese-DEGAS scores contained two subpopulations derived largely from either non-diabetic or T2D donors. The obesity-associated non-diabetic cells were enriched for translation and unfolded protein response genes compared to T2D cells. We selected DLK1 for validation by immunostaining in human pancreas sections from healthy and T2D donors. DLK1 was heterogeneously expressed among β-cells and appeared depleted from T2D islets. In conclusion, DEGAS has the potential to advance our holistic understanding of the β-cell transcriptomic phenotypes, including features that distinguish β-cells in obese non-diabetic or lean T2D states. Future work will expand this approach to additional human islet omics datasets to reveal the complex multicellular interactions driving T2D.

PancrESS - a meta-analysis resource for understanding cell-type specific expression in the human pancreas

BACKGROUND

The human pancreas is composed of specialized cell types producing hormones and enzymes critical to human health. These specialized functions are the result of cell type-specific transcriptional programs which manifest in cell-specific gene expression. Understanding these programs is essential to developing therapies for pancreatic disorders. Transcription in the human pancreas has been widely studied by single-cell RNA technologies, however the diversity of protocols and analysis methods hinders their interpretability in the aggregate.

RESULTS

In this work, we perform a meta-analysis of pancreatic single-cell RNA sequencing data. We present a database for reference transcriptome abundances and cell-type specificity metrics. This database facilitates the identification and definition of marker genes within the pancreas. Additionally, we introduce a versatile tool which is freely available as an R package, and should permit integration into existing workflows. Our tool accepts count data files generated by widely-used single-cell gene expression platforms in their original format, eliminating an additional pre-formatting step. Although we designed it to calculate expression specificity of pancreas cell types, our tool is agnostic to the biological source of count data, extending its applicability to other biological systems.

CONCLUSIONS

Our findings enhance the current understanding of expression specificity within the pancreas, surpassing previous work in terms of scope and detail. Furthermore, our database and tool enable researchers to perform similar calculations in diverse biological systems, expanding the applicability of marker gene identification and facilitating comparative analyses.

Patch-seq: Advances and Biological Applications

Multimodal analysis of gene-expression patterns, electrophysiological properties, and morphological phenotypes at the single-cell/single-nucleus level has been arduous because of the diversity and complexity of neurons. The emergence of Patch-sequencing (Patch-seq) directly links transcriptomics, morphology, and electrophysiology, taking neuroscience research to a multimodal era. In this review, we summarized the development of Patch-seq and recent applications in the cortex, hippocampus, and other nervous systems. Through generating multimodal cell type atlases, targeting specific cell populations, and correlating transcriptomic data with phenotypic information, Patch-seq has provided new insight into outstanding questions in neuroscience. We highlight the challenges and opportunities of Patch-seq in neuroscience and hope to shed new light on future neuroscience research.

Genetic risk converges on regulatory networks mediating early type 2 diabetes

Type 2 diabetes mellitus (T2D), a major cause of worldwide morbidity and mortality, is characterized by dysfunction of insulin-producing pancreatic islet β cells1,2. T2D genome-wide association studies (GWAS) have identified hundreds of signals in non-coding and β cell regulatory genomic regions, but deciphering their biological mechanisms remains challenging3-5. Here, to identify early disease-driving events, we performed traditional and multiplexed pancreatic tissue imaging, sorted-islet cell transcriptomics and islet functional analysis of early-stage T2D and control donors. By integrating diverse modalities, we show that early-stage T2D is characterized by β cell-intrinsic defects that can be proportioned into gene regulatory modules with enrichment in signals of genetic risk. After identifying the β cell hub gene and transcription factor RFX6 within one such module, we demonstrated multiple layers of genetic risk that converge on an RFX6-mediated network to reduce insulin secretion by β cells. RFX6 perturbation in primary human islet cells alters β cell chromatin architecture at regions enriched for T2D GWAS signals, and population-scale genetic analyses causally link genetically predicted reduced RFX6 expression with increased T2D risk. Understanding the molecular mechanisms of complex, systemic diseases necessitates integration of signals from multiple molecules, cells, organs and individuals, and thus we anticipate that this approach will be a useful template to identify and validate key regulatory networks and master hub genes for other diseases or traits using GWAS data.

Differential CpG methylation at Nnat in the early establishment of beta cell heterogeneity

Aims/hypothesis

Beta cells within the pancreatic islet represent a heterogenous population wherein individual sub-groups of cells make distinct contributions to the overall control of insulin secretion. These include a subpopulation of highly-connected 'hub' cells, important for the propagation of intercellular Ca2+ waves. Functional subpopulations have also been demonstrated in human beta cells, with an altered subtype distribution apparent in type 2 diabetes. At present, the molecular mechanisms through which beta cell hierarchy is established are poorly understood. Changes at the level of the epigenome provide one such possibility which we explore here by focussing on the imprinted gene neuronatin (Nnat), which is required for normal insulin synthesis and secretion.

Methods

Single cell RNA-seq datasets were examined using Seurat 4.0 and ClusterProfiler running under R. Transgenic mice expressing eGFP under the control of the Nnat enhancer/promoter regions were generated for fluorescence-activated cell (FAC) sorting of beta cells and downstream analysis of CpG methylation by bisulphite and RNA sequencing, respectively. Animals deleted for the de novo methyltransferase, DNMT3A from the pancreatic progenitor stage were used to explore control of promoter methylation. Proteomics was performed using affinity purification mass spectrometry and Ca2+ dynamics explored by rapid confocal imaging of Cal-520 and Cal-590. Insulin secretion was measured using Homogeneous Time Resolved Fluorescence Imaging.

Results

Nnat mRNA was differentially expressed in a discrete beta cell population in a developmental stage- and DNA methylation (DNMT3A)-dependent manner. Thus, pseudo-time analysis of embryonic data sets demonstrated the early establishment of Nnat-positive and negative subpopulations during embryogenesis. NNAT expression is also restricted to a subset of beta cells across the human islet that is maintained throughout adult life. NNAT+ beta cells also displayed a discrete transcriptome at adult stages, representing a sub-population specialised for insulin production, reminiscent of recently-described "βHI" cells and were diminished in db/db mice. 'Hub' cells were less abundant in the NNAT+ population, consistent with epigenetic control of this functional specialization.

Conclusions/interpretation

These findings demonstrate that differential DNA methylation at Nnat represents a novel means through which beta cell heterogeneity is established during development. We therefore hypothesise that changes in methylation at this locus may thus contribute to a loss of beta cell hierarchy and connectivity, potentially contributing to defective insulin secretion in some forms of diabetes.

VTA dopamine neurons are hyperexcitable in 3xTg-AD mice due to casein kinase 2-dependent SK channel dysfunction

Alzheimer's disease (AD) patients exhibit neuropsychiatric symptoms that extend beyond classical cognitive deficits, suggesting involvement of subcortical areas. Here, we investigated the role of midbrain dopamine (DA) neurons in AD using the amyloid + tau-driven 3xTg-AD mouse model. We found deficits in reward-based operant learning in AD mice, suggesting possible VTA DA neuron dysregulation. Physiological assessment revealed hyperexcitability and disrupted firing in DA neurons caused by reduced activity of small-conductance calcium-activated potassium (SK) channels. RNA sequencing from contents of single patch-clamped DA neurons (Patch-seq) identified up-regulation of the SK channel modulator casein kinase 2 (CK2). Pharmacological inhibition of CK2 restored SK channel activity and normal firing patterns in 3xTg-AD mice. These findings shed light on a complex interplay between neuropsychiatric symptoms and subcortical circuits in AD, paving the way for novel treatment strategies.

Dynamic scRNA-seq of live human pancreatic slices reveals functional endocrine cell neogenesis through an intermediate ducto-acinar stage

Human pancreatic plasticity is implied from multiple single-cell RNA sequencing (scRNA-seq) studies. However, these have been invariably based on static datasets from which fate trajectories can only be inferred using pseudotemporal estimations. Furthermore, the analysis of isolated islets has resulted in a drastic underrepresentation of other cell types, hindering our ability to interrogate exocrine-endocrine interactions. The long-term culture of human pancreatic slices (HPSs) has presented the field with an opportunity to dynamically track tissue plasticity at the single-cell level. Combining datasets from same-donor HPSs at different time points, with or without a known regenerative stimulus (BMP signaling), led to integrated single-cell datasets storing true temporal or treatment-dependent information. This integration revealed population shifts consistent with ductal progenitor activation, blurring of ductal/acinar boundaries, formation of ducto-acinar-endocrine differentiation axes, and detection of transitional insulin-producing cells. This study provides the first longitudinal scRNA-seq analysis of whole human pancreatic tissue, confirming its plasticity in a dynamic fashion.

A link between energy metabolism and plant host adaptation states in the two-spotted spider mite, Tetranychus urticae (Koch)

Energy metabolism is a highly conserved process that balances generation of cellular energy and maintenance of redox homeostasis. It consists of five interconnected pathways: glycolysis, tricarboxylic acid cycle, pentose phosphate, trans-sulfuration, and NAD+ biosynthesis pathways. Environmental stress rewires cellular energy metabolism. Type-2 diabetes is a well-studied energy metabolism rewiring state in human pancreatic β-cells where glucose metabolism is uncoupled from insulin secretion. The two-spotted spider mite, Tetranychus urticae (Koch), exhibits a remarkable ability to adapt to environmental stress. Upon transfer to unfavourable plant hosts, mites experience extreme xenobiotic stress that dramatically affects their survivorship and fecundity. However, within 25 generations, mites adapt to the xenobiotic stress and restore their fitness. Mites' ability to withstand long-term xenobiotic stress raises a question of their energy metabolism states during host adaptation. Here, we compared the transcriptional responses of five energy metabolism pathways between host-adapted and non-adapted mites while using responses in human pancreatic islet donors to model these pathways under stress. We found that non-adapted mites and human pancreatic β-cells responded in a similar manner to host plant transfer and diabetogenic stress respectively, where redox homeostasis maintenance was favoured over energy generation. Remarkably, we found that upon host-adaptation, mite energy metabolic states were restored to normal. These findings suggest that genes involved in energy metabolism can serve as molecular markers for mite host-adaptation.

An Integrated Map of Cell Type-Specific Gene Expression in Pancreatic Islets

Pancreatic islets consist of multiple cell types that produce hormones required for glucose homeostasis, and islet dysfunction is a major factor in type 1 and type 2 diabetes. Numerous studies have assessed transcription across individual cell types using single-cell assays; however, there is no canonical reference of gene expression in islet cell types that is also easily accessible for researchers to query and use in bioinformatics pipelines. Here we present an integrated map of islet cell type-specific gene expression from 192,203 cells from single-cell RNA sequencing of 65 donors without diabetes, donors who were type 1 diabetes autoantibody positive, donors with type 1 diabetes, and donors with type 2 diabetes from the Human Pancreas Analysis Program. We identified 10 distinct cell types, annotated subpopulations of several cell types, and defined cell type-specific marker genes. We tested differential expression within each cell type across disease states and identified 1,701 genes with significant changes in expression, with most changes observed in β-cells from donors with type 1 diabetes. To facilitate user interaction, we provide several single-cell visualization and reference mapping tools, as well as the open-access analytical pipelines used to create this reference. The results will serve as a valuable resource to investigators studying islet biology.

ARTICLE HIGHLIGHTS

Endoplasmic reticulum stress in pancreatic β cells induces incretin desensitization and β-cell dysfunction via ATF4-mediated PDE4D expression

Pancreatic β-cell dysfunction and eventual loss are key steps in the progression of type 2 diabetes (T2D). Endoplasmic reticulum (ER) stress responses, especially those mediated by the protein kinase RNA-like ER kinase and activating transcription factor 4 (PERK-ATF4) pathway, have been implicated in promoting these β-cell pathologies. However, the exact molecular events surrounding the role of the PERK-ATF4 pathway in β-cell dysfunction remain unknown. Here, we report our discovery that ATF4 promotes the expression of PDE4D, which disrupts β-cell function via a downregulation of cAMP signaling. We found that β-cell-specific transgenic expression of ATF4 led to early β-cell dysfunction and loss, a phenotype that resembles accelerated T2D. Expression of ATF4, rather than C/EBP homologous protein (CHOP), promoted PDE4D expression, reduced cAMP signaling, and attenuated responses to incretins and elevated glucose. Furthermore, we found that β-cells of leptin receptor-deficient diabetic (db/db) mice had elevated nuclear localization of ATF4 and PDE4D expression, accompanied by impaired β-cell function. Accordingly, pharmacological inhibition of the ATF4 pathway attenuated PDE4D expression in the islets and promoted incretin-simulated glucose tolerance and insulin secretion in db/db mice. Finally, we found that inhibiting PDE4 activity with selective pharmacological inhibitors improved β-cell function in both db/db mice and β-cell-specific ATF4 transgenic mice. In summary, our results indicate that ER stress causes β-cell failure via ATF4-mediated PDE4D production, suggesting the ATF4-PDE4D pathway could be a therapeutic target for protecting β-cell function during the progression of T2D.NEW & NOTEWORTHY Endoplasmic reticulum stress has been implied to cause multiple β-cell pathologies during the progression of type 2 diabetes (T2D). However, the precise molecular events underlying this remain unknown. Here, we discovered that elevated ATF4 activity, which was seen in T2D β cells, attenuated β-cell proliferation and impaired insulin secretion via PDE4D-mediated downregulation of cAMP signaling. Additionally, we demonstrated that pharmacological inhibition of the ATF4 pathway or PDE4D activity alleviated β-cell dysfunction, suggesting its therapeutic usefulness against T2D.

An integrative single-cell multi-omics profiling of human pancreatic islets identifies T1D associated genes and regulatory signals

Genome wide association studies (GWAS) have identified over 100 signals associated with type 1 diabetes (T1D). However, translating any given T1D GWAS signal into mechanistic insights, including putative causal variants and the context (cell type and cell state) in which they function, has been limited. Here, we present a comprehensive multi-omic integrative analysis of single-cell/nucleus resolution profiles of gene expression and chromatin accessibility in healthy and autoantibody+ (AAB+) human islets, as well as islets under multiple T1D stimulatory conditions. We broadly nominate effector cell types for all T1D GWAS signals. We further nominated higher-resolution contexts, including effector cell types, regulatory elements, and genes for three independent T1D risk variants acting through islet cells within the pancreas at the DLK1/MEG3, RASGRP1, and TOX loci. Subsequently, we created isogenic gene knockouts DLK1-/-, RASGRP1-/-, and TOX-/-, and the corresponding regulatory region knockout, RASGRP1Δ, and DLK1Δ hESCs. Loss of RASGRP1 or DLK1, as well as knockout of the regulatory region of RASGRP1 or DLK1, increased β cell apoptosis. Additionally, pancreatic β cells derived from isogenic hESCs carrying the risk allele of rs3783355A/A exhibited increased β cell death. Finally, RNA-seq and ATAC-seq identified five genes upregulated in both RASGRP1-/- and DLK1-/- β-like cells, four of which are associated with T1D. Together, this work reports an integrative approach for combining single cell multi-omics, GWAS, and isogenic hESC-derived β-like cells to prioritize the T1D associated signals and their underlying context-specific cell types, genes, SNPs, and regulatory elements, to illuminate biological functions and molecular mechanisms.

Hyperglucagonaemia in diabetes: altered amino acid metabolism triggers mTORC1 activation, which drives glucagon production

AIM/HYPOTHESIS

Hyperglycaemia is associated with alpha cell dysfunction, leading to dysregulated glucagon secretion in type 1 and type 2 diabetes; however, the mechanisms involved are still elusive. The nutrient sensor mammalian target of rapamycin complex 1 (mTORC1) plays a major role in the maintenance of alpha cell mass and function. We studied the regulation of alpha cell mTORC1 by nutrients and its role in the development of hyperglucagonaemia in diabetes.

METHODS

Alpha cell mTORC1 activity was assessed by immunostaining for phosphorylation of its downstream target, the ribosomal protein S6, and glucagon, followed by confocal microscopy on pancreatic sections and flow cytometry on dispersed human and mouse islets and the alpha cell line, αTC1-6. Metabolomics and metabolic flux were studied by 13C glucose labelling in 2.8 or 16.7 mmol/l glucose followed by LC-MS analysis. To study the role of mTORC1 in mediating hyperglucagonaemia in diabetes, we generated an inducible alpha cell-specific Rptor knockout in the Akita mouse model of diabetes and tested the effects on glucose tolerance by IPGTT and on glucagon secretion.

RESULTS

mTORC1 activity was increased in alpha cells from diabetic Akita mice in parallel to the development of hyperglycaemia and hyperglucagonaemia (two- to eightfold increase). Acute exposure of mouse and human islets to amino acids stimulated alpha cell mTORC1 (3.5-fold increase), whereas high glucose concentrations inhibited mTORC1 (1.4-fold decrease). The mTORC1 response to glucose was abolished in human and mouse diabetic alpha cells following prolonged islet exposure to high glucose levels, resulting in sustained activation of mTORC1, along with increased glucagon secretion. Metabolomics and metabolic flux analysis showed that exposure to high glucose levels enhanced glycolysis, glucose oxidation and the synthesis of glucose-derived amino acids. In addition, chronic exposure to high glucose levels increased the expression of Slc7a2 and Slc38a4, which encode amino acid transporters, as well as the levels of branched-chain amino acids and methionine cycle metabolites (~1.3-fold increase for both). Finally, conditional Rptor knockout in alpha cells from adult diabetic mice inhibited mTORC1, thereby inhibiting glucagon secretion (~sixfold decrease) and improving diabetes, despite persistent insulin deficiency.

CONCLUSIONS/INTERPRETATION

Alpha cell exposure to hyperglycaemia enhances amino acid synthesis and transport, resulting in sustained activation of mTORC1, thereby increasing glucagon secretion. mTORC1 therefore plays a major role in mediating alpha cell dysfunction in diabetes.

DATA AVAILABILITY

All sequencing data are available from the Gene Expression Omnibus (GEO) repository (accession no. GSE154126; https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE154126 ).

Global proteomics reveals insulin abundance as a marker of human islet homeostasis alterations

AIM

The variation in quality between the human islet samples represents a major problem for research, especially when used as control material. The assays assessing the quality of human islets used in research are non-standardized and limited, with many important parameters not being consistently assessed. High-throughput studies aimed at characterizing the diversity and segregation markers among apparently functionally healthy islet preps are thus a requirement. Here, we designed a pilot study to comprehensively identify the diversity of global proteome signatures and the deviation from normal homeostasis in randomly selected human-isolated islet samples.

METHODS

By using Tandem Mass Tag 16-plex proteomics, we focused on the recurrently observed disparity in the detected insulin abundance between the samples, used it as a segregating parameter, and analyzed the correlated changes in the proteome signature and homeostasis by pathway analysis.

RESULTS

In this pilot study, we showed that insulin protein abundance is a predictor of human islet homeostasis and quality. This parameter is independent of other quality predictors within their acceptable range, thus being able to further stratify islets samples of apparent good quality. Human islets with low amounts of insulin displayed changes in their metabolic and signaling profile, especially in regard to energy homeostasis and cell identity maintenance. We further showed that xenotransplantation into diabetic hosts is not expected to improve the pre-transplantation signature, as it has a negative effect on energy balance, antioxidant activity, and islet cell identity.

CONCLUSIONS

Insulin protein abundance predicts significant changes in human islet homeostasis among random samples of apparently good quality.

The simplified Kirchhoff network model (SKNM): a cell-based reaction-diffusion model of excitable tissue

Cell-based models of excitable tissues offer the advantage of cell-level precision, which cannot be achieved using traditional homogenized electrophysiological models. However, this enhanced accuracy comes at the cost of increased computational demands, necessitating the development of efficient cell-based models. The widely-accepted bidomain model serves as the standard in computational cardiac electrophysiology, and under certain anisotropy ratio conditions, it is well known that it can be reduced to the simpler monodomain model. Recently, the Kirchhoff Network Model (KNM) was developed as a cell-based counterpart to the bidomain model. In this paper, we aim to demonstrate that KNM can be simplified using the same steps employed to derive the monodomain model from the bidomain model. We present the cell-based Simplified Kirchhoff Network Model (SKNM), which produces results closely aligned with those of KNM while requiring significantly less computational resources.

MRPS6 modulates glucose-stimulated insulin secretion in mouse islet cells through mitochondrial unfolded protein response

Lack of efficient insulin secretion from the pancreas can lead to impaired glucose tolerance (IGT), prediabetes, and diabetes. We have previously identified two IGT-associated single nucleotide polymorphisms (SNPs) rs62212118 and rs13052524 located at two overlapping genes: MRPS6 and SLC5A3. In this study, we show that MRPS6 but not SLC5A3 regulates glucose-stimulated insulin secretion (GSIS) in primary human β-cell and a mouse pancreatic insulinoma β-cell line. Data mining and biochemical studies reveal that MRPS6 is positively regulated by the mitochondrial unfolded protein response (UPRmt), but feedback inhibits UPRmt. Disruption of such feedback by MRPS6 knockdown causes UPRmt hyperactivation in high glucose conditions, hence elevated ROS levels, increased apoptosis, and impaired GSIS. Conversely, MRPS6 overexpression reduces UPRmt, mitigates high glucose-induced ROS levels and apoptosis, and enhances GSIS in an ATF5-dependent manner. Consistently, UPRmt up-regulation or down-regulation by modulating ATF5 expression is sufficient to decrease or increase GSIS. The negative role of UPRmt in GSIS is further supported by analysis of public transcriptomic data from murine islets. In all, our studies identify MRPS6 and UPRmt as novel modulators of GSIS and apoptosis in β-cells, contributing to our understanding of the molecular and cellular mechanisms of IGT, prediabetes, and diabetes.

NEUROD1 reinforces endocrine cell fate acquisition in pancreatic development

NEUROD1 is a transcription factor that helps maintain a mature phenotype of pancreatic β cells. Disruption of Neurod1 during pancreatic development causes severe neonatal diabetes; however, the exact role of NEUROD1 in the differentiation programs of endocrine cells is unknown. Here, we report a crucial role of the NEUROD1 regulatory network in endocrine lineage commitment and differentiation. Mechanistically, transcriptome and chromatin landscape analyses demonstrate that Neurod1 inactivation triggers a downregulation of endocrine differentiation transcription factors and upregulation of non-endocrine genes within the Neurod1-deficient endocrine cell population, disturbing endocrine identity acquisition. Neurod1 deficiency altered the H3K27me3 histone modification pattern in promoter regions of differentially expressed genes, which resulted in gene regulatory network changes in the differentiation pathway of endocrine cells, compromising endocrine cell potential, differentiation, and functional properties.

Loss of ATP-sensitive channel expression and function decreases opioid sensitivity in a mouse model of type 2 diabetes

During diabetes, β-cell dysfunction due to loss of potassium channels sensitive to ATP, known as KATP channels occurs progressively over time contributing to hyperglycemia. KATP channels are additionally present in the central and peripheral nervous systems and are downstream targets of opioid receptor signaling. The aim of this study is to investigate if KATP channel expression or activity in the nervous system changes in diabetic mice and if morphine antinociception changes in mice fed a high fat diet (HFD) for 16 weeks compared to controls. Mechanical thresholds were also monitored before and after administration of glyburide or nateglinide, KATP channel antagonists, for four weeks. HFD mice have decreased antinociception to systemic morphine, which is exacerbated after systemic treatment with glyburide or nateglinide. HFD mice also have lower rotarod scores, decreased mobility in an open field test, and lower burrowing behavior compared to their control diet counterparts, which is unaffected by KATP channel antagonist delivery. Expression of KATP channel subunits, Kcnj11 (Kir6.2) and Abcc8 (SUR1), were decreased in the peripheral and central nervous system in HFD mice, which is significantly correlated with baseline paw withdrawal thresholds. Upregulation of SUR1 through an adenovirus delivered intrathecally increased morphine antinociception in HFD mice, whereas Kir6.2 upregulation improved morphine antinociception only marginally. Perspective: This article presents the potential link between KATP channel function and neuropathy during diabetes. There is a need for increased knowledge in how diabetes affects structural and molecular changes in the nervous system to lead to the progression of chronic pain and sensory issues.

Single cell multiomic analysis reveals diabetes-associated β-cell heterogeneity driven by HNF1A

Broad heterogeneity in pancreatic β-cell function and morphology has been widely reported. However, determining which components of this cellular heterogeneity serve a diabetes-relevant function remains challenging. Here, we integrate single-cell transcriptome, single-nuclei chromatin accessibility, and cell-type specific 3D genome profiles from human islets and identify Type II Diabetes (T2D)-associated β-cell heterogeneity at both transcriptomic and epigenomic levels. We develop a computational method to explicitly dissect the intra-donor and inter-donor heterogeneity between single β-cells, which reflect distinct mechanisms of T2D pathogenesis. Integrative transcriptomic and epigenomic analysis identifies HNF1A as a principal driver of intra-donor heterogeneity between β-cells from the same donors; HNF1A expression is also reduced in β-cells from T2D donors. Interestingly, HNF1A activity in single β-cells is significantly associated with lower Na+ currents and we nominate a HNF1A target, FXYD2, as the primary mitigator. Our study demonstrates the value of investigating disease-associated single-cell heterogeneity and provides new insights into the pathogenesis of T2D.

Delineating mouse β-cell identity during lifetime and in diabetes with a single cell atlas

Although multiple pancreatic islet single-cell RNA-sequencing (scRNA-seq) datasets have been generated, a consensus on pancreatic cell states in development, homeostasis and diabetes as well as the value of preclinical animal models is missing. Here, we present an scRNA-seq cross-condition mouse islet atlas (MIA), a curated resource for interactive exploration and computational querying. We integrate over 300,000 cells from nine scRNA-seq datasets consisting of 56 samples, varying in age, sex and diabetes models, including an autoimmune type 1 diabetes model (NOD), a glucotoxicity/lipotoxicity type 2 diabetes model (db/db) and a chemical streptozotocin β-cell ablation model. The β-cell landscape of MIA reveals new cell states during disease progression and cross-publication differences between previously suggested marker genes. We show that β-cells in the streptozotocin model transcriptionally correlate with those in human type 2 diabetes and mouse db/db models, but are less similar to human type 1 diabetes and mouse NOD β-cells. We also report pathways that are shared between β-cells in immature, aged and diabetes models. MIA enables a comprehensive analysis of β-cell responses to different stressors, providing a roadmap for the understanding of β-cell plasticity, compensation and demise.

Islet cells in human type 1 diabetes: from recent advances to novel therapies - a symposium-based roadmap for future research

There is a growing understanding that the early phases of type 1 diabetes (T1D) are characterised by a deleterious dialogue between the pancreatic beta cells and the immune system. This, combined with the urgent need to better translate this growing knowledge into novel therapies, provided the background for the JDRF-DiabetesUK-INNODIA-nPOD symposium entitled 'Islet cells in human T1D: from recent advances to novel therapies', which took place in Stockholm, Sweden, in September 2022. We provide in this article an overview of the main themes addressed in the symposium, pointing to both promising conclusions and key unmet needs that remain to be addressed in order to achieve better approaches to prevent or reverse T1D.

Pancreatic β-cell senescence in diabetes: mechanisms, markers and therapies

Cellular senescence is a response to a wide variety of stressors, including DNA damage, oncogene activation and physiologic aging, and pathologically accelerated senescence contributes to human disease, including diabetes mellitus. Indeed, recent work in this field has demonstrated a role for pancreatic β-cell senescence in the pathogenesis of Type 1 Diabetes, Type 2 Diabetes and monogenic diabetes. Small molecule or genetic targeting of senescent β-cells has shown promise as a novel therapeutic approach for preventing and treating diabetes. Despite these advances, major questions remain around the molecular mechanisms driving senescence in the β-cell, identification of molecular markers that distinguish senescent from non-senescent β-cell subpopulations, and translation of proof-of-concept therapies into novel treatments for diabetes in humans. Here, we summarize the current state of the field of β-cell senescence, highlighting insights from mouse models as well as studies on human islets and β-cells. We identify markers that have been used to detect β-cell senescence to unify future research efforts in this field. We discuss emerging concepts of the natural history of senescence in β-cells, heterogeneity of senescent β-cells subpopulations, role of sex differences in senescent responses, and the consequences of senescence on integrated islet function and microenvironment. As a young and developing field, there remain many open research questions which need to be addressed to move senescence-targeted approaches towards clinical investigation.

Not the second fiddle: α cell development, identity, and function in health and diabetes

Historic and emerging studies provide evidence for the deterioration of pancreatic α cell function and identity in diabetes mellitus. Increased access to human tissue and the availability of more sophisticated molecular technologies have identified key insights into how α cell function and identity are preserved in healthy conditions and how they become dysfunctional in response to stress. These studies have revealed evidence of impaired glucagon secretion, shifts in α cell electrophysiology, changes in α cell mass, dysregulation of α cell transcription, and α-to-β cell conversion prior to and during diabetes. In this review, we outline the current state of research on α cell identity in health and disease. Evidence in model organisms and humans suggests that in addition to β cell dysfunction, diabetes is associated with a fundamental dysregulation of α cell identity. Importantly, epigenetic studies have revealed that α cells retain more poised and open chromatin at key cell-specific and diabetes-dysregulated genes, supporting the model that the inherent epigenetic plasticity of α cells makes them susceptible to the transcriptional changes that potentiate the loss of identity and function seen in diabetes. Thus, additional research into the maintenance of α cell identity and function is critical to fully understanding diabetes. Furthermore, these studies suggest α cells could represent an alternative source of new β cells for diabetes treatment.

Immune cell identifier and classifier (ImmunIC) for single cell transcriptomic readouts

Single cell RNA sequencing has a central role in immune profiling, identifying specific immune cells as disease markers and suggesting therapeutic target genes of immune cells. Immune cell-type annotation from single cell transcriptomics is in high demand for dissecting complex immune signatures from multicellular blood and organ samples. However, accurate cell type assignment from single-cell RNA sequencing data alone is complicated by a high level of gene expression heterogeneity. Many computational methods have been developed to respond to this challenge, but immune cell annotation accuracy is not highly desirable. We present ImmunIC, a simple and robust tool for immune cell identification and classification by combining marker genes with a machine learning method. With over two million immune cells and half-million non-immune cells from 66 single cell RNA sequencing studies, ImmunIC shows 98% accuracy in the identification of immune cells. ImmunIC outperforms existing immune cell classifiers, categorizing into ten immune cell types with 92% accuracy. We determine peripheral blood mononuclear cell compositions of severe COVID-19 cases and healthy controls using previously published single cell transcriptomic data, permitting the identification of immune cell-type specific differential pathways. Our publicly available tool can maximize the utility of single cell RNA profiling by functioning as a stand-alone bioinformatic cell sorter, advancing cell-type specific immune profiling for the discovery of disease-specific immune signatures and therapeutic targets.

Why does the immune system destroy pancreatic β-cells but not α-cells in type 1 diabetes?

A perplexing feature of type 1 diabetes (T1D) is that the immune system destroys pancreatic β-cells but not neighbouring α-cells, even though both β-cells and α-cells are dysfunctional. Dysfunction, however, progresses to death only for β-cells. Recent findings indicate important differences between these two cell types. First, expression of BCL2L1, a key antiapoptotic gene, is higher in α-cells than in β-cells. Second, endoplasmic reticulum (ER) stress-related genes are differentially expressed, with higher expression levels of pro-apoptotic CHOP in β-cells than in α-cells and higher expression levels of HSPA5 (which encodes the protective chaperone BiP) in α-cells than in β-cells. Third, expression of viral recognition and innate immune response genes is higher in α-cells than in β-cells, contributing to the enhanced resistance of α-cells to coxsackievirus infection. Fourth, expression of the immune-inhibitory HLA-E molecule is higher in α-cells than in β-cells. Of note, α-cells are less immunogenic than β-cells, and the CD8+ T cells invading the islets in T1D are reactive to pre-proinsulin but not to glucagon. We suggest that this finding is a result of the enhanced capacity of the α-cell to endure viral infections and ER stress, which enables them to better survive early stressors that can cause cell death and consequently amplify antigen presentation to the immune system. Moreover, the processing of the pre-proglucagon precursor in enteroendocrine cells might favour immune tolerance towards this potential self-antigen compared to pre-proinsulin.

The human α cell in health and disease

In commemoration of 100 years since the discovery of glucagon, we review current knowledge about the human α cell. Alpha cells make up 30-40% of human islet endocrine cells and play a major role in regulating whole-body glucose homeostasis, largely through the direct actions of their main secretory product - glucagon - on peripheral organs. Additionally, glucagon and other secretory products of α cells, namely acetylcholine, glutamate, and glucagon-like peptide-1, have been shown to play an indirect role in the modulation of glucose homeostasis through autocrine and paracrine interactions within the islet. Studies of glucagon's role as a counterregulatory hormone have revealed additional important functions of the α cell, including the regulation of multiple aspects of energy metabolism outside that of glucose. At the molecular level, human α cells are defined by the expression of conserved islet-enriched transcription factors and various enriched signature genes, many of which have currently unknown cellular functions. Despite these common threads, notable heterogeneity exists amongst human α cell gene expression and function. Even greater differences are noted at the inter-species level, underscoring the importance of further study of α cell physiology in the human context. Finally, studies on α cell morphology and function in type 1 and type 2 diabetes, as well as other forms of metabolic stress, reveal a key contribution of α cell dysfunction to dysregulated glucose homeostasis in disease pathogenesis, making targeting the α cell an important focus for improving treatment.

Integrating genetics with single-cell multiomic measurements across disease states identifies mechanisms of beta cell dysfunction in type 2 diabetes

Dysfunctional pancreatic islet beta cells are a hallmark of type 2 diabetes (T2D), but a comprehensive understanding of the underlying mechanisms, including gene dysregulation, is lacking. Here we integrate information from measurements of chromatin accessibility, gene expression and function in single beta cells with genetic association data to nominate disease-causal gene regulatory changes in T2D. Using machine learning on chromatin accessibility data from 34 nondiabetic, pre-T2D and T2D donors, we identify two transcriptionally and functionally distinct beta cell subtypes that undergo an abundance shift during T2D progression. Subtype-defining accessible chromatin is enriched for T2D risk variants, suggesting a causal contribution of subtype identity to T2D. Both beta cell subtypes exhibit activation of a stress-response transcriptional program and functional impairment in T2D, which is probably induced by the T2D-associated metabolic environment. Our findings demonstrate the power of multimodal single-cell measurements combined with machine learning for characterizing mechanisms of complex diseases.

Epigenetic dosage identifies two major and functionally distinct β cell subtypes

The mechanisms that specify and stabilize cell subtypes remain poorly understood. Here, we identify two major subtypes of pancreatic β cells based on histone mark heterogeneity (βHI and βLO). βHI cells exhibit ∼4-fold higher levels of H3K27me3, distinct chromatin organization and compaction, and a specific transcriptional pattern. βHI and βLO cells also differ in size, morphology, cytosolic and nuclear ultrastructure, epigenomes, cell surface marker expression, and function, and can be FACS separated into CD24+ and CD24- fractions. Functionally, βHI cells have increased mitochondrial mass, activity, and insulin secretion in vivo and ex vivo. Partial loss of function indicates that H3K27me3 dosage regulates βHI/βLO ratio in vivo, suggesting that control of β cell subtype identity and ratio is at least partially uncoupled. Both subtypes are conserved in humans, with βHI cells enriched in humans with type 2 diabetes. Thus, epigenetic dosage is a novel regulator of cell subtype specification and identifies two functionally distinct β cell subtypes.

Single-nucleus RNA sequencing of human pancreatic islets identifies novel gene sets and distinguishes β-cell subpopulations with dynamic transcriptome profiles

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) provides valuable insights into human islet cell types and their corresponding stable gene expression profiles. However, this approach requires cell dissociation that complicates its utility in vivo. On the other hand, single-nucleus RNA sequencing (snRNA-seq) has compatibility with frozen samples, elimination of dissociation-induced transcriptional stress responses, and affords enhanced information from intronic sequences that can be leveraged to identify pre-mRNA transcripts.

METHODS

We obtained nuclear preparations from fresh human islet cells and generated snRNA-seq datasets. We compared these datasets to scRNA-seq output obtained from human islet cells from the same donor. We employed snRNA-seq to obtain the transcriptomic profile of human islets engrafted in immunodeficient mice. In both analyses, we included the intronic reads in the snRNA-seq data with the GRCh38-2020-A library.

RESULTS

First, snRNA-seq analysis shows that the top four differentially and selectively expressed genes in human islet endocrine cells in vitro and in vivo are not the canonical genes but a new set of non-canonical gene markers including ZNF385D, TRPM3, LRFN2, PLUT (β-cells); PTPRT, FAP, PDK4, LOXL4 (α-cells); LRFN5, ADARB2, ERBB4, KCNT2 (δ-cells); and CACNA2D3, THSD7A, CNTNAP5, RBFOX3 (γ-cells). Second, by integrating information from scRNA-seq and snRNA-seq of human islet cells, we distinguish three β-cell sub-clusters: an INS pre-mRNA cluster (β3), an intermediate INS mRNA cluster (β2), and an INS mRNA-rich cluster (β1). These display distinct gene expression patterns representing different biological dynamic states both in vitro and in vivo. Interestingly, the INS mRNA-rich cluster (β1) becomes the predominant sub-cluster in vivo.

CONCLUSIONS

In summary, snRNA-seq and pre-mRNA analysis of human islet cells can accurately identify human islet cell populations, subpopulations, and their dynamic transcriptome profile in vivo.