Estimating the transfer rates of bacterial plasmids with an adapted Luria-Delbrück fluctuation analysis.

To increase our basic understanding of the ecology and evolution of conjugative plasmids, we need reliable estimates of their rate of transfer between bacterial cells. Current assays to measure transfer rate are based on deterministic modeling frameworks. However, some cell numbers in these assays can be very small, making estimates that rely on these numbers prone to noise. Here, we take a different approach to estimate plasmid transfer rate, which explicitly embraces this noise. Inspired by the classic fluctuation analysis of Luria and Delbrück, our method is grounded in a stochastic modeling framework. In addition to capturing the random nature of plasmid conjugation, our new methodology, the Luria-Delbrück method ("LDM"), can be used on a diverse set of bacterial systems, including cases for which current approaches are inaccurate. A notable example involves plasmid transfer between different strains or species where the rate that one type of cell donates the plasmid is not equal to the rate at which the other cell type donates. Asymmetry in these rates has the potential to bias or constrain current transfer estimates, thereby limiting our capabilities for estimating transfer in microbial communities. In contrast, the LDM overcomes obstacles of traditional methods by avoiding restrictive assumptions about growth and transfer rates for each population within the assay. Using stochastic simulations and experiments, we show that the LDM has high accuracy and precision for estimation of transfer rates compared to the most widely used methods, which can produce estimates that differ from the LDM estimate by orders of magnitude.

Growth dynamics in naturally progressing chronic lymphocytic leukaemia.

How the genomic features of a patient's cancer relate to individual disease kinetics remains poorly understood. Here we used the indolent growth dynamics of chronic lymphocytic leukaemia (CLL) to analyse the growth rates and corresponding genomic patterns of leukaemia cells from 107 patients with CLL, spanning decades-long disease courses. We found that CLL commonly demonstrates not only exponential expansion but also logistic growth, which is sigmoidal and reaches a certain steady-state level. Each growth pattern was associated with marked differences in genetic composition, the pace of disease progression and the extent of clonal evolution. In a subset of patients, whose serial samples underwent next-generation sequencing, we found that dynamic changes in the disease course of CLL were shaped by the genetic events that were already present in the early slow-growing stages. Finally, by analysing the growth rates of subclones compared with their parental clones, we quantified the growth advantage conferred by putative CLL drivers in vivo.

Clonal interactions in cancer: Integrating quantitative models with experimental and clinical data.

Tumors consist of different genotypically distinct subpopulations-or subclones-of cells. These subclones can influence neighboring clones in a process called "clonal interaction." Conventionally, research on driver mutations in cancer has focused on their cell-autonomous effects that lead to an increase in fitness of the cells containing the driver. Recently, with the advent of improved experimental and computational technologies for investigating tumor heterogeneity and clonal dynamics, new studies have shown the importance of clonal interactions in cancer initiation, progression, and metastasis. In this review we provide an overview of clonal interactions in cancer, discussing key discoveries from a diverse range of approaches to cancer biology research. We discuss common types of clonal interactions, such as cooperation and competition, its mechanisms, and the overall effect on tumorigenesis, with important implications for tumor heterogeneity, resistance to treatment, and tumor suppression. Quantitative models-in coordination with cell culture and animal model experiments-have played a vital role in investigating the nature of clonal interactions and the complex clonal dynamics they generate. We present mathematical and computational models that can be used to represent clonal interactions and provide examples of the roles they have played in identifying and quantifying the strength of clonal interactions in experimental systems. Clonal interactions have proved difficult to observe in clinical data; however, several very recent quantitative approaches enable their detection. We conclude by discussing ways in which researchers can further integrate quantitative methods with experimental and clinical data to elucidate the critical-and often surprising-roles of clonal interactions in human cancers.

Ultra-deep mutational landscape in chronic lymphocytic leukemia uncovers dynamics of resistance to targeted therapies.

