Skip to main content

ORIGINAL RESEARCH article

Front. Med., 07 September 2022
Sec. Pathology
Volume 9 - 2022 | https://doi.org/10.3389/fmed.2022.994467

Computational tumor stroma reaction evaluation led to novel prognosis-associated fibrosis and molecular signature discoveries in high-grade serous ovarian carcinoma

Jun Jiang1† Burak Tekin2† Lin Yuan3† Sebastian Armasu1 Stacey J. Winham1 Ellen L. Goode1 Hongfang Liu4*‡ Yajue Huang2*‡ Ruifeng Guo2*‡ Chen Wang1*‡
  • 1Department of Quantitative Health Sciences, Mayo Clinic, Rochester, MN, United States
  • 2Department of Laboratory Medicine and Pathology, Mayo Clinic, Rochester, MN, United States
  • 3Pathology Center, Shanghai General Hospital, Shanghai, China
  • 4Department of Artificial Intelligence and Informatics, Mayo Clinic, Rochester, MN, United States

Background: As one of the key criteria to differentiate benign vs. malignant tumors in ovarian and other solid cancers, tumor-stroma reaction (TSR) is long observed by pathologists and has been found correlated with patient prognosis. However, paucity of study aims to overcome subjective bias or automate TSR evaluation for enabling association analysis to a large cohort.

Materials and methods: Serving as positive and negative sets of TSR studies, H&E slides of primary tumors of high-grade serous ovarian carcinoma (HGSOC) (n = 291) and serous borderline ovarian tumor (SBOT) (n = 15) were digitally scanned. Three pathologist-defined quantification criteria were used to characterize the extents of TSR. Scores for each criterion were annotated (0/1/2 as none-low/intermediate/high) in the training set consisting of 18,265 H&E patches. Serial of deep learning (DL) models were trained to identify tumor vs. stroma regions and predict TSR scores. After cross-validation and independent validations, the trained models were generalized to the entire HGSOC cohort and correlated with clinical characteristics. In a subset of cases tumor transcriptomes were available, gene- and pathway-level association studies were conducted with TSR scores.

Results: The trained models accurately identified the tumor stroma tissue regions and predicted TSR scores. Within tumor stroma interface region, TSR fibrosis scores were strongly associated with patient prognosis. Cancer signaling aberrations associated 14 KEGG pathways were also found positively correlated with TSR-fibrosis score.

Conclusion: With the aid of DL, TSR evaluation could be generalized to large cohort to enable prognostic association analysis and facilitate discovering novel gene and pathways associated with disease progress.

Introduction

Ovarian cancer (OC) is one of the leading causes of mortality among cancers in women. Pathologically, ovarian cancer is divided into high-grade and low-grade carcinomas, and the high-grade carcinomas can be further classified into various histological subtypes, most commonly serous, endometrioid and clear cell. Among them, high-grade serous ovarian carcinoma (HGSOC) is the prevalent histotype and accounts for vast majority of ovarian cancer associated mortality (1). These patients often present with rapid clinical progression, disseminated peritoneal metastasis, distant metastasis, and resistance to treatments. On the other hand, low-grade ovarian carcinomas usually present with slow-progressing diseases and are associated with much lower mortality but protracted clinical courses. Diagnostically, low-grade ovarian carcinomas can be difficult to distinguish from ovarian borderline tumors with much more indolent clinical behavior, while sometimes may share overlapping histological features with aggressive high-grade carcinomas. While all malignant OCs regardless histology types are treated similarly using platinum-based front-line chemotherapies, different surgical resections and chemotherapy treatments options could be applied to different histologic subtypes. In clinical diagnosis, recognizable histological features play a critical role in differentiating these subtypes.

Although histological diagnosis of HGSOC has been well-established, many studies have shown highly heterogenous clinical courses in these patients (24). Interestingly, pathologists have long observed the high variability of tumor associated stroma reaction in HGSOCs in daily practice (4, 5). Similar to a process of normal wound healing, the tumor-stroma reaction (TSR) in cancer has been associated with increased extracellular matrix and production of growth factors to facilitate recovery growth of injured tissues (6). In ovarian tumors, histopathological examination of tumor-stroma reaction is critical to differentiate low-grade serous carcinoma from serous borderline serous tumor (SBOT), with the latter lacking tumor triggered stroma reaction. More importantly, tumor-stroma reaction has been reported to facilitate tumorigenesis and associated with prognostic differences in many solid cancers such as cholangiocarcinoma, pancreatic cancer, melanoma, and OC (710). Though numerous studies have demonstrated that the interactions between tumor cells and stroma play a critical role in cancer progression and metastasis across multiple cancer types (1113), the association of histological feature of stromal reaction with molecular mechanism is still underexplored. One of major reasons responsible for this gap is the lack of quantitative evaluation of TSR in solid tumors. In daily pathology practice and many research studies, TSR were examined by manually reviewing H&E-stained slides by individual pathologist, which is highly subjective and labor-intensive. Interobserver variability remains a main challenge thus limiting large-scale investigation of TSR. More importantly, evaluation of TSR by pathologist relies on personal experience, while an unbiased quantification becomes unrealistic which may cause heterogenous quality of the TSR scoring data with poorly reproducibility.

With the advancement of digital pathology, there has been substantial interest in exploring the role of quantitative attributes computationally extracted from H&E-stained whole slide images (WSI). Li et al. (14) introduced a digital pathology-based pipeline to early-stage estrogen receptor-positive invasive breast cancers for association analysis. Their results suggested that the orientation disorder of collagen fiber is prognostic for early-stage breast cancer (14). Geessink et al. (15) trained a deep learning model to segment relevant tissue types in rectal cancer histology and subsequently calculate tumor-stroma ratio for intra-tumoral stroma. Their results showed that tumor-stroma ratio is an independent prognosticator in rectal cancer when assessed automatically in user-provided stroma hot-spots (15). Failmezger et al. (16) introduced topological features extraction method to quantify stromal recruitment for immunosuppression in melanoma histology using graph based spatial model. This research revealed that tumors with high stromal clustering and barrier had reduced expression of pathways involved in naïve CD4 signaling, MAPK, and PI3K signaling, and indicated that computational histology-based stromal phenotypes within the tumor microenvironment are significantly associated with prognosis and immune exclusion in melanoma (16).

