Nanometry is widely used in biological sciences to analyze the movement of molecules or molecular assemblies in cells and in vivo. In cardiac muscle, a change in sarcomere length (SL) by a mere ∼100 nm causes a substantial change in contractility, indicating the need for the simultaneous measurement of SL and intracellular Ca2+ concentration ([Ca2+]i) in cardiomyocytes at high spatial and temporal resolution. To accurately analyze the motion of individual sarcomeres with nanometer precision during excitation–contraction coupling, we applied nanometry techniques to primary-cultured rat neonatal cardiomyocytes. First, we developed an experimental system for simultaneous nanoscale analysis of single sarcomere dynamics and [Ca2+]i changes via the expression of AcGFP in Z discs. We found that the averaging of the lengths of sarcomeres along the myocyte, a method generally used in today’s myocardial research, caused marked underestimation of sarcomere lengthening speed because of the superpositioning of different timings for lengthening between sequentially connected sarcomeres. Then, we found that after treatment with ionomycin, neonatal myocytes exhibited spontaneous sarcomeric oscillations (cell-SPOCs) at partial activation with blockage of sarcoplasmic reticulum functions, and the waveform properties were indistinguishable from those obtained in electric field stimulation. The myosin activator omecamtiv mecarbil markedly enhanced Z-disc displacement during cell-SPOC. Finally, we interpreted the present experimental findings in the framework of our mathematical model of SPOCs. The present experimental system has a broad range of application possibilities for unveiling single sarcomere dynamics during excitation–contraction coupling in cardiomyocytes under various settings.

The state of the contractile system of cardiac muscle is regulated by a change in the intracellular Ca2+ concentration ([Ca2+]i), regardless of the development stage (Bers, 2002; Colella et al., 2008 and references therein). It is likewise well established that cardiac contractile performance is highly dependent on sarcomere length (SL) in that a change of merely ∼100 nm causes a substantial change in myocardial contractility (i.e., the Frank-Starling relation; Allen and Kentish, 1985; Fukuda et al., 2001; Katz, 2002; Hanft et al., 2008). This intrinsic nature of cardiac myofilaments requires an accurate measurement of SL to enhance our understanding of myocardial dynamic properties, not only in adults but also in neonates whose cardiomyocytes undergo rapid growth for normal heart development via biochemical/paracrine and mechanical signaling (Jacot et al., 2008, 2010; Rodriguez et al., 2011 and references therein). However, because of the lack of technology to measure single SL at high spatial and temporal resolution, sarcomere dynamics and the impact on myocyte motion have not been systematically investigated, especially in neonatal cardiomyocytes.

Although several studies have been conducted by using various types of cardiac specimens, either neonatal or adult, SL is usually averaged along myofibrils. However, given the notion that [Ca2+]i rises locally in cardiomyocytes (in both adults [Bers, 2002] and neonates [Colella et al., 2008]), it is imperative to investigate how the mechanical properties of myocytes are fine-tuned via single sarcomere dynamics. Recently, we applied quantum dots (QDs) to isolated adult cardiomyocytes and successfully measured the length of a single sarcomere at 30-nm spatial precision at 30 fps (Serizawa et al., 2011). However, given the notion that only a small fraction of QDs were incorporated into living myocytes, even with the use of FuGene, the properties of sarcomeric motion in relation to a change in [Ca2+]i could not be analyzed in detail (Serizawa et al., 2011).

α-Actinin is a 100-kD filamentous actin–binding protein that is concentrated in the dense bodies of stress fibers and smooth muscle cells and in the Z discs of striated muscle cells. Dabiri et al. (1997) first demonstrated green fluorescent protein fusion to the amino terminus of α-actinin in living chick embryonic cardiomyocytes and provided evidence that muscle thick and thin filaments were assembled in different cellular locations to form the sarcomeric structure. Recently, by taking advantage of this method, Hersch et al. (2013) estimated the contractile force of rat neonatal cardiomyocytes.

In the present study, we coupled AcGFP to α-actinin and transfected rat neonatal cardiomyocytes to directly visualize single sarcomeric movements at nanometer precision simultaneously with [Ca2+]i. Results clearly show that (a) the coordinated movements of sarcomeres in series produce motion of a given region of a myocyte; (b) single sarcomeric oscillation waveforms during electric field stimulation are indistinguishable from the spontaneous sarcomeric oscillations (SPOCs) occurring at constant, partial [Ca2+] in skinned myocardium (termed cell-SPOC; Sasaki et al., 2006; Ishiwata et al., 2011; Serizawa et al., 2011); and (c) the myosin activator omecamtiv mecarbil (OM) enhances Z-disc movements during cell-SPOC. Moreover, the preset experimental system was useful for the measurement of SL displacement in adult ventricular myocytes. The molecular mechanisms of cardiac sarcomere dynamics under various conditions, as well as the application possibilities of our experimental system, are discussed.

AcGFP expression in the Z discs of neonatal cardiomyocytes

All of the experiments in this study conformed to the Guide for the Care and Use of Laboratory Animals published by the National Institutes of Health (NIH publication No. 85-23, revised in 1996).

Ventricular myocytes were isolated from 1-d-old Wistar rats and were cultured based on a previous publication (Kojima et al., 2003). pAcGFP-actinin plasmid transfection was performed on the third day of incubation by using Lipofectamine LTX (Invitrogen). pAcGFP-actinin plasmids were obtained by inserting the α-actinin gene into pAcGFP-N1 (Takara Bio Inc.).

Microscopic system

Real-time nanoimaging was performed on AcGFP-expressing neonatal myocytes by using an inverted microscope (IX-70; Olympus) equipped with an electron-multiplying charge coupled device camera (iXon +897; Andor Technology) using a 60× oil immersion objective (numerical aperture of 1.45; Olympus; Oyama et al., 2012). Fluorescence images were obtained from myocytes expressed with α-actinin–AcGFP. When changes in [Ca2+]i were measured simultaneously with sarcomeric motion, the myocytes were bathed in 2 mM Ca2+-HEPES-Tyrode’s solution containing 2 µM Fluo-4-AM (Dojindo) for 20 min at 25°C. Thereafter, experiments were performed.

