Search

Ca 2+ ions play a central role in the stimulus-secretion coupling cascade of pancreatic beta cells. The use of confocal microscopy in conjunction with the acute pancreas tissue slice technique offers valuable insights into changes in the intracellular calcium concentration following stimulation by secretagogues. This allows the study of beta cells on a single cell level, as well as their behavior on a multicellular scale within an intact environment. With the use of advanced analytical tools, this approach offers insight into how single cells contribute to the functional unit of islets of Langerhans and processes underlying insulin secretion. Here we describe a comprehensive protocol for the preparation and utilization of acute pancreas tissue slices in mice, the use of high-resolution confocal microscopy for observation of glucose-stimulated calcium dynamics in beta cells, and the computational analysis for objective evaluation of calcium signals., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

We propose a detailed computational beta cell model that emphasizes the role of anaplerotic metabolism under glucose and glucose-glutamine stimulation. This model goes beyond the traditional focus on mitochondrial oxidative phosphorylation and ATP-sensitive K + channels, highlighting the predominant generation of ATP from phosphoenolpyruvate in the vicinity of K ATP channels. It also underlines the modulatory role of H 2 O 2 as a signaling molecule in the first phase of glucose-stimulated insulin secretion. In the second phase, the model emphasizes the critical role of anaplerotic pathways, activated by glucose stimulation via pyruvate carboxylase and by glutamine via glutamate dehydrogenase. It particularly focuses on the production of NADPH and glutamate as key enhancers of insulin secretion. The predictions of the model are consistent with empirical data, highlighting the complex interplay of metabolic pathways and emphasizing the primary role of glucose and the facilitating role of glutamine in insulin secretion. By delineating these crucial metabolic pathways, the model provides valuable insights into potential therapeutic targets for diabetes., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Within the islets of Langerhans, beta cells orchestrate synchronized insulin secretion, a pivotal aspect of metabolic homeostasis. Despite the inherent heterogeneity and multimodal activity of individual cells, intercellular coupling acts as a homogenizing force, enabling coordinated responses through the propagation of intercellular waves. Disruptions in this coordination are implicated in irregular insulin secretion, a hallmark of diabetes. Recently, innovative approaches, such as integrating multicellular calcium imaging with network analysis, have emerged for a quantitative assessment of the cellular activity in islets. However, different groups use distinct experimental preparations, microscopic techniques, apply different methods to process the measured signals and use various methods to derive functional connectivity patterns. This makes comparisons between findings and their integration into a bigger picture difficult and has led to disputes in functional connectivity interpretations. To address these issues, we present here a systematic analysis of how different approaches influence the network representation of islet activity. Our findings show that the choice of methods used to construct networks is not crucial, although care is needed when combining data from different islets. Conversely, the conclusions drawn from network analysis can be heavily affected by the pre-processing of the time series, the type of the oscillatory component in the signals, and by the experimental preparation. Our tutorial-like investigation aims to resolve interpretational issues, reconcile conflicting views, advance functional implications, and encourage researchers to adopt connectivity analysis. As we conclude, we outline challenges for future research, emphasizing the broader applicability of our conclusions to other tissues exhibiting complex multicellular dynamics., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Šterk et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Tight control of beta cell stimulus-secretion coupling is crucial for maintaining homeostasis of energy-rich nutrients. While glucose serves as a primary regulator of this process, incretins augment beta cell function, partly by enhancing cytosolic [Ca 2+ ] dynamics. However, the details of how precisely they affect beta cell recruitment during activation, their active time, and functional connectivity during plateau activity, and how they influence beta cell deactivation remain to be described. Performing functional multicellular Ca 2+ imaging in acute mouse pancreas tissue slices enabled us to systematically assess the effects of the GLP-1 receptor agonist exendin-4 (Ex-4) simultaneously in many coupled beta cells with high resolution. In otherwise substimulatory glucose, Ex-4 was able to recruit approximately a quarter of beta cells into an active state. Costimulation with Ex-4 and stimulatory glucose shortened the activation delays and accelerated beta cell activation dynamics. More specifically, active time increased faster, and the time required to reach half-maximal activation was effectively halved in the presence of Ex-4. Moreover, the active time and regularity of [Ca 2+ ] IC oscillations increased, especially during the first part of beta cell response. In contrast, subsequent addition of Ex-4 to already active cells did not significantly enhance beta cell activity. Network analyses further confirmed increased connectivity during activation and activity in the presence of Ex-4, with hub cell roles remaining rather stable in both control experiments and experiments with Ex-4. Interestingly, Ex-4 demonstrated a biphasic effect on deactivation, slightly prolonging beta cell activity at physiological concentrations and shortening deactivation delays at supraphysiological concentrations. In sum, costimulation by Ex-4 and glucose increases [Ca 2+ ] IC during beta cell activation and activity, indicating that the effect of incretins may, to an important extent, be explained by enhanced [Ca 2+ ] IC signals. During deactivation, previous incretin stimulation does not critically prolong cellular activity, which corroborates their low risk of hypoglycemia., 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. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Paradiž Leitgeb, Kerčmar, Križančić Bombek, Pohorec, Skelin Klemen, Slak Rupnik, Gosak, Dolenšek and Stožer.)