In this paper, we present a digital pathology-based pipeline that is able to prognosticate patient survival by estimating degree of TSR directly from multiple aspects of digitized H&E images. Specifically, the automated pipeline consists of image processing techniques combined with several machine learning models trained from pathologists’ annotations. Serving as the cores of the pipeline, the trained models were used to identify tumor-associated stroma regions, from which we subsequently predicted TSR scores with H&E images as inputs. As shown in Figure 1, the developed pipeline was applied to our research cohort to establish associations between tissue-level features, prognosis, and molecular pathways of HGSOC.

FIGURE 1
www.frontiersin.org

Figure 1. Overview of our research workflow. (A) Slide scanning and annotation. (B) Tumor-stroma segmentation and TSR score estimation. (C) Tumor stroma interface area identification. (D) Tissue-level feature summarization. (E) Association analysis.

Materials and methods

Cohort summary

Our research cohort consists of 291 HGSOC and 15 serous borderline ovarian tumor (SBOT) cases ascertained at the Mayo Clinic between 1994 and 2009. The SBOT cases should not have significant TSR by definition, therefore served as the negative control in both TSR score prediction and evaluation. As the cohort selection criteria, all the cases had retrievable clinical, molecular, and tissue blocks. Survival data were obtained from the Mayo Clinic Tumor Registry, electronic medical records. Gene expression profiles and histological information, including tissue sites (primary or metastasis) and tumor stage, were collected from the EHR system in Mayo Clinic. The cohort characteristics were summarized in Table 1. All the slides were scanned in Pathology Research Core at Mayo Clinic with a digital whole-slide scanner (Aperio Scanscope XT). To preserve cell and tissue details, the slides were scanned with 40× resolution (pixel size: 0.25 um). Imaging quality was manually checked by histology technicians when the slides were scanned. All patients provided informed consent for use of their tissues and data in research; all protocols were approved by the Mayo Clinic Institutional Review Board.

TABLE 1
www.frontiersin.org

Table 1. Research cohort statistics.

Pathologist-guided image annotation

According to the consensus of three pathologists, three types of histopathological evaluation metrics (sub-scores) were defined as criteria to characterize the extent of TSR: (i) increased fibrosis, characterized by collagen deposition, (ii) increased stromal cellularity due to fibroblastic and/or myofibroblastic proliferation, and (iii) orientation of stromal cells (14). Based on histopathologic examination of H&E-stained slide areas, the sub-TSR scores were assessed as 0 (none/weak), 1 (intermediate) or 2 (strong) (Supplementary Figure 1). The criteria were chosen to reflect the histopathologic changes commonly observed and evaluated in the clinical practice. With regard to the criterion of stromal cellularity, a minimal density of (myo)fibroblasts was assigned a score of 0, while a score of 2 was given if the area occupied by (myo)fibroblasts exceeded the area occupied by the acellular stroma in a given field. As for the fibrosis criterion, minimal deposition of fine collagen fibers with significant fiber spacing was assigned a score of 0, whereas dense collagen deposition with sclerosis was classified as 2. In terms of orientation of stromal cells, a relatively linear, unidirectional orientation was assigned a score of 2, while a haphazard orientation without appreciable directionality was scored as 0.

To create an annotated dataset for model training, five HGSOC slides were randomly selected. Within each slide, five most representative regions of interest (ROIs) were circled by pathologists in tumor-stroma interface regions, which were areas where the borders of the tumor islands came into close proximity with the surrounding stromal areas and the stroma exhibited morphologic characteristics different than the non-neoplastic ovarian stroma (Figure 1A). Two experienced pathologists were invited to annotate three TSR scores using an interactive tool named QuPath (17). In each slide, five most representative ROIs were circled by pathologists in tumor-stroma interface regions. The size of each ROI was at least 1,024*1,024 pixels. Within each ROI, polygons were used to annotate homogeneous regions with the same TSR scores. Sub-regions with the same TSR scores were labeled to the same category (Supplementary Figure 1B).

Dataset preprocessing

Using a framework developed in our previous work (18), TSR annotations were parsed using Groovy script within QuPath and converted into a pair of image and annotation mask for each annotated ROI. For the convenience of visualization, annotation masks were encoded from dark to light R/G/B colors for each TSR scoring metric (Supplementary Figures 1D–F). To extract regular size of images and annotation masks for model training and evaluation, a 256*256 pixel sliding window was applied to the ROIs. Taking full capacity of pathologists’ annotations, the sampling stride was set to 128 for the aim of creating augmented/enlarged dataset. In total, 11,240 image patches with TSR annotations were obtained from HGSOC cases.

Considering that SBOT confirmed cases were free of significant TSR, we proposed to train the TSR prediction model with image samples from SBOT WSIs as negative controls. Thus, image patches from SBOT cases were also prepared and labeled to TSR score zero. We randomly selected five slides from our previous research (18), from which 1,405 image patches from stroma regions were randomly extracted and added into “annotated” dataset. In total, 7,025 image patches were obtained from SBOT cases.

To identify tumor-stroma interface areas and quantify TSR inside these regions, we repurposed TSR annotations for tumor-stroma segmentation modeling. The stroma region is defined as all the tissue region except the tumor region. Within each annotated ROI, overall stroma regions were obtained by merging all three different TSR score regions, while the tumor regions were defined as the remaining tissue regions inside the ROIs. The same sliding-window sampling strategy was used to extract image patches for segmentation modeling.

Tumor-stroma tissue segmentation

In order to identify tumor-stroma interface areas where TSR occurs, tumor and stroma regions were segmented using a deep learning neural network named Mask-RCNN (19), as shown in Figure 1B. Mask -RCNN was preferred in this study as it has been used in many histological image processing tasks (20, 21) and was reported to be more robust than the U-Net for image segmentation (22, 23). The hyperparameters of Mask-RCNN, such as the dimension of convolutional layers (input dimension = 256 × 256), learning rate (lr = 0.01) and RPN anchor scales (RPN = [8, 16, 32, 64, 128]) were modified to adapt to our image segmentation task. Taking annotated dataset described in Data Preprocessing section, images from HGSOC cases were shuffled and divided into training, validation, and testing groups (3,317:1,105:1,110). From the training subset, the input layers of Mask-RCNN took both original images and tumor-stroma multilabel masks as training samples. In the training process, tumor and stroma areas were iteratively proposed by a sub-structure of Mask-RCNN named Region Proposal Network (RPN). Fully connected network layers were concatenated to the forehead layers to identify differences between ground truth (annotations) and proposed segmentation masks. By minimizing the differences, Mask-RCNN was trained to segment tumor and stroma within WSIs. To facilitate training convergence, weights from the model pretrained with Coco dataset (24) were loaded to our model as initial settings and fine-tuned on our training dataset. Specifically, the model trained at the 315th epoch reached the lowest loss in the validation set.

