An algorithm for automated detection,localization and measurement of local calcium signals from camera-based imaging

Local Ca²⁺ transients such as puffs and sparks form the building blocks of cellular Ca²⁺ signaling in numerous cell types. They have traditionally been studied by linescan confocal microscopy, but advances in TIRF microscopy together with improved electron-multiplied CCD (EMCCD) cameras now enable rapid (>500 frames s⁻¹) imaging of subcellular Ca²⁺ signals with high spatial resolution in two dimensions. This approach yields vastly more information (ca. 1 Gb min⁻¹) than linescan imaging, rendering visual identification and analysis of local events imaged both laborious and subject to user bias. Here we describe a routine to rapidly automate identification and analysis of local Ca²⁺ events. This features an intuitive graphical user-interfaces and runs under Matlab and the open-source Python software. The underlying algorithm features spatial and temporal noise filtering to reliably detect even small events in the presence of noisy and fluctuating baselines; localizes sites of Ca²⁺ release with sub-pixel resolution; facilitates user review and editing of data; and outputs time-sequences of fluorescence ratio signals for identified event sites along with Excel-compatible tables listing amplitudes and kinetics of events.     

Imaging single-channel calcium microdomains

The Ca²⁺ microdomains generated around the mouth of open ion channels represent the basic building blocks from which cytosolic Ca²⁺ signals are constructed. Recent improvements in optical imaging techniques now allow these microdomains to be visualized as single channel calcium fluorescence transients (SCCaFTs), providing information about channel properties that was previously accessible only by electrophysiological patch-clamp recordings. We review recent advances in single channel Ca²⁺ imaging methodologies, with emphasis on total internal reflection fluorescence microscopy (TIRFM) as the technique of choice for recording SCCaFTs from voltage- and ligand-gated plasmalemmal ion channels. This technique of ‘optical patch-clamp recording’ is massively parallel, permitting simultaneous imaging of hundreds of channels; provides millisecond resolution of gating kinetics together with sub-micron spatial resolution of channel locations; and is applicable to diverse families of membrane channels that display partial permeability to Ca²⁺ ions.     

Superresolution Localization of Single Functional IP3R Channels Utilizing Ca Flux as a Readout

The subcellular localization of membrane Ca²⁺ channels is crucial for their functioning, but is difficult to study because channels may be distributed more closely than the resolution of conventional microscopy is able to detect. We describe a technique, stochastic channel Ca²⁺ nanoscale resolution (SCCaNR), employing Ca²⁺-sensitive fluorescent dyes to localize stochastic openings and closings of single Ca²⁺-permeable channels within <50 nm, and apply it to examine the clustered arrangement of inositol trisphosphate receptor (IP₃R) channels underlying local Ca²⁺ puffs. Fluorescence signals (blips) arising from single functional IP₃Rs are almost immotile (diffusion coefficient <0.003 μm² s⁻¹), as are puff sites over prolonged periods, suggesting that the architecture of this signaling system is stable and not subject to rapid, dynamic rearrangement. However, rapid stepwise changes in centroid position of fluorescence are evident within the durations of individual puffs. These apparent movements likely result from asynchronous gating of IP₃Rs distributed within clusters that have an overall diameter of ∼400 nm, indicating that the nanoscale architecture of IP₃R clusters is important in shaping local Ca²⁺ signals. We anticipate that SCCaNR will complement superresolution techniques such as PALM and STORM for studies of Ca²⁺ channels as it obviates the need for photoswitchable labels and provides functional as well as spatial information.     

Noise analysis of cytosolic calcium image data

Cellular Ca²⁺ signals are often constrained to cytosolic micro- or nano-domains where stochastic openings of Ca²⁺ channels cause large fluctuations in local Ca²⁺ concentration (Ca²⁺ ‘noise’). With the advent of TIRF microscopy to image the fluorescence of Ca²⁺-sensitive probes from attoliter volumes it has become possible to directly monitor these signals, which closely track the gating of plasmalemmal and ER Ca²⁺-permeable channels. Nevertheless, it is likely that many physiologically important Ca²⁺ signals are too small to resolve as discrete events in fluorescence recordings. By analogy with noise analysis of electrophysiological data, we explore here the use of statistical approaches to detect and analyze such Ca²⁺ noise in images obtained using Ca²⁺-sensitive indicator dyes. We describe two techniques - power spectrum analysis and spatio-temporal correlation - and demonstrate that both effectively identify discrete, spatially localized calcium release events (Ca²⁺ puffs). Moreover, we show they are able to detect localized noise fluctuations in a case where discrete events cannot directly be resolved.     

