Glucose metabolism, β-cell Ca²⁺ flux, and insulin secretion are pulsatile, and have been found to oscillate in both humans and mice (8, 9, 18-20). Because they are interconnected, understanding the factors governing oscillation patterns can inform about the mechanisms that regulate insulin secretion (6, 11, 21, 22). To explore the influence of genetic background on Ca²⁺ oscillations, we measured Ca²⁺ in islets of the eight CC founder strains, that together harbor as much genetic diversity as humans: A/J, C57BL/6J (B6), 129S1/SvlmJ (129), NOD/ShiLtJ (NOD), NZO/HILtJ (NZO), CAST/EiJ (CAST), PWK/PhJ (PWK), and WSB/EiJ (WSB).

All mice were maintained on a Western-style diet (WD) high in fat and sucrose for 16 weeks, prior to isolating their islets for Ca²⁺ imaging with Fura Red, a Ca²⁺-sensitive fluorescent dye (Figure 1A). We measured Ca²⁺ dynamics in response to four conditions: 1) 8 mM glucose (8G); 2) 8 mM glucose + 2 mM glutamine, 0.5 mM leucine, and 1.5 mM alanine (8G/QLA); 3) 8G/QLA + 10 nM glucose-dependent insulinotropic polypeptide (8G/QLA/GIP); and 4) 2 mM glucose (2G) (Figure 1B). There was a high degree of similarity between three of the five classical strains (A/J, B6, 129), which were dominated by slow oscillations (period 2-10 minutes) in 8G and 8G/QLA/GIP, and had relatively fewer islets reach plateau (continuous peak activity without oscillation) in 8G/QLA. Likewise, the wild-derived strains (CAST, WSB, and PWK) closely matched one another, while differing from the classic strains. The wild-derived mouse islets were dominated by fast oscillations (period <2 minutes) in 8G, resulting in plateaus for 8G/QLA and 8G/QLA/GIP.

Two strains stood out from the others. Islets from NOD mice showed characteristics from both the wild-derived and classic strains; slow oscillations in 8G and a sustained plateau in response to 8G/AA and 8G/QLA/GIP with fast oscillations superimposed. The NZO mice also differed from the other classic strains, likely because they were all diabetic (blood glucose >250 mg/dL). Their islets were minimally responsive to 8G, but did respond with a strong pulse in 8G/QLA and Ca²⁺ remained elevated in 8G/QLA/GIP.

Many of the strain differences seen in the male mice were maintained in the females (Supplemental Figure 1A). The classic strains were once again highly similar to one another, as were the wild-derived strains. Furthermore, the NZO females, of which all but one were diabetic, mirrored the behavior of the male islets. One interesting observation that emerged from the female islets is that the NOD females displayed a greater variation in their Ca²⁺ oscillations than the NOD males (Supplemental Figure 2). Some of the islets maintained slow oscillations throughout the various conditions, while some demonstrated fast oscillations and plateaued like the wild-derived strains. Yet others appeared strikingly similar to the islets from diabetic NZO mice, despite none of the NOD mice being diabetic. Finally, the one non-diabetic female NZO displayed oscillatory behavior comparable to that of the other classic strains, with clear, slow oscillations (Supplemental Figure 1B).

An understanding of the mechanisms regulating insulin secretion, including the roles of specific metabolic pathways, ion channels, and hormones, has been derived from the shape and frequency of islet Ca²⁺ oscillations (6, 9, 11, 19, 23-26). To elucidate strain differences in Ca²⁺ dynamics, we focused on six parameters of the Ca²⁺ waveform (Figure 2A): 1) silent duration (the electrically-silent “triggering” phase, also known as “Mito_Cat” (8), which culminates in K_ATP closure and membrane depolarization), 2) peak Ca²⁺ (the top of each oscillation); 3) active duration (the length of time for each Ca²⁺ oscillation measured at half of the peak height, also known the oxidative “secretory” phase, or “Mito_Ox” (8); 4) pulse duration (active duration plus extra time for Ca²⁺ extrusion), 5) period (the length of time between two peaks, and the sum of silent duration and active duration); 6) plateau fraction (the active duration divided by the period, or the fraction of time spent in the active “secretory” phase). We also assessed the spectral density for every islet to extract additional information from complex oscillations where multiple components were visible (Supplemental Figure 3A). We analyzed each trace to determine the top two frequencies contributing to the trace (1^st and 2^nd component frequencies, Supplemental Figure 3B) and their respective contributions (1^st and 2^nd component amplitudes). Because certain features, such as metabolically-driven (slow) and electrically-driven (fast) oscillations have characteristic frequencies (27), extracting the top two frequencies may highlight additional information beyond that previously collected.

