This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Build The Index Between Genome And Mrnas?

Hi. I finished calling SNP from genome resequencing data by using samtools. Now I need to know where the SNP located(5'-UTR, exon, intron, 3'UTR, et al). So is there any softwere I can use to build a index between genome and annotated mRNA so that I can calculate the variation of SNP in genes?

index genome mrna snp variant
• 5,250 views • link
updated by 2184687-1231-83- 2184687-1231-83- • written by Haiping • Haiping

3 answers

Using your SNP coordinates, if your genome is in Ensembl, you could use the Ensembl's variant effect predictor (VEP) tool: http://www.ensembl.org/tools.html

Since you are working on grape, you can use http://plants.ensembl.org/Vitis_vinifera/Info/Index

• 1 views • link
written by 2184687-1231-83- 2184687-1231-83-

Thanks for you reply. Since I worked on grape genomes, there is no information for grape on Ensembl's. I do hava the grape genomes and predicted fasta mRNAs.So is there any softwere that can output a .gff files by using genome and mRNAs? Thanks again!

• 0 views • link
written by Haiping • Haiping

Ok. Thanks. I never found it before. i will try to use VEP tools. Thanks a lot for you help.

• 0 views • link
written by Haiping • Haiping

use bedtools: to know/map the SNPs to the genic regions

Simple example can be: 

intersectBed -a snp.bed -b exons.bed > snpsInExons.bed

you can get the gene/exon or other information from ucsc:

ADDED: For grape genome look annotation files;
you can extract required information from the GFF file provided.

• 1 views • link
updated by Ram Ram • written by Rm Rm

Thank you for you reply. I checked the information from UCSC.Since I worked on grape, there is no information related to plant on UCSC. Is there another website can do this kind of job or is there another sofewere i can used to do it myselfs? Thanks again.

• 0 views • link
written by Haiping • Haiping

edited my answer to add grape annotation information

• 0 views • link
written by Rm Rm

I got it! Thanks a lot. Actually I saw that there is a files but I don't know it's a GFF files I needed. Now I had Gff file and SNP(.txt) files, It's enough to use bedtools to map the snps to the genomes? I never used the bedtools before. So I need to learn it.

• 0 views • link
written by Haiping • Haiping

bed tools is simple to use: requires chr start end other-column (exon infor etc)

• 0 views • link
written by Rm Rm

OK. Thanks a lot! I will try to do it!

• 0 views • link
written by Haiping • Haiping

In www.ensemblgenomes.org, go to Ensembl Plants, and choose your species from the drop-down box or "View list of all Ensembl Plants species":

alt text

Then go to "Manage your data" and choose "Variant Effect Predictor":

alt text

You can now copy+paste your SNPs into the form or upload your SNPs file. You can also use the Perl script for processing large batches:

alt text

• 1 views • link
written by 2184687-1231-83- 2184687-1231-83-

Thanks for your explanation step by step. Now I know how to use it. Acturlly I am just analysis subset of genes and not the whole genomes, so that would be more convenient if I can do the idex myself by using the DNA and mRNAs for genes I analysised.It seems that I can extract the informations I needed from .GFF. I will try to do it. Thanks for all your comments.

• 0 views • link
written by Haiping • Haiping

Log in to answer this question.

More posts like this