Do different turbid media with matched bulk optical properties also exhibit similar polarization properties?

Mueller matrix polarization imaging and quantitative parameters analysis method

Mueller matrix polarization imaging is a new biomedical optical imaging method that can generate both polarization and isotropic intensity images of structures of the biological tissue sample surface. In this paper, a Mueller polarization imaging system in the reflection mode is described for obtaining the Mueller matrix of the specimens. Diattenuation, phase retardation, and depolarization of the specimens are derived by using the conventional Mueller matrix polarization decomposition method and a newly proposed direct method. The results show that the direct method is more convenient and faster than the conventional decomposition method. The polarization parameter combination method is then presented in which any two of the diattenuation, phase retardation, and depolarization parameters are combined, and three new quantitative parameters are defined in order to reveal more detailed anisotropic structures. The images of in vitro samples are presented to demonstrate the capability of the parameters introduced.

Polarization memory rate as a metric to differentiate benign and malignant tissues

Non-invasive optical methods for cancer diagnostics, such as microscopy, spectroscopy, and polarimetry, are rapidly advancing. In this respect, finding new and powerful optical metrics is an indispensable task. Here we introduce polarization memory rate (PMR) as a sensitive metric for optical cancer diagnostics. PMR characterizes the preservation of circularly polarized light relative to linearly polarized light as light propagates in a medium. We hypothesize that because of well-known indicators associated with the morphological changes of cancer cells, like an enlarged nucleus size and higher chromatin density, PMR should be greater for cancerous than for the non-cancerous tissues. A thorough literature review reveals how this difference arises from the anomalous depolarization behaviour of many biological tissues. In physical terms, though most biological tissue primarily exhibits Mie scattering, it typically exhibits Rayleigh depolarization. However, in cancerous tissue the Mie depolarization regime becomes more prominent than Rayleigh. Experimental evidence of this metric is found in a preliminary clinical study using a novel Stokes polarimetry probe. We conducted in vivo measurements of 20 benign, 28 malignant and 59 normal skin sites with a 660 nm laser diode. The median PMR values for cancer vs non-cancer are significantly higher for cancer which supports our hypothesis. The reported fundamental differences in depolarization may persist for other types of cancer and create a conceptual basis for further developments in polarimetry applications for cancer detection.

Discriminating turbid media by scatterer size and scattering coefficient using backscattered linearly and circularly polarized light

The effects of scatterer size and scattering coefficient on backscattered linearly and circularly polarized light are investigated through Stokes polarimetry. High-SNR polarization modulation/synchronous detection measurements are corroborated by polarization-sensitive Monte Carlo simulations. Circular degree of polarization (DOP) is found to be sensitive to scatterer size, but is equivocal at times due to helicity flipping effects; linear DOP appears to be mostly dependent on the medium scattering coefficient. We exploit these trends to generate a DOP_C - DOP_L response surface which clusters turbid samples based on these medium properties. This work may prove useful in biomedicine, for example in noninvasive assessment of epithelial precancer progression.

Controlling the optical pathlength in continuous-wave reflectance spectroscopy using polarization

We investigate potential improvements of continuous-wave diffuse reflectance spectroscopy within highly scattering media by employing polarization gating. Simulations are used to show the extent at which the effective optical pathlength varies in a typical scattering medium as a function of the optical wavelength, the total level of absorption, and the selected polarization channels, including elliptical and circular polarization channels. Experiments then demonstrate that a wavelength dependent polarization gating scheme may reduce the prior knowledge required to solve the problem of chromophore quantification. This is achieved by finding combinations of polarization channels which have similar effective optical pathlengths through the medium at each wavelength.

Evaluating β-amyloidosis progression in Alzheimer's disease with Mueller polarimetry

We applied the wide-field Mueller imaging polarimetry for the screening of formalin-fixed paraffin-embedded samples of mouse brain tissue at different stages of brain β-amyloidosis in Alzheimer's disease (AD). The accumulation of amyloid-beta (Aβ) deposits throughout the brain tissue is one of the key pathological hallmarks observed with the AD progression. We demonstrate that the presence of Aβ plaques influences the properties of backscattered polarized light, in particular, its degree of depolarization. By means of statistical analysis, we demonstrate that the high-order statistical moments of depolarization distributions, acquired with the multi-spectral Mueller imaging polarimetry, can be used as sensitive markers of the growing presence of Aβ plaques. The introduced label-free polarimetric approach has a potential to facilitate the current practice of the histopathology screening in terms of diagnosis accuracy, time and cost efficiency.

Role of scattering and birefringence in phase retardation revealed by locus of Stokes vector on Poincaré sphere

