Detecting ancient admixture in humans using sequence polymorphism data

A debate of long-standing interest in human evolution centers around whether archaic human populations (such as the Neanderthals) have contributed to the modern gene pool. A model of ancient population structure with recent mixing is introduced, and it is determined how much information (i.e., sequence data from how many unlinked nuclear loci) would be necessary to distinguish between different demographic scenarios. It is found that approximately 50-100 loci are necessary if plausible parameter estimates are used. There are not enough data available at the present to support either the "single origin" or the "multiregional" model of modern human evolution. However, this information should be available in a few years.     

Limits to local adaptation in six populations of the annual plant Diodia teres

* Local adaptation is common, but tests for adaptive differentiation frequently compare populations from strongly divergent environments, making it unlikely that any influence of stochastic processes such as drift or mutation on local adaptation will be detected. Here, the hypothesis that local adaptation is more likely to develop when the native environments of populations are more distinct than when they are similar was tested. * A reciprocal transplant experiment including two populations from each of three habitats was conducted to determine the pattern of local adaptation. In addition to testing for local adaptation at the population level, the hypothesis was tested that local adaptation is more common between populations from different habitats than between populations from the same habitat. * Local adaptation was not common, but more evidence was found of local adaptation between populations from different habitats than between populations from the same habitat. Two instances of foreign genotype fitness advantage confirm that stochastic processes such as drift can limit local adaptation. * These results are consistent with the hypothesis that stochastic processes can inhibit local adaptation but are more likely to be overwhelmed by natural selection when populations occur in divergent environments.     

Evolution of dispersal polymorphism and local adaptation of dispersal distance in spatially structured landscapes

Many organisms show polymorphism in dispersal distance strategies. This variation is particularly ecological relevant if it encompasses a functional separation of short‐ (SDD) and long‐distance dispersal (LDD). It remains, however, an open question whether both parts of the dispersal kernel are similarly affected by landscape related selection pressures. We implemented an individual‐based model to analyze the evolution of dispersal traits in fractal landscapes that vary in the proportion of habitat and its spatial configuration. Individuals are parthenogenetic with dispersal distance determined by two alleles on each individual's genome: one allele coding for the probability of global dispersal and one allele coding for the variance σ of a Gaussian local dispersal with mean value zero. Simulations show that mean distances of local dispersal and the probability of global dispersal, increase with increasing habitat availability, but that changes in the habitat's spatial autocorrelation impose opposing selective pressure: local dispersal distances decrease and global dispersal probabilities increase with decreasing spatial autocorrelation of the available habitat. Local adaptation of local dispersal distance emerges in landscapes with less than 70% of clumped habitat. These results demonstrate that long and short distance dispersal evolve separately according to different properties of the landscape. The landscape structure may consequently largely affect the evolution of dispersal distance strategies and the level of dispersal polymorphism.     

Estimating the size of closed populations using inverse multiple-recapture sampling

A log-linear model for estimating the size of a closed population is defined for inverse multiple-recapture sampling with dependent samples. Efficient estimators of the log-linear model parameters and the population size are obtained by the method of minimum chi-square. A chi-square test of the general linear hypothesis regarding the log-linear model parameters is defined.     

Detecting chaotic dynamics of insect populations from long-term survey data

1. Estimates of the Lyapunov exponent, a statistic that measures the sensitive dependence of the dynamic behaviour of a system on its initial conditions, are used to characterize several sets of insect time series.
2. A new method is described to overcome the difficulty of defining the dynamics of an observed, noisy, short ecological time series. This method provides two test statistics for the estimated Lyapunov exponent.
3. This method is applied to forty-six time series comprising six aphid species from five sites and four moth species from six sites. There are few positive Lyapunov exponents and none is sufficiently large to characterize its time series as chaotic.
4. Two methods to estimate the Lyapunov exponent are compared; that based on logarithmically transformed counts yields less variable estimates for highly variable insect data than that based on untransformed counts.     

Major histocompatibility complex polymorphism: dynamics and consequences of parasite-mediated local adaptation in fishes

