Blog

Quantification of collagen fiber properties in alcoholic liver fibrosis using polarization-resolved second harmonic generation microscopy | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  13, Article number: 22100 (2023 ) Cite this article Mineral Wool Rockwool Difference

Quantification of collagen fiber properties in alcoholic liver fibrosis using polarization-resolved second harmonic generation microscopy | Scientific Reports

Liver fibrosis is assessed mainly by conventional staining or second harmonic generation (SHG) microscopy, which can only provide collagen content in fibrotic area. We propose to use polarization-resolved SHG (PR-SHG) microscopy to quantify liver fibrosis in terms of collagen fiber orientation and crystallization. Liver samples obtained from autopsy cases with fibrosis stage of F0–F4 were evaluated with an SHG microscope, and 12 consecutive PR-SHG images were acquired while changing the polarization azimuth angle of the irradiated laser from 0° to 165° in 15° increments using polarizer. The fiber orientation angle (φ) and degree (ρ) of collagen were estimated from the images. The SHG-positive area increased as the fibrosis stage progressed, which was well consistent with Sirius Red staining. The value of φ was random regardless of fibrosis stage. The mean value of ρ (ρ-mean), which represents collagen fiber crystallinity, varied more as fibrosis progressed to stage F3, and converged to a significantly higher value in F4 than in other stages. Spatial dispersion of ρ (ρ-entropy) also showed increased variation in the stage F3 and decreased variation in the stage F4. It was shown that PR-SHG could provide new information on the properties of collagen fibers in human liver fibrosis.

Radu Hristu, Stefan G. Stanciu, ... George A. Stanciu

Danielle Tokarz, Richard Cisek, … Virginijus Beard

Adib Keikhosravi, Michael Shribak, … Kevin W. Eliceiri

When external stress or damage is applied, local inflammation occurs in tissues. Damaged tissue is removed by the inflammatory process, and fibrous tissue, mainly extracellular matrix such as collagen fibers, is formed to fill the space to maintain tissue shape and homeostasis. The excessive production and accumulation of collagen fibers beyond the physiological range is called fibrosis. In the liver, fibrosis is caused by chronic inflammation due to alcohol, viral hepatitis, and autoimmune hepatitis. In recent years, chronic inflammation and fibrosis of the liver caused by excessive fat accumulation in hepatocytes, especially in well-nourished industrialized countries, has become a social and medical problem, and this condition has attracted attention as nonalcoholic fatty liver disease (NAFLD) or nonalcoholic steatohepatitis (NASH)1. In the liver, the clinical outcomes of the patients such as carcinogenesis and liver failure have been reported to correlate with fibrosis stage and quantified collagen content2. Therefore, many basic and clinical studies have been performed to elucidate the mechanism of lever fibrosis. Previously, fibrosis was thought to only progress irreversibly. Following recent advances in antiviral therapy, it has become clear that there is a bidirectional outcome of progression and resolution, depending on the balance between collagen fiber production and degradation, and that reducing the degree of fibrosis reduces the risk of carcinogenesis3. However, the risk of carcinogenesis remains even without fibrosis4, suggesting that there are residual risks other than the amount of fibrosis, such as the characteristics of fibrosis.

Liver fibrosis has traditionally been assessed by histological methods such as Azan staining and Sirius Red staining. Liver fibrosis is pathologically categorized into five stages from F0 to F4 according to the new Inuyama classification or METAVIR score. The evaluation is not classified by a precisely quantified serial coefficient, but by a subjective staircase score based on the approximate percentage of fibrotic areas in the pathological specimen. For more objective evaluation, some attempts have been made to quantify the amount of fibrosis on images. However, it is still difficult to automatically distinguish between areas of physiological and pathological fibrosis, and this technique has not yet been put to practical use. It has also been noted that the nature of fibrosis varies according to the causative disease. The starting point and the pattern of fibrosis differ depending on the primary disease, such as viral hepatitis and alcoholic liver disease (ALD) or NASH. To understand the pathophysiology of fibrosis progression and regression, to accurately evaluate the therapeutic efficacy of antifibrotic agents, and to clarify the relationship between fibrosis and carcinogenesis for accurate assessment of the risk of carcinogenesis, it is essential not only to quantitatively assess the amount and area of fibrosis but also to accurately characterize the collagen fibers involved in fibrosis.

Second harmonic generation (SHG) microscopy has attracted attention as a new label-free method for evaluating fibrous tissue since early 2000s. SHG is a nonlinear optical phenomenon that generates light with a frequency twice that of the irradiated light (i.e., half the wavelength) when a hyperpolarizable non-centrosymmetric material is irradiated with intense pulsed light. In biological samples, collagen generates strong SHG light and is used for label-free observation of fibrotic tissue in skin, tendons, ligaments, and various other tissues5,6,7. Gailhouste et al. showed that SHG is useful in the assessment and scoring of liver fibrosis in hepatitis8. It has also been shown that a combination of SHG and two-photon excitation autofluorescence can be useful in assessing and scoring liver fibrosis due to hepatitis and NAFLD9. Since SHG per se can only assess the amount of fibrosis, more recently the concept of polarization-resolved SHG (PR-SHG) has been proposed. PR-SHG is based on the principle that pulsed laser irradiation with linear polarization parallel to the fiber orientation produces strong SHG light, whereas laser irradiation with orthogonal linear polarization produces weak SHG light10. PR-SHG can obtain new information on collagen fiber orientation as well as collagen content11,12,13. Although the previous studies have proved the efficacy of PR-SHG in analysis of collagen fibers in vitro and in animal models of liver fibrosis, the usefulness of PR-SHG in human liver tissue has not yet been verified. In the present study, we aimed to quantitatively evaluate the properties of collagen fibers in patients with ALD in terms of the orientation and crystallization using PR-SHG microscopy, and verified its potential for future clinical application.