In the testing phase, tumor stroma regions were segmented and saved as multilabel masks for each image patch from the hold-out testing dataset using trained tumor-stroma segmentation models. Three commonly used evaluation metrics for image segmentation tasks (25, 26) were calculated to measure concordance between model prediction and ground truth, including DSC (Dice similarity coefficient), IoU (intersection of union) and AP (averaged precision).

Tumor-stroma reaction scoring modeling

In our work, TSR score estimation was formulated into an image classification problem. In other words, different TSR scores corresponded to different image categories. We employed a commonly used DL network architecture named VGG16 (27) as our image classification model, as this model also achieved an encouraging performance on identifying tumor infiltrating lymphocytes (TIL) (28). To estimate three TSR scores, three VGG16 models (Figure 1B) were trained with the combinational dataset (from both HGSOC and SBOT cases) established in data preprocessing steps. The annotated dataset was divided into training, validation, and testing subsets (12,743:2,247:2,249, 6,000 training samples were from SBOT as negative control). During the training phase, the training dataset was divided into batches (32 samples per batch) to meet the computational resource limitations. The maximum training epoch was set to 30. At the end of each training epoch, training loss was calculated on the validation dataset. To avoid overfitting, the training process was set to stop when the loss variation is less than 10–3 within four epochs. To increase generalizability and avoid bias from different H&E-staining conditions, training image dataset was augmented using linear image transformation, such as rotation and flipping. With the same training strategy, three VGG16 image classification models were trained independently to estimate fibrosis, cellularity, and orientation TSR scores, respectively (Figure 1B). During the testing phase, three TSR scores were estimated for each input image from the hold-out testing dataset using the three trained VGG16 models. The performances of the three models were evaluated by comparing the concordance between model estimation and human annotation.

Extrinsic model evaluation

Before applying our models to the entire research cohort, it was essential to evaluate model performances on an independent dataset (extrinsic) as our DL models were trained and evaluated using annotated images from ROIs (intrinsic). To this end, we developed an interactive evaluation tool (Supplementary Figure 2) for the assessment on a dataset that was independent of annotated ROIs and WSIs. For the sake of pathologists’ convenience, the original image as well as corresponding tumor-stroma segmentation and TSR scores predicted by the trained models were loaded into this evaluation tool. With the original image in the center, eight neighborhood image patches were also shown in the user interface as additional references for pathologists to make an accurate judgment. The concordance (average precision) between pathologists’ evaluations and predictions were recorded as pathologists proceeded reviewing by clicking buttons and checking boxes.

To create an independent dataset for extrinsic evaluation, image patch-level TSR score distribution was calculated for each slide. Slides with ultra-high (TSR = 2) or low (TSR = 1) TSR score ratio were selected to epitomize the performances of our trained models. Unannotated cases were selected based on TSR score ratio distribution, only the upper 5% quantile and low 5% quantile were included for manual evaluation. To limit the number of images to be evaluated, we randomly selected at least 10 but less than 30 images from each slide within tumor stroma interfacing regions.

Applying deep learning models to research cohort

Trained DL models were generalized to the research cohort according to the following procedures.

Tissue detection and patch extraction

To cut down the computational cost, only the tissue regions (foreground of WSIs) were included in the testing phase. Tissue regions were detected within a down-sampled whole slide image (down-sample rate = 128), and then image patches were extracted from tissue regions accordingly. To detect foreground within down-sampled whole slide images, the color space was converted from RGB to LAB, and then threshold method was applied to L channel for tissue detections. As a commonly used foreground detection method for bright field whole slide images, feasibility of this method has been proved in our previous research (18, 29). By mapping the coordinates of pixels within detected tissue regions back to original WSI, image patches (256*256 pixels) were extracted in a tiling manner. With the same threshold-based method, foreground within the original resolution of image patches was detected. If more than 50% of the patch is background, it was excluded from our study. In the end, the extracted images were fed into our trained models for patch-level predictions.

Tumor-stroma interface area identification

To investigate TSR in invasive tumor front, tumor-stroma interface areas were identified based on tumor-stroma segmentation results. By down-sampling (r = 1/128) and stitching patch level results back to their original locations, slide level tissue segmentation (tumor vs. stroma) were reconstructed. Within slide level tumor-stroma segmentation results, tumor-stroma interface areas were identified using serials of image morphological and logical operations (Figure 1C and Supplementary Figure 5). The calculation process can be formulated as follows:

T u m o r c o r e = C ( I T , S )
S t r o m a c o r e = C ( I S , S )
R O I i n t e r f a c e = a n d ( x o r [ D ( S t r o m a c o r e , S ) , E ( S t r o m a c o r e , S ) ] ,
S t r o m a c o r e , D ( T u m o r c o r e , S ) ) ,

in which, IT and IS denotes the tumor and stroma multilabel image from slide level segmentation, respectively. S denotes the structural elements for morphological operations, including closing C(I, S), erosion E(I, S) and dilation D(I, S) (30).

To evaluate the accuracy of this automatic tumor-stroma localization method, the multilabel masks of interface areas and the counterpart WSIs were shown side-by-side and reviewed by our pathologists. Specifically, top five largest connected components were detected as the representative sub-regions for detailed reviews (31). The misidentifications were recorded for quantitative assessment metrics calculation.

Slide level feature summarization

To enable association analysis, TSRs were summarized to abstract slide level descriptors (Figures 1D,E). Since the ROI (tumor-stroma interface area) size varies from case to case, we used mean and standard deviation to denote slide level characteristics. Normalized distributions of TSR scores were calculated by counting TSR scores of each image patch within tumor-stroma interface regions. The entire assembled workflow was generalized to all the slides in our cohort. The summarized features were prepared for association analysis.

Downstream association analysis

The summarized TSR characteristics were associated with clinical and molecular information. In our work, only HGSOCs were included for downstream analysis.

Clinical associations

In HGSOC cases, for each TSR score (Fibrosis, Cellularity, and Orientation), median split was used to divide patients into two groups (i.e., score-high and low) to facilitate categorical comparisons. For univariate and multivariable [adjusted for age, FIGO stage (IV vs. III), and residual tumor after primary debulking surgery] survival analysis, a Cox proportional hazards regression model was used, and hazard ratios (HRs) and associated 95% confidence intervals (CIs) were estimated. All statistical tests were two-sided, and a P-value of less than 0.05 was considered statistically significant.

