Difference between revisions of "Pascal"
(8 intermediate revisions by 2 users not shown) | |||
Line 6: | Line 6: | ||
<date> Jan 25 </date> | <date> Jan 25 </date> | ||
</teaser> | </teaser> | ||
− | |||
== Rigorous gene and pathway analysis of GWAS == | == Rigorous gene and pathway analysis of GWAS == | ||
Line 13: | Line 12: | ||
''Download'' | ''Download'' | ||
− | * '''[http://www2.unil.ch/cbg/images/3/3d/PASCAL.zip Pascal package]''' (Download might take a while because the 1KG-EUR data | + | * '''[http://www2.unil.ch/cbg/images/3/3d/PASCAL.zip Pascal package]''' (Download might take a while because the 1KG-EUR data are included) |
− | * '''[[PascalTestData | Test data]]''' (Additional data that | + | * '''[[PascalTestData | Test data]]''' (Additional data that were used for evaluation in the paper) |
''Reference'' | ''Reference'' | ||
− | * '''Fast and rigorous computation of gene and pathway scores from SNP-based summary statistics'''. ('''[http://regulatorycircuits.org/data/papers/LamparterMarbach2016.pdf PDF]''')<br> Lamparter D*, Marbach D*, | + | * '''Fast and rigorous computation of gene and pathway scores from SNP-based summary statistics'''. ('''[http://regulatorycircuits.org/data/papers/LamparterMarbach2016.pdf PDF]''')<br> Lamparter D*, Marbach D*, Rueedi R, Kutalik Z, and Bergmann S. <br> [http://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1004714 ''PLoS Computational Biology'' 12, e1004714, 2016.] |
− | |||
[[File:PascalFigure1.jpg|500px|left]] | [[File:PascalFigure1.jpg|500px|left]] | ||
'''Figure: Overview of methodology to compute gene and pathway scores''' | '''Figure: Overview of methodology to compute gene and pathway scores''' | ||
Line 26: | Line 24: | ||
('''b''') We use external databases to define gene sets for each reported pathway. We then compute pathway scores by combining the scores of genes that belong to the same pathways, i.e. gene sets. The fast gene scoring method allows us to dynamically recalculate gene scores by aggregating SNP p-values across pathway genes that are in LD and thus cannot be treated independently. This amounts to fusing the genes and computing a new score that takes the full LD structure of the corresponding locus into account. We evaluate pathway enrichment of high-scoring (potentially fused) genes using one of two parameter-free procedures (chi-square or empirical score), avoiding any p-value cutoffs inherent to standard binary enrichment tests. | ('''b''') We use external databases to define gene sets for each reported pathway. We then compute pathway scores by combining the scores of genes that belong to the same pathways, i.e. gene sets. The fast gene scoring method allows us to dynamically recalculate gene scores by aggregating SNP p-values across pathway genes that are in LD and thus cannot be treated independently. This amounts to fusing the genes and computing a new score that takes the full LD structure of the corresponding locus into account. We evaluate pathway enrichment of high-scoring (potentially fused) genes using one of two parameter-free procedures (chi-square or empirical score), avoiding any p-value cutoffs inherent to standard binary enrichment tests. | ||
+ | <br><br> | ||
+ | |||
+ | == Tissue-specific regulatory circuits disrupted in complex disease == | ||
+ | |||
+ | [[File:Pascal_network_analysis.png|260px|right]] | ||
+ | '''The efficiency and accuracy of Pascal opens the door to large-scale analyses''' that would not have been possible with previous tools. For example, summarizing SNP p-values at the level of genes is a crucial step in most network-based GWAS analysis methods. Pascal was key for our recent work, where we integrated '''37 GWAS datasets''' with close to '''400 tissue-specific gene regulatory circuits''' to systematically analyze the inter-connectivity of genes that are perturbed by trait-associated genetic variants. This study showed that disease-associated genetic variants often disturb regulatory modules in cell types or tissues that are highly specific to that disease, giving new insights on disease mechanisms. | ||
+ | |||
+ | ''Reference'' | ||
+ | |||
+ | * '''Tissue-specific regulatory circuits reveal variable modular perturbations across complex diseases'''. ('''[http://compbio.mit.edu/publications/144_Marbach_NatureMethods_16.pdf PDF]''') <br>Marbach D, Lamparter D, Quon G, Kellis M, Kutalik Z, and Bergmann S. <br> [http://www.nature.com/nmeth/journal/v13/n4/abs/nmeth.3799.html ''Nature Methods'', 13, 366-370, 2016]. | ||
+ | **[http://dx.doi.org/10.1038/nrg.2016.36 ''Nature Reviews Genetics'' highlight], [http://regulatorycircuits.org/papers/PR_genes_social_network_en_fr.pdf SIB news], [https://www.unil.ch/getactu/wwwfbm/1457426312508/ UNIL news] | ||
+ | ** [http://regulatorycircuits.org regulatorycircuits.org] |
Latest revision as of 13:39, 13 April 2016
Rigorous gene and pathway analysis of GWAS
Pascal (Pathway scoring algorithm) is an easy-to-use tool for gene scoring and pathway analysis from GWAS results. Pascal uses external data to estimate linkage disequilibrium. Therefore, the user only needs to supply genome wide SNP p-values. Pascal then derives p-values for genes and predefined pathways. Pascal doesn’t use Monte-Carlo simulation to derive gene p-values. This leads to increased speed and accuracy. This speed in the gene scoring is then leveraged to control the false positive rate in pathway scoring. For pathway scoring we implemented and tested enrichment strategies that compared very favorably compared to hypergeometric enrichment. This comparison was done on a large collection of GWAS results giving us confidence to recommend Pascal for downstream analysis of GWAS results. Pascal is mainly written in Java and has been tested on Unix systems and Mac OsX.
Download
- Pascal package (Download might take a while because the 1KG-EUR data are included)
- Test data (Additional data that were used for evaluation in the paper)
Reference
- Fast and rigorous computation of gene and pathway scores from SNP-based summary statistics. (PDF)
Lamparter D*, Marbach D*, Rueedi R, Kutalik Z, and Bergmann S.
PLoS Computational Biology 12, e1004714, 2016.
Figure: Overview of methodology to compute gene and pathway scores
(a) We compute gene scores by aggregating SNP p-values from a GWAS meta-analysis (without the need for individual genotypes), while correcting for linkage disequilibrium (LD) structure. To this end, we use numerical and analytic solutions to compute gene p-values efficiently and accurately given LD information from a reference population. Two options are available: the max and sum of chi-squared statistics, which are based on the most significant SNP and the average association signal across the region, respectively.
(b) We use external databases to define gene sets for each reported pathway. We then compute pathway scores by combining the scores of genes that belong to the same pathways, i.e. gene sets. The fast gene scoring method allows us to dynamically recalculate gene scores by aggregating SNP p-values across pathway genes that are in LD and thus cannot be treated independently. This amounts to fusing the genes and computing a new score that takes the full LD structure of the corresponding locus into account. We evaluate pathway enrichment of high-scoring (potentially fused) genes using one of two parameter-free procedures (chi-square or empirical score), avoiding any p-value cutoffs inherent to standard binary enrichment tests.
Tissue-specific regulatory circuits disrupted in complex disease
The efficiency and accuracy of Pascal opens the door to large-scale analyses that would not have been possible with previous tools. For example, summarizing SNP p-values at the level of genes is a crucial step in most network-based GWAS analysis methods. Pascal was key for our recent work, where we integrated 37 GWAS datasets with close to 400 tissue-specific gene regulatory circuits to systematically analyze the inter-connectivity of genes that are perturbed by trait-associated genetic variants. This study showed that disease-associated genetic variants often disturb regulatory modules in cell types or tissues that are highly specific to that disease, giving new insights on disease mechanisms.
Reference
- Tissue-specific regulatory circuits reveal variable modular perturbations across complex diseases. (PDF)
Marbach D, Lamparter D, Quon G, Kellis M, Kutalik Z, and Bergmann S.
Nature Methods, 13, 366-370, 2016.