Hostname: page-component-cd9895bd7-mkpzs Total loading time: 0 Render date: 2024-12-22T10:06:44.816Z Has data issue: false hasContentIssue false

Fibrosis severity scoring on Sirius red histology with multiple-instance deep learning

Published online by Cambridge University Press:  18 July 2023

Sneha N. Naik*
Affiliation:
ITMAT Data Science Group, NIHR Imperial BRC, Imperial College, London, United Kingdom Heffner Biomedical Imaging Lab, Department of Biomedical Engineering, Columbia University, New York, NY, USA
Roberta Forlano
Affiliation:
Faculty of Medicine, Department of Metabolism, Digestion and Reproduction, Imperial College, London, United Kingdom
Pinelopi Manousou
Affiliation:
Department of Hepatology, Imperial College Healthcare NHS Trust, London, United Kingdom
Robert Goldin
Affiliation:
Section for Pathology, Imperial College, London, United Kingdom
Elsa D. Angelini
Affiliation:
ITMAT Data Science Group, NIHR Imperial BRC, Imperial College, London, United Kingdom Heffner Biomedical Imaging Lab, Department of Biomedical Engineering, Columbia University, New York, NY, USA Faculty of Medicine, Department of Metabolism, Digestion and Reproduction, Imperial College, London, United Kingdom Telecom Paris, Institut Polytechnique de Paris, LTCI, Palaiseau, France
*
Corresponding author: Sneha N. Naik; Email: [email protected]

Abstract

Non-alcoholic fatty liver disease (NAFLD) is now the leading cause of chronic liver disease, affecting approximately 30% of people worldwide. Histopathology reading of fibrosis patterns is crucial to diagnosing NAFLD. In particular, separating mild from severe stages corresponds to a critical transition as it correlates with clinical outcomes. Deep Learning for digitized histopathology whole-slide images (WSIs) can reduce high inter- and intra-rater variability. We demonstrate a novel solution to score fibrosis severity on a retrospective cohort of 152 Sirius-Red WSIs, with fibrosis stage annotated at slide level by an expert pathologist. We exploit multiple instance learning and multiple-inferences to address the sparsity of pathological signs. We achieved an accuracy of $ 78.98\pm 5.86\% $, an F1 score of $ 77.99\pm 5.64\%, $ and an AUC of $ 0.87\pm 0.06 $. These results set new state-of-the-art benchmarks for this application.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2023. Published by Cambridge University Press

Impact Statement

We demonstrate, for the first time, the capabilities of an end-to-end weakly supervised multiple-instance DL pipeline to accurately classify liver fibrosis stages between mild (stage $ \le 2 $ ) and severe (stage $ >2 $ ) on whole-slide images. Such a framework works towards reducing expert annotation workload and human-induced uncertainty.

1. Introduction

1.1. NAFLD and fibrosis scoring

Non-alcoholic fatty liver disease (NAFLD) is currently the leading cause of chronic liver disease worldwide and is typically associated with obesity, insulin resistance, and diabetes(Reference Godoy-Matos, Júnior and Valerio 1 ). In a study of over 8 million people worldwide, it was estimated that the global prevalence of NAFLD is 30% with the highest prevalence in the Middle East and South America(Reference Younossi, Koenig, Abdelatif, Fazel, Henry and Wymer 2 ). Furthermore, in the USA, NAFLD became the first indication for liver transplantation in 2022(Reference Lonardo, Mantovani, Petta, Carraro, Byrne and Targher 3 ). Given the increasing worldwide incidence of obesity, it is likely that NAFLD in the global population will continue to rise in decades to come. Fibrosis stage is the main predictor of liver- and non-liver-related outcomes and fibrosis progression or non-worsening of fibrosis is the end-point in all clinical trials. Therefore, the evaluation of fibrosis progression in clinical trials is of key interest.

1.2. Histopathology of the liver

Liver biopsy is the reference “gold” standard for the diagnosis and staging of NAFLD(Reference Boyd, Cain, Chauhan and Webb 4 ). The liver consists of the portal tracts (containing bile ducts and vascular branches), terminal hepatic venule, hepatocytes, and sinusoids. Histopathologists most commonly use the semi-quantitative Non-alcoholic Steatohepatitis Clinical Research Network (NASH CRN) scoring system to stage the severity of liver disease seen on biopsies(Reference Brunt, Janney, Di Bisceglie, Neuschwander-Tetri and Bacon 5 ). The scoring system (described in detail in(Reference Kleiner, Brunt and Van Natta 6 , Reference Nascimbeni, Ballestri and Machado 7 )) relies on visual signs of four pathological processes(Reference Boyd, Cain, Chauhan and Webb 4 ): Steatosis, Lobular inflammation, Ballooning, and Fibrosis. The presence of the first three components (Steatosis, Ballooning, Lobular inflammation) are combined into an unweighted sum to generate the NAFLD Activity Score (NAS score) which ranges from 0 to 8 and is an indicator of active injury in the liver. Fibrosis is defined as the accumulation of extracellular matrix proteins (e.g., collagen) in the liver(Reference Bataller and Brenner 8 ). It is scored separately from the other features as it is generally considered a result of disease activity rather than a feature of active injury and is less reversible. Fibrosis is graded from 0 to 4 based on the following architectural patterns visible: Stage 0: No visible fibrosis; Stage 1a: Mild perisinusoidal fibrosis; Stage 1b: Moderate-to-severe perisinusoidal fibrosis; Stage 1c: Portal / periportal fibrosis; Stage 2: perisinusoidal and portal / periportal fibrosis; Stage 3: Bridging fibrosis between adjacent portal tracts or between a portal tract and an adjacent hepatic venule(Reference Boyd, Cain, Chauhan and Webb 4 ); Stage 4: Cirrhosis, which can lead to cell necrosis and liver failure.

