Experimental evolution reveals the synergistic genomic mechanisms of adaptation to ocean warming and acidification in a marine copepod

Highly structured populations of copepods at risk to deep-sea mining: Integration of genomic data with demogenetic and biophysical modelling

Copepoda is the most abundant taxon in deep-sea hydrothermal vents, where hard substrate is available. Despite the increasing interest in seafloor massive sulphides exploitation, there have been no population genomic studies conducted on vent meiofauna, which are known to contribute over 50% to metazoan biodiversity at vents. To bridge this knowledge gap, restriction-site-associated DNA sequencing, specifically 2b-RADseq, was used to retrieve thousands of genome-wide single-nucleotide polymorphisms (SNPs) from abundant populations of the vent-obligate copepod Stygiopontius lauensis from the Lau Basin. SNPs were used to investigate population structure, demographic histories and genotype-environment associations at a basin scale. Genetic analyses also helped to evaluate the suitability of tailored larval dispersal models and the parameterization of life-history traits that better fit the population patterns observed in the genomic dataset for the target organism. Highly structured populations were observed on both spatial and temporal scales, with divergence of populations between the north, mid, and south of the basin estimated to have occurred after the creation of the major transform fault dividing the Australian and the Niuafo'ou tectonic plate (350 kya), with relatively recent secondary contact events (<20 kya). Larval dispersal models were able to predict the high levels of structure and the highly asymmetric northward low-level gene flow observed in the genomic data. These results differ from most studies conducted on megafauna in the region, elucidating the need to incorporate smaller size when considering site prospecting for deep-sea exploitation of seafloor massive sulphides, and the creation of area-based management tools to protect areas at risk of local extinction, should mining occur.

Multistress Interplay: Time and Duration of Ocean Acidification Modulate the Toxicity of Mercury and Other Metals

The current understanding of multistress interplay assumes stresses occur in perfect synchrony, but this assumption is rarely met in the natural marine ecosystem. To understand the interplay between nonperfectly overlapped stresses in the ocean, we manipulated a multigenerational experiment (F0-F3) to explore how different temporal scenarios of ocean acidification will affect mercury toxicity in a marine copepod Pseudodiaptomus annandalei. We found that the scenario of past acidification aggravated mercury toxicity but current and persistent acidification mitigated its toxicity. We specifically performed a proteomics analysis for the copepods of F3. The results indicated that current and persistent acidification initiated the energy compensation for development and mercury efflux, whereas past acidification lacked the barrier of H+ and had dysfunction in the detoxification and efflux system, providing a mechanistic understanding of mercury toxicity under different acidification scenarios. Furthermore, we conducted a meta-analysis on marine animals, demonstrating that different acidification scenarios could alter the toxicity of several other metals, despite evidence from nonsynchronous scenarios remaining limited. Our study thus demonstrates that time and duration of ocean acidification modulate mercury toxicity in marine copepods and suggests that future studies should move beyond the oversimplified scenario of perfect synchrony in understanding multistress interaction.

The contribution of gene flow, selection, and genetic drift to five thousand years of human allele frequency change

Genomic time series from experimental evolution studies and ancient DNA datasets offer us a chance to directly observe the interplay of various evolutionary forces. We show how the genome-wide variance in allele frequency change between two time points can be decomposed into the contributions of gene flow, genetic drift, and linked selection. In closed populations, the contribution of linked selection is identifiable because it creates covariances between time intervals, and genetic drift does not. However, repeated gene flow between populations can also produce directionality in allele frequency change, creating covariances. We show how to accurately separate the fraction of variance in allele frequency change due to admixture and linked selection in a population receiving gene flow. We use two human ancient DNA datasets, spanning around 5,000 y, as time transects to quantify the contributions to the genome-wide variance in allele frequency change. We find that a large fraction of genome-wide change is due to gene flow. In both cases, after correcting for known major gene flow events, we do not observe a signal of genome-wide linked selection. Thus despite the known role of selection in shaping long-term polymorphism levels, and an increasing number of examples of strong selection on single loci and polygenic scores from ancient DNA, it appears to be gene flow and drift, and not selection, that are the main determinants of recent genome-wide allele frequency change. Our approach should be applicable to the growing number of contemporary and ancient temporal population genomics datasets.

Developmental temperature, more than long-term evolution, defines thermal tolerance in an estuarine copepod

Climate change is resulting in increasing ocean temperatures and salinity variability, particularly in estuarine environments. Tolerance of temperature and salinity change interact and thus may impact organismal resilience. Populations can respond to multiple stressors in the short-term (i.e., plasticity) or over longer timescales (i.e., adaptation). However, little is known about the short- or long-term effects of elevated temperature on the tolerance of acute temperature and salinity changes. Here, we characterized the response of the near-shore and estuarine copepod, Acartia tonsa, to temperature and salinity stress. Copepods originated from one of two sets of replicated >40 generation-old temperature-adapted lines: ambient (AM, 18°C) and ocean warming (OW, 22°C). Copepods from these lines were subjected to one and three generations at the reciprocal temperature. Copepods from all treatments were then assessed for differences in acute temperature and salinity tolerance. Development (one generation), three generations, and >40 generations of warming increased thermal tolerance compared to Ambient conditions, with development in OW resulting in equal thermal tolerance to three and >40 generations of OW. Strikingly, developmental OW and >40 generations of OW had no effect on low salinity tolerance relative to ambient. By contrast, when environmental salinity was reduced first, copepods had lower thermal tolerances. These results highlight the critical role for plasticity in the copepod climate response and suggest that salinity variability may reduce copepod tolerance to subsequent warming.

