Showing 1–42 of 42 results for author: Fontanari, J F

Solving the prisoner's dilemma trap in Hamilton's model of temporarily formed random groups

Authors: José F. Fontanari, Mauro Santos

Abstract: Explaining the evolution of cooperation in the strong altruism scenario, where a cooperator does not benefit from her contribution to the public goods, is a challenging problem that requires positive assortment among cooperators (i.e., cooperators must tend to associate with other cooperators) or punishment of defectors. The need for these drastic measures stems from the analysis of a group select… ▽ More Explaining the evolution of cooperation in the strong altruism scenario, where a cooperator does not benefit from her contribution to the public goods, is a challenging problem that requires positive assortment among cooperators (i.e., cooperators must tend to associate with other cooperators) or punishment of defectors. The need for these drastic measures stems from the analysis of a group selection model of temporarily formed random groups introduced by Hamilton nearly fifty years ago to describe the fate of altruistic behavior in a population. Challenging conventional wisdom, we show analytically here that strong altruism evolves in Hamilton's original model in the case of biparental sexual reproduction. Moreover, when the cost of cooperation is small and the amplified contribution shared by group members is large, cooperation is the only stable strategy in equilibrium. Thus, our results provide a solution to the `problem of origination' of strong altruism, i.e. how cooperation can take off from an initial low frequency of cooperators. We discuss a possible reassessment of cooperation in cases of viral co-infection, as cooperation may even be favored in situations where the prisoner's dilemma applies. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Revisiting institutional punishment in the $N$-person prisoner's dilemma

Authors: Bianca Y. S. Ishikawa, José F. Fontanari

Abstract: The conflict between individual and collective interests makes fostering cooperation in human societies a challenging task, requiring drastic measures such as the establishment of sanctioning institutions. These institutions are costly because they have to be maintained regardless of the presence or absence of offenders. Here, we propose realistic improvements to the standard $N$-person prisoner's… ▽ More The conflict between individual and collective interests makes fostering cooperation in human societies a challenging task, requiring drastic measures such as the establishment of sanctioning institutions. These institutions are costly because they have to be maintained regardless of the presence or absence of offenders. Here, we propose realistic improvements to the standard $N$-person prisoner's dilemma formulation with institutional punishment by eliminating overpunishment, requiring a minimum number of contributors to establish the sanctioning institution, and sharing the cost among them once this minimum number is reached. In addition, we focus on large groups or communities for which sanctioning institutions are ubiquitous. Using the replicator equation framework for an infinite population, we find that by sufficiently fining players who fail to contribute either to the public good or to the sanctioning institution, a population of contributors immune to invasion by these free riders can be established, provided that the contributors are sufficiently numerous. In a finite population, we use finite-size scaling to show that, for some parameter settings, demographic noise helps to fixate the strategy that contributes to the public good but not to the sanctioning institution even for infinitely large populations when, somewhat counterintuitively, its proportion in the initial population vanishes with a small power of the population size. △ Less

Submitted 13 June, 2024; v1 submitted 9 June, 2024; originally announced June 2024.

Imitation dynamics and the replicator equation

Authors: José F. Fontanari

Abstract: Evolutionary game theory has impacted many fields of research by providing a mathematical framework for studying the evolution and maintenance of social and moral behaviors. This success is owed in large part to the demonstration that the central equation of this theory - the replicator equation - is the deterministic limit of a stochastic imitation (social learning) dynamics. Here we offer an alt… ▽ More Evolutionary game theory has impacted many fields of research by providing a mathematical framework for studying the evolution and maintenance of social and moral behaviors. This success is owed in large part to the demonstration that the central equation of this theory - the replicator equation - is the deterministic limit of a stochastic imitation (social learning) dynamics. Here we offer an alternative elementary proof of this result, which holds for the scenario where players compare their instantaneous (not average) payoffs to decide whether to maintain or change their strategies, and only more successful individuals can be imitated. △ Less

Submitted 23 May, 2024; v1 submitted 31 March, 2024; originally announced April 2024.

Journal ref: Europhysics Letters 146 (2024) 47001

A soluble model for synchronized rhythmic activity in ant colonies

Authors: Pedro M. M. da Silveira, José F. Fontanari

Abstract: Synchronization is one of the most striking instances of collective behavior, occurring in many natural phenomena. For example, in some ant species, ants are inactive within the nest most of the time, but their bursts of activity are highly synchronized and involve the entire nest population. Here we revisit a simulation model that generates this synchronized rhythmic activity through autocatalyti… ▽ More Synchronization is one of the most striking instances of collective behavior, occurring in many natural phenomena. For example, in some ant species, ants are inactive within the nest most of the time, but their bursts of activity are highly synchronized and involve the entire nest population. Here we revisit a simulation model that generates this synchronized rhythmic activity through autocatalytic behavior, i.e., active ants can activate inactive ants, followed by a period of rest. We derive a set of delay differential equations that provide an accurate description of the short- and medium-time simulations for large ant colonies. Analysis of the fixed-point solutions, complemented by numerical integration of the equations, indicates the existence of stable limit-cycle solutions when the rest period is greater than a threshold and the event of spontaneous activation of inactive ants is very unlikely, so that most of the arousal of ants is done by active ants. Furthermore, we argue that the persistent oscillations observed in the simulations for colonies of finite size are due to resonant amplification of demographic noise. △ Less

Submitted 31 May, 2024; v1 submitted 17 February, 2024; originally announced February 2024.

Journal ref: Mathematical Biosciences 375 (2024) 109245

The dynamics of casual groups can keep free-riders at bay

Authors: José F. Fontanari, Mauro Santos

Abstract: Understanding the conditions for maintaining cooperation in groups of unrelated individuals despite the presence of non-cooperative members is a major research topic in contemporary biological, sociological, and economic theory. The $N$-person snowdrift game models the type of social dilemma where cooperative actions are costly, but there is a reward for performing them. We study this game in a sc… ▽ More Understanding the conditions for maintaining cooperation in groups of unrelated individuals despite the presence of non-cooperative members is a major research topic in contemporary biological, sociological, and economic theory. The $N$-person snowdrift game models the type of social dilemma where cooperative actions are costly, but there is a reward for performing them. We study this game in a scenario where players move between play groups following the casual group dynamics, where groups grow by recruiting isolates and shrink by losing individuals who then become isolates. This describes the size distribution of spontaneous human groups and also the formation of sleeping groups in monkeys. We consider three scenarios according to the probability of isolates joining a group. We find that for appropriate choices of the cost-benefit ratio of cooperation and the aggregation-disaggregation ratio in the formation of casual groups, free-riders can be completely eliminated from the population. If individuals are more attracted to large groups, we find that cooperators persist in the population even when the mean group size diverges. We also point out the remarkable similarity between the replicator equation approach to public goods games and the trait group formulation of structured demes. △ Less

Submitted 7 March, 2024; v1 submitted 28 January, 2024; originally announced January 2024.

Journal ref: Mathematical Biosciences 372 (2024) 109188

Evolving division of labor in a response threshold model

Authors: José F. Fontanari, Viviane M. de Oliveira, Paulo R. A. Campos

