A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.

Inferring population history with DIY ABC: a user-friendly approach to approximate Bayesian computation

Genetic data obtained on population samples convey information about their evolutionary history. Inference methods can extract part of this information but they require sophisticated statistical techniques that have been made available to the biologist community (through computer programs) only for simple and standard situations typically involving a small number of samples. We propose here a computer program (DIY ABC) for inference based on approximate Bayesian computation (ABC), in which scenarios can be customized by the user to fit many complex situations involving any number of populations and samples. Such scenarios involve any combination of population divergences, admixtures and population size changes. DIY ABC can be used to compare competing scenarios, estimate parameters for one or more scenarios and compute bias and precision measures for a given scenario and known values of parameters (the current version applies to unlinked microsatellite data). This article describes key methods used in the program and provides its main features. The analysis of one simulated and one real dataset, both with complex evolutionary scenarios, illustrates the main possibilities of DIY ABC.

An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure

There is great interest in using amplified fragment length polymorphism (AFLP) markers because they are inexpensive and easy to produce. It is, therefore, possible to generate a large number of markers that have a wide coverage of species genotnes. Several statistical methods have been proposed to study the genetic structure using AFLP's but they assume Hardy-Weinberg equilibrium and do not estimate the inbreeding coefficient, F-IS. A Bayesian method has been proposed by Holsinger and colleagues that relaxes these simplifying assumptions but we have identified two sources of bias that can influence estimates based on these markers: (i) the use of a uniform prior on ancestral allele frequencies and (ii) the ascertainment bias of AFLP markers. We present a new Bayesian method that avoids these biases by using an implementation based on the approximate Bayesian computation (ABC) algorithm. This new method estimates population-specific F-IS and F-ST values and offers users the possibility of taking into account the criteria for selecting the markers that are used in the analyses. The software is available at our web site (http://www-leca.uif-grenoble.fi-/logiciels.htm). Finally, we provide advice on how to avoid the effects of ascertainment bias.

Extending Zelterman’s approach for robust estimation of population size to zero-truncated clustered data

Estimation of population size with missing zero-class is an important problem that is encountered in epidemiological assessment studies. Fitting a Poisson model to the observed data by the method of maximum likelihood and estimation of the population size based on this fit is an approach that has been widely used for this purpose. In practice, however, the Poisson assumption is seldom satisfied. Zelterman (1988) has proposed a robust estimator for unclustered data that works well in a wide class of distributions applicable for count data. In the work presented here, we extend this estimator to clustered data. The estimator requires fitting a zero-truncated homogeneous Poisson model by maximum likelihood and thereby using a Horvitz-Thompson estimator of population size. This was found to work well, when the data follow the hypothesized homogeneous Poisson model. However, when the true distribution deviates from the hypothesized model, the population size was found to be underestimated. In the search of a more robust estimator, we focused on three models that use all clusters with exactly one case, those clusters with exactly two cases and those with exactly three cases to estimate the probability of the zero-class and thereby use data collected on all the clusters in the Horvitz-Thompson estimator of population size. Loss in efficiency associated with gain in robustness was examined based on a simulation study. As a trade-off between gain in robustness and loss in efficiency, the model that uses data collected on clusters with at most three cases to estimate the probability of the zero-class was found to be preferred in general. In applications, we recommend obtaining estimates from all three models and making a choice considering the estimates from the three models, robustness and the loss in efficiency. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Farm income, population and farmland prices: a relative information approach

This paper uses an entropy-based information approach to determine if farmland values are more closely associated with urban pressure or farm income. The basic question is: how much information on changes in farm real estate values is contained in changes in population versus changes in returns to production agriculture? Results suggest population is informative, but changes in farmland values are more strongly associated with changes in the distribution of returns. However, this relationship is not true for every region nor does it hold over time, as for some regions and time periods changes in population are more informative. Results have policy implications for both equity and efficiency.

Pollinator community responses to the spatial population structure of wild plants: A pan-European approach

