CoChat
Sign Up
Cited by 18

The DOI matches its Crossref record; 2 of the 2 references it cites were found in a registry or an index.

Lusca: FIJI (ImageJ) based tool for automated morphological analysis of cellular and subcellular structures

Šimunić, Iva, Jagečić, Denis, Isaković, Jasmina, Dobrivojević Radmilović, Marina, Mitrečić, Dinko·2024·Scientific Reports
DOI 10.1038/s41598-024-57650-6PMC PMC10978859
Open this paper in CoChat — free account

Abstract

The human body consists of diverse subcellular, cellular and supracellular structures. Neurons possess varying-sized projections that interact with different cellular structures leading to the development of highly complex morphologies. Aiming to enhance image analysis of complex biological forms including neurons using available FIJI (ImageJ) plugins, Lusca, an advanced open-source tool, was developed. Lusca utilizes machine learning for image segmentation with intensity and size thresholds. It performs particle analysis to ascertain parameters such as area/volume, quantity, and intensity, in addition to skeletonization for determining length, branching, and width. Moreover, in conjunction with colocalization measurements, it provides an extensive set of 29 morphometric parameters for both 2D and 3D analysis. This is a significant enhancement compared to other scripts that offer only 5–15 parameters. Consequently, it ensures quicker and more precise quantification by effectively eliminating noise and discerning subtle details. With three times larger execution speed, fewer false positive and negative results, and the capacity to measure various parameters, Lusca surpasses other existing open-source solutions. Its implementation of machine learning-based segmentation facilitates versatile applications for different cell types and biological structures, including mitochondria, fibres, and vessels. Lusca’s automated and precise measurement capability makes it an ideal choice for diverse biological image analyses.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41598-024-57650-6

Cites

