Genomic Predictions and Genome-Wide Association Study of Resistance Against Piscirickettsia salmonis in Coho Salmon (Oncorhynchus kisutch) Using ddRAD Sequencing

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.

GWAS on Imputed Whole-Genome Sequence Variants Reveal Genes Associated with Resistance to Piscirickettsia salmonis in Rainbow Trout (Oncorhynchus mykiss)

Genome-wide association studies (GWAS) allow the identification of associations between genetic variants and important phenotypes in domestic animals, including disease-resistance traits. Whole Genome Sequencing (WGS) data can help increase the resolution and statistical power of association mapping. Here, we conduced GWAS to asses he facultative intracellular bacterium Piscirickettsia salmonis, which affects farmed rainbow trout, Oncorhynchus mykiss, in Chile using imputed genotypes at the sequence level and searched for candidate genes located in genomic regions associated with the trait. A total of 2130 rainbow trout were intraperitoneally challenged with P. salmonis under controlled conditions and genotyped using a 57K single nucleotide polymorphism (SNP) panel. Genotype imputation was performed in all the genotyped animals using WGS data from 102 individuals. A total of 488,979 imputed WGS variants were available in the 2130 individuals after quality control. GWAS revealed genome-wide significant quantitative trait loci (QTL) in Omy02, Omy03, Omy25, Omy26 and Omy27 for time to death and in Omy26 for binary survival. Twenty-four (24) candidate genes associated with P. salmonis resistance were identified, which were mainly related to phagocytosis, innate immune response, inflammation, oxidative response, lipid metabolism and apoptotic process. Our results provide further knowledge on the genetic variants and genes associated with resistance to intracellular bacterial infection in rainbow trout.

Comparing mixed models and Random Forest association tests using naturalGWAS and a Striped Bass SNP dataset

In this study, we used the phenotype simulation package naturalGWAS to test the performance of Zhao's Random Forest method in comparison to an uncorrected Random Forest test, latent factor mixed models (LFMM), genome-wide efficient mixed models (GEMMA), and confounder adjusted linear regression (CATE). We created 400 sets of phenotypes, corresponding to five effect sizes and 2, 5, 15, or 30 causal loci, simulated from two empirical datasets containing SNPs from Striped Bass representing three and 13 populations. All association methods were evaluated for their ability to detect genotype-phenotype associations based on power, false discovery rates, and number of false positives. Genomic inflation was highest for uncorrected Random Forest and LFMM tests and lowest for Gemma and Zhao's Random Forest. All association tests had similar power to detect causal loci, and Zhao's Random Forest had the lowest false discovery rate in all scenarios. To measure the performance of association tests in small datasets with few loci surrounding a causal gene we also ran analyses again after removing causal loci from each dataset. All association tests were only able to find true positives, defined as loci located within 30k bp of a causal locus, in 3%-18% of simulations. In contrast, at least one false positive was found in 17%-44% of simulations. Zhao's Random Forest again identified the fewest false positives of all association tests studied. The ability to test the power of association tests for individual empirical datasets can be an extremely useful first step when designing a GWAS study.

Genotyping-by-Sequencing Reveals the Impact of Restocking on Wild Common Carp Populations of the Southern Caspian Basin

Understanding the population structure and level of genetic diversity of wild populations is fundamental for appropriate stock management and species conservation. The common carp (Cyprinus carpio) is one of the most important bony fish throughout the Southern coastline of the Caspian Sea, but captures of this species have seen a dramatic reduction during the last decade. As a consequence, a restocking program has been put in place to maintain C. carpio populations, but its impact is not clear. In the present study, the population structure and genetic diversity of C. carpio in the Southern Caspian basin was determined using 17,828 single-nucleotide polymorphism (SNP) markers. A total of 117 individuals collected from four different locations in the southern Caspian basin and a farm were genotyped by genotyping-by-sequencing. The overall Fst obtained was 0.04, indicating a low level of differentiation between populations, and most genetic diversity was attributed to within population variation (97%). The low Fst values suggest that frequent migration events between different locations occur, and three migration events were inferred in the present study. However, each population still showed a distinct genetic profile, which allowed distinguishing the origin of the fish. This indicates that the ongoing restocking program is maintaining the differences between populations to some extent. Nonetheless, high inbreeding and low heterozygosity were detected in all populations, suggesting that additional conservation efforts are required to protect C. carpio populations in the Southern coast of the Caspian Sea.

Why Does Piscirickettsia salmonis Break the Immunological Paradigm in Farmed Salmon? Biological Context to Understand the Relative Control of Piscirickettsiosis

Piscirickettsiosis (SRS) has been the most important infectious disease in Chilean salmon farming since the 1980s. It was one of the first to be described, and to date, it continues to be the main infectious cause of mortality. How can we better understand the epidemiological situation of SRS? The catch-all answer is that the Chilean salmon farming industry must fight year after year against a multifactorial disease, and apparently only the environment in Chile seems to favor the presence and persistence of Piscirickettsia salmonis. This is a fastidious, facultative intracellular bacterium that replicates in the host’s own immune cells and antigen-presenting cells and evades the adaptive cell-mediated immune response, which is why the existing vaccines are not effective in controlling it. Therefore, the Chilean salmon farming industry uses a lot of antibiotics—to control SRS—because otherwise, fish health and welfare would be significantly impaired, and a significantly higher volume of biomass would be lost per year. How can the ever-present risk of negative consequences of antibiotic use in salmon farming be balanced with the productive and economic viability of an animal production industry, as well as with the care of the aquatic environment and public health and with the sustainability of the industry? The answer that is easy, but no less true, is that we must know the enemy and how it interacts with its host. Much knowledge has been generated using this line of inquiry, however it remains insufficient. Considering the state-of-the-art summarized in this review, it can be stated that, from the point of view of fish immunology and vaccinology, we are quite far from reaching an effective and long-term solution for the control of SRS. For this reason, the aim of this critical review is to comprehensively discuss the current knowledge on the interaction between the bacteria and the host to promote the generation of more and better measures for the prevention and control of SRS.

Genome-wide analysis of acute low salinity tolerance in the eastern oyster Crassostrea virginica and potential of genomic selection for trait improvement

As the global demand for seafood increases, research into the genetic basis of traits that can increase aquaculture production is critical. The eastern oyster (Crassostrea virginica) is an important aquaculture species along the Atlantic and Gulf Coasts of the United States, but increases in heavy rainfall events expose oysters to acute low salinity conditions, which negatively impact production. Low salinity survival is known to be a moderately heritable trait, but the genetic architecture underlying this trait is still poorly understood. In this study, we used ddRAD sequencing to generate genome-wide single-nucleotide polymorphism (SNP) data for four F₂ families to investigate the genomic regions associated with survival in extreme low salinity (<3). SNP data were also used to assess the feasibility of genomic selection (GS) for improving this trait. Quantitative trait locus (QTL) mapping and combined linkage disequilibrium analysis revealed significant QTL on eastern oyster chromosomes 1 and 7 underlying both survival and day to death in a 36-day experimental challenge. Significant QTL were located in genes related to DNA/RNA function and repair, ion binding and membrane transport, and general response to stress. GS was investigated using Bayesian linear regression models and prediction accuracies ranged from 0.48 to 0.57. Genomic prediction accuracies were largest using the BayesB prior and prediction accuracies did not substantially decrease when SNPs located within the QTL region on Chr1 were removed, suggesting that this trait is controlled by many genes of small effect. Our results suggest that GS will likely be a viable option for improvement of survival in extreme low salinity.

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.

Accuracies of genomic predictions for disease resistance of striped catfish to Edwardsiella ictaluri using artificial intelligence algorithms

Assessments of genomic prediction accuracies using artificial intelligent (AI) algorithms (i.e., machine and deep learning methods) are currently not available or very limited in aquaculture species. The principal aim of this study was to examine the predictive performance of these new methods for disease resistance to Edwardsiella ictaluri in a population of striped catfish Pangasianodon hypophthalmus and to make comparisons with four common methods, i.e., pedigree-based best linear unbiased prediction (PBLUP), genomic-based best linear unbiased prediction (GBLUP), single-step GBLUP (ssGBLUP) and a nonlinear Bayesian approach (notably BayesR). Our analyses using machine learning (i.e., ML-KAML) and deep learning (i.e., DL-MLP and DL-CNN) together with the four common methods (PBLUP, GBLUP, ssGBLUP, and BayesR) were conducted for two main disease resistance traits (i.e., survival status coded as 0 and 1 and survival time, i.e., days that the animals were still alive after the challenge test) in a pedigree consisting of 560 individual animals (490 offspring and 70 parents) genotyped for 14,154 single nucleotide polymorphism (SNPs). The results using 6,470 SNPs after quality control showed that machine learning methods outperformed PBLUP, GBLUP, and ssGBLUP, with the increases in the prediction accuracies for both traits by 9.1–15.4%. However, the prediction accuracies obtained from machine learning methods were comparable to those estimated using BayesR. Imputation of missing genotypes using AlphaFamImpute increased the prediction accuracies by 5.3–19.2% in all the methods and data used. On the other hand, there were insignificant decreases (0.3–5.6%) in the prediction accuracies for both survival status and survival time when multivariate models were used in comparison to univariate analyses. Interestingly, the genomic prediction accuracies based on only highly significant SNPs (P < 0.00001, 318–400 SNPs for survival status and 1,362–1,589 SNPs for survival time) were somewhat lower (0.3–15.6%) than those obtained from the whole set of 6,470 SNPs. In most of our analyses, the accuracies of genomic prediction were somewhat higher for survival time than survival status (0/1 data). It is concluded that although there are prospects for the application of genomic selection to increase disease resistance to E. ictaluri in striped catfish breeding programs, further evaluation of these methods should be made in independent families/populations when more data are accumulated in future generations to avoid possible biases in the genetic parameters estimates and prediction accuracies for the disease-resistant traits studied in this population of striped catfish P. hypophthalmus.

