Validating a molecular clock for nudibranchs-No fossils to the rescue.

Time calibrated phylogenies are typically reconstructed with fossil information but for soft-bodied marine invertebrates that lack hard parts, a fossil record is lacking. In these cases, biogeographic calibrations or the rates of divergence for related taxa are often used. Although nudibranch phylogenies have advanced with the input of molecular data, no study has derived a divergence rate for this diverse group of invertebrates. Here, we use an updated closure date for the Isthmus of Panama (2.8 Ma) to derive the first divergence rates for chromodorid nudibranchs using multigene data from a geminate pair with broad phylogeographic sampling. Examining the species Chromolaichma sedna (Marcus & Marcus, 1967), we uncover deep divergences among eastern Pacific and western Atlantic clades and we erect a new species designation for the latter (Chromolaichma hemera sp. nov.). Next, we discover extensive phylogeographic structure within C. hemera sp. nov. sensu lato, thereby refuting the hypothesis of a recent introduction. Lastly, we derive divergence rates for mitochondrial and nuclear loci that exceed known rates for other gastropods and we highlight significant rate heterogeneity both among markers and taxa. Together, these findings improve understanding of nudibranch systematics and provide rates useful to apply to divergence scenarios in this diverse group.

Geography, environment, and colonization history interact with morph type to shape genomic variation in an Arctic fish.

Polymorphic species are useful models for investigating the evolutionary processes driving diversification. Such processes include colonization history as well as contemporary selection, gene flow, and genetic drift, which can vary between intraspecific morphs as a function of their distinct life histories. The interactive and relative influence of such evolutionary processes on morph differentiation critically informs morph-specific management decisions and our understanding of incipient speciation. We therefore investigated how geographic distance, environmental conditions, and colonization history interacted with morph migratory capacity in the highly polymorphic fish species, Arctic Charr (Salvelinus alpinus). Using an 87 k SNP chip we genetically characterized recently evolved anadromous, resident, and landlocked charr collected from 45 locations across a secondary contact zone of three charr glacial lineages in eastern Canada. A strong pattern of isolation by distance across all populations suggested geographic distance principally shaped genetic structure. Landlocked populations had lower genetic diversities and higher genetic differentiation than anadromous populations. However, effective population size was generally temporally stable in landlocked populations in comparison to anadromous populations. Genetic diversity positively correlated with latitude, potentially indicating southern anadromous populations' vulnerability to climate change and greater introgression between the Arctic and Atlantic glacial lineages in northern Labrador. Local adaptation was suggested by the observation of several environmental variables strongly associating with functionally relevant outlier genes including a region on chromosome AC21 potentially associated with anadromy. Our results demonstrate that gene flow, colonization history, and local adaptation uniquely interact to influence the genetic variation and evolutionary trajectory of populations.

Application of Omics Tools in Designing and Monitoring Marine Protected Areas For a Sustainable Blue Economy.

A key component of the global blue economy strategy is the sustainable extraction of marine resources and conservation of marine environments through networks of marine protected areas (MPAs). Connectivity and representativity are essential factors that underlie successful implementation of MPA networks, which can safeguard biological diversity and ecosystem function, and ultimately support the blue economy strategy by balancing ocean use with conservation. New "big data" omics approaches, including genomics and transcriptomics, are becoming essential tools for the development and maintenance of MPA networks. Current molecular omics techniques, including population-scale genome sequencing, have direct applications for assessing population connectivity and for evaluating how genetic variation is represented within and among MPAs. Effective baseline characterization and long-term, scalable, and comprehensive monitoring are essential for successful MPA management, and omics approaches hold great promise to characterize the full range of marine life, spanning the microbiome to megafauna across a range of environmental conditions (shallow sea to the deep ocean). Omics tools, such as eDNA metabarcoding can provide a cost-effective basis for biodiversity monitoring in large and remote conservation areas. Here we provide an overview of current omics applications for conservation planning and monitoring, with a focus on metabarcoding, metagenomics, and population genomics. Emerging approaches, including whole-genome sequencing, characterization of genomic architecture, epigenomics, and genomic vulnerability to climate change are also reviewed. We demonstrate that the operationalization of omics tools can enhance the design, monitoring, and management of MPAs and thus will play an important role in a modern and comprehensive blue economy strategy.

The Genomic Consistency of the Loss of Anadromy in an Arctic Fish (Salvelinus alpinus).

