Purpose: Preoperative chemoradiotherapy (CRT) is the treatment of choice for rectal cancer (RC), but half of the patients do not respond, suffer unnecessary toxicities, and surgery delays. We aimed to develop a model that could predict a clinically meaningful response to CRT by using formalin-fixed paraffin-embedded (FFPE) biopsies.

Experimental Design: We first carried out an exploratory screening of candidate genes by using SAGE technology to evaluate dynamic changes in the RC transcriptome in selected refractory patients before and after CRT. Next, 53 genes (24 from SAGE and 29 from the literature) were analyzed by qPCR arrays in FFPE initial biopsies from 94 stage II/III RC patients who were preoperatively treated with CRT. Tumor response was defined by using Dworak's tumor regression grade (2–3–4 vs. 0–1). Multivariate Cox methods and stepwise algorithms were applied to generate an optimized predictor of response and outcome.

Results: In the training cohort (57 patients), a 13-gene signature predicted tumor response with 86% accuracy, 87% sensitivity, and 82% specificity. In a testing cohort (37 patients), the model correctly classified 6 of 7 nonresponders, with an overall accuracy of 76%. A signature-based score identified patients with a higher risk of relapse in univariate (3-year disease-free survival 64% vs. 90%, P = 0.001) and multivariate analysis (HR = 4.35 95% CI: 1.2–15.75, P = 0.02), in which it remained the only statistically significant prognostic factor.

Conclusions: A basal 13-gene signature efficiently predicted CRT response and outcome. Multicentric validation by the GEMCAD collaborative group is currently ongoing. If confirmed, the predictor could be used to improve patient selection in RC studies. Clin Cancer Res; 17(12); 4145–54. ©2011 AACR.

Translational Relevance

One of the most pressing questions in rectal cancer (RC) to date is to accurately identify those patients that benefit from preoperative chemoradiotherapy (CRT). Unfortunately, no basal prognostic factors predict response to preoperative CRT. To our knowledge, only 4 gene signatures have been published in RC aiming to predict pathologic response to CRT. None have been shown to predict patients with poor pathologic response and bad prognosis, who are likely to be the best candidates to explore alternative strategies to standard CRT. Herein, we present a basal 13-gene signature, the first one in RC obtained from clinical routinely prepared paraffin embedded tumors by using quantitative PCR that efficiently classifies CRT-induced tumor regression grade (accuracy 86%) and predicts poor responders in training and testing sets of patients. Importantly, a signature-based score remained the only prognostic factor in the multivariate analysis.

Preoperative chemoradiotherapy (CRT) is the standard treatment for stage II/III rectal cancer (RC). Compared with postoperative therapy, it improved local control and tolerability in the CAO/AROAIO 94 trial (1). More recently, the NSABP R-03 trial has reported a significant benefit in 5-year disease-free survival (DFS; ref. 2). Unfortunately, 5-fluorouracil (5FU)-based CRT produces a complete pathologic response in only 8% to 14% of the patients, and distant metastases develop in one-third of the cases after 5 years (1, 3, 4). Furthermore, patients who do not respond to preoperative CRT are exposed to unnecessary toxicities and resection of their primary tumor is delayed. Thus, the ability to predict treatment response and outcome before CRT is of the utmost importance to clinicians.

A number of postsurgical prognostic factors have been proposed (5–7), but CRT response cannot be clinically predicted. The identification of basal markers of resistance biomarkers could offer great help on this regard. Directed strategies that explore individual markers have not yielded clinically validated assays (7–9). Instead, gene expression arrays are appealing to predict the multifactorial event of treatment response (10–13). TaqMan Low Density Arrays (TLDA) can be done from routine formalin-fixed paraffin-embedded (FFPE) biopsies, making this platform of great interest for clinical application, as shown for the Oncotype Dx Assay (Genomic Health Inc.) signature in breast cancer (14), endorsed by ASCO (15) and NCCN guidelines (16). Selecting appropriate genes for analysis seems critical.

In this study, we sought to define a clinically exploitable gene signature that could predict treatment response and outcome in preoperatively treated RC patients. First, we carried out an undirected screening of candidate genes by evaluating dynamic changes in the RC transcriptome, before and after CRT in selected treatment-refractory patients. Serial analysis of gene expression (SAGE) was ideal because it does not require a priori knowledge of the transcripts that are present in cells and it accurately quantitates mRNA levels of the transcriptome, allowing the discovery of new molecular targets (17). It provided a set of 24 genes whose expression was altered after CRT. In the second stage, we evaluated 29 genes selected from the literature and the SAGE-emerging candidates (for a total of 53 genes) by using TLDA cards in FFPE initial biopsies from 94 locally advanced preoperatively treated RC patients. Basal changes in gene expression were correlated with tumor regression grade (TRG) and patient outcome (defined as 3-year disease free survival). A set of 13 genes was identified that could predict both endpoints.

The experimental design is depicted in Figure 1.

Figure 1.

Flow chart of patients. Abbreviations: HCB, Hospital Clinic Barcelona; HULP, Hospital Universitario La Paz.

Figure 1.

Flow chart of patients. Abbreviations: HCB, Hospital Clinic Barcelona; HULP, Hospital Universitario La Paz.

Close modal

Patients and treatment

Undirected screening of candidate genes.

For the initial SAGE experiment, 25 patients from a RC neoadjuvant CRT trial were prospectively recruited (18). Chemotherapy consisted of a preoperative combination of oxaliplatin and raltitrexed (130 mg/m2 and 3 mg/m2, days 1, 21, and 42) on 3 cycles, together with preoperative radiotherapy (RT; 50.4 Gy in 28 fractions). We sought to evaluate dynamic changes in tumor transcriptome after CRT in refractory patients, as defined by lack of downstaging and minimal or null pathologic response according to Dworak's classification of TRG (TRG 0–1; ref. 19), with similar initial clinical stage (cT3N0, as diagnosed by using endorectal ultrasounds and MRI) and adequate fresh tissue samples. Successful libraries were obtained in 3 cases (patients A, B, and C) fulfilling the criteria. Patients A and C developed lung metastases 16 and 54 months after diagnosis, respectively, supporting the scenario of resistant disease.

Directed screening.

