Search PubMed⌕ Search

Biomedical subjects

Mark A Beaumont

Publications and source records attributed to Mark A Beaumont.

11 recordsLinked to original sources

Geographical variation in the clouded leopard, Neofelis nebulosa, reveals two species.

The clouded leopard, Neofelis nebulosa, is an endangered semiarboreal felid with a wide distribution in tropical forests of southern and southeast Asia, including the islands of Sumatra and Borneo in the Indonesian archipelago. In common with many larger animal species, it displays morphological variation within its wide geographical range and is currently regarded as comprising of up to four subspecies. It is widely recognized that taxonomic designation has a major impact on conservation planning and action. Given that the last taxonomic revision was made over 50 years ago, a more detailed examination of geographical variation is needed. We describe here the results of a morphometric analysis of the pelages of 57 clouded leopards sampled throughout the species' range. We conclude that there are two distinct morphological groups, which differ primarily in the size of their cloud markings. These results are supported by a recent genetic analysis. On that basis, we give diagnoses for the distinction of two species, one in mainland Asia (N. nebulosa) and the other in Indonesia (N. diardi). The implications for conservation that arise from this new taxonomic arrangement are discussed.

Animals↗

Adaptation and speciation: what can F(st) tell us?

A useful way of summarizing genetic variability among different populations is through estimates of the inbreeding coefficient, F(st). Several recent studies have tried to use the distribution of estimates of F(st) from individual genetic loci to detect the effects of natural selection. However, the promise of this approach has yet to be fully realized owing to the pervasive dogma that this distribution is highly dependent on demographic history. Here, I review recent theoretical results that indicate that the distribution of estimates of F(st) is generally expected to be robust to the vagaries of demographic history. I suggest that analyses based on it provide a useful first step for identifying candidate genes that might be under selection, and explore the ways in which this information can be used in ecological and evolutionary studies.

Journal Article↗

Identifying adaptive genetic divergence among populations from genome scans.

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of FST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of FST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of FST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

Adaptation, Biological↗

Comparative evaluation of a new effective population size estimator based on approximate bayesian computation.

We describe and evaluate a new estimator of the effective population size (N(e)), a critical parameter in evolutionary and conservation biology. This new "SummStat" N(e) estimator is based upon the use of summary statistics in an approximate Bayesian computation framework to infer N(e). Simulations of a Wright-Fisher population with known N(e) show that the SummStat estimator is useful across a realistic range of individuals and loci sampled, generations between samples, and N(e) values. We also address the paucity of information about the relative performance of N(e) estimators by comparing the SummStat estimator to two recently developed likelihood-based estimators and a traditional moment-based estimator. The SummStat estimator is the least biased of the four estimators compared. In 32 of 36 parameter combinations investigated using initial allele frequencies drawn from a Dirichlet distribution, it has the lowest bias. The relative mean square error (RMSE) of the SummStat estimator was generally intermediate to the others. All of the estimators had RMSE > 1 when small samples (n = 20, five loci) were collected a generation apart. In contrast, when samples were separated by three or more generations and N(e) < or = 50, the SummStat and likelihood-based estimators all had greatly reduced RMSE. Under the conditions simulated, SummStat confidence intervals were more conservative than the likelihood-based estimators and more likely to include true N(e). The greatest strength of the SummStat estimator is its flexible structure. This flexibility allows it to incorporate any potentially informative summary statistic from population genetic data.

Bayes Theorem↗

Estimation of population growth or decline in genetically monitored populations.

This article introduces a new general method for genealogical inference that samples independent genealogical histories using importance sampling (IS) and then samples other parameters with Markov chain Monte Carlo (MCMC). It is then possible to more easily utilize the advantages of importance sampling in a fully Bayesian framework. The method is applied to the problem of estimating recent changes in effective population size from temporally spaced gene frequency data. The method gives the posterior distribution of effective population size at the time of the oldest sample and at the time of the most recent sample, assuming a model of exponential growth or decline during the interval. The effect of changes in number of alleles, number of loci, and sample size on the accuracy of the method is described using test simulations, and it is concluded that these have an approximately equivalent effect. The method is used on three example data sets and problems in interpreting the posterior densities are highlighted and discussed.

Alleles↗

Y genetic data support the Neolithic demic diffusion model.

There still is no general agreement on the origins of the European gene pool, even though Europe has been more thoroughly investigated than any other continent. In particular, there is continuing controversy about the relative contributions of European Palaeolithic hunter-gatherers and of migrant Near Eastern Neolithic farmers, who brought agriculture to Europe. Here, we apply a statistical framework that we have developed to obtain direct estimates of the contribution of these two groups at the time they met. We analyze a large dataset of 22 binary markers from the non-recombining region of the Y chromosome (NRY), by using a genealogical likelihood-based approach. The results reveal a significantly larger genetic contribution from Neolithic farmers than did previous indirect approaches based on the distribution of haplotypes selected by using post hoc criteria. We detect a significant decrease in admixture across the entire range between the Near East and Western Europe. We also argue that local hunter-gatherers contributed less than 30% in the original settlements. This finding leads us to reject a predominantly cultural transmission of agriculture. Instead, we argue that the demic diffusion model introduced by Ammerman and Cavalli-Sforza [Ammerman, A. J. & Cavalli-Sforza, L. L. (1984) The Neolithic Transition and the Genetics of Populations in Europe (Princeton Univ. Press, Princeton)] captures the major features of this dramatic episode in European prehistory.