Images of Ca flux in astrocytes: evidence for spatially distinct sites of Ca release and uptake

In this study, we have developed a mathematical method to derive the Ca²⁺ fluxes underlying agonist-evoked Ca²⁺ waves in cultured rat cortical astrocytes. Astrocytes were stimulated with norepinephrine (100 nM) to evoke Ca²⁺ waves, which were recorded by measuring FIuo-3 fluorescence changes with high spatial and temporal resolution. Normalized fluorescence (ΔF/F) was analyzed in discrete cellular spaces in a series of successive slices along the length of the cell. From these data, Ca²⁺ flux was then calculated using a one dimensional reaction-diffusion equation which utilizes the temporal and spatial derivatives of the fluorescence data and the diffusion coefficient of Ca²⁺ in the cytosol. This method identified distinct sites of positive flux (Ca²⁺ release into the cytosol) and of negative flux (Ca²⁺ removal from cytosol) and showed that in astrocytes, sites of Ca²⁺ release from stores regularly alternate with sites of Ca²⁺ removal from the cytosol. Cross correlation analysis of the two distribution patterns gave positive correlation at 2 μm out of phase and a negative correlation in phase. Thapsigargin-induced Ca²⁺ waves were analyzed to determine if the negative flux was due to Ca²⁺ uptake via thapsigargin-sensitive Ca²⁺ pumps. Negative flux sites were still found under these conditions, suggesting that multiple mechanisms of Ca²⁺ removal from the cytosol may contribute to negative flux sites. This method of calculation of flux may serve as a means to describe the distribution of functional ion channels and pumps participating in cellular Ca²⁺ signalling.     

Dynamic Ca2+ imaging with a simplified lattice light-sheet microscope: A sideways view of subcellular Ca2+ puffs

We describe the construction of a simplified, inexpensive lattice light-sheet microscope, and illustrate its use for imaging subcellular Ca²⁺ puffs evoked by photoreleased i-IP₃ in cultured SH-SY5Y neuroblastoma cells loaded with the Ca²⁺ probe Cal520. The microscope provides sub-micron spatial resolution and enables recording of local Ca²⁺ transients in single-slice mode with a signal-to-noise ratio and temporal resolution (2 ms) at least as good as confocal or total internal reflection microscopy. Signals arising from openings of individual IP₃R channels are clearly resolved, as are stepwise changes in fluorescence reflecting openings and closings of individual channels during puffs. Moreover, by stepping the specimen through the light-sheet, the entire volume of a cell can be scanned within a few hundred ms. The ability to directly visualize a sideways (axial) section through cells directly reveals that IP₃-evoked Ca²⁺ puffs originate at sites in very close (≤a few hundred nm) to the plasma membrane, suggesting they play a specific role in signaling to the membrane.     

Thermal modulation of epicardial Ca2+ dynamics uncovers molecular mechanisms of Ca2+ alternans

Ca²⁺ alternans (Ca-Alts) are alternating beat-to-beat changes in the amplitude of Ca²⁺ transients that frequently occur during tachycardia, ischemia, or hypothermia that can lead to sudden cardiac death. Ca-Alts appear to result from a variation in the amount of Ca²⁺ released from the sarcoplasmic reticulum (SR) between two consecutive heartbeats. This variable Ca²⁺ release has been attributed to the alternation of the action potential duration, delay in the recovery from inactivation of RYR Ca²⁺ release channel (RYR2), or an incomplete Ca²⁺ refilling of the SR. In all three cases, the RYR2 mobilizes less Ca²⁺ from the SR in an alternating manner, thereby generating an alternating profile of the Ca²⁺ transients. We used a new experimental approach, fluorescence local field optical mapping (FLOM), to record at the epicardial layer of an intact heart with subcellular resolution. In conjunction with a local cold finger, a series of images were recorded within an area where the local cooling induced a temperature gradient. Ca-Alts were larger in colder regions and occurred without changes in action potential duration. Analysis of the change in the enthalpy and Q₁₀ of several kinetic processes defining intracellular Ca²⁺ dynamics indicated that the effects of temperature change on the relaxation of intracellular Ca²⁺ transients involved both passive and active mechanisms. The steep temperature dependency of Ca-Alts during tachycardia suggests Ca-Alts are generated by insufficient SERCA-mediated Ca²⁺ uptake into the SR. We found that Ca-Alts are heavily dependent on intra-SR Ca²⁺ and can be promoted through partial pharmacologic inhibition of SERCA2a. Finally, the FLOM experimental approach has the potential to help us understand how arrhythmogenesis correlates with the spatial distribution of metabolically impaired myocytes along the myocardium.     

