Boom-bust population dynamics drive rapid genetic change.

Increasing environmental threats and more extreme environmental perturbations place species at risk of population declines, with associated loss of genetic diversity and evolutionary potential. While theory shows that rapid population declines can cause loss of genetic diversity, populations in some environments, like Australia's arid zone, are repeatedly subject to major population fluctuations yet persist and appear able to maintain genetic diversity. Here, we use repeated population sampling over 13 y and genotype-by-sequencing of 1903 individuals to investigate the genetic consequences of repeated population fluctuations in two small mammals in the Australian arid zone. The sandy inland mouse (Pseudomys hermannsburgensis) experiences marked boom-bust population dynamics in response to the highly variable desert environment. We show that heterozygosity levels declined, and population differentiation (FST) increased, during bust periods when populations became small and isolated, but that heterozygosity was rapidly restored during episodic population booms. In contrast, the lesser hairy-footed dunnart (Sminthopsis youngsoni), a desert marsupial that maintains relatively stable population sizes, showed no linear declines in heterozygosity. These results reveal two contrasting ways in which genetic diversity is maintained in highly variable environments. In one species, diversity is conserved through the maintenance of stable population sizes across time. In the other species, diversity is conserved through rapid genetic mixing during population booms that restores heterozygosity lost during population busts.

Changing patterns of genetic differentiation in the slender wild oat, Avena barbata.

The slender wild oat (Avena barbata) was widely studied in California using allozymes in the 1970s and interpreted as a case of ecotypic adaptation to contrasting moisture environments. However, common garden studies suggested that the moist-associated ("mesic") ecotype had high fitness in both moist and dry habitats, thus predicting an adaptive spread into areas occupied by the dry associated ("xeric") ecotype. To test this prediction, we revisited 100 populations of A. barbata that were screened genetically 40 y ago. As expected, mesic allozyme and morphological markers are much more common than in the 1970s. The less-fit xeric ecotype, while still widespread, has declined markedly in range and frequency. Genotyping by sequencing of modern populations reveals striking genetic uniformity within each of the two ecotypes. In recombinants between the two ecotypes, the mesic allele at a major fitness quantitative trait locus (QTL) shows a high frequency but so do many other genomic regions not identified as fitness QTL. Additional introduced genotypes are diverse and more widespread than in the past, and our results show that these have spread into the former range of the xeric ecotype to an even greater extent than the mesic ecotype has. While these results confirm the prediction of contemporary evolution from common gardens, they also suggest that much of the change has been driven by additional waves of introduced genotypes.

Evidence for continent-wide convergent evolution and stasis throughout 150 y of a biological invasion.

The extent to which evolution can rescue a species from extinction, or facilitate range expansion, depends critically on the rate, duration, and geographical extent of the evolutionary response to natural selection. Adaptive evolution can occur quickly, but the duration and geographical extent of contemporary evolution in natural systems remain poorly studied. This is particularly true for species with large geographical ranges and for timescales that lie between "long-term" field experiments and the fossil record. Here, we introduce the Virtual Common Garden (VCG) to investigate phenotypic evolution in natural history collections while controlling for phenotypic plasticity in response to local growing conditions. Reconstructing 150 y of evolution in Lythrum salicaria (purple loosestrife) as it invaded North America, we analyze phenology measurements of 3,429 herbarium records, reconstruct growing conditions from more than 12 million local temperature records, and validate predictions across three common gardens spanning 10° of latitude. We find that phenological clines have evolved repeatedly throughout the range, during the first century of evolution. Thereafter, the rate of microevolution stalls, recapitulating macroevolutionary stasis observed in the fossil record. Our study demonstrates that preserved specimens are a critical resource for investigating limits to evolution in natural populations. Our results show how natural selection and trade-offs measured in field studies predict adaptive divergence observable in herbarium specimens over 15 decades at a continental scale.

Environmental change unlinks a gene from its trait.