Traditional glass microscopy slides are scanned using a digital slide scanner to generate whole-slide images (WSIs). Advances in slide scanning technology and reduction in the cost of digital storage have accelerated the design of solutions for automated segmentation of disease patterns on WSIs(Reference Hu, Li and Li 9 ), and most recently the scoring of those patterns.

1.3. Deep learning in histopathology

Deep Learning (DL) classifiers are raising significant breakthroughs in WSI segmentation, scoring, and classification(Reference Hu, Li and Li 9 Reference Zhou, Li and Rahaman 11 ). The vast majority of published work focuses on H&E stains for cancer datasets.

WSIs are very large images, typically of size 0.5–1 Gb per WSI. This prevents opening WSIs at full resolution on standard computers, but also loading the whole field of view at full resolution on standard GPUs to train DL architectures. Options commonly used in DL are either working with pixel-level annotations–which is time-consuming to prepare or using multiple instance learning (MIL) with annotations at the WSI level. In a MIL framework, small “tiles” are extracted from WSIs to form “bags” with a single label assigned to all bags from the same WSI and which corresponds to the WSI label (commonly called weak labeling). Features are extracted from each tile (typically using a CNN), and tile embeddings or per-tile predictions are aggregated to generate an overall prediction from the bag(Reference Ilse, Tomczak and Welling 12 ).

MIL frameworks rely on specific strategies for tile selection, feature extraction, and multiple inference aggregation. Common MIL paradigms rely on one training “bag” per image after careful tile selection (e.g., tissue content above a threshold)(Reference Ilse, Tomczak and Welling 12 Reference Li, Li, Gertych, Knudsen, Speier and Arnold 14 ); data augmentations(Reference Ilse, Tomczak and Welling 12 Reference Xie, Muhammad and Vanderbilt 17 ); use of transfer learning from ImageNet via GoogLeNet, InceptionNet, ResNet, and MobileNet(Reference Bejnordi, Veta and Van Diest 18 , Reference Kwok 19 ) for cancer, and AlexNet and ResNet for NAFLD(Reference Arjmand, Christou and Tzallas 20 , Reference Qu, Minacapelli and Tait 21 ); aggregation of tile-level inferences via max-pooling, or of tile-level features via average-pooling, attention-based or RNN based frameworks ahead of WSI-level prediction(Reference Ilse, Tomczak and Welling 12 , Reference Campanella, Hanna and Geneslaw 13 , Reference Lu, Williamson, Chen, Chen, Barbieri and Mahmood 16 , Reference Qu, Minacapelli and Tait 21 , Reference Bilal, Jewsbury and Wang 22 ).

Stain normalization is often used as a pre-processing step. Classic approaches include template color transfer, such as Reinhard stain normalization(Reference Reinhard, Adhikhmin, Gooch and Shirley 23 ), and stain-deconvolution(Reference Ruifrok and Johnston 24 ). More recently, DL architectures were proposed such as sparse auto-encoders(Reference Janowczyk, Basavanhally and Madabhushi 25 ), GANs(Reference Cho, Lim, Choi and Min 26 ), and the recent Stain Standardization Capsule (SSC) architecture(Reference Zheng, Jiang and Zhang 27 ) which outperformed both color augmentation and state-of-the-art normalization methods on three tumor classification tasks on H&E stained datasets.

1.4. Specific challenges

A specific challenge faced by computational methods when scoring liver fibrosis on those Sirius-Red stained WSIs is the presence of natural collagen around the portal tracts and hepatic veins. This causes direct measures of the percentage area of collagen to be similar for different fibrosis stages(Reference Forlano, Mullish and Giannakeas 28 ), and suggests that both locating and quantifying the abnormal fibrosis patterns is crucial.

Histological scoring systems used to manually assess liver biopsy, such as the NASH CRN, use a semi-quantitative approach, based on describing the localization and visual patterns of liver fibrosis, but without measuring the amount of fibrosis. Such a scoring system is very subjective and is subject to significant inter and intra-observer variability. Inter-reader weighted kappa for fibrosis staging in NAFLD has been reported approximately 0.484 from scores of three independent hepatopathologists on 687 H&E-stained biopsies of NASH patients. Intra-reader kappas in the same study were moderately higher at 0.679(Reference Davison, Harrison and Cotter 29 ). This is considered as modest to poor agreement (on a score of −1 for total disagreement to 1 for perfect agreement). Reported kappas vary across datasets: some show fair inter-rater agreement at 0.72–0.79 for H&E and Trichrome datasets with sizes ranging from 40 to 300+ WSIs(Reference Soon, Liu, Leow, Wee, Wei and Sanyal 30 , Reference Taylor-Weiner, Pokkalla and Han 31 ). Others report cases of a two-stage difference in fibrosis severity assigned by two independent pathologists(Reference Gawrieh, Sethunath and Cummings 32 ).

