Submit manuscript...
eISSN: 2378-315X

Biometrics & Biostatistics International Journal

Mini Review Volume 2 Issue 4

A mini review for aggregating analyses for genome-wide association studies (GWAS)

Caiyan Li

Takeda Pharmaceuticals, USA

Correspondence: Caiyan Li, Takeda Pharmaceuticals, USA

Received: April 08, 2015 | Published: April 29, 2015

Citation: Caiyan L. A mini review for aggregating analyses for genome-wide association studies (GWAS). Biom Biostat Int J. 2015;11(2):115-116. DOI: 10.15406/bbij.2015.02.00038

Download PDF

Univariate analysis

Genome-wide Association Studies (GWAS) has emerged as a popular tool for identifying common genetic variants for complex disease. A typical case-control GWAS often involves genotyping hundreds of thousands of Single-nucleotide polymorphism (SNPs), in thousands of disease cases and healthy controls, with the goal of identifying individual SNPs that are associated with the disease (or outcome). GWAS usually consists of a discovery phase, in which SNPs are identified followed by a validation phase, in which the SNPs identified in discovery phase are replicated in a separate study cohort. The standard GWAS method is univariate analysis, in which the phenotype is regressed onto each individual SNP to generate p-value for association measurement. The SNPs are then ranked based on their univariate p-values, and a threshold is set such that any SNP with a p-value below the threshold is identified as significant. Multiplicity adjustment methods, such as bonferroni, false discovery rate (FDR) can be used to set the threshold.

Though the use of univariate analysis has helped to identify many disease-susceptibility variants, and it is easy to comprehend and conduct, there are limitations which pose difficulty in some settings. Due to the large number of tested individual SNPs (e.g., 4 million SNPs), the adjusted threshold for genome-wide significance can be extremely low, e.g. α= 10 7 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbakabeg7aHj abg2da9iaaigdacaaIWaWaaWbaaeqajuaibaGaeyOeI0IaaG4naaaa aaa@3C91@ ,1 which is very difficult to attain. Secondly, univariate-SNP analysis often has poor reproducibility. Many highly-ranked SNPs identified in the discovery phase are false positives due to the low power for detecting SNPs with small effects. Moreover, rare variants (defined as alleles with a frequency of 1%-5%) require a significant large sample size in order to be detected via standard univariate analysis. Since multivariate analysis has the ability to take into account of the correlation between multiple SNPs, as well as help reduce the multiplicity burden (e.g., lower the significance level), it could be advantageous to consider the joint effect of multiple SNPs together. To this end, researchers have established various methodologies to aggregate individual SNPs together as a set and then conduct association analyses, which has aided with the discovery of disease-related variants, especially on rare variants. Among these methods, there are two major categories: burden tests and kernel-based variance component tests. Below, we will review some popular methods within these two categories. 

Burden tests

Burden tests assess the cumulative effects of multiple variants within a genomic region by collapsing or summarizing the variants within a region by a single value, which is then used to test for association with the outcome of interest. Some popular burden tests are briefed below:

Cohort allelic sum test (CAST):2 collapse genotypes across all variants within a specific region into a single dichotomous value for each subject. A subject is coded as 1 if a rare variant is present at any of the variant sites within the region; otherwise as 0. One variation of CAST involves collapsing by counting the number of rare variants within a specific region, instead of dichotomizing.

Combined Multivariate and Collapsing (CMC):3 The CMC method is a unified approach that combines collapsing and multivariate tests. SNPs within the region are divided into subgroups based on predefined criteria, such as allele frequencies, and within each subgroup CART was used to collapse genotypes. A multivariate test (e.g., Hotellings’stest) is then applied to the subgroups for association analysis.