AbstractThe potentially significant genetic consequences associated with the loss of migratory capacity of diadromous fishes that have become landlocked in freshwater are poorly understood. Consistent selective pressures associated with freshwater residency may drive repeated differentiation both between allopatric landlocked and anadromous populations and within landlocked populations (resulting in sympatric morphs). Alternatively, the strong genetic drift anticipated in isolated landlocked populations could hinder consistent adaptation, limiting genetic parallelism. Understanding the degree of genetic parallelism underlying differentiation has implications for both the predictability of evolution and management practices. We employed an 87k single-nucleotide polymorphism (SNP) array to examine the genetic characteristics of landlocked and anadromous Arctic char (Salvelinus alpinus) populations from five drainages within Labrador, Canada. One gene was detected as an outlier between sympatric, size-differentiated morphs in each of two landlocked lakes. While no single locus differentiated all replicate pairs of landlocked and anadromous populations, several SNPs, genes, and paralogs were consistently detected as outliers in at least 70% of these pairwise comparisons. A significant C-score suggested that the amount of shared outlier SNPs across all paired landlocked and anadromous populations was greater than expected by chance. Our results indicate that despite their isolation, selection due to the loss of diadromy may drive consistent genetic responses in landlocked populations.

Harnessing the power of multi-omics data for predicting climate change response.

Predicting how species will respond to future climate change is of central importance in the midst of the global biodiversity crisis, and recent work has demonstrated the utility of population genomics for improving these predictions. Here, we suggest a broadening of the approach to include other types of genomic variants that play an important role in adaptation, like structural (e.g. copy number variants) and epigenetic variants (e.g. DNA methylation). These data could provide additional power for forecasting response, especially in weakly structured or panmictic species. Incorporating structural and epigenetic variation into estimates of climate change vulnerability, or maladaptation, may not only improve prediction power but also provide insight into the molecular mechanisms underpinning species' response to climate change.

