Question: Finding Transcriptional Factor Binding Motifs by Table Browser in UCSC
gravatar for mitaminmog
21 months ago by
mitaminmog0 wrote:

I want to find how many motifs bind to the Transcription Factor of ANO1 gene. Would you please help me how to use Table Browser to find Transcriptional Factor binding motifs for gene?

Thank you very much.

tfbs motifs • 897 views
ADD COMMENTlink modified 21 months ago by Alex Reynolds28k • written 21 months ago by mitaminmog0
gravatar for Alex Reynolds
21 months ago by
Alex Reynolds28k
Seattle, WA USA
Alex Reynolds28k wrote:

To my knowledge, the UCSC browser does not offer tables with transcription factor binding sites.

You could make a set of binding sites with FIMO, the JASPAR MEME or other MEME database files available via the MEME download site, and your genome build of interest (per-chromosome FASTA files from UCSC goldenpath, say), e.g. for hg19:

$ for chr in `seq 1 22` X Y; do echo $chr; wget -qO-$chr.fa.gz | gunzip -c - > hg19.chr$chr.fa; done

The output of a FIMO search includes a GFF file you can convert to BED via gff2bed.

Once you have a sorted BED file of binding sites from gff2bed, you can use bedmap --echo-map or similar to locate and count the number of binding sites which overlap annotations that associate with the ANO1 HGNC description. To get HGNC descriptions and their locations, e.g. for hg19:

$ wget -qO- | gunzip -c | awk -v OFS="\t" '($$9>1){ print $$3,$$5,$$6,$$13,$$9,$$4 }' | sort-bed - > hg19.hgnc.bed

Filter for ANO1 regions:

$ awk '$4=="ANO1"' hg19.hgnc.bed > hg19.hgnc.ANO1.bed

Then, for example:

$ bedmap --echo --echo-map-id-uniq hg19.hgnc.ANO1.bed hg19.jaspar_fimo.bed > answer.bed

To count the number of unique JASPAR motif models that overlap the ANO1 space, for example:

$ wc -l answer.bed
ADD COMMENTlink modified 21 months ago • written 21 months ago by Alex Reynolds28k
Please log in to add an answer.


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