The contribution of gene flow, selection, and genetic drift to five thousand years of human allele frequency change

Genomic time series from experimental evolution studies and ancient DNA datasets offer us a chance to directly observe the interplay of various evolutionary forces. We show how the genome-wide variance in allele frequency change between two time points can be decomposed into the contributions of gene flow, genetic drift, and linked selection. In closed populations, the contribution of linked selection is identifiable because it creates covariances between time intervals, and genetic drift does not. However, repeated gene flow between populations can also produce directionality in allele frequency change, creating covariances. We show how to accurately separate the fraction of variance in allele frequency change due to admixture and linked selection in a population receiving gene flow. We use two human ancient DNA datasets, spanning around 5,000 years, as time transects to quantify the contributions to the genome-wide variance in allele frequency change. We find that a large fraction of genome-wide change is due to gene flow. In both cases, after correcting for known major gene flow events, we do not observe a signal of genome-wide linked selection. Thus despite the known role of selection in shaping long-term polymorphism levels, and an increasing number of examples of strong selection on single loci and polygenic scores from ancient DNA, it appears to be gene flow and drift, and not selection, that are the main determinants of recent genome-wide allele frequency change. Our approach should be applicable to the growing number of contemporary and ancient temporal population genomics datasets.

Comparative Population Transcriptomics Provide New Insight into the Evolutionary History and Adaptive Potential of World Ocean Krill

Genetic variation is instrumental for adaptation to changing environments but it is unclear how it is structured and contributes to adaptation in pelagic species lacking clear barriers to gene flow. Here, we applied comparative genomics to extensive transcriptome datasets from 20 krill species collected across the Atlantic, Indian, Pacific, and Southern Oceans. We compared genetic variation both within and between species to elucidate their evolutionary history and genomic bases of adaptation. We resolved phylogenetic interrelationships and uncovered genomic evidence to elevate the cryptic Euphausia similis var. armata into species. Levels of genetic variation and rates of adaptive protein evolution vary widely. Species endemic to the cold Southern Ocean, such as the Antarctic krill Euphausia superba, showed less genetic variation and lower evolutionary rates than other species. This could suggest a low adaptive potential to rapid climate change. We uncovered hundreds of candidate genes with signatures of adaptive evolution among Antarctic Euphausia but did not observe strong evidence of adaptive convergence with the predominantly Arctic Thysanoessa. We instead identified candidates for cold-adaptation that have also been detected in Antarctic fish, including genes that govern thermal reception such as TrpA1. Our results suggest parallel genetic responses to similar selection pressures across Antarctic taxa and provide new insights into the adaptive potential of important zooplankton already affected by climate change.

Genome architecture underlying salinity adaptation in the invasive copepod Eurytemora affinis species complex: A review

With climate change, habitat salinity is shifting rapidly throughout the globe. In addition, many destructive freshwater invaders are recent immigrants from saline habitats. Recently, populations of the copepod Eurytemora affinis species complex have invaded freshwater habitats multiple times independently from saline estuaries on three continents. This review discusses features of this species complex that could enhance their evolutionary potential during rapid environmental change. Remarkably, across independent freshwater invasions, natural selection has repeatedly favored the same alleles far more than expected. This high degree of parallelism is surprising, given the expectation of nonparallel evolution for polygenic adaptation. Factors such as population structure and the genome architecture underlying critical traits under selection might help drive rapid adaptation and parallel evolution. Given the preponderance of saline-to-freshwater invasions and climate-induced salinity change, the principles found here could provide invaluable insights into mechanisms operating in other systems and the potential for adaptation in a changing planet.

Simultaneous warming and acidification limit population fitness and reveal phenotype costs for a marine copepod

Phenotypic plasticity and evolutionary adaptation allow populations to cope with global change, but limits and costs to adaptation under multiple stressors are insufficiently understood. We reared a foundational copepod species, Acartia hudsonica, under ambient (AM), ocean warming (OW), ocean acidification (OA), and combined ocean warming and acidification (OWA) conditions for 11 generations (approx. 1 year) and measured population fitness (net reproductive rate) derived from six life-history traits (egg production, hatching success, survival, development time, body size and sex ratio). Copepods under OW and OWA exhibited an initial approximately 40% fitness decline relative to AM, but fully recovered within four generations, consistent with an adaptive response and demonstrating synergy between stressors. At generation 11, however, fitness was approximately 24% lower for OWA compared with the AM lineage, consistent with the cost of producing OWA-adapted phenotypes. Fitness of the OWA lineage was not affected by reversal to AM or low food environments, indicating sustained phenotypic plasticity. These results mimic those of a congener, Acartia tonsa, while additionally suggesting that synergistic effects of simultaneous stressors exert costs that limit fitness recovery but can sustain plasticity. Thus, even when closely related species experience similar stressors, species-specific costs shape their unique adaptive responses.

