Search

Beta cells within the pancreatic islets of Langerhans respond to stimulation with coherent oscillations of membrane potential and intracellular calcium concentration that presumably drive the pulsatile exocytosis of insulin. Their rhythmic activity is multimodal, resulting from networked feedback interactions of various oscillatory subsystems, such as the glycolytic, mitochondrial, and electrical/calcium components. How these oscillatory modules interact and affect the collective cellular activity, which is a prerequisite for proper hormone release, is incompletely understood. In the present work, we combined advanced confocal Ca 2+ imaging in fresh mouse pancreas tissue slices with time series analysis and network science approaches to unveil the glucose-dependent characteristics of different oscillatory components on both the intra- and inter-cellular level. Our results reveal an interrelationship between the metabolically driven low-frequency component and the electrically driven high-frequency component, with the latter exhibiting the highest bursting rates around the peaks of the slow component and the lowest around the nadirs. Moreover, the activity, as well as the average synchronicity of the fast component, considerably increased with increasing stimulatory glucose concentration, whereas the stimulation level did not affect any of these parameters in the slow component domain. Remarkably, in both dynamical components, the average correlation decreased similarly with intercellular distance, which implies that intercellular communication affects the synchronicity of both types of oscillations. To explore the intra-islet synchronization patterns in more detail, we constructed functional connectivity maps. The subsequent comparison of network characteristics of different oscillatory components showed more locally clustered and segregated networks of fast oscillatory activity, while the slow oscillations were more global, resulting in several long-range connections and a more cohesive structure. Besides the structural differences, we found a relatively weak relationship between the fast and slow network layer, which suggests that different synchronization mechanisms shape the collective cellular activity in islets, a finding which has to be kept in mind in future studies employing different oscillations for constructing networks., Competing Interests: 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., (Copyright © 2021 Zmazek, Klemen, Markovič, Dolenšek, Marhl, Stožer and Gosak.)

Introduction: The tissue slice technique offers several benefits compared to isolated cells and cell clusters that help us understand the (patho)physiology of several organs in situ . The most prominent features are preserved architecture and function, with intact homotypic and heterotypic interactions between cells in slices. In the pancreas, this technique has been utilized successfully to study acinar and endocrine islet cells. However, it has never been used to investigate ductal function. Since pancreatic ductal epithelial cells (PDECs) play an essential role in the physiology of the pancreas, our aim was to use this technique to study PDEC structure and function in situ . Materials and methods: Eight- to sixteen weeks old C57BL/6 mice were used for preparation of pancreas tissue slices. Low melting point agarose was injected into the common bile duct and the whole organ was extracted. For morphological studies, pieces of tissue were embedded in agarose and cryosectioned to obtain 15 μm thick slices. In order to visualize pancreatic ducts, (i) the Giemsa dye was added to the agarose and visualized using light microscopy or (ii) immunostaining for the cystic fibrosis transmembrane conductance regulator (CFTR) was performed. For functional characterization, agarose-embedded tissue was immediately cut to 140 μm thick tissue slices that were loaded with the cell permeant form of the Oregon Green 488 BAPTA-1 dye and used for confocal calcium imaging. Results: Giemsa staining has shown that the injected agarose reaches the head and body of the pancreas to a greater extent than the tail, without disrupting the tissue architecture. Strong CFTR expression was detected at the apical membranes of PDECs and acinar cells, whereas islet cells were completely negative for CFTR. Stimulation with chenodeoxycholic acid (CDCA, 1 mM) resulted in a robust transient increase in intracellular calcium concentration that was readily visible in >40 ductal cells per slice. Conclusion: Our results confirm that the acutely-isolated pancreas tissue slice technique is suitable for structural and functional investigation of PDECs and their relationship with other cell types, such as acini and endocrine cells in situ . In combination with different genetic, pharmacological or dietary approaches it could become a method of choice in the foreseeable future.

Self-organized critical dynamics is assumed to be an attractive mode of functioning for several real-life systems and entails an emergent activity in which the extent of observables follows a power-law distribution. The hallmarks of criticality have recently been observed in a plethora of biological systems, including beta cell populations within pancreatic islets of Langerhans. In the present study, we systematically explored the mechanisms that drive the critical and supercritical behavior in networks of coupled beta cells under different circumstances by means of experimental and computational approaches. Experimentally, we employed high-speed functional multicellular calcium imaging of fluorescently labeled acute mouse pancreas tissue slices to record calcium signals in a large number of beta cells simultaneously, and with a high spatiotemporal resolution. Our experimental results revealed that the cellular responses to stimulation with glucose are biphasic and glucose-dependent. Under physiological as well as under supraphysiological levels of stimulation, an initial activation phase was followed by a supercritical plateau phase with a high number of global intercellular calcium waves. However, the activation phase displayed fingerprints of critical behavior under lower stimulation levels, with a progressive recruitment of cells and a power-law distribution of calcium wave sizes. On the other hand, the activation phase provoked by pathophysiologically high glucose concentrations, differed considerably and was more rapid, less continuous, and supercritical. To gain a deeper insight into the experimentally observed complex dynamical patterns, we built up a phenomenological model of coupled excitable cells and explored empirically the model's necessities that ensured a good overlap between computational and experimental results. It turned out that such a good agreement between experimental and computational findings was attained when both heterogeneous and stimulus-dependent time lags, variability in excitability levels, as well as a heterogeneous cell-cell coupling were included into the model. Most importantly, since our phenomenological approach involved only a few parameters, it naturally lends itself not only for determining key mechanisms of self-organized criticality at the tissue level, but also points out various features for comprehensive and realistic modeling of different excitable systems in nature.

A coordinated functioning of beta cells within pancreatic islets is mediated by oscillatory membrane depolarization and subsequent changes in cytoplasmic calcium concentration. While gap junctions allow for intraislet information exchange, beta cells within islets form complex syncytia that are intrinsically nonlinear and highly heterogeneous. To study spatiotemporal calcium dynamics within these syncytia, we make use of computational modeling and confocal high-speed functional multicellular imaging. We show that model predictions are in good agreement with experimental data, especially if a high degree of heterogeneity in the intercellular coupling term is assumed. In particular, during the first few minutes after stimulation, the probability distribution of calcium wave sizes is characterized by a power law, thus indicating critical behavior. After this period, the dynamics changes qualitatively such that the number of global intercellular calcium events increases to the point where the behavior becomes supercritical. To better mimic normal in vivo conditions, we compare the described behavior during supraphysiological non-oscillatory stimulation with the behavior during exposure to a slightly lower and oscillatory glucose challenge. In the case of this protocol, we observe only critical behavior in both experiment and model. Our results indicate that the loss of oscillatory changes, along with the rise in plasma glucose observed in diabetes, could be associated with a switch to supercritical calcium dynamics and loss of beta cell functionality.