Genome-wide association study for milk production and conformation traits in Canadian Alpine and Saanen dairy goats.

Increasing the productivity of Canadian dairy goats is critical to the competitiveness of the sector; however, little is known about the underlying genetic architecture of economically important traits in these populations. Consequently, the objectives of this study were as follows: (1) to perform a single-step GWAS for milk production traits (milk, protein, and fat yields, and protein and fat percentages in first and later lactations) and conformation traits (body capacity, dairy character, feet and legs, fore udder, general appearance, rear udder, suspensory ligament, and teats) in the Canadian Alpine and Saanen breeds; and (2) to identify positional and functional candidate genes related to these traits. The data available for analysis included 305-d milk production records for 6,409 Alpine and 3,434 Saanen does in first lactation and 5,827 Alpine and 2,632 Saanen does in later lactations; as well as linear type conformation records for 5,158 Alpine and 2,342 Saanen does. Genotypes were available for 833 Alpine and 874 Saanen animals. Both single-breed and multiple-breed GWAS were performed using single-trait animal models. Positional and functional candidate genes were then identified in downstream analyses. The GWAS identified 189 unique SNP that were significant at the chromosomal level, corresponding to 271 unique positional candidate genes within 50 kb up- and downstream, across breeds and traits. This study provides evidence for the economic importance of several candidate genes (e.g., CSN1S1, CSN2, CSN1S2, CSN3, DGAT1, and ZNF16) in the Canadian Alpine and Saanen populations that have been previously reported in other dairy goat populations. Moreover, several novel positional and functional candidate genes (e.g., RPL8, DCK, and MOB1B) were also identified. Overall, the results of this study have provided greater insight into the genetic architecture of milk production and conformation traits in the Canadian Alpine and Saanen populations. Greater understanding of these traits will help to improve dairy goat breeding programs.

Single- and multiple-breed genomic evaluations for conformation traits in Canadian Alpine and Saanen dairy goats.

Conformation traits are functional traits known to affect longevity, production efficiency, and profitability of dairy goats. However, genetic progress for these traits is expected to be slower than for milk production traits due to the limited number of herds participating in type classification programs, and often lower heritability estimates. Genomic selection substantially accelerates the rate of genetic progress in many species and industries, especially for lowly heritable, difficult, or expensive to measure traits. Therefore, the main objectives of this study were (1) to evaluate the potential benefits of the implementation of single-step genomic evaluations for conformation traits in Canadian Alpine and Saanen dairy goats, and (2) to investigate the effect of the use of single- and multiple-breed training populations. The phenotypes used in this study were linear conformation scores, on a 1-to-9 scale, for 8 traits (i.e., body capacity, dairy character, fore udder, feet and legs, general appearance, rear udder, medial suspensory ligament, and teats) of 5,158 Alpine and 2,342 Saanen does. Genotypes were available for 833 Alpine and 874 Saanen animals. Averaged across all traits, the use of multiple-breed analyses increased validation accuracy for Saanen, and reduced bias of genomically enhanced breeding values (GEBV) for both Alpine and Saanen compared with single-breed analyses. Little benefit was observed from the use of GEBV relative to pedigree-based EBV in terms of validation accuracy and bias, possibly due to limitations in the validation design, but substantial gains of 0.14 to 0.21 (32-50%) were observed in the theoretical accuracy of validation animals when averaged across traits for single- and multiple-breed analyses. Across the whole genotyped population, average gains in theoretical accuracy for GEBV compared with EBV across all traits ranged from 0.15 to 0.17 (32-37%) for Alpine and 0.17 to 0.19 (40-41%) for Saanen, depending on the model used. The largest gains were observed for does without classification records (0.19-0.22 or 50-55%) and bucks without daughter classification records (0.20-0.27 or 57-82%), which have the least information contributing to their traditional EBV. The use of multiple-breed rather than single-breed models was most beneficial for the Saanen breed, which had fewer phenotypic records available for the analyses. These results suggest that the implementation of genomic selection could increase the accuracy of breeding values for conformation traits in Canadian dairy goats.

Single-step genomic evaluation of milk production traits in Canadian Alpine and Saanen dairy goats.

