Evolution of migration under kin selection and local adaptation

We present here a stochastic two-locus, two-habitat model for the evolution of migration with local adaptation and kin selection. One locus determines the migration rate while the other causes local adaptation. We show that the opposing forces of kin competition and local adaptation can lead to the existence of one or two convergence stable migration rates, notably depending on the recombination rate between the two loci. We show that linkage between migration and local adaptation loci has two antagonist effects: when linkage is tight, cost of local adaptation increases, leading to smaller equilibrium migration rates. However, when linkage is tighter, the population structure at the migration locus tends to be very high because of the indirect selection, and thus equilibrium migration rates increases. This result, qualitatively different from results obtained with other models of migration evolution, indicates that ignoring drift or the detail of the genetic architecture may lead to incorrect conclusions.     

Evolution of migration rate in a spatially realistic metapopulation model

We use an individual-based, spatially realistic metapopulation model to study the evolution of migration rate. We first explore the consequences of habitat change in hypothetical patch networks on a regular lattice. If the primary consequence of habitat change is an increase in local extinction risk as a result of decreased local population sizes, migration rate increases. A nonmonotonic response, with migration rate decreasing at high extinction rate, was obtained only by assuming very frequent catastrophes. If the quality of the matrix habitat deteriorates, leading to increased mortality during migration, the evolutionary response is more complex. As long as habitat patch occupancy does not decrease markedly with increased migration mortality, reduced migration rate evolves. However, once mortality becomes so high that empty patches remain uncolonized for a long time, evolution tends to increase migration rate, which may lead to an "evolutionary rescue" in a fragmented landscape. Kin competition has a quantitative effect on the evolution of migration rate in our model, but these patterns in the evolution of migration rate appear to be primarily caused by spatiotemporal variation in fitness and mortality during migration. We apply the model to real habitat patch networks occupied by two checkerspot butterfly (Melitaea) species, for which sufficient data are available to estimate rigorously most of the model parameters. The model-predicted migration rate is not significantly different from the empirically observed one. Regional variation in patch areas and connectivities leads to regional variation in the optimal migration rate, predictions that can be tested empirically.     

Influence of finite-sites mutation, population subdivision and sampling schemes on patterns of nucleotide polymorphism for species with molecular hyperdiversity

