Computational modelling of perivascular-niche dynamics for the optimization of treatment schedules for glioblastoma.

Glioblastoma stem-like cells dynamically transition between a chemoradiation-resistant state and a chemoradiation-sensitive state. However, physical barriers in the tumour microenvironment restrict the delivery of chemotherapy to tumour compartments that are distant from blood vessels. Here, we show that a massively parallel computational model of the spatiotemporal dynamics of the perivascular niche that incorporates glioblastoma stem-like cells and differentiated tumour cells as well as relevant tissue-level phenomena can be used to optimize the administration schedules of concurrent radiation and temozolomide-the standard-of-care treatment for glioblastoma. In mice with platelet-derived growth factor (PDGF)-driven glioblastoma, the model-optimized treatment schedule increased the survival of the animals. For standard radiation fractionation in patients, the model predicts that chemotherapy may be optimally administered about one hour before radiation treatment. Computational models of the spatiotemporal dynamics of the tumour microenvironment could be used to predict tumour responses to a broader range of treatments and to optimize treatment regimens.

Mathematical Modeling Links Pregnancy-Associated Changes and Breast Cancer Risk.

Recent debate has concentrated on the contribution of bad luck to cancer development. The tight correlation between the number of tissue-specific stem cell divisions and cancer risk of the same tissue suggests that bad luck has an important role to play in tumor development, but the full extent of this contribution remains an open question. Improved understanding of the interplay between extrinsic and intrinsic factors at the molecular level is one promising route to identifying the limits on extrinsic control of tumor initiation, which is highly relevant to cancer prevention. Here, we use a simple mathematical model to show that recent data on the variation in numbers of breast epithelial cells with progenitor features due to pregnancy are sufficient to explain the known protective effect of full-term pregnancy in early adulthood for estrogen receptor-positive (ER+) breast cancer later in life. Our work provides a mechanism for this previously ill-understood effect and illuminates the complex influence of extrinsic factors at the molecular level in breast cancer. These findings represent an important contribution to the ongoing research into the role of bad luck in human tumorigenesis. Cancer Res; 77(11); 2800-9. ©2017 AACR.

Most human non-GCIMP glioblastoma subtypes evolve from a common proneural-like precursor glioma.

To understand the relationships between the non-GCIMP glioblastoma (GBM) subgroups, we performed mathematical modeling to predict the temporal sequence of driver events during tumorigenesis. The most common order of evolutionary events is 1) chromosome (chr) 7 gain and chr10 loss, followed by 2) CDKN2A loss and/or TP53 mutation, and 3) alterations canonical for specific subtypes. We then developed a computational methodology to identify drivers of broad copy number changes, identifying PDGFA (chr7) and PTEN (chr10) as driving initial nondisjunction events. These predictions were validated using mouse modeling, showing that PDGFA is sufficient to induce proneural-like gliomas and that additional NF1 loss converts proneural to the mesenchymal subtype. Our findings suggest that most non-GCIMP mesenchymal GBMs arise as, and evolve from, a proneural-like precursor.

Genetic and phenotypic diversity in breast tumor metastases.

Metastatic disease is the main cause of cancer-related mortality due to almost universal therapeutic resistance. Despite its high clinical relevance, our knowledge of how cancer cell populations change during metastatic progression is limited. Here, we investigated intratumor genetic and phenotypic heterogeneity during metastatic progression of breast cancer. We analyzed cellular genotypes and phenotypes at the single cell level by performing immunoFISH in intact tissue sections of distant metastatic tumors from rapid autopsy cases and from primary tumors and matched lymph node metastases collected before systemic therapy. We calculated the Shannon index of intratumor diversity in all cancer cells and within phenotypically distinct cell populations. We found that the extent of intratumor genetic diversity was similar regardless of the chromosomal region analyzed, implying that it may reflect an inherent property of the tumors. We observed that genetic diversity was highest in distant metastases and was generally concordant across lesions within the same patient, whereas treatment-naïve primary tumors and matched lymph node metastases were frequently genetically more divergent. In contrast, cellular phenotypes were more discordant between distant metastases than primary tumors and matched lymph node metastases. Diversity for 8q24 was consistently higher in HER2(+) tumors compared with other subtypes and in metastases of triple-negative tumors relative to primary sites. We conclude that our integrative method that couples ecologic models with experimental data in human tissue samples could be used for the improved prognostication of patients with cancer and for the design of more effective therapies for progressive disease.

Inference of tumor evolution during chemotherapy by computational modeling and in situ analysis of genetic and phenotypic cellular diversity.

Cancer therapy exerts a strong selection pressure that shapes tumor evolution, yet our knowledge of how tumors change during treatment is limited. Here, we report the analysis of cellular heterogeneity for genetic and phenotypic features and their spatial distribution in breast tumors pre- and post-neoadjuvant chemotherapy. We found that intratumor genetic diversity was tumor-subtype specific, and it did not change during treatment in tumors with partial or no response. However, lower pretreatment genetic diversity was significantly associated with pathologic complete response. In contrast, phenotypic diversity was different between pre- and posttreatment samples. We also observed significant changes in the spatial distribution of cells with distinct genetic and phenotypic features. We used these experimental data to develop a stochastic computational model to infer tumor growth patterns and evolutionary dynamics. Our results highlight the importance of integrated analysis of genotypes and phenotypes of single cells in intact tissues to predict tumor evolution.

