Source author record

Simone Pigolotti

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

26 published item(s)

preprint2022arXiv

Coalescent dynamics of planktonic communities

Planktonic communities are extremely diverse and include a vast number of rare species. The dynamics of these rare species is best described by individual-based models. However, individual-based approaches to planktonic diversity face substantial difficulties, due to the large number of individuals required to make realistic predictions. In this paper, we study diversity of planktonic communities by means of a spatial coalescence model, that incorporates transport by oceanic currents. As a main advantage, our approach requires simulating a number of individuals equal to the size of the sample one is interested in, rather than the size of the entire community. By theoretical analysis and simulations, we explore the conditions upon which our coalescence model is equivalent to individual-based dynamics. As an application, we use our model to predict the impact of chaotic advection by oceanic currents on biodiversity. We conclude that the coalescent approach permits to simulate marine microbial communities much more efficiently than with individual-based models.

preprint2022arXiv

Population genetics in microchannels

Spatial constraints such as rigid barriers affect the dynamics of cell populations, potentially altering the course of natural evolution. In this paper, we investigate the population genetics of Escherichia coli proliferating in microchannels with open ends. Our analysis is based on a population model in which reproducing cells shift entire lanes of cells towards the open ends of the channel. The model predicts that diversity is lost very rapidly within lanes, but at a much slower pace among lanes. As a consequence, two mixed, neutral E. coli strains competing in a microchannel must organize into an ordered regular stripe pattern in the course of a few generations. These predictions are in quantitative agreement with our experiments. We also demonstrate that random mutations appearing in the middle of the channel are much more likely to reach fixation than those occurring elsewhere. Our results illustrate fundamental mechanisms of microbial evolution in spatially confined space.

preprint2016arXiv

Competition between fast- and slow-diffusing species in non-homogeneous environments

We study an individual-based model in which two spatially-distributed species, characterized by different diffusivities, compete for resources. We consider three different ecological settings. In the first, diffusing faster has a cost in terms of reproduction rate. In the second case, resources are not uniformly distributed in space. In the third case, the two species are transported by a fluid flow. In all these cases, at varying the parameters, we observe a transition from a regime in which diffusing faster confers an effective selective advantage to one in which it constitutes a disadvantage. We analytically estimate the magnitude of this advantage (or disadvantage) and test it by measuring fixation probabilities in simulations of the individual-based model. Our results provide a framework to quantify evolutionary pressure for increased or decreased dispersal in a given environment.

preprint2016arXiv

Protocols for copying and proofreading in template-assisted polymerization

We discuss how information encoded in a template polymer can be stochastically copied into a copy polymer. We consider four different stochastic copy protocols of increasing complexity, inspired by building blocks of the mRNA translation pathway. In the first protocol, monomer incorporation occurs in a single stochastic transition. We then move to a more elaborate protocol in which an intermediate step can be used for error correction. Finally, we discuss the operating regimes of two kinetic proofreading protocols: one in which proofreading acts from the final copying step, and one in which it acts from an intermediate step. We review known results for these models and, in some cases, extend them to analyze all possible combinations of energetic and kinetic discrimination. We show that, in each of these protocols, only a limited number of these combinations leads to an improvement of the overall copying accuracy.

preprint2015arXiv

Inherent directionality explains the lack of feedback loops in empirical networks

We explore the hypothesis that the relative abundance of feedback loops in many empirical complex networks is severely reduced owing to the presence of an inherent global directionality. Aimed at quantifying this idea, we propose a simple probabilistic model in which a free parameter $γ$ controls the degree of inherent directionality. Upon strengthening such directionality, the model predicts a drastic reduction in the fraction of loops which are also feedback loops. To test this prediction, we extensively enumerated loops and feedback loops in many empirical biological, ecological and socio- technological directed networks. We show that, in almost all cases, empirical networks have a much smaller fraction of feedback loops than network randomizations. Quite remarkably, this empirical finding is quantitatively reproduced, for all loop lengths, by our model by fitting its only parameter $γ$. Moreover, the fitted value of $γ$ correlates quite well with another direct measurement of network directionality, performed by means of a novel algorithm. We conclude that the existence of an inherent network directionality provides a parsimonious quantitative explanation for the observed lack of feedback loops in empirical networks.

