Hi all,
When I tried to retrieve attributes "description" from the ensembl database for Drosophila, it returned NA for all the genes.
dmGeneList <- c("FBgn0000003 ", "FBgn0000008", "FBgn0000014")
ensembl = useMart("ensembl",dataset="dmelanogaster_gene_ensembl")
annoTable <- getBM(
attributes=c("description", "external_gene_name", "flybase_gene_id"),
filters = 'flybase_gene_id',
values = dmGeneList,
mart = ensembl)
description external_gene_name flybase_gene_id
1 NA 7SLRNA:CR32864 FBgn0000003
2 NA a FBgn0000008
3 NA abd-A FBgn0000014
I searched the attributes with the name "description",
searchAttributes(ensembl, pattern = "^description")
It turns out that there are attributes with the same name "description" on different pages.
name description page
5 description Gene description feature_page
137 description Gene description structure
164 description Gene description homologs
2471 description Gene description sequences
I guess this is the reason for the "NA". How can I retrieve the attributes from only the "feature_page"?
Thanks!
Wanpeng
Thank you, James! It's a great idea to do the search first! The "gene description" column seems to be empty in this version. It worked previously. I'm turning to flybase and use the "gene snapshot" table. This annotation package for fly is also very useful. Thank you for sharing!
Wanpeng