Copy number alterations and allelic ratio in relation to recurrence of rectal cancer
BMC Genomics volume 16, Article number: 438 (2015)
In rectal cancer, total mesorectal excision surgery combined with preoperative (chemo)radiotherapy reduces local recurrence rates but does not improve overall patient survival, a result that may be due to the harmful side effects and/or co-morbidity of preoperative treatment. New biomarkers are needed to facilitate identification of rectal cancer patients at high risk for local recurrent disease. This would allow for preoperative (chemo)radiotherapy to be restricted to high-risk patients, thereby reducing overtreatment and allowing personalized treatment protocols. We analyzed genome-wide DNA copy number (CN) and allelic alterations in 112 tumors from preoperatively untreated rectal cancer patients. Sixty-six patients with local and/or distant recurrent disease were compared to matched controls without recurrence. Results were validated in a second cohort of tumors from 95 matched rectal cancer patients. Additionally, we performed a meta-analysis that included 42 studies reporting on CN alterations in colorectal cancer and compared results to our own data.
The genomic profiles in our study were comparable to other rectal cancer studies. Results of the meta-analysis supported the hypothesis that colon cancer and rectal cancer may be distinct disease entities. In our discovery patient study cohort, allelic retention of chromosome 7 was significantly associated with local recurrent disease. Data from the validation cohort were supportive, albeit not statistically significant, of this finding.
We showed that retention of heterozygosity on chromosome 7 may be associated with local recurrence in rectal cancer. Further research is warranted to elucidate the mechanisms and effect of retention of chromosome 7 on the development of local recurrent disease in rectal cancer.
The Dutch Total Mesorectal Excision (TME) trial  changed standard treatment guidelines for rectal cancer patients . This international trial was designed to provide a clinical assessment of whether additional preoperative short-term radiotherapy (pRT) could reduce the number of local recurrences compared to TME surgery alone. Based on the Dutch TME trial, patients with tumor stage 2 (T2) were recommended to receive pRT in addition to TME surgery . The beneficial effect of TME surgery combined with pRT on rectal cancer local recurrence rates was subsequently confirmed in a Dutch population-based study . In addition to pRT, preoperative concurrent chemoradiation treatment (pCCRT) has been studied internationally and proposed as a preoperative treatment for the reduction of local recurrent disease related to rectal cancer. Reductions in local recurrence achieved by pCCRT were similar to those for pRT , and the postoperative complications associated with pRT and pCCRT were also reported to be comparable . In 2013, an EURECCA consensus was published which provided treatment recommendations for rectal cancer aimed at minimizing the differences in rectal cancer treatment regimes in Europe [6,7]. Although both pRT and pCCRT are recommended for specific TNM stages and introduction of these approaches has reduced local recurrence rates markedly, no difference in overall survival has yet been achieved [2,4,8]. This lack of survival benefit is most likely caused by the harmful side effects  and co-morbidity of the preoperative treatment . Recently, a phase III clinical trial named RAPIDO was initiated to assess the survival benefit of pRT followed by full-dose preoperative chemotherapy as an alternative for pCCRT with and without postoperative chemotherapy .
Under current guidelines, preoperative treatment is offered to a broad group of patients, although only around 10% of patients are actually at risk for development of a local recurrence without preoperative treatment . This means that up to 90% of rectal cancer patients are unnecessarily exposed to preoperative treatment . This undesirable situation calls for effective biomarkers that allow for selection of those patients with a high probability for development of a local recurrence after surgery alone. Effective selection for pRT or pCCRT treatment of only those rectal cancer patients with an increased risk of recurrent local disease would reduce overtreatment and allow for more personalized treatment protocols.
Genetic profiling of copy number (CN) alterations in individual chromosomes has previously been recognized as an independent predictor for metastatic relapse of early stage colorectal cancer . The risk of development of a local recurrence in rectal cancer can be deduced from the genetic aberrations found in the primary tumor. In the current study, we were especially interested in the prognostic value of genomic alterations related to local recurrent disease in rectal cancer. The aim of our study was therefore to identify CN alterations and allelic imbalances that could predict risk of recurrence in rectal cancer. Using a genome-wide CN analysis approach, we analyzed tissues from 112 rectal cancer patients enrolled in the TME trial. High-density SNP arrays were used to assess CN alterations and allelic (im)balances of genomic DNA (gDNA) segments [12,13]. As many studies have reported on CN alterations in (colo)rectal cancer, we additionally performed a meta-analysis of 42 studies that reported chromosomal CN alterations (high frequency CN alterations at the level of chromosome arms) and compared these to our current findings. We hypothesized that the CN alterations and/or aberrant allelic ratio patterns of certain gDNA segments might be prognostic for a local recurrence in rectal cancer, and might therefore identify those patients that would benefit most from pRT or pCCRT.
Tumor and matching normal tissues from the Dutch TME trial, which recruited rectal cancer patients in 118 European centers and one Canadian center, were available for analysis. The trial design has been described previously . We selected patients (from the non-radiation treatment arm) with clinically resected TNM stage I-III adenocarcinomas of the rectum and a presentation of local and/or distant recurrent disease at follow-up. To exclude uncontrolled bias, patients with clinically resected TNM stage I-III adenocarcinomas of the rectum presenting without recurrent disease were additionally selected to match those presenting with recurrent disease. Individual control group patients matched patients in the local recurrence group, the distant recurrence group and the local & distant recurrence group. For every patient within a recurrence group (local, distant, or local & distant) a unique match was therefore included. Matching criteria were TNM stage (exact match), CRM involvement (exact match), gender (exact match), and age at surgery (average difference of 2 years; range 0-7 years).
For the discovery phase, we selected tissues from 112 rectal cancer patients for whom fresh-frozen (FF) tumor and normal matched tissues were available. These included samples from patients with a recurrence (all available local N = 10, distant N = 41, and all available local & distant N = 15) and matched control patients without a recurrence (N = 46). Three samples from the discovery cohort were used for technical validation. For validation of results, we selected tumor and normal matched formalin-fixed paraffin-embedded (FFPE) tissue from 95 patients. Patients presenting with a recurrence (all available local, N = 12, distant, N = 24, and all available local & distant, N = 22) were matched to controls (without any recurrence, N = 37). The discovery study cohort (N = 112) and validation cohort (N = 95), which did not overlap, are described in Table 1.
All samples were coded in accordance with national ethical guidelines (“Code for Proper Secondary Use of Human Tissue”, Dutch Federation of Medical Scientific Societies). The use of these specimens was approved by the Medical Ethical Committee of the Leiden University Medical Center (LUMC). Patients were included in The Dutch Total Mesorectal Excision (TME) trial after informed consent was obtained .
Genomic DNA (gDNA) from tumor and corresponding normal tissue (FF or FFPE) was isolated, as described previously [14,15], following macrodissection or laser capture microdissection (LCM). The DNA concentration was measured using the PicoGreen method (Invitrogen, Carlsbad, CA, USA). If necessary, samples were concentrated using a Speedvac (SC110A, ThermoFisher, Waltham, MA, USA) to 10 ng/μl or to a minimum of 15 μl.
CytoSNP arrays and analysis
Hybridization of gDNA to high resolution Illumina Human CytoSNP12v2 arrays (Illumina, San Diego, CA, USA), intensity data extraction and the first quality control steps were performed by ServiceXS (Leiden, the Netherlands). Bioinformatic analysis was performed using the beadarraySNP package in R . After normalization, de-waving and automated segmentation analysis, we obtained CN data and allelic ratio data for all arrays. Thresholds for CN alterations after normalization were set at 0.92 for losses and 1.08 for gains, compared to the sample mean. To designate the allelic ratio groups, DNA segments were divided per individual sample into three classes, retention of heterozygosity, imbalance and loss of heterozygosity (LOH).
A globaltest (R package globaltest [16,17]) was performed on both the continuous CN data and the overall allelic ratio group data separately to determine overall statistical differences between analysis groups, comparing local recurrence group (L-group) versus control group (C-group), distant recurrence group (D-group) versus C-group and local & distant recurrence group (LD-group) versus C-group. A significant global test on the overall CN data or allelic ratio group data resulted in a second global test on individual chromosome arms. If significant, the smaller underlying segments underwent a further global test to identify regions with significant differences between the analysis groups. Multiple testing correction was performed using the Benjamin-Hochberg (BH) method .
For allelic ratio group analyses of chromosome combinations, the overall chromosome status was classified into balanced = 1, imbalanced = 2 and LOH = 3. Overall chromosome status was defined as most abundant allelic ratio group number on the chromosome. In the chromosome combination analyses, the highest number of overall chromosome status was used. Differences between groups were assessed with Fisher’s exact test for count data.
Survival analysis was performed using the Cox proportional hazards model (R-package survival). Survival data were available for 12 years of follow-up. The mean follow-up time of the discovery cohort and the validation cohort were 6.6 years (range 0.07-13.4 years) and 5.7 years (range 0.09-13.6 years), respectively. Overall survival (OS) was defined as the time from surgery until death by any cause. Disease specific survival (DSS) was defined as the time from surgery until death by rectal cancer. Local recurrence-free period (LRFP) was defined as the time from surgery until the discovery of a local recurrence. Distant recurrence-free period (DRFP) was defined as the time from surgery until the discovery of a distant recurrence. Multivariate models included the predetermined clinically important covariates TNM stage, age at surgery, gender and circumferential margin involvement, irrespective of statistical significance.
Dynamic array and analysis
Based on the results from the discovery phase, reference SNP (rs) identification numbers on chromosome 7 and 13 were extracted from the CytoSNP12 array. A search was performed for validated ABI-Taqman SNP assays, based on rs-numbers, in the SNP browser program (Applied Biosystems, Foster City, CA, USA) and 48 SNPs were selected for validation, with 16 on chromosome 7p, 16 on chromosome 7q and 16 on chromosome 13 (Additional file 1 S1). Selection was based on the highest percentage of heterozygosity (at least 40%) as determined in the genotype call for normal tissues (N = 112) in the discovery phase.
Taqman assays targeting the 48 selected SNPs were tested using the 96.96 BioMark Dynamic Array for quantitative Real-Time PCR (Fluidigm Corporation, San Francisco, CA, USA). On the 96.96 BioMark dynamic platform, Taqman SNP targets were assayed in duplicate and samples were assayed in triplicate (technical validation) or duplicate (validation) using 14 cycles of Specific Target Amplification for each sample replicate prior to the qPCR on the array. Assays were performed by ServiceXS. Fluidigm Real-Time PCR analysis software was used to extract cycle threshold (Ct) values, while Fluidigm SNP Genotyping analysis software was used to extract endpoint measurements at cycle 40 and for genotype calling. The linear derivative was used to determine Ct threshold values per SNP assay. Samples with a Ct >30 and missing values were set at Ct = 30. The number of starting DNA molecules per reaction chamber was computed using the following formula: # DNA molecules = 2^(26-Ct). In the 96.96 BioMark dynamic assay, a Ct value of 26 corresponds to one DNA molecule in the reaction chamber. Total sample molecule amount of the VIC (A-allele) and FAM (B-allele) channels was used for sample normalization. The allelic ratio, based on calculated amounts of DNA molecules rather than intensity, can be calculated by dividing the values for the VIC (A-allele) and FAM (B-allele) channels.
Validation of results
Based on quality control, samples with a mean log2 value lower than 4, and patients with non-matching tumor and normal SNP genotypes, were excluded from analyses. The Welch Two Sample T-test was used to assess whether average allelic ratio values along chromosome 7p, chromosome 7q or chromosome 13 were significantly different between the recurrence groups and the control group. Survival analysis for validation of results was performed as described above.
Results and discussion
Description of frequently occurring copy number alterations, and meta-analysis
We assessed CN alterations and chromosomal aberrations in rectal cancers from the Dutch TME trial, selected from the treatment arm that did not receive pRT. To gain a broader perspective on the results of our analysis of chromosomal CN alterations, we performed a meta-analysis of previously published CN studies. A PubMed search (performed on April 7, 2014), using the search criteria described in Additional file 2 S2, yielded 325 published studies. Studies from 2005 to search date and additional studies (selected from references of the selected articles published in 2005 or after) that described older CN alteration studies were assessed for eligibility, which required description of at least 6 tumors and detailed information on CN alterations per chromosome arm. A total of 42 studies that reported data on CN alterations in rectal cancer or colorectal cancer were included [11,19-59]. We combined results from these 42 studies (Table 2) and one additional study published online just after our Pubmed search (Table 2) and used the data to locate high frequency CN alterations in colorectal cancer at the level of chromosome arms (Figure 1A). CN alterations reported in at least 40% of the studies for at least 25% of the study cases were considered to be common CN alterations with high frequencies.
Many included studies pooled colon and rectal tumors and referred to them as colorectal cancer. However, evidence suggests that tumors arising from colon and rectal tissues should be considered as distinct disease entities [60-62]. We therefore compared common CN alterations with high frequencies in rectal versus colon cancers. Twenty studies (including the present study) specifically reported CN alterations in rectal carcinomas (Figure 1B). Frequent gains common in rectal cancer included regions of chromosomes 7p (found in 75% of the 20 rectal cancer studies, with a frequency >25%), 7q (65%), 8q (100%), 13q (100%), 20p (70%) and 20q (85%). Frequent losses common in rectal cancer were found on chromosomes 1p (found in 55% of the 20 rectal cancer studies, with a frequency >25%), 4q (50%), 5q (40%), 8p (80%), 14q (40%), 15q (65%), 17p (95%), 18p (70%), 18q (100%) and 22q (50%). Twelve studies reported CN alterations specific to colon carcinomas (Figure 1C), with frequent gains common in colon cancer on chromosomes 7p (found in 67% of the 12 studies, with a frequency >25%), 7q (75%), 8q (75%), 13q (83%), 20p (42%) and 20q (92%). Frequent losses common in colon cancer were found on chromosomes 8p (found in 42% of the 12 studies, with a frequency >25%), 17p (67%), 18p (67%) and 18q (92%). In contrast to rectal cancer, frequent losses on chromosomes 1p, 4q, 5q, 14q, 15q and 22q were not observed in at least 40% of the colon studies. With the exception of loss at 4q, these chromosome arms showed a difference of >20% between rectal cancer and colon cancer studies, indicating that rectal cancers were associated with a higher frequency of losses on these particular chromosome arms (Figure 1D). In addition, the percentage of studies reporting frequent losses on chromosomes 1q, 8p, 17p, and 21q, and the percentage of studies reporting frequent gains on chromosomes 8q and 20p, was higher in rectal cancer compared to colon cancer. These differences may be explained, in part, by the much higher percentage of microsatellite unstable colon carcinomas in comparison to rectal carcinomas, because microsatellite instability reflects DNA mismatch repair deficiency, in contrast to chromosomal instability characterized by aneuploidy. Where possible, microsatellite unstable tumors were omitted from our comparisons, although many studies did not report microsatellite stability status.
In conclusion, frequent CN alterations common in rectal cancer were located on chromosomes 1p, 4q, 5q, 7, 8, 13q, 14q, 15q, 17p, 18, 20 and 22q, while colon cancer showed a similar pattern but with a lower frequency of gains at chromosomes 8q and 20p, a lower frequency of losses at chromosomes 1q, 8p, 17p and 21q, and an absence of frequent CN losses at chromosomes 1p, 4q, 5q, 14q, 15q and 22q.
Comparison of our study with the meta-analysis
In order to identify CN alterations associated with a recurrence in rectal cancer, we assessed CN alterations in a series of non-preoperatively treated rectal cancers from the Dutch TME trial. Rectal cancers with local and/or distant recurrences were compared with matching control cases without recurrence. In terms of the location of frequent gains, our study cohort was comparable to previously described rectal cancer cohorts (Figure 1B upper panel). In contrast, compared to previous studies we detected a higher number of chromosome arms with frequent losses (Figure 1B lower panel). This might be due to our use of high-density arrays rather than CGH arrays, and to the selection of previously studied CGH locations. Another explanation might be the larger sample size of our study compared to most of the earlier studies, as smaller studies are more prone to higher variance. Additionally, our study differs from other studies due to our use of both laser capture microdissection and macrodissection. This could have influenced results by reducing a potential effect in the analyses.
Genomic abnormalities and local recurrence
Continuous CN values and allelic ratio groups were analyzed in an effort to identify differentially affected genomic regions in rectal cancer patients presenting with and without local and/or distant recurrence. We included all patients with a local recurrence for whom fresh frozen tissue was available (N = 25). Fifteen of these patients also presented with a distant recurrence. Additionally, we included 41 patients presenting with only distant recurrences. All patients with recurrent disease were matched to control samples of patients who did not present with a recurrence during follow-up (N = 46).
The frequency of gains (Additional file 3 S3A) and losses (Additional file 3 S3B) along the length of each chromosome was plotted for each analysis group (control group, local recurrence group, distant recurrence group, and local & distant recurrence group). CN gains were less frequently observed in the local recurrence group compared to the other analysis groups, while CN losses were more frequently observed in patients with a local recurrence. However, CN alterations were not significantly associated with (local) recurrent disease in our study cohort. A similar pattern can be deduced from a study by Diep et al. , with fewer CN gains and more CN losses in local recurrences compared to independent primary colorectal tumor samples. A lower frequency of gains along chromosome 7 and higher frequency of losses along chromosome 18q were especially associated with local recurrence in both studies. A study by Kodeda et al. , comparing locally recurrent and non-recurrent rectal tumors, reported that both CN gains and CN losses were less frequently observed in locally recurrent tumors. Lower percentages of CN gains on chromosomes 8q, 13q and 20q, and a lower frequency of CN losses on chromosome 18q were observed in locally recurrent tumors in both our study cohort and the Kodeda cohort. In addition, we observed a lower frequency of CN gains on chromosomes 9p, 16p, 17q, 19 and 22q in the locally recurrent tumors. In contrast to the Kodeda study, we observed a higher frequency of CN losses on chromosomes 1p and 5q in locally recurrent tumors, in addition to a higher frequency of losses on chromosomes 3p, 9p and 15q. Overall, the results of the Kodeda study were comparable with our study regarding the lower frequency of CN gains, but conflicting on frequency of CN losses. This difference might be a result of the higher number of small deletions identified in our study due to our use of high-density SNP arrays compared to the CGH arrays used by Kodeda et al. However, despite these differences both the Kodeda study  and our study share a common conclusion that none of the above described CN alterations are (statistically) significantly associated with local recurrent disease.
An important feature of the SNP arrays used in our study was that they allow assessment of allelic imbalances, in addition to CN alterations. Using array data, allelic ratio groups were analyzed in order to identify genomic regions that showed allelic (im)balances associated with local and/or distant recurrence. The allelic ratio of DNA segments was classified into three classes: retention of heterozygosity, imbalance and LOH. In contrast to CN alterations, allelic imbalance showed a significant association with local recurrence, while no association was found with distant recurrence or a combination of both. Tumors from patients in the local recurrence group showed overall statistically significant differences in allelic ratios compared to the control samples (p < 0.005). Individual chromosome arms, and underlying segments and sub-segments, were then analyzed. Several chromosome arms showed different allelic ratio groups with (uncorrected) p-values <0.05 (Table 3A). In all cases, the local recurrence group displayed fewer allelic aberrations (allelic imbalance or LOH) than the control group. Statistically significant sub-regions, with p-values <0.05 after adjustment for multiple testing, were identified on chromosome 7. These regions showed almost no LOH (Additional file 4 S4). All tumors from patients presenting with only a local recurrence showed retention of heterozygosity and balanced alleles. In contrast, approximately 50 percent of control group cases showed imbalanced alleles. The percentage of patients with retention, imbalanced alleles or LOH along the length of chromosome 7 is shown in Figure 2. Balanced alleles were more prominent in cases with only a local recurrence, indicating that retention on chromosome 7 might be associated with the development of local recurrent disease. Frequent gains of chromosome 7, with or without allelic imbalances, were most often identified, and true LOH was almost non-existent in our study cohort. This suggests that heterozygosity of chromosome 7 is important for rectal cancer tumorigenesis and warrants further studies on the role of chromosome 7 retention in local recurrences of rectal cancer.
The use of heterozygous SNPs to determine allelic ratios of normal and tumor samples identified statistically significant differences on chromosome 7 between patients presenting with and without only local recurrent disease.
A validation of the association of chromosome 7 retention with local recurrence was performed using the 96.96 BioMark dynamic platform. A technical validation demonstrated that this platform could be used for validation of our results (Additional file 5 S5). We selected all (remaining) patients with a local recurrence for whom FFPE tissue was available (N = 34). Of these patients, 22 also presented with a distant recurrence. Additionally, we included 24 patients presenting with only distant recurrences. All patients with recurrent disease were matched to control patient samples without recurrence in follow-up (N = 37).
Differences in frequency of allelic (im)balances of specific SNPs on chromosome 7 were assessed in FFPE tissues of rectal cancer patients presenting with and without local and/or distant recurrence. Based on the results from the discovery phase, we selected 32 SNPs located on chromosome 7p (N = 16) and chromosome 7q (N = 16) (Additional file 1 S1). As a high percentage of patients presenting with only a local recurrence also showed retention on chromosome 13 (Figure 2), an additional set of 16 SNPs was selected on this chromosome to potentially increase discriminative power. SNPs on both chromosome 13 and chromosome 7 were therefore included to enhance differences between analysis groups. During the discovery phase, the overall chromosome status of combined chromosomes 7 and 13 was significantly different in patients in the local recurrence group compared to the control group (Fisher’s exact test for count data, p-value < 0.0002; Table 3B). Balanced alleles were more prominent in cases with a local recurrence, indicating that retention on chromosome 7, and to a lesser extent on chromosome 13, might be associated with the development of a local recurrence.
Using the technically validated dynamic array approach, differences between the local recurrence group and the control group showed a trend towards significance on chromosome 7p (p-value = 0.07). Retention of heterozygosity was more frequently observed in tumors with local recurrence compared to control group tumors, which is in accordance with findings in the discovery phase (Figure 3). The telomeric region of chromosome 7p showed allelic imbalance or LOH in consistently lower percentages of patients in the local recurrence group compared to the control group. For the centromeric region, the same pattern was observed for both the local recurrence group and for the local & distant recurrence group. Chromosome 13 retention results could not be reproduced.
In brief, validation cohort data were not conclusive but do support the notion that retention of heterozygosity on the telomeric and centromeric regions of chromosome 7 may be associated with local recurrence in rectal cancer.
Prognostic value of CN alterations and allelic aberrations
The clinical prognostic value of continuous CN profiles and allelic ratio group profiles was assessed in the discovery cohort using the Cox proportional hazards model. Multivariate analyses showed a trend in association of chromosome 7p with local recurrence (LRFP), but no associations with any other clinical outcome in our cohort of rectal cancer patients (Additional file 6 S6). In our discovery study cohort, allelic imbalance was associated with death by rectal cancer and local recurrent disease, with the allelic imbalance at chromosome 7p being chiefly responsible for the prognostic effect for local recurrent disease (Additional file 6 S6). These results could not be confirmed with the 48 selected SNPs using the dynamic array (data not shown), and this outcome is reflected in the absence of literature describing associations between CN alterations or allelic ratio profiles and patient survival in rectal cancer patients.
Only one study specifically dedicated to rectal cancer patients reported an association of CN alterations with patient survival or recurrent disease . The study by Doyen et al.  showed (in multivariate analysis) an association of loss at chromosome 8p with worse cancer specific survival (CSS) and the occurrence of metachronous distant metastases (DRFP). Unfortunately, the covariates included in the multivariate analysis where not reported. Four previously published colorectal cancer studies reported associations of CN alterations with patient survival or recurrent disease [11,45,48,64]. A study by De Angelis et al.  showed worse overall survival for patients with losses at chromosomes 1p and 8p in multivariate analyses. Bardi et al.  observed a significant association between loss at chromosome 4 and worse disease-free survival in univariate analyses. In multivariate analysis, loss of chromosome 18 was reported to be associated with worse overall survival. Al-Mulla et al.  reported that loss at chromosomes 4p and 5q was associated with worse disease-free survival in early stage colorectal cancers in multivariate analyses. Allelic imbalance at chromosome 8p was reported by Halling et al.  and associated with overall survival (OS) and time to recurrence (DRFP) in multivariate analysis. Association between allelic imbalance or loss at chromosome 8p and worse clinical outcome was described in three independent articles. However, only one was dedicated specifically to rectal cancer, and our study could not validate the results. The differences in results between these studies and our own data might be due to differences in tumor locations, as we focused on rectal cancer patients alone, whereas four out of five earlier studies included both colon and rectal cancer patients. Additionally, our use of laser capture microdissection along with macrodissection in order to reduce the potential effect of tumor microenvironment on the analyses, and thereby reduction of intratumoral heterogeneity of the rectal tumors, might provide a truer picture of the association between CN alterations and clinical outcome.
Allelic aberrations and LOH in colorectal cancer have been widely investigated, and loss of 18q and 17p are seen as prognostic markers for clinical outcomes in colorectal cancer (reviewed in ). To the best of our knowledge, no study has previously focused specifically on rectal cancers. A study by Choi et al.  showed that higher levels of LOH were significantly associated with tumor location specifically in the rec and the distal portion of the colon. This finding indicates that rectal and colon tumors should be considered distinct disease entities in relation to allelic alterations and LOH in particular. Prognostic indicators identified in studies of colorectal cancer patients cannot be adequately compared with (our) survival data on rectal cancer patients.
Implications for clinical use
Rectal cancer is associated with a high rate of local recurrence, in contrast to colon cancer. The location of the rectum, fixed in the smaller pelvis, provides opportunities for pRT or pCCRT treatment. Although introduction of pRT and pCCRT led to markedly reduced local recurrence rates (from 11% to 5% for pRT ), no difference in overall survival was observed [2,4,8]. This suggests that the majority of rectal cancer patients (over 90%) are currently receiving unnecessary preoperative treatment to reduce the local recurrence risk, when in fact this risk is only relevant for less than 10% of all patients. Identification of patients who are likely to show local recurrent disease could guide decision-making for the preoperative treatment of rectal cancer patients.
Validation data from the present study on allelic ratios of chromosome 7 are supportive (but not conclusive) of our initial finding that retention of heterozygosity on chromosome 7 is associated with local recurrent disease. While these data do not yet provide sufficient grounds for development of a clinically useful platform based on these observations, further research is warranted to elucidate underlying mechanisms. Chromosome 7 harbors many interesting genes, including druggable targets such as the oncogenes EGFR, BRAF and MET, but at present their relation to the retention of chromosome 7 and the development of local recurrent disease in rectal cancer is unclear. Comparison of data on CN alterations and allelic imbalance from our rectal cancer cohort with previously published studies provided support for the hypothesis that colon cancer and rectal cancer may be distinct disease entities, and thus may require stratification in (survival) analyses accordingly.
Kapiteijn E, Marijnen CA, Nagtegaal ID, Putter H, Steup WH, Wiggers T, et al. Preoperative radiotherapy combined with total mesorectal excision for resectable rectal cancer. N Engl J Med. 2001;345:638–46.
van Gijn W, Marijnen CA, Nagtegaal ID, Kranenbarg EM, Putter H, Wiggers T, et al. Preoperative radiotherapy combined with total mesorectal excision for resectable rectal cancer: 12-year follow-up of the multicentre, randomised controlled TME trial. Lancet Oncol. 2011;12:575–82.
den Dulk M, Krijnen P, Marijnen CA, Rutten HJ, van de Poll-Franse LV, Putter H, et al. Improved overall survival for patients with rectal cancer since 1990: the effects of TME surgery and pre-operative radiotherapy. Eur J Cancer. 2008;44:1710–6.
Yeh CH, Chen MF, Lai CH, Huang WS, Lee SP, Chen WC, et al. Comparison of treatment results between surgery alone, preoperative short-course radiotherapy, or long-course concurrent chemoradiotherapy in locally advanced rectal cancer. Int J Clin Oncol. 2012;17:482–90.
Bujko K, Nowacki MP, Kepka L, Oledzki J, Bebenek M, Kryj M, et al. Postoperative complications in patients irradiated pre-operatively for rectal cancer: report of a randomised trial comparing short-term radiotherapy vs chemoradiation. Colorectal Dis. 2005;7:410–6.
Valentini V, Glimelius B, Haustermans K, Marijnen CA, Rodel C, Gambacorta MA, et al. EURECCA consensus conference highlights about rectal cancer clinical management: the radiation oncologist’s expert review. Radiother Oncol. 2014;110:195–8.
van de Velde CJ, Aristei C, Boelens PG, Beets-Tan RG, Blomqvist L, Borras JM, et al. EURECCA colorectal: multidisciplinary mission statement on better care for patients with colon and rectal cancer in Europe. Eur J Cancer. 2013;49:2784–90.
Peeters KC, Marijnen CA, Nagtegaal ID, Kranenbarg EK, Putter H, Wiggers T, et al. The TME trial after a median follow-up of 6 years: increased local control but no survival benefit in irradiated patients with resectable rectal carcinoma. Ann Surg. 2007;246:693–701.
Marijnen CA, van de Velde CJ, Putter H, van den Brink M, Maas CP, Martijn H, et al. Impact of short-term preoperative radiotherapy on health-related quality of life and sexual functioning in primary rectal cancer: report of a multicenter randomized trial. J Clin Oncol. 2005;23:1847–58.
Nilsson PJ, van EB, Hospers GA, Pahlman L, van de Velde CJ, Beets-Tan RG, et al. Short-course radiotherapy followed by neo-adjuvant chemotherapy in locally advanced rectal cancer–the RAPIDO trial. BMC Cancer. 2013;13:279.
Al-Mulla F, Behbehani AI, Bitar MS, Varadharaj G, Going JJ. Genetic profiling of stage I and II colorectal cancer may predict metastatic relapse. Mod Pathol. 2006;19:648–58.
Corver WE, Middeldorp A, ter Haar NT, Jordanova ES, van PM, van ER, et al. Genome-wide allelic state analysis on flow-sorted tumor fractions provides an accurate measure of chromosomal aberrations. Cancer Res. 2008;68:10333–40.
Oosting J, Lips EH, van ER, Eilers PH, Szuhai K, Wijmenga C, et al. High-resolution copy number analysis of paraffin-embedded archival tissue using SNP BeadArrays. Genome Res. 2007;17:368–76.
DNA isolation protocol. http://med.stanford.edu/labs/vanderijn-west/documents/DNAextractionfromTRIZOL_Organicphase_update.doc. 2014.
van Eijk R, Stevens L, Morreau H, van WT. Assessment of a fully automated high-throughput DNA extraction method from formalin-fixed, paraffin-embedded tissue for KRAS, and BRAF somatic mutation analysis. Exp Mol Pathol. 2013;94:121–5.
Goeman JJ, van de Geer SA, de KF, van Houwelingen HC. A global test for groups of genes: testing association with a clinical outcome. Bioinformatics. 2004;20:93–9.
Goeman JJ, Oosting J, Cleton-Jansen AM, Anninga JK, van Houwelingen HC. Testing association of a pathway with survival using gene expression data. Bioinformatics. 2005;21:1950–7.
Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc, Series B. 1995;57:289–300.
Alcock HE, Stephenson TJ, Royds JA, Hammond DW. Analysis of colorectal tumor progression by microdissection and comparative genomic hybridization. Genes Chromosomes Cancer. 2003;37:369–80.
Ghadimi BM, Grade M, Liersch T, Langer C, Siemer A, Fuzesi L, et al. Gain of chromosome 8q23-24 is a predictive marker for lymph node positivity in colorectal cancer. Clin Cancer Res. 2003;9:1808–14.
He QJ, Zeng WF, Sham JS, Xie D, Yang XW, Lin HL, et al. Recurrent genetic alterations in 26 colorectal carcinomas and 21 adenomas from Chinese patients. Cancer Genet Cytogenet. 2003;144:112–8.
Meijer GA, Hermsen MA, Baak JP, van Diest PJ, Meuwissen SG, Belien JA, et al. Progression from colorectal adenoma to carcinoma is associated with non-random chromosomal gains as detected by comparative genomic hybridisation. J Clin Pathol. 1998;51:901–9.
Paredes-Zaglul A, Kang JJ, Essig YP, Mao W, Irby R, Wloch M, et al. Analysis of colorectal cancer by comparative genomic hybridization: evidence for induction of the metastatic phenotype by loss of tumor suppressor genes. Clin Cancer Res. 1998;4:879–86.
Ried T, Knutzen R, Steinbeck R, Blegen H, Schrock E, Heselmeyer K, et al. Comparative genomic hybridization reveals a specific pattern of chromosomal gains and losses during the genesis of colorectal tumors. Genes Chromosomes Cancer. 1996;15:234–45.
Tanami H, Tsuda H, Okabe S, Iwai T, Sugihara K, Imoto I, et al. Involvement of cyclin D3 in liver metastasis of colorectal cancer, revealed by genome-wide copy-number analysis. Lab Invest. 2005;85:1118–29.
Xiao XY, Zhou XY, Yan G, Sun MH, Du X. Chromosomal alteration in Chinese sporadic colorectal carcinomas detected by comparative genomic hybridization. Diagn Mol Pathol. 2007;16:96–103.
Chen Z, Liu Z, Li W, Qu K, Deng K, 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.
Chen Z, Liu Z, Deng X, Warden C, Li W, Garcia-Aguilar J. Chromosomal copy number alterations are associated with persistent lymph node metastasis after chemoradiation in locally advanced rectal cancer. Dis. Colon Rectum. 2012;55:677–85.
Grade M, Becker H, Liersch T, Ried T, Ghadimi BM. Molecular cytogenetics: genomic imbalances in colorectal cancer and their clinical impact. Cell Oncol. 2006;28:71–84.
Grade M, Gaedcke J, Wangsa D, Varma S, Beckmann J, Liersch T. Chromosomal copy number changes of locally advanced rectal cancers treated with preoperative chemoradiotherapy. Cancer Genet Cytogenet. 2009;193:19–28.
Kodeda K, Asting AG, Lonnroth C, Derwinger K, Wettergren Y, Nordgren S, et al. Genomic CGH-assessed structural DNA alterations in rectal carcinoma as related to local recurrence following primary operation for cure. Int J Oncol. 2012;41:1397–404.
Liang JW, Shi ZZ, Zhang TT, Hao JJ, Wang Z, Wang XM, et al. Analysis of genomic aberrations associated with the clinicopathological parameters of rectal cancer by arraybased comparative genomic hybridization. Oncol Rep. 2013;29:1827–34.
Lips EH, de Graaf EJ, Tollenaar RA, van ER, Oosting J, Szuhai K, et al. Single nucleotide polymorphism array analysis of chromosomal instability patterns discriminates rectal adenomas from carcinomas. J Pathol. 2007;212:269–77.
Lips EH, van ER, de Graaf EJ, Doornebosch PG, de Miranda NF, Oosting J, et al. Progression and tumor heterogeneity analysis in early rectal cancer. Clin Cancer Res. 2008;14:772–81.
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.
Shi ZZ, Zhang YM, Shang L, Hao JJ, Zhang TT, Wang BS, et al. Genomic profiling of rectal adenoma and carcinoma by array-based comparative genomic hybridization. BMC Med Genomics. 2012;5:52.
Zhou HT, Shi ZZ, Zhou ZX, Jiang YY, Hao JJ, Zhang TT, Shi F, et al. Genomic changes in rectal adenocarcinoma associated with liver metastasis. Cancer Biomark. 2013;13:281–8.
Grade M, Hormann P, Becker S, Hummon AB, Wangsa D, Varma S, et al. Gene expression profiling reveals a massive, aneuploidy-dependent transcriptional deregulation and distinct differences between lymph node-negative and lymph node-positive colon carcinomas. Cancer Res. 2007;67:41–56.
Korn WM, Yasutake T, Kuo WL, Warren RS, Collins C, Tomita M, et al. Chromosome arm 20q gains and other genomic alterations in colorectal cancer metastatic to liver, as analyzed by comparative genomic hybridization and fluorescence in situ hybridization. Genes Chromosomes Cancer. 1999;25:82–90.
Nakao K, Shibusawa M, Tsunoda A, Yoshizawa H, Murakami M, Kusano M, et al. Genetic changes in primary colorectal cancer by comparative genomic hybridization. Surg Today. 1998;28:567–9.
Schlegel J, Stumm G, Scherthan H, Bocker T, Zirngibl H, Ruschoff J, et al. Comparative genomic in situ hybridization of colon carcinomas with replication error. Cancer Res. 1995;55:6002–5.
Al-Mulla F, Keith WN, Pickford IR, Going JJ, Birnie GD. Comparative genomic hybridization analysis of primary colorectal carcinomas and their synchronous metastases. Genes Chromosomes Cancer. 1999;24:306–14.
Aragane H, Sakakura C, Nakanishi M, Yasuoka R, Fujita Y, Taniguchi H, et al. Chromosomal aberrations in colorectal cancers and liver metastases analyzed by comparative genomic hybridization. Int J Cancer. 2001;94:623–9.
Aust DE, Willenbucher RF, Terdiman JP, Ferrell LD, Chang CG, Moore DH, et al. Chromosomal alterations in ulcerative colitis-related and sporadic colorectal cancers by comparative genomic hybridization. Hum Pathol. 2000;31:109–14.
Bardi G, Fenger C, Johansson B, Mitelman F, Heim S. Tumor karyotype predicts clinical outcome in colorectal cancer patients. J Clin Oncol. 2004;22:2623–34.
Chan TL, Curtis LC, Leung SY, Farrington SM, Ho JW, Chan AS, et al. Early-onset colorectal cancer with stable microsatellite DNA and near-diploid chromosomes. Oncogene. 2001;20:4871–6.
De Angelis PM, Clausen OP, Schjolberg A, Stokke T. Chromosomal gains and losses in primary colorectal carcinomas detected by CGH and their associations with tumour DNA ploidy, genotypes and phenotypes. Br J Cancer. 1999;80:526–35.
De Angelis PM, Stokke T, Beigi M, Mjaland O, Clausen OP. Prognostic significance of recurrent chromosomal aberrations detected by comparative genomic hybridization in sporadic colorectal cancer. Int J Colorectal Dis. 2001;16:38–45.
Diep CB, Teixeira MR, Thorstensen L, Wiig JN, Eknaes M, Nesland JM, et al. Genome characteristics of primary carcinomas, local recurrences, carcinomatoses, and liver metastases from colorectal cancer patients. Mol Cancer. 2004;3:6.
Georgiades IB, Curtis LJ, Morris RM, Bird CC, Wyllie AH. Heterogeneity studies identify a subset of sporadic colorectal cancers without evidence for chromosomal or microsatellite instability. Oncogene. 1999;18:7933–40.
Hermsen M, Postma C, Baak J, Weiss M, Rapallo A, Sciutto A, et al. Colorectal adenoma to carcinoma progression follows multiple pathways of chromosomal instability. Gastroenterology. 2002;123:1109–19.
Knosel T, Petersen S, Schwabe H, Schluns K, Stein U, Schlag PM, et al. Incidence of chromosomal imbalances in advanced colorectal carcinomas and their metastases. Virchows Arch. 2002;440:187–94.
Lagerstedt KK, Kristiansson E, Lonnroth C, Andersson M, Iresjo BM, Gustafsson A, et al. Genes with relevance for early to late progression of colon carcinoma based on combined genomic and transcriptomic information from the same patients. Cancer Inform. 2010;9:79–91.
Nakao K, Shibusawa M, Ishihara A, Yoshizawa H, Tsunoda A, Kusano M, et al. Genetic changes in colorectal carcinoma tumors with liver metastases analyzed by comparative genomic hybridization and DNA ploidy. Cancer. 2001;91:721–6.
Nakao M, Kawauchi S, Uchiyama T, Adachi J, Ito H, Chochi Y, et al. DNA copy number aberrations associated with the clinicopathological features of colorectal cancers: Identification of genomic biomarkers by array-based comparative genomic hybridization. Oncol Rep. 2011;25:1603–11.
Poeaim S, Rerkamnuaychoke B, Jesdapatarakul S, Campiranon A. Chromosome alterations in colorectal cancer in Thai patients. Cancer Genet Cytogenet. 2005;160:152–9.
Rooney PH, Boonsong A, McKay JA, Marsh S, Stevenson DA, Murray GI, et al. Colorectal cancer genomics: evidence for multiple genotypes which influence survival. Br J Cancer. 2001;85:1492–8.
Leslie A, Pratt NR, Gillespie K, Sales M, Kernohan NM, Smith G, et al. Mutations of APC, K-ras, and p53 are associated with specific chromosomal aberrations in colorectal adenocarcinomas. Cancer Res. 2003;63:4656–61.
Nakao K, Mehta KR, Fridlyand J, Moore DH, Jain AN, Lafuente A, et al. High-resolution analysis of DNA copy number alterations in colorectal cancer by array-based comparative genomic hybridization. Carcinogenesis. 2004;25:1345–57.
Birkenkamp-Demtroder K, Olesen SH, Sorensen FB, Laurberg S, Laiho P, Aaltonen LA, et al. Differential gene expression in colon cancer of the caecum versus the sigmoid and rectosigmoid. Gut. 2005;54:374–84.
Komuro K, Tada M, Tamoto E, Kawakami A, Matsunaga A, Teramoto K, et al. Right- and left-sided colorectal cancers display distinct expression profiles and the anatomical stratification allows a high accuracy prediction of lymph node metastasis. J Surg Res. 2005;124:216–24.
Goossens-Beumer IJ, Zeestraten EC, Benard A, Christen T, Reimers MS, Keijzer R, et al. Clinical prognostic value of combined analysis of Aldh1, Survivin, and EpCAM expression in colorectal cancer. Br J Cancer. 2014;110:2935–44.
Doyen J, Letouze E, Marisa L, de RA, Milano G, Etienne-Grimaldi MC, et al. High-resolution analysis of DNA copy number alterations in rectal cancer: correlation with metastasis, survival, and mRNA expression. Strahlenther Onkol. 2014;190:1028–36.
Halling KC, French AJ, McDonnell SK, Burgart LJ, Schaid DJ, Peterson BJ, et al. Microsatellite instability and 8p allelic imbalance in stage B2 and C colorectal cancers. J Natl Cancer Inst. 1999;91:1295–303.
Westra JL, Plukker JT, Buys CH, Hofstra RM. Genetic alterations in locally advanced stage II/III colon cancer: a search for prognostic markers. Clin Colorectal Cancer. 2004;4:252–9.
Choi SW, Lee KJ, Bae YA, Min KO, Kwon MS, Kim KM, et al. Genetic classification of colorectal cancer based on chromosomal loss and microsatellite instability predicts survival. Clin Cancer Res. 2002;8:2311–22.
This study was performed within the framework of CTMM, the Center for Translational Molecular Medicine. DeCoDe project (grant 03O-101). We thank Ronald van Vlierberghe, Rob Keyzer, Melanie Schrumpf, Fahim Behrouz and Matthijs Kuppen for expert technical support.
The authors declare that they have no competing interests.
IJGB drafted the manuscript, participated in all aspects of the study, study design, analysed all data and performed the meta-analysis. JO analysed all data and carried out the statistical analysis. WEC participated in analysis and interpretation of CytoSNP array data. MJFWJ carried out and participated in the design and analysis of the Dynamic array experiments. BJ participated in the design of the Dynamic array experiments and in the coordination of CytoSNP array data acquisition. WVW participated in the design of the Dynamic array experiments and coordination of CytoSNP array studies. ECMZ participated in sample acquisition and sample analysis. CJHVDV participated in the collection of clinical data. HM co-conceived the study, participated in the study design and sample acquisition, and maintained oversight of all aspects of the study. PJKK co-conceived the study, participated in sample acquisition and study coordination and provided supervision on all aspects of the study. TVW co-conceived and designed the study, drafted the manuscript and participated in all aspects of the study. All authors critically read, corrected and approved the final manuscript.
Search criteria for meta-analysis.
Frequency plots for CN gains.
Statistically significant regions in allelic ratio on chromosome 7q.
Technical validation of SNPs.
Overall clinical prognostic value for copy number and allelic ratio.
About this article
Cite this article
Goossens-Beumer, I.J., Oosting, J., Corver, W.E. et al. Copy number alterations and allelic ratio in relation to recurrence of rectal cancer. BMC Genomics 16, 438 (2015). https://doi.org/10.1186/s12864-015-1550-0
- Rectal Cancer
- Total Mesorectal Excision
- Rectal Cancer Patient
- Copy Number Alteration
- Allelic Imbalance