For the quantitative PCR (qPCR) study, 125 additional patients were retrospectively recruited. The patients had clinical stages II and III adenocarcinoma (defined by endoscopic ultrasound and/or MRI) of the rectum within 10 cms from the anal verge and were treated with 2 or more preoperative cycles of chemotherapy, concomitant with 50.4 Gy in 28 fractions. The specified clinical and pathologic variables were recorded (Table 1). Thirty-one patients were excluded because of unsuitable pathologic samples or incomplete gene expression or clinical data. Ninety-four patients were selected for the response prediction statistical analysis. Patients received 4 different chemotherapy regimens as follows: (a) 29 UFT (300 mg/m2 on days 1–14 concurrent with radiotherapy) and leucovorin (250 mg/m2 i.v. on day 1 and 7.5 mg/m2/12 h p.o. on days 2–14; ref. 20; b) 17 5FU 225 mg/m2/d in continuous infusion during RT; 48 patients were treated with concurrent oxaliplatin-based CRT, including 27 (c) receiving raltitrexed and oxaliplatin (3 mg/m2 and 130 mg/m2, respectively, on days 1, 21, and 42), and 21 (d) 2 cycles of standard XELOX followed capecitabine (825 mg/m2/12 h) on radiotherapy days. Follow-up occurred at 3-month intervals for 2 years, 6-month intervals the third year, and yearly thereafter. Evaluations consisted of anamnesis, physical examination, and blood tests. Complete colonoscopy, chest radiography, abdominal ultrasound, and abdominopelvic CT were also scheduled according to local guidelines. Local and distance recurrence was histologically confirmed whenever possible. Alternate criteria included sequential enlargement of a mass or distant nodules in radiologic studies.

Figure 2.

DFS according to the score determined by the 13-gene signature (≤2 high risk, versus >2 low risk).

Figure 2.

DFS according to the score determined by the 13-gene signature (≤2 high risk, versus >2 low risk).

Close modal

Institutional approval from our ethical committee was obtained for the conduct of the study. Patients provided written consent so that their samples and clinical data could be used for investigational purposes.

Response and TRG

Two independent expert surgical pathologists who were blinded to patient outcome evaluated TRG according to the classification of Dworak, as described by Rodel and colleagues (19, 21). Response to CRT was considered as TRG greater than 0–1. The kappa index was used to evaluate concordance between the 2 observers and the pathologists discussed discordant results.

Construction of SAGE libraries and analysis

To identify genes involved in RC treatment response, a previous exploratory SAGE study was done. SAGE libraries were generated from initial tumor biopsies before CRT and postoperative specimens from patients A, B, and C. Tumor specimens were snap frozen in liquid nitrogen at −70°C and stained with hematoxylin and eosin (H&E). Eligible samples were composed of at least 90% of tumor cells. SAGE was done as described in the Supplementary Material. The pairwise comparisons (before and after treatment for each patient) were determined by using Monte Carlo simulation to identify transcripts in which difference was statistically different at P ≤ 0.01 and P ≤ 0.05 (Supplementary Tables S1 and S2). Gene matches for significant tags were manually verified by using SAGE Genie (22). Eligibility of SAGE genes for qPCR analysis was based on statistically significant altered expression, similar regulation trend (up/downregulation) across the 3 patients, and biological plausibility according to the published literature.

RNA extraction and qPCR

FFPE tumor sections from the initial biopsies of the 125 patients were reviewed by an expert pathologist. Ninety-four patients were selected after excluding cases with insufficient tumor cells, incomplete clinical data, or unsuccessful amplification results. More than 75% tumor cells enrichment was ensured, when necessary, by subsequent macrodissection with the use of a safety blade and new confirmatory H&E staining in postoperative biopsies. RNA was extracted from 10 to 15 five-micrometer sections by using the Masterpure RNA Purification Kit (Epicentre Biotechnologies) according to the manufacturer's instructions, detailed in Supplementary Material.

qPCR was done by using the ABI Prims 7900 HT Sequence Detection System (Applied Biosystems), according to manufacturer's instructions. TLDAs design was done as described in the Supplementary Material.

Statistical analysis

Testing for associations between response and clinical variables was done with χ2 test (Yates corrected) or Fisher's exact tests for categorical variables.

The expression of each gene was measured in triplicate, and the average intensity in the microarray was assessed. To reduce variation between microarrays, the intensity values for each sample in each microarray were rescaled by means of a quantile normalization method. Each intensity value was then log-transformed to a base-2 scale. qPCR data results were normalized against the 3 most stable genes (GAPDH, B2M, and PMSB4) before analysis. The stability of the expression of the candidate control genes was determined with genNorm (23).

The total sample size was randomly assigned to either a training cohort (60%) or a testing cohort (40%; Supplementary Fig. S1). Patients were classified into 3 groups (TRG 0–1, TRG2–3, and TRG4). We used a supervised class prediction by stepwise forwards and backwards discriminant analysis method to define a subset of genes predictive of response. This model was applied to the testing cohort to validate its usefulness. Accuracy, sensitivity, specificity, positive predictive values (PPV), and negative predictive values (NPV) were calculated (more detailed in Supplementary Materials).

Two-factor analyses were calculated (F1 and F2, the discriminant functions) for the genes that were selected by discriminant analysis. A score was obtained for every patient with the formula:

$${\rm Score} = \sqrt {F_1^2 + F_2^2}$$

To determine differences in DFS and overall survival (OS; measured from the time of diagnosis), univariate analyses were done according to the Kaplan–Meier method. Comparisons between curves were evaluated with the log-rank test. Cox regression was used for the multivariate analysis. A value of P < 0.05 was considered statistically significant, and all tests were 2-tailed. The proportional hazards assumption for covariates was investigated by examining the scaled Schoenfeld residuals. SPSS 17.0 and SLPUS v. 6. were used for statistical analysis.

SAGE screening of transcriptome changes after CRT

Six SAGE-libraries were generated from patients A, B, and C. Every patient was studied before and after treatment. Pairwise comparisons revealed statistically significant dynamic changes in gene expression when libraries were compared before and after treatment for each patient (P < 0.01). The numbers of tags, unique tags for each library, and dynamic changes are shown in Supplementary Tables S1 and S2. Certain tag-corresponding genes significantly altered after treatment included tissue specific genes such as the trefoil-family factor 3 (TFF3) and the gastrointestinal glutathione peroxidase 2 (GPX2) or transcripts associated with treatment response as glutathione S-transferase pi or with fibrotic response after injury, such as collagen type I alpha 1 (COL1A1) and secreted protein acidic cysteine-rich (SPARC). For qPCR analysis, we focused on 24 genes that were significantly regulated in pairwise comparisons (Supplementary Table S2) in addition to 29 genes previously reported to be involved in colorectal carcinogenesis or treatment response (Supplementary Table S3).