preprint2015arXiv

Numerical and experimental study of the effects of noise on the permutation entropy

We analyze the effects of noise on the permutation entropy of dynamical systems. We take as numerical examples the logistic map and the Rössler system. Upon varying the noise strengthfaster, we find a transition from an almost-deterministic regime, where the permutation entropy grows slower than linearly with the pattern dimension, to a noise-dominated regime, where the permutation entropy grows faster than linearly with the pattern dimension. We perform the same analysis on experimental time-series by considering the stochastic spiking output of a semiconductor laser with optical feedback. Because of the experimental conditions, the dynamics is found to be always in the noise-dominated regime. Nevertheless, the analysis allows to detect regularities of the underlying dynamics. By comparing the results of these three different examples, we discuss the possibility of determining from a time series whether the underlying dynamics is dominated by noise or not.

preprint2015arXiv

Stochasticity enhances the gaining of bet-hedging strategies in contact-process-like dynamics

In biology and ecology, individuals or communities of individuals living in unpredictable environments often alternate between different evolutionary strategies to spread and reduce risks. Such behavior is commonly referred to as "bet-hedging." Long-term survival probabilities and population sizes can be much enhanced by exploiting such hybrid strategies. Here, we study the simplest possible birth-death stochastic model in which individuals can choose among a poor but safe strategy, a better but risky alternative, or a combination of both. We show analytically and computationally that the benefits derived from bet-hedging strategies are much stronger for higher environmental variabilities (large external noise) and/or for small spatial dimensions (large intrinsic noise). These circumstances are typically encountered by living systems, thus providing us with a possible justification for the ubiquitousness of bet-hedging in nature.

preprint2015arXiv

Thermodynamics of Error Correction

Information processing at the molecular scale is limited by thermal fluctuations. This can cause undesired consequences in copying information since thermal noise can lead to errors that can compromise the functionality of the copy. For example, a high error rate during DNA duplication can lead to cell death. Given the importance of accurate copying at the molecular scale, it is fundamental to understand its thermodynamic features. In this paper, we derive a universal expression for the copy error as a function of entropy production and {\cred work dissipated by the system during wrong incorporations}. Its derivation is based on the second law of thermodynamics, hence its validity is independent of the details of the molecular machinery, be it any polymerase or artificial copying device. Using this expression, we find that information can be copied in three different regimes. In two of them, work is dissipated to either increase or decrease the error. In the third regime, the protocol extracts work while correcting errors, reminiscent of a Maxwell demon. As a case study, we apply our framework to study a copy protocol assisted by kinetic proofreading, and show that it can operate in any of these three regimes. We finally show that, for any effective proofreading scheme, error reduction is limited by the chemical driving of the proofreading reaction.

preprint2014arXiv

Analysis of self-overlap reveals trade-offs in plankton swimming trajectories

Movement is a fundamental behaviour of organisms that brings about beneficial encounters with resources and mates, but at the same time exposes the organism to dangerous encounters with predators. The movement patterns adopted by organisms should reflect a balance between these contrasting processes. This trade-off can be hypothesized as being evident in the behaviour of plankton, which inhabit a dilute 3D environment with few refuges or orienting landmarks. We present an analysis of the swimming path geometries based on a volumetric Monte Carlo sampling approach, which is particularly adept at revealing such trade-offs by measuring the self-overlap of the trajectories. Application of this method to experimentally measured trajectories reveals that swimming patterns in copepods are shaped to efficiently explore volumes at small scales, while achieving a large overlap at larger scales. Regularities in the observed trajectories make the transition between these two regimes always sharper than in randomized trajectories or as predicted by random walk theory. Thus real trajectories present a stronger separation between exploration for food and exposure to predators. The specific scale and features of this transition depend on species, gender, and local environmental conditions, pointing at adaptation to state and stage dependent evolutionary trade-offs.

preprint2014arXiv

Selective advantage of diffusing faster

