Transcriptional Profiling of a Patient-Matched Cohort of Glioblastoma (IDH-Wildtype) for Therapeutic Target and Repurposing Drug Identification.

Glioblastoma (GBM) is the most prevalent and aggressive adult brain tumor. Despite multi-modal therapies, GBM recurs, and patients have poor survival (~14 months). Resistance to therapy may originate from a subpopulation of tumor cells identified as glioma-stem cells (GSC), and new treatments are urgently needed to target these. The biology underpinning GBM recurrence was investigated using whole transcriptome profiling of patient-matched initial and recurrent GBM (recGBM). Differential expression analysis identified 147 significant probes. In total, 24 genes were validated using expression data from four public cohorts and the literature. Functional analyses revealed that transcriptional changes to recGBM were dominated by angiogenesis and immune-related processes. The role of MHC class II proteins in antigen presentation and the differentiation, proliferation, and infiltration of immune cells was enriched. These results suggest recGBM would benefit from immunotherapies. The altered gene signature was further analyzed in a connectivity mapping analysis with QUADrATiC software to identify FDA-approved repurposing drugs. Top-ranking target compounds that may be effective against GSC and GBM recurrence were rosiglitazone, nizatidine, pantoprazole, and tolmetin. Our translational bioinformatics pipeline provides an approach to identify target compounds for repurposing that may add clinical benefit in addition to standard therapies against resistant cancers such as GBM.

Particle swarm optimization artificial intelligence technique for gene signature discovery in transcriptomic cohorts.

The development of gene signatures is key for delivering personalized medicine, despite only a few signatures being available for use in the clinic for cancer patients. Gene signature discovery tends to revolve around identifying a single signature. However, it has been shown that various highly predictive signatures can be produced from the same dataset. This study assumes that the presentation of top ranked signatures will allow greater efforts in the selection of gene signatures for validation on external datasets and for their clinical translation. Particle swarm optimization (PSO) is an evolutionary algorithm often used as a search strategy and largely represented as binary PSO (BPSO) in this domain. BPSO, however, fails to produce succinct feature sets for complex optimization problems, thus affecting its overall runtime and optimization performance. Enhanced BPSO (EBPSO) was developed to overcome these shortcomings. Thus, this study will validate unique candidate gene signatures for different underlying biology from EBPSO on transcriptomics cohorts. EBPSO was consistently seen to be as accurate as BPSO with substantially smaller feature signatures and significantly faster runtimes. 100% accuracy was achieved in all but two of the selected data sets. Using clinical transcriptomics cohorts, EBPSO has demonstrated the ability to identify accurate, succinct, and significantly prognostic signatures that are unique from one another. This has been proposed as a promising alternative to overcome the issues regarding traditional single gene signature generation. Interpretation of key genes within the signatures provided biological insights into the associated functions that were well correlated to their cancer type.

A bespoke target selection tool to guide biomarker discovery in tubo-ovarian cancer.

Introduction: Cancers presenting at advanced stages inherently have poor prognosis. High grade serous carcinoma (HGSC) is the most common and aggressive form of tubo-ovarian cancer. Clinical tests to accurately diagnose and monitor this condition are lacking. Hence, development of disease-specific tests are urgently required. Methods: The molecular profile of HGSC during disease progression was investigated in a unique patient cohort. A bespoke data browser was developed to analyse gene expression and DNA methylation datasets for biomarker discovery. The Ovarian Cancer Data Browser (OCDB) is built in C# with a.NET framework using an integrated development environment of Microsoft Visual Studio and fast access files (.faf). The graphical user interface is easy to navigate between four analytical modes (gene expression; methylation; combined gene expression and methylation data; methylation clusters), with a rapid query response time. A user should first define a disease progression trend for prioritising results. Single or multiomics data are then mined to identify probes, genes and methylation clusters that exhibit the desired trend. A unique scoring system based on the percentage change in expression/methylation between disease stages is used. Results are filtered and ranked using weighting and penalties. Results: The OCDB's utility for biomarker discovery is demonstrated with the identified target OSR2. Trends in OSR2 repression and hypermethylation with HGSC disease progression were confirmed in the browser samples and an independent cohort using bioassays. The OSR2 methylation biomarker could discriminate HGSC with high specificity (95%) and sensitivity (93.18%). Conclusions: The OCDB has been refined and validated to be an integral part of a unique biomarker discovery pipeline. It may also be used independently to aid identification of novel targets. It carries the potential to identify further biomarker assays that can reduce type I and II errors within clinical diagnostics.