We first compared the appearances of fibrotic regions indicated by Sirius Red staining and by wide-area SHG imaging. Since this study focused on ALD, the comparison focused on fibrosis around the central vein, which is a characteristic pathological finding of ALD. Figure 1 shows representative images of SHG microscopy and Sirius Red staining of the same region around a central vein in consecutive serial sections of the sample. The SHG-positive areas were highly consistent with positive areas of Sirius Red staining. As shown in Fig. 2a, the percentage of SHG-positive pixels to total pixels in the wide-area SHG image increased as the fibrosis stage progressed, with significantly large SHG-positive areas in F4 stage than in other stages. The total SHG light intensity also increased as the fibrosis stage progressed, and was markedly increased in stage F4 (Fig. 2b). The percentage of Sirius Red-stained area to the total area of the specimen increased with progression of fibrosis stage (Fig. 2c), showing the same trend as for SHG-positive area (Fig. 2a). There was a good correlation of the percentage of SHG-positive area with the percentage of Sirius Red-stained area in the same sample (correlation coefficient [R2] = 0.7763) (Fig. 2d).

Representative SHG and Sirius Red-stained images in each stage of fibrosis. Of two consecutive liver sections, one was imaged with SHG microscopy (left) and the other with Sirius Red staining (right). Fibrosis was compared in the same region around a central vein. Representative photographs are shown for each stage of fibrosis from F0 to F4 (corresponding to panels (a) through (e)). The scale bars represent 200 µm.

Summarized data and correlation of SHG and Sirius Red-stained images. (a) Percentage of SHG-positive area at each fibrosis stage. (b) Total SHG intensity within the observation area at each fibrosis stage. (c) Percentage of Sirius Red-stained area at each fibrosis stage. (d) Scatter plot of SHG-positive area vs. Sirius Red-stained area. †p < 0.01 vs. F0, F1, F2, and F3. ††p < 0.05 vs. F2 and F3. **p < 0.05 vs. F0, F1, and F2. One-way ANOVA with Tukey’s post-hoc test.

Collagen fiber orientation was evaluated by creating color maps and histograms of orientation angle (φ) and degree (ρ) (H(φ) and H(ρ), respectively) as shown in Fig. 3. The value of φ, which indicates the orientation angle of collagen fibers in the observed area, showed wide variation from sample to sample and image to image even within the same fibrosis stage, and we could not find any apparent tendency according to fibrosis stages. On the other hand, the value of ρ, which represents the crystallization of collagen fibers in the observed area, showed that the peak shifted to higher values as the fibrosis progressed.

Representative PR-SHG images and histograms of φ and ρ. Representative PR-SHG images of φ (upper left) and ρ (upper right) and histograms for φ (bottom left) and ρ (bottom right) based on each image obtained from samples of each stage of fibrosis from F0 to F4 (corresponding to panels (a) to (e)). The scale bars represent 10 µm.

Figure 4 summarizes the φ and ρ data. The orientation index of φ (φ-OI) as a measure of variation in collagen fiber orientation revealed no statistical difference by fibrosis stage (Fig. 4a). The mean value of ρ (ρ-mean), which represents the degree of collagen fiber crystallinity, showed a very characteristic change. In stage F0, ρ-mean was approximately 2.1 and the inter-quartile range (IQR) was 0.189, showing little variation. As fibrosis stage progressed from F1 to F2, the average and median values of ρ-mean did not change significantly, whereas the IQR increased slightly to around 0.20. At stage F3, there was little change in the average and median values but the IQR was 0.604, indicating large variation in the data. At stage F4, ρ-mean was 2.5, which was significantly higher than at any other stage, and the IQR was 0.097, showing less variation (Fig. 4b). A similar trend was observed for ρ-entropy, which represents spatial heterogeneity of collagen fiber crystallinity in the observed area. The IQR of ρ-entropy was higher at stage F3 (0.264) than stages F0–F2, although there was no change in the mean or median values. At stage F4, the average value of ρ-entropy (− 3.12) was significantly lower and the IQR (0.075) showed less variation compared with other stages. We further evaluated the relationship between ρ-mean and ρ-entropy of each of the two samples. When ρ-entropy was plotted against ρ-mean, all data were linearly distributed from the upper left (i.e., the region of low ρ-mean and high ρ-entropy) to the lower right (i.e., the region of high ρ-mean and low ρ-entropy) (Fig. 5). Among them, data from F0 samples were concentrated in the upper left of the scatter plot (solid oval), whereas data from F4 samples were concentrated in the lower right of the scatter plot (dotted oval). The data for the F1 to F3 stages were widely distributed among these data.

