Evidence for adaptation from standing genetic variation on an antimicrobial peptide gene in the mussel Mytilus edulis

Genome scans of population differentiation identify candidate loci for adaptation but provide little information on how selection has influenced the genetic structure of these loci. Following a genome scan, we investigated the nature of the selection responsible for the outlying differentiation observed between populations of the marine mussel Mytilus edulis at a leucine/arginine polymorphism (L31R) in the antimicrobial peptide MGD2. We analysed DNA sequence polymorphisms, allele frequencies and population differentiation of polymorphisms closely linked to L31R, and pairwise and third‐order linkage disequilibria. An outlying level of population differentiation was observed at L31R only, while no departure from panmixia was observed at linked loci surrounding L31R, as in most of the genome. Selection therefore seems to affect L31R directly. Three hypotheses can explain the lack of differentiation in the chromosomal region close to L31R: (i) hitchhiking has occurred but migration and recombination subsequently erased the signal, (ii) selection was weak enough and recombination strong enough to limit the hitchhiking effect to a very small chromosomal region or (iii) selection acted on a pre‐existing polymorphism (i.e. standing variation) at linkage equilibrium with its background. Linkage equilibrium was observed between L31R and linked polymorphisms in every population analysed, as expected under the three hypotheses. However, linkage disequilibrium was observed in some populations between pairs of loci located upstream and downstream to L31R, generating a complex pattern of third‐order linkage disequilibria which is best explained by the hypothesis of selection on a pre‐existing polymorphism. We hypothesise that selection could be either balanced, maintaining alleles at different frequencies depending on the pathogen community encountered locally by mussels, or intermittent, resulting in sporadic fluctuations in allele frequency.     

The genomic signature of parallel adaptation from shared genetic variation

Parallel adaptation is common and may often occur from shared genetic variation, but the genomic consequences of this process remain poorly understood. We first use individual‐based simulations to demonstrate that comparisons between populations adapted in parallel to similar environments from shared variation reveal a characteristic genomic signature around a selected locus: a low‐divergence valley centred at the locus and flanked by twin peaks of high divergence. This signature is initiated by the hitchhiking of haplotype tracts differing between derived populations in the broader neighbourhood of the selected locus (driving the high‐divergence twin peaks) and shared haplotype tracts in the tight neighbourhood of the locus (driving the low‐divergence valley). This initial hitchhiking signature is reinforced over time because the selected locus acts as a barrier to gene flow from the source to the derived populations, thus promoting divergence by drift in its close neighbourhood. We next empirically confirm the peak‐valley‐peak signature by combining targeted and RAD sequence data at three candidate adaptation genes in multiple marine (source) and freshwater (derived) populations of threespine stickleback. Finally, we use a genome‐wide screen for the peak‐valley‐peak signature to discover additional genome regions involved in parallel marine‐freshwater divergence. Our findings offer a new explanation for heterogeneous genomic divergence and thus challenge the standard view that peaks in population divergence harbour divergently selected loci and that low‐divergence regions result from balancing selection or localized introgression. We anticipate that genome scans for peak‐valley‐peak divergence signatures will promote the discovery of adaptation genes in other organisms.     

Lack of adaptation from standing genetic variation despite the presence of putatively adaptive alleles in introduced sweet vernal grass (Anthoxanthum odoratum)

Population genetic theory predicts that the availability of appropriate standing genetic variation should facilitate rapid evolution when species are introduced to new environments. However, few tests of rapid evolution have been paired with empirical surveys for the presence of previously identified adaptive genetic variants in natural populations. In this study, we examined local adaptation to soil Al toxicity in the introduced range of sweet vernal grass (Anthoxanthum odoratum), and we genotyped populations for the presence of Al tolerance alleles previously identified at the long‐term ecological Park Grass Experiment (PGE, Harpenden, UK) in the species native range. We found that markers associated with Al tolerance at the PGE were present at appreciable frequency in introduced populations. Despite this, there was no strong evidence of local adaptation to soil Al toxicity among populations. Populations demonstrated significantly different intrinsic root growth rates in the absence of Al. This suggests that selection on correlated root growth traits may constrain the ability of populations to evolve significantly different root growth responses to Al. Our results demonstrate that genotype–phenotype associations may differ substantially between the native and introduced parts of a species range and that adaptive alleles from a native species range may not necessarily promote phenotypic differentiation in the introduced range.     

