Predicting the future of our oceans-Evaluating genomic forecasting approaches in marine species

Climate change is restructuring biodiversity on multiple scales and there is a pressing need to understand the downstream ecological and genomic consequences of this change. Recent advancements in the field of eco-evolutionary genomics have sought to include evolutionary processes in forecasting species' responses to climate change (e.g., genomic offset), but to date, much of this work has focused on terrestrial species. Coastal and offshore species, and the fisheries they support, may be even more vulnerable to climate change than their terrestrial counterparts, warranting a critical appraisal of these approaches in marine systems. First, we synthesize knowledge about the genomic basis of adaptation in marine species, and then we discuss the few examples where genomic forecasting has been applied in marine systems. Next, we identify the key challenges in validating genomic offset estimates in marine species, and we advocate for the inclusion of historical sampling data and hindcasting in the validation phase. Lastly, we describe a workflow to guide marine managers in incorporating these predictions into the decision-making process.

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 87k 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.

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 (F_ST  = 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.

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.

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.

Biodiversity inventories in high gear: DNA barcoding facilitates a rapid biotic survey of a temperate nature reserve

BACKGROUND

Comprehensive biotic surveys, or 'all taxon biodiversity inventories' (ATBI), have traditionally been limited in scale or scope due to the complications surrounding specimen sorting and species identification. To circumvent these issues, several ATBI projects have successfully integrated DNA barcoding into their identification procedures and witnessed acceleration in their surveys and subsequent increase in project scope and scale. The Biodiversity Institute of Ontario partnered with the rare Charitable Research Reserve and delegates of the 6th International Barcode of Life Conference to complete its own rapid, barcode-assisted ATBI of an established land trust in Cambridge, Ontario, Canada.

NEW INFORMATION

The existing species inventory for the rare Charitable Research Reserve was rapidly expanded by integrating a DNA barcoding workflow with two surveying strategies - a comprehensive sampling scheme over four months, followed by a one-day bioblitz involving international taxonomic experts. The two surveys resulted in 25,287 and 3,502 specimens barcoded, respectively, as well as 127 human observations. This barcoded material, all vouchered at the Biodiversity Institute of Ontario collection, covers 14 phyla, 29 classes, 117 orders, and 531 families of animals, plants, fungi, and lichens. Overall, the ATBI documented 1,102 new species records for the nature reserve, expanding the existing long-term inventory by 49%. In addition, 2,793 distinct Barcode Index Numbers (BINs) were assigned to genus or higher level taxonomy, and represent additional species that will be added once their taxonomy is resolved. For the 3,502 specimens, the collection, sequence analysis, taxonomic assignment, data release and manuscript submission by 100+ co-authors all occurred in less than one week. This demonstrates the speed at which barcode-assisted inventories can be completed and the utility that barcoding provides in minimizing and guiding valuable taxonomic specialist time. The final product is more than a comprehensive biotic inventory - it is also a rich dataset of fine-scale occurrence and sequence data, all archived and cross-linked in the major biodiversity data repositories. This model of rapid generation and dissemination of essential biodiversity data could be followed to conduct regional assessments of biodiversity status and change, and potentially be employed for evaluating progress towards the Aichi Targets of the Strategic Plan for Biodiversity 2011-2020.