Potential of genomic selection for improvement of resistance to ostreid herpesvirus in Pacific oyster (Crassostrea gigas)

Intraspecific variations in oyster (Magallana gigas) ploidy does not affect physiological responses to microplastic pollution

Recent advances in genetic manipulation such as triploid breeding and artificial selection, have rapidly emerged as valuable hatchery methodologies for enhancing seafood stocks. The Pacific oyster Magallana gigas is a leading aquaculture species worldwide and key ecosystem engineer that has received particular attention in this field of science. In light of the growing recognition of the ecological effects of intraspecific variation, oyster polyploids provide a valuable opportunity to assess whether intraspecific diversity affects physiological responses to environmental stressors. While the responses of diploid and triploid oysters to climate change have been extensively investigated, research on their sensitivity to environmental pollution remains scarce. Here, we assess whether genotypic (i.e., ploidy) variation within Magallana gigas affects physiological responses to microplastic pollution. We show that diploid and triploid M. gigas have similar clearance rates and ingest similar amounts of microplastics under laboratory-controlled condition. In addition, they exhibited similar heart rates after prolonged exposure to microplastic leachates. Our findings suggest that intraspecific variations within M. gigas ploidy does not affect oyster responses to microplastic pollution. However, regardless of ploidy, our work highlights significant adverse effects of microplastic leachates on the heart rate of M. gigas and provides evidence of microplastic ingestion in the laboratory.

Genome-wide allele frequency studies in Pacific oyster families identify candidate genes for tolerance to ostreid herpesvirus 1 (OsHV-1)

BACKGROUND

Host genetics influences the development of infectious diseases in many agricultural animal species. Identifying genes associated with disease development has the potential to make selective breeding for disease tolerance more likely to succeed through the selection of different genes in diverse signaling pathways. In this study, four families of Pacific oysters (Crassostrea gigas) were identified to be segregating for a quantitative trait locus (QTL) on chromosome 8. This QTL was previously found to be associated with basal antiviral gene expression and survival to ostreid herpesvirus 1 (OsHV-1) mortality events in Tomales Bay, California. Individuals from these four families were phenotyped and genotyped in an attempt to find candidate genes associated with the QTL on chromosome 8.

RESULTS

Genome-wide allele frequencies of oysters from each family prior to being planting in Tomales Bay were compared with the allele frequencies of oysters from respective families that survived an OsHV-1 mortality event. Six significant unique QTL were identified in two families in these genome-wide allele frequency studies, all of which were located on chromosome 8. Three QTL were assigned to candidate genes (ABCA1, PIK3R1, and WBP2) that have been previously associated with antiviral innate immunity in vertebrates.

CONCLUSION

The identification of vertebrate antiviral innate immunity genes as candidate genes involved in molluscan antiviral innate immunity reinforces the similarities between the innate immune systems of these two groups. Causal variant identification in these candidate genes will enable future functional studies of these genes in an effort to better understand their antiviral modes of action.

Epigenetic variations are more substantial than genetic variations in rapid adaptation of oyster to Pacific oyster mortality syndrome

Disease emergence is accelerating with global changes. Understanding by which mechanisms host populations can rapidly adapt will be crucial for management practices. Pacific oyster mortality syndrome (POMS) imposes a substantial and recurrent selective pressure on oyster populations, and rapid adaptation may arise through genetics and epigenetics. In this study, we used (epi)genome-wide association mapping to show that oysters differentially exposed to POMS displayed genetic and epigenetic signatures of selection. Consistent with higher resistance to POMS, the genes targeted included many genes in several pathways related to immunity. By combining correlation, DNA methylation quantitative trait loci, and variance partitioning, we revealed that a third of phenotypic variation was explained by interactions between the genetic and epigenetic information, ~14% by the genome, and up to 25% by the epigenome alone. Similar to genetically based adaptation, epigenetic mechanisms notably governing immune responses can contribute substantially to the rapid adaptation of hosts to emerging infectious diseases.

Evaluation of low-density SNP panels and imputation for cost-effective genomic selection in four aquaculture species

