Recent studies suggest CD133, a surface protein widely used for isolation of colon cancer stem cells, to be associated with tumor angiogenesis and recurrence. We hypothesized that gene expression levels and germline variations in CD133 will predict clinical outcome in patients with metastatic colorectal cancer (mCRC), treated in first-line setting with 5-fluorouracil, oxaliplatin and bevacizumab (BV), and we investigated whether there is a correlation with gene expression levels of CD133, vascular endothelial growth factor (VEGF) and its receptors. We evaluated intra-tumoral gene expression levels by quantitative real-time (RT) PCR from 54 patients and three germline variants of the CD133 gene by PCR-restriction-fragment length polymorphism from 91 patients with genomic DNA. High gene expression levels of CD133 (>7.76) conferred a significantly greater tumor response (RR=86%) than patients with low expression levels (7.76, RR=38%, adjusted P=0.003), independent of VEGF or its receptor gene expression levels. Gene expression levels of CD133 were significantly associated with VEGF and its receptors messenger RNA levels (VEGFR-1 (P<0.01), -2 and -3, P<0.05). Combined analyses of two polymorphisms showed a significant association with progression-free survival (PFS) (18.5 months vs 9.8 months, P=0.004) in a multivariate analysis as an independent prognostic factor for PFS (adjusted P=0.002). These results suggest that CD133 is a predictive marker for standard first-line BV-based treatment in mCRC.