BTK inhibitors, Bcl-2 inhibitors, and other targeted therapies have significantly improved the outcomes of patients with chronic lymphocytic leukemia (CLL). With increased survivorship, monitoring disease and deciphering potential mechanisms of resistance to these agents are critical for devising effective treatment strategies. We used duplex sequencing, a technology that enables detection of mutations at ultra-low allelic frequencies, to identify mutations in five genes associated with drug resistance in CLL and followed their evolution in two patients who received multiple targeted therapies and ultimately developed disease progression on pirtobrutinib. In both patients we detected variants that expanded and reached significant cancer cell fractions (CCF). In patient R001, multiple known resistance mutations in both BTK and PLCG2 appeared following progression on zanubrutinib (BTK p.L528W, p.C481S; PLCG2 S707F, L845F, R665W, and D993H). In contrast, patient R002 developed multiple BTK mutations following acalabrutinib treatment, including known resistance mutations p.C481R, p.T474I and p.C481S. We found that pirtobrutinib was able to suppress, but not completely eradicate, BTK p.C481S mutations in both patients, but other resistance mutations such as mutations in PLCG2 and new BTK mutations increased while the patients were receiving pirtobrutinib. For example, BTK p.L528W in patient R001 increased in frequency more than 1,000-fold (from a CCF of 0.02% to 35%), and the CCF in p.T474I in patient R002 increased from 0.03% to 4.2% (more than 100-fold). Our data illuminate the evolutionary dynamics of resistant clones over the patients' disease course and under selective pressure from different targeted treatments.

Stochastic model for cell population dynamics quantifies homeostasis in colonic crypts and its disruption in early tumorigenesis.

The questions of how healthy colonic crypts maintain their size, and how homeostasis is disrupted by driver mutations, are central to understanding colorectal tumorigenesis. We propose a three-type stochastic branching process, which accounts for stem, transit-amplifying (TA) and fully differentiated (FD) cells, to model the dynamics of cell populations residing in colonic crypts. Our model is simple in its formulation, allowing us to estimate all but one of the model parameters from the literature. Fitting the single remaining parameter, we find that model results agree well with data from healthy human colonic crypts, capturing the considerable variance in population sizes observed experimentally. Importantly, our model predicts a steady-state population in healthy colonic crypts for relevant parameter values. We show that APC and KRAS mutations, the most significant early alterations leading to colorectal cancer, result in increased steady-state populations in mutated crypts, in agreement with experimental results. Finally, our model predicts a simple condition for unbounded growth of cells in a crypt, corresponding to colorectal malignancy. This is predicted to occur when the division rate of TA cells exceeds their differentiation rate, with implications for therapeutic cancer prevention strategies.

Waiting times in a branching process model of colorectal cancer initiation.

We study a multi-stage model for the development of colorectal cancer from initially healthy tissue. The model incorporates a complex sequence of driver gene alterations, some of which result in immediate growth advantage, while others have initially neutral effects. We derive analytic estimates for the sizes of premalignant subpopulations, and use these results to compute the waiting times to premalignant and malignant genotypes. This work contributes to the quantitative understanding of colorectal tumor evolution and the lifetime risk of colorectal cancer.

Individual and environmental correlates of childhood maltreatment and exposure to community violence: Utilizing a latent profile and a multilevel meta-analytic approach.

BACKGROUND: Childhood maltreatment (CM) and exposure to community violence (ETV) are correlated with physical/mental health and psychosocial problems. Typically, CM and ETV are examined separately, by subtypes within category, or collapsed across both into one category of adversity. Consequently, research is limited in identifying subgroups of individuals with different amounts of exposure to both CM and ETV. Accordingly, we lack sufficient understanding of the extent to which problems associated with CM and ETV vary based on the amount (i.e. dose) of exposure to both of these experiences. METHODS: We used 20 samples (28,300 individuals) to estimate person-centered profiles of CM and ETV occurrence and co-occurrence within each sample. An individual data multilevel meta-analytic framework was used to determine the average effect size across samples for different profiles and conditional probability correlations within sociodemographic, neighborhood, health, mental health, and psychosocial domains. RESULTS: The profile characterized by high levels of CM and high levels of ETV correlated with stressful life events, depression and anxiety symptoms, and general indicators of externalizing behaviors. CM predominant profiles were associated with mental health diagnoses and treatment. ETV predominant profiles associated with risk-taking/violent behavior and neighborhood-level disadvantage. However, nuance based on the dose of CM or ETV was evident. CONCLUSIONS: It is important to identify subgroups based on the amount of exposure to CM and ETV. These subgroups have differential relationships with correlates across domains. Greater delineation and description of the lived experience will allow for more precision in addressing the burden of childhood adversity.

