Population differentiation in the open sea: insights from the pelagic copepod Pleuromamma xiphias

Revision of the 'Acanthephyra purpurea' species complex (Crustacea: Decapoda), with an emphasis on species diversification in the Atlantic

We inventoried all nine species of the 'Acanthephyra purpurea' complex, one of the most abundant and cosmopolitan group of mesopelagic shrimps. We used 119 specimens at hand and genetic data for 124 specimens from GenBank and BOLD. Phylogenetic analysis of four genes (COI, 16S, NaK, and enolase) showed that the 'Acanthephyra purpurea' complex is polyphyletic and encompasses two species groups, 'A. purpurea' (mostly Atlantic) and 'A. smithi' (Indo-West Pacific). The 'A. purpurea' species group consists of two major molecular clades A. pelagica and A. kingsleyi - A. purpurea - A. quadrispinosa. Molecular data suggest that hitherto accepted species A. acanthitelsonis, A. pelagica, and A. sica should be considered as synonyms. The Atlantic is inhabited by at least two cryptic genetic lineages of A. pelagica and A. quadrispinosa. Morphological analyses of qualitative and quantitative (900 measurements) characters resulted in a tabular key to species and in a finding of four evolutionary traits. Atlantic species showed various scenarios of diversification visible on mitochondrial gene level, nuclear gene level, and morphological level. We recorded and discussed similar phylogeographic trends in diversification and in distribution of genetic lineages within two different clades: A. pelagica and A. kingsleyi - A. purpurea - A. quadrispinosa.

Morphological and taxonomic diversity of mesozooplankton is an important driver of carbon export fluxes in the ocean

Mesozooplankton is a very diverse group of small animals ranging in size from 0.2 to 20 mm not able to swim against ocean currents. It is a key component of pelagic ecosystems through its roles in the trophic networks and the biological carbon pump. Traditionally studied through microscopes, recent methods have been however developed to rapidly acquire large amounts of data (morphological, molecular) at the individual scale, making it possible to study mesozooplankton using a trait-based approach. Here, combining quantitative imaging with metabarcoding time-series data obtained in the Sargasso Sea at the Bermuda Atlantic Time-series Study (BATS) site, we showed that organisms' transparency might be an important trait to also consider regarding mesozooplankton impact on carbon export, contrary to the common assumption that just size is the master trait directing most mesozooplankton-linked processes. Three distinct communities were defined based on taxonomic composition, and succeeded one another throughout the study period, with changing levels of transparency among the community. A co-occurrences' network was built from metabarcoding data revealing six groups of taxa. These were related to changes in the functioning of the ecosystem and/or in the community's morphology. The importance of Diel Vertical Migration at BATS was confirmed by the existence of a group made of taxa known to be strong migrators. Finally, we assessed if metabarcoding can provide a quantitative approach to biomass and/or abundance of certain taxa. Knowing more about mesozooplankton diversity and its impact on ecosystem functioning would allow to better represent them in biogeochemical models.

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

Reevaluating claims of ecological speciation in Halichoeres bivittatus

Allopatry has traditionally been viewed as the primary driver of speciation in marine taxa, but the geography of the marine environment and the larval dispersal capabilities of many marine organisms render this view somewhat questionable. In marine fishes, one of the earliest and most highly cited empirical examples of ecological speciation with gene flow is the slippery dick wrasse, Halichoeres bivittatus. Evidence for this cryptic or incipient speciation event was primarily in the form of a deep divergence in a single mitochondrial locus between the northern and southern Gulf of Mexico, combined with a finding that these two haplotypes were associated with different habitat types ("tropical" vs. "subtropical") in the Florida Keys and Bermuda, where they overlap. Here, we examine habitat assortment in the Florida Keys using a broader sampling of populations and habitat types than were available for the original study. We find no evidence to support the claim that haplotype frequencies differ between habitat types, and little evidence to support any differences between populations in the Keys. These results undermine claims of ecological speciation with gene flow in Halichoeres bivittatus. Future claims of this type should be supported by multiple lines of evidence that illuminate potential mechanisms and allow researchers to rule out alternative explanations for spatial patterns of genetic differences.

