Abstract
Missense mutations in the polymerase epsilon (POLE) gene have been reported to generate proofreading defects resulting in an ultramutated genome and to sensitize tumors to checkpoint blockade immunotherapy. However, many POLE-mutated tumors do not respond to such treatment. To better understand the link between POLE mutation variants and response to immunotherapy, we prospectively assessed the efficacy of nivolumab in a multicenter clinical trial in patients bearing advanced mismatch repair–proficient POLE-mutated solid tumors. We found that only tumors harboring selective POLE pathogenic mutations in the DNA binding or catalytic site of the exonuclease domain presented high mutational burden with a specific single-base substitution signature, high T-cell infiltrates, and a high response rate to anti–PD-1 monotherapy. This study illustrates how specific DNA repair defects sensitize to immunotherapy. POLE proofreading deficiency represents a novel agnostic biomarker for response to PD-1 checkpoint blockade therapy.
POLE proofreading deficiency leads to high tumor mutational burden with high tumor-infiltrating lymphocytes and predicts anti–PD-1 efficacy in mismatch repair–proficient tumors. Conversely, tumors harboring POLE mutations not affecting proofreading derived no benefit from PD-1 blockade. POLE proofreading deficiency is a new tissue-agnostic biomarker for cancer immunotherapy.
See related video: https://vimeo.com/720727355
This article is highlighted in the In This Issue feature, p. 1397
INTRODUCTION
Polymerase epsilon (POLE) is a major eukaryotic DNA polymerase involved in both DNA replication and repair (1). Its proofreading domain, harboring exonuclease activity, deletes mispaired bases inserted during DNA replication. Mutations that alter the exonuclease domain of POLE and inactivate its editing function lead to the extreme accumulation of missense and nonsense mutations with a very high tumor mutational burden (TMB) with as many as 100 mutations per megabase pair (TMB > 100 mt/Mb; refs. 1–3) with a specific POLE-related mutational single-base substitution (COSMIC SBS10a and 10b) signature >50% (4, 5). A high TMB has been associated with improved outcomes with immune-checkpoint blockade (6, 7), especially in tumors associated with known carcinogens (e.g., tobacco, UV light exposure) and in tumors with DNA repair defects such as mismatch repair deficiency (MMRd; ref. 8).
Defective proofreading POLE (POLE-pd) mutations are rare and predominantly described in mismatch repair–proficient (MMRp) colorectal and endometrial cancers as sporadic and germline events (2, 9). In general, POLE-pd tumors have a good prognosis (2, 10), and as a result, lower prevalence of POLE-pd mutations is expected in the advanced setting. Proofreading defects have been documented for multiple hotspots in the exonuclease domain (D275, P286, S297, N363K, F367, V411, L424, P436, M444, A456, Y458, S459, S461, and A463; refs. 2, 11). However, the functional impact of non-hotspot missense mutations across the POLE gene [i.e., variants of unknown significance (VUS)] in terms of the hypermutant phenotype and response to immunotherapy is incompletely characterized (12, 13).
To define the characteristics of POLE variants across solid tumors and their impact on the TMB and immunogenomic phenotype, we studied genomic databases across solid tumors and initiated a clinical trial following a two-stage Bayesian enrichment design testing the role of PD-1 blockade in POLE-mutated tumors. Cases were identified using a nationwide (France) prospective screening protocol for POLE variants, and proofreading deficiency was assessed by an ad hoc centralized molecular tumor board. Nivolumab (240 mg i.v. once every 2 weeks) was administered until disease progression [progressive disease (PD)], until toxicity, or up to 2 years. The primary endpoint was the objective response rate (ORR) assessed by RECIST v1.1 at 12 weeks. A secondary endpoint was to assess if POLE-pd led to improved outcomes.
RESULTS
Immunogenomic Analyses of POLE Mutations across Solid Tumors
According to the Memorial Sloan Kettering Cancer Center (MSKCC) Integrated Mutation Profiling of Actionable Cancer Targets (IMPACT) molecular screening program, POLE alterations are observed in 3.4% of all solid tumors (Fig. 1A), but POLE-pd variants account for only 0.4% of mutations and 0.1% in the advanced setting. Most POLE variants are benign, VUS, or located outside the exonuclease domain. These POLE VUS or nonpathogenic mutations are most often observed in high-TMB tumors such as melanoma, lung, or microsatellite instability–high (MSI-H) tumors, and likely represent passenger mutations rather than oncogenic drivers (Supplementary Figs. S1A–S1D and S2A–S2D; Supplementary Tables S1–S3). Crossing The Cancer Genome Atlas (TCGA) and IMPACT databases, we confirmed that POLE-pd mutations are mostly observed in endometrial cancer and colorectal cancer but also found in other gastrointestinal tract, glial, urothelial, prostate, and gynecologic cancers (Fig. 1B; Supplementary Tables S1 and S2). Patients with POLE-pd tumors are younger than other patients (Supplementary Table S3), and in gastrointestinal cancers POLE-pd is mainly observed in male patients (Supplementary Table S2C). Occasionally, we observed that MSI-H tumors could present concomitant POLE-pd, with alternative hotspots and amino acid substitutions compared with microsatellite-stable (MSS) tumors (Fig. 1C; Supplementary Fig. S2C and S2D). In MSS tumors, the most common variants are P286R and V411L. POLE-pd is associated with an ultramutated phenotype in both MMRd/MSI-H and MMRp/MSS tumors, with a mutation burden significantly higher than other tumors harboring other POLE alterations (Fig. 1D; Supplementary Tables S1 and S3). Only one POLE-pd variant, L424V, usually germline and associated with polyposis (14), was observed in a low TMB tumor (Fig. 1D). Analysis of samples with multiple POLE alterations showed that all POLE-pd variants were monoallelic without evidence of loss-of-heterozygosity (N = 96). Allelic frequency (AF) of POLE-pd variants showed a positive correlation with TMB, whereas no correlation was observed for nonpathogenic variants (Supplementary Fig. S1C).
To study the specifics of POLE-pd tumors, Pan-Cancer TCGA RNA-sequencing (RNA-seq) analyses of immune infiltrates and related pathways were performed comparing tumors harboring POLE-pd variants or other POLE alterations, accounting for TMB and microsatellite status (Fig. 1E; Supplementary Table S4). While POLE-pd tumors did not display higher cytotoxic T-cell infiltrate than other POLE-mutated tumors, they had the highest Th1 infiltrate and higher PD-1 (PDCD1) expression compared with other MSS tumors. PD-L1 (CD274) expression was decreased compared with MSS TMB-high tumors, a subgroup composed mainly of melanoma and lung cancers.
To further assess the immune landscape of MSS POLE-pd tumors, we investigated the RNA-seq data of the colorectal TCGA cohort, comparing with MSS non-POLE tumors and MSI-H tumors (Fig. 1F; Supplementary Table S5). Our analyses showed that POLE-pd tumors have a high cytotoxic Th1 T-cell infiltrate in line with higher PD-1 (PDCD1) expression. The POLE-pd tumor microenvironment also demonstrates strong upregulation of IFNG, PD-L1 (CD274), and CTLA-4 compared with MSS tumors, which resembles that seen in MSI-H tumors. Altogether, these data suggest that MSS POLE-pd tumors are good candidates for immune-checkpoint blockade.
Nivolumab Efficacy in Patients with Advanced Solid Tumors Harboring POLE Mutations
From January 2018 to December 2020, 61 patients harboring POLE mutations had been screened nationwide and assessed by the ad hoc molecular board (Fig. 2A). Twenty-one patients were enrolled in cohort 6 of Acsé Nivolumab, dedicated to patients harboring POLE mutations and excluding patients with MSI-H tumors or any MSS solid tumor types for which immunotherapy was currently available, such as melanoma and lung cancers. As TMB is not an approved biomarker for immunotherapy in Europe and is not routinely assessed, the molecular board assessed pathogenicity of POLE variants blinded to any TMB results. Twenty patients received at least one dose of nivolumab and were further considered for the efficacy analysis.
Clinical characteristics of included patients are summarized in Fig. 2B, and POLE variants are detailed in Fig. 2C. The median follow-up was 13.1 months (range = 0.5–26.1 months).
An ad hoc committee confirmed the pathogenicity of the POLE variants in the tumors of 12 (57%) patients, including three patients with P286R, three with V411L, two with N363K, two with S459F, one with S297Y, and one with A463V mutations, then classified as POLE-pd (Fig. 2C; Supplementary Table S6). Of note, one N363K variant was germline. Four patients had VUS, and five had POLE proofreading–proficient (POLE-pp) variants (Fig. 2B and C).
Figure 2D and E report, respectively, best response and RECIST target changes at each radiologic assessment. The trial achieved its primary objective with a 38% (7/19) ORR at 12 weeks, which was comprised of seven partial responses (Fig. 2F). Two patients had further complete response as best responses. The 12-week disease control rate (DCR) was 58% (11/19), with four patients maintaining stable disease. One additional patient had stable disease after initial radiologic progression. The safety profile of nivolumab was in accordance with safety data reported in other tumor types (Supplementary Table S7), and one patient stopped nivolumab for serious adverse events.
Responses were observed exclusively in patients with MMRp/MSS colorectal cancer (N = 5) and endometrial cancer (N = 2). In the overall cohort, median progression-free survival (mPFS) and median overall survival (mOS) were 5.4 months (Fig. 2G) and not reached, respectively (Fig. 2H).
Figure 2I displays a nearly complete response observed in a 42-year-old patient with an unresectable transverse colon tumor with a POLE-pd P286R variant. This response allowed resection of the residual lesion and revealed a 4.2-cm-long, pan-parietal mass with a mucoid aspect (Supplementary Fig. S3A and S3B). The pretreatment tumor (Supplementary Fig. S4A–S4D) was rich in tumor-infiltrating lymphocytes (TIL) and CD8+ cells and had low PD-L1 expression. The posttreatment lesion displayed mucoid pools without tumor cells, dissociated by trabeculae of paucicellular fibrosis, and was enriched in CD8+ cells (Supplementary Fig. S4E–S4H). There were also foci of infarctoid necrosis bordered by macrophagic cells (Supplementary Fig. S4I–S4L).
Activity of Nivolumab According to POLE Proofreading Defect
We investigated the impact of the prospectively determined POLE proofreading defect on nivolumab benefit. Whole-exome sequencing (WES) was successfully performed on 14 tumors, providing TMB and SBS signatures for pathogenicity confirmation. Two more patients had TMB results thanks to panel sequencing performed before inclusion in the trial.
Figures 2D and 3 display the overall results stratified by pathogenicity, TMB, and SBS POLE signatures. Sequencing confirmed that all POLE-pd variants were monoallelic, whereas other nonpathogenic mutations could be observed simultaneously.
TMB assessment confirmed that the pathogenicity assessed by the molecular board was concordant in 100% (10/10) of POLE-pd instances (Figs. 2D and 3A), as all displayed high TMB (median of 114 mt/Mb, min: 25, max: 385). In the same line of evidence, 100% (3/3) of the POLE-pp tumors had low TMB (median of 5 mt/Mb, min: 4, max: 9). Only two tumors with VUS were successfully sequenced, with both displaying low TMB (median of 3 mt/Mb, min: 2, max: 4); of note the patient with the V464A variant had low cellularity possibly impairing TMB assessment.
Study of the POLE SBS signature showed that all POLE-pd tumors had a typical POLE mutational profile (100%, 8/8; Figs. 2D and 3B). Among these POLE-pd tumors, we observed a wide variation in the proportion of POLE SBS signature (median 67%, min: 18%, max: 80%). Interestingly, two tumors with POLE-pp variants showed an attenuated POLE SBS signature below 15% (R47W and R446Q). The variants G6R, K425R, and V464A had a POLE SBS signature of 0% but surprisingly displayed MMRd SBS signatures >30%, whereas TMB was low, without any argument for MMRd (IHC and/or MSI score).
Then, using RNA-seq, we assessed the immune infiltrates according to proofreading deficiency. We observed that patients with POLE-pd/VUS variants had higher CD8+ T-cell infiltrates (Supplementary Fig. S5A and S5B) compared with POLE-pp variants. Of note, two patients with VUS (G330R and V464A) and long-lasting control had tumors with high CD8+ T-cell infiltrate.
Among assessable patients with POLE-pd mutations, five patients showed responses, for an ORR of 46% (5/11), and an additional 27% (3/11) of patients achieved stable disease, which yielded a DCR of 73%. Best responses were observed early with long-lasting control (Figs. 2E and 3A). Progressions were observed as best response in patients with endometrial cancer and glioblastoma despite having POLE-pd variants. Concerning the four patients with VUS, two had responses, one had stable disease, and one progressed. In contrast, at 12 weeks, all five patients with POLE-pp variants progressed. According to the Bayesian enrichment design, because no evidence of activity was observed in patients harboring POLE-pp mutations after the first 15 included patients, we subsequently enrolled only patients with POLE-pd variants or VUS.
Higher TMB was associated with higher likelihood of responses (Fig. 3B). The intensity of the SBS POLE signature was able to discriminate responders, stabilization, and progression (mean of 80%, 60%, and 15%, respectively). Correction of TMB by the amount of SBS POLE signature revealed that responders had the strongest POLE-pd phenotypes with extremely high TMB composed mostly of POLE-related mutations. The study of the AF of POLE-pd variants did not show a correlation with response or SBS signature (Supplementary Fig. S6A–S6E), whereas a trend for lower TMB was observed for lower AFs, and an SBS POLE signature <50% was observed only for low allelic variants (Supplementary Fig. S6C). All POLE-pd variants were monoallelic.
Figure 3C and E display PFS and OS results according to TMB and molecular board–assessed pathogenicity. Both TMB and pathogenicity were informative to predict the benefit, and the VUS category seems to present an intermediate benefit. When comparing nonpathogenic and pathogenic/VUS groups, the respective mPFS and mOS were 2.3 versus 9.6 months [HR (PFS) = 0.2; 95% confidence interval (CI) = 0.1–0.7] and 4.7 months versus not reached [HR (OS) = 0.1; 95% CI = 0.02–0.7; Fig. 3D and E). Figure 3F details best responses per variants, pathogenicity class, and related median TMB.
As the main POLE-pd mutations were observed in colorectal cancer, we inquired as to the outcomes of POLE-pd patients compared with MMRd/MSI-H patients in an independent cohort of advanced colorectal cancer exposed to single-agent PD-1 inhibitor (N = 35; Supplementary Fig. S7A–S7C). The observed PFS and OS of POLE-pd MMRp patients when exposed to anti–PD-1 was similar to those of MMRd/MSI-H colorectal cancer patients, confirming the benefit observed in this rare population.
In Silico POLE Mutation Pathogenicity Analyses Reveal Two Tridimensional Hotspot Sites Correlating with Benefit from Nivolumab
To assess the functional impact of POLE variants on proofreading activity, in silico tridimensional analyses combined with DNA affinity models were performed (Fig. 4A; Supplementary Figs. S8 and S9; Supplementary Video S1 and Supplementary Table S8).
The mutations in the exonuclease domain can be spatially and functionally grouped into three categories: (i) residues on the DNA binding side of the exonuclease domain (Fig. 4A, part d); (ii) residues affecting the catalytic core of the exonuclease domain (Fig. 4A, part e); and (iii) surface residues not predicted to hinder exonuclease function (Fig. 4A, part f). Classifying the variants according to their predicted tridimensional site (DNA binding site, catalytic site, or predicted to have a neutral impact as away from the DNA binding or catalytic site; Supplementary Table S7), we assessed the related TMB, POLE SBS signature, and associated benefit derived from nivolumab (Fig. 4B–D). Patients with variants within the catalytic site or the DNA binding site seemed to have a similar benefit (Fig. 4B and C), including two patients with VUS G330R and V464A and long-term benefits who were both recategorized as belonging to the catalytic site thanks to the in silico model. In contrast, in the neutral site group, only one response was observed (ORR = 17%), in one patient with a D435G VUS and a short benefit (Fig. 3A).
TMB was high for all variants within the DNA binding and catalytic sites, whereas TMB was low for all the neutral variants (Fig. 4D; Supplementary Fig. S1D). The POLE SBS signature had similar amplitudes for the DNA binding and catalytic site (Fig. 4D). Excellent PFS and OS were observed for patients with tumors harboring DNA binding and catalytic sites variants (Fig. 4E and F), with long-lasting plateaus. To the contrary, patients with tumors having neutral variants showed no sustained benefit.
DISCUSSION
Here we report the first prospective assessment of anti–PD-1 immunotherapy for patients with advanced tumors harboring predefined POLE mutations. While the primary endpoint was achieved in the overall population, responses were limited to patients with variants affecting proofreading, clustering in the DNA binding or catalytic sites. Patients were included with blinded TMB values, and similar to MMRd/MSI-H, POLE-pd status alone was sufficient to predict high TMB and benefit from anti–PD-1 therapy (15, 16). Strikingly, tumor types that exhibit high rates of MMRd/MSI-H, including colorectal cancer and endometrial cancer, also exhibit the highest relative prevalence of POLE-pd mutations, suggesting that intrinsic DNA repair impairment is a strong biomarker of sensitivity to checkpoint blockade in specific cancer types.
POLE-pd prevalence across solid tumors is extremely low, raising concerns on how to identify them in clinical practice. Large sequencing panels and TMB assessments are not available in clinical routine uniformly worldwide. Our study shows that some tumor types should be considered for POLE testing and provides characterization of variants impairing proofreading. Like MMRd/MSI-H tumors, metastatic MSS samples show low prevalence of POLE-pd because of an excellent prognosis at the early stage. In addition, this trial was supported by the French National Cancer Institute (INCa), which organized and provided guidelines for screening to all national sequencing platforms. Our study highlights how centralized coordinated efforts could help to identify patients with rare molecular alterations in the era of precision oncology.
Our Pan-Cancer analyses, combined with a clinical trial designed to include only patients with tumor types not eligible otherwise for immunotherapy, helped us to better understand how POLE mutations and TMB are intertwined. First, we confirmed that POLE-pd variants causing DNA repair impairment and high POLE-related TMB were exclusively monoallelic missense variants clustering in DNA binding and catalytic domains. Second, we showed that neutral missense mutations could be secondary to other hypermutability processes, as observed in UV light– or tobacco-driven tumors, displaying usually high sensitivity to immune-checkpoint blockade in line with their TMB and PD-L1 expression. Finally, we showed prospectively for the first time that tumors with POLE neutral missense variants in classically low-TMB tumors derived no benefit from PD-1 blockade, providing nuance to previous retrospective reports (12, 17). Altogether, these results highlight the need for careful joint assessment of POLE pathogenicity and TMB if available to avoid unnecessary and potentially harmful treatment by immunotherapy.
Tridimensional assessment of POLE variants seems to provide additional predictive value. Indeed, PFS and OS of patients with POLE DNA binding or catalytic sites mutated tumors were significantly improved compared with patients with tumors mutated in POLE away from these tridimensional sites. Our results also suggest that patients with tumors harboring missense VUS within POLE DNA binding and catalytic sites may derive a benefit from immunotherapy, whereas TMB and mutational signature assessments did not confirm their pathogenicity. Responsive patients with VUS in the Acsé trial could correspond to rare POLE variants with an attenuated phenotype. This hypothesis remains to be confirmed in a larger cohort and in vivo functional studies.
Translational analyses showed that POLE-pd tumors present with a strong T-cell infiltrate, immune-checkpoint expression, and high IFNγ. In this study, we also confirmed that advanced POLE-pd tumors presented with high CD8+ T-cell infiltrates compared with POLE-pp variants. Interestingly, patients with tumors harboring POLE VUS within the DNA binding/catalytic site also had a high CD8+ T-cell infiltrate.
WES allowed us to study the relationship between response to PD-1 blockade and molecular features of POLE-mutated tumors. TMB and pathogenicity assessed by the molecular board were perfectly congruent. All patients benefiting from nivolumab had tumors with high TMB, except from one with low cellularity. Inversely, two patients (19%) with high TMB had progression as best response. COSMIC POLE-related SBS signature assessment revealed that responses correlated with higher proportion of POLE-related SBS mutations, measuring the intensity of the proofreading defect phenotype. The most benefiting patients had high TMB with more than 60% of their mutations attributable to POLE proofreading defect. The study of AFs of POLE-pd variants showed that response could be observed if the POLE variant was clonal (AF ≥10%) or subclonal (AF <10%). Subclonal cases tended to have lower TMB and lower SBS POLE signature, suggesting that late emergence of a POLE-pd variant could lead to an attenuated phenotype. These data are in accordance with recent reports that specific mutational processes are highly immunogenic (18), confirming that the intensity of the POLE proofreading defect drives immunogenicity not only in MMRd/MSI-H tumors (19) but also in a pure MMRp/MSS context.
Altogether, this study showed that POLE-pd mutations predict high activity from anti–PD-1 in MMRp/MSS solid tumors and should be considered as a new tumor-agnostic biomarker for PD1-blockade immunotherapy. POLE-pd mutations lead to TIL-rich ultramutated tumors highly benefiting from immunotherapy. To the contrary, MMRp/MSS tumors with POLE mutations not affecting proofreading derived no benefit in tumor types that are classically considered poorly responsive to PD-1 blockade. This trial supports the idea to extend the concept of an agnostic predictive biomarker for checkpoint blockade of MSI to a larger set of tumors with intrinsic DNA repair impairment, MMRd/MSI-H, and/or POLE-pd.
METHODS
Acsé Nivolumab
Acsé Nivolumab is a multicenter (n = 46 sites) phase II basket trial sponsored by INCa (www.e-cancer.fr) and the French network of comprehensive cancer centers (Unicancer; www.unicancer.fr) investigating in multiple single-arm cohorts the efficacy of anti–PD-1 immunotherapy in rare tumor types. Version 1.1 from January 30, 2017, of the protocol was approved by the Agence Nationale de Sécurité du Médicament on February 10, 2017 (EUDRACT #2016-002257-37), was approved by the ethical committee “Comité de Protection des personnes (CPP) Ouest II (Angers)” on February 16, 2017 (#CPP 2017/02), and declared on ClinicalTrials.gov on January 6, 2017 (NCT03012581). Written informed consent was obtained from all the patients. The POLE cohort of the Acsé Nivolumab trial allowed the inclusion of patients with solid tumors harboring a missense POLE mutation, according to local sequencing methods, prospectively categorized for pathogenicity by an ad hoc centralized molecular board as POLE-pd, non–POLE-pp, or VUS. The main inclusion criteria were Eastern Cooperative Oncology Group Performance Status ≤1, advanced MMRp/MSS tumors confirmed by local IHC and/or PCR, tumors refractory to standard therapies, without prior exposure to immune-checkpoint inhibitors and/or without a primary tumor having an active approval for immunotherapy. Nivolumab (240 mg i.v. once every 2 weeks) was administered until disease progression (PD), until toxicity, or up to 2 years.
The primary endpoint was ORR assessed by RECIST v1.1 at 12 weeks by an independent review committee (IRC). Secondary endpoints included tolerance (Common Terminology Criteria for Adverse Events v4.0), best response, PFS, OS, and subgroup analyses according to proofreading defect.
The provision of baseline tumor material was mandatory per protocol at enrollment. Therefore, all patients included in the Acsé nivolumab trial consented to provide blood and either provide archival tumor material (less than 6 months old) or underwent a de novo tumor biopsy during the screening period for subsequent ancillary analysis. Pretreatment frozen peripheral blood mononuclear cells (PBMC; germinal material) and formalin-fixed, paraffin-embedded (FFPE) tumor blocks (somatic material) were submitted for DNA and RNA extraction, followed by WES and RNA-seq.
Design
The trial used a two-stage Bayesian enrichment design, which allowed for a smaller more informative trial that is specifically tied to decision-making within a drug development program. This process allowed for incremental changes to each cohort based on current data rather than restricting revisions in a trial design using fixed sample sizes (20, 21).
The first stage led to an interim analysis of the primary endpoint (ORR at 12 weeks), which was planned after 10 patients were enrolled in the cohort and completed the first scheduled disease evaluation. Based on an expected minimal ORR of 10% (85% confidence), the cohort continued only if at least one patient presented a radiologic response out of 10. For the second stage, intracohort analyses were performed after every five additional patients accrued (20). These analyses used all data accumulated up to that point and considered potential predictive markers of efficacy such as POLE pathogenicity assessment. Subset analyses allowed identification of subcohorts of poor outcomes that may be proposed for enrollment discontinuation, whereas those with good outcomes could enrich the design (22).
POLE Variant Classification
In the Pan-Cancer analyses, POLE pathogenicity was classified according to Campbell and colleagues (2) and enriched by other well-described pathogenic variants from the literature (Supplementary References). Considered pathogenic missense variants, whatever the amino acid substitution, were D275, P286, S297, N363, F367, V411, L424, P436, M444, A456, Y458, S459, S461, and A463.
Variants observed in the Acsé trial were classified prospectively according to the American College of Medical Genetics and Genomics standards and guidelines for the interpretation of sequence variants (23). A detailed description of the variants included in the trial and their classification is given in Supplementary Table S6.
DNA and RNA Extraction
Genomic DNA was extracted from FFPE tissues using the Maxwell 16 FFPE Plus LEV DNA Purification Kit (Promega) and from blood samples using the QIAamp Mini Kit (Qiagen) following the manufacturers’ instructions. DNA concentration was dosed by a Qubit fluorometer (Invitrogen), absorbances were analyzed by a NanoDrop spectrophotometer (Thermo Fisher Scientific), and quality profiles were assessed on Tapestation (Agilent Technologies).
Total RNA from FFPE tissues was extracted using FormaPure RNA (Beckman Coulter). RNA concentration and absorbances were analyzed with a NanoDrop spectrophotometer (Thermo Fisher Scientific), and quality profiles were assessed on Tapestation (Agilent Technologies).
WES and RNA-seq
Whole exome was captured from a FFPE, and paired constitutional DNA and libraries were prepared with Agilent SureSelect XTHS2 All Exon v8 (42 Mb) reagent kits according to the manufacturer's protocols (Agilent Technologies). Libraries were sequenced on a NovaSeq 6000 Illumina sequencer in 100-bp paired-end reads (Illumina).
For RNA-seq, libraries were prepared with the TruSeq Exome Kit following recommendations (Illumina) and sequenced on a NovaSeq 6000 sequencer in 75-bp paired-end reads.
WES Bioinformatics Analysis
All exome sequencing data were mapped to the hg38 human reference genome with bwa aligner v0.7.15-r1140 (https://github.com/lh3/bwa). Alignments were sorted and duplicates were marked using biobambam v2.0.79 (https://gitlab.com/german.tischler/biobambam2).
Somatic point mutations and short indels were called using the Mutect2 module from GATK v4.1.2 (https://github.com/broadinstitute/gatk). The genome aggregation database gnomAD v2.1.1 was used as a resource to calibrate variant calling model in Mutect2. Also, a panel of 80 normal samples sequenced on the same sequencing machine (NovaSeq6000) was used to remove recurring sequencing artifacts. All variants were annotated using the ensembl-vep v98.3 annotation tool, which integrates COSMIC v89 and dbSNP v152 resources in its cache.
TMB was computed from nonsynonymous variants (point mutations or indels) on the exome target with a variant allele frequency of at least 5% and is reported in mt/Mb.
SBS10 mutation signatures (SBS 10a + 10b corresponding to POLE proofreading defect) and MMRd signatures (4) were computed using SigProfilerSingleSample (COSMIC mutational signature V3; https://github.com/AlexandrovLab/SigProfilerSingleSample).
Translational Analyses and Data Source
Clinicogenomic data were extracted from the Pan-Cancer TCGA and the MSK-IMPACT databases on cBioPortal [Institutional Review Board (IRB) Protocol 12-245; refs. 26–28]. Microsatellite status was determined in both cohorts using the MSI sensor score, considering MSI-H with a score ≥10 and MSS with a score <10. For TCGA, the thresholds of 10 mt/Mb and 100 mt/Mb were used to define hypermutated tumors and ultramutated tumors, respectively, as previously described (2). Bulk RNA-seq data from the TCGA series were downloaded for the COAD/READ cohorts and for available samples harboring POLE alterations across all tumor types.
RNA-seq Bioinformatics Analysis and Immune Deconvolution of Bulk RNA-seq Data
Alignments were performed using STAR (version 2.7.10a) on the hg38 human reference genome. FastQC (version 0.11.9) and MultiQC (version 1.12) were used for quality control.
Analyses were carried out using R system software, version 4.0.3. TCGA expression matrices corresponded to unnormalized counts aggregated by HUGO gene symbol. We normalized these data using an original R script implementing the TPM algorithm (29). We then log2-transformed the TPM normalized data after adding 1e−5 to these data; in the resulting log2 TPM-normalized matrices, we set to zero the negative values. As TPM uses the length of the genes, we used the biomaRt R package to get this information from the Ensembl Homo sapiens gene database (version 105). We then performed immune cell–type deconvolution of the log2 TPM data, according to three algorithms: MCPcounter (30), implemented in the MCPcounter R package (https://github.com/ebecht/MCPcounter); XCELL (31); and EPIC [ref. 32; both xCELL and EPIC implemented in the immunedeconv R package (https://github.com/icbi-lab/immunedeconv)]. Markers related to immune multicellular structures (e.g., tertiary lymphoid structures) or immune functions (e.g., immunosuppression) were obtained from published signatures (33, 34).
Statistical Analyses
Analysis of efficacy endpoints was performed using the IRC assessment of response according to RECIST v1.1. According to parameters, statistical data are expressed in the form of frequencies and percentages (qualitative variables), mean and standard deviation, median and interquartile range (quantitative variables), and survival curve estimates calculated with the Kaplan–Meier method. In translational analyses, statistical tests are mentioned in the legends of the figures and included Fisher exact test, Kruskal–Wallis rank sum test, and Pearson Chi-squared test. Each primary and secondary evaluation criterion is reported with a CI of 95% unless otherwise mentioned. The Cox model was used to assess the HR when comparing subgroups for PFS and OS. Data have been analyzed with the SAS (SAS Inc.) and R (https://www.R-project.org/) software packages.
Concerning the statistical methodology for comparison of PanCancer and COAD/READ groups (Fig. 1E and F), we used one-way ordinary ANOVA tests allowing multiple comparison of MSS POLE pathogenic category with other categories, corrected for multiplicity using Dunnet test. Analyses were performed on MCPcounter, EPIC, or XCELL output for each parameter of interest using GraphPad software version 9. The level of significance displayed in the figures is as follows: ns, nonsignificant, >0.05; *, ≥0.05 and >0.01; **, ≤0.01 and >0.001; ***, ≤0.001 and >0.0001; ****, ≤0.0001. Similar statistical methodology was used for comparisons for Figs. 1D, 3B, and 4D and E.
Tumor IHC
An IHC study using antibodies to CD3, CD8, and PD-L1 was performed on one patient both in pretreatment biopsies and a surgical specimen (35). For PD-L1, slides were incubated with the clone QR1 (Diagomics) on unstained tissue sections (4 μm) using the Leica Bond III platform automated stainer according to the manufacturers’ recommendations. Positive PD-L1 expression was defined as any membranous staining in either tumor cells or in infiltrating inflammatory cells. The percentage of tumor cells demonstrating membranous PD-L1 staining was also scored in 5% increments. T cells in the tumor-associated immune infiltrates were assessed by semiquantitative estimation of CD3-positive (CD3: Thermo Scientific, clone SP7, 1:300 dilution) and CD8-positive (CD8: DAKO, clone C8/144B, 1:30 dilution) cell densities. They were scored as follows: 0, no positive cells; 1, scattered positive cells; 2, moderate number of positive cells; and 3, abundant occurrence of positive cells. In each tumor, the density of T-cell infiltration was examined at the invasive front (cells localized in stroma adjacent to the invasive tumor margin) and in the intratumoral compartment.
MSK Cohort of Patients with Colorectal Cancer Exposed to Immune-Checkpoint Inhibitors
The characteristics of the MSK colorectal cancer cohort have been presented before. In brief, clinical and genetic data from patients followed at MSKCC with colorectal cancer who underwent testing from 2015 through 2020 with the next-generation sequencing targeted MSK-IMPACT panel of somatic mutations were obtained from an internal database (26–28). Patients who underwent MSK-IMPACT testing and had metastatic or unresectable and locally advanced (n = 140) colon or rectal adenocarcinoma and who received at least one immune-checkpoint inhibitor between April 11, 2014, and December 2, 2019, were identified. Clinical–demographic data and cancer characteristics (histology, primary tumor location, staging) as well as prior treatment details were extracted from the electronic medical record and verified by licensed physicians (S.B. Maron and B.H. Diplas). MMRd status and POLE status were determined as reported in the Methods. The study was approved by the MSKCC IRB (2020-013).
Data Availability
A part of the analyses that support the findings of this study are based upon data generated by the TCGA Research Network (https://www.cancer.gov/tcga) and are openly available at https://portal.gdc.cancer.gov/projects.
Other reported data from the Acsé POLE clinical trials and MSK cohort of patients with colorectal cancer exposed to immune-checkpoint inhibitors are available by the corresponding authors upon reasonable request. The data are not publicly available due to information that could compromise the privacy of research participants.
Authors’ Disclosures
B. Rousseau reports nonfinancial support from FoundationOne, grants and personal fees from Neophore, other support from Bristol Myers Squibb and MSKCC, and grants from Swim Across America and the Nuovo-Soldati Foundation during the conduct of the study; personal fees from Bayer, Roche, Servier, and Gilead outside the submitted work; and a patent for methods and composition for cancer immunotherapy issued to MSKCC. R. Cohen reports personal fees and nonfinancial support from Bristol Myers Squibb and MSD Oncology during the conduct of the study, as well as personal fees from Amgen and Pierre Fabre, nonfinancial support from Mylan Medical, and grants from Servier Institute outside the submitted work. V. Simmet reports personal fees from Lilly, Servier, and Sanofi outside the submitted work. P. Augereau reports personal fees and nonfinancial support from AstraZeneca, Pfizer, and GlaxoSmithKline, nonfinancial support from Novartis, and personal fees from Covis outside the submitted work. D. Malka reports grants from Bristol Myers Squibb during the conduct of the study, as well as personal fees and nonfinancial support from Amgen, Bayer, Bristol Myers Squibb, Merck Serono, Merck, Pierre Fabre, Roche, Sanofi, and Servier, and personal fees from Incyte, Viatris, and AstraZeneca outside the submitted work. A. Hollebecque reports grants from Bristol Myers Squibb during the conduct of the study, as well as grants and personal fees from Incyte, grants from AstraZeneca, Bristol Myers Squibb, Boehringer Ingelheim, Janssen Cilag, Novartis, Pfizer, Roche, and Sanofi, and personal fees from Amgen, Bayer, Debiopharm, QED Therapeutics, Eisai, and Lilly outside the submitted work. D. Pouessel reports grants and personal fees from AstraZeneca and Bristol Myers Squibb, grants from Incyte, Roche, Janssen Oncology, and Seagen, grants, personal fees, and nonfinancial support from MSD, personal fees and nonfinancial support from Pfizer, and personal fees from Merck and Ipsen outside the submitted work. C. Gomez-Roca reports other support from Institut National du Cancer during the conduct of the study; grants and personal fees from Roche/Genentech and Bristol Myers Squibb, personal fees from Eisai and Foundation Medicine, and nonfinancial support from MSD, Pierre Fabre, and Macomics outside the submitted work; and is an ESMO Officer. M. Guillet reports personal fees from Bayer and Servier, and nonfinancial support from Sanofi, Roche, and Amgen outside the submitted work. J.-J. Grob reports personal fees from Bristol Myers Squibb, Merck, Novartis, Roche, Philogen, Sanofi, Pierre Fabre, and Amgen outside the submitted work. C. de la Fouchardiere reports personal fees from Ipsen and Incyte, grants and nonfinancial support from Pierre Fabre Oncologie, and personal fees and nonfinancial support from MSD, Servier, Amgen, Bayer, Eisai, Lilly, Daiichi Sankyo, Merck, and Roche outside the submitted work. F. Rolland reports personal fees from Bristol Myers Squibb and MSD during the conduct of the study. S. Hiret reports personal fees and other support from Takeda, Bristol Myers Squibb, Roche, and AstraZeneca outside the submitted work. E. Saada-Bouzid reports personal fees from Bristol Myers Squibb, MSD, and Merck Serono outside the submitted work. O. Bouche reports personal fees from Merck KGaA, Roche, Bayer, AstraZeneca, Grunenthal, MSD, Amgen, Servier, and Pierre Fabre outside the submitted work. T. Andre reports personal fees from Bristol Myers Squibb, AstraZeneca, Amgen, Astellas Pharma, Gritstone Oncology, GlaxoSmithKline, HallioDx, Kaleido Biosciences, MSD, Pierre Fabre, Roche, Seagen, Servier, and Transgene outside the submitted work. S. Oudard reports grants and personal fees from Bristol Myers Squibb, Pfizer, Novartis, Sanofi, and Bayer, and personal fees from Merck, Janssen, and Astellas outside the submitted work. C. Tournigand reports nonfinancial support and other support from Bristol Myers Squibb during the conduct of the study, as well as nonfinancial support and other support from MSD, Sanofi, and Bayer outside the submitted work. J.-C. Soria reports other support from Gritstone Bio, Relay Therapeutics, Amgen, and AstraZeneca outside the submitted work. S. Champiat reports personal fees from Amgen, AstraZeneca, Bristol Myers Squibb, Eisai, Janssen, MSD, Novartis, and Roche, other support (as principal investigator of clinical trials) from AbbVie, Amgen, Cytovation, Eisai, Imcheck Therapeutics, Molecular Partners Ag, MSD, Ose Pharma, Pierre Fabre, Sanofi Aventis, Sotio A.S, and Transgene, other support (as an advisory board member) from Alderaan Biotechnology, Amgen, AstraZeneca, Avacta, Ellipses Pharma, Oncovita, Seagen, and UltraHuman, and other support (travel and congress) from AstraZeneca, MSD, Ose Pharma, Roche, and Sotio outside the submitted work; as part of the Drug Development Department (DITEP) for Gustave Roussy is/has been the principal/subinvestigator of clinical trials for AbbVie, Adaptimmune, Adlai Nortye USA Inc., Aduro Biotech, Agios Pharmaceuticals, Amgen, Argen-X Bvba, Arno Therapeutics, Astex Pharmaceuticals, AstraZeneca Ab, Aveo, Basilea Pharmaceutica International Ltd., Bayer Healthcare Ag, Bbb Technologies Bv, Beigene, BicycleTx Ltd., Bioalliance Pharma, Blueprint Medicines, Boehringer Ingelheim, Boston Pharmaceuticals, Bristol Myers Squibb, Ca, Celgene Corporation, Chugai Pharmaceutical, Cullinan-Apollo, CureVac, Daiichi Sankyo, Debiopharm, Eisai, Eisai Limited, Eli Lilly, Exelixis, Faron Pharmaceuticals Ltd., Forma Tharapeutics, Gamamabs, Genentech, GlaxoSmithKline, H3 Biomedicine, Hoffmann La Roche Ag, Imcheck Therapeutics, Innate Pharma, Institut De Recherche Pierre Fabre, Iris Servier, Iteos Belgium SA, Janssen Cilag, Janssen Research Foundation, Kura Oncology, Kyowa Kirin Pharm. Dev, Lilly France, Loxo Oncology, Lytix Biopharma As, MedImmune, Menarini Ricerche, Merck Sharp & Dohme-Chibret, Merus, Molecular Partners Ag, Nanobiotix, Nektar Therapeutics, Novartis Pharma, Octimet Oncology Nv, Oncoethix, Oncopeptides, Onyx Therapeutics, Orion Pharma, Oryzon Genomics, Ose Pharma, Pfizer, Pharma Mar, Pierre Fabre Médicament, Plexxikon, Roche, Sanofi Aventis, Seattle Genetics, Sotio A.S, Syros Pharmaceuticals, Taiho Pharma, Tesaro, Turning Point Therapeutics, and Xencor; research grants from AstraZeneca, Bristol Myers Squibb, Boehringer Ingelheim, GlaxoSmithKline, INCa, Janssen Cilag, Merck, Novartis, Pfizer, Roche, and Sanofi; and nonfinancial support (drug supplied) from AstraZeneca, Bayer, Bristol Myers Squibb, Boehringer Ingelheim, GlaxoSmithKline, MedImmune, Merck, NH TherAGuiX, Pfizer, and Roche. M.F. Lamendola-Essel reports personal fees from Rockefeller University outside the submitted work. S.B. Maron reports personal fees from Natera, Bicara, Novartis, Daiichi Sankyo, and Basilea, and nonfinancial support from Bayer, Guardant Health, and Genentech outside the submitted work, as well as prior Calithera stock ownership. G. Argiles reports personal fees from Gadeta BV and other support from Amgen outside the submitted work. N.H. Segal reports personal fees from AstraZeneca, GlaxoSmithKline, Novartis, ABL Bio, Revitope, Amgen, Boehringer Ingelheim, and PsiOxus, grants and personal fees from Roche/Genentech, Immunocore, and PureTech Ventures, and grants from Pfizer, Merck, Bristol Myers Squibb, and Regeneron outside the submitted work. A. Cercek reports personal fees from Bayer, Merck, GlaxoSmithKline, Seagen, and Janssen, and grants from Seagen, GlaxoSmithKline, and RGenix outside the submitted work. L.A. Diaz Jr is a member of the board of directors for Jounce Therapeutics; is a compensated consultant for PetDx, Innovatus CP, Se'er, Delfi, Kinnate, and Neophore; is an uncompensated consultant for Merck but has received clinical trial support from Merck; is an inventor of multiple licensed patents related to technology for circulating tumor DNA analyses and mismatch repair deficiency for diagnosis and therapy from Johns Hopkins University (some of these licenses and relationships are associated with equity or royalty payments directly to Johns Hopkins and L.A. Diaz Jr); and holds equity in Jounce Therapeutics, PetDx, Se'er, Delfi, Kinnate, and Neophore. L.A. Diaz Jr divested his equity in Personal Genome Diagnostics to LabCorp in February 2022 and divested his equity in Thrive Earlier Detection to Exact Biosciences in January 2021. His spouse holds equity in Amgen. The terms of all these arrangements are being managed by Johns Hopkins and Memorial Sloan Kettering in accordance with their conflict of interest policies. P. Saintigny reports grants from Illumina, HTG Molecular Diagnostics, AstraZeneca, Roche, Bristol Myers Squibb Foundation, Novartis, and Inivata, and personal fees from Archer and Roche Molecular Diagnostics during the conduct of the study, as well as grants from Omicure outside the submitted work. A. Marabelle reports grants from Ligue contre le Cancer and Bristol Myers Squibb, and nonfinancial support from Institut National du Cancer during the conduct of the study; grants, personal fees, nonfinancial support, and other support from Bristol Myers Squibb, MSD, and Roche/Genentech, personal fees and other support from Pfizer/Merck Serono, other support from AstraZeneca, and grants, personal fees, and other support from Sanofi outside the submitted work; and was the lead investigator of the KN-158 trial that contributed to the approval of pembrolizumab for MSI- and TMB-high tumors. No disclosures were reported by the other authors.One of the Editors-in-Chief is an author on this article. In keeping with the AACR's editorial policy, the peer review of this submission was managed by a member of Cancer Discovery's Board of Scientific Editors, who rendered the final decision concerning acceptability.
Authors’ Contributions
B. Rousseau: Conceptualization, formal analysis, supervision, investigation, methodology, writing–original draft, writing–review and editing. I. Bieche: Conceptualization, formal analysis, supervision, investigation, methodology, writing–original draft, writing–review and editing. E. Pasmant: Conceptualization, formal analysis, supervision, investigation, methodology, writing–original draft, writing–review and editing. N. Hamzaoui: Conceptualization, data curation, formal analysis, supervision, investigation, writing–review and editing. N. Leulliot: Formal analysis, investigation, visualization, writing–review and editing. L. Michon: Data curation, software, investigation, methodology. A. de Reynies: Data curation, software, investigation, visualization, methodology. V. Attignon: Resources, data curation. M.B. Foote: Formal analysis, investigation, writing–review and editing. J. Masliah-Planchon: Formal analysis, investigation, writing–review and editing. M. Svrcek: Formal analysis, investigation, visualization, writing–review and editing. R. Cohen: Investigation, writing–review and editing. V. Simmet: Investigation, writing–review and editing. P. Augereau: Investigation, writing–review and editing. D. Malka: Investigation, writing–review and editing. A. Hollebecque: Investigation, writing–review and editing. D. Pouessel: Investigation, writing–review and editing. C. Gomez-Roca: Investigation, writing–review and editing. R. Guimbaud: Investigation. A. Bruyas: Investigation, writing–review and editing. M. Guillet: Investigation, writing–review and editing. J.-J. Grob: Investigation. M. Duluc: Investigation. S. Cousin: Investigation. C. de la Fouchardiere: Investigation. A. Flechon: Investigation. F. Rolland: Investigation. S. Hiret: Investigation. E. Saada-Bouzid: Investigation. O. Bouche: Investigation. T. Andre: Investigation, writing–review and editing. D. Pannier: Investigation. F. El Hajbi: Investigation. S. Oudard: Investigation. C. Tournigand: Investigation. J.-C. Soria: Investigation. S. Champiat: Conceptualization. D.G. Gerber: Data curation, investigation. D. Stephens: Data curation, investigation. M.F. Lamendola-Essel: Resources, project administration, writing–review and editing. S.B. Maron: Data curation, investigation, writing–review and editing. B.H. Diplas: Data curation, investigation. G. Argiles: Investigation, writing–review and editing. A.R. Krishnan: Data curation, investigation. S. Tabone-Eglinger: Resources, data curation. A. Ferrari: Resources, data curation, software. N.H. Segal: Investigation, writing–review and editing. A. Cercek: Investigation. N. Hoog-Labouret: Supervision. F. Legrand: Supervision. C. Simon: Resources, supervision, project administration. A. Lamrani-Ghaouti: Data curation, supervision, project administration. L.A. Diaz Jr: Supervision, investigation, writing–review and editing. P. Saintigny: Conceptualization, resources, formal analysis, supervision, investigation, methodology, writing–review and editing. S. Chevret: Conceptualization, formal analysis, supervision, methodology, writing–original draft, writing–review and editing. A. Marabelle: Conceptualization, supervision, investigation, methodology, writing–original draft, writing–review and editing.
Acknowledgments
The Acsé Nivolumab trial is sponsored by the French network of comprehensive cancer centers (Unicancer) and the Institut National du Cancer (INCa). The funding of the trial has been supported by Bristol Myers Squibb, the Ligue contre le Cancer (www.ligue-cancer.net), and l'INCa. The supply of nivolumab was kindly provided by Bristol Myers Squibb. Blood PBMCs and tumor FFPE blocks were prospectively collected, stored, and prepared for analysis by the Léon Bérard Cancer Center (Centre de Ressources Biologiques du Centre Léon Bérard, #BB-0033-00050, Lyon, France). The sequencing of tumors (WES and RNA-seq) was part of the ancillary analysis program of Acsé Nivolumab (Acsé Cible) funded by a research grant from the Ligue contre le Cancer. B. Rousseau thanks the Nuovo-Soldati and Swim Across America foundations for salary support. This work was funded in part by the Marie-Josée and Henry R. Kravis Center for Molecular Oncology and NCI Cancer Center Core Grant P30-CA008748. We gratefully acknowledge the members of the Molecular Diagnostics Service in the Department of Pathology of MSKCC. Also, we thank Mr. Daniel Couch, the genomic platform personnel who screened for patients harboring POLE mutations, all the study coordinators, the nurses, and the patients and their families who have actively participated in the POLE cohort of the Acsé Nivolumab trial. This work was supported by Bristol Myers Squibb, Unicancer, the Ligue contre le Cancer, and INCa. M.B. Foote is funded by an NIH T32-CA009512 and an ASCO Young Investigator Award. This work was also supported by funding from the Integrated Cancer Research Site LYriCAN (INCa-DGOS-Inserm_12563) and PIA Institut Convergence Francois Rabelais PLAsCAN (PLASCAN, 17-CONV-0002).