A representative Ca²⁺ dynamic from a female B6 islet is illustrated in Figure 2A. The transition from 8G to 8G/QLA resulted in an increased active duration, yielding a longer period and increased plateau fraction. For an islet that plateaued at the peak, as seen in 8G/QLA (Figure 2B), we computed a plateau-fraction of one, an active and pulse duration of 40 minutes (the measurement time), and a period of zero minutes. An islet that returned to baseline and ceased to oscillate, as seen in 2 mM glucose (Figure 2B), was determined to have a plateau-fraction, active duration, and pulse duration of zero, and a period of 40 minutes. Dissecting the strain-dependence of these key parameters of the Ca²⁺ oscillations is important for identifying underlying mechanisms, as illustrated in Figure 2C. While both traces have a similar active duration (blue bars), trace 1 has a longer period (red bars), resulting in an increased silent duration and a decreased plateau-fraction.

Examples of pathways altering specific components of Ca²⁺ oscillations have previously been established (Figure 2D) (9, 11, 18, 28, 29). For example, when K_ATP channels are pharmacologically closed with tolbutamide, the silent duration is shortened, resulting in increased frequency without a change in pulse shape (upper panel). The addition of glucose leads to increased glucose metabolism and glucokinase (GK) activity (11). The resulting rise in ATP inhibits K_ATP channels (6) and is used as a substrate for additional processes that affect Ca²⁺, such as SERCA pumps (30, 31). Thus, glucose alters both the active and silent durations, resulting in a change in both frequency and shape of the Ca²⁺ oscillations (lower panel).

Average Ca²⁺ is commonly used for analyzing Ca²⁺ dynamics and is frequently assumed to be highly correlated to insulin secretion. To determine whether average Ca²⁺ is predictive of insulin secretion, we performed ex vivo perifusion studies on islets from WSB and 129 male mice, two strains that showed similar average Ca²⁺ (Figure 4B) but exhibited vastly different Ca²⁺ oscillations (Figure 1B). WSB mice had significantly higher insulin secretion in each of the secretory conditions (Figure 3A), suggesting another Ca²⁺ parameter better predicts insulin secretion.

To identify parameters of the Ca²⁺ dynamics most strongly correlated to insulin secretion, we computed the correlation between the Ca²⁺ oscillation parameters and insulin secretion in similar conditions (8G, 8G/QLA, basal) for the same sexes and strains (Figure 4A and Supplemental Figure 4) (32). Consistent with our observations from the perifusion data in the WSB and 129 islets, we found that average Ca²⁺ was not strongly correlated to insulin secretion. Other metrics, such as active duration in 8G, and the silent durations in 8G/QLA, were more highly correlated to insulin secretion. Meanwhile, the 1^st component frequency in 8G from the spectral density analysis was highly correlated with decreased insulin secretion. These metrics were also the most highly correlated with multiple clinical measures in the founder mice, particularly plasma insulin (Figure 4B, Supplemental Figure 5), for which silent duration in 8G/QLA/GIP had the strongest correlations.

Several parameters of the Ca²⁺ oscillatory waveform showed strong strain and sex effects. For example, basal Ca²⁺ (average Ca²⁺ in 2G, Figure 4C) was relatively consistent among the strains, except NZO where it was highest in islets from male mice. For the overall pulse duration (Figure 4D), the NZO mice were once again the highest, followed by CAST and WSB. A noticeable sex-effect was measured for the CAST mice, where male mice had a longer pulse duration than the female mice. The 1^st component frequency (Figure 4E) is driven by the differences observed in the wild-derived strains, for which CAST has the highest frequency, followed by PWK and WSB. Finally, the trend for a sex effect in the classic strains on the silent duration (at 8G, 8G/QLA, and 8G/QLA/GIP) is absent in the NZO and wild-derived mice with the former having greater silent duration in males and the latter frequently having islets plateau in response to these stimuli. These data demonstrate that genetics has a profound influence on key parameters of islet Ca²⁺ oscillations.