Novel deep learning-based solution for identification of prognostic subgroups in liver cancer (Hepatocellular carcinoma).

BACKGROUND: Liver cancer (Hepatocellular carcinoma; HCC) prevalence is increasing and with poor clinical outcome expected it means greater understanding of HCC aetiology is urgently required. This study explored a deep learning solution to detect biologically important features that distinguish prognostic subgroups. A novel architecture of an Artificial Neural Network (ANN) trained with a customised objective function (LRSC) was developed. The ANN should discover new data representations, to detect patient subgroups that are biologically homogenous (clustering loss) and similar in survival (survival loss) while removing noise from the data (reconstruction loss). The model was applied to TCGA-HCC multi-omics data and benchmarked against baseline models that only use a reconstruction objective function (BCE, MSE) for learning. With the baseline models, the new features are then filtered based on survival information and used for clustering patients. Different variants of the customised objective function, incorporating only reconstruction and clustering losses (LRC); and reconstruction and survival losses (LRS) were also evaluated. Robust features consistently detected were compared between models and validated in TCGA and LIRI-JP HCC cohorts. RESULTS: The combined loss (LRSC) discovered highly significant prognostic subgroups (P-value = 1.55E-77) with more accurate sample assignment (Silhouette scores: 0.59-0.7) compared to baseline models (0.18-0.3). All LRSC bottleneck features (N = 100) were significant for survival, compared to only 11-21 for baseline models. Prognostic subgroups were not explained by disease grade or risk factors. Instead LRSC identified robust features including 377 mRNAs, many of which were novel (61.27%) compared to those identified by the other losses. Some 75 mRNAs were prognostic in TCGA, while 29 were prognostic in LIRI-JP also. LRSC also identified 15 robust miRNAs including two novel (hsa-let-7g; hsa-mir-550a-1) and 328 methylation features with 71% being prognostic. Gene-enrichment and Functional Annotation Analysis identified seven pathways differentiating prognostic clusters. CONCLUSIONS: Combining cluster and survival metrics with the reconstruction objective function facilitated superior prognostic subgroup identification. The hybrid model identified more homogeneous clusters that consequently were more biologically meaningful. The novel and prognostic robust features extracted provide additional information to improve our understanding of a complex disease to help reveal its aetiology. Moreover, the gene features identified may have clinical applications as therapeutic targets.

Repurposing FDA approved drugs as radiosensitizers for treating hypoxic prostate cancer.

BACKGROUND: The presence of hypoxia is a poor prognostic factor in prostate cancer and the hypoxic tumor microenvironment promotes radioresistance. There is potential for drug radiotherapy combinations to improve the therapeutic ratio. We aimed to investigate whether hypoxia-associated genes could be used to identify FDA approved drugs for repurposing for the treatment of hypoxic prostate cancer. METHODS: Hypoxia associated genes were identified and used in the connectivity mapping software QUADrATIC to identify FDA approved drugs as candidates for repurposing. Drugs identified were tested in vitro in prostate cancer cell lines (DU145, PC3, LNCAP). Cytotoxicity was investigated using the sulforhodamine B assay and radiosensitization using a clonogenic assay in normoxia and hypoxia. RESULTS: Menadione and gemcitabine had similar cytotoxicity in normoxia and hypoxia in all three cell lines. In DU145 cells, the radiation sensitizer enhancement ratio (SER) of menadione was 1.02 in normoxia and 1.15 in hypoxia. The SER of gemcitabine was 1.27 in normoxia and 1.09 in hypoxia. No radiosensitization was seen in PC3 cells. CONCLUSION: Connectivity mapping can identify FDA approved drugs for potential repurposing that are linked to a radiobiologically relevant phenotype. Gemcitabine and menadione could be further investigated as potential radiosensitizers in prostate cancer.