Molecular associations

Tumor gene expression profiles were measured using Agilent Whole Human Genome 4x44K Expression Arrays and processed as previously described (2, 4). For gene-level association analysis, normalized expression levels of each gene were correlated with each TSR score from the same tumors using Spearman Rank correlation. For over-representation pathway analysis purposes, genes with positive and negative correlations with each TSR score (nominal p-value < 0.05) were analyzed using DAVID bioinformatics tool (32, 33), to reveal pathways statistically enriched in correlated gene sets. False discovery rates were computed to correct for multiple hypothesis testing.

Results

Tumor-stroma segmentation

The developed tumor stroma segmentation model identified the tumor vs. stroma region within both HGSOC and SBOT WSIs (Figure 2A). Figure 2B demonstrates the whole slide level segmentation results by stitching patch-level segmentation results together according to patch locations. Different tissue types were colored according to predicted categories. More examples of slide-level segmentations were shown in Supplementary Figure 3.

FIGURE 2
www.frontiersin.org

Figure 2. Examples of tumor-stroma segmentation and TSR scoring results. (A) Original WSIs, with HGSOC and SBOT each; (B) tumor-stroma segmentation, tumor and stroma were encoded with cyan and yellow; (C) TSR scores measured from three metrics, including fibrosis (Red), cellularity (Green) and orientation (Blue). Each metric was encoded from dark to light color, denoting TSR score from low to high. *For better visualization, TSR scores within all stroma regions were shown, but only the tumor-stroma interface regions were included for analysis.

In the hold-out testing subset (N = 1,110) reserved for segmentation accuracy assessment, DSC, IoU and AP achieved 93.5, 88.65, and 95.34%, respectively (Figure 3A). Moreover, we observed that IoU and DSC dramatically decreased if there was a tissue type misdetection (Figure 3B). Our evaluation also suggested that IoU and DSC were highly correlated to each other (Figure 3C), and AP could be a more suitable metric for measuring our segmentation accuracy. Since the hold-out testing set is from annotated ROIs, our evaluation results suggested that the trained tumor-stroma segmentation model performed well in intrinsic cases.

FIGURE 3
www.frontiersin.org

Figure 3. Tumor-stroma segmentation evaluation. (A) Boxplot of three evaluation metrics, including IoU, AP and DSC. (B) Examples of segmentation. Red arrows point to missed targets in segmentation. (C) Correlation of three evaluation metrics. Each dot represents an image sample. Linear regression was used to calculate correlation.

Based on the criteria mentioned in the methods, 15 unannotated slides were identified, from which 615 image patches were sampled for independent tumor-stroma segmentation and TSR score evaluation. By analyzing the review records from the extrinsic evaluation tool (Supplementary Figure 2), our model achieved 90.6% accuracy, indicating that pathologists generally agreed with our tumor-stroma segmentation performance within independent WSIs. It is noteworthy that our trained model can be applied to the entire research cohort to generate tumor-stroma segmentation across the whole slide for the downstream analytical steps.

Based on our tumor-stroma segmentation results, some WSIs with low stroma tissue areas were identified. By checking the original images of these cases within QuPath, the pathologists confirmed that our tumor-stroma segmentation results were accurate, as the tumor islands occupied the majority of these slides, while the stromal areas consisted of a significant number of adipocytes, necrosis, and/or hemorrhage, with minimal collagenous and cellular stroma [Supplementary Figure 4 case (4)].

Tumor-stroma reaction scoring and evaluations

Using our trained models, three TSR scores were predicted for each patch inside the stroma regions based on the tumor stoma segmentation results. The spatial overview of slide-level TSR was reconstructed by mapping three TSR scores to different colors with ranked saturation (Figure 2C). More examples of slide level segmentations are shown in Supplementary Figure 4. More zoomed in results can be found in our GitHub repository.

Our qualitative results suggested that the predicted TSR scores were not even in stroma regions of all HGSOC. Heterogeneity between regions and cases were high, as shown in Figure 2C and Supplementary Figure 4. Confusion matrices were calculated to quantitatively evaluate model performances on the hold-out testing dataset (N = 2,249). The results indicated that our model achieved accurate TSR score estimation, especially in predicting fibrosis (>90%), as shown in Figure 4A. In the extrinsic evaluation, with an average precision over 82.8%, the results suggested that the trained model can be generalized to our research cohort for objective TSR scoring.

FIGURE 4
www.frontiersin.org

Figure 4. Tumor-stroma reaction (TSR) scoring evaluation. (A) Confusion matrix of intrinsic evaluation. (B) Violine plot for three TSR metrics within HGSOC vs. SBOT. Majority of SBOT images have low TSR score, no matter in which metric.

We also observed false positives (TSR scores > 0) in some region from three SBOT cases, as illustrated Supplementary Figure 4 case (3). By mapping the TSR scores back to the WSIs and observing with high resolution, we identified that these false-positive predictions were presumably due to these slides having mostly non-neoplastic ovarian stroma, which inherently has a relatively cellular and fibrotic composition. Since our model is mainly trained on annotated HGSOC regions, the trained model did not capture texture patterns within normal SBOT regions.

Violin plots illustrated the distributions of TSR score ratios per case. Figure 4B indicates that the majority of image patches had low TSR scores (TSR = 0), regardless of the diagnosis being SBOT or HGSOC. However, compared to SBOT cases, HGSOC cases were more likely to have higher TSR scores (>1), especially for fibrosis score. We observed a significant proportion of image patches from HGSOC cases had fibrosis TSR scores of 1. After checking the training dataset, we confirmed that half (3,339 out of 6,743) of the annotated images of HGSOC cases had moderate fibrosis TSR scores, indicating that the ambiguity of fibrosis scores could be high, and our TSR scoring model was trained to match pathologists’ interpretations.

Identification of tumor-stroma interface regions

Our tumor-stroma interface region identification strategy identified five regions within each testing slide (Supplementary Figure 5). The proposed interface regions were localized and overlayed to the WSIs. According to pathologists’ manual review, 83.4% (207 out of 250) proposed tumor-stroma regions were confirmed to be tumor-stroma interface area. After checking the falsely proposed tumor-stroma interface regions, we found flaws of tumor-stroma segmentation in those regions to be responsible for the failure, indicating that the tumor-stroma interface identification relies heavily on tumor-stroma segmentation results. To facilitate replication of our work, all the code for this paper is public available on GitHub.1 The pretrained models for tumor-stroma segmentation and tumor-stroma reaction prediction can be obtained via contacting authors.