Inferring parameters of cancer evolution in chronic lymphocytic leukemia.

As a cancer develops, its cells accrue new mutations, resulting in a heterogeneous, complex genomic profile. We make use of this heterogeneity to derive simple, analytic estimates of parameters driving carcinogenesis and reconstruct the timeline of selective events following initiation of an individual cancer, where two longitudinal samples are available for sequencing. Using stochastic computer simulations of cancer growth, we show that we can accurately estimate mutation rate, time before and after a driver event occurred, and growth rates of both initiated cancer cells and subsequently appearing subclones. We demonstrate that in order to obtain accurate estimates of mutation rate and timing of events, observed mutation counts should be corrected to account for clonal mutations that occurred after the founding of the tumor, as well as sequencing coverage. Chronic lymphocytic leukemia (CLL), which often does not require treatment for years after diagnosis, presents an optimal system to study the untreated, natural evolution of cancer cell populations. When we apply our methodology to reconstruct the individual evolutionary histories of CLL patients, we find that the parental leukemic clone typically appears within the first fifteen years of life.

Mathematical model of colorectal cancer initiation.

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 premalignant mutational genotypes on the way to colorectal cancer. We parameterize 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 lifetime risk 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 primarily by the fitness effects that they provide, rather than their mutation rates. Our results imply that there may not be significant immune suppression of untreated benign and malignant colorectal lesions.

A mechanistic mathematical model of initiation and malignant transformation in sporadic vestibular schwannoma.

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.

A spatial model predicts that dispersal and cell turnover limit intratumour heterogeneity.

Most cancers in humans are large, measuring centimetres in diameter, and composed of many billions of cells. An equivalent mass of normal cells would be highly heterogeneous as a result of the mutations that occur during each cell division. What is remarkable about cancers is that virtually every neoplastic cell within a large tumour often contains the same core set of genetic alterations, with heterogeneity confined to mutations that emerge late during tumour growth. How such alterations expand within the spatially constrained three-dimensional architecture of a tumour, and come to dominate a large, pre-existing lesion, has been unclear. Here we describe a model for tumour evolution that shows how short-range dispersal and cell turnover can account for rapid cell mixing inside the tumour. We show that even a small selective advantage of a single cell within a large tumour allows the descendants of that cell to replace the precursor mass in a clinically relevant time frame. We also demonstrate that the same mechanisms can be responsible for the rapid onset of resistance to chemotherapy. Our model not only provides insights into spatial and temporal aspects of tumour growth, but also suggests that targeting short-range cellular migratory activity could have marked effects on tumour growth rates.

Mutations driving CLL and their evolution in progression and relapse.

Which genetic alterations drive tumorigenesis and how they evolve over the course of disease and therapy are central questions in cancer biology. Here we identify 44 recurrently mutated genes and 11 recurrent somatic copy number variations through whole-exome sequencing of 538 chronic lymphocytic leukaemia (CLL) and matched germline DNA samples, 278 of which were collected in a prospective clinical trial. These include previously unrecognized putative cancer drivers (RPS15, IKZF3), and collectively identify RNA processing and export, MYC activity, and MAPK signalling as central pathways involved in CLL. Clonality analysis of this large data set further enabled reconstruction of temporal relationships between driver events. Direct comparison between matched pre-treatment and relapse samples from 59 patients demonstrated highly frequent clonal evolution. Thus, large sequencing data sets of clinically informative samples enable the discovery of novel genes associated with cancer, the network of relationships between the driver events, and their impact on disease relapse and clinical outcome.

