Microenvironmental factors in cell segregation and heterogeneity in breast cancer development

We analyzed a quantitative model that describes the epigenetic dynamics during the growth and evolution of an avascular tumor. A gene regulatory network (GRN) formed by a set of ten genes that are believed to play an important role in breast cancer development was kinetically coupled to the microenvironmental agents: glucose, estrogens and oxygen. The dynamics of spontaneous mutations was described by a Yule-Furry master equation whose solution represents the probability that a given cell in the tissue undergoes a certain number of mutations at a given time. We assumed that mutations rate is modified by nutrients spatial gradients. The tumor mass was grown by means of a cellular automata supplemented with a set of reaction diffusion equations that described the transport of the microenvironmental agents. By analyzing the epigenetic states space described by the GRN dynamics, we found three attractors that were identified with the cellular epigenetic states: normal, precancer and cancer. For two-dimensional (2D) and three-dimensional (3D) tumors we calculated the spatial distributions of the following quantities: (i) number of mutations, (ii) mutations of each gene and, (iii) phenotypes. Using estrogens as the principal microenvironmental agent that regulates cells proliferation process, we obtained the tumor shapes for different values of the estrogen consumption and supply rates. It was found that he majority of mutations occurred in cells that were located close to the 2D tumor perimeter or close to the 3D tumor surface. Also It was found that the occurrence of different phenotypes in the tumor are controlled by the levels of estrogen concentration since they can change the individual cell threshold and gene expression levels. All the results were consistently observed for 2D and 3D tumors.