chr1swallace / coloc

Repo for the R package coloc
152 stars 43 forks source link

Colocalization of eQTL and GWAS #114

Open cy6n opened 1 year ago

cy6n commented 1 year ago

Hi Chris,

Thanks for making this fantastic package!

I am working on the colocalization analysis with fine mapping using eQTL and GWAS. And I got some questions when I use this package.

Currently, I have the fine-mapping results of eQTL coming from our collaborators. For example, for one specific gene, we have a credible set including two putative causal variants. And then I did my own fine mapping using SuSIE (susie_rss) for this region (position of this gene +/- 1MB) at GWAS and I got one credible set including 25 putative causal variants.

However, there is only one common/overlap SNP at the fine-mapping result of GWAS and eQTL. Could I use this one specific SNP to do the coloc.abf? And if I got two common/overlap SNPs from the fine-mapping of GWAS and eQTL. Could I use these two SNPs to conduct the coloc.susie by using their LD martix?

Greatly appreciate your help!

Best, Chaojie

mocksu commented 1 year ago

Did you check the genomic coordinates of your datasets? They might be in different versions (e.g. hg37 vs hg38).

chr1swallace commented 1 year ago

coloc evaluates different hypotheses including no causal variants for either trait, so you need statistics for variants that cover the whole region of interest in both traits, not just the intersection of independent fine mapping.

cy6n commented 1 year ago

Did you check the genomic coordinates of your datasets? They might be in different versions (e.g. hg37 vs hg38). Thanks for your reply.

Yes -- they are in same version (hg38)

cy6n commented 1 year ago

coloc evaluates different hypotheses including no causal variants for either trait, so you need statistics for variants that cover the whole region of interest in both traits, not just the intersection of independent fine mapping.

Thanks for your responses and that is very helpful. However, I was wondering why we need to apply the color_susie if coloc evaluates different hypotheses including no causal variants for either trait?

chr1swallace commented 1 year ago

Because coloc.abf doesn't allow for more than one casual variant per trait. Coloc.susie allows for that

https://chr1swallace.github.io


From: cy6n @.> Sent: Thursday, March 9, 2023 9:58:49 PM To: chr1swallace/coloc @.> Cc: Chris Wallace @.>; Comment @.> Subject: Re: [chr1swallace/coloc] Colocalization of eQTL and GWAS (Issue #114)

coloc evaluates different hypotheses including no causal variants for either trait, so you need statistics for variants that cover the whole region of interest in both traits, not just the intersection of independent fine mapping.

Thanks for your responses and that is very helpful. However, I was wondering why we need to apply the color_susie if coloc evaluates different hypotheses including no causal variants for either trait?

— Reply to this email directly, view it on GitHubhttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fchr1swallace%2Fcoloc%2Fissues%2F114%23issuecomment-1462881641&data=05%7C01%7Ccew54%40universityofcambridgecloud.onmicrosoft.com%7C12dfb0aaf4de471f246d08db20e977c0%7C49a50445bdfa4b79ade3547b4f3986e9%7C1%7C0%7C638139959325356916%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=Tbq3Y9bJ6pMkB2x%2FW1ZKs8t%2BpT9lYRH2qcwHq6RMazg%3D&reserved=0, or unsubscribehttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FAAQWR2ATN2DNNADEM2ERVZTW3JHBTANCNFSM6AAAAAAUYHSCGU&data=05%7C01%7Ccew54%40universityofcambridgecloud.onmicrosoft.com%7C12dfb0aaf4de471f246d08db20e977c0%7C49a50445bdfa4b79ade3547b4f3986e9%7C1%7C0%7C638139959325356916%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=buSpzgWizC%2FZIThmF95l7x45CI%2BiOMxaWy14ox1a%2BSM%3D&reserved=0. You are receiving this because you commented.Message ID: @.***>