Such a scoring system is also insufficient, indeed, the absence of a quantitative assessment has important consequences in clinical practice as the extent of liver fibrosis is commonly associated with poor clinical outcomes in these patients. Overall, the lack of an objective and quantitative assessment of histology has an enormous impact on the interpretation of clinical trials, as it might not be sensitive enough to detect small changes in fibrosis pattern. Therefore, in order to better diagnose and monitor liver fibrosis status in clinical trials on NAFLD, there is a strong need for fast and reproducible scoring of WSIs(Reference Nam, Chapiro, Paradis, Seraphin and Kather 33 ).

In this paper, we propose a MIL framework to train a weakly supervised binary fibrosis classifier to distinguish Stages 0–2 from Stages 3–4, considered as a clinically-crucial fibrosis severity transition.

2. Methods

Our proposed method was designed to tackle a binary classification task with the following challenges at hand: small image cohort, only image-level annotations, sparse pathological signs in large tissue samples, significant stain color variability between samples and presence of artifacts inherent to non-curated retrospective clinical data.

2.1. The LiFib cohort

We exploit a retrospective clinical cohort, referred to in this paper as the LiFib cohort, of 152 digitized WSIs from liver biopsy samples performed as part of routine clinical care at St Mary’s Hospital. Liver tissue was collected via needle biopsy using a percutaneous technique under ultrasound guidance. WSIs were recorded with a Hamamatsu NanoZoomer2.0 HT digital slide scanner at 40x magnification using Sirius-Red stain. Each WSI contains two to three cores of tissue. The WSI images in our cohort have a pixel size of $ 0.23\mu m $ at max zoom factor of 40x. The WSIs were annotated at a slide level by an expert pathologist with 35 years of clinical experience. The four-stage fibrosis scale mentioned in the introduction is used, with the following distribution in our cohort: stage 0 = 12%, stage 1 = 23%, stage 2 = 15%, stage 3 = 35%, and stage 4 = 15%. For the binary classification task of separating mild (stage $ \le 2) $ from severe (stage $ >2 $ ) fibrosis, we have an equal number (n = 76) of cases per class.

Our WSIs are affected by common slide-preparation artifacts: air bubbles, knife-slicing irregularity, tissue folds, cracks, crumbling tissue (particularly for severe fibrosis cases), and fixation problems. WSIs also contain red artifacts where a vein was cut with the tissue during extraction (Figure 3a). From visual inspection, blur digitization artifacts affect only a small portion of WSIs. The LiFib dataset was not curated to exclude cases with artifacts, as those are common in routine clinical data.

We manually characterized WSI visual appearance into three categories (Figure 2): (a) green tissue on a grey background (n = 75), (b) green tissue on a blue background (n = 56), (c) yellow tissue on a grey background (n = 21). In all cases, collagen and fibrosis components appear as shades of red. However, in older WSIs, the stain is more faded.

We constructed five independent test sets for five-fold cross-validation, stratifying on fibrosis scores and color categories. For each fold a 70%–10%–20% train-validation-test split was used.

2.2. Bag construction and deep-learning pipelines

An overview of the proposed bag construction and DL pipelines is provided in Figure 1.

Figure 1. Overview of our bag construction and deep-learning pipelines for fibrosis scoring on Sirius Red histopathology WSIs. Two MIL bag construction methods are compared: (1) RED-1, which uses a priori knowledge on red appearance of fibrosis to generate 1 bag per WSI, and (2) RAND-n, which generates $ n $ bags per WSI only based on tissue content. Deep-learning pipeline (best model): SE-ResNet18 initialized with ImageNet pre-trained weights with gated attention aggregation to output a binary label on a given bag.

2.2.1. Tissue segmentation and tile extraction

Tissue segmentation was performed using Ilastik’s Pixel Classification workflow(Reference Berg, Kutra and Kroeger 34 ). Background from foreground separation was trained with manual brush strokes on one WSI per color category, which led to satisfactory segmentation quality on all cases. Tiles were generated on a grid using 2.5 mm square patches at 5× magnification (pixel size 1651 × 1651), with 50% overlap. We excluded tiles from the background using a threshold on tissue content. This threshold was set to 50% except in few cases impacted by tissue crumbling for which it was iteratively decreased to retain a minimum of 10 tiles per WSI.

2.2.2. Bag composition

Bags are composed of a pre-fixed number of $ k $ tiles and are given the original WSI label as a proxy ground-truth label. The majority of existing literature constructs a single bag per WSI that includes from thousands to tens of thousands of tiles with deemed-sufficient tissue content(Reference Ilse, Tomczak and Welling 12 Reference Lu, Williamson, Chen, Chen, Barbieri and Mahmood 16 ). Unlike most cancer focused datasets, the pathological signs in NAFLD are sparse. In the LiFib cohort, for cases of mild-medium fibrosis, a large proportion of tiles generated are uninformative (containing only artifact or an absence of collagen). Consequently, using too many tiles per bag might result in a prohibitively weak learning signal in each bag. We therefore tested two alternative approaches: (1) Dataloader RED-1: we generate 1 bag per WSI using the $ k $ tiles with the highest percentage of red pixel content (proxy for fibrosis tissue) measured before stain normalization. (2) Dataloader RAND-n: we generate $ n=\lceil \frac{m}{k}\rceil $ bags per WSI with $ k $ tiles randomly sampled without replacement from a pool of $ m $ available valid tile candidates. This leads to a variable number of bags per WSI. Within each bag, tiles are downsampled to 224 × 224 pixels using bicubic interpolations and data augmentation was performed using random rotations in range $ \left[-50,+50\right]{}^{\circ} $ , horizontal and vertical flips.

