The genetic consequences of selection in natural populations

Evolutionary rescue in resistance to pesticides

Evolutionary rescue occurs when the genetic evolution of adaptation saves a population from decline or extinction after environmental change. The evolution of resistance to pesticides is a special scenario of abrupt environmental change, where rescue occurs under (very) strong selection for one or a few de novo resistance mutations of large effect. Here, a population genetic model of evolutionary rescue with density-dependent population change is developed, with a focus on deriving results that are important to resistance management. Massive stochastic simulations are used to generate observations, which are accurately predicted using analytical approximations. Key results include the probability density function for the time to resistance and the probability of population extinction. The distribution of resistance times shows a lag period, a narrow peak and a long tail. Surprisingly, the mean time to resistance can increase with the strength of selection because, if a mutation does not occur early on, then its emergence is delayed by the pesticide reducing the population size. The probability of population extinction shows a sharp transition, in that when extinction is possible, it is also highly likely. Consequently, population suppression and (local) eradication can be theoretically achievable goals, as novel strategies to delay resistance evolution.

Cross-pollination in seed-blended refuge and selection for Vip3A resistance in a lepidopteran pest as detected by genomic monitoring

The evolution of pest resistance to management tools reduces productivity and results in economic losses in agricultural systems. To slow its emergence and spread, monitoring and prevention practices are implemented in resistance management programs. Recent work suggests that genomic approaches can identify signs of emerging resistance to aid in resistance management. Here, we empirically examined the sensitivity of genomic monitoring for resistance management in transgenic Bt crops, a globally important agricultural innovation. Whole genome resequencing of wild North American Helicoverpa zea collected from non-expressing refuge and plants expressing Cry1Ab confirmed that resistance-associated signatures of selection were detectable after a single generation of exposure. Upon demonstrating its sensitivity, we applied genomic monitoring to wild H. zea that survived Vip3A exposure resulting from cross-pollination of refuge plants in seed-blended plots. Refuge seed interplanted with transgenic seed exposed H. zea to sublethal doses of Vip3A protein in corn ears and was associated with allele frequency divergence across the genome. Some of the greatest allele frequency divergence occurred in genomic regions adjacent to a previously described candidate gene for Vip3A resistance. Our work highlights the power of genomic monitoring to sensitively detect heritable changes associated with field exposure to Bt toxins and suggests that seed-blended refuge will likely hasten the evolution of resistance to Vip3A in lepidopteran pests.

Adaptation to distinct habitats is maintained by contrasting selection at different life stages in sunflower ecotypes

Conspecific populations living in adjacent but contrasting microenvironments represent excellent systems for studying natural selection. These systems are valuable because gene flow is expected to force genetic homogeneity except at loci experiencing divergent selection. A history of reciprocal transplant and common garden studies in such systems, and a growing number of genomic studies, have contributed to understanding how selection operates in natural populations. While selection can vary across different fitness components and life stages, few studies have investigated how this ultimately affects allele frequencies and the maintenance of divergence between populations. Here, we study two sunflower ecotypes in distinct, adjacent habitats by combining demographic models with genome-wide sequence data to estimate fitness and allele frequency change at multiple life stages. This framework allows us to estimate that only local ecotypes are likely to experience positive population growth (λ > 1) and that the maintenance of divergent adaptation appears to be mediated via habitat- and life stage-specific selection. We identify genetic variation, significantly driven by loci in chromosomal inversions, associated with different life history strategies in neighbouring ecotypes that optimize different fitness components and may contribute to the maintenance of distinct ecotypes.

Large haploblocks underlie rapid adaptation in the invasive weed Ambrosia artemisiifolia

Adaptation is the central feature and leading explanation for the evolutionary diversification of life. Adaptation is also notoriously difficult to study in nature, owing to its complexity and logistically prohibitive timescale. Here, we leverage extensive contemporary and historical collections of Ambrosia artemisiifolia-an aggressively invasive weed and primary cause of pollen-induced hayfever-to track the phenotypic and genetic causes of recent local adaptation across its native and invasive ranges in North America and Europe, respectively. Large haploblocks-indicative of chromosomal inversions-contain a disproportionate share (26%) of genomic regions conferring parallel adaptation to local climates between ranges, are associated with rapidly adapting traits, and exhibit dramatic frequency shifts over space and time. These results highlight the importance of large-effect standing variants in rapid adaptation, which have been critical to A. artemisiifolia's global spread across vast climatic gradients.

Allele frequencies and selection coefficients in locally adapted populations

Understanding the role of natural selection in driving evolutionary change requires accurate estimates of the strength of selection acting at the genetic level in the wild. This is challenging to achieve but may be easier in the case of populations in migration-selection balance. When two populations are at equilibrium under migration-selection balance, there exist loci whose alleles are selected different ways in the two populations. Such loci can be identified from genome sequencing by their high values of F_ST. This raises the question of what is the strength of selection on locally-adaptive alleles. To answer this question we analyse a 1-locus 2-allele model of a population distributed between two niches. We show by simulation of selected cases that the outputs from finite-population models are essentially the same as those from deterministic infinite-population models. We then derive theory for the infinite-population model showing the dependence of selection coefficients on equilibrium allele frequencies, migration rates, dominance and relative population sizes in the two niches. An Excel spreadsheet is provided for the calculation of selection coefficients and their approximate standard errors from observed values of population parameters. We illustrate our results with a worked example, with graphs showing the dependence of selection coefficients on equilibrium allele frequencies, and graphs showing how F_ST depends on the selection coefficients acting on the alleles at a locus. Given the extent of recent progress in ecological genomics, we hope our methods may help those studying migration-selection balance to quantify the advantages conferred by adaptive genes.

