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

Metazoan adaptation to global change relies on selection of standing genetic variation. Determining the extent to which this variation exists in natural populations, particularly for responses to simultaneous stressors, is essential to make accurate predictions for persistence in future conditions. Here, we identified the genetic variation enabling the copepod Acartia tonsa to adapt to experimental ocean warming, acidification, and combined ocean warming and acidification (OWA) over 25 generations of continual selection. Replicate populations showed a consistent polygenic response to each condition, targeting an array of adaptive mechanisms including cellular homeostasis, development, and stress response. We used a genome-wide covariance approach to partition the allelic changes into three categories: selection, drift and replicate-specific selection, and laboratory adaptation responses. The majority of allele frequency change in warming (57%) and OWA (63%) was driven by shared selection pressures across replicates, but this effect was weaker under acidification alone (20%). OWA and warming shared 37% of their response to selection but OWA and acidification shared just 1%, indicating that warming is the dominant driver of selection in OWA. Despite the dominance of warming, the interaction with acidification was still critical as the OWA selection response was highly synergistic with 47% of the allelic selection response unique from either individual treatment. These results disentangle how genomic targets of selection differ between single and multiple stressors and demonstrate the complexity that nonadditive multiple stressors will contribute to predictions of adaptation to complex environmental shifts caused by global change.

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.

Methylmercury accumulation in plankton on the continental margin of the northwest Atlantic Ocean.

Accumulation of monomethylmercury (MMHg) by plankton is a key process influencing concentrations of this toxic mercury species in marine food webs and seafood. We examined bioaccumulation and biomagnification of MMHg in microseston and four size fractions of zooplankton on the continental shelf, slope, and rise of the northwest Atlantic Ocean. The bioaccumulation factor (BAF, L/kg) for MMHg in microseston averaged 10(4.3±0.3) among 21 locations, and concentrations were unrelated to those in colocated, filtered surface water. Instead, concentrations and the BAF of MMHg in microseston were related inversely with total suspended solids in surface water, a proxy for planktonic biomass at these remote locations. MMHg was biomagnified by a factor of 4 from microseston to zooplankton, and both concentrations of MMHg and the fraction of total mercury as MMHg increased with larger size fractions of zooplankton. These results suggest that the initial magnitude of MMHg uptake into pelagic marine food webs is influenced by the degree of primary production in surface waters and propagated up through large zooplankton. Accordingly, biological productivity, in addition to inputs of MMHg to surface waters, must be considered when predicting how MMHg bioaccumulation will vary spatially and temporally in the ocean.

Loss of transcriptional plasticity but sustained adaptive capacity after adaptation to global change conditions in a marine copepod.

Adaptive evolution and phenotypic plasticity will fuel resilience in the geologically unprecedented warming and acidification of the earth's oceans, however, we have much to learn about the interactions and costs of these mechanisms of resilience. Here, using 20 generations of experimental evolution followed by three generations of reciprocal transplants, we investigated the relationship between adaptation and plasticity in the marine copepod, Acartia tonsa, in future global change conditions (high temperature and high CO2). We found parallel adaptation to global change conditions in genes related to stress response, gene expression regulation, actin regulation, developmental processes, and energy production. However, reciprocal transplantation showed that adaptation resulted in a loss of transcriptional plasticity, reduced fecundity, and reduced population growth when global change-adapted animals were returned to ambient conditions or reared in low food conditions. However, after three successive transplant generations, global change-adapted animals were able to match the ambient-adaptive transcriptional profile. Concurrent changes in allele frequencies and erosion of nucleotide diversity suggest that this recovery occurred via adaptation back to ancestral conditions. These results demonstrate that while plasticity facilitated initial survival in global change conditions, it eroded after 20 generations as populations adapted, limiting resilience to new stressors and previously benign environments.

Determining the Advantages, Costs, and Trade-Offs of a Novel Sodium Channel Mutation in the Copepod Acartia hudsonica to Paralytic Shellfish Toxins (PST).

The marine copepod Acartia hudsonica was shown to be adapted to dinoflagellate prey, Alexandrium fundyense, which produce paralytic shellfish toxins (PST). Adaptation to PSTs in other organisms is caused by a mutation in the sodium channel. Recently, a mutation in the sodium channel in A. hudsonica was found. In this study, we rigorously tested for advantages, costs, and trade-offs associated with the mutant isoform of A. hudsonica under toxic and non-toxic conditions. We combined fitness with wild-type: mutant isoform ratio measurements on the same individual copepod to test our hypotheses. All A. hudsonica copepods express both the wild-type and mutant sodium channel isoforms, but in different proportions; some individuals express predominantly mutant (PMI) or wild-type isoforms (PWI), while most individuals express relatively equal amounts of each (EI). There was no consistent pattern of improved performance as a function of toxin dose for egg production rate (EPR), ingestion rate (I), and gross growth efficiency (GGE) for individuals in the PMI group relative to individuals in the PWI expression group. Neither was there any evidence to indicate a fitness benefit to the mutant isoform at intermediate toxin doses. No clear advantage under toxic conditions was associated with the mutation. Using a mixed-diet approach, there was also no observed relationship between individual wild-type: mutant isoform ratios and among expression groups, on both toxic and non-toxic diets, for eggs produced over three days. Lastly, expression of the mutant isoform did not mitigate the negative effects of the toxin. That is, the reductions in EPR from a toxic to non-toxic diet for copepods were independent of expression groups. Overall, the results did not support our hypotheses; the mutant sodium channel isoform does not appear to be related to adaptation to PST in A. hudsonica. Other potential mechanisms responsible for the adaptation are discussed.

No evidence for induction or selection of mutant sodium channel expression in the copepod Acartia husdsonica challenged with the toxic dinoflagellate Alexandrium fundyense.

Some species in the dinoflagellate genus Alexandrium spp. produce a suite of neurotoxins that block sodium channels, known as paralytic shellfish toxins (PST), which have deleterious effects on grazers. Populations of the ubiquitous copepod grazer Acartia hudsonica that have co-occurred with toxic Alexandrium spp. are better adapted than naïve populations. The mechanism of adaptation is currently unknown. We hypothesized that a mutation in the sodium channel could account for the grazer adaptation. We tested two hypotheses: (1) Expression of the mutant sodium channel could be induced by exposure to toxic Alexandrium fundyense; (2) in the absence of induction, selection exerted by toxic A. fundyense would favor copepods that predominantly express the mutant isoform. In the copepod A. hudsonica, both isoforms are expressed in all individuals in varying proportions. Thus, in addition to comparing expression ratios of wild-type to mutant isoforms for individual copepods, we also partitioned copepods into three groups: those that predominantly express the mutant (PMI) isoform, the wild-type (PWI) isoform, or both isoforms approximately equally (EI). There were no differences in isoform expression between individuals that were fed toxic and nontoxic food after three and 6 days; induction of mutant isoform expression did not occur. Furthermore, the hypothesis that mutant isoform expression responds to toxic food was also rejected. That is, no consistent evidence showed that the wild-type to mutant isoform ratios decreased, or that the relative proportion of PMI individuals increased, due to the consumption of toxic food over four generations. However, in the selected line that was continuously exposed to toxic food sources, egg production rate increased, which suggested that adaptation occurred but was unrelated to sodium channel isoform expression.