Open-Source Tools for Volume Estimation of 3D Multicellular Aggregates

Featured Application: The main goal of this work is to provide an overview of open-source tools available for researchers interested Abstract: The volume is one of the most relevant features that deﬁne the treatment of an in vivo tumour. When using cancer 3D in vitro models in pre-clinical studies, it becomes important to evaluate the macroscopic e ﬀ ects of drugs and radiotherapy treatments. Depending on the nature of the 3D in vitro model used, di ﬀ erent open-source solutions can be used for measuring the volume by starting from microscope-acquired images. In this work, we introduced several open-source tools today available for estimating the volume of 3D multicellular aggregates (e.g., spheroids, organoids), also giving hints for deﬁning the “best software” by analysing characteristics of 3D in vitro models and limits of the tools. Finally, using several cancer organoids imaged by a ﬂuorescent microscope, we compared volume estimations obtained with di ﬀ erent tools, besides presenting a new version of the Reconstruction and Visualization from Multiple Sections (ReViMS version 2.0) tool. This work aims to be the reference for researchers interested in estimating the volume of 3D multicellular aggregates through an open-source tool.


Introduction
Commonly, the treatment planning of cancer patients begins with localisation of primary and secondary tumours and quantification of their volume. Depending on these and many other mass(es) parameters, cancer treatments are based on different chemotherapy, radiotherapy, immunotherapy, and chirurgical surgery procedures. Therefore, volume is one of the most relevant features for the treatment planning of an in vivo tumour. Symmetrically, it becomes important to evaluate the effects of drugs and radiotherapy treatments when using 3D in vitro models, among which multicellular The next sections are organised as follows: Section 2 presents a short overview of the tools considered in the experiments, and then detailed in Section 3. Finally, Section 4 summarises the main findings of the work.

Materials and Methods
In this section, a brief description is given for material and each software used in the experiments. Figure 1 shows the Graphical User Interface (GUI) of the different considered software.

ReViSP Version 2.1
ReViSP is an open-source software for estimating the volume of spheroids starting from the binary mask of a single projection [5]. ReViSP is written in MATLAB (The MathWorks, Inc., Natick, MA). Source code and standalone software (for Windows) are publicly available for download at: http://sourceforge.net/p/revisp. ReViSP allows reconstruction of the external surface of a spheroid and estimation of its volume by just loading the binary mask of a single projection. No parameters are required. The binary mask can be easily obtained by manually segmenting a single projection with ImageJ [13]. Otherwise, AnaSP can be used to automatically segment and analyse the original spheroid's image [9]. External surface and volume are then automatically computed by supposing a local circular symmetry of the spheroid around its skeleton. On such base, ReViSP turns out to be an extremely user-friendly software. However, it is worth noting that today several tools similar to ReViSP exist for volume and surface estimation of spheroids [3].

ReViMS Version 2.0
ReViMS is an open-source, user-friendly software for estimating the volume of cell aggregates with irregular shape and organoids (i.e., multicellular aggregates that self-organise into ex vivo mini-organs [2]) by automatically segmenting single slices of z-stacks of fluorescence images [12]. ReViMS is written in MATLAB and has been designed for biologists: it allows reconstruction of the external surface of a cell aggregate and estimation of its volume without any programming skills. Source code and standalone versions (for Windows and Linux) are publicly available for download at: http://sourceforge.net/p/revims. In the first software release (i.e., ReViMS version 1.0), for each slice, the border of the aggregate was defined with a standard histogram-based segmentation method (e.g., Otsu segmentation [14], Triangle segmentation [15], or a combination of them). In order to improve the software applicability, we implemented two new Active Contour segmentation methods that permit handling even blurry fluorescent objects.

