Evaluation of alternative methods for estimating the precision of REML-based estimates of variance components and heritability.

Residual Maximum Likelihood (REML) analysis is the most widely used method to estimate variance components and heritability. This method is based on large sample theory under the assumption that the parameter estimates are asymptotically multivariate normally distributed with covariance matrix given by the inverse of the information matrix. Hence, these sampling variances could be biased if the assumption of asymptotic approximation is incorrect, especially when the sample size is small. Though it is difficult to assess the impact of sample size, an alternative option is to generate a full distribution of the parameters to determine the uncertainty of estimates. In this study, we compared the REML estimates of variance components, heritability and sampling variances of body-weight (BW), body-depth (BD), and condition-factor (K) with those obtained from four sampling-based methods viz., parametric and nonparametric bootstrap, asymptotic sampling and Bayesian estimation. The aim was to understand if a sample size of order 1413 was sufficient to contain adequate information for a reliable asymptotic approximation. The REML solution was close to values obtained from different sampling-based methods indicating that the present sample size was sufficient to estimate reliable genetic variation in different traits with varying heritability. The variance and heritability estimated by a nonparametric bootstrap estimate based on randomization of family effects gave comparable results as evaluated by REML for different traits. Hence, the nonparametric bootstrap estimate can be effectively used to understand whether the sample size is large enough to contain sufficient information under likelihood estimation assumptions.

Genomic selection for white spot syndrome virus resistance in whiteleg shrimp boosts survival under an experimental challenge test.

White spot syndrome virus (WSSV) causes major worldwide losses in shrimp aquaculture. The development of resistant shrimp populations is an attractive option for management of the disease. However, heritability for WSSV resistance is generally low and genetic improvement by conventional selection has been slow. This study was designed to determine the power and accuracy of genomic selection to improve WSSV resistance in Litopenaeus vannamei. Shrimp were experimentally challenged with WSSV and resistance was evaluated as dead or alive (DOA) 23 days after infestation. All shrimp in the challenge test were genotyped for 18,643 single nucleotide polymorphisms. Breeding candidates (G0) were ranked on genomic breeding values for WSSV resistance. Two G1 populations were produced, one from G0 breeders with high and the other with low estimated breeding values. A third population was produced from "random" mating of parent stock. The average survival was 25% in the low, 38% in the random and 51% in the high-genomic breeding value groups. Genomic heritability for DOA (0.41 in G1) was high for this type of trait. The realised genetic gain and high heritability clearly demonstrates large potential for further genetic improvement of WSSV resistance in the evaluated L. vannamei population using genomic selection.

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.

Genomic predictions can accelerate selection for resistance against Piscirickettsia salmonis in Atlantic salmon (Salmo salar).

BACKGROUND: Salmon Rickettsial Syndrome (SRS) caused by Piscirickettsia salmonis is a major disease affecting the Chilean salmon industry. Genomic selection (GS) is a method wherein genome-wide markers and phenotype information of full-sibs are used to predict genomic EBV (GEBV) of selection candidates and is expected to have increased accuracy and response to selection over traditional pedigree based Best Linear Unbiased Prediction (PBLUP). Widely used GS methods such as genomic BLUP (GBLUP), SNPBLUP, Bayes C and Bayesian Lasso may perform differently with respect to accuracy of GEBV prediction. Our aim was to compare the accuracy, in terms of reliability of genome-enabled prediction, from different GS methods with PBLUP for resistance to SRS in an Atlantic salmon breeding program. Number of days to death (DAYS), binary survival status (STATUS) phenotypes, and 50 K SNP array genotypes were obtained from 2601 smolts challenged with P. salmonis. The reliability of different GS methods at different SNP densities with and without pedigree were compared to PBLUP using a five-fold cross validation scheme. RESULTS: Heritability estimated from GS methods was significantly higher than PBLUP. Pearson's correlation between predicted GEBV from PBLUP and GS models ranged from 0.79 to 0.91 and 0.79-0.95 for DAYS and STATUS, respectively. The relative increase in reliability from different GS methods for DAYS and STATUS with 50 K SNP ranged from 8 to 25% and 27-30%, respectively. All GS methods outperformed PBLUP at all marker densities. DAYS and STATUS showed superior reliability over PBLUP even at the lowest marker density of 3 K and 500 SNP, respectively. 20 K SNP showed close to maximal reliability for both traits with little improvement using higher densities. CONCLUSIONS: These results indicate that genomic predictions can accelerate genetic progress for SRS resistance in Atlantic salmon and implementation of this approach will contribute to the control of SRS in Chile. We recommend GBLUP for routine GS evaluation because this method is computationally faster and the results are very similar with other GS methods. The use of lower density SNP or the combination of low density SNP and an imputation strategy may help to reduce genotyping costs without compromising gain in reliability.

The use of genomic information increases the accuracy of breeding value predictions for sea louse (Caligus rogercresseyi) resistance in Atlantic salmon (Salmo salar).

Sea lice infestations caused by Caligus rogercresseyi are a main concern to the salmon farming industry due to associated economic losses. Resistance to this parasite was shown to have low to moderate genetic variation and its genetic architecture was suggested to be polygenic. The aim of this study was to compare accuracies of breeding value predictions obtained with pedigree-based best linear unbiased prediction (P-BLUP) methodology against different genomic prediction approaches: genomic BLUP (G-BLUP), Bayesian Lasso, and Bayes C. To achieve this, 2404 individuals from 118 families were measured for C. rogercresseyi count after a challenge and genotyped using 37 K single nucleotide polymorphisms. Accuracies were assessed using fivefold cross-validation and SNP densities of 0.5, 1, 5, 10, 25 and 37 K. Accuracy of genomic predictions increased with increasing SNP density and was higher than pedigree-based BLUP predictions by up to 22%. Both Bayesian and G-BLUP methods can predict breeding values with higher accuracies than pedigree-based BLUP, however, G-BLUP may be the preferred method because of reduced computation time and ease of implementation. A relatively low marker density (i.e. 10 K) is sufficient for maximal increase in accuracy when using G-BLUP or Bayesian methods for genomic prediction of C. rogercresseyi resistance in Atlantic salmon.