Stochastic dynamics of metastasis formation

Tumor metastasis accounts for the majority of deaths in cancer patients. The metastatic behavior of cancer cells is promoted by mutations in many genes, including activation of oncogenes such as RAS and MYC. Here, we develop a mathematical framework to analyse the dynamics of mutations enabling cells to metastasize. We consider situations in which one mutation is necessary to confer metastatic ability to the cell. We study different population sizes of the main tumor and different somatic fitness values of metastatic cells. We compare mutations that are positively selected in the main tumor with those that are neutral or negatively selected, but faster at forming metastases. We study whether metastatic potential is the property of all (or the majority of) cells in the main tumor or only the property of a small subset. Our theory shows how to calculate the expected number of metastases that are formed by a tumor.     

Somatic selection for and against cancer

In multicellular organisms, cells cooperate within a well-defined developmental program. Cancer is a breakdown of such cooperation: cells mutate to phenotypes of uncoordinated proliferation. We study basic principles of the architecture of solid tissues that influence the rate of cancer initiation. In particular, we explore how somatic selection acts to prevent or to promote cancer. Cells with mutations in oncogenes or tumor suppressor genes often have increased proliferation rates. Somatic selection increases their abundance and thus enhances the risk of cancer. Many potentially harmful mutations, however, increase the probability of triggering apoptosis and, hence, initially lead to cells with reduced net proliferation rates. Such cells are eliminated by somatic selection, which therefore also works to reduce the risk of cancer. We show that a tissue organization into small compartments avoids the rapid spread of mutations in oncogenes and tumor suppressor genes, but promotes genetic instability. In small compartments, genetic instability, which confers a selective disadvantage for the cell, can spread by random drift. If both deleterious and advantageous mutations participate in tumor initiation, then we find an intermediate optimum for the compartment size.     

Growth, metastasis, and invasiveness of Drosophila tumors caused by mutations in specific tumor suppressor genes

 More than 50 genes have been identified in Drosophila by loss-of-function mutations that lead to overgrowth of specific tissues. Loss-of-function mutations in the lethal giant larvae, discs large, or brain tumor genes cause neoplastic overgrowth of larval brains and imaginal discs. In the present study, the growth and metastatic potential of tumors resulting from mutations in these genes were quantified. Overgrown brains and imaginal discs were transplanted into adults and β-galactosidase accumulation was used as a marker to identify donor cells. Mutations in these three genes generated tumors with similar metastatic patterns. For brain tumors, the metastatic index (a measure we defined as the fraction of hosts that acquired secondary tumors normalized for the amount of primary tumor growth) of each of the three mutants was similar. Analysis of cell proliferation in mutant brains suggests that the tumors arise from a population of several hundred cells which represent only 1–2% of the cells in third instar larval brains. For imaginal disc tumors from lethal giant larvae and brain tumor mutants, it is shown for the first time that they can be metastatic and invasive. Primary imaginal disc tumors from lethal giant larvae and brain tumor mutants formed secondary tumors in 43 and 53% of the hosts, respectively, although the secondary tumors were, in general, smaller than the secondary tumors derived from primary brain tumors. Received: 18 August 1997 / Accepted: 16 October 1997     

Population genetics of tumor suppressor genes

Cancer emerges when a single cell receives multiple mutations. For example, the inactivation of both alleles of a tumor suppressor gene (TSG) can imply a net reproductive advantage of the cell and might lead to clonal expansion. In this paper, we calculate the probability as a function of time that a population of cells has generated at least one cell with two inactivated alleles of a TSG. Different kinetic laws hold for small and large populations. The inactivation of the first allele can either be neutral or lead to a selective advantage or disadvantage. The inactivation of the first and of the second allele can occur at equal or different rates. Our calculations provide insights into basic aspects of population genetics determining cancer initiation and progression.     

Encapsulated human hepatocellular carcinoma cells by alginate gel beads as an in vitro metastasis model

