Evidence of recent signatures of selection during domestication in an Atlantic salmon population

Unveiling population structure and selection signatures of riverine and genetically improved rohu, Labeo rohita using genome wide SNPs

BACKGROUND

Captive breeding, along with artificial selection can significantly impact population structure by influencing allele frequencies and driving populations towards specific adaptation. Selective sweeps are powerful forces in shaping genetic variation within populations and can drive rapid spread of beneficial alleles while simultaneously reducing genetic diversity in localized regions of the genome. The present work was undertaken to assess the genetic structure and consequences of artificial selection in 10th generation of genetically improved rohu by comparing with wild populations.

METHODS AND RESULTS

The present study used 11,022 high-quality genome wide SNPs to compare the population genetic structure and signatures of selection between Jayanti rohu population and its wild counterpart. Outlier analysis revealed presence of 14 adaptive SNPs, out of which 5 were classified to be under decisive selection pressure. Notably, Jayanti rohu (JR) displayed 297 private alleles exclusive to its population. Chromosomes 7 and 16 emerged as potential hotspots containing a majority of the identified SNPs. Structure and principal component analysis revealed two distinct clusters, effectively distinguishing the JR and wild rohu populations. Phylogenetic analysis indicated a separate cluster of JR population distant from wild groups.

CONCLUSION

The results of present study shall help in elucidating patterns of genetic variation and characterizing selection signatures associated with captive bred and natural populations of rohu. The genomic resources generated through this work shall be helpful in improving the traceability of selectively bred germplasm for developing future strategies of genetic management.

Genome-wide detection of positive and balancing signatures of selection shared by four domesticated rainbow trout populations (Oncorhynchus mykiss)

BACKGROUND

Evolutionary processes leave footprints along the genome over time. Highly homozygous regions may correspond to positive selection of favorable alleles, while maintenance of heterozygous regions may be due to balancing selection phenomena. We analyzed data from 176 fish from four disconnected domestic rainbow trout populations that were genotyped using a high-density Axiom Trout genotyping 665K single nucleotide polymorphism array, including 20 from the US and 156 from three French lines. Using methods based on runs of homozygosity and extended haplotype homozygosity, we detected signatures of selection in these four populations.

RESULTS

Nine genomic regions that included 253 genes were identified as being under positive selection in all four populations Most were located on chromosome 2 but also on chromosomes 12, 15, 16, and 20. In addition, four heterozygous regions that contain 29 genes that are putatively under balancing selection were also shared by the four populations. These were located on chromosomes 10, 13, and 19. Regardless of the homozygous or heterozygous nature of the regions, in each region, we detected several genes that are highly conserved among vertebrates due to their critical roles in cellular and nuclear organization, embryonic development, or immunity. We identified new candidate genes involved in rainbow trout fitness, as well as 17 genes that were previously identified to be under positive selection, 10 of which in other fishes (auts2, atp1b3, zp4, znf135, igf-1α, brd2, col9a2, mrap2, pbx1, and emilin-3).

CONCLUSIONS

Using material from disconnected populations of different origins allowed us to draw a genome-wide map of signatures of positive selection that are shared between these rainbow trout populations, and to identify several regions that are putatively under balancing selection. These results provide a valuable resource for future investigations of the dynamics of genetic diversity and genome evolution during domestication.

Whole‐genome re‐sequencing provides key genomic insights in farmed Arctic charr ( Salvelinus alpinus ) populations of anadromous and landlocked origin from Scandinavia

Arctic charr (Salvelinus alpinus) is a niche-market high-value species for Nordic aquaculture. Similar to other salmonids, both anadromous and landlocked populations are encountered. Whole-genome re-sequencing (22X coverage) was performed on two farmed populations of anadromous (Sigerfjord; n = 24) and landlocked (Arctic Superior; n = 24) origin from Norway and Sweden respectively. More than 5 million SNPs were used to study their genetic diversity and to scan for selection signatures. The two populations were clearly distinguished through principal component analysis, with the mean fixation index being ~0.12. Furthermore, the levels of genomic inbreeding estimated from runs of homozygosity were 6.23% and 8.66% for the Norwegian and the Swedish population respectively. Biological processes that could be linked to selection pressure associated primarily with the anadromous background and/or secondarily with domestication were suggested. Overall, our study provided insights regarding the genetic composition of two main strains of farmed Arctic charr from Scandinavia. At the same time, ample genomic resources were produced in the magnitude of millions of SNPs that could assist the transition of Nordic Arctic charr farming in the genomics era.

Whole-genome resequencing reveals recent signatures of selection in five populations of largemouth bass ( Micropterus salmoides)

