A quantitative theoretical model for the development of malignancy in ductal carcinoma in situ

Mathematical models and clinical observations have demonstrated that microenvironmental hypoxia and acidosis are important selection factors during the later stages of the somatic evolution of breast cancer. The consequent promotion of constitutive upregulation of glycolysis and resistance to acid-induced cellular toxicity is hypothesized to be critical for the ability of cancer cells to invade host tissue. In this work we developed a 3D fixed lattice cellular automata model to study the role of these two phenotypes in determining morphology and the potential for invasion of ductal carcinoma in situ (DCIS), which in this work is defined as the erosion of a healthy epithelial cell layer and direct contact with the basement membrane. The model was conceived as a 40-cell wide epithelial duct surrounded by blood vessels and composed of a basement membrane and one internal layer of epithelial cells. Our results show that an increment in the order of 8-fold in glucose metabolism and an increase in acid resistance corresponding to pH thresholds of approximately 6.8 and 6.45 for quiescence and death, respectively, are required for the tumor to breach through the layer of healthy epithelial cells and reach the basement membrane as a first step for invasion. Our model also suggests correlations between classic morphologies and different values of hyperglycolytic and acid-resistant phenotypes, indicating that immunohistochemistry studies targeting these genes may improve the predictive power of morphological analyses of biopsies.     

Oncogenic alterations of metabolism

Over seven decades ago, classical biochemical studies showed that tumors have altered metabolic profiles and display high rates of glucose uptake and glycolysis. Although these metabolic changes are not the fundamental defects that cause cancer, they might confer a common advantage on many different types of cancers, which allows the cells to survive and invade. Recent molecular studies have revealed that several of the multiple genetic alterations that cause tumor development directly affect glycolysis, the cellular response to hypoxia and the ability of tumor cells to recruit new blood vessels.     

Effects of Different Tissue Microenvironments on Gene Expression in Breast Cancer Cells

In metastasis, circulating tumor cells penetrate the walls of blood vessels and enter the metastatic target tissue, thereby becoming exposed to novel and relatively unsupportive microenvironments. In the new microenvironments, the tumor cells often remain in a dormant state indefinitely and must adapt before they are able to successfully colonize the tissue. Very little is known about this adaptive process. We studied temporal changes in gene expression when breast cancer cells adapt to survive and grow on brain, bone marrow, and lung tissue maintained in an in vivo culture system, as models of the metastatic colonization of these tissues. We observed the transient activation of genes typically associated with homeostasis and stress during the initial stages of adaptation, followed by the activation of genes that mediate more advanced functions, such as elaboration of cell morphology and cell division, as the cells adapted to thrive in the host tissue microenvironment. We also observed the temporary induction of genes characteristic of the host tissue, which was particularly evident when tumor cells were grown on brain tissue. These early transient gene expression events suggest potential points of therapeutic intervention that are not evident in data from well-established tumors.     

An update on therapeutic opportunities offered by cancer glycolytic metabolism

Almost all invasive cancers, regardless of tissue origin, are characterized by specific modifications of their cellular energy metabolism. In fact, a strong predominance of aerobic glycolysis over oxidative phosphorylation (Warburg effect) is usually associated with aggressive tumour phenotypes. This metabolic shift offers a survival advantage to cancer cells, since they may continue to produce energy and anabolites even when they are exposed to either transient or permanent hypoxic conditions. Moreover, it ensures a high production rate of glycolysis intermediates, useful as building blocks for fast cell proliferation of cancer cells. This peculiar metabolic profile may constitute an ideal target for therapeutic interventions that selectively hit cancer cells with minimal residual systemic toxicity. In this review we provide an update about some of the most recent advances in the discovery of new bioactive molecules that are able to interfere with cancer glycolysis.     

Integrin alpha9 (ITGA9) expression and epigenetic silencing in human breast tumors

Cancer metastasis is the major cause of cancer-associated death. Accordingly, identification of the regulatory mechanisms that control whether or not tumor cells become “directed walkers” is a crucial issue of cancer research. The deregulation of cell migration during cancer progression determines the capacity of tumor cells to escape from the primary tumors and invade adjacent tissues to finally form metastases. The ability to switch from a predominantly oxidative metabolism to glycolysis and the production of lactate even when oxygen is plentiful is a key characteristic of cancer cells. This metabolic switch, known as the Warburg effect, was first described in 1920s, and affected not only tumor cell growth but also tumor cell migration. In this review, we will focus on the recent studies on how cancer cell metabolism affects tumor cell migration and invasion. Understanding the new aspects on molecular mechanisms and signaling pathways controlling tumor cell migration is critical for development of therapeutic strategies for cancer patients.     

