A Robust Actin Filaments Image Analysis Framework

Abstract

The cytoskeleton is a highly dynamical protein network that plays a central role in numerous cellular physiological processes, and is traditionally divided into three components according to its chemical composition, i.e. actin, tubulin and intermediate filament cytoskeletons. Understanding the cytoskeleton dynamics is of prime importance to unveil mechanisms involved in cell adaptation to any stress type. Fluorescence imaging of cytoskeleton structures allows analyzing the impact of mechanical stimulation in the cytoskeleton, but it also imposes additional challenges in the image processing stage, such as the presence of imaging-related artifacts and heavy blurring introduced by (high-throughput) automated scans. However, although there exists a considerable number of image-based analytical tools to address the image processing and analysis, most of them are unfit to cope with the aforementioned challenges. Filamentous structures in images can be considered as a piecewise composition of quasi-straight segments (at least in some finer or coarser scale). Based on this observation, we propose a three-steps actin filaments extraction methodology: (i) first the input image is decomposed into a 'cartoon' part corresponding to the filament structures in the image, and a noise/texture part, (ii) on the 'cartoon' image, we apply a multi-scale line detector coupled with a (iii) quasi-straight filaments merging algorithm for fiber extraction. The proposed robust actin filaments image analysis framework allows extracting individual filaments in the presence of noise, artifacts and heavy blurring. Moreover, it provides numerous parameters such as filaments orientation, position and length, useful for further analysis. Cell image decomposition is relatively under-exploited in biological images processing, and our study shows the benefits it provides when addressing such tasks. Experimental validation was conducted using publicly available datasets, and in osteoblasts grown in two different conditions: static (control) and fluid shear stress. The proposed methodology exhibited higher sensitivity values and similar accuracy compared to state-of-the-art methods.

Abstract

The cytoskeleton is a highly dynamical protein network that plays a central role in numerous cellular physiological processes, and is traditionally divided into three components according to its chemical composition, i.e. actin, tubulin and intermediate filament cytoskeletons. Understanding the cytoskeleton dynamics is of prime importance to unveil mechanisms involved in cell adaptation to any stress type. Fluorescence imaging of cytoskeleton structures allows analyzing the impact of mechanical stimulation in the cytoskeleton, but it also imposes additional challenges in the image processing stage, such as the presence of imaging-related artifacts and heavy blurring introduced by (high-throughput) automated scans. However, although there exists a considerable number of image-based analytical tools to address the image processing and analysis, most of them are unfit to cope with the aforementioned challenges. Filamentous structures in images can be considered as a piecewise composition of quasi-straight segments (at least in some finer or coarser scale). Based on this observation, we propose a three-steps actin filaments extraction methodology: (i) first the input image is decomposed into a 'cartoon' part corresponding to the filament structures in the image, and a noise/texture part, (ii) on the 'cartoon' image, we apply a multi-scale line detector coupled with a (iii) quasi-straight filaments merging algorithm for fiber extraction. The proposed robust actin filaments image analysis framework allows extracting individual filaments in the presence of noise, artifacts and heavy blurring. Moreover, it provides numerous parameters such as filaments orientation, position and length, useful for further analysis. Cell image decomposition is relatively under-exploited in biological images processing, and our study shows the benefits it provides when addressing such tasks. Experimental validation was conducted using publicly available datasets, and in osteoblasts grown in two different conditions: static (control) and fluid shear stress. The proposed methodology exhibited higher sensitivity values and similar accuracy compared to state-of-the-art methods.

Author Summary

We propose a novel actin filaments cytoskeleton analysis framework that allows extracting quasi-straight individual fibers in a robust manner, and provides their respective position, orientation, and length as output. The proposed framework is defined as a three-steps processing sequence, that can explicitly cope with high-throughput imaging related issues, such as noise/artifacts presence and heavy blurring, and can similarly process artifacts-free and well-focused images.

Citation: Alioscha-Perez M, Benadiba C, Goossens K, Kasas S, Dietler G, Willaert R, et al. (2016) A Robust Actin Filaments Image Analysis Framework. PLoS Comput Biol 12(8): e1005063. doi:10.1371/journal.pcbi.1005063

Editor: Niels Volkmann, Burnham Institute for Medical Research, UNITED STATES

Received: November 19, 2015; Accepted: July 15, 2016; Published: August 23, 2016

Copyright: © 2016 Alioscha-Perez et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Simulated data and results can be

Fig 1. Obtained cytoskeleton image can include some degrees of blurring.

(a) Illustration of the cap and basal actin layers. (b) Image of a rat osteoblast (bone producing cell) with phalloidin (red) stained actin cytoskeleton and dapi stained nucleus (blue).

http://dx.doi.org/10.1371/journal.pcbi.1005063.g001

Even in such cases, one can assume that actin fibers are a piecewise composition of quasi-straight segments (at least in some finer or coarser scale), and the directionality of those segments can serve as guide to track individual (blurred and potentially overlapping) filaments. On the contrary, noise and artifacts present in the image are unlikely to exhibit the fibers specific geometry. Thus, quasi-straight clues present in the image can be associated to fibers, since such structures are unlikely to emerge from a random process (Helmholtz principle [22]).

To separate noise/texture from actin fibers, we propose using image decomposition [23-25]. Given an input image f, such approach decompose it as f = u + v, where the image u represents the 'cartoon' information, and v the 'texture' and/or noise image. Among the existing techniques, the authors of [23] proposed a morphological component analysis (MCA) based image separation which constructs a sparse representation of an image and separates the image into morphological components (MCs). This method assumes the 'cartoon' and 'texture' be represented via some known basis functions. We borrow the image decomposition approach of [23] and propose the curvelets model [26] as basis-function for the quasi-linear fiber content, and a wavelet model for the artifacts (texture/noise) content. Our motivation of using the curvelets transform is due to the fact that common denominators among actin filaments structures (as shown in Fig 1b) are, (i) they display anisotropic line-like features, (ii) they basically show behaviors of C2-continuous curves, and (iii) they are relatively smooth in the direction along the filament. Curvelets enables the possibility to directionally analyze an image with different angular resolutions in a single and effective transform [27]. Curvelets take the form of basis elements, which have elongated effective support; i.e., length > width [28]. This method is a good candidate for the detection of anisotropic structures of different lengths [27], such as actin filaments.

In this work, we propose a novel actin filaments cytoskeleton analysis framework that allows extracting quasi-straight individual fibers in a robust manner, and provides their respective position, orientation, and length as output. The proposed framework can explicitly cope with high-throughput imaging related issues, such as noise/artifacts presence and heavy blurring, and can similarly process artifacts-free and well-focused images. During experiments, the proposed model was able to extract a higher number of individual fibers compared to other state-of-the-art models in several cytoskeleton images.

