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.     

Adaptation dynamics and evolutionary rescue under sulfide selection in cyanobacteria: a comparative study between Microcystis aeruginosa and Oscillatoria sp. (cyanobacteria)

Experimental evolution studies using cyanobacteria as model organisms are scarce despite the role of cyanobacteria in the evolution of photosynthesis. Three different experimental evolution approaches have been applied to shed light on the sulfide adaptation process, which played a key role in the evolution of this group. We used a Microcystis aeruginosa sulfide‐sensitive strain, unable to grow above ~0.1 mM, and an Oscillatoria sp. strain, isolated from a sulfureous spa (~0.2 mM total sulfide). First, performing a fluctuation analysis design using the spa waters as selective agent, we proved that M. aeruginosa was able to adapt to this sulfide level by rare spontaneous mutations. Second, applying a ratchet protocol, we tested if the limit of adaptation to sulfide of the two taxa was dependent on their initial sulfide tolerance, finding that M. aeruginosa adapted to 0.4 mM sulfide, and Oscillatoria sp. to ~2 mM sulfide, twice it highest tolerance level. Third, using an evolutionary rescue approach, we observed that both speed of exposure to increasing sulfide concentrations (deterioration rate) and populations’ genetic variation determined the survival of M. aeruginosa at lethal sulfide levels, with a higher dependence on genetic diversity. In conclusion, sulfide adaptation of sensitive cyanobacterial strains is possible by rare spontaneous mutations and the adaptation limits depend on the sulfide level present in strain’s original habitat. The high genetic diversity of a sulfide‐sensitive strain, even at fast environmental deterioration rates, could increase its possibility of survival even to a severe sulfide stress.     

Restricted gene flow and local adaptation highlight the vulnerability of high‐latitude reefs to rapid environmental change

Global climate change poses a serious threat to the future health of coral reef ecosystems. This calls for management strategies that are focused on maximizing the evolutionary potential of coral reefs. Fundamental to this is an accurate understanding of the spatial genetic structure in dominant reef‐building coral species. In this study, we apply a genotyping‐by‐sequencing approach to investigate genome‐wide patterns of genetic diversity, gene flow, and local adaptation in a reef‐building coral, Pocillopora damicornis, across 10 degrees of latitude and a transition from temperate to tropical waters. We identified strong patterns of differentiation and reduced genetic diversity in high‐latitude populations. In addition, genome‐wide scans for selection identified a number of outlier loci putatively under directional selection with homology to proteins previously known to be involved in heat tolerance in corals and associated with processes such as photoprotection, protein degradation, and immunity. This study provides genomic evidence for both restricted gene flow and local adaptation in a widely distributed coral species, and highlights the potential vulnerability of leading‐edge populations to rapid environmental change as they are locally adapted, reproductively isolated, and have reduced levels of genetic diversity.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site‐specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole‐genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree‐based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.     

On the potential strength and consequences for nonrandom gene flow caused by local adaptation in flowering time

Gene flow is generally considered a random process, that is the loci under consideration have no effect on dispersal success. Edelaar and Bolnick (Trends Ecol Evol, 27, 2012 659) recently argued that nonrandom gene flow could exert a significant evolutionary force. It can, for instance, ameliorate the maladaptive effects of immigration into locally adapted populations. I examined the potential strength for nonrandom gene flow for flowering time genes, a trait frequently found to be locally adapted. The idea is that plants that successfully export pollen into a locally adapted resident population will be a genetically biased subset of their natal population – they will have resident‐like flowering times. Reciprocally, recipients will be more migrant‐like than the resident population average. I quantified the potential for biased pollen exchange among three populations along a flowering time cline in Brassica rapa from southern California. A two‐generation line cross experiment demonstrated genetic variance in flowering time, both within and among populations. Calculations based on the variation in individual flowering schedules showed that resident plants with the most migrant‐like flowering times could expect to have up to 10 times more of the their flowers pollinated by immigrant pollen than the least migrant‐like. Further, the mean flowering time of the pollen exporters that have access to resident mates differs by up to 4 weeks from the mean in the exporters’ natal population. The data from these three populations suggest that the bias in gene flow for flowering time cuts the impact on the resident population by as much as half. This implies that when selection is divergent between populations, migrants with the highest mating success tend to be resident‐like in their flowering times, and so, fewer maladaptive alleles will be introduced into the locally adapting gene pool.     

Evolutionary changes in plant tolerance against herbivory through a resurrection experiment

