Accurate Automatic Detection of Densely Distributed Cell Nuclei in 3D Space
The animal brain is the most complex information processing system in living organisms ., To elucidate how real nervous systems perform computations is one of the fundamental goals of neuroscience and systems biology ., The wiring information for neural circuits and visualization of their activity at cellular resolution are required for achieving this goal ., Advances in microscopy techniques in recent years have enabled whole-brain activity imaging of small animals at cellular resolution 1–4 ., The wiring information of all the neurons in the mouse brain can be obtained using recently developed brain-transparentization techniques 5–9 ., Detection of neurons from microscopy images is necessary for optical measurements of neuronal activity or for obtaining wiring information ., Because there are many neurons in the images , methods of automatic neuron detection , rather than manual selection of ROIs ( regions of interest ) , are required and several such methods have been proposed 10 , 11 ., Detection of cells that are distributed in three-dimensional ( 3D ) space is also important in other fields of biology such as embryonic development studies 12–17 ., In these methods , cell nuclei are often labeled by fluorescent probes and used as a marker of a cell ., To identify nuclei in such images , the basic method is blob detection , which for example consists of local peak detection followed by watershed segmentation ., If the cells are sparsely distributed , blob detection methods are powerful techniques for nucleus detection ., However , if two or more cells are close to each other , the blobs are fused , and some cells will be overlooked ., These false negatives may be trivial for the statistics of the cells but may strongly affect individual measurements such as those of neuronal activity ., Overlooking some nuclei should be avoided when subsequent analyses assume that all the cells were detected , for example , when making wiring diagram of neurons or establishing a cell lineage in embryonic development ., Therefore , correct detection of all nuclei from images without false negatives is a fundamental problem in the field of bio-image informatics ., Although many efforts have been made to develop methods that avoid such false negatives , these methods seem to insufficiently overcome the problem ., In the head region of Caenorhabditis elegans , for example , the neuronal nuclei are densely packed and existing methods produce many false negatives , as shown below ., Actually , in the studies of whole-brain activity imaging of C . elegans reported so far , the local peak detection method that can overlook many nuclei was employed 3 , 18 , or the nuclei were manually detected 19 , 20 ., Highly accurate automatic nucleus detection methods should be developed in order to improve the efficiency and accuracy of such image analysis ., Here we propose a highly accurate automatic nucleus detection method for densely distributed cell nuclei in 3D space ., The proposed method is based on newly developed clump splitting method suitable for 3D images and improves the detection of all nuclei in 3D images of neurons of nematodes ., A combination of this approach with a Gaussian mixture fitting algorithm yields highly accurate locations of densely packed nuclei and enables automatic tracking and measuring of these nuclei ., The performance of the proposed method is demonstrated by using various images of densely-packed head neurons of nematodes which was obtained by various types of microscopes ., In this study , we focused on the head neurons of the soil nematode C . elegans , which constitute the major neuronal ensemble of this animal 21 ., All the neuronal nuclei in a worm of strain JN2100 were visualized by the red fluorescent protein mCherry ., The head region of the worm was imaged by a confocal microscope , and we obtained 3D images of 12 animals ( Data 1 , Fig 1A ) ., The shape of the nuclei was roughly ellipsoidal ( Fig 1B ) ., The fluorescence intensity increased toward the centers of the nuclei ( Fig 1D ) ., The typical half-radius of the nuclei was about 1 . 10 μm ( S1 Fig ) ., The distance to the nearest neighboring nucleus was 4 . 30 ± 2 . 13 μm ( mean and standard deviation , S1 Fig ) , suggesting that the neurons are densely distributed in 3D space ., The mean fluorescence intensities differed among neurons by one order of magnitude ( S1 Fig ) , making it difficult to detect a darker nucleus near a bright nucleus ., We first applied conventional blob detection techniques to the 3D image ( Fig 1C–1E ) ., Salt-and-pepper noise and background intensities were removed from the image ., The image was smoothed to avoid over-segmentation ( Fig 1C and 1D ) ., Local intensity peaks in the preprocessed image were detected and used as seeds for 3D seeded grayscale watershed segmentation ., Each segmented region was regarded as a nucleus ( Fig 1E ) ., We found that dark nuclei in high-density regions often escaped detection ., If the dark nucleus was adjacent to a bright nucleus , the fluorescence of the bright nucleus overlapped that of the dark one , and the local intensity peak in the dark nucleus was masked ( Fig 1D ) ., As a result , the seed for the dark nucleus was lost , and the dark nucleus fused with the bright nucleus ( Fig 1E ) ., The rate of false-negative nuclei was 18 . 9% ., In contrast , our proposed method successfully detected and segmented the dark nuclei ( Fig 1F ) ., The shapes of the nuclei are roughly ellipsoidal , and the fluorescence intensity increased toward the centers of the nuclei , suggesting that the intensity of nuclei can be approximated by a mixture of trivariate Gaussian distributions ., The intensities fk of the k-th Gaussian distribution gk at voxel position, x∈R3, can be written as, fk ( x ) =πkgk ( x|μk , Σk ) =πkexp\u2061 ( −12 ( x−μk ) TΣk−1 ( x−μk ) ) ,, where μk and Σk are the mean vector and covariance matrix of gk , respectively , and πk is an intensity scaling factor ., To explain the effect on the curvature , typical bright and dark nuclei were approximated by the Gaussian distribution and are shown in Fig 2 as iso-intensity contour lines ( Fig 2A , 2C and 2E ) and plots of the intensity along the cross section ( Fig 2B , 2D and 2F ) ., When a bright nucleus was near a dark nucleus , the peak intensity of the dark nucleus merged with the tail of the fluorescence intensity distribution of the bright nucleus and no longer formed a peak ., These false negatives can be avoided by using methods for dividing a close pair of objects , or clump splitting ., Such methods have been developed for correct detection of objects in two-dimensional ( 2D ) images 22–26 ., These methods focus on the concavity of the outline of a blob ., The concavity was calculated based on one of or a combination of various measurements such as angle 25 , area 27 , curvature 26 , and distance measurements 24 of the outline ., In these methods , after binarization of the image , concavity was obtained for each point on the outline ., Then the concave points were determined as the local peaks of the concavity ., After determination of concave points , a line connecting a pair of concave points is regarded as the boundary between the objects ., When we regard the outermost contour line in Fig 2E as the outline of the fused blob ( Fig 2G ) , the conventional 2D clump splitting method can be easily applied and two concave points are detected from the fused blob ( Fig 2G , red circles ) ., The blob was divided into two parts by a border line connecting the two points , and the dark nucleus was detected ., In the ideal case in Fig 2E , we obtained necessary and sufficient number of concave points ., In real images , however , we might obtain too many concave points because outlines often contain noise and are not smooth ., However , the number of concave points to choose is unknown because it is hard to know how many nuclei are included in a blob in a real image ., Further , it is not obvious how to find the correct combinations of concave points to be linked if a blob contains three or more objects ., In addition , for 3D images , the concepts of border lines that connect two concave points cannot be naturally expanded to three dimensions , because now we need some extra processes such as connecting groups of concave points in order to form border surfaces ., Even if we regard a 3D image as a stack of 2D images , it is hard to split objects fused in the z direction ( direction of the stacks ) 11 , 27 ., Here we introduce a concept of areas of concavity instead of concave points ( i . e . local peak of concavity ) ., Hereafter we use curvature as a measure of concavity and focus on areas of negative curvature for simplicity and clarity , but other measures such as angle , area , and distance from convex hull may be applicable ., Furthermore we used the iso-intensity contour lines inside the object in addition to the outline of the object ., Near the concave points in Fig 2E , the iso-intensity contour lines have negative curvature; i . e . , they curve in the direction of low intensities ., Negative curvature may be a landmark of the border line because a single Gaussian distribution has positive curvature everywhere ., Actually , the voxels at which an iso-intensity contour line has negative curvature were between two Gaussian distributions ( Fig 2E , area between the broken lines ) ., Once these voxels are removed from the blob , detection of two nuclei should be straightforward ., This approach is different from the classic clump splitting methods in two respects; focusing on area rather than local peak of concavity ( concave points ) , and using iso-intensity contour lines in addition to the outline ., These differences eliminate the need for determining how many concave points should be chosen and for obtaining correct combinations of the concave points because the area of negative curvature will cover the border lines ., Therefore we can use the approach even if a blob contains three or more objects ., In addition , this approach is robust to noise because it does not depend on a single contour line ., Furthermore , this approach can be expanded to 3D images naturally because the 3D area ( i . e . voxels ) of negative curvatures will cover the border surfaces of the 3D objects ., Iso-intensity contour lines in 2D images are parts of iso-intensity contour surfaces in three dimensions ., A point on an iso-intensity surface has two principal curvatures , which can be calculated from the intensities of surrounding voxels ( S2 Text ) 28 ., The smaller of the two principal curvatures is positive at any point in a single Gaussian distribution but is negative around the border of two Gaussian distributions ., Therefore , once voxels that have negative curvature are removed from the blob , two or more nuclei should be detected easily in 3D images ., Thus our approach solves the above problems of the classic clump splitting methods ., We applied the above approach to real 3D images ( Fig 3 ) ., The original images were processed by denoising , background removal , and smoothing to obtain the preprocessed images ., The peak detection algorithm could find only a peak from the bright nucleus , and the blob obtained by watershed segmentation contained both nuclei ., The principal curvatures of the iso-intensity surface were calculated from the preprocessed image ., There were voxels of negative curvature in the area between two nuclei , but the area did not divide the two nuclei completely ., The voxels of negative curvature were removed from the blob , and the blob was distance-transformed; these procedures were followed by 3D watershed segmentation ., Thus , the two nuclei were separated , and the dark nucleus was successfully detected ., After voxels of negative curvature were removed from the blobs , the size of blobs obtained by the second watershed segmentation tended to be smaller than real nuclei , and the distances between the blobs tended to be larger ., To obtain the precise positions and sizes of the nuclei , least squares fitting with a Gaussian mixture was applied to the entire 3D image using a newly developed method ( see Methods ) ., The number of Gaussian distributions and the initial values of the centers of the distributions were derived from the above results ., Repeated application of watershed segmentation may increase over-segmentation ., If the distance between two fitted Gaussian distributions is too small , the two distributions may represent the same nuclei ., In this case , one of the two distributions was removed to avoid over-segmentation , and the fitting procedure was repeated with a single Gaussian distribution ., The proposed method detected 194 out of 198 nuclei in the 3D image ( Fig 4 ) ., Among the four overlooked nuclei , the intensities of two of them were too low to be detected ., The other two had moderate intensities but were adjacent to brighter nuclei ., In these cases , curvature-based clump splitting successfully split the two nuclei ., However , deviations of the brighter nuclei from Gaussian distributions disrupted the fitting of the Gaussian distributions and resulted in misplacement of the Gaussian distributions for the darker nuclei , which were instead fitted to the brighter nuclei ., On the other hand , the proposed method returned 11 false positives ., Two of them resulted from the misplacement of the Gaussian distribution for the darker nuclei described above ., Four of them were not neuronal nuclei but were fluorescence foci intrinsic to the gut ., Three of them were the result of over-segmentation of complex-shaped non-neuronal nuclei ., One of them was mislocalized fluorescence in the cytosol ., The last one was the result of over-segmentation of a large nucleus that was fitted with two Gaussian distributions separated by a distance larger than the cutoff distance ., We compared the performance of the proposed method with five previously published methods for nucleus segmentation ( Fig 5 and Table 1 ) ., Ilastik 29 is based on machine learning techniques and uses image features such as Laplacian of Gaussian ., FARSight 30 is based on graph cut techniques ., RPHC 1 was designed for multi-object tracking problems such as whole-brain activity imaging of C . elegans and uses a numerical optimization-based peak detection technique for object detection ., 3D watershed plugin in ImageJ 31 consists of local peak detection and seeded watershed ., This method is almost the same as the conventional blob detection method used in our proposed method ., CellSegmentation3D 32 uses gradient flow tracking techniques and was developed for clump splitting ., This method has been used in the study of automated nucleus detection and annotation in 3D images of adult C . elegans 33 ., We applied these six methods to 12 animals in Data 1 ( Fig 5 ) and obtained the performance indices ( Table 1 , see Methods ) ., The parameters of each method were optimized for the dataset ., The 3D images in the dataset contains 190 . 92 nuclei on average , based on manual counting ., The proposed method found 96 . 9% of the nuclei and the false negative rate was 3 . 1% , whereas the false negative rate of the other methods were 11 . 2% or more ., The false positive rate of the proposed method was 4 . 9% and that of the other methods ranged from 2 . 1% to 21 . 2% ., The proposed method shows the best performance with both of the well-established indices , F-measure 12 and Accuracy 34 , because of the very low false negative rate and modest false positive rate ., It should be noted that all of the compared methods overlooked more than 10% of nuclei in our dataset ., The reason for this was suggested by the segmentation results , in which almost all of these methods failed to detect the dark nuclei near the bright nuclei and fused them ( Fig 5 , right column ) ., These results suggest that all the compared methods have difficulty in handling 3D images with either large variance of object intensity or dense packing of objects , or both ( S1 Fig ) ., These results clearly indicate that our proposed method detects densely distributed cell nuclei in 3D space with highest accuracy ., Very low false negative rate is the most significant improvement of the proposed method from the other methods , suggesting that the proposed method will improve efficiency and accuracy of image analysis steps drastically ., Because none of the computational image analysis methods is perfect , experimenters should be able to correct any errors they find ., Therefore , a user-friendly graphical user interface ( GUI ) for visualization and correction of the results is required ., We developed a GUI called RoiEdit3D for visualizing the result of the proposed method and correcting it manually ( S2 Fig ) ., Because RoiEdit3D is based on ImageJ/Fiji 35 , 36 in MATLAB through Miji 37 , experimenters can use the familiar interface and tools of ImageJ directly ., Developers can extend the functionality using a favorite framework chosen from various options such as ImageJ macros , Java , MATLAB scripts , and C++ languages ., Interface with downstream analyses should be straightforward because the corrected results are saved in the standard MATLAB data format and can be exported to Microsoft Excel ., Three-dimensional images are shown as trihedral figures using the customized Orthogonal View plugin in ImageJ ( S2 Fig ) ., Fitted Gaussian distributions are shown as ellipsoidal regions of interest ( ROIs ) in each view ., The parameters of the Gaussian distributions are shown in the Customized ROI Manager window in tabular form ., The Customized ROI Manager and trihedral figures are linked , and selected ROIs are highlighted in both windows ., When the parameters of the distributions or the names of nuclei are changed in the Customized ROI Manager window , the corresponding ROIs in the trihedral figures are updated immediately ., Least squares fitting with a Gaussian mixture can be applied after ROIs are manually removed or added ., RoiEdit3D can be used for multi-object tracking ., The fitted Gaussian mixture at a time point is used as an initial value for the mixture at the next time point , and a fitting procedure is executed ( Fig 6A ) ., Additionally , the intensities of nuclei can be obtained as parameters of the fitted Gaussian distributions ., We tried to track and measure the fluorescence intensity of nuclei in real time-lapse 3D images ( Data2 ) ., The animal in the image expressed a calcium indicator , so neural activity during stimulation with the sensory stimulus , sodium chloride , could be measured as changes in the fluorescent intensity ., The proposed nucleus detection method was applied to the first time point in the image and found 194 nuclei out of 198 nuclei ., Seventeen false positives and four false negatives were corrected manually using RoiEdit3D ., Then the nuclei in the time-lapse 3D image were tracked by the proposed method ., Most of the nuclei were successfully tracked ., One or more tracking errors occurred in 27 nuclei during 591 frames , and the success rate was 86 . 4% , which is comparable to that in the previous work 1 ., The tracking process takes 19 . 83 sec per frame ( total 3 . 25 hr ) ., The ASER gustatory neuron was successfully identified and tracked in the time-lapse 3D image by the proposed method ( Fig 6B ) ., The ASER neuron reportedly responds to changes in the sodium chloride concentration 38 , 39 ., We identified a similar response of the ASER neuron using the proposed method ( Fig 6C ) ., This result indicates that the proposed method can be used for multi-object tracking and measuring , which is an essential function for whole-brain activity imaging ., Furthermore the proposed method was utilized to measure the fluorescence intensity of nuclei in time-lapse 2D images ( Data 3 ) ., The proposed nucleus detection method was applied to the image for the first time point ( S3 Fig ) ., Data 3 does not contain images of a highly-localized nuclear marker , and therefore the images of calcium indicator that was weakly localized to the nuclei were used instead ., The proposed method found 7 nuclei out of 9 nuclei ., Six false positives and two false negatives were corrected manually using RoiEdit3D ., Then the nuclei were tracked by the proposed method ., All of the nuclei were successfully tracked during 241 time frames ., The ASER neuron was successfully identified and tracked in the 2D images ., The response of the ASER neuron in the 2D images ( S3 Fig ) is similar to that in the 3D images ., This result indicates that the proposed method can be used for multi-object tracking and measuring of 2D images as well as 3D images ., In this article , we proposed a method that accurately detects neuronal nuclei densely distributed in 3D space ., Our GUI enables visualization and manual correction of the results of automatic detection of nuclei from 3D images as well as 2D images ., Additionally , our GUI successfully tracked and measured multiple objects in time-lapse 2D and 3D images ., Thus , the proposed method can be used as a comprehensive tool for analysis of neuronal activity , including whole-brain activity imaging ., Although the microscopy methods for whole-brain activity imaging of C . elegans have been intensively developed in recent years 3 , 18–20 , computational image analysis methods were underdeveloped ., In these works , the neuronal nuclei in the whole-brain activity imaging data were detected either manually or automatically by peak detection ., Manual detection is most reliable but time- and labor-consuming , whereas the accuracy of the automatic peak detection is relatively low because of overlooking dark nuclei near bright nuclei ., Our proposed method will reduce the difficulty and improve the accuracy ., Furthermore , the numbers of the neuronal nuclei found or tracked in these four works were less than the real number of neuronal nuclei 3 , 18–21 ., The scarcity may be due not only to the experimental limitations such as fluctuation of fluorescent protein expression or low image resolution , but also to the limitations of the image analysis methods that may overlook nuclei ., The proposed method can detect almost all the nuclei in our whole-brain activity imaging data ( Fig 6 ) , suggesting that the proposed method can avoid errors that may be caused by overlooking nuclei , such as erroneous measurements of neural activities and misidentifications of neuron classes ., Thus , our method will be highly useful for the purpose ., Peng and colleagues have intensively developed the computational methods for automatic annotation of cell nuclei in C . elegans 33 , 40 , 41 ., Although their methods successfully annotate cells in many tissue such as body wall muscles and intestine , the methods seem not to be applicable to annotations of head neurons in adult worms , which is highly desired in the field of whole-brain activity imaging 20 ., They pointed out that the positions of neuronal nuclei in adult worms are highly variable 33 and this may be one of the reasons for the difficulty ., The accuracy of detection and segmentation of neuronal nuclei may be another reasons because CellSegmentation3D that was incorporated in their latest annotation framework 33 shows compromised performance in our dataset ( Table 1 , Fig 5 ) ., Our proposed method improves the accuracy of neuronal nucleus detection and will promote developing the automatic annotation methods for the neurons ., It is noteworthy that the method of simultaneous detection and annotation of cells 41 is unique and useful in the studies of C . elegans ., Because the method assigns the positions of reference to the sample image directly and avoid the detection step , the method find cells without overlooking under some conditions , but would not work correctly under the large variation of the numbers or the relative positions of the nuclei , both of which are observed in our dataset ., The optimal method for accurate detection of nuclei will vary depending on the characteristics of the nuclei ., Many conditions such as the visualization method , shape , and distribution of nuclei will affect these characteristics ., In our case , the distributions of the fluorescence intensity of nuclei were similar to Gaussian distributions; thus , we developed an optimal method for such cases ., Even if an original image does not have these characteristics , some preprocessing steps such as applying a Gaussian smoothing filter may enable application of our method to the image ., Although choosing the optimal method and tuning its parameters might be more work than manual identification , the automatic detection method would improve subjectivity and effectivity ., In the field of biology , it is often the case that hundreds or thousands of animals should be analyzed equally well ., In such case , manual detection would be time-consuming and the automatic detection method would be required ., For tracking the nuclei in time lapse images , we can apply the detection method to each time frame separately and then link the detected nuclei between frames ., In this case , some false negatives and false positives would be separately produced for each frame , and they might disrupt the link step , resulting in increase of tracking errors ., On the other hand , in the proposed method , the result of the automatic detection could be corrected manually , resulting in decrease of tracking errors ., The proposed tracking method is a simplistic approach ., Combination with existing excellent tracking methods will likely improve tracking performance of the proposed method ., Cell division and cell death did not occur in our data , but they are fundamental problems in the analysis of embryonic development ., It may be important to improve our method if it is to be applied to these problems so that the method handles such phenomena appropriately ., C . elegans strains JN2100 and JN2101 were used in this study ., Animals were raised on nematode growth medium at 20°C ., E . coli strain OP50 was used as a food source ., We used three datasets in this study ., Data1 and 2 contain ~200 neuronal nuclei , and Data3 contains 9 nuclei ., The positions of the centers of the nuclei were manually corrected by experimental specialists using the proposed GUI ., The blobs of the nuclei were detected by the conventional method ( Steps 1 & 2 ) ., Under-segmented blobs were detected and split in Step, 3 . The precise positions and sizes of the nuclei were obtained in Step, 4 . The names and parameter values of the filters used in the proposed method are shown in S1 Table ., The performance of proposed method for cell detection was compared with five state-of-the-art methods: Ilastik , FARSight , RPHC , 3D watershed plugin in ImageJ , and CellSegmentation3D ., Ilastik is machine learning-based method and required a training data that was created manually ., The parameters of RPHC was the same as the literature 1 ., The parameters of the other methods were optimized based on F-measure and accuracy ., The parallel displacements of the raw 3D images of 12 animals in Data 1 was corrected , and the methods were applied to the images ., Because FARSight crashed during processing , its command line version ( segment_nuclei . exe ) was used 50 ., The input images for FARSight and CellSegmentation3D were converted to 8-bit images because they could not operate with 32-bit grayscale images ., For CellSegmentation3D , because it could not operate with our whole 3D image , the input images were divided and processed separately ., The comparison was performed and the processing time was measured on the same PC as that used for the proposed method ., All the methods other than CellSegmentation3D might be able to utilize multi-threading ., The centroids of the segmented regions obtained by each program were used as the representative points of the objects ., For the proposed method , the means of the fitted Gaussian distributions ( μk ) were used as the representative points ., The Euclid distances of the representative points and manually pointed Ground Truth were obtained ., If a representative point was nearest-neighbor of a point of Ground Truth and vice versa , the object was regarded as a True Positive ., If only the former condition was met , the Ground Truth was regarded as a False Negative ., If only the latter condition was met , the object was regarded as a False Positive ., We obtained the indices of the performance 12 , 34 , 50 as follows:, Truepositiverate=TPGT ,, Falsepositiverate=FPGT ,, Falsenegativerate=FNGT ,, F−measure=2×TP2×TP+FN+FP ,, Accuracy=TPTP+FN+FP ,, where, GT=TP+FN ., GT , TP , FP and FN mean Ground Truth , True Positive , False Positive and False Negative , respectively .
To measure the activity of neurons using whole-brain activity imaging , precise detection of each neuron or its nucleus is required ., In the head region of the nematode C . elegans , the neuronal cell bodies are distributed densely in three-dimensional ( 3D ) space ., However , no existing computational methods of image analysis can separate them with sufficient accuracy ., Here we propose a highly accurate segmentation method based on the curvatures of the iso-intensity surfaces ., To obtain accurate positions of nuclei , we also developed a new procedure for least squares fitting with a Gaussian mixture model ., Combining these methods enables accurate detection of densely distributed cell nuclei in a 3D space ., The proposed method was implemented as a graphical user interface program that allows visualization and correction of the results of automatic detection ., Additionally , the proposed method was applied to time-lapse 3D calcium imaging data , and most of the nuclei in the images were successfully tracked and measured .
To reach the ultimate goal of neuroscience to understanding how each neuron functions in the brain , whole-brain activity imaging techniques with single-cell resolution have been intensively developed ., There are many neurons in the whole-brain images and manual detection of the neurons is very time-consuming ., However , the neurons are often packed densely in the 3D space and existing automatic methods fail to correctly split the clumps ., In fact , in previous reports of whole-brain activity imaging of C . elegans , the number of detected neurons were less than expected ., Such scarcity may be a cause of measurement errors and misidentification of neuron classes ., Here we developed a highly accurate automatic cell detection method for densely-packed cells ., The proposed method successfully detected almost all neurons in whole-brain images of the nematode ., Our method can be used to track multi-objects and enables automatic measurements of the neuronal activities from whole-brain activity imaging data ., We also developed a visualization and correction tool that is helpful for experimenters ., Additionally , the proposed method can be a fundamental technique for other applications such as making wiring diagram of neurons or establishing a cell lineage in embryonic development ., Thus our framework supports effective and accurate bio-image analyses .
invertebrates, fluorescence imaging, engineering and technology, caenorhabditis, neuroscience, animals, animal models, caenorhabditis elegans, human factors engineering, model organisms, neuroimaging, research and analysis methods, computer and information sciences, imaging techniques, animal cells, man-computer interface, calcium imaging, graphical user interface, cellular neuroscience, cell biology, computer architecture, neurons, nematoda, biology and life sciences, cellular types, image analysis, organisms, user interfaces