2.2.3. MIL architectures and multiple inference aggregation

Our baseline architecture uses a ResNet18 pre-trained on ImageNet. Our advanced architecture uses the Squeeze-and-Excite(Reference Hu, Shen and Sun 35 ) module in a pre-trained SE-ResNet18 encoder to enhance robustness of features to artefacts and upweight the importance of pathological signs.

For the aggregation of tile-level features in the MIL inference, we tested max-pooling and a gated attention framework(Reference Ilse, Tomczak and Welling 12 ), shown to outperform max-pooling of tile-level predictions on small datasets(Reference Lu, Williamson, Chen, Chen, Barbieri and Mahmood 16 ).

The max-pooling architecture applies a single linear layer to the encoder outputs $ {\mathbf{x}}_k $ , to generate per-tile predictions $ {y}_k=\mathrm{sigm}\left({\mathbf{Wx}}_k^T\right) $ and assigns to the WSI a label equal to the maximum of per-tile predictions, $ Y=\underset{k}{\max }{y}_k $ . Gated Attention aggregation (Figure 1, Equation 1) first computes attention weights per-tile $ \left({a}_k\right) $ using two linear layers $ \left(\mathbf{U},\mathbf{V}\right) $ . A single WSI-level feature vector $ \left(\mathbf{z}\right) $ is generated from the weighted sum of tile-level features. Finally, a single linear layer $ \left(\mathbf{W}\right) $ returns the final overall WSI label.

(1) $$ {\displaystyle \begin{array}{l}{a}_k=\frac{\exp \left\{{\mathbf{w}}^T\left(\tanh \left({\mathbf{Vx}}_k^T\right)\odot \mathrm{sigm}\left({\mathbf{Ux}}_k^T\right)\right)\right\}}{\sum_{j=1}^K\exp \left\{{\mathbf{w}}^T\left(\tanh \left({\mathbf{Vx}}_j^T\right)\odot \mathrm{sigm}\left({\mathbf{Ux}}_j^T\right)\right)\right\}}\\ {}\mathbf{z}=\sum \limits_{k=1}^n{a}_k{\mathbf{x}}_k\\ {}Y=\mathrm{sigm}\left({\mathbf{Wz}}^T\right)\end{array}} $$

Model training and evaluation procedures are different for our two bag construction methods. (1) Dataloader RED- 1:Using 1 bag per WSI, assigned with the WSI label as proxy ground truth. During training and testing, one inference pass is run to generate the bag-level prediction. (2) Dataloader RAND-n: Using n bags per WSI, all assigned with the WSI label as proxy ground truths. We conduct one inference pass per bag to generate n bag-level predictions. We train using the binary cross-entropy loss between each bag prediction, and the proxy ground truth label. During training and testing, one inference pass is run per bag to generate n bag-level predictions. The final WSI prediction is set to the most severe predicted fibrosis stage across the bag-level predictions (one-wins-all aggregation).

3. Results

In our cohort, we extracted on average $ m=37\pm 19 $ valid tiles per WSI, with $ m\in \left[10\;100\right] $ (12 WSIs with only $ m=10 $ tiles). We tested bag sizes $ k\in \left[5\;20\right] $ for Dataloader RED-1 and we worked with $ n $ = 1 to 10 bags for Dataloader RAND-n for bag size $ k=10 $ .

We report mean and standard deviation values over the 5 folds for Accuracy and F1 scores, to assess models performance and reproducibility across training sets. We also report Recall, Specificity, Precision, F2, and Matthews Correlation Coefficient (MCC) values measured on the test confusion matrix aggregated over the 5 folds.

3.1. Single bag approach: MIL aggregation and bag size effect

We first report results in Table 1 using Dataloader RED-1 (cf. Section 2.2) considered as the most favorable training setup guided with a priori knowledge on disease sign characteristics. The three most important hyperparameters to reach training stability across folds were the encoder learning rate (= $ {10}^{-4} $ ), batch size (=6), and early stopping on the validation set. Other optimally tuned hyperparameters include weight decay (= $ {10}^{-4} $ ) and Adam optimizer.

Table 1. Single-bag MIL architectures performance–pre-trained with ImageNet and using Dataloader RED-1 with bag size $ k=10 $ .

Note. Two MIL aggregation strategies are compared: max-pooling on tile inference or gated attention on tile features. Best values are indicated in bold.

From Table 1, the SE-ResNet18 with gated attention shows the best performance. We tested bag size effects on this architecture in Table 2 with $ k= $ 5 to 20.

