Power Calculation for Testing If Disease is Associated with Marker in a Case-Control Study Using the GeneticsDesign Package

Size: px
Start display at page:

Download "Power Calculation for Testing If Disease is Associated with Marker in a Case-Control Study Using the GeneticsDesign Package"

Transcription

1 Power Calculation for Testing If Disease is Associated with Marker in a Case-Control Study Using the GeneticsDesign Package Weiliang Qiu weiliang.qiu@gmail.com Ross Lazarus ross.lazarus@channing.harvard.edu October 30, Introduction In genetics association studies, investigators are interested in testing the association between disease and DSL (disease susceptibility locus) in a case-control study. However, DSL is usually unknown. Hence what one can do is to test the association between disease and DSL indirectly by testing the association between disease and a marker. The power of the association between disease and DSL depends on the power of the association between disease and the marker, and on the LD (linkage disequilibrium) between disease and the marker. Assume that both DSL and the marker are biallelic and assume additive model. Given the MAFs (minor allele frequencies) for both DSL and the marker, relative risks, r 2 or D between DSL and marker, the numbers of cases and controls, and the significant level for hypothesis testing, the tools in the GeneticsDesign package can calculate the power for the association test that disease is associated with DSL in a case-control study. This information can help investigator to determine appropriate sample sizes in experiment design stage. 2 Examples To call the functions in the R package GeneticsDesign, we first need to load it into R: > library(geneticsdesign) The following is a sample code to get a table of power for different combinations of high risk allele frequency P r(a) and genotype relative risk RR(Aa aa). 1

2 > set1<-seq(from0.1, to0.5, by0.1) > set2<-c(1.25, 1.5, 1.75, 2.0) > len1<-length(set1) > len2<-length(set2) > mat<-matrix(0, nrowlen1, ncollen2) > rownames(mat)<-paste("maf", set1, sep"") > colnames(mat)<-paste("rraa", set2, sep"") > for(i in 1:len1) + { a<-set1[i] + for(j in 1:len2) + { b<-set2[j] + res<-gpc.default(paa,pd0.1,rraa(1+b)/2, RRAAb, Dprime1,pBa, quiett) + mat[i,j]<-res$power + } + } > print(round(mat,3)) RRAA1.25 RRAA1.5 RRAA1.75 RRAA2 MAF MAF MAF MAF MAF A Technical Details Suppose that A is the minor allele for the disease locus; a is the common allele for the disease locus; B is the minor allele for the marker locus; and b is the common allele for the marker locus. We use D to denote the disease status diseased and use D to denote the disease status not diseased. Given the minor allele frequencies P r(a) and P r(b) and Linkage disequilibrium (LD) measure D, we can get haplotype frequencies P r(ab), P r(ab), P r(ab), and P r(ab): where P r(ab) P r(a)p r(b) + D P r(ab) P r(a)p r(b) D P r(ab) P r(a)p r(b) D P r(ab) P r(a)p r(b) + D, D D d max 2

