Harnessing genomics to fast-track genetic improvement in aquaculture.

Aquaculture is the fastest-growing farmed food sector and will soon become the primary source of fish and shellfish for human diets. In contrast to crop and livestock production, aquaculture production is derived from numerous, exceptionally diverse species that are typically in the early stages of domestication. Genetic improvement of production traits via well-designed, managed breeding programmes has great potential to help meet the rising seafood demand driven by human population growth. Supported by continuous advances in sequencing and bioinformatics, genomics is increasingly being applied across the broad range of aquaculture species and at all stages of the domestication process to optimize selective breeding. In the future, combining genomic selection with biotechnological innovations, such as genome editing and surrogate broodstock technologies, may further expedite genetic improvement in aquaculture.

What do the terms resistance, tolerance, and resilience mean in the case of Ostrea edulis infected by the haplosporidian parasite Bonamia ostreae.

The decline of the European flat oyster Ostrea edulis represents a loss to European coastal economies both in terms of food security and by affecting the Good Environmental Status of the marine environment as set out by the European Council's Marine Strategy Framework Directive (2008/56/EC). Restoration of O. edulis habitat is being widely discussed across Europe, addressing key challenges such as the devastating impact of the haplosporidian parasite Bonamia ostreae. The use of resistant, tolerant, or resilient oysters as restoration broodstock has been proposed by restoration practitioners, but the definitions and implications of these superficially familiar terms have yet to be defined and agreed by all stakeholders. This opinion piece considers the challenges of differentiating Bonamia resistance, tolerance, and resilience; challenges which impede the adoption of robust definitions. We argue that, disease-resistance is reduced susceptibility to infection by the parasite, or active suppression of the parasites ability to multiply and proliferate. Disease-tolerance is the retention of fitness and an ability to neutralise the virulence of the parasite. Disease-resilience is the ability to recover from illness and, at population level, tolerance could be interpreted as resilience. We concede that further work is required to resolve practical uncertainty in applying these definitions, and argue for a collaboration of experts to achieve consensus. Failure to act now might result in the future dispersal of this disease into new locations and populations, because robust definitions are important components of regulatory mechanisms that underpin marine management.

A Single-Tube HNB-Based Loop-Mediated Isothermal Amplification for the Robust Detection of the Ostreid herpesvirus 1.

The Ostreid herpesvirus 1 species affects shellfish, contributing significantly to high economic losses during production. To counteract the threat related to mortality, there is a need for the development of novel point-of-care testing (POCT) that can be implemented in aquaculture production to prevent disease outbreaks. In this study, a simple, rapid and specific colorimetric loop-mediated isothermal amplification (LAMP) assay has been developed for the detection of Ostreid herpesvirus1 (OsHV-1) and its variants infecting Crassostrea gigas (C. gigas). The LAMP assay has been optimized to use hydroxynaphthol blue (HNB) for visual colorimetric distinction of positive and negative templates. The effect of an additional Tte UvrD helicase enzyme used in the reaction was also evaluated with an improved reaction time of 10 min. Additionally, this study provides a robust workflow for optimization of primers for uncultured viruses using designed target plasmid when DNA availability is limited.

'Candidatus Aquirickettsiella gammari' (Gammaproteobacteria: Legionellales: Coxiellaceae): A bacterial pathogen of the freshwater crustacean Gammarus fossarum (Malacostraca: Amphipoda).

Invasive and non-native species can pose risks to vulnerable ecosystems by co-introducing bacterial pathogens. Alternatively, co-introduced bacterial pathogens may regulate invasive population size and invasive traits. We describe a novel candidate genus and species of bacteria ('Candidatus Aquirickettsiella gammari') found to infect Gammarus fossarum, from its native range in Poland. The bacterium develops intracellularly within the haemocytes and cells of the musculature, hepatopancreas, connective tissues, nervous system and gonad of the host. The developmental cycle of 'Candidatus Aquirickettsiella gammari' includes an elementary body (496.73 nm ±â€¯37.56 nm in length, and 176.89 nm ±â€¯36.29 nm in width), an elliptical, condensed spherical stage (737.61 nm ±â€¯44.51 nm in length and 300.07 nm ±â€¯44.02 nm in width), a divisional stage, and a spherical initial body (1397.59 nm ±â€¯21.26 nm in diameter). We provide a partial genome for 'Candidatus Aquirickettsiella gammari', which clades phylogenetically alongside environmental 16S rRNA sequences from aquatic habitats, and bacterial symbionts from aquatic isopods (Asellus aquaticus), grouping separately from the Rickettsiella, a genus that includes bacterial pathogens of terrestrial insects and isopods. Increased understanding of the diversity of symbionts carried by G. fossarum identifies those that might regulate host population size, or those that could pose a risk to native species in the invasive range. Identification of 'Candidatus Aquirickettsiella gammari' and its potential for adaptation as a biological control agent is explored.

