Intrinsic islet heterogeneity and gap junction coupling determine spatiotemporal Ca²⁺ wave dynamics

Exploring pancreatic beta-cell subgroups and their connectivity

Functional pancreatic islet beta cells are essential to ensure glucose homeostasis across species from zebrafish to humans. These cells show significant heterogeneity, and emerging studies have revealed that connectivity across a hierarchical network is required for normal insulin release. Here, we discuss current thinking and areas of debate around intra-islet connectivity, cellular hierarchies and potential "controlling" beta-cell populations. We focus on methodologies, including comparisons of different cell preparations as well as in vitro and in vivo approaches to imaging and controlling the activity of human and rodent islet preparations. We also discuss the analytical approaches that can be applied to live-cell data to identify and study critical subgroups of cells with a disproportionate role in control Ca2+ dynamics and thus insulin secretion (such as "first responders", "leaders" and "hubs", as defined by Ca2+ responses to glucose stimulation). Possible mechanisms by which this hierarchy is achieved, its physiological relevance and how its loss may contribute to islet failure in diabetes mellitus are also considered. A glossary of terms and links to computational resources are provided.

Network representation of multicellular activity in pancreatic islets: Technical considerations for functional connectivity analysis

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.

Integrating single-cell transcriptomics with cellular phenotypes: cell morphology, Ca2+ imaging and electrophysiology

I review recent technological advancements in coupling single-cell transcriptomics with cellular phenotypes including morphology, calcium signaling, and electrophysiology. Single-cell RNA sequencing (scRNAseq) has revolutionized cell type classifications by capturing the transcriptional diversity of cells. A new wave of methods to integrate scRNAseq and biophysical measurements is facilitating the linkage of transcriptomic data to cellular function, which provides physiological insight into cellular states. I briefly discuss critical factors of these phenotypical characterizations such as timescales, information content, and analytical tools. Dedicated sections focus on the integration with cell morphology, calcium imaging, and electrophysiology (patch-seq), emphasizing their complementary roles. I discuss their application in elucidating cellular states, refining cell type classifications, and uncovering functional differences in cell subtypes. To illustrate the practical applications and benefits of these methods, I highlight their use in tissues with excitable cell-types such as the brain, pancreatic islets, and the retina. The potential of combining functional phenotyping with spatial transcriptomics for a detailed mapping of cell phenotypes in situ is explored. Finally, I discuss open questions and future perspectives, emphasizing the need for a shift towards broader accessibility through increased throughput.

Conversion of spikers to bursters in pituitary cell networks: Is it better to disperse for maximum exposure or circle the wagons?

The endocrine cells of the pituitary gland are electrically active, and in vivo they form small networks where the bidirectional cell-cell coupling is through gap junctions. Numerous studies of dispersed pituitary cells have shown that typical behaviors are tonic spiking and bursting, the latter being more effective at evoking secretion. In this article, we use mathematical modeling to examine the dynamics of small networks of spiking and bursting pituitary cells. We demonstrate that intrinsic bursting cells are capable of converting intrinsic spikers into bursters, and perform a fast/slow analysis to show why this occurs. We then demonstrate the sensitivity of network dynamics to the placement of bursting cells within the network, and demonstrate strategies that are most effective at maximizing secretion from the population of cells. This study provides insights into the in vivo behavior of cells such as the stress-hormone-secreting pituitary corticotrophs that are switched from spiking to bursting by hypothalamic neurohormones. While much is known about the electrical properties of these cells when isolated from the pituitary, how they behave when part of an electrically coupled network has been largely unstudied.

A primer on modelling pancreatic islets: from models of coupled β-cells to multicellular islet models

Pancreatic islets are mini-organs composed of hundreds or thousands of ɑ, β and δ-cells, which, respectively, secrete glucagon, insulin and somatostatin, key hormones for the regulation of blood glucose. In pancreatic islets, hormone secretion is tightly regulated by both internal and external mechanisms, including electrical communication and paracrine signaling between islet cells. Given its complexity, the experimental study of pancreatic islets has been complemented with computational modeling as a tool to gain a better understanding about how all the mechanisms involved at different levels of organization interact. In this review, we describe how multicellular models of pancreatic cells have evolved from the early models of electrically coupled β-cells to models in which experimentally derived architectures and both electrical and paracrine signals have been considered.

Deconstructing the integrated oscillator model for pancreatic β-cells

Electrical bursting oscillations in the β-cells of pancreatic islets have been a focus of investigation for more than fifty years. This has been aided by mathematical models, which are descendants of the pioneering Chay-Keizer model. This article describes the key biophysical and mathematical elements of this model, and then describes the path forward from there to the Integrated Oscillator Model (IOM). It is both a history and a deconstruction of the IOM that describes the various elements that have been added to the model over time, and the motivation for adding them. Finally, the article is a celebration of the 40th anniversary of the publication of the Chay-Keizer model.

Multilevel synchronization of human β-cells networks

β-cells within the endocrine pancreas are fundamental for glucose, lipid and protein homeostasis. Gap junctions between cells constitute the primary coupling mechanism through which cells synchronize their electrical and metabolic activities. This evidence is still only partially investigated through models and numerical simulations. In this contribution, we explore the effect of combined electrical and metabolic coupling in β-cell clusters using a detailed biophysical model. We add heterogeneity and stochasticity to realistically reproduce β-cell dynamics and study networks mimicking arrangements of β-cells within human pancreatic islets. Model simulations are performed over different couplings and heterogeneities, analyzing emerging synchronization at the membrane potential, calcium, and metabolites levels. To describe network synchronization, we use the formalism of multiplex networks and investigate functional network properties and multiplex synchronization motifs over the structural, electrical, and metabolic layers. Our results show that metabolic coupling can support slow wave propagation in human islets, that combined electrical and metabolic synchronization is realized in small aggregates, and that metabolic long-range correlation is more pronounced with respect to the electrical one.

The effect of forskolin and the role of Epac2A during activation, activity, and deactivation of beta cell networks

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.