Genomic Selection and Genome-wide Association Study for Feed-Efficiency Traits in a Farmed Nile Tilapia (Oreochromis niloticus) Population

Nile tilapia is a key aquaculture species with one of the highest production volumes globally. Genetic improvement of feed efficiency via selective breeding is an important goal, and genomic selection may expedite this process. The aims of this study were to 1) dissect the genetic architecture of feed-efficiency traits in a Nile tilapia breeding population, 2) map the genomic regions associated with these traits and identify candidate genes, 3) evaluate the accuracy of breeding value prediction using genomic data, and 4) assess the impact of the genetic marker density on genomic prediction accuracies. Using an experimental video recording trial, feed conversion ratio (FCR), body weight gain (BWG), residual feed intake (RFI) and feed intake (FI) traits were recorded in 40 full-sibling families from the GIFT (Genetically Improved Farmed Tilapia) Nile tilapia breeding population. Fish were genotyped with a ThermoFisher Axiom 65 K Nile tilapia SNP array. Significant heritabilities, ranging from 0.12 to 0.22, were estimated for all the assessed traits using the genomic relationship matrix. A negative but favourable genetic correlation was found between BWG and the feed-efficiency related traits; -0.60 and -0.63 for FCR and RFI, respectively. While the genome-wide association analyses suggested a polygenic genetic architecture for all the measured traits, there were significant QTL identified for BWG and FI on chromosomes seven and five respectively. Candidate genes previously found to be associated with feed-efficiency traits were located in these QTL regions, including ntrk3a, ghrh and eif4e3. The accuracy of breeding value prediction using the genomic data was up to 34% higher than using pedigree records. A SNP density of approximately 5,000 SNPs was sufficient to achieve similar prediction accuracy as the full genotype data set. Our results highlight the potential of genomic selection to improve feed efficiency traits in Nile tilapia breeding programmes.

Impact of genotypic errors with equal and unequal family contribution on accuracy of genomic prediction in aquaculture using simulation

Genotypic errors, conflict between recorded genotype and the true genotype, can lead to false or biased population genetic parameters. Here, the effect of genotypic errors on accuracy of genomic predictions and genomic relationship matrix are investigated using a simulation study based on population and genomic structure comparable to black tiger prawn, Penaeus monodon. Fifty full-sib families across five generations with phenotypic and genotypic information on 53 K SNPs were simulated. Ten replicates of different scenarios with three heritability estimates, equal and unequal family contributions were generated. Within each scenario, four SNP densities and three genotypic error rates in each SNP density were implemented. Results showed that family contribution did not have a substantial impact on accuracy of predictions across different datasets. In the absence of genotypic errors, 3 K SNP density was found to be efficient in estimating the accuracy, whilst increasing the SNP density from 3 to 20 K resulted in a marginal increase in accuracy of genomic predictions using the current population and genomic parameters. In addition, results showed that the presence of even 10% errors in a 10 and 20 K SNP panel might not have a severe impact on accuracy of predictions. However, below 10 K marker density, even a 5% error can result in lower accuracy of predictions.

Temporal Gene Expression Signature of Plasma Extracellular Vesicles-MicroRNAs from Post-Smolt Coho Salmon Challenged with Piscirickettsia salmonis

Piscirickettsiosis is the most important bacterial disease in the Chilean salmon industry, which has borne major economic losses due to failure to control it. Cells use extracellular vesicles (EVs) as an inter-cellular communicators to deliver several factors (e.g., microRNAs) that may regulate the responses of other cells. However, there is limited knowledge about the identification and characterization of EV-miRNAs in salmonids or the effect of infections on these. In this study, Illumina sequencing technology was used to identify Coho salmon plasma EV-miRNAs upon Piscirickettsia salmonis infection at four different time points. A total of 118 novels and 188 known EV-miRNAs, including key immune teleost miRNAs families (e.g., miR-146, miR-122), were identified. A total of 245 EV-miRNAs were detected as differentially expressed (FDR < 5%) in terms of control, with a clear down-regulation pattern throughout the disease. KEGG enrichment results of EV-miRNAs target genes showed that they were grouped mainly in cellular, stress, inflammation and immune responses. Therefore, it is hypothesized that P. salmonis could potentially benefit from unbalanced modulation response of Coho salmon EV-miRNAs in order to promote a hyper-inflammatory and compromised immune response through the suppression of different key immune host miRNAs during the course of the infection, as indicated by the results of this study.

The Phagosome-Lysosome Fusion Is the Target of a Purified Quillaja saponin Extract (PQSE) in Reducing Infection of Fish Macrophages by the Bacterial Pathogen Piscirickettsia salmonis

Piscirickettsia salmonis, the etiological agent of Piscirickettsiosis, is a Gram-negative and facultative intracellular pathogen that has affected the Chilean salmon industry since 1989. The bacterium is highly aggressive and can survive and replicate within fish macrophages using the Dot/Icm secretion system to evade the host’s immune response and spread systemically. To date, no efficient control measures have been developed for this disease; therefore, the producers use large amounts of antibiotics to control this pathogen. In this frame, this work has focused on evaluating the use of saponins from Quillaja saponaria as a new alternative to control the Piscirickettsiosis. It has been previously reported that purified extract of Q. saponaria (PQSE) displays both antimicrobial activity against pathogenic bacteria and viruses and adjuvant properties. Our results show that PQSE does not present antimicrobial activity against P. salmonis, although it reduces P. salmonis infection in an in vitro model, promoting the phagosome–lysosome fusion. Additionally, we demonstrate that PQSE modulates the expression of IL-12 and IL-10 in infected cells, promoting the immune response against the pathogen and reducing the expression of pathogen virulence genes. These results together strongly argue for specific anti-invasion and anti-intracellular replication effects induced by the PQSE in macrophages.

Optimizing genomic prediction using low-density marker panels for streptococcosis resistance in red tilapia (Oreochromis spp.)

Streptococcosis is a major disease that causes huge economic losses to tilapia farming in Thailand. Breeding for Streptococcosis agalactiae resistant strains using the pedigree BLUP method has proven an effective approach to control the disease in red tilapia, but the accuracy of selection is relatively low. Genomic selection, which is based on genome-wide markers to predict genomic breeding values of selection candidates, provides a powerful approach for accelerating genetic progress and producing permanent gains in the population. We evaluated the implementation of four genomic prediction models, GBLUP, ssGBLUP, BayesB and BayesC, using 19 sets of SNP markers (ranging from 500 to 24 582 SNPs) in 886 fish challenged with S. agalactiae. The accuracy of prediction was estimated using a five-fold cross-validation approach, with 708 and 178 individuals sampled for the training and validation sets respectively. Prediction of the accuracy of each of the models was improved substantially compared with PBLUP (10%) using 1000 informative SNPs. The GBLUP model (65%), which required less computing time, outperformed the remaining models - ssGBLUP (53%), BayesB (47%) and BayesC (42%).

A genomic‐based vision on the genetic diversity and key performance traits in selectively bred Arctic charr ( Salvelinus alpinus )