We study a stochastic spatial model of biological competition in which two species have the same birth and death rates, but different diffusion constants. In the absence of this difference, the model can be considered as an off-lattice version of the Voter model and presents similar coarsening properties. We show that even a relative difference in diffusivity on the order of a few percent may lead to a strong bias in the coarsening process favoring the more agile species. We theoretically quantify this selective advantage and present analytical formulas for the average growth of the fastest species and its fixation probability.

preprint2013arXiv

Kinetic vs. energetic discrimination in biological copying

We study stochastic copying schemes in which discrimination between a right and a wrong match is achieved via different kinetic barriers or different binding energies of the two matches. We demonstrate that, in single-step reactions, the two discrimination mechanisms are strictly alternative and can not be mixed to further reduce the error fraction. Close to the lowest error limit, kinetic discrimination results in a diverging copying velocity and dissipation per copied bit. On the opposite, energetic discrimination reaches its lowest error limit in an adiabatic regime where dissipation and velocity vanish. By analyzing experimentally measured kinetic rates of two DNA polymerases, T7 and Polγ, we argue that one of them operates in the kinetic and the other in the energetic regime. Finally, we show how the two mechanisms can be combined in copying schemes implementing error correction through a proofreading pathway

preprint2013arXiv

Replicator dynamics with turnover of players

We study adaptive dynamics in games where players abandon the population at a given rate, and are replaced by naive players characterized by a prior distribution over the admitted strategies. We demonstrate how such process leads macroscopically to a variant of the replicator equation, with an additional term accounting for player turnover. We study how Nash equilibria and the dynamics of the system are modified by this additional term, for prototypical examples such as the rock-scissor-paper game and different classes of two-action games played between two distinct populations. We conclude by showing how player turnover can account for non-trivial departures from Nash equilibria observed in data from lowest unique bid auctions.

preprint2013arXiv

Species abundances and lifetimes: from neutral to niche-stabilized communities

We study a stochastic community model able to interpolate from a neutral regime to a niche partitioned regime upon varying a single parameter tuning the intensity of niche stabilization, namely the difference between intraspecific and interspecific competition. By means of a self-consistent approach, we obtain an analytical expression for the species abundance distribution, in excellent agreement with stochastic simulations of the model. In the neutral limit, the Fisher log-series is recovered, while upon increasing the stabilization strength the species abundance distribution develops a maximum for species at intermediate abundances, corresponding to the emergence of a carrying capacity. Numerical studies of species extinction-time distribution show that niche-stabilization strongly affects also the dynamical properties of the system by increasing the average species lifetimes, while suppressing their fluctuations. The results are discussed in view of the niche-neutral debate and of their potential relevance to field data.

preprint2012arXiv

Equilibrium strategy and population-size effects in lowest unique bid auctions

In lowest unique bid auctions, $N$ players bid for an item. The winner is whoever places the \emph{lowest} bid, provided that it is also unique. We use a grand canonical approach to derive an analytical expression for the equilibrium distribution of strategies. We then study the properties of the solution as a function of the mean number of players, and compare them with a large dataset of internet auctions. The theory agrees with the data with striking accuracy for small population size $N$, while for larger $N$ a qualitatively different distribution is observed. We interpret this result as the emergence of two different regimes, one in which adaptation is feasible and one in which it is not. Our results question the actual possibility of a large population to adapt and find the optimal strategy when participating in a collective game.

preprint2012arXiv

Growth, competition and cooperation in spatial population genetics

We study an individual based model describing competition in space between two different alleles. Although the model is similar in spirit to classic models of spatial population genetics such as the stepping stone model, here however space is continuous and the total density of competing individuals fluctuates due to demographic stochasticity. By means of analytics and numerical simulations, we study the behavior of fixation probabilities, fixation times, and heterozygosity, in a neutral setting and in cases where the two species can compete or cooperate. By concluding with examples in which individuals are transported by fluid flows, we argue that this model is a natural choice to describe competition in marine environments.

preprint2012arXiv

