This viewpoint paper explores the potential of genomics technology to provide accurate, rapid, and cost efficient observations of the marine environment. The use of such approaches in next generation marine monitoring programs will help achieve the goals of marine legislation implemented world-wide. Genomic methods can yield faster results from monitoring, easier and more reliable taxonomic identification, as well as quicker and better assessment of the environmental status of marine waters. A summary of genomic methods that are ready or show high potential for integration into existing monitoring programs is provided (e.g. qPCR, SNP based methods, DNA barcoding, microarrays, metagenetics, metagenomics, transcriptomics). These approaches are mapped to existing indicators and descriptors and a series of case studies is presented to assess the cost and added value of these molecular techniques in comparison with traditional monitoring systems. Finally, guidelines and recommendations are suggested for how such methods can enter marine monitoring programs in a standardized manner.

The Pacific oyster (Crassostrea gigas) is a bivalve mollusc with vital roles in coastal ecosystems and aquaculture globally. While extensive genomic tools are available for C. gigas, highly contiguous reference genomes are required to support both fundamental and applied research. Herein we report the creation and annotation of a chromosome-level assembly for C. gigas. High-coverage long- and short-read sequence data generated on Pacific Biosciences and Illumina platforms were used to generate an initial assembly, which was then scaffolded into 10 pseudo-chromosomes using both Hi-C sequencing and a high-density linkage map. The assembly has a scaffold N50 of 58.4 Mb and a contig N50 of 1.8 Mb, representing a step advance on the previously published C. gigas assembly. Annotation based on Pacific Biosciences Iso-Seq and Illumina RNA-Seq resulted in identification of ∼30,000 putative protein-coding genes. Annotation of putative repeat elements highlighted an enrichment of Helitron rolling-circle transposable elements, suggesting their potential role in shaping the evolution of the C. gigas genome. This new chromosome-level assembly will be an enabling resource for genetics and genomics studies to support fundamental insight into bivalve biology, as well as for selective breeding of C. gigas in aquaculture.

Pacific oysters are a key aquaculture species globally, and genetic improvement via selective breeding is a major target. Genomic selection has the potential to expedite genetic gain for key target traits of a breeding program, but has not yet been evaluated in oyster. The recent development of SNP arrays for Pacific oyster (Crassostrea gigas) raises the opportunity to test genomic selection strategies for polygenic traits. In this study, a population of 820 oysters (comprising 23 full-sibling families) were genotyped using a medium density SNP array (23 K informative SNPs), and the genetic architecture of growth-related traits [shell height (SH), shell length (SL), and wet weight (WW)] was evaluated. Heritability was estimated to be moderate for the three traits (0.26 ± 0.06 for SH, 0.23 ± 0.06 for SL and 0.35 ± 0.05 for WW), and results of a GWAS indicated that the underlying genetic architecture was polygenic. Genomic prediction approaches were used to estimate breeding values for growth, and compared to pedigree based approaches. The accuracy of the genomic prediction models (GBLUP) outperformed the traditional pedigree approach (PBLUP) by ∼25% for SL and WW, and ∼30% for SH. Further, reduction in SNP marker density had little impact on prediction accuracy, even when density was reduced to a few hundred SNPs. These results suggest that the use of genomic selection in oyster breeding could offer benefits for the selection of breeding candidates to improve complex economic traits at relatively modest cost.

