Proteomic predictors of individualized nutrient-specific insulin secretion in health and disease
================================================================================================
* Jelena Kolic
* WenQing Grace Sun
* Haoning Howard Cen
* Jessica Ewald
* Jason C. Rogalski
* Shugo Sasaki
* Han Sun
* Varsha Rajesh
* Yi Han Xia
* Renata Moravcova
* Søs Skovsø
* Aliya F. Spigelman
* Jocelyn E. Manning Fox
* James Lyon
* Leanne Beet
* Jianguo Xia
* Francis C. Lynn
* Anna L. Gloyn
* Leonard J. Foster
* Patrick E. MacDonald
* James D. Johnson
## ABSTRACT
Population level variation and molecular mechanisms behind insulin secretion in response to carbohydrate, protein, and fat remain uncharacterized despite ramifications for personalized nutrition. Here, we define prototypical insulin secretion dynamics in response to the three macronutrients in islets from 140 cadaveric donors, including those diagnosed with type 2 diabetes. While islets from the majority of donors exhibited the expected relative response magnitudes, with glucose being highest, amino acid moderate, and fatty acid small, 9% of islets stimulated with amino acid and 8% of islets stimulated with fatty acids had larger responses compared with high glucose. We leveraged this insulin response heterogeneity and used transcriptomics and proteomics to identify molecular correlates of specific nutrient responsiveness, as well as those proteins and mRNAs altered in type 2 diabetes. We also examine nutrient-responsiveness in stem cell-derived islet clusters and observe that they have dysregulated fuel sensitivity, which is a hallmark of functionally immature cells. Our study now represents the first comparison of dynamic responses to nutrients and multi-omics analysis in human insulin secreting cells. Responses of different people’s islets to carbohydrate, protein, and fat lay the groundwork for personalized nutrition.
**ONE-SENTENCE SUMMARY** Deep phenotyping and multi-omics reveal individualized nutrient-specific insulin secretion propensity.
## INTRODUCTION
Insulin is released by pancreatic islet beta cells in response to nutrient stimuli to maintain energy homeostasis. The major driver of insulin secretion is glucose. However, proteins and fats may also modulate insulin release and the effects of non-carbohydrate nutrients on insulin secretion remain underexplored. Much of our limited understanding of nutrient-stimulated insulin secretion is extrapolated from rodents, although more recent availability of cadaveric human islets for research purposes has expanded our pre-clinical knowledge. However, current human islet datasets generally only examine a single nutrient stimulus, glucose1–3. A few small studies examined other nutrients4,5, but no large-scale direct comparison of insulin secretion in human islets, stimulated by carbohydrates, proteins and fats has been reported. Understanding of nutrient-induced insulin secretion is important in the context type 2 diabetes and emerging studies linking hyperinsulinemia with cancer6. Indeed, large prospective clinical trials shows broad beneficial effects of diets targeting hyperinsulinemia7.
Individuals respond differently to diets8 and there is high interpersonal variability in postprandial responses to even one macronutrient, glucose9. *Ex vivo studies* measuring insulin secretion from human islets show high variability, only some of which can be explained by donor characteristics or islet isolation parameters2. The concept that insulin-responses to food types or different macronutrients are individualized has not been investigated. No studies have leveraged macronutrient-induced insulin secretion heterogeneity and large-scale multi-omics to elucidate the associated molecular mechanisms.
Here we address these critical knowledge gaps by measuring dynamic insulin secretion in response to three model macronutrient stimuli in islets from non-diabetic and donors with type 2 diabetes, as well as stem cell-derived islets. Our transcriptomic and proteomic analysis revealed molecular signatures of each islet donor and identified distinct clusters of proteins that predict insulin secretion responses to carbohydrates, proteins, and fats. We also report on a subset of donor islets that islets that are relatively hyperresponsive lipids in a way that resembles functionally immature human embryonic stem cell-derived islets. This is the largest human islet dataset that includes both macronutrient-stimulated insulin secretion measurements and multi-omic profiling, coupled with the first side-by-side comparison of nutrient responses and proteomes between human islets and stem cell-derived islet clusters. This resource helps us understand why how individuals’ islets respond differently to sugar, protein, and fat and advocates for greater application of personalized recommendations and treatments for individuals living with diabetes.
## RESULTS
### Prototypical nutrient-stimulated insulin secretion dynamics from human islets
Between 2016 and 2022, we systematically measured insulin secretion in response to carbohydrate (15 mM or 6 mM glucose), amino acid (5 mM leucine), and fat (1.5 mM oleate/palmitate mix) from islets isolated from 140 cadaveric donors2 reflective of the general population (Fig. 1A)10. In most islet donors, we confirmed that carbohydrate was the strongest insulin secretagogue, followed by amino acid and then fat, which only weakly stimulated insulin secretion on average (Figs. 1B, S1A,B). Islets exhibited biphasic insulin secretion in response to glucose11. We found that insulin secretion in response to amino acid is also biphasic, with a distinct 1st phase lasting ∼15 minutes and a sustained 2nd phase lasting the duration of the challenge (Figs. 1B, S1A,B). In contrast, the response to fat, when present, was monophasic (Figs. 1B, S1A,B).
[Fig. 1.](http://medrxiv.org/content/early/2024/03/05/2023.05.24.23290298/F1)
Fig. 1. Nutrient-stimulated insulin secretion, transcriptomes, and proteomes from metabolically healthy and T2D islets.
(**A**) Histogram with trend line of ages for both male (green, solid line) and female (yellow, dashed line) donors. Additional details of islet isolation parameters and donor characteristics are available in Table S1 and will be available on [www.humanislets.com](http://www.humanislets.com) **(B)** Averaged traces of dynamic insulin secretion measurements in response to glucose (15 or 6 mM), leucine (5 mM), 1.5 mM oleate and palmitate (1:1 mixture) or 30 mM KCl in islets isolated from non-diabetic donors (n=123) (*left)* and donors with type 2 diabetes (n=17) *(right)* are shown as indicated. Basal glucose was 3 mM. **(C)** Islet samples from the cadaveric donors were collected for RNAseq. Enriched mRNAs in islets from ND donors (n= 82) are shown in teal and enriched mRNAs in islets from donors with type 2 diabetes (n=8) are shown in red. The top 40 most significant differentially abundant proteins are highlighted by labelling with gene name. **(D)** Shows the same as E, except donor islets were collected for mass-spec based proteomics (n= 118) in ND, and (n=16) in T2D **(E)** *(top)* Venn diagrams show the overlap of differentially expressed mRNAs and abundant proteins, and this is further depicted in the correlation plot *(bottom)* **(F)** Across-gene correlation between RNAs and proteins **(G)** Differentially abundant proteins (with a greater than 0.5 log2 fold change) between the non-diabetic and type 2 diabetic donor islets were connected using STRING and depicted in the context of beta cell signalling. The color of the nodes represents the fold change (ND/T2D), while the thickness of the line around the nodes represents the p-value (adjusted).
We also assessed the role of macronutrient order on insulin secretion, inspired by clinical meal-order studies12,13. In our experiments, prior exposure to high glucose, amino acid or fat did not alter insulin secretion stimulated by moderately elevated 6 mM glucose (Figs. S1A,B). Amino acid on top of 6 mM glucose further increased insulin secretion (Figs. S1A). Fat produced no enhanced response (Figs. S1A,B), in contrast to previous small-scale rodent or human islet observations14,15. Interestingly, insulin secretion stimulated by direct depolarization with KCl was inhibited after prior exposure to lipid in islets from normoglycemic donors (Fig. S1A), perhaps foreshadowing lipotoxic effects on the insulin secretory machinery16. This large dataset provides the first side-by-side response profiles for each of the main three macronutrients in human islets.
### Prototypical nutrient-stimulated insulin secretion dynamics in type 2 diabetes
We next examined the relationship between type 2 diabetes and nutrient-stimulated insulin secretion. Islets from donors with family-reported type 2 diabetes (Table S1) had ∼40% lower insulin secretion in response to 15 mM glucose and ∼35% lower insulin secretion in response to moderate 6 mM glucose (Figs. 1B, S2C,F). Insulin secretion in response to direct depolarization by KCl was reduced by ∼22% but failed to reach statistical significance (Fig. S2I, p=0.07). When examining the kinetics of insulin secretion, we saw a significant delay in time-to-insulin peak in response to high glucose in islets from donors with diabetes (Fig. S1E). Insulin secretion in response to fatty acids was also lowered by ∼55% (Figs. 1B, S2E). However, insulin content, baseline insulin secretion, and insulin secretion in response to 5 mM leucine were not different (Figs. 1B, S2A,B,D). This preserved amino acid-stimulated insulin secretion is consistent with clinical data17 and suggests therapeutic protein intake could be exploited in diabetes management. However, leucine together with 6 mM glucose induced ∼36% (on average) less insulin secretion in donors with diabetes (Fig. S2G), and leucine was not stimulatory on top of 6 mM glucose in these donors (Fig. S1B), emphasizing a need for additional clinical research into context-dependent amino acid-stimulated insulin secretion. Known glucose-lowering medication status of these donors with diabetes did not significantly correlate with any insulin-secretion parameters; however, there was an overall trend with the need for exogenous insulin and lower overall insulin secretory capacity (Fig. S2J-Q).
### Comprehensive transcriptomics and proteomics of islets with and without type 2 diabetes
To better understand the relationship between secretory response and variation in gene expression in islets from donors with and without type 2 diabetes, we performed comprehensive transcriptomic and proteomic analysis from corresponding batches of donor islets. On average, after quality control, we measured over 20,000 mRNAs using RNA sequencing (from 82 ND and 8 T2D donors) and almost 8000 proteins (from 118 ND and 16 T2D donors) using mass-spectrometry. We first used our multi-omic dataset to estimate islet cell composition in our donor preps based on abundance of proteins typically enriched in islet and acinar cells18. The model indicates that islets from donors with T2D have altered islet cell composition, with lower beta cell and higher alpha and delta cell markers, as well as a higher percentage of acinar cell markers (Fig. S3B,D).
Using RNA sequencing of a large subset of human islets, we next showed that 286 mRNAs were enriched in the non-diabetic donors and 146 were enriched in the donors with diabetes (Fig. 1C, Table S2). Importantly, differences were more apparent at the protein level, with 355 proteins significantly more abundant in islets from non-diabetic donors, and 200 proteins more abundant in islets from donors with diabetes (Fig. 1D, Table S3). Despite only finding 48 gene products that were significantly altered at both mRNA and protein levels (Fig. 1 E), individual protein fold changes were mostly consistent with mRNA fold changes (Fig. S4A), and both lists of features were enriched with many of the same pathways (Fig. S4B). We also used an alternative approach to compare the protein and RNA changes in T2D by utilizing the rank-rank hypergeometric overlap test19, which identified ∼3300 proteins and RNAs that were changed in the same direction in T2D (Fig. S5). However, we found that this method was not stringent enough for our purposes since the vast majority of the overlapping proteins and RNAs it identified had negligible fold change difference in T2D (Fig. S5).
Gene products decreased in both expression and abundance in islets from donors with diabetes included the sulfonylurea receptor subunit of the KATP channel (*ABCC8*), the well-known target of oral hypoglycaemic agents, as well as pyruvate carboxylase (*PC*), which is critical for mitochondrial metabolism and glucose-stimulated insulin release20, hedgehog acyltransferase-like (*HHATL*), which negatively regulates protein palmitoylation, a process implicated in type 2 diabetes21 and islet amyloid polypeptide (*IAPP*), a hormone co-secreted with insulin with roles in glycemic control and gastric emptying22 (Fig. 1E). Only four gene products were increased in both expression and abundance in type 2 diabetes: RNA-binding protein with multiple splicing (RBPMS), AP-3 complex subunit sigma-1 (AP3S1) – involved in lysosomal trafficking23, syntaxin binding protein 6 (*STXBP6*) which is thought to limit insulin release by limiting the size of the granule fusion pore24 and beta arrestin 1(*ARRB1*) which is suggested to be involved in GLP-1 stimulated insulin secretion25.
There was overall discordance between islet RNA expression and protein abundance (7609 gene products) from the same donor/isolations, with an R value of 0.5 (Fig. 1F), consistent with other studies of primary tissues26,27. As such, we did not detect significant RNA-protein correlations for 81% of the gene products, which interestingly mostly encoded for mitochondria-related and ribosomal proteins (Fig. S4C,D,E). This discrepancy was likely not due to unreliable RNAseq measurements because we observed a strong correlation between the RNAseq and NanoString analysis of 130 mRNAs (Fig. S4F,G). Because proteins had lower co-efficients of variation (Fig. S4H), were more stable in the face of isolation variables (see below), and may provide better insight into disease phenotype than mRNAs27, we focused on the proteomic data to elucidate key mechanisms. We mapped the protein-protein interaction networks for proteins showing fold-change greater than 0.5 log2 using STRING28, and subsequently assigned them to their intracellular location in a diagram using subcellular location information found in UniProt (Fig. 1G). Consistent with a lower insulin secretory response to stimulatory glucose (Fig. S2C,F), islets from donors with type 2 diabetes had lower abundance of key proteins predicted to be involved in glucose-stimulated insulin secretion. Islet proteins reduced in type 2 diabetes (Fig. 1 G) included mitochondrial proteins arginase 2 (ARG2), N-acetylaspartate synthetase (NAT8L), hydroxyacyl-coenzyme A dehydrogenase, mitochondrial (HADH), protein phosphatase 1K, mitochondrial (PPM1K) and alpha-methylacyl-CoA racemase (AMACR), regulators of glycolysis (6-phosphofructo-2-kinase/fructose-2, 6-bisphosphatase 2 (PFKFB2) and glucose-6-phosphatase 2 (G6PC2), endoplasmic reticulum Ca2+ homeostasis factors ERO1-like protein beta (ERO1B), inactive ubiquitin thioesterase (OTULINL), sarcoplasmic/endoplasmic reticulum calcium ATPase 3 (ATP2A3), and the extracellular Ca2+-sensing receptor (CASR). Proteins that were more abundant in islets from donors with diabetes included those crucial for cell adhesion (tight junction protein ZO-2 (TJP2), integrin alpha-2/beta-1 (ITGA2), Golgi function arylsulfatase L (ARSL), transmembrane BAX inhibitor motif containing 1 (TMBIM1), transcriptional proteins such as TLE family member 3 (TLE3), helicase with zinc finger domain 2 (HELZ2) and RNA-binding protein with multiple splicing (RBPMS), and proteins involved in cytoskeletal reorganization like MARCKS-related protein (MARCKSK1). Collectively, these data identify multiple proteins that correlate with the insufficient glucose-stimulated insulin secretion in islets from donors with diabetes.
### Individuality of nutrient-stimulated insulin secretion
Islets from individual donors exhibited a large range of insulin secretion rates at baseline, and in response to high glucose, moderate glucose, amino acid, fat, and direct depolarization (Fig. 2A). Surprisingly, we observed that some donors had more robust responses to fatty acids than to glucose, challenging the long-standing idea that dietary fats alone have negligible effects on insulin release5,14,29. This degree of response heterogeneity across all macronutrients was not found in C57Bl6J mouse islets, even when including both sexes and a wide range of ages (Fig. 2B).
[Fig. 2.](http://medrxiv.org/content/early/2024/03/05/2023.05.24.23290298/F2)
Fig. 2. Clustering heterogenous insulin responses to macronutrients and donor proteomes.
(**A**) *(left)* Individual traces of dynamic insulin secretion stimulated by glucose (15 mM or 6 mm) or KCl (30 mM). Basal glucose was 3 mM. Average responses from non-diabetic donors are illustrated with solid teal line and average responses from donors with type 2 diabetes are in shown in the dashed salmon line. Floating dot plot inserts illustrate the heterogeneity in insulin AUC for the corresponding section of the perifusion curve (salmon dots illustrate the responses from donors with T2D). *(middle)* illustrates the same as *(left)* except islets were stimulated with 5 mM leucine alone, or in combination with 6 mM glucose as depicted in figure panel. *(right)* illustrates the same as *(left)*, except islets were stimulated with 1.5 mM oleate/palmitate (1:1 mix), or in combination with 6 mM glucose. **(B)** Illustrates the same as (A) except in mouse islets of both sexes (7 males, 9 females, 7-90 weeks of age). Floating box plots illustrate the variance of AUC responses between the mouse islet (grey) and human islet (brown) responses. **(C)** Co-correlation analysis identified 18 distinct modules of proteins whose abundance shows similar patterns. Modules are annotated with KEGG pathways and GO terms. The Pearson correlation coefficients between the modules and the donor meta data or functional data (x-axis) are shown in each rectangle, and the corresponding adjusted p value is shown in the brackets. Positive correlations are shown in shades of red and negative in shades of blue. Significant correlations are highlighted with black boxes. **(D)** Illustrates the main connection in the co-expression network. The WGCNA adjacency matrix was filtered to remove: 1) any protein not in a module, 2) protein-protein adjacency distances less than 0.2, and 3) any protein without any connections after filters 1 and 2. The resulting network contained 1600 protein nodes and 16 209 edges. Each protein node is colored according to its module (same color scheme as C).
We explored the source of this variation by examining known donor characteristics (Table S1) and islet isolation parameters. Cold ischemic time of the pancreas was negatively correlated with insulin secretion in response to direct depolarization (Fig. S6P), but it was not correlated to insulin secretion stimulated by any of the three macronutrients tested. We also found that islets from female donors had lower insulin secretion at 3 mM glucose and 6 mM glucose, with or without fatty acid treatment (Fig. S6A-C). Donor age and BMI were not different between male and female donors in our study (Table S1), suggesting that these differences are a result of biological sex. Somewhat surprisingly, donor BMI had only minor effects on macronutrient-stimulated insulin secretion (Fig. S6D-L). Expectedly, high HbA1c was negatively correlated to insulin secretion stimulated by high glucose and KCl (Fig. S6M-O). These observations suggest that donor characteristics and islet isolation parameters minimally contribute to individualized nutrient responses and that separate factors must drive the observed heterogeneity.
### Co-expression network analysis uncovers protein-networks correlated with islet function
We used co-expression network analysis to obtain an overview of protein-protein relationships within the proteome (and transcriptome), and then analyzed the network along with clinical and functional outcomes to gain insight into the molecular drivers of macronutrient-stimulated insulin secretion response heterogeneity. We identified 18 network modules of highly co-regulated proteins (Figs. 2C, S7) and plotted them as a co-expression network to illustrate the connections between modules (Fig. 2D). All modules were significantly enriched in gene sets (Gene Ontology Cellular Component, Molecular Function, and Biological Process) and pathways (KEGG and Reactome), allowing for their higher-level functional annotation. The overall protein abundances of 9 modules were significantly correlated with islet functional data (Fig. 2C). Several modules were positively correlated with insulin secretion in response to multiple nutrients. These modules (lightcyan, dark green, greenyellow and black) contained proteins with critical roles in cytoskeletal reorganization, mitochondrial metabolism, and insulin processing. Two modules of interest (blue and royalblue) with roles in protein transport and localization, were positively correlated with insulin secretion stimulated by KCl but negatively correlated with insulin secretion stimulated by lipids. Notable proteins in this module include enzymes involved in fatty acid synthesis, acetyl-CoA carboxylase 1 (ACACA) and fatty acid synthase (FASN)30,31, and lipid metabolism glycerol-3-phosphate dehydrogenase **(**GPD1)32. Consistent with the discordance between RNA expression and protein abundance, and suggestive of more dynamic mRNA levels, we did not observe any significantly correlated transcript modules in the 29 network modules detected (Fig. S8). Together, these data demonstrate that these protein network modules in human islets likely underpin their differential responsiveness to stimuli.
### Associations between individual proteins, donor data and islet function
We first correlated the abundance of each mRNA (Fig. S9, Table S4) and protein (Fig. S10, Table S5) identified in our dataset to donor traits, independent of diabetes status. Some of these correlations were driven by disease status rather than islet function (Table S6), so we next reanalyzed our dataset using a linear regression model to adjust for T2D diagnosis (Fig. 3, S11). We focused our attention on the proteome because we found that mRNAs were poorly correlated with islet function (Fig. S9,11) and were more prominently associated with islet culture and isolation variables (Fig. S12, Tables S6-8).
[Fig. 3.](http://medrxiv.org/content/early/2024/03/05/2023.05.24.23290298/F3)
Fig. 3. Regression coefficients (adjusting for T2D status) of individual proteins with donor characteristics and insulin secretory responses.
Volcano plots are shown depicting significant positive (red) and negative (blue) linear regression coefficients of proteins to continuous donor characteristics: **(A)** donor age and **(B)** HbA1c; or functional parameters in response to: **(C)** 3 mM glucose **(D)** 15 mM glucose **(E)** 6 mM glucose **(F)** 5 mM leucine **(G)** 1.5 mM oleate/palmitate (1:1 molar ratio) **(H)** 1.5 mM oleate/palmitate + 6 mM glucose and **(I)** 30 mM KCl. Protein abundances (log2) are depicted by size of circle, and the coefficient of variation (log10) is depicted by color gradient. **(J)** Venn diagram showing the overlap of the number of protein abundances that positively associate with the indicated nutrient stimuli following adjustment for T2D disease status diagnosis **(K)** Heat map depicts the top 50 most positively associated proteins to the indicated stimuli. **(L-M)** Shows the same as (J-K), but for negative associations.
Our analysis suggests that increasing donor age was positively associated to proteins with known (and postulated) roles in lysosomal degradation, cellular senescence, inflammation and mitochondrial dysfunction: palmitoyl-protein thioesterase 1 (PPT1)33 scavenger receptor class B member 2 (SCARB2)34, acid ceramidase (ASAH1)35, purinergic receptor P2X436, and NADPH:adrenodoxin oxidoreductase (FDXR)37 (Fig. 3A, Table S8). We found that the abundances of 161 proteins were associated with HbA1c (Fig. 3B). Interestingly, the most significant positive relationships were seen to proteins with known roles in cancer progression and development: protein YIF1A38, ubiquitin associated protein 1 (UBAP1)39, DDRGK domain-containing protein 140 and syntaxin-5 (STX5)41. This is relevant to reports of poor glycemic control being associated with and increased risk of multiple cancers42.
We also looked for relationships between protein abundance and islet function in response to each nutrient stimulus (Fig. 3C-I). The greatest number of significant associations to protein abundance were seen with insulin secretion induced by lipids (2053 significant associations) (Fig. 3G, Tables S6,8), pointing to somewhat unique mechanisms for this macronutrient as a possible modulator of insulin release. The abundance of 3 proteins was positively associated to insulin secretion stimulated by all three macronutrients as well as direct depolarization by KCl (Fig. 3J,K), suggesting that these proteins play a role in the overall secretory capacity of the beta cell. These pan-stimulus enabling proteins were ADP-Ribosylation Factor GTPase 6 Interacting Protein 6 (ARL6IP6), solute carrier family 25 member 4 (SLC25A4), and secretory carrier membrane protein 4 (SCAMP4). These findings highlight the critical role of mitochondria in nutrient– and depolarization-stimulated insulin secretion43.
Interestingly, alpha cell transcription factor, aristaless related homeobox (ARX)44 was negatively associated with basal insulin secretion (Fig. 3C), insulin secretion stimulated by glucose (Fig. 3E), as well as insulin secretion stimulated by amino acid (Fig. 3F) (Table S8). While these correlations may initially suggest a general decrease in beta cell ability to secrete insulin on account of higher alpha cell mass, we also observe that insulin secretion induced by direct depolarization with KCl was not correlated to ARX abundance (Fig. 3M, Table S8), arguing against such a simple interpretation. Collectively, these data provide unprecedented information on the proteins that associate with islet secretory function, in a general and nutrient-specific way. We have built an online resource where the relationship between any detected protein, donor characteristics, islet isolation parameters and islet function can be displayed ([www.humanislets.com](http://www.humanislets.com)).
### Unique proteome of lipid– or amino acid-hyper-responders
During the course of our multi-nutrient phenotyping studies, we identified a previously un-reported sub-group of donors with unusually large insulin secretory response to lipids in basal glucose conditions. In fact, ∼8% of donors secreted more insulin in response to oleate/palmitate than to 15 mM glucose (Fig. 4A). These donors also secreted more insulin when lipid was present together with 6 mM glucose, but less insulin in response to direct depolarization by KCl (Fig. 4A, Table S9). Overall, lipid hyper-responsive islets came from donors with higher HbA1C, but no significant differences in known islet isolation or culture parameters were seen between the two groups (Table S9).
[Fig. 4.](http://medrxiv.org/content/early/2024/03/05/2023.05.24.23290298/F4)
Fig. 4. Prototypical proteomes and responses of fat and protein hyper-responders.
(**A**) The average insulin secretion response from donors classified as “high lipid responders” is shown in dark blue (n=11) and the average insulin secretion response from donors classified as “low lipid responders” is shown in the light blue (n=129). **(B)** Principal component analysis (PCA) plot of protein abundance in the high (dark blue circles) vs. low (light blue circles) lipid responders. **(C)** Volcano plot showing differential protein abundance between high and low lipid responders. **(F)** The average insulin secretion response from donors classified as “high protein responders” is shown in orange (n=13) and the average insulin secretion response from donors classified as “low protein responders” is shown in yellow (n=127). **(G)** PCA plot of protein abundance in the high (orange circles) vs. low (yellow circles) protein responders. **(H)** Volcano plot showing differential protein abundance between high and low protein responders.
However, proteomic comparison of these two donor groups showed that the fat hyper-responsive islets exhibited a 30% decrease in the abundance of Wnt family member 4 (WNT4), hinting at a less-mature state45 (Fig. 4B,C). Fat hyper-responsive islets also had less protein phosphatase 1 catalytic subunit gamma (PPP1CC), which is involved in glycogen metabolism, hypoxanthine phosphoribosyltransferase 1 (HPRT1), which is involved in purine metabolism, and AKT serine/threonine kinase 1 (AKT1), which is involved in growth factor signaling and survival. On the other hand, lipid hyper-responders exhibited elevated hepatocyte nuclear factor 1-beta (HNF1B), the gene product underlying a rare form of diabetes, maturity onset diabetes of the young (MODY5), sulfotransferase family 1A member 1 (SULT1A1), and platelet derived growth factor receptor alpha (PDGFRA) (Fig. 4B,C). Pathway analysis using multiple databases 46 reveals defects in small-molecule protein modification, AKT signaling and lipid metabolism in donors classified as fat-hyper responders (Fig. S13A). Interestingly, proteins involved in endoplasmic reticulum signal integration and organization are increased in fat-hyper responders (Fig. S13B), which is perhaps suggestive of endoplasmic reticulum dysfunction in these donors47.
We also identified a previously un-reported subset of donors (∼9%) that secreted more insulin in response to leucine than to high glucose (Fig. 4D, Table S11). Amino acid hyper-responders were more likely to been diagnosed with type 2 diabetes, but there was no difference in HbA1C between the high and low responders (Table S10). Amino acid hyper-response was also associated with a longer *in vitro* culture time (Table S10), suggesting that this observed phenotype may be influenced by adaptation to our culture conditions. Only two differentially abundant proteins (Fig. 4E,F) were associated with this phenotype, thymocyte selection-associated high mobility group box protein (TOX; transcriptional regulator) and angiotensinogen (AGT), an essential component in the renin-angiotensin system, with important roles in the endocrine pancreas48. Our work now suggests that some donors’ islet proteomes may be pre-programmed or may have adapted to *in vitro* culture conditions, which enables them to hyper-respond to lipids and/or amino acids.
### Human embryonic stem cell-derived islet-like clusters are lipid-responsive
The differentiation of beta cell surrogates from embryonic stem cells has potential as a diabetes therapy and can be used to model ‘environmentally naïve’ insulin secretion49. The *in vitro* insulin-secretory function of these cells, even in response to supraphysiological levels of glucose, still does not match that of human islets49. There are no reports that directly compare how stem cell-derived islets and primary human islets respond to multiple nutrients49, however recent work has suggested that human mitochondria from pluripotent stem cell-derived islets have a higher sensitivity for non-glucose fuels50. In our hands, INS2AGFP stem cell-derived islet-clusters51 (Fig. 5A,B) exhibited lower levels of typical beta cell enriched markers (Fig. S3C,D) and had ∼1/10 the level of basal insulin release when compared to human islets (Fig. 5C). They showed poor responsiveness to glucose, were unresponsive to amino acid (either alone or in combination with glucose) but did release significant amounts of insulin upon direct depolarization by KCl (Figs. 5C, S14). This suggests that these *in vitro* derived cells have the capacity for regulated insulin secretion but have specific defects in glucose and amino acid responses and may prefer non-glucose substrates50.
[Fig. 5.](http://medrxiv.org/content/early/2024/03/05/2023.05.24.23290298/F5)
Fig. 5. Stem cell-derived islet-like clusters hyper-respond to fat and can be matured with extended culture.
(**A**) Summary of the human embryonic stem cell-derived islet-like clusters differentiation protocol. **(B)** Representative images of unsorted (day 35) Sc-derived islet-like clusters *(left)* and FAC-sorted (day 35) Sc-derived beta cell-like clusters *(right)* **(C)** Averaged traces of dynamic insulin secretion measurements in response to glucose (15 or 6 mM), leucine (5 mM), 1.5 mM oleate and palmitate (1:1 mixture) or 30 mM KCl in stem cell derived islet-like clusters (average of 8 younger-immature and 8 older-maturing preparations). Basal glucose was 3 mM. For comparison, the average human islet dynamic insulin secretion measurement traces are shown in the dashed lines. **(D)** PCA plot of protein abundance in the stem cell-derived islet-like clusters (n=25) (blue circles) compared to donors with type 2 diabetes (salmon circles, n=16) and normoglycemic donors (teal circles) (n=118). **(E)** Volcano plot showing differential protein abundance between stem cell derived islet like clusters (blue circles, n=25) and donors without diabetes (teal circles, n=118). The top 40 most significant differentially abundant proteins are highlighted by labelling with gene name. **(F)** *(left)* Compares average dynamic insulin secretion stimulated by glucose (15 mM or 6 mm) or KCl (30 mM) between younger-immature (dashed brown lines, n=8) and older-maturing (solid brown lines, n=8) clusters. Basal glucose was 3 mM. Floating dot plot inserts illustrate the AUC responses for the corresponding section of the perifusion curve between younger-immature clusters (black circles) and older-maturing clusters (open circles). *(middle)* illustrates the same as *(left)* except islet-like clusters were stimulated with 5 mM leucine alone, or in combination with 6 mM glucose as depicted in figure panel. *(right)* illustrates the same as *(left)*, except islet-like were stimulated with 1.5 mM oleate/palmitate, or in combination with 6 mM glucose. **(G)** PCA plot and **(H)** volcano plot of protein abundance comparing the younger-immature clusters (light green, n=14) and older-maturing clusters (dark green, n=10).
Comparing their proteomes to human islets revealed massive differences; 3196 proteins were significantly more abundant in the non-diabetic human donor islets (Fig. 5D,E, Table S3). Some differences were due to the presence of non-islet exocrine cells trapped in our hand-picked preparations (Figs. 5E, S3). Other proteomic differences were more informative. WNT4, which we found was less abundant in lipid hyper-responders, was also more than 2 times lower in these hESC-derived islet-like clusters than in islets from non-diabetic donors (Table S3). On the other hand, 2553 proteins were more abundant in the stem cell-derived islet clusters (Table S3). Higher levels of this protein in the pancreatic beta cell are thought to limit voltage-dependent calcium channel activity, and thus decrease glucose-stimulated insulin secretion52. The abundance of monocarboxylate transporter 1 (SLC16A1), a beta-cell disallowed gene whose presence is suggestive a neonatal-like immature beta-cell phenotype53, increased over 6-fold in these islet-like clusters. Compellingly, proteins linked to important roles in fatty acid transport and metabolism long-chain fatty acid transport protein 3 (SLC27A3), mitochondrial coenzyme A transporter (SLC25A42), and solute carrier family 22 member 18 (SLC22A18) were increased ∼2-fold.
Extended maturation of stem cell-derived islet-like clusters has been shown to increase the ability of these cells to respond appropriately to a glucose stimulus54,55. We therefore compared nutrient-stimulated insulin secretion from younger-immature (days 21-35) and older-maturing (days 50+) stem cell preparations (Figs. 5F, S14). Nanostring analysis showed increased expression of key maturity markers (MAFA, IAPP and UCN3) in cells older than day 35 (Fig. S15). Basal insulin secretion was 5-fold higher in the older-maturing islet-like cells, and they also showed evidence of glucose-responsiveness, and maintained their sensitivity to fatty acid. At the level of the proteome, 2072 proteins more abundant in the older-maturing clusters and 2066 more abundant in the younger-immature clusters (Fig. 5G,H, Table S3). Proteins with significant roles in glucose-stimulated insulin secretion and insulin synthesis including pyruvate carboxylase (PC), mitochondrial citrate carrier (SLC25A1), high-voltage gated calcium channel subunit (CACNA2D2), glucose transporter 1 (SLC2A1), an integral SNARE protein (VAMP8) and prohormone convertase 2 (PCSK2) were all significantly more abundant in the older-maturing SC-islet clusters. These cells also had a 73% decrease in the abundance of the beta-cell disallowed gene monocarboxylate transporter (SLC16A1), when compared to more immature clusters. These proteomic data now provide a roadmap to produce more faithful human islet surrogates. Although the maturation of stem cell derived islets remains incomplete *in vitro*55,56, and current cell replacement tend to focus on beta-like cells, our dataset can now be used to make better stem cell-derived islets by exploiting the differences seen at the protein levels between these cells and real human islets.
## DISCUSSION
The goal of this study was to compare dynamic insulin secretion responses to carbohydrate, protein, and fat from human islets representative of the population of donors with and without type 2 diabetes, and to define the transcriptomic and proteomic mechanisms that underly the variation in response to each macronutrient. Our study is part of a large multi-institutional human islet deep phenotyping network which generates and links functional and multi-omic datasets for insight into normal and pathological variation in islet function. Our initiative provides unprecedented and unbiased mechanistic detail of the cellular processes associated with insulin secretion in response to each nutrient class, and the first to directly compare functional data to both the proteome and transcriptome. As a part of the Islet Deep Phenotyping Network Resource, this dataset forms the basis on an online platform ([www.humanislets.com](http://www.humanislets.com)) providing open access to islet donor/isolation data and additional phenotyping. Despite working with high-quality research islets, all isolated from the same centre2, isolation parameters can have unavoidable impacts; this was particularly obvious when examining mRNA, which was more sensitive than protein to technical isolation and culture parameters. Along with weak mRNA correlation to protein in our study and others’26,27, this should give pause when inferring protein abundance (or function) from mRNA expression.
Our findings have relevance for personalized therapeutic nutrition and the clinical management of diabetes. For example, while both first and second phase insulin secretion in response to high glucose were blunted in islets from donors with type 2 diabetes, we did not detect significantly impaired amino acid-stimulated insulin secretion, which bolsters the case that protein-rich diets could have therapeutic benefits in patients with type 2 diabetes57 and highlights the need for additional research into amino-acid stimulated insulin secretion.
A recent study sought to identify differences between islets from 17 normal donors and 12 donors with type 2 diabetes, but was only able to quantify 3036 proteins, none of which were reported statistically significant after correcting for multiple comparisons58. Here, with more donors and more than twice the proteomic coverage, we report 555 proteins that are differentially abundant in islets from donors with type 2 diabetes. This allowed us to conduct protein-centric network analysis that was not previously possible. Proteins critical to glucose and mitochondrial metabolism coupling are less abundant in islets from donors with type 2 diabetes, explaining the blunted glucose-stimulated insulin secretory response. Conversely, proteins involved in cytoskeletal actin remodeling are on average more abundant in islets from donors with type 2 diabetes, consistent with the role of filamentous actin limiting glucose-stimulated insulin secretion in pancreatic beta cells59. Identifying the pathways that are perturbed in islets from donors with type 2 diabetes is crucial for better understanding the disease pathology and discovering of novel therapeutic targets.
Previous studies revealed heterogeneity in insulin secretion from human islets in response to glucose, but they do not identify the source of this variation or examine other nutrients1–3. Here, we leveraged the heterogeneity in nutrient-stimulated insulin secretion to identify networks of proteins that drive nutrient-specific responses. Interestingly, we identified more significant correlations between protein and insulin secretion in response to physiologically relevant 6 mM glucose, when compared to 15 mM, calling into question the routine use of such supraphysiological glucose concentrations *in vitro*.
Heterogeneity in insulin secretion in response to the three main classes of macronutrients is reminiscent of the interindividual variability of postprandial glucose responses to carbohydrates9,60,61. It has been suggested that this interindividual variability in glycemic response is in part responsible for the mixed weight loss outcomes following specific diet interventions. A recent randomized clinical trial did not support this idea62, but post-prandial glucose and insulin responses were not measured to confirm the effectiveness of the reported personal diets. On the other hand, a diet specifically targeting hyperinsulinemia was shown to be highly effect in the prevention of multiple diseases7. Multi-omic profiling of blood plasma recently revealed heterogeneity in individual responses to a “healthy lifestyle intervention”63. Combined with our results, these studies provide the rationale for a clinical trial to test the insulin response to standardized macronutrients challenges in the general population of normoglycemic individuals and those living with type 2 diabetes.
One of our most surprising findings was that ∼8% of donor islets showed relative hyperresponsiveness to a fat stimulus at basal glucose. Strikingly, the greatest number of correlations to protein abundance were in response to the fatty acid stimulus, meaning that these islets are quite different from those exhibiting the typical response ratio. We speculate that this lipid responsiveness is related to beta cell immaturity because it is also observed in stem-cell derived immature cells. The magnitude of fatty acid-stimulated insulin secretion has been debated4,14,64, and we are not aware of clinical literature that has reported on individuals with higher lipid-stimulated insulin secretion compared with glucose-stimulated insulin secretion *in vivo*. Studies with sufficient power and participant diversity have not been conducted and reported yet. Our literature search of clinical studies examining the effects of glucose and oral or iv lipid infusions in human participants was only able to find studies with low n numbers that did not provide detailed information on the response distribution (Table S11)65–70, so direct comparisons with our dataset are not possible. We hope our results will encourage clinical researchers to look for response heterogeneity for each macronutrient in heterogenous cohorts >100.
In summary, we present evidence that nutrient-specific insulin secretion is heterogeneous in the general population. In the face of hyperinsulinemia being causal to numerous health conditions, our benchmark study lays the groundwork for future clinical studies aiming to advance the area of personalized therapeutic nutrition.
## DATA AND MATERIALS AVAILABILITY
All data are available, both via a purpose-build web portal and the appropriate public databases. The raw RNA sequencing data has been deposited in the European Genome-phenome Archive (EGA) under study number EGAS00001007241. The raw proteomics data has been deposited to ProteomExchange via MASSive under dataset number PXD045422. This dataset, together with the functional data will also be available for download and browsing on a new portal website [www.humanislets.com](http://www.humanislets.com).
## MATERIALS AND METHODS
### Human islet culture
Islets were isolated from pancreas of cadaveric human donors at the Alberta Diabetes Institute IsletCore and used with approval of the Human Research Ethics Board at the University of Alberta (Pro00013094; Pro00001754), the University of British Columbia (UBC) Clinical Research Ethics Board (H13-01865), the University of Oxford’s Oxford Tropical Research Ethics Committee (OxTREC Ref.2-15), the Oxfordshire Regional Ethics Committee B (REC reference: 09/H0605/2), or by the Stanford Center for Biomedical Ethics (IRB Protocol: 57310). T2D diagnosis was family-declared at the time of organ donation. All families of organ donors provided written informed consent for use in research. Islet isolation was performed according to the detailed methods deposited in the protocols.io repositor 71. Islets were shipped in CMRL media (Thermo Fisher Scientific) overnight from Edmonton to UBC. Upon arrival, islets were immediately purified by handpicking under a stereomicroscope and suspended in RPMI 1640 medium (Thermo Fisher Scientific, Cat#, 11879-020) supplemented with 5.5 mmol glucose, 10% fetal bovine serum (FBS), and 100 units/mL penicillin/streptomycin. Islets were cultured in 10 cm non-adhesive petri dishes (Thermo Fisher Scientific, Cat # FB0875713) at a density of 10-16 islets/cm2 (100-160 IEQ/ cm2) for 24-72 hours prior to experiments to allow for recovery from shipment. Total islet culture time (as well as additional donor and isolation characteristics) are listed in Table 1.
### Deep dynamic phenotyping of human and mouse islets, and stem cell-derived islet-like clusters
Our standard approach 72 compared the response to 6 mM or 15 mM glucose stimulation and direct depolarization with 30 mM KCl. In parallel, we also measured insulin secretion in response to 5 mM leucine (Sigma, Cat. # L8912, dissolved in 1M HCl, then pH adjusted with 1M NaOH) – an essential branched chain amino acid previously shown to increase insulin secretion in human islets5; or a 1:1 mixture of 0.75 mM oleic acid (Sigma, Cat. # 364525) and 0.75 mM palmitic acid (Sigma, Cat. #P5585) at either 3 mM or 6 mM glucose (see figure legends). Fatty acids stock solutions were prepared in 50% ethanol at 65°C for 30 minutes, and then added to a solution of fatty acid free bovine serum albumin (BSA) (Sigma, Cat. # A7030) in a 6:1 molar ratio in a 37 °C water bath for 60 min.
We loaded perifusion columns with either 65 human islets, 100 (hESC)-stem cell-derived islet-like clusters, or 100 mouse islets and perifused them at 0.4 mL/min with 3 mM glucose Krebs-Ringer Modified Buffer (KRB) solution as described previously 72 for 60 minutes to equilibrate the islets to the KRB and flow rate, and then with the indicated conditions. First-phase insulin release was defined as the amount of insulin secreted during the first 15 minutes of 15 mM glucose stimulation, while the remaining 25 min of stimulation were defined as second-phase release. Peak insulin secretion was defined as single point whereby the amount of insulin released during the first 15 minutes of a solution change was the highest. Samples were stored at –20 °C and insulin secretion was quantified using human (Millipore Cat. # HI-14K) or rat (Millipore Cat. # RI-13K) insulin radioimmunoassay kits. All insulin secretion units have been converted to pmol/islet, to allow for direct comparison between species.
### Mouse islet isolation and culture
C57BL/6J mice were purchased from the Jackson Laboratory and housed in the UBC Modified Barrier Facility (temperature-controlled) using protocols approved by the UBC Animal Care Committee and in accordance with international guidelines. The mice were on a 12-hour light/dark cycle with chow diet (LabDiet #5053) and drinking water *ad libitum*. Mouse islets were isolated by ductal collagenase (Type XI; Sigma, Cat. # C7657) injection followed by hand-picking. To allow for recovery post digestion, islets were cultured overnight in islet culture media (RPMI media with 11.1 mM D-glucose supplemented with 10% vol/vol fetal bovine serum (Thermo Fisher Scientific, Cat. #12483020) and 1% vol/vol Penicillin-Streptomycin (Gibco, Cat #15140-148)) at 37°C with 5% CO2.
### Human embryonic stem cell-derived islet-like cluster culture and differentiation
A human embryonic stem cell line, based on the WiCell WA01 line, that contained knock-in add-on of EGFP downstream of the insulin coding sequence 51 was differentiated, as previously described, through a 6-stage protocol to day 21, media components are included in the table below51,73. After the initial 21 days of differentiation, SC-β cells were cultured in CMRL with 5.6 mM glucose (Thermo Fisher Scientific) containing 1% fatty acid-free BSA (Proliant), 1:100 Glutamax, 1:100 NEAA, 1 mM pyruvate, 10 mM HEPES, 1:100 ITS (Thermo Fisher Scientific), 10 μg/ml of heparin sulfate, 1 mM N-acetyl cysteine, 10 μM zinc sulfate, 1.75 μL 2-mercaptoethanol, 2 μM T3 and 0.155 mM ascorbic acid (Sigma Aldrich) until Day 27. Between Day 28-50, aggregates were grown in CMRL with 5.6 mM glucose containing 1% fatty acid-free BSA, 1:100 Glutamax, 1:100 NEAA, 1 mM Pyruvate, 10 mM HEPES, 1:100 ITS, 10 μg/ml heparin sulfate, 1 mM N-acetyl cysteine, 10 μM zinc sulfate, 1.75 μL 2-mercaptoethanol, 10 nM T3, 1:2000 Trace elements A (Cellgro), 1:2000 trace elements B (Cellgro), 1:2000 Lipid Concentrate (Thermo Fisher Scientific) and 0.5 µM ZM447439 (Selleckchem) 55.
Dispersion/reaggregation of the stem cell-derived aggregates was carried out on Day 25 as previously described 74. Briefly, aggregates were dissociated with Accumax (Thermo Fisher Scientific) for 8 minutes at 37 **°**C with occasional shaking. Cells were washed with DPBS, centrifuged (5min/200xg) then resuspended in DPBS with 0.2% FBS and 10µM Y-27632. Cells were filtered through a strainer-capped FACS tube prior to FACS purification of the GFP+ fraction. FAC-sorted cells were collected in Stage 6.2 media containing 10µM Y-27632 and then aggregated in AggreWell 400 24-well plates (STEMCELL Technologies) using the manufacturer’s protocol at 1000 cells per aggregate.
View this table:
[Table1](http://medrxiv.org/content/early/2024/03/05/2023.05.24.23290298/T1)
Media formulations
### Proteomics
#### Library generation
We generated a pooled protein sample of 10 non-diabetic donors (as indicated in Table S1) at a concentration of 1.5ug/uL. Briefly, 100 μl of SDS lysis buffer (4 % SDS, 100 mM Tris, pH = 8) was added to each sample. Samples were then vortexed for 1 minute, boiled at 100°C for 10 minutes, vortexed again for 1 minute and then centrifuged at 4°C at 10,000 RPM for 10 minutes. The supernatants were collected and tested for protein concentration by BCA protein assay (ThermoFisher Scientific, Cat. # PI-23225). Measurements were performed on a Spark plate reader (TECAN). All protein samples were stored at –80°C.
Next, 297 μg of the pooled protein sample was used for further processing. The Trapped Ion Mobility – Time of Flight Mass Spectrometer (TimsTOF Pro; Bruker Daltonics, Germany) was set to Parallel Accumulation-Serial Fragmentation (PASEF) scan mode for DDA acquisition scanning 100 – 1700 m/z with 5 PASEF ramps. The capillary voltage was set to 1800V, drying gas to 3L/min, and drying temperature to 180°C. The MS and MS/MS spectra were acquired from m/z 100 – 1700. As for TIMS setting, ion mobility range (1/k0) was set to 0.70 – 1.35 V·s/cm2, 100ms ramp time and accumulation time (100% duty cycle), and ramp rate of 9.42Hz; this resulted in 0.64s of total cycle time.
Linear precursor repetitions were applied with a target intensity of 21,000 and 2500 intensity threshold. The active exclusion was enabled with a 0.4 min release. The collision energy was ramped linearly as a function of mobility from 27eV at 1/k0 = 0.7 V·s/cm2 to 55eV at 1/k0 = 1.35 V·s/cm2. Isolation widths were set at 2.07 m/z at < 400 m/z and 3.46 m/z at > 1000 m/z. Mass accuracy: error of mass measurement is typically within 3 ppm and is not allowed to exceed 7 ppm. For calibration of ion mobility dimension, the ions of Agilent ESI-Low Tuning Mix ions were selected (m/z [Th], 1/k0 [Th]: 622.0290, 0.9915; 922.0098, 1.1986; 1221.9906, 1.3934). TimsTOF Pro was run with timsControl v. 3.0.0 (Bruker). LC and MS were controlled with HyStar 6.0 (6.0.30.0, Bruker).
Acquired data were then searched using FragPipe75 computational platform (v. 17.1) with MSFragger76,77 (v. 3.4), Philosopher78 (v. 3.8), and EasyPQP (v. 0.1.27) components to build a spectral library. Protein sequence database Homo sapiens from Uniprot (reviewed sequences only; downloaded on January 15, 2021) and common contaminants protein, containing in total 20428 sequences were used, where reversed protein sequences were appended to the original database as decoys. For the MSFragger analysis, both precursor and fragment mass tolerances were set to 50 ppm. Enzyme specificity was set to “trypsin”, with up to 2 missed cleavages allowed.
For each analysis, the MS/MS search results were further processed using Philosopher, where final reports were generated and filtered at 1% protein FDR plus 1% PSM/ion/peptide-level FDR79 for each corresponding PSM.tsv, ion.tsv, and peptide.tsv files. Finally, .tsv files were used as input to EasyPQP for the generation of consensus spectrum libraries. The final spectral library was filtered to 1% protein and 1% peptide-level FDR.
#### Experimental samples
For each donor approximately 300 hand-picked islets were transferred into a 1.5 ml microcentrifuge tube. Culture media was removed, islets were washed twice in 1X PBS (Invitrogen, Cat. #10010049) and cell pellets were immediately flash frozen in liquid nitrogen and then stored at –80°C. Samples were processed and protein was extracted and measured as described above.
Next, 20 μg of each lysed sample was taken for further processing. Reduction of disulfide bonds was done by incubation with DTT (30 mM) for 35 min at 37°C, followed by alkylation by CAA (50 mM) for 20 min at 37°C. Samples were loaded onto 10% Mini-PROTEAN® TGX™ Precast Protein Gels (BioRad) and run for 30 min at 85V. Proteins were visualized by Coomassie blue (stained for 20 min). Lanes were cut out and destained in Ambic:EtOH (60:40), dehydrated and digested with Trypsin, first round 0.32ug/lane for 15 hours at 37°C, second round 0.128ug/lane for 2 hours, at 37°C (0.448ug in total/lane). Digestion was stopped with 10% FA, and samples were extracted with a series of extraction solutions, twice by 50% ACN, 50% of 0.1% TFA, and twice by 80% ACN, 20% of 0.1% TFA. Samples were then concentrated via vacuum centrifugation. Extracted peptide samples were then cleaned up via STAGE-tip purification, briefly: Resolubilized acidified sample was forced through a conditioned and equilibrated column with 8-12 mm of C18 packing, washed with 1% TFA twice, and eluted into clean plates by buffer containing 40% ACN, 0.1% TFA, then dried down.
Each sample was reconstituted in 0.5% ACN, 0.1% formic acid for LC-MS/MS analysis of 150 ng total on-column injections (with n = 3 technical replicates). The digest was separated using NanoElute UHPLC system (Bruker Daltonics) with Aurora Series Gen2 (CSI) analytical column (25 cm x 75 μm 1.6 μm FSC C18, with Gen2 nanoZero and CSI fitting; Ion Opticks, Parkville, Victoria, Australia) heated to 50°C and coupled to a Trapped Ion Mobility – Time of Flight mass spectrometer (timsTOF Pro; Bruker Daltonics, Germany) operated in Data-Independent Acquisition – Parallel Accumulation-Serial Fragmentation (DIA-PASEF) mode. A standard 60 min gradient was run from 2% B to 12% B over 30 min, then to 33% B from 30 to 60 min, then to 95% B over 0.5 min, and held at 95% B for 7.72 min. Before each run, the analytical column was conditioned with 4 column volumes of buffer A. Where buffer A consisted of 0.1% aqueous formic acid and 0.5 % acetonitrile in water, and buffer B consisted of 0.1% formic acid in 99.4 % acetonitrile. The NanoElute thermostat temperature was set at 7°C. The analysis was performed at 0.3 μL/min flow rate.
The TimsTOF Pro was set to Parallel Accumulation-Serial Fragmentation (PASEF) scan mode for DIA acquisition scanning 100 – 1700 m/z. The capillary voltage was set to 1800V, drying gas to 3L/min, and drying temperature to 180°C. The MS1 scan was followed by 17 consecutive PASEF ramps containing 22 non-overlapping 35 m/z isolation windows, covering the m/z range 319.5 – 1089.5 (more information in DIA windows). As for TIMS setting, ion mobility range (1/k0) was set to 0.70 – 1.35 V·s/cm2, 100ms ramp time and accumulation time (100% duty cycle), and ramp rate of 9.42 Hz; this resulted in 1.91s of total cycle time. The collision energy was ramped linearly as a function of mobility from 27eV at 1/k0 = 0.7 V·s/cm2 to 55eV at 1/k0 = 1.35 V·s/cm2. Mass accuracy: error of mass measurement is typically within 3 ppm and is not allowed to exceed 7 ppm. For calibration of ion mobility dimension, the ions of Agilent ESI-Low Tuning Mix ions were selected (m/z [Th], 1/k0 [Th]: 622.0290, 0.9915; 922.0098, 1.1986; 1221.9906, 1.3934). TimsTOF Pro was run with timsControl v. 3.0.0 (Bruker). LC and MS were controlled with HyStar 6.0 (6.0.30.0, Bruker).
Acquired diaPASEF data were then searched using FragPipe 78 computational platform (v. 17.1) with MSFragger 79,80 (v. 3.4), Philosopher 78 (v. 3.8), EasyPQP (v. 0.1.27) and DIA-NN (v. 1.8) to obtain DIA quantification, with use of the spectral library generated from the before mentioned highly fractionated sample. Quantification mode was set to “Any LC (high precision)”. All other settings were left default.
### RNAseq and Nanostring
Islets for RNA-seq were processed in two centers, initially in Oxford (n=49) with later donors being processed at Stanford (n = 47). For those processed in Oxford the methods have been described 81. Briefly, freshly isolated human islets were collected at the Alberta Diabetes Institute IsletCore ([www.isletcore.ca](http://www.isletcore.ca)) in Edmonton, Canada. Freshly isolated islets were processed for RNA and DNA extraction after 1-3 days in culture in CMRL media. RNA was extracted from human islets using Trizol (Ambion, UK or Sigma-Aldrich, Canada). To clean remaining media from the islets, samples were washed three times with phosphate-buffered saline (Sigma-Aldrich, UK). After the final cleaning step 1 mL Trizol was added to the cells. The cells were lysed by pipetting immediately to ensure rapid inhibition of RNase activity and incubated at room temperature for 10 min. Lysates were then transferred to clean 1.5 mL RNase-free centrifuge tubes (Applied Biosystems, UK). RNA quality (RIN score) was determined using an Agilent 2100 Bioanalyser (Agilent, UK), with a RIN score > 6 deemed acceptable for inclusion in the study. Samples were stored at −80 °C prior to sequencing. PolyA selected libraries were prepared from total RNA at the Oxford Genomics Centre using NEBNext ultra directional RNA library prep kit for Illumina with custom 8 bp indexes 82. Libraries were multiplexed (three samples per lane), clustered using TruSeq PE Cluster Kit v3, and paired-end sequenced (75bp) using Illumina TruSeq v3 chemistry on the Illumina HiSeq2000 platform. For the samples processed at Stanford, freshly isolated human islets were collected at the Alberta Diabetes Institute IsletCore in Edmonton, Canada, the same as the Oxford samples. The islets were picked to >95% purity and washed with PBS, then stored in 1 mL of Trizol (Ambion, UK). Islets in trizol were stored at –80 C until they were processed for RNA using the phenol-chloroform extraction method. The extracted RNA was resuspended in water and then dnased to remove contaminating DNA. RNA quality (RIN score) and concentration were determined using an Agilent 2100 Bioanalyzer (Agilent, US), with a RIN score > 5 deemed acceptable for inclusion in the study. Samples were stored at −80 °C prior to sequencing. PolyA selected cDNA libraries were prepared from total RNA and then sequenced using NovaSeq 6000 by Novogene. The raw sequencing data (fastq files) were checked first in regarding of the number of reads (>= 20 million per sample and extra sequencing were performed if necessary) and then in regarding of quality control, including quality score per base and length distribution using fastqc (v0.11.9). Reads were aligned to human genome reference (GRCh38) using STAR (v2.7.9a) with ENSEMBL gene annotations (v101). Gene expression levels were counted using featureCounts (v2.0.1) on exonic reads only. Differential expression was compared using the Wald test in DESeq2 (v1.26.0). *p* values were adjusted using the Benjamini and Hochberg method.
Fifty islets from each donor batch were used to assess islet quality by profiling expression of 132 human islet genes as described 83. Gene expression was measured using nCounter prep kits and nCounter SPRINT profiler according to manufacturer’s protocol (NanoString, USA).
### Co-expression analysis
Analysis was performed with both RNA-seq and proteomics data to assess relationships between ‘modules’ of omics features and donor metadata, including technical islet isolation parameters, clinical metadata, and functional outcomes. There were three main steps: data processing, network construction and functional characterization, and module-donor metadata correlation analysis.
Raw protein abundance matrices (134 donors, 8489 proteins; Uniprot IDs) were filtered to remove proteins with greater than 50% missing values, reducing the number of proteins to 7919. Next, the proteomics data were normalized by sample median, log transformed (base 10), and missing values were imputed with the *missForest* R package 84 using a random forest-based algorithm. The RNA-seq counts matrix (96 donors, 17 673 genes; Entrez IDs) was filtered based on abundance (15% with lowest mean counts removed), normalized with the ‘variance stabilizing transformation’ from the *DeSeq2* R package 85, and then filtered based on variance (15% with lowest variance removed). The RNA-seq data was measured in seven batches; batch effects were adjusted for using the ComBat method in the *sva* R package 86. One of the batches (*n* = 6 samples) had a sequencing depth of about 50% compared to the other samples, resulting in ∼5000 genes with zero counts that prevented effective batch effect adjustment. These six samples were removed from the dataset. Donor metadata included clinical metadata (age, sex, BMI, diabetes diagnosis, HbA1c), technical islet isolation metadata (culture time, cold ischemia time, digestion time, purity, and insulin content), and functional outcomes (AUC values from perifusion experiments after glucose, leucine, and oleate/palmitate exposures). Continuous variables with right-skew distributions were log-transformed (base 10). A constant value was added to the perifusion outcomes that had negative values to ensure all data were > 0 before log-transformation.
Two co-expression networks were constructed using the same methods, one for each omics type, using the *WGCNA* R package 87. First, the soft threshold parameter was selected by choosing the lowest value at which the R2 value of the scale free topology model fit did not substantially increase (threshold parameter = 8 and 20 for proteomics and RNA-seq respectively; selected via visual examination of an elbow plot). Next, pairwise distances between omics features were calculated by computing signed adjacency and topological overlap matrices using the soft threshold parameter and the ‘bicor’ correlation function. Finally, modules of co-expressed (or co-abundant) features were defined by first hierarchically clustering features based on the distance matrix, and then detecting modules of densely connected features using the dynamic tree cut algorithm. Modules were annotated with KEGG pathways and GO terms (biological process, cellular component, and molecular function) by performing overrepresentation analysis with the list of features in each module.
The pairwise partial correlation between module expression/abundance levels and donor metadata, accounting for the technical isolation parameters (purity, culture time, digestion time, and cold ischemia time), were calculated with the partial Pearson correlation function in the *ppcor* R package 88. Module expression/abundance levels were summarized with *WGCNA* ‘eigengenes’, which are the first principal component scores for each module, re-scaled so that positive and negative scores can be interpreted as higher and lower expression/abundance levels respectively. Correlation p-values were adjusted using the Benjamini-Hochberg method, with the FDR set to 5%.
### Multi-omics data analyses
RNAseq data was analyzed by DESeq2 R package (version 1.36.0)85. Six samples exhibiting large batch difference were excluded, leaving a total of 90 individuals (82 ND and 8 T2D) included in the analyses. Genes with low counts (<5 raw counts in more than 50% samples) were removed. Raw counts of the genes were normalized by variance stabilizing transformation (VST). T2D group was contrasted to ND group, including batch as a covariate in the default Wald test. Genes with Benjamini-Hochberg adjusted p-values < 0.05 were considered as differentially expressed.
For proteomics data analyses, proteins that were undetected in more than 50% samples were removed. Benjamini-Hochberg adjusted p-values < 0.05 of student t-test were used to identify differentially expressed proteins between the various identified groups.
To identify proteins associated with islet parameters, linear regression was conducted using log2 protein abundances of each protein as a dependent variable, each islet parameter as an independent variable, and diabetes status as a covariate to adjust for the effect of T2D. The coefficients of islet parameters and their FDR-adjusted p-values were reported. Similarly, to identify RNAs associated with islet parameters, the VST normalized gene counts were used as dependent variables instead. The batches of sequencing were included as additional covariates. Pearson correlations were also conducted instead of linear regression.
To assess the how well RNA expression reflects the abundance of the corresponding protein, we conducted Pearson correlation between proteins and RNAs. For across-gene correlation, mean protein abundances and mean RNA TPMs (transcript per million) were correlated. We used protein signal/protein length as a raw estimate of relative abundances across different proteins and TPM as an estimate of relative RNA levels across different genes. For within-gene across sample correlations, we correlated each gene’s VST normalized RNA counts with log2 protein abundances, and the p values were adjusted by Benjamini-Hochberg. Instead of Pearson correlation, we also checked within-gene association using linear regression adjusting for batches of sequencing, and the results were mostly the same (not shown).
Gene set enrichment analysis (GSEA) was performed on the transcriptomic or proteomics changes (T2D/ND) using KEGG or Hallmark gene set libraries. The direction-signed –log10(p values) were used as rank scores. Over representation analysis (ORA) was performed on genes with positive, negative or non-significant RNA-protein correlations using Gene Ontology – Cellular Component (GO-CC) gene set library. Genes that are commonly detected in RNAseq and proteomics were used as background gene list. Both GSEA and ORA were conducted in clusterProfiler R package89.
To evaluate the concordant changes of RNAs and proteins in T2D, in addition to comparing differentially expressed RNAs and proteins, we applied the threshold-free method called rank-rank hypergeometric overlap test19. RNAs and proteins were ranked by direction-signed –log10(p value) or log2 fold change (T2D/ND). The ranked lists were processed in the online tool ([https://systems.crump.ucla.edu/rankrank/](https://systems.crump.ucla.edu/rankrank/)) and RRHO R package using a step size of 100.
Cell type deconvolution was conducted using the BisqueRNA R package using a marker-based decomposition approach. RNAseq DESeq2 normalized counts (without VST) or protein abundances (without log2 transformation) were used as the input. Islet cell type markers (alpha, beta, gamma, delta) were from van Gurp et al. 90, and acinar cell markers were from Segerstolpe 91. All code can be found at Github repository [https://github.com/hcen/Kolic_macronutrients](https://github.com/hcen/Kolic_macronutrients).
### General statistical analyses and data visualization
Statistical analyses and data presentation for perifusion studies of dynamic insulin secretion were carried out using GraphPad Prism 9 (Graphpad Software, San Diego, CA, USA) or R (v 4.1.1). using Student’s *t*-test for parametric data and Mann-Whitney test for non-parametric data. For all statistical analyses, differences were considered significant if the p-value was less than 0.05. *: p< 0.05; ** p< 0.01; \***| p< 0.001. Data were presented as means ± SEM with individual data points from biological replicates (average of two technical replicates).
## AUTHOR CONTRIBUTIONS
Conceptualization: JK, JX, FCL, ALG, LJF, PEM, JDJ
Methodology: JK, HHC, JE, HS, VR, SS, AFG, JEM, JL, LB, JX, FCL, ALG, LJF, PEM, JDJ
Investigation: JK, RM, JCR, HS, SS, FCL, ALG
Visualization: JK, WGS, HC, JE, YHX, JDJ
Funding acquisition: JX, FCL, ALG, LJF, PEM, JDJ
Supervision: JX, FCL, ALG, LJF, PEM, JDJ
Writing – original draft: JK, JDJ
## COMPETING INTERESTS
ALG declares that her spouse holds stock options in Roche and is an employee of Genentech. Other authors declare that they have no competing interests.
## Supporting information
Supplemental Figures [[supplements/290298_file03.pdf]](pending:yes)
Table S1 [[supplements/290298_file04.xlsx]](pending:yes)
Table S1 [[supplements/290298_file05.xlsx]](pending:yes)
Table S3 [[supplements/290298_file06.xlsx]](pending:yes)
Table S4 [[supplements/290298_file07.xlsx]](pending:yes)
Table S5 [[supplements/290298_file08.xlsx]](pending:yes)
Table S6 [[supplements/290298_file09.xlsx]](pending:yes)
Table S7 [[supplements/290298_file10.xlsx]](pending:yes)
Table S8 [[supplements/290298_file11.xlsx]](pending:yes)
Table S9 [[supplements/290298_file12.xlsx]](pending:yes)
Table S10 [[supplements/290298_file13.xlsx]](pending:yes)
Table S11 [[supplements/290298_file14.xlsx]](pending:yes)
## Data Availability
All data produced in the present study are available upon reasonable request to the authors.
## FUNDING
This work was supported by a Canadian Institutes for Health Research (CIHR) operating grant (168857) to J.D.J., a CIHR Team Grant (ASD-179092/5-SRA-2021-1149-S-B) to P.E.M., F.C.L., J.X., L.J.F., and J.D.J, a CIHR-JDRF Team (ASD-173663/5-SRA-2020-1059-S-B) to F.C.L., P.E.M., and J.D.J, and the JDRF Centre of Excellence at UBC (3-COE-2022-1103-M-B). This work was funded by Wellcome Trust grants to A.L.G at Oxford (095101, 200837, 106130, 203141, NIH (U01-DK105535; U01-DK085545) and by the Wellcome Trust (200837) and NIH (UM-1DK126185) at Stanford. Proteomics infrastructure and analysis used in this study was supported by the UBC Life Sciences Institute, Canada Foundation for Innovation, BC Knowledge Development Fund, and Genome Canada/BC (PRO264). J.K. was individually supported by a Banting Fellowship from CIHR.
## ACKNOWLEDGMENTS
We thank the Human Organ Procurement and Exchange (HOPE) program and Trillium Gift of Life Network (TGLN) for their work procuring human donor pancreas for research. We especially thank the organ donors and their families for their gift in support of diabetes research. The authors acknowledge that work done at UBC and BC Children’s Hospital occurs on the traditional, ancestral, and unceded territories of the Coast Salish peoples, the Sḵwx̱ wú7mesh (Squamish), səl̓ ilwətaɁɬ (Tsleil-Waututh), and xʷməθkʷəy̓ əm (Musqueam) Nations.
## Footnotes
* Methods section has been updated. Multiple new sets of analysis and figures have been added.
* Received May 24, 2023.
* Revision received March 4, 2024.
* Accepted March 5, 2024.
* © 2024, Posted by Cold Spring Harbor Laboratory
The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.
## REFERENCES
1. 1.Kayton, N.S., Poffenberger, G., Henske, J., Dai, C., Thompson, C., Aramandla, R., Shostak, A., Nicholson, W., Brissova, M., Bush, W.S., and Powers, A.C. (2015). Human islet preparations distributed for research exhibit a variety of insulin-secretory profiles. Am J Physiol Endocrinol Metab 308, E592–602. doi:10.1152/ajpendo.00437.2014.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/ajpendo.00437.2014&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=25648831&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
2. 2.Lyon, J., Manning Fox, J.E., Spigelman, A.F., Kim, R., Smith, N., O’Gorman, D., Kin, T., Shapiro, A.M., Rajotte, R.V., and MacDonald, P.E. (2016). Research-Focused Isolation of Human Islets From Donors With and Without Diabetes at the Alberta Diabetes Institute IsletCore. Endocrinology 157, 560–569. doi:10.1210/en.2015-1562.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/en.2015-1562&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=26653569&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
3. 3.Alcazar, O., and Buchwald, P. (2019). Concentration-Dependency and Time Profile of Insulin Secretion: Dynamic Perifusion Studies With Human and Murine Islets. Front Endocrinol (Lausanne) 10, 680. doi:10.3389/fendo.2019.00680.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fendo.2019.00680&link_type=DOI)
4. 4.Cen, J., Sargsyan, E., and Bergsten, P. (2016). Fatty acids stimulate insulin secretion from human pancreatic islets at fasting glucose concentrations via mitochondria-dependent and – independent mechanisms. Nutr Metab (Lond) 13, 59. doi:10.1186/s12986-016-0119-5.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12986-016-0119-5&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=27582778&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
5. 5.Henquin, J.C., Dufrane, D., and Nenquin, M. (2006). Nutrient control of insulin secretion in isolated normal human islets. Diabetes 55, 3470–3477. doi:10.2337/db06-0868.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6MTA6IjU1LzEyLzM0NzAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wMy8wNS8yMDIzLjA1LjI0LjIzMjkwMjk4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)
6. 6.Zhang, A.M.Y., Wellberg, E.A., Kopp, J.L., and Johnson, J.D. (2021). Hyperinsulinemia in Obesity, Inflammation, and Cancer. Diabetes Metab J 45, 285–311. doi:10.4093/dmj.2020.0250.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4093/dmj.2020.0250&link_type=DOI)
7. 7.Wang, P., Song, M., Eliassen, A.H., Wang, M., Fung, T.T., Clinton, S.K., Rimm, E.B., Hu, F.B., Willett, W.C., Tabung, F.K., and Giovannucci, E.L. (2023). Optimal dietary patterns for prevention of chronic disease. Nature Medicine. doi:10.1038/s41591-023-02235-5.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-023-02235-5&link_type=DOI)
8. 8.Berry, S.E., Valdes, A.M., Drew, D.A., Asnicar, F., Mazidi, M., Wolf, J., Capdevila, J., Hadjigeorgiou, G., Davies, R., Al Khatib, H., et al. (2020). Human postprandial responses to food and potential for precision nutrition. Nat Med 26, 964–973. doi:10.1038/s41591-020-0934-0.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0934-0&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=32528151&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
9. 9.Zeevi, D., Korem, T., Zmora, N., Israeli, D., Rothschild, D., Weinberger, A., Ben-Yacov, O., Lador, D., Avnit-Sagi, T., Lotan-Pompan, M., et al. (2015). Personalized Nutrition by Prediction of Glycemic Responses. Cell 163, 1079–1094. doi:10.1016/j.cell.2015.11.001.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2015.11.001&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=26590418&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
10. 10.Canada, S. (2022). Population estimates on July 1st, by age and sex. [https://www150.statcan.gc.ca/t1/tbl1/en/tv.action?pid=1710000501](https://www150.statcan.gc.ca/t1/tbl1/en/tv.action?pid=1710000501).
11. 11.Henquin, J.C., Dufrane, D., Kerr-Conte, J., and Nenquin, M. (2015). Dynamics of glucose-induced insulin secretion in normal human islets. Am J Physiol Endocrinol Metab 309, E640–650. doi:10.1152/ajpendo.00251.2015.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/ajpendo.00251.2015&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=26264556&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
12. 12.Basturk, B., Koc Ozerson, Z., and Yuksel, A. (2021). Evaluation of the Effect of Macronutrients Combination on Blood Sugar Levels in Healthy Individuals. Iran J Public Health 50, 280–287. doi:10.18502/ijph.v50i2.5340.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18502/ijph.v50i2.5340&link_type=DOI)
13. 13.Shukla, A.P., Iliescu, R.G., Thomas, C.E., and Aronne, L.J. (2015). Food Order Has a Significant Impact on Postprandial Glucose and Insulin Levels. Diabetes Care 38, e98–99. doi:10.2337/dc15-0429.
[FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGlhY2FyZSI7czo1OiJyZXNpZCI7czo4OiIzOC83L2U5OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAzLzA1LzIwMjMuMDUuMjQuMjMyOTAyOTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)
14. 14.Gravena, C., Mathias, P.C., and Ashcroft, S.J. (2002). Acute effects of fatty acids on insulin secretion from rat and human islets of Langerhans. J Endocrinol 173, 73–80. doi:10.1677/joe.0.1730073.
[Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoiam9lIjtzOjU6InJlc2lkIjtzOjg6IjE3My8xLzczIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDMvMDUvMjAyMy4wNS4yNC4yMzI5MDI5OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)
15. 15.Ferdaoussi, M., Bergeron, V., Zarrouki, B., Kolic, J., Cantley, J., Fielitz, J., Olson, E.N., Prentki, M., Biden, T., MacDonald, P.E., and Poitout, V. (2012). G protein-coupled receptor (GPR)40-dependent potentiation of insulin secretion in mouse islets is mediated by protein kinase D1. Diabetologia 55, 2682–2692. doi:10.1007/s00125-012-2650-x.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-012-2650-x&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=22820510&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
16. 16.Kreutzberger, A.J.B., Kiessling, V., Doyle, C.A., Schenk, N., Upchurch, C.M., Elmer-Dixon, M., Ward, A.E., Preobraschenski, J., Hussein, S.S., Tomaka, W., et al. (2020). Distinct insulin granule subpopulations implicated in the secretory pathology of diabetes types 1 and 2. Elife 9. doi:10.7554/eLife.62506.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.62506&link_type=DOI)
17. 17.Manders, R.J., Koopman, R., Sluijsmans, W.E., van den Berg, R., Verbeek, K., Saris, W.H., Wagenmakers, A.J., and van Loon, L.J. (2006). Co-ingestion of a protein hydrolysate with or without additional leucine effectively reduces postprandial blood glucose excursions in Type 2 diabetic men. J Nutr 136, 1294–1299. doi:10.1093/jn/136.5.1294.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjEwOiIxMzYvNS8xMjk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDMvMDUvMjAyMy4wNS4yNC4yMzI5MDI5OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)
18. 18.Segerstolpe, A., Palasantza, A., Eliasson, P., Andersson, E.M., Andreasson, A.C., Sun, X., Picelli, S., Sabirsh, A., Clausen, M., Bjursell, M.K., et al. (2016). Single-Cell Transcriptome Profiling of Human Pancreatic Islets in Health and Type 2 Diabetes. Cell Metab 24, 593–607. doi:10.1016/j.cmet.2016.08.020.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2016.08.020&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=27667667&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
19. 19.Plaisier, S.B., Taschereau, R., Wong, J.A., and Graeber, T.G. (2010). Rank-rank hypergeometric overlap: identification of statistically significant overlap between gene-expression signatures. Nucleic Acids Res 38, e169. doi:10.1093/nar/gkq636.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkq636&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=20660011&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
20. 20.Xu, J., Han, J., Long, Y.S., Epstein, P.N., and Liu, Y.Q. (2008). The role of pyruvate carboxylase in insulin secretion and proliferation in rat pancreatic beta cells. Diabetologia 51, 2022–2030. doi:10.1007/s00125-008-1130-9.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-008-1130-9&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=18769905&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259816500012&link_type=ISI)
21. 21.Dong, G., Adak, S., Spyropoulos, G., Zhang, Q., Feng, C., Yin, L., Speck, S.L., Shyr, Z., Morikawa, S., Kitamura, R.A., et al. (2023). Palmitoylation couples insulin hypersecretion with beta cell failure in diabetes. Cell Metab 35, 332–344 e337. doi:10.1016/j.cmet.2022.12.012.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2022.12.012&link_type=DOI)
22. 22.Jaikaran, E.T., and Clark, A. (2001). Islet amyloid and type 2 diabetes: from molecular misfolding to islet pathophysiology. Biochim Biophys Acta 1537, 179–203. doi:10.1016/s0925-4439(01)00078-3.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0925-4439(01)00078-3&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=11731221&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
23. 23.Simpson, F., Peden, A.A., Christopoulou, L., and Robinson, M.S. (1997). Characterization of the adaptor-related protein complex, AP-3. J Cell Biol 137, 835–845. doi:10.1083/jcb.137.4.835.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNiIjtzOjU6InJlc2lkIjtzOjk6IjEzNy80LzgzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAzLzA1LzIwMjMuMDUuMjQuMjMyOTAyOTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)
24. 24.Collins, S.C., Do, H.W., Hastoy, B., Hugill, A., Adam, J., Chibalina, M.V., Galvanovskis, J., Godazgar, M., Lee, S., Goldsworthy, M., et al. (2016). Increased Expression of the Diabetes Gene SOX4 Reduces Insulin Secretion by Impaired Fusion Pore Expansion. Diabetes 65, 1952–1961. doi:10.2337/db15-1489.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjUvNy8xOTUyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDMvMDUvMjAyMy4wNS4yNC4yMzI5MDI5OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)
25. 25.Sladek, R. (2023). Predicting the response to GLP-1 receptor agonists: an unexpected role for beta-arrestin-1. Lancet Diabetes Endocrinol 11, 3–4. doi:10.1016/S2213-8587(22)00357-6.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2213-8587(22)00357-6&link_type=DOI)
26. 26.Wegler, C., Olander, M., Wisniewski, J.R., Lundquist, P., Zettl, K., Asberg, A., Hjelmesaeth, J., Andersson, T.B., and Artursson, P. (2020). Global variability analysis of mRNA and protein concentrations across and within human tissues. NAR Genom Bioinform 2, lqz010. doi:10.1093/nargab/lqz010.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nargab/lqz010&link_type=DOI)
27. 27.Jiang, L., Wang, M., Lin, S., Jian, R., Li, X., Chan, J., Dong, G., Fang, H., Robinson, A.E., Consortium, G.T., and Snyder, M.P. (2020). A Quantitative Proteome Map of the Human Body. Cell 183, 269–283 e219. doi:10.1016/j.cell.2020.08.036.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.08.036&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=32916130&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
28. 28.Szklarczyk, D., Gable, A.L., Nastou, K.C., Lyon, D., Kirsch, R., Pyysalo, S., Doncheva, N.T., Legeay, M., Fang, T., Bork, P., et al. (2021). The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res 49, D605–D612. doi:10.1093/nar/gkaa1074.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkaa1074&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=33237311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
29. 29.Warnotte, C., Nenquin, M., and Henquin, J.C. (1999). Unbound rather than total concentration and saturation rather than unsaturation determine the potency of fatty acids on insulin secretion. Mol Cell Endocrinol 153, 147–153. doi:10.1016/s0303-7207(99)00069-6.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0303-7207(99)00069-6&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=10459862&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000081494900015&link_type=ISI)
30. 30.Cantley, J., Davenport, A., Vetterli, L., Nemes, N.J., Whitworth, P.T., Boslem, E., Thai, L.M., Mellett, N., Meikle, P.J., Hoehn, K.L., et al. (2019). Disruption of beta cell acetyl-CoA carboxylase-1 in mice impairs insulin secretion and beta cell mass. Diabetologia 62, 99–111. doi:10.1007/s00125-018-4743-7.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-018-4743-7&link_type=DOI)
31. 31.Mo, H., Liu, Y., Zhang, M., Qiu, Z., Li, Y., Zhang, Z., Li, Y., and Xu, G. (2022). The Role of Pancreatic Fatty Acid Synthesis in Islet Morphology and Function after Caloric Restriction or Roux-En-Y Gastric Bypass Surgery in Mice. Genes (Basel) 14. doi:10.3390/genes14010005.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/genes14010005&link_type=DOI)
32. 32.Roomp, K., Kristinsson, H., Schvartz, D., Ubhayasekera, K., Sargsyan, E., Manukyan, L., Chowdhury, A., Manell, H., Satagopam, V., Groebe, K., et al. (2017). Combined lipidomic and proteomic analysis of isolated human islets exposed to palmitate reveals time-dependent changes in insulin secretion and lipid metabolism. PLoS One 12, e0176391. doi:10.1371/journal.pone.0176391.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0176391&link_type=DOI)
33. 33.Mondal, A., Appu, A.P., Sadhukhan, T., Bagh, M.B., Previde, R.M., Sadhukhan, S., Stojilkovic, S., Liu, A., and Mukherjee, A.B. (2022). Ppt1-deficiency dysregulates lysosomal Ca(++) homeostasis contributing to pathogenesis in a mouse model of CLN1 disease. J Inherit Metab Dis 45, 635–656. doi:10.1002/jimd.12485.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jimd.12485&link_type=DOI)
34. 34.Zou, Y., Pei, J., Wang, Y., Chen, Q., Sun, M., Kang, L., Zhang, X., Zhang, L., Gao, X., and Lin, Z. (2022). The Deficiency of SCARB2/LIMP-2 Impairs Metabolism via Disrupted mTORC1-Dependent Mitochondrial OXPHOS. Int J Mol Sci 23. doi:10.3390/ijms23158634.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijms23158634&link_type=DOI)
35. 35.Munk, R., Anerillas, C., Rossi, M., Tsitsipatis, D., Martindale, J.L., Herman, A.B., Yang, J.H., Roberts, J.A., Varma, V.R., Pandey, P.R., et al. (2021). Acid ceramidase promotes senescent cell survival. Aging (Albany NY) 13, 15750–15769. doi:10.18632/aging.203170.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18632/aging.203170&link_type=DOI)
36. 36.Suurvali, J., Boudinot, P., Kanellopoulos, J., and Ruutel Boudinot, S. (2017). P2X4: A fast and sensitive purinergic receptor. Biomed J 40, 245–256. doi:10.1016/j.bj.2017.06.010.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bj.2017.06.010&link_type=DOI)
37. 37.Zhang, J., Kong, X., Zhang, Y., Sun, W., Wang, J., Chen, M., and Chen, X. (2020). FDXR regulates TP73 tumor suppressor via IRP2 to modulate aging and tumor suppression. J Pathol 251, 284–296. doi:10.1002/path.5451.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/path.5451&link_type=DOI)
38. 38.Wang, Y., Zhu, H., and Wang, X. (2022). Prognosis and immune infiltration analysis of endoplasmic reticulum stress-related genes in bladder urothelial carcinoma. Front Genet 13, 965100. doi:10.3389/fgene.2022.965100.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fgene.2022.965100&link_type=DOI)
39. 39.Wu, P., Mo, Y., Peng, M., Tang, T., Zhong, Y., Deng, X., Xiong, F., Guo, C., Wu, X., Li, Y., et al. (2020). Emerging role of tumor-related functional peptides encoded by lncRNA and circRNA. Mol Cancer 19, 22. doi:10.1186/s12943-020-1147-3.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12943-020-1147-3&link_type=DOI)
40. 40.Wang, X., Zhou, T., Yang, X., Cao, X., Jin, G., Zhang, P., Guo, J., Rong, K., Li, B., Hu, Y., et al. (2023). DDRGK1 Enhances Osteosarcoma Chemoresistance via Inhibiting KEAP1-Mediated NRF2 Ubiquitination. Adv Sci (Weinh) 10, e2204438. doi:10.1002/advs.202204438.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/advs.202204438&link_type=DOI)
41. 41.Zhang, B., Zhao, Z., Wang, Y., Guo, T., He, M., Guan, G., Peng, P., Cai, J., Zhang, B., Liu, X., and Song, Q. (2023). STX5 Inhibits Hepatocellular Carcinoma Adhesion and Promotes Metastasis by Regulating the PI3K/mTOR Pathway. J Clin Transl Hepatol 11, 572–583. doi:10.14218/JCTH.2022.00276.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.14218/JCTH.2022.00276&link_type=DOI)
42. 42.Hope, C., Robertshaw, A., Cheung, K.L., Idris, I., and English, E. (2016). Relationship between HbA1c and cancer in people with or without diabetes: a systematic review. Diabet Med 33, 1013–1025. doi:10.1111/dme.13031.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/dme.13031&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
43. 43.Jezek, P., Holendova, B., Jaburek, M., Dlaskova, A., and Plecita-Hlavata, L. (2022). Contribution of Mitochondria to Insulin Secretion by Various Secretagogues. Antioxid Redox Signal 36, 920–952. doi:10.1089/ars.2021.0113.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/ars.2021.0113&link_type=DOI)
44. 44.Courtney, M., Gjernes, E., Druelle, N., Ravaud, C., Vieira, A., Ben-Othman, N., Pfeifer, A., Avolio, F., Leuckx, G., Lacas-Gervais, S., et al. (2013). The inactivation of Arx in pancreatic alpha-cells triggers their neogenesis and conversion into functional beta-like cells. PLoS Genet 9, e1003934. doi:10.1371/journal.pgen.1003934.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1003934&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=24204325&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
45. 45.Katsumoto, K., Yennek, S., Chen, C., Silva, L.F.D., Traikov, S., Sever, D., Azad, A., Shan, J., Vainio, S., Ninov, N., et al. (2022). Wnt4 is heterogeneously activated in maturing beta-cells to control calcium signaling, metabolism and function. Nat Commun 13, 6255. doi:10.1038/s41467-022-33841-5.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-022-33841-5&link_type=DOI)
46. 46. Liis Kolberg, U.R., Ivan Kuzmin, Jaak Vilo, Hedi Peterson (2020). gprofiler2 –-an R package for gene list functional enrichment analysis and namespace conversion toolset g:Profiler [version 2; peer review: 2 approved]. F1000Research 2020 9(ELIXIR).
47. 47.Lai, E., Bikopoulos, G., Wheeler, M.B., Rozakis-Adcock, M., and Volchuk, A. (2008). Differential activation of ER stress and apoptosis in response to chronically elevated free fatty acids in pancreatic beta-cells. Am J Physiol Endocrinol Metab 294, E540–550. doi:10.1152/ajpendo.00478.2007.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/ajpendo.00478.2007&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=18198352&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253816000008&link_type=ISI)
48. 48.Leung, P.S. (2003). Pancreatic renin-angiotensin system: a novel target for the potential treatment of pancreatic diseases? JOP 4, 89–91.
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=12629265&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
49. 49.Kolic, J., and Johnson, J.D. (2021). Promises and pitfalls of beta cell-replacement therapies. Nat Metab 3, 1036–1037. doi:10.1038/s42255-021-00433-4.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s42255-021-00433-4&link_type=DOI)
50. 50.Barsby, T., Vahakangas, E., Ustinov, J., Montaser, H., Ibrahim, H., Lithovius, V., Kuuluvainen, E., Chandra, V., Saarimaki-Vire, J., Katajisto, P., et al. (2023). Aberrant metabolite trafficking and fuel sensitivity in human pluripotent stem cell-derived islets. Cell Rep 42, 112970. doi:10.1016/j.celrep.2023.112970.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2023.112970&link_type=DOI)
51. 51.Novakovsky, G., Sasaki, S., Fornes, O., Omur, M.E., Huang, H., Bayly, C.L., Zhang, D., Lim, N., Cherkasov, A., Pavlidis, P., et al. (2023). In silico discovery of small molecules for efficient stem cell differentiation into definitive endoderm. Stem Cell Reports 18, 765–781. doi:10.1016/j.stemcr.2023.01.008.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.stemcr.2023.01.008&link_type=DOI)
52. 52.Jacobson, D.A., and Shyng, S.L. (2020). Ion Channels of the Islets in Type 2 Diabetes. J Mol Biol 432, 1326–1346. doi:10.1016/j.jmb.2019.08.014.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jmb.2019.08.014&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=31473158&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
53. 53.Puri, S., Roy, N., Russ, H.A., Leonhardt, L., French, E.K., Roy, R., Bengtsson, H., Scott, D.K., Stewart, A.F., and Hebrok, M. (2018). Replication confers beta cell immaturity. Nat Commun 9, 485. doi:10.1038/s41467-018-02939-0.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-02939-0&link_type=DOI)
54. 54.Barsby, T., and Otonkoski, T. (2022). Maturation of beta cells: lessons from in vivo and in vitro models. Diabetologia 65, 917–930. doi:10.1007/s00125-022-05672-y.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-022-05672-y&link_type=DOI)
55. 55.Balboa, D., Barsby, T., Lithovius, V., Saarimaki-Vire, J., Omar-Hmeadi, M., Dyachok, O., Montaser, H., Lund, P.E., Yang, M., Ibrahim, H., et al. (2022). Functional, metabolic and transcriptional maturation of human pancreatic islets derived from stem cells. Nat Biotechnol 40, 1042–1055. doi:10.1038/s41587-022-01219-z.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41587-022-01219-z&link_type=DOI)
56. 56.Maxwell, K.G., Kim, M.H., Gale, S.E., and Millman, J.R. (2022). Differential Function and Maturation of Human Stem Cell-Derived Islets After Transplantation. Stem Cells Transl Med 11, 322–331. doi:10.1093/stcltm/szab013.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/stcltm/szab013&link_type=DOI)
57. 57.Malaeb, S., Bakker, C., Chow, L.S., and Bantle, A.E. (2019). High-Protein Diets for Treatment of Type 2 Diabetes Mellitus: A Systematic Review. Adv Nutr 10, 621–633. doi:10.1093/advances/nmz002.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/advances/nmz002&link_type=DOI)
58. 58.Diamanti, K., Cavalli, M., Pereira, M.J., Pan, G., Castillejo-López, C., Kumar, C., Mundt, F., Komorowski, J., Deshmukh, A.S., Mann, M., et al. (2022). Organ-specific metabolic pathways distinguish prediabetes, type 2 diabetes, and normal tissues. Cell Rep Med 3, 100763. doi:10.1016/j.xcrm.2022.100763.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.xcrm.2022.100763&link_type=DOI)
59. 59.Pigeau, G.M., Kolic, J., Ball, B.J., Hoppa, M.B., Wang, Y.W., Ruckle, T., Woo, M., Manning Fox, J.E., and MacDonald, P.E. (2009). Insulin granule recruitment and exocytosis is dependent on p110gamma in insulinoma and human beta-cells. Diabetes 58, 2084–2092. doi:10.2337/db08-1371.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNTgvOS8yMDg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDMvMDUvMjAyMy4wNS4yNC4yMzI5MDI5OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)
60. 60.Vega-Lopez, S., Ausman, L.M., Griffith, J.L., and Lichtenstein, A.H. (2007). Interindividual variability and intra-individual reproducibility of glycemic index values for commercial white bread. Diabetes Care 30, 1412–1417. doi:10.2337/dc06-1598.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGlhY2FyZSI7czo1OiJyZXNpZCI7czo5OiIzMC82LzE0MTIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wMy8wNS8yMDIzLjA1LjI0LjIzMjkwMjk4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)
61. 61.Rein, M., Ben-Yacov, O., Godneva, A., Shilo, S., Zmora, N., Kolobkov, D., Cohen-Dolev, N., Wolf, B.C., Kosower, N., Lotan-Pompan, M., et al. (2022). Effects of personalized diets by prediction of glycemic responses on glycemic control and metabolic health in newly diagnosed T2DM: a randomized dietary intervention pilot trial. BMC Med 20, 56. doi:10.1186/s12916-022-02254-y.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12916-022-02254-y&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
62. 62.Popp, C.J., Hu, L., Kharmats, A.Y., Curran, M., Berube, L., Wang, C., Pompeii, M.L., Illiano, P., St-Jules, D.E., Mottern, M., et al. (2022). Effect of a Personalized Diet to Reduce Postprandial Glycemic Response vs a Low-fat Diet on Weight Loss in Adults With Abnormal Glucose Metabolism and Obesity: A Randomized Clinical Trial. JAMA Netw Open 5, e2233760. doi:10.1001/jamanetworkopen.2022.33760.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2022.33760&link_type=DOI)
63. 63.Watanabe, K., Wilmanski, T., Diener, C., Earls, J.C., Zimmer, A., Lincoln, B., Hadlock, J.J., Lovejoy, J.C., Gibbons, S.M., Magis, A.T., et al. (2023). Multiomic signatures of body mass index identify heterogeneous health phenotypes and responses to a lifestyle intervention. Nat Med. doi:10.1038/s41591-023-02248-0.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-023-02248-0&link_type=DOI)
64. 64. Dos Santos, K.C., Olofsson, C., Cunha, J., Roberts, F., Catrina, S.B., Fex, M., Ekberg, N.R., and Spegel, P. (2022). The impact of macronutrient composition on metabolic regulation: An Islet-Centric view. Acta Physiol (Oxf) 236, e13884. doi:10.1111/apha.13884.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/apha.13884&link_type=DOI)
65. 65.Thiebaud, D., DeFronzo, R.A., Jacot, E., Golay, A., Acheson, K., Maeder, E., Jequier, E., and Felber, J.P. (1982). Effect of long chain triglyceride infusion on glucose metabolism in man. Metabolism 31, 1128–1136. doi:10.1016/0026-0495(82)90163-9.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0026-0495(82)90163-9&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=6752642&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1982PQ13100010&link_type=ISI)
66. 66.Jensen, C.B., Storgaard, H., Holst, J.J., Dela, F., Madsbad, S., and Vaag, A.A. (2003). Insulin secretion and cellular glucose metabolism after prolonged low-grade intralipid infusion in young men. J Clin Endocrinol Metab 88, 2775–2783. doi:10.1210/jc.2002-021430.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2002-021430&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=12788887&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183318200060&link_type=ISI)
67. 67.Trico, D., Mengozzi, A., Baldi, S., Bizzotto, R., Olaniru, O., Toczyska, K., Huang, G.C., Seghieri, M., Frascerra, S., Amiel, S.A., et al. (2022). Lipid-induced glucose intolerance is driven by impaired glucose kinetics and insulin metabolism in healthy individuals. Metabolism 134, 155247. doi:10.1016/j.metabol.2022.155247.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.metabol.2022.155247&link_type=DOI)
68. 68.Storgaard, H., Jensen, C.B., Vaag, A.A., Volund, A., and Madsbad, S. (2003). Insulin secretion after short– and long-term low-grade free fatty acid infusion in men with increased risk of developing type 2 diabetes. Metabolism 52, 885–894. doi:10.1016/s0026-0495(03)00102-1.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0026-0495(03)00102-1&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=12870166&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184197500015&link_type=ISI)
69. 69.Kashyap, S., Belfort, R., Gastaldelli, A., Pratipanawatr, T., Berria, R., Pratipanawatr, W., Bajaj, M., Mandarino, L., DeFronzo, R., and Cusi, K. (2003). A sustained increase in plasma free fatty acids impairs insulin secretion in nondiabetic subjects genetically predisposed to develop type 2 diabetes. Diabetes 52, 2461–2474. doi:10.2337/diabetes.52.10.2461.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6MTA6IjUyLzEwLzI0NjEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wMy8wNS8yMDIzLjA1LjI0LjIzMjkwMjk4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)
70. 70.Michaliszyn, S.F., Bonadonna, R.C., Sjaarda, L.A., Lee, S., Farchoukh, L., and Arslanian, S.A. (2013). beta-Cell lipotoxicity in response to free fatty acid elevation in prepubertal youth: African American versus Caucasian contrast. Diabetes 62, 2917–2922. doi:10.2337/db12-1664.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjIvOC8yOTE3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDMvMDUvMjAyMy4wNS4yNC4yMzI5MDI5OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)
71. 71.Lyon, J., Spigelman, A.F., Manning Fox, J.E., and MacDonald, P.E. (2021). ADI Isletcore protocols for the isolation, assessment and cryopreservation of human pancreatic islets of langerhans for research purposes. Protocols.io 2. dx.doi.org/10.17504/protocols.io.bupanvie.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.17504/protocols.io.bupanvie&link_type=DOI)
72. 72.Kolic, J., Beet, L., Overby, P., Cen, H.H., Panzhinskiy, E., Ure, D.R., Cross, J.L., Huizinga, R.B., and Johnson, J.D. (2020). Differential Effects of Voclosporin and Tacrolimus on Insulin Secretion From Human Islets. Endocrinology 161. doi:10.1210/endocr/bqaa162.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/endocr/bqaa162&link_type=DOI)
73. 73.Castro-Gutierrez, R., Alkanani, A., Mathews, C.E., Michels, A., and Russ, H.A. (2021). Protecting Stem Cell Derived Pancreatic Beta-Like Cells From Diabetogenic T Cell Recognition. Front Endocrinol (Lausanne) 12, 707881. doi:10.3389/fendo.2021.707881.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fendo.2021.707881&link_type=DOI)
74. 74.Mar, S., Filatov, E., Nian, C., Sasaki, S., Zhang, D., and Lynn, F.C. (2023). Tracking insulin– and glucagon-expressing bihormonal cells during differentiation using an INSULIN and GLUCAGON double reporter human embryonic stem cell line. bioRxiv, 2023.2004.2019.537542. doi:10.1101/2023.04.19.537542.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMy4wNC4xOS41Mzc1NDJ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAzLzA1LzIwMjMuMDUuMjQuMjMyOTAyOTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)
75. 75.Demichev, V., Messner, C.B., Vernardis, S.I., Lilley, K.S., and Ralser, M. (2020). DIA-NN: neural networks and interference correction enable deep proteome coverage in high throughput. Nat Methods 17, 41–44. doi:10.1038/s41592-019-0638-x.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41592-019-0638-x&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=31768060&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
76. 76.Kong, A.T., Leprevost, F.V., Avtonomov, D.M., Mellacheruvu, D., and Nesvizhskii, A.I. (2017). MSFragger: ultrafast and comprehensive peptide identification in mass spectrometry-based proteomics. Nat Methods 14, 513–520. doi:10.1038/nmeth.4256.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.4256&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=28394336&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
77. 77.Yu, F., Haynes, S.E., Teo, G.C., Avtonomov, D.M., Polasky, D.A., and Nesvizhskii, A.I. (2020). Fast Quantitative Analysis of timsTOF PASEF Data with MSFragger and IonQuant. Mol Cell Proteomics 19, 1575–1585. doi:10.1074/mcp.TIR120.002048.
[Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoibWNwcm90IjtzOjU6InJlc2lkIjtzOjk6IjE5LzkvMTU3NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAzLzA1LzIwMjMuMDUuMjQuMjMyOTAyOTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)
78. 78.da Veiga Leprevost, F., Haynes, S.E., Avtonomov, D.M., Chang, H.Y., Shanmugam, A.K., Mellacheruvu, D., Kong, A.T., and Nesvizhskii, A.I. (2020). Philosopher: a versatile toolkit for shotgun proteomics data analysis. Nat Methods 17, 869–870. doi:10.1038/s41592-020-0912-y.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41592-020-0912-y&link_type=DOI)
79. 79.Bern, M.W., and Kil, Y.J. (2011). Two-dimensional target decoy strategy for shotgun proteomics. J Proteome Res 10, 5296–5301. doi:10.1021/pr200780j.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/pr200780j&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=22010998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
80. 80.Demichev, V., Szyrwiel, L., Yu, F., Teo, G.C., Rosenberger, G., Niewienda, A., Ludwig, D., Decker, J., Kaspar-Schoenefeld, S., Lilley, K.S., et al. (2022). dia-PASEF data analysis using FragPipe and DIA-NN for deep proteomics of low sample amounts. Nat Commun 13, 3944. doi:10.1038/s41467-022-31492-0.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-022-31492-0&link_type=DOI)
81. 81.van de Bunt, M., Manning Fox, J.E., Dai, X., Barrett, A., Grey, C., Li, L., Bennett, A.J., Johnson, P.R., Rajotte, R.V., Gaulton, K.J., et al. (2015). Transcript Expression Data from Human Islets Links Regulatory Signals from Genome-Wide Association Studies for Type 2 Diabetes and Glycemic Traits to Their Downstream Effectors. PLoS Genet 11, e1005694. doi:10.1371/journal.pgen.1005694.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1005694&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=26624892&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
82. 82.Lamble, S., Batty, E., Attar, M., Buck, D., Bowden, R., Lunter, G., Crook, D., El-Fahmawi, B., and Piazza, P. (2013). Improved workflows for high throughput library preparation using the transposome-based Nextera system. BMC Biotechnol 13, 104. doi:10.1186/1472-6750-13-104.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1472-6750-13-104&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=24256843&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
83. 83.Yoon, J.S., Sasaki, S., Velghe, J., Lee, M.Y.Y., Winata, H., Nian, C., and Lynn, F.C. (2022). Calcium-dependent transcriptional changes in human pancreatic islet cells reveal functional diversity in islet cell subtypes. Diabetologia 65, 1519–1533. doi:10.1007/s00125-022-05718-1.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-022-05718-1&link_type=DOI)
84. 84.Stekhoven, D.J., and Buhlmann, P. (2012). MissForest--non-parametric missing value imputation for mixed-type data. Bioinformatics 28, 112–118. doi:10.1093/bioinformatics/btr597.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr597&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=22039212&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298383800017&link_type=ISI)
85. 85.Love, M.I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15, 550. doi:10.1186/s13059-014-0550-8.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-014-0550-8&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=25516281&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
86. 86. Jeffrey T. Leek, W. Evan Johnson, Hilary S. Parker, E.J.F., and Andrew E. Jaffe, Y.Z., John D. Storey, Leonardo Collado Torres (2022). sva: Surrogate Variable Analysis. R package version 3.46.0.
87. 87.Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9, 559. doi:10.1186/1471-2105-9-559.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-9-559&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=19114008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
88. 88.Kim, S. (2015). ppcor: An R Package for a Fast Calculation to Semi-partial Correlation Coefficients. Commun Stat Appl Methods 22, 665–674. doi:10.5351/CSAM.2015.22.6.665.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.5351/CSAM.2015.22.6.665&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=26688802&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
89. 89.Yu, G., Wang, L.G., Han, Y., and He, Q.Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287. doi:10.1089/omi.2011.0118.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/omi.2011.0118&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=22455463&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)
[Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000303653300007&link_type=ISI)
90. 90.van Gurp, L., Fodoulian, L., Oropeza, D., Furuyama, K., Bru-Tari, E., Vu, A.N., Kaddis, J.S., Rodríguez, I., Thorel, F., and Herrera, P.L. (2022). Generation of human islet cell type-specific identity genesets. Nat Commun 13, 2020. doi:10.1038/s41467-022-29588-8.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-022-29588-8&link_type=DOI)
91. 91.Segerstolpe, Å., Palasantza, A., Eliasson, P., Andersson, E.M., Andréasson, A.C., Sun, X., Picelli, S., Sabirsh, A., Clausen, M., Bjursell, M.K., et al. (2016). Single-Cell Transcriptome Profiling of Human Pancreatic Islets in Health and Type 2 Diabetes. Cell Metab 24, 593–607. doi:10.1016/j.cmet.2016.08.020.
[CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2016.08.020&link_type=DOI)
[PubMed](http://medrxiv.org/lookup/external-ref?access_num=27667667&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F05%2F2023.05.24.23290298.atom)