Showing papers in "arXiv: Populations and Evolution in 2005"

A nonlocal continuum model for biological aggregation

Abstract: We construct a continuum model for biological aggregations in which individuals experience long-range social attraction and short range dispersal. For the case of one spatial dimension, we study the steady states analytically and numerically. There exist strongly nonlinear states with compact support and steep edges that correspond to localized biological aggregations, or clumps. These steady state clumps are approached through a dynamic coarsening process. In the limit of large population size, the clumps approach a constant density swarm with abrupt edges. We use energy arguments to understand the nonlinear selection of clump solutions, and to predict the internal density in the large population limit. The energy result holds in higher dimensions as well, and is demonstrated via numerical simulations in two dimensions.

Fixation of strategies for an evolutionary game in finite populations

Abstract: A stochastic evolutionary dynamics of two strategies given by 2 x 2 matrix games is studied in finite populations. We focus on stochastic properties of fixation: how a strategy represented by a single individual wins over the entire population. The process is discussed in the framework of a random walk with arbitrary hopping rates. The time of fixation is found to be identical for both strategies in any particular game. The asymptotic behavior of the fixation time and fixation probabilities in the large population size limit is also discussed. We show that fixation is fast when there is at least one pure evolutionary stable strategy (ESS) in the infinite population size limit, while fixation is slow when the ESS is the coexistence of the two strategies.

Mimivirus and the emerging concept of "giant" virus

Abstract: The recently discovered Acanthamoeba polyphaga Mimivirus is the largest known DNA virus. Its particle size (>400 nm), genome length (1.2 million bp) and large gene repertoire (911 protein coding genes) blur the established boundaries between viruses and parasitic cellular organisms. In addition, the analysis of its genome sequence identified new types of genes not expected to be seen in a virus, such as aminoacyl-tRNA synthetases and other central components of the translation machinery. In this article, we examine how the finding of a giant virus for the first time overlapping with the world of cellular organisms in terms of size and genome complexity might durably influence the way we look at microbial biodiversity, and force us to fundamentally revise our classification of life forms. We propose to introduce the word "girus" to recognize the intermediate status of these giant DNA viruses, the genome complexity of which make them closer to small parasitic prokaryotes than to regular viruses.

Food Webs: Experts Consuming Families of Experts

Abstract: The question what determines the structure of natural food webs has been listed among the nine most important unanswered questions in ecology. It arises naturally from many problems related to ecosystem stability and resilience. The traditional view is that population-dynamical stability is crucial for understanding the observed structures. But phylogeny (evolutionary history) has also been suggested as the dominant mechanism. Here we show that observed topological features of predatory food webs can be reproduced to unprecedented accuracy by a mechanism taking into account only phylogeny, size constraints, and the heredity of the trophically relevant traits of prey and predators. The analysis reveals a tendency to avoid resource competition rather than apparent competition. In food webs with many parasites this pattern is reversed.

Phase Transitions and Spatio-Temporal Fluctuations in Stochastic Lattice Lotka-Volterra Models

Abstract: We study the general properties of stochastic two-species models for predator-prey competition and coexistence with Lotka-Volterra type interactions defined on a $d$-dimensional lattice. Introducing spatial degrees of freedom and allowing for stochastic fluctuations generically invalidates the classical, deterministic mean-field picture. Already within mean-field theory, however, spatial constraints, modeling locally limited resources, lead to the emergence of a continuous active-to-absorbing state phase transition. Field-theoretic arguments, supported by Monte Carlo simulation results, indicate that this transition, which represents an extinction threshold for the predator population, is governed by the directed percolation universality class. In the active state, where predators and prey coexist, the classical center singularities with associated population cycles are replaced by either nodes or foci. In the vicinity of the stable nodes, the system is characterized by essentially stationary localized clusters of predators in a sea of prey. Near the stable foci, however, the stochastic lattice Lotka-Volterra system displays complex, correlated spatio-temporal patterns of competing activity fronts. Correspondingly, the population densities in our numerical simulations turn out to oscillate irregularly in time, with amplitudes that tend to zero in the thermodynamic limit. Yet in finite systems these oscillatory fluctuations are quite persistent, and their features are determined by the intrinsic interaction rates rather than the initial conditions. We emphasize the robustness of this scenario with respect to various model perturbations.