Summarized data of φ-OI, ρ-mean, and ρ-entropy. The data of φ-OI (a), ρ-mean (b), and ρ-entropy (c) compiled into box-and-whisker diagrams are superimposed on a plot of actual data from F0 to F4. The boxes indicate the inter-quartile range (IQR) and the bars above and below the boxes indicate the range between maximum and minimum values. †p < 0.01 vs. F0, F1, F2, and F3. *p < 0.01 vs. F3. One-way ANOVA with Tukey’s post-hoc test.

Scatterplot showing the distribution of ρ-entropy relative to ρ-mean. White diamonds indicate F0, gray circles indicate F1–F3 (light, medium, and dark grey, respectively), and black diamonds indicate F4 data. The solid and dotted ovals indicate areas where the scatter data accumulates in stages F0 and F4, respectively.

In the present study, we assessed fibrosis due to alcoholic liver injury in human postmortem tissue using SHG microscopy. The major findings of the study were as follows. (1) SHG microscopy could accurately detect liver fibrosis, and there was good correlation of SHG images with fibrosis stage based on Sirius Red staining. (2) We have confirmed for the first time that PR-SHG microscopy could estimate the orientation angle (φ) and degree (ρ) of collagen fibers in human liver tissue. (3) We further evaluated collagen fiber orientation in terms of φ-OI, ρ-mean, and ρ-entropy in human liver tissue. (4) φ-OI showed no obvious correlation with stage of fibrosis, whereas rho-mean and rho-entropy showed data dispersion in the process of liver tissue inflammation caused by alcohol and data convergence in cirrhosis.

In recent years, the development of new imaging techniques such as shear-wave elastography (SWE), magnetic resonance elastography (MRE), and transient elastography (FibroScan®) has been actively pursued for quantitative evaluation of liver fibrosis in a noninvasive manner14, which is becoming increasingly important as the number of patients with NAFLD continues to increase rapidly. Following the development and use of anti-fibrotic agents in clinical medicine15,16, it is also becoming increasingly necessary to distinguish between drug-responsive and drug-resistant fibrosis to clearly define the treatment criteria. SWE and MRE enable noninvasive assessment of fat accumulation and fibrosis in the liver, and their diagnostic value has been increasing in recent years due to enhanced spatial resolution. The FibroScan® is commonly used for noninvasive monitoring of the progression and regression of liver fibrosis17. SWE has also been shown to accurately reflect changes in liver stiffness induced by antiviral treatment in viral hepatitis18. MRE can analyze the elastic and viscous components of the liver independently by applying external vibration with different frequencies, which is helpful in distinguishing fibrosis, inflammation, and congestion19. These new methods have the advantage of noninvasive diagnosis, however, confounding factors and spatial resolution still remain problems20,21. Electron microscopy is a method to observe collagen fibers in more detail. However, electron microscopy is time-consuming to prepare samples, and it is questionable whether it accurately reflects the clinical phenotypes because of the artifacts during sample preparation. On the other hand, SHG microscopy can be used to observe organs immediately after removal as long as a thin slice specimen is prepared, and organs can be observed as they are if a microscope with a reflective configuration is prepared. The assessment of liver fibrosis by PR-SHG would enable the evaluation of fibrosis specific to collagen fibers at the microscopic level and have potential to bridge the gap between macro (such as liver stiffness and serum biomarkers) and micro information (such as collagen fiber properties).

SHG imaging has recently attracted attention as a novel label-free evaluation of fibrotic tissues and collagen fibers. It has been used to visualize and evaluate collagen in soft tissues such as skin, tendons, and ligaments, as well as pathological fibrosis in cancer and liver injury5,6,7,8,9,22. We found high consistency between SHG imaging and Sirius Red staining according to fibrosis stage, which reaffirms the potential of SHG as a method for assessing liver fibrosis. However, the information obtained by SHG microscopy is limited to the amount and distribution of collagen and not up to the superiority over conventional staining evaluation. PR-SHG microscopy additionally provided us the indices of collagen fiber orientation angle (φ) and degree (ρ), which are not available by conventional microscopy or SHG microscopy with normal configuration. Among them, we could not find any statistically significant difference in φ corresponding to fibrosis stages in human liver tissue. As φ represents fiber orientation, it could be reasonable to use PR-SHG microscopy to assess pathological changes in tissues where fibers are oriented uniformly in a certain direction, for example in tendons and ligaments. On the other hand, in fibrosis of parenchymal organs such as the liver, fibrosis occurs to fill in tissue damaged by inflammation and other factors, so the orientation of collagen fibers is not aligned in a certain direction. This might be the reason why little difference in φ depending on fibrosis stage was observed in liver fibrosis in the present study. In contrast, the value of ρ, which represents the crystallinity of collagen fibers, showed characteristic changes according to the stage of fibrosis. The data in F4 stage had higher ρ-mean and lower ρ-entropy than the other stages, suggesting that collagen in F4 fibrosis had heterogenous orientation at fibril level, and also had higher spatial heterogeneity in the organ. In self-assembling collagen gel and in a rat model of liver fibrosis induced by bile duct ligation, the anisotropy of SHG and the ratio of forward to backward scattering of PR-SHG were shown to change, possibly due to an increase in type III collagen relative to type I collagen23,24,25. The previous article also suggested SHG images in breast cancer differed depending on cross-linking between collagen fibers26. The morphological organization of collagen within breast cancer samples depends not only on cross-linking but also on its supramolecular symmetry, that can be assessed using PR-SHG microscopy27. The changes in ρ-mean and ρ-entropy that we observe in the study might represent the changes in collagen type or in the cross-linking status between collagen fibers. Thus, it was significant to confirm that such collagen fiber-related data could be obtained from human liver tissue by PR-SHG for the first time, even though from postmortem samples. In the future, the ratio of type I to type III collagen or the cross-links between collagen fibers could be correlated with the value of ρ-mean and ρ-entropy obtained by PR-SHG to provide a more accurate information of the nature of fibrosis.

