 Proceedings
 Open access
 Published:
Genomewide association analyses of the 15^{th} QTLMAS workshop data using mixed model based single locus regression analysis
BMC Proceedings volumeÂ 6, ArticleÂ number:Â S5 (2012)
Abstract
Background
The mixed model based single locus regression analysis (MMRA) method was used to analyse the common simulated dataset of the 15th QTLMAS workshop to detect potential significant association between single nucleotide polymorphisms (SNPs) and the simulated trait. A Wald chisquared statistic with df =1 was employed as test statistic and the permutation test was performed. For adjusting multiple testing, phenotypic observations were permutated 10,000 times against the genotype and pedigree data to obtain the threshold for declaring genomewide significant SNPs. Linkage disequilibrium (LD) in term of D' between significant SNPs was quantified and LD blocks were defined to indicate quantitative trait loci (QTL) regions.
Results
The estimated heritability of the simulated trait is approximately 0.30. 82 genomewide significant SNPs (P < 0.05) on chromosomes 1, 2 and 3 were detected. Through the LD blocks of the significant SNPs, we confirmed 5 and 1 QTL regions on chromosomes 1 and 3, respectively. No block was detected on chromosome 2, and no significant SNP was detected on chromosomes 4 and 5.
Conclusion
MMRA is a suitable method for detecting additive QTL and a fast method with feasibility of performing permutation test. Using LD blocks can effectively detect QTL regions.
Background
Recently, the highdensity single nucleotide polymorphism (SNP) arrays have been developed for almost all domestic animals, which offer the prerequisite of genomewide association study (GWAS), a more powerful approach for highresolution mapping of loci controlling phenotypic traits in domestic animals [1]. In GWAS, two basic designs of resource population have been widely used for association analysis, one is the casecontrol design with unrelated individuals, and the other is the familybased design with pedigree structure. Corresponding to these two designs, different approaches for association analysis have been proposed. However, there is no clear evidence showing general superiority of one approach over others. In farm animals, family based design is more relevant because of complex pedigree structure in almost all animal populations. In our previous GWAS study [2], we employed a mixed model based single locus regression analysis (MMRA) to test the association between SNPs and milk production traits in dairy cattle. We found this method was more powerful than the TDTbased single locus regression analysis. To further verify its performance in terms of power and type I error, we applied it to the common dataset provided in the 15th QTLMAS workshop.
Methods
The simulated population consisted of 3,220 individuals in two generations. The first generation consisted of 20 sires and 200 dams, which were assumed to be unrelated. Each sire mated with 10 dams and each dam produced 15 progenies, leading to a total of 3,000 individuals in the second generation. Of the 15 progenies of each dam, 10 were phenotyped for a continuous trait. All of the 3,220 individuals were genotyped for 9,990 SNP markers distributed on 5 chromosomes without missing. Each chromosome had a size of 1 Morgan (M) and carried 1,998 evenly distributed SNPs.
Variance component estimation
We applied the software DMU (Version 6, release 5.0) [3] to estimate the variance components of the simulated trait, which would be used in the subsequent association analysis, based on the following model
Where y is the vector of phenotypes of the 2,000 phenotyped individuals, Î¼ is the overall mean, a is the vector of the residual polygenic effect with \mathbf{a}~N\left(0,\mathbf{A}{\mathrm{\xcf\u0192}}_{a}^{2}\right) (where A is the additive genetic relationship matrix and {\mathrm{\xcf\u0192}}_{a}^{2} is the additive genetic variance), Z is the incidence matrix of a, and e is the vector of residual errors with \mathbf{e}~N\left(0,\mathbf{I}{\mathrm{\xcf\u0192}}_{e}^{2}\right) (where I is a unit matrix and {\mathrm{\xcf\u0192}}_{e}^{2} is the residual error variance).
Genotype quality control
We removed the 1,000 progenies without phenotypes off the genotype data, and we calculated the minor allele frequency (MAF) for each SNP for the remained 2,220 individuals (2,000 progenies and 220 parents). We found that 2,879 SNPs were homozygous (MAF = 0) for all the tested individuals and additionally 715 SNPs had a MAF less than 0.03. These SNPs were removed and 6,396 SNPs remained for the subsequent analyses.
Association analysis
The mixed model based single locus analysis [2, 4] was performed based on the following linear mixed model:
where y is the vector of phenotypes of the 2000 phenotyped individuals, Î¼ is the overall mean, Ã— is the vector of the SNP genotype indicators which takes values 0, 1 or 2 corresponding to the three genotypes 11, 12 and 22 (assuming 2 is the allele with a minor frequency), b is the regression coefficient of phenotypes on SNP genotypes (i.e., the substitution effect of the SNP), a is the vector of the residual polygenic effect with \mathbf{a}~N\left(0,\mathbf{A}{\mathrm{\xcf\u0192}}_{a}^{2}\right), Z is the incidence matrix of a, and e is the vector of residual errors with \mathbf{e}~N\left(0,\mathbf{I}{\mathrm{\xcf\u0192}}_{e}^{2}\right).
For each SNP, the estimate of b and the corresponding sampling variances Var\left(\stackrel{\xe2\u02c6\S}{b}\right) can be obtained via mixed model equations (MME), and a Wald chisquared statistic {\hat{b}}^{2}/Var\left(\hat{b}\right) with df =1 was constructed to examine whether the SNP is associated with the trait.
Statistical inference
For the analyses above, the permutation method was adopted to adjust for multiple testing from the number of SNP loci detected. In our method, the phenotypes were permuted 10,000 times against the genotype and pedigree data and the empirical distribution of the Wald chisquared statistic under the null hypothesis (no association existed between any SNP and the trait in genomewide level) was obtained using the largest Wald chisquared statistic value across all SNPs from each permuted dataset. The threshold value for declaring a significant association was determined by choosing the 95th percentile of the empirical distribution, i.e., we declared a significant SNP at a 0.05 genomewide significance level if its raw value of the Wald chisquared statistic was larger than the empirical threshold value.
For the significant SNPs, linkage disequilibrium (LD) in term of D' between them was quantified using Haploview [5] and the LD blocks were defined by the criteria of Gabriel et al. [6] with default parameters.
Results and discussion
Association analysis results
The estimates of {\mathrm{\xcf\u0192}}_{a}^{2} and {\mathrm{\xcf\u0192}}_{e}^{2} are 24.82 and 58.65, respectively, so that the heritability estimate is 0.30 approximately. The profile of the raw p values (from the chidistribution and in terms of log10 p) of all tested SNPs is shown in Figure 1. By using simply Bonferroni correction, we detected 119 significant SNPs of 0.05 genomewide significance level (raw p values < 7.82E6). However, by using permutation test, we detected 82 significant SNPs of 0.05 genomewide significance level (raw p values < 2.31E7) for the simulated trait. The 82 significant SNPs are distributed on chromosomes 1, 2 and 3, i.e., 63 on chromosome 1, 3 on chromosome 2, and 16 on chromosome 3. The significant SNPs on chromosome 1 cover a large interval between 0.15cM and 15.30cM. The simulated SNP with the biggest effect is at 2.90cM (No.58), while the most significant SNP identified is at 3.55cM (No.71). The positions of the 3 significant SNPs on chromosome 2 are 81.90cM, 83.10cM and 95.80cM, respectively. Moreover, 13 of the 16 significant SNPs on chromosome 3 cover an interval between 4.25cM and 5.65cM, and the other 3 are at 3.70cM, 16.10cM and 26.75cM, respectively. No SNPs on chromosome 4 and 5 were found to be associated with the trait significantly. If we set the significant level at 0.01 for the permutation test, the number of significant SNPs reduces to 32, of which 3 are on chromosome 3 and all others on chromosome 1.
To further pinpoint the relationship among the detected SNPs, we analysed the LD levels in terms of D' between the significant SNPs (Figures 2, 3 and 4) for chromosomes 13, respectively. Through the criteria of Gabriel et al. [6] with default parameters in Haploview [5], we defined 5 LD blocks on chromosome 1, which harbour 4 to 10 significant SNPs, and 1 LD block on chromosome 3, which harbour 10 significant SNPs with 6 outside. No block was detected on chromosome 2. The LD patterns show that these significant SNPs links to each other in different LD levels.
Comparison of the significant SNPs with the simulated QTN
On chromosome 1, there is one simulated QTN located at 2.85cM (No.57), which had the largest effect among all simulated QTNs. We detected 63 significant SNPs on this chromosome. However, the true QTN at 2.85cM has a MAF of 0 and was discarded after quality control, and the adjacent SNP at 2.90cM (No.58), which has the largest estimated effect among all significant SNPs, is accordingly considered as the putative QTN. Although a large number of pseudo significant SNPs were identified on this chromosome, the LD levels between the most effective SNP and other 62 significant SNPs (Figure 2) showed that 47 of them are in strong LD (D'>0.5) with it. This suggests that the simulated QTN may be surrogated by a suite of "ghost" QTNs nearby due to high LD level.
On chromosome 2, there are two simulated QTNs in coupling linkage phase located at 81.90cM (No.3638) and 93.75cM (No.3875), respectively. We detected 3 significant SNPs, the first is exactly at 81.90cM, and the second (No.3662, at 83.10cM) is in strong LD (D' = 0.97) with the first one (Figure 3). But the third one No.3916 is at 95.80cM and is 2.05cM away from the second simulated QTN, while the LD level between them is strong (D'=0.69, Figure 3).
On chromosome 3, there are two simulated QTN in repulsion linkage phase located at 5cM (No.4100) and 15cM (No.4300), respectively. However, the first simulated SNP on this chromosome also has a MAF of 0 and was discarded after quality control. Of the 16 significant SNPs detected, 10 are harboured in the LD block covering the interval between 4.75cM and 5.65cM with an average LD level of 0.97 (D'), in which SNP No.4101 is just adjacent to the first true QTN. The second simulated QTN is 1.10cM away from the significant SNP (No.4322) and the LD level between them is strong (D'=0.93).
The one simulated imprinting QTN on chromosome 4, and 2 simulated epistatic QTNs on chromosome 5 were not detected by our analysis. This is because our method dose not account for both imprinting effect and epistatic effect. Our method needs to be further improved to account for interaction effects between SNPs and imprinting effects from parents.
Comparison of the significant SNPs with the those with high effects estimated via Bayesian approaches
To further validate significant SNPs identified herein, we compared the most promising SNPs detected with those with highest effects estimated via Bayesian approaches (BayesA, BayesB and BayesCÏ€) reported in our another analysis on prediction of genomic breeding values for the same data set [7]. Since the results from the three Bayesian approaches are similar and BayesCÏ€ performed best, we only compare with BayesCÏ€ here. Specifically, on chromosome 1, the most effective SNP (No.58) identified by MMRA is exactly the same as that by BayesCÏ€. On chromosome 2, BayesCÏ€ revealed SNP No. 3660 with the largest effect and SNP No.3873 with the second largest effect, which are close to and in strong LD with the significant SNPs No. 3662 and No. 3916 detected by MMRA. On chromosome 3, the two promising SNPs detected by MMRA are No.4101 and No.4322, which are close to the SNP with the largest (No.4092) and the second largest (No.4331) effect estimated by BayesCÏ€, respectively. In all, most of findings herein are largely consistent with those with highest effect estimates via BayesCÏ€. This further demonstrates that the Bayesian approaches (particularly BayesCÏ€) could also sever as tools for QTL mapping, as suggested by Fan et al. [8].
Computing time
All analyses were implemented through Fortran programs and performed on an octalcore Linux Server (Intel Xeon E5504 2.00GHz; 48.00GB RAM). The time needed was about 1.5 minutes for one permutation analysis. The 10,000 permutations were performed through 8 threads, each was assigned 1,250 permutations. So, the total computing time was about 31 hours. This shows that MMRA is a fast method with feasibility of performing a large number of permutations.
Conclusion
Our results herein show that the MMRA method is suitable for detecting additive QTL, and it is a fast method with feasibility of performing permutation test. And the LD region on chromosome 3 can effectively integrate significant SNPs for QTL region detection. However, we detects only one true additive QTN (No.3638), two SNPs (No.58 and No.4101) close to two true additive QTNs (No.57 and No.4100) with many false positives, which remains to be further investigated and the MMRA method needs to be further improved to account for other nonadditive effects.
Abbreviations
 GWAS:

GenomeWide Association Study
 MMRA:

Mixed Model based Single Locus Regression Analysis
 SNP:

Single Nucleotide Polymorphisms
 QTL:

Quantitative Trait Locus
 QTN:

Quantitative Trait Nucleotide
 LD:

Linkage disequilibrium
 M:

Morgan
 MAF:

Minor Allele Frequency
 TDT:

Transmission Disequilibrium Test.
References
Andersson L: Genomewide association analysis in domestic animals: a powerful approach for genetic dissection of trait loci. Genetica. 2009, 136 (2): 341349. 10.1007/s1070900893124.
Jiang L, Liu J, Sun D, Ma P, Ding X, Yu Y, Zhang Q: Genome wide association studies for milk production traits in Chinese Holstein population. Plos One. 2010, 5 (10): e1366110.1371/journal.pone.0013661.
Madsen P, Jensen J: DMU: A user's Guide. A Package for Analysing Multivariate Mixed Models. University of Aarhus, Faculty of Agricultural Sciences, Department of Animal Breeding and Genetics. 2007
Daetwyler HD, Schenkel FS, Sargolzaei M, Robinson JA: A genome scan to detect quantitative trait loci for economically important traits in Holstein cattle using two methods and a dense single nucleotide polymorphism map. J Dairy Sci. 2008, 91 (8): 32253236. 10.3168/jds.20070333.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263265. 10.1093/bioinformatics/bth457.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, et al: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 22252229. 10.1126/science.1069424.
Wang CL, Ma PP, Zhang Z, Ding XD, Liu JF, Fu WX, Weng ZQ, Zhang Q: Comparison of five methods for genomic breeding value estimation for the common dataset of the 15^{th} QTLMAS Workshop. Proceedings of the XVth QTLMAS Workshop: 1920 May 2011; Rennes France. Edited by: Demeure O, Elsen JM, Filangi O, Le Roy P. 2012
Fan B, Onteru SK, Du ZQ, Garrick DJ, Stalder KJ, Rothschild MF: Genomewide association study identifies Loci for body composition and structural soundness traits in pigs. Plos One. 2011, 6 (2): e1472610.1371/journal.pone.0014726.
Acknowledgements
This work was supported by the National High Technology Research and Development Program of China (863 Program 2011AA100302), the National Natural Science Foundation of China (Grant No. 30972092, 30800776), and the Modern Pig Industry Technology System Program of Anhui Province.
This article has been published as part of BMC Proceedings Volume 6 Supplement 2, 2012: Proceedings of the 15th European workshop on QTL mapping and marker assisted selection (QTLMAS). The full contents of the supplement are available online at http://www.biomedcentral.com/bmcproc/supplements/6/S2.
Author information
Authors and Affiliations
Corresponding authors
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
WXF, CLW, XDD, ZZ, PPM and ZQW carried out the analyses and contributed the manuscript. JFL and QZ coordinated the analyses and drafted the writeup. All authors have read and contributed to the final text of the manuscript.
WeiXuan Fu, ChongLong Wang contributed equally to this work.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Fu, WX., Wang, CL., Ding, XD. et al. Genomewide association analyses of the 15^{th} QTLMAS workshop data using mixed model based single locus regression analysis. BMC Proc 6 (Suppl 2), S5 (2012). https://doi.org/10.1186/175365616S2S5
Published:
DOI: https://doi.org/10.1186/175365616S2S5