Increased Levels of Genomic Instability and Mutations in Homologous Recombination Genes in Locally Advanced Rectal Carcinomas

Pre-operative 5-fluoracil-based chemoradiotherapy (nCRT) is the standard treatment for patients with locally advanced rectal cancer (LARC). Patients with pathological complete response (pCR–0% of tumor cells in the surgical specimen after nCRT) have better overall survival and lower risk of recurrence in comparison with incomplete responders (pIR). Predictive biomarkers to be used for new therapeutic strategies and capable of stratifying patients to avoid overtreatment are needed. We evaluated the genomic profiles of 33 pre-treatment LARC biopsies using SNP array and targeted-next generation sequencing (tNGS). Based on the large number of identified genomic alterations, we calculated the genomic instability index (GII) and three homologous recombination deficiency (HRD) scores, which have been reported as impaired DNA repair markers. We observed high GII in our LARC cases, which was confirmed in 165 rectal cancer cases from TCGA. Patients with pCR presented higher GII compared with pIR. Moreover, a negative correlation between GII and the fraction of tumor cells remaining after surgery was observed (ρ = –0.382, P = 0.02). High HRD scores were detected in 61% of LARC, of which 70% were incomplete responders. Using tNGS (105 cancer-related genes, 13 involved in HR and 5 in mismatch repair pathways), we identified 23% of cases with mutations in HR genes, mostly in pIR cases (86% of mutated cases). In agreement, the analysis of the TCGA dataset (N = 145) revealed 21% of tumors with mutations in HR genes. The HRD scores were shown to be predictive of better response to PARP-inhibitors and platinum-based chemotherapy in breast and ovarian cancer. Our results suggest that the same strategy could be applied in a set of LARC patients with HRD. In conclusion, we identified high genomic instability in LARC, which was related to alterations in the HR pathway, especially in pIR. These findings suggest that patients with impaired HRD would clinically benefit from PARP-inhibitors and platinum-based therapy.

Introduction

Several studies have demonstrated that colorectal cancer (CRC) is a heterogeneous disease and that rectal and colon tumors differ in many aspects (1). Among them, risk factors (gender, age, body mass index, dietary association), mutational profiling (most rectal cancers are negative for mutations in mismatch repair genes, while other genes are poorly explored), microsatellite instability (MSI is rare in rectal cancer) and CpG island methylator phenotype (CIMP is 0-low in rectal cancer) have been described (2, 3).

Chromosomal instability (CIN) and MSI are well-recognized mechanisms in CRC. CIN results from high levels of chromosomal mis-segregation, defects on telomere stability and/or DNA damage response (4). Although several approaches to quantify CIN have been described, the thresholds to consider CIN-positive are not standardized (5, 6).

Large DNA copy number alterations (named genomic scars) have been associated with the inability to repair DNA double strand breaks (DSB), which is the most hazardous type of DNA damage. The DSB are mostly repaired by non-homologous end joining (NHEJ) and homologous recombination (HR) repair mechanisms. The ability to predict HR deficiency (HRD), and consequently impaired DNA damage repair, can be determinant of response to chemotherapeutic agents, as reported in breast cancer (7–9). Scores used to capture large-scale transition (LST), telomeric allelic imbalance (tAI), and homologous recombination deficiency-loss of heterozygosity (HRD-LOH) have been used as tools to predict HRD (7–11). Overall, tumors with high scores are very likely to carry a deficient DNA damage repair system, resulting in a better response to platinum-based therapy (10).

The management of patients with locally advanced rectal carcinomas (LARC) with neoadjuvant 5-fluorouracil-based chemo and radiotherapy (nCRT) followed by total mesorectal excision has resulted in reduced locoregional recurrence and death by the disease. Cases achieving pathological complete response (pCR) to nCRT have shown lower rates of local and distant recurrences, and better survival compared to patients with incomplete response (pIR) (12). However, 70–90% of the patients present partial response and ~ 20% show resistance to treatment (13). In addition, the nCRT treatment results in significant morbidity, including long-term side effects and a high rate of clinical complications after surgery (14–18). Recent results from the German CAO/ARO/AIO-04 study showed that adding oxaliplatin to the current nCRT can result in improved disease-free survival (19). However, no molecular characterization was used to select the patients included in this strategy of treatment.

Although few studies described the genomic profile of LARC, the association of nCRT response with chromosomal imbalances and genomic instability has been poorly explored (20, 21). The first study reporting such association used comparative genomic hybridization (CGH) (22). The authors showed higher average of chromosomal imbalances in responders compared with non-responders. Subsequently, array-CGH data revealed higher, but not significant, chromosomal instability in resistant tumors (20, 23). Scores of HRD have not yet been reported in LARC, but an elevated chromosomal missegregation and loss of MRE11 function, a key player of the homologous recombination repair pathway, were associated with better pathological response to chemoradiotherapy in these patients (21).

Based on the current protocol used to treat LARC patients, it is mandatory to identify tumor inherent characteristics capable of stratifying patients according to response to therapy. Moreover, new therapeutic strategies associated with lower morbidity, recurrence rates and better overall survival are needed. Herein, we performed target enrichment for next generation sequencing and genome-wide SNP array in 33 LARC samples carefully selected aiming to verify the predictive value of chromosomal instability based in scores as biomarkers of nCRT response.

