Question: findings start and stop position of a intron
1
gravatar for avari
3.2 years ago by
avari80
Germany
avari80 wrote:

Hello,

I want to find the start and stop position of a particular intron for a human gene and then exact SNPs in my dataset that fall within this window using Plink. Where would be a good place to find this information ?

Thanks!

intron plink position snps gene • 1.1k views
ADD COMMENTlink written 3.2 years ago by avari80
3

I would do following

  • Create a BED file with Introns only (Like this)
  • BED file containing SNP coordiantes
  • Then use bedtools intersect function on these two BED files to find out SNPs falling in intronic regions.
ADD REPLYlink written 3.2 years ago by venu6.2k

Thanks for the instructions, that seems useful but a little complicated! I'm afraid I don’t have a bioinformatics background so I am unfamiliar with some of these tools. I take it there is no quick and dirty way of searching for this in some database? Just looking for the start/stop base pair positions.

ADD REPLYlink written 3.2 years ago by avari80
1

My suggestion would be to get familiar with these tools & UNIX environment. I can't say its not possible to do these tasks without learning those tools (Though I'm a bit biased to the end of 'not possible'), it would be much more easy to do such things if you are familiar.

And take a look at this link: How do I get the position of a SNP when it is located in an intron?. It may help but I am not sure whether its a perfect solution.

ADD REPLYlink written 3.2 years ago by venu6.2k
1

If you have the SNP coordinates and the alleles e.g. 15 2233 2233 A/G (ref/alt), you can use the Variant Effect Predictor to find out where your SNPs map to (intronic regions or not) and what is their effect on proteins, transcripts, regulatory elements, etc. If your species is not in Ensembl, you could go about building a cache from GTF or GFF file..

ADD REPLYlink written 3.2 years ago by Denise - Open Targets5.0k

Thanks, that might work actually as I can extract all the SNPs for the gene in Plink first and then feed them in. Will give it a try. I am specifically interested in the third intron, so hopefully the intron number is specified.

ADD REPLYlink written 3.2 years ago by avari80
1

You can do this without programming using UCSC Genome Browser. Navigate to the region of interest (your intron). Then go to Tools > Table Browser. Under "region", select "position" to limit yourself to that region. Under "group" and "track", select the appropriate SNP track.

ADD REPLYlink written 3.2 years ago by igor8.1k

Thanks, that sounds useful too!

ADD REPLYlink written 3.2 years ago by avari80
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 619 users visited in the last hour