Competitive history shapes rapid evolution in a seasonal climate.

Eco-evolutionary dynamics will play a critical role in determining species' fates as climatic conditions change. Unfortunately, we have little understanding of how rapid evolutionary responses to climate play out when species are embedded in the competitive communities that they inhabit in nature. We tested the effects of rapid evolution in response to interspecific competition on subsequent ecological and evolutionary trajectories in a seasonally changing climate using a field-based evolution experiment with Drosophila melanogaster Populations of D. melanogaster were either exposed, or not exposed, to interspecific competition with an invasive competitor, Zaprionus indianus, over the summer. We then quantified these populations' ecological trajectories (abundances) and evolutionary trajectories (heritable phenotypic change) when exposed to a cooling fall climate. We found that competition with Z. indianus in the summer affected the subsequent evolutionary trajectory of D. melanogaster populations in the fall, after all interspecific competition had ceased. Specifically, flies with a history of interspecific competition evolved under fall conditions to be larger and have lower cold fecundity and faster development than flies without a history of interspecific competition. Surprisingly, this divergent fall evolutionary trajectory occurred in the absence of any detectible effect of the summer competitive environment on phenotypic evolution over the summer or population dynamics in the fall. This study demonstrates that competitive interactions can leave a legacy that shapes evolutionary responses to climate even after competition has ceased, and more broadly, that evolution in response to one selective pressure can fundamentally alter evolution in response to subsequent agents of selection.

Microbiome composition shapes rapid genomic adaptation of Drosophila melanogaster.

Population genomic data has revealed patterns of genetic variation associated with adaptation in many taxa. Yet understanding the adaptive process that drives such patterns is challenging; it requires disentangling the ecological agents of selection, determining the relevant timescales over which evolution occurs, and elucidating the genetic architecture of adaptation. Doing so for the adaptation of hosts to their microbiome is of particular interest with growing recognition of the importance and complexity of host-microbe interactions. Here, we track the pace and genomic architecture of adaptation to an experimental microbiome manipulation in replicate populations of Drosophila melanogaster in field mesocosms. Shifts in microbiome composition altered population dynamics and led to divergence between treatments in allele frequencies, with regions showing strong divergence found on all chromosomes. Moreover, at divergent loci previously associated with adaptation across natural populations, we found that the more common allele in fly populations experimentally enriched for a certain microbial group was also more common in natural populations with high relative abundance of that microbial group. These results suggest that microbiomes may be an agent of selection that shapes the pattern and process of adaptation and, more broadly, that variation in a single ecological factor within a complex environment can drive rapid, polygenic adaptation over short timescales.

The microbiota influences the Drosophila melanogaster life history strategy.

Organisms are locally adapted when members of a population have a fitness advantage in one location relative to conspecifics in other geographies. For example, across latitudinal gradients, some organisms may trade off between traits that maximize fitness components in one, but not both, of somatic maintenance or reproductive output. Latitudinal gradients in life history strategies are traditionally attributed to environmental selection on an animal's genotype, without any consideration of the possible impact of associated microorganisms ("microbiota") on life history traits. Here, we show in Drosophila melanogaster, a key model for studying local adaptation and life history strategy, that excluding the microbiota from definitions of local adaptation is a major shortfall. First, we reveal that an isogenic fly line reared with different bacteria varies the investment in early reproduction versus somatic maintenance. Next, we show that in wild fruit flies, the abundance of these same bacteria was correlated with the latitude and life history strategy of the flies, suggesting geographic specificity of the microbiota composition. Variation in microbiota composition of locally adapted D. melanogaster could be attributed to both the wild environment and host genetic selection. Finally, by eliminating or manipulating the microbiota of fly lines collected across a latitudinal gradient, we reveal that host genotype contributes to latitude-specific life history traits independent of the microbiota and that variation in the microbiota can suppress or reverse the differences between locally adapted fly lines. Together, these findings establish the microbiota composition of a model animal as an essential consideration in local adaptation.