Evolutionary genetic algorithm identifies IL2RB as a potential predictive biomarker for immune-checkpoint therapy in colorectal cancer.

Identifying robust predictive biomarkers to stratify colorectal cancer (CRC) patients based on their response to immune-checkpoint therapy is an area of unmet clinical need. Our evolutionary algorithm Atlas Correlation Explorer (ACE) represents a novel approach for mining The Cancer Genome Atlas (TCGA) data for clinically relevant associations. We deployed ACE to identify candidate predictive biomarkers of response to immune-checkpoint therapy in CRC. We interrogated the colon adenocarcinoma (COAD) gene expression data across nine immune-checkpoints (PDL1, PDCD1, CTLA4, LAG3, TIM3, TIGIT, ICOS, IDO1 and BTLA). IL2RB was identified as the most common gene associated with immune-checkpoint genes in CRC. Using human/murine single-cell RNA-seq data, we demonstrated that IL2RB was expressed predominantly in a subset of T-cells associated with increased immune-checkpoint expression (P < 0.0001). Confirmatory IL2RB immunohistochemistry (IHC) analysis in a large MSI-H colon cancer tissue microarray (TMA; n = 115) revealed sensitive, specific staining of a subset of lymphocytes and a strong association with FOXP3+ lymphocytes (P < 0.0001). IL2RB mRNA positively correlated with three previously-published gene signatures of response to immune-checkpoint therapy (P < 0.0001). Our evolutionary algorithm has identified IL2RB to be extensively linked to immune-checkpoints in CRC; its expression should be investigated for clinical utility as a potential predictive biomarker for CRC patients receiving immune-checkpoint blockade.

Prostate cancer heterogeneity assessment with multi-regional sampling and alignment-free methods.

Combining alignment-free methods for phylogenetic analysis with multi-regional sampling using next-generation sequencing can provide an assessment of intra-patient tumour heterogeneity. From multi-regional sampling divergent branching, we validated two different lesions within a patient's prostate. Where multi-regional sampling has not been used, a single sample from one of these areas could misguide as to which drugs or therapies would best benefit this patient, due to the fact these tumours appear to be genetically different. This application has the power to render, in a fraction of the time used by other approaches, intra-patient heterogeneity and decipher aberrant biomarkers. Another alignment-free method for calling single-nucleotide variants from raw next-generation sequencing samples has determined possible variants and genomic locations that may be able to characterize the differences between the two main branching patterns. Alignment-free approaches have been applied to relevant clinical multi-regional samples and may be considered as a valuable option for comparing and determining heterogeneity to help deliver personalized medicine through more robust efforts in identifying targetable pathways and therapeutic strategies. Our study highlights the application these tools could have on patient-aligned treatment indications.

Immune status is prognostic for poor survival in colorectal cancer patients and is associated with tumour hypoxia.

BACKGROUND: Immunohistochemical quantification of the immune response is prognostic for colorectal cancer (CRC). Here, we evaluate the suitability of alternative immune classifiers on prognosis and assess whether they relate to biological features amenable to targeted therapy. METHODS: Overall survival by immune (CD3, CD4, CD8, CD20 and FOXP3) and immune-checkpoint (ICOS, IDO-1 and PD-L1) biomarkers in independent CRC cohorts was evaluated. Matched mutational and transcriptomic data were interrogated to identify associated biology. RESULTS: Determination of immune-cold tumours by combined low-density cell counts of CD3, CD4 and CD8 immunohistochemistry constituted the best prognosticator across stage II-IV CRC, particularly in patients with stage IV disease (HR 1.98 [95% CI: 1.47-2.67]). These immune-cold CRCs were associated with tumour hypoxia, confirmed using CAIX immunohistochemistry (P = 0.0009), which may mediate disease progression through common biology (KRAS mutations, CRIS-B subtype and SPP1 mRNA overexpression). CONCLUSIONS: Given the significantly poorer survival of immune-cold CRC patients, these data illustrate that assessment of CD4-expressing cells complements low CD3 and CD8 immunohistochemical quantification in the tumour bulk, potentially facilitating immunophenotyping of patient biopsies to predict prognosis. In addition, we found immune-cold CRCs to associate with a difficult-to-treat, poor prognosis hypoxia signature, indicating that these patients may benefit from hypoxia-targeting clinical trials.