As a general rule of optics, spatial precision decreases with an increase in camera speed (frame rate). Therefore, to obtain high precision in the measurement of SL displacement, we imaged sarcomeric motions at 50 fps in the majority of the experiments in this study (Drago et al., 2012).

Experimental procedure

Nanoscale imaging was performed 1 d after pAcGFP-actinin plasmid transfection. In experiments under electric field stimulation, 2 mM Ca2+-HEPES-Tyrode’s solution was used as the bathing solution, and myocytes were electrically stimulated (Serizawa et al., 2011 and references therein). For experiments on cell-SPOC, myocytes were treated with 4 µM ionomycin (Sigma-Aldrich) to increase the Ca2+ permeability of the membranes. Then, they were bathed in Ca-SPOC solution (pCa 5.75 buffered by Ca2+/10 mM EGTA; see the Methods section of the online supplemental material) containing 4 µM thapsigargin (Sigma-Aldrich) and 200 µM ryanodine (Sigma-Aldrich). OM (Selleck) was initially dissolved in ethanol and diluted with 0.6 µM Ca-SPOC solution. The final concentration of ethanol was 0.1%, having no effect on cell-SPOC.

Except for the experiment on spontaneous beating (performed at 27 ± 0.5°C), all the mechanical experiments were performed at 36.0 ± 0.5°C in this study. When changes in [Ca2+]i were measured simultaneously with sarcomeric motion, the myocytes were bathed in the 2 mM Ca2+-HEPES-Tyrode’s solution containing 2 µM Fluo-4-AM for 20 min at 25 ± 0.5°C. Thereafter, the imaging experiments were performed at 36.0 ± 0.5°C.

Data analysis

The intensity profiles of the fluorescence images were obtained by iQ software (Andor Technology), transferred to ImageJ software (National Institutes of Health), and analyzed by using our Microsoft Excel macro, as described in previous studies (Sasaki et al., 2006; Serizawa et al., 2011). In both the electric field stimulation and cell-SPOC experiments, we obtained the following parameters as indices of sarcomeric contraction: SL oscillation frequency, SL amplitude, sarcomere shortening time (velocity), sarcomere lengthening time (velocity), and Z-disc amplitude. The SL amplitude was defined as the SL at peak lengthening minus that at peak shortening. The shortening (lengthening) time was defined as the length of time from the onset to the end of SL shortening (lengthening). The Z-disc amplitude was defined as the distance of the movement of the Z disc during a cycle of contraction. Data obtained from three consecutive contraction cycles were averaged for the aforementioned parameters.

Statistics

Significant differences were assigned using the paired or unpaired t test as appropriate. Data are expressed as means ± SEM unless otherwise noted, with n representing the number of preparations. Linear regression analyses were performed in accordance with the method used in a previous study (Fukuda et al., 2001). Statistical significance was assumed to be P < 0.05. NS indicates P > 0.05.

Online supplemental material

Supplemental text provides more information about the methods used in this paper. Fig. S1 presents image analysis for a rat neonatal cardiomyocyte transfected with the α-actinin–AcGFP plasmid showing the colocalization of α-actinin–AcGFP with the phase-dense Z discs. Fig. S2 shows the fluctuation analysis of the length of a single sarcomere in an α-actinin–AcGFP-expressing neonatal cardiomyocyte in the presence of Fluo-4. Fig. S3 shows changes in SL and fluorescence intensity (FI) in a cardiomyocyte upon electric stimulation in the absence of Fluo-4. Fig. S4 shows α-actinin–AcGFP expression and SL nanometry in an adult cardiomyocyte of the rat. Fig. S5 shows the time course of changes in the lengths of sarcomeres during spontaneous beating in an AcGFP-expressing neonatal cardiomyocyte. Fig. S6 presents histograms showing averaged data versus individual data on amplitude, shortening, and lengthening properties during spontaneous beating. Fig. S7 shows sarcomeric motions during spontaneous beating in a neonatal cardiomyocyte without transfection of the α-actinin–AcGFP plasmid. Fig. S8 shows the relationship of electric stimulation frequency versus time to peak FI in intact neonatal cardiomyocytes. Fig. S9 shows cell-SPOC properties in AcGFP-expressing neonatal cardiomyocytes. Fig. S10 shows waveform analyses of cell-SPOC in the absence and presence of OM. Videos 1 and 2 show neonatal cardiomyocytes expressed with α-actinin–AcGFP in Z discs at rest. Video 3 shows reconstruction of a three-dimensional image of a neonatal cardiomyocyte. Video 4 shows that the FI from a near full range of a cardiomyocyte does not change during the course of sarcomeric contraction. Video 5 shows an adult cardiomyocyte expressed with α-actinin–AcGFP in Z discs at rest. Videos 6 and 7 present neonatal cardiomyocytes showing spontaneous beating. Videos 8 and 9 present neonatal cardiomyocytes showing contractions in response to electric field stimulation. Videos 10–12 present neonatal cardiomyocytes showing cell-SPOC.

First, the α-actinin–AcGFP construct produced fluorescent-labeled protein in Z discs within 24 h after transfection when introduced into neonatal myocytes (as demonstrated in neonatal chick [Dabiri et al., 1997] and rat myocytes [Hersch et al., 2013]; Fig. 1 A). This fusion protein was incorporated into the Z discs of normal cellular structures rich in endogenous α-actinin, and α-actinin–AcGFP colocalized in a banded pattern with the phase-dense Z discs (Fig. S1).