Table 2. Effect of bag sizes $ k $ using SE-ResNet18 (Gated Attention) pre-trained with ImageNet on Dataloader RED-1 with $ n=1 $ bag per WSI.

Note. Other hyperparameters as in Table 1. Best values are indicated in bold.

We can see in Table 2 that bag size can decrease Recall by 20% corresponding to an increase in false negatives. A bag size $ k $ = 10 leads to the best metrics except for Specificity (−4% compared to $ k $ = 15) and was therefore selected for subsequent experiments. These results confirm the tradeoff that exists in MIL between bags too small that can lead to missed information (e.g., here in cases with lots of red artifacts, where thresholding on red pixels did not retain fibrosis signs of interest) and bags too large that can degrade learning capacity.

At this stage, our model suffers from two limitations likely due to some information loss when using Dataloader RED-1: (a) low Recall (64.47%) and (b) large performance variability across folds (10% standard deviation in F1). The subsequent refinement aims to improve these results.

3.2. Gain from multiple inference training

We compare in Table 3 our baseline that uses Dataloader RED-1 with the use of multiple bags with Dataloader RAND-n during testing only or also during training. We use the optimal bag size for Dataloader RED-1 (k = 10). The number of bags $ n $ varies across WSIs, depending on the number of valid tiles with an average of $ m=37\pm 19 $ tiles extracted per WSI and a range of valid tiles per WSI from 10 to 100, hence $ n $ ranging from 1 to 10 for $ k=10 $ . Class balance was maintained when using RAND-n, working with a total of 313 bags for the mild fibrosis class and 318 bags for the severe fibrosis class.

Table 3. Comparing Dataloader RED-1 ( $ n=1 $ ) versus Dataloader RAND-n with bag size $ k=10. $

Note. RAND-n uses multiple inferences with “one wins all” aggregation. All other details as in Table 2. Best values are indicated in bold.

Using Dataloader RAND-n only at testing has a detrimental effect on Specificity and Precision (−7.9% and − 5.7%) which corresponds to an increase in false positives. On the other hand, using Dataloader RAND-n at both training and testing leads to the best performing model (except for −1.3% in Specificity), with clear gains in Accuracy of ( $ +4.6\% $ ), and F1 score ( $ +7.8\% $ ) along with stable or decreased standard deviations across folds. It also returns an increase of +9.5% in Recall, corresponding to fewer false negatives, which could be attributed to the prevention of information missed at inference by utilizing more tiles in a WSI. This model achieves an AUC of $ 0.87\pm 0.06 $ across the five test folds.

3.3. Qualitative evaluation of attention coefficients

Our best model (SE-ResNet18 + Gated Attention + Dataloader RAND-n train/test) outputs attention weights for all tiles inside bags, which can be interpreted as the importance given by the network to a given tile to make the final inference at the bag level. We reviewed tiles with high and low attention weights with three clinical experts. For positive cases, we focused separately on true-positive (TP–Figure 2) and false-positive (FP–Figure 3) cases. We also reviewed whole WSIs for false-negative (FN–Figure 4) cases to analyze the missed pathological signs.

Figure 2. True positive (TP) cases of severe fibrosis illustrated on tiles with the 10 highest attention weights. WSI color categories are indicated under each case. Arrows color-coding: blue = relevant pathological signs of bridging, purple = red artefacts, green = healthy portal tracts. (a–f) Cases with high attention correctly focused on fibrotic signs as well as on artifacts and portal tracts. (g) Case showing higher attention put on portal tracts sliced longitudinally than on fibrotic signs. (h) Case with strong blurring still showing high attention on fibrotic signs.

Figure 3. False positive (FP) cases of severe fibrosis. (a) Red artifacts: The tiles with two highest and two lowest attention weights show a focus on red pixels from a vein cut within the tissue. (b) Borderline case: This case was reconsidered to be a potential severe case by the expert clinician.

Figure 4. False Negative (FN) cases. (a) Faded red stain; (b) Poor-quality biopsy with tissue crumbling; (c,d) Borderline cases with sparse bridging signs shown in zoomed boxes: (c) most red pixels correspond to healthy portal tracts sliced longitudinally; (d) Very thin bridging patterns.

This extensive review leads to the following:

General observations:

  1. 1. Tiles with high attention scores in TP cases are usually of high diagnostic significance with the confirmed presence of fibrotic patterns (cf. Figure 2).

  2. 2. The presence of large red artifacts in TP cases, such as vein tissue in cirrhotic cases biopsied via wires (Figure 2c) did not distract the model to pay high attention to fibrotic patterns.

  3. 3. Some common characteristics seem to explain the majority of cases classified incorrectly by our model, (cf. Figures 3 and 4), detailed in the next paragraph.