Materials and Methods

Patients

Clinical data from all patients with locally advanced rectal cancer admitted at the A. C. Camargo Cancer Center and Barretos Cancer Hospital, São Paulo Brazil, from 2006 to 2015, were carefully revised. Patients diagnosed with LARC by clinical (according to the AJCC 7th ed.) (24) and image evaluations referred to neoadjuvant treatment with continuous infusion of 5-fluouracil or oral capecitabine and radiotherapy (total dose of 50.4 Gy) followed by surgery were included in this study. In addition, the presence of previous cancer and/ or metastases at diagnosis was an exclusion criterion. From this cohort, a large number of patients were excluded due to incomplete clinical data, different treatment strategies or unavailable biopsies. According to these criteria, 33 unrelated patients were included in this study. The biopsy specimens were obtained during colonoscopy prior to pre-operative chemoradiotherapy and evaluated by two specialized pathologists (MDFSB and CSN). Comprehensive clinical, pathological and epidemiological data were obtained from the medical records. The histopathological analysis was revised and a tumor area with more than 70% tumor fraction was selected for molecular evaluation. The Human Research Ethics Committee from both Institutions approved the study (Protocols 1884/14 and 1030/2015, respectively). All patients provided written informed consent prior to sample collection.

Classification of Response

Pathological complete response (pCR; ypT0N0) or pathological incomplete response (pIR) groups were selected according to the absence or presence of reminiscent viable tumor cells in the surgical specimens, respectively. Tumor regression grade (TRG) was scored according to Dworak and Mandard system with modifications (25). The cases were assigned into four categories according to the percentage of viable tumor cells: Grade 0: complete tumor regression (0% of residual tumor cells); Grade I: subtotal tumor regression (< 25% of viable tumor cells); Grade II: partial tumor regression (25–50% of viable tumor cells); Grade III: minimal tumor regression (> 50% of viable tumor cells). Samples with complete tumor regression and showing lymph nodes (ypT0N+) were classified as Grade II (25).

Genomic Profiling

