==== Front Res Sq ResearchSquare Research Square American Journal Experts 37398289 10.21203/rs.3.rs-3043911/v1 10.21203/rs.3.rs-3043911 preprint 1 Article Impact of stretch on sarcomere length variability in isolated fully relaxed rat cardiac myocytes Lookin Oleg Boulali Najlae Université de Montpellier, INSERM, CNRS, CHU Arnaud de Villeneuve Cazorla Olivier Université de Montpellier, INSERM, CNRS, CHU Arnaud de Villeneuve Tombe Pieter Université de Montpellier, INSERM, CNRS, CHU Arnaud de Villeneuve AUTHOR CONTRIBUTIONS OL, NB, OC and PdT contributed to the conception of the study, design of experiments, analysis, and interpretation of the results. NB, OC and PdT contributed to the experimental measurements which were performed at the Montpellier laboratories. OL designed the software for data processing. The manuscript was written by OL, OC and PdT. All authors approved the final version of the manuscript. ✉ pdetombe@gmail.com 14 6 2023 rs.3.rs-3043911https://creativecommons.org/licenses/by/4.0/ This work is licensed under a Creative Commons Attribution 4.0 International License, which allows reusers to distribute, remix, adapt, and build upon the material in any medium or format, so long as attribution is given to the creator. The license allows for commercial use. nihpp-rs3043911v1.pdf The contractility of cardiac muscle is greatly affected by preload via the Frank-Starling Mechanism (FSM). It is based on the preload-dependent activation of sarcomeres – the elementary contractile units in muscle cells. Recent findings show a natural variability in sarcomere length (SL) in resting cardiomyocytes that, moreover, is altered in an actively contracting myocyte. SL variability may contribute to the FSM but it remains unresolved whether the change in the SL variability is regulated by activation process per se or simply by changes in cell stretch, i.e. average SL. To separate the roles of activation and SL, we characterized SL variability in isolated fully relaxed rat ventricular cardiomyocytes (n = 12) subjected to a longitudinal stretch with the carbon fiber (CF) technique. Each cell was tested in three states: without CF attachment (control, no preload), with CF attachment without stretch, and with CF attachment and ~ 10% stretch of initial SL. The cells were imaged by transmitted light microscopy to retrieve and analyze individual SL and SL variability off-line by multiple quantitative measures like coefficient of variation or median absolute deviation. We found that CF attachment without stretch did not affect the extent of SL variability and averaged SL. In stretched myocytes, the averaged SL significantly increased while the SL variability remained unchanged. This result clearly indicates that the non-uniformity of individual SL is not sensitive to the average SL itself in fully relaxed myocytes. We conclude that SL variability per se does not contribute to the FSM in the heart. isolated cardiomyocyte sarcomere sarcomere length variability carbon fibers stretch Centre National de la Recherche Scientifique#IEA00401 Agence Nationale pour la RechercheANR-21-CE14-0058 NIHHL62426 ==== Body pmcINTRODUCTION Sarcomere length (SL) is the fundamental parameter that controls preload-dependent regulation of muscle contraction [6, 7, 14]. The sarcomeric contractile response in cardiac myocyte non-linearly increases in proportion to its length over the physiological range [5]. The in-series arrangement of sarcomeres in myofibrils, therefore, strongly dictates that their individual lengths, simply due to their mechanical arrangement, induce significant mutual effects to each other during mechanical activation [6]. In resting skeletal and cardiac myocytes/myofibrils, SL heterogeneity may be as high as 20% of the average SL [1, 12, 15, 19, 20, 24, 32], and in one report even more [25]. Moreover, SL variability is increased in the actively contracting cardiomyocyte, i.e. in the systolic phase [17, 21] or during steady state activation in skeletal myocytes [3, 12, 15, 16, 25, 26, 27, 30]. However, it remains unresolved whether SL heterogeneity is governed by the (physiological) activation process, or whether it is caused simply by passive mechanical properties of the sarcomere, and thus can be “preadjusted” simply by prestretch (preloading) of the cell. In our previous study we have shown that the SL variability in an isolated mechanically unloaded and actively contracting cardiomyocyte is significantly higher in the end-systolic state at lower average SL vs the end-diastolic state at higher averaged SL [21]. The aim of the present study was to examine whether a physical change in SL itself directly affects SL variability in the absence of the activation process in a cell. We recorded the sarcomere striation signal in transmission light mode with enhanced contrast and high magnification using a confocal microscopy system. Each cell was imaged in the following states: without carbon fiber (CF) attachment (“unattached”), with CF attachment but without stretch (“attached”), and with stretch by CF (“stretched”). The extent of SL variability in each state was thoroughly characterized by multiple statistical measures, both absolute, such as coefficient of variation, inter-quartile range, median absolute deviation, and relative measures. METHODS All investigations conformed to European Parliament Directive 2010/63/EU and approved by the ethics committee “Comité d’éthique pour l’expérimentation animale Languedoc-Roussillon”. Isolated cardiomyocytes Ventricular cardiomyocytes were isolated by enzymatic digestion as previously described [2]. Briefly, male rats (~ 300 g) received a lethal dose of pentobarbital (50 mg/kg). The heart was quickly cannulated via the aorta and perfused with a HEPES physiological solution at 37°C containing type IV collagenase enzyme (1.25 mg/ml, Worthington). After enzymatic digestion, the external Ca2+ concentration was set to 0.1 mM. All subsequent experiments were performed using a modified Tyrode solution (140 mM NaCl, 4 mM KCl, 1 mM MgCl2, 20 mM HEPES, 0.1 mM CaCl2, and 11 mM glucose, pH 7.4) at room temperature without electrical pacing. All solutions contained 20 mM BDM to prevent spontaneous contractile activity. Optical measurements The analysis of intracellular variability in sarcomere lengths in isolated cardiomyocytes was performed using transmission light (TL) measurements. Briefly, myocytes were attached on both ends to 7 μm carbon rods using Myotak according to the previously described methods [8]. The CFs were mounted to hydraulic XYZ manipulators that aided in the carbon rod attachment, as well as gentle SL stretch of the attached myocytes. The cells were mounted in a glass bottomed 35mm plastic petri dish placed on top of an inverted microscope (Nikon, TE3000 Eclipse) equipped with a Nikon 60/1.2 water immersion objective. Transmitted light DIC images were obtained with a confocal microscope (Biorad Radiance 2100) using a 488 nm DPSS laser (Melles-Griot), as a 1024×1024 image with ~ 200 nm spatial resolution in the X-Y plane. Representative images of a cardiomyocyte in the three states – unattached to CFs, attached to CFs but not stretched, and stretched by CFs – are shown in Fig. 1. The stretch was approximately 10% of cell length in the resting state without CF attachment (SL stretched from ~ 1.9 μm to ~ 2.2 μm). Measurements were performed on n = 12 cardiomyocytes obtained from n = 7 hearts. Analysis of variability in sarcomere lengths The analysis of variability in individual sarcomere lengths was implemented by custom-made software EqapAll6 (IDE Delphi Embarcadero). Full-frame transmitted light images were taken in “unattached”, “attached”, and “stretched” states, each with five ROIs to cover as much area of a cell as possible. The lengths of individual sarcomeres were determined using an algorithm based on the analysis of the sarcomere striation profile in each ROI as input data (Fig. 2A). In brief, the algorithm works as follows: 1) input data are analyzed for the presence of non-periodic changes in intensity (typically due to the uneven brightness or image artifacts) which were then removed using low-frequency filtering, 2) the corrected input data were next filtered by moving summing up of the values and the resulting signal is leveled and scaled according to the initial range of intensity (Fig. 2A, bottom part of the panel), 3) next, the adjusted data were analyzed for the presence of local maxima by finding the spatial position of sign reversion of the first derivative of the signal, 4) when that position was found, the local maximum was calculated by parabolic approximation of three neighboring positions of the signal, where the central position is at the point of sign reversion of the first derivative. The X-axis positions of local maxima were then interpreted as Z-disks and were used for calculation of the individual SL as a distance between two neighboring Z-disks (Fig. 2B). The use of a moving sum to smooth the striation profile allowed for more accurate determination of the local maxima by rejecting false, noise- or artifact-related, local peaks; this approach also resulted in much denser SL distribution (Fig. 2C). In addition, the calculated SL values from all 5 ROIs were combined into a single individual SL data set and only values between 1.2 and 2.6 μm were used in further analysis of SL variability. Statistical analysis Statistical analyses were carried out with GraphPad Prism 7.0 (GraphPad Software). To verify that the samples of individual SL data do or do not follow a Gaussian distribution, both D’Agostino-Pearson and Shapiro-Wilk normality tests were applied. Thus, these normality tests were applied to the mean/median SL values and the values of corresponding SL variability measures, all treated as samples. The evaluation of significance in the mean/median SL values as well as the absolute/relative measures of SL variability between “unattached”, “attached” and “stretched” states was implemented by the Friedman test with Dunn’s multiple comparisons test that is appropriate for samples with a non-normal shape of distribution. Differences were considered significant at p < 0.05. Data are presented as mean ± S.D. RESULTS We evaluated whether stretch itself can affect SL variability in the non-stimulated fully relaxed cardiomyocyte in the following states: without CF attachment (“unattached”), with CF attachment but without stretch (“attached”), and with stretch by CF (“stretched”). Only cells measured in all three states were used for the analysis. First, we analyzed how many sets of individual SL values followed a normal distribution and whether this was affected by CF attachment or stretch. Out of total 36 analyzed SL data sets (12 sets per each state), 7 sets (19%) and 6 sets (17%) were non-normal according to the Shapiro-Wilk and D’Agostino & Pearson normality tests, respectively. The majority of SL data sets that did not pass these normality tests were obtained for “unattached” or “attached” states (11 sets), compared to the “stretched” state (2 sets). Next, we analyzed quantitative statistical measures calculated for the individual SL data sets – mean, median, standard deviation, median absolute deviation, etc… – obtained from all cells (n = 12) for normality of their distribution. It was needed for proper selection of significance comparison test between each two states of the cells. Most of the statistical measures passed both D’Agostino & Pearson and Shapiro-Wilk normality tests but the following did not: inter-quartile range (IRQ), median absolute deviation (MAD), IQR-to-median SL value, and MAD-to-median SL value. We therefore used the Friedman test with Dunn’s multiple comparisons test to compare each statistical measure between two states of a cell, i.e. “unattached” vs “attached”, “unattached” vs “stretched”, “attached” vs “stretched”, assuming non-normal distribution of the samples. The mean SL values were as follows: 1.94 ± 0.05 μm in the “unattached” state, 1.94 ± 0.04 μm in the “attached” state, and 2.15 ± 0.07 μm in the “stretched” state. A significant difference in mean SL was found for “stretched” state vs either “unattached” or “attached” state (p = 0.0001 and p = 0.0033, respectively), but not between the two latter states (Fig. 3A). Likewise, median SL values were found to be 1.94 ± 0.06 μm in the “unattached” state, 1.95 ± 0.05 μm in the “attached” state, and 2.20 ± 0.11 μm in “stretched” state, with a significant difference between “stretched” state vs either “unattached” or “attached” state (p = 0.0007 for both comparisons) but not between the two latter states (Fig. 3B). An important finding in non-stimulated fully relaxed isolated cardiomyocytes, imaged by the TL mode, was that the extent of SL variability, as evaluated by several absolute and relative indexes, was neither affected by the CF attachment itself, nor by cell stretch (Fig. 4A–E). This could also be demonstrated by direct comparison of SL distribution plots obtained for these states of the cells, where the width of the plots remains virtually identical regardless the state (Fig. 4F). Also, if we plotted the measure of SL variability (for example, MAD-to-median SL value) against the median SL value itself for all three states, we did not see a functional relationship between the groups of data (Fig. 5A). Finally, we calculated delta values between SL obtained in “unattached” and “attached”/”stretched” states, with further plotting of these delta SL values vs MAD-to-median SL values (Fig. 5B). The relationship for “attached” state was uncertain because the CF attachment did not produce any effect on the average SL. In the “stretched” state, however, there was also non significant relationship between the extent of SL variability and average SL value (e.g. Spearman r value was 0.52, p = 0.09). These findings indicate that SL variability is not sensitive to average sarcomere length per se. DISCUSSION In the present study, we tested the effect of increase in average SL, i.e. the direct effect of preload, on the variability in the individual SLs in rat ventricular cardiomyocytes held in a fully relaxed state. We found that neither the extent of SL variability nor the character of this variability is affected by the preload itself. Therefore, experimentally observed increase in the SL heterogeneity in activated cardiomyocytes [17, 21] or in skeletal muscle [4, 16, 25, 27] appears to be related to the inhomogeneous process, both spatially and in time, of Ca2+ activation of individual sarcomeres, but not to the actual change in average SL. However, our finding was obtained in fully relaxed cardiomyocytes, i.e. with low Ca2+ and BDM uncoupling between actin filament and myosin motors – a situation that may not be present under physiological conditions. Under normal conditions, a quiescent myocyte still experiences a certain level of background interaction between actin and myosin due to the non-disturbed myosin activity and a low but physiologically important diastolic level of [Ca2+]i [9, 22]. This may predetermine an initial extent of SL non-uniformity in a resting (but not fully relaxed) cardiomyocyte. Accordingly, in an activated cell the contractile response of individual sarcomeres may be substantially non-synchronous [18, 37], thereby contributing to increased SL inhomogeneity due to the in-series mechanical communication between strongly activated (contracting) and weakly activated individual sarcomeres that are lengthened by their contracting neighbors [13]. The phenomenon of SL non-uniformity has already been extensively studied in skeletal muscle [4, 11, 33]. SL heterogeneity is substantially increased even if only a single sarcomere is disrupted in the initially untreated myofibril [12] and this is accompanied by diminished myofibrillar active force production [23]. Importantly, SL heterogeneity is even more increased in a muscle that is subjected to mechanical stretch upon activation [29, 30, 31] – these results clearly indicates the importance of dynamic redistribution of individual SLs in activated muscle. In contrast to the multitude of studies in skeletal muscle, the direct effect of sarcomere pre-stretch on SL variability in contracting cardiac myocytes remains to be determined. To elucidate this aspect in detail, future studies on SL variability should include cardiomyocytes that start their contraction from different end-diastolic SL (i.e. varied preloads). A limitation of the present study is the use of non-confocal transmission light (TL) mode of measurements. The images taken in this mode are less accurate in terms of spatial determination of sarcomeres [17, 36]. An important note has been raised in [32], where the authors experimentally verified that the nanometer variations in SL, when obtained by conventional light microscopy, should not be interpreted as a real variation but rather as noise. They suggested reasonable and optimal levels of optical resolution – 100 and 50 nm, respectively – that are necessary to conform to a reliable quantification of SL variability. The use of advanced optical methods such as second-harmonic generation microscopy [10] may help to satisfy the requirements imposed by the recently introduced term “SL nanometry” [34, 35, 38]. Nevertheless, we recently reported hat SL variability measured by confocal microscopy in fluorescently stained (ANEPPS) non-stretched fully relaxed isolated rat cardiomyocytes was nearly the same by magnitude as for the TL mode in the current study (Fig. 2F in [20]). Similar results were also obtained in another study that combined confocal and TL modes of SL measurements [28]. Of note, our recent and current studies employed comparable spatial resolutions of the optical systems. We therefore conclude that, provided the optical settings provide proper spatial resolution, the TL imaging may be at least as good for SL analysis as confocal imaging. ACKNOWLEDGEMENTS We thank Patrice Bideaux for assistance with the rat cell isolation procedures. FUNDING The study was supported by grants from the Centre National de la Recherche Scientifique (CNRS, France, grant #IEA00401 to OC), Agence Nationale pour la Recherche (ANR-21-CE14-0058, NITROSOCARD to OC) and NIH (HL62426 to PdT). DATA AVAILABILITY Data acquired and analyzed during this study are available from the corresponding author upon reasonable request. Figure 1: Representative transmission light images of a cardiomyocyte in three states: “unattached” (no attached carbon fibers, CF), “attached” (CF attached but no stretch applied), and “stretched” (myocyte is stretched by CF). Five regions of interest were used in each image to collect sarcomeric data. Figure 2: The implementation of the custom-made algorithm designed for calculation of individual sarcomere lengths (SL) using a sarcomere striation pattern. (A) Example of transmitted light image of a cardiomyocyte with 5 regions of interests (ROIs) placed over the whole area of the cell. The sarcomere striation profile is obtained from the selected ROI. The algorithm uses raw data (black curves) as input and first removes non-regular components of the data, then it produces (via several steps, see text for details) the adjusted data with significantly decreased impurities in the regular striation pattern (red curved). The local maxima of the adjusted signal correspond to the spatial positions of sarcomeric Z-disks; the distance between two neighboring Z-disks is interpreted as an individual SL. (B) The calculated individual SLs as a function of their X-positions in the striation profile obtained from raw or algorithmically adjusted data. (C) The distributions of SL obtained from raw or algorithmically adjusted data. All data are shown for an individual exemplary cardiac cell. Figure 3: Stretch of resting cardiomyocytes affects sarcomere length (SL) mean and median values. The carbon fiber (CF) technique was used to test each myocyte in three states: without CF attachment (“unattached”), with CF attached but without stretch (“attached”), and with ~10% cell stretch applied (“stretched”). Measurements were made in transmitted light mode for n = 12 isolated cardiomyocytes from N = 7 hearts. (A) Mean SL value, and (B) Median SL value. * Difference between “stretched” state and either “unattached” or “attached” state is significant at p< 0.01 (Friedman’s test with Dunn’s multiple comparisons test). Figure 4: Stretch of resting cardiomyocyte does not affect the variability in individual sarcomere lengths (SL). The measurements were obtained in transmitted light mode for n = 12 isolated cardiomyocytes (each in “unattached”, “attached”, and “stretched” state) from N = 7 hearts. The following absolute and relative quantitative measures of SL variability were analyzed: (A) Coefficient of variation, (B) Inter-quartile range (IQR), (C) Median absolute deviation (MAD), (D) IQR divided by median SL value, and (E) MAD divided by median SL value. No significant differences were observed between “unattached”, “attached” and “stretched” states (except for the inter-quartile range divided by median SL value between “unattached” vs “attached”, indicated by “#” on panel D). (F) Average distribution plots for individual SL values obtained in “unattached”, “attached”, and “stretched” states. The peak of the plots corresponds to the averaged median SL value and the distributions are shown as a % of the corresponding median SL value. Figure 5: Lack of significant relationship between averaged sarcomere length (SL) and extent of SL variability. The measurements were obtained in transmitted light mode for n = 12 isolated cardiomyocytes (each in “unattached”, “attached”, and “stretched” state) from N = 7 hearts. (A) MAD-to-median SL value is plotted as function of average SL for all three states. (B) MAD-to-median SL value is plotted as function of absolute change in average SL (ΔSL) for “attached” and “stretched” states (ΔSL are calculated relatively to “unattached” state). Note that there is no significant relationship between ΔSL and the extent of SL variability (Spearman r= 0.52 with p = 0.09). Conflict of interest None. Competing interests The authors declare no competing interests. ==== Refs References 1. Adkins AN , Fong RM , Dewald JPA , Murray WM (2022) Variability of in vivo sarcomere length measures in the upper limb obtained with second harmonic generation microendoscopy. Front Physiol 12 :817334. 10.3389/fphys.2021.817334 35211028 2. Andre L , Boissière J , Reboul C , Perrier R , Zalvidea S , Meyer G , Thireau J , Tanguy S , Bideaux P , Hayot M , Boucher F , Obert P , Cazorla O , Richard S (2010) Carbon monoxide pollution aggravates ischemic heart failure through oxidative stress pathway. Am J Respir Crit Care Med 181 (6 ):587–595. 10.1164/rccm.200905-0794OC 20019346 3. de Souza Leite F , Minozzo FC , Altman D , Rassier DE (2017) Microfluidic perfusion shows intersarcomere dynamics within single skeletal muscle myofibrils. Proc Natl Acad Sci U S A 114 (33 ):8794–8799. 10.1073/pnas.1700615114 28765372 4. de Souza Leite F , Rassier DE (2020) Sarcomere length nonuniformity and force regulation in myofibrils and sarcomeres. Biophys J 119 :1–6. 10.1016/j.bpj.2020.11.005 32521239 5. de Tombe PP , Mateja RD , Tachampa K , Ait Mou Y , Farman GP , Irving TC (2010) Myofilament length dependent activation. J Mol Cell Cardiol 48 :851–858. 10.1016/j.yjmcc.2009.12.017 20053351 6. de Tombe PP , ter Keurs HEDJ (2016) Cardiac muscle mechanics: sarcomere length matters. J Mol Cell Cardiol 91 :148–150. 10.1016/j.yjmcc.2015.12.006 26678623 7. Dobesh DP , Konhilas JP , de Tombe PP (2002) Cooperative activation in cardiac muscle: impact of sarcomere length. Am J Physiol Heart Circ Physiol 282 :H1055–H1062. 10.1152/ajpheart.00667.2001 11834504 8. Dvornikov AV , Dewan S , Alekhina OV , Pickett FB , de Tombe PP (2014) Novel approaches to determine contractile function of the isolated adult zebrafish ventricular cardiac myocyte. J Physiol 592 (9 ):1949–1956. 10.1113/jphysiol.2014.270678 24591576 9. Eisner DA , Caldwell JL , Trafford AW , Hutchings DC (2020) The control of diastolic calcium in the heart. Basic mechanisms and functional implications. Circ Res 126 : 395–412. 10.1161/CIRCRESAHA.119.315891 31999537 10. Fujita H , Kaneshiro J , Takeda M , Sasaki K , Yamamoto R , Umetsu D , Kuranaga E , Higo S , Kondo T , Asano Y , Sakata Y , Miyagawa S , Watanabe TM (2023) Estimation of crossbridge-state during cardiomyocyte beating using second harmonic generation. Life Sci Alliance 6 (7 ):e202302070. 10.26508/lsa.202302070 37236659 11. Haeger R , de Souza Leite F , Rassier DE (2020) Sarcomere length nonuniformities dictate force production along the descending limb of the force–length relation. Proc R Soc B 287 :20202133. 10.1098/rspb.2020.2133 12. Haeger RM , Rassier DE (2020) Force enhancement after stretch of isolated myofibrils is increased by sarcomere length non-uniformities. Sci Rep 10 :21590. 10.1038/s41598-020-78457-1 33299041 13. Helmes M , Palmer BM (2022) Sarcomere length in the beating heart: Synchronicity is optional. J Gen Physiol 154 (2 ):e202113022. 10.1085/jgp.202113022 35015810 14. Herzog W (2022) What can we learn from single sarcomere and myofibril preparations? Front Physiol 13 :837611. 10.3389/fphys.2022.837611 35574477 15. Johnston K , Jinha A , Herzog W (2016) The role of sarcomere length non-uniformities in residual force enhancement of skeletal muscle myofibrils. R Soc Open Sci 3 :150657. 10.1098/rsos.150657 27069655 16. Johnston K , Moo EK , Jinha A , Herzog W (2019) On sarcomere length stability during isometric contractions before and after active stretching. J Exp Biol 222 :jeb209924. 10.1242/jeb.209924 31704896 17. Kobirumaki-Shimozawa F , Oyama K , Shimozawa T , Mizuno A , Ohki T , Terui T , Minamisawa S , Ishiwata S , Fukuda N (2016) Nano-imaging of the beating mouse heart in vivo: Importance of sarcomere dynamics, as opposed to sarcomere length per se, in the regulation of cardiac function. J Gen Physiol 147 (1 ):53–62. 10.1085/jgp.201511484 26712849 18. Kobirumaki-Shimozawa F , Shimozawa T , Oyama K , Baba S , Li J , Nakanishi T , Terui T , Louch WE , Ishiwata S , Fukuda N (2021) Synchrony of sarcomeric movement regulates left ventricular pump function in the in vivo beating mouse heart. J Gen Physiol 153 (11 ):e202012860. 10.1085/jgp.202012860 34605861 19. Llewellyn ME , Barretto RPJ , Delp SL , Schnitzer MJ (2008) Minimally invasive high-speed imaging of sarcomere contractile dynamics in mice and humans. Nature 454 (7205 ):784–788. 10.1038/nature07104 18600262 20. Lookin O , de Tombe P , Boulali N , Gergely C , Cloitre T , Cazorla O (2023) Cardiomyocyte sarcomere length variability: Membrane fluorescence versus second harmonic generation myosin imaging. J Gen Physiol 155 (4 ):e202213289. 10.1085/jgp.202213289 36695814 21. Lookin O , Khokhlova A , Myachina T , Butova X , Cazorla O , de Tombe P (2022) Contractile state dependent sarcomere length variability in isolated guinea-pig cardiomyocytes. Front Physiol 13 :857471. 10.3389/fphys.2022.857471 35444559 22. Louch WE , Stokke MK , Sjaastad I , Christensen G , Sejersted OM (2012) No rest for the weary: diastolic calcium homeostasis in the normal and failing myocardium. Physiology (Bethesda) 27 (5 ):308–323. 10.1152/physiol.00021.2012 23026754 23. Mendoza AC , Rassier DE (2020) Extraction of thick filaments in individual sarcomeres affects force production by single myofibrils. Biophys J 118 (8 ):1921–1929. 10.1016/j.bpj.2020.03.007 32251620 24. Moo EK , Fortuna R , Sibole SC , Abusara Z , Herzog W (2016) In vivo sarcomere lengths and sarcomere elongations are not uniform across an intact muscle. Front Physiol 7 :187. 10.3389/fphys.2016.00187 27252660 25. Moo EK , Herzog W (2018) Single sarcomere contraction dynamics in a whole muscle. Sci Rep 8 :15235. 10.1038/s41598-018-33658-7 30323321 26. Moo EK , Herzog W (2020) Sarcomere lengths become more uniform over time in intact muscle-tendon unit during isometric contractions. Front Physiol 11 :448. 10.3389/fphys.2020.00448 32477162 27. Moo EK , Leonard TR , Herzog W (2017) In vivo sarcomere lengths become more non-uniform upon activation in intact whole muscle. Front Physiol 8 :1015. 10.3389/fphys.2017.01015 29270135 28. Nance ME , Whitfield JT , Zhu Y , Gibson AK , Hanft LM , Campbell KS , Meininger GA , McDonald KS , Segal SS , Domeier TL (2015) Attenuated sarcomere lengthening of the aged murine left ventricle observed using two-photon fluorescence microscopy. Am J Physiol Heart Circ Physiol 309 :H918–H925. 10.1152/ajpheart.00315.2015 26209054 29. Rassier DE (2012) Residual force enhancement in skeletal muscles: one sarcomere after the other. J Muscle Res Cell Motil 33 :155–165. 10.1007/s10974-012-9308-7 22729612 30. Rassier DE (2017) Sarcomere mechanics in striated muscles: from molecules to sarcomeres to cells. Am J Physiol Cell Physiol 313 (2 ):C134–C145. 10.1152/ajpcell.00050.2017 28539306 31. Rassier DE , Pavlov I (2012) Force produced by isolated sarcomeres and half-sarcomeres after an imposed stretch. Am J Physiol Cell Physiol 302 :C240–C248. 10.1152/ajpcell.00208.2011 21998143 32. Schmidt J , Jinha A , Herzog W (2021) Sarcomere length measurement reliability in single myofibrils. J Biomech 126 :110628. 10.1016/j.jbiomech.2021.110628 34274869 33. Shimamoto Y , Suzuki M , Mikhailenko SV , Yasuda K , Ishiwata S (2009) Inter-sarcomere coordination in muscle revealed through individual sarcomere response to quick stretch. Proc Natl Acad Sci U S A 106 (29 ):11954–9. 10.1073/pnas.0813288106 19515816 34. Shimozawa T , Hirokawa E , Kobirumaki-Shimozawa F , Oyama K , Shintani SA , Terui T , Kushida Y , Tsukamoto S , Fujii T , Ishiwata S , Fukuda N (2017) In vivo cardiac nano-imaging: A new technology for high-precision analyses of sarcomere dynamics in the heart. Prog Biophys Mol Biol 124 :31–40. 10.1016/j.pbiomolbio.2016.09.006 27664770 35. Shintani SA , Oyama K , Kobirumaki-Shimozawa F , Ohki T , Ishiwata S , Fukuda N (2014) Sarcomere length nanometry in rat neonatal cardiomyocytes expressed with a-actinin-AcGFP in Z discs. J Gen Physiol 143 :513–524. 10.1085/jgp.201311118 24638993 36. Telley IA , Denoth J , Stussi E , Pfitzer G , Stehle R (2006) Half-sarcomere dynamics in myofibrils during activation and relaxation studied by tracking fluorescent markers. Biophys J 90 :514–530. 10.1529/biophysj.105.070334 16239326 37. Tsukamoto S , Fujii T , Oyama K , Shintani SA , Shimozawa T , Kobirumaki-Shimozawa F , Ishiwata S , Fukuda N (2016) Simultaneous imaging of local calcium and single sarcomere length in rat neonatal cardiomyocytes using yellow Cameleon-Nano140. J Gen Physiol 148 (4 ):341–355. 10.1085/jgp.201611604 27670899 38. Young KW , Kuo B P-P , O’Connor SM , Radic S , Lieber RL (2017) In vivo sarcomere length measurement in whole muscles during passive stretch and twitch contractions. Biophys J 112 (4 ):805–812. 10.1016/j.bpj.2016.12.046 28256239