Land-use changes can alter the spatial population structure of plant species, which may in turn affect the attractiveness of flower aggregations to different groups of pollinators at different spatial scales. To assess how pollinators respond to spatial heterogeneity of plant distributions and whether honeybees affect visitation by other pollinators we used an extensive data set comprising ten plant species and their flower visitors from five European countries. In particular we tested the hypothesis that the composition of the flower visitor community in terms of visitation frequencies by different pollinator groups were affected by the spatial plant population structure, viz. area and density measures, at a within-population (‘patch’) and among-population (‘population’) scale. We found that patch area and population density were the spatial variables that best explained the variation in visitation frequencies within the pollinator community. Honeybees had higher visitation frequencies in larger patches, while bumblebees and hoverflies had higher visitation frequencies in sparser populations. Solitary bees had higher visitation frequencies in sparser populations and smaller patches. We also tested the hypothesis that honeybees affect the composition of the pollinator community by altering the visitation frequencies of other groups of pollinators. There was a positive relationship between visitation frequencies of honeybees and bumblebees, while the relationship with hoverflies and solitary bees varied (positive, negative and no relationship) depending on the plant species under study. The overall conclusion is that the spatial structure of plant populations affects different groups of pollinators in contrasting ways at both the local (‘patch’) and the larger (‘population’) scales and, that honeybees affect the flower visitation by other pollinator groups in various ways, depending on the plant species under study. These contrasting responses emphasize the need to investigate the entire pollinator community when the effects of landscape change on plant–pollinator interactions are studied.

Palaeobiology of the medieval population of Albano (Rome, Italy): a combined morphological and biomolecular approach

Recent years have seen increased interest in skeletal populations from the Imperial Roman Age in Italy, but much less is known about diet and standards of living in the subsequent medieval period. To fill this gap, we conducted a morphological analysis of human remains from Albano, an Italian town near Rome, as well as a stable isotope analysis of bone collagen to reconstruct diet. The sample was recovered from a Medieval cemetery (1040–1220 cal. yr. BP) located in the gardens of the historical Palazzo Doria Pamphili in Albano. A minimum number of 40 individuals (31 adults and 9 sub-adults) were examined using standard methods. Though the general health status of the population was good, signs of cribra orbitalia and diffuse enthesopathies were noted during the morphological examination. Stable carbon and nitrogen isotope analyses of the bone collagen from 24 adult humans and three faunal bones indicate that the diet of the population may be described as predominantly terrestrial and C3-plant based although the data for some of the individuals suggest a modest consumption of C4-(millet) based or aquatic proteins. No evidence of significant dietary differences between the sexes was found. The comparison of the isotope data from Albano with those from populations recovered in the same region is consistent with a shift from a terrestrial, possibly plant foods-dominated subsistence in the Early Middle Ages to a diet with a higher contribution from animal proteins, both terrestrial and aquatic, in the Later Middle Ages.

Migration and diversity in Roman Britain: a multidisciplinary approach to immigrants in Roman York, England

Previous anthropological investigations at Trentholme Drive, in Roman York identified an unusual amount of cranial variation amongst the inhabitants, with some individuals suggested as having originated from the Middle East or North Africa. The current study investigates the validity of this assessment using modern anthropological methods to assess cranial variation in two groups: The Railway and Trentholme Drive. Strontium and oxygen isotope evidence derived from the dentition of 43 of these individuals was combined with the craniometric data to provide information on possible levels of migration and the range of homelands that may be represented. The results of the craniometric analysis indicated that the majority of the York population had European origins, but that 11% of the Trentholme Drive and 12% of The Railway study samples were likely of African decent. Oxygen analysis identified four incomers, three from areas warmer than the UK and one from a cooler or more continental climate. Although based on a relatively small sample of the overall population at York, this multidisciplinary approach made it possible to identify incomers, both men and women, from across the Empire. Evidence for possible second generation migrants was also suggested. The results confirm the presence of a heterogeneous population resident in York and highlight the diversity, rather than the uniformity, of the population in Roman Britain. Am J Phys Anthropol 140:546-561, 2009. (C) 2009 Wiley-Liss, Inc

