Background Immunogenic cell death (ICD) is a tumor cell death involving both innate and adaptive immune responses. Given published findings that oxaliplatin, but not irinotecan, drives ICD, we investigated whether single nucleotide polymorphisms (SNPs) in the ICD pathway are associated with the efficacy of oxaliplatin-based chemotherapy in metastatic colorectal cancer (mCRC).
Methods Two randomized clinical trials data were analyzed: discovery cohort, FOLFOX/bevacizumab arm (MAVERICC); validation cohort, FOLFOXIRI/bevacizumab arm (TRIBE); and two control cohorts, FOLFIRI/bevacizumab arms (both trials). Genomic DNA extracted from blood samples was genotyped. Ten SNPs in the ICD pathway were tested for associations with clinical outcomes.
Results In total, 648 patients were included. In the discovery cohort, three SNPs were significantly associated with clinical outcomes in univariate analysis: CALR rs1010222 with progression-free survival (G/G vs any A, HR=0.61, 95% CI 0.43–0.88), ANXA1 rs1050305 with overall survival (OS) (A/A vs any G, HR=1.87, 95% CI 1.04–3.35), and LRP1 rs1799986 with OS (C/C vs any T, HR=1.69, 95% CI 1.07–2.70). Multivariate analysis confirmed the trend, but statistical significance was not reached. In the validation cohort, ANXA1 rs1050305, and LRP1 rs1799986 were validated to have the significant associations with clinical outcome. No significant associations of these SNPs were observed in the two control cohorts. Treatment-by-SNP interaction test confirmed the predictive values.
Conclusions The predictive utility of ICD-related SNPs for the efficacy of oxaliplatin-based chemotherapy was demonstrated, warranting further validation studies to be translated into personalized treatment strategies using conventional cytotoxic agents in mCRC.
- gastrointestinal neoplasms
- adaptive immunity
- translational medical research
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
The current standard of care for patients with microsatellite-stable metastatic colorectal cancer (mCRC) consists of the combination of a cytotoxic chemotherapy backbone and a monoclonal antibody based on RAS/BRAF testing.1 2 In clinical practice, there are two choices for cytotoxic agent, namely oxaliplatin or irinotecan in combination with 5-fluorouracil, and two for monoclonal antibodies targeting vascular endothelial growth factor (VEGF) or epidermal growth factor receptor (EGFR). Most research has focused on the efficacy of monoclonal antibodies, leading to the successful identification of predictive and prognostic biomarkers for EGFR inhibitors such as RAS and BRAF mutation status and primary tumor sidedness.3–6 However, there are no predictive biomarkers for the selection of oxaliplatin or irinotecan which have distinct mechanisms of action: oxaliplatin is a platinum compound forming intra-strand and inter-strand cross-links; irinotecan is a topoisomerase I (TOP1) inhibitor trapping TOP1 cleavage complex with DNA replication and transcription. The choice of these two drugs remains an important clinical question for optimizing treatment in individual patients.
Immunogenic cell death (ICD) is a novel concept of tumor cell death that engages both innate and adaptive immune responses, thereby conferring an additional immunogenic antitumor effect to cytotoxic agents.7–9 However, most chemotherapeutics kill tumor cells in a non-immunogenic manner, and only a few cytotoxic agents can prime antitumor immunity by inducing ICD.10 11 The mechanisms of ICD induction are characterized by the ability to induce endoplasmic reticulum stress, which in turn leads to the release of damage-associated molecular patterns (DAMPs) from dying tumor cells and the subsequent activation of pattern recognition receptors (PRRs) of host innate immune cells such as dendritic cells.12 The engagement of innate immune responses cross-primes adaptive responses controlled by cytotoxic T lymphocytes through proinflammatory cytokine release, potentially leading to the continuous elimination of tumor cells.12 13 Oxaliplatin has been recently well-recognized as an ICD inducer, similarly as other drugs such as anthracyclines and bortezomib.7 In a preclinical study, CRC cells treated with oxaliplatin exhibited increased pre-apoptotic exposure of calreticulin (CALR) and post-apoptotic release of high-mobility group box 1 (HMGB1), which are the signals of DAMPs required for ICD induction.14 By contrast, irinotecan lacks the mechanism of action that drives ICD, and thus, its cytotoxic effects depend on non-immunogenic mechanisms.10 15 These findings suggested that genetic variants in the ICD pathway are related to the efficacy of oxaliplatin but not irinotecan, prompting us to test whether these novel biomarkers will be useful for patient selection. Therefore, we investigated the predictive utility of single nucleotide polymorphisms (SNPs) in ICD-related genes encoding essential DAMPs and PRRs using genetic and clinical data from two first-line randomized clinical trials comparing oxaliplatin- (or oxaliplatin-containing) and irinotecan-based treatments in patients with mCRC.
Patient population and study design
This study included patients with mCRC enrolled in two first-line randomized trials: MAVERICC (NCT01374425)16 and TRIBE (NCT00719797).17 In the MAVERICC trial, patients were randomized to treatment with either FOLFOX plus bevacizumab or FOLFIRI plus bevacizumab. In the TRIBE trial, patients were randomized to the FOLFOXIRI plus bevacizumab or FOLFIRI plus bevacizumab arm. Patients without sufficient peripheral whole blood samples for analyzes were excluded from our study. To assess whether the ICD pathway is specifically related to the efficacy of oxaliplatin, we set the cohorts as follows: (1) discovery cohort, FOLFOX plus bevacizumab arm in MAVERICC; (2) validation cohort, FOLFOXIRI plus bevacizumab arm in TRIBE; (3) control cohort 1, FOLFIRI plus bevacizumab arm in MAVERICC; and (4) control cohort 2, FOLFIRI plus bevacizumab arm in TRIBE (figure 1). All patients provided informed consent for molecular research prior to study enrollment.
Genotyping and selecting polymorphisms
Genomic DNA was extracted from peripheral whole blood collected before treatment initiation using a QIAmp Kit (Qiagen, Valencia, California, USA) in accordance with the manufacturer’s protocol (www.qiagen.com). The OncoArray of 530 K SNPs was used for genotyping (Illumina, San Diego, California, USA).18 After genotyping, we planned this exploratory analysis focusing on five genes encoding DAMPs or PRRs related to chemotherapy-induced ICD,12 including three DAMP-encoding genes (CALR, HMGB1, and ANXA1) and two PRR-encoding genes (LRP1 and P2R×7) (online supplementary table S1). Although FPR1, TLR3, and TLR4 are also key PRR-encoding genes, these genes were not included in this study because they were not predictive of oxaliplatin efficacy in patients with CRC in another large biomarker study.19 The candidate SNPs for this study were arbitrarily selected according to the following criteria: (1) minor allele frequency in Caucasians of at least 10% in the Ensemble Genome Browser (https://www.ensembl.org), (2) having potential biological functions based on public databases (https://snpinfo.niehs.nih.gov; https://www.ncbi.nlm.nih.gov), and (3) tag SNPs chosen from HapMap genotype data with an r2 threshold of 0.8 (https://snpinfo.niehs.nih.gov). In total, 10 SNPs were selected, as presented in online supplementary table S2.
Selected SNPs were evaluated for their associations with tumor response, progression-free survival (PFS), and overall survival (OS). The overall response rate (ORR) was calculated as the percentage of patients with either a complete or partial response using the Response Evaluation Criteria in Solid Tumors version 1.1. PFS was defined as the time from randomization to disease progression or death from any cause. OS was defined as the time from randomization to death from any cause. Patients who experienced no events were censored at the last follow-up date. The correlation between each SNP and ORR was examined using the χ2 test. To test the association between each SNP and PFS or OS, univariate and multivariate analyzes using the Cox proportional hazards regression model and the log-rank test were performed. In the multivariate analyzes, the following study-specific adjusted covariates were used: the TRIBE trial included sex, age, Eastern Cooperative Oncology Group (ECOG) performance status (PS), primary tumor site, liver-limited disease, adjuvant chemotherapy, BRAF status, and RAS status; and the MAVERICC trial included ethnicity, sex, age, ECOG PS, primary tumor site, primary tumor resected, number of metastases, and KRAS status. SNPs were coded using a dominant genetic model for the number of variant alleles. To formally assess the predictive value, the treatment-by-SNP interaction was tested within each trial using the HR calculated in the multivariate analyzes. All analyzes were two-sided at a significance level of 0.05 and were performed using SAS V.9.4 software.
In total, 648 patients were included in this study (figure 1). Some characteristics were unbalanced between the trial cohorts. In particular ECOG PS 0, left-sided primary tumor, metastasis in ≤2 organs, and primary tumor resection were more prevalent in the TRIBE trial arms than in the MAVERICC trial arms (online supplementary table S3). The median follow-up and survival time in each cohort are summarized in online supplementary table S4.
Predictive value of ICD-related SNPs
All tested associations between each SNP and clinical outcomes are presented in online supplemental tables S5–8. Although several significant associations were observed in the discovery or validation cohort, few associations were confirmed in control cohort 1 or 2.
In particular, three SNPs displayed significant associations with outcome in the univariate analysis in the discovery cohort (table 1, figure 2). First, patients carrying any A allele of CALR rs1010222 had a better PFS than those with a G/G genotype (HR=0.61, 95% CI 0.43–0.88, p=0.008). Second, patients carrying any G allele of ANXA1 rs1050305 had worse OS than those with an A/A genotype (HR=1.87, 95% CI 1.04–3.35, p=0.03). Finally, patients carrying any T allele of LRP1 rs1799986 had worse OS than those with a C/C genotype (HR=1.69, 95% CI 1.07–2.70, p=0.03). Multivariate analysis confirmed these trends, but none of the associations reached statistical significance. In the validation cohort, two SNPs were confirmed to have significant associations with clinical outcome (table 1, figure 2). Namely, patients carrying any G allele of ANXA1 rs1050305 had worse ORR, PFS (revealed in univariate analysis), and OS (revealed in both univariate and multivariate analyzes) than those with an A/A genotype; and patients carrying any T allele of LRP1 rs1799986 had worse PFS (revealed in multivariate analysis) than those with a C/C genotype, even though this clinical endpoint was not same as that in the discovery cohort. Of note, in the two control cohorts, no significant associations were observed between these SNPs (CALR rs1010222, ANXA1 rs1050305, and LRP1 rs1799986) and clinical outcomes (table 1, figure 2).
The results of the treatment-by-SNP interaction test are presented in table 2. In total, six SNPs exhibited significant interactions with treatment: HMGB1 rs1360485 (PFS in TRIBE), ANXA1 rs1050305 (OS in TRIBE), LRP1 rs1799986 (PFS in TRIBE), LRP1 rs11172113 (both PFS and OS in TRIBE), P2R×7 rs208294 (OS in MAVERICC and OS in TRIBE), and P2R×7 rs1718119 (OS in MAVERICC).
This is the first study to identify genetic variants in ICD-related pathways that may predict the efficacy of oxaliplatin-based first-line chemotherapy in patients with mCRC. No predictive biomarkers for the personalized use of cytotoxic agents in mCRC have been identified so far. DNA damage response genes and the transcriptome-based consensus molecular subtype have been examined in biomarker research, but their predictive value has not been validated.16 20 Gray et al reported that genetic variants of three PRR genes (FPR1, TLR3, and TLR4) were not related to the efficacy of oxaliplatin in patients with CRC in the adjuvant setting (using data from the SCOT trials) and the first-line setting for patients with metastasis (using data from the COIN/COIN-B trials).19 However, DAMPs are also essential tumor-sided factors for ICD induction, whereas PRRs are host-side factors. Therefore, in our study, we extensively assessed the functional variants of both DAMPs and PRRs. According to our hypothesis, ICD-related SNPs were significantly associated with the efficacy of oxaliplatin-based combination. These findings imply the linkage between ICD pathway and oxaliplatin efficacy, which can be coupled with the potential synergy of oxaliplatin and immunotherapy.21 22
Our findings identified ANXA1 rs1050305 and LRP1 rs1799986 as SNPs exhibiting significant associations with the efficacy of oxaliplatin but not irinotecan as confirmed in the discovery and validation cohorts. ANXA1 is one of the DAMPs released from the cytoplasm of dying cells, and it interacts with FPR1, which is expressed on dendritic cells.23 The extracellular accumulation of ANXA1 constitutes a hallmark of ICD.12 ANXA1 rs1050305 is a synonymous SNP located at an exonic splicing enhancer/silencer site that disrupts splicing activity and causes alternative splicing (https://snpinfo.niehs.nih.gov/snpinfo/snpfunc.html). LRP1 is a PRR activated via docking of cell surface-exposed CALR, which is another hallmark of chemotherapy-induced ICD.12 24 LRP1 rs1799986 is a missense SNP that may participate in splicing regulation. This SNP is reported to be associated with susceptibility to Alzheimer’s disease, supporting its functional significance.25 Of note, while the predictive value of ANXA1 rs1050305 was consistently validated on OS, that of LRP1 rs1799986 was observed on inconsistent clinical outcomes between the discovery and validation cohorts: the discovery cohort showed a significant association with OS, and the validation cohort showed that with PFS. This discrepancy might be due to some unidentified factors such as different contribution of second-line treatment between two cohorts, which might affect the associations between biomarker and clinical outcome.
Our study had great strengths in its design, featuring discovery, validation, and control cohorts obtained from two randomized clinical trials. The validated effects of the SNPs on clinical outcomes in the discovery and validation cohorts, lacking of that in the two control cohorts, strongly supports the predictive value of ANXA1 rs1050305 and LRP1 rs1799986 specifically for oxaliplatin efficacy. The multiple significant associations observed in the discovery and validation cohorts but not in the control cohorts further indicate that the ICD pathway plays an important role in the efficacy of oxaliplatin but not irinotecan. Another strength was that the treatment-by-SNP interaction test confirmed the predictive value of ANXA1 rs1050305 and LRP1 rs1799986 and highlighted the predictive utility of 6 out of 10 tested SNPs.26
This study had two key limitations. First, this was a prospective–retrospective study; thus, the results need to be validated in prospective clinical trials. According to published guidelines, one or more validation studies are required to reach a sufficient level of evidence for medical utility.27 28 Second, the validation cohort consisted of patients treated with FOLFOXIRI plus bevacizumab, which is both an oxaliplatin-containing and irinotecan-containing regimen. To validate the different effects of SNPs on oxaliplatin and irinotecan efficacy, it will be better to use a treatment arm containing only oxaliplatin-based regimens as the validation cohort. However, our study provides reliable evidence describing the specific connection of ICD-related SNPs with the efficacy of “oxaliplatin-containing” treatment based on the use of “oxaliplatin-non-containing” cohorts as controls.
In conclusion, our study provides the first evidence that genetic variants in the ICD pathway are significantly associated with the efficacy of oxaliplatin-based, but not irinotecan-based, first-line chemotherapy in patients with mCRC. Further translational studies are warranted to validate our novel findings and realize the personalized use of conventional cytotoxic chemotherapeutic agents.
We thank all patients who contributed to this study.
Contributors HA and H-JL primarily planned, designed, and drafted the manuscript. YX and JM did statistical analyses. All authors made substantial contributions to data collection and drafting the manuscript. All authors read and approved the final manuscript.
Funding This work was supported by the National Cancer Institute [P30CA 014089 to HJL], Gloria Borges WunderGlo Foundation, Dhont Family Foundation, Victoria and Philip Wilson Research Fund, San Pedro Peninsula Cancer Guild, and Daniel Butler Research Fund.
Competing interests None declared.
Patient consent for publication Not required.
Ethics approval All patients provided informed consent for molecular research prior to study enrollment. The study protocol was approved by the Institutional Review Board of each participating institution and was conducted in accordance with the tenets of the Declaration of Helsinki as well as the Good Clinical Practice and REMARK guidelines.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data are available upon reasonable request. De-identified datasets analyzed in the current study are available from the corresponding author on reasonable request.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.