Weighted sum statistic (WSS):4 The WSS method specifically considers the case-control setting and collapse SNPs within a set into a single weighted average of the number of rare alleles for each individual, and then apply Wilcoxon rank sum test to compare between case and control. Specifically, the weight for SNP i is defined as w i = 1 n i q i (1 q i ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbakaadEhada WgaaqcfasaaiaadMgaaKqbagqaaiabg2da9abaaaaaaaaapeWaaSaa a8aabaWdbiaaigdaa8aabaWdbmaakaaapaqaa8qacaWGUbWdamaaBa aajuaibaWdbiaadMgaaKqba+aabeaapeGaamyCa8aadaWgaaqcfasa a8qacaWGPbaajuaGpaqabaWdbiaacIcacaaIXaGaeyOeI0IaamyCa8 aadaWgaaqcfasaa8qacaWGPbaajuaGpaqabaWdbiaacMcaaeqaaaaa aaa@4799@ where q i = m i U +1 2 n i U +2 , m i U MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbakaadghada WgaaqcfasaaiaadMgaaKqbagqaaiabg2da9abaaaaaaaaapeWaaSaa a8aabaWdbiaad2gapaWaa0baaKqbGeaapeGaamyAaaWdaeaapeGaam yvaaaajuaGcqGHRaWkcaaIXaaapaqaa8qacaaIYaGaamOBa8aadaqh aaqcfasaa8qacaWGPbaapaqaa8qacaWGvbaaaKqbakabgUcaRiaaik daaaGaaiilaiaad2gadaqhaaqcfasaaiaadMgaaeaacaWGvbaaaaaa @4A0C@ is the number of rare alleles observed for SNP i in the control group, n i U MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbakaad6gada qhaaqcfasaaiaadMgaaeaacaWGvbaaaaaa@3984@ is the number of control subjects genotypes for SNP i, and n i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbakaad6gada WgaaqcfasaaiaadMgaaKqbagqaaaaa@3937@ is the total number of subjects genotypes for SNP i. The weight is the inverse of the estimated standard deviation of the total number of rare alleles in the sample (including both cases and controls), under the null hypothesis of no frequency differences. It is used to down-weight rare alleles in constructing the weighted-sum score test.

One limitation for all these burden tests is that they implicitly assume that all rare variants influence the phenotype in the same direction and with the same magnitude of effect (after incorporating weights). However, biological mechanism leads us to believe that most variants within a region have little or no effect on phenotype; in addition some variants are protective and some are deleterious, and the magnitude of effects may vary as well. In these situations, employ the above burden tests may introduce substantial noise and largely reduce the power. Thus, researchers proposed methods that can account for the different direction and magnitude of effects for SNPs within a region.    

Sequence kernel association rest (SKAT)5

Kernel-based test methods, such as SKAT, are non-burden tests. Instead of aggregating variants, SKAT aggregates the associations between variants and the phenotype through a kernel matrix and can allow for SNP-SNP interactions, i.e., epistatic effects. SKAT is a flexible, computationally efficient, regression based approach that tests for association between variants in a region (both common and rare) and a phenotype (dichotomous or continuous) while adjusting for covariates. It borrows information from correlated SNPs that are grouped on the basis of prior biological knowledge (e.g., gene or pathway) and hence produce results with improved reproducibility and power. In this review, we assume a population-based case-control GWAS with n independent genotyped subjects. For a given SNP-set containing p SNPs, let Z i1 , Z i2 , Z i3 ,......, Z ip MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbakaadQfada WgaaqcfasaaiaadMgacaaIXaaajuaGbeaacaGGSaGaamOwamaaBaaa juaibaGaamyAaiaaikdaaKqbagqaaiaacYcacaWGAbWaaSbaaKqbGe aacaWGPbGaaG4maaqcfayabaGaaiilaiaac6cacaGGUaGaaiOlaiaa c6cacaGGUaGaaiOlaiaacYcacaWGAbWaaSbaaKqbGeaacaWGPbGaam iCaaqcfayabaaaaa@4B36@ be the genotype values for subject i( i=1,...n ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbakaadMgada qadaqaaiaadMgacqGH9aqpcaaIXaGaaiilaiaac6cacaGGUaGaaiOl aiaad6gaaiaawIcacaGLPaaaaaa@3F58@ . Let denote the case-control status for subject i (1 for case and 0 for control); X i1 , X i2 ,......, X im MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbakaadIfada WgaaqcfasaaiaadMgacaaIXaaajuaGbeaacaGGSaGaamiwamaaBaaa juaibaGaamyAaiaaikdaaKqbagqaaiaacYcacaGGUaGaaiOlaiaac6 cacaGGUaGaaiOlaiaac6cacaGGSaGaamiwamaaBaaajuaibaGaamyA aiaad2gaaKqbagqaaaaa@4716@ denote the covariates (e.g., demographic and clinical variables) that we would like to adjust for. The goal is to test the global null of whether any of the p SNPs is related to the outcome while adjusting for covariates. In evaluating the significance of the joint effect of the SNP-set, a logistic kernel-machine regression model is employed as follows:

logit P( y i =1 )= α 0 + α 1 x i1 ++ α m x im +h( z i1, z i2, , z ip ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaciiBaiaac+gacaGGNbGaamyAaiaadshacaqGGaGaamiuamaa bmaapaqaa8qacaWG5bWdamaaBaaajuaibaWdbiaadMgaaKqba+aabe aapeGaeyypa0JaaGymaaGaayjkaiaawMcaaiabg2da9iabeg7aH9aa daWgaaqcfasaa8qacaaIWaaajuaGpaqabaGaey4kaSYdbiabeg7aH9 aadaWgaaqcfasaa8qacaaIXaaajuaGpaqabaWdbiaadIhapaWaaSba aeaapeGaamyAaKqbGiaaigdaaKqba+aabeaapeGaey4kaSIaeyOjGW Raey4kaSIaeqySde2damaaBaaajuaibaWdbiaad2gaaKqba+aabeaa peGaamiEa8aadaWgaaqcfasaa8qacaWGPbGaamyBaaqcfa4daeqaa8 qacqGHRaWkcaWGObGaaiikaiaadQhapaWaaSbaaeaajuaipeGaamyA aiaaigdajuaGcaGGSaaapaqabaWdbiaadQhapaWaaSbaaeaajuaipe GaamyAaiaaikdajuaGcaGGSaaapaqabaWdbiabgAci8kaacYcacaWG 6bWdamaaBaaajuaibaWdbiaadMgacaWGWbaapaqabaqcfa4dbiaacM caaaa@6D0B@

The SNPs influence the outcome through the general function h(.), which is an arbitrary function that has a form defined only by a positive semidefinite kernel function K(. , .). For subject i, we can fully define using the kernel function K(. , .) as: h( z i1, z i2, , z ip )= i=1 n γ i ' K( Z i , Z i ' ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaamiAamaabmaapaqaa8qacaWG6bWdamaaBaaabaWdbiaadMga caaIXaGaaiilaaWdaeqaa8qacaWG6bWdamaaBaaabaWdbiaadMgaca aIYaGaaiilaaWdaeqaa8qacqGHMacVcaGGSaGaamOEa8aadaWgaaqa a8qacaWGPbGaamiCaaWdaeqaaaWdbiaawIcacaGLPaaacqGH9aqpcq GHris5daqhaaqaaiaadMgacqGH9aqpcaaIXaaabaGaamOBaaaacqaH ZoWzpaWaaSbaaeaapeGaamyAa8aadaahaaqabeaapeGaai4jaaaaa8 aabeaapeGaam4saiaacIcacaWGAbWdamaaBaaabaWdbiaadMgaa8aa beaapeGaaiilaiaadQfapaWaaSbaaeaapeGaamyAa8aadaahaaqabe aapeGaai4jaaaaa8aabeaapeGaaiykaaaa@5887@ , for some γ 1 ,, γ n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaeq4SdC2damaaBaaajqwba+Faa8qacaaIXaaajuaGpaqabaWd biaacYcacqGHMacVcaGGSaGaeq4SdC2damaaBaaajuaibaWdbiaad6 gaaKqba+aabeaaaaa@426D@ . By choosing different kernel functions, we can specify different bases and corresponding models to assess the effects of SNPs. For example, if we define K(. , .) as the linear kernel such that K( Z i , Z i ' )= j=1 p z ij z i ' j MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaam4samaabmaapaqaa8qacaWGAbWdamaaBaaajuaibaWdbiaa dMgaaKqba+aabeaapeGaaiilaiaadQfapaWaaSbaaKqbGeaapeGaam yAaKqba+aadaahaaqcfasabeaapeGaai4jaaaaaKqba+aabeaaa8qa caGLOaGaayzkaaGaeyypa0JaeyyeIu+aa0baaKqbGeaacaWGQbGaey ypa0JaaGymaaqaaiaadchaaaqcfaOaamOEa8aadaWgaaqcfasaa8qa caWGPbGaamOAaaqcfa4daeqaa8qacaWG6bWdamaaBaaajuaibaWdbi aadMgajuaGpaWaaWbaaKqbGeqabaWdbiaacEcaaaGaamOAaaqcfa4d aeqaaaaa@51ED@ then we are assuming the simple linear and logistic model defined by logit P( y i =1 )= α 0 + α 1 x i1 ++ α m x im + β 1 z i1 ++ β p z ip MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaciiBaiaac+gacaGGNbGaamyAaiaadshacaqGGaGaamiuamaa bmaapaqaa8qacaWG5bWdamaaBaaajuaibaWdbiaadMgaaKqba+aabe aapeGaeyypa0JaaGymaaGaayjkaiaawMcaaiabg2da9iabeg7aH9aa daWgaaqcfasaa8qacaaIWaaajuaGpaqabaGaey4kaSYdbiabeg7aH9 aadaWgaaqcfasaa8qacaaIXaaajuaGpaqabaWdbiaadIhapaWaaSba aKqbGeaapeGaamyAaiaaigdaaKqba+aabeaapeGaey4kaSIaeyOjGW Raey4kaSIaeqySde2damaaBaaajuaibaWdbiaad2gaaKqba+aabeaa peGaamiEa8aadaWgaaqcfasaa8qacaWGPbGaamyBaaqcfa4daeqaa8 qacqGHRaWkcqaHYoGypaWaaSbaaKqbGeaapeGaaGymaaqcfa4daeqa a8qacaWG6bWdamaaBaaajuaibaWdbiaadMgacaaIXaaajuaGpaqaba WdbiabgUcaRiabgAci8kabgUcaRiabek7aI9aadaWgaaqcfasaa8qa caWGWbaajuaGpaqabaWdbiaadQhapaWaaSbaaKqbGeaapeGaamyAai aadchaaKqba+aabeaaaaa@6DCD@ The outcome depends on the SNPs solely through the function h(.), thus, in order to test whether there is a true SPN-set effect, the following null hypothesis is considered: H 0 :h( z )=0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaamisa8aadaWgaaqcfasaa8qacaaIWaaajuaGpaqabaWdbiaa cQdacaWGObWaaeWaa8aabaWdbiaadQhaaiaawIcacaGLPaaacqGH9a qpcaaIWaaaaa@3F4E@ against the general alternative. A variance-component score test is developed to test the above hypothesis. To learn more about the technical details, please refer to references. 

The choice of kernel changes the underlying base for the nonparametric function governing the relationship between the outcome and SNPs in the set. Essentially, K(. , .) is a function that projects the genotype data from the original space to another space and then h(.) is modeled linearly in this new space. More intuitively, K( Z i , Z i ' ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaam4samaabmaapaqaa8qacaWGAbWdamaaBaaajuaibaWdbiaa dMgaa8aabeaajuaGpeGaaiilaiaadQfapaWaaSbaaKqbGeaapeGaam yAaKqba+aadaahaaqcfasabeaapeGaai4jaaaaaKqba+aabeaaa8qa caGLOaGaayzkaaaaaa@413A@ can be viewed as a function that measures the similarity between two individuals with regard to the genotypes. Some specific kernels that SKAT considers are weighted linear kernel, K( Z i , Z i ' )= j=1 p w j z ij z i ' j MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaam4samaabmaapaqaa8qacaWGAbWdamaaBaaajuaibaWdbiaa dMgaaKqba+aabeaapeGaaiilaiaadQfapaWaaSbaaKqbGeaapeGaam yAaKqba+aadaahaaqcfasabeaapeGaai4jaaaaaKqba+aabeaaa8qa caGLOaGaayzkaaGaeyypa0JaeyyeIu+aa0baaKqbGeaacaWGQbGaey ypa0JaaGymaaqaaiaadchaaaqcfaOaam4Da8aadaWgaaqcfasaa8qa caWGQbaajuaGpaqabaWdbiaadQhapaWaaSbaaKqbGeaapeGaamyAai aadQgaaKqba+aabeaapeGaamOEa8aadaWgaaqcfasaa8qacaWGPbqc fa4damaaCaaajuaibeqaa8qacaGGNaaaaiaadQgaaKqba+aabeaaaa a@54F3@ which implies that the outcome depends on SNPs in a linear fashion; weighted quadratic kernel, K( Z i , Z i ' )= (1+ j=1 p w j z ij z i ' j ) 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaam4samaabmaapaqaa8qacaWGAbWdamaaBaaajuaibaWdbiaa dMgaaKqba+aabeaapeGaaiilaiaadQfapaWaaSbaaKqbGeaapeGaam yAaKqba+aadaahaaqcfasabeaapeGaai4jaaaaaKqba+aabeaaa8qa caGLOaGaayzkaaGaeyypa0JaaiikaiaaigdacqGHRaWkcqGHris5da qhaaqcfasaaiaadQgacqGH9aqpcaaIXaaabaGaamiCaaaajuaGcaWG 3bWdamaaBaaajuaibaWdbiaadQgaaKqba+aabeaapeGaamOEa8aada Wgaaqcfasaa8qacaWGPbGaamOAaaqcfa4daeqaa8qacaWG6bWdamaa BaaajuaibaWdbiaadMgajuaGpaWaaWbaaKqbGeqabaWdbiaacEcaaa GaamOAaaqcfa4daeqaa8qacaGGPaWdamaaCaaabeqcfasaa8qacaaI Yaaaaaaa@5924@ which assumes that the model depends on the main effects and the quadratic terms of SNPs and first order of SNP-by-SNP interactions; weighted IBS (identical-by-state) kernel K( Z i , Z i ' )= j=1 p w j IBS( z ij , z i ' j ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaam4samaabmaapaqaa8qacaWGAbWdamaaBaaajuaibaWdbiaa dMgaaKqba+aabeaapeGaaiilaiaadQfapaWaaSbaaKqbGeaapeGaam yAaKqba+aadaahaaqcfasabeaapeGaai4jaaaaaKqba+aabeaaa8qa caGLOaGaayzkaaGaeyypa0JaeyyeIu+aa0baaKqbGeaacaWGQbGaey ypa0JaaGymaaqaaiaadchaaaqcfaOaam4Da8aadaWgaaqcfasaa8qa caWGQbaajuaGpaqabaWdbiaadMeacaWGcbGaam4uaiaacIcacaWG6b WdamaaBaaajuaibaWdbiaadMgacaWGQbaajuaGpaqabaWdbiaacYca caWG6bWdamaaBaaajuaibaWdbiaadMgajuaGpaWaaWbaaKqbGeqaba WdbiaacEcaaaGaamOAaaqcfa4daeqaa8qacaGGPaaaaa@5979@ which defines similarity between two individuals as the number of alleles that share IBS. The weighted IBS kernel allows for epistatic effects because it does not assume linearity or interactions of any particular order (e.g., second order). Experiences suggest using the IBS kernel when the number of interacting SNPs within the region is modest. In all these kernels,  is an allele specific weight that controls the relative importance of the  SNP. Without prior information, SKAT suggests to use w j =Beta(MA F j ;1,25) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeWaaOaaa8aabaWdbiaadEhapaWaaSbaaKqbGeaapeGaamOAaaqc fa4daeqaaaWdbeqaaiabg2da9iaadkeacaWGLbGaamiDaiaadggaca GGOaGaamytaiaadgeacaWGgbWdamaaBaaajuaibaWdbiaadQgaaKqb a+aabeaapeGaai4oaiaaigdacaGGSaGaaGOmaiaaiwdacaGGPaaaaa@47D0@ , where MA F j MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaamytaiaadgeacaWGgbWdamaaBaaajuaibaWdbiaadQgaaKqb a+aabeaaaaa@3AF7@ is the minor allele frequency of the j th MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaKqbacbaaaaaaa aapeGaamOAa8aadaahaaqabKqbGeaapeGaamiDaiaadIgaaaaaaa@39DE@ SNP. If prior information is available, weight can be selected to increase (or decrease) the weight for likely functionality.

Although SKAT makes few assumptions about rare-variants effects and provides attractive power, it has some limitations in certain settings. For example, in the setting that a large proportion of the rare variants in a given region are truly causal and influence the phenotype in the same direction, SKAT can be less powerful than burden tests. Thus, SKAT-O6 is proposed, which automatically behaves like the burden test when the burden test is more powerful than SKAT, and behaves like SKAT when SKAT is more powerful than the burden test.

The SKAT research team has developed a powerful R package which implements SKAT and SKAT-O and their Small-Sample Adjustments. The team also developed functions to calculate power and sample size to help design studies to evaluate SNP-set effects. http://www.hsph.harvard.edu/~xlin/software.html

Discussion

In this mini review, we have briefly summarized three-types of methodologies for analyzing SNP data: univariate analysis, burden tests and SKAT (SKAT-O). While univariate analysis lacks of power, it is the classic method that has been used in GWAS to help identify many disease-related SNPs. SKAT has the advantage of taking into account of correlation of a set of SNPs, thus improve power. SKAT-O is a combination of SKAT and burden tests. Personal experiences suggest carrying out both univariate analysis and SKAT (SKAT-O) while dealing with genome-wide association studies. Deep-sequencing will soon generate enormous genetic information for large disease samples and the need for statistical methods to analyze these big data is increasing. More research in this field to search for unified optimal approach which incorporates known biological information is needed.

Acknowledgments

None.

Conflicts of interest

Authors declare that there are no conflicts of interests.

References

Creative Commons Attribution License

©2015 Caiyan. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.