Liver fibrosis was previously thought to be a unidirectional process that worsened with inflammation and accumulation of tissue damage and was generally staged by the new Inuyama classification or METAVIR score, which determines the fibrosis stage based on the amount of visually observed fibrotic tissue. Recent basic research has shown that fibrosis in parenchymal organs such as liver could have both progressive and regressive outcomes, depending on the balance of collagen fiber formation and degradation. Furthermore, because the extent and timing of inflammatory cell infiltration into the tissue is not uniform, inflammation in the tissue is spatiotemporally heterogeneous, which would be expected to result in spatiotemporal heterogeneity in the formation and degradation of fibrosis. In contrast, in cases where fibrosis has progressed significantly to cirrhosis, the processes of collagen fiber formation and degradation are inactivated, and fibrosis is considered to be irreversible and to lose its spatial and temporal heterogeneity. Our experimental data could support this pathophysiology of the fibrotic process in the liver. In the present data, the variation in ρ-mean and ρ-entropy values increased as the fibrosis stage progressed from F1 to F3, and the F4 samples showed less variation in ρ-mean and ρ-entropy. Especially in the case of F3 samples, it was assumed that there was a large variation in the data due to temporal heterogeneity in the progression of fibrosis as well as sample-to-sample variation. In contrast, stage F4 was dominated by complete and stable fibrosis, resulting in less variation in the data. This speculation was also supported by the finding that total SHG intensity was significantly higher in the stage F4 samples than in the others, as it has been suggested that SHG luminosity is higher in thicker, mature collagen fibers28,29. Thus in liver fibrosis that repeatedly progresses and regresses, the PR-SHG might be useful in precisely identifying the pathological state of the patients, and in the future, it is expected to be combined with morphological classifications such as the new Inuyama classification and METAVIR score to enable more detailed classification of liver fibrosis. If such knowledge can be accumulated, it might enable us to predict whether fibrotic regions have a tendency to form and expand or to break down and shrink, whether they still have plasticity or already in irreversible state, and even the responses to anti-fibrotic therapy.

The present study has several limitations. First, it should be noted that the study was validated only in cases of ALD. Pathological conditions that produce liver fibrosis include viral hepatitis and NAFLD or NASH, which has been viewed as a major medical problem in recent years. It should also be noted that this study just focused on cases of ALD used postmortem samples and evaluated only one point for each specimen. However, all the autopsy cases analyzed in this study were untreated patients with respect to alcoholic liver disease, which enabled evaluation of samples in a purely fibrotic state according to stage without any intervention. In addition, since ALD is mainly characterized by fibrosis around the central vein, it was relatively easy to interpret the PR-SHG data by focusing on the area around the central vein for comparison. Nevertheless, it was still difficult to interpret data with large individual differences and spatial heterogeneity among samples and temporal heterogeneity due to the time phase of the pathology, and we have not been able to fully examine how this heterogeneous data reflected the clinical pathology. To more accurately understand the meaning of ρ-mean and ρ-entropy in the pathogenesis of liver fibrosis, it will be necessary to conduct animal experiments in which liver fibrosis can be observed over time in single individuals, and clinical studies in which patients treated with antifibrotic therapy are followed over time. By integrating the various information obtained from these studies, such as biomarkers, imaging findings such as ultrasound and CT, fibrosis classification based on conventional staining methods such as the new Inuyama classification and METAVIR score, quantitative information on the amount of fibrosis by conventional SHG imaging, and quantitative information on the nature of fibrosis by PR-SHG, it would be possible to classify fibrosis in a new and more detailed way that is more relevant to clinical symptoms and treatment indications.