SNP arrays are enabling tools for high-resolution studies of the genetic basis of complex traits in farmed and wild animals. Oysters are of critical importance in many regions from both an ecological and economic perspective, and oyster aquaculture forms a key component of global food security. The aim of our study was to design a combined-species, medium density SNP array for Pacific oyster (Crassostrea gigas) and European flat oyster (Ostrea edulis), and to test the performance of this array on farmed and wild populations from multiple locations, with a focus on European populations. SNP discovery was carried out by whole-genome sequencing (WGS) of pooled genomic DNA samples from eight C. gigas populations, and restriction site-associated DNA sequencing (RAD-Seq) of 11 geographically diverse O. edulis populations. Nearly 12 million candidate SNPs were discovered and filtered based on several criteria, including preference for SNPs segregating in multiple populations and SNPs with monomorphic flanking regions. An Affymetrix Axiom Custom Array was created and tested on a diverse set of samples (n = 219) showing ∼27 K high quality SNPs for C. gigas and ∼11 K high quality SNPs for O. edulis segregating in these populations. A high proportion of SNPs were segregating in each of the populations, and the array was used to detect population structure and levels of linkage disequilibrium (LD). Further testing of the array on three C. gigas nuclear families (n = 165) revealed that the array can be used to clearly distinguish between both families based on identity-by-state (IBS) clustering parental assignment software. This medium density, combined-species array will be publicly available through Affymetrix, and will be applied for genome-wide association and evolutionary genetic studies, and for genomic selection in oyster breeding programs.

Nuclear receptors are a superfamily of transcription factors important in key biological, developmental and reproductive processes. Several of these receptors are ligand- activated and through their ability to bind endogenous and exogenous ligands, are potentially vulnerable to xenobiotics. Molluscs are key ecological species in defining aquatic and terrestrial habitats and are sensitive to xenobiotic compounds in the environment. However, the understanding of nuclear receptor presence, function and xenobiotic disruption in the phylum Mollusca is limited. Here, forty-three nuclear receptor sequences were mined from the genome of the Pacific oyster, Crassostrea gigas. They include members of NR0-NR5 subfamilies, notably lacking any NR6 members. Phylogenetic analyses of the oyster nuclear receptors have been conducted showing the presence of a large novel subfamily group not previously reported, which is named NR1P. Homologues to all previous identified nuclear receptors in other mollusc species have also been determined including the putative heterodimer partner retinoid X receptor, estrogen receptor and estrogen related receptor. C. gigas contains a highly diverse set of nuclear receptors including a novel NR1 group, which provides important information on presence and evolution of this transcription factor superfamily in invertebrates. The Pacific oyster possesses two members of NR3, the sex steroid hormone receptor analogues, of which there are 9 in humans. This provides increasing evidence that steroid ligand specific expansion of this family is deuterostome specific. This new knowledge on divergence and emergence of nuclear receptors in C. gigas provides essential information for studying regulation of molluscan gene expression and the potential effects of xenobiotics.

Ostreid herpesvirus (OsHV) can cause mass mortality events in Pacific oyster aquaculture. While various factors impact on the severity of outbreaks, it is clear that genetic resistance of the host is an important determinant of mortality levels. This raises the possibility of selective breeding strategies to improve the genetic resistance of farmed oyster stocks, thereby contributing to disease control. Traditional selective breeding can be augmented by use of genetic markers, either via marker-assisted or genomic selection. The aim of the current study was to investigate the genetic architecture of resistance to OsHV in Pacific oyster, to identify genomic regions containing putative resistance genes, and to inform the use of genomics to enhance efforts to breed for resistance. To achieve this, a population of ∼1,000 juvenile oysters were experimentally challenged with a virulent form of OsHV, with samples taken from mortalities and survivors for genotyping and qPCR measurement of viral load. The samples were genotyped using a recently-developed SNP array, and the genotype data were used to reconstruct the pedigree. Using these pedigree and genotype data, the first high density linkage map was constructed for Pacific oyster, containing 20,353 SNPs mapped to the ten pairs of chromosomes. Genetic parameters for resistance to OsHV were estimated, indicating a significant but low heritability for the binary trait of survival and also for viral load measures (h2 0.12 – 0.25). A genome-wide association study highlighted a region of linkage group 6 containing a significant QTL affecting host resistance. These results are an important step toward identification of genes underlying resistance to OsHV in oyster, and a step toward applying genomic data to enhance selective breeding for disease resistance in oyster aquaculture.