The hologenome of Daphnia magna reveals possible DNA methylation and microbiome-mediated evolution of the host genome.

Properties that make organisms ideal laboratory models in developmental and medical research are often the ones that also make them less representative of wild relatives. The waterflea Daphnia magna is an exception, by both sharing many properties with established laboratory models and being a keystone species, a sentinel species for assessing water quality, an indicator of environmental change and an established ecotoxicology model. Yet, Daphnia's full potential has not been fully exploited because of the challenges associated with assembling and annotating its gene-rich genome. Here, we present the first hologenome of Daphnia magna, consisting of a chromosomal-level assembly of the D. magna genome and the draft assembly of its metagenome. By sequencing and mapping transcriptomes from exposures to environmental conditions and from developmental morphological landmarks, we expand the previously annotates gene set for this species. We also provide evidence for the potential role of gene-body DNA-methylation as a mutagen mediating genome evolution. For the first time, our study shows that the gut microbes provide resistance to commonly used antibiotics and virulence factors, potentially mediating Daphnia's environmental-driven rapid evolution. Key findings in this study improve our understanding of the contribution of DNA methylation and gut microbiota to genome evolution in response to rapidly changing environments.

Historical exposure to chemicals reduces tolerance to novel chemical stress in Daphnia (waterflea).

Until the last few decades, anthropogenic chemicals used in most production processes have not been comprehensively assessed for their risk and impact on wildlife and humans. They are transported globally and usually end up in the environment as unintentional pollutants, causing long-term adverse effects. Modern toxicology practices typically use acute toxicity tests of unrealistic concentrations of chemicals to determine their safe use, missing pathological effects arising from long-term exposures to environmentally relevant concentrations. Here, we study the transgenerational effect of environmentally relevant concentrations of five chemicals on the priority list of international regulatory frameworks on the keystone species Daphnia magna. We expose Daphnia genotypes resurrected from the sedimentary archive of a lake with a known history of chemical pollution to the five chemicals to understand how historical exposure to chemicals influences adaptive responses to novel chemical stress. We measure within- and transgenerational plasticity in fitness-linked life history traits following exposure of "experienced" and "naive" genotypes to novel chemical stress. As the revived Daphnia originate from the same genetic pool sampled at different times in the past, we are able to quantify the long-term evolutionary impact of chemical pollution by studying genome-wide diversity and identifying functional pathways affected by historical chemical stress. Our results suggest that historical exposure to chemical stress causes reduced genome-wide diversity, leading to lower cross-generational tolerance to novel chemical stress. Lower tolerance is underpinned by reduced gene diversity at detoxification, catabolism and endocrine genes in experienced genotypes. We show that these genes sit within pathways that are conserved and potential chemical targets in other species, including humans.

Daphnia as a Sentinel Species for Environmental Health Protection: A Perspective on Biomonitoring and Bioremediation of Chemical Pollution.

Despite available technology and the knowledge that chemical pollution damages human and ecosystem health, chemical pollution remains rampant, ineffectively monitored, rarely prevented, and only occasionally mitigated. We present a framework that helps address current major challenges in the monitoring and assessment of chemical pollution by broadening the use of the sentinel species Daphnia as a diagnostic agent of water pollution. And where prevention has failed, we propose the application of Daphnia as a bioremediation agent to help reduce hazards from chemical mixtures in the environment. By applying "omics" technologies to Daphnia exposed to real-world ambient chemical mixtures, we show improvements at detecting bioactive components of chemical mixtures, determining the potential effects of untested chemicals within mixtures, and identifying targets of toxicity. We also show that using Daphnia strains that naturally adapted to chemical pollution as removal agents of ambient chemical mixtures can sustainably improve environmental health protection. Expanding the use of Daphnia beyond its current applications in regulatory toxicology has the potential to improve both the assessment and the remediation of environmental pollution.

Adaptive responses along a depth and a latitudinal gradient in the endemic seagrass Posidonia oceanica.

Seagrass meadows provide important ecosystem services and are critical for the survival of the associated invertebrate community. However, they are threatened worldwide by human-driven environmental change. Understanding the seagrasses' potential for adaptation is critical to assess not only their ability to persist under future global change scenarios, but also to assess the persistence of the associated communities. Here we screened a wild population of Posidonia oceanica, an endemic long-lived seagrass in the Mediterranean Sea, for genes that may be target of environmental selection, using an outlier and a genome-wide transcriptome analysis. We identified loci where polymorphism or differential expression was associated with either a latitudinal or a bathymetric gradient, as well as with both gradients in an effort to identify loci associated with temperature and light. We found the candidate genes underlying growth and immunity to be divergent between populations adapted to different latitudes and/or depths, providing evidence for local adaptation. Furthermore, we found evidence of reduced gene flow among populations including adjacent populations. Reduced gene flow, combined with low sexual recombination, small effective population size, and long generation time of P. oceanica raises concerns for the long-term persistence of this species, especially in the face of rapid environmental change driven by human activities.

