Chaetognaths (arrow worms) are an enigmatic group of marine animals whose phylogenetic position remains elusive, in part because they display a mix of developmental and morphological characters ...associated with other groups 1, 2. In particular, it remains unclear whether they are a sister group to protostomes 1, 2, one of the principal animal superclades, or whether they bear a closer relationship with some spiralian phyla 3, 4. Addressing the phylogenetic position of chaetognaths and refining our understanding of relationships among spiralians are essential to fully comprehend character changes during bilaterian evolution 5. To tackle these questions, we generated new transcriptomes for ten chaetognath species, compiling an extensive phylogenomic dataset that maximizes data occupancy and taxonomic representation. We employed inference methods that consider rate and compositional heterogeneity across taxa to avoid limitations of earlier analyses 6. In this way, we greatly improved the resolution of the protostome tree of life. We find that chaetognaths cluster together with rotifers, gnathostomulids, and micrognathozoans within an expanded Gnathifera clade and that this clade is the sister group to other spiralians 7, 8. Our analysis shows that several previously proposed groupings are likely due to systematic error, and we propose a revised organization of Lophotrochozoa with three main clades: Tetraneuralia (mollusks and entoprocts), Lophophorata (brachiopods, phoronids, and ectoprocts), and a third unnamed clade gathering annelids, nemerteans, and platyhelminthes. Consideration of classical morphological, developmental, and genomic characters in light of this topology indicates secondary loss as a fundamental trend in spiralian evolution.
•Chaetognaths are gnathiferans, the sister group to other spiralians•Lophotrochozoans are divided in three main clades, including Tetraneuralia•A new clade is made of annelids, nemerteans, and platyhelminths•Compositional heterogeneity affects the reconstruction of the protostome tree
Marlétaz et al. used new data and refined methods to improve the phylogenetic placement of chaetognaths, an enigmatic group of marine animals. They relate to gnathiferans, a sister group of other spiralians, and include rotifers and gnathostomulids. Spiralians are subdivided in three main clades (including Tetraneuralia and Lophotrochozoa).
Shelled pteropods are widely regarded as bioindicators for ocean acidification, because their fragile aragonite shells are susceptible to increasing ocean acidity. While short-term incubations have ...demonstrated that pteropod calcification is negatively impacted by ocean acidification, we know little about net calcification in response to varying ocean conditions in natural populations. Here, we examine in situ calcification of Limacina helicina pteropods collected from the California Current Ecosystem, a coastal upwelling system with strong spatial gradients in ocean carbonate chemistry, dissolved oxygen and temperature. Depth-averaged pH ranged from 8.03 in warmer offshore waters to 7.77 in cold CO
-rich waters nearshore. Based on high-resolution micro-CT technology, we showed that shell thickness declined by ~ 37% along the upwelling gradient from offshore to nearshore water. Dissolution marks covered only ~ 2% of the shell surface area and were not associated with the observed variation in shell thickness. We thus infer that pteropods make thinner shells where upwelling brings more acidified and colder waters to the surface. Probably the thinner shells do not result from enhanced dissolution, but are due to a decline in calcification. Reduced calcification of pteropods is likely to have major ecological and biogeochemical implications for the cycling of calcium carbonate in the oceans.
Characterization of species diversity of zooplankton is key to understanding, assessing, and predicting the function and future of pelagic ecosystems throughout the global ocean. The marine ...zooplankton assemblage, including only metazoans, is highly diverse and taxonomically complex, with an estimated ~28,000 species of 41 major taxonomic groups. This review provides a comprehensive summary of DNA sequences for the barcode region of mitochondrial cytochrome oxidase I (COI) for identified specimens. The foundation of this summary is the MetaZooGene Barcode Atlas and Database (MZGdb), a new open-access data and metadata portal that is linked to NCBI GenBank and BOLD data repositories. The MZGdb provides enhanced quality control and tools for assembling COI reference sequence databases that are specific to selected taxonomic groups and/or ocean regions, with associated metadata (e.g., collection georeferencing, verification of species identification, molecular protocols), and tools for statistical analysis, mapping, and visualization. To date, over 150,000 COI sequences for ~ 5600 described species of marine metazoan plankton (including holo- and meroplankton) are available via the MZGdb portal. This review uses the MZGdb as a resource for summaries of COI barcode data and metadata for important taxonomic groups of marine zooplankton and selected regions, including the North Atlantic, Arctic, North Pacific, and Southern Oceans. The MZGdb is designed to provide a foundation for analysis of species diversity of marine zooplankton based on DNA barcoding and metabarcoding for assessment of marine ecosystems and rapid detection of the impacts of climate change.
Pteropods, a group of holoplanktonic gastropods, are regarded as bioindicators of the effects of ocean acidification on open ocean ecosystems, because their thin aragonitic shells are susceptible to ...dissolution. While there have been recent efforts to address their capacity for physiological acclimation, it is also important to gain predictive understanding of their ability to adapt to future ocean conditions. However, little is known about the levels of genetic variation and large‐scale population structuring of pteropods, key characteristics enabling local adaptation. We examined the spatial distribution of genetic diversity in the mitochondrial cytochrome c oxidase I (COI) and nuclear 28S gene fragments, as well as shell shape variation, across a latitudinal transect in the Atlantic Ocean (35°N–36°S) for the pteropod Limacina bulimoides. We observed high levels of genetic variability (COI π = 0.034, 28S π = 0.0021) and strong spatial structuring (COI ΦST = 0.230, 28S ΦST = 0.255) across this transect. Based on the congruence of mitochondrial and nuclear differentiation, as well as differences in shell shape, we identified a primary dispersal barrier in the southern Atlantic subtropical gyre (15–18°S). This barrier is maintained despite the presence of expatriates, a gyral current system, and in the absence of any distinct oceanographic gradients in this region, suggesting that reproductive isolation between these populations must be strong. A secondary dispersal barrier supported only by 28S pairwise ΦST comparisons was identified in the equatorial upwelling region (between 15°N and 4°S), which is concordant with barriers observed in other zooplankton species. Both oceanic dispersal barriers were congruent with regions of low abundance reported for a similar basin‐scale transect that was sampled 2 years later. Our finding supports the hypothesis that low abundance indicates areas of suboptimal habitat that result in barriers to gene flow in widely distributed zooplankton species. Such species may in fact consist of several populations or (sub)species that are adapted to local environmental conditions, limiting their potential for adaptive responses to ocean changes. Future analyses of genome‐wide diversity in pteropods could provide further insight into the strength, formation and maintenance of oceanic dispersal barriers.
Shelled pteropods or 'sea butterflies' are proposed as bio‐indicators to monitor impacts of ocean acidification on open ocean ecosystems. Local adaptation is affected by the presence and geographic extent of distinct populations, however, little is known about the large‐scale population structure of pteropods. We thus examined the spatial distribution of genetic and morphometric variation in the species Limacina bulimoides across the Atlantic Ocean. Two dispersal barriers were identified that were congruent with regions of low abundance, suggesting that regions of suboptimal habitat may create barriers to gene flow in widely‐distributed zooplankton species.
Pteropods are planktonic gastropods that are considered as bio-indicators to monitor impacts of ocean acidification on marine ecosystems. In order to gain insight into their adaptive potential to ...future environmental changes, it is critical to use adequate molecular tools to delimit species and population boundaries and to assess their genetic connectivity. We developed a set of target capture probes to investigate genetic variation across their large-sized genome using a population genomics approach. Target capture is less limited by DNA amount and quality than other genome-reduced representation protocols, and has the potential for application on closely related species based on probes designed from one species.
We generated the first draft genome of a pteropod, Limacina bulimoides, resulting in a fragmented assembly of 2.9 Gbp. Using this assembly and a transcriptome as a reference, we designed a set of 2899 genome-wide target capture probes for L. bulimoides. The set of probes includes 2812 single copy nuclear targets, the 28S rDNA sequence, ten mitochondrial genes, 35 candidate biomineralisation genes, and 41 non-coding regions. The capture reaction performed with these probes was highly efficient with 97% of the targets recovered on the focal species. A total of 137,938 single nucleotide polymorphism markers were obtained from the captured sequences across a test panel of nine individuals. The probes set was also tested on four related species: L. trochiformis, L. lesueurii, L. helicina, and Heliconoides inflatus, showing an exponential decrease in capture efficiency with increased genetic distance from the focal species. Sixty-two targets were sufficiently conserved to be recovered consistently across all five species.
The target capture protocol used in this study was effective in capturing genome-wide variation in the focal species L. bulimoides, suitable for population genomic analyses, while providing insights into conserved genomic regions in related species. The present study provides new genomic resources for pteropods and supports the use of target capture-based protocols to efficiently characterise genomic variation in small non-model organisms with large genomes.
Pteropods are a group of planktonic gastropods that are widely regarded as biological indicators for assessing the impacts of ocean acidification. Their aragonitic shells are highly sensitive to ...acute changes in ocean chemistry. However, to gain insight into their potential to adapt to current climate change, we need to accurately reconstruct their evolutionary history and assess their responses to past changes in the Earth’s carbon cycle. Here, we resolve the phylogeny and timing of pteropod evolution with a phylogenomic dataset (2,654 genes) incorporating new data for 21 pteropod species and revised fossil evidence. In agreement with traditional taxonomy, we recovered molecular support for a division between “sea butterflies” (Thecosomata; mucus-web feeders) and “sea angels” (Gymnosomata; active predators). Molecular dating demonstrated that these two lineages diverged in the early Cretaceous, and that all main pteropod clades, including shelled, partially-shelled, and unshelled groups, diverged in the mid- to late Cretaceous. Hence, these clades originated prior to and subsequently survived major global change events, including the Paleocene–Eocene Thermal Maximum (PETM), the closest analog to modern-day ocean acidification and warming. Our findings indicate that planktonic aragonitic calcifiers have shown resilience to perturbations in the Earth’s carbon cycle over evolutionary timescales.
Pteropods are a widespread group of holoplanktonic gastropod molluscs and are uniquely suitable for study of long-term evolutionary processes in the open ocean because they are the only living ...metazoan plankton with a good fossil record. Pteropods have been proposed as bioindicators to monitor the impacts of ocean acidification and in consequence have attracted considerable research interest, however, a robust evolutionary framework for the group is still lacking. Here we reconstruct their phylogenetic relationships and examine the evolutionary history of pteropods based on combined analyses of Cytochrome Oxidase I, 28S, and 18S ribosomal rRNA sequences and a molecular clock calibrated using fossils and the estimated timing of the formation of the Isthmus of Panama. Euthecosomes with uncoiled shells were monophyletic with Creseis as the earliest diverging lineage, estimated at 41-38 million years ago (mya). The coiled euthecosomes (Limacina, Heliconoides, Thielea) were not monophyletic contrary to the accepted morphology-based taxonomy; however, due to their high rate heterogeneity no firm conclusions can be drawn. We found strong support for monophyly of most euthecosome genera, but Clio appeared as a polyphyletic group, and Diacavolinia grouped within Cavolinia, making the latter genus paraphyletic. The highest evolutionary rates were observed in Heliconoides inflatus and Limacina bulimoides for both 28S and 18S partitions. Using a fossil-calibrated phylogeny that sets the first occurrence of coiled euthecosomes at 79-66 mya, we estimate that uncoiled euthecosomes evolved 51-42 mya and that most extant uncoiled genera originated 40-15 mya. These findings are congruent with a molecular clock analysis using the Isthmus of Panama formation as an independent calibration. Although not all phylogenetic relationships could be resolved based on three molecular markers, this study provides a useful resource to study pteropod diversity and provides general insight into the processes that generate and maintain their diversity in the open ocean.
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.
The aragonite shelled, planktonic gastropod family Atlantidae (shelled heteropods) is likely to be one of the first groups to be impacted by imminent ocean changes, including ocean warming and ocean ...acidification. With a fossil record spanning at least 100 Ma, atlantids have experienced and survived global-scale ocean changes and extinction events in the past. However, the diversification patterns and tempo of evolution in this family are largely unknown.
Based on a concatenated maximum likelihood phylogeny of three genes (cytochrome c oxidase subunit 1 mitochondrial DNA, 28S and 18S ribosomal rRNA) we show that the three extant genera of the family Atlantidae, Atlanta, Protatlanta and Oxygyrus, form monophyletic groups. The genus Atlanta is split into two groups, one exhibiting smaller, well ornamented shells, and the other having larger, less ornamented shells. The fossil record, in combination with a fossil-calibrated phylogeny, suggests that large scale atlantid extinction was accompanied by considerable and rapid diversification over the last 25 Ma, potentially driven by vicariance events.
Now confronted with a rapidly changing modern ocean, the ability of atlantids to survive past global change crises gives some optimism that they may be able to persist through the Anthropocene.
The existence of multiple independently derived populations in landlocked marine lakes provides an opportunity for fundamental research into the role of isolation in population divergence and ...speciation in marine taxa. Marine lakes are landlocked water bodies that maintain a marine character through narrow submarine connections to the sea and could be regarded as the marine equivalents of terrestrial islands. The sponge Suberites diversicolor (Porifera: Demospongiae: Suberitidae) is typical of marine lake habitats in the Indo-Australian Archipelago. Four molecular markers (two mitochondrial and two nuclear) were employed to study genetic structure of populations within and between marine lakes in Indonesia and three coastal locations in Indonesia, Singapore and Australia. Within populations of S. diversicolor two strongly divergent lineages (A & B) (COI: p = 0.4% and ITS: p = 7.3%) were found, that may constitute cryptic species. Lineage A only occurred in Kakaban lake (East Kalimantan), while lineage B was present in all sampled populations. Within lineage B, we found low levels of genetic diversity in lakes, though there was spatial genetic population structuring. The Australian population is genetically differentiated from the Indonesian populations. Within Indonesia we did not record an East-West barrier, which has frequently been reported for other marine invertebrates. Kakaban lake is the largest and most isolated marine lake in Indonesia and contains the highest genetic diversity with genetic variants not observed elsewhere. Kakaban lake may be an area where multiple putative refugia populations have come into secondary contact, resulting in high levels of genetic diversity and a high number of endemic species.