Population dynamics in compressible flows

Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible turbulence is relevant for population dynamics in the ocean and we consider cases both where the velocity field is turbulent and when it is static. Furthermore, we investigate populations in terms of a continuos density field and when the populations are treated via discrete particles. In the last case we focus on the competition and fixation of one species compared to another

preprint2012arXiv

Quality Control System Response to Stochastic Growth of Amyloid Fibrils

We introduce a stochastic model describing aggregation of misfolded proteins and degradation by the protein quality control system in a single cell. In analogy with existing literature, aggregates can grow, nucleate and fragment stochastically. We assume that the quality control system acts as an enzyme that can degrade aggregates at different stages of the growth process, with an efficiency that decreases with the size of the aggregate. We show how this stochastic dynamics, depending on the parameter choice, leads to two qualitatively different behaviors: a homeostatic state, where the quality control system is stable and aggregates of large sizes are not formed, and an oscillatory state, where the quality control system periodically breaks down, allowing for the formation of large aggregates. We discuss how these periodic breakdowns may constitute a mechanism for the sporadic development of neurodegenerative diseases.

preprint2012arXiv

Switching between oscillations and homeostasis in competing negative and positive feedback motifs

We analyze a class of network motifs in which a short, two-node positive feed- back motif is inserted in a three-node negative feedback loop. We demonstrate that such networks can undergo a bifurcation to a state where a stable fixed point and a stable limit cycle coexist. At the bifurcation point the period of the oscillations diverges. Further, intrinsic noise can make the system switch between oscillatory state and the stationary state spontaneously. We find that this switching also occurs in previous models of circadian clocks that use this combination of positive and negative feedback. Our results suggest that real- life circadian systems may need specific regulation to prevent or minimize such switching events.

preprint2012arXiv

What ecological factors shape species-area curves in neutral models?

Understanding factors that shape biodiversity and species coexistence across scales is of utmost importance in ecology, both theoretically and for conservation policies. Species-area relationships (SARs), measuring how the number of observed species increases upon enlarging the sampled area, constitute a convenient tool for quantifying the spatial structure of biodiversity. While general features of species-area curves are quite universal across ecosystems, some quantitative aspects can change significantly. Several attempts have been made to link these variations to ecological forces. Within the framework of spatially explicit neutral models, here we scrutinize the effect of varying the local population size (i.e. the number of individuals per site) and the level of habitat saturation (allowing for empty sites). We conclude that species-area curves become shallower when the local population size increases, while habitat saturation, unless strongly violated, plays a marginal role. Our findings provide a plausible explanation of why SARs for microorganisms are flatter than those for larger organisms.

preprint2011arXiv

Population genetics in compressible flows

We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We then study three examples of compressible advecting fields: a shell model of turbulence, a sinusoidal velocity field and a linear velocity sink. In all cases, advection leads to a striking drop in the fixation time, as well as a large reduction in the global carrying capacity. Despite localization on convergence zones, one species goes extinct much more rapidly than in well-mixed populations. For a weak harmonic potential, one finds a bimodal distribution of fixation times. The long-lived states in this case are demixed configurations with a single boundary, whose location depends on the fitness advantage.

preprint2010arXiv

Coexistence and invasibility in a two-species competition model with habitat-preference

The outcome of competition among species is influenced by the spatial distribution of species and effects such as demographic stochasticity, immigration fluxes, and the existence of preferred habitats. We introduce an individual-based model describing the competition of two species and incorporating all the above ingredients. We find that the presence of habitat preference --- generating spatial niches --- strongly stabilizes the coexistence of the two species. Eliminating habitat preference --- neutral dynamics --- the model generates patterns, such as distribution of population sizes, practically identical to those obtained in the presence of habitat preference, provided an higher immigration rate is considered. Notwithstanding the similarity in the population distribution, we show that invasibility properties depend on habitat preference in a non-trivial way. In particular, the neutral model results results more invasible or less invasible depending on whether the comparison is made at equal immigration rate or at equal distribution of population size, respectively. We discuss the relevance of these results for the interpretation of invasibility experiments and the species occupancy of preferred habitats.

preprint2010arXiv