Both theoretical and empirical works have highlighted the difference in the evolutionary implications of host resistance and tolerance against their enemies. However, it has been difficult to show evolutionary changes in host defences in natural populations; thus, evaluating theoretical predictions of simultaneous evolution of defences remains a challenge. We studied the evolutionary changes in traits related to resistance and tolerance against herbivory in a natural plant population using seeds from two collections made in a period of 20 years. In a common garden experiment, we compared defensive traits of ancestral (1987) and descendant (2007) subpopulations of the annual plant Datura stramonium that shows genetic variation for tolerance and to which the specialist herbivore Lema daturaphila is locally adapted. We also examined the effects of different plant genotypes on the herbivore for testing the plant genetic variation in resistance. Based on the response to the contemporary herbivore populations, results revealed a nonsignificant response in plant resistance traits (herbivore consumption, foliar trichomes and tropane alkaloids), but a significant one in tolerance. The survival of herbivores in laboratory experiments depended on the plant genotype, which suggests genetic variation in plant resistance. Although we cannot identify the selective agent for the change nor exclude genetic drift, the results are consistent with the expectation that when resistance fails to control herbivory, tolerance should play a more important role in the evolution of the interaction.     

Signatures of local adaptation along environmental gradients in a range‐expanding damselfly (Ischnura elegans)

Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range‐expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site‐associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V‐ATPase) and visual processes (long‐wavelength‐sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate‐mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients.     

Genotype‐by‐environment interactions due to antibiotic resistance and adaptation in Escherichia coli

Mutations that are beneficial in one environment can have different fitness effects in other environments. In the context of antibiotic resistance, the resulting genotype‐by‐environment interactions potentially make selection on resistance unpredictable in heterogeneous environments. Furthermore, resistant bacteria frequently fix additional mutations during evolution in the absence of antibiotics. How do these two types of mutations interact to determine the bacterial phenotype across different environments? To address this, I used Escherichia coli as a model system, measuring the effects of nine different rifampicin resistance mutations on bacterial growth in 31 antibiotic‐free environments. I did this both before and after approximately 200 generations of experimental evolution in antibiotic‐free conditions (LB medium), and did the same for the antibiotic‐sensitive wild type after adaptation to the same environment. The following results were observed: (i) bacteria with and without costly resistance mutations adapted to experimental conditions and reached similar levels of competitive fitness; (ii) rifampicin resistance mutations and adaptation to LB both indirectly altered growth in other environments; and (iii) resistant‐evolved genotypes were more phenotypically different from the ancestor and from each other than resistant‐nonevolved and sensitive‐evolved genotypes. This suggests genotype‐by‐environment interactions generated by antibiotic resistance mutations, observed previously in short‐term experiments, are more pronounced after adaptation to other types of environmental variation, making it difficult to predict long‐term selection on resistance mutations from fitness effects in a single environment.     

A recent local sweep at the PHYA locus in the Northern European Spiterstulen population of Arabidopsis lyrata

Northern and central European Arabidopsis lyrata ssp. petraea populations are locally adapted to prevailing climatic conditions through differences in timing of life history events. The timing of flowering and, in perennials, the timing of growth cessation influence fitness. Phytochrome A may have an important role in regulating these life history traits as it perceives changes in daylength. We asked whether PHYA has contributed to local adaptation to the northern conditions in A. l. petraea. To search for signals of directional selection at the PHYA locus, we resequenced PHYA and 9 short fragments around PHYA from a 57‐kb region from a German (Plech) and a Norwegian (Spiterstulen) population and compared patterns of differentiation and diversity to a set of 19 reference loci around the genome. First, we found that the populations were highly differentiated: there were three nonsynonymous fixed differences at the PHYA locus, which was in stark contrast with the total four fixed differences in the 19 reference loci. Compatible with a sweep hypothesis, variation was almost completely removed from the 9.4‐kb region around PHYA in the northern Spiterstulen population. The overall level of linkage disequilibrium (LD) was higher in Spiterstulen, but there was no LD across the PHYA locus in the population, which is also a known consequence of a selective sweep. The sweep has likely occurred after the last glacial maximum, which suggests that it has contributed to adaptation to the northern conditions.     

Dynamics of genomic change during evolutionary rescue in the seed beetle Callosobruchus maculatus

Rapid adaptation can prevent extinction when populations are exposed to extremely marginal or stressful environments. Factors that affect the likelihood of evolutionary rescue from extinction have been identified, but much less is known about the evolutionary dynamics (e.g., rates and patterns of allele frequency change) and genomic basis of successful rescue, particularly in multicellular organisms. We conducted an evolve‐and‐resequence experiment to investigate the dynamics of evolutionary rescue at the genetic level in the cowpea seed beetle, Callosobruchus maculatus, when it is experimentally shifted to a stressful host plant, lentil. Low survival (~1%) at the onset of the experiment caused population decline. But adaptive evolution quickly rescued the population, with survival rates climbing to 69% by the F5 generation and 90% by the F10 generation. Population genomic data showed that rescue likely was caused by rapid evolutionary change at multiple loci, with many alleles fixing or nearly fixing within five generations of selection on lentil. Selection on these loci was only moderately consistent in time, but parallel evolutionary changes were evident in sublines formed after the lentil line had passed through a bottleneck. By comparing estimates of selection and genomic change on lentil across five independent C. maculatus lines (the new lentil‐adapted line, three long‐established lines and one case of failed evolutionary rescue), we found that adaptation on lentil occurred via somewhat idiosyncratic evolutionary changes. Overall, our results suggest that evolutionary rescue in this system can be caused by very strong selection on multiple loci driving rapid and pronounced genomic change.     

Identification of a genetic network for an ecologically relevant behavioural phenotype in Drosophila melanogaster

Pupation site choice of Drosophila third‐instar larvae is critical for the survival of individuals, as pupae are exposed to various biotic and abiotic dangers while immobilized during the 3–4 days of metamorphosis. This singular behavioural choice is sensitive to both environmental and genetic factors. Here, we developed a high‐throughput phenotyping approach to assay the variation in pupation height in Drosophila melanogaster, while controlling for possibly confounding factors. We find substantial variation of mean pupation height among sampled natural stocks and we show that the Drosophila Genetic Reference Panel (DGRP) reflects this variation. Using the DGRP stocks for genome‐wide association (GWA) mapping, 16 loci involved in determining pupation height could be resolved. The candidate genes in these loci are enriched for high expression in the larval central nervous system. A genetic network could be constructed from the candidate loci, which places scribble (scrib) at the centre, plus other genes known to be involved in nervous system development, such as Epidermal growth factor receptor (Egfr) and p53. Using gene disruption lines, we could functionally validate several of the initially identified loci, as well as additional loci predicted from network analysis. Our study shows that the combination of high‐throughput phenotyping with a genetic analysis of variation captured from the wild can be used to approach the genetic dissection of an environmentally relevant behavioural phenotype.     

Riverscape genomics of a threatened fish across a hydroclimatically heterogeneous river basin

Understanding how natural selection generates and maintains adaptive genetic diversity in heterogeneous environments is key to predicting the evolutionary response of populations to rapid environmental change. Detecting selection in complex spatial environments remains challenging, especially for threatened species where the effects of strong genetic drift may overwhelm signatures of selection. We carried out a basinwide riverscape genomic analysis in the threatened southern pygmy perch (Nannoperca australis), an ecological specialist with low dispersal potential. High‐resolution environmental data and 5162 high‐quality filtered SNPs were used to clarify spatial population structure and to assess footprints of selection associated with a steep hydroclimatic gradient and with human disturbance across the naturally and anthropogenically fragmented Murray–Darling Basin (Australia). Our approach included F_ST outlier tests to define neutral loci, and a combination of spatially explicit genotype–environment association analyses to identify candidate adaptive loci while controlling for the effects of landscape structure and shared population history. We found low levels of genetic diversity and strong neutral population structure consistent with expectations based on spatial stream hierarchy and life history. In contrast, variables related to precipitation and temperature appeared as the most important environmental surrogates for putatively adaptive genetic variation at both regional and local scales. Human disturbance also influenced the variation in candidate loci for adaptation, but only at a local scale. Our study contributes to understanding of adaptive evolution along naturally and anthropogenically fragmented ecosystems. It also offers a tangible example of the potential contributions of landscape genomics for informing in situ and ex situ conservation management of biodiversity.     

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host–pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1′665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen‐mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host–pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.     

Population genomics fits the bill: genetics of adaptive beak variation in Darwin's finches

Darwin's finches are an iconic case of adaptive radiation. The size and shape of their beaks are key adaptive traits related to trophic niche that vary among species and evolve rapidly when the food supply changes. Building on recent studies, a paper in this issue of Molecular Ecology (Chaves et al. 2016 ) investigates the genomic basis of beak size variation in sympatric populations of three species of ground finch (Geospiza) by performing a Genome‐wide association study using RAD‐seq data. The authors find that variation in a small number of markers can explain a substantial proportion of variation in beak size. Some of these markers are in genomic regions that have previously been implicated in beak size variation in Darwin's finches, whereas other markers have not, suggesting both conservation and divergence in the genetic basis of morphological evolution. Overall, the study confirms that loci of large effect are involved in beak size variation, which helps to explain the high heritability and rapid response to selection of this trait. The independent identification of regions containing HMGA2 and DLK1 loci in a GWAS makes them prime targets for functional studies. The study also shows that under the right conditions, RAD‐seq can be a viable alternative to genome sequencing for GWAS in wild vertebrate populations.