Relationship between purebred and crossbred parameters

Selection for social genetic effects in purebreds increases growth in crossbreds

Background

Average daily gain (ADG) in pigs is affected by both direct and social genetic effects (SGE). However, selection for SGE in purebreds has not conclusively been shown to improve ADG in crossbreds, and it is unknown whether SGE in purebreds are equal to those in crossbreds. Moreover, SGE may reflect dominance related behaviour, which is affected by the variation in body weight within a group. Therefore, we hypothesized that (a) there is a positive effect of parent average SGE estimated in purebred pigs on phenotypic ADG in crossbred offspring, and (b) there is an interaction between SGE on ADG and standard deviation in starting weight of pigs within the group. We also hypothesized that (c) social genetic variance for ADG exists in crossbred pigs, and (d) there is a favourable genetic correlation between SGE on ADG in purebred and crossbred pigs.

Results

We found a statistically significant interaction between the standard deviation in starting weight and SGE within groups, and conditioning on the mean standard deviation in starting weight, we found a favourable regression coefficient (0.37 ± 0.21) of ADG in crossbreds on SGE in purebreds. Variances for SGE were small in both Landrace (L) and Yorkshire (Y), and higher for SGE in both the dam and sire component of crossbred YL. The genetic correlations between SGE in purebreds and the dam or sire component of SGE in crossbreds were also favourable (0.52 ± 0.48 and 0.34 ± 0.42, respectively), although not significantly different from 0.

Conclusions

We confirmed that there is a positive effect of SGE estimated using purebred information on phenotypic ADG in crossbreds, and that the largest effect is achieved when the within-group variation in starting weight is small. Our results indicate that social genetic variance in crossbreds exists and that there is a favourable genetic correlation between social genetic effects in purebreds and crossbreds. Collectively, our results indicate that selection for SGE on ADG in purebreds in a nucleus farm environment with little competition for resources can improve ADG in crossbreds in a commercial environment.

Optimizing genomic reference populations to improve crossbred performance

Background

In pig and poultry breeding, the objective is to improve the performance of crossbred production animals, while selection takes place in the purebred parent lines. One way to achieve this is to use genomic prediction with a crossbred reference population. A crossbred reference population benefits from expressing the breeding goal trait but suffers from a lower genetic relatedness with the purebred selection candidates than a purebred reference population. Our aim was to investigate the benefit of using a crossbred reference population for genomic prediction of crossbred performance for: (1) different levels of relatedness between the crossbred reference population and purebred selection candidates, (2) different levels of the purebred-crossbred correlation, and (3) different reference population sizes. We simulated a crossbred breeding program with 0, 1 or 2 multiplication steps to generate the crossbreds, and compared the accuracy of genomic prediction of crossbred performance in one generation using either a purebred or a crossbred reference population. For each scenario, we investigated the empirical accuracy based on simulation and the predicted accuracy based on the estimated effective number of independent chromosome segments between the reference animals and selection candidates.

Results

When the purebred-crossbred correlation was 0.75, the accuracy was highest for a two-way crossbred reference population but similar for purebred and four-way crossbred reference populations, for all reference population sizes. When the purebred-crossbred correlation was 0.5, a purebred reference population always resulted in the lowest accuracy. Among the different crossbred reference populations, the accuracy was slightly lower when more multiplication steps were used to create the crossbreds. In general, the benefit of crossbred reference populations increased when the size of the reference population increased. All predicted accuracies overestimated their corresponding empirical accuracies, but the different scenarios were ranked accurately when the reference population was large.

Conclusions

The benefit of a crossbred reference population becomes larger when the crossbred population is more related to the purebred selection candidates, when the purebred-crossbred correlation is lower, and when the reference population is larger. The purebred-crossbred correlation and reference population size interact with each other with respect to their impact on the accuracy of genomic estimated breeding values.

Genetic parameters of milk production traits in response to a short once-daily milking period in crossbred Holstein × Normande dairy cows

