Running title: CT-based subtypes of pancreatic ductal adenocarcinoma. Keywords Statement of translational relevanceThere are currently limited methods to stratify patients with pancreatic ductal adenocarcinoma (PDAC) into prognostic groups based on disease biology. We focused on the morphological features of the disease on computed tomography (CT) scans and correlated these features with genomic, pathological, and clinical data. We found that tumors with a distinct interface in relation to the surrounding parenchyma (called high delta tumors) had more aggressive biological features than tumors without a distinct interface (called low delta tumors). Patients with high delta tumors were more likely to develop early distant metastasis and die faster than those with low delta tumors. Mathematical modeling suggested that stromal elements strongly influenced the morphological patterns seen on CT scans. These findings indicate that a universally available diagnostic test can be used to interrogate the biology of this deadly disease, providing a means to stratify patients at diagnosis and aiding the design of future clinical trials.

Feedback regulation of cell lineage progression plays an important role in tissue size homeostasis, but whether such feedback also plays an important role in tissue morphogenesis has yet to be explored. Here we use mathematical modeling to show that a particular feedback architecture in which both positive and negative diffusible signals act on stem and/or progenitor cells leads to the appearance of bistable or bi-modal growth behaviors, ultrasensitivity to external growth cues, local growth-driven budding, self-sustaining elongation, and the triggering of self-organization in the form of lamellar fingers. Such behaviors arise not through regulation of cell cycle speeds, but through the control of stem or progenitor self-renewal. Even though the spatial patterns that arise in this setting are the result of interactions between diffusible factors with antagonistic effects, morphogenesis is not the consequence of Turing-type instabilities.

Mutational activation of the BRAF proto-oncogene in melanocytes reliably produces benign nevi (pigmented ‘moles’), yet the same change is the most common driver mutation in melanoma. The reason nevi stop growing, and do not progress to melanoma, is widely attributed to a cell-autonomous process of ‘oncogene-induced senescence’. Using a mouse model of Braf-driven nevus formation, analyzing both proliferative dynamics and single-cell gene expression, we found no evidence that nevus cells are senescent, either compared with other skin cells, or other melanocytes. We also found that nevus size distributions could not be fit by any simple cell-autonomous model of growth arrest, yet were easily fit by models based on collective cell behavior, for example in which arresting cells release an arrest-promoting factor. We suggest that nevus growth arrest is more likely related to the cell interactions that mediate size control in normal tissues, than to any cell-autonomous, ‘oncogene-induced’ program of senescence.

Glioblastoma (GBM), the most aggressive brain tumor in human patients, is decidedly heterogeneous and highly vascularized. Glioma stem/initiating cells (GSC) are found to play a crucial role by increasing cancer aggressiveness and promoting resistance to therapy. Recently, crosstalk between GSC and vascular endothelial cells has been shown to significantly promote GSC self-renewal and tumor progression. Further, GSC also transdifferentiate into bona-fide vascular endothelial cells (GEC), which inherit mutations present in GSC and are resistant to traditional anti-angiogenic therapies. Here we use 3D mathematical modeling to investigate GBM progression and response to therapy. The model predicted that GSC drive invasive fingering and that GEC spontaneously form a network within the hypoxic core, consistent with published experimental findings. Standard-of-care treatments using DNA-targeted therapy (radiation/chemo) together with anti-angiogenic therapies, reduced GBM tumor size but increased invasiveness. Anti-GEC treatments blocked the GEC support of GSC and reduced tumor size but led to increased invasiveness. Anti-GSC therapies that promote differentiation or disturb the stem cell niche effectively reduced tumor invasiveness and size, but were ultimately limited in reducing tumor size because GEC maintain GSC. Our study suggests that a combinatorial regimen targeting the vasculature, GSC, and GEC, using drugs already approved by the FDA, can reduce both tumor size and invasiveness and could lead to tumor eradication.

We develop a three-dimensional multispecies mathematical model to simulate the growth of colon cancer organoids containing stem, progenitor and terminally differentiated cells, as a model of early (prevascular) tumor growth. Stem cells (SCs) secrete short-range self-renewal promoters (e.g., Wnt) and their long-range inhibitors (e.g., Dkk) and proliferate slowly. Committed progenitor (CP) cells proliferate more rapidly and differentiate to produce post-mitotic terminally differentiated cells that release differentiation promoters, forming negative feedback loops on SC and CP self-renewal. We demonstrate that SCs play a central role in normal and cancer colon organoids. Spatial patterning of the SC self-renewal promoter gives rise to SC clusters, which mimic stem cell niches, around the organoid surface, and drive the development of invasive fingers. We also study the effects of externally applied signaling factors. Applying bone morphogenic proteins, which inhibit SC and CP self-renewal, reduces invasiveness and organoid size. Applying hepatocyte growth factor, which enhances SC self-renewal, produces larger sizes and enhances finger development at low concentrations but suppresses fingers at high concentrations. These results are consistent with recent experiments on colon organoids. Because many cancers are hierarchically organized and are subject to feedback regulation similar to that in normal tissues, our results suggest that in cancer, control of cancer stem cell self-renewal should influence the size and shape in similar ways, thereby opening the door to novel therapies.

