Designing eco-evolutionary experiments for restoration projects: Opportunities and constraints revealed during stickleback introductions.

Eco-evolutionary experiments are typically conducted in semi-unnatural controlled settings, such as mesocosms; yet inferences about how evolution and ecology interact in the real world would surely benefit from experiments in natural uncontrolled settings. Opportunities for such experiments are rare but do arise in the context of restoration ecology-where different "types" of a given species can be introduced into different "replicate" locations. Designing such experiments requires wrestling with consequential questions. (Q1) Which specific "types" of a focal species should be introduced to the restoration location? (Q2) How many sources of each type should be used-and should they be mixed together? (Q3) Which specific source populations should be used? (Q4) Which type(s) or population(s) should be introduced into which restoration sites? We recently grappled with these questions when designing an eco-evolutionary experiment with threespine stickleback (Gasterosteus aculeatus) introduced into nine small lakes and ponds on the Kenai Peninsula in Alaska that required restoration. After considering the options at length, we decided to use benthic versus limnetic ecotypes (Q1) to create a mixed group of colonists from four source populations of each ecotype (Q2), where ecotypes were identified based on trophic morphology (Q3), and were then introduced into nine restoration lakes scaled by lake size (Q4). We hope that outlining the alternatives and resulting choices will make the rationales clear for future studies leveraging our experiment, while also proving useful for investigators considering similar experiments in the future.

Freshwater fishes maintain multi-trait phenotypic stability across an environmental gradient in aqueous calcium.

Reductions in a limiting nutrient might be expected to necessitate compromises in the functional traits that depend on that nutrient; yet populations existing in locations with low levels of such nutrients often do not show the expected degradation of functional traits. Indeed, logperch (Percina caprodes), pumpkinseed sunfish (Lepomis gibbosus) and yellow perch (Perca flavescens) residing in low-calcium water in the Upper St. Lawrence River were all previously found to maintain levels of scale calcium comparable to those of conspecific populations in high-calcium water. Yet it remains possible that the maintenance of one functional trait (i.e., scale calcium) under nutrient-limited (i.e., low calcium) conditions could come at the expense of maintaining other functional traits that depend on the same nutrient. The present study therefore examines other calcium-dependent traits, specifically skeletal element sizes and bone densities in the same fish species in the same area. Using radiographs of 101 fish from the three species across four locations (two in high-calcium water and two in low-calcium water), this new work documents multi-trait "homeostasis" along the gradient of water calcium. That is, no effect of calcium regime (low-calcium vs. high-calcium) was detected on any of the measured variables. Further, effect sizes for the skeletal traits were very low - lower even than effect sizes previously documented for scale calcium. These results thus show that native fishes maintain phenotypic stability across a suite of functional traits linked to calcium regulation, perhaps pointing to an "organismal-level homeostasis" scenario rather than a "trait-level homeostasis" scenario.

Dimensionality and Modularity of Adaptive Variation: Divergence in Threespine Stickleback from Diverse Environments.

AbstractPopulations are subjected to diverse environmental conditions that affect fitness and induce evolutionary or plastic responses, resulting in phenotypic divergence. Some authors contend that such divergence is concentrated along a single major axis of trait covariance even if that axis does not lead populations directly toward a fitness optimum. Other authors argue that divergence can occur readily along many phenotype axes at the same time. We use populations of threespine stickleback (Gasterosteus aculeatus) from 14 lakes with contrasting ecological conditions to find some resolution along the continuum between these two extremes. Unlike many previous studies, we included several functional suites of traits (defensive, swimming, trophic) potentially subject to different sources of selection. We find that populations exhibit dimensionality of divergence that is high enough to preclude a history of constraint along a single axis-both for divergence in multivariate mean trait values and for the structure of trait covariances. Dimensionality varied among trait suites and were strongly influenced by the inclusion of specific traits, and integration of trait suites varied between populations. We leverage this variation into new insights about the process of divergence and suggest that similar analyses could increase understanding of other adaptive radiations.

Current water quality guidelines across North America and Europe do not protect lakes from salinization.