Tumor-stroma reaction clinical and molecular associations

All three TSR scores were significantly elevated in HGSOC cases vs. SBOT cases (p < 0.001, Figure 5A). Moreover, in HGSOC cases, higher fibrosis score (>median) was significantly associated with worse survival (p = 0.02; Figure 5B), and the prognostic association remained significant (p = 0.04; Figure 5C) after multivariate adjusting for other established prognostic factors (age at diagnosis, stage, and residual tumor after surgical debulking). In order to gain further insight into possible molecular mechanisms associated with each TSR score, gene-level correlations were computed between mRNA level of each gene and TSR score from the same tumors; and significant associations were found in two correlations: (1) correlation between fibrosis and molecular findings, and (2) correlation between cellularity and molecular findings (Figure 5D). Further genetic analysis suggested different molecular bases between the three TSR scores. Through pathway enrichment analysis, genes positively correlated with TSR-fibrosis score were found to be enriched in 14 KEGG pathways [FDR (false discovery rate) < 5%], which are mostly associated with cancer signaling aberrations. On the other hand, genes positively correlated with the TSR-orientation score were enriched in 79 KEGG pathways, with leading significant pathways implicated with immune response (Supplementary Table 1). In contrast, genes having positive correlations with TSR-cellularity score were only significantly enriched with one KEGG pathway (hsa01100: Metabolic pathways; FDR = 0.04). Detailed molecular association results including gene- and pathway-level results were shown in Supplementary Table 1.

FIGURE 5
www.frontiersin.org

Figure 5. Prognosis and molecular associations of fibrosis score. (A) Tumor-stroma reaction (TSR) score boxplots for HGSOC and SBOT groups. (B) Overall survival differences between fibrosis high vs. low. (C) The prognostic association for other established prognostic factors (age at diagnosis, stage and residual tumor after surgical debulking). (D) Correlation between fibrosis/cellularity/orientation and molecular findings.

Discussion

In this study focusing on digital analysis of reactions between tumor and stroma, combined with a critical pathology review using subjective scoring systems, we demonstrated highly concordant computational prediction based on VGG16 DL structure with training annotations and independent validations by multiple pathologists in HGSOC. Conventionally, TSR is often lumped as an overall subjective assessment by pathologists. Herein we further dissected it into three essential aspects of TSR and examined their individual and combined significance by digital detection. With a series of digital detection and quantification procedures including tumor-stroma interface area detections, the trained DL model has been successfully generalized to a large OC cohort from a single institution consisting of nearly 300 patients with long-term clinical follow-up and tumor transcriptome data. Interestingly, among the three aspects of TSR, the data revealed significant prognosis association only with fibrosis score. This is the first study demonstrating the outstanding significance of fibrosis over other TSR pathological features in HGSOC, indicating these features did not carry the equal weight regarding clinical significance. At the design phase of this study, the orientation score was selected as one of the three quantification criteria, mainly to further characterize the fibroblastic and/or myofibroblastic proliferation. In fact, collagen fiber organization has been associated with prognosis in breast cancer in the literature and DL approaches have been employed to quantify this as a histomorphometric feature (14, 34). However, our team observed this criterion to be a highly subjective one. In addition, the orientation score did not reveal any significant prognosis associations in our transcriptome association analysis, calling the usefulness of this criterion with regard to ovarian cancer TSR assessment into question. These observations warrant further study on individual components of the TSR, as well as aberrant gene- and pathway-level activities associated with different digital TSR scores.

Of note, the design of the TSR scoring in our study was not specific or limited to HGSOC, and the same histological principle can be applied to in majority types of solid cancers. Therefore, the digital platform developed by our study can be potentially generalized and applied to various tumor types. These findings highlight potentials of powerful DL approaches to generalize digital pathology-based predictions for large-scale translational research and enable molecular discoveries to better understand tumorigenesis and cancer progression.

To consolidate our discoveries, we considered including other pre-existing cases into our research cohort. However, we found that images scanned at different times could be dramatically different in hue even if they were from the same patient, same institution (Mayo Clinic) and shared the same staining and image acquisition protocol. This inconsistency may be due to multiple technical variables, including scanner settings and/or age of the H&E-stained slides. Since these batch effects in pathology image data could be hidden variables in deep learning digital pathology that compromise the accuracy of classification systems (35, 36), we opted for not including our previous HGSOC data into this research. Many previous studies introduced color normalization methods to minimize staining inconsistencies (29, 37). Though it is hard to measure the preservation of diagnostic information after image transformation, many integrative studies investigating cancer subtype classification and prognosis association achieved optimistic performances by introducing image normalization (38, 39). From this point of view, color normalization could be beneficial for assorted research cohort from miscellaneous data sources, especially from multiple institutions. Meanwhile, we also noticed that some investigators improved their model performance by synthesizing images using generative adversarial network (GAN) (40, 41), which could be another potential way to enhance the generalizability of our TSR estimation model.

Although our tumor-stroma interface region detection relies on patch-level tumor-stroma segmentation, the strategy we introduced (Supplementary Figure 5) could partially offset this limitation. To generate a slide-level overview of tissue context (tumor vs. stroma) for image morphological manipulations, patch-level tumor-stroma segmentation results were down-sampled and stitched back to their original locations. In this process, the tissue type (tumor or stroma) in slide-level was determined by the dominated component of patch-level segmentation. In other words, for tumor-stroma interface region detection task, tumor-stroma segmentation results were not required to achieve pixel-level accuracy.

We observed that our TSR scoring models highlighted some regions with a potentially high TSR in five SBOT cases, for example, case (3) in Supplementary Figure 4. The main reason contributing to this flaw is that our models were not trained to differentiate normal vs. abnormal ovarian stroma. We anticipate that our TSR scoring pipeline can achieve a better estimation if models can be trained using extra normal vs. abnormal ovarian stroma annotations. Meanwhile, we acknowledge that using TSR as the sole measurement is not enough to describe the complex tumor micro-environment (TME). It has been reported that TILs can also be assessed with the aid of digital pathology in advanced-stage, HPV-negative head and neck tumors (42). We will introduce more interpretable measurements for pathology image metadata summarization, which will bring more opportunities for novel discoveries. To integrate cellular level features into large cohort analysis, we also plan to introduce more advanced cell segmentation modules to our workflow for better cell level representations (43). As reported in our previous work (44), over- and under-segmentations lead to inaccurate downstream analysis impute to erroneous features calculated based on them.