3 and That is, d max { min [P r(a)p r(b), P r(a)p r(b)], if D > 0, max [ P r(a)p r(b), P r(a)p r(b)], if D < 0. D P r(ab) P r(a)p r(b). Note that D > 0 means P r(ab) > P r(a)p r(b), i.e., the probability of occuring the haplotype AB is higher than the probability that the haplotype AB occurs merely by chance. D < 0 means P r(ab) < P r(a)p r(b), i.e., the probability of occuring the haplotype AB is smaller than the probability that the haplotype AB occurs merely by chance. In other words, given the same sample size, we can observe haplotype AB more often when D > 0 than when D < 0. Hence the power of association test that marker is associated with disease is larger when D > 0 than when D < 0. Hence, we assume that D > 0. D can also be rewritten as D P r(b A)P r(a) P r(a)p r(b) [P r(b A) P r(b)] P r(a). Hence D > 0 is equivalent to P r(b A) > P r(a) which indicates positive association between minor allele A in disease locus and minor allele B in marker locus. The relative risks are RR(AA aa) P r(d AA) P r(d aa) RR(Aa aa) P r(d Aa) P r(d aa) The disease prevalence can be rewritten as P r(d AA)P r(aa) + P r(d Aa)P r(aa) + P r(d aa)p r(aa). Dividing both sides by P r(d aa), we get Hence P r(d aa) P r(d AA) P r(d Aa) P r(d aa) P r(aa) + P r(aa) + P r(d aa) P r(d aa) P r(d aa) P r(aa) RR(AA aa) [P r(a)] 2 + RR(Aa aa)2p r(a)p r(a) + [P r(a)] 2. P r(d aa) RR(AA aa) [P r(a)] 2 + RR(Aa aa)2p r(a)p r(a) + [P r(a)] 2 P r(d Aa) RR(Aa aa)p r(d aa) P r(d AA) RR(AA aa)p r(d aa) 3

4 Once we obtain the penetrances of disease locus (P r(d aa), P r(d Aa), and P r(d AA)), we can get the sampling probabilities: P r(d, BB) P r(bb D) P r(d, BB, AA) + P r(d, BB, Aa) + P r(d, BB, aa) P r(d BB, AA)P r(bb, AA) + P r(d BB, Aa)P r(bb, Aa) + P r(d BB, aa)p r(bb, aa) P r(bb D) P r(d, Bb) P r(d, Bb, AA) + P r(d, Bb, Aa) + P r(d, Bb, aa) P r(d Bb, AA)P r(bb, AA) + P r(d Bb, Aa)P r(bb, Aa) + P r(d Bb, aa)p r(bb, aa) P r(d, bb) P r(bb D) P r(d, bb, AA) + P r(d, bb, Aa) + P r(d, bb, aa) P r(d bb, AA)P r(bb, AA) + P r(d bb, Aa)P r(bb, Aa) + P r(d bb, aa)p r(bb, aa) We assume that P r(d BB, AA) P r(d AA), P r(d BB, Aa) P r(d Aa), P r(d BB, aa) P r(d aa). We also have P r(bb, AA) [P r(ab)] 2 P r(bb, Aa) 2P r(ab)p r(ab) P r(bb, aa) [P r(ab)] 2 P r(bb, AA) 2P r(ab)p r(ab) P r(bb, Aa) 2 [P r(ab) P r(ab) + P r(ab)p r(ab)] P r(bb, aa) 2P r(ab)p r(ab) P r(bb, AA) [P r(ab)] 2 P r(bb, Aa) 2P r(ab)p r(ab) P r(bb, aa) [P r(ab)] 2 4

5 Hence P r(d BB, AA)P r(bb, AA) + P r(d BB, Aa)P r(bb, Aa) + P r(d BB, aa)p r(bb, aa) P r(bb D) P r(d AA) [P r(ab)]2 + P r(d Aa) [2P r(ab)p r(ab)] + P r(d aa) [P r(ab)] 2 P r(d Bb, AA)P r(bb, AA) + P r(d Bb, Aa)P r(bb, Aa) + P r(d Bb, aa)p r(bb, aa) P r(bb D) P r(d AA) [2P r(ab)p r(ab)] + P r(d Aa) {2 [P r(ab)p r(ab) + P r(ab)p r(ab)]} P r(d aa) [2P r(ab)p r(ab)] + P r(d bb, AA)P r(bb, AA) + P r(d bb, Aa)P r(bb, Aa) + P r(d bb, aa)p r(bb, aa) P r(bb D) P r(d AA) [P r(ab)]2 + P r(d Aa) [2P r(ab)p r(ab)] + P r(d aa) [P r(ab)] 2 To calculate the sampling probabilities P r(bb D), P r(bb D), and P r(bb D), we can apply Bayesian rules again: P r(bb D) P r( D BB)P r(bb) P r( D) P r(bb D) P r( D Bb)P r(bb) P r( D) P r(bb D) P r( D bb)p r(bb) P r( D) [1 P r(d BB)][P r(b)]2 1 [1 P r(d Bb)]2P r(b)p r(b) 1 [1 P r(d bb)][p r(b)]2 1 and P r(d BB) P r(d Bb) P r(d bb) P r(bb D) P r(bb D) P r(bb) [P r(b)] 2 P r(bb D) P r(bb D) P r(bb) [2P r(b)p r(b)] P r(bb D) P r(bb D) P r(bb) [P r(b)] 2. 5

6 The expected allele frequencies are The counts for cases P r(b D) 2n casep r(bb D) + n case P r(bb D) 2n case P r(bb D) + P r(bb D)/2 P r(b D) 2n casep r(bb D) + n case P r(bb D) 2n case P r(bb D) + P r(bb D)/2 P r(b D) P r(bb D) + P r(bb D)/2 P r(b D) P r(bb D) + P r(bb D)/2. n 2c n cases P r(bb D), n 1c n cases P r(bb D), n 0c n cases P r(bb D). The counts for controls n 2n n controls P r(bb D), n 1n n controls P r(bb D), n 0n n controls P r(bb D). Table 1: Counts for cases and controls marker genotype cases controls BB n 2c n 2n Bb n 1c n 1n bb n 0c n 0n total n cases n controls Odds ratios are Define OR(BB bb) n 2cn 0n n 2n n 0c P r(bb D)P r(bb D) P r(bb D)P r(bb D) OR(Bb bb) n 1cn 0n n 1n n 0c P r(bb D)P r(bb D) P r(bb D)P r(bb D) U 3 ( S x i N r i R ) N s i, i1 6

7 where We have R n 0c + n 1c + n 2c, S n 0n + n 1n + n 2n, N R + S.x 1 2, x 2 1, x 3 0, r 1 n 0c, r 2 n 1c, r 3 n 2c, s 1 n 0n, s 2 n 1n, s 3 n 2n. V V ar(u) RS N N 3 ( 3 3 ) 2 x 2 i n i x i n i. i1 i1 The Linear Trend Test Statistic 1 Z 2 T ( U V ) 2 asymptotically follows a χ 2 1 distribution under the null hypothesis. Under alternative hypothesis, Z 2 T asymptotically follows a non-central chi-square distribution χ2 1,λ T with non-centrality parameter where λ T RS [ 3 i1 x i(p 1i p 0i ) ] 2 3 i1 x2 i (Rp 0i + Sp 1i ) [ 3 i1 x i(rp 0i + Sp 1i ] 2 /N, x 1 2, x 2 1, x 3 0, n 1 n 2c + n 2n, n 2 n 1c + n 1n, n 3 n 0c + n 0n, p 01 P r(bb D), p 02 P r(bb D), p 03 P r(bb D), p 11 P r(bb D), p 12 P r(bb D), p 13 P r(bb D). Given significant level α, the power 1 β satisfies: References P r(z 2 T > c H 0 ) α P r(z 2 T > c H a ) 1 β. [1] Armitage, P. Tests for linear trends in proportions and frequencies. Biometrics, 11, , [2] Cochran, W. G. Some methods for strengthening the common chi-squared tests. Biometrics, 10, ,

8 [3] Gordon D. and Finch S. J. and Nothnagel M. and Ott J. Power and sample size calculations for case-control genetic association tests when errors are present: application to single nucleotide polymorphisms. Hum. Hered., 54:22-33, [4] Gordon D. and Haynes C. and Blumenfeld J. and Finch S. J. PAWE-3D: visualizing Power for Association With Error in case/control genetic studies of complex traits. Bioinformatics, 21: , [5] Purcell S. and Cherny S. S. and Sham P. C. Genetic Power Calculator: design of linkage and association genetic mapping studies of complex trait Bioinformatics, 19(1): , [6] Sham P. Statistics in Human Genetics. Arnold Applications of Statistics,

BST227 Introduction to Statistical Genetics. Lecture 4: Introduction to linkage and association analysis

BST227 Introduction to Statistical Genetics. Lecture 4: Introduction to linkage and association analysis BST227 Introduction to Statistical Genetics Lecture 4: Introduction to linkage and association analysis 1 Housekeeping Homework #1 due today Homework #2 posted (due Monday) Lab at 5:30PM today (FXB G13)

More information

For more information about how to cite these materials visit

For more information about how to cite these materials visit Author(s): Kerby Shedden, Ph.D., 2010 License: Unless otherwise noted, this material is made available under the terms of the Creative Commons Attribution Share Alike 3.0 License: http://creativecommons.org/licenses/by-sa/3.0/

More information

New Enhancements: GWAS Workflows with SVS

New Enhancements: GWAS Workflows with SVS New Enhancements: GWAS Workflows with SVS August 9 th, 2017 Gabe Rudy VP Product & Engineering 20 most promising Biotech Technology Providers Top 10 Analytics Solution Providers Hype Cycle for Life sciences

More information

Transmission Disequilibrium Test in GWAS

Transmission Disequilibrium Test in GWAS Department of Computer Science Brown University, Providence sorin@cs.brown.edu November 10, 2010 Outline 1 Outline 2 3 4 The transmission/disequilibrium test (TDT) was intro- duced several years ago by

More information

Statistical Tests for X Chromosome Association Study. with Simulations. Jian Wang July 10, 2012

Statistical Tests for X Chromosome Association Study. with Simulations. Jian Wang July 10, 2012 Statistical Tests for X Chromosome Association Study with Simulations Jian Wang July 10, 2012 Statistical Tests Zheng G, et al. 2007. Testing association for markers on the X chromosome. Genetic Epidemiology

More information

CS2220 Introduction to Computational Biology

CS2220 Introduction to Computational Biology CS2220 Introduction to Computational Biology WEEK 8: GENOME-WIDE ASSOCIATION STUDIES (GWAS) 1 Dr. Mengling FENG Institute for Infocomm Research Massachusetts Institute of Technology mfeng@mit.edu PLANS

More information

Non-parametric methods for linkage analysis

Non-parametric methods for linkage analysis BIOSTT516 Statistical Methods in Genetic Epidemiology utumn 005 Non-parametric methods for linkage analysis To this point, we have discussed model-based linkage analyses. These require one to specify a

More information

Statistical power and significance testing in large-scale genetic studies

Statistical power and significance testing in large-scale genetic studies STUDY DESIGNS Statistical power and significance testing in large-scale genetic studies Pak C. Sham 1 and Shaun M. Purcell 2,3 Abstract Significance testing was developed as an objective method for summarizing

More information

Human population sub-structure and genetic association studies

Human population sub-structure and genetic association studies Human population sub-structure and genetic association studies Stephanie A. Santorico, Ph.D. Department of Mathematical & Statistical Sciences Stephanie.Santorico@ucdenver.edu Global Similarity Map from

More information

Ascertainment Through Family History of Disease Often Decreases the Power of Family-based Association Studies

Ascertainment Through Family History of Disease Often Decreases the Power of Family-based Association Studies Behav Genet (2007) 37:631 636 DOI 17/s10519-007-9149-0 ORIGINAL PAPER Ascertainment Through Family History of Disease Often Decreases the Power of Family-based Association Studies Manuel A. R. Ferreira

More information

Imaging Genetics: Heritability, Linkage & Association

Imaging Genetics: Heritability, Linkage & Association Imaging Genetics: Heritability, Linkage & Association David C. Glahn, PhD Olin Neuropsychiatry Research Center & Department of Psychiatry, Yale University July 17, 2011 Memory Activation & APOE ε4 Risk

More information

Indian Journal of Nephrology Indian J Nephrol 2001;11: 88-97

Indian Journal of Nephrology Indian J Nephrol 2001;11: 88-97 88 Indian Journal of Nephrology Indian J Nephrol 2001;11: 88-97 ARTICLE HLA gene and haplotype frequency in renal transplant recipients and donors of Uttar Pradesh (North India) S Agrawal, AK Singh, RK

More information

Roadmap. Inbreeding How inbred is a population? What are the consequences of inbreeding?

Roadmap. Inbreeding How inbred is a population? What are the consequences of inbreeding? 1 Roadmap Quantitative traits What kinds of variation can selection work on? How much will a population respond to selection? Heritability How can response be restored? Inbreeding How inbred is a population?

More information

Genetics and Pharmacogenetics in Human Complex Disorders (Example of Bipolar Disorder)

Genetics and Pharmacogenetics in Human Complex Disorders (Example of Bipolar Disorder) Genetics and Pharmacogenetics in Human Complex Disorders (Example of Bipolar Disorder) September 14, 2012 Chun Xu M.D, M.Sc, Ph.D. Assistant professor Texas Tech University Health Sciences Center Paul

More information

Assessing Accuracy of Genotype Imputation in American Indians

Assessing Accuracy of Genotype Imputation in American Indians Assessing Accuracy of Genotype Imputation in American Indians Alka Malhotra*, Sayuko Kobes, Clifton Bogardus, William C. Knowler, Leslie J. Baier, Robert L. Hanson Phoenix Epidemiology and Clinical Research

More information

Nonparametric Linkage Analysis. Nonparametric Linkage Analysis

Nonparametric Linkage Analysis. Nonparametric Linkage Analysis Limitations of Parametric Linkage Analysis We previously discued parametric linkage analysis Genetic model for the disease must be specified: allele frequency parameters and penetrance parameters Lod scores

More information

IN SILICO EVALUATION OF DNA-POOLED ALLELOTYPING VERSUS INDIVIDUAL GENOTYPING FOR GENOME-WIDE ASSOCIATION STUDIES OF COMPLEX DISEASE.

IN SILICO EVALUATION OF DNA-POOLED ALLELOTYPING VERSUS INDIVIDUAL GENOTYPING FOR GENOME-WIDE ASSOCIATION STUDIES OF COMPLEX DISEASE. IN SILICO EVALUATION OF DNA-POOLED ALLELOTYPING VERSUS INDIVIDUAL GENOTYPING FOR GENOME-WIDE ASSOCIATION STUDIES OF COMPLEX DISEASE By Siddharth Pratap Thesis Submitted to the Faculty of the Graduate School

More information

Lecture 6: Linkage analysis in medical genetics

Lecture 6: Linkage analysis in medical genetics Lecture 6: Linkage analysis in medical genetics Magnus Dehli Vigeland NORBIS course, 8 th 12 th of January 2018, Oslo Approaches to genetic mapping of disease Multifactorial disease Monogenic disease Syke

More information

An Introduction to Quantitative Genetics I. Heather A Lawson Advanced Genetics Spring2018

An Introduction to Quantitative Genetics I. Heather A Lawson Advanced Genetics Spring2018 An Introduction to Quantitative Genetics I Heather A Lawson Advanced Genetics Spring2018 Outline What is Quantitative Genetics? Genotypic Values and Genetic Effects Heritability Linkage Disequilibrium

More information

Replacing IBS with IBD: The MLS Method. Biostatistics 666 Lecture 15

Replacing IBS with IBD: The MLS Method. Biostatistics 666 Lecture 15 Replacing IBS with IBD: The MLS Method Biostatistics 666 Lecture 5 Previous Lecture Analysis of Affected Relative Pairs Test for Increased Sharing at Marker Expected Amount of IBS Sharing Previous Lecture:

More information

Chapter 2. Linkage Analysis. JenniferH.BarrettandM.DawnTeare. Abstract. 1. Introduction

Chapter 2. Linkage Analysis. JenniferH.BarrettandM.DawnTeare. Abstract. 1. Introduction Chapter 2 Linkage Analysis JenniferH.BarrettandM.DawnTeare Abstract Linkage analysis is used to map genetic loci using observations on relatives. It can be applied to both major gene disorders (parametric

More information

A Comparison of Sample Size and Power in Case-Only Association Studies of Gene-Environment Interaction

A Comparison of Sample Size and Power in Case-Only Association Studies of Gene-Environment Interaction American Journal of Epidemiology ª The Author 2010. Published by Oxford University Press on behalf of the Johns Hopkins Bloomberg School of Public Health. This is an Open Access article distributed under

More information

Single SNP/Gene Analysis. Typical Results of GWAS Analysis (Single SNP Approach) Typical Results of GWAS Analysis (Single SNP Approach)

Single SNP/Gene Analysis. Typical Results of GWAS Analysis (Single SNP Approach) Typical Results of GWAS Analysis (Single SNP Approach) High-Throughput Sequencing Course Gene-Set Analysis Biostatistics and Bioinformatics Summer 28 Section Introduction What is Gene Set Analysis? Many names for gene set analysis: Pathway analysis Gene set

More information

Haplotypes of VKORC1, NQO1 and GGCX, their effect on activity levels of vitamin K-dependent coagulation factors, and the risk of venous thrombosis

Haplotypes of VKORC1, NQO1 and GGCX, their effect on activity levels of vitamin K-dependent coagulation factors, and the risk of venous thrombosis Haplotypes of VKORC1, NQO1 and GGCX, their effect on activity levels of vitamin K-dependent coagulation factors, and the risk of venous thrombosis Haplotypes of VKORC1, NQO1 and GGCX, their effect on activity

More information

Transmission Disequilibrium Methods for Family-Based Studies Daniel J. Schaid Technical Report #72 July, 2004

Transmission Disequilibrium Methods for Family-Based Studies Daniel J. Schaid Technical Report #72 July, 2004 Transmission Disequilibrium Methods for Family-Based Studies Daniel J. Schaid Technical Report #72 July, 2004 Correspondence to: Daniel J. Schaid, Ph.D., Harwick 775, Division of Biostatistics Mayo Clinic/Foundation,

More information

Association mapping (qualitative) Association scan, quantitative. Office hours Wednesday 3-4pm 304A Stanley Hall. Association scan, qualitative

Association mapping (qualitative) Association scan, quantitative. Office hours Wednesday 3-4pm 304A Stanley Hall. Association scan, qualitative Association mapping (qualitative) Office hours Wednesday 3-4pm 304A Stanley Hall Fig. 11.26 Association scan, qualitative Association scan, quantitative osteoarthritis controls χ 2 test C s G s 141 47

More information

Heritability and genetic correlations explained by common SNPs for MetS traits. Shashaank Vattikuti, Juen Guo and Carson Chow LBM/NIDDK

Heritability and genetic correlations explained by common SNPs for MetS traits. Shashaank Vattikuti, Juen Guo and Carson Chow LBM/NIDDK Heritability and genetic correlations explained by common SNPs for MetS traits Shashaank Vattikuti, Juen Guo and Carson Chow LBM/NIDDK The Genomewide Association Study. Manolio TA. N Engl J Med 2010;363:166-176.

More information

Association-heterogeneity mapping identifies an Asian-specific association of the GTF2I locus with rheumatoid arthritis

Association-heterogeneity mapping identifies an Asian-specific association of the GTF2I locus with rheumatoid arthritis Supplementary Material Association-heterogeneity mapping identifies an Asian-specific association of the GTF2I locus with rheumatoid arthritis Kwangwoo Kim 1,, So-Young Bang 1,, Katsunori Ikari 2,3, Dae

More information

Decomposition of the Genotypic Value

Decomposition of the Genotypic Value Decomposition of the Genotypic Value 1 / 17 Partitioning of Phenotypic Values We introduced the general model of Y = G + E in the first lecture, where Y is the phenotypic value, G is the genotypic value,

More information

Using ancestry estimates as tools to better understand group or individual differences in disease risk or disease outcomes

Using ancestry estimates as tools to better understand group or individual differences in disease risk or disease outcomes Using ancestry estimates as tools to better understand group or individual differences in disease risk or disease outcomes Jill Barnholtz-Sloan, Ph.D. Assistant Professor Case Comprehensive Cancer Center

More information

Regression-based Linkage Analysis Methods

Regression-based Linkage Analysis Methods Chapter 1 Regression-based Linkage Analysis Methods Tao Wang and Robert C. Elston Department of Epidemiology and Biostatistics Case Western Reserve University, Wolstein Research Building 2103 Cornell Road,

More information

Genetic association analysis incorporating intermediate phenotypes information for complex diseases

Genetic association analysis incorporating intermediate phenotypes information for complex diseases University of Iowa Iowa Research Online Theses and Dissertations Fall 2011 Genetic association analysis incorporating intermediate phenotypes information for complex diseases Yafang Li University of Iowa

More information

Performing. linkage analysis using MERLIN

Performing. linkage analysis using MERLIN Performing linkage analysis using MERLIN David Duffy Queensland Institute of Medical Research Brisbane, Australia Overview MERLIN and associated programs Error checking Parametric linkage analysis Nonparametric

More information

Quantitative Genetic Methods to Dissect Heterogeneity in Complex Traits

Quantitative Genetic Methods to Dissect Heterogeneity in Complex Traits Virginia Commonwealth University VCU Scholars Compass Theses and Dissertations Graduate School 2012 Quantitative Genetic Methods to Dissect Heterogeneity in Complex Traits T. Bernard Bigdeli Virginia Commonwealth

More information

Mendelian Inheritance. Jurg Ott Columbia and Rockefeller Universities New York

Mendelian Inheritance. Jurg Ott Columbia and Rockefeller Universities New York Mendelian Inheritance Jurg Ott Columbia and Rockefeller Universities New York Genes Mendelian Inheritance Gregor Mendel, monk in a monastery in Brünn (now Brno in Czech Republic): Breeding experiments

More information

Association of XRCC5 polymorphisms with COPD and COPD-related phenotypes in the Han Chinese population: a case-control cohort study

Association of XRCC5 polymorphisms with COPD and COPD-related phenotypes in the Han Chinese population: a case-control cohort study Association of XRCC5 polymorphisms with COPD and COPD-related phenotypes in the Han Chinese population: a case-control cohort study B. Wang 1 *, J. Yang 1 *, J. Xiao 1, B. Liang 1, H.X. Zhou 1, Z. Su 1,

More information

Introduction to the Genetics of Complex Disease

Introduction to the Genetics of Complex Disease Introduction to the Genetics of Complex Disease Jeremiah M. Scharf, MD, PhD Departments of Neurology, Psychiatry and Center for Human Genetic Research Massachusetts General Hospital Breakthroughs in Genome

More information

Introduction to linkage and family based designs to study the genetic epidemiology of complex traits. Harold Snieder

Introduction to linkage and family based designs to study the genetic epidemiology of complex traits. Harold Snieder Introduction to linkage and family based designs to study the genetic epidemiology of complex traits Harold Snieder Overview of presentation Designs: population vs. family based Mendelian vs. complex diseases/traits

More information

Supplementary Figure 1: Attenuation of association signals after conditioning for the lead SNP. a) attenuation of association signal at the 9p22.

Supplementary Figure 1: Attenuation of association signals after conditioning for the lead SNP. a) attenuation of association signal at the 9p22. Supplementary Figure 1: Attenuation of association signals after conditioning for the lead SNP. a) attenuation of association signal at the 9p22.32 PCOS locus after conditioning for the lead SNP rs10993397;

