Maša Skelin Klemen1
Jurij Dolenšek1,2
Lidija Križančić Bombek1
Viljem Pohorec1
Marko Gosak1,2,3
Marjan Slak Rupnik1,3,4
Andraž Stožer1*- 1Institute of Physiology, Faculty of Medicine, University of Maribor, Maribor, Slovenia
- 2Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia
- 3Alma Mater Europaea, European Center Maribor, Maribor, Slovenia
- 4Center for Physiology and Pharmacology, Medical University of Vienna, Vienna, Austria
Beta cells couple stimulation by glucose with insulin secretion and impairments in this coupling play a central role in diabetes mellitus. Cyclic adenosine monophosphate (cAMP) amplifies stimulus-secretion coupling via protein kinase A and guanine nucleotide exchange protein 2 (Epac2A). With the present research, we aimed to clarify the influence of cAMP-elevating diterpene forskolin on cytoplasmic calcium dynamics and intercellular network activity, which are two of the crucial elements of normal beta cell stimulus-secretion coupling, and the role of Epac2A under normal and stimulated conditions. To this end, we performed functional multicellular calcium imaging of beta cells in mouse pancreas tissue slices after stimulation with glucose and forskolin in wild-type and Epac2A knock-out mice. Forskolin evoked calcium signals in otherwise substimulatory glucose and beta cells from Epac2A knock-out mice displayed a faster activation. During the plateau phase, beta cells from Epac2A knock-out mice displayed a slightly higher active time in response to glucose compared with wild-type littermates, and stimulation with forskolin increased the active time via an increase in oscillation frequency and a decrease in oscillation duration in both Epac2A knock-out and wild-type mice. Functional network properties during stimulation with glucose did not differ in Epac2A knock-out mice, but the presence of Epac2A was crucial for the protective effect of stimulation with forskolin in preventing a decline in beta cell functional connectivity with time. Finally, stimulation with forskolin prolonged beta cell activity during deactivation, especially in Epac2A knock-out mice.
Introduction
Insulin serves as an important regulator of nutrient homeostasis and a lack of its effects, either due to beta cell failure or resistance in peripheral tissues, or both, results in diabetes mellitus, which affects more than 537 million people worldwide (1). 90% of all diabetes cases are type 2 diabetes mellitus (T2DM), which is characterized by impaired insulin secretion, insulin resistance, and obesity (2).
Pancreatic beta cells act as sensors for blood glucose changes and respond with altering insulin secretion via a process termed stimulus-secretion coupling, by activation of the triggering and amplifying pathways (3, 4). The triggering pathway is initiated by glucose entry into beta cells through glucose transporters (GLUT), followed by glucose metabolism and intracellular ATP concentration increase, which results in decreased open probability of the KATP channels. The subsequent decrease in potassium efflux depolarizes the cell membrane, leading to an increase in permeability of VDCC allowing Ca2+ ions influx into the cytoplasm (5). Being the central intracellular messenger, Ca2+ activates exocytosis of insulin-containing granules (3, 6, 7).
The triggering pathway is augmented and fine-tuned by amplifying pathways which start with either (i) an increase in [Ca2+]IC (8), (ii) GPCRs activation (9), or (iii) phospholipase C (PLC) activation (10). The result of either [Ca2+]IC elevation or GPCR activation is enhanced adenylyl cyclase activity, which increases the concentration of cAMP. Consequently, cAMP enhances fusion of insulin-containing vesicles with plasma membrane via both a protein kinase A (PKA)-dependent pathway and activation of the guanine nucleotide exchange protein 2 (Epac2A) (11–13), whose role is not fully elucidated yet.
When blood glucose levels are high, PKA-dependent amplification mostly increases the second phase of insulin secretion (14) by inhibiting KATP channel activity in an ADP-dependent manner (15), as well as inhibits K+ currents through voltage dependent KV1.4 channels (16, 17). Further enhancement of Ca2+ signaling can be mediated through VDCCs (18–20) or Ca2+ mobilization from internal stores through both IP3 and ryanodine receptors (18–22). In addition, PKA increases the number of insulin-containing granules that are highly sensitive to Ca2+ (23, 24), the mobility and replenishment of the readily releasable pool (RRP) of insulin-containing granules (11), and the overall sensitivity of the secretory machinery to [Ca2+]IC (25).
Along with influences on PKA, cAMP affects insulin secretion via Epac2A (13, 26–28), which has been implicated in both the first and the second phase of insulin secretion (13, 14). While the Epac1 isoform is ubiquitous, Epac2A is located mainly in the neural, endocrine, and neuroendocrine tissue (27, 29), and seemingly predominant in mouse islets (30). The effects of Epac2A on insulin secretion are mediated in a myriad of ways. It has been proposed that GLP-1-induced mitochondrial ATP synthesis, which contributes to closure of KATP channels, is partially mediated by Epac2A (31). Furthermore, Epac2A contributes to increased activity of the enzyme glucokinase, increasing the production of mitochondrial metabolism substrates (32). Under the influence of Epac2A, KATP channels exhibit increased ATP sensitivity, facilitating membrane depolarization (33, 34). Epac2A also influences [Ca2+]IC dynamics by RYR2 sensitization, thus encouraging Ca2+ release from the endoplasmic reticulum Ca2+ stores upon stimulation with increased [Ca2+]IC via the Rap protein and PLC-ϵ (35), a process designated as calcium-induced calcium release (36, 37). Epac2A may also augment Ca2+ influx by influencing VDCCs via Rim2α (38, 39). Downstream from [Ca2+]IC, Epac2A increases the size of the readily releasable pool of granules through Epac2/Rap1 signaling. It also influences docking and priming of granules when complexed with the Rab3-interacting molecule Rim2α (39), and stimulates insulin granule acidification, an important step in granule priming by regulating Cl- influx (12). Previous research on mice lacking Epac2A (Epac2A KO) has demonstrated a propensity for obesity (40), a diminished first phase of insulin secretion (13), and impaired glucose-stimulated insulin secretion under diet-induced insulin resistance (41).
Besides exhibiting complex and intertwined intracellular signaling pathways, beta cells are intrinsically highly heterogeneous and operate in a multicellular environment. They communicate with each other through different mechanisms, whereby the electrical coupling through gap-junctions is recognized as the main synchronizing agent that facilitates the spreading of depolarization and Ca2+ waves across the islets (42–48). This enables beta cells to operate in a coordinated manner, which was found essential for the normal control of hormone secretion (3, 49–53). However, due to their multifaced cell-to-cell variability, the resulting multicellular activity is not completely synchronized and manifests itself in the form of non-stationary waves, guided by the intrinsic islet heterogeneity (54–59). In recent years, graph-theoretical approaches have been used as a powerful methodology to quantify this non-trivial spatiotemporal behavior within the islets (48, 60–63). By this means, functional connectivity networks are typically constructed based on statistical similarity of Ca2+ traces and previous research has shown that the beta cells form architectures with high communication capacities that are locally segregated into functional subcompartments (63–65). Moreover, beta cell networks are highly heterogeneous and some very well-connected cells, which were identified as metabolically highly active, are believed to be crucial for routing information between cells (66–69). However, the role of different signaling pathways in the coordination of the collective cellular rhythms in islets, and particularly the role of cAMP signaling are incompletely understood. Previous studies suggest that intracellular cAMP concentration influences gap-junctional coupling between beta cells (70–74). Recently, PKA and Epac2A were reported to increase Cx36 gap-junctional coupling by different means (51, 75). Noteworthy, cAMP serves as a strong modulator of intercellular communication and the effect can be controlled by either arm of the incretin signaling cascade (76–79). Therefore, it remains of great interest to explore in more detail how neurohormonal amplifying pathways affect multicellular beta cell activity.
In this study, we simultaneously recorded [Ca2+]IC oscillations in a large number of beta cells in acute pancreas tissue slices with single-cell resolution to assess the role of Epac2A knockout during activation, activity, and deactivation, as well as the effects on intercellular coupling assessed by functional connectivity analyses. We further analyzed the effects of the cAMP-elevating diterpene forskolin on [Ca2+]IC and determine the contribution of Epac2A in mediating these effects. To this end, we performed experiments using forskolin in otherwise substimulatory as well as in stimulatory glucose conditions, in both Epac2A KO mice and their WT littermates.
Materials and methods
Ethics statement
The study was conducted as per the protocol approved by the Administration for Food Safety, Veterinary Sector and Plant Protection of the Republic of Slovenia (permit numbers: U3440-1-61/2009/2 and U34401-12/2015/3) and complied with all national and European recommendations relating to care and work with experimental animals to minimize animal discomfort.
Animals, tissue slice preparation and dye loading
Epac2A KO mice (Rapgef4tm1.1Sse) were generated by Professor Susumu Seino of Kobe University as described previously (13). In short, a loxP site was inserted into exon 1 and a floxed neo cassette was inserted downstream of exon 1. Cre-mediated recombination was used to delete a portion of exon 1. The absence of protein product in pancreatic islet cells was confirmed by RT-PCR analysis. Epac2A HET mice with C57BL/6 background were maintained to obtain Epac2A WT and Epac2A KO littermates. The study was conducted on 5 Epac2A WT and 5 Epac2A KO mice.
Acute pancreas tissue slices were obtained from 12-24 weeks old Epac2A KO mice and their WT littermates of either sex were kept in individually ventilated cages (Allentown LLC, USA) with a 12:12 hours light: dark schedule, as described previously (43, 80). In brief, mice were sacrificed using high concentration of CO2, and the abdominal cavity was accessed via laparotomy. The pancreas was injected with low-melting point 1.9% agarose (Lonza, USA) kept at 40°C dissolved in extracellular solution (ECS, consisting of (in mM) 125 NaCl, 26 NaHCO3, 6 glucose, 6 lactic acid, 3 myo-inositol, 2.5 KCl, 2 Na-pyruvate, 2 CaCl2, 1.25 NaH2PO4, 1 MgCl2, 0.5 ascorbic acid) with access through the proximal common bile duct, which was clamped distally at the major duodenal papilla. Following injection, the pancreas was cooled with ice-cold ECS and extracted from the abdominal cavity. A vibratome (VT 1000 S, Leica) was used to create 140 µm thick tissue slices, which were collected at room temperature in HEPES-buffered saline (HBS, consisting of (in mM) 150 NaCl, 10 HEPES, 6 glucose, 5 KCl, 2 CaCl2, 1 MgCl2; titrated to pH=7.4 using 1 M NaOH). The slices were incubated in the dye-loading solution (6 µM Oregon Green 488 BAPTA-1 AM (OGB-1, Invitrogen), 0.03% Pluronic F-127 (w/v), and 0.12% dimethylsulphoxide (v/v) dissolved in HBS) for 50 minutes at RT. All chemicals were obtained from Sigma-Aldrich (St. Louis, Missouri, USA) unless otherwise specified.
Stimulation protocol and calcium imaging
Individual tissue slices were exposed to single square pulse-like stimulation per islet in a perifusion system containing carbogenated ECS at 37°C. Regardless of the number of islets observed on an individual tissue slice, each slice was stimulated with only a single stimulatory protocol. In control experiments, slices were exposed to 12 mM glucose for 25 minutes, followed by incubation in a solution with substimulatory 6 mM glucose concentration until apparent cessation of activity. In experiments using forskolin slices were exposed to 6 mM or 12 mM glucose for 5 and 15 min, respectively, followed by the addition of 10 µM forskolin for 10 min and the subsequent substimulatory 6 mM glucose concentration again until the deactivation of cells. Calcium imaging was performed on a Leica TCS SP5 DMI6000 CS inverted confocal system (20X HC PL APO water/oil immersion objective, NA 0.7) and a Leica TCS SP5 AOBS Tandem II upright confocal system (20x HCX APO L water immersion objective, NA 1.0). Acquisition frequency was set to 1-2 Hz at 512 x 512 pixels. For a more precise quantification of [Ca2+]IC oscillations in the phases of sustained activity, a resolution of 10 Hz at 256 x 256 pixels was maintained throughout the whole stimulation protocol. The laser power was adapted to avoid photobleaching and prolong the maximum time of recording, while maintaining a satisfactory signal strength. Cells were imaged at the approximate depth of 15 µm or more to avoid recording from cells at the surface. The optical section thickness was kept at about 4 µm to assure recordings from a single cell only. OGB-1 was excited by an argon 488 nm laser line and emitted fluorescence was detected by Leica HyD hybrid detector in the range of 500-700 nm (all from Leica Microsystems, Germany), as described previously (43, 46, 81).
Processing of Ca2+ signals and beta cell functional network analysis
Individual ROIs were selected manually with respect to cell morphology using a high-resolution image, maximal projection, or frame average. Time series were exported using custom-made software (ImageFiltering, copyright Denis Špelič). Custom-made Matlab and Phyton scripts were used for subsequent off-line analysis. Fluorescence signals were expressed as F/F0, the ratio of the fluorescence signal (F) at a certain time point of the experiment relative to the initial fluorescence (F0). A combination of exponential and linear fit was applied to correct the time series data for photobleaching. Time series with distorted signals due to movement artefacts were removed from further analysis. Also excluded were all series with evident non-beta cell-like patterns, such as low-glucose-active and high-glucose-inactive cells, or constantly active or sporadically active with various oscillation frequencies, suggestive of other islet cell types (82–86). Beginnings of [Ca2+]IC increases and decreases after exposure to high and low glucose were used to manually select activation and deactivation times, respectively.
A zero-lag digital filter was used on Ca2+ traces to extract the classical signaling parameters in the phase of sustained activity, with the purpose of removing noise and possible low frequency baseline variations. Typical cutoff frequencies were determined empirically in ranges of 0.02-0.04 Hz and 1-2 Hz for the lower and upper boundary, respectively, as described previously (62, 87). This was followed by time series binarization so that the values from the onset to the end of individual oscillations were 1, and values between the oscillations were 0. Specifically, this process involved detecting the onsets, peaks, and endings of oscillations in a given time series. The onsets were determined by analyzing the derivative (slopes) of the time series, identifying significant changes, i.e., beginning of Ca2+ increases. The peaks represent the maximum signal intensity within each oscillation. The endings were defined as the points where the signal intensity decays to 50% of the maximum value in a given oscillation. These steps are repeated throughout the whole time series, i.e., for all oscillations. Binarized time series data was then used to characterize the oscillatory activity based on the following parameters: (i) average frequencies, (ii) average durations of oscillations, and (iii) relative active times of individual cells. The latter represents the average fraction of time that cells spend in an active state with increased Ca2+ and serves as a suitable measure of the cell’s overall activity (81). The variability of the activation delays was measured by calculating delays relative to the first responding cells in an islet (any-cell-first-responder delays) or to account for the inter-islet differences, dividing the data with the respective islet median delay (relative any-cell-first-responder delays).
To assess the collective beta cell activity in each islet, we generated functional connectivity networks. By this means, nodes represent individual beta cells, and their locations correspond to physical positions of cells in tissue slices. Connections between node pairs were created based on the temporal similarity of the measured Ca2+ dynamics, as evaluated through the calculation of the correlation coefficient, as explained previously (65). We have used variable thresholds to extract the connectivity matrix, so that the average node degree in each islet, i.e., the average number of connections per cell <k>, in the first part in the sustained activity phase (12 mM glucose only) was <k>=8 and the same connectivity threshold was then used in the second part in the phase of sustained activity (either 12 mM glucose or 12 mM glucose + 10 µM forskolin) in the same islet. In this way a comparison between different islets was facilitated as well as a direct evaluation of the changes between both intervals. The extracted functional networks were analyzed with common network metrics. Specifically, the average correlation coefficient along with the average node degree were used to evaluate the average level of beta cell synchronicity in the given islet. Modularity was used to characterize the level of functional segregation, i.e., the extent of partitioning into smaller subpopulations. It measures how well the nodes in a network can be clustered into communities and is computed as the fraction of the edges that exist within communities minus the expected fraction of edges that would exist within communities if the network were random. The values of modularity span between 0 and 1, where 0 means that the nodes in the network are not clustered at all, while a value of 1 means that the nodes are perfectly divided into functional submodules. Moreover, we computed the relative largest component to quantify the level of the network’s functional integration. The largest component of a network is the connected subgraph with the most nodes and is computed by starting with a single node and then iteratively adding all of the nodes that are connected to it. This process is repeated until no more nodes can be added. If the largest component is very small (i.e., close to zero), then the network is likely to be fragmented and difficult to traverse, whilst high values (i.e., close to 1) indicate a very cohesive and integrated architecture. For details, see (63, 65).
Results
A typical beta cell response to glucose stimulation features three subsequent phases of [Ca2+]IC activity, i.e., the activation, the plateau, and the deactivation phase (Figure 1 (43, 81)). To enhance the reader’s comprehension of the recorded signals, an additional time-lapse video was included (Video S1). Below, we present the results according to this order of the three response phases.
Figure 1 A typical beta cell response to glucose and forskolin. (A) 12 mM glucose induced a transient increase in [Ca2+]IC, followed by [Ca2+]IC oscillations during the plateau phase that ceased after stimulus withdrawal. Responses of 4 cells from the same WT islet are shown. Activation (Delay_A) and deactivation (Delay_D) delays were defined as indicated in the figure; (B) [Ca2+]IC activity after addition of 10 µM forskolin in 4 cells from the same WT islet. Rectangles indicate intervals used in subsequent analysis of the plateau phase, the results of which are shown in indicated figures. Note that the prolonged stimulation with glucose only (A) served as a control for possible time effects, since forskolin was added to glucose after a stable plateau was achieved in response to glucose (B); (C) Responses to 12 mM glucose of 4 cells from the same KO islet are shown; (D) [Ca2+]IC activity after addition of 10 µM forskolin in 4 cells from the same KO islet. The analysis methodology was the same as demonstrated in (A, B).
The effect of forskolin and the role of Epac2A in activation of beta cells
Beta cells were inactive in 6 mM glucose in both WT and Epac2A KO mice (Figure 1). We tested whether stimulation by forskolin produced a shift in glucose sensitivity, and whether this effect was Epac2A-dependent (Figures 2, S1). 10 µM forskolin added to 6 mM glucose activated beta cells in a qualitatively similar manner as 12 mM glucose. In quantitative terms, the delays between the stimulus onset and the activation of cells after the forskolin stimulus were smaller in KO mice (243 s vs. 693 s, delay_A). Surprisingly, the absence of Epac2A affected the activation properties in a similar manner even during stimulation with glucose only, as the activation delays were shorter in Epac2A KO mice (112 s vs. 168 s, Figure 2), suggesting a possible inhibitory role of Epac2A during the activation phase at both basal conditions and during forskolin-stimulated conditions.
Figure 2 The effect of forskolin and the role of Epac2A during activation of beta cells from WT and KO mice. (A) The activation delays (Delay_A) after stimulation with either 6 mM glucose + 10 µM forskolin or 12 mM glucose in WT and Epac2A KO mice. Note that the 6 mM glucose failed to elicit a response in beta cells. 1st quartile/median/3rd quartile (Q1/M/Q3, in seconds): 151/693/868 (6 mM glucose + 10 µM forskolin in WTs), 144/243/659 (6 mM glucose + 10 µM forskolin in KOs), 127/168/243 (12 mM glucose in WTs), and 81/112/142 (12 mM glucose in KOs); (B) The any-cell-first-responder delays after stimulation with either 6 mM glucose + 10 µm forskolin or 12 mM glucose in WT and Epac2A KO mice. Q1/M/Q3 (in seconds): 126/278/533 (6 mM glucose + 10 µM forskolin in WTs), 70/178/518 (6 mM glucose + 10 µM forskolin in KOs), 17/35/58 (12 mM glucose in WTs), and 22/48/69 (12 mM glucose in KOs); (C) Relative any-cell-first-responder delays after stimulation with either 6 mM glucose + 10 µm forskolin or 12 mM glucose. Q1/M/Q3: 0.32/0.58/0.90 (6 mM glucose + 10 µM forskolin in WTs), 0.31/0.59/0.95 (6 mM glucose + 10 µM forskolin in KOs), 0.009/0.21/0.40 (12 mM glucose in WTs), and 0.24/0.47/0.69 (12 mM glucose in KO). Data pooled from the following number of cells/islets: 388/6 (6 mM glucose + 10 µM forskolin in WTs), 428/8 (6 mM glucose + 10 µM forskolin in KOs), 1373/24 (12 mM glucose in WTs), 1375/22 (12 mM glucose in KOs). Data were analyzed using Mann-Whitney U test, p values are indicated on graphs.
We have shown previously that the glucose activation of beta cells was relatively heterogeneous (Figure 1), and that the differences in activation delays were strongly glucose-dependent (81). In this vein, we assessed heterogeneity during activation by calculating delays of every cell relative to the first-responding cell in an islet (i.e., any-cell-first-responder delays, delay_A1, Figure 2B). The lack of Epac2A enlarged the any-cell-first-responder delay after stimulation with 12 mM glucose by 26% (48 s in KO vs. 35 s in WT mice), in contrast to a decrease by 36% after stimulation with forskolin plus 6 mM glucose (178 s in KO vs. 278 s in WT mice). Importantly, to account for the inter-islet differences and the shorter activation in Epac2A KO mice, the data were normalized with the respective islet median delays (Figure 2C), the Epac2A deletion increased activation heterogeneity after stimulation with 12 mM glucose (0.47 in KO and 0.21 in WT mice) and did not have a significant effect after stimulation with forskolin plus 6 mM glucose. In sum, Epac2A seems to contribute to a less heterogeneous activation during stimulation with glucose.
The effect of forskolin and the role of Epac2A during the plateau phase
We further explored the effect of stimulation by forskolin and the specific role of Epac2A during the plateau phase by analyzing the classical and network functional parameters.
Classical functional parameters
12 mM glucose evoked repetitive fast [Ca2+]IC oscillations (Figure 1). Relative active time of a cell (defined as the percentage of oscillatory activity that the cells spent at a [Ca2+]IC > 50% of the oscillation amplitude, was 41% at a frequency of 0.033 Hz and an oscillation duration of 10.7 s (Figures 3A, C, E). Epac2A deletion raised the active time by 12% (AT = 46%), by increasing both oscillation duration (11.7 s) and frequency (0.039 Hz).
Figure 3 The effect of forskolin and the role of Epac2A during the plateau phase of response to glucose. (A) The oscillation frequency after stimulation with 12 mM glucose. 1st quartile/median/3rd quartile (Q1/M/Q3) in Hz: 0.025/0.033/0.042 (in WTs) and 0.036/0.039/0.044 (in KOs); (B) The oscillation frequency after either prolonged stimulation with 12 mM or addition of 10 µm forskolin. Q1/M/Q3 in Hz: 0.031/0.033/0.042 (12 mM glucose in WTs), 0.045/0.058/0.069 (12 mM glucose + 10 µM forskolin in WTs), 0.025/0.033/0.056 (12 mM glucose in KOs) and 0.047/0.056/0.064 (12 mM glucose + 10 µM forskolin in KOs); (C) The oscillation duration after stimulation with 12 mM glucose. Q1/M/Q3 in seconds: 9.7/10.7/13.4 (in WTs) and 10.4/11.7/12.7 (in KOs); (D) The oscillation duration after either prolonged stimulation with 12 mM or addition of 10 µm forskolin. Q1/M/Q3 in seconds: 10.8/12.1/13.3 (12 mM glucose in WTs), 7.9/8.5/9.9 (12 mM glucose + 10 µM forskolin in WTs), 9.2/10.9/18.5 (12 mM glucose in KOs), and 9.2/9.9/10.6 (12 mM glucose + 10 µM forskolin in KOs); (E) The relative active time after stimulation with 12 mM glucose. Q1/M/Q3: 0.33/0.41/0.48 (in WTs) and 0.42/0.46/0.51 (in KOs); (F) The relative active time after either prolonged stimulation with 12 mM or addition of 10 µm forskolin. Q1/M/Q3: 0.36/0.43/0.49 (12 mM glucose in WTs), 0.44/0.55/0.59 (12 mM glucose + 10 µM forskolin in WTs), 0.39/0.47/0.53 (12 mM glucose in KOs) and 0.49/0.55/0.62 (12 mM glucose + 10 µM forskolin in KOs). Data pooled from the following number of cells/islets: 687/11 (12 mM glucose in WTs), 656/11 (12 mM glucose + 10 µM forskolin in WTs), 548/7 (12 mM glucose in KOs), and 631/9 (12 mM glucose + 10 µM forskolin in KOs). Data were analyzed using Mann-Whitney U test (for 2 samples) or one-way ANOVA on ranks (Kruskal-Wallis test) followed by Dunn’s multiple comparisons test for more than 2 samples), p values are indicated on graphs.
Next, we tested the effect of forskolin. In WT mice, 10 µM forskolin added to the stimulatory 12 mM glucose raised the active time further by 28% (from 43% to 55%, Figures 3F, S3C). There was a considerable drop in oscillation duration by 30% (8.5 s vs. 12.1 s, Figures 3D, S3B), but active time nevertheless rose due to a large increase in frequency by 75% (0.058 Hz vs. 0.033 Hz, Figures 3B, S3A), corroborating the importance of analyzing the active time rather than the frequency or the durations only (88). In Epac2A KO mice the frequency also increased, oscillation duration decreased, and the active time increased to a similar extent and a similar absolute value as in WT mice (increase by 17%, i.e. from 47% to 55% Figures 3F, S3). These changes in frequency, duration and AT were not due to prolonged exposure to the stimulating glucose concentration alone, as significantly smaller changes were observed between the two time intervals analyzed in cells stimulated with high glucose only (Figures S2, S3). To provide a more detailed insight into how individual islets responded to glucose only and forskolin, we additionally present in Figure S3 how active time, duration, and frequency have changed between intervals, separately for each individual islet.
These results seem to suggest that the absence of Epac2A brings about a minor rise in active time during the response to 12 mM glucose and that pharmacological stimulation with forskolin raises the active time through a large increase in oscillation frequency despite a concomitant decrease in oscillation duration. This modulation of [Ca2+]IC oscillations is preserved in Epac2A KO mice, which implies that it is most probably largely Epac2A-independent.
Functional network parameters
To characterize the collective beta cell [Ca2+]IC activity, we constructed functional connectivity networks, as described in the Methods section. First, we present characteristic beta cell networks in islets from WT and Epac2A KO mice subjected to 12 mM glucose stimulation only (Figures 4A, C) and in islets from both types of mice which were subsequently stimulated with 10 µM forskolin (Figures 4B, D). Connectivity maps are presented separately for all four scenarios. After prolonged stimulation with 12 mM glucose only, the networks became much sparser in both WT and KO mice. In contrast, in the second protocol where 10 µM forskolin was added to 12 mM glucose, the networks remained integral, and the number of connections stayed approximately the same. These qualitative observations in the exemplary islets indicate that forskolin rather profoundly affected the collective beta cell behavior by preventing a decline in or even increasing intercellular connectivity, probably to a comparable extent in both types of mice.
Figure 4 The effect of forskolin on functional beta cell connectivity networks and the role of Epac2A. (A-D) Characteristic functional network representations of islets from WT and Epac2A KO mice stimulated either with 12 mM glucose or with 12 mM glucose + 10 µM forskolin. Nodes signify positions of beta cells within an islet, connections stand for functional associations in Ca2+ activity, and colors of cells denote different communities. In all four panels the left networks correspond to the first interval in the plateau phase, whilst the right networks represent connectivity patterns in the second interval of sustained activity [i.e., either prolonged glucose stimulation (A, C) or stimulation with glucose + forskolin (B, D)]; (E-K) Synchronization and network metrics for the average pooled data from all islets in the given group: average node degree (E), average correlation coefficient (F, G), modularity (H, I), relative largest component (J, K). In (F, H, J) the values in the first part of the plateau phase (glucose only) are shown, whereas in (E, G, I, K) the values corresponding to the second part (either glucose only or glucose+forskolin) are presented. Individual dots represent the average values in individual islets, whereas the horizontal lines denote the median value.
To gain a more general and quantitative insight, we show in Figures 4E-K different metrics of synchronous beta cell dynamics and network parameters pooled from all islets. We first present the average node degrees, i.e., the average number of functional connections. Since we used an adaptive thresholding approach to set the average degree in the first part in the plateau phase to <k>=8 in all islets, we show only the values in the second interval. Results in Figure 4E reveal that in both WT and KO mice the average node degree was profoundly lower (median k≈2) after prolonged stimulation with 12 mM glucose only. In the stimulation protocol with forskolin, the node degree in the second part of the plateau phase was on average higher in both types of mice when compared to stimulation with glucose only, with median values k≈10 and k≈5 for the islets from WT and KO mice, respectively. However, the difference between the control and the forskolin protocol was significant only for WT mice.
Next, we show the average correlation coefficients (Figures 4F, G), modularity values (Figures 4H, I), and relative largest components (Figures 4J, K) of beta cell networks. For these metrics, we present first their absolute values in the first part of the plateau phase, separately for the islets from WT and KO mice (Figures 4F, H, J). None of these parameters were found to differ between the islets from the WT and KO mice, indicating that under stimulation with 12 mM glucose only, Epac2A does not notably affect either the average synchronicity or the beta cell network’s functional integrity and segregation patterns. We proceeded with examining the differences between the glucose stimulation only protocol and the forskolin protocol by comparing the parameters in the second part of the plateau phase. In the islets which were subsequently stimulated with 10 µM forskolin, the average correlation was higher when compared to prolonged 12 mM glucose stimulation only. Furthermore, the networks were less segregated, and the level of functional integration was greater, as reflected by lower values of modularity (Figure 4I) and higher values of the relative largest component (Figure 4K), respectively. However, while all the above differences were significant in islets from WT mice, the effect of forskolin was less well-pronounced in islets from Epac2A KO mice and not statistically significant for any of the parameters, suggesting that Epac2A probably importantly contributes to the greater functional integration of cell networks during stimulation by forskolin. To provide a more detailed insight, we additionally present in Figure S4 how the beta cell synchronization and different network parameters have changed between intervals, separately for each individual islet.
The effect of forskolin and the role of Epac2A during deactivation of beta cells
After cessation of stimulation with 12 mM glucose, there was a delay before beta cells became inactive (Figures 5A, S5). This delay did not differ between WT and KO mice. Additional stimulation by forskolin significantly prolonged the activity of cells (in WTs by 22%, from 348 s to 426 s, Figure 5A). This prolongation was not Epac2A-dependent, as it was even significantly longer in KO mice (from 297 to 820s, Figure 5A). Thus, at least under conditions of additional stimulation by forskolin, Epac2A seems to restrict the activity of beta cells upon removal of the stimulus. In accordance with our previous studies, comparing deactivation delays between different cells (Figure 5B) demonstrated relatively large intercellular differences (43, 64, 81). Stimulation by forskolin enlarged the first-cell-any-cell delays (in WTs from 87 s to 214 s, Figure 5B). Although this effect of forskolin did not depend on Epac2A in absolute terms (214 s in WT and 210 s in KO mice, Figure 5B), taking into account the inter-islet variability of the median delays and the absolutely longer delays in KOs revealed that forskolin increased the relative heterogeneity only in WTs (from 0.26 to 0.51, Figure 5C). On the other hand, normalization with median values of delays also showed that the relative heterogeneity may be slightly higher in KOs after stimulation with glucose only (0.21 vs. 0.34, Figure 5C).
Figure 5 The effect of forskolin and the role of Epac2A during deactivation of beta cells from WT and KO mice. (A) The deactivation delays after cessation of stimulation with 12 mM glucose in the presence or absence of 10 µm forskolin. 1st quartile/median/3rd quartile (Q1/M/Q3), in seconds): 197/348/490 (12 mM glucose in WTs), 303/426/586 (12 mM glucose + 10 µM forskolin in WTs), 211/297/479 (12 mM glucose in KOs), and 477/820/861 (12 mM glucose + 10 µM forskolin in KOs); (B) The deactivation delays of cells after 1st deactivated cell after cessation of stimulation with 12 mM glucose in the presence or absence of 10 µm forskolin. Q1/M/Q3 (in seconds): 39/87/158 (12 mM glucose in WTs), 97/214/370 (12 mM glucose + 10 µM forskolin in WTs), 45/83/303 (12 mM glucose in KOs), and 153/210/293 (12 mM glucose + 10 µM forskolin in KOs); (C) Relative deactivation delays of cells after 1st deactivated cell after cessation of stimulation with 12 mM glucose in the presence or absence of 10 µm forskolin. Q1/M/Q3 (in seconds): 0.16/0.26/0.54 (12 mM glucose in WTs), 0.28/0.51/0.68 (12 mM glucose + 10 µM forskolin in WTs), 0.17/0.34/0.66 (12 mM glucose in KOs), and 0.24/0.32/0.58 (12 mM glucose + 10 µM forskolin in KOs). Data pooled from the following number of cells/islets: 678/13 (12 mM glucose in WTs), 640/10 (12 mM glucose + 10 µM forskolin in WT), 567/10 (12 mM glucose in KOs), 422/7 (12 mM glucose + 10 µM forskolin in KOs). Data were analyzed using one-way ANOVA on ranks (Kruskal-Wallis test) followed by Dunn’s multiple comparisons test, p values are indicated on graphs.
Taken together, these results suggest that forskolin increases beta cell activity not only by making the activation more likely and by raising the active time during the plateau phase, but also by prolonging beta cell activity upon removal of the stimulus, extending their activity into an otherwise silent period. This effect of forskolin on deactivation seems to be Epac2A-independent and under normal conditions, Epac2A may even restrict it. Forskolin also made deactivation more heterogenous in absolute terms, but given the huge prolongation of deactivation in KOs, this increase in heterogeneity was relatively larger in WTs. On the other hand, after stimulation with glucose only, the deactivation was in relative terms slightly more heterogenous in KOs.
Discussion
The principal aims of this study were to describe in detail the possible role of Epac2A and the effects of forskolin on [Ca2+]IC oscillations in beta cells during activation, activity, and deactivation, in terms of classical physiological parameters and by using network analyses. First, it should be noted that the differences between Epac2A KO mice and their WT littermates observed in our study are probably not the consequence of altered physiological and biochemical parameters, since previous studies have shown that Epac2A KO mice are generally healthy, with no noticeable physiological abnormalities. Their body weights and food intake are similar as in WT mice, the fasting blood glucose levels and levels of insulin secretion are also similar, and they have normal glucose tolerance and insulin sensitivity (40, 41, 89). On the other hand, the same studies revealed that Epac2A KO mice have larger adipocytes and higher amounts of white adipose tissue, their plasma leptin level is increased, and adiponectin level is decreased already at 7 weeks of age compared with WT mice. Furthermore, these mice are more susceptible to developing obesity when fed a high fat diet. Leptin resistance test performed in vivo showed suppressed hypothalamic leptin signaling (40, 89). Since Epac2A is expressed in hypothalamus and heart, but neither in white and brown adipose tissue nor in skeletal muscles, the defect in hypothalamic leptin signaling likely accounts for the disrupted physiological and biochemical parameters observed in Epac2A KO mice when exposed to high metabolic demands (40, 41). Considering the above evidence, we strongly believe that the differences between Epac2A KO mice and their WT littermates observed in our study are at least for the most part due to intrinsic differences in the neurohormonal amplifying pathway of insulin secretion in beta cells due to the lack of Epac2A.
Second, we showed previously that pancreatic beta cells within acute mouse pancreas tissue slices are typically silent in 6 mM glucose and respond to stimulatory glucose concentrations above 7 mM with a transient increase of [Ca2+]IC followed by fast [Ca2+]IC oscillations (43, 57, 81, 88). In the present work, we demonstrated that beta cells from both Epac2A KO and WT mice responded to the otherwise substimulatory glucose (6 mM) after stimulation by forskolin. This indicates that a rise in cAMP due to adenylyl cyclase activation can affect [Ca2+]IC, probably in a PKA-dependent and Epac2A-independent manner, as already described before (15, 17, 90). Several studies demonstrated that this could involve an increased cytosolic ATP/ADP ratio, either by accelerated ATP synthesis, decreased ATP degradation, or by redistribution of ATP between different intracellular pools (91), decreased KATP conductance (92), inhibition of K+ currents through voltage-dependent KV1.4 channels (17, 90), enhancement of Ca2+ signals through VDCC and/or Ca2+ mobilization from internal stores via RyR or IP3R (19–22). More studies are required to identify the molecular targets responsible for the above effects of cAMP during activation and quantify their contributions, as well as to determine possible cAMP-independent effects of forskolin. Although forskolin is a potent adenylyl cyclase activator that successfully evoked oscillations of [Ca2+]IC in 6 mM glucose, the effect of forskolin in 6 mM glucose seems to be significantly weaker compared to the stimulation with 12 mM glucose, since we observed a much shorter activation delay in 12 mM glucose. This supports previously published data showing progressively shorter activation delay with increasing glucose concentrations (43, 57, 81, 88), a phenomenon most probably resulting from the amount of metabolized glucose needed to activate beta cells or direct cAMP production (93, 94). Furthermore, while in 12 mM glucose most beta cells within an islet responded to stimulation, when stimulated with 6 mM glucose + forskolin, only a minority of beta cells became active (data not shown), demonstrating large differences in their metabolic activity and therefore their sensitivity to glucose also under stimulation by forskolin. Since in our protocol beta cells were exposed to forskolin only for 10 min, we wish to point out the possibility that the majority of beta cells within an islet could eventually be recruited and would respond to prolonged stimulation with forskolin, as seen during exposure to lower stimulatory concentrations of glucose (81, 88). Interestingly, beta cells from Epac2A KO mice responded to both protocols with a shorter activation delay compared to their WT littermates (Figure 2). As mentioned before, cAMP can affect [Ca2+]IC in a PKA-dependent (15, 17, 90) or an Epac2A-dependent manner (19–22, 33). In Epac2A KO cells, all the available cAMP produced either by forskolin-enhanced adenylyl cyclase activity in substimulatory glucose conditions or endogenously synthesized in beta cells after a high glucose load, will act only through PKA-dependent pathways, increasing for instance the L-type VDCC activity after phosphorylation with PKA (19, 20, 95). Our results point to the possibility that the role of Epac2A may be at least partly inhibitory during activation, that in Epac2A KO mice there is a compensatory upregulation of the stimulatory PKA-dependent mechanisms, or a combination of both. However, further studies are required to confirm and further clarify this mechanism. The activation delays among islets as well as among beta cells within the same islet were very heterogeneous, as observed from the delays between the first-responding cell and the others from the same islet in both protocols (Figures 2, S1). Most previous studies confirm this, showing at least some degree of heterogeneity between cells during their activation (72, 81, 88, 93, 96–101). Here we showed that in both WT and KO cells, the heterogeneity of delays decreased with faster activation, enabling beta cells a more homogeneous response to a stronger stimulation (i.e., 12 mM glucose vs. 6 mM glucose + 10 µM forskolin). Interestingly, the relative heterogeneity of activation delays was higher in beta cells from Epac2A KO mice in the high glucose regime (Figure 2C), suggesting the involvement of Epac2A in the coordination of activation under a high glucose load. The latter could be mediated through Cx36 gap junction coupling in an Epac2A-dependent manner to overcome the extensive intrinsic heterogeneity present in beta cells and to ensure a more coordinated response to glucose (78, 81). This view is also consistent with the above finding that the responses are faster in Epac2A KO mice, since in more weakly coupled syncytia, the intrinsically more sensitive cells can escape the inhibition from less responsive cells (57, 102, 103). Most importantly, our findings regarding network and deactivation properties also argue for the role of Epac2A in intercellular coupling (see below) and are consistent with another recent report (51), but deserve to be explored further in the future.
Next, the plateau phase of response after stimulation with high glucose consists of repetitive fast [Ca2+]IC oscillations that reflect bursts of membrane potential depolarizations (46, 97, 104, 105) or Ca2+ release from intracellular Ca2+ stores (106). Epac2A ablation increased the beta cell active time to a minor extent (Figure 3) that is most probably biologically irrelevant, since the Epac2 KO animals display normal insulin and glucose levels, as discussed above. In our hands, the absence of Epac2A did not prevent the forskolin-mediated rise in frequency of oscillations and active time (Figure 3). A similar effect was observed for oscillations of membrane potential in GLP-1, forskolin, or glucagon-activated cells under high glucose conditions (92, 107–109). Interestingly, the cells that were previously less active increased their oscillation frequency the most, while the duration of oscillations decreased slightly in the majority of cells (Figure S2). This response to forskolin clearly shows that the activation of the neurohormonal amplifying pathway through cAMP involves different molecular mechanisms compared to glucose stimulation, which results in longer oscillation duration at the cost of reduced oscillation frequency or in increased frequency, but without a significant drop in oscillation duration (48, 62, 81, 88). The potentiation of [Ca2+]IC oscillations under the influence of forskolin is probably exerted independently of KATP channels, since GLP-1 is able to increase insulin secretion also in Kir6.2 KO mice that lack functional KATP channels (110) and probably involves increased Ca2+ influx through L-type VDCCs, phosphorylated by PKA (18) as well as PKA- and Epac2A-dependent mobilization of Ca2+ from internal stores (19–22, 111). In addition to the possible direct effects of cAMP on [Ca2+]IC, there is evidence that the amplifying effect of cAMP could be explained by a larger gap junctional conductance between beta cells (70). In beta cell syncytium, electrical coupling through Cx36 plays a crucial role in regulating coordinated glucose induced Ca2+ oscillations, thereby determining the dynamics of insulin secretion (75, 112, 113). In the present paper, we demonstrated that under stimulatory glucose levels, electrical coupling results in well-coordinated Ca2+ oscillations, but prolonged exposure to glucose causes an Epac2A-independent beta cell desynchronization, lower network integrity and higher segregation (Figure 4), as described previously (114). On the other hand, forskolin can prevent this decline in network function, but in an at least partly Epac2A-dependent manner, since in Epac2A KO mice, the improvement in network parameters was only partial. Disrupted cellular communication was described before in metabolically overloaded cells (48, 51, 67, 72, 115), while the GLP-1R agonist Exendin-4 increased Cx36 coupling and improved Ca2+ oscillation coordination (51). Several mechanisms have been proposed to explain how cAMP regulates Cx36, either by changing Cx36 gene expression, increasing Cx36 coupling, or by changing distribution of Cx36 on the cell membrane (75, 116, 117). In the retina, cAMP has been shown to regulate gap junction coupling in a PKA-dependent manner (118), with no effect on trafficking or changing the distribution of Cx36 on the plasma membrane (76, 119), while in myocardial cells, PKA is responsible for opening of Cx36 (120). On the other hand, in neurons, regulation of Cx36 is Epac2A-mediated (121). In pancreatic beta cells, both PKA and Epac2A seem to be responsible for Cx36 regulation. PKA was proposed as a channel gating regulator, while Epac2A probably influences Cx36 coupling via slower mechanisms, such as trafficking, assembly, or turnover (51), similarly as in rat myocardial cells (77, 78). Overall, our data suggest that the presence of Epac2A is not critical for basal beta cell functional network integrity upon stimulation with glucose, as no notable differences in beta cell synchronicity and network connectivity were observed when comparing the behavior of control and Epac2A KO islets. However, when the islets were additionally stimulated with forskolin, cells from Epac2A KO mice failed to fully exploit the positive effects of elevated cAMP on beta cell network activity. This might be a consequence of the Epac2A-related deficiency of the connexon trafficking pathways (51, 75, 79, 122). However, additional studies are required to determine the effects of Epac2A deficiency on connexon trafficking in beta cells and to elucidate whether this is the main factor that affects the beta cell network dynamics when cAMP is increased.
Finally, after cessation of stimulation, the oscillatory activity in beta cells gradually stops, and [Ca2+]IC returns to baseline. Proper deactivation is especially important to prevent hypoglycemia (102, 123). Like activation, the deactivation phase is also glucose-dependent (81, 88). After decreasing glucose from 12 mM to the substimulatory level (6 mM), beta cells deactivated with 2-3 times greater time lags compared to the activation phase, with no differences between WT and Epac2A KO mice (Figure 5). Again, considerable heterogeneity was observed among cells, which is in accordance with previously published studies (43, 64, 81, 88). When beta cells were additionally stimulated by forskolin, the deactivation delay was prolonged significantly, similarly to what can be seen under an extremely high glucose load (81, 88), and this prolongation was especially pronounced in cells lacking Epac2A. Longer and more heterogenous deactivation delays probably indicate a greater degree of activation during the preceding stimulation period, not only through the triggering pathway, but also through the neurohormonal amplifying pathway and probably elevated cAMP levels. It is reasonable to speculate that following strong stimulation, every beta cell will need longer to decrease the concentration of the triggering and amplifying secondary messengers below the stimulatory level. To speculate even further, we believe that slight changes in Cx36 coupling in Epac2A KO mice could at least partly account for the faster and more heterogeneous activation, decreased responsiveness of network parameters to cAMP, and also the longer deactivation delay, since stronger intercellular coupling during deactivation is expected to bring about a stronger hyperpolarizing influence from the already deactivated cells on the still active cells. This hypothesis is a good starting point for future work that should also test the translational relevance of our findings in mouse models of diabetes mellitus and for the human islets, and aim at measuring cAMP dynamics simultaneously with [Ca2+]IC. With regard to the latter, a concentration of forskolin equal to the concentration used in this study tripled cAMP levels in isolated MIN6 cells (from ~2,5 to 7.5 ng/mg total protein (124)) and rat beta cells (from ~ 0,7 to 2,1 pmol/ml (125)). Experimental quantification of cAMP levels would also enable us to demonstrate explicitly that forskolin indeed elevated cAMP. To our knowledge, there are currently two approaches to measure cAMP levels: lysis of isolated islets (124, 125) and transfection with biosensor-encoding adenoviruses (126). While the first approach is not feasible in the tissue slice (lysis of the tissue slice would include contamination with other cell types, such as acinar and ductal cells), the second approach requires a cultivation period that can critically alter the normal glucose response (127). Based on previously published results, we can thus only assume that cAMP and [Ca2+]IC dynamics are tightly correlated (126, 128). We also wish to point out that the most distal part in the stimulus-secretion coupling cascade, i.e., insulin secretion, was not addressed in this study. However, we and others showed previously that cAMP increases the sensitivity of insulin granule fusion or their availability (11, 23–25). The positive Epac2A-independent effect on [Ca2+]IC dynamics demonstrated in this study is therefore presumably accompanied also by the effect on the distal steps of the stimulus-secretion coupling cascade. In the future, simultaneous or complementary measurements of secretion dynamics could help us understand the respective contributions of the effects of cAMP on [Ca2+]IC and increased sensitivity of the exocytotic apparatus at the level of insulin secretion and determine whether they are additive or synergistic.
Conclusions
Forskolin was able to weakly activate beta cells exposed to substimulatory glucose and beta cells from Epac2A KO mice responded faster than beta cells from WT littermates to both high glucose and to forskolin added to low glucose. During the plateau phase, forskolin added to high glucose resulted in increased oscillation frequency and relative active time, with well-coordinated activity among beta cells and this response did not critically depend on Epac2A. Prolonged exposure to glucose caused Epac2A-independent beta cell desynchronization, lower network integrity, and higher segregation, while activation of the neurohormonal amplifying pathway prevented this decline in network function in an at least partly Epac2A-dependent manner. In the end, following stimulation with high glucose and forskolin, beta cells deactivated more slowly and the prolongation of activity into the otherwise already silent period was especially well pronounced in beta cells from Epac2A KO mice. Taken together, our results suggest that especially under conditions of stimulated cAMP production, Epac2A may play a role in coordinating beta cell collective activity, with its absence resulting in earlier activation, weaker functional coupling during activity, and later deactivation.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://doi.org/10.6084/m9.figshare.23723700.v2.
Ethics statement
The animal study was approved by Administration for Food Safety, Veterinary Sector and Plant Protection of the Republic of Slovenia. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
Conceptualization, MSK, AS, and MSR; methodology, MSK, JD, MG and AS; software, JD and MG; validation, MSK and AS; formal analysis, MSK, JD and MG; investigation, MSK, JD, LK and VP; resources, MSK, JD and MG; data curation, MSK, JD and MG; writing—original draft preparation. MSK, JD, LK, VP, MG and AS; writing—review and editing, AS and MSK; visualization, MSK and MG; supervision, AS and MSR; project administration, MSK, AS and MSR; funding acquisition, AS and MSR. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by SLOVENIAN RESEARCH AGENCY, grant number P3-0396, I0-0029, N3-0048, N3-0133, N3-170, J3-3077, and J3-9289.
Acknowledgments
We thank Rudi Mlakar for his excellent technical assistance. We also thank Professor Susumo Seino, Kobe University Graduate School of Medicine, Kobe, Japan and Professor Martina Schmidt, University of Groningen, Groningen, the Netherlands, for providing Epac2A KO mice to establish our colony.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2023.1225486/full#supplementary-material
Supplementary Video 1 | A time series illustrating the response of pancreatic islet cells to 12 mM glucose. Combined movie of average calcium signal (left) and F/F0 (right). F0 was defined as the first 100 frames where no activity was detected. Scale bar = 50 µm. Time stamps were added to the upper right corner.
References
2. Weyer C, Bogardus C, Mott DM, Pratley RE. The natural history of insulin secretory dysfunction and insulin resistance in the pathogenesis of type 2 diabetes mellitus. J Clin Invest (1999) 104:787–94. doi: 10.1172/JCI7231
3. Skelin Klemen M, Dolenšek J, Slak Rupnik M, Stožer A. The triggering pathway to insulin secretion: Functional similarities and differences between the human and the mouse β cells and their translational relevance. Islets (2017) 9:109–39. doi: 10.1080/19382014.2017.1342022
4. Seino S, Sugawara K, Yokoi N, Takahashi H. beta-Cell signalling and insulin secretagogues: A path for improved diabetes therapy. Diabetes Obes Metab (2017) 19(Suppl 1):22–9. doi: 10.1111/dom.12995
5. Ashcroft FM, Ashcroft SJH. Mechanism of insulin secretion. In: Ashcroft FM, Ashcroft SJH, editors. Insulin: molecular biology to pathology. Oxford, UK: IRL Press at Oxford University Press (1992). p. 97–150.
6. Rorsman P, Renström E. Insulin granule dynamics in pancreatic beta cells. Diabetologia (2003) 46:1029–45. doi: 10.1007/s00125-003-1153-1
7. Rorsman P, Ashcroft FM. Pancreatic β-cell electrical activity and insulin secretion: of mice and men. Physiol Rev (2018) 98:117–214. doi: 10.1152/physrev.00008.2017
8. Dou H, Wang C, Wu X, Yao L, Zhang X, Teng S, et al. Calcium influx activates adenylyl cyclase 8 for sustained insulin secretion in rat pancreatic beta cells. Diabetologia (2015) 58:324–33. doi: 10.1007/s00125-014-3437-z
9. Tengholm A, Gylfe E. cAMP signalling in insulin and glucagon secretion. Diabetes Obes Metab (2017) 19(Suppl 1):42–53. doi: 10.1111/dom.12993
10. Trexler AJ, Taraska JW. Regulation of insulin exocytosis by calcium-dependent protein kinase C in beta cells. Cell Calcium (2017) 67:1–10. doi: 10.1016/j.ceca.2017.07.008
11. Renström E, Eliasson L, Rorsman P. Protein kinase A-dependent and -independent stimulation of exocytosis by cAMP in mouse pancreatic B-cells. J Physiol (1997) 502:105–18. doi: 10.1111/j.1469-7793.1997.105bl.x
12. Eliasson L, Ma X, Renstrom E, Barg S, Berggren P-O, Galvanovskis J, et al. SUR1 regulates PKA-independent cAMP-induced granule priming in mouse pancreatic B-cells. J Gen Physiol (2003) 121:181–97. doi: 10.1085/jgp.20028707
13. Shibasaki T, Takahashi H, Miki T, Sunaga Y, Matsumura K, Yamanaka M, et al. Essential role of Epac2/Rap1 signaling in regulation of insulin granule dynamics by cAMP. Proc Natl Acad Sci (2007) 104:19333–8. doi: 10.1073/pnas.0707054104
14. Henquin JC, Nenquin M. Activators of PKA and Epac distinctly influence insulin secretion and cytosolic Ca2+ in female mouse islets stimulated by glucose and tolbutamide. Endocrinology (2014) 155:3274–87. doi: 10.1210/en.2014-1247
15. Light PE, Manning Fox JE, Riedel MJ, Wheeler MB. Glucagon-like peptide-1 inhibits pancreatic ATP-sensitive potassium channels via a protein kinase A- and ADP-dependent mechanism. Mol Endocrinol (2002) 16:2135–44. doi: 10.1210/me.2002-0084
16. MacDonald PE, Salapatek AMF, Wheeler MB. Temperature and redox state dependence of native Kv2.1 currents in rat pancreatic beta-cells. J Physiol (2003) 546:647–53. doi: 10.1113/jphysiol.2002.035709
17. Kim SJ, Choi WS, Han JS, Warnock G, Fedida D, McIntosh CH. A novel mechanism for the suppression of a voltage-gated potassium channel by glucose-dependent insulinotropic polypeptide: protein kinase A-dependent endocytosis. J Biol Chem (2005) 280:28692–700. doi: 10.1074/jbc.M504913200
18. Gromada J, Holst JJ, Rorsman P. Cellular regulation of islet hormone secretion by the incretin hormone glucagon-like peptide 1. Pflügers Arch (1998) 435:583–94. doi: 10.1007/s004240050558
19. Ammala C, Ashcroft FM, Rorsman P. Calcium-independent potentiation of insulin release by cyclic AMP in single [beta]-cells. Nature (1993) 363:356–8. doi: 10.1038/363356a0
20. Kanno T, Suga S, Wu J, Kimura M, Wakui M. Intracellular cAMP potentiates voltage-dependent activation of L-type Ca2+ channels in rat islet beta-cells. Pflugers Archiv: Eur J Physiol (1998) 435:578–80. doi: 10.1007/s004240050556
21. Chepurny OG, Kelley GG, Dzhura I, Leech CA, Roe MW, Dzhura E, et al. PKA-dependent potentiation of glucose-stimulated insulin secretion by Epac activator 8-pCPT-2’-O-Me-cAMP-AM in human islets of Langerhans. Am J Physiol Endocrinol Metab (2010) 298:E622–33. doi: 10.1152/ajpendo.00630.2009
22. Kang G, Chepurny OG, Rindler MJ, Collis L, Chepurny Z, Li W-h, et al. A cAMP and Ca2+ coincidence detector in support of Ca2+-induced Ca2+ release in mouse pancreatic β cells. J Physiol (2005) 566:173–88. doi: 10.1113/jphysiol.2005.087510
23. Wan QF, Dong Y, Yang H, Lou X, Ding J, Xu T. Protein kinase activation increases insulin secretion by sensitizing the secretory machinery to Ca2+. J Gen Physiol (2004) 124:653–62. doi: 10.1085/jgp.200409082
24. Dolensek J, Skelin M, Rupnik MS. Calcium dependencies of regulated exocytosis in different endocrine cells. Physiol Res (2011) 60:S29–38. doi: 10.33549/physiolres.932176
25. Skelin M, Rupnik M. cAMP increases the sensitivity of exocytosis to Ca2+ primarily through protein kinase A in mouse pancreatic beta cells. Cell Calcium (2011) 49:89–99. doi: 10.1016/j.ceca.2010.12.005
26. Holz GG. Epac: A new cAMP-binding protein in support of glucagon-like peptide-1 receptor-mediated signal transduction in the pancreatic beta-cell. Diabetes (2004) 53:5–13. doi: 10.2337/diabetes.53.1.5
27. Ozaki N, Shibasaki T, Kashima Y, Miki T, Takahashi K, Ueno H, et al. cAMP-GEFII is a direct target of cAMP in regulated exocytosis. Nat Cell Biol (2000) 2:805–11. doi: 10.1038/35041046
28. Seino S, Shibasaki T. PKA-dependent and PKA-independent pathways for cAMP-regulated exocytosis. Physiol Rev (2005) 85:1303–42. doi: 10.1152/physrev.00001.2005
29. Kawasaki H, Springett GM, Mochizuki N, Toki S, Nakaya M, Matsuda M, et al. A family of cAMP-binding proteins that directly activate Rap1. Science (1998) 282:2275–9. doi: 10.1126/science.282.5397.2275
30. Kelley GG, Chepurny OG, Schwede F, Genieser HG, Leech CA, Roe MW, et al. Glucose-dependent potentiation of mouse islet insulin secretion by Epac activator 8-pCPT-2’-O-Me-cAMP-AM. Islets (2009) 1:260–5. doi: 10.4161/isl.1.3.9645
31. Tsuboi T, Rutter GA. Insulin secretion by ‘kiss-and-run’ exocytosis in clonal pancreatic islet beta-cells. Biochem Soc Trans (2003) 31:833–6. doi: 10.1042/bst0310833
32. Park J-H, Kim S-J, Park S-H, Son D-G, Bae J-H, Kim HK, et al. Glucagon-like peptide-1 enhances glucokinase activity in pancreatic β-cells through the association of epac2 with rim2 and rab3A. Endocrinology (2012) 153(2):574–82. doi: 10.1210/en.2011-0259
33. Kang G, Chepurny OG, Malester B, Rindler MJ, Rehmann H, Bos JL, et al. cAMP sensor Epac as a determinant of ATP-sensitive potassium channel activity in human pancreatic beta cells and rat INS-1 cells. J Physiol (2006) 595–609. doi: 10.1113/jphysiol.2006.107391
34. Kang G, Leech CA, Chepurny OG, Coetzee WA, Holz GG. Role of the cAMP sensor Epac as a determinant of KATP channel ATP sensitivity in human pancreatic β-cells and rat INS-1 cells. J Physiol (2008) 586(Pt 5):1307–19. doi: 10.1113/jphysiol.2007.143818
35. Dzhura I, Chepurny OG, Leech CA, Roe MW, Dzhura E, Xu X, et al. Phospholipase C-ϵ links Epac2 activation to the potentiation of glucose-stimulated insulin secretion from mouse islets of Langerhans. Islets (2011) 3:121–8. doi: 10.4161/isl.3.3.15507
36. Kang G, Chepurny OG, Holz GG. cAMP-regulated guanine nucleotide exchange factor II (Epac2) mediates Ca2+-induced Ca2+ release in INS-1 pancreatic beta-cells. J Physiol (2001) 536:375–85. doi: 10.1111/j.1469-7793.2001.0375c.xd
37. Kang G, Joseph JW, Chepurny OG, Monaco M, Wheeler MB, Bos JL, et al. Epac-selective cAMP analog 8-pCPT-2′-O-me-cAMP as a stimulus for ca2+-induced ca2+ Release and exocytosis in pancreatic β-cells. J Biol Chem (2003) 278:8279–85. doi: 10.1074/jbc.M211682200
38. Shibasaki T, Sunaga Y, Fujimoto K, Kashima Y, Seino S. Interaction of ATP sensor, cAMP sensor, ca2+ Sensor, and voltage-dependent ca2+ Channel in insulin granule exocytosis. J Biol Chem (2004) 279(9):7956–61. doi: 10.1074/jbc.M309068200
39. Yasuda T, Shibasaki T, Minami K, Takahashi H, Mizoguchi A, Uriu Y, et al. Rim2alpha determines docking and priming states in insulin granule exocytosis. Cell Metab (2010) 12:117–29. doi: 10.1016/j.cmet.2010.05.017
40. Hwang M, Go Y, Park JH, Shin SK, Song SE, Oh BC, et al. Epac2a-null mice exhibit obesity-prone nature more susceptible to leptin resistance. Int J Obes (2017) 41:279–88. doi: 10.1038/ijo.2016.208
41. Song WJ, Mondal P, Li Y, Lee SE, Hussain MA. Pancreatic beta-cell response to increased metabolic demand and to pharmacologic secretagogues requires EPAC2A. Diabetes (2013) 62:2796–807. doi: 10.2337/db12-1394
42. Benninger RK, Zhang M, Head WS, Satin LS, Piston DW. Gap junction coupling and calcium waves in the pancreatic islet. Biophys J (2008) 95:5048–61. doi: 10.1529/biophysj.108.140863
43. Stožer A, Dolenšek J, Rupnik MS. Glucose-stimulated calcium dynamics in islets of langerhans in acute mouse pancreas tissue slices. PloS One (2013) 8:e54638. doi: 10.1371/journal.pone.0054638
44. Bosco D, Haefliger J-A, Meda P. Connexins: key mediators of endocrine function. Physiol Rev (2011) 91:1393–445. doi: 10.1152/physrev.00027.2010
45. Šterk M, Dolenšek J, Bombek LK, Markovič R, Zakelšek D, Perc M, et al. Assessing the origin and velocity of Ca2+ waves in three-dimensional tissue: Insights from a mathematical model and confocal imaging in mouse pancreas tissue slices. Commun Nonlinear Sci Numerical Simulation (2021) 93:105495. doi: 10.1016/j.cnsns.2020.105495
46. Dolenšek J, Stožer A, Skelin Klemen M, Miller EW, Slak Rupnik M. The relationship between membrane potential and calcium dynamics in glucose-stimulated beta cell syncytium in acute mouse pancreas tissue slices. PloS One (2013) 8:e82374. doi: 10.1371/journal.pone.0082374
47. Stozer A, Hojs R, Dolensek J. Beta cell functional adaptation and dysfunction in insulin resistance and the role of chronic kidney disease. Nephron (2019) 143:33–7. doi: 10.1159/000495665
48. Gosak M, Yan-Do R, Lin H, MacDonald PE, Stozer A. Ca2+ Oscillations, waves, and networks in islets from human donors with and without type 2 diabetes. Diabetes (2022) 71:2584–96. doi: 10.2337/db22-0004
49. Bavamian S, Klee P, Britan A, Populaire C, Caille D, Cancela J, et al. Islet-cell-to-cell communication as basis for normal insulin secretion. Diabetes Obes Metab (2007) 9:118–32. doi: 10.1111/j.1463-1326.2007.00780.x
50. Benninger RK, Piston DW. Cellular communication and heterogeneity in pancreatic islet insulin secretion dynamics. Trends Endocrinol metabolism: TEM (2014) 25(8):399–406. doi: 10.1016/j.tem.2014.02.005
51. Farnsworth NL, Walter R, Piscopio RA, Schleicher WE, Benninger RKP. Exendin-4 overcomes cytokine-induced decreases in gap junction coupling via protein kinase A and Epac2 in mouse and human islets. J Physiol (2019) 597:431–47. doi: 10.1113/JP276106
52. Zavala E, Wedgwood KCA, Voliotis M, Tabak J, Spiga F, Lightman SL, et al. Mathematical modelling of endocrine systems. Trends Endocrinol Metab (2019) 30:244–57. doi: 10.1016/j.tem.2019.01.008
53. Corezola do Amaral ME, Kravets V, Dwulet JM, Farnsworth NL, Piscopio R, Schleicher WE, et al. Caloric restriction recovers impaired β-cell-β-cell gap junction coupling, calcium oscillation coordination, and insulin secretion in prediabetic mice. Am J Physiol-Endocrinol Metab (2020) 319:E709–20. doi: 10.1152/ajpendo.00132.2020
54. Benninger RK, Hutchens T, Head WS, McCaughey MJ, Zhang M, Le Marchand SJ, et al. Intrinsic islet heterogeneity and gap junction coupling determine spatiotemporal Ca(2)(+) wave dynamics. Biophys J (2014) 107:2723–33. doi: 10.1016/j.bpj.2014.10.048
55. Gosak M, Stozer A, Markovic R, Dolensek J, Perc M, Rupnik MS, et al. Critical and supercritical spatiotemporal calcium dynamics in beta cells. Front Physiol (2017) 8:1106. doi: 10.3389/fphys.2017.01106
56. Westacott MJ, Ludin NWF, Benninger RKP. Spatially Organized beta-Cell Subpopulations Control Electrical Dynamics across Islets of Langerhans. Biophys J (2017) 113:1093–108. doi: 10.1016/j.bpj.2017.07.021
57. Stozer A, Markovic R, Dolensek J, Perc M, Marhl M, Slak Rupnik M, et al. Heterogeneity and delayed activation as hallmarks of self-organization and criticality in excitable tissue. Front Physiol (2019) 10:869. doi: 10.3389/fphys.2019.00869
58. Kravets V, Dwulet JM, Schleicher WE, Hodson DJ, Davis AM, Piscopio RA, et al. Functional architecture of the pancreatic islets: First responder cells drive the first-phase [Ca2+] response. bioRxiv (2020) 2020:12.22.424082. doi: 10.1101/2020.12.22.424082
59. Benninger RKP, Kravets V. The physiological role of β-cell heterogeneity in pancreatic islet function. Nat Rev Endocrinol (2022) 18:9–22. doi: 10.1038/s41574-021-00568-0
60. Stožer A, Gosak M, Dolenšek J, Perc M, Marhl M, Rupnik MS, et al. Functional connectivity in islets of langerhans from mouse pancreas tissue slices. PloS Comput Biol (2013) 9:e1002923. doi: 10.1371/journal.pcbi.1002923
61. Salem V, Silva LD, Suba K, Georgiadou E, Neda Mousavy Gharavy S, Akhtar N, et al. Leader beta-cells coordinate Ca(2+) dynamics across pancreatic islets in vivo. Nat Metab (2019) 1:615–29. doi: 10.1038/s42255-019-0075-2
62. Zmazek J, Klemen MS, Markovič R, Dolenšek J, Marhl M, Stožer A, et al. Assessing different temporal scales of calcium dynamics in networks of beta cell populations. Front Physiol (2021) 12:612233–3. doi: 10.3389/fphys.2021.612233
63. Stozer A, Sterk M, Paradiz Leitgeb E, Markovic R, Skelin Klemen M, Ellis CE, et al. From isles of konigsberg to islets of langerhans: examining the function of the endocrine pancreas through network science. Front Endocrinol (2022) 13:922640. doi: 10.3389/fendo.2022.922640
64. Markovic R, Stozer A, Gosak M, Dolensek J, Marhl M, Rupnik MS. Progressive glucose stimulation of islet beta cells reveals a transition from segregated to integrated modular functional connectivity patterns. Sci Rep (2015) 5:7845. doi: 10.1038/srep07845
65. Gosak M, Markovic R, Dolensek J, Slak Rupnik M, Marhl M, Stozer A, et al. Network science of biological systems at different scales: A review. Phys Life Rev (2018) 24:118–35. doi: 10.1016/j.plrev.2017.11.003
66. Gosak M, Stozer A, Markovic R, Dolensek J, Marhl M, Rupnik MS, et al. The relationship between node degree and dissipation rate in networks of diffusively coupled oscillators and its significance for pancreatic beta cells. Chaos (Woodbury N.Y.) (2015) 25:073115. doi: 10.1063/1.4926673
67. Johnston NR, Mitchell RK, Haythorne E, Pessoa MP, Semplici F, Ferrer J, et al. Beta cell hubs dictate pancreatic islet responses to glucose. Cell Metab (2016) 24:389–401. doi: 10.1016/j.cmet.2016.06.020
68. Rutter GA, Hodson DJ, Chabosseau P, Haythorne E, Pullen TJ, Leclerc I. Local and regional control of calcium dynamics in the pancreatic islet. Diabetes Obes Metab (2017) 19(Suppl 1):30–41. doi: 10.1111/dom.12990
69. Sterk M, Dolensek J, Skelin Klemen M, Krizancic Bombek L, Paradiz Leitgeb E, Kercmar J, et al. Functional characteristics of hub and wave-initiator cells in beta cell networks. Biophys J (2023) 122(5):784–801. doi: 10.1016/j.bpj.2023.01.039
70. Mears D, Sheppard NF, Atwater I, Rojas E. Magnitude and modulation of pancreatic beta-cell gap junction electrical conductance in-situ. J Membrane Biol (1995) 146:163–76. doi: 10.1007/BF00238006
71. Allagnat F, Martin D, Condorelli DF, Waeber G, Haefliger J-A. Glucose represses connexin36 in insulin-secreting cells. J Cell Sci (2005) 118:5335–44. doi: 10.1242/jcs.02600
72. Hodson DJ, Mitchell RK, Bellomo EA, Sun G, Vinet L, Meda P, et al. Lipotoxicity disrupts incretin-regulated human beta cell connectivity. J Clin Invest (2013) 123:4182–94. doi: 10.1172/JCI68459
73. Haefliger J-A, Martin D, Favre D, Petremand Y, Mazzolai L, Abderrahmani A, et al. Reduction of connexin36 content by ICER-1 contributes to insulin-secreting cells apoptosis induced by oxidized LDL particles. PloS One (2013) 8:e55198. doi: 10.1371/journal.pone.0055198
74. Farnsworth NL, Hemmati A, Pozzoli M, Benninger RK. Fluorescence recovery after photobleaching reveals regulation and distribution of connexin36 gap junction coupling within mouse islets of Langerhans. J Physiol (2014) 592:4431–46. doi: 10.1113/jphysiol.2014.276733
75. Farnsworth NL, Benninger RK. New insights into the role of connexins in pancreatic islet function and diabetes. FEBS Lett (2014) 588:1278–87. doi: 10.1016/j.febslet.2014.02.035
76. Kothmann WW, Massey SC, O’Brien J. Dopamine-stimulated dephosphorylation of connexin 36 mediates AII amacrine cell uncoupling. J Neurosci (2009) 29:14903–11. doi: 10.1523/JNEUROSCI.3436-09.2009
77. Lee T-M, Lin S-Z, Chang N-C. Both PKA and Epac pathways mediate N-acetylcysteine-induced Connexin43 preservation in rats with myocardial infarction. PloS One (2013) 8:e71878. doi: 10.1371/journal.pone.0071878
78. Somekawa S, Fukuhara S, Nakaoka Y, Fujita H, Saito Y, Mochizuki N. Enhanced functional gap junction neoformation by protein kinase A-dependent and Epac-dependent signals downstream of cAMP in cardiac myocytes. Circ Res (2005) 97:655–62. doi: 10.1161/01.RES.0000183880.49270.f9
79. Stožer A, Paradiž Leitgeb E, Pohorec V, Dolenšek J, Križančić Bombek L, Gosak M, et al. The role of cAMP in beta cell stimulus-secretion and intercellular coupling. Cells (2021) 10:1658. doi: 10.3390/cells10071658
80. Speier S, Rupnik M. A novel approach to in situ characterization of pancreatic ß-cells. Pflügers Archiv Eur J Physiol (2003) 446:553–8. doi: 10.1007/s00424-003-1097-9
81. Stožer A, Skelin Klemen M, Gosak M, Križančić Bombek L, Pohorec V, Slak Rupnik M, et al. Glucose-dependent activation, activity, and deactivation of beta cell networks in acute mouse pancreas tissue slices. Am J Physiol-Endocrinol Metab (2021) 321:E305–23. doi: 10.1152/ajpendo.00043.2021
82. Nadal A, Quesada I, Soria B. Homologous and heterologous asynchronicity between identified α-, β- and δ-cells within intact islets of Langerhans in the mouse. J Physiol (1999) 517:85–93. doi: 10.1111/j.1469-7793.1999.0085z.x
83. Asada N, Shibuya I, Iwanaga T, Niwa K, Kanno T. Identification of alpha- and beta-cells in intact isolated islets of Langerhans by their characteristic cytoplasmic Ca2+ concentration dynamics and immunocytochemical staining. Diabetes (1998) 47:751–7. doi: 10.2337/diabetes.47.5.751
84. Quesada I, Nadal A, Soria B. Different effects of tolbutamide and diazoxide in alpha, beta-, and delta-cells within intact islets of Langerhans. Diabetes (1999) 48:2390–7. doi: 10.2337/diabetes.48.12.2390
85. Shuai H, Xu Y, Yu Q, Gylfe E, Tengholm A. Fluorescent protein vectors for pancreatic islet cell identification in live-cell imaging. Pflugers Archiv: Eur J Physiol (2016) 468:1765–77. doi: 10.1007/s00424-016-1864-z
86. Stožer A, Dolenšek J, Križančić Bombek L, Pohorec V, Slak Rupnik M, Klemen MS. Confocal laser scanning microscopy of calcium dynamics in acute mouse pancreatic tissue slices. JoVE (2021):e62293. doi: 10.3791/62293
87. Stožer A, Dolenšek J, Križančić Bombek L, Gosak M, Skelin Klemen M. Calcium imaging, Imaging Modalities for Biological and Preclinical Research: A Compendium Vol. 1. Bristol, UK: IOP Publishing (2021) p. I.1.e–1-I.1.e-13.
88. Pohorec V, Krizancic Bombek L, Skelin Klemen M, Dolensek J, Stozer A. Glucose-stimulated calcium dynamics in beta cells from male C57BL/6J, C57BL/6N, and NMRI mice: A comparison of activation, activity, and deactivation properties in tissue slices. Front Endocrinol (2022) 13:867663. doi: 10.3389/fendo.2022.867663
89. Song SE, Shin SK, Park SY, Hwang IS, Im SS, Bae JH, et al. Epac2a-knockout mice are resistant to dexamethasone-induced skeletal muscle atrophy and short-term cold stress. BMB Rep (2018) 51:39–44. doi: 10.5483/BMBRep.2018.51.1.132
90. MacDonald PE, Salapatek AM, Wheeler MB. Glucagon-like peptide-1 receptor activation antagonizes voltage-dependent repolarizing K(+) currents in beta-cells: a possible glucose-dependent insulinotropic mechanism. Diabetes (2002) 51(Suppl 3):S443–7. doi: 10.2337/diabetes.51.2007.S443
91. Hodson DJ, Tarasov AI, Gimeno Brias S, Mitchell RK, Johnston NR, Haghollahi S, et al. Incretin-modulated beta cell energetics in intact islets of Langerhans. Mol Endocrinol (2014) 28(6):me.2014–1038. doi: 10.1210/me.2014-1038
92. Gromada J, Bokvist K, Ding WG, Holst JJ, Nielsen JH, Rorsman P. Glucagon-like peptide 1 (7-36) amide stimulates exocytosis in human pancreatic beta-cells by both proximal and distal regulatory steps in stimulus-secretion coupling. Diabetes (1998) 47:57–65. doi: 10.2337/diab.47.1.57
93. Irles E, Neco P, Lluesma M, Villar-Pazos S, Santos-Silva JC, Vettorazzi JF, et al. Enhanced glucose-induced intracellular signaling promotes insulin hypersecretion: pancreatic beta-cell functional adaptations in a model of genetic obesity and prediabetes. Mol Cell Endocrinol (2015) 404:46–55. doi: 10.1016/j.mce.2015.01.033
94. Sluga N, Krizancic Bombek L, Kercmar J, Sarikas S, Postic S, Pfabe J, et al. Physiological levels of adrenaline fail to stop pancreatic beta cell activity at unphysiologically high glucose levels. Front Endocrinol (2022) 13:1013697. doi: 10.3389/fendo.2022.1013697
95. Britsch S, Krippeit-Drews P, Lang F, Gregor M, Drews G. Glucagon-like peptide-1 modulates Ca2+ current but not K+ATP current in intact mouse pancreatic B-cells. Biochem Biophys Res Commun (1995) 207:33–9. doi: 10.1006/bbrc.1995.1149
96. Do OH, Low JT, Gaisano HY, Thorn P. The secretory deficit in islets from db/db mice is mainly due to a loss of responding beta cells. Diabetologia (2014) 57:1400–9. doi: 10.1007/s00125-014-3226-8
97. Gilon P, Henquin JC. Influence of membrane potential changes on cytoplasmic Ca2+ concentration in an electrically excitable cell, the insulin-secreting pancreatic B-cell. J Biol Chem (1992) 267:20713–20. doi: 10.1016/S0021-9258(19)36744-4
98. MacDonald PE, Rorsman P. Oscillations, intercellular coupling, and insulin secretion in pancreatic β Cells. PloS Biol (2006) 4:e49. doi: 10.1371/journal.pbio.0040049
99. Jonkers FC, Henquin J-C. Measurements of cytoplasmic ca2+ in islet cell clusters show that glucose rapidly recruits β-cells and gradually increases the individual cell response. Diabetes (2001) 50:540–50. doi: 10.2337/diabetes.50.3.540
100. Gonzalez A, Merino B, Marroqui L, Neco P, Alonso-Magdalena P, Caballero-Garrido E, et al. Insulin hypersecretion in islets from diet-induced hyperinsulinemic obese female mice is associated with several functional adaptations in individual beta-cells. Endocrinology (2013) 154:3515–24. doi: 10.1210/en.2013-1424
101. Bertuzzi F, Davalli AM, Nano R, Socci C, Codazzi F, Fesce R, et al. Mechanisms of coordination of Ca2+ signals in pancreatic islet cells. Diabetes (1999) 48:1971–8. doi: 10.2337/diabetes.48.10.1971
102. Speier S, Gjinovci A, Charollais A, Meda P, Rupnik M. Cx36-mediated coupling reduces β-cell heterogeneity, confines the stimulating glucose concentration range, and affects insulin release kinetics. Diabetes (2007) 56:1078–86. doi: 10.2337/db06-0232
103. Satin LS, Zhang Q, Rorsman P. “Take me to your leader”: an electrophysiological appraisal of the role of hub cells in pancreatic islets. Diabetes (2020) 69:830–6. doi: 10.2337/dbi19-0012
104. Santos RM, Rosario LM, Nadal A, Garcia-Sancho J, Soria B, Valdeolmillos M. Widespread synchronous Ca oscillations due to bursting electrical activity in single pancreatic islets. Pflügers Archiv Eur J Physiol (1991) 418:417–22. doi: 10.1007/BF00550880
105. Bergsten P. Role of oscillations in membrane potential, cytoplasmic ca2+, and metabolism for plasma insulin oscillations. Diabetes (2002) 51:S171–6. doi: 10.2337/diabetes.51.2007.S171
106. Postic S, Sarikas S, Pfabe J, Pohorec V, Krizancic Bombek L, Sluga N, et al. High-resolution analysis of the cytosolic Ca(2+) events in beta cell collectives in situ. Am J Physiol Endocrinol Metab (2023) 324:E42–55. doi: 10.1152/ajpendo.00165.2022
107. Gromada J, Ding W-G, Barg S, Renström E, Rorsman P. Multisite regulation of insulin secretion by cAMP-increasing agonists: evidence that glucagon-like peptide 1 and glucagon act via distinct receptors. Pflügers Archiv Eur J Physiol (1997) 434:515–24. doi: 10.1007/s004240050431
108. Ikeuchi M, Cook DL. Glucagon and forskolin have dual effects upon islet cell electrical activity. Life Sci (1984) 35:685–91. doi: 10.1016/0024-3205(84)90264-9
109. Henquin JC, Schmeer W, Meissner HP. Forskolin, an activator of adenylate cyclase, increases CA2+-dependent electrical activity induced by glucose in mouse pancreatic B cells. Endocrinology (1983) 112:2218–20. doi: 10.1210/endo-112-6-2218
110. Miki T, Minami K, Shinozaki H, Matsumura K, Saraya A, Ikeda H, et al. Distinct effects of glucose-dependent insulinotropic polypeptide and glucagon-like peptide-1 on insulin secretion and gut motility. Diabetes (2005) 54:1056–63. doi: 10.2337/diabetes.54.4.1056
111. Dzhura I, Chepurny OG, Kelley GG, Leech CA, Roe MW, Dzhura E, et al. Epac2-dependent mobilization of intracellular Ca2+ by glucagon-like peptide-1 receptor agonist exendin-4 is disrupted in β-cells of phospholipase C-ϵ knockout mice. J Physiol (2010) 588:4871–89. doi: 10.1113/jphysiol.2010.198424
112. Meda P. The in vivo beta-to-beta-cell chat room: connexin connections matter. Diabetes (2012) 61:1656–8. doi: 10.2337/db12-0336
113. Ravier MA, Güldenagel M, Charollais A, Gjinovci A, Caille D, Söhl G, et al. Loss of connexin36 channels alters β-cell coupling, islet synchronization of glucose-induced ca2+ and insulin oscillations, and basal insulin release. Diabetes (2005) 54:1798–807. doi: 10.2337/diabetes.54.6.1798
114. Šterk M, Križančić Bombek L, Skelin Klemen M, Slak Rupnik M, Marhl M, Stožer A, et al. NMDA receptor inhibition increases, synchronizes, and stabilizes the collective pancreatic beta cell activity: Insights through multilayer network analysis. PloS Comput Biol (2021) 17:e1009002. doi: 10.1371/journal.pcbi.1009002
115. MIranda JG, Schleicher WE, Wells KL, Ramirez DG, Landgrave SP, Benninger RKP. Dynamic changes in β-cell [Ca(2+)] regulate NFAT activation, gene transcription, and islet gap junction communication. Mol Metab (2022) 57:101430. doi: 10.1016/j.molmet.2021.101430
116. Segretain D, Falk MM. Regulation of connexin biosynthesis, assembly, gap junction formation, and removal. Biochim Biophys Acta (2004) 1662:3–21. doi: 10.1016/j.bbamem.2004.01.007
117. Laird DW. Life cycle of connexins in health and disease. Biochem J (2006) 394:527–43. doi: 10.1042/BJ20051922
118. Urschel S, Hoher T, Schubert T, Alev C, Sohl G, Worsdorfer P, et al. Protein kinase A-mediated phosphorylation of connexin36 in mouse retina results in decreased gap junctional communication between AII amacrine cells. J Biol Chem (2006) 281:33163–71. doi: 10.1074/jbc.M606396200
119. Ivanova E, Yee CW, Sagdullaev BT. Increased phosphorylation of Cx36 gap junctions in the AII amacrine cells of RD retina. Front Cell Neurosci (2015) 9:390. doi: 10.3389/fncel.2015.00390
120. Pidoux G, Tasken K. Anchored PKA as a gatekeeper for gap junctions. Commun Integr Biol (2015) 8:e1057361. doi: 10.1080/19420889.2015.1057361
121. Li X, Lynn BD, Nagy JI. The effector and scaffolding proteins AF6 and MUPP1 interact with connexin36 and localize at gap junctions that form electrical synapses in rodent brain. Eur J Neurosci (2012) 35:166–81. doi: 10.1111/j.1460-9568.2011.07947.x
122. Ishii M, Mueller I Fau - Nakajima T, Nakajima T Fau - Pasquale EB, Pasquale Eb Fau - Ogawa K, Ogawa K. EphB signaling inhibits gap junctional intercellular communication and synchronized contraction in cultured cardiomyocytes. Basic Res Cardiol (2011) 106(6):1057–68. doi: 10.1007/s00395-011-0219-3
123. Rupnik M. All together now: exocytose or fail. Islets (2009) 1:78–80. doi: 10.4161/isl.1.1.8757
124. Zhang F, Tzanakakis ES. Optogenetic regulation of insulin secretion in pancreatic beta-cells. Sci Rep (2017) 7:9357. doi: 10.1038/s41598-017-09937-0
125. McQuaid TS, Saleh MC, Joseph JW, Gyulkhandanyan A, Manning-Fox JE, MacLellan JD, et al. cAMP-mediated signaling norMalizes glucose-stimulated insulin secretion in uncoupling protein-2 overexpressing beta-cells. J Endocrinol (2006) 190:669–80. doi: 10.1677/joe.1.06723
126. Tian G, Sandler S, Gylfe E, Tengholm A. Glucose- and hormone-induced cAMP oscillations in α- and β-cells within intact pancreatic islets. Diabetes (2011) 60:1535–43. doi: 10.2337/db10-1087
127. Qadir MMF, Álvarez-Cubela S, Weitz J, Panzer JK, Klein D, Moreno-Hernández Y, et al. Long-term culture of human pancreatic slices as a model to study real-time islet regeneration. Nat Commun (2020) 11:3265–5. doi: 10.1038/s41467-020-17040-8
Keywords: pancreas, tissue slices, beta cells, calcium imaging, amplifying pathway, forskolin, Epac2A KO, intercellular network
Citation: Skelin Klemen M, Dolenšek J, Križančić Bombek L, Pohorec V, Gosak M, Slak Rupnik M and Stožer A (2023) The effect of forskolin and the role of Epac2A during activation, activity, and deactivation of beta cell networks. Front. Endocrinol. 14:1225486. doi: 10.3389/fendo.2023.1225486
Received: 19 May 2023; Accepted: 08 August 2023;
Published: 28 August 2023.
Edited by:
Marcia Hiriart, Universidad Nacional Autonoma de Mexico, MexicoReviewed by:
David Hodson, University of Oxford, United KingdomKazuki Harada, The University of Tokyo, Japan
Copyright © 2023 Skelin Klemen, Dolenšek, Križančić Bombek, Pohorec, Gosak, Slak Rupnik and Stožer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Andraž Stožer, YW5kcmF6LnN0b3plckB1bS5zaQ==