Effects of patch size and density on flower visitation and seed set of wild plants: a pan-European approach

1. Habitat fragmentation can affect pollinator and plant population structure in terms of species composition, abundance, area covered and density of flowering plants. This, in turn, may affect pollinator visitation frequency, pollen deposition, seed set and plant fitness. 2. A reduction in the quantity of flower visits can be coupled with a reduction in the quality of pollination service and hence the plants’ overall reproductive success and long-term survival. Understanding the relationship between plant population size and⁄ or isolation and pollination limitation is of fundamental importance for plant conservation. 3. Weexamined flower visitation and seed set of 10 different plant species fromfive European countries to investigate the general effects of plant populations size and density, both within (patch level) and between populations (population level), on seed set and pollination limitation. 4. Wefound evidence that the effects of area and density of flowering plant assemblages were generally more pronounced at the patch level than at the population level. We also found that patch and population level together influenced flower visitation and seed set, and the latter increased with increasing patch area and density, but this effect was only apparent in small populations. 5. Synthesis. By using an extensive pan-European data set on flower visitation and seed set we have identified a general pattern in the interplay between the attractiveness of flowering plant patches for pollinators and density dependence of flower visitation, and also a strong plant species-specific response to habitat fragmentation effects. This can guide efforts to conserve plant–pollinator interactions, ecosystem functioning and plant fitness in fragmented habitats.

Comparison of survival rates of captive-reared and wild-bred Mauritius kestrels (Falco punctatus) in a re-introduced population

Re-introduction is a technique widely used in the conservation of threatened bird species. With advances in aviculture the use of captive-produced individuals as the release stock is becoming more commonplace, and ideally, survival of captive-produced, released individuals should be no different from their wild-bred counterparts. During the late 1980s the Critically Endangered Mauritius kestrel (Falco punctatus) was successfully re-introduced into the Bambous mountain range, Mauritius, some 30 years after its local extinction. Between 1987 and 2001 the developing population was closely monitored enabling us to construct re-sighting histories for 88 released and 284 wild-bred kestrels. We used age-structured models in the survival analysis software program MARK to determine if an individual's origin influenced its subsequent survival. Our analysis indicated no compelling evidence for reduced survival among juvenile captive-reared and released individuals, relative to their wild-bred counterparts, across the majority of cohorts and only limited evidence of a cohort-specific effect. This study illustrates that despite the lack of a formal experimental approach it is still feasible to conduct an assessment of re-introduction outcomes and techniques.

Quantification of ruminal Clostridium proteoclasticum by real-time PCR using a molecular beacon approach

Aims: All members of the ruminal Butyrivibrio group convert linoleic acid (cis-9,cis-12-18 : 2) via conjugated 18 : 2 metabolites (mainly cis-9,trans-11-18 : 2, conjugated linoleic acid) to vaccenic acid (trans-11-18 : 1), but only members of a small branch, which includes Clostridium proteoclasticum, of this heterogeneous group further reduce vaccenic acid to stearic acid (18 : 0, SA). The aims of this study were to develop a real-time polymerase chain reaction (PCR) assay that would detect and quantify these key SA producers and to use this method to detect diet-associated changes in their populations in ruminal digesta of lactating cows. Materials and Results: The use of primers targeting the 16S rRNA gene of Cl. proteoclasticum was not sufficiently specific when only binding dyes were used for detection in real-time PCR. Their sequences were too similar to some nonproducing strains. A molecular beacon probe was designed specifically to detect and quantify the 16S rRNA genes of the Cl. proteoclasticum subgroup. The probe was characterized by its melting curve and validated using five SA-producing and ten nonproducing Butyrivibrio-like strains and 13 other common ruminal bacteria. Analysis of ruminal digesta collected from dairy cows fed different proportions of starch and fibre indicated a Cl. proteoclasticum population of 2-9% of the eubacterial community. The influence of diet on numbers of these bacteria was less than variations between individual cows. Conclusion: A molecular beacon approach in qPCR enables the detection of Cl. proteoclasticum in ruminal digesta. Their numbers are highly variable between individual animals. Signifance and Impact of the Study: SA producers are fundamental to the flow of polyunsaturated fatty acid and vaccenic acid from the rumen. The method described here enabled preliminary information to be obtained about the size of this population. Further application of the method to digesta samples from cows fed diets of more variable composition should enable us to understand how to control these bacteria in order to enhance the nutritional characteristics of ruminant-derived foods, including milk and beef.