More information

Supplementary Online Content

Supplementary Online Content Supplementary Online Content Hartwig FP, Borges MC, Lessa Horta B, Bowden J, Davey Smith G. Inflammatory biomarkers and risk of schizophrenia: a 2-sample mendelian randomization study. JAMA Psychiatry.

More information

Handling Immunogenetic Data Managing and Validating HLA Data

Handling Immunogenetic Data Managing and Validating HLA Data Handling Immunogenetic Data Managing and Validating HLA Data Steven J. Mack PhD Children s Hospital Oakland Research Institute 16 th IHIW & Joint Conference Sunday 3 June, 2012 Overview 1. Master Analytical

More information

GENETIC LINKAGE ANALYSIS

GENETIC LINKAGE ANALYSIS Atlas of Genetics and Cytogenetics in Oncology and Haematology GENETIC LINKAGE ANALYSIS * I- Recombination fraction II- Definition of the "lod score" of a family III- Test for linkage IV- Estimation of

More information

Genome-wide association studies (case/control and family-based) Heather J. Cordell, Institute of Genetic Medicine Newcastle University, UK

Genome-wide association studies (case/control and family-based) Heather J. Cordell, Institute of Genetic Medicine Newcastle University, UK Genome-wide association studies (case/control and family-based) Heather J. Cordell, Institute of Genetic Medicine Newcastle University, UK GWAS For the last 8 years, genome-wide association studies (GWAS)