Voltage-Activated Elementary Calcium Release Events in Isolated Mouse Skeletal Muscle Fibers

The elementary Ca²⁺-release events underlying voltage-activated myoplasmic Ca²⁺ transients in mammalian muscle remain elusive. Here, we looked for such events in confocal line-scan (x,t) images of fluo-3 fluorescence taken from isolated adult mouse skeletal muscle fibers held under voltage-clamp conditions. In response to step depolarizations, spatially segregated fluorescence signals could be detected that were riding on a global increase in fluorescence. These discrete signals were separated using digital filtering in the spatial domain; mean values for their spatial half-width and amplitude were 1.99 ± 0.09 μm and 0.16 ± 0.005 ΔF/F ₀ (n = 151), respectively. Under control conditions, the duration of the events was limited by the pulse duration. In contrast, in the presence of maurocalcine, a scorpion toxin suspected to disrupt the process of repolarization-induced ryanodine receptor (RyR) closure, events uninterrupted by the end of the pulse were readily detected. Overall results establish these voltage-activated low-amplitude local Ca²⁺ signals as inherent components of the physiological Ca²⁺-release process of mammalian muscle and suggest that they result from the opening of either one RyR or a coherently operating group of RyRs, under the control of the plasma membrane polarization.     

Automated Detection and Analysis of Ca2+ Sparks in x-y Image Stacks Using a Thresholding Algorithm Implemented within the Open-Source Image Analysis Platform ImageJ

Previous studies have used analysis of Ca²⁺ sparks extensively to investigate both normal and pathological Ca²⁺ regulation in cardiac myocytes. The great majority of these studies used line-scan confocal imaging. In part, this is because the development of open-source software for automatic detection of Ca²⁺ sparks in line-scan images has greatly simplified data analysis. A disadvantage of line-scan imaging is that data are collected from a single row of pixels, representing only a small fraction of the cell, and in many instances x-y confocal imaging is preferable. However, the limited availability of software for Ca²⁺ spark analysis in two-dimensional x-y image stacks presents an obstacle to its wider application. This study describes the development and characterization of software to enable automatic detection and analysis of Ca²⁺ sparks within x-y image stacks, implemented as a plugin within the open-source image analysis platform ImageJ. The program includes methods to enable precise identification of cells within confocal fluorescence images, compensation for changes in background fluorescence, and options that allow exclusion of events based on spatial characteristics.     

Automated Detection and Analysis of Ca Sparks in x-y Image Stacks Using a Thresholding Algorithm Implemented within the Open-Source Image Analysis Platform ImageJ

Previous studies have used analysis of Ca²⁺ sparks extensively to investigate both normal and pathological Ca²⁺ regulation in cardiac myocytes. The great majority of these studies used line-scan confocal imaging. In part, this is because the development of open-source software for automatic detection of Ca²⁺ sparks in line-scan images has greatly simplified data analysis. A disadvantage of line-scan imaging is that data are collected from a single row of pixels, representing only a small fraction of the cell, and in many instances x-y confocal imaging is preferable. However, the limited availability of software for Ca²⁺ spark analysis in two-dimensional x-y image stacks presents an obstacle to its wider application. This study describes the development and characterization of software to enable automatic detection and analysis of Ca²⁺ sparks within x-y image stacks, implemented as a plugin within the open-source image analysis platform ImageJ. The program includes methods to enable precise identification of cells within confocal fluorescence images, compensation for changes in background fluorescence, and options that allow exclusion of events based on spatial characteristics.     

Combining Voltage and Calcium Imaging from Neuronal Dendrites