Despite its potential utility for predicting cows' milk yield responses to once-daily milking (ODM), the genetic basis of cow milk trait responses to ODM has been scarcely if ever described in the literature, especially for short ODM periods. This study set out to (1) estimate the genetic determinism of milk yield and composition during a 3-wk ODM period, (2) estimate the genetic determinism of milk yield responses (i.e., milk yield loss upon switching cows to ODM and milk yield recovery upon switching them back to twice-daily milking; TDM), and (3) seek predictors of milk yield responses to ODM, in particular using the first day of ODM. Our trial used 430 crossbred Holstein × Normande cows and comprised 3 successive periods: 1 wk of TDM (control), 3 wk of ODM, and 2 wk of TDM. Implementing ODM for 3 wk reduced milk yield by 27.5% on average, and after resuming TDM cows recovered on average 57% of the milk lost. Heritability estimates in the TDM control period and 3-wk ODM period were, respectively, 0.41 and 0.35 for milk yield, 0.66 and 0.61 for milk fat content, 0.60 and 0.80 for milk protein content, 0.66 and 0.36 for milk lactose content, and 0.20 and 0.15 for milk somatic cell score content. Milk yield and composition during 3-wk ODM and TDM periods were genetically close (within-trait genetic correlations between experimental periods all exceeding 0.80) but were genetically closer within the same milking frequency. Heritabilities of milk yield loss observed upon switching cows to ODM (0.39 and 0.34 for milk yield loss in kg/d and %, respectively) were moderate and similar to milk yield heritabilities. Milk yield recovery (kg/d) upon resuming TDM was a trait of high heritability (0.63). Because they are easy to measure, TDM milk yield and composition and milk yield responses on the first day of ODM were investigated as predictors of milk yield responses to a 3-wk ODM to easily detect animals that are well adapted to ODM. Twice-daily milking milk yield and composition were found to be partly genetically correlated with milk yield responses but not closely enough for practical application. With genetic correlations of 0.98 and 0.96 with 3-wk ODM milk yield losses (in kg/d and %, respectively), milk yield losses on the first day of ODM proved to be more accurate in predicting milk yield responses on longer term ODM than TDM milk yield.

Genomic preselection with genotyping-by-sequencing increases performance of commercial oil palm hybrid crosses

Background

There is great potential for the genetic improvement of oil palm yield. Traditional progeny tests allow accurate selection but limit the number of individuals evaluated. Genomic selection (GS) could overcome this constraint. We estimated the accuracy of GS prediction of seven oil yield components using A × B hybrid progeny tests with almost 500 crosses for training and 200 crosses for independent validation. Genotyping-by-sequencing (GBS) yielded +5000 single nucleotide polymorphisms (SNPs) on the parents of the crosses. The genomic best linear unbiased prediction method gave genomic predictions using the SNPs of the training and validation sets and the phenotypes of the training crosses. The practical impact was illustrated by quantifying the additional bunch production of the crosses selected in the validation experiment if genomic preselection had been applied in the parental populations before progeny tests.

Results

We found that prediction accuracies for cross values plateaued at 500 to 2000 SNPs, with high (0.73) or low (0.28) values depending on traits. Similar results were obtained when parental breeding values were predicted. GS was able to capture genetic differences within parental families, requiring at least 2000 SNPs with less than 5% missing data, imputed using pedigrees. Genomic preselection could have increased the selected hybrids bunch production by more than 10%.

Conclusions

Finally, preselection for yield components using GBS is the first possible application of GS in oil palm. This will increase selection intensity, thus improving the performance of commercial hybrids. Further research is required to increase the benefits from GS, which should revolutionize oil palm breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4179-3) contains supplementary material, which is available to authorized users.

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.

Genetic evaluation for three-way crossbreeding

BACKGROUND

Commercial pig producers generally use a terminal crossbreeding system with three breeds. Many pig breeding organisations have started to use genomic selection for which genetic evaluation is often done by applying single-step methods for which the pedigree-based additive genetic relationship matrix is replaced by a combined relationship matrix based on both marker genotypes and pedigree. Genomic selection is implemented for purebreds, but it also offers opportunities for incorporating information from crossbreds and selecting for crossbred performance. However, models for genetic evaluation for the three-way crossbreeding system have not been developed.

RESULTS

Four-variate models for three-way terminal crossbreeding are presented in which the first three variables contain the records for the three pure breeds and the fourth variable contains the records for the three-way crossbreds. For purebred animals, the models provide breeding values for both purebred and crossbred performances. Heterogeneity of genetic architecture between breeds and genotype by environment interactions are modelled through genetic correlations between these breeding values. Specification of the additive genetic relationships is essential for these models and can be defined either within populations or across populations. Based on these two types of additive genetic relationships, both pedigree-based, marker-based and combined relationships based on both pedigree and marker information are presented. All these models for three-way crossbreeding can be formulated using Kronecker matrix products and therefore fitted using Henderson's mixed model equations and standard animal breeding software.

CONCLUSIONS