State of the Art

Computer vision and image-analytical tools are essential in order to study the biology of cells [29]. Several approaches exist for the analysis of filamentous structures [19-21, 30-36], consisting in one or several sequential processing steps: pre-processing, filaments network segmentation, and individual fibers extraction.

The filaments analysis problem can be solved in many different ways, as illustrated in Fig 2a. These solutions can be roughly organized in at least three main categories (middle layers in Fig 2a and 2b): based on filaments directionality, based on filaments network, and based on single-filaments extraction. The processing path (from the image to output) will naturally determine the computational burden, accuracy and the amount of information to be extracted from the image. The strategy we followed is depicted in Fig 2b.

Fig 2. Different analysis strategies will determine the amount of information to be extracted from the image.

(a) Possible fibers analysis strategies. (b) The strategy followed involve three sequential stages: directional filaments filtering and enhancement, filaments network segmentation and individual fibers extraction.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g002

In the first category, the filaments information is directly extracted from the pre-processed image [30, 35], sidestepping the foreground/background segmentation step. Segmentation-related errors were avoided at the cost of limiting the model applicability to directionality-related analysis [37] only. Networks-related approaches extract the filaments networks [19] after filtering, allowing to complement orientation information with motion and filaments density analysis [38]. One category of methods for the identification of fiber networks uses template matching [39], where prior knowledge about the target is incorporated into 3D (or 2D) templates. These template-based detection methods are more selective and impose a similarity function. In [40], actin filaments in cryo-ET data sets have been segmented using a stochastic template-based search, which combines a genetic algorithm and a bidirectional expansion strategy. Template matching methods does not resolve the possible template overlap when tracing converging filament branches, thus network junctions could be left undetected. Such problem has been solved in [41] where an active contour based method allowed increasing the segmentation robustness by incorporating prior information about the filament shape. Such approach explicitly models the linear nature of filaments.

The overall strategy of the network-based approaches is to extract the whole curvilinear network at once. In contrast, our work aims at providing information at the single-fiber level, where individual filament segments [20, 21, 32] are extracted. This allows to perform further analysis [11, 16, 42-46] taking into account fibers position, orientation, and length. However, errors introduced during the different processing stages can accumulate and, unless appropriate computational steps are taken, later analysis (statistical and other) can be compromised.

Directional filtering and edges enhancement.

Decomposing an image into meaningful components is an important and challenging inverse problem in computer vision. Image filtering/denoising is an example of image decomposition, which separates the image into signal parts and noise parts. Laplace filters [47], steerable (directional) Gaussian [48], or a combination of them [21] have been used in the literature. More sophisticated techniques, such as Coherence-enhancement filters [49] are well designed to enhance flow-like structures and merge those with (small) discontinuities in a proper way, but they are unable to tackle blurring satisfactorily. In [44], the Radon transform has been used in the context of filaments detection, as a pre-processing step applied to the input image.

Under the perspective that we are addressing a content-specific computer vision application, i.e. separating filamentous structures from other image content, in this work we consider image decomposition as the process of separating an image in conceptually and theoretically different components [23-25]. In these methods a typical assumption is made that the given image is a linear mixture of several source images of more coherent origin. For instance, a nature image is linearly mixed by three coherent components, 'cartoon' structure image, 'texture' image and 'noise' image. The different components require content-specific modelling and representations, and motivate an analysis in two parallel channels. The structure part for 1D feature detection (edges, ridges, etc.), segmentation, object recognition and shape analysis, while the texture component can be used for surface analysis, etc. In the proposed framework, we first perform a sparse multi-source separation strategy [23, 50, 51] to separate the image content into two main parts (sources): the artifacts related content (considered as texture and/or noise) and the filaments related content. Edge enhancement and filament extraction is then applied on the image containing the filaments. We make use of a curvelets model [27] for the quasi-linear fiber content representation, and a wavelet model for the artifacts content. As indicated above, curvelets are ideal candidates for multi-scale edges and directional representations compared to other models.

Network binarization.

The network segmentation requires separating those pixels belonging to the network from the rest, in a process known as binarization. Standard binarization techniques includes global thresholding, such as Otsu's method [52], and local adaptive thresholding based neighborhood means; a combination of both is often used [21, 53]. Instead, taking into account the multi-scale directional nature of the fibers, we performed a multi-scale line segmentation similar to [34] but combined with adaptive local thresholding. The results provided the actin filament network but not yet the individual fibers.

Fibers extraction.

Individual filaments extraction is based on line segment detectors. Hough-based detection approaches exhibit a high computational burden, while light models [32] can miss line segments whose probability is near the decision boundary. Extraction strategies such as the one proposed in [21] are efficient, but they fail to detect several heavy-blurred filaments. We propose a segments extraction algorithm that detects overlapping segments based on a multi-scale directional line response. Quasi-straight line segments of fixed length are extracted at some specific scale, and iteratively merged in order to obtain the final individual fibers.

Materials and Methods

In this work, we propose a sequential three stage processing framework as illustrated in Fig 3. It involves a multi-source filaments separation for the first stage, a multi-scale line detection for the second stage, and a filaments segments merging algorithm for the last stage. The output of each processing step serve as input to the next step, detailed in the following sections.

Fig 3. The proposed framework is based on a specific three stage sequence of processing steps.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g003

Cells Culture

MC3T3-E1 cells, established as an osteoblastic cell line, were provided from Sigma (99072810). Osteoblasts were grown in alpha-modified minimal essential medium ([alpha] -MEM; Life Technologies) supplemented with 10% fetal bovine serum (Life Technologies), 1% L-glutamin (Life Technologies) and 1% penicillin/streptomycin (Life Technologies) in an incubator at 5% CO2 heated at 37°C. The medium was changed twice weekly, and the cells were subcultured into 75 cm2 culture flasks by detaching them gently after a brief PBS rinsing step followed by Trypsin treatment once the cells were reaching subconfluency. For mechanical stimulation, MC3T3-E1 cells were plated into Ibidi M-slide device (M-Slide I 0.4 Luer, ibiTreat from Ibidi) at a concentration of 1.5 × 104 cells/mL. After overnight culture, the medium was replaced by [alpha] -MEM (powder exempt of Bicarbonate, Life Technologies) 10% FBS 1% L-glutamin 1% penicillin/streptomycin 25 mM HEPES (Life Technologies) adjusted to pH 7.4 for 4 hours in the incubator without CO2 at 37°C.

Fluid Shear Stress Induction

The shear stress was induced by a fluidic system composed of two containers and a pump. This system, depicted in Fig 4a, permits a gravity driven constant flow of culture medium in the chamber containing the osteoblasts.

Fig 4. Fluid shear stress induction.