Routine implementation of genomic information for guiding selection decisions is not yet common in the majority of aquaculture species. Reduced representation sequencing approaches offer a cost‐effective solution for obtaining genome‐wide information in species with a limited availability of genomic resources. In the current study, we implemented double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) on an Arctic charr strain with the longest known history of selection (approximately 40 years) aiming to improve selection decisions. In total, 1730 animals reared at four different farms in Sweden and spanning from year classes 2013–2017 were genotyped using ddRAD‐seq. Approximately 5000 single nucleotide polymorphisms (SNPs) were identified, genetic diversity‐related metrics were estimated, and genome‐wide association studies (GWAS) for body length at different time points and age of sexual maturation were conducted. Low genetic differentiation amongst animals from the different farms was observed based on both the results from pairwise F_st values and principal component analysis (PCA). The existence of associations was investigated between the mean genome‐wide heterozygosity of each full‐sib family (year class 2017) and the corresponding inbreeding coefficient or survival to the eyed stage. A moderate correlation (−0.33) was estimated between the mean observed heterozygosity of each full‐sib family and the corresponding inbreeding coefficient, while no linear association was obtained with the survival to the eyed stage. GWAS did not detect loci with major effect for any of the studied traits. However, genomic regions explaining more than 1% of the additive genetic variance for either studied traits were suggested across 14 different chromosomes. Overall, key insights valuable for future selection decisions of Arctic charr have been obtained, suggesting ddRAD as an attractive genotyping platform for obtaining genome‐wide information in a cost‐effective manner.

Genome-wide association study of cyanogenic glycosides, proline, sugars, and pigments in Eucalyptus cladocalyx after 18 consecutive dry summers

Natural variation of cyanogenic glycosides, soluble sugars, proline, and nondestructive optical sensing of pigments (chlorophyll, flavonols, and anthocyanins) was examined in ex situ natural populations of Eucalyptus cladocalyx F. Muell. grown under dry environmental conditions in the southern Atacama Desert, Chile. After 18 consecutive dry seasons, considerable plant-to-plant phenotypic variation for all the traits was observed in the field. For example, leaf hydrogen cyanide (HCN) concentrations varied from 0 (two acyanogenic individuals) to 1.54 mg cyanide g^-1 DW. Subsequent genome-wide association study revealed associations with several genes with a known function in plants. HCN content was associated robustly with genes encoding Cytochrome P450 proteins, and with genes involved in the detoxification mechanism of HCN in cells (β-cyanoalanine synthase and cyanoalanine nitrilase). Another important finding was that sugars, proline, and pigment content were linked to genes involved in transport, biosynthesis, and/or catabolism. Estimates of genomic heritability (based on haplotypes) ranged between 0.46 and 0.84 (HCN and proline content, respectively). Proline and soluble sugars had the highest predictive ability of genomic prediction models (PA = 0.65 and PA = 0.71, respectively). PA values for HCN content and flavonols were relatively moderate, with estimates ranging from 0.44 to 0.50. These findings provide new understanding on the genetic architecture of cyanogenic capacity, and other key complex traits in cyanogenic E. cladocalyx.

Increased accuracy of genomic predictions for growth under chronic thermal stress in rainbow trout by prioritizing variants from GWAS using imputed sequence data

Through imputation of genotypes, genome‐wide association study (GWAS) and genomic prediction (GP) using whole‐genome sequencing (WGS) data are cost‐efficient and feasible in aquaculture breeding schemes. The objective was to dissect the genetic architecture of growth traits under chronic heat stress in rainbow trout (Oncorhynchus mykiss) and to assess the accuracy of GP based on imputed WGS and different preselected single nucleotide polymorphism (SNP) arrays. A total of 192 and 764 fish challenged to a heat stress experiment for 62 days were genotyped using a customized 1 K and 26 K SNP panels, respectively, and then, genotype imputation was performed from a low‐density chip to WGS using 102 parents (36 males and 66 females) as the reference population. Imputed WGS data were used to perform GWAS and test GP accuracy under different preselected SNP scenarios. Heritability was estimated for body weight (BW), body length (BL) and average daily gain (ADG). Estimates using imputed WGS data ranged from 0.33 ± 0.05 to 0.55 ± 0.05 for growth traits under chronic heat stress. GWAS revealed that the top five cumulatively SNPs explained a maximum of 0.94%, 0.86% and 0.51% of genetic variance for BW, BL and ADG, respectively. Some important functional candidate genes associated with growth‐related traits were found among the most important SNPs, including signal transducer and activator of transcription 5B and 3 (STAT5B and STAT3, respectively) and cytokine‐inducible SH2‐containing protein (CISH). WGS data resulted in a slight increase in prediction accuracy compared with pedigree‐based method, whereas preselected SNPs based on the top GWAS hits improved prediction accuracies, with values ranging from 1.2 to 13.3%. Our results support the evidence of the polygenic nature of growth traits when measured under heat stress. The accuracies of GP can be improved using preselected variants from GWAS, and the use of WGS marginally increases prediction accuracy.

Exploring genetic resistance to infectious salmon anaemia virus in Atlantic salmon by genome-wide association and RNA sequencing

Background

Infectious Salmonid Anaemia Virus (ISAV) causes a notifiable disease that poses a large threat for Atlantic salmon (Salmo salar) aquaculture worldwide. There is no fully effective treatment or vaccine, and therefore selective breeding to increase resistance to ISAV is a promising avenue for disease prevention. Genomic selection and potentially genome editing can be applied to enhance host resistance, and these approaches benefit from improved knowledge of the genetic and functional basis of the target trait. The aim of this study was to characterise the genetic architecture of resistance to ISAV in a commercial Atlantic salmon population and study its underlying functional genomic basis using RNA Sequencing.

Results

A total of 2833 Atlantic salmon parr belonging to 194 families were exposed to ISAV in a cohabitation challenge in which cumulative mortality reached 63% over 55 days. A total of 1353 animals were genotyped using a 55 K SNP array, and the estimate of heritability for the trait of binary survival was 0.13–0.33 (pedigree-genomic). A genome-wide association analysis confirmed that resistance to ISAV was a polygenic trait, albeit a genomic region in chromosome Ssa13 was significantly associated with resistance and explained 3% of the genetic variance. RNA sequencing of the heart of 16 infected (7 and 14 days post infection) and 8 control fish highlighted 4927 and 2437 differentially expressed genes at 7 and 14 days post infection respectively. The complement and coagulation pathway was down-regulated in infected fish, while several metabolic pathways were up-regulated. The interferon pathway showed little evidence of up-regulation at 7 days post infection but was mildly activated at 14 days, suggesting a potential crosstalk between host and virus. Comparison of the transcriptomic response of fish with high and low breeding values for resistance highlighted TRIM25 as being up-regulated in resistant fish.

Conclusions

ISAV resistance shows moderate heritability with a polygenic architecture, but a significant QTL was detected on chromosome 13. A mild up-regulation of the interferon pathway characterises the response to the virus in heart samples from this population of Atlantic salmon, and candidate genes showing differential expression between samples with high and low breeding values for resistance were identified.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07671-6.

Detection of selection signatures in farmed coho salmon (Oncorhynchus kisutch) using dense genome-wide information

Animal domestication and artificial selection give rise to gradual changes at the genomic level in populations. Subsequent footprints of selection, known as selection signatures or selective sweeps, have been traced in the genomes of many animal livestock species by exploiting variation in linkage disequilibrium patterns and/or reduction of genetic diversity. Domestication of most aquatic species is recent in comparison with land animals, and salmonids are one of the most important fish species in aquaculture. Coho salmon (Oncorhynchus kisutch), cultivated primarily in Chile, has been subjected to breeding programs to improve growth, disease resistance traits, and flesh color. This study aimed to identify selection signatures that may be involved in adaptation to culture conditions and traits of productive interest. To do so, individuals of two domestic populations cultured in Chile were genotyped with 200 thousand SNPs, and analyses were conducted using iHS, XP-EHH and CLR. Several signatures of selection on different chromosomal regions were detected across both populations. Some of the identified regions under selection contained genes such anapc2, alad, chp2 and myn, which have been previously associated with body weight in Atlantic salmon, or sec24d and robo1, which have been associated with resistance to Piscirickettsia salmonis in coho salmon. Findings in our study can contribute to an integrated genome-wide map of selection signatures, to help identify the genetic mechanisms of phenotypic diversity in coho salmon.

Genetic Dissection of a Precocious Phenotype in Male Tiger Pufferfish (Takifugu rubripes) using Genotyping by Random Amplicon Sequencing, Direct (GRAS-Di)

The novel non-targeted PCR-based genotyping system, namely Genotyping by Random Amplicon Sequencing, Direct (GRAS-Di), is characterized by the simplicity in library construction and robustness against DNA degradation and is expected to facilitate advancements in genetics, in both basic and applied sciences. In this study, we tested the utility of GRAS-Di for genetic analysis in a cultured population of the tiger pufferfish Takifugu rubripes. The genetic analyses included family structure analysis, genetic map construction, and quantitative trait locus (QTL) analysis for the male precocious phenotype using a population consisting of four full-sib families derived from a genetically precocious line. An average of 4.7 million raw reads were obtained from 198 fish. Trimmed reads were mapped onto a Fugu reference genome for genotyping, and 21,938 putative single-nucleotide polymorphisms (SNPs) were obtained. These 22 K SNPs accurately resolved the sibship and parent-offspring pairs. A fine-scale linkage map (total size: 1,949 cM; average interval: 1.75 cM) was constructed from 1,423 effective SNPs, for which the allele inheritance patterns were known. QTL analysis detected a significant locus for testes weight on Chr_14 and three suggestive loci on Chr_1, Chr_8, and Chr_19. The significant QTL was shared by body length and body weight. The effect of each QTL was small (phenotypic variation explained, PVE: 3.1-5.9%), suggesting that the precociousness seen in the cultured pufferfish is polygenic. Taken together, these results indicate that GRAS-Di is a practical genotyping tool for aquaculture species and applicable for molecular breeding programs, such as marker-assisted selection and genomic selection.