Genomic selection can accelerate genetic progress in aquaculture breeding programmes, particularly for traits measured on siblings of selection candidates. However, it is not widely implemented in most aquaculture species, and remains expensive due to high genotyping costs. Genotype imputation is a promising strategy that can reduce genotyping costs and facilitate the broader uptake of genomic selection in aquaculture breeding programmes. Genotype imputation can predict ungenotyped SNPs in populations genotyped at a low-density (LD), using a reference population genotyped at a high-density (HD). In this study, we used datasets of four aquaculture species (Atlantic salmon, turbot, common carp and Pacific oyster), phenotyped for different traits, to investigate the efficacy of genotype imputation for cost-effective genomic selection. The four datasets had been genotyped at HD, and eight LD panels (300-6,000 SNPs) were generated in silico. SNPs were selected to be: i) evenly distributed according to physical position ii) selected to minimise the linkage disequilibrium between adjacent SNPs or iii) randomly selected. Imputation was performed with three different software packages (AlphaImpute2, FImpute v.3 and findhap v.4). The results revealed that FImpute v.3 was faster and achieved higher imputation accuracies. Imputation accuracy increased with increasing panel density for both SNP selection methods, reaching correlations greater than 0.95 in the three fish species and 0.80 in Pacific oyster. In terms of genomic prediction accuracy, the LD and the imputed panels performed similarly, reaching values very close to the HD panels, except in the pacific oyster dataset, where the LD panel performed better than the imputed panel. In the fish species, when LD panels were used for genomic prediction without imputation, selection of markers based on either physical or genetic distance (instead of randomly) resulted in a high prediction accuracy, whereas imputation achieved near maximal prediction accuracy independently of the LD panel, showing higher reliability. Our results suggests that, in fish species, well-selected LD panels may achieve near maximal genomic selection prediction accuracy, and that the addition of imputation will result in maximal accuracy independently of the LD panel. These strategies represent effective and affordable methods to incorporate genomic selection into most aquaculture settings.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Inactivated ostreid herpesvirus-1 induces an innate immune response in the Pacific oyster, Crassostrea gigas, hemocytes

Infectious diseases are a major constraint to the expansion of shellfish production worldwide. Pacific oyster mortality syndrome (POMS), a polymicrobial disease triggered by the Ostreid herpesvirus-1 (OsHV-1), has devastated the global Pacific oyster (Crassostrea gigas) aquaculture industry. Recent ground-breaking research revealed that C. gigas possess an immune memory, capable of adaption, which improves the immune response upon a second exposure to a pathogen. This paradigm shift opens the door for developing 'vaccines' to improve shellfish survival during disease outbreaks. In the present study, we developed an in-vitro assay using hemocytes - the main effectors of the C. gigas immune system - collected from juvenile oysters susceptible to OsHV-1. The potency of multiple antigen preparations (e.g., chemically and physically inactivated OsHV-1, viral DNA, and protein extracts) to stimulate an immune response in hemocytes was evaluated using flow cytometry and droplet digital PCR to measure immune-related subcellular functions and gene expression, respectively. The immune response to the different antigens was benchmarked against that of hemocytes treated with Poly (I:C). We identified 10 antigen preparations capable of inducing immune stimulation in hemocytes (ROS production and positively expressed immune- related genes) after 1 h of exposure, without causing cytotoxicity. These findings are significant, as they evidence the potential for priming the innate immunity of oysters using viral antigens, which may enable cost-effective therapeutic treatment to mitigate OsHV-1/POMS. Further testing of these antigen preparations using an in-vivo infection model is essential to validate promising candidate pseudo-vaccines.

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.

ADAR-Editing during Ostreid Herpesvirus 1 Infection in Crassostrea gigas: Facts and Limitations

Ostreid herpesvirus-1 (OsHV-1) RNAs are enzymatically modified by A-to-I conversions during the infection of Crassostrea gigas. The increase of ADAR1 expression and hyper-editing activity parallel to OsHV-1 RNAs suggests a functional connection between dsRNA editing and antiviral responses. We analyzed 87 RNA-seq data sets from immuno-primed, resistant, and susceptible oysters exposed to OsHV-1 to compare the ADAR hyper-editing levels on host and viral transcripts and trace hyper-editing on the oyster genes. Host RNAs were more hyper-edited than viral RNAs, despite the increased editing of viral RNAs in late infection phases. A set of genes, representing ∼0.5% of the oyster transcriptome and including several tripartite motif-containing sequences, were constantly hyper-edited. Conversely, we identified genes involved in antiviral response, miRNA maturation, and epigenetic regulation that were hyper-edited in specific conditions only. Despite technical and biological bottlenecks that hamper the understanding of the bivalve “RNA editome,” available tools and technologies can be adapted to bivalve mollusks.

