Comparing genomic prediction accuracy from purebred, crossbred and combined purebred and crossbred reference populations in sheep

Combined genomic evaluation of Merino and Dohne Merino Australian sheep populations

BACKGROUND

The Dohne Merino sheep was introduced to Australia from South Africa in the 1990s. It was primarily used in crosses with the Merino breed sheep to improve on attributes such as reproduction and carcass composition. Since then, this breed has continued to expand in Australia but the number of genotyped and phenotyped purebred individuals remains low, calling into question the accuracy of genomic selection. The Australian Merino, on the other hand, has a substantial reference population in a separate genomic evaluation (MERINOSELECT). Combining these resources could fast track the impact of genomic selection on the smaller breed, but the efficacy of this needs to be investigated. This study was based on a dataset of 53,663 genotypes and more than 2 million phenotypes. Its main objectives were (1) to characterize the genetic structure of Merino and Dohne Merino breeds, (2) to investigate the utility of combining their evaluations in terms of quality of predictions, and (3) to compare several methods of genetic grouping. We used the 'LR-method' (Linear Regression) for these assessments.

RESULTS

We found very low Fst values (below 0.048) between the different Merino lines and Dohne breed considered in our study, indicating very low genetic differentiation. Principal component analysis revealed three distinct groups, identified as purebred Merino, purebred Dohne, and crossbred animals. Considering the whole population in the reference led to the best quality of predictions and the largest increase in accuracy (from 'LR-method') from pedigree to genomic-based evaluations: 0.18, 0.14 and 0.16 for yearling fibre diameter (YFD), yearling greasy fleece weight (YGFW) and yearling liveweight (YWT), respectively. Combined genomic evaluations showed higher accuracies than the evaluation based on the Dohne reference only (accuracies increased by 0.16, 0.06 and 0.07 for YFD, YGFW, and YWT, respectively). For the combined genomic evaluations, metafounder models were more accurate than Unknown Parent Groups models (accuracies increased by 0.04, 0.04 and 0.06 for YFD, YGFW and YWT, respectively).

CONCLUSIONS

We found promising results for the future transition of the Dohne breed from pedigree to genomic selection. A combined genomic evaluation, with the MERINOSELECT evaluation in addition to using metafounders, is expected to enhance the quality of genomic predictions for the Dohne Merino breed.

Assessment of genomic breeding values and their accuracies for carcass traits in Jeju Black cattle using whole-genome SNP chip panels

The objective of the present study was to evaluate the breeding value and accuracy of genomic estimated breeding values (GEBVs) of carcass traits in Jeju Black cattle (JBC) using Hanwoo steers and JBC as a reference population using the single-trait animal model. Our research included genotype and phenotype information on 19,154 Hanwoo steers with 1097 JBC acting as the reference population. Likewise, the test population consisted of 418 genotyped JBC individuals with no phenotypic records for those carcass traits. For estimating the accuracy of GEBV, we divided the entire population into three groups. Hanwoo and JBC make up the first group; Hanwoo and JBC, who has both the genotype and phenotypic records, are referred to as the reference (training) population, and JBC, who lacks phenotypic information is referred to as the test (validation) population. The second group consists of the JBC (without phenotype) as the test population and Hanwoo as a reference population with phenotype and genotypic data. The only JBCs in the third group are those who have genotypic and phenotypic data on them as a reference population but no phenotypic data on them as a test population. The single-trait animal model was used in all three groups for statistical purposes. The reference populations estimated heritabilities for carcass weight (CWT), eye muscle area (EMA), backfat thickness (BF), and marbling score (MS) as 0.30, 0.26, 0.26, and 0.34 for the Hanwoo steer and 0.42, 0.27, 0.26, and 0.48 for JBC. The average accuracy for carcass traits in Group 1 was 0.80 for the Hanwoo and JBC reference population compared with 0.73 for the JBC test population. Although the average accuracy for carcass traits in Group 2 was 0.80, it was 0.80 for the Hanwoo reference population and only 0.56 for the JBC test population. The average accuracy for the JBC reference and test populations was 0.68 and 0.50, respectively, when they were included in the accuracy comparison without the Hanwoo reference population. Groups 1 and 2 used Hanwoo as reference population, which led to a better average accuracy; however, Group 3 only used the JBC reference and test population, which led to a lower average accuracy. This might be due to the fact that Group 3 used a smaller reference size than the group that came before it and that the genetic makeup of the Hanwoo and JBC breeds differed. The GEBV accuracy for MS was higher than that of other traits across all three analysis groups, followed by CWT, EMA, and BF, which may be partially explained by the MS traits' higher heritability. This study suggests that in order to achieve more accuracy, a large reference population particular to a breed should be established. Therefore, to increase the accuracy of GEBV prediction and the genetic benefit from genomic selection in JBC, individual reference breeds, and large populations are required.

Effect of subdivision of the Lacaune dairy sheep breed on the accuracy of genomic prediction