2 of 2 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Pipeline for image analysis of neural projections and bodies. Blue squares represent image segmentation, while grey, orange, green, red, violet, and yellow represent different options for image analysis. These include colocalization, neural bodies, neural projections, length and branching, width, and area, number, and intensity, respectively. Step 1 shows the enlarged input images of neurons and nuclei stained with MAP2 (red), SMI312 (green), and DAPI (cyan). ROIs are segmented from the background with TWS in step 2, following intensity, area, and circularity thresholding to obtain the segmented image in step 3. “Neural bodies” analysis includes the calculation of area, number, and intensity for both nuclei and soma images. Neural bodies image is acquired by using the Boolean operator “AND” on previously dilated nuclei image and neuron image from step 3 to avoid false signals. The area, number, and intensity of both neural bodies and nuclei are calculated with the particle analyser after redirection to the corresponding input images (orange square, step 4). Analysis option “Neural projections” includes the calculation of area, number, intensity, length, branching, and width. Particle analysis was performed again to obtain area, number, and intensity results (green square, step 4). Length and branching are calculated after forming skeletons in step 5, while width calculation further involves the transformation of the step 5 image into a 32-bit image. Following thresholding and deduction of 255, values of NaN for the background pixels, and 0 for skeletons, are obtained. Simultaneously, using the Local Thickness on the image from step 3, accurate width dimensions are achieved (step 6). The images from step 5 and step 6 are added to obtain the image in step 7 from which width results are calculated. The separation of these options into distinct squares (red, violet, and yellow) facilitates better user comprehension and easier application, especially considering the versatility of Lusca in analysing various biological objects beyond neurons.
Figure 1 Pipeline for image analysis of neural projections and bodies. Blue squares represent image segmentation, while grey, orange, green, red, violet, and yellow represent different options for image analysis. These include colocalization, neural bodies, neural projections, length and branching, width, and area, number, and intensity, respectively. Step 1 shows the enlarged input images of neurons and nuclei stained with MAP2 (red), SMI312 (green), and DAPI (cyan). ROIs are segmented from the background with TWS in step 2, following intensity, area, and circularity thresholding to obtain the segmented image in step 3. “Neural bodies” analysis includes the calculation of area, number, and intensity for both nuclei and soma images. Neural bodies image is acquired by using the Boolean operator “AND” on previously dilated nuclei image and neuron image from step 3 to avoid false signals. The area, number, and intensity of both neural bodies and nuclei are calculated with the particle analyser after redirection to the corresponding input images (orange square, step 4). Analysis option “Neural projections” includes the calculation of area, number, intensity, length, branching, and width. Particle analysis was performed again to obtain area, number, and intensity results (green square, step 4). Length and branching are calculated after forming skeletons in step 5, while width calculation further involves the transformation of the step 5 image into a 32-bit image. Following thresholding and deduction of 255, values of NaN for the background pixels, and 0 for skeletons, are obtained. Simultaneously, using the Local Thickness on the image from step 3, accurate width dimensions are achieved (step 6). The images from step 5 and step 6 are added to obtain the image in step 7 from which width results are calculated. The separation of these options into distinct squares (red, violet, and yellow) facilitates better user comprehension and easier application, especially considering the versatility of Lusca in analysing various biological objects beyond neurons.CC BY · PubMed Central
Figure 2 Qualitative and quantitative comparison of different scripts for neuron analysis with Lusca and manual tracing on 2D high and low-quality-stained and 3D images. Comparison between open-source programs, Lusca, and manual tracing performed on MAP2 high and low-quality-stained input images. (a) Qualitative comparison of Lusca, NeurphologyJ, and Neurite Tracer final output images for each image stain quality. Datasets (neuronal bodies and nuclear count, neurite length and width, and neuron area/volume) generated by each open-source program, Lusca and manual measurements subjected to linear regression for 2D (b) high and (c) low-quality-stained images, as well as (d) 3D images. The equations for the lines of best fit and the coefficients of determination are presented in the figure key.
Figure 2 Qualitative and quantitative comparison of different scripts for neuron analysis with Lusca and manual tracing on 2D high and low-quality-stained and 3D images. Comparison between open-source programs, Lusca, and manual tracing performed on MAP2 high and low-quality-stained input images. (a) Qualitative comparison of Lusca, NeurphologyJ, and Neurite Tracer final output images for each image stain quality. Datasets (neuronal bodies and nuclear count, neurite length and width, and neuron area/volume) generated by each open-source program, Lusca and manual measurements subjected to linear regression for 2D (b) high and (c) low-quality-stained images, as well as (d) 3D images. The equations for the lines of best fit and the coefficients of determination are presented in the figure key.CC BY · PubMed Central
Figure 3 False positive and negative operational classifications and measurements rate comparison between NeurphologyJ, NeuriteTracer and Lusca. False negative signal defined as an unrecognized visible signal by the macro. False positive signal characterised by the macro but lacking an actual signal. Recognised false positive signals for neural projections: (a) background noise mistaken for neurite, while for neural bodies these include (b) thick neurite and (c) background noise misplaced for neural body. False negative signals for neural projections: (d) neurites recognised as background and (e) thick neurite misplaced to neural bodies, and for neural bodies (f) two connected neural bodies and (g) neural body recognised as background. For further comparison of ImageJ/FIJI macros, false positive and negative signals measured for each script on high and low-quality-stained MAP2 neuron images: false positive neurite length rate (h), false negative neurite length rate (i), false positive count rate of neuronal bodies (j), false negative count rate of neuronal bodies (k).
Figure 3 False positive and negative operational classifications and measurements rate comparison between NeurphologyJ, NeuriteTracer and Lusca. False negative signal defined as an unrecognized visible signal by the macro. False positive signal characterised by the macro but lacking an actual signal. Recognised false positive signals for neural projections: (a) background noise mistaken for neurite, while for neural bodies these include (b) thick neurite and (c) background noise misplaced for neural body. False negative signals for neural projections: (d) neurites recognised as background and (e) thick neurite misplaced to neural bodies, and for neural bodies (f) two connected neural bodies and (g) neural body recognised as background. For further comparison of ImageJ/FIJI macros, false positive and negative signals measured for each script on high and low-quality-stained MAP2 neuron images: false positive neurite length rate (h), false negative neurite length rate (i), false positive count rate of neuronal bodies (j), false negative count rate of neuronal bodies (k).CC BY · PubMed Central
Figure 4 Lusca pipeline applied for morphological 3D magnetic resonance angiography (MRA) stack analysis. For 3D perception the anatomical planes (a) coronal, (b) sagittal and (c) transversal are shown. Blue, yellow, red, and violet squares represent image segmentation, area, number and intensity, length and branching, and width measurements respectively. The anatomical landmarks of MRA stacks used to standardize the volume and position of the maximum intensity projection ROIs (Optional step, grey square). Vessels on the input images in step 1 segmented from the background with TWS in step 2 and, after area and intensity thresholding, the step 3 segmented image is acquired. In step 4, after redirection of segmented image to the input image, the area, number, and intensity are obtained. For length and branching measurements, in step 3 the image undergoes skeletonization and analysis. Width is calculated by applying Local thickness mask on the image from step 3 to get step 6 image. Simultaneously, the threshold is applied to the 32-bit image from step 5, and after subtracting 255 values NaN for the background and 0 for skeleton pixels are obtained. Images from steps 5 and 6 are added to get step 7 image that serves for width calculations.
Figure 4 Lusca pipeline applied for morphological 3D magnetic resonance angiography (MRA) stack analysis. For 3D perception the anatomical planes (a) coronal, (b) sagittal and (c) transversal are shown. Blue, yellow, red, and violet squares represent image segmentation, area, number and intensity, length and branching, and width measurements respectively. The anatomical landmarks of MRA stacks used to standardize the volume and position of the maximum intensity projection ROIs (Optional step, grey square). Vessels on the input images in step 1 segmented from the background with TWS in step 2 and, after area and intensity thresholding, the step 3 segmented image is acquired. In step 4, after redirection of segmented image to the input image, the area, number, and intensity are obtained. For length and branching measurements, in step 3 the image undergoes skeletonization and analysis. Width is calculated by applying Local thickness mask on the image from step 3 to get step 6 image. Simultaneously, the threshold is applied to the 32-bit image from step 5, and after subtracting 255 values NaN for the background and 0 for skeleton pixels are obtained. Images from steps 5 and 6 are added to get step 7 image that serves for width calculations.CC BY · PubMed Central
Figure 5 Lusca pipeline for morphological analysis of mitochondria. Blue, red, and yellow squares represent image segmentation, length and branching, and area, count and intensity measurements respectively. Mitochondria immunocytochemistry images in (a) normoxic and (b) after oxygen–glucose deprivation treatment were stained with Tomm20 antibody. Mitochondria (step 1) are segmented from the background with TWS in step 2 and after area, intensity and circularity thresholding, the step 3 image is acquired. In step 4 after redirecting the input image the result area, count and intensity are obtained. For further length and branching measurements, the step 3 image is skeletonized and analysed (step 5).
Figure 5 Lusca pipeline for morphological analysis of mitochondria. Blue, red, and yellow squares represent image segmentation, length and branching, and area, count and intensity measurements respectively. Mitochondria immunocytochemistry images in (a) normoxic and (b) after oxygen–glucose deprivation treatment were stained with Tomm20 antibody. Mitochondria (step 1) are segmented from the background with TWS in step 2 and after area, intensity and circularity thresholding, the step 3 image is acquired. In step 4 after redirecting the input image the result area, count and intensity are obtained. For further length and branching measurements, the step 3 image is skeletonized and analysed (step 5).CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 2 references