Pancreatic beta cells are coupled excitable oscillators that synchronize their activity via different communication pathways. Their oscillatory activity manifests itself on multiple timescales and consists of bursting electrical activity, subsequent oscillations in the intracellular Ca^{2+}, as well as oscillations in metabolism and exocytosis. The coordination of the intricate activity on the multicellular level plays a key role in the regulation of physiological pulsatile insulin secretion and is incompletely understood. In this paper, we investigate theoretically the principles that give rise to the synchronized activity of beta cell populations by building up a phenomenological multicellular model that incorporates the basic features of beta cell dynamics. Specifically, the model is composed of coupled slow and fast oscillatory units that reflect metabolic processes and electrical activity, respectively. Using a realistic description of the intercellular interactions, we study how the combination of electrical and metabolic coupling generates collective rhythmicity and shapes functional beta cell networks. It turns out that while electrical coupling solely can synchronize the responses, the addition of metabolic interactions further enhances coordination, the spatial range of interactions increases the number of connections in the functional beta cell networks, and ensures a better consistency with experimental findings. Moreover, our computational results provide additional insights into the relationship between beta cell heterogeneity, their activity profiles, and functional connectivity, supplementing thereby recent experimental results on endocrine networks.

Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

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., 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 © 2023 Skelin Klemen, Dolenšek, Križančić Bombek, Pohorec, Gosak, Slak Rupnik and Stožer.)

We study collective failures in biologically realistic networks that consist of coupled excitable units. The networks have broad-scale degree distribution, high modularity, and small-world properties, while the excitable dynamics is determined by the paradigmatic FitzHugh-Nagumo model. We consider different coupling strengths, bifurcation distances, and various aging scenarios as potential culprits of collective failure. We find that for intermediate coupling strengths, the network remains globally active the longest if the high-degree nodes are first targets for inactivation. This agrees well with previously published results, which showed that oscillatory networks can be highly fragile to the targeted inactivation of low-degree nodes, especially under weak coupling. However, we also show that the most efficient strategy to enact collective failure does not only non-monotonically depend on the coupling strength, but it also depends on the distance from the bifurcation point to the oscillatory behavior of individual excitable units. Altogether, we provide a comprehensive account of determinants of collective failure in excitable networks, and we hope this will prove useful for better understanding breakdowns in systems that are subject to such dynamics., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

Islets of Langerhans operate as multicellular networks in which several hundred β cells work in synchrony to produce secretory pulses of insulin, a hormone crucial for controlling metabolic homeostasis. Their collective rhythmic activity is facilitated by gap junctional coupling and affected by their functional heterogeneity, but the details of this robust and coordinated behavior are still not fully understood. Recent advances in multicellular imaging and optogenetic and photopharmacological strategies, as well as in network science, have led to the discovery of specialized β cell subpopulations that were suggested to critically determine the collective dynamics in the islets. In particular hubs, i.e., β cells with many functional connections, are believed to significantly enhance communication capacities of the intercellular network and facilitate an efficient spreading of intercellular Ca 2+ waves, whereas wave-initiator cells trigger intercellular signals in their cohorts. Here, we determined Ca 2+ signaling characteristics of these two β cell subpopulations and the relationship between them by means of functional multicellular Ca 2+ imaging in mouse pancreatic tissue slices in combination with methods of complex network theory. We constructed network layers based on individual Ca 2+ waves to identify wave initiators, and functional correlation-based networks to detect hubs. We found that both cell types exhibit a higher-than-average active time under both physiological and supraphysiological glucose concentrations, but also that they differ significantly in many other functional characteristics. Specifically, Ca 2+ oscillations in hubs are more regular, and their role appears to be much more stable over time than for initiator cells. Moreover, in contrast to wave initiators, hubs transmit intercellular signals faster than other cells, which implies a stronger intercellular coupling. Our research indicates that hubs and wave-initiator cell subpopulations are both natural features of healthy pancreatic islets, but their functional roles in principle do not overlap and should thus not be considered equal., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Pancreatic islets are highly interconnected structures that produce pulses of insulin and other hormones, maintaining normal homeostasis of glucose and other nutrients. Normal stimulus-secretion and intercellular coupling are essential to regulated secretory responses, and these hallmarks are known to be altered in diabetes. In the current study, we used calcium imaging of isolated human islets to assess their collective behavior. The activity occurred in the form of calcium oscillations, was synchronized across different regions of islets through calcium waves, and was glucose dependent: higher glucose enhanced the activity, elicited a greater proportion of global calcium waves, and led to denser and less fragmented functional networks. Hub regions were identified in stimulatory conditions, and they were characterized by long active times. Moreover, calcium waves were found to be initiated in different subregions and the roles of initiators and hubs did not overlap. In type 2 diabetes, glucose dependence was retained, but reduced activity, locally restricted waves, and more segregated networks were detected compared with control islets. Interestingly, hub regions seemed to suffer the most by losing a disproportionately large fraction of connections. These changes affected islets from donors with diabetes in a heterogeneous manner., (© 2022 by the American Diabetes Association.)

Our core body temperature is held around [Formula: see text]C by an effective internal thermoregulatory system. However, various clinical scenarios have a more favorable outcome under external temperature regulation. Therapeutic hypothermia, for example, was found beneficial for the outcome of resuscitated cardiac arrest patients due to its protection against cerebral ischemia. Nonetheless, practice shows that outcomes of targeted temperature management vary considerably in dependence on individual tissue damage levels and differences in therapeutic strategies and protocols. Here, we address these differences in detail by means of computational modeling. We develop a multi-segment and multi-node thermoregulatory model that takes into account details related to specific post-cardiac arrest-related conditions, such as thermal imbalances due to sedation and anesthesia, increased metabolic rates induced by inflammatory processes, and various external cooling techniques. In our simulations, we track the evolution of the body temperature in patients subjected to post-resuscitation care, with particular emphasis on temperature regulation via an esophageal heat transfer device, on the examination of the alternative gastric cooling with ice slurry, and on how anesthesia and the level of inflammatory response influence thermal behavior. Our research provides a better understanding of the heat transfer processes and therapies used in post-cardiac arrest patients., (© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Technological advances in imaging techniques and biometric data acquisition have enabled us to apply methods of network science to study the morphology and structural design of organelles, organs, and tissues, as well as the coordinated interactions among them that yield a healthy physiology at the level of whole organisms. We here review research dedicated to these advances, in particular focusing on networks between cells, the topology of multicellular structures, neural interactions, fluid transportation networks, and anatomical networks. The percolation of blood vessels, structural connectivity within the brain, the porous structure of bones, and relations between different anatomical parts of the human body are just some of the examples that we explore in detail. We argue and show that the models, methods, and algorithms developed in the realm of network science are ushering in a new era of network-based inquiry into the morphology and structural design of living systems in the broadest possible terms. We also emphasize that the need and applicability of this research is likely to increase significantly in the years to come due to the rapid progress made in the development of bioartificial substitutes and tissue engineering., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Author(s). Published by Elsevier B.V. All rights reserved.)