Genomic selection was deployed in Lacaune dairy breed in 2015. Lacaune population split in 1972 into 2 breeding companies with associated flocks, and there have been very few exchanges of animals between the subpopulations, leading to divergence of the 2 subpopulations. In spite of that, there is a joint genomic prediction. The objective of this study is to understand how this structuring affects prediction accuracy. We analyzed all the data available from Lacaune breeding program for milk yield: around 6 million phenotypes, 2 million animals in the pedigree and more than 29,000 genotyped animals, including 3,434 and 2,868 AI rams for each company. To consider missing pedigree, we set up genetic groups using the theory of metafounders. First, we studied the pedigree and genomic structures of the 2 subpopulations calculating Fst, evolution of average pedigree relationships across time and principal components analysis of genomic relationships. In a second part, we compared the reliability between different scenarios: an evaluation with a single reference population (Alone), an evaluation with a joint reference population (Together) and an evaluation of one subpopulation based on the reference population of the other group (Indirect). The low Fst value (0.02) reveals that the 2 subpopulations are still genetically close. Nevertheless, a low and constant average relationship between the animals of the 2 subpopulations confirms the absence of recent connections between them. We can see with principal component analysis results that even if they are close, they diverge over time. Finally, we observe small gains in accuracy of Together versus Alone, in spite of whereas doubling the reference population size in Together. These gains vary across years and subpopulations: less than 0.08 (0.46 to 0.54; ratio of accuracy for the partial and whole evaluations-corresponding to the greatest change in this ratio for breeding company 1, observed for the cohort 2016) for one subpopulation and between 0.03 (0.55 to 0.58) and 0.17 (0.48 to 0.65) for the other. To conclude, the 2 subpopulations remain close enough genetically so that their combined evaluation is advantageous, even if only slightly.

Effects of the Breeding Strategy Beef-on-Dairy at Animal, Farm and Sector Levels

The decline in farm revenue due to volatile milk prices has led to an increase in the use of beef semen in dairy herds. While this strategy ("Beef-on-dairy" (BoD)) can have economic benefits, it can also lead to unintended consequences affecting animal welfare. Semen sale trends from breeding organizations depict increasing sales of beef semen across the globe. Calves born from such breeding strategies can perform better when compared to purebred dairy calves, especially in terms of meat quality and growth traits. The Beef-on-dairy strategy can lead to unintentional negative impacts including an increase in gestation length, and increased dystocia and stillbirth rates. Studies in this regard have found the highest gestation length for Limousin crossbred calves followed by calves from the Angus breed. This increase in gestation length can lead to economic losses ranging from 3 to 5 US$ per animal for each additional day. In terms of the growth performance of crossbred animals, literature studies are inconclusive due to the vast differences in farming structure across the regions. But almost all the studies agree regarding improvement in the meat quality in terms of color, fiber type, and intra-muscular fat content for crossbred animals. Utilization of genomic selection, and development of specialized Beef-on-dairy indexes for the sires, can be a viable strategy to make selection easier for the farmers.

Clonal diploid and autopolyploid breeding strategies to harness heterosis: insights from stochastic simulation

KEY MESSAGE

Reciprocal recurrent selection sometimes increases genetic gain per unit cost in clonal diploids with heterosis due to dominance, but it typically does not benefit autopolyploids. Breeding can change the dominance as well as additive genetic value of populations, thus utilizing heterosis. A common hybrid breeding strategy is reciprocal recurrent selection (RRS), in which parents of hybrids are typically recycled within pools based on general combining ability. However, the relative performances of RRS and other breeding strategies have not been thoroughly compared. RRS can have relatively increased costs and longer cycle lengths, but these are sometimes outweighed by its ability to harness heterosis due to dominance. Here, we used stochastic simulation to compare genetic gain per unit cost of RRS, terminal crossing, recurrent selection on breeding value, and recurrent selection on cross performance considering different amounts of population heterosis due to dominance, relative cycle lengths, time horizons, estimation methods, selection intensities, and ploidy levels. In diploids with phenotypic selection at high intensity, whether RRS was the optimal breeding strategy depended on the initial population heterosis. However, in diploids with rapid-cycling genomic selection at high intensity, RRS was the optimal breeding strategy after 50 years over almost all amounts of initial population heterosis under the study assumptions. Diploid RRS required more population heterosis to outperform other strategies as its relative cycle length increased and as selection intensity and time horizon decreased. The optimal strategy depended on selection intensity, a proxy for inbreeding rate. Use of diploid fully inbred parents vs. outbred parents with RRS typically did not affect genetic gain. In autopolyploids, RRS typically did not outperform one-pool strategies regardless of the initial population heterosis.

Progress and opportunities through use of genomics in animal production

The rearing of farmed animals is a vital component of global food production systems, but its impact on the environment, human health, animal welfare, and biodiversity is being increasingly challenged. Developments in genetic and genomic technologies have had a key role in improving the productivity of farmed animals for decades. Advances in genome sequencing, annotation, and editing offer a means not only to continue that trend, but also, when combined with advanced data collection, analytics, cloud computing, appropriate infrastructure, and regulation, to take precision livestock farming (PLF) and conservation to an advanced level. Such an approach could generate substantial additional benefits in terms of reducing use of resources, health treatments, and environmental impact, while also improving animal health and welfare.

Genotyping and phenotyping strategies for genetic improvement of meat quality and carcass composition in swine

Background

Meat quality and composition traits have become valuable in modern pork production; however, genetic improvement has been slow due to high phenotyping costs. Combining genomic information with multi-trait indirect selection based on cheaper indicator traits is an alternative for continued cost-effective genetic improvement.

Methods