Pancreatic β-cell heterogeneity in adult human islets and stem cell-derived islets

Recent studies reported that pancreatic β-cells are heterogeneous in terms of their transcriptional profiles and their abilities for insulin secretion. Sub-populations of pancreatic β-cells have been identified based on the functionality and expression of specific surface markers. Under diabetes condition, β-cell identity is altered leading to different β-cell sub-populations. Furthermore, cell-cell contact between β-cells and other endocrine cells within the islet play an important role in regulating insulin secretion. This highlights the significance of generating a cell product derived from stem cells containing β-cells along with other major islet cells for treating patients with diabetes, instead of transplanting a purified population of β-cells. Another key question is how close in terms of heterogeneity are the islet cells derived from stem cells? In this review, we summarize the heterogeneity in islet cells of the adult pancreas and those generated from stem cells. In addition, we highlight the significance of this heterogeneity in health and disease conditions and how this can be used to design a stem cell-derived product for diabetes cell therapy.

Role of Delta/Notch-like EGF-related receptor in blood glucose homeostasis

Cell-cell interactions are necessary for optimal endocrine functions in the pancreas. β-cells, characterized by the expression and secretion of the hormone insulin, are a major constituent of functional micro-organs in the pancreas known as islets of Langerhans. Cell-cell contacts between β-cells are required to regulate insulin production and glucose-stimulated insulin secretion, which are key determinants of blood glucose homeostasis. Contact-dependent interactions between β-cells are mediated by gap junctions and cell adhesion molecules such as E-cadherin and N-CAM. Recent genome-wide studies have implicated Delta/Notch-like EGF-related receptor (Dner) as a potential susceptibility locus for Type 2 Diabetes in humans. DNER is a transmembrane protein and a proposed Notch ligand. DNER has been implicated in neuron-glia development and cell-cell interactions. Studies herein demonstrate that DNER is expressed in β-cells with an onset during early postnatal life and sustained throughout adulthood in mice. DNER loss in adult β-cells in mice (β-Dner cKO mice) disrupted islet architecture and decreased the expression of N-CAM and E-cadherin. β-Dner cKO mice also exhibited impaired glucose tolerance, defects in glucose- and KCl-induced insulin secretion, and decreased insulin sensitivity. Together, these studies suggest that DNER plays a crucial role in mediating islet cell-cell interactions and glucose homeostasis.

Spatial distribution of heterogeneity as a modulator of collective dynamics in pancreatic beta-cell networks and beyond

We study the impact of spatial distribution of heterogeneity on collective dynamics in gap-junction coupled beta-cell networks comprised on cells from two populations that differ in their intrinsic excitability. Initially, these populations are uniformly and randomly distributed throughout the networks. We develop and apply an iterative algorithm for perturbing the arrangement of the network such that cells from the same population are increasingly likely to be adjacent to one another. We find that the global input strength, or network drive, necessary to transition the network from a state of quiescence to a state of synchronised and oscillatory activity decreases as network sortedness increases. Moreover, for weak coupling, we find that regimes of partial synchronisation and wave propagation arise, which depend both on network drive and network sortedness. We then demonstrate the utility of this algorithm for studying the distribution of heterogeneity in general networks, for which we use Watts-Strogatz networks as a case study. This work highlights the importance of heterogeneity in node dynamics in establishing collective rhythms in complex, excitable networks and has implications for a wide range of real-world systems that exhibit such heterogeneity.

Functional characteristics of hub and wave-initiator cells in β cell networks

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 Ca2+ waves, whereas wave-initiator cells trigger intercellular signals in their cohorts. Here, we determined Ca2+ signaling characteristics of these two β cell subpopulations and the relationship between them by means of functional multicellular Ca2+ imaging in mouse pancreatic tissue slices in combination with methods of complex network theory. We constructed network layers based on individual Ca2+ 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, Ca2+ 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.

Determinants and dynamics of pancreatic islet architecture

The islets of Langerhans are highly organized structures that have species-specific, three-dimensional tissue architecture. Islet architecture is critical for proper hormone secretion in response to nutritional stimuli. Islet architecture is disrupted in all types of diabetes mellitus and in cadaveric islets for transplantation during isolation, culture, and perfusion, limiting patient outcomes. Moreover, recapitulating native islet architecture remains a key challenge for in vitro generation of islets from stem cells. In this review, we discuss work that has led to the current understanding of determinants of pancreatic islet architecture, and how this architecture is maintained or disrupted during tissue remodeling in response to normal and pathological metabolic changes. We further discuss both empirical and modeling data that highlight the importance of islet architecture for islet function.

Coordination of pancreatic islet rhythmic activity by delayed negative feedback

Secretion of insulin from the pancreas is pulsatile, driven by intrinsic oscillations within individual islets of Langerhans. The secretions are coordinated among the many islets distributed throughout the pancreas producing a synchronized rhythm in vivo that is essential for maintaining normal glucose levels. One hypothesized mechanism for the coordination of islet activity is negative feedback, whereby sequestration of glucose in response to elevated insulin leads to a reduction in the blood glucose level that is sensed by the islet population. This global signal of glucose then coordinates the individual islets. In this study, we tested how this coordination mechanism is affected by time delays in the negative feedback, using a microfluidic system to monitor Ca2+ levels in a small population of islets and implementing glucose control through a negative feedback system. We found that islet synchronization occurs even with time delays in the feedback of up to 7 min. We also found that a second, slower closed-loop oscillation period is produced during delayed feedback in which islet oscillations are clustered into episodes. The period of this second oscillatory mode increases with the time delay and appears to be a second stable behavior that coexists with the faster synchronized oscillation. The general conclusion is that islet coordination through negative feedback is a viable means of islet coordination that is robust to delays in the timing of the feedback, and could complement other potential coordination mechanisms such as entrainment by pancreatic ganglia.NEW & NOTEWORTHY Insulin secretion from islets of Langerhans is rhythmic, and these rhythms are coordinated to produce oscillatory plasma insulin levels. Using a combination of microfluidics and computational modeling, we demonstrate that coordination can occur through negative feedback of the type provided by the liver, even if that feedback is delayed by several minutes. We also demonstrate that a second, slower, mode of oscillations can occur when feedback is delayed where faster oscillations are grouped into episodes.

