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.

A breeding program for Nile tilapia in Brazil: Results from nine generations of selection to increase the growth rate in cages.

Variance components and heritabilities for daily weight gain (DWG) were estimated for Nile tilapia farmed in cages across nine generations (G1-G9) of selection in a breeding program in Brazil. DWG was measured in 16,272 accumulated tagged animals representing 535 full- and half-sib families of Nile tilapia under cage farming. The additive genetic variance showed a slight variation (0.051-0.066), and heritability estimates ranged from 0.20 to 0.33. The common environmental effect accounted for a higher proportion of the total variance in DWG, especially in the last generations (6%-24%). A genetic trend based on all data available showed a substantial increase in the DWG (about 3.3% per generation) of Nile tilapia across nine generations of selection. Furthermore, our results demonstrate ample scope for further genetic improvement.

Multi-trait GWAS using imputed high-density genotypes from whole-genome sequencing identifies genes associated with body traits in Nile tilapia.

BACKGROUND: Body traits are generally controlled by several genes in vertebrates (i.e. polygenes), which in turn make them difficult to identify through association mapping. Increasing the power of association studies by combining approaches such as genotype imputation and multi-trait analysis improves the ability to detect quantitative trait loci associated with polygenic traits, such as body traits. RESULTS: A multi-trait genome-wide association study (mtGWAS) was performed to identify quantitative trait loci (QTL) and genes associated with body traits in Nile tilapia (Oreochromis niloticus) using genotypes imputed to whole-genome sequences (WGS). To increase the statistical power of mtGWAS for the detection of genetic associations, summary statistics from single-trait genome-wide association studies (stGWAS) for eight different body traits recorded in 1309 animals were used. The mtGWAS increased the statistical power from the original sample size from 13 to 44%, depending on the trait analyzed. The better resolution of the WGS data, combined with the increased power of the mtGWAS approach, allowed the detection of significant markers which were not previously found in the stGWAS. Some of the lead single nucleotide polymorphisms (SNPs) were found within important functional candidate genes previously associated with growth-related traits in other terrestrial species. For instance, we identified SNP within the α1,6-fucosyltransferase (FUT8), solute carrier family 4 member 2 (SLC4A2), A disintegrin and metalloproteinase with thrombospondin motifs 9 (ADAMTS9) and heart development protein with EGF like domains 1 (HEG1) genes, which have been associated with average daily gain in sheep, osteopetrosis in cattle, chest size in goats, and growth and meat quality in sheep, respectively. CONCLUSIONS: The high-resolution mtGWAS presented here allowed the identification of significant SNPs, linked to strong functional candidate genes, associated with body traits in Nile tilapia. These results provide further insights about the genetic variants and genes underlying body trait variation in cichlid fish with high accuracy and strong statistical support.

Effects of genotype × environment interaction on the estimation of genetic parameters and gains in Nile tilapia.

Nile tilapia (Oreochromis niloticus) is the major fish species produced in Brazil, a country with a vast territory and great climate diversity. This study assessed the effects of the genotype × environment interaction on heritability estimates and selection responses in Nile tilapia (Tilamax strain) cultivated in earthen ponds and net cages. The weight at harvest, trunk length, and head percentage of 4400 individuals were determined. Trait heritabilities were higher in pond fish (0.27-0.52) than in caged fish (0.09-0.33). Genetic correlations between farming systems were lower than 0.5 for the three traits. The rank position of the top 10 families differed according to the environment, as did the response to direct and indirect selection. The results revealed significant genotype × environment effects on the heritability of Nile tilapia farmed under different systems.

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

Estimates of Autozygosity Through Runs of Homozygosity in Farmed Coho Salmon.