In Atlantic salmon (Salmo salar), sea age is a major life history trait governed by a sex-specific trade-off between reproductive success and survival. In this issue of Molecular Ecology, Besnier et al. (Molecular Ecology, 2023) found evidence to suggest that the disassociation between sea age and major effect loci, including the previously identified candidate genes vgll3 and six6, may be related to the recently observed trend towards slower growth and later maturation. These results are of importance because they challenge the prevailing view that evolution moves in a slow shuffle, and they provide a pertinent example of how an optimal phenotype can change due to growth-driven plasticity and lead to contemporary molecular and phenotypic evolution.

Population genomic structure of a widespread, urban-dwelling mammal: The eastern grey squirrel (Sciurus carolinensis).

Urbanization is a persistent and widespread driver of global environmental change, potentially shaping evolutionary processes due to genetic drift and reduced gene flow in cities induced by habitat fragmentation and small population sizes. We tested this prediction for the eastern grey squirrel (Sciurus carolinensis), a common and conspicuous forest-dwelling rodent, by obtaining 44K SNPs using reduced representation sequencing (ddRAD) for 403 individuals sampled across the species' native range in eastern North America. We observed moderate levels of genetic diversity, low levels of inbreeding, and only a modest signal of isolation-by-distance. Clustering and migration analyses show that estimated levels of migration and genetic connectivity were higher than expected across cities and forested areas, specifically within the eastern portion of the species' range dominated by urbanization, and genetic connectivity was less than expected within the western range where the landscape is fragmented by agriculture. Landscape genetic methods revealed greater gene flow among individual squirrels in forested regions, which likely provide abundant food and shelter for squirrels. Although gene flow appears to be higher in areas with more tree cover, only slight discontinuities in gene flow suggest eastern grey squirrels have maintained connected populations across urban areas in all but the most heavily fragmented agricultural landscapes. Our results suggest urbanization shapes biological evolution in wildlife species depending strongly on the composition and habitability of the landscape matrix surrounding urban areas.

Spatial variation in genomic signatures of local adaptation during the cane toad invasion of Australia.

Adaptive evolution can facilitate species' range expansions across environmentally heterogeneous landscapes. However, serial founder effects can limit the efficacy of selection, and the evolution of increased dispersal during range expansions may result in gene flow swamping local adaptation. Here, we study how genetic drift, gene flow and selection interact during the cane toad's (Rhinella marina) invasion across the heterogeneous landscape of Australia. Following its introduction in 1935, the cane toad colonised eastern Australia and established several stable range edges. The ongoing, more rapid range expansion in north-central Australia has occurred concomitant with an evolved increase in dispersal capacity. Using reduced representation genomic data of Australian cane toads from the expansion front and from two areas of their established range, we test the hypothesis that high gene flow constrains local adaptation at the expansion front relative to established areas. Genetic analyses indicate the three study areas are genetically distinct but show similar levels of allelic richness, heterozygosity and inbreeding. Markedly higher gene flow or recency of colonisation at the expansion front have likely hindered local adaptation at the time of sampling, as indicated by reduced slopes of genetic-environment associations (GEAs) estimated using a novel application of geographically weighted regression that accounts for allele surfing; GEA slopes are significantly steeper in established parts of the range. Our work bolsters evidence supporting adaptation of invasive species post-introduction and adds novel evidence for differing strengths of evolutionary forces among geographic areas with different invasion histories.

Temporal genomics help in deciphering neutral and adaptive patterns in the contemporary evolution of kelp populations.

The impact of climate change on populations will be contingent upon their contemporary adaptive evolution. In this study, we investigated the contemporary evolution of 4 populations of the cold-water kelp Laminaria digitata by analyzing their spatial and temporal genomic variations using ddRAD-sequencing. These populations were sampled from the center to the southern margin of its north-eastern Atlantic distribution at 2 time points, spanning at least 2 generations. Through genome scans for local adaptation at a single time point, we identified candidate loci that showed clinal variation correlated with changes in sea surface temperature (SST) along latitudinal gradients. This finding suggests that SST may drive the adaptive response of these kelp populations, although factors such as species' demographic history should also be considered. Additionally, we performed a simulation approach to distinguish the effect of selection from genetic drift in allele frequency changes over time. This enabled the detection of loci in the southernmost population that exhibited temporal differentiation beyond what would be expected from genetic drift alone: these are candidate loci which could have evolved under selection over time. In contrast, we did not detect any outlier locus based on temporal differentiation in the population from the North Sea, which also displayed low and decreasing levels of genetic diversity. The diverse evolutionary scenarios observed among populations can be attributed to variations in the prevalence of selection relative to genetic drift across different environments. Therefore, our study highlights the potential of temporal genomics to offer valuable insights into the contemporary evolution of marine foundation species facing climate change.