Largemouth bass ( Micropterus salmoides) is an economically important fish species in North America, Europe, and China. Various genetic improvement programs and domestication processes have modified its genome sequence through selective pressure, leaving nucleotide signals that can be detected at the genomic level. In this study, we sequenced 149 largemouth bass fish, including protospecies (imported from the US) and improved breeds (four domestic breeding populations from China). We detected genomic regions harboring certain genes associated with improved traits, which may be useful molecular markers for practical domestication, breeding, and selection. Subsequent analyses of genetic diversity and population structure revealed that the improved breeds have undergone more rigorous genetic changes. Through selective signal analysis, we identified hundreds of putative selective sweep regions in each largemouth bass line. Interestingly, we predicted 103 putative candidate genes potentially subjected to selection, including several associated with growth (p sst1 and grb10), early development ( klf9, sp4, and sp8), and immune traits ( pkn2, sept2, bcl6, and ripk2). These candidate genes represent potential genomic landmarks that could be used to improve important traits of biological and commercial interest. In summary, this study provides a genome-wide map of genetic variations and selection footprints in largemouth bass, which may benefit genetic studies and accelerate genetic improvement of this economically important fish.

Genomic signatures of artificial selection in fecundity of Pacific white shrimp, Penaeus vannamei

Penaeusvannamei is the most important economic shrimp in the world. Many selective breeding programs are carried out to improve its production and performance traits. Although significant differences in the reproductive ability of female P. vannamei under artificial breeding conditions have been reported, the genome-wide adaption of the reproductive ability of domesticated female P. vannamei is less investigated. In this study, whole-genome analysis was performed along with pooled DNA sequencing on two fecundity separated bulks, high fecundity bulk (HB), and low fecundity bulk (LB). Each bulk contained 30 individuals from 3 commercial populations. A sequencing depth of >30× was achieved for each bulk, leading to the identification of 625,181 and 629,748 single nucleotide polymorphisms (SNPs) in HB and LB, respectively. Fixation index (Fst) combined with p ratio allowed for the identification of 145 selective sweep regions, with a sequence length of 14.5 Mb, accounting for 0.59% of the genome. Among the 145 selective sweep regions, a total of 64,046 SNPs were identified, and further verification was performed by genotyping 50 candidate SNPs on 60 samples from the offspring of the three populations. Furthermore, 121 genes were screened from the sweep regions. GO annotation and KEGG enrichment analyses showed that partial genes were essential for fecundity regulation. This study provides important information for in-depth investigation of genomic characteristics for long-term selective breeding on the fecundity of female P. vannamei and will also be important for genome-assisted breeding of P. vannamei in the future.

Distinct early life stage gene expression effects of hybridization among European and North American farmed and wild Atlantic salmon populations

Due to multi-generation domestication selection, farmed and wild Atlantic salmon diverge genetically, which raises concerns about potential genetic interactions among escaped farmed and wild populations and disruption of local adaptation through introgression. When farmed strains of distant geographic origin are used, it is unknown whether the genetic consequences posed by escaped farmed fish will be greater than if more locally derived strains are used. Quantifying gene transcript expression differences among divergent farmed, wild and F₁ hybrids under controlled conditions is one of the ways to explore the consequences of hybridization. We compared the transcriptomes of fry at the end of yolk sac absorption of a European (EO) farmed ("StofnFiskur", Norwegian strain), a North American (NA) farmed (Saint John River, NB strain), a Newfoundland (NF) wild population with EO ancestry, and related F₁ hybrids using 44K microarrays. Our findings indicate that the wild population showed greater transcriptome differences from the EO farmed strain than that of the NA farmed strain. We also found the largest differences in global gene expression between the two farmed strains. We detected the fewest differentially expressed transcripts between F₁ hybrids and domesticated/wild maternal strains. We also found that the differentially expressed genes between cross types over-represented GO terms associated with metabolism, development, growth, immune response, and redox homeostasis processes. These findings suggest that the interbreeding of escaped EO/NA farmed and NF wild population would alter gene transcription, and the consequences of hybridization would be greater from escaped EO farmed than NA farmed salmon, resulting in potential effects on the wild populations.

Genomic selection signatures in farmed Colossoma macropomum from tropical and subtropical regions in South America

Tambaqui or cachama (Colossoma macropomum) is one of the most important neotropical freshwater fish used for aquaculture in South America, and its production is concentrated at low latitudes (close to the Equator, 0°), where the water temperature is warm. Therefore, understanding how selection shapes genetic variations and structure in farmed populations is of paramount importance in evolutionary biology. High‐throughput sequencing to generate genome‐wide data for fish species allows for elucidating the genomic basis of adaptation to local or farmed conditions and uncovering genes that control the phenotypes of interest. The present study aimed to detect genomic selection signatures and analyze the genetic variability in farmed populations of tambaqui in South America using single‐nucleotide polymorphism (SNP) markers obtained with double‐digest restriction site‐associated DNA sequencing. Initially, 199 samples of tambaqui farmed populations from different locations (located in Brazil, Colombia, and Peru), a wild population (Amazon River, Brazil), and the base population of a breeding program (Aquaculture Center, CAUNESP, Jaboticabal, SP, Brazil) were genotyped. Observed and expected heterozygosity was 0.231–0.350 and 0.288–0.360, respectively. Significant genetic differentiation was observed using global FST analyses of SNP loci (FST = 0.064, p < 0.050). Farmed populations from Colombia and Peru that differentiated from the Brazilian populations formed distinct groups. Several regions, particularly those harboring the genes of significance to aquaculture, were identified to be under positive selection, suggesting local adaptation to stress under different farming conditions and management practices. Studies aimed at improving the knowledge of genomics of tambaqui farmed populations are essential for aquaculture to gain deeper insights into the evolutionary history of these fish and provide resources for the establishment of breeding programs.