Self-consistent method for density estimation

The estimation of a density profile from experimental data points is a challenging problem, usually tackled by plotting a histogram. Prior assumptions on the nature of the density, from its smoothness to the specification of its form, allow the design of more accurate estimation procedures, such as Maximum Likelihood. Our aim is to construct a procedure that makes no explicit assumptions, but still providing an accurate estimate of the density. We introduce the self-consistent estimate: the power spectrum of a candidate density is given, and an estimation procedure is constructed on the assumption, to be released \emph{a posteriori}, that the candidate is correct. The self-consistent estimate is defined as a prior candidate density that precisely reproduces itself. Our main result is to derive the exact expression of the self-consistent estimate for any given dataset, and to study its properties. Applications of the method require neither priors on the form of the density nor the subjective choice of parameters. A cutoff frequency, akin to a bin size or a kernel bandwidth, emerges naturally from the derivation. We apply the self-consistent estimate to artificial data generated from various distributions and show that it reaches the theoretical limit for the scaling of the square error with the dataset size.

preprint2009arXiv

How Gaussian competition leads to lumpy or uniform species distributions

A central model in theoretical ecology considers the competition of a range of species for a broad spectrum of resources. Recent studies have shown that essentially two different outcomes are possible. Either the species surviving competition are more or less uniformly distributed over the resource spectrum, or their distribution is 'lumped' (or 'clumped'), consisting of clusters of species with similar resource use that are separated by gaps in resource space. Which of these outcomes will occur crucially depends on the competition kernel, which reflects the shape of the resource utilization pattern of the competing species. Most models considered in the literature assume a Gaussian competition kernel. This is unfortunate, since predictions based on such a Gaussian assumption are not robust. In fact, Gaussian kernels are a border case scenario, and slight deviations from this function can lead to either uniform or lumped species distributions. Here we illustrate the non-robustness of the Gaussian assumption by simulating different implementations of the standard competition model with constant carrying capacity. In this scenario, lumped species distributions can come about by secondary ecological or evolutionary mechanisms or by details of the numerical implementation of the model. We analyze the origin of this sensitivity and discuss it in the context of recent applications of the model.

preprint2006arXiv

Absorbing processes in Richardson diffusion: analytical results

We consider the recently addressed problem of a passive particle (a predator), being the center of a ``sphere of interception'' of radius $R$ and able to absorb other passive particles (the preys) entering into the sphere. Assuming that all the particles are advected by a turbulent flow and that, in particular, the Richardson equation properly describes the relative dispersion, we calculate an analytical expression for the flux into the sphere as a function of time, assuming an initial constant density of preys outside the sphere. In the same framework, we show that the distribution of times of first passage into the sphere has a $t^{-5/2}$ power law tail, seen in contrast to the $t^{-3/2}$ appearing in standard 3D diffusion. We also discuss the correction due to the integral length scale on the results in the stationary case.

preprint2005arXiv

Boltzmann entropy and chaos in a large assembly of weakly interacting systems

We introduce a high dimensional symplectic map, modeling a large system consisting of weakly interacting chaotic subsystems, as a toy model to analyze the interplay between single-particle chaotic dynamics and particles interactions in thermodynamic systems. We study the growth with time of the Boltzmann entropy, S_B, in this system as a function of the coarse graining resolution. We show that a characteristic scale emerges, and that the behavior of S_B vs t, at variance with the Gibbs entropy, does not depend on the coarse graining resolution, as far as it is finer than this scale. The interaction among particles is crucial to achieve this result, while the rate of entropy growth depends essentially on the single-particle chaotic dynamics (for t not too small). It is possible to interpret the basic features of the dynamics in terms of a suitable Markov approximation.

preprint2004arXiv

A Stochastic Model for the Species Abundance Problem in an Ecological Community

We propose a model based on coupled multiplicative stochastic processes to understand the dynamics of competing species in an ecosystem. This process can be conveniently described by a Fokker-Planck equation. We provide an analytical expression for the marginalized stationary distribution. Our solution is found in excellent agreement with numerical simulations and compares rather well with observational data from tropical forests.