As shown in the plot profile in Fig. 1 A, SL averaged ∼1.95 µm under the resting condition (1.965 ± 0.021 µm; from 10 myocytes × 3 sarcomeres). This value is similar to that reported previously in rat neonatal myocytes (Hersch et al., 2013), and it is slightly longer than that reported by others in isolated adult rat and mouse myocytes (i.e., ∼1.70–1.85 µm; Bub et al., 2010; King et al., 2011 and references therein; compare with Fig. S4). The differential value of SL is likely to primarily underlie the different experimental settings. Namely, in the present study, neonatal cardiomyocytes were cultured, as opposed to adult cardiomyocytes isolated from the heart in these previous studies. Also, it is likely that the differential titin (connectin) expression profile in neonatal versus adult cardiomyocytes may in part underlie the differential SL values. Indeed, Lahmers et al. (2004) reported that the titin expression profile in rats differs in adults versus neonates, with stiff N2B titin predominantly expressed in adults in contrast to more compliant neonatal isoforms (e.g., N2BA isoforms) expressed in neonates, with titin-based passive force ∼0 mN/mm2 at SLs up to ∼2.0 µm. Given the widely accepted role of titin as a molecular ruler to determine the SL during myofibrillogenesis in striated muscles (Gregorio et al., 1999; Udaka et al., 2008 and references therein), the slightly longer SL observed in the present study under the isometric condition may be partly caused by N2BA titin isoforms predominantly expressed in neonatal myocytes.

The transfection efficacy of AcGFP expression in Z discs (i.e., the ratio of AcGFP-expressing myocytes [determined by fluorescence observation] to all myocytes attached to the glass plate [determined by phase-contrast observation]) was 38.14 ± 5.47% (number of observations, 17, at the image size of 138 × 138 µm with a 60× oil immersion objective).

Nanometry is a powerful technique that is widely used in single-molecule biophysics by many researchers for things like actomyosin interaction (Uemura et al., 2004; Komori et al., 2009; Kaya and Higuchi, 2010) and was recently used in molecular imaging in cells (Yoo et al., 2008; Watanabe et al., 2013) and even in vivo (Tada et al., 2007; Gonda et al., 2010) in various bioscience fields. In the present study, AcGFP was clearly expressed in Z discs in neonatal cardiomyocytes with relative ease via lipofection, providing a suitable experimental system for nanometry on a single SL displacement measurement (Fig. 1 B). Under the present experimental condition, the SL displacement precision was 3 nm at 50 fps (Fig. 1 B and Video 1). Likewise, it was 8 nm at 50 fps in the presence of the Ca2+ indicator Fluo-4 (Fig. S2 and Video 2). The present method (i.e., SL nanometry) provides high accuracy of the positioning of the AcGFP-expressed Z disc by fitting the 3 pixels, including the peak point (with the width of 0.81 µm), thereby minimizing the influence of the fluorescence of the neighboring Z discs (Fig. 1 B). Given a relatively short range of SL in cardiac preparations compared with that in skeletal muscle preparations, it can be considered that the present analysis provides a powerful tool in quantitating the SL displacement in cardiac muscle.

By moving the stage of the microscope (i.e., at 0.496 µm/s) in the Z direction (Video 3), the thickness of the myocyte in Fig. 1 A was estimated as ∼6 µm. However, considering the bleed over of the AcGFP fluorescence signal, this will be an overestimated value. To accurately measure the thickness of the myocyte cultured under the present setting, it will be necessary to measure the distance between the fluorescence beads (or nanoparticles) with a known diameter adhered to the top and bottom of the myocyte.

It is worth noting that because of the similar wavelengths of AcGFP (peak: 505 nm) and Fluo-4 (peak: 518 nm), the present experimental system enables simultaneous imaging of the SL displacement and [Ca2+]i changes without the use of double light-path microscopy. As shown in Fig. S3 (also see Video 4), the FI obtained from the near full range of a cardiomyocyte did not change during the course of sarcomeric contraction (from ∼2.0 to ∼1.8 µm) in the absence of Fluo-4 during spontaneous beating, indicating that there was no interference of AcGFP and Fluo-4 fluorescence.

Clear striation patterns were likewise observed under confocal microscopy in isolated adult cardiomyocytes of the rat by using the adenovirus vector system, enabling us to measure single SL displacement with 15-nm precision at 30 fps (Fig. S4 and Video 5). It should be noted that this value is approximately twofold greater than that obtained in our previous study using QDs conjugated with anti–α-actinin antibody (i.e., 30-nm precision at 30 fps; Serizawa et al., 2011), indicating the usefulness of our experimental system in the analysis of sarcomere dynamics in isolated adult cardiomyocytes. In the present study, however, we systematically investigated the single sarcomere dynamics in cultured neonatal cardiomyocytes, taking into consideration the following three merits. First, experiments can be performed under the isometric condition (hence the results can be simply extrapolated to the previous findings obtained from skinned fibers/myocytes of various cardiac specimens). Second, a very high precision (∼3 nm) in the measurement of SL displacement can be obtained with no need for a confocal system because of a thin layer of myofibrils (compare with Fig. 1 A). Third, AcGFP can be routinely expressed with a success rate of ∼40% in Z discs via lipofection (i.e., no need for advanced gene expression techniques). Future research will need to thoroughly investigate the sarcomere dynamics in both healthy and diseased adult cardiomyocytes at nanometer precision.

