Rapid screening for phenotype-genotype associations by linear transformations of genomic evaluations

BMC Bioinformatics

Jose L Gualdrón Duarte²  Rodolfo JC Cantet²  Ronald O Bates¹  Catherine W Ernst¹  Nancy E Raney¹  Juan P Steibel³ 

[1] Department of Animal Science, Michigan State University, East Lansing, MI, USA

[2] Departamento de Producción Animal, Facultad de Agronomía, UBA-CONICET, Buenos Aires, Argentina

[3] Department of Fisheries and Wildlife, Michigan State University, East Lansing, MI, USA

关键词: Pig genotype;    Marker variance;    Genome wide association;   

Others  :  1087546
DOI  :  10.1186/1471-2105-15-246

 received in 2014-02-03, accepted in 2014-07-07,  发布年份 2014

【 摘 要 】

Background

Currently, association studies are analysed using statistical mixed models, with marker effects estimated by a linear transformation of genomic breeding values. The variances of marker effects are needed when performing the tests of association. However, approaches used to estimate the parameters rely on a prior variance or on a constant estimate of the additive variance. Alternatively, we propose a standardized test of association using the variance of each marker effect, which generally differ among each other. Random breeding values from a mixed model including fixed effects and a genomic covariance matrix are linearly transformed to estimate the marker effects.

Results

The standardized test was neither conservative nor liberal with respect to type I error rate (false-positives), compared to a similar test using Predictor Error Variance, a method that was too conservative. Furthermore, genomic predictions are solved efficiently by the procedure, and the p-values are virtually identical to those calculated from tests for one marker effect at a time. Moreover, the standardized test reduces computing time and memory requirements.

The following steps are used to locate genome segments displaying strong association. The marker with the highest − log(p-value) in each chromosome is selected, and the segment is expanded one Mb upstream and one Mb downstream of the marker. A genomic matrix is calculated using the information from those markers only, which is used as the variance-covariance of the segment effects in a model that also includes fixed effects and random genomic breeding values. The likelihood ratio is then calculated to test for the effect in every chromosome against a reduced model with fixed effects and genomic breeding values. In a case study with pigs, a significant segment from chromosome 6 explained 11% of total genetic variance.

Conclusions

The standardized test of marker effects using their own variance helps in detecting specific genomic regions involved in the additive variance, and in reducing false positives. Moreover, genome scanning of candidate segments can be used in meta-analyses of genome-wide association studies, as it enables the detection of specific genome regions that affect an economically relevant trait when using multiple populations.

【 授权许可】

   
2014 Gualdrón Duarte et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

【 参考文献 】

