Similar to the question at https://support.bioconductor.org/p/107839/, it appears "refGene" has been removed from the supportedUCSCtracks("mm10") table in GenomicFeatures 1.40.1 and 1.42.1 (if not earlier) but the track is present for mm39. We use mm10's refGene to build our lab's annotation.
It may have been fixed already in the devel version since this code seems to add it manually in makeTxDbFromUCSC.R in lines 253-269 that were added 8 months ago, but I'm not sure how this lines up with those two release versions of GenomicFeatures (https://github.com/Bioconductor/GenomicFeatures/blob/master/R/makeTxDbFromUCSC.R):
if (!(genome %in% c("hg38", "hg19"))) {
## Keep only existing tracks.
ans <- ans[ans$track %in% names(genome_tracknames), , drop=FALSE]
rownames(ans) <- NULL
## Associate subtrack "UCSC RefSeq" to table "refGene" for a few
## genome builds.
if (genome %in% c("mm10", "rn6",
"bosTau9", "danRer10", "danRer11",
"ce11", "dm6", "galGal6", "panTro6",
"rheMac10", "sacCer3"))
{
ans_subtrack <- ans[ , "subtrack"]
ans_subtrack[ans[ , "tablename"] == "refGene"] <- "UCSC RefSeq"
ans[ , "subtrack"] <- ans_subtrack
}
}
Is it possible to restore mm10 refGene to the release version or should we just update to the devel version now? Many thanks!
Session Info:
R version 4.0.4 (2021-02-15)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS Catalina 10.15.6
Matrix products: default BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib LAPACK: /Library/Frameworks/R.framework/Versions/4.0/Resources/lib/libRlapack.dylib
locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages: [1] stats4 parallel grid stats graphics grDevices utils datasets methods base
other attached packages:
[1] ggrepel_0.9.1 scales_1.1.1
[3] broom_0.7.5 readxl_1.3.1
[5] dplyr_1.0.4 plyr_1.8.6
[7] data.table_1.14.0 GenomicAlignments_1.26.0
[9] Rsamtools_2.6.0 BSgenome_1.58.0
[11] rtracklayer_1.50.0 AnnotationHub_2.22.0
[13] BiocFileCache_1.14.0 dbplyr_2.1.0
[15] VennDiagram_1.6.20 futile.logger_1.4.3
[17] RColorBrewer_1.1-2 survival_3.2-7
[19] TCGAbiolinks_2.18.0 RMariaDB_1.1.0
[21] biomaRt_2.46.3 goseq_1.42.0
[23] geneLenDataBase_1.26.0 BiasedUrn_1.07
[25] GO.db_3.12.1 beanplot_1.2
[27] IlluminaHumanMethylation450kmanifest_0.4.0 minfi_1.36.0
[29] bumphunter_1.32.0 locfit_1.5-9.4
[31] iterators_1.0.13 foreach_1.5.1
[33] Biostrings_2.58.0 XVector_0.30.0
[35] SummarizedExperiment_1.20.0 MatrixGenerics_1.2.1
[37] matrixStats_0.58.0 FDb.InfiniumMethylation.hg19_2.2.0
[39] org.Hs.eg.db_3.12.0 TxDb.Hsapiens.UCSC.hg19.knownGene_3.2.2
[41] GenomicFeatures_1.42.1 AnnotationDbi_1.52.0
[43] Biobase_2.50.0 GenomicRanges_1.42.0
[45] GenomeInfoDb_1.26.2 IRanges_2.24.1
[47] S4Vectors_0.28.1 BiocGenerics_0.36.0
[49] seqLogo_1.56.0 gplots_3.1.1
[51] fastcluster_1.1.25 Rcpp_1.0.6
[53] mgcv_1.8-34 nlme_3.1-152
[55] forcats_0.5.1 stringr_1.4.0
[57] purrr_0.3.4 readr_1.4.0
[59] tidyr_1.1.2 tibble_3.0.6
[61] ggplot2_3.3.3 tidyverse_1.3.0
[63] BiocManager_1.30.10
Huh. In a new session
makeTxDbFromUCSC("mm10","refGene")
does work for me. One of the vagaries of R, I guess. Thanks for the quick response!