Skip to main content

A study of gene expression markers for predictive significance for bevacizumab benefit in patients with metastatic colon cancer: a translational research study of the Hellenic Cooperative Oncology Group (HeCOG)



Bevacizumab, an antibody neutralizing Vascular Endothelial Growth Factor (VEGF), is licensed for the management of patients with advanced colon cancer. However, tumor biomarkers identifying the molecular tumor subsets most amenable to angiogenesis modulation are lacking.


We profiled expession of 24526 genes by means of whole genome 24 K DASL (c-DNA-mediated, Annealing, Selection and Ligation) arrays, (Illumina, CA) in 16 bevacizumab-treated patients with advanced colon cancer (Test set). Genes with correlation to 8-month Progression-free status were studied by means of qPCR in two independent colon cancer cohorts: 49 patients treated with bevacizumab + chemotherapy (Bevacizumab qPCR set) and 72 patients treated with chemotherapy only (Control qPCR set). Endpoints were best tumor response before metastasectomy (ORR) and progression-free survival (PFS).


Five genes were significantly correlated to 8-month progression-free status in the Test set: overexpression of KLF12 and downregulation of AGR2, ALDH6A1, MCM5, TFF2. In the two independent datasets, irinotecan- or oxaliplatin-based chemotherapy was administered as first-line treatment and metastasectomies were subsequently applied in 8-14% of patients. No prognostically significant gene classifier encompassing all five genes could be validated in the Bevacizumab or Control qPCR sets. The complex gene expression profile of all-low tumor (ALDH6A1 + TFF2 + MCM5) was strongly associated with ORR in the Bevacizumab qPCR set (ORR 85.7%, p = 0.007), but not in the Control set (ORR 36.4%, p = 0.747). The Odds Ratio for response for the all-low tumor (ALDH6A1 + TFF2 + MCM5) profile versus any other ALDH6A1 + TFF2 + MCM5 profile was 15 (p = 0.018) in the Bevacizumab qPCR set but only 0.72 (p = 0.63) in the Control set. The tumor expression profile of (KLF12-high + TFF2-low) was significantly associated with PFS only in the Bevacizumab qPCR set: bevacizumab-treated patients with (KLF12-high + TFF2-low) tumors had superior PFS (median 14 months, 95% CI 2-21) compared to patients with any other (KLF12 + TFF2) expression profile (median PFS 7 months, 95% CI 5-10, p = 0.021). The Hazard Ratio for disease progression for (KLF12-high + TFF2-low) versus any other KLF12 + TFF2 expression profile was 2.92 (p = 0.03) in the Validation and 1.29 (p = 0.39) in the Control set.


Our «three-stage» hypothesis-generating study failed to validate the prognostic significance of a five-gene classifier in mCRC patients. Exploratory analyses suggest two gene signatures that are potentially associated with bevazicumab benefit in patients with advanced colon cancer.

Peer Review reports


The high cost of targeted therapies as well as their conceptual definition as «targeting» specific molecular aberrations mandate the use of biomarkers in modern oncology practice [1]. Biomarkers are tumor and host characteristics that either define the natural course of a malignancy irrespective of therapy (prognostic) or the probability of patient benefit from a therapy administered (predictive) [2]. Although both are clinically relevant, less progress has been made in the field of the latter.

Angiogenesis is the process of new blood vessel formation and is pivotal for tumor growth, invasion and metastases [3]. Bevacizumab, a humanized monoclonal antibody that binds and neutralizes one of the main effectors of malignant angiogenesis, the Vascular Endothelial Growth Factor (VEGF), has been licensed for the treatment of patients with metastatic colorectal cancer combined with chemotherapy [4]. However, the rather modest improvement in response and survival outcomes achieved indicate the rich tumor heterogeneity and the probability that only a subset of tumors are amenable to VEGF modulation. The molecular characterization of tumors responsive to bevacizumab remains the Holy Grail for a worlwide community of investigators.

Genomic technologies are being widely used to study tumors at the molecular level. Since the extraction of RNA from formalin-fixed, paraffin-embedded (FFPE) tumor tissue has been optimized, microarray-based multigene expression profiling platforms have been developed for the identification of molecular signatures associated with various tumor characteristics [5]. The larger scale availabilty and more straightforward feasibility of performing quantitative PCR (qPCR) assays commonly led to attempts to adapt microarray signatures to qPCR methodologies.

In this study, we used a microarray platform to profile the expression of 24526 genes in a test set of 16 patients with metastatic colon cancer treated with bevacizumab, aiming to identify a select set of genes associated with superior outcome on bevacizumab. We then studied the expression of these genes using qPCR in an independent set of patients who received bevacizumab and in a control set of patients who were treated with chemotherapy only, in order to confirm their significance and to dissect their potential predictive from prognostic utility.


