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

Genetic diversity patterns in farmed rainbow trout (Oncorhynchus mykiss) populations using genome-wide SNP and haplotype data

Rainbow trout is one of the most popular aquaculture species worldwide, with a long history of domestication. However, limited information exists about the genetic diversity of farmed rainbow trout populations globally, with most available reports relying on low-throughput genotyping technologies. Notably, no information exists about the genetic diversity status of farmed rainbow trout in Sweden. Double-digest restriction-site-associated DNA sequencing was performed on more than 500 broodfish from two leading producers in Sweden and from the country's national breeding program. Following the detection of single nucleotide polymorphisms (SNPs), genetic diversity was studied by using either individual SNPs (n = 8680; one SNP retained per 300 bp sequence reads) or through SNP haplotypes (n = 20 558; all SNPs retained in 300 bp sequence reads). Similar amounts of genetic diversity were found amongst the three populations when individual SNPs were used. Furthermore, principal component analysis and discriminant analysis of principal components suggested two genetic clusters with the two industry populations grouped together. Genetic differentiation based on the FST fixation index was ~0.01 between the industry populations and ~0.05 when those were compared with the breeding program. Preliminary estimates of effective population size (Ne ) and inbreeding (based on runs of homozygosity; FROH ) were similar amongst the three populations (Ne  ≈ 50-80; median FROH  ≈ 0.11). Finally, the haplotype-based analysis suggested that animals from the breeding program had higher shared coancestry levels than those from the other two populations. Overall, our study provides novel insights into the genetic diversity and structure of Sweden's three main farmed rainbow trout populations, which could guide their future management.

Sex-Associated SNP Confirmation of Sex-Reversed Male Farmed Japanese Flounder Paralichthys olivaceus

Female Japanese flounder Paralichthys olivaceus grow more rapidly than the male. The goal of all-female commercial production requires an efficient method of genetic sex identification. We conducted genome-wide association analysis of female and male farmed Japanese flounder (n = 24 per phenotypic sex) and found all regions of chromosome 24 to be significantly associated with phenotypic sex, suggesting it as the sex chromosome. Genetic sex was identified based on single nucleotide polymorphisms (SNP) on chromosome 24 (n = 3568) using multidimensional scaling analysis, and individuals were clearly separated according to sex by the first dimension. The 61 SNPs most highly associated with sex were selected, and an amplicon-based SNP panel was developed. This was used to determine genetic sex of 39 females and 40 males. Eleven phenotypic males were assigned as female with XX genotype, suggesting sex reversal. Genetic sex was also assessed based on the indel of the amh gene promoter, which is the major candidate sex gene of Japanese flounder. We found four SNPs perfectly associated with genotypic sex in the sex-associated SNP panel, one of which was located in exon 2 of the amh gene. Along with the indel of the amh gene promoter, the sex-associated SNP panel will be of value in identifying genetic sex of farmed Japanese flounder. Molecular sexing will facilitate all-female production by breeding sex-reversed males.

Analysis of Autozygosity Using Whole-Genome Sequence Data of Full-Sib Families in Pikeperch (Sander lucioperca)

Pikeperch (Sander lucioperca) has emerged as a high value species to the aquaculture industry. However, its farming techniques are at an early stage and its production is often performed without a selective breeding program, potentially leading to high levels of inbreeding. In this study, we identified and characterized autozygosity based on genome-wide runs of homozygosity (ROH) on a sample of parental and offspring individuals, determined effective population size (N_e), and assessed relatedness among parental individuals. A mean of 2,235 ± 526 and 1,841 ± 363 ROH segments per individual, resulting in a mean inbreeding coefficient of 0.33 ± 0.06 and 0.25 ± 0.06 were estimated for the progeny and parents, respectively. N_e was about 12 until four generations ago and at most 106 for 63 generations in the past, with varying genetic relatedness amongst the parents. This study shows the importance of genomic information when family relationships are unknown and the need of selective breeding programs for reproductive management decisions in the aquaculture industry.

Estimates of recent and historical effective population size in turbot, seabream, seabass and carp selective breeding programmes

Background

The high fecundity of fish species allows intense selection to be practised and therefore leads to fast genetic gains. Based on this, numerous selective breeding programmes have been started in Europe in the last decades, but in general, little is known about how the base populations of breeders have been built. Such knowledge is important because base populations can be created from very few individuals, which can lead to small effective population sizes and associated reductions in genetic variability. In this study, we used genomic information that was recently made available for turbot (Scophthalmus maximus), gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and common carp (Cyprinus carpio) to obtain accurate estimates of the effective size for commercial populations.