Defining the molecular evolution of extrauterine high grade serous carcinoma.

OBJECTIVE: High grade serous carcinoma (HGSC) is the most common and most aggressive, subtype of epithelial ovarian cancer. It presents as advanced stage disease with poor prognosis. Recent pathological evidence strongly suggests HGSC arises from the fallopian tube via the precursor lesion; serous tubal intraepithelial carcinoma (STIC). However, further definition of the molecular evolution of HGSC has major implications for both clinical management and research. This study aims to more clearly define the molecular pathogenesis of HGSC. METHODS: Six cases of HGSC were identified at the Northern Ireland Gynaecological Cancer Centre (NIGCC) that each contained ovarian HGSC (HGSC), omental HGSC (OMT), STIC, normal fallopian tube epithelium (FTE) and normal ovarian surface epithelium (OSE). The relevant formalin-fixed paraffin embedded (FFPE) tissue samples were retrieved from the pathology archive via the Northern Ireland Biobank following attaining ethical approval (NIB11:005). Full microarray-based gene expression profiling was performed on the cohort. The resulting data was analysed bioinformatically and the results were validated in a HGSC-specific in-vitro model. RESULTS: The carcinogenesis of HGSC was investigated and showed the molecular profile of HGSC to be more closely related to normal FTE than OSE. STIC lesions also clustered closely with HGSC, indicating a common molecular origin. CONCLUSION: This study provides strong evidence suggesting that extrauterine HGSC arises from the fimbria of the distal fallopian tube. Furthermore, several potential pathways were identified which could be targeted by novel therapies for HGSC. These findings have significant translational relevance for both primary prevention and clinical management of the disease.

NUQA: Estimating Cancer Spatial and Temporal Heterogeneity and Evolution through Alignment-Free Methods.

Longitudinal next-generation sequencing of cancer patient samples has enhanced our understanding of the evolution and progression of various cancers. As a result, and due to our increasing knowledge of heterogeneity, such sampling is becoming increasingly common in research and clinical trial sample collections. Traditionally, the evolutionary analysis of these cohorts involves the use of an aligner followed by subsequent stringent downstream analyses. However, this can lead to large levels of information loss due to the vast mutational landscape that characterizes tumor samples. Here, we propose an alignment-free approach for sequence comparison-a well-established approach in a range of biological applications including typical phylogenetic classification. Such methods could be used to compare information collated in raw sequence files to allow an unsupervised assessment of the evolutionary trajectory of patient genomic profiles. In order to highlight this utility in cancer research we have applied our alignment-free approach using a previously established metric, Jensen-Shannon divergence, and a metric novel to this area, Hellinger distance, to two longitudinal cancer patient cohorts in glioma and clear cell renal cell carcinoma using our software, NUQA. We hypothesize that this approach has the potential to reveal novel information about the heterogeneity and evolutionary trajectory of spatiotemporal tumor samples, potentially revealing early events in tumorigenesis and the origins of metastases and recurrences. Key words: alignment-free, Hellinger distance, exome-seq, evolution, phylogenetics, longitudinal.

A Novel Role for Cathepsin S as a Potential Biomarker in Triple Negative Breast Cancer.