IMPORTANCE Ostreid herpesvirus-1 (OsHV-1) is a harmful pathogen of bivalve species, such as oysters. However, knowledge is lacking about host–virus interactions at the molecular level, hampering the possibility of a correct management of viral outbreaks and related massive mortalities. Notably, OsHV-1 transcripts are massively modified by host RNA editing enzyme during infection, resulting in multiple A-to-I variations along RNAs assuming double-strand conformations. The impact of these modifications on host transcripts is, however, not completely clear. Analyzing RNA-seq data of oysters infected with OsHV-1, we revealed that ∼0.5% of the oyster transcriptome is always enzymatically modified by ADAR, whereas genes involved in antiviral response, miRNA maturation, and epigenetic regulation were hyper-edited in specific conditions only. Despite our results, relevant technical bottlenecks impair an accurate quantification of RNA editing events, making necessary an approach specifically dedicated to the progressive understanding of oyster “RNA editome.”

Genomic Selection in Aquaculture Species

To date, genomic prediction has been conducted in about 20 aquaculture species, with a preference for intra-family genomic selection (GS). For every trait under GS, the increase in accuracy obtained by genomic estimated breeding values instead of classical pedigree-based estimation of breeding values is very important in aquaculture species ranging from 15% to 89% for growth traits, and from 0% to 567% for disease resistance. Although the implementation of GS in aquaculture is of little additional investment in breeding programs already implementing sib testing on pedigree, the deployment of GS remains sparse, but could be boosted by adaptation of cost-effective imputation from low-density panels. Moreover, GS could help to anticipate the effect of climate change by improving sustainability-related traits such as production yield (e.g., carcass or fillet yields), feed efficiency or disease resistance, and by improving resistance to environmental variation (tolerance to temperature or salinity variation). This chapter synthesized the literature in applications of GS in finfish, crustaceans and molluscs aquaculture in the present and future breeding programs.

Genomic Prediction for Whole Weight, Body Shape, Meat Yield, and Color Traits in the Portuguese Oyster Crassostrea angulata

Genetic improvement for quality traits, especially color and meat yield, has been limited in aquaculture because the assessment of these traits requires that the animals be slaughtered first. Genotyping technologies do, however, provide an opportunity to improve the selection efficiency for these traits. The main purpose of this study is to assess the potential for using genomic information to improve meat yield (soft tissue weight and condition index), body shape (cup and fan ratios), color (shell and mantle), and whole weight traits at harvest in the Portuguese oyster, Crassostrea angulata. The study consisted of 647 oysters: 188 oysters from 57 full-sib families from the first generation and 459 oysters from 33 full-sib families from the second generation. The number per family ranged from two to eight oysters for the first and 12–15 oysters for the second generation. After quality control, a set of 13,048 markers were analyzed to estimate the genetic parameters (heritability and genetic correlation) and predictive accuracy of the genomic selection for these traits. The multi-locus mixed model analysis indicated high estimates of heritability for meat yield traits: 0.43 for soft tissue weight and 0.77 for condition index. The estimated genomic heritabilities were 0.45 for whole weight, 0.24 for cup ratio, and 0.33 for fan ratio and ranged from 0.14 to 0.54 for color traits. The genetic correlations among whole weight, meat yield, and body shape traits were favorably positive, suggesting that the selection for whole weight would have beneficial effects on meat yield and body shape traits. Of paramount importance is the fact that the genomic prediction showed moderate to high accuracy for the traits studied (0.38–0.92). Therefore, there are good prospects to improve whole weight, meat yield, body shape, and color traits using genomic information. A multi-trait selection program using the genomic information can boost the genetic gain and minimize inbreeding in the long-term for this population.

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'.

Susceptibility variation to the main pathogens of Crassostrea gigas at the larval, spat and juvenile stages using unselected and selected oysters to OsHV-1 and/or V. aestuarianus