Evolution of butterfly seasonal plasticity driven by climate change varies across life stages.

Photoperiod is a common cue for seasonal plasticity and phenology, but climate change can create cue-environment mismatches for organisms that rely on it. Evolution could potentially correct these mismatches, but phenology often depends on multiple plastic decisions made during different life stages and seasons that may evolve separately. For example, Pararge aegeria (Speckled wood butterfly) has photoperiod-cued seasonal life history plasticity in two different life stages: larval development time and pupal diapause. We tested for climate change-associated evolution of this plasticity by replicating common garden experiments conducted on two Swedish populations 30 years ago. We found evidence for evolutionary change in the contemporary larval reaction norm-although these changes differed between populations-but no evidence for evolution of the pupal reaction norm. This variation in evolution across life stages demonstrates the need to consider how climate change affects the whole life cycle to understand its impacts on phenology.

Eco-evolution from deep time to contemporary dynamics: The role of timescales and rate modulators.

Eco-evolutionary dynamics, or eco-evolution for short, are often thought to involve rapid demography (ecology) and equally rapid heritable phenotypic changes (evolution) leading to novel, emergent system behaviours. We argue that this focus on contemporary dynamics is too narrow: Eco-evolution should be extended, first, beyond pure demography to include all environmental dimensions and, second, to include slow eco-evolution which unfolds over thousands or millions of years. This extension allows us to conceptualise biological systems as occupying a two-dimensional time space along axes that capture the speed of ecology and evolution. Using Hutchinson's analogy: Time is the 'theatre' in which ecology and evolution are two interacting 'players'. Eco-evolutionary systems are therefore dynamic: We identify modulators of ecological and evolutionary rates, like temperature or sensitivity to mutation, which can change the speed of ecology and evolution, and hence impact eco-evolution. Environmental change may synchronise the speed of ecology and evolution via these rate modulators, increasing the occurrence of eco-evolution and emergent system behaviours. This represents substantial challenges for prediction, especially in the context of global change. Our perspective attempts to integrate ecology and evolution across disciplines, from gene-regulatory networks to geomorphology and across timescales, from today to deep time.

Resurrecting urban sunflowers: Phenotypic and molecular changes between antecedent and modern populations separated by 36 years.

Resurrection experiments provide a unique opportunity to evaluate phenotypic and molecular evolution in response to environmental challenges. To understand the evolution of urban populations of Helianthus annuus, we compared plants from 36-year-old antecedent seed collections to modern seed collections from the same area using molecular and quantitative genetic approaches. We found 200 differentially expressed transcripts between antecedent and modern groups, and transcript expression was generally higher in modern samples as compared to antecedent samples. Admixture analysis indicated gene flow from domesticated to modern populations over time. After a greenhouse refresher generation, one antecedent-modern population pair was grown under two water availability (well-watered and drought) and temperature (ambient and elevated by 2.8°C) conditions reflecting historical and contemporary climates. Overall, 78% (7 out of 9) of traits differed between the antecedent and modern populations, with modern individuals displaying some trait changes that are coherent with climate changes expectations and some trait changes in the direction of crop varieties. Phenotypic selection analysis showed that modern trait values were often favoured by selection, especially in environmental treatments resembling modern conditions. Trait heritability in the antecedent population was five times as high as in the modern population, on average. In addition, phenotypic plasticity for some traits, such as flowering phenology, was present in the antecedent population but absent in the modern population. The combination of phenotypic and molecular information suggests that evolution has been influenced by crop-wild introgression, adaptive processes and drift. We discuss these results in the context of continued evolution in response to anthropogenic factors.

