How To Build The Index Between Genome And Mrnas?
3
2
Entering edit mode
11.3 years ago
Haiping ▴ 110

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 • 2.8k views
ADD COMMENT
3
Entering edit mode
11.3 years ago

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

ADD COMMENT
1
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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!

ADD REPLY
3
Entering edit mode
11.3 years ago
Rm 8.2k

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.

ADD COMMENT
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

edited my answer to add grape annotation information

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode
11.3 years ago

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

ADD COMMENT
1
Entering edit mode

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.

ADD REPLY

Login before adding your answer.

Traffic: 1233 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6