The importance of information on relatives for the prediction of genomic breeding values and the implications for the makeup of reference data sets in livestock breeding schemes
© Clark et al; licensee BioMed Central Ltd. 2012
Received: 19 July 2011
Accepted: 9 February 2012
Published: 9 February 2012
The theory of genomic selection is based on the prediction of the effects of genetic markers in linkage disequilibrium with quantitative trait loci. However, genomic selection also relies on relationships between individuals to accurately predict genetic value. This study aimed to examine the importance of information on relatives versus that of unrelated or more distantly related individuals on the estimation of genomic breeding values.
Simulated and real data were used to examine the effects of various degrees of relationship on the accuracy of genomic selection. Genomic Best Linear Unbiased Prediction (gBLUP) was compared to two pedigree based BLUP methods, one with a shallow one generation pedigree and the other with a deep ten generation pedigree. The accuracy of estimated breeding values for different groups of selection candidates that had varying degrees of relationships to a reference data set of 1750 animals was investigated.
The gBLUP method predicted breeding values more accurately than BLUP. The most accurate breeding values were estimated using gBLUP for closely related animals. Similarly, the pedigree based BLUP methods were also accurate for closely related animals, however when the pedigree based BLUP methods were used to predict unrelated animals, the accuracy was close to zero. In contrast, gBLUP breeding values, for animals that had no pedigree relationship with animals in the reference data set, allowed substantial accuracy.
An animal's relationship to the reference data set is an important factor for the accuracy of genomic predictions. Animals that share a close relationship to the reference data set had the highest accuracy from genomic predictions. However a baseline accuracy that is driven by the reference data set size and the overall population effective population size enables gBLUP to estimate a breeding value for unrelated animals within a population (breed), using information previously ignored by pedigree based BLUP methods.
Genomic selection (GS) is a method that uses genomic information to estimate breeding values and rank selection candidates in livestock breeding programs. It has become widely used in some livestock industries e.g. dairy cattle and pig improvement programs. Initial studies on genomic evaluation have suggested that GS predicts the effects of markers in linkage disequilibrium (LD) with quantitative trait loci (QTL). This implies that accurate predictions of breeding value may persist for several generations, allowing for: 1) a reduced number of phenotypic measurements in each generation  and; 2) the possibility of accurate predictions across different breeds provided sufficient marker density . Habier et al.  proposed that genomic predictions also rely on the genetic relationships between individuals with phenotypic records, usually known as the reference data set, and those whose breeding value is to be predicted [4, 5]. The following question arises: does an animal that has its breeding value predicted from genomic information require relatives in a reference data set?
The reference data set is used to gain information on important phenotypes and genotypes so that genomic estimated breeding values (GEBV) can be highly accurate for selection candidates. The makeup and size of this data set, combined with the methods used to predict the breeding value, govern the accuracy achieved in many breeding schemes . Due to the cost of measuring genotypes and phenotypes on large numbers of individuals, it has been suggested that using a specially selected reference data set may be a cost effective way of gaining the economic advantage presented by genomic selection, especially for species such as beef cattle and sheep that do not have the nucleus structure of the dairy cattle, pig and poultry industries [7, 8].
Various methods are used to predict breeding values from genomic data. These range from variable selection methods such as BayesB, which allows only a small number of loci to have an effect, some of them potentially large, to gBLUP, which assumes equal variance across all loci . Empirical evidence across livestock populations has shown that in many cases these methods obtain very similar accuracies of the estimated breeding value . This suggests that additive genetic variation for many traits is controlled by many genes with a small effect, somewhat like Fisher's (1918)  'infinitesimal model'.
The gBLUP method to estimate genomic breeding values has been widely described [5, 10, 12, 13]. This method uses genomic information in the form of a genomic relationship matrix (GRM) that defines the additive genetic covariance between animals . The GRM then replaces the pedigree-based numerator relationship matrix (NRM) in the traditional BLUP equations. The GRM is expected to give a more accurate estimate of the covariance between individuals, however, it is important to understand how much accuracy is gained from improved measures of covariance among known relatives and how much is gained from information on distant 'relatives' previously ignored via the pedigree method.
Goddard  and Daetwyler et al.  presented methods to predict the accuracy of genomic selection. These methods predict accuracy based on the effective population size, the number of animals in the reference data set, the heritability and the effective number of chromosome segments segregating within the population. Hayes et al.  extended the approach by Goddard  to consider information from relatives and illustrated that the GRM uses information on true relationships, which can vary around the expected additive genetic relationship predicted from pedigree. This study will compare the accuracy of an animal's breeding value that has a strong pedigree relationship with a reference data set with that of an animal that is essentially unrelated to the reference data set, and discuss the effect of these relationships on the design of reference data sets used in genomic selection breeding schemes.
To examine the effect of relationships between animals in the test and reference data sets, we used both computer simulation and real phenotypic data from the Australian Sheep CRC Information Nucleus Flock (INF). The INF animals are located at eight sites across Australia and managed by CRC partner organizations, including The University of New England, NSW Department of Primary Industries, Victorian Department of Primary Industries, South Australian Research and Development Institute and The Department of Agriculture of Western Australia. The experimental data in this paper were obtained according to protocols approved by the Animal Ethics committees of these organizations.
Genotype simulations were conducted using the Markovian Coalescence Simulator (MaCS)  to simulate 2 000 base haplotypes, with an effective population size (N e ) of 100. As described in Clark et al. , thirty chromosomes each with base haplotypes of a 100 cM region (1·108 base pairs) were simulated, with a per site mutation rate of 2.5·10-8. The total number of SNP segregating on the genome was approximately 1 670 000 (SNP sequence). Sixty thousand SNP markers and 10 000 QTL were randomly selected from the SNP sequence in the base generation to be used in the genomic analysis (9428 QTL segregating in the final generation). Therefore each SNP had a 3% chance of being used as a marker and a 0.5% chance of being used as a QTL. The additive effect of each QTL was drawn from a gamma distribution with a shape and scale of 0.4 and 1.66 respectively  and had a 50% chance of being positive or negative.
The base population haplotypes were randomly allocated to 80 base males and 2 000 base female animals of a simulated population structure, with 10 subsequent generations receiving these haplotypes via mendelian inheritance, allowing recombination to occur according to genetic distance, i.e. 1% recombination per cM. The population was simulated for 10 generations and each generation contained 4000 animals, half male and half female. Eighty males were randomly selected in each generation and each male was randomly mated to 25 females, which each had two offspring per generation. Only breeding animals were allocated breeding values and phenotypes.
where β j is the additive effect of QTL genotype (j) and Q kj is the QTL genotype at locus j which is coded as 0, 1, or 2 and is the number of copies of the QTL that an individual (k) carries. Trait phenotypes were simulated based on a heritability (h2) of 0.3 and all other variation in phenotype was due to a random environmental effect drawn from a normal distribution with variance .
Three reference data sets of 1 750 animals were formed for the simulation study. Reference data set 1 (closely related) consisted of animals from the 10th generation. This reference data set was constructed such that animals in the test data set had 20 half sibs in the reference data set. Reference data set 2 (distantly related) also consisted of animals in the 10th generation but there were no close relationships between animals in the test and reference data sets. However there were some second degree relationships (1st cousins) between the two data sets. The final reference data set (unrelated) consisted of females from generation 1 and resulted in a very low or zero relationship between the two data sets. The accuracy of prediction was assessed in the test data set which consisted of 250 animals from the 10th generation and the average correlation between TBV and estimated breeding value (genomic or pedigree based) was calculated over 10 replicates of the simulation study.
where y is a vector of phenotypes, μ is the mean, 1n is a vector of 1s, Z is a design matrix allocating records to breeding values, g is a vector of breeding values for animals in the reference set and the test set and e is a vector of normal deviates with variance . Furthermore where G is the genomic relationship matrix (GRM), and is the genetic variance for this model. The GRM (G) was formed using the method as defined by VanRaden 
Traditional best linear unbiased prediction (BLUP) was also performed, using a deep (BLUP-D), 10 generation pedigree or a shallow, single generation pedigree (BLUP-S). Traditional BLUP ignores genomic data and relies on information from ancestors using a numerator relationship matrix (A) based on pedigree information. This method uses the same model as gBLUP (above) however with the vector of additive genetic values g replaced by a, with where A is the numerator relationship matrix and is the additive genetic variance. Variance components for both BLUP methods were estimated with ASREML  and the model solutions yielded estimated breeding values.
Merino sheep phenotypic data
The reference data set consisted of phenotypic and genotypic records for the Merino sheep from the Australian Sheep Cooperative Research Centre information nucleus flock (INF) . The traits ultrasound scanned eye muscle depth (EMD; 1781 animals) and live weight at ultrasound scanning (SC_WT; 1743 animals) were evaluated. Scanned EMD is used to estimate the size of the rib-eye muscle, which produces high value cuts of meat, and SC_WT is highly correlated to an animal's weight at post weaning. Animals in the INF were sired by rams from the wider Merino population; these sires were chosen to maximize the connectedness with the Australian sheep flock by sampling artificial insemination sires from a wide range of sheep breeders.
The test data set consisted of a population of Australian Merino industry sires with highly accurate Australian sheep breeding values (ASBV). Information about ASBV definitions can be found at the following website maintained by Australian Wool Innovation Ltd and Meat and Livestock Australia . The industry sires were divided into closely, distantly and unrelated groups based on their pedigree relationship the animals in the INF flock. The maximum relationship of an animal in the test data set with an animal in the reference data set ranged from 0.125 to 0.5 (no progeny included) for the 48 closely related test animals, from 0 to 0.125 for the 60 distantly related test animals, and the 53 unrelated test animals shared no pedigree relationship to the reference data set.
All animals were genotyped using the Illumina 50 K ovine SNP chip (Illumina Inc., San Diego, CA, USA), which includes 54 977 SNP. Following the genotyping procedures, quality control measures were applied to all SNP as follows: SNP were removed if they had a call rate of less than 95%, a GC score (proportion of guanine-cytosine pairs) of less than 0.6, a minor allele frequency of less than 0.01, a SNP heterozygosity of greater than 3 s.d. from the mean (mean heterozygosity, 0.374; s.d., 0.129), were not in Hardy-Weinberg equilibrium (a P-value cut-off of 1·10-15), had no genome location or were in greater than 0.99 LD with another SNP on the chip . After these quality control measures were applied, 48 640 SNP were used. Missing genotypes were imputed using fast PHASE .
The following fixed effects were fitted in both trait models: sex, birth type, rearing type, age of dam, contemporary group (birth year • birth month • site • management group) and age-at-trait recording, SC_WT was fitted in the analysis of EMD.
where X is a design matrix relating the fixed effects (as described above) to each animal and b is a vector of fixed effects. Genetic evaluation was undertaken using BLUP-D, using a deep pedigree of 7277 animals that ranged from one to eight generations in length (depending on the individual), and gBLUP as defined above. Variance components for all methods were again estimated using ASREML  and the model solutions yielded estimated breeding values.
Validation and accuracy
The empirical accuracy (r(cor)) for the Merino data analysis was evaluated as the Pearson product-moment correlation between the GEBV and a progeny test ASBV of the animals in the test data set. The empirical accuracy may be an underestimate of the 'real' accuracy because the ASBV accuracies are below 1. The validation sires had an ASBV accuracy greater than 0.5 and the mean accuracy was 0.85 for EMD and 0.9 for SC_WT. The ASBVs used were calculated such that they included no information from animals in the reference data set. The Merino sheep population is highly heterogeneous and can be divided into strains defined as fine, medium and strong wool types. Correlations between GEBV and ASBV were calculated after accounting for the effect of strain.
The empirical accuracy of the breeding values estimated in the test set, for the simulation example was defined as the correlation between the true and estimated breeding value. The accuracy was also estimated for each individual as: where; PEV is the prediction error variance estimated using elements from the mixed model equations, G ii is the diagonal of the GRM for animal i and is substituted for A ii in traditional BLUP, is the additive genetic variance. Furthermore, where; C ii is the diagonal of inverse of the coefficient matrix for animal i and is the residual variance (See Appendix 1).
To determine the effect of an individual's relationship to the reference data set on the accuracy of genomic predictions, a range of comparisons were made between varying definitions of relatedness and an individual's GEBV accuracy using r(PEV). Four measures of genomic relatedness were considered: a) An animal's mean relationship with the reference data set; b) its maximum relationship; c) its mean top 10 relationships and d) its mean top 100 relationships.
Empirical accuracy1 (± S.E.)2 using genomic and pedigree based methods in simulated data.
Relationship to reference [pedigree relationship]
There was no significant difference in accuracy between BLUP-S and BLUP-D when the animals in the test and reference data sets had a close relationship. However, when a shallow pedigree was used, and animals in the test and reference data sets were distantly related or unrelated, all breeding values estimated using BLUP-S were zero. In contrast, BLUP-D predicted a breeding value with a significant accuracy when the reference and test data sets shared a distant relationship and accuracy reduced to close to zero when animals in the reference and test data sets were unrelated.
Merino sheep data analysis
Empirical1 (r(cor)) and estimated accuracy2 (r(pev)) using genomic and pedigree based methods for the Merino EMD data.
Relationship to reference [Maximum pedigree relationship]
Empirical1 (r(cor)) and estimated accuracy2 (r(pev)) using genomic and pedigree based methods for the Merino SC_WT data.
Relationship to reference [Maximum pedigree relationship]
For the EMD example, the estimated and empirical accuracies were very similar when using gBLUP. However for the SC_WT scenarios, there was a large difference between the estimated accuracy and the empirical accuracy.
The relationship between the animals in the test and reference data sets has an effect on the accuracy of genomic predictions. Close relationships between the two data sets' result in the highest accuracy for GEBV. Similar results were predicted by Hayes et al.  and observed by Habier et al. [3, 4] for populations that share a close relationship. However, breeding values that are predicted for closely related animals using the traditional pedigree-based BLUP approach also achieve high accuracy. The current study has shown that when there is a distant relationship between the animals in the test and reference data sets, gBLUP is still able to predict an animal's breeding value with some accuracy. Furthermore, when the animals are unrelated by pedigree or when the pedigree relationships are low, gBLUP can use information from distant relatives to maintain a proportion of accuracy of the GEBV.
The information gathered from only distantly related animals enabled an estimate of breeding value to be made with some accuracy. However, when relatives were included in the reference data set, the importance of information on distantly related animals may be reduced. Selection index theory shows that when information on closely related animals is available, more weight is placed on this information and therefore information from distantly related animals becomes less important. Although the importance of information from distant relatives is reduced, this extra information, which is not used in pedigree-based methods, enables gBLUP to achieve a higher accuracy of the EBV. The inclusion of information on relatives improves the accuracy of the predicted breeding values.
Where ρ = (1 + a+ 2√a), with a = 1+2*λ/N and N is the number of animals in the reference,
λ = σ2e /σ2u where σ2e is the residual variance and σ2u is the genetic variance at a single locus and is estimated by σ2u = h2/M e ·k where M e = 2N e L and is the effective number of chromosome segments, h2 is the heritability and k = 1/log(2N e ). For the simulation example N = 1750, N e = 100, h2 = 0.3 and L = 30. Then k = 0.189, λ = 3773.8, a = 5.31 and consequently the accuracy for an individual with no phenotype was equal to 0.36. Similarly, the alternative method described by Daetwyler et al.  results in a predicted accuracy of 0.28 (details not shown). The predicted accuracies resulting from either method were similar to the baseline accuracy in our study achieved by gBLUP in unrelated individuals (0.34). In the theoretical prediction methods, there is some ambiguity about the approximation of M e [5, 23], with proposed values equal to: a) 2N e L/ln( 4N e L); b) 4N e L and c) 2N e L. Using  for each of these values results in predicted accuracies of a) 0.74 b) 0.27 and c) 0.36. Consequently 2N e L appears to be the most appropriate variable for baseline accuracy in our simulation example. For the Merino sheep data, with an estimated N e of approximately 1,000 , the expected accuracy was 0.15 and lower than that achieved by gBLUP for EMD (0.28) and for SC_WT (0.18). This increase for gBLUP in the real data is possibly due to extra information from animals that shared a genomic relationship but were unknown in the pedigree, or the estimation of N e may have been affected by heterogeneity of the breed, which really consists of several sub-populations.
Accuracy estimated using the prediction error variance of the mixed model equations (r(PEV)) was shown to be a good approximation of empirical accuracy for the simulation example. Estimated and empirical accuracies were also very similar when using gBLUP for the EMD example. However, some differences between r(PEV)and empirical accuracy were observed for both, BLUP-D and gBLUP in real data in the case of SC_WT. In the simulation example, the empirical accuracy was the correlation between the TBV and EBV (or GEBV), whereas in the Merino data example, the empirical accuracy was the correlation between the ASBV and EBV (or GEBV). The ASBVs are progeny test estimates and have some prediction error associated with them. The empirical accuracy was also likely to be affected by sampling because of the small size of each test data set (50-60 animals). Furthermore, unlike the simulation data, where all animals were linked by a true pedigree, many Merino animals in the unrelated test set had no direct pedigree relationships with the reference data set and therefore only zero breeding values were estimated for these animals. In contrast, in the case of missing pedigree, gBLUP could use genomic relationship information and a more accurate breeding value was estimated for all animals in the test set.
Another complexity in our real data example is the heterogeneity of the Merino sheep population, as it consists of many sub-populations. In routine ASBV analyses, this population structure is accounted for using pedigree information and genetic groups based on individual flock data. When correlating GEBV and ASBV, we accounted for sub-population effects by assigning sires to groups of "fine wool", "medium wool" and "strong wool". Empirical accuracies for SC_WT were clearly affected by correcting for the sub-population structure, which may explain why there are some differences between r( PEV)and r(cor) for this trait. The corrections had little to no effect on empirical accuracy for EMD. Note that EMD was corrected for SC_WT and this may have removed some of the sub-population effects on EMD.
The makeup of reference data sets is an important factor for the design of genomic evaluation systems to enable additional genetic gain from genomic selection at the lowest cost. This is especially true for beef cattle and sheep breeding programs that do not have a distinct nucleus tier. We have shown that genomic predictions are more accurate when animals are related to the reference data set; however substantial baseline accuracy can be achieved for all animals in the population. To achieve this, the reference data set will need to include a large number of animals that cover the genetic diversity of the given population (breed). It may be important to include animals that are expected to contribute more to the future gene pool in that breed but these contributions need to be balanced by contributions to genetic diversity .
The optimal size of the reference data set will depend on N e of the given population; populations with higher N e may need a larger reference data set so that suitable baseline accuracies can be achieved. If the baseline accuracy is low (large N e and small reference data set size) the contribution of relatives' information will be larger, however this information from relatives is only limited to closely related individuals and will not last over many generations.
The relationship between animals in the reference and test data sets affects the accuracy of predicting breeding values using gBLUP. When there is a close relationship between the animals in the reference and test data sets, gBLUP can estimate breeding values with a high accuracy. When there is only a distant relationship between the animals in test and reference data sets, gBLUP can still estimate a breeding value with some accuracy. This baseline accuracy depends on the effective population size and the size of the reference data set, and should be carefully considered when designing a reference data set for a breeding program.
Appendix 1: Accuracy estimated using the PEV of the mixed model equations weighted by genomic relationships
where is the residual variance and C ii is the diagonal of the inverse of the coefficient matrix, furthermore:
The authors acknowledge the Cooperative Research Centre for Sheep Industry Innovation for funding. We thank Daniel Brown for providing validation ram ASBV excluding INF data, Klint Gore for managing the phenotypic INF database, Cedric Gondro for performing part of the genotype quality control for the INF data and the many staff involved at the INF sites across Australia. We would also like to thank Ben Hayes, Brian Kinghorn, Andrew Swan and Bruce Tier for their valued comments and discussions regarding this research.
- Muir WM: Comparison of genomic and traditional BLUP-estimated breeding value accuracy and selection response under alternative trait and genomic parameters. J Anim Breed Genet. 2007, 124: 342-355. 10.1111/j.1439-0388.2007.00700.x.View ArticlePubMed
- Goddard ME, Hayes BJ, McPartlan H, Chamberlain AJ: Can the same genetic markers be used in multiple breeds?. Proceedings of the 8th World Congress on Genetics Applied to Livestock Production: August 13-18. 2006, , Brazil. CD-ROM communication no. 22-16
- Habier D, Fernando RL, Dekkers JCM: The impact of genetic relationship information on genome-assisted breeding values. Genetics. 2007, 177: 2389-2397.PubMed CentralPubMed
- Habier D, Tetens J, Seefried FR, 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
- Hayes BJ, Visscher PM, Goddard ME: Increased accuracy of artificial selection by using the realized relationship matrix. Genet Res. 2009, 91: 47-60. 10.1017/S0016672308009981.View Article
- 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
- Goddard ME, Hayes BJ: Genomic Selection. J Anim Breed Genet. 2007, 124: 323-330. 10.1111/j.1439-0388.2007.00702.x.View ArticlePubMed
- 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
- Meuwissen THE, Hayes BJ, Goddard ME: Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001, 157: 1819-1829.PubMed CentralPubMed
- Moser G, Tier B, Crump RE, Khatkar MS, Raadsma HW: A comparison of five methods to predict genomic breeding values of dairy bulls from genome-wide SNP markers. Genet Sel Evol. 2009, 41: 56-10.1186/1297-9686-41-56.PubMed CentralView ArticlePubMed
- Fisher RA: The correlation between relatives on the supposition of mendelian inheritance. Trans R Soc Edin. 1918, 52: 399-433.View Article
- VanRaden PM, Van Tassell CP, Wiggans GR, Sonstegard TS, Schnabel RD, Taylor JF, Schenkel F: 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
- Hayes BJ, Bowman PJ, Chamberlain AC, 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
- Nejati-Javaremi A, Smith C, Gibson JP: Effect of total allelic relationship on accuracy of evaluation and response to selection. J Anim Sci. 1997, 75: 1738-1745.PubMed
- 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
- Chen GK, Marjoram P, Wall JD: Fast and flexible simulation of DNA sequence data. Genome Res. 2009, 19: 136-142.PubMed CentralView ArticlePubMed
- Clark SA, Hickey JM, van der Werf JHJ: Different models of genetic variation and their effect on genomic evaluation. Genet Sel Evol. 2011, 43: 18-10.1186/1297-9686-43-18.PubMed CentralView ArticlePubMed
- VanRaden PM: Efficient methods to compute genomic predictions. J Dairy Sci. 2008, 91: 4414-4423. 10.3168/jds.2007-0980.View ArticlePubMed
- Gilmour AR, Gogel BJ, Cullis BR, Thompson R: ASReml User Guide Release 3.0. Hemel Hempstead: VSN International Ltd;. 2009
- 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
- AWI, MLA: Gain from genetics. 2010, Australian Wool and Innovation, Meat and Livestock Australia, [http://www.makingmorefromsheep.com.au/gain-from-genetics/index.htm]
- 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
- Goddard ME, Hayes BJ, Meuwissen THE: Genomic selection in livestock populations. Genet Res. 2010, 92: 413-421. 10.1017/S0016672310000613.View Article
- Kijas JW, Townley D, Dalrymple BP, Heaton MP, Maddox JF, McGrath A, Wilson P, Ingersoll RG, McCulloch R, McWilliam S, Tang D, McEwan J, Cockett N, Oddy VH, Nicholas FW, Raadsma H: A genome wide survey of SNP variation reveals the genetic structure of sheep breeds. PLoS One. 2009, 4 (3): e4668-10.1371/journal.pone.0004668.PubMed CentralView ArticlePubMed
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.