Clinicopathologic data for qPCR study

A total of 94 patients were analyzed (77 patients from Hospital La Paz, Madrid and 17 from Hospital Clinic of Barcelona). A flow diagram of the patients is depicted in Figure 1. Clinicopathologic data are shown in Table 1. After a median follow-up of 46 months (range 6–127 months), local relapse occurred in 6 patients (7%) and distant metastases occurred in 22 (23%). The 3-year DFS was 73% (95% CI: 64–82) and 3-year OS was 90% (95% CI: 84–96). TRG 0–1, 2, 3, and 4 were seen in 18, 53, 11, and 12 patients, respectively, with good concordance between pathologists (Kappa = 0.7).

Table 1.

Clinicopathologic characteristics

CharacteristicAll, n (%)TRG 01, n (%)TRG 2/3, n (%)TRG 4, n (%)P
Total 94 (100) 18 (19.1) 64 (68.1) 12 (12.8)  
Gender 
 Male 54 (57) 11 (61) 37 (57) 6 (50) 0.83 
 Female 40 (43) 7 (39) 27 (42) 6 (50)  
Age, y 
 ≤65 51 (54) 9 (50) 33 (52) 9 (75) 0.30 
 >65 43 (46) 9 (50) 31 (48) 3 (25)  
ECOG 
 0 51 (54) 8 (44) 36 (56) 7 (58) 0.84 
 1 40 (43) 9 (50) 26 (41) 5 (42)  
 2 3 (3) 1 (6) 2 (3)  
Concurrent chemotherapy 
 Fluoropyrimidines 46 (49) 14 (78) 27 (42) 5 (42) 0.02 
 Oxaliplatin based 48 (51) 4 (22) 37 (58) 7 (58)  
CEA 
 ≤2.5 ng/mL 34 (36) 8 (44) 18 (28) 8 (67) 0.06 
 >2.5 ng/mL 53 (56) 10 (56) 39 (61) 4 (33)  
 Unknown 7 (7) 0 (0) 7 (11) (0)  
Clinical stage (AJCC TNM) 
 I 4 (4) 0 (0) 3 (5) 1 (8) 0.73 
 II 49 (52) 11 (61) 33 (52) 5 (42)  
 III 41 (43) 7 (39) 28 (44) 6 (50)  
Pathologic stage (AJCC TNM) 
 0 12 (13) 0 (0) 0 (0) 12 (100) <.01 
 I 24 (26) 4 (22) 20 (31) 0 (0)  
 II 40 (43) 7 (39) 33 (52) 0 (0)  
 III 18 (19) 7 (39) 11 (17) 0 (0)  
Score 
 >2 32 (34) 9 (50) 17 (27) 6 (50) 0.08 
 ≤2 62 (66) 9 (50) 47 (73) 6 (50)  
Relapse 
 No 66 (70) 11 (61) 44 (69) 11 (92) 0.43 
 Local 6 (6) 2 (11) 4 (6) 0 (0)  
 Distant 22 (23) 5 (28) 16 (25) 1 (8)  
Death      
 No 80 (85) 14 (78) 54 (84) 12 (100) 0.24 
 Yes 14 (15) 4 (22) 10 (16) 0 (0)  
CharacteristicAll, n (%)TRG 01, n (%)TRG 2/3, n (%)TRG 4, n (%)P
Total 94 (100) 18 (19.1) 64 (68.1) 12 (12.8)  
Gender 
 Male 54 (57) 11 (61) 37 (57) 6 (50) 0.83 
 Female 40 (43) 7 (39) 27 (42) 6 (50)  
Age, y 
 ≤65 51 (54) 9 (50) 33 (52) 9 (75) 0.30 
 >65 43 (46) 9 (50) 31 (48) 3 (25)  
ECOG 
 0 51 (54) 8 (44) 36 (56) 7 (58) 0.84 
 1 40 (43) 9 (50) 26 (41) 5 (42)  
 2 3 (3) 1 (6) 2 (3)  
Concurrent chemotherapy 
 Fluoropyrimidines 46 (49) 14 (78) 27 (42) 5 (42) 0.02 
 Oxaliplatin based 48 (51) 4 (22) 37 (58) 7 (58)  
CEA 
 ≤2.5 ng/mL 34 (36) 8 (44) 18 (28) 8 (67) 0.06 
 >2.5 ng/mL 53 (56) 10 (56) 39 (61) 4 (33)  
 Unknown 7 (7) 0 (0) 7 (11) (0)  
Clinical stage (AJCC TNM) 
 I 4 (4) 0 (0) 3 (5) 1 (8) 0.73 
 II 49 (52) 11 (61) 33 (52) 5 (42)  
 III 41 (43) 7 (39) 28 (44) 6 (50)  
Pathologic stage (AJCC TNM) 
 0 12 (13) 0 (0) 0 (0) 12 (100) <.01 
 I 24 (26) 4 (22) 20 (31) 0 (0)  
 II 40 (43) 7 (39) 33 (52) 0 (0)  
 III 18 (19) 7 (39) 11 (17) 0 (0)  
Score 
 >2 32 (34) 9 (50) 17 (27) 6 (50) 0.08 
 ≤2 62 (66) 9 (50) 47 (73) 6 (50)  
Relapse 
 No 66 (70) 11 (61) 44 (69) 11 (92) 0.43 
 Local 6 (6) 2 (11) 4 (6) 0 (0)  
 Distant 22 (23) 5 (28) 16 (25) 1 (8)  
Death      
 No 80 (85) 14 (78) 54 (84) 12 (100) 0.24 
 Yes 14 (15) 4 (22) 10 (16) 0 (0)  

Abbreviations: AJCC, American Joint Committee on Cancer; ECOG, Eastern Cooperative Oncology Group; UFT-LV, UFT and leucovorin.

Development of a response predictor

From the 53 genes, 13 were selected by discriminant analysis to build a predictor of response (Table 2). In the training cohort, the model correctly classified 86% of the cases (49 of 57). Two nonresponsive patients (TRG 0/1) and 6 responsive cases (TRG 2, 3, and 4) were misclassified. The test correctly classified 87% of responders (sensitivity) and 82% of nonresponders (specificity). The PPV was 95% and the NPV was 60%. In the testing cohort, the model correctly classified 6 of 7 nonresponders (specificity = 86%) and 22 of 30 responders (sensitivity = 73%), with an overall accuracy of 76% (Supplementary Fig. S1) and a PPV of 95%.