Animals↗

Likelihood-based estimation of the effective population size using temporal changes in allele frequencies: a genealogical approach.

A new genetic estimator of the effective population size (N(e)) is introduced. This likelihood-based (LB) estimator uses two temporally spaced genetic samples of individuals from a population. We compared its performance to that of the classical F-statistic-based N(e) estimator (N(eFk)) by using data from simulated populations with known N(e) and real populations. The new likelihood-based estimator (N(eLB)) showed narrower credible intervals and greater accuracy than (N(eFk)) when genetic drift was strong, but performed only slightly better when genetic drift was relatively weak. When drift was strong (e.g., N(e) = 20 for five generations), as few as approximately 10 loci (heterozygosity of 0.6; samples of 30 individuals) are sufficient to consistently achieve credible intervals with an upper limit <50 using the LB method. In contrast, approximately 20 loci are required for the same precision when using the classical F-statistic approach. The N(eLB) estimator is much improved over the classical method when there are many rare alleles. It will be especially useful in conservation biology because it less often overestimates N(e) than does N(eLB) and thus is less likely to erroneously suggest that a population is large and has a low extinction risk.

Animals↗

Approximate Bayesian computation in population genetics.

We propose a new method for approximate Bayesian statistical inference on the basis of summary statistics. The method is suited to complex problems that arise in population genetics, extending ideas developed in this setting by earlier authors. Properties of the posterior distribution of a parameter, such as its mean or density curve, are approximated without explicit likelihood calculations. This is achieved by fitting a local-linear regression of simulated parameter values on simulated summary statistics, and then substituting the observed summary statistics into the regression equation. The method combines many of the advantages of Bayesian statistical inference with the computational efficiency of methods based on summary statistics. A key advantage of the method is that the nuisance parameters are automatically integrated out in the simulation step, so that the large numbers of nuisance parameters that arise in population genetics problems can be handled without difficulty. Simulation results indicate computational and statistical efficiency that compares favorably with those of alternative methods previously proposed in the literature. We also compare the relative efficiency of inferences obtained using methods based on summary statistics with those obtained directly from the data using MCMC.

Bayes Theorem↗

Genetic evidence for long-term population decline in a savannah-dwelling primate: inferences from a hierarchical bayesian model.

The purpose of this study was to test for evidence that savannah baboons (Papio cynocephalus) underwent a population expansion in concert with a hypothesized expansion of African human and chimpanzee populations during the late Pleistocene. The rationale is that any type of environmental event sufficient to cause simultaneous population expansions in African humans and chimpanzees would also be expected to affect other codistributed mammals. To test for genetic evidence of population expansion or contraction, we performed a coalescent analysis of multilocus microsatellite data using a hierarchical Bayesian model. Markov chain Monte Carlo (MCMC) simulations were used to estimate the posterior probability density of demographic and genealogical parameters. The model was designed to allow interlocus variation in mutational and demographic parameters, which made it possible to detect aberrant patterns of variation at individual loci that could result from heterogeneity in mutational dynamics or from the effects of selection at linked sites. Results of the MCMC simulations were consistent with zero variance in demographic parameters among loci, but there was evidence for a 10- to 20-fold difference in mutation rate between the most slowly and most rapidly evolving loci. Results of the model provided strong evidence that savannah baboons have undergone a long-term historical decline in population size. The mode of the highest posterior density for the joint distribution of current and ancestral population size indicated a roughly eightfold contraction over the past 1,000 to 250,000 years. These results indicate that savannah baboons apparently did not share a common demographic history with other codistributed primate species.

Animals↗

Testing for genetic evidence of population expansion and contraction: an empirical analysis of microsatellite DNA variation using a hierarchical Bayesian model.

The role of past climatic change in shaping the distributions of tropical rain forest vertebrates is central to long-standing hypotheses about the legacy of the Quaternary ice ages. One approach to testing such hypotheses is to use genetic data to infer the demographic history of codistributed species. Population genetic theory that relates the structure of allelic genealogies to historical changes in effective population size can be used to detect a past history of demographic expansion or contraction. The fruit bats Cynopterus sphinx and C. brachyotis (Chiroptera: Pteropodidae) exhibit markedly different distribution patterns across the Indomalayan region and therefore represent an exemplary species pair to use for such tests. The purpose of this study was to test alternative hypotheses about historical patterns of demographic expansion and contraction in C. sphinx and C. brachyotis using a coalescent-based analysis of microsatellite variation. Specifically, we used a hierarchical Bayesian model based on Markov chain Monte Carlo simulations to estimate the posterior distribution of genealogical and demographic parameters. The results revealed strong evidence for population contraction in both species. Evidence for a population contraction in C. brachyotis was expected on the basis of biogeographic considerations. However, similar evidence for population contraction in C. sphinx does not support the hypothesis that this species underwent a pronounced range expansion during the late Quaternary. Genetic evidence for population decline may reflect the consequences of habitat destruction on a more recent time scale.

Animals↗