Experimental validation of an extended Jones matrix calculus model to study the 3D structural orientation of the collagen fibers in articular cartilage using polarization-sensitive optical coherence tomography

Optical anisotropy composition of benign and malignant prostate tissues revealed by Mueller-matrix imaging

A Mueller matrix imaging approach is employed to disclose the three-dimensional composition framework of optical anisotropy within cancerous biotissues. Visualized by the Mueller matrix technique spatial architecture of optical anisotropy of tissues is characterised by high-order statistical moments. Thus, quantitative analysis of the spatial distribution of optical anisotropy, such as linear and circular birefringence and dichroism, is revealed by using high-order statistical moments, enabling definitively discriminate prostate adenoma and carcinoma. The developed approach provides greater (>90%) accuracy of diagnostic achieved by using either the 3-rd or 4-th order statistical moments of the linear anisotropy parameters. Noticeable difference is observed between prostate adenoma and carcinoma tissue samples in terms of the extinction coefficient and the degree of depolarisation. Juxtaposition to other optical diagnostic modalities demonstrates the greater accuracy of the approach described herein, paving the way for its wider application in cancer diagnosis and tissue characterization.

Differentiation of Human GBM From Non-GBM Brain Tissue With Polarization Imaging Technique

As for optical techniques, it is difficult for the 5-aminolevulinic (5-ALA) fluorescence guidance technique to completely detect glioma due to residual cells in the blind area and the dead angle of vision under microscopy. The purpose of this research is to characterize different microstructural information and optical properties of formalin-soaked unstained glioblastoma (GBM) and non-GBM tissue with the polarization imaging technique (PIT), and provide a novel method to detect GBM during surgery. In this paper, a 3×3 Mueller matrix polarization experimental system in backscattering mode was built to detect the GBM and non-GBM tissue bulk. The Mueller matrix decomposition and transformation parameters of GBM and non-GBM tissue were calculated and analyzed, and showed that parameters (1−Δ) and t are good indicators for distinguishing GBM from non-GBM tissues. Furthermore, the central moment coefficients (CMCs) of the frequency distribution histogram (FDH) were also calculated and used to distinguish the cancerous tissues. The results of the experiments confirmed the feasibility of PIT applied in the clinic to detect glioma, laying the foundation for the subsequent non-invasive, non-staining glioma detection.

Mueller matrix polarimetry and polar decomposition of articular cartilage imaged in reflectance

Articular cartilage birefringence relates to zonal architecture primarily of type II collagen, which has been assessed extensively in transmission, through thin tissue sections, to evaluate cartilage repair and degeneration. Mueller matrix imaging of articular cartilage in reflection is of potential utility for non-destructive imaging in clinical and research applications. Therefore, such an imaging system was constructed to measure laser reflectance signals, calibrated, and tested with optical standards. Polar decomposition was chosen as a method to extract fundamental optical parameters from the experimental Mueller matrices, with performance confirmed by simulations. Adult bovine articular cartilage from the patellofemoral groove was found to have ∼0.93 radians retardance, low diattenuation of ∼0.2, and moderately high depolarization of 0.66. Simulations showed that variation in depolarization drives inaccuracy of depolarization and retardance maps derived by polar decomposition. These results create a basis for further investigation of the clinical utility of polarized signals from knee tissue and suggest potential approaches for improving the accuracy of polar decomposition maps.

Embossed topographic depolarisation maps of biological tissues with different morphological structures

Layered topographic maps of the depolarisation due to diffuse biological tissues are produced using a polarisation-holographic Mueller matrix method approach. Histological sections of myocardial tissue with a spatially structured optically anisotropic fibrillar network, and parenchymal liver tissue with a polycrystalline island structure are successfully mapped. The topography of the myocardium maps relates to the scattering multiplicity within the volume and the specific morphological structures of the biological crystallite networks. The overall depolarisation map is a convolution of the effects of these two factors. Parenchymal liver tissues behave broadly similarly, but the different biological structures present cause the degree of scattering multiplicity to increase more rapidly with increasing phase. Through statistical analysis, the dependences of the magnitudes of the first to fourth order statistical moments are determined. These moments characterise the changing distributions of the depolarisation values through the volume of biological tissues with different morphological structures. Parenchymal liver tissue depolarisation maps are characterised by larger mean and variance, and less skewness and kurtosis, compared to the distributions for the myocardium. This work demonstrates that a polarisation-holographic Mueller matrix method can be applied to the assessment of the 3D morphology of biological tissues, with applications in disease diagnosis.

Polarization-sensitive optical coherence tomography with a conical beam scan for the investigation of birefringence and collagen alignment in the human cervix

