Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 05 February 2024
Sec. Interdisciplinary Physics
This article is part of the Research Topic Biophotonics and Nanophotonics: Theory and Applications View all articles

Mueller matrix polarization interferometry of optically anisotropic architectonics of biological tissue object fields: the fundamental and applied aspects

Alexander Ushenko,
Alexander Ushenko1,2*Alexander DubolazovAlexander Dubolazov2Jun Zheng
Jun Zheng1*Oxana BakunOxana Bakun3Mykhaylo GorskyMykhaylo Gorsky2Yuriy UshenkoYuriy Ushenko4Olexandra LitvinenkoOlexandra Litvinenko3Ivan GordeyIvan Gordey4Chen ZheboChen Zhebo1Valeriy SklyarchukValeriy Sklyarchuk2
  • 1Photoelectric Information Center, Taizhou Institute of Zhejiang University, Taizhou, China
  • 2Optics and Publishing Department, Chernivtsi National University, Chernivtsi, Ukraine
  • 3Department of Forensic Medicine and Medical Jurisprudence, Bucovinian State Medical University, Chernivtsi, Ukraine
  • 4Computer Science Department, Chernivtsi National University, Chernivtsi, Ukraine

The presented results concern the diagnostic application of vector-parameter polarization mapping and polarization-interference phase scanning applied to layered azimuths and ellipticity polarization maps obtained with the help of digital holographic reproduction for phase-inhomogeneous tissue layers in the female reproductive system. The differential diagnosis of the pathological changes in optically anisotropic polycrystalline tissue components, specifically “benign (fibroids) and precancerous (endometriosis),” was investigated. All studies within the representative groups of the histological samples were conducted using circularly polarized laser radiation, ensuring azimuthal invariance and polarization measurement reliability. Integral and layered maps and polarization azimuths and ellipticity random variable distribution histograms were obtained from the microscopic images of histological sections from the female reproductive system tumors. Systematized tables present statistical moments which characterize azimuths and ellipticity polarization maps for fibroids and endometriosis tissues. It is shown that the statistical parameters (diagnostic markers) are the most sensitive to benign and precancerous changes in the female reproductive system tissues. The results of information analysis, including the accuracy determination of the diagnostic vector-parameter polarization and polarization-interference methods for detecting and differentiating the samples of fibroids and precancerous endometriosis tissues in the female reproductive system, are presented.

1 Introduction

Enhancing the capacity to detect cancer in the female reproductive system, particularly in the early preclinical stages such as endometriosis, is a crucial aspect of addressing the global challenge of “sustainable development.” Endometriosis involves the growth of tissues resembling the uterine inner layer outside the uterus, leading to pain and potential infertility [1]. Worldwide, endometriosis affects approximately 10% (190 million) of women and girls of reproductive age [2]. It is a chronic disease, with symptoms such as severe, interfering pain during menstruation, intercourse, defecation, and/or urination; chronic pelvic pain; abdominal bloating; nausea; increased fatigue; and sometimes, also depression, anxiety disorder, and infertility. Because of the variety of symptoms of endometriosis, it is not easy for health professionals to diagnose the disease, and many people who suffer from it know little about it. Therefore, it sometimes takes quite a long time from the onset of symptoms to diagnosis [3]. Currently, there is no universally recognized treatment for endometriosis, and in most cases, treatment is limited to the elimination of symptoms [4].

In recent years, modern optical–physical methods of oncological diagnostics have progressed. A significant achievement in this field was the development of polarization-sensitive optical coherence tomography (PSOCT) techniques which have become promising in this direction. The basic principles and achievements of the PSOCT of biological objects over the past 25 years have been consistently set out in a series of review articles [57]. This laser polarization-interference technique not only makes it possible to obtain layered images from small depths (up to 2 mm) of biological tissues [8, 9] but also provides the possibility of obtaining coordinate distributions of the Jones and Mueller matrix elements and anisotropy maps. This was most effectively demonstrated in studies of transparent eye tissues [1012].

PSOCT functionality has been extended with the studies of layered polarization and matrix images of radiation which was scattered (depolarized) on the biological tissues of various morphological structures and physiological states [8]. In particular, it was possible to accurately diagnose fibrosis and differentiate tumor sites with low fibrosis [13], along with benign and malignant tumors of the larynx [14]. However, the sensitivity and resolution of such systems was limited due to the distortion effect produced by the high level of depolarized laser speckle background. This leads to low contrast for polarizing images which are obtained by such tissue-layer illumination.

Based on this, the creation and development of high-resolution (0.5 μm–1 μm) polarization-interference microscopy systems, which functionally complement PSOCT techniques, is relevant. This technique, combined with digital holographic phase scanning of a multiple-scattered object field, can not only significantly increase the contrast of polarization images but also provide a Mueller matrix reconstruction of the optical anisotropy parameters for the samples, such as histological biological tissue sections.

To address the aforementioned challenges, further fundamental and applied development of optical diagnostic methods is needed. The desired results can be achieved by integrating modern laser polarimetry techniques [1521] based on Mueller matrix formalism [2229] and digital holography [3033], which provide exhaustive information about the layered optical anisotropic hierarchical structure of biological tissues. The promising initial results have been demonstrated in the high-precision (∼95%) comprehensive early detection of prostate cancer using layered Mueller matrix mapping [32].

However, the obtained results remain somewhat empirical. Currently, there is no information about the possibilities with this method for the samples with different optical thicknesses, morphological structures, and various pathologies. From a physical point of view, the task is to form relationships between the structure of polycrystalline networks of biological tissues and the layered distribution of polarization parameters (azimuth and ellipticity).

Our research will concentrate on advancing a novel and logically complementary PSOCT method that combines Mueller matrix interferometry with the microscopic imaging of the histological sections which are made from diverse biological tissues. These sections exhibit various morphological structures and physiological states. Methodologically, this approach relies on a comprehensive Mueller matrix model to characterize the morphological structure of biological tissues. In this model, the biological tissues are conceptualized as amorphous and polycrystalline structures [20, 21]. The polarization distributions of azimuths and ellipticity in the microscopic image of biological tissues are shaped by the polycrystalline component. From a crystal-optical perspective, this component is regarded as a combination of the structural mechanisms (such as linear birefringence and dichroism in fibrillar networks composed of diverse protein molecules, such as collagen, myosin, and elastin) and chiral mechanisms (involving circular birefringence and dichroism in optically active protein complexes with molecular chiral rings) [3436].