SIGNIFICANCE

Biological tissues are typically characterized by high anisotropic scattering and may also exhibit linear form birefringence. Both scattering and birefringence bias the phase shift between transverse electric field components of polarized light. These phase alterations are associated with particular structural malformations in the tissue. In fact, the majority of polarization-based techniques are unable to distinguish the nature of the phase shift induced by birefringence or scattering of light.

AIM

We explore the distinct contributions of scattering and birefringence in the phase retardation of circularly polarized light propagated in turbid tissue-like scattering medium.

APPROACH

The circularly polarized light in frame of Stokes polarimetry approach is used for the screening of biotissue phantoms and chicken skin samples. The change of optical properties in chicken skin is accomplished by optical clearing, which reduces scattering, and mechanical stretch, which induces birefringence. The change of optical properties of skin tissue is confirmed by spectrophotometric measurements and second-harmonic generation imaging.

RESULTS

The contributions of scattering and birefringence in the phase retardation of circularly polarized light propagated in biological tissues are distinguished by the locus of the Stokes vector mapped on the Poincaré sphere. The phase retardation of circularly polarized light due to scattering alterations is assessed. The value of birefringence in chicken skin is estimated as 0.3  ×  10^-3, which agrees with alternative studies. The change of birefringence of skin tissue due to mechanical stretch in the order of 10^-6 is detected.

CONCLUSIONS

While the polarimetric parameters on their own do not allow distinguishing the contributions of scattering and birefringence, the resultant Stokes vector trajectory on the Poincaré sphere reveals the role of scattering and birefringence in the total phase retardation. The described approach, applied independently or in combination with Mueller polarimetry, can be beneficial for the advanced characterization of various types of malformations within biological tissues.

Designing phantoms to accurately replicate circular depolarization in biological scattering media

We introduce an iterative method for designing optical phantoms that are able to replicate the depolarization profiles of various target media, including colloidal suspensions of Intralipid, bovine milk, and ex vivo samples of ovine kidney cortex tissue. The designed phantoms comprise spherical scattering particles with fine-tuned size distributions and are capable of simultaneously reproducing spatially resolved intensity measurements and depolarization measurements of target media when illuminated with circularly polarized light.

Optical phantoms for biomedical polarimetry: a review

Calibration, quantification, and standardization of the polarimetric instrumentation, as well as interpretation and understanding of the obtained data, require the development and use of well-calibrated phantoms and standards. We reviewed the status of tissue phantoms for a variety of applications in polarimetry; more than 500 papers are considered. We divided the phantoms into five groups according to their origin (biological/nonbiological) and fundamental polarimetric properties of retardation, depolarization, and diattenuation. We found that, while biological media are generally depolarizing, retarding, and diattenuating, only one of all the phantoms reviewed incorporated all these properties, and few considered at least combined retardation and depolarization. Samples derived from biological tissue, such as tendon and muscle, remain extremely popular to quickly ascertain a polarimetric system, but do not provide quantifiable results aside from relative direction of their principal optical axis. Microspheres suspensions are the most utilized phantoms for depolarization, and combined with theoretical models can offer true quantification of depolarization or degree of polarization. There is a real paucity of birefringent phantoms despite the retardance being one of the most interesting parameters measurable with polarization techniques. Therefore, future work should be directed at generating truly reliable and repeatable phantoms for this metric determination. Diattenuating phantoms are rare and application-specific. Given that diattenuation is considered to be low in most biological tissues, the lack of such phantoms is seen as less problematic. The heterogeneity of the phantoms reviewed points to a critical need for standardization in this field. Ultimately, all research groups involved in polarimetric studies and instruments development would benefit from sharing a limited set of standardized polarimetric phantoms, as is done earlier in the round robin investigations in ellipsometry.

Characterizing the depolarization of circularly polarized light in turbid scattering media

We investigate the effectiveness of various bulk optical parameters in characterizing the degree of circular polarization (DOCP) of light diffusely reflected from scattering media. It is demonstrated that the traditional set of bulk optical parameters (namely, the scattering and absorption coefficients and the scattering asymmetry parameter) fail to characterize the observed depolarization. However, we find that there exists an additional parameter connected to the circular polarization memory phenomenon that consistently relates to observations, even in media with widely varying refractive indices and particle size distributions. This relationship is demonstrated using both Monte Carlo simulations and a new method for designing microsphere-based phantom media, which contain carefully controlled particle size distributions and depolarization characteristics.

Polarization-resolved second harmonic generation imaging of human ovarian cancer