Among forensic autopsies performed at the Office of the Tokyo Metropolitan Coroner between 2016 and 2019, we analyzed liver samples obtained from patients whose medical histories confirmed their constant alcohol consumption before death (Table 1). The degree of fibrosis in each case was assessed by Dr. Kimura, the forensic scientist in charge, using formalin-fixed paraffin-embedded (FFPE) specimens of the liver, according to the new Inuyama classification. No personal information was linked to the specimen number, and use of the specimens was limited to histological evaluation. Ethical review was performed and approved by the Tokyo Metropolitan Medical Examiner's Office for these uses (No. 2019-2). All methods were carried out in accordance with relevant guidelines and regulations. This study used only already existing samples for which autopsies have been performed in the past, and all samples was not identifiable to any specific individual. Therefore, we did not obtain informed consent, and the information about the research was made available to all bereaved families by the postings and the website of Tokyo Medical Examiner’s Office (https://www.hokeniryo.metro.tokyo.lg.jp/kansatsu/kenkyuu/index.html), and the opportunity for the families to refuse the research being conducted was ensured.

Of two consecutive FFPE sections with 5-µm thickness obtained per case, one was stained with Sirius Red for conventional microscopy and the other was used for SHG microscopy without any staining. For Sirius Red staining, specimens were immersed in a 3:100 mixture of 1% Sirius Red (Muto Pure Chemical, Tokyo, Japan) and saturated picric acid (Wangysong Solution A, Muto Pure Chemical) for 10 min after deparaffinization. The specimens were then dehydrated with 100% ethanol and xylene and coverslipped.

For bright field microscopic analysis of Sirius Red-stained specimens, we used a BZ-X700 automated microscopy system (Keyence Corporation, Osaka, Japan) controlled by a software BZ-X Viewer. Briefly, an image of normal resolution (960 × 720 pixels) was acquired through a 4 × objective lens (CFI Plan Apo Lambda 4X, Nikon, Inc., Tokyo, Japan; magnification = 4; numerical aperture (NA) = 0.20; working distance = 20 mm), and multiple images were combined by tiling to create a macroscopic image of the entire specimen. The images were then stored in a tagged image file format and analyzed using BZ-X Analyzer software ver.1.3.0.3. The hybrid cell counting program was applied using bright field mode and region masking with one-color extraction based on RGB color mode. The analyzed area was set around the periphery of the entire sample within a polygonal region of interest (ROI), and color extraction was performed by setting the cursor on the Sirius Red-stained area inside the ROI and adjusting it so that the entire Sirius Red-stained area was masked within a tolerance of 35–45. The automated measurement outputs the area of the entire sample and the area of the Sirius-Red-stained area, and the percentage of the Sirius Red-stained area to the sample area was calculated from these values.

To perform SHG imaging and PR-SHG imaging, we used a custom-built laser scanning SHG microscope based on a wavelength-tunable femtosecond optical parametric oscillator (InSight DeepSee, Spectra-Physics, Inc., Tokyo, Japan; tuning range = 680–1300 nm; pulse duration = ~ 110 fs; repetition rate = 80 MHz) and an inverted optical microscope (Ti2-U, Nikon, Inc.). The central wavelength and the laser power at the sample were set to be 800 nm and 10–20 mW, respectively. The focal spot was scanned two-dimensionally onto the sample using a pair of galvanometer mirrors (GMs), a pair of relay lenses, and an objective lens (CFI Plan Apo Lambda 60XC, Nikon, Inc.; magnification = 60; NA = 0.95; working distance = 110–210 μm). Forward-propagated SHG light was collected via the condenser lens (TI-C-LWD, Nikon, Inc.; NA = 0.52) and then separated from the excitation laser light by a dichroic mirror (FF705-Di01, Semrock, Inc., Rochester, NY, the US) and an optical band-pass filter (FBH405-10, Thorlabs, Inc., Newton, NJ, the US). Finally, the SHG light was detected by a photon-counting photomultiplier (H8259-01, Hamamatsu Photonics, K.K., Hamamatsu, Japan) connected to a pulse counter. Using the above setup, SHG images of a 137 × 137 μm region, composed of 128 × 128 pixels, were acquired at a rate of 1 image/s. To further expand the imaging region, we scanned the sample position horizontally or vertically at intervals of 137 μm using a stepping-motor-driven translation stage every time an SHG image was acquired. Finally, a large-area SHG image of size 1.37 × 1.37 mm was obtained, corresponding to 1280 × 1280 pixels, by stitching 100 SHG images together in a matrix of 10 rows and 10 columns.

Groups of PR-SHG images were acquired using an additional polarization control unit based on a half-waveplate (WPH05M-808, Thorlabs, Inc.) and quarter-waveplate (WPQ05M-808, Thorlabs, Inc.) installed before the GMs. The input linear polarization angles were rotated uniformly from 0° to 165° in 15° steps. All linear polarization states had a good polarization extinction ratio (> 100) at the back pupil of the objective lens.

PR-SHG analysis exploits the difference in SHG light intensity in SHG images acquired while rotating the input linear polarization angle. Following the method used in previous studies12,30,31,32,33,34, we retrieved the fiber orientation angle φ and SHG anisotropic parameter ρ in specimens at single-pixel resolution. Briefly, when linearly polarized laser light is incident at angle α on collagen fiber, the resultant SHG intensity I2ω measured without an analyzer in front of the detector can be written as