Cathepsin S (CTSS) has previously been implicated in a number of cancer types, where it is associated with poor clinical features and outcome. To date, patient outcome in breast cancer has not been examined with respect to this protease. Here, we carried out immunohistochemical (IHC) staining of CTSS using a breast cancer tissue microarray in patients who received adjuvant therapy. We scored CTSS expression in the epithelial and stromal compartments and evaluated the association of CTSS expression with matched clinical outcome data. We observed differences in outcome based on CTSS expression, with stromal-derived CTSS expression correlating with a poor outcome and epithelial CTSS expression associated with an improved outcome. Further subtype characterisation revealed high epithelial CTSS expression in TNBC patients with improved outcome, which remained consistent across two independent TMA cohorts. Further in silico gene expression analysis, using both in-house and publicly available datasets, confirmed these observations and suggested high CTSS expression may also be beneficial to outcome in ER-/HER2+ cancer. Furthermore, high CTSS expression was associated with the BL1 Lehmann subgroup, which is characterised by defects in DNA damage repair pathways and correlates with improved outcome. Finally, analysis of matching IHC analysis reveals an increased M1 (tumour destructive) polarisation in macrophage in patients exhibiting high epithelial CTSS expression. In conclusion, our observations suggest epithelial CTSS expression may be prognostic of improved outcome in TNBC. Improved outcome observed with HER2+ at the gene expression level furthermore suggests CTSS may be prognostic of improved outcome in ER- cancers as a whole. Lastly, from the context of these patients receiving adjuvant therapy and as a result of its association with BL1 subgroup CTSS may be elevated in patients with defects in DNA damage repair pathways, indicating it may be predictive of tumour sensitivity to DNA damaging agents.

ACE: A Workbench Using Evolutionary Genetic Algorithms for Analyzing Association in TCGA.

Modern methods of acquiring molecular data have improved rapidly in recent years, making it easier for researchers to collect large volumes of information. However, this has increased the challenge of recognizing interesting patterns within the data. Atlas Correlation Explorer (ACE) is a user-friendly workbench for seeking associations between attributes in The Cancer Genome Atlas (TCGA) database. It allows any combination of clinical and genomic data streams to be searched using an evolutionary algorithm approach. To showcase ACE, we assessed which RNA sequencing transcripts were associated with estrogen receptor (ESR1) in the TCGA breast cancer cohort. The analysis revealed already well-established associations with XBP1 and FOXA1, but also identified a strong association with CT62, a potential immunotherapeutic target with few previous associations with breast cancer. In conclusion, ACE can produce results for very large searches in a short time and will serve as an increasingly useful tool for biomarker discovery in the big data era. SIGNIFICANCE: ACE uses an evolutionary algorithm approach to perform large searches for associations between any combinations of data in the TCGA database.

Impact of Variable RNA-Sequencing Depth on Gene Expression Signatures and Target Compound Robustness: Case Study Examining Brain Tumor (Glioma) Disease Progression.

PURPOSE: Gene expression profiling can uncover biologic mechanisms underlying disease and is important in drug development. RNA sequencing (RNA-seq) is routinely used to assess gene expression, but costs remain high. Sample multiplexing reduces RNAseq costs; however, multiplexed samples have lower cDNA sequencing depth, which can hinder accurate differential gene expression detection. The impact of sequencing depth alteration on RNA-seq-based downstream analyses such as gene expression connectivity mapping is not known, where this method is used to identify potential therapeutic compounds for repurposing. METHODS: In this study, published RNA-seq profiles from patients with brain tumor (glioma) were assembled into two disease progression gene signature contrasts for astrocytoma. Available treatments for glioma have limited effectiveness, rendering this a disease of poor clinical outcome. Gene signatures were subsampled to simulate sequencing alterations and analyzed in connectivity mapping to investigate target compound robustness. RESULTS: Data loss to gene signatures led to the loss, gain, and consistent identification of significant connections. The most accurate gene signature contrast with consistent patient gene expression profiles was more resilient to data loss and identified robust target compounds. Target compounds lost included candidate compounds of potential clinical utility in glioma (eg, suramin, dasatinib). Lost connections may have been linked to low-abundance genes in the gene signature that closely characterized the disease phenotype. Consistently identified connections may have been related to highly expressed abundant genes that were ever-present in gene signatures, despite data reductions. Potential noise surrounding findings included false-positive connections that were gained as a result of gene signature modification with data loss. CONCLUSION: Findings highlight the necessity for gene signature accuracy for connectivity mapping, which should improve the clinical utility of future target compound discoveries.