Remodeling of the extracellular matrix in human ovarian cancer can be manifested in increased collagen concentration, changes in alignment within fibrils/fibers and/or up-regulation of different collagen isoforms. We used pixel-based second harmonic generation (SHG) polarization microscopy analyses to probe these molecular changes in human ovarian tissues [normal stroma, benign tumors, and high-grade serous (HGS) tumors] by: (i) determination of the α-helical pitch angle via the single-axis molecular model, (ii) collagen alignment within fibrils via SHG anisotropy, and (iii) chirality via SHG circular dichroism (SHG-CD). Pixel approaches are required due to the complex structure of the matrix that lacks a high degree of fiber alignment. The largest differences in the helical pitch angle were between normal stroma and benign tumors, consistent with gene expression showing the Col III isoform is up-regulated in the latter. The data were not consistent with up-regulation of Col III in HGS tumors as previous reports have suggested. The different tissues also displayed differing SHG anisotropies and SHG-CD responses, consistent with either Col III incorporation or randomization of Col I alignment within benign and malignant tumors. Additionally, the high-grade tumors displayed higher collagen concentration, where this desmoplasia is consistent with the higher fiber density in these tissues. These results collectively indicate that the fibril assemblies are distinct in all tissues, where these differences likely result from the synthesis of collagen rather than remodeling of existing collagen. Importantly, these analyses are label-free and interrogate subresolution collagen structure on intact tissues, without the need for conventional structural biology tools.

Assessment of tissue polarimetric properties using Stokes polarimetric imaging with circularly polarized illumination

Tissue-depolarization and linear-retardance are the main polarization characteristics of interest for bulk tissue characterization, and are normally interpreted from Mueller polarimetry. Stokes polarimetry can be conducted using simpler instrumentation and in a shorter time. Here, we use Stokes polarimetric imaging with circularly polarized illumination to assess the circular-depolarization and linear-retardance properties of tissue. Results obtained were compared with Mueller polarimetry in transmission and reflection geometry, respectively. It is found that circular-depolarization obtained from these 2 methods is very similar in both geometries, and that linear-retardance is highly quantitatively similar for transmission geometry and qualitatively similar for reflection geometry. The majority of tissue circular-depolarization and linear-retardance image information (represented by local image contrast features) obtained from Mueller polarimetry is well preserved from Stokes polarimetry in both geometries. These findings can be referred to for further understanding tissue Stokes polarimetric data, and for further application of Stokes polarimetry under the circumstances where short acquisition time or low optical system complexity is a priority, such as polarimetric endoscopy and microscopy.

Wavelength-Dependent Second Harmonic Generation Circular Dichroism for Differentiation of Col I and Col III Isoforms in Stromal Models of Ovarian Cancer Based on Intrinsic Chirality Differences

Extensive remodeling of the extracellular matrix (ECM) occurs in many epithelial cancers. For example, in ovarian cancer, upregulation of collagen isoform type III has been linked to invasive forms of the disease, and this change may be a potential biomarker. To examine this possibility, we implemented wavelength-dependent second harmonic generation circular dichroism (SHG-CD) imaging microscopy to quantitatively determine changes in chirality in ECM models comprised of different Col I/Col III composition. In these models, Col III was varied between 0 and 40%, and we found increasing Col III results in reduced net chirality, consistent with structural biology studies of Col I and III in tissues where the isoforms comingle in the same fibrils. We further examined the wavelength dependence of the SHG-CD to both optimize the response and gain insight into the underlying mechanism. We found using shorter SHG excitation wavelengths resulted in increased SHG-CD sensitivity, where this is consistent with the electric-dipole-coupled oscillator model suggested previously for the nonlinear chirality response from thin films. Moreover, the sensitivity is further consistent with the wavelength dependency of SHG intensity fit to a two-state model of the two-photon absorption in collagen. We also provide experimental calibration protocols to implement the SHG-CD modality on a laser scanning microscope. We last suggest that the technique has broad applicability in probing a wide range of diseased states with changes in collagen molecular structure.

Polarization influence on reflectance measurements in the spatial frequency domain

In this work, we quantify the influence of crossed polarizers on reflectance measurements in the spatial frequency domain. The use of crossed polarizers is a very common approach for suppression of specular surface reflections. However, measurements are typically evaluated using a non-polarized scalar theory. The consequences of this discrepancy are the focus of our study, and we also quantify the related errors of the derived optical properties. We used polarized Monte Carlo simulations for forward calculation of the reflectance from different samples. The samples' scatterers are assumed to be spherical, allowing for the calculation of the scattering functions by Mie theory. From the forward calculations, the reduced scattering coefficient [Formula: see text] and the absorption coefficient μa were derived by means of a scalar theory, as commonly used. Here, we use the analytical solution of the scalar radiative transfer equation. With this evaluation approach, which does not consider polarization, we found large errors in [Formula: see text] and μa in the range of 25% and above. Furthermore, we investigated the applicability of the use of a reference measurement to reduce these errors as suggested in literature. We found that this method is not able to generally improve the accuracy of measurements in the spatial frequency domain. Our general recommendation is to apply a polarized theory when using crossed polarizers.