The genetic basis of evolution, adaptation and speciation

A primary question in biology concerns the genetic basis of the evolution of novel traits, often in response to environmental changes, and how this can subsequently cause species isolation. This topic was the focus of the symposium on the Genetics of Speciation and Evolution at the annual meeting of the Canadian Society for Ecology and Evolution, held in Banff in May 2011. The presentations revealed some of the rapid advances being made in understanding the genetic basis of adaptation and speciation, as well as the elegant interplay between an organism's genetic complement and the environment that organism experiences.     

Phenotypic variation and covariation indicate high evolvability of acoustic communication in crickets

Studying the genetic architecture of sexual traits provides insight into the rate and direction at which traits can respond to selection. Traits associated with few loci and limited genetic and phenotypic constraints tend to evolve at high rates typically observed for secondary sexual characters. Here, we examined the genetic architecture of song traits and female song preferences in the field crickets Gryllus rubens and Gryllus texensis. Song and preference data were collected from both species and interspecific F1 and F2 hybrids. We first analysed phenotypic variation to examine interspecific differentiation and trait distributions in parental and hybrid generations. Then, the relative contribution of additive and additive‐dominance variation was estimated. Finally, phenotypic variance–covariance ( P ) matrices were estimated to evaluate the multivariate phenotype available for selection. Song traits and preferences had unimodal trait distributions, and hybrid offspring were intermediate with respect to the parents. We uncovered additive and dominance variation in song traits and preferences. For two song traits, we found evidence for X‐linked inheritance. On the one hand, the observed genetic architecture does not suggest rapid divergence, although sex linkage may have allowed for somewhat higher evolutionary rates. On the other hand, P matrices revealed that multivariate variation in song traits aligned with major dimensions in song preferences, suggesting a strong selection response. We also found strong covariance between the main traits that are sexually selected and traits that are not directly selected by females, providing an explanation for the striking multivariate divergence in male calling songs despite limited divergence in female preferences.     

The probability of parallel evolution

Abstract How often will natural selection drive parallel evolution at the DNA sequence level? More precisely, what is the probability that selection will cause two populations that live in identical environments to substitute the same beneficial mutation? Here I show that, under fairly general conditions, the answer is simple: if a wild‐type sequence can mutate to n different beneficial mutations, replicate populations will on average fix the same mutation with probability P= 2/(n + 1). This probability, which is derived using extreme value theory, is independent of most biological details, including the length of the gene in question and the precise distribution of fitness effects among alleles. I conclude that the probability of parallel evolution under natural selection is nearly twice as large as that under neutrality.     

The probability of genetic parallelism and convergence in natural populations

Genomic and genetic methods allow investigation of how frequently the same genes are used by different populations during adaptive evolution, yielding insights into the predictability of evolution at the genetic level. We estimated the probability of gene reuse in parallel and convergent phenotypic evolution in nature using data from published studies. The estimates are surprisingly high, with mean probabilities of 0.32 for genetic mapping studies and 0.55 for candidate gene studies. The probability declines with increasing age of the common ancestor of compared taxa, from about 0.8 for young nodes to 0.1–0.4 for the oldest nodes in our study. Probability of gene reuse is higher when populations begin from the same ancestor (genetic parallelism) than when they begin from divergent ancestors (genetic convergence). Our estimates are broadly consistent with genomic estimates of gene reuse during repeated adaptation to similar environments, but most genomic studies lack data on phenotypic traits affected. Frequent reuse of the same genes during repeated phenotypic evolution suggests that strong biases and constraints affect adaptive evolution, resulting in changes at a relatively small subset of available genes. Declines in the probability of gene reuse with increasing age suggest that these biases diverge with time.     

Influence of voltine ecotype and geographic distance on genetic and haplotype variation in the Asian corn borer

Diapause is an adaptive dormancy strategy by which arthropods endure extended periods of adverse climatic conditions. Seasonal variation in larval diapause initiation and duration in Ostrinia furnacalis may influence adult mating generation number (voltinism) across different local environments. The degree to which voltine ecotype, geographic distance, or other ecological factors influence O. furnacalis population genetic structure remains uncertain. Genetic differentiation was estimated between voltine ecotypes collected from 8 locations. Mitochondrial haplotypes were significantly different between historically allopatric univoltine and bivoltine locations, but confounded by a strong correlation with geographic distance. In contrast, single nucleotide polymorphism (SNP) genotypes show low but significant levels of variation and a lack of influence of geographic distance between allopatric voltine locations. Regardless, 11 of 257 SNP loci were predicted to be under selection, suggesting population genetic homogenization except at loci proximal to factors putatively under selection. These findings provide evidence of haplotype divergent voltine ecotypes that may be maintained in allopatric and sympatric areas despite relatively high rates of nuclear gene flow, yet influence of voltinism on maintenance of observed haplotype divergence remains unresolved.     