French commercial hatcheries are massively producing Crassostrea gigas selected for their higher resistance to OsHV-1, and soon should also implement selection for increasing resistance to Vibrio aestuarianus. The first objective of this study was to optimize the breeding programs for dual resistance to OsHV-1 and V. aestuarianus to determine the earliest life stage for which oysters are able to develop disease resistance. Wild stocks and selected families were tested using experimental infections by both pathogens at the larval, spat and juvenile stages. Oyster families could be evaluated for OsHV-1 as soon as the larval stage by a bath method, but this only highlighted the most resistant families; those that showed the highest resistance to V. aestuarianus could be determined using the cohabitation method at the juvenile stage. The second objective of this study was to determine if selection to increase/decrease the resistance to OsHV-1 and V. aestuarianus could have an impact on other major pathogens currently detected in hatchery at the larval stage, and in nursery and field at the spat/juveniles stages (V. coralliilyticus, V. crassostreae, V. tasmaniensis, V. neptunius, V. europaeus, V. harveyi, V. chagasi). No relationship was found between mortality caused by V. aestuarianus/OsHV-1 and the mortality caused by the other virulent bacterial strains tested regardless the stages, except between OsHV-1 and V. tasmaniensis at the juvenile stage. Finally, miscellaneous findings were evidenced such as (1) bath for bacterial challenges was not adapted for spat, (2) the main pathogens at the larval stage were OsHV-1 and V. coralliilyticus using bath, while it was V. coralliilyticus, V. europaeus, and V. neptunius at the juvenile stage by injection, and (4) variation in mortality was observed among families/wild controls for all pathogens at larval and juvenile stages, except for V. harveyi for larvae.

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.

Prediction Accuracies of Genomic Selection for Nine Commercially Important Traits in the Portuguese Oyster (Crassostrea angulata) Using DArT-Seq Technology

Genomic selection has been widely used in terrestrial animals but has had limited application in aquaculture due to relatively high genotyping costs. Genomic information has an important role in improving the prediction accuracy of breeding values, especially for traits that are difficult or expensive to measure. The purposes of this study were to (i) further evaluate the use of genomic information to improve prediction accuracies of breeding values from, (ii) compare different prediction methods (BayesA, BayesCπ and GBLUP) on prediction accuracies in our field data, and (iii) investigate the effects of different SNP marker densities on prediction accuracies of traits in the Portuguese oyster (Crassostrea angulata). The traits studied are all of economic importance and included morphometric traits (shell length, shell width, shell depth, shell weight), edibility traits (tenderness, taste, moisture content), and disease traits (Polydora sp. and Marteilioides chungmuensis). A total of 18,849 single nucleotide polymorphisms were obtained from genotyping by sequencing and used to estimate genetic parameters (heritability and genetic correlation) and the prediction accuracy of genomic selection for these traits. Multi-locus mixed model analysis indicated high estimates of heritability for edibility traits; 0.44 for moisture content, 0.59 for taste, and 0.72 for tenderness. The morphometric traits, shell length, shell width, shell depth and shell weight had estimated genomic heritabilities ranging from 0.28 to 0.55. The genomic heritabilities were relatively low for the disease related traits: Polydora sp. prevalence (0.11) and M. chungmuensis (0.10). Genomic correlations between whole weight and other morphometric traits were from moderate to high and positive (0.58–0.90). However, unfavourably positive genomic correlations were observed between whole weight and the disease traits (0.35–0.37). The genomic best linear unbiased prediction method (GBLUP) showed slightly higher accuracy for the traits studied (0.240–0.794) compared with both BayesA and BayesCπ methods but these differences were not significant. In addition, there is a large potential for using low-density SNP markers for genomic selection in this population at a number of 3000 SNPs. Therefore, there is the prospect to improve morphometric, edibility and disease related traits using genomic information in this species.

Differential Mortality and High Viral Load in Naive Pacific Oyster Families Exposed to OsHV-1 Suggests Tolerance Rather than Resistance to Infection