Hepatocellular carcinoma (HCC) is the most common primary liver cancer and often forms metastases, which are the most important prognostic factors. For further elucidation of the mechanism underlying the progression and metastasis of HCC, a culture system mimicking the in vivo tumor microenvironment is needed. In this study, we investigated the metastatic ability of HCC cells cultured within alginate gel (ALG) beads. In the culture system, HCC cells formed spheroids by proliferation and maintained in nuclear abnormalities. The gene and protein expression of metastasis-related molecules was increased in ALG beads, compared with the traditional adhesion culture. Furthermore, several gene expression levels in ALG bead culture system were even closer to liver cancer tissues. More importantly, in vitro invasion assay showed that the invasion cells derived from ALG beads was 7.8-fold higher than adhesion cells. Our results indicated that the in vitro three-dimensional (3D) model based on ALG beads increased metastatic ability compared with adhesion culture, even partly mimicked the in vivo tumor tissues. Moreover, due to the controllable preparation conditions, steady characteristics and production at large-scale, the 3D ALG bead model would become an important tool used in the high-throughput screening of anti-metastasis drugs and the metastatic mechanism research.     

Role of platelets and breast cancer stem cells in metastasis

The high mortality rate of breast cancer is mainly caused by the metastatic ability of cancer cells, resistance to chemotherapy and radiotherapy, and tumor regression capacity. In recent years, it has been shown that the presence of breast cancer stem cells is closely associated with the migration and metastatic ability of cancer cells, as well as with their resistance to chemotherapy and radiotherapy. The tumor microenvironment is one of the main molecular factors involved in cancer and metastatic processes development, in this sense it is interesting to study the role of platelets, one of the main communicator cells in the human body which are activated by the signals they receive from the microenvironment and can generate more than one response. Platelets can ingest and release RNA, proteins, cytokines and growth factors. After the platelets interact with the tumor microenvironment, they are called "tumor-educated platelets." Tumor-educated platelets transport material from the tumor microenvironment to sites adjacent to the tumor, thus helping to create microenvironments conducive for the development of primary and metastatic tumors. It has been observed that the clone capable of carrying out the metastatic process is a cancer cell with stem cell characteristics. Cancer stem cells go through a series of processes, including epithelial-mesenchymal transition, intravasation into blood vessels, movement through blood vessels, extravasation at the site of the establishment of a metastatic focus, and site colonization. Tumor-educated platelets support all these processes.     

Genetic instability and clonal expansion

Inactivation of tumor suppressor genes can lead to clonal expansion. We study the evolutionary dynamics of this process and calculate the probability that inactivation of a tumor suppressor gene is preceded by mutations in genes that confer genetic instability. Unstable cells might have a slower rate of clonal expansion than stable cells because of an increased probability of generating lethal mutations or inducing apoptosis. We show that the different growth rates of genetically stable and unstable cells during clonal expansion represent, in general, only a small disadvantage for genetic instability. The intuitive reason for this conclusion is that robust clonal expansion, where cellular birth rates are significantly greater than death rates, occurs on a much faster time scale than waiting for those mutations that allow clonal expansion. Moreover, in special cases where clonal expansion is very slow, genetically unstable cells have a higher probability to accumulate additional mutations during clonal expansion that confer a selective advantage. Clonal expansion represents a major disadvantage for genetic instability only when inactivation of the tumor suppressor gene leads to a very small increase of the cellular reproductive rate that is cancelled by the increased mortality of unstable cells.     

More is less: Inactivation and deletion events and the search for tumor suppressor genes