Notes from a snail island: Littorinid evolution and adaptation

The most successful study systems are built on a foundation of decades of research on the basic biology, ecology and life history of the organisms in question. Combined with new technologies, this can provide a formidable means to address important issues in evolutionary biology and molecular ecology. Littorinid marine snails are a good example of this, with a rich literature on their taxonomy, speciation, thermal tolerance and behavioural adaptations. In August 2017, an international meeting on Littorinid evolution was held at the Tjärnö Marine Research Laboratory in Western Sweden. In this meeting review, I provide a summary of some of the exciting work on parallel evolution, sexual selection and adaptation to environmental stress presented there. I argue that newly available genomic resources present an opportunity for integrating the traditionally divergent fields of speciation and environmental adaptation in Littorinid research.     

Genomic evidence for parallel adaptation to cities

Urban evolutionary biology is the study of rapid evolutionary change in response to humans and our uses of land to support city dwellers. Because cities are relatively modern additions to the natural world, research on urban evolution tends to focus on microevolutionary change that has happened across a few to many hundreds of generations. These questions still fall under the broad purview of evolutionary ecology. However, the severity, rapidity and replication of environmental changes that drive evolution in this context make it worthy of specific attention. Urban evolution provides the opportunity to study the earliest stages of evolution in a context that is scientifically interesting and societally important. The newness of urban populations and their proximity to natural populations also creates challenges when trying to detect population genetic change. In a From the Cover article in this issue of Molecular Ecology, Mueller et al. use whole genome resequencing data to address some of these challenges while exploring genetic changes associated with urbanization in three replicate urban‐rural burrowing owl (Athene cunicularia) populations. Combining multiple approaches across these sample sites Mueller et al. find evidence for selection on genes whose function is related to synapses, neuron projections, brain connectivity and cognitive function in general. That selection was parallel suggests that phenotypes related to brain processes were probably particularly important for urban adaptation.     

Assessing footprints of selection in commercial Atlantic salmon populations using microsatellite data

Relatively large rates of response to traits of economic importance have been observed in different selection experiments in salmon. Several QTL have been mapped in the salmon genome, explaining unprecedented levels of phenotypic variation. Owing to the relatively large selection intensity, individual loci may be indirectly selected, leaving molecular footprints of selection, together with increased inbreeding, as its likely relatives will share the selected loci. We used population differentiation and levels of linkage disequilibrium in chromosomes known to be harbouring QTL for body weight, infectious pancreatic necrosis resistance and infectious salmon anaemia resistance to assess the recent selection history at the genomic level in Atlantic salmon. The results clearly suggest that the marker SSA0343BSFU on chromosome 3 (body weight QTL) showed strong evidence of directional selection. It is intriguing that this marker is physically mapped to a region near the coding sequence of DVL2 , making it an ideal candidate gene to explain the rapid evolutionary response of this chromosome to selection for growth in Salmo salar. Weak evidence of diversifying selection was observed in the QTL associated with infectious pancreatic necrosis and infectious salmon anaemia resistance. Overall, this study showed that artificial selection has produced important changes in the Atlantic salmon genome, validating QTL in commercial salmon populations used for production purposes according to the recent selection history.     

The evolution of genetic canalization under fluctuating selection

Abstract.— If the direction of selection changes from generation to generation, the ability to respond to selection is maladaptive: the response to selection in one generation leads to reduced fitness in the next. Because the response is determined by the amount of genetic variance expressed at the phenotypic level, rapidly fluctuating selection should favor modifier genes that reduce the phenotypic effect of alleles segregating at structural loci underlying the trait. Such reduction in phenotypic expression of genetic variation has been named "genetic canalization." I support this argument with a series of two- and multilocus models with alternating linear selection and Gaussian selection with fluctuating optimum. A canalizing modifier gene affects the fitness of its carriers in three ways: (1) it reduces the phenotypic consequences of genetic response to previous selection; (2) it reduces the genetic response to selection, which is manifested as linkage disequilibrium between the modifier and structural loci; and (3) it reduces the phenotypic variance. The first two effects reduce fitness under directional selection sustained for several generations, but improve fitness when the direction of selection has just been reversed. The net effect tends to favor a canalizing modifier under rapidly fluctuating selection regimes (period of eight generations or less). The third effect improves fitness of the modifier allele if the fitness function is convex and reduces it if the function is concave. Under fluctuating Gaussian selection, the population is more likely to experience the concave portion of the fitness function when selection is stronger. Therefore, only weak to moderately strong fluctuating Gaussian selection favors genetic canalization. This paper considerably broadens the conditions that favor genetic canalization, which so far has only been postulated to evolve under long-term stabilizing selection.     