By measuring the phase retardance of a cervical extracellular matrix, our in-house polarization-sensitive optical coherence tomography (PS-OCT) was shown to be capable of (1) mapping the distribution of collagen fibers in the non-gravid cervix, (2) accurately determining birefringence, and (3) measuring the distinctive depolarization of the cervical tissue. A conical beam scan strategy was also employed to explore the 3D orientation of the collagen fibers in the cervix by interrogating the samples with an incident light at 45° and successive azimuthal rotations of 0-360°. Our results confirmed previous observations by X-ray diffraction, suggesting that in the non-gravid human cervix collagen fibers adjacent to the endocervical canal and in the outermost areas tend to arrange in a longitudinal fashion whereas in the middle area they are oriented circumferentially. PS-OCT can assess the microstructure of the human cervical collagen in vitro and holds the potential to help us better understand cervical remodeling prior to birth pending the development of an in vivo probe.

Measurement of biaxial optical birefringence in articular cartilage

Optical birefringence is a valuable parameter in evaluating collagen organization and assessing collagen degeneration in articular cartilage. A uniaxial birefringent model is implied when using existing methods such as polarized light microscopy in characterizing birefringence in cartilage. However, some studies suggest the existence of a sheet-like collagen organization in articular cartilage, which requires a biaxial birefringence model to describe. In this study, we applied a multi-incident birefringence measurement procedure to investigate the biaxial birefringence in articular cartilage. The results supported the existence of a small yet significant biaxial birefringence effect, which was in agreement with the expectation from a sheet-like collagen organization.

Slope-based segmentation of articular cartilage using polarization-sensitive optical coherence tomography phase retardation image

A segmentation method based on phase retardation measurements from polarization-sensitive optical coherence tomography (PS-OCT) is developed to differentiate the structural zones of articular cartilage. The organization of collagen matrix in articular cartilage varies over the different structural zones, generating different tissue birefringence. Analyzing the slope of the accumulated phase retardation at different depths can detect the variation in tissue birefringence and be used to segment the structural zones. The method is validated on phantoms composed of layers of different materials. Articular cartilage samples from adult swine are segmented with the method. The characteristics in each segmented zone are also examined by histology and high-resolution second-harmonic generation imaging, showing distinctive properties that match with the anatomical structure of articular cartilage. The segmentation algorithm is also applied on PS-OCT images acquired at multiple illumination angles, where the angular dependence of tissue birefringence in the deep zone is detected. This method offers a noninvasive imaging approach to differentiating the structural zones of articular cartilage, as well as a quantification approach based on the phase retardation measurements of PS-OCT. This method has great potential in studying depth-related progression of cartilage degeneration.

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.

Optic axis mapping with catheter-based polarization-sensitive optical coherence tomography

Birefringence offers an intrinsic contrast mechanism related to the microstructure and arrangement of fibrillary tissue components. Here we present a reconstruction strategy to recover not only the scalar amount of birefringence but also its optic axis orientation as a function of depth in tissue from measurements with catheter-based polarization sensitive optical coherence tomography. A polarization symmetry constraint, intrinsic to imaging in the backscatter direction, facilitates the required compensation for wavelength-dependent transmission through system elements, the rotating catheter, and overlying tissue layers. Applied to intravascular imaging of coronary atherosclerosis in human patients, the optic axis affords refined interpretation of plaque architecture.

Mapping 3D fiber orientation in tissue using dual-angle optical polarization tractography

Optical polarization tractography (OPT) has recently been applied to map fiber organization in the heart, skeletal muscle, and arterial vessel wall with high resolution. The fiber orientation measured in OPT represents the 2D projected fiber angle in a plane that is perpendicular to the incident light. We report here a dual-angle extension of the OPT technology to measure the actual 3D fiber orientation in tissue. This method was first verified by imaging the murine extensor digitorum muscle placed at various known orientations in space. The accuracy of the method was further studied by analyzing the 3D fiber orientation of the mouse tibialis anterior muscle. Finally we showed that dual-angle OPT successfully revealed the unique 3D "arcade" fiber structure in the bovine articular cartilage.

Polarization-sensitive optical coherence tomography-based imaging, parameterization, and quantification of human cartilage degeneration

Polarization-sensitive optical coherence tomography (PS-OCT) is a light-based, high-resolution, real-time, noninvasive, and nondestructive imaging modality yielding quasimicroscopic cross-sectional images of cartilage. As yet, comprehensive parameterization and quantification of birefringence and tissue properties have not been performed on human cartilage. PS-OCT and algorithm-based image analysis were used to objectively grade human cartilage degeneration in terms of surface irregularity, tissue homogeneity, signal attenuation, as well as birefringence coefficient and band width, height, depth, and number. Degeneration-dependent changes were noted for the former three parameters exclusively, thereby questioning the diagnostic value of PS-OCT in the assessment of human cartilage degeneration.