Herbaria reveal cost of the Green Revolution

Rapid weed evolution is exposed by genome sequencing of natural history collections.

Utility and challenges of using whole‐genome resequencing to detect emerging insect and mite resistance in agroecosystems

Arthropods that invade agricultural ecosystems systematically evolve resistance to the control measures used against them, and this remains a significant and ongoing challenge for sustainable food production systems. Early detection of resistance evolution could prompt remedial action to slow the spread of resistance alleles in the landscape. Historical approaches used to detect emerging resistance included phenotypic monitoring of agricultural pest populations, as well as monitoring of allele frequency changes at one or a few candidate pesticide resistance genes. In this article, I discuss the successes and limitations of these traditional monitoring approaches and then consider whether whole‐genome scanning could be applied to samples collected from agroecosystems over time for resistance monitoring. I examine the qualities of agroecosystems that could impact application of this approach to pesticide resistance monitoring and describe a recent retrospective analysis where genome scanning successfully detected an oligogenic response to selection by pesticides years prior to pest management failure. I conclude by considering areas of further study that will shed light on the feasibility of applying whole‐genome scanning for resistance risk monitoring in agricultural pest species.

The short-term, genome-wide effects of indirect selection deserve study: A response to Charlesworth and Jensen (2022)

We recently published a paper quantifying the genome-wide consequences of natural selection, including the effects of indirect selection due to the correlation of genetic regions (neutral or selected) with directly selected regions (Gompert et al., 2022). In their critique of our paper, Charlesworth and Jensen (2022) make two main points: (i) indirect selection is equivalent to hitchhiking and thus well documented (i.e., our results are not novel) and (ii) that we do not demonstrate the source of linkage disequilibrium (LD) between SNPs and the Mel-Stripe locus in the Timema cristinae experiment we analyse. As we discuss in detail below, neither of these are substantial criticisms of our work.

Freshwater Colonization, Adaptation, and Genomic Divergence in Threespine Stickleback

The Threespine Stickleback is ancestrally a marine fish, but many marine populations breed in fresh water (i.e., are anadromous), facilitating their colonization of isolated freshwater habitats a few years after they form. Repeated adaptation to fresh water during at least 10 My and continuing today has led to Threespine Stickleback becoming a premier system to study rapid adaptation. Anadromous and freshwater stickleback breed in sympatry and may hybridize, resulting in introgression of freshwater-adaptive alleles into anadromous populations, where they are maintained at low frequencies as ancient standing genetic variation. Anadromous stickleback have accumulated hundreds of freshwater-adaptive alleles that are disbursed as few loci per marine individual and provide the basis for adaptation when they colonize fresh water. Recent whole-lake experiments in lakes around Cook Inlet, Alaska have revealed how astonishingly rapid and repeatable this process is, with the frequency of 40% of the identified freshwater-adaptive alleles increasing from negligible (∼1%) in the marine founder to ≥50% within ten generations in fresh water, and freshwater phenotypes evolving accordingly. These high rates of genomic and phenotypic evolution imply very intense directional selection on phenotypes of heterozygotes. Sexual recombination rapidly assembles freshwater-adaptive alleles that originated in different founders into multilocus freshwater haplotypes, and regions important for adaptation to freshwater have suppressed recombination that keeps advantageous alleles linked within large haploblocks. These large haploblocks are also older and appear to have accumulated linked advantageous mutations. The contemporary evolution of Threespine Stickleback has provided broadly applicable insights into the mechanisms that facilitate rapid adaptation.

Faster-haplodiploid evolution under divergence-with-gene-flow: simulations and empirical data from pine-feeding hymenopterans

Although haplodiploidy is widespread in nature, the evolutionary consequences of this mode of reproduction are not well characterized. Here, we examine how genome-wide hemizygosity and a lack of recombination in haploid males affects genomic differentiation in populations that diverge via natural selection while experiencing gene flow. First, we simulated diploid and haplodiploid "genomes" (500-kb loci) evolving under an isolation-with-migration model with mutation, drift, selection, migration, and recombination; and examined differentiation at neutral sites both tightly and loosely linked to a divergently selected site. So long as there is divergent selection and migration, sex-limited hemizygosity and recombination cause elevated differentiation (i.e., produce a "faster-haplodiploid effect") in haplodiploid populations relative to otherwise equivalent diploid populations, for both recessive and codominant mutations. Second, we used genome-wide SNP data to model divergence history and describe patterns of genomic differentiation between sympatric populations of Neodiprion lecontei and N. pinetum, a pair of pine sawfly species (order: Hymenoptera; family: Diprionidae) that are specialized on different pine hosts. These analyses support a history of continuous gene exchange throughout divergence and reveal a pattern of heterogeneous genomic differentiation that is consistent with divergent selection on many unlinked loci. Third, using simulations of haplodiploid and diploid populations evolving according to the estimated divergence history of N. lecontei and N. pinetum, we found that divergent selection would lead to higher differentiation in haplodiploids. Based on these results, we hypothesize that haplodiploids undergo divergence-with-gene-flow and sympatric speciation more readily than diploids.

