<?xml version="1.0" ?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.3 20210610//EN"  "JATS-archivearticle1-mathml3.dtd"><article xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="1.3" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">elife</journal-id>
<journal-id journal-id-type="publisher-id">eLife</journal-id>
<journal-title-group>
<journal-title>eLife</journal-title>
</journal-title-group>
<issn publication-format="electronic" pub-type="epub">2050-084X</issn>
<publisher>
<publisher-name>eLife Sciences Publications, Ltd</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">87065</article-id>
<article-id pub-id-type="doi">10.7554/eLife.87065</article-id>
<article-id pub-id-type="doi" specific-use="version">10.7554/eLife.87065.2</article-id>
<article-version-alternatives>
<article-version article-version-type="publication-state">reviewed preprint</article-version>
<article-version article-version-type="preprint-version">1.2</article-version>
</article-version-alternatives>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Cell Biology</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Independent regulation of Z-lines and M-lines during sarcomere assembly in cardiac myocytes revealed by the automatic image analysis software sarcApp</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Neininger-Castro</surname>
<given-names>Abigail C.</given-names>
</name>
<degrees>PhD</degrees>
<xref ref-type="aff" rid="a1">1</xref>
<xref ref-type="corresp" rid="cor1">*</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Hayes</surname>
<given-names>James B.</given-names>
<suffix>Jr.</suffix></name>
<xref ref-type="aff" rid="a1">1</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Sanchez</surname>
<given-names>Zachary C.</given-names>
</name>
<xref ref-type="aff" rid="a1">1</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Taneja</surname>
<given-names>Nilay</given-names>
</name>
<degrees>PhD</degrees>
<xref ref-type="aff" rid="a1">1</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fenix</surname>
<given-names>Aidan M.</given-names>
</name>
<degrees>PhD</degrees>
<xref ref-type="aff" rid="a1">1</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Moparthi</surname>
<given-names>Satish</given-names>
</name>
<degrees>PhD</degrees>
<xref ref-type="aff" rid="a2">2</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Vassilopoulos</surname>
<given-names>Stéphane</given-names>
</name>
<degrees>PhD</degrees>
<xref ref-type="aff" rid="a2">2</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Burnette</surname>
<given-names>Dylan T.</given-names>
</name>
<degrees>PhD</degrees>
<xref ref-type="aff" rid="a1">1</xref>
<xref ref-type="corresp" rid="cor1">*</xref>
</contrib>
<aff id="a1"><label>1</label><institution>Department of Cell and Developmental Biology, Vanderbilt University School of Medicine Basic Sciences</institution>, Nashville, TN</aff>
<aff id="a2"><label>2</label><institution>Sorbonne Université, INSERM, Institut de Myologie, Centre de Recherche en Myologie</institution>, Paris, <country>France</country></aff>
</contrib-group>
<contrib-group content-type="section">
<contrib contrib-type="editor">
<name>
<surname>Dubois</surname>
<given-names>Nicole C</given-names>
</name>
<role>Reviewing Editor</role>
<aff>
<institution-wrap>
<institution>Icahn School of Medicine at Mount Sinai</institution>
</institution-wrap>
<city>New York</city>
<country>United States of America</country>
</aff>
</contrib>
<contrib contrib-type="senior_editor">
<name>
<surname>Stainier</surname>
<given-names>Didier YR</given-names>
</name>
<role>Senior Editor</role>
<aff>
<institution-wrap>
<institution>Max Planck Institute for Heart and Lung Research</institution>
</institution-wrap>
<city>Bad Nauheim</city>
<country>Germany</country>
</aff>
</contrib>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>*</label>To whom correspondence should be addressed: <email>dylan.burnette@vanderbilt.edu</email> and <email>abbieneininger@gmail.com</email></corresp>
</author-notes>
<pub-date date-type="original-publication" iso-8601-date="2023-05-15">
<day>15</day>
<month>05</month>
<year>2023</year>
</pub-date>
<pub-date date-type="update" iso-8601-date="2023-09-04">
<day>04</day>
<month>09</month>
<year>2023</year>
</pub-date>
<volume>12</volume>
<elocation-id>RP87065</elocation-id>
<history>
<date date-type="sent-for-review" iso-8601-date="2023-03-02">
<day>02</day>
<month>03</month>
<year>2023</year>
</date>
</history>
<pub-history>
<event>
<event-desc>Preprint posted</event-desc>
<date date-type="preprint" iso-8601-date="2023-01-12">
<day>12</day>
<month>01</month>
<year>2023</year>
</date>
<self-uri content-type="preprint" xlink:href="https://doi.org/10.1101/2023.01.11.523681"/>
</event>
<event>
<event-desc>Reviewed preprint v1</event-desc>
<date date-type="reviewed-preprint" iso-8601-date="2023-05-15">
<day>15</day>
<month>05</month>
<year>2023</year>
</date>
<self-uri content-type="reviewed-preprint" xlink:href="https://doi.org/10.7554/eLife.87065.1"/>
<self-uri content-type="editor-report" xlink:href="https://doi.org/10.7554/eLife.87065.1.sa3">eLife assessment</self-uri>
<self-uri content-type="referee-report" xlink:href="https://doi.org/10.7554/eLife.87065.1.sa2">Reviewer #1 (Public Review):</self-uri>
<self-uri content-type="referee-report" xlink:href="https://doi.org/10.7554/eLife.87065.1.sa1">Reviewer #2 (Public Review):</self-uri>
<self-uri content-type="referee-report" xlink:href="https://doi.org/10.7554/eLife.87065.1.sa0">Reviewer #3 (Public Review):</self-uri>
</event>
</pub-history>
<permissions>
<copyright-statement>© 2023, Neininger-Castro et al</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Neininger-Castro et al</copyright-holder>
<ali:free_to_read/>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<ali:license_ref>https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="elife-preprint-87065-v2.pdf"/>
<abstract>
<title>Abstract</title>
<p>Sarcomeres are the basic contractile units within cardiac myocytes, and the collective shortening of sarcomeres aligned along myofibrils generates the force driving the heartbeat. The alignment of the individual sarcomeres is important for proper force generation, and misaligned sarcomeres are associated with diseases including cardiomyopathies and COVID-19. The actin bundling protein, α-actinin-2, localizes to the “Z-Bodies” of sarcomere precursors and the “Z-Lines” of sarcomeres, and has been used previously to assess sarcomere assembly and maintenance. Previous measurements of α-actinin-2 organization have been largely accomplished manually, which is time-consuming and has hampered research progress. Here, we introduce sarcApp, an image analysis tool that quantifies several components of the cardiac sarcomere and their alignment in muscle cells and tissue. We first developed sarcApp to utilize deep learning-based segmentation and real space quantification to measure α-actinin-2 structures and determine the organization of both precursors and sarcomeres/myofibrils. We then expanded sarcApp to analyze “M-Lines” using the localization of myomesin and a protein that connects the Z-Lines to the M-Line (titin). sarcApp produces 33 distinct measurements per cell and 24 per myofibril that allow for precise quantification of changes in sarcomeres, myofibrils, and their precursors. We validated this system with perturbations to sarcomere assembly. We found perturbations that affected Z-Lines and M-Lines differently, suggesting that they may be regulated independently during sarcomere assembly.</p>
</abstract>

