Pseudotime Ordering of Single Human β-Cells Reveals States of Insulin Production and Unfolded Protein Response

Single Nucleus Multiome Analysis Reveals Early Inflammatory Response to High-Fat Diet in Mouse Pancreatic Islets

In periods of sustained hyper-nutrition, pancreatic β-cells undergo functional compensation through transcriptional upregulation of gene programs driving insulin secretion. This adaptation is essential for maintaining systemic glucose homeostasis and metabolic health. Using single nuclei multiomics, we have mapped the early transcriptional compensation mechanisms in murine islets of Langerhans exposed to high-fat diet (HFD) for one and three weeks. We show that β-cells exhibit the largest transcriptional response to HFD, characterized by early activation of proinflammatory eRegulons and downregulation of β-cell identity genes, particularly in a distinct subset of β-cells. Our observations translate to humans, as we observe an increase in the inflammatory gene signatures in human β-cells in pre-diabetes and diabetes. Collectively, these observations point to cellular cross-talk through proinflammatory signaling as a central and early driver of β-cell dysfunction that limits the compensatory capacity of β-cells, which is closely linked to the development of diabetes.

Signal transduction pathways controlling Ins2 gene activity and beta cell state transitions

Pancreatic β cells exist in low and high insulin gene activity states that are dynamic on a scale of hours to days. Here, we used live 3D imaging, mass spectrometry proteomics, and targeted perturbations of β cell signaling to comprehensively investigate Ins2(GFP)HIGH and Ins2(GFP)LOW β cell states. We identified the two Ins2 gene activity states in intact isolated islets and showed that cells in the same state were more likely to be nearer to each other. We report the proteomes of pure β cells to a depth of 5555 proteins and show that β cells with high Ins2 gene activity had reduced β cell immaturity factors, as well as increased translation. We identified activators of cAMP signaling (GLP1, IBMX) as powerful drivers of Ins2(GFP)LOW to Ins2(GFP)HIGH transitions. Okadaic acid and cyclosporine A had the opposite effects. This study provides new insight into the proteomic profiles and regulation of β cell states.

Islet single-cell transcriptomic profiling during obesity-induced beta cell expansion in female mice

Targeting beta cell proliferation is an appealing approach to restore glucose control in type 1 diabetes. However, the underlying mechanisms of beta cell proliferation remain incompletely understood, limiting identification of new therapeutic targets. Obesity is a naturally occurring process that potently induces human and rodent beta cell replication, representing an ideal model to study mechanisms of beta cell proliferation. We showed previously acute whole-body Lepr gene deletion in adult mice induces obesity and massive beta cell expansion. Here, using single-cell transcriptomics with female Lepr KO islets, we identified distinct populations of beta cells undergoing unfolded protein response (UPR), stress resolution, and cell cycle progression. Lepr KO beta cells undergoing UPR markedly increased chaperone protein, ribosomal biogenesis, and cell cycle transcriptional programs that were enriched for Xbp1 and Myc target genes. Our findings suggest a coordinated transcriptional mechanism involving Xbp1 and Myc to alleviate UPR and stimulate beta cell proliferation in obese female mice.

scCobra allows contrastive cell embedding learning with domain adaptation for single cell data integration and harmonization

The rapid advancement of single-cell technologies has created an urgent need for effective methods to integrate and harmonize single-cell data. Technical and biological variations across studies complicate data integration, while conventional tools often struggle with reliance on gene expression distribution assumptions and over-correction. Here, we present scCobra, a deep generative neural network designed to overcome these challenges through contrastive learning with domain adaptation. scCobra effectively mitigates batch effects, minimizes over-correction, and ensures biologically meaningful data integration without assuming specific gene expression distributions. It enables online label transfer across datasets with batch effects, allowing continuous integration of new data without retraining. Additionally, scCobra supports batch effect simulation, advanced multi-omic integration, and scalable processing of large datasets. By integrating and harmonizing datasets from similar studies, scCobra expands the available data for investigating specific biological problems, improving cross-study comparability, and revealing insights that may be obscured in isolated datasets.

SPOCK2 controls the proliferation and function of immature pancreatic β-cells through MMP2

Human pluripotent stem cell-derived β-cells (SC-β-cells) represent an alternative cell source for transplantation in diabetic patients. Although mitogens could in theory be used to expand β-cells, adult β-cells very rarely replicate. In contrast, newly formed β-cells, including SC-β-cells, display higher proliferative capacity and distinct transcriptional and functional profiles. Through bidirectional expression modulation and single-cell RNA-seq, we identified SPOCK2, an ECM protein, as an inhibitor of immature β-cell proliferation. Human β-cells lacking SPOCK2 presented elevated MMP2 expression and activity, leading to β-integrin-FAK-c-JUN pathway activation. Treatment with the MMP2 protein resulted in pronounced short- and long-term SC-β-cell expansion, significantly increasing glucose-stimulated insulin secretion in vitro and in vivo. These findings suggest that SPOCK2 mediates fetal β-cell proliferation and maturation. In summary, we identified a molecular mechanism that specifically regulates SC-β-cell proliferation and function, highlighting a unique signaling milieu of SC-β-cells with promise for the robust derivation of fully functional cells for transplantation.

Harnessing beta cell regeneration biology for diabetes therapy

The pandemic scale of diabetes mellitus is alarming, its complications remain devastating, and current treatments still pose a major burden on those affected and on the healthcare system as a whole. As the disease emanates from the destruction or dysfunction of insulin-producing pancreatic β-cells, a real cure requires their restoration and protection. An attractive strategy is to regenerate β-cells directly within the pancreas; however, while several approaches for β-cell regeneration have been proposed in the past, clinical translation has proven challenging. This review scrutinizes recent findings in β-cell regeneration and discusses their potential clinical implementation. Hereby, we aim to delineate a path for innovative, targeted therapies to help shift from 'caring for' to 'curing' diabetes.

Vemurafenib inhibits the replication of diabetogenic enteroviruses in intestinal epithelial and pancreatic beta cells