Genotyping Strategies Using ddRAD Sequencing in Farmed Arctic Charr (Salvelinus alpinus)

Incorporation of genomic technologies into fish breeding programs is a modern reality, promising substantial advances regarding the accuracy of selection, monitoring the genetic diversity and pedigree record verification. Single nucleotide polymorphism (SNP) arrays are the most commonly used genomic tool, but the investments required make them unsustainable for emerging species, such as Arctic charr (Salvelinus alpinus), where production volume is low. The requirement to genotype a large number of animals for breeding practices necessitates cost effective genotyping approaches. In the current study, we used double digest restriction site-associated DNA (ddRAD) sequencing of either high or low coverage to genotype Arctic charr from the Swedish national breeding program and performed analytical procedures to assess their utility in a range of tasks. SNPs were identified and used for deciphering the genetic structure of the studied population, estimating genomic relationships and implementing an association study for growth-related traits. Missing information and underestimation of heterozygosity in the low coverage set were limiting factors in genetic diversity and genomic relationship analyses, where high coverage performed notably better. On the other hand, the high coverage dataset proved to be valuable when it comes to identifying loci that are associated with phenotypic traits of interest. In general, both genotyping strategies offer sustainable alternatives to hybridization-based genotyping platforms and show potential for applications in aquaculture selective breeding.

First genomic prediction and genome‐wide association for complex growth‐related traits in Rock Bream (Oplegnathus fasciatus)

Rock Bream (Oplegnathus fasciatus) is an important aquaculture species for offshore cage aquaculture and fish stocking of marine ranching in East Asia. Genomic selection has the potential to expedite genetic gain for the key target traits of a breeding program, but has not yet been evaluated in Oplegnathus. The purposes of the present study were to explore the performance of genomic selection to improve breeding value accuracy through real data analyses using six statistical models and to carry out genome‐wide association studies (GWAS) to dissect the genetic architecture of economically vital growth‐related traits (body weight, total length, and body depth) in the O. fasciatus population. After quality control, genotypes for 16,162 SNPs were acquired for 455 fish. Heritability was estimated to be moderate for the three traits (0.38 for BW, 0.33 for TL, and 0.24 for BD), and results of GWAS indicated that the underlying genetic architecture was polygenic. Six statistic models (GBLUP, BayesA, BayesB, BayesC, Bayesian Ridge‐Regression, and Bayesian LASSO) showed similar performance for the predictability of genomic estimated breeding value (GEBV). The low SNP density (around 1 K selected SNP based on GWAS) is sufficient for accurate prediction on the breeding value for the three growth‐related traits in the current studied population, which will provide a good compromise between genotyping costs and predictability in such standard breeding populations advanced. These consequences illustrate that the employment of genomic selection in O. fasciatus breeding could provide advantages for the selection of breeding candidates to facilitate complex economic growth traits.

Investigating mechanisms underlying genetic resistance to Salmon Rickettsial Syndrome in Atlantic salmon using RNA sequencing

Background

Salmon Rickettsial Syndrome (SRS), caused by Piscirickettsia salmonis, is one of the primary causes of morbidity and mortality in Atlantic salmon aquaculture, particularly in Chile. Host resistance is a heritable trait, and functional genomic studies have highlighted genes and pathways important in the response of salmon to the bacteria. However, the functional mechanisms underpinning genetic resistance are not yet well understood. In the current study, a large population of salmon pre-smolts were challenged with P. salmonis, with mortality levels recorded and samples taken for genotyping. In parallel, head kidney and liver samples were taken from animals of the same population with high and low genomic breeding values for resistance, and used for RNA-Sequencing to compare their transcriptome profile both pre and post infection.

Results

A significant and moderate heritability (h² = 0.43) was shown for the trait of binary survival. Genome-wide association analyses using 38 K imputed SNP genotypes across 2265 animals highlighted that resistance is a polygenic trait. Several thousand genes were identified as differentially expressed between controls and infected samples, and enriched pathways related to the host immune response were highlighted. In addition, several networks with significant correlation with SRS resistance breeding values were identified, suggesting their involvement in mediating genetic resistance. These included apoptosis, cytoskeletal organisation, and the inflammasome.

Conclusions

While resistance to SRS is a polygenic trait, this study has highlighted several relevant networks and genes that are likely to play a role in mediating genetic resistance. These genes may be future targets for functional studies, including genome editing, to further elucidate their role underpinning genetic variation in host resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07443-2.

Genome-Wide Marker Analysis for Traits of Economic Importance in Asian Seabass Lates calcarifer

To date, it is not known whether animal breeding values in Asian seabass (Lates calcarifer) can be estimated using single nucleotide polymorphisms (SNPs) generated from new high-throughput genotyping by sequencing platforms. The principal aim of the present study was to assess the genomic prediction accuracy for growth traits, survival, cannibalism, and disease resistance against Streptococcus iniae in this species L. calcarifer. Additionally, this study attempted to identify markers associated with the five traits studied as well as to understand if the genotype data can be used to estimate genetic parameters for these complex traits. The genomic best linear unbiased prediction (gBLUP) method was used to analyze 11,084 SNPs and showed that the prediction accuracies for growth traits (weight and length) were high (0.67–0.75). By contrast, these estimates for survival were low (0.25). Multi-locus mixed model analyses identified four SNPs significantly associated with body weight (p < 5 × 10−8 or −log10 p ≥ 5). There were, however, no significant associations detected for other traits. Similarly, the SNP heritability was moderate, while the estimates for other traits were approximated to zero and not significant. Genetic correlations between body weight and standard length were close to unity. Collectively, the results obtained from this study suggest that genotyping by sequencing platforms can provide informative DNA markers to conduct genome-wide association analysis, estimation of genetic parameters, and evaluation of genomic prediction accuracy for complex traits in Asian seabass.

Community Parameters and Genome-Wide RAD-Seq Loci of Ceratothoa oestroides Imply Its Transfer between Farmed European Sea Bass and Wild Farm-Aggregating Fish

Wild fish assemblages that aggregate within commercial marine aquaculture sites for feeding and shelter have been considered as a primary source of pathogenic parasites vectored to farmed fish maintained in net pens at an elevated density. In order to evaluate whether Ceratothoa oestroides (Isopoda, Cymothoidae), a generalist and pestilent isopod that is frequently found in Adriatic and Greek stocks of farmed European sea bass (Dicentrarchus labrax), transfers between wild and farmed fish, a RAD-Seq (restriction-site-associated DNA sequencing)-mediated genetic screening approach was employed. The double-digest RAD-Seq of 310 C. oestroides specimens collected from farmed European sea bass (138) and different wild farm-aggregating fish (172) identified 313 robust SNPs that evidenced a close genetic relatedness between the “wild” and “farmed” genotypes. ddRAD-Seq proved to be an effective method for detecting the discrete genetic structuring of C. oestroides and genotype intermixing between two populations. The parasite prevalence in the farmed sea bass was 1.02%, with a mean intensity of 2.0 and mean abundance of 0.02, while in the wild fish, the prevalence was 8.1%; the mean intensity, 1.81; and the mean abundance, 0.15. Such differences are likely a consequence of human interventions during the farmed fish’s rearing cycle that, nevertheless, did not affect the transfer of C. oestroides.

Temporal genome-wide DNA methylation signature of post-smolt Pacific salmon challenged with Piscirickettsia salmonis

Piscirickettsiosis is the most important bacterial disease in the Chilean salmon industry, which has sorted several efforts to its control, generating enormous economic losses. Epigenetic alterations, such as DNA methylation, can play a relevant role in the modulation of the metazoans response to pathogens. Bacterial disease may activate global and local immune responses generating intricate responses with significant biological impact in the host. However, it is scarcely understood how bacterial infections influence fish epigenetic alterations. In the present study, we utilized Pacific salmon and Piscirickettsiosis as model, to gain understanding into the dynamics of DNA methylation among fish-bacterial infection interactions. A genome-wide analysis of DNA methylation patterns in female spleen tissue of Pacific salmon was achieved by reduced representation bisulphite sequencing from a time course design. We determined 2,251, 1,918, and 2,516 differentially methylated regions DMRs among infected and control Pacific salmon in 1 dpi, 5 dpi, and 15 dpi, respectively. The mean methylation difference per DMR among control and infected groups was of ~35%, with an oscillatory pattern of hypo, hyper, and hypomethylation across the disease. DMCs, among the control and infected group, showed that they were statistically enriched in intergenic regions and depleted in exons. Functional annotation of the DMR genes demonstrated three KEGG principal categories, associated directly with the host response to pathogens infections. Our results provide the first evidence of epigenetic variation in fish provoked by bacterial infection and demonstrate that this variation can be modulated across the disease.

