Intra-Tumoral Heterogeneity in Metastatic Potential and Survival Signaling between Iso-Clonal HCT116 and HCT116b Human Colon Carcinoma Cell Lines

Background Colorectal cancer (CRC) metastasis is a leading cause of cancer-related deaths in the United States. The molecular mechanisms underlying this complex, multi-step pathway are yet to be completely elucidated. Recent reports have stressed the importance of intra-tumoral heterogeneity in the development of a metastatic phenotype. The purpose of this study was to characterize the intra-tumoral phenotypic heterogeneity between two iso-clonal human colon cancer sublines HCT116 and HCT116b on their ability to undergo metastatic colonization and survive under growth factor deprivation stress (GFDS). Materials and Methods HCT116 and HCT116b cells were transfected with green fluorescence protein and subcutaneously injected into BALB/c nude male mice. Once xenografts were established, they were excised and orthotopically implanted into other male BALB/c nude mice using microsurgical techniques. Animal tissues were studied for metastases using histochemical techniques. Microarray analysis was performed to generate gene signatures associated with each subline. In vitro assessment of growth factor signaling pathway was performed under GFDS for 3 and 5 days. Results Both HCT116 and HCT116b iso-clonal variants demonstrated 100% primary tumor growth, invasion and peritoneal spread. However, HCT116 was highly metastatic with 68% metastasis observed in liver and/or lungs compared to 4% in HCT116b. Microarray analysis revealed an upregulation of survival and metastatic genes in HCT116 cells compared to HCT116b cells. In vitro analysis showed that HCT116 upregulated survival and migratory signaling proteins and downregulated apoptotic agents under GFDS. However, HCT116b cells effectively showed the opposite response under stress inducing cell death. Conclusions We demonstrate the importance of clonal variation in determining metastatic potential of colorectal cancer cells using the HCT116/HCT116b iso-clonal variants in an orthotopic metastatic mouse model. Determination of clonal heterogeneity in patient tumors can serve as useful tools to identify clinically relevant biomarkers for diagnostic and therapeutic assessment of metastatic colorectal cancer.


Introduction
Colorectal cancer (CRC) is a major contributor of cancerrelated deaths in the United States [1].Metastasis to distant organ sites significantly affects the mortality rate resulting from the disease [1,2].The Cancer Genome Atlas Network has recently reported the multi-dimensional genomic changes associated with CRC with the goal to provide deeper insight into the pathophysiology of CRC to identify potential therapeutic targets [3].Identification and characterization of novel molecular targets for CRC metastasis is a pressing need since to date there are no effective anti-metastatic therapies available.Tumor cells display remarkable clonal heterogeneity due to both genetic and nongenetic influences [4].Several clinically important phenotypes including the ability to undergo metastatic colonization have been attributed to this clonal variance [4].Therefore understanding the extent of difference between these clonal variants is crucial for effectively targeting these clones.The characterization using in vivo models of the clonal heterogeneity arising from a single patient's colon tumor is still discrete.
Numerous in vitro and in vivo assays have been developed to study CRC progression.However, none of these techniques are effective in recapitulating the multi-step dissemination process.We have developed an orthotopic mouse CRC metastasis model that can quantitatively and qualitatively reproduce the metastatic phenotype of the human disease to the liver and/or lungs in an in vivo setting [1,5,6,7,8,9].Earlier work from our lab has characterized several human colon carcinoma cell lines using the orthotopic model [1,6,10].In this study, we compared the iso-clonal human colon carcinoma cell lines HCT116 and HCT116b isolated from the same patient primary colon carcinoma [11].Previously, we have demonstrated in vitro that HCT cells are growth factorindependent [12].In contrast, HCT116b cells are growth factordependent subcompartment of the malignant HCT116 cells [12].These isogenic cell lines demonstrate the clonal variance associated with malignant progression within tumors in vivo [6,12].
Comparison of HCT116 and HCT116b cells in vivo showed a significant difference in their ability to form metastatic deposits.Microarray analysis comparing the primary colon carcinoma arising from the two iso-clonal variants revealed striking differences in their gene signature.In vitro analysis of the two iso-clonal sublines revealed differences in cell survival and motility signaling under GFDS conditions.

