KING Tutorial: Gene Mapping
KING is a toolset to explore genotype data from a genome-wide association study (GWAS) or a sequencing project.
KING can be used to map genes for diseases and complex traits.
GENERAL INPUT FILES
The input files need to be in PLINK binary format.
Besides the standard PLINK binary format (ex.fam, ex.bim, ex.bed), two other files can be specified,
including ex.phe for phenotypes and ex.cov for covariates.
KING searches for all 5 files automatically even though only one file (ex.bed) needs to be specified in the command line.
GENOME-WIDE ASSOCIATION SCAN
Two association analyses are available in KING at the moment: TDT and multi-trait score test (mtscore). Examples are:
prompt> king -b ex.bed --tdt
prompt> king -b ex.bed --cov, --mtscore --maxP 5E-8 --invnorm
--tdt implements the well-known Transmission/Disequilibrium Test for family data that consist of parent-affected child trios.
--lmm implements a linear mixed models for association between a SNP and a quantitative trait.
The LMM implementation here is quite efficient, especially for a lot of traits, e.g., in eQTL/pQTL/meQTL/mQTL analysis,
where association needs to be examined exhaustively between each of the 10,000s of traits and 100,0000s of SNPs.
The example above applies an inverse normal transformation to all traits prior to the GWAS scan.
OTHER PARAMETERS
The following parameters can also be specified:
--prefix specifies the name of the file that stores GWAS scan results. "king" is used as default.
--cpus specifies the number of CPU cores to be used in the parallel computing. If not specified, the default number is half of the total number of (logical) cores.
--invnorm carries out inverse normal transformation for quantitative traits prior to association analysis.
--trait specifies the trait names to be analyzed in the association analysis.
--covariate specifies the covariate names to be adjusted in the association analysis.
REFERENCE
Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, Chen WM
(2010) Robust relationship inference in genome-wide association studies.
Bioinformatics 26(22):2867-2873
[Abstract]
[PDF][Citations]
======================================
Last updated: October 10, 2022 by Wei-Min Chen
|