We present a freely available workflow built for rapid exploration and accurate analysis of cellular bodies in specific cell compartments in fluorescence microscopy images. This user-friendly workflow is designed on the open-source software Icy and also uses ImageJ functionalities. The pipeline is affordable without knowledge in image analysis.
The last decade has been characterized by breakthroughs in fluorescence microscopy techniques illustrated by spatial resolution improvement but also in live-cell imaging and high-throughput microscopy techniques. This led to a constant increase in the amount and complexity of the microscopy data for a single experiment. Because manual analysis of microscopy data is very time consuming, subjective, and prohibits quantitative analyses, automation of bioimage analysis is becoming almost unavoidable. We built an informatics workflow called Substructure Analyzer to fully automate signal analysis in bioimages from fluorescent microscopy. This workflow is developed on the user-friendly open-source platform Icy and is completed by functionalities from ImageJ. It includes the pre-processing of images to improve the signal to noise ratio, the individual segmentation of cells (detection of cell boundaries) and the detection/quantification of cell bodies enriched in specific cell compartments. The main advantage of this workflow is to propose complex bio-imaging functionalities to users without image analysis expertise through a user-friendly interface. Moreover, it is highly modular and adapted to several issues from the characterization of nuclear/cytoplasmic translocation to the comparative analysis of different cell bodies in different cellular sub-structures. The functionality of this workflow is illustrated through the study of the Cajal (coiled) Bodies under oxidative stress (OS) conditions. Data from fluorescence microscopy show that their integrity in human cells is impacted a few hours after the induction of OS. This effect is characterized by a decrease of coilin nucleation into characteristic Cajal Bodies, associated with a nucleoplasmic redistribution of coilin into an increased number of smaller foci. The central role of coilin in the exchange between CB components and the surrounding nucleoplasm suggests that OS induced redistribution of coilin could affect the composition and the functionality of Cajal Bodies.
Light microscopy and, more particularly, fluorescence microscopy are robust and versatile techniques commonly used in biological sciences. They give access to the precise localization of various biomolecules like proteins or RNA through their specific fluorescent labeling. The last decade has been characterized by rapid advances in microscopy and imaging technologies as evidenced by the 2014 Nobel Prize in Chemistry awarding Eric Betzig, Stefan W. Hell and William E. Moerner for the development of super-resolved fluorescence microscopy (SRFM)1. SFRM bypasses the diffraction limit of traditional optical microscopy to bring it into the nanodimension. Improvement in techniques like live-imaging or high throughput screening approaches also increases the amount and the complexity of the data to treat for each experiment. Most of the time, researchers are faced with high heterogeneous populations of cells and want to analyze phenotypes at the single-cell level.
Initially, analyses such as foci counting were performed by eye, which is preferred by some researchers since it provides full visual control over the counting process. However, manual analysis of such data is too time consuming, leads to variability between observers, and does not give access to more complex features so that computer-assisted approaches are becoming widely used and almost unavoidable2. Bioimage informatics methods substantially increase the efficiency of data analysis and are free of the unavoidable operator subjectivity and potential bias of the manual counting analysis. The increased demand in this field and the improvement of computer power led to the development of a large number of image analysis platforms. Some of them are freely available and give access to various tools to perform analysis with personal computers. A classification of open access tools has been recently established3 and presents Icy4 as a powerful software combining usability and functionality. Moreover, Icy has the advantage of communicating with ImageJ.
For users without image analysis expertise, the main obstacles are to choose the appropriate tool according to the problematic and correctly tune parameters that are often not well understood. Moreover, setup times are often long. Icy proposes a user-friendly point-and-click interface named “Protocols” to develop workflow by combining some plugins found within an exhaustive collection4. The flexible modular design and the point-and-click interface make setting up an analysis feasible for non-programmers. Here we present a workflow called Substructure Analyzer, developed in Icy’s interface, whose function is to analyze fluorescent signals in specific cellular compartments and measure different features as brightness, foci number, foci size, and spatial distribution. This workflow addresses several issues such as quantification of signal translocation, analysis of transfected cells expressing a fluorescent reporter, or analysis of foci from different cellular substructures in individual cells. It allows the simultaneous processing of multiple images, and output results are exported to a tab-delimited worksheet that can be opened in commonly used spreadsheet programs.
The Substructure Analyzer pipeline is presented in Figure 1. First, all the images contained in a specified folder are pre-processed to improve their signal to noise ratio. This step increases the efficiency of the following steps and decreases the running time. Then, the Regions of Interest (ROIs), corresponding to the image areas where the fluorescent signal should be detected, are identified and segmented. Finally, the fluorescent signal is analyzed, and results are exported into a tab-delimited worksheet.
Object segmentation (detection of boundaries) is the most challenging step in image analysis, and its efficiency determines the accuracy of the resulting cell measurements. The first objects identified in an image (called primary objects) are often nuclei from DNA-stained images (DAPI or Hoechst staining), although primary objects can also be whole cells, beads, speckles, tumors, or whatever stained objects. In most biological images, cells or nuclei touch each other or overlap causing the simple and fast algorithms to fail. To date, no universal algorithm can perform perfect segmentation of all objects, mostly because their characteristics (size, shape, or texture) modulate the efficiency of segmentation5. The segmentation tools commonly distributed with microscopy software (such as the MetaMorph Imaging Software by Molecular Devices6, or the NIS-Elements Advances Research software by Nikon7) are generally based on standard techniques such as correlation matching, thresholding, or morphological operations. Although efficient in basic systems, these overgeneralized methods rapidly present limitations when used in more challenging and specific contexts. Indeed, segmentation is highly sensitive to experimental parameters such as cell type, cell density, or biomarkers, and frequently requires repeated adjustment for a large data set. The Substructure Analyzer workflow integrates both simple and more sophisticated algorithms to propose different alternatives adapted to image complexity and user needs. It notably proposes the marker-based watershed algorithm8 for highly clustered objects. The efficiency of this segmentation method relies on the selection of individual markers on each object. These markers are manually chosen most of the time to get correct parameters for full segmentation, which is highly time consuming when users face a high number of objects. Substructure Analyzer proposes an automatic detection of these markers, providing a highly efficient segmentation process. Segmentation is, most of the time, the limiting step of image analysis and can considerably modify the processing time depending on the resolution of the image, the number of objects per image, and the level of clustering of objects. Typical pipelines require a few seconds to 5 minutes per image on a standard desktop computer. Analysis of more complex images can require a more powerful computer and some basic knowledge in image analysis.
The flexibility and functionality of this workflow are illustrated with various examples in the representative results. The advantages of this workflow are notably displayed through the study of nuclear substructures under oxidative stress (OS) conditions. OS corresponds to an imbalance of the redox homeostasis in favor of oxidants and is associated with high levels of reactive oxygen species (ROS). Since ROS act as signaling molecules, changes in their concentration and subcellular localization affect positively or negatively a myriad of pathways and networks that regulate physiological functions, including signal transduction, repair mechanisms, gene expression, cell death, and proliferation9,10. OS is thus directly involved in various pathologies (neurodegenerative and cardiovascular diseases, cancers, diabetes, etc.), but also cellular aging. Therefore, deciphering the consequences of OS on the human cell’s organization and function constitutes a crucial step in the understanding of the roles of OS in the onset and development of human pathologies. It has been established that OS regulates gene expression by modulating transcription through several transcription factors (p53, Nrf2, FOXO3A)11, but also by affecting the regulation of several co- and post-transcriptional processes such as alternative splicing (AS) of pre-RNAs12,13,14. Alternative splicing of primary coding and non-coding transcripts is an essential mechanism that increases the encoding capacity of the genomes by producing transcript isoforms. AS is performed by a huge ribonucleoprotein complex called spliceosome, containing almost 300 proteins and 5 U-rich small nuclear RNAs (UsnRNAs)15. Spliceosome assembly and AS are tightly controlled in cells and some steps of the spliceosome maturation occur within membrane-less nuclear compartments named Cajal Bodies. These nuclear substructures are characterized by the dynamic nature of their structure and their composition, which are mainly conducted by multivalent interactions of their RNA and protein components with the coilin protein. Analysis of thousands of cells with the Substructure Analyzer workflow allowed characterization of never described effects of OS on Cajal Bodies. Indeed, obtained data suggest that OS modifies the nucleation of Cajal Bodies, inducing a nucleoplasmic redistribution of the coilin protein into numerous smaller nuclear foci. Such a change of the structure of Cajal Bodies might affect the maturation of the spliceosome and participate in AS modulation by OS.
NOTE: User-friendly tutorials are available on Icy’s website http://icy.bioimageanalysis.org.
1. Download Icy and the Substructure Analyzer protocol
2. Opening the protocol
3. Interacting with the workflow on Icy
NOTE: Each block or box is numbered and has a specific rank within the workflow (Figure 2b). By clicking on this number, the closest possible position to the first is assigned to the selected block/box then the position of the other blocks/boxes is re-organized. Respect the right order of the blocks when preparing the workflow. For example, Spot Detector block needs pre-defined ROIs so that Segmentation blocks have to run before Spot Detector blocks. Do not modify the position of boxes. Do not use “.” in the image’s name.
4. Merging of the channels of an image
5. Segmentation of the regions of interest
NOTE: Substructure Analyzer integrates both simple and more sophisticated algorithms to propose different alternatives adapted to image complexity and user needs.
6. Fluorescent signal detection and analysis
7. Run the protocol
All the described analyses have been performed on a standard laptop (64-bit, quad-core processor at 2.80 GHz with 16 GB random-access memory (RAM)) working with the 64-bit version of Java. Random-access memory is an important parameter to consider, depending on the amount and the resolution of images to analyze. Using the 32-bit version of Java limits the memory to about 1300 MB, which could be unsuitable for big data analysis, whereas the 64-bit version allows increasing the memory allocated to Icy. Figure 3 reports the time needed for segmentation for different types of images and different resolutions. It confirms that high resolution substantially increases the time of primary object segmentation.
The data presented in the next paragraphs demonstrate that the Substructure Analyzer workflow can be used to solve most of the common problematics encountered in cellular and molecular biology (cell counting, foci counting and analysis, translocation analysis).
Measurement in numerous cells of the precise proportion of molecules concentrated within a given compartment or the change of localization across subcellular domains can be a very cumbersome and error-prone task, especially when it is performed manually. Figure 4 illustrates the ability of the workflow to rapidly and precisely quantify the well-described nuclear translocation of the transcription factor NFκB in response to TNFα stimulation. Images used for the analysis were graciously compiled by Ilya Ravkin (http://www.ravkin.net/SBS/Image-Library.htm) and are publicly available at the Broad Bioimage Benchmark Collection16. This set contains images of MCF7 and A549 cell lines treated with increasing concentrations of TNFα. The analysis has been performed on more than 40,000 cells and 96 images (48 images for each channel) from the MCF7 cell line (Figure 4a). The whole analysis (from the import of images to data saving) took 26 minutes. Each image corresponds to a specific TNFα concentration. The nuclear/cytoplasmic proportions of NFκB were homogeneous across all cells for a given image. For this reason, for each image, the values of all nuclear or cytoplasmic pixels were summed to calculate the nuclear and cytoplasmic proportions of NFκB, and individualization of touching nuclei/cytoplasms during their segmentation was not required. DAPI images were processed in the Segmentation C: Clustered objects with convex shapes block to segment nuclei and the green channel was used to delineate the cytoplasms with the Segmentation E: Clustered cytoplasm block (Figure 4b). The Fluorescence Analysis D: Global Translocation block was used to export the sum of nuclear and cytoplasmic pixel values. The obtained data are represented in a dose-response curve (Figure 4c) and illustrate the increase of NFκB nuclear proportion after stimulation with increasing concentrations of TNFα.
The workflow not only analyzes the global intensity signal in different sub-cellular compartments but can also be used to detect foci and extract specific information about their features. Figure 5a illustrates the detection of P-bodies in individual cells by localizing the enhancer of mRNA decapping 4 (EDC4) protein. The Segmentation A: Non clustered objects block was used to segment nuclei that present a low level of clustering, the Segmentation E: Clustered cytoplasm block to delineate the cytoplasms and the Fluorescence Analysis C: 2 channels in two compartments to detect EDC4 foci. Segmented objects are identified (Segmented nuclei_0, Segmented cytoplasms_0), and multiple information is collected into the corresponding worksheet (Segmented nuclei Information and Segmented cytoplasms Information) like the EDC4 mean intensity or the number/size of detected bodies in each ROI (Figure 5b). Before being analyzed, detected bodies can be preliminarily filtered according to their area, which can be useful to exclude objects below the resolution power of the microscope. To estimate the size of conserved objects, an additional area threshold can be introduced. Areas of all the detected bodies are reported in dedicated worksheets (Segmented nuclei_Distribution foci size, Segmented cytoplasms_Distribution foci size). In this example, we highlight the analysis of two cells (Segmented nuclei_0 and_1 and Segmented cytoplasm_0 and_1), where cytoplasmic EDC4 foci are detected. Note that even though the signal of the EDC4 protein is detected in both the nuclear and cytoplasmic compartments, only the cytoplasmic foci correspond to P-bodies. The nuclear signal most likely results from the cross-reactivity of the antibody used to perform immunofluorescence experiments. The size in pixels of each of the foci is given.
Then, as the effect of oxidative stress (OS) on Cajal Bodies was still unknown, we took advantage of the versatility of the workflow to study it during time kinetics (2 to 20 hours after OS induction with 500 µM H2O2) (Figure 6). Cajal Bodies are dynamic structures nucleating and dissolving within the nucleoplasm under the control of specific parameters like the transcription rate or cell cycle progression. Cajal Bodies were visualized by localizing their main structural and functional component, the coilin protein. Information about the number and the size of Cajal bodies were collected by studying 2300 individual cells from high-resolution images (3840X3072) containing 3 channels: DAPI (blue), 53BP1 (green) and coilin (red) (Figure 6a). Full processing took less than 1 h. Because nuclei presented convex shape and some of them were clustered, the block Segmentation C: Clustered objects with convex shapes was chosen to perform the segmentation. As OS is known to induce DNA double-strand breaks (DDSBs), the p53-binding protein 1 (53BP1), known to be an essential effector of DNA damage signaling pathways, was used as a marker to discriminate between stressed and non-stressed cells. Indeed, in the absence of DNA damages, 53BP1 is homogeneously distributed within the nucleoplasm, whereas following DNA damages, it concentrates on DDSBs, which form easily distinguishable foci in microscopy. The number and the size of both 53BP1 and coilin nuclear foci in each cell were analyzed by using the block Fluorescence Analysis B: 2 Channels in the same compartment. Firstly, the analysis of 53BP1 data helped to determine for each time point of the kinetics the best threshold to classify cells according to their stress status. Among the whole kinetics, OS induces a significant increase of the 53BP1 foci number with an 11-fold increase at 2, 4, and 8 h after the induction of OS (Figure S1). This effect is less pronounced at 6 h (6.5-fold) because of a higher initial level of DSB in non-stressed cells. Even after 20 h, a sustained increase remains (almost 6-fold). These data reflect the efficiency of the treatment to induce the OS and establish 53BP1 as an efficient stress-marker for both early and late observations in microscopy. On the other hand, a small proportion of stressed-cells with a low level of DDSBs was observed at each time point of the kinetics, suggesting that cells are not homogeneously stressed, enforcing the importance of using a stress marker in stress experiments. Based on the ROC analysis, an overall threshold of 17 53BP1 foci was selected to discriminate between stressed and unstressed cells throughout the kinetics (Figure S1).
Then the number and size of nuclear foci of coilin were analyzed according to the cell’s stress status. A significant increase in the number of foci was observed 2, 4, and 6 h after stress induction (Figure 6b). The most sustained effect is observed at 2 h, the number of cells having more than 10 foci increasing from 0% in the unstressed cells to 75% in the stressed cells. Moreover, more than 25% of stressed cells of that time point had more than 20 coilin foci. This effect progressively decreases until 8 h. At 20 h, a small increase of the first and the third quartiles is observed, but data also show that 84% and 80% of cells present less than 8 foci in non-stressed and stressed cells, respectively. These data show that OS also has late effects on parameters controlling the nucleation of Cajal Bodies, which might be different from early effects. Interestingly, careful analysis of the features of the coilin nuclear foci revealed that the observed increase of the number of coilin foci associates with a decrease in their size (Figure 6c). For instance, 2 h after OS induction, the proportion of coilin foci with an area below 0.2 µm2 increases from 26% to 64%. This effect decreases until 8 h, but, unlike Cajal Bodies’ number, no significant change is observed at 20 h. This might reflect that early and late nucleoplasmic redistribution of Cajal Bodies are different events.
Altogether these data strongly suggest that OS changes the nucleation power of Cajal Bodies, inducing a nucleoplasmic redistribution into numerous smaller nuclear foci. Since nucleation of Cajal Bodies is driven by their protein and RNA components, this suggests that the OS effect might change the composition of Cajal Bodies. Within Cajal Bodies, the coilin protein interacts with several different protein partners such as components of the SMN complex and Sm proteins. These interactions often involved phosphodomains of coilin, and one can imagine that a modification of the structure of Cajal Bodies could be linked to changes in the phosphorylation status of coilin. Moreover, recent works demonstrated that Cajal Bodies are not randomly localized within the nucleus, and can affect gene expression through proximal association with specific gene loci17. Given all these facts, it would not be surprising if an effect on Cajal bodies' functionality would accompany changes in their structure. From these results, we can also ask if such Cajal Bodies’ remodeling is a passive consequence of OS or participates in the response by interacting with specific gene loci, for example.
To test if a change in coilin expression could alter its localization and change the nucleation of Cajal Bodies, we overexpressed an exogenous GFP-coilin fusion protein. Nuclei were not clustered, so segmentation was performed with the block Segmentation A: Non clustered objects. Then, to precisely quantify the number and size of the coilin foci (red signal, channel 1) according to the level of GFP-coilin (green signal, channel 2) overexpression, the block Fluorescence Analysis B: 2 Channels in the same compartment was used. The level of GFP-coilin was reflected by the mean intensity of the GFP signal in individual nucleus (Figure 7a). In cells with medium or high level of GFP-coilin overexpression (GFP intensity higher than 10), the number of Cajal Bodies per nucleus significantly increases with some nuclei displaying up to 21 foci (Figure 7b). We also noticed that the brightness (mean intensity) of Cajal Bodies considerably increases in parallel with the overexpression level leading to intensity saturation. Such saturated regions might thus mask the presence of smaller and fainter foci. In cells overexpressing a medium or high level of GFP-coilin, the size of Cajal Bodies also significantly increases, the median values rising from 1 to 2 µm2 (Figure 7c). Moreover, more than 25% of cells with a high GFP-coilin expression level present foci with an area larger than 2.5 µm2, whereas it never exceeds 0.8 µm2 in non-transfected cells. In conclusion, GFP-Coilin overexpression leads to a significant increase in both the number and size of Cajal Bodies. Since OS increases the number of Cajal Bodies but reduces their size, these data might reflect that OS effect on the structure of Cajal Bodies is most probably induced by a change of their composition rather than by an effect on the cellular amount of coilin.
Figure 1: The Substructure Analyzer pipeline
The workflow is developed in Icy, an open community platform for bioimage informatics, and performs automatic analysis of multiple fluorescence microscopy images. Firstly, multi-channel images are automatically loaded within the protocol and pre-processed to improve, if necessary, the signal to noise ratio and to remove imaging artifacts. Then, image segmentation isolates regions of interest (ROIs) from the background. Several methods of segmentation are available depending on the level of clustering and the nature of the objects of interest. Segmented objects are saved with a specific descriptor (for example, Image name_Nucleus_1) in a specific folder and can be used/re-used for subsequent analysis. Fluorescent signals like spots are then analyzed within the ROIs and multiple features (location, size, shape, intensity, texture, spot number, and size) are exported into an automatically created spreadsheet. All the measured features, like the number of detected spots, are reported to the descriptor of the corresponding ROI. Please click here to view a larger version of this figure.
Figure 2: Graphical interface of the Icy workflow
(a) The workflow is composed of 13 general blocks, which can be grouped according to their general function: Select Folder (block 1) allow access to images; Merge Channels (block 2) is used to merge different channels of images. Blocks 3 to 8 are dedicated to object segmentation, each one being adapted for a specific context: Segmentation A: Non clustered objects, Segmentation B: Poorly clustered objects, Segmentation C: Clustered objects with convex shapes, Segmentation D: Clustered objects with irregular shapes (block 7, works in association with block 6), Segmentation E: Clustered cytoplasms. Blocks 9 to 11 are used to count/analyze foci within subcellular compartments: Fluorescence Analysis A: 1 Channel, Fluorescence Analysis B: 2 Channels in the same compartment and Fluorescence Analysis C: 2 Channels in two compartments. Blocks 12 and 13 are adapted for the analysis of translocation events: Fluorescence Analysis D: Global Translocation and Fluorescence Analysis E: Individual Cell Translocation. (b) Global presentation of a block architecture. Please click here to view a larger version of this figure.
Figure 3: Time needed for segmentation as a function of image resolution and object clustering level
(a) Graphical illustration of the time needed to perform segmentation depending on image resolution and clustering level. Two different segmentation blocks have been tested: Segmentation A (for non-clustered objects) and Segmentation D (for clustered objects with irregular shapes). The time of processing (in seconds per image) is plotted as a function of image resolution (number of pixels). (b) Example of images analyzed with either Segmentation A or Segmentation D block. DAPI staining and segmented nuclei are indicated. Please click here to view a larger version of this figure.
Figure 4: Workflow analysis of an NFκB nuclear translocation assay in response to TNFα stimulation
(a) Example of images of human MCF7 cells (human breast adenocarcinoma cell line) treated with TNFα from the BBBC014 image set (Broad Bioimage Benchmark Collection)12. Localization of the transcription factor NFκB has been studied in cells treated with 12 increasing concentrations of TNFα (from 0 to 10-7g/mL). For each concentration, 4 replicates were done. DAPI images (blue channel) were processed in the Segmentation C: Clustered objects with convex shapes block to segment nuclei, and then NFκB staining with FITC (green channel) was used to delineate the cytoplasms with the Segmentation E: Clustered cytoplasm block. The Fluorescence Analysis D: Global Translocation block was used to export the sum of nuclear and cytoplasmic pixel values. (b) Segmented nuclei (grey) and cytoplasms (white) are used to determine the total intensity of the NFκB fluorescent signal in each compartment. (c) The dose-response curve obtained from data analysis illustrates the increase of NFκB nuclear proportion as TNFα concentration increases. Images are in an 8-bit BMP format, and the image size is 1360 x 1024 pixels. For each concentration, the standard deviation between the 4 replicates is calculated and illustrated as error bars. Please click here to view a larger version of this figure.
Figure 5: Quantification of EDC4 foci (P-bodies) within the cytoplasm in individual cells
(a) Co-staining with DAPI (blue channel) and Phalloidin (red channel) allows the segmentation and identification of cell nuclei and F-actin, respectively. The EDC4 protein (green channel) is used as a marker for the detection of P-bodies known to be exclusively cytoplasmic. Scale bar, 10 µm. (b) Foci are counted, and several features (here their area in pixels) are exported to the Results spreadsheet structured in several sheets. Two sheets are dedicated to the analysis of nuclei and two to the analysis of cytoplasms. Each raw embeds the information of a specific ROI. In Nucleus/Cytoplasm Information sheets, the characteristics of segmented objects are exported (ROI ID, size, and mean intensity), followed by the number of detected foci within each ROI. If necessary, a size threshold (100 pixels in this example) can be added, and the number of foci below and above this threshold is reported in the two last columns. In Nucleus/Cytoplasm_Distribution of foci size sheets, the area (in pixels) of all the detected foci is reported in the raw of the corresponding ROI. In this example, we highlight the analysis of two cells (Segmented cytoplasm_0 and_1), where cytoplasmic EDC4 foci have been detected. The size in pixels of each of the foci is also given. Please click here to view a larger version of this figure.
Figure 6: Kinetics of coilin and 53BP1 nucleoplasmic distribution after oxidative stress
The analysis of the number and size of foci was based on more than 110 individual cells for each time point from 3 independent kinetics. (a) Immunofluorescence detection of coilin and 53BP1 foci in HeLa S3 cells throughout treatment kinetics with 500 µM H2O2. 53BP1 concentrates on DNA Double-Strand Break sites and is used to assess the efficiency of the treatment on each cell. Coilin is enriched in Cajal Bodies. Cells were fixed 2, 4, 6, 8, and 20 h after oxidative stress induction and co-stained with anti-coilin (red channel) and anti-53BP1 (green channel) antibodies. Scale bar, 10 µm. (b) The number of coilin nuclear foci after oxidative stress. Results are shown as box plots, where the central mark is the median, the lower and upper edges of the box are the 25th and 75th percentiles. In cells treated with H2O2, a significant increase inthe number of coilin foci is detected and is maximal after 2 and 4 h of treatment. Wilcoxon – Mann Whitney test analysis demonstrates a significant difference between untreated and H2O2 treated cells (*, p<0.05; **, p<0.01; ***, p<0.001). (c) The proportion of cells in function of the coilin foci area (µm2). Most foci have a smaller area at 2 h and 4 h time points as compared to the control. Please click here to view a larger version of this figure.
Figure 7: Study of Cajal Bodies’ nucleation in cells overexpressing a GFP-coilin fusion protein
(a) HeLa S3 cells were transiently transfected in biological triplicates with 500 ng of a plasmid expressing a GFP-Coilin fusion protein using a standard procedure following the manufacturer’s recommendations. After 48 h, cells were fixed and stained with an antibody against coilin. DAPI (blue channel) permits to segment nuclei. For these analyses, more than 100 cells were analyzed. GFP signal reflects ectopic coilin-GFP protein (green channel) localization, whereas the coilin signal (red channel) corresponds to both endogenous and exogenous protein localization. For each nucleus, the mean intensity of the GFP signal reflects the level of coilin expression. According to this parameter, the features of Coilin foci have been analyzed using the block Fluorescence Analysis B: 2 Channels in the same compartment. Scale bar, 10 µm. (b) Relationship between the number of coilin foci and mean GFP intensity per nucleus. Results are shown as box plots, where the central mark is the median, the lower and upper edge of the box are the 25th and 75th percentiles. Wilcoxon – Mann Whitney test analysis demonstrates a significant increase of the number of foci for a GFP intensity >10 (*, p<0.05; **, p<0.01; ***, p<0.001). (c) Relationship between the area of coilin foci and mean GFP intensity per nucleus. Results are shown as box plots, where the central mark is the median, the lower and upper edge of the box are the 25th and 75th percentiles. Wilcoxon – Mann Whitney test analysis demonstrates a significant increase of the area of coilin foci for a GFP intensity >10 (*, p<0.05; **, p<0.01; ***, p<0.001). Please click here to view a larger version of this figure.
Figure S1: Kinetics of 53BP1 nucleoplasmic distribution after oxidative stress
(a) Quantification of 53BP1 nuclear foci in untreated or H2O2 treated cells after different incubation times. Results are shown as box plots, where the central mark is the median, the lower and upper edge of the box are the 25th and 75th percentiles. Wilcoxon – Mann Whitney test analysis demonstrates a significant difference between untreated and H2O2 treated cells (*, p<0.05; **, p<0.01; ***, p<0.001). (b) For each time point of the kinetics, a ROC (Receiver Operating Characteristic) curve has been established to determine the best threshold to discriminate between stressed and unstressed cells. The parameters of each test are reported in the table (Sensitivity, specificity, TP: true positive, TN: true negative, FP: false positive, FN: false negative). From these data, a global threshold of 17 53BP1 foci was determined to discriminate stressed from non-stressed cells through the kinetics. Please click here to download this file.
An increasing number of free software tools are available for the analysis of fluorescence cell images. Users must correctly choose the adequate software according to the complexity of their problematic, to their knowledge in image processing, and to the time they want to spend in their analysis. Icy, CellProfiler, or ImageJ/Fiji are powerful tools combining both usability and functionality3. Icy is a stand-alone tool that presents a clear graphical user interface (GUI), and notably its “Protocols” point and click interface through which workflows can be easily designed or manipulated4. The functionality of this software is enhanced by supporting and utilizing ImageJ plug-ins, and a lot of documentation is available thanks to its large community of users. We used this strong combination of usability and functionality of the Icy software to develop the Substructure Analyzer workflow. Its main objective is to propose an automated solution to perform a full analysis of cell fluorescent signal from multiple images. The analysis encompasses image pre-processing, object segmentation, and fluorescent signal analysis.
Several protocols are gratefully shared on the Icy website and propose to use Icy functionalities to detect and analyze fluorescent signals like cellular foci. However, some of these protocols only process one image per run, do not export the results in a specific file or analyze a single channel. Others do not contain preliminary segmentation to determine regions of interest (ROIs) or propose simple algorithms like “HK-means” that are not suited for the segmentation of highly clustered objects. Substructure Analyzer automates all the steps of the image analysis from image pre-processing to the segmentation of objects and the detection/analysis of fluorescent signals. The protocol also proposes simple or more complex methods for object segmentation and export of the data (names of the segmented objects, foci analysis) is realized in optimized outputs adapted to the problematic. CellProfiler also proposes powerful pipelines composed of modules that encapsulate functionalities18,19. The available pipelines are well adapted to specific problems. Icy workflows are a network of different boxes, each of them performing a specific task. Through the “Protocols” interface, boxes composing the network are intuitive and easy to manipulate. Substructure Analyzer is adapted to several contexts like the simple merging of channels, the quantification of fluorescent signals distribution between two subcellular compartments, the analysis of foci in one or several channels from one or two cellular compartments in individual cells. Several displays allow visualizing the intermediate results during each run to control the processing. Moreover, Icy presents icon bars that group the methods by topic and from which functionalities found in the workflow can be separately manipulated to manually test specific parameters on some images before using them on a whole image set.
Like in the bioimage analysis field, the most critical step of this protocol is the object segmentation. The protocol proposes the segmentation of primary objects, which are most of the time cell nuclei, and also contains another block which task is to segment a second type of objects within the cell. Nuclei segmentation can be difficult when the objects are highly clustered so that they touch or overlap together. Different alternatives for the segmentation of primary objects are contained in the protocol, according to the shape of the objects and the level of clustering even if to date no universal algorithm success in the full segmentation of highly clustered objects. The segmentation process is mainly limited by the determination of correct parameters that could successfully segment a part of the objects but lead to under- or over-segmentation of another part. The user would have to perform different runs with different parameters to get a final correct segmentation, especially for clustered objects presenting heterogeneous and non-convex shapes.
This protocol is mostly limited by the segmentation step, which could be highly time consuming and need powerful computer configuration for the most complex cases. More specifically, the more numerous or complex the images are (clustered objects to segment, a high number of pixels in the images), the more random-access memory (RAM) the user will require. Globally, for any problematic, the user will need to run the protocol on a 64 bits Java Runtime Environment (freely available on Java website) and to use a 64 bits operating system (OS) to increase the available memory used by Icy (memory is limited to 1300 MB for 32 bits JRE). Because of the script, the protocol does not work correctly on a Mac computer. Moreover, the number of parameters increases with the complexity of the problem and will require more knowledge in image analysis. Concerning the functionality criteria, this protocol has not yet been adapted for the analysis of stacked images (time stacks, z-stacks) even if Icy efficiently performs analysis on this type of data.
Future updates will be realized notably to improve the segmentation of complex clustered objects. Additional blocks will be developed to adapt the segmentation to clustered cellular structures with irregular non-convex shapes. We will also adapt the workflow for the analysis of time and z-stacks (for live imaging and 3D data).
The authors have nothing to disclose.
G.H. was supported by a graduate fellowship from the Ministère Délégué à la Recherche et aux Technologies. L.H. was supported by a graduate fellowship from the Institut de Cancérologie de Lorraine (ICL), whereas Q.T. was supported by a public grant overseen by the French National Research Agency (ANR) as part of the second “Investissements d’Avenir” program FIGHT-HF (reference: ANR-15-RHU4570004). This work was funded by CNRS and Université de Lorraine (UMR 7365).
16% Formaldehyde solution (w/v) methanol free | Thermo Fisher Scientific | 28908 | to fix the cells |
Alexa Fluor 488 of goat anti-rabbit | Thermo Fisher Scientific | A-11008 | fluorescent secondary antibody |
Alexa Fluor 555 of goat anti-mouse | Thermo Fisher Scientific | A-21425 | fluorescent secondary antibody |
Alexa Fluor 555 Phalloidin | Thermo Fisher Scientific | A34055 | fluorescent secondary antibody |
Bovine serum albumin standard (BSA) | euromedex | 04-100-812-E | |
DMEM | Sigma-Aldrich | D5796-500ml | cell culture medium |
Duolink In Situ Mounting Medium with DAPI | Sigma-Aldrich | DUO82040-5ML | mounting medium |
Human: HeLa S3 cells | IGBMC, Strasbourg, France | cell line used to perform the experiments | |
Hydrogen peroxide solution 30% (H2O2) | Sigma-Aldrich | H1009-100ml | used as a stressing agent |
Lipofectamine 2000 Reagent | Thermo Fisher Scientific | 11668-019 | transfection reagent |
Mouse monoclonal anti-coilin | abcam | ab11822 | Coilin-specific antibody |
Nikon Optiphot-2 fluorescence microscope | Nikon | epifluoresecence microscope | |
Opti-MEM I Reduced Serum Medium | Thermo Fisher Scientific | 31985062 | transfection medium |
PBS pH 7.4 (10x) | gibco | 70011-036 | to wash the cells |
Rabbit polyclonal anti-53BP1 | Thermo Fisher Scientific | PA1-16565 | 53BP1-specific antibody |
Rabbit polyclonal anti-EDC4 | Sigma-Aldrich | SAB4200114 | EDC4-specific antibody |
Triton X-100 | Roth | 6683 | to permeabilize the cells |