surviveR: a flexible shiny application for patient survival analysis

Kaplan-Meier (KM) survival analyses based on complex patient categorization due to the burgeoning volumes of genomic, molecular and phenotypic data, are an increasingly important aspect of the biomedical researcher's toolkit. Commercial statistics and graphing packages for such analyses are functionally limited, whereas open-source tools have a high barrier-to-entry in terms of understanding of methodologies and computational expertise. We developed surviveR to address this unmet need for a survival analysis tool that can enable users with limited computational expertise to conduct routine but complex analyses. surviveR is a cloud-based Shiny application, that addresses our identified unmet need for an easy-to-use web-based tool that can plot and analyse survival based datasets. Integrated customization options allows a user with limited computational expertise to easily filter patients to enable custom cohort generation, automatically calculate log-rank test and Cox hazard ratios. Continuous datasets can be integrated, such as RNA or protein expression measurements which can be then used as categories for survival plotting. We further demonstrate the utility through exemplifying its application to a clinically relevant colorectal cancer patient dataset. surviveR is a cloud-based web application available at https://generatr.qub.ac.uk/app/surviveR , that can be used by non-experts users to perform complex custom survival analysis.

classifieR a flexible interactive cloud-application for functional annotation of cancer transcriptomes

BACKGROUND

Transcriptionally informed predictions are increasingly important for sub-typing cancer patients, understanding underlying biology and to inform novel treatment strategies. For instance, colorectal cancers (CRCs) can be classified into four CRC consensus molecular subgroups (CMS) or five intrinsic (CRIS) sub-types that have prognostic and predictive value. Breast cancer (BRCA) has five PAM50 molecular subgroups with similar value, and the OncotypeDX test provides transcriptomic based clinically actionable treatment-risk stratification. However, assigning samples to these subtypes and other transcriptionally inferred predictions is time consuming and requires significant bioinformatics experience. There is no "universal" method of using data from diverse assay/sequencing platforms to provide subgroup classification using the established classifier sets of genes (CMS, CRIS, PAM50, OncotypeDX), nor one which in provides additional useful functional annotations such as cellular composition, single-sample Gene Set Enrichment Analysis, or prediction of transcription factor activity.

RESULTS

To address this bottleneck, we developed classifieR, an easy-to-use R-Shiny based web application that supports flexible rapid single sample annotation of transcriptional profiles derived from cancer patient samples form diverse platforms. We demonstrate the utility of the " classifieR" framework to applications focused on the analysis of transcriptional profiles from colorectal (classifieRc) and breast (classifieRb). Samples are annotated with disease relevant transcriptional subgroups (CMS/CRIS sub-types in classifieRc and PAM50/inferred OncotypeDX in classifieRb), estimation of cellular composition using MCP-counter and xCell, single-sample Gene Set Enrichment Analysis (ssGSEA) and transcription factor activity predictions with Discriminant Regulon Expression Analysis (DoRothEA).

CONCLUSIONS

classifieR provides a framework which enables labs without access to a dedicated bioinformation can get information on the molecular makeup of their samples, providing an insight into patient prognosis, druggability and also as a tool for analysis and discovery. Applications are hosted online at https://generatr.qub.ac.uk/app/classifieRc and https://generatr.qub.ac.uk/app/classifieRb after signing up for an account on https://generatr.qub.ac.uk .

Clinical Positioning of the IAP Antagonist Tolinapant (ASTX660) in Colorectal Cancer

Inhibitors of apoptosis proteins (IAPs) are intracellular proteins, with important roles in regulating cell death, inflammation, and immunity. Here, we examined the clinical and therapeutic relevance of IAPs in colorectal cancer. We found that elevated expression of cIAP1 and cIAP2 (but not XIAP) significantly correlated with poor prognosis in patients with microsatellite stable (MSS) stage III colorectal cancer treated with 5-fluorouracil (5FU)-based adjuvant chemotherapy, suggesting their involvement in promoting chemoresistance. A novel IAP antagonist tolinapant (ASTX660) potently and rapidly downregulated cIAP1 in colorectal cancer models, demonstrating its robust on-target efficacy. In cells co-cultured with TNFα to mimic an inflammatory tumor microenvironment, tolinapant induced caspase-8-dependent apoptosis in colorectal cancer cell line models; however, the extent of apoptosis was limited because of inhibition by the caspase-8 paralogs FLIP and, unexpectedly, caspase-10. Importantly, tolinapant-induced apoptosis was augmented by FOLFOX in human colorectal cancer and murine organoid models in vitro and in vivo, due (at least in part) to FOLFOX-induced downregulation of class I histone deacetylases (HDAC), leading to acetylation of the FLIP-binding partner Ku70 and downregulation of FLIP. Moreover, the effects of FOLFOX could be phenocopied using the clinically relevant class I HDAC inhibitor, entinostat, which also induced acetylation of Ku70 and FLIP downregulation. Further analyses revealed that caspase-8 knockout RIPK3-positive colorectal cancer models were sensitive to tolinapant-induced necroptosis, an effect that could be exploited in caspase-8-proficient models using the clinically relevant caspase inhibitor emricasan. Our study provides evidence for immediate clinical exploration of tolinapant in combination with FOLFOX in poor prognosis MSS colorectal cancer with elevated cIAP1/2 expression.