Why are deleterious mutations maintained in selfing populations? An analysis of the effects of early- and late-acting mutations by a two-locus two-allele model

Frequencies of deleterious mutations are higher than expected in many plants. Here, by developing a two-locus two-allele model, I examine the effects of differential timing of the expression of deleterious mutations (two-stage effects) on the maintenance of mutations. I assume early- and late-acting loci to distinguish whether maintenance of mutations in populations with high selfing rates is explained better by two-stage effects of single mutations, or by separate mutations in both early- and late-acting loci. I found that, when ovules are overproduced, the stable frequency of early-acting mutations is higher if mutations also occur in a late-acting locus than if a late-acting mutation is lacking. The stable frequency of late-acting mutations is higher if mutations also occur in an early-acting locus than if an early-acting mutation is lacking. Selective interference does not account for these results because analyses in which the number of loci subject to mutations is equalized are included. Overproduction of ovules has little effect on maintenance if either early- or late-acting mutations are lacking, whereas when ovules are not overproduced, the two-stage effect does not enhance the maintenance of mutations. Hence, mutations occurring in both loci coupled with overproduction of ovules enhances the maintenance of mutations in populations with high selfing rates. The detailed mechanisms underlying the two-stage effect were also analyzed.

Natural selection drives genome-wide evolution via chance genetic associations

Understanding selection's impact on the genome is a major theme in biology. Functionally neutral genetic regions can be affected indirectly by natural selection, via their statistical association with genes under direct selection. The genomic extent of such indirect selection, particularly across loci not physically linked to those under direct selection, remains poorly understood, as does the time scale at which indirect selection occurs. Here, we use field experiments and genomic data in stick insects, deer mice and stickleback fish to show that widespread statistical associations with genes known to affect fitness cause many genetic loci across the genome to be impacted indirectly by selection. This includes regions physically distant from those directly under selection. Then, focusing on the stick insect system, we show that statistical associations between SNPs and other unknown, causal variants result in additional indirect selection in general and specifically within genomic regions of physically linked loci. This widespread indirect selection necessarily makes aspects of evolution more predictable. Thus, natural selection combines with chance genetic associations to affect genome-wide evolution across linked and unlinked loci and even in modest-sized populations. This process has implications for the application of evolutionary principles in basic and applied science.

Reviewing the consequences of genetic purging on the success of rescue programs

Genetic rescue is increasingly considered a promising and underused conservation strategy to reduce inbreeding depression and restore genetic diversity in endangered populations, but the empirical evidence supporting its application is limited to a few generations. Here we discuss on the light of theory the role of inbreeding depression arising from partially recessive deleterious mutations and of genetic purging as main determinants of the medium to long-term success of rescue programs. This role depends on two main predictions: (1) The inbreeding load hidden in populations with a long stable demography increases with the effective population size; and (2) After a population shrinks, purging tends to remove its (partially) recessive deleterious alleles, a process that is slower but more efficient for large populations than for small ones. We also carry out computer simulations to investigate the impact of genetic purging on the medium to long term success of genetic rescue programs. For some scenarios, it is found that hybrid vigor followed by purging will lead to sustained successful rescue. However, there may be specific situations where the recipient population is so small that it cannot purge the inbreeding load introduced by migrants, which would lead to increased fitness inbreeding depression and extinction risk in the medium to long term. In such cases, the risk is expected to be higher if migrants came from a large non-purged population with high inbreeding load, particularly after the accumulation of the stochastic effects ascribed to repeated occasional migration events. Therefore, under the specific deleterious recessive mutation model considered, we conclude that additional caution should be taken in rescue programs. Unless the endangered population harbors some distinctive genetic singularity whose conservation is a main concern, restoration by continuous stable gene flow should be considered, whenever feasible, as it reduces the extinction risk compared to repeated occasional migration and can also allow recolonization events.

Long-term exhaustion of the inbreeding load in Drosophila melanogaster

Inbreeding depression, the decline in fitness of inbred individuals, is a ubiquitous phenomenon of great relevance in evolutionary biology and in the fields of animal and plant breeding and conservation. Inbreeding depression is due to the expression of recessive deleterious alleles that are concealed in heterozygous state in noninbred individuals, the so-called inbreeding load. Genetic purging reduces inbreeding depression by removing these alleles when expressed in homozygosis due to inbreeding. It is generally thought that fast inbreeding (such as that generated by full-sib mating lines) removes only highly deleterious recessive alleles, while slow inbreeding can also remove mildly deleterious ones. However, a question remains regarding which proportion of the inbreeding load can be removed by purging under slow inbreeding in moderately large populations. We report results of two long-term slow inbreeding Drosophila experiments (125-234 generations), each using a large population and a number of derived lines with effective sizes about 1000 and 50, respectively. The inbreeding load was virtually exhausted after more than one hundred generations in large populations and between a few tens and over one hundred generations in the lines. This result is not expected from genetic drift alone, and is in agreement with the theoretical purging predictions. Computer simulations suggest that these results are consistent with a model of relatively few deleterious mutations of large homozygous effects and partially recessive gene action.

Evolution of resistance under alternative models of selective interference