Genomic Signatures of Domestication Selection in the Australasian Snapper (Chrysophrys auratus)

Domestication of teleost fish is a recent development, and in most cases started less than 50 years ago. Shedding light on the genomic changes in key economic traits during the domestication process can provide crucial insights into the evolutionary processes involved and help inform selective breeding programmes. Here we report on the recent domestication of a native marine teleost species in New Zealand, the Australasian snapper (Chrysophrys auratus). Specifically, we use genome-wide data from a three-generation pedigree of this species to uncover genetic signatures of domestication selection for growth. Genotyping-By-Sequencing (GBS) was used to generate genome-wide SNP data from a three-generation pedigree to calculate generation-wide averages of FST between every generation pair. The level of differentiation between generations was further investigated using ADMIXTURE analysis and Principal Component Analysis (PCA). After that, genome scans using Bayescan, LFMM and XP-EHH were applied to identify SNP variants under putative selection following selection for growth. Finally, genes near candidate SNP variants were annotated to gain functional insights. Analysis showed that between generations FST values slightly increased as generational time increased. The extent of these changes was small, and both ADMIXTURE analysis and PCA were unable to form clear clusters. Genome scans revealed a number of SNP outliers, indicative of selection, of which a small number overlapped across analyses methods and populations. Genes of interest within proximity of putative selective SNPs were related to biological functions, and revealed an association with growth, immunity, neural development and behaviour, and tumour repression. Even though few genes overlapped between outlier SNP methods, gene functionalities showed greater overlap between methods. While the genetic changes observed were small in most cases, a number of outlier SNPs could be identified, of which some were found by more than one method. Multiple outlier SNPs appeared to be predominately linked to gene functionalities that modulate growth and survival. Ultimately, the results help to shed light on the genomic changes occurring during the early stages of domestication selection in teleost fish species such as snapper, and will provide useful candidates for the ongoing selective breeding in the future of this and related species.

Detection of selection signatures in the genome of a farmed population of anadromous rainbow trout (Oncorhynchus mykiss)

Domestication processes and artificial selection are likely to leave signatures that can be detected at a molecular level in farmed rainbow trout (Oncorhynchus mykiss). These signatures of selection are genomic regions that contain functional genetic variants conferring a higher fitness to their bearers. We genotyped 749 rainbow trout from a commercial population using a rainbow trout Axiom 57 K SNP array panel and identified putative genomic regions under selection using the pcadapt, Composite Likelihood Ratio (CLR) and Integrated Haplotype Score (iHS) methods. After applying quality-control pipelines and statistical analyses, we detected 12, 96 and 16 SNPs putatively under selection, associated with 96, 781 and 115 candidate genes, respectively. Several of these candidate genes were associated with growth, early development, reproduction, behavior and immune system traits. In addition, some of the SNPs were found in interesting regions located in autosomal inversions on Omy05 and Omy20. These findings could represent a genome-wide map of selection signatures in farmed rainbow trout and could be important in explaining domestication and selection for genetic traits of commercial interest.

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.

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.

Toward Genome-Based Selection in Asian Seabass: What Can We Learn From Other Food Fishes and Farm Animals?

Due to the steadily increasing need for seafood and the plateauing output of fisheries, more fish need to be produced by aquaculture production. In parallel with the improvement of farming methods, elite food fish lines with superior traits for production must be generated by selection programs that utilize cutting-edge tools of genomics. The purpose of this review is to provide a historical overview and status report of a selection program performed on a catadromous predator, the Asian seabass (Lates calcarifer, Bloch 1790) that can change its sex during its lifetime. We describe the practices of wet lab, farm and lab in detail by focusing onto the foundations and achievements of the program. In addition to the approaches used for selection, our review also provides an inventory of genetic/genomic platforms and technologies developed to (i) provide current and future support for the selection process; and (ii) improve our understanding of the biology of the species. Approaches used for the improvement of terrestrial farm animals are used as examples and references, as those processes are far ahead of the ones used in aquaculture and thus they might help those working on fish to select the best possible options and avoid potential pitfalls.