Rapid evolution of novel biotic interactions in the UK Brown Argus butterfly uses genomic variation from across its geographical range.

Understanding the rate and extent to which populations can adapt to novel environments at their ecological margins is fundamental to predicting the persistence of biological communities during ongoing and rapid global change. Recent range expansion in response to climate change in the UK butterfly Aricia agestis is associated with the evolution of novel interactions with a larval food plant, and the loss of its ability to use an ancestral host species. Using ddRAD analysis of 61,210 variable SNPs from 261 females from throughout the UK range of this species, we identify genomic regions at multiple chromosomes that are associated with evolutionary responses, and their association with demographic history and ecological variation. Gene flow appears widespread throughout the range, despite the apparently fragmented nature of the habitats used by this species. Patterns of haplotype variation between selected and neutral genomic regions suggest that evolution associated with climate adaptation is polygenic, resulting from the independent spread of alleles throughout the established range of this species, rather than the colonization of pre-adapted genotypes from coastal populations. These data suggest that rapid responses to climate change do not depend on the availability of pre-adapted genotypes. Instead, the evolution of novel forms of biotic interaction in A. agestis has occurred during range expansion, through the assembly of novel genotypes from alleles from multiple localities.

Do chromosome rearrangements fix by genetic drift or natural selection? Insights from Brenthis butterflies.

Large-scale chromosome rearrangements, such as fissions and fusions, are a common feature of eukaryote evolution. They can have considerable influence on the evolution of populations, yet it remains unclear exactly how rearrangements become established and eventually fix. Rearrangements could fix by genetic drift if they are weakly deleterious or neutral, or they may instead be favoured by positive natural selection. Here, we compare genome assemblies of three closely related Brenthis butterfly species and characterize a complex history of fission and fusion rearrangements. An inferred demographic history of these species suggests that rearrangements became fixed in populations with large long-term effective size (Ne ), consistent with rearrangements being selectively neutral or only very weakly underdominant. Using a recently developed analytic framework for characterizing hard selective sweeps, we find that chromosome fusions are not enriched for evidence of past sweeps compared to other regions of the genome. Nonetheless, we do infer a strong and recent selective sweep around one chromosome fusion in the B. daphne genome. Our results suggest that rearrangements in these species likely have weak absolute fitness effects and fix by genetic drift. However, one putative selective sweep raises the possibility that natural selection may sometimes play a role in the fixation of chromosome fusions.

Keep your cool: Overwintering physiology in response to urbanization in the acorn ant, Temnothorax curvispinosus.

Winter presents a challenge for survival, yet temperate ectotherms have remarkable physiological adaptations to cope with low-temperature conditions. Under recent climate change, rather than strictly relaxing pressure on overwintering survival, warmer winters can instead disrupt these low-temperature trait-environment associations, with negative consequences for populations. While there is increasing evidence of physiological adaptation to contemporary warming during the growing season, the effects of winter warming on physiological traits are less clear. To address this knowledge gap, we performed a common garden experiment using relatively warm-adapted versus cold-adapted populations of the acorn ant, Temnothorax curvispinosus, sampled across an urban heat island gradient, to explore the effects of winter conditions on plasticity and evolution of physiological traits. We found no evidence of evolutionary divergence in chill coma recovery nor in metabolic rate at either of two test temperatures (4 and 10 °C). Although we found the expected plastic response of increased metabolic rate under the 10 °C acute test temperature as compared with the 4 °C test temperature, this plastic response, (i.e., the acute thermal sensitivity of metabolic rate), was not different across populations. Surprisingly, we found that winter-acclimated urban ant populations exhibited higher heat tolerance compared with rural ant populations, and that the magnitude of divergence was comparable to that observed among growing-season acclimated ants. Finally, we found no evidence of differences between populations with respect to changes in colony size from the beginning to the end of the overwintering experiment. Together, these findings indicate that despite the evolution of higher heat tolerance that is often accompanied by losses in low-temperature tolerance, urban acorn ants have retained several components of low-temperature physiological performance when assessed under ecologically relevant overwintering conditions. Our study suggests the importance of measuring physiological traits under seasonally-relevant conditions to understand the causes and consequences of evolutionary responses to contemporary warming.