Patients with chemonaive metastatic colon cancer who received first-line standardized chemotherapy protocols with or without bevacizumab between 2005 and 2009 in oncology centers affiliated with the Hellenic Cooperative Oncology Group (HeCOG), consented for the research use of their biologic material. FFPE blocks were fully annotated with clinicopathologic characteristics. The translational research protocol was approved by the Scientific Commitee, Papageorgiou Hospital, Thessaloniki, 185/8-10-2013. The patient sets consisted of three cohorts: a) the Test set (N = 16, patients treated with chemotherapy and bevacizumab) in which FFPE microarray analysis was performed in order to identify candidate genes predictive of bevacizumab benefit, b) the Bevacizumab qPCR set (N = 49, an independent cohort of patients treated with chemotherapy and bevacizumab) and the Control qPCR set (N = 72, patients treated with chemotherapy without bevacizumab). All patients had a performance status of 0-1. In the latter two independent sets, expression of the selected genes with possible predictive/prognostic significance was quantified by means of qPCR (Additional file 1: Figure S1 for REMARK diagram).

For RNA extraction from FFPE tumors, H&E sections were histologically reviewed and areas containing >50% tumor cells were marked; these were macrodissected from serial unstained sections at 8um after deparaffinization and submitted for RNA extraction with the RNeasy FFPE Kit (Qiagen, Hilden, D). For the Test Set, two series of RNA samples were prepared; one of these was submitted for Illumina profiling. RNA samples from the Test, Bevacizumab qPCR and the Control qPCR Sets were processed for reverse transcription and first strand cDNA synthesis with the Superscript III and random hexamers (Invitrogen/Life Technologies). All reagents and systems were used according to the instructions of the manufacturers. cDNAs were normalized at 25 ng/ul and stored at -20°C until use.

Microarray methodology and analysis

We performed global gene expression profiling on the 16 patients of our test set using whole genome DASL (cDNA-mediated, Annealing, Selection, and Ligation) arrays (Illumina, CA), covering more than 24,000 transcripts. This technology overcomes the challenges of profiling partially degraded RNA, often extracted from FFPE samples and provides high-quality gene expression data [6]. We isolated 250 ng of total RNA in a concentration of 25 ng/μl, as required by Expression Analysis Inc (Durham, NC). The A260/A280 ratio of each RNA specimen exceeded 1.6. Outlier exclusion was based on the percent present call of the samples; detection rate >12000 transcripts. Microarray experiments were carried out at Expression Analysis Inc (Durham, NC) according to the manufacturer’s recommendations. The microarray data have been submitted to Gene Expression Omnibus as study GSE53127 and can be viewed at:

5-gene predictor validation with qPCR

The 5-gene predictor was evaluated on cDNA samples from the Bevacizumab qPCR and Control sets with qPCR and hydrolysis probes (TaqMan® MGB probes, Applied Biosystems/Life Technologies). The following premade assays were selected for amplicons matching the regions targeted by corresponding Illumina probes (assay ID, NM-reference, exon spanning, location, size): AGR2 Hs00356521_m1 (NM_006408.3, ex 7-8, 665, 69 bp); ALDH6A1 Hs00194421_m1 (NM_005589.2, ex 11-12, 1607, 131 bp); KLF12 Hs00273134_m1 (NM_007249.4, ex 6-7, 1089, 100 bp); MCM5 Hs01052142_m1 (NM_006739.3, ex 16-17, 2197, 70 bp); and, TFF2 Hs00989207_m1 (NM_005423.4, ex 3-4, 520, 68 bp). For normalization and relative expression assessment, 3 premade TaqMan® MGB assays for endogenous control transcripts were used: #4333767 F for GUSB; Hs00183533_m1 for IPO8; and Hs00427620_m1 for TBP. Samples were run in duplicates, in 10ul reactions (2ul cDNA template per reaction) in an ABI7900HT real time PCR system under default conditions. A commercially available reference RNA derived from multiple transformed cell lines (TaqMan® Control Total RNA, cat. no 4307281, Applied Biosystems) was applied in multiple positions in each run as positive control and for inter-run evaluation of PCR assay efficiency. No-template controls were also included. Samples were run in duplicates, at least in two metachronous runs. To obtain linear Relative Quantification (RQ) values, relative expression was assessed as (40-dCT), as previously described, whereby dCT (or deltaCT) was calculated as (average target CT) – (average endogenous control CT) from all eligible measurements [7]. Samples were considered eligible for analysis when (a) both endogenous control CTs in duplicates were <36 and when duplicate dCT’s for the same sample within the same run were <0.75. The efficiency of all assays was considered as comparable, since the difference between inter-run RQ values for the reference RNA sample was <1 for all assays. Upon testing for target RQ value compliance per sample with each endogenous control, TPB yielded the most unstable results and was thus not included in the final assessment of RQ values.

Statistical analysis