Development of a SNP linkage map and genome-wide association study for resistance to Aeromonas hydrophila in pacu (Piaractus mesopotamicus)

BACKGROUND

Pacu (Piaractus mesopotamicus) is one of the most important Neotropical aquaculture species from South America. Disease outbreaks caused by Aeromonas hydrophila infection have been considered significant contributors to the declining levels of pacu production. The current implementation of genomic selection for disease resistance has been adopted as a powerful strategy for improvement in fish species. This study aimed to investigate the genetic architecture of resistance to A. hydrophila in pacu via Genome-Wide Association Study (GWAS), the identification of suggestive Quantitative Trait Loci (QTLs) and putative genes associated with this trait. The genetic data were obtained from 381 juvenile individuals belonging to 14 full-sibling families. An experimental challenge was performed to gain access to the levels of genetic variation for resistance against the bacteria using the following trait definitions: binary test survival (TS) and time of death (TD).

RESULTS

The analyses of genetic parameters estimated moderate heritability (h2) for both resistance traits: 0.20 (± 0.09) for TS and 0.35 (± 0.15) for TD. A linkage map for pacu was developed to enable the GWAS, resulting in 27 linkage groups (LGs) with 17,453 mapped Single Nucleotide Polymorphisms (SNPs). The length of the LGs varied from 79.95 (LG14) to 137.01 (LG1) cM, with a total map length of 2755.60 cM. GWAS identified 22 putative QTLs associated to A. hydrophila resistance. They were distributed into 17 LGs, and were considered suggestive genomic regions explaining > 1% of the additive genetic variance (AGV) for the trait. Several candidate genes related to immune response were located close to the suggestive QTLs, such as tbk1, trim16, Il12rb2 and lyz2.

CONCLUSION

This study describes the development of the first medium density linkage map for pacu, which will be used as a framework to study relevant traits to the production of this species. In addition, the resistance to A. hydrophila was found to be moderately heritable but with a polygenic architecture suggesting that genomic selection, instead of marker assisted selection, might be useful for efficiently improving resistance to one of the most problematic diseases that affects the South American aquaculture.

Genomics to accelerate genetic improvement in tilapia

Selective breeding of tilapia populations started in the early 1990s and over the past three decades tilapia has become one of the most important farmed freshwater species, being produced in more than 125 countries around the globe. Although genome assemblies have been available since 2011, most of the tilapia industry still depends on classical selection techniques using mass spawning or pedigree information to select for growth traits with reported genetic gains of up to 20% per generation. The involvement of international breeding companies and research institutions has resulted in the rapid development and application of genomic resources in the last few years. GWAS and genomic selection are expected to contribute to uncovering the genetic variants involved in economically relevant traits and increasing the genetic gain in selective breeding programs, respectively. Developments over the next few years will probably focus on achieving a deep understanding of genetic architecture of complex traits, as well as accelerating genetic progress in the selection for growth-, quality- and robustness-related traits. Novel phenotyping technologies (i.e. phenomics), lower-cost whole-genome sequencing approaches, functional genomics and gene editing tools will be crucial in future developments for the improvement of tilapia aquaculture.

Combining Multiple Approaches and Models to Dissect the Genetic Architecture of Resistance to Infections in Fish

Infectious diseases represent a major threat for the sustainable development of fish farming. Efficient vaccines are not available against all diseases, and growing antibiotics resistance limits the use of antimicrobial drugs in aquaculture. It is therefore important to understand the basis of fish natural resistance to infections to help genetic selection and to develop new approaches against infectious diseases. However, the identification of the main mechanisms determining the resistance or susceptibility of a host to a pathogenic microbe is challenging, integrating the complexity of the variation of host genetics, the variability of pathogens, and their capacity of fast evolution and adaptation. Multiple approaches have been used for this purpose: (i) genetic approaches, QTL (quantitative trait loci) mapping or GWAS (genome-wide association study) analysis, to dissect the genetic architecture of disease resistance, and (ii) transcriptomics and functional assays to link the genetic constitution of a fish to the molecular mechanisms involved in its interactions with pathogens. To date, many studies in a wide range of fish species have investigated the genetic determinism of resistance to many diseases using QTL mapping or GWAS analyses. A few of these studies pointed mainly toward adaptive mechanisms of resistance/susceptibility to infections; others pointed toward innate or intrinsic mechanisms. However, in the majority of studies, underlying mechanisms remain unknown. By comparing gene expression profiles between resistant and susceptible genetic backgrounds, transcriptomics studies have contributed to build a framework of gene pathways determining fish responsiveness to a number of pathogens. Adding functional assays to expression and genetic approaches has led to a better understanding of resistance mechanisms in some cases. The development of knock-out approaches will complement these analyses and help to validate putative candidate genes critical for resistance to infections. In this review, we highlight fish isogenic lines as a unique biological material to unravel the complexity of host response to different pathogens. In the future, combining multiple approaches will lead to a better understanding of the dynamics of interaction between the pathogen and the host immune response, and contribute to the identification of potential targets of selection for improved resistance.

Application of Low Coverage Genotyping by Sequencing in Selectively Bred Arctic Charr (Salvelinus alpinus)

Arctic charr (Salvelinus alpinus) is a species of high economic value for the aquaculture industry, and of high ecological value due to its Holarctic distribution in both marine and freshwater environments. Novel genome sequencing approaches enable the study of population and quantitative genetic parameters even on species with limited or no prior genomic resources. Low coverage genotyping by sequencing (GBS) was applied in a selected strain of Arctic charr in Sweden originating from a landlocked freshwater population. For the needs of the current study, animals from year classes 2013 (171 animals, parental population) and 2017 (759 animals; 13 full sib families) were used as a template for identifying genome wide single nucleotide polymorphisms (SNPs). GBS libraries were constructed using the PstI and MspI restriction enzymes. Approximately 14.5K SNPs passed quality control and were used for estimating a genomic relationship matrix. Thereafter a wide range of analyses were conducted in order to gain insights regarding genetic diversity and investigate the efficiency of the genomic information for parentage assignment and breeding value estimation. Heterozygosity estimates for both year classes suggested a slight excess of heterozygotes. Furthermore, FST estimates among the families of year class 2017 ranged between 0.009 - 0.066. Principal components analysis (PCA) and discriminant analysis of principal components (DAPC) were applied aiming to identify the existence of genetic clusters among the studied population. Results obtained were in accordance with pedigree records allowing the identification of individual families. Additionally, DNA parentage verification was performed, with results in accordance with the pedigree records with the exception of a putative dam where full sib genotypes suggested a potential recording error. Breeding value estimation for juvenile growth through the usage of the estimated genomic relationship matrix clearly outperformed the pedigree equivalent in terms of prediction accuracy (0.51 opposed to 0.31). Overall, low coverage GBS has proven to be a cost-effective genotyping platform that is expected to boost the selection efficiency of the Arctic charr breeding program.

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.

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

In genomic selection (GS), genome-wide SNP markers are used to generate genomic estimated breeding values for selection candidates. The application of GS in shellfish looks promising and has the potential to help in dealing with one of the main issues currently affecting Pacific oyster production worldwide, which is the 'summer mortality syndrome'. This causes periodic mass mortality in farms worldwide and has mainly been attributed to a specific variant of the ostreid herpesvirus (OsHV-1). In the current study, we evaluated the potential of genomic selection for host resistance to OsHV-1 in Pacific oysters, and compared it with pedigree-based approaches. An OsHV-1 disease challenge was performed using an immersion-based virus exposure treatment for oysters for 7 days. A total of 768 samples were genotyped using the medium-density SNP array for oysters. A GWAS was performed for the survival trait using a GBLUP approach in blupf90 software. Heritability ranged from 0.25 ± 0.05 to 0.37 ± 0.05 (mean ± SE) based on pedigree and genomic information respectively. Genomic prediction was more accurate than pedigree prediction, and SNP density reduction had little impact on prediction accuracy until marker densities dropped below approximately 500 SNPs. This demonstrates the potential for GS in Pacific oyster breeding programmes, and importantly, demonstrates that a low number of SNPs might suffice to obtain accurate genomic estimated breeding values, thus potentially making the implementation of GS more cost effective.

Genetic Analysis of QTL for Resistance to Maize Lethal Necrosis in Multiple Mapping Populations