Abstract: The response threshold model explains the emergence of division of labor (i.e., task specialization) in an unstructured population by assuming that the individuals have different propensities to work on different tasks. The incentive to attend to a particular task increases when the task is left unattended and decreases when individuals work on it. Here we derive mean-field equations for the stimu… ▽ More The response threshold model explains the emergence of division of labor (i.e., task specialization) in an unstructured population by assuming that the individuals have different propensities to work on different tasks. The incentive to attend to a particular task increases when the task is left unattended and decreases when individuals work on it. Here we derive mean-field equations for the stimulus dynamics and show that they exhibit complex attractors through period-doubling bifurcation cascades when the noise disrupting the thresholds is small. In addition, we show how the fixed threshold can be set to ensure specialization in both the transient and equilibrium regimes of the stimulus dynamics. However, a complete explanation of the emergence of division of labor requires that we address the question of where the threshold variation comes from, starting from a homogeneous population. We then study a structured population scenario, where the population is divided into a large number of independent groups of equal size, and the fitness of a group is proportional to the weighted mean work performed on the tasks during a fixed period of time. Using a winner-take-all strategy to model group competition and assuming an initial homogeneous metapopulation, we find that a substantial fraction of workers specialize in each task, without the need to penalize task switching. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Journal ref: Ecological Complexity 58 (2024) 101083

The synergy between two threats: disinformation and Covid-19

Authors: Henrique A. Tórtura, José F. Fontanari

Abstract: The breakdown of trusted sources of information is probably one of the most serious problems today, since in the absence of a common ground, it will be impossible to address the problems that trouble our contemporary world. The Covid-19 pandemic is just a recent situation where the lack of agreed stances has led to failure and hopelessness. In fact, disinformation surrounding the Covid-19 has been… ▽ More The breakdown of trusted sources of information is probably one of the most serious problems today, since in the absence of a common ground, it will be impossible to address the problems that trouble our contemporary world. The Covid-19 pandemic is just a recent situation where the lack of agreed stances has led to failure and hopelessness. In fact, disinformation surrounding the Covid-19 has been a distinctive feature of this pandemic since its very beginning and has hampered what is perhaps the most important initiative to prevent the spread of the coronavirus, viz., an effective communication between scientifically-minded health authorities and the general public. To investigate how disinformation threatens epistemic security, here we propose and solve analytically an evolutionary game-theoretic model where the individuals must accurately estimate some property of their hazardous environment. They can either explore the environment or copy the estimate from another individual, who may display a distorted version of its estimate. We find that the exploration-only strategy is optimal when the environment is relatively safe and the individuals are not reliable. In this doomsday scenario, disinformation erodes trust and suppresses the ability of the individuals to share information with one another. △ Less

Submitted 24 September, 2022; v1 submitted 16 May, 2022; originally announced May 2022.

Journal ref: Mathematical Models and Methods in Applied Sciences 32 (2022) 2077-2097

Local adaptation, phenotypic plasticity, and species coexistence

Authors: José F. Fontanari, Margarida Matos, Mauro Santos