Molecular hyperdiversity has been documented in viruses, prokaryotes and eukaryotes. Such organisms undermine the assumptions of the infinite-sites mutational model, because multiple mutational events at a site comprise a non-negligible portion of polymorphisms. Moreover, different sampling schemes of individuals from species with subdivided populations can profoundly influence resulting patterns and interpretations of molecular variation. Inspired by molecular hyperdiversity in the nematode Caenorhabditis sp. 5, which exhibits average pairwise differences among synonymous sites of >5% as well as modest population structure, we investigated via coalescent simulation the joint effects of a finite-sites mutation (FSM) process and population subdivision on the variant frequency spectrum. From many demes interconnected through a stepping-stone migration model, we constructed local samples from a single deme, pooled samples from several demes and scattered samples of a single individual from numerous demes. Compared with a single panmictic population at equilibrium, we find that high population mutation rates induce a deficit of rare variants (positive Tajima's D) under a FSM model. Population structure also induces such a skew for local samples when migration is high and for pooled samples when migration is low. Contrasts of sampling schemes for C. sp. 5 imply high mutational input coupled with high migration. We propose that joint analysis of local, pooled and scattered samples for species with subdivided populations provides a means of improving inference of demographic history, by virtue of the partially distinct patterns of polymorphism that manifest when sequences are analyzed according to differing sampling schemes.     

Effects of environmental heterogeneity on victim-exploiter coevolution

We study victim-exploiter coevolution in a spatially heterogeneous island model. In each species, fitness consequences of between-species interactions are controlled by a single haploid diallelic locus. Our emphasis is on the conditions for the maintenance of genetic variation, the dynamic patterns observed, the extent of local adaptation and genetic differentiation between different demes, and on how different parameters (such as the strength and heterogeneity in selection, migration rates, and the number of sites) affect the dynamic and static behavior of the system. We show that under spatially homogeneous selection the maintenance of genetic variation is possible through asynchronous nonlinear dynamics where the allele frequencies in a majority of demes quickly synchronize but the rest do not. Spatially heterogeneous selection can maintain genetic variation even if migration rates are maximal. This happens in an oscillatory way. Genetic variation is most likely to be maintained at high levels if the heterogeneity in selection is large. If there are some restrictions on migration, genetic variation can be maintained at a stable equilibrium. This behavior is most likely at intermediate migration rates. In this case, the system can exhibit high spatial subdivision as measured by F(ST) values but relatively low local adaptation.     

Towards a correct description of zooplankton feeding in models: Taking into account food-mediated unsynchronized vertical migration

Complex nature of foraging behaviour of zooplankton makes it difficult to describe adequately zooplankton grazing in models with vertical space. In mean-field models (based on systems of PDEs or coupled ODEs), zooplankton feeding at a given depth is normally computed as the product of the local functional response and the zooplankton density at this depth. Such simplification is often at odds with field observations which show the absence of clear relationship between intake rates of organisms and the ambient food density. The observed discrepancy is generic and is often caused by fast non-synchronous vertical migration of organisms with different nutrition status. In this paper, we suggest a simple way of incorporating unsynchronized short-term vertical migration of zooplankton into the mean-field modelling framework. We compute grazing of zooplankton in each layer depending on feeding activity of organisms in the layer. We take into account grazing impact of animals which are in the active phase of foraging cycle at the given moment of time but neglect the impact of animals which are in the non-active phase of the cycle (e.g. digesting food). Unsynchronized vertical migration determines the vertical distribution of actively feeding animals in layers depending on vertical distribution of food. In this paper, we compare two generic plankton models: (i) a model based on ‘classical’ grazing approach and (ii) a model incorporating food-mediated unsynchronized vertical migration of zooplankton. We show that including unsynchronized food-mediated migration would make the behaviour of a plankton model more realistic. This would imply a significant enhancement of ecosystem's stability and some additional mechanisms of regulation of algal blooms. In the system with food-mediated unsynchronized vertical migration, the control of phytoplankton by herbivorous becomes possible even for very large concentrations of nutrients in the water (formally, when the system's carrying capacity tends to infinity).     

Using underdominance to bi-stably transform local populations

Underdominance refers to natural selection against individuals with a heterozygous genotype. Here, we analyze a single-locus underdominant system of two large local populations that exchange individuals at a certain migration rate. The system can be characterized by fixed points in the joint allele frequency space. We address the conditions under which underdominance can be applied to transform a local population that is receiving wildtype immigrants from another population. In a single population, underdominance has the benefit of complete removal of genetically modified alleles (reversibility) and coexistence is not stable. The two population system that exchanges migrants can result in internal stable states, where coexistence is maintained, but with additional release of wildtype individuals the system can be reversed to a fully wildtype state. This property is critically controlled by the migration rate. We approximate the critical minimum frequency required to result in a stable population transformation. We also concentrate on the destabilizing effects of fitness and migration rate asymmetry. Practical implications of our results are discussed in the context of utilizing underdominance to genetically modify wild populations. This is of importance especially for genetic pest management strategies, where locally stable and potentially reversible transformations of populations of disease vector species are of interest.     

Population viscosity can promote the evolution of altruistic sterile helpers and eusociality

Because it increases relatedness between interacting individuals, population viscosity has been proposed to favour the evolution of altruistic helping. However, because it increases local competition between relatives, population viscosity may also act as a brake for the evolution of helping behaviours. In simple models, the kin selected fecundity benefits of helping are exactly cancelled out by the cost of increased competition between relatives when helping occurs after dispersal. This result has lead to the widespread view, especially among people working with social organisms, that special conditions are required for the evolution of altruism. Here, we re-examine this result by constructing a simple population genetic model where we analyse whether the evolution of a sterile worker caste (i.e. an extreme case of altruism) can be selected for by limited dispersal. We show that a sterile worker caste can be selected for even under the simplest life-cycle assumptions. This has relevant consequences for our understanding of the evolution of altruism in social organisms, as many social insects are characterized by limited dispersal and significant genetic population structure.     