Dying to Survive-The p53 Paradox

The p53 tumour suppressor is best known for its canonical role as "guardian of the genome", activating cell cycle arrest and DNA repair in response to DNA damage which, if irreparable or sustained, triggers activation of cell death. However, despite an enormous amount of work identifying the breadth of the gene regulatory networks activated directly and indirectly in response to p53 activation, how p53 activation results in different cell fates in response to different stress signals in homeostasis and in response to p53 activating anti-cancer treatments remains relatively poorly understood. This is likely due to the complex interaction between cell death mechanisms in which p53 has been activated, their neighbouring stressed or unstressed cells and the local stromal and immune microenvironment in which they reside. In this review, we evaluate our understanding of the burgeoning number of cell death pathways affected by p53 activation and how these may paradoxically suppress cell death to ensure tissue integrity and organismal survival. We also discuss how these functions may be advantageous to tumours that maintain wild-type p53, the understanding of which may provide novel opportunity to enhance treatment efficacy.

Development of a protein signature to enable clinical positioning of IAP inhibitors in colorectal cancer

Resistance to chemotherapy-induced cell death is a major barrier to effective treatment of solid tumours such as colorectal cancer, CRC. Herein, we present a study aimed at developing a proteomics-based predictor of response to standard-of-care (SoC) chemotherapy in combination with antagonists of IAPs (inhibitors of apoptosis proteins), which have been implicated as mediators of drug resistance in CRC. We quantified the absolute expression of 19 key apoptotic proteins at baseline in a panel of 12 CRC cell lines representative of the genetic diversity seen in this disease to identify which proteins promote resistance or sensitivity to a model IAP antagonist [birinapant (Bir)] alone and in combination with SoC chemotherapy (5FU plus oxaliplatin). Quantitative western blotting demonstrated heterogeneous expression of IAP interactome proteins across the CRC cell line panel, and cell death analyses revealed a widely varied response to Bir/chemotherapy combinations. Baseline protein expression of cIAP1, caspase-8 and RIPK1 expression robustly correlated with response to Bir/chemotherapy combinations. Classifying cell lines into 'responsive', 'intermediate' and 'resistant' groups and using linear discriminant analysis (LDA) enabled the identification of a 12-protein signature that separated responders to Bir/chemotherapy combinations in the CRC cell line panel with 100% accuracy. Moreover, the LDA model was able to predict response accurately when cells were cocultured with Tumour necrosis factor-alpha to mimic a pro-inflammatory tumour microenvironment. Thus, our study provides the starting point for a proteomics-based companion diagnostic that predicts response to IAP antagonist/SoC chemotherapy combinations in CRC.

Molecular subtype-specific responses of colon cancer cells to the SMAC mimetic Birinapant

