Methodology article | Open | Published:
An integrated enhancement and reconstruction strategy for the quantitative extraction of actin stress fibers from fluorescence micrographs
BMC Bioinformaticsvolume 18, Article number: 268 (2017)
The stress fibers are prominent organization of actin filaments that perform important functions in cellular processes such as migration, polarization, and traction force generation, and whose collective organization reflects the physiological and mechanical activities of the cells. Easily visualized by fluorescence microscopy, the stress fibers are widely used as qualitative descriptors of cell phenotypes. However, due to the complexity of the stress fibers and the presence of other actin-containing cellular features, images of stress fibers are relatively challenging to quantitatively analyze using previously developed approaches, requiring significant user intervention. This poses a challenge for the automation of their detection, segmentation, and quantitative analysis.
Here we describe an open-source software package, SFEX (Stress Fiber Extractor), which is geared for efficient enhancement, segmentation, and analysis of actin stress fibers in adherent tissue culture cells. Our method made use of a carefully chosen image filtering technique to enhance filamentous structures, effectively facilitating the detection and segmentation of stress fibers by binary thresholding. We subdivided the skeletons of stress fiber traces into piecewise-linear fragments, and used a set of geometric criteria to reconstruct the stress fiber networks by pairing appropriate fiber fragments. Our strategy enables the trajectory of a majority of stress fibers within the cells to be comprehensively extracted. We also present a method for quantifying the dimensions of the stress fibers using an image gradient-based approach. We determine the optimal parameter space using sensitivity analysis, and demonstrate the utility of our approach by analyzing actin stress fibers in cells cultured on various micropattern substrates.
We present an open-source graphically-interfaced computational tool for the extraction and quantification of stress fibers in adherent cells with minimal user input. This facilitates the automated extraction of actin stress fibers from fluorescence images. We highlight their potential uses by analyzing images of cells with shapes constrained by fibronectin micropatterns. The method we reported here could serve as the first step in the detection and characterization of the spatial properties of actin stress fibers to enable further detailed morphological analysis.
The stress fibers are prominent assemblies of filamentous actin (F-actin) commonly observed in adherent tissue culture cells. Often considered to be the parallels of the contractile sarcomeric units of muscles , each stress fiber arises from higher-order organization of >10-30 F-actin filaments, numerous actin cross-linking proteins, and non-muscle myosin II molecular motors [1,2,3,4,5,6,7,8]. The stress fibers generate substantial mechanical forces that power cellular contraction against the extracellular matrix [9,10,11,12]. Meanwhile, the formation, stability, dynamics, and morphology of the stress fibers are highly regulated by mechanical and biochemical cues [13,14,15,16,17,18,19,20,21,22]. For instance, upregulation of contractility, actin polymerization, and matrix adhesion promote the formation and thickening of stress fibers, whereas cell relaxation, the inhibition of contractility, and actin cytoskeletal disruption lead to their disassembly and disintegration [15, 17, 18, 20].
A typical adherent cell contains an ensemble of stress fibers that span between adhesion sites or interconnect with one another across the cells, forming an integrated contractility apparatus that plays central roles in morphodynamic programmes such as migration, adhesion, and polarization [21,24,25,, 23–26]. Stress fiber organization therefore underpins important cellular behaviors involved in both normal and pathological processes including developmental morphogenesis and cancer metastasis [27,28,29,30]. Since the stress fibers can be readily visualized both in living or fixed cells, by fluorescence microscopy using F-actin targeting fluorophores [31,32,33], the architecture of the stress fiber network has long been recognized as a key phenotypic reporter of cellular physiology . Nevertheless, although such images may encode valuable information on cellular signaling and mechanobiological states, in many studies the analysis of the stress fiber network architecture were often restricted to qualitative descriptions, in significant part due to the limited availability of appropriate methods for quantitative extraction and analysis of salient features of the stress fiber networks.
Stress fibers typically are observed as networks of numerous elongated filaments. The recognition of filamentous image features has been extensively explored in fields such as geospatial informatics, neurosciences and astrophysics [35,36,37]. However, due to large variations in imaging methods, feature complexity, image resolution, and noise level, methods developed for a given type of curvilinear structures may not be directly applicable to others. For stress fibers, several approaches have previously been developed for their characterization [12,39,40,41,42,43,, 38–44]. For example, order parameters analysis has been used to describe the aggregate image texture and orientation, without explicit treatment of each discrete stress fibers , thus avoiding the challenging task of detecting and segmenting individual filaments. Alternatively, a simulation-based approach can be used to study the stress fiber networks based on Finite Element analysis of idealized cellular architectures [12,41,42,43,44,, 40–45]. Likewise, stress fiber networks can be treated as a micrograph-based linear superposition of filaments such that relevant coefficients can be solved by linear optimization . However, major limitations of these approaches are their inability to extract empirical characteristics such as the dimension, density, and interactions between filaments, which are of key biologically relevance in the study of actin cytoskeletal organization.
We note that an approach that enables the extraction of individual stress fibers is potentially highly beneficial, particularly as this would permit direct correlation between experiments and theoretical models, a key step towards quantitative and predictive understanding of the underlying mechanisms . However, existing methods for discrete filaments extraction have been parametrized for sparsely distributed filaments in cell periphery, cytoskeleton networks polymerized in vitro, or super-resolution microscopy images where improved resolution permits visual distinction of filaments [46,47,48,49,50,51,52,53,54,55]. To our knowledge, a computational method for the identification and complete extraction of the stress fibers in fluorescence micrographs of cells has not been available.
While the stress fibers are often the most prominent F-actin-containing cellular features, a number of technical factors pose significant challenges for their automated extraction. These include the complex organization of the filaments, such as filament intersection and convergence, and the presence of numerous F-actin-containing structures which may appear as bright puncta or as indistinct background intensity which reduces the local contrast of the stress fibers (Fig. 1b, c, Additional file 1: Figure S1). In this study we present a computational strategy to address these challenges, implemented as a package called SFEX (Stress Fibers Extractor) (see Additional files 2 and 3). In brief, this involves two major steps, linear structure enhancement and stress fiber reconstruction (Fig. 1a). In the first, neighborhood-based enhancement methods, line filter transform (LFT) and orientation filter transform (OFT) , were applied to the raw fluorescence image to selectively enhance the contrast of linear objects against different shape profiles, allowing detection and segmentation by binary thresholding. Subsequently, minimal linear filament fragments which represent the centerlines of detected stress fibers were generated from the skeletonized binary images. These were then recombined to reconstruct the traces of individual stress fibers. Carried out iteratively, this process permits extraction of the majority of stress fibers in cells, thus allowing the spatial attributes of both individual fibers and their collective architecture to be determined. Altogether, SFEX enables the automated extraction of the fiber networks from fluorescence micrographs, and thus may facilitate large-scale quantitative analysis of actin stress fiber profiles for dissecting molecular mechanisms or in high-throughput screening applications.
Cell Culture and specimen preparation
Human Osteosarcoma Cells (U2OS) were obtained from the American Type Culture Collection (ATCC, Manassus, VA) and cultured in McCoy’s 5A media (Gibco) supplemented with 10% heat-inactivated fetal bovine serum (Gibco), 1X glutaMAX (Gibco) and 1% penicillin/streptomycin (Gibco) and were maintained in the incubator at 37 °C with 5% CO2. Cells were seeded on a Starter’s CYTOOchipTM in a 35 mm dish at the density of 75,000 cells/ml, and allowed to adhere for 20 mins before incubation. After 30 mins of incubation, unattached cells were removed by triple rinsing with DPBS (Dulbecco’s Phosphate Buffered Saline). The attached cells were allowed to spread for 3 hours in cell culture media before fixation. Cells were fixed with 4% Paraformaldehyde (Electron Microscopy Science) in PBS (Phosphate Buffered Saline), permeabilized with 0.2% Triton X-100 (Sigma) and stained with Alexa Flour 568 Phalloidin (Life Technologies) overnight. The samples were then mounted on a glass slide with PBS as imaging buffer and sealed by vaseline-lanolin-paraffin mixture  for TIRF imaging.
Total Internal Reflection Fluorescence Microscopy (TIRFM) imaging
The specimens were imaged by Nikon Eclipse Ti-E inverted microscope with motorized total internal reflection fluorescence (TIRF) illuminator. The microscope is equipped with a sCMOS camera (Orca Flash 4.0, Hamamatsu) and a 405/488/561/647 TIRF Laser Dichroic filter (Chroma Technologies). Single cells were acquired under TIRF mode with a 60X oil-immersion objective (NA 1.49 Apo TIRF). Fluorophores were excited at 30% intensity of a 60 mW 561 nm laser.
Image enhancement by line and orientation filter transform
For LFT, at each (x, y) pixel we defined a neighborhood of radius r, within which linear features are to be assessed (Fig. 2b). A line segment of length 2r centered at each pixel is rotated stepwise with an angle θ between -90° and 90° (Fig. 2b). The direction along which the accumulated image intensity is the largest is designated the preferred orientation, θ max . As defined by Eq. 1 and 2, this process is repeated for all pixels to generate two image maps: the intensity map (L intensity ), where each entry is the mean pixel value along the preferential direction of that pixel, and the orientation map (L orientation ), which contains the preferred direction at each pixel.
The LFT step above serves to enhance linear features using only the image intensity information. However, in addition to the intensity, linear structures can also be recognized by considering the preferred directions of both the base pixel and its neighbors. To incorporate this information, for each pixel we performed a second filter transform, OFT, to explore whether the neighboring pixels along θ max have similar preferential directions. As calculated by Eqs. 3–8, these criteria thus assigned the probability score for a pixel being on a filamentous structure.
Reconstruction of the stress fiber traces
The significantly increased contrast (Fig. 2f, g) due to LFT and OFT facilitates the use of a segmentation threshold to extract the stress fiber networks. We used Otsu’s method  to determine the initial threshold level. Subsequently the binarized image was skeletonized to determine the centerlines of the filaments. However, we noted that a typical stress fiber is rarely an isolated linear structure, but instead is usually associated with numerous small fibrils branching off to the sides (Fig. 1b, c, red, green and blue arrows). This characteristic gives rise to highly branched skeletons, especially in the arc regions, which impedes simple identification of the centerlines of actin bundles.
To facilitate the tracing of the appropriate stress fiber centerlines, we partitioned the skeleton into separate unbranched linear segments by removing the junction regions whose local 8-connected neighborhood contains more than three filament pixels. This results in a pool of piecewise-linear ‘filament fragments’. Subsequently, we applied a series of geometric constraints to group together the fragments that best capture the stress fibers traces in the original image. We then calculated the propagation direction for each terminus of the fragments, defined as the orientation pointing from the center of mass of the fragment to the tip itself (Fig. 3a). Next, a fan-shaped sector region is generated, anchored at each tip (Fig. 3b). A local search is then performed within each search fan to locate potential tips for fragment connection. Finally, the reconstruction of each stress fiber is carried out by combining every pair of tips, i and j, that satisfies three geometric constraints: (1) similarity, (2) proximity, (3) continuity, as defined by Eqs. 9, 10, 11 (Fig. 3c–e).
with the definition of variables listed in Table 1.
However, in many occasions, one or more bridging fragments that themselves satisfying these three criteria can be found in the search fan (Fig. 3f). To address this, the fourth criterion is applied to check for the existence of such inner fragments and whether they are eligible to bridge the two fragments being investigated, defined as follows (variables definition in Table 1):
For each iteration, it is possible that multiple fragment termini may satisfy the aforementioned four conditions with respect to the base terminus. To select the optimal partner terminus for combination, we therefore introduced a scoring system to compute the priority of all termini based on how they satisfy the similarity and continuity criteria as shown in Eq. 16 (variables definition in Table 1).
The termini pair with the lowest score are then assigned to the same stress fibers, while the unpaired fragment termini are designated as the termini of stress fibers. The process is iterated for all fragments and, when completed, yields the most probable centerline trace of each stress fiber within the networks. Using a desktop workstation, the entire calculation for a typical image (512 × 512 pixels or 55 × 55 μm) was completed in less than two minutes. Subsequently, the centerline trace for each stress fiber can be used for further analysis of the stress fibers in lieu of the entire image, thereby helping to reduce the data set dimension. For example, the orientation of the stress fibers can be calculated using the local information centered around each pixel along the centerline.
Generation of synthetic images and quantification of extraction accuracy
To systematically evaluate the capability of our method in extracting stress fiber filaments, sensitivity analyses were performed to benchmark the dependence of our method on the values of key parameters and the magnitude of background noise. As shown in Fig. 4a, we generated synthetic images using simple geometric patterns of curves and straight lines as the ground-truth, approximating the two types of stress fibers shapes commonly observed in cells. The ground truth skeletons were convoluted with a Gaussian filter with a width comparable to the resolving power of our TIRF microscope (σ = 72 nm, image pixel size 108 nm), and corrupted by Gaussian noise added at varying Peak Signal-to-Noise Ratio (PSNR) (Fig. 4b, Additional file 4: Figure S2). To quantify the performance of our method in detecting individual fibers (Fig. 4c), we determined the false positives and false negatives ratios, defined as 1-(M/D) and 1-(M/G), respectively, where D is the number of computer-identified filaments, G is the number of ground truth filaments and M denotes the number of matches between D and G .
Automated determination of stress fiber widths
The widths of the stress fibers are known to be strongly correlated with both cell contractility or activities of actin regulatory pathways [1,61,, 23, 60–62], and thus can serve as a useful read-out of cell mechanical properties. Although this can be calculated interactively from line profile of stress fibers, due to the large numbers of stress fibers per cells, we sought to automate this process. For each thick stress fiber trace (Additional file 5: Figure S3A, synthetic image example shown), we first calculated the Euclidean distance map (Figure S3D, synthetic image example shown). Since each stress fiber can be considered as an open curve, each distance level generated by connecting pixels with the same distance value appears as a closed loop enveloping the stress fiber. Based on the original fluorescence micrographs, we then calculated the mean of the gradients along each distance level. The mean gradients as a function of the distances from the stress fibers were then obtained, and the distance level with the highest mean image gradient is identified, with the average width of the stress fiber defined as twice this value (Additional file 5: Figure S3G-J).
Extracting information from secondary actin stress fiber networks
The actin stress fibers exhibit significant variations in their sizes and intensity. Thus, while the prominent primary stress fibers can be easily detected, lower intensity secondary fibers that often form complex networks are more challenging to segment accurately (Fig. 5a-c). Nevertheless, as these types of stress fibers can comprise a significant portion of F-actin structures in cells, we also implemented a method to glean quantitative information from such networks. Following the extraction of the high-intensity stress fiber networks from the cell images, Otsu thresholding is used to segment cell areas devoid of thick stress fibers (Fig. 5a-c). The integrated intensity of pixels within this region can then be calculated, for example, as a function of the distance from the cell edge (Fig. 5d, black region). This provides a useful metric for the density distribution of the actin networks as a function of the cell morphology, particularly in case of cells on micropattern, described further below, where the cell edge can be used as a common spatial frame of reference.
Statistical analysis and computation
The computational routines were programmed in MATLAB (Release R2015a, Natick, MA). All computations were performed on Windows 7 workstation (Intel(R) Xeon(R) CPU E5-2640 v3 @ 2.60GHz; RAM 192GB; 64-bit OS). Statistics and graphing were performed in MATLAB.
Results and Discussion
Fluorescence imaging and analysis of actin stress fibers
The spatial distribution of the actin stress fibers is closely coupled to the morphology of the cells. However, with conventional cell culture methods, cell shapes for a given population become highly heterogeneous, complicating a systematic study of the spatial organization of actin structures. One of the methods to regularize cell morphology is the use of micropatterns to constrain cell adhesion to within areas printed with extracellular matrix proteins such as fibronectin (insets 1-3 in Fig. 1b), with the outlying areas coated with non-adhesive materials [31, 63, 64]. As seen in Fig. 1b, human osteosarcoma (U2OS) cells cultured on such fibronectin-micropatterned cover glasses (CYTOOChips, Cytoo Inc.) were highly restricted to geometric forms such as Y-, crossbow-, and disc-shaped patterns, exhibiting significant uniformity across the population (Fig. 1b, Additional file 1: Figure S1). Using Alexa Fluor 568 phalloidin to label F-actin, diffraction-limited fluorescence images of stress fiber organization can be obtained by total internal reflection fluorescence microscopy (TIRFM) imaging.
As seen in Fig. 1b-c and Additional file 1: Figure S1, while the stress fibers are often the most prominent F-actin containing features, fluorescence micrographs of F-actin usually contain numerous other image features that pose significant challenges for the automated extraction of stress fibers. These include the complex organization of the filaments, such as filament intersection and convergence (Fig. 1b, c, Additional file 1: Figure S1), numerous F-actin-containing cellular structures which may appear as bright puncta (Fig. 1b, c, Additional file 1: Figure S1), and nebulous background intensity which reduces the local contrast of the stress fibers (Fig. 1b, c, Additional file 1: Figure S1). Thus, while the qualitative patterns of stress fiber organization can be readily recognized by visual inspection, the automated extraction of the fibers against feature clutters and high background noises remains a difficult task. To address these, we made use of anisotropic image enhancement methods to accentuate the fibrous structure of interest. Anisotropic image enhancement has been widely used for recognizing cellular structures such as cytoskeleton and membrane, with prior knowledge about the shapes of interest [56,66,, 65–67]. We found that robust anisotropic enhancement can be achieved using the LFT and OFT methods which are relatively easier to implement compared to other approaches. LFT and OFT electively highlight all component pixels along filamentous structures, while suppressing non-linear features . As shown in Fig. 2, LFT and OFT enhance the relative contrast of filamentous features and also suppress the intensity of non-filamentous structures, therefore enabling the use of a simple binarization threshold to extract features of interest. The binarized image of the stress fiber networks is then skeletonized and regions of filament junctions are removed to generate unbranched linear fragments. Subsequently, the filaments are reconstructed by using four geometric criteria for filament fragments recombination, as illustrated in Fig. 3. These computational steps, and subsequent quantitative analysis of the networks, can be performed using the included software package, SFEX (Stress Fibers Extractor), described in the Supplementary Information.
Assessing the performance of stress fiber reconstruction
Since the architectures of stress fiber networks differ drastically between different cells even under the same condition, the optimal parameter set for fiber extraction also varies from cell to cell. To aid users in estimating the appropriate input parameter range, particularly against image noise, we performed sensitivity analyses by systematically varying a given pair of key parameters, keeping the rest fixed, and used SFEX to extract the stress fibers from synthetic images as shown in Fig. 4a–c. The results of the analyses were then scored against the ground-truth and the errors were visualized as the heat maps shown in Fig. 4d–i.
We first evaluated how the neighborhood radius r of the LFT/OFT enhancement together with the noise level of the original images affected the accuracy of stress fiber detection. As shown in Fig. 4d and g, both the neighborhood radius r and image quality are strong determinants of the tracing accuracy. The highest accuracy (low false positives and false negatives ratio, blue regions in Fig. 4d and g) was found to correspond to r = ~1.2 μm, which represents the optimal balance, as the large values of r may introduce distortion to curved fibers and hence a low tracing accuracy, while the smaller values of r may not provide sufficient enhancement for the fibers relative to noises. Also, these results suggest that our method appear to perform reliably for image quality exceeding PSNR ~ 20, with rapid deterioration with the increase in image noise (see Additional file 4: Figure S2 for examples of noise levels).
We next explored how the search radius parameter for the fragment reconstruction affected fiber reconstruction accuracy. Surprisingly, we observed that beyond the lower limit of search radius ~16 pixels which corresponds to ~1.7 μm, a relatively wide range of search radius can be used to obtain high reconstruction accuracy, as seen in Fig. 4e and h. This permissiveness suggests that the detection accuracy is not significantly affected as long as the search region encloses enough potential fragment termini for recombination analysis. Also this indicates that our algorithm for termini pairing is probably sufficiently robust that the optimal pairing can be determined regardless of the size of the search fan. Likewise, the performance of our method is relatively resilient with respect to the search fan angle beyond ~20 degree, as shown in Fig. 4f and i, suggesting that a relatively large search sector should be used to avoid potential mismatch of termini pair. Similar to the case of the OFT radius above, the accuracy degrades rapidly once image quality is reduced below PSNR ~20. Here, the filament extraction error likely originates from noise-induced error in the filament propagation direction (Fig. 3b). Interestingly, we note that a narrow region of relatively good accuracy can be observed at ~ PSNR 12-14. As can be seen from Additional file 4: Figure S2, the images are significantly degraded by noise in that range, such that the filaments are extracted in multiple small segments. In this case, the bridging filament criterion (Fig. 3f) would be invoked, which may help improve the accuracy to a certain extent. However, as the zone where this effect takes place is relatively narrow, it is advisable to obtain image quality exceeding PSNR ~ 20 for accurate analysis of the stress fibers.
We also performed sensitivity analysis for the stress fiber width determination. For this, we used a ground truth image consisting of a curved fiber (Additional file 5: Figure S3A and Additional file 6: Figure S4A). The synthetic images were corrupted by noise (Additional file 6: Figure S4B-G) and the width of the filaments calculated using the distance map-based method described above. As shown in Additional file 5: Figure S3B-C, E-F, the width of the fiber obtained from the analysis appears to capture well the width of the fiber across a wide range of noise level, (Additional file 6: Figure S4H), suggesting that our distance map gradient method is highly resistant to image noise.
Quantitative analysis of geometrically confined cells
In recent years, image-based morphometric profiling of cells has been an active area of investigation both in basic cell biological studies and in theranostic developments [68, 69]. In particular, the invention of extracellular matrix micropatterning technologies which standardize cell shapes [31, 63, 64] has greatly facilitated image processing and statistical analysis. Fluorescence micrographs of numerous cells with similar geometrically well-defined shapes can be rapidly acquired by automated microscopes, and ensemble-averaged to yield maps of how specific molecules are stereotypically organized [63, 70, 71]. However, the effect of shape confinement on actin stress fibers organization has not been systematically investigated, despite their dramatic differences between cell shapes (Fig. 6a-f, Additional file 1: Figure S1), in part due to the difficulties of extracting stress fibers from the images. Thus in this study we used our method to explore how the micropatterns may influence actin stress fibers organization.
For this work, we made use of three different micropatterns which are designed to elicit different cell morphological responses, but with comparable total cell areas: Y-shaped, which promotes strong adhesion and prominent stress fibers at cell edges spanning between each vertex (Fig. 1b, left); disc-shaped, which keeps the cells unpolarized, thus devoid of thick stress fibers (Fig. 1b, right), featuring instead circumferential network of secondary stress fibers that connects with adhesion sites near the cell boundary at semi-regular intervals; and crossbow-shaped, which induces clear front-back polarization, with prominent stress fibers connecting the vertices, converging on the cell rear, and intricate networks of secondary actin stress fibers in the frontal arc (Fig. 1b, middle).
To overlay actin images of cells on a given micropattern, we segmented the cells and determined the centers of mass (COMs) of the cell region (Additional file 7: Figure S5A-C), and then performed translation and rotation to align their COMs (Additional file 7: Figure S5D-F). The ensemble image is then calculated from the mean of the aligned images. As shown in Fig. 6a, b, d, e, we observed that the prominent primary stress fibers appear to be nearly uniform between different cells on similar patterns, while the organization of secondary actin stress fibers, especially toward the cell front in polarized cell, is much more complicated. Thus, although micropatterns help standardize the overall cell morphology, significant cell-to-cell variation in discrete structures such as actin stress fibers are still observed and thus spatial information on stress fibers organization is significantly obscured by simple ensemble-averaging analysis.
Architecture of actin stress fiber networks as a function of cell shapes
Since the architecture of the stress fiber networks can be analyzed on a cell-to-cell basis using SFEX, we next demonstrated the utility of our approach in quantifying actin stress fibers in cells on various micropatterns. We performed stress fiber analysis on these cells and outlined below some important metrics useful for differentiating between their distinct actin stress fiber organization phenotypes.
Stress fiber networks, especially in protrusive regions, play important roles in cell polarity and migration. The mesh of secondary actin stress fibers forms ‘band’- like lamella structure [21, 72] parallel to the protrusive edges of the cell (Fig. 6b, c, e, f), which are believed to be the zones of contractility generation by myosin II motors. We therefore first investigated whether the size of the lamella zone is dependent on micropattern-induced cell polarity. We compared actin organization in unpolarized disc pattern and polarized crossbow pattern, plotting the density of F-actin (coverage ratio of F-actin versus total cell area) as a function of the distance from cell edge, using the curved edge of the cell as the spatial frame of reference. As shown in Fig. 6k-l, F-actin densities peak at distance ~1 μm behind the cell edge, then monotonously decrease toward the cell interior. A similar pattern is observed for both unpolarized (disc) and polarized (crossbow), suggesting that the overall radial organization of the actin cytoskeleton may be comparable, regardless of the extent of cell polarization.
We next investigated the orientation distribution of individual secondary actin stress fibers, as these may be another useful indicator of network organization. We skeletonized the network and removed the junction regions to generate linearized fragments and calculated the orientation of F-actin for each skeleton pixel using a 15-by-15-pixel window (Fig. 5e). Our analysis revealed a non-uniform arrangement of actin orientations in polarized cells (cross-bow) as shown in Fig. 6i relative to non-polarized (disc) cells (Fig. 6j, p = 8.4 x 10-5, Kolmogorov-Smirnov test). Such difference in the azimuthal organization of the actin networks likely reflects the shape constraint imposed by the cell edge, since such secondary actin stress fibers are believed to arise from the lamellipodial front (i.e. the arc) , and thus may originate uniformly in disc-shaped cells, but with orientation preference in crossbow-shaped cells. Thus, while the actin networks exhibit comparable radial organization, polarity appears to manifest in terms of azimuthal organizational bias. Future studies that systematically vary the dimension or the symmetry of the pattern may thus be useful in studying how mechanical constraints are sensed and propagated into actin organization differences and vice versa.
Finally, we note that while the width of the stress fibers are generally considered to be mechanosensitive, dependent on factors such as substrate stiffness as well as myosin-II dependent contractility [1,61,, 23, 60–62], such relationship is thus far considered to be largely qualitative [23, 62], and a predictive relationship between stress fiber attributes and cell mechanics has not been established. Towards addressing this question, we explored whether the dimensions of the thick stress fibers are dependent on the size of non-adhesive regions in micropatterns, that is, their lengths. As shown in Fig. 6a, b, actin filament forms thick bundles spanning the entire borders of the non-adhesive regions. The average total length of stress fibers on Y-shaped pattern is ~1.6 times longer than those on crossbow-shaped pattern (Fig. 6g, p < 2x10-21, t-test). Interestingly, the width of the stress fibers on the latter is also ~1.6 times larger than on the former (Fig. 6h, p < 5x10-8, T test) indicative of an inverse proportional relationship between the width and length of stress fibers, which may be interesting to investigate further by variation of the pattern size and in conjunction with pharmacological perturbations.
We proposed an integrated computational strategy for the enhancement, detection, and analysis of actin stress fibers from fluorescence images of F-actin in adherent tissue culture cells. The associated software package, SFEX, is an open-source and highly customizable tool equipped with graphical user interface. We implemented anisotropic neighborhood-based filtering to enhance linear features corresponding to the stress fibers relative to other undesirable structures. This enables efficient segmentation and skeletonization of stress fibers to yield the centerline traces. Using a set of geometric criteria, the stress fiber skeletons were regrouped from piecewise-linear fragments to reconstruct the entire networks. We evaluated the robustness of our approach and defined the parameter range for optimal performance. Using cells cultured on well-defined micropattern substrates, we analyzed their stress fiber organization to highlight their dependence on cell shapes. Our approach can be performed with minimal user input and thus should be amenable for processing large datasets. This could be useful both for time-lapse movies where the analysis of large number of frames by self-consistent criteria are required, or for high-throughput applications whereby cells are subject to external mechanical or biochemical perturbations.
Total Internal Reflection Fluorescence Microscope
Line Filter Transform
Orientation Filter Transform
Stress Fiber Extractor
Human Osteosarcoma Cells
Dulbecco’s Phosphate Buffered Saline
Phosphate Buffered Saline
Peak Signal-to-Noise Ratio
center(s) of mass
Burridge K, Wittchen ES. The tension mounts: stress fibers as force-generating mechanotransducers. J Cell Biol. 2013;200(1):9–19.
Cramer LP, Siebert M, Mitchison TJ. Identification of novel graded polarity actin filament bundles in locomoting heart fibroblasts: implications for the generation of motile force. J Cell Biol. 1997;136(6):1287–305.
Abercrombie M, Heaysman JE, Pegrum SM. The locomotion of fibroblasts in culture. IV. Electron microscopy of the leading lamella. Exp Cell Res. 1971;67(2):359–67.
Buckley IK, Porter KR. Cytoplasmic fibrils in living cultured cells. A light and electron microscope study. Protoplasma. 1967;64(4):349–80.
Weber K, Groeschel-Stewart U. Antibody to myosin: the specific visualization of myosin-containing filaments in nonmuscle cells. Proc Natl Acad Sci U S A. 1974;71(11):4561–4.
Lazarides E, Burridge K. Alpha-actinin: immunofluorescent localization of a muscle structural protein in nonmuscle cells. Cell. 1975;6(3):289–98.
Adams JC. Formation of stable microspikes containing actin and the 55 kDa actin bundling protein, fascin, is a consequence of cell adhesion to thrombospondin-1: implications for the anti-adhesive activities of thrombospondin-1. J Cell Sci. 1995;108(Pt 5):1977–90.
Wang K, Ash JF, Singer SJ. Filamin, a new high-molecular-weight protein found in smooth muscle and non-muscle cells. Proc Natl Acad Sci U S A. 1975;72(11):4483–6.
Sawada Y, Tamada M, Dubin-Thaler BJ, Cherniavskaya O, Sakai R, Tanaka S, Sheetz MP. Force sensing by mechanical extension of the Src family kinase substrate p130Cas. Cell. 2006;127(5):1015–26.
Kanchanawong P, Shtengel G, Pasapera AM, Ramko EB, Davidson MW, Hess HF, Waterman CM. Nanoscale architecture of integrin-based cell adhesions. Nature. 2010;468(7323):580–84.
Burridge K, Guilluy C. Focal adhesions, stress fibers and mechanical tension. Exp Cell Res. 2016;343(1):14–20.
Walcott S, Sun SX. A mechanical model of actin stress fiber formation and substrate elasticity sensing in adherent cells. Proc Natl Acad Sci U S A. 2010;107(17):7757–62.
Ridley AJ, Hall A. The small GTP-binding protein rho regulates the assembly of focal adhesions and actin stress fibers in response to growth factors. Cell. 1992;70(3):389–99.
Maekawa M, Ishizaki T, Boku S, Watanabe N, Fujita A, Iwamatsu A, Obinata T, Ohashi K, Mizuno K, Narumiya S. Signaling from Rho to the actin cytoskeleton through protein kinases ROCK and LIM-kinase. Science. 1999;285(5429):895–8.
Pollard TD. Regulation of actin filament assembly by Arp2/3 complex and formins. Annu Rev Biophys Biomol Struct. 2007;36:451–77.
Bugyi B, Carlier MF. Control of actin filament treadmilling in cell motility. Annu Rev Biophys. 2010;39:449–70.
Bershadsky AD, Ballestrem C, Carramusa L, Zilberman Y, Gilquin B, Khochbin S, Alexandrova AY, Verkhovsky AB, Shemesh T, Kozlov MM. Assembly and mechanosensory function of focal adhesions: experiments and models. Eur J Cell Biol. 2006;85(3-4):165–73.
Chrzanowska-Wodnicka M, Burridge K. Rho-stimulated contractility drives the formation of stress fibers and focal adhesions. J Cell Biol. 1996;133(6):1403–15.
Smith MA, Blankman E, Gardel ML, Luettjohann L, Waterman CM, Beckerle MC. A zyxin-mediated mechanism for actin stress fiber maintenance and repair. Dev Cell. 2010;19(3):365–76.
Discher DE, Janmey P, Wang YL. Tissue cells feel and respond to the stiffness of their substrate. Science. 2005;310(5751):1139–43.
Geiger B, Spatz JP, Bershadsky AD. Environmental sensing through focal adhesions. Nat Rev Mol Cell Biol. 2009;10(1):21–33.
Riveline D, Zamir E, Balaban NQ, Schwarz US, Ishizaki T, Narumiya S, Kam Z, Geiger B, Bershadsky AD. Focal contacts as mechanosensors: externally applied local mechanical force induces growth of focal contacts by an mDia1-dependent and ROCK-independent mechanism. J Cell Biol. 2001;153(6):1175–86.
Hotulainen P, Lappalainen P. Stress fibers are generated by two distinct actin assembly mechanisms in motile cells. J Cell Biol. 2006;173(3):383–94.
Chen WT. Mechanism of retraction of the trailing edge during fibroblast movement. J Cell Biol. 1981;90(1):187–200.
Parsons JT, Horwitz AR, Schwartz MA. Cell adhesion: integrating cytoskeletal dynamics and cellular tension. Nat Rev Mol Cell Biol. 2010;11(9):633–43.
Prager-Khoutorsky M, Lichtenstein A, Krishnan R, Rajendran K, Mayo A, Kam Z, Geiger B, Bershadsky AD. Fibroblast polarization is a matrix-rigidity-dependent process controlled by focal adhesion mechanosensing. Nat Cell Biol. 2011;13(12):1457–65.
Gilbert SF. Developmental Biology. 7th ed. Sunderland, MA. Sinauer Associates, Inc.; 2003.
Ridley AJ, Schwartz MA, Burridge K, Firtel RA, Ginsberg MH, Borisy G, Parsons JT, Horwitz AR. Cell migration: integrating signals from front to back. Science. 2003;302(5651):1704–9.
Tojkander S, Gateva G, Lappalainen P. Actin stress fibers--assembly, dynamics and biological roles. Nat Cell Biol. 2012;125(Pt 8):1855–64.
Pollard TD, Cooper JA. Actin, a central player in cell shape and movement. Science. 2009;326(5957):1208–12.
Tee YH, Shemesh T, Thiagarajan V, Hariadi RF, Anderson KL, Page C, Volkmann N, Hanein D, Sivaramakrishnan S, Kozlov MM, et al. Cellular chirality arising from the self-organization of the actin cytoskeleton. Nat Cell Biol. 2015;17(4):445–57.
Vallenius T. Actin stress fibre subtypes in mesenchymal-migrating cells. Open Bio. 2013;3(6):130001.
Skau CT, Plotnikov SV, Doyle AD, Waterman CM. Inverted formin 2 in focal adhesions promotes dorsal stress fiber and fibrillar adhesion formation to drive extracellular matrix assembly. Proc Natl Acad Sci U S A. 2015;112(19):E2447–2456.
Blanchoin L, Boujemaa-Paterski R, Sykes C, Plastino J. Actin dynamics, architecture, and mechanics in cell motility. Physiol Rev. 2014;94(1):235–63.
Breuer D, Nikoloski Z. DeFiNe: an optimisation-based method for robust disentangling of filamentous networks. Sci Rep. 2015;5:18267.
Helmstaedter M. Cellular-resolution connectomics: challenges of dense neural circuit reconstruction. Nat Methods. 2013;10(6):501–7.
De J, Cheng L, Zhang X, Lin F, Li H, Ong KH, Yu W, Yu Y, Ahmed S. A graph-theoretical approach for tracing filamentary structures in neuronal and retinal images. IEEE Trans Med Imaging. 2016;35(1):257–72.
Gupta M, Sarangi BR, Deschamps J, Nematbakhsh Y, Callan-Jones A, Margadant F, Mege RM, Lim CT, Voituriez R, Ladoux B. Adaptive rheology and ordering of cell cytoskeleton govern matrix rigidity sensing. Nat Commun. 2015;6:7525.
Rigort A, Gunther D, Hegerl R, Baum D, Weber B, Prohaska S, Medalia O, Baumeister W, Hege HC. Automated segmentation of electron tomograms for a quantitative description of actin filament networks. J Struct Biol. 2012;177(1):135–44.
Wood ST, Dean BC, Dean D. A linear programming approach to reconstructing subcellular structures from confocal images for automated generation of representative 3D cellular models. Med Image Anal. 2013;17(3):337–47.
Ingber DE. Tensegrity I. Cell structure and hierarchical systems biology. J Cell Sci. 2003;116(Pt 7):1157–73.
Unnikrishnan GU, Unnikrishnan VU, Reddy JN. Constitutive material modeling of cell: a micromechanics approach. J Biomech Eng. 2007;129(3):315–23.
Peeters EA, Oomens CW, Bouten CV, Bader DL, Baaijens FP. Mechanical and failure properties of single attached cells under compression. J Biomech. 2005;38(8):1685–93.
Karcher H, Lammerding J, Huang H, Lee RT, Kamm RD, Kaazempur-Mofrad MR. A three-dimensional viscoelastic model for cell deformation with experimental verification. Biophys J. 2003;85(5):3336–49.
Gouget CL, Hwang Y, Barakat AI. Model of cellular mechanotransduction via actin stress fibers. Biomech Model Mechanobiol. 2016;15(2):331–44.
Yang Q, Karpikov A, Toomre D, Duncan JS. 3-D reconstruction of microtubules from multi-angle total internal reflection fluorescence microscopy using Bayesian framework. IEEE Trans Image Process. 2011;20(8):2248–59.
Li HS, Shen T, Huang XL. Actin filament segmentation using dynamic programming. Lect Notes Comput Sc. 2011;6801:411–23.
Li HS, Shen TA, Vavylonis D, Huang XL. Actin filament segmentation using spatiotemporal active-surface and active-contour models. Med Image Comput Comput Assist Interv. 2010;6361(Pt I):86–94.
Li HS, Shen T, Smith MB, Fujiwara I, Vavylonis D, Huang XL. Automated actin filament segmentation, tracking and tip elongation measurements based on open active contour models. 2009 Ieee International Symposium on Biomedical Imaging: From Nano to Macro, Vols 1 and 2. 2009. p. 1302–5.
Li HS, Shen T, Vavylonis D, Huang XL. Actin filament tracking based on particle filters and stretching open active contour models. Medical Image Computing and Computer-Assisted Intervention - Miccai 2009, Pt Ii, Proceedings. 2009;5762:673–81.
Xu T, Vavylonis D, Huang XL. 3D actin network centerline extraction with multiple active contours. Med Image Anal. 2014;18(2):272–84.
Sargin ME, Altinok A, Kiris E, Feinstein SC, Wilson L, Rose K, Manjunath BS: Tracing microtubules in live cell images. 2007 4th Ieee International Symposium on Biomedical Imaging : Macro to Nano, Vols 1-3 2007:296-299
Smith MB, Li H, Shen T, Huang X, Yusuf E, Vavylonis D. Segmentation and tracking of cytoskeletal filaments using open active contours. Cytoskeleton. 2010;67(11):693–705.
Xu T, Li H, Shen T, Ojkic N, Vavylonis D, Huang X. Extraction and analysis of actin networks based on open active contour models. InBiomedical Imaging. From Nano to Macro, 2011 IEEE International Symposium. 2011. p. 1334–40. IEEE.
Zhang Z, Nishimura Y, Kanchanawong P. Extracting microtubule networks from superresolution single-molecule localization microscopy data. Mol Biol Cell. 2017;28(2):333–45.
Sandberg K, Brega M. Segmentation of thin structures in electron micrographs using orientation fields. J Struct Biol. 2007;157(2):403–15.
Shin WD, Fischer RS, Kanchanawong P, Kim Y, Lim J, Meyers KA, Nishimura Y, Plotnikov SV, Thievessen I, Yarar D, et al. A Versatile, Multicolor Total Internal Reflection Fluorescence and Spinning Disk Confocal Microscope System for High-Resolution Live Cell Imaging. In: Goldman RD, Swedlow JR, Spector DL, editors. Live Cell Imaging: A Laboratory Manual. 2nd ed. Cold Spring Harbor, New York: Cold Spring Harbor Laboratory Press; 2010.
Otsu N. A threshold selection method from gray-level histograms. IEEE Trans Syst Man Cybern. 1979;9(1):62–6.
Matov A, Applegate K, Kumar P, Thoma C, Krek W, Danuser G, Wittmann T. Analysis of microtubule dynamic instability using a plus-end growth marker. Nat Methods. 2010;7(9):761–8.
Maraldi M, Valero C, Garikipati K. A computational study of stress fiber-focal adhesion dynamics governing cell contractility. Biophys J. 2014;106(9):1890–901.
Peterson LJ, Rajfur Z, Maddox AS, Freel CD, Chen Y, Edlund M, Otey C, Burridge K. Simultaneous stretching and contraction of stress fibers in vivo. Mol Biol Cell. 2004;15(7):3497–508.
Stachowiak MR, Smith MA, Blankman E, Chapin LM, Balcioglu HE, Wang S, Beckerle MC, O’Shaughnessy B. A mechanical-biochemical feedback loop regulates remodeling in the actin cytoskeleton. Proc Natl Acad Sci U S A. 2014;111(49):17528–33.
Thery M, Racine V, Piel M, Pepin A, Dimitrov A, Chen Y, Sibarita JB, Bornens M. Anisotropy of cell adhesive microenvironment governs cell internal organization and orientation of polarity. Proc Natl Acad Sci U S A. 2006;103(52):19771–6.
Wei SC, Fattet L, Tsai JH, Guo Y, Pai VH, Majeski HE, Chen AC, Sah RL, Taylor SS, Engler AJ, et al. Matrix stiffness drives epithelial-mesenchymal transition and tumour metastasis through a TWIST1-G3BP2 mechanotransduction pathway. Nat Cell Biol. 2015;17(5):678–88.
Jiang M, Ji Q, McEwen BF. Enhancement of microtubules in em tomography. InBiomedical Imaging. Nano to Macro, 2004. IEEE International Symposium. 2004. p. 1123–26. IEEE.
Neumann MH, von Sachs R. Wavelet thresholding in anisotropic function classes and application to adaptive estimation of evolutionary spectra. 1997:38–76.
Mosaliganti KR, Noche RR, Xiong F, Swinburne IA, Megason SG. ACME: automated cell morphology extractor for comprehensive reconstruction of cell membranes. PLoS Comput Biol. 2012;8(12):e1002780.
Nikitina L, Ahammer H, Blaschitz A, Gismondi A, Glasner A, Schimek MG, Dohr G, Sedlmayr P. A new method for morphometric analysis of tissue distribution of mobile cells in relation to immobile tissue structures. PLoS One. 2011;6(3):e15086.
Yu H, Lim KP, Xiong S, Tan LP, Shim W. Functional morphometric analysis in cellular behaviors: shape and size matter. Adv Healthc Mater. 2013;2(9):1188–97.
von Erlach TC, Hedegaard MA, Stevens MM. High resolution Raman spectroscopy mapping of stem cell micropatterns. Analyst. 2015;140(6):1798–803.
Tseng Q, Wang I, Duchemin-Pelletier E, Azioune A, Carpi N, Gao J, Filhol O, Piel M, Thery M, Balland M. A new micropatterning method of soft substrates reveals that different tumorigenic signals can promote or reduce cell contraction levels. Lab Chip. 2011;11(13):2231–40.
Ponti A, Machacek M, Gupton SL, Waterman-Storer CM, Danuser G. Two distinct actin networks drive the protrusion of migrating cells. Science. 2004;305(5691):1782.
Burnette DT, Manley S, Sengupta P, Sougrat R, Davidson MW, Kachar B, Lippincott-Schwartz J. A role for actin arcs in the leading-edge advance of migrating cells. Nat Cell Biol. 2011;13(4):371–81.
We thank the staffs of the microscopy, computing, and open lab core facilities in Mechanobiology Institute.
Singapore National Research Foundation: NRF Competitive Research Programme (NRF2012NRF-CRP-001-084); Singapore National Research Foundation: NRF Fellowship (NRF-NRFF-2011-04); Mechanobiology Institute (Seed Funding and Graduate Scholarship); Singapore Ministry of Education (MOE-T2-1-124).
Availability of data and materials
Other datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Z.Z. and P.K. conceived and designed the study. X.S. performed the cell biological and microscopy experiments. Z.Z. developed the software and performed data analysis. Z.Z. and P.K. wrote the manuscript with inputs from all authors. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Actin stress fibers architectures of U2OS cells on micropatterns. Inverse contrast TIRFM images of F-actin in U2OS cells plated on Y- (top row), disk- (middle row) and crossbow-shaped (bottom row) micropatterns. Scale bar, 5 μm (TIF 4029 kb).
SupplementarySoftware_SFEX_BINF-D-16-00942. Stress Fiber Extractor (SFEX) 1.0. This file contains a test image and all the source code associated with the method proposed in this study (DOCX 6571 kb).
SupplementarySoftware_SFEX User-Manual-BINF-D-16-00942. SFEX 1.0 User Manual. This file is a step-by-step introduction about the usage of SFEX (The software in Additional file 2) (ZIP 633 kb).
Synthetic images used for assessing filament reconstruction accuracy. Synthetic ground truth images with introduced noise at different levels (TIF 6847 kb).
Sensitivity analysis of filament width measurement. A) Ground truth image. B, C) 3D visualization of images with low (B) and high (C) noise. Projected 2D images are shown below their 3D view. Detected filament contour is highlighted in black. D) Distance map based on ground truth image (A). Colorbar, distance to the centerline of ground truth filament. E, F) Enlarged view of regions highlighted by blue arrows in (B) and (C). Detected filament contour is highlighted in black. G) Stress fibers in a U2OS cell plated on Y-shaped micropattern. Contours of ventral stress fibers are highlighted in red. Scale bar, 5 μm. H-J) Enlarged views of regions indicated by red arrows in (G). Scale bar, 1 μm (TIF 13861 kb).
Images for sensitivity analysis of filament width measurement. A) Synthetic ground truth image. B-G) Synthetic image with introduced noise at different levels. H) Detected filament length as a function of image noise (TIF 4838 kb).
Alignment of cells on micropatterns. A-C) Overlay of cell contours before alignment for Y, crossbow, and disc patterns, respectively. Crosses denote the center-of-mass (COM) of the cell regions. D-F) COM-based overlay of aligned cell contours. n = 10 for each micropattern (TIF 943 kb).