• [1]Crossa J, Pérez P, de los Campos G, Mahuku G, Dreisigacker S, Magorokosho C: Genomic selection and prediction in plant breeding. J Crop Improv 2011, 25:239-261.
• [2]Goddard ME, Hayes BJ: Mapping genes for complex traits in domestic animals and their use in breeding programmes. Nat Rev Genet 2009, 10:381-391.
• [3]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.
• [4]Hayes BJ, Pryce J, Chamberlain AJ, Bowman PJ, Goddard ME: Genetic architecture of complex traits and accuracy of genomic prediction: coat colour, milk-fat percentage, and type in Holstein cattle as contrasting model traits. PLoS Genet 2010, 6:e1001139.
• [5]Kumar S, Garrick DJ, Bink MC, Whitworth C, Chagné D, Volz RK: Novel genomic approaches unravel genetic architecture of complex traits in apple. BMC Genomics 2013, 14:393.
• [6]Zhou X, Stephens M: Genome-wide efficient mixed-model analysis for association studies. Nat Genet 2012, 44:821-824.
• [7]Meuwissen TH, Hayes BJ, Goddard ME: Prediction of total genetic value using genome-wide dense marker maps. Genetics 2001, 157:1819-1829.
• [8]Garrick DJ: Equivalent mixed model equations for genomic selection. J Bone Miner Res 2007, 90(Suppl):376. (Abstr.)
• [9]Strandén I, Garrick DJ: Technical note: derivation of equivalent computing algorithms for genomic predictions and reliabilities of animal merit. J Dairy Sci 2009, 92:2971-2975.
• [10]Edwards DB, Ernst CW, Tempelman RJ, Rosa GJM, Raney NE, Hoge MD, Bates RO: Quantitative trait loci mapping in an F2 Duroc x Pietrain resource population: I. Growth traits. J Anim Sci 2008, 86:241-253.
• [11]Edwards DB, Ernst CW, Raney NE, Doumit ME, Hoge MD, Bates RO: Quantitative trait locus mapping in an F2 Duroc x Pietrain resource population: II. Carcass and meat quality traits. J Anim Sci 2008, 86:254-266.
• [12]Ramos AM, Crooijmans RP, Affara NA, Amaral AJ, Archibald AL, Beever JE, Bendixen C, Churcher C, Clark R, Dehais P, Hansen MS, Hedegaard J, Hu Z-L, Kerstens HH, Law AS, Megens H-J, Milan D, Nonneman DJ, Rohrer GA, Rothschild MF, Smith TPL, Schnabel RD, Van Tassell CP, Taylor JF, Wiedmann RT, Schook LB, Groenen MA: Design of a high density SNP genotyping assay in the pig using SNPs identified and characterized by next generation sequencing technology. PLoS One 2009, 4:e6524.
• [13]Badke YM, Bates RO, Ernst CW, Schwab C, Fix J, Van Tassell CP, Steibel JP: Methods of tagSNP selection and other variables affecting imputation accuracy in swine. BMC Genet 2013, 14:8.
• [14]Badke YM, Bates RO, Ernst CW, Schwab C, Steibel JP: Estimation of linkage disequilibrium in four US pig breeds. BMC Genomics 2012, 13:24.
• [15]Gualdrón Duarte JL, Bates RO, Ernst CW, Raney NE, Cantet RJC, Steibel JP: Genotype imputation accuracy in a F2 pig population using high density and low density SNP panels. BMC Genet 2013, 14:38.
• [16]Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007, 81:559-575.
• [17]Hickey JM, Kinghorn BP, Tier B, van der Werf JH, Cleveland MA: A phasing and imputation method for pedigreed populations that results in a single-stage genomic evaluation. Genet Sel Evol 2012, 44:9.
• [18]VanRaden PM: Efficient methods to compute genomic predictions. J Dairy Sci 2008, 91:4414-4423.
• [19]Clifford D, McCullagh P: The regress function. R News 2006, 6:6-10.
• [20]McClure MC, Ramey HR, Rolf MM, McKay SD, Decker JE, Chapple RH, Kim JW, Taxis TM, Weaber RL, Schnabel RD, Taylor JF: Genome-wide association analysis for quantitative trait loci influencing Warner-Bratzler shear force in five taurine cattle breeds. Anim Genet 2012, 43:662-673.
• [21]Wang H, Misztal I, Aguilar I, Legarra A, Muir WM: Genome-wide association mapping including phenotypes from relatives without genotypes. Genet Res (Camb) 2012, 94:73-83.
• [22]Henderson C: Applications of Linear Models in Animal Breeding. Guelph: University of Guelph; 1984.
• [23]Yu J, Pressoir G, Briggs WH, Vroh Bi I, Yamasaki M, Doebley JF, McMullen MD, Gaut BS, Nielsen DM, Holland JB, Kresovich S, Buckler ES: A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet 2006, 38:203-208.
• [24]Kang HM, Zaitlen NA, Wade CM, Kirby A, Heckerman D, Daly MJ, Eskin E: Efficient control of population structure in model organism association mapping. Genetics 2008, 178:1709-1723.
• [25]Endelman JB: Ridge regression and other Kernels for genomic selection with R package rrBLUP. Plant Genome J 2011, 4:250.
• [26]Rangkasenee N, Murani E, Brunner RM, Schellander K, Cinar MU, Luther H, Hofer A, Stoll M, Witten A, Ponsuksili S, Wimmers K: Genome-wide association identifies TBX5 as candidate gene for Osteochondrosis providing a functional link to cartilage perfusion as initial factor. Front Genet 2013, 4:78.
• [27]Do DN, Ostersen T, Strathe AB, Mark T, Jensen J, Kadarmideen HN: Genome-wide association and systems genetic analyses of residual feed intake, daily feed consumption, backfat and weight gain in pigs. BMC Genet 2014, 15:27.
• [28]Fan Y, Xing Y, Zhang Z, Ai H, Ouyang Z, Ouyang J, Yang M, Li P, Chen Y, Gao J, Li L, Huang L, Ren J: A further look at porcine chromosome 7 reveals VRTN variants associated with vertebral number in Chinese and Western pigs. PLoS One 2013, 8:e62534.
• [29]Liang K-Y, Self SG: On the asymptotic behaviour of the Pseudolikelihood Ratio Test Statistic. J R Stat Soc Ser B 1996, 58:785-796.
• [30]Self SG, Liang K-Y: Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. J Am Stat Assoc 1987, 82:605-610.
• [31]Vaughan LK, Divers J, Padilla M, Redden DT, Hemant K, Pomp D, Allison DB: The use of plasmodes as a supplement to simulations: a simple example evaluating individual admixture estimation methodologies. Comput Stat Data Anal 2009, 53:1755-1766.
• [32]Klein RJ, Zeiss C, Chew EY, Tsai J, Sackler RS, Haynes C, Henning AK, Sangiovanni JP, Mane SM, Susan T, Bracken MB, Ferris FL, Ott J, Barnstable C, Hoh J: Complement factor H Polymorphism in age-related macular degeneration. Science (80-) 2006, 308:385-389.
• [33]Legarra A, Aguilar I, Misztal I: A relationship matrix including full pedigree and genomic information. J Dairy Sci 2009, 92:4656-4663.
• [34]Sun X, Fernando RL, Garrick DJ, Dekkers JCM: An iterative approach for efficient calculation of breed- ing values and genome-wide association analysis using weighted genomic BLUP. J Anim Sci 2011, 89(E–Suppl 2):e11.
• [35]Falconer D, Mackay T: Introduction to quantitative genetics. New York: Longman; 1996.
• [36]Choi I, Steibel JP, Bates RO, Raney NE, Rumph JM, Ernst CW: Application of alternative models to identify QTL for growth traits in an F2 Duroc x Pietrain pig resource population. BMC Genet 2010, 11:97.
• [37]Fan B, Onteru SK, Du Z-Q, Garrick DJ, Stalder KJ, Rothschild MF: Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs. PLoS One 2011, 6:e14726.
• [38]Switonski M, Stachowiak M, Cieslak J, Bartz M, Grzes M: Genetics of fat tissue accumulation in pigs: a comparative approach. J Appl Genet 2010, 51:153-168.