The use of multiple pesticides or drugs can lead to a simultaneous selection pressure for resistance alleles at different loci. Models of resistance evolution focus on how this can delay the spread of resistance through a population, but often neglect how this can also reduce the probability that a resistance allele spreads. This neglected factor has been studied in a parallel literature as selective interference. Models of interference use alternative constructions of fitness, where selection coefficients from different loci either add or multiply. Although these are equivalent under weak selection, the two constructions make alternative predictions under the strong selection that characterizes resistance evolution. Here, simulations are used to examine the effects of interference on the probability of fixation and time to fixation of a new and strongly beneficial mutation in the presence of another strongly beneficial allele with variable starting frequency. The results from simulations show a complicated pattern of effects. The key result is that, under multiplicativity, the presence of the strongly beneficial allele leads to a small reduction in the probability of fixation for the new beneficial mutation up to ~10%, and a negligible increase in the average time to fixation up to ~2%, whereas under additivity, the effect is more substantial at up to ~50% for the probability of fixation and ~100% for the average time to fixation. Consequently, the effect of interference is only an important feature of resistance evolution under additivity. Current evidence from studies of experimental evolution provides widespread support for the basic features of additivity, which suggests that interference may afford resistance a different pattern of evolution than other adaptations: rather than the gradual and simultaneous selection of many alleles with small effects, the rapid evolution of resistance may involve the sequential selection of alleles with large effects.

The Antimicrobial Resistance Crisis: An Inadvertent, Unfortunate but Nevertheless Informative Experiment in Evolutionary Biology

The global rise of antimicrobial resistance (AMR) phenotypes is an exemplar for rapid evolutionary response. Resistance arises as a consequence of humanity’s widespread and largely indiscriminate use of antimicrobial compounds. However, some features of this crisis remain perplexing. The remarkably widespread and rapid rise of diverse, novel and effective resistance phenotypes is in stark contrast to the apparent paucity of antimicrobial producers in the global microbiota. From the viewpoint of evolutionary theory, it should be possible to use selection coefficients to examine these phenomena. In this work we introduce an elaboration on the selection coefficient s termed selective efficiency, considering the genetic, metabolic, ecological and evolutionary impacts that accompany selective phenotypes. We then demonstrate the utility of the selective efficiency concept using AMR and antimicrobial production phenotypes as ‘worked examples’ of the concept. In accomplishing this objective, we also put forward cogent hypotheses to explain currently puzzling aspects of the AMR crisis. Finally, we extend the selective efficiency concept into a consideration of the ongoing management of the AMR crisis.

Predicting future from past: The genomic basis of recurrent and rapid stickleback evolution

Similar forms often evolve repeatedly in nature, raising long-standing questions about the underlying mechanisms. Here, we use repeated evolution in stickleback to identify a large set of genomic loci that change recurrently during colonization of freshwater habitats by marine fish. The same loci used repeatedly in extant populations also show rapid allele frequency changes when new freshwater populations are experimentally established from marine ancestors. Marked genotypic and phenotypic changes arise within 5 years, facilitated by standing genetic variation and linkage between adaptive regions. Both the speed and location of changes can be predicted using empirical observations of recurrence in natural populations or fundamental genomic features like allelic age, recombination rates, density of divergent loci, and overlap with mapped traits. A composite model trained on these stickleback features can also predict the location of key evolutionary loci in Darwin's finches, suggesting that similar features are important for evolution across diverse taxa.

Preferred reporting items for systematic reviews and meta-analyses in ecology and evolutionary biology: a PRISMA extension

Since the early 1990s, ecologists and evolutionary biologists have aggregated primary research using meta-analytic methods to understand ecological and evolutionary phenomena. Meta-analyses can resolve long-standing disputes, dispel spurious claims, and generate new research questions. At their worst, however, meta-analysis publications are wolves in sheep's clothing: subjective with biased conclusions, hidden under coats of objective authority. Conclusions can be rendered unreliable by inappropriate statistical methods, problems with the methods used to select primary research, or problems within the primary research itself. Because of these risks, meta-analyses are increasingly conducted as part of systematic reviews, which use structured, transparent, and reproducible methods to collate and summarise evidence. For readers to determine whether the conclusions from a systematic review or meta-analysis should be trusted - and to be able to build upon the review - authors need to report what they did, why they did it, and what they found. Complete, transparent, and reproducible reporting is measured by 'reporting quality'. To assess perceptions and standards of reporting quality of systematic reviews and meta-analyses published in ecology and evolutionary biology, we surveyed 208 researchers with relevant experience (as authors, reviewers, or editors), and conducted detailed evaluations of 102 systematic review and meta-analysis papers published between 2010 and 2019. Reporting quality was far below optimal and approximately normally distributed. Measured reporting quality was lower than what the community perceived, particularly for the systematic review methods required to measure trustworthiness. The minority of assessed papers that referenced a guideline (~16%) showed substantially higher reporting quality than average, and surveyed researchers showed interest in using a reporting guideline to improve reporting quality. The leading guideline for improving reporting quality of systematic reviews is the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) statement. Here we unveil an extension of PRISMA to serve the meta-analysis community in ecology and evolutionary biology: PRISMA-EcoEvo (version 1.0). PRISMA-EcoEvo is a checklist of 27 main items that, when applicable, should be reported in systematic review and meta-analysis publications summarising primary research in ecology and evolutionary biology. In this explanation and elaboration document, we provide guidance for authors, reviewers, and editors, with explanations for each item on the checklist, including supplementary examples from published papers. Authors can consult this PRISMA-EcoEvo guideline both in the planning and writing stages of a systematic review and meta-analysis, to increase reporting quality of submitted manuscripts. Reviewers and editors can use the checklist to assess reporting quality in the manuscripts they review. Overall, PRISMA-EcoEvo is a resource for the ecology and evolutionary biology community to facilitate transparent and comprehensively reported systematic reviews and meta-analyses.