Active Contour Segmentation
Cell objects analysed in fluorescence typically generate images with blurry contours. Consequently, the definition of the cell borders becomes more subjective. Histogram-based segmentation methods provide a rule for defining foreground limits. However, slight variations of the threshold value are liable to generate strongly different segmentations outcomes. Spatial-based segmentation methods are in general more robust in defining the border of blurry objects.
In order to extend the segmentation opportunities provided by ReViMS with methods with better performances in case of blurry objects, we included two Active Contour methods: an edge-based segmentation method, well-known in literature as Geodesic Active Contour [16], and a region-based energy method, typically called Chan-Vese model [17]. Both methods belong to the class of "snakes" methods.
The classical "snakes" approaches work by deforming an initial contour towards the boundary of the object to be detected [18]. The deformation is obtained by iteratively minimising a functional designed so that its local minimum is obtained at the boundary of the object. Basically, to stop the curve evolution the classical "snakes" models rely on an edge-function depending on the image gradient.
The Geodesic Active Contour is a connection between classical energy-based deformable contours and geometric curve evolution ones. In particular, a new term is introduced to the curve evolution models to attract the deforming curve to the boundary, thus improving the detection of boundaries with large differences in their gradient. The result is an Active Contour approach independent of object geometry and topology. The Chan-Vese model is a further extension of Geodesic Active Contour: whilst classical "snakes" models can detect only objects with edges defined by gradient, in the Chan-Vese model, the stopping process is based on a different metric [19], not related to the gradient of the image. Accordingly, it can detect even contours without gradient, for instance, objects with smooth or discontinuous boundaries.