Shifts in survival and reproduction after chronic warming enhance the potential of a marine copepod to persist under extreme heat events

The study of a species' thermal tolerance and vital rates responses provides useful metrics to characterize its vulnerability to ocean warming. Under prolonged thermal stress, plastic and adaptive processes can adjust the physiology of organisms. Yet it is uncertain whether the species can expand their upper thermal limits to cope with rapid and extreme changes in environmental temperature. In this study, we reared the marine copepod Paracartia grani at control (19°C) and warmer conditions (25°C) for >18 generations and assessed their survival and fecundity under short-term exposure to a range of temperatures (11-34°C). After multigenerational warming, the upper tolerance to acute exposure (24 h) increased by 1-1.3°C, although this enhancement decreased to 0.3-0.8°C after longer thermal stress (7 days). Warm-reared copepods were smaller and produced significantly fewer offspring at the optimum temperature. No shift in the thermal breadth of the reproductive response was observed. Yet the fecundity rates of the warm-reared copepods in the upper thermal range were up to 21-fold higher than the control. Our results show that chronic warming improved tolerance to stress temperatures and fecundity of P. grani, therefore, enhancing its chances to persist under extreme heat events.

Phenotypic Plasticity Drives Seasonal Thermal Tolerance in a Baltic Copepod

Seasonal changes in environmental conditions require substantial physiological responses for population persistence. Phenotypic plasticity is a common mechanism to tolerate these changes, but for organisms with short generation times rapid adaptation may also be a contributing factor. Here, we aimed to disentangle the impacts of adaptation from phenotypic plasticity on thermal tolerance of the calanoid copepod Acartia hudsonica collected throughout spring and summer of a single year. We used a common garden (11 °C and 18 °C) design to determine the relative impacts of plasticity versus adaptation. Acartia hudsonica were collected from five time points across the season and thermal tolerance was determined using critical thermal maximum (CTmax) followed by additional measurements after one generation of common garden. As sea surface temperature increased through the season, field collected individuals showed corresponding increases in thermal tolerance but decreases in body size. Despite different thermal tolerances of wild collections, common garden animals did not differ in CTmax within thermal treatments. Instead, there was evidence of phenotypic plasticity where higher temperatures were tolerated by the 18 °C versus the 11 °C treatment animals across all collections. Acclimation also had significant effects on body size, with higher temperatures resulting in smaller individuals, consistent with the temperature size rule. Therefore, the differences in thermal tolerance and body size observed in field collected A. hudsonica were likely driven by plasticity rather than adaptation. However, the observed decrease in body size suggests that nutrient availability and ecosystem functioning could be impacted if temperatures consistently increase with no change in copepod abundance. This is the first record of A. hudsonica in the Baltic Sea known to the authors.

Copepod life history evolution under high- and low-food regimes

Copepods play a critical role in the carbon cycle of the planet - they mediate the sequestration of carbon into the deep ocean and are the trophic link between phytoplankton and marine food webs. Global change stressors that decrease copepod productivity create the potential for catastrophic positive feedback loops. Accordingly, a growing list of studies examine the evolutionary capacity of copepods to adapt to the two primary stressors associated with global change: warmer temperatures and lower pH. But the evolutionary capacity of copepods to adapt to changing food regimes, the third major stressor associated with global change, remains unknown. We used experimental evolution to explore how a 10-fold difference in food availability affects life history evolution in the copepod, Tisbe sp. over 2 years, and spanning 30+ generations. Different food regimes evoked evolutionary responses across the entire copepod life history: we observed evolution in body size, size-fecundity relationships and offspring investment strategies. Our results suggest that changes to food regimes reshape life histories and that cryptic evolution in traits such as body size is likely. We demonstrate that evolution in response to changes in ocean productivity will alter consumer life histories and may distort trophic links in marine foodchains. Evolution in response to changing phytoplankton productivity may alter the efficacy of the global carbon pump in ways that have not been anticipated until now.

Population-specific vulnerability to ocean change in a multistressor environment

Variation in environmental conditions across a species' range can alter their responses to environmental change through local adaptation and acclimation. Evolutionary responses, however, may be challenged in ecosystems with tightly coupled environmental conditions, where changes in the covariance of environmental factors may make it more difficult for species to adapt to global change. Here, we conduct a 3-month-long mesocosm experiment and find evidence for local adaptation/acclimation in populations of red sea urchins, Mesocentrotus franciscanus, to multiple environmental drivers. Moreover, populations differ in their response to projected concurrent changes in pH, temperature, and dissolved oxygen. Our results highlight the potential for local adaptation/acclimation to multivariate environmental regimes but suggest that thresholds in responses to a single environmental variable, such as temperature, may be more important than changes to environmental covariance. Therefore, identifying physiological thresholds in key environmental drivers may be particularly useful for preserving biodiversity and ecosystem functioning.