Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long‐term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene–temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene–environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.     

Spatial and temporal population genetic structure of four northeastern Pacific littorinid gastropods: the effect of mode of larval development on variation at one mitochondrial and two nuclear DNA markers

We investigated the effect of development mode on the spatial and temporal population genetic structure of four littorinid gastropod species. Snails were collected from the same three sites on the west coast of Vancouver Island, Canada in 1997 and again in 2007. DNA sequences were obtained for one mitochondrial gene, cytochrome b ( Cyt b ), and for up to two nuclear genes, heat shock cognate 70 ( HSC70 ) and aminopeptidase N intron ( APN54 ). We found that the mean level of genetic diversity and long-term effective population sizes ( N _e) were significantly greater for two species, Littorina scutulata and L. plena , that had a planktotrophic larval stage than for two species, Littorina sitkana and L. subrotundata , that laid benthic egg masses which hatched directly into crawl-away juveniles. Predictably, two poorly dispersing species, L. sitkana and L. subrotundata , showed significant spatial genetic structure at an 11- to 65-km geographical scale that was not observed in the two planktotrophic species. Conversely, the two planktotrophic species had more temporal genetic structure over a 10-year interval than did the two direct-developing species and showed highly significant temporal structure for spatially pooled samples. The greater temporal genetic variation of the two planktotrophic species may have been caused by their high fecundity, high larval dispersal, and low but spatially correlated early survivorship. The sweepstakes-like reproductive success of the planktotrophic species could allow a few related females to populate hundreds of kilometres of coastline and may explain their substantially larger temporal genetic variance but lower spatial genetic variance relative to the direct-developing species.     

A likelihood-based comparison of population histories in a parasitoid guild

Little is known about the stability of trophic relationships in complex natural communities over evolutionary timescales. Here, we use sequence data from 18 nuclear loci to reconstruct and compare the intraspecific histories of major Pleistocene refugial populations in the Middle East, the Balkans and Iberia in a guild of four Chalcid parasitoids (Cecidostiba fungosa, Cecidostiba semifascia, Hobbya stenonota and Mesopolobus amaenus) all attacking Cynipid oak galls. We develop a likelihood method to numerically estimate models of divergence between three populations from multilocus data. We investigate the power of this framework on simulated data, and—using triplet alignments of intronic loci—quantify the support for all possible divergence relationships between refugial populations in the four parasitoids. Although an East to West order of population divergence has highest support in all but one species, we cannot rule out alternative population tree topologies. Comparing the estimated times of population splits between species, we find that one species, M. amaenus, has a significantly older history than the rest of the guild and must have arrived in central Europe at least one glacial cycle prior to other guild members. This suggests that although all four species may share a common origin in the East, they expanded westwards into Europe at different times.     

Investigating the evolutionary history of the Pacific Northwest mesic forest ecosystem: hypothesis testing within a comparative phylogeographic framework

We examine the evolution of mesic forest ecosystems in the Pacific Northwest of North America using a statistical phylogeography approach in four animal and two plant lineages. Three a priori hypotheses, which explain the disjunction in the mesic forest ecosystem with either recent dispersal or ancient vicariance, are tested with phylogenetic and coalescent methods. We find strong support in three amphibian lineages (Ascaphus spp., and Dicampton spp., and Plethodon vandykei and P. idahoensis) for deep divergence between coastal and inland populations, as predicted by the ancient vicariance hypothesis. Unlike the amphibians, the disjunction in other Pacific Northwest lineages is likely due to recent dispersal along a northern route. Topological and population divergence tests support the northern dispersal hypothesis in the water vole (Microtus richardsoni) and northern dispersal has some support in both the dusky willow (Salix melanopsis) and whitebark pine (Pinus albicaulis). These analyses demonstrate that genetic data sampled from across an ecosystem can provide insight into the evolution of ecological communities and suggest that the advantages of a statistical phylogeographic approach are most pronounced in comparisons across multiple taxa in a particular ecosystem. Genetic patterns in organisms as diverse as willows and salamanders can be used to test general regional hypotheses, providing a consistent metric for comparison among members of an ecosystem with disparate life-history traits.     

基于高通量测序的群体基因组学——植物病原真菌研究新方向

群体基因组学能够从全基因组水平揭示种群结构与进化、物种形成、适应性机制等。随着高通量技术的不断发展,基因组测序成本不断降低,大规模测序已成为可能。近几年被全基因组测序的真菌数量迅速增加,极大地促进了真菌群体基因组学的发展,加深了人们对植物病原真菌起源、遗传多样性、选择作用、致病性、毒力因子、杀菌剂抗药性、寄主专化型等生物学特性的认识。本文简要介绍了植物病原真菌的全基因测序以及比较基因组学的研究进展,重点综述了基于高通量测序的病原真菌群体基因组学的最新研究动态。群体基因组学将成为植物病原真菌一个新的研究方向。     

Did Late Pleistocene climate change result in parallel genetic structure and demographic bottlenecks in sympatric Central African crocodiles,Mecistops and Osteolaemus?