Ca2+ Oscillations, Waves, and Networks in Islets From Human Donors With and Without Type 2 Diabetes

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.

Functional architecture of pancreatic islets identifies a population of first responder cells that drive the first-phase calcium response

Insulin-secreting β-cells are functionally heterogeneous. Whether there exist cells driving the first-phase calcium response in individual islets, has not been examined. We examine “first responder” cells, defined by the earliest [Ca²⁺] response during first-phase [Ca²⁺] elevation, distinct from previously identified “hub” and “leader” cells. We used islets isolated from Mip-Cre^ER; Rosa-Stop-Lox-Stop-GCamP6s mice (β-GCamP6s) that show β-cell-specific GCamP6s expression following tamoxifen-induced CreER-mediated recombination. First responder cells showed characteristics of high membrane excitability and lower electrical coupling to their neighbors. The first-phase response time of β-cells in the islet was spatially organized, dependent on the cell’s distance to the first responder cell, and consistent over time up to approximately 24 h. When first responder cells were laser ablated, the first-phase [Ca²⁺] was slowed down, diminished, and discoordinated compared to random cell ablation. Cells that were next earliest to respond often took over the role of the first responder upon ablation. In summary, we discover and characterize a distinct first responder β-cell state, critical for the islet first-phase response to glucose.

A hallmark of the early stages in diabetes is the disruption of the first-phase peak of insulin secretion. This study reveals a state of beta cells that drives the first-phase calcium response (a precursor to insulin secretion), challenging the existing paradigm that first-phase calcium response to glucose is not organized within the islet.

Pancreatic α and β cells are globally phase-locked

The Ca²⁺ modulated pulsatile glucagon and insulin secretions by pancreatic α and β cells play a crucial role in glucose homeostasis. However, how α and β cells coordinate to produce various Ca²⁺ oscillation patterns is still elusive. Using a microfluidic device and transgenic mice, we recorded Ca²⁺ signals from islet α and β cells, and observed heterogeneous Ca²⁺ oscillation patterns intrinsic to each islet. After a brief period of glucose stimulation, α and β cells’ oscillations were globally phase-locked. While the activation of α cells displayed a fixed time delay of ~20 s to that of β cells, β cells activated with a tunable period. Moreover, islet α cell number correlated with oscillation frequency. We built a mathematical model of islet Ca²⁺ oscillation incorporating paracrine interactions, which quantitatively agreed with the experimental data. Our study highlights the importance of cell-cell interaction in generating stable but tunable islet oscillation patterns.

The Ca²⁺ modulated pulsatile glucagon and insulin secretions by pancreatic α and β cells are critical in glucose homeostasis. Here the authors show that the Ca²⁺ oscillations of α and β cells are phase-locked, and that the oscillation pattern is tuned by paracrine interactions between α and β cells.

From Isles of Königsberg to Islets of Langerhans: Examining the Function of the Endocrine Pancreas Through Network Science

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.

Dynamic changes in β-cell [Ca2+] regulate NFAT activation, gene transcription, and islet gap junction communication

OBJECTIVE

Diabetes occurs because of insufficient insulin secretion due to β-cell dysfunction within the islet of Langerhans. Elevated glucose levels trigger β-cell membrane depolarization, action potential generation, and slow sustained free-Ca2+ ([Ca2+]) oscillations, which trigger insulin release. Nuclear factor of activated T-cell (NFAT) is a transcription factor, which is regulated by the increases in [Ca2+] and calceineurin (CaN) activation. NFAT regulation links cell activity with gene transcription in many systems and regulates proliferation and insulin granule biogenesis within the β-cell. However, the link between the regulation of β-cell electrical activity and oscillatory [Ca2+] dynamics with NFAT activation and downstream transcription is poorly understood. Here, we tested whether dynamic changes to β-cell electrical activity and [Ca2+] regulate NFAT activation and downstream transcription.

METHODS

In cell lines, mouse islets, and human islets, including those from donors with type 2 diabetes, we applied both agonists/antagonists of ion channels together with optogenetics to modulate β-cell electrical activity. We measured the dynamics of [Ca2+] and NFAT activation as well as performed whole transcriptome and functional analyses.

RESULTS

Both glucose-induced membrane depolarization and optogenetic stimulation triggered NFAT activation as well as increased the transcription of NFAT targets and intermediate early genes (IEGs). Importantly, slow, sustained [Ca2+] oscillation conditions led to NFAT activation and downstream transcription. In contrast, in human islets from donors with type2 diabetes, NFAT activation by glucose was diminished, but rescued upon pharmacological stimulation of electrical activity. NFAT activation regulated GJD2 expression and increased Cx36 gap junction permeability upon elevated oscillatory [Ca2+] dynamics. However, it is unclear if NFAT directly binds the GJD2 gene to regulate expression.

CONCLUSIONS

This study provides an insight into the specific patterns of electrical activity that regulate NFAT activation, gene transcription, and islet function. In addition, it provides information on how these factors are disrupted in diabetes.

The physiological role of β-cell heterogeneity in pancreatic islet function

Endocrine cells within the pancreatic islets of Langerhans are heterogeneous in terms of transcriptional profile, protein expression and the regulation of hormone release. Even though this heterogeneity has long been appreciated, only within the past 5 years have detailed molecular analyses led to an improved understanding of its basis. Although we are beginning to recognize why some subpopulations of endocrine cells are phenotypically different to others, arguably the most important consideration is how this heterogeneity affects the regulation of hormone release to control the homeostasis of glucose and other energy-rich nutrients. The focus of this Review is the description of how endocrine cell heterogeneity (and principally that of insulin-secreting β-cells) affects the regulation of hormone secretion within the islets of Langerhans. This discussion includes an overview of the functional characteristics of the different islet cell subpopulations and describes how they can communicate to influence islet function under basal and glucose-stimulated conditions. We further discuss how changes to the specific islet cell subpopulations or their numbers might underlie islet dysfunction in type 2 diabetes mellitus. We conclude with a discussion of several key open questions regarding the physiological role of islet cell heterogeneity.

