A protocol to adopt the mixture model by Zhu et al. for the analysis of population stratification on the data with missing allele calls

doi:10.1038/nprot.2008.129

Method Article

A protocol to adopt the mixture model by Zhu et al. for the analysis of population stratification on the data with missing allele calls

https://doi.org/10.1038/nprot.2008.129

This protocol has been posted on Protocol Exchange, an open repository of community-contributed protocols sponsored by Nature Portfolio. These protocols are posted directly on the Protocol Exchange by authors and are made freely available to the scientific community for use and comment.

Version 1

posted

You are reading this latest protocol version

Genetics

population stratification

genome wide association study

principal component analysis

expectation-maximization algorithm

There are two kinds of applications of principal component analysis (PCA) to analyze population substructures of genetic polymorphism data. One application is for an individual covariance matrix, and the other application is for a marker covariance matrix. The former method is already implemented in EIGENSTRAT [1]; the latter method, however, is not common because it cannot be applied, if data include missing typing data (allele call). Here, we describe some modification of a Mixture Model [2] (MM), so that it can handle data with missing allele calls (we call it a compensated mixture model (CMM) protocol). MM applies PCA to a marker covariance matrix before applying the normal-distribution mixture model.

Genotype data file on markers (e.g. SNPs in our GWAS on gastric cancer), which were selected so that the marker loci would be independent each other (an example of such selection criteria is given below for the analysis shown in Table 1 and Figure 1).
CMM program module (please contact us if you want to use our in-house software which is written by C++)

The calculation procedures for CMM are as follows:

Calculate allele frequencies for each locus.
Sample genotype randomly based on the allele frequencies at the missing-data loci for each of the subjects showing missing allele calls of the loci.
Calculate M times M marker covariance matrix (M is the number of marker loci).
Calculate eigenvectors up to the 3rd or 4th largest eigenvalues of the covariance matrix.
Calculate Bayesian information criterions (BICs) of the principle components, assuming K normal-distributions mixture models (K corresponds to the number of subpopulations).
Count the inferred subpopulation number K based on minimum BIC.
Iterate the above steps from 2 to 6 (we iterated this procedure 200 times in our paper).

The result on the 5,197 SNP typing data on the Chinese and Japanese population of the HapMap project (SNPs were selected by the following criteria: physical distances among the SNPs are more than 500kbp, minor allele frequency more than 3%, and missing genotype call rate less than 5%) are shown in Table 1 and Figure 1.

[1] Price, A.L., Patterson, N.J., Plenge, R.M., Weinblatt, M.E., Shadick, N.A. & Reich, D. Principal components analysis corrects for stratification in genome-wide association studies, Nature Genetics 38, 904-909 (2006).

[2] Zhu, X., Zhang, S., Zhao, H. & Cooper, R.S. Association mapping, using a mixture model for complex traits. Genet. Epidemiol. 23, 181-196 (2002).

This work was supported in Japan by the program for promotion of Fundamental Studies in Health Sciences of the National Institute of Biomedical Innovation (NiBio).

Download PDF

Version 1

posted

You are reading this latest protocol version

A protocol to adopt the mixture model by Zhu et al. for the analysis of population stratification on the data with missing allele calls

Status:

Version 1

Figures

Introduction

Equipment

Procedure

References

Acknowledgements

Associated Publications

Status:

Version 1

Privacy Policy

Terms of Service

A protocol to adopt the mixture model by Zhu et al. for the analysis of population stratification on the data with missing allele calls

Status:

Version 1

Figures

Introduction

Equipment

Procedure

References

Acknowledgements

Associated Publications

Status:

Version 1

Privacy Policy

Terms of Service

Manage Cookie Preferences