Going now into more details on errors patterns:

  1. 4. False-positive WSIs generally fall into two categories: (a) misleading large red artifacts such as a vein cut lengthwise (Figure 3a) shown in all tiles with high attention (i.e., focus on a single large region in the WSI); (b) borderline cases reclassified as a severe case after reviewing the highest attention weighted tiles (Figure 3b).

  2. 5. False-negative WSIs have four key characteristics: (a) faded stain and loss of color contrast on old WSIs in color categories green/grey and green/blue (Figure 4a) which would require higher magnification detect fibrosis; (b) poor-quality biopsy due to tissue crumbling in advanced fibrosis cases (Figure 4b); (c) cases with very delicate fibrosis or very sparse pathological signs (Figure 4c,d). In such cases, high attention was correctly assigned to tiles with healthy tissue. Case in Figure 4c corresponds to a WSI that is overall healthy except for a single, borderline bridge in the zoomed box. It is interesting to note on this case that red structures from the lumen of blood vessels cut longitudinally at branching points were not misinterpreted as pathological signs. Case in Figure 4d corresponds to a WSI where fibrosis bridging is only a few cells thick; There was also a case where a large piece of portal tract was surrounded by faint fibrosis; (d) WSIs in color category yellow/grey are more often misclassified as healthy tissue, likely due to the small sample size.

4. Conclusions

In this study, we have proposed an end-to-end weakly-supervised multiple-instance DL framework that can distinguish severe from mild fibrosis on histopathological WSIs from a retrospective clinical dataset stained with Sirius Red. To the best of our knowledge, this is the first study applying DL on Sirius-Red stained WSIs to stage fibrosis.

We designed our method to handle multiple challenges in this new application: small image cohort, sparse signs on large tissue samples, large stain variability and presence of visual artifacts, a priori knowledge on the red color of pathological signs of interest (e.g., bridging).

Our key contributions are as follows: We first demonstrated the benefit of using two advanced concepts (squeeze-excitation encoder and gated attention pooling) when training on one bag per WSI composed of tiles with high red-pixel content. These first experiments lead to an average baseline Accuracy of $ 74.32\pm 5.38\% $ and an average F1 score of $ 70.18\pm 10\% $ .

We further demonstrated the benefit of switching to multiple inferences on multiple bags per WSI for our problem. Our final model has an average accuracy of $ 78.98\pm 5.86\% $ , and an average F1 score of $ 77.99\pm 5.64\% $ and an AUC of $ 0.87\pm 0.06 $ on the LiFib dataset. Finally, we showed some clinical interest in reviewing tiles with high attention weights, which lead to one case reclassified as severe.

Our results compare well with the performance reported on fibrosis staging using a different stain and more annotations(Reference Gawrieh, Sethunath and Cummings 32 ). Our results are also on par with an NAFLD scoring solution trained with WSI-level annotations on H&E and Trichrome stained biopsies(Reference Qu, Minacapelli and Tait 21 ).

We acknowledge three limitations in our work. First, we relied on a single expert for ground truth annotations, while visual fibrosis stage scoring on WSIs has high inter- and intra- clinical expert variability. Second, GPU constraints limited batch and bag size ranges that could be tested. Finally, we conducted minimal pre-processing on our clinical dataset.

Based on an extensive visual review of our results with three expert clinicians, we have identified crucial points to be addressed in future work such as (a) training with multiple zoom factors, (b) forcing diversity of localization of tiles in each bag to avoid domination effect of large red artifacts, (c) pre-identifying cases with crumbling tissue and adapting the training procedure to give more importance to such challenging cases. Future investigations could also incorporate additional preprocessing steps(Reference Taylor-Weiner, Pokkalla and Han 31 ), stain normalization on Sirius-Red stained images(Reference Zheng, Jiang and Zhang 27 ), and extend the model to multi-class classification. We could also consider multimodal learning, since reasonable NAFLD diagnosis performance has been reported using DL on non-imaging data(Reference Sorino, Campanella and Bonfiglio 36 , Reference Li, Wang, Zhang, Zhang and Jiao 37 ).

Replication code and pretrained models are available at https://github.com/s-n-naik/Deep-Learning-for-Sirius-Red-Stained-Histopatholological-Images.

Authorship contribution

R.G. and E.D.A. conceptualized the study. S.N.N., R.F., P.M., and R.G. curated data. Investigation and data visualisation was performed by S.N.N. and E.D.A. supervised the project. E.D.A. administered the project and R.F., P.M., R.G., and E.D.A. provided resources. S.N.N. wrote the original draft and E.D.A. reviewed and edited the draft. All authors approved the final submitted draft.

Competing interest

The authors declare no competing interests exist.

Funding statement

This study was partially funded by the NIHR Imperial Biomedical Research Centre (BRC). The views expressed are those of the author(s) and not necessarily those of the NIHR or the Department of Health and Social Care.

Data availability statement

The replication code is available at https://github.com/s-n-naik/Deep-Learning-for-Sirius-Red-Stained-Histopatholological-Images.

References

