ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Reciprocal cooperation between unrelated rats depends on cost to donor and benefit to recipient

Background Although evolutionary models of cooperation build on the intuition that costs of the donor and benefits to the receiver are the most general fundamental parameters, it is largely unknown how they affect the decision of animals to cooperate with an unrelated social partner. Here we test experimentally whether costs to the donor and need of the receiver decide about the amount of help provided by unrelated rats in an iterated prisoner's dilemma game. Results Fourteen unrelated Norway rats were alternately presented to a cooperative or defective partner for whom they could provide food via a mechanical apparatus. Direct costs for this task and the need of the receiver were manipulated in two separate experiments. Rats provided more food to cooperative partners than to defectors (direct reciprocity). The propensity to discriminate between helpful and non-helpful social partners was contingent on costs: An experimentally increased resistance in one Newton steps to pull food for the social partner reduced the help provided to defectors more strongly than the help returned to cooperators. Furthermore, test rats provided more help to hungry receivers that were light or in poor condition, which might suggest empathy, whereas this relationship was inverse when experimental partners were satiated. Conclusions In a prisoner's dilemma situation rats seem to take effect of own costs and potential benefits to a receiver when deciding about helping a social partner, which confirms the predictions of reciprocal cooperation. Thus, factors that had been believed to be largely confined to human social behaviour apparently influence the behaviour of other social animals as well, despite widespread scepticism. Therefore our results shed new light on the biological basis of reciprocity.

Explosive speciation rates and unusual species richness in haplochromine cichlid fishes: Effects of sexual selection

Ancient lakes are often unusually species rich, mostly as a result of radiation and species-flock formation having taken place in only one or a few of many taxa present. Understanding why some taxa radiate and others do not is at the heart of understanding biodiversity. In this chapter I discuss possible explanations for disproportionally large species numbers in some cichlid fish lineages in East African Great Lakes: the halochromine cichlid fishes in Lakes Victoria and Malawi. I show that speciation rates in this group are higher than in any other lacustrine fish radiation. Against this background, I review hypotheses put forward to explain diversity in cichlid species flocks. The evolution of species diversity requires three processes: speciation, ecological radiation and anatomical diversification, and it is wrong to consider hypotheses that are relevant to different processes as alternatives to each other. The African cichlid species flocks show unusually high ecological species packing in several phylogenetic groups and unusually high speciation rates in haplochromines. Therefore, it maybe concluded that at least two evolutionary models are required to explain the difference between cichlid diversity and other fish diversity in East African Lakes: one for speciation in haplochromines and one for coexistence. Subsequently I review work on speciation in haplochromines, and in particular studies aimed at testing the hypothesis of speciation by sexual selection. Haplochromines have a polygynous mating system, conducive to sexual selection, but other polygynous cichlids are not particularly species rich. This suggests that more than just strong sexual selection is required to explain haplochromine species richness. Recent palaeoecological evidence undermines the previously popular hypotheses that explained the species richness of Lake Victoria in terms of speciation under varying natural or sexual selection regimes in satellite lakes or in isolated lake basins. I summarize experimental and comparative studies, which provide evidence for two mechanisms of sympatric speciation by disruptive sexual selection on polymorphic coloration. Such modes of speciation may explain (i) the high speciation rates in colour polymorphic lineages of haplochromine cichlids under conditions where colour variation is visible in clear water, and (ii) in combination with factors that affect population survival, the unusual species richness in haplochromine species flocks. I argue that sexual selection, if disruptive, can accelerate the pace of adaptive radiation because the resultant genetic population fragmentation allows a much increased rate of differential response to disruptive natural selection. Hence, the ecological pattern of diversity resembles that produced by disruptive natural selection, with the difference that disruptive sexual selection continues to cause (gross) speciation even after niche space is saturated. This may explain the unusually high numbers of very closely related and ecologically similar species in haplochromine species flocks. The role of disruptive sexual selection is twofold: it not only causes speciation, but also maintains reproductive isolation in sympatry between species that have evolved in sympatry or allopatry. Therefore, the maintenance of diversity in species flocks that originated through sexual selection depends on the persistence of the selection regime within the environmental signal space under which that diversity evolved.

Population genomic tests of models of adaptive radiation in Lake Victoria Region cichlid fish

