scholarly journals Robust estimation of recent effective population size from number of independent origins in soft sweeps

2018 ◽  
Author(s):  
Bhavin S. Khatri ◽  
Austin Burt

Estimating recent effective population size is of great importance in characterising and predicting the evolution of natural populations. Methods based on nucleotide diversity may underestimate current day effective population sizes due to historical bottlenecks, whilst methods that reconstruct demographic history typically only detect long-term variations. However, soft selective sweeps, which leave a fingerprint of mutational history by recurrent mutations on independent haplotype backgrounds, holds promise of an estimate more representative of recent population history. Here we present a simple and robust method of estimation based only on knowledge of the number of independent recurrent origins and the current frequency of the beneficial allele in a population sample, independent of the strength of selection and age of the mutation. Using a forward time theoretical framework, we show the mean number of origins is a function of θ = 2Nμ and current allele frequency, through a simple equation, and the distribution is approximately Poisson. This estimate is robust to whether mutants pre-existed before selection arose, and is equally accurate for diploid populations with incomplete dominance. For fast (e.g., seasonal) demographic changes compared to time scale for fixation of the mutant allele, and for moderate peak-to-trough ratios, we show our constant population size estimate can be used to bound the maximum and minimum population size. Applied to the Vgsc gene of Anopheles gambiae, we estimate an effective population size of roughly 6 × 107, and including seasonal demographic oscillations, a minimum effective population size greater than 6 × 106 and a maximum less than 3 × 109.

2019 ◽  
Vol 36 (9) ◽  
pp. 2040-2052 ◽  
Author(s):  
Bhavin S Khatri ◽  
Austin Burt

Abstract Estimating recent effective population size is of great importance in characterizing and predicting the evolution of natural populations. Methods based on nucleotide diversity may underestimate current day effective population sizes due to historical bottlenecks, whereas methods that reconstruct demographic history typically only detect long-term variations. However, soft selective sweeps, which leave a fingerprint of mutational history by recurrent mutations on independent haplotype backgrounds, holds promise of an estimate more representative of recent population history. Here, we present a simple and robust method of estimation based only on knowledge of the number of independent recurrent origins and the current frequency of the beneficial allele in a population sample, independent of the strength of selection and age of the mutation. Using a forward-time theoretical framework, we show the mean number of origins is a function of θ=2Nμ and current allele frequency, through a simple equation, and the distribution is approximately Poisson. This estimate is robust to whether mutants preexisted before selection arose and is equally accurate for diploid populations with incomplete dominance. For fast (e.g., seasonal) demographic changes compared with time scale for fixation of the mutant allele, and for moderate peak-to-trough ratios, we show our constant population size estimate can be used to bound the maximum and minimum population size. Applied to the Vgsc gene of Anopheles gambiae, we estimate an effective population size of roughly 6×107, and including seasonal demographic oscillations, a minimum effective population size >3×107, and a maximum <6×109, suggesting a mean ∼109.


Author(s):  
Andy Foote ◽  
Rebecca Hooper ◽  
Alana Alexander ◽  
Robin Baird ◽  
Charles Baker ◽  
...  

Runs of homozygosity (ROH) occur when offspring inherit haplotypes that are identical by descent from each parent. Length distributions of ROH are informative about population history; specifically the probability of inbreeding mediated by mating system and/or population demography. Here, we investigate whether variation in killer whale (Orcinus orca) demographic history is reflected in genome-wide heterozygosity and ROH length distributions, using a global dataset of 26 genomes representative of geographic and ecotypic variation in this species, and two F1 admixed individuals with Pacific-Atlantic parentage. We first reconstruct demographic history for each population as changes in effective population size through time using the pairwise sequential Markovian coalescent (PSMC) method. We find a subset of populations declined in effective population size during the Late Pleistocene, while others had more stable demography. Genomes inferred to have undergone ancestral declines in effective population size, were autozygous at hundreds of short ROH (<1Mb), reflecting high background relatedness due to coalescence of haplotypes deep within the pedigree. In contrast, longer and therefore younger ROH (>1.5 Mb) were found in low latitude populations and populations of known conservation concern, including a Scottish population, for which 37.8% of the autosomes comprised of ROH >1.5 Mb in length.


2021 ◽  
Vol 78 (2) ◽  
Author(s):  
Błażej Wójkiewicz ◽  
Andrzewj Lewandowski ◽  
Weronika B. Żukowska ◽  
Monika Litkowiec ◽  
Witold Wachowiak