More information

Introduction to Genetics and Genomics

Introduction to Genetics and Genomics 2016 Introduction to enetics and enomics 3. ssociation Studies ggibson.gt@gmail.com http://www.cig.gatech.edu Outline eneral overview of association studies Sample results hree steps to WS: primary scan,

More information

Lecture 1 Mendelian Inheritance

Lecture 1 Mendelian Inheritance Genes Mendelian Inheritance Lecture 1 Mendelian Inheritance Jurg Ott Gregor Mendel, monk in a monastery in Brünn (now Brno in Czech Republic): Breeding experiments with the garden pea: Flower color and

More information

The Loss of Heterozygosity (LOH) Algorithm in Genotyping Console 2.0

The Loss of Heterozygosity (LOH) Algorithm in Genotyping Console 2.0 The Loss of Heterozygosity (LOH) Algorithm in Genotyping Console 2.0 Introduction Loss of erozygosity (LOH) represents the loss of allelic differences. The SNP markers on the SNP Array 6.0 can be used

More information

Association between interleukin-17a polymorphism and coronary artery disease susceptibility in the Chinese Han population

Association between interleukin-17a polymorphism and coronary artery disease susceptibility in the Chinese Han population Association between interleukin-17a polymorphism and coronary artery disease susceptibility in the Chinese Han population G.B. Su, X.L. Guo, X.C. Liu, Q.T. Cui and C.Y. Zhou Department of Cardiothoracic