Biodiversity in model ecosystems, I: Coexistence conditions for competing species

Abstract: This is the first of two papers where we discuss the limits imposed by competition to the biodiversity of species communities. In this first paper we study the coexistence of competing species at the fixed point of population dynamic equations. For many simple models, this imposes a limit on the width of the productivity distribution, which is more severe the more diverse the ecosystem is (Chesson, 1994). Here we review and generalize this analysis, beyond the ``mean-field''-like approximation of the competition matrix used in previous works, and extend it to structured food webs. In all cases analysed, we obtain qualitatively similar relations between biodiversity and competition: the narrower the productivity distribution is, the more species can stably coexist. We discuss how this result, considered together with environmental fluctuations, limits the maximal biodiversity that a trophic level can host.

The identifiability of tree topology for phylogenetic models, including covarion and mixture models

Abstract: For a model of molecular evolution to be useful for phylogenetic inference, the topology of evolutionary trees must be identifiable. That is, from a joint distribution the model predicts, it must be possible to recover the tree parameter. We establish tree identifiability for a number of phylogenetic models, including a covarion model and a variety of mixture models with a limited number of classes. The proof is based on the introduction of a more general model, allowing more states at internal nodes of the tree than at leaves, and the study of the algebraic variety formed by the joint distributions to which it gives rise. Tree identifiability is first established for this general model through the use of certain phylogenetic invariants.

Adaptation in simple and complex fitness landscapes

Abstract: This is an introductory review of deterministic mutation-selection models for asexual populations (i.e., quasispecies theory) and related topics. First, the basic concepts of fitness, mutations, and sequence space are introduced. Different types of mutation-selection dynamics are defined and their relation to problems of statistical physics are outlined. Then the stationary population distribution in simple, single peak fitness landscapes is discussed at length, with particular emphasis on the error threshold phenomenon. Extensions of the theory covering e.g. epistatic interactions, diploid organisms, semiconservative replication and time-dependent fitness peaks are briefly described. A further section is devoted to randomly rugged fitness landscapes, which may display fitness correlations of various degree as well as extended neutral networks. The final two sections address evolutionary dynamics in both rugged and smooth fitness landscapes, and provide a brief overview of pertinent experiments.

Epidemics in small world networks

Abstract: For many infectious diseases, a small-world network on an underlying regular lattice is a suitable simplified model for the contact structure of the host population. It is well known that the contact network, described in this setting by a single parameter, the small-world parameter $p$, plays an important role both in the short term and in the long term dynamics of epidemic spread. We have studied the effect of the network structure on models of immune for life diseases and found that in addition to the reduction of the effective transmission rate, through the screening of infectives, spatial correlations may strongly enhance the stochastic fluctuations. As a consequence, time series of unforced Susceptible-Exposed-Infected-Recovered (SEIR) models provide patterns of recurrent epidemics with realistic amplitudes, suggesting that these models together with complex networks of contacts are the key ingredients to describe the prevaccination dynamical patterns of diseases such as measles and pertussis. We have also studied the role of the host contact strucuture in pathogen antigenic variation, through its effect on the final outcome of an invasion by a viral strain of a population where a very similar virus is endemic. Similar viral strains are modelled by the same infection and reinfection parameters, and by a given degree of cross immunity that represents the antigenic distance between the competing strains. We have found, somewhat surprisingly, that clustering on the network decreases the potential to sustain pathogen diversity.

Maximum Likelihood Jukes-Cantor Triplets: Analytic Solutions

Abstract: Complex systems of polynomial equations have to be set up and solved algebraically in order to obtain analytic solutions for maximum likelihood on phylogenetic trees. This has restricted the types of systems previously resolved to the simplest models - three and four taxa under a molecular clock, with just two state characters. In this work we give, for the first time, analytic solutions for a family of trees with four state characters, like normal DNA or RNA. The model of substitution we use is the Jukes-Cantor model, and the trees are on three taxa under molecular clock, namely rooted triplets. We employ a number of approaches and tools to solve this system: Spectral methods (Hadamard conjugation), a new representation of variables (the path-set spectrum), and algebraic geometry tools (the resultant of two polynomials). All these, combined with heavy application of computer algebra packages (Maple), let us derive the desired solution.