Deep tissue volume imaging of birefringence through fibre-optic needle probes for the delineation of breast tumour

Identifying tumour margins during breast-conserving surgeries is a persistent challenge. We have previously developed miniature needle probes that could enable intraoperative volume imaging with optical coherence tomography. In many situations, however, scattering contrast alone is insufficient to clearly identify and delineate malignant regions. Additional polarization-sensitive measurements provide the means to assess birefringence, which is elevated in oriented collagen fibres and may offer an intrinsic biomarker to differentiate tumour from benign tissue. Here, we performed polarization-sensitive optical coherence tomography through miniature imaging needles and developed an algorithm to efficiently reconstruct images of the depth-resolved tissue birefringence free of artefacts. First ex vivo imaging of breast tumour samples revealed excellent contrast between lowly birefringent malignant regions, and stromal tissue, which is rich in oriented collagen and exhibits higher birefringence, as confirmed with co-located histology. The ability to clearly differentiate between tumour and uninvolved stroma based on intrinsic contrast could prove decisive for the intraoperative assessment of tumour margins.

Polarized light interaction with tissues

This tutorial-review introduces the fundamentals of polarized light interaction with biological tissues and presents some of the recent key polarization optical methods that have made possible the quantitative studies essential for biomedical diagnostics. Tissue structures and the corresponding models showing linear and circular birefringence, dichroism, and chirality are analyzed. As the basis for a quantitative description of the interaction of polarized light with tissues, the theory of polarization transfer in a random medium is used. This theory employs the modified transfer equation for Stokes parameters to predict the polarization properties of single- and multiple-scattered optical fields. The near-order of scatterers in tissues is accounted for to provide an adequate description of tissue polarization properties. Biomedical diagnostic techniques based on polarized light detection, including polarization imaging and spectroscopy, amplitude and intensity light scattering matrix measurements, and polarization-sensitive optical coherence tomography are described. Examples of biomedical applications of these techniques for early diagnostics of cataracts, detection of precancer, and prediction of skin disease are presented. The substantial reduction of light scattering multiplicity at tissue optical clearing that leads to a lesser influence of scattering on the measured intrinsic polarization properties of the tissue and allows for more precise quantification of these properties is demonstrated.

3D Human cartilage surface characterization by optical coherence tomography

Early diagnosis and treatment of cartilage degeneration is of high clinical interest. Loss of surface integrity is considered one of the earliest and most reliable signs of degeneration, but cannot currently be evaluated objectively. Optical Coherence Tomography (OCT) is an arthroscopically available light-based non-destructive real-time imaging technology that allows imaging at micrometre resolutions to millimetre depths. As OCT-based surface evaluation standards remain to be defined, the present study investigated the diagnostic potential of 3D surface profile parameters in the comprehensive evaluation of cartilage degeneration. To this end, 45 cartilage samples of different degenerative grades were obtained from total knee replacements (2 males, 10 females; mean age 63.8 years), cut to standard size and imaged using a spectral-domain OCT device (Thorlabs, Germany). 3D OCT datasets of 8  ×  8, 4  ×  4 and 1  ×  1 mm (width  ×  length) were obtained and pre-processed (image adjustments, morphological filtering). Subsequent automated surface identification algorithms were used to obtain the 3D primary profiles, which were then filtered and processed using established algorithms employing ISO standards. The 3D surface profile thus obtained was used to calculate a set of 21 3D surface profile parameters, i.e. height (e.g. Sa), functional (e.g. Sk), hybrid (e.g. Sdq) and segmentation-related parameters (e.g. Spd). Samples underwent reference histological assessment according to the Degenerative Joint Disease classification. Statistical analyses included calculation of Spearman's rho and assessment of inter-group differences using the Kruskal Wallis test. Overall, the majority of 3D surface profile parameters revealed significant degeneration-dependent differences and correlations with the exception of severe end-stage degeneration and were of distinct diagnostic value in the assessment of surface integrity. None of the 3D surface profile parameters investigated were capable of reliably differentiating healthy from early-degenerative cartilage, while scan area sizes considerably affected parameter values. In conclusion, cartilage surface integrity may be adequately assessed by 3D surface profile parameters, which should be used in combination for the comprehensive and thorough evaluation and overall improved diagnostic performance. OCT- and image-based surface assessment could become a valuable adjunct tool to standard arthroscopy.

Detecting early stage osteoarthritis by optical coherence tomography?