Models for genetic evaluation in the three-way crossbreeding system are presented. They provide estimated breeding values for both purebred and crossbred performances, and can use pedigree-based or marker-based relationships, or combined relationships based on both pedigree and marker information. This provides a framework that allows information from three-way crossbred animals to be incorporated into a genetic evaluation system.

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

Background

The accuracy of genomic prediction depends largely on the number of animals with phenotypes and genotypes. In some industries, such as sheep and beef cattle, data are often available from a mixture of breeds, multiple strains within a breed or from crossbred animals. The objective of this study was to compare the accuracy of genomic prediction for several economically important traits in sheep when using data from purebreds, crossbreds or a combination of those in a reference population.

Methods

The reference populations were purebred Merinos, crossbreds of Border Leicester (BL), Poll Dorset (PD) or White Suffolk (WS) with Merinos and combinations of purebred and crossbred animals. Genomic breeding values (GBV) were calculated based on genomic best linear unbiased prediction (GBLUP), using a genomic relationship matrix calculated based on 48 599 Ovine SNP (single nucleotide polymorphisms) genotypes. The accuracy of GBV was assessed in a group of purebred industry sires based on the correlation coefficient between GBV and accurate estimated breeding values based on progeny records.

Results

The accuracy of GBV for Merino sires increased with a larger purebred Merino reference population, but decreased when a large purebred Merino reference population was augmented with records from crossbred animals. The GBV accuracy for BL, PD and WS breeds based on crossbred data was the same or tended to decrease when more purebred Merinos were added to the crossbred reference population. The prediction accuracy for a particular breed was close to zero when the reference population did not contain any haplotypes of the target breed, except for some low accuracies that were obtained when predicting PD from WS and vice versa.

Conclusions

This study demonstrates that crossbred animals can be used for genomic prediction of purebred animals using 50 k SNP marker density and GBLUP, but crossbred data provided lower accuracy than purebred data. Including data from distant breeds in a reference population had a neutral to slightly negative effect on the accuracy of genomic prediction. Accounting for differences in marker allele frequencies between breeds had only a small effect on the accuracy of genomic prediction from crossbred or combined crossbred and purebred reference populations.

Genomic evaluation of both purebred and crossbred performances

Background

For a two-breed crossbreeding system, Wei and van der Werf presented a model for genetic evaluation using information from both purebred and crossbred animals. The model provides breeding values for both purebred and crossbred performances. Genomic evaluation incorporates marker genotypes into a genetic evaluation system. Among popular methods are the so-called single-step methods, in which marker genotypes are incorporated into a traditional animal model by using a combined relationship matrix that extends the marker-based relationship matrix to non-genotyped animals. However, a single-step method for genomic evaluation of both purebred and crossbred performances has not been developed yet.

Results

An extension of the Wei and van der Werf model that incorporates genomic information is presented. The extension consists of four steps: (1) the Wei van der Werf model is reformulated using two partial relationship matrices for the two breeds; (2) marker-based partial relationship matrices are constructed; (3) marker-based partial relationship matrices are adjusted to be compatible to pedigree-based partial relationship matrices and (4) combined partial relationship matrices are constructed using information from both pedigree and marker genotypes. The extension of the Wei van der Werf model can be implemented using software that allows inverse covariance matrices in sparse format as input.

Conclusions

A method for genomic evaluation of both purebred and crossbred performances was developed for a two-breed crossbreeding system. The method allows information from crossbred animals to be incorporated in a coherent manner for such crossbreeding systems.

Genetic parameters of the piglet mortality traits stillborn, weak at birth, starvation, crushing, and miscellaneous in crossbred pigs

This study aimed to estimate genetic parameters for the mortality causes stillborn, weak at birth, starvation, crushing, and miscellaneous in crossbred piglets produced by crossbred dams. Data were collected in a single Danish commercial herd from October 2006 to July 2008 and consisted of 34,194 piglets (2,152 litters), which originated from 195 Danish Duroc sires and 955 crossbreds between Danish Landrace and Danish Yorkshire dams. Of the 34,194 piglets born, 11.5% were stillborn, 4.2% were crushed by the sow, 2.7% died due to starvation, 2.3% were weak at birth, and 2.2% died of miscellaneous causes before weaning. The first 4 mentioned causes were analyzed multivariately using a generalized linear mixed model with a probit link function, including the genetic effect of both sire and dam. Heritabilities based on the sire component ranged between 0.08 for stillborn and 0.21 for starvation whereas heritabilities based on the dam component ranged between 0.01 for miscellaneous and 0.24 for stillborn, indicating that reducing piglet mortality through genetic selection is possible. The expected observed responses to selection would, however, be low. The genetic correlations between mortality traits based on the sire component ranged from -0.05 between stillborn and starvation to 0.35 between stillborn and weak at birth whereas genetic correlations based on the dam component ranged from -0.11 between weak at birth and starvation to 0.76 between crushing and starvation. There seemed to be a favorable relationship between the 2 traits stillborn and weak at birth and between crushing and starvation, implying that care should be taken with correct recordings of mortality causes. The genetic correlation precision was rather low, and if nonadditive effects are not accounted for, there may be unexpected correlated responses among the different mortality causes in the crossbred mortalities.