Evolution of thermal tolerance in multifarious environments.

Species extinction rates are many times greater than the direst predictions made two decades ago by environmentalists, largely because of human impact. Major concerns are associated with the predicted higher recurrence and severity of extreme events, such as heat waves. Although tolerance to these extreme events is instrumental to species survival, little is known whether and how it evolves in natural populations, and to what extent it is affected by other environmental stressors. Here, we study physiological and molecular mechanisms of thermal tolerance over evolutionary times in multifarious environments. Using the practice of "resurrection ecology" on the keystone grazer Daphnia magna, we quantified genetic and plastic differences in physiological and molecular traits linked to thermal tolerance in historical and modern genotypes of the same population. This population experienced an increase in average temperature and occurrence of heat waves, in addition to dramatic changes in water chemistry, over five decades. On genotypes resurrected across the five decades, we measured plastic and genetic differences in CTmax , body size, Hb content and differential expression of four heat shock proteins after exposure to temperature as single stress and in combination with food levels and insecticide loads. We observed evolution of the critical thermal maximum and plastic response in body size, HSP expression and Hb content over time in a warming only scenario. Molecular and physiological responses to extreme temperature in multifarious environments were not predictable from the response to warming alone. Underestimating the effect of multiple stressors on thermal tolerance can lead to wrong estimates of species evolvability and persistence.

Early transcriptional response pathways in Daphnia magna are coordinated in networks of crustacean-specific genes.

Natural habitats are exposed to an increasing number of environmental stressors that cause important ecological consequences. However, the multifarious nature of environmental change, the strength and the relative timing of each stressor largely limit our understanding of biological responses to environmental change. In particular, early response to unpredictable environmental change, critical to survival and fitness in later life stages, is largely uncharacterized. Here, we characterize the early transcriptional response of the keystone species Daphnia magna to twelve environmental perturbations, including biotic and abiotic stressors. We first perform a differential expression analysis aimed at identifying differential regulation of individual genes in response to stress. This preliminary analysis revealed that a few individual genes were responsive to environmental perturbations and they were modulated in a stressor and genotype-specific manner. Given the limited number of differentially regulated genes, we were unable to identify pathways involved in stress response. Hence, to gain a better understanding of the genetic and functional foundation of tolerance to multiple environmental stressors, we leveraged the correlative nature of networks and performed a weighted gene co-expression network analysis. We discovered that approximately one-third of the Daphnia genes, enriched for metabolism, cell signalling and general stress response, drives transcriptional early response to environmental stress and it is shared among genetic backgrounds. This initial response is followed by a genotype- and/or condition-specific transcriptional response with a strong genotype-by-environment interaction. Intriguingly, genotype- and condition-specific transcriptional response is found in genes not conserved beyond crustaceans, suggesting niche-specific adaptation.

Bisulfite Sequencing with Daphnia Highlights a Role for Epigenetics in Regulating Stress Response to Microcystis through Preferential Differential Methylation of Serine and Threonine Amino Acids.

Little is known about the influence that environmental stressors may have on genome-wide methylation patterns, and to what extent epigenetics may be involved in environmental stress response. Yet, studies of methylation patterns under stress could provide crucial insights on stress response and toxicity pathways. Here, we focus on genome-wide methylation patterns in the microcrustacean Daphnia magna, a model organism in ecotoxicology and risk assessment, exposed to the toxic cyanobacterium Microcystis aeruginosa. Bisulfite sequencing of exposed and control animals highlighted differential methylation patterns in Daphnia upon exposure to Microcystis primarily in exonic regions. These patterns are enriched for serine/threonine amino acid codons and genes related to protein synthesis, transport and degradation. Furthermore, we observed that genes with differential methylation corresponded well with genes susceptible to alternative splicing in response to Microcystis stress. Overall, our results suggest a complex mechanistic response in Daphnia characterized by interactions between DNA methylation and gene regulation mechanisms. These results underscore that DNA methylation is modulated by environmental stress and can also be an integral part of the toxicity response in our study species.