Maize lethal necrosis (MLN) occurs when maize chlorotic mottle virus (MCMV) and sugarcane mosaic virus (SCMV) co-infect maize plant. Yield loss of up to 100% can be experienced under severe infections. Identification and validation of genomic regions and their flanking markers can facilitate marker assisted breeding for resistance to MLN. To understand the status of previously identified quantitative trait loci (QTL)in diverse genetic background, F₃ progenies derived from seven bi-parental populations were genotyped using 500 selected kompetitive allele specific PCR (KASP) SNPs. The F₃ progenies were evaluated under artificial MLN inoculation for three seasons. Phenotypic analyses revealed significant variability (P ≤ 0.01) among genotypes for responses to MLN infections, with high heritability estimates (0.62 to 0.82) for MLN disease severity and AUDPC values. Linkage mapping and joint linkage association mapping revealed at least seven major QTL (qMLN3_130 and qMLN3_142, qMLN5_190 and qMLN5_202, qMLN6_85 and qMLN6_157 qMLN8_10 and qMLN9_142) spread across the 7-biparetal populations, for resistance to MLN infections and were consistent with those reported previously. The seven QTL appeared to be stable across genetic backgrounds and across environments. Therefore, these QTL could be useful for marker assisted breeding for resistance to MLN.

Utilizing field collected insects for next generation sequencing: Effects of sampling, storage, and DNA extraction methods

DNA sequencing technologies continue to advance the biological sciences, expanding opportunities for genomic studies of non-model organisms for basic and applied questions. Despite these opportunities, many next generation sequencing protocols have been developed assuming a substantial quantity of high molecular weight DNA (>100 ng), which can be difficult to obtain for many study systems. In particular, the ability to sequence field-collected specimens that exhibit varying levels of DNA degradation remains largely unexplored. In this study we investigate the influence of five traditional insect capture and curation methods on Double-Digest Restriction Enzyme Associated DNA (ddRAD) sequencing success for three wild bee species. We sequenced a total of 105 specimens (between 7-13 specimens per species and treatment). We additionally investigated how different DNA quality metrics (including pre-sequence concentration and contamination) predicted downstream sequencing success, and also compared two DNA extraction methods. We report successful library preparation for all specimens, with all treatments and extraction methods producing enough highly reliable loci for population genetic analyses. Although results varied between species, we found that specimens collected by net sampling directly into 100% EtOH, or by passive trapping followed by 100% EtOH storage before pinning tended to produce higher quality ddRAD assemblies, likely as a result of rapid specimen desiccation. Surprisingly, we found that specimens preserved in propylene glycol during field sampling exhibited lower-quality assemblies. We provide recommendations for each treatment, extraction method, and DNA quality assessment, and further encourage researchers to consider utilizing a wider variety of specimens for genomic analyses.

Genomic prediction for disease resistance to Hepatopancreatic parvovirus and growth, carcass and quality traits in Banana shrimp Fenneropenaeus merguiensis

Conventional genetic improvement of disease resistance in aquatic animal species involves challenge tests or using qPCR to quantify viral load that is costly, time-consuming and causing biosecurity concerns. Recent developments in high throughput next generation genome sequencing platforms such as genotyping by sequencing (GBS) have opened new possibilities for improving disease traits based on DNA information. The principal aim of this study was thus to examine potential application of genomic selection to improve resistance to hepatopancreatic parvovirus (HPV) in banana shrimp Fenneropenaeus merguiensis. Specifically, we used a total of 9472 single nucleotide polymorphisms (SNPs) developed de novo from GBS platforms to assess accuracy of genomic prediction for HPV resistance and growth, carcass and quality-related traits in this white shrimp species. Our multi-locus mixed model analysis showed moderate heritabilities for HPV resistance (h² = 0.46) and other traits studied (0.10 to 0.55). Genetic correlations of HPV titre with growth and carcass traits, estimated using SNPs markers, were negative (i.e., favourable), suggesting that selection for improved growth and carcass traits may have increased HPV resistance (i.e., reduced HPV titre). More importantly, our gBLUP model demonstrated that the accuracy of gBLUP prediction was moderate for HPV disease resistance (0.46). The genomic prediction accuracy was somewhat greater for growth and carcass related traits especially for body weight (0.76) and meat or tail weight (0.77). On the other hand, the prediction accuracy was from 0.25 to 0.41 for quality traits (raw and cooked colour and flesh streaks). Collectively, it is concluded that there are prospects to apply genomic selection in the genetic improvement for increased disease resistance, carcass and quality-related traits in this population of banana shrimp F. merguiensis.

Single-Step Genome-Wide Association Study for Resistance to Piscirickettsia salmonis in Rainbow Trout (Oncorhynchus mykiss)

One of the main pathogens affecting rainbow trout (Oncorhynchus mykiss) farming is the facultative intracellular bacteria Piscirickettsia salmonis Current treatments, such as antibiotics and vaccines, have not had the expected effectiveness in field conditions. Genetic improvement by means of selection for resistance is proposed as a viable alternative for control. Genomic information can be used to identify the genomic regions associated with resistance and enhance the genetic evaluation methods to speed up the genetic improvement for the trait. The objectives of this study were to i) identify the genomic regions associated with resistance to P. salmonis; and ii) identify candidate genes associated with the trait in rainbow trout. We experimentally challenged 2,130 rainbow trout with P. salmonis and genotyped them with a 57 K single nucleotide polymorphism (SNP) array. Resistance to P. salmonis was defined as time to death (TD) and as binary survival (BS). Significant heritabilities were estimated for TD and BS (0.48 ± 0.04 and 0.34 ± 0.04, respectively). A total of 2,047 fish and 26,068 SNPs passed quality control for samples and genotypes. Using a single-step genome wide association analysis (ssGWAS) we identified four genomic regions explaining over 1% of the genetic variance for TD and three for BS. Interestingly, the same genomic region located on Omy27 was found to explain the highest proportion of genetic variance for both traits (2.4 and 1.5% for TD and BS, respectively). The identified SNP in this region is located within an exon of a gene related with actin cytoskeletal organization, a protein exploited by P. salmonis during infection. Other important candidate genes identified are related with innate immune response and oxidative stress. The moderate heritability values estimated in the present study show it is possible to improve resistance to P. salmonis through artificial selection in the rainbow trout population studied here. Furthermore, our results suggest a polygenic genetic architecture for the trait and provide novel insights into the candidate genes underpinning resistance to P. salmonis in O. mykiss.

Multiple Selection Signatures in Farmed Atlantic Salmon Adapted to Different Environments Across Hemispheres

Domestication of Atlantic salmon started approximately 40 years ago, using artificial selection through genetic improvement programs. Selection is likely to have imposed distinctive signatures on the salmon genome, which are often characterized by high genetic differentiation across population and/or reduction in genetic diversity in regions associated to traits under selection. The identification of such selection signatures may give insights into the candidate genomic regions of biological and commercial interest. Here, we used three complementary statistics to detect selection signatures, two haplotype-based (iHS and XP-EHH), and one F_ST-based method (BayeScan) among four populations of Atlantic salmon with a common genetic origin. Several regions were identified for these techniques that harbored genes, such as kind1 and chp2, which have been associated with growth-related traits or the kcnb2 gene related to immune system in Atlantic salmon, making them particularly relevant in the context of aquaculture. Our results provide candidate genes to inform the evolutionary and biological mechanisms controlling complex selected traits in Atlantic salmon.

First Genome-wide Association Analysis for Growth Traits in the Largest Coral Reef-Dwelling Bony Fishes, the Giant Grouper (Epinephelus lanceolatus)

The giant grouper, Epinephelus lanceolatus, is the largest coral reef-dwelling bony fish species. However, despite extremely fast growth performance and the considerable economic importance in this species, its genetic regulation of growth remains unknown. Here, we performed the first genome-wide association study (GWAS) for five growth traits in 289 giant groupers using 42,323 single nucleotide polymorphisms (SNPs) obtained by genotyping-by-sequencing (GBS). We identified a total of 36 growth-related SNPs, of which 11 SNPs reached a genome-wide significance level. The phenotypic variance explained by these SNPs varied from 7.09% for body height to 18.42% for body length. Moreover, 22 quantitative trait loci (QTLs) for growth traits, including nine significant QTLs and 13 suggestive QTLs, were found on multiple chromosomes. Interestingly, the QTL (LG17: 6934451) was shared between body weight and body height, while two significant QTLs (LG7: 22596399 and LG15: 11877836) for body length were consistent with the associated regions of total length at the genome-wide suggestive level. Eight potential candidate genes close to the associated SNPs were selected for expression analysis, of which four genes (phosphatidylinositol transfer protein cytoplasmic 1, protein tyrosine phosphatase receptor type E, alpha/beta hydrolase domain-containing protein 17C, and vascular endothelial growth factor A-A) were differentially expressed and involved in metabolism, development, response stress, etc. This study improves our understanding of the complex genetic architecture of growth in the giant grouper. The results contribute to the selective breeding of grouper species and the conservation of coral reef fishes.