The sockeye salmon genome, transcriptome, and analyses identifying population defining regions of the genome

Sockeye salmon (Oncorhynchus nerka) is a commercially and culturally important species to the people that live along the northern Pacific Ocean coast. There are two main sockeye salmon ecotypes—the ocean-going (anadromous) ecotype and the fresh-water ecotype known as kokanee. The goal of this study was to better understand the population structure of sockeye salmon and identify possible genomic differences among populations and between the two ecotypes. In pursuit of this goal, we generated the first reference sockeye salmon genome assembly and an RNA-seq transcriptome data set to better annotate features of the assembly. Resequenced whole-genomes of 140 sockeye salmon and kokanee were analyzed to understand population structure and identify genomic differences between ecotypes. Three distinct geographic and genetic groups were identified from analyses of the resequencing data. Nucleotide variants in an immunoglobulin heavy chain variable gene cluster on chromosome 26 were found to differentiate the northwestern group from the southern and upper Columbia River groups. Several candidate genes were found to be associated with the kokanee ecotype. Many of these genes were related to ammonia tolerance or vision. Finally, the sex chromosomes of this species were better characterized, and an alternative sex-determination mechanism was identified in a subset of upper Columbia River kokanee.

Comparing Genomic Signatures of Selection Between the Abbassa Strain and Eight Wild Populations of Nile Tilapia (Oreochromis niloticus) in Egypt

Domestication to captive rearing conditions, along with targeted selective breeding have genetic consequences that vary from those in wild environments. Nile tilapia (Oreochromis niloticus) is one of the most translocated and farmed aquaculture species globally, farmed throughout Asia, North and South America, and its African native range. In Egypt, a breeding program established the Abbassa Strain of Nile tilapia (AS) in 2002 based on local broodstock sourced from the Nile River. The AS has been intensively selected for growth and has gone through genetic bottlenecks which have likely shifted levels and composition of genetic diversity within the strain. Consequently, there are questions on the possible genetic impact AS escapees may have on endemic populations of Nile tilapia. However, to date there have been no genetic studies comparing genetic changes in the domesticated AS to local wild populations. This study used 9,827 genome-wide SNPs to investigate population genetic structure and signatures of selection in the AS (generations 9–11) and eight wild Nile tilapia populations from Egypt. SNP analyses identified two major genetic clusters (captive and wild populations), with wild populations showing evidence of isolation-by-distance among the Nile Delta and upstream riverine populations. Between genetic clusters, approximately 6.9% of SNPs were identified as outliers with outliers identified on all 22 O. niloticus chromosomes. A lack of localized outlier clustering on the genome suggests that no genes of major effect were presently detected. The AS has retained high levels of genetic diversity (H_{o_All} = 0.21 ± 0.01; H_{e_All} = 0.23 ± 0.01) when compared to wild populations (H_{o_All} = 0.18 ± 0.01; H_{e_All} = 0.17 ± 0.01) after 11 years of domestication and selective breeding. Additionally, 565 SNPs were unique within the AS line. While these private SNPs may be due to domestication signals or founder effects, it is suspected that introgression with blue tilapia (Oreochromis aureus) has occurred. This study highlights the importance of understanding the effects of domestication in addition to wild population structure to inform future management and dissemination decisions. Furthermore, by conducting a baseline genetic study of wild populations prior to the dissemination of a domestic line, the effects of aquaculture on these populations can be monitored over time.

The Effects of Quantitative Trait Architecture on Detection Power in Short-Term Artificial Selection Experiments

Evolve and resequence (E&R) experiments, in which artificial selection is imposed on organisms in a controlled environment, are becoming an increasingly accessible tool for studying the genetic basis of adaptation. Previous work has assessed how different experimental design parameters affect the power to detect the quantitative trait loci (QTL) that underlie adaptive responses in such experiments, but so far there has been little exploration of how this power varies with the genetic architecture of the evolving traits. In this study, we use forward simulation to build a more realistic model of an E&R experiment in which a quantitative polygenic trait experiences a short, but strong, episode of truncation selection. We study the expected power for QTL detection in such an experiment and how this power is influenced by different aspects of trait architecture, including the number of QTL affecting the trait, their starting frequencies, effect sizes, clustering along a chromosome, dominance, and epistasis patterns. We show that all of these parameters can affect allele frequency dynamics at the QTL and linked loci in complex and often unintuitive ways, and thus influence our power to detect them. One consequence of this is that existing detection methods based on models of independent selective sweeps at individual QTL often have lower detection power than a simple measurement of allele frequency differences before and after selection. Our findings highlight the importance of taking trait architecture into account when designing and interpreting studies of molecular adaptation with temporal data. We provide a customizable modeling framework that will enable researchers to easily simulate E&R experiments with different trait architectures and parameters tuned to their specific study system, allowing for assessment of expected detection power and optimization of experimental design.