Abstract Context Black poplar (Populus nigra L.) is a keystone species of European riparian ecosystems that has been negatively impacted by riverside urbanization for centuries. Consequently, it has become an endangered tree species in many European countries. The establishment of a suitable rescue plan of the remaining black poplar forest stands requires a preliminary knowledge about the distribution of genetic variation among species populations. However, for some parts of the P. nigra distribution in Europe, the genetic resources and demographic history remain poorly recognized. Aims Here, we present the first study on identifying and characterizing the genetic resources of black poplar from the Oder valley in Poland. This study (1) assessed the genetic variability and effective population size of populations and (2) examined whether gene flow is limited by distance or there is a single migrant pool along the studied river system. Methods A total of 582 poplar trees derived from nine black poplar populations were investigated with nuclear microsatellite markers. Results (1) The allelic richness and heterozygosity level were high and comparable between populations. (2) The genetic structure of the studied poplar stands was not homogenous. (3) The signatures of past bottlenecks were detected. Conclusion Our study (1) provides evidence for genetic substructuring of natural black poplar populations from the studied river catchment, which is not a frequent phenomenon reported for this species in Europe, and (2) indicates which poplar stands may serve as new genetic conservation units (GCUs) of this species in Europe. Key message The genetic resources of black poplar in the Oder River valley are still substantial compared to those reported for rivers in Western Europe. On the other hand, clear signals of isolation by distance and genetic erosion reflected in small effective population sizes and high spatial genetic structure of the analyzed populations were detected. Based on these findings, we recommend the in situ and ex situ conservation strategies for conserving and restoring the genetic resources of black poplar populations in this strongly transformed by human river valley ecosystem.


2010 ◽  
Vol 107 (5) ◽  
pp. 2147-2152 ◽  
Author(s):  
Chad D. Huff ◽  
Jinchuan Xing ◽  
Alan R. Rogers ◽  
David Witherspoon ◽  
Lynn B. Jorde

The genealogies of different genetic loci vary in depth. The deeper the genealogy, the greater the chance that it will include a rare event, such as the insertion of a mobile element. Therefore, the genealogy of a region that contains a mobile element is on average older than that of the rest of the genome. In a simple demographic model, the expected time to most recent common ancestor (TMRCA) is doubled if a rare insertion is present. We test this expectation by examining single nucleotide polymorphisms around polymorphic Alu insertions from two completely sequenced human genomes. The estimated TMRCA for regions containing a polymorphic insertion is two times larger than the genomic average (P < <10−30), as predicted. Because genealogies that contain polymorphic mobile elements are old, they are shaped largely by the forces of ancient population history and are insensitive to recent demographic events, such as bottlenecks and expansions. Remarkably, the information in just two human DNA sequences provides substantial information about ancient human population size. By comparing the likelihood of various demographic models, we estimate that the effective population size of human ancestors living before 1.2 million years ago was 18,500, and we can reject all models where the ancient effective population size was larger than 26,000. This result implies an unusually small population for a species spread across the entire Old World, particularly in light of the effective population sizes of chimpanzees (21,000) and gorillas (25,000), which each inhabit only one part of a single continent.


Genetics ◽  
1997 ◽  
Vol 147 (3) ◽  
pp. 1249-1258 ◽  
Author(s):  
Loren M Miller ◽  
Anne R Kapuscinski

Effective population size (Ne) of a natural fish population was estimated from temporal changes in allele frequencies at seven microsatellite loci. Use of a historical collection of fish scales made it possible to increase the precision of estimates by increasing the time interval between samples and to use an equation developed for discrete generations without correcting for demographic parameters. Estimates of Ne for the time intervals 1961–1977 and 1977–1993 were 35 and 72, respectively. For the entire interval, 1961–1993, the estimate of Ne was 48 when based on a weighted mean derived from the above two estimates or 125 when calculated from 1961 and 1993 samples only. Corresponding ratios of effective size to adult census size ranged from 0.03 to 0.14. An Ne of 48 over a 32-year period would imply that this population lost as much as 8% of its heterozygosity in that time. Results suggest the potential for using genetic methods based on microsatellite loci data to compare historical trends in Ne with population dynamic parameters. Such comparisons will help to evaluate the relationship between genetic diversity and long-term persistence of natural populations.


2017 ◽  
Author(s):  
Erik M. Volz ◽  
Xavier Didelot

