- Open Access
Complete response to neoadjuvant chemoradiotherapy in rectal cancer is associated with RAS/AKT mutations and high tumour mutational burden
Radiation Oncology volume 16, Article number: 129 (2021)
Pathological complete response (pathCR) in rectal cancer is beneficial, as up to 75% of patients do not experience regrowth of the primary tumour, but it is poorly understood. We hypothesised that the changes seen in the pre-treatment biopsies of pathCR but not seen in residual tumour after chemoradiotherapy were the determinants of responsiveness.
Two groups of patients with either complete response (pathCR group, N = 24) or no response (poor response group, N = 24) were retrieved. Pre-treatment biopsies of cancers from these patients underwent high read depth amplicon sequencing for a targeted panel, exome sequencing, methylation profiling and immunohistochemistry for DNA repair pathway proteins.
Twenty four patients who underwent pathCR and twenty-four who underwent poor response underwent molecular characterisation. Patients in the pathCR group had significantly higher tumour mutational burden and neoantigen load, frequent copy number alterations but fewer structural variants and enrichment for driver mutations in the PI3K/AKT/mTOR signalling pathway. There were no significant differences in tumour heterogeneity as measured by MATH score. Methylation analysis demonstrated enrichment for hypomethyation in the PI3K/AKT/mTOR signalling pathway.
The phenomenon of pathCR in rectal cancer may be related to immunovisibility caused by a high tumour mutational burden phenotype. Potential therapy resistance mechanisms involve the PI3K/AKT/mTOR signalling pathway, but tumour heterogeneity does not seem to play a role in resistance.
Rectal cancer is a common malignancy , with approximately 11,000 cases per year in the UK . Treatment typically consists of excisional surgery  with neoadjuvant therapy if the cancer is locally advanced, consisting of either short course radiotherapy  (25 Gy in 5 fractions over 1 week) with surgery the following week, or long course neo-adjuvant chemoradiotherapy (nCRT)  (45–50 Gy in 25 fractions over 5 weeks, with synchronous iv 5-FU or oral capecitabine) with surgery 6–10 weeks later. The former regimen usually demonstrates little, if any tumour regression, but if it does occur, is associated with more favourable prognosis, and the latter regimen can lead to significant tumour shrinkage and downstaging, with pathological complete response (pathCR, defined as complete regression of tumour in the resection specimen) observed in approximately 10–15% of patients . Multiple investigators have shown higher rates of response [7, 8] with higher radiotherapy doses, with maximum pCR rates of 25–30%. The current standard of care within the United Kingdom for locally advanced rectal cancer is long course chemoradiotherapy, however the recent development of the concept of “total” neo-adjuvant therapy (TNT), whereby consolidation chemotherapy is given after chemoradiotherapy, increased pCR rates from 15 to 25% in the consolidation group, acting as an additional therapeutic option in these patients which may lead to its introduction as standard of care . Short course radiotherapy with a delay to surgery , or with neoadjuvant chemotherapy in the RAPIDO trial , has also been shown to be of benefit and demonstrates the critical nature of the sequencing and timing of treatment in order to maximise response.
Clinical complete response (CCR) is defined as the absence of tumour on imaging and/or clinical examination, but does not definitively exclude residual tumour, which requires a resection specimen to confirm. CCR is well correlated with pathCR and, this could allow routine use of a watch and wait strategy. In addition to its role as an indicator of potential cure  (defined as > 5 years recurrence free), path CR could be used to delay excisional surgery  or allowing organ preservation surgery , such as TEMS or TAMIS .
The molecular drivers of pathCR are unclear [16, 17], but they are thought to relate to factors that promote radiotherapy- and chemotherapy-related tumour cell killing. For example, rectal cancers exist within a low oxygen tension environment  leading to intrinsic resistance. As radiotherapy induces the formation of oxygen derived free radicals, which cause tumour cell death by direct DNA damage, a low oxygen environment leads to fewer available oxygen free radicals, leading to lower cell death. Another hypothesised mechanism underlying differences pathCR is variation in DNA repair. Ionising radiation  induces a DNA double strand break (DSB), which is then repaired either by homologous recombination (HR), where a sister chromatid is used to repair the defect, or non-homologous end joining (NHEJ), where a complex repair mechanism including BRCA, ATM, ERCC5 and others contribute to a DNA repair complex that re-joins the damaged segments of DNA. There is evidence from previous studies that aberrant functioning of the NHEJ pathway is associated with a longer survival, presumably as a consequence of better response in radiotherapy , however the data concerning the role of NHEJ and the response to radiotherapy is unclear, with studies giving conflicting reports of its relationship with radiation response . Lal et al.  also showed that the immune context of a tumour depended on genomic factors, with KRAS mutation, CMS2 or CMS3 classification all being independently associated with a reduced immune infiltration. This has implications for responsiveness to therapy, as a large proportion (40%) of colorectal cancer has KRAS mutation, and immunovisibility of the tumour is essential for a good response to neoadjuvant therapy.
Tumour heterogeneity undoubtedly also plays a role in determining pathCR [16, 23], as does immunogenicity caused by the formation and expression of clonal neoantigens. Rectal cancers have little hypermutation and therefore rates of immunovisibility are low. Whilst Akiyoshi et al.  have recently shown that levels of clonal neoantigens are higher in patients undergoing a good response to neoadjuvant chemoradiotherapy, suggesting that hypermutation is important, a previous integrated molecular analysis of rectal cancer  found no key genomic features that correlated with resistance. This may be because of the significant heterogeneity of the datasets used in the analysis.
Owing to the uncertainty surrounding the precise mechanisms of sensitivity of rectal cancer to chemoradiotherapy, we aimed to study the phenomenon of pathCR in rectal cancer. We hypothesised that the genomic changes responsible for the phenomenon of pathCR would be seen uniquely in the pre-treatment biopsies of those undergoing a complete response, and not seen at all in the residual post-treatment specimens of patients who had undergone neo-adjuvant chemoradiotherapy and therefore would represent two divergent opposites to identify the phenomenon.
A prospective database of all patients undergoing neoadjuvant chemoradiotherapy for rectal cancer was used to identify patients. Ethical approval was obtained from the NW Research Ethics committee (ref 15/NW/0079). Patients who underwent long course chemoradiotherapy and achieved pathCR were identified, as determined by complete regression of the tumour, with absolutely no tumour cells remaining (Mandard grade 1/ TRG 1), on examination of the specimen by a Consultant Histopathologist, as opposed to minimal residual disease where several cells were allowable. All patients underwent long course chemoradiotherapy with either oral capecitabine (825 mg/m2) or infusional 5-flourouracil during a radiotherapy course of 45 Gy in 25 fractions over 35 days. Pre-treatment stage and post-treatment response (at 6 weeks after finishing treatment) was assessed using magnetic resonance imaging. Resection of the primary tumour then occurred as soon as possible after the 6 week MRI scan. The histopathology archives were then searched to find the pre-treatment endoscopic biopsies of these patients for downstream analysis. This was defined as the “pathCR cohort”. A randomly selected second cohort of rectal cancer patients was identified, who had no response to treatment or progressed whilst on treatment, as defined both by post-treatment MRI and Mandard grade 4 or 5. This was defined as the “non-responder cohort”. For the validation cohort, patients from work stream 3 of the Stratification in Colorectal cancer (S-CORT) were utilised, which represent a cohort of patients undergoing long course chemoradiotherapy for whom pre-treatment biopsy material was available. For germline DNA, representative normal tissue from the proximal resection margin of the discovery cohort was obtained for genomic analysis.
Formalin fixed blocks were retrieved for these patients and cut into 4uM sections on frosted glass slides for needle macrodissection and immunohistochemistry. Pre-treatment specimens were all endoscopically obtained biopsies, whereas post-treatment specimens consisted of tumour blocks selected by a histopathologist as representative of the tumour. This consisted of the block whereby the H&E section showed maximum tumour content. A representative H&E section was used to target tumour cells for macrodissection to maximise tumour content for DNA extraction using a modified protocol of the Qiagen DNEasy kit (Qiagen). Eluted DNA was quantified using Nanodrop spectrophotometry (for contaminants) and Qubit fluorimetry (for concentration).
This was performed on 4uM slides as previously described using a Leica BondMax autostainer. IHC was performed against yH2AX (Abcam ref ab26350), ATM (ab36810), Ku70/80 (ab53126), MLH1 (ab92312), MSH2 (ab52266) and MSH6 (ab14204) and slides were then scanned on a Leica slide imaging platform. IHC was scored by the QuPath system . Briefly, a single section was zoomed to 10 × view of a representative area of epithelium/tumour, calibration was performed according to the QuPath manual, either nuclei or membranous counting was set depending on the antibody and auto counting of DAB stained positive/negative cells was carried out. Staining was reported as a percentages of positive/negative cells.
Extracted genomic DNA was used for targeted amplicon resequencing using the Fluidigm 48 × 48 Access array. PCR amplicons covering APC, KRAS, BRAF, NRAS, FBXW7 and SOX7 (primer sequences available on request) were designed using Primer3 . Primer specificity was checked using PrimerBLAST and UCSC in-silico PCR. Briefly, 20 ng of FFPE DNA was injected into the Access Array system with PCR primers and thermal cycled according to manufacturer’s specifications. Amplicons were then ligated to Illumina sequencing indexes & adapters and pooled and sequenced on an Illumina MiSeq to an average read depth of > 1000 × using a 100 bp PE sequencing strategy. For FFPE exome sequencing, a custom modification of the Illumina TruSeq Exome hybridisation kit was used. At least 300 ng of FFPE-derived DNA was prepared with the following modifications: Firstly, no size selection was performed after end repair and DNA fragments were amplified with 12 cycles of PCR. Enrichment was performed using a bead ratio of 0.8, then samples were combined into pools of 3 plex for coding exome (TruSeq exome, 45mb in size) probe hybridisation and subsequent clean up. 10 cycles of amplification were performed to enrich the final libraries which were then pooled into a final 12 plex library. Sequencing was performed on an Illumina NextSeq. For the S-CORT WS3 samples, a custom panel consisting of 61 oncogenes was sequenced to 500 × using Agilent SureSelect bait capture at the Wellcome Trust Sanger Institute. Methylation was interrogated using the Illumina HumanMethylation 450 array system. Between 100 and 250 ng of FFPE DNA was processed using the Illumina FFPE restore kit, and then hybridised to the HumanMethylation 450 array following the manufacturer’s instructions. Slides were washed and scanned on an Illumina iScan scanner.
Imaging mass cytometry
Slide staining and CyTOF were performed on FFPE sections using methods as previously described . CyTOF IMC data were analysed using an image-processing pipeline as described (https://github.com/BodenmillerGroup/imctools). Ilastik generated Probability probability maps and raw multi-channel files from each region of interest were analysed using CellProfiler and the Cytomapper R package (v3.12) [27, 28].
For the amplicon resequencing, FASTQ files were trimmed (Trimgalore) and aligned to the GRCh38 reference genome using a standard pipeline using the BWA (v0.7.17-r1188) aligner . Mutation calling was performed using FreeBayes (v1.0.0) . For the samples from the SCORT consortium, alignment using BWA was carried out to the GRCh38 reference genome and mutation calling performed with Caveman (v1.14.0) and Pindel (v1.0) [31, 32]. For the exome sequencing analysis, FASTQ files were trimmed and aligned to the GRCh38 reference genome using an exome sequencing pipeline Isaac v4 aligner , Manta SV caller (v1.6.0)  and Canvas CNV (Canvas 126.96.36.1993 + master)  caller. Enrichment was determined by comparison to the Illumina TruSeq exome v1.2 BED file. Mutation calling was performed using Strelka 2  and MuTect2 in tumour-normal subtraction mode. For both amplicon sequencing and exome analysis, significantly mutated genes were identified using Intogen (v1.0) , MutSigCV (v2.1)  and dNdScv (v.0.1.0) . Mutational signatures were estimated using the MutationalPatterns (v3.0.1)  R package. Tumour mutational burden was estimating by annotation of the combined, overlapping variant calls from Mutect/Strelka2 to filter to non-synonymous variants (nsV). NsV were then divided by the TruSeq Exome panel size (45.3 Mb) to give a figure in mutations/mb.
For structural variant (SV) calls, Manta was used in tumour/normal mode on exome sequencing data and lists of structural variants outputted to VCF file. For copy number calls (CNV), Canvas was used in tumour normal mode and calls outputted as VCF files.
Clonality was determined by running superFreq (v1.0) , a cancer specific tumour exome caller and river plots were producing from this package. For neoantigen calls, tumour-normal subtracted VCF files produced with Strelka, GT fields were added via conversion of the SGT field (using a custom script), annotated with the variant effect predictor (VEP), filtered for indels and then analysed using PVacTools v2.0.0 against MHC Class I binding predictions . Neoantigens with a “best” median IC of < 50 nmol L−1 were counted as binding neoantigens. HLA typing for each patient was determined with HLA-LA* (v1.0.1) on germline exome sequencing data .
For methylation, IDAT files were imported into R/Bioconductor and analysed on the ChAMP (v2.20.1) pipeline  using standard settings. iDAT files were imported, standardised to beta methylation values, then underwent SWAN normalisation. Normalised beta-values underwent differential methylation analysis using limma, DMR were called using dmrLasso. Pathway analyses were conducted in GProfiler2 .
For statistical analysis, Stata 15.1 was used, and all distributions plotted to check for a Gaussian distribution. Unpaired t-testing was set with a significance threshold of 0.05, and any missing values were imputed using the Stata impute commands.
In total, there were 48 patients with samples available in the study (Fig. 1). In the pathCR group, 24 pre-treatment biopsies were available. In the post-treatment group there were 24 post-treatment samples that had a poor response to chemoradiotherapy as defined by histological tumour regression grade. The S-CORT WS3 cohort had 231 pre-treatment biopsies available for analysis with variable levels of response. Matched germline DNA was available for all samples apart from the S-CORT WS3 cohort.
All available samples successfully underwent amplicon sequencing and exome sequencing. For the amplicon sequencing average read depth was 1020 × (range 357–5210 ×). For the whole exome sequencing the average read depth was 99 × for tumour samples and 55 × for control normal. For methylation arrays all samples hybridised successfully.
In the amplicon sequencing group, using treatment response as a classifier to select top significantly mutated genes correlated with response as defined by MutSigCV,were PIK3CA (p = 2.4 × 10−9), FBXW7 (p = 8.68 × 10−4) and PTEN (p = 1.02 × 10−3). Neither dNdScV nor Intogen can classify mutations by co-variates therefore a Fisher exact test was used to compare the two groups (pathCR vs. non-responder) for the output of each caller. Using dNdScV, the top significantly mutated genes were AKT1 (pmis = 2.01 × 10−3), FBXW7 (pmis = 5.90 × 10−3), FAM123B (pmis = 2.48 × 10−2) and POLE (pmis = 1.32 × 10−2). Gene centric analysis with Intogen demonstrated recurrent mutations in PIK3CA, PTEN, FBXW7 and POLE. Pathway analysis of the genes present in the pathCR group but not the non-responder group demonstrated 14 pathways (Table 1) being significantly over-represented in the dataset including hsa05210 (Colorectal Cancer, p = 1.88 × 10−16, q = 1.62 × 10−14), hsa05222 (Small cell lung cancer, p = 8.43 × 10−12, q = 3.63 × 10−10), and hsa04150 (mTOR pathway, p = 4.99 × 10−10, q = 1.43 × 10−8), all of which contain genes involved in the mTOR/AKT pathway.
For the exome sequencing data, analysis by MutSigCV using pathCR as a covariate revealed 1,412 genes that were significantly mutated (Additional file 1). The top genes were HIVEP3, HS6ST3, KIAA1671, LRRC4C and ROBO2. This was due to the preponderance of hypermutant samples within the pathway which lead to select of non-driver genes, and filtering would have removed samples. Pathway analysis by GProfiler demonstrated no enrichment of any KEGG pathways. Analysis by Intogen using revealed demonstrated 1620 genes significantly mutated (Additional file 1 as determined either by OncoDriveClust or OncoDriveFM. The top five genes were CDC27, CTBP2, IGSF3, PABPC3 and ZNF432. Pathway analysis with Intogen demonstrated over-representation of focal adhesion (hsa04510, p = 2.25 × 10−135, q = 5.82 × 10−133), which contains within it the MAPK/PIK3K and Wnt signalling pathways; as the top rated pathway (Additional file 1). Analysis by dNdScV demonstrated 209 genes significantly mutated (p < 0.05), the top five of which were ZNF717, MUC3A, APC, OR4C5 and KRAS.
Mutational signature analysis (version 3 single base substitution) was performed on the exome sequencing samples from the pathCR group. The top ranked signatures were signatures three, five and thirty. Signature three is postulated to be due to defective base repair due to faulty homologous recombination, signature five is due to the effects of transcription coupled nucleotide repair and signature 30 is due to a defect in base-excision repair due to inactivating mutations in NTHL1 (Additional file 1).
Tumour mutation burden (TMB) was significantly higher in pre-treatment samples from pathCR patients than non-responders (exome data, medianpathCR = 38.28 muts/mb, range 15.93–86.95 v medianNR 7.27 muts/mb, range 3.08–47.3, p = 0.02 Wilcoxon, Fig. 2).
In order to validate the finding that tumour mutational burden seemed to correlate with response, we analysed data from the SCORT WS3 “Grampian” cohort consisting of pre-treatment biopsies (n = 231) from patients undergoing long course chemoradiotherapy. Total numbers of coding mutations were significantly (p = 0.036) greater in the “responder” (n = 35, mean 9.2 mutations) vs. “non-responder” (n = 196, mean 7.8 mutations) groups.
Structural variants and copy number variation
Copy number estimation was performed on all tumour: normal exome pairs successfully (Fig. 3). Three biopsy samples showed an extremely complex pattern of copy number gain and loss (with modal CN of 67, 52 and 66 respectively), which corresponded with complete pathological response to chemoradiotherapy. The median copy number of the pathCR group was 50 (IQR 45–66) and of the non-responder group was 46 (IQR 44–47, Wilcoxon Ranked sums p = 0.01). No recurrent copy number alterations were observed. In the poor response group, the copy number seen in the pre-treatment biopsy (median 45, IQR 44–49) did not significantly vary with the post treatment specimen (median 46, IQR 44–47, Wilcoxon Ranked sums p = 0.94).
Structural variant (SV) calling was performed on all tumour: normal exome pairs successfully. The median number of structural variants in the pre-treatment biopsies of the pathCR group was 4 (IQR 0–10) vs. 52 in non-responders group (IQR 18–80, p < 0.001, Fig. 2). There was no significant difference in the numbers of SV between the pre-treatment biopsies of the non-responders and the post-treatment specimens.
Clonal evolution analysis
In the seven cases from whole trios of normal, pre-treatment biopsy and post-treatment specimens were available (only in the poor response group), clonal evolution analysis with superFreq was carried out. This revealed that subsequent to radiation therapy there was an increase in clonal diversity (Fig. 4) in all patients, suggesting that radiation therapy drove the generation of an increase in clonal diversity. For both the pathCR group and the non-responder group, we calculated MATH score, a measure of tumour heterogeneity for all samples in the pre-treatment biopsies of all samples, finding that there was a median MATH score of 4.7 in non-responders and 4.2 in pathCR group (Mann Whitney p = 0.5036), suggesting that tumour heterogeneity did not play a role in treatment resistance (Fig. 5).
In order to understand if there were epigenetic determinants of response, we compared the methylome of the pre-treatment biopsies of the pathCR group to the tumour of the no response group. This demonstrated 1853 differentially methylated positions and Go:Profiler analysis of these differentially methylated genes demonstrated that the WikiPathways WP306 (Focal adhesion p = 1.16 × 10−3), WP185 (Integrin mediated focal adhesion p = 2.42 × 10−3) and WP3932 (PI3K-Akt-MTOR pathway p = 3.31 × 10−2) were enriched in the pathCR group for demethylation, suggesting potential activation of this pathway. Pathway WP306 contains the PI3K/Akt/MTOR signalling genes as well as WP3932 and represents a general enrichment of this pathway.
Analysis of differentially methylated regions (contiguous blocks of differentially methylated CpGs > 10 in number, DMR) 347 differentially methylated regions. HOMER annotation of top 100 DMRs and then analysis by g:Profiler revealed only enrichment for methylation in the Human Phenotype pathway HP0000104 Renal Agenesis (p = 2.82 × 10−2), which contains Wnt signalling and DNA repair genes (ATRX) and the GO:MF term Class II MHC binding (p = 1.77 × 10−2).
Neoantigen prediction using pVacTools revealed a median of 78 neoantigens per sample (range 9–383) of which the pathCR group had a median 91 neoantigens (IQR 78–154) and the non-responder group a median 74 neoantigens (IQR 28–89, Wilcoxon rank sum p = 0.034).
DSB pathway and mismatch repair deficiency
Semi-quantitative analysis using QuPath was carried out on white light DAB stained images of the expression of yH2AX, Ku70/80, ATM, MLH1, MSH2 and MSH6 respectively. For the pathCR group the pre-treatment biopsies only were studied and for the non-responder group the post-treatment specimens.
For DNA repair proteins, pathCR cases showed significant differences compared to the non-responder group in the DNA repair associated proteins γH2AX (89% cells expressing vs. 78% cells expressing, p = 0.007), ATM (80% vs. 69%, p = 0.04) and Ku70/80 (88% vs. 29%, p < 0.001), suggesting that a deficiency of normal DNA double strand break proteins was associated with resistance to treatment.
In the mismatch repair associated proteins, pathCR group showed significant differences compared to the non-responder group in MLH1 (70% cells expressing vs. 84% cells expressing, p = 0.001) and MSH6 (76% vs. 88%, p < 0.001). However MSH2 expression (33% vs. 41%, p = 0.10) was not significant. Especially for the mismatch repair proteins, despite there being statistically significant differences, these were only a few percentage points different and little complete loss of expression was seen, suggesting that the role of these proteins is uncertain in this context.
Imaging mass cytometry analysis
A 40 marker image mass cytometry analysis was performed, and cell counting using CellProfiler determined that with 44 regions of interest across the responders vs. non-responders there were significantly more (proportions as part of total immune cells compartment 7.98% vs. 3.06%, p = 0.0021) CD8 + T-cells infiltrating in tumours with a complete response to radiotherapy than those with a poor response.
In this study, we have demonstrated a number of interesting findings comparing pre-treatment biopsies and post treatment specimens in patients undergoing (nCRT) for rectal cancer.
The most important finding is that in patients achieving a complete or near complete response to nCRT is that the pre-treatment tumour has a high tumour mutational burden . This is because hypermutation within the tumour due to an intrinsic DNA repair defect and leads to increased presentation of neoantigens because of indel and frameshift mutations . This increases immunovisibility and we hypothesise that the increased immunovisibility, coupled with activation of the immune system by irradiation, causes cGAS/STING activation which has been shown to lead to a Type I interferon response  leading to migration of immune cells and enhanced regression. Also, it suggests that these patients may benefit from neo-adjuvant anti-PD1/PD-L1 or anti-CTLA4 immunotherapy as high TMB (defined as > 10 muts/mb) has been shown to be correlated with responsiveness to these agents . Another intriguing possibility is that genotoxic therapy, coupled with radiation therapy could be delivered as part of neoadjuvant treatment in the tumour, possibly increasing neoantigen burden and making more patients suitable for immunotherapy . Our sample size of pathological complete responders is relatively small, however we deliberately chose tumour regression where absolutely no cells remained, which is a very rare phenomenon, compared to the phenomenon of “minimal residual disease”, but we believed that this would give a stronger biological signal. We found that this was replicated in a larger cohort (SCORT WS3) although this cohort only underwent limited amplicon sequencing which can serve at best as a proxy for tumour mutational burden.
We have also demonstrated that there is enrichment for mutations in the mTOR/AKT signalling pathways, specifically PIK3CA but also as a general trend towards mutations and epigenetic changes in this pathway. PIK3CA is a member of the mTOR signalling pathway and makes up the alpha subunit of the PI3K protein. mTOR signalling has previously been highlighted as being of possible relevance in radiosensitivity  as a cellular marker of stress in the tumours of patients undergoing chemoradiotherapy. PIK3CA signals through the mTORC1/mTORC2 and exerts its downstream effect on AKT . Targeted agents for PIK3CA (apitosilib ), mTORC1/2 (visusertib ) and AKT (MK2206 ) exist and are at various stages of clinical development. We suggest that these agents may be utilised as part of a neoadjuvant therapy strategy in order to increase response rates.
FBXW7, a gene previously implicated in cell cycle control by ubiquitination? of cyclin-E1  was also significantly enriched for mutations within it in this study. Zhang et al. demonstrated that FBXW7 also had a role in non-homologous end joining [56, 57], being a binding partner in the complex that repairs damage caused by double strand breaks. Mutations in FBXW7 may affect its ability to participate in NHEJ and therefore increase radiosensitivity. Mutational signature analysis of the exome sequencing dataset also showed multiple mutational signatures consistent with enrichment in DNA repair, specifically faulty homologous recombination (HR) and nucleotide excision repair. The finding that both PIK3CA and FBXW7 mutations are both enriched in pre-treatment biopsies and found in post-treatment specimens from cancers that do not respond to neoadjuvant therapy means that they may act as biomarkers for lack of response. In oesophageal cancer, induction of novel mutations in cancer driver genes in post treatment specimens has been observed after chemoradiotherapy treatment and our findings may reflect this.
Tumour heterogeneity is a significant problem across all tumours due to drivers that may cause a differential response to therapy because of mutational clonality. Our results show, unsurprisingly, that treatment of samples with neoadjuvant chemoradiotherapy causes an increase in clonal diversity, which may be a driver in the phenomenon of radiation resistance. However, we did not demonstrate any difference between responders and non-responders in terms of MATH score as a measure of tumour heterogeneity, unlike other papers . As we only sampled a single region of tumour, it is also possible the observed lack of difference in heterogeneity is purely due to chance.
The immunohistochemical analysis of tumour samples showed increased yH2AX and ATM expression was associated with response, as well as increased expression of the mismatch repair proteins MLH1 & MSH6. We found this puzzling, as we would have expected the opposite to be true, i.e. loss of expression was needed for response, especially in ATM and yH2AX. However, we hypothesise that these samples may have had regions of loss of normal response and these disappeared as a consequence of response to radiotherapy, or that, is as typical tumour biopsies usually have low tumour content.
Clearly, the best way to understand these defects and investigate them further would be to build a cellular model of rectal cancer in order to modulate these pathways in order to measure responsiveness . Current cell lines have a bias towards their micro environment and although provide reasonable models of single pathway alterations lack the fidelity to measure therapy response when modulated. The ideal model would be an organoid based rectal cancer therapy model  as this provides both the 3D structure (enabling cell/cell communication a more representative element of intra-tumoural hypoxia) and more accurate response to therapy, as well as the ability to evolve and resist therapy and the ability to co-culture with other cells in the microenvironment such as T-cells and fibroblasts.
Therefore, we suggest that based on these findings, a number of factors contribute to the response to neoadjuvant chemoradiotherapy: hypermutation leading to increased neoantigen presentation; enrichment in defects in the mTOR signalling pathway; hypoxia regulated by miR-21-5p and an increase in clonal diversity. Our findings agree with those of Akiyoshi et al.  in that increase in neoantigen diversity correlated with response. Kamran  however, found no clear molecular defects that predisposed to radioresistance. This could be due to the fact that their analysis was geared towards pathways of treatment resistance rather than those that lead to radiosensitivity.
Our findings suggest a number of new therapeutic avenues for increasing responsiveness to chemoradiotherapy in rectal cancer. We plan to further study these in complex 3D models such as organoids in order to understand whether they will increase response rates when appropriately targeted.
Availability data and materials
All data will be uploaded to the SRA upon acceptance of this manuscript.
Pathological complete response
AKR thymoma pathway
Mammalian target of rapamycin
Complete clinical response
Transanal endoscopic microsurgery
Transanal minimally invasive surgery
Double strand break
Non homologous end joining
Tumour regression grade
Formalin fixed, paraffin embedded
Human leucocyte antigen
Major histocompatibility complex
Copy number variant
Mutant allele tumour heterogeneity
Akiyoshi T, Tanaka N, Kiyotani K, Gotoh O, Yamamoto N, Oba K, et al. Immunogenomic profiles associated with response to neoadjuvant chemoradiotherapy in patients with rectal cancer. Br J Surg. 2019;106(10):1381–92.
Glynne-Jones R, Wyrwicz L, Tiret E, Brown G, Rodel C, Cervantes A, et al. Rectal cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2017;28(suppl_4):iv22–iv40.
Prete FP, Pezzolla A, Prete F, Testini M, Marzaioli R, Patriti A, et al. Robotic versus laparoscopic minimally invasive surgery for rectal cancer: a systematic review and meta-analysis of randomized controlled trials. Ann Surg. 2018;267(6):1034–46.
Ngan SY, Burmeister B, Fisher RJ, Solomon M, Goldstein D, Joseph D, et al. Randomized trial of short-course radiotherapy versus long-course chemoradiation comparing rates of local recurrence in patients with T3 rectal cancer: Trans-Tasman Radiation Oncology Group trial 01.04. J Clin Oncol. 2012;30(31):3827–33.
Bujko K, Wyrwicz L, Rutkowski A, Malinowska M, Pietrzak L, Krynski J, et al. Long-course oxaliplatin-based preoperative chemoradiation versus 5 × 5 Gy and consolidation chemotherapy for cT4 or fixed cT3 rectal cancer: results of a randomized phase III study. Ann Oncol. 2016;27(5):834–42.
Dossa F, Chesney TR, Acuna SA, Baxter NN. A watch-and-wait approach for locally advanced rectal cancer after a clinical complete response following neoadjuvant chemoradiation: a systematic review and meta-analysis. Lancet Gastroenterol Hepatol. 2017;2(7):501–13.
Habr-Gama A, Sabbaga J, Gama-Rodrigues J, Sao Juliao GP, Proscurshim I, Bailao Aguilar P, et al. Watch and wait approach following extended neoadjuvant chemoradiation for distal rectal cancer: are we getting closer to anal cancer management? Dis Colon Rectum. 2013;56(10):1109–17.
Sanghera P, Wong DW, McConkey CC, Geh JI, Hartley A. Chemoradiotherapy for rectal cancer: an updated analysis of factors affecting pathological response. Clin Oncol (R Coll Radiol). 2008;20(2):176–83.
Fokas E, Allgauer M, Polat B, Klautke G, Grabenbauer GG, Fietkau R, et al. Randomized phase II trial of chemoradiotherapy plus induction or consolidation chemotherapy as total neoadjuvant therapy for locally advanced rectal cancer: CAO/ARO/AIO-12. J Clin Oncol. 2019;37(34):3212–22.
Pettersson D, Holm T, Iversen H, Blomqvist L, Glimelius B, Martling A. Preoperative short-course radiotherapy with delayed surgery in primary rectal cancer. Br J Surg. 2012;99(4):577–83.
Bahadoer RR, Dijkstra EA, van Etten B, Marijnen CAM, Putter H, Kranenbarg EM, et al. Short-course radiotherapy followed by chemotherapy before total mesorectal excision (TME) versus preoperative chemoradiotherapy, TME, and optional adjuvant chemotherapy in locally advanced rectal cancer (RAPIDO): a randomised, open-label, phase 3 trial. Lancet Oncol. 2021;22(1):29–42.
Cercek A, Goodman KA, Hajj C, Weisberger E, Segal NH, Reidy-Lagunes DL, et al. Neoadjuvant chemotherapy first, followed by chemoradiation and then surgery, in the management of locally advanced rectal cancer. J Natl Compr Canc Netw. 2014;12(4):513–9.
van der Valk MJM, Hilling DE, Bastiaannet E, Meershoek-Klein Kranenbarg E, Beets GL, Figueiredo NL, et al. Long-term outcomes of clinical complete responders after neoadjuvant treatment for rectal cancer in the International Watch & Wait Database (IWWD): an international multicentre registry study. Lancet. 2018;391(10139):2537–45.
Chadi SA, Malcomson L, Ensor J, Riley RD, Vaccaro CA, Rossi GL, et al. Factors affecting local regrowth after watch and wait for patients with a clinical complete response following chemoradiotherapy in rectal cancer (InterCoRe consortium): an individual participant data meta-analysis. Lancet Gastroenterol Hepatol. 2018;3(12):825–36.
Rombouts AJM, Al-Najami I, Abbott NL, Appelt A, Baatrup G, Bach S, et al. Can we Save the rectum by watchful waiting or TransAnal microsurgery following (chemo) Radiotherapy versus Total mesorectal excision for early REctal Cancer (STAR-TREC study)? Protocol for a multicentre, randomised feasibility study. BMJ Open. 2017;7(12):e019474.
Kamran SC, Lennerz JK, Margolis CA, Liu D, Reardon B, Wankowicz SA, et al. Integrative molecular characterization of resistance to neoadjuvant chemoradiation in rectal cancer. Clin Cancer Res. 2019;25(18):5561–71.
Agostini M, Janssen KP, Kim IJ, D’Angelo E, Pizzini S, Zangrando A, et al. An integrative approach for the identification of prognostic and predictive biomarkers in rectal cancer. Oncotarget. 2015;6(32):32561–74.
Rasheed S, Harris AL, Tekkis PP, Turley H, Silver A, McDonald PJ, et al. Hypoxia-inducible factor-1alpha and -2alpha are expressed in most rectal cancers but only hypoxia-inducible factor-1alpha is associated with prognosis. Br J Cancer. 2009;100(10):1666–73.
Sak A, Stuschke M. Use of gammaH2AX and other biomarkers of double-strand breaks during radiotherapy. Semin Radiat Oncol. 2010;20(4):223–31.
Beggs AD, Domingo E, McGregor M, Presz M, Johnstone E, Midgley R, et al. Loss of expression of the double strand break repair protein ATM is associated with worse prognosis in colorectal cancer and loss of Ku70 expression is associated with CIN. Oncotarget. 2012;3(11):1348–55.
Trenner A, Sartori AA. Harnessing DNA double-strand break repair for cancer treatment. Front Oncol. 2019;9:1388.
Lal N, White BS, Goussous G, Pickles O, Mason MJ, Beggs AD, et al. KRAS mutation and consensus molecular subtypes 2 and 3 are independently associated with reduced immune infiltration and reactivity in colorectal cancer. Clin Cancer Res. 2018;24(1):224–33.
Blank A, Roberts DE, 2nd, Dawson H, Zlobec I, Lugli A. Tumor heterogeneity in primary colorectal cancer and corresponding metastases. Does the apple fall far from the tree? Front Med (Lausanne). 2018;5:234.
Bankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, et al. QuPath: open source software for digital pathology image analysis. Sci Rep. 2017;7(1):16878.
Koressaar T, Lepamets M, Kaplinski L, Raime K, Andreson R, Remm M. Primer3_masker: integrating masking of template sequence with primer design software. Bioinformatics. 2018;34(11):1937–8.
Ijsselsteijn ME, van der Breggen R, Farina Sarasqueta A, Koning F, de Miranda N. A 40-marker panel for high dimensional characterization of cancer immune microenvironments by imaging mass cytometry. Front Immunol. 2019;10:2534.
Eling N, Damond N, Hoch T, Bodenmiller B. Cytomapper: an R/bioconductor package for visualisation of highly multiplexed imaging data. Bioinformatics. 2020.
Kamentsky L, Jones TR, Fraser A, Bray MA, Logan DJ, Madden KL, et al. Improved structure, function and compatibility for Cell Profiler: modular high-throughput image analysis software. Bioinformatics. 2011;27(8):1179–80.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. arXiv preprint arXiv:1207.3907arXiv:1207.3907. 2012.
Jones D, Raine KM, Davies H, Tarpey PS, Butler AP, Teague JW, et al. cgpCaVEManWrapper: simple execution of CaVEMan in order to detect somatic single nucleotide variants in NGS data. Curr Protoc Bioinform. 2016;56:15 0 1- 0 8.
Ye K, Schulz MH, Long Q, Apweiler R, Ning Z. Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics. 2009;25(21):2865–71.
Raczy C, Petrovski R, Saunders CT, Chorny I, Kruglyak S, Margulies EH, et al. Isaac: ultra-fast whole-genome secondary analysis on Illumina sequencing platforms. Bioinformatics. 2013;29(16):2041–3.
Chen X, Schulz-Trieglaff O, Shaw R, Barnes B, Schlesinger F, Kallberg M, et al. Manta: rapid detection of structural variants and indels for germline and cancer sequencing applications. Bioinformatics. 2016;32(8):1220–2.
Roller E, Ivakhno S, Lee S, Royce T, Tanner S. Canvas: versatile and scalable detection of copy number variants. Bioinformatics. 2016;32(15):2375–7.
Saunders CT, Wong WS, Swamy S, Becq J, Murray LJ, Cheetham RK. Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs. Bioinformatics. 2012;28(14):1811–7.
Gonzalez-Perez A, Perez-Llamas C, Deu-Pons J, Tamborero D, Schroeder MP, Jene-Sanz A, et al. IntOGen-mutations identifies cancer drivers across tumor types. Nat Methods. 2013;10(11):1081–2.
Lawrence MS, Stojanov P, Polak P, Kryukov GV, Cibulskis K, Sivachenko A, et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature. 2013;499(7457):214–8.
Martincorena I, Raine KM, Gerstung M, Dawson KJ, Haase K, Van Loo P, et al. Universal patterns of selection in cancer and somatic tissues. Cell. 2018;173(7):1823.
Blokzijl F, Janssen R, van Boxtel R, Cuppen E. MutationalPatterns: comprehensive genome-wide analysis of mutational processes. Genome Med. 2018;10(1):33.
Flensburg C, Sargeant T, Oshlack A, Majewski I. SuperFreq: integrated mutation detection and clonal tracking in cancer. BioRxiv. 2018:380097.
Hundal J, Carreno BM, Petti AA, Linette GP, Griffith OL, Mardis ER, et al. pVAC-Seq: a genome-guided in silico approach to identifying tumor neoantigens. Genome Med. 2016;8(1):11.
Dilthey AT, Mentzer AJ, Carapito R, Cutland C, Cereb N, Madhi SA, et al. HLA*LA - HLA typing from linearly projected graph alignments. Bioinformatics. 2019.
Tian Y, Morris TJ, Webster AP, Yang Z, Beck S, Feber A, et al. ChAMP: updated methylation analysis pipeline for Illumina BeadChips. Bioinformatics. 2017;33(24):3982–4.
Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, et al. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 2019;47(W1):W191–8.
Rospo G, Lorenzato A, Amirouchene-Angelozzi N, Magri A, Cancelliere C, Corti G, et al. Evolving neoantigen profiles in colorectal cancers with DNA repair defects. Genome Med. 2019;11(1):42.
Goodman AM, Kato S, Bazhenova L, Patel SP, Frampton GM, Miller V, et al. Tumor mutational burden as an independent predictor of response to immunotherapy in diverse cancers. Mol Cancer Ther. 2017;16(11):2598–608.
Deng L, Liang H, Xu M, Yang X, Burnette B, Arina A, et al. STING-dependent cytosolic DNA sensing promotes radiation-induced type I interferon-dependent antitumor immunity in immunogenic tumors. Immunity. 2014;41(5):843–52.
Swift LH, Golsteyn RM. Genotoxic anti-cancer agents and their relationship to DNA damage, mitosis, and checkpoint adaptation in proliferating cancer cells. Int J Mol Sci. 2014;15(3):3403–31.
Chang L, Graham PH, Hao J, Ni J, Bucci J, Cozzi PJ, et al. PI3K/Akt/mTOR pathway inhibitors enhance radiosensitivity in radioresistant prostate cancer cells through inducing apoptosis, reducing autophagy, suppressing NHEJ and HR repair pathways. Cell Death Dis. 2014;5:e1437.
Kennedy AL, Morton JP, Manoharan I, Nelson DM, Jamieson NB, Pawlikowski JS, et al. Activation of the PIK3CA/AKT pathway suppresses senescence induced by an activated RAS oncogene to promote tumorigenesis. Mol Cell. 2011;42(1):36–49.
Dolly SO, Wagner AJ, Bendell JC, Kindler HL, Krug LM, Seiwert TY, et al. Phase I study of apitolisib (GDC-0980), dual phosphatidylinositol-3-kinase and mammalian target of rapamycin kinase inhibitor, in patients with advanced solid tumors. Clin Cancer Res. 2016;22(12):2874–84.
Basu B, Krebs MG, Sundar R, Wilson RH, Spicer J, Jones R, et al. Vistusertib (dual m-TORC1/2 inhibitor) in combination with paclitaxel in patients with high-grade serous ovarian and squamous non-small-cell lung cancer. Ann Oncol. 2018;29(9):1918–25.
Narayan RS, Fedrigo CA, Brands E, Dik R, Stalpers LJ, Baumert BG, et al. The allosteric AKT inhibitor MK2206 shows a synergistic interaction with chemotherapy and radiotherapy in glioblastoma spheroid cultures. BMC Cancer. 2017;17(1):204.
Yeh CH, Bellon M, Nicot C. FBXW7: a critical tumor suppressor of human cancers. Mol Cancer. 2018;17(1):115.
Zhang Q, Karnak D, Tan M, Lawrence TS, Morgan MA, Sun Y. FBXW7 facilitates nonhomologous end-joining via K63-linked polyubiquitylation of XRCC4. Mol Cell. 2016;61(3):419–33.
Zhang Q, Mady ASA, Ma Y, Ryan C, Lawrence TS, Nikolovska-Coleska Z, et al. The WD40 domain of FBXW7 is a poly(ADP-ribose)-binding domain that mediates the early DNA damage response. Nucleic Acids Res. 2019;47(8):4039–53.
Noorbakhsh J, Kim H, Namburi S, Chuang JH. Distribution-based measures of tumor heterogeneity are sensitive to mutation calling and lack strong clinical predictive power. Sci Rep. 2018;8(1):11445.
Ganesh K, Wu C, O’Rourke KP, Szeglin BC, Zheng Y, Sauve CG, et al. A rectal cancer organoid platform to study individual responses to chemoradiation. Nat Med. 2019;25(10):1607–14.
Nagle PW, Coppes RP. Current and Future perspectives of the use of organoids in radiobiology. Cells. 2020;9(12).
We acknowledge the support of the Genomics Birmingham Sequencing Facility (of which ADB is academic lead).
This study was supported by grants from the Academy of Medical Sciences and the Wellcome Trust (Ref 102732/Z/13/Z). ADB is currently supported by a Cancer Research UK Advanced Clinician Scientist award (Ref C31641/A23923).
Ethical approval and consent to participate
NW Research Ethics committee (15/NW/0079). All patients gave informed consent to participate.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Supplementary data from mutational analysis (MutSigCV, Intogen, Intogen Pathways and COSMIC mutational signatures).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Stockton, J.D., Tee, L., Whalley, C. et al. Complete response to neoadjuvant chemoradiotherapy in rectal cancer is associated with RAS/AKT mutations and high tumour mutational burden. Radiat Oncol 16, 129 (2021). https://doi.org/10.1186/s13014-021-01853-y