Distribution of genetic variation underlying adult migration timing in steelhead of the Columbia River basin

Fish migrations are energetically costly, especially when moving between freshwater and saltwater, but are a viable strategy for Pacific salmon and trout (Oncorhynchus spp.) due to the advantageous resources available at various life stages. Anadromous steelhead (O. mykiss) migrate vast distances and exhibit variation for adult migration phenotypes that have a genetic basis at candidate genes known as greb1L and rock1. We examined the distribution of genetic variation at 13 candidate markers spanning greb1L, intergenic, and rock1 regions versus 226 neutral markers for 113 populations (n = 9,471) of steelhead from inland and coastal lineages in the Columbia River. Patterns of population structure with neutral markers reflected genetic similarity by geographic region as demonstrated in previous studies, but candidate markers clustered populations by genetic variation associated with adult migration timing. Mature alleles for late migration had the highest frequency overall in steelhead populations throughout the Columbia River, with only 9 of 113 populations that had a higher frequency of premature alleles for early migration. While a single haplotype block was evident for the coastal lineage, we identified multiple haplotype blocks for the inland lineage. The inland lineage had one haplotype block that corresponded to candidate markers within the greb1L gene and immediately upstream in the intergenic region, and the second block only contained candidate markers from the intergenic region. Haplotype frequencies had similar patterns of geographic distribution as single markers, but there were distinct differences in frequency between the two haplotype blocks for the inland lineage. This may represent multiple recombination events that differed between lineages where phenotypic differences exist between freshwater entry versus arrival timing as indicated by Micheletti et al. (2018a). Redundancy analyses were used to model environmental effects on allelic frequencies of candidate markers, and significant variables were migration distance, temperature, isothermality, and annual precipitation. This study improves our understanding of the spatial distribution of genetic variation underlying adult migration timing in steelhead as well as associated environmental factors and has direct conservation and management implications.     

The evolution of genetic architecture. I. Diversification of genetic backgrounds by genetic drift

The genetic architecture of a phenotype plays a critical role in determining phenotypic evolution through its effects on patterns of genetic variation. Genetic architecture is often considered to be constant in evolutionary quantitative genetic models. However, genetic architecture may be variable and itself evolve when there are dominance and epistatic interactions among alleles at the same and different loci, respectively. The evolution of genetic architecture by genetic drift is examined here by testing the breeding value of four standard inbred mouse strains mated across a set of 26 related recombinant quasi-inbred (RqI) lines generated from the intercross of the Large (LG/J) and Small (SM/J) inbred mouse strains. Phenotypes of interest include age-specific body weights, growth, and adult body composition. If the genetic architecture of these traits has differentiated by genetic drift during the production of the RqI strains, we should observe interactions between tester strain and RqI strain. The breeding values of the tester strains will change relative to one another depending on which RqI strain they are crossed to. The study included an average of 15.1 offspring per cross, over a total of 100 different crosses. Multivariate and univariate analyses of variance indicate that there is strongly significant interaction for all traits. Interaction is more pronounced in males than in females and accounted for an average of about 40% of the explained variation in males and 30% in females. These results indicate that the genetic architecture of these traits has differentiated by genetic drift in the RqI strains since their isolation from a common founder population. Further analysis indicates that this differentiation results in changes in the order of tester strain effects so that common patterns of selection in these differentiated populations could result in the fixation of different alleles.     

Comparative landscape genetics reveals the evolution of viviparity reduces genetic connectivity in fire salamanders