Ionome and elemental transport kinetics shaped by parallel evolution in threespine stickleback.

Evidence that organisms evolve rapidly enough to alter ecological dynamics necessitates investigation of the reciprocal links between ecology and evolution. Data that link genotype to phenotype to ecology are needed to understand both the process and ecological consequences of rapid evolution. Here, we quantified the suite of elements in individuals (i.e., ionome) and differences in the fluxes of key nutrients across populations of threespine stickleback. We find that allelic variation associated with freshwater adaptation that controls bony plating is associated with changes in the ionome and nutrient recycling. More broadly, we find that adaptation of marine stickleback to freshwater conditions shifts the ionomes of natural populations and populations raised in common gardens. In both cases ionomic divergence between populations was primarily driven by differences in trace elements rather than elements typically associated with bone. These findings demonstrate the utility of ecological stoichiometry and the importance of ionome-wide data in understanding eco-evolutionary dynamics.

Parallel changes in gut microbiome composition and function during colonization, local adaptation and ecological speciation.

The processes of local adaptation and ecological speciation are often strongly shaped by biotic interactions such as competition and predation. One of the strongest lines of evidence that biotic interactions drive evolution comes from the repeated divergence of lineages in association with repeated changes in the community of interacting species. Yet relatively little is known about the repeatability of changes in gut microbial communities and their role in adaptation and divergence of host populations in nature. Here we use three cases of rapid, parallel adaptation and speciation in freshwater threespine stickleback to test for parallel changes in associated gut microbiomes. We find that features of the gut microbial communities have shifted repeatedly in the same direction in association with parallel divergence and speciation of stickleback hosts. These results suggest that changes to gut microbiomes can occur rapidly and predictably in conjunction with host evolution, and that host-microbe interactions might play an important role in host adaptation and diversification.

Detecting the macroevolutionary signal of species interactions.

Species interactions lie at the heart of many theories of macroevolution, from adaptive radiation to the Red Queen. Although some theories describe the imprint that interactions will have over long timescales, we are still missing a comprehensive understanding of the effects of interactions on macroevolution. Current research shows strong evidence for the impact of interactions on macroevolutionary patterns of trait evolution and diversification, yet many macroevolutionary studies have only a tenuous relationship to ecological studies of interactions over shorter timescales. We review current research in this area, highlighting approaches that explicitly model species interactions and connect them to broad-scale macroevolutionary patterns. We also suggest that progress has been made by taking an integrative interdisciplinary look at individual clades. We focus on African cichlids as a case study of how this approach can be fruitful. Overall, although the evidence for species interactions shaping macroevolution is strong, further work using integrative and model-based approaches is needed to spur progress towards understanding the complex dynamics that structure communities over time and space.

Differential predation alters pigmentation in threespine stickleback (Gasterosteus aculeatus).

Animal pigmentation plays a key role in many biological interactions, including courtship and predator avoidance. Sympatric benthic and limnetic ecotypes of threespine stickleback (Gasterosteus aculeatus) exhibit divergent pigment patterns. To test whether differential predation by cutthroat trout contributes to the differences in pigmentation seen between the ecotypes, we used a within-generation selection experiment on F2 benthic-limnetic hybrids. After 10 months of differential selection, we compared the pigmentation of fish under trout predation to control fish not exposed to trout predation. We found that stickleback exhibited more lateral barring in ponds with trout predation. Ponds with trout were also less turbid, and a greater degree of barring was negatively correlated with the magnitude of turbidity across pond replicates. A more benthic diet, a proxy for habitat use, was also correlated with greater lateral barring and green dorsal pigmentation. These patterns suggest that differential exposure to cutthroat trout predation may explain the divergence in body pigmentation between benthic and limnetic ecotypes.

Responses to simulated winter conditions differ between threespine stickleback ecotypes.

