Figures
Abstract
Although stochasticity in oceanographic conditions is known to be an important driver of temporal genetic change in many marine species, little is known about whether genetically distinct plankton populations can persist in open ocean habitats. A prior study demonstrated significant population genetic structure among oceanic gyres in the mesopelagic copepod Haloptilus longicornis in both the Atlantic and Pacific Oceans, and we hypothesized that populations within each gyre represent distinct gene pools that persist over time. We tested this expectation through basin-scale sampling across the Atlantic Ocean in 2010 and 2012. Using both mitochondrial (mtCOII) and microsatellite markers (7 loci), we show that the genetic composition of populations was stable across two years in both the northern and southern subtropical gyres. Genetic variation in this species was partitioned among ocean gyres (FCT = 0.285, P < 0.0001 for mtCOII, FCT = 0.013, P < 0.0001 for microsatellites), suggesting strong spatial population structure, but no significant partitioning was found among sampling years. This temporal persistence of population structure across a large geographic scale was coupled with chaotic genetic patchiness at smaller spatial scales, but the magnitude of genetic differentiation was an order of magnitude lower at these smaller scales. Our results demonstrate that genetically distinct plankton populations persist over time in highly-dispersive open ocean habitats, and this is the first study to rigorously test for temporal stability of large scale population structure in the plankton.
Citation: Goetze E, Andrews KR, Peijnenburg KTCA, Portner E, Norton EL (2015) Temporal Stability of Genetic Structure in a Mesopelagic Copepod. PLoS ONE 10(8): e0136087. https://doi.org/10.1371/journal.pone.0136087
Editor: Hans G. Dam, University of Connecticut, UNITED STATES
Received: April 4, 2015; Accepted: July 29, 2015; Published: August 24, 2015
Copyright: © 2015 Goetze et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: mtCOII sequence data are available in GenBank under accession numbers KR872026—KR872295 and KC713636—KC713781. Raw microsatellite genotypes are available as Supporting Information. Oceanographic data from both Atlantic Meridional Transect cruises are available through the British Oceanographic Data Center.
Funding: This work was supported by National Science Foundation grants OCE-1029478 and OCE-1338959 to E. Goetze, NWO-VENI grant 863.08.024 to K.T.C.A. Peijnenburg, as well as National Geographic Society/Waitt grant W119-10 to E. Goetze. Plankton collections were partially supported by the UK Natural Environment Research Council National Capability funding to Plymouth Marine Laboratory and the National Oceanography Centre, Southampton, through the Atlantic Meridional Transect (AMT) Programme.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Understanding the ecology and evolution of marine species requires knowledge of the extent to which conspecific populations are genetically differentiated across space, and whether these spatial patterns are stable through time. At large geographic scales, marine species with large and stable population sizes, broad species ranges, long-lived planktonic larvae, and high fecundity are expected to exhibit temporal stability in allele frequencies coupled with weak or no genetic structure among populations (e.g., [1]), due to the combined effects of low genetic drift and high migration among sites. Although the expectation of weak population differentiation in such ‘high gene flow’ species is commonly observed (e.g., [2–6]), temporal stability in allele frequencies has been documented in very few cases [7,8]. In contrast, there is an extensive literature describing unexpectedly high levels of non-geographic differentiation among samples (e.g., [9–11]), termed chaotic genetic patchiness [12,13], due to genetic variation at fine temporal and/or spatial scales. In benthic species with meroplanktonic larvae, this fine-scale pattern appears to be largely driven by high variance in reproductive success, coupled with stochasticity in larval survivorship and transport due to variation in oceanographic conditions (e.g., [14–16]). Similar patterns also have been reported for pelagic marine fish. For example, in the European eel, results previously considered to represent spatial structuring among populations have since been suggested to reflect temporal heterogeneity [17,18], with both isolation by spawning cohorts as well as variance in reproductive success within cohorts contributing to differentiation among samples.
Marine holozooplankton, ocean drifters during their entire lifecycle, arguably represent the end of a continuum for marine animals: They are characterized by exceptionally large population size and among the greatest dispersal potential of all marine species. Because of these characteristics, we expect to observe little or no genetic differentiation among populations, and high temporal stability in allele frequencies at large spatial scales. However, despite their very high capacity for dispersal in the open sea, many holozooplankton species are strongly differentiated among populations at ocean gyre and basin spatial scales [19]. For example, populations in the northern and southern subtropical gyres of the Atlantic, Pacific and Indian Oceans exhibit strong genetic differentiation in the cosmopolitan copepods Eucalanus spinifer, Eucalanus hyalinus, and Pleuromamma xiphias [20,21], and the temperate-boreal North Atlantic species Meganyctiphanes norwegica (euphausid) and Sagitta setosa (chaetognath) also have genetically distinct populations across water masses and/or across European coastal seas [22–25]. Yet virtually nothing is known regarding the stability of spatial genetic structure in marine zooplankton, because temporally-replicated sampling remains very rare. Important exceptions include work on two krill species [25,26] and a chaetognath [27], in which temporally-replicated samples in certain areas showed that the genetic composition was stable over time. However, since sampling in these studies was not replicated over the entire spatial scale of study, it has never been explicitly tested whether the pattern of spatial structuring was temporally stable.
In this study, we examine the temporal stability of spatial population genetic structure across oceanic habitats in the mesopelagic copepod Haloptilus longicornis. This common upper-mesopelagic species has been shown in prior work to have genetically distinct populations within subtropical gyres of both the North and South Pacific and the North and South Atlantic Oceans [28,29]. Spatial patterns of abundance in the Atlantic Ocean for this species indicate that subtropical gyres are the preferred habitat, and that although the species is present in equatorial waters, this region serves as a biophysical barrier to migration among gyre populations [30]. Previous analyses based on both mtDNA and nuclear microsatellite markers indicate that the nominal species H. longicornis is composed of two morphologically cryptic species [28]. Throughout the remainder of the paper we refer to these as H. longicornis species 1 and H. longicornis species 2 (or sp. 1 & sp. 2). Both species are circumglobal in distribution, and occur sympatrically across subtropical and tropical waters. This study focuses on understanding temporal genetic patterns within H. longicornis species 1, the more abundant of the two species.
Given prior observations of significant population structure in H. longicornis, this system presents an opportunity to test whether genetically distinct populations of marine zooplankton persist in open ocean habitats. In this study, we addressed the questions: (1) Are populations of H. longicornis species 1 sampled in the northern and southern subtropical gyres of the Atlantic Ocean in 2010 genetically indistinguishable from populations found in these same ocean regions in 2012? and (2) Is the spatial population structure observed across ocean gyres temporally stable, and present during both sampling years? To answer these questions, we used material from repeat basin-scale transects in the Atlantic Ocean, and data from both mitochondrial and microsatellite markers to assess population differentiation across time and space in this system. This study is the first to rigorously test for temporal stability of population structure in oceanic zooplankton.
Materials and Methods
Plankton collections and oceanographic data
Bulk plankton samples were collected on the Atlantic Meridional Transect (AMT) cruises in 2010 and 2012 (AMT20, AMT22; Fig 1, Table 1). Permits were not required for these collections, and the work did not involve endangered or protected species. On the 2010 cruise, plankton collections were made by vertical casts of a 0.57-m diameter bongo net (200 μm, 200–0 m depths). On the 2012 cruise, a 0.71-m diameter bongo net (200 μm mesh) was towed obliquely between on average 336 m and the surface. This study includes material from 12 stations on the 2010 cruise and 15 stations on the 2012 cruise (Fig 1, Table 1), with collection locations and dates as listed in Table 1. Bulk plankton was preserved immediately in 100% ethyl alcohol, alcohol was changed within 24 hours of collection, and samples were stored at -20°C. Adult females of Haloptilus longicornis were identified to species following Bradford-Grieve [31]. Note that H. longicornis species 1 and H. longicornis species 2 cannot yet be distinguished morphologically, and both species were included in DNA sequencing and microsatellite genotyping (Table 1). DNA was extracted from individual specimens using the DNeasy Blood & Tissue Kit (QIAGEN), with modifications to the manufacturer’s protocol as described in Norton and Goetze [29].
Symbols for sites included in mtDNA and microsatellite analyses as shown in the legend (and as listed in Table 1).
MtDNA marker and sequence analysis
A 546 bp fragment of the mitochondrial gene cytochrome oxidase subunit II (mtCOII) was amplified in polymerase chain reaction (PCR) using primers COII_F6 (5’–GTC TAC AGG ATG CAA ACT CC – 3’) and COII_R9 (5’–AGA GCA TTG CCC AAA CCT GA – 3’; [29]). Amplification conditions and preparation of PCR products for sequencing were as described in Norton and Goetze [29], with Sanger sequencing of both forward and reverse strands. Forward and reverse sequences were aligned (Geneious v5.5.3), checked for errors, and unique haplotypes were identified using FaBox v1.41. MtCOII sequences were obtained from 292 and 377 animals collected in 2010 and 2012, respectively (Table 1). Of these, a total of 199 and 273 were of H. longicornis species 1 (Table 1). Animals were identified to H. longicornis species 1 based on placement in mitochondrial lineage 1 or into Microsatellite Cluster 1 by Factorial Correspondence Analysis (FCA, see below), with complete concordance observed between these two approaches in the placement of animals to species (as described in [28]). All mtCOII sequences from 2012 are first reported on here; data from 2010 were included in [29] and [28]. A sequence alignment was created including only individuals from H. longicornis species 1 (using MUSCLE, within Geneious v7.1.7, [32]. A phylogenetic tree for all haplotypes was inferred under maximum likelihood (ML) using MEGA v6.06 [33], and the Tamura and Nei (+ G) substitution model. This ML tree was converted to a haplotype genealogy and plotted in Haploviewer (http://www.cibiv.at/~greg/haploviewer).
Microsatellite markers and genotyping
Individuals were genotyped at 7 microsatellite markers that were developed for H. longicornis [28,34]: HALOMS027, HALOMS032, HALOMS064, HALOMS066, HALOMS086, HALOMS091, and HALOMS175. Primer sequences and PCR conditions were as reported in [28]. Microsatellite loci were amplified in 10 μl multiplex PCRs containing 1X Type-it Multiplex PCR Master Mix (Qiagen) and 2 μM of each primer. PCR products were genotyped using an ABI3730 Genetic Analyzer, and scoring of microsatellite chromatograms was conducted using GENEMAPPER v4.0. A total of 1058 individuals were genotyped (Table 1). Specimens were assigned to H. longicornis species 1 or H. longicornis species 2 based on factorial correspondence analysis (FCA) using GENETIX v4.05 [35], as described in [28]. A total of 840 individuals were assigned to H. longicornis species 1, and were included in all subsequent analyses (Table 1). From 28 to 60 individuals of H. longicornis species 1 were genotyped from each site, with medians of 37 and 47 individuals in 2010 and 2012, respectively. All data from 2012 are first reported here; data from 2010 were included in [28].
Population genetic analyses
For H. longicornis species 1, deviations from Hardy-Weinberg equilibrium (HWE) and linkage disequilibrium were examined using ARLEQUIN v3.5.1.3 and GENEPOP v4.2 for all microsatellite loci [36–38]. We tested for the presence of null alleles in microsatellite data using MICROCHECKER v2.2.3 [39], and estimated null allele frequencies and calculated population pairwise FST values with correction for null alleles in FreeNA [40]. Microsatellite genetic diversity indices of observed and expected heterozygosity, average alleles per locus, and allele richness were calculated in GENETIX v4.05 and FSTAT [35,41]. Pairwise FST values were calculated among all sample sites using both microsatellite and mtCOII data, as a measure of population subdivision across samples (ARLEQUIN v3.5.1.3, [38]). Significance was assessed following correction for multiple comparisons using the false discovery rate (FDR, [42,43]). Pairwise ΦST values also were calculated for the mtCOII data. We identified the nucleotide substitution model that best fit our mtCOII data using the Akaike Information Criterion, as implemented in jModelTest v2.1.4 [44], and the K81 or three-parameter model was selected as the best model (TPM3uf+G). The Tamura and Nei substitution model, which was the closest available model in Arlequin, was used to calculate pairwise and global ΦST values, and to estimate genetic diversity at each site. Hierarchical Analyses of Molecular Variance (AMOVA) based on FST were carried out to partition the genetic variance across both space (ocean gyres) and time (sampling years), for both marker types. In these analyses, we tested for population structure under the following groupings: with samples stratified by (1) northern and southern subtropical gyres (2 gyres), and (2) across two sampling years (2010, 2012). Global FST values were estimated using non-hierarchical AMOVAs among all samples, as well as among subsets of the data across ocean gyres and sampling years. Significance was tested with 10,000 permutations of genotypes or haplotypes among populations. Principal coordinate analysis (PCA) plots of linearized pairwise FST values based on both mtCOII and microsatellite data were used to visualize spatial and temporal genetic differentiation among samples. Population structure was further examined using a Bayesian clustering method implemented in STRUCTURE [45,46] for microsatellite loci. We used admixture and correlated allele frequency models, with a burn-in of 105 steps followed by 106 steps, with and without using sampling location as a prior. We ran these analyses for each of the 2010 and 2012 datasets using K = 1 to K = 10, and for the dataset of combined years using K = 1 to K = 20. We ran three separate replicates for each K to investigate consistency of Pr(X|K). The true K was evaluated by visual inspection of barplots and comparing Pr(X|K) across K values.
Results
Genetic Diversity
A total of 58 mtCOII haplotypes occurred among the 472 H. longicornis sp. 1 animals sequenced, with an average of 7 haplotypes observed at each sampling site (S1 Table). Average nucleotide diversity across all samples was 0.00553, with a range from 0.00333–0.00774. No difference was observed in the number of haplotypes, haplotype diversity, or nucleotide diversity between samples collected in 2010 and 2012 (Mann-Whitney rank sum or t-tests, P >> 0.05 in all cases). However, there were significant differences between the northern and southern subtropical gyres in the number of mtCOII haplotypes (Mann-Whitney rank sum, P = 0.006), haplotype diversity (t-test, P < 0.001), and nucleotide diversity (Mann-Whitney rank sum, P = 0.021), with higher diversity observed in the southern gyre.
Our seven microsatellite markers were moderately polymorphic, with the average number of alleles per locus ranging from 4.9 to 6.9 across all samples (S2 Table). The average observed and expected heterozygosities across all loci ranged across 0.28–0.51 and 0.44–0.58, respectively. Only one of 148 locus-by-population comparisons showed deviation from HWE after Bonferroni correction. MICROCHECKER found 4 of 7 microsatellite markers exhibiting null alleles, with null allele frequency at these markers ranging from 0.12 to 0.34. However, as noted by Andrews et al. [28] for these same microsatellite markers, when pairwise and global FST values are calculated both with and without correction for null alleles, higher values are observed in the corrected FST values. For this reason, we consider it conservative to present these values as uncorrected for null alleles, as this is least likely to identify spurious, elevated structure among populations.
Population structure
There was significant genetic differentiation among samples at both mtCOII (global FST = 0.175, P < 0.0001, global ΦST = 0.048, P = 0.0001) and microsatellite markers (global FST = 0.01, P < 0.0001, Table 2). The dominant pattern in the data was strong spatial differentiation among populations in the northern and southern subtropical gyres. In hierarchical AMOVA analyses, significant amounts of genetic variation were partitioned among the two ocean gyres (FCT = 0.285, P < 0.0001 for mtCOII, FCT = 0.013, P < 0.0001 for microsatellites), while no significant results were found among the two sampling years (Table 2). These results held for analyses including all samples, as well as for analyses with samples collected only in either 2010 or 2012 and for samples collected only within each gyre across both years. Results of non-hierarchical AMOVA analyses of both mtCOII and microsatellite data also were significant when they included samples distributed across both ocean gyres, either including both years (global mtCOII FST = 0.175, P < 0.0001, global msat FST = 0.010, P < 0.0001), 2012 only (global msat FST = 0.013, P < 0.0001, global mtCOII FST = 0.163, P < 0.0001), or 2010 only for the mtCOII data (global mtCOII FST = 0.201, P < 0.0001; Table 2). All among-gyre pairwise FST values based on mtCOII were highly significant, irrespective of sampling year (P < 0.0001; S3 Table), with values ranging up to 0.407. The majority of among-gyre pairwise ΦST comparisons also were significant, although pairwise ΦST values were lower than FST values for the same sample comparisons (S3 Table). Among-gyre values also were more often significant in microsatellite pairwise FST (46 of 99 comparisons), than were within-gyre comparisons (5 of 91 comparisons significant; S4 Table). Haplotype genealogies illustrated broad mtCOII haplotype sharing among years (Fig 2A), but with spatial separation among northern and southern gyres (Fig 2B). However, there was no deep phylogenetic division between gyres, with several haplotypes that were restricted to the North Atlantic only one or a few mutational steps away from haplotypes that were unique to or shared with the South Atlantic. The dominant mtCOII haplotypes differed between the northern and southern subtropical gyres, with haplotype 3 (blue, Fig 3A) dominant in the northern gyre and haplotype 6 (white) dominant in the southern gyre (across both years; Fig 3A). Haplotype dominance within each gyre was stable across the two sampling years. Principal coordinate analysis (PCA) of linearized pairwise FST values among all samples separated populations primarily along an axis defined by space rather than time, with 94.8% (mtCOII) and 55.3% (msat) of the variation explained by the first principal coordinate (PCo1; Fig 4). For both mtCOII and microsatellite marker types, there was no notable separation of samples by sampling year (2010, 2012). Bayesian clustering analyses using microsatellites indicated K = 1 for all datasets (2010, 2012, and both years combined); however, this type of analysis lacks power for levels of divergence comparable to our microsatellite dataset [46]. In sum, at large spatial scales, population structure among northern and southern subtropical gyres is the dominant genetic pattern.
Genealogies are plotted to illustrate the extent of haplotype sharing across (A) sampling years 2010, 2012, and (B) ocean gyres (northern and southern subtropical gyres).
(A) mtCOII haplotype frequencies within each subtropical gyre, sampled in 2010 (top) and 2012 (bottom). Each pie represents a population sample, with PopID as listed in Table 1 and shown in Fig 1 (number inside the pie). Each color is a unique haplotype; black indicates all haplotypes occurring only once in the sample. Population samples with N > 20 were included. (B) Biogeographic provinces across the Atlantic Ocean shown in section plots of (top) Chlorophyll a, and (bottom) seawater temperature (0–300 m), from the 2012 AMT22 cruise.
1 subtropical gyre populations shown in principal coordinate analysis plots, based on linearized pairwise FST values among sampling sites. (A) Results from mtCOII sequence data, shown for 15 sampling sites in 2010 and 2012. (B) Results from 7 microsatellite loci, for 20 sampling sites in 2010 and 2012. In both plots, red indicates north gyre and black indicates south gyre sampling sites, and the circle and square symbols are from 2010 and the diamond and triangle symbols are from 2012. Note the absence of separation in temporal samples within each gyre.
At the within-gyre scale, a small number of sample comparisons also were significant, suggesting chaotic genetic patchiness at fine geographic scales. One pairwise ΦST value based on mtCOII sequence data was significant in a within-gyre comparison (corrected for multiple comparisons, S3 Table). Pairwise FST values based on microsatellite data were significant at within-gyre scales in three cases, with two cases involving a comparison between years and one case involving a comparison among sites, both of which were sampled in 2012 (S4 Table). Sample sizes were high at these sites, with between 37 and 60 animals genotyped from these locations.
Discussion
Although stochastic processes are well known to underlie temporal genetic change in many marine species [1,47,48], very few studies have examined temporal variation in the marine holoplankton. Because holoplanktonic species are similar to many other marine species in being characterized by large population size, high dispersal capability, and broad geographic range, we might also expect them to exhibit temporal stability in allele frequencies and weak population genetic structure due to the combined effects of low genetic drift and high migration among sites. However, given that the pelagic habitat and the planktonic animals it contains are in constant motion, there also is uncertainty regarding whether the genetic composition of plankton in any given ocean region would be stable over time. Results reported here are the first to document that genetically distinct populations of planktonic species persist over time within subtropical gyres, even though these populations are continuously undergoing advection and diffusion in the open ocean. Our results are particularly interesting because they appear counter-intuitive given the highly dispersive nature of the pelagic habitat.
Our results suggest that temporal genetic change in zooplankton populations is heavily influenced by the spatial scale of study. At large spatial scales, we observed temporal stability in the common mesopelagic copepod Haloptilus longicornis species 1. The strongest signal in our data is of genetic differentiation between northern and southern subtropical gyres, irrespective of sampling year. The significant and strong spatial sub-structuring observed among northern and southern gyres in this species in 2010 was also present in 2012, with the same haplotypes and alleles dominant in each gyre across the two years. Samples collected in 2012 were not consistently differentiated from those collected in the same ocean gyre in 2010, and there was no significant partitioning of genetic variation among years in this species (hierarchical AMOVA, Table 2). The two-year time span between samples likely corresponds to on the order of 25–30 generations for this species [49], given the seawater temperature in their core subtropical distribution, which should be sufficient to detect stochastic variation in the genetic composition of cohorts or generations. Similar results of an absence of differentiation among samples were found in two other studies examining the temporal component of population structure in marine zooplankton [25–27]. In contrast to H. longicornis, these species have temperate-boreal (Meganyctiphanes norvegica, Sagitta setosa) or circum-antarctic (Euphausia superba) distribution patterns and temporal stability was only examined for a subset of geographic samples. The shared observation of genetically distinct populations that persist through time in the open sea suggests that temporal stability in large-scale population structure may not be uncommon for holozooplankton in a variety of habitats. Our results support our initial expectations that plankton populations within subtropical gyres represent distinct gene pools that persist over time.
Our results, along with previous studies, also indicate that temporal stability of population genetic structure at large spatial scales is coupled with chaotic genetic patchiness, or weak and ephemeral genetic differentiation, at fine spatial scales. Marine zooplankton have highly patchy abundance and distribution at the submesoscale (1–10 km; [50,51]), and chaotic genetic patchiness may occur on spatial scales of meters to kilometers, or temporal scales of days to months. Statistical heterogeneity in the genetic composition of zooplankton samples on the mesoscale and submesoscale has been observed in a number of prior studies, although the interpretation of this pattern has been varied. Bortolotto et al. [26] included temporal replicates of krill samples from the Ross Sea (5 years) and South Georgia (2 years) and concluded that although there was a single panmictic population of krill on a large geographic scale (Southern Ocean), there was also evidence for subtle, though not statistically significant, spatial and temporal heterogeneity. Other examples include spatial genetic patchiness in the copepod Metridia pacifica populations across mesoscale coastal upwelling filaments and eddies in the California Current [52,53], as well as outlier samples that show high pairwise FST values unrelated to geographic distance or ocean habitat among sites in several other studies (Eucalanus hyalinus, Pleuromamma xiphias, [20,21]). A number of other studies report unexplained genetic heterogeneity at various spatial or temporal scales in zooplankton species (e.g., [54–59]); however, in some cases, these observations were confounded by sampling artifacts and low sample size (spurious, significant FST). Genetic patchiness also can be seen in the data reported here for Haloptilus longicornis. There are population samples within both the northern and southern gyres that show weak but significant differentiation relative to other samples within the same gyre (both marker types, S3 & S4 Tables). As in other marine species, this patchiness can include either a temporal or spatial component (or both; S3 & S4 Tables). There is little understanding of the mechanisms that underlie this small-scale genetic patchiness in marine zooplankton, but in addition to variance in reproductive success (and clonal reproduction; [60]), they may include biological or bio-physical aggregation of animals (e.g., swarms, thin-layers, vertical migration; [61,62]), and selective mortality across the life history. In sum, at smaller spatial scales, zooplankton may exhibit genetic patchiness in either temporal or spatial domains, as is well-described in many other highly-dispersing marine species (e.g., [10,17,18,63–65]).
Operationally, our results suggest that the broad-scale genetic patterns of population structure within holozooplankton species can be captured with limited temporal sampling given sufficient sampling intensity (sample size). This result is reassuring, as the majority of prior studies have used samples from a wide range of sampling dates, and interpreted the patterns observed exclusively in a spatial context (e.g., [20,21,28,30,66–71]). Given spatially extensive, in many cases circumglobal, range distributions for marine zooplankton, this approach has been necessary to collecting material across the distributional range of the species. However, it has always been recognized that the lack of temporally-stratified sampling resulted in poor understanding of the extent to which our interpretation of population structure in marine holoplankton was confounded by genetic variation across both space and time (e.g., [25–27,72]). Results reported here are therefore important in documenting that genetic structure occurring over large spatial scales is temporally persistent (Table 2). If this is true for other species, then combining temporally-stratified samples to examine spatial structure is unlikely to impede interpretation of the larger-scale population structure in these species.
In summary, this study is the first to demonstrate that both the genetic composition of plankton populations within oceanic gyres and the spatial patterns of population sub-structure across ocean gyres are persistent through time. This large-scale pattern of significant and temporally-stable genetic structure for marine zooplankton is coupled with chaotic genetic patchiness at finer spatial scales (in this case, within subtropical gyres). Our results demonstrate temporal stability in the genetic structure of a ‘high dispersal’ marine species, and are particularly important given the dispersive nature of the pelagic habitat and very high vagility of planktonic species. Additional work is needed to understand the extent and underlying causes of fine-scale genetic patchiness in the plankton, as well as to interpret its biological impact on species. Confirmation of the observations of temporally stable genetic structure over large spatial scales from other zooplankton species also would lend greater support to the conclusion that these patterns are the norm for the holozooplankton.
Supporting Information
S1 File. Microsatellite genotypes for all specimens, with sample ID and cruise and station as listed for each specimen.
Please contact the corresponding author regarding subsequent use of these data (egoetze@hawaii.edu).
https://doi.org/10.1371/journal.pone.0136087.s001
(CSV)
S1 Table. Mitochondrial cytochrome oxidase subunit II (mtCOII) summary statistics and diversity indices for population samples, including only specimens of Haloptilus longicornis species 1.
Collection sites for plankton samples included in this study, from (A) the 2010 AMT cruise (AMT20), and (B) the 2012 AMT cruise (AMT22). Pop ID = the population identifier referred to throughout the manuscript; Station = the cruise and station number of each sample; N = the number of adult females included; H = the number of haplotypes; h = haplotype diversity; π = nucleotide diversity.
https://doi.org/10.1371/journal.pone.0136087.s002
(PDF)
S2 Table. Microsatellite diversity indices for population samples, including only specimens of Haloptilus longicornis species 1.
Collection sites for plankton samples included in this study, from (A) the 2010 AMT cruise (AMT20), and (B) the 2012 AMT cruise (AMT22). Pop ID = the population identifier referred to throughout the manuscript; Station = the cruise and station number of each sample; N = the number of adult females included. Ho and He are the observed and expected heterozygosity averaged across all microsatellite loci for each population. k = average number of alleles per locus.
https://doi.org/10.1371/journal.pone.0136087.s003
(PDF)
S3 Table. Pairwise FST (below diagonal) and ΦST (above diagonal) values between sample sites in the Atlantic Ocean for Haloptilus longicornis sp 1, based on mtCOII data.
Samples are grouped by ocean gyres, with sample site numbers as in Tables 1 & 2. Material from both AMT20 (2010) and AMT22 (2012) are included. Significant values (P < 0.05) are shaded grey, bold indicates significance following correction for multiple comparisons (FDR).
https://doi.org/10.1371/journal.pone.0136087.s004
(PDF)
S4 Table. Pairwise FST values between sample sites in the Atlantic Ocean for Haloptilus longicornis sp 1, based on microsatellite data.
Samples are grouped by ocean gyres, with sample site numbers as in Tables 1 & 2. Material from both AMT20 (2010) and AMT22 (2012) are included. Significant values (P < 0.05) are shaded grey, bold indicates significance following correction for multiple comparisons (FDR).
https://doi.org/10.1371/journal.pone.0136087.s005
(PDF)
Acknowledgments
We thank R. Harmer for field assistance in 2010, and L. Van Woudenberg for help in the laboratory. We thank M. Iacchei and P. Meirmans for insightful comments and suggestions. Three anonymous reviewers provided helpful comments that improved the manuscript. This is contribution # 263 of the Atlantic Meridional Transect (AMT) Programme.
Author Contributions
Conceived and designed the experiments: EG KRA. Performed the experiments: EG KRA EP ELN. Analyzed the data: EG KRA KTCAP. Contributed reagents/materials/analysis tools: EG. Wrote the paper: EG KRA KTCAP EP ELN.
References
- 1.
Toonen R, Grosberg RK (2011) Causes of chaos: spatial and temporal genetic heterogeneity in the intertidal anomuran crab Petrolisthes cintipes. In: Koenemann S, Held C, Schubart C, editors. Phylogeography and Population Genetics in Crustacea: CRC Press. pp. 75–107.
- 2. Provan J, Beatty GE, Keating SL, Maggs CA, Savidge G (2009) High dispersal potential has maintained long-term population stability in the North Atlantic copepod Calanus finmarchicus. Proceedings of the Royal Society B-Biological Sciences 276: 301–307.
- 3. Cote CL, Gagnaire PA, Bourret V, Verreault G, Castonguay M, Bernatchez L (2013) Population genetics of the American eel (Anguilla rostrata): FST = 0 and North Atlantic Oscillation effects on demographic fluctuations of a panmictic species. Molecular Ecology 22: 1763–1776. pmid:23216918
- 4. Horne JB, van Herwerden L (2013) Long-term panmixia in a cosmopolitan Indo-Pacific coral reef fish and a nebulous genetic boundary with its broadly sympatric sister species. Journal of Evolutionary Biology 26: 783–799. pmid:23305496
- 5. Albrecht GT, Valentin AE, Hundertmark KJ, Hardy SM (2014) Panmixia in Alaskan Populations of the Snow Crab Chionoecetes Opilio (Malacostraca: Decapoda) in the Bering, Chukchi, and Beaufort Seas. Journal of Crustacean Biology 34: 31–39.
- 6. Reisser CMO, Bell JJ, Gardner JPA (2014) Correlation between pelagic larval duration and realised dispersal: long-distance genetic connectivity between northern New Zealand and the Kermadec Islands archipelago. Marine Biology 161: 297–312.
- 7. Hedgecock D, Barber PH, Edmands S (2007) Genetic approaches to measuring connectivity. Oceanography 20: 70–79.
- 8. Burton RS (1997) Genetic evidence for long term persistence of marine invertebrate populations in an ephemeral environment. Evolution 51: 993–998.
- 9. Pujolar JM, Maes GE, Volckaert FAM (2006) Genetic patchiness among recruits in the European eel Anguilla anguilla. Mar Ecol Prog Ser 307: 209–217.
- 10. Lee HJ, Boulding EG (2009) Spatial and temporal population genetic structure of four northeastern Pacific littorinid gastropods: the effect of mode of larval development on variation at one mitochondrial and two nuclear markers. Molecular Ecology 18: 2165–2184. pmid:19344352
- 11. Hogan JD, Thiessen RJ, Heath DD (2010) Variability in connectivity indicated by chaotic genetic patchiness within and among populations of a marine fish. Marine Ecology Progress Series 417: 263–U289.
- 12. Johnson MS, Black R (1982) Chaotic genetic patchiness in an intertidal limpet, Siphonaria sp. Marine Biology 70: 157–164.
- 13. Johnson MS, Black R (1984) Pattern beneath the chaos: the effect of recruitment on genetic patchiness in an intertidal limpet. Evolution 38: 1371–1383.
- 14. Broquet T, Viard F, Yearsley JM (2013) Genetic Drift and Collective Dispersal Can Result in Chaotic Genetic Patchiness. Evolution 67: 1660–1675. pmid:23730760
- 15. Iacchei M, Ben-Horin T, Selkoe KA, Bird CE, Garcia-Rodriguez FJ, Toonen RJ (2013) Combined analyses of kinship and FST suggest potential drivers of chaotic genetic patchiness in high gene-flow populations. Molecular Ecology 22: 3476–3494. pmid:23802550
- 16. Selkoe KA, Watson JR, White C, Ben Horin T, Iacchei M, Mitarai S, et al. (2010) Taking the chaos out of genetic patchiness: seascape genetics reveals ecological and oceanographic drivers of genetic patterns in three temperate reef species. Molecular Ecology 19: 3708–3726. pmid:20723063
- 17. Dannewitz J, Maes GE, Johansson L, Wickstrom H, Volckaert FAM, Jarvi T (2005) Panmixia in the European eel: a matter of time… Proceedings of the Royal Society B-Biological Sciences 272: 1129–1137.
- 18. Maes GE, Pujolar JM, Hellemans B, Volckaert FAM (2006) Evidence for isolation by time in the European eel (Anguilla anguilla L.). Molecular Ecology 15: 2095–2107. pmid:16780427
- 19. Peijnenburg KTCA, Goetze E (2013) High evolutionary potential of marine zooplankton. Ecology & Evolution 3: 2765–2781.
- 20. Goetze E (2005) Global population genetic structure and biogeography of the oceanic copepods, Eucalanus hyalinus and E. spinifer. Evolution 59: 2378–2398. pmid:16396179
- 21. Goetze E (2011) Population differentiation in the open sea: Insights from the pelagic copepod Pleuromamma xiphias. Integrative and Comparative Biology 51: 580–597. pmid:21940778
- 22. Peijnenburg KTCA, Fauvelot C, Breeuwer JAJ, Menken SBJ (2006) Spatial and temporal genetic structure of the planktonic Sagitta setosa (Chaetognatha) in European seas as revealed by mitochondrial and nuclear DNA markers. Molecular Ecology 15: 3319–3338. pmid:16968273
- 23. Peijnenburg KTCA, Breeuwer JAJ, Pierrot-Bults AC, Menken SBJ (2004) Phylogeography of the planktonic chaetognath Sagitta setosa reveals isolation in European seas. Evolution 58: 1472–1487. pmid:15341150
- 24. Patarnello T, Papetti C, Zane L (2010) Genetics of northern krill (Meganyctiphanes norvergica Sars). Advances in Marine Biology 57: 41–57. pmid:20955888
- 25. Papetti C, Zane L, Bortolotto E, Bucklin A, Patarnello T (2005) Genetic differentiation and local temporal stability of population structure in the euphausiid Meganyctiphanes norvegica. Mar Ecol Prog Ser 289: 225–235.
- 26. Bortolotto E, Bucklin A, Mezzavilla M, Zane L, Patarnello T (2011) Gone with the currents: lack of genetic differentiation at the circum-continental scale in Antarctic krill Euphausia superba. BMC Genetics 12: 32. pmid:21486439
- 27. Peijnenburg KTCA, Fauvelot C, Breeuwer AJ, Menken SBJ (2006) Spatial and temporal genetic structure of the planktonic Sagitta setosa (Chaetognatha) in European seas as revealed by mitochondrial and nuclear DNA markers. Molecular Ecology 15: 3319–3338. pmid:16968273
- 28. Andrews KR, Norton E, Fernandez-Silva I, Portner E, Goetze E (2014) Multilocus evidence for globally-distributed cryptic species and distinct populations across ocean gyres in a mesopelagic copepod. Molecular Ecology 23: 5462–5479. pmid:25283587
- 29. Norton EL, Goetze E (2013) Equatorial dispersal barriers and limited population connectivity among oceans in a planktonic copepod. Limnology and Oceanography 58: 1581–1596.
- 30.
Norton E (2013) Empirical and biophysical modeling studies of dispersal barriers for marine plankton [Masters]. Honolulu, Hawaii: University of Hawaii at Manoa. 105 p.
- 31.
Bradford-Grieve J (1999) The Marine Fauna of New Zealand: Pelagic Calanoid Copepoda: Bathypontiidae, Arietellidae, Augaptilidae, Heterorhabdidae, Lucicutiidae, Metridinidae, Phyllopodidae, Centropagidae, Pseudodiaptomidae, Temoridae, Candaciidae, Pontellidae, Sulcanidae, Acartiidae, Tortanidae; Gordon DP, editor. Wellington, New Zealand: National Institute of Water and Atmospheric Research. 268 p.
- 32. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research 32: 1792–1797. pmid:15034147
- 33. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Molecular Biology and Evolution 30: 2725–2729. pmid:24132122
- 34. Fernandez-Silva I, Whitney J, Wainwright B, Andrews KR, Ylitalo-Ward H, Bowen B, et al. (2013) Microsatellites for Next-Generation Ecologists: A Post-Sequencing Bioinformatics Pipeline. PloS One 8.
- 35.
Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F (1996–2004) GENETIX 4.05, logiciel sous Windows pour la génétique des populations. Montpellier, France: Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II.
- 36. Raymond M, Rousset F (1995) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Heredity 86: 248–249.
- 37. Excoffier L, Laval LG, Schneider S (2005) Arlequin ver 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online 1: 47–50.
- 38. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources 10: 564–567. pmid:21565059
- 39. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P (2004) MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes 4: 535–538.
- 40. Chapuis MP, Estoup A (2007) Microsatellite null alleles and estimation of population differentiation. Molecular Biology and Evolution 24: 621–631. pmid:17150975
- 41. Goudet J (1995) FSTAT (Version 1.2): A Computer Program to Calculate F-Statistics. J Hered 86: 485–486.
- 42. Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate—a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological 57: 289–300.
- 43. Narum SR (2006) Beyond Bonferroni: Less conservative analyses for conservation genetics. Conservation Genetics 7: 783–787.
- 44. Darriba D, Taboada GL, Doallo R, Posada D (2012) jModelTest 2: more models, new heuristics and parallel computing. Nature Methods 9: 772–772.
- 45. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959. pmid:10835412
- 46. Hubisz MJ, Falush D, Stephens M, Pritchard JK (2009) Inferring weak population structure with the assistance of sample group information. Molecular Ecology Resources 9: 1322–1332. pmid:21564903
- 47. Barshis DJ, Sotka EE, Kelly RP, Sivasundar A, Menge BA, Barth JA, et al. (2011) Coastal upwelling is linked to temporal genetic variability in the acorn barnacle Balanus glandula. Mar Ecol Prog Ser 439: 139–150.
- 48. Selkoe KA, Gaines SD, Caselle JE, Warner RR (2006) Current shifts and kin aggregation explain genetic patchiness in fish recruits. Ecology 87: 3082–3094. pmid:17249233
- 49.
Mauchline J (1998) The biology of calanoid copepods. San Diego: Academic Press. vii, 710 p.
- 50. Folt CL, Burns JM (1999) Biological drivers of zooplankton patchiness. Trends in Ecology & Evolution 14: 300–305.
- 51.
Haury LR, McGowan JA, Weibe PH (1978) Patterns and processes in the time-space scales of plankton distributions. In: Steele JH, editor. Spatial Pattern in Plankton Communities. Plenum, N.Y. pp. 277–327.
- 52. Bucklin A, Rienecker MM, Mooers CNK (1989) Genetic tracers of zooplankton transport in cosatal filaments off northern California. Journal of Geophysical Research 94: 8277–8288.
- 53. Bucklin A (1991) Population genetic responses of the planktonic copepod Metridia pacifica to a coastal eddy in the California Current. Journal of Geophysical Research 96: 14799–14808.
- 54. Bucklin A, Smolenack SB, Bentley AM, Wiebe P (1997) Gene flow patterns of the Euphausiid, Meganyctiphanes norvegica, in the NW Atlantic based on mtDNA sequences for cytochrome b and cytochrome oxidase I. Journal of Plankton Research 19: 1763–1781.
- 55. Bucklin A, Wiebe PH, Smolenack SB, Copley NJ, Clarke ME (2002) Integrated biochemical, molecular genetic, and bioacoustical analysis of mesoscale variability of the euphausiid Nematoscelis difficilis in the California Current. Deep-Sea Research Part I-Oceanographic Research Papers 49: 437–462.
- 56. Bucklin A, Wiebe PH (1986) Genetic Heterogeneity in Euphausiid Populations—Euphausia krohnii and Nematoscelis megalops in North-Atlantic Slope Water. Limnology and Oceanography 31: 1346–1352.
- 57. Thuesen EV, Numachi K, Nemoto T (1993) Genetic Variation in the Planktonic Chaetognaths Parasagitta elegans and Eukrohnia hamata. Marine Ecology-Progress Series 101: 243–251.
- 58. Jarman SN, Elliott NG, McMinn A (2002) Genetic differentiation in the Antarctic coastal krill Euphausia crystallorophias. Heredity 88: 280–287. pmid:11920136
- 59. Jarman SN, Nicol S (2002) Sources of variance in studies of krill population genetics. CCAMLR Science 9: 107–116.
- 60. Aglieri G, Papetti C, Zane L, Milisenda G, Boero F, Piraino S (2014) Genetic patchiness in the holoplanktonic jellyfish Pelagia noctiluca (Scyphozoa, Cnidaria). PLoS One 9.
- 61. Benoit-Bird KJ, McManus MA (2012) Bottom-up regulation of a pelagic community through spatial aggregations. Biology Letters 8: 813–816. pmid:22552636
- 62. Ambler JW (2002) Zooplankton swarms: charactersitics, proximal cues and proposed advantages. Hydrobiologia 480: 155–164.
- 63. Muths D, Jollivet D, Gentil F, Davoult D (2009) Large-scale genetic patchiness among NE Atlantic populations of the brittle star Ophiothrix fragilis. Aquatic Biology 5: 117–132.
- 64. Lee HJ, Boulding EG (2007) Mitochondrial DNA variation in space and time in the northeastern Pacific gastropod, Littorina keenae. Molecular Ecology 16: 3084–3103. pmid:17651189
- 65. Hedgecock D, Pudovkin AI (2011) Sweepstakes reproductive success in highly fecund marine fish and shellfish: A review and commentary. Bulletin of Marine Science 87: 971–1002.
- 66. Yebra L, Bonnet D, Harris RP, Lindeque PK, Peijnenburg KTCA (2011) Barriers in the pelagic: population structuring of Calanus helgolandicus and C. euxinus in European waters. Marine Ecology Progress Series 428: 135–149.
- 67. Unal E, Frost BW, Armbrust EV, Kideys AE (2006) Phylogeography of Calanus helgolandicus and the Black Sea copepod Calanus euxinus, with notes on Pseudocalanus elongatus (Copepoda: Calanoida). Deep-Sea Research II 53: 1961–1975.
- 68. Peijnenburg KTCA, van Haastrecht EK, Fauvelot C (2005) Present-day genetic composition suggests contrasting demographic histories of two dominant chaetognaths of the North-East Atlantic, Sagitta elegans and S setosa. Marine Biology 147: 1279–1289.
- 69. Nelson RJ, Carmack EC, McLaughlin FA, Cooper GA (2009) Penetration of Pacific zooplankton into the western Arctic Ocean tracked with molecular population genetics. Mar Ecol Prog Ser 381: 129–138.
- 70. Chen G, Hare MP (2011) Cryptic diversity and comparative phylogeography of the estuarine copepod Acartia tonsa on the US Atlantic coast. Molecular Ecology 20: 2425–2441. pmid:21521392
- 71. Goetze E (2010) Integrated molecular and morphological biogeography of the calanoid copepod family Eucalanidae. Deep-Sea Research II 57: 2110–2129.
- 72. Batta-Lona PG, Bucklin A, Weibe PH, Patarnello T, Copley NJ (2011) Population genetic variation of the Southern Ocean krill, Euphausia superba, in the Western Antarctic Peninsula region based on mitochondrial single nucleotide polymorphisms (SNP). Deep-Sea Research II 58: 1652–1661.