Human-induced salinization caused by the use of road deicing salts, agricultural practices, mining operations, and climate change is a major threat to the biodiversity and functioning of freshwater ecosystems. Yet, it is unclear if freshwater ecosystems are protected from salinization by current water quality guidelines. Leveraging an experimental network of land-based and in-lake mesocosms across North America and Europe, we tested how salinization-indicated as elevated chloride (Cl-) concentration-will affect lake food webs and if two of the lowest Cl- thresholds found globally are sufficient to protect these food webs. Our results indicated that salinization will cause substantial zooplankton mortality at the lowest Cl- thresholds established in Canada (120 mg Cl-/L) and the United States (230 mg Cl-/L) and throughout Europe where Cl- thresholds are generally higher. For instance, at 73% of our study sites, Cl- concentrations that caused a ≥50% reduction in cladoceran abundance were at or below Cl- thresholds in Canada, in the United States, and throughout Europe. Similar trends occurred for copepod and rotifer zooplankton. The loss of zooplankton triggered a cascading effect causing an increase in phytoplankton biomass at 47% of study sites. Such changes in lake food webs could alter nutrient cycling and water clarity and trigger declines in fish production. Current Cl- thresholds across North America and Europe clearly do not adequately protect lake food webs. Water quality guidelines should be developed where they do not exist, and there is an urgent need to reassess existing guidelines to protect lake ecosystems from human-induced salinization.

Neutral and adaptive drivers of genomic change in introduced brook trout (Salvelinus fontinalis) populations revealed by pooled sequencing.

Understanding the drivers of successful species invasions is important for conserving native biodiversity and for mitigating the economic impacts of introduced species. However, whole-genome resolution investigations of the underlying contributions of neutral and adaptive genetic variation in successful introductions are rare. Increased propagule pressure should result in greater neutral genetic variation, while environmental differences should elicit selective pressures on introduced populations, leading to adaptive differentiation. We investigated neutral and adaptive variation among nine introduced brook trout (Salvelinus fontinalis) populations using whole-genome pooled sequencing. The populations inhabit isolated alpine lakes in western Canada and descend from a common source, with an average of ~19 (range of 7-41) generations since introduction. We found some evidence of bottlenecks without recovery, no strong evidence of purifying selection, and little support that varying propagule pressure or differences in local environments shaped observed neutral genetic variation differences. Putative adaptive loci analysis revealed nonconvergent patterns of adaptive differentiation among lakes with minimal putatively adaptive loci (0.001%-0.15%) that did not correspond with tested environmental variables. Our results suggest that (i) introduction success is not always strongly influenced by genetic load; (ii) observed differentiation among introduced populations can be idiosyncratic, population-specific, or stochastic; and (iii) conservatively, in some introduced species, colonization barriers may be overcome by support through one aspect of propagule pressure or benign environmental conditions.

Integrating physiology and environmental dynamics to operationalize environmental DNA (eDNA) as a means to monitor freshwater macro-organism abundance.

Research has demonstrated consistent positive correlations between organism abundance and absolute environmental DNA (eDNA) concentrations. Robust correlations in laboratory experiments indicate strong functional links, suggesting the potential for eDNA to monitor organism abundance in nature. However, correlations between absolute eDNA concentrations and organism abundance in nature tend to be weaker because myriad biotic and abiotic factors influence steady-state eDNA concentrations, decoupling its direct functional link with abundance. Additional technical challenges can also weaken correlations between relative organism abundance and relative eDNA data derived from metabarcoding. Future research must account for these factors to improve the inference of organism abundance from eDNA, including integrating the effects of organism physiology on eDNA production, eDNA dynamics in lentic/lotic systems, and key environmental parameters that impact estimated steady-state concentrations. Additionally, it is critical to manage expectations surrounding the accuracy and precision that eDNA can provide - eDNA, for example, cannot provide abundance estimates comparable to intensively managed freshwater fisheries that enumerate every individual fish. Recent developments, however, are encouraging. Current methods could provide meaningful information regarding qualitative conservation thresholds and emergent research has demonstrated that eDNA concentrations in natural ecosystems can provide rough quantitative estimates of abundance, particularly when models integrate physiology and/or eDNA dynamics. Operationalizing eDNA to infer abundance will probably require more than simple correlations with organism biomass/density. Nevertheless, the future is promising - models that integrate eDNA dynamics in nature could represent an effective means to infer abundance, particularly when traditional methods are considered too "costly" or difficult to obtain.

The evolutionary ecology of fatty-acid variation: Implications for consumer adaptation and diversification.