Methods

Restriction-site associated DNA sequencing data were used to estimate current and historical effective population sizes. We used a novel method that considers the linkage disequilibrium spectrum for the whole range of genetic distances between all pairs of single nucleotide polymorphisms (SNPs), and thus accounts for potential fluctuations in population size over time.

Results

Our results show that the current effective population size for these populations is small (equal to or less than 50 fish), potentially putting the sustainability of the breeding programmes at risk. We have also detected important drops in effective population size about five to nine generations ago, most likely as a result of domestication and the start of selective breeding programmes for these species in Europe.

Conclusions

Our findings highlight the need to broaden the genetic composition of the base populations from which selection programmes start, and suggest that measures designed to increase effective population size within all farmed populations analysed here should be implemented in order to manage genetic variability and ensure the sustainability of the breeding programmes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00680-9.

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.

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.

Development and Validation of an Open Access SNP Array for Nile Tilapia (Oreochromis niloticus)

Tilapia are among the most important farmed fish species worldwide, and are fundamental for the food security of many developing countries. Several genetically improved Nile tilapia (Oreochromis niloticus) strains exist, such as the iconic Genetically Improved Farmed Tilapia (GIFT), and breeding programs typically follow classical pedigree-based selection. The use of genome-wide single-nucleotide polymorphism (SNP) data can enable an understanding of the genetic architecture of economically important traits and the acceleration of genetic gain via genomic selection. Due to the global importance and diversity of Nile tilapia, an open access SNP array would be beneficial for aquaculture research and production. In the current study, a ∼65K SNP array was designed based on SNPs discovered from whole-genome sequence data from a GIFT breeding nucleus population and the overlap with SNP datasets from wild fish populations and several other farmed Nile tilapia strains. The SNP array was applied to clearly distinguish between different tilapia populations across Asia and Africa, with at least ∼30,000 SNPs segregating in each of the diverse population samples tested. It is anticipated that this SNP array will be an enabling tool for population genetics and tilapia breeding research, facilitating consistency and comparison of results across studies.

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.

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.

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.

Genome-wide estimates of genetic diversity, inbreeding and effective size of experimental and commercial rainbow trout lines undergoing selective breeding

BACKGROUND

Selective breeding is a relatively recent practice in aquaculture species compared to terrestrial livestock. Nevertheless, the genetic variability of farmed salmonid lines, which have been selected for several generations, should be assessed. Indeed, a significant decrease in genetic variability due to high selection intensity could have occurred, potentially jeopardizing the long-term genetic progress as well as the adaptive capacities of populations facing change(s) in the environment. Thus, it is important to evaluate the impact of selection practices on genetic diversity to limit future inbreeding. The current study presents an analysis of genetic diversity within and between six French rainbow trout (Oncorhynchus mykiss) experimental or commercial lines based on a medium-density single nucleotide polymorphism (SNP) chip and various molecular genetic indicators: fixation index (FST), linkage disequilibrium (LD), effective population size (Ne) and inbreeding coefficient derived from runs of homozygosity (ROH).

RESULTS

Our results showed a moderate level of genetic differentiation between selected lines (FST ranging from 0.08 to 0.15). LD declined rapidly over the first 100 kb, but then remained quite high at long distances, leading to low estimates of Ne in the last generation ranging from 24 to 68 depending on the line and methodology considered. These results were consistent with inbreeding estimates that varied from 10.0% in an unselected experimental line to 19.5% in a commercial line, and which are clearly higher than corresponding estimates in ruminants or pigs. In addition, strong variations in LD and inbreeding were observed along the genome that may be due to differences in local rates of recombination or due to key genes that tended to have fixed favorable alleles for domestication or production.

CONCLUSIONS

This is the first report on ROH for any aquaculture species. Inbreeding appeared to be moderate to high in the six French rainbow trout lines, due to founder effects at the start of the breeding programs, but also likely to sweepstakes reproductive success in addition to selection for the selected lines. Efficient management of inbreeding is a major goal in breeding programs to ensure that populations can adapt to future breeding objectives and SNP information can be used to manage the rate at which inbreeding builds up in the fish genome.

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.