Abiotic factors can act as barriers to colonization and drive local adaptation. During colonization, organisms may cope with changes in abiotic factors using existing phenotypic plasticity, but the role of phenotypic plasticity in assisting or hindering the process of local adaptation remains unclear. To address these questions, we explore the role of winter conditions in driving divergence during freshwater colonization and the effects of plasticity on local adaptation in ancestral marine and derived freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We found that freshwater-resident stickleback had greater tolerance of acute exposure to low temperatures than marine stickleback, but these differences were abolished after acclimation to simulated winter conditions (9L:15D photoperiod at 4 °C). Plasma chloride levels differed between the ecotypes, but showed a similar degree of plasticity between ecotypes. Gene expression of the epithelial calcium channel (ECaC) differed between ecotypes, with the freshwater ecotype demonstrating substantially greater expression than the marine ecotype, but there was no plasticity in this trait under these conditions in either ecotype. In contrast, growth (assessed as final mass) and the expression of an isoform of the electroneutral Na(+)/H(+) exchanger (NHE3) exhibited substantial change with temperature in the marine ecotype that was not observed in the freshwater ecotype under the conditions tested here, which is consistent with evolution of these traits by a process such as genetic assimilation. These data demonstrate substantial divergence in many of these traits between freshwater and marine stickleback, but also illustrate the complexity of possible relationships between plasticity and local adaptation.

Investment in boney defensive traits alters organismal stoichiometry and excretion in fish.

Understanding how trait diversification alters ecosystem processes is an important goal for ecological and evolutionary studies. Ecological stoichiometry provides a framework for predicting how traits affect ecosystem function. The growth rate hypothesis of ecological stoichiometry links growth and phosphorus (P) body composition in taxa where nucleic acids are a significant pool of body P. In vertebrates, however, most of the P is bound within bone, and organisms with boney structures can vary in terms of the relative contributions of bones to body composition. Threespine stickleback populations have substantial variation in boney armour plating. Shaped by natural selection, this variation provides a model system to study the links between evolution of bone content, elemental body composition, and P excretion. We measure carbon:nitrogen:P body composition from stickleback populations that vary in armour phenotype. We develop a mechanistic mass-balance model to explore factors affecting P excretion, and measure P excretion from two populations with contrasting armour phenotypes. Completely armoured morphs have higher body %P but excrete more P per unit body mass than other morphs. The model suggests that such differences are driven by phenotypic differences in P intake as well as body %P composition. Our results show that while investment in boney traits alters the elemental composition of vertebrate bodies, excretion rates depend on how acquisition and assimilation traits covary with boney trait investment. These results also provide a stoichiometric hypothesis to explain the repeated loss of boney armour in threespine sticklebacks upon colonizing freshwater ecosystems.

Adaptive genetic variation mediates bottom-up and top-down control in an aquatic ecosystem.

Research in eco-evolutionary dynamics and community genetics has demonstrated that variation within a species can have strong impacts on associated communities and ecosystem processes. Yet, these studies have centred around individual focal species and at single trophic levels, ignoring the role of phenotypic variation in multiple taxa within an ecosystem. Given the ubiquitous nature of local adaptation, and thus intraspecific variation, we sought to understand how combinations of intraspecific variation in multiple species within an ecosystem impacts its ecology. Using two species that co-occur and demonstrate adaptation to their natal environments, black cottonwood (Populus trichocarpa) and three-spined stickleback (Gasterosteus aculeatus), we investigated the effects of intraspecific phenotypic variation on both top-down and bottom-up forces using a large-scale aquatic mesocosm experiment. Black cottonwood genotypes exhibit genetic variation in their productivity and consequently their leaf litter subsidies to the aquatic system, which mediates the strength of top-down effects from stickleback on prey abundances. Abundances of four common invertebrate prey species and available phosphorous, the most critically limiting nutrient in freshwater systems, are dictated by the interaction between genetic variation in cottonwood productivity and stickleback morphology. These interactive effects fit with ecological theory on the relationship between productivity and top-down control and are comparable in strength to the effects of predator addition. Our results illustrate that intraspecific variation, which can evolve rapidly, is an under-appreciated driver of community structure and ecosystem function, demonstrating that a multi-trophic perspective is essential to understanding the role of evolution in structuring ecological patterns.