Another imperfection of our study is way we summarize predicted TSR score to slide level for association analysis. For the sake of simplicity, we employed simple statistics and assigned equal weights to patch level TSR predictions. However, the relative importance of tissue regions contributing to the diagnosis could be dramatically different depending on tissue context. In this study, we introduced tumor-stroma interface area identification methods that were aimed at mimicking pathologists’ diagnoses. This strategy is simple and works well in most cases; however, it highly depends on the tumor-stroma segmentation accuracy. We noticed that some studies proposed to introduce multi-resolution analysis for capturing subtle tissue features within different WSI scales (45, 46). Attention based deep neural networks (47) are also tangible options for locating diagnostically relevant regions and assign those regions with higher weights in automatic analysis. We will consider these solutions to fill the gaps in our current work.

Further limitations of the current study include that all the samples were from a single institution, requiring further validations in external cohorts. H&E imaging-based digital pathology studies as such may be also affected by paraffin block preservation protocols and digital scanning parameter settings, which could lead to model differences when generalizing to WSI samples collected and scanned following different protocols. Computational developments and evaluations will be made to address these challenges.

Conclusion

Our developed system achieved encouraging performances in tissue segmentation and TSR score predictions and generalized successfully to a large single-institution OC cohort, resulting in novel discoveries of clinical prognosis associations and molecular findings implicated in different TSR scores.

Data availability statement

The original contributions presented in this study are publicly available. The source code can be found here: https://github.com/smujiang/TumorStromaReaction. The pre-trained model can be obtained by contacting authors.

Ethics statement

The studies involving human participants were reviewed and approved by Mayo Clinic Institutional Review Board. The patients/participants provided their written informed consent to participate in this study. Written informed consent was not obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

RG and YH conceived of the presented histological concept. BT, LY, and RG defined data annotation criteria, performed annotation, and validated model performance. JJ and CW conceived and performed majority of experiments and manuscript drafting. SA conducted cohort characteristic statistics. SW contributed to the final version of the manuscript. HL provided funding resources and supervised the deep learning parts of this project. EG also supervised the funding of this work. All authors discussed the results and contributed to the final manuscript.

Funding

This research was partially supported by the Mayo Clinic Ovarian Cancer SPORE (P50 CA136393), National Institutes of Health’s National Institute (R01EB019403), Mayo Clinic Center for Individualized Medicine, and Mayo Clinic Department of Laboratory Medicine and Pathology, Division of Anatomic Pathology Intramural Research Fund.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2022.994467/full#supplementary-material

Supplementary Figure 1 | Tumor-stroma reaction (TSR) annotation. (A) Original image within ROI selected for annotation. (B) Annotated ROIs. Polygons were used to label regions with different TSR scores. (C) Legend of three TSR score measurements. (D–F) Parsed annotations. TSR scores were encoded into R/G/B colors to represent three measurements (fibrosis, cellularity, and orientation), respectively; Panels (G,H) are two zoom in examples. Panel (G) was annotated as Fibrosis = 2, Cellularity = 1, Orientation = 1; Panel (H) was annotated as Fibrosis = 1, Cellularity = 2, Orientation = 2.

Supplementary Figure 2 | Interactive tool for extrinsic evaluation. Source code available in our GitHub. Buttons and checkboxes on the right are clickable, pathologists’ interactions were recorded for extrinsic evaluation.

Supplementary Figure 3 | Extra examples of tumor-stroma segmentation results, including five HGSOCs, five SBOTs and their tumor stroma segmentation results.

Supplementary Figure 4 | Extra examples (two HGSOC and two SBOT) of TSR scoring results. TSR scores measured with fibrosis (Red), cellularity (Green), and orientation (Blue). From dark to light, TSR scores were encoded into R/G/B colors. *For better visualization, TSR scores within all stroma regions were shown, but only the tumor-stroma interface regions were included for analysis.

Supplementary Figure 5 | Tumor-stroma interface area identification. (A) Original WSI and tumor-stroma segmentation results. (B) Morphological and logical operations were conducted on tumor-stroma segmentation for localizing tumor-stroma interface regions. (C) Proposed ROIs (red rectangles) for TSR score summarization.

Supplementary Table 1 | Detailed molecular association results, including top 10 genes positively/negatively associated with TSR-Fibrosis score, and pathways positively associated with Fibrosis and Orientation scores, respectively.

Footnotes

  1. ^ https://github.com/smujiang/TumorStromaReaction

References

1. Lisio M-A, Fu L, Goyeneche A, Gao Z-H, Telleria C. High-grade serous ovarian cancer: basic sciences, clinical and therapeutic standpoints. Int J Mol Sci. (2019) 20:952. doi: 10.3390/ijms20040952

CrossRef Full Text | Google Scholar

2. Wang C, Armasu SM, Kalli KR, Maurer MJ, Heinzen EP, Keeney GL, et al. Pooled clustering of high-grade serous ovarian cancer gene expression leads to novel consensus subtypes associated with survival and surgical outcomes. Clin Cancer Res. (2017) 23:4077–85. doi: 10.1158/1078-0432.CCR-17-0246

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Talhouk A, George J, Wang C, Budden T, Tan TZ, Chiu DS, et al. Development and validation of the gene expression predictor of high-grade serous ovarian carcinoma molecular subTYPE (PrOTYPE). Clin Cancer Res. (2020) 26:5411–23. doi: 10.1158/1557-3265.OVCA19-A03

CrossRef Full Text | Google Scholar

4. Konecny GE, Wang C, Hamidi H, Winterhoff B, Kalli KR, Dering J, et al. Prognostic and therapeutic relevance of molecular subtypes in high-grade serous ovarian cancer. J Natl Cancer Inst. (2014) 106:dju249. doi: 10.1093/jnci/dju249

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Murakami R, Matsumura N, Mandai M, Yoshihara K, Tanabe H, Nakai H, et al. Establishment of a novel histopathological classification of high-grade serous ovarian carcinoma correlated with prognostically distinct gene expression subtypes. Am J Pathol. (2016) 186:1103–13. doi: 10.1016/j.ajpath.2015.12.029

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ohtani H. Stromal reaction in cancer tissue: pathophysiologic significance of the expression of matrix-degrading enzymes in relation to matrix turnover and immune/inflammatory reactions. Pathol Int. (1998) 48:1–9. doi: 10.1111/j.1440-1827.1998.tb03820.x

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sirica AE, Gores GJ. Desmoplastic stroma and cholangiocarcinoma: clinical implications and therapeutic targeting. Hepatology (Baltimore, Md). (2014) 59:2397. doi: 10.1002/hep.26762

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wang LM, Silva MA, D’Costa Z, Bockelmann R, Soonawalla Z, Liu S, et al. The prognostic role of desmoplastic stroma in pancreatic ductal adenocarcinoma. Oncotarget. (2016) 7:4183. doi: 10.18632/oncotarget.6770

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Busam KJ, Mujumdar U, Hummer AJ, Nobrega J, Hawkins WG, Coit DG, et al. Cutaneous desmoplastic melanoma: reappraisal of morphologic heterogeneity and prognostic factors. Am J Surg Pathol. (2004) 28:1518–25. doi: 10.1097/01.pas.0000141391.91677.a4