Colorectal cancer is a molecularly heterogeneous disease. Responses to genotoxic chemotherapy in the adjuvant or palliative setting vary greatly between patients, and colorectal cancer cells often resist chemotherapy by evading apoptosis. Antagonists of an inhibitor of apoptosis proteins (IAPs) can restore defective apoptosis signaling by degrading cIAP1 and cIAP2 proteins and by inhibition of XIAP. Due to the multiple molecular mechanisms-of-action of these targets, responses to IAP antagonist may differ between molecularly distinct colon cancer cells. In this study, responses to the IAP antagonist Birinapant and oxaliplatin/5-fluorouracil (5-FU) were investigated in 14 colon cancer cell lines, representing the consensus molecular subtypes (CMS). Treatment with Birinapant alone did not result in a substantial increase in apoptotic cells in this cell line panel. Annexin-V/PI assays quantified by flow cytometry and high-content screening showed that Birinapant increased responses of CMS1 and partially CMS3 cell lines to oxaliplatin/5-FU, whereas CMS2 cells were not effectively sensitized. FRET-based imaging of caspase-8 and -3 activation validated these differences at the single-cell level, with CMS1 cells displaying sustained activation of caspase-8-like activity during Birinapant and oxaliplatin/5-FU co-treatment, ultimately activating the intrinsic mitochondrial apoptosis pathway. In CMS2 cell lines, Birinapant exhibited synergistic effects in combination with TNFα, suggesting that Birinapant can restore extrinsic apoptosis signaling in the context of inflammatory signals in this subtype. To explore this further, we co-cultured CMS2 and CMS1 colon cancer cells with peripheral blood mononuclear cells. We observed increased cell death during Birinapant single treatment in these co-cultures, which was abrogated by anti-TNFα-neutralizing antibodies. Collectively, our study demonstrates that IAP inhibition is a promising modulator of response to oxaliplatin/5-FU in colorectal cancers of the CMS1 subtype, and may show promise as in the CMS2 subtype, suggesting that molecular subtyping may aid as a patient stratification tool for IAP antagonists in this disease.

Abstract 1386: 5-FU-induces PD-L1 in the absence of p53 via ATM dependent activation of STAT3

Background: 5-Fluorouracil (5-FU) is an antimetabolite DNA damaging chemotherapeutics that forms the backbone for treatment of colorectal cancer (CRC) and currently it is the most commonly used anti-cancer drug. Despite the central role played by p53 in regulating DNA-damage response, its role in regulating response to 5-FU remains unclear.

Methods: To better understand how p53 status differentially affects response to 5-FU, we conducted a detailed phenotypic analysis of the effects of 5FU alone or in combination with Oxaliplatin on cell cycle, DNA damage, cell death, DNA repair and consequent signaling across a panel of p53 isogenic, null and mutant CRC cell lines.

Results: This revealed that while 5-FU induce both p53 dependent and independent cell death, the effects on cell cycle are profoundly different, with p53 deficient cells accumulating in S-phase with prolongated double strand breaks. This is not observed in p53 proficient cells, which are protected from these effects by induction of the p53 target gene p21/CDKN1A limiting cell cycle and enforcing suppression of cell cycle genes. Notably in both p53 or p21 deficient cells 5-FU results in sustained DNA damage-signaling through ATR and ATM in p53 deficient cells which preferentially arrest in S-phase. Importantly, these cells exhibit significant induction of cell surface expression of programmed death ligand-1 (PD-L1), which is enhanced by ATR inhibition concomitant with increased ATM activation. Further analysis indicates that transcriptional induction of PD-L1 is mediated through STAT3 in an interferon independent manner.

Conclusion: This work adds significantly to our understanding of how p53 status impacts 5FU induced cell cycle arrest and response in CRC patients and have important implications for understanding 5-FU response, particularly in future combinations with immune checkpoint inhibitors.