Godoy-Matos, AF, Júnior, WSS & Valerio, CM (2020) NAFLD as a continuum: from obesity to metabolic syndrome and diabetes. Diabetol Metab Syndr 12, 120.CrossRefGoogle ScholarPubMed
Younossi, ZM, Koenig, AB, Abdelatif, D, Fazel, Y, Henry, L & Wymer, M (2016) Global epidemiology of nonalcoholic fatty liver disease—meta-analytic assessment of prevalence, incidence, and outcomes. J Hepatol 64, 7384.CrossRefGoogle ScholarPubMed
Lonardo, A, Mantovani, A, Petta, S, Carraro, A, Byrne, CD & Targher, G (2022) Metabolic mechanisms for and treatment of NAFLD or NASH occurring after liver transplantation. Nat Rev Endocrinol 18(10), 638650.CrossRefGoogle ScholarPubMed
Boyd, A, Cain, O, Chauhan, A & Webb, GJ (2020) Medical liver biopsy: background, indications, procedure and histopathology. Frontline Gastroenterol 11, 4047.CrossRefGoogle ScholarPubMed
Brunt, EM, Janney, CG, Di Bisceglie, AM, Neuschwander-Tetri, BA & Bacon, BR (1999) Nonalcoholic steatohepatitis: a proposal for grading and staging the histological lesions. Am J Gastroenterol 94, 24672474.CrossRefGoogle ScholarPubMed
Kleiner, DE, Brunt, EM, Van Natta, M, et al. (2005) Design and validation of a histological scoring system for nonalcoholic fatty liver disease. J Hepatol 41, 13131321.CrossRefGoogle ScholarPubMed
Nascimbeni, F, Ballestri, S, Machado, MV, et al. (2018) Clinical relevance of liver histopathology and different histological classifications of NASH in adults. Expert Rev Gastroenterol Hepatol 12, 351367.CrossRefGoogle ScholarPubMed
Bataller, R & Brenner, DA (2005) Liver fibrosis. J Clin Invest 115, 209218.CrossRefGoogle ScholarPubMed
Hu, W, Li, X, Li, C, et al. (2023 ) A state-of-the-art survey of artificial neural networks for whole-slide image analysis: from popular convolutional neural networks to potential visual transformers. Comput Biol Med 161, 107034.CrossRefGoogle ScholarPubMed
Srinidhi, CL, Ciga, O & Martel, AL (2020) Deep neural network models for computational histopathology: a survey. Med Image Anal 67, 101813.CrossRefGoogle ScholarPubMed
Zhou, X, Li, C, Rahaman, MM, et al. (2020) A comprehensive review for breast histopathology image analysis using classical and deep neural networks. IEEE Access 8, 9093190956.CrossRefGoogle Scholar
Ilse, M, Tomczak, J & Welling, M (2018) Attention-based deep multiple instance learning. In Proceedings of the 35th International Conference on Machine Learning, PMLR 80, pp. 21272136.Google Scholar
Campanella, G, Hanna, MG, Geneslaw, L, et al. (2019) Clinical-grade computational pathology using weakly supervised deep learning on whole slide images. Nat Med 25, 13011309.CrossRefGoogle ScholarPubMed
Li, J, Li, W, Gertych, A, Knudsen, BS, Speier, W & Arnold, CW (2019) An attention-based multi-resolution model for prostate whole slide imageclassification and localization. Preprint, arXiv:1905.13208.Google Scholar
Li, M, Wu, L, Wiliem, A, Zhao, K, Zhang, T & Lovell, B (2019) Deep instance-level hard negative mining model for histopathology images. In Medical Image Computing and Computer Assisted Intervention (MICCAI 2019), pp. 514522 [Shen, D et al., editors], Lecture Notes in Computer Science, 11764. Cham: Springer.CrossRefGoogle Scholar
Lu, MY, Williamson, DF, Chen, TY, Chen, RJ, Barbieri, M & Mahmood, F (2021) Data-efficient and weakly supervised computational pathology on whole-slide images. Nat Biomed Eng 5, 555570.CrossRefGoogle ScholarPubMed
Xie, C, Muhammad, H, Vanderbilt, CM, et al. (2020) Beyond classification: whole slide tissue histopathology analysis by end-to-end part learning. In Proceedings of the Third Conference on Medical Imaging with Deep Learning, PMLR 121, pp. 843856.Google Scholar
Bejnordi, BE, Veta, M, Van Diest, PJ, et al. (2017) Diagnostic assessment of deep learning algorithms for detection of lymph node metastases in women with breast cancer. JAMA 318, 21992210.CrossRefGoogle Scholar
Kwok, S (2018) Multiclass classification of breast cancer in whole-slide images. In Proceedings of Image Analysis and Recognition (ICIAR), pp. 931940. Cham: Springer.CrossRefGoogle Scholar
Arjmand, A, Christou, V, Tzallas, AT, et al. (2020) Transfer learning versus custom CNN architectures in NAFLD biopsy images. In 2020 43rd International Conference on Telecommunications and Signal Processing (TSP), pp. 480483. IEEE.CrossRefGoogle Scholar
Qu, H, Minacapelli, CD, Tait, C, et al. (2021) Training of computational algorithms to predict NAFLD activity score and fibrosis stage from liver histopathology slides. Comput Methods Programs Biomed 207, 106153.CrossRefGoogle ScholarPubMed
Bilal, M, Jewsbury, R, Wang, R, et al. (2023) An aggregation of aggregation methods in computational pathology. Med Image Anal 88, 102885.CrossRefGoogle ScholarPubMed
Reinhard, E, Adhikhmin, M, Gooch, B & Shirley, P (2001) Color transfer between images. IEEE Comput Graph Appl 21, 3441.CrossRefGoogle Scholar
Ruifrok, AC & Johnston, DA (2001) Quantification of histochemical staining by color deconvolution. Anal Quant Cytol Histol 23, 291299.Google ScholarPubMed
Janowczyk, A, Basavanhally, A & Madabhushi, A (2017) Stain normalization using sparse autoencoders (StaNoSA): application to digital pathology. Comput Med Imaging Graph 57, 5061.CrossRefGoogle ScholarPubMed
Cho, H, Lim, S, Choi, G & Min, H (2017) Neural stain-style transfer learning using GAN for histopathological images. arXiv:1710.08543.Google Scholar
Zheng, Y, Jiang, Z, Zhang, H, et al. (2020) Stain standardization capsule for application-driven histopathological image normalization. IEEE J Biomed Health Inform 25, 337347.CrossRefGoogle Scholar
Forlano, R, Mullish, BH, Giannakeas, N, et al. (2020) High-throughput, machine learning–based quantification of steatosis, inflammation, ballooning, and fibrosis in biopsies from patients with nonalcoholic fatty liver disease. J Clin Gastroenterol Hepatol 18, 20812090.CrossRefGoogle ScholarPubMed
Davison, BA, Harrison, SA, Cotter, G, et al. (2020) Suboptimal reliability of liver biopsy evaluation has implications for randomized clinical trials. J Hepatol 73, 13221332.CrossRefGoogle ScholarPubMed
Soon, GS, Liu, F, Leow, WQ, Wee, A, Wei, L & Sanyal, AJ (2022) Artificial intelligence improves pathologist agreement for fibrosis scores in nonalcoholic steatohepatitis patients. Clin Gastroenterol Hepatol 21, 19401942.e3.CrossRefGoogle ScholarPubMed
Taylor-Weiner, A, Pokkalla, H, Han, L, et al. (2021) A machine learning approach enables quantitative measurement of liver histology and disease monitoring in NASH. J Hepatol 74, 133147.CrossRefGoogle ScholarPubMed
Gawrieh, S, Sethunath, D, Cummings, OW, et al. (2020) Automated quantification and architectural pattern detection of hepatic fibrosis in NAFLD. Ann Diagn Pathol 47, 151518.CrossRefGoogle ScholarPubMed
Nam, D, Chapiro, J, Paradis, V, Seraphin, TP & Kather, JN (2022) Artificial intelligence in liver diseases: improving diagnostics, prognostics and response prediction. JHEP Rep 4, 100443.CrossRefGoogle ScholarPubMed
Berg, S, Kutra, D, Kroeger, T, et al. (2019) ilastik: interactive machine learning for (bio)image analysis. Nat Methods 16, 12261232.CrossRefGoogle ScholarPubMed
Hu, J, Shen, L & Sun, G (2018) Squeeze-and-excitation networks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pp. 71327141. IEEE.CrossRefGoogle Scholar
Sorino, P, Campanella, A, Bonfiglio, C, et al. (2021) Development and validation of a neural network for NAFLD diagnosis. Sci Rep 11(1), 113.CrossRefGoogle ScholarPubMed
Li, Y, Wang, X, Zhang, J, Zhang, S & Jiao, J (2022) Applications of artificial intelligence (AI) in researches on non-alcoholic fatty liver disease (NAFLD): a systematic review. Rev Endocr Metab Disord 23, 387400.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Overview of our bag construction and deep-learning pipelines for fibrosis scoring on Sirius Red histopathology WSIs. Two MIL bag construction methods are compared: (1) RED-1, which uses a priori knowledge on red appearance of fibrosis to generate 1 bag per WSI, and (2) RAND-n, which generates $ n $ bags per WSI only based on tissue content. Deep-learning pipeline (best model): SE-ResNet18 initialized with ImageNet pre-trained weights with gated attention aggregation to output a binary label on a given bag.