Temporal genetic stability in natural populations of the waterflea Daphnia magna in response to strong selection pressure.

Studies monitoring changes in genetic diversity and composition through time allow a unique understanding of evolutionary dynamics and persistence of natural populations. However, such studies are often limited to species with short generation times that can be propagated in the laboratory or few exceptional cases in the wild. Species that produce dormant stages provide powerful models for the reconstruction of evolutionary dynamics in the natural environment. A remaining open question is to what extent dormant egg banks are an unbiased representation of populations and hence of the species' evolutionary potential, especially in the presence of strong environmental selection. We address this key question using the water flea Daphnia magna, which produces dormant stages that accumulate in biological archives over time. We assess temporal genetic stability in three biological archives, previously used in resurrection ecology studies showing adaptive evolutionary responses to rapid environmental change. We show that neutral genetic diversity does not decline with the age of the population and it is maintained in the presence of strong selection. In addition, by comparing temporal genetic stability in hatched and unhatched populations from the same biological archive, we show that dormant egg banks can be consulted to obtain a reliable measure of genetic diversity over time, at least in the multidecadal time frame studied here. The stability of neutral genetic diversity through time is likely mediated by the buffering effect of the resting egg bank.

An integrated spatio-temporal view of riverine biodiversity using environmental DNA metabarcoding.

Anthropogenically forced changes in global freshwater biodiversity demand more efficient monitoring approaches. Consequently, environmental DNA (eDNA) analysis is enabling ecosystem-scale biodiversity assessment, yet the appropriate spatio-temporal resolution of robust biodiversity assessment remains ambiguous. Here, using intensive, spatio-temporal eDNA sampling across space (five rivers in Europe and North America, with an upper range of 20-35 km between samples), time (19 timepoints between 2017 and 2018) and environmental conditions (river flow, pH, conductivity, temperature and rainfall), we characterise the resolution at which information on diversity across the animal kingdom can be gathered from rivers using eDNA. In space, beta diversity was mainly dictated by turnover, on a scale of tens of kilometres, highlighting that diversity measures are not confounded by eDNA from upstream. Fish communities showed nested assemblages along some rivers, coinciding with habitat use. Across time, seasonal life history events, including salmon and eel migration, were detected. Finally, effects of environmental conditions were taxon-specific, reflecting habitat filtering of communities rather than effects on DNA molecules. We conclude that riverine eDNA metabarcoding can measure biodiversity at spatio-temporal scales relevant to species and community ecology, demonstrating its utility in delivering insights into river community ecology during a time of environmental change.

The role of selection in driving landscape genomic structure of the waterflea Daphnia magna.

The combined analysis of neutral and adaptive genetic variation is crucial to reconstruct the processes driving population genetic structure in the wild. However, such combined analysis is challenging because of the complex interaction among neutral and selective processes in the landscape. Overcoming this level of complexity requires an unbiased search for the evidence of selection in the genomes of populations sampled from their natural habitats and the identification of demographic processes that lead to present-day populations' genetic structure. Ecological model species with a suite of genomic tools and well-understood ecologies are best suited to resolve this complexity and elucidate the role of selective and demographic processes in the landscape genomic structure of natural populations. Here we investigate the water flea Daphnia magna, an emerging model system in genomics and a renowned ecological model system. We infer past and recent demographic processes by contrasting patterns of local and regional neutral genetic diversity at markers with different mutation rates. We assess the role of the environment in driving genetic variation in our study system by identifying correlates between biotic and abiotic variables naturally occurring in the landscape and patterns of neutral and adaptive genetic variation. Our results indicate that selection plays a major role in determining the population genomic structure of D. magna. First, environmental selection directly impacts genetic variation at loci hitchhiking with genes under selection. Second, priority effects enhanced by local genetic adaptation (cf. monopolization) affect neutral genetic variation by reducing gene flow among populations and genetic diversity within populations.

Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization.