The mid‐Holocene has had profound demographic impacts on wildlife on the African continent, although there is little known about the impacts on species from Central Africa. Understanding the impacts of climate change on codistributed species can enhance our understanding of ecosystem dynamics and for formulating restoration objectives. We took a multigenome comparative approach to examine the phylogeographic structure of two poorly known Central African crocodile species—Mecistops sp. aff. cataphractus and Osteolaemus tetraspis. In addition, we conducted coalescent‐based demographic reconstructions to test the hypothesis that population decline was driven by climate change since the Last Glacial Maximum, vs. more recent anthropogenic pressures. Using a hierarchical Bayesian model to reconstruct demographic history, we show that both species had dramatic declines (>97%) in effective population size in the 'period following the Last Glacial Maximum 1,500–18,000 YBP. Identification of genetic structuring showed both species have similar regional structure corresponding to major geological features (i.e., hydrologic basin) and that small observed differences between them are best explained by the differences in their ecology and the likely impact that climate change had on their habitat needs. Our results support our hypothesis that climatic effects, presumably on forest and wetland habitat, had a congruent negative impact on both species.     

Genetic panmixia and demographic dependence across the North Atlantic in the deep-sea fish,blue hake (Antimora rostrata)

The efficient investment of resources and effort into conservation strategies depends on the accurate identification of management units. At the same time, understanding the processes by which population structure evolves requires an understanding of the conditions under which panmixia may exist. Here, we study a species with an unusual, apparently sex-biased pattern of distribution, and test the hypothesis that distribution processes associated with this pattern (for example, congregating at a single dominant spawning site or periodic mixing during reproduction) could lead to panmixia over a large geographic range. Using 13 microsatellite markers, we compared 393 blue hake (Antimora rostrata) from 11 sample sites across a geographic range of over 3000 km, and found no evidence of population structure. We estimated current effective population size and found it to be large (∼15 000) across the sampled area. In addition, we used simulation models to test expectations about demographic correlation among populations and our ability to detect relevant levels of gene flow. All data were consistent with the interpretation of long-range panmixia.     

From where did the Western honeybee (Apis mellifera) originate?

The native range of the honeybee Apis mellifera encompasses Europe, Africa, and the Middle East, whereas the nine other species of Apis are found exclusively in Asia. It is therefore commonly assumed that A. mellifera arose in Asia and expanded into Europe and Africa. However, other hypotheses for the origin of A. mellifera have also been proposed based on phylogenetic trees constructed from genetic markers. In particular, an analysis based on >1000 single-nucleotide polymorphism markers placed the root of the tree of A. mellifera subspecies among samples from Africa, suggestive of an out-of-Africa expansion. Here, we re-evaluate the evidence for this and other hypotheses by testing the robustness of the tree topology to different tree-building methods and by removing specimens with a potentially hybrid background. These analyses do not unequivocally place the root of the tree of A. mellifera subspecies within Africa, and are potentially consistent with a variety of hypotheses for honeybee evolution, including an expansion out of Asia. Our analyses also support high divergence between western and eastern European populations of A. mellifera, suggesting they are likely derived from two distinct colonization routes, although the sources of these expansions are still unclear.     

How to infer reliable diploid genotypes from NGS or traditional sequence data: from basic probability to experimental optimization

The use of diploid sequence markers is still challenging despite the good quality of the information they provide. There is a common problem to all sequencing approaches [traditional cloning and sequencing of PCR amplicons as well as next-generation sequencing (NGS)]: when no variation is found within the sequences from a given individual, homozygozity can never be asserted with certainty. As a consequence, sequence data from diploid markers are mostly analysed at the population (not the individual level) particularly in animal studies. This study aims at contributing to solve this. Using the Bayes theorem and the binomial law, useful results are derived, among which: (i) the number of sequence reads per individual (or sequencing depth) which is required to ensure, at a given probability threshold, that some heterozygotes are not considered erroneously as homozygotes, as a function of the observed heterozygozity (H(o) ) of the locus in the population; (ii) a way of estimating H(o) from low coverage NGS data; (iii) a way of testing the null hypothesis that a genetic marker corresponds to a single and diploid locus, in the absence of data from controlled crosses; (iv) strategies for characterizing sequence genotypes in populations minimizing the average number of sequence reads per individual; (v) a rationale to decide which are the variations that one needs to consider along the sequence, as a function of the sequencing depth affordable, the level of polymorphism desired and the risk of sequencing error. For traditional sequencing technology, optimal strategies appear surprisingly different from the usual empirical ones. The average number of sequence reads required to obtain 99% of fully determined genotypes never exceeds six, this value corresponding to the worst situation when H(o) equals 0.6. This threshold value of H(o) is strikingly stable when the tolerated proportion of nonfully resolved genotypes varies in a reasonable range. These results do not rely on the Hardy-Weinberg equilibrium assumption or on diallelism of nucleotidic sites.     

Pushing north one bottleneck at a time: site frequency spectra tell the history of Sitka spruce

Reconstructing the history of populations is a longstanding goal of molecular ecologists. In addition to a better understanding of the past, it is hoped that this knowledge would also facilitate predictions regarding species’ responses to future events such as climate change. The traditional way of doing this is through the fossil record, but these historical records are often incomplete. Inferring historical demography from patterns of nucleotide variability can help to fill these gaps. In this issue of Molecular Ecology, Holliday et al. (2010) glimpse into the demographic past of Sitka spruce, Picea sitchensis, an economically and ecologically important species native to northwestern United States and Canada, by examining the site frequency spectrum (SFS) of 153 loci in six populations covering the species entire range.