It is well known that spontaneous, periodic releases of Ca2+ occur from the SR in cultured neonatal myocytes under physiological conditions, subsequently inducing contractions with large sarcomeric amplitudes (Jacot et al., 2008, 2010; Rodriguez et al., 2011; Hersch et al., 2013). Because the advantage of the present method is to track the individual SL instead of measuring the mean value, we conducted simultaneous imaging of single sarcomeres and [Ca2+]i in α-actinin–AcGFP-expressing myocytes during spontaneous beating (at ∼1 Hz; see Fig. 2 A for six consecutive beats; also see Video 6). The striations of AcGFP fluorescence were clearly seen along the myocyte even after the loading of the Ca2+ indicator Fluo-4 during the course of beating, with Fluo-4 fluorescence increasing only in the SL shortening phase. We then conducted in-depth SL analyses on the second and third beats (Fig. 2 B). It was found that sarcomeres first shortened by ∼0.3 µm in response to a rise in [Ca2+]i (phase 1), followed by quick lengthening (phase 2) and subsequent shortening to the original SL (phase 3). Although phase 1 was observed in all sarcomeres within the myocyte, phase 2 (and the subsequent phase 3) was observed in a time-dependent and SL-dependent manner (Fig. S5), independent of [Ca2+]i (see Fig. 2 B, top, a kymograph, for similar [Ca2+]i changes in various sarcomeres). Namely, sarcomeres with lengths longer than the mean value (∼1.95 µm in this myocyte) were lengthened in a consecutive manner, and shorter sarcomeres did not exhibit the lengthening phase, coupled presumably with a decrease in load on myofibrils. Our analysis revealed that the averaging of SL measurement caused marked underestimation of the lengthening speed as a result of delineation of the timing for lengthening between sequentially connecting sarcomeres (which occurred in a time-dependent and length-dependent manner). The histograms in Fig. S6 show that mean values differ from the individual values on shortening and lengthening properties and amplitude.

Here, it may be suggested that some sarcomeres are overstretched under the present experimental setting (i.e., up to ∼2.8 µm; Fig. 2, A and B). However, Opitz et al. (2004) demonstrated in rats that N2BA titin-expressing fetal/neonatal myofibrils produce passive force ∼15 times lower than that of adult myofibrils, and fetal/neonatal myofibrils can be reversibly elongated up to 3.5–4.0 µm (compare with 2.5–3.0 µm in N2B titin-expressing adult myofibrils). It is therefore likely that depending on the condition, sarcomeres can be elongated up to ∼2.8 µm in cultured neonatal cardiomyocytes of the rat, as observed in the present study. Future research will clarify the molecular mechanism of the propagation of lengthening/shortening of sarcomeres along myofibrils during spontaneous beating.

As evident by the finding of nonuniform sarcomere motion even within the same cell, it is reasonable to consider that our SL nanometry has broad potential to reveal slight but fundamental changes in sarcomere dynamics in neonatal as well as adult (Fig. S4) cardiomyocytes under various conditions. We consider that SL nanometry will enable nanoscale measurements of the displacement of the distance between dye-labeled T tubules in the isolated heart, not only when the heart is at rest (Bub et al., 2010; Serizawa et al., 2011; Botcherby et al., 2013; Inoue et al., 2013) but also in motion (e.g., in the beating heart in vivo).

The magnitude of SL shortening during spontaneous beating in myocytes with no AcGFP expression was similar to that observed in AcGFP-expressing myocytes (i.e., ∼0.3 µm; Fig. S7 and Video 7). This demonstrates that AcGFP expression in Z discs does not hinder sarcomeric contractile properties.

Recently, Llewellyn et al. (2008) successfully imaged sarcomere motion with the second harmonic generation (SHG) method. SHG allows for label-free imaging of the A band of the sarcomere at a subsarcomere precision and thus appears to be versatile for the investigation of sarcomere dynamics in various muscle types. However, the present method using the α-actinin–AcGFP expression provides (a) higher precision in the analysis of SL displacement (3 and 15 nm in neonatal and adult cardiomyocytes, respectively, as compared with 20–50 nm by using SHG; Boulesteix et al., 2004; Llewellyn et al., 2008); (b) nanoscale analysis of individual sarcomere dynamics without the use of smoothing techniques (which are frequently used in SHG; Llewellyn et al., 2008); (c) sharper fluorescence images based on the smaller width of Z discs than that of the A band; and (d) steady fluorescence images independent of the viewing angle (as compared with SHG; Plotnikov et al., 2006). In addition, our system allows for the simultaneous visualization of individual sarcomere dynamics and [Ca2+]i changes at one wavelength excitation (with the precision of SL displacement of 8 nm in the presence of the Ca2+-sensitive dye Fluo-4). Therefore, we consider that the present method using α-actinin–AcGFP expression is a more suitable choice than SHG for the analysis of excitation–contraction coupling in living cardiomyocytes at the individual sarcomere level.

Next, to investigate whether SL nanometry is indeed useful in analyzing sarcomere dynamics under conditions with electric field stimulation, i.e., an experimentation method frequently conducted in cardiac physiology, we simultaneously recorded the motion of single sarcomeres and changes in [Ca2+]i at various frequencies. Fig. 3 A shows a typical Fluo-4–loaded AcGFP-expressed myocyte at rest and at peak shortening and peak lengthening. Clear striation patterns at both phases indicate that our SL measurement system is of high precision regardless of the level of [Ca2+]i. Fig. 3 B represents typical recordings showing changes in SL and Fluo-4–dependent FI at 4 Hz (Video 8; also see Video 9 for stimulation at 2 Hz in the same myocyte). All the sarcomeres started to shorten (lengthen) almost simultaneously in response to a rise (decrease) in FI. By analyzing the waveform properties of single sarcomeres as a function of stimulation frequency, we found that the lengthening time decreased with little change in shortening time, resulting in a significant difference between shortening time and lengthening time at 3, 4, and 5 Hz (shortening time > lengthening time; P < 0.05; Fig. 3 C). However, the amplitude was unaffected by a change in stimulation frequency, which is consistent with the finding in a previous study with rat adult cardiomyocytes placed under the isometric condition by using carbon fibers (Peterson et al., 2013).