Empirical population genetic studies have been dominated by a neutralist view, according to which gene flow and drift are the main forces driving population genetic structure in nature. The neutralist view in essence describes a process of isolation by dispersal limitation (IBDL) that generally leads to a pattern of isolation by distance (IBD). Recently, however, conceptual frameworks have been put forward that view local genetic adaptation as an important driver of population genetic structure. Isolation by adaptation (IBA) and monopolization (M) posit that gene flow among natural populations is reduced as a consequence of local genetic adaptation. IBA stresses that effective gene flow is reduced among habitats that show dissimilar ecological characteristics, leading to a pattern of isolation by environment. In monopolization, local genetic adaptation of initial colonizing genotypes results in a reduction in gene flow that fosters the persistence of founder effects. Here, we relate these different processes driving landscape genetic structure to patterns of IBD and isolation by environment (IBE). We propose a method to detect whether IBDL, IBA and M shape genetic differentiation in natural landscapes by studying patterns of variation at neutral and non-neutral markers as well as at ecologically relevant traits. Finally, we reinterpret a representative number of studies from the recent literature by associating patterns to processes and identify patterns associated with local genetic adaptation to be as common as IBDL in structuring regional genetic variation of populations in the wild. Our results point to the importance of quantifying environmental gradients and incorporating ecology in the analysis of population genetics.

100 years of anthropogenic impact causes changes in freshwater functional biodiversity.

Despite efforts from scientists and regulators, biodiversity is declining at an alarming rate. Unless we find transformative solutions to preserve biodiversity, future generations may not be able to enjoy nature's services. We have developed a conceptual framework that establishes the links between biodiversity dynamics and abiotic change through time and space using artificial intelligence. Here, we apply this framework to a freshwater ecosystem with a known history of human impact and study 100 years of community-level biodiversity, climate change and chemical pollution trends. We apply explainable network models with multimodal learning to community-level functional biodiversity measured with multilocus metabarcoding, to establish correlations with biocides and climate change records. We observed that the freshwater community assemblage and functionality changed over time without returning to its original state, even if the lake partially recovered in recent times. Insecticides and fungicides, combined with extreme temperature events and precipitation, explained up to 90% of the functional biodiversity changes. The community-level biodiversity approach used here reliably explained freshwater ecosystem shifts. These shifts were not observed when using traditional quality indices (e.g. Trophic Diatom Index). Our study advocates the use of high-throughput systemic approaches on long-term trends over species-focused ecological surveys to identify the environmental factors that cause loss of biodiversity and disrupt ecosystem functions.

Over long periods of time, environmental changes ­ such as chemical pollution and climate change ­ affect the diversity of organisms that live in an ecosystem, known as 'biodiversity'. Understanding the impact of these changes is challenging because they can happen slowly, their effect is only measurable after years, and historical records are limited. This can make it difficult to determine when specific changes happened, what might have driven them and what impact they might be having. One way to measure changes in biodiversity over time is by analysing traces of DNA shed by organisms. Plants, animals, and bacteria living in lakes leave behind genetic material that gets trapped and buried in the sediment at the bottom of lakes. Similarly, biocides ­ substances used to kill or control populations of living organisms ­ that run-off into lakes leach into the sediment and can be measured years later. Therefore, this sediment holds a record of life and environmental impacts in the lake over past centuries. Eastwood, Zhou et al. wanted to understand the relationship between environmental changes (such as the use of biocides and climate change) and shifts in lake biodiversity. To do so, the researchers studied a lake community that had experienced major environmental impacts over the last century (including nutrient pollution, chemical pollution and climate change), but which appeared to improve over the last few years of the 20^th century. Using machine learning to find connections over time between biodiversity and non-living environmental changes, Eastwood, Zhou et al. showed that, despite apparent recovery in water quality, the biodiversity of the lake was not restored to its original state. A combination of climate factors (such as rainfall levels and extreme temperatures) and biocide application (particularly insecticides and fungicides) explained up to 90% of the biodiversity changes that occurred in the lake. These changes had not been identified before using traditional techniques. The functional roles microorganisms played in the ecosystem (such as degradation and nitrogen metabolism) were also altered, suggesting that loss of biodiversity may lead to loss of ecosystem functions. The findings described by Eastwood, Zhou et al. can be used by environmental regulators to identify species or ecosystems at risk from environmental change and prioritise them for intervention. The approach can also be used to identify which chemicals pose the greatest threat to biodiversity. Additionally, the use of environmental DNA from sediment can provide rich historical biodiversity data, which can be used to train artificial intelligence-based models to improve predictions of how ecosystems will respond to complex environmental changes.

Harnessing water fleas for water reclamation: A nature-based tertiary wastewater treatment technology.