Cell Culture
HCT116 and HCT116b sublines were isolated from a primary tissue culture of a single human colon carcinoma as described by Brattain et al. [11].The two cell lines were adapted to growth in serum-free medium [11,13,14] consisting of McCoy's 5A medium (Sigma) supplemented with amino acids, pyruvate, and antibiotics (designated SM) containing the growth factors transferrin (4 mg/ ml; Sigma), insulin (20 mg/ml; Sigma) and EGF (5 ng/ml; Collaborative Research) [15].

Green Fluorescence Protein (GFP) Transfection
Packaging cells, 293 GP (Clontech, Mountain View, CA), were co-transfected with a plasmid encoding VSVG envelope protein and a retroviral vector encoding GFP and the G418 resistance gene using FuGene (Invitrogen, Carlsbad, CA).The viruses were collected 48 h later and used to infect HCT116 and HCT116b cells.After 48 h, the infected HCT116 and HCT116b cells were selected by treatment with G418 for 5 d.This resulted in a stable transfection.

Animals
Male athymic BALB/c nude mice between 4-6 weeks of age were acquired from the National Cancer Institute (NCI).Housing for these animals was maintained in a HEPA-filtrated environment within sterilized cages.All animals were subjected to a daily 12 hr light/12-hr dark cycle.All animal procedures were conducted with approval of and in compliance with the Roswell Park Cancer Institute Institutional Animal Care and Use Committee (IACUC).

Orthotopic Implantation and Imaging
Exponentially growing GFP-labeled HCT116 and HCT116b cells (5610 6 cells) were inoculated subcutaneously onto the dorsal surfaces of separate BALB/c nude male mice.Once the xenografts were established, measuring approximately 500 mm3, they were excised and minced into 1 mm3 pieces.Two of these pieces were then orthotopically implanted into the colons of other male BALB/c nude mice as previously reported.Briefly, for operative procedures, animals were anesthetized with isoflurane inhalation.A 1 cm laparotomy was performed, and both the cecum and ascending colon were exteriorized.Using 76 magnification and microsurgical techniques, the serosa was disrupted in two different locations.The 1 mm3 pieces of xenograft were sub-serosally implanted using a 9-0 nylon suture at the two points of serosal disruption.The bowel was then returned to the peritoneal cavity and the abdomen was closed with interrupted 5-0 vicryl sutures.

Imaging
Starting at 1 wk post-implantation, animals were anesthetized with a 1:1 mixture of ketamine (10 mg/mL) and xylazine (1 mg/ mL) by intraperitoneal injection (0.01 mL/mg) and weekly GFP fluorescence imaging was performed for up to 7 wk.Specifically, GFP fluorescence imaging was performed using a light box illuminated by fiberoptic lighting at 470 nm (Illumatool BLS; Lightools Research, Encinitas, CA) using a Retiga EXi color CCD camera (QImaging, Burnaby BC, Canada).High-resolution images were captured directly using a MS-Windows based PC.Images were visually optimized for contrast and brightness using commercial software (Adobe Photoshop, CS2; Adobe, San Jose, CA).Excitation of GFP in the light box facilitated identification of primary and metastatic disease by direct near-real time visualization of fluorescence in live animals.

H and E Staining
At approximately 6-8 weeks post-implantation, the animals were euthanized.The colon, including any primary tumor, was harvested in addition to the liver, lungs, and heart.These organs were explanted, imaged, and immediately placed in 10% neutral buffered formalin fixative for 24 h.The tissues were then processed and embedded in paraffin.Slides were then cut for hematoxylin and eosin (H and E).