Tumor suppressor genes are frequently inactivated in cancer by large‐scale deletion events or epigenetic silencing, and experimental demonstration of such inactivation has historically been considered as support for assigning tumor suppressive function to a given gene. However, the discovery of a number of chromosomal domains wherein large deletions naturally occur at frequencies up to 100 times the average for the genome as a whole leads us to reevaluate the significance of sporadic deletions found within genes associated with these hotspots. Similarly, our recent demonstration that epigenetic chromatin silencing frequently spreads in cancer cells from gene‐poor into gene‐rich regions with apparent indifference to the gene content of the affected domain raises questions about the pertinence of inactivation as a criterion for ascribing tumor suppressor function to a given gene. We suggest that a number of putative suppressor genes for which inactivation and/or deletion events have been documented may simply be victims of collateral damage when these events occur, and the implication that these genes are being selected against during cancer progression should in some cases be reassessed. J. Cell. Biochem. 110: 281–287, 2010. © 2010 Wiley‐Liss, Inc.     

Translational approaches using metastasis suppressor genes

Cancer metastasis is a significant contributor to breast cancer patient morbidity and mortality. In order to develop new anti-metastatic therapies, we need to understand the biological and biochemical mechanisms of metastasis. Toward these efforts, we and others have studied metastasis suppressor genes, which halt metastasis in vivo without affecting primary tumor growth. The first metastasis suppressor gene identified was nm23, also known as NDP kinase. Nm23 represents the most widely validated metastasis suppressor gene, based on transfection and knock-out mouse strategies. The biochemical mechanism of metastasis suppression via Nm23 is unknown and likely complex. Two potential mechanisms include binding proteins and a histidine kinase activity. Elevation of Nm23 expression in micrometastatic tumor cells may constitute a translational strategy for the limitation of metastatic colonization in high risk cancer patients. To date, medroxyprogesterone acetate (MPA) has been identified as a candidate compound for clinical testing.     

p53 and PCNA Expression in Malignant Melanomas of the Head and Neck

Mutation in the p53 tumor suppressor gene is the most common genetic alteration in human cancer. As in mutant p53 the protein is stabilised and the half-life is extended, it becomes detectable by immunohistological staining. p53 immunoreactivity thus seems to be a potential biomarker for the assessment of the oncogenic potential of malignant melanomas. In 103 tissue sections of primary and metastatic malignant melanomas of the head and neck detectable levels of p53 were only found in 3 of the primary tumors and in none of the metastases. At the same time the proliferation status of the malignant melanoma lesions was determined using the cell cycle specific antibody PCNA. 55 primary and metastatic tumors were stained with a PCNA-MAb to determine the proliferation activity of the tumors. The results of our immunohistochemical investigation suggest that immunoreactivity of p53 cannot be used to determine the malignant potential of melanomas in the head and neck. PCNA staining showed that the majority of the tumors and metastases were proliferating rapidly.     

Differential expression of metabolic genes in tumor and stromal components of primary and metastatic loci in pancreatic adenocarcinoma

Background

Pancreatic cancer is the fourth leading cause of cancer related deaths in the United States with a five-year survival rate of 6%. It is characterized by extremely aggressive tumor growth rate and high incidence of metastasis. One of the most common and profound biochemical phenotypes of animal and human cancer cells is their ability to metabolize glucose at high rates, even under aerobic conditions. However, the contribution of metabolic interrelationships between tumor cells and cells of the surrounding microenvironment to the progression of cancer is not well understood. We evaluated differential expression of metabolic genes and, hence, metabolic pathways in primary tumor and metastases of patients with pancreatic adenocarcinoma.

Methods and Findings

We analyzed the metabolic gene (those involved in glycolysis, tri-carboxylic acid pathway, pentose-phosphate pathway and fatty acid metabolism) expression profiles of primary and metastatic lesions from pancreatic cancer patients by gene expression arrays. We observed two principal results: genes that were upregulated in primary and most of the metastatic lesions; and genes that were upregulated only in specific metastatic lesions in a site-specific manner. Immunohistochemical (IHC) analyses of several metabolic gene products confirmed the gene expression patterns at the protein level. The IHC analyses also revealed differential tumor and stromal expression patterns of metabolic enzymes that were correlated with the metastasis sites.

Conclusions

Here, we present the first comprehensive studies that establish differential metabolic status of tumor and stromal components and elevation of aerobic glycolysis gene expression in pancreatic cancer.     

Transformation effector and suppressor genes.

