We present a model of hepatic colorectal metastases which represents monoclonal cell lines double-labeled by luciferase and tdTomato. These cells form liver metastasis in varying numbers and patterns similar to those observed in patients. Using in vivo and ex vivo luminescent and fluorescent imaging we determine the growth kinetics and clonogenic frequency of tumor cells colonizing liver. Molecular profiling detected stable expressional differences between clones consistent with their phenotypes. The data indicate that clinically relevant phenotypes of liver metastases can be modeled in vivo.