Data from an ongoing breeding program were used in this study. Phenotypic and genomic information was collected on three-way crossbred and purebred Duroc animals belonging to 28 half-sib families. We applied different methods to assess the value of using purebred and crossbred information (both genomic and phenotypic) to predict expensive-to-record traits measured on crossbred individuals. Estimation of multi-trait variance components set the basis for comparing the different scenarios, together with a fourfold cross-validation approach to validate the phenotyping schemes under four genotyping strategies.

Results

The benefit of including genomic information for multi-trait prediction depended on the breeding goal trait, the indicator traits included, and the source of genomic information. While some traits benefitted significantly from genotyping crossbreds (e.g., loin intramuscular fat content, backfat depth, and belly weight), multi-trait prediction was advantageous for some traits even in the absence of genomic information (e.g., loin muscle weight, subjective color, and subjective firmness).

Conclusions

Our results show the value of using different sources of phenotypic and genomic information. For most of the traits studied, including crossbred genomic information was more beneficial than performing multi-trait prediction. Thus, we recommend including crossbred individuals in the reference population when these are phenotyped for the breeding objective.

Genomic evaluation with multibreed and crossbred data

•

We found low accuracy of genomic evaluation of crossbreds based on purebred data.

•

We found higher accuracy for crossbreds based on crossbred data.

•

Use of putative sequence variants had a small impact on genomic accuracy.

Several types of multibreed genomic evaluation are in use. These include evaluation of crossbreds based on purebred SNP effects, joint evaluation of all purebreds and crossbreds with a single additive effect, and treating each purebred and crossbred group as a separate trait. Additionally, putative quantitative trait nucleotides can be exploited to increase the accuracy of prediction. Existing studies indicate that the prediction of crossbreds based on purebred data has low accuracy, that a joint evaluation can potentially provide accurate evaluations for crossbreds but could lower accuracy for purebreds compared with single-breed evaluations, and that the use of putative quantitative trait nucleotides only marginally increases the accuracy. One hypothesis is that genomic selection is based on estimation of clusters of independent chromosome segments. Subsequently, predicting a particular group type would require a reference population of the same type, and crosses with same breed percentage but different type (F₁ vs. F₂) would, at best, use separate reference populations. The genomic selection of multibreed population is still an active research topic.

Genomic prediction in a numerically small breed population using prioritized genetic markers from whole-genome sequence data

The objective of this study was to investigate the accuracy of genomic prediction of body weight and eating quality traits in a numerically small sheep population (Dorper sheep). Prediction was based on a large multi-breed/admixed reference population and using (a) 50k or 500k single nucleotide polymorphism (SNP) genotypes, (b) imputed whole-genome sequencing data (~31 million), (c) selected SNPs from whole genome sequence data and (d) 50k SNP genotypes plus selected SNPs from whole-genome sequence data. Furthermore, the impact of using a breed-adjusted genomic relationship matrix on accuracy of genomic breeding value was assessed. The selection of genetic variants was based on an association study performed on imputed whole-genome sequence data in an independent population, which was chosen either randomly from the base population or according to higher genetic proximity to the target population. Genomic prediction was based on genomic best linear unbiased prediction (GBLUP), and the accuracy of genomic prediction was assessed according to the correlation between genomic breeding value and corrected phenotypes divided by the square root of trait heritability. The accuracy of genomic prediction was between 0.20 and 0.30 across different traits based on common 50k SNP genotypes, which improved on average by 0.06 (absolute value) on average based on using prioritized genetic markers from whole-genome sequence data. Using prioritized genetic markers from a genetically more related GWAS population resulted in slightly higher prediction accuracy (0.02 absolute value) compared to genetic markers derived from a random GWAS population. Using high-density SNP genotypes or imputed whole-genome sequence data in GBLUP showed almost no improvement in genomic prediction accuracy however, accounting for different marker allele frequencies in reference population according to a breed-adjusted GRM resulted to on average 0.024 (absolute value) increase in accuracy of genomic prediction.

Within-breed selection is sufficient to improve terminal crossbred beef marbling: a review of reciprocal recurrent genomic selection

Reciprocal recurrent selection is the selection of purebreds for crossbred performance and takes advantage of additive and non-additive variance by using pedigreed progeny performance records. Developed in maize, the adoption of this approach in livestock breeding has been limited to the pork and poultry industries; genomic selection may facilitate its extension into the beef industry by replacing pedigree. The literature regarding the relative importance of additive versus non-additive variance and reciprocal recurrent genomic selection models was reviewed. The potential for using reciprocal recurrent genomic selection in a terminal Wagyu × Angus cross scenario was examined. Non-additive variance is more important for fitness traits and accounts for a small proportion of variance related to production traits such as marbling. In general, reciprocal recurrent selection was not significantly better at improving performance of crossbreds than was traditional selection within parental breeds using only additive variance in the studies examined. Simulation studies showed benefits of including dominance or breed-specific allele effects in prediction models but advantages were small as more realistic simulations were examined. On the basis of the evidence, it is likely that in a terminal two-way cross-beef scenario utilising Wagyu sires and Angus dams, where selection emphasis is on marbling, selection of purebreds on the basis of additive variance will allow substantial progress to be realised.

Breed-adjusted genomic relationship matrices as a method to account for population stratification in multibreed populations of tropically adapted beef heifers

