Hello!
I'm trying to identify gene names for the list of SNP (GWAS result) and have some inconsistencies. For some SNPs gene names were identified using BioMart package in R (by BP - if it is located inside the gene and #chr - the same as gene). For these SNPs I make a plot in LocusZoom and get a list of genes. And not always gene which was identified with BioMart is included in the list of genes from LocusZoom.
Could you please explain, what is the reason of this inconsistencies? And can I consider genes identified by BP and chr# using BioMart correct?
Thank you!
•
written
by
alioniya •
1 answer
Two possible things to consider:
What reference genome does LocusZoom use? What reference genome did you use with biomaRt? What reference genome did you call your variants on? If any of these don't match, you're going to be looking at completely different genomic regions.
What gene set are LocusZoom using? BiomaRt will give you the latest genes from GENCODE. If LocusZoom are using a lower quality or older gene set, then they may give you less accurate gene mapping.
Log in to answer this question.
-
Mapping SNP rsIds/ positions to genes coding for proteins
written by Harsh •
-
How to identify carriers from a SNps listing?
written by itziarlamiquiz •
-
How to get all genes and their annotations within a LD block?
written by b.g.tamang •
-
Ensembl IDs as input in ReactomePA package in R
written by sainikeshav6 •
-
Generate SNP/indel annotation
written by oussriradnane •
-
using biomaRt in R to annotate SNPs to Genes within 1000kb
written by Sheila
-
Identification of SNPs in regulatory regions
written by valerie
-
Getting the function of genes based on gene names
written by GabrielMontenegro
-
Plot with locuszoom and chr:pos SNPs
written by alesssia
-
What is the first exon in negative strand genes?
written by Zakaria Benmounah •