The value of crossbreeding information in selection of poultry under a dominance model

A stochastic simulation was used to investigate the value of crossbreeding information in a two line crossbreeding system in poultry under different genetic scenarios. Populations consisted of 25 sires and 250 dams, which were mated to produce both purebred and crossbred progeny. The next generation parents were selected either based on purebred information (PLS) (sib means, and own performance for females) or additional crossbred sib means were included (CCPS). The trait under selection was controlled by 20 loci with varying degree of dominance. Pure lines differed in initial allele frequencies. If the trait was controlled by loci with partial dominance, little or no extra benefit was obtained from including crossbred information over the pure line information. Under complete dominance and overdominance CCPS outperformed PLS. As a practical rule, CCPS is better than PLS if the ratio between dominance variance and total genetic variance is around 0·3 or higher. In this case the most probable cause of the dominance variation is loci with full or overdominance.

Maximizing genetic gain for the sire line of a crossbreeding scheme utilizing both purebred and crossbred information

A selection index procedure which utilizes both, purebred and crossbred information was developed for the sire line of a three-path crossbreeding scheme in pigs, to predict response to best linear unbiased prediction (BLUP) selection with an animal model. Purebred and crossbred performance were treated as correlated traits. The breeding goal was crossbred performance but methods can be applied to other goals. A hierarchical mating structure was used. Sires were mated to purebred dams to generate replacements and to F^ from the dam line to generate fattening pigs. Generations were discrete, inbreeding was ignored. The selection index included purebred and crossbred phenotypic information of the current generation and estimated breeding values for purebred and crossbred performance of parents and mates of parents from the previous generation. Reduction of genetic variance due to linkage disequilibrium and reduction of selection intensity due to finite population size and due to correlated index values was accounted for. Selection was undertaken until asymptotic responses were reached. The index was used to optimize the number of selected parents per generation and the number of offspring tested per litter, and to make inferences on the value of crossbred information when the breeding goal was crossbred performance. It was optimal to test a maximum number of offspring per litter, mainly due to increased female selection intensities. Maximum response reductions due to linkage disequilibrium and correlated index values were 32% and 29% respectively. Correcting for correlated index values changed ranking of breeding schemes. Benefit of crossbred information was largest when the genetic correlation between purebred and crossbred performance was low. Due to high correlations between index values in that case, the optimum number of selected sires increased considerably when crossbred information was included.

Maximizing genetic response in crossbreds using both purebred and crossbred information