More information

Supplementary Methods. 1. Cancer Genetic Markers of Susceptibility (CGEMS) Prostate Cancer Genome-Wide Association Scan

Supplementary Methods. 1. Cancer Genetic Markers of Susceptibility (CGEMS) Prostate Cancer Genome-Wide Association Scan Supplementary Methods 1. Cancer Genetic Markers of Susceptibility (CGEMS) Prostate Cancer Genome-Wide Association Scan The CGEMS data portal provides public access to summary results for approximately

More information

Effects of age-at-diagnosis and duration of diabetes on GADA and IA-2A positivity

Effects of age-at-diagnosis and duration of diabetes on GADA and IA-2A positivity Effects of age-at-diagnosis and duration of diabetes on GADA and IA-2A positivity Duration of diabetes was inversely correlated with age-at-diagnosis (ρ=-0.13). However, as backward stepwise regression

More information

Statistical Genetics : Gene Mappin g through Linkag e and Associatio n

Statistical Genetics : Gene Mappin g through Linkag e and Associatio n Statistical Genetics : Gene Mappin g through Linkag e and Associatio n Benjamin M Neale Manuel AR Ferreira Sarah E Medlan d Danielle Posthuma About the editors List of contributors Preface Acknowledgements

More information

MBG* Animal Breeding Methods Fall Final Exam

MBG* Animal Breeding Methods Fall Final Exam MBG*4030 - Animal Breeding Methods Fall 2007 - Final Exam 1 Problem Questions Mick Dundee used his financial resources to purchase the Now That s A Croc crocodile farm that had been operating for a number

More information

(b) What is the allele frequency of the b allele in the new merged population on the island?