Enteroviruses, which infect via the gut, have been implicated in type 1 diabetes (T1D) development. Prolonged faecal shedding of enterovirus has been associated with islet autoimmunity. Additionally, enteroviral proteins and viral RNA have been detected in the pancreatic islets of individuals with recent-onset T1D, implicating their possible role in beta cell destruction. Despite this, no approved antiviral drugs currently exist that specifically target enterovirus infections for utilisation in disease interventions. Drug repurposing allows for the discovery of new clinical uses for existing drugs and can expedite drug discovery. Previously, the cancer drug Vemurafenib demonstrated unprecedented antiviral activity against several enteroviruses. In the present study, we assessed the efficacy of Vemurafenib and an analogue thereof in preventing infection or reducing the replication of enteroviruses associated with T1D. We tested Vemurafenib in intestinal epithelial cells (IECs) and insulin-producing beta cells. Additionally, we established a protocol for infecting human stem cell-derived islets (SC-islets) and used Vemurafenib and its analogue in this model. Our studies revealed that Vemurafenib exhibited strong antiviral properties in IECs and a beta cell line. The antiviral effect was also seen with the Vemurafenib analogue. SC-islets expressed the viral receptors CAR and DAF, with their highest expression in insulin- and glucagon-positive cells, respectively. SC-islets were successfully infected by CVBs and the antiviral activity of Vemurafenib and its analogue was confirmed in most SC-islet batches. In summary, our observations suggest that Vemurafenib and its analogue warrant further exploration as potential antiviral agents for the treatment of enterovirus-induced diseases, including T1D.

Heterogeneous enhancer states orchestrate β cell responses to metabolic stress

Obesity-induced β cell dysfunction contributes to the onset of type 2 diabetes. Nevertheless, elucidating epigenetic mechanisms underlying islet dysfunction at single cell level remains challenging. Here we profile single-nuclei RNA along with enhancer marks H3K4me1 or H3K27ac in islets from lean or obese mice. Our study identifies distinct gene signatures and enhancer states correlating with β cell dysfunction trajectory. Intriguingly, while many metabolic stress-induced genes exhibit concordant changes in both H3K4me1 and H3K27ac at their enhancers, expression changes of specific subsets are solely attributable to either H3K4me1 or H3K27ac dynamics. Remarkably, a subset of H3K4me1+H3K27ac- primed enhancers prevalent in lean β cells and occupied by FoxA2 are largely absent after metabolic stress. Lastly, cell-cell communication analysis identified the nerve growth factor (NGF) as protective paracrine signaling for β cells through repressing ER stress. In summary, our findings define the heterogeneous enhancer responses to metabolic challenges in individual β cells.

Analysis of the heterogeneity and complexity of murine extraorbital lacrimal gland via single-cell RNA sequencing

PURPOSE

The lacrimal gland is essential for maintaining ocular surface health and avoiding external damage by secreting an aqueous layer of the tear film. However, a healthy lacrimal gland's inventory of cell types and heterogeneity remains understudied.

METHODS

Here, 10X Genome-based single-cell RNA sequencing was used to generate an unbiased classification of cellular diversity in the extraorbital lacrimal gland (ELG) of C57BL/6J mice. From 43,850 high-quality cells, we produced an atlas of cell heterogeneity and defined cell types using classic marker genes. The possible functions of these cells were analyzed through bioinformatics analysis. Additionally, the CellChat was employed for a preliminary analysis of the cell-cell communication network in the ELG.

RESULTS

Over 37 subclasses of cells were identified, including seven types of glandular epithelial cells, three types of fibroblasts, ten types of myeloid-derived immune cells, at least eleven types of lymphoid-derived immune cells, and five types of vascular-associated cell subsets. The cell-cell communication network analysis revealed that fibroblasts and immune cells play a pivotal role in the dense intercellular communication network within the mouse ELG.

CONCLUSIONS

This study provides a comprehensive transcriptome atlas and related database of the mouse ELG.

Disruption of perinatal myeloid niches impacts the aging clock of pancreatic β cells

Perinatal expansion of pancreatic β cells is critical to metabolic adaptation. Yet, mechanisms surveying the fidelity by which proliferative events generate functional β cell pools remain unknown. We have previously identified a CCR2+ myeloid niche required for peri-natal β cell replication, with β cells dynamically responding to loss and repopulation of these myeloid cells with growth arrest and rebound expansion, respectively. Here, using a timed single-cell RNA-sequencing approach, we show that transient disruption of perinatal CCR2+ macrophages change islet β cell repertoires in young mice to resemble those of aged mice. Gene expression profiling and functional assays disclose prominent mitochondrial defects in β cells coupled to impaired redox states, NAD depletion, and DNA damage, leading to accelerated islets' dysfunction with age. These findings reveal an unexpected vulnerability of mitochondrial β cells' bioenergetics to the disruption of perinatal CCR2+ macrophages, implicating these cells in surveying early in life both the size and energy homeostasis of β cells populations.

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.

Identifying cell states in single-cell RNA-seq data at statistically maximal resolution

Single-cell RNA sequencing (scRNA-seq) has become a popular experimental method to study variation of gene expression within a population of cells. However, obtaining an accurate picture of the diversity of distinct gene expression states that are present in a given dataset is highly challenging because of the sparsity of the scRNA-seq data and its inhomogeneous measurement noise properties. Although a vast number of different methods is applied in the literature for clustering cells into subsets with 'similar' expression profiles, these methods generally lack rigorously specified objectives, involve multiple complex layers of normalization, filtering, feature selection, dimensionality-reduction, employ ad hoc measures of distance or similarity between cells, often ignore the known measurement noise properties of scRNA-seq measurements, and include a large number of tunable parameters. Consequently, it is virtually impossible to assign concrete biophysical meaning to the clusterings that result from these methods. Here we address the following problem: Given raw unique molecule identifier (UMI) counts of an scRNA-seq dataset, partition the cells into subsets such that the gene expression states of the cells in each subset are statistically indistinguishable, and each subset corresponds to a distinct gene expression state. That is, we aim to partition cells so as to maximally reduce the complexity of the dataset without removing any of its meaningful structure. We show that, given the known measurement noise structure of scRNA-seq data, this problem is mathematically well-defined and derive its unique solution from first principles. We have implemented this solution in a tool called Cellstates which operates directly on the raw data and automatically determines the optimal partition and cluster number, with zero tunable parameters. We show that, on synthetic datasets, Cellstates almost perfectly recovers optimal partitions. On real data, Cellstates robustly identifies subtle substructure within groups of cells that are traditionally annotated as a common cell type. Moreover, we show that the diversity of gene expression states that Cellstates identifies systematically depends on the tissue of origin and not on technical features of the experiments such as the total number of cells and total UMI count per cell. In addition to the Cellstates tool we also provide a small toolbox of software to place the identified cellstates into a hierarchical tree of higher-order clusters, to identify the most important differentially expressed genes at each branch of this hierarchy, and to visualize these results.

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.