Genome-Wide Association Analysis for Resistance to Infectious Pancreatic Necrosis Virus Identifies Candidate Genes Involved in Viral Replication and Immune Response in Rainbow Trout (Oncorhynchus mykiss)

Infectious pancreatic necrosis (IPN) is a viral disease with considerable negative impact on the rainbow trout (Oncorhynchus mykiss) aquaculture industry. The aim of the present work was to detect genomic regions that explain resistance to infectious pancreatic necrosis virus (IPNV) in rainbow trout. A total of 2,278 fish from 58 full-sib families were challenged with IPNV and 768 individuals were genotyped (488 resistant and 280 susceptible), using a 57K SNP panel Axiom, Affymetrix. A genome-wide association study (GWAS) was performed using the phenotypes time to death (TD) and binary survival (BS), along with the genotypes of the challenged fish using a Bayesian model (Bayes C). Heritabilities for resistance to IPNV estimated using genomic information, were 0.53 and 0.82 for TD and BS, respectively. The Bayesian GWAS detected a SNP located on chromosome 5 explaining 19% of the genetic variance for TD. The proximity of Sentrin-specific protease 5 (SENP5) to this SNP makes it a candidate gene for resistance against IPNV. In case of BS, a SNP located on chromosome 23 was detected explaining 9% of the genetic variance. However, the moderate-low proportion of variance explained by the detected marker leads to the conclusion that the incorporation of all genomic information, through genomic selection, would be the most appropriate approach to accelerate genetic progress for the improvement of resistance against IPNV in rainbow trout.

Genome-Wide Patterns of Population Structure and Linkage Disequilibrium in Farmed Nile Tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus) is one of the most produced farmed fish in the world and represents an important source of protein for human consumption. Farmed Nile tilapia populations are increasingly based on genetically improved stocks, which have been established from admixed populations. To date, there is scarce information about the population genomics of farmed Nile tilapia, assessed by dense single nucleotide polymorphism (SNP) panels. The patterns of linkage disequilibrium (LD) may affect the success of genome-wide association studies (GWAS) and genomic selection (GS), and also provide key information about demographic history of farmed Nile tilapia populations. The objectives of this study were to provide further knowledge about the population structure and LD patterns, as well as, estimate the effective population size (N e ) for three farmed Nile tilapia populations, one from Brazil (POP A) and two from Costa Rica (POP B and POP C). A total of 55 individuals from each population, were genotyped using a 50K SNP panel selected from a whole-genome sequencing (WGS) experiment. The first two principal components explained about 20% of the total variation and clearly differentiated between the three populations. Population genetic structure analysis showed evidence of admixture, especially for POP C. The contemporary N e estimated, based on LD values, ranged from 78 to 159. No differences were observed in the LD decay among populations, with a rapid decrease of r 2 with increasing inter-marker distance. Average r 2 between adjacent SNP pairs ranged from 0.19 to 0.03 for both POP A and C, and 0.20 to 0.03 f or POP B. Based on the number of independent chromosome segments in the Nile tilapia genome, at least 9.4, 7.6, and 4.6K SNPs for POP A, POP B, and POP C respectively, are required for the implementation of GS in the present farmed Nile tilapia populations.

Association Mapping Based on a Common-Garden Migration Experiment Reveals Candidate Genes for Migration Tendency in Brown Trout

A better understanding of the environmental and genetic contribution to migratory behavior and the evolution of traits linked to migration is crucial for fish conservation and fisheries management. Up to date, a few genes with unequivocal influence on the adoption of alternative migration strategies have been identified in salmonids. Here, we used a common garden set-up to measure individual migration distances of generally highly polymorphic brown trout Salmo trutta from two populations. Fish from the assumedly resident population showed clearly shorter migration distances than the fish from the assumed migratory population at the ages of 2 and 3 years. By using two alternative analytical pipelines with 22186 and 18264 SNPs obtained through RAD-sequencing, we searched for associations between individual migration distance, and both called genotypes and genotype probabilities. None of the SNPs showed statistically significant individual effects on migration after correction for multiple testing. By choosing a less stringent threshold, defined as an overlap of the top 0.1% SNPs identified by the analytical pipelines, GAPIT and Angsd, we identified eight candidate genes that are potentially linked to individual migration distance. While our results demonstrate large individual and population level differences in migration distances, the detected genetic associations were weak suggesting that migration traits likely have multigenic control.

Genome-Wide Association Study and Cost-Efficient Genomic Predictions for Growth and Fillet Yield in Nile Tilapia (Oreochromis niloticus)

Fillet yield (FY) and harvest weight (HW) are economically important traits in Nile tilapia production. Genetic improvement of these traits, especially for FY, are lacking, due to the absence of efficient methods to measure the traits without sacrificing fish and the use of information from relatives to selection. However, genomic information could be used by genomic selection to improve traits that are difficult to measure directly in selection candidates, as in the case of FY. The objectives of this study were: (i) to perform genome-wide association studies (GWAS) to dissect the genetic architecture of FY and HW, (ii) to evaluate the accuracy of genotype imputation and (iii) to assess the accuracy of genomic selection using true and imputed low-density (LD) single nucleotide polymorphism (SNP) panels to determine a cost-effective strategy for practical implementation of genomic information in tilapia breeding programs. The data set consisted of 5,866 phenotyped animals and 1,238 genotyped animals (108 parents and 1,130 offspring) using a 50K SNP panel. The GWAS were performed using all genotyped and phenotyped animals. The genotyped imputation was performed from LD panels (LD0.5K, LD1K and LD3K) to high-density panel (HD), using information from parents and 20% of offspring in the reference set and the remaining 80% in the validation set. In addition, we tested the accuracy of genomic selection using true and imputed genotypes comparing the accuracy obtained from pedigree-based best linear unbiased prediction (PBLUP) and genomic predictions. The results from GWAS supports evidence of the polygenic nature of FY and HW. The accuracy of imputation ranged from 0.90 to 0.98 for LD0.5K and LD3K, respectively. The accuracy of genomic prediction outperformed the estimated breeding value from PBLUP. The use of imputation for genomic selection resulted in an increased relative accuracy independent of the trait and LD panel analyzed. The present results suggest that genotype imputation could be a cost-effective strategy for genomic selection in Nile tilapia breeding programs.

Comparative Genomic Analysis of Three Salmonid Species Identifies Functional Candidate Genes Involved in Resistance to the Intracellular Bacterium Piscirickettsia salmonis

Piscirickettsia salmonis is the etiologic agent of salmon rickettsial syndrome (SRS) and is responsible for considerable economic losses in salmon aquaculture. The bacterium affects coho salmon (CS; Oncorhynchus kisutch), Atlantic salmon (AS; Salmo salar), and rainbow trout (RT; Oncorhynchus mykiss) in several countries, including Norway, Canada, Scotland, Ireland, and Chile. We used Bayesian genome-wide association study analyses to investigate the genetic architecture of resistance to P. salmonis in farmed populations of these species. Resistance to SRS was defined as the number of days to death and as binary survival (BS). A total of 828 CS, 2130 RT, and 2601 AS individuals were phenotyped and then genotyped using double-digest restriction site-associated DNA sequencing and 57K and 50K Affymetrix® Axiom® single nucleotide polymorphism (SNP) panels, respectively. Both traits of SRS resistance in CS and RT appeared to be under oligogenic control. In AS, there was evidence of polygenic control of SRS resistance. To identify candidate genes associated with resistance, we applied a comparative genomics approach in which we systematically explored the complete set of genes adjacent to SNPs, which explained more than 1% of the genetic variance of resistance in each salmonid species (533 genes in total). Thus, genes were classified based on the following criteria: i) shared function of their protein domains among species, ii) shared orthology among species, iii) proximity to the SNP explaining the highest proportion of the genetic variance, and iv) presence in more than one genomic region explaining more than 1% of the genetic variance within species. Our results allowed us to identify 120 candidate genes belonging to at least one of the four criteria described above. Of these, 21 of them were part of at least two of the criteria defined above and are suggested to be strong functional candidates influencing P. salmonis resistance. These genes are related to diverse biological processes, such as kinase activity, GTP hydrolysis, helicase activity, lipid metabolism, cytoskeletal dynamics, inflammation, and innate immune response, which seem essential in the host response against P. salmonis infection. These results provide fundamental knowledge on the potential functional genes underpinning resistance against P. salmonis in three salmonid species.

Optimizing Genomic Prediction of Host Resistance to Koi Herpesvirus Disease in Carp