Time to peak FI significantly decreased as frequency increased (Fig. S8), as previously demonstrated in adult mouse cardiomyocytes (Werdich et al., 2008). It is likewise important to note that diastolic FI significantly increased in a frequency-dependent manner (P < 0.05 at 3 Hz and higher; Fig. 3 D; Serizawa et al., 2011), despite the fact that these experiments were conducted at 36 ± 0.5°C. Thus, given high heartbeat frequencies in neonatal rats (∼240 bpm on day 0 and ∼300 bpm on day 1; Smotherman et al., 1991), it is considered that diastolic [Ca2+]i is elevated in physiological settings, presumably as a result of insufficient SR Ca2+ pump activities (Serizawa et al., 2011 and references therein).

We previously reported that by using QD-treated isolated adult cardiomyocytes, a sarcomeric contraction cycle consisted of shortening followed by slow lengthening at a low frequency of 1 Hz, and a phase shift occurred upon an increase in stimulation frequency, in that the cycle consisted of shortening followed by quick lengthening at physiologically high frequencies (Serizawa et al., 2011). However, in the present study, shortening time was longer over the range of frequencies tested (Fig. 3 C). We consider that the difference in the waveform at a low frequency in the previous study (Serizawa et al., 2011) and the present study results primarily from the difference in the experimental condition (isotonic vs. isometric), rather than from the difference in the preparation (adult vs. neonatal).

Then, to further investigate the capability of SL nanometry under conditions in which the highest spatial resolution is obtained (i.e., with no interference of [Ca2+]i changes), we analyzed sarcomeric waveform properties during SPOCs in ionomycin (Kauffman et al., 1980; Morgan and Jacob, 1994)-treated neonatal cardiomyocytes with SR functions blocked by ryanodine (Sutko et al., 1997) and thapsigargin (Davidson and Varhol, 1995; Treiman et al., 1998). Fig. 4 A shows an epi-illumination image of an α-actinin–AcGFP-expressing cardiomyocyte (left) and the waveforms of SPOC at two different time phases (right; Video 10; also see Video 11 for cell-SPOC with a low frequency of ∼1 Hz). Despite different magnitudes of Z-disc movement (at time phases 1 and 2), the amplitude of sarcomeric oscillation was similar (Fig. 4 A). As previously reported in studies on adult cardiac muscles of various animals (Sasaki et al., 2006; Ishiwata et al., 2011; Serizawa et al., 2011), cell-SPOC waves consisted of slow shortening followed by quick lengthening. In contrast to the situation in skinned adult cardiomyocytes in which the SPOC frequencies are relatively constant (and varies according to animal species; Sasaki et al., 2006; Ishiwata et al., 2011; Serizawa et al., 2011), relatively large cell-to-cell variations were observed (i.e., from ∼1 to ∼4.5 Hz; 2.08 ± 0.11 Hz on average from 31 myocytes) in ionomycin-treated neonatal myocytes, with the amplitude inversely correlated with frequency (Fig. 4 B). Likewise, the shortening time/lengthening time ratio decreased with an increase in frequency (Fig. 4 C) as a result of a frequency-dependent decrease in shortening time with no significant change in lengthening time (Fig. S9). It should be stressed that at a physiologically relevant frequency of 4 Hz, the values in the relationships in Fig. 4 (B and C) became statistically indistinguishable between electric field stimulation and cell-SPOC, suggesting that the sawtooth waveform is intrinsic to sarcomeres themselves, with and without changes in [Ca2+] (Sasaki et al., 2006; Ishiwata et al., 2011). When converted to velocity, it was found that shortening velocity increased and lengthening velocity decreased, both in a linear fashion, with an increase in frequency (Fig. 4 D). The present findings suggest that sarcomeric auto-oscillatory properties may operate physiologically upon relaxation by facilitating relaxation (lengthening) to adjacent sarcomeres (Stehle et al., 2006). However, care should be taken in the interpretation of the data; because the stiffness of the glass is higher than that of cardiac tissue, the sarcomere dynamics observed in the present study may reflect the properties of a cardiomyocyte attached to the glass surface rather than that under physiological situations (i.e., in the heart).

Here, it should be noted that although cell-SPOC occurs at partial activation (pCa 5.75 in the present study) and its waveform properties are similar to those obtained under electric field stimulation (compare with Fig. 4), cell-SPOC should be regarded as a useful experimental system that reveals important characteristics of sarcomere dynamics (e.g., shortening and lengthening properties, oscillation frequency, and the magnitude of contraction [Z-disc amplitude]) in living cardiomyocytes. Future studies are needed to conduct SL nanometry in the beating heart in vivo to investigate whether or not sarcomeric auto-oscillations indeed occur under physiological conditions at various developmental stages.

Recently, we constructed a mathematical model of SPOC (Ishiwata et al., 2011; Sato et al., 2011, 2013). In the model analysis, we demonstrated that an increase in the myosin attachment rate (α) results in an increase in amplitude and a decrease in frequency via an increase in shortening time. We therefore consider that a variance of regional [Ca2+]i around myofilaments in neonatal myocytes (and thus a variance in α; see Ishiwata et al., 2011; Sato et al., 2011, 2013 for effects of Ca2+ on the myosin attachment rate), albeit by a small magnitude, modulates the properties of cell-SPOC. Consistent with this interpretation, at frequencies <1.5 Hz, the values of SL at peak shortening became <1.6 µm (i.e., shorter than the A-band length; Fig. 4 E), indicating that substantial actomyosin-based force was generated upon shortening. With an increase in frequency, maximal SL decreased and minimal SL increased, both in a linear fashion, which is consistent with the notion that internal load along myofibrils is reduced as frequency increases (coupled with a decrease in α).

