Next Article in Journal
Prognostic Impact of let-7e MicroRNA and Its Target Genes in Localized High-Risk Intestinal GIST: A Spanish Group for Research on Sarcoma (GEIS) Study
Next Article in Special Issue
Patient-Derived Xenograft and Organoid Models for Precision Medicine Targeting of the Tumour Microenvironment in Head and Neck Cancer
Previous Article in Journal
FGFR4 Inhibitor BLU9931 Attenuates Pancreatic Cancer Cell Proliferation and Invasion While Inducing Senescence: Evidence for Senolytic Therapy Potential in Pancreatic Cancer
Previous Article in Special Issue
Fast, In Vivo Model for Drug-Response Prediction in Patients with B-Cell Precursor Acute Lymphoblastic Leukemia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Implementing Systems Modelling and Molecular Imaging to Predict the Efficacy of BCL-2 Inhibition in Colorectal Cancer Patient-Derived Xenograft Models

by
Alice C. O’Farrell
1,†,
Monika A. Jarzabek
1,†,
Andreas U. Lindner
2,
Steven Carberry
2,
Emer Conroy
3,
Ian S. Miller
1,
Kate Connor
1,
Liam Shiels
1,
Eugenia R. Zanella
4,5,
Federico Lucantoni
2,
Adam Lafferty
1,
Kieron White
1,
Mariangela Meyer Villamandos
2,
Patrick Dicker
6,
William M. Gallagher
3,
Simon A. Keek
7,
Sebastian Sanduleanu
7,
Philippe Lambin
7,8,
Henry C. Woodruff
7,8,
Andrea Bertotti
4,5,
Livio Trusolino
4,5,
Annette T. Byrne
1,3,‡ and
Jochen H. M. Prehn
2,*,‡
add Show full author list remove Hide full author list
1
Precision Cancer Medicine Group, Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, D02 YN77 Dublin, Ireland
2
Department of Physiology and Medical Physics, Centre for Systems Medicine, Royal College of Surgeons in Ireland, D02 YN77 Dublin, Ireland
3
UCD School of Biomolecular and Biomedical Science, UCD Conway Institute, University College Dublin, Belfield, D04 V1W8 Dublin, Ireland
4
Candiolo Cancer Institute—FPO IRCCS, Candiolo, 10060 Torino, Italy
5
Department of Oncology, University of Torino, Candiolo, 10060 Torino, Italy
6
Department of Epidemiology and Public Health Medicine, Royal College of Surgeons in Ireland, D02 YN77 Dublin, Ireland
7
The D-Lab, Department of Precision Medicine, GROW—School for Oncology, Maastricht University, 6229 ER Maastricht, The Netherlands
8
Department of Radiology and Nuclear Imaging, Maastricht University Medical Center, 6229 HX Maastricht, The Netherlands
*
Author to whom correspondence should be addressed.
These authors contributed equally to the work.
Joint senior authors.
Cancers 2020, 12(10), 2978; https://doi.org/10.3390/cancers12102978
Submission received: 16 September 2020 / Accepted: 9 October 2020 / Published: 14 October 2020
(This article belongs to the Special Issue Patient-Derived Xenograft-Models in Cancer Research)

Abstract

:

Simple Summary

Drugs that sensitise tumours to chemotherapy by enhancing cell death signalling are of significant clinical interest. However, it is challenging to determine which colorectal cancer patients may benefit from such sensitisers. The ability to predict this would be advantageous. Here we show that protein profiling combined with mathematical modelling identifies responsive tumours. Using our modelling method, we predicted the effect of adding a sensitizer drug to chemotherapy in two patient-derived colorectal tumours. We grew the tumours in mice, treated animals with these drugs and performed PET/CT imaging. The predicted “sensitive” tumours were smaller when the sensitising drug was added to chemotherapy whilst it did not further reduce tumour size in “non-sensitive” tumours, thus validating our prediction. PET imaging also supported our predictions. CT analysis (radiomics) revealed features that distinguished the two tumours. This was the first application of radiomic analyses to PDX derived CT data.

Abstract

Resistance to chemotherapy often results from dysfunctional apoptosis, however multiple proteins with overlapping functions regulate this pathway. We sought to determine whether an extensively validated, deterministic apoptosis systems model, ‘DR_MOMP’, could be used as a stratification tool for the apoptosis sensitiser and BCL-2 antagonist, ABT-199 in patient-derived xenograft (PDX) models of colorectal cancer (CRC). Through quantitative profiling of BCL-2 family proteins, we identified two PDX models which were predicted by DR_MOMP to be sufficiently sensitive to 5-fluorouracil (5-FU)-based chemotherapy (CRC0344), or less responsive to chemotherapy but sensitised by ABT-199 (CRC0076). Treatment with ABT-199 significantly improved responses of CRC0076 PDXs to 5-FU-based chemotherapy, but showed no sensitisation in CRC0344 PDXs, as predicted from systems modelling. 18F-Fluorodeoxyglucose positron emission tomography/computed tomography (18F-FDG-PET/CT) scans were performed to investigate possible early biomarkers of response. In CRC0076, a significant post-treatment decrease in mean standard uptake value was indeed evident only in the combination treatment group. Radiomic CT feature analysis of pre-treatment images in CRC0076 and CRC0344 PDXs identified features which could phenotypically discriminate between models, but were not predictive of treatment responses. Collectively our data indicate that systems modelling may identify metastatic (m)CRC patients benefitting from ABT-199, and that 18F-FDG-PET could independently support such predictions.

1. Introduction

Several chemotherapeutic agents kill cancer cells by initiating apoptosis. Dysfunctional apoptosis is recognised as a key factor in malignant progression and development of chemotherapy resistance in cancer [1,2]. The mitochondrial apoptosis pathway plays an important role in chemotherapy-induced cell death [1]. BCL-2 family proteins are master regulators of mitochondrial apoptosis, and act on mitochondria to induce, or inhibit, mitochondrial outer membrane permeabilization (MOMP). MOMP is achieved through the activation and homo-oligomerisation of two channel-forming, pro-apoptotic multi-domain proteins within the B-cell lymphoma 2 (BCL-2) family, bcl-2-like protein 4 (BAX) and Bcl-2 homologous antagonist/killer (BAK). Monomeric activated BAX and BAK are sequestered and neutralised by a number of anti-apoptotic multi-domain proteins within the BCL-2 family. The anti-apoptotic function of these proteins are antagonised by a third set of proteins within the family, the Bcl-2 homology-3 (BH3)-only proteins. These are activated by transcription or post-translational modification in response to cellular stress [3,4]. Thus, whether cells undergo mitochondrial apoptosis or not is regulated by a complex network of BCL-2 protein interactions.
Small molecule apoptosis sensitisers, such as BCL-2 antagonists restore or induce cell death signalling in cancer cells. There are now several selective BCL-2 antagonists targeting anti-apoptotic BCL-2 family proteins, with the selective BCL-2 antagonist ABT-199 (venetoclax) now approved for the treatment of acute myeloid leukaemia, chronic and small lymphocytic lymphoma patients who received at least one prior therapy [5]. Unlike many haematological cancers, which often overexpress a single anti-apoptotic BCL-2 protein, solid tumours routinely express several pro- and anti-apoptotic proteins [6]. In the setting of solid tumours, it is therefore more challenging to determine which cancer type and which individual patients may best benefit from a given BCL-2 antagonist, highlighting the need for novel systems-based patient stratification tools.
We previously developed and extensively validated a computational tool, DR_MOMP, which determines the sensitivity of cancer cells to undergo apoptosis [7,8,9]. DR_MOMP is a systems model based on ordinary differential equations (ODEs). It incorporates BCL-2 protein levels of five BCL-2 family members BAK, BAX, BCL-2, BCL-2-like 1-long (BCLxL), and myeloid leukaemia cell differentiation protein (MCL-1), measured by quantitative Western blotting or reverse phase protein array analysis, their interaction kinetics, as well as translocation and degradation kinetics of these proteins to calculate the amount of BH3-only proteins (‘stress dose’) required to induce the process of MOMP [7,9]. DR_MOMP hence delivers for each tumour/cell type a numerical value (‘stress dose’) that is a prognostic marker of clinical outcome in colorectal cancer (CRC) [7,8,9], and capable of predicting the response of CRC [7] and breast cancer cells [10] to genotoxic therapy. Here, we extended the application of DR_MOMP to include the effect of ABT-199 on the apoptosis signalling network, and to determine whether the deterministic model could also be used as a stratification tool for the BCL-2 antagonist ABT-199 in in vivo CRC patient derived xenograft (PDX) models (Figure 1). A secondary, hypothesis generating-aspect of the research question was to investigate whether tumour glucose metabolism and/or radiomics features derived from preclinical 18F-Flourodeoxyglucose positron emission tomography/computed tomography (18F-FDG-PET/CT) supported model predictions.

2. Results

2.1. In Vivo PDX Response to ABT-199 Combined with FOLFOX Validates DR_MOMP Predictions