Context Beef cattle breeds in Australia can broadly be broken up into two subspecies, namely, Bos indicus and Bos taurus. Due to the time since divergence between the subspecies, it is likely that mutations affecting quantitative traits have developed independently in each. Aims We hypothesise that this will affect the prediction accuracy of genomic selection of admixed and composite populations that include both ancestral subspecies. Our study investigates methods to quantify population stratification in a multibreed population of tropically adapted heifers, with the aim of improving prediction accuracy of genomic selection for reproductive maturity score. Methods We used genotypes and reproductive maturity phenotypes from 3695 tropically adapted heifers from three purebred populations, namely, Brahman, Santa Gertrudis and Droughtmaster. Two of these breeds, Santa Gertrudis and Droughtmaster, are stabilised composites of varying B. indicus × B. taurus ancestry, and the third breed, Brahman, has predominately B. indicus ancestry. Genotypes were imputed to three marker-panel densities and population stratification was accounted for in genomic relationship matrices by using breed-specific allele frequencies when calculating the genomic relationships among animals. Prediction accuracy and bias were determined using a five-fold cross validation of randomly selected multibreed cohorts. Key Results Our results showed that the use of breed-adjusted genomic relationship matrices did not improve either prediction accuracy or bias for a lowly heritable trait such as reproductive maturity score. However, using breed-adjusted genomic relationship matrices allowed the capture of a higher proportion of additive genetic effects when estimating variance components. Conclusions These findings suggest that, despite seeing no improvement in prediction accuracy, it may still be beneficial to use breed-adjusted genomic relationship matrices in multibreed populations to improve the estimation of variance components. Implications As such, genomic evaluations using breed-adjusted genomic relationship matrices may be beneficial in multibreed populations.

Comparison of prediction accuracy for genomic estimated breeding value using the reference pig population of single-breed and admixed-breed

This study was performed to increase the accuracy of genomic estimated breeding value (GEBV) predictions for domestic pigs using single-breed and admixed reference populations (single-breed of Berkshire pigs [BS] with cross breed of Korean native pigs and Landrace pigs [CB]). The principal component analysis (PCA), linkage disequilibrium (LD), and genome-wide association study (GWAS) were performed to analyze the population structure prior to genomic prediction. Reference and test population data sets were randomly sampled 10 times each and precision accuracy was analyzed according to the size of the reference population (100, 200, 300, or 400 animals). For the BS population, prediction accuracy was higher for all economically important traits with larger reference population size. Prediction accuracy was ranged from −0.05 to 0.003, for all traits except carcass weight (CWT), when CB was used as the reference population and BS as the test. The accuracy of CB for backfat thickness (BF) and shear force (SF) using admixed population as reference increased with reference population size, while the results for CWT and muscle pH at 24 hours after slaughter (pH) were equivocal with respect to the relationship between accuracy and reference population size, although overall accuracy was similar to that using the BS as the reference.

Genomic prediction based on selected variants from imputed whole-genome sequence data in Australian sheep populations

Background

Whole-genome sequence (WGS) data could contain information on genetic variants at or in high linkage disequilibrium with causative mutations that underlie the genetic variation of polygenic traits. Thus far, genomic prediction accuracy has shown limited increase when using such information in dairy cattle studies, in which one or few breeds with limited diversity predominate. The objective of our study was to evaluate the accuracy of genomic prediction in a multi-breed Australian sheep population of relatively less related target individuals, when using information on imputed WGS genotypes.

Methods

Between 9626 and 26,657 animals with phenotypes were available for nine economically important sheep production traits and all had WGS imputed genotypes. About 30% of the data were used to discover predictive single nucleotide polymorphism (SNPs) based on a genome-wide association study (GWAS) and the remaining data were used for training and validation of genomic prediction. Prediction accuracy using selected variants from imputed sequence data was compared to that using a standard array of 50k SNP genotypes, thereby comparing genomic best linear prediction (GBLUP) and Bayesian methods (BayesR/BayesRC). Accuracy of genomic prediction was evaluated in two independent populations that were each lowly related to the training set, one being purebred Merino and the other crossbred Border Leicester x Merino sheep.

Results

A substantial improvement in prediction accuracy was observed when selected sequence variants were fitted alongside 50k genotypes as a separate variance component in GBLUP (2GBLUP) or in Bayesian analysis as a separate category of SNPs (BayesRC). From an average accuracy of 0.27 in both validation sets for the 50k array, the average absolute increase in accuracy across traits with 2GBLUP was 0.083 and 0.073 for purebred and crossbred animals, respectively, whereas with BayesRC it was 0.102 and 0.087. The average gain in accuracy was smaller when selected sequence variants were treated in the same category as 50k SNPs. Very little improvement over 50k prediction was observed when using all WGS variants.

Conclusions

Accuracy of genomic prediction in diverse sheep populations increased substantially by using variants selected from whole-genome sequence data based on an independent multi-breed GWAS, when compared to genomic prediction using standard 50K genotypes.

Validation of genomic predictions for body weight in broilers using crossbred information and considering breed-of-origin of alleles

Background