β-Cell Pathophysiology: A Review of Advanced Optical Microscopy Applications

β-cells convert glucose (input) resulting in the controlled release of insulin (output), which in turn has the role to maintain glucose homeostasis. β-cell function is regulated by a complex interplay between the metabolic processing of the input, its transformation into second-messenger signals, and final mobilization of insulin-containing granules towards secretion of the output. Failure at any level in this process marks β-cell dysfunction in diabetes, thus making β-cells obvious potential targets for therapeutic purposes. Addressing quantitatively β-cell (dys)function at the molecular level in living samples requires probing simultaneously the spatial and temporal dimensions at the proper resolution. To this aim, an increasing amount of research efforts are exploiting the potentiality of biophysical techniques. In particular, using excitation light in the visible/infrared range, a number of optical-microscopy-based approaches have been tailored to the study of β-cell-(dys)function at the molecular level, either in label-free mode (i.e., exploiting intrinsic autofluorescence of cells) or by the use of organic/genetically-encoded fluorescent probes. Here, relevant examples from the literature are reviewed and discussed. Based on this, new potential lines of development in the field are drawn.

Microtubules regulate pancreatic β-cell heterogeneity via spatiotemporal control of insulin secretion hot spots

Heterogeneity of glucose-stimulated insulin secretion (GSIS) in pancreatic islets is physiologically important but poorly understood. Here, we utilize mouse islets to determine how microtubules (MTs) affect secretion toward the vascular extracellular matrix at single cell and subcellular levels. Our data indicate that MT stability in the β-cell population is heterogenous, and that GSIS is suppressed in cells with highly stable MTs. Consistently, MT hyper-stabilization prevents, and MT depolymerization promotes the capacity of single β-cell for GSIS. Analysis of spatiotemporal patterns of secretion events shows that MT depolymerization activates otherwise dormant β-cells via initiation of secretion clusters (hot spots). MT depolymerization also enhances secretion from individual cells, introducing both additional clusters and scattered events. Interestingly, without MTs, the timing of clustered secretion is dysregulated, extending the first phase of GSIS and causing oversecretion. In contrast, glucose-induced Ca²⁺ influx was not affected by MT depolymerization yet required for secretion under these conditions, indicating that MT-dependent regulation of secretion hot spots acts in parallel with Ca²⁺ signaling. Our findings uncover a novel MT function in tuning insulin secretion hot spots, which leads to accurately measured and timed response to glucose stimuli and promotes functional β-cell heterogeneity.

Glucose-induced [Ca2+]i oscillations in β cells are composed of trains of spikes within a subplasmalemmal microdomain

Electrical activity and oscillations of cytosolic Ca²⁺ concentrations ([Ca²⁺]_i) that trigger insulin release in response to glucose are key functions of pancreatic β cells. Although oscillatory Ca²⁺ signals have been intensively studied in β cells, their lower frequency did not match that of electrical activity. In addition, the measured peak [Ca²⁺]_i did not reach levels that are typically required by synaptotagmins to elicit the release of insulin-containing vesicles in live-cell experiments. We therefore sought to resolve the Ca²⁺ dynamics in the subplasmalemmal microdomain that is critical for triggering fast exocytosis. Applying total internal reflection fluorescence (TIRF) microscopy in insulin-producing INS-1E and primary mouse β cells, we resolved extraordinary fast trains of Ca²⁺ spiking (frequency > 3 s^-1) in response to glucose exposure. Using a low-affinity [Ca²⁺]_i indicator dye, we provide experimental evidence that Ca²⁺ spikes reach low micromolar apparent concentrations in the vicinity of the plasma membrane_. Analysis of Ca²⁺ spikes evoked by repeated depolarization for 10 ms closely matched the Ca²⁺ dynamics observed upon glucose application. To our knowledge, this is the first study that experimentally demonstrates Ca²⁺ spikes in β cells with velocities that resemble those of bursting or continuously appearing trains of action potentials (APs) in non-patched cells.

Glucose-dependent activation, activity, and deactivation of beta cell networks in acute mouse pancreas tissue slices

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.

Reduced synchroneity of intra-islet Ca2+ oscillations in vivo in Robo-deficient β cells

The spatial architecture of the islets of Langerhans is hypothesized to facilitate synchronized insulin secretion among β cells, yet testing this in vivo in the intact pancreas is challenging. Robo βKO mice, in which the genes Robo1 and Robo2 are deleted selectively in β cells, provide a unique model of altered islet spatial architecture without loss of β cell differentiation or islet damage from diabetes. Combining Robo βKO mice with intravital microscopy, we show here that Robo βKO islets have reduced synchronized intra-islet Ca2+ oscillations among β cells in vivo. We provide evidence that this loss is not due to a β cell-intrinsic function of Robo, mis-expression or mis-localization of Cx36 gap junctions, or changes in islet vascularization or innervation, suggesting that the islet architecture itself is required for synchronized Ca2+ oscillations. These results have implications for understanding structure-function relationships in the islets during progression to diabetes as well as engineering islets from stem cells.

The Role of cAMP in Beta Cell Stimulus-Secretion and Intercellular Coupling

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 Ca2+ 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.

Intercellular Communication in the Islet of Langerhans in Health and Disease