Fitness maps to a large-effect locus in introduced stickleback populations

Mutations of small effect underlie most adaptation to new environments, but beneficial variants with large fitness effects are expected to contribute under certain conditions. Genes and genomic regions having large effects on phenotypic differences between populations are known from numerous taxa, but fitness effect sizes have rarely been estimated. We mapped fitness over a generation in an F2 intercross between a marine and a lake stickleback population introduced to a freshwater pond. A quantitative trait locus map of the number of surviving offspring per F2 female detected a single, large-effect locus near Ectodysplasin (Eda), a gene having an ancient freshwater allele causing reduced bony armor and other changes. F2 females homozygous for the freshwater allele had twice the number of surviving offspring as homozygotes for the marine allele, producing a large selection coefficient, s = 0.50 ± 0.09 SE. Correspondingly, the frequency of the freshwater allele increased from 0.50 in F2 mothers to 0.58 in surviving offspring. We compare these results to allele frequency changes at the Eda gene in an Alaskan lake population colonized by marine stickleback in the 1980s. The frequency of the freshwater Eda allele rose steadily over multiple generations and reached 95% within 20 y, yielding a similar estimate of selection, s = 0.49 ± 0.05, but a different degree of dominance. These findings are consistent with other studies suggesting strong selection on this gene (and/or linked genes) in fresh water. Selection on ancient genetic variants carried by colonizing ancestors is likely to increase the prevalence of large-effect fitness variants in adaptive evolution.

Natural selection maintains species despite frequent hybridization in the desert shrub Encelia

Natural selection is an important driver of genetic and phenotypic differentiation between species. For species in which potential gene flow is high but realized gene flow is low, adaptation via natural selection may be a particularly important force maintaining species. For a recent radiation of New World desert shrubs (Encelia: Asteraceae), we use fine-scale geographic sampling and population genomics to determine patterns of gene flow across two hybrid zones formed between two independent pairs of species with parapatric distributions. After finding evidence for extremely strong selection at both hybrid zones, we use a combination of field experiments, high-resolution imaging, and physiological measurements to determine the ecological basis for selection at one of the hybrid zones. Our results identify multiple ecological mechanisms of selection (drought, salinity, herbivory, and burial) that together are sufficient to maintain species boundaries despite high rates of hybridization. Given that multiple pairs of Encelia species hybridize at ecologically divergent parapatric boundaries, such mechanisms may maintain species boundaries throughout Encelia.

Drift, selection and adaptive variation in small populations of a threatened rattlesnake

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long- and short-term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre-date the effects of recent drift, and that functional variation in these loci persists despite small short-term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of "drift debt," a nonequilibrium state where present-day levels of variation overestimate the amount of functional genetic diversity present in future populations.

Root and shoot variation in relation to potential intermittent drought adaptation of Mesoamerican wild common bean (Phaseolus vulgaris L.)

BACKGROUND

Wild crop relatives have been potentially subjected to stresses on an evolutionary time scale prior to domestication. Among these stresses, drought is one of the main factors limiting crop productivity and its impact is likely to increase under current scenarios of global climate change. We sought to determine to what extent wild common bean (Phaseolus vulgaris) exhibited adaptation to drought stress, whether this potential adaptation is dependent on the climatic conditions of the location of origin of individual populations, and to what extent domesticated common bean reflects potential drought adaptation.

METHODS

An extensive and diverse set of wild beans from across Mesoamerica, along with a set of reference Mesoamerican domesticated cultivars, were evaluated for root and shoot traits related to drought adaptation. A water deficit experiment was conducted by growing each genotype in a long transparent tube in greenhouse conditions so that root growth, in addition to shoot growth, could be monitored.

RESULTS

Phenotypic and landscape genomic analyses, based on single-nucleotide polymorphisms, suggested that beans originating from central and north-west Mexico and Oaxaca, in the driest parts of their distribution, produced more biomass and were deeper-rooted. Nevertheless, deeper rooting was correlated with less root biomass production relative to total biomass. Compared with wild types, domesticated types showed a stronger reduction and delay in growth and development in response to drought stress. Specific genomic regions were associated with root depth, biomass productivity and drought response, some of which showed signals of selection and were previously related to productivity and drought tolerance.

CONCLUSIONS

The drought tolerance of wild beans consists in its stronger ability, compared with domesticated types, to continue growth in spite of water-limited conditions. This study is the first to relate bean response to drought to environment of origin for a diverse selection of wild beans. It provides information that needs to be corroborated in crosses between wild and domesticated beans to make it applicable to breeding programmes.

Natural selection on the Arabidopsis thaliana genome in present and future climates