Drivers and Cascading Ecological Consequences of Gambusia affinis Trait Variation.

AbstractPhenotypic trait differences among populations can shape ecological outcomes for communities and ecosystems. However, few studies have mechanistically linked heritable and plastic components of trait variation to generalizable processes of ecology, such as trophic cascades. Here, we assess morphological and behavioral trait variation in nine populations of common garden-reared western mosquitofish (Gambusia affinis) from three distinct ancestral predator environments (three populations per environment), each reared in the presence and absence of predator cues. We then use a pond mesocosm experiment to examine the ecological consequences of mosquitofish trait variation and density variation. Our results show significant among-population trait variation, but this variation was generally unrelated to ancestral predator environment. When traits did vary congruently with respect to ancestral predator environment, this trait variation was driven by gene-by-environment interactions. Variation in several mosquitofish traits altered the cascading effects of mosquitofish on zooplankton and primary producers, but the effect of any given trait was typically weaker than that of density. We note that the relatively stronger ecological effects of density may mask the effects of traits in some systems. Our example here shows that trait variation can be highly noncongruent with respect to a perceived selective agent, phenotypic change is a product of complex interactions between genes and the environment, and numerous interacting phenotypes generate significant but potentially cryptic cascading ecological change.

Inconsistent evolution and growth-survival tradeoffs in Gambusia affinis.

Growth-survival tradeoffs may be a generalizable mechanism influencing trajectories of prey evolution. Here, we investigate evolutionary contributions to growth and survival in western mosquitofish (Gambusia affinis) from 10 populations from high- and low-predation ancestral environments. We assess (i) the degree to which evolutionary components of growth and survival are consistent or inconsistent across populations within ancestral predation environments, and (ii) whether growth and survival trade off at the population level. We measure growth and survival on groups of common-reared mosquitofish in pond mesocosms. We find that evolution of growth is consistent, with fish from low-predation ancestral environments showing higher growth, while the evolution of survival is inconsistent, with significant population-level divergence unrelated to ancestral predation environment. Such inconsistency prevents a growth-survival tradeoff across populations. Thus, the generalizability of contemporary evolution probably depends on local context of evolutionary tradeoffs, and a continued focus on singular selective agents (e.g. predators) without such local context will impede insights into generalizable evolutionary patterns.

Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for antimalarial intervention deployment.