As per our previous studies that demonstrated the ability of DR_MOMP to accurately predict responses of breast cancer cells to BCL-2 antagonists when combined with chemotherapy [10], we started our investigations by confirming that DR_MOMP can predict the sensitisation of CRC cell lines to genotoxic therapeutics by the BCL-2 antagonist ABT-199. We determined BCL-2 protein profiles in DLD-1 and HCT 116 colon cancer cells by quantitative Western Blotting (Figures S1 and S2A, Table S1). We next calculated the ‘stress dose’ required to induce mitochondrial apoptosis in response to genotoxic chemotherapy by employing DR_MOMP systems modelling. DR_MOMP calculated an equivalent apoptosis sensitivity in HCT-116 and DLD-1 cells (Figure S2B, Table S2). Upon treatment with 5-FU/OX, as predicted by the model, DLD-1 and HCT-116 cells did not show a significant difference in mean cell survival (t-tests, adjusted Holms p > 0.49; statistic performed on data presented in Figure S2C,D). However, DR_MOMP analysis indicated that HCT-116 were predicted to be more efficiently sensitised to chemotherapy by ABT-199 than DLD-1 cells. ABT-199 concentration ranges were chosen to reflect therapeutic plasma concentrations previously been shown to reach approximately 5 µM [11]. To validate these model predictions, we performed a 6 × 6 dose matrix assay using a wide range of 5-fluorouracil (5-FU)/oxaliplatin (OX) and ABT-199 concentrations (Figure S2C,D). Data were also processed using isobologram analysis to evaluate drug interactions (Figure S2E). Cell survival analysis confirmed the greater synergistic effect in HCT-116 cells at high ABT-199 concentrations, in line with DR_MOMP predictions.
To explore the above confirmed ability of DR_MOMP to predict responses to BCL-2 antagonists in a pre-clinical in vivo setting, we first profiled BCL-2 proteins in several metastatic CRC PDX tissue samples that were all originally derived from metastatic (m)CRC patient liver metastases, established by Bertotti et al. [12,13]. We identified two PDXs, CRC0344 and CRC0076, which differed in their predicted responses to genotoxic chemotherapy, and their predicted sensitisation by ABT-199. BCL-2 family protein levels in PDX samples were quantified using cell lines with known absolute BCL-2 protein levels as standards, and were in line with previous studies in CRC tissue samples [7,8]. Data are presented in Table S1 and Figure 2A. Western blots are provided in Figure S3. Using these BCL-2 protein profiles, we calculated the intrinsic sensitivity of the two models to genotoxic stress caused by 5-FU either with or without ABT-199, varying the ABT-199 doses between 0 to 10 μM (Figure 2B, Table S2). In the absence of ABT-199, the calculated DR_MOMP stress dose of CRC0076 was 4271 nM, indicating resistance to 5-FU/OX treatment in this model, which is in line with previous findings [7]. CRC0344 showed a lower DR_MOMP stress doses of 2601 nM (Figure 2B) indicating a higher sensitivity to 5-FU/OX treatment. In the presence of a 10.0 µM ABT-199 dose, the model predicted that the DR_MOMP stress dose was reduced in both models. Of note, the model predicted that the stress dose in CRC0076 was reduced to 1506 nM, thus rendering it potentially more vulnerable to 5-FU/OX treatment upon the introduction of ABT-199. BCL-2 concentrations were 5772 nM and 1980 nM for PDX models CRC0076 and CRC0344 respectively (Figure 2A). Overall, PDX CRC0076 tumours were predicted to be more resistant to genotoxic stress than CRC0344, with the addition of 10 µM ABT-199 sensitising this model. Conversely, CRC0344 was sufficiently sensitive to 5-FU/OX alone.
In order to interrogate the effect of FOLFOX, ABT-199, and combination treatment, we performed a blinded in vivo efficacy study using the two mCRC PDX models CRC0076 and CRC0344. These models were chosen based on previously confirmed engraftment rates (>90%) and their differing DR_MOMP-predicted values. Models (n = 13–14 mice/group) were treated with ABT-199 either with or without FOLFOX (5-FU + folinic acid + OX) [replicating the clinical regimen]. Tumour growth was measured using calipers over 4 weeks (Figure 2C,D). In CRC0076, following 28 days of treatment, combination therapy significantly decreased average tumour size (p = 0.02) compared to vehicle. However, there was no significant difference between FOLFOX alone and vehicle-treated tumours (p = 0.08) (Figure 2C). Conversely, in CRC0344 following 28 days of treatment, both FOLFOX alone and combination therapy decreased average tumour size compared to vehicle (p = 0.0013 and p = 0.0001 respectively). Thus, the addition of ABT-199 conferred no additional benefit over FOLFOX alone in this model (Figure 2D). These data are in line with DR_MOMP predictions where CRC0344 was predicted to require a substantially lower genotoxic stress dose than CRC0076, thus making it more sensitive to treatment with FOLFOX alone (Figure 2B). Conversely, CRC0076 was sensitised by the addition of ABT-199.

2.2. ABT-199 Treatment In Vivo Leads to Decrease of p53 Upregulated Modulator of Apoptosis (PUMA) Levels in the CRC0076 PDX Model

On treatment day 28, c.4 h after the final dose of ABT-199, tumours taken from three representative mice from each treatment group were analysed using immunohistochemistry (IHC) and Western blot (This sample size is deemed appropriate for a pharmacodynamic study as previously outlined by Workman et al. [14]). Tissue was used to study the effects of ABT-199 in the two PDX models.
In CRC0076, we found lower levels of p53 upregulated modulator of apoptosis (PUMA) in the three samples treated with ABT-199 alone when compared to vehicle (Figure 3A, p = 0.049). Moreover, the combination treated CRC0076 tumours also showed reduced PUMA levels (Figure 3A, p = 0.02). BCL2-Interacting Mediator of cell death (BIM) levels also tended to be lowered by ABT-199 treatment, but this effect was more variable than the effect observed with PUMA. In contrast, we did not observe an effect of ABT-199 on PUMA or BIM levels in the FOLFOX sensitive CRC0344 model (Figure 3B). Full Western blots are presented in Figure S4. Histological examination of CRC0076 (ANOVA p = 0.53) and CRC0344 PDX tumours (ANOVA p = 0.63) revealed no statistically significant change between the treatment groups in cleaved caspase-3 levels after 28 days of treatment (Figure 3C,D). Our data suggest a numerical trend towards decreased cell proliferation (Ki-67 staining) in CRC0076 (which expresses high levels of BCL-2) following ABT-199 treatment (Figure S5).

2.3. Decreased Tumour Glucose Uptake as an Early Biomarker of Response to FOLFOX and ABT-199 Combination Therapy

We next sought to determine whether variability in tumour glucose uptake and metabolism was an early response biomarker in the predicted combination responder, CRC0076, and the predicted FOLFOX alone responder, CRC0344. For each model, six mice from each treatment group were randomly selected to undergo 18F-FDG-PET/CT.
To determine the effect of treatment on glucose uptake and metabolism (indicated by tumour standard uptake value (SUV) calculations) in the PDX tumours, pre- and post-treatment scans (taken 2 weeks after treatment initiation) were performed in the same animals (representative images shown in Figure 4A). When comparing pre- and post-treatment SUVs, a significant post-treatment decrease in glucose uptake and metabolism was only evident (decreased SUVmean) in the combination group of the CRC0076 model (p = 0.01). In the CRC0344 model no significant change in post-treatment glucose uptake and metabolism was observed in any treatment group (Figure 4B).

2.4. Radiomic CT Feature Analysis in DR_MOMP Predicted Combination-Only Responder (CRC0076) and FOLFOX Alone Responder (CRC0344) PDX Models

To further mine the available PET/CT imaging data (Section 2.3, Figure 4), we sought to interrogate available CT data using a radiomic pipeline. This hypothesis-generating study aimed to assess whether radiomic features extracted from pre-treatment images could predict end-study (4- weeks) treatment outcome in the combination-only responder (CRC0076) and/or the FOLFOX alone responder (CRC0344).
CT data from pre-treatment and 2 weeks post-treatment mouse tumour scans underwent additional radiomic analysis (workflow and image segmentation described in Figure 5A–C); Firstly, the relationship between CT-derived radiomic features and the underlying phenotype of DR_MOMP predicted combination therapy responder (CRC0076) and FOLFOX alone responder (CRC0344) models was studied. Subsequently, the baseline CT features that were able to discriminate between models in univariate analysis were interrogated in scans taken at 2-weeks post treatment to assess their ability to predict end study outcome. For each PDX (CRC0076 and CRC0344) n = 36 pre-treatment images were used as a training cohort, whilst n = 9 pre-treatment images were retained as a validation cohort. The training cohort was used both to establish a radiomics pipeline (not previously established for PDX-derived CT data) and to identify potential radiomic biomarker candidates for further analysis.
Ultimately, seven features best discriminated between the DR_MOMP predicted combination therapy response (CRC0076) and the FOLFOX alone responder (CRC0344) (Figure 5D). Using the n = 36 training cohort, a radiomic random forest classifier was trained to separate the two PDXs using these seven identified features. The random forest classifier was then tested to determine if the two PDXs in the validation cohort of pre-treatment scans (n = 9) could be successfully separated. Performance of the radiomics classifier in the validation dataset was shown to have an area under the curve (AUC) of 0.80 (95% CI 0.408–1, Figure 5E).
We next assessed whether any of the seven CT-derived radiomic features could predict end-study (4-weeks) treatment outcome in the combination-only responder (CRC0076) and/or the FOLFOX alone responder (CRC0344). For all PDX treatment groups (vehicle, ABT-199, FOLFOX, combination in CRC0076 and CRC0344), the value of each feature at 2 weeks was compared to the pre-treatment values to determine treatment effect on the radiomics feature after 14 days. As before, the post-treatment CRC PDX CT images (n = 45) were segmented and radiomics features extracted. Data is presented in Figure S6.
In CRC0076 tumours treated with ABT-199, the Wavelet.HHL_ngtdm_Busyness feature intensity was significantly higher (p = 0.04526) after 2 weeks of treatment compared to its pre-treatment intensity (Figure S6A). Nevertheless, this did not predict significant treatment response in the ABT-199 cohort over 4 weeks (Figure 2C). All other significant changes in radiomic feature intensity observed after 2 weeks of treatment in either PDXs were seen across all cohorts including vehicle (Figure S6C,D,F,H,N). Ultimately, after two weeks of treatment, none of the 7 identified radiomics features predicted treatment efficacy after 4 weeks in either the combination sensitive CRC0076 tumour or the FOLFOX sensitive CRC0344.

3. Discussion