Comparative Analysis of the Transcriptome and Distribution of Putative SNPs in Two Rainbow Trout (Oncorhynchus mykiss) Breeding Strains by Using Next-Generation Sequencing

Selective breeding can significantly improve the establishment of sustainable and profitable aquaculture fish farming. For rainbow trout (Oncorhynchus mykiss), one of the main aquaculture coldwater species in Europe, a variety of selected hatchery strains are commercially available. In this study, we investigated the genetic variation between the local Born strain, selected for survival, and the commercially available Silver Steelhead strain, selected for growth. We sequenced the transcriptome of six tissues (gills, head kidney, heart, liver, spleen, and white muscle) from eight healthy individuals per strain, using RNA-seq technology to identify strain-specific gene-expression patterns and single nucleotide polymorphisms (SNPs). In total, 1760 annotated genes were differentially expressed across all tissues. Pathway analysis assigned them to different gene networks. We also identified a set of SNPs, which are heterozygous for one of the two breeding strains: 1229 of which represent polymorphisms over all tissues and individuals. Our data indicate a strong genetic differentiation between Born and Silver Steelhead trout, despite the relatively short time of evolutionary separation of the two breeding strains. The results most likely reflect their specifically adapted genotypes and might contribute to the understanding of differences regarding their robustness toward high stress and pathogenic challenge described in former studies.

Whole genome re-sequencing reveals recent signatures of selection in three strains of farmed Nile tilapia (Oreochromis niloticus)

Nile tilapia belongs to the second most cultivated group of fish in the world, mainly because of its favorable characteristics for production. Genetic improvement programs and domestication process of Nile tilapia may have modified the genome through selective pressure, leaving signals that can be detected at the molecular level. In this work, signatures of selection were identified using genome-wide SNP data, by two haplotype-based (iHS and Rsb) and one F_ST based method. Whole-genome re-sequencing of 326 individuals from three strains (A, B and C) of farmed tilapia maintained in Brazil and Costa Rica was carried out using Illumina HiSeq 2500 technology. After applying conventional SNP-calling and quality-control filters, ~ 1.3 M high-quality SNPs were inferred and used as input for the iHS, Rsb and F_ST based methods. We detected several candidate genes putatively subjected to selection in each strain. A considerable number of these genes are associated with growth (e.g. NCAPG, KLF3, TBC1D1, TTN), early development (e.g. FGFR3, PFKFB3), and immunity traits (e.g. NLRC3, PIGR, MAP1S). These candidate genes represent putative genomic landmarks that could be associated to traits of biological and commercial interest in farmed Nile tilapia.

Relative genomic impacts of translocation history, hatchery practices, and farm selection in Pacific oyster Crassostrea gigas throughout the Northern Hemisphere

Pacific oyster Crassostrea gigas, endemic to coastal Asia, has been translocated globally throughout the past century, resulting in self-sustaining introduced populations (naturalized). Oyster aquaculture industries in many parts of the world depend on commercially available seed (hatchery-farmed) or naturalized/wild oysters to move onto a farm (naturalized-farmed). It is therefore important to understand genetic variation among populations and farm types. Here, we genotype naturalized/wild populations from France, Japan, China, and most extensively in coastal British Columbia, Canada. We also genotype cultured populations from throughout the Northern Hemisphere to compare with naturalized populations. In total, 16,942 markers were identified using double-digest RAD-sequencing in 182 naturalized, 112 hatchery-farmed, and 72 naturalized-farmed oysters (n = 366). Consistent with previous studies, very low genetic differentiation was observed around Vancouver Island (mean F _ST = 0.0019) and low differentiation between countries in the Japan-Canada-France historical translocation lineage (France-Canada F _ST = 0.0024; Japan-Canada F _ST = 0.0060). Chinese populations were more differentiated (China-Japan F _ST = 0.0241). Hatchery-propagated populations had higher interindividual relatedness suggesting family structure. Within-population inbreeding was not detected on farms, but nucleotide diversity and polymorphism rate were lower in one farm population. Moving oysters from nature onto farms did not result in strong within-generation selection. Private alleles at substantial frequency were identified in several hatchery populations grown in BC, suggesting nonlocal origins. Tests of selection identified outlier loci consistent with selective differences associated with domestication, in some cases consistently identified in multiple farms. Top outlier candidates were nearby genes involved in calcium signaling and calmodulin activity. Implications of potential introgression from hatchery-farmed oysters depend on whether naturalized populations are valued as a locally adapted resource or as an introduced, invasive species. Given the value of the industry in BC and the challenges the industry faces (e.g., climate change, crop losses, biotic stressors), this remains an important question.

Changed Patterns of Genomic Variation Following Recent Domestication: Selection Sweeps in Farmed Atlantic Salmon