肿瘤生长的能量代谢特点及其临床应用

肿瘤细胞与人体正常细胞在代谢上有些不同,这主要体现在能量代谢和物质代谢上。肿瘤细胞能量代谢的特点表现在活跃地摄取葡萄糖和谷胺酰胺,进行有氧糖酵解(Warburg效应)。这种看上去很不经济的能量供给方式对肿瘤细胞却是必需的,它既为肿瘤细胞的不断生长提供能量,也为它们提供了生物合成的原料。肿瘤不同的代谢方式既是挑战也是机遇,弄清肿瘤细胞的代谢机制,对肿瘤早期诊断和靶向治疗具有重要意义。     

Distinct mechanisms of tumor invasion and metastasis

Most cancer deaths are caused by metastasis rather than the primary tumor. Cancer cells invade normal tissue as epithelial sheets or single cells by inducing expression of programs characteristic of developmental processes. Depending on their tissue of origin, cancer cells subsequently spread to distinct target organs where they seed secondary tumors (metastasis). Recent experimental evidence suggests that metastasis requires changes not only in cancer cells but also in the tumor microenvironment and in the metastatic target site. For example, a premetastatic niche is formed in target organs that attract cancer cells. Understanding the distinct mechanisms used by cancer cells to form metastasis will enable better patient evaluation and the design of innovative therapeutic approaches.     

Origins and evolution of cell phenotypes in breast tumors

This study presents a stochastic model that correlates genomic instability with tumor formation. The model describes the time- and space-variant volumetric concentrations of cancer cells of various phenotypes in a breast tumor. The cells of epithelial origin in the cancerous breast tissue are classified into four different phenotypes, normal epithelial cells and the grade 1, grade 2 and grade 3 cancer cell types with increasing potential for growth and invasion. Equations governing the time course of volumetric concentrations of cell phenotypes are derived by using the principle of conservation of mass. Cell migration into and from the stroma is taken into account. The transformations between cell phenotypes are due to genetic inheritance and chromosome aberrations. These transformations are assumed to be stochastic functions of the local cell concentration. The simulations of the model for planar geometry replicate the shapes of human breast tumors and capture the time history of tumor growth in animal models. Simulations point to transformation of tumor cell population from heterogeneous compositions to a single phenotype at advanced stages of invasive tumors. Systematic variations of model parameters in the computations indicate the important roles the migration capacity, proliferation rate, and phenotype transition probability play in tumor growth. The model developed provides realistic simulations for standard breast cancer therapies and can be used in the optimization studies of chemotherapy, radiotherapy, hormone therapy and emerging individualized therapies for cancer.     

The early pregnancy placenta foreshadows DNA methylation alterations of solid tumors

The placenta relies on phenotypes that are characteristic of cancer to successfully implant the embryo in the uterus during early pregnancy. Notably, it has to invade its host tissues, promote angiogenesis—while surviving hypoxia—, and escape the immune system. Similarities in DNA methylation patterns between the placenta and cancers suggest that common epigenetic mechanisms may be involved in regulating these behaviors. We show here that megabase-scale patterns of hypomethylation distinguish first from third trimester chorionic villi in the placenta, and that these patterns mirror those that distinguish many tumors from corresponding normal tissues. We confirmed these findings in villous cytotrophoblasts isolated from the placenta and identified a time window at the end of the first trimester, when these cells come into contact with maternal blood, as the likely time period for the methylome alterations. Furthermore, the large genomic regions affected by these patterns of hypomethylation encompass genes involved in pathways related to epithelial-mesenchymal transition, immune response, and inflammation. Analyses of expression profiles corresponding to genes in these hypomethylated regions in colon adenocarcinoma tumors point to networks of differentially expressed genes previously implicated in carcinogenesis and placentogenesis, where nuclear factor kappa B is a key hub. Taken together, our results suggest the existence of epigenetic switches involving large-scale changes of methylation in the placenta during pregnancy and in tumors during neoplastic transformation. The characterization of such epigenetic switches might lead to the identification of biomarkers and drug targets in oncology as well as in obstetrics and gynecology.     

How does a cell anchor and invade an organ?

In multicellular organisms, most cells are confined to a particular tissue. However, some cells invade organs during normal development and in diseases (e.g., angiogenesis and cancer). Recent studies reveal a fascinating step-by-step process in which specific vulval cells induce and attract a single gonadal cell to invade an epithelial tubular organ in order to connect the uterus to the vulva in C. elegans.     

Invadopodia and podosomes in tumor invasion

Cell migration through the extracellular matrix (ECM) is necessary for cancer cells to invade adjacent tissues and metastasize to an organ distant from primary tumors. Highly invasive carcinoma cells form ECM-degrading membrane protrusions called invadopodia. Tumor-associated macrophages have been shown to promote the migratory phenotypes of carcinoma cells, and macrophages are known to form podosomes, similar structures to invadopodia. However, the role of invadopodia and podosomes in vivo remains to be determined. In this paper, we propose a model for possible functions and interactions of invadopodia and podosomes in tumor invasion, based on observations that macrophage podosomes degrade ECM and that podosome formation is regulated by colony-stimulating factor-1 signaling.     

Competition and natural selection in a mathematical model of cancer

A malignant tumor is a dynamic amalgamation of various cell phenotypes, both cancerous (parenchyma) and healthy (stroma). These diverse cells compete over resources as well as cooperate to maintain tumor viability. Therefore, tumors are both an ecological community and an integrated tissue. An understanding of how natural selection operates in this unique ecological context should expose unappreciated vulnerabilities shared by all cancers. In this study I address natural selection’s role in tumor evolution by developing and exploring a mathematical model of a heterogenous primary neoplasm. The model is a system of nonlinear ordinary differential equations tracking the mass of up to two different parenchyma cell types, the mass of vascular endothelial cells from which new tumor blood vessels are built and the total length of tumor microvessels. Results predict the possibility of a hypertumor—a focus of aggressively reproducing parenchyma cells that invade and destroy part or all of the tumor, perhaps before it becomes a clinical entity. If this phenomenon occurs, then we should see examples of tumors that develop an aggressive histology but are paradoxically prone to extinction. Neuroblastoma, a common childhood cancer, may sometimes fit this pattern. In addition, this model suggests that parenchyma cell diversity can be maintained by a tissue-like integration of cells specialized to provide different services.     

Multiparametric Classification Links Tumor Microenvironments with Tumor Cell Phenotype

While it has been established that a number of microenvironment components can affect the likelihood of metastasis, the link between microenvironment and tumor cell phenotypes is poorly understood. Here we have examined microenvironment control over two different tumor cell motility phenotypes required for metastasis. By high-resolution multiphoton microscopy of mammary carcinoma in mice, we detected two phenotypes of motile tumor cells, different in locomotion speed. Only slower tumor cells exhibited protrusions with molecular, morphological, and functional characteristics associated with invadopodia. Each region in the primary tumor exhibited either fast- or slow-locomotion. To understand how the tumor microenvironment controls invadopodium formation and tumor cell locomotion, we systematically analyzed components of the microenvironment previously associated with cell invasion and migration. No single microenvironmental property was able to predict the locations of tumor cell phenotypes in the tumor if used in isolation or combined linearly. To solve this, we utilized the support vector machine (SVM) algorithm to classify phenotypes in a nonlinear fashion. This approach identified conditions that promoted either motility phenotype. We then demonstrated that varying one of the conditions may change tumor cell behavior only in a context-dependent manner. In addition, to establish the link between phenotypes and cell fates, we photoconverted and monitored the fate of tumor cells in different microenvironments, finding that only tumor cells in the invadopodium-rich microenvironments degraded extracellular matrix (ECM) and disseminated. The number of invadopodia positively correlated with degradation, while the inhibiting metalloproteases eliminated degradation and lung metastasis, consistent with a direct link among invadopodia, ECM degradation, and metastasis. We have detected and characterized two phenotypes of motile tumor cells in vivo, which occurred in spatially distinct microenvironments of primary tumors. We show how machine-learning analysis can classify heterogeneous microenvironments in vivo to enable prediction of motility phenotypes and tumor cell fate. The ability to predict the locations of tumor cell behavior leading to metastasis in breast cancer models may lead towards understanding the heterogeneity of response to treatment.     

Salmonella typhimurium specifically chemotax and proliferate in heterogeneous tumor tissue in vitro

Multi-drug resistance greatly limits the efficacy of conventional blood-born chemotherapeutics, which have limited ability to penetrate tumor tissue and are ineffective at killing quiescent cells far from tumor vasculature. Nonpathogenic, motile bacteria can overcome both of theses limitations. We hypothesize that the accumulation of S. typhimurium in tumors is controlled by two mechanisms: (1) chemotaxis towards compounds produced by quiescent cancer cells and (2) preferential growth within tumor tissue. We tested this hypothesis by quantifying the relative contributions of these mechanisms using the tumor cylindroid model, which mimics the microenvironments of in vivo tumors. Time-lapse fluorescence microscopy was used to measure the accumulation of GFP-labeled S. typhimurium into cylindroids of different size. Cylindroids larger than 500 microm in diameter contain quiescent cells, whereas cylindroids smaller than 500 microm do not. Spatio-temporal profiles of bacterial concentration were fit to a mathematical model to calculate two parameters that describe bacterial interaction with tumors: intratumoral bacterial growth, M, and intratumoral bacterial chemoattraction, K. It was observed that S. typhimurium is attracted to cylindroids and accumulate at long time points in the central region of large cylindroids. Both intratumoral bacterial growth and chemotaxis were significantly greater in large cylindroids, suggesting that quiescent cells secrete bacterial chemoattractants and the presence of necrotic and quiescent cells enable S. typhimurium to replicate in tumor tissue. In this study, several mechanisms of S. typhimurium accumulation in solid tumors have been quantified, which we believe is an important step in the development of bacterial-based therapeutics to target tumor quiescence.     

Cancer stem cells in glioblastoma — molecular signaling and therapeutic targeting

Glioblastomas (GBMs) are highly lethal primary brain tumors. Despite current therapeutic advances in other solid cancers, the treatment of these malignant gliomas remains essentially palliative. GBMs are extremely resistant to conventional radiation and chemotherapies. We and others have demonstrated that a highly tumorigenic subpopulation of cancer cells called GBM stem cells (GSCs) promotes therapeutic resistance. We also found that GSCs stimulate tumor angiogenesis by expressing elevated levels of VEGF and contribute to tumor growth, which has been translated into a useful therapeutic strategy in the treatment of recurrent or progressive GBMs. Furthermore, stem cell-like cancer cells (cancer stem cells) have been shown to promote metastasis. Although GBMs rarely metastasize beyond the central nervous system, these highly infiltrative cancers often invade into normal brain tissues preventing surgical resection, and GSCs display an aggressive invasive phenotype. These studies suggest that targeting GSCs may effectively reduce tumor recurrence and significantly improve GBM treatment. Recent studies indicate that cancer stem cells share core signaling pathways with normal somatic or embryonic stem cells, but also display critical distinctions that provide important clues into useful therapeutic targets. In this review, we summarize the current understanding and advances in glioma stem cell research, and discuss potential targeting strategies for future development of anti-GSC therapies.     

Evolutionary game theory elucidates the role of glycolysis in glioma progression and invasion

Abstract. Objectives: Tumour progression has been described as a sequence of traits or phenotypes that cells have to acquire if the neoplasm is to become an invasive and malignant cancer. Although genetic mutations that lead to these phenotypes are random, the process by which some of these mutations become successful and cells spread is influenced by tumour microenvironment and the presence of other cell phenotypes. It is thus likely that some phenotypes that are essential in tumour progression will emerge in the tumour population only with prior presence of other different phenotypes. Materials and methods: In this study, we use evolutionary game theory to analyse the interactions between three different tumour cell phenotypes defined by autonomous growth, anaerobic glycolysis, and cancer cell invasion. The model allows us to understand certain specific aspects of glioma progression such as the emergence of diffuse tumour cell invasion in low‐grade tumours. Results: We have found that the invasive phenotype is more likely to evolve after appearance of the glycolytic phenotype which would explain the ubiquitous presence of invasive growth in malignant tumours. Conclusions: The result suggests that therapies, which increase the fitness cost of switching to anaerobic glycolysis, might decrease probability of the emergence of more invasive phenotypes.     

Recurrent patterns of DNA copy number alterations in tumors reflect metabolic selection pressures

Copy number alteration (CNA) profiling of human tumors has revealed recurrent patterns of DNA amplifications and deletions across diverse cancer types. These patterns are suggestive of conserved selection pressures during tumor evolution but cannot be fully explained by known oncogenes and tumor suppressor genes. Using a pan‐cancer analysis of CNA data from patient tumors and experimental systems, here we show that principal component analysis‐defined CNA signatures are predictive of glycolytic phenotypes, including ¹⁸F‐fluorodeoxy‐glucose (FDG) avidity of patient tumors, and increased proliferation. The primary CNA signature is enriched for p53 mutations and is associated with glycolysis through coordinate amplification of glycolytic genes and other cancer‐linked metabolic enzymes. A pan‐cancer and cross‐species comparison of CNAs highlighted 26 consistently altered DNA regions, containing 11 enzymes in the glycolysis pathway in addition to known cancer‐driving genes. Furthermore, exogenous expression of hexokinase and enolase enzymes in an experimental immortalization system altered the subsequent copy number status of the corresponding endogenous loci, supporting the hypothesis that these metabolic genes act as drivers within the conserved CNA amplification regions. Taken together, these results demonstrate that metabolic stress acts as a selective pressure underlying the recurrent CNAs observed in human tumors, and further cast genomic instability as an enabling event in tumorigenesis and metabolic evolution.