Empty gene names from ensembl query using biomart
Here is the code:
require(biomaRt)
ensembl = useMart("ensembl", ,dataset="hsapiens_gene_ensembl")
# filters = listFilters(ensembl)
# filters[1:5, ]
# attributes = listAttributes(ensembl)
# ensembl_attr = cbind(attributes$name, attributes$description)
# head(ensembl_attr, 15)
getBM(c("hgnc_symbol", "start_position", "end_position"), filters=c("chromosome_name", "start", "end"), values=list(16, 1100000, 1250000), mart=ensembl)
Result:
hgnc_symbol start_position end_position
1 1103280 1105461
2 1111627 1113399
3 TPSB2 1227272 1230184
4 1206560 1207124
5 1223639 1224143
6 1156976 1157974
7 1159548 1160176
8 CACNA1H 1153241 1221771
9 TPSG1 1221651 1225257
10 TPSAB1 1240696 1242554
11 1148224 1148754
As you can see, some gene names are missing, why is that?
• 3,444 views
•
link
1 answer
It's because those genes have not got HGNC names. They rather have a clone based name. You may want to include Clone based VEGA gene name and/or Clone based VEGA transcript name in your Attributes.
• 0 views
•
link
Log in to answer this question.