Molecular profiling of human mammary gland links breast cancer risk to a p27(+) cell population with progenitor characteristics.

Early full-term pregnancy is one of the most effective natural protections against breast cancer. To investigate this effect, we have characterized the global gene expression and epigenetic profiles of multiple cell types from normal breast tissue of nulliparous and parous women and carriers of BRCA1 or BRCA2 mutations. We found significant differences in CD44(+) progenitor cells, where the levels of many stem cell-related genes and pathways, including the cell-cycle regulator p27, are lower in parous women without BRCA1/BRCA2 mutations. We also noted a significant reduction in the frequency of CD44(+)p27(+) cells in parous women and showed, using explant cultures, that parity-related signaling pathways play a role in regulating the number of p27(+) cells and their proliferation. Our results suggest that pathways controlling p27(+) mammary epithelial cells and the numbers of these cells relate to breast cancer risk and can be explored for cancer risk assessment and prevention.

A mathematical methodology for determining the temporal order of pathway alterations arising during gliomagenesis.

Human cancer is caused by the accumulation of genetic alterations in cells. Of special importance are changes that occur early during malignant transformation because they may result in oncogene addiction and thus represent promising targets for therapeutic intervention. We have previously described a computational approach, called Retracing the Evolutionary Steps in Cancer (RESIC), to determine the temporal sequence of genetic alterations during tumorigenesis from cross-sectional genomic data of tumors at their fully transformed stage. Since alterations within a set of genes belonging to a particular signaling pathway may have similar or equivalent effects, we applied a pathway-based systems biology approach to the RESIC methodology. This method was used to determine whether alterations of specific pathways develop early or late during malignant transformation. When applied to primary glioblastoma (GBM) copy number data from The Cancer Genome Atlas (TCGA) project, RESIC identified a temporal order of pathway alterations consistent with the order of events in secondary GBMs. We then further subdivided the samples into the four main GBM subtypes and determined the relative contributions of each subtype to the overall results: we found that the overall ordering applied for the proneural subtype but differed for mesenchymal samples. The temporal sequence of events could not be identified for neural and classical subtypes, possibly due to a limited number of samples. Moreover, for samples of the proneural subtype, we detected two distinct temporal sequences of events: (i) RAS pathway activation was followed by TP53 inactivation and finally PI3K2 activation, and (ii) RAS activation preceded only AKT activation. This extension of the RESIC methodology provides an evolutionary mathematical approach to identify the temporal sequence of pathway changes driving tumorigenesis and may be useful in guiding the understanding of signaling rearrangements in cancer development.

The probable cell of origin of NF1- and PDGF-driven glioblastomas.

Primary glioblastomas are subdivided into several molecular subtypes. There is an ongoing debate over the cell of origin for these tumor types where some suggest a progenitor while others argue for a stem cell origin. Even within the same molecular subgroup, and using lineage tracing in mouse models, different groups have reached different conclusions. We addressed this problem from a combined mathematical modeling and experimental standpoint. We designed a novel mathematical framework to identify the most likely cells of origin of two glioma subtypes. Our mathematical model of the unperturbed in vivo system predicts that if a genetic event contributing to tumor initiation imparts symmetric self-renewing cell division (such as PDGF overexpression), then the cell of origin is a transit amplifier. Otherwise, the initiating mutations arise in stem cells. The mathematical framework was validated with the RCAS/tv-a system of somatic gene transfer in mice. We demonstrated that PDGF-induced gliomas can be derived from GFAP-expressing cells of the subventricular zone or the cortex (reactive astrocytes), thus validating the predictions of our mathematical model. This interdisciplinary approach allowed us to determine the likelihood that individual cell types serve as the cells of origin of gliomas in an unperturbed system.

A mathematical framework to determine the temporal sequence of somatic genetic events in cancer.

Human cancer is caused by the accumulation of genetic alterations in cells. Of special importance are changes that occur early during malignant transformation because they may result in oncogene addiction and represent promising targets for therapeutic intervention. Here we describe a computational approach, called Retracing the Evolutionary Steps in Cancer (RESIC), to deduce the temporal sequence of genetic events during tumorigenesis from cross-sectional genomic data of tumors at their fully transformed stage. When applied to a dataset of 70 advanced colorectal cancers, our algorithm accurately predicts the sequence of APC, KRAS, and TP53 mutations previously defined by analyzing tumors at different stages of colon cancer formation. We further validate the method with glioblastoma and leukemia sample data and then apply it to complex integrated genomics databases, finding that high-level EGFR amplification appears to be a late event in primary glioblastomas. RESIC represents the first evolutionary mathematical approach to identify the temporal sequence of mutations driving tumorigenesis and may be useful to guide the validation of candidate genes emerging from cancer genome surveys.