The ability to monitor membrane potential (V _m) and calcium (Ca²⁺) transients at multiple locations on the same neuron can facilitate further progress in our understanding of neuronal function. Here we describe a method to combine V _m and Ca²⁺ imaging using styryl voltage sensitive dyes and Fura type UV-excitable Ca²⁺ indicators. In all cases V _m optical signals are linear with membrane potential changes, but the calibration of optical signals on an absolute scale is presently possible only in some neurons. The interpretation of Ca²⁺ optical signals depends on the indicator Ca²⁺ buffering capacity relative to the cell endogenous buffering capacity. In hippocampal CA1 pyramidal neurons, loaded with JPW-3028 and 300 μM Bis-Fura-2, V _m optical signals cannot be calibrated and the physiological Ca²⁺ dynamics are compromised by the presence of the indicator. Nevertheless, at each individual site, relative changes in V _m and Ca²⁺ fluorescence signals under different conditions can provide meaningful new information on local dendritic integration. In cerebellar Purkinje neurons, loaded with JPW-1114 and 1 mM Fura-FF, V _m optical signals can be calibrated in terms of mV and Ca²⁺ optical signals quantitatively reveal the physiological changes in free Ca²⁺. Using these two examples, the method is explained in detail.     

Distinct Contributions of Orai1 and TRPC1 to Agonist-Induced [Ca2+]i Signals Determine Specificity of Ca2+-Dependent Gene Expression

Regulation of critical cellular functions, including Ca²⁺-dependent gene expression, is determined by the temporal and spatial aspects of agonist-induced Ca²⁺ signals. Stimulation of cells with physiological concentrations of agonists trigger increases [Ca²⁺]_i due to intracellular Ca²⁺ release and Ca²⁺ influx. While Orai1-STIM1 channels account for agonist-stimulated [Ca²⁺]_i increase as well as activation of NFAT in cells such as lymphocytes, RBL and mast cells, both Orai1-STIM1 and TRPC1-STIM1 channels contribute to [Ca²⁺]_i increases in human submandibular gland (HSG) cells. However, only Orai1-mediated Ca²⁺ entry regulates the activation of NFAT in HSG cells. Since both TRPC1 and Orai1 are activated following internal Ca²⁺ store depletion in these cells, it is not clear how the cells decode individual Ca²⁺ signals generated by the two channels for the regulation of specific cellular functions. Here we have examined the contributions of Orai1 and TRPC1 to carbachol (CCh)-induced [Ca²⁺]_i signals and activation of NFAT in single cells. We report that Orai1-mediated Ca²⁺ entry generates [Ca²⁺]_i oscillations at different [CCh], ranging from very low to high. In contrast, TRPC1-mediated Ca²⁺ entry generates sustained [Ca²⁺]_i elevation at high [CCh] and contributes to frequency of [Ca²⁺]_i oscillations at lower [agonist]. More importantly, the two channels are coupled to activation of distinct Ca²⁺ dependent gene expression pathways, consistent with the different patterns of [Ca²⁺]_i signals mediated by them. Nuclear translocation of NFAT and NFAT-dependent gene expression display “all-or-none” activation that is exclusively driven by local [Ca²⁺]_i generated by Orai1, independent of global [Ca²⁺]_i changes or TRPC1-mediated Ca²⁺ entry. In contrast, Ca²⁺ entry via TRPC1 primarily regulates NFκB-mediated gene expression. Together, these findings reveal that Orai1 and TRPC1 mediate distinct local and global Ca²⁺ signals following agonist stimulation of cells, which determine the functional specificity of the channels in activating different Ca²⁺-dependent gene expression pathways.     

Automatic Detection and Classification of Ca2+ Release Events in Line- and Frame-Scan Images

