Volume 6 Supplement 2
Proceedings of the 15th European workshop on QTL mapping and marker assisted selection (QTLMAS)
Genomewide association analyses of the 15^{th} QTLMAS workshop data using mixed model based single locus regression analysis
 WeiXuan Fu†^{1},
 ChongLong Wang†^{1, 2},
 XiangDong Ding^{1},
 Zhe Zhang^{1},
 PeiPei Ma^{1},
 ZiQing Weng^{1},
 JianFeng Liu^{1}Email author and
 Qin Zhang^{1}Email author
DOI: 10.1186/175365616S2S5
© Fu et al.; licensee BioMed Central Ltd. 2012
Published: 21 May 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
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}{\sigma}_{a}^{2}\right)$ (where A is the additive genetic relationship matrix and ${\sigma}_{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}{\sigma}_{e}^{2}\right)$ (where I is a unit matrix and ${\sigma}_{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
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}{\sigma}_{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}{\sigma}_{e}^{2}\right)$.
For each SNP, the estimate of b and the corresponding sampling variances $Var\left(\stackrel{\wedge}{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
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.
Notes
List of abbreviation used
 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.
Declarations
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.
Authors’ Affiliations
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.View ArticlePubMed
 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.PubMed CentralView ArticlePubMed
 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.View ArticlePubMed
 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.View ArticlePubMed
 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.View ArticlePubMed
 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.PubMed CentralView ArticlePubMed
Copyright
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.