290 likes | 646 Views
Admixture Mapping. Qunyuan Zhang Division of Statistical Genomics GEMS Course M21-621 Computational Statistical Genetics March 25, 2010. Three Mapping Strategies. Linkage Analysis (linkage): genotype & phenotype data from family (or families)
E N D
Admixture Mapping Qunyuan Zhang Division of Statistical Genomics GEMS Course M21-621 Computational Statistical Genetics March 25, 2010
Three Mapping Strategies Linkage Analysis (linkage): genotype & phenotype data from family (or families) Association Scan (LD): genotype & phenotype data from population(s) or families Admixture Mapping (LD): genotype data from admixed and ancestral populations, phenotype data from admixed populations (1) Ancestry-phenotype association mapping (2) Ancestry info for population structure control
Ancestral Population 2 Caucasians Ancestral Population 1 Africans Admixture Information (Ancestry Analysis) Admixture Mapping Admixed Population African Americans Genetic Admixture
Rationale of Admixture Mapping If a disease has some genetic factors, and the disease gene frequency in pop 2 is higher than in pop 1. After the admixture of pop 1 and 2, the diseased individuals in admixed generations will carry disease genes/alleles that have more ancestry from pop 2 than from pop 1. If a marker is linked with disease genes, because of linkage disequilibrium, the diseased individuals will also carry the marker copies that have more ancestry from pop 2 than from pop 1. Inversely, if we find a marker/locus whose ancestry from pop 2 in diseased group is significantly different from that in non-diseased group, we consider this marker/locus to be linked with (or a part of ) disease gene.
Advantages of Admixture Mapping Admixed population has more genetic variation and polymorphism than relatively pure ancestral populations. Admixture produces new LD in admixed population. Compared with ancestral populations, shorter genetic history of admixture population keeps more LD (long genetic history will destroy LD), In admixed population, LD could be detected for relatively loose linkage. Ancestry information can be used to control population stratification caused by genetic admixture. According to simulation, admixture mapping demonstrates higher power than regular methods, needs less sample size. Flexible design: case-control or case-only, qualitative or quantitative traits, no need of pedigree information
Ancestry Proportion of genetic materials descending from each founding population Population level : population admixture proportion Individual level: individual admixture proportion Individual-locus level: locus-specific ancestry
Two Ways of Using Ancestral Info. Individual Ancestry (IA) can be used as a genetic background covariate for population structure control Phenotype= a + b * Genotype + c * IA + Error Locus-specific Ancestry (LSA) can be directly used to detect association (admixture mapping) Phenotype=a + b * LSA
Individual Ancestry (IA) Estimation using MLE • G: Observed genotypes of admixed and ancestral populations • Q: Allelic frequencies in ancestral populations • P : Individual Ancestry to be estimated • Goal: obtain P that maximizes Pr(G|P,Q) • Assign prior values for Q (randomly or estimated from ancestral population genotype data) & P (randomly) • Compute P(i) by solving • Compute Q(i) by solving • Iterate Steps 1 and 2 until convergence. • Tang et al. Genetic Epidemiology, 2005(28): 289–301
ObservedG:genotypes of admixed and ancestral populations • UnknownZ :admixed individuals’ locus specific ancestries from ancestral populations • Problem: How to estimateZ ? • Maximum Likelihood Estimate(MLE): • How to obtain aZthat maximizes Pr(G|Z) ? • Zis a huge space of parameters, in which search is difficult for likelihood method. • Bayesian and Markov Chain Monte Carlo (MCMC) methods • Assume ancestral population numberK • Define prior distribution Pr(Z) underK • Use MCMC to sample from posterior distribution Pr(Z|G) = Pr(Z)∙ Pr(G|Z) • Average over large number of MCMC samples to obtain estimate ofZ • Falush et al. Genetics, 2003(164):1567–1587 Locus-specific Ancestry Estimationusing MCMC
Software STRUCTUREFalush D, Stephens M, Pritchard JK (2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164:1567–1587. ADMIXMAPHoggart CJ, Parra EJ, Shriver MD, Bonilla C, Kittles RA, Clayton DG, McKeigue PM (2003) Control of confounding of genetic associations in stratified populations. Am J Hum Genet 72:1492–1504. ANCESTRYMAPPatterson N, Hattangadi N, Lane B, Lohmueller KE, Hafler DA, Oksenberg JR, Hauser SL, Smith MW, O’Brien SJ, Altshuler D, Daly MJ, Reich D (2004) Methods for high-density admixture mapping of disease genes. Am J Hum Genet 74:979–1000
References D.C.Rife. Populations of hybrid origin as source material for the detection of linkage. Am.J.Hum.Genet. 1954, (6):26-33 R.Chakraborty et al. Adimixture as a tool for finding linked genes and detecting that difference from allelic association between loci. Proc.Natl.Acad.Sci. 1988,Vol.85:9119-9123 N. Risch. Mapping genes for complex disease using association studies with recently admixed populations. Am.J.Hum.Genet.Suppl. 1992, 51:13 … P.M.McKeigue. Prospects for admixture mapping of complex traits. Am.J.Hum.Genet. 2005, Vol.76:1-7 X.Zhu et al. Admixture mapping for hypertention loci with genome-scan markers. Nature Genetics. 2005,Vol.37(2): 177-181 Q Zhang et al. Genome-wide admixture mapping for coronary artery calcification in African Americans: the NHLBI Family Heart Study. Genet Epidemiol. 2008 Apr;32(3):264-72.
Mean=0.22 Std Dev=0.1003 (MIC) Marker Information Content (MIC ) DistributionUsed for Simulation (300 Loci) Freqency of allele k at locus i in Caucasians Freqency of allelek at locusi in Africans Allele number of locus i
400 microsatellite markers Average distance 10 cM Coronary and aortic artery calcium (CAC) African Americans 622 Subjects from 211 families Admixture Mapping CAC Loci calcified plaque Quantified by CT
Data Samples 1672 subjects from 3 populations: 622 African Americans (211 families) from FHS-SCAN 893 Caucasians (320 families) from FHS-SCAN 157 Africans (unrelated) from Marshfield Center Genotypes 302 microsatellite Loci of all subjects Average marker distance 11.9cM Phenotype Coronary and aortic artery calcium (CAC) of 622 African Americans, BLOM transformation
Statisticl Procedure Step 1 Randomly draw one subject from each family to create a sample of 688 unrelated subjects which comprises : 211 African Americans from 211 families (FHS-SCAN) 320 whites from 320 families (FHS-SCAN) 157 unrelated Africans (Marshfield Center) Step 2 Ancestry estimation, STRUCTURE 2.1 Step 3 Ancestry-CAC association analysis, regress 211 African Americans’ CAC scores on their locus-specific ancestries from Africans. Step 4 Repeat step1~step3 (100 times), obtain the average p-value of each locus Step 5 For each locus: permutation test on average p-value Number of random permutations: 10000
RESULTSAncestry Analysis at Individual Level Individual Ancestry Distribution of 622 African Americans Ancestry-from-Africans: average 77.96% (3.1%~96.9%)
RESULTSAncestry Analysis at Individual-locus Level Distribution of Locus-specific Ancestries from Africans An Example African American Ancestry from Africans 302 Microsatellite Loci ordered by chromosome and position from Chrom. 1 (4.22cM) to Chrom. 23 (104.83cM)
-log(p value) of Markers on Chromosome 4 GATA42H02
-log(p value) of Markers on Chromosome 6 GATA64D02