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

Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycemia, beta cell glucotoxicity, and ultimately type 2 diabetes (T2D). In this study, we sought to explore the effects of hyperglycemia on human pancreatic islet (HPI) gene expression by exposing HPIs from two donors to low (2.8mM) and high (15.0mM) glucose concentrations over 24 hours, assaying the transcriptome at seven time points using single-cell RNA sequencing (scRNA-seq). We modeled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Across all cell types, we identified 1,528 genes associated with time, 1,185 genes associated with glucose exposure, and 845 genes associated with interaction effects between time and glucose. We clustered differentially expressed genes across cell types and found 347 modules of genes with similar expression patterns across time and glucose conditions, including two beta cell modules enriched in genes associated with T2D. Finally, by integrating genomic features from this study and genetic summary statistics for T2D and related traits, we nominate 363 candidate effector genes that may underlie genetic associations for T2D and related traits.

[1]  Xiaochen Bo,et al.  clusterProfiler 4.0: A universal enrichment tool for interpreting omics data , 2021, Innovation.

[2]  C. Langefeld,et al.  A practical solution to pseudoreplication bias in single-cell studies , 2021, Nature Communications.

[3]  Tao Yang,et al.  δ-Cells: The Neighborhood Watch in the Islet Community , 2021, Biology.

[4]  Trevor Hastie,et al.  Genetics of 35 blood and urine biomarkers in the UK Biobank , 2020, Nature Genetics.

[5]  Joshua D. Campbell,et al.  Celda: a Bayesian model to perform co-clustering of genes into modules and cells into subpopulations using single-cell RNA-seq data , 2020, bioRxiv.

[6]  Jacob C. Ulirsch,et al.  Leveraging polygenic enrichments of gene features to predict genes underlying complex traits and diseases , 2020, Nature Genetics.

[7]  T. Pers,et al.  Genetic mapping of etiologic brain cell types for obesity , 2020, eLife.

[8]  A. Gloyn,et al.  Insights into pancreatic islet cell dysfunction from type 2 diabetes mellitus genetics , 2020, Nature Reviews Endocrinology.

[9]  D. Chicco,et al.  The advantages of the Matthews correlation coefficient (MCC) over F1 score and accuracy in binary classification evaluation , 2020, BMC Genomics.

[10]  Kerstin B. Meyer,et al.  BBKNN: fast batch alignment of single cell transcriptomes , 2019, Bioinform..

[11]  Joshua D. Campbell,et al.  Decontamination of ambient RNA in single-cell RNA-seq with DecontX , 2019, Genome Biology.

[12]  Allon M Klein,et al.  Scrublet: Computational Identification of Cell Doublets in Single-Cell Transcriptomic Data. , 2019, Cell systems.

[13]  T. Wolfsberg,et al.  Single-cell transcriptomics from human pancreatic islets: sample preparation matters. , 2019, Biology methods & protocols.

[14]  Vincent A. Traag,et al.  From Louvain to Leiden: guaranteeing well-connected communities , 2018, Scientific Reports.

[15]  G. Rutter,et al.  The α-cell in diabetes mellitus , 2018, Nature Reviews Endocrinology.

[16]  Anthony J. Payne,et al.  Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps , 2018, Nature Genetics.

[17]  Fabian J Theis,et al.  SCANPY: large-scale single-cell gene expression data analysis , 2018, Genome Biology.

[18]  C. Ling,et al.  The effects of high glucose exposure on global gene expression and DNA methylation in human pancreatic islets , 2017, Molecular and Cellular Endocrinology.

[19]  L. Groop,et al.  Glucose-Induced Changes in Gene Expression in Human Pancreatic Islets: Causes or Consequences of Chronic Hyperglycemia , 2017, Diabetes.

[20]  Amy K. Schmid,et al.  Clustering gene expression time series data using an infinite Gaussian process mixture model , 2017, bioRxiv.

[21]  K. Kaestner,et al.  PAX6 maintains &bgr; cell identity by repressing genes of alternative islet cell types , 2017, The Journal of clinical investigation.

[22]  J. K. Oeser,et al.  G6PC2 Modulates Fasting Blood Glucose In Male Mice in Response to Stress. , 2016, Endocrinology.

[23]  S. O’Rahilly,et al.  PCSK1 Mutations and Human Endocrinopathies: From Obesity to Gastrointestinal Disorders. , 2016, Endocrine reviews.