Urbanisation, population growth, and climate change have put unprecedented pressure on water resources, leading to a global water crisis and the need for water reuse. However, water reuse is unsafe unless persistent chemical pollutants are removed from reclaimed water. State-of-the-art technologies for the reduction of persistent chemical pollutants in wastewater typically impose high operational and energy costs and potentially generate toxic by-products (e.g., bromate from ozonation). Nature-base solutions are preferred to these technologies for their lower environmental impact. However, so far, bio-based tertiary wastewater treatments have been inefficient for industrial-scale applications. Moreover, they often demand significant financial investment and large infrastructure, undermining sustainability objectives. Here, we present a scalable, low-cost, low-carbon, and retrofittable nature-inspired solution to remove persistent chemical pollutants (pharmaceutical, pesticides and industrial chemicals). We showed Daphnia's removal efficiency of individual chemicals and chemicals from wastewater at laboratory scale ranging between 50 % for PFOS and 90 % for diclofenac. We validated the removal efficiency of diclofenac at prototype scale, showing sustained performance over four weeks in outdoor seminatural conditions. A techno-commercial analysis on the Daphnia-based technology suggested several technical, commercial and sustainability advantages over established and emerging treatments at comparable removal efficiency, benchmarked on available data on individual chemicals. Further testing of the technology is underway in open flow environments holding real wastewater. The technology has the potential to improve the quality of wastewater effluent, meeting requirements to produce water appropriate for reuse in irrigation, industrial application, and household use. By preventing persistent chemicals from entering waterways, this technology has the potential to maximise the shift to clean growth, enabling water reuse, reducing resource depletion and preventing environmental pollution.

Genomic signature of natural and anthropogenic stress in wild populations of the waterflea Daphnia magna: validation in space, time and experimental evolution.

Natural populations are confronted with multiple selection pressures resulting in a mosaic of environmental stressors at the landscape level. Identifying the genetic underpinning of adaptation to these complex selection environments and assigning causes of natural selection within multidimensional selection regimes in the wild is challenging. The water flea Daphnia is a renowned ecological model system with its well-documented ecology, the possibility to analyse subfossil dormant egg banks and the short generation time allowing an experimental evolution approach. Capitalizing on the strengths of this model system, we here link candidate genome regions to three selection pressures, known to induce micro-evolutionary responses in Daphnia magna: fish predation, parasitism and land use. Using a genome scan approach in space, time and experimental evolution trials, we provide solid evidence of selection at the genome level under well-characterized environmental gradients in the wild and identify candidate genes linked to the three environmental stressors. Our study reveals differential selection at the genome level in Daphnia populations and provides evidence for repeatable patterns of local adaptation in a geographic mosaic of environmental stressors fuelled by standing genetic variation. Our results imply high evolutionary potential of local populations, which is relevant to understand the dynamics of trait changes in natural populations and their impact on community and ecosystem responses through eco-evolutionary feedbacks.

The Time Machine framework: monitoring and prediction of biodiversity loss.

Transdisciplinary solutions are needed to achieve the sustainability of ecosystem services for future generations. We propose a framework to identify the causes of ecosystem function loss and to forecast the future of ecosystem services under different climate and pollution scenarios. The framework (i) applies an artificial intelligence (AI) time-series analysis to identify relationships among environmental change, biodiversity dynamics and ecosystem functions; (ii) validates relationships between loss of biodiversity and environmental change in fabricated ecosystems; and (iii) forecasts the likely future of ecosystem services and their socioeconomic impact under different pollution and climate scenarios. We illustrate the framework by applying it to watersheds, and provide system-level approaches that enable natural capital restoration by associating multidecadal biodiversity changes to chemical pollution.

Single nucleotide polymorphism discovery from expressed sequence tags in the waterflea Daphnia magna.

BACKGROUND: Daphnia (Crustacea: Cladocera) plays a central role in standing aquatic ecosystems, has a well known ecology and is widely used in population studies and environmental risk assessments. Daphnia magna is, especially in Europe, intensively used to study stress responses of natural populations to pollutants, climate change, and antagonistic interactions with predators and parasites, which have all been demonstrated to induce micro-evolutionary and adaptive responses. Although its ecology and evolutionary biology is intensively studied, little is known on the functional genomics underpinning of phenotypic responses to environmental stressors. The aim of the present study was to find genes expressed in presence of environmental stressors, and target such genes for single nucleotide polymorphic (SNP) marker development. RESULTS: We developed three expressed sequence tag (EST) libraries using clonal lineages of D. magna exposed to ecological stressors, namely fish predation, parasite infection and pesticide exposure. We used these newly developed ESTs and other Daphnia ESTs retrieved from NCBI GeneBank to mine for SNP markers targeting synonymous as well as non synonymous genetic variation. We validate the developed SNPs in six natural populations of D. magna distributed at regional scale. CONCLUSIONS: A large proportion (47%) of the produced ESTs are Daphnia lineage specific genes, which are potentially involved in responses to environmental stress rather than to general cellular functions and metabolic activities, or reflect the arthropod's aquatic lifestyle. The characterization of genes expressed under stress and the validation of their SNPs for population genetic study is important for identifying ecologically responsive genes in D. magna.