To explore relationships between Ca²⁺ oscillations and islet proteins, we took advantage of our whole islet proteomic survey from the eight founder strains (17). To identify proteins that may underly the strain-differences in Ca²⁺ oscillations, we computed the correlation between islet protein abundance and Ca²⁺ dynamics across all mice used in our study (Figure 5 and Supplemental Figure 6). Our previous survey of islet proteomics included both sexes for all strains, except NZO males, resulting in a quantitative measure of 4054 proteins (17). Figure 5A illustrates a heatmap of the correlation between islet proteins and several parameters of Ca²⁺ oscillations. Unsupervised clustering was used to show that groups of proteins showed strong positive or negative correlation to a given Ca²⁺ parameter, yielding distinct correlation architecture. For example, proteins highly correlated to the 8G 1^st component frequency, tended to also be strongly anticorrelated to the silent duration conditions, which were very similar to one another. The active and pulse durations for 8G had nearly identical correlation structure. Additionally, the conditions with the fewest highly correlated proteins were the average Ca²⁺ measures for 8G, 8G/QLA, and 8G/QLA/GIP, and the structure for these was largely inverted from the active duration conditions. Finally, despite the differences in the overall correlations between the different metrics, there were proteins that did overlap (e.g. the block of proteins with high correlation to both 8G A_D and 8G/QLA S_D) suggesting that while there were clusters of distinct proteins/pathways for any given metric some proteins may modify more than one metric.

Among the 4054 islet proteins, 363 had high absolute correlation coefficients (r > |0.5|) to 3 or more of the parameters our data suggest most strongly correlate to insulin secretion and plasma insulin (Basal Ca²⁺, 8G A_D, 8G P_D, 8G/QLA S_D, 8G S_D, 8G/QLA/GIP S_D). Interestingly, of the proteins correlated to these traits, many have been previously implicated in islet biology, including PCSK1, GCK, SUR1, GLUT2, PDX1, and GLP-1 (28-36). Notably, the highly correlated proteins enriched for tissues, pathways, and transcription factors that support their role in insulin secretion (Figure 5B, C, D, Enrichr links in Supplemental Table 1). For instance, proteins highly anti-correlated to active duration in 8G were enriched for components of oxidative metabolism and had their gene promoters enrich for binding to the islet transcription factor MAFA (33). These enrichment data provide a framework for discovering new genes of interest for their role in islet function.

The data presented in Figure 5A illustrates the correlation between islet proteins and Ca²⁺ dynamics. Importantly, a protein strongly correlated to Ca²⁺ does not necessarily reflect a causal relationship, i.e., a change in protein abundance may or may not cause a change in the Ca²⁺ signal. To take our analysis beyond correlation, we integrated our data with human GWAS of glycemia-related traits.

