A Computational Modeling Approach for Deriving Biomarkers to Predict Cancer Risk in Premalignant Disease

Evolution-Informed Strategies for Combating Drug Resistance in Cancer

The ever-changing nature of cancer poses the most difficult challenge oncologists face today. Cancer's remarkable adaptability has inspired many to work toward understanding the evolutionary dynamics that underlie this disease in hopes of learning new ways to fight it. Eco-evolutionary dynamics of a tumor are not accounted for in most standard treatment regimens, but exploiting them would help us combat treatment-resistant effectively. Here, we outline several notable efforts to exploit these dynamics and circumvent drug resistance in cancer.

Capturing the spatial and temporal dynamics of tumor stroma for on-chip optimization of microenvironmental targeting nanomedicine

Malignant cells grow in a complex microenvironment that plays a key role in cancer progression. The "dynamic reciprocity" existing between cancer cells and their microenvironment is involved in cancer differentiation, proliferation, invasion, metastasis, and drug response. Therefore, understanding the molecular mechanisms underlying the crosstalk between cancer cells and their surrounding tissue (i.e., tumor stroma) and how this interplay affects the disease progression is fundamental to design and validate novel nanotherapeutic approaches. As an important regulator of tumor progression, metastasis and therapy resistance, the extracellular matrix of tumors, the acellular component of the tumor microenvironment, has been identified as very promising target of anticancer treatment, revolutionizing the traditional therapeutic paradigm that sees the neoplastic cells as the preferential objective to fight cancer. To design and to validate such a target therapy, advanced 3D preclinical models are necessary to correctly mimic the complex, dynamic and heterogeneous tumor microenvironment. In addition, the recent advancement in microfluidic technology allows fine-tuning and controlling microenvironmental parameters in tissue-on-chip devices in order to emulate the in vivo conditions. In this review, after a brief description of the origin of tumor microenvironment heterogeneity, some examples of nanomedicine approaches targeting the tumor microenvironment have been reported. Further, how advanced 3D bioengineered tumor models coupled with a microfluidic device can improve the design and testing of anti-cancer nanomedicine targeting the tumor microenvironment has been discussed. We highlight that the presence of a dynamic extracellular matrix, able to capture the spatiotemporal heterogeneity of tumor stroma, is an indispensable requisite for tumor-on-chip model and nanomedicine testing.

Dietary stilbenes as modulators of specific miRNAs in prostate cancer

Accumulated experimental data have suggested that natural plant products may be effective miRNA-modulating chemopreventive and therapeutic agents. Dietary polyphenols such as flavonoids, stilbenes, and lignans, among others, have been intensively studied for their miRNA-mediated cardioprotective, antioxidant, anti-inflammatory and anticancer properties. The aim of this review is to outline known stilbene-regulated miRNAs in cancer, with a special focus on the interplay between various miRNAs and MTA1 signaling in prostate cancer. MTA1 is an epigenetic reader and an oncogenic transcription factor that is overexpressed in advanced prostate cancer and metastasis. Not surprisingly, miRNAs that are linked to MTA1 affect cancer progression and the metastatic potential of cells. Studies led to the identification of MTA1-associated pro-oncogenic miRNAs, which are regulated by stilbenes such as resveratrol and pterostilbene. Specifically, it has been shown that inhibition of the activity of the MTA1 regulated oncogenic miR-17 family of miRNAs, miR-22, and miR-34a by stilbenes leads to inhibition of prostatic hyperplasia and tumor progression in mice and reduction of proliferation, survival and invasion of prostate cancer cells in vitro. Taken together, these findings implicate the use of resveratrol and its analogs as an attractive miRNA-mediated chemopreventive and therapeutic strategy in prostate cancer and the use of circulating miRNAs as potential predictive biomarkers for clinical development.

Dietary Pterostilbene for MTA1-Targeted Interception in High-Risk Premalignant Prostate Cancer

