Source author record

Tibor Antal

Tibor Antal appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2020arXiv

Genetic composition of an exponentially growing cell population

We study a simple model of DNA evolution in a growing population of cells. Each cell contains a nucleotide sequence which randomly mutates at cell division. Cells divide according to a branching process. Following typical parameter values in bacteria and cancer cell populations, we take the mutation rate to zero and the final number of cells to infinity. We prove that almost every site (entry of the nucleotide sequence) is mutated in only a finite number of cells, and these numbers are independent across sites. However independence breaks down for the rare sites which are mutated in a positive fraction of the population. The model is free from the popular but disputed infinite sites assumption. Violations of the infinite sites assumption are widespread while their impact on mutation frequencies is negligible at the scale of population fractions. Some results are generalised to allow for cell death, selection, and site-specific mutation rates. For illustration we estimate mutation rates in a lung adenocarcinoma.

preprint2016arXiv

Universal asymptotic clone size distribution for general population growth

Deterministically growing (wild-type) populations which seed stochastically developing mutant clones have found an expanding number of applications from microbial populations to cancer. The special case of exponential wild-type population growth, usually termed the Luria-Delbrück or Lea-Coulson model, is often assumed but seldom realistic. In this article we generalise this model to different types of wild-type population growth, with mutants evolving as a birth-death branching process. Our focus is on the size distribution of clones - that is the number of progeny of a founder mutant - which can be mapped to the total number of mutants. Exact expressions are derived for exponential, power-law and logistic population growth. Additionally for a large class of population growth we prove that the long time limit of the clone size distribution has a general two-parameter form, whose tail decays as a power-law. Considering metastases in cancer as the mutant clones, upon analysing a data-set of their size distribution, we indeed find that a power-law tail is more likely than an exponential one.

preprint2014arXiv

Mutant number distribution in an exponentially growing population

We present an explicit solution to a classic model of cell-population growth introduced by Luria and Delbrueck 70 years ago to study the emergence of mutations in bacterial populations. In this model a wild-type population is assumed to grow exponentially in a deterministic fashion. Proportional to the wild-type population size, mutants arrive randomly and initiate new sub-populations of mutants that grows stochastically according to a supercritical birth and death process. We give an exact expression for the generating function of the total number of mutants at a given wild type population size. We present a simple expression for the probability of finding no mutants, and a recursion formula for the probability of finding a given number of mutants. In the "large population-small mutation"-limit we recover recent results of Kessler and Levin for a fully stochastic version of the process.

preprint2013arXiv

Spatial evolution of tumors with successive driver mutations

We study the spatial evolutionary dynamics of solid tumors as they obtain additional driver mutations. We start with a cancer clone that expands uniformly in three dimensions giving rise to a spherical shape. We assume that cell division occurs on the surface of the growing tumor. Each cell division has a chance to give rise to a mutation that activates an additional driver gene. The resulting clone has an enhanced growth rate, which generates a local ensemble of faster growing cells, thereby distorting the spherical shape of the tumor. We derive analytic formulas for the geometric boundary that separates the original cancer clone from the new mutant as well as the expanding frontier of the new mutant. The total number of original cancer cells converges to a constant as time goes to infinity, because this clone becomes enveloped by mutants. We derive formulas for the abundance and diversity of additional driver mutations as function of time. Our model is semi-deterministic: the spatial growth of the various cancer clones follows deterministic equations, but the arrival of a new mutant is a stochastic event.

preprint2012arXiv

Molecular Spiders on the Plane

Synthetic bio-molecular spiders with "legs" made of single-stranded segments of DNA can move on a surface covered by single-stranded segments of DNA called substrates when the substrate DNA is complementary to the leg DNA. If the motion of a spider does not affect the substrates, the spider behaves asymptotically as a random walk. We study the diffusion coefficient and the number of visited sites for spiders moving on the square lattice with a substrate in each lattice site. The spider's legs hop to nearest-neighbor sites with the constraint that the distance between any two legs cannot exceed a maximal span. We establish analytic results for bipedal spiders, and investigate multileg spiders numerically. In experimental realizations legs usually convert substrates into products (visited sites). The binding of legs to products is weaker, so the hopping rate from the substrates is smaller. This makes the problem non-Markovian and we investigate it numerically. We demonstrate the emergence of a counter-intuitive behavior - the more spiders are slowed down on unvisited sites, the more motile they become.

preprint2012arXiv

Outbreak size distributions in epidemics with multiple stages

Multiple-type branching processes that model the spread of infectious diseases are investigated. In these stochastic processes, the disease goes through multiple stages before it eventually disappears. We mostly focus on the critical multistage Susceptible-Infected-Recovered (SIR) infection process. In the infinite population limit, we compute the outbreak size distributions and show that asymptotic results apply to more general multiple-type critical branching processes. Finally using heuristic arguments and simulations we establish scaling laws for a multistage SIR model in a finite population.

preprint2011arXiv

Exact solution of a two-type branching process: Models of tumor progression