For each Ca²⁺ parameter, we focused on those proteins in the tails of the correlation histogram where r > |0.5| (e.g. Figure 5B). We identified human homologues for 3073 proteins that were correlated to Ca²⁺ in either direction for at least one of our parameters of interest. We then searched the Type 2 Diabetes Knowledge Portal (https://t2d.hugeamp.org/) for SNPs that are associated with one or more glycemia-related traits (see Table 2) with a P-value < 10^-8, and are located within ±100 kbp of the homologous gene (e.g. COBLL1, Figure 6A), or in regions that contact the gene’s promoter region (determined using human islet promoter-capture HiC data (37)) as illustrated by ACP1 (Figure 6B). This yielded a list of 647 human genes strongly associated with diabetes-related SNPs. Among these genes, 478 were not previously associated with insulin secretion (see Methods), suggesting they may have understudied roles in islet function (Figure 7A, B; Table 3, Table 4). In summary, our approach leverages the genetic diversity of the eight CC founder strains and human GWAS for diabetes-related traits to highlight genes that may play a novel role in islet function and relate to diabetes risk.

Importantly, the relevance of these tools depends on the mouse model in which they are used to validate the roles of the candidate protein. To aid in the selection of mouse strain, we provide a resource with proteomic and Ca²⁺ data (Figure 7C, D, E; https://data-viz.it.wisc.edu/FounderCalciumStudy/, https://doi.org/10.5061/dryad.j0zpc86jc, https://rstudio.it.wisc.edu/FounderCalciumStudy). This will enable the user to identify proteins correlated to proteins or traits of interest, and from there, identify which strain(s) may be most appropriate for studies of their protein of interest. In the examples illustrated in Figure 7C, the user queries for the strain/sex distribution of the protein GALNS, which shows a high negative correlation to multiple traits including the active duration time in 8G/QLA. Strains at the extremes of this trait are also extremes regarding GALNS protein abundance. Strains with high abundance (AJ, for example (yellow arrow)) would be ideal models for inhibition or knockout, whereas CAST mice (green arrow) could be a comparison strain for validating the role of the protein, as they express much less GALNS. Users can also query for the proteins or calcium parameters and see their correlation to the other proteins and/or calcium parameters (Figure 7D). Finally, the user can see which traits or proteins show the strongest strain, sex, or sex-by-strain effects using the options in the volcano plot (Figure 7E). Together, these and other tools in the resource will allow researchers to explore their traits or proteins of interest as well as determine the appropriate model systems and conditions that may best interrogate their experimental questions of interest.

While the development and progression of T2D is potentiated by environmental factors, an estimated 50% of disease risk is driven by genetic factors (1, 3, 38). Therefore, to study the genetic variation contributing to T2D, we took advantage of the genetic diversity contained within the eight CC founder strains. These mice collectively contain a level of genetic diversity mirroring that seen in humans, making them an excellent experimental platform to link genetics with altered islet function (14, 15). We demonstrate that they also vary in their Ca²⁺ response to various insulin secretagogues, supporting the use of these mice to identify novel genes involved in regulating islet biology.

Variations in Ca²⁺ dynamics are highly complex, reflecting changes in metabolism, extra-islet signaling, and Ca²⁺ itself (8). We therefore selected stimulatory conditions to assess each of these components in islets of the eight mouse strains. 8 mM glucose was first used to survey glycolytic responses, because we have observed that several strains reliably oscillate at this glucose concentration. Furthermore, this glucose concentration remains close to the stimulatory threshold, thus reducing the possibility of oscillations plateauing if islet Ca²⁺ responses were left-shifted in any strains (23, 39, 40). We then added QLA as fuel to engage mitochondrial metabolism and paracrine signaling from α-cells, providing a survey of α-to-β cell communication in the islet (6, 7, 41, 42). Finally, we used GIP to interrogate the islet incretin responses and the cAMP amplification pathway (41), before returning to a low glucose concentration that enabled us to establish baseline Ca²⁺ levels.

The variation in Ca²⁺ response to these conditions can be better understood by examining the multitude of pathways regulating Ca²⁺ dynamics. As mentioned previously, altering ionic pathways involved in regulating Ca²⁺, such as K_ATP channels, has a very different effect on Ca²⁺ oscillations compared to altering glycolysis. It is important to further dissect these pathways, as changing specific components of glucose metabolism can elicit different effects. For example, activating glucokinase (GK, which is rate limiting for glycolysis) and activating pyruvate kinase (PK, an ATP-generating enzyme which directly controls K_ATP channel closure in the final step of glycolysis) both reduce the silent phase duration by accelerating K_ATP channel closure. GK and PK activation also alter the active duration and the plateau fraction, however they do so in opposite directions (6). Activating GK increases the active duration and oscillation period, while activating PK decreases those same parameters via Ca²⁺ extrusion (7). The incretin hormones GLP-1 and GIP reduce silent duration and oscillation period, most likely due to their ability to activate Epac2 and sensitize K_ATP channels to ATP-dependent closure (43, 44). Thus, PK and GLP1 have a common target (i.e. K_ATP), and therefore a similar effect on silent duration. These examples illustrate the benefit of analyzing multiple Ca²⁺ parameters for understanding pathways of interest.

The importance of analyzing a variety of Ca²⁺ parameters is further supported by the insulin secretion measurements in the male WSB and 129 mice. While average Ca²⁺ is a common metric used to predict insulin secretion, relying on only this metric would suggest that the two strains secrete insulin similarly. However, the WSB mice secreted significantly more insulin in 8G, 8G/QLA, and 8G/QLA/GIP (Figure 3). Based on our correlation analysis between Ca²⁺ parameters and insulin secretion across each sex and strain, active duration, and pulse duration in 8G more accurately predicted insulin secretion and may be highly informative when used with other data (Figure 4). This is similar to results published by some other groups, suggesting that average Ca²⁺ does not correlate well with insulin secretion (45).

Notably, several of the strains appeared to cluster together with similar responses. One such group is composed of three classical strains (A/J, B6, and 129), which had similar waveforms that were dominated by slow oscillations. These differed from a second group containing the wild-derived strains (CAST, WSB, and PWK) which closely matched one another and were dominated by faster oscillations. The classic strains have been highly inbred (>150+ generations) and descend from related common ancestors, the “fancy mice.” They also have extremely low genetic diversity, with 97% of their genomes explained by fewer than 10 haplotypes (3, 46, 47). By contrast, wild-derived strains are each independent in their parental origin, inbred for far fewer generations than the classic strains (20 vs. >150), and include significant contributions from other subspecies of Mus musculus than the predominant subspecies (M. m. domesticus) in the traditional strains, particularly CAST (M. m. castaneous) and PWK (M. m. musculus) (3, 46, 47). It is thus unsurprising that the two primary Ca²⁺ response clusters were composed of the classic and wild-derived strains. Multiple loci have already been linked to islet dysfunction and differential metabolic homeostasis in the classic strains (3). Our work here highlights the promise in using wild-derived strains to unmask previously underappreciated islet phenomena, something we and others have previously shown (17, 48, 49).

While considered one of the classic strains, the NOD mice differed from the two primary clusters noted above. They displayed a combination of features from both groups and had a high degree of inter-islet variability, especially the female mice. NOD share common ancestors with the Swiss-Webster mice, which do not share parental origin with the other classic strains (46) and also display a “mixed” phenotype consisting of islet Ca²⁺ oscillations in response to glucose, with both slow and fast components (23).

Additionally, while all NOD were normoglycemic, a heterogeneous response was observed in islets from female, but not male, NOD mice. Female NOD mice are known to develop islet immune infiltration and subsequent autoimmune diabetes whereas males are largely protected from this (50). Male NOD islets were largely consistent in their Ca²⁺ waveforms. In the females, however, a high degree of heterogeneity in responses was observed across the female’s islets (Supplemental Figure 2). For any NOD female, some islets resembled those from NOD males in their clear oscillations, others largely lacked oscillatory behavior other than a strong pulse in response to 8G/QLA, and still others had an intermediate response. These observations may reflect varying degrees of dysfunction in the NOD female islets as the mice progress to diabetes, though we cannot say whether this results from variation in beta-cell intrinsic defects or islet immune infiltration.

The NZO mice also varied from the two clusters previously discussed. Male, and all but one of the female NZO mice were diabetic. Islets from diabetic mice had reduced amplitude and oscillatory behavior, other than a single pulse in 8G/QLA. This pattern is similar to the patterns observed in many of the NOD female islets. On the other hand, islets of the one non-diabetic female NZO mouse demonstrated clear, slow oscillations (Supplemental Figure 1B), which was surprising given reports of low K_ATP abundance due to Abcc8 mutations in the NZO (51) and the strong role of K_ATP channels in regulating islet Ca²⁺(11, 32, 33, 52). While not of the same lineage as the NOD, the NZO do exhibit some autoimmune infiltration in the pancreas, and the marked difference between the non-diabetic and diabetic NZOs, along with the variation in female NOD islet responses, further suggests that intra-islet variability for the NOD mice may be the result of disease progression.

Understanding the genetic variation driving islet responses in the founders may be informative beyond these specific strains. Screens in the DO mice and similar outbred populations can track SNPs associated with trait variation to their parental inbred strain of origin. Our previous genetic screen for drivers of islet function observed that many of the quantitative trait loci (QTL) appearing for ex vivo islet traits had effects driven by the SNPs from the wild-derived strains as opposed to the traditional inbred groups (16). The QTL mediated by Zfp148, which drives Ca²⁺ and insulin secretion phenotypes in β-cells (40), also had strong strain effects from wild-derived strains (16).

Previous studies of islet Ca²⁺ have largely been confined to a handful of strains, and many studies by individual labs tend to use the strains linked with their projects. While this does include a few outbred stocks (e.g. NMRI (53-55), CD-1 (9, 56, 57)), direct comparisons of these to the traditional inbred lines are rare (53), and studies of specific genes often use traditional inbred lines as wild-derived lines do not respond as well to the conventional assistive reproductive technologies required for genome editing and transgenesis (58, 59).

One study, comparing the outbred (NMRI) stock to the C57BL6/J and C57BL/6N strains (53), found that the NMRI displayed significantly lower frequencies than the C57 lines, particularly in physiological glucose ranges and had similar active periods. While highly informative, there were important differences between these studies and our studies here. Of note, the studies were done in acute slice culture, in only one sex, and the frequencies detected did not resemble the (at least for the C57 lines) the slow oscillations observed in isolated islets from these inbred strains (e.g. Figure 1, (6, 40)).

One limitation of our current study is that the association between islet proteins and Ca²⁺ waveforms is correlative and therefore cannot distinguish proteins that are causal for the differences in islet Ca²⁺ between strains from proteins that change because of these differences. One approach to discriminate cause from effect, and establish the relevance to humans, is to identify whether genes encoding human orthologues of these proteins are associated with glycemic traits in human GWAS. SNPs for glycemic traits (Table 2), particularly those involving insulin, suggest that alterations in these proteins may impart disease risk, which is less likely for proteins that do not play critical regulatory roles. Thus, filtering our candidates for glycemic trait associations in human GWAS, while not definitive, suggests a likely causal role for these proteins in mediating differences in islet Ca²⁺ and insulin secretion among the different mice. Integrating human GWAS data with the proteins most correlated to Ca²⁺ dynamics nominated ∼650 protein candidates, of which approximately a third have been previously shown to have roles in islet biology. These include well established drivers of insulin secretion; e.g. SUR1, GLUT2, and GNAS. Other previously unknown candidates show promise for validation, as they are already targets of small molecule compounds (e.g. ACP1 and others, (60-65)), are secreted (e.g. COBLL1 and others, (66-73)), or have been knocked out in mice, resulting in metabolic phenotypes (Figure 7B, Table 3, Table 4 (74)).

In addition to the candidate regulators with potential relevance to human islet biology, we provide a user-friendly interface to our data where users can determine whether their gene of interest has a potential regulatory role in islets. Multiple inferences regarding the roles of specific pathways are possible via analysis of Ca²⁺ oscillations in islets (6, 8, 9, 19), and our protein correlation data provides a resource to identify which parameter most closely correlates to a number of Ca²⁺ traits. Additionally, it highlights strain/sex outliers for a given trait or gene product, which can be used to select which strain/sex is best to explore that gene’s role (e.g. Figure 7C). Newer technologies in reproductive assistance, transgenesis, and gene editing, together with more accurate genome sequencing and single mutations conferring docility, are quickly making utilization of the wild-derived mice more practical (58, 59, 75-77). As many of the QTL identified in DO-based studies often have strong driver SNPs from the wild-derived strains, a further understanding of which experimental questions might be best addressed by use of these strains will be important.

We have previously provided user-friendly web interfaces that allow searches of gene expression as a function of diet (WD vs chow, (78, 79)) and background (e.g. BTBR and B6, (78, 79)), correlation and QTL scans in F2 intercrosses of these mice (80), and where these may align with QTL in our DO studies (16). Many of our candidates are strongly altered by diet and have strong correlations in the F2 data for certain clinical traits including insulin and glycemic parameters (80). Here, we provide the correlation data for islet proteins against multiple parameters describing islet Ca²⁺ responses between strains (DOI: https://doi.org/10.5061/dryad.j0zpc86jc, https://data-viz.it.wisc.edu/FounderCalciumStudy/, https://github.com/byandell/FounderCalciumStudy, https://rstudio.it.wisc.edu/FounderCalciumStudy). These will enable researchers to better identify proteins or parameters of interest as well as appropriate background strains with which to determine the functions of these proteins.

All general chemicals, amino acids, BSA, DMSO, glucose, glucose-dependent insulinotropic polypeptide (GIP, G2269), cOmplete Mini EDTA-free Protease Inhibitor Cocktail Tablets (11836170001), and heat-inactivated FBS (12306C) were purchased from Sigma Aldrich. RPMI 1640 base medium (11-875-093), antibiotic–antimycotic solutions (15240112), NP-40 Alternative (492016), Fura Red Ca²⁺ imaging dye (F3020), DiR (D12731), and agarose (BP1356-500) were purchased from ThermoFisher. Glass-bottomed culture dishes were ordered from Mattek (P35G-0-14-C). Fura Red stocks were prepared at 5 mM concentrations in DMSO, aliquoted into light-shielded tubes, and stored at –20°C until day of use (5 μM final concentration). DiR was prepared in DMSO at 2 mg/ mL, aliquoted to light-shielded tubes, and stored at 4°C until use. All imaging solutions were prepared in a bicarbonate/HEPES-buffered imaging medium (formula in Table 1). Amino acids were prepared as 100× stock in the biocarbonate/HEPES-buffered imaging medium, aliquoted into 1.5 mL tubes, and frozen at –20°C until day of use. Aliquots of GIP stock were prepared at 100 μM in water and kept at - 20°C until day of use.

Animal care and experimental protocols were approved by the University of Wisconsin-Madison Animal Care and Use Committee. Most strains (B6, AJ, 129, NOD, PWK, and WSB) were bred inhouse, although two strains (CAST and NZO) were purchased from Jackson Laboratory (Bar Harbor, ME). All mice were fed a high-fat, high-sucrose Western-style diet (WD, consisting of 44.6% kcal fat, 34% carbohydrate, and 17.3% protein) from Envigo Teklad (TD.08811) beginning at 4 weeks and continuing until sacrifice (aged ∼19-20 weeks for all strains except the NZO males). The NZO males were sacrificed at 12 weeks of age owing to complications from severe diabetes. For each strain, 3-7 males and females from at least 2 litters were analyzed. Animals were sacrificed by cervical dislocation prior to islet isolation.

Fasting blood glucose and insulin levels were measured in mice at 19 weeks of age, except for the NZO males which were measured at 12 weeks of age. Glucose was analyzed by the glucose oxidase method using a commercially available kit (TR15221, Thermo Fisher Scientific), and insulin was measured by radioimmunoassay (RIA; SRI-13K, Millipore).

Islets were isolated as previously described (81) and incubated in recovery medium (RPMI 1640, 11.1 mM glucose, 1% antibiotic/antimycotic, 10% FBS) overnight at 37°C and 5% CO₂. Islets were then incubated with Fura Red (5 μM in recovery medium) at 37°C for 45 minutes. Imaging dishes were created from glass-bottomed 10 cm² dishes that had been filled with agarose. A channel with a central well was cut into the agarose with expanded ports on either side of the well for inflow and outflow lines. Prior to loading the chambers were perfused with the initial imaging solution (8 mM glucose in imaging medium). Islets were then loaded into these dishes. The imaging chamber was placed on a 37°C– heated microscope stage (Tokai Hit TIZ) of a Nikon A1R-Si+ confocal microscope. All solution reservoirs were kept in a 37°C water bath. Solutions were perfused through the chamber at 0.25 mL/min, with constant flow controlled by a Fluigent MCFS-EZ and M-switch valve assembly (Fluigent). The scope was integrated with a Nikon Eclipse-Ti Inverted scope and equipped with a Nikon CFI Apochromat Lambda D 10x/0.45 objective (Nikon Instruments), fluorescence spectral detector, and multiple laser lines (Nikon LU-NV laser unit; 405, 440, 488, 514, 561, 640nm). Bound dye was excited with the 405nm laser and the spectral detector’s variable filter was set to 620-690nm. The free dye was excited with the 488nm laser and the variable filter collected from 640-690nm. Images were collected at 1 frame/sec at 6 second intervals. Each islet was considered a region of interest for further analysis. ROI intensity was collected by NIS Elements and exported for further analysis. All microscopy was performed at the University of Wisconsin-Madison Biochemistry Optical Core, which was established with support from the University of Wisconsin-Madison Department of Biochemistry Endowment.

Isolated islets were kept in RPMI-based medium (see above) overnight prior to perifusion, which was performed as previously described, with minor modifications (40, 82). Islets were equilibrated in 2 mM glucose for 55 minutes, after which 100 μL fractions were collected every minute with the perifusion solutions set at a flow rate of 100 μL/min. All solutions and islet chambers were kept at 37°C. After the final fraction was collected, islet chambers were disconnected, inverted, and flushed with 2 mL of NP-40 Alternative lysis buffer containing protease inhibitors for islet insulin extraction.

Insulin in each perifusion fraction and islet insulin content were determined using a custom assay, as previously described (17).

Trace segments for each solution condition were analyzed using Matlab and R. Traces were detrended using custom R scripts and Graphpad PRISM. Custom Matlab scripts (7) (https://github.com/hrfoster/Merrins-Lab-Matlab-Scripts, also stored on Dryad https://doi.org/10.5281/zenodo.6540721) determined oscillation peak amplitude, pulse duration, active duration (the time when Ca²⁺ is above 50% peak amplitude), silent duration (the difference between period and active duration), plateau fraction (the fraction of overall time per pulse spent in the active duration), pulse period and other parameters. Spectral density deconvolution for the trace segments to determine principal frequencies was done using R. Animal averages for the different parameters defined by Matlab and R were computed and graphed using custom R scripts. Figures were created using CorelDraw and Biorender.com. All R scripts and the citations for the relevant packages used to generate them are available via Dryad (https://doi.org/10.5061/dryad.j0zpc86jc).

Correlation analysis was performed using the imaging data measurements and our published islet protein abundance data, ex vivo static insulin secretion measurements, and in vivo measurements made in a separate cohort of mice on the WD from the same strains and sexes used in these studies (17). For each imaging parameter or previously published measurement, the Z-score was calculated using the formula z = (x-μ) / σ where z is the Z-score, x is the animal average for that trait given the strain and sex, μ is the average of all animals’ values for that trait, and σ is the standard deviation for all animals’ values for that trait. Z-scores were computed in R and excel for the imaging parameters and the previously published (17) islet proteomic, ex vivo secretion, and in vivo measurements.

Correlation coefficients between the Z-score values of the imaging parameters and Z-scores of the previously published protein abundance, islet secretion, and in vivo traits were computed in Excel using the CORREL function. The equation used for this function is:

Where X and Y are the Z-scores for the correlated traits/parameters, ẋ is the population average for trait X and ẏ is the population average for trait Y. Traits were considered highly correlated if absolute value for their Z-score correlation coefficients was ≥ 0.5.

Proteins highly correlated or anticorrelated to imaging parameters were further analyzed using pathway enrichment and presence of human GWAS SNPs. Briefly, for a given parameter, pathway analysis for the highly correlated or anti-correlated proteins to that parameter was done using Enrichr (83, 84). Enrichr links for the subsets of proteins highly correlated to specific calcium parameters are provided in Supplemental Table 1, which is stored on Dryad (https://doi.org/10.5061/dryad.j0zpc86jc).

For GWAS analysis, human orthologues for genes encoding the previously measured islet proteins were identified using BioMart (85). For highly correlated proteins, the protein was deemed of human interest if its orthologue had SNPs for glycemia-related traits (see table 2) either along the gene body, within +/- 100 kbp of the gene start or end, or if any region in the gene body was connected to regions with SNPs by chromatin looping. SNPs were queried using Lunaris tool of the Common Metabolic Diseases Knowledge Portal (cmdkp.org). Chromatin loop anchor points for the relevant gene orthologues were identified using previously published human islet promoter-capture HiC data (37) and the alignment between these anchor loops and orthologues of interest was done using R scripts.

For those proteins having orthologues with SNPs via this analysis, we conducted further literature searches using Pubmed, Google Scholar, ChEMBL (62, 63, 86), canSAR (61), Uniprot (66), Tabula Muris (87), and the Human Protein Atlas (67, 68), and other resources (69, 88, 89) to determine tissue expression and identify any prior roles in islet biology. Figures for the relevant protein examples were created using GraphPad Prism, CorelDraw, and the WashU Epigenome Browser (90).

We further narrowed this list by searching each of the genes and their aliases in the PubMed, Google Scholar, and Google Search Engines along with “insulin secretion.” This allowed us to identify which genes have a known role in altering the insulin secretory pathway, and which genes may be understudied (Figure 7 A).

A web resource was created to explore the islet calcium and proteomic data and their relationships (https://data-viz.it.wisc.edu/FounderCalciumStudy, https://rstudio.it.wisc.edu/FounderCalciumStudy). This resource sits on an RStudio/Connect server (see https://posit.co/). It enables the user to select traits from the calcium and protein datasets to plot by strain, sex, and calcium parameters. Distinct mice were assayed for calcium and protein. Individual strains can be selected on the main menu using the checkboxes, or all strains (default) can be viewed.

The different datasets available in the main menu are:

The calcium data have three stimulatory conditions (8G, 8G/QLA, and 8G/QLA/GIP) that are displayed together for each calcium parameter. The proteomic data (protein) are displayed for each identified peptide. In rare cases of multiple peptides per gene, both gene symbol and peptide identifier (PP number) are included (e.g. Pkm_PP_1521 for the M1 isoform of the protein PKM). Desired proteins can be selected simultaneously with desired calcium parameters for correlation analysis and paired display by both datasets. The basal elements retained from the calcium data include the Average Calcium measurement for 2mM glucose. Proteomic data were log10-transformed. All traits were transformed into normal scores, keeping the sample mean and variance the same.

Scatter plots display data across sex and calcium conditions. When plotting calcium against protein or basal traits, means by strain and sex are used, as the two experiments used different mice. Correlation of selected traits with all other traits in the resource use Pearson correlation on pairwise-complete data. The user can order traits by their significance or by their correlation to other selected traits.

Statistical modeling terms include strain, sex, and the strain:sex interaction, plus additional terms for comparing calcium condition with respect to strain and sex. Users can view volcano plots displaying deviation of term effects, measured as the standard deviation (SD = square root of mean square error) divided by the raw SD for that trait, against their significance (p-value after adjusting for all other model terms, presented on -log10 scale). In addition to the terms, a composite “signal” captures the combined effect of terms strain:condition + strain:condition:sex using a general F test computation.

All data handling and web app construction for the resource was performed using R scripts in publicly available GitHub repositories, with specifics for the calcium study at https://github.com/byandell/FounderCalciumStudy and the general purpose analysis and web deployment package at https://github.com/byandell/foundr.

For the islet perifusion insulin measurements, statistics were determined in GraphPad Prism. Fractional secretion area-under-the-curve (AUC) was determined using Prism and differences in AUCs analyzed using post-tests following 2-way ANOVA for the indicated trace segments. Islet total insulins between strains were compared using a two-tailed Student’s t-test with Welch’s correction.

All R scripts and raw data are available via Dryad (https://doi.org/10.5061/dryad.j0zpc86jc). Code for the web resource is available via Github (https://github.com/byandell/FounderCalciumStudy). Image files are available upon request.

All protocols were approved by the University of Wisconsin-Madison IACUC (Protocol A005821-R01)