AbstractNon-parametric population genetic modeling provides a simple and flexible approach for studying demographic history and epidemic dynamics using pathogen sequence data. Existing Bayesian approaches are premised on stationary stochastic processes which may provide an unrealistic prior for epidemic histories which feature extended period of exponential growth or decline. We show that non-parametric models defined in terms of the growth rate of the effective population size can provide a more realistic prior for epidemic history. We propose a non-parametric autoregressive model on the growth rate as a prior for effective population size, which corresponds to the dynamics expected under many epidemic situations. We demonstrate the use of this model within a Bayesian phylodynamic inference framework. Our method correctly reconstructs trends of epidemic growth and decline from pathogen genealogies even when genealogical data is sparse and conventional skyline estimators erroneously predict stable population size. We also propose a regression approach for relating growth rates of pathogen effective population size and time-varying variables that may impact the replicative fitness of a pathogen. The model is applied to real data from rabies virus and Staphylococcus aureus epidemics. We find a close correspondence between the estimated growth rates of a lineage of methicillin-resistant S. aureus and population-level prescription rates of β-lactam antibiotics. The new models are implemented in an open source R package called skygrowth which is available at https://mrc-ide.github.io/skygrowth/.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Kai Song ◽  
Bin Gao ◽  
Peter Halvarsson ◽  
Yun Fang ◽  
Siegfried Klaus ◽  
...  

Abstract Background The boreal forest is one of the largest biomes on earth, supporting thousands of species. The global climate fluctuations in the Quaternary, especially the ice ages, had a significant influence on the distribution of boreal forest, as well as the divergence and evolution of species inhabiting this biome. To understand the possible effects of on-going and future climate change it would be useful to reconstruct past population size changes and relate such to climatic events in the past. We sequenced the genomes of 32 individuals from two forest inhabiting bird species, Hazel Grouse (Tetrastes bonasia) and Chinese Grouse (T. sewerzowi) and three representatives of two outgroup species from Europe and China. Results We estimated the divergence time of Chinese Grouse and Hazel Grouse to 1.76 (0.46–3.37) MYA. The demographic history of different populations in these two sibling species was reconstructed, and showed that peaks and bottlenecks of effective population size occurred at different times for the two species. The northern Qilian population of Chinese Grouse became separated from the rest of the species residing in the south approximately 250,000 years ago and have since then showed consistently lower effective population size than the southern population. The Chinese Hazel Grouse population had a higher effective population size at the peak of the Last Glacial Period (approx. 300,000 years ago) than the European population. Both species have decreased recently and now have low effective population sizes. Conclusions Combined with the uplift history and reconstructed climate change during the Quaternary, our results support that cold-adapted grouse species diverged in response to changes in the distribution of palaeo-boreal forest and the formation of the Loess Plateau. The combined effects of climate change and an increased human pressure impose major threats to the survival and conservation of both species.


2017 ◽  
Author(s):  
John Hawks

AbstractHuman populations have a complex history of introgression and of changing population size. Human genetic variation has been affected by both these processes, so that inference of past population size depends upon the pattern of gene flow and introgression among past populations. One remarkable aspect of human population history as inferred from genetics is a consistent “wave” of larger effective population size, prior to the bottlenecks and expansions of the last 100,000 years. Here I carry out a series of simulations to investigate how introgression and gene flow from genetically divergent ancestral populations affect the inference of ancestral effective population size. Both introgression and gene flow from an extinct, genetically divergent population consistently produce a wave in the history of inferred effective population size. The time and amplitude of the wave reflect the time of origin of the genetically divergent ancestral populations and the strength of introgression or gene flow. These results demonstrate that even small fractions of introgression or gene flow from ancient populations may have large effects on the inference of effective population size.


2019 ◽  
Author(s):  
Xi Wang ◽  
Carolina Bernhardsson ◽  
Pär K. Ingvarsson

AbstractUnder the neutral theory, species with larger effective population sizes are expected to harbour higher genetic diversity. However, across a wide variety of organisms, the range of genetic diversity is orders of magnitude more narrow than the range of effective population size. This observation has become known as Lewontin’s paradox and although aspects of this phenomenon have been extensively studied, the underlying causes for the paradox remain unclear. Norway spruce (Picea abies) is a widely distributed conifer species across the northern hemisphere and it consequently plays a major role in European forestry. Here, we use whole-genome re-sequencing data from 35 individuals to perform population genomic analyses in P. abies in an effort to understand what drives genome-wide patterns of variation in this species. Despite having a very wide geographic distribution and an enormous current population size, our analyses find that genetic diversity of P.abies is low across a number of populations (p=0.005-0.006). To assess the reasons for the low levels of genetic diversity, we infer the demographic history of the species and find that it is characterised by several re-occurring bottlenecks with concomitant decreases in effective population size can, at least partly, provide an explanation for low polymorphism we observe in P. abies. Further analyses suggest that recurrent natural selection, both purifying and positive selection, can also contribute to the loss of genetic diversity in Norway spruce by reducing genetic diversity at linked sites. Finally, the overall low mutation rates seen in conifers can also help explain the low genetic diversity maintained in Norway spruce.


Sign in / Sign up

Export Citation Format

Share Document