Parasitism is a common form of life and represents a strong selective pressure for host organisms. In response to this evolutionary pressure, vertebrates have developed genetically coded defences such as the major histocompatibility complex (MHC). Mechanisms of parasite-mediated selection not only maintain outstanding polymorphism in these genes but have also been proposed to further promote host population divergence and ultimately speciation because it can drive evolution of local adaptation in which MHC genes play a crucial role. This review first highlights the dynamics and complexity of parasite-mediated selection in natural systems, which not only depends on dominating parasite strategies and on the taxonomic diversity of the parasite community but also includes the differences in parasite communities between habitats and niches, creating divergent selection on locally adapted populations. Then the different ways in which MHC genes potentially allow vertebrates to respond to these dynamics and to adapt locally are outlined. Finally, it is proposed that varying selection strength in time and space may lead to variation in the strength of precopulatory reproductive isolation which has evolved to maintain local adaptation.     

Effective and efficient data sampling using bitmap indices

With growing computational capabilities of parallel machines, scientific simulations are being performed at finer spatial and temporal scales, leading to a data explosion. The growing sizes are making it extremely hard to store, manage, disseminate, analyze, and visualize these datasets, especially as neither the memory capacity of parallel machines, memory access speeds, nor disk bandwidths are increasing at the same rate as the computing power. Sampling can be an effective technique to address the above challenges, but it is extremely important to ensure that dataset characteristics are preserved, and the loss of accuracy is within acceptable levels. In this paper, we address the data explosion problems by developing a novel sampling approach, and implementing it in a flexible system that supports server-side sampling and data subsetting. We observe that to allow subsetting over scientific datasets, data repositories are likely to use an indexing technique. Among these techniques, we see that bitmap indexing can not only effectively support subsetting over scientific datasets, but can also help create samples that preserve both value and spatial distributions over scientific datasets. We have developed algorithms for using bitmap indices to sample datasets. We have also shown how only a small amount of additional metadata stored with bitvectors can help assess loss of accuracy with a particular subsampling level. Some of the other properties of this novel approach include: (1) sampling can be flexibly applied to a subset of the original dataset, which may be specified using a value-based and/or a dimension-based subsetting predicate, and (2) no data reorganization is needed, once bitmap indices have been generated. We have extensively evaluated our method with different types of datasets and applications, and demonstrated the effectiveness of our approach.     

Evidence for local adaptation in neighbouring Daphnia populations: a laboratory transplant experiment

1. We tested whether two neighbouring Daphnia galeata populations (from Lake Blankaart and Fish Pond), separated only by 5 m of land and with the occasional exchange of water were genetically differentiated in allozyme markers and life history traits. Allozyme electrophoresis revealed that the populations differed in allelic as well as in genotypic composition. 2. In a laboratory transplant experiment, in which animals of four clones of each of the populations were raised in the water of both ponds, survival in Blankaart water was high for both the Blankaart and Fish Pond clones, whereas survival in Fish Pond water was high for the Fish Pond clones, but low for the Blankaart clones. 3. Fish Pond clones produced fewer neonates than Blankaart clones when cultured in Blankaart water. High egg mortality was observed for animals that were raised in Blankaart water, and this egg mortality was higher for Fish Pond clones than for Blankaart clones. 4. Our results provide evidence for genetic differentiation between Daphnia populations inhabiting neighbouring water bodies and suggest local adaptation to environmental conditions other than direct predation.     

Parameters adaptation in the populations models

Ecology-evolutionary models of low dimensions were developed on the basis of competitive selection criteria. Dynamics of variables (number of individuals) and the search of evolutionary-stable values of parameters (biological characterictics of populations) were monitored in the suggested models. If the environmental temperature is changing periodically, the average (a) and width (d) of temperature tolerance range appears to be the important parameters. By model experiments it was established that stable values of temperature (a), favorable for development of highly specialized algae (d is low) were close to minimum and maximum of temperature curve. And for the low specialized algae (d is high) this values were close to the average temperature of environment. In a similar manner, a set of evolutionally stable parameters (a, d) was established for either of the two interacted populations (competitors and "predator-prey"). The hypotheses concerning it's geometric structure and the process of coevolution is formulated.     