Analysis of Ca²⁺ signals obtained in various cell types (i.e., cardiomyocytes) is always a tradeoff between acquisition speed and signal/noise ratio of the fluorescence signal. This becomes especially apparent during fast two- or three-dimensional confocal imaging when local intracellular fluorescence signals originating from Ca²⁺ release from intracellular Ca²⁺ stores (e.g., sarcoplasmic reticulum) need to be examined. Mathematical methods have been developed to remedy a high noise level by fitting each pixel with a transient function to “denoise” the image. So far, current available analytical approaches have been impaired by a number of constraints (e.g., inability to fit local, concurrent, and consecutive events) and the limited ability to customize implementation. Here, we suggest a, to our knowledge, novel approach for detailed analysis of subcellular micro-Ca²⁺ events based on pixel-by-pixel denoising of confocal frame- and line-scan images. The algorithm enables spatiotemporally overlapping events (e.g., a Ca²⁺ spark occurring during the decaying phase of a Ca²⁺ wave) to be extracted so that various types of Ca²⁺ events can be detected at a pixel time level of precision. The method allows a nonconstant baseline to be estimated for each pixel, foregoing the need to subtract fluorescence background or apply self-ratio methods before image analysis. Furthermore, by using a clustering algorithm, identified single-pixel events are grouped into “physiologically relevant” Ca²⁺ signaling events spanning multiple pixels (sparks, waves, puffs, transients, etc.), from which spatiotemporal event parameters (e.g., full duration at half maximal amplitude, full width at half maximal amplitude, amplitude, wave speed, rise, and decay times) can be easily extracted. The method was implemented with cross-platform open source software, providing a comprehensive and easy-to-use graphical user interface enabling rapid line-scan images and rapid frame-scan image sequences (up to 150 frames/s) to be analyzed and repetitive Ca²⁺ events (Ca²⁺ sparks and Ca²⁺ puffs) originating from clusters of Ca²⁺ release channels located in the sarcoplasmic reticulum membrane (ryanodine receptors and inositol 1,4,5-trisphosphate receptors) of isolated cardiomyocytes to be examined with a high level of precision.     

Examination of the Effects of Heterogeneous Organization of RyR Clusters,Myofibrils and Mitochondria on Ca2+ Release Patterns in Cardiomyocytes

Spatio-temporal dynamics of intracellular calcium, [Ca²⁺]_i, regulate the contractile function of cardiac muscle cells. Measuring [Ca²⁺]_i flux is central to the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease. However, current imaging techniques are limited in the spatial resolution to which changes in [Ca²⁺]_i can be detected. Using spatial point process statistics techniques we developed a novel method to simulate the spatial distribution of RyR clusters, which act as the major mediators of contractile Ca²⁺ release, upon a physiologically-realistic cellular landscape composed of tightly-packed mitochondria and myofibrils. We applied this method to computationally combine confocal-scale (~ 200 nm) data of RyR clusters with 3D electron microscopy data (~ 30 nm) of myofibrils and mitochondria, both collected from adult rat left ventricular myocytes. Using this hybrid-scale spatial model, we simulated reaction-diffusion of [Ca²⁺]_i during the rising phase of the transient (first 30 ms after initiation). At 30 ms, the average peak of the simulated [Ca²⁺]_i transient and of the simulated fluorescence intensity signal, F/F₀, reached values similar to that found in the literature ([Ca²⁺]_i ≈1 μM; F/F₀≈5.5). However, our model predicted the variation in [Ca²⁺]_i to be between 0.3 and 12.7 μM (~3 to 100 fold from resting value of 0.1 μM) and the corresponding F/F₀ signal ranging from 3 to 9.5. We demonstrate in this study that: (i) heterogeneities in the [Ca²⁺]_i transient are due not only to heterogeneous distribution and clustering of mitochondria; (ii) but also to heterogeneous local densities of RyR clusters. Further, we show that: (iii) these structure-induced heterogeneities in [Ca²⁺]_i can appear in line scan data. Finally, using our unique method for generating RyR cluster distributions, we demonstrate the robustness in the [Ca²⁺]_i transient to differences in RyR cluster distributions measured between rat and human cardiomyocytes.     

Factors Determining the Recruitment of Inositol Trisphosphate Receptor Channels During Calcium Puffs

Puffs are localized, transient elevations in cytosolic Ca²⁺ that serve both as the building blocks of global cellular Ca²⁺ signals and as local signals in their own right. They arise from clustered inositol 1,4,5-trisphosphate receptor/channels (IP₃Rs), whose openings are coordinated by Ca²⁺-induced Ca²⁺ release (CICR). We utilized total internal reflection fluorescence imaging of Ca²⁺ signals in neuroblastoma cells with single-channel resolution to elucidate the mechanisms determining the triggering, amplitudes, kinetics, and spatial spread of puffs. We find that any given channel in a cluster has a mean probability of ∼66% of opening following opening of an initial “trigger” channel, and the probability of puff triggering thus increases steeply with increasing number of channels in a cluster (cluster size). Mean puff amplitudes scale with cluster size, but individual amplitudes vary widely, even at sites of similar cluster size, displaying similar proportions of events involving any given number of the channels in the cluster. Stochastic variation in numbers of Ca²⁺-inhibited IP₃Rs likely contributes to the variability of amplitudes of repeated puffs at a site but the amplitudes of successive puffs were uncorrelated, even though we observed statistical correlations between interpuff intervals and puff amplitudes. Initial puffs evoked following photorelease of IP₃—which would not be subject to earlier Ca²⁺-inhibition—also showed wide variability, indicating that mechanisms such as stochastic variation in IP₃ binding and channel recruitment by CICR further determine puff amplitudes. The mean termination time of puffs lengthened with increasing puff amplitude size, consistent with independent closings of channels after a given mean open time, but we found no correlation of termination time with cluster size independent of puff amplitude. The spatial extent of puffs increased with their amplitude, and puffs of similar size were of similar width, independent of cluster size.     