Analyses of the microarray data were performed using BRB-ArrayTools Software developed by Dr. Richard Simon and BRB-ArrayTools Development Team [8]. After quantile normalization of the samples, we excluded one fourth of the genes showing minimal variation across our dataset. In order to assess gene expression profiles predictive of bevacizumab benefit, we utilized Compound Covariate Predictor, Diagonal Linear Discriminant Analysis, Nearest Neighbor Classification, Support Vector Machines with linear kernel and Bayesian Compound Covariate Predictor. These algorithms incorporate genes differentially expressed among different classes as assessed by the random variance t-test. Evaluation of the predictive value of these methods was based on Leave-One-Out-Cross-Validation. The 8-month Progression-Free status was used as endpoint and surrogate marker of bevacizumab benefit.

For all the markers the median (50th percentile) were examined as possible threshold for prognostic significance categorizing the gene expression levels into high versus low. The expression of five genes was examined for correlation to the following parameters as endpoints: a) Objective response rate (ORR- Best response to therapy; complete or partial response), prior to any metastasectomy, b) Progression-free survival (PFS), calculated from the initiation of first-line therapy to disease progression, death or last follow-up, whichever occurs first.

ORR was chosen as an objective, easy to measure endpoint which is not confounded by the potentially curative resection of metastases in some patients. On the other hand, PFS was used as a survival endpoint in order to investigate the possible association of genes with survival, without impact on tumor regression rates.

The Fisher’s exact test was used to examine possible associations between gene expressions with the overall response rate (ORR), while odds ratios were also calculated in order to measure the association. Time-to-event distributions were estimated using the Kaplan-Meier curves. The log-rank test and Cox’s proportional hazards models was used to examine the univariate prognostic significance of the markers for PFS. For all univariate tests the significance level was set at α = 0.05. Multivariate analysis included clinical parameters and gene expression profiles that were significant in the univariate setting. The SAS software was used for statistical analysis (SAS for Windows, version 9.3, SAS Institute Inc., Cary, NC, USA).


Microarray identification of candidate genes in the Test set

The sixteen patients in the Test set had metastatic colon cancer treated with FOLFIRI + bevacizumab (n = 13) or CapOx + bevacizumab (n = 3). Three patients had achieved a complete remission and nine a partial remission (objective response rate, ORR 75.0%). Ten patients (62.6%) remained progression-free for at least eight months. Patient and tumor characteristics are shown in Table 1. The 8-month Progression-Free status was used as the endpoint for examining the association of genome-wide gene expression with it as a surrogate marker of bevacizumab activity. In this test set, we developed gene expression models using different algorithms to predict which patients would progress within 8 months. The optimal predictor comprised of five genes, differentially expressed between patients who had progressed within 8 months and those who remained progression free at that time point, at a significance level of 0.00005 (random variance t-test). Prediction accuracies fluctuated between 89–94%, based on the five different algorithms. Tumor tissue samples from patients with unfavorable PFS status were found to overexpress four out of the five genes (AGR2, ALDH6A1, TFF2, MCM5) and underexpress KLF12. Information on the five-gene model can be found in Table 1 and Figure 1.

Table 1 Patient demographics and Gene selection in Test set
Figure 1
figure 1

Gene clustering in microarray analysis of test set according to 8-month progression-free endpoint. Hierarchical clustering was performed using the euclidean distance and the average linkage algorithm. UPPER ROW-PD8m: Progressive disease during first 8 months from treatment start, noPD8m: Progression-free during 8 months from treatment start. LOWER ROW- Patient tumor samples (n = 16, two patients ineligible). VERTICAL LEFT: Expression of five genes. Red denotes overexpression, Green denotes underexpression.

Application of the 5-gene predictor with qPCR in the Bevacizumab qPCR and Control sets

Application of the 5 genes of the predictor in the Test set yielded significant associations with PFS in the predicted direction for individual markers AGR2, ALDH6A1, MCM5 and TFF2 (favorable PFS for patients with tumors expressing low RQ values, with low at median cut-off) but none for KLF12 (Additional file 2: Figure S2). Predetermined clustering of RQ values for the same genes according to the predictor pattern was not possible for the 16 tumors in this set, since no groups applicable for statistics could be formed. Therefore, the Test set could not be used for the evaluation of the predictor with qPCR.


In the Bevacizumab qPCR set, forty-nine patients with metastatic colon cancer, of a median age of 60.1 years, received mostly irinotecan- or oxaliplatin-based regimens combined with bevacizumab as first-line therapy. An objective response was seen in 36.7% of them, while 51.0% were free from progression for at least 8 months. Seven patients (14.3%) underwent potentially curative metastasectomies following therapy-induced cytoreduction. In the control set, seventy-two advanced colon cancer patients of a median age of 64.5 years received mostly irinotecan-based chemotherapy (58.3%) or oxaliplatin-based regimens as first-line therapy. No bevacizumab was administered. An objective response was seen in 43.1% of them, while 54.2% were free from progression for at least 8 months. Six patients (8.3%) underwent potentially curative metastasectomies following therapy-induced cytoreduction. Most characteristics of the Bevacizumab qPCR and Control sets were matched and are shown in Table 2.