Unusual mitochondrial DNA polymorphism in two local populations of blue tit Parus caeruleus

Mitochondrial DNA (mtDNA) from 25 blue tits Parus caeruleus sampled from two populations of the Grenoble region (France) was assayed for polymorphism with 17 restriction endonucleases. Nine genotypes were found. Several mtDNA genotypes were also analysed by amplification via the polymerase chain reaction (PCR) and direct sequencing of 903 bp of the cytochrome b gene. The mtDNA polymorphism is greater in P. caeruleus than in other comparable bird species and results from the presence of two clearly differentiated mitochondrial lineages. Using the data of restriction polymorphism, the mean sequence divergence between individuals of the two lineages is 1.23%. Therefore, P. caeruleus should fall into the category II of phylogeographic pattern sensu Avise et al. (1987): discontinuous mtDNA genotypes which co-occur in the same region. P. caeruleus, like humans and other mobile species with high gene flow, seems to have lost its geographic structure in terms of mtDNA phylogeny. This unusual mitochondrial polymorphism can be explained by the recent admixture of two long-term isolated populations. This could be accounted for by two different scenarios. One assumes a simultaneous post-glacial colonization of the Grenoble region by two isolated European populations of P. caeruleus. Alternatively, hybridization between P. caeruleus and P. cyanus could have caused the observed pattern of mtDNA variation.     

Testing the role of phenotypic plasticity for local adaptation: growth and development in time-constrained Rana temporaria populations

Phenotypic plasticity can be important for local adaptation, because it enables individuals to survive in a novel environment until genetic changes have been accumulated by genetic accommodation. By analysing the relationship between development rate and growth rate, it can be determined whether plasticity in life-history traits is caused by changed physiology or behaviour. We extended this to examine whether plasticity had been aiding local adaptation, by investigating whether the plastic response had been fixed in locally adapted populations. Tadpoles from island populations of Rana temporaria, locally adapted to different pool-drying regimes, were monitored in a common garden. Individual differences in development rate were caused by different foraging efficiency. However, developmental plasticity was physiologically mediated by trading off growth against development rate. Surprisingly, plasticity has not aided local adaptation to time-stressed environments, because local adaptation was not caused by genetic assimilation but on selection on the standing genetic variation in development time.     

Inference of parasite local adaptation using two different fitness components

Estimating parasite fitness is central to studies aiming to understand parasite evolution. Theoretical models generally use the basic reproductive rate R(0) to express fitness, yet it is very difficult to quantify R(0) empirically and experimental studies often use fitness components such as infection intensity or infectivity as substitutes. These surrogate measures may be biased in several ways. We assessed local adaptation of the microsporidium Ordospora colligata to its host, the crustacean Daphnia magna using two different parasite fitness components: infection persistence over several host generations in experimental populations and infection intensity in individual hosts. We argue that infection persistence is a close estimator of R(0), whereas infection intensity measures only a component of it. Both measures show a pattern that is consistent with parasite local adaptation and they correlate positively. However, several inconsistencies between them suggest that infection intensity may at times provide an inadequate estimate of parasite fitness.     

Detecting bottlenecks using BOTTLENECK 1.2.02 in wild populations: the importance of the microsatellite structure

Reduced, or bottlenecked, populations are more prone to adverse events. Thus, the detection of genetic bottleneck signatures in wildlife is an important issue for conservation. BOTTLENECK 1.2.02 is a software commonly used for detecting genetic characteristics of past bottlenecks. Here we test the efficiency with which this software detects bottlenecks in two koala populations of known history. The sign test performed well for both populations, particularly under the infinite alleles model for mutation. This suggests this model could be the more realistic for marsupial microsatellites than other mutation models. Under the allele frequency distribution test, the two populations falsely appeared to be at mutation/drift equilibrium. However, this test could detect the bottleneck when only imperfect repeat microsatellites were included in the analysis. We thus recommend further investigation of imperfect repeat microsatellites, which could be more powerful for bottleneck detection. These results underline the cautious approach researchers and conservationists should take when studying the past of unknown populations.     