TUNEL and Ki67 Staining
Serial sections were cut to complement the H and E sections and were stained with an IgG1 rabbit polyclonal antibody for Ki-67 (Dako North America, Inc., Carpinteria, CA).Slides from paraffin embedded tissue blocks were stained according to the Apotag (Oncor, Gaithersburg, MD) terminal nucleotidyl transferase-mediated nick end labeling (TUNEL) assay kit.The Ki67 and TUNEL assay protocols have been described in details in our previous publications [1,6,7,8,9].

Apoptosis Assays
Apoptosis assays were performed to determine the response of HCT116 and HCT116b cells to growth factor deprivation stress (GFDS) [16].Cells were seeded in 24-well plates and allowed to grow to 80% confluent.The cells were then changed to starving medium (SM) for 2 or 4 days.Apoptosis assays were then performed using a DNA fragmentation ELISA kit as described in the manufacturer's protocol (Invitrogen).

Transwell Motility Assays
Transwell motility assays were performed utilizing 8mm pore, 6.5 mm polycarbonate transwell filters (Corning Costar Corp.).After trypsinizing the cells, single cell suspensions were seeded in supplemental McCoy's 5A medium in the absence of growth factors onto the upper surface of the filters and allowed to migrate towards McCoy's 5A medium with 10% FBS.After 18 hr incubation, MTT was added to the medium.The cells on the upper surface of the filter were removed with a cotton swab, and the cells that had migrated to the underside of the filter were visualized under the microscope and dissolved in DMSO.Absorbance was read at 570 nm.

Western Blot Analysis
Cells were lysed in a Tris-HCl based buffer containing 0.5% NP-40 and appropriate protease and phosphate inhibitors as described previously [1,15].Protein concentration was determined by bicinhoninic acid assay (Pierce).

Gene Expression Analysis
The microarray dataset comparing the HCT116 and HCT116b primary colon tumors has been deposited at the Gene Expression Omnibus (GEO).The accession number to the dataset is GSE44381.

Bioinformatics Analysis
Bioinformatics analysis was performed using the ingenuity pathway analysis (IPA) software tool obtained from the UNMC bioinformatics core facility.

High variance in metastatic rates between HCT116 and HCT116b cells
We compared the clonal variation in metastatic capability between HCT116 cells and its iso-clonal variant HCT116b by utilizing the colonic orthotopic implantation of subcutaneously grown xenografts.The HCT116 animals were euthanized between 6-8 weeks.The HCT 116b were sacrificed on day 63 (9 weeks).The HCT116 animals had to be put down earlier because of tumor burden and animal cachexia in accordance with IACUC regulations.As shown in Table 1, both cell lines gives rise to an invasive primary tumor in 100% of the animals implanted as demonstrated by invasive sheet observed during the immunohistochemical evaluation of the HCT116 and HCT116b primary tumors (Figure 1A).Strikingly, however, the cell lines vary dramatically with respect to their metastatic potential.HCT116 implanted mice generated histologically confirmed metastasis to liver and/or lungs in 68% animals (Figure 1B).HCT116 cells colonized to the liver in 21/40 (52%) and to the lungs in 19/40 (48%) of mice.In contrast, HCT116b cells showed a liver and/or lungs metastatic rate of only about 8%.HCT116b cells colonized to the liver in 2/49 (4%) and to the lungs in 2/49 (4%) of mice (Table 1).

Increased TUNEL staining observed in HCT116b primary tumors
We next characterized the HCT116 and HCT116b primary colonic tumors for their effect on cell survival and proliferation.As indicated in Figure 2, HCT116 showed an increased cell survival signaling as reflected by lower TUNEL staining rates compared to HCT116b tumors indicating that the repression of metastatic colonization in HCT116b is associated with repression of cell survival signaling in vivo.Interestingly, no change in Ki67 staining was observed between HCT116 and HCT116b primary tumors indicating no difference in proliferation rate in vivo between these clonal variants (data not shown).