Diffusion and Home Range Parameters from Rodent Population Measurements in Panama

Abstract: Simple random walk considerations are used to interpret rodent population data collected in Hantavirus-related investigations in Panama regarding the short-tailed cane mouse, \emph{Zygodontomys brevicauda}. The diffusion constant of mice is evaluated to be of the order of (and larger than) 200 meters squared per day. The investigation also shows that the rodent mean square displacement saturates in time, indicating the existence of a spatial scale which could, in principle, be the home range of the rodents. This home range is concluded to be of the order of 70 meters. Theoretical analysis is provided for interpreting animal movement data in terms of an interplay of the home ranges, the diffusion constant, and the size of the grid used to monitor the movement. The study gives impetus to a substantial modification of existing theory of the spread of the Hantavirus epidemic which has been based on simple diffusive motion of the rodents, and additionally emphasizes the importance for developing more accurate techniques for the measurement of rodent movement.

Biodiversity in model ecosystems, II: Species assembly and food web structure

Abstract: This is the second of two papers dedicated to the relationship between population models of competition and biodiversity. Here we consider species assembly models where the population dynamics is kept far from fixed points through the continuous introduction of new species, and generalize to such models thecoexistence condition derived for systems at the fixed point. The ecological overlap between species with shared preys, that we define here, provides a quantitative measure of the effective interspecies competition and of the trophic network topology. We obtain distributions of the overlap from simulations of a new model based both on immigration and speciation, and show that they are in good agreement with those measured for three large natural food webs. As discussed in the first paper, rapid environmental fluctuations, interacting with the condition for coexistence of competing species, limit the maximal biodiversity that a trophic level can host. This horizontal limitation to biodiversity is here combined with either dissipation of energy or growth of fluctuations, which in our model limit the length of food webs in the vertical direction. These ingredients yield an effective model of food webs that produce a biodiversity profile with a maximum at an intermediate trophic level, in agreement with field studies.

Evolution as context-driven actualization of potential: Toward an interdisciplinary theory of change of state

Abstract: It is increasingly evident that there is more to biological evolution than natural selection; moreover, the concept of evolution is not limited to biology. We propose an integrative framework for characterizing how entities evolve, in which evolution is viewed as a process of context-driven actualization of potential (CAP). Processes of change differ according to the degree of nondeterminism, and the degree to which they are sensitive to, internalize, and depend upon a particular context. The approach enables us to embed phenomena across disciplines into a broad conceptual framework. We give examples of insights into physics, biology, culture and cognition that derive from this unifying framework.

The Tangled nature model with inheritance and constraint: Evolutionary ecology restricted by a conserved resource

Abstract: We study a version of the Tangled Nature model of evolutionary ecology redefined in a phenotype space where mutants have properties correlated to their parents. The model has individual-based dynamics whilst incorporating species scale competitive constraints and a system scale resource constraint. Multiple species arise that coexist in a species interaction network with evolving global properties. Both the mean interaction strength and the network connectance increase relative to the null system as mutualism becomes more extensive. From a study of the dependence of average degree on the resource level we extract the diversity-connectance relationship which conforms to the hyperbolic form seen in field data. This is adjudged to arise as a consequence of the evolutionary pressure to achieve positive interactions. The network degree distributions conform more strongly to exponential than to the null binomial distributions in all cases. This effect is believed to be caused by correlations in the reproductive process. We also study how resource availability influences the phenotypical lifetime distribution which is approximately of power law form. We observe that the mean lifetime is inversely related to the resource level.

A minimal stochastic model for influenza evolution

Abstract: We introduce and discuss a minimal individual-based model for influenza dynamics. The model takes into account the effects of specific immunization against viral strains, but also infectivity randomness and the presence of a short-lived strain transcending immunity recently suggested in the literature. We show by simulations that the resulting model exhibits substitution of viral strains along the years, but that their divergence remains bounded. We also show that dropping any of these features results in a drastically different behavior, leading either to the extinction of the disease, to the proliferation of the viral strains, or to their divergence.