New Invasive Nemertean Species (Cephalothrix Simula) in England with High Levels of Tetrodotoxin and a Microbiome Linked to Toxin Metabolism.

The marine nemertean Cephalothrix simula originates from the Pacific Ocean but in recent years has been discovered in northern Europe. The species has been associated with high levels of the marine neurotoxin Tetrodotoxin, traditionally associated with Pufferfish Poisoning. This study reports the first discovery of two organisms of C. simula in the UK, showing the geographical extent of this species is wider than originally described. Species identification was initially conducted morphologically, with confirmation by Cox 1 DNA sequencing. 16S gene sequencing enabled the taxonomic assignment of the microbiome, showing the prevalence of a large number of bacterial genera previously associated with TTX production including Alteromonas, Vibrio and Pseudomonas. LC-MS/MS analysis of the nemertean tissue revealed the presence of multiple analogues of TTX, dominated by the parent TTX, with a total toxin concentration quantified at 54 µg TTX per g of tissue. Pseudomonas luteola isolated from C. simula, together with Vibrio alginolyticus from the native nemertean Tubulanus annulatus, were cultured at low temperature and both found to contain TTX. Overall, this paper confirms the high toxicity of a newly discovered invasive nemertean species with links to toxin-producing marine bacteria and the potential risk to human safety. Further work is required to assess the geographical extent and toxicity range of C. simula along the UK coast in order to properly gauge the potential impacts on the environment and human safety.

Dynamics of nuclear receptor gene expression during Pacific oyster development.

BACKGROUND: Nuclear receptors are a highly conserved set of ligand binding transcription factors, with essential roles regulating aspects of vertebrate and invertebrate biology alike. Current understanding of nuclear receptor regulated gene expression in invertebrates remains sparse, limiting our ability to elucidate gene function and the conservation of developmental processes across phyla. Here, we studied nuclear receptor expression in the early life stages of the Pacific oyster, Crassostrea gigas, to identify at which specific key stages nuclear receptors are expressed RESULTS: We used quantitative RT-PCR to determine the expression profiles of 34 nuclear receptors, revealing three developmental key stages, during which nuclear receptor expression is dynamically regulated: embryogenesis, mid development from gastrulation to trochophore larva, and late larval development prior to metamorphosis. Clustering of nuclear receptor expression patterns demonstrated that transcriptional regulation was not directly related to gene phylogeny, suggesting closely related genes may have distinct functions. Expression of gene homologs of vertebrate retinoid receptors suggests participation in organogenesis and shell-formation, as they are highly expressed at the gastrulation and trochophore larval initial shell formation stages. The ecdysone receptor homolog showed high expression just before larval settlement, suggesting a potential role in metamorphosis. CONCLUSION: Throughout early oyster development nuclear receptors exhibited highly dynamic expression profiles, which were not confined by gene phylogeny. These results provide fundamental information on the presence of nuclear receptors during key developmental stages, which aids elucidation of their function in the developmental process. This understanding is essential as ligand sensing nuclear receptors can be disrupted by xenobiotics, a mode of action through which anthropogenic environmental pollutants have been found to mediate effects.

The nuclear receptor gene family in the Pacific oyster, Crassostrea gigas, contains a novel subfamily group.

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.

A chromosome-scale Mytilus edulis genome assembly for aquaculture, marine ecology, and evolution.

The smooth-shelled blue mussel, Mytilus edulis is part of the Mytilus species complex, encompassing at least three putative species: M. edulis, Mytilus galloprovincialis, and Mytilus trossulus. These three species occur on both sides of the Atlantic and hybridize in nature, and both M. edulis and M. galloprovincialis are important aquaculture species. They are also invasive species in many parts of the world. Here, we present a chromosome-level assembly of M. edulis. We used a combination of PacBio sequencing and Dovetail's Omni-C technology to generate an assembly with 14 long scaffolds containing 94% of the predicted length of the M. edulis genome (1.6 out of 1.7 Gb). Assembly statistics were as follows: total length = 1.65 Gb, N50 = 116 Mb, L50 = 7, and L90 = 13. BUSCO analysis showed 92.55% eukaryote BUSCOs identified. AB-Initio annotation using RNA-seq from mantle, gills, muscle, and foot predicted 47,128 genes. These gene models were combined with IsoSeq validation resulting in 45,379 full CDS protein sequences and 129,708 isoforms. Using GBS and shotgun sequencing, we also sequenced several eastern Canadian populations of Mytilus to characterize single-nucleotide as well as structural variance. This high-quality genome for M. edulis provides a platform to develop tools that can be used in breeding, molecular ecology and evolution to address questions of both commercial and environmental perspectives.