(b) What is the allele frequency of the b allele in the new merged population on the island? 2005 7.03 Problem Set 6 KEY Due before 5 PM on WEDNESDAY, November 23, 2005. Turn answers in to the box outside of 68-120. PLEASE WRITE YOUR ANSWERS ON THIS PRINTOUT. 1. Two populations (Population One

More information

Research Article Power Estimation for Gene-Longevity Association Analysis Using Concordant Twins

Research Article Power Estimation for Gene-Longevity Association Analysis Using Concordant Twins Genetics Research International, Article ID 154204, 8 pages http://dx.doi.org/10.1155/2014/154204 Research Article Power Estimation for Gene-Longevity Association Analysis Using Concordant Twins Qihua

More information

Passing It On. QUESTION: How are inherited characteristics passed from parent to offspring? toothpicks - red and green

Passing It On. QUESTION: How are inherited characteristics passed from parent to offspring? toothpicks - red and green Passing It On QUESTION: How are inherited characteristics passed from parent to offspring? MATERIALS: allele sacks allele sets (7) gumdrops - green and black map pins marshmallows push pins - green and

More information

Alzheimer Disease and Complex Segregation Analysis p.1/29

Alzheimer Disease and Complex Segregation Analysis p.1/29 Alzheimer Disease and Complex Segregation Analysis Amanda Halladay Dalhousie University Alzheimer Disease and Complex Segregation Analysis p.1/29 Outline Background Information on Alzheimer Disease Alzheimer

More information

5/2/18. After this class students should be able to: Stephanie Moon, Ph.D. - GWAS. How do we distinguish Mendelian from non-mendelian traits?

5/2/18. After this class students should be able to: Stephanie Moon, Ph.D. - GWAS. How do we distinguish Mendelian from non-mendelian traits? corebio II - genetics: WED 25 April 2018. 2018 Stephanie Moon, Ph.D. - GWAS After this class students should be able to: 1. Compare and contrast methods used to discover the genetic basis of traits or

More information

Sponsored document from The Journal of Allergy and Clinical Immunology

Sponsored document from The Journal of Allergy and Clinical Immunology Sponsored document from The Journal of Allergy and Clinical Immunology Chromosome 11q13.5 variant associated with childhood eczema: An effect supplementary to filaggrin mutations Gráinne M. O'Regan a,

More information

University of Groningen. Metabolic risk in people with psychotic disorders Bruins, Jojanneke

University of Groningen. Metabolic risk in people with psychotic disorders Bruins, Jojanneke University of Groningen Metabolic risk in people with psychotic disorders Bruins, Jojanneke IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from

More information

Supplementary Figure 1. Principal components analysis of European ancestry in the African American, Native Hawaiian and Latino populations.

Supplementary Figure 1. Principal components analysis of European ancestry in the African American, Native Hawaiian and Latino populations. Supplementary Figure. Principal components analysis of European ancestry in the African American, Native Hawaiian and Latino populations. a Eigenvector 2.5..5.5. African Americans European Americans e

More information

The Role of Host: Genetic Variation

The Role of Host: Genetic Variation The Role of Host: Genetic Variation Patrick J. Stover, PhD The Janet and Gordon Lankton Professor of Nutrition Director, Division of Nutritional Sciences, Cornell University Dietary Requirements are Complex

More information

A Bayesian approach for constructing genetic maps when markers are miscoded

A Bayesian approach for constructing genetic maps when markers are miscoded Genet. Sel. Evol. 34 (2002) 353 369 353 INRA, EDP Sciences, 2002 DOI: 10.1051/gse:2002012 Original article A Bayesian approach for constructing genetic maps when markers are miscoded Guilherme J.M. ROSA

More information

Department of Psychiatry, Henan Mental Hospital, The Second Affiliated Hospital of Xinxiang Medical University, Xinxiang, China

Department of Psychiatry, Henan Mental Hospital, The Second Affiliated Hospital of Xinxiang Medical University, Xinxiang, China Association study of suppressor with morphogenetic effect on genitalia protein 6 (SMG6) polymorphisms and schizophrenia symptoms in the Han Chinese population Hongyan Yu 1,, Yongfeng Yang 1,2,3,, Wenqiang

More information

Genome-wide association studies for human narcolepsy and other complex diseases

Genome-wide association studies for human narcolepsy and other complex diseases ETHZ-JST Joint WS Sept. 15-16, 2008 Genome-wide association studies for human narcolepsy and other complex diseases Katsushi Tokunaga Department of Human Genetics Graduate School of Medicine The University

More information

Statistical questions for statistical methods

Statistical questions for statistical methods Statistical questions for statistical methods Unpaired (two-sample) t-test DECIDE: Does the numerical outcome have a relationship with the categorical explanatory variable? Is the mean of the outcome the

More information

Antigen Presentation to T lymphocytes

Antigen Presentation to T lymphocytes Antigen Presentation to T lymphocytes Immunology 441 Lectures 6 & 7 Chapter 6 October 10 & 12, 2016 Jessica Hamerman jhamerman@benaroyaresearch.org Office hours by arrangement Antigen processing: How are

More information

Title: Pinpointing resilience in Bipolar Disorder

Title: Pinpointing resilience in Bipolar Disorder Title: Pinpointing resilience in Bipolar Disorder 1. AIM OF THE RESEARCH AND BRIEF BACKGROUND Bipolar disorder (BD) is a mood disorder characterised by episodes of depression and mania. It ranks as one

More information

ASSOCIATION OF KCNJ1 VARIATION WITH CHANGE IN FASTING GLUCOSE AND NEW ONSET DIABETES DURING HCTZ TREATMENT

ASSOCIATION OF KCNJ1 VARIATION WITH CHANGE IN FASTING GLUCOSE AND NEW ONSET DIABETES DURING HCTZ TREATMENT ONLINE SUPPLEMENT ASSOCIATION OF KCNJ1 VARIATION WITH CHANGE IN FASTING GLUCOSE AND NEW ONSET DIABETES DURING HCTZ TREATMENT Jason H Karnes, PharmD 1, Caitrin W McDonough, PhD 1, Yan Gong, PhD 1, Teresa

More information

MENDELIAN GENETICS. MENDEL RULE AND LAWS Please read and make sure you understand the following instructions and knowledge before you go on.

MENDELIAN GENETICS. MENDEL RULE AND LAWS Please read and make sure you understand the following instructions and knowledge before you go on. MENDELIAN GENETICS Objectives Upon completion of this lab, students should: 1. Understand the principles and terms used in Mendelian genetics. 2. Know how to complete a Punnett square to estimate phenotypic

More information

Supplementary Figures

Supplementary Figures Supplementary Figures Supplementary Fig 1. Comparison of sub-samples on the first two principal components of genetic variation. TheBritishsampleisplottedwithredpoints.The sub-samples of the diverse sample

More information

Correlations between the COMT gene rs4680 polymorphism and susceptibility to ovarian cancer

Correlations between the COMT gene rs4680 polymorphism and susceptibility to ovarian cancer Correlations between the COMT gene rs4680 polymorphism and susceptibility to ovarian cancer W. Pan 1 and H. Liao 2 1 Department of Obstetrics and Gynecology, Huangshi Central Hospital of Hubei Province

More information

Selection at one locus with many alleles, fertility selection, and sexual selection

Selection at one locus with many alleles, fertility selection, and sexual selection Selection at one locus with many alleles, fertility selection, and sexual selection Introduction It s easy to extend the Hardy-Weinberg principle to multiple alleles at a single locus. In fact, we already

More information

Global variation in copy number in the human genome

Global variation in copy number in the human genome Global variation in copy number in the human genome Redon et. al. Nature 444:444-454 (2006) 12.03.2007 Tarmo Puurand Study 270 individuals (HapMap collection) Affymetrix 500K Whole Genome TilePath (WGTP)

More information

Animal Science USAMV Iaşi; University of Agricultural Sciences and Veterinary Medicine Ion Ionescu de la Brad Iaşi ; Institute of Life Sciences at

Animal Science USAMV Iaşi; University of Agricultural Sciences and Veterinary Medicine Ion Ionescu de la Brad Iaşi ; Institute of Life Sciences at PhD THESIS ABSTRACT The PhD thesis entitled Genetic polymorphism of milk proteins and the relationships between their genetic variants with quantitative and qualitative traits in some cattle populations

More information

Relationship between genetic polymorphisms in the DRD5 gene and paranoid schizophrenia in northern Han Chinese

Relationship between genetic polymorphisms in the DRD5 gene and paranoid schizophrenia in northern Han Chinese Relationship between genetic polymorphisms in the DRD5 gene and paranoid schizophrenia in northern Han Chinese Y. Zhao, M. Ding, H. Pang, X.M. Xu and B.J. Wang Department of Forensic Medicine, China Medical

More information

SNPrints: Defining SNP signatures for prediction of onset in complex diseases

SNPrints: Defining SNP signatures for prediction of onset in complex diseases SNPrints: Defining SNP signatures for prediction of onset in complex diseases Linda Liu, Biomedical Informatics, Stanford University Daniel Newburger, Biomedical Informatics, Stanford University Grace

More information

Mendelian Genetics using Fast Plants Report due Sept. 15/16. Readings: Mendelian genetics: Hartwell Chapter 2 pp , Chapter 5 pp

Mendelian Genetics using Fast Plants Report due Sept. 15/16. Readings: Mendelian genetics: Hartwell Chapter 2 pp , Chapter 5 pp 1 Biology 423L Sept. 1/2 Mendelian Genetics using Fast Plants Report due Sept. 15/16. Readings: Mendelian genetics: Hartwell Chapter 2 pp. 13-27, Chapter 5 pp. 127-130. FASTPLANTS: Williams et al. (1986)

More information

Hellinger Optimal Criterion and HPA- Optimum Designs for Model Discrimination and Probability-Based Optimality

Hellinger Optimal Criterion and HPA- Optimum Designs for Model Discrimination and Probability-Based Optimality International Journal of Mathematics and Statistics Invention (IJMSI) E-ISSN: 2321 4767 P-ISSN: 2321-4759 Volume 5 Issue 2 February. 2017 PP-71-77 Hellinger Optimal Criterion and HPA- Optimum Designs for

More information

Genetic variants on 17q21 are associated with asthma in a Han Chinese population

Genetic variants on 17q21 are associated with asthma in a Han Chinese population Genetic variants on 17q21 are associated with asthma in a Han Chinese population F.-X. Li 1 *, J.-Y. Tan 2 *, X.-X. Yang 1, Y.-S. Wu 1, D. Wu 3 and M. Li 1 1 School of Biotechnology, Southern Medical University,

More information

Population Genetics 4: Assortative mating

Population Genetics 4: Assortative mating opulation Genetics 4: Assortative mating Mating system Random ositive assortment Negative assortment Inbreeding Mate choice is independent of both phenotype and genotype Mate choice is based on similarity

More information

Original Article. C18orf1 located on chromosome 18p11.2 may confer susceptibility to schizophrenia

Original Article. C18orf1 located on chromosome 18p11.2 may confer susceptibility to schizophrenia J Med Dent Sci 2003; 50: 225 229 Original Article C18orf1 located on chromosome 18p11.2 may confer susceptibility to schizophrenia Mika Kikuchi 1,2, Kazuo Yamada 1, Tomoko Toyota 1,2 and Takeo Yoshikawa

More information

Supplementary Online Content

Supplementary Online Content Supplementary Online Content Lotta LA, Stewart ID, Sharp SJ, et al. Association of genetically enhanced lipoprotein lipase mediated lipolysis and low-density lipoprotein cholesterol lowering alleles with

More information

Nature Genetics: doi: /ng Supplementary Figure 1. Study design.

Nature Genetics: doi: /ng Supplementary Figure 1. Study design. Supplementary Figure 1 Study design. Leukopenia was classified as early when it occurred within the first 8 weeks of thiopurine therapy and as late when it occurred more than 8 weeks after the start of

More information

Unequal Numbers of Judges per Subject

Unequal Numbers of Judges per Subject The Reliability of Dichotomous Judgments: Unequal Numbers of Judges per Subject Joseph L. Fleiss Columbia University and New York State Psychiatric Institute Jack Cuzick Columbia University Consider a

More information

The Making of the Fittest: Natural Selection in Humans

The Making of the Fittest: Natural Selection in Humans INTRODUCTION MENDELIAN GENETICS, PROBABILITY, PEDIGREES, AND CHI-SQUARE STATISTICS Hemoglobin is a protein found in red blood cells (RBCs) that transports oxygen throughout the body. The hemoglobin protein

More information

Tutorial on Genome-Wide Association Studies

Tutorial on Genome-Wide Association Studies Tutorial on Genome-Wide Association Studies Assistant Professor Institute for Computational Biology Department of Epidemiology and Biostatistics Case Western Reserve University Acknowledgements Dana Crawford

More information

Myoglobin A79G polymorphism association with exercise-induced skeletal muscle damage

Myoglobin A79G polymorphism association with exercise-induced skeletal muscle damage Myoglobin A79G polymorphism association with exercise-induced skeletal muscle damage T. Cui and M.S. Jiang College of Physical Education, Shandong University of Finance and Economics, Ji nan, Shandong,

More information

Combining Different Marker Densities in Genomic Evaluation

Combining Different Marker Densities in Genomic Evaluation Combining Different Marker Densities in Genomic Evaluation 1, Jeff O Connell 2, George Wiggans 1, Kent Weigel 3 1 Animal Improvement Programs Lab, USDA, Beltsville, MD, USA 2 University of Maryland School

More information

Association of Single Nucleotide Polymorphisms (SNPs) in CCR6, TAGAP and TNFAIP3 with Rheumatoid Arthritis in African Americans

Association of Single Nucleotide Polymorphisms (SNPs) in CCR6, TAGAP and TNFAIP3 with Rheumatoid Arthritis in African Americans Association of Single Nucleotide Polymorphisms (SNPs) in CCR6, TAGAP and TNFAIP3 with Rheumatoid Arthritis in African Americans Elizabeth A. Perkins, University of Alabama at Birmingham Dawn Landis, University

More information

TADA: Analyzing De Novo, Transmission and Case-Control Sequencing Data

TADA: Analyzing De Novo, Transmission and Case-Control Sequencing Data TADA: Analyzing De Novo, Transmission and Case-Control Sequencing Data Each person inherits mutations from parents, some of which may predispose the person to certain diseases. Meanwhile, new mutations

More information

Drosophila melanogaster. Introduction. Drosophila melanogaster is a kind of flies fruit fly that is widely used in genetic

Drosophila melanogaster. Introduction. Drosophila melanogaster is a kind of flies fruit fly that is widely used in genetic Jessie Tran Mrs. Lajoie Honors Biology Date of Experiment: 4 May 2015 Due Date: 12 May 2015 Determining the Inheritance Patterns of Purple Eyes, Lobe Eyes, and Yellow Body Genes of Drosophila melanogaster

More information

Systems of Mating: Systems of Mating:

Systems of Mating: Systems of Mating: 8/29/2 Systems of Mating: the rules by which pairs of gametes are chosen from the local gene pool to be united in a zygote with respect to a particular locus or genetic system. Systems of Mating: A deme

More information

Mendelian & Complex Traits. Quantitative Imaging Genomics. Genetics Terminology 2. Genetics Terminology 1. Human Genome. Genetics Terminology 3

Mendelian & Complex Traits. Quantitative Imaging Genomics. Genetics Terminology 2. Genetics Terminology 1. Human Genome. Genetics Terminology 3 Mendelian & Complex Traits Quantitative Imaging Genomics David C. Glahn, PhD Olin Neuropsychiatry Research Center & Department of Psychiatry, Yale University July, 010 Mendelian Trait A trait influenced

More information

The association between TCM syndromes and SCAP polymorphisms in subjects with non-alcoholic fatty liver disease

The association between TCM syndromes and SCAP polymorphisms in subjects with non-alcoholic fatty liver disease The association between TCM syndromes and SCAP polymorphisms in subjects with non-alcoholic fatty liver disease Shanshan Sun, Tao Wu, Miao Wang, Wei Li, Lin Wang, Songhua He, Huafeng Wei, Haiyan Song,

More information

COMBINING PHENOTYPIC AND GENOTYPIC DATA TO DISCOVER MULTIPLE DISEASE GENES. Hannu Toivonen, Saara Hyvönen, Petteri Sevon

COMBINING PHENOTYPIC AND GENOTYPIC DATA TO DISCOVER MULTIPLE DISEASE GENES. Hannu Toivonen, Saara Hyvönen, Petteri Sevon COMBINING PHENOTYPIC AND GENOTYPIC DATA TO DISCOVER MULTIPLE DISEASE GENES Hannu Toivonen, Saara Hyvönen, Petteri Sevon Department of Computer Science, University of Helsinki P.O.Box 68, FI 4 University

More information

MULTIFACTORIAL DISEASES. MG L-10 July 7 th 2014

MULTIFACTORIAL DISEASES. MG L-10 July 7 th 2014 MULTIFACTORIAL DISEASES MG L-10 July 7 th 2014 Genetic Diseases Unifactorial Chromosomal Multifactorial AD Numerical AR Structural X-linked Microdeletions Mitochondrial Spectrum of Alterations in DNA Sequence

More information

Rare Variant Burden Tests. Biostatistics 666

Rare Variant Burden Tests. Biostatistics 666 Rare Variant Burden Tests Biostatistics 666 Last Lecture Analysis of Short Read Sequence Data Low pass sequencing approaches Modeling haplotype sharing between individuals allows accurate variant calls

More information

PopGen4: Assortative mating

PopGen4: Assortative mating opgen4: Assortative mating Introduction Although random mating is the most important system of mating in many natural populations, non-random mating can also be an important mating system in some populations.

More information

Using Imputed Genotypes for Relative Risk Estimation in Case-Parent Studies

Using Imputed Genotypes for Relative Risk Estimation in Case-Parent Studies American Journal of Epidemiology Published by Oxford University Press on behalf of the Johns Hopkins Bloomberg School of Public Health 011. Vol. 173, No. 5 DOI: 10.1093/aje/kwq363 Advance Access publication:

More information

The Association Design and a Continuous Phenotype

The Association Design and a Continuous Phenotype PSYC 5102: Association Design & Continuous Phenotypes (4/4/07) 1 The Association Design and a Continuous Phenotype The purpose of this note is to demonstrate how to perform a population-based association

More information

Evaluating the role of the 620W allele of protein tyrosine phosphatase PTPN22 in Crohn s disease and multiple sclerosis

Evaluating the role of the 620W allele of protein tyrosine phosphatase PTPN22 in Crohn s disease and multiple sclerosis (2006) 14, 317 321 & 2006 Nature Publishing Group All rights reserved 1018-4813/06 $30.00 www.nature.com/ejhg ARTICLE Evaluating the role of the 620W allele of protein tyrosine phosphatase PTPN22 in Crohn

More information