Study of the optimum haplotype length to build genomic relationship matrices.
BACKGROUND
As genomic data becomes more abundant, genomic prediction is more routinely used to estimate breeding values. In genomic prediction, the relationship matrix ([Formula: see text]), which is traditionally used in genetic evaluations is replaced by the genomic relationship matrix ([Formula: see text]). This paper considers alternative ways of building relationship matrices either using single markers or haplotypes of different lengths. We compared the prediction accuracies and log-likelihoods when using these alternative relationship matrices and the traditional [Formula: see text] matrix, for real and simulated data.


METHODS
For real data, we built relationship matrices using 50k genotype data for a population of Brahman cattle to analyze three traits: scrotal circumference (SC), age at puberty (AGECL) and weight at first corpus luteum (WTCL). Haplotypes were phased with hsphase and imputed with BEAGLE. The relationship matrices were built using three methods based on haplotypes of different lengths. The log-likelihood was considered to define the optimum haplotype lengths for each trait and each haplotype-based relationship matrix.


RESULTS
Based on simulated data, we showed that the inverse of [Formula: see text] matrix and the inverse of the haplotype relationship matrices for methods using one-single nucleotide polymorphism (SNP) phased haplotypes provided coefficients of determination (R2) close to 1, although the estimated genetic variances differed across methods. Using real data and multiple SNPs in the haplotype segments to build the relationship matrices provided better results than the [Formula: see text] matrix based on one-SNP haplotypes. However, the optimal haplotype length to achieve the highest log-likelihood depended on the method used and the trait. The optimal haplotype length (7 to 8 SNPs) was similar for SC and AGECL. One of the haplotype-based methods achieved the largest increase in log-likelihood for SC, i.e. from -1330 when using [Formula: see text] to -1325 when using haplotypes with eight SNPs.


CONCLUSIONS
Building the relationship matrix by using haplotypes that comprise multiple SNPs will increase the accuracy of estimated breeding values. However, the optimum haplotype length that shows the correct relationship among individuals for each trait can be derived from the data.
The Centre for Genetic Analysis and Applications, School of Environmental and Rural Science, University of New England, Armidale, Australia. mferdosi@myune.edu.au.|Animal Genetics and Breeding Unit, University of New England, Armidale, Australia. mferdosi@myune.edu.au.
Mohammad H
MH
Ferdosi
Cobb-Vantress, Siloam Springs, AR, USA.|CSIRO Agriculture Flagship, FD McMaster Laboratory Chiswick, Armidale, Australia.
John
J
Henshall
Animal Genetics and Breeding Unit, University of New England, Armidale, Australia. btier@une.edu.au.
Bruce
B
Tier
J Dairy Sci. 2008 Nov;91(11):4414-23
18946147
Nat Genet. 2010 Jul;42(7):565-9
20562875
J Anim Sci. 2012 May;90(5):1398-410
22100599
Genet Sel Evol. 2011 Jul 12;43:1-10
21749687
Biometrics. 1975 Jun;31(2):423-47
1174616
Genetics. 2008 Jan;178(1):553-61
18202394
Genet Sel Evol. 2011 Jan 05;43:1
21208445
Genet Sel Evol. 2011 Jun 26;43:25
21703021
Genet Sel Evol. 2014 Aug 01;46:49
25158793
J Anim Breed Genet. 2009 Feb;126(1):3-13
19207924
J Anim Breed Genet. 2006 Aug;123(4):218-23
16882088
J Anim Breed Genet. 2013 Aug;130(4):259-69
23855628
J Anim Sci. 1997 Jul;75(7):1738-45
9222829
Animal. 2010 Feb;4(2):157-64
22443868
Genetics. 2013 Jul;194(3):597-607
23640517
Theor Popul Biol. 2016 Feb;107:26-30
26341159
Genet Sel Evol. 2012 Aug 31;44:28
22937985
PLoS Genet. 2010 Nov 04;6(11):e1001188
21079680
Am J Hum Genet. 2007 Nov;81(5):1084-97
17924348
BMC Bioinformatics. 2014 Jun 07;15:172
24906803
Genet Sel Evol. 2014 Feb 04;46:11
24495596