Figure 1

Table 1. Single-bag MIL architectures performance–pre-trained with ImageNet and using Dataloader RED-1 with bag size $ k=10 $.

Figure 2

Table 2. Effect of bag sizes $ k $ using SE-ResNet18 (Gated Attention) pre-trained with ImageNet on Dataloader RED-1 with $ n=1 $ bag per WSI.

Figure 3

Table 3. Comparing Dataloader RED-1 ($ n=1 $) versus Dataloader RAND-n with bag size $ k=10. $

Figure 4

Figure 2. True positive (TP) cases of severe fibrosis illustrated on tiles with the 10 highest attention weights. WSI color categories are indicated under each case. Arrows color-coding: blue = relevant pathological signs of bridging, purple = red artefacts, green = healthy portal tracts. (a–f) Cases with high attention correctly focused on fibrotic signs as well as on artifacts and portal tracts. (g) Case showing higher attention put on portal tracts sliced longitudinally than on fibrotic signs. (h) Case with strong blurring still showing high attention on fibrotic signs.

Figure 5

Figure 3. False positive (FP) cases of severe fibrosis. (a) Red artifacts: The tiles with two highest and two lowest attention weights show a focus on red pixels from a vein cut within the tissue. (b) Borderline case: This case was reconsidered to be a potential severe case by the expert clinician.

Figure 6

Figure 4. False Negative (FN) cases. (a) Faded red stain; (b) Poor-quality biopsy with tissue crumbling; (c,d) Borderline cases with sparse bridging signs shown in zoomed boxes: (c) most red pixels correspond to healthy portal tracts sliced longitudinally; (d) Very thin bridging patterns.