Bcl-xL as a poor prognostic biomarker and predictor of response to adjuvant chemotherapy specifically in BRAF-mutant stage II and III colon cancer.

Purpose: BRAF mutation occurs in 8-15% of colon cancers (CC), and is associated with poor prognosis in metastatic disease. Compared to wild-type BRAF (BRAFWT) disease, stage II/III CC patients with BRAF mutant (BRAFMT) tumors have shorter overall survival after relapse; however, time-to-relapse is not significantly different. The aim of this investigation was to identify, and validate, novel predictors of relapse of stage II/III BRAFMT CC. Experimental design: We used gene expression data from a cohort of 460 patients (GSE39582) to perform a supervised classification analysis based on risk-of-relapse within BRAFMT stage II/III CC, to identify transcriptomic biomarkers associated with prognosis within this genotype. These findings were validated using immunohistochemistry in an independent population-based cohort of Stage II/III CC (n = 691), applying Cox proportional hazards analysis to determine associations with survival. Results: High gene expression levels of Bcl-xL, a key regulator of apoptosis, were associated with increased risk of relapse, specifically in BRAFMT tumors (HR = 8.3, 95% CI 1.7-41.7), but not KRASMT/BRAFWT or KRASWT/BRAFWT tumors. High Bcl-xL protein expression in BRAFMT, untreated, stage II/III CC was confirmed to be associated with an increased risk of death in an independent cohort (HR = 12.13, 95% CI 2.49-59.13). Additionally, BRAFMT tumors with high levels of Bcl-xL protein expression appeared to benefit from adjuvant chemotherapy (P for interaction = 0.006), indicating the potential predictive value of Bcl-xL expression in this setting. Conclusions: These findings provide evidence that Bcl-xL gene and/or protein expression identifies a poor prognostic subgroup of BRAFMT stage II/III CC patients, who may benefit from adjuvant chemotherapy.

Depletion of DNMT1 in differentiated human cells highlights key classes of sensitive genes and an interplay with polycomb repression.

BACKGROUND: DNA methylation plays a vital role in the cell, but loss-of-function mutations of the maintenance methyltransferase DNMT1 in normal human cells are lethal, precluding target identification, and existing hypomorphic lines are tumour cells. We generated instead a hypomorphic series in normal hTERT-immortalised fibroblasts using stably integrated short hairpin RNA. RESULTS: Approximately two-thirds of sites showed demethylation as expected, with one-third showing hypermethylation, and targets were shared between the three independently derived lines. Enrichment analysis indicated significant losses at promoters and gene bodies with four gene classes most affected: (1) protocadherins, which are key to neural cell identity; (2) genes involved in fat homoeostasis/body mass determination; (3) olfactory receptors and (4) cancer/testis antigen (CTA) genes. Overall effects on transcription were relatively small in these fibroblasts, but CTA genes showed robust derepression. Comparison with siRNA-treated cells indicated that shRNA lines show substantial remethylation over time. Regions showing persistent hypomethylation in the shRNA lines were associated with polycomb repression and were derepressed on addition of an EZH2 inhibitor. Persistent hypermethylation in shRNA lines was, in contrast, associated with poised promoters. CONCLUSIONS: We have assessed for the first time the effects of chronic depletion of DNMT1 in an untransformed, differentiated human cell type. Our results suggest polycomb marking blocks remethylation and indicate the sensitivity of key neural, adipose and cancer-associated genes to loss of maintenance methylation activity.

