ABSTRACT
Cell polarity involves the compartmentalization of the cell cortex. The establishment of cortical compartments arises from the spatial bias in the activity and concentration of cortical proteins. The mechanistic dissection of cell polarity requires the accurate detection of dynamic changes in cortical proteins, but the fluctuations of cell shape and the inhomogeneous distributions of cortical proteins greatly complicate the quantitative extraction of their global and local changes during cell polarization. To address these problems, we introduce an open-source software package, ImaEdge, which automates the segmentation of the cortex from time-lapse movies, and enables quantitative extraction of cortical protein intensities. We demonstrate that ImaEdge enables efficient and rigorous analysis of the dynamic evolution of cortical PAR proteins during Caenorhabditis elegans embryogenesis. It is also capable of accurate tracking of varying levels of transgene expression and discontinuous signals of the actomyosin cytoskeleton during multiple rounds of cell division. ImaEdge provides a unique resource for quantitative studies of cortical polarization, with the potential for application to many types of polarized cells.
This article has an associated First Person interview with the first authors of the paper.
INTRODUCTION
The cell cortex consists of a thin layer of cytoplasmic proteins and the actomyosin cytoskeleton organized beneath the plasma membrane (Salbreux et al., 2012). The cortex can be compartmentalized via the spatially biased dynamic behaviors of the actomyosin cytoskeleton and cortical proteins (Bray and White, 1988; Clark et al., 2014). These cortical proteins can diffuse laterally within the cortex and exchange between the cortex and the cytoplasm. The polarization of these cortical compartments, in turn, controls various aspects of metazoan development. For example, the establishment of apical and basolateral cortical domains in epithelial cells allows for directional vesicle transport (Rodriguez-Boulan and Macara, 2014). Axons and dendrites also arise from asymmetrical segregation of cortical proteins in neuronal cells, allowing for the establishment of neural circuits (Knoblich, 2010).
The understanding of mechanisms underlying cortical polarization has long been a goal in the study of cell polarity. Robust quantitative analysis of protein dynamics at the cortex is integral towards this objective, but has remained technically challenging. While sophisticated imaging techniques have enabled cell and developmental biologists to attain higher-resolution images of cells, the dynamic process of cell polarization has not been sufficiently described. This is largely due to a lack of development of unbiased quantitative analyses of dynamic protein behaviors at the cortex (Myers, 2012). The spatiotemporal dynamics of cortical proteins are commonly represented either by a gallery of consecutive frames from time-lapse movies or by a kymograph with straight lines drawn at multiple locations of interest. Although they are simple to perform, these approaches are far from able to capture the full extent of information available in time-lapse movies. There are several computational tools that can track the position of cell edges and cortical proteins, but they are mostly optimized for cultured cells, which are typically labeled by overexpressed bright fluorescent probes with a largely uniform distribution (Tables S1 and S2) (Barry et al., 2015; Biro et al., 2013; Bosgraaf and Van Haastert, 2010; Dimopoulos et al., 2014; Dormann et al., 2002; Held et al., 2010; Hodneland et al., 2013; Machacek et al., 2009; Tsygankov et al., 2014; Wang et al., 2010). It has been a challenge to apply the existing tools to study cortical proteins in multicellular transgenic organisms such as developing embryos, mainly because of the weak fluorescence signals expressed in these systems, together with the complex and fluctuating cell contours. Additionally, although dedicated plasma membrane probes, such as GFP-tagged PH domain, can be used to aid in cortex segmentation (Balla and Varnai, 2009), the use of such probe takes up one valuable spectral channel, restricting the choice of fluorophores and constraining multi-color multi-protein analysis.
To overcome the aforementioned problems, we sought to develop a simple but robust method that can precisely locate the cell edge contours using the inhomogeneous distribution of fluorescence intensities of a polarity protein itself. This is an important prerequisite for investigating certain proteins that exhibit polarized distribution along the cortex, such as cortical polarity regulators and the actomyosin cytoskeleton (Goldstein and Macara, 2007; St Johnston and Ahringer, 2010). In the case of these proteins, segmentation using a single threshold value cannot capture the proper cell edges for the entire duration of cell polarization, during which the intensities and patterns of cortical proteins will change. These problems could be addressed in principle by using segmentation methods based on the progressive evolution of an active contour that is defined either by the minimization of a customized energy function or by parameter-free methods that track the local information from the cell edge (Ersoy et al., 2008a,b; Li et al., 2009; Nejati Javaremi et al., 2013; Song et al., 2017; Wu et al., 2015). These methods in practice involve relatively complex segmentation algorithms with multiple key input parameters, requiring significant user intervention to analyze large-scale datasets from long time-lapse movies, such as those obtained during cell polarization. Therefore, the important segmentation steps tend to be labor intensive, time consuming and error prone.
Here, we present a computational approach, implemented as a software called ImaEdge (IMAge analysis of cell EDGE dynamics), which addresses the problems highlighted above. Instead of requiring a global threshold for defining the cortex, ImaEdge utilizes an adaptive gradient method to automatically determine the relevant parameters for all frames of a time-lapse movie, using the fluorescence intensity of the cortical protein of interest itself (Fig. 1). Firstly, the images are radially segmented using adaptive thresholds. The optimal cell contour is computed in accordance with the maximum average image intensity gradient without the creation of a complex energy function. Secondly, to detect the thickness of the cortical region, the cell contour is propagated inward until the maximum average intensity gradient is located. Thirdly, to obtain a global view of the dynamics and spatiotemporal evolution of cortical protein intensity, the cortex is sampled at specified intervals, and the temporal evolution is calculated by closest-sub-region linkage between adjacent frames. We demonstrate the capability of ImaEdge by quantifying the spatiotemporal dynamics of two cell polarity proteins, PAR-6 and PAR-2, and non-muscle myosin-II heavy chain, NMY-2, during cell polarization in Caenorhabditis elegans embryos. ImaEdge was able to trace the local and global changes of cortical protein concentrations from one-cell to two-cell embryos without the need for dedicated cortical markers. ImaEdge is equipped with a graphical interface for ease of use by cell and developmental biologists [software user manual available on GitHub (http://github.com/KanchanawongLab)], thus providing a unique resource for a quantitative study of cortical polarization in various polarized cell types.
RESULTS
Monitoring cortical polarization in C. elegans zygotes
To develop a method to automatically determine the cell cortex region in every frame of a time-lapse movie, we used a classic example of cell polarization, the C. elegans zygote. The C. elegans zygote begins to polarize shortly after the entry of the sperm-donated centrosome (Cowan and Hyman, 2004, 2006; Hamill et al., 2002; O'Connell et al., 2000; Schumacher et al., 1998). Cortical actomyosin and partitioning defective (PAR) proteins become segregated into two distinct cortical domains along the anterior–posterior axis (Kemphues, 2000). The actomyosin network and the anterior PAR proteins (two PDZ-domain proteins PAR-3 and PAR-6, and atypical protein kinase C) initially localize throughout the cortex and then become segregated to the anterior cortex during polarization (Etemad-Moghadam et al., 1995; Hung and Kemphues, 1999; Tabuse et al., 1998). The posterior PAR proteins (PAR-1 kinase, PAR-2, and LGL-1) translocate from the cytoplasm onto the posterior cortex (Beatty et al., 2010; Boyd et al., 1996; Guo and Kemphues, 1995; Hoege et al., 2010). Extensive overlap between these two cortical domains is prevented by reciprocal inhibitory interactions between the anterior and the posterior PAR proteins (Cuenca et al., 2003; Hao et al., 2006; Hoege et al., 2010; Motegi et al., 2011). Since PAR-6 and PAR-2 have commonly been used as proxies for the anterior and the posterior PAR protein complexes, respectively (Cuenca et al., 2003; Goehring et al., 2011), we monitored mCherry::PAR-6 and GFP::PAR-2 in live zygotes (Fig. S1 and Movie 1).
Adaptive radial segmentation of cell edges
The non-uniform distribution of cortical PAR proteins poses a major challenge for the accurate and unbiased segmentation of the cell cortex (Fig. 2A–C). Therefore, ImaEdge was designed to automate these processes while providing intuitive visual guides to allow a biologist full control of the process. ImaEdge first enhances the intensity contrast of the loaded image (MATLAB command: imadjust) as an initial processing step before subsequent segmentation. We implemented an adaptive radial segmentation algorithm to determine the contour of the cell edge in each frame of the time-lapse movies (Fig. 2D–I). The image of the first frame (either single-channel, or merged multi-channel) was binarized using Otsu's method (Otsu, 1979), keeping the largest 8-connected binary region as the initial binary mask for the cell of interest, while discarding smaller regions. This yields the initial mask that captures the majority of the cell region, albeit with irregular edges that do not generally conform to the actual cell edge contour. Next, to obtain the accurate cell contour, the center of mass of the cell of interest was calculated. This position was used to radially sub-divide the region of the cell of interest into multiple sectors (typically 10) (Fig. 2D). A small margin (typically 35 pixels or 9.3 µm) was added to each sector to include a small portion of the out-of-cell background regions for accurate estimates of the background intensity. For each sector, a binary threshold between the cell and the background was calculated using Otsu's method (Fig. 2E–G). This serves as the starting estimate for the cell contour within a given sector. Subsequently, the largest 8-connected objects in all sectors were then merged into the whole-cell binary mask. The merged mask was median-filtered (typically a 10×10 pixels median filter) to smoothen the cell contour, especially to eliminate disjointed edges at the boundaries between sectors (Fig. 2H,I). The size of the median filter determines the extent to which the cell contour will be smoothened. For example, a larger median filter focuses more on the global cell morphology, while a smaller one better preserves local context along the cell edge. To find the optimal binary thresholds for these sectors, varying scaling factors within a customized range (typically from 50% to 120% in 2% increments) were applied to the initial Otsu threshold. To segment the cortical region, the intensity gradient map was calculated by applying a Sobel–Feldman filter to the contrast-enhanced image. The averaged image gradient with a specific contour was then calculated.
The optimal scaling factor was determined as the value which yields the maximum average image gradient. This was then multiplied with the Otsu threshold to generate the optimal threshold levels for all sectors. The same procedure is then applied to each time-lapse frame, thereby automatically segmenting the cell contour for the entire movie.
Quantification of cortical protein spatiotemporal dynamics
Once the cell edge contours were defined, the thickness of the cortex was determined by computing a distance map relative to the segmented cell edge (Fig. 3 and Fig. 4A,E). The distance levels were propagated inward, allowing the calculation of the average image gradient for each distance level (Fig. 3B,C). The minimum distance that fully covers the first and second maximum spatial derivatives of the image gradient was then selected as the outer and inner boundaries of the cell cortex (Fig. 3D–G). Since cortical thickness may vary in different frames, this process was then iterated for all image frames, thereby automatically capturing the cortical width dynamics of the entire cell (Fig. 3H).
To quantify the spatiotemporal dynamics of fluorescence-tagged proteins at the cortex, the entire cortical region was subdivided into a user-defined number of sampling windows (typically 100). This was done by placing uniformly spaced sampling points along the cell contour, which are then connected to their closest partner points along the inner boundary (Fig. 4A–H; Movie 2). The nearest box-oriented alignment was performed to align the leftmost sampled window to the closest window in every two consecutive frames. This was done by checking the Euclidean distance between the centers of two sampling boxes. The leftmost window was used as the starting point, and the other boxes were aligned consecutively in a clockwise fashion. ImaEdge continued searching for the closest window in the next frame and repeated the alignment until all windows were accurately assigned to all frames (Fig. 4D,H; Movies 2 and 3). The mean fluorescence intensity in each sampling window can then be calculated to represent the local intensity of proteins of interest at the cortex.
Recasting the spatiotemporal dynamics of the cortex into a rectilinear coordinate allows for direct comparison across multiple cells (Machacek et al., 2009). The spatiotemporal dynamics of cortical proteins can be represented in a 2D heat map, where positions of each sampling window are plotted against their respective frame numbers (Fig. 4I,J). For example, the entry value in the rth row and cth column of the 2D heat map is the mean intensity of the rth sampling box in the cth frame. To approximate the continuous form of the spatiotemporal dynamics of cortical proteins, the heat maps can also be smoothened with a 2D Gaussian filter with customizable standard deviations. The average intensity within a group of windows can also be calculated to represent the changes in the intensity of the cortical protein (Fig. 4K).
Sensitivity to various signal-to-noise levels
Live imaging of long biological processes is often constrained by the fusion probe expression levels, the brightness and photostability of fluorophores, and the phototoxic susceptibility of the cells of interest. These factors often limit the achievable image quality. To determine the viable signal-to-noise level for robust quantitative analysis, we next evaluated the sensitivity of ImaEdge to image noise. We used the cell contours and the detected cortex from the merged images of GFP::PAR-2 and mCherry::PAR-6 channels from a C. elegans zygote obtained under typical experimental conditions as the ground-truth information for error measurement. We determined the signal-to-noise ratio (SNR) to be ∼17 dB from this test movie (Fig. 5A). Subsequently, we synthesized noise-added images at varying levels (SNR from 17 dB to 10 dB) by introducing Gaussian noise to the ground-truth movie (Fig. 5B; Fig. S2). We then assessed the accuracy of cell edge segmentation during the adaptive radial thresholding. Although strong noise intensity can be clearly observed in low-SNR images, the overlay of identified cell contours from synthesized noisy images did not show large deviations from the original cell contours. The errors, calculated as the piecewise minimum distances between segmented and ground-truth contour pixels along the ground-truth cell contour (Fig. 5C,D), were generally within 1 µm. The average root-mean-square (RMS) error of a segmented cell contour at a specific noise level can also be visualized as an error map (Fig. 5E). These analyses indicate that segmentation by ImaEdge is relatively resistant to noise over a large range of noise levels (SNR from 12 dB to 17 dB).
However, we noticed small fluctuations in the very low-SNR images (SNR from 10 dB to 11 dB). This was probably due to the intrinsic inconsistency in quality over the frames of the ground-truth movies. To understand the general trend of how image noise influences our segmentation accuracy, we determined the average RMS error over all frames (Fig. 5F). As expected, the segmentation error increased with higher noise levels. However, small standard deviations for movies with SNR from 13 dB to 17 dB were observed, indicating that the inconsistencies in movie quality do have minor influences on accurate cell segmentations.
We next investigated how the quantification of protein spatiotemporal dynamics is affected by image noise. Errors were estimated using the differences in the centroid positions of the sampling windows under various noise levels. This metric therefore takes into account both the errors in the cell segmentation and cortex detection steps. As expected, the decrease in SNR translates to an increase in sampling error, with larger variations over all noise levels compared with the previous analysis of segmentation sensitivity (Fig. 5G,H). Nevertheless, the overall RMS error of cortex sampling remained at a low level, typically within 1 µm. Taken together, our sensitivity analysis highlights the robustness of ImaEdge in processing and analyzing movies of developing C. elegans zygotes over a wide range of noise levels.
Robust cortical segmentation and measurement of cortical proteins in C. elegans zygotes
We applied ImaEdge to analyze the dynamics of mCherry::PAR-6 and GFP::PAR-2 at the cortex of C. elegans zygotes (Fig. 6A–C). The contours of the cell edge were determined from images that overlaid both GFP and mCherry channels. The defined cortical region was divided into 100 sampling windows along the entire cortex (Fig. 6A–C). Consistent with previous observations (Cuenca et al., 2003; Goehring et al., 2011), the 2D heat map depicts mutual exclusivity of mCherry::PAR-6 and GFP::PAR-2 throughout the first mitotic division (Fig. 6A–C). At the first ‘establishment’ phase after symmetry breaking, mCherry::PAR-6 was excluded from the posterior polar cortex, allowing GFP::PAR-2 to be loaded onto it. The size of the mCherry::PAR-6 and GFP::PAR-2 cortical domains gradually regressed and expanded, respectively. ImaEdge consistently captured the progressive increase in GFP::PAR-2 intensity within the posterior cortical domain until the initial over-expansion of the growing cortical GFP::PAR-2 domain was stalled and then retracted (Fig. 6A–C). In contrast, cortical intensity of mCherry::PAR-6 was generally constant during the establishment phase (Fig. 6A–C). During the following ‘maintenance’ phase, the size of the two cortical domains reached steady state (Fig. 6A–C). These analyses demonstrated that ImaEdge is not only capable of tracing the cortex, but is also able to integrate the spatiotemporal dynamics of cortical polarity proteins during polarization.
We further assessed the capability of ImaEdge in monitoring the dynamics of mCherry::PAR-6 and GFP::PAR-2 in zygotes defective in cell polarization (Fig. 6D–I). SPD-5 is a scaffold of the pericentriolar material (PCM) and is essential for centrosome maturation during establishment phase in zygotes (Hamill et al., 2002; Woodruff et al., 2015). Depletion of SPD-5 by RNA-mediated interference (RNAi) caused loss or significant delay in segregation of cortical PAR proteins (Cowan and Hyman, 2004; Hamill et al., 2002). Analysis by ImaEdge with default parameters optimized for wild-type zygotes reflected significant delay in the exclusion of mCherry::PAR-6 and loading of GFP::PAR-2 at the posterior cortex (Fig. 6D–F). The expansion of the posterior domain in the spd-5(RNAi) zygote was slower than that in the control, resulting in a smaller size of the posterior domain at the steady state (Fig. 6D–F). CDC-42 is a Rho-type small GTPase that controls the complex of PAR-6 and PKC-3 during cell polarization in zygotes (Gotta et al., 2001; Kay and Hunter, 2001). CDC-42 is required to maintain PAR-6 at the anterior domain during the maintenance phase (Aceto et al., 2006; Motegi and Sugimoto, 2006) and is essential to exclude PAR-2 from the cortex (Rodriguez et al., 2017; Schonegg and Hyman, 2006). Analysis of cdc-42(RNAi) zygotes by ImaEdge revealed uniform cortical distribution of GFP::PAR-2 during the establishment phase, while mCherry::PAR-6 was progressively segregated into the anterior cortical domain (Fig. 6G-I). Four out of five cdc-42(RNAi) embryos eventually excluded GFP::PAR-2 from the anterior cortical domain where mCherry::PAR-6 became enriched (Fig. 6G–I). Thereafter, mCherry::PAR-6 disappeared from the cortex during the maintenance phase (Fig. 6G–I). Failure in cortical exclusion of PAR-2 during the establishment phase supports a model in which CDC-42 is essential to activate PKC-3 in the anterior cortical domain (Rodriguez et al., 2017). Cortical exclusion of PAR-2 from the anterior cortex during the transition from establishment to maintenance phases suggests the presence of CDC-42-independent mechanism(s) for activation of PKC-3 at the cortex. Based on these observations, ImaEdge is suitable for the spatiotemporal analysis of cortical polarity proteins, not only in wild-type zygotes but also in a wide variety of zygotes that are defective in the process of cell polarization.
Sensitivity to fluorescence levels, discontinuous cortical signals and excessive shape changes in C. elegans zygotes
To test the sensitivity of ImaEdge under different experimental conditions, we analyzed various levels of GFP::PAR-2 in zygotes (Fig. 7A–G and Movies 4–6). The levels of gfp::par-2 transgene expression were modified via the codon-usage optimization method (Redemann et al., 2011). Synthetic transgenes with a higher usage of adapted codons results in faster translation in the germline, thus increasing the level of GFP::PAR-2 protein in early embryos. We imaged zygotes from transgenic animals expressing three synthetic gfp::par-2 transgenes with codon adaptation index (CAI) values of 0.6, 0.4 and 0.26 (Fig. 7A). ImaEdge was able to define the cell edge contour and to extract the spatiotemporal dynamics of cortical GFP::PAR-2 even using the images from embryos expressing lower levels of GFP::PAR-2 (Fig. 7A–G). Consistent with a previous report (Goehring et al., 2011), the boundary between the anterior and the posterior domain was shifted toward the anterior pole in zygotes expressing higher levels of GFP::PAR-2, while conversely, the boundary was shifted toward the posterior pole in zygotes expressing lower levels of GFP::PAR-2 (Fig. 7A–D). Analysis by ImaEdge revealed that the lower levels of GFP::PAR-2 significantly retarded the expansion of the posterior cortical domain, while the higher level of GFP::PAR-2 caused efficient expansion of the posterior cortical domain (Fig. 7A–D). An increase in the level of GFP::PAR-2 changed the kinetics of cortical loading of GFP::PAR-2 within the posterior cortex (Fig. 7B–G) and blocked the stalling-and-retraction behavior of the posterior cortical domain (Fig. 7B–D). These observations and analyses by ImaEdge support a model in which PAR-2 stimulates the expansion of the posterior cortical domain.
We next assessed the capability of ImaEdge in quantifying discontinuous signals of cortical non-muscle myosin-II heavy chain (NMY-2; Fig. 7H–J and Movie 7). Cortical NMY-2 appears as numerous foci, which undergo dynamic turnover and move toward the anterior cortex during polarity establishment (Munro et al., 2004) (Fig. 7H). Despite the discontinuous intensity of cortical NMY-2::mKate, ImaEdge achieved robust segmentation of the cell edge contour just by using the images from the NMY-2::mKate channel (Fig. 7H–J). The analysis of cortical NMY-2::mKate revealed that NMY-2 foci dynamically appeared and disappeared during the establishment and the maintenance phases, respectively (Fig. 7J). The timing of the disappearance of NMY-2::mKate occurred simultaneously with the retraction of the posterior GFP::PAR-2 domain (Fig. 7I,J). These analyses demonstrated that ImaEdge can be effectively and reliably applied to the analysis of a wide variety of cortical proteins in polarizing zygotes.
We further tested the ability of ImaEdge to analyze the cortex in zygotes that exhibit excessive cell shape changes. Depletion of the RhoGAP RGA-3 causes hyper-contraction of cortical actomyosin because of a failure in the inactivation of the small GTPase RHO-1 (Schmutz et al., 2007; Schonegg et al., 2007). Although ImaEdge was able to consistently track both cortical PAR-6::GFP and cortical NMY-2::mKate in zygotes showing excessive ruffling, it was unable to recognize the deeper cortical invaginations, where cortical PAR-6 and NMY-2::mKate signals penetrated into the cytoplasm (Fig. S3). Similarly, ImaEdge was not able to recognize cortical signals at cleavage furrows and subsequent cell–cell contact sites (see Fig. 8).
Monitoring polarization of PAR proteins in C. elegans two-cell embryos
We next applied ImaEdge to analyze the polarization of PAR proteins during the second cell division by extracting cortical information along the surface of two-cell stage embryos. In two-cell stage embryos, the anterior cell (AB) is destined to divide symmetrically, resulting in equal inheritance of PAR proteins by two daughter cells, ABa and ABp (Fig. 8A and Movie 8). In contrast, the posterior cell (P1) is engaged to undergo asymmetric cell division, resulting in asymmetric inheritance of PAR-6 and PAR-2 by two daughter cells, EMS and P2, respectively (Fig. 8A). Analysis by ImaEdge was able to automatically detect dynamic processes in the segregation of mCherry::PAR-6 to the anterovertical cortex, and that of GFP::PAR-2 to the posterior polar cortex in P1 cells (Fig. 8B–E). In contrast to zygotes, P1 cells initially showed GFP::PAR-2 localizing throughout the cortex (Fig. 8A,B). During prophase, GFP::PAR-2 was excluded from the anterovertical cortex, where mCherry::PAR-6 became enriched (Fig. 8A–C). The size and intensity of the posteriorly localized cortical GFP::PAR-2 domain was stable until prometaphase, when it started to polarize towards one side of the lateral cortex (Fig. 8B,D). Subsequently, the GFP::PAR-2 domain was inherited by the posterior-positioned daughter cell, P2.
We next tested whether ImaEdge was sensitive enough to detect the abnormal behavior of GFP::PAR-2 in P1 cells. We performed an RNAi-based screen to look for genes that, when depleted, do not affect PAR-2 polarization in zygotes, but whose absence would compromise PAR-2 behavior in P1 cells. Among the candidates tested, we found that the depletion of zyg-1 fulfilled the above requirements. ZYG-1 is a member of the Polo-like kinase 4 family and is essential for centriole replication. P1 cells depleted of ZYG-1 failed to divide, because of the absence of bipolar spindles during the second cell division. Zygotes treated with zyg-1(RNAi) exhibited cortical PAR domains with normal intensities and size during the first cell division (Fig. 8F and Movie 9). Analysis by ImaEdge showed that the distribution of GFP::PAR-2 at the posterior cortex was maintained in the P1 cell until prometaphase, when it started to disappear from the posterior cortex or remained localized throughout the entire P1 cell cortex (Fig. 8G–J). Given that P1 cells polarize in a cell-autonomous manner (independent of the anterior cell, AB) (Priess and Thomson, 1987), these observations suggest that the formation of a bipolar mitotic spindle is essential for the establishment of a polarized PAR-2 domain in P1 cells. These analyses demonstrated the versatility of ImaEdge in quantifying cortical dynamics of polarity proteins not only at the one-cell stage, but also in multicellular-stage embryos.
DISCUSSION
In this study, we have introduced an image analysis platform, ImaEdge, which allows for efficient quantitative analysis of cortical proteins in living embryos. ImaEdge enables the automatic detection of cell edge contours and accurate extraction of the cell cortex from time-lapse movies. Profiles depicting the fluorescence intensities of cortical proteins can be sampled using a finite number of boxes, allowing the spatiotemporal dynamics of cortical proteins to be represented as 2D heat maps. As pioneered earlier (Machacek et al., 2009), the heat map of cortical protein signals can also be averaged for statistical analysis. The radial adaptive segmentation method we propose here is capable of robust extraction of the entire cell cortex, even those with inhomogeneous distribution of cortical proteins in C. elegans embryos. ImaEdge comes with a set of default parameters for automatic segmentation of the cortex that can analyze both wild-type and polarity-defective mutant embryos. The current version of ImaEdge is unable to handle multiple cells. However, this step can be tackled in a relatively straightforward way by detecting multiple ROIs and analyzing them separately. In the meantime, although ImaEdge currently fails to differentiate between cells in multicellular embryos, it is still able to extract the cortical information along the surface of the embryos.
As described in Table S3, users may choose to adjust any parameters to suit their systems of interest, such as the number of sector sub-regions, the size of median filter, the number of the sampling windows and whether to use radial adaptive segmentation or simple segmentation (faster for simple datasets). In terms of the ease of further development, ImaEdge adopts simple algorithms that can be readily understood and modified by researchers in the fields of cell biology and developmental biology. Image segmentation of both the cell boundary and the cortex by ImaEdge requires only rough ranges of a couple of parameters (Table S3), because the program can scan these parameters and define the most suitable parameter values for each frame in time-lapse movies. Although the optimal parameters may be variable for different frames, the same parameter scopes can be generally used to process all frames.
ImaEdge is currently optimized for time-series images of C. elegans embryos, whose motility and cell volume changes are generally negligible (except during cytokinesis). Therefore, the constant radial sub-division mask can be determined from the first frame and then applied for all frames. Our method of adaptive radial segmentation can be readily extended to cases involving cell migration and extensive cell shape/volume changes by invoking a radially partitioned mask adaptively constructed for the cell image of each frame. Based on the extracted contours and their corresponding distance maps, protein quantities at the cortex can be sampled in a similar fashion and aligned to multiple references, such as center of mass, nuclear position or protruding cell edges.
The segmentation process used by ImaEdge is dependent on the local contrast between the cell and the background, as calculated using Otsu's method. Thus, there are currently two limitations in detection of the highly-convoluted cortex and the multicellular case where the cortex from neighboring cells are in close apposition. Ultimately, these limitations stem from the resolution limit of the imaging system. In the first case, regions of the highly convoluted cortex with high radius of curvature are not well resolved, and thus these errors inevitably propagate into any subsequent processing. In the second case, the difficulty stems from resolving one cortex from another, since the cell-to-cell gap is well below the diffraction limit. These limitations are clearly demonstrated in examples from Fig. 8 whereby ImaEdge can reveal the distinct behaviors of PAR proteins in two-cell-stage embryos between wild-type and zyg-1(RNAi) conditions, but is unable to recognize their distribution at the contact between AB and P1 blastomeres. Similarly, ImaEdge performs well on wild-type embryos, but is unable to detect the deeper cortical invaginations in hyper-contractile rga-3(RNAi) embryos. Possible solutions to these problems may require super-resolution microscopy approaches (Bertocchi et al., 2013). With current diffraction-limited microscopy, we therefore recommend users to carefully examine cortical regions defined by ImaEdge, when deeper cortical invaginations and furrows are prominent in the target cell(s).
The current version of ImaEdge is designed to analyze two-dimensional time-series images. It might be possible to extend an analogous approach to analyze higher-dimensional datasets, such as 4D images (three-dimensional time-series images). The 3D ROI can be radially divided into a finite number of cone-shaped sub-regions for adaptive segmentation. The 3D cell cortex determined using a 3D image gradient can then be sampled by connecting a set of equally spaced control points at the segmented cell surface to their closest points at the inner surface of cell cortex. These will inevitably entail additional complexity both in terms of the algorithms and visual display of the analysis results. However, at present a more acute challenge lies in the acquisition of such datasets with sufficient duration and spatiotemporal resolution while preserving sample viability. Light-sheet microscopy may address this need but this is not yet widely available (Chen et al., 2014).
ImaEdge is fully open-source and equipped with a user-friendly graphical interface, thereby providing a ready-made platform for the cell and developmental biology community to view, modify and customize their analysis as desired [ImaEdge source code and user manual can be found on GitHub (http://github.com/KanchanawongLab)]. For example, ImaEdge can be integrated with existing analyses, such as cell-marking methods (Murray et al., 2006; Ooi et al., 2006), imaging of cell-lineage tracing (Bao et al., 2006) or other analytical software (Boyle et al., 2006; Du et al., 2014; Santella et al., 2010, 2014). This will allow for the integrated and rigorous quantification of cortical protein dynamics in specific target cells over longer periods of embryogenesis such as those observed through in toto imaging (Keller et al., 2008; McMahon et al., 2008; Udan et al., 2014). We anticipate that ImaEdge will be instrumental in bringing about a deep quantitative understanding of the dynamic process of cell polarization, not only in C. elegans embryos, but also in a wide variety of polarized cells. ImaEdge thus provides a tool to link quantitative molecular information at the cortex to theoretical models, which marks a significant step forward for quantitative and predictive interpretation of cell polarization.
MATERIALS AND METHODS
Worm strains and live imaging
The C. elegans strains used in this study are MOT119: unc-119(ed3) III; ddIs238[gfp::par-2 re-coded, CAI=0.26]; par-2(ok1723); temIs17[mCherry::par-6], MOT121: unc-119(ed3) III; temIs7[gfp::par-2 no codon optimization, CAI=0.26]; par-2(ok1723); temIs17[mCherry::par-6], MOT118: unc-119(ed3) III; par-2(ok1723); ddIs25[gfp::par-2 re-coded, CAI, 0.41]]; temIs17[mCherry::par-6], MOT120: unc-119(ed3) III; ddIs239[gfp::par-2 re-coded, CAI=0.60]; par-2(ok1723); temIs17[mCherry::par-6], MOT295: nmy-2(cp52[nmy-2::mkate]) I; unc-119(ed3) III; pkc-3(it309[pkc3::gfp]). All strains were maintained at 25°C according to standard procedures (Brenner, 1974). RNAi experiments were performed by the feeding method. In short, L4440-based RNAi clones were transformed into E. coli HT115 cells and the transformants were inoculated into liquid LB with 500 µg/ml carbenicillin and grown at 37°C. Then, 100 µl of successfully transformed E. coli liquid culture was seeded onto nematode growth medium plates with 1 mM IPTG. The bacterial lawn was allowed to grow overnight at room temperature. Worms at the L1 stage and L4 stage were transferred to feeding plates for RNAi against zyg-1 and other genes, respectively, and incubated at 25°C.
C. elegans embryos were isolated from gravid hermaphrodites into egg salt buffer, placed onto coverslips, and inverted on slides with 20 µm monodisperse polystyrene beads (Bangs Labs) as a spacer. All embryos in this imaging condition completed their first mitotic cell cycle and underwent successful cytokinesis. Embryos were observed at 25°C with a CFI Plan Apochromat 60× NA1.4 oil-immersion lens on a Nikon Ti-E motorized upright microscope outfitted with a CSU-X1 spinning-disc confocal system (Yokogawa Electric) with LaserStack 491 and 561 solid-state diode lasers (Intelligent Imaging Innovation). Images with pixel size of 266.7 nm were acquired with a Photometrics Evolve512 camera (Photometrics) controlled by Metamorph (Intelligent Imaging Innovation) every 2 s using 250 ms exposure for GFP fusions and 250 ms for mKate and mCherry fusions at 20% power on the 491 and 561 lasers and 1×1 binning in the camera.
Image analysis
where A is the original image, Gx and Gy are gradient maps along the x- and y-axis, G refers to the map of the magnitude of image gradient, and Gaverage is the mean gradient along a contour containing n pixels. For the segmentation of cell contour and cortex, we adopt the automation function with multicore computing provided by ImaEdge. To smoothen the final heat map for each channel, we applied a Gaussian filter with standard deviations of 5 pixels (or 5 frames) along the x-axis and 1 pixel (or 1 sampling window) along the y-axis, to the original heatmap.
ImaEdge was coded in MATLAB (Release R2015a). All computations were performed on a Windows 7 workstation [Intel(R) Xeon(R) CPU E5-2640 v3 @ 2.60 GHz; RAM 192GB; 64-bit OS].
Acknowledgements
We thank Ken Kemphues (Cornell University), Bob Goldstein and Dan Dickinson (University of North Carolina), Nathan Goehring (The Francis Crick Institute) and the Caenorhabditis Genetic Center for sharing C. elegans strains, and core facilities in the MBI and TLL for their assistance.
Footnotes
Funding
We acknowledge funding support from the National Research Foundation Singapore (NRF) under the NRF Fellowship scheme to F.M. (NRF_NRFF2012-08) and P.K. (NRF-NRFF-2011-04), the NRF Competitive Research Programme to P.K. (NRF2012NRF-CRP-001-084), and the Ministry of Education - Singapore Academic Research Fund Tier 2 to P.K. (MOE-T2-1-124). F.M. was also supported by the Strategic Japan-Singapore Cooperative Research Program by the Japan Science and Technology Agency and the Singapore Agency for Science, Technology, and Research (1514324022). Z.Z. is supported by the MBI Graduate Scholarship.
Data availability
The source code can be downloaded from GitHub (http://github.com/KanchanawongLab). The compiled standalone software packages for Windows and Mac systems are also available from the corresponding authors on reasonable request.
References
Competing interests
The authors declare no competing or financial interests.