Laboratory Replication of Ostreid Herpes Virus (OsHV-1) Using Pacific Oyster Tissue Explants.

Pacific oysters (Crassostrea or Magallana gigas) are one of the most economically important aquaculture species globally. Over the past two decades, ostreid herpesvirus (OsHV-1) has become a major pathogen of cultured Pacific oysters, resulting in widespread mortality with a global distribution. Experimental use of OsHV-1 is challenging for many reasons, including both complexity of host-pathogen dynamics and a lack of functioning model systems. The goal of this study was to improve the tools available for working with OsHV-1 in both whole animals and in tissue explants established from oysters maintained in controlled laboratory conditions. Tissue explants were taken from oysters originating from two different sources that have different levels of mortality in experimental OsHV-1 infections and were exposed to OsHV-1. A whole-animal infection experiment was run concurrently as a comparison. Quantitative PCR and electron microscopy were used to confirm that the explants were capable of replicating OsHV-1. Furthermore, the quantitative PCR results suggest that the source of the oysters was significant in determining the outcome of infection in the explants, supporting the validity of the explant model for OsHV-1 infection. This tissue explant approach for studying OsHV-1 allows for the control of confounding factors in the disease outcome that is not possible in whole-animal experiments, providing a new tool for the study of OsHV-1 in Pacific oysters.

Cryopreservation and storage of mussel (Mytilus spp.) haemocytes for latent analysis by the Comet assay.

Estuarine and coastal habitats are known to be polluted by a range of chemical contaminants from both industrial and domestic sources. Blue mussels (Mytilus spp.), which inhabit these areas, are widely used as bio-indicators in eco-toxicological studies, because of their sedentary nature and their ability to bio-accumulate contaminants. The analysis of DNA damage in mussel haemocytes is a valuable tool for biomonitoring but sampling issues related to storage, handling and transportation have often limited its application in large-scale monitoring programmes. This study uses a trial and error method to evaluate and validate a suitable protocol for cryopreservation of mussel haemocytes, thereby allowing material collected in the field to be analysed later under controlled laboratory conditions. Three different cell-culture media, i.e. Leibovitz-15, Hank's balanced salt solution and mussel physiological saline, along with four different cryoprotectants, i.e. dimethyl sulphoxide (10% and 20%), 1,2-propanediol (10%), ethylene glycol (10%) and glycerol (10%) were tested to assess their suitability for cryopreservation of mussel haemocytes for analysis in the comet assay. Experimental studies where mussel haemocytes were also exposed to UV radiation or benzo(a)pyrene were conducted in order to mimic environmental stresses and to verify the effectiveness of newly defined cryopreservation protocols. The comet assay was used to demonstrate that mussel haemocytes could be preserved at cryogenic temperatures for a month without altering levels of DNA damage, which could possibly be used for lab or field studies where time constraints or facilities do not allow instant analysis.

Oxidative DNA damage may not mediate Ni-induced genotoxicity in marine mussels: assessment of genotoxic biomarkers and transcriptional responses of key stress genes.

