Introduction
Fibrosis is described as the biological evidence of an error in the dynamic process of cellular repair1. Liver cirrhosis is the final pathological condition resulting from persistent damage in the liver, in which healthy tissue is replaced by scar tissue, which causes a decompensation and alteration in its functioning and architecture, often irreversible. Cirrhosis represents an important cause of morbidity and mortality in the world population2. In some countries of Latin America, it occupies between the fifth and sixth place as a cause of general death. In Mexico, the mortality due to liver cirrhosis in 2013 was positioned as cause number five3. Although this condition is considered a problem that implies the loss of global health, there is not enough, accurate, and updated information on the epidemiological scope and, moreover, on the precise clinical characteristics of its evolution and reversion processes, thus limiting the possibilities of more adequate prevention, intervention, and assistance4.
To contribute to the study of this pathology, multiple trials have been conducted using animal models. One of the methods frequently found in the literature, due to its simplicity and relative reproducibility, is the use of substances that generate a significant activation of the hepatic stellar cells, which causes an increase in the production of extracellular matrix, inflammation, and necrosis5. Some of the components used for these models are, for example, dimethylnitrosamine, galactosamine, and carbon tetrachloride (CCl4). An important limitation in the study of liver cirrhosis is the subjective evaluation of the progression of fibrosis, which in experimental trials is commonly estimated from the quantification of serum markers and anatomopathological studies of the biopsy6-8. In the latter case, a pathology specialist performs qualitative assessments of the structural characteristics of the tissue. The inter- and intra-observer variability of the diagnoses, once the biopsy is taken, is also criticized and represents an important factor of variability9. Although there are different classification systems for liver fibrosis, one of the most frequently used is the Metavir score10.
Several works have been aimed to the development of digital image analysis methods for the quantification of hepatic fibrosis11-15. Dahab et al.11. developed an automated method for the digital quantification of fibrosis in liver images of rats making correlation studies on a qualitative conventional scale. The automation used is a restriction, since the experience of a skilled observer is disregarded and the evaluations depend mostly on the quality and staining of the histological specimen. En Goodman et al.12 presented a study to evaluate an anti-fibrosis agent using a commercial automated system for the acquisition of biopsy images, with a processing in color RGB and grayscales using cross-polarized filters. As a consequence of tissue fragmentation and insufficient length and/or width, 144 cases of the 389 patients had one or both biopsies discarded for analysis. They reported a correlation coefficient of 0.49 between phases of fibrosis in the Ishak score. Although results indicate a degree of average correlation, they showed the importance of applying this type of techniques to be used in future therapeutic trials of agents to inhibit the progression of fibrosis. Zaitoun et al.13 present the quantification of periportal, pericellular, and septal fibrosis; as well as steatosis, by stereological and morphometric techniques, in biopsies of patients with chronic hepatitis C and patients with alcoholic liver disease (ALD). They found a relationship between the quantification of fibrosis by means of the commercial system Profit 5.2™, and the modified histological score of Ishak. However, in patients with cirrhosis, the amount of fibrosis varied between 20% and 74%. Significant differences were found in the amount of fibrosis between patients with chronic hepatitis C and those with ALD.
In this work, we present a methodology for the quantification of fibrosis in liver images based on morphometric processing of the color represented in the LAB color space of the Commission International de l’Eclairage. In the LAB space, it is possible to represent more colors with respect to the RGB scale; in addition to, better emulate the non-linear response of the human eye. The color representation is given in a three-dimensional space of the components of luminosity (L), transition of red-green (A), and transition of yellow-cyan (B). In this space, the color of the fibrosis selected in the image is represented as a point in the three-dimensional space L*A*B. Thus, the final obtaining of the color of the fibrosis is calculated from the selection of three zones identified as fibrosis, where each of the LAB components is averaged to find a resulting color, which is used as the selection value of each pixel of the image in which color matches. The percentage of fibrosis is calculated by counting the pixels that match the selected value with respect to the total number of pixels in the image. We present the correlation between the fibrosis quantification in a murine model and the histopathological evaluation, hence validating the principle of the proposed methodology.
Materials and methods
We conducted a descriptive, longitudinal, prospective, experimental, analytical, and correlational study.
Animal model of induction of fibrosis by CCl4
Fibrosis was induced in male Wistar rats of 6 weeks of age, weighing between 190 and 220 g. Four groups of four rats were randomly formed, in which CCl4 was injected intraperitoneally in olive oil, in an initial dose of 250 µl, at a concentration of 33%, twice a week16,17. After individuals gained 20% of their initial weight, the induction dose was increased weekly to 1 µl/g. At the end of the induction time (Table 1), individuals were sacrificed with an overdose of sodium pentobarbital (Pisabental, Pisa agropecuaria®) intraperitoneally and a total hepatectomy performed.
Liver samples and imaging
A cut was systematically taken from the right hepatic medial lobe of all the hepatectomized livers and stained with Masson trichrome. Using an optical microscope, the slides were photographed at a ×10, to obtain ten images of each one, controlling that the field of vision completely corresponded to tissue.
Histopathological evaluation of the images in Metavir score
Three histopathological evaluations were carried out in a blinded and independent manner of each image using the Metavir score, described in table 2. The final evaluation was taken as the most concordant, thus resulting in evaluations with complete agreement, namely 3/3; or evaluations with one non-concordant evaluator, namely, two-third. When the disparity was complete (none of the three evaluators agreed on the evaluation), the opinion of a fourth was requested to perform the conclusive evaluation based on 2/4 coincidences.
Morphometric quantification of fibrosis
The quantification of fibrosis was performed in the 160 images taken from the liver tissue. The quantification was performed by the same experienced user, trained in the visual identification of fibrosis and its histological characteristics. The quantification was performed by selecting fibrous tissue from three different areas of the image as distant as possible, to be more representative of possible color differences in the staining.
Statistical analysis
The concordance between the evaluators for each staging was calculated. The data obtained from the quantification of the fibrosis percentage (FP) were expressed as mean and standard deviation (SD) and to determine statistical differences between the induction groups, an ANOVA test was applied; the significance was determined when p < 0.05. A Spearman’s rho test was used to establish the correlation between the FP variable and the histopathological evaluation. The statistical analysis was carried out with the program SPSS™, v.15.
Results
As previously described, a histological cut stained with Masson trichrome was obtained from each of the 16 individuals and photographed in ten different areas, thus acquiring a total of 160 images.
Histopathological evaluation
The range of distribution of fibrosis evaluation in each induction group (T1-T4) is shown in figure 1. The Metavir score evaluation for the 40 images for each induction group overlaps with other fibrosis phases that did not correspond to the expected phase according to the induction time. The concordance in the histopathological evaluation of fibrosis is illustrated in figure 2A. The frequency of the evaluations of the 160 analyzed images is shown in figure 2B according to the Metavir evaluation groups being: 46 in F4, 33 in F3, 40 in F2, 41 in F1, and no cases in F0.
Quantification of fibrosis
Using the morphometric algorithm, the FP in the 160 images was calculated. The FP was grouped according to the evaluated Metavir phase, the box plot is shown in figure 3. Table 3 shows the descriptive statistics for the FP for each group according to the Metavir evaluation. The results of the ANOVA test show significant statistical differences for the Metavir evaluation groups (p = 0.0001). To specify which groups differ from each other, a post hoc analysis was carried out and the statistical differences were found to be significant in all cases (p = 0.0001). On the other hand, Spearman’s rho test indicates a strong correlation between the percentage of fibrosis and histopathological diagnosis (p = 0.0001, r = 0.9).
Discussion
Among the alternatives of digital analysis for the images applied in the evaluation of hepatic fibrosis, it was found that some of the methods have important limitations, such as poor accuracy to detect perisinusoidal fibrosis and/or considerable coefficients of variation to evaluate cirrhosis. In most of the works, the image processing was done in the RGB color space; this feature mainly restricts the range of colors available. In addition, a full automation for the acquisition and processing of the images compromises the evaluation of defects and artifacts in tissue specimens. Many of the software were designed as a general purpose tool, and therefore, not specific for this pathology yielding important errors derived from the particularities in the hepatic architecture.
In this work, we present a morphometric method for the quantification of hepatic fibrosis, emphasizing its applicability as an auxiliary tool in the evaluation and diagnosis of this pathology. As shown in figure 1, most of the evaluated images in the induction Groups T1, T2, and T4 were in the expected phases F1, F2, and F4, respectively. The dispersion of the liver fibrosis data suggested that the induction times did not generate a predominant stage of fibrosis but rather a dispersion of stages that are more evident in the T3 group in which 47.5% of the images evaluated corresponded to F4 and only 30% to F3 which was the expected phase. For the 160 images analyzed, there was a concordance of the three evaluators (3/3) in 70%. However, 48 complex cases were also found, in which at least one of the observers did not coincide with the rest. It is also observed that for F1 the staging with a 3/3 agreement was achieved in 58.5% of the cases, for F2 in 67.5%, and for F3 in 66.7%, while F4 had the highest number of cases totally concordant with 84.8%. This analysis suggests that the inter-observer variability is greater in cases with incipient (F1) or moderate fibrosis (F2 and F3). The fibrosis evaluation was distributed as follows: the most populated group was the group with F4 evaluation, with 28.8% of the images, followed by F1 and F2 with 25.6% and 25%, respectively. Group F3 is the least numerous representing 20.6% of cases. Despite having groups with a non-homogeneous population, it can be observed that the difference between the most populated and the least populated barely represents 8% of the total cases.
The results of the quantification of the FP range from 0.21% to 7.75%. In figure 3, it is possible to notice an increase in FP as the fibrosis phase is more advanced with a non-linear trend. A strong correlation was obtained between these two variables. Numerically, the average of the FP for each Metavir stage (Table 3) has an increase of 1.34% when going from F1 to F2, 1.98% in F2 to F3, and 0.85% in the transition from F3 to F4. Using this measurement and the SD for each group, FP ranges were calculated for each Metavir phase; this is illustrated in the diagram of figure 4. The groups are distinguished from each other except for the overlap between Groups F3 and F4; however, the ANOVA test and its post hoc indicate that all the groups are statistically different from each other.
The method of induction of fibrosis by CCl4 was strict and carefully monitored to minimize errors; however, the animal model revealed a heterogeneous distribution of the disorder, within the same individual and between individuals of the same group corresponding to the same induction group: T1, T2, T3y, and T4. This is evidenced in the individual-by-individual analysis for each group, as shown in table 4, where, at least one individual was evaluated in a different Metavir stage than that expected, (this is to say), T1 was expected to correspond to F1 and so on to T4/F4. Especially, an ethereal behavior is observed in the individual B of the Group T1, which evolved to F3 phase, even when it was submitted only to 4 weeks of induction. The T3 group was the most seriously affected since only one out of four individuals was evaluated in expected phase, F3.
Induction group | Individual | Metavir | FP | ||
---|---|---|---|---|---|
Evaluation | Range | (%) | (%) | ||
T1 | A | F1 | 1 | 0.8 | 0.5 |
B | F3 | 2 | 2.8 | 0.8 | |
C | F1 | 1 | 1.4 | 0.3 | |
D | F1 | 2 | 0.9 | 0.9 | |
T2 | A | F2 | 1 | 1.5 | 0.6 |
B | F2 | 1 | 1.5 | 0.8 | |
C | F1 | 1 | 0.7 | 0.3 | |
D | F2 | 2 | 2.7 | 1.0 | |
T3 | A | F2 | 2 | 2.7 | 1.1 |
B | F4 | 1 | 3.8 | 0.7 | |
C | F3 | 1 | 4.5 | 1.0 | |
D | F4 | 1 | 4.6 | 1.1 | |
T4 | A | F4 | 0 | 5.6 | 1.5 |
B | F4 | 0 | 5.1 | 0.9 | |
C | F3 | 1 | 4.3 | 1.2 | |
D | F4 | 1 | 5.2 | 1.0 |
FP: fibrosis percentage.
The dispersion of fibrosis within the same individual can be seen through the Metavir and FP evaluation. For the first case, 12.5% of the individuals showed a rank equal to 0, that is, only in two of the 16 rats, all the observations were evaluated in the same Metavir phase. On the other hand, 25% of the individuals of the model have a rank in Metavir evaluations equal to two that means that in four rats, up to three different phases were found in the same individual. Regarding the FP, it was observed that the SD, of the quantification per individual goes from 0.3% to 1.5%. The average SD of the individuals staged in F1 is equal to 0.5%, 0.86% for F2, 1% for F3, and 1.0% for F4. The evaluation of these two methodologies suggests that the use of the Metavir score, yield more robust diagnoses for F4, since observations are based on shape characteristics, namely, the observation of complete nodules, regardless of the amount of fibrosis, or the thickness of the septa that form those nodules. This same discussion could justify the dispersion of FP that was observed in Group T4 and that is greater (> 0.9 for all individuals) than in the other induction groups.
A qualitative evaluation such as anatomo-histopathological assessment does not allow complete precision in the appraisal of liver damage once liver cirrhosis has been diagnosed. A more precise quantification of fibrosis would facilitate a distinction in the staging of liver cirrhosis, where a quantitative measurement is fundamental to evaluate the effect of treatments for the control or reversion of the pathology. Therefore, the quantification of FP by digital morphometry can increase the resolution for a complete evaluation of liver fibrosis, especially in intermediate stages.
Conclusions
In conclusion, this study confirms the use of HepaScan as a quantitative tool to evaluate the amount of fibrosis in liver biopsy images in a murine model of hepatic fibrosis. The morphometric quantification in terms of the fibrosis percentage FP had a high correlation with the histopathological evaluation of images of histological sections of liver biopsies. The quantification of the FP by digital morphometry will potentiate its use in the evaluation of the effect of treatments for the control or reversion of the pathology.