High evolutionary potential of marine zooplankton

Abstract Open ocean zooplankton often have been viewed as slowly evolving species that have limited capacity to respond adaptively to changing ocean conditions. Hence, attention has focused on the ecological responses of zooplankton to current global change, including range shifts and changing phenology. Here, we argue that zooplankton also are well poised for evolutionary responses to global change. We present theoretical arguments that suggest plankton species may respond rapidly to selection on mildly beneficial mutations due to exceptionally large population size, and consider the circumstantial evidence that supports our inference that selection may be particularly important for these species. We also review all primary population genetic studies of open ocean zooplankton and show that genetic isolation can be achieved at the scale of gyre systems in open ocean habitats (100s to 1000s of km). Furthermore, population genetic structure often varies across planktonic taxa, and appears to be linked to the particular ecological requirements of the organism. In combination, these characteristics should facilitate adaptive evolution to distinct oceanographic habitats in the plankton. We conclude that marine zooplankton may be capable of rapid evolutionary as well as ecological responses to changing ocean conditions, and discuss the implications of this view. We further suggest two priority areas for future research to test our hypothesis of high evolutionary potential in open ocean zooplankton, which will require (1) assessing how pervasive selection is in driving population divergence and (2) rigorously quantifying the spatial and temporal scales of population differentiation in the open ocean. Recent attention has focused on the ecological responses of open ocean zooplankton to current global change, including range shifts and changing phenology. Here, we argue that marine zooplankton also are well poised for evolutionary responses to global change.


Abstract
Open ocean zooplankton often have been viewed as slowly evolving species that have limited capacity to respond adaptively to changing ocean conditions. Hence, attention has focused on the ecological responses of zooplankton to current global change, including range shifts and changing phenology. Here, we argue that zooplankton also are well poised for evolutionary responses to global change. We present theoretical arguments that suggest plankton species may respond rapidly to selection on mildly beneficial mutations due to exceptionally large population size, and consider the circumstantial evidence that supports our inference that selection may be particularly important for these species. We also review all primary population genetic studies of open ocean zooplankton and show that genetic isolation can be achieved at the scale of gyre systems in open ocean habitats (100s to 1000s of km). Furthermore, population genetic structure often varies across planktonic taxa, and appears to be linked to the particular ecological requirements of the organism. In combination, these characteristics should facilitate adaptive evolution to distinct oceanographic habitats in the plankton. We conclude that marine zooplankton may be capable of rapid evolutionary as well as ecological responses to changing ocean conditions, and discuss the implications of this view. We further suggest two priority areas for future research to test our hypothesis of high evolutionary potential in open ocean zooplankton, which will require (1) assessing how pervasive selection is in driving population divergence and (2) rigorously quantifying the spatial and temporal scales of population differentiation in the open ocean.