where K is the proportional constant depending on an experimental condition such as laser power, and φ denotes the angle of the main axis of the fiber with respect to the x (horizontal)-axis on the SHG image. α is also defined as the angle between the incident laser polarization and the x-axis on the SHG image. Here, we assumed that the collagen fiber lay in the xy plane and that the laser light propagated in the z-direction in Cartesian coordinates xyz for the local frame. In this formalism, we also assumed that the collagen fiber possessed cylindrical symmetry (C∞). In addition, Kleinman symmetry can be applied because of the difference between the photon energy of SHG light and the first electronic transition in the collagen molecule10. Therefore, the second-order nonlinear optical susceptibility tensor χ(2) has only two independent nonvanishing tensorial components: χ(2)xxx and χ(2)xyy. ρ is the ratio of these remaining components (ρ = χ(2)xxx/χ(2)xyy), and we considered that this parameter was inversely proportional to the degree of collagen orientation and/or the structural maturity of collagen fiber. To obtain an analytic solution of these values, we applied linear least square (LLS) fitting34 to the polarization-dependent SHG images at single-pixel resolution and reconstructed φ and ρ images of the same pixel size as the SHG image. Due to the fast-fitting algorism in the LLS method, the total image processing time for an image of 128 × 128 pixels was within 1 min.

To obtain quantitative values from the φ and ρ images, we created histograms of the φ and ρ values in each image and calculated three parameters, φ-OI, ρ-mean, and ρ-entropy, as previously described35,36,37.

From the φ image, we calculated the orientation index (OI), which reflects the percentage of fibers with orientation angle parallel to the main orientation angle (φmax) in the φ-histogram (H(φ)), defined as

We then calculated the mean ρ value of the ρ images (ρmean) and the statistical entropy (S), which was defined as

is the probability of encountering a pixel that has the value ρ. Statistical entropy expresses the randomness of the ρ value in the ρ-histogram (H(ρ)), with lower S values corresponding to a more uniform ρ-distribution. By definition of the calculation, a lower value of ρ-mean indicates a more crystallized fiber, and a lower value of ρ-entropy indicates a higher spatial heterogeneity of fiber crystallinity.

The data are expressed as mean ± standard error of the mean. Differences between group averages were tested by a one-way analysis of variance with a post-hoc test (Tukey’s honestly significant difference) and values of p < 0.05 were considered significant. Statistical analysis was performed using R software ver.3.6.338.

Data obtained from the experiments are available from the corresponding author upon reasonable request.

Le, M. H. et al. Global incidence of non-alcoholic fatty liver disease: A systematic review and meta-analysis of 63 studies and 1,201,807 persons. J. Hepatol. https://doi.org/10.1016/j.jhep.2023.03.040 (2023).

Calvaruso, V. et al. Computer-assisted image analysis of liver collagen: Relationship to Ishak scoring and hepatic venous pressure gradient. Hepatology 49, 1236–1244. https://doi.org/10.1002/hep.22745 (2009).

Wanless, I. R., Nakashima, E. & Sherman, M. Regression of human cirrhosis. Morphologic features and the genesis of incomplete septal cirrhosis. Arch. Pathol. Lab. Med. 124, 1599–1607. https://doi.org/10.5858/2000-124-1599-rohc (2000).

Article  CAS  PubMed  Google Scholar 

Calvaruso, V. et al. Liver collagen proportionate area predicts decompensation in patients with recurrent hepatitis C virus cirrhosis after liver transplantation. J. Gastroenterol. Hepatol. 27, 1227–1232. https://doi.org/10.1111/j.1440-1746.2012.07136.x (2012).

Medyukhina, A. et al. Automated classification of healthy and keloidal collagen patterns based on processing of SHG images of human skin. J. Biophotonics 4, 627–636. https://doi.org/10.1002/jbio.201100028 (2011).

Article  CAS  PubMed  Google Scholar 

Hase, E. et al. Evaluation of the histological and mechanical features of tendon healing in a rabbit model with the use of second-harmonic-generation imaging and tensile testing. Bone Joint Res. 5, 577–585. https://doi.org/10.1302/2046-3758.511.Bjr-2016-0162.R1 (2016).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Lee, W., Rahman, M. H., Kersh, M. & Toussaint, K. Application of quantitative second-harmonic generation microscopy to posterior cruciate ligament for crimp analysis studies. J. Biomed. Opt. 22, 046009 (2017).

Gailhouste, L. et al. Fibrillar collagen scoring by second harmonic microscopy: A new tool in the assessment of liver fibrosis. J. Hepatol. 52, 398–406. https://doi.org/10.1016/j.jhep.2009.12.009 (2010).

Article  CAS  PubMed  Google Scholar 

Chang, P. E. et al. Second harmonic generation microscopy provides accurate automated staging of liver fibrosis in patients with non-alcoholic fatty liver disease. PLoS One 13, e0199166. https://doi.org/10.1371/journal.pone.0199166 (2018).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Yew, E. Y. S. & Sheppard, C. J. R. Effects of axial field components on second harmonic generation microscopy. Opt. Express 14, 1167–1174. https://doi.org/10.1364/OE.14.001167 (2006).