Nickel (Ni) is a known carcinogenic and mutagenic compound and an important contaminant of aquatic environments. Ni toxicity and its potential impact on aquatic organisms are, however, not well understood. This study used an integrated approach to evaluate genotoxic effects, tissue-specific accumulation and transcriptional profiles of key genes in mussels, Mytilus galloprovincialis, exposed to a range of concentrations of Ni. The genotoxic effects assessed were total and oxidative DNA damage (DNA strand breaks measured using the enzyme modified comet assay), and induction of micronuclei (MN; clastogenic and/or aneugenic effects) using haemocytes as the target cells. Six genes (pgp, mt10, mt20, sod, hsp70 and gst) were selected for transcriptional analysis in the gills based on their key role in the stress response. Following exposure to sublethal concentrations of Ni (0-3600µgL(-1)) for 5 days, mussel haemocytes showed significant genotoxicity at >1800µgL(-1) (4-fold increase for DNA strand breaks and 3-fold increase for MN induction). There was no significant difference between buffer (control) and enzyme treatments which target oxidised DNA bases (formamidopyrimidine glycosylase or endonuclease IIII). This suggested that, in haemocytes, oxidative DNA damage is not a major mechanism for Ni-induced genotoxicity. The expression of mt20 and gst genes in gill was up-regulated at genotoxic concentrations, whilst pgp expression was markedly up-regulated, particularly at 18µgL(-1) Ni (19-fold increase). Pearson's correlation analysis revealed significant associations between % tail DNA and MN induction in haemocytes (r=0.88, p<0.05), and between Ni accumulation in foot (r=0.47, p<0.05) and digestive gland (r=0.41, p<0.05) and induction of MN in the haemocytes. Our results are the first to suggest that Ni-induced genotoxicity in mussel haemocytes may not be a result of oxidative DNA damage, and that multixenobiotic resistance (MXR) may play an important role in Ni detoxification in this species.

SNP discovery and genetic structure in blue mussel species using low coverage sequencing and a medium density 60 K SNP-array.

Blue mussels from the genus Mytilus are an abundant component of the benthic community, found in the high latitude habitats. These foundation species are relevant to the aquaculture industry, with over 2 million tonnes produced globally each year. Mussels withstand a wide range of environmental conditions and species from the Mytilus edulis complex readily hybridize in regions where their distributions overlap. Significant effort has been made to investigate the consequences of environmental stress on mussel physiology, reproductive isolation, and local adaptation. Yet our understanding on the genomic mechanisms underlying such processes remains limited. In this study, we developed a multi species medium-density 60 K SNP-array including four species of the Mytilus genus. SNPs included in the platform were called from 138 mussels from 23 globally distributed mussel populations, sequenced using a whole-genome low coverage approach. The array contains polymorphic SNPs which capture the genetic diversity present in mussel populations thriving across a gradient of environmental conditions (~59 K SNPs) and a set of published and validated SNPs informative for species identification and for diagnosis of transmissible cancer (610 SNPs). The array will allow the consistent genotyping of individuals, facilitating the investigation of ecological and evolutionary processes in these taxa. The applications of this array extend to shellfish aquaculture, contributing to the optimization of this industry via genomic selection of blue mussels, parentage assignment, inbreeding assessment and traceability. Further applications such as genome wide association studies (GWAS) for key production traits and those related to environmental resilience are especially relevant to safeguard aquaculture production under climate change.

Evaluation of interactive effects of phosphorus-32 and copper on marine and freshwater bivalve mollusks.

PURPOSE: Contaminants seldom occur in isolation in the aquatic environment. While pollution of coastal and inland water bodies has received considerable attention to date, there is limited information on potential interactive effects between radionuclides and metals. Whether by accidental or controlled release, such contaminants co-exist in aquatic ecosystems and can pose an enhanced threat to biota. Using a range of biological responses, the study aimed to evaluate relative interactive effects on representative freshwater and marine bivalve species. METHODS: An integrated, multi-biomarker approach was adopted to investigate response to copper (Cu, 18 µg L-1), a known environmentally relevant genotoxic metal and differing concentrations of phosphorus-32 (32P; 0.1 and 1 mGy d-1), alone and in combination in marine (Mytilus galloprovincialis) and freshwater (Dreissena polymorpha) mussels. Genetic and molecular biomarkers were determined post-exposure and included DNA damage (as measured by the comet assay), micronuclei (MN) formation, γ-H2AX foci induction and the expression of key stress-related genes (i.e. hsp70/90, sod, cat, gst). RESULTS: Overall, using a tissue-specific (i.e. gill and digestive gland) approach, genotoxic response was reflective of exposures where Cu had a slight additive effect on 32P-induced damage across the species (but not all), cell types and dose rates. Multivariate analysis found significant correlations between comet and γ-H2AX assays, across both the tissues. Transcriptional expression of selected genes were generally unaltered in response to contaminant exposures, independent of species or tissues. CONCLUSIONS: Our study is the first to explore the interactive effects of ionizing radiation (IR) and Cu on two bivalve species representing two ecological habitats. The complexity of IR-metal interactions demonstrate that extrapolation of findings obtained from single stressor studies into field conditions could be misrepresentative of real-world environments. In turn, environmental protective strategies deemed suitable in protecting biota from a single, isolated stressor may not be wholly adequate.

Two parallel chromosome-level reference genomes to support restoration and aquaculture of European flat oyster Ostrea edulis.