The response of amino acid frequencies to directional mutation pressure in mitochondrial genome sequences is related to the physical properties of the amino acids and to the structure of the genetic code

Abstract: The frequencies of A, C, G and T in mitochondrial DNA vary among species due to unequal rates of mutation between the bases. The frequencies of bases at four-fold degenerate sites respond directly to mutation pressure. At 1st and 2nd positions, selection reduces the degree of frequency variation. Using a simple evolutionary model, we show that 1st position sites are less constrained by selection than 2nd position sites, and therefore that the frequencies of bases at 1st position are more responsive to mutation pressure than those at 2nd position. We define a similarity measure between amino acids that is a function of 8 measured physical properties. We define a proximity measure for each amino acid, which is the average similarity between an amino acid and all others that are accessible via single point mutations in the genetic code. We also define a responsiveness for each amino acid, which measures how rapidly an amino acid frequency changes as a result of mutation pressure acting on the base frequencies. There is a strong correlation between responsiveness and proximity, and both these quantities are also correlated with the mutability of amino acids estimated from the mtREV substitution rate matrix. We also consider the variation of base frequencies between strands and between genes on a strand. These trends are consistent with the patterns expected from analysis of the variation among genomes

Evolution of Variance in Offspring Number: the Effects of Population Size and Migration

Abstract: It was shown by Gillespie (1974) that if two genotypes produce the same average number of offspring on but have a different variance associated within each generation, the genotype with a lower variance will have a higher effective fitness. Specifically, the effective fitness is W(effective)=w-var/N, where w is the mean fitness, var is the variance in offspring number, and N is the total population size. The model also predicts that if a strategy has a higher arithmetic mean fitness and a higher variance than the competitor, the outcome of selection will depend on the population size (with larger population sizes favoring the high variance, high mean genotype). This suggests that for metapopulations with large numbers of (relatively) small demes, a strategy with lower variance and lower mean may be favored if the migration rate is low while higher migration rates (consistent with a larger effective population size) favor the opposite strategy. Individual based simulation confirms that this is indeed the case for an island model of migration, though the effect of migration differs greatly depending on whether migration precedes or follows selection. It is noted in the appendix that while Gillespie 1974 does seem to be heuristically accurate, it is not clear that the definition of effective fitness follows from his derivation.

Pair Approximation Models for Disease Spread

Abstract: We consider a Susceptible-Infective-Recovered (SIR) model, where the mechanism for the renewal of susceptibles is demographic, on a ring with next nearest neighbour interactions, and a family of correlated pair approximations (CPA), parametrized by a measure of the relative contributions of loops and open triplets of the sites involved in the infection process. We have found that the phase diagram of the CPA, at fixed coordination number, changes qualitatively as the relative weight of the loops increases, from the phase diagram of the uncorrelated pair approximation to phase diagrams typical of one-dimensional systems. In addition, we have performed computer simulations of the same model and shown that while the CPA with a constant correlation parameter cannot describe the global behaviour of the model, a reasonable description of the endemic equilibria as well as of the phase diagram may be obtained by allowing the parameter to depend on the demographic rate.

Giant viruses in the oceans : the 4th Algal Virus Workshop

Abstract: Giant double-stranded DNA viruses (such as record breaking Acanthamoeba polyphaga Mimivirus), with particle sizes of 0.2 to 0.6 micron, genomes of 300 kbp to 1.200 kbp, and commensurate complex gene contents, constitute an evolutionary mystery. They challenge the common vision of viruses, traditionally seen as highly streamlined genomes optimally fitted to the smallest possible -filterable- package. Such giant viruses are now discovered in increasing numbers through the systematic sampling of ocean waters as well as freshwater aquatic environments, where they play a significant role in controlling phyto- and bacterio- plankton populations. The 4th algal virus workshop showed that the study of these ecologically important viruses is now massively entering the genomic era, promising a better understanding of their diversity and, hopefully, some insights on their origin and the evolutionary forces that shaped their genomes.

Fluctuations in models of biological macroevolution