Abstract: Understanding the mechanisms of species coexistence has always been a fundamental topic in ecology. Classical theory predicts that interspecific competition may select for traits that stabilize niche differences, although recent work shows that this is not strictly necessary. Here we ask whether adaptive phenotypic plasticity could allow species coexistence (i.e., some stability at an equilibrium… ▽ More Understanding the mechanisms of species coexistence has always been a fundamental topic in ecology. Classical theory predicts that interspecific competition may select for traits that stabilize niche differences, although recent work shows that this is not strictly necessary. Here we ask whether adaptive phenotypic plasticity could allow species coexistence (i.e., some stability at an equilibrium point) without ecological differentiation in habitat use. We used individual-based stochastic simulations defining a landscape composed of spatially uncorrelated or autocorrelated environmental patches, where two species with the same competitive strategies, not able to coexist without some form of phenotypic plasticity, expanded their ranges in the absence of a competition-colonization trade-off (a well-studied mechanism for species diversity). Each patch is characterized by a random environmental value that determines the optimal phenotype of its occupants. In such a scenario, only local adaptation and gene flow (migration) may interact to promote genetic variation and coexistence in the metapopulation. Results show that a competitively inferior species with adaptive phenotypic plasticity can coexist in a same patch with a competitively superior, non-plastic species, provided the migration rates and variances of the patches' environmental values are sufficiently large. △ Less

Submitted 10 March, 2023; v1 submitted 29 April, 2022; originally announced May 2022.

Journal ref: Front. Ecol. Evol. 11 (2023) 1077374

The surprising little effectiveness of cooperative algorithms in parallel problem solving

Authors: Sandro M. Reia, Larissa F. Aquino, José F. Fontanari

Abstract: Biological and cultural inspired optimization algorithms are nowadays part of the basic toolkit of a great many research domains. By mimicking processes in nature and animal societies, these general-purpose search algorithms promise to deliver optimal or near-optimal solutions using hardly any information on the optimization problems they are set to tackle. Here we study the performances of a cult… ▽ More Biological and cultural inspired optimization algorithms are nowadays part of the basic toolkit of a great many research domains. By mimicking processes in nature and animal societies, these general-purpose search algorithms promise to deliver optimal or near-optimal solutions using hardly any information on the optimization problems they are set to tackle. Here we study the performances of a cultural-inspired algorithm -- the imitative learning search -- as well as of asexual and sexual variants of evolutionary algorithms in finding the global maxima of NK-fitness landscapes. The main performance measure is the total number of agent updates required by the algorithms to find those global maxima and the baseline performance, which establishes the effectiveness of the cooperative algorithms, is set by the blind search in which the agents explore the problem space (binary strings) by flipping bits at random. We find that even for smooth landscapes that exhibit a single maximum, the evolutionary algorithms do not perform much better than the blind search due to the stochastic effects of the genetic roulette. The imitative learning is immune to this effect thanks to the deterministic choice of the fittest string in the population, which is used as a model for imitation. The tradeoff is that for rugged landscapes the imitative learning search is more prone to be trapped in local maxima than the evolutionary algorithms. In fact, in the case of rugged landscapes with a mild density of local maxima, the blind search either beats or matches the cooperative algorithms regardless of whether the task is to find the global maximum or to find the fittest state within a given runtime. △ Less

Submitted 29 May, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Journal ref: Eur. Phys. J. B (2020) 93: 140

Influence of technological progress and renewability on the sustainability of ecosystem engineers populations

Authors: Guilherme M. Lopes, José F. Fontanari

Abstract: Overpopulation and environmental degradation due to inadequate resource-use are outcomes of human's ecosystem engineering that has profoundly modified the world's landscape. Despite the age-old concern that unchecked population and economic growth may be unsustainable, the prospect of societal collapse remains contentious today. Contrasting with the usual approach to modeling human-nature interact… ▽ More Overpopulation and environmental degradation due to inadequate resource-use are outcomes of human's ecosystem engineering that has profoundly modified the world's landscape. Despite the age-old concern that unchecked population and economic growth may be unsustainable, the prospect of societal collapse remains contentious today. Contrasting with the usual approach to modeling human-nature interactions, which are based on the Lotka-Volterra predator-prey model with humans as the predators and nature as the prey, here we address this issue using a discrete-time population dynamics model of ecosystem engineers. The growth of the population of engineers is modeled by the Beverton-Holt equation with a density-dependent carrying capacity that is proportional to the number of usable habitats. These habitats (e.g., farms) are the products of the work of the individuals on the virgin habitats (e.g., native forests), hence the denomination engineers of ecosystems to those agents. The human-made habitats decay into degraded habitats, which eventually regenerate into virgin habitats. For slow regeneration resources, we find that the dynamics is dominated by cycles of prosperity and collapse, in which the population reaches vanishing small densities. However, increase of the efficiency of the engineers to explore the resources eliminates the dangerous cyclical patterns of feast and famine and leads to a stable equilibrium that balances population growth and resource availability. This finding supports the viewpoint of growth optimists that technological progress may avoid collapse. △ Less

Submitted 28 December, 2018; originally announced December 2018.

Journal ref: Mathematical Biosciences and Engineering 16, 3450-3464 ( 2019)

Mobility helps problem-solving systems to avoid Groupthink

Authors: Paulo F. Gomes, Sandro M. Reia, Francisco A. Rodrigues, José F. Fontanari

Abstract: Groupthink occurs when everyone in a group starts thinking alike, as when people put unlimited faith in a leader. Avoiding this phenomenon is a ubiquitous challenge to problem-solving enterprises and typical countermeasures involve the mobility of group members. Here we use an agent-based model of imitative learning to study the influence of the mobility of the agents on the time they require to f… ▽ More Groupthink occurs when everyone in a group starts thinking alike, as when people put unlimited faith in a leader. Avoiding this phenomenon is a ubiquitous challenge to problem-solving enterprises and typical countermeasures involve the mobility of group members. Here we use an agent-based model of imitative learning to study the influence of the mobility of the agents on the time they require to find the global maxima of NK-fitness landscapes. The agents cooperate by exchanging information on their fitness and use this information to copy the fittest agent in their influence neighborhoods, which are determined by face-to-face interaction networks. The influence neighborhoods are variable since the agents perform random walks in a two-dimensional space. We find that mobility is slightly harmful for solving easy problems, i.e. problems that do not exhibit suboptimal solutions or local maxima. For difficult problems, however, mobility can prevent the imitative search being trapped in suboptimal solutions and guarantees a better performance than the independent search for any system size. △ Less

Submitted 7 January, 2019; v1 submitted 7 August, 2018; originally announced August 2018.

Journal ref: Phys. Rev. E 99, 032301 (2019)

Predictability of the imitative learning trajectories

Authors: Paulo R. A. Campos, José F. Fontanari

Abstract: The fitness landscape metaphor plays a central role on the modeling of optimizing principles in many research fields, ranging from evolutionary biology, where it was first introduced, to management research. Here we consider the ensemble of trajectories of the imitative learning search, in which agents exchange information on their fitness and imitate the fittest agent in the population aiming at… ▽ More The fitness landscape metaphor plays a central role on the modeling of optimizing principles in many research fields, ranging from evolutionary biology, where it was first introduced, to management research. Here we consider the ensemble of trajectories of the imitative learning search, in which agents exchange information on their fitness and imitate the fittest agent in the population aiming at reaching the global maximum of the fitness landscape. We assess the degree to which the starting and ending points determine the learning trajectories using two measures, namely, the predictability that yields the probability that two randomly chosen trajectories are the same, and the mean path divergence that gauges the dissimilarity between two learning trajectories. We find that the predictability is greater in rugged landscapes than in smooth ones. The mean path divergence, however, is strongly affected by the search parameters -- population size and imitation propensity -- that obliterate the influence of the underlying landscape. The learning trajectories become more deterministic, in the sense that there are fewer distinct trajectories and those trajectories are more similar to each other, with increasing population size and imitation propensity. In addition, we find that the roughness of the learning trajectories, which measures the deviation from additivity of the fitness function, is always greater than the roughness estimated over the entire fitness landscape. △ Less

Submitted 15 November, 2018; v1 submitted 12 July, 2018; originally announced July 2018.

Journal ref: J. Stat. Mech. (2019) 013501

The collapse of ecosystem engineer populations

Authors: José F. Fontanari

Abstract: Humans are the ultimate ecosystem engineers who have profoundly transformed the world's landscapes in order to enhance their survival. Somewhat paradoxically, however, sometimes the unforeseen effect of this ecosystem engineering is the very collapse of the population it intended to protect. Here we use a spatial version of a standard population dynamics model of ecosystem engineers to study the c… ▽ More Humans are the ultimate ecosystem engineers who have profoundly transformed the world's landscapes in order to enhance their survival. Somewhat paradoxically, however, sometimes the unforeseen effect of this ecosystem engineering is the very collapse of the population it intended to protect. Here we use a spatial version of a standard population dynamics model of ecosystem engineers to study the colonization of unexplored virgin territories by a small settlement of engineers. We find that during the expansion phase the population density reaches values much higher than those the environment can support in the equilibrium situation. When the colonization front reaches the boundary of the available space, the population density plunges sharply and attains its equilibrium value. The collapse takes place without warning and happens just after the population reaches its peak number. We conclude that overpopulation and the consequent collapse of an expanding population of ecosystem engineers is a natural consequence of the nonlinear feedback between the population and environment variables. △ Less

Submitted 26 November, 2017; originally announced November 2017.

Journal ref: Mathematics (2018) 6: 9

The revival of the Baldwin Effect

Authors: José F. Fontanari, Mauro Santos

Abstract: The idea that a genetically fixed behavior evolved from the once differential learning ability of individuals that performed the behavior is known as the Baldwin effect. A highly influential paper [Hinton G.E., Nowlan S.J., 1987. How learning can guide evolution. Complex Syst. 1, 495-502] claimed that this effect can be observed in silico, but here we argue that what was actually shown is that the… ▽ More The idea that a genetically fixed behavior evolved from the once differential learning ability of individuals that performed the behavior is known as the Baldwin effect. A highly influential paper [Hinton G.E., Nowlan S.J., 1987. How learning can guide evolution. Complex Syst. 1, 495-502] claimed that this effect can be observed in silico, but here we argue that what was actually shown is that the learning ability is easily selected for. Then we demonstrate the Baldwin effect to happen in the in silico scenario by estimating the probability and waiting times for the learned behavior to become innate. Depending on parameter values, we find that learning can increase the chance of fixation of the learned behavior by several orders of magnitude compared with the non-learning situation. △ Less

Submitted 9 August, 2017; v1 submitted 27 February, 2017; originally announced February 2017.

Journal ref: Eur. Phys. J. B (2017) 90: 186

The spatial dynamics of ecosystem engineers

Authors: Caroline Franco, José F. Fontanari

Abstract: The changes on abiotic features of ecosystems have rarely been taken into account by population dynamics models, which typically focus on trophic and competitive interactions between species. However, understanding the population dynamics of organisms that must modify their habitats in order to survive, the so-called ecosystem engineers, requires the explicit incorporation of abiotic interactions… ▽ More The changes on abiotic features of ecosystems have rarely been taken into account by population dynamics models, which typically focus on trophic and competitive interactions between species. However, understanding the population dynamics of organisms that must modify their habitats in order to survive, the so-called ecosystem engineers, requires the explicit incorporation of abiotic interactions in the models. Here we study a model of ecosystem engineers that is discrete both in space and time, and where the engineers and their habitats are arranged in patches fixed to the sites of regular lattices. The growth of the engineer population is modeled by Ricker equation with a density-dependent carrying capacity that is given by the number of modified habitats. A diffusive dispersal stage ensures that a fraction of the engineers move from their birth patches to neighboring patches. We find that dispersal influences the metapopulation dynamics only in the case that the local or single-patch dynamics exhibits chaotic behavior. In that case, it can suppress the chaotic behavior and avoid extinctions in the regime of large intrinsic growth rate of the population. △ Less

Submitted 10 August, 2017; v1 submitted 28 November, 2016; originally announced November 2016.

Journal ref: Mathematical Biosciences 292, 76-85 (2017)

Phenotypic Plasticity, the Baldwin Effect, and the Speeding up of Evolution: the Computational Roots of an Illusion

Authors: Mauro Santos, Eörs Szathmáry, José F. Fontanari

Abstract: An increasing number of dissident voices claim that the standard neo-Darwinian view of genes as 'leaders' and phenotypes as 'followers' during the process of adaptive evolution should be turned on its head. This idea is older than the rediscovery of Mendel's laws of inheritance and has been given several names before its final 'Baldwin effect' label. A condition for this effect is that environment… ▽ More An increasing number of dissident voices claim that the standard neo-Darwinian view of genes as 'leaders' and phenotypes as 'followers' during the process of adaptive evolution should be turned on its head. This idea is older than the rediscovery of Mendel's laws of inheritance and has been given several names before its final 'Baldwin effect' label. A condition for this effect is that environmentally induced variation such as phenotypic plasticity or learning is crucial for the initial establishment of a population. This gives the necessary time for natural selection to act on genetic variation and the adaptive trait can be eventually encoded in the genotype. An influential paper published in the late 1980s showed the Baldwin effect to happen in computer simulations, and claimed that it was crucial to solve a difficult adaptive task. This generated much excitement among scholars in various disciplines that regard neo-Darwinian accounts to explain the evolutionary emergence of high-order phenotypic traits such as consciousness or language almost hopeless. Here, we use analytical and computational approaches to show that a standard population genetics treatment can easily crack what the scientific community has granted as an unsolvable adaptive problem without learning. The Baldwin effect is once again in need of convincing theoretical foundations. △ Less

Submitted 25 November, 2014; originally announced November 2014.

Journal ref: Journal of Theoretical Biology 371, 127-136 (2015)

Effect of migration in a diffusion model for template coexistence in protocells

Authors: Jose F. Fontanari, Maurizio Serva

Abstract: The compartmentalization of distinct templates in protocells and the exchange of templates between them (migration) are key elements of a modern scenario for prebiotic evolution. Here we use the diffusion approximation of population genetics to study analytically the steady-state properties of such prebiotic scenario. The coexistence of distinct template types inside a protocell is achieved by a s… ▽ More The compartmentalization of distinct templates in protocells and the exchange of templates between them (migration) are key elements of a modern scenario for prebiotic evolution. Here we use the diffusion approximation of population genetics to study analytically the steady-state properties of such prebiotic scenario. The coexistence of distinct template types inside a protocell is achieved by a selective pressure at the protocell level (group selection) favoring protocells with a mixed template composition. In the degenerate case, where the templates have the same replication rate, we find that a vanishingly small migration rate suffices to eliminate the segregation effect of random drift and so to promote coexistence. In the non-degenerate case, a small migration rate greatly boosts coexistence as compared with the situation where there is no migration. However, increase of the migration rate beyond a critical value leads to the complete dominance of the more efficient template type (homogeneous regime). In this case, we find a continuous phase transition separating the homogeneous and the coexistence regimes, with the order parameter vanishing linearly with the distance to the transition point. △ Less

Submitted 15 February, 2014; originally announced February 2014.

Journal ref: Bulletin of Mathematical Biology 76 (2014) 654-672

Nonlinear group survival in Kimura's model for the evolution of altruism

Authors: Jose F. Fontanari, Maurizio Serva

Abstract: Establishing the conditions that guarantee the spreading or the sustenance of altruistic traits in a population is the main goal of intergroup selection models. Of particular interest is the balance of the parameters associated to group size, migration and group survival against the selective advantage of the non-altruistic individuals. Here we use Kimura's diffusion model of intergroup selection… ▽ More Establishing the conditions that guarantee the spreading or the sustenance of altruistic traits in a population is the main goal of intergroup selection models. Of particular interest is the balance of the parameters associated to group size, migration and group survival against the selective advantage of the non-altruistic individuals. Here we use Kimura's diffusion model of intergroup selection to determine those conditions in the case the group survival probability is a nonlinear non-decreasing function of the proportion of altruists in a group. In the case this function is linear, there are two possible steady states which correspond to the non-altruistic and the altruistic phases. At the discontinuous transition line separating these phases there is a non-ergodic coexistence phase. For a continuous concave survival function, we find an ergodic coexistence phase that occupies a finite region of the parameter space in between the altruistic and the non-altruistic phases, and is separated from these phases by continuous transition lines. For a convex survival function, the coexistence phase disappears altogether but a bistable phase appears for which the choice of the initial condition determines whether the evolutionary dynamics leads to the altruistic or the non-altruistic steady state. △ Less

Submitted 24 January, 2014; v1 submitted 23 July, 2013; originally announced July 2013.

Journal ref: Mathematical Biosciences 249 (2014), 18-26

Reinforcement and inference in cross-situational word learning

Authors: Paulo F. C. Tilles, Jose F. Fontanari

Abstract: Cross-situational word learning is based on the notion that a learner can determine the referent of a word by finding something in common across many observed uses of that word. Here we propose an adaptive learning algorithm that contains a parameter that controls the strength of the reinforcement applied to associations between concurrent words and referents, and a parameter that regulates infere… ▽ More Cross-situational word learning is based on the notion that a learner can determine the referent of a word by finding something in common across many observed uses of that word. Here we propose an adaptive learning algorithm that contains a parameter that controls the strength of the reinforcement applied to associations between concurrent words and referents, and a parameter that regulates inference, which includes built-in biases, such as mutual exclusivity, and information of past learning events. By adjusting these parameters so that the model predictions agree with data from representative experiments on cross-situational word learning, we were able to explain the learning strategies adopted by the participants of those experiments in terms of a trade-off between reinforcement and inference. These strategies can vary wildly depending on the conditions of the experiments. For instance, for fast mapping experiments (i.e., the correct referent could, in principle, be inferred in a single observation) inference is prevalent, whereas for segregated contextual diversity experiments (i.e., the referents are separated in groups and are exhibited with members of their groups only) reinforcement is predominant. Other experiments are explained with more balanced doses of reinforcement and inference. △ Less

Submitted 15 July, 2013; originally announced July 2013.

Journal ref: Front. Behav. Neurosci. 7:163 (2013)

Solvable model for template coexistence in protocells

Authors: J. F. Fontanari, M. Serva

Abstract: Compartmentalization of self-replicating molecules (templates) in protocells is a necessary step towards the evolution of modern cells. However, coexistence between distinct template types inside a protocell can be achieved only if there is a selective pressure favoring protocells with a mixed template composition. Here we study analytically a group selection model for the coexistence between two… ▽ More Compartmentalization of self-replicating molecules (templates) in protocells is a necessary step towards the evolution of modern cells. However, coexistence between distinct template types inside a protocell can be achieved only if there is a selective pressure favoring protocells with a mixed template composition. Here we study analytically a group selection model for the coexistence between two template types using the diffusion approximation of population genetics. The model combines competition at the template and protocell levels as well as genetic drift inside protocells. At the steady state, we find a continuous phase transition separating the coexistence and segregation regimes, with the order parameter vanishing linearly with the distance to the critical point. In addition, we derive explicit analytical expressions for the critical steady-state probability density of protocell compositions. △ Less

Submitted 26 January, 2013; v1 submitted 17 October, 2012; originally announced October 2012.

Journal ref: Europhysics Letters, 101 (2013) 38006

Minimal model of associative learning for cross-situational lexicon acquisition

Authors: Paulo F. C. Tilles, Jose F. Fontanari

Abstract: An explanation for the acquisition of word-object mappings is the associative learning in a cross-situational scenario. Here we present analytical results of the performance of a simple associative learning algorithm for acquiring a one-to-one mapping between $N$ objects and $N$ words based solely on the co-occurrence between objects and words. In particular, a learning trial in our learning scena… ▽ More An explanation for the acquisition of word-object mappings is the associative learning in a cross-situational scenario. Here we present analytical results of the performance of a simple associative learning algorithm for acquiring a one-to-one mapping between $N$ objects and $N$ words based solely on the co-occurrence between objects and words. In particular, a learning trial in our learning scenario consists of the presentation of $C + 1 < N$ objects together with a target word, which refers to one of the objects in the context. We find that the learning times are distributed exponentially and the learning rates are given by $\ln{[\frac{N(N-1)}{C + (N-1)^{2}}]}$ in the case the $N$ target words are sampled randomly and by $\frac{1}{N} \ln [\frac{N-1}{C}] $ in the case they follow a deterministic presentation sequence. This learning performance is much superior to those exhibited by humans and more realistic learning algorithms in cross-situational experiments. We show that introduction of discrimination limitations using Weber's law and forgetting reduce the performance of the associative algorithm to the human level. △ Less

Submitted 17 December, 2012; v1 submitted 6 April, 2012; originally announced April 2012.

Journal ref: J. Math. Psych. 56, 396-403 (2012)

A structural model of emotions of cognitive dissonances

Authors: Jose F. Fontanari, Marie-Claude Bonniot-Cabanac, Michel Cabanac, Leonid I. Perlovsky

Abstract: Cognitive dissonance is the stress that comes from holding two conflicting thoughts simultaneously in the mind, usually arising when people are asked to choose between two detrimental or two beneficial options. In view of the well-established role of emotions in decision making, here we investigate whether the conventional structural models used to represent the relationships among basic emotions,… ▽ More Cognitive dissonance is the stress that comes from holding two conflicting thoughts simultaneously in the mind, usually arising when people are asked to choose between two detrimental or two beneficial options. In view of the well-established role of emotions in decision making, here we investigate whether the conventional structural models used to represent the relationships among basic emotions, such as the Circumplex model of affect, can describe the emotions of cognitive dissonance as well. We presented a questionnaire to 34 anonymous participants, where each question described a decision to be made among two conflicting motivations and asked the participants to rate analogically the pleasantness and the intensity of the experienced emotion. We found that the results were compatible with the predictions of the Circumplex model for basic emotions. △ Less

Submitted 28 February, 2012; originally announced February 2012.

Journal ref: Neural Networks 32 (2012) 57-64

The optimization and shock waves in evolution dynamics

Authors: David B. Saakian, Jose F. Fontanari

Abstract: We consider the optimal dynamics in the infinite population evolution models with general symmetric fitness landscape. The search of optimal evolution trajectories are complicated due to sharp transitions (like shock waves) in evolution dynamics with smooth fitness landscapes, which exist even in case of popular quadratic fitness. We found exact analytical solutions for discontinuous dynamics at… ▽ More We consider the optimal dynamics in the infinite population evolution models with general symmetric fitness landscape. The search of optimal evolution trajectories are complicated due to sharp transitions (like shock waves) in evolution dynamics with smooth fitness landscapes, which exist even in case of popular quadratic fitness. We found exact analytical solutions for discontinuous dynamics at the large genome length limit. We found the optimal mutation rates for the fixed fitness landscape. The single peak fitness landscape gives the fastest dynamics to send the vast majority of the population from the initial sequence to the neighborhood of the final sequence. △ Less

Submitted 17 August, 2009; originally announced August 2009.

Comments: 6 pages

The information capacity of hypercycles

Authors: Daniel A. M. M. Silvestre, Jose F. Fontanari

Abstract: Hypercycles are information integration systems which are thought to overcome the information crisis of prebiotic evolution by ensuring the coexistence of several short templates. For imperfect template replication, we derive a simple expression for the maximum number of distinct templates $n_m$ that can coexist in a hypercycle and show that it is a decreasing function of the length $L$ of the t… ▽ More Hypercycles are information integration systems which are thought to overcome the information crisis of prebiotic evolution by ensuring the coexistence of several short templates. For imperfect template replication, we derive a simple expression for the maximum number of distinct templates $n_m$ that can coexist in a hypercycle and show that it is a decreasing function of the length $L$ of the templates. In the case of high replication accuracy we find that the product $n_m L$ tends to a constant value, limiting thus the information content of the hypercycle. Template coexistence is achieved either as a stationary equilibrium (stable fixed point) or a stable periodic orbit in which the total concentration of functional templates is nonzero. For the hypercycle system studied here we find numerical evidence that the existence of an unstable fixed point is a necessary condition for the presence of periodic orbits. △ Less

Submitted 12 June, 2008; originally announced June 2008.

Comments: 4 pages, 1 figure, submitted to J. Theor. Biol

Package models and the information crisis of prebiotic evolution

Authors: Daniel A. M. M. Silvestre, José F. Fontanari

Abstract: The coexistence between different types of templates has been the choice solution to the information crisis of prebiotic evolution, triggered by the finding that a single RNA-like template cannot carry enough information to code for any useful replicase. In principle, confining $d$ distinct templates of length $L$ in a package or protocell, whose survival depends on the coexistence of the templa… ▽ More The coexistence between different types of templates has been the choice solution to the information crisis of prebiotic evolution, triggered by the finding that a single RNA-like template cannot carry enough information to code for any useful replicase. In principle, confining $d$ distinct templates of length $L$ in a package or protocell, whose survival depends on the coexistence of the templates it holds in, could resolve this crisis provided that $d$ is made sufficiently large. Here we review the prototypical package model of Niesert et al. 1981 which guarantees the greatest possible region of viability of the protocell population, and show that this model, and hence the entire package approach, does not resolve the information crisis. This is so because to secure survival the total information content of the protocell, $Ld$, must tend to a constant value that depends only on the spontaneous error rate per nucleotide of the template replication mechanism. As a result, an increase of $d$ must be followed by a decrease of $L$ to ensure the protocell viability, so that the net information gain is null. △ Less

Submitted 17 October, 2007; originally announced October 2007.

Comments: 11 pages, two columns, 11 figures, submitted to J. Theor. Biol

Preservation of information in a prebiotic package model

Authors: Daniel A. M. M. Silvestre, Jose F. Fontanari

Abstract: The coexistence between different informational molecules has been the preferred mode to circumvent the limitation posed by imperfect replication on the amount of information stored by each of these molecules. Here we reexamine a classic package model in which distinct information carriers or templates are forced to coexist within vesicles, which in turn can proliferate freely through binary div… ▽ More The coexistence between different informational molecules has been the preferred mode to circumvent the limitation posed by imperfect replication on the amount of information stored by each of these molecules. Here we reexamine a classic package model in which distinct information carriers or templates are forced to coexist within vesicles, which in turn can proliferate freely through binary division. The combined dynamics of vesicles and templates is described by a multitype branching process which allows us to write equations for the average number of the different types of vesicles as well as for their extinction probabilities. The threshold phenomenon associated to the extinction of the vesicle population is studied quantitatively using finite-size scaling techniques. We conclude that the resultant coexistence is too frail in the presence of parasites and so confinement of templates in vesicles without an explicit mechanism of cooperation does not resolve the information crisis of prebiotic evolution. △ Less

Submitted 19 April, 2007; v1 submitted 25 December, 2006; originally announced December 2006.

Comments: 9 pages, 8 figures, accepted version, to be published in PRE

Revisiting the nonequilibrium phase transition of the triplet-creation model

Authors: Giovano O. Cardozo, Jose F. Fontanari

Abstract: The nonequilibrium phase transition in the triplet-creation model is investigated using critical spreading and the conservative diffusive contact process. The results support the claim that at high enough diffusion the phase transition becomes discontinuous. As the diffusion probability increases the critical exponents change continuously from the ordinary directed percolation (DP) class to the… ▽ More The nonequilibrium phase transition in the triplet-creation model is investigated using critical spreading and the conservative diffusive contact process. The results support the claim that at high enough diffusion the phase transition becomes discontinuous. As the diffusion probability increases the critical exponents change continuously from the ordinary directed percolation (DP) class to the compact directed percolation (CDP). The fractal dimension of the critical cluster, however, switches abruptly between those two universality classes. Strong crossover effects in both methods make it difficult, if not impossible, to establish the exact location of the tricritical point. △ Less

Submitted 26 July, 2006; originally announced July 2006.

Comments: 7 pages, 12 figures

Template coexistence in prebiotic vesicle models

Authors: Daniel G. M. Silvestre, Jose F. Fontanari

Abstract: The coexistence of distinct templates is a common feature of the diverse proposals advanced to resolve the information crisis of prebiotic evolution. However, achieving robust template coexistence turned out to be such a difficult demand that only a class of models, the so-called package models, seems to have met it so far. Here we apply Wright's Island formulation of group selection to study th… ▽ More The coexistence of distinct templates is a common feature of the diverse proposals advanced to resolve the information crisis of prebiotic evolution. However, achieving robust template coexistence turned out to be such a difficult demand that only a class of models, the so-called package models, seems to have met it so far. Here we apply Wright's Island formulation of group selection to study the conditions for the coexistence of two distinct template types confined in packages (vesicles) of finite capacity. In particular, we show how selection acting at the level of the vesicles can neutralize the pressures towards the fixation of any one of the template types (random drift) and of the type with higher replication rate (deterministic competition). We give emphasis to the role of the distinct generation times of templates and vesicles as yet another obstacle to coexistence. △ Less

Submitted 30 November, 2006; v1 submitted 25 May, 2005; originally announced May 2005.

Comments: 7 pages, 8 figures

Journal ref: Eur. Phys. J. B 47, 423-429 (2005)

Multifractal analysis of DNA walks and trails

Authors: Alexandre Rosas, Edvaldo Nogueira Jr., José F. Fontanari

Abstract: The characterization of the long-range order and fractal properties of DNA sequences has proved a difficult though highly rewarding task due mainly to the mosaic character of DNA consisting of many interwoven patches of various lengths with different nucleotide constitutions. We apply here a recently proposed generalization of the detrended fluctuation analysis method to show that the DNA walk c… ▽ More The characterization of the long-range order and fractal properties of DNA sequences has proved a difficult though highly rewarding task due mainly to the mosaic character of DNA consisting of many interwoven patches of various lengths with different nucleotide constitutions. We apply here a recently proposed generalization of the detrended fluctuation analysis method to show that the DNA walk construction, in which the DNA sequence is viewed as a time series, exhibits a monofractal structure regardless of the existence of local trends in the series. In addition, we point out that the monofractal structure of the DNA walks carries over to an apparently alternative graphical construction given by the projection of the DNA walk into the $d$ spatial coordinates, termed DNA trails. In particular, we calculate the fractal dimension $D_t$ of the DNA trails using a well-known result of fractal theory linking $D_t$ to the Hurst exponent $H$ of the corresponding DNA walk. Comparison with estimates obtained by the standard box-counting method allows the evaluation of both finite-length and local trends effects. △ Less

Submitted 9 October, 2002; v1 submitted 17 September, 2002; originally announced September 2002.

Comments: 1 - We have changed the title of the paper, which reads now ``Multifractal analysis of DNA walks and trails''. 2 - We have clarified the meanings of $s(n)$, $s_ν(n)$, $m$, etc.. 3 - We have now pointed out in all places that the DNA walks exhibit a monofractal structure regardless of the existence of local trends in the series. 6 pages, 5 figures

Complementarity and diversity in a soluble model ecosystem

Authors: Viviane M. de Oliveira, J. F. Fontanari

Abstract: Complementarity among species with different traits is one of the basic processes affecting biodiversity, defined as the number of species in the ecosystem. We present here a soluble model ecosystem in which the species are characterized by binary traits and their pairwise interactions follow a complementarity principle. Manipulation of the species composition, and so the study of its effects on… ▽ More Complementarity among species with different traits is one of the basic processes affecting biodiversity, defined as the number of species in the ecosystem. We present here a soluble model ecosystem in which the species are characterized by binary traits and their pairwise interactions follow a complementarity principle. Manipulation of the species composition, and so the study of its effects on the species diversity is achieved through the introduction of a bias parameter favoring one of the traits. Using statistical mechanics tools we find explicit expressions for the allowed values of the equilibrium species concentrations in terms of the control parameters of the model. △ Less

Submitted 19 March, 2002; originally announced March 2002.

Phase-Transitions in a model for the formation of herpes simplex ulcers

Authors: Claudia Pio Ferreira, Jose Fernando Fontanari, Rita M. Zorzenon dos Santos

Abstract: The critical properties of a cellular automaton model describing the spreading of infection of the Herpes Simplex Virus in corneal tissue are investigated through the dynamic Monte Carlo method. The model takes into account different cell susceptibilities to the viral infection, as suggested by experimental findings. In a two-dimensional square lattice, the sites are associated to two distinct t… ▽ More The critical properties of a cellular automaton model describing the spreading of infection of the Herpes Simplex Virus in corneal tissue are investigated through the dynamic Monte Carlo method. The model takes into account different cell susceptibilities to the viral infection, as suggested by experimental findings. In a two-dimensional square lattice, the sites are associated to two distinct types of cells, namely, permissive and resistant to the infection. While a permissive cell becomes infected in the presence of a single infected cell in its neighborhood, a resistant cell needs to be surrounded by at least R>1 infected or dead cells in order to become infected. The infection is followed by the death of the cells resulting in ulcers whose forms may be dendritic (self-limited clusters) or amoeboid (percolating clusters) depending on the degree of resistance R of the resistant cells as well as on the density of permissive cells in the healthy tissue. We show that a phase transition between these two regimes occurs only for R>=5 and, in addition, that the phase-transition is in the universality class of the ordinary percolation. △ Less

Submitted 13 June, 2001; originally announced June 2001.

Nonequilibrium phase transitions in a model for the origin of life

Authors: Claudia P. Ferreira, J. F. Fontanari

Abstract: The requisites for the persistence of small colonies of self-replicating molecules living in a two-dimensional lattice are investigated analytically in the infinite diffusion or mean-field limit and through Monte Carlo simulations in the position-fixed or contact process limit. The molecules are modeled by hipercyclic replicators A which are capable of replicating via binary fission A + E --2A w… ▽ More The requisites for the persistence of small colonies of self-replicating molecules living in a two-dimensional lattice are investigated analytically in the infinite diffusion or mean-field limit and through Monte Carlo simulations in the position-fixed or contact process limit. The molecules are modeled by hipercyclic replicators A which are capable of replicating via binary fission A + E --2A with production rates s as well as via catalytically assisted replication 2A + E --3A with rate c. In addition, a molecule can degrade into its source materials E with rate $γ$. In the asymptotic regime the population can be characterized by the presence (active phase) and the absence (empty phase) of replicators in the lattice. In both diffusion regimes, we find that for small values of the ratio c/$γ$ these phases are separated by a second-order phase transition which is in the universality class of the directed percolation, while for small values s/$gamma$ the phase transition is of first order. Furthermore, we illustrate the suitability of the dynamic Monte Carlo method, which is based on the analysis of the spreading behaviour of a few active cells in the center of an otherwise infinite empty lattice, to adress the problem of emergence of replicators. Rather surprisingly, we show that this method allows an unambiguous identification of the order of the nonequilibrium phase transition. △ Less

Submitted 24 April, 2001; originally announced April 2001.

Extinctions in the random replicator model

Authors: Viviane M. de Oliveira, J. F. Fontanari

Abstract: The statistical properties of an ecosystem composed of species interacting via pairwise, random interactions and deterministic, concentration limiting self-interaction are studied analytically with tools of equilibrium statistical mechanics of disordered systems. Emphasis is given to the effects of externally induced extinction of a fixed fraction of species at the outset of the evolutionary pro… ▽ More The statistical properties of an ecosystem composed of species interacting via pairwise, random interactions and deterministic, concentration limiting self-interaction are studied analytically with tools of equilibrium statistical mechanics of disordered systems. Emphasis is given to the effects of externally induced extinction of a fixed fraction of species at the outset of the evolutionary process. The manner the ecosystem copes with the initial extinction event depends on the degree of competition among the species as well as on the strength of that event. For instance, in the regime of high competition the ecosystem diversity, given by the fraction of surviving species, is practically insensitive to the strengt of the initial extinction provided it is not too large, while in the less competitive regime the diversity decreases linearly with the size of the event. In the case of large extinction events we find that no further biotic extinctions take place and, furthermore, that rare species become very unlikely to be found in the ecosystem at equilibrium. In addition, we show that the reciprocal of the Edwards-Anderson order parameter yields a good measure of the diversity of the model ecosystem. △ Less

Submitted 24 April, 2001; originally announced April 2001.

Group selection models in prebiotic evolution

Authors: D. Alves, P. R. A. Campos, A. T. C. Silva, J. F. Fontanari

Abstract: The evolution of enzyme production is studied analytically using ideas of the group selection theory for the evolution of altruistic behavior. In particular, we argue that the mathematical formulation of Wilson's structured deme model ({\it The Evolution of Populations and Communities}, Benjamin/Cumings, Menlo Park, 1980) is a mean-field approach in which the actual environment that a particular… ▽ More The evolution of enzyme production is studied analytically using ideas of the group selection theory for the evolution of altruistic behavior. In particular, we argue that the mathematical formulation of Wilson's structured deme model ({\it The Evolution of Populations and Communities}, Benjamin/Cumings, Menlo Park, 1980) is a mean-field approach in which the actual environment that a particular individual experiences is replaced by an {\it average} environment. That formalism is further developed so as to avoid the mean-field approximation and then applied to the problem of enzyme production in the prebiotic context, where the enzyme producer molecules play the altruists role while the molecules that benefit from the catalyst without paying its production cost play the non-altruists role. The effects of synergism (i.e., division of labor) as well as of mutations are also considered and the results of the equilibrium analysis are summarized in phase diagrams showing the regions of the space of parameters where the altruistic, non-altruistic and the coexistence regimes are stable. In general, those regions are delimitated by discontinuous transition lines which end at critical points. △ Less

Submitted 11 September, 2000; originally announced September 2000.

Comments: 22 pages, 10 figures

Random replicators with high-order interactions

Authors: Viviane M. de Oliveira, J. F. Fontanari

Abstract: We use tools of the equilibrium statistical mechanics of disordered systems to study analytically the statistical properties of an ecosystem composed of N species interacting via random, Gaussian interactions of order p >= 2, and deterministic self-interactions u <= 0. We show that for nonzero u the effect of increasing the order of the interactions is to make the system more cooperative, in the… ▽ More We use tools of the equilibrium statistical mechanics of disordered systems to study analytically the statistical properties of an ecosystem composed of N species interacting via random, Gaussian interactions of order p >= 2, and deterministic self-interactions u <= 0. We show that for nonzero u the effect of increasing the order of the interactions is to make the system more cooperative, in the sense that the fraction of extinct species is greatly reduced. Furthermore, we find that for p > 2 there is a threshold value which gives a lower bound to the concentration of the surviving species, preventing then the existence of rare species and, consequently, increasing the robustness of the ecosystem to external perturbations. △ Less

Submitted 4 May, 2000; originally announced May 2000.

Comments: 7 pages, 4 Postscript figures

On the structure pf genealogical trees in the presence of selection

Authors: P. R. A. Campos, M. T. Sonoda, J. F. Fontanari

Abstract: We investigate through numerical simulations the effect of selection on two summary statistics for nucleotide variation in a sample of two genes from a population of N asexually reproducing haploid individuals. One is the mean time since two individuals had their most recent common ancestor ($\bar{T_s}$), and the other is the mean number of nucleotide differences between two genes in the sample… ▽ More We investigate through numerical simulations the effect of selection on two summary statistics for nucleotide variation in a sample of two genes from a population of N asexually reproducing haploid individuals. One is the mean time since two individuals had their most recent common ancestor ($\bar{T_s}$), and the other is the mean number of nucleotide differences between two genes in the sample ($\bar{d_s}$). In the case of diminishing epistasis, in which the deleterious effect of a new mutation is attenuated, we find that the scale of $\bar{d_s}$ with the population size depends on the mutation rate, leading then to the onset of a sharp threshold phenomenon as N becomes large. △ Less

Submitted 1 November, 1999; originally announced October 1999.

Comments: 6 pages

Error Propagation in the Hypercycle

Authors: P. R. A. Campos, J. F. Fontanari, P. F. Stadler

Abstract: We study analytically the steady-state regime of a network of n error-prone self-replicating templates forming an asymmetric hypercycle and its error tail. We show that the existence of a master template with a higher non-catalyzed self-replicative productivity, a, than the error tail ensures the stability of chains in which m<n-1 templates coexist with the master species. The stability of these… ▽ More We study analytically the steady-state regime of a network of n error-prone self-replicating templates forming an asymmetric hypercycle and its error tail. We show that the existence of a master template with a higher non-catalyzed self-replicative productivity, a, than the error tail ensures the stability of chains in which m<n-1 templates coexist with the master species. The stability of these chains against the error tail is guaranteed for catalytic coupling strengths (K) of order of a. We find that the hypercycle becomes more stable than the chains only for K of order of a2. Furthermore, we show that the minimal replication accuracy per template needed to maintain the hypercycle, the so-called error threshold, vanishes like sqrt(n/K) for large K and n<=4. △ Less

Submitted 29 July, 1999; originally announced July 1999.

Journal ref: Phys. Rev. E 61: 2996-3002 (2000)

Stochastic group selection model for the evolution of altruism

Authors: Ana T. C. Silva, J. F. Fontanari

Abstract: We study numerically and analytically a stochastic group selection model in which a population of asexually reproducing individuals, each of which can be either altruist or non-altruist, is subdivided into $M$ reproductively isolated groups (demes) of size $N$. The cost associated with being altruistic is modelled by assigning the fitness $1- τ$, with $τ\in [0,1]$, to the altruists and the fitne… ▽ More We study numerically and analytically a stochastic group selection model in which a population of asexually reproducing individuals, each of which can be either altruist or non-altruist, is subdivided into $M$ reproductively isolated groups (demes) of size $N$. The cost associated with being altruistic is modelled by assigning the fitness $1- τ$, with $τ\in [0,1]$, to the altruists and the fitness 1 to the non-altruists. In the case that the altruistic disadvantage $τ$ is not too large, we show that the finite $M$ fluctuations are small and practically do not alter the deterministic results obtained for $M \to \infty$. However, for large $τ$ these fluctuations greatly increase the instability of the altruistic demes to mutations. These results may be relevant to the dynamics of parasite-host systems and, in particular, to explain the importance of mutation in the evolution of parasite virulence. △ Less

Submitted 25 January, 1999; originally announced January 1999.

Comments: 12 pages, 7 figures

Finite-size scaling of the error threshold transition in finite population

Authors: P. R. A. Campos, J. F. Fontanari

Abstract: The error threshold transition in a stochastic (i.e. finite population) version of the quasispecies model of molecular evolution is studied using finite-size scaling. For the single-sharp-peak replication landscape, the deterministic model exhibits a first-order transition at $Q=Q_c=1/a$, where $% Q$ is the probability of exact replication of a molecule of length $L \to \infty$, and $a$ is the s… ▽ More The error threshold transition in a stochastic (i.e. finite population) version of the quasispecies model of molecular evolution is studied using finite-size scaling. For the single-sharp-peak replication landscape, the deterministic model exhibits a first-order transition at $Q=Q_c=1/a$, where $% Q$ is the probability of exact replication of a molecule of length $L \to \infty$, and $a$ is the selective advantage of the master string. For sufficiently large population size, $N$, we show that in the critical region the characteristic time for the vanishing of the master strings from the population is described very well by the scaling assumption $τ= N^{1/2} f_a \left [ \left (Q - Q_c) N^{1/2} \right ] $, where $f_a$ is an $a$-dependent scaling function. △ Less

Submitted 15 September, 1998; originally announced September 1998.

Comments: 8 pages, 3 ps figures. submitted to J. Phys. A

Error threshold in finite populations

Authors: Domingos Alves, J. F. Fontanari

Abstract: A simple analytical framework to study the molecular quasispecies evolution of finite populations is proposed, in which the population is assumed to be a random combination of the constiyuent molecules in each generation,i.e., linkage disequilibrium at the population level is neglected. In particular, for the single-sharp-peak replication landscape we investigate the dependence of the error thre… ▽ More A simple analytical framework to study the molecular quasispecies evolution of finite populations is proposed, in which the population is assumed to be a random combination of the constiyuent molecules in each generation,i.e., linkage disequilibrium at the population level is neglected. In particular, for the single-sharp-peak replication landscape we investigate the dependence of the error threshold on the population size and find that the replication accuracy at threshold increases linearly with the reciprocal of the population size for sufficiently large populations. Furthermore, in the deterministic limit our formulation yields the exact steady-state of the quasispecies model, indicating then the population composition is a random combination of the molecules. △ Less

Submitted 31 March, 1998; v1 submitted 23 September, 1997; originally announced September 1997.

Comments: 14 pages and 4 figures

Journal ref: Phys. Rev. E, vol.57, number 6, june 1998

Error threshold in the evolution of diploid organisms

Authors: Domingos Alves, J. F. Fontanari

Abstract: The effects of error propagation in the reproduction of diploid organisms are studied within the populational genetics framework of the quasispecies model. The dependence of the error threshold on the dominance parameter is fully investigated. In particular, it is shown that dominance can protect the wild-type alleles from the error catastrophe. The analysis is restricted to a diploid analogue o… ▽ More The effects of error propagation in the reproduction of diploid organisms are studied within the populational genetics framework of the quasispecies model. The dependence of the error threshold on the dominance parameter is fully investigated. In particular, it is shown that dominance can protect the wild-type alleles from the error catastrophe. The analysis is restricted to a diploid analogue of the single-peaked landscape. △ Less

Submitted 9 December, 1996; originally announced December 1996.

Comments: 9 pages, 4 Postscript figures. Submitted to J. Phy. A: Mat. and Gen

Journal ref: J. Phys. A: Math and Gen., v. 30, (1997), p. 2601-2607

A Population Genetic Approach to the Quasispecies Model

Authors: Domingos Alves, Jose Fernando Fontanari

Abstract: A population genetics formulation of Eigen's molecular quasispecies model is proposed and several simple replication landscapes are investigated analytically. Our results show a remarcable similarity to those obtained with the original kinetics formulation of the quasispecies model. However, due to the simplicity of our approach, the space of the parameters that define the model can be explored.… ▽ More A population genetics formulation of Eigen's molecular quasispecies model is proposed and several simple replication landscapes are investigated analytically. Our results show a remarcable similarity to those obtained with the original kinetics formulation of the quasispecies model. However, due to the simplicity of our approach, the space of the parameters that define the model can be explored. In particular, for the simgle-sharp-peak landscape our analysis yelds some interesting predictions such as the existence of a maximum peak height and a mini- mum molecule length for the onset of the error threshold transition. △ Less

Submitted 31 May, 1996; v1 submitted 24 May, 1996; originally announced May 1996.

Comments: 16 pages, 4 Postscript figures. Submited to Phy. Rev. E