Table 2 Patient demographics in Bevacizumab qPCR and Control sets


The expression of five selected genes individually as well as the complex expression profiles of gene combinations were studied for association with objective response rate. Only the expression profile of ALDH6A1 + TFF2 + MCM5 correlated strongly with response to bevacizumab (see Table 3). Among patients harboring tumors with (ALDH6A1 + TFF2 + MCM5)-all low gene expression, 85.7% responded to bevacizumab chemoimmunotherapy in the Bevacizumab qPCR set versus only 28.6% of those with any other (ALDH6A1 + TFF2 + MCM5) expression profile (p = 0.007). In the Control set, only 36.4% of patients harboring (ALDH6A1 + TFF2 + MCM5)-all low tumors responded to chemotherapy-only, compared to a rather similar ORR of 44.3% in patients with any other (ALDH6A1 + TFF2 + MCM5) expression profile (p = 0.747). The Odds Ratio for Response in patients with (ALDH6A1 + TFF2 + MCM5)-all low tumors compared to any other (ALDH6A1 + TFF2 + MCM5) profile was 0.72 (p = 0.63) in the Control set, but was 15.00 (p = 0.0168) in the Bevacizumab qPCR set. The distribution of patients with tumor regression, stable and progressive disease according to the (ALDH6A1 + TFF2 + MCM5) expression profile is visualized in Figure 2 for the Bevacizumab qPCR set. Of note, low expression of the MCM5 gene also correlated significantly, though less strongly, with ORR (MCM5-low tumors: Bevacizumab qPCR set, ORR 52.0%, p = 0.038 - Control set, ORR 38.9%, p = 0.63).

Table 3 Gene expession associated with patient outcomes
Figure 2
figure 2

Objective response rate by (ALDH6A1 + TFF2 + MCM5) gene expression in Bevacizumab qPCR set.


When we examined the association of expression of our five genes with PFS, a different picture emerged. In the Bevacizumab qPCR set, patients with (KLF12-high + TFF2-low) tumors had superior PFS (median 14 months, 95% CI 2-21) compared to patients with any other (KLF12 + TFF2) expression profile (median PFS 7 months, 95% CI 5-10, p = 0.021, Figure 3). On the contrary, patients not treated with bevacizumab in the Control set had a median PFS of 11 months (95% CI 8-13) when harboring (KLF12-high + TFF2-low) tumors, not significantly different from the median PFS of 8 months (95% CI 7-10) of patients with any other (KLF12 + TFF2) expression profile (p = 0.38). The Hazard Ratio for risk of progression for any other (KLF12 + TFF2) expression profile versus the reference category of (KLF12-high + TFF2-low) tumors was 2.92 (p = 0.03) in the presence of bevacizumab and 1.29 (p = 0.39) in the absence of bevacizumab.

Figure 3
figure 3

Progression-free survival by KLF12 + TFF2 gene expression. a) Control set. b) Bevacizumab qPCR set.

Of note, the profile of (ALDH6A1 + TFF2 + MCM5) expression did not show any significant association with PFS, while that of (KLF12 + TFF2) no association with tumor response (ORR). Finally, the profile of all-low (ALDH6A1 + TFF2) tumor gene expression showed a marginally significant association with superior survival in the chemotherapy-only Control set, but not in the Bevacizumab qPCR set (Table 3). Detailed information on all examined qPCR gene expression profiles in all sets are given in Additional file 3: Table S1 and Additional file 4: Table S2.

Multivariate analysis

A multivariate analysis model incorporating performance status, gender, primary site in the colon, type of first line chemotherapy, age, occurrence of metastasectomy and gene expression profiles significant in univariate analyses was applied in the total population of Bevacizumab qPCR and Control sets. No single gene or complex multigene tumor expression profile showed any independent significant association with either ORR or PFS. Parameters with independent prognostic significance for PFS were poor performance status (PS 2-3 vs 0-1, HR for disease progression 18.2, p = 0.0004) and occurrence of metastasectomy (HR for disease progression 0.04, p = 0.0015).


Our microarray-based exploratory study identified the expression of five genes with significant correlation to the probability of disease control from bevacizumab therapy beyond the 8-month benchmark duration. KLF12 (Kruppel-like factor 12) is located at 13q22, encoding a member of the Kruppel-like zinc finger protein family of transcription factors. It represses the expression of the Activator protein-2 alpha (AP-2 alpha) gene, an important regulator of vertebrate development and carcinogenesis, by binding its promoter [9]. As a transcription factor, overexpression of KLF12 in endometrial cancer cell lines significantly repressed proliferation and secretion of pro-survival factors such as insulin-like growth factor binding protein-1 [10]. On the other hand, KLF12 was shown to induce cell proliferation, angiogenesis and invasion in gastric cancer cell lines and clinical samples [11]. High KLF12 gene expession may constitute a surrogate marker of tumors with high proliferative, angiogenic and migratory potential, amenable to VEGF blockade.