Extracellular pH has the potential to affect various aspects of the pancreatic beta cell function. To explain this effect, a number of mechanisms was proposed involving both extracellular and intracellular targets and pathways. Here, we focus on reassessing the influence of extracellular pH on glucose-dependent beta cell activation and collective activity in physiological conditions. To this end we employed mouse pancreatic tissue slices to perform high-temporally resolved functional imaging of cytosolic Ca 2+ oscillations. We investigated the effect of either physiological H + excess or depletion on the activation properties as well as on the collective activity of beta cell in an islet. Our results indicate that lowered pH invokes activation of a subset of beta cells in substimulatory glucose concentrations, enhances the average activity of beta cells, and alters the beta cell network properties in an islet. The enhanced average activity of beta cells was determined indirectly utilizing cytosolic Ca 2+ imaging, while direct measuring of insulin secretion confirmed that this enhanced activity is accompanied by a higher insulin release. Furthermore, reduced functional connectivity and higher functional segregation at lower pH, both signs of a reduced intercellular communication, do not necessary result in an impaired insulin release., 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 © 2022 Postić, Gosak, Tsai, Pfabe, Sarikas, Stožer, Korošak, Yang and Slak Rupnik.)

Islets of Langerhans are multicellular microorgans located in the pancreas that play a central role in whole-body energy homeostasis. Through secretion of insulin and other hormones they regulate postprandial storage and interprandial usage of energy-rich nutrients. In these clusters of hormone-secreting endocrine cells, intricate cell-cell communication is essential for proper function. Electrical coupling between the insulin-secreting beta cells through gap junctions composed of connexin36 is particularly important, as it provides the required, most important, basis for coordinated responses of the beta cell population. The increasing evidence that gap-junctional communication and its modulation are vital to well-regulated secretion of insulin has stimulated immense interest in how subpopulations of heterogeneous beta cells are functionally arranged throughout the islets and how they mediate intercellular signals. In the last decade, several novel techniques have been proposed to assess cooperation between cells in islets, including the prosperous combination of multicellular imaging and network science. In the present contribution, we review recent advances related to the application of complex network approaches to uncover the functional connectivity patterns among cells within the islets. We first provide an accessible introduction to the basic principles of network theory, enumerating the measures characterizing the intercellular interactions and quantifying the functional integration and segregation of a multicellular system. Then we describe methodological approaches to construct functional beta cell networks, point out possible pitfalls, and specify the functional implications of beta cell network examinations. We continue by highlighting the recent findings obtained through advanced multicellular imaging techniques supported by network-based analyses, giving special emphasis to the current developments in both mouse and human islets, as well as outlining challenges offered by the multilayer network formalism in exploring the collective activity of islet cell populations. Finally, we emphasize that the combination of these imaging techniques and network-based analyses does not only represent an innovative concept that can be used to describe and interpret the physiology of islets, but also provides fertile ground for delineating normal from pathological function and for quantifying the changes in islet communication networks associated with the development of diabetes mellitus., 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 © 2022 Stožer, Šterk, Paradiž Leitgeb, Markovič, Skelin Klemen, Ellis, Križančić Bombek, Dolenšek, MacDonald and Gosak.)

The physiology and pathophysiology of the exocrine pancreas are in close connection to changes in intra-cellular Ca2+ concentration. Most of our knowledge is based on in vitro experiments on acinar cells or acini enzymatically isolated from their surroundings, which can alter their structure, physiology, and limit our understanding. Due to these limitations, the acute pancreas tissue slice technique was introduced almost two decades ago as a complementary approach to assess the morphology and physiology of both the endocrine and exocrine pancreas in a more conserved in situ setting. In this study, we extend previous work to functional multicellular calcium imaging on acinar cells in tissue slices. The viability and morphological characteristics of acinar cells within the tissue slice were assessed using the LIVE/DEAD assay, transmission electron microscopy, and immunofluorescence imaging. The main aim of our study was to characterize the responses of acinar cells to stimulation with acetylcholine and compare them with responses to cerulein in pancreatic tissue slices, with special emphasis on inter-cellular and inter-acinar heterogeneity and coupling. To this end, calcium imaging was performed employing confocal microscopy during stimulation with a wide range of acetylcholine concentrations and selected concentrations of cerulein. We show that various calcium oscillation parameters depend monotonically on the stimulus concentration and that the activity is rather well synchronized within acini, but not between acini. The acute pancreas tissue slice represents a viable and reliable experimental approach for the evaluation of both intra- and inter-cellular signaling characteristics of acinar cell calcium dynamics. It can be utilized to assess many cells simultaneously with a high spatiotemporal resolution, thus providing an efficient and high-yield platform for future studies of normal acinar cell biology, pathophysiology, and screening pharmacological substances., Competing Interests: The authors have declared that no competing interests exist.

The firing patterns of neuronal populations often exhibit emergent collective oscillations, which can display substantial regularity even though the dynamics of individual elements is very stochastic. One of the many phenomena that is often studied in this context is coherence resonance, where additional noise leads to improved regularity of spiking activity in neurons. In this work, we investigate how the coherence resonance phenomenon manifests itself in populations of excitatory and inhibitory neurons. In our simulations, we use the coupled FitzHugh-Nagumo oscillators in the excitable regime and in the presence of neuronal noise. Formally, our model is based on the concept of a two-layered network, where one layer contains inhibitory neurons, the other excitatory neurons, and the interlayer connections represent heterotypic interactions. The neuronal activity is simulated in realistic coupling schemes in which neurons within each layer are connected with undirected connections, whereas neurons of different types are connected with directed interlayer connections. In this setting, we investigate how different neurophysiological determinants affect the coherence resonance. Specifically, we focus on the proportion of inhibitory neurons, the proportion of excitatory interlayer axons, and the architecture of interlayer connections between inhibitory and excitatory neurons. Our results reveal that the regularity of simulated neural activity can be increased by a stronger damping of the excitatory layer. This can be accomplished with a higher proportion of inhibitory neurons, a higher fraction of inhibitory interlayer axons, a stronger coupling between inhibitory axons, or by a heterogeneous configuration of interlayer connections. Our approach of modeling multilayered neuronal networks in combination with stochastic dynamics offers a novel perspective on how the neural architecture can affect neural information processing and provide possible applications in designing networks of artificial neural circuits to optimize their function via noise-induced phenomena., 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 © 2022 Ristič and Gosak.)

Many details of glucose-stimulated intracellular calcium changes in β cells during activation, activity, and deactivation, as well as their concentration-dependence, remain to be analyzed. Classical physiological experiments indicated that in islets, functional differences between individual cells are largely attenuated, but recent findings suggest considerable intercellular heterogeneity, with some cells possibly coordinating the collective responses. To address the above with an emphasis on heterogeneity and describing the relations between classical physiological and functional network properties, we performed functional multicellular calcium imaging in mouse pancreas tissue slices over a wide range of glucose concentrations. During activation, delays to activation of cells and any-cell-to-first-responder delays are shortened, and the sizes of simultaneously responding clusters increased with increasing glucose concentrations. Exactly the opposite characterized deactivation. The frequency of fast calcium oscillations during activity increased with increasing glucose up to 12 mM glucose concentration, beyond which oscillation duration became longer, resulting in a homogenous increase in active time. In terms of functional connectivity, islets progressed from a very segregated network to a single large functional unit with increasing glucose concentration. A comparison between classical physiological and network parameters revealed that the first-responders during activation had longer active times during plateau and the most active cells during the plateau tended to deactivate later. Cells with the most functional connections tended to activate sooner, have longer active times, and deactivate later. Our findings provide a common ground for recent differing views on β cell heterogeneity and an important baseline for future studies of stimulus-secretion and intercellular coupling. NEW & NOTEWORTHY We assessed concentration-dependence in coupled β cells, degree of functional heterogeneity, and uncovered possible specialized subpopulations during the different phases of the response to glucose at the level of many individual cells. To this aim, we combined acute mouse pancreas tissue slices with functional multicellular calcium imaging over a wide range from threshold (7 mM) and physiological (8 and 9 mM) to supraphysiological (12 and 16 mM) glucose concentrations, classical physiological, and advanced network analyses.

Pancreatic beta cells secrete insulin in response to stimulation with glucose and other nutrients, and impaired insulin secretion plays a central role in development of diabetes mellitus. Pharmacological management of diabetes includes various antidiabetic drugs, including incretins. The incretin hormones, glucagon-like peptide-1 and gastric inhibitory polypeptide, potentiate glucose-stimulated insulin secretion by binding to G protein-coupled receptors, resulting in stimulation of adenylate cyclase and production of the secondary messenger cAMP, which exerts its intracellular effects through activation of protein kinase A or the guanine nucleotide exchange protein 2A. The molecular mechanisms behind these two downstream signaling arms are still not fully elucidated and involve many steps in the stimulus-secretion coupling cascade, ranging from the proximal regulation of ion channel activity to the central Ca 2+ signal and the most distal exocytosis. In addition to modifying intracellular coupling, the effect of cAMP on insulin secretion could also be at least partly explained by the impact on intercellular coupling. In this review, we systematically describe the possible roles of cAMP at these intra- and inter-cellular signaling nodes, keeping in mind the relevance for the whole organism and translation to humans.

We propose and study an epidemiological model on a social network that takes into account heterogeneity of the population and different vaccination strategies. In particular, we study how the COVID-19 epidemics evolves and how it is contained by different vaccination scenarios by taking into account data showing that older people, as well as individuals with comorbidities and poor metabolic health, and people coming from economically depressed areas with lower quality of life in general, are more likely to develop severe COVID-19 symptoms, and quicker loss of immunity and are therefore more prone to reinfection. Our results reveal that the structure and the spatial arrangement of subpopulations are important epidemiological determinants. In a healthier society the disease spreads more rapidly but the consequences are less disastrous as in a society with more prevalent chronic comorbidities. If individuals with poor health are segregated within one community, the epidemic outcome is less favorable. Moreover, we show that, contrary to currently widely adopted vaccination policies, prioritizing elderly and other higher-risk groups is beneficial only if the supply of vaccine is high. If, however, the vaccination availability is limited, and if the demographic distribution across the social network is homogeneous, better epidemic outcomes are achieved if healthy people are vaccinated first. Only when higher-risk groups are segregated, like in elderly homes, their prioritization will lead to lower COVID-19 related deaths. Accordingly, young and healthy individuals should view vaccine uptake as not only protecting them, but perhaps even more so protecting the more vulnerable socio-demographic groups., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 The Author(s).)

NMDA receptors promote repolarization in pancreatic beta cells and thereby reduce glucose-stimulated insulin secretion. Therefore, NMDA receptors are a potential therapeutic target for diabetes. While the mechanism of NMDA receptor inhibition in beta cells is rather well understood at the molecular level, its possible effects on the collective cellular activity have not been addressed to date, even though proper insulin secretion patterns result from well-synchronized beta cell behavior. The latter is enabled by strong intercellular connectivity, which governs propagating calcium waves across the islets and makes the heterogeneous beta cell population work in synchrony. Since a disrupted collective activity is an important and possibly early contributor to impaired insulin secretion and glucose intolerance, it is of utmost importance to understand possible effects of NMDA receptor inhibition on beta cell functional connectivity. To address this issue, we combined confocal functional multicellular calcium imaging in mouse tissue slices with network science approaches. Our results revealed that NMDA receptor inhibition increases, synchronizes, and stabilizes beta cell activity without affecting the velocity or size of calcium waves. To explore intercellular interactions more precisely, we made use of the multilayer network formalism by regarding each calcium wave as an individual network layer, with weighted directed connections portraying the intercellular propagation. NMDA receptor inhibition stabilized both the role of wave initiators and the course of waves. The findings obtained with the experimental antagonist of NMDA receptors, MK-801, were additionally validated with dextrorphan, the active metabolite of the approved drug dextromethorphan, as well as with experiments on NMDA receptor KO mice. In sum, our results provide additional and new evidence for a possible role of NMDA receptor inhibition in treatment of type 2 diabetes and introduce the multilayer network paradigm as a general strategy to examine effects of drugs on connectivity in multicellular systems., Competing Interests: The authors have declared that no competing interests exist.

Ca 2+ homeostasis and signaling disturbances are associated with lens pathophysiology and are involved in cataract formation. Here, we explored the spatiotemporal changes in Ca 2+ signaling in lens epithelial cells (LECs) upon local mechanical stimulation, to better understand the LECs' intercellular communication and its association with cataractogenesis. We were interested in if the progression of the cataract affects the Ca 2+ signaling and if modifications of the Ca 2+ homeostasis in LECs are associated with different cataract types. Experiments were done on the human postoperative anterior lens capsule (LC) preparations consisting of the monolayer of LECs on the basement membrane. Our findings revealed that the Ca 2+ signal spreads radially from the stimulation point and that the amplitude of Ca 2+ transients decreases with increasing distance. It is noteworthy that a comparison of signaling characteristics with respect to the degree of cataract progression revealed that, in LCs from more developed cataracts, the Ca 2+ wave propagates faster and the amplitudes of Ca 2+ signals are lower, while their durations are longer. No differences were identified when comparing LCs with regard to the cataract type. Moreover, experiments with Apyrase have revealed that the Ca 2+ signals are not affected by ATP-dependent paracrine communication. Our results indicated that cataract progression is associated with modifications in Ca 2+ signaling in LECs, suggesting the functional importance of altered Ca 2+ signaling of LECs in cataractogenesis.

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.)

Social distancing is an effective strategy to mitigate the impact of infectious diseases. If sick or healthy, or both, predominantly socially distance, the epidemic curve flattens. Contact reductions may occur for different reasons during a pandemic including health-related mobility loss (severity of symptoms), duty of care for a member of a high-risk group, and forced quarantine. Other decisions to reduce contacts are of a more voluntary nature. In particular, sick people reduce contacts consciously to avoid infecting others, and healthy individuals reduce contacts in order to stay healthy. We use game theory to formalize the interaction of voluntary social distancing in a partially infected population. This improves the behavioral micro-foundations of epidemiological models, and predicts differential social distancing rates dependent on health status. The model's key predictions in terms of comparative statics are derived, which concern changes and interactions between social distancing behaviors of sick and healthy. We fit the relevant parameters for endogenous social distancing to an epidemiological model with evidence from influenza waves to provide a benchmark for an epidemic curve with endogenous social distancing. Our results suggest that spreading similar in peak and case numbers to what partial immobilization of the population produces, yet quicker to pass, could occur endogenously. Going forward, eventual social distancing orders and lockdown policies should be benchmarked against more realistic epidemic models that take endogenous social distancing into account, rather than be driven by static, and therefore unrealistic, estimates for social mixing that intrinsically overestimate spreading.

Type 2 diabetes mellitus is a complex multifactorial disease of epidemic proportions. It involves genetic and lifestyle factors that lead to dysregulations in hormone secretion and metabolic homeostasis. Accumulating evidence indicates that altered mitochondrial structure, function, and particularly bioenergetics of cells in different tissues have a central role in the pathogenesis of type 2 diabetes mellitus. In the present study, we explore how mitochondrial dysfunction impairs the coupling between metabolism and exocytosis in the pancreatic alpha and beta cells. We demonstrate that reduced mitochondrial ATP production is linked with the observed defects in insulin and glucagon secretion by utilizing computational modeling approach. Specifically, a 30-40% reduction in alpha cells' mitochondrial function leads to a pathological shift of glucagon secretion, characterized by oversecretion at high glucose concentrations and insufficient secretion in hypoglycemia. In beta cells, the impaired mitochondrial energy metabolism is accompanied by reduced insulin secretion at all glucose levels, but the differences, compared to a normal beta cell, are the most pronounced in hyperglycemia. These findings improve our understanding of metabolic pathways and mitochondrial bioenergetics in the pathology of type 2 diabetes mellitus and might help drive the development of innovative therapies to treat various metabolic diseases.

Thermoregulation plays a vital role in homeostasis. Many species of animals as well as humans have evolved various physiological mechanisms for body temperature control, which are characteristically flexible and enable a fine-tuned spatial and temporal regulation of body temperature in different environmental conditions and circumstances. Human beings normally maintain a core body temperature at around 37°C, and maintenance of this relatively high temperature is critical for survival. Therefore, principles of thermoregulatory control have also important clinical implications. Infections can cause the body temperature to rise internally and several diseases can cause a dysfunction of thermoregulatory mechanisms. Moreover, the utilization of thermotherapies in treating various diseases has been known for thousands of years with a recent resurgence of interest. An increasing amount of research suggests that targeted temperature management is of paramount importance to patient outcomes in certain clinical scenarios. We provide a concise summary of the basic concepts of thermoregulation. Emphasis is given to the principles of thermoregulation in humans in basic pathological states and to targeted temperature management strategies in the clinical environment, with special attention on therapeutic hypothermia in postcardiac arrest patients. Finally, the discussion is focused on the potential offered by computational thermophysiological models for predicting thermal responses of patients in various clinical circumstances, for proposing new perspectives in the design of novel thermal therapies, and to optimize targeted temperature management strategies. This article is categorized under: Cardiovascular Diseases > Cardiovascular Diseases>Computational Models Cardiovascular Diseases > Cardiovascular Diseases>Environmental Factors Cardiovascular Diseases > Cardiovascular Diseases>Biomedical Engineering., (© 2020 Wiley Periodicals LLC.)

The human brain is increasingly seen as a dynamic neural system, the function of which relies on a diverse set of connections between brain regions. To assess these complex dynamical interactions, formalism of complex networks was suggested as one of the most promising tools to offer new insight into the brain's structural and functional organization, with a potential also for clinical implications. Irrespective of the brain mapping technique, modern network approaches have revealed fundamental aspects of normal brain-network organization, such as small-world and scale-free patterns, hierarchical modularity, and the presence of hubs. Moreover, the utility of these approaches, to gain a better understanding of neurological diseases, is of great interest. In the present contribution, we first describe the basic network measures and how the brain networks are constructed on the basis of brain activity data in order to introduce clinical neurologists to this new theoretical paradigm. We then demonstrate how network formalism can be used to detect changes in EEG-based functional connectivity patterns in six paediatric patients with childhood absence epilepsy. Notably, our results do not only indicate enhanced synchronicity during epileptic episodes but also reveal specific spatial changes in the electrical activity of the brain. We argue that the network-based evaluation of functional brain networks can provide clinicians with more detailed insight into the activity of a pathological brain and can also be regarded as a support for objective diagnosis and treatment for various neurological diseases.

If interaction partners in social dilemma games are not selected randomly from the population but are instead determined by a network of contacts, it has far reaching consequences for the evolutionary dynamics. Selecting partners randomly leads to a well-mixed population, where pattern formation is essentially impossible. This rules out important mechanisms that can facilitate cooperation, most notably network reciprocity. In contrast, if interactions are determined by a lattice or a network, then the population is said to be structured, where cooperators can form compact clusters that protect them from invading defectors. Between these two extremes, however, there is ample middle ground that can be brought about by the consideration of temporal networks, mobility, or other coevolutionary processes. The question that we here seek to answer is, when does mixing on a lattice actually lead to well-mixed conditions? To that effect, we use the public goods game on a square lattice, and we consider nearest-neighbor and random mixing with different frequencies, as well as a mix of both mixing protocols. Not surprisingly, we find that nearest-neighbor mixing requires a higher frequency than random mixing to arrive at the well-mixed limit. The differences between the two mixing protocols are most expressed at intermediate mixing frequencies, whilst at very low and very high mixing frequencies the two almost converge. We also find a near universal exponential growth of the average size of cooperator clusters as their fraction increases from zero to one, regardless of whether this increase is due to increasing the multiplication factor of the public goods, decreasing the frequency of mixing, or gradually shifting the mixing from random to nearest neighbors.

We present a mathematical model of the energy-driven metabolic switch for glucagon and insulin secretion from pancreatic alpha and beta cells, respectively. The energy status related to hormone secretion is studied for various glucose concentrations. Additionally, the physiological response is studied with regards to the presence of other metabolites, particularly the free-fatty acids. At low glucose, the ATP production in alpha cells is high due to free-fatty acids oxidation in mitochondria, which enables glucagon secretion. When the glucose concentration is elevated above the threshold value, the glucagon secretion is switched off due to the contribution of glycolytic ATP production, representing an "anaerobic switch". On the other hand, during hypoglycemia, the ATP production in beta cells is low, reflecting a "waiting state" for glucose as the main metabolite. When glucose is elevated above the threshold value, the oxidative fate of glucose in mitochondria is the main source of energy required for effective insulin secretion, i.e. the "aerobic switch". Our results show the importance of well-regulated and fine-tuned energetic processes in pancreatic alpha and beta cells required for efficient hormone secretion and hence effective blood glucose regulation. These energetic processes have to be appropriately switched on and off based on the sensing of different metabolites by alpha and beta cells. Our computational results indicate that disturbances in cell energetics (e.g. mitochondrial dysfunction), and dysfunctional metabolite sensing and distribution throughout the cell might be related to pathologies such as metabolic syndrome and diabetes., (Copyright © 2020. Published by Elsevier Ltd.)

Purpose: To determine whether the psychological state of patients with head and neck cancer (HCN) is associated with their nutritional status., Methods: In 40 patients with locally advanced HNC treated with definitive or adjuvant (chemo)radiotherapy, psychological and nutritional status were assessed before treatment, at its completion and 3 months' post-therapy. Psychosocial distress was measured using the Hospital Anxiety and Depression Scale questionnaire (HADS-A, HADS-D), whereas the nutritional status was evaluated using standard methods (Nutritional Risk Screening Tool 2002, anthropometric data, dynamometry and laboratory tests) and with a bioelectrical impedance analysis parameter phase angle (PA)., Results: Before treatment, more patients were screened positive for anxiety than at treatment completion (p = 0.037) or 3 months' post-therapy (p = 0.083). Depression prevalence was non-significantly higher at the end and after therapy. Compared to the baseline, more cachectic patients and a reduction of PA values were found at successive assessments. Anxiety was more often recorded among malnourished/cachectic patients (assessment 1, p = 0.017; assessment 2, p = 0.020) who were also found more frequently depressed (assessment 2, p = 0.045; assessment 3, p = 0.023). Significantly higher PA values were measured in patients without distress determined at 3 months' post-therapy by the HADS-A (p = 0.027)., Conclusion: The association between the psychological and nutritional status found in this pilot study and the options for intervention warrants further clarification in a larger prospective trial.

Type 2 diabetes mellitus (T2DM) has been associated with insulin resistance and the failure of β-cells to produce and secrete enough insulin as the disease progresses. However, clinical treatments based solely on insulin secretion and action have had limited success. The focus is therefore shifting towards α-cells, in particular to the dysregulated secretion of glucagon. Our qualitative electron-microscopy-based observations gave an indication that mitochondria in α-cells are altered in Western-diet-induced T2DM. In particular, α-cells extracted from mouse pancreatic tissue showed a lower density of mitochondria, a less expressed matrix and a lower number of cristae. These deformities in mitochondrial ultrastructure imply a decreased efficiency in mitochondrial ATP production, which prompted us to theoretically explore and clarify one of the most challenging problems associated with T2DM, namely the lack of glucagon secretion in hypoglycaemia and its oversecretion at high blood glucose concentrations. To this purpose, we constructed a novel computational model that links α-cell metabolism with their electrical activity and glucagon secretion. Our results show that defective mitochondrial metabolism in α-cells can account for dysregulated glucagon secretion in T2DM, thus improving our understanding of T2DM pathophysiology and indicating possibilities for new clinical treatments., Competing Interests: M.P. was a member of the Royal Society Open Science Editorial Board at the time of submission and publication of this paper, but was not involved in the assessment of the manuscript., (© 2020 The Authors.)

Signatures of criticality, such as power law scaling of observables, have been empirically found in a plethora of real-life settings, including biological systems. The presence of critical states is believed to have many functional advantages and is associated with optimal operational abilities. Typically, critical dynamics arises in the proximity of phase transition points between absorbing disordered states (subcriticality) and ordered active regimes (supercriticality) and requires a high degree of fine tuning to emerge, which is unlikely to occur in real biological systems. In the present study we propose a rather simple, and biologically relevant mechanism that profoundly expands the critical-like region. In particular, by means of numerical simulation we show that incorporating spatial heterogeneities into the square lattice of map-based excitable oscillators broadens the parameter space in which the distribution of excitation wave sizes follows closely a power law. Most importantly, this behavior is only observed if the spatial profile exhibits intermediate-sized patches with similar excitability levels, whereas for large and small spatial clusters only marginal widening of the critical state is detected. Furthermore, it turned out that the presence of spatial disorder in general amplifies the size of excitation waves, whereby the relatively highest contributions are observed in the proximity of the critical point. We argue that the reported mechanism is of particular importance for excitable systems with local interactions between individual elements.

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.

Network science is today established as a backbone for description of structure and function of various physical, chemical, biological, technological, and social systems. Here we review recent advances in the study of complex biological systems that were inspired and enabled by methods of network science. First, we present research highlights ranging from determination of the molecular interaction network within a cell to studies of architectural and functional properties of brain networks and biological transportation networks. Second, we focus on synergies between network science and data analysis, which enable us to determine functional connectivity patterns in multicellular systems. Until now, this intermediate scale of biological organization received the least attention from the network perspective. As an example, we review the methodology for the extraction of functional beta cell networks in pancreatic islets of Langerhans by means of advanced imaging techniques. Third, we concentrate on the emerging field of multilayer networks and review the first endeavors and novel perspectives offered by this framework in exploring biological complexity. We conclude by outlining challenges and directions for future research that encompass utilization of the multilayer network formalism in exploring intercellular communication patterns in tissues, and we advocate for network science being one of the key pillars for assessing physiological function of complex biological systems-from organelles to organs-in health and disease., (Copyright © 2017 Elsevier B.V. All rights reserved.)

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.

SNAP-25 is a protein of the core SNARE complex mediating stimulus-dependent release of insulin from pancreatic β cells. The protein exists as two alternatively spliced isoforms, SNAP-25a and SNAP-25b, differing in 9 out of 206 amino acids, yet their specific roles in pancreatic β cells remain unclear. We explored the effect of SNAP-25b-deficiency on glucose-stimulated insulin release in islets and found increased secretion both in vivo and in vitro. However, slow photo-release of caged Ca 2+ in β cells within pancreatic slices showed no significant differences in Ca 2+ -sensitivity, amplitude or rate of exocytosis between SNAP-25b-deficient and wild-type littermates. Therefore, we next investigated if Ca 2+ handling was affected in glucose-stimulated β cells using intracellular Ca 2+ -imaging and found premature activation and delayed termination of [Ca 2+ ] i elevations. These findings were accompanied by less synchronized Ca 2+ -oscillations and hence more segregated functional β cell networks in SNAP-25b-deficient mice. Islet gross morphology and architecture were maintained in mutant mice, although sex specific compensatory changes were observed. Thus, our study proposes that SNAP-25b in pancreatic β cells, except for participating in the core SNARE complex, is necessary for accurate regulation of Ca 2+ -dynamics.

Quantitative analysis of the vascular network anatomy is critical for the understanding of the vasculature structure and function. In this study, we have combined microcomputed tomography (microCT) and computational analysis to provide quantitative three-dimensional geometrical and topological characterization of the normal kidney vasculature, and to investigate how 2 core genes of the Wnt/planar cell polarity, Frizzled4 and Frizzled6, affect vascular network morphogenesis. Experiments were performed on frizzled4 (Fzd4-/-) and frizzled6 (Fzd6-/-) deleted mice and littermate controls (WT) perfused with a contrast medium after euthanasia and exsanguination. The kidneys were scanned with a high-resolution (16 μm) microCT imaging system, followed by 3D reconstruction of the arterial vasculature. Computational treatment includes decomposition of 3D networks based on Diameter-Defined Strahler Order (DDSO). We have calculated quantitative (i) Global scale parameters, such as the volume of the vasculature and its fractal dimension (ii) Structural parameters depending on the DDSO hierarchical levels such as hierarchical ordering, diameter, length and branching angles of the vessel segments, and (iii) Functional parameters such as estimated resistance to blood flow alongside the vascular tree and average density of terminal arterioles. In normal kidneys, fractal dimension was 2.07±0.11 (n = 7), and was significantly lower in Fzd4-/- (1.71±0.04; n = 4), and Fzd6-/- (1.54±0.09; n = 3) kidneys. The DDSO number was 5 in WT and Fzd4-/-, and only 4 in Fzd6-/-. Scaling characteristics such as diameter and length of vessel segments were altered in mutants, whereas bifurcation angles were not different from WT. Fzd4 and Fzd6 deletion increased vessel resistance, calculated using the Hagen-Poiseuille equation, for each DDSO, and decreased the density and the homogeneity of the distal vessel segments. Our results show that our methodology is suitable for 3D quantitative characterization of vascular networks, and that Fzd4 and Fzd6 genes have a deep patterning effect on arterial vessel morphogenesis that may determine its functional efficiency.

In this work we investigated how modifications of the Ca2+ homeostasis in anterior lens epithelial cells (LECs) are associated with different types of cataract (cortical or nuclear) and how the progression of the cataract (mild or moderate) affects the Ca2+ signaling. We systematically analyzed different aspects of intra- and inter-cellular Ca2+ signaling in the human LECs, which are attached to surgically isolated lens capsule (LC), obtained during cataract surgery. We monitored the temporal and spatial changes in intracellular Ca2+ concentration after stimulation with acetylcholine by means of Fura-2 fluorescence captured with an inverted microscope. In our analysis we compared the features of Ca2+ signals in individual cells, synchronized activations, spatio-temporal grouping and the nature of intercellular communication between LECs. The latter was assessed by using the methodologies of the complex network theory. Our results point out that at the level of individual cells there are no significant differences when comparing the features of the signals with regard either to the type or the stage of the cataract. On the other hand, noticeable differences are observed at the multicellular level, despite inter-capsule variability. LCs associated with more developed cataracts were found to exhibit a slower collective response to stimulation, a less pronounced spatio-temporal clustering of LECs with similar signaling characteristics. The reconstructed intercellular networks were found to be sparser and more segregated than in LCs associated with mild cataracts. Moreover, we show that spontaneously active LECs often operate in localized groups with quite well aligned Ca2+ activity. The presence of spontaneous activity was also found to affect the stimulated Ca2+ responses of individual cells. Our findings indicate that the cataract progression entails the impairment of intercellular signaling thereby suggesting the functional importance of altered Ca2+ signaling of LECs in cataractogenesis.

Beta cells in the pancreatic islets of Langerhans are precise biological sensors for glucose and play a central role in balancing the organism between catabolic and anabolic needs. A hallmark of the beta cell response to glucose are oscillatory changes of membrane potential that are tightly coupled with oscillatory changes in intracellular calcium concentration which, in turn, elicit oscillations of insulin secretion. Both membrane potential and calcium changes spread from one beta cell to the other in a wave-like manner. In order to assess the properties of the abovementioned responses to physiological and pathological stimuli, the main challenge remains how to effectively measure membrane potential and calcium changes at the same time with high spatial and temporal resolution, and also in as many cells as possible. To date, the most wide-spread approach has employed the electrophysiological patch-clamp method to monitor membrane potential changes. Inherently, this technique has many advantages, such as a direct contact with the cell and a high temporal resolution. However, it allows one to assess information from a single cell only. In some instances, this technique has been used in conjunction with CCD camera-based imaging, offering the opportunity to simultaneously monitor membrane potential and calcium changes, but not in the same cells and not with a reliable cellular or subcellular spatial resolution. Recently, a novel family of highly-sensitive membrane potential reporter dyes in combination with high temporal and spatial confocal calcium imaging allows for simultaneously detecting membrane potential and calcium changes in many cells at a time. Since the signals yielded from both types of reporter dyes are inherently noisy, we have developed complex methods of data denoising that permit for visualization and pixel-wise analysis of signals. Combining the experimental approach of high-resolution imaging with the advanced analysis of noisy data enables novel physiological insights and reassessment of current concepts in unprecedented detail.