$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Nuclear volume quantification in tumor spheroids in 3D
Automated 3D microscopy (z-stacks) enables scientists to visualize complex multicellular systems such as organoids. To quantify morphological and fluorescence signal properties at the single-cell level, cell segmentation in 3D is often required. The following example demonstrates how Cell-ACDC can segment the nuclei in organoids containing thousands of cells from the nuclear staining channel and quantify their volume (data from Ref.9). Segmentation was performed using a custom Cellpose model that was trained and published in Ref.9. The trained Cellpose model can be used directly in Cell-ACDC by providing the weights file path in the GUI when selecting the model parameters for Cellpose v2. The segmentation was performed using the second module of Cell-ACDC to batch-process all images (Figure 1A, "Segment and track" module). Next, the result was visualized in the third module GUI (Figure 1A, "Visualize and correct" module). This module has been optimized to visualize thousands of single objects with multiple annotation options (e.g., contours, overlaid segmentation masks, or text IDs). After the measurements from the 3D masks were computed, all the available measurements were documented directly in the GUI. They are accessible by going to the top menu bar (Figure 6, "Menu bar"), selecting the Measurements menu, and then Set measurements…. The pop-up dialogue lets users obtain measurement-specific information from the info buttons and select which measurements to save. For the representative results in Figure 10, "cell_vol_fl_3D" was used, which is the volume of each object calculated by multiplying the total voxels in the object by the pixel size squared and the voxel depth. These properties are either automatically extracted from the microscopy raw file or provided by the user. Computation of the measurements from this GUI requires loading the raw images. Thus, to streamline the process and enable batch processing, the measurements can also be computed from the Utilities menu (Figure 1B, "Utilities"), by going to the sub-menu Measurements and then Compute measurements for one or more experiments. Finally, the nuclear volume distribution was plotted as a representative result (Figure 10). This analysis reveals a significant fraction of small nuclei, likely due to segmentation artifacts. They could easily be removed by filtering out small objects from the segmentation mask. At the same time, the very large nuclei could be due to merged nuclei during segmentation. It is always recommended to plot the volume distribution of objects (e.g., single cells) to identify artifacts and extract additional biological information about cell size.
Quantification of time-lapse microscopy data
With time-lapse microscopy, cellular dynamics can be directly observed at the single-cell level. Besides cell segmentation, extracting temporal dynamics requires additional analysis, including cell tracking and cell pedigree annotation. Due to the interdependence of these analysis stages, errors introduced early in the pipeline can propagate to later steps. As a result, ongoing visualization and correction of segmentation, tracking, and annotation errors are necessary. The following demonstrates that Cell-ACDC is suitable for tackling these tasks. Two datasets of two different model organisms were selected: 1) budding yeast (strain DCY001-1 from Ref.35, where the two histone H2B proteins, Htb1 and Htb2, are tagged with mCitrine), and 2) mouse embryonic stem cells (mESCs, data from Ref.22).
These two datasets highlight two annotation modes available in Cell-ACDC: asymmetric and symmetric (i.e., symmetric cytokinesis, or "normal") cell division. As the mode of division differs, the two organisms require a different tracking and annotation framework.
In asymmetric division, the mother cell forms a bud that grows and eventually separates to become a daughter cell. After division, the mother cell retains its original cell ID, and its generation number increases by one, while the daughter cell is assigned a new cell ID and its generation number is set to one. The budding phase corresponds to the S/G2/M phases of the cell cycle and is annotated as such in Cell-ACDC. These annotation choices help address typical biological questions related to the cell cycle in budding yeast.
In the case of "symmetric" division (e.g., mammalian cells), the mother cell divides into two daughter cells. The mother cell, i.e., its ID, disappears upon division, and the two daughter cells receive new IDs. Additionally, the generation number of the daughter cells is increased by one relative to the mother cell. Cell-ACDC also keeps track of the parent ID, the root ID (the original ancestor cell at the beginning of a lineage), and the sister ID.
For both annotation modes, an innovative framework was developed to correct annotation errors, where the correction is automatically propagated to all past and future relevant time points. The visualization, annotation, and correction were performed in the third module GUI (Figure 1A, "Visualize and correct" module and Figure 6). To segment and track the cells in dataset 1, the model YeaZ_v226 was applied to the phase contrast channel, while for the nuclear (histone) channel, the StarDist25 model was used. Then, using the utility Tracking and lineage > Track and/or count sub-cellular objects (Figure 1B, Utilities menu bar), Cell-ACDC assigned each nucleus the Cell ID of its corresponding cell, ensuring consistency between the tables generated from cell and nucleus masks.
For dataset 2, the DeepSea22 segmentation model was used. All three models are already available in Cell-ACDC, showcasing the advantage of integrating several segmentation models into the software.
Once segmentation and tracking errors were corrected, cell pedigrees were annotated, numerical features were computed, and downstream analysis was carried out. For dataset 1, the column TaYFP_amount_autoBkgr and the number of nuclei segmented (Figure 11A) were plotted against time. "TaYFP" is the name of the nuclear channel. "amount_autoBkgr" is a proxy for the total cellular protein amount extracted from epifluorescence images36. It is calculated as the difference between the mean fluorescence intensity in each cell mask and the background median, multiplied by the cell area (in pixels). Here, the background median is calculated from all the pixels that are not segmented as cells. As expected, the H2B amount starts increasing at bud emergence (Figure 11A-ii) and reaches a constant value before nuclear division (Figure 11A-iii). This is an important quality control in histone protein homeostasis, as the amount of histone proteins is expected to be cell cycle dependent. Additionally, plotting the number of nuclei over time confirms that histone protein amounts reach the maximum roughly around nuclear division.
For dataset 2, the cell area over time for a selected cell undergoing cell division was plotted. As expected, the cell area increases until reaching a maximum value (Figure 11B-i). Then, it decreases until cell division (Figure 11B-ii) as the cell contracts. Finally, the cycle restarts for the two daughter cells. This is another recommended analysis since checking cell size changes over the cell cycle is essential to confirm that cells are growing and dividing as expected (or not in the case of specific mutants).

Figure 1: Cell-ACDC modules. (A) Overview of the 4 main modules that can be launched from the main Cell-ACDC launcher. After segmenting, tracking, and annotating microscopy data, numerical features can be computed from either the third module ("Visualise and correct"), or from (B) the Utilities menu on the top menu bar. The "Utilities" are the routines that can run automatically on multiple datasets without user input. Besides computing the measurements, other utilities include the concatenation of multiple output tables into a single table, tracking sub-cellular objects, and image pre-processing. Cell-ACDC supports 2D, 3D (z-stack or time-lapse), and 4D data (z-stacks over time), with any number of additional channels. The output table with the numerical features can then be used for downstream analysis and biological discovery (Figure 10 and Figure 11). To this purpose, Jupyter notebooks on the Cell-ACDC GitHub page are available, which include examples of plots that can be obtained from the output table. Please click here to view a larger version of this figure.

Figure 2: Cell-ACDC decision flowchart. A flowchart outlining the module to use depending on the dataset type and analysis requirements. Please click here to view a larger version of this figure.

Figure 3: Download example data. (A) Open Welcome Guide. (B) Download example data required to replicate the protocol. Please click here to view a larger version of this figure.

Figure 4: Load data for data prep. (A) Load data into the data prep GUI. (B) Select channel to load. (C) Edit and confirm image metadata. Please click here to view a larger version of this figure.

Figure 5: Run data prep process. (A) Start the process. (B) Position ROIs (for cropping) and background ROIs. Please click here to view a larger version of this figure.