CrossRef Full Text | Google Scholar

10. Davidson B, Trope CG, Reich R. The role of the tumor stroma in ovarian cancer. Front Oncol. (2014) 4:104. doi: 10.3389/fonc.2014.00104

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Bremnes RM, Dønnem T, Al-Saad S, Al-Shibli K, Andersen S, Sirera R, et al. The role of tumor stroma in cancer progression and prognosis: emphasis on carcinoma-associated fibroblasts and non-small cell lung cancer. J Thorac Oncol. (2011) 6:209–17. doi: 10.1097/JTO.0b013e3181f8a1bd

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Freeman MR, Li Q, Chung LW. Can stroma reaction predict cancer lethality? Clin Cancer Res. (2013) 19:4905–7. doi: 10.1158/1078-0432.CCR-13-1694

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Ueno H, Ishiguro M, Nakatani E, Ishikawa T, Uetake H, Murotani K, et al. Prognostic value of desmoplastic reaction characterisation in stage II colon cancer: prospective validation in a Phase 3 study (SACURA Trial). Br J Cancer. (2021) 124:1088–97. doi: 10.1038/s41416-020-01222-8

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Li HJ, Bera K, Toro P, Fu PF, Zhang ZL, Lu C, et al. Collagen fiber orientation disorder from H&E images is prognostic for early stage breast cancer: clinical trial validation. Npj Breast Cancer. (2021) 7:104. doi: 10.1038/s41523-021-00310-z

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Geessink OGF, Baidoshvili A, Klaase JM, Ehteshami Bejnordi B, Litjens GJS, van Pelt GW, et al. Computer aided quantification of intratumoral stroma yields an independent prognosticator in rectal cancer. Cell Oncol. (2019) 42:331–41. doi: 10.1007/s13402-019-00429-z

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Failmezger H, Muralidhar S, Rullan A, de Andrea CE, Sahai E, Yuan Y. Topological tumor graphs: a graph-based spatial model to infer stromal recruitment for immunosuppression in melanoma histology. Cancer Res. (2020) 80:1199–209. doi: 10.1158/0008-5472.CAN-19-2268

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Bankhead P, Loughrey MB, Fernández JA, Dombrowski Y, McArt DG, Dunne PD, et al. QuPath: open source software for digital pathology image analysis. Sci Rep. (2017) 7:1–7. doi: 10.1038/s41598-017-17204-5

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Jiang J, Tekin B, Guo R, Liu H, Huang Y, Wang C. Digital pathology-based study of cell-and tissue-level morphologic features in serous borderline ovarian tumor and high-grade serous ovarian cancer. J Pathol Inform. (2021) 12:24. doi: 10.4103/jpi.jpi_76_20

CrossRef Full Text | Google Scholar

19. He K, Gkioxari G, Dollár P, Girshick R. Mask r-cnn. Proceedings of the 2017 IEEE International Conference on Computer Vision (ICCV). Venice: (2017). p. 2961–9. doi: 10.1109/ICCV.2017.322

CrossRef Full Text | Google Scholar

20. Wang S, Rong R, Yang DM, Fujimoto J, Yan S, Cai L, et al. Computational staining of pathology images to study the tumor microenvironment in lung cancer. Cancer Res. (2020) 80:2056–66. doi: 10.1158/0008-5472.CAN-19-1629

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Mulay S, Ram K, Sivaprakasam M, Vinekar A. Early detection of retinopathy of prematurity stage using deep learning approach. Proceedings of the SPIE 10950, Medical Imaging 2019: Computer-Aided Diagnosis. San Diego, CA: SPIE (2019). p. 758–64. doi: 10.1117/12.2512719

CrossRef Full Text | Google Scholar

22. Quoc TTP, Linh TT, Minh TNT. Comparing U-Net convolutional network with mask R-CNN in agricultural area segmentation on satellite images. Proceedings of the 2020 7th NAFOSTED Conference on Information and Computer Science (NICS). Ho Chi Minh: IEEE (2020). p. 124–9. doi: 10.1109/NICS51282.2020.9335856

CrossRef Full Text | Google Scholar

23. Durkee MS, Abraham R, Ai J, Fuhrman JD, Clark MR, Giger ML. Comparing Mask R-CNN and U-Net architectures for robust automatic segmentation of immune cells in immunofluorescence images of Lupus Nephritis biopsies. Proceedings of the Imaging, Manipulation, and Analysis of Biomolecules, Cells, and Tissue. Bellingham, WA: SPIE (2021). p. 109–15. doi: 10.1117/12.2577785

CrossRef Full Text | Google Scholar

24. Abdulla W. Mask R-CNN for Object Detection and Instance Segmentation on Keras and TensorFlow. San Francisco, CA: GitHub repository (2017).

Google Scholar

25. Jiang J, Wu Y, Huang M, Yang W, Chen W, Feng Q. 3D brain tumor segmentation in multimodal MR images based on learning population-and patient-specific feature sets. Comput Med Imaging Graph. (2013) 37:512–21. doi: 10.1016/j.compmedimag.2013.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zhou D, Fang J, Song X, Guan C, Yin J, Dai Y, et al. Iou loss for 2d/3d object detection. Proceedings of the 2019 International Conference on 3D Vision (3DV). Quebec City, QC: IEEE (2019). p. 85–94. doi: 10.1109/3DV.2019.00019

CrossRef Full Text | Google Scholar

27. Simonyan K, Zisserman A. Very deep convolutional networks for large-scale image recognition. arXiv [Preprint] (2014): doi: 10.48550/arXiv.1409.1556

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Saltz J, Gupta R, Hou L, Kurc T, Singh P, Nguyen V, et al. Spatial organization and molecular correlation of tumor-infiltrating lymphocytes using deep learning on pathology images. Cell Rep. (2018) 23:181.–193.