Pacific oysters, Crassostrea gigas, are one of the most productive aquaculture species in the world. However, they are threatened by the spread of Ostreid herpesvirus-1 (OsHV-1) and its microvariants (collectively "µvars"), which cause mass mortalities in all life stages of Pacific oysters globally. Breeding programs have been successful in reducing mortality due to OsHV-1 variants following viral outbreaks; however, an OsHV-1-resistant oyster line does not yet exist in the United States (US), and it is unknown how OsHV-1 µvars will affect US oyster populations compared to the current variant, which is similar to the OsHV-1 reference, found in Tomales Bay, CA. The goals of this study were to investigate the resistance of C. gigas juveniles produced by the Molluscan Broodstock Program (MBP) to three variants of OsHV-1: a California reference OsHV-1, an Australian µvar, and a French µvar. This is the first study to directly compare OsHV-1 µvars to a non-µvar. The survival probability of oysters exposed to the French (FRA) or Australian (AUS) µvar was significantly lower (43% and 71%, respectively) than to the reference variant and controls (96%). No oyster family demonstrated resistance to all three OsHV-1 variants, and many surviving oysters contained high copy numbers of viral DNA (mean ~3.53 × 108). These results indicate that the introduction of OsHV-1 µvars could have substantial effects on US Pacific oyster aquaculture if truly resistant lines are not achieved, and highlight the need to consider resistance to infection in addition to survival as traits in breeding programs to reduce the risk of the spread of OsHV-1 variants.

Functional Insights From the Evolutionary Diversification of Big Defensins

Big defensins are antimicrobial polypeptides believed to be the ancestors of β-defensins, the most evolutionary conserved family of host defense peptides (HDPs) in vertebrates. Nevertheless, big defensins underwent several independent gene loss events during animal evolution, being only retained in a limited number of phylogenetically distant invertebrates. Here, we explore the evolutionary history of this fascinating HDP family and investigate its patchy distribution in extant metazoans. We highlight the presence of big defensins in various classes of lophotrochozoans, as well as in a few arthropods and basal chordates (amphioxus), mostly adapted to life in marine environments. Bivalve mollusks often display an expanded repertoire of big defensin sequences, which appear to be the product of independent lineage-specific gene tandem duplications, followed by a rapid molecular diversification of newly acquired gene copies. This ongoing evolutionary process could underpin the simultaneous presence of canonical big defensins and non-canonical (β-defensin-like) sequences in some species. The big defensin genes of mussels and oysters, two species target of in-depth studies, are subjected to gene presence/absence variation (PAV), i.e., they can be present or absent in the genomes of different individuals. Moreover, big defensins follow different patterns of gene expression within a given species and respond differently to microbial challenges, suggesting functional divergence. Consistently, current structural data show that big defensin sequence diversity affects the 3D structure and biophysical properties of these polypeptides. We discuss here the role of the N-terminal hydrophobic domain, lost during evolution toward β-defensins, in the big defensin stability to high salt concentrations and its mechanism of action. Finally, we discuss the potential of big defensins as markers for animal health and for the nature-based design of novel therapeutics active at high salt concentrations.

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.

Genomic Prediction Using Low Density Marker Panels in Aquaculture: Performance Across Species, Traits, and Genotyping Platforms

Genomic selection increases the rate of genetic gain in breeding programs, which results in significant cumulative improvements in commercially important traits such as disease resistance. Genomic selection currently relies on collecting genome-wide genotype data accross a large number of individuals, which requires substantial economic investment. However, global aquaculture production predominantly occurs in small and medium sized enterprises for whom this technology can be prohibitively expensive. For genomic selection to benefit these aquaculture sectors, more cost-efficient genotyping is necessary. In this study the utility of low and medium density SNP panels (ranging from 100 to 9,000 SNPs) to accurately predict breeding values was tested and compared in four aquaculture datasets with different characteristics (species, genome size, genotyping platform, family number and size, total population size, and target trait). The traits show heritabilities between 0.19–0.49, and genomic prediction accuracies using the full density panel of 0.55–0.87. A consistent pattern of genomic prediction accuracy was observed across species with little or no accuracy reduction until SNP density was reduced below 1,000 SNPs (prediction accuracies of 0.44–0.75). Below this SNP density, heritability estimates and genomic prediction accuracies tended to be lower and more variable (93% of maximum accuracy achieved with 1,000 SNPs, 89% with 500 SNPs, and 70% with 100 SNPs). A notable drop in accuracy was observed between 200 SNP panels (0.44–0.75) and 100 SNP panels (0.39–0.66). Now that a multitude of studies have highlighted the benefits of genomic over pedigree-based prediction of breeding values in aquaculture species, the results of the current study highlight that these benefits can be achieved at lower SNP densities and at lower cost, raising the possibility of a broader application of genetic improvement in smaller and more fragmented aquaculture settings.