Pig and poultry breeding programs aim at improving crossbred (CB) performance. Selection response may be suboptimal if only purebred (PB) performance is used to compute genomic estimated breeding values (GEBV) because the genetic correlation between PB and CB performance (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc) is often lower than 1. Thus, it may be beneficial to use information on both PB and CB performance. In addition, the accuracy of GEBV of PB animals for CB performance may improve when the breed-of-origin of alleles (BOA) is considered in the genomic relationship matrix (GRM). Thus, our aim was to compare scenarios where GEBV are computed and validated by using (1) either CB offspring averages or individual CB records for validation, (2) either a PB or CB reference population, and (3) a GRM that either accounts for or ignores BOA in the CB individuals. For this purpose, we used data on body weight measured at around 7 (BW7) or 35 (BW35) days in PB and CB broiler chickens and evaluated the accuracy of GEBV based on the correlation GEBV with phenotypes in the validation population (validation correlation).

Results

With validation on CB offspring averages, the validation correlation of GEBV of PB animals for CB performance was lower with a CB reference population than with a PB reference population for BW35 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc = 0.96), and about equal for BW7 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc = 0.80) when BOA was ignored. However, with validation on individual CB records, the validation correlation was higher with a CB reference population for both traits. The use of a GRM that took BOA into account increased the validation correlation for BW7 but reduced it for BW35.

Conclusions

We argue that the benefit of using a CB reference population for genomic prediction of PB animals for CB performance should be assessed either by validation on CB offspring averages, or by validation on individual CB records while using a GRM that accounts for BOA in the CB individuals. With this recommendation in mind, our results show that the accuracy of GEBV of PB animals for CB performance was equal to or higher with a CB reference population than with a PB reference population for a trait with an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc of 0.8, but lower for a trait with an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc of 0.96. In addition, taking BOA into account was beneficial for a trait with an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc of 0.8 but not for a trait with an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc of 0.96.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0481-7) contains supplementary material, which is available to authorized users.

Genomic selection on breeding time in a wild bird population

Artificial selection experiments are a powerful tool in evolutionary biology. Selecting individuals based on multimarker genotypes (genomic selection) has several advantages over phenotype-based selection but has, so far, seen very limited use outside animal and plant breeding. Genomic selection depends on the markers tagging the causal loci that underlie the selected trait. Because the number of necessary markers depends, among other factors, on effective population size, genomic selection may be in practice not feasible in wild populations as most wild populations have much higher effective population sizes than domesticated populations. However, the current possibilities of cost-effective high-throughput genotyping could overcome this limitation and thereby make it possible to apply genomic selection also in wild populations. Using a unique dataset of about 2000 wild great tits (Parus major), a small passerine bird, genotyped on a 650 k SNP chip we calculated genomic breeding values for egg-laying date using the so-called GBLUP approach. In this approach, the pedigree-based relatedness matrix of an "animal model," a special form of the mixed model, is replaced by a marker-based relatedness matrix. Using the marker-based relatedness matrix, the model seemed better able to disentangle genetic and permanent environmental effects. We calculated the accuracy of genomic breeding values by correlating them to the phenotypes of individuals whose phenotypes were excluded from the analysis when estimating the genomic breeding values. The obtained accuracy was about 0.20, with very little effect of the used genomic relatedness estimator but a strong effect of the number of SNPs. The obtained accuracy is lower than typically seen in domesticated species but considerable for a trait with low heritability (∼0.2) as avian breeding time. Our results show that genomic selection is possible also in wild populations with potentially many applications, which we discuss here.

Using a very low-density SNP panel for genomic selection in a breeding program for sheep

Background

Building an efficient reference population for genomic selection is an issue when the recorded population is small and phenotypes are poorly informed, which is often the case in sheep breeding programs. Using stochastic simulation, we evaluated a genomic design based on a reference population with medium-density genotypes [around 45 K single nucleotide polymorphisms (SNPs)] of dams that were imputed from very low-density genotypes (≤ 1000 SNPs).

Methods

A population under selection for a maternal trait was simulated using real genotypes. Genetic gains realized from classical selection and genomic selection designs were compared. Genomic selection scenarios that differed in reference population structure (whether or not dams were included in the reference) and genotype quality (medium-density or imputed to medium-density from very low-density) were evaluated.

Results

The genomic design increased genetic gain by 26% when the reference population was based on sire medium-density genotypes and by 54% when the reference population included both sire and dam medium-density genotypes. When medium-density genotypes of male candidates and dams were replaced by imputed genotypes from very low-density SNP genotypes (1000 SNPs), the increase in gain was 22% for the sire reference population and 42% for the sire and dam reference population. The rate of increase in inbreeding was lower (from − 20 to − 34%) for the genomic design than for the classical design regardless of the genomic scenario.

Conclusions

We show that very low-density genotypes of male candidates and dams combined with an imputation process result in a substantial increase in genetic gain for small sheep breeding programs.

Genomic evaluation for a three-way crossbreeding system considering breed-of-origin of alleles

BACKGROUND