Prostate cancer remains one of the most prevalent cancers in aging men. Active surveillance subpopulation of patients with prostate cancer includes men with varying cancer risk categories of precancerous disease due to prostatic intraepithelial neoplasia (PIN) heterogeneity. Identifying molecular alterations associated with PIN can provide preventable measures through finding novel pharmacologic targets for cancer interception. Targeted nutritional interception may prove to be the most appropriate chemoprevention for intermediate- and high-risk active surveillance patients. Here, we have generated two prostate-specific transgenic mouse models, one overexpressing MTA1 (R26MTA1 ) and the other overexpressing MTA1 on the background of Pten heterozygosity (R26MTA1 ; Pten+/f ), in which we examined the potential chemopreventive efficacy of dietary pterostilbene. We show that MTA1 promotes neoplastic transformation of prostate epithelial cells by activating cell proliferation and survival, leading to PIN development. Moreover, MTA1 cooperates with PTEN deficiency to accelerate PIN development by increasing cell proliferation and MTA1-associated signaling. Further, we show that mice fed with a pterostilbene-supplemented diet exhibited more favorable histopathology with decreased severity and number of PIN foci accompanied by reduced proliferation, angiogenesis, and inflammation concomitant to reduction in MTA1 and MTA1-associated CyclinD1, Notch2, and oncogenic miR-34a and miR-22 levels. PREVENTION RELEVANCE: Developing novel interceptive strategies for prostate cancer chemoprevention is a paramount goal in clinical oncology. We offer preclinical evidence for the potential of pterostilbene as a promising natural agent for MTA1-targeted interceptive strategy in future cancer prevention trials towards protecting select patients with prostate cancer under active surveillance from developing cancer.

Identifying key questions in the ecology and evolution of cancer

The application of evolutionary and ecological principles to cancer prevention and treatment, as well as recognizing cancer as a selection force in nature, has gained impetus over the last 50 years. Following the initial theoretical approaches that combined knowledge from interdisciplinary fields, it became clear that using the eco-evolutionary framework is of key importance to understand cancer. We are now at a pivotal point where accumulating evidence starts to steer the future directions of the discipline and allows us to underpin the key challenges that remain to be addressed. Here, we aim to assess current advancements in the field and to suggest future directions for research. First, we summarize cancer research areas that, so far, have assimilated ecological and evolutionary principles into their approaches and illustrate their key importance. Then, we assembled 33 experts and identified 84 key questions, organized around nine major themes, to pave the foundations for research to come. We highlight the urgent need for broadening the portfolio of research directions to stimulate novel approaches at the interface of oncology and ecological and evolutionary sciences. We conclude that progressive and efficient cross-disciplinary collaborations that draw on the expertise of the fields of ecology, evolution and cancer are essential in order to efficiently address current and future questions about cancer.

Inferring Tumor Proliferative Organization from Phylogenetic Tree Measures in a Computational Model

We use a computational modeling approach to explore whether it is possible to infer a solid tumor's cellular proliferative hierarchy under the assumptions of the cancer stem cell hypothesis and neutral evolution. We work towards inferring the symmetric division probability for cancer stem cells, since this is believed to be a key driver of progression and therapeutic response. Motivated by the advent of multiregion sampling and resulting opportunities to infer tumor evolutionary history, we focus on a suite of statistical measures of the phylogenetic trees resulting from the tumor's evolution in different regions of parameter space and through time. We find strikingly different patterns in these measures for changing symmetric division probability which hinge on the inclusion of spatial constraints. These results give us a starting point to begin stratifying tumors by this biological parameter and also generate a number of actionable clinical and biological hypotheses regarding changes during therapy, and through tumor evolutionary time. [Cancer; evolution; phylogenetics.].

Optimal Timing for Cancer Screening and Adaptive Surveillance Using Mathematical Modeling

Cancer screening and early detection efforts have been partially successful in reducing incidence and mortality, but many improvements are needed. Although current medical practice is informed by epidemiologic studies and experts, the decisions for guidelines are ultimately ad hoc. We propose here that quantitative optimization of protocols can potentially increase screening success and reduce overdiagnosis. Mathematical modeling of the stochastic process of cancer evolution can be used to derive and optimize the timing of clinical screens so that the probability is maximal that a patient is screened within a certain "window of opportunity" for intervention when early cancer development may be observable. Alternative to a strictly empirical approach or microsimulations of a multitude of possible scenarios, biologically based mechanistic modeling can be used for predicting when best to screen and begin adaptive surveillance. We introduce a methodology for optimizing screening, assessing potential risks, and quantifying associated costs to healthcare using multiscale models. As a case study in Barrett's esophagus, these methods were applied for a model of esophageal adenocarcinoma that was previously calibrated to U.S. cancer registry data. Optimal screening ages for patients with symptomatic gastroesophageal reflux disease were older (58 for men and 64 for women) than what is currently recommended (age > 50 years). These ages are in a cost-effective range to start screening and were independently validated by data used in current guidelines. Collectively, our framework captures critical aspects of cancer evolution within patients with Barrett's esophagus for a more personalized screening design. SIGNIFICANCE: This study demonstrates how mathematical modeling of cancer evolution can be used to optimize screening regimes, with the added potential to improve surveillance regimes. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/81/4/1123/F1.large.jpg.