Research works printed but not linked: 30

  1. Titus, A. D. J. et al. Hypobaric hypoxia-induced dendritic atrophy of hippocampal neurons is associated with cognitive impairment in adult rats. Neuroscience 145, 265–278 (2007). · DOI
  2. York, E. M., LeDue, J. M., Bernier, L.-P. & MacVicar, B. A. 3DMorph automatic analysis of microglial morphology in three dimensions from ex vivo and in vivo imaging. eNeuro 5(6), ENEURO.0266–18.2018. (2018). · DOI
  3. Abramoff, M. D., Magalhaes, P. J. & Ram, S. J. Image processing with ImageJ. Biophoton. Int. 11, 36–42 (2004).
  4. Carpenter, A. E. et al. Cell Profiler: Image analysis software for identifying and quantifying cell phenotypes. Genome Biol. 7, R100 (2006). · DOI
  5. Meijering, E. et al. Design and validation of a tool for neurite tracing and analysis in fluorescence microscopy images. Cytometry. A 58, 167–176 (2004). · DOI
  6. Myatt, D. R., Hadlington, T., Ascoli, G. A. & Nasuto, S. J. Neuromantic—from semi-manual to semi-automatic reconstruction of neuron morphology. Front. Neuroinform. 6, 4 (2012). · DOI
  7. Peng, H., Ruan, Z., Long, F., Simpson, J. H. & Myers, E. W. V3D enables real-time 3D visualization and quantitative analysis of large-scale biological image data sets. Nat. Biotechnol. 28, 348–353 (2010). · DOI · Paper page
  8. Peng, H., Long, F., Zhao, T. & Myers, E. Proof-editing is the bottleneck of 3D neuron reconstruction: The problem and solutions. Neuroinformatics 9, 103–105 (2011). · DOI
  9. Rodriguez, A., Ehlenberger, D. B., Hof, P. R. & Wearne, S. L. Rayburst sampling, an algorithm for automated three-dimensional shape analysis from laser scanning microscopy images. Nat. Protoc. 1, 2152–2161 (2006). · DOI
  10. Pool, M., Thiemann, J., Bar-Or, A. & Fournier, A. E. NeuriteTracer: A novel ImageJ plugin for automated quantification of neurite outgrowth. J. Neurosci. Methods 168, 134–139 (2008). · DOI
  11. Ho, S.-Y. et al. NeurphologyJ: An automatic neuronal morphology quantification method and its application in pharmacological discovery. BMC Bioinform. 12, 230 (2011). · DOI
  12. Boulan, B. et al. AutoNeuriteJ: An ImageJ plugin for measurement and classification of neuritic extensions. PLoS ONE 15, e0234529 (2020). · DOI
  13. Arganda-Carreras, I. et al. Trainable weka segmentation: A machine learning tool for microscopy pixel classification. Bioinformatics 33, 2424–2426 (2017). · DOI · Paper page
  14. Bolte, S. & Cordelières, F. P. A guided tour into subcellular colocalization analysis in light microscopy. J. Microsc. 224, 213–232 (2006). · DOI
  15. Cruz-Matías, I. et al. Sphericity and roundness computation for particles using the extreme vertices model. J. Comput. Sci. 30, 28–40 (2019). · DOI
  16. Arganda-Carreras, I., Fernández-González, R., Muñoz-Barrutia, A. & Ortiz-De-Solorzano, C. 3D reconstruction of histological sections: Application to mammary gland tissue. Microsc. Res. Tech. 73, 1019–1029 (2010). · DOI
  17. Dougherty, R. & Kunzelmann, K. Computing local thickness of 3D structures with ImageJ. Microsc. Microanal. 13, 1678–1679 (2007). · DOI
  18. Dunn, K. W., Kamocka, M. M. & McDonald, J. H. A practical guide to evaluating colocalization in biological microscopy. Am. J. Physiol. Cell Physiol. 300, C723–C742 (2011). · DOI
  19. Encarnacion-Rivera, L., Foltz, S., Hartzell, H. C. & Choo, H. Myosoft: An automated muscle histology analysis tool using machine learning algorithm utilizing FIJI/ImageJ software. PLoS ONE 15, e0229041 (2020). · DOI
  20. Justić, H. et al. Redefining the Koizumi model of mouse cerebral ischemia: A comparative longitudinal study of cerebral and retinal ischemia in the Koizumi and Longa middle cerebral artery occlusion models. J. Cereb. Blood Flow Metab Off. J. Int. Soc. Cereb. Blood Flow Metab. 42, 2080–2094 (2022). · DOI
  21. Jagečić, D., Petrović, D. J., Šimunić, I., Isaković, J. & Mitrečić, D. The oxygen and glucose deprivation of immature cells of the nervous system exerts distinct effects on mitochondria, mitophagy, and autophagy, depending on the cells’ differentiation stage differentiation stage. Brain Sci. 13, 910 (2023). · DOI
  22. Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012). · DOI · Paper page
  23. van Ginneken, B. Fifty years of computer analysis in chest imaging: Rule-based, machine learning, deep learning. Radiol. Phys. Technol. 10, 23–32 (2017). · DOI
  24. Deo, R. C. Machine learning in medicine. Circulation 132, 1920–1930 (2015). · DOI · Paper page
  25. Gurcan, M. N. et al. Histopathological image analysis: A review. IEEE Rev. Biomed. Eng. 2, 147–171 (2009). · DOI · Paper page
  26. Patel, L., Shukla, T., Huang, X., Ussery, D. W. & Wang, S. Machine learning methods in drug discovery. Molecules 25, 5277 (2020). · DOI · Paper page
  27. Koizumi, J. et al. Experimental studies of ischemic brain edema. Jpn. J. Stroke 8, 1–8 (1986). · DOI · Paper page
  28. Longa, E. Z., Weinstein, P. R., Carlson, S. & Cummins, R. Reversible middle cerebral artery occlusion without craniectomy in rats. Stroke 20, 84–91 (1989). · DOI · Paper page
  29. Stirling, D. R. et al. Cell Profiler 4: Improvements in speed, utility and usability. BMC Bioinform. 22, 433 (2021). · DOI
  30. Pani, G. et al. MorphoNeuroNet: An automated method for dense neurite network analysis. Cytom. Part A J. Int. Soc. Anal. Cytol. 85, 188–199 (2014). · DOI

Cited by

Sorted by citations

See all 18

Metadata

Authors
Šimunić, I., Jagečić, D., Isaković, J., Dobrivojević Radmilović, M., & Mitrečić, D.
Venue
Scientific Reports
Publisher
Springer Science and Business Media LLC
Published
2024-03-28
Volume
14
Issue
1
ISSN
2045-2322
DOI
10.1038/s41598-024-57650-6
PMC
PMC10978859
Cited by
18
References
2
Language
en
Update notices
correction (2024-04-19)
Licence
Funders
Hrvatska Zaklada za Znanost (UIP-2017-05-8082, IP-2022-10-4656)

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review

How to cite

Cite CoChat for the checks on this page, and the paper for its own findings.

This page
CoChat. (2026). Citation check: Lusca: FIJI (ImageJ) based tool for automated morphological analysis of cellular and subcellular structures. https://app.cochat.ai/papers/10.1038/s41598-024-57650-6
The paper
Šimunić, I., Jagečić, D., Isaković, J., Dobrivojević Radmilović, M., & Mitrečić, D. (2024). Lusca: FIJI (ImageJ) based tool for automated morphological analysis of cellular and subcellular structures. Scientific Reports. https://doi.org/10.1038/s41598-024-57650-6