The amorphous component of biological tissue consists of fats, lipids, biochemical acids, and other optically inactive molecules. Its impact on the object field parameters is characterized by the absorption and attenuation of optical radiation. The proposed Mueller matrix interferometry technique is a technological synthesis which incorporates the Mach–Zehnder polarization-interference platform, digital holographic reconstruction, and step-by-step phase scanning to analyze complex amplitude distributions in the microscopic images of biological tissue histological sections.

Consequently, our approach introduces a novel capability for experimentally implementing the phase selection of laser radiation components with varying scattering multiplicities. This allows for the isolation of an informatively relevant “low-fold” or single-scattered polarization component in microscopic images. The distributions of polarization parameters (azimuths and ellipticity) in this component are maximally and unambiguously correlated with the structural parameters of the polycrystalline component in the native biological tissue section. The anticipated outcomes include a substantial enhancement in method sensitivity and the provision of contemporary diagnostic insights in the pathological changes within depolarizing biological layers’ polycrystalline components.

Our research is concentrated on advancing and experimentally validating Mueller matrix polarization-interferometry methods for characterizing multiple scattered fields in diffusive biological tissue layers within the female reproductive system. The objective is to investigate the potential of the layered digital holographic selection of components with varying scattering multiplicities in the polarization-inhomogeneous field. This leads to the enhanced accuracy of the differential diagnosis between benign conditions, such as fibroids and precancerous tissues, specifically endometriosis, within the female reproductive system.

2 Materials and methods

In this section, we will briefly review (without reducing the completeness of the analysis) the main theoretical relationships within the framework of the linear birefringence (LB) approximation. It describes the processes of polarization structure formation in the cases of single- and multiple-radiation scattering on diffusive tissue layers. This allows us to calculate the amount of object laser field components formed by the histological sections of the female reproductive system with benign and precancerous changes.

2.1 Object field Stokes polarimetry

To obtain reliable (azimuthally invariant) information within the representative histological sections of the tissue biopsy samples from the female reproductive system, illumination is performed using right-circularly (⨂) polarized beams with an azimuth angle α0. The Stokes vector of this polarized beam is given by the following expression [1519] (Eq. 1):

VS0=1001.(1)

2.1.1 “Single-scattering” interaction

The polarization properties of an optically anisotropic (n, the index of linear birefringence, LB) protein fibril (with a geometric size d and spatial orientation of the optical axis ρ) in the volume of the tissue sample from the female reproductive system are represented by a Mueller matrix birefringence operator [16, 19] (Eq. 2):

Wj=10000ω22ω23ω240ω32ω33ω340ω42ω43ω44j=
=10000cos22ρ+sin22ρcosδcos2ρsin2ρ1cosδsin2ρsinδ0cos2ρsin2ρ1cosδsin22ρ+cos22ρcosδcos2ρsinδ0sin2ρsinδcos2ρsinδcosδj.(2)

Here, δ=2πλΔnd is the phase shift between the linearly orthogonally polarized components of the laser beam amplitude, where λ is the wavelength.

The process of single-scattering transformation through the polarization structure of the probing beam S0 by a local (ith) fibril (ρj;δj) is described by the following matrix equation:

VSi*αi,βi=WiVS0.(3)

As a result, a wave, which is formed as a result of a local act of scattering, can be described with the following azimuth αi and ellipticityβi of polarization (Eqs 4, 5):

αj=0.5arctanVS3*VS2*=0.5arctancotan2ρjqρj,(4)
βj=0.5arcinVS4*VS1*=0.5arcincosδjgδj.(5)

Here, VSi=1;2;3;4* represents the parameters of the Stokes vector of the singly scattered laser beam VSi*αj,βj. Thus, the singly scattered component in the object field of the protein fibrillar network in the tissue of the female reproductive system constitutes two polarization distributions:

• “Azimuthal-orientational” distribution denoted by Λαj;ρj

• “Phase” distribution denoted by Ψβj;δj

2.1.2 “n-fold” interaction

For multiple interactions with the protein fibrillar network in the female reproductive system, matrix Eq. 3 for a circularly polarized laser probe VS0 takes the following form [20, 21] (Eqs 68):

VSjnαjn;βjn=WnWn1W2W1VS0.(6)

This optical scenario generates an ensemble of the following random values of azimuth αjn and ellipticity βjn:

αjn=0.5arctanj=1nVS3j*j=1nVS2j*j=1nqρj,(7)
βjn=0.5arcsinj=1nVS4j*j=1nVS1j*j=1ngδj.(8)

As a result, a polarization-variant component of the diffuse field will be formed with a different distribution of the azimuthal Λnj=1nqρj and ellipticity values Ψnj=1ngδj of polarization.

Thus, the polarization structure Ρnαn;βn for the laser field, which was multiply scattered in the tissue of the female reproductive system, can be represented as the superposition of “orientational” and “phase” components (Eq. 9).

Ρnαn;βn=Λnj=1nqρj+Ψnj=1ngδj.(9)

The coherence of laser radiation facilitates an alternative analytical description. It is focused on the amplitude aspects of the processes shaping the probing field within the polarization structure of the diffuse layers in female reproductive tissue.

2.1.3 Amplitude analysis

For coherent laser fields, there exists a direct relationship between the parameters of the Stokes vector and the orthogonal components (Uxj and Uyj) of complex amplitudes [33]. Based on this, the previously obtained expressions (4 and 5) for the polarization parameters can be rewritten as follows (Eqs 10, 11):

αjρj=0.5arctanUxjUyj*+Uxj*UyjUxjUxj*UyjUyj*,(10)
βjδj=0.5arciniUxjUyj*Uxj*UyjUxjUxj*+UyjUyj*.(11)