Genomic prediction of purebred animals for crossbred performance can be based on a model that estimates effects of single nucleotide polymorphisms (SNPs) in purebreds on crossbred performance. For crossbred performance, SNP effects might be breed-specific due to differences between breeds in allele frequencies and linkage disequilibrium patterns between SNPs and quantitative trait loci. Accurately tracing the breed-of-origin of alleles (BOA) in three-way crosses is possible with a recently developed procedure called BOA. A model that accounts for breed-specific SNP effects (BOA model), has never been tested empirically on a three-way crossbreeding scheme. Therefore, the objectives of this study were to evaluate the estimates of variance components and the predictive accuracy of the BOA model compared to models in which SNP effects for crossbred performance were assumed to be the same across breeds, using either breed-specific allele frequencies ([Formula: see text] model) or allele frequencies averaged across breeds ([Formula: see text] model). In this study, we used data from purebred and three-way crossbred pigs on average daily gain (ADG), back fat thickness (BF), and loin depth (LD).

RESULTS

Estimates of variance components for crossbred performance from the BOA model were mostly similar to estimates from models [Formula: see text] and [Formula: see text]. Heritabilities for crossbred performance ranged from 0.24 to 0.46 between traits. Genetic correlations between purebred and crossbred performance ([Formula: see text]) across breeds ranged from 0.30 to 0.62 for ADG and from 0.53 to 0.74 for BF and LD. For ADG, prediction accuracies of the BOA model were higher than those of the [Formula: see text] and [Formula: see text] models, with significantly higher accuracies only for one maternal breed. For BF and LD, prediction accuracies of models [Formula: see text] and [Formula: see text] were higher than those of the BOA model, with no significant differences. Across all traits, models [Formula: see text] and [Formula: see text] yielded similar predictions.

CONCLUSIONS

The BOA model yielded a higher prediction accuracy for ADG in one maternal breed, which had the lowest [Formula: see text] (0.30). Using the BOA model was especially relevant for traits with a low [Formula: see text]. In all other cases, the use of crossbred information in models [Formula: see text] and [Formula: see text], does not jeopardize predictions and these models are more easily implemented than the BOA model.

Multiple-trait QTL mapping and genomic prediction for wool traits in sheep

BACKGROUND

The application of genomic selection to sheep breeding could lead to substantial increases in profitability of wool production due to the availability of accurate breeding values from single nucleotide polymorphism (SNP) data. Several key traits determine the value of wool and influence a sheep's susceptibility to fleece rot and fly strike. Our aim was to predict genomic estimated breeding values (GEBV) and to compare three methods of combining information across traits to map polymorphisms that affect these traits.

METHODS

GEBV for 5726 Merino and Merino crossbred sheep were calculated using BayesR and genomic best linear unbiased prediction (GBLUP) with real and imputed 510,174 SNPs for 22 traits (at yearling and adult ages) including wool production and quality, and breech conformation traits that are associated with susceptibility to fly strike. Accuracies of these GEBV were assessed using fivefold cross-validation. We also devised and compared three approximate multi-trait analyses to map pleiotropic quantitative trait loci (QTL): a multi-trait genome-wide association study and two multi-trait methods that use the output from BayesR analyses. One BayesR method used local GEBV for each trait, while the other used the posterior probabilities that a SNP had an effect on each trait.

RESULTS

BayesR and GBLUP resulted in similar average GEBV accuracies across traits (~0.22). BayesR accuracies were highest for wool yield and fibre diameter (>0.40) and lowest for skin quality and dag score (<0.10). Generally, accuracy was higher for traits with larger reference populations and higher heritability. In total, the three multi-trait analyses identified 206 putative QTL, of which 20 were common to the three analyses. The two BayesR multi-trait approaches mapped QTL in a more defined manner than the multi-trait GWAS. We identified genes with known effects on hair growth (i.e. FGF5, STAT3, KRT86, and ALX4) near SNPs with pleiotropic effects on wool traits.

CONCLUSIONS

The mean accuracy of genomic prediction across wool traits was around 0.22. The three multi-trait analyses identified 206 putative QTL across the ovine genome. Detailed phenotypic information helped to identify likely candidate genes.

Prediction of genomic breeding values for growth, carcass and meat quality traits in a multi-breed sheep population using a HD SNP chip

Background

New Zealand has some unique Terminal Sire composite sheep breeds, which were developed in the last three decades to meet commercial needs. These composite breeds were developed based on crossing various Terminal Sire and Maternal breeds and, therefore, present high genetic diversity compared to other sheep breeds. Their breeding programs are focused on improving carcass and meat quality traits. There is an interest from the industry to implement genomic selection in this population to increase the rates of genetic gain. Therefore, the main objectives of this study were to determine the accuracy of predicted genomic breeding values for various growth, carcass and meat quality traits using a HD SNP chip and to evaluate alternative genomic relationship matrices, validation designs and genomic prediction scenarios. A large multi-breed population (n = 14,845) was genotyped with the HD SNP chip (600 K) and phenotypes were collected for a variety of traits.

Results

The average observed accuracies (± SD) for traits measured in the live animal, carcass, and, meat quality traits ranged from 0.18 ± 0.07 to 0.33 ± 0.10, 0.28 ± 0.09 to 0.55 ± 0.05 and 0.21 ± 0.07 to 0.36 ± 0.08, respectively, depending on the scenario/method used in the genomic predictions. When accounting for population stratification by adjusting for 2, 4 or 6 principal components (PCs) the observed accuracies of molecular breeding values (mBVs) decreased or kept constant for all traits. The mBVs observed accuracies when fitting both G and A matrices were similar to fitting only G matrix. The lowest accuracies were observed for k-means cross-validation and forward validation performed within each k-means cluster.