TFF2 (Trefoil Factor 2) is located at 21q22.3 and encodes a stable secretory protein, member of the Trefoil family, expressed in gastrointestinal mucosa. Their functions are not defined, but they may protect the mucosa from insults [12]. Breast, pancreas and bile duct cancer cell line experiments suggested that TFF2 expression induces cell migration via Platelet-Activation Receptor 4 (PAR4) and Panc1 activation, as well as mitosis via EGFR/MAPK axis signaling [1317]. On the other hand, TFF2 was shown to possess anti-inflammatory properties and to undergo promoter methylation during gastric cancer progression, data pointing to a tumor-suppressive function [18, 19]. ALDH6A1 (aldehyde dehydrogenase 6 family, member A1) was mapped at 14q24.3. The encoded protein is a mitochondrial methylmalonate semialdehyde dehydrogenase that plays a role in the valine and pyrimidine catabolic pathways. This protein catalyzes the irreversible oxidative decarboxylation of malonate and methylmalonate semialdehydes to acetyl- and propionyl-CoA [20, 21]. Despite its regulatory role in mitochondrial energy production and DNA catabolism, no studies examined its putative contribution to cancer homeostasis to date.

MCM5 (minichromosome maintenance complex component 5) at 22q13.1 encodes for a member of the MCM family of chromatin-binding proteins that stimulates cell transition from G0 to G1/S phase of the cell cycle and actively participates in cell cycle regulation [22, 23]. Data from clinical and preclinical models of skin, esophageal, bladder and gastrointestinal carcinomas further confirm the proliferative, migratory and cell cycle activating properties of the MCM5 protein [2427]. Low MCM5 gene expression could mark tumors with low proliferation, abnormal vasculature and hypoxia, the profile most amenable to vessel normalization and cell kill by bevacizumab + chemotherapy. Of note, the most well studied pro-oncogenic gene, AGR2 (Anterior Gradient 2, at 7p21.3) which has identified oncogenic functions such as attenuation of endoplasmic reticulum stress, transition from G0 to G1 phase, inhibition of cell senescence and association with tumor stage, was not found to correlate with either response or progression-free survival in the validation control [28].

We selected a «three-stage» design for our experiment, which should be viewed as hypothesis-generating, rather than proof of principle. First, we used data mining in order to identify genes with potential association with bevacizumab benefit from a test set of 16 patients for whom genome-wide gene expression was studied in a microarray platform. Second, we tested the predictive performance of the 5 genes in the microarray predictor with qPCR, a method more convenient and realistic for clinical practice, in the test set and in an independent Bevacizumab qPCR set of 49 patients who had been treated with bevacizumab. Third, we used the same qPCR approach in order to examine the prognostic significance of the predictor genes in a matched control set of 72 patients who received first-line chemotherapy, but not bevacizumab. We chose to include patients who had metastasectomy after chemotherapy in both cohorts, despite introducing a positive bias: surgical resection of metastases could alter the natural course of disease and be a confounding factor in our search for a predictor of bevacizumab benefit. However, excluding metastasectomy cases would introduce a negative bias, as it is likely that patients led to potentially curative metastasectomy would be the ones with major cytoreduction and disease control from chemotherapy + bevacizumab. Accordingly, we used two metrics for clinical benefit: best tumor response before metastasectomy, which is not confounded by the latter, and progression-free survival, which is more sensitive than overall survival but potentially influenced by metastasectomy. Of note, the incidence of metastasectomy was not significantly different in the Bevacizumab qPCR and Control cohorts. We failed to identify any qPCR gene signature that could incorporate all five genes identified in the Test set, consequently our Bevacizumab qPCR cohort should not be viewed as a «validation» cohort but as an exploratory cohort for the study of a new qPCR signature consisting of some of the preselected five genes.

The conflicting function of genes studied reported by other investigators maybe due to differences in cancer types, tumor microenviroment, expression of multiple other modulating biomolecules, disease stage as well as study design and experimental methodologies. They constitute the interpretation of observed associations of genes with bevacizumab activity extremely difficult, especially in view of lack of constistency when various benefit metrics are examined (response rate, PFS) and the absence of independent significance in multivariate analyses. The inability of our response-predictive qPCR profiles to impact on PFS and vice versa surely raises concerns about the validity of our findings. Still, the decoupling of ORR from PFS benefit may be due to the impact of metastasectomies on PFS [29]. It could also be explained by the commonly reported discrete biologies underlying the phenomena of tumor regression and of control of tumor proliferation, invasion and virulence, especially when anti-angiogenic therapies are administered [30, 31]. Technical differences regarding RNA (microarrays) and mRNA (qPCR) sequence detection could also account for the inability to successfully recapitulate the entire 5-gene microarray predictor with qPCR. In view of the multiple analyses performed (see Additional file 3: Table S1 and Additional file 4: Table S2), it is not safe to conclude that the association of our qPCR profiles with clinical benefit from bevacizumab may be reflecting important functional roles of these genes or establish them as surrogate markers of genetic subsets of tumors responsive to anti-angiogenesis. These associations could simply constitute random findings and the data should only be viewed as hypothesis-generating.

