Noise-Free Visualization of Microscopic Calcium Signaling by Pixel-Wise FittingNovelty and Significance
Rationale: Our insights into physiological and pathophysiological cardiac excitation-contraction coupling has greatly benefited from significant advancement in optical technologies such as high-speed confocal microscopy. This has pushed pixel dwell times into the time domain of nanoseconds, resulting in low signal-to-noise ratios, which have limited data analysis and interpretation.
Objective: Line scan imaging has been and still is dominant in high speed confocal recording. It allows analysis only of a small fraction of a cell's cross section (1.5%), but the appreciation of spatiotemporal fine details of excitation-contraction coupling is instrumental for the further understanding of pathological mechanisms. We aim to provide a novel analysis tool to extract otherwise hidden fine details in cardiac excitation-contraction coupling from high-speed 2-dimensional confocal image series.
Methods and Results: We demonstrate that high-speed 2-dimensional confocal data (150 frames/s) can be analyzed quantitatively by a pixel-wise fitting approach, using a mathematical formalism to phenomenologically describe local calcium transients. Such an approach produces virtually noise-free fluorescence data originating from minute volumes (0.025 femtoliter) and allows extraction of detailed and most importantly quantitative and mechanistically novel information on microscopic calcium signaling and excitation-contraction coupling in a robust manner.
Conclusions: Pixel-wise fitting provides novel insights into cardiac excitation-contraction coupling. Specifically, it revealed microscopic calcium alternans on the level of individual coupling sites. Microscopic calcium alternans is an early precursor of cellular alternans and as such will shed more light onto this mechanism leading to cardiac arrhythmia.
- excitation-contraction coupling
- calcium-induced calcium release
- noise-free image sequences
- arrhythmogenic precursor
- Ca2+ transients
High-resolution live-cell imaging of subcellular signaling events has greatly fostered our understanding of cardiac physiology and pathophysiology.1 In this, a major driving force was the advancement in laser-scanning technology that enabled ultrafast confocal imaging (>100 frames/s) without compromising the optical resolution.2 Although researchers are technically able to follow fast subcellular signaling, such as excitation-contraction coupling (ECC) in cardiac myocytes,3 decreasing the single pixel dwell times and concomitant decreases in the signal-to-noise ratio (SNR) have limited the progress. However, scientific interest is not limited to the occurrence of Ca2+ induced Ca2+ release (CICR) per se but also focuses on where, how much, and how fast the Ca2+ increases occur because the signaling information is encoded in all of these properties.4
Although cardiac Ca2+ transients appear homogeneous under physiological conditions, they arise from the coordinated activity of many microscopic ECC units,3 comprising voltage operated Ca2+ channels in the plasma membrane communicating to CICR channels, the ryanodine receptors (RyR), in the membrane of the Ca2+ storing organelle, the sarcoplasmic reticulum (SR).5 The reliable and repetitive fast communication between these two proteins ensures robust ECC and effective contractility of the myocyte eventually determining the mechanical performance of the entire heart.3 In pathophysiological conditions, such as hypertrophy and heart failure, a partial loss of this coordination has been identified as a major contributor to these diseases.6–8 The tool of choice to investigate local ECC failures occurring in cardiac pathologies is high-speed confocal imaging. Up to now, this field has been dominated by studies using line scan approaches, but line scan images only reflect around 1.5% of the entire accessible confocal cross section albeit imaging can be performed at very high acquisition rates. In addition, line scanning relies on choosing a particular line for sampling and hoping for its representativeness. A better approach to analyze subcellular ECC failures is high-speed imaging of the entire cell cross section. Unfortunately, single pixel dwell times are down to a few 10 nanoseconds when acquiring 2D images at frame rates exceeding 100 Hz, and images are thus concomitant with very low SNR. Although the analysis of the “where” of ECC might still be achieved using data with low SNR, quantitative analysis, such as determining “how much” and “how fast,” requires data with high SNR. Image quality could be improved by generic algorithms that attempt to reduce the noise in the image data by smoothing, wavelet denoising or other filtering.9–11 However, this unavoidably changes the amplitude and/or kinetic properties of recorded signals and often also functionally reduces the spatial resolution.
We introduce an analytic approach to investigate the raw data by phenomenologically describing the time course of local Ca2+ signals with very few and basic assumptions and subsequent pixel-wise fitting of such equations to the recorded data. We will provide virtually pixel-noise free images and movies of cardiac Ca2+ transients and will utilize the information now available to reveal subcellular details of ECC.
An expanded Methods section is available in the Online Data Supplement.
Adult ventricular myocytes from Wistar rats were isolated as previously described.12
High-Resolution, High-Speed Confocal Imaging
Isolated myocytes were loaded with fluo-4 am as previously described13 (1.0 μmol/L for 30 minutes with 15 minutes for deesterification) and imaged with a resonant confocal scanner (Leica TCS SP5, Leica Microsystems GmbH, Wetzlar, Germany) at a frame rate of approximately 150 Hz. For details of solution composition and optical setup as well as for data preparation and calculation of Ca2+ concentrations, please refer to the Online Supplement.
Pixel-wise fitting was performed in MatLab with the “fminsearch” function, which uses the simplex search method of Lagarias et al.14
Pixel-wise Ca2+ transient data were fitted to the following 3 functions, as shown in Online Figure I:(1)
In Equation (1), a gaussian function was used to describe the Ca2+ increase, and an exponential decay function was used for the Ca2+ removal process, where μ denotes the time point of maximal Ca2+ upstroke velocity, σ denotes the half duration at half maximum of the Ca2+ increase, and τ describes the Ca2+ removal process.
In Equation (2), an exponential function was used to describe the Ca2+ increase process, and an exponential decay function was used for the Ca2+ removal process, where μ denotes the onset of the Ca2+ upstroke, τin denotes the time constant of the Ca2+ upstroke, and τout denotes the time constants of Ca2+ decay.
In Equation (3), an exponential function was used to describe the Ca2+ increase process, and a double exponential decay function was used for the Ca2+ removal process, where μ denotes the onset of the Ca2+ upstroke, τin denotes the time constant of the Ca2+ upstroke, and τout1 and τout2 denote the time constants of Ca2+ decay. This is the final equation presented in Figures 1 through 6.(4)
where Ir denotes the recorded Ca2+ concentration (or intensity), and Ic denotes the calculated Ca2+ concentration from the above functions.
2. Preparation of Initial Parameters
The initial values of the temporal parameters (μ, σ, τin, τout or τout1, and τout2) were derived by first fitting the global Ca2+ data to Equations (1),(2), or (3). The initial values of the amplitudes (A) for each pixel were derived from the change in the maximum pixel amplitude. For Equation 3, A1 and A2 were calculated based on the A1/A2 ratio from the global fitting. The initial values of the baseline (B) for each pixel Ca2+ transient were computed either from the average of 10 to 20 time points before the Ca2+ upstroke or from a set concentration (100 nmol/L).
3. Step-by-Step Pixel-Wise Fitting
Pixel-wise fitting was separated into 3 steps (compare Figure 1F). In step 1, the single-pixel Ca2+ transient data were corrected for the baseline, and the amplitude variables [A for Equations (1) and (2), and A1 and A2 for Equation (3)] were fixed, while the temporal parameters were fitted freely using the start values from the global fitting (see step 2, Preparation of Initial Parameters). In step 2, the amplitudes and the basal Ca2+ were fitted from the start values of the Ca2+ transients from the global fitting (see step 2, Preparation of Initial Parameters), while fixing the temporal parameters from step 1. Finally, all temporal parameters were refitted with the obtained pixel amplitudes and baseline values from step 2 as fixed parameters. During step 3, the reconstructed Ca2+ images were calculated with all of the parameters and combined into a single-pixel noise-free image stack.
For correlation analysis and plotting of the data, please refer to the Online Data Supplement.
Simulation of “Ideal” Calcium Transients and Noise Generation
After 3 days in culture, rat ventricular myocytes15 were imaged as described above and additionally stained with 1 μmol/L Di-8-ANEPPS for membrane visualization. Imaging was performed in Tyrode solution containing 5 mmol/L Ca2+ and 5 μmol/L blebbistatin for mechanical uncoupling. The image stacks were denoised based on Poisson unbiased risk estimator (PURE)11 and then subjected to pixel-wise fitting. The reconstructed Ca2+ transient was used as the input signal in the noise evaluation. Noise was added as a convolution of photon statistics Poisson noise and gaussian read-out and amplification noise. We calculated noise levels for different gains as previously described.16
A Novel Approach to Generate Pixel Noise-Free Image Sequences
Figure 1A depicts typical raw images from confocal recordings of a cardiomyocyte during the onset of an electrically evoked Ca2+ transient that were acquired at a frame rate of 146 Hz. It is obvious from the individual images (Figure 1A) and from the plot of the single-pixel fluorescence over time (Figure 1B, blue dots) that single-pixel data are extremely noisy (signal coefficient of variation: 43% at baseline Ca2+). Within such data, several sources of noise limit the detailed analysis of the spatiotemporal aspects of Ca2+ signaling and render interpretation difficult. We designed a fitting algorithm to describe the time course of Ca2+ transients at each pixel purely phenomenologically (without the concept of a mechanistically modeling) (Figure 1E). The optimization procedure of this fitting algorithm is shown in Online Figure I. The descriptive equation (bottom of Figure 1E) comprised a lag phase (Figure 1E, blue), a mono-exponential upstroke (Figure 1E, red), and a biexponential decay phase (Figure 1E, green). When we applied this equation to fit the global Ca2+ transient (Figure 1C), a good representation of the time course became apparent (Figure 1C, minimal residuals in the lower panel). Thus, we used the same approach to extract the underlying Ca2+ transient from single-pixel data (Figure 1B, left). The work flow is detailed in Figure 1F. To determine the appropriate starting values for the pixel-wise fitting (PWF) process, we initially used global fluorescence data (preprocessing in Figure 1F). Thereafter, these parameters were refined in consecutive fitting iterations (each allowing only a subset of variables to run freely) and finalized in a fitting step during which all parameters could run freely (pixel-wise fitting in Figure 1F). At the end of this multistep fitting process, we obtained a single equation describing the time course of the pixel-based Ca2+ transients. This information could subsequently be used to reconstruct the global and local Ca2+ transients (Online Video I). To further illustrate the gain in quality of our approximation, we reconstructed the entire Ca2+ transient image stack and compared it with the raw data (see surface representations of the pseudoline scan images in Figure 1D and Online Video I). The spatial distributions of several quantitative parameters can be extracted exclusively from the PWF data; that is, the detailed spatial distributions of the amplitude (Figure 2A), of the CICR duration (image in Figure 2B), and of the averaged τ for the Ca2+ removal (Figure 2C). A power spectral analysis of the CICR duration showed the characteristic frequency peak at 0.56 μm−1, the sarcomeric spatial frequency (curve in Figure 2B) corresponding to a sarcomere length of 1.78 μm. Considering our spatial resolution in x/y of 0.28 μm, these data strongly supporting the quality of our approach. To quantify the extraction of information in the PWF data relative to the raw data, we analyzed the correlations between consecutive Ca2+ transients in both datasets (Figure 2D). Although there was a rather low correlation and a wide spread in the raw transients (Figure 2D, left panel), we found a high correlation and a tight distribution for the reconstructed transients (Figure 2D; right panel). To quantify the quality of the fit for each individual pixel, we used a PWF-extracted standard deviation of the residual divided by the square root of the intensity (image in Figure 2E) that we determined as “fitting deviation.” Compared with the goodness of fit, R2, this number is rather independent of the intensity and the associated noise, since in our setting noise is a recording parameter rather than a fitting parameter. We regarded single-pixel fittings as “failing” when their residuals exceeded 3 times the standard deviation of the gaussian distribution of the histogram (plot in Figure 2E). In our examples, the number of failing pixels amounted to approximately 0.2% of the total number of pixels (Figure 2F). If these were isolated pixels and their number did not exceed 0.5%, we assumed an appropriate fitting approach for the Ca2+ transient and neglected those failing pixels.
Graded Noise Analysis, Algorithm Robustness, and Reproducibility
To analyze the sensitivity of the algorithm for noise, we performed a graded noise analysis. For this and also to create “worse-case scenarios” with respect to inhomogeneous ECC, we designed an “ideal cell” based on a cardiac myocyte 3 days in culture that was characterized by a reduced T-tubular system.15 We simultaneously visualized the plasma membrane by Di-8-ANEPPS staining (Figure 3A). From Ca2+ transients obtained from this cell, we extracted noise-free input data (Figure 3B). Based on different gain values as described in the Methods section, we added corresponding noise levels to the input transient. This is illustrated as the peak SNR in Figure 3C (left row). With the images containing graded noise we performed pixel wise fitting and the results are presented in Figure 3C, middle and right image columns for amplitude and upstroke, respectively. From these data, we extracted the median of the fitting deviation and plotted it against the peak SNR (Figure 3D). This plot revealed that although the fit was still stable with SNR values below 2, the fitting deviation increased with decreasing SNR. To translate this into the particular fitting parameter, we performed correlation analysis of the fitted amplitude and upstroke 1/τ for the different SNR values relative to the input data (Figure 3E). When considering the amplitude, the correlation was close to 1 for SNR values even lower than 2 (0.954 at a SNR of 1.93). Since the upstroke phase was fitted from only very few points, the correlation levels for the 1/τ distributions did not approach the high values found for the amplitude. However, the images in Figure 3C (right image row) revealed that Ca2+-release sites (sites of fast Ca2+-increase) could still be identified at a gain of 100 (SNR of 2.7). The fitting results for the upstroke parameter would largely benefit from faster image acquisition technologies (see Discussion). For a more theoretical consideration of graded noise evaluation based on amplitude and upstroke gradients, please refer to Online Figure II.
Beside these technical considerations, we addressed the robustness and reproducibility of the algorithm on repetitive transients measured under standard experimental conditions. For healthy cardiac myocytes, it is believed that during steady-state conditions, consecutive electrically evoked Ca2+ transients display similar Ca2+ release patterns.3 We made use of this to verify the reproducibility and robustness of our PWF approach. Figure 4 summarizes the results of such a typical experiment. A rat ventricular myocyte was stimulated at a constant frequency (2 Hz) for approximately 5 minutes to obtain steady-state conditions and thereafter 5 consecutive Ca2+ transients were recorded at 146 Hz. We applied our PWF approach to each of the 5 transients individually and compared the amplitude (Figure 4B and 4D) as well as the Ca2+ upstroke distribution (Figure 4C and 4E) for all transients. Although patterns appear rather similar between transients, we calculated the 2D correlation coefficient for the pairs given in Figure 4F. The result depicted that the correlation between the amplitudes (filled circles) and Ca2+ upstrokes (black triangles) was very high and consistent even between signals that were 5 transients apart.
Microscopic Alternans as a Precursor of Macroscopic Alternans
PWF offers important mechanistic insights that have not been accessible before. To support this, we investigated a phenomenon called Ca2+ transient alternans. This is the cellular equivalent of T-wave alternans in the ECG that is associated with a plethora of disease situations. It reflects alternans of the action potential (AP) repolarization by Ca2+ transients alternating between high and low amplitudes.17,18 Although it is known that a stepwise increase in the stimulation frequency often provokes the occurrence of Ca2+ alternans,19 the process of the onset of alternans remains obscure. In Figures 5 and 6, we compared manifested Ca2+ alternans (right column) with the period preceding these macroscopic alternans (left column) in rat ventricular myocytes. For the latter conditions (Figure 5A through 5C, left column), no obvious changes of the global Ca2+ transients occurred relative to the “healthy” situation apart from a decrease in the global amplitude (Figure 4A, 4B, and 4D). In Figure 5A through 5C (right column), the presence of macroscopic alternans is obvious. In the nonalternans condition, the sites of fast increases were distributed evenly in the cell; however, in the macroscopic alternans condition, a restricted response was evident in only part of the myocytes (Figure 6A and 6B). Despite this difference, we could find coupling sites (Figure 5A, yellow arrows) that surprisingly displayed alternating amplitudes of microscopic Ca2+ transients (Figure 5D). The red traces represent the local Ca2+ transients as a result of PWF. From this, we concluded that despite the different behavior of global or macroscopic Ca2+ transients, the alternans could be identified based on the level of individual coupling sites well before the onset of macroscopic alternans by PWF. Therefore, we refer to the behavior of those Ca2+ transients preceding the macroscopic alternans as microscopic alternans. Quantification of the alternating behavior was achieved by calculating the correlation coefficient between the amplitude distribution images (as depicted in Figure 5A) of the first and all successive transients (Figure 5F). Although the correlation was slightly but significantly changed between the consecutive and alternating transients during microscopic alternans (marked in green; Figure 5G), the macroscopic alternans resulted in an alternation between large positive and negative values for the image correlation coefficient (marked in red, Figure 5F and 5G). This analysis provided further evidence that the behavior of the Ca2+ release machinery was already altered during the period of microscopic alternans. For macroscopic alternans, the negative image correlation coefficient (Figure 5G) pointed to predominantly alternating release sites. In analogy to (Figure 2D), we analyzed the correlations between consecutive Ca2+ transients in raw and PWF datasets and plotted these correlations for 3 transients in (Figure 5H). Although it can be seen that for the raw data almost no differences were apparent between microscopic (left) and macroscopic alternans (right), these scatter plots of the Ca2+ signals depict a rather tight correlation for the microscopic alternans (lower image row in Figure 5H, left column). During macroscopic alternans for consecutive transients such a tight correlation was completely lost (1st and 3rd images in the lower right row, Figure 5H). The two “wings” representing the spatially alternating Ca2+ transients are obvious. The scatterplots between alternating transients (here transient 1 and transient 3; middle image) were indistinguishable from the microscopic situation. These data again exemplify the power of PWF in revealing subcellular details of ECC that were concealed in the noise of the raw images.
Furthermore, we analyzed the Ca2+ upstroke velocity as a measure for the onset of the CICR (compare Figure 6). This parameter that could only be extracted from high speed confocal data in a spatially resolved manner by using PWF. For microscopic alternans, the correlation of the increased velocity (see green circles in Figure 6C) plateaued at values of approximately 0.5 following an initial dip. For the upstroke velocity, the correlation coefficients between the consecutive transients and every second transient (alternating; Figure 6D) were as expected, based on the macroscopic alternans. For the microscopic alternans, the results were slightly more complex. The correlation coefficient was strongly decreased for both the consecutive and alternating transients when investigating the upstroke velocity (green circles in Figure 6D). These data strongly supported our notion that ECC was already flawed in the period preceding global Ca2+ alternans.
In the last decade, optical imaging techniques have pushed the boundary for fast confocal imaging into areas in which photon-statistical limitations have become an increasingly difficult obstacle for the interpretation of the data. Here, often signal averaging does not appear to be appropriate because of the individuality of the signal properties (see, eg, Ca2+ alternans). In cellular Ca2+ signaling, the properties of individual events such as fundamental and elementary Ca2+ signals including Ca2+ blips, quarks, puffs, and sparks play an important role in cellular signaling.20 Traditionally such fast events were recorded with line scans. With a comparable technology as used in this report, this approach allowed much higher acquisition rates in the range of kHz.6,7 This is achieved by sacrificing spatial information since the observation is limited to ≈1.5% of the entire cell cross section. There are applications where the high temporal resolution might be necessary but there are also fundamental questions where a better appreciation of spatial properties (at least from a 2D image) is compulsory. This holds, for example, true for pathological conditions for which characteristic uncoupling between plasma membrane Ca2+ channels and RyR had been suggested (orphaned RyR).8 Very fast confocal 2D (ie, x/y over t) scanning is accompanied by single pixel dwell times of a few 10 nanoseconds leading to single pixel SNR that prohibit researchers from studying Ca2+ signals on this level in greater details. Morad et al performed pioneering work on fast 2D scans of isolated cardiac myocytes,21,22 and, although they already reached a significant gain of knowledge, they experienced the problems described above. They partially circumvented the situation by functionally uncoupling Ca2+ release with 5 mmol/L EGTA intracellularly. For rat atrial myocytes, we have used fast acousto-optic based scanning (120-Hz image rate, 10-kHz line scan rate) without such interventions and revealed distinct but predetermined ECC sites.23
As an alternative to signal averaging, filtering and smoothing of data offers excellent possibilities to increase the SNR,10 but mostly such approaches lead to compromises in the signal amplitude and/or time course. These are both signal properties that are of utmost importance for cardiac Ca2+ signaling.
We introduce a novel and advantageous approach to overcome such technical limitations: PWF of the time course of Ca2+ signals in fast 2D confocal image series. By using an illustrative formalism to phenomenologically describe the time course of cardiac Ca2+ transients, we were able to generate and visualize virtually noise-free single pixel data (compare Figure 1D) of fast cellular and subcellular Ca2+ signals. Ca2+ transients in striated muscle cells and neurons are among the fastest subcellular Ca2+ signals currently known. We therefore validated the PWF approach on Ca2+ transients in electrically paced cardiac myocytes (Figure 2). We extensively tested the algorithm for noise stability (Figure 3 and Online Figure II) and the approach for reproducibility and robustness (Figure 4), and the results are satisfying.
The frame rate of approximately 150 Hz as used in this paper can be regarded as a lower limit to resolve (and fit) details of the ECC such as the Ca2+-upstroke. Recordings at higher acquisition speeds are desirable, especially when acquiring at physiological temperatures of 37°C. There are imaging technologies available, such as acousto-optic scanning21,23 or slit scanning that offer kHz frame rates.24 Furthermore with the accessibility of sCMOS sensors,25 camera based confocal scanners will also be able to approach similar frame rates. Data acquired with these faster recording techniques will benefit from the pixel wise fitting algorithm even more.
However, our approach in its present stage has limitations. It only works for one Ca2+-transient per pixel within a given time interval. For a beating cardiac myocyte such a requirement is sensible and leads to results as presented throughout Figures 1 to 6. If we face repetitive events, such as Ca2+ sparks and Ca2+ waves, if Ca2+ is released in 2 or even more bursts or if Ca2+ release events occur on top of Ca2+ transients, the algorithm in its present form will fail. Therefore, it is necessary to verify these prerequirements before applying PWF to imaging data.
At the moment, the algorithm does not account for movement of the cell. However, the fitting procedure copes well with cellular contraction. Nonetheless, such contractions will unavoidably lead to a minute spatial blurring, but future sophisticated image-by-image registration techniques might overcome such shortcomings. Nevertheless, we minimized contraction artifacts by recording cells under steady-state pacing conditions where contractions are minimal due to postrest behavior. In addition, we chose cells that tightly attached to the inelastic surface of the cover slip. Nonetheless, we do not see an obstacle in applying ECC uncouplers such as blebbistatin.26
Microscopic Versus Macroscopic Alternans
On the level of the ECG, T-wave alternans reflects alternans of action potential repolarization on the single-cell level and alternating Ca2+ transients between high and low amplitudes.17,18 Despite the fact that experimental as well as modeling studies suggest that Ca2+ alternans arises from rhythmic sensitivity changes of the RyR resulting in rhythmic changes of the Ca2+ transient amplitude, experimental evidence usually relies on performing line scans along a single line of the cardiac myocytes17 or cardiac tissue27 to investigate that phenomenon. Such recordings unfortunately only represent a small fraction of the cell's cytosolic volume but analysis of fast 2D confocal recordings has been impossible due to limited SNR.
Although analysis of spatially isolated Ca2+ signals (eg, Ca2+ sparks)15,28 or atrial ECC29 has been possible, the beating ventricular myocyte represents a particular challenge. Because ECC sites are numerous and Ca2+ signals evolving from such sites are distributed in 3 dimensions throughout the entire cell, the resulting Ca2+ signals tend to fuse in time and space. This makes identification of fine details of ECC rather challenging. Various reports have attempted to increase the SNR and to isolate ECC signals by introducing very high concentrations of Ca2+ indicators and/or additional Ca2+ buffers (both in the mM range).22,30 The experimental settings in these reports resulted in high fluorescence signals and in an artificial isolation of Ca2+ release signals. In contrast, PWF can be applied in naive cells that had been loaded with a limited amount of indicator to ensure the lowest amount of Ca2+ buffer possible.
Our algorithm showed its advantages in studying cardiac alternans on 2D image sequences (Figures 5 and 6 and Online Video II). For the first time, we could demonstrate that cardiac ECC was already flawed in the periods preceding macroscopic alternans (a period in which microscopic alternans was revealed). Quantitative information about the subtle details of ECC in cardiac myocytes during the onset of alternans was provided. Using PWF, we were able to extend our investigations not only onto “obvious” parameters of ECC such as the amplitude of local Ca2+ signals but also onto additional important determinants of ECC such as upstroke velocity, an indicator for the goodness and tightness of L-type Ca2+ channel–RyR interactions. In periods preceding global Ca2+ alternans, we could reveal minute alterations of ECC: microscopic alternans. Further investigations into the coordination and accumulation of such local ECC-flaws might enlighten mechanistic details of the procedures contributing to global Ca2+ alternans, electric alternans, and eventually arrhythmogenic events.
Thus, PWF not only allows diminishing the signal noise but also enables researchers to derive quantitative information important to characterize and understand mechanisms underlying subcellular Ca2+ transients.
Extendibility of the Presented Approach
Despite the fact that we tested and validated the PWF approach on cardiac Ca2+ signaling, the applicability of PWF extends well beyond the cardiac research field into all areas of Ca2+ signaling for which the low SNR limits the interpretation of signals, for example, neuronal signal transduction31 and areas of study in which descriptive formalisms for the cellular Ca2+ transients do exist. Technological developments, such as the scientific CMOS cameras, will improve imaging frame rates so that they are on the order of thousands of frames per second.25 In this case, PWF might also be used to reconstruct the time course of individual APs from optical recordings using fluorescent potential sensors.32,33 Therefore, the PWF of APs may not only give results similar to those of classic electrophysiological experiments but also might add the ability to spatially resolve optical AP recordings, possibly yielding new insight into the distribution and propagation of APs.
Sources of Funding
This work was supported by the German Research Foundation (KFO196 and GK1326), the Federal Ministry for Education and Research (CordiLux), and Saarland University Funds (ZFK and HOMFOR).
We thank Prof V. Helms (Bioinformatics, Saarbrücken) and C. Götze (Arivis GmbH, Rostock) for proofreading the manuscript and for helpful discussions.
In April 2012, the average time from submission to first decision for all original research papers submitted to Circulation Research was 12.79 days.
The online-only Data Supplement is available with this article at http://circres.ahajournals.org/lookup/suppl/doi:10.1161/CIRCRESAHA.112.266403/-/DC1.
Non-standard Abbreviations and Acronyms
- action potential
- calcium-induced calcium release
- excitation-contraction coupling
- Poisson unbiased risk estimator
- pixel-wise fitting
- ryanodine receptor
- scientific complementary metal-oxide-semiconductor
- sarcoplasmic reticulum
- signal-to-noise ratio
- Received February 7, 2012.
- Revision received May 10, 2012.
- Accepted May 14, 2012.
- © 2012 American Heart Association, Inc.
- Shorte Spencer L,
- Frischknecht F
- Kaestner L,
- Lipp P
- Bootmann MD,
- Berridge MJ,
- Putney JW,
- Roderick HL
- Fearnley CJ,
- Roderick HL,
- Bootmann MD
- Heinzel FR,
- Bito V,
- Biesmans L,
- Wu M,
- Detre E,
- Wegner von F,
- Claus P,
- Dymarkowski S,
- Maes F,
- Bogaert J,
- Rademakers F,
- D'hooge J,
- Sipido K
- Gomez AM,
- Valdivia HH,
- Cheng H,
- Lederer MR,
- Santana LF,
- Cannell MB,
- McCune SA,
- Altschuld RA,
- Lederer WJ
- Song L,
- Sobie E,
- McCulle S,
- Lederer W,
- Balke C,
- Cheng H
- Kaestner L,
- Scholz A,
- Hammer K,
- Vecerdea A,
- Ruppenthal S,
- Lipp P
- Delpretti S,
- Luisier F,
- Ramani S,
- Blu T
- Weiss JN,
- Nivala M,
- Garfinkel A,
- Qu Z
- Aistrup GL,
- Kelly JE,
- Kapur S,
- Kowalczyk M,
- Sysman-Wolpin I,
- Kadish AH,
- Wasserstrom JA
- Cheng H,
- Lederer WJ
- Cleemann L,
- Wang W,
- Morad M
- Kaestner L
- Kaestner L,
- Lipp P
- Jung G
- Kaestner L,
- Tian Q,
- Lipp P
Novelty and Significance
What Is Known?
Fast confocal imaging of beating cardiac myocytes is traditionally performed only in one spatial dimension thereby restricting representativeness and spatial information.
Ultrafast 2-dimensional confocal data are noisy, and noise-removing procedures inevitably alter the data and its information content.
Calcium alternans in single cardiac myocytes correlate of T-wave alternans in the ECG.
What New Information Does This Article Contribute?
This new method for image processing allows noise removal from microscopic imaging of cardiac myocytes producing noise-free “movies” of beating cells.
Very locally occurring alternans (“microscopic alternans”) are precursors of previously known alternans (“macroscopic alternans”).
Confocal line scans, performed almost exclusively for fast confocal calcium imaging, cover only 1.5% of the entire cross section of cardiac myocytes. We introduce a novel method that allows segregation of “real” calcium transients from image noise in fast confocal 2-dimensional image series, which was could not be analyzed in detail previously due to high noise levels. Using this method to investigate steps preceding calcium alternans, we identified “microscopic alternans” on the level of individual excitation-contraction coupling units. We regard “microscopic alternans” as important mechanistic precursors of cellular or “macroscopic alternans” in the development of arrhythmic behavior in cardiac myocytes. Our novel method is a valuable tool that allows deeper insights and therefore a better understanding of the molecular operation of excitation-coupling units in cardiac excitation-contraction coupling. The application of pixel-wise fitting might also be extended to other important phenomena such as spatially resolved readouts of action potentials by optical imaging.