Tumors are highly dynamic ecosystems in which diverse cancer cell subpopulations compete for space and resources. These complex, often non-linear interactions govern continuous spatial and temporal changes in the size and phenotypic properties of these subpopulations. Because intra-tumoral blood flow is often chaotic, competition for resources may be a critical selection factor in progression and prognosis. Here, we quantify resource competition using 3D spheroid cultures with MDA-MB-231 and MCF-7 breast cancer cells. We hypothesized that MCF-7 cells, which primarily rely on efficient aerobic glucose metabolism, would dominate the population under normal pH and low glucose conditions; and MDA-MB-231 cells, which exhibit high levels of glycolytic metabolism, would dominate under low pH and high glucose conditions. In spheroids with single populations, MCF-7 cells exhibited equal or superior intrinsic growth rates (density-independent measure of success) and carrying capacities (density-dependent measure of success) when compared to MDA-MB-231 cells under all pH and nutrient conditions. Despite these advantages, when grown together, MCF-7 cells do not always outcompete MDA-MB-231 cells. MDA-MB-231 cells outcompete MCF-7 cells in low glucose conditions and coexistence is achieved in low pH conditions. Under all conditions, MDA-MB-231 has a stronger competitive effect (frequency-dependent interaction) on MCF-7 cells than vice-versa. This, and the inability of growth rate or carrying capacity when grown individually to predict the outcome of competition, suggests a reliance on frequency-dependent interactions and the need for competition assays. We frame these results in a game-theoretic (frequency-dependent) model of cancer cell interactions and conclude that competition assays can demonstrate critical density-independent, density-dependent and frequency-dependent interactions that likely contribute to in vivo outcomes.

We identify critical conserved and mutated genes through a theoretical model linking a gene’s fitness contribution to its observed mutational frequency in a clinical cohort. “Passenger” gene mutations do not alter fitness and have mutational frequencies determined by gene size and the mutation rate. Driver mutations, which increase fitness (and proliferation), are observed more frequently than expected. Non-synonymous mutations in essential genes reduce fitness and are eliminated by natural selection resulting in lower prevalence than expected. We apply this “evolutionary triage” principle to TCGA data from EGFR-mutant, KRAS-mutant, and NEK (non-EGFR/KRAS) lung adenocarcinomas. We find frequent overlap of evolutionarily selected non-synonymous gene mutations among the subtypes suggesting enrichment for adaptations to common local tissue selection forces. Overlap of conserved genes in the LUAD subtypes is rare suggesting negative evolutionary selection is strongly dependent on initiating mutational events during carcinogenesis. Highly expressed genes are more likely to be conserved and significant changes in expression (>20% increased/decreased) are common in genes with evolutionarily selected mutations but not in conserved genes. EGFR-mut cancers have fewer average mutations (89) than KRAS-mut (228) and NEK (313). Subtype-specific variation in conserved and mutated genes identify critical molecular components in cell signaling, extracellular matrix remodeling, and membrane transporters. These findings demonstrate subtype-specific patterns of co-adaptations between the defining driver mutation and somatically conserved genes as well as novel insights into epigenetic versus genetic contributions to cancer evolution.

Tumors are highly dynamic ecosystems in which diverse cancer cell subpopulations compete for space and resources. These complex, often non-linear interactions govern continuous spatial and temporal changes in the size and phenotypic properties of these subpopulations. Because intra-tumoral blood flow is often chaotic, competition for resources may be a critical selection factor in progression and prognosis. Here, we quantify resource competition using 3D spheroid cultures with MDA-MB-231 and MCF-7 breast cancer cells. We hypothesized that MCF-7 cells, which primarily rely on efficient aerobic glucose metabolism, would dominate the population under normal pH and low glucose conditions; and MDA-MB-231 cells, which exhibit high levels of glycolytic metabolism, would dominate under low pH and high glucose conditions. In spheroids with single populations, MCF-7 cells exhibited equal or superior intrinsic growth rates (density-independent measure of success) and carrying capacities (density-dependent measure of success) when compared to MDA-MB-231 cells under all pH and nutrient conditions. Despite these advantages, when grown together, MCF-7 cells do not always outcompete MDA-MB-231 cells. MDA-MB-231 cells outcompete MCF-7 cells in low glucose conditions and coexistence is achieved in low pH conditions. Under all conditions, MDA-MB-231 has a stronger competitive effect (frequency-dependent interaction) on MCF-7 cells than vice-versa. This, and the inability of growth rate or carrying capacity when grown individually to predict the outcome of competition, suggests a reliance on frequency-dependent interactions and the need for competition assays. We frame these results in a game-theoretic (frequency-dependent) model of cancer cell interactions and conclude that competition assays can demonstrate critical density-independent, density-dependent and frequency-dependent interactions that likely contribute to in vivo outcomes. Highlights: Demonstrate how mixed-culture spheroids must be used to characterize competition between two cancer cell lines. Competition alters growth dynamics of cancer cells. Competition growth models can be used to quantify density-independent, density-dependent and frequency-dependent effects on competition. Competition affects tumor progression and structure, making it key to understanding tumor development and evolution.