This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What is the best way to retrieve gene name from variant mutation?

Hi all! I was wondering if any of you may suggest a clean a fast way to retrieve gene name from a list of variants that look like this:

CHR_POS_REF_ALT > chr1_33333_A_C

Of course I can generate a file with POS, CHR, REF, ALT and as I am missing some information for a standard VCF file, I was wondering if "VariantAnnotation" from Bioconductor is the correct tool to use.

mutation name gene variant
• 3,300 views • link
written by Lila M Lila M

You only really need POS, CHR, REF, and ALT for a VCF file. 8 columns are needed but everything else can be a dot (i.e., .).

• 0 views • link
written by lethalfang • lethalfang

1 answer 

tabix tabix_indexed.gtf.gz ` echo "chr1_2077466_A_C" | awk -F '_' '{printf("%s:%s-%s\n",$1,$2,$2);}'` |\
awk -F '\t' '($3=="gene")' | cut -f 9 | tr ";" "\n" |\
grep -w gene_name | cut -d '"' -f 2 |\
sort | uniq | paste -sd,
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thank you for your feedback! I've never used tabix before, am I right it I assume tabix_indexed.gtf.gz cames with the software?

• 0 views • link
written by Lila M Lila M

ah but I don't have the gtf file

• 0 views • link
written by Lila M Lila M

If this is for human (or another model organism) you should be able to get a file based on genome build.

• 0 views • link
written by GenoMax GenoMax

Thank you both for your feedback! what I am doing so far is:

gunzip  Homo_sapiens.GRCh38.109.gtf.gz | awk '$1 ~ /^#/ {print $0;next} {print $0 | "sort -k1,1 -k4,4n -k5,5n"}' | bgzip > Homo_sapiens.GRCh38.109_sorted.gtf.gz

tabix Homo_sapiens.GRCh38.109_sorted.gtf.gz

tabix Homo_sapiens.GRCh38.109_sorted.gtf.gz  ` echo "chr1_1000000_G_A" | awk -F '_' '{printf("%s:%s-%s\n",$1,$2,$2);}'` |awk -F '\t' '($3=="gene")' | cut -f 9 | tr ";" "\n" |grep -w gene_name | cut -d '"' -f 2 |sort | uniq | paste -sd,

But unfortunately it raises this

usage: paste [-s] [-d delimiters] file ...

Any idea what is causing it?

• 0 views • link
written by Lila M Lila M

Does not work but does not produce any error either. I am not understanding tabix followed by the echo. It is possible the biostar formatting may have eliminated something in your command let me check.

Update: I got a copy of the original command from your answer and it produced no output.

• 0 views • link
written by GenoMax GenoMax

Pierre Lindenbaum Is the command you posted as answer correct/complete? Can you please check? I can't get it to work. What is the expected output?

• 0 views • link
written by GenoMax GenoMax

Hello again, did we reach any conclusion about this? I really appreciate it :) Thank you!

• 0 views • link
written by Lila M Lila M

Log in to answer this question.

More posts like this