Article  ADS  PubMed  Google Scholar 

Chen, X., Nadiarynkh, O., Plotnikov, S. & Campagnola, P. J. Second harmonic generation microscopy for quantitative analysis of collagen fibrillar structure. Nat. Protoc. 7, 654–669. https://doi.org/10.1038/nprot.2012.009 (2012).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Rouède, D. et al. Determination of extracellular matrix collagen fibril architectures and pathological remodeling by polarization dependent second harmonic microscopy. Sci. Rep. 7, 12197. https://doi.org/10.1038/s41598-017-12398-0 (2017).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Mercatelli, R. et al. Morpho-mechanics of human collagen superstructures revealed by all-optical correlative micro-spectroscopies. Commun. Biol. 26(2), 117. https://doi.org/10.1038/s42003-019-0357-y (2019).

Tamaki, N., Kurosaki, M., Huang, D. Q. & Loomba, R. Noninvasive assessment of liver fibrosis and its clinical significance in nonalcoholic fatty liver disease. Hepatol. Res. 52, 497–507. https://doi.org/10.1111/hepr.13764 (2022).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Caligiuri, A., Gentilini, A., Pastore, M., Gitto, S. & Marra, F. Cellular and molecular mechanisms underlying liver fibrosis regression. Cells 10, 2759 (2021).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Tan, Z. et al. Liver fibrosis: Therapeutic targets and advances in drug therapy. Front. Cell Dev. Biol. https://doi.org/10.3389/fcell.2021.730176 (2021).

Article  PubMed  PubMed Central  Google Scholar 

Castera, L., Forns, X. & Alberti, A. Non-invasive evaluation of liver fibrosis using transient elastography. J. Hepatol. 48, 835–847. https://doi.org/10.1016/j.jhep.2008.02.008 (2008).

Tada, T. et al. Improvement of liver stiffness in patients with hepatitis C virus infection who received direct-acting antiviral therapy and achieved sustained virological response. J. Gastroenterol. Hepatol. 32, 1982–1988. https://doi.org/10.1111/jgh.13788 (2017).

Article  CAS  PubMed  Google Scholar 

Yin, M. et al. Distinguishing between hepatic inflammation and fibrosis with MR elastography. Radiology 284, 694–705. https://doi.org/10.1148/radiol.2017160622 (2017).

Kennedy, P. et al. Quantitative elastography methods in liver disease: Current evidence and future directions. Radiology 286, 738–763. https://doi.org/10.1148/radiol.2018170601 (2018).

Castéra, L. et al. Pitfalls of liver stiffness measurement: A 5-year prospective study of 13,369 examinations. Hepatology 51, 828–835. https://doi.org/10.1002/hep.23425 (2010).

Keikhosravi, A., Bredfeldt, J. S., Sagar, A. K. & Eliceiri, K. W. In Methods in Cell Biology Vol. 123 (eds Waters, J. C. & Wittman, T.) 531–546 (Academic Press, 2014).

Karsdal, M. A. et al. Fibrogenesis assessed by serological type III collagen formation identifies patients with progressive liver fibrosis and responders to a potential antifibrotic therapy. Am. J. Physiol. Gastrointest. Liver Physiol. 311, G1009–G1017. https://doi.org/10.1152/ajpgi.00283.2016 (2016).

Tilbury, K., Lien, C.-H., Chen, S.-J. & Campagnola, P. J. Differentiation of Col I and Col III isoforms in stromal models of ovarian cancer by analysis of second harmonic generation polarization and emission directionality. Biophys. J. 106, 354–365. https://doi.org/10.1016/j.bpj.2013.10.044 (2014).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Lin, J., Pan, S., Zheng, W. & Huang, Z. Polarization-resolved second-harmonic generation imaging for liver fibrosis assessment without labeling. Appl. Phys. Lett. 103, 173701. https://doi.org/10.1063/1.4826516 (2013).

Article  ADS  CAS  Google Scholar 

Levental, K. R. et al. Matrix crosslinking forces tumor progression by enhancing integrin signaling. Cell 139, 891–906. https://doi.org/10.1016/j.cell.2009.10.027 (2009).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Mercatelli, R., Triulzi, T., Pavone, F. S., Orlandi, R. & Cicchi, R. Collagen ultrastructural symmetry and its malignant alterations in human breast cancer revealed by polarization-resolved second-harmonic generation microscopy. J. Biophotonics 13(8), e202000159. https://doi.org/10.1002/jbio.202000159 (2020).

Article  CAS  PubMed  Google Scholar 

Campagnola, P. J. & Dong, C.-Y. Second harmonic generation microscopy: Principles and applications to disease diagnosis. Laser Photonics Rev. 5, 13–26. https://doi.org/10.1002/lpor.200910024 (2011).

Article  ADS  CAS  Google Scholar 