Conclusions

The accuracies observed in this study support the feasibility of genomic selection for growth, carcass and meat quality traits in New Zealand Terminal Sire breeds using the Ovine HD SNP chip. There was a clear advantage on using a mixed training population instead of performing analyzes per genomic clusters. In order to perform genomic predictions per breed group, genotyping more animals is recommended to increase the size of the training population within each group and the genetic relationship between training and validation populations. The different scenarios evaluated in this study will help geneticists and breeders to make wiser decisions in their breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0476-8) contains supplementary material, which is available to authorized users.

Crossbreed evaluations in single-step genomic best linear unbiased predictor using adjusted realized relationship matrices

Combining purebreed and crossbreed information is beneficial for genetic evaluation of some livestock species. Genetic evaluations can use relationships based on genomic information, relying on allele frequencies that are breed specific. Single-step genomic BLUP (ssGBLUP) does not account for different allele frequencies, which could limit the genetic gain in crossbreed evaluations. In this study, we tested the performance of different breed-specific genomic relationship matrices () in ssGBLUP for crossbreed evaluations; we also tested the importance of genotyping crossbred animals. Genotypes were available for purebreeds (AA and BB) and crossbreeds (F) in simulated and real pig populations. The number of genotyped animals was, on average, 4,315 for the simulated population and 15,798 for the real population. Cross-validation was performed on 1,200 and 3,117 F animals in the simulated and real populations, respectively. Simulated scenarios were under no artificial selection, mass selection, or BLUP selection. Two genomic relationship matrices were constructed based on breed-specific allele frequencies: 1) , a genomic relationship matrix centered by breed-specific allele frequencies, and 2) , a genomic relationship matrix centered and scaled by breed-specific allele frequencies. All (the across-breed genomic relationship matrix), , and were also tuned to account for selective genotyping. Using breed-specific allele frequencies reduced the number of negative relationships between 2 purebreeds, pulling the average closer to 0, as in the pedigree-based relationship matrix. For simulated populations that included mass selection, genomic EBV (GEBV) in F, when using and , were, on average, 10% more accurate than ; however, after tuning to account for selective genotyping, provided the same accuracy as for breed-specific genomic relationship matrices. For the real population, accuracies for litter size in F were 0.62 for , , and , and tuning had no impact on accuracy, except for , which was 1 percentage point less accurate. Accuracy of GEBV for number of stillborns in F1 was 0.5 for all tested genomic relationship matrices with no changes after tuning. We observed that genotyping F increased accuracies of GEBV for the same animals by up to 39% compared with having genotypes for only AA and BB. In crossbreed evaluations, accounting for breed-specific allele frequencies promoted changes in G that were not influential enough to improve accuracy of GEBV. Therefore, the best performance of ssGBLUP for crossbreed evaluations requires genotypes for pure- and crossbreeds and no breed-specific adjustments in the realized relationship matrix.

Genomic evaluation by including dominance effects and inbreeding depression for purebred and crossbred performance with an application in pigs

BACKGROUND

Improved performance of crossbred animals is partly due to heterosis. One of the major genetic bases of heterosis is dominance, but it is seldom used in pedigree-based genetic evaluation of livestock. Recently, a trivariate genomic best linear unbiased prediction (GBLUP) model including dominance was developed, which can distinguish purebreds from crossbred animals explicitly. The objectives of this study were: (1) methodological, to show that inclusion of marker-based inbreeding accounts for directional dominance and inbreeding depression in purebred and crossbred animals, to revisit variance components of additive and dominance genetic effects using this model, and to develop marker-based estimators of genetic correlations between purebred and crossbred animals and of correlations of allele substitution effects between breeds; (2) to evaluate the impact of accounting for dominance effects and inbreeding depression on predictive ability for total number of piglets born (TNB) in a pig dataset composed of two purebred populations and their crossbreds. We also developed an equivalent model that makes the estimation of variance components tractable.

RESULTS

For TNB in Danish Landrace and Yorkshire populations and their reciprocal crosses, the estimated proportions of dominance genetic variance to additive genetic variance ranged from 5 to 11%. Genetic correlations between breeding values for purebred and crossbred performances for TNB ranged from 0.79 to 0.95 for Landrace and from 0.43 to 0.54 for Yorkshire across models. The estimated correlation of allele substitution effects between Landrace and Yorkshire was low for purebred performances, but high for crossbred performances. Predictive ability for crossbred animals was similar with or without dominance. The inbreeding depression effect increased predictive ability and the estimated inbreeding depression parameter was more negative for Landrace than for Yorkshire animals and was in between for crossbred animals.

CONCLUSIONS

Methodological developments led to closed-form estimators of inbreeding depression, variance components and correlations that can be easily interpreted in a quantitative genetics context. Our results confirm that genetic correlations of breeding values between purebred and crossbred performances within breed are positive and moderate. Inclusion of dominance in the GBLUP model does not improve predictive ability for crossbred animals, whereas inclusion of inbreeding depression does.

Genomic prediction of crossbred performance based on purebred Landrace and Yorkshire data using a dominance model

Background