Testing a 'genes-to-ecosystems' approach to understanding aquatic-terrestrial linkages.

A 'genes-to-ecosystems' approach has been proposed as a novel avenue for integrating the consequences of intraspecific genetic variation with the underlying genetic architecture of a species to shed light on the relationships among hierarchies of ecological organization (genes → individuals → communities → ecosystems). However, attempts to identify genes with major effect on the structure of communities and/or ecosystem processes have been limited and a comprehensive test of this approach has yet to emerge. Here, we present an interdisciplinary field study that integrated a common garden containing different genotypes of a dominant, riparian tree, Populus trichocarpa, and aquatic mesocosms to determine how intraspecific variation in leaf litter alters both terrestrial and aquatic communities and ecosystem functioning. Moreover, we incorporate data from extensive trait screening and genome-wide association studies estimating the heritability and genes associated with litter characteristics. We found that tree genotypes varied considerably in the quality and production of leaf litter, which contributed to variation in phytoplankton abundances, as well as nutrient dynamics and light availability in aquatic mesocosms. These 'after-life' effects of litter from different genotypes were comparable to the responses of terrestrial communities associated with the living foliage. We found that multiple litter traits corresponding with aquatic community and ecosystem responses differed in their heritability. Moreover, the underlying genetic architecture of these traits was complex, and many genes contributed only a small proportion to phenotypic variation. Our results provide further evidence that genetic variation is a key component of aquatic-terrestrial linkages, but challenge the ability to predict community or ecosystem responses based on the actions of one or a few genes.

Transgenerational epigenetic inheritance increases trait variation but is not adaptive.

Understanding processes that can produce adaptive phenotypic shifts in response to rapid environmental change is critical to reducing biodiversity loss. The ubiquity of environmentally induced epigenetic marks has led to speculation that epigenetic inheritance could potentially enhance population persistence in response to environmental change. Yet, the magnitude and fitness consequences of epigenetic marks carried beyond maternal inheritance are largely unknown. Here, we tested how transgenerational epigenetic inheritance (TEI) shapes the phenotypic response of Daphnia clones to the environmental stressor Microcystis. We split individuals from each of eight genotypes into exposure and control treatments (F0 generation) and tracked the fitness of their descendants to the F3 generation. We found transgenerational epigenetic exposure to Microcystis led to reduced rates of survival and individual growth and no consistent effect on offspring production. Increase in trait variance in the F3 relative to F0 generations suggests potential for heritable bet hedging driven by TEI, which could impact population dynamics. Our findings are counter to the working hypothesis that TEI is a generally adaptive mechanism likely to prevent extinction for populations inhabiting rapidly changing environments.

Cyanotoxins accumulate in Lake St. Clair fish yet their fillets are safe to eat.

Consuming fish exposed to cyanobacterial harmful algal blooms (HABs) may be a major route of microcystin toxin exposure to humans. However, it remains unknown whether fish can accumulate and retain microcystins temporally in waterbodies with recurring seasonal HABs, particularly before and after a HAB event when fishing is active. We conducted a field study on Largemouth Bass, Northern Pike, Smallmouth Bass, Rock Bass, Walleye, White Bass, and Yellow Perch to assess the human health risks to microcystin toxicity via fish consumption. We collected 124 fish in 2016 and 2018 from Lake St. Clair, a large freshwater ecosystem in the North American Great Lakes that is actively fished pre- and post-HAB periods. Muscles were analyzed using the 2-methyl-3-methoxy-4-phenylbutyric acid (MMPB) Lemieux Oxidation method for total microcystins, which was used to perform a human health risk assessment for comparison against fish consumption advisory benchmarks available for Lake St. Clair. From this collection 35 fish livers were additionally extracted to confirm the presence of microcystins. Microcystins were detected in all livers at widely varying concentrations (1-1500 ng g-1 ww), suggesting HABs are an underappreciated and pervasive stressor to fish populations. Conversely, microcystin levels were consistently low in muscles (0-15 ng g-1 ww) and presented negligible risk, empirically supporting that fillets may be safely consumed before and after HAB events following fish consumption advisories.