Circular polarization memory in polydisperse scattering media

We investigate the survival of circularly polarized light in random scattering media. The surprising persistence of this form of polarization has a known dependence on the size and refractive index of scattering particles, however a general description regarding polydisperse media is lacking. Through analysis of Mie theory, we present a means of calculating the magnitude of circular polarization memory in complex media, with total generality in the distribution of particle sizes and refractive indices. Quantification of this memory effect enables an alternate pathway toward recovering particle size distribution, based on measurements of diffusing circularly polarized light.

Probing microstructural information of anisotropic scattering media using rotation-independent polarization parameters

Polarization parameters contain rich information on the micro- and macro-structure of scattering media. However, many of these parameters are sensitive to the spatial orientation of anisotropic media, and may not effectively reveal the microstructural information. In this paper, we take polarization images of different textile samples at different azimuth angles. The results demonstrate that the rotation insensitive polarization parameters from rotating linear polarization imaging and Mueller matrix transformation methods can be used to distinguish the characteristic features of different textile samples. Further examinations using both experiments and Monte Carlo simulations reveal that the residue rotation dependence in these polarization parameters is due to the oblique incidence illumination. This study shows that such rotation independent parameters are potentially capable of quantitatively classifying anisotropic samples, such as textiles or biological tissues.

Monte Carlo model of the depolarization of backscattered linearly polarized light in the sub-diffusion regime

We present a predictive model of the depolarization ratio of backscattered linearly polarized light from spatially continuous refractive index media that is applicable to the sub-diffusion regime of light scattering. Using Monte Carlo simulations, we derived a simple relationship between the depolarization ratio and both the sample optical properties and illumination-collection geometry. Our model was validated on tissue simulating phantoms and found to be in good agreement. We further show the utility of this model by demonstrating its use for measuring the depolarization length from biological tissue in vivo. We expect our results to aid in the interpretation of the depolarization ratio from sub-diffusive reflectance measurements.

Ex vivo assessment of carbon tetrachloride (CCl(4))-induced chronic injury using polarized light spectroscopy

The liver performs various functions, such as the production and detoxification of chemicals; therefore, it is susceptible to hepatotoxins such as carbon tetrachloride (CCl4), which causes chronic injury. Thus, assessment of injury and its status of severity are of prime importance. Current work reports an ex vivo study for probing the severance of hepatic injury induced by CCl4 with polarized light over the spectral range 400-800 nm. Different concentrations of CCl4 were used to induce varying severity of hepatic injury in a rat model. Linear retardance, depolarization rates, and diagonal Mueller matrix elements (m22, m33, and m44), were successfully used as the distinguishing criterion for normal and different liver injuries. Our results show that linear retardance for injured liver samples with lower doses of CCl4 tends to increase when compared with normal liver samples, while samples injured at higher doses of CCl4 offer almost no retardance. Total, linear, and circular depolarizations follow decreasing trends with increased liver injury severity over the entire investigated wavelength range. Linear polarization states were observed to be better maintained as compared to circular polarization states for all samples. Furthermore, numerical values of diagonal elements of the experimentally measured Mueller matrix also increase with increasing doses of CCl4. Liver fibroses, change in transport albedo, and the relative refractive index of the extracellular matrix caused by CCl4 are responsible for the observed differences. These results will provide a pathway to gauge the severity of injury caused by toxic chemicals.

Analytical light reflectance models for overlapping illumination and collection area geometries

Several biomedical applications, such as detection of dysplasia, require selective interrogation of superficial tissue structures less than a few hundred micrometers thick. Techniques and methods have been developed to limit the penetration depth of light in tissue, including the design of systems such as fiber-optic probes that have overlapping illumination and collection areas on the tissue surface. For such geometries, the diffusion approximation to the light-transport equation typically does not apply, and as a result there is no general model to extract tissue optical properties from reflectance measurements. In the current study, we employ Monte Carlo (MC) simulations to develop simple and compact analytical models for the light reflectance from these overlapping geometries. These models incorporate the size of the illumination and collection areas, the collection angle, the polarization of the incident light, and the optical properties of the sample. Moreover, these MC simulations use the Whittle-Matérn model to describe scattering from spatially continuous refractive index media such as tissue, which is more general than models based on the conventionally used Henyey-Greenstein model. We validated these models on tissue-simulating phantoms. The models developed herein will facilitate the extraction of optical properties and aid in the design of optical systems employing overlapping illumination and collection areas, including fiber-optic probes for in vivo tissue diagnosis.

