[Bioc-devel] Converting gene ids to GRanges - ensembl centric TxDb missing for human
Bhagwat, Aditya
Ad|ty@@Bh@gw@t @end|ng |rom mp|-bn@mpg@de
Tue Oct 15 14:15:50 CEST 2019
Dear BioC devel,
I want to convert geneids to GRanges by doing:
GenomicFeatures::genes(txdb)[geneids]
Works wonderfully for mouse, with entrezgene as well ensemblgene-centric TxDbs:
txdb <- TxDb.Mmusculus.UCSC.mm10.knownGene::TxDb.Mmusculus.UCSC.mm10.knownGene
GenomicFeatures::genes(txdb)[c('100009600', '99889', '99982')]
txdb <- TxDb.Mmusculus.UCSC.mm10.ensGene::TxDb.Mmusculus.UCSC.mm10.ensGene
GenomicFeatures::genes(txdb)[c('ENSMUSG00000000001', 'ENSMUSG00000000003')]
For human, hower, ensembl-centric TxDbs seem to be missing:
txdb <- TxDb.Hsapiens.UCSC.hg38.knownGene::TxDb.Hsapiens.UCSC.hg38.knownGene
GenomicFeatures::genes(txdb)[c('1', '10', '100')]
# No TxDb.Hsapiens.UCSC.hg38.ensGene::TxDb.Hsapiens.UCSC.hg38.ensGene
Has this been a (perhaps recent) design choice to no longer offer the ensemble-centric TxDbs?
(The larger context of this question is the development of multicrispr (https://gitlab.gwdg.de/loosolab/software/multicrispr))
Thankyou for feedback!
Aditya
[[alternative HTML version deleted]]
More information about the Bioc-devel
mailing list