Blood glucose homeostasis requires proper function of pancreatic islets, which secrete insulin, glucagon, and somatostatin from the β-, α-, and δ-cells, respectively. Each islet cell type is equipped with intrinsic mechanisms for glucose sensing and secretory actions, but these intrinsic mechanisms alone cannot explain the observed secretory profiles from intact islets. Regulation of secretion involves interconnected mechanisms among and between islet cell types. Islet cells lose their normal functional signatures and secretory behaviors upon dispersal as compared to intact islets and in vivo. In dispersed islet cells, the glucose response of insulin secretion is attenuated from that seen from whole islets, coordinated oscillations in membrane potential and intracellular Ca2+ activity, as well as the two-phase insulin secretion profile, are missing, and glucagon secretion displays higher basal secretion profile and a reverse glucose-dependent response from that of intact islets. These observations highlight the critical roles of intercellular communication within the pancreatic islet, and how these communication pathways are crucial for proper hormonal and nonhormonal secretion and glucose homeostasis. Further, misregulated secretions of islet secretory products that arise from defective intercellular islet communication are implicated in diabetes. Intercellular communication within the islet environment comprises multiple mechanisms, including electrical synapses from gap junctional coupling, paracrine interactions among neighboring cells, and direct cell-to-cell contacts in the form of juxtacrine signaling. In this article, we describe the various mechanisms that contribute to proper islet function for each islet cell type and how intercellular islet communications are coordinated among the same and different islet cell types. © 2021 American Physiological Society. Compr Physiol 11:2191-2225, 2021.

Hubs, diversity, and synchronization in FitzHugh-Nagumo oscillator networks: Resonance effects and biophysical implications

Using the FitzHugh-Nagumo equations to represent the oscillatory electrical behavior of β-cells, we develop a coupled oscillator network model with cubic lattice topology, showing that the emergence of pacemakers or hubs in the system can be viewed as a natural consequence of oscillator population diversity. The optimal hub to nonhub ratio is determined by the position of the diversity-induced resonance maximum for a given set of FitzHugh-Nagumo equation parameters and is predicted by the model to be in a range that is fully consistent with experimental observations. The model also suggests that hubs in a β-cell network should have the ability to "switch on" and "off" their pacemaker function. As a consequence, their relative amount in the population can vary in order to ensure an optimal oscillatory performance of the network in response to environmental changes, such as variations of an external stimulus.

The effects of beta-cell mass and function, intercellular coupling, and islet synchrony on [Formula: see text] dynamics

Type 2 diabetes (T2D) is a challenging metabolic disorder characterized by a substantial loss of [Formula: see text]-cell mass and alteration of [Formula: see text]-cell function in the islets of Langerhans, disrupting insulin secretion and glucose homeostasis. The mechanisms for deficiency in [Formula: see text]-cell mass and function during the hyperglycemia development and T2D pathogenesis are complex. To study the relative contribution of [Formula: see text]-cell mass to [Formula: see text]-cell function in T2D, we make use of a comprehensive electrophysiological model of human [Formula: see text]-cell clusters. We find that defect in [Formula: see text]-cell mass causes a functional decline in single [Formula: see text]-cell, impairment in intra-islet synchrony, and changes in the form of oscillatory patterns of membrane potential and intracellular [Formula: see text] concentration, which can lead to changes in insulin secretion dynamics and in insulin levels. The model demonstrates a good correspondence between suppression of synchronizing electrical activity and published experimental measurements. We then compare the role of gap junction-mediated electrical coupling with both [Formula: see text]-cell synchronization and metabolic coupling in the behavior of [Formula: see text] concentration dynamics within human islets. Our results indicate that inter-[Formula: see text]-cellular electrical coupling depicts a more important factor in shaping the physiological regulation of islet function and in human T2D. We further predict that varying the whole-cell conductance of delayed rectifier [Formula: see text] channels modifies oscillatory activity patterns of [Formula: see text]-cell population lacking intercellular coupling, which significantly affect [Formula: see text] concentration and insulin secretion.

Small subpopulations of β-cells do not drive islet oscillatory [Ca2+] dynamics via gap junction communication

The islets of Langerhans exist as multicellular networks that regulate blood glucose levels. The majority of cells in the islet are excitable, insulin-producing β-cells that are electrically coupled via gap junction channels. β-cells are known to display heterogeneous functionality. However, due to gap junction coupling, β-cells show coordinated [Ca²⁺] oscillations when stimulated with glucose, and global quiescence when unstimulated. Small subpopulations of highly functional β-cells have been suggested to control [Ca²⁺] dynamics across the islet. When these populations were targeted by optogenetic silencing or photoablation, [Ca²⁺] dynamics across the islet were largely disrupted. In this study, we investigated the theoretical basis of these experiments and how small populations can disproportionality control islet [Ca²⁺] dynamics. Using a multicellular islet model, we generated normal, skewed or bimodal distributions of β-cell heterogeneity. We examined how islet [Ca²⁺] dynamics were disrupted when cells were targeted via hyperpolarization or populations were removed; to mimic optogenetic silencing or photoablation, respectively. Targeted cell populations were chosen based on characteristics linked to functional subpopulation, including metabolic rate of glucose oxidation or [Ca²⁺] oscillation frequency. Islets were susceptible to marked suppression of [Ca²⁺] when ~10% of cells with high metabolic activity were hyperpolarized; where hyperpolarizing cells with normal metabolic activity had little effect. However, when highly metabolic cells were removed from the model, [Ca²⁺] oscillations remained. Similarly, when ~10% of cells with either the highest frequency or earliest elevations in [Ca²⁺] were removed from the islet, the [Ca²⁺] oscillation frequency remained largely unchanged. Overall, these results indicate small populations of β-cells with either increased metabolic activity or increased frequency are unable to disproportionately control islet-wide [Ca²⁺] via gap junction coupling. Therefore, we need to reconsider the physiological basis for such small β-cell populations or the mechanism by which they may be acting to control normal islet function.

Many biological systems can be studied using network theory. How heterogeneous cell subpopulations come together to create complex multicellular behavior is of great value in understanding function and dysfunction in tissues. The pancreatic islet of Langerhans is a highly coupled structure that is important for maintaining blood glucose homeostasis. β-cell electrical activity is coordinated via gap junction communication. The function of the insulin-producing β-cell within the islet is disrupted in diabetes. As such, to understand the causes of islet dysfunction we need to understand how different cells within the islet contribute to its overall function via gap junction coupling. Using a computational model of β-cell electrophysiology, we investigated how small highly functional β-cell populations within the islet contribute to its function. We found that when small populations with greater functionality were introduced into the islet, they displayed signatures of this enhanced functionality. However, when these cells were removed, the islet, retained near-normal function. Thus, in a highly coupled system, such as an islet, the heterogeneity of cells allows small subpopulations to be dispensable, and thus their absence is unable to disrupt the larger cellular network. These findings can be applied to other electrical systems that have heterogeneous cell populations.

