2018
DOI: 10.1101/350355
|View full text |Cite
Preprint
|
Sign up to set email alerts
|

FunSPU: a versatile and adaptive multiple functional annotation-based association test of whole-genome sequencing data

Abstract: Despite ongoing large-scale population-based whole-genome sequencing (WGS) projects such as the NIH NHLBI TOPMed program and the NHGRI Genome Sequencing Program, WGS-based association analysis of complex traits remains a tremendous challenge due to the large number of rare variants, many of which are non-trait-associated neutral variants. External biological knowledge, such as functional annotations based on ENCODE, may be helpful in distinguishing causal rare variants from neutral ones; however, each function… Show more

Help me understand this report
View published versions

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
2
1

Citation Types

0
13
0

Year Published

2019
2019
2020
2020

Publication Types

Select...
6

Relationship

2
4

Authors

Journals

citations
Cited by 10 publications
(13 citation statements)
references
References 35 publications
(21 reference statements)
0
13
0
Order By: Relevance
“…First, as transcriptome-wide association studies (TWAS) ( Gamazon et al, 2015 ; Gusev et al, 2016 ) that incorporate eQTL-derived weights into a weighted Sum test ( Xu et al, 2017 ) to both improve statistical power and enhance biological interpretation, our proposed method can incorporate eQTL-derived weights into the test statistics of iSPU( γ ) and aiSPU. Also, some other functional weights ( He et al, 2017 ; Ma and Wei, 2019 ) can be equally applied. We expect that integrating functional genomic information will improve power and gain insights into the mechanisms of complex traits.…”
Section: Discussionmentioning
confidence: 99%
“…First, as transcriptome-wide association studies (TWAS) ( Gamazon et al, 2015 ; Gusev et al, 2016 ) that incorporate eQTL-derived weights into a weighted Sum test ( Xu et al, 2017 ) to both improve statistical power and enhance biological interpretation, our proposed method can incorporate eQTL-derived weights into the test statistics of iSPU( γ ) and aiSPU. Also, some other functional weights ( He et al, 2017 ; Ma and Wei, 2019 ) can be equally applied. We expect that integrating functional genomic information will improve power and gain insights into the mechanisms of complex traits.…”
Section: Discussionmentioning
confidence: 99%
“…In our data example, gene B4GALNT2 was discovered by the aSPU-meta test for the AA subjects in the presence of high between-study heterogeneity and a large number of likely neutral RVs; in contrast, it was missed by competing methods. Although we did not focus on the use of variantspecific or study-specific weights, appropriate incorporation of such weights may increase the statistical power and interpretation, for example, using the functional weights (He, Xu, Lee, & Ionita-Laza, 2017;Ma & Wei, 2019;Zhan & Liu, 2015). We also note that, similar to other competing tests for RVs, one potential drawback of the current implementation of the aSPU-meta test using Monte-Carlo simulations is its use of the asymptotic normal distribution of the score vector, which may not hold for extremely unbalanced case-control studies, and other resampling methods like permutation-based ones may need to be adopted for more accurate p value calculations (Dey et al, 2019).…”
Section: Discussionmentioning
confidence: 99%
“…The inner sum of the weighted score function with power γ1 can be regarded as the aGE test statistic for weighted genotype in each tissue, and they are normalized to the power of 1γ1 before taking a second power γ2. Based on our previous investigations (Kwak & Pan, 2016; Ma & Wei, 2019; Pan et al, 2014) and simulation study to be detailed later, we would recommend using Γ1={1,2,3,4,5,6} and Γ2={1,2,4}, which suffice to strike a balance between statistical power and computational complexity. When γ1=1, the GEw(γ1=1,) test is equivalent to the adaptive weighted burden test, assuming that the interaction effects of SNPs are proportional to the external weights, that is the full model for the m th tissue is reduced to h(μi)=XiβX+βeEi+jqβjGij+α(m)jqwj(m)Sij. Note that formulation () avoids the potential problem of an inflated Type 1 error rate in model ().…”
Section: Methodsmentioning
confidence: 99%