Evolution in the Open Sea
The oceans are changing on a global scale and, in some cases, at rates greatly exceeding those observed in the historical and recent geological record (e.g., Pelejero et al. 2010). Holoplankton, the organisms that spend their entire life cycle in the open water column, are particularly good indicators of climate change (Hays et al. 2005), and show the most dramatic range shifts of any organisms reported in either terrestrial or marine environments (e.g., Beaugrand et al. 2002Beaugrand et al. , 2009Burrows et al. 2011). Marine zooplankton serve as key links in the food web between primary producers and higher trophic levels (e.g., fish, micronekton), and also are important mediators of biogeochemical fluxes in the ocean. Marine zooplankton are a phylogenetically diverse group (see Fig. 1) including representatives from 12 animal phyla (Angel 1993). Most taxa are diploid and sexual species, but some notable exceptions include members of the phyla Cnidaria and Urochordata, which alternate between asexual and sexual phases of their life cycle (e.g., scyphozoans, salps).
Although we know very little about the evolutionary potential of open ocean zooplankton, most authors have explicitly or implicitly assumed that zooplankton will show limited evolutionary responses to climate change (e.g., Helaou€ et Reygondeau and Beaugrand 2010;Stegert et al. 2010; but see Dam 2013). Hence, much attention has been focused on the observations of ecological responses to climate change in these species, for example, on their changing species' distributions and phenology (e.g., Hays et al. 2005;Richardson 2008; Ji et al. 2010). Here, we consider evolutionary responses, that is, those that result in genetic changes in a population, for example, in response to selective pressure. Although the number of studies that have rigorously tested for genetic adaptation in marine zooplankton is small, and limited to estuarine and coastal taxa that are amenable to laboratory experimentation (reviewed in Dam 2013), unequivocal evidence exists for genetic adaptation in marine zooplankton. One example is that grazer populations with a history of exposure to toxic algal blooms have significantly higher fitness when challenged with toxic prey than those with little or no such history of exposure (Colin and Dam 2007). There is also an extensive literature providing evidence for adaptive evolution in freshwater zooplankton, for example, using the water flea Daphnia as a model system (e.g., Critescu et al. 2012;Orsini et al. 2012). Our review focuses on the open ocean zooplankton, the majority of which cannot be cultured in the laboratory. We argue that evolutionary responses to global change are important to consider for these taxa, and explore indirect methods for studying evolution in these oceanic species.
Three primary arguments have traditionally supported the perception that open ocean zooplankton have limited capacity to evolve in comparison to other terrestrial and marine species (e.g., van der Spoel and Pierrot-Bults 1979;Angel 1992;Palumbi 1992;Grosberg and Cunningham 2001). First, low species diversity in the plankton has been interpreted as evidence for depressed speciation rates in comparison to other organisms. The pelagic habitat lacks obvious geographic isolating barriers that would be necessary for speciation in allopatry, the most common geographic mode of speciation (Pierrot-Bults and van der Spoel 1979;Angel 1993;van der Spoel 1994;Norris 2000;Coyne and Orr 2004). Although molecular studies are revealing cryptic species across the spectrum of pelagic animal phyla (Dawson and Jacobs 2001;Darling and Wade 2008;Morard et al. 2009;Goetze 2010a;Jennings et al. 2010;Miyamoto et al. 2010;Ortman et al. 2010), true estimates of species numbers are unlikely to be orders of magnitude higher than current estimates. A recent estimate of global marine diversity reported that~226,000 eukaryotic marine species are described, of a total of 0.7-1.0 million marine species (Appeltans et al. 2012). Of these, only a small fraction are DNA barcoded , and molecular methods are estimated to add tens of thousands, rather than hundreds of thousands, of species to the currently accepted species list (Appeltans et al. 2012). Thus, we can accept that species diversity in the plankton is relatively low. Second, marine zooplankton have among the largest (effective) population sizes of any organisms on Earth (e.g., Lynch et al. 1991;Bucklin and Wiebe 1998;Avise 2000;Goetze 2005;Peijnenburg et al. 2005), and therefore genetic drift is expected to be ineffective at changing allele frequencies within these populations. Third and finally, marine plankton are envisioned to be high-dispersal species (Ekman 1953;Angel 1993;van der Spoel 1994), with extensive migration among populations limiting their capacity for local adaptation. Plankton have been likened to airborne spores or wind-dispersed seeds that can drift almost anywhere in the ocean (Norris 2000), with their biogeographic ranges limited only by their ability to find suitable habitat for the establishment of new populations (Norris and de Vargas 2000;Sexton and Norris 2008). Collectively, these three arguments have supported a persistent view of low evolutionary potential for marine zooplankton, in comparison to other marine and terrestrial groups. However, a number of recent observations and insights suggest that we may have overlooked important processes driving the evolution of open ocean zooplankton. Here, we propose the idea that selection may be a dominant driver of marine zooplankton evolution based on theoretical insights and circumstantial evidence that selection may be widespread in these species. We also review the 46 primary population genetic and phylogeographic studies that focus on population-level differentiation in marine zooplankton (55 taxa; Table 1), and show that  genetic isolation often is observed across distinct pelagic biomes even in these entirely planktonic taxa.