DNA was extracted from frozen tumor specimens using phenol-chloroform-isoamyl alcohol (25:24:1) solution. Genomic imbalances were assessed using the CytoScan HD array (Affymetrix, Santa Clara, CA, USA), according to the manufacturer recommendations. The scanning (Scanner 3000 7G, Affymetrix) and [.CEL] files were generated (Affymetrix^® GeneChip^® Command Console^® Software v.4.0) followed by data analysis Chromosome Analysis Suite (ChAS v.3.1, Affymetrix). Significant alterations were considered with at least 25 probes altered for losses, 50 for gains, and a minimum of 5 Mb for cnLOH (copy-neutral loss of heterozygosity), as previously described (26). Regions with common variants (>1% of the population study) identified using the Database of Genomic Variants (DGV, http://dgv.tcag.ca/dgv/app/home, updated in May 2016) and Affymetrix Database of Variants (composed by 2,421 health individuals evaluated with the Cytoscan HD array) were excluded from subsequent analyses.

Genomic Instability Index and Homologous Recombination Deficiency Scores

The fraction of the genome affected by copy-number changes (Genomic Instability Index, GII) was calculated for each individual case, as previously described (27). The median GII values were compared according to therapy response (pCR vs. pIR, and TRG). Considering that the fraction of tumor cells remaining in the surgical specimen after nCRT is a surrogate of treatment response, we evaluated the correlation between the GII and the % of viable tumor cells identified after surgery.

The homologous recombination deficiency scores (LST, tAI, and HRD-LOH) were calculated as previously described (28). In summary, LST indicates the number of chromosomal breaks between adjacent regions of at least 10 Mb (high: >15 in diploid tumors and > 20 in polyploidy tumors); tAI is referred as the number of subtelomeric regions with allelic imbalance that start beyond the centromere and extended to the telomere region (high: > median value). HRD-LOH indicates the number of LOH regions larger than 15 Mb and shorter than the whole chromosome (high: >10).

Target Enrichment—Next Generation Sequencing (tNGS)

Targeted NGS was performed using a 105 cancer-related genes panel (all exons, 3′UTR and 5′UTR) including 13 genes involved in the HR pathway and five MMR genes (SureSelectXT Custom Panel, Agilent Technologies, Inc., Santa Clara, CA) (Table S1). The libraries were prepared using SureSelectQXT Library Prep Kit (Agilent) according to manufacturer's instructions and sequenced on NextSeq 550 (Illumina, San Diego, CA). For each patient, 2 μL of genomic DNA (25 ng/μL) was used for enzymatic fragmentation. Library amplification was performed using Herculase II Fusion DNA Polymerase (Agilent) and the PCR product was purified using the Agencourt AMPureXP purification bead system (Beckman Coulter; Pasadena, CA). The targeted DNA was captured using streptavidin-coated magnetic beads (Dynabeads MyOne Streptavidin T1; Thermo Fisher Scientific, Waltham, MA) followed by indexing (SureSelectQXT P7 and P5 dual indexing primers). The analysis of amplified indexed library DNA was performed using High Sensitivity D1000 ScreenTape (on Agilent TapeStation). Two samples were excluded (DNA quantity and quality) and 31 were multiplexed into 1.4 pM pool and loaded onto the NextSeq 550 (Illumina).

tNGS Data Analysis

The raw sequencing reads demultiplexing and FASTQ files were obtained, with the resulting reads being mapped to the UCSC human genome reference build 19 using the BWA alignment algorithm (29). Variant calling and quality filtering were performed with GATK (30). Variant annotation was primarily done using ANNOVAR (31), GATK and SnpEff (32). Ingenuity software (Qiagen, Hilden, Germany) was used to classify the variants according to the American College of Medical Genetics and Genomics (ACMG) (33). The pipeline for variant filtering is represented in Figure S1. Shortly, benign or likely benign variants (ACMG) and those with frequency > 0.01 in the GAD (34) and/or in the ExAC (35) were excluded. The remaining variants were manually curated using the Genome Browse software (Golden Helix Inc., Bozeman, MT). The variants mapped in homopolymer regions were excluded and a new filtering was applied to exclude variants of uncertain significance (VUS, classified by ACMG) with CADD score < 3.0.

Cross-Validation Analysis

Genome wide SNP6 segmented somatic copy number alterations data from 165 rectal cancers were retrieved from TCGA (available at Xena Browser, accessed on March 2018: https://tcga.xenahubs.net/download/TCGA.READ.sampleMap/SNP6_nocnv_genomicSegment.gz). Segments were filtered (log2 ratio of 0.2) and GII was calculated as described above (27). The mutational profile of 145 of 165 rectal cancer cases were available (cBioPortal; accessed on April 2018: http://www.cbioportal.org/) and used to evaluate the mutational status of 23 MMR genes (KEGG: hsa03430) and 41 HR genes (KEGG: hsa03440) in association to GII.

Statistical Analysis

Statistical analyses were performed using the SPSS (SPSS 24 Inc, Chicago, IL, USA). Mann-Whitney U and Fisher exact test were applied to investigate association among the clinical parameters, GII, HRD-scores and nCRT response. The non-parametric Kruskall-Wallis test was applied to compare multiple groups. The Spearman correlation coefficient (ρ) was used to assess the relationship between two variables. A comparison of segmented CNAs data between groups according to the response to nCRT was performed using the CoNVaQ web tool (36). Computation of empirical p-values for permutation testing was performed for finding a matching region of the same length or longer, when the size of each group is preserved, but the samples are randomly distributed among groups (36).

Results

Patient Characteristics and Response to nCRT

Clinical parameters and response to nCRT of patients included in this study are summarized in Table 1. Eleven (33%) patients achieved pCR (or TRG Grade 0) and 22 (67%), pIR (TRG 1: 10, TRG 2: 9, and TRG 3: 3). The analysis of clinical features and pathological response showed no statistical significance (data not shown).

TABLE 1

Table 1. Clinical and histopathological characteristics of 33 locally advanced rectal cancer patients included in this study.

Copy Number Alterations and Chromosomal Instability Analysis

A total of 1,442 CNAs was detected including 880 gains, 442 losses and 120 cnLOH (596 in 11 pCR; 846 in 22 pIR), ranging from 1 to 161 per case. Gains of 7p, 8q, 13q, and 20q and losses of 8p, 14, 17p, and 18 were frequently detected. The statistical model of CoNVaQ (36) pointed out 45 significant chromosomal regions mapped at chromosomes 4, 8, 19, 20, and 21 (p < 0.05) associated with therapeutic response (Table S2). Loss of chromosome 21 (start 35020505—end 40248388) was found only in pIR cases (36.4%) (permutation p-value > 0.29).

A similar GII distribution was found comparing our LARC cases (N = 33, median = 0.358) with the TCGA dataset (N = 165, median = 0.284) (P > 0.05). The GII median value was used as threshold (N = 198, median = 0.286). Eighteen (55%; 11/22 pIR and 7/11 pCR, Figure 2) and 81 (49%) rectal tumors from our dataset and TCGA, respectively, presented high GII levels. Although statistically non-significant, pCR presented higher median GII (0.475) compared to incomplete responders (0.294). Similar analysis using the TRG revealed an opposite trend between the median GII and TRG. A significant difference was found grouping TRG 0+1 and TRG 2+3 (P = 0.043) (Figure 1A). A negative correlation between GII and the fraction of tumor cells remaining after surgery was observed (ρ = −0.382, P = 0.02). Information for TCGA cases regarding the response to nCRT was not available, preventing further comparisons.

FIGURE 1

Figure 1. (A) Distribution of the Genomic Instability Index (GII) and (B) the number of mutations found in rectal cancer samples according to pathological response (pCR: complete response or pIR: incomplete) to neoadjuvant chemoradiotherapy. Cases were also categorized into tumor regression grade (TRG) and grouped in TRG 0 + 1 and TRG 2 + 3. (C) Homologous recombination deficiency scores (LST, LOH, and tAI) and GII of samples with mutation in DNA repair genes. (D) GII of each TCGA and internal cases carrying mutation in DNA damage repair pathways (homologous recombination and mismatch repair) genes. Samples carrying mutation in the same gene are shown connected by a line, additional genes are shown below each sample. Samples with >1000 mutations are indicated, and specific genes were not taken into consideration for these cases. Five of 14 cases from TCGA with mutation in ATM presented >1000 mutations, and nine of them are represented.

FIGURE 2

Figure 2. Proportion of locally advanced rectal cases with pathological complete response (pCR) or incomplete response (pIR) to neoadjuvant therapy according to the genomic features evaluated in this study. Patients were divided into high or low levels of Genomic Instability Index (GII), elevated values of at least one of the three scores of Homologous Recombination Repair Deficiency (HRD scores), and mutation in genes related to DNA damage repair (mismatch repair-MMR or HR repair-HRR).

Homologous Recombination Deficiency Scores

The HRD score on median values of tAI, LST and HRD-LOH were 14, 6, and 2, respectively. Using median values > 14.0 as cut-off for tAI, 14 of 22 pIR and 4 of 11 pCR presented high scores. Two pCR cases showed high LST (>15 or >20 for diploid or polyploid tumors) and five cases (3pIR and 2pCR) presented HRD-LOH scores higher than 10. Potential deficiency in HR pathway (at least one high score) was found in 55% of pCR and 64% of pIR cases (Figure 2). No significant differences were found according to the nCRT and clinical features (Table 2, Figure S2).

TABLE 2

Table 2. The number of cases with high Genomic Instability Index (GII) and predicted homologous recombination deficiency (HRD) based on three scores (tAI, LOH, and LST).

The Spearman correlation coefficient between each pair of scores was significant for LST and tAI (ρ = 0.720, p < 0.001) scores. Significant correlation between GII and tAI was detected for pIR cases (ρ = 0.470, p = 0.02), while no correlation was observed between the scores and GII for pCR cases. HRD scores and GII for each case are depicted in Table S3.

Mutational Profile of Rectal Cancer

Targeted-NGS analysis of 31 cases revealed 161 unique variants in 51 genes, with an average of five mutations per case in pCR, and six in pIR (Figure 1B, Table S4). Thirteen genes were mutated regardless of response to therapy, including TP53 (84%), APC (81%) and KRAS (45%). Eight genes were mutated exclusively in pCR cases and 30 genes in pIR. The frequency of mutated genes was similar to the TCGA data (Fisher's, P > 0.05), except for five genes with higher frequency in our cases (FLT1, FLT4, MMP1, RAD51C, WNT1) (Figure S3).

Seven cases (23%) presented HR genes mutation, including ATM (2pIR), BRIP1 (1pIR), POLD1 (1pIR), RAD51B (1pCR), and RAD51C (2pIR). Four cases (13%) harbored MMR gene mutations: MLH1 (1pCR), MLH3 (1pIR) and MSH6 (2pIR) (Figure 2). Similarly, the TCGA dataset (N = 145) presented 21% of tumors with HR gene mutations and 10% with MMR gene mutations. Four cases with HR gene mutations also showed high HRD scores (Figure 1C).

Next generation sequencing data from 145 rectal cancer retrieved from TCGA were used to evaluate the mutational status of 23 MMR genes (KEGG: hsa03430) and 41 HR genes (KEGG: hsa03440) in association to GII. Cases with GII higher than the threshold (median of 0.286) presented mutations almost exclusively in HR genes. Only one case carried MMR mutation (RFC1). Low GII values (median = 0.01) were detected in cases presenting more than 1,000 mutations. Six of 14 ATM mutated cases had high GII (median = 0.35), while five, that were also hypermutated cases (>1,000 including MMR genes), presented low GII (median = 0.01). The GII values from TCGA data and our internal dataset according to the mutation status of HR/MMR genes are depicted in Figure 1D.

Correlation of Genomic Alterations and Response to nCRT

The genomic variables (GII, HRD scores and mutational profile) were categorized in high and low or mutated and wild type (WT) to compare pCR with pIR cases (Figure 2). This analysis revealed slight difference between the two groups. The pCR cases presented increased proportion of high GII, but lower rate of cases with at least one high HRD score and mutations in MMR and HR genes compared with pIR group.

Discussion

Genomic alterations and response to therapy in LARC are insufficiently understood and explored. In our study, LARC cases presented a high level of chromosomal instability recapitulating the most common alterations previously reported in literature (20, 23, 37). In addition, loss of chromosome 21, where RUNX1 is mapped, was detected exclusively in pIR cases (36.4%). This tumor suppressor encodes a transcription factor downregulated in several cancers (38).

We calculated the medium values of the genomic instability index (GII) based on the findings obtained in the cohort of 165 rectal cancer from TCGA and our set of 33 LARC (median GII of 0.286). Comparing with the pan-cancer analysis of genomic scar signatures, this is higher than described in prostate adenocarcinoma (0.081) and lower than breast (0.435) and ovarian serous carcinoma (0.574) (TCGA data) (28).

To investigate the presence of impaired DNA damage repair, three scores used to predict HRD were evaluated by capturing tAI, HRD-LOH and LST (7, 10, 11). Interestingly, high median scores were detected for tAI (14), but low for LST (6) and HRD-LOH (2). Marquard et al. (28) evaluated HRD scores in 15 different tumor types, including colorectal adenocarcinomas. Based on the average of these three scores, colon cancer was ranked at 10th position while our LARC was at 7th position (below breast cancer). High HRD scores were associated to increased sensitivity to platinum-based chemotherapy and PARP-inhibitors in breast and ovarian cancer (9). A considerable subset of our cases presented high HRD scores (61%; of which 70% were pIR) suggesting that these patients could benefit from platinum-based therapy.

The correlation among the three HRD scores revealed a positive significant value between tAI and LST (ρ = 0.711). In agreement, Marquard et al. showed a correlation between these two scores (ρ > 0.7 in 8 of 15 cancer types) and that HRD-LOH presented lower correlation with tAI and LST in 13 of 15 tumors (28). Overall, no significant correlation between GII and these three scores was found. However, GII was positively correlated with tAI (ρ = 0.470, p = 0.02) in pIR cases, suggesting that GII in rectal cancer quantify different aspects of CNA, being probably associated with HRD. Alternative mechanisms could be involved in pCR cases. In fact, the correlation between GII and HRD scores in ovarian cancer was much lower than the median value for all cancer types and more similar to our findings (28).

Interestingly, mutations in genes involved in HR repair were detected in 23% of our LARC cases and in 21% of TCGA cohort, respectively. Six of the seven cases with mutations in one of these genes (ATM, BRIP1, POLD1, RAD51B, RAD51C) presented incomplete response to therapy. We detected an elevated tAI score and GII in two cases presenting RAD51C mutation. In cancer cells, mutated RAD51C recruits non-homologous end joining (NHEJ) proteins, which correlates to higher genomic instability due to error-prone repair. Deficiency of RAD51C confers enhanced sensitivity to Poly (ADP-ribose) Polymerase 1 (PARP1) inhibitors in combination with radiation, which can be used as an alternative synergistic treatment approach for 5-FU resistant tumors (39). The synthetic lethality achieved with PARP1 inhibitors was also observed in pancreatic ductal adenocarcinoma cells with ATM mutation. Two LARC cases with pIR presented ATM mutations. ATM deficiency promotes genomic instability by blunting DNA damage repair and leading to complex genomic aberration such as chromothripsis (40), which we detected in one ATM mutated case (chromosome 20).

The mutational profiling of 145 cases from TCGA revealed ATM mutation in 10% of the cases; 43% of them presented high GII. Nonetheless, this correlation was stronger when excluding cases with more than 1,000 mutations (78% with high GII). All five ATM mutated cases that also carried more than 1,000 mutations, including in MMR genes, had very low GII. Mutations found in MMR and HR genes in such cases can be a consequence of their hypermutated phenotype and should not be considered for an association with GII.

Although the initial results of clinical trials evaluating the addition of oxaliplatin to the conventional 5-FU based nCRT regimen have shown no additional benefits (41), an increased disease-free survival was observed after long-term follow-up (19). Compared with older patients (aged ≥ 70 years), those under 60 years of age presented a reduced number of recurrences, metastasis and better overall survival (42). However, the HR pathway was not assessed, and this analysis could be predictive of response to oxaliplatin.

No association between the number of mutations or a specific mutated gene and the complete response to nCRT was observed. The most frequently mutated genes were TP53, APC, and KRAS. Wild type TP53 was associated with a good response to nCRT (Risk Ratio = 1.30, 95% CI = 1.14–1.49, P < 0.001) in a meta-analysis with 1,830 rectal cancer patients from 30 studies (43). However, the lack of standard treatment regimen, the definition of methods, their quantification, and outcome may lead to bias in the meta-analysis. The predictive power of KRAS in the nCRT response was also reported, in which codon 13 mutation showed association to a worse response to therapy (44, 45). Furthermore, the small number of LARC evaluated in several studies, including our own, make the use of these mutations as predictive biomarkers still controversial.

Among the genes mutated only in our pCR cases, two members of the human ATP-binding cassete (ABC) transporter superfamily present compelling evidence of a role in chemoresistance. ABCG2 is expressed in gastrointestinal tissue and its inhibition was shown to increase the sensitivity to 5-FU in vitro and in vivo (46). Inhibition of ABCC1 expression has been reported as leading to 5-FU enhanced response in CRC cell lines (47). This gene is transcriptionally regulated by MYC, which was also found mutated in our pCR case (48). In breast cancer and glioblastoma cell lines, the downregulation of MYC protein resulted in suppression of ABCC1, being a putative mechanism to sensitize tumors to chemotherapy (48). Although neoadjuvant treatment information from the TCGA cohort is unavailable, a comparison with our data revealed similar frequencies of these ABC gene mutations (ABCC1: 3% in our data and 6%: TCGA; ABCG2: 3% in both).

Although we presented new and relevant data, a set of limitations could be pointed out including the small sample size, a challenge also faced by the previous large-scale genomic studies in LARC (20, 22, 23). Furthermore, the use of single pre-treatment biopsies precludes the evaluation of intra-tumoral heterogeneity, already described in a set of rectal tumors (49). As tumors can be an assembly of different subclones, intra-tumor heterogeneity can also play a role in tumor response to therapy. However, the implication of heterogeneity in the prediction and prognosis of rectal cancer has yet to be explored.

In conclusion, we evaluated the genomic profile of LARC in association with the response to nCRT and quantified the observed high chromosomal instability into HRD scores and GII. Similar to breast and ovarian cancer, the HRD scores were high in LARC and further studies are necessary to evaluate its association with the sensitivity to platinum-based chemotherapy and PARP-inhibitors. The assessment of HRD can be used to select patients for neoadjuvant chemotherapy strategy. A set of tumors with high genomic instability (GII) can be explained by mutations in genes involved in DNA repair by homologous recombination. Alteration in these genes was observed in 23% of our cases and in 21% of the TCGA dataset, revealing promising results as new targets for therapy alone or in combination with the current nCRT.

Ethics Statement

This study was carried out in accordance with the recommendations of Human Research Ethics Committees from A. C. Camargo Cancer Center and Barretos Cancer Hospital, São Paulo Brazil. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by both Committees (Protocols 1884/14 and 1030/2015, respectively).

Author Contributions

SR and SA conceived and designed the study. MB and CS-N performed histopathological evaluation. SA, BC, LC, and SR selected the cases and obtained the clinical data. LC conducted the experiments. LC, SL, JB, AP, and MA analyzed the data. LC and SR wrote and edited the manuscript. All authors read and approved the final version of the manuscript.

Funding

National Institute of Science and Technology in Oncogenomics—INCITO: São Paulo Research Foundation—FAPESP (grant number 2008/57887-9) and National Council for Scientific and Technological Development—CNPq (grant number 573589/08-9); Danish Colorectal Cancer Center South, Denmark; FAPESP (grant numbers 2014/06323-9, 2015/25803-4).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors would like to thank the patients who agreed to participate in this study, Mateus Camargo de Barros Filho and Fabio Albuquerque Marchi for their contribution during the study. We also thank the Biobank and Department of Pathology of A.C. Camargo Cancer Center and Barretos Cancer Hospital (São Paulo, Brazil).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.00395/full#supplementary-material

References

1. Slattery ML, Curtin K, Wolff RK, Boucher KM, Sweeney C, Edwards S, et al. A comparison of colon and rectal somatic DNA alterations. Dis Colon Rectum. (2009) 52:1304–11. doi: 10.1007/DCR.0b013e3181a0e5df

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Kapiteijn E, Liefers GJ, Los LC, Kranenbarg EK, Hermans J, Tollenaar RA, et al. Mechanisms of oncogenesis in colon versus rectal cancer. J Pathol. (2001) 195:171–8. doi: 10.1002/path.918

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Fenebro E, Halvarsson B, Baldetorp B, Nilbert M. Predominance of CIN versus MSI in the development of rectal cancer at young age. BMC Cancer. (2009) 2:25. doi: 10.1186/1471-2407-2-25

CrossRef Full Text | Google Scholar

4. Pino MS, Chung DC. The chromosomal instability pathway in colon cancer. Gastroenterology. (2010) 138:2059–72. doi: 10.1053/j.gastro.2009.12.065

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Watanabe T, Kobunai T, Yamamoto Y, Matsuda K, Ishihara S, Nozawa K, et al. Chromosomal instability (CIN) phenotype, CIN high or CIN low, predicts survival for colorectal cancer. J Clin Oncol. (2012) 30:2256–64. doi: 10.1200/JCO.2011.38.6490

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Burrell RA, Mcclelland SE, Endesfelder D, Groth P, Weller MC, Shaikh N, et al. Replication stress links structural and numerical cancer chromosomal instability. Nature. (2013) 494:492–6. doi: 10.1038/nature11935

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Abkevich V, Timms KM, Hennessy BT, Potter J, Carey MS, Meyer LA, et al. Patterns of genomic loss of heterozygosity predict homologous recombination repair defects in epithelial ovarian cancer. Br J Cancer. (2012) 107:1776–82. doi: 10.1038/bjc.2012.451

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wang ZC, Birkbak NJ, Culhane AC, Drapkin R, Fatima A, Tian R, et al. Profiles of genomic instability in high-grade serous ovarian cancer predict treatment outcome. Clin Cancer Res. (2012) 18:5806–15. doi: 10.1158/1078-0432.CCR-12-0857

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Watkins JA, Irshad S, Grigoriadis A, Tutt AN. Genomic scars as biomarkers of homologous recombination deficiency and drug response in breast and ovarian cancers. Breast Cancer Res. (2014) 16:211. doi: 10.1186/bcr3670

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Birkbak NJ, Wang ZC, Kim JY, Eklund AC, Li Q, Tian R, et al. Telomeric allelic imbalance indicates defective DNA repair and sensitivity to DNA-damaging agents. Cancer Discov. (2012) 2:366–75. doi: 10.1158/2159-8290.CD-11-0206

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Popova T, Manie E, Rieunier G, Caux-Moncoutier V, Tirapo C, Dubois T, et al. Ploidy and large-scale genomic instability consistently identify basal-like breast carcinomas with BRCA1/2 inactivation. Cancer Res. (2012) 72:5454–62. doi: 10.1158/0008-5472.CAN-12-1470

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zorcolo L, Rosman AS, Restivo A, Pisano M, Nigri GR, Fancellu A, et al. Complete pathologic response after combined modality treatment for rectal cancer and long-term survival: a meta-analysis. Ann Surg Oncol. (2012) 19:2822–32. doi: 10.1245/s10434-011-2209-y

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Fokas E, Liersch T, Fietkau R, Hohenberger W, Beissbarth T, Hess C, et al. Tumor regression grading after preoperative chemoradiotherapy for locally advanced rectal carcinoma revisited: updated results of the CAO/ARO/AIO-94 trial. J Clin Oncol. (2014) 32:1554–62. doi: 10.1200/JCO.2013.54.3769

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Loos M, Quentmeier P, Schuster T, Nitsche U, Gertler R, Keerl A, et al. Effect of preoperative radio(chemo)therapy on long-term functional outcome in rectal cancer patients: a systematic review and meta-analysis. Ann Surg Oncol. (2013) 20:1816–28. doi: 10.1245/s10434-012-2827-z

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Tevis SE, Kohlnhofer BM, Stringfield S, Foley EF, Harms BA, Heise CP, et al. Postoperative complications in patients with rectal cancer are associated with delays in chemotherapy that lead to worse disease-free and overall survival. Dis Colon Rectum. (2013) 56:1339–48. doi: 10.1097/DCR.0b013e3182a857eb

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Giandomenico F, Gavaruzzi T, Lotto L, Del Bianco P, Barina A, Perin A, et al. Quality of life after surgery for rectal cancer: a systematic review of comparisons with the general population. Expert Rev Gastroenterol Hepatol. (2015) 9:1227–42. doi: 10.1586/17474124.2015.1070667

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Jiang DM, Raissouni S, Mercer J, Kumar A, Goodwin R, Heng DY, et al. Clinical outcomes of elderly patients receiving neoadjuvant chemoradiation for locally advanced rectal cancer. Ann Oncol. (2015) 26:2102–6. doi: 10.1093/annonc/mdv331

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Hupkens BJP, Martens MH, Stoot JH, Berbee M, Melenhorst J, Beets-Tan RG, et al. Quality of life in rectal cancer patients after chemoradiation: watch-and-wait policy versus standard resection - a matched-controlled study. Dis Colon Rectum. (2017) 60:1032–40. doi: 10.1097/DCR.0000000000000862

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Rödel C, Graeven U, Fietkau R, Hohenberger W, Hothorn T, Arnold D, et al. Oxaliplatin added to fluorouracil-based preoperative chemoradiotherapy and postoperative chemotherapy of locally advanced rectal cancer (the German CAO/ARO/AIO-04 study): final results of the multicentre, open-label, randomised, phase 3 trial. Lancet Oncol. (2015) 16:979–89. doi: 10.1016/S1470-2045(15)00159-X

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Molinari C, Ballardini M, Teodorani N, Giannini M, Zoli W, Emiliani E, et al. Genomic alterations in rectal tumors and response to neoadjuvant chemoradiotherapy: an exploratory study. Radiat Oncol. (2011) 6:161. doi: 10.1186/1748-717X-6-161

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zaki BI, Suriawinata AA, Eastman AR, Garner KM, Bakhoum SF. Chromosomal instability portends superior response of rectal adenocarcinoma to chemoradiation therapy. Cancer. (2014) 120:1733–42. doi: 10.1002/cncr.28656

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Grade M, Gaedcke J, Wangsa D, Varma S, Beckmann J, Liersch T, et al. Chromosomal copy number changes of locally advanced rectal cancers treated with preoperative chemoradiotherapy. Cancer Genet Cytogenet. (2009) 193:19–28. doi: 10.1016/j.cancergencyto.2009.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Chen Z, Liu Z, Li W, Qu K, Deng X, Varma MG, et al. Chromosomal copy number alterations are associated with tumor response to chemoradiation in locally advanced rectal cancer. Genes Chromosomes Cancer. (2011) 50:689–99. doi: 10.1002/gcc.20891

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Sobin L, Gospodarowicz M, Wittekind C. TNM Classification of Malignant Tumours. 7th ed. Oxford: Wiley-Blackwell (2011).

Google Scholar

25. Santos MD, Silva C, Rocha A, Matos E, Nogueira C, Lopes C. Prognostic value of mandard and dworak tumor regression grading in rectal cancer: study of a single tertiary center. ISRN Surg. (2014) 2014:310542. doi: 10.1155/2014/310542

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Villacis RA, Miranda PM, Gomy I, Santos EM, Carraro DM, Achatz MI, et al. Contribution of rare germline copy number variations and common susceptibility loci in Lynch syndrome patients negative for mutations in the mismatch repair genes. Int J Cancer. (2016) 138:1928–35. doi: 10.1002/ijc.29948

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Liu C, Srihari S, Lal S, Gautier B, Simpson PT, Khanna KK, et al. Personalised pathway analysis reveals association between DNA repair pathway dysregulation and chromosomal instability in sporadic breast cancer. Mol Oncol. (2016) 10:179–93. doi: 10.1016/j.molonc.2015.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Marquard AM, Eklund AC, Joshi T, Krzystanek M, Favero F, Wang ZC, et al. Pan-cancer analysis of genomic scar signatures associated with homologous recombination deficiency suggests novel indications for existing cancer drugs. Biomark Res. (2015) 3:9. doi: 10.1186/s40364-015-0033-4

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. (2009) 25:1754–60. doi: 10.1093/bioinformatics/btp324

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Depristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. (2011) 43:491–8. doi: 10.1038/ng.806

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. (2010) 38:e164. doi: 10.1093/nar/gkq603

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Cingolani P, Platts A, Wang Le L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. (2012) 6:80–92. doi: 10.4161/fly.19695

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. (2015) 17:405–24. doi: 10.1038/gim.2015.30

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Becker KG, Barnes KC, Bright TJ, Wang SA. The genetic association database. Nat Genet. (2004) 36:431–2. doi: 10.1038/ng0504-431

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. (2016) 536:285–91. doi: 10.1038/nature19057

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Larsen SJ, Do Canto LM, Rogatto SR, Baumbach J. CoNVaQ: a web tool for copy number variation-based association studies. BMC Genomics. (2018) 19:369. doi: 10.1186/s12864-018-4732-8

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Gonzalez-Gonzalez M, Garcia J, Alcazar JA, Gutierrez ML, Gonzalez LM, Bengoechea O, et al. Association between the cytogenetic profile of tumor cells and response to preoperative radiochemotherapy in locally advanced rectal cancer. Medicine. (2014) 93:e153. doi: 10.1097/MD.0000000000000153

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Chuang LS, Ito K, Ito Y. RUNX family: Regulation and diversification of roles through interacting proteins. Int J Cancer. (2013) 132:1260–71. doi: 10.1002/ijc.27964

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Somyajit K, Mishra A, Jameei A, Nagaraju G. Enhanced non-homologous end joining contributes toward synthetic lethality of pathological RAD51C mutants with poly (ADP-ribose) polymerase. Carcinogenesis. (2015) 36:13–24. doi: 10.1093/carcin/bgu211

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Perkhofer L, Schmitt A, Romero Carrasco MC, Ihle M, Hampp S, Ruess DA, et al. ATM deficiency generating genomic instability sensitizes pancreatic ductal adenocarcinoma cells to therapy-induced DNA damage. Cancer Res. (2017) 77:5576–90. doi: 10.1158/0008-5472.CAN-17-0634

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Aschele C, Cionini L, Lonardi S, Pinto C, Cordio S, Rosati G, et al. Primary tumor response to preoperative chemoradiation with or without oxaliplatin in locally advanced rectal cancer: pathologic results of the STAR-01 randomized phase III trial. J Clin Oncol. (2011) 29:2773–80. doi: 10.1200/JCO.2010.34.4911

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Hofheinz RD, Arnold D, Fokas E, Kaufmann M, Hothorn T, Folprecht G, et al. Impact of age on the efficacy of oxaliplatin in the preoperative chemoradiotherapy and adjuvant chemotherapy of rectal cancer: a post hoc analysis of the CAO/ARO/AIO-04 phase III trial. Ann Oncol. (2018) 29:1793–9. doi: 10.1093/annonc/mdy205

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Chen MB, Wu XY, Yu R, Li C, Wang LQ, Shen W, et al. P53 status as a predictive biomarker for patients receiving neoadjuvant radiation-based treatment: a meta-analysis in rectal cancer. PLoS ONE. (2012) 7:e45388. doi: 10.1371/journal.pone.0045388

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Duldulao MP, Lee W, Nelson RA, Li W, Chen Z, Kim J, et al. Mutations in specific codons of the KRAS oncogene are associated with variable resistance to neoadjuvant chemoradiation therapy in patients with rectal adenocarcinoma. Ann Surg Oncol. (2013) 20:2166–71. doi: 10.1245/s10434-013-2910-0

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Lee JW, Lee JH, Shim BY, Kim SH, Chung MJ, Kye BH, et al. KRAS mutation status is not a predictor for tumor response and survival in rectal cancer patients who received preoperative radiotherapy with 5-fluoropyrimidine followed by curative surgery. Medicine. (2015) 94:e1284. doi: 10.1097/MD.0000000000001284

CrossRef Full Text | Google Scholar

46. Ma L, Liu T, Jin Y, Wei J, Yang Y, Zhang H. ABCG2 is required for self-renewal and chemoresistance of CD133-positive human colorectal cancer cells. Tumour Biol. (2016) 37:12889–96. doi: 10.1007/s13277-016-5209-5

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Chen M, Li D, Gong N, Wu H, Su C, Xie C, et al. miR-133b down-regulates ABCC1 and enhances the sensitivity of CRC to anti-tumor drugs. Oncotarget. (2017) 8:52983–94. doi: 10.18632/oncotarget.17677

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Porro A, Haber M, Diolaiti D, Iraci N, Henderson M, Gherardi S, et al. Direct and coordinate regulation of ATP-binding cassette transporter genes by Myc factors generates specific transcription signatures that significantly affect the chemoresistance phenotype of cancer cells. J Biol Chem. (2010) 285:19532–43. doi: 10.1074/jbc.M109.078584

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Hardiman KM, Ulintz PJ, Kuick RD, Hovelson DH, Gates CM, Bhasi A, et al. Intra-tumor genetic heterogeneity in rectal cancer. Lab Invest. (2016) 96:4–15. doi: 10.1038/labinvest.2015.131

PubMed Abstract | CrossRef Full Text | Google Scholar