NMDA receptor inhibition increases, synchronizes, and stabilizes the collective pancreatic beta cell activity: Insights through multilayer network analysis

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.

Flipping the switch on the hub cell: Islet desynchronization through cell silencing

Pancreatic β cells, responsible for secreting insulin into the bloodstream and maintaining glucose homeostasis, are organized in the islets of Langerhans as clusters of electrically coupled cells. Gap junctions, connecting neighboring cells, coordinate the behavior of the islet, leading to the synchronized oscillations in the intracellular calcium and insulin secretion in healthy islets. Recent experimental work has shown that silencing special hub cells can lead to a disruption in the coordinated behavior, calling into question the democratic paradigm of islet insulin secretion with more or less equal input from each β cell. Islets were shown to have scale-free functional connectivity and a hub cell whose silencing would lead to a loss of functional connectivity and activity in the islet. A mechanistic model representing the electrical and calcium dynamics of β cells during insulin secretion was applied to a network of cells connected by gap junctions to test the hypothesis of hub cells. Functional connectivity networks were built from the simulated calcium traces, with some networks classified as scale-free, confirming experimental results. Potential hub cells were identified using previously defined centrality measures, but silencing them was unable to desynchronize the islet. Instead, switch cells, which were able to turn off the activity of the islet but were not highly functionally connected, were found via systematically silencing each cell in the network.

Assessing Different Temporal Scales of Calcium Dynamics in Networks of Beta Cell Populations

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 Ca2+ 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.

A Practical Guide to Rodent Islet Isolation and Assessment Revisited

Insufficient insulin secretion is a key component of both type 1 and type 2 diabetes. Since insulin is released by the islets of Langerhans, obtaining viable and functional islets is critical for research and transplantation. The effective and efficient isolation of these small islands of endocrine cells from the sea of exocrine tissue that is the rest of the pancreas is not necessarily simple or quick. Choosing and administering the digestive enzyme, separation of the islets from acinar tissue, and culture of islets are all things that must be considered. The purpose of this review is to provide a history of the development of islet isolation procedures and to serve as a practical guide to rodent islet research for newcomers to islet biology. We discuss key elements of mouse islet isolation including choosing collagenase, the digestion process, purification of islets using a density gradient, and islet culture conditions. In addition, this paper reviews techniques for assessing islet viability and function such as visual assessment, glucose-stimulated insulin secretion and intracellular calcium measurements. A detailed protocol is provided that describes a common method our laboratory uses to obtain viable and functional mouse islets for in vitro study. This review thus provides a strong foundation for successful procurement and purification of high-quality mouse islets for research purposes.

Importance of Both Imprinted Genes and Functional Heterogeneity in Pancreatic Beta Cells: Is There a Link?

Diabetes mellitus now affects more than 400 million individuals worldwide, with significant impacts on the lives of those affected and associated socio-economic costs. Although defects in insulin secretion underlie all forms of the disease, the molecular mechanisms which drive them are still poorly understood. Subsets of specialised beta cells have, in recent years, been suggested to play critical roles in "pacing" overall islet activity. The molecular nature of these cells, the means through which their identity is established and the changes which may contribute to their functional demise and "loss of influence" in both type 1 and type 2 diabetes are largely unknown. Genomic imprinting involves the selective silencing of one of the two parental alleles through DNA methylation and modified imprinted gene expression is involved in a number of diseases. Loss of expression, or loss of imprinting, can be shown in mouse models to lead to defects in beta cell function and abnormal insulin secretion. In the present review we survey the evidence that altered expression of imprinted genes contribute to loss of beta cell function, the importance of beta cell heterogeneity in normal and disease states, and hypothesise whether there is a direct link between the two.

β Cell dysfunction during progression of metabolic syndrome to type 2 diabetes

In a society where physical activity is limited and food supply is abundant, metabolic diseases are becoming a serious epidemic. Metabolic syndrome (MetS) represents a cluster of metabolically related symptoms such as obesity, hypertension, dyslipidemia, and carbohydrate intolerance, and significantly increases type 2 diabetes mellitus risk. Insulin resistance and hyperinsulinemia are consistent characteristics of MetS, but which of these features is the initiating insult is still widely debated. Regardless, both of these conditions trigger adverse responses from the pancreatic β cell, which is responsible for producing, storing, and releasing insulin to maintain glucose homeostasis. The observation that the degree of β cell dysfunction correlates with the severity of MetS highlights the need to better understand β cell dysfunction in the development of MetS. This Review focuses on the current understanding from rodent and human studies of the progression of β cell responses during the development of MetS, as well as recent findings addressing the complexity of β cell identity and heterogeneity within the islet during disease progression. The differential responses observed in β cells together with the heterogeneity in disease phenotypes within the patient population emphasize the need to better understand the mechanisms behind β cell adaptation, identity, and dysfunction in MetS.

Rapid lipolytic oscillations in ex vivo adipose tissue explants revealed through microfluidic droplet sampling at high temporal resolution