A novel transcriptomic approach to identify candidate genes for grain quality traits in wheat

A novel methodology is described in which transcriptomics is combined with the measurement of bread-making quality and other agronomic traits for wheat genotypes grown in different environments (wet and cool or hot and dry conditions) to identify transcripts associated with these traits. Seven doubled haploid lines from the Spark x Rialto mapping population were selected to be matched for development and known alleles affecting quality. These were grown in polytunnels with different environments applied 14 days post-anthesis, and the whole experiment was repeated over 2 years. Transcriptomics using the wheat Affymetrix chip was carried out on whole caryopsis samples at two stages during grain filling. Transcript abundance was correlated with the traits for approximately 400 transcripts. About 30 of these were selected as being of most interest, and markers were derived from them and mapped using the population. Expression was identified as being under cis control for 11 of these and under trans control for 18. These transcripts are candidates for involvement in the biological processes which underlie genotypic variation in these traits.

Population genetic models can be used to study the evolution of the interacting behaviours of parents and their progeny

An example of the evolution of the interacting behaviours of parents and progeny is studied using iterative equations linking the frequencies of the gametes produced by the progeny to the frequencies of the gametes in the parental generation. This population genetics approach shows that a model in which both behaviours are determined by a single locus can lead to a stable equilibrium in which the two behaviours continue to segregate. A model in which the behaviours are determined by genes at two separate loci leads eventually to fixation of the alleles at both loci but this can take many generations of selection. Models of the type described in this paper will be needed to understand the evolution of complex behaviour when genomic or experimental information is available about the genetic determinants of behaviour and the selective values of different genomes. (c) 2007 Elsevier Inc. All rights reserved.

Towards a population ecology of stressed environments: the effects of zinc on the springtail Folsomia candida

1. To understand population dynamics in stressed environments it is necessary to join together two classical lines of research. Population responses to environmental stress have been studied at low density in life table response experiments. These show how the population's growth rate (pgr) at low density varies in relation to levels of stress. Population responses to density, on the other hand, are based on examination of the relationship between pgr and population density. 2. The joint effects of stress and density on pgr can be pictured as a contour map in which pgr varies with stress and density in the same way that the height of land above sea level varies with latitude and longitude. Here a microcosm experiment is reported that compared the joint effects of zinc and population density on the pgr of the springtail Folsomia candida (Collembola). 3. Our experiments allowed the plotting of a complete map of the effects of density and a stressor on pgr. Particularly important was the position of the pgr= 0 contour, which suggested that carrying capacity varied little with zinc concentration until toxic levels were reached. 4. This prediction accords well with observations of population abundance in the field. The method also allowed us to demonstrate, simultaneously, hormesis, toxicity, an Allee effect and density dependence. 5. The mechanisms responsible for these phenomena are discussed. As zinc is an essential trace element the initial increase in pgr is probably a consequence of dietary zinc deficiency. The Allee effect may be attributed to productivity of the environment increasing with density at low density. Density dependence is a result of food limitation. 6. Synthesis and applications. We illustrate a novel solution based on mapping a population's growth rate in relation to stress and population density. Our method allows us to demonstrate, simultaneously, hormesis, toxicity, an Allee effect and density dependence in an important ecological indicator species. We hope that the approach followed here will prove to have general applicability enabling predictions of field abundance to be made from estimates of the joint effects of the stressors and density on population growth rate.