you're a BOSS! thanks this is exactly what i was looking for
So I'm kind of tired of always using these online conversions that have a limit for how long the input list is . . . .
Is there anywhere where I can download a file (like through UCSC Table Browser or something) to get every single transcript, gene, and gene symbol in mm10
In this format:
ENSMUSTxxxxx [tab] ENSMUGxxxxx [tab] Upf1
ENSMUSTxxxxx [tab] ENSMUGxxxxx [tab] Upf2
ENSMUSTxxxxx [tab] ENSMUGxxxxx [tab] Upf3a
ENSMUSTxxxxx [tab] ENSMUGxxxxx [tab] Upf3b
ENSMUSTxxxxx [tab] ENSMUGxxxxx [tab] Smg1
5 answers
Yes, this is quite easy using UCSC Table Browser or the UCSC public MySQL server.
Using Table Browser, fill in the fields so as they look like this (you may want to enter a file name):
Then, click "get output" and link to the ensemblToGeneName table, so as the fields look like this:
Click "get output" again; here are the first few lines of output:
#mm10.ensGene.name mm10.ensGene.name2 mm10.ensemblToGeneName.value
ENSMUST00000086465 ENSMUSG00000042429 Adora1
ENSMUST00000038191 ENSMUSG00000042429 Adora1
ENSMUST00000169927 ENSMUSG00000042429 Adora1
ENSMUST00000132064 ENSMUSG00000025909 Sntg1
ENSMUST00000140295 ENSMUSG00000025909 Sntg1
ENSMUST00000140302 ENSMUSG00000025909 Sntg1
ENSMUST00000115484 ENSMUSG00000025909 Sntg1
ENSMUST00000135046 ENSMUSG00000025909 Sntg1
ENSMUST00000115488 ENSMUSG00000025909 Sntg1
In case you are comfortable with command line then you can try Neilfws's solution on command line.
mysql --user=genome -N --host=genome-mysql.cse.ucsc.edu -A -D mm10 -e "select name,name2 from ensGene" > Gene1_table
mysql --user=genome -N --host=genome-mysql.cse.ucsc.edu -A -D mm10 -e "select name,value from mm10.ensemblToGeneName" > Gene2_table
paste Gene1_table Gene2_table > mm10_ensembl.txt
Check out the AnnotationHub package in R/Bioconductor. This way you can easily download and access within R all sorts of annotation in just a few lines of code. See the below presentation from the recent CSAMA 15 workshop for some more detail:
These two short YouTube clips are also a good place to start:
Cheers,
Phil
You can do that directly from the Ensembl fasta files, e.g from here. After download, do:
awk '{if ($1 ~ /^>/ ) print}' <(gzcat Homo_sapiens.GRCh38.cdna.all.fa.gz) \
| awk -F " " 'OFS="\t" {print $1, $4, $7}' \
| awk 'OFS="\t" {gsub(">","");gsub("gene:","");gsub("gene_symbol:",""); print}' > outout.tsv
Log in to answer this question.