MetaZooGene Atlas and Database: Reference Sequences for Marine Ecosystems.

The MetaZooGene Atlas and Database (MZGdb; https://metazoogene.org/mzgdb/ ) is an open-access data and metadata portal synchronized with the NCBI GenBank and BOLD data repositories. The MZGdb includes sequences for genes used for the classification and identification of marine organisms based on DNA barcoding and metabarcoding. The focus of the MZGdb is biodiversity of marine ecosystems, including phytoplankton and microbes, zooplankton and invertebrates, fish, and other marine vertebrates (pinnipeds, cetaceans, and sea turtles). DNA sequences currently included are mitochondrial cytochrome oxidase I (COI), 12S, and 16S rRNA, and nuclear 18S and 28S rRNA. The MZGdb provides data and mapping tools for assembling and downloading compilations of reference sequence data that are specific to selected genes, taxonomic groups, and/or ocean regions. An additional feature of the MZGdb is the Atlas which summarizes data coverage and proportional completeness based on statistics of species with available sequences versus species commonly found in each ocean region.This chapter is a collaborative effort of the Scientific Committee for Ocean Research (SCOR) Working Group WG157: MetaZooGene: Toward a new global view of marine zooplankton biodiversity based on DNA metabarcoding and reference DNA sequence databases ( https://metazoogene.org ).

Transition from stromatolite to thrombolite fabric: potential role for reticulopodial protists in lake microbialites of a Proterozoic ecosystem analog.

Prior observations suggest that foraminiferan protists use their reticulopodia (anastomosing pseudopodia) to alter sediment fabric by disrupting laminations of subtidal marine stromatolites, erasing the layered structures in an experimental setting. Because microbialites and foraminifera are found in non-marine settings, we hypothesized that foraminifera living in lakes could also disrupt layered microbialite fabric. With this aim and using a variety of multidisciplinary approaches, we conducted field surveys and an experiment on microbialites from Green Lake (GL; Fayetteville, New York State, United States), which has been studied as a Proterozoic ecosystem analog. The lake is meromictic and alkaline, receiving calcium sulfate-rich water in the monimolimnion; it supports a well-developed carbonate platform that provides access to living and relict microbialites. The living microbialites grow from early spring to autumn, forming a laminated mat at their surface (top ~5 mm), but a clotted or massive structure exists at depth (> ~ 1 cm). We observed a morphotype of "naked" foraminiferan-like protist in samples from GL microbialites and sediments; thus, considered the possibility of freshwater foraminiferan impact on microbialite fabric. Results of an experiment that seeded the cultured freshwater foraminifer Haplomyxa saranae onto the GL microbialite surface indicates via micro-CT scanning and anisotropy analysis that the introduced foraminifer impacted uppermost microbialite layering (n = 3 cores); those cores with an added inhibitor lacked changes in anisotropy for two of those three cores. Thus, it remains plausible that the much smaller, relatively common, native free-form reticulate protist, which we identified as Chlamydomyxa labyrinthuloides, can disrupt microbialite fabrics on sub-millimeter scales. Our observations do not exclude contributions of other possible causal factors.

Salpa genome and developmental transcriptome analyses reveal molecular flexibility enabling reproductive success in a rapidly changing environment.

Ocean warming favors pelagic tunicates, such as salps, that exhibit increasingly frequent and rapid population blooms, impacting trophic dynamics and composition and human marine-dependent activities. Salp blooms are a result of their successful reproductive life history, alternating seasonally between asexual and sexual protogynous (i.e. sequential) hermaphroditic stages. While predicting future salp bloom frequency and intensity relies on an understanding of the transitions during the sexual stage from female through parturition and subsequent sex change to male, these transitions have not been explored at the molecular level. Here we report the development of the first complete genome of S. thompsoni and the North Atlantic sister species S. aspera. Genome and comparative analyses reveal an abundance of repeats and G-quadruplex (G4) motifs, a highly stable secondary structure, distributed throughout both salp genomes, a feature shared with other tunicates that perform alternating sexual-asexual reproductive strategies. Transcriptional analyses across sexual reproductive stages for S. thompsoni revealed genes associated with male sex differentiation and spermatogenesis are expressed as early as birth and before parturition, inconsistent with previous descriptions of sequential sexual differentiation in salps. Our findings suggest salp are poised for reproductive success at birth, increasing the potential for bloom formation as ocean temperatures rise.

Genetics redraws pelagic biogeography of Calanus.

Planktonic copepods of the genus Calanus play a central role in North Atlantic/Arctic marine food webs. Here, using molecular markers, we redrew the distributional ranges of Calanus species inhabiting the North Atlantic and Arctic Oceans and revealed much wider and more broadly overlapping distributions than previously described. The Arctic shelf species, C. glacialis, dominated the zooplankton assemblage of many Norwegian fjords, where only C. finmarchicus has been reported previously. In these fjords, high occurrences of the Arctic species C. hyperboreus were also found. Molecular markers revealed that the most common method of species identification, prosome length, cannot reliably discriminate the species in Norwegian fjords. Differences in degree of genetic differentiation among fjord populations of the two species suggested that C. glacialis is a more permanent resident of the fjords than C. finmarchicus We found no evidence of hybridization between the species. Our results indicate a critical need for the wider use of molecular markers to reliably identify and discriminate these morphologically similar copepod species, which serve as important indicators of climate responses.

Rapid Evolutionary Rates and Unique Genomic Signatures Discovered in the First Reference Genome for the Southern Ocean Salp, Salpa thompsoni (Urochordata, Thaliacea).

A preliminary genome sequence has been assembled for the Southern Ocean salp, Salpa thompsoni (Urochordata, Thaliacea). Despite the ecological importance of this species in Antarctic pelagic food webs and its potential role as an indicator of changing Southern Ocean ecosystems in response to climate change, no genomic resources are available for S. thompsoni or any closely related urochordate species. Using a multiple-platform, multiple-individual approach, we have produced a 318,767,936-bp genome sequence, covering >50% of the estimated 602 Mb (±173 Mb) genome size for S. thompsoni Using a nonredundant set of predicted proteins, >50% (16,823) of sequences showed significant homology to known proteins and ∼38% (12,151) of the total protein predictions were associated with Gene Ontology functional information. We have generated 109,958 SNP variant and 9,782 indel predictions for this species, serving as a resource for future phylogenomic and population genetic studies. Comparing the salp genome to available assemblies for four other urochordates, Botryllus schlosseri, Ciona intestinalis, Ciona savignyi and Oikopleura dioica, we found that S. thompsoni shares the previously estimated rapid rates of evolution for these species. High mutation rates are thus independent of genome size, suggesting that rates of evolution >1.5 times that observed for vertebrates are a broad taxonomic characteristic of urochordates. Tests for positive selection implemented in PAML revealed a small number of genes with sites undergoing rapid evolution, including genes involved in ribosome biogenesis and metabolic and immune process that may be reflective of both adaptation to polar, planktonic environments as well as the complex life history of the salps. Finally, we performed an initial survey of small RNAs, revealing the presence of known, conserved miRNAs, as well as novel miRNA genes; unique piRNAs; and mature miRNA signatures for varying developmental stages. Collectively, these resources provide a genomic foundation supporting S. thompsoni as a model species for further examination of the exceptional rates and patterns of genomic evolution shown by urochordates. Additionally, genomic data will allow for the development of molecular indicators of key life history events and processes and afford new understandings and predictions of impacts of climate change on this key species of Antarctic pelagic ecosystems.

Phylogeography and connectivity of the Pseudocalanus (Copepoda: Calanoida) species complex in the eastern North Pacific and the Pacific Arctic Region.

The genus Pseudocalanus (Copepoda, Calanoida) is among the most numerically dominant copepods in eastern North Pacific and Pacific-Arctic waters. We compared population connectivity and phylogeography based on DNA sequence variation for a portion of the mitochondrial cytochrome oxidase I gene for four Pseudocalanus species with differing biogeographical ranges within these ocean regions. Genetic analyses were linked to characterization of biological and physical environmental variables for each sampled region. Haplotype diversity was higher for the temperate species (Pseudocalanus mimus and Pseudocalanus newmani) than for the Arctic species (Pseudocalanus acuspes and Pseudocalanus minutus). Genetic differentiation among populations at regional scales was observed for all species, except P. minutus. The program Migrate-N tested the likelihood of alternative models of directional gene flow between sampled populations in relation to oceanographic features. Model results estimated predominantly northward gene flow from the Gulf of Alaska to the Beaufort Sea for P. newmani. Model scenarios that allowed bidirectional gene flow between sampled populations gave the best Bayesian predictions for P. acuspes, P. mimus and P. minutus. Under current warming trends, biogeographical boundaries and barriers for Pseudocalanus species may shift, allowing habitat range expansion or contraction and resulting in altered population connectivity between Arctic and sub-Arctic populations.

Identification, Discrimination, and Discovery of Species of Marine Planktonic Ostracods Using DNA Barcodes.

The Ostracoda (Crustacea; Class Ostracoda) is a diverse, frequently abundant, and ecologically important component of the marine zooplankton assemblage. There are more than 200 described species of marine planktonic ostracods, many of which (especially conspecific species) can be identified only by microscopic examination and dissection of fragile morphological characters. Given the complexity of species identification and increasing lack of expert taxonomists, DNA barcodes (short DNA sequences for species discrimination and identification) are particularly useful and necessary. Results are reported from analysis of 210 specimens of 78 species of marine planktonic ostracods, including two novel species, and 51 species for which barcodes have not been previously published. Specimens were collected during 2006 to 2008 from the Atlantic, Indian, and Southern Oceans, Greenland Sea and Gulf of Alaska. Samples were collected from surface to 5,000 m using various collection devices. DNA sequence variation was analyzed for a 598 base-pair region of the mitochondrial cytochrome oxidase subunit I (COI) gene. Kimura-2-Parameter (K2P) genetic distances within described species (mean = 0.010 ± 0.017 SD) were significantly smaller than between species (0.260 + 0.080), excluding eight taxa hypothesized to comprise cryptic species due to morphological variation (especially different size forms) and/or collection from different geographic regions. These taxa showed similar K2P distance values within (0.014 + 0.026) and between (0.221 ± 0.068) species. All K2P distances > 0.1 resulted from comparisons between identified or cryptic species, with no overlap between intra- and interspecific genetic distances. A Neighbor Joining tree resolved nearly all described species analyzed, with multiple sequences forming monophyletic clusters with high bootstrap values (typically 99%). Based on taxonomically and geographically extensive sampling and analysis (albeit with small sample sizes), the COI barcode region was shown to be a valuable character for discrimination, recognition, identification, and discovery of species of marine planktonic ostracods.

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.

In silico characterization of the insect diapause-associated protein couch potato (CPO) in Calanus finmarchicus (Crustacea: Copepoda).

Couch potato (CPO) is an RNA-binding protein involved in the regulation of nervous system development and adult diapause in insects. Within insects, this protein is highly conserved, yet it has not been identified in another large arthropod group, the Crustacea. Here, functional genomics was used to identify putative CPO homologs in the copepod Calanus finmarchicus, a planktonic crustacean that undergoes seasonal diapause. In silico mining of expressed sequence tag (EST) and 454 pyrosequencing data resulted in the identification of two full-length CPO proteins, each 205 amino acids long. The two C. finmarchicus CPOs (Calfi-CPO I and II) are identical in sequence with the exception of three amino acids, and are predicted to possess a single RNA recognition motif (RRM). Sequence comparison of the two Calfi-CPOs with those of insects shows high levels of amino acid conservation, particularly in their RRMs. Using the C. finmarchicus sequences as queries, ESTs encoding partial CPOs were identified from two other crustaceans, the parasitic copepod Lernaeocera branchialis and shrimp Penaeus monodon. Surprisingly, no convincing CPO-encoding transcripts were identified from crustacean species with very large (>100,000) EST datasets (e.g. Litopenaeus vannamei, Daphnia pulex and Lepeophtheirus salmonis), suggesting that CPO transcript/protein may be expressed at very low levels or absent in some crustaceans. RNA-Seq data suggested stage-specific expression of CPO in C. finmarchicus, with few transcripts present in eggs (which represent mixed embryonic stages) and adults, and high levels in nauplii and copepodites; stages exhibiting high CPO expression are consistent with a role for it in neuronal development.

Multi-gene analysis reveals a lack of genetic divergence between Calanus agulhensis and C. sinicus (Copepoda; Calanoida).

The discrimination and taxonomic identification of marine species continues to pose a challenge despite the growing number of diagnostic metrics and approaches. This study examined the genetic relationship between two sibling species of the genus Calanus (Crustacea; Copepoda; Calanidae), C. agulhensis and C. sinicus, using a multi-gene analysis. DNA sequences were determined for portions of the mitochondrial cytochrome c oxidase I (mtCOI); nuclear citrate synthase (CS), and large subunit (28S) rRNA genes for specimens collected from the Sea of Japan and North East (NE) Pacific Ocean for C. sinicus and from the Benguela Current and Agulhas Bank, off South Africa, for C. agulhensis. For mtCOI, C. sinicus and C. agulhensis showed similar levels of haplotype diversity (H(d) = 0.695 and 0.660, respectively) and nucleotide diversity (π = 0.003 and 0.002, respectively). Pairwise F(ST) distances for mtCOI were significant only between C. agulhensis collected from the Agulhas and two C. sinicus populations: the Sea of Japan (F(ST) = 0.152, p<0.01) and NE Pacific (F(ST) = 0.228, p<0.005). Between the species, F(ST) distances were low for both mtCOI (F(ST) = 0.083, p = 0.003) and CS (F(ST) = 0.050, p = 0.021). Large subunit (28S) rRNA showed no variation between the species. Our results provide evidence of the lack of genetic distinction of C. sinicus and C. agulhensis, raise questions of whether C. agulhensis warrants status as a distinct species, and indicate the clear need for more intensive and extensive ecological and genetic analysis.

Molecular systematic of three species of Oithona (Copepoda, Cyclopoida) from the Atlantic Ocean: comparative analysis using 28S rDNA.

Species of Oithona (Copepoda, Cyclopoida) are highly abundant, ecologically important, and widely distributed throughout the world oceans. Although there are valid and detailed descriptions of the species, routine species identifications remain challenging due to their small size, subtle morphological diagnostic traits, and the description of geographic forms or varieties. This study examined three species of Oithona (O. similis, O. atlantica and O. nana) occurring in the Argentine sector of the South Atlantic Ocean based on DNA sequence variation of a 575 base-pair region of 28S rDNA, with comparative analysis of these species from other North and South Atlantic regions. DNA sequence variation clearly resolved and discriminated the species, and revealed low levels of intraspecific variation among North and South Atlantic populations of each species. The 28S rDNA region was thus shown to provide an accurate and reliable means of identifying the species throughout the sampled domain. Analysis of 28S rDNA variation for additional species collected throughout the global ocean will be useful to accurately characterize biogeographical distributions of the species and to examine phylogenetic relationships among them.

Functional genomics resources for the North Atlantic copepod, Calanus finmarchicus: EST database and physiological microarray.

The copepod, Calanus finmarchicus is a keystone species for the North Atlantic. Because of recent changes in the geographic distribution of this species, there are questions as to how this organism responds physiologically to environmental cues. Molecular techniques allow for examination and new understanding of these physiological changes. Here, we describe the development of a microarray for high-throughput studies of the physiological ecology of C. finmarchicus. An EST database was generated for this species using a normalized cDNA library derived from adult and sub-adult individuals. Sequence data were clustered into contigs and annotated using Blastx. Target transcripts were selected, and unique, 50 base-pair, oligomer probes were generated for 995 genes. Blast2GO processing provided detailed information on gene function. The selected targets included broad representation of biological processes, cellular components, and molecular functions. The microarray was tested in two sets of comparisons: adult females maintained at different food concentrations and field-caught sub-adults showing differences in lipid storage. Up-regulated and down-regulated transcripts were identified for both comparisons. Only a small subset of the genes up-regulated in low food individuals were also up-regulated in lipid-poor animals; no overlap was seen between the genes down-regulated in the two comparisons.

Gone with the currents: lack of genetic differentiation at the circum-continental scale in the Antarctic krill Euphausia superba.

BACKGROUND: Southern Ocean fauna represent a significant amount of global biodiversity, whose origin may be linked to glacial cycles determining local extinction/eradication with ice advance, survival of refugee populations and post-glacial re-colonization. This pattern implies high potential for differentiation in benthic shelf species with limited dispersal, yet consequences for pelagic organisms are less clear. The present study investigates levels of genetic variation and population structure of the Antarctic krill Euphausia superba using mitochondrial DNA and EST-linked microsatellite markers for an unprecedentedly comprehensive sampling of its populations over a circum-Antarctic range. RESULTS: MtDNA (ND1) sequences and EST-linked microsatellite markers indicated no clear sign of genetic structure among populations over large geographic scales, despite considerable power to detect differences inferred from forward-time simulations. Based on ND1, few instances of genetic heterogeneity, not significant after correction for multiple tests, were detected between geographic or temporal samples. Neutrality tests and mismatch distribution based on mtDNA sequences revealed strong evidence of past population expansion. Significant positive values of the parameter g (a measure of population growth) were obtained from microsatellite markers using a coalescent-based genealogical method and suggested a recent start (60,000 - 40,000 years ago) for the expansion. CONCLUSIONS: The results provide evidence of lack of genetic heterogeneity of Antarctic krill at large geographic scales and unequivocal support for recent population expansion. Lack of genetic structuring likely reflects the tight link between krill and circum-Antarctic ocean currents and is consistent with the hypothesis that differentiation processes in Antarctic species are largely influenced by dispersal potential, whereas small-scale spatial and temporal differentiation might be due to local conditions leading to genetic patchiness. The signal of recent population growth suggests differential impact of glacial cycles on pelagic Antarctic species, which experienced population expansion during glaciations with increased available habitat, versus sedentary benthic shelf species. EST-linked microsatellites provide new perspectives to complement the results based on mtDNA and suggest that data-mining of EST libraries will be a useful approach to facilitate use of microsatellites for additional species.

DNA barcoding of marine metazoa.

More than 230,000 known species representing 31 metazoan phyla populate the world's oceans. Perhaps another 1,000,000 or more species remain to be discovered. There is reason for concern that species extinctions may out-pace discovery, especially in diverse and endangered marine habitats such as coral reefs. DNA barcodes (i.e., short DNA sequences for species recognition and discrimination) are useful tools to accelerate species-level analysis of marine biodiversity and to facilitate conservation efforts. This review focuses on the usual barcode region for metazoans: a approximately 648 base-pair region of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Barcodes have also been used for population genetic and phylogeographic analysis, identification of prey in gut contents, detection of invasive species, forensics, and seafood safety. More controversially, barcodes have been used to delimit species boundaries, reveal cryptic species, and discover new species. Emerging frontiers are the use of barcodes for rapid and increasingly automated biodiversity assessment by high-throughput sequencing, including environmental barcoding and the use of barcodes to detect species for which formal identification or scientific naming may never be possible.

Molecular phylogeny of the Calanoida (Crustacea: Copepoda).

The order Calanoida includes some of the most successful planktonic groups in both marine and freshwater environments. Due to the morphological complexity of the taxonomic characters in this group, subdivision and phylogenies have been complex and problematic. This study establishes a multi-gene molecular phylogeny of the calanoid copepods based upon small (18S) and large (28S) subunits of nuclear ribosomal RNA genes and mitochondrial encoded cytochrome b and cytochrome c oxidase subunit-I genes, including 29 families from 7 superfamilies of the order. This analysis is more comprehensive than earlier studies in terms of number of families, range of molecular markers, and breadth of taxonomic levels resolved. Patterns of divergence of ribosomal RNA genes are shown to be significantly heterogeneous among superfamilies, providing a likely explanation for disparate results of previous studies. The multi-gene phylogeny recovers a monophyletic Calanoida, as well as the superfamilies Augaptiloidea, Centropagoidea, Bathypontioidea, Eucalanoidea, Spinocalanoidea and Clausocalanoidea. The phylogeny largely agrees with previously-published morphological phylogenies, including e.g., enlargement of the Bathypontioidea to include the Fosshageniidae.

Barcoding of arrow worms (Phylum Chaetognatha) from three oceans: genetic diversity and evolution within an enigmatic phylum.

Arrow worms (Phylum Chaetognatha) are abundant planktonic organisms and important predators in many food webs; yet, the classification and evolutionary relationships among chaetognath species remain poorly understood. A seemingly simple body plan is underlain by subtle variation in morphological details, obscuring the affinities of species within the phylum. Many species achieve near global distributions, spanning the same latitudinal bands in all ocean basins, while others present disjunct ranges, in some cases with the same species apparently found at both poles. To better understand how these complex evolutionary and geographic variables are reflected in the species makeup of chaetognaths, we analyze DNA barcodes of the mitochondrial cytochrome oxidase c subunit I (COI) gene, from 52 specimens of 14 species of chaetognaths collected mainly from the Atlantic Ocean. Barcoding analysis was highly successful at discriminating described species of chaetognaths across the phylum, and revealed little geographical structure. This barcode analysis reveals hitherto unseen genetic variation among species of arrow worms, and provides insight into some species relationships of this enigmatic group.