Genetic diversity and novel lineages in the cosmopolitan copepod Pleuromamma abdominalis in the Southeast Pacific

Across boundary currents, zooplankton are subject to strong oceanographic gradients and hence strong selective pressures. How such gradients interact with the speciation process of pelagic organisms is still poorly understood in the open ocean realm. Here we report on genetic diversity within the pelagic copepod Pleuromamma abdominalis in the poorly known Southeast Pacific region, with samples spanning an ocean gradient from coastal upwelling to the oligotrophic South Pacific Subtropical Gyre. We assessed variation in fragments of the mitochondrial (mt) genes cytochrome c oxidase subunit I (COI) and Cytochrome b as well as in the nuclear internal transcribed spacer (ITS) region and 28 S rRNA. Phylogenetic analyses revealed the presence of 8 divergent lineages occurring across the gradient with genetic distances in the range of 0.036 and 0.44 (mt genes), and GMYC species delimitation methods support their inference as distinct (undescribed) species. Genetic lineages occurring across the zonal gradient showed strong genetic structuring, with the presence of at least two new lineages within the coastal upwelling zone, revealing an unexpectedly high level of endemism within the Humboldt Current System. Multivariate analyses found strong correlation between genetic variation and surface chlorophyll-a and salinity, suggesting an important role for hydrographic gradients in maintaining genetic diversity. However, the presence of cryptic lineages within the upwelling zone cannot be easily accounted for by environmental heterogeneity and poses challenging questions for understanding the speciation process for oceanic zooplankton.

Assessing species boundaries in the open sea: an integrative taxonomic approach to the pteropod genus Diacavolinia

To track changes in pelagic biodiversity in response to climate change, it is essential to accurately define species boundaries. Shelled pteropods are a group of holoplanktonic gastropods that have been proposed as bio-indicators because of their vulnerability to ocean acidification. A particularly suitable, yet challenging group for integrative taxonomy is the pteropod genus Diacavolinia, which has a circumglobal distribution and is the most species-rich pteropod genus, with 24 described species. We assessed species boundaries in this genus, with inferences based on geometric morphometric analyses of shell-shape variation, genetic (cytochrome c oxidase subunit I, 28S rDNA sequences) and geographic data. We found support for a total of 13 species worldwide, with observations of 706 museum and 263 freshly collected specimens across a global collection of material, including holo‐ and paratype specimens for 14 species. In the Atlantic Ocean, two species are well supported, in contrast to the eight currently described, and in the Indo‐Pacific we found a maximum of 11 species, partially merging 13 of the described species. Distributions of these revised species are congruent with well-known biogeographic provinces. Combining varied datasets in an integrative framework may be suitable for many diverse taxa and is an important first step to predicting species-specific responses to global change.

Regional heterogeneity impacts gene expression in the subarctic zooplankter Neocalanus flemingeri in the northern Gulf of Alaska

Marine pelagic species are being increasingly challenged by environmental change. Their ability to persist will depend on their capacity for physiological acclimatization. Little is known about limits of physiological plasticity in key species at the base of the food web. Here we investigate the capacity for acclimatization in the copepod Neocalanus flemingeri, which inhabits the Gulf of Alaska, a heterogeneous and highly seasonal environment. RNA-Seq analysis of field-collected pre-adults identified large regional differences in expression of genes involved in metabolic and developmental processes and response to stressors. We found that lipid synthesis genes were up-regulated in individuals from Prince William Sound and down-regulated in the Gulf of Alaska. Up-regulation of lipid catabolic genes in offshore individuals suggests they are experiencing nutritional deficits. The expression differences demonstrate physiological plasticity in response to a steep gradient in food availability. Our transcriptional analysis reveals mechanisms of acclimatization that likely contribute to the observed resilience of this population.

Large-scale ocean connectivity and planktonic body size