Indeed, in parallel with this scenario, another process takes place: the interference summation of differently polarized partial coherent waves, leading to the formation of another polarization-variant component in the object field which is scattered by the female reproductive diffuse tissue section.

2.1.4 Interference interaction

For the orthogonal components Uxj and Uyj of the complex amplitudes U1 and U2, respectively, of two partially singly scattered coherent waves at a local point in the object field, the following interference equations can be written (Eqs 12, 13):

Uxj=Ux1+Ux2j=Ux1+Ux2+2Ux1Ux2cosφx12j,(12)
Uyj=Uy1+Uy2j=Uy1+Uy2+2Uy1Uy2cosφy12j,(13)

where Uxj;Uyj are the modules of complex amplitudes and φx12 and φy12 are the phase shifts between (Ux1;Ux2 and Uy1;Uy2, respectively.

The following equations can be written for the resulting orthogonal components of amplitudes Ux and Uy due to the “n-fold” interaction in the scattering process of the laser probe with optical heterogeneities (Eq. 14):

Ux=j=1nUxj;Uy=j=1nUyj;φxy=i=1nφxji=1nφyj.(14)

The interference addition of two phase-shifted orthogonal components Ux and Uy by a phase difference φxy results in an elliptically polarized wave [24] (Eq. 15):

X2Ux2+Y2Uy22XYUxUycosφxy=sin2φxy,(15)

with the following “interference” values (I) of azimuth αI and ellipticity βI:

αI=0.5arcsinsin2UyUx/1+tan2δxycos22UyUx,(16)
βI=0.5arctantanδxysin2UyUx/1+tan2δxycos22UyUx.(17)

As we can see, because of the secondary interference of differently polarized coherent partial waves, a polarization-inhomogeneous component (denoted as ΛIαjn and ΨIβjn) is also formed (according to Eqs. 16, 17). An object field is produced by the scattering on biological tissues from the female reproductive sphere.

2.1.5 Resulting field

Thus, the polarization structure of the laser field, scattered by the female reproductive sphere tissues, can be represented as the superposition of the following components:

Φα;β=Λαj;ρj+Ψβj;δj+Λnj=1nqρj+Ψnj=1ngδj+ΛIαjn;φj++ΨIβjn;φj.(18)

3 Experimental setup and measurement methodology

The Stokеs-polarimetric mapping scheme based on the Mach–Zehnder interferometer is a generalization of the polarization-interferometry scheme [3033], which is shown in Figure 1.

FIGURE 1
www.frontiersin.org

FIGURE 1. Optical scheme for polarization-interference mapping of the Stokes vector parameters.1: He–Ne laser; 2: collimator—“O”; 3 and 11: beam splitters—“BS”; 4 and 5: mirrors—“M”; 6, 9, and 13: polarizers—“P”; 7 and 10: quarter wave plates—“QP”; 8: object; 12: polarization objective—“O”; 14: digital camera—“CCD”; and 15: personal computer—“PC.”

The parallel (Ø=2000μm) illuminating (“probing”) beam of He–Ne (λ=0.6328μm) laser (1), formed by the spatial-frequency filter (2), with a 50% beam splitter (3) is divided into “object” and “reference” ones. The “object” beam with the help of a rotating mirror (5) is directed through the polarizing filter (6–7), which forms a right-circularly polarized beam [relation (1)]. Next, such an “object” beam is sent in the direction of the female biological tissue histological section layer (8). The polarization-inhomogeneous image of the biological tissue histological sections (8) is projected by the strain-free objective (12) into the digital camera (14) plane.

The “reference” beam is directed by the mirror (4) through the polarization filter (9–10), which forms its right-circular state of polarization. The resulting beam is sent into the polarization-inhomogeneous image plane histological sections (8), which forms the microscopic image. Next, a circularly polarized “reference” wave using a beam splitter (11) and a polarizing lens (12) is superimposed on a microscopic image of a native histological section of female biological tissue. As a result, an interference pattern is formed, the coordinate intensity distribution of which is recorded using a digital camera (14) through a polarizer (13).

Before carrying out the measurements of biological tissues, the experimental device passed the metrological certification with the introduction of model objects (“clean air,” “linear polarizer,” and “phase plates 0.25λ” and “0.5λ”). For each type of model object, 50 measurements of ellipticity were carried out. Within each series of measurements, the average (mathematical expectation) and the standard deviation of the measured values were determined according to the standard methodology of probability theory. As a result, the final value of the standard deviation was β=0.0003rad.

The methodology, for a layer-by-layer measurement of the object field Stokes vector-parameter (polarization maps αm,n and βm,n) distributions, is using complex amplitudes Ex and Ey from a digital holographic reconstruction, followed by the phase scanning θk of an object field. It is presented in [3033]. However, detailed information is not provided in this work. For a better understanding of the further discussion, we provide a brief overview of the 3D digital holographic scanning method.

4 The method of 3D polarimetry phase scanning of the object field

1. Using polarizing filters (6–7) and (9–10), circular polarization () is sequentially formed in the “irradiating” (Ir) and “reference” (Re) parallel laser beams, IrRe.

2. For the circular polarization () state, two partial interference patterns are recorded through the polarizer–analyzer (14) with the orientation of the transmission plane at angles Ω=00;Ω=900.

3. The analytical processing of female biological tissue histological section microscopic interference images was carried out using the digital Fourier transform FTυ,ν [3033] (Eq. 19):

FTx;yυ,ν=1M×Nm=0M1n=0N1Ix,yΩ=00;900m,nexpi2πm×υM+n×νN,(19)

where {IxΩ=00m,n=UxΩ=00Ux*Ω=00m,n;IyΩ=900m,n=UyΩ=900Uy*Ω=900m,n;Ux,y are the orthogonal components of complex amplitude for different orientations Ω=00;Ω=900; * denotes the complex conjugation operation; υ,ν are the spatial frequencies; and m=1120;n=960 are the quantity of pixels of the CCD camera.

4. The results of the digital Fourier transform (relation (18)) are used to obtain complex amplitude distributions according to the following algorithms (Eq. 20, 21):

U00UxΩ=00,(20)
U900UyΩ=900expiφxφy.(21)

5. By the means of the stepwise (Δθ)-phase (θk) scanning of the reconstructed complex amplitude (relations (20) and (21)) field using algorithms (10)–(11), we obtain polarization parameters αθk,m,n and βθk,m,n coordinate distributions.

6. The resulting set of polarization maps pαθk,m,n;βθ,m,n was analyzed in a statistical approach using the following algorithms to calculate the mean (Z1), variance (Z2), skewness (Z3), and kurtosis (Z4) [20, 21] (Eq. 22):

Z1=1Kj=1Kpj;Z2=1Kj=1Kp2j;Z3=1Z231Kj=1Kp3j;Z4=1Z241Kj=1Kp4j,(22)

where K is the CCD pixel quantity.

7. By reducing the magnitude of the phase parameter θk, the “boundary” value θk* is determined to isolate the biological object field single-scattered component, starting from which the aggregate Zi=1;2;3;4αθk*,m,nconst and Zi=1;2;3;4βθk*,m,nconst.

5 Object of investigations

The objects of the study were native histological sections of the endometrium, which were prepared according to the standard procedure using a microtome with rapid freezing.

Optically thin (group 1) and optically thick (group 2) samples of histological sections were created, including specimens with benign and precancerous changes in the female reproductive tissue. These samples were used to conduct a comparative analysis of the diagnostic efficiency of the Mueller matrix method under the conditions of scattering with varying multiplicities.

The optical-geometric parameters of the samples are presented in Table 1.

TABLE 1
www.frontiersin.org

TABLE 1. Female reproductive system histological section optical and geometric parameters.

The measurement of the uterine tissue sample extinction coefficient was carried out according to the standard procedure of light attenuation measurement [37] using an integral light-scattering sphere [38, 39]. The sample preparation procedure was conducted in accordance with the principles of the Declaration of Helsinki and in compliance with the International Conference on Harmonization, Good Clinical Practice and local regulatory requirements. The study was reviewed and approved by the appropriate independent ethics committees.

Within groups 1 and 2, two representative subgroups of the histological section samples were formed:

1. Benign samples (myoma): control group 1.1 (14 samples) and control group 2.1 (14 samples)

2. Premalignant samples (endometriosis): experimental group 1.2 (14 samples) and experimental group 2.2 (14 samples)

To ensure the statistical significance of the Mueller matrix polarization-interferometry method, we set the confidence interval value as p0.05.

To determine the representativeness of each sample (the number of samples in each of the groups 1.1, 1.2, 2.1, and 2.2), we calculated the standard deviation σ2 of each of the statistical moments, Z_(i = 1; 2; 3; 4) (n), until its value satisfied the condition σ20.025. This condition was achieved when the number of samples in each group was n=14.

From a morphological point of view, benign and premalignant tumors are characterized by different polycrystalline structures. Benign tumors (myoma) have a collagen and myosin fibrous network. Premalignant tumors (endometriosis) have a “developed” fibrous network of the connective tissue component, which is oriented in specific directions and scales [40]. Similar polarization manifestations of the polycrystalline newly formed fibrillar growth networks of other types of biological tissues have been studied using laser polarimetry in a number of other publications [20, 21, 3436].

The combined selection, for studying two types of uterine tissue samples, allows for the comparative investigation of their effectiveness under different scattering conditions by two methods: Stokes polarimetry and phased selection of polarization components (Λαj;ρj+Ψβj;δj+Λnj=1nqρj+Ψnj=1ngδj+ΛIαjn;φj+ΨIβjn;φj) of the diffuse field and the extraction of the minimally scattered (most diagnostically important) component Λαj;ρj+Ψβj;δj.

Obtaining optically thin histological sections with controlled parameters of single scattering using a microtome is a complex and not always achievable task. It is further justified that polarization-interference scanning of the biological layer object fields (Equation 18) can use optically thick samples.

The results, obtained from the methods of Stokes-polarimetric mapping and polarization-interference layered-phase scanning of object field research sample tissues from the female reproductive system, were analyzed. The informational analysis involves the use of operational characteristics from evidence-based medicine [41].

6 Experimental results and discussion

This section of the article provides a systematic representation of the experimental results along with their discussion, focusing on the structure of azimuths αm,n and ellipticity βm,n maps within the integral, layered, and polarization-inhomogeneous object fields of diffusive samples from the female reproductive system (see Table 1).

6.1 Integral- and layered-polarization azimuth maps

Figure 2 shows the following: integral-polarization azimuth αm,n map fragments (1) and (2); layered phase-polarization azimuth αθk,m,n maps for θk=π/4 fragments (5) and (6) and for θk=π/8 fragments (9) and (10); histogram Gα fragments (3) and (4); histogram Gθk=π/4,α fragments (7) and (8); and Gθk=π/8,α fragments (11) and (12).

FIGURE 2
www.frontiersin.org

FIGURE 2. Coordinate and statistical structure of azimuth-polarization maps for optically thin (τ0.01) myoma and precancerous endometrium histological sections, left and right columns, respectively (explanation in the text).

The analysis of the obtained data revealed the following:

1. Coordinate inhomogeneity and individual statistical structure of the integral (Figure 2; fragments (1) and (2)) and layered phase (Figure 2; fragments (5) and (6) and (9) and (10)) azimuth polarization maps of the myoma and endometrium histological sections

2. Asymmetric histograms of the integral (Gα; Figure 2; fragments (3) and (4)) and layered phase (Gθk=π/4,α,Gθk=π/8,α; Figure 2; fragments (7) and (8) and (11) and (12))

The statistical analysis results of the polarization azimuth maps αm,n and αθk,m,n are presented in Table 2.

TABLE 2
www.frontiersin.org

TABLE 2. Balanced accuracy threshold levels.

It shows the following:

• All the first- and fourth-order statistical moments that characterize the histograms α and Gθk=π/4,α,Gθk=π/8,α are different from zero, Zi=1,2,3,4α,θk=π/4,α,θk=π/8,α0

• The higher-order statistical moments that characterize the asymmetry and excess of distributions Gα і Gθk=π/4,α,Gθk=π/8,α are the most sensitive to changes in the polarization azimuth map

• Value change scenarios of the first- and fourth-order statistical moments, which characterize the polarization azimuth maps of myoma (M) and endometrium (E) histological sections, Zi=1,2М,α,θk<Zi=1,2Е,α,θkZi=3,4М,α,θk>Zi=3,4Е,α,θk

From a physical point of view, the obtained results can be attributed to the pathological alterations in the connective tissue component. These alterations lead to an increase in the “neo-formed” structural anisotropy (the relationship between (2) and (5)) of collagen networks through the geometric and orientational growth of protein fibers. As a consequence, there is an enlargement in the distribution of the orientations of the optical axes ρm,nmax of fibrillar crystals, along with their birefringence δm,n=2πλΔndm,n.

The aforementioned precancerous pathological changes will lead to an increase in the mean Z1 (α,ρ,δ) and variance Z2 (α,ρ,δ) of the polarization azimuth fluctuations in the corresponding carcinoma biopsy histological section microscopic image. At the same time, the magnitude of the third- and fourth-order statistical moments that characterize the asymmetry and excess of the distributions αm,n will decrease (relation (22)), i.e., Z3,4 (α,ρ,δ) (Table 3).

TABLE 3
www.frontiersin.org

TABLE 3. Statistical moments of the first and fourth orders characterizing the layer-by-layer-polarization azimuth maps.

The results of the information analysis (in relation to (23)–(25)) revealed a set of diagnostic markers and the following accuracy levels in the benign and precancerous condition differential diagnoses:

• Integrated-polarization azimuth maps: satisfactory level AcZ4α=82.1%

• Layered-polarization azimuth maps: good level AcZ3,4α,θk=π/4,π/8=85.7%89.2%

6.2 Integral- and layered-polarization ellipticity maps

Figure 3 presents the investigation results of the benign myoma and precancerous endometriosis object fields (notations are similar to those in Figure 2).

FIGURE 3
www.frontiersin.org

FIGURE 3. Coordinate and statistical structure of layered ellipticity-polarization maps for optically thin (τ0.01) myoma and precancerous endometrium histological sections, left and right columns, respectively (explanation is given in the text).

The analysis of the obtained data revealed a united heterogeneous topographic and coordinate structure of the experimentally obtained integral- (Figure 3; fragments (1) and (2)) and layer-by-layer-polarization ellipticity maps (Figure 3; fragments (5) and (6) and (9) and (10)).

Table 4 illustrates the statistical analysis results of the βm,n and βθk,m,n maps.

TABLE 4
www.frontiersin.org

TABLE 4. First- and fourth-order statistical moment characterization of the layer-by-layer-polarization azimuth maps for optically thick histological sections.

Here, we found the following:

• Distinct from the normal and individual statistical distributions of the polarization parameters, Zi=1,2,3,4β,θk=π/4,β,θk=π/8,β0.

• Diagnostic markers of the benign and precancerous changes: higher-order statistical points that characterize the asymmetry and excess of the distributions Gβ and Gθk=π/4,β,Gθk=π/8,β.

• Precancerous changes in the polycrystalline component of the endometrial samples (group 2) are accompanied by an increase in the mean Z1 (β,ρ,δ) and variance Z2 (β,ρ,δ). At the same time, the magnitude of the third- and fourth-order statistical moments decreases, Z3,4 (β,ρ,δ) and Zi=1,2М,β,θk<Zi=1,2Е,β,θk;Zi=3,4М,β,θk>Zi=3,4Е,β,θk.

Information analysis results revealed a significant increase in the accuracy of the differential diagnosis of the benign and precancerous conditions in the female reproductive sphere:

• The integral-polarization ellipticity maps show a good level of accuracy: AcZ3,4β=85.7%

• The layer-wise phase-polarization ellipticity maps demonstrate the following accuracy levels: good level AcZ3β,θk=π/4=89.2%, very good level AcZ4β,θk=π/8=91.7%, and very good level AcZ3,4β,θk=π/8=92.8%

• Starting from θkπ/8, the balanced accuracy remains practically unchanged, AcZ3,4β,θk=π/8const

From a physical point of view, the discovered regularities can be attributed to the fact that for θkπ/8, the “single-scattering” regime occurs within the volume of optically thin biological layers. As a result, the maximum possible polarimetric differential diagnosis accuracy of the female reproductive system benign and precancerous states is achieved.

6.3 Optically thick-layer polarization azimuth maps

Figure 4 displays the maps and histograms illustrating the distribution of polarization azimuth values αm,n in the microscopic images of the histological sections from diffuse benign uterine fibroids and precancerous endometriosis tissues. The notations are consistent with those employed in Figures 2, 3.

FIGURE 4
www.frontiersin.org

FIGURE 4. Coordinate and statistical structure of the layer-by-layer-polarization azimuth maps for optically thick (τ>0.01) myoma and precancerous endometrium histological sections, left and right columns, respectively (explanation is given in the text).

The analysis of the obtained results revealed, similar to the studies on optically thin layers, topographic and coordinate inhomogeneity in the integral (Figure 4; fragments (1) and (2))- and layer-wise phase (Figure 4; fragments (5) and (6) and (9) and (10))-polarization azimuth maps in benign and precancerous uterine tissues of optically thick-histological section microscopic images.

The research results on the diagnostic sensitivity of integral- and phase-polarization azimuth mapping in optically thick-histological section microscopic images of both groups showed predominant secondary interference mechanisms’ influence. These mechanisms form the elliptically polarized diffuse component of the object field in almost all phase planes (except for θkπ/8).

As a result, unsatisfactory levels of diagnostic sensitivity for integral-polarization azimuth maps were obtained (Table 2), AcZ3,4α<80% .

Layer-by-layer phase-polarization azimuth maps illustrate the improvement in accuracy to a satisfactory level, AcZ3,4α,θk=π/8=85.7%.

6.4 Optically thick-layer polarization ellipticity maps

Figure 5 presents the polarization maps βm,n of uterine tissue histological sections.

FIGURE 5
www.frontiersin.org

FIGURE 5. Coordinate and statistical structure of layer-by-layer-polarization ellipticity maps for optically thick (τ>0.01) myoma and precancerous endometrium histological sections, left and right columns, respectively (explanation is given in the text).

The analysis of the obtained results revealed a sufficiently similar, albeit complex, topographic and coordinate structure of polarization ellipticity maps (Figure 5; fragments (1), (2), (5), and (6) and (9) and (10)).

Table 5 illustrates the statistical analysis results of the polarization ellipticity maps βm,n and βθi,m,n.

TABLE 5
www.frontiersin.org

TABLE 5. Statistical moments of the first- and fourth-order characterization of the layer-by-layer ellipticity-polarization maps for optically thick histological sections.

Within the framework of the statistical approach for the analysis of experimental data, we found the following:

• Different from the normal statistical distributions of polarization ellipticity maps, Zi=1,2,3,4β,θk=π/4,β,θk=π/8,β0

• Diagnostic markers are higher-order statistical moments that characterize the asymmetry and excess of distributions Gβ and Gθk=π/4,β,Gθk=π/8,β

We already noted that the pathological changes in the polycrystalline component of diffuse endometrial samples from group 2 are accompanied by opposite scenarios of the transformation of the values of statistical markers, showing an increase in the average Z1 (β,ρ,δ) and variance Z2 (β,ρ,δ) and, conversely, a decrease of the third- and fourth-order statistical moments Z3;4 (β,ρ,δ).

A comparative analysis of the layer-by-layer maps βθi,m,n revealed that as the value of θk decreases, there is a consistent dispersion Z2θk;β decrease (up to three times) (Table 5). At the same time, the value of Z3,4 (β,ρ,δ increases up to three to four times.

A further reduction of the parameter θ<π8 of the object field phase section practically does not lead to changes in the statistical structure of the polarization maps. This fact indicates that for such phase planes, the most “diagnostically favorable” mode of single scattering in the volume of uterine tissue diffuse layers is realized. Therefore, it is in the phase plane θ=π8 that we observe the greatest difference (40%–50%) between the statistical moments Z3;4θk;β, which characterize the distributions of Ψβj;δj in object fields of the samples of benign and precancerous uterine tissues.

As a result, for female reproductive sphere-diffuse tissue samples in precancerous conditions, the accuracy of the balanced differential diagnosis in the phase plane θ=π8 reaches the same level as that for optically thin samples (Table 4), AcZ3,4β,θk=π/8=92.8%.

7 Conclusion

1. The experimental results of the diagnostic applications, employing vector-parametric polarization mapping and polarization-interference phase scanning methods on object fields of female reproductive sphere tissues, were presented and physically analyzed. This approach involves a layer-by-layer digital polarization map holographic reproduction, contributing to the differential diagnosis of benign (myoma) and precancerous (endometriosis) conditions.

2. For all types of pathology, it is an established complex, and for the individual statistics of integral- and layered-polarization maps, all statistical moments of the first and fourth orders, which characterize the histograms Gα,β and Gθk=π/4,α,β,Gθk=π/8,α,β, are different from zero, Zi=1,2,3,4α,β,θk=π/4,α,β,θk=π/8,α,β0.

3. The following maximum precancerous condition differential diagnosis accuracy levels using optically thin histological sections were revealed:

3.1. Integral-polarization azimuth maps: satisfactory AcZ4α=82.1%

3.2. Layer-by-layer phase-polarization azimuth maps: good AcZ3,4α,θk=π/4,π/8=85.7%89.2%

3.3. Integral-polarization ellipticity maps: good AcZ3,4β=85.7%

3.4. Layer-by-layer phase-polarization ellipticity maps: good AcZ3β,θk=π/4=89.2%, very good AcZ4β,θk=π/8=91.7%, and equivalent to very good AcZ3,4β,θk=π/8=92.8%

4. For optically thick histological sections, the obtained results are given as follows:

4.1. The diagnostic sensitivity for integral-polarization azimuth maps was unsatisfactory, AcZ3,4α<80%

4.2. Layer-wise phase-polarization azimuth maps demonstrated an improvement in the accuracy level to a satisfactory level, AcZ3,4α,θk=π/8=85.7%

5. The precancerous states of the female reproductive system diffuse tissue sample detection accuracy for polarization ellipticity maps in the phase plane θ=π8 reached a level similar to those of optically thin samples, AcZ4β,θk=π/8=92.8%.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by the Ethics Committee of the Bureau of Forensic Medicine of the Chernivtsi Region. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin because histologically verified preparations from the archive of the Bureau of Forensic Medicine of the Chernivtsi Region were used as the experimental samples.

Author contributions

AU: project administration, writing–original draft, and writing–review and editing. AD: investigation and writing–original draft. JZ: conceptualization, methodology, supervision, and writing–review and editing. OB: conceptualization, formal analysis, investigation, and writing–review and editing. MG: software, visualization, and writing–review and editing. YU: data curation, investigation, and writing–original draft. OL: formal analysis, investigation, and writing–review and editing. IG: investigation and writing–review and editing. CZ: data curation and writing–review and editing. VS: writing–review and editing, formal Analysis, visualization.

Funding

The authors declare that financial support was received for the research, authorship, and/or publication of this article. The authors acknowledge support from the National Research Foundation of Ukraine, Project 2022.01/0034, and the Scholarship of the Verkhovna Rada of Ukraine for Young Scientists and Doctors of Science.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. World Health Organisation. Classification of diseases (ICD) (2022). Available at: https://www.who.int/standards/classifications/classification-of-diseases.

Google Scholar

2. Zondervan KT, Becker CM, Missmer SA. Endometriosis. New Engl J Med (2020) 382(13):1244–56. doi:10.1056/nejmra1810764

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Agarwal SK, Chapron C, Giudice LC, Laufer MR, Leyland N, Missmer SA, et al. Clinical diagnosis of endometriosis: a call to action. Am J Obstet Gynecol (2019) 220(4):354.e1–354.e12. doi:10.1016/j.ajog.2018.12.039

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Johnson NP, Hummelshoj L, Abrao MS, Adamson GD, Allaire C, Amelung V, et al. Consensus on current management of endometriosis. Hum Reprod (2013) 28(6):1552–68. doi:10.1093/humrep/det050

PubMed Abstract | CrossRef Full Text | Google Scholar

5. de Boer JF, Milner TE. Review of polarization sensitive optical coherence tomography and Stokes vector determination. J Biomed Opt (2002) 7(3):359. doi:10.1117/1.1483879

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Spandana KU, Krishna KM, Mazumder N. Polarization-resolved Stokes-Mueller imaging: a review of technology and applications. Lasers Med Sci (2019) 34(7):1283–93. doi:10.1007/s10103-019-02752-1

PubMed Abstract | CrossRef Full Text | Google Scholar

7. He Y, Li K, Li W, Qiu Y, Li D, Wang C, et al. Polarization coherency matrix tomography. J Biophotonics (2023) 16(9): e202300093. doi:10.1002/jbio.202300093

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Yang D, Yuan Z, Hu M, Liang Y. Zebrafish brain and skull imaging based on polarization-sensitive optical coherence tomography. J Biophotonics (2022) 15(12): e202200112. doi:10.1002/jbio.202200112

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Park JE, Zhou X, Kwon DY, Kim SW, Lee H, Jung MJ, et al. Application of polarization sensitive-optical coherence tomography to the assessment of phase retardation in subpleural cancer in rabbits. Tissue Eng Regenerative Med (2021) 18(1):61–9. doi:10.1007/s13770-020-00318-9

CrossRef Full Text | Google Scholar

10. Willemse J, Maximilian GOG, Verbraak FD, de Boer JF. In vivo 3D determination of peripapillary scleral and retinal layer architecture using polarization-sensitive optical coherence tomography. Translational Vis Sci Tech (2020) 9(11):21–1. doi:10.1167/tvst.9.11.21

CrossRef Full Text | Google Scholar

11. Baumann B, Augustin M, Lichtenegger A, Harper DJ, Muck M, Eugui P, et al. Polarization-sensitive optical coherence tomography imaging of the anterior mouse eye. J Biomed Opt (2018) 23(08):1. doi:10.1117/1.jbo.23.8.086005

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Hariri LP, Adams DC, Applegate MB, Miller AJ, Roop BW, Villiger M, et al. Distinguishing tumor from associated fibrosis to increase diagnostic biopsy yield with polarization-sensitive optical coherence tomography. Clin Cancer Res (2019) 25(17):5242–9. doi:10.1158/1078-0432.ccr-19-0566

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Tao K, Sun K, Ding Z, Ma Y, Kuang H, Zhao H, et al. Catheter-based polarization sensitive optical coherence tomography using similar mueller matrix method. IEEE Trans Biomed Eng (2020) 67(1):60–8. doi:10.1109/tbme.2019.2908031

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Burns JA, Kim KH, deBoer JF, Rox Anderson R, Zeitels SM. Polarization-sensitive optical coherence tomography imaging of benign and malignant laryngeal lesions. Otolaryngol Head Neck Surg (2011) 145(1):91–9. doi:10.1177/0194599811403078

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Tuchin VV. Tissue optics: light scattering methods and instruments for medical diagnosis. Bellingham, Washington, USA: Spie Press (2015).

Google Scholar

16. Ghosh N. Tissue polarimetry: concepts, challenges, applications, and outlook. J Biomed Opt (2011) 16(11):110801. doi:10.1117/1.3652896

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Jacques SL. Polarized light imaging of biological tissues. In: Handbook of biomedical Optics2. Boca-Raton: CRC Press (2011). p. 649–69.

Google Scholar

18. Layden D, Ghosh N, Vitkin IA. Quantitative polarimetry for tissue characterization and diagnosis. In: Advanced biophotonics: tissue optical sectioning. Boca-Raton: CRC Press (2013). p. 73–108.

CrossRef Full Text | Google Scholar

19. Vitkin A, Ghosh N, de Martino A. Tissue polarimetry. In: DL Andrews, editor. Photonics: scientific foundations, technology and applications. Goboken. John Wiley and Sons, Ltd (2015). p. 239–321.

CrossRef Full Text | Google Scholar

20. Angelsky OV, Ushenko YA, Ushenko A. Statistical, correlation, and topological approaches in diagnostics of the structure and physiological state of birefringent biological tissues. In: VV Tuchin, editor. Handbook of photonics for biomedical science. London: CRC Press, Taylor and Francis Publishing (2010). p. 283–322.

CrossRef Full Text | Google Scholar

21. Ushenko YA, Boychuk TM, Bachynsky VT. Diagnostics of structure and physiological state of birefringent biological tissues: statistical, correlation and topological approaches. In: VV Tuchin, editor. Handbook of coherent-domain optical methods: biomedical diagnostics, environmental monitoring, and materials science. New York: Springer New York (2013). p. 107–48.

CrossRef Full Text | Google Scholar

22. Gil JJ, Ossikovski R. Polarized light and the mueller matrix approach. Boca-Raton: CRC Press (2022).

Google Scholar

23. Arteaga O, Kahr B. Mueller matrix polarimetry of bianisotropic materials [Invited]. J Opt Soc America B-optical Phys (2019) 36(8):F72–2. doi:10.1364/josab.36.000f72

CrossRef Full Text | Google Scholar

24. Lee HR, Lotz C, Kai Groeber Becker F, Dembski S, Novikova T. Digital histology of tissue with Mueller microscopy and FastDBSCAN. Appl Opt (2022) 61(32):9616. doi:10.1364/ao.473095

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Novikova T, Kim M, Lee HR, Ossikovski R, Malfait-Jobart A, Lamarque D. Optical diagnosis of gastric tissue biopsies with Mueller microscopy and statistical analysis. J Eur Opt Soc Rapid Publications (2022) 18(2):10. doi:10.1051/jeos/2022011

CrossRef Full Text | Google Scholar

26. Lee HR, Li P, Sang T, Lotz C, Groeber-Becker F, Dembski S, et al. Digital histology with Mueller microscopy: how to mitigate an impact of tissue cut thickness fluctuations. J Biomed Opt (2019) 24(07):076004. doi:10.1117/1.JBO.24.7.076004

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Li P, Lee HR, Chandel S, Lotz C, Groeber-Becker FK, Dembski S, et al. Analysis of tissue microstructure with Mueller microscopy: logarithmic decomposition and Monte Carlo modeling. J Biomed Opt (2020) 25(01):1. doi:10.1117/1.jbo.25.1.015002

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lee HR, Sang T, Li P, Lotz C, Groeber-Becker F, Dembski S, et al. Mueller microscopy of anisotropic scattering media: theory and experiments. France: Centre pour la Communication Scientifique Directe (2018).

Google Scholar

29. Ma H, He H, Ramella-Roman JC. Mueller matrix microscopy. In: JC Ramella-Roman, and T Novikova, editors. Polarized light in biomedical imaging and sensing. Cham: Springer (2023).

CrossRef Full Text | Google Scholar

30. Peyvasteh M, Tryfonyuk L, Ushenko V, Syvokorovskaya A, Dubolazov AV, Vanchulyak O, et al. 3D Mueller-matrix-based azimuthal invariant tomography of polycrystalline structure within benign and malignant soft-tissue tumours. Laser Phys Lett (2020) 17(11):115606–6. doi:10.1088/1612-202x/abbee0

CrossRef Full Text | Google Scholar

31. Ushenko V, Hogan B, Dubolazov A, Piavchenko G, Kuznetsov S, Ushenko A, et al. 3D Mueller matrix mapping of layered distributions of depolarisation degree for analysis of prostate adenoma and carcinoma diffuse tissues. Scientific Rep (2021)(1) 11. doi:10.1038/s41598-021-83986-4

CrossRef Full Text | Google Scholar

32. Yu U, Sdobnov A, Dubolazov A, Grytsiuk M, Bykov A, Meglinski I, et al. Stokes-correlometry analysis of biological tissues with polycrystalline structure. IEEE J Selected Top Quan Electro (2019) 25(1):1–12. doi:10.1109/jstqe.2018.2865443

CrossRef Full Text | Google Scholar

33. Ushenko VA, Hogan BT, Dubolazov A, Grechina AV, Boronikhina TV, Gorsky M, et al. Embossed topographic depolarisation maps of biological tissues with different morphological structures. Scientific Rep (2021)(1) 11. doi:10.1038/s41598-021-83017-2

CrossRef Full Text | Google Scholar

34. Ushenko VA, Sdobnov A, Mishalov WD, Dubolazov AV, Олар OV, Bachinskyi VT, et al. Biomedical applications of Jones-matrix tomography to polycrystalline films of biological fluids. J Innovative Opt Health Sci (2019) 12(06):12. doi:10.1142/s1793545819500172

CrossRef Full Text | Google Scholar

35. Ushenko AG, Dubolazov OV, Ushenko VA, Oyu N, Olar OV. Fourier polarimetry of human skin in the tasks of differentiation of benign and malignant formations. Appl Opt (2016) 55(12):B56. doi:10.1364/ao.55.000b56

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Peyvasteh M, Dubolazov A, Popov A, Ushenko A, Meglinski I, Meglinski I. Two-point Stokes vector diagnostic approach for characterization of optically anisotropic biological tissues. J Phys D (2020) 53(39):395401–1. doi:10.1088/1361-6463/ab9571

CrossRef Full Text | Google Scholar

37. Marchesini R, Bertoni A, Andreola S, Melloni E, Sichirollo AE. Extinction and absorption coefficients and scattering phase functions of human tissues in vitro. Appl Opt (1989) 28(12):2318–24. doi:10.1364/ao.28.002318

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Edwards DK, Gier JT, Nelson KE, Roddick RD. Integrating sphere for imperfectly diffuse samples. J Opt Soc America (1961) 51(11):1279–9. doi:10.1364/josa.51.001279

CrossRef Full Text | Google Scholar

39. Goodman JW. Laser speckle and related phenomena. In: JC Dainty, editor. Topics in applied Physics. Heidelberg: Springer-Verlag Berlin (1975). p. 9–75.

CrossRef Full Text | Google Scholar

40. Rogers PAW, D’Hooghe TM, Fazleabas A, Gargett CE, Giudice LC, Montgomery GW, et al. Priorities for endometriosis research: recommendations from an international consensus workshop. Reprod Sci (2009) 16(4):335–46. doi:10.1177/1933719108330568

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Robinson SP. Principles of forensic medicine. Cambridge: Cambridge University Press (1996).

Google Scholar

Keywords: polarization, interference, holography, microscopic image, optical anisotropy, statistical moments, balanced accuracy, female reproductive system

Citation: Ushenko A, Dubolazov A, Zheng J, Bakun O, Gorsky M, Ushenko Y, Litvinenko O, Gordey I, Zhebo C and Sklyarchuk V (2024) Mueller matrix polarization interferometry of optically anisotropic architectonics of biological tissue object fields: the fundamental and applied aspects. Front. Phys. 11:1302254. doi: 10.3389/fphy.2023.1302254

Received: 26 September 2023; Accepted: 29 December 2023;
Published: 05 February 2024.

Edited by:

Wei Zhao, Northwest University, China

Reviewed by:

V. N. Du Le, The University of Alabama in Huntsville, United States
Radu Hristu, Polytechnic University of Bucharest, Romania

Copyright © 2024 Ushenko, Dubolazov, Zheng, Bakun, Gorsky, Ushenko, Litvinenko, Gordey, Zhebo and Sklyarchuk. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Alexander Ushenko, o.ushenko@chnu.edu.ua; Jun Zheng, dbzj@netease.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.