Osteoarthritis (OA) is the most common chronic disease of our joints, manifested by a dynamically increasing degeneration of hyaline articular cartilage (AC). While currently no therapy can reverse this process, the few available treatment options are hampered by the inability of early diagnosis. Loss of cartilage surface, or extracellular matrix (ECM), integrity is considered the earliest sign of OA. Despite the increasing number of imaging modalities surprisingly few imaging biomarkers exist. In this narrative review, recent developments in optical coherence tomography are critically evaluated for their potential to assess different aspects of AC quality as biomarkers of OA. Special attention is paid to imaging surface irregularities, ECM organization and the evaluation of posttraumatic injuries by light-based modalities.

What can biophotonics tell us about the 3D microstructure of articular cartilage?

Connective tissues such as articular cartilage have been the subject of study using novel optical techniques almost since the invention of polarized light microscopy (PLM). Early studies of polarized light micrographs were the main evidential basis for the establishment of quantitative models of articular cartilage collagen structure by Benninghoff and others. Even now, state of the art optical techniques including quantitative polarized light microscopy (qPLM), optical coherence tomography (OCT), polarization-sensitive optical coherence tomography (PS-OCT), second harmonic generation (SHG) microscopy, Fourier-transform infrared (FTIR) microscopy, Raman and optical hyperspectral reflectance and fluorescence imaging are providing new insights into articular cartilage structure from the nanoscale through to the mesoscale. New insights are promised by emerging modalities such as optical elastography. This short review highlights some key recent results from modern optical techniques.

Bayesian maximum likelihood estimator of phase retardation for quantitative polarization-sensitive optical coherence tomography

This paper presents the theory and numerical implementation of a maximum likelihood estimator for local phase retardation (i.e., birefringence) measured using Jones-matrix-based polarization sensitive optical coherence tomography. Previous studies have shown conventional mean estimations of phase retardation and birefringence are significantly biased in the presence of system noise. Our estimator design is based on a Bayes' rule that relates the distributions of the measured birefringence under a particular true birefringence and the true birefringence under a particular measured birefringence. We used a Monte-Carlo method to calculate the likelihood function that describes the relationship between the distributions and numerically implement the estimator. Our numerical and experimental results show that the proposed estimator was asymptotically unbiased even with low signal-to-noise ratio and/or for the true phase retardations close to the edge of the measurement range. The estimator revealed detailed clinical features when applied to the in vivo anterior human eye.

Conical scan polarization-sensitive optical coherence tomography

We report on a new articular cartilage imaging technique with potential for clinical arthroscopic use, by supplementing the variable-incidence-angle polarization-sensitive optical coherence tomography method previously developed by us with a conical beam scan protocol. The technique is validated on bovine tendon by comparing experimental data with simulated data generated using the extended Jones matrix calculus. A unique capability of this new optical technique is that it can locate the "brushing direction" of collagen fibers in articular cartilage, which is structural information that extends beyond established methods such as split-line photography or birefringent fast-axis measurement in that it is uniquely defined over the full azimuthal-angle range of (-π, + π). The mapping of this direction over the cartilage surface may offer insights into the optimal design of tissue-engineering scaffolds for cartilage repair.

Spectral binning for mitigation of polarization mode dispersion artifacts in catheter-based optical frequency domain imaging

Polarization mode dispersion (PMD) has been recognized as a significant barrier to sensitive and reproducible birefringence measurements with fiber-based, polarization-sensitive optical coherence tomography systems. Here, we present a signal processing strategy that reconstructs the local retardation robustly in the presence of system PMD. The algorithm uses a spectral binning approach to limit the detrimental impact of system PMD and benefits from the final averaging of the PMD-corrected retardation vectors of the spectral bins. The algorithm was validated with numerical simulations and experimental measurements of a rubber phantom. When applied to the imaging of human cadaveric coronary arteries, the algorithm was found to yield a substantial improvement in the reconstructed birefringence maps.

Comparison of different metrics for analysis and visualization in spectroscopic optical coherence tomography

Spectroscopic Optical Coherence Tomography (S-OCT) extracts depth resolved spectra that are inherently available from OCT signals. The back scattered spectra contain useful functional information regarding the sample, since the light is altered by wavelength dependent absorption and scattering caused by chromophores and structures of the sample. Two aspects dominate the performance of S-OCT: (1) the spectral analysis processing method used to obtain the spatially-resolved spectroscopic information and (2) the metrics used to visualize and interpret relevant sample features. In this work, we focus on the second aspect, where we will compare established and novel metrics for S-OCT. These concepts include the adaptation of methods known from multispectral imaging and modern signal processing approaches such as pattern recognition. To compare the performance of the metrics in a quantitative manner, we use phantoms with microsphere scatterers of different sizes that are below the system's resolution and therefore cannot be differentiated using intensity based OCT images. We show that the analysis of the spectral features can clearly separate areas with different scattering properties in multi-layer phantoms. Finally, we demonstrate the performance of our approach for contrast enhancement in bovine articular cartilage.