Prospective patient stratification into robust cancer-cell intrinsic subtypes from colorectal cancer biopsies.

Colorectal cancer (CRC) biopsies underpin accurate diagnosis, but are also relevant for patient stratification in molecularly-guided clinical trials. The consensus molecular subtypes (CMSs) and colorectal cancer intrinsic subtypes (CRISs) transcriptional signatures have potential clinical utility for improving prognostic/predictive patient assignment. However, their ability to provide robust classification, particularly in pretreatment biopsies from multiple regions or at different time points, remains untested. In this study, we undertook a comprehensive assessment of the robustness of CRC transcriptional signatures, including CRIS and CMS, using a range of tumour sampling methodologies currently employed in clinical and translational research. These include analyses using (i) laser-capture microdissected CRC tissue, (ii) eight publically available rectal cancer biopsy data sets (n = 543), (iii) serial biopsies (from AXEBeam trial, NCT00828672; n = 10), (iv) multi-regional biopsies from colon tumours (n = 29 biopsies, n = 7 tumours), and (v) pretreatment biopsies from the phase II rectal cancer trial COPERNCIUS (NCT01263171; n = 44). Compared to previous results obtained using CRC resection material, we demonstrate that CMS classification in biopsy tissue is significantly less capable of reliably classifying patient subtype (43% unknown in biopsy versus 13% unknown in resections, p = 0.0001). In contrast, there was no significant difference in classification rate between biopsies and resections when using the CRIS classifier. Additionally, we demonstrated that CRIS provides significantly better spatially- and temporally- robust classification of molecular subtypes in CRC primary tumour tissue compared to CMS (p = 0.003 and p = 0.02, respectively). These findings have potential to inform ongoing biopsy-based patient stratification in CRC, enabling robust and stable assignment of patients into clinically-informative arms of prospective multi-arm, multi-stage clinical trials. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Integrated tumor identification and automated scoring minimizes pathologist involvement and provides new insights to key biomarkers in breast cancer.