Polarization sensitive optical coherence tomography of melanin provides intrinsic contrast based on depolarization

Polarization sensitive optical coherence tomography (PS-OCT) is a functional extension of OCT. In addition to imaging based on tissue reflectivity, PS-OCT also enables depth-resolved mapping of sample polarization properties such as phase-retardation, birefringent axis orientation, Stokes vectors, and degree of polarization uniformity (DOPU). In this study, PS-OCT was used to investigate the polarization properties of melanin. In-vitro measurements in samples with varying melanin concentrations revealed polarization scrambling, i.e. depolarization of backscattered light. Polarization scrambling in the PS-OCT images was more pronounced for higher melanin concentrations and correlated with the concentration of the melanin granules in the phantoms. Moreover, in-vivo PS-OCT was performed in the retinas of normal subjects and individuals with albinism. Unlike in the normal eye, polarization scrambling in the retinal pigment epithelium (RPE) was less pronounced or even not observable in PS-OCT images of albinos. These results indicate that the depolarizing appearance of pigmented structures like, for instance, the RPE is likely to be caused by the melanin granules contained in these cells.

Tomographic imaging with polarized light

We report three-dimensional tomographic reconstruction of optical parameters for the mesoscopic light scattering regime from experimentally obtained datasets by using polarized light. We present a numerically inexpensive approximation to the radiative transfer equation governing the polarized light transport. This approximation is employed in the reconstruction algorithm, which computes two optical parameters by using parallel and perpendicular polarizations of transmitted light. Datasets were obtained by imaging a scattering phantom embedding highly absorbing inclusions. Reconstruction results are presented and discussed.

Prehistological evaluation of benign and malignant pigmented skin lesions with optical computed tomography

Discrimination of benign and malignant melanocytic lesions is a major issue in clinical dermatology. Assessment of the thickness of melanoma is critical for prognosis and treatment selection. We aimed to evaluate a novel optical computed tomography (optical-CT) system as a tool for three-dimensional (3-D) imaging of melanocytic lesions and its ability to discriminate benign from malignant melanocytic lesions while simultaneously determining the thickness of invasive melanoma. Seventeen melanocytic lesions, one hemangioma, and normal skin were assessed immediately after their excision by optical-CT and subsequently underwent histopathological examination. Tomographic reconstructions were performed with a back-propagation algorithm calculating a 3-D map of the total attenuation coefficient (AC). There was a statistically significant difference between melanomas, dysplastic nevi, and non-dysplastic nevi, as indicated by Kruskal-Wallis test. Median AC values were higher for melanomas compared with dysplastic and non-dysplastic nevi. No statistically significant difference was observed when thickness values obtained by optical-CT were compared with histological thickness using a Wilcoxon sighed rank test. Our results suggest that optical-CT can be important for the immediate prehistological evaluation of biopsies, assisting the physician for a rapid assessment of malignancy and of the thickness of a melanocytic lesion.

Quantitative correlation between light depolarization and transport albedo of various porcine tissues

We present a quantitative study of depolarization in biological tissues and correlate it with measured optical properties (reduced scattering and absorption coefficients). Polarized light imaging was used to examine optically thick samples of both isotropic (liver, kidney cortex, and brain) and anisotropic (cardiac muscle, loin muscle, and tendon) pig tissues in transmission and reflection geometries. Depolarization (total, linear, and circular), as derived from polar decomposition of the measured tissue Mueller matrix, was shown to be related to the measured optical properties. We observed that depolarization increases with the transport albedo for isotropic and anisotropic tissues, independent of measurement geometry. For anisotropic tissues, depolarization was higher compared to isotropic tissues of similar transport albedo, indicating birefringence-caused depolarization effects. For tissues with large transport albedos (greater than ~0.97), backscattering geometry was preferred over transmission due to its greater retention of light polarization; this was not the case for tissues with lower transport albedo. Preferential preservation of linearly polarized light over circularly polarized light was seen in all tissue types and all measurement geometries, implying the dominance of Rayleigh-like scattering. The tabulated polarization properties of different tissue types and their links to bulk optical properties should prove useful in future polarimetric tissue characterization and imaging studies.