The characterization of runs of homozygosity (ROH), using high-density single nucleotide polymorphisms (SNPs) allows inferences to be made about the past demographic history of animal populations and the genomic ROH has become a common approach to characterize the inbreeding. We aimed to analyze and characterize ROH patterns and compare different genomic and pedigree-based methods to estimate the inbreeding coefficient in two pure lines (POP A and B) and one recently admixed line (POP C) of coho salmon (Oncorhynchus kisutch) breeding nuclei, genotyped using a 200 K Affymetrix Axiom® myDesign Custom SNP Array. A large number and greater mean length of ROH were found for the two "pure" lines and the recently admixed line (POP C) showed the lowest number and smaller mean length of ROH. The ROH analysis for different length classes suggests that all three coho salmon lines the genome is largely composed of a high number of short segments (<4 Mb), and for POP C no segment >16 Mb was found. A high variable number of ROH, mean length and inbreeding values across chromosomes; positively the consequence of artificial selection. Pedigree-based inbreeding values tended to underestimate genomic-based inbreeding levels, which in turn varied depending on the method used for estimation. The high positive correlations between different genomic-based inbreeding coefficients suggest that they are consistent and may be more accurate than pedigree-based methods, given that they capture information from past and more recent demographic events, even when there are no pedigree records available.

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.

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.

Heritability of hip dysplasia: Preliminary results for German Shepherd dogs in Brazil.

The heritability of canine hip dysplasia in German Shepherd dogs was estimated using Bayesian methods. Data on hip score and status of 1632 dogs born from 1990 to 2013 were provided by the Brazilian Society of German Shepherd Breeders. Heritability estimates (mean ± standard deviation) were 0.1979 ± 0.058 for hip scores and 0.187 ± 0.055 for hip status. We observed no phenotypic trends and a small rate of genetic trend (0.52%) according to the year of birth, probably because of ineffective phenotypic selection and absence of genetic selection. The heritability estimates in this study can be used to achieve effective selective breeding and genetic gains, which in turn can result in improvements in dog health and welfare.

Fine Mapping Using Whole-Genome Sequencing Confirms Anti-Müllerian Hormone as a Major Gene for Sex Determination in Farmed Nile Tilapia (Oreochromis niloticus L.).

Nile tilapia (Oreochromis niloticus) is one of the most cultivated and economically important species in world aquaculture. Intensive production promotes the use of monosex animals, due to an important dimorphism that favors male growth. Currently, the main mechanism to obtain all-male populations is the use of hormones in feeding during larval and fry phases. Identifying genomic regions associated with sex determination in Nile tilapia is a research topic of great interest. The objective of this study was to identify genomic variants associated with sex determination in three commercial populations of Nile tilapia. Whole-genome sequencing of 326 individuals was performed, and a total of 2.4 million high-quality bi-allelic single nucleotide polymorphisms (SNPs) were identified after quality control. A genome-wide association study (GWAS) was conducted to identify markers associated with the binary sex trait (males = 1; females = 0). A mixed logistic regression GWAS model was fitted and a genome-wide significant signal comprising 36 SNPs, spanning a genomic region of 536 kb in chromosome 23 was identified. Ten out of these 36 genetic variants intercept the anti-Müllerian (Amh) hormone gene. Other significant SNPs were located in the neighboring Amh gene region. This gene has been strongly associated with sex determination in several vertebrate species, playing an essential role in the differentiation of male and female reproductive tissue in early stages of development. This finding provides useful information to better understand the genetic mechanisms underlying sex determination in Nile tilapia.

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

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 2 was 0.10. The Average r 2 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 2 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 2 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.

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.

Single-step genomic evaluation improves accuracy of breeding value predictions for resistance to infectious pancreatic necrosis virus in rainbow trout.

The aim of this study was to compare the accuracy of breeding values (EBVs) predicted using the traditional pedigree based Best Linear Unbiased Prediction (PBLUP) and the single-step genomic Best Linear Unbiased Prediction (ssGBLUP) for resistance against infectious pancreatic necrosis virus (IPNV) in rainbow trout. A total of 2278 animals were challenged against IPNV and 768 individuals were genotyped using a 57 K single nucleotide polymorphism array for rainbow trout. Accuracies for both methods were assessed using five-fold cross-validation. The heritabilities were higher for PBLUP compared to ssGBLUP. The ssGBLUP accuracies outperformed PBLUP in 7 and 11% for days to death and binary survival, respectively. The ssGBLUP could be an alternative approach to improve the accuracy of breeding values for resistance against infectious pancreatic necrosis virus in rainbow trout, using information from genotyped and non-genotyped animals.

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

