Search

Background A vestibular schwannoma (VS) is a relatively rare, benign tumour of the eighth cranial nerve, often involving alterations to the gene NF2. Previous mathematical models of schwannoma incidence have not attempted to account for alterations in specific genes, and could not distinguish between nonsense mutations and loss of heterozygosity (LOH). Methods Here, we present a mechanistic approach to modelling initiation and malignant transformation in schwannoma. Each parameter is associated with a specific gene or mechanism operative in Schwann cells, and can be determined by combining incidence data with empirical frequencies of pathogenic variants and LOH. Results This results in new estimates for the base-pair mutation rate u = 4.48 × 10−10 and the rate of LOH = 2.03 × 10−6/yr in Schwann cells. In addition to new parameter estimates, we extend the approach to estimate the risk of both spontaneous and radiation-induced malignant transformation. Discussion We conclude that radiotherapy is likely to have a negligible excess risk of malignancy for sporadic VS, with a possible exception of rapidly growing tumours.

Quantifying evolutionary dynamics of cancer initiation and progression can provide insights into more effective strategies of early detection and treatment. Here we develop a mathematical model of colorectal cancer initiation through inactivation of two tumor suppressor genes and activation of one oncogene, accounting for the well-known path to colorectal cancer through loss of tumor suppressors APC and TP53, and gain of the KRAS oncogene. In the model, we allow mutations to occur in any order, leading to a complex network of incomplete mutational genotypes on the way to colorectal cancer. We parametrize the model using experimentally measured parameter values, many of them only recently available, and compare its predictions to epidemiological data on colorectal cancer incidence. We find that the reported incidence of colorectal cancer can be recovered using a mathematical model of colorectal cancer initiation together with experimentally measured mutation rates in colorectal tissues and proliferation rates of premalignant lesions. We demonstrate that the order of driver events in colorectal cancer is determined by the combined effect of the rates at which driver genes are mutated and the fitness effects they provide. Our results imply that there may not be significant immune suppression of untreated benign and malignant colorectal lesions.

Recently available cancer sequencing data have revealed a complex view of the cancer genome containing a multitude of mutations, including drivers responsible for cancer progression and neutral passengers. Measuring selection in cancer and distinguishing drivers from passengers have important implications for development of novel treatment strategies. It has recently been argued that a third of cancers are evolving neutrally, as their mutational frequency spectrum follows a 1/f power law expected from neutral evolution in a particular intermediate frequency range. We study a stochastic model of cancer evolution and derive a formula for the probability distribution of the cancer cell frequency of a subclonal driver, demonstrating that driver frequency is biased towards 0 and 1. We show that it is difficult to capture a driver mutation at an intermediate frequency, and thus the calling of neutrality due to a lack of such driver will significantly overestimate the number of neutrally evolving tumors. Our approach provides quantification of the validity of the 1/f statistic across the entire range of relevant parameter values. We also show that our conclusions remain valid for non-exponential models: spatial 3d model and sigmoidal growth, relevant for early- and late stages of cancer growth., Author summary Darwinian evolution in cancer is responsible for the emergence of malignant traits in initially benign tumors. As tumor cells divide, they accumulate new mutations and while most of them are “passengers” which do not confer any selective growth advantage, “driver” mutations endow cells with traits that contribute to cancer spread. Identifying driver mutations that are under selection in cancer can point to new targets for cancer therapeutics and open new avenues for personalized cancer treatment. It has recently been argued that the presence or absence of selection in cancer can be deduced from deviation of mutant allele frequencies from 1/f power law in an intermediate frequency range. Using a stochastic mathematical model of cancer evolution we derive a formula for the frequency of a subclonal driver and show that frequencies of cancer drivers are biased towards 0 and 1; thus most mutations will inevitably appear to be either neutral (frequency ≈ 0) or clonal (frequency ≈ 1) despite very different levels of selection. Consequently, the proposed 1/f statistic will significantly overestimate the number of cancers deemed to be evolving neutrally. Our work quantifies the validity of the proposed neutral evolution statistic across the entire range of relevant parameter values.

Scientific Reports 6: Article number: 39511; published online: 22 December 2016; updated: 20 October 2017. A Supplementary Information data file and Supplementary Mathematica notebook file were omitted from the original version of this Article. This has been corrected in the PDF and HTML version of the Article.

One of the hallmarks of cancer is the accumulation of driver mutations which increase the net reproductive rate of cancer cells and allow them to spread. This process has been studied in mathematical models of well mixed populations, and in computer simulations of three-dimensional spatial models. But the computational complexity of these more realistic, spatial models makes it difficult to simulate realistically large and clinically detectable solid tumours. Here we describe an exactly solvable mathematical model of a tumour featuring replication, mutation and local migration of cancer cells. The model predicts a quasi-exponential growth of large tumours even if different fragments of the tumour grow sub-exponentially due to nutrient and space limitations. The model reproduces clinically observed tumour growth times using biologically plausible rates for cell birth, death, and migration rates. We also show that the expected number of accumulated driver mutations increases exponentially in time if the average fitness gain per driver is constant, and that it reaches a plateau if the gains decrease over time. We discuss the realism of the underlying assumptions and possible extensions of the model.