Genomics in the ecological arena.

This meeting report presents the cutting-edge research that is developing around the waterflea Daphnia, an emerging model system in environmental genomics. Daphnia has been a model species in ecology, toxicology and evolution for many years and is supported by a large community of ecologists, evolutionary biologists and ecotoxicologists. Thanks to new advances in genomics and transciptomics and to the sustained efforts of the Daphnia Genomics Consortium (DGC), Daphnia is also rapidly developing as a model system in environmental genomics. Advances in this emerging field were presented at the DGC 2010, held for the first time in a European University. During the meeting, a plethora of elegant studies were presented on the mechanisms of responses to environmental challenges using recently developed genomic tools. The DGC 2010 is a concrete example of the new trends in ecology and evolution. The times are mature for the application of innovative genomic and transcriptomic tools for studies of environmental genomics in non-model organisms.

Protocol for assay of transposase accessible chromatin sequencing in non-model species.

The assay for transposase accessible chromatin (ATAC-seq) is a method for mapping genome-wide chromatin accessibility. Coupled with high-throughput sequencing, it enables integrative epigenomics analyses. ATAC-seq requires direct access to cell nuclei, a major challenge in non-model species such as small invertebrates, whose soft tissue is surrounded by a protective exoskeleton. Here, we present modifications of the ATAC-seq protocol for applications in small crustaceans, extending applications to non-model species. For complete information on the use and execution of this protocol, please refer to Buenrostro et al. (2013).

Evolutionary mechanisms underpinning fitness response to multiple stressors in Daphnia.

Multiple stressors linked to anthropogenic activities can influence how organisms adapt and evolve. So far, a consensus on how multiple stressors drive adaptive trajectories in natural populations has not been reached. Some meta-analysis reports show predominance of additive effects of stressors on ecological endpoints (e.g., fecundity, mortality), whereas others show synergistic effects more frequently. Moreover, it is unclear what mechanisms of adaptation underpin responses to complex environments. Here, we use populations of Daphnia magna resurrected from different times in the past to investigate mechanisms of adaptation to multiple stressors and to understand how historical exposure to environmental stress shapes adaptive responses of modern populations. Using common garden experiments on resurrected modern and historical populations, we investigate (1) whether exposure to one stress results in higher tolerance to a second stressor; (2) the mechanisms of adaptation underpinning long-term evolution to multistress (genetic evolution, plasticity, evolution of plasticity); and (3) the interaction effects of multiple stressors on fitness (synergism, antagonism, additivity). We measure the combined impact of different levels of resource availability (algae) and biocides on fitness-linked life-history traits and interpret these results in light of historical environmental exposures. We show that exposure to one stressor can alter tolerance to second stressors and that the interaction effect depends on the severity of either stressor. We also show that mechanisms of adaptation underpinning phenotypic evolution significantly differ in single-stress and multistress scenarios. These adaptive responses are driven largely by synergistic effects on fecundity and size at maturity, and additive effects on age at maturity. Exposure to multiple stressors shifts the trade-offs among fitness-linked life-history traits, with a stronger effect on Daphnia populations when low-resource availability and high biocide levels are experienced. Our study indicates that mitigation interventions based on single-stress analysis may not capture realistic threats.

Extensive standing genetic variation from a small number of founders enables rapid adaptation in Daphnia.

We lack a thorough understanding of the origin and maintenance of standing genetic variation that enables rapid evolutionary responses of natural populations. Whole genome sequencing of a resurrected Daphnia population shows that standing genetic variation in over 500 genes follows an evolutionary trajectory that parallels the pronounced and rapid adaptive evolution of multiple traits in response to predator-driven natural selection and its subsequent relaxation. Genetic variation carried by only five founding individuals from the regional genotype pool is shown to suffice at enabling the observed evolution. Our results provide insight on how natural populations can acquire the genomic variation, through colonization by a few regional genotypes, that fuels rapid evolution in response to strong selection pressures. While these evolutionary responses in our study population involved hundreds of genes, we observed no evidence of genetic erosion.