As one of the most salient features of light, polarization can be used to develop imaging techniques capable of probing the structural and optical properties of media.12.–3 Polarization imaging methods have the ability to measure the microstructural information nondestructively in tissues, cells, and other specimens.4,5 Polarization techniques can especially help suppress the multiscattered photons from deep tissues, and thus improve the imaging contrast of superficial layers of samples.6 For this reason, in the past two decades, several techniques such as the degree of polarization and difference polarization have been applied to the detection of human skin cancers in vivo.7,8 Recently, as comprehensive descriptions of polarization properties, the Mueller matrix imaging techniques have been regarded as potential methods of cancerous tissues’ diagnosis9 and preliminarily applied to the detection of colon cancer,1011.–12 thyroid cancer,13 cervical cancer,14,15 and so on.1617.–18 Through calculations using the Mueller matrix polar decomposition (MMPD) and Mueller matrix transformation (MMT), groups of quantitative polarization imaging parameters with explicit physics meanings can be obtained.19,20
Nowadays the increasing cancer incidence rate is becoming one of the biggest threats to human health. Among all types of cancers, liver cancer ranks in the top five in frequency with an estimated number of over 700,000 new cases a year all over the world. In high-risk countries including China, liver cancer can arise before the age of 20 years.21,22 According to the mortality rate in the United States, liver cancer will become the third leading cause of cancer-related death by 2030.23 In China, there are over 350,000 new cases of liver cancer each year, accounting for more than 50% of those in the world.24 During the development of liver hepatitis to cancer, fibrosis often evolves as part of a healing process in response to liver damage, resulting in cirrhosis of liver tissues.25 The diagnosis of liver fibrosis is usually based on histological findings of liver biopsy observed by optical microscope.26 This procedure needs stained tissue sections and experienced pathologists to result in accurate evaluations.27 Recently, the second-harmonic generation (SHG) technique has been applied for quantitative assessment of fibrosis in chronic hepatitis B patients.28 The results and analysis have shown that the fibrosis scoring accuracy can be improved by SHG images and imaging processing methods. On the other hand, as an optical technique with simple structure and low cost, the polarization imaging can also be suitable for the detection of fibrous structures. Compared with the traditional optical microscope and SHG imaging, the Mueller matrix microscope and transformed parameters such as MMPD and MMT may provide more information on unstained fibrous structures whose orders of alignment are not so good. Therefore, the transformed Mueller matrix parameters can possibly be used as good indicators for pathological diagnosis of liver cancer. Recently, in order to extract microstructural information of pathological tissue slices more effectively, we designed a Mueller matrix microscope by adding both the polarization states generator and analyzer (PSG and PSA) to an ordinary commercial optical microscope.25
In this paper, based on the methods proposed in a previous study,25 we take it one step further to quantitatively facilitate the fibrosis detection and scorings of pathological liver tissue samples in different stages from cirrhosis to cancer using the Mueller matrix microscope. We measure and analyze the MMPD and MMT parameters for the fibrotic liver tissue samples in different stages. Then we conduct the Monte Carlo (MC) simulations based on the sphere birefringence model to examine in detail the influence of structural changes in different fibrosis stages on the imaging parameters. Both the experimental and MC simulated results indicate that the polarized light microscope and transformed Mueller matrix parameters have a good application prospect in liver cancer diagnosis and scoring.
Methods and Materials
Figure 1 shows the polarization Mueller matrix microscope used in this study. The microscope was transformed from a commercial transmission microscope (Yuexian L2050, Guangzhou, China) by adding the PSG and PSA. The light source is an LED (3 W, 632 nm, ). The polarization states of the incident light are controlled by the PSG consisting of a set of a linear polarizer (P1, extinction ratio 500:1, Daheng Optics, China) and a quarter-wave plate (R1, Daheng Optics). The polarized light passes through the tissue sample on the stage and the objective lens, and then the PSA, consisting of another set of a quarter-wave plate (R2, Daheng Optics) and a linear polarizer (P2, extinction ratio 500:1, Daheng Optics). Both the PSG and PSA have been designed to be compact modules, as shown in Fig. 1, to be incorporated into the microscope. The emitted light is collected by a 12-bit CCD camera (QImaging 74-0107A, Canada) to produce the polarization images of the sample. For the Mueller matrix measurement, 30 images with specific incident and output polarization states are obtained. Using 30 images, we can estimate the Fourier coefficients needed to compute the Mueller matrix elements.
During the experiments, the polarizers (P1 and P2) are fixed in the 0-deg direction. The two quarter-wave plates (R1 and R2) are rotated harmonically to generate the polarized light. In addition, the rotation rates of the quarter-wave plates R1 and R2 are . Therefore, the intensity of light can be given by the following equation:29. More details for this Mueller matrix imaging method can be found in Refs. 30 and 31. To calibrate the microscope, we measured some standard samples including air, a quarter-wave plate, and a polarizer. We compensated the errors due to nonideal quarter-wave plates and polarizers as described in Ref. 32. The experimental results testified that the maximum errors for the absolute values of the Mueller matrix elements are about 0.01.
Mueller Matrix Polar Decomposition and Mueller Matrix Transformation Parameters
Mueller matrices contain abundant structural and optical information of samples. However, since the Mueller matrix elements lack clear physics meanings and relations to certain microstructures, there are huge difficulties in applying Mueller matrix polarimetry to biomedical studies, especially for quantitative analysis. Recently, several approaches to transform Mueller matrix elements to quantitative parameters with explicit physics meanings have been proposed. In this paper, we adopt the parameters of both MMPD and MMT techniques to analyze the microscopic imaging results of the liver tissue samples. The MMPD method is developed by Lu and Chipman19 and has been widely applied to biomedical studies.1415.16.17.–18 It is based on the three main interactions between the polarized light and media: diattenuation (), retardation (), and depolarization (). In MMPD, a Mueller matrix () can be decomposed and the linear retardance () and its orientation angle () can be derived by the following equations:33
In previous studies, by fitting the Mueller matrix elements into certain trigonometric functions, we proposed the MMT method. The MMT method can provide a set of parameters which is insensitive to the azimuth angle for the backscattering imaging of bulk tissues20 or transmission imaging of thin tissues with limited scattering such as the pathological sections of liver tissues used in this paper.25 Therefore, the MMT parameters and suitable for transmission microscopic imaging can be acquired by the following equations:
Experiments and simulations have shown that the MMT parameter and MMPD parameter are good indicators of the retardance of the media, while the MMT parameter and MMPD parameter are related to the orientation directions of the aligned fibrous structures.25,3435.–36 The MMPD parameters are slightly more sensitive to the fibers, whereas the MMT parameters can be calculated more easily and quickly.37 In our recent studies, we have found that, for the thin tissue slices with anisotropic fibrous structures, the birefringence effect plays the dominant role for the polarization imaging contrast mechanism;37 therefore, for the cancerous liver tissue sections, the values of parameters and and the distribution of parameters and may be used to reflect the density and orientation of the fibrous structure.
Liver Tissue Samples
Polarization techniques can effectively improve the imaging quality and contrast of superficial tissues. Since most cancers occur in the superficial epithelium throughout the body in the early stage,38 polarization imaging techniques have been regarded as potential methods of cancerous tissues diagnosis. In recent studies, we have found that the cancerous liver tissues contain abundant fibrous structures, which can be detected by polarization parameters.25 In this study, we choose some human liver cirrhosis and cancer tissues with fibrosis as the samples for the polarized Mueller matrix microscope. As shown in Fig. 2, the samples are nonstained, dewaxing sections of pathological liver tissue slices, which are prepared and provided by Shenzhen Sixth People’s (Nanshan) Hospital. For the liver cancer tissues, there are always inflammatory reactions accompanied by fibrosis formations. In addition, the pathological liver tissues of different stages have different proportions of fibrous structures, which can be used as the indicators for stage scorings.39 We choose four -thick slices of human liver tissues in different fibrosis stages (F1 to F4) as the imaging samples, as shown in Figs. 2(a)–2(d). For comparisons, the corresponding -thick hematoxylin–eosin (H–E) stained slices shown as Figs. 3(a)–3(d) are also provided. This work was approved by the Ethics Committee of the Shenzhen Sixth People’s (Nanshan) Hospital.
Results and Discussion
Microscopic Imaging Results of Liver Tissue Samples
Figure 2(a) shows the intensity image of the -thick nonstained slice of the liver tissue with fibrosis in F1 stage under objective observation. Figures 4(a)–4(d) are the images of the MMPD and MMT parameters , , , and . Previous studies have shown that and are related to the retardance and and are related to the angle orientation of the fibrous structures. From the intensity image shown as Fig. 2(a), the fibrous structures around the liver cells can hardly be discriminated. However, from the images of the parameters and shown in Fig. 4, we can see that there are blurred circular regions with relatively high values, indicating the existence of fibers. As mentioned above, the fibrosis process is closely related to the inflammatory reaction. It can be observed from Figs. 4(a) and 4(b) that the values of parameters and for the fibrotic circular regions are larger than the other parts, which means that the liver tissues with fibrosis have more prominent retardance than the normal tissues. Moreover, we can see that the lower right regions have the largest values of and , indicating the densest fibers are distributed here. By implementing the image analysis and processing method used by Dubreuil et al.,39 it is found that the largest values of parameters and of the entire image are 0.41 and 0.38, respectively, and the proportion of fibrous distribution is about 6%. In addition to the density, from the images of the parameters and shown as Figs. 4(c) and 4(d), the alignment orientations of the circular aligned fibers can be roughly revealed. The microscopic imaging results indicate that in F1 stage there are fibers existing in the liver tissues, but the fibrosis is not particularly evident, which is in accordance with the pathological description.
For comparisons, we then apply the Mueller matrix microscopy to the liver tissue with fibrosis in the F2 stage, whose intensity image under objective observation is shown as Fig. 2(b). Figures 5(a)–5(d) are the images of the MMPD and MMT parameters , , , and . Some circularly aligned fibrous structures around the liver cells can be vaguely seen in Fig. 2(b). Figures 5(a) and 5(b) confirm that the values of parameters and for the fibrotic circular regions are larger than the other parts, and the lower right regions have the largest values, meaning that the fibrosis degree of these regions is probably more obvious. We also notice that compared with the liver tissue in the F1 stage, the retardance of the sample in the F2 stage can be more prominent. Calculations show that the largest values of parameters and of the entire image are 0.78 and 0.75, respectively, and the proportion of fibrous distribution is about 10%. From the images of the parameters and shown as Figs. 5(c) and 5(d), the orientations of the fibrous structures can be observed: the structural alignments for the upper left and lower right fibrotic regions are much clearer than the other parts. The results indicate that in the F2 stage the fibrosis in liver tissue becomes evident in some parts, which can also be observed in the pathological H–E stained section.
Figure 2(c) shows the intensity image of the liver tissue with cirrhosis in F3 stage under objective observation. Figures 6(a)–6(d) are the images of the MMPD and MMT parameters , , , and . We can clearly see that there are some circularly distributed fibers in Fig. 2(c). The imaging results shown in Figs. 6(a) and 6(b) demonstrate that the parameters and for the fibrotic circular regions have larger values than the other parts, indicating the existence of well-ordered fibers. In addition, there are some regions with relatively lower values surrounding the circular regions, meaning that the distribution range of fibrosis in the F3 stage may become wider than those in the F1 and F2 stages. Calculations show that the largest values of parameters and of the entire image are 0.82 and 0.80, respectively, and the proportion of fibrous distribution is about 12%. It can also be observed from the images of parameters and shown as Figs. 6(c) and 6(d) that the circularly aligned fibrotic regions have very clear distribution orientations. The results shown in Fig. 6 indicate that compared with the F1 and F2 stages, in the F3 stage the fibers in liver tissues are with a wider distribution and the fibrosis becomes more obvious in some regions.
Finally, we choose the -thick nonstained slice of cancerous liver tissue in the F4 stage, whose intensity image under objective observation is shown as Fig. 2(d). In both the unstained Fig. 2(d) and H–E stained Fig. 3(d), we can see that there are abundant fibrous structures around the liver cells. Figures 7(a) and 7(b) show that the values of parameters and for the cancerous liver sample are much larger than the samples in the F1 to F3 stages. Moreover, the boundaries of the fibrotic regions are clear, demonstrating that the fibrosis degree of these regions is very high. Calculations show that the largest values of parameters and of the entire image are 0.96 and 0.95, respectively, and the proportion of fibrous distribution is about 16%. The images of parameters and shown as Figs. 7(c) and 7(d) also reveal that there are a large portion of fibers with distinct orientations in the cancerous liver tissue. The results shown in Figs. 4–7 testify that: (a) both the microscopic MMPD and MMT parameters can be used as the potential indicators for the fibrosis process in pathological liver tissue samples. (b) As the development of cancer-induced inflammatory reaction increases, the proportion of fibrous structures in liver tissues becomes larger. In the F4 stage, the fibers are much denser than in the F1 to F3 stages, which can be clearly observed using the Mueller matrix microscope.
Quantitative Analysis of Imaging Results
For a more detailed analysis of the characteristic variations of parameters , , , and in different stages, we make quantitative studies in this section. First, from Fig. 8, which shows the average values over the entire images of the parameters and in different fibrosis stages, we can see that both the parameters increase as the stage changes from F1 to F4. This is because a higher development stage means a larger amount of fibers in the liver tissues, leading to a more prominent retardance effect. In addition, it is also shown that from F1 to F4, the mean value of parameter changes more prominently than that of parameter , indicating that parameter may be more sensitive to the variations of fibrous structures. Although more statistical studies are still needed, the preliminary results demonstrate that the parameters and may be used as quantitative tools for the detection and scorings of liver cirrhosis and cancer. Second, from Fig. 9, which shows the distribution histograms of the parameter in different stages, we can see that: (a) The values of parameter have a relatively uniform distribution in the F1 stage. (b) The distribution of the parameter begins to converge to certain values in the F2 and F3 stages. (c) In the F4 stage, the convergence of parameter becomes prominent, showing abundant fibrous structures distributed in a relatively narrow range of orientations. The characteristic distribution feature of parameter is very similar to parameter . The results shown in Fig. 9 indicate that, although more statistical analyses are still needed, the distribution feature of the parameter or may offer information to help determine the fibrosis degree. Since the pathological liver tissues in different stages have different proportions of fibrous microstructures, the results shown in Figs. 8 and 9 confirm that the MMPD and MMT parameters can provide additional quantitative information helpful for the accurate fibrosis scorings of liver cirrhosis and cancers.
In order to interpret the experimental results above, we carry on MC simulations based on the sphere birefringence model to analyze the relationship between the parameters and the fibrosis process of the liver tissues.40 In the MC simulations, there are two types of spherical scatterers, the large ones ( in diameter, scattering coefficient is ) represent the cell nuclei and the small ones ( in diameter, scattering coefficient is ) represent the organelles. According to the liver tissue samples, the thickness of the medium is and the refractive indices of the interstitial medium and scatterers are 1.33 and 1.45, respectively.37 As mentioned in Sec. 2.2, for the thin tissue slices with anisotropic fibrous structures, the birefringence effect plays the dominant role for the polarization imaging contrast mechanism.37 In addition, the experimental results have shown that from stage F1 to F4 not only the values of retardance of fibrotic regions, but also the proportion of fibrous distribution increase. Therefore, to quantitatively analyze the fibrosis process during the F1 to F4 stages, in the MC simulations, both the increases of birefringence value and proportion of fibrous distribution are considered. As shown in Fig. 10, in the simulations, the value of is increased from 0.001 to 0.0025 and the optical axis is along the -axis direction. Meanwhile, the proportion of fibrous distribution is changed from 6%, 10%, 12%, to 16% according to the experimental results. The simulated results shown in Fig. 10 demonstrate that the values of parameters and increase as the birefringence increases, which are consistent with the experimental observations. It should be pointed out that in Fig. 10, the absolute values of the parameters and are slightly lower than those in Fig. 8. In the MC simulation, there are several constant parameters used such as the thickness of the medium and the diameter of the scatterers. Actually, these parameters can be varied in different regions of the tissue samples. Since such variations are hard to estimate, we use approximately constant values of these parameters, which may influence the results of the MC simulations. It can also be observed that when changes from 0.001 to 0.0025, the mean value of parameter changes about 0.02, while the mean value of parameter changes only about 0.01, meaning that the MMPD parameter is slightly more sensitive to the variations of . Figure 11 shows that both the parameters and can reflect the orientation of the optical axis or the alignment direction of the fibers. It should be pointed out that in MC simulations, when the value of is changed from 0.001 to 0.0025, or the proportion of fibrous distribution is changed from 6% to 16%, the curves shown in Fig. 11 almost remain the same, confirming that the parameters and are insensitive to the value of .
In summary, the experimental and MC simulated results demonstrate that the polarization Mueller matrix microscope can provide quantitative information of the pathological changes of liver tissues in different fibrosis stages. During the development of liver hepatitis and cancers, the inflammatory reactions can result in fibrosis and cirrhosis, which can be reflected by the values of the MMPD and MMT parameters. Although both the MMPD and MMT methods can be used for the detection of fibrous microstructures, they have different advantages: the MMPD parameters are slightly more sensitive to the fibers, whereas the MMT parameters are easier and quicker to compute.37 In this study, when a Mueller matrix is known, the calculation time of MMT parameters is 0.02 s, while the calculation time of MMPD parameters is 60.22 s with Intel(R) Core(TM) i7-3770 CPU and MATLAB 2014. It is shown that the Mueller matrix microscope and transformed parameters can be applied as potential tools to liver cancer stage scoring.
In this paper, we apply the Mueller matrix microscope to pathological liver cirrhosis and cancer tissue samples in the F1 to F4 fibrosis stages. We also adopt both the MMPD and MMT parameters for the quantitative analysis of the liver tissues. The microscopic imaging results demonstrate that the parameters , , , and can be used as tools for the detection of fibrous liver tissues. Since the pathological liver tissues in different stages have different proportions of fibrous microstructures, the experimental results show that the MMPD and MMT parameters can provide additional quantitative information helpful for the fibrosis detection and accurate scorings of liver cirrhosis and cancers. We conduct MC simulations based on the sphere birefringence model to analyze the relationship between the parameters and the fibrosis process of the liver tissues. Both the experimental and MC simulated results indicate that the polarized light microscope and transformed Mueller matrix parameters have a good application prospect in liver cancer diagnosis. In addition, through some more comprehensive transformation process of the Mueller matrix elements, we may obtain more quantitative parameters for the extraction of the intrinsic microstructural characteristic features of liver cancerous tissues. In future studies, statistical analysis of more fibrotic liver tissue samples is still needed to establish reliable relationship between the pathological structural features and the polarization imaging parameters. Also, analyzing methods such as the central moment parameters can be adopted to extract more quantitative indicators to facilitate fibrosis detections and scorings of liver cirrhosis and cancer tissues.
This work has been supported by National Natural Science Foundation of China (NSFC) (Grant Nos. 11174178, 11374179, 61205199, and 61405102) and Science and Technology Project of Shenzhen (Grant Nos. CXZZ20140509172959978 and GJHZ20150316160614844).