Multicompartment modeling of protein shedding kinetics during vascularized tumor growth

Identification of protein biomarkers for cancer diagnosis and prognosis remains a critical unmet clinical need. A major reason is that the dynamic relationship between proliferating and necrotic cell populations during vascularized tumor growth, and the associated extra- and intra-cellular protein outflux from these populations into blood circulation remains poorly understood. Complementary to experimental efforts, mathematical approaches have been employed to effectively simulate the kinetics of detectable surface proteins (e.g., CA-125) shed into the bloodstream. However, existing models can be difficult to tune and may be unable to capture the dynamics of non-extracellular proteins, such as those shed from necrotic and apoptosing cells. The models may also fail to account for intra-tumoral spatial and microenvironmental heterogeneity. We present a new multi-compartment model to simulate heterogeneously vascularized growing tumors and the corresponding protein outflux. Model parameters can be tuned from histology data, including relative vascular volume, mean vessel diameter, and distance from vasculature to necrotic tissue. The model enables evaluating the difference in shedding rates between extra- and non-extracellular proteins from viable and necrosing cells as a function of heterogeneous vascularization. Simulation results indicate that under certain conditions it is possible for non-extracellular proteins to have superior outflux relative to extracellular proteins. This work contributes towards the goal of cancer biomarker identification by enabling simulation of protein shedding kinetics based on tumor tissue-specific characteristics. Ultimately, we anticipate that models like the one introduced herein will enable examining origins and circulating dynamics of candidate biomarkers, thus facilitating marker selection for validation studies.

When, why and how tumour clonal diversity predicts survival

The utility of intratumour heterogeneity as a prognostic biomarker is the subject of ongoing clinical investigation. However, the relationship between this marker and its clinical impact is mediated by an evolutionary process that is not well understood. Here, we employ a spatial computational model of tumour evolution to assess when, why and how intratumour heterogeneity can be used to forecast tumour growth rate and progression-free survival. We identify three conditions that can lead to a positive correlation between clonal diversity and subsequent growth rate: diversity is measured early in tumour development; selective sweeps are rare; and/or tumours vary in the rate at which they acquire driver mutations. Opposite conditions typically lead to negative correlation. In cohorts of tumours with diverse evolutionary parameters, we find that clonal diversity is a reliable predictor of both growth rate and progression-free survival. We thus offer explanations-grounded in evolutionary theory-for empirical findings in various cancers, including survival analyses reported in the recent TRACERx Renal study of clear-cell renal cell carcinoma. Our work informs the search for new prognostic biomarkers and contributes to the development of predictive oncology.

The prognostic value of a Methylome-based Malignancy Density Scoring System to predict recurrence risk in early-stage Lung Adenocarcinoma

Current NCCN guidelines do not recommend the use of adjuvant chemotherapy for stage IA lung adenocarcinoma patients with R0 surgery. However, 25% to 40% of patients with stage IA disease experience recurrence. Stratifying patients according to the recurrence risk may tailor adjuvant therapy and surveillance imaging for those with a higher risk. However, prognostic markers are often identified by comparing high-risk and low-risk cases which might introduce bias due to the widespread interpatient heterogeneity. Here, we developed a scoring system quantifying the degree of field cancerization in adjacent normal tissues and revealed its association with disease-free survival (DFS).

Methods: We recruited a cohort of 44 patients with resected stage IA lung adenocarcinoma who did not receive adjuvant therapy. Both tumor and adjacent normal tissues were obtained from each patient and subjected to capture-based targeted genomic and epigenomic profiling. A novel methylome-based scoring system namely malignancy density ratio (MD ratio) was developed based on 39 patients by comparing tumor and corresponding adjacent normal tissues of each patient. A MD score was then obtained by Wald statistics. The correlations of MD ratio, MD score, and genomic features with clinical outcome were investigated.