</article-meta>
<notes>
<notes notes-type="competing-interest-statement">
<title>Competing Interest Statement</title><p>The authors have declared no competing interest.</p></notes>
<fn-group content-type="summary-of-updates">
<title>Summary of Updates:</title>
<fn fn-type="update"><p>This revision includes suggested changes by reviewers for eLife.</p></fn>
</fn-group>
<fn-group content-type="external-links">
<fn fn-type="dataset"><p>
<ext-link ext-link-type="uri" xlink:href="https://github.com/abbieneininger/sarcApp">https://github.com/abbieneininger/sarcApp</ext-link>
</p></fn>
</fn-group>
</notes>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>The sarcomere is the fundamental unit of contraction in striated muscle<sup><xref ref-type="bibr" rid="c1">1</xref></sup>. Sarcomeres individually contract to generate a force, and sarcomeres organize and align within myofibrils to augment that force. These forces enable skeletal muscle to coordinate organism-level locomotion and cardiac muscle to drive the heartbeat<sup><xref ref-type="bibr" rid="c2">2</xref>,<xref ref-type="bibr" rid="c3">3</xref></sup>. Sarcomere structure and alignment is critical to muscle function, and dysfunctional sarcomeres and myofibrils have been implicated as causes of disease<sup><xref ref-type="bibr" rid="c4">4</xref>-<xref ref-type="bibr" rid="c6">6</xref></sup>.</p>
<p>Since the mid-20<sup>th</sup> century, the question of how sarcomeres are assembled and subsequently organized into higher-order myofibrils has been a much-debated topic, with multiple separate lines of research emerging to support a small host of potential models<sup><xref ref-type="bibr" rid="c7">7</xref>-<xref ref-type="bibr" rid="c9">9</xref></sup>. Generally, the models of sarcomere assembly that are backed by the most data posit that sarcomeres arise from either i) a “parts-wise” assembly, where precursors of unique sarcomeric regions (e.g., the thin (actin) or thick (myosin) filament) are first assembled separately, then “stitched” together to generate the final structure (i.e., Stitching Model) or ii) a more direct assembly, where a precursor forms that resembles a sarcomere but has a subset of non-sarcomeric proteins that are later replaced by sarcomeric proteins (i.e., Pre-myofibril Model)<sup><xref ref-type="bibr" rid="c7">7</xref>,<xref ref-type="bibr" rid="c10">10</xref></sup>. Recently, our group has shown that sarcomeres descend directly from <underline>m</underline>uscle <underline>s</underline>tress fibers (MSFs), with our data supporting a subset of predictions made by both models, but not necessarily favoring either model<sup><xref ref-type="bibr" rid="c11">11</xref></sup>. All models of sarcomere assembly are currently limited in their testability because the field lacks a tractable workflow for high throughput analysis.</p>
<p>Recently, our lab has shown that human induced pluripotent stem cell-derived cardiac myocytes (hiCMs) can be harnessed as a tractable model for routine imaging of sarcomere assembly in live CMs<sup><xref ref-type="bibr" rid="c11">11</xref>,<xref ref-type="bibr" rid="c12">12</xref></sup>. In our assay, hiCMs assemble sarcomeres de novo, and directly from MSFs, within 24 hours. Despite having powerful tools for imaging, our analysis workflow is still hampered by technical bottlenecks related to i) image binarization and ii) manual quantification. Binarization via classical methods relying on image pixel intensity (e.g., Otsu’s method) is often performed manually—slowly—and can become muddled by a multitude of factors, including i) high-intensity artifacts related to nonspecific localization of antibodies, ii) low-intensity staining due to high background levels/noise, and/or iii) multi-valent interactions by the target protein itself occluding the specific structure-of-interest<sup><xref ref-type="bibr" rid="c13">13</xref>-<xref ref-type="bibr" rid="c15">15</xref></sup>. Meanwhile, we find that even well-binarized, high-resolution images can require several hours of manual quantification per image.</p>
<p>Deep learning has emerged as a powerful potential solution for rapid, accurate binarization of complex greyscale images. Deep learning-based binarization can be accomplished using a framework called a U-Net<sup><xref ref-type="bibr" rid="c16">16</xref></sup>. A U-Net trains a deep learning model to convert images into only background (noise) and foreground (signal) by matching images to expert-annotated “ground truth” binaries – i.e., images manually annotated by the user. Compared to manually performed binarization, U-Net offers the advantages of superior accuracy and speed, but at a high upfront cost, since implementation of a U-Net requires a user with advanced background in mathematics, experience with machine learning, and proficiency in writing software code<sup><xref ref-type="bibr" rid="c16">16</xref>,<xref ref-type="bibr" rid="c17">17</xref></sup>.</p>
<p>While a sarcomere contains 100+ proteins, the field has often resorted to “proxying” sarcomere assembly by staining, often exclusively, for α-actinin-2<sup><xref ref-type="bibr" rid="c11">11</xref>,<xref ref-type="bibr" rid="c18">18</xref>-<xref ref-type="bibr" rid="c22">22</xref></sup>. α-actinin-2 marks the sarcomeric Z-Lines, which border the sarcomeric contractile machinery<sup><xref ref-type="bibr" rid="c23">23</xref>,<xref ref-type="bibr" rid="c24">24</xref></sup>. In the last decade, several methods for automatic quantification of Z-Lines have emerged, with the most widely used methods relying on calculations made in frequency space (i.e., by using a Fourier transformation of the image)<sup><xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c21">21</xref></sup>. While frequency space calculations can detect repeating patterns created by adjacent Z-Lines, current methods yield no information about Z-Lines in real space (e.g., number, size, spacing, or organization within myofibrils), about sarcomere precursors, or about any sarcomeric components other than Z-Lines. Moreover, it is not yet clear quite how the assembly of Z-Lines themselves relates either spatially or temporally to the assembly of other sarcomere components, or vice versa.</p>
<p>Here, we present a high throughput, automated, and non-biased quantification scheme for sarcomere assembly. Our approach is two-pronged and includes a unique method for both 1) rapid, accurate binarization of images that can consistently separate true signal from noise and structures-of-interest from those that are not-of-interest and 2) automated quantification of relevant sarcomeric and pre-sarcomeric structures within the binarized images. Towards 1), we have developed “yoU-Net”, a U-Net-based framework with a user-friendly graphical user interface (GUI) that enables a user with little to no experience to binarize images with the power and flexibility afforded by deep learning. We demonstrate that yoU-Net can accurately binarize images of CM sarcomeres and sarcomere precursors using multiple unique sarcomere stains. Towards 2), we have developed “sarcApp”, a software that automatically annotates and calculates 50+ unique descriptive outputs of sarcomeres and sarcomere precursors in real space using binarized images. Using yoU-Net and sarcApp, we demonstrate that some perturbations to sarcomere assembly preferentially affect specific sarcomere components (e.g., M-Lines) more than others (e.g., Z-Lines). Our data emphasize the need for more high throughput studies into sarcomere assembly that includes multiple sarcomere components, and altogether point to a disjointed or modular sarcomere assembly program within CMs that is more nuanced than was previously believed. The tools presented herein will facilitate future high throughput studies to develop a more mature and complex model of sarcomere assembly.</p>
</sec>
<sec id="s2">
<title>Results</title>
<sec id="s2a">
<title>yoU-Net enables user-friendly binarization of dynamic grayscale images</title>
<p>Our overarching goal was to develop a method for fast and accurate quantification of sarcomeres within hiCMs using our sarcomere assembly assay. We faced two primary bottlenecks in our workflow, the first of which was a slow and cumbersome image binarization process. Binarization facilitates image quantification by converting the complex, raw image of a biological specimen into a far simpler image containing only background (non-relevant information) and signal (relevant information). A well-binarized image should retain high fidelity to the original image but contain only the relevant structures of interest superimposed over a blank background. Individual pixels within a binary image are stripped of dynamic information present within the raw image and are re-assigned to represent only either background or signal – binary images are therefore often represented in black (background) and white (signal).</p>
<p>Traditional methods of binarization fall short for images of sarcomeres and sarcomere precursors in hiCMs due to the dynamic range of pixel intensities present within the raw images. In images of hiCMs 24 hours after plating, sarcomere precursors and sarcomeres are each stained by α-actinin-2 (<xref rid="fig1" ref-type="fig">Figure 1A</xref>), but the Z-Lines of the sarcomeres are markedly brighter than the Z-Bodies of the precursor structures (i.e., muscle stress fibers (MSFs)). Classical, intensity-based binarization methods force a user to make a choice in this scenario – specifically, whether to accurately binarize based on gray-levels of the sarcomere precursors or of the sarcomeres themselves (without the consistent option for both). Classical Otsu-based thresholding (using FIJI) of a representative cell in <xref rid="fig1" ref-type="fig">Figure 1A</xref> shows that binarizing this cell according to Z-Lines oversaturates Z-Bodies into an inseparable clump (<xref rid="fig1" ref-type="fig">Figure 1A</xref>, arrowhead), leaving a user unable to accurately quantify both structures from this single binary.</p>
<fig id="fig1" position="float" fig-type="figure">
<label>Figure 1:</label>
<caption><title>yoU-Net follows a U-Net architecture to binarize immunofluorescence images.</title>
<p>A) Representative images of α-actinin-2 and actin filaments (phalloidin) in a hiCM. The α-actinin-2 binary was thresholded in FIJI using Otsu’s method. Orange arrowhead denotes Z-Bodies. B) yoU-Net architecture from input image to output binary. Details can be found in the Supplement and in Figure S1. C) yoU-Net-generated binary of α-actinin-2, predicted using the trained U-Net described in <xref ref-type="fig" rid="fig1">Figure 1B</xref>. Orange arrowhead: Z-Bodies. D) Model of muscle stress fibers (MSFs) and myofibrils during sarcomere formation. Black arrow denotes direction of MSF translocation as α-actinin-2-positive Z-Bodies elongate and coalesce to form Z-Lines.</p></caption>
<graphic xlink:href="523681v2_fig1.tif" mimetype="image" mime-subtype="tiff"/>
</fig>
<p>Our lab was well-positioned to approach this problem by constructing a deep-learning-based model called a U-Net<sup><xref ref-type="bibr" rid="c16">16</xref></sup>. A U-Net enables a user to train a deep-learning model to automatically binarize images with high precision and speed. Model “training” involves matching raw images to user-annotated binaries – called “ground truth” binaries – using a complex set of mathematical convolution/downsampling and deconvolution/upsampling steps. The architectural layout of a U-Net is represented by <xref rid="fig1" ref-type="fig">Figure 1B</xref>.</p>
<p>We constructed our own U-Net and trained it using 1000 training steps (100 epochs/steps of 10 iterations each) of hiCMs stained for α-actinin-2. After training, our U-Net can reliably and automatically predict high fidelity binaries of hiCM α-actinin-2 stains that are resistant to artifacts related to differences in pixel intensity and can accurately separate individual Z-Bodies (<xref rid="fig1" ref-type="fig">Figure 1C</xref>, arrowhead). Recognizing that other labs would benefit from the flexibility and accuracy of our U-Net, we also wrote software encoding for a graphical user interface (GUI) to accompany our U-Net in the hopes of making it more user-friendly and broadly accessible. Our U-Net and GUI together are known as “yoU-Net” and are available, open source with details in the Supplement (Figure S1). We used yoU-Net to binarize images throughout the remainder of this study.</p>
</sec>
<sec id="s2b">
<title>sarcApp quantifies muscle stress fibers, sarcomeres, and myofibrils</title>
<p>The next bottleneck within our workflow was the issue of manual quantification, which in some cases can require several hours for a single CM. A meaningful quantitative output of sarcomere assembly would discriminate between sarcomeres and sarcomere precursors (muscle stress fibers (MSFs)), and produce readouts for each. Sarcomere assembly involves the formation MSFs near the cell edge that give rise to sarcomeres typically closer to the cell center (<xref rid="fig1" ref-type="fig">Figure 1D</xref>). Within MSFs are punctate α-actinin-2-containing structures termed “Z-Bodies”, which gradually elongate/concatenate over time and eventually transition into the Z-Lines of sarcomeres (<xref rid="fig1" ref-type="fig">Figure 1D</xref>). While the Z-Body to Z-Line transition exists along a continuum, it was necessary for us to define a binary transition point on the basis of α-actinin-2 structure length when a Z-Body “becomes” a Z-Line. We tasked individual lab members with independently annotating 10 images each to determine which structures constituted a Z-Body and which were Z-Lines. Together, we defined a potential Z-Line to be any α-actinin-2-containing structure with a length &gt;1.4µm, with all other structures being potential Z-Bodies. Because many organisms have variable sarcomere spacing and structure, this variable along with all others discussed (lengths, minimum number of structures, etc.) are customizable and can be easily re-defined by the user for their unique model system<sup><xref ref-type="bibr" rid="c25">25</xref></sup>.</p>
<p>Using our predefined size criterion for Z-Lines and Z-Bodies, we next desired a method capable of automated detection, measurement, and descriptive output of both structures. Towards this goal, we developed sarcApp, a software-based code that classifies each α-actinin-2 structure within a binary image as a Z-Line (<xref rid="fig2" ref-type="fig">Figure 2A-B</xref>) or Z-Body (<xref rid="fig2" ref-type="fig">Figure 2C-D</xref>) and further assigns each to either a myofibril (<xref rid="fig2" ref-type="fig">Figure 2B</xref>) or MSF (<xref rid="fig2" ref-type="fig">Figure 2D</xref>). All structures identified as potential Z-Lines are paired to other potential Z-Lines in the cell based on shape, orientation, and location, then assigned to myofibril based on orientation to other Z-Lines (<xref rid="fig2" ref-type="fig">Figure 2A-B</xref>). Z-Lines which belong to a myofibril are considered “confirmed”. This eliminates singular α-actinin-2-positive adhesions. Next, each potential Z-Body is paired to other Z-Bodies based on shape and location, then assigned to an MSF based on orientation to other nearby Z-Bodies (<xref rid="fig2" ref-type="fig">Figure 2C-D</xref>). Like with Z-Lines, only Z-Bodies which belong to an MSF are considered “confirmed”. After detecting, annotating/assigning, and measuring all potential Z-Lines and Z-Bodies within a given cell, sarcApp automatically compiles and generates three spreadsheets summarizing the data on a per-cell, per-myofibril, and per-MSF basis. Cell spreadsheets generated include myofibrils per cell (<xref rid="fig2" ref-type="fig">Figure 2E</xref>), Z-Lines per cell (<xref rid="fig2" ref-type="fig">Figure 2F</xref>), Z-Line lengths (<xref rid="fig2" ref-type="fig">Figure 2G</xref>), myofibril persistence lengths (<xref rid="fig2" ref-type="fig">Figure 2H</xref>), MSFs per cell (<xref rid="fig2" ref-type="fig">Figure 2I</xref>), Z-Bodies per cell (<xref rid="fig2" ref-type="fig">Figure 2J</xref>), Z-Body lengths (<xref rid="fig2" ref-type="fig">Figure 2K</xref>), and MSF persistence lengths (<xref rid="fig2" ref-type="fig">Figure 2L</xref>). Myofibril spreadsheets include analysis of myofibril organization within the cell by calculating the angle of each myofibril relative to the nearest cell edge (Figure M-O), as well as number of Z-Lines, average Z-Line spacing, persistence length, absolute angle of the myofibril long axis (considering the image edges as the X and Y axes), average Z-Line length, distance from the center point of the myofibril to the nearest parallel cell edge, absolute angle of the nearest cell edge, and relative angle of the myofibril to the edge (a measure of parallelism); all on a per individual myofibril basis. MSF spreadsheets include number of Z-Bodies, average Z-Body spacing, and persistence length for each MSF (Figure S3).</p>
<fig id="fig2" position="float" fig-type="figure">
<label>Figure 2:</label>
<caption><title>Quantifying sarcomere and myofibril organization using α-actinin-2 binaries.</title>
<p>A) Z-Lines and myofibrils identified by sarcApp. Each line denotes a Z-Line, and each different color represents a different myofibril. B) Quantification scheme for myofibrils and Z-Lines. Details can be found in Figure S2. C) Z-Bodies identified by sarcApp. Each red circle denotes a Z-Body. D) Quantification scheme for MSFs and Z-Bodies. Details can be found in Figure S2. E) Distribution of myofibrils per hiCM plated for 24 hours (N=188 cells; 4 biological replicates). F) Distributions of Z-Lines per hiCMs from <xref ref-type="fig" rid="fig2">Figure 2E</xref>. G) Distribution of average Z-Line length per cell from <xref ref-type="fig" rid="fig2">Figure 2F</xref> (N=104 cells). H) Distribution of average myofibril persistence lengths per cell from 2G (N=104 cells; quantification details can be found in Figure S2). I) Distribution of MSFs per cell from <xref ref-type="fig" rid="fig2">Figure 2E</xref>. J) Distribution of Z-Bodies per cell from <xref ref-type="fig" rid="fig2">Figure 2E</xref>. K) Distribution of average Z-Body length per cell from <xref ref-type="fig" rid="fig2">Figure 2E</xref>. L) Distribution of average MSF persistence lengths per cell from <xref ref-type="fig" rid="fig2">Figure 2E</xref>. M) Myofibril long axes identified by sarcApp. N) Quantification scheme for myofibril angle relative to edge. Briefly, the closest edge segment to the myofibril long axis (perpendicularly) is used as the reference angle, and the numerical output is the difference between the myofibril long axis angle and the reference edge angle. O) Distribution of myofibril orientation relative to cell edge in the same cells as <xref ref-type="fig" rid="fig2">Figure 2E</xref>. (N=188 cells, 1217 myofibrils). Note that most myofibrils are relatively parallel to the edge in hiCMs plated for 24 hours.</p></caption>
<graphic xlink:href="523681v2_fig2.tif" mimetype="image" mime-subtype="tiff"/>
</fig>
<p>sarcApp must take certain liberties to parametrize and predict MSF/myofibrillar structures within CMs. For example, Z-Lines are assigned to individual myofibrils in part by inter-Z-Line distance (which we define as being &lt;3 µm apart using centroid-to-centroid distance) but also by Z-Line long-axis angle (which we define as being within 30 degrees, accounting for cell curvature). These measurements are customizable, as some species have varied inter-Z-Line distance or varied cell curvature. We increase the computational speed of sarcApp by measuring Z-Line spacing using the distance between each Z-Line’s centroid as opposed to parallel spacing, which can produce slightly different results (Figure S2A). sarcApp defines a myofibril as a structure with a linear collection of 4 or more Z-Lines (i.e., 3 sarcomeres). Potential Z-Lines that are not assigned to a myofibril are not quantified, as they could represent adhesions, which are also α-actinin-2-positive<sup><xref ref-type="bibr" rid="c22">22</xref></sup>. Any structure with a long axis &lt;1.4 µm is defined as a potential Z-Body. Z-Bodies are assigned to MSFs with the criteria that the Z-Bodies are 1) &lt;3 µm apart and 2) have centroids that form a line within 30 degrees parallel to the edge, indicative of a typical Z-Body organization with MSF actin filaments<sup><xref ref-type="bibr" rid="c11">11</xref></sup>.</p>
</sec>
<sec id="s2c">
<title>Blebbistatin treatment reduces Z-Line assembly in hiCMs</title>
<p>To validate the functionality and accuracy of sarcApp, we paired our lab’s sarcomere assembly assay<sup><xref ref-type="bibr" rid="c11">11</xref>,<xref ref-type="bibr" rid="c22">22</xref></sup> with a well-established inhibitor of sarcomere assembly: Blebbistatin, a pan-myosin II inhibitor<sup><xref ref-type="bibr" rid="c26">26</xref></sup>. In brief, we trypsinized and re-plated human induced pluripotent stem cell-derived cardiac myocytes (hiCMs). During this process, the hiCMs lose and re-form their sarcomeres, allowing us to monitor sarcomere assembly. Because this assembly begins at the cell edge and moves rearward toward the cell center, a snapshot of a hiCM after 24 hours of spreading (once sarcomeres have begun to form and myofibrils begin to align) contains information on early, mid, and late sarcomere formation at the leading edge, behind the edge in the lamella, and toward the center of the cell, respectively<sup><xref ref-type="bibr" rid="c11">11</xref></sup>.</p>
<p>Blebbistatin inhibits myosin II-based contractility and Blebbistatin-treated hiCMs exhibit reduced Z-Line assembly in our assay (<xref rid="fig3" ref-type="fig">Figure 3A</xref>, S3-5), consistent with other reports<sup><xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c27">27</xref></sup>. Platinum-replica electron microscopy (EM) shows that Blebbistatin-treated hiCMs do not appear to have well-defined Z-Lines and have a disordered overall orientation of MSFs/myofibrils compared to controls (<xref rid="fig3" ref-type="fig">Figure 3B</xref>). As measured by sarcApp, hiCMs treated with Blebbistatin exhibited fewer myofibrils per cell (<xref rid="fig3" ref-type="fig">Figure 3C</xref>), fewer Z-Lines per cell (<xref rid="fig3" ref-type="fig">Figure 3D</xref>), decreased myofibril persistence length (<xref rid="fig3" ref-type="fig">Figure 3E</xref>), decreased Z-Line length (<xref rid="fig3" ref-type="fig">Figure 3F</xref>), and decreased overall size of all α-actinin-2-positive puncta/structures (<xref rid="fig3" ref-type="fig">Figure 3G</xref>). sarcApp also reported that, while the myofibrils of control hiCMs are oriented ∼parallel to the nearest cell edge (<xref rid="fig3" ref-type="fig">Figure 3H</xref>; similar to <xref rid="fig2" ref-type="fig">Figure 2M-O</xref>), the myofibrils of Blebbistatin-treated hiCMs are more randomly oriented (<xref rid="fig3" ref-type="fig">Figure 3H</xref>), suggesting a disorganized assembly process. Altogether, these data generated by sarcApp are consistent with a wide range of previous studies showing that Blebbistatin inhibits sarcomere assembly<sup><xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c27">27</xref></sup>.</p>
<fig id="fig3" position="float" fig-type="figure">
<label>Figure 3:</label>
<caption><title>Blebbistatin treatment ablates Z-Line formation in hiCMs</title>
<p>A) Representative images and insets of α-actinin-2 and F-actin in hiCMs treated with DMSO, 50 µM Blebbistatin, or 100 µM Blebbistatin. B) Representative platinum replica EM image of a control hiCM and a 50 µM Blebbistatin-treated hiCM. Arrows indicate an elongated Z-Line in the DMSO-treated hiCM and a Z-Body in the Blebbistatin-treated hiCM. C) Myofibrils per cell in hiCMs (N=118 DMSO cells, 108 50 µM Blebbistatin cells, and 93 100 µM Blebbistatin cells; 4 biological replicates). D) Z-Lines per cell in hiCMs from <xref ref-type="fig" rid="fig3">Figure 3C</xref>. E) Average myofibril persistence length per cell in hiCMs from <xref ref-type="fig" rid="fig3">Figure 3C</xref> (N=104 DMSO control cells, 45 50 µM Blebbistatin cells, and 45 100 µM Blebbistatin cells). F) Average Z-Line length per cell from <xref ref-type="fig" rid="fig3">Figure 3E</xref>. G) Average size of all α-actinin-2-positive puncta in hiCMs from <xref ref-type="fig" rid="fig3">Figure 3C</xref>. H) Myofibril orientation relative to the cell edge segment closest to the myofibril center, perpendicularly. N=4 biological replicates, 1217 DMSO control myofibrils, 385 50 µM Blebbistatin myofibrils, 220 100 µM Blebbistatin myofibrils.</p></caption>
<graphic xlink:href="523681v2_fig3.tif" mimetype="image" mime-subtype="tiff"/>
</fig>
</sec>
<sec id="s2d">
<title>sarcApp quantifies titin in hiCMs</title>
<p>Our next goal was to introduce additional quantitative modes of analysis to sarcApp that could be used concurrently with α-actinin-2-based quantification of Z-line assembly. Our focus was on other core components of the sarcomere which, like α-actinin-2, are thought to play an active role in assembly and can be stained and visualized using commercially available antibodies. Our group has shown previously that the protein titin can be used as a spatiotemporal indicator of myofibril maturation status in hiCMs<sup><xref ref-type="bibr" rid="c22">22</xref></sup>. Antibodies to the titin I-band localize to both MSFs and sarcomeres in hiCMs (<xref rid="fig4" ref-type="fig">Figure 4A</xref>). Titin forms ringlike structures around the Z-Bodies of MSFs that are closer to the apparent sarcomere transition point (<xref rid="fig4" ref-type="fig">Figure 4A</xref>). <xref rid="fig4" ref-type="fig">Figure 4B</xref> shows our current model for how titin is oriented around Z-bodies during earlier stages of sarcomere assembly, forming rings. Likely, due to both antibody localization to the titin I-Band region and the ability of titin to bind both actin and the α-actinin-2 N-terminus, titin is oriented with the N-terminus toward the center of the Z-Body and the C-terminus radially oriented outward (<xref rid="fig4" ref-type="fig">Figure 4C</xref>)<sup><xref ref-type="bibr" rid="c28">28</xref>,<xref ref-type="bibr" rid="c29">29</xref></sup>. This is a topic of current research in our lab.</p>
<fig id="fig4" position="float" fig-type="figure">
<label>Figure 4:</label>
<caption><title>sarcApp uses titin binaries to identify myofibrils and precursor ring structures</title>
<p>A) Representative image of titin and α-actinin-2 in a control hiCM. Arrow: an α-actinin-2-positive Z-Body with titin localized in a ring around it. B) Model of titin localization during sarcomere formation. C) Titin doublets identified by sarcApp. Each line denotes a doublet with titin localized, and each color is a myofibril. E) Quantification scheme for myofibrils and titin doublets. Details can be found in Figure S2 and the Supplemental Methods. E) Titin precursor rings identified by sarcApp (red). F) Quantification scheme for titin precursor rings. Details can be found in Figure S2 and the Supplemental Methods.</p></caption>
<graphic xlink:href="523681v2_fig4.tif" mimetype="image" mime-subtype="tiff"/>
</fig>
<p>yoU-Net-generated binaries of titin stains after training with annotated ground-truth binaries accurately recapitulate titin rings (MSFs) and titin doublets (sarcomeres) in hiCMs (<xref rid="fig4" ref-type="fig">Figure 4A</xref>, S1D). We parametrized each titin-based structure into a quantitative scheme and have outfitted sarcApp with a modality for measuring titin. To determine the minimum length required for a potential titin doublet to be considered as part of a sarcomere, we co-stained hiCMs with both titin and α-actinin-2, and found that on average, titin doublets are 0.3 µm longer than the Z-Lines with which they associate. Thus, we set a lower boundary of titin doublet length at 1.7 µm (Z-lines: 1.4 µm). Closer inspection of titin images revealed that, in several cases, the titin I-Band signal “wraps” around the radial tip of Z-Lines, reminiscent of titin rings around Z-Bodies.</p>
<p>sarcApp accurately distinguishes between titin doublets (<xref rid="fig4" ref-type="fig">Figure 4C</xref>) and titin rings (<xref rid="fig4" ref-type="fig">Figure 4E</xref>) and produces readouts for number of myofibrils doublets per cell, myofibril persistence length, doublet length, number of rings, ring diameter, ring aspect ratio, distance of doublets from the edge, and distance of rings from the edge. (<xref rid="fig4" ref-type="fig">Figures 4C-F</xref>). sarcApp is additionally equipped to incorporate a cell edge co-stain (e.g., actin or non-muscle myosin IIA or IIB) alongside titin to measure myofibril orientation and myofibril distance from edge since titin stains do not label the edge (similar to <xref rid="fig2" ref-type="fig">Figure 2 M-O</xref>).</p>
</sec>
<sec id="s2e">
<title>Blebbistatin affects myofibril orientation and doublet length, but not the periodicity of titin structures</title>
<p>We wished to use our sarcomere formation assay in the presence of Blebbistatin to quantify titin structure alignment using sarcApp. As before, hiCMs were re-plated in the presence of Blebbistatin for 24 hours, then fixed and stained for titin as well as phalloidin (F-actin) to aid in the visualization of MSFs. While F-actin in Blebbistatin-treated hiCMs is mostly punctate, titin rings and doublets can still be observed by eye around punctate actin with some degree of organization and periodicity (<xref rid="fig5" ref-type="fig">Figure 5A</xref>). sarcApp-dependent quantification of myofibrils using titin revealed no significant differences in the numbers of myofibrils in hiCMs replated in moderate (50 µM) or high (100 µM) Blebbistatin (<xref rid="fig5" ref-type="fig">Figure 5B</xref>), but found fewer and shorter titin doublets in high Blebbistatin (<xref rid="fig5" ref-type="fig">Figure 5C+D</xref>). sarcApp-based quantification of titin precursor rings within MSFs revealed no difference in absolute number of rings per cell with Blebbistatin treatment (<xref rid="fig5" ref-type="fig">Figure 5E</xref>) but rings that were detected were significantly more rounded as indicated by a decreased aspect ratio (<xref rid="fig5" ref-type="fig">Figure 5F</xref>). sarcApp also detected that myofibrils as quantified by titin were mostly parallel to the edge in control cells but became more randomly aligned relative to the edge in Blebbistatin (<xref rid="fig5" ref-type="fig">Figure 5G</xref>), consistent with sarcApp-dependent quantification of α-actinin-2. In summary, Blebbistatin treatment results in disarrayed and shorter titin doublets that are still capable of aligning into periodic myofibrillar structures.</p>
<fig id="fig5" position="float" fig-type="figure">
<label>Figure 5:</label>
<caption><title>Blebbistatin affects myofibril orientation and the morphology of titin structures</title>
<p>A)Representative images of titin and F-actin in hiCMs treated with DMSO, 50 µM Blebbistatin, and 100 µM Blebbistatin. B)Myofibrils per cell in hiCMs (N=3 biological replicates, 107 DMSO cells, 95 50 µM Blebbistatin cells, and 84 100 µM Blebbistatin cells). C) Doublets per cell in hiCMs from <xref ref-type="fig" rid="fig5">Figure 5B</xref>. D) Doublet length per cell in hiCMs from <xref ref-type="fig" rid="fig5">Figure 5B</xref>: (N= 58 DMSO cells, 32 50 µM Blebbistatin cells, and 21 100 µM Blebbistatin cells. E) Rings per cell in hiCMs from <xref ref-type="fig" rid="fig5">Figure 5B</xref>. F) Average ring aspect ratio per cell in hiCMs from <xref ref-type="fig" rid="fig5">Figure 5B</xref> (N= 105 DMSO control cells, 95 50 µM Blebbistatin cells, and 78 100 µM Blebbistatin cells. G) Myofibril orientation (N=292 DMSO control myofibrils, 151 50 µM Blebbistatin myofibrils, and 69 100 µM Blebbistatin myofibrils</p></caption>
<graphic xlink:href="523681v2_fig5.tif" mimetype="image" mime-subtype="tiff"/>
</fig>
</sec>
<sec id="s2f">
<title>sarcApp quantifies myomesin, a component of the M-Line, in hiCMs</title>
<p>Having equipped sarcApp to quantify two structures associated with Z-Line assembly in hiCMs, our next goal was to introduce a modality for M-Line assembly. The M-Line lies at the midline of the stack of myosin II filaments in the sarcomere and demarcates the midpoint between two Z-Lines (<xref rid="fig6" ref-type="fig">Figure 6A</xref>). M-Lines in hiCMs can be stained using antibodies to the protein myomesin. Myomesin stains closely resemble α-actinin-2 Z-Line stains although myomesin is apparently absent from Z-Bodies (<xref rid="fig6" ref-type="fig">Figure 6B</xref>). Like α-actinin-2 and titin stains, yoU-Net can be used to predict myomesin binaries with high accuracy (<xref rid="fig6" ref-type="fig">Figure 6B</xref>, Figure S1E)</p>
<fig id="fig6" position="float" fig-type="figure">
<label>Figure 6:</label>
<caption><title>sarcApp uses myomesin binaries to identify myofibrils and M-Lines in hiCMs</title>
<p>A) Schematic showing myomesin localization at the M-Line. B) Representative image of myomesin and F-actin in a hiCM. The myomesin binary was predicted using yoU-net as described in the Supplemental Methods. C) Myofibrils and M-Lines identified by sarcApp. Each line denotes an M-Line, and each color represents a myofibril. D) Quantification scheme for myofibrils and M-Lines. Details can be found in Figure S2. E) Representative images of myomesin and F-actin in hiCMs treated with DMSO, 50 µM Blebbistatin, and 100 µM Blebbistatin. F) Myofibrils per cell in hiCMs (N=3 biological replicates; 112 DMSO cells, 90 50 µM Blebbistatin cells, and 89 100 µM Blebbistatin cells). G) M-Lines per cell in hiCMs from <xref ref-type="fig" rid="fig6">Figure 6F</xref>. H) Average M-Line length per cell in hiCMs from <xref ref-type="fig" rid="fig6">Figure 6F</xref> (N=97 DMSO control cells, 16 50 µM Blebbistatin cells, and 13 100 µM Blebbistatin cells). I) Myofibril orientation (N=752 DMSO control myofibrils, 49 50 µM Blebbistatin myofibrils, and 60 100 µM Blebbistatin myofibrils).</p></caption>
<graphic xlink:href="523681v2_fig6.tif" mimetype="image" mime-subtype="tiff"/>
</fig>
<p>We equipped sarcApp to detect and measure myomesin M-Lines from yoU-Net-generated binaries (<xref rid="fig6" ref-type="fig">Figure 6C</xref>) based on a scheme that resembles the one designed for Z-Lines (<xref rid="fig6" ref-type="fig">Figure 6D</xref>). For myomesin, the M-Line must &gt;1.4 µm and adjacent to at least two other M-lines to be considered part of a myofibril (i.e., at least three M-lines per myofibril, as it could be inferred that such a myofibril would have at least four Z-lines). sarcApp-defined readouts using myomesin include myofibrils per cell, M-Lines per cell, myofibril persistence length, and M-Line lengths. As with titin, sarcApp is additionally equipped to incorporate a cell edge co-stain (e.g., actin or non-muscle myosin IIA or IIB) alongside myomesin to measure myofibril orientation and myofibril distance from cell edge since myomesin antibodies do not stain the edge (see <xref rid="fig2" ref-type="fig">Figure 2M-O</xref>). To validate sarcApp as a tool to measure M-Lines, we stained myomesin in hiCMs exposed to Blebbistatin. Representative images in <xref rid="fig6" ref-type="fig">Figure 6E</xref> show that control hiCMs exhibit elongated M-Lines while Blebbistatin-treated M-Lines appear truncated and punctate-like. sarcApp quantification matched observations made by eye, detecting fewer myofibrils per cell (<xref rid="fig6" ref-type="fig">Figure 6F</xref>) and M-Lines per cell (<xref rid="fig6" ref-type="fig">Figure 6G</xref>) in Blebbistatin-treated hiCMs, as well as shortened M-Line length (<xref rid="fig6" ref-type="fig">Figure 6H</xref>). Co-staining with actin enabled sarcApp quantification of myofibril orientation, revealing that myofibril M-Lines are disorganized relative to the cell edge (<xref rid="fig6" ref-type="fig">Figure 6I</xref>, see <xref rid="fig2" ref-type="fig">Figure 2 M-O</xref>) akin to the Blebbistatin-induced phenotype seen in Z-Lines (<xref rid="fig3" ref-type="fig">Figure 3H</xref>).</p>
</sec>
<sec id="s2g">
<title>Neither α nor β cardiac myosin II alone is required for Z-Line assembly</title>
<p>The sarcomere A-band is the region of thick and thin filament overlap (<xref rid="fig7" ref-type="fig">Figure 7A</xref>). The thick filament contains stacks of muscle myosins that produce the forces of contraction<sup><xref ref-type="bibr" rid="c18">18</xref></sup>. Previous studies probing the role of muscle myosins in sarcomere assembly have produced conflicting results<sup><xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c30">30</xref></sup>. With sarcApp capable of measuring both M-Lines and Z-Lines, we were well-poised to ask if and/or how muscle myosins influence the specific assembly of both M-Lines and Z-Lines in hiCMs. Cardiac sarcomeres contain both α myosin II (<italic>MYH6</italic>) and β myosin II (<italic>MYH7</italic>)<sup><xref ref-type="bibr" rid="c18">18</xref>,<xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c31">31</xref></sup>; we exposed hiCMs to one of two unique sequences targeting either <italic>MYH-</italic> or <italic>MYH7</italic> alongside a non-targeting siRNA (control) then re-plated, fixed, and stained hiCMs with antibodies to either α-actinin-2 (Z-Lines), titin (Z-Lines), or myomesin (M-Lines) (<xref rid="fig7" ref-type="fig">Figure 7B-C, H-I</xref>).</p>
<fig id="fig7" position="float" fig-type="figure">
<label>Figure 7:</label>
<caption><title>Knockdown of α or β cardiac myosin II reduces but does not eliminate sarcomeres</title>
<p>A) Schematic showing cardiac myosin localization in a sarcomere. B) Representative western blot showing α cardiac myosin (MYH6) knockdown in hiCMs. C) Representative images of α-actinin-2, titin, and myomesin in siControl hiCMs and α cardiac myosin (MYH6) knockdown hiCMs. D) Number of Z-Lines per cell in hiCMs in two independent groups of siControl (scramble)-treated hiCMs and two separate MYH6 siRNA-treated hiCMs (sequences 1 and 2) (N=3 biological replicates, 81 siCon cells and 78 siMYH6 (1) cells, and 88 siCon cells and 63 siMYH6 (2) cells). E) Average Z-Line length per cell in hiCMs from <xref ref-type="fig" rid="fig7">Figure 7D</xref> (N=68 siCon cells and 75 siMYH6 (1) cells, and 83 siCon cells and 62 siMYH6 (2) cells). F) Average doublet length per cell in hiCMs in siCon (scramble)-treated hiCMs and two MYH6 siRNA sequences (1 and 2) (71 siCon cells, 49 siMYH6 (1) cells, and 50 siMYH6 (2) cells. G) Average M-Line length per hiCMs (85 siCon cells, 78 siMYH6 (1) cells, and 64 siMYH6 (2) cells). H) Representative western blot showing β cardiac myosin (MYH7) knockdown in hiCMs. I) Representative images of α-actinin-2, titin, and myomesin in β cardiac myosin (MYH7) knockdown hiCMs. J) Number of Z-Lines per cell in hiCMs in two independent groups of siCon (scramble)-treated hiCMs and two separate MYH7 siRNA-treated hiCMs (pools 1 and 2) (N= 86 siCon cells and 66 siMYH7 (1) cells, and 97 siCon cells and 62 siMYH7 (2) cells). L) Average Z-Line length per cell in hiCMs from <xref ref-type="fig" rid="fig7">Figure 7J</xref> (81 siCon cells and 63 siMYH7 (1) cells, and 81 siCon cells and 59 siMYH7 (2) cells). M) Number of rings per hiCM (N=3 biological replicates, 94 siCon cells, 72 siMYH7 (1) cells, and 66 siMYH7 (2) cells). N) Average M-Line length per cell in hiCMs (N= 87 siCon cells, 65 siMYH7 (1) cells, and 62 siMYH7 (2) cells).</p></caption>
<graphic xlink:href="523681v2_fig7.tif" mimetype="image" mime-subtype="tiff"/>
</fig>
<p>After 24 hours, we observed by eye that both the Z-Lines and M-Lines of hiCMs depleted of α myosin II (<italic>MYH6</italic>) appeared to largely resemble those of controls (<xref rid="fig7" ref-type="fig">Figure 7C</xref>). sarcApp-dependent quantification of each stain determined that α myosin II knockdown did not alter length of α-actinin-2 Z-lines, titin doublets, or M-lines, though one of the two siRNAs did result in increased numbers of Z-lines (<xref rid="fig7" ref-type="fig">Figure 7D-G</xref>). As with α myosin II, stains of β myosin II-depleted hiCMs (<italic>MYH7</italic>) also appeared to resemble those of controls by eye (<xref rid="fig7" ref-type="fig">Figure 7H-I</xref>) – however, sarcApp-based quantification revealed that hiCMs depleted of β myosin II assembled more abundant, but shorter Z-lines than controls (<xref rid="fig7" ref-type="fig">Figure 7J-K</xref>). Further, while sarcApp detected no change in titin doublet lengths (<xref rid="fig7" ref-type="fig">Fig. 7L</xref>), M-lines were detected as being significantly shorter than controls following β myosin II knockdown (<xref rid="fig7" ref-type="fig">Figure 7M</xref>). Thus, we conclude that neither α nor β myosin II is strictly required for Z- or M-line assembly in hiCMs but that β myosin II likely influences Z-line concatenation. This is not entirely surprising, as Geach and colleagues showed similar results in Xenopus embryos<sup><xref ref-type="bibr" rid="c32">32</xref></sup>.</p>
</sec>
<sec id="s2h">
<title>Myomesin is required for A-Band but not Z-Line assembly</title>
<p>Given our data suggesting that Z-Line and M-Line assembly could be uncoupled (<xref rid="fig7" ref-type="fig">Figure 7</xref>), we next asked if directly disrupting M-Line assembly impacts Z-Line assembly. We disrupted M-Line assembly by exposing hiCMs to one of two types of myomesin<italic>-</italic>targeting siRNAs (gene name: <italic>MYOM</italic>) alongside a non-targeting siRNA (control) then re-plated, fixed, and stained hiCMs with antibodies to either α-actinin-2 or titin (<xref rid="fig8" ref-type="fig">Figure 8A-B</xref>).</p>
<fig id="fig8" position="float" fig-type="figure">
<label>Figure 8:</label>
<caption><title>Myomesin knockdown alters titin and cardiac myosin II localization, but not α-actinin-2</title>
<p>A) Representative images of sarcomeric proteins α-actinin-2 and titin in myomesin (MYOM) knockdown hiCMs. B) Representative western blot and quantification showing MYOM knockdown in hiCMs, N=3 biological replicates. C) Number of Z-Lines per cell in siCon (scramble)-treated hiCMs and two separate MYOM siRNA-treated hiCMs (sequences 1 and 2). N=3 biological replicates, 132 siCon cells, 105 siMYOM (1) cells, and 103 siMYOM (2) cells. D) Average Z-Line length per cell in hiCMs from <xref ref-type="fig" rid="fig8">Figure 8C</xref>: N=4 biological replicates, 117 siCon cells, 104 siMYOM (1) cells, and 92 siMYOM (2) cells (only cells with myofibrils were quantified for <xref ref-type="fig" rid="fig8">Figure 8D</xref>). E) Titin doublets per cell in hiCMs using experimental treatments in <xref ref-type="fig" rid="fig8">Figure 8C</xref>. N=3 biological replicates, 117 siCon cells, 90 siMYOM (1) cells, and 100 siMYOM (2) cells. F) Average doublet length per cell in hiCMs from <xref ref-type="fig" rid="fig8">Figure 8E</xref>: N=3 biological replicates, 74 siCon cells, 68 siMYOM (1) cells, and 70 siMYOM (2) cells. (Only cells with myofibrils were quantified for <xref ref-type="fig" rid="fig8">Figure 8F</xref>). G) Rings per cell in hiCMs from <xref ref-type="fig" rid="fig8">Figure 8E</xref>. H) Representative images of MYH7 and F-actin in siCon, siMYOM (1), and siMYOM (2) hiCMs. I) Length of β cardiac myosin stacks in siCon, siMYOM (1), and siMYOM (2) hiCMs</p></caption>
<graphic xlink:href="523681v2_fig8.tif" mimetype="image" mime-subtype="tiff"/>
</fig>
<p>After 24 hours, we observed by eye that both the α-actinin-2 and titin stains of myomesin-depleted hiCMs resemble those of controls (<xref rid="fig8" ref-type="fig">Figure 8A</xref>). sarcApp-dependent quantification of α-actinin-2 detected no difference in Z-Lines per cell in myomesin-depleted hiCMs compared to controls (<xref rid="fig8" ref-type="fig">Figure 8C</xref>), suggesting Z-Lines assemble independently of myomesin – however, sarcApp did detect a decrease in Z-Line length for one myomesin-targeting siRNA (<xref rid="fig8" ref-type="fig">Figure 8D</xref>, S14). Meanwhile sarcApp-dependent quantification of titin detected no significant change in the number of titin doublets or doublet lengths (<xref rid="fig8" ref-type="fig">Figure 8E-F</xref>). Interestingly, upon myomesin knockdown, the number of titin precursor rings significantly increase, warranting future investigation (<xref rid="fig8" ref-type="fig">Figure 8G</xref>).</p>
<p>We next wanted to assess the influence of myomesin knockdown on M-Line/A-band assembly. Because we were unable to visualize M-Lines in myomesin-depleted hiCMs (data not shown), we stained the A-band using antibodies to the β myosin II motor domain, observing that A-Bands in myomesin-depleted hiCMs appear to be less organized (<xref rid="fig8" ref-type="fig">Figure 8H</xref>) with shorter β myosin II stacks (<xref rid="fig8" ref-type="fig">Figure 8I</xref>). These data suggest that myomesin regulates A-band/M-Line assembly and myofibril maturation, but not Z-Line assembly. Altogether, these data in tandem with data from <xref rid="fig7" ref-type="fig">Figure 7</xref> suggest that the assembly of the sarcomere A-Band/M-Line and Z-Line are decoupled and/or are regulated by mechanisms that are at least partially exclusive to one another.</p>
</sec>
</sec>
<sec id="s3">
<title>Discussion</title>
<p>Here, we introduce sarcApp and yoU-Net, two open-source software packages to binarize microscope images and quantify sarcomere components. We first demonstrate deep learning-based image binarization of three distinct sarcomeric proteins – α-actinin-2, titin, and myomesin – using models trained by a framework we developed called yoU-Net, based off the ubiquitously-used U-Net<sup><xref ref-type="bibr" rid="c16">16</xref></sup>. Each trained model is provided and will function as demonstrated herein based on the use of the antibodies as stated in the Methods. yoU-Net is designed and equipped to be functionally versatile which we hope other researchers can harness to generate additional trained models of proteins not studied here.</p>
<p>The high fidelity of yoU-Net-generated binaries allowed us to quantify several geometric features of the sarcomere in real space as opposed to frequency space. Automatic quantification of sarcomere/myofibril structure with sarcApp offers multiple advantages over traditional, manual quantification including speed, objectivity, and consistency. sarcApp is capable of quantifying 33 measurements per cell and 24 measurements per myofibril at a rate that is orders of magnitude faster than manual quantification, which besides being slower also requires complete attention from the user. sarcApp also removes the potential for user bias/error from quantification and can detect quantitative differences in cells that are not readily apparent by eye (<xref rid="fig7" ref-type="fig">Figures 7</xref>, <xref rid="fig8" ref-type="fig">8</xref>) – we expect this will drive discovery of novel mechanisms in cells that would have remained otherwise undetectable. sarcApp as an automatic tool must operate within a rigid quantitative framework and functions optimally when structures-of-interest have shapes that are consistent and definable by simple geometric measurements, however, to increase usability we have enabled users to self-define parameters using an easy-to-navigate GUI.</p>
<p>Other than basic structural readouts for sarcomere/myofibril precursors (Z-Bodies, titin rings, etc.) and sarcomeres/myofibrils (Z-Lines, titin doublets, M-Lines, etc.), sarcApp is also equipped to proxy sarcomere/myofibril organization. Myofibril organization correlates with function <italic>in vivo</italic> and poor organization (disarray) is often cited as causative for myocardial disease<sup><xref ref-type="bibr" rid="c4">4</xref>-<xref ref-type="bibr" rid="c6">6</xref></sup>. sarcApp reports organization using a novel metric reported here known as myofibril orientation (<xref rid="fig2" ref-type="fig">Figure 2 M-O</xref>). Orientation is reported relative to the cell membrane by comparing the angle of the myofibril long axis to the nearest cell edge; such a metric was chosen because we routinely observe hiCMs assemble myofibrils with a long axis that is ∼parallel (i.e., 0 degrees) to the cell edge <italic>in vitro</italic> and also because tissue sections of myocytes <italic>in vivo</italic> show myofibrillar arrays parallel to the long axis of the membrane<sup><xref ref-type="bibr" rid="c11">11</xref>,<xref ref-type="bibr" rid="c22">22</xref></sup>. We generated this measure of orientation to be more transparent and informative than what is currently available in the hopes of reaching a broader audience<sup><xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c21">21</xref></sup>.</p>
<p>Previous work from our lab and others have included as many cells as reasonably possible in analyses<sup><xref ref-type="bibr" rid="c11">11</xref>,<xref ref-type="bibr" rid="c22">22</xref></sup>, requiring technical expertise in (at minimum) plating cells, sample preparation, microscopy, and image analysis. While microscopy has historically been the rate-limiting step in such analyses, modern microscopes can now rapidly image hundreds to thousands of cells at high resolution. The true bottleneck towards understanding sarcomere assembly now lies for our lab in image analysis/quantification. We summed the quantified cell totals of two recent landmark papers in the field and found the reported totals to be 385<sup><xref ref-type="bibr" rid="c19">19</xref></sup> and 785<sup><xref ref-type="bibr" rid="c11">11</xref></sup>. sarcApp has allowed us to present in this paper multi-dimensional outputs from 3,452 single cells with outputs that we feel represent the data clearly at both basic glance and upon close inspection. Much like yoU-net, we have encoded sarcApp into a graphical user interface (GUI) equipped with a user manual that details the steps one needs to use sarcApp, irrespective of coding expertise.</p>
<p>We validated sarcApp using the pan-myosin II inhibitor Blebbistatin<sup><xref ref-type="bibr" rid="c26">26</xref></sup>. Blebbistatin exposure results in predictable, dose-dependent inhibition of sarcomere assembly<sup><xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c27">27</xref></sup>. sarcApp-generated measurements of α-actinin-2 in the presence or absence of Blebbistatin were consistent with observations made by previous groups<sup><xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c27">27</xref></sup> (<xref rid="fig3" ref-type="fig">Figure 3</xref>). Meanwhile, several other outputs related to titin and myomesin first reported here by sarcApp suggest that periodic, pseudomyofibrillar structures characteristic to both the Z-Line and M-Line can at least partially self-assemble in the absence of myosin II motor-based contractility (<xref rid="fig5" ref-type="fig">Figures 5</xref>-<xref ref-type="fig" rid="fig6">6</xref>). These data indicate other mechanisms beyond myosin II-based motor-based contractility may independently contribute to the assembly and ordered arrangement of myofibrils in hiCMs. Specifically with titin, we also observe the phenomenon that Blebbistatin disrupts the higher-order organization of titin-based structures with respect to the cell membrane but does not affect titin periodicity within the structure. These data suggest that unique mechanistic tiers of organization exist within the cell at the sub- and super-structural level that become uncoupled upon relieving the cell of myosin II-dependent contractility.</p>
<p>We find that neither α nor β cardiac myosin II knockdown had any measurable impact on the assembly of Z-Lines (<xref rid="fig7" ref-type="fig">Figure 7</xref>). Attempts to determine if Z-Lines assemble in the simultaneous absence of both myosins were unsuccessful in our hands. Such an experiment, in addition to the Blebbistatin experiments reported herein and elsewhere<sup><xref ref-type="bibr" rid="c19">19</xref>,<xref ref-type="bibr" rid="c27">27</xref></sup>, could resolve whether myosin II-dependent contractility alone is sufficient to explain the role of myosin II in Z-Line assembly. Meanwhile, β cardiac myosin II knockdown did result in hiCMs with shorter M-Lines, indicating that either contractility or scaffolding by β, but not α myosin II, is required for M-Line maturation (<xref rid="fig7" ref-type="fig">Figure 7G, M</xref>). Knockdown of myomesin similarly did not impact Z-Line assembly but did result in shortening of the thick filament A-band (<xref rid="fig8" ref-type="fig">Figure 8H-I</xref>). Our ability to perturb the M-Line/A-band but not the Z-Line suggests that Z- and M-Line assembly are regulated at least in part by distinct mechanistic axes, as previously hypothesized<sup><xref ref-type="bibr" rid="c33">33</xref>-<xref ref-type="bibr" rid="c35">35</xref></sup>. It should be noted here that solely using Z-Lines as a proxy for sarcomere status would have led us to errantly conclude that β cardiac myosin II and myomesin both play no or at least minor roles in sarcomere assembly. Future studies examining assembly should therefore consider Z-Line assembly to be distinct from, or at least not wholly representative of, assembly of other sarcomere components.</p>
<p>Here we used exclusively fixed-cell immunofluorescence for high-throughput quantification of sarcomere assembly. We eliminated as much bias as possible by selecting cells for quantification solely using a nuclear DAPI stain and are therefore confident that data presented herein fully captures the true biological variability present within whole hiCM cultures. A future goal of the lab is to further adapt sarcApp for live cell quantification. Tracking specific structures across time would be a necessity for a “live” version of sarcApp, a hurdle the lab is currently working to overcome. We also designed this initial version of sarcApp to quantify in 2D rather than in 3D space. As we have published previously, sarcomeres/myofibrils assemble on the dorsal (top) surface of an hiCM plated on glass for up to 48 hours after plating, and we consider them within that timeframe to exist largely within a two-dimensional plane<sup><xref ref-type="bibr" rid="c11">11</xref>,<xref ref-type="bibr" rid="c22">22</xref></sup>. Breaking into three dimensions will require more complex mathematics and produce more convoluted outputs, which could limit adoption by other researchers, but is a long-term goal of the lab nonetheless.</p>
<p>In summary, sarcApp and yoU-Net provide frameworks for automatic binarization, segmentation, and quantitative analysis of sarcomeric proteins in striated muscle cells. We anticipate these software packages will be useful for lower-throughput single-cell analyses of sarcomere structure and/or assembly as well as for high-throughput screens to analyze several hundreds to thousands of cells at a time. In this way we hope to facilitate the entry of more researchers into the field of cardiac cell biology and we continue to introduce novel quantification metrics and additional sarcomere proteins in the future through continued maintenance and management of sarcApp as an open-source community project.</p>
</sec>
<sec id="s4">
<title>Methods</title>
<p><bold>SarcApp and yoU-net can be found at the following address: <ext-link ext-link-type="uri" xlink:href="https://github.com/abbieneininger/sarcApp">https://github.com/abbieneininger/sarcApp</ext-link></bold></p>
<table-wrap id="tbl1" orientation="portrait" position="float">
<label>Table 1:</label>
<caption><title>Key Resources</title></caption>
<graphic xlink:href="523681v2_tbl1.tif" mimetype="image" mime-subtype="tiff"/>
<graphic xlink:href="523681v2_tbl1a.tif" mimetype="image" mime-subtype="tiff"/>
</table-wrap>
<sec id="s4a">
<title>The Supplemental User Manual and sarcApp Software</title>
<p>Details on sarcApp and yoU-Net download and usage documentation are in the supplement and at <ext-link ext-link-type="uri" xlink:href="http://www.github.com/abbieneininger/sarcApp">www.github.com/abbieneininger/sarcApp</ext-link></p>
</sec>
<sec id="s4b">
<title>Cell Culture and Authentication</title>
<p>Human iPSC-derived cardiac myocytes (CMM-100-012-000.5, Cellular Dynamics, Madison, WI) were cultured as per manufacturer’s instructions in proprietary manufacturer-provided cardiac myocyte maintenance medium in polystyrene 96-well cell culture plates. Cells were maintained at 37°C and 5% CO<sub>2</sub>. For re-plating hiCMs onto glass substrates, cells were washed 2 times with 100 µL 1x PBS with no Ca<sup>2+</sup>/Mg<sup>2+</sup> (PBS*, 70011044, Gibco, Grand Island, NY). PBS* was completely removed from hiCMs and 40 µL 0.1% Trypsin-EDTA with no phenol red (15400054, Gibco, Grand Island, NY) was added to hiCMs and incubated at 37°C for 2 minutes. Following incubation, the cells were washed 3 times with trypsin, the plate rotated 180 degrees, and washed another 3 times. Trypsinization was then quenched by adding 120 µL of culture media and total cell mixture was pipetted into a 1.5 mL Eppendorf tube. Cells were centrifuged at 200xg for 3 minutes, and the supernatant was aspirated. The cell pellet was re-suspended in 200 µL of culture media and plated on 35 mm dishes with a 10 mm glass bottom (D35-10-1.5-N; CellVis, Sunnydale, CA) pre-coated with 10 µg/mL fibronectin (354008, Corning) for 1 hr at 37°C.</p>
</sec>
<sec id="s4c">
<title>Antibodies</title>
<p>Alexa Fluor 488-phalloidin (A12379), Alexa Fluor 568-phalloidin (A12380), and Alexa Fluor 647-phalloidin (A22287) were purchased from Invitrogen. Alexa Fluor 488-goat anti-mouse (A11029), Alexa Fluor 488-goat anti-rabbit (A11034), Alexa Fluor 568-goat-anti-rabbit (A11011), Alexa Fluor 568-goat anti-mouse (A11004), Alexa Fluor 647-goat-anti-mouse (A32728), and Alexa Fluor 647-goat-anti-rabbit (A32733) (1:100) antibodies were purchased from Life Technologies (Grand Island, NY).</p>
<p>The Titin (9D10) antibody, Myomesin (MYOM) antibody, and MYH7 (A4.591) (all 1:2) antibody were purchased from the Developmental Studies Hybridoma Bank (University of Iowa). Mouse anti-α-actinin-2 (1:200, A7811) was purchased from Sigma Aldrich. Rabbit anti-MYH6 and MYH7 for western blotting (1:500) were purchased from ProteinTech (22281-1-AP, 22280-1-AP).</p>
</sec>
<sec id="s4d">
<title>Chemicals</title>
<p>Blebbistatin was purchased from Sigma (B0560) and reconstituted to 10 mM in DMSO (Sigma 276855).</p>
</sec>
<sec id="s4e">
<title>Fixation and Immunostaining</title>
<p>Cells were fixed with 4% paraformaldehyde (PFA) in PBS at room temperature for 20 min and then permeabilized for 5 min with 1% Triton X-100/4% PFA in PBS. For actin visualization, phalloidin Alexa 488 or Alexa 568 in 1× PBS (15 μl of stock phalloidin per 200 μl of PBS) was used for 2 hours at room temperature. For immunofluorescence experiments, cells were blocked in 5% bovine serum albumin (BSA) in PBS for 20 min, followed by antibody incubations.</p>
<p>For visualizing titin and myomesin, a live cell extraction was performed to remove cytoplasmic background<sup><xref ref-type="bibr" rid="c36">36</xref></sup>. A cytoskeleton-stabilizing live-cell extraction buffer was made fresh containing 2 mL of stock solution (500 mM 1,4-piperazinediethanesulfonic acid, 25 mM ethylene glycol tetra acetic acid, 25 mM MgCl<sub>2</sub>), 4 mL of 10% polyoxyethylene glycol (PEG; 35,000 molecular weight), 4 ml H2O, and 100 μL of Triton X-100, 10 μM paclitaxel, and 10 μM phalloidin. Cells were treated with this extraction buffer for 1 min, followed by a 1-min wash with wash buffer (extraction buffer without PEG or Triton X-100). Cells were then fixed with 4% PFA for 20 min, followed by antibody labeling. VectaShield with DAPI (H-1200, Vector Laboratories Inc., Burlingame, CA) was used for mounting.</p>
</sec>
<sec id="s4f">
<title>Protein Knockdown</title>
<p>Knockdowns for MYOM were performed using single siRNAs from GE Dharmacon: one to the UTR of MYOM and one to the coding sequence (CDS). Knockdowns for MYH6 were performed using one of two siRNAs for MYH6: one to the UTR and one to the CDS. purchased from GE Dharmacon. Knockdowns for MYH7 were performed using one of two siRNA pools: one of two siRNAs and one of four. Experiments were performed in 96-well culture plates, using the Lipofectamine RNAiMAX reagent and instructions provided by the manufacturer (ThermoFisher, LMRNA015). Following knockdown, cells were re-plated onto glass substrates for 24 hours and fixed for immunofluorescence or lysed for western blotting.</p>
</sec>
<sec id="s4g">
<title>Western Blotting</title>
<p>Gel samples were prepared by mixing cell lysates with LDS sample buffer (Life Technologies, NP0007) and sample reducing buffer (Life Technologies, NP0009) and boiled at 95°C for 5 min. Samples were resolved on Bolt 4–12% gradient Bis-Tris gels (Life Technologies, NW04120BOX). Protein bands were blotted onto a nylon membrane (Millipore). Blots were blocked using 5% nonfat dry milk (NFDM, Research Products International, Mt. Prospect, IL, M17200) in Tris-buffered saline with Tween-20 (TBST, 10X TBS from Corning 46-012-CM, Tween-20 from Sigma P7949). Antibody incubations were also performed in 5% NFDM in TBST. Blots were developed using the Immobilon Chemiluminescence Kit (Millipore, WBKLS0500).</p>
</sec>
<sec id="s4h">
<title>Fluorescence Microscopy</title>
<p>Imaging was performed using either Spinning Disk confocal microscopy or instant Structured Illumination Microscopy (iSIM). The spinning disk images were taken on a Nikon Spinning Disk confocal microscope equipped with Apo TIRF Oil 100X 1.49 NA objective and a Photometrics Prime 95B cMOS monochrome camera, provided by the Nikon Center of Excellence, Vanderbilt University. Images were deconvolved post-acquisition using the FIJI Microvolution software plugin (Microvolution, Cupertino, CA). iSIM imaging was performed with a Visitech iSIM using a Nikon SR HP Apo TIRF 100x oil immersion objective (model number MRD01997) at 1X zoom with NA=1.49. Images were captured using a Hamamatsu ORCA-Fusion Digital CMOS camera (model C14440-20UP) with a 0.1 µm axial step size. Images were deconvolved using Microvolution software (Cupertino, CA) installed in FIJI (Fiji Is Just ImageJ) over 20 iterations.</p>
</sec>
<sec id="s4i">
<title>Platinum Replica Transmission EM of Live Cell Extracted Cells</title>
<p>Adherent plasma membranes from cultured cardimyocytes grown on glass coverslips were detergent extracted. Cells were treated with extraction buffer (2 mL stock buffer (5X Stock Buffer: 500 mM 1,4-piperazinediethanesulfonic acid, 25 mM ethylene glycol tetraacetic acid, 25 mM MgCl2, pH’d and kept at 4C), 4 mL 10% PEG (35,000 MW), 4 mL milliQ H2O, 100 uL of TritonX-100, 10 uM nocodazole, and 10 uM phalloidin) for 30 min, followed by a 1 min wash with wash buffer (2 mL stock buffer, 8 mL milliQ H2O, 10 uM nocodazole, 10 uM phalloidin), followed by fixation (2% PFA, 2% glutaraldehyde) for 20 min. Extracted cells were further sequentially treated with 0.5% OsO4, 1% tannic acid, and 1% uranyl acetate before graded ethanol dehydration and hexamethyldisilazane (HMDS) substitution (LFG Distribution, France). Dried samples were then rotary shadowed with 2 nm of platinum (sputtering) and 4-6 nm of carbon (carbon thread evaporation) using an ACE600 metal coater (Leica Microsystems, Germany). The resultant platinum replica was floated off the glass with hydrofluoric acid (5%), washed several times on distilled water, and picked up on 200 mesh formvar/carbon-coated EM grids. The grids were mounted in a eucentric side-entry goniometer stage of a transmission electron microscope operated at 120 kV (JEOL, Japan), and images were recorded with a Xarosa digital camera (EM-SIS, Germany). Images were processed in Adobe Photoshop to adjust brightness and contrast and presented in inverted contrast.</p>
</sec>
<sec id="s4j">
<title>Statistical Analyses</title>
<p>Analyses comparing three groups (one control and two treatment groups: <xref rid="fig3" ref-type="fig">Figures 3C-G</xref>, <xref ref-type="fig" rid="fig5">5B-F</xref>, <xref ref-type="fig" rid="fig6">6F-H</xref>, <xref ref-type="fig" rid="fig7">7D-H, K-O</xref>, <xref ref-type="fig" rid="fig8">8C-G</xref>, S3A-L, S4A-L, S5A-L, S6A-N, S7A-G, S9A-N, S10A-G, S12A-N, S13A-G, S14A-L, S15A-N) were calculated using a one-way ANOVA. If significant, a post-hoc Tukey test (<xref rid="fig3" ref-type="fig">Figures 3C-F</xref>, <xref ref-type="fig" rid="fig5">5D</xref>, <xref ref-type="fig" rid="fig6">6F-H</xref>, <xref ref-type="fig" rid="fig7">7E, K-M</xref>, <xref ref-type="fig" rid="fig8">8D-G, G</xref>, S4 C, D, F, G, K, S5A-D, S6 C, D, K, M, S7A-D, S9N, S10F, G, S12E, F, N, S13D, F, G, S14A, D, K, S15F, I, L, N) was done. Analyses comparing two groups were calculated using a two-tailed unpaired student’s t test (Figures S8A-L, S11A-L). Graphs were made using SuperPlots<sup><xref ref-type="bibr" rid="c37">37</xref></sup>.</p>
</sec>
</sec>
<sec id="d1e1311" sec-type="supplementary-material">
<title>Supporting information</title>
<supplementary-material id="d1e1434">
<label>Supplemental Figure</label>
<media xlink:href="supplements/523681_file03.pdf"/>
</supplementary-material>
</sec>
</body>
<back>
<ack>
<title>Acknowledgements</title>
<p>We would like to thank the instructors, TAs, and fellow classmates at the Marine Biosciences Laboratory course Deep Learning for Biological Microscopy, especially Jan Funke (Janelia), Dagmar Kainmueller (MDC Berlin), Zachary Whiddon (University of Louisville), and William Patton (Janelia) for assistance on early versions of yoU-Net. We would also like to thank Bryan Millis in the Vanderbilt Biophotonics Center and Kari Seedle in the Vanderbilt Nikon Center for Excellence and the Vanderbilt Center for Imaging Shared Resources (CISR) for experimental and imaging assistance. We thank Vanderbilt’s Program in Developmental Biology, Microtubules and Motors Club and Molecular Biophysics Training Program for project feedback and discussion.</p>
</ack>
<sec id="s5">
<title>Funding</title>
<p>This work was supported by Vanderbilt University T32 5T32HD007502-25 to ACN and ZCS, Vanderbilt University T32 5T32GM008320-32 to JBH, American Heart Predoctoral Fellowship 836090 to JBH and 18PRE33960551 to NT, NIH NIGMS R35 GM125028 to DTB, Vanderbilt University R25 5R25GM062549-18 and American Heart Predoctoral Fellowship 1070985 to ZCS, and F31 HL136081 to AMF.</p>
</sec>
<sec id="s6">
<title>Author Contributions</title>
<p>Conceptualization (ACNC, NT, AMF, DTB), formal analysis (ACNC), funding acquisition (ACNC, DTB, JBH, ZCS, AMF, NT), platinum replica EM imaging (SM, SV), validation (ACNC), investigation (ACNC, JBH, ZCS, NT, AMF), visualization (ACNC), methodology (ACNC), writing-original draft (ACNC, DTB), writing-review and editing (all authors).</p>
</sec>
<ref-list>
<title>References</title>
<ref id="c1"><label>1</label><mixed-citation publication-type="journal"><string-name><surname>Au</surname>, <given-names>Y.</given-names></string-name> <article-title>The muscle ultrastructure: a structural perspective of the sarcomere</article-title>. <source>Cell Mol Life Sci</source> <volume>61</volume>, <fpage>3016</fpage>–<lpage>3033</lpage> (<year>2004</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1007/s00018-004-4282-x">https://doi.org:10.1007/s00018-004-4282-x</ext-link></mixed-citation></ref>
<ref id="c2"><label>2</label><mixed-citation publication-type="journal"><string-name><surname>Hersch</surname>, <given-names>N.</given-names></string-name> <etal>et al.</etal> <article-title>The constant beat: cardiomyocytes adapt their forces by equal contraction upon environmental stiffening</article-title>. <source>Biol Open</source> <volume>2</volume>, <fpage>351</fpage>–<lpage>361</lpage> (<year>2013</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1242/bio.20133830">https://doi.org:10.1242/bio.20133830</ext-link></mixed-citation></ref>
<ref id="c3"><label>3</label><mixed-citation publication-type="journal"><string-name><surname>Wang</surname>, <given-names>Z.</given-names></string-name> <etal>et al.</etal> <article-title>The molecular basis for sarcomere organization in vertebrate skeletal muscle</article-title>. <source>Cell</source> <volume>184</volume>, <fpage>2135</fpage>–<lpage>2150 e2113</lpage> (<year>2021</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1016/j.cell.2021.02.047">https://doi.org:10.1016/j.cell.2021.02.047</ext-link></mixed-citation></ref>
<ref id="c4"><label>4</label><mixed-citation publication-type="journal"><string-name><surname>Eschenhagen</surname>, <given-names>T.</given-names></string-name>, <string-name><surname>Mummery</surname>, <given-names>C.</given-names></string-name> &amp; <string-name><surname>Knollmann</surname>, <given-names>B. C.</given-names></string-name> <article-title>Modelling sarcomeric cardiomyopathies in the dish: from human heart samples to iPSC cardiomyocytes</article-title>. <source>Cardiovasc Res</source> <volume>105</volume>, <fpage>424</fpage>–<lpage>438</lpage> (<year>2015</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1093/cvr/cvv017">https://doi.org:10.1093/cvr/cvv017</ext-link></mixed-citation></ref>
<ref id="c5"><label>5</label><mixed-citation publication-type="journal"><string-name><surname>Harvey</surname>, <given-names>P. A.</given-names></string-name> &amp; <string-name><surname>Leinwand</surname>, <given-names>L. A.</given-names></string-name> <article-title>The cell biology of disease: cellular mechanisms of cardiomyopathy</article-title>. <source>J Cell Biol</source> <volume>194</volume>, <fpage>355</fpage>–<lpage>365</lpage> (<year>2011</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1083/jcb.201101100">https://doi.org:10.1083/jcb.201101100</ext-link></mixed-citation></ref>
<ref id="c6"><label>6</label><mixed-citation publication-type="journal"><string-name><surname>Olsson</surname>, <given-names>M. C.</given-names></string-name>, <string-name><surname>Palmer</surname>, <given-names>B. M.</given-names></string-name>, <string-name><surname>Stauffer</surname>, <given-names>B. L.</given-names></string-name>, <string-name><surname>Leinwand</surname>, <given-names>L. A.</given-names></string-name> &amp; <string-name><surname>Moore</surname>, <given-names>R. L.</given-names></string-name> <article-title>Morphological and functional alterations in ventricular myocytes from male transgenic mice with hypertrophic cardiomyopathy</article-title>. <source>Circ Res</source> <volume>94</volume>, <fpage>201</fpage>–<lpage>207</lpage> (<year>2004</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1161/01.RES.0000111521.40760.18">https://doi.org:10.1161/01.RES.0000111521.40760.18</ext-link></mixed-citation></ref>
<ref id="c7"><label>7</label><mixed-citation publication-type="journal"><string-name><surname>Dlugosz</surname>, <given-names>A. A.</given-names></string-name>, <string-name><surname>Antin</surname>, <given-names>P. B.</given-names></string-name>, <string-name><surname>Nachmias</surname>, <given-names>V. T.</given-names></string-name> &amp; <string-name><surname>Holtzer</surname>, <given-names>H.</given-names></string-name> <article-title>The relationship between stress fiber-like structures and nascent myofibrils in cultured cardiac myocytes</article-title>. <source>J Cell Biol</source> <volume>99</volume>, <fpage>2268</fpage>–<lpage>2278</lpage> (<year>1984</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1083/jcb.99.6.2268">https://doi.org:10.1083/jcb.99.6.2268</ext-link></mixed-citation></ref>
<ref id="c8"><label>8</label><mixed-citation publication-type="journal"><string-name><surname>Sanger</surname>, <given-names>J. W.</given-names></string-name>, <string-name><surname>Mittal</surname>, <given-names>B.</given-names></string-name> &amp; <string-name><surname>Sanger</surname>, <given-names>J. M.</given-names></string-name> <article-title>Formation of myofibrils in spreading chick cardiac myocytes</article-title>. <source>Cell Motil</source> <volume>4</volume>, <fpage>405</fpage>–<lpage>416</lpage> (<year>1984</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1002/cm.970040602">https://doi.org:10.1002/cm.970040602</ext-link></mixed-citation></ref>
<ref id="c9"><label>9</label><mixed-citation publication-type="journal"><string-name><surname>McKenna</surname>, <given-names>N. M.</given-names></string-name>, <string-name><surname>Meigs</surname>, <given-names>J. B.</given-names></string-name> &amp; <string-name><surname>Wang</surname>, <given-names>Y. L.</given-names></string-name> <article-title>Exchangeability of alpha-actinin in living cardiac fibroblasts and muscle cells</article-title>. <source>J Cell Biol</source> <volume>101</volume>, <fpage>2223</fpage>–<lpage>2232</lpage> (<year>1985</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1083/jcb.101.6.2223">https://doi.org:10.1083/jcb.101.6.2223</ext-link></mixed-citation></ref>
<ref id="c10"><label>10</label><mixed-citation publication-type="journal"><string-name><surname>Rhee</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Sanger</surname>, <given-names>J. M.</given-names></string-name> &amp; <string-name><surname>Sanger</surname>, <given-names>J. W.</given-names></string-name> <article-title>The premyofibril: evidence for its role in myofibrillogenesis</article-title>. <source>Cell Motil Cytoskeleton</source> <volume>28</volume>, <fpage>1</fpage>–<lpage>24</lpage> (<year>1994</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1002/cm.970280102">https://doi.org:10.1002/cm.970280102</ext-link></mixed-citation></ref>
<ref id="c11"><label>11</label><mixed-citation publication-type="journal"><string-name><surname>Fenix</surname>, <given-names>A. M.</given-names></string-name> <etal>et al.</etal> <article-title>Muscle-specific stress fibers give rise to sarcomeres in cardiomyocytes</article-title>. <source>Elife</source> <volume>7</volume> (<year>2018</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.7554/eLife.42144">https://doi.org:10.7554/eLife.42144</ext-link></mixed-citation></ref>
<ref id="c12"><label>12</label><mixed-citation publication-type="web"><string-name><surname>Fenix</surname>, <given-names>A. M.</given-names></string-name> <etal>et al.</etal> <article-title>Expansion and concatenation of non-muscle myosin IIA filaments drive cellular contractile system formation during interphase and mitosis</article-title>. <source>Mol Biol Cell</source> (<year>2016</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1091/mbc.E15-10-0725">https://doi.org:10.1091/mbc.E15-10-0725</ext-link></mixed-citation></ref>
<ref id="c13"><label>13</label><mixed-citation publication-type="journal"><string-name><surname>Ostu</surname>, <given-names>N.</given-names></string-name> <article-title>A Threshold Selection Method from Gray-Level Histograms</article-title>. <source>IEEE Transactions on Systems, Man, and Cybernetics</source> <volume>9</volume>, <fpage>62</fpage>–<lpage>66</lpage> (<year>1979</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1109/TSMC.1979.4310076">https://doi.org:10.1109/TSMC.1979.4310076</ext-link></mixed-citation></ref>
<ref id="c14"><label>14</label><mixed-citation publication-type="journal"><string-name><surname>Sezgin</surname>, <given-names>M. S.</given-names></string-name>, B. <article-title>Survey over image thresholding techniques and quantitative performance evaluation</article-title>. <source>Journal of Electronic Imaging</source> <volume>13</volume> (<year>2004</year>). https://doi.org: <pub-id pub-id-type="doi">10.1117/1.1631315</pub-id></mixed-citation></ref>
<ref id="c15"><label>15</label><mixed-citation publication-type="journal"><string-name><surname>Uchida</surname>, <given-names>S.</given-names></string-name> <article-title>Image processing and recognition for biological images</article-title>. <source>Dev Growth Differ</source> <volume>55</volume>, <fpage>523</fpage>–<lpage>549</lpage> (<year>2013</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1111/dgd.12054">https://doi.org:10.1111/dgd.12054</ext-link></mixed-citation></ref>
<ref id="c16"><label>16</label><mixed-citation publication-type="web"><string-name><surname>Ronneberger</surname>, <given-names>O. F.</given-names></string-name>, <string-name><given-names>P.</given-names>; <surname>Brox</surname></string-name>, <article-title>T. U-net: Convolutional networks for biomedical image segmentation</article-title>. <source>International Conference no Medical image computing and computer-assisted intervention</source>, <fpage>234</fpage>–<lpage>241</lpage> (<year>2015</year>). https://doi.org: <pub-id pub-id-type="doi">10.1007/978-3-319-24574-4_28</pub-id></mixed-citation></ref>
<ref id="c17"><label>17</label><mixed-citation publication-type="journal"><string-name><surname>Yin</surname>, <given-names>X. X.</given-names></string-name>, <string-name><surname>Sun</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Fu</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Lu</surname>, <given-names>R.</given-names></string-name> &amp; <string-name><surname>Zhang</surname>, <given-names>Y.</given-names></string-name> <article-title>U-Net-Based Medical Image Segmentation</article-title>. <source>J Healthc Eng</source> <volume>2022</volume>, <fpage>4189781</fpage> (<year>2022</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1155/2022/4189781">https://doi.org:10.1155/2022/4189781</ext-link></mixed-citation></ref>
<ref id="c18"><label>18</label><mixed-citation publication-type="journal"><string-name><surname>Henderson</surname>, <given-names>C. A.</given-names></string-name>, <string-name><surname>Gomez</surname>, <given-names>C. G.</given-names></string-name>, <string-name><surname>Novak</surname>, <given-names>S. M.</given-names></string-name>, <string-name><surname>Mi-Mi</surname>, <given-names>L.</given-names></string-name> &amp; <string-name><surname>Gregorio</surname>, <given-names>C. C.</given-names></string-name> <article-title>Overview of the Muscle Cytoskeleton</article-title>. <source>Compr Physiol</source> <volume>7</volume>, <fpage>891</fpage>–<lpage>944</lpage> (<year>2017</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1002/cphy.c160033">https://doi.org:10.1002/cphy.c160033</ext-link></mixed-citation></ref>
<ref id="c19"><label>19</label><mixed-citation publication-type="journal"><string-name><surname>Chopra</surname>, <given-names>A.</given-names></string-name> <etal>et al.</etal> <article-title>Force Generation via beta-Cardiac Myosin, Titin, and alpha-Actinin Drives Cardiac Sarcomere Assembly from Cell-Matrix Adhesions</article-title>. <source>Dev Cell</source> <volume>44</volume>, <fpage>87</fpage>–<lpage>96 e85</lpage> (<year>2018</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1016/j.devcel.2017.12.012">https://doi.org:10.1016/j.devcel.2017.12.012</ext-link></mixed-citation></ref>
<ref id="c20"><label>20</label><mixed-citation publication-type="journal"><string-name><surname>Ehler</surname>, <given-names>E.</given-names></string-name>, <string-name><surname>Rothen</surname>, <given-names>B. M.</given-names></string-name>, <string-name><surname>Hammerle</surname>, <given-names>S. P.</given-names></string-name>, <string-name><surname>Komiyama</surname>, <given-names>M.</given-names></string-name> &amp; <string-name><surname>Perriard</surname>, <given-names>J. C.</given-names></string-name> <article-title>Myofibrillogenesis in the developing chicken heart: assembly of Z-disk, M-line and the thick filaments</article-title>. <source>J Cell Sci</source> <volume>112</volume> (<issue>Pt10</issue>), <fpage>1529</fpage>–<lpage>1539</lpage> (<year>1999</year>).</mixed-citation></ref>
<ref id="c21"><label>21</label><mixed-citation publication-type="journal"><string-name><surname>Hinson</surname>, <given-names>J. T.</given-names></string-name> <etal>et al.</etal> <article-title>HEART DISEASE. Titin mutations in iPS cells define sarcomere insufficiency as a cause of dilated cardiomyopathy</article-title>. <source>Science</source> <volume>349</volume>, <fpage>982</fpage>–<lpage>986</lpage> (<year>2015</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1126/science.aaa5458">https://doi.org:10.1126/science.aaa5458</ext-link></mixed-citation></ref>
<ref id="c22"><label>22</label><mixed-citation publication-type="journal"><string-name><surname>Taneja</surname>, <given-names>N.</given-names></string-name>, <string-name><surname>Neininger</surname>, <given-names>A. C.</given-names></string-name> &amp; <string-name><surname>Burnette</surname>, <given-names>D. T.</given-names></string-name> <article-title>Coupling to substrate adhesions drives the maturation of muscle stress fibers into myofibrils within cardiomyocytes</article-title>. <source>Mol Biol Cell</source> <volume>31</volume>, <fpage>1273</fpage>–<lpage>1288</lpage> (<year>2020</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1091/mbc.E19-11-0652">https://doi.org:10.1091/mbc.E19-11-0652</ext-link></mixed-citation></ref>
<ref id="c23"><label>23</label><mixed-citation publication-type="journal"><string-name><surname>Murphy</surname>, <given-names>A. C.</given-names></string-name> &amp; <string-name><surname>Young</surname>, <given-names>P. W.</given-names></string-name> <article-title>The actinin family of actin cross-linking proteins - a genetic perspective</article-title>. <source>Cell Biosci</source> <volume>5</volume>, <fpage>49</fpage> (<year>2015</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1186/s13578-015-0029-7">https://doi.org:10.1186/s13578-015-0029-7</ext-link></mixed-citation></ref>
<ref id="c24"><label>24</label><mixed-citation publication-type="journal"><string-name><surname>Ribeiro Ede</surname>, <given-names>A.</given-names>, <suffix>Jr</suffix></string-name>. <etal>et al.</etal> <article-title>The structure and regulation of human muscle alpha-actinin</article-title>. <source>Cell</source> <volume>159</volume>, <fpage>1447</fpage>–<lpage>1460</lpage> (<year>2014</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1016/j.cell.2014.10.056">https://doi.org:10.1016/j.cell.2014.10.056</ext-link></mixed-citation></ref>
<ref id="c25"><label>25</label><mixed-citation publication-type="journal"><string-name><surname>Reedy</surname>, <given-names>M. C.</given-names></string-name> &amp; <string-name><surname>Beall</surname>, <given-names>C.</given-names></string-name> <article-title>Ultrastructure of developing flight muscle in Drosophila. I. Assembly of myofibrils</article-title>. <source>Dev Biol</source> <volume>160</volume>, <fpage>443</fpage>–<lpage>465</lpage> (<year>1993</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1006/dbio.1993.1320">https://doi.org:10.1006/dbio.1993.1320</ext-link></mixed-citation></ref>
<ref id="c26"><label>26</label><mixed-citation publication-type="journal"><string-name><surname>Limouze</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Straight</surname>, <given-names>A. F.</given-names></string-name>, <string-name><surname>Mitchison</surname>, <given-names>T.</given-names></string-name> &amp; <string-name><surname>Sellers</surname>, <given-names>J. R.</given-names></string-name> <article-title>Specificity of blebbistatin, an inhibitor of myosin II</article-title>. <source>J Muscle Res Cell Motil</source> <volume>25</volume>, <fpage>337</fpage>–<lpage>341</lpage> (<year>2004</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1007/s10974-004-6060-7">https://doi.org:10.1007/s10974-004-6060-7</ext-link></mixed-citation></ref>
<ref id="c27"><label>27</label><mixed-citation publication-type="journal"><string-name><surname>Liu</surname>, <given-names>H.</given-names></string-name> <etal>et al.</etal> <article-title>Disassembly of myofibrils and potential imbalanced forces on Z-discs in cultured adult cardiomyocytes</article-title>. <source>Cytoskeleton (Hoboken)</source> <volume>73</volume>, <fpage>246</fpage>–<lpage>257</lpage> (<year>2016</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1002/cm.21298">https://doi.org:10.1002/cm.21298</ext-link></mixed-citation></ref>
<ref id="c28"><label>28</label><mixed-citation publication-type="journal"><string-name><surname>Anderson</surname>, <given-names>B. R.</given-names></string-name> &amp; <string-name><surname>Granzier</surname>, <given-names>H. L.</given-names></string-name> <article-title>Titin-based tension in the cardiac sarcomere: molecular origin and physiological adaptations</article-title>. <source>Prog Biophys Mol Biol</source> <volume>110</volume>, <fpage>204</fpage>–<lpage>217</lpage> (<year>2012</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1016/j.pbiomolbio.2012.08.003">https://doi.org:10.1016/j.pbiomolbio.2012.08.003</ext-link></mixed-citation></ref>
<ref id="c29"><label>29</label><mixed-citation publication-type="journal"><string-name><surname>Herzog</surname>, <given-names>W.</given-names></string-name> <article-title>The multiple roles of titin in muscle contraction and force production</article-title>. <source>Biophys Rev</source> <volume>10</volume>, <fpage>1187</fpage>–<lpage>1199</lpage> (<year>2018</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1007/s12551-017-0395-y">https://doi.org:10.1007/s12551-017-0395-y</ext-link></mixed-citation></ref>
<ref id="c30"><label>30</label><mixed-citation publication-type="journal"><string-name><surname>Yang</surname>, <given-names>K. C.</given-names></string-name> <etal>et al.</etal> <article-title>Novel Adult-Onset Systolic Cardiomyopathy Due to MYH7 E848G Mutation in Patient-Derived Induced Pluripotent Stem Cells</article-title>. <source>JACC Basic Transl Sci</source> <volume>3</volume>, <fpage>728</fpage>–<lpage>740</lpage> (<year>2018</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1016/j.jacbts.2018.08.008">https://doi.org:10.1016/j.jacbts.2018.08.008</ext-link></mixed-citation></ref>
<ref id="c31"><label>31</label><mixed-citation publication-type="journal"><string-name><surname>Lyons</surname>, <given-names>G. E.</given-names></string-name>, <string-name><surname>Schiaffino</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Sassoon</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Barton</surname>, <given-names>P.</given-names></string-name> &amp; <string-name><surname>Buckingham</surname>, <given-names>M.</given-names></string-name> <article-title>Developmental regulation of myosin gene expression in mouse cardiac muscle</article-title>. <source>J Cell Biol</source> <volume>111</volume>, <fpage>2427</fpage>–<lpage>2436</lpage> (<year>1990</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1083/jcb.111.6.2427">https://doi.org:10.1083/jcb.111.6.2427</ext-link></mixed-citation></ref>
<ref id="c32"><label>32</label><mixed-citation publication-type="journal"><string-name><surname>Geach</surname>, <given-names>T. J.</given-names></string-name>, <string-name><surname>Hirst</surname>, <given-names>E. M.</given-names></string-name> &amp; <string-name><surname>Zimmerman</surname>, <given-names>L. B.</given-names></string-name> <article-title>Contractile activity is required for Z-disc sarcomere maturation in vivo</article-title>. <source>Genesis</source> <volume>53</volume>, <fpage>299</fpage>–<lpage>307</lpage> (<year>2015</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1002/dvg.22851">https://doi.org:10.1002/dvg.22851</ext-link></mixed-citation></ref>
<ref id="c33"><label>33</label><mixed-citation publication-type="journal"><string-name><surname>Rudy</surname>, <given-names>D. E.</given-names></string-name>, <string-name><surname>Yatskievych</surname>, <given-names>T. A.</given-names></string-name>, <string-name><surname>Antin</surname>, <given-names>P. B.</given-names></string-name> &amp; <string-name><surname>Gregorio</surname>, <given-names>C. C.</given-names></string-name> <article-title>Assembly of thick, thin, and titin filaments in chick precardiac explants</article-title>. <source>Dev Dyn</source> <volume>221</volume>, <fpage>61</fpage>–<lpage>71</lpage> (<year>2001</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1002/dvdy.1125">https://doi.org:10.1002/dvdy.1125</ext-link></mixed-citation></ref>
<ref id="c34"><label>34</label><mixed-citation publication-type="journal"><string-name><surname>Schafer</surname>, <given-names>D. A.</given-names></string-name>, <string-name><surname>Hug</surname>, <given-names>C.</given-names></string-name> &amp; <string-name><surname>Cooper</surname>, <given-names>J. A.</given-names></string-name> <article-title>Inhibition of CapZ during myofibrillogenesis alters assembly of actin filaments</article-title>. <source>J Cell Biol</source> <volume>128</volume>, <fpage>61</fpage>–<lpage>70</lpage> (<year>1995</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1083/jcb.128.1.61">https://doi.org:10.1083/jcb.128.1.61</ext-link></mixed-citation></ref>
<ref id="c35"><label>35</label><mixed-citation publication-type="journal"><string-name><surname>Beall</surname>, <given-names>C. J.</given-names></string-name>, <string-name><surname>Sepanski</surname>, <given-names>M. A.</given-names></string-name> &amp; <string-name><surname>Fyrberg</surname>, <given-names>E. A.</given-names></string-name> <article-title>Genetic dissection of Drosophila myofibril formation: effects of actin and myosin heavy chain null alleles</article-title>. <source>Genes Dev</source> <volume>3</volume>, <fpage>131</fpage>–<lpage>140</lpage> (<year>1989</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1101/gad.3.2.131">https://doi.org:10.1101/gad.3.2.131</ext-link></mixed-citation></ref>
<ref id="c36"><label>36</label><mixed-citation publication-type="journal"><string-name><surname>Svitkina</surname>, <given-names>T.</given-names></string-name> <article-title>Imaging cytoskeleton components by electron microscopy</article-title>. <source>Methods Mol Biol</source> <volume>586</volume>, <fpage>187</fpage>–<lpage>206</lpage> (<year>2009</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1007/978-1-60761-376-3_10">https://doi.org:10.1007/978-1-60761-376-3_10</ext-link></mixed-citation></ref>
<ref id="c37"><label>37</label><mixed-citation publication-type="journal"><string-name><surname>Lord</surname>, <given-names>S. J.</given-names></string-name>, <string-name><surname>Velle</surname>, <given-names>K. B.</given-names></string-name>, <string-name><surname>Mullins</surname>, <given-names>R. D.</given-names></string-name> &amp; <string-name><surname>Fritz-Laylin</surname>, <given-names>L. K.</given-names></string-name> <article-title>SuperPlots: Communicating reproducibility and variability in cell biology</article-title>. <source>J Cell Biol</source> <volume>219</volume> (<year>2020</year>). <ext-link ext-link-type="uri" xlink:href="https://doi.org:10.1083/jcb.202001064">https://doi.org:10.1083/jcb.202001064</ext-link></mixed-citation></ref>
</ref-list>
</back>
<sub-article id="sa0" article-type="editor-report">
<front-stub>
<article-id pub-id-type="doi">10.7554/eLife.87065.2.sa3</article-id>
<title-group>
<article-title>eLife Assessment</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Dubois</surname>
<given-names>Nicole C</given-names>
</name>
<role specific-use="editor">Reviewing Editor</role>
<aff>
<institution-wrap>
<institution>Icahn School of Medicine at Mount Sinai</institution>
</institution-wrap>
<city>New York</city>
<country>United States of America</country>
</aff>
</contrib>
</contrib-group>
<kwd-group kwd-group-type="evidence-strength">
<kwd>Solid</kwd>
</kwd-group>
<kwd-group kwd-group-type="claim-importance">
<kwd>Useful</kwd>
</kwd-group>
</front-stub>
<body>
<p>This manuscript describes a <bold>useful</bold> tool for quantitative assessment of sarcomere structures in healthy and perturbed cardiomyocytes grown in vitro. The work is <bold>solid</bold>, and the methods, data and analyses broadly support the claims with only minor weaknesses. The tool will be relevant to biologists working on and interested in obtaining quantitative information on sarcomere structure, function and development.</p>
</body>
</sub-article>
<sub-article id="sa1" article-type="referee-report">
<front-stub>
<article-id pub-id-type="doi">10.7554/eLife.87065.2.sa2</article-id>
<title-group>
<article-title>Reviewer #1 (Public Review):</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<anonymous/>
<role specific-use="referee">Reviewer</role>
</contrib>
</contrib-group>
</front-stub>
<body>
<p>This manuscript by Neininger-Castro and colleagues presents a novel automatic image analysis method for assessing sarcomeres, the basic units of myofibrils and validates this tool in a couple of experimental approaches that interfere with sarcomere assembly in iPSC-cardiomyocytes (iPSC-CM).</p>
<p>Automatic quantification of sarcomeres is definitely something that is useful to the field. I am surprised that there is no reference in the manuscript to SarcTrack, published by Toepfer and colleagues in 2019 (PMID 30700234), which has exactly the same purpose. The advantage of the image analysis software presented in the current manuscript appears to me to be that it can cover both mature sarcomeres and nascent sarcomeres in premyofibrils effectively.</p>
</body>
</sub-article>
<sub-article id="sa2" article-type="referee-report">
<front-stub>
<article-id pub-id-type="doi">10.7554/eLife.87065.2.sa1</article-id>
<title-group>
<article-title>Reviewer #2 (Public Review):</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<anonymous/>
<role specific-use="referee">Reviewer</role>
</contrib>
</contrib-group>
</front-stub>
<body>
<p>Neininger-Castro et al report on their original study entitled &quot;Independent regulation of Z-lines and M-lines during sarcomere assembly in cardiac myocytes revealed by the automatic image analysis software sarcApp&quot;, In this study, the research team developed two software, yoU-Net and sarcApp, that provide new binarization and sarcomere quantification methods. The authors further utilized human induced pluripotent stem cell-derived cardiomyocytes (hiCMs) as their model to verify their software by staining multiple sarcomeric components with and without the treatment of Blebbistatin, a known myosin II activity inhibitor. With the treatment of different Blebbistatin concentrations, the morphology of sarcomeric proteins was disturbed. These disrupted sarcomeric structures were further quantified using sarcApp and the quantification data supported the phenotype. The authors further investigated the roles of muscle myosins in sarcomere assembly by knocking down MYH6, MYH7, or MYOM in hiCMs. The knockdown of these genes did not affect Z-line assembly yet the knockdown of MYOM affected M-line assembly. The authors demonstrated that different muscle myosins participate in sarcomere assembly in different manners.</p>
</body>
</sub-article>
<sub-article id="sa3" article-type="referee-report">
<front-stub>
<article-id pub-id-type="doi">10.7554/eLife.87065.2.sa0</article-id>
<title-group>
<article-title>Reviewer #3 (Public Review):</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<anonymous/>
<role specific-use="referee">Reviewer</role>
</contrib>
</contrib-group>
</front-stub>
<body>
<p>Neininger-Castro and colleagues developed software tools for the quantification of sarcomeres and sarcomere-precursor features in immunostained human induced pluripotent stem cell-derived cardiac myocytes (hiCMs). In the first part they used a deep-learning- based model called a U-Net to construct and train a network for binarization of immunostained cardiomyocyte images. They also wrote graphical user interface (GUI) software that will assist other labs to use this approach and made it publicly available. They did not compare their approach to existing ones, but example from one image suggests their binarization tool outperforms Otsu thresholding binarization.</p>
<p>In the second part they developed a software tool called sarcApp that classifies sarcomere structures in the binarized image as a Z-Line or Z-Body and assigns each to either a myofibril or to stress fibers. The tools can then automatically count and measure multiple features (33 per cell and 24 per myofibril) and report them on a per-cell, per-myofibril, and per- stress fiber basis.</p>
<p>To test the tools they used Blebbistatin to inhibit sarcomere assembly and showed that the sarcApp tool could capture changes in multiple features such as fewer myofibrils, fewer Z-Lines, decreased myofibril persistence, decreased Z-Line length and altered myofibril orientation in the Blebbistatin treated cells. With some changes the tool was also shown to quantify sarcomeres in titin and myomesin stained cardiomyocytes.</p>
<p>Finally they used sarcApp to quantify the changes in sarcomere assembly after siRNA mediated knockout of MYH7, MYH7, or MYOM. The analysis indicates that neither MYH6 nor MYH7 knockdown perturbed the assembly of Z- or M-lines, and that knockdown of MYOM perturbed the A-band/M-Line but not the Z-Line assembly according to features captured by the sarcApp tool.</p>
<p>Overall the authors developed and made publicly available an excellent software tool that will be very useful for labs that are interested in studying sarcomere assembly. Multiple features that are difficult to measure or count manually can be automatically measured by the software quickly and accurately.</p>
<p>There are however some remaining questions about these tools:</p>
<p>
1. The binarization tool which is tailored to sarcomere image binarization appears promising but was not systematically compared with existing approaches. Example from one cell suggests it outperforms Otsu's binarization approach.</p>
<p>
2. How robust is the tool? The tool was tested on images from one type of cardiomyocytes (hiCMs) taken from one lab using Nikon Spinning Disk confocal microscope equipped with Apo TIRF Oil 100X 1.49 NA objective or instant Structured Illumination Microscopy (iSIM), using deconvolution (Microvolution software) and in a specific magnification. It remains to be seen whether the tool would be equally effective with images taken with other microscopy systems, with other cardiomyocytes (chick or neonatal rat), with different magnifications, live imaging, etc. The authors state that this approach is also useful in other situations, but the data is not included in this manuscript.</p>
<p>
3. The tool was developed for evaluation of sarcomere assembly. The authors show that for this application it can detect the perturbation by Blebbistatin, or knockdown of sarcomeric genes. It remains to be seen if this tool is also useful for assessment of sarcomere structure for other questions beside sarcomere assembly and in other sarcomere pathologies.</p>
</body>
</sub-article>
<sub-article id="sa4" article-type="author-comment">
<front-stub>
<article-id pub-id-type="doi">10.7554/eLife.87065.2.sa4</article-id>
<title-group>
<article-title>Author Response</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Neininger-Castro</surname>
<given-names>Abigail C.</given-names>
</name>
<role specific-use="author">Author</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Hayes</surname>
<given-names>James B.</given-names>
<suffix>Jr.</suffix></name>
<role specific-use="author">Author</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Sanchez</surname>
<given-names>Zachary C.</given-names>
</name>
<role specific-use="author">Author</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Taneja</surname>
<given-names>Nilay</given-names>
</name>
<role specific-use="author">Author</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fenix</surname>
<given-names>Aidan M.</given-names>
</name>
<role specific-use="author">Author</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Moparthi</surname>
<given-names>Satish</given-names>
</name>
<role specific-use="author">Author</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Vassilopoulos</surname>
<given-names>Stéphane</given-names>
</name>
<role specific-use="author">Author</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Burnette</surname>
<given-names>Dylan T.</given-names>
</name>
<role specific-use="author">Author</role>
</contrib>
</contrib-group>
</front-stub>
<body>
<p>The following is the authors’ response to the original reviews.</p>
<disp-quote content-type="editor-comment">
<p><bold>Reviewer #1 (Public Review):</bold></p>
<p>This manuscript by Neininger-Castro and colleagues presents a novel automatic image analysis method for assessing sarcomeres, the basic units of myofibrils and validates this tool in a couple of experimental approaches that interfere with sarcomere assembly in iPSCcardiomyocytes (iPSC-CM).</p>
<p>Automatic quantification of sarcomeres is definitely something that is useful to the field. I am surprised that there is no reference in the manuscript to SarcTrack, published by Toepfer and colleagues in 2019 (PMID 30700234), which has exactly the same purpose. The advantage of the image analysis software presented in the current manuscript appears to me to be that it can cover both mature sarcomeres and nascent sarcomeres in premyofibrils effectively.</p>
</disp-quote>
<p>We whole-heartedly disagree that SarcTrack has the exact same purpose as sarcApp. sarcApp measures more than the frequency of actinin2 images, and can measure real-space quantifications of actinin, myomesin, and titin, which has not been done before in this way. However, SarcTrack is an interesting method that we hope many researchers find helpful in their research. SarcTrack is a particle tracker that outputs the dimensions of the objects found, but does not distinguish between Z-Lines and other actinin2-positive structures (Z-Bodies, adhesions). It also does not group these structures into higher order structures such as myofibrils and muscle stress fibers.</p>
<disp-quote content-type="editor-comment">
<p>When going through the manuscript there were a few issues that should be addressed in a revised version of the manuscript:</p>
<p>1. I am a bit puzzled that they took 1.4 um length as a cutoff length for a mature A-band in their quantifications, since the consensus in the field for thick filament length seems to be 1.6 um?</p>
</disp-quote>
<p>We use 1.4 µm as a cutoff length for the length of a Z-Line rather than the A-Band. We believe the reviewer is referring to the width of the A-Band perpendicular to the Z-lines, which is indeed 1.6 µm. However, we are referring to the length of the Z-Lines, which can span anywhere from
1.4 µm to up to 10 or more µm. Thank you for allowing us to make the clarification.</p>
<disp-quote content-type="editor-comment">
<p>1. When doing the knockdown for alpha and beta-myosin heavy chain, respectively, why did they not also do a Western blot for the &quot;other&quot; isoform as well (Figure 7)? We know that iPSCCM express a mixture, so the relatively mild phenotype that they observe in single knockdown experiments may well be due to concomitant upregulation of the expression of the other isoform. In my point of view this should be checked.</p>
</disp-quote>
<p>It is likely that in the single knockdown experiments the other isoform is upregulated, which is why we were careful in stating that neither muscle myosin alone is required for sarcomere formation. We do agree this would be an interesting experiment to check beyond the scope of this manuscript.</p>
<disp-quote content-type="editor-comment">
<p>1. There seems to be a disconnect between the images for myomesin knockdown shown in Figure 8H and the quantification shown in Figure 8I, which makes me wonder whether the image shown in H middle (MYOM1 (1) KD), where the beta-myosin doublets do not seem to be much affected is really representative?</p>
</disp-quote>
<p>The image shown in the middle of H is representative of the mean length of beta-myosin doublets in MYOM1 (1) KD hiCMs. While the beta-myosin doublets are still present and organized, they are significantly shorter. In the zoomed out image, you can appreciate much shorter arrays of beta-myosin doublets that, while extending across the entire cell, are thinner than control cells.</p>
<disp-quote content-type="editor-comment">
<p><bold>Reviewer #2 (Public Review):</bold></p>
<p>Neininger-Castro et al report on their original study entitled &quot;Independent regulation of Z-lines and M-lines during sarcomere assembly in cardiac myocytes revealed by the automatic image analysis software sarcApp&quot;, In this study, the research team developed two software, yoU-Net and sarcApp, that provide new binarization and sarcomere quantification methods. The authors further utilized human induced pluripotent stem cell-derived cardiomyocytes (hiCMs) as their model to verify their software by staining multiple sarcomeric components with and without the treatment of Blebbistatin, a known myosin II activity inhibitor. With the treatment of different Blebbistatin concentrations, the morphology of sarcomeric proteins was disturbed. These disrupted sarcomeric structures were further quantified using sarcApp and the quantification data supported the phenotype. The authors further investigated the roles of muscle myosins in sarcomere assembly by knocking down MYH6, MYH7, or MYOM in hiCMs. The knockdown of these genes did not affect Z-line assembly yet the knockdown of MYOM affected M-line assembly. The authors demonstrated that different muscle myosins participate in sarcomere assembly in different manners.</p>
<p><bold>Reviewer #3 (Public Review):</bold></p>
<p>Neininger-Castro and colleagues developed software tools for the quantification of sarcomeres and sarcomere-precursor features in immunostained human induced pluripotent stem cellderived cardiac myocytes (hiCMs). In the first part they used a deep-learning- based model called a U-Net to construct and train a network for binarization of immunostained cardiomyocyte images. They also wrote graphical user interface (GUI) software that will assist other labs in using this approach and made it publicly available. They did not compare their approach to existing ones, but an example from one image suggests their binarization tool outperforms Otsu thresholding binarization.</p>
<p>In the second part they developed a software tool called sarcApp that classifies sarcomere structures in the binarized image as a Z-Line or Z-Body and assigns each to either a myofibril or to stress fibers. The tools can then automatically count and measure multiple features (33 per cell and 24 per myofibril) and report them on a per-cell, per-myofibril, and per- stress fiber basis.</p>
<p>To test the tools they used Blebbistatin to inhibit sarcomere assembly and showed that the sarcApp tool could capture changes in multiple features such as fewer myofibrils, fewer Z-Lines, decreased myofibril persistence, decreased Z-Line length and altered myofibril orientation in the Blebbistatin treated cells. With some changes the tool was also shown to quantify sarcomeres in titin and myomesin stained cardiomyocytes.</p>
<p>Finally they used sarcApp to quantify the changes in sarcomere assembly after siRNA mediated knockout of MYH7, MYH7, or MYOM. The analysis indicates that neither MYH6 nor MYH7 knockdown perturbed the assembly of Z- or M-lines, and that knockdown of MYOM perturbed the A-band/M-Line but not the Z-Line assembly according to features captured by the sarcApp tool.</p>
<p>Overall the authors developed and made publicly available an excellent software tool that will be very useful for labs that are interested in studying sarcomere assembly. Multiple features that are difficult to measure or count manually can be automatically measured by the software quickly and accurately.</p>
<p>There are however some remaining questions about these tools:</p>
<p>1. The binarization tool which is tailored to sarcomere image binarization appears promising but was not systematically compared with existing approaches.</p>
</disp-quote>
<p>We compared it with the existing approach we used previously in the lab, which was Otsu’s method for binarization. We are not aware of several other binarization approaches to compare to, other than using other machine learning techniques that are less advanced than a U-Net, the current standard in image-to-image translation.</p>
<disp-quote content-type="editor-comment">
<p>1. How robust is the tool? The tool was tested on images from one type of cardiomyocytes
(hiCMs) taken from one lab using Nikon Spinning Disk confocal microscope equipped with Apo TIRF Oil 100X 1.49 NA objective or instant Structured Illumination Microscopy (iSIM), using deconvolution (Microvolution software) and in a specific magnification. It remains to be seen whether the tool would be equally effective with images taken with other microscopy systems, with other cardiomyocytes (chick or neonatal rat), with different magnifications, live imaging, etc.</p>
</disp-quote>
<p>We tested the software with several magnifications, with live imaging, and with other tissues. We did not include the information in the manuscript because the data we tested the software with is for future manuscripts studying different aspects of sarcomere formation and maintenance. sarcApp reliably identifies Z-Lines and sarcomeres with deconvolved widefield fluorescence images of hiCMs and frozen human tissue, and are currently using it to measure zebrafish data for another study. Further, it works for live imaging with an actinin2-GFP (or similar) label. For the titin quantification, we would recommend using only 60-100X magnification, as the titin structures (doublets and rings) are not resolvable at lower magnifications.</p>
<disp-quote content-type="editor-comment">
<p>1. The tool was developed for evaluation of sarcomere assembly. The authors show that for this application it can detect the perturbation by Blebbistatin, or knockdown of sarcomeric genes. It remains to be seen if this tool is also useful for assessment of sarcomere structure for other questions beside sarcomere assembly and in other sarcomere pathologies.</p>
</disp-quote>
<p>While this is beyond the scope of this specific methods paper, we welcome other researchers to use our software for other questions in other pathologies. We are currently doing the same for other manuscripts from our lab.</p>
<disp-quote content-type="editor-comment">
<p><bold>Reviewer #1 (Recommendations For The Authors):</bold></p>
<p>1)&quot;alpha-actinin..., which border the sarcomeric contractile machinery (thin and thick filaments); Z-lines do NOT border thick filaments in a relaxed sarcomere</p>
</disp-quote>
<p>We have removed “(thin and thick filaments)” from the text.</p>
<disp-quote content-type="editor-comment">
<p>1. myomesin targeting siRNAs (gene name MYOM): there are actually three genes encoding for myomesin family members, specify, which one was targeted (I am assuming MYOM1).</p>
</disp-quote>
<p>Thank you for the clarification: we do target MYOM1</p>
<disp-quote content-type="editor-comment">
<p>1. I am not surprised that they found not many mature Z-lines in the absence of both sarcomeric myosins; a similar codependence of assembly of mature Z-discs and the presence of functional thick filaments was previously shown by Geach and colleagues in 2015 (PMID 25845369)</p>
</disp-quote>
<p>Thank you for sharing this manuscript: we have added a reference to it in our study.</p>
<disp-quote content-type="editor-comment">
<p><bold>Reviewer #2 (Recommendations For The Authors):</bold></p>
<p>This work offers the possibility to gain more insights into the process of sarcomere assembly through the advancement in sarcomeric or myofibril structure analyses. However, some clarifications are needed from the authors, please see below for the comments.</p>
<p>1. It is recommended that the authors include the time points for replating and harvesting hiCMs. After replating, the cardiomyocytes require at least three to four days for sarcomeric structures to reform. If the hiCMs were fixed before sarcomere assembly had completed, the staining of sarcomeric proteins including ACTN2 and titin could be compromised and it is difficult to tell if the phenotypes observed were consequences of drug treatments or knockdown of sarcomeric genes or simply because the replating hiCMs were fixed before their sarcomeric structures had fully regrown. It is also recommended that the authors replate hiCMs at a fixed time point to avoid discrepancies in the data.</p>
</disp-quote>
<p>Cardiomyocytes do not require three to four days for sarcomeric structures to re-form, and indeed only require 24 hours, with the first sarcomeres typically appearing at ~6 hours. We and others have published several studies demonstrating this (Fenix et al., eLIfe 2018, Taneja, Neininger and Burnette MBoC 2020, Chen et al. Nature Methods, 2022). While sarcomeres continue to develop and turn over after this time, our lab is interested in the beginning steps of sarcomerogenesis rather than the turnover of mature structures.</p>
<disp-quote content-type="editor-comment">
<p>1. The sarcApp automatically identifies Z-lines and Z-bodies; however, is there an option for the users to set their own thresholds? Some users may select different criterions when quantifying sarcomeres. Moreover, the Z-lines and Z-bodies identified by the software are not always accurate. Can the users modify the list manually in an unbiased way. If this function is not available, the authors may consider adding this function to their software. sarcApp measures Zline and Z-bodies length but does not measure Z-line and Z-bodies width, but sometimes it is also necessary to measure the width.</p>
</disp-quote>
<p>Absolutely, users can modify the thresholds to identify Z-Lines and Z-Bodies. There is not a way for users to modify the list in an unbiased way per se, as editing the list of Z-Lines and Z-Bodies based on non-mathematical measurements is inherently biased, but the user is free to add in other Z-Lines and Z-Bodies as they wish. In this context, “manually” and “unbiased” is mutually exclusive.</p>
<disp-quote content-type="editor-comment">
<p>1. It is recommended that the authors include the original images beside the sarcomeric structures identified by sarcApp (Figure 2A, 2C, 4C-F and more). It would be easier to compare the original Z-lines and Z-bodies with those identified by the software.</p>
</disp-quote>
<p>We have added these in Author response image 1.</p>
<fig id="sa4fig1">
<label>Author response image 1.</label>
<caption>
<title>Uncropped images and merges from Figures 2, 4 and 6, respectively.</title>
</caption>
<graphic mime-subtype="jpg" xlink:href="elife-87065-sa4-fig1.jpg" mimetype="image"/>
</fig>
<disp-quote content-type="editor-comment">
<p>1. The M-line length quantification data in Figure 3G, 5F, and 6H showed different colored-dots labeling n1 to n3, but the authors did not discuss the significance of these symbols.</p>
</disp-quote>
<p>We are not sure what the reviewer means by this statement: there is no significance of the different colored dots other than to mark the biological replicate shown. These graphs were created using SuperPlots, which was not stated in the original methods. It has now been added to the Statistical Analysis section.</p>
<disp-quote content-type="editor-comment">
<p>1. Can the authors elaborate more on the reasons why they treated Blebbistatin at concentrations of 50µM and 100µM. Previous studies showed that 25µM of Blebbistatin was sufficient to delay the transformation of cardiomyocytes (PMID 27072942). Can the authors also comment on why they selected 6 hours, 12 hours, and 24 hours post replating for drug treatment. Moreover, the drug treatment at different time points was only done on ACTN2 but not titin or myomesin.</p>
</disp-quote>
<p>We selected 6, 12, and 24 hours for actinin2 to show the time course of sarcomere formation and to show that sarcomeres are developed by 24 hours, as also mentioned above. We are interested in future studies of the time course of titin and myomesin over time, and are working on it in the lab.</p>
<p>We chose 50 and 100 µM Blebbistatin as these completely blocked sarcomere assembly whereas treatment with 25 µM did not. This manuscript is a methods paper that aims to validate sarcApp and show how it could be used. We did not intend for it to be a comprehensive study of how different concentrations of blebbistatin affects sarcomere assembly.</p>
<p>We are also unsure what the reviewer means by “transformation of cardiomyocytes”. The manuscript with the PMID of 27072942 does not address this issue. The paper is a “review and analyze readmission data for patients who received a continuous flow left ventricular assist device (LVAD)”. We assume the reviewer is referring to differentiation. The model system we developed and published in eLife in 2018 does not use differentiating iPSC cardiac myocytes. The hiCMs we use are terminally differentiated but still immature, as they are more transcriptionally similar to primary fetal myocytes. As such, they do not maintain their sarcomeres when they removed from the 96 well and plated onto a glass coverslip for highresolution microscopy. These assemble sarcomeres within 24 hours with the sarcomeres forming close to the dorsal membrane and then rearrange overtime (e.g., moving from the top of the cell to the bottom) (Fenix et al., eLife 2018). With that said, we do agree with the reviewer that a study of sarcomere assembly in the context of cardiac myocyte differentiation would be a fascinating direction for future studies, and we think sarcApp could facilitate such studies.</p>
<disp-quote content-type="editor-comment">
<p>1. The authors mentioned that the myofibrils of Z-line, titin, and M-line were randomly oriented after Blebbistatin treatments. The myofibrils were randomly oriented for titin and M-line. However, the orientation of Z-line after 50µM Blebbistatin treatment was not necessarily random, only the orientation after 100µM Blebbistatin treatment was randomized. The authors might consider changing bar graph to other types of charts if the orientation was really randomized after quantification.</p>
</disp-quote>
<p>We find that the bar chart is the most informative to us, but users can consider other types of charts in their analyses.</p>
<disp-quote content-type="editor-comment">
<p>1. It is recommended that the authors include images staining ACTN2 at lower magnifications (Figure 1A, 1C). With current images, it is true that yoU-Net can separate Z-lines from Z-bodies yet it is difficult to tell if yoU-Net can still distinguish Z-lines from Z-bodies with larger images or it only applies to a small portion of the image.</p>
</disp-quote>
<p>The yoU-Net can distinguish Z-Lines from Z-Bodies with images of any size, as image size (height vs. width in pixels) does not affect how binarization occurs. During binarization, the only pixel requirement is that the width and height are divisible by 8 (for downsampling purposes). Usually this is not the case with raw images, so the image borders are slightly cropped to make them usable. In terms of resolution, we recommend using 60X-100X objectives on confocal or superresolution data for the clearest results. We have, however, successfully binarized deconvolved widefield images at 100X as well.</p>
<disp-quote content-type="editor-comment">
<p>1. The authors mentioned that the knockdown of MYH7 did not affect Z-lines and M-lines; however, the structures of ACTN2, myomesin, and titin appeared more organized as compared to those in control.</p>
</disp-quote>
<p>We agree that the sarcomeres and myofibrils look slightly more organized, and did mean to state that the knockdown did not negatively affect Z-Lines and M-Lines and have updated the manuscript to be more accurate.</p>
<disp-quote content-type="editor-comment">
<p>1. Please provide the merge images for Fig. 4D, 4E, 6B</p>
</disp-quote>
<p>The merge images for Fig. 4D, 4E, and 6B are included with the original images requested above (point 3)</p>
<disp-quote content-type="editor-comment">
<p>1. In the text, they described&quot; &quot;antibodies to the titin I-band localize to both MSFs and sarcomeres in hiCMs (Figure 4A). Titin forms ring-like structures around the Z-Bodies of MSFs that are closer to the apparent sarcomere transition point (Figure 4A)&quot; However, based on the antibody information they provided, it is not explicitly recognized for N-or C-terminus TITIN. Please provide TTN N-terminus or TTN-C terminus co-stainings with ACTN2 antibody to understand which part of TTN together with ACTN2 forms a Z-Body.</p>
</disp-quote>
<p>The TTN antibody is an N-terminal antibody localizing to the I-Band region of sarcomeres. We agree with the reviewer that a more thorough study of titin will be of interest and we are currently undertaking such a study. However, this is a methods paper presenting a tool. While some of the data we present does point to mechanistic hypotheses, it is beyond the scope of this study to fully characterize titin during sarcomere assembly.</p>
<disp-quote content-type="editor-comment">
<p>1. TITIN doublet was used to indicate a sarcomere in Fig. 4C-D. Moreover, they also used another combination (myomesin and F-ACTIN) to label a sarcomere in Fig. 6D. Can they compare the difference between these two methods or by using these two methods (TITIN doublet) and (myomesin and F-ACTIN), how is the average length of sarcomere? Will the sarcomere length be the same?</p>
</disp-quote>
<p>We noted in the manuscript that due to the organization of titin doublets (wrapping around the ends of Z-Lines) that the average titin doublet will be approximately 0.3 um longer than the ZLine. We did not expect to see a difference in lengths of myomesin M-Lines and mature actinin2 Z-Lines and indeed do not see major differences in the average lengths (between 2.0 and 2.5 um in 24 hour control cells)</p>
<disp-quote content-type="editor-comment">
<p>1. They used siRNA method to knockdown MYH6, MYH7 and MYOM and concluded that the knockdown of these genes did not affect the Z-line assembly. Even though they showed very nice knockdown efficiency of these proteins, they should (1) co-stain MYH6/TITIN/actinin2 and MYH6/ myomesin /actinin2 for Fig. 7C. (2) MYH7/TITIN/actinin2 and MYH7/ myomesin /actinin2 for Fig. 7I. (3) MYOM1/TITIN/actinin2 and MYOM2/TITIN/actinin2 for Fig. 8A. (4) MYH7/MYOM1 and MYH7/MYOM2 for Fig. 8H to make sure the cells they measured were truly knockdownpositive cells,</p>
</disp-quote>
<p>The antibodies for alpha and beta myosin are not very efficient for immunofluorescence, and work best for western blots. We decided also to choose a random subset of the cells on the dish to be sure to eliminate any risk of cherry-picking. While imaging cells on the dish, we looked only at the DAPI nuclear channel and selected 50 cells minimum per dish with only this channel, then imaged the other channels.</p>
<disp-quote content-type="editor-comment">
<p>Minor comments:</p>
<p>1. Well-organized sarcomere structure on DMSO treated cells in Fig.5A and Fig. 6A, but it was disarray in Fig. S3M. Why?</p>
</disp-quote>
<p>Figure S3 shows hiCMs that have only been allowed to spread for 6 hours, which have not formed mature sarcomeres yet, hence the disarray.</p>
<disp-quote content-type="editor-comment">
<p>1. Fig 1A, Fig2B: please label the name of the antibody, not the actin filament</p>
</disp-quote>
<p>We used phalloidin labelling here, which marks actin filaments. We have updated the figure legends to be more clear. Thank you!</p>
<disp-quote content-type="editor-comment">
<p>1. Fig. 7I: actinin2 instead of actinin</p>
</disp-quote>
<p>Thank you for catching this! We have fixed it.</p>
<disp-quote content-type="editor-comment">
<p><bold>Reviewer #3 (Recommendations For The Authors):</bold></p>
<p>Testing the app using images shot by other microscopy systems, magnifications, and cardiomyocytes from other species, as noted in the public review above, should make the app even more wildly useful.</p>
<p>A more formal head-to-head comparison with other approaches will be more convincing in showing the new tool is superior</p>
<p>I also think that a more detailed protocol for using the app will help other investigators.</p>
<p>The app counts and measures many features, but it is not always clear how and using what algorithm these are measured. Including these details in a protocol or even as comments in the code will be very helpful for others.</p>
</disp-quote>
<p>The protocol found on the public GitHub for the app will help other investigators to download, use, and understand the application. We have received contact from researchers who have been able to use the application without assistance from us, which is a good sign that the application is user-friendly and that the online protocol is sufficient.</p>
</body>
</sub-article>
</article>