Google Scholar

29. Jiang J, Prodduturi N, Chen D, Gu Q, Flotte T, Feng Q, et al. Image-to-image translation for automatic ink removal in whole slide images. J Med Imaging. (2020) 7:057502. doi: 10.1117/1.JMI.7.5.057502

CrossRef Full Text | Google Scholar

30. Gil JY, Kimmel R. Efficient dilation, erosion, opening, and closing algorithms. IEEE Trans Pattern Analy Mach Intellig. (2002) 24:1606–17. doi: 10.1109/TPAMI.2002.1114852

CrossRef Full Text | Google Scholar

31. Dundar MM, Badve S, Bilgin G, Raykar V, Jain R, Sertel O, et al. Computerized classification of intraductal breast lesions using histopathological images. IEEE Trans Biomed Eng. (2011) 58:1977–84. doi: 10.1109/TBME.2011.2110648

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protocols. (2009) 4:44–57. doi: 10.1038/nprot.2008.211

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. (2009) 37:1–13. doi: 10.1093/nar/gkn923

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Bredfeldt JS, Liu Y, Conklin MW, Keely PJ, Mackie TR, Eliceiri KW. Automated quantification of aligned collagen for human breast carcinoma prognosis. J Pathol Inform. (2014) 5:28. doi: 10.4103/2153-3539.139707

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Fei T, Zhang T, Shi W, Yu T. Mitigating the adverse impact of batch effects in sample pattern detection. Bioinformatics. (2018) 34:2634–41. doi: 10.1093/bioinformatics/bty117

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Schmitt M, Maron RC, Hekler A, Stenzinger A, Hauschild A, Weichenthal M, et al. Hidden variables in deep learning digital pathology and their potential to cause batch effects: prediction model study. J Med Intern Res. (2021) 23:e23436. doi: 10.2196/23436

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zheng Y, Jiang Z, Zhang H, Xie F, Shi J, Xue C. Adaptive color deconvolution for histological WSI normalization. Comput Methods Prog Biomed. (2019) 170:107–20. doi: 10.1016/j.cmpb.2019.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Boschman J, Farahani H, Darbandsari A, Ahmadvand P, Van Spankeren A, Farnell D, et al. The utility of color normalization for AI-based diagnosis of hematoxylin and eosin-stained pathology images. J Pathol. (2022) 256:15–24. doi: 10.1002/path.5797

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Van Eycke Y-R, Allard J, Salmon I, Debeir O, Decaestecker C. Image processing in digital pathology: an opportunity to solve inter-batch variability of immunohistochemical staining. Sci Rep. (2017) 7:42964. doi: 10.1038/srep42964

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Niazi MKK, Parwani AV, Gurcan MN. Digital pathology and artificial intelligence. Lancet Oncol. (2019) 20:e253–61. doi: 10.1016/S1470-2045(19)30154-8

CrossRef Full Text | Google Scholar

41. Hou L, Agarwal A, Samaras D, Kurc TM, Gupta RR, Saltz JH. Robust histopathology image analysis: to label or to synthesize? Proceedings of the 2019 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR). Long Beach, CA: (2019). p. 8533–42. doi: 10.1109/CVPR.2019.00873

PubMed Abstract | CrossRef Full Text | Google Scholar

42. de Ruiter EJ, de Roest RH, Brakenhoff RH, Leemans CR, de Bree R, Terhaard CH, et al. Digital pathology-aided assessment of tumor-infiltrating T lymphocytes in advanced stage, HPV-negative head and neck tumors. Cancer Immunol Immunother. (2020) 69:581–91. doi: 10.1007/s00262-020-02481-3

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Weigert M, Schmidt U, Haase R, Sugawara K, Myers G “Star-convex polyhedra for 3d object detection and segmentation in microscopy,” in Proceedings of the IEEE/CVF Winter Conference on Applications of Computer Vision (Los Alamitos, CA: IEEE Computer Society) (2020). 3666–73. doi: 10.1109/WACV45572.2020.9093435

CrossRef Full Text | Google Scholar

44. Jiang J, Tekin B, Guo R, Liu H, Huang Y, Wang C. Digital pathology-based study of cell-and tissue-level morphologic features in serous borderline ovarian tumor and high-grade serous ovarian cancer. arXiv [Preprint] (2020). doi: 10.48550/arXiv.2008.12479

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Van Rijthoven M, Balkenhol M, Siliòa K, Van Der Laak J, Ciompi F. HookNet: multi-resolution convolutional neural networks for semantic segmentation in histopathology whole-slide images. Med Image Analy. (2021) 68:101890. doi: 10.1016/j.media.2020.101890

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Marini N, Otálora S, Ciompi F, Silvello G, Marchesin S, Vatrano S, et al. Multi-Scale task multiple instance learning for the classification of digital pathology images with global annotations. Proc Mach Learn Res. (2021) 156:170–81.

Google Scholar

47. Sornapudi S, Stanley RJ, Stoecker WV, Long R, Xue Z, Zuna R, et al. DeepCIN: attention-based cervical histology image classification with sequential feature modeling for pathologist-level accuracy. J Pathol Inform. (2020) 11:40. doi: 10.4103/jpi.jpi_50_20

CrossRef Full Text | Google Scholar

Keywords: tumor-stroma reaction, high-grade serous ovarian carcinoma, digital pathology, prognostic fibrosis, molecular signature

Citation: Jiang J, Tekin B, Yuan L, Armasu S, Winham SJ, Goode EL, Liu H, Huang Y, Guo R and Wang C (2022) Computational tumor stroma reaction evaluation led to novel prognosis-associated fibrosis and molecular signature discoveries in high-grade serous ovarian carcinoma. Front. Med. 9:994467. doi: 10.3389/fmed.2022.994467

Received: 14 July 2022; Accepted: 15 August 2022;
Published: 07 September 2022.

Edited by:

Jun Cheng, Shenzhen University, China

Reviewed by:

Shidan Wang, University of Texas Southwestern Medical Center, United States
Yao Wu, Children’s National Hospital, United States

Copyright © 2022 Jiang, Tekin, Yuan, Armasu, Winham, Goode, Liu, Huang, Guo and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chen Wang, wang.chen@mayo.edu; Hongfang Liu, liu.hongfang@mayo.edu; Yajue Huang, huang.yajue@mayo.edu; Ruifeng Guo, guo.ruifeng@mayo.edu

These authors have contributed equally to this work and share first authorship

These authors share last authorship

Download