Accuracy of pedigree and genomic predictions of carcass and novel meat quality traits in multi-breed sheep data assessed by cross-validation
© Daetwyler et al.; licensee BioMed Central Ltd. 2012
Received: 5 July 2012
Accepted: 31 October 2012
Published: 12 November 2012
Genomic predictions can be applied early in life without impacting selection candidates. This is especially useful for meat quality traits in sheep. Carcass and novel meat quality traits were predicted in a multi-breed sheep population that included Merino, Border Leicester, Polled Dorset and White Suffolk sheep and their crosses.
Prediction of breeding values by best linear unbiased prediction (BLUP) based on pedigree information was compared to prediction based on genomic BLUP (GBLUP) and a Bayesian prediction method (BayesR). Cross-validation of predictions across sire families was used to evaluate the accuracy of predictions based on the correlation of predicted and observed values and the regression of observed on predicted values was used to evaluate bias of methods. Accuracies and regression coefficients were calculated using either phenotypes or adjusted phenotypes as observed variables.
Results and conclusions
Genomic methods increased the accuracy of predicted breeding values to on average 0.2 across traits (range 0.07 to 0.31), compared to an average accuracy of 0.09 for pedigree-based BLUP. However, for some traits with smaller reference population size, there was no increase in accuracy or it was small. No clear differences in accuracy were observed between GBLUP and BayesR. The regression of phenotypes on breeding values was close to 1 for all methods, indicating little bias, except for GBLUP and adjusted phenotypes (regression = 0.78). Accuracies calculated with adjusted (for fixed effects) phenotypes were less variable than accuracies based on unadjusted phenotypes, indicating that fixed effects influence the latter. Increasing the reference population size increased accuracy, indicating that adding more records will be beneficial. For the Merino, Polled Dorset and White Suffolk breeds, accuracies were greater than for the Border Leicester breed due to the smaller sample size and limited across-breed prediction. BayesR detected only a few large marker effects but one region on chromosome 6 was associated with large effects for several traits. Cross-validation produced very similar variability of accuracy and regression coefficients for BLUP, GBLUP and BayesR, showing that this variability is not a property of genomic methods alone. Our results show that genomic selection for novel difficult-to-measure traits is a feasible strategy to achieve increased genetic gain.
Sheep meat production is increasing and replacing wool production as the primary product of the Australian sheep industry. Improving growth traits through selection for increased live and a carcass weight is an important driver of profitability. Providing consistently high-quality meat is also essential to maintain high consumer acceptance and depends on several carcass quality criteria, such as intra-muscular fat, shear force and Omega-3 content[1, 2]. Genomic selection is applied in an ever growing number of livestock species, e.g.[4–8], and could increase economic returns from lamb production. Genomic selection can be applied at a number of well-known entry points of breeding schemes to increase genetic progress. Genomic estimated breeding values (GEBV) can be obtained for selection candidates at a young age before phenotypic information is available and be used to increase accuracy of selection and shorten generation intervals. This is useful for traits measured later in life, such as adult greasy fleece weight and reproduction and in cases when phenotypic evaluation involves invasive or destructive approaches such as for carcass composition and meat quality, which are traditionally measured on the relatives of selection candidates.
Estimating GEBV for selection candidates requires a reference population with both marker genotypes and phenotypes. Because selection candidates often lack phenotypic records, the predictive performance of GEBV can be assessed either with a set of validation individuals that have highly accurate EBV, e.g. sires with many progeny[4, 10], or by cross-validation, e.g.[6, 11–13]. Both validation methods require that the validation population and the potential selection candidates have a similar genetic make-up, such that the accuracies obtained for the selection candidates will reflect those calculated using validation individuals. In particular, the validation and selection individuals should have similar relationships to the reference population[14–16].
For difficult-to-measure and novel traits, individuals with highly accurate EBV often do not exist. Thus, in such cases, cross-validation is applied. In the cross-validation approach, the reference population is divided into a number of subsets and each subset is predicted using a reference population that excludes this particular subset. The method used to divide the data has been shown to affect prediction accuracy, e.g.[6, 17, 18]. Choosing fully random subsets is the simplest implementation but this ignores data structures, for example presence of sire half-sib groups. Several studies have divided subsets randomly with constraints on data structure, such as age, family, and relatedness. Another consideration is the size of the subsets used for cross-validation. The larger the subset, the smaller the sampling variance of the correlation between predicted and observed variables is expected to be[11, 19]. However, larger subsets decrease the size of the reference population, resulting in a trade-off between the size of subsets and the accuracy achieved.
The utility of applying genomic prediction must be evaluated against what would be achieved with non-genomic approaches, such as best linear unbiased prediction (BLUP) using pedigree. Cross-validation studies have compared accuracies of EBV using traditional pedigree methods and accuracies of GEBV[6, 18] but these comparisons have not been made for multi-breed livestock data. Another point to consider is the phenotype used to estimate accuracies in cross-validation studies. Most studies correlate with phenotypes but the accuracies resulting from these comparisons may be affected by fixed effects that are often included in the prediction models.
The aim of this study was to predict GEBV for several carcass and novel meat quality traits in a multi-breed sheep population. A previous study in the same population investigated how much of the accuracy of GEBV could be attributed to population structure. Here, an across sire family cross-validation scheme was used to estimate accuracies of GEBV in several sheep breeds and their crosses. GEBV were obtained with three methods: BLUP, genomic BLUP (GBLUP) and BayesR. In addition, accuracies were calculated based on phenotypes or adjusted phenotypes and with or without adding a polygenic effect.
Summary statistics and heritabilities ( h 2 ) for the phenotypic traits analysed for the two data sets (CRC and SG)
The following traits were analysed and phenotypic information for these traits is provided in Table1. Carcass eye muscle depth (EMD, mm), carcass fat depth at site C (FAT, mm, depth of fat at maximum EMD), hot carcass weight (HCWT, kg), dressing percentage (DRESS, %), calculated as the ratio of HCWT to pre-slaughter weight, intra-muscular fat (IMF, %), iron content of wet muscle tissue (IRON, mg/kg ), and the concentration of omega 3 fatty acid compounds eicosapentaenoic acid (EPA, mg/100g) and docosapentaenoic acid (DPA, mg/100g). Lean meat yield (LMY, %) was estimated on the CRC animals by a combination of other carcass traits and validated by computed tomography (CT) scanning. On the SG animals, LMY was computed as the ratio of HCWT and actual lean meat after bone-out. To account for these differences in methodology, LMY was standardised (mean = 0, standard deviation (SD) = 1) within the CRC and SG datasets before the datasets were merged.
All animals were genotyped using the Illumina 50K ovine SNP chip, containing 54 977 single nucleotide polymorphisms (SNP) (Illumina Inc., San Diego, USA). After applying the following quality control measures, 48 599 SNP were retained: SNP were removed if the call rate was less than 95%, if the Illumina Gentrain score was less than 0.6, if the minor allele frequency was less than 0.01, if the SNP was not in Hardy-Weinberg equilibrium (a P-value cut-off of 1×10–15), if the genome location was unknown or if the SNP showed complete linkage disequilibrium (r2 > 0.99) with another SNP on the chip. Data for a genotyped animal were removed if the genotype call rate was less than 0.9 for that animal or if the animal’s mean heterozygosity was higher than 0.5, indicating sample contamination. The genotype database was built over a number of years, missing genotypes were initially imputed using fastPHASE and more recently, missing genotypes were imputed using Beagle, after this program became available.
where y is a vector of phenotypic records, X, and Z 1 are design matrices relating the fixed and random effects to the phenotype, Q is a matrix containing breed proportions for each animal, derived from pedigree information, μ is the mean, b is a vector of fixed effects, a is a vector of random additive polygenic effects, q is a vector of random breed effects, fitted as partial regressions, and e is the vector of residuals. The following distributions were assumed: a ~ N (0,σ a 2A), q ~ N (0,σ q 2I), and e ~ N (0,σ e 2I), where A is the numerator relationship matrix, σ a 2 is the additive variance, σ q 2 is the variance of breed effects, and σ e 2 is the residual variance. The base model included the following fixed effects: sex, birth type, rearing type, contemporary group (birth year × site × slaughter group), and age at trait recording. Age of the dam was fitted only for CRC data. HCWT was included as a fixed covariate for all traits except for DRESS and LMY. The size of the relevant pedigree was 16 985 individuals. The phenotypes (y) were restricted to genotyped animals to make a fair comparison to the genomic prediction methods.
where Z 2 is a design matrix, g is a vector of random additive genomic effects distributed as N (0,σ g 2G), σ g 2 is the genomic variance, G is the genomic relationship matrix. SNP with allele frequencies less than 0.005 were removed from the calculation of G to improve numerical stability. Phenotypes, rather than de-regressed estimated breeding values, were used to ensure independence of reference and validation sets. If all phenotypes were used to calculate breeding values, then the accuracy of predicting an animal without a phenotype would be overestimated, because phenotypes of validation animals contributed to the reference pedigree breeding values.
where W is a design matrix relating adjusted phenotypes to random marker effects (m). BayesR is described in more detail in. Briefly, marker variances can come from distributions with variances σ12= 0, σ22= 0.0001σ g 2, σ32= 0.001σ g 2, or σ42= 0.01σ g 2, and starting values for σ g 2 were from GBLUP analysis. The prior for the proportion of markers in each distribution was drawn from a Dirichlet distribution. Priors for other parameters were chosen as in Erbe et al.. Ten parallel chains of 50 000 iterations (20 000 burn-in) were run for each subset.
Posterior means of marker effects of BayesR resulting from post burn-in chains were averaged across chains and replicates and then standardised by dividing them by the standard deviation of the adjusted phenotypes (SD). SNP with effects greater than 0.005 SD were (arbitrarily) chosen and potential candidate genes were searched for onhttp://www.livestockgenomics.csiro.au/cgi-bin/gbrowse/oarv2.0/ using a 1 Mb interval with 0.5 Mb on each side of the SNP. The probability of the effect of a SNP being in the largest distribution (σ42) was also investigated.
Performance of genomic prediction was evaluated using cross-validation. It is unlikely that potential selection candidates in the Australian sheep population have full or half-sibs in the reference population. Thus, entire sire families in the CRC dataset were randomly chosen and combined into subsets of approximately 500 individuals (CRC subsets). Thus, genomic predictions were evaluated across sire families and larger reference populations had more subsets. The SG dataset was not divided into subsets but was added to each reference population. The performance of predictions was not evaluated for the SG data, because this population is not expected to be representative of the general sheep population. Genomic predictions were calculated for each CRC subset, with the reference set consisting of all other CRC subsets and the SG dataset. Accuracy was evaluated in each validation subset as the Pearson correlation of genomic predicted breeding values () or genomic plus polygenic predicted breeding values ( +) with either phenotypes (y) or adjusted phenotypes (y*). Accuracies were divided by from model 1 to adjust for the upper limit of accuracy of a phenotype/residual. The bias of breeding values (both and +) was calculated as the regression of phenotypes or adjusted phenotypes on predicted breeding values. Accuracy and bias were calculated for the whole validation subset and for each subdivision of each subset by the sire breeds Merino (MER, effective population size, N e ~ 850), Border Leicester (BL, N e ~ 250), Polled Dorset (PD, N e ~ 300), and White Suffolk (WS, N e ~ 300).
Genetic relatedness of validation animals with the reference population was calculated for each subset as the mean of the squared relationships between validation and reference animals and mean of the top 10 relationships. Other studies have concluded that these measures are more predictive of accuracy than mean relationships, hence the latter was not reported[15, 16].
Results and discussion
Table S3 [see Additional file1: Table S3] contains the mean genetic relationships of validation with reference animals, calculated as the mean of the top 10 genomic relationships for each individual. Small differences in relatedness between breed groups and between traits were observed, ranging from 0.102 for DRESS to 0.168 for IRON, both in the Merino breed. No clear relationships of the mean genetic relationship with the achieved accuracy were found, for several reasons. First, the sampling variances of the correlations between GEBV and phenotypes were too large due to the small size of the validation sets [see Additional file1: Table S4]. Secondly, the genomic relationship matrix used here was based on the original Yang et al. implementation, which does not adjust for breed (base) allele frequencies. While these scaling issues are not likely to decrease predictive performance substantially, they will numerically affect mean relationships within a breed and may limit the possibility of finding a relationship between magnitude of mean relationships and within-breed accuracies. One possible solution would be to scale allele frequencies within breeds to their respective breed base allele frequencies before calculating the relationship matrix.
In the cross-validation design applied here, sires were chosen randomly and all their progeny were assigned to subsets. This prevented the upward bias of accuracies that would result from within-family prediction when half sib families are randomly split between reference and validation datasets. The accuracies obtained with our approach are expected to better reflect what would be achieved across a range of industry selection candidates with varying degrees of relationships to the reference animals. A further complication in our study was that the reference and validation populations were mostly made up of crossbreds, yet potential selection candidates in the industry are likely purebred individuals. Dividing the validation sets by sire breed groups was used to approximate the accuracy of purebred selection candidates. Because all animals have a large Merino component, the accuracies obtained with the Border Leicester, Polled Dorset and White Suffolk validation sets (sire breed groups) are not strictly equivalent to the accuracy which would be obtained with purebred animals. However, in the absence of purebred individuals with carcass data, this represents the best possible approximation. However, while the selection of breeding stock takes place among purebred individuals, commercial stock results from crosses between terminal and maternal sheep breeds.
Another aim of this study was to compare the variability of cross-validated EBV accuracies across subsets from pedigree-based BLUP and genomic methods. No large differences were observed, aside from the increase in accuracy using genomic data. In addition, accuracies from pedigree-based BLUP and genomic methods had very similar standard errors, indicating that cross-validation accuracies are just as variable across subsets for pedigree-based BLUP.
Information of SNP with greater than 0.5 SD effects, including all genes present with 0.5 Mb on either side
TSC21, EIF2AK3, RPIA, IGK, PSD4, IL1RN, IL1F10, IL1F5
SLC20A1, CHCHD5, POLR1B, TTL, NCK2, Augurin, cDNA FLJ78230
PPM1K, ABCG2, PKD2, SPP1, MEPE, IBSP, LAP3, FAM182A, DCAF16, NCAPG, LCORL
CRYABA1, NYFIP2, TAOK1, GIT1, ANDKRD13B, SSH2, EFCAB5, NSRP1, SLC6A4, BLMH, CPD
Chymotrypsinogen A, BCAR1, CFDP1, CFDP2, TMEM170A, BVDV, ADAT1, KARS, CNTNAP4
HEATR3, PAPD5, ADCY7, BRD7, NKD1, NX20, NOD2, CYLD
Potential candidate genes within a 1 Mb interval with 0.5 Mb on each side of the SNP are also presented in Table2. One region on chromosome 6 contained SNP with estimated effects ranging from 0.0060 to 0.0106 SD for FAT, IMF, DRESS and LMY. Genes in this region included ATP-binding cassette sub-family G member 2 (ABCG2) and Polycystin-2 (PKD2), which have been reported as having been under selection in an analysis of a large number of sheep breeds. ABCG2, a gene involved in ATP binding, has been found to contain a causative mutation that affects milk yield and composition in dairy cattle and has also been investigated as a candidate gene for facial eczema in sheep. Another potential gene of interest is NCK2 protein. This gene is close to SNP OAR3_64213489 (DRESS 0.0052SD) and codes for an adaptor protein that associates with tyrosine-phosphorylated growth-factor receptors. ARF GTPase-activating protein (GIT1), which is close to SNP OAR11_21345650_X (IMF, 0.0058SD), codes for a GTPase-activating protein that is possibly involved in vesicle trafficking, adhesions and cytoskeletal organisation.
The main benefit of genomic selection for carcass and novel meat quality traits is that it is not necessary to sacrifice valuable selection candidates for testing, and GEBV can be obtained early in life. Genomic predictions can be trained within a set of industry representative individuals and then applied in the general sheep population. In its current form, this process is implemented in a centralised approach in sheep, in which test animals are housed in information nucleus flocks. However, data for training could also be collected during slaughter of industry stock, which could substantially increase the number of records. One advantage of an information nucleus is that animals are well identified and similarly managed, and fixed effects are fully recorded. Including industry records would require further investment in recording and tracking of animal production and movement and the development of uniform standards for measurement, sampling and testing at slaughter facilities.
One application of genomic selection is the prediction of accurate breeding values of juveniles without phenotypic records. This allows for significant shortening of generation intervals. In addition, some sheep breeders use juvenile in-vitro fertilised embryo transfer (JIVET), which consists of harvesting immature oocytes from 6 to 8 week old ewe lambs and implanting these into sexually mature individuals after in vitro fertilisation. The combination of genomic selection and JIVET could be a powerful tool to increase genetic gain for novel meat traits. For example, lines with high omega-3 content or superior eating quality could be developed. The increase in genetic gain resulting from genomic selection would need to be combined with a strategy to limit a reduction in genetic diversity, such as using optimised contributions and mating schemes[37–40].
The accuracies of EBV obtained with genomic methods were not substantially higher than accuracies of EBV obtained with pedigree-based BLUP. Given the large reference population size, one could have expected a larger increase. However, in our case, many breeds contributed to the reference population and the limited increase in accuracy could be explained by small contributions obtained from across-breed prediction, which has been found to be very low in this population. Increasing marker density, either through a high-density SNP chip or through next-generation sequencing, could increase accuracies from across-breed prediction. Simulation studies with simple genetic architectures have shown that using sequence data can be beneficial. In contrast, in a perhaps under-powered study using empirical Drosophila sequence data, no increase in genomic prediction accuracy was observed compared to the use of lower SNP densities. In dairy cattle, there is some evidence that across-breed prediction can be increased when using a high-density array. A higher density SNP chip could potentially be more beneficial in sheep than in cattle, because the N e is greater for most sheep breeds than in dairy cattle. In Holstein cattle, 80% of the genetic variance was captured by the 50k bovine SNP chip. Using the same methodology between 30 to 55% of the genetic variance was captured by the 50k ovine SNP chip in Merino sheep, depending on trait (results not shown). Thus, increasing marker density may result in substantial increases in accuracy, both within and across sheep breeds. Currently, the implementation of high-density arrays, and potentially sequence data, is accomplished using a two-step approach. First, reference populations that are genotyped at medium density (e.g. 50 000 SNP) are imputed up to higher density, using a sample of individuals that is genotyped at the higher density. Second, the imputed reference population is used for genomic prediction.
Genomic predictions for meat quality traits in sheep are potentially valuable because they can be applied early in life and do not require potential selection candidates to be sacrificed. In a large multi-breed sheep dataset, genomic prediction resulted in greater accuracies of EBV than pedigree-based BLUP, but for some traits the increase in accuracy was small. Accuracy increased as reference population size increased and the accuracy was greater for the Merino, Polled Dorset and White Suffolk breeds than for the Border Leicester breed. The latter result is explained in part by the lower proportion of Border Leicester sheep in the reference population. It also suggests that across-breed prediction is limited with the 50k SNP chip. The methods GBLUP and BayesR produced very similar accuracies of GEBV, with a mean accuracy of approximately 0.2 across traits. Few markers with large effects were discovered but one region on chromosome 6 was associated with large effects for several traits. Validation correlations of GEBV with phenotypes adjusted for estimates of fixed effects were less variable than correlations with unadjusted phenotypes, and there was little evidence of bias in the GEBV. The general behaviour of cross-validation accuracies was very similar for pedigree-based BLUP, GBLUP and BayesR. In conclusion, genomic breeding values can provide a powerful tool to increase genetic progress in sheep, especially when combined with reproductive technologies.
The authors gratefully acknowledge funding from the Cooperative Research Centre for Sheep Industry Innovation, Meat and Livestock Australia, and Australian Wool Innovation Ltd. We thank Klint Gore and Ken Geenty for managing the CRC information nucleus database, Cedric Gondro for performing part of the genotype quality control for the CRC data and the many staff involved at the CRC and SG sites across Australia. We thank the reviewers for their constructive comments.
- Rowe JB: The Australian sheep industry - undergoing transformation. Anim Prod Sci. 2010, 50: 991-997. 10.1071/AN10142.View Article
- Pethick D, Banks RG, Hales J, Ross JR: Australian prime lamb - a vision for 2020. Int J Sheep Wool Sci. 2006, 54: 66-73.
- Meuwissen THE, Hayes BJ, Goddard ME: Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001, 157: 1819-1829.PubMed CentralPubMed
- Daetwyler HD, Hickey JM, Henshall JM, Dominik S, Gredler B, van der Werf JHJ, Hayes BJ: Accuracy of estimated genomic breeding values for wool and meat traits in a multi-breed sheep population. Anim Prod Sci. 2010, 50: 1004-1010. 10.1071/AN10096.View Article
- Wolc A, Stricker C, Arango J, Settar P, Fulton JE, O'Sullivan NP, Preisinger R, Habier D, Fernando R, Garrick DJ, Lamont SJ, Dekkers JC: Breeding value prediction for production traits in layer chickens using pedigree or genomic relationships in a reduced animal model. Genet Sel Evol. 2011, 43: 5-10.1186/1297-9686-43-5.PubMed CentralView ArticlePubMed
- Saatchi M, McClure MC, McKay SD, Rolf MM, Kim J, Decker JE, Taxis TM, Chapple RH, Ramey HR, Northcutt SL, Bauck S, Woodward B, Dekkers JC, Fernando RL, Schnabel RD, Garrick DJ, Taylor JF: Accuracies of genomic breeding values in American Angus beef cattle using K-means clustering for cross-validation. Genet Sel Evol. 2011, 43: 40-10.1186/1297-9686-43-40.PubMed CentralView ArticlePubMed
- Hayes BJ, Bowman PJ, Chamberlain AJ, Goddard ME: Invited review: Genomic selection in dairy cattle: progress and challenges. J Dairy Sci. 2009, 92: 433-443. 10.3168/jds.2008-1646.View ArticlePubMed
- Lund MS, de Ross APW, de Vries AG, Druet T, Ducrocq V, Fritz S, Guillaume F, Guldbrandtsen B, Liu Z, Reents R, Schrooten C, Seefried F, Su G: A common reference population from four European Holstein populations increases reliability of genomic predictions. Genet Sel Evol. 2011, 43: 43-10.1186/1297-9686-43-43.PubMed CentralView ArticlePubMed
- Banks RG, van der Werf JHJ: Economic evaluation of whole genome selection, using meat sheep as a case study. Proceedings of the 18th Conference of the Association for the Advancement of Animal Breeding and Genetics: 28 September - 1 October 2009; Barossa Valley. 2009, AAABG Distributors, Armidale, Australia, 430-433.
- VanRaden PM, Van Tassell CP, Wiggans GR, Sonstegard TS, Schnabel RD, Taylor JF, Schenkel FS: Invited review: Reliability of genomic predictions for North American Holstein bulls. J Dairy Sci. 2009, 92: 16-24. 10.3168/jds.2008-1514.View ArticlePubMed
- Erbe M, Pimentel ECG, Sharifi AR, Simianer H: Assessment of cross-validation strategies for genomic prediction in cattle. 9th World Congress of Genetics Applied to Livestock Production: 1–6 August 2009; Leipzig. 2010, Gesellschaft für Tierzuchtwissenschaften e. V, Giessen, Germany
- Legarra A, Robert-Granié C, Manfredi E, Elsen JM: Performance of genomic selection in mice. Genetics. 2008, 180: 611-618. 10.1534/genetics.108.088575.PubMed CentralView ArticlePubMed
- Pryce JE, Arias J, Bowman PJ, Davis SR, Macdonald KA, Waghorn GC, Wales WJ, Williams YJ, Spelman RJ, Hayes BJ: Accuracy of genomic predictions of residual feed intake and 250-day body weight in growing heifers using 625,000 single nucleotide polymorphism markers. J Dairy Sci. 2012, 95: 2108-2119. 10.3168/jds.2011-4628.View ArticlePubMed
- Habier D, Tetens J, Seefried F-R, Lichtner P, Thaller G: The impact of genetic relationship information on genomic breeding values in German Holstein cattle. Genet Sel Evol. 2010, 42: 5-10.1186/1297-9686-42-5.PubMed CentralView ArticlePubMed
- Clark SA, Hickey JM, Daetwyler HD, Van der Werf JHJ: The importance of information on relatives for the prediction of genomic breeding values and implications for the makeup of reference populations in livestock breeding schemes. Genet Sel Evol. 2012, 44: 4-10.1186/1297-9686-44-4.PubMed CentralView ArticlePubMed
- Pszczola M, Strabel T, Mulder HA, Calus MP: Reliability of direct genomic values for animals with different relationships within and to the reference population. J Dairy Sci. 2012, 95: 389-400. 10.3168/jds.2011-4338.View ArticlePubMed
- Luan T, Woolliams JA, Lien S, Kent M, Svendsen M, Meuwissen THE: The accuracy of genomic selection in Norwegian Red cattle assessed by cross-validation. Genetics. 2009, 183: 1119-1126. 10.1534/genetics.109.107391.PubMed CentralView ArticlePubMed
- Lee SH, van der Werf JHJ, Hayes BJ, Goddard ME, Visscher PM: Predicting unobserved phenotypes for complex traits from whole-genome SNP data. PLoS Genet. 2008, 4: e1000231-10.1371/journal.pgen.1000231.PubMed CentralView ArticlePubMed
- Fisher RA: Frequency distribution of the values of the correlation coefficient in samples from an indefinitely large population. Biometrika. 1915, 10: 507-521.
- Henderson CR: Best linear unbiased estimation and prediction under a selection model. Biometrics. 1975, 31: 423-447. 10.2307/2529430.View ArticlePubMed
- Daetwyler HD, Kemper KE, van der Werf JHJ, Hayes BJ: Components of the accuracy of genomic prediction in a multi-breed sheep population. JAnim Sci. 2012, 90: 3375-3384.
- van der Werf JHJ, Kinghorn BP, Banks RG: Design and role of an information nucleus in sheep breeding programs. Anim Prod Sci. 2010, 50: 998-1003. 10.1071/AN10151.View Article
- White JD, Allingham PG, Gorman CM, Emery DL, Hynd P, Owens J, Bell A, Siddell J, Harper G, Hayes BJ, Daetwyler HD, Usmar J, Goddard ME, Henshall JM, Dominik S, Brewer H, van der Werf JHJ, Nicholas FW, Warner R, Hofmyer C, Longhurst T, Fisher T, Swan P, Forage R, Oddy VH: Design and phenotyping procedures for recording wool, skin, parasite resistance, growth, carcass yield and quality traits of the SheepGENOMICS mapping flock. Anim Prod Sci. 2012, 52: 157-171. 10.1071/AN11085.View Article
- Gardner GE, Williams A, Siddell J, Ball AJ, Mortimer S, Jacob RH, Pearce KL, Hocking Edwards JE, Rowe JB, Pethick DW: Using Australian sheep breeding values to increase lean meat yield percentage. Anim Prod Sci. 2010, 50: 1098-1106. 10.1071/AN10144.View Article
- Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006, 78: 629-644. 10.1086/502802.PubMed CentralView ArticlePubMed
- Browning BL, Browning SR: A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. AmJ Hum Genet. 2009, 84: 210-223. 10.1016/j.ajhg.2009.01.005.View Article
- Gilmour AR, Gogel B, Cullis BR, Thompson R: 2009 ASReml user guide release 3.0. 2009, VSN International Ltd, Hemel Hempstead
- Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR, Madden PA, Heath AC, Martin NG, Montgomery GW, Goddard ME, Visscher PM: Common SNPs explain a large proportion of the heritability for human height. Nat Genet. 2010, 42: 565-569. 10.1038/ng.608.PubMed CentralView ArticlePubMed
- Erbe M, Hayes BJ, Matukumalli LK, Goswani S, Bowman PJ, Reich CM, Mason BA, Goddard ME: Improving accuracy of genomic predictions within and between dairy cattle breeds with high density SNP panels. J Dairy Sci. 2012, 95: 4114-4129. 10.3168/jds.2011-5019.View ArticlePubMed
- Kijas JW, Lenstra JA, Hayes B, Boitard S, Porto Neto LR, San Cristobal M, Servin B, McCulloch R, Whan V, Gietzen K, Paiva S, Barendse W, Ciani E, Raadsma H, McEwan J, Dalrymple B: International Sheep Genomics Consortium: Genome-wide analysis of the World's sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol. 2012, 10: e1001258-10.1371/journal.pbio.1001258.PubMed CentralView ArticlePubMed
- Daetwyler HD, Pong-Wong R, Villanueva B, Woolliams JA: The impact of genetic architecture on genome-wide evaluation methods. Genetics. 2010, 185: 1021-1031. 10.1534/genetics.110.116855.PubMed CentralView ArticlePubMed
- Daetwyler HD, Villanueva B, Woolliams JA: Accuracy of predicting the genetic risk of disease using a genome-wide approach. PLoS ONE. 2008, 3: e3395-10.1371/journal.pone.0003395.PubMed CentralView ArticlePubMed
- Goddard ME: Genomic selection: prediction of accuracy and maximisation of long term response. Genetica. 2009, 136: 245-257. 10.1007/s10709-008-9308-0.View ArticlePubMed
- Cohen-Zinder M, Seroussi E, Larkin DM, Loor JJ, Everts-van der Wind A, Lee JH, Drackley JK, Band MR, Hernandez AG, Shani M, Lewin HA, Weller JI, Ron M: Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle. Genome Res. 2005, 15: 936-944. 10.1101/gr.3806705.PubMed CentralView ArticlePubMed
- Duncan EJ, Dodds KG, Henry HM, Thompson MP, Phua SH: Cloning, mapping and association studies of the ovine ABCG2 gene with facial eczema disease in sheep. Anim Genet. 2007, 38: 126-131. 10.1111/j.1365-2052.2006.01557.x.View ArticlePubMed
- O'Brien JK, Catt SL, Ireland KA, Maxwell WM, Evans G: In vitro and in vivo developmental capacity of oocytes from prepubertal and adult sheep. Theriogenology. 1997, 47: 1433-1443. 10.1016/S0093-691X(97)00134-9.View ArticlePubMed
- Meuwissen THE: Maximizing the response of selection with a predefined rate of inbreeding. J Anim Sci. 1997, 75: 934-940.PubMed
- Grundy B, Villanueva B, Woolliams JA: Dynamic selection procedures for constrained inbreeding and their consequences for pedigree development. Genet Res. 1998, 72: 159-168. 10.1017/S0016672398003474.View Article
- Sonesson AK, Woolliams JA, Meuwissen TH: Genomic selection requires genomic control of inbreeding. Genet Sel Evol. 2012, 44: 27-10.1186/1297-9686-44-27.PubMed CentralView ArticlePubMed
- Pryce JE, Hayes BJ, Goddard ME: Novel strategies to minimize progeny inbreeding while maximizing genetic gain using genomic information. J Dairy Sci. 2012, 95: 377-388. 10.3168/jds.2011-4254.View ArticlePubMed
- Meuwissen T, Goddard M: Accurate prediction of genetic values for complex traits by whole-genome resequencing. Genetics. 2010, 185: 623-631. 10.1534/genetics.110.116590.PubMed CentralView ArticlePubMed
- Ober U, Ayroles JF, Stone EA, Richards S, Zhu D, Gibbs RA, Stricker C, Gianola D, Schlather M, Mackay TFC, Simianer H: Using whole-genome sequence data to predict quantitative trait phenotypes in Drosophila melanogaster. PLoS Genet. 2012, 8: e1002685-10.1371/journal.pgen.1002685.PubMed CentralView ArticlePubMed
- Daetwyler HD: Genome-wide evaluation of populations. 2009, PhD thesis. Wageningen University, ISBN: 978-90-8585-528-6
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.