Digital Image Analysis of Picrosirius Red Staining: A Robust Method for Multi-Organ Fibrosis Quantification and Characterization

Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Associated Data

GUID: DD30F16B-8F1C-4086-8BF8-6BD83FC8B0D2

Abstract

Current understanding of fibrosis remains incomplete despite the increasing burden of related diseases. Preclinical models are used to dissect the pathogenesis and dynamics of fibrosis, and to evaluate anti-fibrotic therapies. These studies require objective and accurate measurements of fibrosis. Existing histological quantification methods are operator-dependent, organ-specific, and/or need advanced equipment. Therefore, we developed a robust, minimally operator-dependent, and tissue-transposable digital method for fibrosis quantification. The proposed method involves a novel algorithm for more specific and more sensitive detection of collagen fibers stained by picrosirius red (PSR), a computer-assisted segmentation of histological structures, and a new automated morphological classification of fibers according to their compactness. The new algorithm proved more accurate than classical filtering using principal color component (red-green-blue; RGB) for PSR detection. We applied this new method on established mouse models of liver, lung, and kidney fibrosis and demonstrated its validity by evidencing topological collagen accumulation in relevant histological compartments. Our data also showed an overall accumulation of compact fibers concomitant with worsening fibrosis and evidenced topological changes in fiber compactness proper to each model. In conclusion, we describe here a robust digital method for fibrosis analysis allowing accurate quantification, pattern recognition, and multi-organ comparisons useful to understand fibrosis dynamics.

Keywords: fibrosis, picrosirius red, digital analysis, collagen proportionate area, whole section, region-of-interest, fibrosis pattern

1. Introduction

Fibrosis is a pathological manifestation of excessive or inappropriate wound healing resulting in collagen deposition that may affect the function of diverse organs [1]. So far, no effective anti-fibrotic therapy is available while the burden of fibrotic diseases is increasing and intense research on preclinical models is underway. Therefore, accurate and reliable quantification of collagen deposits and their characterization are necessary to stage and map the natural history of diseases, to study the pathogenesis of fibrosis, to identify therapeutic targets, to evaluate anti-fibrotic strategies or to validate non-invasive diagnostic tools.

Collagen biochemical assays (e.g., hydroxyproline measurements) on total tissue homogenate provide a global estimation of collagen levels within a tissue specimen. However, this approach faces multiples limitations: (1) techniques are somehow painstaking (hydroxyproline content or sodium dodecyl sulfate-polyacrylamide gel electrophoresis followed by densitometry analysis) or require purification steps to limit overestimation (Sircol) [2], carrying the risk of missing a significant component in discarded fractions; (2) they include both induced and structural collagen; (3) histological compartments cannot be distinguished, complicating data interpretation and hiding subtle changes [3].

Alternatively, histopathological parameters based on collagen staining on tissue sections have been validated in many organs and disease conditions for the grading and the staging of fibrosis or for follow-up purposes. In this approach, the first critical step is to provide a specific signal. The most common techniques to visualize collagen deposits are the Masson’s trichrome and picrosirius red (PSR) (Sirius red F3B, Direct red 80). Both dyes possess good affinity for collagen fibers owe to multiple base-acid interactions, with PSR having a higher binding affinity for fibrillary collagen because of its elongated chemical structure, as confirmed by immunodetection [4,5,6]. PSR staining offers numerous advantages: besides being inexpensive, PSR solution is stable for years, the staining protocol is fast and simple, consistent in terms of sensitivity and hue [5], and provides a contrasted staining that enables a more reliable evaluation of fibrosis than with the Masson’s trichrome [7,8,9]. Several protocols have been described for PSR staining, the most common being with yellow counterstaining by picric acid only which provides a good signal-to-noise ratio [4,10], with the use of additional dyes for specific applications [11,12] or combined with immunohistochemistry but with some limitations [13].

The staining interpretation represents the second critical step. It requires (1) the accurate and reliable detection of the signal over noise, (2) the identification of the relevant histological structures, and (3) a meaningful quantification. This can be achieved by skilled pathologists or using digital image analysis (DIA), both methods having intrinsic strengths and limitations. On the one hand, trained human brains can perform incredibly fast for completing scoring systems (i.e., Ishak, Metavir, or Ashcroft systems), but these human resources are usually limited and inter- and even intra-individual variations are a well-known source of error [14,15]. Furthermore, these scoring systems are by essence descriptive or only semi-quantitative, and comprehensive information from whole tissue PSR-stained sections remains under-exploited. On the other hand, modern DIA has benefited of the onset of artificial intelligence. Advanced algorithms have been developed to recognize some specific fibrosis patterns [16,17,18], but these methods require software that can be expensive and proper training. Furthermore, they are not intended for multi-organ comparison since their transposition would require specific adaptation and optimization. This can be solved by computer-assisted recognition of tissue architecture followed by automated quantification of fibrosis on large tissue samples [19,20,21,22], thereby outclassing visual scoring on selected fields of view.

Outcomes generated by DIA can then be expressed as the area occupied by collagen out of the entire examined tissue section (referred to as collagen proportionate area, CPA). This accessible method is in use for fibrosis quantification in the liver [8,23,24,25,26,27], kidney [9,28,29], or lung [26,30]. Additionally, DIA combined with advanced microscopy (second harmonic generation) allows the morphometric segmentation of collagen fibers and hence discriminating aggregated and distributed fibers [31,32]. Nevertheless, to the best of our knowledge, no transposable method to quantify fibrosis in multiple organs and to characterize fibrosis dynamics on PSR-stained sections imaged in brightfield has been reported.

The present DIA method was developed to offer (1) a robust and versatile quantification that enables the regionalized analysis of fibrosis in diverse tissues for multi-organ or inter-study comparison and (2) objective parameters to characterize the pattern of collagen deposits by the comprehensive analysis on large tissue specimens.

2. Materials and Methods

2.1. Murine Models of Fibrosis

We used in the present study recognized models of liver, kidney or lung fibrosis. All male mice were kept under standard laboratory conditions in 12:12 h light-dark cycles, and food and water were supplied ad libitum. Animal procedures were in accordance with the National Research Council Guide for the Care and Use of Laboratory Animals, in agreement with ARRIVE guidelines 2.0, and were approved by the Institutional Review Board of UCLouvain (liver model: 2012/UCL/MD/026; lung model: 2018/UCL/MD/012; kidney model: 2018/UCL/MD/015). In every model, control animals were compared with animals exposed to the fibrosis agent and sacrificed at two time points corresponding to intermediate and severe fibrosis. We analyzed n = 5 animals for each group and time point.

2.1.1. Liver Fibrosis

Carbon tetrachloride (CCl4)-induced liver fibrosis is a classical and robust model of bridging fibrosis in rodent. CCl4 is specifically metabolized in ~50% of hepatocytes in the center of the liver lobule to induce an acute CYP2E1-dependent production of highly reactive tricholoromethyl free radicals, causing centro-lobular necrosis that triggers a wound healing response. When the insult is repeated, successive rounds of wound healing occur, causing peri-central fibrosis and progressively forming fibrotic bridges between vascular structures dissecting the parenchyma [33,34]. Fifteen 7-week-old C57BL/6J male mice were purchased from Janvier (Le-Genest-St-Ile, France). Ten mice received an intraperitoneal injection of CCl4 (800 μL/kg body weight in corn oil) three times a week to induce fibrosis. Among them, five mice were treated for 2 weeks to induce incomplete central bridges (named hereafter 2w CCl4 group) and five mice were injected for 7 weeks to yield severe bridging fibrosis (7w CCl4 group) [33,34]. The remaining five mice received the vehicle only (controls). At the time of sacrifice, 48 h after the last injection, mice where anesthetized. The caudate and the left lateral liver lobes were dissected out, fixed with 4% formaldehyde for 24 h and embedded in paraffin.

2.1.2. Lung Fibrosis

Pulmonary fibrosis is a well-known side effect of bleomycin treatment because the lung expresses extremely low levels of bleomycin hydrolase. In the presence of iron and oxygen, non-hydrolyzed bleomycin generates reactive oxygen species. Subsequent cytotoxic lipid peroxides and oxidized proteins cause tissue injury and parenchymal fibrosis [35]. Fifteen 8-week-old C57BL/6J female mice were purchased from Janvier (Le-Genest-St-Ile, France) and divided in three groups of five animals. In the first group, the mice received a single endotracheal instillation with 0.02 IU bleomycin in 50 µL saline and were sacrificed 28 days after instillation. The unique instillation of bleomycin (named after BLM-acute) is a model known to induce a rapid but transient fibrosis. In the second group, the mice underwent repeated endotracheal instillation with 0.02 IU once-a-week for 8 weeks to promote severe and persistent fibrosis [35,36] (BLM-chronic). The control group was instilled with saline solution only. At the time of sacrifice, mice were anesthetized by an intraperitoneal injection of 20 mg sodium pentobarbital (Certa, Braine-l’Alleud, Belgium). One lung was dissected out, fixed with 4% formaldehyde for 24 h and embedded in paraffin.

2.1.3. Kidney Fibrosis

Kidney fibrosis is a primary dose-limiting side effect of antineoplastic treatment with cisplatin. Cisplatin is predominantly eliminated by the kidney, accumulating in the cells lining the proximal and more distal tubules. Intracellular accumulation of cisplatin generates reactive oxygen species, cytotoxicity, and DNA damage, causing a cortical fibrosis sparing the glomeruli [37]. Kidney fibrosis was generated using a well-established model of repeated low-dose cisplatin, as previously described [38]. Fifteen male FVB/n mice aged 8–10 weeks (Charles River, Brussels, Belgium) were administered a weekly intraperitoneal injection of low-dose cisplatin (7 mg/kg, Sigma-Aldrich, Overijse, Belgium) for 3 or 4 weeks (named after 3w CIS and 4w CIS groups) vs. normal saline only (controls). At the time of sacrifice, mice where anesthetized, one kidney was dissected out, fixed with 4% formaldehyde for 24 h and embedded in paraffin.

2.2. Tissue Preparation and Standardized Protocol for Picrosirius Red Staining

Protocol for PSR staining was performed as described [4,5] with recent optimization [10]. Briefly, 5 µm paraffin sections [20] were first deparaffinized and rehydrated. Slides were immersed in 1.0% phosphomolybdic acid solution for 2 min for pH acclimatization and yellow counterstaining [10], rinsed in water and then incubated in a saturated aqueous picric acid solution containing 0.1% Direct red 80 (Sigma #365548) for 2 h at room temperature. After a brief 2 min wash in 0.01 N HCL and an additional wash in water, slides were dehydrated and mounted with a Dako automated coverslipper. Tissue sections were digitalized using a SCN400 slide scanner (Leica Biosystems, Wetzlar, Germany) at 20× magnification.

2.3. Analysis Workflow

Whole tissue sections were analyzed using the image analysis tool Author™ version 2017.2 (Visiopharm, Hørsholm, Denmark) in order to (1) detect the tissue and assist the definition of histological compartments, (2) detect PSR staining, and (3) quantify this staining. These three steps are described hereafter.

2.3.1. Tissue Recognition and Computer-Assisted Segmentation of Histological Compartments

In all models, the tissue section was automatically defined. Artifacts and regions-of-exclusion (ROEs, listed in Table 1 ) were then automatically discarded and manual correction applied when required (i.e., incomplete or incorrect automatic detection). Total tissue area was quantified. Then, computer-assisted delineation of regions-of-interest (ROIs) was performed for each tissue model. Details and illustrations are provided in Supplementary Material 1. Briefly,

For liver sections, connective tissue was classified as perivascular ROI when it directly surrounded a lumen (with or without erythrocytes) or as bridges when located in the parenchyma. The remaining tissue was considered as parenchyma.

For the lung, large connective tissue surrounding air ducts was discarded since its presence was highly variable between sections, introducing a bias. PSR staining was referred as peri-air ducts or perivascular ROIs when found at a maximum distance of 25 µm from air ducts (bronchi and bronchioles) or from blood vessels. The remaining tissue was considered as parenchyma.

For the kidney, the inner medulla, the outer medulla, and the cortex were manually defined based on morphological criteria, as described [39]. The cortex was further separated divided into large cortical vessels and vessel-free cortex.

Table 1

Characteristics of regions-of-interest (ROIs) and regions-of-exclusion (ROEs).

Total Tissue Area (mm 2 )Regions-of-Interest (ROIs)Regions-of-Exclusion
(ROEs)
Mean (±SD)
Relative Area (%) p-Val
Mean (±SD)
ControlOrgans
(n = 5)
Fibrotic Organs 1
(n = 5)
Liver 43.3 (13.9)Perivascular tissue1.20 (0.32)1.92 (0.70)0.007Capsule
Bridges0.04 (0.02)2.79 (0.30)Large lumens
Parenchyma98.76 (0.30)95.29 (0.84)
Lung 32.77 (6.18)Peri-air ducts tissue4.22 (0.39)4.24 (0.81)0.950Pleura
Perivascular tissue2.94 (0.78)2.93 (0.45)0.992Dense peri-air duct connective tissue (variable) 2
Parenchyma92.85 (0.95)92.83 (1.01)0.973Large lumens & arteries
Kidney 17.63 (1.73)Inner medulla4.46 (1.15)6.49 (3.07)0.204Capsule
Outer medulla10.88 (2.32)13.03 (2.61)0.205Hilum
Cortex84.66 (1.78)80.48 (2.81)0.022Large lumens

1 Last time point of treatment (liver: 7w CCl4; lung: BLM-chronic; kidney: 4w CIS). 2 Strongly stained fibers surrounding very large air ducts were included in the total tissue area but not considered as a ROI. This region represented 5.02% (1.42) and 6.46% (2.12) of the total tissue area (p > 0.05).

ROIs and ROEs are illustrated in Figure 1 . The respective proportion of the area of each ROI vs. total tissue area was calculated as

Relative ROI area = Σ A r e a s e l e c t e d R O I A r e a t o t a l t i s s u e × 100 ( i n % )

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g001.jpg

Illustration of fibrosis development in the murine models of liver (AC), lung (DF) and kidney (GI) fibrosis. Representative samples of control (A), 2w CCl4 (B), and 7w CCl4 (C) livers; control (D), BLM-acute (E), and BLM-chronic (F) lungs; control (G), 3w CIS (H), and 4w CIS (I) kidneys. Global view of entire histological sections (scale bar = 1 mm) and higher magnification of specific regions (inserts).

2.3.2. Optimized Algorithm to Detect PSR-Stained Collagen Deposits and Fiber Segmentation

Two detection methods for PSR staining were designed and illustrated in the Supplement Material 2. The first method (PSRRGB) extracted the red pixels of red-green-blue (RGB) filtering (RGB-G filter). For the second method optimized for increased selectivity and specificity (PSROPT), images were preprocessed using a combination of RGB-G and H-AEC filters to generate an image of the noise filter. This noise was subtracted to red–green contrast filter providing a robust detection of the true signal. Chromaticity red filter completed this detection thanks to its high resolution. Based on this optimized detection, PSR-stained collagen fragments were subsequently segmented (PSRMORF) by watershed and object separation and then classified according to their intensity values and density into strongly stained and plain segments (compact fibers) or weakly stained and loose segments (scattered fibers), in a close manner as defined elsewhere [31]. For each detection method, one threshold was arbitrarily chosen on stained vs. non-stained regions of the three organs to allow multi-organ comparison, as explained in Supplementary Material 3.

2.3.3. PSR-Stained Collagen Deposits Quantification Parameters

Collagen proportionate area (CPA) value is a morphometric measurement of signal density. Values were calculated as conventionally defined [9,23,24,25,26,28,29,30]. When performed on the entire tissue section, the measure was here referred to as in toto CPA [Equation (2)]. When tissue is segmented in multiple ROIs, additional information becomes available such as the distribution of all collagen among the different ROIs, defined here as collagen proportion [CP, Equation (3)], and the proportion of collagen detected in a specific ROI within total tissue, defined here as ROI CPA [Equation (4)]. These measurements were calculated as follows:

I n t o t o CPA = s t a i n e d a r e a t o t a l t i s s u e a r e a t o t a l t i s s u e × 100 ( i n % )

Collagen proportion = Σ s t a i n e d a r e a s e l e c t e d R O I s t a i n e d a r e a t o t a l t i s s u e × 100 ( i n % )

ROI CPA = Σ s t a i n e d a r e a s e l e c t e d R O I a r e a t o t a l t i s s u e × 100 ( i n % )

The same parameters were applied for the quantification of compact and scattered fibers detected with the PSRMORF method. Quantifications were performed at maximal resolution and raw values were expressed as µm 2 .

The analysis workflow is illustrated in Figure 2 .

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g002.jpg

Analysis workflow for in toto vs. regionalized quantification using the three developed detection algorithms: RGB method (PSRRGB), optimized method (PSROPT) or morphological fiber classification (PSRMORF). ROI, region-of-interest; ROE, region-of-exclusion; CPA, collagen proportionate area; CP, collagen proportion.

2.4. Statistics

Statistical analyses were performed using GraphPad Prism 8 on Windows10. One-way ANOVA were performed with Fisher’s post-hoc test. For correlation analysis, coefficient of determination (goodness-of-fit, R 2 ) was computed with a confidence interval of 95%. For the fiber classification, inter-group comparisons were achieved by one-way ANOVA, while intra-group comparisons were computed by paired (CP) or unpaired (ROI CPA) T-tests. The respective control group for each model served for the multiple comparisons. p-value

3. Results

3.1. Visual Examination of the Histological Sections

As expected, specimen examination revealed that CCl4 exposure causes a peri-centrilobular liver fibrosis characterized by thickening of the vascular wall and progressive development of centro-central bridges in a time-dependent manner ( Figure 3 A). In the lungs, instillation of bleomycin induced alveolar hyperplasia (more extensive in the chronic model), patchy alveolar fibrosis, and inflammatory and/or fibroblastic foci often localized close to blood vessels ( Figure 3 B). Cisplatin injections led to the development of a time-dependent fibrosis around proximal and distal tubules in the inner layer of the renal cortex, as well as in the inner and outer medulla ( Figure 3 C).

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g003.jpg

Computer-assisted segmentation of histological structures in the three studied organs. An illustration is given before (left) and after (right) semi-automated segmentation in (A) a 7w CCl4 liver, (B) a BLM-chronic lung, and (C) a 4w CIS kidney. Details in Supplementary Material 1. Scale bar = 1 mm.

3.2. A Robust and Tissue-Transposable Algorithm to Detect PSR Staining in Diverse Fibrotic Organs

The first objective of the present study was to robustly detect PSR staining with an algorithm applicable to diverse organs. In first intention, we applied an RGB-based filtering to uncover PSR-stained pixels. As shown in Figure 4 B, there is an apparent lack of sensitivity as fibers were incompletely detected (blue arrows). Moreover, despite a low detection threshold, non-PSR stained pixels were recognized as staining ( Figure 4 B, black arrow). These errors were more frequent in the livers than in lungs or kidneys.

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g004.jpg

Detection methods of discrete liver fibrosis stained in PSR. (A) Native image. (B) PSRRGB detection method, using RGB-G filter (insert f1). Staining detection is highlighted in red. Blue arrows: undetected pixels (false negative); black arrow: unsaturated pixels (false positive). (C) PSROPT detection method using optimized filters (inserts f2 and f3). Staining detection is highlighted in red. (D) Optimized detection is used for morphological segmentation (based on fiber density and staining intensity) with algorithm PSRMORF. Compact fibers are labeled in red and scattered fibers are labeled in blue. Refer to Supplementary Materials 2 and 3A for the filters, the detailed procedure, and the dynamic ranges of sensitivity. Scale bar = 50 µm.

To improve the specificity and sensitivity of the detection, we developed an optimized algorithm based on a more robust filtering of stained pixels (PSROPT method). This preprocessing increased the specificity for PSR staining vs. background, allowing lower detection threshold and more dynamic range, thereby also improving sensitivity of the detection ( Figure 4 C and Supplementary Material 2D).

The specificity and sensitivity of PSRRGB and PSROPT detection methods were assessed by comparing in toto CPA values between organs at baseline and between the different time points within the same model. In accordance with the visual appraisal of staining at baseline (see Figure 3 ), in toto PSROPT CPA were ~4× higher in the control lungs than in the livers and the kidneys ( Figure 5 B and Supplementary Material 3B), and this difference was lower with the PSRRGB method ( Figure 5 A). Comparing in toto CPA in the treated groups vs. their respective controls, PSRRGB method evidenced fibrosis only at the latest time point of treatment in livers and kidneys and with chronic bleomycin treatment in the lungs ( Figure 5 A). By contrast, fibrosis at the early time point and progressive fibrosis were evidenced in the three organs using PSROPT ( Figure 5 B).

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g005.jpg

In toto CPA analyses using PSRRGB (A) vs. PSROPT (B) filtering method for signal detection. Quantification performed on entire tissue sections (ROEs excluded). One-way ANOVA with Fisher’s post-hoc test. a: p < 0.050, b: p < 0.010, c: p < 0.001 vs. controls; *: p < 0.050, $: p < 0.001 last time point vs. intermediate time point.

Correlation between the two methods was modest with a correlation coefficient weaker for the liver sections (0.499) as compared with the lungs (0.647) and the kidneys (0.693) (Supplementary Material 3C), tentatively explained by a higher noise in the liver fibrosis model.

These results suggest that the PSRRGB detection method is not sensitive enough to highlight collagen deposits accumulation during fibrosis onset in these models, and this limitation could be circumvented, at least partially, by the PSROPT method.

3.3. Regionalized Analysis Revealed Topical Fibrosis Changes

The second aim was to evaluate whether regionalized analysis would allow a better characterization of fibrosis. To meet this goal, histologically relevant ROIs were defined for each organ ( Figure 1 and Supplementary Material 1). The list of ROIs and ROEs and the relative area of ROIs are given in Table 1 . PSR staining was then quantified in each ROIs using the PSROPT detection method. Results were expressed as CP (distribution of collagen among different ROIs) and ROI CPA (percentage of collagen detected in a specific ROI within total tissue).

3.3.1. Liver Fibrosis

While in normal livers, parenchyma represented more than 95% of the tissue sections ( Table 1 ), PSR staining in this region only corresponded to ~10% of the total collagen amount ( Table 2 ). The greatest CP was detected in the perivascular areas. This proportion decreased proportionally as the bridges formed ( Table 2 ). Although the ROI CPA moderately increased around vascular structures in fibrotic livers (~2 times), this increase was much larger and time-dependent in bridges (~70 times at 2 weeks and ~130 times at 7 weeks). Collagen accumulation in the parenchyma and outside the fibrous septa was only detectable in the severe fibrotic condition ( Figure 6 ).

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g006.jpg

Liver ROI CPA calculated for each ROI using the PSROPT detection method. Mean (±SD) and one-way ANOVA with Fisher’s post-hoc test. a: p < 0.050, b: p < 0.010, c: p < 0.001 vs. controls; $: p < 0.001 vs. middle group.

Table 2

Liver collagen proportion per ROI. Mean percentages (±SD).

Collagen Proportion (%)Controls2w CCl4p-Val7w CCl4p-Val
Perivascular85.71 (11.10)64.22 (7.71)0.00645.06 (11.40)
Bridges1.89 (1.94)32.09 (7.72)46.27 (8.44)
Parenchyma12.40 (9.24)3.69 (1.24)0.0378.67 (4.04)0.334
Total100%100% 100%

3.3.2. Lung Fibrosis

As shown in Table 1 , the distribution of the ROIs (peri-air ducts, perivascular tissue, and parenchyma) was similar in the normal and bleomycin-treated tissues. In the healthy lung, the greatest CP was detected in the parenchyma (~75%), where it further increased after single or repeated bleomycin instillation ( Table 3 ). The ROI CPA around air ducts was unchanged after fibrosis induction, but bleomycin caused a significant enrichment of PSR staining in the perivascular region ( Figure 7 ).

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g007.jpg

Lung ROI CPA calculated for each ROI using the PSROPT detection method. Mean (±SD) and one-way ANOVA with Fisher’s post-hoc test. b: p < 0.010, c: p < 0.001 vs. controls.

Table 3

Lung collagen proportion per ROI. Mean percentages (±SD).

Collagen Proportion (%)ControlsBLM-Acutep-ValBLM-Chronicp-Val
Perivascular5.46 (1.58)4.89 (1.84)0.5565.10 (0.80)0.709
Peri-air ducts19.64 (5.45)8.66 (2.61)0.0019.81 (1.43)
Parenchyma74.91 (6.39)86.44 (4.33)0.00285.09 (1.74)0.004
Total100%100% 100%

3.3.3. Kidney Fibrosis

Under normal condition, the cortex was the dominant ROI (85% of the tissue) ( Table 1 ) and the most important source of PSR signal (~85%); a large part of it was attributed to vessels ( Table 4 ). Cisplatin treatment induced a reduction of the relative cortical area ( Table 1 ), a time-dependent increase in the CP in the cortex (not attributed to vessels) and an elevation of the CP within the outer medulla ( Table 4 ). The ROI CPA was increased in the inner medulla, outer medulla and cortex upon cisplatin administration ( Figure 8 A), with the largest accumulation in the (vessel-free) cortex ( Figure 8 B). The results of in toto CPA were validated independently, by using a previously described DIA method for PSR analysis [40] (Supplementary Material 4).

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g008.jpg

Kidney ROI CPA calculated using the PSROPT detection method for each ROI (A) and in sub-analyzing the cortex (B). Mean (±SD) and one-way ANOVA with Fisher’s post-hoc test. b: p < 0.010, c: p < 0.001 vs. controls; *: p < 0.050, #: p < 0.010, $: p < 0.001 last time point vs. intermediate time point.

Table 4

Kidney collagen proportion per ROI. Mean percentages (±SD). Cortex was sub-analyzed to distinguish the parenchyma (1) and the large cortical vessels (2) .

Collagen Proportion (%)Controls3w CISp-Val4w CISp-Val
Inner medulla10.29 (7.95)15.83 (3.78)0.14714.54 (4.28)0.257
Outer medulla4.94 (1.07)13.06 (3.49)4.98 (1.62)0.979
Cortex 1,2 84.78 (8.59)71.10 (4.02)0.00480.48 (4.60)0.286
Total100%100% 100%
1 Cortex (minus vessels)36.07 (8.00)54.90 (4.57)0.00165.33 (8.51)
2 Large cortical vessels48.70 (10.52)16.21 (7.87)15.15 (6.46)

3.4. Morphological Parameters for Fibrosis Classification

Our third aim was to provide a tool for the automated categorization of collagen deposits based on the morphological organization of the fibers, as proposed. Using the same detection filter than PSROPT, a new algorithm (PSRMORF) was developed to automatically segment the PSR fibers by watershed and separate objects, and quantify collagen fibers defined as compact or scattered based on the mean signal intensity of each segment ( Figure 4 D) (details given in Supplement Material 2B,D).

3.4.1. In Toto Analysis of Fiber Morphology

Because the fiber classification by PSRMORF was obtained from PSROPT algorithm, the sum of the two types of fibers calculated in toto is equal to the values described above ( Figure 5 bottom panels). Compact fibers were dominant in control livers and kidneys, representing respectively ~75% and 60% of total collagen ( Figure 9 A,C), while the normal lung was proportionally richer in scattered fibers, ( Figure 9 B), at least after exclusion of the dense connective tissue surrounding air ducts as detailed above. Fibrosis in the liver was characterized by a proportional increase in both compact and scattered fibers. In contrast, bleomycin and cisplatin caused a preferential accumulation in scattered fibers in lungs and kidneys.

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g009.jpg

In toto analysis of compact and scattered fibers using PSRMORF detection method, performed in total tissue for (A) the livers, (B) the lungs and (C) the kidneys. Mean (±SD) and one-way ANOVA with Fisher’s post-hoc test. a: p < 0.050, b: p < 0.010, c: p < 0.001 vs. controls; #: p < 0.010, $: p < 0.001 last time point vs. intermediate time point.

3.4.2. Regionalized Analysis of Fiber Morphology

Spatiotemporal fibrotic changes were investigated by calculating the evolution of the relative proportion of compact and scattered fibers in each ROI ( Table 5 ) and of their abundance in the tissue ( Figure 10 ).

An external file that holds a picture, illustration, etc. Object name is biomolecules-10-01585-g010.jpg

Regionalized analysis of compact and scattered fibers with PSRMORF detection method. Results are shown in the liver for the perivascular region (A), bridges (B), and parenchyma (C), in the lung for the perivascular region (D), the peri-air duct region (E), and the parenchyma (F), and in the kidney for the inner medulla (G), the outer medulla (H) and the cortex (vessels included) (I). Mean (±SD). One-way ANOVA with Fisher’s post-hoc test. a: p < 0.050, b: p < 0.010, c: p < 0.001 vs. controls; *: p < 0.050, #: p < 0.010, $: p < 0.001 last time point vs. intermediate time point.

Table 5

Collagen proportion (%) per fiber type in specific histological compartments. Mean percentages (±SD).

Compact/Scattered Fibers Relative Area
Regions-of-Interest
(ROIs)
Fiber ClassControls
(n = 5)
Treatment Time-Point 1
(n = 5)
p-ValTreatment Time-Point 2
(n = 5)
p-Val
Liver Perivascular tissuecompact81.23 (7.82)87.41 (3.70)0.10086.00 (3.92)0.194
scattered18.77 (7.82)12.59 (3.70)14.00 (3.92)
Bridgescompact66.78 (9.28)71.29 (4.61)0.32867.00 (6.27)0.961
scattered32.22 (9.28)28.71 (4.61)33.00 (6.27
Parenchymacompact23.15 (3.32)18.90 (3.55)0.27223.49 (8.87)0.928
scattered76.85 (3.32)81.10 (3.55)76.51 (8.87)
Lung Peri-air ducts tissuecompact59.61 (7.03)41.53 (5.10)48.07 (5.71)0.010
scattered40.39 (7.03)58.47 (5.10)51.93 (5.71)
Perivascular tissuecompact45.53 (7.47)36.44 (5.40)0.02642.81 (3.32)0.461
scattered54.47 (7.47)63.56 (5.40)57.19 (3.32)
Parenchymacompact36.32 (6.14)33.36 (3.64)0.29841.45 (2.06)0.083
scattered63.68 (6.14)66.64 (3.64)58.55 (2.06)
Kidney Inner medullacompact6.04 (4.40)7.20 (2.67)0.6969.95 (6.04)0.203
scattered93.96 (4.40)92.80 (2.67)90.05 (6.04)
Outer medullacompact10.81 (7.42)18.70 (10.41)0.15415.48 (6.15)0.385
scattered89.19 (7.42)81.30 (10.41)84.52 (6.15)
Cortexcompact68.47 (3.48)39.48 (7.63)42.25 (6.01)
scattered31.53 (3.48)60.52 (7.63)54.75 (6.01)

In livers, the compact fibers mostly located in the perivascular connective tissue and the bridges, while scattered fibers were predominant in the parenchyma, irrespectively to the fibrotic status ( Table 5 ). The modest increase in collagen amount detected in the perivascular area of CCl4-treated livers was attributed to a local enrichment in compact fibers ( Figure 10 A). The bridges showed gradual increase in the abundance of both fiber types upon CCl4 treatment, with compact fibers being twice predominant over the scattered fibers ( Figure 10 B). Collagen fibers that accumulated in the parenchyma after 7 weeks of treatment displayed mostly a scattered phenotype ( Figure 10 C).

In normal lungs, all compartments comprised both compact and scattered fibers. Upon fibrosis induction, higher proportion of scattered fibers was detected around the air ducts (both dosing regimen) and the vessels (BLM-acute only) ( Table 5 ). In the parenchyma, bleomycin did not change the fiber proportion (two-thirds scattered, one-third compact) ( Table 5 ), presumably because both fiber types accumulated in this ROI ( Figure 10 F). The two regimens increased the abundance of scattered fibers around vessels whereas a net increase in compact fibers was detected only in the chronic model ( Figure 10 D). The area surrounding air ducts remained essentially unaffected ( Figure 10 E).

In normal kidneys, the inner and outer medullas were mostly made up of scattered fibers, while the cortex was enriched in compact fibers. Cisplatin treatment did not change these proportions in the medulla but reversed the balance in the cortex ( Table 5 ). Looking at the regional abundance per fiber type of fiber, cisplatin increased collagen deposits in all compartments, especially in the form of scattered fibers. Compact fibers also accumulated in all ROIs of cisplatin-treated kidneys compared to controls, particularly in the cortex in a time-dependent manner ( Figure 10 G–I), pointing to the cortex as a primary site of dense interstitial fibrosis.

4. Discussion

Whole slide imaging and modern DIA at high resolution enlarge the possibilities for accurate quantification. In the present study, we exploited DIA to address the said limitations of analyzing PSR stained sections. We studied three organs in which assessment of fibrosis is of clinical relevance. For this, we developed algorithms to (1) focus the analysis on parenchymal compartments and exclude physiologically collagen-rich structures and (2) quantitatively describe the morphological organization of PSR-stained collagen fibers, here defined as compact or scattered.

4.1. Novelty, Advantages, and Limitations of the Proposed Method

In line with previous comparisons of red-component (RGB-G) vs. color-specific threshold for PSR staining analysis [22], we found the RGB-based method relatively poorly sensitive, and particularly for detecting early or faint fibrosis. The PSROPT algorithm was thus optimized to be more sensitive and more specific to PSR staining, and minimally dependent on the thresholding. Similar to [20] on severe cardiac fibrosis, our method exploited the red-green contrast filter for PSR detection. However, thanks to the addition of other preprocess filters, the dynamic range of detection of PSROPT was not only (i) more robust to naturally occurring intra- and inter-studies variations, but also (ii) ampler than that of PSRRGB method (illustrated in the Supplementary Material 2). Indeed, the use of hue-saturation-value (from which derives chromaticity filtering) can explain the gain in robustness to staining variations and optimal detection of faint signal [41].

The analysis of the entire tissue (in toto results) is a fast, simple and reliable tool for global detection of accumulated collagen. However, thanks to the regionalized quantifications, we can reveal additional changes, unsuspected or hidden in the global analysis. We provide thereby a tool to (i) evaluate whether fibrosis occurs in a regionalized manner or not, (ii) quantitatively document the relative distribution of collagen according to the histological structure, and (iii) quantify local collagen accumulation in the specific ROI(s). It is thus now possible to finely characterize fibrosis specificities in various organs using the same methodology (as discussed below).

Finally, our method offers a morphometric analysis of the organization of collagen deposits. Beside collagen content, automated fiber segmentation enables to distinguish small intertwined fibers from large clumps, each collagen segment being classified according to the mean intensity value of all the pixels that compose it. Thus, the segmentation in ROIs and the classification in fiber types offer a refined description of fibrosis for assessing spatiotemporal changes in a finely tuned manner [21].

Our method has some limitations:

Intrinsically, the analysis of tissue section is an incomplete representation of the three-dimensional structures.

Despite the computer-assisted recognition of the structures, the validation of the ROIs remains essential. This manual correction is a non-negligible source of error.

The outcomes are defined by an arbitrary-defined threshold, based on visual evaluation. The developed PSROPT algorithm limits the impact of adjusting threshold. Robustly standardized histological procedures are still essential before quantitative DIA.

The classification into compact and scattered fibers was not based on the collagen nature or effective physical parameters but on the rationale estimation that denser or thicker collagen fibrils contain more basic residues binding the dye and are therefore more intensively stained.

4.2. Biological Results

4.2.1. Liver Fibrosis

First, the progressive 2- and 3-fold accumulation of collagen deposits found after 2 weeks or 7 weeks exposure to CCl4 was in line with existing literature [42,43]. In accordance with previous analyses in the same model [44,45], our regionalized analysis located most of the accumulated collagen in the bridges. This local enrichment of collagen in the bridges could be attributed to both compact (predominant) and scattered fibers. This agrees with earlier reports showing a relative increase in aggregated fibers with extended exposure to CCl4 [44]. Bridges have been shown to be rich in type I, III and IV collagen [42,45], with a prominent role of type III and IV during fibrosis progression [19,45]. Since PSR indifferently stains these three types of collagens [5,6,46], our results should reflect any change in the sum of these different fibers. Besides bridges, we also evidenced larger amount of collagen deposits in the parenchyma after long-term treatment (7 weeks), potentially corresponding to small or early forming deposits. Usually not integrated in the pathological descriptors, this information obtained by regionalized analysis however indicates a worsening of the disease. Together, our data confirm the progressive local accumulation of collagen in the bridges in response to CCl4 in a time-dependent manner, starting as scattered fibers and evolving to a more compact phenotype. Whether the distribution of scattered vs. compact fibers in the bridges and parenchyma reflects the enrichment in one type of collagen and/or reflects disease progression or resolution remains to be elucidated.

We can conclude that CCl4 hepatotoxicity is characterized by the accumulation of dense collagen fibers around vessels, and a strong induction of the bridging fibrosis that progresses with exposure time.

4.2.2. Lung Fibrosis

Bleomycin, irrespectively of the instillation protocol, induced a significant increase in collagen. The histological aspect of bleomycin-treated lungs fairly recapitulates the observations of the fundamental study by Degryse and Lawson [36]. The authors noted patchy areas of parenchymal fibrosis in the acute model, while the chronic model was characterized by coexistence of normal lung, dense scar, active remodeling with interstitial collagen deposition, and alveolar hyperplasia. In our specimens, we saw no dense scar, while we noted foci located close to vessels in a variable proportion between samples. The regionalized analysis pointed to a predominant parenchymal fibrosis in bleomycin-treated lungs, irrespective of the administration protocol but with a greater extend in the repetitive model. This interstitial fibrosis was a mixture of compact and loose fibers, with an increasing role of the compact fibers in the chronic model, suggesting a more severe status after repeated instillation [47]. Lung fibrotic regions have been shown to be rich in type IV collagen [48], which does not organize in dense beams of fibers, possibly explaining the predominance of scattered fibers in our study. Our results also suggest a local effect in the direct vicinity of blood vessels. First, we noted an accumulation of fibro-inflammatory foci close to the vasculature. Second, the two-classes segmentation of fibers showed a perivascular enrichment in scattered fibers in all treated lungs and an accumulation of dense fibers in this region only in chronic-treated lungs. This suggests that fibrosis induced by chronic exposure to bleomycin is more compact, perhaps more mature, and/or less subject to resorption than after a unique instillation [36]. Our results may support the concept of vascular alterations in the initiation of interstitial fibrosis and disease progression, as already proposed [49,50]. Interestingly, the formation of vessel-related structures seen by medical tomography in fibrotic patients can predict the likelihood of decline of the forced vital capacity [51] underlying the key role of vasculature in the disease. Bleomycin-induced fibrosis causes bronchial dilation [52], a phenomenon also seen in patients. In our study, while the collagen amount at the margin of air ducts remained unaffected in fibrotic samples, an increase proportion of scattered fibers was noted in this ROI. It is possible that fibers are mechanically stretched, and therefore appear less dense.

In conclusion, bleomycin induces a parenchymal lung fibrosis, and our data support a pathological action on the vasculature.

4.2.3. Kidney Fibrosis

In patients, chronic use of cisplatin causes significant tubular atrophy and interstitial fibrosis with relative glomerular sparing. The repeated low-dose model used here mimics this condition more closely than the standard models (with unique and higher dose) in terms of effect on kidney function, kidney damage (interstitial cortical fibrosis) and survival rate [38].

The samples used here recapitulated the 3-fold increase in collagen deposition after 4 weeks of cisplatin as documented [38].

However, our sensitive method allowed to detect fibrosis as early as 3 weeks after the initiation of the cisplatin treatment. At this time-point, collagen accumulated in the cortex, essentially in the form of scattered fibers. One week later as fibrosis progressed, fibers were more compact, with no impact on large cortical vessels. These observations support the current understanding of cortical uptake and tubular alteration by cisplatin leading to fibrosis in these regions, despite vasculature injury [37]. Both the inner and outer medulla displayed collagen accumulation, in line with recent untargeted metabolomics data showing higher sensitivity of the medulla to cisplatin than the cortex [53]. This increase was characterized by a strong accumulation of scattered fibers, but also of compact fibers, suggesting fibrosis maturation similarly to the cortex.

4.3. Perspectives for Digital Pathology and Research

We demonstrate here the relevance and reproducibility of analyzing the local collagen accumulation in specific tissue compartments, a fortiori in the form of compact or scattered fibers with accessible material (paraffin sections), using a standard method (PSR staining) and a standard imaging system (brightfield slide scanner). Nevertheless, to go further in this characterization, complementary techniques are needed.

First, as illustrated in this work, an unbiased definition of the structures of interest for the different organs is needed, especially owing the possibility to fully automatize this key step by artificial intelligence [17,54,55].

Second, PSR stained sections can then be visualized with brightfield microscope but also under circular polarized light (although polarization data remain matter to debate for their interpretation) or fluorescent illumination, the latest possibly being as sensitive as the second harmonic generation for collagen detection [56,57,58]. Therefore, the validation of our morphometric classification of the compact vs. scattered fibers may be addressed by one of these relatively accessible methods in the future. Alternatively, label-free second-harmonic generation/two-photon excited fluorescence microscopy may address this question. Indeed, when combined with deep learning, this is used to reveal the spatiotemporal changes in liver fibrosis [27,31,32,44,59,60]. This technology is promising since the recent marketing of a slide scanner using this microscopy method for clinical applications [61]. Furthermore, another type of advance microscopy (staining-free fluorescence lifetime imaging) is proposed to distinguish type I vs. type III collagen [62], but this claim still needs validation by complementary approaches such as immunohistology. With such technologies, additional morphometric features can be investigated such as the fiber length, diameter, orientation, alignment, anisotropy, or shape factors (i.e., waviness) which have been associated with pathologies, developing the concept of pathology-associated collagen signatures [63,64]. Such classification could ultimately serve to stratify the degree of fiber density and compactness using conventional microscopy.

5. Conclusions

The golden rules for a reliable method are simplicity (easy to handle, low manpower need, etc.) and reproducibility within and between studies. Our new digital method to detect collagen fibers using PSR staining in brightfield is simple, robust, efficient, and tissue-transposable, allowing inter-studies and multi-organ comparison. The regionalized analysis offers a spatiotemporal characterization of collagen deposit organization and the proposed method of fiber classification enables pattern recognition to better describe fibrosis characteristics. Beside the useful place of PSR for accurate diagnosis, the present work implicates that the analysis opportunities offered by this staining are often underexploited.

Acknowledgments

The authors thank Natacha Feza-Bingi, Yousof Yakoub, Yvette Cnops, and Sebastien Druart for their expertise in animal models. We are indebted to Michele de Beuckelaer for her help in histological preparation and to Aurélie Daumerie for her advice with Visiopharm.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-273X/10/11/1585/s1. Supplementary Material 1: Computer assisted organ-specific definition of the regions-of-interest. Supplementary Material 2: Detection filters and threshold determination for the RGB-based or the PSR-optimized methods. Supplementary Material 3: RGB-based and the PSR-optimized methods used for multi-organ comparison in liver, lung, and kidney (illustration and correlation). Supplementary Material 4: Validation of the developed PSR-optimized method vs. previously described PSR-detection method on murine kidney samples.

Author Contributions

Conceptualization, G.E.C. and C.B.; Methodology, G.E.C. and C.B.; Software, G.E.C.; Validation, C.B., I.L., A.F., G.S., J.M., O.D., and F.H.; Formal analysis, G.E.C. and G.S.; Data curation, G.E.C. and C.B.; Writing—original draft preparation, G.E.C.; Writing—review and editing, all authors; Supervision, C.B., I.L., O.D., J.M., and F.H. All authors have read and agreed to the published version of the manuscript.

Funding

G.E.C. is supported by the Fonds de Recherche Clinique (Saint-Luc); J.M. is supported by the Fonds National pour la Recherche Scientifique and by the Fondation Saint-Luc; G.S. and O.D. are funded by the Swiss National Science Foundation (310030_189044). No other funding to disclose.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Rosenbloom J., Macarak E., Piera-Velazquez S., Jimenez S.A. Human Fibrotic Diseases: Current Challenges in Fibrosis Research. Methods Mol. Biol. 2017; 1627 :1–23. doi: 10.1007/978-1-4939-7113-8_1. [PubMed] [CrossRef] [Google Scholar]

2. Kliment C.R., Englert J.M., Crum L.P., Oury T.D. A novel method for accurate collagen and biochemical assessment of pulmonary tissue utilizing one animal. Int. J. Clin. Exp. Pathol. 2011; 4 :349–355. [PMC free article] [PubMed] [Google Scholar]

3. Coentro J.Q., Capella-Monsonis H., Graceffa V., Wu Z., Mullen A.M., Raghunath M., Zeugolis D.I. Collagen Quantification in Tissue Specimens. Methods Mol. Biol. 2017; 1627 :341–350. doi: 10.1007/978-1-4939-7113-8_22. [PubMed] [CrossRef] [Google Scholar]

4. Junqueira L.C., Bignolas G., Brentani R.R. Picrosirius staining plus polarization microscopy, a specific method for collagen detection in tissue sections. Histochem. J. 1979; 11 :447–455. doi: 10.1007/BF01002772. [PubMed] [CrossRef] [Google Scholar]

5. Dapson R.W., Fagan C., Kiernan J.A., Wickersham T.W. Certification procedures for sirius red F3B (CI 35780, Direct red 80) Biotech. Histochem. 2011; 86 :133–139. doi: 10.3109/10520295.2011.570277. [PubMed] [CrossRef] [Google Scholar]

6. Lattouf R., Younes R., Lutomski D., Naaman N., Godeau G., Senni K., Changotade S. Picrosirius red staining: A useful tool to appraise collagen networks in normal and pathological tissues. J. Histochem. Cytochem. 2014; 62 :751–758. doi: 10.1369/0022155414545787. [PubMed] [CrossRef] [Google Scholar]

7. Whittaker P., Kloner R.A., Boughner D.R., Pickering J.G. Quantitative assessment of myocardial collagen with picrosirius red staining and circularly polarized light. Basic Res. Cardiol. 1994; 89 :397–410. doi: 10.1007/BF00788278. [PubMed] [CrossRef] [Google Scholar]

8. Huang Y., de Boer W.B., Adams L.A., MacQuillan G., Rossi E., Rigby P., Raftopoulos S.C., Bulsara M., Jeffrey G.P. Image analysis of liver collagen using sirius red is more accurate and correlates better with serum fibrosis markers than trichrome. Liver Int. 2013; 33 :1249–1256. doi: 10.1111/liv.12184. [PubMed] [CrossRef] [Google Scholar]

9. Street J.M., Souza A.C., Alvarez-Prats A., Horino T., Hu X., Yuen P.S., Star R.A. Automated quantification of renal fibrosis with Sirius Red and polarization contrast microscopy. Physiol. Rep. 2014; 2 doi: 10.14814/phy2.12088. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

10. Hatori M., Moriya S., Fujimori M., Kobayashi S., Ikota H., Shirabe K., Yokoo H., Kimura K., Saio M. Phosphomolybdic Acid Prevents Nonspecific Nuclear Staining by Picrosirius Red but Is Converted to Molybdenum Blue by Blue Light. J. Histochem. Cytochem. 2020; 68 :621–634. doi: 10.1369/0022155420942620. [PubMed] [CrossRef] [Google Scholar]

11. Segnani C., Ippolito C., Antonioli L., Pellegrini C., Blandizzi C., Dolfi A., Bernardini N. Histochemical Detection of Collagen Fibers by Sirius Red/Fast Green Is More Sensitive than van Gieson or Sirius Red Alone in Normal and Inflamed Rat Colon. PLoS ONE. 2015; 10 :e0144630. doi: 10.1371/journal.pone.0144630. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

12. van Steenberghe M., Schubert T., Gerelli S., Bouzin C., Guiot Y., Xhema D., Bollen X., Abdelhamid K., Gianello P. Porcine pulmonary valve decellularization with NaOH-based vs detergent process: Preliminary in vitro and in vivo assessments. J. Cardiothorac. Surg. 2018; 13 :34. doi: 10.1186/s13019-018-0720-y. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

13. Gadd V.L. Combining immunodetection with histochemical techniques: The effect of heat-induced antigen retrieval on picro-Sirius red staining. J. Histochem. Cytochem. 2014; 62 :902–906. doi: 10.1369/0022155414553667. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

14. Rousselet M.C., Michalak S., Dupre F., Croue A., Bedossa P., Saint-Andre J.P., Cales P., Hepatitis N. Sources of variability in histological scoring of chronic viral hepatitis. Hepatology. 2005; 41 :257–264. doi: 10.1002/hep.20535. [PubMed] [CrossRef] [Google Scholar]

15. Arjmand A., Tsipouras M.G., Tzallas A.T., Forlano R., Manousou P., Giannakeas N. Quantification of Liver Fibrosis—A Comparative Study. Appl. Sci. 2020; 10 :447. doi: 10.3390/app10020447. [CrossRef] [Google Scholar]

16. Yu Y., Wang J., Ng C.W., Ma Y., Mo S., Fong E.L.S., Xing J., Song Z., Xie Y., Si K., et al. Deep learning enables automated scoring of liver fibrosis stages. Sci. Rep. 2018; 8 :16016. doi: 10.1038/s41598-018-34300-2. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

17. Seger S., Stritt M., Vezzali E., Nayler O., Hess P., Groenen P.M.A., Stalder A.K. A fully automated image analysis method to quantify lung fibrosis in the bleomycin-induced rat model. PLoS ONE. 2018; 13 :e0193057. doi: 10.1371/journal.pone.0193057. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

18. Kolachalama V.B., Singh P., Lin C.Q., Mun D., Belghasem M.E., Henderson J.M., Francis J.M., Salant D.J., Chitalia V.C. Association of Pathological Fibrosis With Renal Survival Using Deep Neural Networks. Kidney Int. Rep. 2018; 3 :464–475. doi: 10.1016/j.ekir.2017.11.002. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

19. Vassiliadis E., Larsen D.V., Clausen R.E., Veidal S.S., Barascuk N., Larsen L., Simonsen H., Silvestre T.S., Hansen C., Overgaard T., et al. Measurement of CO3–610, a potential liver biomarker derived from matrix metalloproteinase-9 degradation of collagen type iii, in a rat model of reversible carbon-tetrachloride-induced fibrosis. Biomark. Insights. 2011; 6 :49–58. doi: 10.4137/BMI.S6347. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

20. Schipke J., Brandenberger C., Rajces A., Manninger M., Alogna A., Post H., Muhlfeld C. Assessment of cardiac fibrosis: A morphometric method comparison for collagen quantification. J. Appl. Physiol. 2017; 122 :1019–1030. doi: 10.1152/japplphysiol.00987.2016. [PubMed] [CrossRef] [Google Scholar]

21. Flores Molina M., Fabre T., Cleret-Buhot A., Soucy G., Meunier L., Abdelnabi M.N., Belforte N., Turcotte S., Shoukry N.H. Visualization, Quantification, and Mapping of Immune Cell Populations in the Tumor Microenvironment. J. Vis. Exp. 2020:e60740. doi: 10.3791/60740. [PubMed] [CrossRef] [Google Scholar]

22. Miles C., Westaby J., Ster I.C., Asimaki A., Boardman P., Joshi A., Papadakis M., Sharma S., Behr E.R., Sheppard M.N. Morphometric characterization of collagen and fat in normal ventricular myocardium. Cardiovasc. Pathol. 2020; 48 :107224. doi: 10.1016/j.carpath.2020.107224. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

23. Calvaruso V., Burroughs A.K., Standish R., Manousou P., Grillo F., Leandro G., Maimone S., Pleguezuelo M., Xirouchakis I., Guerrini G.P., et al. Computer-assisted image analysis of liver collagen: Relationship to Ishak scoring and hepatic venous pressure gradient. Hepatology. 2009; 49 :1236–1244. doi: 10.1002/hep.22745. [PubMed] [CrossRef] [Google Scholar]

24. Hall A.R., Tsochatzis E., Morris R., Burroughs A.K., Dhillon A.P. Sample size requirement for digital image analysis of collagen proportionate area in cirrhotic livers. Histopathology. 2013; 62 :421–430. doi: 10.1111/his.12010. [PubMed] [CrossRef] [Google Scholar]

25. Tsochatzis E.A., Bosch J., Burroughs A.K. Liver cirrhosis. Lancet. 2014; 383 :1749–1761. doi: 10.1016/S0140-6736(14)60121-5. [PubMed] [CrossRef] [Google Scholar]

26. Chen H.H., Waghorn P.A., Wei L., Tapias L.F., Schuhle D.T., Rotile N.J., Jones C.M., Looby R.J., Zhao G., Elliott J.M., et al. Molecular imaging of oxidized collagen quantifies pulmonary and hepatic fibrogenesis. JCI Insight. 2017; 2 doi: 10.1172/jci.insight.91506. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

27. Xu S., Wang Y., Tai D.C.S., Wang S., Cheng C.L., Peng Q., Yan J., Chen Y., Sun J., Liang X., et al. qFibrosis: A fully-quantitative innovative method incorporating histological features to facilitate accurate fibrosis scoring in animal model and chronic hepatitis B patients. J. Hepatol. 2014; 61 :260–269. doi: 10.1016/j.jhep.2014.02.015. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

28. Farris A.B., Alpers C.E. What is the best way to measure renal fibrosis?: A pathologist’s perspective. Kidney Int. Suppl. 2014; 4 :9–15. doi: 10.1038/kisup.2014.3. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

29. Reed C.L., Sharma S.G., Prieto S.P., Muldoon T.J. Quantitative Analysis of Renal Fibrosis Using a Colorimetric System; Proceedings of the 2018 IEEE Conference on Multimedia Information Processing and Retrieval (MIPR); Miami, FL, USA. 10–12 April 2018; pp. 178–183. [CrossRef] [Google Scholar]

30. Akam E.A., Abston E., Rotile N.J., Slattery H.R., Zhou I.Y., Lanuti M., Caravan P. Improving the reactivity of hydrazine-bearing MRI probes for in vivo imaging of lung fibrogenesis. Chem. Sci. 2020; 11 :224–231. doi: 10.1039/C9SC04821A. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

31. Wang B., Sun Y., Zhou J., Wu X., Chen S., Wu S., Liu H., Wang T., Ou X., Jia J., et al. Advanced septa size quantitation determines the evaluation of histological fibrosis outcome in chronic hepatitis B patients. Mod. Pathol. 2018; 31 :1567–1577. doi: 10.1038/s41379-018-0059-x. [PubMed] [CrossRef] [Google Scholar]

32. Liu F., Goh G.B., Tiniakos D., Wee A., Leow W.Q., Zhao J.M., Rao H.Y., Wang X.X., Wang Q., Wan W.K., et al. qFIBS: An Automated Technique for Quantitative Evaluation of Fibrosis, Inflammation, Ballooning, and Steatosis in Patients With Nonalcoholic Steatohepatitis. Hepatology. 2020; 71 :1953–1966. doi: 10.1002/hep.30986. [PubMed] [CrossRef] [Google Scholar]

33. Weber L.W., Boll M., Stampfl A. Hepatotoxicity and mechanism of action of haloalkanes: Carbon tetrachloride as a toxicological model. Crit. Rev. Toxicol. 2003; 33 :105–136. doi: 10.1080/713611034. [PubMed] [CrossRef] [Google Scholar]

34. Delire B., Starkel P., Leclercq I. Animal Models for Fibrotic Liver Diseases: What We Have, What We Need, and What Is under Development. J. Clin. Transl. Hepatol. 2015; 3 :53–66. doi: 10.14218/JCTH.2014.00035. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

35. Liu T., De Los Santos F.G., Phan S.H. The Bleomycin Model of Pulmonary Fibrosis. Methods Mol. Biol. 2017; 1627 :27–42. doi: 10.1007/978-1-4939-7113-8_2. [PubMed] [CrossRef] [Google Scholar]

36. Degryse A.L., Lawson W.E. Progress toward improving animal models for idiopathic pulmonary fibrosis. Am. J. Med. Sci. 2011; 341 :444–449. doi: 10.1097/MAJ.0b013e31821aa000. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

37. Manohar S., Leung N. Cisplatin nephrotoxicity: A review of the literature. J. Nephrol. 2018; 31 :15–25. doi: 10.1007/s40620-017-0392-z. [PubMed] [CrossRef] [Google Scholar]

38. Sharp C.N., Doll M.A., Dupre T.V., Shah P.P., Subathra M., Siow D., Arteel G.E., Megyesi J., Beverly L.J., Siskind L.J. Repeated administration of low-dose cisplatin in mice induces fibrosis. Am. J. Physiol. Ren. Physiol. 2016; 310 :F560–F568. doi: 10.1152/ajprenal.00512.2015. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

39. Jouret F., Bernard A., Hermans C., Dom G., Terryn S., Leal T., Lebecque P., Cassiman J.J., Scholte B.J., de Jonge H.R., et al. Cystic fibrosis is associated with a defect in apical receptor-mediated endocytosis in mouse and human kidney. J. Am. Soc. Nephrol. 2007; 18 :707–718. doi: 10.1681/ASN.2006030269. [PubMed] [CrossRef] [Google Scholar]

40. Bienaime F., Muorah M., Yammine L., Burtin M., Nguyen C., Baron W., Garbay S., Viau A., Broueilh M., Blanc T., et al. Stat3 Controls Tubulointerstitial Communication during CKD. J. Am. Soc. Nephrol. 2016; 27 :3690–3705. doi: 10.1681/ASN.2015091014. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

41. Yabusaki K., Faits T., McMullen E., Figueiredo J.L., Aikawa M., Aikawa E. A novel quantitative approach for eliminating sample-to-sample variation using a hue saturation value analysis program. PLoS ONE. 2014; 9 :e89627. doi: 10.1371/journal.pone.0089627. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

42. Chobert M.N., Couchie D., Fourcot A., Zafrani E.S., Laperche Y., Mavier P., Brouillet A. Liver precursor cells increase hepatic fibrosis induced by chronic carbon tetrachloride intoxication in rats. Lab. Investig. 2012; 92 :135–150. doi: 10.1038/labinvest.2011.143. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

43. Delire B., Henriet P., Lemoine P., Leclercq I.A., Starkel P. Chronic liver injury promotes hepatocarcinoma cell seeding and growth, associated with infiltration by macrophages. Cancer Sci. 2018; 109 :2141–2152. doi: 10.1111/cas.13628. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

44. Sun W., Chang S., Tai D.C., Tan N., Xiao G., Tang H., Yu H. Nonlinear optical microscopy: Use of second harmonic generation and two-photon microscopy for automated quantitative liver fibrosis studies. J. Biomed. Opt. 2008; 13 :064010. doi: 10.1117/1.3041159. [PubMed] [CrossRef] [Google Scholar]

45. Veidal S.S., Karsdal M.A., Nawrocki A., Larsen M.R., Dai Y., Zheng Q., Hagglund P., Vainer B., Skjot-Arkil H., Leeming D.J. Assessment of proteolytic degradation of the basement membrane: A fragment of type IV collagen as a biochemical marker for liver fibrosis. Fibrogenesis Tissue Repair. 2011; 4 :22. doi: 10.1186/1755-1536-4-22. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

46. Puchtler H., Sweat F. Histochemical Specifity of Staining Methods for Connective Tissue Fibers: Resorcin-Fuchsin and Van Gieson’s Picro-Fuchsin. Z. Zellforch. Microsk. Anat. Histochem. 1964; 79 :24–34. doi: 10.1007/BF00304175. [PubMed] [CrossRef] [Google Scholar]

47. Tashiro J., Rubio G.A., Limper A.H., Williams K., Elliot S.J., Ninou I., Aidinis V., Tzouvelekis A., Glassberg M.K. Exploring Animal Models That Resemble Idiopathic Pulmonary Fibrosis. Front. Med. 2017; 4 :118. doi: 10.3389/fmed.2017.00118. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

48. Morales-Nebreda L.I., Rogel M.R., Eisenberg J.L., Hamill K.J., Soberanes S., Nigdelioglu R., Chi M., Cho T., Radigan K.A., Ridge K.M., et al. Lung-specific loss of alpha3 laminin worsens bleomycin-induced pulmonary fibrosis. Am. J. Respir. Cell Mol. Biol. 2015; 52 :503–512. doi: 10.1165/rcmb.2014-0057OC. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

49. Colombat M., Mal H., Groussard O., Capron F., Thabut G., Jebrak G., Brugiere O., Dauriat G., Castier Y., Leseche G., et al. Pulmonary vascular lesions in end-stage idiopathic pulmonary fibrosis: Histopathologic study on lung explant specimens and correlations with pulmonary hemodynamics. Hum. Pathol. 2007; 38 :60–65. doi: 10.1016/j.humpath.2006.06.007. [PubMed] [CrossRef] [Google Scholar]

50. Parra E.R., Teodoro W.R., Velosa A.P., de Oliveira C.C., Yoshinari N.H., Capelozzi V.L. Interstitial and vascular type V collagen morphologic disorganization in usual interstitial pneumonia. J. Histochem. Cytochem. 2006; 54 :1315–1325. doi: 10.1369/jhc.6A6969.2006. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

51. Jacob J., Bartholmai B.J., Rajagopalan S., van Moorsel C.H.M., van Es H.W., van Beek F.T., Struik M.H.L., Kokosi M., Egashira R., Brun A.L., et al. Predicting Outcomes in Idiopathic Pulmonary Fibrosis Using Automated Computed Tomographic Analysis. Am. J. Respir. Crit. Care Med. 2018; 198 :767–776. doi: 10.1164/rccm.201711-2174OC. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

52. Ruscitti F., Ravanetti F., Essers J., Ridwan Y., Belenkov S., Vos W., Ferreira F., KleinJan A., van Heijningen P., Van Holsbeke C., et al. Longitudinal assessment of bleomycin-induced lung fibrosis by Micro-CT correlates with histological evaluation in mice. Multidiscip. Respir. Med. 2017; 12 :8. doi: 10.1186/s40248-017-0089-0. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

53. Zhang P., Chen J.Q., Huang W.Q., Li W., Huang Y., Zhang Z.J., Xu F.G. Renal Medulla is More Sensitive to Cisplatin than Cortex Revealed by Untargeted Mass Spectrometry-Based Metabolomics in Rats. Sci. Rep. 2017; 7 :44804. doi: 10.1038/srep44804. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

54. Tey W.K., Kuang Y.C., Ooi M.P., Khoo J.J. Automated quantification of renal interstitial fibrosis for computer-aided diagnosis: A comprehensive tissue structure segmentation method. Comput. Methods Programs Biomed. 2018; 155 :109–120. doi: 10.1016/j.cmpb.2017.12.004. [PubMed] [CrossRef] [Google Scholar]

55. Heinemann F., Birk G., Schoenberger T., Stierstorfer B. Deep neural network based histological scoring of lung fibrosis and inflammation in the mouse model system. PLoS ONE. 2018; 13 :e0202708. doi: 10.1371/journal.pone.0202708. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

56. Vogel B., Siebert H., Hofmann U., Frantz S. Determination of collagen content within picrosirius red stained paraffin-embedded tissue sections using fluorescence microscopy. MethodsX. 2015; 2 :124–134. doi: 10.1016/j.mex.2015.02.007. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

57. Drifka C.R., Loeffler A.G., Mathewson K., Mehta G., Keikhosravi A., Liu Y., Lemancik S., Ricke W.A., Weber S.M., Kao W.J., et al. Comparison of Picrosirius Red Staining With Second Harmonic Generation Imaging for the Quantification of Clinically Relevant Collagen Fiber Features in Histopathology Samples. J. Histochem. Cytochem. 2016; 64 :519–529. doi: 10.1369/0022155416659249. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

58. Wegner K.A., Keikhosravi A., Eliceiri K.W., Vezina C.M. Fluorescence of Picrosirius Red Multiplexed With Immunohistochemistry for the Quantitative Assessment of Collagen in Tissue Sections. J. Histochem. Cytochem. 2017; 65 :479–490. doi: 10.1369/0022155417718541. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

59. Vuillemin N., Mahou P., Debarre D., Gacoin T., Tharaux P.L., Schanne-Klein M.C., Supatto W., Beaurepaire E. Efficient second-harmonic imaging of collagen in histological slides using Bessel beam excitation. Sci. Rep. 2016; 6 :29863. doi: 10.1038/srep29863. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

60. Liu F., Chen L., Rao H.Y., Teng X., Ren Y.Y., Lu Y.Q., Zhang W., Wu N., Liu F.F., Wei L. Automated evaluation of liver fibrosis in thioacetamide, carbon tetrachloride, and bile duct ligation rodent models using second-harmonic generation/two-photon excited fluorescence microscopy. Lab. Investig. 2017; 97 :84–92. doi: 10.1038/labinvest.2016.128. [PubMed] [CrossRef] [Google Scholar]

61. Royce S., Erdem C., Tan A., Jaffar J., Westall G., Bourke J., Samuel C., Widdop R. Assessment of collagen in human idiopathic pulmonary fibrosis using second harmonics on Genesis 200. Eur. Respir. J. 2018; 52 doi: 10.1183/13993003.congress-2018.pa4803. [CrossRef] [Google Scholar]

62. Ranjit S., Dvornikov A., Stakic M., Hong S.H., Levi M., Evans R.M., Gratton E. Imaging Fibrosis and Separating Collagens using Second Harmonic Generation and Phasor Approach to Fluorescence Lifetime Imaging. Sci. Rep. 2015; 5 :13378. doi: 10.1038/srep13378. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

63. Qaiser T., Pugh M., Margielewska S., Hollows R., Murray P., Rajpoot N. Digital Pathology. Volume 11435. Springer; Cham, Switzerland: 2019. Digital Tumor-Collagen Proximity Signature Predicts Survival in Diffuse Large B-Cell Lymphoma; pp. 163–171. ECDP 2019. Lecture Notes in Computer Science. [CrossRef] [Google Scholar]

64. Liu Y., Keikhosravi A., Mehta G.S., Drifka C.R., Eliceiri K.W. Methods for Quantifying Fibrillar Collagen Alignment. Methods Mol. Biol. 2017; 1627 :429–451. doi: 10.1007/978-1-4939-7113-8_28. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

Articles from Biomolecules are provided here courtesy of Multidisciplinary Digital Publishing Institute (MDPI)