Figure 6: Third module GUI to visualize and correct results. Screenshot of the third module GUI ("Visualise and correct" in Figure 1A) with highlighted items. Note that most of the buttons on the toolbars have a tooltip (accessible by hovering the mouse cursor over the button) explaining how to use that specific function. The mode selector can be used to switch between 5 modes: "Viewer", "Segmentation and Tracking", "Cell cycle analysis" (for asymmetrically dividing cells), "Normal division: Lineage tree" (for symmetrically dividing cells, e.g., mammalian cells), and "Custom annotations". Note that a toolbar or menu bar is often present in the other GUIs (e.g., "Data pre-processing" module, Figure 1). The Edit toolbar contains all the functions that can be used to edit and correct segmentation and tracking errors (e.g., brush, eraser, edit ID, etc.). The loaded image is displayed in a two-panel view, which is helpful when different annotation options are needed (e.g., cell cycle info on the left image and IDs on the right image). The right image can also be toggled off (right-click on the image and de-select Show mirrored image). Each image panel includes a LUT slider on the side to quickly adjust intensity levels. By right-clicking on the LUT control, the user can select different color maps for the intensity images. Additionally, on the right-hand side, there is a LUT selector for the color of the segmentation labels to overlay on the intensity images (annotation option called Segm. masks). On the left-side of the annotation options for the left image, there are additional toggles to control some settings, such as auto-save, font size, etc. Please click here to view a larger version of this figure.

Figure 7: Visualize pre-processing in the main GUI. (A) Open pre-processing dialogue. (B) Pre-processing dialogue with parameters used in Protocol initialized. Please click here to view a larger version of this figure.

Figure 8: Visualize segmentation output in the main GUI. (A) Select a segmentation model. (B) Set up segmentation parameters. Please click here to view a larger version of this figure.

Figure 9: Folder structure required by Cell-ACDC. To work with Cell-ACDC, the data must be arranged in a specific folder structure. While Cell-ACDC provides a module to automatically generate this structure, it is important to understand what this structure should look like. First, all files must be inside a folder called Images. Next, they all need to start with the same name, the so-called "basename_". The minimum set of required files is a single-channel TIFF file (2D, 3D z-stack, or 3D+time) and a CSV file ending with "_metadata.csv". This file must be a table with two columns, the first column called Description and the second column called values, and it should contain at least entries for SizeT, and SizeZ for the number of frames and z-slices, respectively. If an image file does not have z-slices, SizeZ must be set to 1. The same is true for no time-lapse images, where SizeT must be 1. Being a CSV file, an entry is a single line of Description,value, separated by a comma, e.g., SizeZ,1. For multiple channels, one TIFF file per channel must be generated. The folder Images must then be placed inside a folder called "Position_1". Multiple positions are allowed, and they need to be named with a consecutive number. When loading data into any of the Cell-ACDC modules, the user can either select a specific Position folder or the entire experiment folder. Please click here to view a larger version of this figure.

Figure 10: 3D quantification of tumor organoids. Screenshot of a representative tumor organoid (data from9) loaded into the third module GUI of Cell-ACDC (left), example z-slices with red contours highlighting the segmentation masks (center), and histogram of the cell volume distribution computed from the 3D segmentation masks (right). Please click here to view a larger version of this figure.

Figure 11: Quantification of time-lapse microscopy data. (A) Screenshot of time-lapse microscopy data of budding yeast cells loaded into the third module GUI of Cell-ACDC (left), and histone H2B protein amount quantification over time in a representative cell cycle (right). The zoomed images show the example cell and its bud (white arrows) at the start of the cell cycle (i), at bud emergence (ii), and at nuclear division (iii). Data is taken from Chatzitheodoridou et al.35 (B) Screenshot of time-lapse microscopy data of mouse embryonic stem cells loaded into the third module GUI of Cell-ACDC (left) and cell area (µm2) plotted as a function of time of a representative cell undergoing cell division. The zoomed images show the example cell and its daughters at maximum cell area before division (i), division into two daughter cells (ii), and the last analyzed frame after division (iii). Data is taken from Zargari et al.22,33. Please click here to view a larger version of this figure.