Our understanding of adipose tissue biology has steadily evolved. While structural and energy storage functionalities have been in the forefront, a key endocrine role for adipocytes was revealed only over the last few decades. In contrast to the wealth of information on dynamic function of other endocrine tissues, few studies have focused on dynamic adipose tissue function or on tool development toward that end. Here, we apply our unique droplet-based microfluidic devices to culture, perfuse, and sample secretions from primary murine epididymal white adipose tissue (eWAT), and from predifferentiated clusters of 3T3-L1 adipocytes. Through automated control, oil-segmented aqueous droplets (∼2.6 nL) were sampled from tissue or cells at 3.5 second temporal resolution (including sample and reference droplets), with integrated enzyme assays enabling real-time quantification of glycerol (down to 1.9 fmol per droplet). This high resolution revealed previously unreported oscillations in secreted glycerol at frequencies of 0.2 to 2.0 min^-1 (∼30-300 s periods) present in the primary tissue but not in clustered cells. Low-level bursts (∼50 fmol) released in basal conditions were contrasted with larger bursts (∼300 fmol) during stimulation. Further, both fold changes and burst magnitudes were decreased in eWAT of aged and obese mice. These results, combined with immunostaining and photobleaching analyses, suggest that gap-junctional coupling or nerve cell innervation within the intact ex vivo tissue explants play important roles in this apparent tissue-level, lipolytic synchronization. High-resolution, quantitative sampling by droplet microfluidics thus permitted unique biological information to be observed, giving an analytical framework poised for future studies of dynamic oscillatory function of adipose and other tissues.

Heterogeneous alpha-cell population modeling of glucose-induced inhibition of electrical activity

Glucagon release from the pancreatic alpha-cells is regulated by glucose, but the underlying mechanisms are far from understood. It is known that the alpha-cell population is very heterogeneous, but - compared to the insulin-secreting beta-cells - the consequences of this cell-to-cell variation are much less studied. Since the alpha-cells are not electrically coupled, large differences in the single cell responses are to be expected, and this variation may contribute to the confusion regarding the mechanisms of glucose-induced suppression of glucagon release. Using mathematical modeling of alpha-cells with realistic cell-to-cell parameter variation based on recent experimental results, we show that the simulated alpha-cells exhibit great diversity in their electrophysiological behavior. To robustly reproduce experimental recordings from alpha-cell exposed to a rise in glucose levels, we must assume that both intrinsic mechanisms and paracrine signals contribute to glucose-induced changes in electrical activity. Our simulations suggest that the sum of different electrophysiological responses due to alpha-cell heterogeneity is involved in glucose-suppressed glucagon secretion, and that more than one mechanism contribute to control the alpha-cell populations' behavior. Finally, we apply regression analysis to our synthetic alpha-cell population to infer which membrane currents influence electrical activity in alpha-cells at different glucose levels. The results from such statistical modeling suggest possible disturbances underlying defect regulation of alpha-cell electrical behavior in diabetics. Thus, although alpha-cells appear to be inherently complex and heterogeneous as reflected in published data, realistic modeling of the alpha-cells at the population level provides insight into the mechanisms of glucagon release.

Metabolic and Functional Heterogeneity in Pancreatic β Cells

Metabolic and secretory heterogeneity are fundamental properties of pancreatic islet β cells. Emerging data suggest that stable differences in the transcriptome and proteome of individual cells may create cellular hierarchies, which, in turn, establish coordinated functional networks. These networks appear to govern the secretory activity of the whole islet and be affected in some forms of diabetes mellitus. Functional imaging, for example, of intracellular calcium dynamics, has led to the demonstration of "small worlds" behavior, and the identification of highly connected "hub" (or "leader") cells and of follower populations subservient to them. Subsequent inactivation of members of either population, for example, using optogenetic approaches or photoablation, has confirmed the importance of hub cells as possible pacemakers. Hub cells appear to be enriched for the glucose phosphorylating enzyme glucokinase and for genes encoding other enzymes involved in glucose metabolism compared to follower cells. Recent findings have shown the relevance of cellular hierarchy in islets from multiple species including human, mouse and fish, and shown that it is preserved in vivo in the context of the fully vascularized and innervated islet. Importantly, connectivity is impaired by insults, which mimic the diabetic milieu, including high glucose and/or fatty levels, and by the ablation of genes associated with type 2 diabetes risk in genome-wide association studies. We discuss here the evidence for the existence of these networks and their failure in disease settings. We also briefly survey the challenges in understanding their properties.

Biophysical modeling of β-cells networks: Realistic architectures and heterogeneity effects

The β-cells dynamics is the regulator of insulin secretion in the pancreas, and its investigation is a central aspect in designing effective treatment strategies for diabetes. Despite great efforts, much is still unknown about the complex organization of such endocrine cells and realistic mathematical modeling represents a useful tool to elucidate key aspects of glucose control in humans. In this contribution, we study the human β-cells collective behaviour, by modeling their electric and metabolic coupling in a cluster, of size and architecture similar to human islets of Langerhans. We focus on the effect of coupling on various dynamics regimes observed in the islets, that are spiking and bursting on multiple timescales. In particular, we test the effect of hubs, that are highly glucose-sensitive β-cells, on the overall network dynamics, observing different modulation depending on the timescale of the dynamics. By properly taking into account the role of cells heterogeneity, recently emerged, our model effectively describes the effect of hubs on the synchronization of the islet response and the correlation of β-cells activity.

Heterogeneity and Delayed Activation as Hallmarks of Self-Organization and Criticality in Excitable Tissue

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.

Beta-cell excitability and excitability-driven diabetes in adult Zebrafish islets

Islet β-cell membrane excitability is a well-established regulator of mammalian insulin secretion, and defects in β-cell excitability are linked to multiple forms of diabetes. Evolutionary conservation of islet excitability in lower organisms is largely unexplored. Here we show that adult zebrafish islet calcium levels rise in response to elevated extracellular [glucose], with similar concentration-response relationship to mammalian β-cells. However, zebrafish islet calcium transients are nor well coupled, with a shallower glucose-dependence of cytoplasmic calcium concentration. We have also generated transgenic zebrafish that conditionally express gain-of-function mutations in ATP-sensitive K+ channels (KATP -GOF) in β-cells. Following induction, these fish become profoundly diabetic, paralleling features of mammalian diabetes resulting from equivalent mutations. KATP -GOF fish become severely hyperglycemic, with slowed growth, and their islets lose glucose-induced calcium responses. These results indicate that, although lacking tight cell-cell coupling of intracellular Ca2+ , adult zebrafish islets recapitulate similar excitability-driven β-cell glucose responsiveness to those in mammals, and exhibit profound susceptibility to diabetes as a result of inexcitability. While illustrating evolutionary conservation of islet excitability in lower vertebrates, these results also provide important validation of zebrafish as a suitable animal model in which to identify modulators of islet excitability and diabetes.

