CC BY
6
Practical Aspects of Using to Assess the Morphology of Biolo;
DOI: 10.17691/stm2023.15.3.04 Received October 24, 2022
A.M. Ignatova, DSc, Senior Researcher1; Researcher2; M.A. Zemlyanova, MD, DSc, Associate Professor, Chief Researcher2; Professor, Department of Environmental Protection3;
O.B. Naimark, DSc, Professor, Head of the Laboratory1;
N.V. Zaitseva, MD, DSc, Academician of the Russian Academy of Sciences, Scientific Advis
11nstitute of Continuous Media Mechanics of the Ural Branch of the Russian Academy of Science
a Branch of the Perm Federal Research Center, Ural Branch of the Russian Academy of Sciences,
13A Lenina St., Perm, 614990, Russia;
2Federal Scientific Center for Medical and Preventive Health Risk Management Technologies,
82 Monastyrskaya St., Perm, 614045, Russia;
3Perm National Research Polytechnic University, 29 Komsomolsky Prospekt, Perm, 614990, Russia
The aim of the study is to identify practical aspects of using multifractal formalism to assess the morphology of biological tissues.
Materials and Methods. The objects of the study were histological images of lung tissues of Wistar rats without pathology and with detected pathological changes, obtained at 50*, 100*, 200* magnifications. Image processing was carried out using the ImageJ/Fiji universal software. The multifractal spectrum of the images, processed to obtain a linear contour, was calculated with the use of FracLac — a module for ImageJ. This module was used to determine the scaling exponent (the function of the Renyi exponent, x(q)) and the singularity spectrum itself.
Results. The singularity spectra for tissues with no pathology have signs of multifractality. The image spectrum of tissue with pathology is shifted to the left relative to the spectrum characteristic of tissue without pathology. A decrease in the spectral height in the presence of pathology indicates a "simplification" of the alveolar pattern, which is presumably associated with the presence of widespread vasculitis, since it causes areas of hemorrhage to appear on the image; this leads to leveling the contour of the alveolar pattern, reducing the surface area of the alveoli and emerging areas inflamed by erythrocytes. At lower magnification, images with pathology lose signs of multifractality.
Conclusion. Correct results of assessing multifractal spectra of histological images can be achieved at 200* magnification and preprocessing to obtain linear contours. Significant differences between the morphological structure of lung tissues with and without pathology are observed when comparing the height, width, and position of the spectrum relative to the origin.
Key words: image analysis; histology; lung tissue; fractal analysis; alveolar pattern.
How to cite: Ignatova A.M., Zemlyanova M.A., Naimark O.B., Zaitseva N.V. Practical aspects of using multifractal formalism to assess the morphology of biological tissues. Sovremennye tehnologii v medicine 2023; 15(3): 35, https://doi.org/10.17691/stm2023.15.3.04
This is an open access article under the CC BY 4.0 license (https://creativecommons.org/licenses/by/4.0/).
Introduction
In terms of classical concepts, histological analysis does not imply obtaining calculable indicators, however, in recent years, in addition to it, image analysis is increasingly used to obtain quantitative characteristics [1-3]. Quantitative characteristics provide the comparison of the results of a large number of experiments and their application for modeling [4-5]. When quantifying the morphology of biological tissues, the problem of determining evaluation criteria arises. Currently, various parameters for the analysis of histological images, in particular, a degree of tissue differentiation based on the Gleason grading system, are used in scientific and clinical practice [6]. This approach and others close to it are
applicable only to the assessment of specific pathologies, such as neoplasms, and not suitable for diagnosing and classifying other changes. In a number of works [7-9], in investigation of non-tumor pathologies, the criteria for studying images of biological tissues by means of image analysis are discussed. It is noted that fractal dimension is applicable to assess the organization of neurons, and indicators of dendritic geometry are used to describe the alveolar pattern.
The structure of histological images often cannot be described by a single fractal set. It is a complex of several fractals superimposed on each other, each of which has its own dimension [10]. Such a multifractal nature of histological images is due to the fact that biological processes develop according to power laws
Corresponding author: Anna M. Ignatova, e-mail: iampstu@yandex.ru
Multifractal Formalism for Assessing the Morphology of Biological Tissues CTM j 2023 j vol.
15 j No.3 35
and, therefore, their phase portraits are characterized by a non-uniform distribution of a certain measure, which can be described using not one, but several scales [7]. That is why a promising method for obtaining image characteristics is multifractal formalism [11-13].
The characteristics of the multifractal spectrum of biological tissue images are proposed to be used as a universal evaluation approach [8].
Direct calculation of the multifractal spectrum of histological images, as a rule, is ineffective due to the presence of artifacts that can lead to an unreliable result, for this reason the use of the multifractal formalism requires specialized image creation. The search for methodological approaches to the use of multifractal image analysis in relation to the assessment of the morphology of biological tissues and objects is relevant.
The aim of the study is to identify practical aspects of using multifractal formalism to assess the morphology of biological tissues.
Materials and Methods
The objects of the study were histological images of lung tissues of Wistar rats without pathology and with detected pathological changes (Figure 1), obtained at 50*, 100*, 200* magnifications. Histological material was obtained from experiments, the conditions of which met the requirements of the European Convention for the Protection of Vertebrate Animals used for Experimental or Other Scientific Purposes (ETS No.123), and the Ethics Committee of the Federal Scientific Center for Medical and Preventive Health Risk Management Technologies (Perm, Russia).
The study examined images of lung tissues with diagnosed hyperplasia of lymphoid tissue associated
with the bronchial wall with an admixture of eosinophils and widespread vasculitis, interstitial pneumonia, and alveolar pulmonary edema.
The images were obtained using microscopy of lung tissue specimens. For their making, the organs were fixed in 10% aqueous neutral formalin solution and then dehydrated in ascending alcohols, then soaked in chloroform and paraffin, and poured with Histomix homogenized paraffin medium (BioVitrum, Russia). 4-^m thick histosections were obtained using a JUNG SM2000R sledge microtome (Leica Microsystems, Germany) and stained with hematoxylin and eosin. Inspection and imaging of the morphological structure of tissues was performed using an Axiostar light-optical microscope (Carl Zeiss, Germany).
It is noted in [14] that one of the possible means of minimizing the effect of artifacts in the image analysis is binarization, edge detection and its skeletonization. Image processing was performed using ImageJ/Fiji universal software (open-source software, developed by Wayne Rasband, National Institutes of Health, USA).
Before making a multifractal analysis, the images were converted to grayscale black and white, and then to two-color ones by the binarization method with a threshold value of 220 (Figure 2 (a)).
The value of lacunarity was estimated for the binary images. This indicator is a measure of the heterogeneity of filling the space with a pattern, which in this case is an alveolar pattern. An increase in the lacunarity value indicates an increase in the number of chaotic dips in the pattern.
In this study, the measure of lacunarity (L) is a change in image density when scanning it with a grid with cells of various sizes (square cells ranging in size from 1 to
Figure 1. Lung tissue of Wistar rats:
(a) without pathology; (b) with hyperplasia of the lymphoid tissue associated with the bronchial wall and widespread vasculitis; (c) with alveolar pulmonary edema; (d) with interstitial pneumonia; ><200
////////////////////^^
36 CTM I 2023 I vol. 15 j No.3 A.M. Ignatova, M.A. Zemlyanova, O.B. Naimark, N.V. Zaitseva
Figure 2. Lung tissue images of Wistar rats prepared for analysis:
(a) the result of binarization; (b) the result of edge detection; (c) the result of contour skeletonization
0.1
60 pixels were used). The following formula is used to calculate lacunarity:
L=(a/M)2,
where a is the standard deviation of the number of pixels (for a binary image) of a fractal aggregate in grid cells of a given size; m is the average value of the mass of the aggregate or the average number of pixels in cells of a given size.
Lacunarity was determined by the slope of the regression line in the coordinates of e/L (e is the grid cell size in fractions of the raster size) and by the range of L values for £=0.01-0.6. When graph plotting, the lacunarity was calculated by the formula
L=(a/M)2+1, where a is used to plot the graph uniformly at zero standard deviation (a=0).
The edge detection of the resulting binary images was performed using the outline function (Figure 2 (b)), and then its skeletonization was carried out (Figure 2 (c)). As a result of processing, the images of lung tissues presented images of a double linear contour of the alveolar pattern, which were used to obtain a multifractal spectrum. To calculate the multifractal spectrum of images processed to obtain a linear contour, we used the module for ImageJ — FracLac. Using this module, the scaling exponent (the function of the Renyi exponent, T(q)) and the singularity spectrum were determined. The following parameters were determined for the singularity spectrum: the width and height of the spectrum; the width and height of the right and left branches [15].
Results and Discussion
The results of assessing the lacunarity of the alveolar pattern are shown in Figure 3 and in Table 1. Figure 4 shows singularity spectra characteristic of processed
2,2 2,0 1.8 1.6 1.4 1.2 1.0
o
o • . 4
8... m ■ • a • a
---- ........ ____ l < ........ vT-a ■ ■ 8 —c
• A A • • 4» ----i; % Ä ■iTT^ ---V
• i « • • % 1 A fi t r- f 4 ■ A vi
0.2
0.3
0.4
0.5
0.6
Without pathology, *50 Without pathology, *200 With pathology, x100 Alveolar edema, *200 Linear (without pathology, *50) Linear (without pathology, *200) Linear (with pathology, *100) Linear (alveolar edema, *200)
Without pathology, *100 With pathology, *50 With pathology, *200 Interstitial pneumonia, *200 Linear (without pathology, *100) ■ Linear (with pathology, *50) i Linear (with pathology, *200) . Linear (interstitial pneumonia, *200)
Figure 3. Image lacunarity functions
Table 1
Lacunarity characteristics of the alveolar pattern
Lacunarity characteristics
Image characteristics Coefficient of the slope of linear regression Minimum value Maximum value
Without pathology
*50 -0.26 1.0 1.83
*100 -0.29 1.01 1.40
x200 -0.66 1.02 1.66
With pathology
*50 0.04 1.0 1.61
*100 0.01 1.0 1.85
x200 -0.71 1.0 1.69
Alveolar edema, *200 -1.08 1.03 2.11
Interstitial pneumonia, *200 -0.69 1.0 1.78
Multifractal Formalism for Assessing the Morphology of Biological Tissues CTM j 2023 j vol.
15 j No.3 37
in the range of £=0.4^-0.5 is apparently associated with the physiological characteristics of the tissues, namely with the appearance of blood vessels in the alveolar pattern. To the greatest extent, the value of lacunarity identifies alveolar edema. For images of tissues with pathology at 50* and 100* magnification, the lacunarity index is close to zero, which, along with the linear nature of the dependence of the scaling exponents for these images (Figure 5), is a sign of monofractality. At 200* magnification, the singularity spectrum for tissues with pathology, on the contrary, has pronounced signs of multifractality, but the characteristic of the spectrum differs from similar indicators characteristic of tissues without morphological pathologies. In the presence of pathology, the height of the spectrum was lower by 0.2; the spectrum width was more by 0.3; at the same time, the ratio between the widths of the right and left branches remained similar to that observed for tissues without pathologies.
For tissues with no pathology detected, the singularity spectra have signs of multifractality, regardless of the magnification at which the images were obtained. However, with an increase in the degree of magnification of histological images the parameters of the spectra change, namely, the width of the spectrum uniformly increases and its height decreases; at the same time, the ratios of the parameters of the right and left branches are preserved. At 200* magnification, the height of the right and left branches slightly decreases relative to similar figures observed at 100* magnification. These changes in the spectra with increasing magnification correspond to a decrease in the number of alveoli in sight and a more distinct contour of the alveolar pattern.
On the whole, the image spectrum of tissue with pathology is shifted to the left relative to the spectrum characteristic of tissue without pathology. The latter, together with a histological analysis, indicates the growing dominance of lymphoid tissue, which is characteristic of its pathological hyperplasia.
A decrease in the width of the singularity spectrum, which is characteristic of images with pathology, indicates a more uniform distribution of the structural components of such a tissue than tissues without pathology, i.e. the periodicity (rhythm) of the pattern of the alveolar pattern is disturbed; in relation to the lung tissue, this indicates a reduction in the alveolar lumens. A decrease in the height of the spectrum in the presence of pathology indicates a "simplification" of the alveolar pattern, which is presumably associated with the pneumonia, *200 0.84 1.57 0.51 0.70 0.33 0.56 presence of widespread vasculitis,
2.0 1.9 1.8 1.7 1.6 1.5 1.4 1.3 1.2 1,1
1.0 0.9 0.8 07 0.6 0.5
1
D0=1.65
Do=1.84; DH.81
D0=1.78; Dq=1.73
Dq=1.5S; Dq=1.46 0o=1.5S; Dq=1.4S
Do=1.71; Dq=1.64
20 1.30 1.40 1.50 1.60 1.70 1.80 1.90 2.00 2.10 2.20 2.30 2.40 2.50
a
Without pathology, *50 ----Without pathology, *100
Without pathology, *200 - With pathology, *50
With pathology, *100 -With pathology, *200
Alveolar edema, *200 ...... Interstitial pneumonia, *200
Figure 4. Image singularity spectra
D0 — monofractal dimension, Dq — informational dimension
images of histological lung tissues; the spectra characteristics are given in Table 2.
For tissues with pathology, the lacunarity degree of the alveolar pattern was lower than that for tissues without pathology. A local increase in the lacunarity
Table 2
Multifractal spectra characteristics of histological images
Spectrum characteristics
Image characteristics Spectrum width, amax-amin Spectrum height, 1W Width of the right branch, amax"°0 Height of the right branch, f(a0)-f(amax) Width of the left branch, a0-amin Height of the left branch, f(a0)-f(amin)
Without pathology
*50 0.71 1.83 0.44 0.95 0.27 0.71
*100 0.93 1.77 0.61 1.20 1.32 0.89
*200 1.01 1.70 0.63 1.04 0.38 0.77
With pathology
*50 Monofractal pattern
*100 Monofractal pattern
*200 0.68 1.49 0.44 0.64 0.24 0.29
Alveolar edema, *200 0.83 1.57 0.44 0.65 0.38 0.52
Interstitial
//////////////////////^^^^
38 CTM I 2023 I vol. 15 j No.3 A.M. Ignatova, M.A. Zemlyanova, O.B. Naimark, N.V. Zaitseva
20 15 10 5 О -5 -10 -15 -20 -25 -30
since it causes areas of hemorrhage to come in sight on the image. This leads to the alignment of the contour of the alveolar pattern, the surface area of the alveoli is reduced, and there appear areas filled with red blood cells.
Most notably, at lower magnification, images with pathology lose signs of multifractality. Presumably, the absence of signs of multifractality is also associated with smoothing the contours of the pattern due to the general swelling and focal distribution of erythrocytes. In this regard, the characteristic details are revealed only at higher magnification.
The singularity spectra of images of specified pathological manifestations, namely alveolar edema and interstitial pneumonia, at 200* magnification are close to the spectrum characteristic of pathological changes in general. They have the same height and almost equal width, which, however, is greater than for pathological changes in general, which is associated with a violation of the alternation of alveolar spaces in the pattern characteristic of lung tissues. The spectrum observed in alveolar edema is most shifted to the left, which confirms the assumption of the dominance of a pathological component in the tissue structure.
Conclusion
Methodological approaches are proposed for using multifractal formalism to assess the morphology of biological tissues. It has been established that the correct results of evaluating the multifractal spectra of histological images can be achieved at 200* magnification and preprocessing to obtain linear contours. Evaluation criteria for the quantitative characteristic of biological tissues by the method of multifractal analysis are the statistical characteristics of the singularity spectrum. Significant differences between the morphological structure of lung tissues with and without pathology are recorded when comparing the height, width and position of the spectrum relative to the origin. Lacunarity is not a universal evaluation criterion.
Study funding. The work was carried out without additional funding.
Conflicts of interest. There are no conflicts of interest.
References
1. Huang H.K. Biomedical image processing. Crit Rev Bioeng 1981; 5(3): 185-271.
2. Chung Y., Shin S., Shim H., Sohn J.Y., Lee D.E., Lee H.,
-6
-4
-Without pathology, *50 . Without pathology, *200 -With pathology, *100 . Alveolar edema, *200
Without pathology, *100 With pathology, *50 With pathology, *200 Interstitial pneumonia, *200
Figure 5. Scaling exponents of images
Eom H.S., Kim K.G., Kong S.Y. Development of an automated image analyzer for microvessel density measurement in bone marrow biopsies. Ann Lab Med 2020; 40(4): 312-316, https:// doi.org/10.3343/alm.2020.40.4.312.
3. Zhang Y., Xie M., Xue R., Tang Q., Zhu X., Wang J., Yang H., Ma C. A novel cell morphology analyzer application in head and neck cancer. Int J Gen Med 2021; 14: 9307-9314, https://doi.org/10.2147/ijgm.s341420.
4. Albers J., Pacile S., Markus M.A., Wiart M., Vande Velde G., Tromba G., Dullin C. X-ray-based 3D virtual histology-adding the next dimension to histological analysis. Mol Imaging Biol 2018; 20(5): 732-741, https://doi.org/10.1007/s11307-018-1246-3.
5. Sokolova N.A., Orel V.E., Gusynin A.V., Selezneva A.A., Kolesnik S.V. Algorithm for image computerization analysis of histological specimens. Problemi informacijnih tehnologij 2012; 11: 129-139.
6. Gorban N.A., Kudaibergenova A.G. Current views of the Gleason grading system. Onkourologia 2010; 6(1): 69-75.
7. Shupletsov Yu.V. Direct calculation of the multifractal spectrum for images of biomedical drugs. Izvestia Rossijskogo gosudarstvennogo pedagogiceskogo universiteta im. A.I. Gercena 2015; 173: 88-96.
8. Ampilova N.B., Soloviev I.P., Shupletzov Y.V. Multifractal spectrum as a classification sign for biomedical preparations images. In: Proceedings of 7th International Conference on Communications, Electromagnetics and Medical Applications. Athens: School of Electrical and Computer Engineering; 2013; p. 98-101.
9. Majeed H., Sridharan S., Mir M., Ma L., Min E., Jung W., Popescu G. Quantitative phase imaging for medical diagnosis. J Biophotonics 2017; 10(2): 177-205, https://doi.org/10.1002/ jbio.201600113.
10. Naimark O.B., Nikityuk A.S., Nebogatikov V.O., Grishko V.V. Assessment of nonlinear dynamics of cellular structure damage as a promising method of personalized oncodiagnostics. Al'manah kliniceskoj mediciny 2018; 46(8): 742-747, https://doi.org/10.18786/2072-0505-201846-8-742-747.
q
Multifractal Formalism for Assessing the Morphology of Biological Tissues
CTM | 2023 | vol. 15 | No.3 39
11. Gelashvili D.B., ludin D.I., Rozenberg G.S., Yakimov V.N., Solntsev L.A. Fraktaly i mul'tifraktaly v bioekologii [Fractals and multifractals in bioecology]. Nizhny Novgorod: Izdatel'stvo Nizhegorodskogo gosudarstvennogo universiteta im. N.I. Lobachevskogo; 2013; 370 p.
12. Bozhokin S.V., Parshin D.A. Fraktaly i mul'tifraktaly [Fractals and multifractals]. Moscow-Izhevsk: NITs; 2001; 128 p.
13. Mandelbrot B. Fraktal'naya geometriya prirody [Fractal
geometry of nature]. Moscow: Institut komp'yuternykh issledovaniy; 2002; 656 p.
14. Zaitseva N.V., Zemlyanova M.A., Ignatova A.M., Stepankov M.S. Morphological changes in lung tissues of mice caused by exposure to nano-sized particles of nickel oxide. Nanotechnol Russ 2018; 13(7-8): 393-399.
15. Ampilova N., Soloviev I. Application of fractal and multifractal analysis algorithms to image segmentation and classification. WSEAS Trans Biol Biomed 2016; 13(3): 14-21.
40 СТМ j 2023 j vol. 15 j No.3
A.M. Ignatova, M.A. Zemlyanova, O.B. Naimark, N.V. Zaitseva