The introduction of wild Atlantic salmon into captivity, and their subsequent artificial selection for production traits, has caused phenotypic differences between domesticated fish and their wild counterparts. Identification of regions of the genome underling these changes offers the promise of characterizing the early biological consequences of domestication. In the current study, we sequenced a population of farmed European Atlantic salmon and compared the observed patterns of SNP variation to those found in conspecific wild populations. This identified 139 genomic regions that contained significantly elevated SNP homozygosity in farmed fish when compared to their wild counterparts. The most extreme was adjacent to versican, a gene involved in control of neural crest cell migration. To control for false positive signals, a second and independent dataset of farmed and wild European Atlantic salmon was assessed using the same methodology. A total of 81 outlier regions detected in the first dataset showed significantly reduced homozygosity within the second one, strongly suggesting the genomic regions identified are enriched for true selection sweeps. Examination of the associated genes identified a number previously characterized as targets of selection in other domestic species and that have roles in development, behavior and olfactory system. These include arcvf, sema6, errb4, id2-like, and 6n1-like genes. Finally, we searched for evidence of parallel sweeps using a farmed population of North American origin. This failed to detect a convincing overlap to the putative sweeps present in European populations, suggesting the factors that drive patterns of variation under domestication and early artificial selection were largely independent. This is the first analysis on domestication of aquaculture species exploiting whole-genome sequence data and resulted in the identification of sweeps common to multiple independent populations of farmed European Atlantic salmon.

Genetic Differentiation in Hatchery and Stocked Populations of Sea Trout in the Southern Baltic: Selection Evidence at SNP Loci

The impacts and interactions between hatchery-bred fish and wild fish populations has been a topic of active investigation in recent decades. In some instances, the benefits of stocking can be overshadowed by negative effects such as genetic introgression with natural populations, loss of genetic diversity, and dilution of local adaptations. Methods that facilitate the identification of stocked fish enable us to estimate not only the effectiveness of stocking but also the level of natural reproduction and the degree of hybridization. The longest Baltic river, the Vistula, also has the second highest discharge. Historically, it hosted numerous populations of the anadromous form of brown trout (sea trout); however, dam construction has since interfered with and reduced spawning migration to a rate that is much lower than before. Reduced spawning has resulted in a population collapse and a negative flow-on effect on commercial catches. In response, Poland (along with many other Baltic countries) initiated an intensive stocking program which continues today and which sees the average annual release of 700,000 smolts. As a consequence, today’s main-river and inshore catches come from stock-enhanced populations. High-throughput single-nucleotide polymorphism (SNP) genotyping was performed on samples of sea trout from southern Baltic populations; results suggest that a significant portion of the sea trout catches in the Vistula mouth region have direct hatchery origin and indicate the presence of Pomeranian specimens. SNP loci identified as outliers indicate a potential selection pressure that may be related with effects of hatchery breeding and mixing with natural populations. The brown trout SNP array applied in this study showed high effectiveness not only for population differentiation, but more importantly, it emerged as a sensitive tool to provide evidence of detection selection.

High-Throughput Single Nucleotide Polymorphism (SNP) Discovery and Validation Through Whole-Genome Resequencing in Nile Tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus) is the second most important farmed fish in the world and a sustainable source of protein for human consumption. Several genetic improvement programs have been established for this species in the world. Currently, the estimation of genetic merit of breeders is typically based on genealogical and phenotypic information. Genome-wide information can be exploited to efficiently incorporate traits that are difficult to measure into the breeding goal. Thus, single nucleotide polymorphisms (SNPs) are required to investigate phenotype-genotype associations and determine the genomic basis of economically important traits. We performed de novo SNP discovery in three different populations of farmed Nile tilapia. A total of 29.9 million non-redundant SNPs were identified through Illumina (HiSeq 2500) whole-genome resequencing of 326 individual samples. After applying several filtering steps, including removing SNP based on genotype and site quality, presence of Mendelian errors, and non-unique position in the genome, a total of 50,000 high-quality SNPs were selected for the development of a custom Illumina BeadChip SNP panel. These SNPs were highly informative in the three populations analyzed showing between 43,869 (94%) and 46,139 (99%) SNPs in Hardy-Weinberg Equilibrium; 37,843 (76%) and 45,171(90%) SNPs with a minor allele frequency (MAF) higher than 0.05; and 43,450 (87%) and 46,570 (93%) SNPs with a MAF higher than 0.01. The 50K SNP panel developed in the current work will be useful for the dissection of economically relevant traits, enhancing breeding programs through genomic selection, as well as supporting genetic studies in farmed populations of Nile tilapia using dense genome-wide information.

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.

Atlantic salmon (Salmo salar L.) genetics in the 21st century: taking leaps forward in aquaculture and biological understanding