Although the search for a validated gene signature predictive for bevacizumab benefit did not bear fruit, some findings are consistently reported. Gene expression profiling studies in bevacizumab-treated patients with glioblastomas, breast and colon cancer identified predictive gene signatures with little or no gene overlap which possessed however a common repertoire of the gene functional ontologies implicated: cell proliferation, mitochondrial energy production, lipid metabolism, migration/invasion, hypoxia regulation and immune response [3235]. The genes identified here are also characterized by the functional roles above. Brauer et al suggested that a genetic profile predictive for benefit from anti-angiogenesis may be independent of tumor primary, while Fiebig at al could not assign a known function in 59% of the 35 genes predicting for bevacizumab benefit in colorectal cancer xenografts and clinical samples [33, 34]. Hu et al reported that although a change in multigene expession in 21 bevacizumab-treated glioblastoma patients correlated to outcome, they could not identify a baseline gene expression signature with prognostic significance. In our case, only gene expression data at baseline were available [36].


To conclude, we identified two distinct qPCR gene expression profiles correlating with response (low expression of ALDH6A1 + TFF2 + MCM5) and with PFS (KLF12-high + TFF2-low) in advanced colon cancer patients managed with bevacizumab and chemotherapy. Despite our three-cohort experimental design, the moderate sample size, the plethora of variables under study and of analyses performed preclude us from establishing a predictive utillity for these genetic profiles before further validation in independent cohorts.


  1. Blay JY, Lacombe D, Meunier F, Stupp R: Personalised medicine in oncology: questions for the next 20 years. Lancet Oncol. 2012, 13: 448-449. 10.1016/S1470-2045(12)70156-0.

    Article  PubMed  Google Scholar 

  2. Heckman-Stoddard BM: Oncology biomarkers: discovery, validation, and clinical use. Semin Oncol Nurs. 2012, 28: 93-98. 10.1016/j.soncn.2012.03.003.

    Article  PubMed  Google Scholar 

  3. Carmeliet P, Jain RK: Angiogenesis in cancer and other diseases. Nature. 2000, 407: 249-257. 10.1038/35025220.

    Article  CAS  PubMed  Google Scholar 

  4. Hayes DF: Bevacizumab treatment for solid tumors: boon or bust?. JAMA. 2011, 305: 506-508. 10.1001/jama.2011.57.

    Article  CAS  PubMed  Google Scholar 

  5. Khosravi Shahi P, Perez MG: The application of Microarrays technology in clinical oncology. An Med Interna. 2006, 23: 255-256.

    CAS  PubMed  Google Scholar 

  6. April C, Klotzle B, Royce T, et al: Whole-genome gene expression profiling of formalin-fixed, paraffin-embedded tissue samples. PLoS One. 2009, 4: e8162-10.1371/journal.pone.0008162.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Hennig G, Gehrmann M, Stropp U, et al: Automated extraction of DNA and RNA from a single formalin-fixed paraffin-embedded tissue section for analysis of both single-nucleotide polymorphisms and mRNA expression. Clin Chem. 2010, 56: 1845-1853. 10.1373/clinchem.2010.151233.

    Article  CAS  PubMed  Google Scholar 

  8. Simon R, Lam A, Li MC, et al: Analysis of gene expression data using BRB-ArrayTools. Cancer Inform. 2007, 3: 11-17.

    PubMed  PubMed Central  Google Scholar 

  9. Zhu CH, Huang Y, Broman MT, Domann FE: Expression of AP-2 alpha in SV40 immortalized human lung fibroblasts is associated with a distinct pattern of cytosine methylation in the AP-2 alpha promoter. Biochim Biophys Acta. 2001, 1519: 85-91. 10.1016/S0167-4781(01)00227-5.

    Article  CAS  PubMed  Google Scholar 

  10. Shen X, Hu Y, Jiang Y, et al: Kruppel-like factor 12 negatively regulates human endometrial stromal cell decidualization. Biochem Biophys Res Commun. 2013, 433: 11-17. 10.1016/j.bbrc.2013.02.078.

    Article  CAS  PubMed  Google Scholar 

  11. Nakamura Y, Migita T, Hosoda F, et al: Kruppel-like factor 12 plays a significant role in poorly differentiated gastric cancer progression. Int J Cancer. 2009, 125: 1859-1867. 10.1002/ijc.24538.

    Article  CAS  PubMed  Google Scholar 

  12. Hanisch FG, Ragge H, Kalinski T, et al: Human gastric TFF2 peptide contains an N-linked fucosylated N, N’-diacetyllactosediamine (LacdiNAc) oligosaccharide. Glycobiology. 2013, 23: 2-11. 10.1093/glycob/cws131.

    Article  CAS  PubMed  Google Scholar 

  13. May FE, Semple JI, Prest SJ, Westley BR: Expression and motogenic activity of TFF2 in human breast cancer cells. Peptides. 2004, 25: 865-872. 10.1016/j.peptides.2003.12.024.

    Article  CAS  PubMed  Google Scholar 

  14. Guppy NJ, El-Bahrawy ME, Kocher HM, et al: Trefoil factor family peptides in normal and diseased human pancreas. Pancreas. 2012, 41: 888-896. 10.1097/MPA.0b013e31823c9ec5.

    Article  CAS  PubMed  Google Scholar 

  15. Kamlua S, Patrakitkomjorn S, Jearanaikoon P, et al: A novel TFF2 splice variant (EX2TFF2) correlates with longer overall survival time in cholangiocarcinoma. Oncol Rep. 2012, 27: 1207-1212.

    CAS  PubMed  Google Scholar 

  16. Zhang Y, Yu G, Wang Y, et al: Activation of protease-activated receptor (PAR) 1 by frog trefoil factor (TFF) 2 and PAR4 by human TFF2. Cell Mol Life Sci. 2011, 68: 3771-3780. 10.1007/s00018-011-0678-6.

    Article  CAS  PubMed  Google Scholar 

  17. Kosriwong K, Menheniott TR, Giraud AS, et al: Trefoil factors: tumor progression markers and mitogens via EGFR/MAPK activation in cholangiocarcinoma. World J Gastroenterol. 2011, 17: 1631-1641. 10.3748/wjg.v17.i12.1631.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Peterson AJ, Menheniott TR, O’Connor L, et al: Helicobacter pylori infection promotes methylation and silencing of trefoil factor 2, leading to gastric tumor development in mice and humans. Gastroenterology. 2010, 139: 2005-2017. 10.1053/j.gastro.2010.08.043.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Wills-Karp M, Rani R, Dienger K, et al: Trefoil factor 2 rapidly induces interleukin 33 to promote type 2 immunity during allergic asthma and hookworm infection. J Exp Med. 2012, 209: 607-622. 10.1084/jem.20110079.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Sass JO, Walter M, Shield JP, et al: 3-Hydroxyisobutyrate aciduria and mutations in the ALDH6A1 gene coding for methylmalonate semialdehyde dehydrogenase. J Inherit Metab Dis. 2012, 35: 437-442. 10.1007/s10545-011-9381-x.

    Article  CAS  PubMed  Google Scholar 

  21. Chambliss KL, Gray RG, Rylance G, et al: Molecular characterization of methylmalonate semialdehyde dehydrogenase deficiency. J Inherit Metab Dis. 2000, 23: 497-504. 10.1023/A:1005616315087.

    Article  CAS  PubMed  Google Scholar 

  22. Snyder M, He W, Zhang JJ: The DNA replication factor MCM5 is essential for Stat1-mediated transcriptional activation. Proc Natl Acad Sci USA. 2005, 102: 14539-14544. 10.1073/pnas.0507479102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Snyder M, Huang XY, Zhang JJ: The minichromosome maintenance proteins 2-7 (MCM2-7) are necessary for RNA polymerase II (Pol II)-mediated transcription. J Biol Chem. 2009, 284: 13466-13472. 10.1074/jbc.M809471200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Liu H, Takeuchi S, Moroi Y, et al: Expression of minichromosome maintenance 5 protein in proliferative and malignant skin diseases. Int J Dermatol. 2007, 46: 1171-1176. 10.1111/j.1365-4632.2007.03335.x.

    Article  CAS  PubMed  Google Scholar 

  25. Going JJ, Keith WN, Neilson L, et al: Aberrant expression of minichromosome maintenance proteins 2 and 5, and Ki-67 in dysplastic squamous oesophageal epithelium and Barrett’s mucosa. Gut. 2002, 50: 373-377. 10.1136/gut.50.3.373.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Giaginis C, Giagini A, Tsourouflis G, et al: MCM-2 and MCM-5 expression in gastric adenocarcinoma: clinical significance and comparison with Ki-67 proliferative marker. Dig Dis Sci. 2011, 56: 777-785. 10.1007/s10620-010-1348-5.

    Article  CAS  PubMed  Google Scholar 

  27. Brems-Eskildsen AS, Zieger K, Toldbod H, et al: Prediction and diagnosis of bladder cancer recurrence based on urinary content of hTERT, SENP1, PPP1CA, and MCM5 transcripts. BMC Cancer. 2010, 10: 646-10.1186/1471-2407-10-646.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Chevet E, Fessart D, Delom F, et al: Emerging roles for the pro-oncogenic anterior gradient-2 in cancer development. Oncogene. 2013, 32: 2499-2509. 10.1038/onc.2012.346.

    Article  CAS  PubMed  Google Scholar 

  29. Adam R: The importance of visceral metastasectomy in colorectal cancer. Ann Oncol. 2000, 11 (Suppl 3): 29-36. 10.1093/annonc/11.suppl_3.29.

    Article  PubMed  Google Scholar 

  30. Valera VA, Walter BA, Yokoyama N, et al: Prognostic groups in colorectal carcinoma patients based on tumor cell proliferation and classification and regression tree (CART) survival analysis. Ann Surg Oncol. 2007, 14: 34-40.

    Article  PubMed  Google Scholar 

  31. Smaglo BG, Hwang J: Continuum of care with anti-angiogenic therapies in metastatic colorectal cancer. J Gastrointest Oncol. 2013, 4: 299-307.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Piao Y, Liang J, Holmes L, et al: Acquired Resistance to Anti-VEGF Therapy in Glioblastoma Is Associated with a Mesenchymal Transition. Clin Cancer Res. 2013, 19: 4392-4403. 10.1158/1078-0432.CCR-12-1557.

    Article  CAS  PubMed  Google Scholar 

  33. Fiebig HH, Vuaroqueaux V, Korrat A, et al: Predictive gene signatures for bevacizumab and cetuximab as well as cytotoxic agents. Int J Clin Pharmacol Ther. 2012, 50: 70-71. 10.5414/CPP50070.

    Article  CAS  PubMed  Google Scholar 

  34. Brauer MJ, Zhuang G, Schmidt M, et al: Identification and analysis of in vivo VEGF downstream markers link VEGF pathway activity with efficacy of anti-VEGF therapies. Clin Cancer Res. 2013, 19: 3681-3692. 10.1158/1078-0432.CCR-12-3635.

    Article  CAS  PubMed  Google Scholar 

  35. Yang SX, Steinberg SM, Nguyen D, et al: Gene expression profile and angiogenic marker correlates with response to neoadjuvant bevacizumab followed by bevacizumab plus chemotherapy in breast cancer. Clin Cancer Res. 2008, 14: 5893-5899. 10.1158/1078-0432.CCR-07-4762.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Hu YL, DeLay M, Jahangiri A, et al: Hypoxia-induced autophagy promotes tumor cell survival and adaptation to antiangiogenic treatment in glioblastoma. Cancer Res. 2012, 72: 1773-1783. 10.1158/0008-5472.CAN-11-3831.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references