Microarray analysis profiling gene signatures associated with HCT116 and HCT116b tumors
Next we sought to determine the differences in gene expression between HCT116 and HCT116b primary colon carcinoma tumor samples.Transcription profiles of the samples were generated using the Affymetrix HGU133plus2.0genechips.A heat map dendrogram generated and ranked using the 2-fold up-or downregulation cut-off is shown in the Figure S1.A summary of genes differentially regulated between the two iso-clonal tumors have also been shown in Table S1.Gene ontology analysis was performed to delineate the differences in gene expression between different cellular compartments.As shown in Figure S2, 26% of all genes were differentially regulated in the cytoplasm.There was a 23% difference in gene expression profiles in the plasma membrane followed by 17% in the nucleus and 10% in the extracellular space.We next performed the Ingenuity Pathway Analysis (IPA) to categorize the difference in gene signatures between HCT116 and HCT116b tumors according to their molecular participation in cell survival and metastasis.Strikingly, several genes involved in providing cancer cells with the survival advantage were observed to be downregulated in HCT116b cells.Some of the genes that were differentially regulated in HCT116 and HCT116b cells have been shown in Figure 3.Among these are growth factor ligands IGF1, EGF, TGFb and FGF that are significantly repressed in HCT116b.Similarly, genes involved in conferring metastatic ability, such as MMPs, fibronectin 1 (FN1) and vimentin (VIM) are transcriptionally repressed in HCT116b.In contrast, cyclin dependent kinase inhibitors, such as CDKN2A and CDKN2C were upregulated in HCT116b cells.
In vitro analysis of cell death markers in HCT116 and HCT116b cells under GFDS Previously, we examined the effects of nutrient and growth factor deprivation in vitro in HCT116 and HCT116b human colon carcinoma cells [12].Under conditions of nutrient deprivationmediated growth arrest, HCT116b cells required nutrients and growth factor supplementation for the reinitiation of DNA synthesis.However, HCT116 only required nutrients and showed a transcriptional upregulation of TGFa during stress [12].Therefore the combined intrinsic ability of HCT116 cells to survive under stress and concomitantly re-enter cell cycle in response to nutrient exposure without exogenous growth factors would seem to provide the malignant cells with a substantial growth advantage in situations of limited blood supply and during the process of metastasis.Based on these prior observations and our current gene signatures obtained by microarray analysis, we hypothesized that the downstream survival pathways associated with malignant progression in CRC would be differentially regulated in HCT116 and HCT116b cell under stress conditions.
We first performed a DNA fragmentation assay to examine the level of cell death in these cells under nutrient and growth factor deprivation stress (GFDS) [12,16].Confluent cells in culture were used as control (reference point) and compared to quiescent (Q) cells that have been deprived of exogenous growth factors for 2 and 4 days respectively (termed as Q2 and Q4 respectively).HCT116 cells did not show any increase in apoptosis under GFDS for 4 days.However, HCT116b cells showed approximately 2-fold increase in cell death by Q4 (Figure 4A).We next performed western blot analysis of HCT116 and HCT116b cells under GFDS conditions to investigate the protein expression of different apoptotic markers.Interestingly, we observed a stable expression of total-PARP and -caspase-3 in HCT116 cells.In contrast, HCT116b cells showed a time-dependent increase in cleaved-PARP and -caspase-3 indicative of increased apoptosis.Similarly, we observed phosphorylation of anti-apoptotic marker pBad in HCT116 cells.However, pBad expression progressively reduced under GFDS in HCT116b cells (Figure 4B).

In vitro analysis of cell survival signaling pathways in HCT116 and HCT116b cells under GFDS
We next determined the level of phosphorylation under GFDS of major signal transduction pathways such as the mitogen activated protein kinase (MAPK), phosphatidylinositol 3-kinase (PI 3 K) and phospholipase Cc (PLCc) pathways which control proliferation and cell survival in colon cancer [8,16,17,18].As shown in Figure 5, we observed a reduction in phosphorylation of AKT, mTOR, ERK and PLCc proteins in HCT116b cells under GFDS.However, HCT116 cells continued to demonstrate high phosphorylation activity of these oncogenes indicative of its growth-factor independent cell survival machinery to counter stress.