Through the lens of evolution, climate change is an agent of natural selection that forces populations to change and adapt, or face extinction. However, current assessments of the risk of biodiversity associated with climate change¹ do not typically take into account how natural selection influences populations differently depending on their genetic makeup². Here we make use of the extensive genome information that is available for Arabidopsis thaliana and measure how manipulation of the amount of rainfall affected the fitness of 517 natural Arabidopsis lines that were grown in Spain and Germany. This allowed us to directly infer selection along the genome³. Natural selection was particularly strong in the hot-dry location in Spain, where 63% of lines were killed and where natural selection substantially changed the frequency of approximately 5% of all genome-wide variants. A significant portion of this climate-driven natural selection of variants was predictable from signatures of local adaptation (R² = 29-52%), as genetic variants that were found in geographical areas with climates more similar to the experimental sites were positively selected. Field-validated predictions across the species range indicated that Mediterranean and western Siberian populations-at the edges of the environmental limits of this species-currently experience the strongest climate-driven selection. With more frequent droughts and rising temperatures in Europe⁴, we forecast an increase in directional natural selection moving northwards from the southern end of Europe, putting many native A. thaliana populations at evolutionary risk.

Understanding Maladaptation by Uniting Ecological and Evolutionary Perspectives

Evolutionary biologists have long trained their sights on adaptation, focusing on the power of natural selection to produce relative fitness advantages while often ignoring changes in absolute fitness. Ecologists generally have taken a different tack, focusing on changes in abundance and ranges that reflect absolute fitness while often ignoring relative fitness. Uniting these perspectives, we articulate various causes of relative and absolute maladaptation and review numerous examples of their occurrence. This review indicates that maladaptation is reasonably common from both perspectives, yet often in contrasting ways. That is, maladaptation can appear strong from a relative fitness perspective, yet populations can be growing in abundance. Conversely, resident individuals can appear locally adapted (relative to nonresident individuals) yet be declining in abundance. Understanding and interpreting these disconnects between relative and absolute maladaptation, as well as the cases of agreement, is increasingly critical in the face of accelerating human-mediated environmental change. We therefore present a framework for studying maladaptation, focusing in particular on the relationship between absolute and relative fitness, thereby drawing together evolutionary and ecological perspectives. The unification of these ecological and evolutionary perspectives has the potential to bring together previously disjunct research areas while addressing key conceptual issues and specific practical problems.

How diverse can rare species be on the margins of genera distribution?

Different genetic patterns have been demonstrated for narrowly distributed taxa, many of them linking rarity to evolutionary history. Quite a few species in young genera are endemics and have several populations that present low variability, sometimes attributed to geographical isolation or dispersion processes. Assessing the genetic diversity and structure of such species may be important for protecting them and understanding their diversification history. In this study, we used microsatellite markers and plastid sequences to characterize the levels of genetic variation and population structure of two endemic and restricted species that grow in isolated areas on the margin of the distribution of their respective genera. Plastid and nuclear diversities were very low and weakly structured in their populations. Evolutionary scenarios for both species are compatible with open-field expansions during the Pleistocene interglacial periods and genetic variability supports founder effects to explain diversification. At present, both species are suffering from habitat loss and changes in the environment can lead these species towards extinction.

Estimating the relative fitness of escaped farmed salmon offspring in the wild and modelling the consequences of invasion for wild populations

Throughout their native range, wild Atlantic salmon populations are threatened by hybridization and introgression with escapees from net-pen salmon aquaculture. Although domestic-wild hybrid offspring have shown reduced fitness in laboratory and field experiments, consequential impacts on population abundance and genetic integrity remain difficult to predict in the field, in part because the strength of selection against domestic offspring is often unknown and context-dependent. Here, we follow a single large escape event of farmed Atlantic salmon in southern Newfoundland and monitor changes in the in-river proportions of hybrids and feral individuals over time using genetically based hybrid identification. Over a three-year period following the escape, the overall proportion of wild parr increased consistently (total wild proportion of 71.6%, 75.1% and 87.5% each year, respectively), with subsequent declines in feral (genetically pure farmed individuals originating from escaped, farmed adults) and hybrid parr. We quantify the strength of selection against parr of aquaculture ancestry and explore the genetic and demographic consequences for populations in the region. Within-cohort changes in the relative proportions of feral and F1 parr suggest reduced relative survival compared to wild individuals over the first (0.15 and 0.81 for feral and F1, respectively) and second years of life (0.26, 0.83). These relative survivorship estimates were used to inform an individual-based salmon eco-genetic model to project changes in adult abundance and overall allele frequency across three invasion scenarios ranging from short-term to long-term invasion and three relative survival scenarios. Modelling results indicate that total population abundance and time to recovery were greatly affected by relative survivorship and predict significant declines in wild population abundance under continued large escape events and calculated survivorship. Overall, this work demonstrates the importance of estimating the strength of selection against domestic offspring in the wild to predict the long-term impact of farmed salmon escape events on wild populations.

Linking a mutation to survival in wild mice

Adaptive evolution in new or changing environments can be difficult to predict because the functional connections between genotype, phenotype, and fitness are complex. Here, we make these explicit connections by combining field and laboratory experiments in wild mice. We first directly estimate natural selection on pigmentation traits and an underlying pigment locus, Agouti, by using experimental enclosures of mice on different soil colors. Next, we show how a mutation in Agouti associated with survival causes lighter coat color through changes in its protein binding properties. Together, our findings demonstrate how a sequence variant alters phenotype and then reveal the ensuing ecological consequences that drive changes in population allele frequency, thereby illuminating the process of evolution by natural selection.

Evolution as an ecosystem process: insights from genomics