Adaptive radiation is usually thought to be associated with speciation, but the evolution of intraspecific polymorphisms without speciation is also possible. The radiation of cichlid fish in Lake Victoria (LV) is perhaps the most impressive example of a recent rapid adaptive radiation, with 600+ very young species. Key questions about its origin remain poorly characterized, such as the importance of speciation versus polymorphism, whether species persist on evolutionary time scales, and if speciation happens more commonly in small isolated or in large connected populations. We used 320 individuals from 105 putative species from Lakes Victoria, Edward, Kivu, Albert, Nabugabo and Saka, in a radiation-wide amplified fragment length polymorphism (AFLP) genome scan to address some of these questions. We demonstrate pervasive signatures of speciation supporting the classical model of adaptive radiation associated with speciation. A positive relationship between the age of lakes and the average genomic differentiation of their species, and a significant fraction of molecular variance explained by above-species level taxonomy suggest the persistence of species on evolutionary time scales, with radiation through sequential speciation rather than a single starburst. Finally the large gene diversity retained from colonization to individual species in every radiation suggests large effective population sizes and makes speciation in small geographical isolates unlikely.

Evolutionary demography of long-lived monocarpic perennials: a time-lagged integral projection model

1. The evolution of flowering strategies (when and at what size to flower) in monocarpic perennials is determined by balancing current reproduction with expected future reproduction, and these are largely determined by size-specific patterns of growth and survival. However, because of the difficulty in following long-lived individuals throughout their lives, this theory has largely been tested using short-lived species (< 5 years). 2. Here, we tested this theory using the long-lived monocarpic perennial Campanula thyrsoides which can live up to 16 years. We used a novel approach that combined permanent plot and herb chronology data from a 3-year field study to parameterize and validate integral projection models (IPMs). 3. Similar to other monocarpic species, the rosette leaves of C. thyrsoides wither over winter and so size cannot be measured in the year of flowering. We therefore extended the existing IPM framework to incorporate an additional time delay that arises because flowering demography must be predicted from rosette size in the year before flowering. 4. We found that all main demographic functions (growth, survival probability, flowering probability and fecundity) were strongly size-dependent and there was a pronounced threshold size of flowering. There was good agreement between the predicted distribution of flowering ages obtained from the IPMs and that estimated in the field. Mostly, there was good agreement between the IPM predictions and the direct quantitative field measurements regarding the demographic parameters lambda, R-0 and T. We therefore conclude that the model captures the main demographic features of the field populations. 5. Elasticity analysis indicated that changes in the survival and growth function had the largest effect (c. 80%) on lambda and this was considerably larger than in short-lived monocarps. We found only weak selection pressure operating on the observed flowering strategy which was close to the predicted evolutionary stable strategy. 6. Synthesis. The extended IPM accurately described the demography of a long-lived monocarpic perennial using data collected over a relatively short period. We could show that the evolution of flowering strategies in short- and long-lived monocarps seem to follow the same general rules but with a longevity-related emphasis on survival over fecundity.

Nuclear and plastid markers reveal the persistence of genetic identity: A new perspective on the evolutionary history of Petunia exserta

Recently divergent species that can hybridize are ideal models for investigating the genetic exchanges that can occur while preserving the species boundaries. Petunia exserta is an endemic species from a very limited and specific area that grows exclusively in rocky shelters. These shaded spots are an inhospitable habitat for all other Petunia species, including the closely related and widely distributed species P. axillaris. Individuals with intermediate morphologic characteristics have been found near the rocky shelters and were believed to be putative hybrids between P. exserta and P. axillaris, suggesting a situation where Petunia exserta is losing its genetic identity. In the current study, we analyzed the plastid intergenic spacers trnS/trnG and trnH/psbA and six nuclear CAPS markers in a large sampling design of both species to understand the evolutionary process occurring in this biological system. Bayesian clustering methods, cpDNA haplotype networks, genetic diversity statistics, and coalescence-based analyses support a scenario where hybridization occurs while two genetic clusters corresponding to two species are maintained. Our results reinforce the importance of coupling differentially inherited markers with an extensive geographic sample to assess the evolutionary dynamics of recently diverged species that can hybridize. (C) 2013 Elsevier Inc. All rights reserved.