Much has been learned about the molecular basis of cancer from the study of the dominantly acting viral and cellular oncogenes and their normal progenitors, the proto-oncogenes. More recent studies have resulted in the isolation and characterization of several genes prototypic of a second class of cancer genes. Whereas oncogenes act to promote the growth of cells, members of this latter class of genes act to inhibit cellular growth and are believed to contribute to the tumorigenic phenotype only when their activities are absent. This new class of cancer genes is referred to by a number of different names including; anti-oncogenes, recessive oncogenes, growth suppressor genes, tumor suppressor genes and emerogenes. Although only a few of these cancer genes have been identified, to date, it is likely that many additional genes of this class await identification. A third class of genes, necessary for the development of the cancer phenotype, is comprised of the transformation effector genes. These are normal cellular genes that encode proteins that cooperate with or activate oncogene functions and thereby induce the development of the neoplastic phenotype. The inactivation of transformation effector functions would therefore inhibit the ability of certain dominantly acting oncogenes to transform cells. The approaches outlined here describe functional assays for the isolation and molecular characterization of transformation effector and suppressor genes.     

Murine tumor suppressor models

Tumor suppressor genes have been shown to be necessary for proper maintenance of cell growth control. Inactivation of these genes in the germline of humans is linked to inherited cancer predisposition. Moreover, sporadically arising human tumors often have somatic mutations in tumor suppressor genes. During the past few years, advances in molecular and cellular biology have led to the creation of animal models that have germline mutations of various tumor suppressor genes. Such mice potentially represent important animal models for familial cancer predisposition syndromes, and the study of the tumorigenesis process has been greatly assisted by their development. Such models have also demonstrated the importance of tumor suppressor function in embryonic development. In this review, we describe mice with inactivated germline tumor suppressor genes that are genetically analogous to 10 different inherited cancer syndromes in humans. We describe the variable usefulness of the mutant mice as models for human disease.     

Chromothripsis is a common mechanism driving genomic rearrangements in primary and metastatic colorectal cancer

Background

Structural rearrangements form a major class of somatic variation in cancer genomes. Local chromosome shattering, termed chromothripsis, is a mechanism proposed to be the cause of clustered chromosomal rearrangements and was recently described to occur in a small percentage of tumors. The significance of these clusters for tumor development or metastatic spread is largely unclear.

Results

We used genome-wide long mate-pair sequencing and SNP array profiling to reveal that chromothripsis is a widespread phenomenon in primary colorectal cancer and metastases. We find large and small chromothripsis events in nearly every colorectal tumor sample and show that several breakpoints of chromothripsis clusters and isolated rearrangements affect cancer genes, including NOTCH2, EXO1 and MLL3. We complemented the structural variation studies by sequencing the coding regions of a cancer exome in all colorectal tumor samples and found somatic mutations in 24 genes, including APC, KRAS, SMAD4 and PIK3CA. A pairwise comparison of somatic variations in primary and metastatic samples indicated that many chromothripsis clusters, isolated rearrangements and point mutations are exclusively present in either the primary tumor or the metastasis and may affect cancer genes in a lesion-specific manner.

Conclusions

We conclude that chromothripsis is a prevalent mechanism driving structural rearrangements in colorectal cancer and show that a complex interplay between point mutations, simple copy number changes and chromothripsis events drive colorectal tumor development and metastasis.     

Data on the Recurrence of Breast Tumors Fit a Model in which Dormant Cells are Subject to Slow Attrition but Can Randomly Awaken to Become Malignant