Studies of insecticide resistance provide insights into the capacity of populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on pyrethroid insecticides, primarily in long lasting insecticidal nets (LLINs). Resistance in the major malaria vectors has increased in concert with the expansion of LLIN distributions. Identifying genetic mechanisms underlying high-level resistance is crucial for the development and deployment of resistance-breaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000g) data we identified a very recent selective sweep in mosquitoes from Uganda which localized to a cluster of cytochrome P450 genes. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar-like transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border, with stepwise replacement of the double-mutant (Zanzibar-like TE and Cyp6p4-236 M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolize pyrethroids and is strongly predictive of resistance to pyrethroids most notably deltamethrin. Importantly, there was increased mortality in mosquitoes carrying the triple-mutation when exposed to nets cotreated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid cotreated LLINs across African countries.

Wild sunflower goes viral: Citizen science and comparative genomics allow tracking the origin and establishment of invasive sunflower in the Levant.

Globalization and intensified volume of trade and transport around the world are accelerating the rate of biological invasions. It is therefore increasingly important to understand the processes through which invasive species colonize new habitats, often to the detriment of native flora. The initial steps of an invasion are particularly critical, as the introduced species relies on limited genetic diversity to adapt to a new environment. However, our understanding of this critical stage of the invasion is currently limited. We used a citizen science approach and social media to survey the distribution of invasive sunflower in Israel. We then sampled and sequenced a representative collection and compared it with available genomic data sets of North American wild sunflower, landraces and cultivars. We show that invasive wild sunflower is rapidly establishing throughout Israel, probably from a single, recent introduction from Texas, while maintaining high genetic diversity through ongoing gene flow. Since its introduction, invasive sunflower has spread quickly to most regions, and differentiation was detected despite extensive gene flow between clusters. Our findings suggest that rapid spread followed by continuous gene flow between diverging populations can serve as an efficient mechanism for maintaining sufficient genetic diversity at the early stages of invasion, promoting rapid adaptation and establishment in the new territory.

The pace of modern life, revisited.

Wild populations must continuously respond to environmental changes or they risk extinction. Those responses can be measured as phenotypic rates of change, which can allow us to predict contemporary adaptive responses, some of which are evolutionary. About two decades ago, a database of phenotypic rates of change in wild populations was compiled. Since then, researchers have used (and expanded) this database to examine phenotypic responses to specific types of human disturbance. Here, we update the database by adding 5675 new estimates of phenotypic change. Using this newer version of the data base, now containing 7338 estimates of phenotypic change, we revisit the conclusions of four published articles. We then synthesize the expanded database to compare rates of change across different types of human disturbance. Analyses of this expanded database suggest that: (i) a small absolute difference in rates of change exists between human disturbed and natural populations, (ii) harvesting by humans results in higher rates of change than other types of disturbance, (iii) introduced populations have increased rates of change, and (iv) body size does not increase through time. Thus, findings from earlier analyses have largely held-up in analyses of our new database that encompass a much larger breadth of species, traits, and human disturbances. Lastly, we use new analyses to explore how various types of human disturbances affect rates of phenotypic change, and we call for this database to serve as a steppingstone for further analyses to understand patterns of contemporary phenotypic change.

Hybridization decreases native cutthroat trout reproductive fitness.

Examining natural selection in wild populations is challenging, but crucial to understanding many ecological and evolutionary processes. Additionally, in hybridizing populations, natural selection may be an important determinant of the eventual outcome of hybridization. We characterized several components of relative fitness in hybridizing populations of Yellowstone cutthroat trout and rainbow trout in an effort to better understand the prolonged persistence of both parental species despite predictions of extirpation. Thousands of genomic loci enabled precise quantification of hybrid status in adult and subsequent juvenile generations; a subset of those data also identified parent-offspring relationships. We used linear models and simulations to assess the effects of ancestry on reproductive output and mate choice decisions. We found a relatively low number of late-stage (F3+) hybrids and an excess of F2 juveniles relative to the adult generation in one location, which suggests the presence of hybrid breakdown decreasing the fitness of F2+ hybrids later in life. Assessments of reproductive output showed that Yellowstone cutthroat trout are more likely to successfully reproduce and produce slightly more offspring than their rainbow trout and hybrid counterparts. Mate choice appeared to be largely random, though we did find statistical support for slight female preference for males of similar ancestry. Together, these results show that native Yellowstone cutthroat trout are able to outperform rainbow trout in terms of reproduction and suggest that management action to exclude rainbow trout from spawning locations may bolster the now-rare Yellowstone cutthroat trout.

Rapid evolution generates synergism between multiple stressors: Linking theory and an evolution experiment.

Global change encompasses many co-occurring anthropogenic stressors. Understanding the interactions between these multiple stressors, whether they be additive, antagonistic or synergistic, is critical for ecosystem managers when prioritizing which stressors to mitigate in the face of global change. While such interactions between stressors appear prevalent, it remains unclear if and how these interactions change over time, as the majority of multiple-stressor studies rarely span multiple generations of study organisms. Although meta-analyses have reported some intriguing temporal trends in stressor interactions, for example that synergism may take time to emerge, the mechanistic basis for such observations is unknown. In this study, by analysing data from an evolution experiment with the rotifer Brachionus calyciflorus (~35 generations and 31,320 observations), we show that adaptation to multiple stressors shifts stressor interactions towards synergism. We show that trade-offs, where populations cannot optimally perform multiple tasks (i.e. adapting to multiple stressors), generate this bias towards synergism. We also show that removal of stressors from evolved populations does not necessarily increase fitness and that there is variation in the evolutionary trajectories of populations that experienced the same stressor regimes. Our results highlight outstanding questions at the interface between evolution and global change biology, and illustrate the importance of considering rapid adaptation when managing or restoring ecosystems subjected to multiple stressors under global change.