Non-coding region of WGS
0
0
Entering edit mode
3.8 years ago
FL512 ▴ 20

Dear all,

I am new to work on WGS/WES analysis. I have VCF files obtained from WGS already and I would like to focus on SNPs located on non-coding regions such as promoter, enhancer, and 5'/3'-UTRs, or if possible, all non-coding regions I would like to analyze.

  1. Is there any database, reference sequence (txt, bed whatever), or even script available to extract SNPs of non-coding region from my VCF files?
  2. If anyone has already posted and figured it out, please kindly let me know where the post is.

I do not want to analyze by using biased approach like focusing of the 5k -20kb upstream of genes of interest. Rather, I want to do a global analysis, therefore, I am struggling...

Thank you very much in advance. Any kind of your comments/suggestion should be helpful.

non-coding region WGS SNP • 1.5k views
ADD COMMENT
2
Entering edit mode

Extract the non-UTR exons in BED format, use bedtools complement to get the complement of the entire genome and the exons, that is the non-coding part. Intersect that file with your VCF, these are your non-coding variants, bedtools intersect is probably what you want.

ADD REPLY
1
Entering edit mode

Adding to this answer you can get the non-coding region part from your reference GTF file and then you can use bedtools intersect to get your region.

ADD REPLY
0
Entering edit mode

Thank you for letting me know. By the way, how am I able to extract non-coding regions from my reference file, which I was struggling for a couple of days. I thought I had to download the non-UTR exons from UCSC genome browser. Anyway, thank you for your help!

ADD REPLY
1
Entering edit mode

Hello, GTF/GFF file has a column called 'type' from which you filter the biotype you need. gff file specification under mRNA you will find the 'intron' attributes. I hope this helps.

ADD REPLY
0
Entering edit mode

Thank you for your help!

ADD REPLY
0
Entering edit mode

UPDATE: for someone like me, this is how I downloaded the region of interests. Bed File With Introns Only

ADD REPLY
0
Entering edit mode

Thank you very much. That was also I was thinking of but I did not know how to do it because of the lacking of knowledge and experience. I will keep you updated. Thank you again.

ADD REPLY

Login before adding your answer.

Traffic: 1293 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