This volume of Evolutionary Applications sees the publication of two genomes for the European native flat oyster Ostrea edulis, a species of significant evolutionary, ecological and commercial value. Each is a highly contiguous chromosome-level assembly from individuals of different genetic backgrounds, which have been benchmarked against one another. This situation has resulted from the serendipitous discovery that two independent research groups were both deep into the process of building, annotating and investigating separately produced assemblies. Due to constraints with funder requirements and the need to recognize early career researchers for their work, alongside the technical challenge of integrating assemblies from two very different genomes, there was limited capacity to merge the sequences into one publication at the stage of discovery. This issue is likely to become very common over the next few years until the technologies for working with multiple genomes at once, for example, graph genomes, become commonplace in nonmodel species. Consequently, both of our teams have decided to collaborate rather than compete, recognizing the benefit to copublishing two separate genome resources for the research community, each with distinct scientific investigations, and working collaboratively to benchmark the assemblies.

Chromosome-level reference genome for European flat oyster (Ostrea edulis L.).

The European flat oyster (Ostrea edulis L.) is a bivalve naturally distributed across Europe, which was an integral part of human diets for centuries, until anthropogenic activities and disease outbreaks severely reduced wild populations. Despite a growing interest in genetic applications to support population management and aquaculture, a reference genome for this species is lacking to date. Here, we report a chromosome-level assembly and annotation for the European Flat oyster genome, generated using Oxford Nanopore, Illumina, Dovetail OmniC™ proximity ligation and RNA sequencing. A contig assembly (N50: 2.38 Mb) was scaffolded into the expected karyotype of 10 pseudochromosomes. The final assembly is 935.13 Mb, with a scaffold-N50 of 95.56 Mb, with a predicted repeat landscape dominated by unclassified elements specific to O. edulis. The assembly was verified for accuracy and completeness using multiple approaches, including a novel linkage map built with ddRAD-Seq technology, comprising 4016 SNPs from four full-sib families (eight parents and 163 F1 offspring). Annotation of the genome integrating multitissue transcriptome data, comparative protein evidence and ab-initio gene prediction identified 35,699 protein-coding genes. Chromosome-level synteny was demonstrated against multiple high-quality bivalve genome assemblies, including an O. edulis genome generated independently for a French O. edulis individual. Comparative genomics was used to characterize gene family expansions during Ostrea evolution that potentially facilitated adaptation. This new reference genome for European flat oyster will enable high-resolution genomics in support of conservation and aquaculture initiatives, and improves our understanding of bivalve genome evolution.

A single genomic region involving a putative chromosome rearrangement in flat oyster (Ostrea edulis) is associated with differential host resilience to the parasite Bonamia ostreae.

European flat oyster (Ostrea edulis) is an ecologically and economically important marine bivalve, that has been severely affected by the intracellular parasite Bonamia ostreae. In this study, a flat oyster SNP array (~14,000 SNPs) was used to validate previously reported outlier loci for divergent selection associated with B. ostreae exposure in the Northeast Atlantic Area. A total of 134 wild and hatchery individuals from the North Sea, collected in naïve (NV) and long-term affected (LTA) areas, were analysed. Genetic diversity and differentiation were related to the sampling origin (wild vs. hatchery) when using neutral markers, and to bonamiosis status (NV vs. LTA) when using outlier loci for divergent selection. Two genetic clusters appeared intermingled in all sampling locations when using outlier loci, and their frequency was associated with their bonamiosis status. When both clusters were compared, outlier data sets showed high genetic divergence (F ST > 0.25) unlike neutral loci (F ST not ≠ 0). Moreover, the cluster associated with LTA samples showed much higher genetic diversity and significant heterozygote excess with outlier loci, but not with neutral data. Most outliers mapped on chromosome 8 (OE-C8) of the flat oyster genome, supporting a main genomic region underlying resilience to bonamiosis. Furthermore, differentially expressed genes previously reported between NV and LTA strains showed higher mapping density on OE-C8. A range of relevant immune functions were specifically enriched among genes annotated on OE-C8, providing hypotheses for resilience mechanisms to an intracellular parasite. The results suggest that marker-assisted selection could be applied to breed resilient strains of O. edulis to bonamiosis, if lower parasite load and/or higher viability of the LTA genetic cluster following B. ostreae infection is demonstrated.

Genome-Wide Association and Genomic Prediction of Growth Traits in the European Flat Oyster (Ostrea edulis).

