Thank you for your help!
Hello,
I'm new here and and I need help please!
I need to download a list of all human genes with their respective gene symbol | chromosome | strand | transcription start site | Txen | and Ensembl gene name,
Actually I'm using UCSC table and I get something like this:
#hg38.knownGene.name hg38.knownGene.chrom hg38.knownGene.strand hg38.knownGene.txStart hg38.knownGene.txEnd hg38.kgXref.geneSymbol
uc001aaa.3 chr1 + 11873 14409 DDX11L1
And I want to know if it's possible to include ensembl gene symbol with UCSC table or with another method
Thank you in advance
Cherif
4 answers
Unfortunately UCSC table browser doesn't have Ensembl gene track for hg38. But they do have it for hg19 and the below command should work for you.
mysql \
--user=genome \
-N \
--host=genome-mysql.cse.ucsc.edu \
-A \
-D hg19 \
-e "select ensGene.name, name2, chrom, strand, txStart, txEnd, value from ensGene, ensemblToGeneName where ensGene.name = ensemblToGeneName.name" > \
output.txt
You can try the same command with hg38 but you will have to choose other gene models such as refseq or ucsc.
Thank you very much!!! Yep it works!!
Follow these steps: http://kandurilab.org/bioinformatics/biostars/UCSC_tableBrowser_annotation.pdf
Here's an easy way to do it from UCSC's table browser:
- In the table browser, select Ensembl Genes as your track
- Under output format, choose "selected fields from primary and related tables"
- Add your output file name, if you want one (otherwise, it will print to the browser)
- Click "get output"
- On the next page, you will get to choose your fields.
- Under linked tables, check
ensemblToGeneName, then press "Allow selection from checked tables" - The page will refresh, and you should have a new table called
hg19.ensemblToGeneName - Check of name, chrom, strand, txStart, txEnd, and name2 in hg19.ensGene (or any fields you'd like)
- In
hg19.ensemblToGeneName, check "value", which has the description "alternate gene name" - Press "get output"
If you did it right, you should get a table that looks a bit like this (I took this chunk from chr1:100,000,000-150,000,000):
#hg19.ensGene.name hg19.ensGene.chrom hg19.ensGene.strand hg19.ensGene.txStart hg19.ensGene.txEnd hg19.ensGene.name2 hg19.ensemblToGeneName.value
ENST00000263174 chr1 + 100111498 100160097 ENSG00000099260 PALMD
ENST00000605497 chr1 + 100111748 100155633 ENSG00000099260 PALMD
ENST00000605613 chr1 + 100133135 100135379 ENSG00000099260 PALMD
ENST00000496843 chr1 + 100148821 100160097 ENSG00000099260 PALMD
ENST00000434734 chr1 + 100163797 100164734 ENSG00000223656 HMGB3P10
Ashutosh noticed what I missed -- hg38 does not have Ensembl annotations in UCSC yet. Is there a reason you're choosing hg38 and not hg19?
Actually I can use Hg19, but maybe I will need it soon (with Hg38)
ps: Thanks a lot! Yes I get it.
Just one more question please
When select Ensembl Genes as my track , I get "204941" genes and their respective txStart and ..
and when I select ucsc genes as my track , I get just 82961 => about the half.
Is it normal?
Those are two different genome annotations--
The UCSC gene track is described here. It is a set of genes taken from RefSeq, GenBank, CCDS, Rfam, and the tRNA genes track.
I couldn't find a similarly clear description for Ensembl, but this is a good start. It seems they rely on deposited mRNAs and protein sequences in public databases. That might mean that their curation is a bit more relaxed than RefSeq, CCDS, etc.
FWIW, whenever I do annotation, I've generally relied on CCDS and RefSeq.
Log in to answer this question.
Hi! Do you tried Biomart?
I'm about to try (until now I do not know how proceed :/)