Global patterns of planktonic diversity are mainly determined by the dispersal of propagules with ocean currents. However, the role that abundance and body size play in determining spatial patterns of diversity remains unclear. Here we analyse spatial community structure - β-diversity - for several planktonic and nektonic organisms from prokaryotes to small mesopelagic fishes collected during the Malaspina 2010 Expedition. β-diversity was compared to surface ocean transit times derived from a global circulation model, revealing a significant negative relationship that is stronger than environmental differences. Estimated dispersal scales for different groups show a negative correlation with body size, where less abundant large-bodied communities have significantly shorter dispersal scales and larger species spatial turnover rates than more abundant small-bodied plankton. Our results confirm that the dispersal scale of planktonic and micro-nektonic organisms is determined by local abundance, which scales with body size, ultimately setting global spatial patterns of diversity.

Body size is hypothesised to be a major determinant of β-diversity in passively-dispersing marine organisms. Here, Villarino et al. show that plankton body size determines rates of dispersal along marine currents, with shorter dispersal and higher species spatial turnover in larger organisms.

Comparative phylogeography of the ocean planet

Understanding how geography, oceanography, and climate have ultimately shaped marine biodiversity requires aligning the distributions of genetic diversity across multiple taxa. Here, we examine phylogeographic partitions in the sea against a backdrop of biogeographic provinces defined by taxonomy, endemism, and species composition. The taxonomic identities used to define biogeographic provinces are routinely accompanied by diagnostic genetic differences between sister species, indicating interspecific concordance between biogeography and phylogeography. In cases where individual species are distributed across two or more biogeographic provinces, shifts in genotype frequencies often align with biogeographic boundaries, providing intraspecific concordance between biogeography and phylogeography. Here, we provide examples of comparative phylogeography from (i) tropical seas that host the highest marine biodiversity, (ii) temperate seas with high productivity but volatile coastlines, (iii) migratory marine fauna, and (iv) plankton that are the most abundant eukaryotes on earth. Tropical and temperate zones both show impacts of glacial cycles, the former primarily through changing sea levels, and the latter through coastal habitat disruption. The general concordance between biogeography and phylogeography indicates that the population-level genetic divergences observed between provinces are a starting point for macroevolutionary divergences between species. However, isolation between provinces does not account for all marine biodiversity; the remainder arises through alternative pathways, such as ecological speciation and parapatric (semiisolated) divergences within provinces and biodiversity hotspots.

Box Jellyfish Alatina alata Has a Circumtropical Distribution

Species of the box jellyfish (Cubozoa) genus Alatina are notorious for their sting along the beaches of several localities of the Atlantic and Pacific. These species include Alatina alata on the Caribbean Island of Bonaire (the Netherlands), A. moseri in Hawaii, and A. mordens in Australia. Most cubozoans inhabit coastal waters, but Alatina is unusual in that specimens have also been collected in the open ocean at great depths. Alatina is notable in that populations form monthly aggregations for spermcast mating in conjunction with the lunar cycle. Nominal species are difficult to differentiate morphologically, and it has been unclear whether they are distinct or a single species with worldwide distribution. Here we report the results of a population genetic study, using nuclear and mitochondrial sequence data from four geographical localities. Our analyses revealed a general lack of geographic structure among Alatina populations, and slight though significant isolation by distance. These data corroborate morphological and behavioral similarities observed in the geographically disparate localities, and indicate the presence of a single, pantropically distributed species, Alatina alata. While repeated, human-mediated introductions of A. alata could explain the patterns we have observed, it seems more likely that genetic metapopulation cohesion is maintained via dispersal through the swimming medusa stage, and perhaps via dispersal of encysted planulae, which are described here for the first time in Alatina.

Transoceanic dispersal and cryptic diversity in a cosmopolitan rafting nudibranch