Averaging Different 2D/3D Segmentations
One of the major problems when assessing quality performances of segmentations tools is the lack of a robust Ground Truth for the considered dataset. Typically, manual segmentation by an expert annotator can be accounted as "Ground Truth". However, the inter-and intra-rater reliability in the segmentation procedure is a well-known problem in several fields, ranging from radiology to histology [20]. Different experts clearly differ in segmenting the same object, and the variability is empathised in case of blurry images such as those typical of fluorescence microscopy.
In order to obtain a robust Ground Truth, improving and comparing 2D/3D manual segmentations obtained by different experts becomes fundamental, but these are challenging tasks. Practically, today there is not a standard procedure for "averaging" different 2D curves and 3D surfaces, even excluding outliers. Logically, the term "averaging" is not meant in the strict analytical sense because there is not a standard definition of average between curves.
In this scenario, the expectation-maximisation algorithm for Simultaneous Truth and Performance Level Estimation (STAPLE) is an interesting offer for "averaging" 2D curves [21]. The algorithm considers a collection of binary segmentations and computes a probabilistic estimate of the true segmentation and a measure of the performance level represented by each segmentation. In particular, the probabilistic estimate of the true segmentation is formed by estimating an optimal combination of the segmentations, weighting each segmentation on the estimated performance level. Besides proving the quality of their approach on several datasets, the authors also provide an open-source version of the algorithm for several programming languages, including a MATLAB version working with 2D binary segmentations (https://bit.ly/2Su6jLJ).
In order to "average" different 2D curves and 3D surfaces created by different human annotators to obtain a more robust Ground Truth to assess quality performances of automatic segmentation tools, we have included in ReViMS version 2.0 an extended version of STAPLE working with 3D segmentations. Briefly, the 3D "averaging" segmentation approach now included in ReViMS can be summarised as follows:

1.
Calculation, for each 3D object k 1 segmented by a first expert annotator, the list of the overlapping objects k j (j = 2, . . . , n; n = number of annotators) in the segmentations created by the other experts; 2.
Selection, among the listed, of the objects k j with volume spatially overlapping more than 50% with k 1 ; 3.
Selection, among the objects k j with volume spatially overlapping more than 50%, of the 2D z-sections which overlap more than 50% with the corresponding z-section of k 1 ; 4.
"Averaging" of the segmentations exploiting the standard STAPLE algorithm, by considering one-by-one each selected z-section.
Source-code and standalone versions of ReViMS version 2.0, including the above described 3D "averaging" segmentation algorithm derived from STAPLE, are available at the website: http: //sourceforge.net/p/revims.

MITK Version 2018.04
The Medical Imaging Interaction Toolkit (MITK) is a free open-source software system for development of interactive medical image processing tools. Standalone versions (for Windows, Linux and Mac), C++ code, documentation, user manuals and video tutorials are distributed at: http://mitk.org/. MITK combines the Insight Toolkit (ITK) and the Visualization Toolkit (VTK) with an application framework [22]. Precisely, MITK is an extension of VTK and ITK: it virtually re-uses anything from VTK and ITK and implements functions that they do not support. For instance, it offers multiple, consistent and synchronised orthogonal 2D-and 3D-views of the same data using labels and colours helping the user to recognise and analyse Region of Interests (ROIs).
MITK was originally developed as a common framework for Ph.D. students in the Division of Medical and Biological Informatics (MBI) at the German Cancer Research Center, with the aim of supporting the development of leading-edge medical imaging software with a high degree of interaction [23]. MITK and the MITK applications are not certified medical products. Accordingly, they must not be used in patient care and may be used in a research setting only. However, today MITK is applied worldwide in Radiology and Nuclear Medicine research projects to segment anatomical and pathological structures from CT, MR, PET and SPECT images, and thanks to its versatility, it can be also used to segment multicellular aggregates, even at even single-cell level, imaged with confocal/multiphoton/LSFM instruments. The currently available version of MITK is the version 2018.04.
In particular, the Segmentation Plugin enables for manually and (semi-)automatically segmenting ROIs, then post-processable using several available tools, such as 2D and 3D Interpolation. The 2D Interpolation creates suggestions for a segmentation whenever there is a slice clear of segmentations (i.e., there will be no suggestion in slices with even only a single segmented pixel) and with neighbouring segmented slices (these do not need to be direct neighbours but could also be a couple of slices distant). Interpolated suggestions are displayed with a different colour than manual segmentations are, until the user accepts them as part of the segmentation. The 3D Interpolation is a pure mathematical interpolation, disregarding of image radiometric information. It is activated by default when using the Manual Segmentation tools. From the second segmented contour onwards, the surface of the segmented area will be interpolated based on the given contour information. By each further contour, the interpolation result will be improved. Accordingly, segmentations of sparse contours in arbitrarily oriented planes is the best approach to achieve an optimal interpolation result. Alternatively, MITK also provides several additional plugins for automatic (e.g., such as Otsu segmentation method [14]) and manually-adjustable thresholding, for obtaining initial segmentations.

OpenSegSPIM Version 1.1
OpenSegSPIM is an open-source and user-friendly 3D automatic quantitative analysis tool for confocal/multiphoton/LSFM data [11]. Standalone versions (for Windows and Mac), MATLAB code, documentation, user manuals and video tutorials are distributed at: http://opensegspim.weebly.com/. OpenSegSPIM assumes no prior knowledge of image processing or programming and is designed to easily segment nuclei/cells from confocal/multiphoton/LSFM image stacks and extract single-cell quantitative relevant information, such as volume, sphericity, relative distance and intensity of the segmented objects.
The software leads the user through a sequence of simple pre-processing steps to obtain single-cell segmentations. A direct visual feedback is provided for the outcome of each step, and the quality of the segmentation depends on a few input parameters easily set using the GUI: the diameter of the objects of interest (set by manually drawing a few lines on the objects) and the range of intensity values considered as foreground (set using an intuitive interactive intensity adjustment tool). OpenSegSPIM also includes: (a) a sub-cellular segmentation tool, to analyse subcellular compartments such as the cell membrane; (b) a simple post-processing tool for manually editing the segmentations obtained; and (c) an automatic batch process of time series and numerous datasets.
Despite OpenSegSPIM being originally developed as a tool for single-cell/nucleus segmentation and analysis, it can be used to estimate the volume of an entire cancer multicellular aggregate by just considering the aggregate as a single object, or by summing the volume of the single composing cells. The currently available version of OpenSegSPIM is version 1.1.

Cancer Multicellular Aggregates and LSFM
In order to assess the accuracy of volume instances estimated with ReViSP, ReViMS, MITK, and OpenSegSPIM, we used the three cancer multicellular aggregates already considered in [12]. These multicellular aggregates were generated using A549 cells (lung adenocarcinoma human cell line, American Type Culture Collection, Rockville, MD, USA). A549 cells were cultured growing in adhesion in Kaighn's modification of Ham's F-12 medium (F12K, ATCC), supplemented with 10% fetal bovine serum (FBS, EuroClone, Milan, Italy), 1% penicillin/streptomycin (GE Healthcare, Milan, Italy) and 2% ampho-tericin B (Euroclone). A549 cells were expanded and maintained as a monolayer in incubator at 5% CO 2 humidified atmosphere at 37 • C and checked periodically for mycoplasma contamination using the MycoAlertTM Mycoplasma Detection Kit (Lonza, Basel, Switzerland). Today, many systems are available for generating 3D cell cultures [1]. They can be divided into the scaffold-free and scaffold-based platforms [24]. The A549 aggregates used in this work were generated with a scaffold-free rotatory cell culture system, the RCCS-8DQ bioreactor (Synthecon Inc., Houston, TX, USA) [25]. The rotator bases were put in a humidified 37 • C, 5% CO 2 incubator, and connected to power supplies installed on the exterior side of the incubator. All activities were performed in sterile conditions. A single cell suspension of 10 6 cells/mL was placed in a 50-mL bioreactor rotating chamber at an initial speed of 12 rpm. The speed was increased over time to avoid aggregate sedimentation within the culture vessels. The culture medium was refreshed every two days. After five days, the aggregates reached an approximate diameter of 300 µm and were thus transferred (one aggregate/well) to sterile 96-well low-attachment culture plates (Corning Inc., Corning, NY, USA), with each well previously filled with 100 µl of fresh culture medium. The plates were incubated at 5% CO 2 humidified atmosphere at 37 • C. The A549 aggregates were then fixed using a glutaraldehyde solution, which provides a homogeneous and stable auto-fluorescent signal in the green range of wavelengths, useful to clearly visualise the 3D surface of the aggregates. Finally, the aggregates were imaged using a Zeiss Lightsheet Z.1 with a filter at 561 µm and a 20×/1.0 NA detection optics and two-sided 10×/0.2 NA illumination optics. The z-stacks were acquired by considering the sampling Nyquist criterion to compute the z distance between subsequent frames, and finally obtaining z-stacks of around 500, 1920 × 1920 pixel resolution, 8-bit grey level images. Precisely, for Multicellular Aggregate 1 (hereafter, MA1) 262 sections were acquired, 492 sections for MA2, and 525 for MA3. Figure 2 shows the 3D representation of MA1, MA2, and MA3. It is worth noting that MA2 is characterised by a spherical shape, practically, it is a spheroid.  Figure 2 shows the 3D representation of MA1, MA2, and MA3. It is worth noting that MA2 is characterised by a spherical shape, practically, it is a spheroid.

Results and Discussion
In order to provide guidelines for researchers that must decide a tool for estimating the volume of 3D in vitro models, we compared volume measurements obtained with ReViSP, ReViMS, MITK, and OpenSegSPIM, using cancer multicellular aggregates imaged with a LSFM. First of all, we analysed the manual segmentations provided by an expert operator segmenting two times the same images, to evaluate the Ground Truth robustness when considering fluorescence datasets. Then, in a second experiment, we showed the limitations of tools designed for spheroids when applied to multicellular aggregates volume estimation. Finally, we compared segmentations automatically obtained with ReViMS, MITK, and OpenSegSPIM, before giving hints for defining the "best software" according to the characteristics of the used 3D in vitro model.

Ground Truth Robustness When Considering Fluorescence Datasets
Objects bordering Ground Truth is necessary when analysing the accuracy of automatic segmentation methods. However, the unreliability of Ground Truth is a well-known issue when considering fluorescence datasets and approaches for "averaging" different experts' segmentations, whilst removing outliers are required to improve the Ground Truth robustness.
In these experiments, we asked an expert microscopist to manually segment slice-by-slice the images of the three considered datasets. Precisely, we asked the operator to re-segment the images after a significant amount of time (i.e., one year from the first time), so to avoid remembering the contours drawn in the first place. In this way, the segmentation was essentially done from scratch without any prior information. Figure 3 reports the segmentations created by the operator for a couple of random slides of MA1 (reported in red and yellow in Figure 3c,e). The mismatch between the segmentations is more evident in boundary slices (e.g., green plane in Figure 3a), where the fluorescent signal from out-of-focus regions is more emphasised. Quantitatively, we measured the percentage absolute difference (AD%) between the number of pixels labelled as foreground in the two different segmentation stages. Average AD% maintained around 5%, this showing good repeatability of the operator when segmenting. In particular, for MA1 AD% was 4.1%, 3.5% for MA2, and 5.3% for MA3.
In order to further improve Ground Truth robustness, we used the 3D "averaging" module included in ReViMS to combine the two segmentation datasets obtained by the operator. The new segmentation dataset was then used as reference Ground Truth for the following experiments performed in this work.

Results and Discussion
In order to provide guidelines for researchers that must decide a tool for estimating the volume of 3D in vitro models, we compared volume measurements obtained with ReViSP, ReViMS, MITK, and OpenSegSPIM, using cancer multicellular aggregates imaged with a LSFM. First of all, we analysed the manual segmentations provided by an expert operator segmenting two times the same images, to evaluate the Ground Truth robustness when considering fluorescence datasets. Then, in a second experiment, we showed the limitations of tools designed for spheroids when applied to multicellular aggregates volume estimation. Finally, we compared segmentations automatically obtained with ReViMS, MITK, and OpenSegSPIM, before giving hints for defining the "best software" according to the characteristics of the used 3D in vitro model.

Ground Truth Robustness When Considering Fluorescence Datasets
Objects bordering Ground Truth is necessary when analysing the accuracy of automatic segmentation methods. However, the unreliability of Ground Truth is a well-known issue when considering fluorescence datasets and approaches for "averaging" different experts' segmentations, whilst removing outliers are required to improve the Ground Truth robustness.
In these experiments, we asked an expert microscopist to manually segment slice-by-slice the images of the three considered datasets. Precisely, we asked the operator to re-segment the images after a significant amount of time (i.e., one year from the first time), so to avoid remembering the contours drawn in the first place. In this way, the segmentation was essentially done from scratch without any prior information. Figure 3 reports the segmentations created by the operator for a couple of random slides of MA1 (reported in red and yellow in Figure 3c,e). The mismatch between the segmentations is more evident in boundary slices (e.g., green plane in Figure 3a), where the fluorescent signal from out-of-focus regions is more emphasised. Quantitatively, we measured the percentage absolute difference (AD%) between the number of pixels labelled as foreground in the two different segmentation stages. Average AD% maintained around 5%, this showing good repeatability of the operator when segmenting. In particular, for MA1 AD% was 4.1%, 3.5% for MA2, and 5.3% for MA3.
In order to further improve Ground Truth robustness, we used the 3D "averaging" module included in ReViMS to combine the two segmentation datasets obtained by the operator. The new segmentation dataset was then used as reference Ground Truth for the following experiments performed in this work.

Tools for Estimating the Volume of Spheroids Used for Multicellular Aggregates
The tools for estimating the volume of spheroids assume an approximate spherical symmetry of the multicellular aggregates to be measured. Most of the time, maximum diameter M and its orthogonal diameter P, or the average diameter D are estimated by using ellipsoid's (Equation (1)) [26] or sphere's (Equation (2)) [27] formula to approximate the mass volume: Recently, ReViSP was introduced to estimate the volume just assuming a local circular symmetry of the spheroid instead of a global spherical symmetry. However, none of these methods allow for accurate 3D representation and reliable volume data collection in the case of multicellular aggregates, where it is not possible to infer the 3D shape of the object. For instance, Figure 3f reports the 3D surface obtained by ReViSP when analysing MA3, a multicellular aggregate with irregular shape. It results in being quite different from the original surface of MA3, strongly elongated along the z-axis (Figure 2).
In order to quantitatively compare the volume estimations of ReViSP and ReViMS, we first computed a binary version of the maximum intensity projection of MA1, MA2 and MA3 with ImageJ [13] and we reconstructed the volumes by ReViSP. Then, by computing AD%, we compared the volume instances obtained with ReViMS and ReViMS with Ground Truth. In particular, the new implemented Chan-Vese model was used for segmenting the images with ReViMS. Figure 3g shows the 3D surface obtained by ReViMS for MA3. Table 1 reports the volume measurements obtained using ReViSP and ReViMS, together with Ground Truth volume and AD%. The average AD% for ReViSP (i.e., 36,6) resulted 40 times higher than the average AD% for ReViMS (i.e., 0,8). Furthermore,

Tools for Estimating the Volume of Spheroids Used for Multicellular Aggregates
The tools for estimating the volume of spheroids assume an approximate spherical symmetry of the multicellular aggregates to be measured. Most of the time, maximum diameter M and its orthogonal diameter P, or the average diameter D are estimated by using ellipsoid's (Equation (1)) [26] or sphere's (Equation (2)) [27] formula to approximate the mass volume: Recently, ReViSP was introduced to estimate the volume just assuming a local circular symmetry of the spheroid instead of a global spherical symmetry. However, none of these methods allow for accurate 3D representation and reliable volume data collection in the case of multicellular aggregates, where it is not possible to infer the 3D shape of the object. For instance, Figure 3f reports the 3D surface obtained by ReViSP when analysing MA3, a multicellular aggregate with irregular shape. It results in being quite different from the original surface of MA3, strongly elongated along the z-axis (Figure 2).
In order to quantitatively compare the volume estimations of ReViSP and ReViMS, we first computed a binary version of the maximum intensity projection of MA1, MA2 and MA3 with ImageJ [13] and we reconstructed the volumes by ReViSP. Then, by computing AD%, we compared the volume instances obtained with ReViMS and ReViMS with Ground Truth. In particular, the new implemented Chan-Vese model was used for segmenting the images with ReViMS. Figure 3g shows the 3D surface obtained by ReViMS for MA3. Table 1 reports the volume measurements obtained using ReViSP and ReViMS, together with Ground Truth volume and AD%. The average AD% for ReViSP (i.e., 36,6) resulted 40 times higher than the average AD% for ReViMS (i.e., 0,8). Furthermore, it is worth noting that ReViSP provided acceptable performance only for MA2 (i.e., 6,0), which is a multicellular aggregate with spherical shape, practically, a spheroid [2].

Tools for Estimating the Volume of Multicellular Aggregates
In order to compare the accuracy of ReViMS with other tools for segmenting cell aggregates, we segmented MA1, MA2, and MA3 also with MITK and OpenSegSPIM. First of all, we asked the same operator that previously obtained the Ground Truth, to use MITK for manually re-segmenting five distributed slices for each orthogonal plane (e.g., axial, sagittal, coronal plane), for a total of 15 slices segmented for each considered multicellular aggregate. Then, the MITK 3D Interpolation plugin was used to compute the 3D surface and volume of MA1, MA2 and MA3. On the other hand, no prior data were required to estimate the volumes with OpenSegSPIM, except that the original image datasets of the multicellular aggregates (approximately, 1 GB each dataset) resulted in being too big to be loaded in OpenSegSPIM and a 0.1 x-y resize factor was needed. In this case, the final volume for each multicellular aggregate was obtained by summing the volume of the single cells automatically segmented by OpenSegSPIM and multiplying it for a resize factor of 10. Figure 3h,i show the final MA3 surface obtained by MITK and OpenSegSPIM, respectively. Table 2 reports the volume measurements obtained using ReViMS, MITK and OpenSegSPIM, together with AD%. The average AD% for MITK and OpenSegSPIM resulted in being 1,1 and 9,6, respectively. Thus, while MITK resulted in being perfectly comparable with ReViMS (average AD% = 0,8), the OpenSegSPIM accordance was slight smaller. It is worth noting that MITK is the only tool working also in case of images containing multiple objects, as images containing more spheroids, but to use MITK is time-consuming because it always requires initial manual segmentations of the object of interest.  Figure 4 reports a flow-chart intended to help biologists, and researchers in general, to define an appropriate open-source tool to estimate the volume of multicellular aggregates according to the characteristics of the used 3D in vitro model; Table 3 summarises the main features of the suggested tools, according to the template table proposed in [28]. If the cell aggregate has a nearly spherical shape (i.e., it is a spheroid) and there is no interest in considering the single volumes of the cells composing the aggregate, ReViSP is the suggested open-source tool, as it just requires one single projection of the spheroid (i.e., a 2D image), easily attainable with a standard brightfield, phase-contrast, or DIC microscope without need of fluorescence staining. When, in the interest of estimating the total volume of the cell aggregate also considering statistics about the volumes of the composing single cells, a 3D image is required, which is easily attainable with a confocal/multiphoton/LSFM. In this case, the tool we suggest is OpenSegSPIM or MITK. In case the used 3D in vitro model is not a spheroid (i.e., it is characterised by an irregular shape) and the researcher is only interested in the total volume of the aggregate, ReViMS is the tool we suggest using.

Defining the Volume-Estimation Tool
Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 13 used 3D in vitro model is not a spheroid (i.e., it is characterised by an irregular shape) and the researcher is only interested in the total volume of the aggregate, ReViMS is the tool we suggest using.

Conclusions
This work provides an extensive overview of open-source tools today available for researchers interested in estimating the volume of multicellular 3D-oids, besides presenting a new version of ReViMS (http://sourceforge.net/p/revims), a user-friendly software tool specifically designed to estimate the volume of organoids and, in general, multicellular aggregates with irregular shape.
Using multicellular cancer aggregates imaged by a light-sheet fluorescent microscope, we compared different open-source tools for measuring the culture volume by starting from microscope-acquired images.
First of all, we proved that even an expert operator provides slightly different segmentations if asked to segment two times the same blurry fluorescence images. Accordingly, in case of fluorescent datasets, the Ground Truth of the object's border is characterised by limited robustness and algorithms for combining different segmentations whilst removing outliers are needed. Then, depending on the geometry of the 3D in vitro model considered (i.e., spherical or irregular shape), solutions based on a single projection (e.g., ReViSP) or suites analysing different sections (e.g., ReViMS and MITK) can be used to estimate the volume with different accuracy. In particular, we proved that considering a single projection works just in case of spheroids (i.e., multicellular aggregates with a spherical shape) and very inaccurate measurements are obtained in case of general cell aggregates with irregular shape. Semi-automatic approaches like MITK generally provide good results also in case of cell aggregates, but the accuracy of the measures strongly depend on the number of initial manually segmented images. However, segmentation errors affect the accuracy of the final measurement. Eventually, software like OpenSegSPIM, which automatically outlines the single cells of multicellular aggregates, can be used to obtain the volume of the entire culture by summing the single terms. However, the fairness of estimation logically depends on the extent by which cells overlap. ReViMS represents today the best trade-off between segmentation accuracy and computational time required in case of volume estimation of multicellular aggregates characterised by an irregular shape.
All the images and segmentations considered in this work, useful for instance for testing other volume estimation tools, are freely provided at: www.filippopiccinini.it/FP/2019FP-LSFM.zip.