Examining the Relationship between Pre-Malignant Breast Lesions, Carcinogenesis and Tumor Evolution in the Mammary Epithelium Using an Agent-Based Model

Introduction Breast cancer, the product of numerous rare mutational events that occur over an extended time period, presents numerous challenges to investigators interested in studying the transformation from normal breast epithelium to malignancy using traditional laboratory methods, particularly with respect to characterizing transitional and pre-malignant states. Dynamic computational modeling can provide insight into these pathophysiological dynamics, and as such we use a previously validated agent-based computational model of the mammary epithelium (the DEABM) to investigate the probabilistic mechanisms by which normal populations of ductal cells could transform into states replicating features of both pre-malignant breast lesions and a diverse set of breast cancer subtypes. Methods The DEABM consists of simulated cellular populations governed by algorithms based on accepted and previously published cellular mechanisms. Cells respond to hormones, undergo mitosis, apoptosis and cellular differentiation. Heritable mutations to 12 genes prominently implicated in breast cancer are acquired via a probabilistic mechanism. 3000 simulations of the 40-year period of menstrual cycling were run in wild-type (WT) and BRCA1-mutated groups. Simulations were analyzed by development of hyperplastic states, incidence of malignancy, hormone receptor and HER-2 status, frequency of mutation to particular genes, and whether mutations were early events in carcinogenesis. Results Cancer incidence in WT (2.6%) and BRCA1-mutated (45.9%) populations closely matched published epidemiologic rates. Hormone receptor expression profiles in both WT and BRCA groups also closely matched epidemiologic data. Hyperplastic populations carried more mutations than normal populations and mutations were similar to early mutations found in ER+ tumors (telomerase, E-cadherin, TGFB, RUNX3, p < .01). ER-tumors carried significantly more mutations and carried more early mutations in BRCA1, c-MYC and genes associated with epithelial-mesenchymal transition. Conclusions The DEABM generates diverse tumors that express tumor markers consistent with epidemiologic data. The DEABM also generates non-invasive, hyperplastic populations, analogous to atypia or ductal carcinoma in situ (DCIS), via mutations to genes known to be present in hyperplastic lesions and as early mutations in breast cancers. The results demonstrate that agent-based models are well-suited to studying tumor evolution through stages of carcinogenesis and have the potential to be used to develop prevention and treatment strategies.