Modeling CAR T-Cell Therapy with Patient Preconditioning.

The Federal Drug Administration approved the first Chimeric Antigen Receptor T-cell (CAR T-cell) therapies for the treatment of several blood cancers in 2017, and efforts are underway to broaden CAR T technology to address other cancer types. Standard treatment protocols incorporate a preconditioning regimen of lymphodepleting chemotherapy prior to CAR T-cell infusion. However, the connection between preconditioning regimens and patient outcomes is still not fully understood. Optimizing patient preconditioning plans and reducing the CAR T-cell dose necessary for achieving remission could make therapy safer. In this paper, we test treatment regimens consisting of sequential administration of chemotherapy and CAR T-cell therapy on a system of differential equations that models the tumor-immune interaction. We use numerical simulations of treatment plans from within the scope of current medical practice to assess the effect of preconditioning plans on the success of CAR T-cell therapy. Model results affirm clinical observations that preconditioning can be crucial for most patients, not just to reduce side effects, but to even achieve remission at all. We demonstrate that preconditioning plans using the same CAR T-cell dose and the same total concentration of chemotherapy can lead to different patient outcomes due to different delivery schedules. Results from sensitivity analysis of the model parameters suggest that making small improvements in the effectiveness of CAR T-cells in attacking cancer cells will significantly reduce the minimum dose required for successful treatment. Our modeling framework represents a starting point for evaluating the efficacy of patient preconditioning in the context of CAR T-cell therapy.

On measuring selection in cancer from subclonal mutation frequencies.

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.

The molecular evolution of acquired resistance to targeted EGFR blockade in colorectal cancers.

Colorectal tumours that are wild type for KRAS are often sensitive to EGFR blockade, but almost always develop resistance within several months of initiating therapy. The mechanisms underlying this acquired resistance to anti-EGFR antibodies are largely unknown. This situation is in marked contrast to that of small-molecule targeted agents, such as inhibitors of ABL, EGFR, BRAF and MEK, in which mutations in the genes encoding the protein targets render the tumours resistant to the effects of the drugs. The simplest hypothesis to account for the development of resistance to EGFR blockade is that rare cells with KRAS mutations pre-exist at low levels in tumours with ostensibly wild-type KRAS genes. Although this hypothesis would seem readily testable, there is no evidence in pre-clinical models to support it, nor is there data from patients. To test this hypothesis, we determined whether mutant KRAS DNA could be detected in the circulation of 28 patients receiving monotherapy with panitumumab, a therapeutic anti-EGFR antibody. We found that 9 out of 24 (38%) patients whose tumours were initially KRAS wild type developed detectable mutations in KRAS in their sera, three of which developed multiple different KRAS mutations. The appearance of these mutations was very consistent, generally occurring between 5 and 6 months following treatment. Mathematical modelling indicated that the mutations were present in expanded subclones before the initiation of panitumumab treatment. These results suggest that the emergence of KRAS mutations is a mediator of acquired resistance to EGFR blockade and that these mutations can be detected in a non-invasive manner. They explain why solid tumours develop resistance to targeted therapies in a highly reproducible fashion.

Quantifying Clonal and Subclonal Passenger Mutations in Cancer Evolution.

The vast majority of mutations in the exome of cancer cells are passengers, which do not affect the reproductive rate of the cell. Passengers can provide important information about the evolutionary history of an individual cancer, and serve as a molecular clock. Passengers can also become targets for immunotherapy or confer resistance to treatment. We study the stochastic expansion of a population of cancer cells describing the growth of primary tumors or metastatic lesions. We first analyze the process by looking forward in time and calculate the fixation probabilities and frequencies of successive passenger mutations ordered by their time of appearance. We compute the likelihood of specific evolutionary trees, thereby informing the phylogenetic reconstruction of cancer evolution in individual patients. Next, we derive results looking backward in time: for a given subclonal mutation we estimate the number of cancer cells that were present at the time when that mutation arose. We derive exact formulas for the expected numbers of subclonal mutations of any frequency. Fitting this formula to cancer sequencing data leads to an estimate for the ratio of birth and death rates of cancer cells during the early stages of clonal expansion.