The nutritional diversity of resources can affect the adaptive evolution of consumer metabolism and consumer diversification. The omega-3 long-chain polyunsaturated fatty acids eicosapentaenoic acid (EPA; 20:5n-3) and docosahexaenoic acid (DHA; 22:6n-3) have a high potential to affect consumer fitness, through their widespread effects on reproduction, growth and survival. However, few studies consider the evolution of fatty acid metabolism within an ecological context. In this review, we first document the extensive diversity in both primary producer and consumer fatty acid distributions amongst major ecosystems, between habitats and amongst species within habitats. We highlight some of the key nutritional contrasts that can shape behavioural and/or metabolic adaptation in consumers, discussing how consumers can evolve in response to the spatial, seasonal and community-level variation of resource quality. We propose a hierarchical trait-based approach for studying the evolution of consumers' metabolic networks and review the evolutionary genetic mechanisms underpinning consumer adaptation to EPA and DHA distributions. In doing so, we consider how the metabolic traits of consumers are hierarchically structured, from cell membrane function to maternal investment, and have strongly environment-dependent expression. Finally, we conclude with an outlook on how studying the metabolic adaptation of consumers within the context of nutritional landscapes can open up new opportunities for understanding evolutionary diversification.

Phenotypic stability in scalar calcium of freshwater fish across a wide range of aqueous calcium availability in nature.

Spatial environmental gradients can promote adaptive differences among conspecific populations as a result of local adaptation or phenotypic plasticity. Such divergence can be opposed by various constraints, including gene flow, limited genetic variation, temporal fluctuations, or developmental constraints. We focus on the constraint that can be imposed when some populations are found in locations characterized by low levels of an essential nutrient. We use scales of wild fish to investigate phenotypic effects of spatial variation in a potentially limiting nutrient-calcium. If scale calcium (we use "scalar" calcium for consistency with the physiology literature) simply reflects environmental calcium availability, we expect higher levels of scalar calcium in fish from calcium-rich water, compared to fish from calcium-poor water. To consider this "passive response" scenario, we analyzed scalar calcium concentrations from three native fish species (Lepomis gibbosus, Percina caprodes, and Perca flavescens) collected at multiple sites across a dissolved calcium gradient in the Upper St. Lawrence River. Contradicting the "passive response" scenario, we did not detect strong or consistent relationships between scalar calcium and water calcium. Instead, for a given proportional increase in water calcium across the wide environmental gradient, the corresponding proportional change in scalar calcium was much smaller. We thus favor the alternative "active homeostasis" scenario, wherein fish from calcium-poor water are better able to uptake, mobilize, and deposit calcium than are fish from calcium-rich water. We further highlight the importance of studying functional traits, such as scales, in their natural setting as opposed to only laboratory studies.

A continuum of genetic mixing for conservation management along the (mal)adaptation spectrum: A comment on Hoffmann et al.

When restoring gene flow for conservation management, genetic variation should be viewed along a continuum of genetic divergence between donor and recipient populations. On the one hand, maintaining local adaptation (low divergence between donors and recipients) can enhance conservation success in the short term. On the other hand, reducing local adaptation in the short term by increasing genetic diversity (high divergence between some donors and recipients) might have better long-term success in the face of changing environmental conditions. Both Hoffman et al. (2020) and a paper we previously published in a Special Issue on Maladaptation in Applied Conservation (Derry et al., 2019) provide frameworks and syntheses for how best to apply conservation strategies in light of genetic variation and adaptation. A key difference between these two studies was that whereas Derry et al. (2019) performed a quantitative meta-analysis, Hoffman et al. (2020) relied on case studies and theoretical considerations, yielding slightly different conclusions. We here provide a summary of the two studies and contrast of the main similarities and differences between them, while highlighting terminology used to describe and explain main concepts.

Environmental RNA: A Revolution in Ecological Resolution?

Current advancements in environmental RNA (eRNA) exploit its relatively fast turnover rate relative to environmental DNA (eDNA) to assess 'metabolically active' or temporally/spatially recent community diversity. However, this focus significantly underutilizes the trove of potential ecological information encrypted in eRNA. Here, we argue for pushing beyond current species-level eDNA detection capabilities by using eRNA to detect any organisms with unique eRNA profiles, potentially including different life-history stages, sexes, or even specific phenotypes within a species. We also discuss the future of eRNA as a means of assessing the physiological status of organisms and the ecological health of populations and communities, reflecting ecosystem-level conditions. We posit that eRNA has the potential to significantly improve the resolution of organism detection, biological monitoring, and biomonitoring applications in ecology.

Freshwater zooplankton metapopulations and metacommunities respond differently to environmental and spatial variation.

