Multi-omics profiling of living human pancreatic islet donors reveals heterogeneous beta cell trajectories toward type 2 diabetes

Existing studies do not sufficiently describe the molecular changes of pancreatic islet beta cells leading to their deficient insulin secretion in type 2 diabetes (T2D). Here we address this deficiency with a comprehensive multi-omics analysis of metabolically profiled pancreatectomized living human donors stratified along the glycemic continuum from normoglycemia to T2D. Islet pools isolated from surgical samples by laser-capture microdissection had remarkably heterogeneous transcriptomic and proteomic profiles in diabetics, but not in non-diabetic controls. Transcriptomics analysis of this unique cohort revealed islet genes already dysregulated in prediabetic individuals with impaired glucose tolerance. Our findings demonstrate a progressive but disharmonic remodeling of mature beta cells, challenging current hypotheses of linear trajectories toward precursor or trans-differentiation stages in T2D. Further, integration of islet transcriptomics and pre-operative blood plasma lipidomics data enabled us to define the relative importance of gene co-expression modules and lipids positively or negatively associated with HbA1c levels, pointing to potential prognostic markers.

[1]  Jürgen Cox,et al.  MaxQuant Software for Ion Mobility Enhanced Shotgun Proteomics* , 2020, Molecular & Cellular Proteomics.

[2]  Monther Alhamdoosh,et al.  RNA-seq analysis is easy as 1-2-3 with limma, Glimma and edgeR , 2016, F1000Research.

[3]  Age K. Smilde,et al.  Real-life metabolomics data analysis : how to deal with complex data ? , 2010 .

[4]  A. Murphy,et al.  RNA Sequencing of Single Human Islet Cells Reveals Type 2 Diabetes Genes. , 2016, Cell metabolism.

[5]  Christian Gieger,et al.  Novel biomarkers for pre-diabetes identified by metabolomics , 2012, Molecular systems biology.

[6]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[7]  D. Zozulinska-Ziolkiewicz,et al.  209-LB: Is the Abdominal Aorta Diameter a Novel Marker of Diabetes Risk in Elderly Women? , 2019, Diabetes.

[8]  Roman Fischer,et al.  MaxQuant Software for Ion Mobility Enhanced Shotgun Proteomics* , 2019, Molecular & Cellular Proteomics.

[9]  Y. Kuroda,et al.  Activation of Macrophage-Associated Molecules after Brain Death in Islets , 2003, Cell transplantation.

[10]  S. Wingett,et al.  FastQ Screen: A tool for multi-genome mapping and quality control. , 2018, F1000Research.

[11]  Steve Horvath,et al.  WGCNA: an R package for weighted correlation network analysis , 2008, BMC Bioinformatics.

[12]  Tommi Suvitaival,et al.  Lipidome as a predictive tool in progression to type 2 diabetes in Finnish men. , 2018, Metabolism: clinical and experimental.

[13]  K. Kaestner,et al.  Virgin Beta Cells Persist throughout Life at a Neogenic Niche within Pancreatic Islets. , 2017, Cell metabolism.

[14]  R. Jewkes,et al.  Perceptions and Experiences of Research Participants on Gender-Based Violence Community Based Survey: Implications for Ethical Guidelines , 2012, PloS one.

[15]  B. Thorens.  GLUT2, glucose sensing and glucose homeostasis , 2015, Diabetologia.

[16]  Colin N. Dewey,et al.  RNA-Seq gene expression estimation with read mapping uncertainty , 2009, Bioinform..

[17]  M. A. Surma,et al.  An automated shotgun lipidomics platform for high throughput, comprehensive, and quantitative analysis of blood plasma intact lipids , 2015, European journal of lipid science and technology : EJLST.

[18]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[19]  M. Pędziwiatr,et al.  Type 2 Diabetes Remission 5 Years After Laparoscopic Sleeve Gastrectomy: Multicenter Cohort Study , 2020, Obesity Surgery.

[20]  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.

[21]  P. Froguel,et al.  Improved protocol for laser microdissection of human pancreatic islets from surgical specimens. , 2013, Journal of visualized experiments : JoVE.

[22]  D. Rutledge,et al.  Iterative weighting of multiblock data in the orthogonal partial least squares framework. , 2014, Analytica chimica acta.

[23]  Julien Boccard,et al.  A consensus orthogonal partial least squares discriminant analysis (OPLS-DA) strategy for multiblock Omics data fusion. , 2013, Analytica chimica acta.

[24]  Jürgen Cox,et al.  1D and 2D annotation enrichment: a statistical method integrating quantitative proteomics with complementary high-throughput data , 2012, BMC Bioinformatics.

[25]  Anton J. Enright,et al.  Kraken: A set of tools for quality control and analysis of high-throughput sequence data , 2013, Methods.

[26]  R. Sladek,et al.  Analysis of Beta-Cell Gene Expression Reveals Inflammatory Signaling and Evidence of Dedifferentiation following Human Islet Isolation and Culture , 2012, PloS one.

[27]  A. Shevchenko,et al.  Detection of Independent Associations of Plasma Lipidomic Parameters with Insulin Sensitivity Indices Using Data Mining Methodology , 2016, PloS one.