We would like to acknowledge Mrs Zoi Alexopoulou and Mrs Maria Moschoni for their contribution to data analysis (ZA) and logistical support (MM).

Author information

Authors and Affiliations


Corresponding author

Correspondence to George Pentheroudakis.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contribution

GP contributed to conception and coordination of the study, analysis and interpretation of data, writing of the manuscript. VK contributed to design of RT-PCR experiments, conception of the study, analysis and interpretation of data. EF contributed to design of microarray experiments, analysis and interpretation of data. GK contributed to to analysis and interpretation of data. GF contributed to conception of the study, analysis and interpretation of data. All other authors contributed to analysis and interpretation of data. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Figure S1: REMARK diagram. (DOC 27 KB)


Additional file 2: Figure S2: Performance of all 5 genes of the microarray predictor in the test, validation and control sets with respect to patient PFS. AGR2, ALDH6A1 and MCM2 were consistent with the predictor in the test set but not in the validation set; instead, 2 out of 3 genes were associated with longer PFS in the non-bevacizumab treated cohort. TFF2 was consistent with the predictor in the test and validation sets. KLF12 was the only gene that could not be validated in the test set with qPCR; high transcript levels of this gene were, however, showed a trend for better outcome in the validation set. Based on these findings and upon failure to transfer the entire 5-gene signature into a single qPCR profile, KLF12 and TFF2 RQ values were profiled for assessing their possible value in predicting PFS upon bevacizumab treatment. (PPT 2 MB)


Additional file 3: Table S1: Univariate Cox regression for each qPCR gene expression and their combinations among dataset groups in terms of PFS. (DOC 56 KB)


Additional file 4: Table S2: Fisher’s exact test for each qPCR gene expression and their combinations among dataset groups in terms of ORR. (DOC 178 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Pentheroudakis, G., Kotoula, V., Fountzilas, E. et al. A study of gene expression markers for predictive significance for bevacizumab benefit in patients with metastatic colon cancer: a translational research study of the Hellenic Cooperative Oncology Group (HeCOG). BMC Cancer 14, 111 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: