|J Pathol Inform 2015,
Automated morphometry provides accurate and reproducible virtual staging of liver fibrosis in chronic hepatitis C
Paul Calès1, Julien Chaigneau2, Gilles Hunault2, Sophie Michalak3, Christine Cavaro-Menard4, Jean-Baptiste Fasquel4, Sandrine Bertrais2, Marie-Christine Rousselet3
1 HIFIH Laboratory, Unité Propre de Recherche de l'Enseignement Supérieur 3859, Sructure Fédérative de Recherche 4208; Department of Liver Gastroenterology, LUNAM University, Angers, France
2 HIFIH Laboratory, Unité Propre de Recherche de l'Enseignement Supérieur 3859, Sructure Fédérative de Recherche 4208, LUNAM University, Angers, France
3 HIFIH Laboratory, Unité Propre de Recherche de l'Enseignement Supérieur 3859, Sructure Fédérative de Recherche 4208; Department of Cellular and Tissue Pathology, CHU Angers, LUNAM University, Angers, France
4 LARIS Laboratory, Unité Propre de Recherche de l'Enseignement Supérieur 7315, LUNAM University, Angers, France
|Date of Submission||01-Sep-2014|
|Date of Acceptance||16-Mar-2015|
|Date of Web Publication||28-May-2015|
HIFIH Laboratory, Unité Propre de Recherche de l'Enseignement Supérieur 3859, Sructure Fédérative de Recherche 4208; Department of Cellular and Tissue Pathology, CHU Angers, LUNAM University, Angers
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background: Liver fibrosis staging provides prognostic value, although hampered by observer variability. We used digital analysis to develop diagnostic morphometric scores for significant fibrosis, cirrhosis and fibrosis staging in chronic hepatitis C. Materials and Methods: We automated the measurement of 44 classical and new morphometric descriptors. The reference was histological METAVIR fibrosis (F) staging (F0 to F4) on liver biopsies. The derivation population included 416 patients and liver biopsies ≥20 mm-length. Two validation population included 438 patients. Results: In the derivation population, the area under the receiver operating characteristic (AUROC) for clinically significant fibrosis (F stage ≥2) of a logistic score combining 5 new descriptors (stellar fibrosis area, edge linearity, bridge thickness, bridge number, nodularity) was 0.957. The AUROC for cirrhosis of 6 new descriptors (edge linearity, nodularity, portal stellar fibrosis area, portal distance, granularity, fragmentation) was 0.994. Predicted METAVIR F staging combining 8 morphometric descriptors agreed well with METAVIR F staging by pathologists: k = 0.868. Morphometric score of clinically significant fibrosis had a higher correlation with porto-septal fibrosis area (rs = 0.835) than METAVIR F staging (rs = 0.756, P < 0.001) and the same correlations with fibrosis biomarkers, e.g., serum hyaluronate: rs = 0.484 versus rs = 0.476 for METAVIR F (P = 0.862). In the validation population, the AUROCs of clinically significant fibrosis and cirrhosis scores were, respectively: 0.893 and 0.993 in 153 patients (biopsy < 20 mm); 0.955 and 0.994 in 285 patients (biopsy ≥ 20 mm). The three morphometric diagnoses agreed with consensus expert reference as well as or better than diagnoses by first-line pathologists in 285 patients, respectively: significant fibrosis: 0.733 versus 0.733 (k), cirrhosis: 0.900 versus 0.827, METAVIR F: 0.881 versus 0.865. Conclusion: The new automated morphometric scores provide reproducible and accurate diagnoses of fibrosis stages via "virtual expert pathologist."
Keywords: Cirrhosis, digital histopathology, image analysis, liver biopsy, METAVIR staging
|How to cite this article:|
Calès P, Chaigneau J, Hunault G, Michalak S, Cavaro-Menard C, Fasquel JB, Bertrais S, Rousselet MC. Automated morphometry provides accurate and reproducible virtual staging of liver fibrosis in chronic hepatitis C. J Pathol Inform 2015;6:20
|How to cite this URL:|
Calès P, Chaigneau J, Hunault G, Michalak S, Cavaro-Menard C, Fasquel JB, Bertrais S, Rousselet MC. Automated morphometry provides accurate and reproducible virtual staging of liver fibrosis in chronic hepatitis C. J Pathol Inform [serial online] 2015 [cited 2022 Jan 25];6:20. Available from: https://www.jpathinformatics.org/text.asp?2015/6/1/20/157782
| Introduction|| |
All chronic liver diseases are characterized by the development of liver fibrosis, which includes the growth of fibrous tissue around hepatocytes (perisinusoidal fibrosis), the enlargement of portal tracts and the development of abnormal connections between portal or venular structures, named bridging or septal fibrosis. These events lead to liver architecture modifications and in turn to major complications responsible for increased mortality.  Determination of the severity of hepatic fibrosis, especially the diagnosis of bridging fibrosis and cirrhosis, is important for management. Currently, the reference method for staging liver fibrosis is the microscopic examination of a liver specimen, usually obtained by liver biopsy, by a pathologist who yields his conclusion with a semi-quantitative score, such as the METAVIR classification.  However, fibrosis staging by such a conventional liver biopsy examination is limited by poor inter/intra-observer reproducibility.  This limit can be partially circumvented by an expert reading, although a recent study showed that single experts may not reach excellent reproducibility and that an expert panel seemed to be the best reading reference.  Thus, working with centers of expertise would be interesting in liver pathology but carries with it the problem of pathologist availability. In addition, expert pathologists are lacking in some countries.
Another approach to liver evaluation is to quantify fibrosis characteristics with morphometric methods and measurements.  Area of fibrosis (also referred as collagen proportionate area) is the main descriptor, but studies have also looked at fibrosis fractal dimension  and the perimeter and size of collagenous elements.  These published techniques provide quantitative information but correlate moderately with reference pathological staging.
Our objective in the present study was to develop a computer-assisted fibrosis assessment method that correlates well with expert staging and could, therefore, circumvent observer variability or stand in for expert pathologist when they are unavailable. We first developed several new fibrosis-related morphometric descriptors that quantitatively describe the information evaluated by pathologists when they stage fibrosis. Second, we automated these measurements to provide automated scorings for the diagnosis of significant fibrosis, cirrhosis and METAVIR fibrosis stages. This proof of concept study was performed in patients with chronic viral hepatitis C.
| Materials and methods|| |
The study design is summarized in [Figure 1]. Briefly, there were four steps (the first two being chronologically independent): Construction of morphometric descriptors of fibrosis and related liver lesions by an engineer guided by a pathologist on digital images; classical fibrosis METAVIR staging independently performed by pathologists; statistical combination of morphometric descriptors allowing the statistical construction of morphometric scores targeted on fibrosis staging; and finally statistical development of morphometric diagnoses reflecting histological diagnoses.
|Figure 1: Overall study design. Construction of morphometric diagnoses in 4 steps|
Click here to view
Objectives and Definitions
The observed fibrosis stage by pathologist is called F hereafter. The main objective was to provide automated scorings for the diagnosis of significant fibrosis and cirrhosis. Significant fibrosis was considered when METAVIR stages were F ≥2. The dependent (predicted) variable was thus METAVIR stages F2 + F3 + F4 versus F0 + F1. Significant means clinically significant since long-term survival is decreased by these stages (F ≥2), and consequently anti-viral drugs are indicated. Cirrhosis was considered when METAVIR stage was = F4. The dependent variable was thus METAVIR stages F4 versus F0 + F1 + F2 + F3. The independent variables (predictors) were fine descriptors of fibrosis obtained by morphometry. The 44 descriptors measured were called morphometric descriptors. Statistical function provided logistic scores ranging from 0 to 1. This roughly means from 0% to 100% probability of the diagnostic target. The two logistic scores obtained were called morphometric scores and more precisely (morphometric) significant fibrosis score and cirrhosis score.
The automated diagnosis of METAVIR fibrosis stages was a secondary objective since the statistical functions for multinomial targets like METAVIR staging are less performant than those for binary target like significant fibrosis or cirrhosis. The discriminant function directly provided the most probable F stage that was called predicted F stage. Finally, diagnosis of significant fibrosis by significant fibrosis score, diagnosis of cirrhosis by cirrhosis score and diagnosis of Metavir F stage by predicted F stage were all grouped as morphometric diagnoses.
The two included populations had a liver biopsy scored with METAVIR fibrosis staging. The largest population with adequate (≥20 mm or cirrhosis [F4]) liver biopsy was used as the derivation population for morphometric diagnostics validated in other populations. Ultimately, 834 patients were included in the study. All patients had single-etiology, compensated, chronic viral hepatitis C without causal or anti-fibrotic treatment during the 6 months preceding the inclusion. Initial patient groups were originally recruited for studies focused on the evaluation of noninvasive liver fibrosis tests. The only additional selection criterion was the availability of a liver specimen for centralized morphometric analysis. Informed consent was obtained from patients in all studies after Institutional Review Board approval.
This population included 416 patients with a digitized liver biopsy ≥20 mm in METAVIR stages F0 to F3 or cirrhosis (liver biopsy length ≥12.1 mm) provided by a previously published population number 1  (but excluding the Fibrostar population). The F stage distribution was: 18 F0, 169 F1, 116 F2, 59 F3, and 54 F4.
First, 153 patients with liver biopsy <20 mm were available in the population number 1. This population included patients with cirrhosis explaining that the sum of the two groups in this population (n = 569) was superior to the total number of patients included (n = 549). Second, 285 other patients with liver biopsy ≥20 mm were recruited in the population number 2 from the Fibrostar study.  Thus, the two validation population included 438 patients.
All liver biopsy were performed by needle (usually suction technique, 1.6 mm diameter) intercostally. Specimens were formalin-fixed and paraffin-embedded.
The METAVIR fibrosis stages were as originally defined:  F0: No fibrosis, F1: Stellar enlargement of portal tracts without septa, F2: Enlargement of portal tracts with rare septa, F3: Numerous septa without cirrhosis, F4: Cirrhosis. In all population, METAVIR fibrosis staging was performed centrally by the same two senior liver pathologists from the Angers tertiary center with a consensus rereading in cases of initial disagreement. These pathologists (MCR, SM) were experienced in the reproducibility conditions of Metavir staging.  In the Fibrostar study, first-line diagnoses (i.e. F stages as initially assessed by the pathologists at the tertiary centers) were also available.
Digitized image analysis
Conditions for liver biopsy specimens were the same in all populations, especially as concerns length measurement and staining. 5 mm-liver sections were stained with picrosirius red without hematoxylin counterstain as the target of interest was fibrosis only. Morphometry was centrally performed by a single engineer (JC) in the Angers center. The liver specimen slides were fully scanned with an Aperio digital slide scanner (Scanscope CS System, Aperio Technologies, Vista CA, USA) image processor that provided high quality images with a maximum scanning area capacity of 120,000 × 50,000 pixels (resolution of 0.5 mm/pixel at magnification ×20). After selecting the scanning area [Figure 2]a, a digital selected image was created [Figure 2]b. Images were compressed using JPEG2000 software (quality Q = 70) in order to reduce the file size (for example, an image can have a size reduction from 1 GB to 30 MB). An automated preprocessing algorithm was used to clean the background and the tissue to analyze [Figure 2]c. When necessary, the engineer manually removed tissue artefacts or inappropriate stained structures such as folds, dust, large vascular or biliary lumens and liver capsule, which yielded a fully cleaned digital image [Figure 2]d. Further image analysis of the whole liver section included three manipulations. First, the digital color image composed of three channels red, green and blue, was decomposed to work on the green channel only as it offers the best contrast to visualize fibrosis. Only fibrosis pixels were taken into account, thanks to the preprocessing algorithm and the optional manual cleaning step that deleted artefacts or inappropriate structures. Second, this graded-intensity green image was thresholded according to an automated process: , a fuzzy generalized classification process allowed for the merging of pixel intensities into three classes (fibrosis, steatosis or healthy tissue) using the minimization of an original energy function. This produced a binary black-and-white image where fibrosis appeared in black, and all other structures appeared in white. The third step is the quantitative measure of the black space occupied by fibrosis. An example of the different steps is shown in [Figure 3]. The present innovation consisted in the distinction of different elementary lesions related to fibrosis (new morphometric descriptors) in addition to classical morphometry.
|Figure 2: Steps of the cleaning process of digital liver specimen image. (a) Glimpse of the whole slide in order to select the scanning area of a picrosirius red stained section of the whole liver specimen (green rectangle). (b) First step on the scanned liver specimen: The light-grey background and small artefacts (arrows) have to be removed by the automated process. (c) Second step: The liver specimen still displays some inappropriate areas for fibrosis measure (circles). These are, on higher magnification, c1: Liver capsule with adipose tissue; c2 and c3: Torn fragments at the edge of the specimen. (d) Final image after manual cleaning|
Click here to view
|Figure 3: Steps of the automated analysis of a digital liver specimen for the determination of the morphometric descriptors in a Metavir F3 stage. (a) Picrosirius red stained digital section after cleaning steps. (b) Green channel of the digital liver specimen used to threshold whole fibrosis. (c) Binary image after automated threshold of whole fibrosis (black pixels). (d) Automated detection of porto-septal fibrosis. (e) Automated detection of perisinusoidal fibrosis. (f) Application of the porto-septal mask to measure the granularity (granules in black pixels). (g) Automated detection of fibrous bridges. (h) Automated detection of stellar fibrosis. (i) and (j) (higher magnification). Final image with morphometric descriptors: Edge linearity percentage is highlighted in dark blue pixels, stellar fibrosis in red pixels, porto-septal fibrosis in yellow pixels, perisinusoidal fibrosis in black pixels, and bridges between portal tracts in green pixels|
Click here to view
[Table 1] summarizes the classical morphometric descriptors.
|Table 1: Patterns related to liver fibrosis measured by automated liver morphometry|
Click here to view
Area and fractal dimension of fibrosis and steatosis
Area of whole fibrosis [Figure 4]a and area of steatosis measurement has been described elsewhere.  The "box-counting" method  has been extensively used for measuring the fractal dimension of many histological objects as a complexity index. This box-counting method provided the fractal dimension of Kolmogorov: Fractal dimension of whole fibrosis or steatosis.
|Figure 4: Binary images after automated thresholding of picrosirius red-stained liver specimens. (a) Whole fibrosis in a F2 specimen. (b) Porto-septal fibrosis. (c) Perisinusoidal fibrosis. (d) Stellar fibrosis (arrows) expands a portal tract. (e) Fibrous bridge (arrow). (f) Porto-septal fibrosis in cirrhosis. (g) Application of a porto-septal mask to measure granularity. (h) Nodules detection: Circular fibrosis around nodules (1) and count of nodules (2). (i) Edge linearity: Straight (1) and irregular (2) edges|
Click here to view
Area and fractal dimension of porto-septal and perisinusoidal fibrosis
Definition assumptions: Porto-septal fibrosis was defined as fibrosis including at least one of the two following patterns: The portal fibrosis with stellar fine lobular fibrous expansions closely related to portal spaces defining F1; dense fibrous connexions (bridges or septa) between portal spaces or between portal spaces and centrilobular veins defining F2 or F3. Perisinusoidal fibrosis was defined as fibrosis around hepatocytes and around small centrilobular veins; details are provided in Supplement and in a previous study.  Area of perisinusoidal fibrosis was easily measured by the difference between the areas of whole fibrosis and of porto-septal fibrosis.
The automated mask detection of porto-septal fibrosis [Figure 4]b allowed us to distinguish perisinusoidal fibrosis [Figure 4]c from the whole fibrosis. We thus measured the area and fractal dimension of porto-septal fibrosis, the area and fractal dimension of perisinusoidal fibrosis, the area of lobular perisinusoidal fibrosis and the ratio of perisinusoidal fibrosis area among the whole fibrosis; details on fibrosis metrics are available elsewhere.  The automated mask detection of porto-septal fibrosis and the new descriptors described were measured at magnification ×5 since they do not need high-resolution images (details in the Supplement [Additional file 1]).
New morphometric descriptors
[Table 1] summarizes the new morphometric descriptors; technical details on the measurement of the following lesions and characteristics are provided in the Supplement.
The main feature differentiating F0 from higher stages is the presence of stellar fibrosis [Figure 4]d, that is, small, fine fibrils branching from the dense porto-septal regions. Therefore, we automated the quantitation of this stellar fibrosis on the whole liver biopsy specimen surface (whole area of stellar fibrosis) on the surface of porto-septal regions (portal area of stellar fibrosis) and on the surface of lobular regions (lobular area of stellar fibrosis). We also measured the mean area of stellar fibrosis per porto-septal region (mean portal area of stellar fibrosis).
We measured the number of porto-septal regions; the mean area of porto-septal regions (mean area of porto-septal fibrosis) and the average distance between the portal tracts (mean portal distance). This last parameter reflects either a normal regular distribution of portal tracts or an architectural distortion with some portal tracts closer to others. In advanced fibrosis, portal tracts embedded in septa are automatically detected as a unique porto-septal region and consequently the mean portal distance is calculated as zero.
Bridging (or septal) fibrosis [Figure 4]e occurs in stages F2 to F4. A bridge was defined with morphometry as a structure lying between two thick elements. There is no literature reference to distinguish F2 from F3 stages with strict quantitative criteria; pathologists subjectively choose between a few septa (F2) and numerous septa (F3). The process of automated scoring required a quantitative definition. Guided by the personal experience of the participating pathologists, we considered that the F3 or F4 stage has been reached when the ratio of bridges was >50%. Therefore, we measured the number of bridges, the ratio of bridges in the porto-septal regions (portal ratio of bridges), the area of fibrosis in the bridges (area of bridging fibrosis) and the mean bridge thickness obtained from their perimeter (mean bridge perimeter) and surface (mean bridge area). Liver architecture is modified in high F stages by many fibrous bridges [Figure 4]f. These bridges isolate parenchymal areas among lobules that we called granules [Figure 4]g. Thus, we measured the mean granularity percentage quantifying the liver structure disruption by bridging fibrosis.
Cirrhosis is mainly characterized by the formation of nodules, which we defined in image analysis as circular and non-fibrotic (lacking fibrous septa) areas surrounded by fibrosis [Figure 4]h1. We measured the number of nodules [Figure 4]h2 and the mean percentage of fibrosis around nodules (mean nodularity percentage).
Liver biopsy specimens may be fragmented. This is fibrosis stage dependent and occurs particularly in cirrhosis where there may be several small (≤5 mm 2 ) fragments. We measured the number of fragments and the fragmentation index; this latter defined as the ratio between the surface of detected small fragments and the total liver biopsy specimen surface.
Liver biopsy specimens change in shape according to F stages: Edges are rather straight in low stages (F0, F1, F2) [Figure 4]i1 but tend to become more curved and irregular in high stages (F3, F4) [Figure 4]i2]. Thus, we automated the measurement of the fractal dimension of edges and the edge linearity percentage.
Liver specimen dimensions
Pathologists usually consider that a representative liver biopsy specimen must be ≥15 or ≥20 mm long. We automated the measurement of liver specimen length, which required specific rules especially for twisted specimens. We also measured liver specimen perimeter and liver specimen surface.
The performance of morphometric descriptors depends on specimen staining quality (picrosirius red in the present study). We measured stain intensity on the three color components of the image (red, green, blue) and obtained a mean color intensity index by averaging all pixel intensities for, respectively, the whole specimen, fibrosis and parenchyma. The term "color intensity" was used to refer to this evaluation of staining quality. The contrast between the mean color intensity of fibrosis and either the liver biopsy specimen or the parenchyma was also recorded.
Noninvasive Fibrosis Tests
The two following blood tests were calculated: Fibrotest,  FibroMeter V2G .  Liver stiffness  evaluation by elastometry (Fibroscan, Echosens, Paris, France) and the combined test Elasto-FibroMeter 2G were also performed.
Quantitative variables were expressed as mean ± standard deviation. The diagnostic accuracy of each morphometric score was expressed as the area under the receiver operating characteristic (AUROC) and the overall accuracy (rate of well-classified patients according to F stages). Data were reported according to STARD statements  and analyzed on an intention to diagnose the basis.
Independent predictors of binary diagnostic targets (significant fibrosis and cirrhosis) were determined by binary logistic regression. Independent predictors of F stages were determined by a linear discriminant analysis, providing a probability score from 0 to 1 for each F stage. Independent predictors of liver stiffness were determined by multiple linear regression.
The main statistical analyses were performed by professional statisticians (GH, SB) using SPSS version 18.0 (IBM, Armonk, NY, USA), SAS 9.1 (SAS Institute Inc., Cary, NC, USA) and R software. 
| Results|| |
The main characteristics of the population are depicted in [Table 2].
Development of Morphometric Scores
Significant fibrosis diagnosis
In the derivation population, the best logistic model provided a morphometric significant fibrosis score combining 5 independent morphometric descriptors: Whole area of stellar fibrosis (P < 0.0001 in multivariate analysis), edge linearity percentage (P < 0.0001), mean bridge thickness (P = 0.0018), number of bridges (P = 0.0032) and mean nodularity percentage (P = 0.0293); descriptor rank is that of stepwise analysis as indicated in [Table 1]. Examples of score determination are shown in [Figure 5]. The value of morphometric score, varying from 0 to 1, is proportional to METAVIR F stage [Figure 6], e.g. predicted F1 values are (median, 1 st and 3 rd quartiles): 0.011 (0.003-0.042), predicted F2: 0.102 (0.038-0.269). With an a priori threshold score at 0.5, the AUROC was 0.957 (95% confidence interval: 0.940-0.973) and 87.3% (84.1-90.5) of the patients were correctly classified for significant fibrosis diagnosis. [Table 3] provides the overall performance results.
|Figure 5: Examples of digital sections of Metavir stages: F1 (a), F2 (b), F3 (c), and F4 (d). First, liver specimens stained with picrosirius red (panels 1) then complete automated morphometric analysis (panels 2 and 3). They show the final images with the five morphometric descriptors of the morphometric significant fibrosis score: (a2) edge linearity percentage in dark blue pixels, stellar fibrosis in red pixels, porto-septal fibrosis in yellow pixels, perisinusoidal fibrosis in black pixels; (b2) one fibrous bridge between portal tracts is highlighted in green pixels; (c2) porto-septal fibrosis with many bridges in green pixels; (d2) nodules surrounded by >80% of fibrosis are highlighted in cyan pixels. Higher magnification of morphometric details of the four Metavir stages (a3, b3, c3 and d3)|
Click here to view
|Figure 6: Relationship between morphometric scores and METAVIR F staging. Population number 1 (549 patients). Box plots (median, interquartile range and extremes) of morphometric scores against METAVIR F staging|
Click here to view
|Table 3: Diagnostic accuracy for significant fibrosis and cirrhosis by respective morphometric scores|
Click here to view
The best logistic model provided a morphometric cirrhosis score including 6 independent morphometric descriptors: Edge linearity percentage (P = 0.0058), mean nodularity percentage (P = 0.0029), portal area of stellar fibrosis (P < 0.0001), mean portal distance (P = 0.0145), mean granularity percentage (P = 0.0003) and fragmentation index (P = 0.0012). The morphometric score threshold was selected here a posteriori at 0.1567 to maximize sensitivity with minimal deterioration of specificity. Thus, sensitivity was 100% and ensuing specificity was 96.4%. The AUROC was 0.994 (0.989-0.999), and 96.9% (95.2-98.5) of the patients were correctly classified for cirrhosis diagnosis. [Table 3] provides the overall performance results.
METAVIR F stage diagnosis
The final discriminant model included 8 descriptors of predicted F stage: Mean granularity percentage (P < 0.0001), edge linearity percentage (P < 0.0001), portal ratio of bridges (P < 0.0001), fractal dimension of perisinusoidal fibrosis (P < 0.0001), ratio of perisinusoidal fibrosis area (P < 0.0001), fractal dimension of porto-septal fibrosis (P < 0.0001), fragmentation index (P = 0.0001) and mean nodularity percentage (P = 0.0016). Agreement between observed METAVIR F (by pathologists) and predicted F stage (by automated morphometry) was very good according to the weighted kappa index: 0.868 (0.844-0.891). Thus, there was no significant discordance (difference ≥2 F stages) between observed METAVIR F and predicted F stage. The rate of correctly classified patients for F stage diagnosis was, overall: 68.5% (64.0-72.3); and per METAVIR stage: F0: 40.0%, F1: 74.7%, F2: 61.7%, F3: 55.9%, F4: 82.0%.
Validation of Morphometric Descriptors and Scores
The metrical definition of the new morphometric descriptors was validated by a good correlation with METAVIR F, quantitative classical measures of fibrosis (whole and porto-septal fibrosis area) and noninvasive liver fibrosis tests, especially those quite independent of liver pathology (liver stiffness, hyaluronate, alpha2-macroglobulin) [details in supplement, [Table S1]. In addition, the metrical definition of porto-septal fibrosis area was validated by a higher correlation (P < 0.0001) with METAVIR F stages (rs : 0.756) than did whole area of fibrosis (rs : 0.598) [Table 4].
|Table 4: Correlations (Spearman's rho) between fibrosis descriptors. Derivation population number 1 (416 patients)|
Click here to view
The 13 independent morphometric descriptors for the 3 previous diagnostic targets showed a significant and regular progression (increase or decrease) as a function of METAVIR F stages [Figure 7] except for portal area of stellar fibrosis [Table 5].
|Figure 7: Relationship between main composite descriptors (with comparable % units) of morphometric scores and METAVIR F staging. Population number 1 (549 patients). Box plots (median, interquartile range and extremes) of composite descriptors against METAVIR F staging|
Click here to view
|Table 5: Values (mean±SD) and correlation of the 13 independent morphometric descriptors as a function of METAVIR F stages. Derivation population number 1 (416 patients)|
Click here to view
The metrical definition of other morphometric descriptors was indirectly validated by the very significant increase in diagnostic accuracy of new morphometric scores versus classical morphometric descriptors [Table 6].
|Table 6: Comparison of diagnostic accuracy (AUROC) of new morphometric scores and classical morphometric descriptors for the two main diagnostic targets. Derivation population number1 (416 patients)|
Click here to view
Fibrosis morphometric scores were validated by high correlations with other liver fibrosis descriptors in the derivation population, especially the histological references: Area of porto-septal fibrosis (reference independent from their construction) and METAVIR F stage [Table 4]. Morphometric significant fibrosis score was the score the most correlated with area of porto-septal fibrosis (rs : 0.835) and METAVIR F (rs : 0.851).
Comparison with METAVIR F
Significant fibrosis score had a higher correlation with area of porto-septal fibrosis (rs : 0.835) than METAVIR F (rs : 0.756, P < 0.001). Significant fibrosis score allowed morphometry to reach the same high correlations with noninvasive fibrosis tests than did METAVIR F, e.g. rs with serum hyaluronate: 0.484 versus 0.476 (P = 0.862), respectively, [Table S1] whereas the correlations of other morphometric descriptors were significantly different from those of METAVIR F.
The plots of morphometric significant fibrosis and cirrhosis scores versus reference F stages showed a marked progression thus an excellent F stage discrimination by morphometric scores [Figure 6].
The scatter plot of significant fibrosis score versus cirrhosis score showed a perfect complementarity with a progression of cirrhosis score from 0 to 1 only when significant fibrosis score reached 1 [Figure 8].
|Figure 8: Scatter plot of significant fibrosis score versus cirrhosis score as a function of METAVIR F stages. Bold line indicates the nonlinear regression by LOWESS. Derivation population number 1 (416 patients)|
Click here to view
Accuracies for the 3 morphometric diagnoses were validated in one different population (number 2). AUROCs of morphometric significant fibrosis and cirrhosis scores were, respectively: 0.955 and 0.994. [Table 3] presents the overall results for significant fibrosis and cirrhosis diagnosis. Considering the predicted F stage, there was globally an expected (optimism bias deletion), but small, decrease in accuracy. Thus, the overall rate of correctly classified patients was 66.3% for F stage diagnosis.
We analyzed the effect of liver biopsy length on diagnostic performance. AUROCs of morphometric significant fibrosis and cirrhosis scores were, respectively: 0.893 and 0.993 in 153 patients with biopsy <20 mm; the decrease in accuracy when liver biopsy length was <20 mm was thus minor [Table 3].
Reproducibility and cleaning effect
Intra-investigator reproducibility of morphometric scores and section cleaning were separately evaluated in 20 patients (details in Supplement). Intra-investigator reproducibility was excellent with intra-class correlation coefficient ≥0.993 for the three scores [Table S2]. Cleaning had no statistically significant effect on morphometric scores. However, at the individual level there was a clinically significant impact in about 1 out 20 patients (details in Supplement).
The clinical interest (as well as metrical definition) of morphometric descriptors was evaluated by taking liver stiffness as an independent reference in the derivation population. By multivariate analysis including METAVIR F and activity, classical morphometry and the 13 main new independent descriptors [Table 1], liver stiffness was independently predicted (overall R²: 0.502) by porto-septal fibrosis area (R²: 0.349), METAVIR F ( R²: 0.069), fractal dimension of porto-septal fibrosis ( R²: 0.029), area of steatosis ( R²: 0.012), mean granularity percentage ( R²: 0.006), fragmentation index ( R²: 0.007), mean portal distance ( R²: 0.006), fractal dimension of perisinusoidal fibrosis ( R²: 0.007), portal area of stellar fibrosis ( R²: 0.012) and whole area of stellar fibrosis ( R²: 0.005). Thus, 86.3% of the information were brought by morphometry (among which 14.4% by new descriptors) and 13.7% by METAVIR F.
In the validation population, the agreement between histological F reference (expert consensus) and the three morphometric diagnoses was equal or superior to the agreement between histological reference (expert consensus) and the histological diagnosis by first-line pathologists from tertiary centers, respectively, morphometry versus first-line pathologist: Significant fibrosis: 0.733 versus 0.733 (k), cirrhosis: 0.900 versus 0.827 (k) and METAVIR F staging: 0.881 versus 0.865 (weighted k).
| Discussion|| |
To date, morphometry has been used mainly to measure the surface (area) and subsequently the fractal dimension of liver lesions. Initially, morphometry was based on the manual determination of regions of interest  before some procedures became automated.  In the present study, we used a highly automated technique. Dioguardi et al. standardized metrical evaluation of the geometric properties of the parenchyma, inflammation, fibrosis, and alterations in liver tissue tectonics, but its application was not evaluated. 
The present work describes a method for automating the pathological staging of liver fibrosis in chronic viral hepatitis C based on a histopathological classification frequently used by pathologists. Among the 44 morphometric descriptors, 13 had independent value for the three diagnostic targets (significant fibrosis, cirrhosis, F stage). Most were expected as they reflect patterns linked to METAVIR F staging; however, the most accurate and constant was a new morphometric descriptor, the edge linearity percentage, which underlines the importance of considering indirect signs of liver fibrosis, such as the external deformation of the liver biopsy specimen. Other descriptors were unexpected like those reflecting perisinusoidal fibrosis.
We used consensus expert diagnosis as a reference to construct morphometric scores for the binary diagnosis. The morphometric scores of significant fibrosis and cirrhosis only included new morphometric descriptors, which again emphasizes the originality of the present study. Moreover, the morphometric diagnosis of F stages was based on a different statistical method. This was the sole diagnosis to include several relatively classical descriptors, such as area or fractal dimension of fibrosis.
It should be noted that the binary morphometric diagnosis of significant fibrosis and cirrhosis was based on a numerical score from 0 to 1. These morphometric scores can be used as quantitative descriptors of the fibrosis degree intrinsically. Thus, morphometric significant fibrosis score had the highest correlations with other fibrosis descriptors, especially excellent correlations with area of porto-septal fibrosis or METAVIR F stages. With their precision and reproducibility, morphometric scores might be more sensitive in detecting changes than semi-quantitative scorings, from a statistical point of view. Thus, it has recently been suggested that morphometry of liver fibrosis, with a simple measurement such as area of fibrosis,  could be the best prognosis predictor. Moreover, liver morphometry might be a sensitive descriptor for predicting or monitoring liver fibrosis progression. , The present new morphometric scores and classical morphometry will need to be compared in these domains. Morphometric significant fibrosis and cirrhosis scores were very complementary [Figure 8]. This complementarity has recently been observed with corresponding blood tests for significant fibrosis and cirrhosis.  Thus, cirrhosis score can well describe the fibrosis progression within cirrhosis stage whereas significant fibrosis score is the basic score. Indeed, significant fibrosis score was the most accurate score with a significantly higher correlation with porto-septal fibrosis area than METAVIR F staging and the only morphometric diagnosis reaching high correlations with noninvasive fibrosis tests, not significantly different from METAVIR F staging correlations.
To summarize, based on expert consensus in a large, well-defined database, we created algorithms providing automated morphometric diagnoses reflecting that made by expert pathologists. In this sense, this new technique may be considered as a "virtual expert pathologist", providing reproducible and accurate automated diagnoses.
Although we used a consensus expert reading, the main limit to our study is the reference based on classical histological fibrosis staging, which comprises two sources of variability: The observer and the specimen size/location. Specimen size has been evaluated in the present work, but location bias will be evaluated in further study. Liver biopsy should be considered as a best standard but not a gold standard  or an imperfect gold standard.  It has already been shown that blood tests, when constructed on liver biopsy, might offer better prognostication than their reference for construction.  Comparatively, it remains to be determined if morphometric scores provide better prognostication than histological staging.
Our material included liver biopsy with a pathological diagnosis of cirrhosis irrespective of the liver specimen length. Recently, sample size requirements for the digital analysis of cirrhosis was evaluated and appeared dependent of liver disease etiology.  Furthermore, no decompensated cirrhosis was included in our study. Morphometry of decompensated cirrhosis might improve cirrhosis prognostication  and should be compared to a prognostic histological score, e.g. the Laennec score,  several items of which were included in the present morphometric scores.
The main factors limiting the reliability of image analysis are thresholding, artefacts, staining quality and the representativeness of whole liver fibrosis. A classical limiting factor is thresholding for black-and-white binarization but, in our technique, thresholding was automated. , We observed that cleaning of artefacts had nonstatistically significant effect on morphometric scores, but they induced a large individual variation in about 1 out 20 patients. Since this technical prerequisite is useful, a further study aimed at the automated detection of artefacts is ongoing. The lack of standardization in staining techniques may be another limit. Color intensity parameters may reduce variability, but this needs further testing. In the present study, none of the color intensity parameters influenced diagnosis. This putative limit might be circumvented by nonstaining techniques.  Staining performed in other centers will be evaluated, and optionally our algorithms will be adapted. Finally, our results show that the present new morphometric scores are less sensitive to specimen length than area of fibrosis. 
Another limit is the interest of developing morphology compared to that of noninvasive tests, especially in chronic hepatitis C. First, it should be underlined that a weakness of classical histological fibrosis assessment is observer variability, furthermore present in all areas of diagnosis in medicine like hepatology,  especially in clinical practice.  However, this inconvenience was only marginally present in the present automated morphometric diagnosis. Thus, the application of reliable and consistent imaging techniques might renew interest in liver biopsy utilization, particularly when issues of fibrosis and staging are paramount. Moreover, virtual liver biopsy , may become available in the near future, for which automated morphometry would be a natural partner.
Images of the area of portal or perisinusoidal fibrosis were quantitatively validated by experts.  Images of the new morphometric descriptors were validated qualitatively by experts and quantitatively by correlations with other fibrosis descriptors and accuracies for pathological diagnoses or liver stiffness. Other validation steps are needed: Other reference pathologists, relationship with other staging systems like Ishak's staging, other population and clinical outcomes.
| Conclusion|| |
The morphometric descriptors presented here describe the full range of main lesions present in chronic viral hepatitis. The morphometric scores provide accurate qualitative diagnosis of significant fibrosis, cirrhosis and METAVIR fibrosis classification. They can be also used as quantitative descriptors useful for precise staging and follow-up by experienced physicians (e.g., pathologists, hepatologists) for diagnosis, prognosis and drug trials. The results of the present study open a door toward "virtual expert pathology." However, as fibrosis staging is only a part of liver biopsy evaluation, classical histopathology remains a cornerstone in liver diseases.
| Acknowledgments|| |
We thank all participants from the SNIFF 12 study:  Jérτme Boursier, Frédéric Oberti, Isabelle Hubert Fouchard; and from the Fibrostar study supported by ANRS:  Jean Pierre Zarski was the principal investigator with Nathalie Sturm as the principal pathologist. Kevin Erwin provided writing assistance for English proofreading.
| References|| |
Poynard T, Bedossa P, Opolon P. Natural history of liver fibrosis progression in patients with chronic hepatitis C. The OBSVIRC, METAVIR, CLINIVIR, and DOSVIRC groups. Lancet 1997;349:825-32.
Intraobserver and interobserver variations in liver biopsy interpretation in patients with chronic hepatitis C. The French METAVIR Cooperative Study Group. Hepatology 1994;20:15-20.
Rousselet MC, Michalak S, Dupré F, Croué A, Bedossa P, Saint-André JP, et al.
Sources of variability in histological scoring of chronic viral hepatitis. Hepatology 2005;41:257-64.
Boursier J, Bertrais S, Oberti F, Gallois Y, Fouchard-Hubert I, Rousselet MC, et al.
Comparison of accuracy of fibrosis degree classifications by liver biopsy and non-invasive tests in chronic hepatitis C. BMC Gastroenterol 2011;11:132.
Germani G, Burroughs AK, Dhillon AP. The relationship between liver disease stage and liver fibrosis: A tangled web. Histopathology 2010;57:773-84.
Dioguardi N, Grizzi F, Bossi P, Roncalli M. Fractal and spectral dimension analysis of liver fibrosis in needle biopsy specimens. Anal Quant Cytol Histol 1999;21:262-6.
Dioguardi N, Grizzi F, Fiamengo B, Russo C. Metrically measuring liver biopsy: A chronic hepatitis B and C computer-aided morphologic description. World J Gastroenterol 2008;14:7335-44.
Boursier J, de Ledinghen V, Zarski JP, Fouchard-Hubert I, Gallois Y, Oberti F, et al.
Comparison of eight diagnostic algorithms for liver fibrosis in hepatitis C: New algorithms are more precise and entirely noninvasive. Hepatology 2012;55:58-67.
Zarski JP, Sturm N, Guechot J, Paris A, Zafrani ES, Asselah T, et al.
Comparison of nine blood tests and transient elastography for liver fibrosis in chronic hepatitis C: The ANRS HCEP-23 study. J Hepatol 2012;56:55-62.
Roullier V, Cavaro-Ménard C, Guillaume C, Aubé C. Fuzzy algorithms to extract vacuoles of steatosis on liver histological color images. Conf Proc IEEE Eng Med Biol Soc 2007;2007:5575-8.
Sandrini J, Boursier J, Chaigneau J, Sturm N, Zarski JP, Le Bail B, et al.
Quantification of portal-bridging fibrosis area more accurately reflects fibrosis stage and liver stiffness than whole fibrosis or perisinusoidal fibrosis areas in chronic hepatitis C. Mod Pathol 2014;27:1035-45.
Boursier J, Chaigneau J, Roullier V, Lainé F, Sandrini J, Michalak S, et al.
Steatosis degree, measured by morphometry, is linked to other liver lesions and metabolic syndrome components in patients with NAFLD. Eur J Gastroenterol Hepatol 2011;23:974-81.
Moal F, Chappard D, Wang J, Vuillemin E, Michalak-Provost S, Rousselet MC, et al.
Fractal dimension can distinguish models and pharmacologic changes in liver fibrosis in rats. Hepatology 2002;36:840-9.
Castéra L, Vergniol J, Foucher J, Le Bail B, Chanteloup E, Haaser M, et al.
Prospective comparison of transient elastography, Fibrotest, APRI, and liver biopsy for the assessment of fibrosis in chronic hepatitis C. Gastroenterology 2005;128:343-50.
Leroy V, Halfon P, Bacq Y, Boursier J, Rousselet MC, Bourlière M, et al.
Diagnostic accuracy, reproducibility and robustness of fibrosis blood tests in chronic hepatitis C: A meta-analysis with individual data. Clin Biochem 2008;41:1368-76.
Calès P, Boursier J, Ducancelle A, Oberti F, Hubert I, Hunault G, et al.
Improved fibrosis staging by elastometry and blood test in chronic hepatitis C. Liver Int 2014;34:907-17.
Bossuyt PM, Reitsma JB, Bruns DE, Gatsonis CA, Glasziou PP, Irwig LM, et al
. The STARD statement for reporting studies of diagnostic accuracy: Explanation and elaboration. Clin Chem 2003;49:7-18.
Team RDC. R: A Language and Environment for Statistical Computing; 2011. Available from: http://www.R-project.org/.
Pilette C, Rousselet MC, Bedossa P, Chappard D, Oberti F, Rifflet H, et al.
Histopathological evaluation of liver fibrosis: Quantitative image analysis vs semi-quantitative scores. Comparison with serum markers. J Hepatol 1998;28:439-46.
Calvaruso V, Dhillon AP, Tsochatzis E, Manousou P, Grillo F, Germani G, et al.
Liver collagen proportionate area predicts decompensation in patients with recurrent hepatitis C virus cirrhosis after liver transplantation. J Gastroenterol Hepatol 2012;27:1227-32.
Calès P, Zarski JP, Chapplain JM, Bertrais S, Sturm N, Michelet C, et al.
Fibrosis progression under maintenance interferon in hepatitis C is better detected by blood test than liver morphometry. J Viral Hepat 2012;19:e143-53.
Manousou P, Burroughs AK, Tsochatzis E, Isgro G, Hall A, Green A, et al.
Digital image analysis of collagen assessment of progression of fibrosis in recurrent HCV after liver transplantation. J Hepatol 2013;58:962-8.
Boursier J, Brochard C, Bertrais S, Michalak S, Gallois Y, Fouchard-Hubert I, et al.
Combination of blood tests for significant fibrosis and cirrhosis improves the assessment of liver-prognosis in chronic hepatitis C. Aliment Pharmacol Ther 2014;40:178-88.
Bedossa P, Carrat F. Liver biopsy: The best, not the gold standard. J Hepatol 2009;50:1-3.
Poynard T, Ingiliz P, Elkrief L, Munteanu M, Lebray P, Morra R, et al.
Concordance in a world without a gold standard: A new non-invasive methodology for improving accuracy of fibrosis markers. PLoS One 2008;3:e3857.
Ngo Y, Munteanu M, Messous D, Charlotte F, Imbert-Bismut F, Thabut D, et al.
A prospective analysis of the prognostic value of biomarkers (FibroTest) in patients with chronic hepatitis C. Clin Chem 2006;52:1887-96.
Hall AR, Tsochatzis E, Morris R, Burroughs AK, Dhillon AP. Sample size requirement for digital image analysis of collagen proportionate area in cirrhotic livers. Histopathology 2013;62:421-30.
Kim SU, Oh HJ, Wanless IR, Lee S, Han KH, Park YN. The Laennec staging system for histological sub-classification of cirrhosis is useful for stratification of prognosis in patients with liver cirrhosis. J Hepatol 2012;57:556-63.
Gailhouste L, Le Grand Y, Odin C, Guyader D, Turlin B, Ezan F, et al.
Fibrillar collagen scoring by second harmonic microscopy: A new tool in the assessment of liver fibrosis. J Hepatol 2010;52:398-406.
Bedossa P, Dargère D, Paradis V. Sampling variability of liver fibrosis in chronic hepatitis C. Hepatology 2003;38:1449-57.
Winkfield B, Aubé C, Burtin P, Calès P. Inter-observer and intra-observer variability in hepatology. Eur J Gastroenterol Hepatol 2003;15:959-66.
Campo-Ruiz V, Lauwers GY, Anderson RR, Delgado-Baeza E, González S. In vivo
and ex vivo
virtual biopsy of the liver with near-infrared, reflectance confocal microscopy. Mod Pathol 2005;18:290-300.
Fabila D, de la Rosa JM, Stolik S, Moreno E, Suárez-Álvarez K, López-Navarrete G, et al. In vivo
assessment of liver fibrosis using diffuse reflectance and fluorescence spectroscopy: A proof of concept. Photodiagnosis Photodyn Ther 2012;9:376-82.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8]
[Table 1], [Table 2], [Table 3], [Table 4], [Table 5], [Table 6]
|This article has been cited by|
||Autophagy of hepatic stellate cell induced by Clonorchis sinensis
| ||Bao Zheng, Zhiyan Gao, Liumei Liang, Yunyu Lu, Yongting Kong, Wanting Chen, Keying Lin, Wanqi Chen, Jingying Mai, Yanwen Li, Changling Ma |
| ||Molecular Biology Reports. 2021; |
|[Pubmed] | [DOI]|
||Liver fibrosis staging by computed tomography: prospective randomized multicentric evaluation of image analyses
| ||Carole Vitellius, Anita Paisant, Adrien Lannes, Julien Chaigneau, Frédéric Oberti, Jérôme Lebigot, Isabelle Fouchard, Jérôme Boursier, Pascal David, Christophe Aubé, Paul Calès, the CDMIR group |
| ||Clinics and Research in Hepatology and Gastroenterology. 2021; : 101797 |
|[Pubmed] | [DOI]|
||Fractal analysis tools for early assessment of liver inflammation induced by chronic consumption of linseed, palm and sunflower oils
| ||Dejan Opric, Anica D. Stankovich, Aleksandra Nenadovic, Sanjin Kovacevic, Danilo D. Obradovic, Silvio de Luka, Jelena Neovic-Ostojic, Jelena Milain, Andjelija . Ilic, Alexander M. Trbovich |
| ||Biomedical Signal Processing and Control. 2020; 61: 101959 |
|[Pubmed] | [DOI]|
||Elastogram quality assessment score in vibration-controlled transient elastography: Diagnostic performance compared to digital morphometric analysis of liver biopsy in chronic hepatitis C
| ||L. C. Mendes,P. A. Ferreira,N. Miotto,L. Zanaga,E. S. L. Gonçales,M. N. Pedro,M. S. Lazarini,F. L. G. Júnior,R. S. B. Stucchi,A. G. Vigani |
| ||Journal of Viral Hepatitis. 2017; |
|[Pubmed] | [DOI]|
||Evaluating segmental liver function using T1 mapping on Gd-EOB-DTPA-enhanced MRI with a 3.0 Tesla
| ||Zhi-Peng Zhou,Li-Ling Long,Wei-Jia Qiu,Ge Cheng,Li-Juan Huang,Teng-Fei Yang,Zhong-Kui Huang |
| ||BMC Medical Imaging. 2017; 17(1) |
|[Pubmed] | [DOI]|
||Repeatability, Reproducibility, Separative Power and Subjectivity of Different Fish Morphometric Analysis Methods
| ||Péter Takács,Zoltán Vitál,Árpád Ferincz,Ádám Staszny,Nicholas S Foulkes |
| ||PLOS ONE. 2016; 11(6): e0157890 |
|[Pubmed] | [DOI]|
||Management of chronic hepatitis C in 2017
| ||L. von Köckritz,J.-F. Dufour |
| ||Hämostaseologie. 2016; 37(1) |
|[Pubmed] | [DOI]|
||A high throughput approach for analysis of cell nuclear deformability at single cell level
| ||Menekse Ermis,Derya Akkaynak,Pu Chen,Utkan Demirci,Vasif Hasirci |
| ||Scientific Reports. 2016; 6: 36917 |
|[Pubmed] | [DOI]|