Piscirickettsia salmonis is one of the main infectious diseases affecting coho salmon (Oncorhynchus kisutch) farming, and current treatments have been ineffective for the control of this disease. Genetic improvement for P. salmonis resistance has been proposed as a feasible alternative for the control of this infectious disease in farmed fish. Genotyping by sequencing (GBS) strategies allow genotyping of hundreds of individuals with thousands of single nucleotide polymorphisms (SNPs), which can be used to perform genome wide association studies (GWAS) and predict genetic values using genome-wide information. We used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic architecture of resistance against P. salmonis in a farmed coho salmon population and to identify molecular markers associated with the trait. We also evaluated genomic selection (GS) models in order to determine the potential to accelerate the genetic improvement of this trait by means of using genome-wide molecular information. A total of 764 individuals from 33 full-sib families (17 highly resistant and 16 highly susceptible) were experimentally challenged against P. salmonis and their genotypes were assayed using ddRAD sequencing. A total of 9,389 SNPs markers were identified in the population. These markers were used to test genomic selection models and compare different GWAS methodologies for resistance measured as day of death (DD) and binary survival (BIN). Genomic selection models showed higher accuracies than the traditional pedigree-based best linear unbiased prediction (PBLUP) method, for both DD and BIN. The models showed an improvement of up to 95% and 155% respectively over PBLUP. One SNP related with B-cell development was identified as a potential functional candidate associated with resistance to P. salmonis defined as DD.

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.

Genomic Prediction Accuracy for Resistance Against Piscirickettsia salmonis in Farmed Rainbow Trout.

Salmonid rickettsial syndrome (SRS), caused by the intracellular bacterium Piscirickettsia salmonis, is one of the main diseases affecting rainbow trout (Oncorhynchus mykiss) farming. To accelerate genetic progress, genomic selection methods can be used as an effective approach to control the disease. The aims of this study were: (i) to compare the accuracy of estimated breeding values using pedigree-based best linear unbiased prediction (PBLUP) with genomic BLUP (GBLUP), single-step GBLUP (ssGBLUP), Bayes C, and Bayesian Lasso (LASSO); and (ii) to test the accuracy of genomic prediction and PBLUP using different marker densities (0.5, 3, 10, 20, and 27 K) for resistance against P. salmonis in rainbow trout. Phenotypes were recorded as number of days to death (DD) and binary survival (BS) from 2416 fish challenged with P. salmonis A total of 1934 fish were genotyped using a 57 K single-nucleotide polymorphism (SNP) array. All genomic prediction methods achieved higher accuracies than PBLUP. The relative increase in accuracy for different genomic models ranged from 28 to 41% for both DD and BS at 27 K SNP. Between different genomic models, the highest relative increase in accuracy was obtained with Bayes C (∼40%), where 3 K SNP was enough to achieve a similar accuracy to that of the 27 K SNP for both traits. For resistance against P. salmonis in rainbow trout, we showed that genomic predictions using GBLUP, ssGBLUP, Bayes C, and LASSO can increase accuracy compared with PBLUP. Moreover, it is possible to use relatively low-density SNP panels for genomic prediction without compromising accuracy predictions for resistance against P. salmonis in rainbow trout.

Genetic parameters in female reproductive traits of Nile tilapia (Oreochromis niloticus).

Genetic parameters for reproductive traits in female Nile tilapia were estimated in this study using Bayesian inference method. The data set presented information from 451 Nile tilapia females that were evaluated at two different places in Maringá - Paraná - Brazil (hapas of 1 and 10 m³) and at one location in Alfenas - Minas Gerais - Brazil. A spawning examination was conducted once a week from October 2012 to March 2013 for a total of 23 weeks of evaluation. Single-trait analyses for spawning success, multiple spawning, spawning frequency, and volume of eggs/female were performed by using the software MTGSAM Threshold. The heritability estimates were 0.14, 0.16, 0.53, and 0.08 for spawning success, multiple spawning, spawning frequency and volume of eggs/female, respectively, indicating it is possible to achieve a substantial genetic gain using these reproductive traits as selection criteria.