In this study we have investigated whether the BCL-2 systems model DR_MOMP could be exploited as a predictive tool for response to ABT-199 treatment in mCRC patients. In line with previous studies [9,16], the model successfully identified tumour cells resistant to 5-FU/OX treatment. To note, in a previous study using a panel of colon cancer cell lines (n = 6) we showed a correlation between DR_MOMP scores and sensitivity to genotoxic drugs 5-FU/OX ([7], Figure 3). Here, we demonstrate that DR_MOMP can identify mCRC tumours that benefit from the addition of ABT-199 to 5-FU/OX-based genotoxic therapy. ABT-199 was approved by the Food and Drug Administration (FDA) in 2016 for treatment of chronic lymphocytic leukaemia and is currently incorporated in more than 150 ongoing clinical trials [17]. However, there are limited trials in solid tumours and none specifically related to CRC, which is surprising considering the well-known relationship between BCL-2 family protein expression to both prognosis and resistance to treatment in CRC [7,17,18,19,20]. CRC tumours show high heterogeneity in their expression of multiple pro- and anti-apoptotic BCL-2 family proteins [7,9]. Tools such as DR_MOMP may, as demonstrated here, identify individual mCRC patients who might benefit from BCL-2 antagonist therapy in future clinical studies.
In vitro, DR_MOMP predicted for DLD-1 and HCT-116 CRC cell lines similar responses to ABT-199 and FOLFOX despite differences in their BCL-2 and BCL-xL levels. In both DLD-1 and HCT-116 cells, the combination of ABT-199 with low concentrations of genotoxic drugs (0.2 µM OX plus 1 µM 5-FU) did not result in more cell death when compared to cultures treated with ABT-199 only (Figure S2C,D), and combination of ABT-199 with higher concentrations of OX plus 5FU (5.0 µM and 25 µM) were required to induce higher levels of cell death. Both cell lines showed a numerically homologous response in their in vitro decrease in cell viability and in silico decrease in DR_MOMP score (Figure S2B–D), however the possibility also remains that low concentrations of OX plus 5-FU were protecting cells from ABT-199 induced death by activating protective stress responses induced by DNA damage.
In the PDX models, which also showed varying BCL-2 family protein levels, while DR_MOMP predicted a differential response to ABT-199 and FOLFOX for CRC0076 and CRC0344, an n-fold greater BCL-2 concentration in one tumour did not translate to an n-fold higher requirement for BCL-2 antagonists (Figure 2A,B). Previous observations have suggested that it may be sufficient to quantify BCL-2 protein levels as stratification tool for ABT-199 [21]. However, our findings in both cell and tumour models, reinforce the assertion that individual BCL-2 family protein levels alone may be insufficient as a sole predictive biomarkers for ABT-199 response in CRC patients. This highlights the power of systems approaches to estimate drug responsiveness and further supports the role of systems-informed biomarkers in this setting.
To validate DR_MOMP predictions in an in vivo, clinically relevant setting, two PDX models were studied. CRC0076 was identified by DR_MOMP as having higher resistance to 5-FU/OX (thus rendering it potentially more susceptible to sensitisation by ABT-199) compared to CRC0344 (Figure 2B). It is noteworthy that CRC0076 is faster growing than CRC0344, reflecting clinical differences in tumour growth rates as evidenced in patients. As predicted, a significant in vivo response to FOLFOX alone was elicited in PDX CRC0344, which was not significantly different to the response to ABT-199 + FOLFOX combination. In contrast, within PDX CRC0076 a significant response was only evident upon treatment with the ABT-199 + FOLFOX combination. The DR_MOMP prediction for which PDX would benefit from the addition of ABT-199 to FOLFOX was thus validated in our in vivo model. As ABT-199 is not an approved drug for the treatment of CRC (and therefore was not used to treat the two patients from which PDXs were derived), we were unable to establish if DR_MOMP was also predictive of patient outcome. Nevertheless, we have previously shown that DR_MOMP can predict response to adjuvant chemotherapy and neoadjuvant chemoradiation therapy in colon and rectal cancer patients. Moreover, DR_MOMP can also be used to calculate a predicted response of patients to BCL-2 antagonists [7,8,9]. Future studies to further determine the translational relevance of DR_MOMP are now warranted (e.g., large-scale PDX co-clinical trials). Ultimately, future implementation of this technology in the clinic will require a tailored clinical trial.
In growing tumours, the pro-apoptotic BH3-only protein PUMA is normally sequestered by BCL-2 and stabilised. In PDX CRC0076, which required ABT-199 for sensitisation, we observed significantly reduced levels of PUMA in ABT-199 treated animals, suggesting that it was released from BCL-2 (which is particularly high in this model) and degraded through the proteasome after ABT-199 treatment [22]. These experiments therefore suggest efficient target engagement in vivo. BIM protein levels are also subject to rapid proteasomal degradation [23], but were not significantly altered by ABT-199, although there was variability between the different groups. In PDX CRC0344, PUMA protein levels were not significantly altered by ABT-199 treatment, suggesting that any PUMA induced is functional (but rapidly degraded) and not sequestered by BCL-2. BIM levels were likewise not altered by the treatments.
Although histological examination of CRC0076 and CRC0344 PDX tumours revealed no significant change between the treatment groups in cleaved caspase-3 levels after 28 days of treatment, it is important to note that cleaved, active caspase-3 is likely only detectable in tissues for minutes to hours [24]. Moreover, BCL-2 family mediated mitochondrial membrane permeabilisation controls both caspase-dependent and caspase-independent cell death pathways such as necroptosis, autophagic cell death and necrosis (reviewed in [25]). Ki-67 staining showed a trend towards reduced cell proliferation in the high BCL-2 expressing CRC0076 PDXs upon ABT-199 treatment. These data suggest that BCL-2 might also control or correlate with cell proliferation, and that ABT-199 inhibits cell proliferation when administrated as a monotherapy (as also seen in ER-positive breast cancer [26]). This effect may be linked to the ability of BCL-2 to improve bioenergetics in cancer cells, which in turn would be inhibited by Bcl-2 and Bcl-xL antagonists [10]. However, we did not observe an inhibitory effect on cell proliferation in the combined treatment group. Here, the induction of apoptosis may have activated a compensatory effect on cell proliferation, a process that is known to occur both physiologically and in the cancer setting (reviewed in [27]).
We next implemented an 18F-FDG-PET imaging approach to determine if early changes in tumour glucose uptake and metabolism could predict treatment outcome. The rationale here was that in addition to controlling the process of MOMP during apoptosis, BCL-2 family proteins also regulate the bioenergetics status of cells and may act directly on the mitochondrial respiratory chain as discussed above [28,29,30,31,32]. Thus, we hypothesised that 18F-FDG-PET may predict response to therapy at an earlier time point than that seen using other clinical parameters (e.g., change in tumour size). Interestingly after two weeks of treatment, only CRC0076 tumours exposed to combination therapy showed a significant decrease in SUV mean (indicative of overall glucose metabolism). SUVmax (generated using the most intense 4 voxels per volume of interest (VOI) was also investigated but provided no additional insight). In this model, tumours in the combination group were significantly smaller than those treated with FOLFOX or ABT-199 alone. These data suggest that in tumours characterised by high BCL-2 levels, and predicted by DR_MOMP to respond favourably to the addition of ABT-199 to FOLFOX, 18F-FDG-PET may be an early response biomarker. Corroborating these results, a recent phase 1b study of ABT-199 combined with tamoxifen in ER and BCL-2–positive metastatic breast cancer treatment responses were pre-empted by metabolic responses using 18F-FDG-PET [26]. Importantly, this was the first clinical study to evaluate ABT-199 in a solid tumour. It remains unclear as to why glucose metabolism was not reduced in the CRC0344 model despite both FOLFOX alone and the combination tumours being significantly smaller than untreated controls. One possible explanation is that non-apoptotic functions of BCL-2 family proteins include effects on glucose metabolism [33,34]. In Figure 2A CRC0076 is seen to have elevated BCL-2 but lower BCL-xL levels compared to CRC0344. Moreover, the effect on PUMA and BIM following treatment was variable between the two models (Figure 3). It is possible that differential expression of other BCL-2 family members could underpin variance in glucose metabolism across the two models. Further mechanistic studies are warranted.
Next, CT data from pre-treatment and 2 weeks post-treatment PDX mouse tumour scans from DR_MOMP predicted combination-only responder (CRC0076) and FOLFOX alone responder (CRC0344) PDX models underwent additional radiomic analysis. Radiomics is a multistep process in which radiographic features based on shape, pixel intensities, and texture are extracted both from clinical and pre-clinical radiological images [35,36,37]. Radiomic analysis of CT data was undertaken due to the clinical relevance of CT in the CRC patient setting. In addition to providing novel methods of tumour classification and staging, radiomic analyses may also reveal critical information underpinning the biology of tumours [38]. Here, our CT radiomics pipeline identified 7 radiographic features, extracted from pre-treatment PDX CT images, that when combined into a random forest model were able to classify PDX CRC0076 from PDX CRC0344 (Figure 5D, AUC 0.8 95% CI 0.408–1). These hypothesis-generating data suggest that the model can discriminate the two PDXs based on radiomics features. Interestingly, 6 of 7 of these features were wavelet or multiscale textural features. Such features have previously been employed in patient tumour classification. Specifically, Chaddad et al. used wavelets to classify and discriminate whole slide colon cancer images into different grades of CRC with 93.33% accuracy [39]. This suggests that a CRC PDX CT radiomics pipeline is likely to have useful translational potential. Further validation to understand the predictive potential of wavelet or multiscale textural CT radiomic features in interventional CRC PDX population [40] and clinical trial datasets are warranted.
We subsequently investigated if any of the identified 7 features could individually predict end-study treatment outcome in the DR_MOMP predicted combination-only (CRC0076) responder or FOLFOX alone responder (CRC0344) models. Due to the limitation of treatment group size (n = 5–6), modelling was not undertaken on the post-treatment scan feature sets. Rather, we compared the values of each radiomics feature, extracted from 2 week post-treatment images, to pre-treatment values (Figure S6). In CRC0076 PDX tumours treated with ABT-199, after 2 weeks of treatment the intensity of feature Wavelet.HHL_ngtdm_Busyness was significantly higher (p = 0.04526) than at baseline (Figure S6A). This observation suggests that treatment with ABT-199 alone alters the radiographic phenotype of the tumour. However this feature change was not represented by a change in tumour volume at this time point. Moreover, this change did not precede a significant treatment response after 4 weeks (Figure 2C). All other significant changes in radiomic feature intensity observed after 2 weeks of treatment in either PDX were seen across all cohorts including vehicle (Figure S6C,D,F,H,N). As feature intensities changed equivalently across cohorts, these changes are likely associated with inherent phenotypic changes in PDX tumours over time, including tumour volume. Future studies should consider the implication of tumour size; Specifically, comparing tumours of a similar volume would decrease “noise” and improve comparisons.
A limitation of the current radiomics study is the relatively small number of available PET/CT images compared with previous reports [41]. Here, radiomics was performed using limited available data from a preclinical study powered in the first-instance for drug efficacy comparisons. Moreover, test-retest analyses were not performed. Nevertheless, to the best of our knowledge we have for the first time established a CT radiomics pipeline for application in CRC PDX studies. Such analyses will likely become increasingly important as the PDX field evolves, with an expanding focus on the use of orthotopic models requiring clinically relevant imaging [40,42,43]. Indeed, the translational relevance of pre-clinical radiomic ‘discovery’ studies in animal models is now becoming evident. For example, Zinn et al. recently employed radiomic analyses on glioblastoma magnetic resonance data sets from mice and humans, and showed that selected imaging features are conserved across species [44]. In the current context, future work employing CT scan data from large scale population-based CRC PDX trials is now warranted. Increasing the number of CT images available for analysis could identify more robust predictive radiomics signatures having utility as predictive biomarkers.