(a) Schematics of the fluidic system we used to exert a constant shear stress on osteoblasts. 1: Ibidi culture chamber with osteoblasts attached to the lower plate (not depicted). 2: Container located under the culture chamber. 3: Peristaltic pump. 4: Container located above the culture chamber. Its position above the culture chamber permits to adjust the flow rate. 5: Microscope and camera (not shown). (b) Image of the culturing chamber showing the border of the micro-fluidic chip, and the direction of the shear stress flow at [asymptotically =] 80°.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g004

The flow rate was adjusted by setting the container (see Fig 4a nb 4) at the appropriate height above the osteoblasts. In our case, it provided a shear stress between 9 and 12 dyn/cm2. This value corresponds roughly to the one to which endothelial cells are exposed to in our arteries. The direction of the shear stress flow is estimated around [asymptotically =] 80° as illustrated in Fig 4b.

Actin Staining of Osteoblasts

After stress exposure for 4 hours, the cells were rinsed in phosphate saline buffer (PBS) and chemically fixed with 4% paraformaldehyde (PFA) for 15 minutes. Before immunostaining the cells were rinsed twice in PBS. Eventually, the fixed cells were incubated with PBS 0.2% Triton X100 for 20 minutes, and exposed to Alexa Fluor 568 phalloidin (Molecular Probes) for 1 hour at room temperature. Finally the cells were rinsed again with PBS before a quick water rinse and coverage with cover glass (by using polyvinylalcohol (Sigma)). A static osteoblast culture, in the same ibidi device, was performed and actin stained as control.

Osteoblasts Imaging

Imaging was performed with an inverted Axiovert 200M system with a 40x Plan-Neofluor (Carl Zeiss, Oberkochen, Germany). We used a motorized platform (MS-2000, Applied Scientific Instrumentation, with NanoDrive controller; Mad City Labs, Madison, WI) to scan the sample. Acquisition was performed using a CoolSnap HQ2 camera (PhotoMetrics, Tucson, AZ) and the Multi-Dimensional Acquisition module of the software Metamorph (Molecular Devices, Sunnyvale, CA).

Image Decomposition

The problem of separating an image into different semantic constituent parts, 'textures' and 'cartoon', can be addressed by several approaches. Variational calculus [24, 25], and sparse multi-source separation [23] are among the most popular methods and their success varies depending on the image nature. In this work, we follow the Starck's [23] image decomposition which is based on the Basis-Pursuit denoising (BPDN) algorithm. The basic idea behind this algorithm is to choose two appropriate dictionaries, one for the representation of 'textures', and the other for the 'cartoon' parts. Both dictionaries are to be designed such that each leads to sparse representations over the images it is serving, while yielding non-sparse representations on the other content type.

In our problem, we model an imaged (actin) cytoskeleton f as a combination of two sources plus some additive noise: (1) being va the background artifacts related content, and uf the filamentous elements; the noise is represented by [eta]. In [23] the sparse source separation problem, for the image f, has been defined as: (2) with [gamma] > 0, [delta] ≥ 0 and p ∈ {0, 1}. The obtained sparse coefficients [Lambda]f provides the fibers-related content uf = Df [Lambda]f, referred to as fibers image in the rest of the paper, and [Lambda]a provides the artifacts related content ua = Da [Lambda]a, referred to as artifacts image; the norm p ∈ {0, 1} determines the reguralizer type and the parameter [gamma] regulates the solution coefficients. The reconstructed image is an approximation of f involving artifacts-related and fiber-related dictionaries, Da and Df respectively; the reminder is usually related to noise [eta].

The definition of the dictionaries Da and Df is very much related to the nature of the different contents present in the image, therefore we considered different dictionaries based on fast transforms. For the artifacts related content we used an undecimated wavelet transform (UDWT) for modeling the dictionary Da, whereas for modeling the fibers dictionary Df, we used the curvelets transform. The latter being, as discussed in the introduction section, a natural choice for modeling curvilinear structures, while the wavelets transform allows modeling artifacts present in the images. Note that, in [23] the image component has been modelled by ridgelets and the texture (artefact) component by a Discrete Cosine Transform (DCT). Ridgelets has been also used in [50] to represent global lines in images.

Our implementation of the image decomposition is based on the MCALab library provided in [54], running a maximum of 100 iterations. The parameter [gamma] (in Eq 2) was linearly decremented during the iterations and initialized accordingly. For all reported experiments, we set p = 0 (in Eq 2), namely a [cursive l]0-norm for the model definition, and [delta] = 3. Results of the the image decomposition are illustrated in Fig 5. We refer the reader to [50] for more details on the solution algorithm, including initialization and updates of the [gamma] parameter (in Eq 2), and to the Supporting Information (SN. 1.1 in S1 Text) for a discussion on the impact of the different parameters values on the decomposition results.

Fig 5. The decomposition process of the original image f into the artifacts image va, the fibers image uf and noise .

(a) Image f. (b) Artifacts image va. (c) Fibers image uf. (d) Reminder noise . Solution obtained after 100 iterations with [delta] = 3.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g005

Filaments enhancement.

Once the fibers related image (uf) has been properly separated from other artifacts, the filaments are enhanced in order to make the fibers more visible with respect to the background. Likewise in [21], we used a sequence of filters to enhance the filaments, consisting in: a Gaussian filter (a convolution operator that uses a kernel representing the shape of a Gaussian, parameterized by [sigma] the standard deviation of the distribution); followed by a Laplace filter, with the following kernel: (3) where is a parameter that we usually set as 0 < [Beta] ≤ 10.0. In most of our experiments we set [Beta] = 10. The result of such operator is a sharpened image highlighting the edges. For the third preprocessing step, we apply a directional Gaussian filter, with a filter size of r = 2[left ceiling]3[sigma]dg[right ceiling] (next larger integer) [21] with [sigma]dg = 10.0 in most of our experiments. Fig 6 illustrates the results of this processing chain, denoted as uE, applied on the fibers image uf.

Fig 6. Fibers enhancement of the fibers image uf provides sharpened filaments network.

(a) Fibers image uf after sparse multi-source separation. (b) Filaments enhanced image, uE, obtained with [sigma] = 1.0, [Beta] = 10.0 and [sigma]dg = 10.0, for the Gaussian, Laplace and directional Gaussian, respectively.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g006

Multi-scale Line Segmentation

For the filaments network segmentation, we opted for a multi-scale line detector based on a structural element of different orientations and widths, representing the scales. The multi-scale line detector basically analyze each pixel's neighborhood of uE at different scales by evaluating (according to a score) if such pixel is part of a line of certain width. Such evaluation is performed by computing a 'line response' for the width evaluation, but also for a length evaluation that will be used in the next processing step for individual lines segments detection, as illustrated in Fig 7.

Fig 7. Linear response for width and length evaluation at certain orientation.

(a) Width evaluation. (b) Length evaluation.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g007

The line response takes place at each pixel in the image, in a (discrete) set of possible orientations between 0 and 180 degrees, considering linear elements of size s ∈ [1, W], with W the expected fibers width (this parameter is related to the magnification of the image). Then, the final multi-scale line response of each pixel (x, y) of the enhanced fibers image, uE, is provided as [34]: (4) where R(.), applied on the fibers enhanced image uE, provides a score indicating how likely is certain point to be the center of a line passing though it with a certain width and a certain direction, by analyzing the neighboring points. The reader is referred to [34] for the details on how R(.) is computed.

In our implementation, based on the source code provided by [34], we retained the advantages of the multi-scale linear response image, uG, and at the same time introduced some modifications to obtain the final binary image, uB of the fibers-related pixels. For the latter, we applied on the uG gray-scale image a local thresholding algorithm to separate background and fibers pixels. For this purpose, we made use of the Wellner's adaptive thresholding [55] were a median filtering provides an estimated local threshold. The final retained local threshold is a percentage b of the estimated threshold. The smaller b is the more line segments candidates are retained in the final binary image, representing a set of edge segments, which are clean, contiguous, i.e., connected, chains of edge pixels. Fig 8 illustrates the output of this step. Note that the obtained edges have a width larger than 1 pixel.

Fig 8. The binarization step provides the segmented filaments network.

(a) Multi-scale linear response uG. (b) Binary image uB obtained from uG. Width parameters set as W = 4.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g008

Line segmentation.

From the obtained binary image we perform a line segmentation step that provides a list of straight line segments with fixed-length, denoted by (F, [theta]

). Each line segment is represented by its parametric equation. The basic idea of the algorithm is to walk over the pixels in sequence, and fit lines to the pixels using the least squares line fitting method until the error exceeds a certain threshold, i.e. 1 pixel error. When the error exceeds this threshold, we generate a new line segment. The line segment is kept if its length is (at least) equal to L, a predefined parameter. In most of our experiments this parameter was set to L = 30.

Filament Segment Merging--Fibers Extraction

Several complex approaches of line grouping, such as the on in [56], have been proposed within the computer vision community. However, for our purpose of fibers extraction, we developed a simple algorithm capable of iteratively extracting continuous linear segments (denoted as filament segments) by connecting the fixed-length segments extracted in the previous step. We first associate to a given filament segment all overlapping fixed-length segments oriented in the same direction (same [theta]

). By repeating this process, all the overlapping segments with the same orientation will be combined into a longer straight-line segment. We then connect segments according to their orientation difference up to a 'curvature' threshold [theta]

< T[theta]

. In addition, when merging (i.e. connecting) the k-th and the i-th segments (Fk[N-ARY UNION]Fi), we discard all the pixels that are beyond the connection (intersection) point. The above described procedure is detailed in Algorithm (1), and illustrated in Fig 9.

Fig 9. Fixed length segments are merged if they overlap and their angle difference is not too large [theta]

< T[theta]

.

(a) Filament partitioned into fixed-length segments. (b) Overlapping and angle difference between segments.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g009

Algorithm 1 Fibers extraction

[theta]

i: Indicates the orientation of extracted fiber Fi

[delta](Fi): Indicates the length of extracted fiber Fi

Input: T[theta]

as angle merging tolerance

(F, [theta]

)[arrow left] List of fixed-length extracted segments tuples (Line segmentation section)

for [xi] = 0...T[theta]

0

for i = 1... do

OverlapEndPoints(Fk, Fi)}

if S [not equal] ∅ then

k [arrow left] S(1)

end if

end for

end for

Output: list of individual fibers (F, [theta]

)

Note: Fk[N-ARY UNION]Fi combines the segments Fk and Fi discarding all the pixels located beyond the connection (intersection) point.

Results and Discussion

Fig 10 illustrates the results obtained at each step of the proposed framework (Fig 3). The robustness of the framework is qualitatively illustrated using a defocused image in Fig 10a and properly focused image in Fig 10e. In the defocused image of Fig 10a, only the lower part of the images exhibit clear (visually) filaments, and in other parts of the image the filaments are barely visible. As result of the image decomposition and multi-scale line segmentation, we obtain most of the filaments (Fig 10d). More details on the different steps and the impact of their parameters stetting are given in S1 Text.

Fig 10. Individual filaments are extracted from the binarized image.

(a) Input defocused image f. (b) Fibers image uf. (c) Enhanced fibers image uE. (d) Fibers from uB. (e) Input well focused image x. (f) Fiber image uf. (g) Enhanced fibers image uE. (h) Fibers from uB. For each processing stage, the parameters were set as in Figs 5, 6 and 8, respectively.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g010

In order to qualitatively and quantitatively assess the performances of the proposed framework, we considered two datasets. The first dataset, from [21], consist in a benchmark of 10 simulated images (referred to as S1, ..., S10) and one real specimen cytoskeleton image (referred to as B2). The second dataset consist in 68 images of osteoblast cells (MC3T3-E1 cell line from mouse) grown in the two different conditions, static (control) and fluid shear stress, as described in section Materials and Methods.

Qualitative Validation

Cytoskeleton image B2.

From the data provided by [21], we selected the image of Fig 11a (referred to as B2). This challenging image has been selected as it exhibits highly blurred regions, with curvilinear filaments in addition to quasi-straight ones. Fig 11 illustrate the extracted filaments obtained using the proposed framework, compared to the method of [21]. A zoom-in of a very blurred image area is also given in Fig 12, where filaments indicated by the green arrows were missed by the method of [21], whereas the proposed framework succeeded to detect them.

Fig 11. The proposed framework extract a higher number of filaments.

(a) Original image B2. (b) Method of [21]. (c) Proposed framework with W = 4 and L = 30.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g011

Fig 12. Zoomed area highlighting blurred filaments detected by the proposed framework.

(a) Zoomed image B2. (b) Method of [21]. (c) Proposed framework.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g012

The parameters W and L play an important role in reducing false positives. Note that results in Fig 11c exhibit some false positives, while in Fig 13c most of the filaments have been extracted with less number of false positives.

Fig 13. Filaments tracing.

(a) Original image B2. (b) Method of [21]. (c) Proposed framework with W = 8 and L = 60.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g013

Peri-nuclear zone.

The central role of actin cap in cells mechanosensing and mechanotransduction has been recently discovered [11, 12, 57, 58]. Thus, a proper filaments extraction on this area can be important to detect significant changes as result of mechanical stimulation. Images of actin filaments in the peri-nuclear zone are likely to exhibit some degree of blurring. Therefore, we assessed the framework performance in the peri-nuclear area of osteoblasts grown under fluid shear stress conditions. Fig 14 illustrates some of the extracted fibers using our method, compared to the results obtained using the method of [21].

Fig 14. Extracted fibers in the peri-nuclear area.

(a) Method of [21]. (b) Proposed framework.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g014

As can be seen, compared to [21], the proposed framework is able to find most of the visible filaments. Although both methods detected several filaments in the bottom-left cell, only the proposed approach detected most of the filaments located in the peri-nuclear area of the top-right cell. Indeed, thanks to the image decomposition approach, several blurred quasi-straight clues were associated as part of filaments by the proposed framework.

Quantitative Evaluation

Simulated filament networks.

In order to quantitatively measure the model performance, we considered a set of 10 simulated images available in the benchmark database of [21], which consist in a set of filaments of different widths and lengths, corrupted with noise.

From the experimental results we computed the model accuracy (Acc), sensitivity (Sn) and specificity (Sp), defined as: (5) (6) (7) where TP, TN, FP, FN stands for:

* TP: pixels correctly identified as part of the filaments network

* FP: pixels wrongly identified as part of the filaments network

* TN: pixels correctly identified as not being part of the filaments network

* FN: pixels wrongly identified as not being part of the filaments network

Results of the method in [21] were publicly available for the benchmark database, thus we only had to run experiments for the proposed method; Fig 15 illustrates both results in one of the simulated images of [21].

Fig 15. Most of errors accounted in the proposed framework consist of fibers longer than they should.

(a) Simulated image (S1). (b) Method of [21]. (c) Proposed framework.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g015

The models accuracy and sensitivity are shown in Fig 16 for the 10 simulated images. As can be seen, both models perform similarly in terms of accuracy, although the proposed framework exhibited a little bit higher accuracy in 6 of the 10 images. In the proposed method, most of the incurred errors consist in false positive detections related to elongated fibers, and to a lesser extent some undetected thin filaments (false negatives).

Fig 16. The models exhibit similar accuracy but a higher sensitivity, in 10 simulated images; Stars corresponds to method in [21] and Circles to the proposed framework results.

(a) Accuracy. (b) Sensitivity.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g016

However, the proposed model has a considerably higher sensitivity and a relatively similar (though lower) specificity, as detailed in Tables 1 and 2. This means that our model detected a higher number of filament pixels incurring in some false positive detections (mostly elongated fibers), while the model in [21] incurred in the same number of errors by failing to detect filaments.

Table 1. Sensitivity (Sn) comparison.

Top-row: the method of [21]; Bottom-row: the proposed framework (PF).

http://dx.doi.org/10.1371/journal.pcbi.1005063.t001

Table 2. Specificity (Sp) comparison of method [21] in top-row and the proposed framework (PF) in bottom-row.

http://dx.doi.org/10.1371/journal.pcbi.1005063.t002

With respect to the specificity, there are no major differences, mostly due to the large number of background (negative) pixels present in the image; differences in sensitivity are significant.

Visually annotated osteoblasts images.

To asses the suitability of the proposed framework for filament's orientation analysis, a set of osteoblats images/sub-images were visually annotated into upper-left oriented, upper-right oriented, or other orientation, according to the orientation of their fibers. From the annotated images, we randomly selected 18 images, 9 showing upper-left oriented fibers, denoted as Left-Set, and the other 9 showing upper-right fibers, denoted as Right-Set; Fig 17 depicts some images from each set. As it can be noticed, some of the selected images are heavily blurred, which can be handled by our framework (see Fig 10a-10d) thanks to the proposed image decomposition approach.

Fig 17. Visually annotated osteoblasts for orientation validation.

(a-d) Images with fibers oriented to the left. (e-h) Images with fibers with fibers oriented to the right.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g017

The normalized angular distribution is depicted in Fig 18. It is computed by counting the number of detected individual fibers oriented toward each discrete direction, and dividing them by the maximum value. As it can be seen, the estimated fiber's orientations correspond to the visually observed orientations. This confirms the suitability of the proposed framework for filament's orientation analysis.

Fig 18. Normalized angular distribution of the Left-Set, and Right-Set, respectively, considering the horizontal axis of the image as reference.

(a) Left-oriented fibers. (b) Right-oriented fibers.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g018

Effect of shear stress on actin filaments in osteoblasts.

In order to demonstrate the effectiveness of the proposed framework, the analysis of orientation, fibers quantification and length was performed in two population of osteoblasts grown in static conditions and under fluid shear stress (as described in Section Fluid shear stress induction). Application of shear stress to the cells is expected to cause reorientation of the actin stress fibers in a way that the filaments become nearly aligned with the shear stress flow direction [59, 60]. In our experiments, the osteoblast population grown under fluid shear stress contained a total of 32 images. The resultant filaments angular distribution over all the 32 images subjected to shear flow is illustrated in Fig 19. Fig 20 illustrates two osteoblast images and the detected fibers, in static and fluid- shear stress conditions.

Fig 19. Normalized angular distribution of fibers grown in different stress conditions, taking as reference the horizontal axis of the image, and considering all images of each population.

(a) Osteoblasts grown in static conditions. (b) Osteoblasts grown under fluid shear stress.

http://dx.doi.org/10.1371/journal.pcbi.1005063.g019

Fig 20. The proposed framework can properly extract highly overlapping fibers.

(a) and (d) Osteoblast image. (b) and (e) Extracted fibers located near the nucleus. (c) and (f) Detected filaments network. Top row: Osteoblasts grown in normal conditions; Bottom row: Osteoblasts grown under fluid shear stress. The arrows depict the shear stress flow direction around [asymptotically =] 80° (Fig 4b).

http://dx.doi.org/10.1371/journal.pcbi.1005063.g020

To highlight the complexity of extracting the filaments in a cytoskeleton, we selected the image (Fig 20d) among the 32 images taken in the shear stress experiment. The results confirm that only a few filaments are oriented diagonally to the flow direction, and most of the stress fibers are aligned with the flow direction (see Fig 19e). The angular distribution of fibers orientation (with mean M = 90.3 and variance [sigma]2 = 1334.90) follows the almost vertical flow direction imposed to the shear stress cell culture (see Fig 4b), thus validating the accuracy of the method. The vast majority of the filaments located near the nucleus (Fig 20e) are well aligned with the flow. A smaller number of filaments, properly detected by our method, exhibited an orientation nearly diagonal with respect to the flow. In addition, most of these filaments are located in the peri-nuclear region thus likely related to the actin cap, suggesting their association to some strain avoidance response [11].