Objective In glioblastoma, the crosstalk between vascular endothelial cells (VECs) and glioma stem cells (GSCs) has been shown to enhance tumor growth. We propose a multiscale mathematical model to study this mechanism, explore tumor growth under various initial and microenvironmental conditions and investigate the effects of blocking this crosstalk. Methods We develop a hybrid continuum-discrete model of highly organized, vascularized tumors. VEC-GSC crosstalk is modeled via VEGF production by tumor cells and by secretion of soluble factors by VECs that promote GSC self-renewal and proliferation. Results VEC-GSC crosstalk increases both tumor size and GSC fraction by enhancing GSC activity and neovascular development. VEGF promotes vessel formation, and larger VEGF sources typically increase vessel numbers, which enhances tumor growth and stabilizes the tumor shape. Increasing the initial GSC fraction has a similar effect. Partially disrupting the crosstalk by blocking VEC secretion of GSC promoters reduces tumor size but does not increase invasiveness, which is in contrast to anti-angiogenic therapies, which reduce tumor size but may significantly increase tumor invasiveness. Significance Multiscale modeling supports the targeting of VEC-GSC crosstalk as a promising approach for cancer therapy.

Blebs are pressure-driven cell protrusions implicated in cellular functions such as cell division, apoptosis, and cell motility, including motility of protease-inhibited cancer cells. Because of their mechanical nature, blebs inform us about general cell-surface mechanics, including membrane dynamics, pressure propagation throughout the cytoplasm, and the architecture and dynamics of the actin cortex. Mathematical models including detailed fluid dynamics have previously been used to understand bleb expansion. Here, we develop mathematical models in two and three dimensions on longer timescales that recapitulate the full bleb life cycle, including both expansion and healing by cortex reformation, in terms of experimentally accessible biophysical parameters such as myosin contractility, osmotic pressure, and turnover of actin and ezrin. The model provides conditions under which blebbing occurs, and naturally gives rise to traveling blebs. The model predicts conditions under which blebs travel or remain stationary, as well as the bleb traveling velocity, a quantity that has remained elusive in previous models. As previous studies have used blebs as reporters of membrane tension and pressure dynamics within the cell, we have used our system to investigate various pressure equilibration models and dynamic, nonuniform membrane tension to account for the shape of a traveling bleb. We also find that traveling blebs tend to expand in all directions unless otherwise constrained. One possible constraint could be provided by spatial heterogeneity in, for example, adhesion density.

Experiments on tumor spheroids have shown that compressive stress from their environment can reversibly decrease tumor expansion rates and final sizes. Stress release experiments show that nonuniform anisotropic elastic stresses can be distributed throughout. The elastic stresses are maintained by structural proteins and adhesive molecules, and can be actively relaxed by a variety of biophysical processes. In this paper, we present a new continuum model to investigate how the growth-induced elastic stresses and active stress relaxation, in conjunction with cell size control feedback machinery, regulate the cell density and stress distributions within growing tumors as well as the tumor sizes in the presence of external physical confinement and gradients of growth-promoting chemical fields. We introduce an adaptive reference map that relates the current position with the reference position but adapts to the current position in the Eulerian frame (lab coordinates) via relaxation. This type of stress relaxation is similar to but simpler than the classical Maxwell model of viscoelasticity in its formulation. By fitting the model to experimental data from two independent studies of tumor spheroid growth and their cell density distributions, treating the tumors as incompressible, neo-Hookean elastic materials, we find that the rates of stress relaxation of tumor tissues can be comparable to volumetric growth rates. Our study provides insight on how the biophysical properties of the tumor and host microenvironment, mechanical feedback control and diffusion-limited differential growth act in concert to regulate spatial patterns of stress and growth. When the tumor is stiffer than the host, our model predicts tumors are more able to change their size and mechanical state autonomously, which may help to explain why increased tumor stiffness is an established hallmark of malignant tumors.