Results: Patients with a high-risk MD ratio showed a significantly shorter postsurgical DFS compared with those with a low-risk MD ratio (HR=4.47, P=0.01). The MD ratio was not associated with T stage (P=1), tumor cell fraction (P=0.748) nor inflammatory status (p=0.548). Patients with a high-risk MD score also demonstrated an inferior DFS (HR=4.69, P=0.039). In addition, multivariate analysis revealed EGFR 19 del (HR=5.39, P=0.012) and MD score (HR= 7.90, P=0.01) were independent prognostic markers.

Conclusion: The novel methylome-based scoring system, developed by comparing the signatures between tumor and corresponding adjacent normal tissues of individual patients, largely minimizes the bias of interpatient heterogeneity and reveals a robust prognostic value in patients with resected lung adenocarcinoma.

Precancer in ulcerative colitis: the role of the field effect and its clinical implications

Cumulative evidence indicates that a significant proportion of cancer evolution may occur before the development of histological abnormalities. While recent improvements in DNA sequencing technology have begun to reveal the presence of these early preneoplastic clones, the concept of 'premalignant field' was already introduced by Slaughter more than half a century ago. Also referred to as 'field effect', 'field defect' or 'field cancerization', these terms describe the phenomenon by which molecular alterations develop in normal-appearing tissue and expand to form premalignant patches with the potential to progress to dysplasia and cancer. Field effects have been well-characterized in ulcerative colitis, an inflammatory bowel disease that increases the risk of colorectal cancer. The study of the molecular alterations that define these fields is informative of mechanisms of tumor initiation and progression and has provided potential targets for early cancer detection. Herein, we summarize the current knowledge about the molecular alterations that comprise the field effect in ulcerative colitis and the clinical utility of these fields for cancer screening and prevention.

An evolutionary perspective on field cancerization

Tumorigenesis begins long before the growth of a clinically detectable lesion and, indeed, even before any of the usual morphological correlates of pre-malignancy are recognizable. Field cancerization, which is the replacement of the normal cell population by a cancer-primed cell population that may show no morphological change, is now recognized to underlie the development of many types of cancer, including the common carcinomas of the lung, colon, skin, prostate and bladder. Field cancerization is the consequence of the evolution of somatic cells in the body that results in cells that carry some but not all phenotypes required for malignancy. Here, we review the evidence of field cancerization across organs and examine the biological mechanisms that drive the evolutionary process that results in field creation. We discuss the clinical implications, principally, how measurements of the cancerized field could improve cancer risk prediction in patients with pre-malignant disease.

Evolution of Premalignant Disease

Where does cancer come from? Although the cell-of-origin is difficult to pinpoint, cancer clones harbor information about their clonal ancestries. In an effort to find cells before they evolve into a life-threatening cancer, physicians currently diagnose premalignant diseases at frequencies that substantially exceed those of clinical cancers. Cancer risk prediction relies on our ability to distinguish between which premalignant features will lead to cancer mortality and which are characteristic of inconsequential disease. Here, we review the evolution of cancer from premalignant disease, and discuss the concept that even phenotypically normal cell progenies inherently gain more malignant potential with age. We describe the hurdles of prognosticating cancer risk in premalignant disease by making reference to the underlying continuous and multivariate natures of genotypes and phenotypes and the particular challenge inherent in defining a cell lineage as "cancerized."

Spatial Measures of Genetic Heterogeneity During Carcinogenesis

In this work we explore the temporal dynamics of spatial heterogeneity during the process of tumorigenesis from healthy tissue. We utilize a spatial stochastic model of mutation accumulation and clonal expansion in a structured tissue to describe this process. Under a two-step tumorigenesis model, we first derive estimates of a non-spatial measure of diversity: Simpson's Index, which is the probability that two individuals sampled at random from the population are identical, in the premalignant population. We next analyze two new measures of spatial population heterogeneity. In particular we study the typical length scale of genetic heterogeneity during the carcinogenesis process and estimate the extent of a surrounding premalignant clone given a clinical observation of a premalignant point biopsy. This evolutionary framework contributes to a growing literature focused on developing a better understanding of the spatial population dynamics of cancer initiation and progression.