Predictions of Taylor's power law, density dependence and pink noise from a neutrally modeled time series

There has recently been increasing interest in neutral models of biodiversity and their ability to reproduce the patterns observed in nature, such as species abundance distributions. Here we investigate the ability of a neutral model to predict phenomena observed in single-population time series, a study complementary to most existing work that concentrates on snapshots in time of the whole community. We consider tests for density dependence, the dominant frequencies of population fluctuation (spectral density) and a relationship between the mean and variance of a fluctuating population (Taylor's power law). We simulated an archipelago model of a set of interconnected local communities with variable mortality rate, migration rate, speciation rate, size of local community and number of local communities. Our spectral analysis showed ‘pink noise’: a departure from a standard random walk dynamics in favor of the higher frequency fluctuations which is partly consistent with empirical data. We detected density dependence in local community time series but not in metacommunity time series. The slope of the Taylor's power law in the model was similar to the slopes observed in natural populations, but the fit to the power law was worse. Our observations of pink noise and density dependence can be attributed to the presence of an upper limit to community sizes and to the effect of migration which distorts temporal autocorrelation in local time series. We conclude that some of the phenomena observed in natural time series can emerge from neutral processes, as a result of random zero-sum birth, death and migration. This suggests the neutral model would be a parsimonious null model for future studies of time series data.     

Range expansions of sexual versus asexual organisms: Effects of reproductive assurance and migration load

It is generally considered that sexual organisms show faster evolutionary adaptation than asexual organisms because sexuals can accumulate adaptive mutations through recombination. Yet, empirical evidence often shows that the geographic range size of sexual species is narrower than that of closely related asexual species, which may seem as if asexuals can adapt to more varied environments. Two potential explanations for this apparent contradiction considered by the existing theory are reproduction assurance and migration load. Here, we consider both reproductive assurance and migration load within a single model to comparatively examine their effects on range expansions of sexuals and asexuals across an environmental gradient. The model shows that higher dispersal propensity decreases sexuals' disadvantage in reproductive assurance while increasing their disadvantage in migration load. Moreover, lower mutation rate constrains adaptation more strongly in asexuals than in sexuals. Thus, high dispersal propensity and high mutation rates promote that asexuals have wider range sizes than sexuals. Intriguingly, our model reveals that sexuals can have wider geographic range sizes than asexuals under low dispersal propensity and low mutation rates, a pattern consistent with a few exceptional empirical cases. Combining reproductive assurance and migration load provides a useful perspective to better understand the relationships between species' mating systems and their geographic ranges.     

The effects of the mating system on the evolution of migration in a spatially heterogeneous population

Summary Verbal explanations for the evolution of migration and dispersal often invoke inbreeding depression as an important force. Experimental work on plant populations indicates that while inbreeding depression may favor increased migration rates, adaptation to local environments may reduce the advantage to migrants. We formalize and test this hypothesis using a two-locus genetic model that incorporates lowered fitness in offspring produced by self-fertilization, and habitat differentiation. We also use the model to address questions about the general theory of genetic modifiers and the modifier reduction principle. We find that even under conditions when migration would increase the mean fitness of a population, migration may not be favored. This result is due to the associations that develop between genotypes at a locus subject to overdominant selection and at a neutral locus controlling the migration rate. Thus, it appears that, in this model, the forces of local adaptation, which favor a reduction in the migration rate, overwhelm those of inbreeding depression, which may favor dispersal.     

Development of an in-stream migration model for Gammarus pulex L. (Crustacea, Amphipoda) as a tool in river restoration management

Ecological models can act as interesting tools to support decision-making in river restoration management. In particular models which are capable of predicting the habitat requirements of species are of considerable importance to ensure that the planned actions have the desired effects on the aquatic ecosystem. To this end, Artificial Neural Network (ANN) models were tested and optimized for the prediction of the habitat suitability for Gammarus pulex, a relevant indicator species in water quality assessment. Although ANN models are in general quite robust with a rather high predictive reliability, the model performance had to be increased with regard to simulations for river restoration management. In particular, it has been shown that spatial and temporal expert-rules could possibly be included. Migration dynamics of downstream drift and upstream migration of the organisms and migration barriers along the river (weirs, culverted river sections,␣...) might indeed deliver important additional information on the effectiveness of the restoration plans, and also on the timing of the expected effects. In this context, an additional in-stream migration model for Gammarus pulex was developed. This migration model, implemented in a Geographical Information System (GIS), has been used to simulate a practical river restoration scenario for a river in Flanders, Belgium. The case study illustrated that the removal of a weir, at a particular site, resulted in the improvement of the habitat suitability for Gammarus pulex. The ANN models predicted that after restoration the habitat was suitable again for Gammarus pulex. The migration model indicated that the restored parts of the river would be recolonized within about 2 months. In this way, decision makers can have an idea whether and when a restoration option will have a desired effect.     

Competition and species coexistence in a metapopulation model: Can fast asymmetric migration reverse the outcome of competition in a homogeneous environment?

We investigate whether asymmetric fast migration can modify the predictions of classical competition theory and, in particular revert species dominance. We consider a model of two species competing for an implicit resource on a habitat divided into two patches. Both patches are connected through constant migration rates and in each patch local dynamics are driven by a Lotka-Volterra competition system.Local competition is asymmetric with the same superior competitor in both patches. Migration is asymmetric, species dependent and fast in comparison to local competitive interactions. The species and patches are taken to be otherwise similar: in both patches we assume the same carrying capacities for both species, and the same growth rates and pair-wise competition coefficients for each species.We show that global dynamics can be described by a classical Lotka-Volterra competition model. We found that by modifying the ratio of intraspecific migration rates for both species all possible combinations of global species relative dominance can be achieved. We find specific conditions for which the local superior competitor is globally excluded. This is to our knowledge the first study showing that fast asymmetric migration can lead to inferior competitor dominance in a homogeneous environment. We conclude that disparity of temporal scales between migration and local dynamics may have important consequences for the maintenance of biodiversity in spatially structured populations.     

Altruism in viscous populations — an inclusive fitness model

Summary A viscous population (Hamilton, 1964) is one in which the movement of organisms from their place of birth is relatively slow. This viscosity has two important effects: one is that local interactions tend to be among relatives, and the other is that competition for resources tends to be among relatives. The first effect tends to promote and the second to oppose the evolution of altruistic behaviour. In a simulation model of Wilsonet al. (1992) these two factors appear to exactly balance one another, thus opposing the evolution of local altruistic behaviour. Here I show, with an inclusive fitness model, that the same result holds in a patchstructured population.     

集合群落及其在河流鱼类群聚研究中的应用

集合群落(metacommunity)是指多个潜在相互作用的物种通过它们之间的扩散而连接在一起的一组局域群落,目前已成为斑块生境下生物群落结构、格局和动态的重要理论基础之一。斑块动态、物种排序、群体效应和中性模型等4种理论模型,可用于解释不同情形下集合群落内物种的迁移状况,描述集合群落的动态。可采用群落结构或生态学机制等途径,来阐述所研究的群落是属于哪一种特定的集合群落类型。集合群落可用于研究河流鱼类群聚,解释鱼类的群落结构等问题。另外本文还结合我国水域生态环境及水生生物现状,对今后集合群落的研究作了展望。     

Heterosis increases the effective migration rate

Individuals coming from the same subpopulation are more likely to share deleterious mutations at any given locus than hybrids formed between parents from different populations. Offspring of migrants therefore may experience heterosis and have higher fitness than resident individuals. This will, in turn, result in the immigrant alleles being present in higher frequencies than predicted from neutral expectations and thus a higher effective migration rate. In this paper we derive a formula to calculate the effective migration rate in the presence of heterosis. It is shown that the effect of heterosis on the migration rate can be substantial when fitness reduction within local populations is severe. The effect will be more pronounced in species with relatively short map lengths. Furthermore the heterosis effect will be highly variable throughout the genome, with the largest effect seen near selected genes and in regions of high gene density.     

Asymmetry in species regional dispersal ability and the neutral theory

The neutral assumption that individuals of either the same or different species share exactly the same birth, death, migration, and speciation probabilities is fundamental yet controversial to the neutral theory. Several theoretical studies have demonstrated that a slight difference in species per capita birth or death rates can have a profound consequence on species coexistence and community structure. Whether asymmetry in migration, a vital demographic parameter in the neutral model, plays an important role in community assembly still remains unknown. In this paper, we relaxed the ecological equivalence assumption of the neutral model by introducing differences into species regional dispersal ability. We investigated the effect of asymmetric dispersal on the neutral local community structure. We found that per capita asymmetric dispersal among species could reduce species richness of the local community and result in deviations of species abundance distributions from those predicted by the neutral model. But the effect was moderate compared with that of asymmetries in birth or death rates, unless very large asymmetries in dispersal were assumed. A large difference in species dispersal ability, if there is, can overwhelm the role of random drift and make local community dynamics deterministic. In this case, species with higher regional dispersal abilities tended to dominate in the local community. However, the species abundance distribution of the local community under asymmetric dispersal could be well fitted by the neutral model, but the neutral model generally underestimated the fundamental biodiversity number but overestimated the migration rate in such communities.     

Metapopulation structure of a seed-predator weevil and its host plant in arms race coevolution

Although the importance of gene flow in the geographic structuring of host-parasite interactions has been well discussed, little is known about how dispersal drives the spatial dynamics of other types of coevolutionary interactions in nature. We evaluated the roles of gene flow in the geographically structured processes of a predator-prey arms race involving a seed-predatory weevil with a long mouthpart and its host camellia plant with a thick fruit coat. Molecular genetic analyses showed that both weevil and camellia populations were structured at a spatial scale of several kilometers. Importantly, the spatial pattern of the migration of weevils, but not that of camellias, imposed significant effects on the geographic configuration of the levels of coevolutionary escalation. This result suggests that even if migration is limited in one species (camellia), local coevolution with the other species that migrates between neighboring localities (weevil) can reduce the interpopulation difference in the local adaptive optima of the former species. Thus, gene flow of a species potentially homogenizes the local biological environments provided by the species and thereby promotes the evolutionary convergence of its coevolving counterparts. Consequently, by focusing on coevolutionary interactions in natural communities, "indirect" effects of gene flow on the adaptive divergence of organisms could be identified.     

Guidance of primordial germ cell migration

Primordial germ cells (PGCs), the progenitors of the gametes, migrate from the position where they are specified towards the region where the gonad develops. To reach their target, the PGCs obtain directional cues from cells positioned along their migration path. One such cue, the chemokine SDF-1, has recently been found to be critical for proper PGC migration in zebrafish and in mice. In Drosophila, too, a molecule that is structurally related to chemokine receptors and is important for PGC migration has been identified. The ability to visualize chemokine-guided migration at a high resolution in vivo in these model organisms provides a unique opportunity to study this process, which is relevant for many events in normal development and disease.     

Modeling patterns of zooplankton diel vertical migration

Realized predation pressure, defined as the product of predationpressure and light intensity, expresses the mortality pressuredue to visual predation. The part of realized predation pressurewhich is sensed by organisms is here considered to be relatedto food level and temperature. This partly realized predationpressure is referred to as sensed predation pressure. We proposea possible control mechanism of diel vertical migration (DVM):organisms move vertically following the minimum change in sensedpredation pressure. To investigate this assumption, we presenta math ematical model of DVM. We assume that when predatorsare present, the food level is above a minimal level, and temperatureis higher than the tolerance of organisms to growth, prey organismsundertake DVM following the minimum change in sensed predationpressure. We examine how patterns of migration may be affectedby changes in water clarity, predation pressure, food leveland temperature. This work supports the assumption that minimizingchanges in sensed predation pressure can explain the wide variationin the vertical profile of zooplankton.     

Island models and the coalescent process

Using a coalescent approach, we derive several classical results and extend them to more general models. We find that the classic result for constant population size and constant migration rates holds in models with varying population size and varying migration rates with the obvious substitution of effective population size and mean migration fraction. In addition, the relationship of a 'local' F _ST to local gene flow is derived. This result may be useful for analysing gene flow in a regional subset of a large global population, using only data from the regional subset.