Atlantic salmon (Salmo salar L.) is among the most iconic and economically important fish species and was the first member of Salmonidae to have a high-quality reference genome assembly published. Advances in genomics have become increasingly central to the genetic improvement of farmed Atlantic salmon as well as conservation of wild salmon stocks. The salmon genome has also been pivotal in shaping our understanding of the evolutionary and functional consequences arising from an ancestral whole-genome duplication event characterising all Salmonidae members. Here, we provide a review of the current status of Atlantic salmon genetics and genomics, focussed on progress made from genome-wide research aimed at improving aquaculture production and enhancing understanding of salmonid ecology, physiology and evolution. We present our views on the future direction of salmon genomics, including the role of emerging technologies (e.g. genome editing) in elucidating genetic features that underpin functional variation in traits of commercial and evolutionary importance.

Comparing genomic signatures of domestication in two Atlantic salmon (Salmo salar L.) populations with different geographical origins

Selective breeding and genetic improvement have left detectable signatures on the genomes of domestic species. The elucidation of such signatures is fundamental for detecting genomic regions of biological relevance to domestication and improving management practices. In aquaculture, domestication was carried out independently in different locations worldwide, which provides opportunities to study the parallel effects of domestication on the genome of individuals that have been selected for similar traits. In this study, we aimed to detect potential genomic signatures of domestication in two independent pairs of wild/domesticated Atlantic salmon populations of Canadian and Scottish origins, respectively. Putative genomic regions under divergent selection were investigated using a 200K SNP array by combining three different statistical methods based either on allele frequencies (LFMM, Bayescan) or haplotype differentiation (Rsb). We identified 337 and 270 SNPs potentially under divergent selection in wild and hatchery populations of Canadian and Scottish origins, respectively. We observed little overlap between results obtained from different statistical methods, highlighting the need to test complementary approaches for detecting a broad range of genomic footprints of selection. The vast majority of the outliers detected were population-specific but we found four candidate genes that were shared between the populations. We propose that these candidate genes may play a role in the parallel process of domestication. Overall, our results suggest that genetic drift may have override the effect of artificial selection and/or point toward a different genetic basis underlying the expression of similar traits in different domesticated strains. Finally, it is likely that domestication may predominantly target polygenic traits (e.g., growth) such that its genomic impact might be more difficult to detect with methods assuming selective sweeps.

Development and Validation of 58K SNP-Array and High-Density Linkage Map in Nile Tilapia (O. niloticus)

Despite being the second most important aquaculture species in the world accounting for 7.4% of global production in 2015, tilapia aquaculture has lacked genomic tools like SNP-arrays and high-density linkage maps to improve selection accuracy and accelerate genetic progress. In this paper, we describe the development of a genotyping array containing more than 58,000 SNPs for Nile tilapia (Oreochromis niloticus). SNPs were identified from whole genome resequencing of 32 individuals from the commercial population of the Genomar strain, and were selected for the SNP-array based on polymorphic information content and physical distribution across the genome using the Orenil1.1 genome assembly as reference sequence. SNP-performance was evaluated by genotyping 4991 individuals, including 689 offspring belonging to 41 full-sib families, which revealed high-quality genotype data for 43,588 SNPs. A preliminary genetic linkage map was constructed using Lepmap2 which in turn was integrated with information from the O_niloticus_UMD1 genome assembly to produce an integrated physical and genetic linkage map comprising 40,186 SNPs distributed across 22 linkage groups (LGs). Around one-third of the LGs showed a different recombination rate between sexes, with the female being greater than the male map by a factor of 1.2 (1632.9 to 1359.6 cM, respectively), with most LGs displaying a sigmoid recombination profile. Finally, the sex-determining locus was mapped to position 40.53 cM on LG23, in the vicinity of the anti-Müllerian hormone (amh) gene. These new resources has the potential to greatly influence and improve the genetic gain when applying genomic selection and surpass the difficulties of efficient selection for invasively measured traits in Nile tilapia.

Comparisons of reproductive function and fatty acid fillet quality between triploid and diploid farm Atlantic salmon (Salmo salar)

Triploidy could prevent escaped farm salmon breeding in the wild, while also improving nutrient quality within farmed fillets. Despite these potential advantages, triploid Atlantic salmon have not been widely used in aquaculture, and their reproductive function has yet to be fully evaluated. Here, we compare reproductive function and fillet composition between triploid and diploid farm salmon under standard aquaculture rearing conditions. We show that female triploids are sterile and do not develop gonads. By contrast, males produce large numbers of motile spermatozoa capable of fertilizing wild salmon eggs. However, compared with diploids, reproductive development and survival rates of eggs fertilized by triploid males were significantly reduced, with less than 1% of eggs sired by triploid males reaching late-eyed stages of development. Analyses of fillets showed that total lipid and fatty acid quantities were significantly lower in triploid than in diploid Atlantic salmon fillets. However, when fatty acids were normalized to total lipid content, triploid fillets had significantly higher relative levels of important omega-3 long-chain polyunsaturated fatty acids. Our results show that: (i) escaped triploid farm salmon are very unlikely to reproduce in the wild and (ii) if able to match diploid fillet lipid content, triploid farm salmon could achieve better fillet quality in terms of essential fatty acids.