In vitro analysis of cell motility related signaling pathways in HCT116 and HCT116b cells under GFDS
Invasion and migration are the two initial steps that facilitate the multi-step cascade of dissemination of cancer cells to distant organ sites [19].We determined the migratory properties of HCT116 and HCT116b cells under conditions of GFDS for 5 days using the transwell motility assay as described in the materials and methods section.Under GFDS, the HCT116b showed approximately 90% reduction in motility compared to HCT116 cells (Figure 6A).To identify the contributions of oncogenic signaling molecules in regulating motility of these cells, we performed western blot analysis and analyzed for several biochemical markers associated with invasion and migration of cancer cells.RhoA, ROCK1 and Rac1 proteins are GTPases that regulate F-actin and thereby regulate cancer cell migration and have been implicated in metastasis of several types of cancer including CRC [19].As described in Figure 6B, RhoA, Rac1 and ROCK1 protein expression decreased in HCT116b cells under GFDS.In contrast, under stress conditions there was an increase in the protein expression of these GTPases in HCT116 cells.We further observed an increase in the phosphorylation of FAK and paxillin in HCT116 cells under GFDS that are also implicated in CRC progression and metastasis [20].

Discussion
In this study, we explored the phenotypic heterogeneity of clonal variants isolated from the same patient primary colon carcinoma using an orthotopic mouse metastasis model for studying CRC metastasis.As evident from the recent review by Marusyk et al [4], there are both genetic and non-genetic causes of phenotypic heterogeneity of cancer cells.The intra-tumoral heterogeneity has substantial implication in the prognosis, subtype categorization, risk of metastasis and the development of therapeutic resistance [4].Inter-tumoral heterogeneity has been  well studied and has provided several model systems to study the initiation and progression of cancer.However, the drawback of these model systems is that the comparison of non-isogenic cell lines has the potential to generate erroneous conclusions since the biological and growth regulatory properties may reflect individual variations of such cells [12].We utilized the iso-clonal human colon carcinoma cell lines HCT116 and HCT116b to study the intra-tumoral heterogeneity in colon cancer cells.These cell lines were isolated and developed into primary tissue cultures directly from a single primary tumor, thus providing a unique opportunity to assess iso-clonal variations in genetic heterogeneity in an in vivo metastatic setting.Initial studies comparing these two clones showed differences in their cellular morphology [11].The HCT116b cells were differentiated and weakly tumorigenic in comparison with the HCT116 cell lines [11].HCT116 cells have been extensively characterized and demonstrated exogenous growth factors independence for DNA synthesis and growth.Awwad et al [12] made the striking observation that endogenously produced growth factor activation in colon cancer cells lead to their exogenous growth factor independence.These cells can survive under nutrient and growth factor deprivation in comparison to their iso-clonal variants that are growth-factor dependent.This report demonstrated that the development of growth factor independent phenotype represents clonal progression within the tumor in vivo.
This study investigated the metastatic properties of the two isoclones utilizing an in vivo mouse metastasis model developed in our laboratory.This model utilizes orthotopic xenograft transplants from human CRC cell lines that recapitulate the nature of metastatic pattern of the disease in CRC patients in that liver and lung metastases are observed after the growth of primary tumors in the colon [1,6,7,8,9,10].The HCT116 cells were found to be highly metastatic with a metastatic rate of 68% in liver and/or lungs.The iso-clonal weakly tumorigenic and differentiated cell line HCT116b was found to be poorly metastatic with a metastatic colonization rate of about 4%.We observed high variability in the primary tumor burden between the two experimental groups (HCT116 and HCT116b) as well as within the same group receiving orthotopically implanted xenograft tumors (HCT116 or HCT116b).Due to this inter-and intra-tumor variability in tumor burden, we performed Ki-67 staining to measure the cell proliferation rate between several HCT116 and HCT116b primary tumors in order to address the possibility that higher level of metastasis might be a consequence of increased primary tumor burden.Interestingly, HCT116 and HCT116b implanted animals had no difference in their rate of cell proliferation in the primary tumors, irrespective of the tumor burden as indicated by Ki-67 staining.However, HCT116 primary tumors have a significantly lower level of TUNEL staining indicative of low apoptotic index, and enhanced cell survival capabilities compared to HCT116b primary tumors.
In this work, we made the novel observation that both iso-clonal variants HCT116 and HCT116b cells showed 100% primary tumor invasion but different rates of metastases.Further analysis of the HCT116 and HCT116b primary colonic tumors showed similar rates of cell proliferation as demonstrated by Ki-67 staining.In contrast, these tumors showed clonal variability in their cell survival capabilities.HCT116b primary tumors showed increased TUNEL staining, which is indicative of cell death compared to HCT116 primary tumors.These findings indicate that the HCT116 and HCT116b cells have developed variability in their ability to form progressively growing colonies at the distal sites rather than modulating their dissemination from the primary tumors sites.This result also suggested that the low rate of distant metastasis in HCT116b cells might be related to stress-induced cell death in the foreign microenvironment for growth of the colon cancer cells in the liver and/or lungs.As evident from the in vitro data presented in Figures 4-6, induction of stress by growth factor deprivation in HCT116b cells led to increased cell death through the downregulation of several key oncogenic molecules that are known to participate in the regulation of aberrant cell survival and metastasis.Ki-67 and TUNEL staining has been used extensively in numerous clinical and pre-clinical studies.Ki-67 has been clinically relevant having predictive and prognostic value.In contrast, TUNEL staining, which is indicative of apoptotic index has been extensively linked to patient survival.The results obtained in this study are significant because they raise the clinical importance of cell survival signaling that could be playing a major role in the dissemination of colon cancer cells.Therefore, successful targeting of these pathways would be an effective strategy to target the treatment of metastasis.
We assessed the gene signatures associated with these clones using microarray and bioinformatics analysis.As represented in the heat map dendrogram (Figure S1), as expected numerous genes were differentially regulated between HCT116 and HCT116b cells.Several of these genes are markers associated with increased stemness, angiogenesis, invasion-migration, growth, survival and metastasis of the HCT116 clones.We focused our attention specifically to the genes affecting the cell survival and metastasis functions of CRC.The linkage of aberrant cell survival signaling for various steps in metastasis has been established with particular emphasis on the importance of aberrant cell survival to successful colonization at the metastatic site [21].Therefore, the understanding of mechanisms that govern cell survival fate of highly metastatic cells would lead to the understanding of a new paradigm for control of metastatic potential and would provide the basis for developing novel strategies for the treatment of metastases.Several growth factor signaling pathways were differentially regulated between the HCT116 and HCT116b cells analyzed under GFDS conditions for 3 to 5 days.An important criterion selected for this study was the selection of up to Day 5 under quiescence state (Q5) as the time at which the experiments were performed.Previous work by Wang et al [22] has shown that at this stage, these cells show a minimal baseline incorporation of 3 H thymidine yet retain greater than 95% viability.Basal 3 H thymidine incorporation is minimal because nutrients from the growth factor free culture medium have been exhausted. 3H thymidine incorporation is restored by addition of fresh nutrients to these cells without any exogenous growth factors [23].Restriction of nutrients past 5 days results in release of the cells from the culture plate and death.Therefore, we posited that operationally, 5 days post-confluency should be a time when the cell's endogenous cell survival and cell cycle re-entry capabilities are maximized.Since these functions are closely aligned with malignancy it is reasonable to project that the endogenous signaling mechanisms exhibited by these cells at this point are relevant to the maintenance of the metastatic phenotype of these cells.The PI3K/AKT, MAPK and PLC pathways are implicated in cell survival and metastasis of cancer cells.In this study, we demonstrate for the first time, the variation of activation of these signaling nodes under GFDS comparing the two iso-clonal sublines.We further demonstrated the increased expression of proteins involved in invasion-migration of CRC cells as a demonstration of the cancer cells extraordinary ability to evade stress mechanisms under anchorage-independent conditions.
The subject of intra-tumoral phenotypic heterogeneity in metastatic dissemination remains mostly unexplored.Marusyk et al [4] documented several distinct sources of phenotypic heterogeneity within tumor cell population.It has been reported that such variation in tumor population are an integration of multiple genetic and non-genetic factors.Previously, the association of intra-tumor heterogeneity with poor prognosis has been demonstrated in esophageal and breast cancer [24,25].In this study, we  demonstrate the relevance of intra-tumoral heterogeneity in colorectal cancer metastasis using the orthotopic metastasis mouse model.The highly metastatic HCT116 clonal variant showed remarkable heterogeneity in gene signatures associated with cell survival and metastatic capabilities in comparison to the poorly metastatic HCT116b clones.The HCT116 cells displayed resistance to nutrient stress by upregulating the protein expression of several key oncogenes associated with aberrant cell survival and metastasis of CRC.In contrast, the poorly metastatic isoclonal variant HCT116b cells downregulated these proteins under similar GFDS conditions demonstrating strikingly opposite response to induce cell death under survival-associated stress.

