$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
The protocol described above has been performed (Figure 1) on 20-month-old male UM-HET3 mice. The mice were administered either vehicle (Control) or a chemical reprogramming cocktail (Treated) using model 2004 osmotic minipumps53. Livers were then sectioned and imaged on an Airyscan2 confocal microscope before processing the z-stacks using the analysis pipeline.
The 3D fluorescence images were processed and analyzed using ZEISS arivis Pro Software 4.1.2. Before segmentation of mitochondria, the particle enhancement algorithm was applied to denoise the images and remove background, out-of-focus light, autofluorescence, and/or fluorescence originating from non-specific binding. The particle enhancement selectively amplifies the signal from a given diameter. After denoising, individual mitochondria were segmented in 3D using the Blob Finder algorithm, which detects objects based on local intensity maxima within a defined size range. Finally, object segmentation reduces each 2D slice to a binary image.
To illustrate the segmentation workflow, Figure 2 presents representative optical slices at each stage of the processing pipeline for Control and Treated samples. The raw fluorescence image shows the unprocessed Tom20 signal, which contains a mixture of specific mitochondrial fluorescence and unwanted background, out-of-focus light, and autofluorescence. After applying particle enhancement denoising, the background signal was substantially reduced, and the remaining fluorescence is confined to discrete, well-defined structures (Figure 2). This results in a cleaner representation of mitochondrial boundaries that facilitates accurate segmentation.
The output of the Blob Finder segmentation is shown in green. To confirm that the z-stacks are accurately segmented, the denoised 2D slices with the associated segmented binary images were overlaid (Figure 2). If the segmentation process is sufficiently optimized, strong and consistent co-localization of mitochondrial fluorescence with the binary objects should be observed. If the segmentation process is over-segmenting (i.e., splitting individual mitochondria into multiple objects) or under-segmenting (i.e., combining multiple mitochondria into one), adjust the segmentation parameters. Visual inspection of the insets confirms high-quality segmentation: detected objects co-localize precisely with the denoised Tom20 fluorescence signal, with no evidence of over-segmentation, object oversizing, inappropriate merging of spatially adjacent mitochondria into single objects, or generation of spurious objects in regions without any Tom20 signal. This level of segmentation fidelity was consistently achieved across both Control and Treated samples, demonstrating that the pipeline parameters optimized in ZEISS arivis software are robust and generalizable across the experimental groups. Together, these results validate the segmentation approach and support the reliability of the morphological measurements derived from it.
The segmented objects can be measured for mitochondrial parameters such as volume, sphericity, and surface area (Figure 3). To account for the large number of mitochondria detected per biological sample and the non-independence of individual organelles within a cell, per-sample median values were calculated and used as the unit of statistical comparison. Group differences were assessed using the Wilcoxon test. Treated samples exhibited significantly higher mitochondrial volume and surface area compared to controls, indicating that chemical reprogramming induces substantial remodeling of mitochondrial morphology.
Lastly, the segmented mitochondria can be classified into different subtypes54,55 (Figure 4). The elbow method56 was used to determine the optimal number of clusters for k-means clustering57,58 (Figure 4A). Then, using k = 4 clusters, random subsamples of mitochondria for Control and Treated groups were plotted via uniform manifold approximation and projection (UMAP) dimensionality reduction (Figure 4B), colored by both mitochondrial subtype (left) and treatment group (right). Finally, the proportions of each mitochondrial subtype (Figure 4C) for the Control and Treated groups were plotted. In the Treatment group, there were fewer fragmented mitochondria but no significant change in the proportion of networked or elongated mitochondria.

Figure 1: Overview of protocol. Mice are treated with drugs or a vehicle for up to 6 weeks using subcutaneous osmotic minipump implantation (1). Following treatment, the mouse tissues are fixed, infiltrated with sucrose, embedded in OCT Compound, and cryo-sectioned (2). The tissue sections are then stained with antibodies against Tom20 and secondary antibodies labeled with Alexa Fluor 568 (3). Super-resolution z-stacks are collected on an Airyscan2 point scanning confocal microscope (4), and the denoised processed z-stacks are segmented in arivis to quantify mitochondrial morphology (5). Please click here to view a larger version of this figure.

Figure 2: Automated pipeline for mitochondrial segmentation and morphology analysis.
Representative single optical slices from Control (upper panels) and Treated (lower panels) samples. The Tom20 immunofluorescence signal is shown in gray (left). The same image after denoising is displayed in gray, followed by the segmentation output highlighting detected mitochondrial objects in green. The merged image shows the overlap between the denoised signal and segmented objects. Insets show higher magnification views of the indicated regions, illustrating the effects of denoising and the accuracy of segmentation. Please click here to view a larger version of this figure.

Figure 3: Mitochondrial morphology in the Control and Treated groups. Box plots show the median ± interquartile range of per-sample median values (n = 10 per group). Each dot represents one biological sample, summarized as the median across all mitochondria detected within that sample. Volume and surface area are displayed on a log10 scale. Groups were compared using Wilcoxon rank-sum tests. Significance levels: ns p ≥ 0.05, **p < 0.01, ***p < 0.001. Please click here to view a larger version of this figure.

Figure 4: Unsupervised morphometric classification of mitochondrial subtypes. (A) Elbow plot showing the total within-cluster sum of squares (WSS) as a function of the number of clusters (k), calculated on a random subsample of 10,000 mitochondria. The inflection point at k = 4 was selected as the optimal number of clusters for downstream analysis. (B) (Left) UMAP dimensionality reduction of mitochondrial morphometric features (n = 20,000 randomly sampled mitochondria), colored by assigned subtype. Four morphologically distinct subtypes were identified: fragmented (small, high sphericity), elongated (long, low sphericity), networked (large, irregular), and intermediate. (Right) UMAP projection colored by experimental group (Control vs. Treated) showing the distribution of each group across the morphological space. (C) Per-sample proportions of each mitochondrial subtype in each experimental group. Each datapoint represents one biological sample (n = 10 per group). Box plots show the interquartile range and median, with whiskers extending to 1.5 times the interquartile range. Subtype proportions were compared between groups using the Wilcoxon rank-sum test, and p-values were corrected for multiple comparisons using the Benjamini-Hochberg method. ns p ≥ 0.05, ****p < 0.0001. Please click here to view a larger version of this figure.
Supplementary File 1: Example arivis pipeline. Please click here to download this file.
Supplementary File 2: R script for quantifying mitochondrial morphology and subtype populations.Please click here to download this file.