Factors Determining the Recruitment of Inositol Trisphosphate Receptor Channels During Calcium Puffs

Puffs are localized, transient elevations in cytosolic Ca²⁺ that serve both as the building blocks of global cellular Ca²⁺ signals and as local signals in their own right. They arise from clustered inositol 1,4,5-trisphosphate receptor/channels (IP₃Rs), whose openings are coordinated by Ca²⁺-induced Ca²⁺ release (CICR). We utilized total internal reflection fluorescence imaging of Ca²⁺ signals in neuroblastoma cells with single-channel resolution to elucidate the mechanisms determining the triggering, amplitudes, kinetics, and spatial spread of puffs. We find that any given channel in a cluster has a mean probability of ∼66% of opening following opening of an initial “trigger” channel, and the probability of puff triggering thus increases steeply with increasing number of channels in a cluster (cluster size). Mean puff amplitudes scale with cluster size, but individual amplitudes vary widely, even at sites of similar cluster size, displaying similar proportions of events involving any given number of the channels in the cluster. Stochastic variation in numbers of Ca²⁺-inhibited IP₃Rs likely contributes to the variability of amplitudes of repeated puffs at a site but the amplitudes of successive puffs were uncorrelated, even though we observed statistical correlations between interpuff intervals and puff amplitudes. Initial puffs evoked following photorelease of IP₃—which would not be subject to earlier Ca²⁺-inhibition—also showed wide variability, indicating that mechanisms such as stochastic variation in IP₃ binding and channel recruitment by CICR further determine puff amplitudes. The mean termination time of puffs lengthened with increasing puff amplitude size, consistent with independent closings of channels after a given mean open time, but we found no correlation of termination time with cluster size independent of puff amplitude. The spatial extent of puffs increased with their amplitude, and puffs of similar size were of similar width, independent of cluster size.     

Smooth muscle gap-junctions allow propagation of intercellular Ca2+ waves and vasoconstriction due to Ca2+ based action potentials in rat mesenteric resistance arteries

The role of vascular gap junctions in the conduction of intercellular Ca²⁺ and vasoconstriction along small resistance arteries is not entirely understood. Some depolarizing agents trigger conducted vasoconstriction while others only evoke a local depolarization. Here we use a novel technique to investigate the temporal and spatial relationship between intercellular Ca²⁺ signals generated by smooth muscle action potentials (APs) and vasoconstriction in mesenteric resistance arteries (MA). Pulses of exogenous KCl to depolarize the downstream end (T1) of a 3 mm long artery increased intracellular Ca²⁺ associated with vasoconstriction. The spatial spread and amplitude of both depended on the duration of the pulse, with only a restricted non-conducting vasoconstriction to a 1 s pulse. While blocking smooth muscle cell (SMC) K⁺ channels with TEA and activating L-type voltage-gated Ca²⁺ channels (VGCCs) with BayK 8644 spread was dramatically facilitated, so the 1 s pulse evoked intercellular Ca²⁺ waves and vasoconstriction that spread along an entire artery segment 3000 μm long. Ca²⁺ waves spread as nifedipine-sensitive Ca²⁺ spikes due to SMC action potentials, and evoked vasoconstriction. Both intercellular Ca²⁺ and vasoconstriction spread at circa 3 mm s⁻¹ and were independent of the endothelium. The spread but not the generation of Ca²⁺ spikes was reversibly blocked by the gap junction inhibitor 18β-GA. Thus, smooth muscle gap junctions enable depolarization to spread along resistance arteries, and once regenerative Ca²⁺-based APs occur, spread along the entire length of an artery followed by widespread vasoconstriction.