<?xml-model href='http://www.tei-c.org/release/xml/tei/custom/schema/relaxng/tei_all.rng' schematypens='http://relaxng.org/ns/structure/1.0'?><TEI xmlns="http://www.tei-c.org/ns/1.0">
	<teiHeader>
		<fileDesc>
			<titleStmt><title level='a'>Subcellular context-specific tuning of actomyosin ring contractility within a common cytoplasm</title></titleStmt>
			<publicationStmt>
				<publisher>bioRxiv</publisher>
				<date>08/09/2024</date>
			</publicationStmt>
			<sourceDesc>
				<bibl> 
					<idno type="par_id">10577406</idno>
					<idno type="doi">10.1101/2024.08.08.607200</idno>
					
					<author>John B Linehan</author><author>Alexandra Zampetaki</author><author>Michael E Werner</author><author>Bryan Heck</author><author>Paul S Maddox</author><author>Sebastian Fürthauer</author><author>Amy S Maddox</author>
				</bibl>
			</sourceDesc>
		</fileDesc>
		<profileDesc>
			<abstract><ab><![CDATA[<title>Summary</title> <p>The non-muscle actomyosin cytoskeleton generates contractile force through the dynamic rearrangement of its constituent parts. Actomyosin rings are a specialization of the non-muscle actomyosin cytoskeleton that drive cell shape changes during division, wound healing, and other events. Contractile rings throughout phylogeny and in a range of cellular contexts are built from conserved components including non-muscle myosin II (NMMII), actin filaments (F-actin), and crosslinking proteins. However, it is unknown whether diverse actomyosin rings close via a single unifying mechanism. To explore how contractile forces are generated by actomyosin rings, we studied three instances of ring closure within the common cytoplasm of the<italic>C. elegans</italic>oogenic germline: mitotic cytokinesis of germline stem cells (GSCs), apoptosis of meiotic compartments, and cellularization of oocytes. We found that each ring type closed with unique kinetics, protein density and abundance dynamics. These measurements suggested that the mechanism of contractile force generation varied across the subcellular contexts. Next, we formulated a physical model that related the forces generated by filament-filament interactions to the material properties of these rings that dictate the kinetics of their closure. Using this framework, we related the density of conserved cytoskeletal proteins anillin and NMMII to the kinematics of ring closure. We fitted model rings to in situ measurements to estimate parameters that are currently experimentally inaccessible, such as the asymmetric distribution of protein along the length of F-actin, which occurs naturally due to differences in the dimensions of the crosslinker and NMMII filaments. Our work predicted that the role of NMMII varies across these ring types, due in part to its distribution along F-actin and motoring. Our model also predicted that the degree of contractility and the impact of ring material properties on contractility differs among ring types.</p>]]></ab></abstract>
		</profileDesc>
	</teiHeader>
	<text><body xmlns="http://www.tei-c.org/ns/1.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink">
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Introduction</head><p>The non-muscle actomyosin cytoskeleton dynamically remodels to drive numerous cell biological and developmental processes <ref type="bibr">(Bear and Haugh, 2014;</ref><ref type="bibr">Martin, 2020)</ref>. The abilities of an actomyosin network to rearrange and to generate and resist forces are determined by the interactions among its molecular components: flexible actin filaments (F-actin), bipolar arrays of myosin motors (non-muscle myosin II; NMMII), and non-motor F-actin crosslinkers including anillin, plastin/fimbrin, filamin, and others. Actomyosin rings are specializations of the plasma membrane-associated cortical cytoskeleton that exist in many contexts throughout animal and fungal cell physiology <ref type="bibr">(Green et al., 2012;</ref><ref type="bibr">Schwayer et al., 2016;</ref><ref type="bibr">Cheffings et al., 2016;</ref><ref type="bibr">Pollard and Wu, 2010)</ref> . Constricting actomyosin rings drive cell expulsion from monolayers, healing of multi-or subcellular wounds, enucleation, and cytokinesis <ref type="bibr">(Schwayer et al., 2016)</ref>. Stable actomyosin rings can maintain cytoplasmic connections among compartments, such as in germline syncytia (see below; <ref type="bibr">(Haglund et al., 2011))</ref>.</p><p>Cytokinesis, when an actomyosin ring separates the two daughter cells during division, is the context in which the functions of many cytoskeletal proteins including NMMII and anillin scaffold proteins have been defined <ref type="bibr">(Green et al., 2012;</ref><ref type="bibr">Kamasaki et al., 2007;</ref><ref type="bibr">Pollard and Wu, 2010;</ref><ref type="bibr">Pollard and Wu, 2010;</ref><ref type="bibr">Schwayer et al., 2016;</ref><ref type="bibr">O'Shaughnessy and Thiyagarajan, 2018)</ref>. Anillins are multidomain proteins that bind the plasma membrane and structural and regulatory elements of the actomyosin cytoskeleton <ref type="bibr">(Maddox et al., 2005;</ref><ref type="bibr">D'Avino, 2009;</ref><ref type="bibr">Piekny and Maddox, 2010;</ref><ref type="bibr">Zhang and Maddox, 2010)</ref>. They contribute to the organization and therefore effectiveness of actomyosin contractile networks <ref type="bibr">(Piekny and Maddox, 2010;</ref><ref type="bibr">D'Avino, 2009;</ref><ref type="bibr">Ku&#269;era et al., 2021)</ref>. NMMII is a motor protein that forms bipolar filamentous aggregates that crosslink and slide F-actin <ref type="bibr">(Glotzer, 2005;</ref><ref type="bibr">Os&#243;rio et al., 2019;</ref><ref type="bibr">Weeds and Lowey, 1971;</ref><ref type="bibr">Niederman and Pollard, 1975)</ref>. NMMII is critical for cytokinesis in most animal and fungal cell types <ref type="bibr">(Mabuchi and Okuno, 1977;</ref><ref type="bibr">De Lozanne and Spudich, 1987;</ref><ref type="bibr">Kitayama et al., 1997;</ref><ref type="bibr">Straight et al., 2003;</ref><ref type="bibr">Ma et al., 2012;</ref><ref type="bibr">Pollard, 2017)</ref> via its abilities to build stresses by motoring along and crosslinking F-actin <ref type="bibr">(Fang et al., 2010;</ref><ref type="bibr">Ma et al., 2012;</ref><ref type="bibr">Palani et al., 2017;</ref><ref type="bibr">Povea-Cabello et al., 2017;</ref><ref type="bibr">Wang et al., 2020;</ref><ref type="bibr">Os&#243;rio et al., 2019)</ref>. NMMII is dispensable for cytokinesis in budding yeast and slime molds due to the dominant mechanical contributions of cell wall deposition and cell motility, respectively <ref type="bibr">(Wloka and Bi, 2012;</ref><ref type="bibr">Lord et al., 2005;</ref><ref type="bibr">Mendes Pinto et al., 2012)</ref>.</p><p>Insight into the mechanisms of actomyosin ring constriction has also come primarily from studying cytokinesis, specifically by measuring closure speed. For example, closure speed scales with ring starting size during C. elegans early embryogenesis and vulval precursor cell division, as well as in the filamentous fungus N. crassa, such that larger cytokinetic rings close faster than smaller cytokinetic rings and duration is relatively conserved <ref type="bibr">(Carvalho et al., 2009;</ref><ref type="bibr">Calvert et al., 2011;</ref><ref type="bibr">Lan et al., 2019;</ref><ref type="bibr">Bourdages et al., 2014;</ref><ref type="bibr">Khaliullin et al., 2018;</ref><ref type="bibr">Lee et al., 2018)</ref>. These and other studies support the idea that ring closure arises from constriction of contractile units arranged in series around the ring <ref type="bibr">(Capco and Bement, 1991;</ref><ref type="bibr">Carvalho et al., 2009)</ref>. By contrast, the distinct kinetics and molecular requirements by different cell types in the same embryo suggest that many different mechanisms exist for contractility <ref type="bibr">(Davies et al., 2018;</ref><ref type="bibr">Ozugergin et al., 2022;</ref><ref type="bibr">reviewed in</ref>  <ref type="bibr">(Ozugergin and Piekny, 2022)</ref>). Thus, while many aspects of actomyosin ring structure and kinetics are widely shared, there exists precedent for mechanistic variations among cell types.</p><p>Various approaches to the phenomenological modeling of ring closure have quantitatively recapitulated closure dynamics of cytokinetic rings in C. elegans, S. pombe (fission yeast), and S. cerevisiae (budding yeast) (reviewed in <ref type="bibr">(Cortes et al., 2018;</ref><ref type="bibr">Pollard, 2014)</ref>). A subset of these models utilize continuum dynamics to relate cytoskeletal protein density to the kinematics of ring closure <ref type="bibr">(Zumdieck et al., 2007;</ref><ref type="bibr">Turlier et al., 2014;</ref><ref type="bibr">Sain et al., 2015)</ref>. These models described the relationship between the ring material properties and dynamics of closure without including the role of motor and crosslinker density within the ring and lack the ability to evaluate the behavior of individual component proteins. By contrast, fully discretized physical agent-based models evaluate ring closure from the perspective of the ring's component proteins. These types of models capture ring closure by simulating each of the components within the ring, and evolving those components position over time <ref type="bibr">(Cortes et al. 2022;</ref><ref type="bibr">Nguyen et al. 2018;</ref><ref type="bibr">Vavylonis et al. 2008)</ref>. A recently developed continuum mechanics model based in active gel theory relates the microscopic interactions of individual filaments to the material properties of a cytoskeletal network <ref type="bibr">(F&#252;rthauer et al., 2019</ref><ref type="bibr">(F&#252;rthauer et al., , 2021;;</ref><ref type="bibr">Foster et al., 2022)</ref>.</p><p>Actomyosin rings are a conserved feature of germlines which, throughout phylogeny, have a complex syncytial architecture wherein many nucleus-containing compartments are interconnected via cytoplasmic bridges <ref type="bibr">(Koch and King, 1969;</ref><ref type="bibr">Mahajan-Miklos and Cooley, 1994;</ref><ref type="bibr">Kumar and Elkouby, 2023;</ref><ref type="bibr">Pepling and Spradling, 1998)</ref>. In the C. elegans syncytial oogenic germline, hundreds of actomyosin rings connect nucleus-containing compartments to a central core of common cytoplasm, known as the rachis <ref type="bibr">(Hirsh et al., 1976;</ref><ref type="bibr">Hall et al., 1999)</ref>. Compartments bear such rings, called rachis bridges, throughout their lifetime. After remaining stably open for tens of hours, a rachis bridge can close if the associated compartment undergoes apoptosis, as do an estimated 50 percent of compartments <ref type="bibr">(Gartner et al., 2008)</ref>. For compartments that do not undergo apoptosis and instead become oocytes, rachis bridges close when nascent oocytes sever from the syncytium <ref type="bibr">(Hall et al., 1999;</ref><ref type="bibr">McCarter et al., 1999)</ref>. Since rings throughout the germline syncytium reside in a common cytoplasm and are comprised of the same conserved cytoskeletal components, it could be hypothesized that they operate via a shared mechanism.</p><p>We explored the general principles of and variations on non-muscle contractility by quantitatively comparing actomyosin ring closure in germline stem cell mitosis, apoptosis and cellularization. We found that these different actomyosin rings within the C. elegans germline exhibit distinct kinetics, as well as profiles of retention of anillin (ANI-1) and NMMII (NMY-2). We utilized a physical framework that relates molecular scale cytoskeletal interactions to the material properties of the ring. Our model suggested that ring closure speed scales not with starting size but with instantaneous size and depends on the material properties of rings, which we found to be dynamic throughout closure and unique to each subcellular context within a large syncytial cell. Our work predicts that differential activity or behavior of NMMII distinguishes the dynamics of these rings from one another.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Methods</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Strain maintenance</head><p>MDX40 C. elegans strain in which mNeonGreen::ANI-1 and mKate2::NMY-2 were expressed from their endogenous loci (ani-1 (mon7[mNeonGreen^3xFlag::ani-1])III x nmy-2(cp52[nmy-2::mKate2+LoxP unc-119(+)LoxP])I; unc-119(ed3) III) was maintained at 20 o Celsius using standard methods <ref type="bibr">(Rehain-Bell et al., 2017)</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Bacteria-impregnated hydrogel -Mounting Method</head><p>Worms were transferred via worm pick into a 1.5 mL tube containing 100 microliters of 0.06% levamisole in S Basal for 10 minutes. 90 -95 microliters of 0.06% levamisole were removed, leaving the worms in 5 -10 mL 0.06% levamisole.</p><p>Photoactivatable hydrogel solution was prepared as described <ref type="bibr">(Burnett et al., 2018)</ref>. 2-Hydroxy-4'-(2-hydroxyethoxy)-2 methylpropiophenome (Sigma Aldrich) photoactivator was suspended in S Basal (final concentrations of all components) to create a 0.001% stock ("photoactivator suspension"), and stored at room temperature, covered in aluminum foil. A 10% suspension of Poly Ethelyne Glycol (Acrylate)2 (PEGDA) (VWR, BioMatrix) was prepared in photoactivator suspension. 2 mL of OP50 liquid culture were spun at 5000 rpm for 10 minutes to pellet the bacteria. The bacteria pellet was resuspended into the photo-activator+PEGDA solution.</p><p>300 microliters of photoactivator+PEGDA+bacteria suspension was transferred first to the wormcontaining tube, and then to a cell culture dish (4-chambered cover glass system #1.5 High-performance cover glass, CellVis, product number: C4-1.5H-N). Worms were pushed to the glass bottom of the dish using an eyelash tool. The dish was then carefully placed on a Speed Light Platinum Gel Documentation system, covered with aluminum foil, and exposed to UV light for 35 seconds to harden the hydrogel. 250 microliters of OP50 liquid culture were pipetted onto the cured gel containing worms.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Fluorescence Imaging</head><p>Imaging was performed using a Nikon A1R Confocal Fluorescence Microscope (Nikon) with 1.27 NA 60x water immersion or a 1.41 NA 60x oil immersion lens and a Gallium arsenide phosphide photomultiplier tube (GaAsP PMT) using Nikon Elements Software and Nikon Perfect Focus.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Ring Dynamic Measurements</head><p>Actomyosin ring closure kinetics were measured by cropping the 4-D fluorescence imaging stack around the ring of interest and then performing a sum projection of the focal planes spanned by the ring. Cropped movies were registered by first manually tracking the center of the actomyosin ring using FIJI <ref type="bibr">(Schindelin et al., 2012)</ref>, and then registering the movies so that the actin ring was in the center of each frame using a custom image registration program written in MATLAB. Rings in registered movies were manually annotated using a custom script in MATLAB. A line was traced around the circumference of the ring as it closed. The pixels covered by the line and the fluorescence intensity within the pixel were recorded. The product of the line length (number of pixels) and pixel size is used as the value for circumference. The total protein was calculated as the sum of the fluorescence intensity of every pixel along the circumference of the ring. The protein density was calculated as the average of the pixels ' fluorescence intensity.</p><p>Measurements of actomyosin ring circumference were likely underestimates of the true circumference since rings were approximated to lie within an optical section (focal plane). The true measurement is within an approximation of n&#8730;2 of the measured value, where n was the number of pixels that intersected the actomyosin ring.</p><p>Measurements of cytoskeletal ring circumference mNeonGreen::ANI-1 and mCherry-NMY-2 CRISPR (LP229) [nmy-2(cp52[nmy-2::mkate2 + LoxP unc-119(+) LoxP]) I; unc-119 (ed3) III] fluorescence intensity were collected for each actomyosin ring throughout closure for as long as accurate measurements could be made. The accuracy of actomyosin ring circumference measurements decayed as the ring completed closure in all ring types as the diminishing opening in actomyosin rings became diffraction-limited. Thus, ring circumference time series were truncated and underestimated the duration of ring closure time. Circumference time series data were smoothed by an averaging filter computed over 5 time points.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Photo-Bleaching Correction</head><p>A photobleaching correction was performed to accurately estimate the change in fluorescence intensity of actomyosin rings over prolonged fluorescence imaging. The decrease in fluorescence intensity was measured in time-lapse image stacks using the Nikon Elements Software package. The measurements of fluorescence intensity were read into MatLab, and a lab-written script was used to correct for fluorescence intensity loss due to photobleaching. The reduction in fluorescence signal was quantified within two regions of interest (ROI); the background, and the dynamic region of interest (where actomyosin rings were measured). The average fluorescence intensity at each frame within each region of interest was tabulated. The background fluorescence intensity was subtracted away from fluorescence intensity values measured from actomyosin rings and the region of interest. An exponential function was fit to the natural logarithm transformed average intensity of the dynamic region of interest to approximate the change in fluorescence intensity due to photobleaching. The ratio method was used to correct photobleaching <ref type="bibr">(Miura, 2020)</ref>. The corrected fluorescence intensity value was calculated,</p><p>Where &#119868;&#119868; &#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888; (&#119905;&#119905;) is the bleach corrected fluorescence intensity at time t, &#119868;&#119868; &#119898;&#119898;&#119888;&#119888;&#119898;&#119898;&#119898;&#119898;&#119898;&#119898;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888; (&#119905;&#119905;) is the measured fluorescence intensity at time t, &#119868;&#119868; &#119887;&#119887;&#119898;&#119898;&#119888;&#119888;&#119887;&#119887;&#119887;&#119887;&#119888;&#119888;&#119888;&#119888;&#119898;&#119898;&#119887;&#119887;&#119888;&#119888; (&#119905;&#119905;) is the fluorescence intensity of the background, and &#119868;&#119868; &#119891;&#119891;&#119891;&#119891;&#119888;&#119888; (&#119905;&#119905;) is the exponential fit of the fluorescence intensity in the dynamic range ROI.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Kymograph and Montage Preparation</head><p>Kymographs were generated from registered actomyosin ring timelapse movies (see Ring Dynamic Measurement) in Image J using the PlugIn KymographBuilder (Plugins -&gt; Kymograph -&gt; Kymograph Builder). Montages were generated using the Fiji tool Make Kymographs (Image -&gt; Stacks -&gt; Make Montage).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Population circumference and protein fluorescence intensity time alignment</head><p>Actomyosin rings within each group were aligned to the final value of each time series; that value corresponds to the last measurement of actomyosin ring circumference made before closure. The time series data for each instance of actomyosin ring closure were truncated to the duration of the shortest acquisition, so that the first time point prior to closure of the population average contained a majority of replicate measurements. The mean value of both circumference and fluorescence intensity at a given time were calculated, along with the standard error. The standard error is calculated as, &#119904;&#119904;&#119905;&#119905;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904; = &#120590;&#120590; &#8730;&#119887;&#119887; , where &#120590;&#120590; is the standard deviation of the measurements at time &#119905;&#119905;, relative to the last point measured, and &#119904;&#119904; the number of measurements at time point &#119905;&#119905; relative to the last point measured.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Determination of the start of actomyosin ring closure</head><p>Actomyosin ring closure onset was determined by plotting the circumference time series and visually determining the time point after which the circumference curve was clearly and consistently decreasing in value.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Average rate of actomyosin ring closure calculation</head><p>The actomyosin ring average circumference change, or closure rate, was calculated as,</p><p>Where &#119889;&#119889;(&#119905;&#119905; = 0) is the circumference at the onset of actomyosin ring closure, &#119889;&#119889;(&#119905;&#119905; = &#119904;&#119904;&#119904;&#119904;&#119904;&#119904;) is the last measured circumference, and &#119879;&#119879; is the time between the onset of closure and the last measured value of circumference in the time series.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Inferred time to complete closure</head><p>The inferred time to complete actomyosin ring closure, assuming that closure rate was constant, is calculated as,</p><p>Where &#119889;&#119889;(&#119905;&#119905; = 0) was the actomyosin ring circumference at the onset of closure, &#119888;&#119888;&#119862;&#119862; &#119888;&#119888;&#119888;&#119888; was the average rate of closure, and &#119879;&#119879; &#119888;&#119888;&#119888;&#119888;&#119898;&#119898;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888;&#119888; is the time in minutes to complete closure.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Actomyosin ring protein content, and protein density time-series alignment to percentage closure</head><p>Actomyosin ring protein content and density were aligned by percentage closure. Percentage closure was calculated as,</p><p>Where &#119889;&#119889;(&#119905;&#119905;) is the actomyosin ring circumference at time &#119905;&#119905; , &#119889;&#119889;(&#119905;&#119905; = 0) is the actomyosin ring circumference at the onset of closure. The circumferential and fluorescence intensity time series (of both ANI-1 and NMY-2) had the same length, so that the index values of elements in each measurement category corresponded to the same point in time. The sum fluorescence intensity (protein content of every actomyosin ring was then sorted into bins (groups) by percentage closure. The bins contained the associated protein content and density values of 5% increments of closure. The average and standard error of each bin, that contained protein content values organized by percentage closure, was taken.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Actomyosin ring size and rate of ingression scaling analysis</head><p>Quantitative relationships between cytokinetic ring circumference at onset of closure and average rate of closure have been reported <ref type="bibr">(Carvalho et al., 2009;</ref><ref type="bibr">Bourdages et al., 2014)</ref>. This scaling law takes the form,</p><p>where &#119888;&#119888;&#119862;&#119862; &#119888;&#119888;&#119888;&#119888; is the average change in circumference throughout closure, &#119889;&#119889; &#119891;&#119891; is the circumference at the onset of closure, &#119898;&#119898; is a proportionality constant with units 1 &#119898;&#119898; , and &#119887;&#119887; is the intercept. In C. elegans blastomeres, cytokinetic ring closure speed scales with ring size as &#119898;&#119898; = 0.0038 and &#119887;&#119887; = 0 <ref type="bibr">(Carvalho et al., 2009)</ref>. When not only blastomeres but also C. elegans vulval precursor cells are considered, this relationship is &#119898;&#119898; = 0.0042, &#119887;&#119887; = -0.02 <ref type="bibr">(Bourdages et al., 2014)</ref>. For each cellular context of actomyosin ring closure in the germline, a range of predicted closure rates was determined using the relation from <ref type="bibr">(Bourdages et al., 2014)</ref>. The lower bound was determined using the mean value and subtracting away the standard error for the value of initial circumference. The upper bound was determined using the mean value plus the standard error for the group initial circumference.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Differential regulation of ANI-1 and NMY-2 turnover rates</head><p>We generated a parameterized model to quantitatively express the hypothesis that the turnover rate of protein on each of the actomyosin rings across cellular contexts of closure are the same. If the regulation of turnover rate is different, then the parameter values would be significantly different across cellular contexts. The hypothesis is expressed quantitatively as the change in fluorescence intensity within the rings,</p><p>Where &#119896;&#119896; &#119888;&#119888;&#119887;&#119887; is the rate of new protein that associates to the ring per unit time and length, and &#119896;&#119896; &#119888;&#119888;&#119891;&#119891;&#119891;&#119891; is the fraction of the total protein in the ring that dissociates, per unit time. The product &#119896;&#119896; &#119888;&#119888;&#119887;&#119887; &#119888;&#119888;(&#119905;&#119905;) is the total rate of new protein incorporated into the ring, and the product &#119896;&#119896; &#119888;&#119888;&#119891;&#119891;&#119891;&#119891; &#119868;&#119868;(&#119905;&#119905;) is the total rate of protein lost from the ring. Integration of the rate of change in fluorescence intensity equation provides an estimate of the fluorescence intensity (total protein content) at time t,</p><p>(3)</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Model Parameterization using Bayesian statistics</head><p>Bayesian analysis was used to determine whether the parameterized model (equation 3) provided an adequate description of the observed fluorescence intensity time-series <ref type="bibr">(Neal, 1993;</ref><ref type="bibr">Hoffman and Gelman, 2011;</ref><ref type="bibr">Chib and Greenberg, 1995)</ref>. The Markov Chain Monte Carlo method with the No U-turn sampling (NUTS) algorithm in the python package PYMC was used to generate distributions for the parameter values &#119896;&#119896; &#119888;&#119888;&#119887;&#119887; , &#119896;&#119896; &#119888;&#119888;&#119891;&#119891;&#119891;&#119891; , for each cellular context of actomyosin ring closure <ref type="bibr">(Patil et al., 2010)</ref>. The initial protein fluorescence intensity is a spurious parameter since the intensity time-series were normalized to the initial value. The PYMC Bayesian analysis generated the joint likelihood that parameter combinations describe the measured fluorescence intensity time-series for ANI-1 and NMY-2 given the model. The likelihoods of the parameter combinations, and the associated distributions of parameter values, provide a way of comparing the most likely parameter combinations for each context of actomyosin ring closure to one another. The model evaluated here took as input the actomyosin ring circumference and estimated the values of the parameters &#119896;&#119896; &#119888;&#119888;&#119887;&#119887; , &#119896;&#119896; &#119888;&#119888;&#119891;&#119891;&#119891;&#119891; using equation ( <ref type="formula">3</ref>).</p><p>The MCMC method with the NUTS algorithm was used as above to determine the distribution of parameters values for &#120572;&#120572; &#119909;&#119909; that described the observed ring closure curves (Figures <ref type="figure">1</ref>, <ref type="figure">5</ref>). A separate analysis was performed, as above, to determine joint likelihood of the parameters &#120572;&#120572; &#119898;&#119898; and &#119881;&#119881; || (Figure <ref type="figure">5</ref>). The rings studied here were aligned to the last time point measured in each trajectory (see Population circumference and protein fluorescence intensity time alignment). The variance among measurements was higher near the onset of closure in apoptotic and cellularizing rachis bridge closure. To account for this, we truncated the first 13 time points from the cellularization data, and the first 10 time points from the apoptosis data (as indicated along the x-axis of Figure <ref type="figure">6</ref>).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Figures</head><p>Graphs, and boxplots were generated in MatLab. All figures were assembled in Adobe Illustrator.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Numerical Programming</head><p>Computer programs were written in Matlab and Jupyter Notebook for Python. Source code can be found on Github under the username JackLinehan in the repository Cellular-Context-Specific-Tuningof-Actomyosin-Ring-Contractility.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Results</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Long-term fluorescence imaging of actomyosin ring closure in the C. elegans oogenic germline</head><p>Cytoskeletal rings are enrichments of F-actin, NMMII, anillin, and diverse other cytoskeletal proteins. These actomyosin rings play a role in a range of cellular processes in addition to canonical cytokinesis. To study contractile force generation in actomyosin rings, we leveraged the co-existence of multiple actomyosin ring closure events in a single syncytial cell: those closing in mitotic, apoptotic, and cellularizing germline compartments in adult C. elegans hermaphrodites (Supplemental Movie 1). We first studied the closure of the germline stem cell (GSC) cytokinetic ring during mitosis in the stem cell niche using fluorescently labeled conserved cytoskeletal components anillin 1 (ANI-1) (Figure <ref type="figure">1A blue</ref>, <ref type="figure">B</ref>, <ref type="figure">B'</ref>) and NMMII heavy chain (NMY-2) (Figure <ref type="figure">1</ref> C, C'), which enrich on the germline rachis lining, including dynamic and stable rings. GSC cytokinetic ring closure commenced as soon as enrichment of ANI-1 and NMY-2 was apparent (Figure1 B, C) (Figure <ref type="figure">1D</ref>, Supplemental Movies 2, 3).</p><p>When germline compartments undergo apoptosis, as occurs in approximately 50% of compartments during pachytene of meiosis I prophase (Figure <ref type="figure">1A</ref>, red) <ref type="bibr">(Gumienny et al., 1999)</ref>, the rachis bridge closes, separating the compartment from the common cytoplasm. Despite the importance of these rings and the knowledge that they bear many of the same components as other contractile rings, the mechanism of their closure is not known. We next used the same methods to study the closure of actomyosin rings that initially connect apoptotic compartments to the common cytoplasm (Figure <ref type="figure">E -G</ref>). The actomyosin rings associated with apoptotic compartments began as anillin-and NMMII-enriched rachis bridges. After remaining at a roughly constant size for many minutes, they closed at a roughly constant rate (Figure <ref type="figure">1</ref> G, Supplemental Movies 4, 5).</p><p>In the short arm of the germline, proximal to the spermatheca, nascent oocytes cellularize from the syncytium via closure of their rachis bridge (Figure <ref type="figure">1A</ref>, green). We also studied rachis bridge closure during the cellularization of nascent oocytes using fluorescently labeled ANI-1 (Figure <ref type="figure">1</ref> H, H') to track circumference (Figure <ref type="figure">1J</ref>) and ANI-1 and NMY-2 abundance (Figure <ref type="figure">1</ref> I, I')(Figure <ref type="figure">1J</ref>, Supplemental Movies 6, 7). As for rings on apoptotic compartments, the rachis bridges that drove cellularization had remained enriched for anillin and NMMII and stably open for some time. Cellularizing rings closed at a roughly constant rate.</p><p>In sum, immobilized live, intact adult hermaphrodite C. elegans allowed us to observe and measure three distinct types of closing rings in the syncytial oogenic germline.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Kinetics of germline actomyosin ring closure</head><p>Insights into the mechanisms of closure by cytoskeletal remodeling and contractile force generation can be gained from measurements of rings 'starting size, closure rate, and closure duration <ref type="bibr">(Carvalho et al., 2009;</ref><ref type="bibr">Calvert et al., 2011)</ref>. Towards our goal of exploring and comparing closure mechanisms of germline rings, we measured ring size over time in the three contexts (Figure <ref type="figure">1 D</ref>, <ref type="figure">G</ref>, <ref type="figure">J</ref>). We first compared the initial circumference of rings among groups and found that GSC cytokinetic rings were the smallest (13.0 &#177; 1.2&#120583;&#120583;&#119898;&#119898;; mean plus or minus standard error; Figure <ref type="figure">2A</ref>), apoptotic germ cell compartment actomyosin rings were approximately 1.4 times larger than GSC cytokinetic rings (17 &#177; 1.6&#120583;&#120583;&#119898;&#119898;; Figure <ref type="figure">2A</ref>), and cellularizing germ cell compartment actomyosin rings were roughly 1.5 times larger than GSC cytokinetic rings (21 &#177; 1.8&#120583;&#120583;&#119898;&#119898;; Fig. <ref type="figure">2A</ref>).</p><p>The average closure rate was the fastest, at 1.46 &#177; 0.1&#120583;&#120583;&#119898;&#119898;/&#119898;&#119898;&#119898;&#119898;&#119904;&#119904;&#119898;&#119898;&#119905;&#119905;&#119904;&#119904;, in GSC cytokinetic rings (Figure <ref type="figure">2A'</ref>). The apoptotic actomyosin ring average closure rate was 0.40 &#177; 0.05&#120583;&#120583;&#119898;&#119898;/&#119898;&#119898;&#119898;&#119898;&#119904;&#119904;&#119898;&#119898;&#119905;&#119905;&#119904;&#119904; (Figure <ref type="figure">3A'</ref>). Cellularizing rings closed at an average rate of 0.44 &#177; 0.1&#120583;&#120583;&#119898;&#119898;/&#119898;&#119898;&#119898;&#119898;&#119904;&#119904;&#119898;&#119898;&#119905;&#119905;&#119904;&#119904; (Figure <ref type="figure">2A'</ref>). We next inferred duration of ring closure by dividing the initial circumference by the average ring closure rate. The average inferred time to close in GSC cytokinetic rings was 9.0 &#177; 0.70 minutes (Figure <ref type="figure">2A''</ref>). The average time for the actomyosin ring to close in apoptotic germ cell compartments was 37.0 &#177; 3.0 minutes (Figure <ref type="figure">2A'')</ref>. The average inferred time to close in cellularizing actomyosin rings was 52.0 &#177; 8.0 minutes (Figure <ref type="figure">2A''</ref>). Thus, the speeds and calculated closure times varied among the three contexts in the germline.</p><p>Cytokinetic ring kinematics have previously been explained by a scaling law relating the initial circumference to closure rate <ref type="bibr">(Carvalho et al., 2009;</ref><ref type="bibr">Calvert et al., 2011;</ref><ref type="bibr">Bourdages et al., 2014;</ref><ref type="bibr">Calvert et al., 2011)</ref>. We explored whether starting ring size was sufficient to predict closure speed for the three C. elegans germline ring types (Figure <ref type="figure">2</ref> B, B')(see Materials and Methods). The measured rate of GSC cytokinetic ring closure indeed fell within the range of closure speed predicted by starting size. This finding suggested that rings in at least three contexts of cytokinesis share mechanisms. By contrast, the measured rates of ring closure for apoptotic or cellularizing rachis bridges were substantially lower than what was predicted by the relationship between initial ring circumference and rate of closure for other C. elegans cell types (Figure <ref type="figure">3B'</ref>). Thus, rachis bridges either close via distinct mechanisms, or close via the same shared mechanism as do cytokinetic rings, but experience greater forces opposing constriction.</p><p>Despite the compatibility of the scaling argument for cytokinesis in GSCs and other cell types, the relationship between the initial ring circumference and closure rate derived from measurements of blastomeres overestimated the measured GSC cytokinetic ring closure rates by roughly a factor of two (predicted rate of closure 0.045 &#120583;&#120583;&#119898;&#119898; &#119898;&#119898; ; measured rate of closure 0.024 &#120583;&#120583;&#119898;&#119898; &#119898;&#119898; ) <ref type="bibr">(Carvalho et al., 2009)</ref>. An adjusted relationship that includes both embryonic and vulval precursor cell cytokinesis data predicts that GSC rings close at 0.03 &#120583;&#120583;&#119898;&#119898; &#119898;&#119898; , which is closer to the measured value of 0.024 &#120583;&#120583;&#119898;&#119898; &#119898;&#119898; (Figure <ref type="figure">2B'</ref>) <ref type="bibr">(Bourdages et al., 2014)</ref>. The earlier quantification of the relationship between ring speed and closure rate implied a direct scaling, while the latter utilized a linear relationship. In total, this analysis revealed that the kinetics of GSC cytokinesis differ from rachis bridge closure in both apoptotic germline compartments and cellularizing nascent oocytes. The finding that germline rachis bridges closed more slowly than predicted by the scaling relationship between initial circumference and closure speed suggests that: 1) forces resisting ring closure are higher in the rachis lining than in GSCs, or 2) rachis bridges and cytokinetic rings contract via different mechanisms.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Time independent ANI-1 and NMY-2 evolution distinguishes the GSC cytokinetic ring from the rachis bridges</head><p>The dynamics of the abundance of ring components over the course of ring closure varies among components and across organisms; compaction and retention dominate some systems, whereas protein density remains relatively constant in other cell types <ref type="bibr">(Wu and Pollard, 2005;</ref><ref type="bibr">Mendes Pinto et al., 2012;</ref><ref type="bibr">Khaliullin et al., 2018;</ref><ref type="bibr">Okada et al., 2021)</ref>. In silico modeling demonstrated that protein abundance dynamics (recruitment, retention, and loss) of key structural components strongly influences the kinetics of actomyosin ring closure <ref type="bibr">(Cortes et al., 2022)</ref>. We reasoned that differences in protein abundance over the course of ring closure among the three ring types would be evidence that the ring types close via unique mechanisms. By contrast, highly similar fold changes in ring component abundance among ring types would suggest a consistent ring-intrinsic mechanism of closure.</p><p>To study the evolution of the three germline ring types independent of the kinematic differences reported above, we measured the fold change of the abundance of the conserved cytoskeletal proteins ANI-1 and NMY-2 as a function of percentage closure. To interpret our protein abundance measurements, we calculated density change for various theoretical regimes of recruitment, retention, and loss (Figure <ref type="figure">4</ref> A, A'). If the net amount of protein is constant throughout closure, we would expect a 5-fold increase in density when the ring is 80% closed (Figure <ref type="figure">4</ref> A, A'; purple). If protein is lost at a rate that is exactly equal to the rate of closure, then the density would remain constant (Figure <ref type="figure">A</ref>, <ref type="figure">A'</ref>; pink). A fold change in density greater than 5x indicates a net gain in the ring, and a fold change that is between 1x and 5x corresponds to a partial loss of protein.</p><p>Normalized fluorescence density and abundance of tagged ANI-1 and NMY-2 varied among the three ring types. GSC cytokinetic rings experienced a greater fold change in density of ANI-1 and NMY-2 than the rachis bridges, which exhibited similar fold changes (Figure <ref type="figure">3</ref> B', C'). ANI-1 abundance initially increased during GSC cytokinesis, and around 50% closure began to decrease. In rachis bridges, ANI-1 was lost throughout closure (Figure <ref type="figure">3B</ref>). NMY-2 abundance increased in GSC cytokinetic rings and decreased in both rachis bridges (Figure <ref type="figure">3C</ref>). Our measurements of ring protein abundance suggest that these dynamics are determined by ring type or cellular context; such a difference in contractile regulation may underlie kinetic differences.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Time dependent ANI-1 and NMY-2 abundance dynamics differ among subcellular contexts of actomyosin ring closure</head><p>The total abundance and average density of ANI-1 and NMY-2 diverged among ring types as they closed (Figure <ref type="figure">3</ref>). Given that the time to complete ring closure varied several fold across subcellular contexts, we next asked how protein abundance varied in time. Furthermore, we leveraged these timeresolved abundance measurements to estimate turnover rates for each protein in each ring type. We used a generalized model for the change in protein content over time:</p><p>where material incorporation was the product of &#119896;&#119896; &#119888;&#119888;&#119887;&#119887; (the rate of protein added per unit space per unit time) and the current circumference of the ring &#119889;&#119889;(&#119905;&#119905;) and material loss was the product of &#119896;&#119896; &#119888;&#119888;&#119891;&#119891;&#119891;&#119891; (the fraction of ring protein shed per unit of time) and the current protein content, &#119868;&#119868;(&#119905;&#119905;) (Figure <ref type="figure">4A</ref>). We calculated the likelihood that the measured fluorescence intensity time-series data were generated by the linear model given a normally distributed error. The Markov Chain Monte Carlo (MCMC) method was used to determine the distribution of parameters &#119896;&#119896; &#119888;&#119888;&#119887;&#119887; and &#119896;&#119896; &#119888;&#119888;&#119891;&#119891;&#119891;&#119891; that fit each average fluorescence intensity over time curve. As expected, the predicted fluorescence intensity time series calculated from the mean parameter values for &#119896;&#119896; &#119888;&#119888;&#119887;&#119887; and &#119896;&#119896; &#119888;&#119888;&#119891;&#119891;&#119891;&#119891; (Figure <ref type="figure">4 B</ref>, <ref type="figure">C</ref>) agreed well with our measured fluorescence intensity time series. The distributions of the parameter value &#119896;&#119896; &#119888;&#119888;&#119887;&#119887; and &#119896;&#119896; &#119888;&#119888;&#119891;&#119891;&#119891;&#119891; were well separated (Figure <ref type="figure">4</ref> B', C'). These results suggest that the regulation of both ANI-1 and NMY-2 differs among the three ring types studied here. Together, our findings indicate that context-dependent actomyosin ring kinetics arise from differences in mechanisms of contractility, independent of any difference that may exist in external force opposing constriction.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Theory for actomyosin rings</head><p>The similarities in ring makeup (all contain F-actin, NMMII, and non-motor crosslinkers including anillin (as reviewed in <ref type="bibr">(Schwayer et al., 2016)</ref>) suggested that a unified theoretical framework could help us explore the differences in their mechanisms of contractility. We used a theoretical physical framework that relates molecular scale interactions to filament network material properties <ref type="bibr">(F&#252;rthauer et al., 2021;</ref><ref type="bibr">F&#252;rthauer et al., 2019)</ref>. The theory derives contractility kinetics from component abundance, beginning with a description of the forces acting on actin-like filaments filament &#119898;&#119898; due to interactions with neighboring filament &#119895;&#119895; via connections made by crosslinkers representing non-motor crosslinkers including anillin, and motors representing NMMII (Figure <ref type="figure">5B</ref>).</p><p>is the force on actin filament &#119898;&#119898; due to interactions via crosslinkers, and</p><p>is the force due to interactions with neighboring filaments via motors. The coefficients &#120574;&#120574; &#119909;&#119909; , &#120574;&#120574; &#119898;&#119898; are the friction constants of crosslinkers &#119909;&#119909; and motors &#119898;&#119898;. The terms &#119907;&#119907; &#8594; &#119891;&#119891; , &#119907;&#119907; &#8594; &#119894;&#119894; describe the velocity of the respective filaments. &#119888;&#119888; &#119909;&#119909; &#65533;&#119904;&#119904; &#119891;&#119891; , &#119904;&#119904; &#119894;&#119894; &#65533;, &#119888;&#119888; &#119898;&#119898; &#65533;&#119904;&#119904; &#119891;&#119891; , &#119904;&#119904; &#119894;&#119894; &#65533; is the density of crosslinkers and motors bound between points &#119904;&#119904; &#119891;&#119891; on filament &#119898;&#119898; and point &#119904;&#119904; &#119894;&#119894; on filament &#119895;&#119895;. The parameters &#119901;&#119901; &#8594; &#119891;&#119891; , &#119901;&#119901; &#8594; &#119894;&#119894; are the unit vectors describing the direction filaments are oriented (Figure <ref type="figure">5A</ref>). &#119881;&#119881; || is the motor stepping speed, taken to be 150 &#119887;&#119887;&#119898;&#119898; &#119898;&#119898; <ref type="bibr">(Nagy et al., 2013)</ref>.</p><p>As previously shown <ref type="bibr">(F&#252;rthauer et al., 2021)</ref>, we can without loss of generality, approximate the crosslink and motor densities by the linear functions of &#119904;&#119904; &#119891;&#119891; , &#119904;&#119904; &#119894;&#119894; , such that &#119888;&#119888; &#119909;&#119909; &#65533;&#119904;&#119904; &#119891;&#119891; , &#119904;&#119904; &#119894;&#119894; &#65533; = &#119888;&#119888; 0 &#119909;&#119909; + &#120572;&#120572; &#119909;&#119909; &#119888;&#119888; 0 &#119909;&#119909; &#65533;&#119904;&#119904; &#119891;&#119891; + &#119904;&#119904; &#119894;&#119894; &#65533;&#119871;&#119871; -1 (6)</p><p>where &#119888;&#119888; 0 &#119909;&#119909;,&#119898;&#119898; is the total amount of crosslinker or motor along the length of the filament, &#120572;&#120572; &#119909;&#119909; &#119888;&#119888; 0 &#119909;&#119909; and &#120572;&#120572; &#119898;&#119898; &#119888;&#119888; 0 &#119898;&#119898; quantifies the spatial asymmetry of the crosslinking and motor proteins 'distribution along the length of filaments (Figure <ref type="figure">5B</ref>). Starting from the above molecular detail equations, and following the coarse graining procedure in NJP, we derived the macroscopic material properties of the network. The network stress obeys</p><p>where &#119907;&#119907; &#8594; is the mean network velocity and &#119868;&#119868; is the identify matrix. &#120578;&#120578; &#119898;&#119898; is the viscosity, given by</p><p>where &#119871;&#119871; = 500&#119904;&#119904;&#119898;&#119898; is the length of the filaments <ref type="bibr">(Wu and Pollard, 2005;</ref><ref type="bibr">Courtemanche et al., 2016;</ref><ref type="bibr">Cortes et al., 2022)</ref>, the crosslinker and motor bridge filaments with a distance &#119877;&#119877; = 13&#119904;&#119904;&#119898;&#119898;, and &#120588;&#120588; is the density of F-actin. The active stress in the network is</p><p>Where we assumed that filaments in the ring would be aligned along the ring.</p><p>The stress and force balance equations are applied to a ring parameterized by the angle &#120579;&#120579;, such that the arc length &#119904;&#119904; is &#119904;&#119904; = &#119904;&#119904;&#120579;&#120579;; where r is the ring radius. We assume that the cytoskeletal network is a thin ribbon with thickness &#8462; (Figure <ref type="figure">5A</ref>). Next, we define the normal vector &#119904;&#119904; ^&#119911;&#119911; and the tangent vector &#119904;&#119904; ^&#119898;&#119898; (Figure <ref type="figure">5A</ref>). We now write the force balance &#120571;&#120571; &#8901; &#120564;&#120564; = 0 for this geometry given the thin shell approximation &#8462; &#8810; &#119904;&#119904;,</p><p>for the tangential direction force balance equation and</p><p>is the radial direction force balance equation. &#119881;&#119881; &#119898;&#119898; is the tangential velocity along the ring circumference opposed by external friction &#120574;&#120574; due to interactions with the plasma membrane and actomyosin cortex. &#119881;&#119881; &#119911;&#119911; is the normal velocity which induces changes in the radius &#119904;&#119904;. We assume that &#120574;&#120574; is sufficiently large so that &#119881;&#119881; &#119898;&#119898; = 0. The dynamics of the ring system are found by deriving the solution to equation ( <ref type="formula">12</ref>),</p><p>where &#119881;&#119881; &#119911;&#119911; = &#119888;&#119888;&#119888;&#119888; &#119888;&#119888;&#119888;&#119888; is the radial velocity at a point along the circumference of the ring. By implementing equations ( <ref type="formula">9</ref>) and ( <ref type="formula">10</ref>) in ( <ref type="formula">13</ref>) we have</p><p>where</p><p>. Equation ( <ref type="formula">14</ref>) relates the changes in protein abundance within rings to the kinematics of ring closure in the absence of an external force.</p><p>We note that the derivation of the model presented in Equations 4-14 assumed isotropic actin filament orientations in the ring, for simplicity. An alternative hypothesis, that the cytoskeleton has strong nematic alignment within the ring, leads to the same final form of Equation <ref type="formula">14</ref>, but</p><p>. For the following, we used &#119860;&#119860; 0 from the isotropic case. Importantly, none of the conclusions of this work depended on this choice (derivation not shown).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>The asymmetric distribution of cytoskeletal components along actin-like filaments is sufficient to recapitulate observed ring closure trajectories</head><p>To understand what microscopic details of network architecture determine the kinetics of actomyosin ring closure, we considered the distribution of the crosslinker and the motor along the length of actin-like filaments represented in the theory. We first assessed a minimal model wherein either the crosslinker or motor was asymmetrically distributed along the length of filaments. To focus our current work on the effects of the abundance of motor and crosslinker, we set the value for friction (g; proportion bound) &#120574;&#120574; &#119909;&#119909; = &#120574;&#120574; &#119898;&#119898; . The equation of motion for a single asymmetry a was given by</p><p>We first considered the case in which the motor is equally distributed along the length of filaments &#120572;&#120572; &#119898;&#119898; &#119888;&#119888; 0 &#119898;&#119898; = 0, while the distribution of the crosslinker was biased towards the filament plus end &#120572;&#120572; &#119909;&#119909; &#119888;&#119888; 0 &#119909;&#119909; &gt; 0. To determine the extent to which a crosslinking protein would need to be asymmetrically distributed to cause closure, we used</p><p>2 (equation 15). We used the MCMC method to determine the distribution of parameter values of &#120572;&#120572; that related the measured protein density dynamics to the observed ring closure curves (Figure <ref type="figure">5A</ref>; black: model predicted curves using the mean &#120572;&#120572; value). For germline stem cell cytokinesis, &#120572;&#120572; &#119909;&#119909; = 0.033 &#177; 0.0006 allowed a good fit to the observed ring closure kinetics (Figure <ref type="figure">6A</ref>, <ref type="figure">A'</ref>). For apoptotic or cellularizing actomyosin ring closure, &#120572;&#120572; &#119909;&#119909; = 0.011 &#177; 0.0003, and &#120572;&#120572; &#119909;&#119909; = 0.005 &#177; 0.0002 (&#119898;&#119898;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904; &#177; &#119904;&#119904;&#119905;&#119905;&#119904;&#119904;) provided sufficient estimates, respectively (Figure <ref type="figure">6 A</ref>, <ref type="figure">A'</ref>). Thus, the MCMC method revealed that measured ring kinetics could be fit well with unique &#120572;&#120572; values for each of the three ring types, providing further support for the idea that the three ring types close via distinct mechanisms.</p><p>Next, we considered the case that the crosslinker was equally distributed along the length of filaments &#120572;&#120572; &#119909;&#119909; &#119888;&#119888; 0 &#119909;&#119909; = 0, while the distribution of the motor was biased towards one end of the filament.</p><p>These conditions resulted in equation ( <ref type="formula">15</ref>) having the form,</p><p>2 , for the ring to close &#120572;&#120572; &#119898;&#119898; &#119888;&#119888; 0 &#119898;&#119898; &lt; 0. Since we had defined &#120574;&#120574; &#119909;&#119909; = 1, the absolute value of this expression was the same as was evaluated above for (&#120572;&#120572; &#119909;&#119909; &#8800; 0, &#120572;&#120572; &#119898;&#119898; = 0), but of opposite sign. The minimal model yields &#120572;&#120572; &lt; 0 when the motor is asymmetrically distributed along F-actin, and &#120572;&#120572; &gt; 0 when the crosslinker is considered asymmetric; the absolute value of the distribution of &#120572;&#120572; was the same.</p><p>In summary, these results validate the relationships among individual component behavior, the material properties of the ring, and the kinetics of closure modeled in equations 4-14 as sufficient to recapitulate the dynamics of actomyosin ring closure.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>The material properties of actomyosin rings change throughout closure</head><p>Active stress and viscosity are material properties of actomyosin networks that dictate the kinetics of their contractility. These important characteristics are difficult to directly measure experimentally, but they can be estimated using our theoretical framework. To understand how the changing protein densities dictate the material properties of these rings, we calculated the active stress and viscosity for each ring type using the mean value of the asymmetry parameter distributions using equations ( <ref type="formula">9</ref>) and (10) (Figure <ref type="figure">6A</ref>', evaluated at &#120572;&#120572; &#119898;&#119898; = 0). (Because asymmetric distribution of either the non-motor crosslinker or the motor was sufficient to fit the model to our kinetics measurements (see Fig. <ref type="figure">6A</ref>), the results of our exploration of the material properties of the ring, are equivalent when either &#120572;&#120572; &#119909;&#119909; &#119888;&#119888; 0 &#119909;&#119909; = 0; &#120572;&#120572; &#119898;&#119898; &#119888;&#119888; 0 &#119898;&#119898; &#8800; 0 or &#120572;&#120572; &#119898;&#119898; &#119888;&#119888; 0 &#119898;&#119898; = 0; &#120572;&#120572; &#119909;&#119909; &#119888;&#119888; 0 &#119909;&#119909; &#8800; 0.) Since not only the starting state, but also in the dynamics of kinetics and protein abundance changed over the course of ring closure, we calculated how active stress and viscosity evolved over time. Active stress and viscosity grew in each of the cellular contexts of closure, with cytokinesis growing 10-fold more than apoptosis, and 50-fold greater than in cellularization (Figure <ref type="figure">6B</ref>, <ref type="figure">B'</ref>). These findings revealed that the material properties of germline rings during closure were dynamic, and unique to each cellular context.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Actomyosin ring contractility exhibits unique dynamics across subcellular contexts</head><p>Our physics-based framework revealed that closure speed is due to two distinct parameters, a ring's size and its material properties (</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>&#120578;&#120578; &#119904;&#119904;</head><p>). Both active stresses and viscosity can be dynamic, so the dynamics of ring closure could result from the constant or dynamically evolving value of one or both of these factors. To calculate and compare contractility, independent of ring size, for our rings, we next determined the contractility, defined as</p><p>, from equation ( <ref type="formula">13</ref>) for each ring type. We found that the GSC cytokinetic ring experienced a ~10% reduction in contractility during closure (Figure <ref type="figure">6</ref> C; blue). Apoptotic compartment rings experienced a similar, but more gradual, reduction in contractility (Figure <ref type="figure">6</ref> C; red). By contrast, in cellularization rings, contractility remained constant throughout closure (Figure <ref type="figure">6C</ref>; green). Together, these results showed that the mechanism of ring closure was the same across the three cellular contexts. The dynamics of closure diverged amongst groups due to the cellular context specific regulation of protein abundance within the rings.</p><p>The magnitude of the active stress may vary across subcellular contexts due to differences in force production by the motor Above, we estimated contractility using measured molecular component densities and a single free parameter, &#120572;&#120572;, the asymmetry of protein distribution along the length of the F-actin like filaments (equation 15; Figure <ref type="figure">6A</ref>, <ref type="figure">A'</ref>). For that work, we held another key parameter constant: the velocity of NMMII motoring (&#119881;&#119881; || ). While good estimates of &#119881;&#119881; || exist <ref type="bibr">(Nagy et al., 2013;</ref><ref type="bibr">Cuda et al., 1997)</ref>, NMMII motoring speed varies across conditions including load <ref type="bibr">(Lord et al., 2005;</ref><ref type="bibr">Wloka and Bi, 2012;</ref><ref type="bibr">Stam et al., 2015;</ref><ref type="bibr">Guo and Guilford, 2006)</ref>. We hypothesized that &#119881;&#119881; || could vary among germline ring types.</p><p>We tested this hypothesis using the MCMC method to calculate the motor protein distribution asymmetry (&#120572;&#120572;) and the motor speed &#119881;&#119881; || . We found that &#120572;&#120572; &#119898;&#119898; , &#119881;&#119881; || , and their product |&#120572;&#120572; &#119898;&#119898; | * &#119881;&#119881; || varied among the subcellular contexts (Figure <ref type="figure">6 D</ref>, <ref type="figure">D'</ref>). |&#120572;&#120572; &#119898;&#119898; | and &#119881;&#119881; || were inversely related (Figure <ref type="figure">6D</ref>), indicating that these two parameters provide the means of tuning motor-based force production. For a given value of |&#120572;&#120572; &#119898;&#119898; |, &#119881;&#119881; || was always higher, and vice versa, for the GSC rings than for either type of rachis bridge, suggesting that NMMII contributes more to the closure of GSC rings than rachis bridges. Together, these results indicate that the motoring speed, and contribution to contractility, of NMMII varies across subcellular contexts.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Discussion</head><p>To explore how contractility emerges from the cytoskeletal interactions within non-muscle actomyosin networks, we studied three instances of actomyosin ring closure within a large syncytial cell. We observed actomyosin rings closing in germline stem cell cytokinesis, the closure of apoptotic compartments, and the cellularization of nascent oocytes. We found that the kinetics of ring closure varied amongst the three groups, and that the conserved cytoskeletal proteins anillin and non-muscle myosin II were differentially retained across contexts. This indicated that the rings closed via unique strategies. We generated a physical model that related protein density dynamics and closure kinetics to the material properties of these rings. This framework allowed us to quantitatively predict the behavior of ring components and compare the predictions across the three rings. Our model revealed that the material properties of actomyosin rings are dynamic throughout closure, and predicted that the role of myosin motoring varies across contexts (Figure <ref type="figure">6</ref>). Together, our measurements and physical modelling NMMII ring closure indicate that ring closure results from cellular context-specific tuning of ring contractility (</p><p>) via regulation of the density of component proteins within rings and the asymmetric distribution of crosslinking and motor proteins along F-actin. Our data suggest that the role of NMMII varies amongst these contexts, potentially due to differences in its ring recruitment kinetics.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Physical modeling and observations of ring closure motivate testable hypothesis of ring contractility</head><p>Our assumption that NMMII abundance scales with the abundance of F-actin, as has been observed <ref type="bibr">(Mendes Pinto et al., 2012;</ref><ref type="bibr">Cortes et al., 2022;</ref><ref type="bibr">Wu and Pollard, 2005;</ref><ref type="bibr">Okada et al., 2021)</ref>, largely dictated the relationship between the ring active stress and viscosity (Figure <ref type="figure">6</ref> B-B'; equations ( <ref type="formula">9</ref>) and ( <ref type="formula">10</ref>)). Additionally, we calculated inferred turnover rates from our measured protein (NMMII and anillin) abundance time series for each context of ring closure. These data indicated that the differential retention of both ANI-1 and NMY-2 across subcellular contexts could be due to cellular context specific turnover of these proteins within rings. Our results motivate several measurements including that of Factin abundance throughout ring closure, and of the turnover (or fluorescence recovery after photobleaching (FRAP) of conserved cytoskeletal proteins. Interestingly, our work indicated that the GSC cytokinetic ring generated more active stress than either of the rachis bridges. Laser ablation of rings can be used to measure recoil velocity and compare the relative magnitude of ring intrinsic forces during closure among ring types. In addition, our modeling work indicated that NMMII motoring contributed more to the closure of the GSC cytokinetic ring than to the closure of rachis bridges. Genetically tuning NMMII motor activity (e.g. via the use of a C. elegans strain expressing temperature-sensitive NMMII mutant) could help test whether GSC cytokinesis is more sensitive to temperature shifting than apoptosis and cellularization.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>The magnitude of contractile stress generated by the rings varied across subcellular contexts</head><p>The rate of ring closure in GSC cytokinesis was accurately predicted using a size-speed scaling relationship with both embryonic and larval somatic C. elegans cell types, indicating that the mechanics of cytokinesis are broadly conserved across cell type and size at least within this species. Interestingly, the same relationship failed to predict the observed closure rates of the two types of rachis bridges. We leveraged our model to test the published explanation for size-speed scaling: the "contractile unit hypothesis" that speed is a function of the number of cytoskeletal segments arranged in series around the ring, corresponding to the starting circumference <ref type="bibr">(Capco and Bement, 1991;</ref><ref type="bibr">Carvalho et al., 2009)</ref>. We found that across the three ring types the rate of closure was dependent not on the initial circumference but on instantaneous size &#119904;&#119904;, the ring radius, and that the proportionality constant was determined by the material properties of the actomyosin ring (&#119898;&#119898; = &#120564;&#120564; &#119860;&#119860; &#120578;&#120578; &#119878;&#119878; , equations 1, and 13) rather than the starting amount of its constituent components. The independent representation of material properties and ring size was necessary for our framework (Equation <ref type="formula">13</ref>) to recapitulate observed kinetics. We can recover the size scaling relationship using equation 15 if larger cytokinetic rings have 1) the same parameter values and 2) similar protein kinetics, as we observed in the GSC cytokinetic ring. The relationship then for ring size and closure rate using our framework is</p><p>Thus, in keeping with the contractile unit or "horse-power" hypothesis, the dimensions of a contractile array are important, but material properties, which vary not only among ring types but over time within a ring type, are also key to explaining contractile speed (Figure <ref type="figure">6B</ref>).</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>The asymmetric distribution of NMMII along F-actin and its motoring tune ring contractility</head><p>Non-motor crosslinkers including anillin can generate contractile forces in reconstituted F-actin networks but only when the level of filament overlap is low <ref type="bibr">(Ku&#269;era et al., 2021)</ref>; in cytokinetic rings Factin is thought to be maximally overlapped <ref type="bibr">(Swulius et al., 2018;</ref><ref type="bibr">reviewed in (Mangione and Gould, 2019)</ref>). Therefore, we assumed that the non-motor crosslinker contributes minimally to contractility. In vivo and in our model, NMMII generates contractile stresses via motoring and contributes to network viscosity by crosslinking F-actin <ref type="bibr">(Wang et al., 2020)</ref>. The contribution of NMMII motoring to ring closure in vivo varies across organismal contexts <ref type="bibr">(Lord et al., 2005;</ref><ref type="bibr">Os&#243;rio et al., 2019;</ref><ref type="bibr">Wloka and Bi, 2012)</ref> and could vary among the three subcellular contexts we explored. Since the magnitude of the force generated by NMMII depends on both its polarized distribution along filaments and its motor activity (Equation <ref type="formula">5</ref>), it could be that either the asymmetric distribution, degree of filament sliding, or a combination of both contributed to tuning the contractility of rings across cellular contexts.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>The asymmetric distribution of crosslinking and motor protein along F-actin is a physical inevitability</head><p>Our work indicates that the asymmetric distribution of crosslinking and motor protein results in contractile stresses within the ring. The interactions between cytoskeletal proteins are known to directly impact network architecture (reviewed in <ref type="bibr">(Kadzik et al., 2020)</ref>). Asymmetry likely naturally occurs given that NMMII forms bundles that span larger distances between actin filaments than many individual crosslinking proteins <ref type="bibr">(Sinard et al., 1989)</ref>. For example, the crosslinking proteins fascin and &#120572;&#120572; -actinin segregate due to differences in size, resulting in mutually exclusive binding of F-actin due to the space between filaments within bundles <ref type="bibr">(Winkelman et al., 2016;</ref><ref type="bibr">Li et al., 2015)</ref>. Thus, it is likely that cytoskeletal proteins occupy mutually exclusive (asymmetric) zones along actin filaments within actomyosin rings.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>During closure, rings become increasingly viscous, which reduces contractility</head><p>Our measurements and modeling suggest that as the density of molecular components increases, the ring becomes increasingly viscous, and therefore less contractile (Figure <ref type="figure">6 B-B'</ref>, <ref type="figure">C</ref>). Thus, the model predicts that the GSC cytokinetic ring and apoptotic compartment ring slow near to 80% closure, as is observed in the cytokinetic ring of the C. elegans zygote <ref type="bibr">(Carvalho et al., 2009;</ref><ref type="bibr">Dorn et al., 2016;</ref><ref type="bibr">Khaliullin et al., 2018;</ref><ref type="bibr">Cortes et al., 2022)</ref>. In contrast, contractility of rachis bridges on cellularizing oocytes was constant throughout closure and was associated with a lower measured density of NMMII and inferred motoring speed. The abundances of ANI-1 and NMY-2 in cellularizing rings trended similarly throughout closure as the homologous proteins' abundances in budding yeast, where cell wall deposition contributes to ring closure. Cellularization occurs on a tapered conical surface; this geometry may contribute to the closure of rachis bridges during cellularization.</p><p>In sum, our theoretical model implementing in vivo measurements reveals that different rings in a syncytium close via the same mechanics but with distinct regulation of protein component density, abundance, and distribution. Our framework makes many testable hypotheses about how protein behavior and abundance dynamics, material properties, and the kinetics of actomyosin ring closure are linked. oocytes (I). C', F', I') Kymographs generated from non-muscle myosin II imaging. D, G, J) Actomyosin ring circumferential time series, &#119898;&#119898;&#119904;&#119904;&#119904;&#119904;&#119904;&#119904; &#177; &#119904;&#119904;&#119904;&#119904;; each actomyosin ring trajectory was aligned to the final point in its trajectory. D) GSC cytokinesis, G) apoptotic compartment rachis bridge closure, J) cellularizing nascent oocyte rachis bridge closure.     </p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>References</head></div></body>
		</text>
</TEI>