4. Materials and Methods

4.1. Quantitative Western Blotting

Western blot analysis was performed as previously described [8]. Briefly, tumour specimen pellets (tumours were halved following removal and all material was used for blotting) were directly homogenised and mixed with lysis buffer containing 0.5 mmol/L Tris(hydroxymethyl)aminomethane—hydrogen chloride (Tris-HCl) (pH 6.8), 10% glycerin (w/v), 2% sodium dodecyl sulfate (SDS) (w/v), and protease and phosphatase inhibitor cocktails (Sigma Aldrich, Gillingham, UK). Cell pellets were directly mixed with the lysis buffer. Proteins were transferred onto nitrocellulose membranes using semi-dry transfer. Primary antibodies for BAX, BCL-xL (Cell Signaling Technology, London, UK), BIM (Abcam, Cambridge, UK), BCL-2, MCL-1, BAK, and PUMA (Santa Cruz, Heidelberg, Germany) were used. A 1:5000 dilution was used for β-actin (Sigma Aldrich) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) (Abcam) proteins as loadings control. After overnight incubation, peroxidase-conjugated secondary antibodies (Chemicon, France or Jackson ImmunoResearch, Ely, UK) at a 1:5000 dilution were added. Images were captured using a LAS-3000 Imager equipped with a cooled 12-bit digital charge-coupled device camera (FUJIFILM UK Ltd. Systems, Sheffield, UK). For BCL-2 family protein quantification HeLa cell extracts (20 µg) were blotted alongside samples for standardisation and quantification. For densitometry analysis, the intensity of protein bands within tumour samples were normalised to the loading control after background subtraction using ImageJ software (National Institutes of Health, Bethesda, MD, USA). Absolute protein levels were calculated by normalising values to HeLa cells. Absolute protein concentrations in HeLa cells were determined with recombinant proteins previously [7].

4.2. DR_MOMP Calculation

Cells’ sensitivity to undergo apoptosis was calculated using the ordinary differential equation-based systems model DR_MOMP [7]. Absolute protein levels of BAK, BAX, BCL-2, BCL-xL, and MCL-1 were used to calculate the genotoxic stress dose that induces MOMP (‘stress dose’). ABT-199’s binding affinities were modelled with dissociation constants 0.01 nM for BCL-2, 48.00 nM for BCL-xL and 444.00 nM for MCL-1 [10,21]. Sensitivity of CRC cell lines and PDX models to undergo apoptosis was determined without and with different doses of ABT-199 ranging from 0 to 10 μM. The ABT-199 dose was applied at the same time and for the same duration as the modelled genotoxic stress (12 h). DR_MOMP was executed with Matlab R2007b (V.7.5.0.342; The MathWorks, Natick, MA, USA). Where p53 status was known (DLD-1 cells), DR_MOMP calculations accounted for deficient Puma and Noxa induction [7].

4.3. Chemicals

ABT-199 was kindly provided by AbbVie (North Chicago, IL, USA). Oxaliplatin (OX) and Fluorouracil (5-FU) were obtained from Sigma-Aldrich (St. Louis, MO, USA) when used for in vitro studies. OX (TEVA UK Ltd., Castleford, UK), 5-FU (ITH PHARMA, London, UK) and folinic acid (TEVA UK Ltd., Castleford, UK) were kindly provided by Beaumont Hospital Dublin when used for in vivo studies (FOLFOX treatment regimen).

4.4. Cell Lines

Authenticated human CRC (HCT-116 and DLD-1) cell lines were provided by the School of Medicine, Dentistry and Biomedical Sciences, Queen’s University Belfast, United Kingdom. Furthermore, all cell lines were Mycoplasma tested regularly. Cells were cultured in RPMI 1640 medium (Sigma-Aldrich). All cell lines were maintained at 37 °C in a humidified atmosphere supplied with 5% CO2.

4.5. In Vitro Toxicity Assays and Synergy Calculations

Acid phosphatase assay was used to measure cell viability based on the conversion of 4-Nitrophenyl phosphate (pNPP) to p-nitrophenol by cytosolic acid phosphatase [45]. pNPP disodium salt hexahydrate was purchased from Sigma-Aldrich. Cells were grown in a 96 well plate at a density of 1.5 × 104 cells per well and treated with increasing concentrations of 5-FU/OX (with a concentration ratio of 1:5 between OX and 5-FU) in combination with increasing concentrations of ABT-199. After 48 h treatment medium was removed and each well was washed once with 200 μL of 1× phosphate buffered saline (PBS). To each well, 100 μL of assay buffer (0.1 M sodium acetate at pH 5.0, 0.1% Triton X-100, and 7.25 mM p-nitrophenyl phosphate) was added. The plates were then incubated at 37 °C for 2 h. The reaction was finally stopped with the addition of 50 μL and colour development was assayed at 405 nm using a Clariostar plate reader (BMG Labtech, Offenburg, Germany), and applying a path length correction of 200 µL. The non-enzymatic hydrolysis of the pNPP substrate was also determined by including wells with the assay buffer and without any cells. An excel template was used to calculate the fraction affected from viability percentages and the results were analysed with the web version of Chalice Analyzer (Horizon Discovery, Waterbeach, UK) to calculate isobolograms.

4.6. Animals

Four- to six-week old (c. 20–25 g) female non-obese diabetic/severe combined immunodeficient (NOD/SCID) mice were purchased from Charles River Laboratories (Cambridge, UK). All animal procedures were approved by the Health Products Regulatory Authority (HPRA) (#AE18982-P099) and the University College Dublin Animal Research Ethics Committee (AREC) (#AREC-16-11).

4.7. PDX Efficacy Study

PDX CRC0076 and CRC0344 were selected for in vivo studies based on model predictions and known tumorigenicity of these models in vivo. The researcher performing the in vivo studies was blinded to the predicted DR_MOMP predictions. PDX implantation and expansion was performed as previously described [46].
Mice with established tumours (n = 13–14/group/PDX) average volume 360 mm3) were assigned into groups ensuring that an even distribution of tumour sizes was represented in each group and that group averages were equal. Groups were subsequently treated for 4 weeks with either (i) ABT-199 (oral gavage, once daily at 100 mg/kg, dissolved in 60% phosal, 50 propylene glycol (PG), 30% polyethylene glycol (PEG) 400 and 10% ethanol), (ii) FOLFOX (all constituent drugs delivered IP, once weekly (on day 3 of each cycle), 5-FU (40 mg/kg in PBS) + Folinic Acid (13.4 mg/kg in PBS) followed two hours later by OX (2.4 mg/kg in 5% glucose/water (v/v)), (iii) a combination of both regimens (ABT-199 + FOLFOX), or (iv) all diluent vehicle solutions, as described. Tumour size was evaluated three times per week by caliper and volume calculated using the formula Tumour Volume = 4/3π × (d/2)2 × D/2, where d is the minor tumour axis and D is the major tumour axis.

4.8. Immunohistochemistry (IHC)

Subcutaneous tumours were excised, rinsed twice in Dulbecco’s-(D)PBS and fixed in 4% formaldehyde for 48 h. Tissues were embedded in paraffin and 5 μm thick sections were cut. Routine hematoxylin and eosin (H&E) staining was performed to facilitate histological evaluation. Standard immunohistochemistry methods were used to detect cleaved caspase-3 on 3 tumours from each treatment group. Briefly, antigen retrieval was performed using pH 6.8 sodium citrate buffer. Primary antibodies (Rabbit polyclonal antibody to cleaved caspase-3, 1:200 dilution or Mouse IgG antibody to Ki-67, 1:150 dilution [both Cell Signaling Technology, London, UK]) were incubated for 1 hr at room temperature, and staining completed using the Lab Vision UltraVision Large Volume Detection System (Fisher Scientific, Dublin, Ireland). All sections were scored blindly and semi-quantitatively by two independent observers. The Allred Score Method based on proportion of positive cells (scored on a scale of 0–5) and staining intensity (scored on a scale of 0–3) was used for cleaved caspase-3. The proportion and intensity were then summed and averaged to produce a total score [15]. For Ki-67 analysis, slides were imaged and colour deconvolution was performed using ImageJ. The Ki-67 proliferation index was determined by calculating Ki-67 stained cells as a percentage of total cells.

4.9. 18F-FDG-PET/CT Study

Mice were fasted for 6 to 12 h prior to imaging and water provided ad libitum. 18F-FDGPET/CT scans was performed on a LabPET4 Triumph scanner (Trifoil Imaging, Chatsworth, CA, USA). PET system calibration was carried out by scanning a phantom scan of known activity as measured in a dose calibrator (Capintec, Florham, NJ, USA). A baseline scan was performed prior to treatment starting and a second scan was performed after 2 weeks of treatment. For imaging, mice were placed in a warming box for 20 min followed by intravenous (IV) tail vein injection under isoflurane anaesthesia, of 8.9 megabecquerels (MBq) (+/− 0.8 MBq) 18F-FDG diluted with saline to a volume of 200 µL. Following injection of 18F-FDG, mice remained under anaesthesia for a 1-h uptake time, with appropriate heating provided to maintain body temperature. A CT scan was taken immediately prior to a 15 min PET scan. Cone beam CT images (focal spot size, 33 mm) were acquired over ~5 min with 512 projections at a geometric magnification equal to 2. X-ray tube voltage and current were 50 peak kilovoltage (kVp) and 350 µA respectively and exposure was 232 ms. CT images were reconstructed with Cobra (Exxim Computing version 7.3.4, Pleasanton, CA, USA)) to a voxel size of 0.1 × 0.1 × 0.1 mm3. PET images were re-constructed with a 3-dimensional (3D) MLEM (Maximum Likelihood Expectation Maximisation) reconstruction algorithm (15 iterations, voxel size of 0.5 × 0.5 × 0.6 m3) using the LabPET software (version 1.14.0, Trifoil Imaging) PMOD (Version 3.208, PMOD Technologies, Zürich, Switzerland) software was used to draw VOIs around the entire tumour and SUV was calculated using the following equation: SUV = A/(ID/BW), where A is the activity in the tumour per unit volume, ID/BW is the injected dose per bodyweight and all activities are decay corrected to injection time. SUVmean is the mean SUV value of all voxels included in the VOI and is a measure of glucose uptake and metabolism (i.e., overall glucose metabolism) within the VOI, in this case in the PDX tumours. SUVmax is the average value of the four most intense voxels from each VOI. PET images were co-registered and resliced to CT data for illustration purposes or for guidance as required when delineating CT tumour volumes.

