scholarly journals Demographic history and divergence of sibling grouse species inferred from whole genome sequencing reveal past effects of climate change

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.

2019 ◽  
Vol 286 (1911) ◽  
pp. 20190304 ◽  
Author(s):  
Balaji Chattopadhyay ◽  
Kritika M. Garg ◽  
Rajasri Ray ◽  
Frank E. Rheindt

Over the last approximately 2.6 Myr, Earth's climate has been dominated by cyclical ice ages that have profoundly affected species' population sizes, but the impact of impending anthropogenic climate change on species’ extinction potential remains a worrying problem. We investigated 11 bat species from different taxonomic, ecological and geographical backgrounds using combined information from palaeoclimatic habitat reconstructions and genomes to analyse biotic impacts of historic climate change. We discover tightly correlated fluctuations between species' historic distribution and effective population size, identify frugivores as particularly susceptible to global warming, pinpoint large insectivores as having overall low effective population size and flag the onset of the Holocene (approx. 10–12 000 years ago) as the period with the generally lowest effective population sizes across the last approximately 1 Myr. Our study shows that combining genomic and palaeoclimatological approaches reveals effects of climatic shifts on genetic diversity and may help predict impacts of future climate change.


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.


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/.


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.


2015 ◽  
Author(s):  
Eva E Deinum ◽  
Daniel L Halligan ◽  
Rob W Ness ◽  
Yao-Hua Zhang ◽  
Lin Cong ◽  
...  

The brown rat, Rattus norvegicus, is both a notorious pest and a frequently used model in biomedical research. By analysing genome sequences of 12 wild-caught brown rats from their ancestral range in NE China, along with the sequence of a black rat, R. rattus, we investigate the selective and demographic forces shaping variation in the genome. We estimate that the recent effective population size (N_e) of this species = 1.24 x 10^5, based on silent site diversity. We compare patterns of diversity in these genomes with patterns in multiple genome sequences of the house mouse Mus musculus castaneus), which has a much larger N_e. This reveals an important role for variation in the strength of genetic drift in mammalian genome evolution. By a Pairwise Sequentially Markovian Coalescent (PSMC) analysis of demographic history, we infer that there has been a recent population size bottleneck in wild rats, which we date to approximately 20,000 years ago. Consistent with this, wild rat populations have experienced an increased flux of mildly deleterious mutations, which segregate at higher frequencies in protein-coding genes and conserved noncoding elements (CNEs). This leads to negative estimates of the rate of adaptive evolution (alpha) in proteins and CNEs, a result which we discuss in relation to the strongly positive estimates observed in wild house mice. As a consequence of the population bottleneck, wild rats also show a markedly slower decay of linkage disequilibrium with physical distance than wild house mice.


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.


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.


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.


2017 ◽  
Author(s):  
Aaron P. Ragsdale ◽  
Ryan N. Gutenkunst

AbstractPopulation demographic history may be learned from contemporary genetic variation data. Methods based on aggregating the statistics of many single loci into an allele frequency spectrum (AFS) have proven powerful, but such methods ignore potentially informative patterns of linkage disequilibrium (LD) between neighboring loci. To leverage such patterns, we developed a composite-likelihood framework for inferring demographic history from aggregated statistics of pairs of loci. Using this framework, we show that two-locus statistics are indeed more sensitive to demographic history than single-locus statistics such as the AFS. In particular, two-locus statistics escape the notorious confounding of depth and duration of a bottleneck, and they provide a means to estimate effective population size based on the recombination rather than mutation rate. We applied our approach to a Zambian population of Drosophila melanogaster. Notably, using both single– and two-locus statistics, we found substantially lower estimates of effective population size than previous works. Together, our results demonstrate the broad potential for two-locus statistics to enable powerful population genetic inference.


Sign in / Sign up

Export Citation Format

Share Document