Theory predicts that population genetic structure and metacommunity structure are linked by the common processes of drift and migration, but how population genetic structure and metacommunity structure are related in nature is still unknown. Deeper understanding of the processes influencing both genetic and community diversity is vital for better predicting how environmental change will impact biodiversity patterns. We examined how crustacean zooplankton and rotifer species' metapopulation genetic structure and metacommunities respond to environmental and spatial variation both within and across four regions of boreal Canada. Metapopulation and metacommunity variation partitioning results were compared within and across the four regions. Metapopulations and metacommunities responded differently to environmental variation and spatial structure both within and across regions, as metapopulations were influenced by different environmental variables compared to metacommunities. At larger spatial scales both metapopulations and metacommunities exhibited greater spatial and environmental structuring, again responding to a different subset of environmental variables. Our findings suggest that even though both genetic and species diversity are linked by the same processes, regional variation in environmental characteristics and spatial structure influence resulting biodiversity patterns differently. To date, no other empirical research has explored relationships between entire metapopulation and metacommunity assemblages at large regional spatial scales.

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.

Causes of maladaptation.

Evolutionary biologists tend to approach the study of the natural world within a framework of adaptation, inspired perhaps by the power of natural selection to produce fitness advantages that drive population persistence and biological diversity. In contrast, evolution has rarely been studied through the lens of adaptation's complement, maladaptation. This contrast is surprising because maladaptation is a prevalent feature of evolution: population trait values are rarely distributed optimally; local populations often have lower fitness than imported ones; populations decline; and local and global extinctions are common. Yet we lack a general framework for understanding maladaptation; for instance in terms of distribution, severity, and dynamics. Similar uncertainties apply to the causes of maladaptation. We suggest that incorporating maladaptation-based perspectives into evolutionary biology would facilitate better understanding of the natural world. Approaches within a maladaptation framework might be especially profitable in applied evolution contexts - where reductions in fitness are common. Toward advancing a more balanced study of evolution, here we present a conceptual framework describing causes of maladaptation. As the introductory article for a Special Feature on maladaptation, we also summarize the studies in this Issue, highlighting the causes of maladaptation in each study. We hope that our framework and the papers in this Special Issue will help catalyze the study of maladaptation in applied evolution, supporting greater understanding of evolutionary dynamics in our rapidly changing world.

Global patterns and drivers of ecosystem functioning in rivers and riparian zones.

River ecosystems receive and process vast quantities of terrestrial organic carbon, the fate of which depends strongly on microbial activity. Variation in and controls of processing rates, however, are poorly characterized at the global scale. In response, we used a peer-sourced research network and a highly standardized carbon processing assay to conduct a global-scale field experiment in greater than 1000 river and riparian sites. We found that Earth's biomes have distinct carbon processing signatures. Slow processing is evident across latitudes, whereas rapid rates are restricted to lower latitudes. Both the mean rate and variability decline with latitude, suggesting temperature constraints toward the poles and greater roles for other environmental drivers (e.g., nutrient loading) toward the equator. These results and data set the stage for unprecedented "next-generation biomonitoring" by establishing baselines to help quantify environmental impacts to the functioning of ecosystems at a global scale.

Climate alters intraspecific variation in copepod effect traits through pond food webs.

Essential fatty acids (EFAs) are primarily generated by phytoplankton in aquatic ecosystems, and can limit the growth, development, and reproduction of higher consumers. Among the most critical of the EFAs are highly unsaturated fatty acids (HUFAs), which are only produced by certain groups of phytoplankton. Changing environmental conditions can alter phytoplankton community and fatty acid composition and affect the HUFA content of higher trophic levels. Almost no research has addressed intraspecific variation in HUFAs in zooplankton, nor intraspecific relationships of HUFAs with body size and fecundity. This is despite that intraspecific variation in HUFAs can exceed interspecific variation and that intraspecific trait variation in body size and fecundity is increasingly recognized to have an important role in food web ecology (effect traits). Our study addressed the relative influences of abiotic selection and food web effects associated with climate change on intraspecific differences and interrelationships between HUFA content, body size, and fecundity of freshwater copepods. We applied structural equation modeling and regression analyses to intraspecific variation in a dominant calanoid copepod, Leptodiatomus minutus, among a series of shallow north-temperate ponds. Climate-driven diurnal temperature fluctuations favored the coexistence of diversity of phytoplankton groups with different temperature optima and nutritive quality. This resulted in unexpected positive relationships between temperature, copepod DHA content and body size. Temperature correlated positively with diatom biovolume, and mediated relationships between copepod HUFA content and body size, and between copepod body size and fecundity. The presence of brook trout further accentuated these positive effects in warm ponds, likely through nutrient cycling and stimulation of phytoplankton resources. Climate change may have previously unrecognized positive effects on freshwater copepod DHA content, body size, and fecundity in the small, shallow bodies of inland waters that are commonly found in north-temperate landscapes.