Lake Erie fish safe to eat yet afflicted by algal hepatotoxins.

Toxic harmful algal blooms (HABs) pose serious threats to human health and instances of wildlife death have been documented across taxa. However, the extent of toxicological impacts on wildlife species is largely unresolved, raising uncertainty about the repercussions of increasingly severe HABs on the biodiversity and functioning of aquatic ecosystems. Here, we conducted a field study to assess human health risks from consuming fish caught across all stages of a HAB and to determine the pervasiveness of potentially harmful levels of the cosmopolitan toxin microcystin on fish populations. We collected 190 fish in 2015 and 2017 from Lake Erie, a large freshwater ecosystem that is highly productive for fisheries and is an epicenter of HABs and microcystin toxicity events. Fish muscles and livers were analyzed for total microcystins, which was used to conduct a human health risk assessment for comparison against fish consumption advisory benchmarks available for Lake Erie. We found microcystins pose low risks to human health from fillet consumption (mean 1.80 ng g-1 ww) but substantial risks to fish health and recruitment from liver concentrations measured well before and after seasonal bloom events (mean 460.13 ng g-1 ww). Our findings indicate HABs are a previously underappreciated but pervasive threat to fish populations.

Intraspecific genetic variation is critical to robust toxicological predictions of aquatic contaminants.

Environmental risk assessment is a critical tool for protecting aquatic life and its effectiveness is predicated on predicting how natural populations respond to contaminants. Yet, routine toxicity testing typically examines only one genotype, which may render risk assessments inaccurate as populations are most often composed of genetically distinct individuals. To determine the importance of intraspecific variation in the translation of toxicity testing to populations, we quantified the magnitude of genetic variation within 20 Daphnia magna clones derived from one lake using whole genome sequencing and phenotypic assays. We repeated these assays across two exposure levels of microcystins, a cosmopolitan and lethal aquatic contaminant produced by harmful algal blooms. We found considerable intraspecific genetic variation in survival, growth, and reproduction, which was amplified by microcystins exposure. Finally, using simulations we demonstrate that the common practice of employing a single genotype to calculate toxicity tolerance failed to produce an estimate within the 95% confidence interval over half of the time. These results illuminate the importance of incorporating intraspecific genetic variation into toxicity testing to reliably predict how natural populations will respond to aquatic contaminants.

Direct observation of adaptive tracking on ecological time scales in Drosophila.

Direct observation of evolution in response to natural environmental change can resolve fundamental questions about adaptation, including its pace, temporal dynamics, and underlying phenotypic and genomic architecture. We tracked the evolution of fitness-associated phenotypes and allele frequencies genome-wide in 10 replicate field populations of Drosophila melanogaster over 10 generations from summer to late fall. Adaptation was evident over each sampling interval (one to four generations), with exceptionally rapid phenotypic adaptation and large allele frequency shifts at many independent loci. The direction and basis of the adaptive response shifted repeatedly over time, consistent with the action of strong and rapidly fluctuating selection. Overall, we found clear phenotypic and genomic evidence of adaptive tracking occurring contemporaneously with environmental change, thus demonstrating the temporally dynamic nature of adaptation.

Genetics of adaptation: Experimental test of a biotic mechanism driving divergence in traits and genes.

