Associating Genomic Variations with Phenotypes

download report

Transcript Associating Genomic Variations with Phenotypes

Associating Genomic
Variations with
Phenotypes
Model comparison, rare variants,
and analysis pipeline
Qunyuan Zhang
Division of Statistical Genomics & Genome Institute
Washington University School of Medicine
1
Data & Question
i
1
Y
y1
X
x11
x12
... x1m
2 y2
x21
... ... ... ...
x22
...
... x2 m
... ...
n
xn 2
..
yn
Phenotypes
(quantitative,
categorical)
xn1
xnm
Genotypes:
SNP
Insertion
Deletion
Duplication
Inversion
Translocation
…
Relationship
between X and Y ?
2
Linkage & Association
i
1
x11
X
q12
x12
...
2 y2
x21
... ... ... ...
q22
...
x22
...
...
...
n
Y
y1
yn
xn1
qn 2
Genotypes
xn 2 ...
Association: (Y,X)
Linkage: (Y,Q)
Phenotype
Q is unobservable
r1 Q r2
Putative QTL
3
A Fixed-effect Mixture Model For Linkage
Commonly used in plant genetics
P1
X
P2
SNP A
SNP B
F1
r1 Q r2
F2
3
f ( yi )   P(Q j | X i , r )
j 1
1

j
 1 yi   j 2 
exp  (
) 
 2 

2
j


n
L(Y )   f ( yi )
i 1
4
A Variance-component Model For Linkage
Commonly used in human genetics
SNP A
SNP B
ΔQ
r 1 Q r2
 1

T
1
L(Y ) 
exp  (Y   ) V (Y   ) 
n/2
1/ 2
(2 ) | V |
 2

1
V  Cov(Y )  ΔQ  Δg  I
2
Q
QTL IBD
matrix
2
g
Background
IBD matrix
2
e
Diagonal unit
matrix
5
Variance-component Model
= Random-effect Linear Model
Y  μ  ZQ γ Q  Z g γ g  e
MVN (0, Δg )
MVN (0, ΔQ )
2
Q
2
g
V  ΔQ  Δg  I
2
Q
2
g
Random
effects
N (0,  e2 )
2
e
 1

T
1
L(Y ) 
exp  (Y   ) V (Y   ) 
n/2
1/ 2
(2 ) | V |
 2

1
6
From Linkage to Association
QTL effect(s)
Y  μ  ZQ γ Q  Z g γ g  e
Y  μ  Xβ  Zg γ g  e
marker
effect(s)
Linkage model
Family-based
association model
V  Δg  I
2
g
2
e
fixed
effect(s)
 1

T
1
L(Y ) 
exp  (Y    X ) V (Y    X ) 
n/2
1/ 2
(2 ) | V |
 2

1
7
A Simple Association Model
For Unrelated Subjects
Y  μ  Xβ  e
L(Y ) 
1
(2 ) n / 2 | V |1/ 2
n

i 1
1

e
V  I
2
e
 1

T
1
exp  (Y    X ) V (Y    X ) 
 2

 1 yi    X 2 
exp   (
) 
e
2
 2

8
Covariate(s):
Adjusting For Confounder(s)
Y  μ  Xβ  XC βC  e
Observed confounders: age, sex etc.
Hidden confounders: population structure
Population structure can be estimated by:
-PCA
-Clustering
-Admixture/ancestry
9
Modeling Hidden Genetic Correlation
Between Subjects
marker
fixed
effect(s)
covariate
fixed
effect(s)
Genetic
background
random effects
Y  μ  Xβ  XCβC  Z g γ g  e
V  Δg  I
2
g
2
e
Family data, pedigree => IBD matrix
Population data, hidden, marker data => IBS matrix
10
Modeling Rare Variants
Y  μ  Xβ  XCβC  Z g γ g  e
Common variants, tested individually, H0: β1=0. One p-value per variant
Y  μ  1 X1  ...
Rare variants, tested as an entire group (burden test), usually by gene
H0: β1= β2=…=βk=0 . One p-value per group of variants
Y  μ  1 X1  2 X 2  ... k X k  ...
Incorporated with variable selection, with loose criteria
β can be treated as random effects, variance components
test, can be weighted by prior information
11
Collapsing Model
Y  μ  1 X1  2 X 2  ... k X k  ...
Y  μ   X   ...
Collapsing multiple variables into one
subject X 1
X2
X3
X
1
2
0
0
0
1
0
0