Abstract: Fluctuations in diversity and extinction sizes are discussed and compared for two different, individual-based models of biological coevolution. Both models display power-law distributions for various quantities of evolutionary interest, such as the lifetimes of individual species, the quiet periods between evolutionary upheavals larger than a given cutoff, and the sizes of extinction events. Time series of the diversity and measures of the size of extinctions give rise to flicker noise. Surprisingly, the power-law behaviors of the probability densities of quiet periods in the two models differ, while the distributions of the lifetimes of individual species are the same.

Evolutionary origin of power-laws in Biochemical Reaction Network; embedding abundance distribution into topology

Abstract: The evolutionary origin of universal statistics in biochemical reaction network is studied, to explain the power-law distribution of reaction links and the power-law distributions of chemical abundances. Using cell models with catalytic reaction network, we find evidence that the power-law distribution in abundances of chemicals emerges by the selection of cells with higher growth speeds. Through the further evolution, this inhomogeneity in chemical abundances is shown to be embedded in the distribution of links, leading to the power-law distribution. These findings provide novel insights into the nature of network evolution in living cells.

On the Genealogy of Populations: Trees, Branches and Offsprings

Abstract: We consider a neutral haploid population whose generations are not overlapping and whose size is large and constantly of $N$ individuals. Any generation is replaced by a new one and any individual has a single parent. We do not choose the stochastic rule which assigns the number of offsprings to any individual since results do not depend on the details of the dynamics, and, as a consequence, the model is parameter free. The genealogical tree is very complex, and distances between individuals (number of generations from the common ancestor) are distributed according to probability density which remains random in the thermodynamic limit (large population). We give a theoretical and numerical description of this distribution and we also consider the dynamical aspects of the problem describing the time evolution of the maximum and mean distances in a single population.

Simulation of Rapoport's rule for latitudinal species spread

Abstract: Rapoport's rule claims that latitudinal ranges of plant and animal species are generally smaller at low than at high latitudes. However, doubts as to the generality of the rule have been expressed, because studies providing evidence against the rule are more numerous than those in support of it. In groups for which support has been provided, the trend of increasing latitudinal ranges with latitude is restricted to or at least most distinct at high latitudes, suggesting that the effect may be a local phenomenon, for example the result of glaciations. Here we test the rule using two models, a simple one-dimensional one with a fixed number of animals expanding in a northern or southerly direction only, and the evolutionary/ecological Chowdhury model using birth, ageing, death, mutation, speciation, prey-predator relations and food levels. Simulations with both models gave results contradicting Rapoport's rule. In the first, latitudinal ranges were roughly independent of latitude, in the second, latitudinal ranges were greatest at low latitudes, as also shown empirically for some well studied groups of animals.

Divergent evolution paths of different genetic families in the Penna model

Abstract: We present some results of simulations of population growth and evolution, using the standard asexual Penna model, with individuals characterized by a string of bits representing a genome containing some possible mutations. After about 20000 simulation steps, when only a few genetic families are still present from among rich variety of families at the beginning of the simulation game, strong peaks in mutation distribution functions are observed. This known effect is due to evolution rules with hereditary mechanism. The birth and death balance in the simulation game also leads to elimination of families specified by different genomes. Number of families $G(t)$ versus time $t$ follow the power law, $G \propto t^n$. Our results show the power coefficient exponent $n$ is changing as the time goes. Starting from about --1, smoothly achieves about --2 after hundreds of steps, and finally has semi-smooth transition to 0, when only one family exists in the environment. This is in contrast with constant $n$ about --1 as found, for example, in \cite{bib:evolution}. We suspect that this discrepancy may be due to two different time scales in simulations - initial stages follow the $n\approx-1$ law, yet for large number of simulation steps we get $n\approx-2$, providing random initial population was sufficiently big to allow for still reliable statistical analysis. The $n\approx-1$ evolution stage seems to be associated with the Verhulst mechanism of population elimination due to the limited environmental capacity - when the standard evolution rules were modified, we observed a plateau ($n=0$) in the power law in short time scale, again followed by $n\approx -2$ law for longer times. The modified model uses birth rate controlled by the current population instead of the standard Verhulst death factor.

Evolutionary ecology in-silico: Does mathematical modelling help in understanding the "generic" trends?