The filaments of the static (control) cells population (38 images) seems to cover the whole angular spectrum (with mean M = 68.37 and variance [sigma]2 = 3071.86), albeit a preference towards some specific angular orientations (Fig 19a) can be seen. This preference could reflect the orientation of the cells after hours of static growth in the channel.

The shear stress experiment exhibited a small variance value ([sigma]2 = 1334.90) in the fibers orientation, which indicates that the vast majority of fibers are oriented towards the same direction. On the contrary, the variance of the static/control experiment ([sigma]2 = 3071.86) is more than twice bigger compared to the shear stress experiment, evidencing that fibers are oriented towards several angular directions.

Conclusions

The cytoskeleton plays an important role in numerous physiological and pathological processes, its morphological characteristics are therefore of prime importance to understand numerous basic cellular phenomena, such as cellular adaptation to physical or chemical stress. However, cytoskeleton quantification and analysis is far from being straightforward, and sophisticated algorithms are required to fulfill the task. In this work we present a processing framework that efficiently detects cytoskeletal fibers and quantifies its morphological characteristics such as the number of filaments it contains, their length and orientation. The proposed model was tested on images of osteoblasts cultivated in shear stress and static (control) conditions. The detection of highly oriented actin fibers in shear stress cultivated cells corroborates what one would expect in such a condition, i.e. an alignment of the fibers with the direction of the flow. In addition, our algorithm successfully detected the perinuclear actin cap, a structure difficult to detect by other state-of-the-art methods, and that seems to play an important role in mechano- transduction [11, 58]. It was also shown that the model separated filaments and imaging-related artifacts very efficiently, even in the presence of heavy blurring, a step that endowed the model with high sensitive detection capabilities. The proposed framework can be extended to extract 3D meshwork of actin filaments. Indeed, after the image decomposition and filament enhancements steps, a method such as the multiple Stretching Open Active Contours (SOACs) of [61] could be used for fibers tracing.

Supporting Information

S1 Fig. Effect of number of iteration on image decomposition results.

(a) Fibers image uf after 10 iterations. (b) Fibers image uf after 100 iterations. (c) Fibers image uf after 300 iterations. Those elements within the image that does not exhibit a filamentous geometry are progressively removed from the filaments component and put into the non-filament part.

(TIF)

S2 Fig. Synthetic image.

(a) Ground truth image u. (b) Generated Artifacts v. (c) Obtained synthetic image f = u + v + [eta] with [eta] a Gaussian noise [sigma] = 0.04.

(TIF)

S3 Fig. Image decomposition of the synthetic image of S2c Fig, using 100 iterations and [gamma] = 3.

(a) Fibers image uf. (b) Artifacts image va. (d) Reminder noise with estimated noise level [sigma] = 0.07.

(TIF)

S4 Fig. Different value of the filaments enhancement parameters (Gaussian filter, [sigma], Laplace filter, [Beta], and linear Gaussian, [sigma]dg) provides deferent filament saliencies.

(a) [sigma] = 0.5, [Beta] = 1.0 and [sigma]dg = 2.0. (b) [sigma] = 1.0, [Beta] = 5.0 and [sigma]dg = 5.0. (c)[sigma] = 1.0, [Beta] = 10.0 and [sigma]dg = 10.0. For visualization purposes, the brightness and contrast was regulated in exactly the same quantity on the three images.

(TIF)

S5 Fig. Filament tracing.

(a) Filaments enhancement of original image. (b) Filaments enhancement after image decomposition. (c) Multi-scale linear response (W = 2) followed by binarization step b = 10. (d) The line segmentation stage keeps only quasi-straight segments of a minimum length L = 30, discarding the others. (e) W = 2, b = 0.1 and L = 30. (f) W = 4, b = 0.1 and L = 30.

(TIF)

S6 Fig. Straight line segments obtained with T[theta]

= 0.

(a) Binary image. (b) Individual filaments. Only the 100 longest filaments are displayed. The different colors depict different orientations.

(TIF)

S7 Fig. Extracted filaments with T[theta]

= 2.

(a) Binary image. (b) Individual filaments. The different colors depict different orientations.

(TIF)

S1 Text. Supplementary Notes.

(PDF)

Author Contributions

1. Conceived and designed the experiments: MAP CB KG SK GD RW HS.

2. Performed the experiments: MAP CB KG.

3. Analyzed the data: MAP HS.

4. Contributed reagents/materials/analysis tools: MAP CB KG SK GD RW HS.

5. Wrote the paper: MAP CB KG SK GD RW HS.

References

Cai S, Pestic-Dragovich L, O'Donnell ME, Wang N, Ingber D, Elson E, et al. Regulation of cytoskeletal mechanics and cell growth by myosin light chain phosphorylation. The American journal of physiology-Cell Physiology. 1998;275(5):C1349-C1356

Bunnell TM, Burbach BJ, Shimizu Y, Ervasti JM. [Beta]-Actin specifically controls cell growth, migration, and the G-actin pool. Molecular Biology of the Cell. 2011;22(21): 4047-4058. doi: 10.1091/mbc.E11-06-0582. pmid:21900491

Sivaramakrishnan S, Spudich JA. Coupled myosin VI motors facilitate unidirectional movement on an F-actin network. The Journal of cell biology. 2009;187(1): 53-60. doi: 10.1083/jcb.200906133. pmid:19786577

Tseng Q, Duchemin-Pelletier E, Deshiere A, Balland M, Guillou H, et al. Spatial organization of the extracellular matrix regulates cell-cell junction positioning. Proceedings of the National Academy of Sciences of the United States of America. 2012;109(5): 1506-1511. doi: 10.1073/pnas.1106377109. pmid:22307605

Stanisavljevic J, Loubat-Casanovas J. Snail1-expressing fibroblasts in the tumor microenvironment display mechanical properties that support metastasis. Cancer research. 2015; 75(2): 284-295. doi: 10.1158/0008-5472.CAN-14-1903. pmid:25488750

McBeath R, Pirone DM, Nelson CM, Bhadriraju K, Chen CS. Cell Shape, Cytoskeletal Tension, and RhoA Regulate Stem Cell Lineage Commitment. Developmental Cell. 2004;6(4): 483-495. doi: 10.1016/S1534-5807(04)00075-9. pmid:15068789

Treiser MD, Yang EH, Gordonov S, Cohen DM, Androulakis IP, Kohn J, et al. Cytoskeleton-based forecasting of stem cell lineage fates. Proceedings of the National Academy of Sciences of the United States of America. 2010;107(2): 610-615. doi: 10.1073/pnas.0909597107. pmid:20080726

Kilian Ka, Bugarija B, Lahn BT, Mrksich M. Geometric cues for directing the differentiation of mesenchymal stem cells. Proceedings of the National Academy of Sciences of the United States of America. 2010;107(11): 4872-4877. doi: 10.1073/pnas.0903269107. pmid:20194780