Disrupted Endoplasmic Reticulum Ca2+ Handling: A Harβinger of β-Cell Failure

The β-cell workload increases in the setting of insulin resistance and reduced β-cell mass, which occurs in type 2 and type 1 diabetes, respectively. The prolonged elevation of insulin production and secretion during the pathogenesis of diabetes results in β-cell ER stress. The depletion of β-cell Ca2+ER during ER stress activates the unfolded protein response, leading to β-cell dysfunction. Ca2+ER is involved in many pathways that are critical to β-cell function, such as protein processing, tuning organelle and cytosolic Ca2+ handling, and modulating lipid homeostasis. Mutations that promote β-cell ER stress and deplete Ca2+ER stores are associated with or cause diabetes (e.g., mutations in ryanodine receptors and insulin). Thus, improving β-cell Ca2+ER handling and reducing ER stress under diabetogenic conditions could preserve β-cell function and delay or prevent the onset of diabetes. This review focuses on how mechanisms that control β-cell Ca2+ER are perturbed during the pathogenesis of diabetes and contribute to β-cell failure.

Endoplasmic reticulum stress in pancreatic β-cell dysfunction: The potential therapeutic role of dietary flavonoids

Diabetes mellitus (DM) is a global health burden that is characterized by the loss or dysfunction of pancreatic β-cells. In pancreatic β-cells, endoplasmic reticulum (ER) stress is a fact of life that contributes to β-cell loss or dysfunction. Despite recent advances in research, the existing treatment approaches such as lifestyle modification and use of conventional therapeutics could not prevent the loss or dysfunction of pancreatic β-cells to abrogate the disease progression. Therefore, targeting ER stress and the consequent unfolded protein response (UPR) in pancreatic β-cells may be a potential therapeutic strategy for diabetes treatment. Dietary phytochemicals have therapeutic applications in human health owing to their broad spectrum of biochemical and pharmacological activities. Flavonoids, which are commonly obtained from fruits and vegetables worldwide, have shown promising prospects in alleviating ER stress. Dietary flavonoids including quercetin, kaempferol, myricetin, isorhamnetin, fisetin, icariin, apigenin, apigetrin, vitexin, baicalein, baicalin, nobiletin hesperidin, naringenin, epigallocatechin 3-O-gallate hesperidin (EGCG), tectorigenin, liquiritigenin, and acacetin have shown inhibitory effects on ER stress in pancreatic β-cells. Dietary flavonoids modulate ER stress signaling components, chaperone proteins, transcription factors, oxidative stress, autophagy, apoptosis, and inflammatory responses to exert their pharmacological effects on pancreatic β-cells ER stress. This review focuses on the role of dietary flavonoids as potential therapeutic adjuvants in preserving pancreatic β-cells from ER stress. Highlights of the underlying mechanisms of action are also presented as well as possible strategies for clinical translation in the management of DM.

NKX2-2 based nuclei sorting on frozen human archival pancreas enables the enrichment of islet endocrine populations for single-nucleus RNA sequencing

BACKGROUND

Current approaches to profile the single-cell transcriptomics of human pancreatic endocrine cells almost exclusively rely on freshly isolated islets. However, human islets are limited in availability. Furthermore, the extensive processing steps during islet isolation and subsequent single cell dissolution might alter gene expressions. In this work, we report the development of a single-nucleus RNA sequencing (snRNA-seq) approach with targeted islet cell enrichment for endocrine-population focused transcriptomic profiling using frozen archival pancreatic tissues without islet isolation.

RESULTS

We cross-compared five nuclei isolation protocols and selected the citric acid method as the best strategy to isolate nuclei with high RNA integrity and low cytoplasmic contamination from frozen archival human pancreata. We innovated fluorescence-activated nuclei sorting based on the positive signal of NKX2-2 antibody to enrich nuclei of the endocrine population from the entire nuclei pool of the pancreas. Our sample preparation procedure generated high-quality single-nucleus gene-expression libraries while preserving the endocrine population diversity. In comparison with single-cell RNA sequencing (scRNA-seq) library generated with live cells from freshly isolated human islets, the snRNA-seq library displayed comparable endocrine cellular composition and cell type signature gene expression. However, between these two types of libraries, differential enrichments of transcripts belonging to different functional classes could be observed.

CONCLUSIONS

Our work fills a technological gap and helps to unleash frozen archival pancreatic tissues for molecular profiling targeting the endocrine population. This study opens doors to retrospective mappings of endocrine cell dynamics in pancreatic tissues of complex histopathology. We expect that our protocol is applicable to enrich nuclei for transcriptomics studies from various populations in different types of frozen archival tissues.

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.

LDHB contributes to the regulation of lactate levels and basal insulin secretion in human pancreatic β cells

Using 13C6 glucose labeling coupled to gas chromatography-mass spectrometry and 2D 1H-13C heteronuclear single quantum coherence NMR spectroscopy, we have obtained a comparative high-resolution map of glucose fate underpinning β cell function. In both mouse and human islets, the contribution of glucose to the tricarboxylic acid (TCA) cycle is similar. Pyruvate fueling of the TCA cycle is primarily mediated by the activity of pyruvate dehydrogenase, with lower flux through pyruvate carboxylase. While the conversion of pyruvate to lactate by lactate dehydrogenase (LDH) can be detected in islets of both species, lactate accumulation is 6-fold higher in human islets. Human islets express LDH, with low-moderate LDHA expression and β cell-specific LDHB expression. LDHB inhibition amplifies LDHA-dependent lactate generation in mouse and human β cells and increases basal insulin release. Lastly, cis-instrument Mendelian randomization shows that low LDHB expression levels correlate with elevated fasting insulin in humans. Thus, LDHB limits lactate generation in β cells to maintain appropriate insulin release.

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.

ApoA-I Protects Pancreatic β-Cells From Cholesterol-Induced Mitochondrial Damage and Restores Their Ability to Secrete Insulin