Table 2.

13-Gene signature

GeneGene originUnigene numberExpression in responders compared with nonrespondersDescription
ALDH1A1 Hs.76392 Overexpressed Aldehyde dehydrogenase 1 family, member A1 
BAK1 Hs.485139 Downregulated BCL2-antagonist/killer 1 
CDKN1 Hs.370771 Overexpressed Cyclin-dependent kinase inhibitor 1A (p21, Cip1) 
FOS Hs.25647 Overexpressed V-fos FBJ murine osteosarcoma viral oncogene homolog 
MLH1 Hs.195364 Downregulated MutL homolog 1, colon cancer, nonpolyposis type 2 (E. coli
RELB Hs.654402 Overexpressed V-rel reticuloendotheliosis viral oncogene homolog B 
STAT3 Hs.463059 Overexpressed Signal transducer and activator of transcription 3 (acute-phase response factor) 
TYMS Hs.592338 Downregulated Thymidylate synthetase 
CKB Hs.173724 Downregulated Creatine kinase, brain 
GPX2 Hs.2704 Downregulated Glutathione peroxidase 2 (gastrointestinal) 
HIG2 Hs.710088 Downregulated Hypoxia-inducible protein 2 
PH-4 Hs.654944 Downregulated Hypoxia-inducible factor prolyl 4-hydroxylase 
TFF3 Hs.82961 Overexpressed Trefoil factor 3 (intestinal) 
GeneGene originUnigene numberExpression in responders compared with nonrespondersDescription
ALDH1A1 Hs.76392 Overexpressed Aldehyde dehydrogenase 1 family, member A1 
BAK1 Hs.485139 Downregulated BCL2-antagonist/killer 1 
CDKN1 Hs.370771 Overexpressed Cyclin-dependent kinase inhibitor 1A (p21, Cip1) 
FOS Hs.25647 Overexpressed V-fos FBJ murine osteosarcoma viral oncogene homolog 
MLH1 Hs.195364 Downregulated MutL homolog 1, colon cancer, nonpolyposis type 2 (E. coli
RELB Hs.654402 Overexpressed V-rel reticuloendotheliosis viral oncogene homolog B 
STAT3 Hs.463059 Overexpressed Signal transducer and activator of transcription 3 (acute-phase response factor) 
TYMS Hs.592338 Downregulated Thymidylate synthetase 
CKB Hs.173724 Downregulated Creatine kinase, brain 
GPX2 Hs.2704 Downregulated Glutathione peroxidase 2 (gastrointestinal) 
HIG2 Hs.710088 Downregulated Hypoxia-inducible protein 2 
PH-4 Hs.654944 Downregulated Hypoxia-inducible factor prolyl 4-hydroxylase 
TFF3 Hs.82961 Overexpressed Trefoil factor 3 (intestinal) 

1: Bibliography; 2: SAGE analysis.

Score for DFS

The 13-gene score was applied to the entire cohort of patients and was inversely correlated with DFS. Prognostic groups were defined as low (>2) and high risk (≤2) (Figure 2). Clinical variables, except for basal carcinoembryonic antigen level (CEA), did not differ between risks groups. On univariate analysis, the risk score was significantly correlated with DFS along with pathologic variables [AJCC Stage, ypT (post–CRT-yatrogenic pathologic T) and ypN (post–CRT-yatrogenic pathologic N; Table 3)]. Given the heterogeneity of the sample, we investigated whether the type of chemotherapy received impacted the scoring results. When stratifying by treatment, the score remained statistically significant (P = 0.03 and P = 0.01 for fluoropyrimidines alone and oxaliplatin-based CRT, respectively; Supplementary Fig. S2). Univariate stratified analysis for other factors, such as ypT and ypN, confirmed the independent prognostic value of the score for DFS (Supplementary Figs. S3 and S4). Interestingly, the score was such a strong predictor of DFS that it remained the only statistically significant factor when Cox multivariate analysis was used, adjusted for age, ypT, ypN, CEA, and tumor response (HR = 4.35 95% CI: 1.20–15.75; Table 4).

Table 3.

Univariate analysis for DFS

FactorNo. of patients3-year DFS, % (95% CI)PHR (95% CI)P
All 94 73 (64–82)    
Age, y 
 ≤65 51 70 (57–83) 0.52 1.28 (0.6–2.7) 0.52 
 >65 43 77 (64–90)    
ypT 
 T0 14 92 (76–100) 0.15 2.03 (1.1–3.5)a 0.01 
 T1 100 (–)    
 T2 21 74 (55–94)    
 T3 53 65 (51–78)    
 T4 100 (–)    
ypN 
 N0 76 79 (69–88) .01 0.39 (0.18–0.85) 0.02 
 N+ 18 53 (29–77)    
AJCC TNM pathologic stage 
 0 12 90 (71–100) 0.06 1.75 (1.1–2.7)b 0.01 
 I 24 78 (61–95)    
 II 40 75 (61–90)    
 III 18 53 (29–77)    
Response to treatment 
 No (TRG 0/1) 18 63 (39–86) 0.41 0.41 
 Yes (TRG 2,3,4) 76 75 (65–85)  0.7 (0.3–1.6)  
CEA level 
 0–2.5 ng/mL 34 84 (71–97) 0.15 0.54 (0.22–1.28) 0.16 
 >2.5 ng/mL 53 69 (56–82)    
Risk score 
 ≤2 62 64 (52–77) <0.01   
 >2 32 90 (79–100)  0.185 (0.06–0.61) <0.01 
FactorNo. of patients3-year DFS, % (95% CI)PHR (95% CI)P
All 94 73 (64–82)    
Age, y 
 ≤65 51 70 (57–83) 0.52 1.28 (0.6–2.7) 0.52 
 >65 43 77 (64–90)    
ypT 
 T0 14 92 (76–100) 0.15 2.03 (1.1–3.5)a 0.01 
 T1 100 (–)    
 T2 21 74 (55–94)    
 T3 53 65 (51–78)    
 T4 100 (–)    
ypN 
 N0 76 79 (69–88) .01 0.39 (0.18–0.85) 0.02 
 N+ 18 53 (29–77)    
AJCC TNM pathologic stage 
 0 12 90 (71–100) 0.06 1.75 (1.1–2.7)b 0.01 
 I 24 78 (61–95)    
 II 40 75 (61–90)    
 III 18 53 (29–77)    
Response to treatment 
 No (TRG 0/1) 18 63 (39–86) 0.41 0.41 
 Yes (TRG 2,3,4) 76 75 (65–85)  0.7 (0.3–1.6)  
CEA level 
 0–2.5 ng/mL 34 84 (71–97) 0.15 0.54 (0.22–1.28) 0.16 
 >2.5 ng/mL 53 69 (56–82)    
Risk score 
 ≤2 62 64 (52–77) <0.01   
 >2 32 90 (79–100)  0.185 (0.06–0.61) <0.01 

aReference value for HR: ypT increment.

bReference value for HR: AJCC TNM pathologic stage increment.

Table 4.

Multivariate analysis for DFS

HR (95% CI)P
Age (≤65 vs. >65) 1.25 (0.55–2.82) 0.59 
CEA (0–2.5 vs. >2.5 ng/mL) 0.88 (0.35–2.24) 0.80 
ypT (T increment) 1.60 (0.96–2.69) 0.07 
ypN (negative vs. positive) 0.65 (0.27–1.58) 0.58 
Response to treatment (yes vs. no)a 0.76 (0.29–1.98) 0.59 
Risk score (≤2 vs. >2) 4.35 (1.20–15.75) 0.02 
HR (95% CI)P
Age (≤65 vs. >65) 1.25 (0.55–2.82) 0.59 
CEA (0–2.5 vs. >2.5 ng/mL) 0.88 (0.35–2.24) 0.80 
ypT (T increment) 1.60 (0.96–2.69) 0.07 
ypN (negative vs. positive) 0.65 (0.27–1.58) 0.58 
Response to treatment (yes vs. no)a 0.76 (0.29–1.98) 0.59 
Risk score (≤2 vs. >2) 4.35 (1.20–15.75) 0.02 

aResponse to treatment assessed as TRG 2, 3, 4.

Treatment and tumor-related factors determined after surgery, primarily CRM and nodal status, have been associated with RC recurrence (5). However, individualized treatment selection for patients with RC relies on pretreatment clinical variables that do not predict the efficacy of neoadjuvant therapy. This study used a sequential genomic strategy to identify a basal 13-gene signature that correctly predicted RC preoperative CRT response in 86% of patients in a training set and 6 of 7 nonresponders in a confirmatory testing cohort. Also, it predicted 2 groups of patients with clear differences in DFS.

This predictor model offers some advantages over other published tumor gene signatures that evaluate response to preoperative RT (24) or combined CRT (10, 12, 13, 25) in terms of clinical feasibility, relevance of the surrogate marker, and prognostic ability. This is the only signature obtained from small amounts of routinely prepared FFPE endoscopic samples by using a qPCR platform. In addition, the endpoint of tumor downstaging used by Ghadimi and colleagues (10) is subjected to numerous evaluation pitfalls. In contrast, TRG provides direct biological information about cytotoxicity, clearly beyond downstaging (26), and it seems to correlate with prognosis in both RC (21, 27–30) and other tumors (31–35). Furthermore, the 13-gene signature is the only one that predicts Dworak's TRG 0–1, which is the group containing the poorest responders. These patients likely have a bad prognosis and would be good candidates for alternative strategies to standard CRT, such as new drugs, doses, and time strategies. TRG 0–1 accounted for 25% of the patients in the largest data set from a randomized trial (21) and 19.1% of the patients in our series. Conversely, identifying patients with some grade of response (95% PPV in the training and testing cohorts) would be very important to further explore the promising strategy of induction chemotherapy that has been recently supported by the GEMCAD GCR3 randomized phase II trial (36). It would exclude those patients who do not benefit from CRT and prevent them from delaying surgical definitive treatment. Because CRT is standard treatment for stage II/III RC, the clinical interest of predicting Dworak's TRG 1–2–3 versus TRG 4 or patients with less than 10% tumor viable cells, as reported by Kim (12) and Rimkus (13), respectively, is less evident. In contrast, a new kinase profile has shown correlation between basal kinase activity and response to CRT in good, intermediate, and poor responders, although it was only validated among good and intermediate responders (25). Besides, these signatures have not been associated with prognosis. Our signature prognostic score remained the only statistically significant factor, even above postsurgical ypT and ypN, following multivariate analysis (HR = 4.35 95% CI: 1.2–15.75, P = 0.02). In our series, the TRG predictive signature-derived score, but not the TRG classification, correlated with outcome. A conceivable explanation for this phenomenon could be that the molecular portrait retrieves deeper informative biological aspects beyond TRG. In fact, a prognostic versus merely predictive capability cannot be excluded.

Limitations of our study include the small number of patients comprising the training and testing sets (57 and 37, respectively) and the heterogeneous treatment schedules. However, this study was larger than other RT/CRT predictive signatures published in RC. All patients received combined CRT with at least 2 cycles of chemotherapy. Therefore, it is not possible to determine whether their response was the result of either treatment component or the combined therapy. Nonetheless, combined CRT is the usual strategy, and the score identified significant DFS differences in both 5FU and oxaliplatin-based CRT groups (P = 0.03 and P = 0.015, respectively). Of note, recent data from randomized trials do not support the use of oxaliplatin during concurrent preoperative CRT in RC (37, 38).

In 2003, Crane and colleagues underlined an interest in identifying altered marker expression in response to preoperative treatment (39). We present the first global “dynamic” transcriptome study. Importantly, this study was done in CRT refractory patients by using SAGE, a powerful genomic tool. This labor-intensive method is very quantitative and is not restricted to a finite number of gene sequences localized to a chip. Although the exploratory study was restricted to 3 selected patients (cT3N0, TRG 0–1), the “intrapatient” changes in the Monte Carlo analysis provided plausible biological data. The majority of regulated genes were involved in oncogenic functions, and one-fourth of the observed changes (regulation of 24 genes) were shared among the 3 cases. These data highlighted interesting candidates for prognostic and functional evaluations. Moreover, regulated tags included an interesting 10% of expressed sequence tags.

The CRT-induced fibrosis in the surgical specimens was mirrored at the molecular level by upregulation of genes involved in the healing process. These genes included COL1A1, the key enzyme in the synthesis of collagens, PH-4 (40), and the TFF3 protein, which is involved in mucous restitution (41). The production of free radicals and hypoxia resulted in the regulation of GPX2 (42), HIG-2 (43), and angiopoietin-like 4 (44). Other genes that were consistently regulated have roles in cell-cycle proliferation and apoptosis and include CCNB1IP1 (45) and RCC1 (46). EEF2, an essential factor for protein synthesis that plays a role in gastrointestinal carcinogenesis (47), and WARS, a gene associated with poor prognosis in colorectal cancer (48), were also identified. Genes participating in the metastatic process were identified, including MMP14 (49), integrin B4 (50), and SPARC, a key controller of cell–matrix interactions that is also associated with tumor angiogenesis and metastasis (51, 52).

The SAGE screen of candidate genes added important value to the development of the gene predictor set. SAGE provided 5 of the 13 predictive genes (CKB, GPX2, HIG-2, PH-4, and TFF3) and a relatively small number of genes was required as compared with the other RC gene signatures [13 vs. 33 (ref. 11), 42 (ref. 13), 54 (ref. 10), 86 (ref. 25), and 95 (12)]. In general, genes in the signature suggest cell defense against CRT injury. This is exemplified by the downregulation of cell proliferation and apoptosis moderators, including the antiapoptotic BAK1 (53) gene and the antioxidant GPX2 (42), and the overexpression of the growth inhibitor CDKN1 (54) in responders. Consistent with most reports, thymidylate synthetase was overexpressed in nonresponders (55). Notably, only GPX2 has been previously reported as a component of other RC signatures (24). The lack of consistency among RC signatures is indicative of the complexity of resistance pathways in response to treatment injury.

In summary, by using a sequential genomic approach, we have developed a clinically exploitable gene signature that eventually may help select RC patients for customized treatment. Our results need to be confirmed and retrospective validation in well-controlled prospective clinical trials is currently ongoing by GEMCAD.

No potential conflicts of interest were disclosed.

This work has been funded by a FIS PI021094 grant. P. Cejas was supported by Asociación Española Contra el Cáncer (AECC) programa PAO 2010. V.M. Garcia was supported by Fundación Para la Investigación Biomédica del Hospital La Paz (FIBHULP) grant: REX 09.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Sauer
R
,
Becker
H
,
Hohenberger
W
,
Rödel
C
,
Wittekind
C
,
Fietkau
R
, et al
Preoperative versus postoperative chemoradiotherapy for rectal cancer
.
N Engl J Med
2004
;
351
:
1731
40
.
2.
Roh
MS
,
Colangelo
LH
,
O'Connell
MJ
,
Yothers
G
,
Deutsch
M
,
Allegra
CJ
, et al
Preoperative multimodality therapy improves disease-free survival in patients with carcinoma of the rectum: NSABP R-03
.
J Clin Oncol
2009
;
27
:
5124
30
.
3.
Bosset
JF
,
Collette
L
,
Calais
G
,
Mineur
L
,
Maingon
P
,
Radosevic-Jelic
L
, et al
Chemotherapy with preoperative radiotherapy in rectal cancer
.
N Engl J Med
2006
;
355
:
1114
23
.
4.
Gérard
JP
,
Conroy
T
,
Bonnetain
F
,
Bouché
O
,
Chapet
O
,
Closon-Dejardin
MT
, et al
Preoperative radiotherapy with or without concurrent fluorouracil and leucovorin in T3–4 rectal cancers: results of FFCD 9203
.
J Clin Oncol
2006
;
24
:
4620
5
.
5.
Nagtegaal
ID
,
Gosens
MJ
,
Marijnen
CA
,
Rutten
HJ
,
van d
V
,
van Krieken
JH
. 
Combinations of tumor and treatment parameters are more discriminative for prognosis than the present TNM system in rectal cancer
.
J Clin Oncol
2007
;
25
:
1647
50
.
6.
Quirke
P
,
Steele
R
,
Monson
J
,
Grieve
R
,
Khanna
S
,
Couture
J
, et al
Effect of the plane of surgery achieved on local recurrence in patients with operable rectal cancer: a prospective study using data from the MRC CR07 and NCIC-CTG CO16 randomised clinical trial
.
Lancet
2009
;
373
:
821
8
.
7.
Smith
FM
,
Reynolds
JV
,
Miller
N
,
Stephens
RB
,
Kennedy
MJ
. 
Pathological and molecular predictors of the response of rectal cancer to neoadjuvant radiochemotherapy
.
Eur J Surg Oncol
2006
;
32
:
55
64
.
8.
Johnston
PG
. 
Prognostic markers of local relapse in rectal cancer: are we any further forward?
J Clin Oncol
2006
;
24
:
4049
50
.
9.
Bertolini
F
,
Bengala
C
,
Losi
L
,
Pagano
M
,
Iachetta
F
,
Dealis
C
, et al
Prognostic and predictive value of baseline and posttreatment molecular marker expression in locally advanced rectal cancer treated with neoadjuvant chemoradiotherapy
.
Int J Radiat Oncol Biol Phys
2007
;
68
:
1455
61
.
10.
Ghadimi
BM
,
Grade
M
,
Difilippantonio
MJ
,
Varma
S
,
Simon
R
,
Montagna
C
, et al
Effectiveness of gene expression profiling for response prediction of rectal adenocarcinomas to preoperative chemoradiotherapy
.
J Clin Oncol
2005
;
23
:
1826
38
.
11.
Watanabe
T
,
Komuro
Y
,
Kiyomatsu
T
,
Kanazawa
T
,
Kazama
Y
,
Tanaka
J
, et al
Prediction of sensitivity of rectal cancer cells in response to preoperative radiotherapy by DNA microarray analysis of gene expression profiles
.
Cancer Res
2006
;
66
:
3370
4
.
12.
Kim
IJ
,
Lim
SB
,
Kang
HC
,
Chang
HJ
,
Ahn
SA
,
Park
HW
, et al
Microarray gene expression profiling for predicting complete response to preoperative chemoradiotherapy in patients with advanced rectal cancer
.
Dis Colon Rectum
2007
;
50
:
1342
53
.
13.
Rimkus
C
,
Friederichs
J
,
Boulesteix
AL
,
Theisen
J
,
Mages
J
,
Becker
K
, et al
Microarray-based prediction of tumor response to neoadjuvant radiochemotherapy of patients with locally advanced rectal cancer
.
Clin Gastroenterol Hepatol
2008
;
6
:
53
61
.
14.
Paik
S
,
Shak
S
,
Tang
G
,
Kim
C
,
Baker
J
,
Cronin
M
, et al
A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer
.
N Engl J Med
2004
;
351
:
2817
26
.
15.
Harris
L
,
Fritsche
H
,
Mennel
R
,
Norton
L
,
Ravdin
P
,
Taube
S
, et al
American Society of Clinical Oncology 2007 update of recommendations for the use of tumor markers in breast cancer
.
J Clin Oncol
2007
;
25
:
5287
312
.
16.
NCCN Clinical practice guidelines in oncology breast cancer
(
V.2.2010
). 
2010
.
17.
Velculescu
VE
,
Zhang
L
,
Vogelstein
B
,
Kinzler
KW
. 
Serial analysis of gene expression
.
Science
1995
;
270
:
484
7
.
18.
Casado
E
,
De Castro
J
,
Castelo
B
,
Sanchez
B
,
Feliu
J
,
Escribano
A
, et al
Phase II study of neoadjuvant treatment of rectal cancer with oxaliplatin, raltitrexed and radiotherapy
.
J Clin Oncol
22
:
14s
2004
(
suppl: abstr 3746
).
19.
Dworak
O
,
Keilholz
L
,
Hoffmann
A
. 
Pathological features of rectal cancer after preoperative radiochemotherapy
.
Int J Colorectal Dis
1997
;
12
:
19
23
.
20.
Feliu
J
,
Calvilio
J
,
Escribano
A
,
de Castro
J
,
Sánchez
ME
,
Mata
A
, et al
Neoadjuvant therapy of rectal carcinoma with UFT-leucovorin plus radiotherapy
.
Ann Oncol
2002
;
13
:
730
6
.
21.
Rödel
C
,
Martus
P
,
Papadoupolos
T
,
Füzesi
L
,
Klimpfinger
M
,
Fietkau
R
, et al
Prognostic significance of tumor regression after preoperative chemoradiotherapy for rectal cancer
.
J Clin Oncol
2005
;
23
:
8688
96
.
22.
Lal
A
,
Lash
AE
,
Altschul
SF
,
Velculescu
V
,
Zhang
L
,
McLendon
RE
, et al
A public database for gene expression in human cancers
.
Cancer Res
1999
;
59
:
5403
7
.
23.
Vandesompele
J
,
De Preter
K
,
Pattyn
F
,
Poppe
B
,
Van Roy
N
,
De Paepe
A
, et al
Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes
.
Genome Biol
2002
;
3
:
34.1
34.11
.
24.
Watanabe
T
,
Komuro
Y
,
Kiyomatsu
T
,
Kanazawa
T
,
Kazama
Y
,
Tanaka
J
, et al
Prediction of sensitivity of rectal cancer cells in response to preoperative radiotherapy by DNA microarray analysis of gene expression profiles
.
Cancer Res
2006
;
66
:
3370
4
.
25.
Folkvord
S
,
Flatmark
K
,
Dueland
S
,
de Wijn
R
,
Grøholt
KK
,
Hole
KH
, et al
Prediction of response to preoperative chemoradiotherapy in rectal cancer by multiplex kinase activity profiling
.
Int J Radiat Oncol Biol Phys
2010
;
78
:
555
62
.
26.
Vironen
J
,
Juhola
M
,
Kairaluoma
M
,
Jantunen
I
,
Kellokumpu
I
. 
Tumour regression grading in the evaluation of tumour response after different preoperative radiotherapy treatments for rectal carcinoma
.
Int J Colorectal Dis
2005
;
20
:
440
5
.
27.
Bouzourene
H
,
Bosman
FT
,
Seelentag
W
,
Matter
M
,
Coucke
P
. 
Importance of tumor regression assessment in predicting the outcome in patients with locally advanced rectal carcinoma who are treated with preoperative radiotherapy
.
Cancer
2002
;
94
:
1121
30
.
28.
Losi
L
,
Luppi
G
,
Gavioli
M
,
Iachetta
F
,
Bertolini
F
,
D'Amico
R
, et al
Prognostic value of Dworak grade of regression (GR) in patients with rectal carcinoma treated with preoperative radiochemotherapy
.
Int J Colorectal Dis
2006
;
21
:
645
51
.
29.
Lindebjerg
J
,
Spindler
KL
,
Ploen
J
,
Jakobsen
A
. 
The prognostic value of lymph node metastases and tumour regression grade in rectal cancer patients treated with long-course preoperative chemoradiotherapy
.
Colorectal Dis
2009
;
11
:
264
9
.
30.
Guillem
JG
,
Chessin
DB
,
Cohen
AM
,
Shia
J
,
Mazumdar
M
,
Enker
W
, et al
Long-term oncologic outcome following preoperative combined modality therapy and total mesorectal excision of locally advanced rectal cancer
.
Ann Surg
2005
;
241
:
829
36
.
31.
Mandard
AM
,
Dalibard
F
,
Mandard
JC
,
Marnay
J
,
Henry-Amar
M
,
Petiot
JF
, et al
Pathologic assessment of tumor regression after preoperative chemoradiotherapy of esophageal carcinoma. Clinicopathologic correlations
.
Cancer
1994
;
73
:
2680
6
.
32.
Becker
K
,
Mueller
JD
,
Schulmacher
C
,
Ott
K
,
Fink
U
,
Busch
R
, et al
Histomorphology and grading of regression in gastric carcinoma treated with neoadjuvant chemotherapy
.
Cancer
2003
;
98
:
1521
30
.
33.
Lowy
AM
,
Mansfield
PF
,
Leach
SD
,
Pazdur
R
,
Dumas
P
,
Ajani
JA
. 
Response to neoadjuvant chemotherapy best predicts survival after curative resection of gastric cancer
.
Ann Surg
1999
;
229
:
303
8
.
34.
Rödel
C
,
Grabenbauer
GG
,
Kühn
R
,
Papadopoulos
T
,
Dunst
J
,
Meyer
M
, et al
Combined-modality treatment and selective organ preservation in invasive bladder cancer: long-term results
.
J Clin Oncol
2002
;
20
:
3061
71
.
35.
Braun
OM
,
Neumeister
B
,
Popp
W
,
Scherrer
R
,
Dobrowsky
E
,
Dobrowsky
W
, et al
Histologic tumor regression grades in squamous cell carcinoma of the head and neck after preoperative radiochemotherapy
.
Cancer
1989
;
63
:
1097
100
.
36.
Fernández-Martos
C
,
Pericay
C
,
Aparicio
J
,
Salud
A
,
Safont
M
,
Massuti
B
, et al
Phase II, randomized study of concomitant chemoradiotherapy followed by surgery and adjuvant capecitabine plus oxaliplatin (CAPOX) compared with induction CAPOX followed by concomitant chemoradiotherapy and surgery in magnetic resonance imaging-defined, locally advanced rectal cancer: Grupo cancer de recto 3 study
.
J Clin Oncol
2010
;
28
:
859
65
.
37.
Gérard
JP
,
Azria
D
,
Gourgou-Bourgade
S
,
Martel-Laffay
I
,
Hennequin
C
,
Etienne
PL
, et al
Comparison of two neoadjuvant chemoradiotherapy regimens for locally advanced rectal cancer: results of the phase III trial ACCORD 12/0405-Prodige 2
.
J Clin Oncol
2010
;
28
:
1638
44
.
38.
Aschele
C
,
Pinto
C
,
Cordio
G
,
Rosati
G
,
Tagliagambe
A
,
Artale
S
, et al
Preoperative fluoruracil (FU)-based chemoradiation with and without weekly oxaliplatin in locally advanced rectal cancer: pathologic response analysis of the Studio Terapia Adiuvante Retto (STAR)-01 randomized phase III trial
.
J Clin Oncol
2009
;
27
.
39.
Crane
CH
,
Thames
HD
,
Hamilton
SR
. 
Will identifying or targeting altered marker expression in response to cytotoxic therapy be of prognostic or therapeutic value?
J Clin Oncol
2003
;
21
:
3381
2
.
40.
Myllyharju
J
. 
Prolyl 4-hydroxylases, key enzymes in the synthesis of collagens and regulation of the response to hypoxia, and their roles as treatment targets
.
Ann Med
2008
;
40
:
402
17
.
41.
Podolsky
DK
,
Lynch-Devaney
K
,
Stow
JL
,
Oates
P
,
Murgue
B
,
DeBeaumont
M
, et al
Identification of human intestinal trefoil factor. Goblet cell-specific expression of a peptide targeted for apical secretion
.
J Biol Chem
1993
;
268
:
6694
702
.
42.
Brigelius-Flohe
R
. 
Tissue-specific functions of individual glutathione peroxidases
.
Free Radic Biol Med
1999
;
27
:
951
65
.
43.
Denko
N
,
Schindler
C
,
Koong
A
,
Laderoute
K
,
Green
C
,
Giaccia
A
. 
Epigenetic regulation of gene expression in cervical cancer cells by the tumor microenvironment
.
Clin Cancer Res
2000
;
6
:
480
7
.
44.
Le Jan
S
,
Amy
C
,
Cazes
A
,
Monnot
C
,
Lamandé
N
,
Favier
J
, et al
Angiopoietin-like 4 is a proangiogenic factor produced during ischemia and in conventional renal cell carcinoma
.
Am J Pathol
2003
;
162
:
1521
8
.
45.
Grönholm
M
,
Muranen
T
,
Toby
GG
,
Utermark
T
,
Hanemann
CO
,
Golemis
EA
, et al
A functional association between merlin and HEI10, a cell cycle regulator
.
Oncogene
2006
;
25
:
4389
98
.
46.
Ohtsubo
M
,
Okazaki
H
,
Nishimoto
T
. 
The RCC1 protein, a regulator for the onset of chromosome condensation locates in the nucleus and binds to DNA
.
J Cell Biol
1989
;
109
:
1389
97
.
47.
Nakamura
J
,
Aoyagi
S
,
Nanchi
I
,
Nakatsuka
S
,
Hirata
E
,
Shibata
S
, et al
Overexpression of eukaryotic elongation factor eEF2 in gastrointestinal cancers and its involvement in G2/M progression in the cell cycle
.
Int J Oncol
2009
;
34
:
1181
9
.
48.
Ghanipour
A
,
Jirstrom
K
,
Ponten
F
,
Glimelius
B
,
Pahlman
L
,
Birgisson
H
. 
The prognostic significance of tryptophanyl-tRNA synthetase in colorectal cancer
.
Cancer Epidemiol Biomarkers Prev
2009
;
18
:
2949
56
.
49.
Fisher
KE
,
Sacharidou
A
,
Stratman
AN
,
Mayo
AM
,
Fisher
SB
,
Mahan
RD
, et al
MT1-MMP- and Cdc42-dependent signaling co-regulate cell invasion and tunnel formation in 3D collagen matrices
.
J Cell Sci
2009
;
122
:
4558
69
.
50.
Robertson
JH
,
Yang
SY
,
Winslet
MC
,
Seifalian
AM
. 
Functional blocking of specific integrins inhibit colonic cancer migration
.
Clin Exp Metastasis
2009
;
26
:
769
80
.
51.
Le
BB
,
Faouzi
S
,
Boussarie
L
,
Guirouilh
J
,
Blanc
JF
,
Carles
J
, et al
Osteonectin/SPARC is overexpressed in human hepatocellular carcinoma
.
J Pathol
1999
;
189
:
46
52
.
52.
Porte
H
,
Chastre
E
,
Prevot
S
,
Nordlinger
B
,
Empereur
S
,
Basset
P
, et al
Neoplastic progression of human colorectal cancer is associated with overexpression of the stromelysin-3 and BM-40/SPARC genes
.
Int J Cancer
1995
;
64
:
70
5
.
53.
Narita
M
,
Shimizu
S
,
Ito
T
,
Chittenden
T
,
Lutz
RJ
,
Matsuda
H
, et al
Bax interacts with the permeability transition pore to induce permeability transition and cytochrome c release in isolated mitochondria
.
Proc Natl Acad Sci U S A
1998
;
95
:
14681
6
.
54.
Schreiber
M
,
Muller
WJ
,
Singh
G
,
Graham
FL
. 
Comparison of the effectiveness of adenovirus vectors expressing cyclin kinase inhibitors p16INK4A, p18INK4C, p19INK4D, p21(WAF1/CIP1) and p27KIP1 in inducing cell cycle arrest, apoptosis and inhibition of tumorigenicity
.
Oncogene
1999
;
18
:
1663
76
.
55.
Popat
S
,
Matakidou
A
,
Houlston
RS
. 
Thymidylate synthase expression and prognosis in colorectal cancer: a systematic review and meta-analysis
.
J Clin Oncol
2004
;
22
:
529
36
.