Genetic progression and the waiting time to cancer

Cancer results from genetic alterations that disturb the normal cooperative behavior of cells. Recent high-throughput genomic studies of cancer cells have shown that the mutational landscape of cancer is complex and that individual cancers may evolve through mutations in as many as 20 different cancer-associated genes. We use data published by Sjoblom et al. (2006) to develop a new mathematical model for the somatic evolution of colorectal cancers. We employ the Wright-Fisher process for exploring the basic parameters of this evolutionary process and derive an analytical approximation for the expected waiting time to the cancer phenotype. Our results highlight the relative importance of selection over both the size of the cell population at risk and the mutation rate. The model predicts that the observed genetic diversity of cancer genomes can arise under a normal mutation rate if the average selective advantage per mutation is on the order of 1%. Increased mutation rates due to genetic instability would allow even smaller selective advantages during tumorigenesis. The complexity of cancer progression thus can be understood as the result of multiple sequential mutations, each of which has a relatively small but positive effect on net cell growth.

💡 Research Summary

The paper tackles the fundamental question of how many genetic alterations and how much time are required for a normal cell to become a malignant cancer cell. Using the extensive mutational data from Sjöblom et al. (2006), which identified mutations in roughly twenty cancer‑associated genes in colorectal tumors, the authors construct a stochastic evolutionary model based on the classic Wright–Fisher process. In this framework, a fixed population of N cells (representing the epithelial stem cell pool of the colon, on the order of 10⁸–10⁹ cells) reproduces in discrete generations. Each replication event carries a per‑gene mutation probability μ, and each newly acquired driver mutation confers a selective advantage s that increases the probability that the mutant lineage will be chosen as a parent in the next generation.

The model assumes that driver mutations act independently, each providing the same average fitness increment, and that cancer is defined as the first appearance of a cell that has accumulated a threshold number k of such driver mutations (typically k ≈ 5–7, consistent with empirical observations). By running extensive Monte‑Carlo simulations across a wide range of (N, μ, s) parameter sets, the authors measure the expected waiting time T(k) until a cell with k drivers emerges. They then derive an analytical approximation by treating the sequential acquisition of mutations as a series of rare events in a Markov chain. The resulting expression,

\