Genomic selection (GS) is increasingly applied in breeding programs of major aquaculture species, enabling improved prediction accuracy and genetic gain compared to pedigree-based approaches. Koi Herpesvirus disease (KHVD) is notifiable by the World Organization for Animal Health and the European Union, causing major economic losses to carp production. GS has potential to breed carp with improved resistance to KHVD, thereby contributing to disease control. In the current study, Restriction-site Associated DNA sequencing (RAD-seq) was applied on a population of 1,425 common carp juveniles which had been challenged with Koi herpes virus, followed by sampling of survivors and mortalities. GS was tested on a wide range of scenarios by varying both SNP densities and the genetic relationships between training and validation sets. The accuracy of correctly identifying KHVD resistant animals using GS was between 8 and 18% higher than pedigree best linear unbiased predictor (pBLUP) depending on the tested scenario. Furthermore, minor decreases in prediction accuracy were observed with decreased SNP density. However, the genetic relationship between the training and validation sets was a key factor in the efficacy of genomic prediction of KHVD resistance in carp, with substantially lower prediction accuracy when the relationships between the training and validation sets did not contain close relatives.

Whole Genome Linkage Disequilibrium and Effective Population Size in a Coho Salmon (Oncorhynchus kisutch) Breeding Population Using a High-Density SNP Array

The estimation of linkage disequilibrium between molecular markers within a population is critical when establishing the minimum number of markers required for association studies, genomic selection, and inferring historical events influencing different populations. This work aimed to evaluate the extent and decay of linkage disequilibrium in a coho salmon breeding population using a high-density SNP array. Linkage disequilibrium was estimated between a total of 93,502 SNPs found in 64 individuals (33 dams and 31 sires) from the breeding population. The markers encompass all 30 coho salmon chromosomes and comprise 1,684.62 Mb of the genome. The average density of markers per chromosome ranged from 48.31 to 66 per 1 Mb. The minor allele frequency averaged 0.26 (with a range from 0.22 to 0.27). The overall average linkage disequilibrium among SNPs pairs measured as r² was 0.10. The Average r² value decreased with increasing physical distance, with values ranging from 0.21 to 0.07 at a distance lower than 1 kb and up to 10 Mb, respectively. An r² threshold of 0.2 was reached at distance of approximately 40 Kb. Chromosomes Okis05, Okis15 and Okis28 showed high levels of linkage disequilibrium (>0.20 at distances lower than 1 Mb). Average r² values were lower than 0.15 for all chromosomes at distances greater than 4 Mb. An effective population size of 43 was estimated for the population 10 generations ago, and 325, for 139 generations ago. Based on the effective number of chromosome segments, we suggest that at least 74,000 SNPs would be necessary for an association mapping study and genomic predictions. Therefore, the SNP panel used allowed us to capture high-resolution information in the farmed coho salmon population. Furthermore, based on the contemporary N_e, a new mate allocation strategy is suggested to increase the effective population size.

High-Density Linkage Map and QTLs for Growth in Snapper (Chrysophrys auratus)

Characterizing the genetic variation underlying phenotypic traits is a central objective in biological research. This research has been hampered in the past by the limited genomic resources available for most non-model species. However, recent advances in sequencing technologies and related genotyping methods are rapidly changing this. Here we report the use of genome-wide SNP data from the ecologically and commercially important marine fish species Chrysophrys auratus (snapper) to 1) construct the first linkage map for this species, 2) scan for growth QTL, and 3) search for putative candidate genes in the surrounding QTL regions. The newly constructed linkage map contained ∼11K SNP markers and is one of the densest maps to date in the fish family Sparidae. Comparisons with genome scaffolds of the recently assembled snapper genome indicated that marker placement was mostly consistent between the scaffolds and linkage map (R = 0.7), but that at fine scales (< 5 cM) some precision limitations occurred. Of the 24 linkage groups, which likely reflect the 24 chromosomes of this species, three were found to contain QTL with genome-wide significance for growth-related traits. A scan of 13 candidate growth genes located the growth hormone, myogenin, and parvalbumin genes within 5.3, 9.6, and 25.0 cM of these QTL, respectively. The linkage map and QTL found in this study will advance the investigation of genome structure and aquaculture breeding efforts in this and related species.

Genomic Selection in Aquaculture: Application, Limitations and Opportunities With Special Reference to Marine Shrimp and Pearl Oysters

Within aquaculture industries, selection based on genomic information (genomic selection) has the profound potential to change genetic improvement programs and production systems. Genomic selection exploits the use of realized genomic relationships among individuals and information from genome-wide markers in close linkage disequilibrium with genes of biological and economic importance. We discuss the technical advances, practical requirements, and commercial applications that have made genomic selection feasible in a range of aquaculture industries, with a particular focus on molluscs (pearl oysters, Pinctada maxima) and marine shrimp (Litopenaeus vannamei and Penaeus monodon). The use of low-cost genome sequencing has enabled cost-effective genotyping on a large scale and is of particular value for species without a reference genome or access to commercial genotyping arrays. We highlight the pitfalls and offer the solutions to the genotyping by sequencing approach and the building of appropriate genetic resources to undertake genomic selection from first-hand experience. We describe the potential to capture large-scale commercial phenotypes based on image analysis and artificial intelligence through machine learning, as inputs for calculation of genomic breeding values. The application of genomic selection over traditional aquatic breeding programs offers significant advantages through being able to accurately predict complex polygenic traits including disease resistance; increasing rates of genetic gain; minimizing inbreeding; and negating potential limiting effects of genotype by environment interactions. Further practical advantages of genomic selection through the use of large-scale communal mating and rearing systems are highlighted, as well as presenting rate-limiting steps that impact on attaining maximum benefits from adopting genomic selection. Genomic selection is now at the tipping point where commercial applications can be readily adopted and offer significant short- and long-term solutions to sustainable and profitable aquaculture industries.

Genome-wide association analysis for body weight identifies candidate genes related to development and metabolism in rainbow trout (Oncorhynchus mykiss)

Growth is one of the most important traits from both a physiological and economic perspective in aquaculture species. Thus, identifying the genomic regions and genes underpinning genetic variation for this trait is of particular interest in several fish species, including rainbow trout. In this work, we perform a genome-wide association study (GWAS) to identify the genomic regions associated with body weight at tagging (BWT) and at 18 months (BW18M) using a dense SNP panel (57 k) and 4596 genotyped rainbow trout from 105 full-sib families belonging to a Chilean breeding population. Analysis was performed by means of single-step GBLUP approach. Genetic variance explained by 20 adjacent SNP windows across the whole genome is reported. To further explore candidate genes, we focused on windows that explained the highest proportion of genetic variance in the top 10 chromosomes for each trait. The main window from the top 10 chromosomes was explored by BLAST using the first and last SNP position of each window to determine the target nucleotide sequence. As expected, the percentage of genetic variance explained by windows was relatively low, due to the polygenic nature of body weight. The most important genomic region for BWT and BW18M were located on chromosomes 15 and 24 and they explained 2.14% and 3.02% of the genetic variance for each trait, respectively. Candidate genes including several growth factors, genes involved in development of skeletal muscle and bone tissue and nutrient metabolism were identified within the associated regions for both traits BWT and BW18M. These results indicate that body weight is polygenic in nature in rainbow trout, with the most important loci explaining as much as 3% of the genetic variance for the trait. The genes identified here represent good candidates for further functional validation to uncover biological mechanisms underlying variation for growth in rainbow trout.

Population Genomic Structure and Genome-Wide Linkage Disequilibrium in Farmed Atlantic Salmon (Salmo salar L.) Using Dense SNP Genotypes

Chilean Farmed Atlantic salmon (Salmo salar) populations were established with individuals of both European and North American origins. These populations are expected to be highly genetically differentiated due to evolutionary history and poor gene flow between ancestral populations from different continents. The extent and decay of linkage disequilibrium (LD) among single nucleotide polymorphism (SNP) impacts the implementation of genome-wide association studies and genomic selection and provides relevant information about demographic processes of fish populations. We assessed the population structure and characterized the extent and decay of LD in three Chilean commercial populations of Atlantic salmon with North American (NAM), Scottish (SCO), and Norwegian (NOR) origin. A total of 123 animals were genotyped using a 159 K SNP Axiom® myDesignTM Genotyping Array. A total of 32 K SNP markers, representing the common SNPs along the three populations after quality control were used. The principal component analysis explained 78.9% of the genetic diversity between populations, clearly discriminating between populations of North American and European origin, and also between European populations. NAM had the lowest effective population size, followed by SCO and NOR. Large differences in the LD decay were observed between populations of North American and European origin. An r 2 threshold of 0.2 was estimated for marker pairs separated by 7,800, 64, and 50 kb in the NAM, SCO, and NOR populations, respectively. In this study we show that this SNP panel can be used to detect association between markers and traits of interests and also to capture high-resolution information for genome-enabled predictions. Also, we suggest the feasibility to achieve similar prediction accuracies using a smaller SNP data set for the NAM population, compared with samples with European origin which would need a higher density SNP array.