We successfully modeled the recurrence of tumors in breast cancer patients, assuming that:(i) A breast cancer patient is likely to have some circulating metastatic cells, even after initial surgery. (ii) These metastatic cells are dormant. (iii) The dormant cells are subject to attrition by the body’s immune system, or by random apoptosis or senescence.(iv) Recurrence suppressor mechanisms exist. (v) When such genes are disabled by random mutations, the dormant metastatic cell is activated, and will develop to a cancer recurrence. The model was also fitted to data on the survival of pancreatic cancer patients. The time course of cancer recurrence in a group of poor prognosis breast cancer patients could not be linked to the over- (or under-) expression of any gene in the primary tumors from which the recurrent tumors derived. Thus, the recurrence of the tumor in breast cancer patients appears to be a random event. Inasmuch as the kinetics of cancer recurrence in published data sets closely follows the model found for the appearance of sporadic retinoblastoma, tumor recurrence could be triggered by mutations in awakening-suppressor mechanisms. The retinoblastoma tumor suppressor gene was identified by tracing its occurrence in familial retinoblastoma pedigrees. Will it be possible to track the postulated cancer recurrence, awakening suppressor gene(s) in early recurrence breast cancer patients?     

Sine oculis homeobox homolog 1 promotes α5β1-mediated invasive migration and metastasis of cervical cancer cells

Sine oculis homeobox homolog 1 (SIX1) has been supposed to be correlated with the metastasis and poor prognosis of several malignancies. However, the effect of SIX1 on the metastatic phenotype of tumor cells and the underlying mechanisms were still unclear to date. Here we report that SIX1 can promote α5β1-mediated metastatic capability of cervical cancer cells. SIX1 promoted the expression of α5β1 integrin to enhance the adhesion capacity of tumor cells in vitro and tumor cell arrest in circulation in vivo. Moreover, higher expression of SIX1 in tumor cells resulted in the increased production of active MMP-2 and MMP-9, up-regulation of anti-apoptotic genes (BCL-XL and BCL2) and down-regulation of pro-apoptotic genes (BIM and BAX), thus promoting the invasive migration and anoikis-resistance of tumor cells. Importantly, blocking α5β1 abrogated the regulatory effect of SIX1 on the expression of these genes, and also abolished the promotional effect of SIX1 on invasive capability of tumor cells. Furthermore, knock-down of α5 could abolish the promoting effect of SIX1 on the development of metastatic lesions in both experimental and spontaneous metastasis model. Therefore, by up-regulating α5β1 expression, SIX1 not only promoted the adhesion capacity, but also augmented ECM-α5β1-mediated regulation of gene expression to enhance the metastatic potential of cervical cancer cells. These results suggest that SIX1/α5β1 might be considered as valuable marker for metastatic potential of cervical cancer cells, or a therapeutic target in cervical cancer treatment.     

Regulation of breast cancer-induced bone lesions by β-catenin protein signaling

Breast cancer patients have an extremely high rate of bone metastases. Morphological analyses of the bones in most of the patients have revealed the mixed bone lesions, comprising both osteolytic and osteoblastic elements. β-Catenin plays a key role in both embryonic skeletogenesis and postnatal bone regeneration. Although this pathway is also involved in many bone malignancy, such as osteosarcoma and prostate cancer-induced bone metastases, its regulation of breast cancer bone metastases remains unknown. Here, we provide evidence that the β-catenin signaling pathway has a significant impact on the bone lesion phenotype. In this study, we established a novel mouse model of mixed bone lesions using intratibial injection of TM40D-MB cells, a breast cancer cell line that is highly metastatic to bone. We found that both upstream and downstream molecules of the β-catenin pathway are up-regulated in TM40D-MB cells compared with non-bone metastatic TM40D cells. TM40D-MB cells also have a higher T cell factor (TCF) reporter activity than TM40D cells. Inactivation of β-catenin in TM40D-MB cells through expression of a dominant negative TCF4 not only increases osteoclast differentiation in a tumor-bone co-culture system and enhances osteolytic bone destruction in mice, but also inhibits osteoblast differentiation. Surprisingly, although tumor cells overexpressing β-catenin did induce a slight increase of osteoblast differentiation in vitro, these cells display a minimal effect on osteoblastic bone formation in mice. These data collectively demonstrate that β-catenin acts as an important determinant in mixed bone lesions, especially in controlling osteoblastic effect within tumor-harboring bone environment.