Comparing methods to normalize insulin secretion shows the process may not be needed

Glucose-stimulated insulin secretion (GSIS) is a well-accepted method to investigate the physiological and pathophysiological function of islets. However, there is little consensus about which method is best for normalizing and presenting GSIS data. In this study, we evaluated the sufficiency of islet area, total protein, total DNA, and total insulin content as parameters to normalize GSIS data. First, we tested if there is a linear correlation between each parameter and the number of islets (10, 20, 30, and 40 islets). Islet area, total protein, and insulin content produced excellent linear correlations with islet number (R2 >0.9 for each) from the same islet material. Insulin secretion in 11mM glucose also correlated reasonably well for islet area (R2=0.69), protein (R2=0.49), and insulin content (R2=0.58). DNA content was difficult to reliably measure and was excluded from additional comparisons. We next measured GSIS for 18 replicates of 20 islets each, measuring 3mM and 11mM glucose to calculate the stimulation index and to compare each normalization parameter. Using these similar islet masses for each replicate, none of the parameters produced linear correlations with GSIS (R2<0.05), suggesting that inherent differences in GSIS dominate small differences in islet mass. We conclude that when comparing GSIS for islets of reasonably similar size (<50% variance), normalization does not improve the representation of GSIS data. Normalization may be beneficial when substantial differences in islet mass are involved. In such situations, we suggest that using islet cross-sectional area is superior to other commonly used techniques for normalizing GSIS data.

Transcriptional Heterogeneity of Beta Cells in the Intact Pancreas

Pancreatic beta cells have been shown to be heterogeneous at multiple levels. However, spatially interrogating transcriptional heterogeneity in the intact tissue has been challenging. Here, we developed an optimized protocol for single-molecule transcript imaging in the intact pancreas and used it to identify a sub-population of "extreme" beta cells with elevated mRNA levels of insulin and other secretory genes. Extreme beta cells contain higher ribosomal and proinsulin content but lower levels of insulin protein in fasted states, suggesting they may be tuned for basal insulin secretion. They exhibit a distinctive intra-cellular polarization pattern, with elevated mRNA concentrations in an apical ER-enriched compartment, distinct from the localization of nascent and mature proteins. The proportion of extreme cells increases in db/db diabetic mice, potentially facilitating the required increase in basal insulin. Our results thus highlight a sub-population of beta cells that may carry distinct functional roles along physiological and pathological timescales.

The Difference δ-Cells Make in Glucose Control

The role of beta and α-cells to glucose control are established, but the physiological role of δ-cells is poorly understood. Delta-cells are ideally positioned within pancreatic islets to modulate insulin and glucagon secretion at their source. We review the evidence for a negative feedback loop between delta and β-cells that determines the blood glucose set point and suggest that local δ-cell-mediated feedback stabilizes glycemic control.

Regulation of islet glucagon secretion: Beyond calcium

The islet of Langerhans plays a key role in glucose homeostasis through regulated secretion of the hormones insulin and glucagon. Islet research has focused on the insulin-secreting β-cells, even though aberrant glucagon secretion from α-cells also contributes to the aetiology of diabetes. Despite its importance, the mechanisms controlling glucagon secretion remain controversial. Proper α-cell function requires the islet milieu, where β- and δ-cells drive and constrain α-cell dynamics. The response of glucagon to glucose is similar between isolated islets and that measured in vivo, so it appears that the glucose dependence requires only islet-intrinsic factors and not input from blood flow or the nervous system. Elevated intracellular free Ca²⁺ is needed for α-cell exocytosis, but interpreting Ca²⁺ data is tricky since it is heterogeneous among α-cells at all physiological glucose levels. Total Ca²⁺ activity in α-cells increases slightly with glucose, so Ca²⁺ may serve a permissive, rather than regulatory, role in glucagon secretion. On the other hand, cAMP is a more promising candidate for controlling glucagon secretion and is itself driven by paracrine signalling from β- and δ-cells. Another pathway, juxtacrine signalling through the α-cell EphA receptors, stimulated by β-cell ephrin ligands, leads to a tonic inhibition of glucagon secretion. We discuss potential combinations of Ca²⁺ , cAMP, paracrine and juxtacrine factors in the regulation of glucagon secretion, focusing on recent data in the literature that might unify the field towards a quantitative understanding of α-cell function.

The role of beta cell heterogeneity in islet function and insulin release

It is becoming increasingly apparent that not all insulin-secreting beta cells are equal. Subtle differences exist at the transcriptomic and protein expression levels, with repercussions for beta cell survival/proliferation, calcium signalling and insulin release. Notably, beta cell heterogeneity displays plasticity during development, metabolic stress and type 2 diabetes mellitus (T2DM). Thus, heterogeneity or lack thereof may be an important contributor to beta cell failure during T2DM in both rodents and humans. The present review will discuss the molecular and cellular features of beta cell heterogeneity at both the single-cell and islet level, explore how this influences islet function and insulin release and look into the alterations that may occur during obesity and T2DM.

Imaging endocrinology in animal models of endocrine disease

Endocrine organs secrete a variety of hormones involved in the regulation of a multitude of body functions. Although pancreatic islets were discovered at the turn of the 19th century, other endocrine glands remained commonly described as diffuse endocrine systems. Over the last two decades, development of new imaging techniques and genetically-modified animals with cell-specific fluorescent tags or specific hormone deficiencies have enabled in vivo imaging of endocrine organs and revealed intricate endocrine cell network structures and plasticity. Overall, these new tools have revolutionized our understanding of endocrine function. The overarching aim of this Review is to describe the current mechanistic understanding that has emerged from imaging studies of endocrine cell network structure/function relationships in animal models, with a particular emphasis on the pituitary gland and the endocrine pancreas.