Prowse PD, Elliott CG, Hutter J, Hamilton DW. Inhibition of Rac and ROCK signalling influence osteoblast adhesion, differentiation and mineralization on titanium topographies. PloS One. 2013;8(3). doi: 10.1371/journal.pone.0058898. pmid:23505566

Kim HY, Davidson LA. Punctuated actin contractions during convergent extension and their permissive regulation by the non-canonical Wnt-signaling pathway. Journal of cell science. 2011;124(4): 635-646. doi: 10.1242/jcs.067579. pmid:21266466

Tamiello C, Bouten CVC, Baaijens FPT. Competition between cap and basal actin fiber orientation in cells subjected to contact guidance and cyclic strain. Scientific Reports. 2015;5. doi: 10.1038/srep08752. pmid:25736393

Tojkander S, Gateva G, Lappalainen P. Actin stress fibers--assembly, dynamics and biological roles. Journal of Cell Science. 2012;125(8): 1855-1864. doi: 10.1242/jcs.098087. pmid:22544950

Jackson WM, Jaasma MJ, Tang RY, Keaveny TM. Mechanical loading by fluid shear is sufficient to alter the cytoskeletal composition of osteoblastic cells. American journal of physiology-Cell physiology. 2008;295(4): C1007-C1015. doi: 10.1152/ajpcell.00509.2007. pmid:18701650

Gardinier JD, Majumdar S, Duncan RL, Wang L. Cyclic Hydraulic Pressure and Fluid Flow differentially modulate cytoskeleton re-organization in MC3T3 Osteoblasts. Cellular and molecular bioengineering. 2010;2(1): 133-143. doi: 10.1007/s12195-008-0038-2

Nabavi N, Khandani A, Camirand A, Harrison RE. Effects of microgravity on osteoclast bone resorption and osteoblast cytoskeletal organization and adhesion. Bone. 2011;49(5): 965-974. doi: 10.1016/j.bone.2011.07.036. pmid:21839189

Dai Z, Wu F, Chen J, Xu H, Wang H, Guo F, et al. Actin microfilament mediates osteoblast Cbfa1 responsiveness to BMP2 under simulated microgravity. PloS One. 2013;8(5). doi: 10.1371/journal.pone.0063661

Kaunas R, Usami S, Chien S. Regulation of stretch-induced JNK activation by stress fiber orientation. Cellular Signalling. 2006;18(11): 1924-1931. doi: 10.1016/j.cellsig.2006.02.008. pmid:16581230

Zhao L, Sang C, Yang C, Zhuang F. Effects of stress fiber contractility on uniaxial stretch guiding mitosis orientation and stress fiber alignment. Journal of Biomechanics. 2011;44(13): 2388-2394. doi: 10.1016/j.jbiomech.2011.06.033. pmid:21767844

Herberich G, Würflinger T, Sechi A, Windoffer R, Leube R, Aach T. Fluorescence microscopic imaging and image analysis of the cytoskeleton. In: Asilomar Conference on Signals, Systems and Computers; 2010. p. 1359-1363

Elosegui-Artola A, Jorge-Peñas A, Moreno-Arotzena O, Oregi A, Lasa M, et al. Image Analysis for the Quantitative Comparison of Stress Fibers and Focal Adhesions. PLoS One. 2014;9(9). doi: 10.1371/journal.pone.0107393. pmid:25269086

Eltzner B, Wollnik C, Gottschlich C, Huckemann S, Rehfeldt F. The Filament Sensor for Near Real-Time Detection of Cytoskeletal Fiber Structures. Plos One. 2015;10(5). doi: 10.1371/journal.pone.0126346. pmid:25996921

Desolneux A, Moisan L, Morel J-M. From Gestalt Theory to Image Analysis: a probabilistic approach. Springer Science & Bussiness Media; 2008(34). doi: 10.1007/978-0-387-74378-3

Starck JL, Elad M, Donoho D. Image decomposition via the combination of sparse representation and a variational approach. IEEE Transactions on Image Processing. 2005;14(10): 1570-1582. doi: 10.1109/TIP.2005.852206. pmid:16238062

Vese L., Osher S. Modeling textures with total variation minimization and oscillating patterns in image processing. Journal of Scientific Computing. 2003;19:553-577. doi: 10.1023/A:1025384832106

Aujol F., Gilboa G., Chan T., Osher S. Structure-texture image decomposition-modeling, algorithms and parameter selection. Int. Journal Computer Vision. 2006;67(1):111-136. doi: 10.1007/s11263-006-4331-z

Candés E. J., Demanet L., Donoho D., Ying L. Fast discrete curvelet transforms. Multiscale Model. Simul. 2006;5(3):861-899. doi: 10.1137/05064182X

Candés E. J., Donoho D. L. New tight frames of curvelets and optimal representations of objects with piecewise C2 singularities. Communications on Pure and Applied Mathematics. 2004;57(2):219-266. doi: 10.1002/cpa.10116

Akakin H.C., Gurcan M.N. Content-based microscopic image retrieval system for multi-image queries. IEEE Trans. Inform. Technol. Biomed. 2012;16:758-769. doi: 10.1109/TITB.2012.2185829

Danuser G. Computer vision in cell biology. Cell. 2011;147(5): 973-978. doi: 10.1016/j.cell.2011.11.001. pmid:22118455

Marquez JP. Fourier analysis and automated measurement of cell and fiber angular orientation distributions. International Journal of Solids and Structures. 2006;43(21): 6413-6423. doi: 10.1016/j.ijsolstr.2005.11.003

Matschegewski C, Staehlke S, Birkholz H, Lange R, Beck U, Engel K, et al. Automatic Actin Filament Quantification of Osteoblasts and Their Morphometric Analysis on Microtextured Silicon-Titanium Arrays. Materials. 2012;5(12): 1176-1195. doi: 10.3390/ma5071176

Rafael GVG, Jérémie J, Jean-Michel M, Gregory R. LSD: a Line Segment Detector. Image Processing On Line. 2012;2: 35-55. doi: 10.5201/ipol.2012.gjmr-lsd

Weichsel J, Urban E, Small JV, Schwarz US. Reconstructing the orientation distribution of actin filaments in the lamellipodium of migrating keratocytes from electron microscopy tomography data. Cytometry Part A: the journal of the International Society for Analytical Cytology. 2012;81(6): 496-507. doi: 10.1002/cyto.a.22050

Nguyen UTV, Bhuiyan A, Park LA, Ramamohanarao K. An effective retinal blood vessel segmentation method using multi-scale line detection. Pattern Recognition. 2013;46(3): 703-715. doi: 10.1016/j.patcog.2012.08.009