Genomic basis of deep-water adaptation in Arctic Charr (Salvelinus alpinus) morphs.

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (FST  = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

Mimicry and mitonuclear discordance in nudibranchs: New insights from exon capture phylogenomics.

Phylogenetic inference and species delimitation can be challenging in taxonomic groups that have recently radiated and where introgression produces conflicting gene trees, especially when species delimitation has traditionally relied on mitochondrial data and color pattern. Chromodoris, a genus of colorful and toxic nudibranch in the Indo-Pacific, has been shown to have extraordinary cryptic diversity and mimicry, and has recently radiated, ultimately complicating species delimitation. In these cases, additional genome-wide data can help improve phylogenetic resolution and provide important insights about evolutionary history. Here, we employ a transcriptome-based exon capture approach to resolve Chromodoris phylogeny with data from 2,925 exons and 1,630 genes, derived from 15 nudibranch transcriptomes. We show that some previously identified mimics instead show mitonuclear discordance, likely deriving from introgression or mitochondrial capture, but we confirm one "pure" mimic in Western Australia. Sister-species relationships and species-level entities were recovered with high support in both concatenated maximum likelihood (ML) and summary coalescent phylogenies, but the ML topologies were highly variable while the coalescent topologies were consistent across datasets. Our work also demonstrates the broad phylogenetic utility of 149 genes that were previously identified from eupulmonate gastropods. This study is one of the first to (a) demonstrate the efficacy of exon capture for recovering relationships among recently radiated invertebrate taxa, (b) employ genome-wide nuclear markers to test mimicry hypotheses in nudibranchs and (c) provide evidence for introgression and mitochondrial capture in nudibranchs.

Limited genetic parallelism underlies recent, repeated incipient speciation in geographically proximate populations of an Arctic fish (Salvelinus alpinus).

The genetic underpinnings of incipient speciation, including the genomic mechanisms which contribute to morphological and ecological differentiation and reproductive isolation, remain poorly understood. The repeated evolution of consistently, phenotypically distinct morphs of Arctic Charr (Salvelinus alpinus) within the Quaternary period offer an ideal model to study the repeatability of evolution at the genomic level. Sympatric morphs of Arctic Charr are found across this species' circumpolar distribution. However, the specific genetic mechanisms driving this morph differentiation are largely unknown despite the cultural and economic importance of the anadromous morph. We used a newly designed 87k SNP chip to investigate the character and consistency of the genomic differences among sympatric morphs within three recently deglaciated and geographically proximate lakes in Labrador, Canada. We found genetically distinct small and large morph Arctic Charr in all three lakes consistent with resident and anadromous morphs, respectively. A degree of reproductive isolation among sympatric morphs is likely given genome-wide distributions of outlier SNPs and high genome-wide FST s. Across all lakes, outlier SNPs were largely nonoverlapping suggesting a lack of genetic parallelism driving morph differentiation. Alternatively, several genes and paralogous copies of the same gene consistently differentiated morphs across multiple lakes suggesting their importance to the manifestation of morphs. Our results confirm the utility of Arctic Charr as a model for investigating the predictability of evolution and support the importance of both genetic parallelism and nonparallelism to the incipient speciation of Arctic Charr morphs.

Resolving fine-scale population structure and fishery exploitation using sequenced microsatellites in a northern fish.

The resiliency of populations and species to environmental change is dependent on the maintenance of genetic diversity, and as such, quantifying diversity is central to combating ongoing widespread reductions in biodiversity. With the advent of next-generation sequencing, several methods now exist for resolving fine-scale population structure, but the comparative performance of these methods for genetic assignment has rarely been tested. Here, we evaluate the performance of sequenced microsatellites and a single nucleotide polymorphism (SNP) array to resolve fine-scale population structure in a critically important salmonid in north eastern Canada, Arctic Charr (Salvelinus alpinus). We also assess the utility of sequenced microsatellites for fisheries applications by quantifying the spatial scales of movement and exploitation through genetic assignment of fishery samples to rivers of origin and comparing these results with a 29-year tagging dataset. Self-assignment and simulation-based analyses of 111 genome-wide microsatellite loci and 500 informative SNPs from 28 populations of Arctic Charr in north-eastern Canada identified largely river-specific genetic structure. Despite large differences (~4X) in the number of loci surveyed between panels, mean self-assignment accuracy was similar with the microsatellite loci and the SNP panel (>90%). Subsequent analysis of 996 fishery-collected samples using the microsatellite panel revealed that larger rivers contribute greater numbers of individuals to the fishery and that coastal fisheries largely exploit individuals originating from nearby rivers, corroborating results from traditional tagging experiments. Our results demonstrate the efficacy of sequence-based microsatellite genotyping to advance understanding of fine-scale population structure and harvest composition in northern and understudied species.

A reference library for Canadian invertebrates with 1.5 million barcodes, voucher specimens, and DNA samples.

The reliable taxonomic identification of organisms through DNA sequence data requires a well parameterized library of curated reference sequences. However, it is estimated that just 15% of described animal species are represented in public sequence repositories. To begin to address this deficiency, we provide DNA barcodes for 1,500,003 animal specimens collected from 23 terrestrial and aquatic ecozones at sites across Canada, a nation that comprises 7% of the planet's land surface. In total, 14 phyla, 43 classes, 163 orders, 1123 families, 6186 genera, and 64,264 Barcode Index Numbers (BINs; a proxy for species) are represented. Species-level taxonomy was available for 38% of the specimens, but higher proportions were assigned to a genus (69.5%) and a family (99.9%). Voucher specimens and DNA extracts are archived at the Centre for Biodiversity Genomics where they are available for further research. The corresponding sequence and taxonomic data can be accessed through the Barcode of Life Data System, GenBank, the Global Biodiversity Information Facility, and the Global Genome Biodiversity Network Data Portal.

A newly discovered radiation of endoparasitic gastropods and their coevolution with asteroid hosts in Antarctica.

BACKGROUND: Marine invertebrates are abundant and diverse on the continental shelf in Antarctica, but little is known about their parasitic counterparts. Endoparasites are especially understudied because they often possess highly modified body plans that pose problems for their identification. Asterophila, a genus of endoparasitic gastropod in the family Eulimidae, forms cysts in the arms and central discs of asteroid sea stars. There are currently four known species in this genus, one of which has been described from the Antarctic Peninsula (A. perknasteri). This study employs molecular and morphological data to investigate the diversity of Asterophila in Antarctica and explore cophylogenetic patterns between host and parasite. RESULTS: A maximum-likelihood phylogeny of Asterophila and subsequent species-delimitation analysis uncovered nine well-supported putative species, eight of which are new to science. Most Asterophila species were found on a single host species, but four species were found on multiple hosts from one or two closely related genera, showing phylogenetic conservatism of host use. Both distance-based and event-based cophylogenetic analyses uncovered a strong signal of coevolution in this system, but most associations were explained by non-cospeciation events. DISCUSSION: The prevalence of duplication and host-switching events in Asterophila and its asteroid hosts suggests that synchronous evolution may be rare even in obligate endoparasitic systems. The apparent restricted distribution of Asterophila from around the Scotia Arc may be an artefact of concentrated sampling in the area and a low obvious prevalence of infection. Given the richness of parasites on a global scale, their role in promoting host diversification, and the threat of their loss through coextinction, future work should continue to investigate parasite diversity and coevolution in vulnerable ecosystems.

Flexible colour patterns obscure identification and mimicry in Indo-Pacific Chromodoris nudibranchs (Gastropoda: Chromodorididae).

Chromodoris is a genus of colourful nudibranchs that feed on sponges and is found across the Indo-Pacific. While this was once the most diverse chromodorid genus, recent work has shown that the genus should be restricted to a monophyletic lineage that contains only 22 species, all of which exhibit black pigmentation and planar spawning behaviour. Earlier phylogenies of this group are poorly resolved and thus additional work is needed to clarify species boundaries within Chromodoris. This study presents a maximum-likelihood phylogeny based on mitochondrial loci (COI, 16S) for 345 Chromodoris specimens, including data from 323 new specimens and 22 from GenBank, from across the Indo-Pacific. Species hypotheses and phylogenetic analysis uncovered 39 taxa in total containing 18 undescribed species, with only five of 39 taxa showing stable colour patterns and distinct morphotypes. This study also presents the first evidence for regional mimicry in this genus, with C. colemani and C. joshi displaying geographically-based variation in colour patterns which appear to match locally abundant congenerics, highlighting the flexibility of these colour patterns in Chromodoris nudibranchs. The current phylogeny contains short branch lengths, polytomies and poor support at interior nodes, which is indicative of a recent radiation. As such, future work will employ a transcriptome-based exon capture approach for resolving the phylogeny of this group. In all, this study included 21 of the 22 described species in the Chromodoris sensu stricto group with broad sampling coverage from across the Indo-Pacific, constituting the most comprehensive sampling of this group to date. This work highlights several cases of undocumented diversity, ultimately expanding our knowledge of species boundaries in this group, while also demonstrating the limitations of colour patterns for species identification in this genus.

Exploring Canadian Echinoderm Diversity through DNA Barcodes.

DNA barcoding has proven an effective tool for species identification in varied groups of marine invertebrates including crustaceans, molluscs, polychaetes and echinoderms. In this study, we further validate its utility by analyzing almost half of the 300 species of Echinodermata known from Canadian waters. COI sequences from 999 specimens were assigned to 145 BINs. In most cases, species discrimination was straightforward due to the large difference (25-fold) between mean intra- (0.48%) and inter- (12.0%) specific divergence. Six species were flagged for further taxonomic investigation because specimens assigned to them fell into two or three discrete sequence clusters. The potential influence of larval dispersal capacity and glacial events on patterns of genetic diversity is discussed for 19 trans-oceanic species. Although additional research is needed to clarify biogeographic patterns and resolve taxonomic questions, this study represents an important step in the assembly of a DNA barcode library for all Canadian echinoderms, a valuable resource for future biosurveillance programs.

Patterns of DNA barcode variation in Canadian marine molluscs.

BACKGROUND: Molluscs are the most diverse marine phylum and this high diversity has resulted in considerable taxonomic problems. Because the number of species in Canadian oceans remains uncertain, there is a need to incorporate molecular methods into species identifications. A 648 base pair segment of the cytochrome c oxidase subunit I gene has proven useful for the identification and discovery of species in many animal lineages. While the utility of DNA barcoding in molluscs has been demonstrated in other studies, this is the first effort to construct a DNA barcode registry for marine molluscs across such a large geographic area. METHODOLOGY/PRINCIPAL FINDINGS: This study examines patterns of DNA barcode variation in 227 species of Canadian marine molluscs. Intraspecific sequence divergences ranged from 0-26.4% and a barcode gap existed for most taxa. Eleven cases of relatively deep (>2%) intraspecific divergence were detected, suggesting the possible presence of overlooked species. Structural variation was detected in COI with indels found in 37 species, mostly bivalves. Some indels were present in divergent lineages, primarily in the region of the first external loop, suggesting certain areas are hotspots for change. Lastly, mean GC content varied substantially among orders (24.5%-46.5%), and showed a significant positive correlation with nearest neighbour distances. CONCLUSIONS/SIGNIFICANCE: DNA barcoding is an effective tool for the identification of Canadian marine molluscs and for revealing possible cases of overlooked species. Some species with deep intraspecific divergence showed a biogeographic partition between lineages on the Atlantic, Arctic and Pacific coasts, suggesting the role of Pleistocene glaciations in the subdivision of their populations. Indels were prevalent in the barcode region of the COI gene in bivalves and gastropods. This study highlights the efficacy of DNA barcoding for providing insights into sequence variation across a broad taxonomic group on a large geographic scale.