The aeolid nudibranch Fiona pinnata (Eschscholtz, 1831) is an obligate rafter that occurs exclusively on macroalgal rafts and other floating substrata, and has a seemingly cosmopolitan marine distribution. Mitochondrial (mtDNA) and nuclear DNA sequence data were generated from specimens collected worldwide to test for global connectivity in this species. Phylogeographic analyses revealed three divergent mtDNA lineages, two of which were abundant and widespread. One of these lineages has an apparent circumequatorial distribution, whereas the other may have an antitropical distribution within the Pacific Ocean. Low genetic divergences within each lineage suggest that rafting can mediate dispersal across transoceanic scales. A third, highly divergent, lineage was detected in a single Indonesian specimen. Broadly concordant phylogeographic relationships were detected for the nuclear ITS1 region, with distinct tropical versus antitropical lineages observed. The substantial genetic divergences and largely allopatric distributions observed among the F. pinnata lineages suggest that they represent a species complex.

Temporal Stability of Genetic Structure in a Mesopelagic Copepod

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 (F_CT = 0.285, P < 0.0001 for mtCOII, F_CT = 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.

Multilocus evidence for globally distributed cryptic species and distinct populations across ocean gyres in a mesopelagic copepod

Zooplanktonic taxa have a greater number of distinct populations and species than might be predicted based on their large population sizes and open-ocean habitat, which lacks obvious physical barriers to dispersal and gene flow. To gain insight into the evolutionary mechanisms driving genetic diversification in zooplankton, we developed eight microsatellite markers to examine the population structure of an abundant, globally distributed mesopelagic copepod, Haloptilus longicornis, at 18 sample sites across the Atlantic and Pacific Oceans (n = 761). When comparing our microsatellite results with those of a prior study that used a mtDNA marker (mtCOII, n = 1059, 43 sample sites), we unexpectedly found evidence for the presence of a cryptic species pair. These species were globally distributed and apparently sympatric, and were separated by relatively weak genetic divergence (reciprocally monophyletic mtCOII lineages 1.6% divergent; microsatellite FST ranging from 0.28 to 0.88 across loci, P < 0.00001). Using both mtDNA and microsatellite data for the most common of the two species (n = 669 for microsatellites, n = 572 for mtDNA), we also found evidence for allopatric barriers to gene flow within species, with distinct populations separated by continental landmasses and equatorial waters in both the Atlantic and Pacific Ocean basins. Our study shows that oceanic barriers to gene flow can act as a mechanism promoting allopatric diversification in holoplanktonic taxa, despite the high potential dispersal abilities and pelagic habitat for these species.

Global biogeography and evolution of Cuvierina pteropods

BACKGROUND

Shelled pteropods are planktonic gastropods that are potentially good indicators of the effects of ocean acidification. They also have high potential for the study of zooplankton evolution because they are metazoan plankton with a good fossil record. We investigated phenotypic and genetic variation in pteropods belonging to the genus Cuvierina in relation to their biogeographic distribution across the world's oceans. We aimed to assess species boundaries and to reconstruct their evolutionary history.

RESULTS

We distinguished six morphotypes based on geometric morphometric analyses of shells from 926 museum and 113 fresh specimens. These morphotypes have distinct geographic distributions across the Atlantic, Pacific and Indian oceans, and belong to three major genetic clades based on COI and 28S DNA sequence data. Using a fossil-calibrated phylogeny, we estimated that these clades separated in the Late Oligocene and Early to Middle Miocene. We found evidence for ecological differentiation among all morphotypes based on ecological niche modelling with sea surface temperature, salinity and phytoplankton biomass as primary determinants. Across all analyses, we found highly congruent patterns of differentiation suggesting species level divergences between morphotypes. However, we also found distinct morphotypes (e.g. in the Atlantic Ocean) that were ecologically, but not genetically differentiated.

CONCLUSIONS

Given the distinct ecological and phenotypic specializations found among both described and undescribed Cuvierina taxa, they may not respond equally to future ocean changes and may not be equally sensitive to ocean acidification. Our findings support the view that ecological differentiation may be an important driving force in the speciation of zooplankton.

Phylogenetic information content of Copepoda ribosomal DNA repeat units: ITS1 and ITS2 impact

The utility of various regions of the ribosomal repeat unit for phylogenetic analysis was examined in 16 species representing four families, nine genera, and two orders of the subclass Copepoda (Crustacea). Fragments approximately 2000 bp in length containing the ribosomal DNA (rDNA) 18S and 28S gene fragments, the 5.8S gene, and the internal transcribed spacer regions I and II (ITS1 and ITS2) were amplified and analyzed. The DAMBE (Data Analysis in Molecular Biology and Evolution) software was used to analyze the saturation of nucleotide substitutions; this test revealed the suitability of both the 28S gene fragment and the ITS1/ITS2 rDNA regions for the reconstruction of phylogenetic trees. Distance (minimum evolution) and probabilistic (maximum likelihood, Bayesian) analyses of the data revealed that the 28S rDNA and the ITS1 and ITS2 regions are informative markers for inferring phylogenetic relationships among families of copepods and within the Cyclopidae family and associated genera. Split-graph analysis of concatenated ITS1/ITS2 rDNA regions of cyclopoid copepods suggested that the Mesocyclops, Thermocyclops, and Macrocyclops genera share complex evolutionary relationships. This study revealed that the ITS1 and ITS2 regions potentially represent different phylogenetic signals.

DNA barcoding of marine copepods: assessment of analytical approaches to species identification

More than 2,500 species of copepods (Class Maxillopoda; Subclass Copepoda) occur in the marine planktonic environment. The exceptional morphological conservation of the group, with numerous sibling species groups, makes the identification of species challenging, even for expert taxonomists. Molecular approaches to species identification have allowed rapid detection, discrimination, and identification of species based on DNA sequencing of single specimens and environmental samples. Despite the recent development of diverse genetic and genomic markers, the barcode region of the mitochondrial cytochrome c oxidase subunit I (COI) gene remains a useful and - in some cases - unequaled diagnostic character for species-level identification of copepods. This study reports 800 new barcode sequences for 63 copepod species not included in any previous study and examines the reliability and resolution of diverse statistical approaches to species identification based upon a dataset of 1,381 barcode sequences for 195 copepod species. We explore the impact of missing data (i.e., species not represented in the barcode database) on the accuracy and reliability of species identifications. Among the tested approaches, the best close match analysis resulted in accurate identification of all individuals to species, with no errors (false positives), and out-performed automated tree-based or BLAST based analyses. This comparative analysis yields new understanding of the strengths and weaknesses of DNA barcoding and confirms the value of DNA barcodes for species identification of copepods, including both individual specimens and bulk samples. Continued integrative morphological-molecular taxonomic analysis is needed to produce a taxonomically-comprehensive database of barcode sequences for all species of marine copepods.

Phylogeography of the tropical planktonic foraminifera lineage globigerinella reveals isolation inconsistent with passive dispersal by ocean currents

Morphologically defined species of marine plankton often harbor a considerable level of cryptic diversity. Since many morphospecies show cosmopolitan distribution, an understanding of biogeographic and evolutionary processes at the level of genetic diversity requires global sampling. We use a database of 387 single-specimen sequences of the SSU rDNA of the planktonic foraminifera Globigerinella as a model to assess the biogeographic and phylogenetic distributions of cryptic diversity in marine microplankton on a global scale. Our data confirm the existence of multiple, well isolated genetic lineages. An analysis of their abundance and distribution indicates that our sampling is likely to approximate the actual total diversity. Unexpectedly, we observe an uneven allocation of cryptic diversity among the phylogenetic lineages. We show that this pattern is neither an artifact of sampling intensity nor a function of lineage age. Instead, we argue that it reflects an ongoing speciation process in one of the three major lineages. Surprisingly, four of the six genetic types in the hyperdiverse lineage are biogeographically restricted to the Indopacific. Their mutual co-occurrence and their hierarchical phylogenetic structure provide no evidence for an origin through sudden habitat fragmentation and their limitation to the Indopacific challenges the view of a global gene flow within the warm-water provinces. This phenomenon shows that passive dispersal is not sufficient to describe the distribution of plankton diversity. Rather, these organisms show differentiated distribution patterns shaped by species interactions and reflecting phylogenetic contingency with unique histories of diversification rates.

High evolutionary potential of marine zooplankton

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.