article: To survive starvation, the bacterium Bacillus subtilis forms durable spores. The initial step of sporulation is asymmetric cell division, leading to a large mother-cell and a small forespore compartment. After division is completed and the dividing septum is thinned, the mother cell engulfs the forespore in a slow process based on cell-wall degradation and synthesis. However, recently a new cell-wall independent mechanism was shown to significantly contribute, which can even lead to fast engulfment in 60 of the cases when the cell wall is completely removed. In this backup mechanism, strong ligand-receptor binding between mother-cell protein SpoIIIAH and forespore-protein SpoIIQ leads to zipper-like engulfment, but quantitative understanding is missing. In our work, we combined fluorescence image analysis and stochastic Langevin simulations of the fluctuating membrane to investigate the origin of fast bistable engulfment in absence of the cell wall. Our cell morphologies compare favorably with experimental time-lapse microscopy, with engulfment sensitive to the number of SpoIIQ-SpoIIIAH bonds in a threshold-like manner. By systematic exploration of model parameters, we predict regions of osmotic pressure and membrane-surface tension that produce successful engulfment. Indeed, decreasing the medium osmolarity in experiments prevents engulfment in line with our predictions. Forespore engulfment may thus not only be an ideal model system to study decision-making in single cells, but its biophysical principles are likely applicable to engulfment in other cell types, e. g. during phagocytosis in eukaryotes.
  To survive starvation the Gram-positive bacterium Bacillus subtilis develops durable spores among other survival strategies [1]. During sporulation, bacteria go through a costly developmental process under limited energy resources. The initial morphological step of sporulation is asymmetric cell division, resulting in a large mother-cell and a small forespore compartment [2]. Subsequently, the dividing septum is largely degraded and the mother-cell membrane moves around the forespore. This membrane movement is similar to phagocytosis whereby immune cells clear our bodies from pathogens and other particles [3], [4]. Finally, the engulfed forespore matures into a spore and the mother cell lyzes for its release. The origin of the engulfment force has been a topic of current research [5]–[11]. Cell-wall degradation and new cell-wall deposition were shown to play a significant role in advancing the mother-cell membrane leading edge. Strikingly, when the cell wall is enzymatically removed engulfment still occurs, surprisingly taking only 1–2 min compared to 45 min with the cell wall (see Fig. 1, Movie S1) [8]. Furthermore, engulfment is successful in 60 of cells while the remaining 40 retract. This observation raises questions on the origin of bistability and decision-making in relatively simple systems under severe energy limitations. In the absence of the cell wall, migration of the mother-cell membrane around the forespore depends on the two membrane proteins that bind each other with high affinity [12], constituting a backup mechanism under severe perturbations [8], [13]: SpoIIQ expressed in the forespore and SpoIIIAH expressed in the mother cell [8], [14], [15] (see Fig. 1A–D). To facilitate engulfment a physical mechanism similar to a Brownian ratchet was proposed [8]. Specifically, thermal fluctuations move the leading membrane edge forward, thus establishing new SpoIIQ-SpoIIIAH bonds that prevent backward membrane movement. One striking feature, however, is that the membrane cup surrounding the forespore is very thin (Fig. 1A, top). This either indicates a fast nonequilibrium mechanism for engulfment or additional forces that produce high membrane curvatures around the cup' s neck region. Even though modeling of similar processes such as membrane budding and phagocytosis helped us understand the role of physical constraints on engulfment [16]–[19], quantitative modeling of forespore engulfment as a fundamental development process is still missing. Here, using image analysis, Langevin simulations and simple analytical approaches we show that fast forespore engulfment in the absence of the cell wall occurs below 1 min, consistent with out-of-equilibrium dynamics driven by strong SpoIIQ-SpoIIIAH binding. Furthermore, we find physical parameter regimes responsible for bistable engulfment, including the number of bonds necessary for threshold-like engulfment and suitable osmotic pressures. The former prediction matches previously published data, while we successfully tested the latter with time-lapse microscopy. Hence, our model makes testable predictions on the measurable physical parameters leading to fast, energy-efficient engulfment. Forespore engulfment in the absence of the cell wall is thus an ideal system to study phagocytosis-like processes and decision-making in single-cell organisms.
  To better understand the process of engulfment in the absence of the cell wall, we analyzed the volume and surface area of sporulating cells treated with cell-wall removal enzyme (lysozyme) from previously published data [8]. For this purpose we used the semi-automated image-analysis software JFilament [20] (see Materials and Methods). Briefly, JFilament software allows assisted manipulation of active contours that comply with the bright linear structures of the images. Therefore, membranes that were fluorescently stained with FM 4–64 were easily tracked over time. Information about cell-membrane position in the medial focal plane was used to calculate volume and surface area assuming rotational symmetry around the axis connecting mother and forespore center of mass. Upon cell-wall removal, a drastic volume loss of 35 in engulfing mother cells was observed (see Fig 2A). The onset of volume loss for each cell was set to 0 min, and surface area and engulfment measurements were aligned in time based on this time point. Similar analyses showed that no changes in forespore volume were observed (Fig. 2A). However, surface-area analysis during engulfment did reveal a minor reduction of the mother-cell surface by 5–15 (Fig. 2B), depending on the assumed shapes of the progressing membrane in the image analysis. To correlate the drastic volume loss with the onset of engulfment we created kymographs along the fluorescently labeled forespore membranes of the engulfing cells (Fig. 2C and D). Extracting from the kymographs the percentage of forespore engulfment over time revealed that engulfment in the absence of the cell wall occurs on a time scale of 1–2 min as opposed to 45 min with the cell wall [8], [10]. Calculating cross-correlations we found significant anticorrelations between volume and engulfment, and surface area and engulfment (Fig. 2D). No time delays were detected (the minimum of the cross-correlation function is at 0 min) at a sampling rate of 1. 3 frame/min (for details see Materials and Methods). However, under the same experimental conditions no volume or surface-area losses were observed in cells lacking SpoIIQ or SpoIIIAH zipper-molecules after cell-wall removal (see Fig. S1). Therefore, this led us to conclude that the drastic volume loss may play a facilitating role in engulfment. During fast engulfment the mother-cell membrane drastically changes shape. Based on fluorescence intensity measurements, a very thin and tight cup was proposed at the late stage of engulfment 20 min after cell-wall removal (Fig. 3A) [8]. However, this tight cup amounts to very high membrane curvatures. Therefore, we wondered how the thin cup with high membrane-bending energy may emerge from an initial broad cup, as observed in cells prior to cell-wall removal (Fig. 1A, top). To answer this question, we analyzed fluorescence images in two channels: fluorescent FM 4–64 (red) that uniformly binds to all membranes exposed to the medium, and SpoIIIJ-GFP (green), a labeled protein that is only expressed in the mother cell where it is uniformly recruited to the cell membrane [8]. In Fig. 3A the tight-cup model predicts four membrane folds at the mother-forespore boundary in the FM 4–64 channel and three membrane folds in the GFP channel. Likewise, the broad cup model predicts two membrane folds in the FM 4–64 channel and a single membrane in the GFP channel. To quantify the transition from a broad to the proposed thin cup we measured the average intensity of the mother-forespore boundary () and mother-cell membrane intensity far from the boundary (, see Fig. 3A). Ratio is then used to determine the number of membrane folds on the mother-forespore interface. Relative boundary analysis of FM 4–64 intensity revealed that during late-stage engulfment cup shape conforms in between a broad and tight cup with 3. 5 membrane folds on the mother-forespore interface. Similarly, analysis of the noisy GFP channel (Fig. 3C, and Movie S2) suggests that a broad cup indeed undergoes a transition towards a thin cup through significant morphological changes at the onset of engulfment. To quantitatively explain the observed volume loss, membrane morphological changes, and fast bistable engulfment, we implemented a biophysical model of forespore engulfment using Langevin dynamics. We hypothesized that thin cups form due to fast nonadiabatic engulfment away from equilibrium. If, however, it turns out that the engulfment dynamics are not fast enough under biophysical membrane constraints, then additional forces would need to be postulated to produce such high-curvature membrane features. To test our hypothesis we used Langevin dynamics that account for out-of equilibrium processes. In our model the 3D mother-cell membrane is represented by a string of beads assuming rotational symmetry around the -axis, while the forespore membrane was modeled as a hard sphere (Fig. 4). Indeed, experiments show negligible deformation of the forespore during engulfment (Fig. 2A and B). Specifically, the Langevin dynamic equation of the bead at position is given by: (1) where each bead at the position = (,) represents a ribbon of width and length shown in Fig. 4A. The left-hand side of Eq. 1 depends on the drag coefficient [21], with is the effective medium viscosity (see Text S1). On the right-hand side of Eq. 1 we have contributions of membrane bending, stochastic thermal fluctuations, zipper-molecule binding, surface tension, and osmotic pressure. For each term we give a brief description, while the detailed model equations and analytical derivations can be found in Materials and Methods and Text S1, respectively. The membrane-bending force () restores the curved membrane to the equilibrium flat configuration. The stochastic term () is used to simulate thermal fluctuations with an effective temperature representing a driving force for the leading-edge ratchet movement. The amplitude of the membrane thermal fluctuations is chosen to be 15 nm (see Fig. S2) as typically observed for lipid bilayers [22] and red-blood cells [23]. The zipper-molecule force () accounts for the high binding affinity between SpoIIQ and SpoIIIAH [15]. This zipper mechanism is a strong driving force of the membrane leading edge. For the zipper-protein surface density () we initially chose the maximum possible value that corresponds to a single molecule per 100 nm2 of membrane (based on 10 nm for the size of the protein [14], [15] and assumed dense packing). The surface-tension force () is characterized by linear () and nonlinear () surface tensions (see Materials and Methods) [16]. The force due to osmotic pressure () is characterized by the pressure difference () between inner and outer medium. While the surface-tension term causes membrane contraction, the pressure-difference term produces volume expansion. In thermal equilibrium these two forces balance each other. An example of an implemented simulation is shown in Fig. 4B. To validate our simulations we considered a spherical mother cell prior to engulfment, allowing us to quantify membrane fluctuations in thermal equilibrium. Specifically, we used the Langevin equation without the zipper term on the right-hand side of Eq. 1. Obtained numerical results were compared with analytical results of thermal membrane fluctuations [24], [25] and fluctuation spectra [26] (Fig. S3). This validation showed that our 3D model of the mother cell has indeed appropriate membrane biophysical properties. For further details see Materials and Methods. Using the model governed by Eq. 1 we numerically simulated forespore engulfment in real time (Fig. 5, Movie S3). At the beginning of the simulations mother cell and forespore have spherical shape and intersect at a single point of their perimeters as shown in Fig. 4A. To better understand contributions of linear surface tension () and pressure difference () on engulfment we varied them while keeping other parameters constant (for parameters see Text S1). We note that the explored surface tension (100 pN/µm) is smaller than the experimentally observed rupture tension (20 nN/µm) [27], [28] and that explored pressure differences (1000 Pa) are suitable for an osmotically balanced medium (see Materials and Methods). In Fig. 5A we show simulation snapshots at s. Simulations that reached full engulfment earlier than 5 s were terminated and last snapshots were displayed. In Fig. 5B the white dashed line separates successful engulfment (575 Pa) from retraction (575 Pa). In the region of successful engulfment with 200 Pa and 30 pN/m we observed engulfment with thin, tight cups as experimentally observed (Fig. 5A and B). In the region between retraction and thin cups we also observed successful engulfment but with adiabatic broad cups resulting in almost spherical mother cells. In the region of thin cups engulfment is fast, taking only 1–2 s which is even faster than 1–2 min from experiments (Fig. 2). We attributed this time discrepancy to the limiting factor of cell-wall removal in the experimental setup, since residues of the cell wall can prevent significant membrane fluctuations in our simulations, therefore delaying engulfment (Fig. S4). Indeed, movies of engulfment show elongated cells even minutes after addition of lysozyme [8]. For example, in Fig. 1A, top, at min (the onset of volume loss) the engulfing mother cell still has elongated shape even though this is min after lysozyme treatment. To explore possible simulated engulfment times we varied the kinematic parameter that represents the effective medium viscosity. This parameter sets the time scale of engulfment but does not influence the morphology of the engulfing cup (Fig. S5). We found that even for extremely high, experimentally observed, simulated engulfment times are still about an order of magnitude smaller than the experimentally observed engulfment times (Fig. S5). To better understand physical changes of the simulated mother cell during engulfment, we plotted volume and surface area relative to their initial values (Fig. 5C and D). In the parameter region of successful engulfment we observed a mother-cell volume loss of 0. 2 µm3, which corresponds to the forespore volume. However, the experimentally observed mother-cell volume loss (Fig. 2A) is about 2 times higher. Furthermore, simulated surface areas increase in this region, while experimentally observed surface areas decrease slightly (Fig. 2B). This discrepancy may be explained by effects not included in our simulations, e. g. cytosol leakage (see Discussion). Experiments showed that engulfment depends critically on the number of expressed SpoIIQ zipper protein in a threshold-like fashion (Fig. 6) [8]; for SpoIIQ expression levels below a critical value mother-cell membrane retracts while for expression levels above the critical value the forespore is successfully engulfed (Fig. 6C). For wild-type cells with naturally occurring expression-level variation this presumably leads to the observed 60/40 bistable engulfment outcome [8]. To numerically determine the number of SpoIIQ proteins required for successful engulfment, we varied the protein-surface concentration () for model parameters = 50 pN/µm and = 500 Pa (Fig. 6A and B). Since in experiments the expression of SpoIIQ was reduced compared to wild-type cells, we assumed that the pool of SpoIIIAH molecules is large compared to SpoIIQ. Therefore, the limiting factor in engulfment is attributed to the surface concentration of SpoIIQ proteins. In Fig. 6A successful engulfment occurs for SpoIIQ surface densities above 3350 µm−2. Since forespore surface area is 2 µm2 (Fig. 2B), this produces a lower bound of 6700 SpoIIQ molecules in the forespore necessary for engulfment. To further explore the role of physical constraints on the number of critical zipper proteins that lead to successful engulfment, we scanned the parameter space for successful engulfment in Fig. 5B and found lower bounds on the number of molecules indispensable for engulfment (Fig. 6D). This lower bound can be as low as 120 SpoIIQ molecules for 30 pN/µm and 200 Pa and increases to 7200 SpoIIQ molecules for 550 Pa. Therefore, our simulations predict that under certain osmotic conditions and surface tensions, engulfment can occur with only 100 SpoIIQ molecules. Our model predicts that high osmotic pressure differences across the mother-cell membrane lead to swelling, reduced membrane fluctuations and membrane retraction, therefore preventing the completion of engulfment (Fig. 5A and B). To assess the validity of these predictions, we followed the progression of engulfment in absence of the cell wall for different osmotic conditions by time-lapse fluorescence microscopy (see Materials and Methods). Briefly, we resuspended sporulating cells in SMM buffer, either with 0. 5 M sucrose to support protoplast engulfment (Fig. 1A, Movie S1) [8] or without sucrose, leading to an increase in the osmotic pressure difference between the cytoplasm and the extracellular medium. We then stained cell membranes with FM 4–64, added lysozyme to remove the cell wall, and performed time-lapse microscopy as the cell wall was degraded (Fig. 7A). In the presence of 0. 5 M sucrose, more than 50 of the cells completed engulfment as previously reported (Fig. 7B, Movie S1; [8]). However, in absence of sucrose the engulfing membrane retracted for almost all cells, and less than 2 of them completed engulfment (Fig. 7A and B, Movie S4). Furthermore, membrane retraction was accompanied by an increase in mother-cell volume and surface area of 30 and 5, respectively (Fig. 7D and E). Hence, our data confirm the predictions of our biophysical model on the response to changes in medium osmolarity.
  In this work we presented image analysis and modeling of forespore engulfment in the absence of the cell wall. Image analysis showed that engulfment occurs extremely fast (1–2 min compared to 45 min with the cell wall), accompanied by a drastic volume loss (35) of the mother cell. During engulfment the initial broad cup dynamically changes to a thin cup, forming high-curvature membrane folds at the intersection between mother cell and forespore. Using Langevin simulations we showed that a Brownian ratchet model reproduces fast, out-of-equilibrium engulfment. Additionally, we numerically determined regions of engulfment and retraction, and predicted the number of SpoIIQ molecules necessary for successful engulfment. Similar, out-of-equilibrium Brownian ratchet mechanisms were previously used to explain molecular-motor directional movement [29], Lysteria motility from actin comets [30], [31], filapodia protrusion [32], and unidirectional movement of other microscopic objects [33]. Our model makes a number of predictions. The phase diagram in Fig. 5B predicts engulfment success and mother-cell morphology for a wide range of surface tensions () and pressure differences (). For example, for a given SpoIIQ surface density, high surface tension restricts the engulfment region while high osmotic pressure prevents engulfment. To test this prediction we increased the osmotic pressure difference by lowering the osmolarity of the suspended buffer. Decreased osmolarity indeed caused mother-cell swelling and stopped engulfment in line with our predictions (Fig. 7). Another prediction, such as the need for an excess of mother-cell membrane, might be tested by controlling the production of FapR, the major lipid homeostasis regulator [34], [35]. There are a number of model limitations, which raise interesting issues. Based on our current model there is no stabilizing force for maintaining high curvatures (thin cups) once engulfment is accomplished. Therefore, the restoring constrain forces () and bending force () will flatten high curvatures within 1 s on the length scale of 1 µm (see Fig. S3C and Eq. 11). To estimate the force necessary to prevent membrane flattening we measured restoring forces in the neck region after engulfment is completed (Fig. S6). Typical radial forces were about 10 pN. Three factors can contribute to experimentally observed “snowman-like” shapes at late stages of engulfment (Fig. 1A, top): residues of the mother cell wall, potential outer/inner membrane binding, and lipid/protein sorting. These factors could be investigated experimentally. First, fluorescent cell-wall labeling [36], [37] could rule out possible remnants of cell wall after lysozyme treatment, which may preserve high membrane curvatures (see also Fig. S4). Second, thin cups could also form if an unknown cohesive factor binds to outer and inner cup membrane thus preventing their separation. Third, membrane lipids that localize to high negative curvatures [38]–[40] together with curvature sensing proteins such as SpoVM [41] could contribute to the formation of stable structures preventing membrane flattening. To better explore this stabilizing mechanism we simulated engulfment with membranes that have positive and negative intrinsic curvatures (Fig. S7). We found that successful engulfment proceeds in range (−120 20) µm−1. Interestingly, high negative curvatures prevent membrane flattening by forming tight and thin cups in the neck region typical for the snowman shape (Fig. S7, bottom, left). Therefore, the stability of tight cups after engulfment completion should be a topic of future experimental investigation, leading to better insights and advancement in modeling. Our current model fell short in reproducing the exact volume loss and surface-area conservation. To numerically find parameter regimes that lead to thin cups and observed volume and surface-area changes, we constructed a new albeit less physically grounded model with explicit volume and surface-area constraints (see Fig. S8 and Text S1). Using this model we determined the parameter region in which experimentally observed volume and surface-area changes occur similar to Fig. 5. As a result, engulfment occurs when either volume or surface area is not conserved as previously proposed [42] (see Movie S3). Volume loss may contribute to engulfment by effectively decreasing the surface tension, therefore boosting excess membrane for forespore engulfment. However, by close inspection of sporulating cells, we observed leakage of the cytosol in one of the cells after membrane removal at the onset of volume loss (Fig. S9). This direct loss of cytosol was not explicitly included in our models. One plausible interpretation of volume loss can be attributed to the hypertonic solution of suspended buffer containing 0. 5 M of sucrose and producing an osmotic pressure of 12 atm [8]. Although this pressure is comparable to osmotic pressures inside of bacteria to balance pressures across the membrane, water leakage from bacteria and/or partial lyses after cell-wall removal is expected [43], [44]. Therefore, future experiments may help address this issue. An important future goal is the theoretical understanding of forespore engulfment in the presence of the cell wall, including membrane fission as the last stage of engulfment [45]. High turgor pressure and constraints from the cell wall and septum must provide difficult constraints for the engulfing mother cell, partially explaining the long engulfment time. In the presence of the cell wall, it has been proposed that peptidoglycan hydrolysis and new cell-wall deposition play major roles in the leading-edge membrane movement around the forespore [5]–[7], [9]–[11]. First, membrane proteins SpoIID, SpoIIM, and SpoIIP (DMP) form a complex and localize to the leading edge of the moving membrane. Since SpoIID and SpoIIP degrade peptidoglycans [5], [46] and play an important role for thinning the septum, it has been proposed that the DMP complex is a processive motor for membrane advancement. However, a mechanistic description of this motor is still missing. Second, new cell-wall deposition at the leading edge may provide an additional motor-like mechanism for membrane movement [10], [47]. A similar mechanism was proposed for cytokinesis of fission yeast Schizosaccharomyces pombe, where polymerizing septum fibrils contribute to inward septum ingression. Using Brownian ratchet modeling of this process it was estimated that a single -glucan fibril can exert polymerization force of 10 pN [48]. Together, future modeling will produce a better understanding of the complicated process of forespore engulfment and decipher the contributions from each mechanism towards the total force at the leading membrane edge. In conclusion, our quantitative model of engulfment in absence of the cell wall provided first insights into mother-cell morphologies, such as cup shape, engulfment dynamics, and bistability. Due to simplicity from cell-wall removal, energy limitations, and absence of cytoskeletal cortex, forespore engulfment could present a minimal system for studying bistability and decision-making. Interestingly, bistability and commitment to engulf were previously identified in phagocytosis [18], [49], thus showing similarities to forespore engulfment. This similarity becomes enhanced, if we speculate that the SpoIIQ-SpoIIIAH backup mechanism may instead be the original core mechanism, which may have evolved before the more complex DMP-based mechanism. In fact, the C-terminal domain of SpoIIIAH is homologous to YscJ/FliF protein family forming multimeric rings in type-III secretion system and flagella motors [50], pointing towards an ancient mechanism. General biophysical principles may also apply to other types of engulfment including the penetration of the red blood cell by the malaria parasite [51].
  We used the semi-automated active contour software JFilament [20] available as ImageJ plugin to extract the membrane position over time. All movies analyzed were previously published in [8]. The information about membrane positions obtained from medial focal plane is used to obtain 3D volume and surface area by assuming rotational symmetry around the axis connecting center of mass of mother cell and forespore. Kymographs as in Fig. 2C were created by collecting intensities along the forespore contours using JFilament. Subsequently, pixel angles were determined using pixel position relative to the mother-forespore frame as defined in inset of Fig. 2C. To test the image analysis method, which was used to estimate the number of membrane folds in the neck region of cup (Fig. 3), we measured fluorescence intensities along known single () and double () membranes (Fig. S10). This analysis produced = 2. 3 0. 6 as expected. To compare two signals and that are given at discrete time points, we calculated the cross-correlation function: (2) Here, and are the average signals, and is the total number of discrete time points. A stochastic Langevin equation is used to simulate mother-cell membrane dynamics (see Eq. 1). Simulation time step and distance between neighboring beads were s and = 10 nm, respectively. When the distance between two neighboring beads exceeded the equilibrium distance by 25 the whole contour was rebeaded using a linear-interpolation method [20]. Additionally, membrane-excluding volume is also implemented; whenever the distance between two beads that are not nearest neighbors was less than a repelling radial force of 50 pN was applied to both beads. Analytical details of the model of Eq. 1 are explained in the following sections. To validate our biophysical membrane model, we quantified thermal fluctuations from simulations and compared this with analytical results from [24]–[26] (see Fig. S3). Fourier mode with wave number is calculated as in [26]: (10) with the average membrane radius. Fourier modes are collected during first 10 s of simulations. Autocorrrelation function () of each Fourier mode is calculated and fitted to exponential function () [56], where is the relaxation time for each Fourier mode that satisfies analytical expression [24], [25], [57]: (11) with the reduced membrane tension, the medium viscosity, and the bending stiffness. To validate membrane shapes in thermal equilibrium, we collected 6000 simulated membrane contours and calculated the variance for each Fourier mode, also known as dimensionless spectrum. The analytical expression for dimensionless spectrum of planer membranes is given by [26]: (12) Bacillus subtilis PY79 sporulation was induced by resuspension at 37°C similar to [58], except that the bacteria were grown in 25 LB prior to resuspension rather than in CH medium. Samples were taken two hours and forty-five minutes after resuspension, spun at 9000 rpm for 10 s, and resuspended either in 25 µl of SMM buffer (0. 5 M sucrose, 20 mM maleic acid, 20 mM MgCl2, pH 6. 5) or in the same buffer without sucrose. 10 µl of the resuspended culture were placed on a poly-L-lysine-treated coverslip and mixed with lysozyme and FM 4–64 (final concentrations 1 mg/ml and 5 µg/ml, respectively). Pictures were taken at room temperature, every 45 seconds for one hour, using an Applied Precision optical sectioning microscope equipped with a Photometrics CoolsnapHQ2 camera. Images were deconvolved and analyzed with SoftWoRx version 5. 5 (Applied Precision) and ImageJ.
summary: When the bacterium B. subtilis runs out of food, it undergoes a fundamental development process by which it forms durable spores. Sporulation is initiated by asymmetric cell division after which the larger mother cell engulfs the smaller forespore, followed by spore maturation and release. This survival strategy is so robust that engulfment even proceeds when cells are deprived of their protective cell wall. Under these severe perturbations, 60 of the mother cells still engulf their forespores in only 10 of the normal engulfment time, while the remaining 40 of mother cells withdraw from engulfment. This all-or-none outcome of engulfment suggests decision-making, which was recently also identified in other types of engulfment, e. g. during phagocytosis when immune cells engulf and destroy pathogens. Here, we developed a biophysical model to explain fast bistable forespore engulfment in absence of the cell wall and energy sources. Our discovered principles may prove very general, thus predicting key ingredients of successful engulfment across all kingdoms of life.
keywords: bacteriology
  bacterial physiology
  microbial physiology
  biology and life sciences
  bacterial spores
  microbiology
  cell mechanics
  biophysics
  biomechanics
  biophysical simulations
title: Bistable Forespore Engulfment in Bacillus subtilis by a Zipper Mechanism in Absence of the Cell Wall