[24]  P. Linsley,et al.  MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data , 2015, Genome Biology.

[25]  L. Velloso,et al.  Delta cell death in the islet of Langerhans and the progression from normal glucose tolerance to type 2 diabetes in non-human primates (baboon, Papio hamadryas) , 2015, Diabetologia.

[26]  Emma Ahlqvist,et al.  Identification of novel genes for glucose metabolism based upon expression pattern in human islets and effect on insulin secretion and glycemia. , 2015, Human molecular genetics.

[27]  M. Fornage,et al.  Contribution of common non-synonymous variants in PCSK1 to body mass index variation and risk of obesity: a systematic review and meta-analysis with evidence from up to 331 175 individuals. , 2015, Human molecular genetics.

[28]  P. Franks,et al.  The association of common variants in PCSK1 with obesity: a HuGE review and meta-analysis. , 2014, American journal of epidemiology.

[29]  Emmanuelle Gouillart,et al.  scikit-image: image processing in Python , 2014, PeerJ.

[30]  G. Gyapay,et al.  A nonsense loss-of-function mutation in PCSK1 contributes to dominantly inherited human obesity , 2014, International Journal of Obesity.

[31]  Thomas Meitinger,et al.  Loss-of-function mutations in SLC30A8 protect against type 2 diabetes , 2014, Nature Genetics.

[32]  Chandler J. Faulman,et al.  G6PC2: A Negative Regulator of Basal Glucose-Stimulated Insulin Secretion , 2013, Diabetes.

[33]  C. Poisson,et al.  Pax6 is crucial for β-cell function, insulin biosynthesis, and glucose-induced insulin secretion. , 2012, Molecular endocrinology.

[34]  F. Ashcroft,et al.  Diabetes Mellitus and the β Cell: The Last Ten Years , 2012, Cell.

[35]  David E. Irwin,et al.  Finding a "Kneedle" in a Haystack: Detecting Knee Points in System Behavior , 2011, 2011 31st International Conference on Distributed Computing Systems Workshops.

[36]  L. Welschen,et al.  Combined effects of single-nucleotide polymorphisms in GCK, GCKR, G6PC2 and MTNR1B on fasting plasma glucose and type 2 diabetes risk , 2009, Diabetologia.

[37]  Douglas Steinley,et al.  Stability analysis in K-means clustering. , 2008, The British journal of mathematical and statistical psychology.

[38]  Jean Tichet,et al.  A Polymorphism Within the G6PC2 Gene Is Associated with Fasting Plasma Glucose Levels , 2008, Science.

[39]  Delbert Dueck,et al.  Clustering by Passing Messages Between Data Points , 2007, Science.

[40]  R. Scharfmann,et al.  Activating mutations in the ABCC8 gene in neonatal diabetes mellitus. , 2006, The New England journal of medicine.

[41]  Thomas Lengauer,et al.  A new measure for functional similarity of gene products based on Gene Ontology , 2006, BMC Bioinformatics.

[42]  F. Ashcroft,et al.  A heterozygous activating mutation in the sulphonylurea receptor SUR1 (ABCC8) causes neonatal diabetes. , 2006, Human molecular genetics.

[43]  P. Sham,et al.  A note on the calculation of empirical P values from Monte Carlo procedures. , 2002, American journal of human genetics.

[44]  R. Cattell.  The Scree Test For The Number Of Factors. , 1966, Multivariate behavioral research.

[45]  Stephen C. J. Parker,et al.  Integrative analysis of gene expression, DNA methylation, physiological traits, and genetic variation in human skeletal muscle. , 2019, Proceedings of the National Academy of Sciences of the United States of America.

[46]  Stephen Lynch,et al.  Image Processing with Python , 2018 .

[47]  Karin Kraft,et al.  [Type-2 diabetes]. , 2010, MMW Fortschritte der Medizin.

[48]  N. Kaiser,et al.  Glucotoxicity and beta-cell failure in type 2 diabetes mellitus. , 2003, Journal of pediatric endocrinology & metabolism : JPEM.

[49]  Masatsugu Hori,et al.  PAX6 mutation as a genetic factor common to aniridia and glucose intolerance. , 2002, Diabetes.

[50]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[51]  G. Abecasis,et al.  Supporting Online Material Materials and Methods Figs. S1 to S8 Tables S1 to S10 References a Genome-wide Association Study of Type 2 Diabetes in Finns Detects Multiple Susceptibility Variants , 2022 .