BACKGROUND

High cholesterol levels in pancreatic β-cells cause oxidative stress and decrease insulin secretion. β-cells can internalize apo (apolipoprotein) A-I, which increases insulin secretion. This study asks whether internalization of apoA-I improves β-cell insulin secretion by reducing oxidative stress.

METHODS

Ins-1E cells were cholesterol-loaded by incubation with cholesterol-methyl-β-cyclodextrin. Insulin secretion in the presence of 2.8 or 25 mmol/L glucose was quantified by radioimmunoassay. Internalization of fluorescently labeled apoA-I by β-cells was monitored by flow cytometry. The effects of apoA-I internalization on β-cell gene expression were evaluated by RNA sequencing. ApoA-I-binding partners on the β-cell surface were identified by mass spectrometry. Mitochondrial oxidative stress was quantified in β-cells and isolated islets with MitoSOX and confocal microscopy.

RESULTS

An F1-ATPase β-subunit on the β-cell surface was identified as the main apoA-I-binding partner. β-cell internalization of apoA-I was time-, concentration-, temperature-, cholesterol-, and F1-ATPase β-subunit-dependent. β-cells with internalized apoA-I (apoA-I+ cells) had higher cholesterol and cell surface F1-ATPase β-subunit levels than β-cells without internalized apoA-I (apoA-I- cells). The internalized apoA-I colocalized with mitochondria and was associated with reduced oxidative stress and increased insulin secretion. The IF1 (ATPase inhibitory factor 1) attenuated apoA-I internalization and increased oxidative stress in Ins-1E β-cells and isolated mouse islets. Differentially expressed genes in apoA-I+ and apoA-I- Ins-1E cells were related to protein synthesis, the unfolded protein response, insulin secretion, and mitochondrial function.

CONCLUSIONS

These results establish that β-cells are functionally heterogeneous, and apoA-I restores insulin secretion in β-cells with elevated cholesterol levels by improving mitochondrial redox balance.

Repression of latent NF-κB enhancers by PDX1 regulates β cell functional heterogeneity

Interactions between lineage-determining and activity-dependent transcription factors determine single-cell identity and function within multicellular tissues through incompletely known mechanisms. By assembling a single-cell atlas of chromatin state within human islets, we identified β cell subtypes governed by either high or low activity of the lineage-determining factor pancreatic duodenal homeobox-1 (PDX1). β cells with reduced PDX1 activity displayed increased chromatin accessibility at latent nuclear factor κB (NF-κB) enhancers. Pdx1 hypomorphic mice exhibited de-repression of NF-κB and impaired glucose tolerance at night. Three-dimensional analyses in tandem with chromatin immunoprecipitation (ChIP) sequencing revealed that PDX1 silences NF-κB at circadian and inflammatory enhancers through long-range chromatin contacts involving SIN3A. Conversely, Bmal1 ablation in β cells disrupted genome-wide PDX1 and NF-κB DNA binding. Finally, antagonizing the interleukin (IL)-1β receptor, an NF-κB target, improved insulin secretion in Pdx1 hypomorphic islets. Our studies reveal functional subtypes of single β cells defined by a gradient in PDX1 activity and identify NF-κB as a target for insulinotropic therapy.

Transcriptional activation of the Myc gene by glucose in β-cells requires a ChREBP-dependent 3-D chromatin interaction between the Myc and Pvt1 genes

OBJECTIVE

All forms of diabetes result from insufficient functional β-cell mass. Thus, achieving the therapeutic goal of expanding β-cell mass requires a better mechanistic understanding of how β-cells proliferate. Glucose is a natural β-cell mitogen that mediates its effects in part through the glucose-responsive transcription factor, carbohydrate response element binding protein (ChREBP) and the anabolic transcription factor, MYC. However, mechanistic details by which glucose activates Myc at the transcriptional level are poorly understood.

METHODS

Here, siRNA was used to test the role of ChREBP in the glucose response of MYC, ChIP and ChIPseq to identify potential regulatory binding sites, chromatin conformation capture to identify DNA/DNA interactions, and an adenovirus was constructed to expresses x-dCas9 and an sgRNA that specifically disrupts the recruitment of ChREBP to a specific targeted ChoRE.

RESULTS

We found that ChREBP is essential for glucose-mediated transcriptional induction of Myc, and for increases in Myc mRNA and protein abundance. Further, ChIPseq revealed that the carbohydrate response element (ChoRE) nearest to the Myc transcriptional start site (TSS) is immediately upstream of the gene encoding the lncRNA, Pvt1, 60,000 bp downstream of the Myc gene. Chromatin Conformation Capture (3C) confirmed a glucose-dependent interaction between these two sites. Transduction with an adenovirus expressing x-dCas9 and an sgRNA specifically targeting the highly conserved Pvt1 ChoRE, attenuates ChREBP recruitment, decreases Myc-Pvt1 DNA/DNA interaction, and decreases expression of the Pvt1 and Myc genes in response to glucose. Importantly, isolated and dispersed rat islet cells transduced with the ChoRE-disrupting adenovirus also display specific decreases in ChREBP-dependent, glucose-mediated expression of Pvt1 and Myc, as well as decreased glucose-stimulated β-cell proliferation.

CONCLUSIONS

The mitogenic glucose response of Myc is mediated via glucose-dependent recruitment of ChREBP to the promoter of the Pvt1 gene and subsequent DNA looping with the Myc promoter.

JOINTLY: interpretable joint clustering of single-cell transcriptomes

Single-cell and single-nucleus RNA-sequencing (sxRNA-seq) is increasingly being used to characterise the transcriptomic state of cell types at homeostasis, during development and in disease. However, this is a challenging task, as biological effects can be masked by technical variation. Here, we present JOINTLY, an algorithm enabling joint clustering of sxRNA-seq datasets across batches. JOINTLY performs on par or better than state-of-the-art batch integration methods in clustering tasks and outperforms other intrinsically interpretable methods. We demonstrate that JOINTLY is robust against over-correction while retaining subtle cell state differences between biological conditions and highlight how the interpretation of JOINTLY can be used to annotate cell types and identify active signalling programs across cell types and pseudo-time. Finally, we use JOINTLY to construct a reference atlas of white adipose tissue (WATLAS), an expandable and comprehensive community resource, in which we describe four adipocyte subpopulations and map compositional changes in obesity and between depots.

Beta cell lipotoxicity in the development of type 2 diabetes: the need for species-specific understanding

The propensity to develop type 2 diabetes (T2D) is known to have both environmental and hereditary components. In those with a genetic predisposition to T2D, it is widely believed that elevated concentrations of circulatory long-chain fatty acids (LC-FFA) significantly contribute towards the demise of insulin-producing pancreatic β-cells - the fundamental feature of the development of T2D. Over 25 years of research support that LC-FFA are deleterious to β-cells, through a process termed lipotoxicity. However, the work underpinning the theory of β-cell lipotoxicity is mostly based on rodent studies. Doubts have been raised as to whether lipotoxicity also occurs in humans. In this review, we examine the evidence, both in vivo and in vitro, for the pathogenic effects of LC-FFA on β-cell viability and function in humans, highlighting key species differences. In this way, we aim to uncover the role of lipotoxicity in the human pathogenesis of T2D and motivate the need for species-specific understanding.

Automatic quality control of single-cell and single-nucleus RNA-seq using valiDrops

Single-cell and single-nucleus RNA-sequencing (sxRNA-seq) measures gene expression in individual cells or nuclei enabling comprehensive characterization of cell types and states. However, isolation of cells or nuclei for sxRNA-seq releases contaminating RNA, which can distort biological signals, through, for example, cell damage and transcript leakage. Thus, identifying barcodes containing high-quality cells or nuclei is a critical analytical step in the processing of sxRNA-seq data. Here, we present valiDrops, an automated method to identify high-quality barcodes and flag dead cells. In valiDrops, barcodes are initially filtered using data-adaptive thresholding on community-standard quality metrics, and subsequently, valiDrops uses a novel clustering-based approach to identify barcodes with distinct biological signals. We benchmark valiDrops and show that biological signals from cell types and states are more distinct, easier to separate and more consistent after filtering by valiDrops compared to existing tools. Finally, we show that valiDrops can predict and flag dead cells with high accuracy. This novel classifier can further improve data quality or be used to identify dead cells to interrogate the biology of cell death. Thus, valiDrops is an effective and easy-to-use method to improve data quality and biological interpretation. Our method is openly available as an R package at www.github.com/madsen-lab/valiDrops.

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.

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

Interference with the HNF4-dependent gene regulatory network diminishes endoplasmic reticulum stress in hepatocytes

BACKGROUND

In all eukaryotic cell types, the unfolded protein response (UPR) upregulates factors that promote protein folding and misfolded protein clearance to help alleviate endoplasmic reticulum (ER) stress. Yet, ER stress in the liver is uniquely accompanied by the suppression of metabolic genes, the coordination and purpose of which are largely unknown.

METHODS

Here, we combined in silico machine learning, in vivo liver-specific deletion of the master regulator of hepatocyte differentiation HNF4α, and in vitro manipulation of hepatocyte differentiation state to determine how the UPR regulates hepatocyte identity and toward what end.

RESULTS

Machine learning identified a cluster of correlated genes that were profoundly suppressed by persistent ER stress in the liver. These genes, which encode diverse functions including metabolism, coagulation, drug detoxification, and bile synthesis, are likely targets of the master regulator of hepatocyte differentiation HNF4α. The response of these genes to ER stress was phenocopied by liver-specific deletion of HNF4α. Strikingly, while deletion of HNF4α exacerbated liver injury in response to an ER stress challenge, it also diminished UPR activation and partially preserved ER ultrastructure, suggesting attenuated ER stress. Conversely, pharmacological maintenance of hepatocyte identity in vitro enhanced sensitivity to stress.

CONCLUSIONS

Together, our findings suggest that the UPR regulates hepatocyte identity through HNF4α to protect ER homeostasis even at the expense of liver function.

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.

Reversing pancreatic β-cell dedifferentiation in the treatment of type 2 diabetes

The maintenance of glucose homeostasis is fundamental for survival and health. Diabetes develops when glucose homeostasis fails. Type 2 diabetes (T2D) is characterized by insulin resistance and pancreatic β-cell failure. The failure of β-cells to compensate for insulin resistance results in hyperglycemia, which in turn drives altered lipid metabolism and β-cell failure. Thus, insulin secretion by pancreatic β-cells is a primary component of glucose homeostasis. Impaired β-cell function and reduced β-cell mass are found in diabetes. Both features stem from a failure to maintain β-cell identity, which causes β-cells to dedifferentiate into nonfunctional endocrine progenitor-like cells or to trans-differentiate into other endocrine cell types. In this regard, one of the key issues in achieving disease modification is how to reestablish β-cell identity. In this review, we focus on the causes and implications of β-cell failure, as well as its potential reversibility as a T2D 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.

Transfer learning enables predictions in network biology

Mapping gene networks requires large amounts of transcriptomic data to learn the connections between genes, which impedes discoveries in settings with limited data, including rare diseases and diseases affecting clinically inaccessible tissues. Recently, transfer learning has revolutionized fields such as natural language understanding1,2 and computer vision3 by leveraging deep learning models pretrained on large-scale general datasets that can then be fine-tuned towards a vast array of downstream tasks with limited task-specific data. Here, we developed a context-aware, attention-based deep learning model, Geneformer, pretrained on a large-scale corpus of about 30 million single-cell transcriptomes to enable context-specific predictions in settings with limited data in network biology. During pretraining, Geneformer gained a fundamental understanding of network dynamics, encoding network hierarchy in the attention weights of the model in a completely self-supervised manner. Fine-tuning towards a diverse panel of downstream tasks relevant to chromatin and network dynamics using limited task-specific data demonstrated that Geneformer consistently boosted predictive accuracy. Applied to disease modelling with limited patient data, Geneformer identified candidate therapeutic targets for cardiomyopathy. Overall, Geneformer represents a pretrained deep learning model from which fine-tuning towards a broad range of downstream applications can be pursued to accelerate discovery of key network regulators and candidate therapeutic targets.

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.

Senescence: a double-edged sword in beta-cell health and failure?

Cellular senescence is a complex process marked by permanent cell-cycle arrest in response to a variety of stressors, and acts as a safeguard against the proliferation of damaged cells. Senescence is not only a key process underlying aging and development of many diseases, but has also been shown to play a vital role in embryogenesis as well as tissue regeneration and repair. In context of the pancreatic beta-cells, that are essential for maintaining glucose homeostasis, replicative senescence is responsible for the age-related decline in regenerative capacity. Stress induced premature senescence is also a key early event underlying beta-cell failure in both type 1 and type 2 diabetes. Targeting senescence has therefore emerged as a promising therapeutic avenue for diabetes. However, the molecular mechanisms that mediate the induction of beta-cell senescence in response to various stressors remain unclear. Nor do we know if senescence plays any role during beta-cell growth and development. In this perspective, we discuss the significance of senescence in beta-cell homeostasis and pathology and highlight emerging directions in this area that warrant our attention.

Understanding cell fate acquisition in stem-cell-derived pancreatic islets using single-cell multiome-inferred regulomes

Pancreatic islet cells derived from human pluripotent stem cells hold great promise for modeling and treating diabetes. Differences between stem-cell-derived and primary islets remain, but molecular insights to inform improvements are limited. Here, we acquire single-cell transcriptomes and accessible chromatin profiles during in vitro islet differentiation and pancreas from childhood and adult donors for comparison. We delineate major cell types, define their regulomes, and describe spatiotemporal gene regulatory relationships between transcription factors. CDX2 emerged as a regulator of enterochromaffin-like cells, which we show resemble a transient, previously unrecognized, serotonin-producing pre-β cell population in fetal pancreas, arguing against a proposed non-pancreatic origin. Furthermore, we observe insufficient activation of signal-dependent transcriptional programs during in vitro β cell maturation and identify sex hormones as drivers of β cell proliferation in childhood. Altogether, our analysis provides a comprehensive understanding of cell fate acquisition in stem-cell-derived islets and a framework for manipulating cell identities and maturity.

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.

Stomach-derived human insulin-secreting organoids restore glucose homeostasis

Gut stem cells are accessible by biopsy and propagate robustly in culture, offering an invaluable resource for autologous cell therapies. Insulin-producing cells can be induced in mouse gut, but it has not been possible to generate abundant and durable insulin-secreting cells from human gut tissues to evaluate their potential as a cell therapy for diabetes. Here we describe a protocol to differentiate cultured human gastric stem cells into pancreatic islet-like organoids containing gastric insulin-secreting (GINS) cells that resemble β-cells in molecular hallmarks and function. Sequential activation of the inducing factors NGN3 and PDX1-MAFA led human gastric stem cells onto a distinctive differentiation path, including a SOX4High endocrine and GalaninHigh GINS precursor, before adopting β-cell identity, at efficiencies close to 70%. GINS organoids acquired glucose-stimulated insulin secretion in 10 days and restored glucose homeostasis for over 100 days in diabetic mice after transplantation, providing proof of concept for a promising approach to treat diabetes.

Pancreatic islet cell type-specific transcriptomic changes during pregnancy and postpartum

Pancreatic β-cell mass expands during pregnancy and regresses in the postpartum period in conjunction with dynamic metabolic demands on maternal glucose homeostasis. To understand transcriptional changes driving these adaptations in β-cells and other islet cell types, we performed single-cell RNA sequencing on islets from virgin, late gestation, and early postpartum mice. We identified transcriptional signatures unique to gestation and the postpartum in β-cells, including induction of the AP-1 transcription factor subunits and other genes involved in the immediate-early response (IEGs). In addition, we found pregnancy and postpartum-induced changes differed within each endocrine cell type, and in endothelial cells and antigen-presenting cells within islets. Together, our data reveal insights into cell type-specific transcriptional changes responsible for adaptations by islet cells to pregnancy and their resolution postpartum.

Insulin secretion deficits in a Prader-Willi syndrome β-cell model are associated with a concerted downregulation of multiple endoplasmic reticulum chaperones

Prader-Willi syndrome (PWS) is a multisystem disorder with neurobehavioral, metabolic, and hormonal phenotypes, caused by loss of expression of a paternally-expressed imprinted gene cluster. Prior evidence from a PWS mouse model identified abnormal pancreatic islet development with retention of aged insulin and deficient insulin secretion. To determine the collective roles of PWS genes in β-cell biology, we used genome-editing to generate isogenic, clonal INS-1 insulinoma lines having 3.16 Mb deletions of the silent, maternal- (control) and active, paternal-allele (PWS). PWS β-cells demonstrated a significant cell autonomous reduction in basal and glucose-stimulated insulin secretion. Further, proteomic analyses revealed reduced levels of cellular and secreted hormones, including all insulin peptides and amylin, concomitant with reduction of at least ten endoplasmic reticulum (ER) chaperones, including GRP78 and GRP94. Critically, differentially expressed genes identified by whole transcriptome studies included reductions in levels of mRNAs encoding these secreted peptides and the group of ER chaperones. In contrast to the dosage compensation previously seen for ER chaperones in Grp78 or Grp94 gene knockouts or knockdown, compensation is precluded by the stress-independent deficiency of ER chaperones in PWS β-cells. Consistent with reduced ER chaperones levels, PWS INS-1 β-cells are more sensitive to ER stress, leading to earlier activation of all three arms of the unfolded protein response. Combined, the findings suggest that a chronic shortage of ER chaperones in PWS β-cells leads to a deficiency of protein folding and/or delay in ER transit of insulin and other cargo. In summary, our results illuminate the pathophysiological basis of pancreatic β-cell hormone deficits in PWS, with evolutionary implications for the multigenic PWS-domain, and indicate that PWS-imprinted genes coordinate concerted regulation of ER chaperone biosynthesis and β-cell secretory pathway function.

A beta cell subset with enhanced insulin secretion and glucose metabolism is reduced in type 2 diabetes

The pancreatic islets are composed of discrete hormone-producing cells that orchestrate systemic glucose homeostasis. Here we identify subsets of beta cells using a single-cell transcriptomic approach. One subset of beta cells marked by high CD63 expression is enriched for the expression of mitochondrial metabolism genes and exhibits higher mitochondrial respiration compared with CD63lo beta cells. Human and murine pseudo-islets derived from CD63hi beta cells demonstrate enhanced glucose-stimulated insulin secretion compared with pseudo-islets from CD63lo beta cells. We show that CD63hi beta cells are diminished in mouse models of and in humans with type 2 diabetes. Finally, transplantation of pseudo-islets generated from CD63hi but not CD63lo beta cells into diabetic mice restores glucose homeostasis. These findings suggest that loss of a specific subset of beta cells may lead to diabetes. Strategies to reconstitute or maintain CD63hi beta cells may represent a potential anti-diabetic therapy.

Single-cell RNA sequencing combined with single-cell proteomics identifies the metabolic adaptation of islet cell subpopulations to high-fat diet in mice

AIMS/HYPOTHESIS

Islets have complex heterogeneity and subpopulations. Cell surface markers representing alpha, beta and delta cell subpopulations are urgently needed for investigations to explore the compositional changes of each subpopulation in obesity progress and diabetes onset, and the adaptation mechanism of islet metabolism induced by a high-fat diet (HFD).

METHODS

Single-cell RNA sequencing (scRNA-seq) was applied to identify alpha, beta and delta cell subpopulation markers in an HFD-induced mouse model of glucose intolerance. Flow cytometry and immunostaining were used to sort and assess the proportion of each subpopulation. Single-cell proteomics was performed on sorted cells, and the functional status of each alpha, beta and delta cell subpopulation in glucose intolerance was deeply elucidated based on protein expression.

RESULTS

A total of 33,999 cells were analysed by scRNA-seq and clustered into eight populations, including alpha, beta and delta cells. For alpha cells, scRNA-seq revealed that the Ace2^low subpopulation had downregulated expression of genes related to alpha cell function and upregulated expression of genes associated with beta cell characteristics in comparison with the Ace2^high subpopulation. The impaired function and increased fragility of ACE2^low alpha cells exposure to HFD was further suggested by single-cell proteomics. As for beta cells, the CD81^high subpopulation may indicate an immature signature of beta cells compared with the CD81^low subpopulation, which had robust function. We also found differential expression of Slc2a2 in delta cells and a potentially stronger cellular function and metabolism in GLUT2^low delta cells than GLUT2^high delta cells. Moreover, an increased proportion of ACE2^low alpha cells and CD81^low beta cells, with a constant proportion of GLUT2^low delta cells, were observed in HFD-induced glucose intolerance.

CONCLUSIONS/INTERPRETATION

We identified ACE2, CD81 and GLUT2 as surface markers to distinguish, respectively, alpha, beta and delta cell subpopulations with heterogeneous maturation and function. The changes in the proportion and functional status of islet endocrine subpopulations reflect the metabolic adaptation of islets to high-fat stress, which weakened the function of alpha cells and enhanced the function of beta and delta cells to bring about glycaemic homeostasis. Our findings provide a fundamental resource for exploring the mechanisms maintaining each islet endocrine subpopulation's fate and function in health and disease.

DATA AVAILABILITY

The scRNA-seq analysis datasets from the current study are available in the Gene Expression Omnibus (GEO) repository under the accession number GSE203376.

Integrated transcriptomics contrasts fatty acid metabolism with hypoxia response in β-cell subpopulations associated with glycemic control

BACKGROUND

Understanding how heterogeneous β-cell function impacts diabetes is imperative for therapy development. Standard single-cell RNA sequencing analysis illuminates some factors driving heterogeneity, but new strategies are required to enhance information capture.

RESULTS

We integrate pancreatic islet single-cell and bulk RNA sequencing data to identify β-cell subpopulations based on gene expression and characterize genetic networks associated with β-cell function in obese SM/J mice. We identify β-cell subpopulations associated with basal insulin secretion, hypoxia response, cell polarity, and stress response. Network analysis associates fatty acid metabolism and basal insulin secretion with hyperglycemic-obesity, while expression of Pdyn and hypoxia response is associated with normoglycemic-obesity.

CONCLUSIONS

By integrating single-cell and bulk islet transcriptomes, our study explores β-cell heterogeneity and identifies novel subpopulations and genetic pathways associated with β-cell function in obesity.

Sex differences in islet stress responses support female β cell resilience

OBJECTIVE

Pancreatic β cells play a key role in maintaining glucose homeostasis; dysfunction of this critical cell type causes type 2 diabetes (T2D). Emerging evidence points to sex differences in β cells, but few studies have examined male-female differences in β cell stress responses and resilience across multiple contexts, including diabetes. Here, we address the need for high-quality information on sex differences in β cell and islet gene expression and function using both human and rodent samples.

METHODS

In humans, we compared β cell gene expression and insulin secretion in donors with T2D to non-diabetic donors in both males and females. In mice, we generated a well-powered islet RNAseq dataset from 20-week-old male and female siblings with similar insulin sensitivity. Our unbiased gene expression analysis pointed to a sex difference in the endoplasmic reticulum (ER) stress response. Based on this analysis, we hypothesized female islets would be more resilient to ER stress than male islets. To test this, we subjected islets isolated from age-matched male and female mice to thapsigargin treatment and monitored protein synthesis, cell death, and β cell insulin production and secretion. Transcriptomic and proteomic analyses were used to characterize sex differences in islet responses to ER stress.

RESULTS

Our single-cell analysis of human β cells revealed sex-specific changes to gene expression and function in T2D, correlating with more robust insulin secretion in human islets isolated from female donors with T2D compared to male donors with T2D. In mice, RNA sequencing revealed differential enrichment of unfolded protein response pathway-associated genes, where female islets showed higher expression of genes linked with protein synthesis, folding, and processing. This differential expression was physiologically significant, as islets isolated from female mice were more resilient to ER stress induction with thapsigargin. Specifically, female islets showed a greater ability to maintain glucose-stimulated insulin production and secretion during ER stress compared with males.

CONCLUSIONS

Our data demonstrate sex differences in β cell gene expression in both humans and mice, and that female β cells show a greater ability to maintain glucose-stimulated insulin secretion across multiple physiological and pathological contexts.

Cell Replacement Therapy for Type 1 Diabetes Patients: Potential Mechanisms Leading to Stem-Cell-Derived Pancreatic β-Cell Loss upon Transplant

Cell replacement therapy using stem-cell-derived insulin-producing β-like cells (sBCs) has been proposed as a practical cure for patients with type one diabetes (T1D). sBCs can correct diabetes in preclinical animal models, demonstrating the promise of this stem cell-based approach. However, in vivo studies have demonstrated that most sBCs, similarly to cadaveric human islets, are lost upon transplantation due to ischemia and other unknown mechanisms. Hence, there is a critical knowledge gap in the current field concerning the fate of sBCs upon engraftment. Here we review, discuss effects, and propose additional potential mechanisms that could contribute toward β-cell loss in vivo. We summarize and highlight some of the literature on phenotypic loss in β-cells under both steady, stressed, and diseased diabetic conditions. Specifically, we focus on β-cell death, dedifferentiation into progenitors, trans-differentiation into other hormone-expressing cells, and/or interconversion into less functional β-cell subtypes as potential mechanisms. While current cell replacement therapy efforts employing sBCs carry great promise as an abundant cell source, addressing the somewhat neglected aspect of β-cell loss in vivo will further accelerate sBC transplantation as a promising therapeutic modality that could significantly enhance the life quality of T1D patients.

Pancreatic Islet Cells Response to IFNγ Relies on Their Spatial Location within an Islet

Type 1 diabetes (T1D) is an auto-immune disease characterized by the progressive destruction of insulin-producing pancreatic beta cells. While beta cells are the target of the immune attack, the other islet endocrine cells, namely the alpha and delta cells, can also be affected by the inflammatory milieu. Here, using a flow cytometry-based strategy, we compared the impact of IFNγ, one of the main cytokines involved in T1D, on the three endocrine cell subsets isolated from C57BL/6 mouse islets. RNA-seq analyses revealed that alpha and delta cells exposed in vitro to IFNγ display a transcriptomic profile very similar to that of beta cells, with an increased expression of inflammation key genes such as MHC class I molecules, the CXCL10 chemokine and the programmed death-ligand 1 (PD-L1), three hallmarks of IFNγ signaling. Interestingly, at low IFNγ concentration, we observed two beta cell populations (responders and non-responders) based on PD-L1 protein expression. Our data indicate that this differential sensitivity relies on the location of the cells within the islet rather than on the existence of two different beta cells subsets. The same findings were corroborated by the in vivo analysis of pancreatic islets from the non-obese diabetic mouse model of T1D, showing more intense PD-L1 staining on endocrine cells close to immune infiltrate. Collectively, our work demonstrates that alpha and delta cells are as sensitive as beta cells to IFNγ, and suggests a gradual diffusion of the cytokine into an islet. These observations provide novel insights into the in situ inflammatory processes occurring in T1D progression.

Biology-inspired data-driven quality control for scientific discovery in single-cell transcriptomics

BACKGROUND

Quality control (QC) of cells, a critical first step in single-cell RNA sequencing data analysis, has largely relied on arbitrarily fixed data-agnostic thresholds applied to QC metrics such as gene complexity and fraction of reads mapping to mitochondrial genes. The few existing data-driven approaches perform QC at the level of samples or studies without accounting for biological variation.

RESULTS

We first demonstrate that QC metrics vary with both tissue and cell types across technologies, study conditions, and species. We then propose data-driven QC (ddqc), an unsupervised adaptive QC framework to perform flexible and data-driven QC at the level of cell types while retaining critical biological insights and improved power for downstream analysis. ddqc applies an adaptive threshold based on the median absolute deviation on four QC metrics (gene and UMI complexity, fraction of reads mapping to mitochondrial and ribosomal genes). ddqc retains over a third more cells when compared to conventional data-agnostic QC filters. Finally, we show that ddqc recovers biologically meaningful trends in gradation of gene complexity among cell types that can help answer questions of biological interest such as which cell types express the least and most number of transcripts overall, and ribosomal transcripts specifically.

CONCLUSIONS

ddqc retains cell types such as metabolically active parenchymal cells and specialized cells such as neutrophils which are often lost by conventional QC. Taken together, our work proposes a revised paradigm to quality filtering best practices-iterative QC, providing a data-driven QC framework compatible with observed biological diversity.

Dynamic Ins2 Gene Activity Defines β-Cell Maturity States

Transcriptional and functional cellular specialization has been described for insulin-secreting β-cells of the endocrine pancreas. However, it is not clear whether β-cell heterogeneity is stable or reflects dynamic cellular states. We investigated the temporal kinetics of endogenous insulin gene activity using live cell imaging, with complementary experiments using FACS and single-cell RNA sequencing, in β-cells from Ins2GFP knockin mice. In vivo staining and FACS analysis of islets from Ins2GFP mice confirmed that at a given moment, ∼25% of β-cells exhibited significantly higher activity at the evolutionarily conserved insulin gene, Ins2. Live cell imaging over days captured Ins2 gene activity dynamics in single β-cells. Autocorrelation analysis revealed a subset of oscillating cells, with mean oscillation periods of 17 h. Increased glucose concentrations stimulated more cells to oscillate and resulted in higher average Ins2 gene activity per cell. Single-cell RNA sequencing showed that Ins2(GFP)HIGH β-cells were enriched for markers of β-cell maturity. Ins2(GFP)HIGH β-cells were also significantly less viable at all glucose concentrations and in the context of endoplasmic reticulum stress. Collectively, our results demonstrate that the heterogeneity of insulin production, observed in mouse and human β-cells, can be accounted for by dynamic states of insulin gene activity.

Non-coding variants disrupting a tissue-specific regulatory element in HK1 cause congenital hyperinsulinism

Gene expression is tightly regulated, with many genes exhibiting cell-specific silencing when their protein product would disrupt normal cellular function1. This silencing is largely controlled by non-coding elements, and their disruption might cause human disease2. We performed gene-agnostic screening of the non-coding regions to discover new molecular causes of congenital hyperinsulinism. This identified 14 non-coding de novo variants affecting a 42-bp conserved region encompassed by a regulatory element in intron 2 of the hexokinase 1 gene (HK1). HK1 is widely expressed across all tissues except in the liver and pancreatic beta cells and is thus termed a 'disallowed gene' in these specific tissues. We demonstrated that the variants result in a loss of repression of HK1 in pancreatic beta cells, thereby causing insulin secretion and congenital hyperinsulinism. Using epigenomic data accessed from public repositories, we demonstrated that these variants reside within a regulatory region that we determine to be critical for cell-specific silencing. Importantly, this has revealed a disease mechanism for non-coding variants that cause inappropriate expression of a disallowed gene.