Finally, to explore whether SL nanometry can be used for the analysis of pharmacological actions of cardiac drugs, we investigated the effects of the actomyosin activator OM (Malik et al., 2011) on cell-SPOC. It has been reported that OM specifically targets and activates myocardial ATPase and improves energy utilization (Malik et al., 2011). In the present study, 0.6 µM OM did not affect amplitude of sarcomeric oscillations; rather, it significantly enhanced the Z-disc motion along myofibrils and shifted the relationship of SL amplitude versus Z-disc motion upward (Fig. 5 A and Video 12). The enhanced Z-disc amplitude likely results from an increase in actomyosin-based force along myofibrils and can be used as an index of contractility (as predicted by our SPOC model; Sato et al., 2011). OM increased (decreased) shortening time (shortening velocity; Fig. 5 B), but it did not significantly affect the lengthening properties (Fig. 5 C), thereby decreasing the cell-SPOC frequency (i.e., 1.63 ± 0.16 Hz on average from eight myocytes; P < 0.05 compared with the value obtained in the absence of OM, 2.08 ± 0.11 Hz). These experimental findings can be qualitatively explained by our SPOC model, which takes into account the “molecular friction” (i.e., the friction between actin and myosin; η [∼1/β, where β is the crossbridge detachment rate] in our SPOC model; Sato et al., 2011) in addition to α (myosin attachment rate). Table 1 summarizes the possible effects of OM on α and η and ensuing changes in cell-SPOC properties. If OM, as an actomyosin activator, increases α and decreases β, then η will be increased, resulting in a decrease in shortening velocity (and hence frequency), as found in the present study (Fig. 5 and Fig. S10). Our model likewise predicts that an increase in η decreases amplitude (Sato et al., 2011), and, accordingly, this effect will offset the increase in amplitude coupled with an increase in α (see Fig. 5 A for no significant effect of OM on cell-SPOC amplitude).

In conclusion, we developed a novel experimental system with cardiomyocytes allowing for nanoimaging of sarcomeric dynamic properties. This system is useful for analyzing cardiac excitation–contraction coupling at nanometer precision and for determining the changes in health and disease during development in future studies.

We thank Dr. Yuta Shimamoto for critical reading of the manuscript.

This work was supported in part by Grants-in-Aid for Scientific Research (B) (to N. Fukuda), Challenging Exploratory Research (to N. Fukuda and F. Kobirumaki-Shimozawa), Scientific Research on Innovative Areas (to N. Fukuda), and Specially Promoted Research and Scientific Research (S) (to S. Ishiwata) from the Ministry of Education, Culture, Sports, Science and Technology of Japan. This study was also supported in part by the Sasagawa Scientific Research Grant from the Japan Science Society (to S.A. Shintani).

The authors declare no competing financial interests.

Richard L. Moss served as editor.