The European flat oyster (Ostrea edulis) is a bivalve mollusc that was once widely distributed across Europe and represented an important food resource for humans for centuries. Populations of O. edulis experienced a severe decline across their biogeographic range mainly due to overexploitation and disease outbreaks. To restore the economic and ecological benefits of European flat oyster populations, extensive protection and restoration efforts are in place within Europe. In line with the increasing interest in supporting restoration and oyster farming through the breeding of stocks with enhanced performance, the present study aimed to evaluate the potential of genomic selection for improving growth traits in a European flat oyster population obtained from successive mass-spawning events. Four growth-related traits were evaluated: total weight (TW), shell height (SH), shell width (SW) and shell length (SL). The heritability of the growth traits was in the low-moderate range, with estimates of 0.45, 0.37, 0.22, and 0.32 for TW, SH, SW and SL, respectively. A genome-wide association analysis revealed a largely polygenic architecture for the four growth traits, with two distinct QTLs detected on chromosome 4. To investigate whether genomic selection can be implemented in flat oyster breeding at a reduced cost, the utility of low-density SNP panels was assessed. Genomic prediction accuracies using the full density panel were high (> 0.83 for all traits). The evaluation of the effect of reducing the number of markers used to predict genomic breeding values revealed that similar selection accuracies could be achieved for all traits with 2K SNPs as for a full panel containing 4,577 SNPs. Only slight reductions in accuracies were observed at the lowest SNP density tested (i.e., 100 SNPs), likely due to a high relatedness between individuals being included in the training and validation sets during cross-validation. Overall, our results suggest that the genetic improvement of growth traits in oysters is feasible. Nevertheless, and although low-density SNP panels appear as a promising strategy for applying GS at a reduced cost, additional populations with different degrees of genetic relatedness should be assessed to derive estimates of prediction accuracies to be expected in practical breeding programmes.

Ancestral Physical Stress and Later Immune Gene Family Expansions Shaped Bivalve Mollusc Evolution.

Bivalve molluscs comprise 20,000 species occupying a wide diversity of marine habitats. As filter feeders and detritivores they act as ecosystem engineers clarifying water, creating reefs, and protecting coastlines. The global decline of natural oyster reefs has led to increased restoration efforts in recent years. Bivalves also play an important role in global food security contributing to >20% of worldwide aquaculture production. Despite this importance, relatively little is known about bivalve evolutionary adaptation strategies. Difficulties previously associated with highly heterozygous and repetitive regions of bivalve genomes have been overcome by long-read sequencing, enabling the generation of accurate bivalve assemblies. With these resources we have analyzed the genomes of 32 species representing each molluscan class, including 15 bivalve species, to identify gene families that have undergone expansion during bivalve evolution. Gene family expansions across bivalve genomes occur at the point of evolutionary pressures. We uncovered two key factors that shape bivalve evolutionary history: expansion of bivalvia into environmental niches with high stress followed by later exposure to specific pathogenic pressures. The conserved expansion of protein recycling gene families we found across bivalvia is mirrored by adaptations to a sedentary lifestyle seen in plants. These results reflect the ability of bivalves to tolerate high levels of environmental stress and constant exposure to pathogens as filter feeders. The increasing availability of accurate genome assemblies will provide greater resolution to these analyses allowing further points of evolutionary pressure to become clear in other understudied taxa and potentially different populations of a single species.

A chromosome-level genome assembly for the Pacific oyster Crassostrea gigas.

BACKGROUND: 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. FINDINGS: 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. CONCLUSIONS: 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.

Potential of genomic technologies to improve disease resistance in molluscan aquaculture.

Molluscan aquaculture is a major contributor to global seafood production, but is hampered by infectious disease outbreaks that can cause serious economic losses. Selective breeding has been widely used to improve disease resistance in major agricultural and aquaculture species, and has clear potential in molluscs, albeit its commercial application remains at a formative stage. Advances in genomic technologies, especially the development of cost-efficient genomic selection, have the potential to accelerate genetic improvement. However, tailored approaches are required owing to the distinctive reproductive and life cycle characteristics of molluscan species. Transgenesis and genome editing, in particular CRISPR/Cas systems, have been successfully trialled in molluscs and may further understanding and improvement of genetic resistance to disease through targeted changes to the host genome. Whole-organism genome editing is achievable on a much greater scale compared to other farmed species, making genome-wide CRISPR screening approaches plausible. This review discusses the current state and future potential of selective breeding, genomic tools and genome editing approaches to understand and improve host resistance to infectious disease in molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.