A genome scan for selection signatures comparing farmed Atlantic salmon with two wild populations: Testing colocalization among outlier markers, candidate genes, and quantitative trait loci for production traits

Comparative genome scans can be used to identify chromosome regions, but not traits, that are putatively under selection. Identification of targeted traits may be more likely in recently domesticated populations under strong artificial selection for increased production. We used a North American Atlantic salmon 6K SNP dataset to locate genome regions of an aquaculture strain (Saint John River) that were highly diverged from that of its putative wild founder population (Tobique River). First, admixed individuals with partial European ancestry were detected using STRUCTURE and removed from the dataset. Outlier loci were then identified as those showing extreme differentiation between the aquaculture population and the founder population. All Arlequin methods identified an overlapping subset of 17 outlier loci, three of which were also identified by BayeScan. Many outlier loci were near candidate genes and some were near published quantitative trait loci (QTLs) for growth, appetite, maturity, or disease resistance. Parallel comparisons using a wild, nonfounder population (Stewiacke River) yielded only one overlapping outlier locus as well as a known maturity QTL. We conclude that genome scans comparing a recently domesticated strain with its wild founder population can facilitate identification of candidate genes for traits known to have been under strong artificial selection.

Plasticity in growth of farmed and wild Atlantic salmon: is the increased growth rate of farmed salmon caused by evolutionary adaptations to the commercial diet?

Background

Domestication of Atlantic salmon for commercial aquaculture has resulted in farmed salmon displaying substantially higher growth rates than wild salmon under farming conditions. In contrast, growth differences between farmed and wild salmon are much smaller when compared in the wild. The mechanisms underlying this contrast between environments remain largely unknown. It is possible that farmed salmon have adapted to the high-energy pellets developed specifically for aquaculture, contributing to inflated growth differences when fed on this diet. We studied growth and survival of 15 families of farmed, wild and F1 hybrid salmon fed three contrasting diets under hatchery conditions; a commercial salmon pellet diet, a commercial carp pellet diet, and a mixed natural diet consisting of preserved invertebrates commonly found in Norwegian rivers.

Results

For all groups, despite equal numbers of calories presented by all diets, overall growth reductions as high 68 and 83%, relative to the salmon diet was observed in the carp and natural diet treatments, respectively. Farmed salmon outgrew hybrid (intermediate) and wild salmon in all treatments. The relative growth difference between wild and farmed fish was highest in the carp diet (1: 2.1), intermediate in the salmon diet (1:1.9) and lowest in the natural diet (1:1.6). However, this trend was non-significant, and all groups displayed similar growth reaction norms and plasticity towards differing diets across the treatments.

Conclusions

No indication of genetic-based adaptation to the form or nutritional content of commercial salmon diets was detected in the farmed salmon. Therefore, we conclude that diet alone, at least in the absence of other environmental stressors, is not the primary cause for the large contrast in growth differences between farmed and wild salmon in the hatchery and wild. Additionally, we conclude that genetically-increased appetite is likely to be the primary reason why farmed salmon display higher growth rates than wild salmon when fed ad lib rations under hatchery conditions. Our results contribute towards an understanding of the potential genetic changes that have occurred in farmed salmon in response to domestication, and the potential mechanisms underpinning genetic and ecological interactions between farmed escapees and wild salmonids.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0841-7) contains supplementary material, which is available to authorized users.

Genomewide single nucleotide polymorphism discovery in Atlantic salmon (Salmo salar): validation in wild and farmed American and European populations

A considerable number of single nucleotide polymorphisms (SNPs) are required to elucidate genotype-phenotype associations and determine the molecular basis of important traits. In this work, we carried out de novo SNP discovery accounting for both genome duplication and genetic variation from American and European salmon populations. A total of 9 736 473 nonredundant SNPs were identified across a set of 20 fish by whole-genome sequencing. After applying six bioinformatic filtering steps, 200 K SNPs were selected to develop an Affymetrix Axiom(®) myDesign Custom Array. This array was used to genotype 480 fish representing wild and farmed salmon from Europe, North America and Chile. A total of 159 099 (79.6%) SNPs were validated as high quality based on clustering properties. A total of 151 509 validated SNPs showed a unique position in the genome. When comparing these SNPs against 238 572 markers currently available in two other Atlantic salmon arrays, only 4.6% of the SNP overlapped with the panel developed in this study. This novel high-density SNP panel will be very useful for the dissection of economically and ecologically relevant traits, enhancing breeding programmes through genomic selection as well as supporting genetic studies in both wild and farmed populations of Atlantic salmon using high-resolution genomewide information.