Selection as a Dominant Evolutionary Force
One important misinterpretation that has supported the view of low evolutionary potential in marine zooplankton is the relative unimportance of genetic drift in influencing allele frequencies in large populations, according to standard population genetic models (Crow and Kimura 1970). However, relatively little attention has been paid to the reverse side of this coin, namely that selection is highly effective in large populations. The reasons for the higher efficacy of selection are twofold. First, more adaptive mutations occur in large populations, simply because there are more individuals to undergo mutation (e.g., Lynch et al. 1991;Barton 2010). Furthermore, beneficial mutations are common enough in large populations to arise both recurrently and on independent genetic backgrounds, increasing the chance that they arrive at the right time, place, and on the appropriate genetic background (Pennings and Hermisson 2006a,b;Ralph and Coop 2010). Second, selection is more effective because the stochastic effect of drift is smaller in large populations (e.g., Gillespie 1999Gillespie , 2000Gillespie , 2001Charlesworth 2009). A simple theoretical model (Fig. 2) demonstrates how even very mild, immeasurably small, selection pressures (s = 10 À9 to 10 À15 ) can quickly drive beneficial mutations to fixation at population sizes typical of marine zooplankton (N e = 10 7 to 10 15 ; Fig. 2). Figure 2 also shows that with increasingly large populations, the selection coefficients that have a substantial effect with respect to substitution rate become progressively smaller. Therefore, the assumption of neutrality will be more likely violated in species with large populations. Even though there is still active debate as to whether most of the genetic variation within and between species is selectively neutral (Kimura 1968(Kimura , 1983, or whether a large proportion of the variation is subject to selection (Gillespie 1991;Hahn 2008;Wares 2010), we argue here that selection may be a dominant force in the evolution of open ocean plankton, because they are at the extreme end of the scale in terms of population size of marine organisms. The effect of population size is illustrated by genome-wide studies in terrestrial model organisms that reported substantial evidence of adaptive evolution in Drosophila and Escherichia coli (large populations), whereas genetic variation in Homo sapiens and Arabidopsis (small populations) conformed to a background selection model, with a large number of segregating deleterious polymorphisms (reviewed in Hahn 2008). Barton (2010) and Karasov et al. (2010) have also argued that adaptation is not mutation limited in very large pop- ulations, and that rapid appearance of adaptive alleles can enable fast evolution, for example, adaptation to insecticide in Drosophila within 50 years, or~1000 generations. Thus far, no studies of open ocean zooplankton have directly addressed the question as to how pervasive selection is in driving evolutionary change in natural populations (but see Dam 2013 and Sanford and Kelly 2011 for selected examples from coastal zooplankton and benthic invertebrates with meroplanktonic larvae). One interesting conclusion from Sanford and Kelly (2011) was that species with planktonic dispersal comprise a surprisingly high percentage (66%) of the marine invertebrates known or suspected of exhibiting local adaptation. Several studies of pelagic marine fish have indicated a much more important role for selection than was previously thought (e.g., Hauser and Carvalho 2008;Gaggiotti et al. 2009;Bradbury et al. 2010; but see McCusker and Bentzen 2010), and selection experiments with the coccolithophore Emiliania huxleyi have shown that adaptive evolution to ocean acidification can occur within~500 generations (Lohbeck et al. 2012).
Three lines of evidence suggest that selection may be a particularly important driver of evolution in open ocean zooplankton. First, if selection works efficiently in large populations, we would expect the assumption of neutrality to be violated often in studies of marine zooplankton. Indeed, we do commonly see significant deviations from neutrality across a broad array of planktonic taxa, with 16 of 21 studies reporting nonneutral evolution based on Tajima's D, Fu and Li's, or Rozas's R 2 tests (Table 1). However, it is well known that these commonly applied neutrality tests (e.g., Tajima's D test, Tajima 1989) are sensitive to fluctuations in population size (e.g., Simonsen et al. 1995;Nielsen 2005), and significant results have often been interpreted in this context. The causative factors that underlie departures from neutrality therefore remain largely unknown, but could indicate that selection commonly influences the extent and distribution of genetic variation in marine zooplankton populations (see also Wares 2010). Broader investigation of these patterns is warranted, and we call for consistent inclusion of rigorous tests of neutrality (see e.g., Nielsen [2005] for a review) as a standard component of data analyses, such that it will be possible to assess how commonly the neutrality assumption is violated across taxa and loci. The McDonald-Kreitman (MK) test (Mcdonald and Kreitman 1991) is particularly informative in cases where multiple protein-coding sequences of related species are available, because this test is robust to demographic assumptions (Nielsen 2001). The MK test is based on the prediction of neutral theory that the ratio of replacement (nonsynonymous) to silent (synonymous) fixed differences observed between species should equal the ratio of replacement to silent polymorphisms within species. A significant excess of replacement fixed differences relative to silent changes is interpreted as evidence for adaptive evolution. To our knowledge, only one study of open ocean zooplankton has applied this test , and it reported significant evidence for selection acting on the mitochondrial cytochrome oxidase II gene in two chaetognath species.
A second line of evidence suggesting that selection may be an important driver of evolution in species with large populations is the finding of large discrepancies between census and effective population sizes, which are commonplace for marine species (Grant and Bowen 1998;Avise 2000;Hauser and Carvalho 2008;Portnoy et al. 2009). Such discrepancies would be expected under widespread deviations from the neutral model. Studies of marine zooplankton that have contrasted census and effective population sizes, as estimated from abundance and genetic data, respectively, have reported dramatic differences ranging from 10 8 to 10 10 fold in chaetognaths , copepods (Bucklin and Wiebe 1998), and krill (Zane et al. 1998; but see Zane et al. 2000, reported a threefold difference). All of these studies used  (Crow and Kimura 1970) showing that substitution rate is sensitive to small selection coefficients in large populations. For slightly beneficial mutations with selection coefficient (s), the fixation probability (P) can be approximated by: where N is diploid effective population size. When s converges to 0 (i.e., mutations are neutral) P is 1/2N and as s grows larger P becomes approximately 2s. For simplicity we assume that the substitution rate can be described as the number of mutations arising times the fixation probability. The substitution rate, relative to the neutral substitution rate, is plotted as a function of effective population size for various immeasurably small selection coefficients (ranging from 10 À9 to 10 À15 ). This model ignores clonal interference, that is, competition between beneficial mutations, which is expected to slow down the response to selection in asexual species (e.g., Gerrish and Lenski 1998).
mitochondrial DNA to estimate effective population sizes, which, if under selection, would be expected to result in lower estimates of N e relative to 'true' N e , as estimated from neutral markers. Recent work has also shown that a key prediction of the neutral theory, namely that species with large populations should have more genetic diversity than species with small populations, does not hold for mitochondrial DNA (Bazin et al. 2006), and this observation was explained by the influence of pervasive selection on the mitochondrial genome (see also Meiklejohn et al. 2007;Galtier et al. 2009;Wares 2010). Similarly for the nuclear genome, levels of genetic diversity across the tree of life do not scale with species abundances, which is inconsistent with a neutral model of evolution (Lynch 2006;Hahn 2008). Third, and finally, as the strength and type of selection will differ between unlinked genetic markers, we would expect to see large differences in the degree of population structure detected by different loci if selection is an important driver of genetic differentiation (see also Fig. 2 and Piganeau et al. 2011). Very few studies of zooplankton population structure have incorporated multiple unlinked genetic markers (Table 1), making it difficult to evaluate how common such discrepancies are. One study that contrasted mitochondrial and nuclear variation in a planktonic chaetognath reported large differences in the degree of structure detected by the two marker types (Peijnenburg et al. 2006). The historical heavy reliance on mtDNA markers in plankton population genetic studies has limited our ability to detect selection, and further tests of the ideas outlined above will become possible as the field moves toward genome-wide data and direct comparisons across loci and taxa. We do not yet have a representative view of the amounts of genomic diversity present in plankton populations and it remains to be tested whether selection models (and what types of model) or classical models of genetic drift are most appropriate to describe molecular evolution in zooplankton.

Genetic Isolation in Open Ocean Habitats
One important supporting argument for the idea that zooplankton are slowly evolving is that high dispersal rates limit their capacity to evolve adaptively in response to spatially varying selection (e.g., Grosberg and Cunningham 2001). But what is known regarding dispersal and connectivity among populations of open ocean zooplankton?
Early genetic research on the population structures of planktonic species often did support the view of nearly unlimited dispersal. Initially, workers focused on ecologically important, keystone species, and usually found genetic homogeneity among samples from distant areas of the distributional range. For example, Euphausia superba, the keystone euphausiid of the Antarctic pelagic ecosystem, has been a target species for population genetic studies dating back to the mid 1980s (e.g., Table 1; MacDonald et al. 1986;Fevolden and Scheppenheim 1989;Bortolotto et al. 2011). Although a number of studies reported weak but significant population structure (Fevolden and Ayala 1981;Zane et al. 1998;Batta-Lona et al. 2011), these patterns disappeared in more rigorous work with higher sample size and higher resolution genetic markers (Fevolden and Scheppenheim 1989;Bortolotto et al. 2011). Based on current evidence, it appears that this species is panmictic throughout its range. Note, however, that with even more powerful data (e.g., thousands of Single Nucleotide Polymorphism (SNP) loci accessed using Next Generation Sequencing technologies) it may well be that the null hypothesis of genetic homogeneity is ultimately rejected. Studies on Calanus finmarchicus, the dominant planktonic copepod in the boreal North Atlantic, followed a similar trajectory, with a number of reports of weak but significant differentiation between gyre systems and at basin spatial scales (Table 1; Bucklin and Kocher 1996;Bucklin et al. 2000a;Unal and Bucklin 2010). However, using a powerful combination of mitochondrial and nuclear microsatellite loci, Provan et al. (2009) demonstrated that the species lacks population genetic structure across the North Atlantic (Provan et al. 2009). A number of studies on other holoplanktonic organisms also have hinted at high dispersal rates, including a notable series of studies on planktonic foraminifers that found identical genotypes across different ocean basins worldwide (de Vargas et al. 1999(de Vargas et al. , 2001(de Vargas et al. , 2002(de Vargas et al. , 2004; note that isolation is also seen in this group, e.g., Aurahs et al. 2009;Ujiie et al. 2012;Seears et al. 2012). In combination, this body of work has demonstrated that there probably are some holoplanktonic species that are panmictic throughout their distributional range. However, we argue here that the population structure of these hyper-abundant species may not be typical of the broader planktonic fauna, and an early focus on these species has inappropriately instilled the sense that holoplankton have universally high dispersal capacity.
The emerging patterns from genetic results on a broader array of species suggest a more complex and nuanced picture of dispersal and connectivity among populations of open ocean zooplankton. First, there is strong evidence for spatial genetic structure and limited gene flow among populations in a number of plankton species, with distributions ranging from coastal seas to open ocean habitats. For example, planktonic chaetognaths and copepods have been shown to have restricted dispersal in the coastal ocean (Peijnenburg et al. 2004(Peijnenburg et al. , 2006Papadopoulos et al. 2005;Unal et al. 2006 Chen and Hare 2011;Yebra et al. 2011), with rare multilocus studies providing the most compelling evidence for limited gene flow among European basins (Peijnenburg et al. 2006;Reusch et al. 2010). Truly oceanic species also have been found to harbor a number of genetically very distinct populations (Goetze 2005(Goetze , 2011Papetti et al. 2005;Nelson et al. 2009), with sample sizes sufficiently high to characterize the distribution of genetic variation at mesoscale (10s to 100s of km) to ocean basin spatial scales (>1000 km). If spatial genetic structure is interpreted in light of gene flow, these studies imply a much more limited dispersal range than expected for holoplanktonic organisms in open ocean habitats (e.g., Goetze 2011). Nonmetazoan holoplankton also exhibit strong spatial genetic structure in a number of cases, implying that dispersal may be more limited than initially expected across a broad diversity of planktonic life (e.g., Darling et al. 2007;McCauley et al. 2009;Rynearson et al. 2009;Casteleyn et al. 2010;Ujiie et al. 2012;Whittaker et al. 2012).
A second emerging insight is that among zooplankton species with circumglobal or cosmopolitan biogeographic ranges, spatial genetic structure often occurs at the scale of pelagic biomes. For example, studies on the copepods Eucalanus hyalinus, E. spinifer, Pleuromamma xiphias, and Haloptilus longicornis all report genetic homogeneity among sample sites within subtropical gyres, with strong and significant population genetic breaks occurring between subtropical gyres and ocean basins (Goetze 2005(Goetze , 2011Norton and Goetze in press). These results imply limited dispersal across equatorial waters in all of these species, irrespective of whether or not the species is antitropical in distribution (Eucalanus species are antitropical, Pleuromamma and Haloptilus spp. are not). In these cases (and others), barriers to dispersal in the open sea occur at oceanographic regions with sharp transitions in both biotic and abiotic properties (Thornhill et al. 2008;Goetze 2011;Norton and Goetze in press), and are likely bio-physical in nature.
Third and finally, increasingly it appears that population genetic structure is species-specific in marine zooplankton, with closely related species often showing very different patterns of differentiation across the same pelagic habitats (Goetze 2005;Blanco-Bercial et al. 2011;Chen and Hare 2011). These results suggest that there are differences between species in realized dispersal, which may be linked to species-specific ecological or habitat requirements. For example, Chen and Hare (2011) documented salinity differences in the preferred habitat of three highly divergent lineages within Acartia tonsa that likely played a role in the development of phylogeographic structure within each of these cryptic species. The genetic lineage found in estuarine habitats showed high phylogeographic divergence among coastal embayments (US Atlantic coast), while the most oceanic lineage entirely lacked population genetic structure across the same sampling sites. Similar ecophysiological mechanisms likely control dispersal among populations of open ocean zooplankton, though there are no studies that have examined this for truly oceanic taxa. What constitutes optimal habitat will vary across species, with the nature of dispersal barriers highly dependent on the ecological niche of the organism. Therefore, it is more appropriate to consider population genetic structure a species-specific trait that varies across holoplanktonic taxa, as does any other ecological or life history trait, than to view the zooplankton as monotypically 'high dispersal' species.
This emerging view of the complexity of dispersal in marine zooplankton extends a number of observations that have been made in other marine species to the pelagic habitat. Panmixia was once expected for all marine species with large population size and a planktonic phase of the life history (e.g., Palumbi 1992Palumbi , 1994. This idea now has been discarded as overly simplistic for both marine fish and benthic invertebrates (with meroplanktonic larvae), with extensive empirical studies having shown 'closed' marine populations and adaptive responses to human-induced and natural environmental change (e.g., Conover et al. 2006;Hauser and Carvalho 2008;Cowen and Sponaugle 2009;Hellberg 2009;Marshall et al. 2010;Sanford and Kelly 2011). There is now considerable evidence that many open ocean zooplankton species also demonstrate population genetic differentiation across distinct pelagic habitats, an observation that has important implications regarding their capacity to respond to local selective forces.

New Views on Marine Zooplankton Evolution
Given the findings and arguments outlined above (on selection and genetic isolation), we hypothesize that open ocean zooplankton have large adaptive potential and may be capable of strong evolutionary responses (genetic adaptation) to changing ocean conditions. These responses could be rapid, occurring on ecological time scales (Hairston et al. 2005;Schoener 2011), and are thus important to consider in discussions about responses of the ocean's biota to contemporary climate change (Davis et al. 2005;Parmesan 2006). Significant findings of spatial population structuring in zooplankton (Table 1) imply that selection could drive fitness toward a phenotypic optimum for ocean conditions experienced within a particular oceanographic habitat (e.g., a subtropical gyre), rather than across the range of environmental conditions that occur throughout the distribution of any particular species (often cir-cumglobal in subtropical and tropical waters). This is important, as evolution can then trend toward specialization on specific ocean habitats rather than toward phenotypic plasticity or ecological generalism. In addition, marine zooplankton have short generation times, which can play a critical role in the rate of evolutionary change (e.g., Thomas et al. 2010), and have been shown capable of very rapid evolution in peripheral marine environments (Dawson and Hamner 2005) and during habitat invasions (Lee et al. 2011, and references therein). These observations suggest that zooplankton can show rapid responses to selection under appropriate conditions.
One implication of our hypothesis is that marine zooplankton may have a more limited capacity to track suitable habitat under changing ocean conditions than is currently assumed (e.g., Provan et al. 2009;Ji et al. 2010;Reygondeau and Beaugrand 2010). Results from the continuous plankton recorder (CPR) survey have convincingly demonstrated that range shifts are occurring at the northern boundary of the distribution of many temperate and boreal North Atlantic species in response to warming ocean conditions (Beaugrand et al. 2002Beaugrand 2003;. However, many of the species within the CPR record have distributions that extend across a range of pelagic biomes (e.g., subtropical gyres, boundary current, and equatorial upwelling ecosystems), and nothing is currently known regarding dispersal within and among different areas of their distribution. These species may in fact consist of a global collection of semi-isolated populations that respond individually to the changing distribution of their preferred habitat. Also, because the presence and strength of barriers to dispersal are expected to be species-specific for marine zooplankton (see above), the capacity of a species to track suitable habitat may vary across taxa and may lead to shifting community structure and trophic mismatches in pelagic foodwebs (e.g., as observed in the North Sea, Reygondeau and Beaugrand 2010).

Priorities for Future Research
Currently, very little is known about adaptation of zooplankton to local oceanographic conditions. Testing our hypothesis of high evolutionary potential in open ocean zooplankton will require research in at least two primary areas. First, the question of how pervasive natural selection is in driving population divergence in zooplankton needs to be addressed. Second, the spatial and temporal scales over which differentiation occurs in the open ocean need to be rigorously quantified. Next Generation sequencing (NGS) technologies have made possible the rapid generation of large-scale sequence data from nonmodel organisms at reasonable cost. We expect that these technologies (reviewed elsewhere in e.g., ) and in particular Restriction-site associated DNA sequencing (RADSeq, Davey and Blaxter 2011) will enable the field to move away from single marker studies (some fine examples of this approach include Emerson et al. 2010;Hohenlohe et al. 2010;Wagner et al. 2013).
Indirect approaches can be used to address the question as to how pervasive selection is in driving the evolution of open ocean zooplankton (taxa that typically cannot be cultured in the lab). One cost-effective approach to examining adaptive variation in natural populations is to study functionally important traits related to, for instance, morphology or life history. Differentiation at these quantitative traits is typically expressed as Q ST or P ST (e.g., Meril€ a and Crnokrak 2001; McKay and Latta 2002; Leinonen et al. 2008) and can be compared to patterns of differentiation at neutral genetic markers, which are typically expressed as F ST . If Q ST > F ST , diversifying selection is invoked, in which different phenotypes are favored in different populations. Whereas if Q ST < F ST , this suggests stabilizing selection, and the same phenotype is favored across populations. Another indirect approach would be to use genome scans to contrast adaptive and neutral variation across the genome (e.g., reviewed in . By comparing divergent populations, it is possible to identify loci showing signatures of selection using, for instance, outlier analyses, cline analyses, or tests of selection based on ratios of neutral and adaptive changes (reviewed in e.g., Nielsen 2005;Nielsen et al. 2009;Rice et al. 2011). A third indirect approach is the candidate gene approach, in which specific loci (genes) of known functional significance are sequenced. A nice example is provided by Larmuseau et al. (2009) in which a sensory gene, rhodopsin, was sequenced for sand goby individuals from different populations in Europe. They found that sequences of this gene did not cluster according to geographical or historical proximity, but according to the general photic conditions of the habitat of the fish. Zooplankton populations sampled along naturally occurring gradients of temperature, salinity, or acidification, which change with latitude and/or depth, could be examined for sequence variation at neutral as well as targeted genetic loci. An example of a genome-wide approach can be found in Bradbury et al. (2010), in which a survey of single nucleotide polymorphisms in Atlantic cod was used to pinpoint 40 loci for which allele frequencies showed parallel temperature-associated clines in the eastern and western north Atlantic. Note that the effect of selection is inferred indirectly using these methods and additional experiments (e.g., laboratory, mapping, and/or gene expression studies) would be necessary to link genetic loci with specific targets of selection (see Feder and Mitchell-Olds' [2003] perspective on the emerging field of evolutionary and ecological functional genomics [EEFG] which seeks to understand the functional basis of evolutionary forces shaping ecologically important traits in natural biological communities). Once several comparative genome-wide datasets become available from natural populations of planktonic taxa, we can begin to draw more general conclusions about the levels of intraspecific and intragenomic genetic variability and the relative importance of selection.
In order to rigorously quantify the spatial and temporal scales over which differentiation occurs in the open ocean, we need to substantially increase the types and amounts of data that are collected. First, it is critically important to be able to compare genetic results across marker types in order to obtain rigorous estimates of dispersal through Bayesian, Maximum Likelihood, or coalescent analytical techniques (e.g., Pritchard et al. 2000;Hey and Nielsen 2004;Zhang 2008;Beerli and Palczewski 2010;Hey 2010;Nielsen et al. 2012). The historical dominance of single marker studies has limited our capacity to distinguish demographic forces (migration, changing population size) from selection. Second, limited sampling in terms of the number of individuals sampled per population (Table 1), but also across space and time, has resulted in low scientific rigor in some studies. Revealing the amounts of genetic variation present within populations is particularly important, because this is a critical component to assessing evolutionary potential (e.g., Chevin et al. 2010). In addition, temporal variability needs to be more comprehensively addressed in future work, as nearly nothing is known regarding the stability of zooplankton population structure through time (but see Peijnenburg et al. 2006 andBortolotto et al. 2011 as rare examples of studies with time-stratified sampling). Third, and finally, resolving spatial patterns of gene flow for a greater diversity of species is also paramount to understanding the nature of dispersal barriers for marine zooplankton. Current evidence suggests that realized dispersal for holoplanktonic species may be linked to the ecological niche of the organism, and this idea needs to be tested with community-wide data on population genetic structures of a diverse panel of species with distinct pelagic niches. Ideally, such species should be collected on the same research cruises so that comparable spatial and temporal scales are sampled. Another approach would be to compare population genetic structures of related species with different life history or ecological characteristics as is common in the benthic invertebrates research community.
Finally, a wide range of questions on the fundamental basis of evolutionary change in marine zooplankton remain to be addressed. We hope that the ideas outlined here, in combination with the suite of high-throughput genomics techniques now available for development of genome-wide markers in nonmodel species will catalyze research in this field.