Inferring admixture proportions from molecular data: extension to any number of parental populations

The relative contribution of two parental populations to a hybrid group (the admixture proportions) can be estimated using not only the frequencies of different alleles, but also the degree of molecular divergence between them. In this paper, we extend this possibility to the case of any number of parental populations. The newly derived multiparental estimator is tested by Monte Carlo simulations and by generating artificial hybrid groups by pooling mtDNA samples from human populations. The general properties (including the variance) of the two-parental estimator seem to be retained by the multiparental estimator. When mixed human populations are considered and hypervariable single-locus data are analyzed (mtDNA control region), errors in the estimated contributions appear reasonably low only when highly differentiated parental populations are involved. Finally, the method applied to the hybrid Canary Island population points to a much lower female contribution from Spain than has previously been estimated.     

Detecting local adaptation in a natural plant-pathogen metapopulation: a laboratory vs. field transplant approach

Antagonistic coevolution between hosts and parasites in spatially structured populations can result in local adaptation of parasites. Traditionally parasite local adaptation has been investigated in field transplant experiments or in the laboratory under a constant environment. Despite the conceptual importance of local adaptation in studies of (co)evolution, to date no study has provided a comparative analysis of these two methods. Here, using information on pathogen population dynamics, I tested local adaptation of the specialist phytopathogen, Podosphaera plantaginis, to its host, Plantago lanceolata at three different spatial scales: sympatric host population, sympatric host metapopulation and allopatric host metapopulations. The experiment was carried out as a field transplant experiment with greenhouse-reared host plants from these three different origins introduced into four pathogen populations. In contrast to results of an earlier study performed with these same host and parasite populations under laboratory conditions, I did not find any evidence for parasite local adaptation. For interactions governed by strain-specific resistance, field studies may not be sensitive enough to detect mean parasite population virulence. Given that parasite transmission potential may be mediated by the abiotic environment and genotype-by-environment interactions, I suggest that relevant environmental variation should be incorporated into laboratory studies of parasite local adaptation.     

Prediction of hybrid performance in maize using molecular markers and joint analyses of hybrids and parental inbreds

The identification of superior hybrids is important for the success of a hybrid breeding program. However, field evaluation of all possible crosses among inbred lines requires extremely large resources. Therefore, efforts have been made to predict hybrid performance (HP) by using field data of related genotypes and molecular markers. In the present study, the main objective was to assess the usefulness of pedigree information in combination with the covariance between general combining ability (GCA) and per se performance of parental lines for HP prediction. In addition, we compared the prediction efficiency of AFLP and SSR marker data, estimated marker effects separately for reciprocal allelic configurations (among heterotic groups) of heterozygous marker loci in hybrids, and imputed missing AFLP marker data for marker-based HP prediction. Unbalanced field data of 400 maize dent × flint hybrids from 9 factorials and of 79 inbred parents were subjected to joint analyses with mixed linear models. The inbreds were genotyped with 910 AFLP and 256 SSR markers. Efficiency of prediction (R ²) was estimated by cross-validation for hybrids having no or one parent evaluated in testcrosses. Best linear unbiased prediction of GCA and specific combining ability resulted in the highest efficiencies for HP prediction for both traits (R ² = 0.6–0.9), if pedigree and line per se data were used. However, without such data, HP for grain yield was more efficiently predicted using molecular markers. The additional modifications of the marker-based approaches had no clear effect. Our study showed the high potential of joint analyses of hybrids and parental inbred lines for the prediction of performance of untested hybrids.     

Mathematical properties of Fst between admixed populations and their parental source populations

We consider the properties of the F(st) measure of genetic divergence between an admixed population and its parental source populations. Among all possible populations admixed among an arbitrary set of parental populations, we show that the value of F(st) between an admixed population and a specific source population is maximized when the admixed population is simply the most distant of the other source populations. For the case with only two parental populations, as a function of the admixture fraction, we further demonstrate that this F(st) value is monotonic and convex, so that F(st) is informative about the admixture fraction. We illustrate our results using example human population-genetic data, showing how they provide a framework in which to interpret the features of F(st) in admixed populations.