Sampo J, Takalo J, Siltanen S, Miettinen A, Lassas M, Timonen J. Curvelet-based method for orientation estimation of particles from optical images. Optical Engineering. 2014;53(3): 033109-033109. doi: 10.1117/1.OE.53.3.033109

Moller B, Piltz E, Nadine B. Quantification of Actin Structures Using Unsupervised Pattern Analysis Techniques. In: International Conference on Pattern Recognition; 2014. p. 3251-3256

Weichsel J, Schwarz US. Two competing orientation patterns explain experimentally observed anomalies in growing actin networks. Proceedings of the National Academy of Sciences of the United States of America. 2010;107(14): 6304-6309. doi: 10.1073/pnas.0913730107. pmid:20308581

Moch M, Herberich G, Aach T, Leube RE, Windoffer R. Measuring the regulation of keratin filament network dynamics. Proceedings of the National Academy of Sciences of the United States of America. 2013;110(26): 10664-10669. doi: 10.1073/pnas.1306020110. pmid:23757496

Rigort A., Günther D., Hegerl R., Baum D., Weber B., Prohaska S., Medalia O., Baumeister W., Hege H.C. Automated segmentation of electron tomograms for a quantitative description of actin filament networks. J. Struc. Biol. 2012;177: 135-144 doi: 10.1016/j.jsb.2011.08.012

Rusu M., Starosolski Z., Wahle M., Rigort A, Wriggers W. Automated tracing of filaments in 3D electron tomography reconstructions using Sculptor and Situs. J. Struc. Biol. 2012;178: 121-28. doi: 10.1016/j.jsb.2012.03.001

Xu T., Vavylonis D., Tsai F.C., Koenderink G.H., Nie W., Yusuf E., Lee I.J., Wu J.Q., Huang X. SOAX: a software for quantification of 3D biopolymer networks. Sci. Rep. 2015;5: 9081. doi: 10.1038/srep09081

Schaub S, Meister JJ, Verkhovsky AB. Analysis of actin filament network organization in lamellipodia by comparing experimental and simulated images. Journal of cell science. 2007;120(8): 1491-1500. doi: 10.1242/jcs.03379. pmid:17401113

Kim DH, Khatau SB, Feng Y, Walcott S, Sun SX, Longmore GD, et al. Actin cap associated focal adhesions and their distinct role in cellular mechanosensing. Scientific Reports. 2012;2. doi: 10.1038/srep00555

Winkler C, Vinzenz M, Small JV, Schmeiser C. Actin filament tracking in electron tomograms of negatively stained lamellipodia using the localized radon transform. Journal of Structural Biology. 2012;178(1): 19-28. doi: 10.1016/j.jsb.2012.02.011. pmid:22387240

Koestler SA, Steffen A, Nemethova M, Winterhoff M, Luo N, Holleboom JM, et al. Arp2/3 complex is essential for actin network treadmilling as well as for targeting of capping protein and cofilin. Molecular biology of the cell. 2013;24(18): 2861-2875. doi: 10.1091/mbc.E12-12-0857. pmid:23885122

Soiné JRD, Brand CA, Stricker J, Oakes PW, Gardel ML, Schwarz US. Model-based Traction Force Microscopy Reveals Differential Tension in Cellular Actin Bundles. PLOS Computational Biology. 2015;11(3). doi: 10.1371/journal.pcbi.1004076. pmid:25748431

van Vliet LJ, Young IT, Beckers GL. A nonlinear laplace operator as edge detector in noisy images. Computer Vision, Graphics, and Image Processing. 1989;45(2): 167-195. doi: 10.1016/0734-189X(89)90131-X

Freeman WT, Adelson EH. The design and use of steerable filters. IEEE Transactions on Pattern Analysis & Machine Intelligence. 1991; 13: 891-906. doi: 10.1109/34.93808

Weickert J. Coherence-enhancing shock filters. DAGM-Symposium. 2003; 1-8. doi: 10.1007/978-3-540-45243-0_1

Fadili JM, Starck JL, Elad M, Donoho DL. MCALab: Reproducible research in signal and image decomposition and inpainting. Computing in Science and Engineering. 2010; 12(1): 44-63. doi: 10.1109/MCSE.2010.14

Tschumperlé D, Fadili J, Bentolila Y. Wire structure pattern extraction and tracking from X-ray images of composite mechanisms. In: Computer Vision and Pattern Recognition; 2010. 2: 2461-2466 doi: 10.1109/cvpr.2006.335

Otsu N. A Threshold Selection Method from Gray-Level Histograms. Automatica. 1975;11: 23-27

Zemel A, Rehfeldt F, Brown aEX, Discher DE, Safran SA. Optimal matrix rigidity for stress-fibre polarization in stem cells. Nature Physics. 2010;6(6): 468-473. doi: 10.1038/nphys1613. pmid:20563235

Signal and Image Decomposition and Inpainting v12.0. https://fadili.users.greyc.fr/demos/WaveRestore/downloads/mcalab/Download.html. (Accessed 30 September 2015)

Wellner PD. Adaptive thresholding for the DigitalDesk. Xerox. 1993;EPC1993-110

Katarztis A., Sahli H., Pizurica V., Cornelis J. A Model-Based Approach to the Automatic Extraction of Linear Features from Airborne Images. IEEE Trans. on Remote Sensing and Geosciences. 2001;39(3):2073-2079. doi: 10.1109/36.951102

Khatau SB, Hale CM, Stewart-Hutchinson PJ, Patel MS, Stewart CL, et al. A perinuclear actin cap regulates nuclear shape. Proceedings of the National Academy of Sciences of the United States of America. 2009;106: 19017-19022. doi: 10.1073/pnas.0908686106. pmid:19850871

Kim DH, Chambliss AB, Wirtz D. The multi-faceted role of the actin cap in cellular mechanosensation and mechanotransduction. Soft Matter. 2013;9: 5516-5523. doi: 10.1039/c3sm50798j

Yoshigi M, Hoffman LM, Jensen CC, Yost HJ, Beckerle MC.J Mechanical force mobilizes zyxin from focal adhesions to actin filaments and regulates cytoskeletal reinforcement. Cell Biol. 2005;171(2): 209-215. doi: 10.1083/jcb.200505018

Azuma N., Akasaka N., Kito H., Ikeda M., Gahtan V., Sasajima T., Sumpio B. E. Role of p38 MAP kinase in endothelial cell alignment induced by fluid shear stress. American Journal of Physiology-Heart and Circulatory Physiology. 2001;280(1):H189-H197. pmid:11123233

Xu T., Vavylonis D., Huang X. 3D actin network centerline extraction with multiple active contours. Medical Image Analysis. 2014;18(2):272-284. doi: 10.1016/j.media.2013.10.015. pmid:24316442