Figure 1 .
Figure 1.Intra-tumoral heterogeneity in the ability to develop metastatic colonization.(A) H&E section demonstrating normal colonic area and primary colon carcinoma (magnification 406).(B) 4 to 6 week old male BALB/c nude mice were orthotopically implanted with HCT116 or HCT116b tumor xenografts.GFP imaging revealed primary tumor growth in both HCT116 and HCT116b animals.However, distant metastasis was only observed in HCT116 implanted animals by week 6.doi:10.1371/journal.pone.0060299.g001

Figure 2 .
Figure 2. Increased cell survival associated with HCT116 cells.Comparison of primary tumor sections of HCT116 and HCT116b mice by TUNEL staining as mentioned in Materials and Methods to determine their apoptotic rates.HCT116b tumor tissues had high TUNEL staining indicative of cell death.However, HCT116 tumor tissues showed negligible response to TUNEL staining indicating high rate of cell survival.doi:10.1371/journal.pone.0060299.g002

Figure 3 .
Figure 3.Comparison of genes signatures associated with metastasis signaling between HCT116 and HCT116b microarray datasets by Ingenuity Pathway Analysis (IPA) tool.The molecules highlighted in green are upregulated in the highly metastatic HCT116 cells and downregulated in poorly metastatic iso-clonal HCT116b cells.The molecules marked in red are upregulated in HCT116b cells compared to HCT116 cells.doi:10.1371/journal.pone.0060299.g003

Figure 4 .
Figure 4. Cell death response under growth factor deprivation stress (GFDS).(A) HCT116 cells are resistant to GFDS-induced cell death as determined by DNA fragmentation assay.However, HCT116b iso-clonal cells induce cell death under GFDS.(B) HCT116b cells induce cell death by time-dependent increase in PARP and caspase 3 cleavages and dephosphorylation of anti-apoptotic pBad protein.doi:10.1371/journal.pone.0060299.g004

Figure 5 .Figure 6 .
Figure 5. Cell survival signaling response under growth factor deprivation stress (GFDS).HCT116 cells upregulated AKT, MAPK and PLC signaling under GFDS conditions to counter stress.In contrast, HCT116b cells downregulated these signaling molecules by dephosphorylating and thereby inactivating AKT, ERK and PLC gamma to induced cell death.doi:10.1371/journal.pone.0060299.g005

Table 1 .
Results of orthotopic implantation of HCT116 and HCT116b tumor xenografts demonstrating the primary invasion and metastases based on histologic evaluation. doi:10.1371/journal.pone.0060299.t001