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

Single-cell sequencing to multi-omics: technologies and applications

Cells, as the fundamental units of life, contain multidimensional spatiotemporal information. Single-cell RNA sequencing (scRNA-seq) is revolutionizing biomedical science by analyzing cellular state and intercellular heterogeneity. Undoubtedly, single-cell transcriptomics has emerged as one of the most vibrant research fields today. With the optimization and innovation of single-cell sequencing technologies, the intricate multidimensional details concealed within cells are gradually unveiled. The combination of scRNA-seq and other multi-omics is at the forefront of the single-cell field. This involves simultaneously measuring various omics data within individual cells, expanding our understanding across a broader spectrum of dimensions. Single-cell multi-omics precisely captures the multidimensional aspects of single-cell transcriptomes, immune repertoire, spatial information, temporal information, epitopes, and other omics in diverse spatiotemporal contexts. In addition to depicting the cell atlas of normal or diseased tissues, it also provides a cornerstone for studying cell differentiation and development patterns, disease heterogeneity, drug resistance mechanisms, and treatment strategies. Herein, we review traditional single-cell sequencing technologies and outline the latest advancements in single-cell multi-omics. We summarize the current status and challenges of applying single-cell multi-omics technologies to biological research and clinical applications. Finally, we discuss the limitations and challenges of single-cell multi-omics and potential strategies to address them.

Exploring the design of clinical research studies on the efficacy mechanisms in type 2 diabetes mellitus

This review examines the complexities of Type 2 Diabetes Mellitus (T2DM), focusing on the critical role of integrating omics technologies with traditional experimental methods. It underscores the advancements in understanding the genetic diversity of T2DM and emphasizes the evolution towards personalized treatment modalities. The paper analyzes a variety of omics approaches, including genomics, methylation, transcriptomics, proteomics, metabolomics, and intestinal microbiomics, delineating their substantial contributions to deciphering the multifaceted mechanisms underlying T2DM. Furthermore, the review highlights the indispensable role of non-omics experimental techniques in comprehending and managing T2DM, advocating for their integration in the development of tailored medicine and precision treatment strategies. By identifying existing research gaps and suggesting future research trajectories, the review underscores the necessity for a comprehensive, multidisciplinary approach. This approach synergistically combines clinical insights with cutting-edge biotechnologies, aiming to refine the management and therapeutic interventions of T2DM, and ultimately enhancing patient outcomes. This synthesis of knowledge and methodologies paves the way for innovative advancements in T2DM research, fostering a deeper understanding and more effective treatment of this complex condition.

Exploring pancreatic beta-cell subgroups and their connectivity

Functional pancreatic islet beta cells are essential to ensure glucose homeostasis across species from zebrafish to humans. These cells show significant heterogeneity, and emerging studies have revealed that connectivity across a hierarchical network is required for normal insulin release. Here, we discuss current thinking and areas of debate around intra-islet connectivity, cellular hierarchies and potential "controlling" beta-cell populations. We focus on methodologies, including comparisons of different cell preparations as well as in vitro and in vivo approaches to imaging and controlling the activity of human and rodent islet preparations. We also discuss the analytical approaches that can be applied to live-cell data to identify and study critical subgroups of cells with a disproportionate role in control Ca2+ dynamics and thus insulin secretion (such as "first responders", "leaders" and "hubs", as defined by Ca2+ responses to glucose stimulation). Possible mechanisms by which this hierarchy is achieved, its physiological relevance and how its loss may contribute to islet failure in diabetes mellitus are also considered. A glossary of terms and links to computational resources are provided.

Single cell multiome profiling of pancreatic islets reveals physiological changes in cell type-specific regulation associated with diabetes risk

Physiological variability in pancreatic cell type gene regulation and the impact on diabetes risk is poorly understood. In this study we mapped gene regulation in pancreatic cell types using single cell multiomic (joint RNA-seq and ATAC-seq) profiling in 28 non-diabetic donors in combination with single cell data from 35 non-diabetic donors in the Human Pancreas Analysis Program. We identified widespread associations with age, sex, BMI, and HbA1c, where gene regulatory responses were highly cell type- and phenotype-specific. In beta cells, donor age associated with hypoxia, apoptosis, unfolded protein response, and external signal-dependent transcriptional regulators, while HbA1c associated with inflammatory responses and gender with chromatin organization. We identified 10.8K loci where genetic variants were QTLs for cis regulatory element (cRE) accessibility, including 20% with lineage- or cell type-specific effects which disrupted distinct transcription factor motifs. Type 2 diabetes and glycemic trait associated variants were enriched in both phenotype- and QTL-associated beta cell cREs, whereas type 1 diabetes showed limited enrichment. Variants at 226 diabetes and glycemic trait loci were QTLs in beta and other cell types, including 40 that were statistically colocalized, and annotating target genes of colocalized QTLs revealed genes with putatively novel roles in disease. Our findings reveal diverse responses of pancreatic cell types to phenotype and genotype in physiology, and identify pathways, networks, and genes through which physiology impacts diabetes risk.

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.

Gluco-regulation & type 2 diabetes: entrenched misconceptions updated to new governing principles for gold standard management

Our understanding of type 2 diabetes (T2D) has evolved dramatically. Advances have upended entrenched dogmas pertaining to the onset and progression of T2D, beliefs that have prevailed from the early era of diabetes research-and continue to populate our medical textbooks and continuing medical education materials. This review article highlights key insights that lend new governing principles for gold standard management of T2D. From the historical context upon which old beliefs arose to new findings, this article outlines evidence and perspectives on beta cell function, the underlying defects in glucoregulation, the remediable nature of T2D, and, the rationale supporting the shift to complication-centric prescribing. Practical approaches translate this rectified understanding of T2D into strategies that fill gaps in current management practices of prediabetes through late type 2 diabetes.

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.

Subcellular Feature-Based Classification of α and β Cells Using Soft X-ray Tomography

The dysfunction of α and β cells in pancreatic islets can lead to diabetes. Many questions remain on the subcellular organization of islet cells during the progression of disease. Existing three-dimensional cellular mapping approaches face challenges such as time-intensive sample sectioning and subjective cellular identification. To address these challenges, we have developed a subcellular feature-based classification approach, which allows us to identify α and β cells and quantify their subcellular structural characteristics using soft X-ray tomography (SXT). We observed significant differences in whole-cell morphological and organelle statistics between the two cell types. Additionally, we characterize subtle biophysical differences between individual insulin and glucagon vesicles by analyzing vesicle size and molecular density distributions, which were not previously possible using other methods. These sub-vesicular parameters enable us to predict cell types systematically using supervised machine learning. We also visualize distinct vesicle and cell subtypes using Uniform Manifold Approximation and Projection (UMAP) embeddings, which provides us with an innovative approach to explore structural heterogeneity in islet cells. This methodology presents an innovative approach for tracking biologically meaningful heterogeneity in cells that can be applied to any cellular system.

Deconvolution of polygenic risk score in single cells unravels cellular and molecular heterogeneity of complex human diseases

Polygenic risk scores (PRSs) are commonly used for predicting an individual's genetic risk of complex diseases. Yet, their implication for disease pathogenesis remains largely limited. Here, we introduce scPRS, a geometric deep learning model that constructs single-cell-resolved PRS leveraging reference single-cell chromatin accessibility profiling data to enhance biological discovery as well as disease prediction. Real-world applications across multiple complex diseases, including type 2 diabetes (T2D), hypertrophic cardiomyopathy (HCM), and Alzheimer's disease (AD), showcase the superior prediction power of scPRS compared to traditional PRS methods. Importantly, scPRS not only predicts disease risk but also uncovers disease-relevant cells, such as hormone-high alpha and beta cells for T2D, cardiomyocytes and pericytes for HCM, and astrocytes, microglia and oligodendrocyte progenitor cells for AD. Facilitated by a layered multi-omic analysis, scPRS further identifies cell-type-specific genetic underpinnings, linking disease-associated genetic variants to gene regulation within corresponding cell types. We substantiate the disease relevance of scPRS-prioritized HCM genes and demonstrate that the suppression of these genes in HCM cardiomyocytes is rescued by Mavacamten treatment. Additionally, we establish a novel microglia-specific regulatory relationship between the AD risk variant rs7922621 and its target genes ANXA11 and TSPAN14. We further illustrate the detrimental effects of suppressing these two genes on microglia phagocytosis. Our work provides a multi-tasking, interpretable framework for precise disease prediction and systematic investigation of the genetic, cellular, and molecular basis of complex diseases, laying the methodological foundation for single-cell genetics.

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.

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.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

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 eight genomic protective alleles for mitochondrial diabetes by Kinship-graph convolutional network

AIMS

Nearly 85% of maternally inherited diabetes and deafness (MIDD) are caused by the m.3243A>G mutation in the mitochondrial DNA. However, the clinical phenotypes of MIDD may also be influenced by the nuclear genome, this study aimed to investigate nuclear genome variants that influence clinical phenotypes associated with m.3243A>G mutation in MIDD based on whole-genome sequencing of the patients belonging to pedigrees.

MATERIALS AND METHODS

We analyzed a whole-genome sequencing (WGS) dataset from blood samples of 38 MIDD patients with the m.3243A > G mutation belonging to 10 pedigrees, by developing a Kinship-graph convolutional network approach, called Ki-GCN, integrated with the conventional genome-wide association study (GWAS) methods.

RESULTS

We identified eight protective alleles in the nuclear genome that have protective effects against the onset of MIDD, related deafness, and also type 2 diabetes. Based on these eight protective alleles, we constructed an effective logistic regression model to predict the early or late onset of MIDD patients.

CONCLUSIONS

There are protective alleles in the nuclear genome that are associated with the onset-age of MIDD patients and might also provide protective effects on the deafness derived from MIDD patients.

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