Digital image analysis (DIA) is becoming central to the quantitative evaluation of tissue biomarkers for discovery, diagnosis and therapeutic selection for the delivery of precision medicine. In this study, automated DIA using a new purpose-built software platform (QuPath) is applied to a cohort of 293 breast cancer patients to score five biomarkers in tissue microarrays (TMAs): ER, PR, HER2, Ki67 and p53. This software is able to measure IHC expression following fully automated tumor recognition in the same immunohistochemical (IHC)-stained tissue section, as part of a rapid workflow to ensure objectivity and accelerate biomarker analysis. The digital scores produced by QuPath were compared with manual scores by a pathologist and shown to have a good level of concordance in all cases (Cohen's κ>0.6), and almost perfect agreement for the clinically relevant biomarkers ER, PR and HER2 (κ>0.86). To assess prognostic value, cutoff thresholds could be applied to both manual and automated scores using the QuPath software, and survival analysis performed for 5-year overall survival. DIA was shown to be capable of replicating the statistically significant stratification of patients achieved using manual scoring across all biomarkers (P<0.01, log-rank test). Furthermore, the image analysis scores were shown to consistently lead to statistical significance across a wide range of potential cutoff thresholds, indicating the robustness of the method, and identify sub-populations of cases exhibiting different expression patterns within the p53 and Ki67 data sets that warrant further investigation. These findings have demonstrated QuPath's suitability for fast, reproducible, high-throughput TMA analysis across a range of important biomarkers. This was achieved using our tumor recognition algorithms for IHC-stained sections, trained interactively without the need for any additional tumor recognition markers, for example, cytokeratin, to obtain greater insight into the relationship between biomarker expression and clinical outcome applicable to a range of cancer types.

Automated Tumour Recognition and Digital Pathology Scoring Unravels New Role for PD-L1 in Predicting Good Outcome in ER-/HER2+ Breast Cancer.

The role of PD-L1 as a prognostic and predictive biomarker is an area of great interest. However, there is a lack of consensus on how to deliver PD-L1 as a clinical biomarker. At the heart of this conundrum is the subjective scoring of PD-L1 IHC in most studies to date. Current standard scoring systems involve separation of epithelial and inflammatory cells and find clinical significance in different percentages of expression, e.g., above or below 1%. Clearly, an objective, reproducible and accurate approach to PD-L1 scoring would bring a degree of necessary consistency to this landscape. Using a systematic comparison of technologies and the application of QuPath, a digital pathology platform, we show that high PD-L1 expression is associated with improved clinical outcome in Triple Negative breast cancer in the context of standard of care (SoC) chemotherapy, consistent with previous findings. In addition, we demonstrate for the first time that high PD-L1 expression is also associated with better outcome in ER- disease as a whole including HER2+ breast cancer. We demonstrate the influence of antibody choice on quantification and clinical impact with the Ventana antibody (SP142) providing the most robust assay in our hands. Through sampling different regions of the tumour, we show that tumour rich regions display the greatest range of PD-L1 expression and this has the most clinical significance compared to stroma and lymphoid rich areas. Furthermore, we observe that both inflammatory and epithelial PD-L1 expression are associated with improved survival in the context of chemotherapy. Moreover, as seen with PD-L1 inhibitor studies, a low threshold of PD-L1 expression stratifies patient outcome. This emphasises the importance of using digital pathology and precise biomarker quantitation to achieve accurate and reproducible scores that can discriminate low PD-L1 expression.

QuPath: Open source software for digital pathology image analysis.

QuPath is new bioimage analysis software designed to meet the growing need for a user-friendly, extensible, open-source solution for digital pathology and whole slide image analysis. In addition to offering a comprehensive panel of tumor identification and high-throughput biomarker evaluation tools, QuPath provides researchers with powerful batch-processing and scripting functionality, and an extensible platform with which to develop and share new algorithms to analyze complex tissue images. Furthermore, QuPath's flexible design makes it suitable for a wide range of additional image analysis applications across biomedical research.

IHC-based subcellular quantification provides new insights into prognostic relevance of FLIP and procaspase-8 in non-small-cell lung cancer.

In this study, we developed an image analysis algorithm for quantification of two potential apoptotic biomarkers in non-small-cell lung cancer (NSCLC): FLIP and procaspase-8. Immunohistochemical expression of FLIP and procaspase-8 in 184 NSCLC tumors were assessed. Individual patient cores were segmented and classified as tumor and stroma using the Definiens Tissue Studio. Subsequently, chromogenic expression of each biomarker was measured separately in the nucleus and cytoplasm and reported as a quantitative histological score. The software package pROC was applied to define biomarker thresholds. Cox proportional hazards analysis was applied to generate hazard ratios (HRs) and associated 95% CI for survival. High cytoplasmic expression of tumoral (but not stromal) FLIP was associated with a 2.5-fold increased risk of death in lung adenocarcinoma patients, even when adjusted for known confounders (HR 2.47, 95% CI 1.14-5.35). Neither nuclear nor cytoplasmic tumoral procaspase-8 expression was associated with overall survival in lung adenocarcinoma patients; however, there was a significant trend (P for trend=0.03) for patients with adenocarcinomas with both high cytoplasmic FLIP and high cytoplasmic procaspase-8 to have a multiplicative increased risk of death. Notably, high stromal nuclear procaspase-8 expression was associated with a reduced risk of death in lung adenocarcinoma patients (adjusted HR 0.31, 95% CI 0.15-0.66). On further examination, the cells with high nuclear procaspase-8 were found to be of lymphoid origin, suggesting that the better prognosis of patients with tumors with high stromal nuclear procaspase-8 is related to immune infiltration, a known favorable prognostic factor. No significant associations were detected in analysis of lung squamous cell carcinoma patients. Our results suggest that cytoplasmic expression of FLIP in the tumor and nuclear expression of procaspase-8 in the stroma are prognostically relevant in non-small-cell adenocarcinomas but not in squamous cell carcinomas of the lung.