4.10. Radiomic Analysis of CRC PDX CT Images

Following reconstruction of image volumes, 45 tumours from n = 2 PDX models (CRC0334 and CRC0076) were delineated by three independent researchers. ROI (region of interest) were generated and segmentation carried out on CT scans of PDX tumours using a semi-automatic method using ITK-Snap [47] and PMOD software (Version 3.208, PMOD Technologies). Radiomic feature extraction was performed with PyRadiomics V2.7.7 (Harvard Medical School, Boston, MA, USA) [48]. A total of 833 features were extracted from each pre-treatment CRC PDX CT image, including first order, shape, Gray Level Size Zone Matrix (GLSZM), and Gray Level Run Length Matrix (GLRLM). The feature data set of the 45 pre-treatment CRC PDX images was next split into two cohorts, training and validation, in a 80:20 ratio (n = 36 and n = 9 respectively, ensuring a 50:50 spilt of the two PDXs (CRC0076 and CRC0344) in each cohort). Feature selection and model training were performed only on the training cohort to prevent information leakage. Features were selected in a three-step approach; All features related to shape (in PDX subcutaneous models, shape features are highly conserved) were removed from further analysis. Inter-correlated features (determined in the training dataset by pairwise spearman correlation (>0.85) and removing the feature with the highest correlation to the remaining features) were likewise removed, yielding 71 features. To further reduce feature dimensionality for this small training dataset a decision tree-based recursive feature elimination method was applied (Figure 5C) and recursive feature elimination was performed using the recursive feature elimination (RFE) function (caret package) with 3 repeats and a resampling method of repeated CV [49]. The performance of the selected number of features was evaluated on the validation dataset via AUC for the ROC. All statistical analysis, model building, and model evaluation were performed in R (version 3.5.2; http://www.r-project.org). The R code used to develop the radiomic classifier to separate PDX CRC0076 and PDX CRC0344 can be found in Supplementary Document 1.

4.11. Statistical Analyses

Analyses were performed using R (Version 3.4.4, The R Foundation, Vienna, Austria) unless otherwise stated. One-way analysis of variance (ANOVA) with Tukey’s ‘Honest Significant Difference’ (Tukey-HSD) post-hoc tests were used to prove for differences in mean values in vitro and for PUMA and BIM protein level, cleaved caspase-3 and Ki-67 IHC data. p-values were not adjusted for multiple comparison among multiple ANOVAs. For in vivo PET analysis, multiple paired t-tests were employed to compare raw PET values for each animal to its own pre-treatment measurement. p-values were adjusted for multiple comparison using the Holm method. SAS Software Version 9.3 (SAS Institute Inc., Cary, NC, USA) was employed for in vivo data management and its statistics, including linear mixed models to model CRC0076 and CRC0344 PDX tumour growth separately. Mouse identifier was included as a random effect and time effects were included as a continuous covariate, in addition to treatment, as a factor. Model diagnostics included evaluation of leverage, studentised residuals and best linear unbiased predictions of mouse effects. All pairwise treatment comparisons were performed and corrected for multiple testing (Bonferroni correction). The nominal 5% level of significance was used and therefore each of the 6 pairwise comparisons were evaluated at the 0.8% level. Throughout, p values ≤ 0.05 were considered statistically significant. For CT radiomic feature analysis, Wilcoxon rank sum test was applied to compare radiomic features between treatment groups. All radiomic feature statistical analysis was performed using R.

5. Conclusions

We demonstrate that the system model DR_MOMP may have potential to effectively predict CRC patients for whom the addition of ABT-199 to a FOLFOX regimen would be advantageous. Tumour SUVmean derived from 18F-FDG-PET imaging may represent an early response biomarker. Finally, to the best of our knowledge we have for the first time established a CT radiomics pipeline for application in CRC PDX studies. Such analyses will likely become increasingly important as the PDX field evolves, with an expanding translational focus on the use of orthotopic models supported by clinically relevant imaging. Our findings warrant further validation; firstly in an expanded PDX “population” cohort [12] employing a wider range of DR_MOMP stress-dose predictions. These analyses should be followed by prediction testing in an appropriate clinical setting.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/10/2978/s1, Table S1: BCL-2 protein profiles by quantitative Western Blotting, Table S2: ‘Stress dose’ calculated using DR_MOMP. Figure S1: Western blots used for quantification of BCL-2 family proteins for DR_MOMP calculations in DLDL-1 cells, Figure S2: A synergistic effect of ABT-199 and 5-FU/OX combination observed in CRC cell lines, Figure S3: Western blots used for quantification of BCL-2 family proteins for DR_MOMP calculations in CRC0076 and CRC-344, Figure S4: Full BIM and PUMA Western blots undertaken in post-mortem CRC0076 and CRC0344 tumours, Figure S5: Ki-67 analysis in post-mortem CRC0076 and CRC0344 tumours, Figure S6: Evaluation of radiomics features as early biomarkers of response to FOLFOX and ABT-199 combination therapy in PDXs CRC0076 and CRC0344. Supplementary Document 1: R code used to develop the radiomic classifier to separate PDX CRC0076 and PDX CRC0344.

Author Contributions

Conceptualization, A.T.B. and J.H.M.P.; Data curation, A.C.O., M.A.J., A.U.L. and P.D.; Formal analysis, A.C.O., M.A.J., A.U.L., S.C., E.C., I.S.M., K.C., L.S., F.L., M.M.V., P.D., S.A.K. and S.S.; Funding acquisition,W.M.G., L.T., A.T.B. and J.H.M.P.; Investigation, A.C.O., M.A.J., S.C., E.C., I.S.M., L.S., F.L., A.L., K.W. and M.M.V.; Methodology, A.T.B. and J.H.M.P.; Resources, E.R.Z., W.M.G., P.L., H.C.W., A.B., L.T., A.T.B. and J.H.M.P.; Supervision, W.M.G., A.B., L.T., A.T.B. and J.H.M.P.; Visualization, A.C.O., M.A.J., A.U.L., A.T.B. and J.H.M.P.; Writing—original draft, A.C.O., M.A.J., A.U.L., S.C., E.C., I.S.M., L.S., F.L., A.L., M.M.V., P.D., A.T.B. and J.H.M.P.; Writing—review & editing, A.C.O., M.A.J., A.U.L., I.S.M., K.C., E.R.Z., F.L., K.W., W.M.G., S.A.K., S.S., P.L., H.C.W., A.B., L.T., A.T.B. and J.H.M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a grant from the Health Research Board, Ireland (HRA-POR-2014-547) to J.H.M.P. and A.T.B., and by Investigator grants from Science Foundation Ireland (13/IA/1881; 14/IA/2582) to J.H.M.P. A.T.B. received funding from Science Foundation Ireland (SFI) Career Development Award ‘Coloforetell’ (grant agreement number: 13/CDA/2183). L.T., A.T.B. and J.H.M.P. are further supported by the European Union’s Horizon 2020 Health Research and Innovation award ‘Colossus’ (grant agreement number: 754923). A.T.B., I.S.M., L.T. and A.B. are members of the EurOPDX Consortium, and receive funding from the European Union’s Horizon 2020 research and innovation programme, grant agreement no. #731105 (EurOPDX Research Infrastructure, www.europdx.eu). W.M.G. and E.C. are supported by the Irish Cancer Society Collaborative Cancer Research Centre BREAST-PREDICT (CCRC13GAL), the Science Foundation Ireland Investigator Programme OPTi-PREDICT (grant code 15/IA/3104) and the Science Foundation Ireland Strategic Partnership Programme Precision Oncology Ireland (grant code 18/SPP/3522). A.T.B., I.S.M., W.M.G. and E.C. are further supported by, and this publication has emanated from, research supported in part by RCSI, UCD and NUIG pursuant to a grant from Science Foundation Ireland (SFI) under grant number 18/RI/5759 and a grant from Science Foundation Ireland (SFI) and the European Regional Development Fund (ERDF) under grant number 13/RC/2073.” H.C.W. is supported by the Dutch Cancer Society (KWF Kankerbestrijding), Project number 12085/2018–2 and TRANSCAN Joint Transnational Call 2016 (JTC2016 “CLEARLY”- no UM 2017–8295).

Acknowledgments

We thank AbbVie (Dublin, Ireland) for providing ABT-199 and M2i Limited (Dublin, Ireland) for supplying [18F]FDG. We thank members of Livio Trusolino’s group (IRCCS, University of Torino) who provided PDX models. Further, we thank our colleagues from the Centre for Cancer Research and Cell Biology (Queen’s University Belfast) for providing DLD-1 and HCT-116 cells used in this study. We thank Orna Bacon (Beaumont Hospital, Dublin) for providing FOLFOX.

Conflicts of Interest

H.C.W. has (minority) shares in the company Oncoradiomics. W.M.G. is a co-founder, shareholder (minority) and Chief Scientific Officer at OncoMark Limited. All other authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [Green Version]
  2. Kaufmann, S.H.; Earnshaw, W.C. Induction of apoptosis by cancer chemotherapy. Exp. Cell Res. 2000, 256, 42–49. [Google Scholar] [CrossRef]
  3. Strasser, A.; Cory, S.; Adams, J.M. Deciphering the rules of programmed cell death to improve therapy of cancer and other diseases. EMBO J. 2011, 30, 3667–3683. [Google Scholar] [CrossRef]
  4. Tait, S.W.; Green, D.R. Mitochondria and cell death: outer membrane permeabilization and beyond. Nat. Rev. Mol. Cell Biol. 2010, 11, 621–632. [Google Scholar] [CrossRef]
  5. Food and Drug Administration (FDA). Prescribing Information Reference ID 4275193, NDA 208573/S-004 and S-005. Available online: https://www.accessdata.fda.gov/drugsatfda_docs/appletter/2018/208573Orig1s004s005ltr.pdf (accessed on 1 July 2019).
  6. Vogler, M. Targeting BCL2-Proteins for the Treatment of Solid Tumours. Adv. Med. 2014, 2014, 943648. [Google Scholar] [CrossRef] [Green Version]
  7. Lindner, A.U.; Concannon, C.G.; Boukes, G.J.; Cannon, M.D.; Llambi, F.; Ryan, D.; Boland, K.; Kehoe, J.; McNamara, D.A.; Murray, F.; et al. Systems analysis of BCL2 protein family interactions establishes a model to predict responses to chemotherapy. Cancer Res. 2013, 73, 519–528. [Google Scholar] [CrossRef] [Green Version]
  8. Flanagan, L.; Lindner, A.U.; de Chaumont, C.; Kehoe, J.; Fay, J.; Bacon, O.; Toomey, S.; Huber, H.J.; Hennessy, B.T.; Kay, E.W.; et al. BCL2 protein signalling determines acute responses to neoadjuvant chemoradiotherapy in rectal cancer. J. Mol. Med. Berl. 2015, 93, 315–326. [Google Scholar] [CrossRef]
  9. Lindner, A.U.; Salvucci, M.; Morgan, C.; Monsefi, N.; Resler, A.J.; Cremona, M.; Curry, S.; Toomey, S.; O’Byrne, R.; Bacon, O.; et al. BCL-2 system analysis identifies high-risk colorectal cancer patients. Gut 2017, 66, 2141–2148. [Google Scholar] [CrossRef]
  10. Lucantoni, F.; Lindner, A.U.; O’Donovan, N.; Dussmann, H.; Prehn, J.H.M. Systems modeling accurately predicts responses to genotoxic agents and their synergism with BCL-2 inhibitors in triple negative breast cancer cells. Cell Death Dis. 2018, 9, 42. [Google Scholar] [CrossRef]
  11. Weiss, J.; Gajek, T.; Kohler, B.C.; Haefeli, W.E. Venetoclax (ABT-199) Might Act as a Perpetrator in Pharmacokinetic Drug-Drug Interactions. Pharmaceutics 2016, 8, 5. [Google Scholar] [CrossRef] [Green Version]
  12. Bertotti, A.; Migliardi, G.; Galimi, F.; Sassi, F.; Torti, D.; Isella, C.; Cora, D.; Di Nicolantonio, F.; Buscarino, M.; Petti, C.; et al. A molecularly annotated platform of patient-derived xenografts (“xenopatients”) identifies HER2 as an effective therapeutic target in cetuximab-resistant colorectal cancer. Cancer Discov. 2011, 1, 508–523. [Google Scholar] [CrossRef] [Green Version]
  13. Zanella, E.R.; Galimi, F.; Sassi, F.; Migliardi, G.; Cottino, F.; Leto, S.M.; Lupo, B.; Erriquez, J.; Isella, C.; Comoglio, P.M.; et al. IGF2 is an actionable target that identifies a distinct subpopulation of colorectal cancer patients with marginal response to anti-EGFR therapies. Sci. Transl. Med. 2015, 7, 272ra12. [Google Scholar] [CrossRef]
  14. Workman, P.; Aboagye, E.O.; Balkwill, F.; Balmain, A.; Bruder, G.; Chaplin, D.J.; Double, J.A.; Everitt, J.; Farningham, D.A.; Glennie, M.J.; et al. Guidelines for the welfare and use of animals in cancer research. Br. J. Cancer 2010, 102, 1555–1577. [Google Scholar] [CrossRef] [Green Version]
  15. Qureshi, A.; Pervez, S. Allred scoring for ER reporting and it’s impact in clearly distinguishing ER negative from ER positive breast cancers. J. Pak. Med. Assoc. 2010, 60, 350–353. [Google Scholar]
  16. Salvucci, M.; Wurstle, M.L.; Morgan, C.; Curry, S.; Cremona, M.; Lindner, A.U.; Bacon, O.; Resler, A.J.; Murphy, A.C.; O’Byrne, R.; et al. A Stepwise Integrated Approach to Personalized Risk Predictions in Stage III Colorectal Cancer. Clin. Cancer Res. 2017, 23, 1200–1212. [Google Scholar] [CrossRef] [Green Version]
  17. Hu, T.; Li, Z.; Gao, C.Y.; Cho, C.H. Mechanisms of drug resistance in colon cancer and its therapeutic strategies. World J. Gastroenterol. 2016, 22, 6876–6889. [Google Scholar] [CrossRef]
  18. Hector, S.; Prehn, J.H. Apoptosis signaling proteins as prognostic biomarkers in colorectal cancer: A review. Biochim. Biophys. Acta 2009, 1795, 117–129. [Google Scholar] [CrossRef]
  19. Violette, S.; Poulain, L.; Dussaulx, E.; Pepin, D.; Faussat, A.M.; Chambaz, J.; Lacorte, J.M.; Staedel, C.; Lesuffleur, T. Resistance of colon cancer cells to long-term 5-fluorouracil exposure is correlated to the relative level of Bcl-2 and Bcl-X(L) in addition to Bax and p53 status. Int. J. Cancer 2002, 98, 498–504. [Google Scholar] [CrossRef]
  20. Koehler, B.C.; Scherr, A.L.; Lorenz, S.; Urbanik, T.; Kautz, N.; Elssner, C.; Welte, S.; Bermejo, J.L.; Jager, D.; Schulze-Bergkamen, H. Beyond cell death—antiapoptotic Bcl-2 proteins regulate migration and invasion of colorectal cancer cells in vitro. PLoS ONE 2013, 8, e76446. [Google Scholar] [CrossRef]
  21. Souers, A.J.; Leverson, J.D.; Boghaert, E.R.; Ackler, S.L.; Catron, N.D.; Chen, J.; Dayton, B.D.; Ding, H.; Enschede, S.H.; Fairbrother, W.J.; et al. ABT-199, a potent and selective BCL-2 inhibitor, achieves antitumor activity while sparing platelets. Nat. Med. 2013, 19, 202–208. [Google Scholar] [CrossRef]
  22. Fricker, M.; O’Prey, J.; Tolkovsky, A.M.; Ryan, K.M. Phosphorylation of Puma modulates its apoptotic function by regulating protein stability. Cell Death Dis. 2010, 1, e59. [Google Scholar] [CrossRef] [Green Version]
  23. Meller, R.; Cameron, J.A.; Torrey, D.J.; Clayton, C.E.; Ordonez, A.N.; Henshall, D.C.; Minami, M.; Schindler, C.K.; Saugstad, J.A.; Simon, R.P. Rapid degradation of Bim by the ubiquitin-proteasome pathway mediates short-term ischemic tolerance in cultured neurons. J. Biol. Chem. 2006, 281, 7429–7436. [Google Scholar] [CrossRef] [Green Version]
  24. Rehm, M.; Dussmann, H.; Janicke, R.U.; Tavare, J.M.; Kogel, D.; Prehn, J.H. Single-cell fluorescence resonance energy transfer analysis demonstrates that caspase activation during apoptosis is a rapid process. Role of caspase-3. J. Biol. Chem. 2002, 277, 24506–24514. [Google Scholar] [CrossRef] [Green Version]
  25. Bock, F.J.; Tait, S.W.G. Mitochondria as multifaceted regulators of cell death. Nat. Rev. Mol. Cell Biol. 2020, 21, 85–100. [Google Scholar] [CrossRef]
  26. Lok, S.W.; Whittle, J.R.; Vaillant, F.; Teh, C.E.; Lo, L.L.; Policheni, A.N.; Bergin, A.R.T.; Desai, J.; Ftouni, S.; Gandolfo, L.C.; et al. A Phase Ib Dose-Escalation and Expansion Study of the BCL2 Inhibitor Venetoclax Combined with Tamoxifen in ER and BCL2-Positive Metastatic Breast Cancer. Cancer Discov. 2019, 9, 354–369. [Google Scholar] [CrossRef]
  27. Boland, K.; Flanagan, L.; Prehn, J.H. Paracrine control of tissue regeneration and cell proliferation by Caspase-3. Cell Death Dis. 2013, 4, e725. [Google Scholar] [CrossRef] [Green Version]
  28. Hardwick, J.M.; Soane, L. Multiple functions of BCL-2 family proteins. Cold Spring Harb. Perspect. Biol. 2013, 5. [Google Scholar] [CrossRef] [Green Version]
  29. Chen, Y.B.; Aon, M.A.; Hsu, Y.T.; Soane, L.; Teng, X.; McCaffery, J.M.; Cheng, W.C.; Qi, B.; Li, H.; Alavian, K.N.; et al. Bcl-xL regulates mitochondrial energetics by stabilizing the inner membrane potential. J. Cell Biol. 2011, 195, 263–276. [Google Scholar] [CrossRef]
  30. Williams, A.; Hayashi, T.; Wolozny, D.; Yin, B.; Su, T.C.; Betenbaugh, M.J.; Su, T.P. The non-apoptotic action of Bcl-xL: regulating Ca(2+) signaling and bioenergetics at the ER-mitochondrion interface. J. Bioenerg. Biomembr. 2016, 48, 211–225. [Google Scholar] [CrossRef]
  31. Gross, A. BCL-2 family proteins as regulators of mitochondria metabolism. Biochim. Biophys. Acta 2016, 1857, 1243–1246. [Google Scholar] [CrossRef]
  32. Lucantoni, F.; Dussmann, H.; Llorente-Folch, I.; Prehn, J.H.M. BCL2 and BCL(X)L selective inhibitors decrease mitochondrial ATP production in breast cancer cells and are synthetically lethal when combined with 2-deoxy-D-glucose. Oncotarget 2018, 9, 26046–26063. [Google Scholar] [CrossRef] [Green Version]
  33. Zhao, Y.; Coloff, J.L.; Ferguson, E.C.; Jacobs, S.R.; Cui, K.; Rathmell, J.C. Glucose metabolism attenuates p53 and Puma-dependent cell death upon growth factor deprivation. J. Biol. Chem. 2008, 283, 36344–36353. [Google Scholar] [CrossRef] [Green Version]
  34. Gross, A.; Katz, S.G. Non-apoptotic functions of BCL-2 family proteins. Cell Death Differ. 2017, 24, 1348–1358. [Google Scholar] [CrossRef]
  35. Rogers, W.; Thulasi Seetha, S.; Refaee, T.A.G.; Lieverse, R.I.Y.; Granzier, R.W.Y.; Ibrahim, A.; Keek, S.A.; Sanduleanu, S.; Primakov, S.P.; Beuque, M.P.L.; et al. Radiomics: from qualitative to quantitative imaging. Br. J. Radiol. 2020, 93, 20190948. [Google Scholar] [CrossRef]
  36. Gillies, R.J.; Kinahan, P.E.; Hricak, H. Radiomics: Images Are More than Pictures, They Are Data. Radiology 2016, 278, 563–577. [Google Scholar] [CrossRef] [Green Version]
  37. Lambin, P.; Rios-Velazquez, E.; Leijenaar, R.; Carvalho, S.; van Stiphout, R.G.; Granton, P.; Zegers, C.M.; Gillies, R.; Boellard, R.; Dekker, A.; et al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur. J. Cancer 2012, 48, 441–446. [Google Scholar] [CrossRef] [Green Version]
  38. Sanduleanu, S.; Woodruff, H.C.; de Jong, E.E.C.; van Timmeren, J.E.; Jochems, A.; Dubois, L.; Lambin, P. Tracking tumor biology with radiomics: A systematic review utilizing a radiomics quality score. Radiother. Oncol. 2018, 127, 349–360. [Google Scholar] [CrossRef]
  39. Chaddad, A.; Daniel, P.; Niazi, T. Radiomics Evaluation of Histological Heterogeneity Using Multiscale Textures Derived From 3D Wavelet Transformation of Multispectral Images. Front. Oncol. 2018, 8, 96. [Google Scholar] [CrossRef]
  40. Byrne, A.T.; Alferez, D.G.; Amant, F.; Annibali, D.; Arribas, J.; Biankin, A.V.; Bruna, A.; Budinska, E.; Caldas, C.; Chang, D.K.; et al. Interrogating open issues in cancer precision medicine with patient-derived xenografts. Nat. Rev. Cancer 2017, 17, 254–268. [Google Scholar] [CrossRef]
  41. Wu, G.; Woodruff, H.C.; Sanduleanu, S.; Refaee, T.; Jochems, A.; Leijenaar, R.; Gietema, H.; Shen, J.; Wang, R.; Xiong, J.; et al. Preoperative CT-based radiomics combined with intraoperative frozen section is predictive of invasive adenocarcinoma in pulmonary nodules: a multicenter study. Eur. Radiol. 2020, 30, 2680–2691. [Google Scholar] [CrossRef] [Green Version]
  42. Panth, K.M.; Leijenaar, R.T.; Carvalho, S.; Lieuwes, N.G.; Yaromina, A.; Dubois, L.; Lambin, P. Is there a causal relationship between genetic changes and radiomics-based image features? An in vivo preclinical experiment with doxycycline inducible GADD34 tumor cells. Radiother. Oncol. 2015, 116, 462–466. [Google Scholar] [CrossRef] [Green Version]
  43. Lambin, P.; Leijenaar, R.T.H.; Deist, T.M.; Peerlings, J.; de Jong, E.E.C.; van Timmeren, J.; Sanduleanu, S.; Larue, R.; Even, A.J.G.; Jochems, A.; et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat. Rev. Clin. Oncol. 2017, 14, 749–762. [Google Scholar] [CrossRef]
  44. Zinn, P.O.; Singh, S.K.; Kotrotsou, A.; Hassan, I.; Thomas, G.; Luedi, M.M.; Elakkad, A.; Elshafeey, N.; Idris, T.; Mosley, J.; et al. A Coclinical Radiogenomic Validation Study: Conserved Magnetic Resonance Radiomic Appearance of Periostin-Expressing Glioblastoma in Patients and Xenograft Models. Clin. Cancer Res. 2018, 24, 6288–6299. [Google Scholar] [CrossRef] [Green Version]
  45. Yang, T.T.; Sinai, P.; Kain, S.R. An acid phosphatase assay for quantifying the growth of adherent and nonadherent cells. Anal. Biochem. 1996, 241, 103–108. [Google Scholar] [CrossRef]
  46. Galimi, F.; Torti, D.; Sassi, F.; Isella, C.; Cora, D.; Gastaldi, S.; Ribero, D.; Muratore, A.; Massucco, P.; Siatis, D.; et al. Genetic and expression analysis of MET, MACC1, and HGF in metastatic colorectal cancer: Response to met inhibition in patient xenografts and pathologic correlations. Clin. Cancer Res. 2011, 17, 3146–3156. [Google Scholar] [CrossRef] [Green Version]
  47. Yushkevich, P.A.; Piven, J.; Hazlett, H.C.; Smith, R.G.; Ho, S.; Gee, J.C.; Gerig, G. User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. Neuroimage 2006, 31, 1116–1128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Van Griethuysen, J.J.M.; Fedorov, A.; Parmar, C.; Hosny, A.; Aucoin, N.; Narayan, V.; Beets-Tan, R.G.H.; Fillion-Robin, J.C.; Pieper, S.; Aerts, H. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer Res. 2017, 77, e104–e107. [Google Scholar] [CrossRef] [Green Version]
  49. Kuhn, M. Building Predictive Models in R Using the caret Package. J. Stat. Softw. 2008, 28, 26. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure 1. Schematic representation of the experimental workflow. The sensitivity of tumour cells to undergo apoptosis was calculated after quantitative protein profiling using the deterministic systems model DR_MOMP [7]. First, protein levels were calculated by quantitative Western blotting in proteins extracted from colorectal cancer (CRC) patient-derived xenograft (PDX) tumours, using cell lines as standards in which absolute BCL-2 protein levels were previously determined. Protein profiles then served as input into the deterministic systems model, DR_MOMP. The model calculates a ‘stress dose’ (pro-apoptotic ‘BH3 only protein dose’) required to induce mitochondrial outer membrane permeabilization (MOMP), the process leading to mitochondrial apoptosis. Kinetics of the BCL-2 antagonist ABT-199 were included in the DR_MOMP model [10] to determine sensitisation of the apoptosis signalling network by the BCL-2 antagonist.
Figure 1. Schematic representation of the experimental workflow. The sensitivity of tumour cells to undergo apoptosis was calculated after quantitative protein profiling using the deterministic systems model DR_MOMP [7]. First, protein levels were calculated by quantitative Western blotting in proteins extracted from colorectal cancer (CRC) patient-derived xenograft (PDX) tumours, using cell lines as standards in which absolute BCL-2 protein levels were previously determined. Protein profiles then served as input into the deterministic systems model, DR_MOMP. The model calculates a ‘stress dose’ (pro-apoptotic ‘BH3 only protein dose’) required to induce mitochondrial outer membrane permeabilization (MOMP), the process leading to mitochondrial apoptosis. Kinetics of the BCL-2 antagonist ABT-199 were included in the DR_MOMP model [10] to determine sensitisation of the apoptosis signalling network by the BCL-2 antagonist.
Cancers 12 02978 g001
Figure 2. DR_MOMP successfully predicated the sensitivity of two CRC PDXs to treatment with FOLFOX and ABT-199 in combination. (A) Two mCRC PDXs (CRC0344 and CRC0076) were assessed for BCL-2 family proteins levels by quantitative Western blot (Figure S3). (B) The susceptibility of each model to undergo MOMP in response to genotoxic chemotherapy alone or in combination treatment with ABT-199 was calculated using DR_MOMP. (C,D) Anti-tumour activity of ABT-199 alone and in combination with FOLFOX chemotherapy in the two models was assessed in vivo. Mice bearing either a CRC0076 or CRC0344 tumour were treated with vehicle control, ABT-199 100 mg/kg, FOLFOX or the combination ABT-199 + FOLFOX; n = 13–14/group, error bars = s.e.m., * p ≤ 0.05. The effect of treatment on absolute tumour volumes of the predicted combination only responder CRC0076 and the predicted FOLFOX alone responder CRC0344 PDX models are shown respectively (C,D).
Figure 2. DR_MOMP successfully predicated the sensitivity of two CRC PDXs to treatment with FOLFOX and ABT-199 in combination. (A) Two mCRC PDXs (CRC0344 and CRC0076) were assessed for BCL-2 family proteins levels by quantitative Western blot (Figure S3). (B) The susceptibility of each model to undergo MOMP in response to genotoxic chemotherapy alone or in combination treatment with ABT-199 was calculated using DR_MOMP. (C,D) Anti-tumour activity of ABT-199 alone and in combination with FOLFOX chemotherapy in the two models was assessed in vivo. Mice bearing either a CRC0076 or CRC0344 tumour were treated with vehicle control, ABT-199 100 mg/kg, FOLFOX or the combination ABT-199 + FOLFOX; n = 13–14/group, error bars = s.e.m., * p ≤ 0.05. The effect of treatment on absolute tumour volumes of the predicted combination only responder CRC0076 and the predicted FOLFOX alone responder CRC0344 PDX models are shown respectively (C,D).
Cancers 12 02978 g002
Figure 3. Decrease in p53 upregulated modulator of apoptosis (PUMA) levels but no effect on BCL2-Interacting Mediator of cell death (BIM) or cleaved caspase-3, following ABT-199 treatment in CRC0076 PDX model. (A,B) Western blot analysis comparing the levels of PUMA and BIM in CRC0076 (A) and CRC0344 (B) PDX models. β-actin or glyceraldehyde 3-phosphate dehydrogenase (GAPDH) (see Figure S4) were used as loading controls. Error bars represent mean +/− 2 SD, * p ≤ 0.05 (ANOVA, Tukey post-hoc). (C,D) Analysis of cleaved caspase-3 levels in CRC0076 (C) and CRC0344 (D) PDX models. Representative images of H&E (magnification 10×, scale bar 100 μm) and cleaved caspase-3 (magnification 5×, scale bar 200 μm) are shown for CRC0076 and CRC0344 in (C,D) respectively. The Allred method [15] was used to quantify the extent and intensity of staining, error bars represent mean +/− s.e.m (n = 3 tumours per treatment, n = 2–3 sections per tumour).
Figure 3. Decrease in p53 upregulated modulator of apoptosis (PUMA) levels but no effect on BCL2-Interacting Mediator of cell death (BIM) or cleaved caspase-3, following ABT-199 treatment in CRC0076 PDX model. (A,B) Western blot analysis comparing the levels of PUMA and BIM in CRC0076 (A) and CRC0344 (B) PDX models. β-actin or glyceraldehyde 3-phosphate dehydrogenase (GAPDH) (see Figure S4) were used as loading controls. Error bars represent mean +/− 2 SD, * p ≤ 0.05 (ANOVA, Tukey post-hoc). (C,D) Analysis of cleaved caspase-3 levels in CRC0076 (C) and CRC0344 (D) PDX models. Representative images of H&E (magnification 10×, scale bar 100 μm) and cleaved caspase-3 (magnification 5×, scale bar 200 μm) are shown for CRC0076 and CRC0344 in (C,D) respectively. The Allred method [15] was used to quantify the extent and intensity of staining, error bars represent mean +/− s.e.m (n = 3 tumours per treatment, n = 2–3 sections per tumour).
Cancers 12 02978 g003
Figure 4. Analysis of glucose uptake and metabolism in the predicted combination responder CRC0076 and the predicted FOLFOX alone responder CRC0344 PDX models. For imaging, n = 6 mice/group/PDX (at pre-treatment time point) received 18F-FDG via tail vein injection, with image acquisition taking place 1 h after injection. Representative fused PET/CT images of CRC0076 and CRC0344 PDXs are shown (the same mice at pre- and on-treatment time points). White arrows indicate the tumours, scale bar indicates standard uptake values (SUVs) (A). N.B. Animals were imaged in a prone position. The coronal plane shown is positioned at the centre of the tumour which is located posterior to the bladder. Where the coronal plane intersects with the bladder, urinary excretion of FDG is observed as an intense signal. Comparison of relative change SUVmean at week 2 of treatment normalised to pre-treatment for CRC0076 and CRC0344 PDXs is presented. Error bars = s.e.m., * p ≤ 0.05 (paired t-tests, adjusted Holms) (B).
Figure 4. Analysis of glucose uptake and metabolism in the predicted combination responder CRC0076 and the predicted FOLFOX alone responder CRC0344 PDX models. For imaging, n = 6 mice/group/PDX (at pre-treatment time point) received 18F-FDG via tail vein injection, with image acquisition taking place 1 h after injection. Representative fused PET/CT images of CRC0076 and CRC0344 PDXs are shown (the same mice at pre- and on-treatment time points). White arrows indicate the tumours, scale bar indicates standard uptake values (SUVs) (A). N.B. Animals were imaged in a prone position. The coronal plane shown is positioned at the centre of the tumour which is located posterior to the bladder. Where the coronal plane intersects with the bladder, urinary excretion of FDG is observed as an intense signal. Comparison of relative change SUVmean at week 2 of treatment normalised to pre-treatment for CRC0076 and CRC0344 PDXs is presented. Error bars = s.e.m., * p ≤ 0.05 (paired t-tests, adjusted Holms) (B).
Cancers 12 02978 g004
Figure 5. Radiomic analysis of pre-treatment CT images of CRC0076 and CRC0344 can successfully distinguish between the PDXs using seven unique radiographic features. (A) Visual Representation of PDX radiomic analysis workflow. n = 45 baseline CT scans were partitioned into a training dataset (n = 36) and validation dataset (n = 9) before Recursive Feature Elimination was performed on the training dataset. (B) Representative CT axial slice of a mouse implanted with CRC PDX CRC0344. Region of interest (ROI, in red) was drawn semi-automatically on every slice containing tumour to enable segmentations of the tumour for radiomic feature extraction. (C) 3D representation of the mask used to segment the tumour. (D) Output of the recursive feature elimination (RFE) method. Using an RFE algorithm the ideal number of features was identified as 7 for radiomic model construction. (E): Graph depicting the result of the receiver operating characteristic (ROC) area under the curve (AUC) analysis of the 7 features identified by the model in a test cohort (n = 9). AUC is 0.8 (95% Confidence interval 0.408–1).
Figure 5. Radiomic analysis of pre-treatment CT images of CRC0076 and CRC0344 can successfully distinguish between the PDXs using seven unique radiographic features. (A) Visual Representation of PDX radiomic analysis workflow. n = 45 baseline CT scans were partitioned into a training dataset (n = 36) and validation dataset (n = 9) before Recursive Feature Elimination was performed on the training dataset. (B) Representative CT axial slice of a mouse implanted with CRC PDX CRC0344. Region of interest (ROI, in red) was drawn semi-automatically on every slice containing tumour to enable segmentations of the tumour for radiomic feature extraction. (C) 3D representation of the mask used to segment the tumour. (D) Output of the recursive feature elimination (RFE) method. Using an RFE algorithm the ideal number of features was identified as 7 for radiomic model construction. (E): Graph depicting the result of the receiver operating characteristic (ROC) area under the curve (AUC) analysis of the 7 features identified by the model in a test cohort (n = 9). AUC is 0.8 (95% Confidence interval 0.408–1).
Cancers 12 02978 g005

Share and Cite

MDPI and ACS Style

O’Farrell, A.C.; Jarzabek, M.A.; Lindner, A.U.; Carberry, S.; Conroy, E.; Miller, I.S.; Connor, K.; Shiels, L.; Zanella, E.R.; Lucantoni, F.; et al. Implementing Systems Modelling and Molecular Imaging to Predict the Efficacy of BCL-2 Inhibition in Colorectal Cancer Patient-Derived Xenograft Models. Cancers 2020, 12, 2978. https://doi.org/10.3390/cancers12102978

AMA Style

O’Farrell AC, Jarzabek MA, Lindner AU, Carberry S, Conroy E, Miller IS, Connor K, Shiels L, Zanella ER, Lucantoni F, et al. Implementing Systems Modelling and Molecular Imaging to Predict the Efficacy of BCL-2 Inhibition in Colorectal Cancer Patient-Derived Xenograft Models. Cancers. 2020; 12(10):2978. https://doi.org/10.3390/cancers12102978

Chicago/Turabian Style

O’Farrell, Alice C., Monika A. Jarzabek, Andreas U. Lindner, Steven Carberry, Emer Conroy, Ian S. Miller, Kate Connor, Liam Shiels, Eugenia R. Zanella, Federico Lucantoni, and et al. 2020. "Implementing Systems Modelling and Molecular Imaging to Predict the Efficacy of BCL-2 Inhibition in Colorectal Cancer Patient-Derived Xenograft Models" Cancers 12, no. 10: 2978. https://doi.org/10.3390/cancers12102978

APA Style

O’Farrell, A. C., Jarzabek, M. A., Lindner, A. U., Carberry, S., Conroy, E., Miller, I. S., Connor, K., Shiels, L., Zanella, E. R., Lucantoni, F., Lafferty, A., White, K., Meyer Villamandos, M., Dicker, P., Gallagher, W. M., Keek, S. A., Sanduleanu, S., Lambin, P., Woodruff, H. C., ... Prehn, J. H. M. (2020). Implementing Systems Modelling and Molecular Imaging to Predict the Efficacy of BCL-2 Inhibition in Colorectal Cancer Patient-Derived Xenograft Models. Cancers, 12(10), 2978. https://doi.org/10.3390/cancers12102978

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop