BED file from GSEA gene set

Question

[2138493o]

Hi all!

Sorry if this is a very simple question but I was wondering what is the easiest way to create a BED file of genes of interest from a GSEA gene set? For context, I would like to use computeMatrix (followed by plotHeatmap) of the deeptools toolkit to plot ChIP-seq signal within the region of genes within a specific GSEA gene set. Thank you in advance everyone!

[ATpoint]

What exactly is the question? Do you need code to make a BED file, and if so based on what? Gene symbols, coordinates of peaks? Please add details.

[2138493o]

I would like to use deeptoolsMatrix to plot ChIP signal within the TSS of my genes of interest; my genes of interest being PRC2 targets. So if I have a set of genes from GSEA (e.g. https://www.gsea-msigdb.org/gsea/msigdb/cards/BENPORATH_PRC2_TARGETS) which provides a list of all the genes and their gene symbols and Entrez ID's, is there anyway I can quickly use this information to create a BED file containing all the genes and their coordinates?

[Marco Pannone]

Download an annotation file (.gff3 file type, for example) and merge the two data frames by column in R, for example performing an inner_join.