The mitochondrial network is extremely complex, making it very challenging to analyze. A novel MATLAB tool analyzes live confocal imaged mitochondria in timelapse images but results in a large output volume requiring individual manual attention. To address this issue, a routine optimization was developed, allowing for speedy file analysis.
The complex mitochondrial network makes it very challenging to segment, follow, and analyze live cells. MATLAB tools allow the analysis of mitochondria in timelapse files, considerably simplifying and speeding up the process of image processing. Nonetheless, existing tools produce a large output volume, requiring individual manual attention, and basic experimental setups have an output of thousands of files, each requiring extensive and time-consuming handling.
To address these issues, a routine optimization was developed, in both MATLAB code and live-script forms, allowing for swift file analysis and significantly reducing document reading and data processing. With a speed of 100 files/min, the optimization allows an overall rapid analysis. The optimization achieves the results output by averaging frame-specific data for individual mitochondria throughout time frames, analyzing data in a defined manner, consistent with those output from existing tools. Live confocal imaging was performed using the dye tetramethylrhodamine methyl ester, and the routine optimization was validated by treating neuronal cells with retinoic acid receptor (RAR) agonists, whose effects on neuronal mitochondria are established in the literature. The results were consistent with the literature and allowed further characterization of mitochondrial network behavior in response to isoform-specific RAR modulation.
This new methodology allowed rapid and validated characterization of whole-neuron mitochondria network, but it also allows for differentiation between axon and cell body mitochondria, an essential feature to apply in the neuroscience field. Moreover, this protocol can be applied to experiments using fast-acting treatments, allowing the imaging of the same cells before and after treatments, transcending the field of neuroscience.
Cellular mitochondria sit at the center of all physiological states, and a thorough understanding of their homeostasis (mitostasis) and behavior is paramount to assist in identifying pharmacological treatment for a wide range of illnesses, including cancer and Alzheimer's disease1,2.
Mitochondria play crucial cellular roles in energy homeostasis, ATP generation, calcium buffering, and ROS regulation, and mitostasis is essential for maintaining protein homeostasis as molecular chaperones are energy-dependent3. These require a constant and dynamic network modulation and adaptation to efficiently meet cellular needs, and mitochondria transport is regulated by different signaling pathways; previous work has described one such pathway, that of retinoic acid receptors (RARs)4,5. Retinoic acid (RA) promotes axonal and neurite outgrowth via RAR activation. In mouse primary cortical neurons, activation of RAR-β encourages mitochondrial growth, speed, and mobility in the neurite6.
Considering the mitochondrial network adaptability and dynamics, the possibility of assessing mitostasis in "real time" is essential not only for investigating energy homeostasis, but also for proteostasis, cellular health, proliferation, or signaling. A commonly used method for evaluating mitostasis relies on confocal microscopy after highlighting mitochondria using a fluorescent dye or marker, as well as a specific microscopy setup allowing temperature and/or CO2 regulation7. This type of experimental setup entails that one experimental replicate be performed at a time. In addition to experimental repetition of different treatments, it should be considered that most experiments should have their technical replicates (where more than one position is imaged per plate), with a series of focal planes (z-stacks) being recorded in a series of time points. Thus, an experimental design with three repetitions of one control and two treatments, with five imaging positions per plate, and 15 time points, results in 225 stacks to be processed. Classically, videos of live mitochondria were analyzed by plotting kymographs, which would be individually analyzed8, in a time-consuming process requiring extensive manual input, even when relying on computer tools.
An algorithm was recently described9 that allows automated segmentation and tracking of mitochondria in live-cell 2-D and 3-D time-lapse files. Other quantification techniques are available, and all have their limitations10. Mitometer, an automated open-source application, is particularly adequate for time lapse and mitochondria dynamics analysis, requiring low user input. This application has a series of advantages over other existing MATLAB-based tools, namely allowing the automatic processing of individual TIF stacks, using up to 13 different parameters, particularly interesting for neurosciences, as it differentiates between peri- and tele-nuclear mitochondria.
However, for an experiment such as the above described, these 13 parameters applied to 225 stacks result in 2,925 individual output files. These require four individual computer inputs, which sum up to over 10,000 manual inputs being required to download all output files. For large experimental designs, this results in a needlessly extremely time-consuming analysis of each file and data integration. Herein we present a routine optimization that allows swift file analysis, greatly reducing document reading and data processing, analyzing data in a defined manner, consistent with the output from existing tools.
NOTE: This protocol has two main steps: a wet lab step, involving cell culture and live confocal microscopy to obtain images of live mitochondria (Figure 1) and an in silico step to analyze obtained images (Figure 2). For automated data analysis of 3D live imaged mitochondria, the MATLAB application Mitometer was used as provided by Lefebvre et al.9. The Routine optimization is written in MATLAB. The software, updated versions and processing ImageJ Macros are freely available online through GitHub, at https://github.com/JoseJoaoMV/Routine_Optimization_Mitometer_APP_MATLAB.
1. Live microscopy
Figure 1: Experimental protocol. SH-SY5Y cells were differentiated and treated with retinoids. (A) TMRM was used to live image healthy mitochondria in treated cells using a confocal microscope, capturing a time lapse z-stack of five visual fields. (B) Mitometer application MATLAB automatically segments and analyzes mitochondria images. In addition to analyzing, this software automatically discriminates mitochondria according to nuclear proximity. Blue dots are mitochondrial initial positions; red dots are final positions. Scale bar = 30 µm. Abbreviation: TMRM = tetramethylrhodamine, methyl ester. Please click here to view a larger version of this figure.
2. Image analysis
Figure 2: Routine optimization. (A) Representative code of the routine optimization. (B) Routine optimization Live-Script. (C) Routine optimization workflow. (D) Routine optimization result validation: representative image of mitochondria in untreated cells (left panel), treated with atRA (10-7 M, 72 h, middle panel), and treated with RAR antagonist BMS493 (10-7 M, 72 h, right panel), imaged after incubation with TMRM (20 nM, 45 min incubation). Scale bar = 30 µm. (E) TMRM Intensity in cell body mitochondria. Significant decrease with all-trans retinoic acid treatment (atRA, 10-7 M, 72 h), compared with control (p=0,0062), not observed when treated with RAR antagonist (BMS493, 10-7 M, 72 hours). Five cells were quantified from each of three repetitions per condition. Please click here to view a larger version of this figure.
To enhance and accelerate the analysis of output files in .txt format, a routine optimization was coded that reads data consistent with Mitometer .txt output files, with columns representing a frame and lines representing identified mitochondria. The routine optimization produces data in a single value per parameter by averaging the frames for each identified mitochondria and then averaging the results of all mitochondria per visual field. The developed routine reads files from folders numbered from 1 upwards. The Live Script Routine Optimization outputs a single spreadsheet file in a table format. In this output, columns represent analyzed parameters (e.g., "Major Axis Length"; "Intensity") and lines represent visual fields of analytic variables (e.g., "Control").
Previously published results describe mitochondrial membrane potential in the cell body of primary neuronal cultures to decrease and axonal mitochondria movement to increase after RARβ activation6.
Similar treatments were performed in differentiated neuroblastoma SH-SY5Y cells, treated with the retinoic acid receptor agonists and antagonists for 72 h (Figure 3). Collected data were plotted and analyzed using nondirectional, 2-tailed, 2-sample, equal variance Student's t-test. Comparisons were made, in the output spreadsheet file, between appropriate groups, with α = 0.05.
Figure 3: Regulation of mitochondrial homeostasis by isoform-specific retinoid signaling. (A) Representative image of mitochondria after treatment (top), respective surface plot (middle), and automated segmentation (bottom). Blue dots are mitochondrial initial positions; red dots are final positions. Scale bars = 30 µm. (B) Heatmap summarizing all variations found to mitochondria parameters; Significant variance was found between the cell body and neurite (two-way ANOVA, p=0.0158), with significant influence of RAR isoform-specific modulation in all mitochondria (two-way ANOVA, p=0.0082). (C) Mitochondrial Length – a significant decrease was observed in cells treated with AM580 (p=0.0179). (D) Mitochondrial surface area – significant decreases were observed in cells treated with AM580 (p=0.000406) and with BMS493 (p=3.01 × 10-8). (E) TMRM intensity, normalized for mitochondrial volume – significant decreases were found in cells treated with RA (p=0.00621) and Ch55 (p=0.000542). * p < 0,05; ** p < 0,01; # p < 0,0001. Five cells were quantified from each of three repetitions per condition. Please click here to view a larger version of this figure.
Analysis reveals that treatment of differentiated SH-SY5Y cells with RARα agonist AM580 results in a decrease in average mitochondria length and surface area; this effect was not found when treating with agonists for other isoforms other than exclusively RARα, but there was an interesting decrease of 35.42 ± 0.5% in the mitochondrial surface area following treatment with RAR pan-antagonist BMS493 (p = 3.01 x 10-8). In contrast, retinoids appear to have an opposite effect in terms of TMRM intensity, which relates to mitochondria membrane polarization11: while treatment with RARα agonist appears to not have a significant effect in TMRM intensity, treatment with RAR pan-agonist atRA results in a dramatic decrease of 54.82 ± 18.01% (p=0.00621). This decrease is also found following treatment with RARα/β agonist CH55 (28.99 ± 4.97% decrease, p=0.000542), and possibly also following treatment with RARβ specific agonist CD2314 (37.01 ± 28.96% decrease, p=0.09134). Importantly, this method differentiated between axonal mitochondria and the ones in the cell body, allowing the study of isoform-specific RAR stimulus and mitochondrial modulation.
Figure 4: Minimal photobleaching throughout imaging protocol. (A) Z-stack timelapses were processed in FIJI, and z-project projections of average intensities were exported for all time points. (B) A z-axis profile of the same processing was plotted according to time. (C) Photobleaching quantification for experimental setup. No significant changes were observed (p = 0.7607; paired t-test for first and last frame mean signal intensity). Please click here to view a larger version of this figure.
Live cell imaging produces large files that require serious computing processing, but even the most recent tools require extensive manual input to process. This routine optimization is focused on simplifying the process of mitochondria analysis on the Mitometer because this tool presents a very good balance between user input and data output. A comprehensive comparison between different tools for mitochondria image analysis has previously been reviewed10. While other pipelines are more focused on analyzing mitochondria networks and cluster mass or analyzing membrane potential variation, this new methodology here presented allows rapid and validated characterization of whole-cell mitochondria network, also allowing for differentiation between axon and cell body mitochondria, an essential feature to apply in the neuroscience field.
The MATLAB application Mitometer9 analyses mitochondria in image series: diffuse background is subtracted from every time-frame and z-plane in the series, which are then convolved with a Gaussian kernel to eliminate high-frequency noise, followed by an intensity threshold resulting in a mask of the segmented mitochondria. The ideal settings maximize the median number of mitochondria identified while minimizing the fluctuation in mitochondria number and area across adjacent temporal frames of the picture, with a mitochondrion being allocated to its track from the previous frame for translational motion.
Differentiated SH-SY5Y cells were used as a neuronal model for experimental validation of the routine optimization. This human cell line is a homogeneous neuroblast-like cell line expressing several neuronal-like features, such as enzyme activity, receptors, or neurofilaments, that easily proliferates in culture. This model allows the experimentation of human-derived cells without the associated ethical concerns and with much lower costs than using primary cultures12. Undifferentiated SH-SY5Y are proliferative, with short processes; differentiation of these cells with sequential retinoic acid and BDNF treatment stops proliferation and promotes neurite elongation, providing a useful in vitro neuronal model13.
All-trans retinoic acid (10-7 M) was used as a pan-RAR agonist; AM580 (10-7 M) is a RARα agonist; CD2314 (10-7 M) is a RARβ agonist; Ch55 (10-7 M) is a RARα and RARβ agonist; BMS493 (10-7 M) is a RAR pan-antagonist and was used as a pharmacological control. This methodology was validated using a similar model that has recently been described: the activation of RAR in neuronal cells regulates mitochondria homeostasis in the neuron6. Similarly, results obtained using this optimized routine are consistent with the literature, showing significant alteration in several mitochondrial parameters (mitochondria length, surface area, and mitochondria membrane potential), discriminating mitochondria in the neurite from those in the cell body (Figure 3). Mitometer can automatically identify and separate mitochondria according to distance from the nucleus (Figure 1). This is a feature from the original application9 and was used as-is.
This optimization allowed quick processing of imaging files, significantly reducing operator input, document reading, data processing (to a speed of 100 files per minute), reducing bias, and increasing experimental blinding. With this experimental setup, each of the frames takes approximately 1 min to complete a cycle; shorter intervals can be acquired by diminishing the number of z-planes captured (possibly accompanied by increasing pinhole aperture) or the number of visual fields; longer intervals can be accomplished by selecting a delay period before captured.
An experimental setup of five positions with 8 z-planes and 15 frames takes approximately 19 min per cell culture dish, resulting in minimal photobleaching (Figure 4). The main experimental limitation with mitochondria live imaging is finding the balance between enough confluence for cells to be healthy and sparsity to allow discriminating between cells, and, particularly, neurites. If cells are plated too confluently in the glass-bottom dishes, their proliferation during the differentiation process results in overlapping cells and the formation of neurite networks, difficulty in imaging individual cells and neurites and identifying transport direction, confounding mitochondria network mapping. In terms of processing, MATLAB does not have the same computational power as other languages, such as Python, but MATLAB is particularly good at signal processing, making it ideal for mitochondria imaging experiments.
This protocol also allows imaging before and after acute treatment with a liquid solution. To do so, the microscope incubation chamber must be carefully opened to give access to the glass-bottom dish for treatment application. This works best with larger volumes (>100 µL), as the turbulence treatment application can distribute it throughout the preparation, while smaller volumes would require agitation of the preparation, possibly shifting its orientation relative to the objective/plate holder. Should no such drift occur, the positions recorded in the software would relate to the same cells imaged before treatment, and a new imaging loop can be initiated. However, it should be considered that this variation would increase photobleaching, and compromises must be made in initial laser intensity. Moreover, this same protocol can be applied to other cellular models, such as neuronal primary cultures6 or even other types of cells14, but it would require optimization of the microscopy section if imaging more mobile cells or for longer imaging periods.
The authors have nothing to disclose.
Image acquisition was performed in the LiM facility of iBiMED, a node of PPBI (Portuguese Platform of BioImaging): POCI-01-0145-FEDER-022122. This work was supported by FCT (EXPL/BTM-SAL/0902/2021) LCF (CI21-00276), a grant to DT from the Fundação para a Ciência e Tecnologia of the Ministério da Educação e Ciência (2020.02006.CEECIND), a grant from ATG-The Gabba Alumni Association to VP, and the Institute for Biomedicine-iBiMED, University of Aveiro.
AM580 | Sigma-Aldrich | A8843 | |
BDNF | Thermo-Fisher | RP8642 | |
BMS493 | Tocris Bioscience | 3409 | |
CD2314 | Tocris Bioscience | 3824 | |
Ch55 | Tocris Bioscience | 2020 | |
Foetal Bovine Serum | Thermo-Fisher | 10270106 | |
GraphPad Prism v4.0 | GraphPad Software, La Jolla | n/a | |
Ham’s F12 Nutrient Mix | Thermo-Fisher | 21765029 | |
MATLAB R2022a | MathWorks | n/a | |
Minimal Essential Medium | Thermo-Fisher | 31095 | |
Nunc Glass Bottom Dishes | Thermo-Fisher | 150680 | |
Phosphate Buffer Saline Solution | Thermo-Fisher | 28372 | |
Retinoic acid | Sigma-Aldrich | R2625 | |
TMRM | Thermo-Fisher | T668 | |
Zeiss LSM 510 | Carl Zeiss | n/a | Equiped with live-cell imaging culture chamber and 63x oil immersion objective |