1
1
3
1
0
0
1
12
Weighted Sum Model
Y  μ  1 X1  2 X 2  ... k X k  ...
k
Y  μ    ( w j X j )  ...
j 1
Y  μ   S  ...
subject X 1
1
2
3
X3
X2
S
w1  0.2 w1  0.5 w1  0.3
0.0
0
0
0
0.8
0
1
1
1
0
0
0.2
13
Weighting Variants
 Base
on allele frequency, continuous or binary(0,1) weight,
variable threshold;
 Based on function annotation/prediction;
 Based on sequencing quality (coverage, mapping quality,
genotyping quality, validated or not etc.);
 Data-driven, using both genotype and phenotype data,
learning weights (including effect directions) from data,
requiring permutation test;
 Any combination …
Grouping Variants
By gene
By transcript
By gene set / pathway
……
By exon
By protein domain
14
Modeling More Data Types
Generalized Linear (Mixed) Model
g (Y )  μ  Xβ  ... e
Link function
For binary Y, logistic model
 P(Y  1) 

g (Y )  logit (Y )  log
 1  P(Y  0) 
exp(μ  Xβ  ...  e)
P(Y  1) 
exp(μ  Xβ  ...  e)  1
15
Longitudinal Data (quantitative)
Time
Fixed effect, time as covariate
Repeated measures, random effect, correlation within subjects
16
Longitudinal Data (binary)
Time
Linear model, time as covariate
Survival analysis, CoxPH model etc.
17
Tools
SAS Procedures
REG, LOGISTIC, GENMOD, MIXED,
HPMIXED, GLIMMIX, PHREG/LIFETEST
R Functions/Packages
lm (), glm()
gee, nlme, kinship2/coxme, lme4, survival
Other Programs
SOLAR, MMAP, EMMA, EMMAX, SKAT
18
Pipeline
Input (data + options)
Job generating/submitting module
Job number controlling module
LSF bsub
job1
job2
…..
Job
N
Options.jobi => self-programmed modules (SAS, R,…)
Options.jobi => external program modules (MMAP, SKAT,..)
Result
1
…..
Result
2
Result
N
Job status monitoring module (all done ?)
Yes
no
Result summarizing module
Wait …
19
gwas.sh options.gwa
#!/bin/sh
OPFILE=$1
...
…
Pheno
Bmi
YES
Obes
YES
HD
Age
Sex
…
type
qt
covar
age,sex
program analysis
SASGLM mixed
ql
NA
SASGLM gee
ql
…
…
age
SASGLM gee
Program
SASGLM
GSTAT
MMAP
language
SAS
R
C
location
Maintainer
/dsg1/code/sas/glm.sas Q.Zhang
/dsg1/code/R/gstat.R
Q.Zhang
/dsg1/code/sas/mmap.sh J. Czajkowski
…
run
NO
[DATA]
database=SAS
genotype_dir=/dsg1/gwas/fhsgeno
genotype_file=
phenotype_file=fhs100
markerinfo_file=mapall
marker_selection=MAF>0.01
pedigree_file=pediall
subjectID=subject
pedgreeID=famid
markername=snp
…
[ANALYSIS]
phenolist_file=
pheno_list=bmi/qt
covariates=
program=SASGLM
analysis=mixed
[OUTPUT]
output_dir=/dsguser/qunyuan/fhs/bmi
output_file=
output_replace=no
[RUN]
clusterjobname=bmimixed
memsize=1000M
maxjobn=300
…
20
Thanks !
21