Hase, E., Matsubara, O., Minamikawa, T., Sato, K. & Yasui, T. In situ time-series monitoring of collagen fibers produced by standing-cultured osteoblasts using a second-harmonic-generation microscope. Appl. Opt. 55, 3261–3267. https://doi.org/10.1364/AO.55.003261 (2016).

Article  ADS  CAS  PubMed  Google Scholar 

Stoller, P., Reiser, K. M., Celliers, P. M. & Rubenchik, A. M. Polarization-modulated second harmonic generation in collagen. Biophys. J. 82, 3330–3342. https://doi.org/10.1016/S0006-3495(02)75673-7 (2002).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Gusachenko, I., Tran, V., Houssen, Y. G., Allain, J.-M. & Schanne-Klein, M.-C. Polarization-resolved second-harmonic generation in tendon upon mechanical stretching. Biophys. J. 102, 2220–2229. https://doi.org/10.1016/j.bpj.2012.03.068 (2012).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Latour, G., Gusachenko, I., Kowalczuk, L., Lamarre, I. & Schanne-Klein, M.-C. In vivo structural imaging of the cornea by polarization-resolved second harmonic microscopy. Biomed. Opt. Express 3, 1–15. https://doi.org/10.1364/BOE.3.000001 (2012).

Gusachenko, I., Latour, G. & Schanne-Klein, M.-C. Polarization-resolved second harmonic microscopy in anisotropic thick tissues. Opt. Express 18, 19339–19352. https://doi.org/10.1364/OE.18.019339 (2010).

Article  ADS  CAS  PubMed  Google Scholar 

Rouède, D., Bellanger, J.-J., Bomo, J., Baffet, G. & Tiaho, F. Linear least square (LLS) method for pixel-resolution analysis of polarization dependent SHG images of collagen fibrils. Opt. Express 23, 13309–13319. https://doi.org/10.1364/OE.23.013309 (2015).

Article  ADS  CAS  PubMed  Google Scholar 

Kim, A., Lakshman, N. & Petroll, W. M. Quantitative assessment of local collagen matrix remodeling in 3-D culture: The role of Rho kinase. Exp. Cell Res. 312, 3683–3692. https://doi.org/10.1016/j.yexcr.2006.08.009 (2006).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Bayan, C., Levitt, J. M., Miller, E., Kaplan, D. & Georgakoudi, I. Fully automated, quantitative, noninvasive assessment of collagen fiber content and organization in thick collagen gels. J. Appl. Phys. 105, 102042. https://doi.org/10.1063/1.3116626 (2009).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Bancelin, S. et al. Ex vivo multiscale quantitation of skin biomechanics in wild-type and genetically-modified mice using multiphoton microscopy. Sci. Rep. 5, 17635. https://doi.org/10.1038/srep17635 (2015).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

R: A language and environment for statistical computing (R Foundation for Statistical Computing, 2020).

The study was financially supported by JSPS KAKENHI (Grant Nos. JP22K07024 [K.T.], JP21K19913 [H.T.], and JP21K18038 [E.H.]).

These authors contributed equally: Saya Matsuzaki, Eiji Hase and Hiroki Takanari.

Department of Radiology, Institute of Biomedical Sciences, Tokushima University Graduate School of Medicine, Tokushima, Japan

Saya Matsuzaki & Masafumi Harada

Division of Interdisciplinary Research for Medicine and Photonics, Institute of Post-LED Photonics, Tokushima University, Tokushima, Japan

Eiji Hase, Hiroki Takanari, Yuri Hayashi, Yusaku Hayashi, Haruto Oshikata, Takeo Minamikawa & Koichi Tsuneyama

Department of Legal Medicine, Nihon University School of Medicine, Tokyo, Japan

Tokushima University Faculty of Medicine, Tokushima, Japan

Yuri Hayashi, Yusaku Hayashi & Haruto Oshikata

Tokyo Medical Examiner’s Office, Tokyo, Japan

Department of Pathology and Laboratory Medicine, Institute of Biomedical Sciences, Tokushima University Graduate School of Medicine, 3-18-15, Kuramoto, Tokushima, 770-8503, Japan

Mayuko Ichimura-Shimizu & Koichi Tsuneyama

Division of Next-Generation Photonics, Institute of Post-LED Photonics, Tokushima University, Tokushima, Japan

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

S.M., E.H., H.T. and K.T. wrote the main manuscript text and figures. Y.H., Y.H. and H.O. helped to measure SHG samples. T.M. and T.Y. created the SHG device. S.K. collected, diagnosed, and provided measured autopsy specimens. M.I.-S., T.Y. and M.H. assisted in the preparation of the manuscript and provided various advice on additional experiments. All authors reviewed the manuscript.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Matsuzaki, S., Hase, E., Takanari, H. et al. Quantification of collagen fiber properties in alcoholic liver fibrosis using polarization-resolved second harmonic generation microscopy. Sci Rep 13, 22100 (2023). https://doi.org/10.1038/s41598-023-48887-8

DOI: https://doi.org/10.1038/s41598-023-48887-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Quantification of collagen fiber properties in alcoholic liver fibrosis using polarization-resolved second harmonic generation microscopy | Scientific Reports

Basalt Composite Fiber Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.