The genes underlying adaptations are becoming known, yet the causes of selection on genes-a key step in the study of the genetics of adaptation-remains uncertain. We address this issue experimentally in a threespine stickleback species pair showing exaggerated divergence in bony defensive armor in association with competition-driven character displacement. We used semi-natural ponds to test the role of a native predator in causing divergent evolution of armor and two known underlying genes. Predator presence/absence altered selection on dorsal spines and allele frequencies at the Msx2a gene across a generation. Evolutionary trajectories of alleles at a second gene, Pitx1, and the pelvic spine trait it controls, were more variable. Our experiment demonstrates how manipulation of putative selective agents helps to identify causes of evolutionary divergence at key genes, rule out phenotypic plasticity as a sole determinant of phenotypic differences, and eliminate reliance on fitness surrogates. Divergence of predation regimes in sympatric stickleback is associated with coevolution in response to resource competition, implying a cascade of biotic interactions driving species divergence. We suggest that as divergence proceeds, an increasing number of biotic interactions generate divergent selection, causing more evolution in turn. In this way, biotic adaptation perpetuates species divergence through time during adaptive radiation in an expanding number of traits and genes.

What genomic data can reveal about eco-evolutionary dynamics.

Recognition that evolution operates on the same timescale as ecological processes has motivated growing interest in eco-evolutionary dynamics. Nonetheless, generating sufficient data to test predictions about eco-evolutionary dynamics has proved challenging, particularly in natural contexts. Here we argue that genomic data can be integrated into the study of eco-evolutionary dynamics in ways that deepen our understanding of the interplay between ecology and evolution. Specifically, we outline five major questions in the study of eco-evolutionary dynamics for which genomic data may provide answers. Although genomic data alone will not be sufficient to resolve these challenges, integrating genomic data can provide a more mechanistic understanding of the causes of phenotypic change, help elucidate the mechanisms driving eco-evolutionary dynamics, and lead to more accurate evolutionary predictions of eco-evolutionary dynamics in nature.

Low temperature and low salinity drive putatively adaptive growth differences in populations of threespine stickleback.

Colonisation can expose organisms to novel combinations of abiotic and biotic factors and drive adaptive divergence. Yet, studies investigating the interactive effects of multiple abiotic factors on the evolution of physiological traits remain rare. Here we examine the effects of low salinity, low temperature, and their interaction on the growth of three North American populations of threespine stickleback (Gasterosteus aculeatus). In north-temperate freshwater habitats, stickleback populations experience a combination of low salinity and low winter temperatures that are not experienced by the ancestral marine and anadromous populations. Here we show that both salinity and temperature, and their interaction, have stronger negative effects on marine and anadromous populations than a freshwater population. Freshwater stickleback showed only a ~20% reduction in specific growth rate when exposed to 4 °C, while marine and anadromous stickleback showed sharp declines (82% and 74% respectively) under these conditions. The modest decreases in growth in freshwater stickleback in fresh water in the cold strongly suggest that this population has the capacity for physiological compensation to offset the negative thermodynamic effects of low temperature on growth. These results are suggestive of adaptive evolution in response to the interactive effects of low salinity and low temperature during freshwater colonisation.

Evosystem Services: Rapid Evolution and the Provision of Ecosystem Services.

Evolution is recognized as the source of all organisms, and hence many ecosystem services. However, the role that contemporary evolution might play in maintaining and enhancing specific ecosystem services has largely been overlooked. Recent advances at the interface of ecology and evolution have demonstrated how contemporary evolution can shape ecological communities and ecosystem functions. We propose a definition and quantitative criteria to study how rapid evolution affects ecosystem services (here termed contemporary evosystem services) and present plausible scenarios where such services might exist. We advocate for the direct measurement of contemporary evosystem services to improve understanding of how changing environments will alter resource availability and human well-being, and highlight the potential utility of managing rapid evolution for future ecosystem services.