Genomic evaluations are routine in most plant and livestock breeding programs but are used infrequently in dairy goat breeding schemes. In this context, the purpose of this study was to investigate the use of the single-step genomic BLUP method for predicting genomic breeding values for milk production traits (milk, protein, and fat yields; protein and fat percentages) in Canadian Alpine and Saanen dairy goats. There were 6,409 and 12,236 Alpine records and 3,434 and 5,008 Saanen records for each trait in first and later lactations, respectively, and a total of 1,707 genotyped animals (833 Alpine and 874 Saanen). Two validation approaches were used, forward validation (i.e., animals born after 2013 with an average estimated breeding value accuracy from the full data set ≥0.50) and forward cross-validation (i.e., subsets of all animals included in the forward validation were used in successive replications). The forward cross-validation approach resulted in similar validation accuracies (0.55 to 0.66 versus 0.54 to 0.61) and biases (-0.01 to -0.07 versus -0.03 to 0.11) to the forward validation when averaged across traits. Additionally, both single and multiple-breed analyses were compared, and similar average accuracies and biases were observed across traits. However, there was a small gain in accuracy from the use of multiple-breed models for the Saanen breed. A small gain in validation accuracy for genomically enhanced estimated breeding values (GEBV) relative to pedigree-based estimated breeding values (EBV) was observed across traits for the Alpine breed, but not for the Saanen breed, possibly due to limitations in the validation design, heritability of the traits evaluated, and size of the training populations. Trait-specific gains in theoretical accuracy of GEBV relative to EBV for the validation animals ranged from 17 to 31% in Alpine and 35 to 55% in Saanen, using the cross-validation approach. The GEBV predicted from the full data set were 12 to 16% more accurate than EBV for genotyped animals, but no gains were observed for nongenotyped animals. The largest gains were found for does without lactation records (35-41%) and bucks without daughter records (46-54%), and consequently, the implementation of genomic selection in the Canadian dairy goat population would be expected to increase selection accuracy for young breeding candidates. Overall, this study represents the first step toward implementation of genomic selection in Canadian dairy goat populations.

The value of incorporating carcass trait phenotypes in terminal sire selection indexes to improve carcass weight and quality of heavy lambs.

Genetic selection for carcass traits is paramount to maximize the profitability and long-term sustainability of any meat-producing livestock species. The main objectives of this research were to evaluate the efficiency of indicator traits for the genetic improvement of lamb carcass traits and to determine the value of including carcass traits into terminal sire selection indexes for the Canadian sheep industry. The carcass traits included hot carcass weight (HCW), fat depth at the GR site (FATGR) and average carcass conformation score (AVGCONF), and were measured on heavy lambs (slaughter age less than 365 days and HCW greater than 16.3 kg) in commercial abattoirs. Growth traits were found to be moderately efficient indicator traits for the genetic improvement of HCW but selection on ultrasound traits was necessary to substantially improve the carcass quality traits (FATGR and AVGCONF). Economic selection indexes were designed by adding various combinations of carcass traits into the Canadian Sheep Genetic Evaluation System terminal indexes. Records measured on individuals and progeny were assumed to be the sources of information for live animal and carcass traits, respectively. The changes in index accuracy, efficiency and expected correlated response were used to assess the value of their inclusion. HCW was found to have a large economic value, and its inclusion into terminal selection indexes was expected to substantially increase their accuracy (0.08-0.12 points) and efficiency (20%-30%). However, further including FATGR (measured 110 mm from the carcass midline over the 12th rib) and AVGCONF had little impact on the accuracy (≤0.03) and efficiency (1%-7%) of the proposed indexes. Thus, the inclusion of carcass traits into the existing terminal selection indexes could be beneficial for the genetic improvement of HCW, but further research is needed to determine optimal methods of increasing carcass fatness and muscularity.

A genetic evaluation of growth, ultrasound, and carcass traits at alternative slaughter endpoints in crossbred heavy lambs.

Genetic parameters were estimated for growth, ultrasound, and carcass traits in a Canadian crossbred heavy lamb population. Traits analyzed included birth, weaning, post-weaning, and ultrasound scanning weights; pre- and post-weaning average daily gain; ultrasonically measured eye muscle and fat depths; hot carcass weight; fat depth at the GR site (110 mm from the midline on the 12th rib); carcass conformation scores; saleable meat yield; price grid value; and total carcass value. The impact of three alternative slaughter endpoints (slaughter age, carcass weight, and carcass fatness) on genetic parameter estimates was also evaluated. In general, carcass traits were found to be moderately heritable, with heritability estimates ranging from 0.17 ± 0.02 for hot carcass weight at a constant slaughter age to 0.34 ± 0.02 for average carcass conformation score at a constant carcass weight. Heritability estimates were similar when observations were adjusted to alternative slaughter endpoints, but for some traits, phenotypic variance and genetic correlation estimates differed. Genetic correlations between carcass traits and growth and ultrasound traits were typically favorable. Ultrasonically measured eye muscle depth and fat depth were found to be moderately to strongly positively correlated with hot carcass weight (0.33 ± 0.15 to 0.71 ± 0.19) and fat depth at the GR site (0.38 ± 0.14 to 0.74 ± 0.12), respectively, reaffirming the usefulness of selection on ultrasound traits to improve carcass yield and quality. Genetic correlations among carcass traits were generally favorable, with the exception of moderate unfavorable positive genetic correlations between fat depth at the GR site and primal cut carcass conformation scores (0.31 ± 0.05 to 0.60 ± 0.05). Overall, the results of this research suggest that there is potential to improve carcass yield and quality through genetic selection and provides the population-specific genetic parameter estimates needed for the genetic evaluation of carcass traits in the Canadian sheep population. Nevertheless, the optimal endpoint for carcass trait genetic evaluations will need to be further investigated, considering both the current findings and additional information on production practices in the industry.