A combined crossbred and purebred selection (CCPS) method, i.e. using crossbred and purebred information, was proposed to achieve genetic response in crossbred animals. Selection index theory was applied to establish a CCPS index. The CCPS was compared with pure-line selection (PLS) and crossbred selection (CS) methods. The genetic correlation between purebred and crossbred performance (rpc) and crossbred heritability (hc2) are crucial factors in the comparison. The CCPS is always better than PLS or CS when a fixed number of purebred progeny is tested. With a fixed total number of purebred and crossbred tested progeny, CCPS is only worse than PLS for very high values of rpc (>0·8). Superiority of CCPS over PLS increases and over CS decreases with decreasing rpc. The larger hc2 is, relative to purebred heritability (hc2 the more response CS and CCPS will achieve. The robustness of CCPS against inappropriate assumptions on rpc and hc2 values was investigated. The expected response is always an overestimate, and the actual response is smaller than the optimal response when rpc is assumed one but the true rpc is smaller. The difference between actual and optimal response increases as rpc decreases but it is small for large rpc values (e.g. <3% for rpc >0·7). The expected response is smaller than the actual response when rpc is large and hc2> hp2 Finally, the actual response to CCPS is larger than the optimal response to PLS for positive values for rpc. The main conclusions are: (1) CCPS method is optimal for obtaining genetic response in crossbreds; and (2) CCPS with inappropriate assumptions on rpc and hc2 values (e.g. recognizing crossbreds as purebreds) achieves more genetic response than PLS for common values of rpc and crossbred heritability.

Genetic correlations between two strains of Durocs and crossbreds from differing production environments for slaughter traits

The aim of this study was to estimate the genetic correlations between 2 purebred Duroc pig populations (P1 and P2) and their terminal crossbreds [C1 = P1 x (Landrace x Large White) and C2 = P2 x (Landrace x Large White)] raised in different production environments. The traits analyzed were backfat (BF), muscle depth (MD), BW at slaughter (WGT), and weight per day of age (WDA). Data sets from P1, P2, C1, and C2 included 26,674, 8,266, 16,806, and 12,350 animals, respectively. Two-trait models (nucleus and commercial crossbreds) for each group included fixed (contemporary group, sex, weight, and age), random additive (animal for P1 and P2 and sire for C1 and C2), random litter, and random dam (C1 and C2 only) effects. Heritability estimates (+/-SE) for BF were 0.46 +/- 0.04, 0.38 +/- 0.02, 0.32 +/- 0.02, and 0.33 +/- 0.02 for P1, P2, C1, and C2, respectively. Heritability estimates for MD were 0.31 +/- 0.01, 0.23 +/- 0.02, 0.19 +/- 0.01, and 0.12 +/- 0.01 for P1, P2, C1, and C2, respectively. The estimates for WGT and WDA were 0.31 +/- 0.01, 0.21 +/- 0.02, 0.16 +/- 0.01, and 0.18 +/- 0.01 and 0.32 +/- 0.01, 0.22 +/- 0.02, 0.16 +/- 0.01, and 0.19 +/- 0.01, respectively. Genetic correlations between purebreds and crossbreds for BF were 0.83 +/- 0.09 (P1 x C1) and 0.89 +/- 0.05 (P2 x C2), for MD 0.78 +/- 0.05 (P1 x C1) and 0.80 +/- 0.08 (P2 x C2). For WGT and WDA, the correlations were 0.53 +/- 0.08 (P1 x C1), 0.80 +/- 0.10 (P2 x C2), and 0.60 +/- 0.07 (P1 x C1) and 0.79 +/- 0.09 (P2 x C2), respectively. (Co)variances in crossbreds were adjusted to a live BW scale. Compared with purebreds, the genetic variances in crossbreds were lower, and the residual variances were greater. Sire variances in crossbreds were approximately 20 to 30% of the animal variances in purebreds for BF and MD but were 13 to 25% for WGT and WDA. The efficiency of purebred selection on crossbreds, assessed by EBV prediction weights, ranged from 0.43 to 0.91 for line 1 and 0.70 to 0.92 for line 2. When nucleus and commercial environments differ substantially, the efficiency of selection varies by line and traits, and selection strategies that include crossbred data from typical production environments may therefore be desirable.

Theory for modelling means and covariances in a two-breed population with dominance inheritance

This paper presents theory and methods to compute genotypic means and covariances in a two breed population under dominance inheritance, assuming multiple unlinked loci. It is shown that the genotypic mean is a linear function of five location parameters and that the genotypic covariance between relatives is a linear function of 25 dispersion parameters. Recursive procedures are given to compute the necessary identity coefficients. In the absence of inbreeding, the number of parameters for the mean is reduced from five to three and the number for the covariance is reduced from 25 to 12. In a two-breed population, for traits exhibiting dominance, the theory presented here can be used to obtain genetic evaluations by best linear unbiased prediction and to estimate genetic parameters by maximum likelihood.

Covariance between relatives in multibreed populations: additive model

Covariance between relatives in a multibreed population was derived for an additive model with multiple unlinked loci. An efficient algorithm to compute the inverse of the additive genetic covariance matrix is given. For an additive model, the variance for a crossbred individual is a function of the additive variances for the pure breeds, the covariance between parents, and segregation variances. Provided that the variance of a crossbred individual is computed as presented here, the covariance between crossbred relatives can be computed using formulae for purebred populations. For additive traits the inverse of the genotypic covariance matrix given here can be used both to obtain genetic evaluations by best linear unbiased prediction and to estimate genetic parameters by maximum likelihood in multibreed populations. For nonadditive traits, the procedure currently used to analyze multibreed data can be improved using the theory presented here to compute additive covariances together with a suitable approximation for nonadditive covariances.