In pig breeding, selection is usually carried out in purebred populations, although the final goal is to improve crossbred performance. Genomic selection can be used to select purebred parental lines for crossbred performance. Dominance is the likely genetic basis of heterosis and explicitly including dominance in the genomic selection model may be an advantage when selecting purebreds for crossbred performance. Our objectives were two-fold: (1) to compare the predictive ability of genomic prediction models with additive or additive plus dominance effects, when the validation criterion is crossbred performance; and (2) to compare the use of two pure line reference populations to a single combined reference population.

Methods

We used data on litter size in the first parity from two pure pig lines (Landrace and Yorkshire) and their reciprocal crosses. Training was performed (1) separately on pure Landrace (2085) and Yorkshire (2145) sows and (2) the two combined pure lines (4230), which were genotyped for 38 k single nucleotide polymorphisms (SNPs). Prediction accuracy was measured as the correlation between genomic estimated breeding values (GEBV) of pure line boars and mean corrected crossbred-progeny performance, divided by the average accuracy of mean-progeny performance. We evaluated a model with additive effects only (MA) and a model with both additive and dominance effects (MAD). Two types of GEBV were computed: GEBV for purebred performance (GEBV) based on either the MA or MAD models, and GEBV for crossbred performance (GEBV-C) based on the MAD. GEBV-C were calculated based on SNP allele frequencies of genotyped animals in the opposite line.

Results

Compared to MA, MAD improved prediction accuracy for both lines. For MAD, GEBV-C improved prediction accuracy compared to GEBV. For Landrace (Yorkshire) boars, prediction accuracies were equal to 0.11 (0.32) for GEBV based on MA, and 0.13 (0.34) and 0.14 (0.36) for GEBV and GEBV-C based on MAD, respectively. Combining animals from both lines into a single reference population yielded higher accuracies than training on each pure line separately. In conclusion, the use of a dominance model increased the accuracy of genomic predictions of crossbred performance based on purebred data.

Electronic supplementary material

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

Genomic correlation: harnessing the benefit of combining two unrelated populations for genomic selection

BACKGROUND

The success of genomic selection in animal breeding hinges on the availability of a large reference population on which genomic-based predictions of additive genetic or breeding values are built. Here, we explore the benefit of combining two unrelated populations into a single reference population.

METHODS

The datasets consisted of 1829 Brahman and 1973 Tropical Composite cattle with measurements on five phenotypes relevant to tropical adaptation and genotypes for 71,726 genome-wide single nucleotide polymorphisms (SNPs). The underlying genomic correlation for the same phenotype across the two breeds was explored on the basis of consistent linkage disequilibrium (LD) phase and marker effects in both breeds.

RESULTS

The proportion of genetic variance explained by the entire set of SNPs ranged from 37.5 to 57.6 %. Estimated genomic correlations were drastically affected by the process used to select SNPs and went from near 0 to more than 0.80 for most traits when using the set of SNPs with significant effects and the same LD phase in the two breeds. We found that, by carefully selecting the subset of SNPs, the missing heritability can be largely recovered and accuracies in genomic predictions can be improved six-fold. However, the increases in accuracy might come at the expense of large biases.

CONCLUSIONS

Our results offer hope for the effective implementation of genomic selection schemes in situations where the number of breeds is large, the sample size within any single breed is small and the breeding objective includes many phenotypes.

Strategies for genotype imputation in composite beef cattle

Background

Genotype imputation has been used to increase genomic information, allow more animals in genome-wide analyses, and reduce genotyping costs. In Brazilian beef cattle production, many animals are resulting from crossbreeding and such an event may alter linkage disequilibrium patterns. Thus, the challenge is to obtain accurately imputed genotypes in crossbred animals. The objective of this study was to evaluate the best fitting and most accurate imputation strategy on the MA genetic group (the progeny of a Charolais sire mated with crossbred Canchim X Zebu cows) and Canchim cattle. The data set contained 400 animals (born between 1999 and 2005) genotyped with the Illumina BovineHD panel. Imputation accuracy of genotypes from the Illumina-Bovine3K (3K), Illumina-BovineLD (6K), GeneSeek-Genomic-Profiler (GGP) BeefLD (GGP9K), GGP-IndicusLD (GGP20Ki), Illumina-BovineSNP50 (50K), GGP-IndicusHD (GGP75Ki), and GGP-BeefHD (GGP80K) to Illumina-BovineHD (HD) SNP panels were investigated. Seven scenarios for reference and target populations were tested; the animals were grouped according with birth year (S1), genetic groups (S2 and S3), genetic groups and birth year (S4 and S5), gender (S6), and gender and birth year (S7). Analyses were performed using FImpute and BEAGLE software and computation run-time was recorded. Genotype imputation accuracy was measured by concordance rate (CR) and allelic R square (R²).

Results

The highest imputation accuracy scenario consisted of a reference population with males and females and a target population with young females. Among the SNP panels in the tested scenarios, from the 50K, GGP75Ki and GGP80K were the most adequate to impute to HD in Canchim cattle. FImpute reduced computation run-time to impute genotypes from 20 to 100 times when compared to BEAGLE.

Conclusion

The genotyping panels possessing at least 50 thousands markers are suitable for genotype imputation to HD with acceptable accuracy. The FImpute algorithm demonstrated a higher efficiency of imputed markers, especially in lower density panels. These considerations may assist to increase genotypic information, reduce genotyping costs, and aid in genomic selection evaluations in crossbred animals.

Electronic supplementary material

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