Allen
D.G.
,
Kentish
J.C.
.
1985
.
The cellular basis of the length-tension relation in cardiac muscle
.
J. Mol. Cell. Cardiol.
17
:
821
840
.
Bers
D.M.
2002
.
Cardiac excitation-contraction coupling
.
Nature.
415
:
198
205
.
Botcherby
E.J.
,
Corbett
A.
,
Burton
R.A.
,
Smith
C.W.
,
Bollensdorff
C.
,
Booth
M.J.
,
Kohl
P.
,
Wilson
T.
,
Bub
G.
.
2013
.
Fast measurement of sarcomere length and cell orientation in Langendorff-perfused hearts using remote focusing microscopy
.
Circ. Res.
113
:
863
870
.
Boulesteix
T.
,
Beaurepaire
E.
,
Sauviat
M.P.
,
Schanne-Klein
M.C.
.
2004
.
Second-harmonic microscopy of unstained living cardiac myocytes: measurements of sarcomere length with 20-nm accuracy
.
Opt. Lett.
29
:
2031
2033
.
Bub
G.
,
Camelliti
P.
,
Bollensdorff
C.
,
Stuckey
D.J.
,
Picton
G.
,
Burton
R.A.
,
Clarke
K.
,
Kohl
P.
.
2010
.
Measurement and analysis of sarcomere length in rat cardiomyocytes in situ and in vitro
.
Am. J. Physiol. Heart Circ. Physiol.
298
:
H1616
H1625
.
Colella
M.
,
Grisan
F.
,
Robert
V.
,
Turner
J.D.
,
Thomas
A.P.
,
Pozzan
T.
.
2008
.
Ca2+ oscillation frequency decoding in cardiac cell hypertrophy: Role of calcineurin/NFAT as Ca2+ signal integrators
.
Proc. Natl. Acad. Sci. USA.
105
:
2859
2864
.
Dabiri
G.A.
,
Turnacioglu
K.K.
,
Sanger
J.M.
,
Sanger
J.W.
.
1997
.
Myofibrillogenesis visualized in living embryonic cardiomyocytes
.
Proc. Natl. Acad. Sci. USA.
94
:
9493
9498
.
Davidson
G.A.
,
Varhol
R.J.
.
1995
.
Kinetics of thapsigargin-Ca2+-ATPase (sarcoplasmic reticulum) interaction reveals a two-step binding mechanism and picomolar inhibition
.
J. Biol. Chem.
270
:
11731
11734
.
Drago
I.
,
De Stefani
D.
,
Rizzuto
R.
,
Pozzan
T.
.
2012
.
Mitochondrial Ca2+ uptake contributes to buffering cytoplasmic Ca2+ peaks in cardiomyocytes
.
Proc. Natl. Acad. Sci. USA.
109
:
12986
12991
.
Fukuda
N.
,
Sasaki
D.
,
Ishiwata
S.
,
Kurihara
S.
.
2001
.
Length Dependence of Tension Generation in Rat Skinned Cardiac Muscle: Role of Titin in the Frank-Starling Mechanism of the Heart
.
Circulation.
104
:
1639
1645
.
Gonda
K.
,
Watanabe
T.M.
,
Ohuchi
N.
,
Higuchi
H.
.
2010
.
In vivo nano-imaging of membrane dynamics in metastatic tumor cells using quantum dots
.
J. Biol. Chem.
285
:
2750
2757
.
Gregorio
C.C.
,
Granzier
H.
,
Sorimachi
H.
,
Labeit
S.
.
1999
.
Muscle assembly: a titanic achievement?
Curr. Opin. Cell Biol.
11
:
18
25
.
Hanft
L.M.
,
Korte
F.S.
,
McDonald
K.S.
.
2008
.
Cardiac function and modulation of sarcomeric function by length
.
Cardiovasc. Res.
77
:
627
636
.
Hersch
N.
,
Wolters
B.
,
Dreissen
G.
,
Springer
R.
,
Kirchgeßner
N.
,
Merkel
R.
,
Hoffmann
B.
.
2013
.
The constant beat: cardiomyocytes adapt their forces by equal contraction upon environmental stiffening
.
Biol. Open.
2
:
351
361
.
Inoue
T.
,
Kobirumaki-Shimozawa
F.
,
Kagemoto
T.
,
Fujii
T.
,
Terui
T.
,
Kusakari
Y.
,
Hongo
K.
,
Morimoto
S.
,
Ohtsuki
I.
,
Hashimoto
K.
,
Fukuda
N.
.
2013
.
Depressed Frank-Starling mechanism in the left ventricular muscle of the knock-in mouse model of dilated cardiomyopathy with troponin T deletion mutation ΔK210
.
J. Mol. Cell. Cardiol.
63
:
69
78
.
Ishiwata
S.
,
Shimamoto
Y.
,
Fukuda
N.
.
2011
.
Contractile system of muscle as an auto-oscillator
.
Prog. Biophys. Mol. Biol.
105
:
187
198
.
Jacot
J.G.
,
McCulloch
A.D.
,
Omens
J.H.
.
2008
.
Substrate stiffness affects the functional maturation of neonatal rat ventricular myocytes
.
Biophys. J.
95
:
3479
3487
.
Jacot
J.G.
,
Kita-Matsuo
H.
,
Wei
K.A.
,
Chen
H.S.
,
Omens
J.H.
,
Mercola
M.
,
McCulloch
A.D.
.
2010
.
Cardiac myocyte force development during differentiation and maturation
.
Ann. N.Y. Acad. Sci.
1188
:
121
127
.
Katz
A.M.
2002
.
Ernest Henry Starling, his predecessors, and the “Law of the Heart”
.
Circulation.
106
:
2986
2992
.
Kauffman
R.F.
,
Taylor
R.W.
,
Pfeiffer
D.R.
.
1980
.
Cation transport and specificity of ionomycin. Comparison with ionophore A23187 in rat liver mitochondria
.
J. Biol. Chem.
255
:
2735
2739
.
Kaya
M.
,
Higuchi
H.
.
2010
.
Nonlinear elasticity and an 8-nm working stroke of single myosin molecules in myofilaments
.
Science.
329
:
686
689
.
King
N.M.
,
Methawasin
M.
,
Nedrud
J.
,
Harrell
N.
,
Chung
C.S.
,
Helmes
M.
,
Granzier
H.
.
2011
.
Mouse intact cardiac myocyte mechanics: cross-bridge and titin-based stress in unactivated cells
.
J. Gen. Physiol.
137
:
81
91
.
Kojima
K.
,
Moriguchi
H.
,
Hattori
A.
,
Kaneko
T.
,
Yasuda
K.
.
2003
.
Two-dimensional network formation of cardiac myocytes in agar microculture chip with 1480 nm infrared laser photo-thermal etching
.
Lab Chip.
3
:
292
296
.
Komori
T.
,
Nishikawa
S.
,
Ariga
T.
,
Iwane
A.H.
,
Yanagida
T.
.
2009
.
Simultaneous measurement of nucleotide occupancy and mechanical displacement in myosin-V, a processive molecular motor
.
Biophys. J.
96
:
L4
L6
.
Lahmers
S.
,
Wu
Y.
,
Call
D.R.
,
Labeit
S.
,
Granzier
H.
.
2004
.
Developmental control of titin isoform expression and passive stiffness in fetal and neonatal myocardium
.
Circ. Res.
94
:
505
513
.
Llewellyn
M.E.
,
Barretto
R.P.
,
Delp
S.L.
,
Schnitzer
M.J.
.
2008
.
Minimally invasive high-speed imaging of sarcomere contractile dynamics in mice and humans
.
Nature.
454
:
784
788
.
Luther
P.K.
,
Bennett
P.M.
,
Knupp
C.
,
Craig
R.
,
Padrón
R.
,
Harris
S.P.
,
Patel
J.
,
Moss
R.L.
.
2008
.
Understanding the organisation and role of myosin binding protein C in normal striated muscle by comparison with MyBP-C knockout cardiac muscle
.
J. Mol. Biol.
384
:
60
72
.
Malik
F.I.
,
Hartman
J.J.
,
Elias
K.A.
,
Morgan
B.P.
,
Rodriguez
H.
,
Brejc
K.
,
Anderson
R.L.
,
Sueoka
S.H.
,
Lee
K.H.
,
Finer
J.T.
et al
.
2011
.
Cardiac myosin activation: a potential therapeutic approach for systolic heart failure
.
Science.
331
:
1439
1443
.
Morgan
A.J.
,
Jacob
R.
.
1994
.
Ionomycin enhances Ca2+ influx by stimulating store-regulated cation entry and not by a direct action at the plasma membrane
.
Biochem. J.
300
:
665
672
.
Opitz
C.A.
,
Leake
M.C.
,
Makarenko
I.
,
Benes
V.
,
Linke
W.A.
.
2004
.
Developmentally regulated switching of titin size alters myofibrillar stiffness in the perinatal heart
.
Circ. Res.
94
:
967
975
.
Oyama
K.
,
Mizuno
A.
,
Shintani
S.A.
,
Itoh
H.
,
Serizawa
T.
,
Fukuda
N.
,
Suzuki
M.
,
Ishiwata
S.
.
2012
.
Microscopic heat pulses induce contraction of cardiomyocytes without calcium transients
.
Biochem. Biophys. Res. Commun.
417
:
607
612
.
Peterson
P.
,
Kalda
M.
,
Vendelin
M.
.
2013
.
Real-time determination of sarcomere length of a single cardiomyocyte during contraction
.
Am. J. Physiol. Cell Physiol.
304
:
C519
C531
.
Plotnikov
S.V.
,
Millard
A.C.
,
Campagnola
P.J.
,
Mohler
W.A.
.
2006
.
Characterization of the myosin-based source for second-harmonic generation from muscle sarcomeres
.
Biophys. J.
90
:
693
703
.
Rodriguez
A.G.
,
Han
S.J.
,
Regnier
M.
,
Sniadecki
N.J.
.
2011
.
Substrate stiffness increases twitch power of neonatal cardiomyocytes in correlation with changes in myofibril structure and intracellular calcium
.
Biophys. J.
101
:
2455
2464
.
Sasaki
D.
,
Fukuda
N.
,
Ishiwata
S.
.
2006
.
Myocardial sarcomeres spontaneously oscillate with the period of heartbeat under physiological conditions
.
Biochem. Biophys. Res. Commun.
343
:
1146
1152
.
Sato
K.
,
Ohtaki
M.
,
Shimamoto
Y.
,
Ishiwata
S.
.
2011
.
A theory on auto-oscillation and contraction in striated muscle
.
Prog. Biophys. Mol. Biol.
105
:
199
207
.
Sato
K.
,
Kuramoto
Y.
,
Ohtaki
M.
,
Shimamoto
Y.
,
Ishiwata
S.
.
2013
.
Locally and Globally Coupled Oscillators in Muscle
.
Phys. Rev. Lett.
111
:
108104
108109
.
Serizawa
T.
,
Terui
T.
,
Kagemoto
T.
,
Mizuno
A.
,
Shimozawa
T.
,
Kobirumaki
F.
,
Ishiwata
S.
,
Kurihara
S.
,
Fukuda
N.
.
2011
.
Real-time measurement of the length of a single sarcomere in rat ventricular myocytes: a novel analysis with quantum dots
.
Am. J. Physiol. Cell Physiol.
301
:
C1116
C1127
.
Smotherman
W.P.
,
Robinson
S.R.
,
Ronca
A.E.
,
Alberts
J.R.
,
Hepper
P.G.
.
1991
.
Heart rate response of the rat fetus and neonate to a chemosensory stimulus
.
Physiol. Behav.
50
:
47
52
.
Stehle
R.
,
Solzin
J.
,
Iorga
B.
,
Gomez
D.
,
Blaudeck
N.
,
Pfitzer
G.
.
2006
.
Mechanical properties of sarcomeres during cardiac myofibrillar relaxation: stretch-induced cross-bridge detachment contributes to early diastolic filling
.
J. Muscle Res. Cell Motil.
27
:
423
434
.
Sutko
J.L.
,
Airey
J.A.
,
Welch
W.
,
Ruest
L.
.
1997
.
The pharmacology of ryanodine and related compounds
.
Pharmacol. Rev.
49
:
53
98
.
Tada
H.
,
Higuchi
H.
,
Wanatabe
T.M.
,
Ohuchi
N.
.
2007
.
In vivo real-time tracking of single quantum dots conjugated with monoclonal anti-HER2 antibody in tumors of mice
.
Cancer Res.
67
:
1138
1144
.
Treiman
M.
,
Caspersen
C.
,
Christensen
S.B.
.
1998
.
A tool coming of age: thapsigargin as an inhibitor of sarco-endoplasmic reticulum Ca2+-ATPases
.
Trends Pharmacol. Sci.
19
:
131
135
.
Udaka
J.
,
Ohmori
S.
,
Terui
T.
,
Ohtsuki
I.
,
Ishiwata
S.
,
Kurihara
S.
,
Fukuda
N.
.
2008
.
Disuse-induced preferential loss of the giant protein titin depresses muscle performance via abnormal sarcomeric organization
.
J. Gen. Physiol.
131
:
33
41
.
Uemura
S.
,
Higuchi
H.
,
Olivares
A.O.
,
De La Cruz
E.M.
,
Ishiwata
S.
.
2004
.
Mechanochemical coupling of two substeps in a single myosin V motor
.
Nat. Struct. Mol. Biol.
11
:
877
883
.
Watanabe
T.M.
,
Fujii
F.
,
Jin
T.
,
Umemoto
E.
,
Miyasaka
M.
,
Fujita
H.
,
Yanagida
T.
.
2013
.
Four-dimensional spatial nanometry of single particles in living cells using polarized quantum rods
.
Biophys. J.
105
:
555
564
.
Werdich
A.A.
,
Lima
E.A.
,
Dzhura
I.
,
Singh
M.V.
,
Li
J.
,
Anderson
M.E.
,
Baudenbacher
F.J.
.
2008
.
Differential effects of phospholamban and Ca2+/calmodulin-dependent kinase II on [Ca2+]i transients in cardiac myocytes at physiological stimulation frequencies
.
Am. J. Physiol. Heart Circ. Physiol.
294
:
H2352
H2362
.
Yoo
J.
,
Kambara
T.
,
Gonda
K.
,
Higuchi
H.
.
2008
.
Intracellular imaging of targeted proteins labeled with quantum dots
.
Exp. Cell Res.
314
:
3563
3569
.

Abbreviations used in this paper:
FI

fluorescence intensity

OM

omecamtiv mecarbil

QD

quantum dot

SHG

second harmonic generation

SL

sarcomere length

SPOC

spontaneous sarcomeric oscillation

This article is distributed under the terms of an Attribution–Noncommercial–Share Alike–No Mirror Sites license for the first six months after the publication date (see http://www.rupress.org/terms). After six months it is available under a Creative Commons License (Attribution–Noncommercial–Share Alike 3.0 Unported license, as described at http://creativecommons.org/licenses/by-nc-sa/3.0/).

Supplementary data