Citation Format: Tamas Sessler, Fiammetta Falcone, Peter Gallagher, Timothy Wright, Kienan Savage, Daniel B. Longley, Simon S. McDade. 5-FU-induces PD-L1 in the absence of p53 via ATM dependent activation of STAT3 [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1386.

The pseudo-caspase FLIP(L) regulates cell fate following p53 activation

p53 is the most frequently mutated, well-studied tumor-suppressor gene, yet the molecular basis of the switch from p53-induced cell-cycle arrest to apoptosis remains poorly understood. Using a combination of transcriptomics and functional genomics, we unexpectedly identified a nodal role for the caspase-8 paralog and only human pseudo-caspase, FLIP(L), in regulating this switch. Moreover, we identify FLIP(L) as a direct p53 transcriptional target gene that is rapidly up-regulated in response to Nutlin-3A, an MDM2 inhibitor that potently activates p53. Genetically or pharmacologically inhibiting expression of FLIP(L) using siRNA or entinostat (a clinically relevant class-I HDAC inhibitor) efficiently promoted apoptosis in colorectal cancer cells in response to Nutlin-3A, which otherwise predominantly induced cell-cycle arrest. Enhanced apoptosis was also observed when entinostat was combined with clinically relevant, p53-activating chemotherapy in vitro, and this translated into enhanced in vivo efficacy. Mechanistically, FLIP(L) inhibited p53-induced apoptosis by blocking activation of caspase-8 by the TRAIL-R2/DR5 death receptor; notably, this activation was not dependent on receptor engagement by its ligand, TRAIL. In the absence of caspase-8, another of its paralogs, caspase-10 (also transcriptionally up-regulated by p53), induced apoptosis in Nutlin-3A-treated, FLIP(L)-depleted cells, albeit to a lesser extent than in caspase-8-proficient cells. FLIP(L) depletion also modulated transcription of canonical p53 target genes, suppressing p53-induced expression of the cell-cycle regulator p21 and enhancing p53-induced up-regulation of proapoptotic PUMA. Thus, even in the absence of caspase-8/10, FLIP(L) silencing promoted p53-induced apoptosis by enhancing PUMA expression. Thus, we report unexpected, therapeutically relevant roles for FLIP(L) in determining cell fate following p53 activation.

The SCFSkp2 ubiquitin ligase complex modulates TRAIL-R2-induced apoptosis by regulating FLIP(L)

TRAIL-R2 (DR5) is a clinically-relevant therapeutic target and a key target for immune effector cells. Herein, we identify a novel interaction between TRAIL-R2 and the Skp1-Cullin-1-F-box (SCF) Cullin-Ring E3 Ubiquitin Ligase complex containing Skp2 (SCF^Skp2). We find that SCF^Skp2 can interact with both TRAIL-R2’s pre-ligand association complex (PLAC) and ligand-activated death-inducing signalling complex (DISC). Moreover, Cullin-1 interacts with TRAIL-R2 in its active NEDDylated form. Inhibiting Cullin-1’s DISC recruitment using the NEDDylation inhibitor MLN4924 (Pevonedistat) or siRNA increased apoptosis induction in response to TRAIL. This correlated with enhanced levels of the caspase-8 regulator FLIP at the TRAIL-R2 DISC, particularly the long splice form, FLIP(L). We subsequently found that FLIP(L) (but not FLIP(S), caspase-8, nor the other core DISC component FADD) interacts with Cullin-1 and Skp2. Importantly, this interaction is enhanced when FLIP(L) is in its DISC-associated, C-terminally truncated p43-form. Prevention of FLIP(L) processing to its p43-form stabilises the protein, suggesting that by enhancing its interaction with SCF^Skp2, cleavage to the p43-form is a critical step in FLIP(L) turnover. In support of this, we found that silencing any of the components of the SCF^Skp2 complex inhibits FLIP ubiquitination, while overexpressing Cullin-1/Skp2 enhances its ubiquitination in a NEDDylation-dependent manner. DISC recruitment of TRAF2, previously identified as an E3 ligase for caspase-8 at the DISC, was also enhanced when Cullin-1’s recruitment was inhibited, although its interaction with Cullin-1 was found to be mediated indirectly via FLIP(L). Notably, the interaction of p43-FLIP(L) with Cullin-1 disrupts its ability to interact with FADD, caspase-8 and TRAF2. Collectively, our results suggest that processing of FLIP(L) to p43-FLIP(L) at the TRAIL-R2 DISC enhances its interaction with co-localised SCF^Skp2, leading to disruption of p43-FLIP(L)’s interactions with other DISC components and promoting its ubiquitination and degradation, thereby modulating TRAIL-R2-mediated apoptosis.

FLIP(L): the pseudo-caspase

Possessing structural homology with their active enzyme counterparts but lacking catalytic activity, pseudoenzymes have been identified for all major enzyme groups. Caspases are a family of cysteine‐dependent aspartate‐directed proteases that play essential roles in regulating cell death and inflammation. Here, we discuss the only human pseudo‐caspase, FLIP(L), a paralog of the apoptosis‐initiating caspases, caspase‐8 and caspase‐10. FLIP(L) has been shown to play a key role in regulating the processing and activity of caspase‐8, thereby modulating apoptotic signaling mediated by death receptors (such as TRAIL‐R1/R2), TNF receptor‐1 (TNFR1), and Toll‐like receptors. In this review, these canonical roles of FLIP(L) are discussed. Additionally, a range of nonclassical pseudoenzyme roles are described, in which FLIP(L) functions independently of caspase‐8. These nonclassical pseudoenzyme functions enable FLIP(L) to play key roles in the regulation of a wide range of biological processes beyond its canonical roles as a modulator of cell death.

A revised model of TRAIL-R2 DISC assembly explains how FLIP(L) can inhibit or promote apoptosis

The long FLIP splice form FLIP(L) can act as both an inhibitor and promoter of caspase‐8 at death‐inducing signalling complexes (DISCs) formed by death receptors such as TRAIL‐R2 and related intracellular complexes such as the ripoptosome. Herein, we describe a revised DISC assembly model that explains how FLIP(L) can have these opposite effects by defining the stoichiometry (with respect to caspase‐8) at which it converts from being anti‐ to pro‐apoptotic at the DISC. We also show that in the complete absence of FLIP(L), procaspase‐8 activation at the TRAIL‐R2 DISC has significantly slower kinetics, although ultimately the extent of apoptosis is significantly greater. This revised model of DISC assembly also explains why FLIP's recruitment to the TRAIL‐R2 DISC is impaired in the absence of caspase‐8 despite showing that it can interact with the DISC adaptor protein FADD and why the short FLIP splice form FLIP(S) is the more potent inhibitor of DISC‐mediated apoptosis.

Structural determinants of DISC function: new insights into death receptor-mediated apoptosis signalling

Death receptors are members of the tumour necrosis factor (TNF) receptor superfamily characterised by an ~80 amino acid long alpha-helical fold, termed the death domain (DD). Death receptors diversified during early vertebrate evolution indicating that the DD fold has plasticity and specificity that can be easily adjusted to attain additional functions. Eight members of the death receptor family have been identified in humans, which can be divided into four structurally homologous groups or clades, namely: the p75(NTR) clade (consisting of ectodysplasin A receptor, death receptor 6 (DR6) and p75 neurotrophin (NTR) receptor); the tumour necrosis factor receptor 1 clade (TNFR1 and DR3), the CD95 clade (CD95/FAS) and the TNF-related apoptosis-inducing ligand receptor (TRAILR) clade (TRAILR1 and TRAILR2). Receptors in the same clade participate in similar processes indicating that structural diversification enabled functional specialisation. On the surface of nearly all human cells multiple death receptors are expressed, enabling the cell to respond to a plethora of external signals. Activation of different death receptors converges on the activation of three main signal transduction pathways: nuclear factor-κB-mediated differentiation or inflammation, mitogen-associated protein kinase-mediated stress response and caspase-mediated apoptosis. While the ability to induce cell death is true for nearly all DRs, the FAS and TRAILR clades have specialised in inducing cell death. Here we summarise recent discoveries about the molecular regulation and structural requirements of apoptosis induction by death receptors and discuss how this information can be used to better explain the biological functions, similarities and distinguishing features of death receptors.

Resistance to TRAIL in non-transformed cells is due to multiple redundant pathways

Tumour necrosis factor-related apoptosis-inducing ligand (TRAIL) is a cytokine and a selective inducer of apoptosis in a range of tumour cells, but not in normal, untransformed cells. A large number of chemotherapeutics as well as biological agents are being tested for their potential to sensitise resistant tumour cells to TRAIL as a means to broaden the range of tumours treatable with TRAIL. However, because of the incomplete understanding of the mechanism(s) underlying TRAIL resistance in non-malignant cells, it is unpredictable whether the effect of these sensitisers will be restricted to tumour cells or they would also sensitise non-transformed cells causing unwanted toxicity. In this study, we carried out a systematic analysis of the mechanisms driving TRAIL resistance in non-transformed cells. We found that cellular FLICE-like inhibitory protein, anti-apoptotic B-cell lymphoma 2 proteins, and X-linked inhibitor of apoptosis protein were independently able to provide resistance to TRAIL. Deficiency of only one of these proteins was not sufficient to elicit TRAIL sensitivity, demonstrating that in non-transformed cells multiple pathways control TRAIL resistance and they act in a redundant manner. This is contrary to the resistance mechanisms found in tumour cell types, many of them tend to rely on a single mechanism of resistance. Supporting this notion we found that 76% of TRAIL-resistant cell lines (13 out of 17) expressed only one of the above-identified anti-apoptotic proteins at a high level (≥1.2-fold higher than the mean expression across all cell lines). Furthermore, inhibition or knockdown of the single overexpressed protein in these tumour cells was sufficient to trigger TRAIL sensitivity. Therefore, the redundancy in resistance pathways in non-transformed cells may offer a safe therapeutic window for TRAIL-based combination therapies where selective sensitisation of the tumour to TRAIL can be achieved by targeting the single non-redundant resistance pathway.