[28]  Marco Y. Hein,et al.  The Perseus computational platform for comprehensive analysis of (prote)omics data , 2016, Nature Methods.

[29]  J. George,et al.  Single-cell transcriptomes identify human islet cell signatures and reveal cell-type–specific expression changes in type 2 diabetes , 2017, Genome research.

[30]  Simon Andrews,et al.  FastQ Screen: A tool for multi-genome mapping and quality control , 2018, F1000Research.

[31]  F. Ashcroft,et al.  Diabetes causes marked inhibition of mitochondrial metabolism in pancreatic β-cells , 2019, Nature Communications.

[32]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[33]  Benjamin S. Aribisala,et al.  Reversal of type 2 diabetes: normalisation of beta cell function in association with decreased pancreas and liver triacylglycerol , 2011, Diabetologia.

[34]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[35]  Ioannis Xenarios,et al.  Plasma Dihydroceramides Are Diabetes Susceptibility Biomarker Candidates in Mice and Humans. , 2017, Cell reports.

[36]  Marco Y. Hein,et al.  A “Proteomic Ruler” for Protein Copy Number and Concentration Estimation without Spike-in Standards* , 2014, Molecular & Cellular Proteomics.

[37]  Wei Li,et al.  RSeQC: quality control of RNA-seq experiments , 2012, Bioinform..

[38]  C. Nichols,et al.  Pancreatic β cell dedifferentiation in diabetes and redifferentiation following insulin therapy. , 2014, Cell metabolism.

[39]  J. Shaw,et al.  Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: results from the International Diabetes Federation Diabetes Atlas, 9th edition. , 2019, Diabetes research and clinical practice.

[40]  D. Sgroi,et al.  Evidence of stress in β cells obtained with laser capture microdissection from pancreases of brain dead donors , 2017, Islets.

[41]  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.

[42]  M. Solimena,et al.  The Expression of Aldolase B in Islets Is Negatively Associated With Insulin Secretion in Humans , 2018, The Journal of clinical endocrinology and metabolism.

[43]  I. Xenarios,et al.  Laser capture microdissection of human pancreatic islets reveals novel eQTLs associated with type 2 diabetes , 2019, Molecular metabolism.

[44]  M. Mann,et al.  MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification , 2008, Nature Biotechnology.

[45]  Paul Theodor Pyl,et al.  HTSeq—a Python framework to work with high-throughput sequencing data , 2014, bioRxiv.

[46]  R. DeFronzo,et al.  Plasma Ceramides Are Elevated in Obese Subjects With Type 2 Diabetes and Correlate With the Severity of Insulin Resistance , 2009, Diabetes.

[47]  Johan Trygg,et al.  K-OPLS package: Kernel-based orthogonal projections to latent structures for prediction and interpretation in feature space , 2008, BMC Bioinformatics.

[48]  Jun Hu,et al.  OSA: a fast and accurate alignment tool for RNA-Seq , 2012, Bioinform..

[49]  D. Aust,et al.  Metabolically phenotyped pancreatectomized patients as living donors for the study of islets in health and diabetes , 2019, Molecular metabolism.

[50]  P. Marchetti,et al.  A call for improved reporting of human islet characteristics in research articles , 2018, Diabetologia.

[51]  Philippe Froguel,et al.  Decreased STARD10 Expression Is Associated with Defective Insulin Secretion in Humans and Mice , 2017, American journal of human genetics.

[52]  Marco Y. Hein,et al.  Accurate Proteome-wide Label-free Quantification by Delayed Normalization and Maximal Peptide Ratio Extraction, Termed MaxLFQ * , 2014, Molecular & Cellular Proteomics.

[53]  Guangchuang Yu,et al.  clusterProfiler: an R package for comparing biological themes among gene clusters. , 2012, Omics : a journal of integrative biology.

[54]  D. Accili,et al.  Evidence of β-Cell Dedifferentiation in Human Type 2 Diabetes. , 2016, The Journal of clinical endocrinology and metabolism.

[55]  G. Weir.  Glucolipotoxicity, β-Cells, and Diabetes: The Emperor Has No Clothes , 2020, Diabetes.

[56]  C. Talchai,et al.  Pancreatic β Cell Dedifferentiation as a Mechanism of Diabetic β Cell Failure , 2012, Cell.

[57]  J. Weitz,et al.  Dysfunction of Persisting β Cells Is a Key Feature of Early Type 2 Diabetes Pathogenesis. , 2020, Cell reports.

[58]  M. Mann,et al.  Minimal, encapsulated proteomic-sample processing applied to copy-number estimation in eukaryotic cells , 2014, Nature Methods.

[59]  Melvin A. Park,et al.  Online Parallel Accumulation–Serial Fragmentation (PASEF) with a Novel Trapped Ion Mobility Mass Spectrometer* , 2018, Molecular & Cellular Proteomics.

[60]  F. Ashcroft,et al.  New insights into KATP channel gene mutations and neonatal diabetes mellitus , 2020, Nature Reviews Endocrinology.

[61]  E. Bonifacio,et al.  Systems biology of the IMIDIA biobank from organ donors and pancreatectomised patients defines a novel transcriptomic signature of islets from individuals with type 2 diabetes , 2017, Diabetologia.