Get gene promoter coordinates
0
1
Entering edit mode
9 months ago
sabin ▴ 40

Hello everyone, I'm trying to find a way to get the coordinates of the promoter region of a gene (if possible, of a list of genes). So let's suppose that i have the gene BRCA2 and I know that for genome version hg38 this gene has coordinates Chr13: 32,315,086-32,400,268. Now what i would like to get is the the start and end coordinate of the promoter region. Is there a way to get such info? I've already tried with biomart-ensembl (but i managed to get only the promoter sequence and not the coordinates) and plus i can do it only one gene at a time and, since i have a long list, that would take too long. Moreover, I tried in R with the function promoter() but it doesn't return results for all the genes in the list, but only for some of them.

Do you know any way to get these coordinates? Any additional info is very appreciated! Thank you

gene promoter coordinates hg38 • 558 views
0
Entering edit mode

The promoter is right upstream of the gene start coordinate (TSS to be precise), people often use something like 200-500bp upstream as a proxy, isn't that sufficient for your purpose?

0
Entering edit mode

Hi! Yes either 500bp or 2k bp are ok. Do you know how to get the coordinates?

0
Entering edit mode

From which package is the promoter() function? The FANTOM project surveyed enhancer and promoter regions genome-wide - perhaps you can use their data: https://fantom.gsc.riken.jp/5/

1
Entering edit mode

Hi Kevin, the function prom promoter() is from the package GenomicFeatures.