Evolution is a fundamental ecosystem process. The study of genomic variation of organisms can not only improve our understanding of evolutionary processes, but also of contemporary and future ecosystem dynamics. We argue that integrative research between the fields of genomics and ecosystem ecology could generate new insights. Specifically, studies of biodiversity and ecosystem functioning, evolutionary rescue, and eco-evolutionary dynamics could all benefit from information about variation in genome structure and the genetic architecture of traits, whereas genomic studies could benefit from information about the ecological context of evolutionary dynamics. We propose new ways to help link research on functional genomic diversity with (reciprocal) interactions between phenotypic evolution and ecosystem change. Despite numerous challenges, we anticipate that the wealth of genomic data being collected on natural populations will improve our understanding of ecosystems.

Seasonally fluctuating selection can maintain polymorphism at many loci via segregation lift

Most natural populations are affected by seasonal changes in temperature, rainfall, or resource availability. Seasonally fluctuating selection could potentially make a large contribution to maintaining genetic polymorphism in populations. However, previous theory suggests that the conditions for multilocus polymorphism are restrictive. Here, we explore a more general class of models with multilocus seasonally fluctuating selection in diploids. In these models, the multilocus genotype is mapped to fitness in two steps. The first mapping is additive across loci and accounts for the relative contributions of heterozygous and homozygous loci-that is, dominance. The second step uses a nonlinear fitness function to account for the strength of selection and epistasis. Using mathematical analysis and individual-based simulations, we show that stable polymorphism at many loci is possible if currently favored alleles are sufficiently dominant. This general mechanism, which we call "segregation lift," requires seasonal changes in dominance, a phenomenon that may arise naturally in situations with antagonistic pleiotropy and seasonal changes in the relative importance of traits for fitness. Segregation lift works best under diminishing-returns epistasis, is not affected by problems of genetic load, and is robust to differences in parameters across loci and seasons. Under segregation lift, loci can exhibit conspicuous seasonal allele-frequency fluctuations, but often fluctuations may be small and hard to detect. An important direction for future work is to formally test for segregation lift in empirical data and to quantify its contribution to maintaining genetic variation in natural populations.

Detecting positive selection in the genome

Population geneticists have long sought to understand the contribution of natural selection to molecular evolution. A variety of approaches have been proposed that use population genetics theory to quantify the rate and strength of positive selection acting in a species’ genome. In this review we discuss methods that use patterns of between-species nucleotide divergence and within-species diversity to estimate positive selection parameters from population genomic data. We also discuss recently proposed methods to detect positive selection from a population’s haplotype structure. The application of these tests has resulted in the detection of pervasive adaptive molecular evolution in multiple species.

Epigenetics in natural animal populations

Phenotypic plasticity is an important mechanism for populations to buffer themselves from environmental change. While it has long been appreciated that natural populations possess genetic variation in the extent of plasticity, a surge of recent evidence suggests that epigenetic variation could also play an important role in shaping phenotypic responses. Compared with genetic variation, epigenetic variation is more likely to have higher spontaneous rates of mutation and a more sensitive reaction to environmental inputs. In our review, we first provide an overview of recent studies on epigenetically encoded thermal plasticity in animals to illustrate environmentally-mediated epigenetic effects within and across generations. Second, we discuss the role of epigenetic effects during adaptation by exploring population epigenetics in natural animal populations. Finally, we evaluate the evolutionary potential of epigenetic variation depending on its autonomy from genetic variation and its transgenerational stability. Although many of the causal links between epigenetic variation and phenotypic plasticity remain elusive, new data has explored the role of epigenetic variation in facilitating evolution in natural populations. This recent progress in ecological epigenetics will be helpful for generating predictive models of the capacity of organisms to adapt to changing climates.

Integrating viability and fecundity selection to illuminate the adaptive nature of genetic clines

Genetically based trait variation across environmental gradients can reflect adaptation to local environments. However, natural populations that appear well-adapted often exhibit directional, not stabilizing, selection on ecologically relevant traits. Temporal variation in the direction of selection could lead to stabilizing selection across multiple episodes of selection, which might be overlooked in short-term studies that evaluate relationships of traits and fitness under only one set of conditions. Furthermore, nonrandom mortality prior to trait expression can bias inferences about trait evolution if viability selection opposes fecundity selection. Here, we leveraged fitness and trait data to test whether phenotypic clines are genetically based and adaptive, whether temporal variation in climate imposes stabilizing selection, and whether viability selection acts on adult phenotypes. We monitored transplants of the subalpine perennial forb, Boechera stricta (Brassicaceae), in common gardens at two elevations over 2-3 years that differed in drought intensity. We quantified viability, and fecundity fitness components for four heritable traits: specific leaf area, integrated water-use efficiency, height at first flower, and flowering phenology. Our results indicate that genetic clines are maintained by selection, but their expression is context dependent, as they do not emerge in all environments. Moreover, selection varied spatially and temporally. Stabilizing selection was most pronounced when we integrated data across years. Finally, viability selection prior to trait expression targeted adult phenotypes (age and size at flowering). Indeed, viability selection for delayed flowering opposed fecundity selection for accelerated flowering; this result demonstrates that neglecting to account for viability selection could lead to inaccurate conclusions that populations are maladapted. Our results suggest that reconciling clinal trait variation with selection requires data collected across multiple spatial scales, time frames, and life-history stages.