Abstract: Motivated by the results of recent laboratory experiments (Yoshida et al. Nature, 424, 303-306 (2003)) as well as many earlier field observations that evolutionary changes can take place in ecosystems over relatively short ecological time scales, several ``unified'' mathematical models of evolutionary ecology have been developed over the last few years with the aim of describing the statistical properties of data related to the evolution of ecosystems. Moreover, because of the availability of sufficiently fast computers, it has become possible to carry out detailed computer simulations of these models. For the sake of completeness and to put these recent developments in the proper perspective, we begin with a brief summary of some older models of ecological phenomena and evolutionary processes. However, the main aim of this article is to review critically these ``unified'' models, particularly those published in the physics literature, in simple language that makes the new theories accessible to wider audience.

Ubiquity of synonymity: almost all large binary trees are not uniquely identified by their spectra or their immanantal polynomials

Abstract: There are several common ways to encode a tree as a matrix, such as the adjacency matrix, the Laplacian matrix (that is, the infinitesimal generator of the natural random walk), and the matrix of pairwise distances between leaves. Such representations involve a specific labeling of the vertices or at least the leaves, and so it is natural to attempt to identify trees by some feature of the associated matrices that is invariant under relabeling. An obvious candidate is the spectrum of eigenvalues (or, equivalently, the characteristic polynomial). We show for any of these choices of matrix that the fraction of binary trees with a unique spectrum goes to zero as the number of leaves goes to infinity. We investigate the rate of convergence of the above fraction to zero using numerical methods. For the adjacency and Laplacian matrices, we show that that the {\em a priori} more informative immanantal polynomials have no greater power to distinguish between trees.

Modeling the dynamics of inhomogeneous natural rotifer populations under toxicant exposure

Abstract: Most population models assume that individuals within a given population are identical, that is, the fundamental role of variation is ignored. Here we develop a general approach to modeling heterogeneous populations with discrete evolutionary time step. The theory is applied to population dynamics of natural rotifer populations. We show that under particular conditions the behavior of the inhomogeneous model possesses complex transition regimes, which depends both on the mean and the variance of the initial parameter distribution and the final state of the population depends on the least possible value from the domain of the parameter. The question of persistence of the population is discussed.

Speciational view of macroevolution: are micro and macroevolution decoupled?

Abstract: We introduce a simple computational model that, with a microscopic dynamics driven by natural selection and mutation alone, allows the description of true speciation events. A statistical analysis of the so generated evolutionary tree captures realistic features showing power laws for frequency distributions in time and size. Albeit these successful predictions, the difficulty in obtaining punctuated dynamics with mass extinctions suggests the necessity of decoupling micro and macro-evolutionary mechanisms in agreement with some ideas of Gould's and Eldredge's theory of punctuated equilibrium.

Sharp gene pool transition in a population affected by phenotype-based selective hunting

Abstract: We use a microscopic model of population dynamics, a modified version of the well known Penna model, to study some aspects of microevolution. This research is motivated by recent reports on the effect of selective hunting on the gene pool of bighorn sheep living in the Ram Mountain region, in Canada. Our model finds a sharp transition in the structure of the gene pool as some threshold for the number of animals hunted is reached.

The Mediterranean deep-sea fauna: historical evolution, bathymetric variations and geographical changes

Abstract: The deep-water fauna of the Mediterranean is characterized by an absence of distinctive characteristics and by a relative impoverishment. Both are a result of events after the Messinian salinity crisis (Late Miocene). The three main classes of phenomena involved in producing or recording these effects are analysed and discussed: - Historical: Sequential faunal changes during the Pliocene and thereafter in particular those during the Quaternary glaciations and still in progress. - Bathymetric: Changes in the vertical aspects of the Bathyal and Abyssal zones that took place under peculiar conditions, i.e. homothermy, a relative oligotrophy, the barrier of the Gibraltar sill, and water mass movement. The deeper the habitat of a species in the Mediterranean, the more extensive is its distribution elsewhere. - Geographical: There are strong affinities and relationships between Mediterranean and Atlantic faunas. Endemic species remain a biogeographical problem. Species always become smaller in size eastward where they occupy a progressively deeper habitat. Thus, the existing deep Mediterranean Sea appears to be younger than any other deep-sea constituent of the World Ocean.