Timing and heterogeneity of mutations associated with drug resistance in metastatic cancers.

Targeted therapies provide an exciting new approach to combat human cancer. The immediate effect is a dramatic reduction in disease burden, but in most cases, the tumor returns as a consequence of resistance. Various mechanisms for the evolution of resistance have been implicated, including mutation of target genes and activation of other drivers. There is increasing evidence that the reason for failure of many targeted treatments is a small preexisting subpopulation of resistant cells; however, little is known about the genetic composition of this resistant subpopulation. Using the novel approach of ordering the resistant subclones according to their time of appearance, here we describe the full spectrum of resistance mutations present in a metastatic lesion. We calculate the expected and median number of cells in each resistant subclone. Surprisingly, the ratio of the medians of successive resistant clones is independent of any parameter in our model; for example, the median of the second clone divided by the median of the first is √2-1. We find that most radiographically detectable lesions harbor at least 10 resistant subclones. Our predictions are in agreement with clinical data on the relative sizes of resistant subclones obtained from liquid biopsies of colorectal cancer patients treated with epidermal growth factor receptor (EGFR) blockade. Our theory quantifies the genetic heterogeneity of resistance that exists before treatment and provides information to design treatment strategies that aim to control resistance.

Tumour and immune cell dynamics explain the PSA bounce after prostate cancer brachytherapy.

BACKGROUND: Interstitial brachytherapy for localised prostate cancer may be followed by transient increases in prostate-specific antigen (PSA) that resolve without therapy. Such PSA bounces may be associated with an improved outcome but often cause alarm in the patient and physician, and have defied explanation. METHODS: We developed a mathematical model to capture the interactions between the tumour, radiation and anti-tumour immune response. The model was fitted to data from a large cohort of patients treated exclusively with interstitial brachytherapy. Immunohistological analysis for T-cell infiltration within the same tumours was also performed. RESULTS: Our minimal model captures well the dynamics of the tumour after therapy, and suggests that a strong anti-tumour immune response coupled with the therapeutic effect of radiation on the tumour is responsible for the PSA bounce. Patients who experience a PSA bounce had a higher density of CD3 and CD8 cells within the tumour that likely contribute to the PSA bounce and the overall better outcomes observed. CONCLUSIONS: Our observations provide a novel and unifying explanation for the PSA bounce in patients with early prostate cancer and also have implications for the use of immune-based therapies in such patients to improve outcomes.

Distant metastasis occurs late during the genetic evolution of pancreatic cancer.

Metastasis, the dissemination and growth of neoplastic cells in an organ distinct from that in which they originated, is the most common cause of death in cancer patients. This is particularly true for pancreatic cancers, where most patients are diagnosed with metastatic disease and few show a sustained response to chemotherapy or radiation therapy. Whether the dismal prognosis of patients with pancreatic cancer compared to patients with other types of cancer is a result of late diagnosis or early dissemination of disease to distant organs is not known. Here we rely on data generated by sequencing the genomes of seven pancreatic cancer metastases to evaluate the clonal relationships among primary and metastatic cancers. We find that clonal populations that give rise to distant metastases are represented within the primary carcinoma, but these clones are genetically evolved from the original parental, non-metastatic clone. Thus, genetic heterogeneity of metastases reflects that within the primary carcinoma. A quantitative analysis of the timing of the genetic evolution of pancreatic cancer was performed, indicating at least a decade between the occurrence of the initiating mutation and the birth of the parental, non-metastatic founder cell. At least five more years are required for the acquisition of metastatic ability and patients die an average of two years thereafter. These data provide novel insights into the genetic features underlying pancreatic cancer progression and define a broad time window of opportunity for early detection to prevent deaths from metastatic disease.