An explicit solution for a general two-type birth-death branching process with one way mutation is presented. This continuous time process mimics the evolution of resistance to treatment, or the onset of an extra driver mutation during tumor progression. We obtain the exact generating function of the process at arbitrary times, and derive various large time scaling limits. In the simultaneous small mutation rate and large time scaling limit, the distribution of the mutant cells develops some atypical properties, including a power law tail and diverging average.

preprint2010arXiv

First Passage Properties of the Polya Urn Process

We study first passage statistics of the Polya urn model. In this random process, the urn contains two types of balls. In each step, one ball is drawn randomly from the urn, and subsequently placed back into the urn together with an additional ball of the same type. We derive the probability G_n that the two types of balls are equal in number, for the first time, when there is a total of 2n balls. This first passage probability decays algebraically, G_n ~ n^{-2}, when n is large. We also derive the probability that a tie ever happens. This probability is between zero and one, so that a tie may occur in some realizations but not in others. The likelihood of a tie is appreciable only if the initial difference in the number balls is of the order of the square-root of the total number of balls.

preprint2009arXiv

Accumulation of driver and passenger mutations during tumor progression

Major efforts to sequence cancer genomes are now occurring throughout the world. Though the emerging data from these studies are illuminating, their reconciliation with epidemiologic and clinical observations poses a major challenge. In the current study, we provide a novel mathematical model that begins to address this challenge. We model tumors as a discrete time branching process that starts with a single driver mutation and proceeds as each new driver mutation leads to a slightly increased rate of clonal expansion. Using the model, we observe tremendous variation in the rate of tumor development - providing an understanding of the heterogeneity in tumor sizes and development times that have been observed by epidemiologists and clinicians. Furthermore, the model provides a simple formula for the number of driver mutations as a function of the total number of mutations in the tumor. Finally, when applied to recent experimental data, the model allows us to calculate, for the first time, the actual selective advantage provided by typical somatic mutations in human tumors in situ. This selective advantage is surprisingly small, 0.005 +- 0.0005, and has major implications for experimental cancer research.

preprint2009arXiv

Exact solution of a two-type branching process: Clone size distribution in cell division kinetics

We study a two-type branching process which provides excellent description of experimental data on cell dynamics in skin tissue (Clayton et al., 2007). The model involves only a single type of progenitor cell, and does not require support from a self-renewed population of stem cells. The progenitor cells divide and may differentiate into post-mitotic cells. We derive an exact solution of this model in terms of generating functions for the total number of cells, and for the number of cells of different types. We also deduce large time asymptotic behaviors drawing on our exact results, and on an independent diffusion approximation.

preprint2009arXiv

Galaxy distribution and extreme value statistics

We consider the conditional galaxy density around each galaxy, and study its fluctuations in the newest samples of the Sloan Digital Sky Survey Data Release 7. Over a large range of scales, both the average conditional density and its variance show a nontrivial scaling behavior, which resembles to criticality. The density depends, for 10 < r < 80 Mpc/h, only weakly (logarithmically) on the system size. Correspondingly, we find that the density fluctuations follow the Gumbel distribution of extreme value statistics. This distribution is clearly distinguishable from a Gaussian distribution, which would arise for a homogeneous spatial galaxy configuration. We also point out similarities between the galaxy distribution and critical systems of statistical physics.

preprint2007arXiv

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.

preprint2007arXiv

Guiding-fields for phase-separation: Controlling Liesegang patterns

Liesegang patterns emerge from precipitation processes and may be used to build bulk structures at submicron lengthscales. Thus they have significant potential for technological applications provided adequate methods of control can be devised. Here we describe a simple, physically realizable pattern-control based on the notion of driven precipitation, meaning that the phase-separation is governed by a guiding field such as, for example, a temperature or a pH field. The phase-separation is modeled through a non-autonomous Cahn-Hilliard equation whose spinodal is determined by the evolving guiding field. Control over the dynamics of the spinodal gives control over the velocity of the instability front which separates the stable and unstable regions of the system. Since the wavelength of the pattern is largely determined by this velocity, the distance between successive precipitation bands becomes controllable. We demonstrate the above ideas by numerical studies of a 1D system with diffusive guiding field. We find that the results can be accurately described by employing a linear stability analysis (pulled-front theory) for determining the velocity -- local-wavelength relationship. From the perspective of the Liesegang theory, our results indicate that the so-called revert patterns may be naturally generated by diffusive guiding fields.

preprint2002arXiv

On the Role of External Constraints in a Spatially Extended Evolutionary Prisoner's Dilemma Game

We study the emergency of mutual cooperation in evolutionary prisoner's dilemma games when the players are located on a square lattice. The players can choose one of the three strategies: cooperation (C), defection (D) or "tit for tat" (T), and their total payoffs come from games with the nearest neighbors. During the random sequential updates the players adopt one of their neighboring strategies if the chosen neighbor has higher payoff. We compare the effect of two types of external constraints added to the Darwinian evolutionary processes. In both cases the strategy of a randomly chosen player is replaced with probability P by another strategy. In the first case, the strategy is replaced by a randomly chosen one among the two others, while in the second case the new strategy is always C. Using generalized mean-field approximations and Monte Carlo simulations the strategy concentrations are evaluated in the stationary state for different strength of external constraints characterized by the probability P.