Entering edit mode
Why does Homo.sapiens import TxDb.Hsapiens.UCSC.hg19.knownGene instead of TxDb.Hsapiens.UCSC.hg38.knownGene?
> library(Homo.sapiens)
Loading required package: AnnotationDbi
Loading required package: stats4
Loading required package: BiocGenerics
Loading required package: parallel
Attaching package: ‘BiocGenerics’
The following objects are masked from ‘package:parallel’:
clusterApply, clusterApplyLB, clusterCall, clusterEvalQ, clusterExport, clusterMap, parApply, parCapply, parLapply, parLapplyLB, parRapply, parSapply, parSapplyLB
The following objects are masked from ‘package:stats’:
IQR, mad, xtabs
The following objects are masked from ‘package:base’:
Filter, Find, Map, Position, Reduce, anyDuplicated, append, as.data.frame, as.vector, cbind, colnames, do.call, duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted, lapply, lengths, mapply, match, mget, order, paste, pmax,
pmax.int, pmin, pmin.int, rank, rbind, rownames, sapply, setdiff, sort, table, tapply, union, unique, unlist, unsplit
Loading required package: Biobase
Welcome to Bioconductor
Vignettes contain introductory material; view with 'browseVignettes()'. To cite Bioconductor, see 'citation("Biobase")', and for packages 'citation("pkgname")'.
Loading required package: IRanges
Loading required package: S4Vectors
Loading required package: OrganismDbi
Loading required package: GenomicFeatures
Loading required package: GenomeInfoDb
Loading required package: GenomicRanges
Loading required package: GO.db
Loading required package: DBI
Loading required package: org.Hs.eg.db
Loading required package: TxDb.Hsapiens.UCSC.hg19.knownGene
Now getting the GODb Object directly
Now getting the OrgDb Object directly
Now getting the TxDb Object directly
> sessionInfo()
R version 3.2.2 (2015-08-14)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Red Hat Enterprise Linux Server release 6.7 (Santiago)
locale:
[1] LC_CTYPE=en_GB.UTF-8 LC_NUMERIC=C LC_TIME=en_GB.UTF-8 LC_COLLATE=C LC_MONETARY=en_GB.UTF-8 LC_MESSAGES=C LC_PAPER=en_GB.UTF-8 LC_NAME=C LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_GB.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] parallel stats4 stats graphics grDevices utils datasets methods base
other attached packages:
[1] Homo.sapiens_1.3.1 TxDb.Hsapiens.UCSC.hg19.knownGene_3.2.2 org.Hs.eg.db_3.2.3 GO.db_3.2.2 RSQLite_1.0.0 DBI_0.3.1
[7] OrganismDbi_1.12.1 GenomicFeatures_1.22.13 GenomicRanges_1.22.4 GenomeInfoDb_1.6.3 AnnotationDbi_1.32.3 IRanges_2.4.8
[13] S4Vectors_0.8.11 Biobase_2.30.0 BiocGenerics_0.16.1 BiocInstaller_1.20.1 vimcom_1.2-7.2 colorout_1.0-2
[19] setwidth_1.0-4
loaded via a namespace (and not attached):
[1] Rcpp_0.12.3 futile.logger_1.4.1 XVector_0.10.0 bitops_1.0-6 futile.options_1.0.0 tools_3.2.2 zlibbioc_1.16.0 biomaRt_2.26.1 graph_1.48.0 Rgraphviz_2.14.0
[11] rtracklayer_1.30.3 CellNOptR_1.16.0 Biostrings_2.38.4 grid_3.2.2 XML_3.98-1.4 RBGL_1.46.0 BiocParallel_1.4.3 lambda.r_1.1.7 Rsamtools_1.22.0 GenomicAlignments_1.6.3
[21] SummarizedExperiment_1.0.2 colorspace_1.2-6 RCurl_1.95-4.8
Thank you,

Also note that you can use the
TxDbaccessor to change theTxDbthat you are usingWhich is sort of a pain, but in the end is only two extra lines of code.
Nice, thanks!
Thanks Valerie,
I just have one comment: since the OrganismDb packages do not refer to a specific genome build with their names, I would expect them to use the current one (i.e.: GRCh38). I don't see any reason to use GRCh37/hg19 these days.
Cheers,