Evolutionary changes in reproductive mode may affect co‐evolving traits, such as dispersal, although this subject remains largely underexplored. The shift from aquatic oviparous or larviparous reproduction to terrestrial viviparous reproduction in some amphibians entails skipping the aquatic larval stage and, thus, greater independence from water. Accordingly, amphibians exhibiting terrestrial viviparous reproduction may potentially disperse across a wider variety of suboptimal habitats and increase population connectivity in fragmented landscapes compared to aquatic‐breeding species. We investigated this hypothesis in the fire salamander (Salamandra salamandra), which exhibits both aquatic‐ (larviparity) and terrestrial‐breeding (viviparity) strategies. We genotyped 426 larviparous and 360 viviparous adult salamanders for 13 microsatellite loci and sequenced a mitochondrial marker for 133 larviparous and 119 viviparous individuals to compare population connectivity and landscape resistance to gene flow within a landscape genetics framework. Contrary to our predictions, viviparous populations exhibited greater differentiation and reduced genetic connectivity compared to larviparous populations. Landscape genetic analyses indicate viviparity may be partially responsible for this pattern, as water courses comprised a significant barrier only in viviparous salamanders, probably due to their fully terrestrial life cycle. Agricultural areas and, to a lesser extent, topography also decreased genetic connectivity in both larviparous and viviparous populations. This study is one of very few to explicitly demonstrate the evolution of a derived reproductive mode affects patterns of genetic connectivity. Our findings open avenues for future research to better understand the eco‐evolutionary implications underlying the emergence of terrestrial reproduction in amphibians.     

A practical introduction to Random Forest for genetic association studies in ecology and evolution

Large genomic studies are becoming increasingly common with advances in sequencing technology, and our ability to understand how genomic variation influences phenotypic variation between individuals has never been greater. The exploration of such relationships first requires the identification of associations between molecular markers and phenotypes. Here, we explore the use of Random Forest (RF), a powerful machine‐learning algorithm, in genomic studies to discern loci underlying both discrete and quantitative traits, particularly when studying wild or nonmodel organisms. RF is becoming increasingly used in ecological and population genetics because, unlike traditional methods, it can efficiently analyse thousands of loci simultaneously and account for nonadditive interactions. However, understanding both the power and limitations of Random Forest is important for its proper implementation and the interpretation of results. We therefore provide a practical introduction to the algorithm and its use for identifying associations between molecular markers and phenotypes, discussing such topics as data limitations, algorithm initiation and optimization, as well as interpretation. We also provide short R tutorials as examples, with the aim of providing a guide to the implementation of the algorithm. Topics discussed here are intended to serve as an entry point for molecular ecologists interested in employing Random Forest to identify trait associations in genomic data sets.     

The signature of positive selection on standing genetic variation

Considerable interest is focused on the use of polymorphism data to identify regions of the genome that underlie recent adaptations. These searches are guided by a simple model of positive selection, in which a mutation is favored as soon as it arises. This assumption may not be realistic, as environmental changes and range expansions may lead previously neutral or deleterious alleles to become beneficial. We examine what effect this mode of selection has on patterns of variation at linked neutral sites by implementing a new coalescent model of positive directional selection on standing variation. In this model, a neutral allele arises and drifts in the population, then at frequency f becomes beneficial, and eventually reaches fixation. Depending on the value of f, this scenario can lead to a large variance in allele frequency spectra and in levels of linkage disequilibrium at linked, neutral sites. In particular, for intermediate f, the beneficial substitution often leads to a loss of rare alleles--a pattern that differs markedly from the signature of directional selection currently relied on by researchers. These findings highlight the importance of an accurate characterization of the effects of positive selection, if we are to reliably identify recent adaptations from polymorphism data.     

The maintenance of heritable variation through social competition

The paradoxical persistence of heritable variation for fitness-related traits is an evolutionary conundrum that remains a preeminent problem in evolutionary biology. Here we describe a simple mechanism in which social competition results in the evolutionary maintenance of heritable variation for fitness related traits. We demonstrate this mechanism using a genetic model with two primary assumptions: the expression of a trait depends upon success in social competition for limited resources; and competitive success of a genotype depends on the genotypes that it competes against. We find that such social competition generates heritable (additive) genetic variation for "competition-dependent" traits. This heritable variation is not eroded by continuous directional selection because, rather than leading to fixation of favored alleles, selection leads instead to allele frequency cycling due to the concerted coevolution of the social environment with the effects of alleles. Our results provide a mechanism for the maintenance of heritable variation in natural populations and suggest an area for research into the importance of competition in the genetic architecture of fitness related traits.