Oxic water column methanogenesis as a major component of aquatic CH4 fluxes.

Methanogenesis has traditionally been assumed to occur only in anoxic environments, yet there is mounting, albeit indirect, evidence of methane (CH4) production in oxic marine and freshwaters. Here we present the first direct, ecosystem-scale demonstration of methanogenesis in oxic lake waters. This methanogenesis appears to be driven by acetoclastic production, and is closely linked to algal dynamics. We show that oxic water methanogenesis is a significant component of the overall CH4 budget in a small, shallow lake, and provide evidence that this pathway may be the main CH4 source in large, deep lakes and open oceans. Our results challenge the current global understanding of aquatic CH4 dynamics, and suggest a hitherto unestablished link between pelagic CH4 emissions and surface-water primary production. This link may be particularly sensitive to widespread and increasing human influences on aquatic ecosystem primary productivity.

Ecological linkages between community and genetic diversity in zooplankton among boreal shield lakes.

Ecological linkages between species diversity in communities and genetic diversity in populations have potential to influence the assembly of communities in habitats recovering from human disturbance, but few studies have attempted to synthesize relationships between these levels of biological organization, especially for locally adapted species. No such studies have been done in freshwater ecosystems despite the plethora of environmental stressors plaguing aquatic communities around the world. We present the first study to test (1) whether diversity and dissimilarity among communities and populations of a locally adapted species are correlated and (2) whether communities and population haplotypes respond differently to environmental selection and spatial structure of habitats. We used a fragment of mitochondrial DNA (mtDNA) belonging to the gene cytochrome oxidase subunit I (COI) as a neutral tag to discriminate among different population haplotype variants. In boreal lakes with different histories of exposure to anthropogenic acidification, diversity and dissimilarity metrics for crustacean zooplankton communities and locally adapted populations of an abundant and broadly distributed calanoid copepod species, Leptodiaptomus minutus, did not correlate. This discord was likely because zooplankton communities responded more strongly to acidity and acidity-related environmental variables than spatial structure of lakes, whereas the distribution of L. minutus haplotypes was more strongly governed by spatial structure of lakes than environmental selection. Although spatial structure was the dominant driver of haplotype structure among L. minutus lake populations, there were similarities in the types of environmental variables that influenced the distributions of species in communities and haplotypes in populations. How haplotype diversity among populations relates to community diversity depends on the relative influence of spatial structure of habitats and selection at each of these scales of biological organization.

Adaptive reversals in acid tolerance in copepods from lakes recovering from historical stress.

Anthropogenic habitat disturbance can often lead to rapid evolution of environmental tolerances in taxa that are able to withstand the stressor. What we do not understand, however, is how species respond when the stressor no longer exists, especially across landscapes and over a considerable length of time. Once anthropogenic disturbance is removed and if there is an ecological trade-off associated with local adaptation to such an historical stressor, then evolutionary theory would predict evolutionary reversals. On the Boreal Shield, tens of thousands of lakes acidified as a result of SO2 emissions, but many of these lakes are undergoing chemical recovery as a consequence of reduced emissions. We investigated the adaptive consequences of disturbance and recovery to zooplankton living in these lakes by asking (1) if contemporary evolution of acid tolerance had arisen among Leptodiaptomus minutus copepod populations in multiple circum-neutral lakes with and without historical acidification, (2) if L. minutus populations were adaptively responding to reversals in selection in historically acidified lakes that had recovered to pH 6.0 for at least 6-8 years, and (3) if there was a fitness trade-off for L. minutus individuals with high acid tolerance at circum-neutral pH. L. minutus populations had higher acid tolerances in circum-neutral lakes with a history of acidification than in local and distant lakes that were never acidified. However, copepods in circum-neutral acid-recovering lakes were less acid-tolerant than were copepods in lakes with longer recovery time. This adaptive reversal in acid tolerance of L. minutus populations following lake recovery was supported by the results of a laboratory experiment that indicated a fitness trade-off in copepods with high acid tolerances at circum-neutral pH. These responses appear to have a genetic basis and suggest that L. minutus is highly adaptive to changes in environmental conditions. Therefore, restoration managers should focus on removing environmental stressors, and adaptable species will be able to reverse evolutionary responses to environmental disturbance in the years following recovery.