Host-selected mutations converging on a global regulator drive an adaptive leap towards symbiosis in bacteria

Host immune and physical barriers protect against pathogens but also impede the establishment of essential symbiotic partnerships. To reveal mechanisms by which beneficial organisms adapt to circumvent host defenses, we experimentally evolved ecologically distinct bioluminescent Vibrio fischeri by colonization and growth within the light organs of the squid Euprymna scolopes. Serial squid passaging of bacteria produced eight distinct mutations in the binK sensor kinase gene, which conferred an exceptional selective advantage that could be demonstrated through both empirical and theoretical analysis. Squid-adaptive binK alleles promoted colonization and immune evasion that were mediated by cell-associated matrices including symbiotic polysaccharide (Syp) and cellulose. binK variation also altered quorum sensing, raising the threshold for luminescence induction. Preexisting coordinated regulation of symbiosis traits by BinK presented an efficient solution where altered BinK function was the key to unlock multiple colonization barriers. These results identify a genetic basis for microbial adaptability and underscore the importance of hosts as selective agents that shape emergent symbiont populations.

Advances in ecological genomics in forest trees and applications to genetic resources conservation and breeding

Forest trees are an unparalleled group of organisms in their combined ecological, economic and societal importance. With widespread distributions, predominantly random mating systems and large population sizes, most tree species harbour extensive genetic variation both within and among populations. At the same time, demographic processes associated with Pleistocene climate oscillations and land-use change have affected contemporary range-wide diversity and may impinge on the potential for future adaptation. Understanding how these adaptive and neutral processes have shaped the genomes of trees species is therefore central to their management and conservation. As for many other taxa, the advent of high-throughput sequencing methods is expected to yield an understanding of the interplay between the genome and environment at a level of detail and depth not possible only a few years ago. An international conference entitled 'Genomics and Forest Tree Genetics' was held in May 2016, in Arcachon (France), and brought together forest geneticists with a wide range of research interests to disseminate recent efforts that leverage contemporary genomic tools to probe the population, quantitative and evolutionary genomics of trees. An important goal of the conference was to discuss how such data can be applied to both genome-enabled breeding and the conservation of forest genetic resources under land use and climate change. Here, we report discoveries presented at the meeting and discuss how the ecological genomic toolkit can be used to address both basic and applied questions in tree biology.

On the maintenance of genetic variation and adaptation to environmental change: considerations from population genomics in fishes

The first goal of this paper was to overview modern approaches to local adaptation, with a focus on the use of population genomics data to detect signals of natural selection in fishes. Several mechanisms are discussed that may enhance the maintenance of genetic variation and evolutionary potential, which have been overlooked and should be considered in future theoretical development and predictive models: the prevalence of soft sweeps, polygenic basis of adaptation, balancing selection and transient polymorphisms, parallel evolution, as well as epigenetic variation. Research on fish population genomics has provided ample evidence for local adaptation at the genome level. Pervasive adaptive evolution, however, seems to almost never involve the fixation of beneficial alleles. Instead, adaptation apparently proceeds most commonly by soft sweeps entailing shifts in frequencies of alleles being shared between differentially adapted populations. One obvious factor contributing to the maintenance of standing genetic variation in the face of selective pressures is that adaptive phenotypic traits are most often highly polygenic, and consequently the response to selection should derive mostly from allelic co-variances among causative loci rather than pronounced allele frequency changes. Balancing selection in its various forms may also play an important role in maintaining adaptive genetic variation and the evolutionary potential of species to cope with environmental change. A large body of literature on fishes also shows that repeated evolution of adaptive phenotypes is a ubiquitous evolutionary phenomenon that seems to occur most often via different genetic solutions, further adding to the potential options of species to cope with a changing environment. Moreover, a paradox is emerging from recent fish studies whereby populations of highly reduced effective population sizes and impoverished genetic diversity can apparently retain their adaptive potential in some circumstances. Although more empirical support is needed, several recent studies suggest that epigenetic variation could account for this apparent paradox. Therefore, epigenetic variation should be fully integrated with considerations pertaining to role of soft sweeps, polygenic and balancing selection, as well as repeated adaptation involving different genetic basis towards improving models predicting the evolutionary potential of species to cope with a changing world.

Genome scans for divergent selection in natural populations of the widespread hardwood species Eucalyptus grandis (Myrtaceae) using microsatellites

Identification of loci or genes under natural selection is important for both understanding the genetic basis of local adaptation and practical applications, and genome scans provide a powerful means for such identification purposes. In this study, genome-wide simple sequence repeats markers (SSRs) were used to scan for molecular footprints of divergent selection in Eucalyptus grandis, a hardwood species occurring widely in costal areas from 32° S to 16° S in Australia. High population diversity levels and weak population structure were detected with putatively neutral genomic SSRs. Using three FST outlier detection methods, a total of 58 outlying SSRs were collectively identified as loci under divergent selection against three non-correlated climatic variables, namely, mean annual temperature, isothermality and annual precipitation. Using a spatial analysis method, nine significant associations were revealed between FST outlier allele frequencies and climatic variables, involving seven alleles from five SSR loci. Of the five significant SSRs, two (EUCeSSR1044 and Embra394) contained alleles of putative genes with known functional importance for response to climatic factors. Our study presents critical information on the population diversity and structure of the important woody species E. grandis and provides insight into the adaptive responses of perennial trees to climatic variations.