EpiTxDb 0.99.5
To create an EpiTxDb object a number of different functions are available.
The most univeral functions are makeEpiTxDb and makeEpiTxDbFromGRanges.
makeEpiTxDb uses four data.frames as input, whereas makeEpiTxDbFromGRanges
is a wrapper for information available as a GRanges object.
The other functions are makeEpiTxDbFromRMBase and makeEpiTxDbFromtRNAdb,
which are aimed to make data available from the RMBase v2.0 database
(Xuan et al. 2017; Sun et al. 2015) or the tRNAdb
(Jühling et al. 2009; Sprinzl and Vassilenko 2005). However, before creating your
EpiTxDb objects, have a look at the already available resources for
H. sapiens, M. musculus and S. cerevisiae.
Additional metadata can be provided as separate data.frame for all functions.
The data.frame must have two columns name and value.
library(GenomicRanges)
library(EpiTxDb)makeEpiTxDb and makeEpiTxDbFromGRangesThe creation of an etdb object is quite easy starting with a GRanges object.
gr <- GRanges(seqnames = "test",
ranges = IRanges::IRanges(1,1),
strand = "+",
DataFrame(mod_id = 1L,
mod_type = "Am",
mod_name = "Am_1"))
etdb <- makeEpiTxDbFromGRanges(gr, metadata = data.frame(name = "test",
value = "Yes"))## Creating EpiTxDb object ... doneetdb## EpiTxDb object:
## # Db type: EpiTxDb
## # Supporting package: EpiTxDb
## # test: Yes
## # Nb of modifications: 1
## # Db created by: EpiTxDb package from Bioconductor
## # Creation time: 2020-04-08 03:42:56 -0400 (Wed, 08 Apr 2020)
## # EpiTxDb version at creation time: 0.99.5
## # RSQLite version at creation time: 2.2.0
## # DBSCHEMAVERSION: 1.0metadata(etdb)Additional data can be provided via the metadata columns of the GRanges
object. For supported columns have a look at ?makeEpiTxDb or
?makeEpiTxDbFromGRanges.
makeEpiTxDbFromtRNAdbThe information of the tRNAdb can be accessed via the tRNAdbImport package
using the RNA database. As a result a ModRNAStringSet object is returned from
which the modifications can be extracted using separate().
The only input require is a valid organism name returned by
listAvailableOrganismsFromtRNAdb().
etdb <- makeEpiTxDbFromtRNAdb("Saccharomyces cerevisiae")## Loading data from tRNAdb ...## Assembling data ...## Creating EpiTxDb object ... doneetdb## EpiTxDb object:
## # Db type: EpiTxDb
## # Supporting package: EpiTxDb
## # Nb of modifications: 557
## # Db created by: EpiTxDb package from Bioconductor
## # Creation time: 2020-04-08 03:43:25 -0400 (Wed, 08 Apr 2020)
## # EpiTxDb version at creation time: 0.99.5
## # RSQLite version at creation time: 2.2.0
## # DBSCHEMAVERSION: 1.0For additional information have a look at ?makeEpiTxDbFromtRNAdb. The result
returned from the tRNAdb is also available as GRanges object, if
gettRNAdbDataAsGRanges() is used.
makeEpiTxDbFromRMBaseAnalogous to the example above makeEpiTxDbFromRMBase() will download the data
from the RMBase v2.0. Three inputs are required, organism, genome and
modtype, which have to valid bia the functions
listAvailableOrganismsFromRMBase(), .listAvailableGenomesFromRMBase() and
listAvailableModFromRMBase.
etdb <- makeEpiTxDbFromRMBase(organism = "Saccharomyces cerevisiae",
genome = "sacCer3",
modtype = "m1A")Internally, the files are cached using the BiocFileCache package and passed
to makeEpiTxDbFromRMBaseFiles(), which can also be used with locally stored
files. The resuls for creating the EpiTxDb class are processed from these
files via the getRMBaseDataAsGRanges() function.
sessionInfo()## R version 4.0.0 alpha (2020-03-31 r78116)
## Platform: x86_64-w64-mingw32/x64 (64-bit)
## Running under: Windows Server 2012 R2 x64 (build 9600)
##
## Matrix products: default
##
## locale:
## [1] LC_COLLATE=C
## [2] LC_CTYPE=English_United States.1252
## [3] LC_MONETARY=English_United States.1252
## [4] LC_NUMERIC=C
## [5] LC_TIME=English_United States.1252
##
## attached base packages:
## [1] parallel stats4 stats graphics grDevices utils datasets
## [8] methods base
##
## other attached packages:
## [1] EpiTxDb_0.99.5 Modstrings_1.3.12 Biostrings_2.55.7
## [4] XVector_0.27.2 AnnotationDbi_1.49.1 Biobase_2.47.3
## [7] GenomicRanges_1.39.3 GenomeInfoDb_1.23.16 IRanges_2.21.8
## [10] S4Vectors_0.25.15 BiocGenerics_0.33.3 BiocStyle_2.15.6
##
## loaded via a namespace (and not attached):
## [1] bitops_1.0-6 matrixStats_0.56.0
## [3] assertive.models_0.0-2 bit64_0.9-7
## [5] progress_1.2.2 httr_1.4.1
## [7] assertive.datetimes_0.0-2 tools_4.0.0
## [9] R6_2.4.1 colorspace_1.4-1
## [11] DBI_1.1.0 assertive.data_0.0-3
## [13] assertive.reflection_0.0-4 tidyselect_1.0.0
## [15] prettyunits_1.1.1 bit_1.1-15.2
## [17] curl_4.3 compiler_4.0.0
## [19] cli_2.0.2 assertive.properties_0.0-4
## [21] xml2_1.3.0 DelayedArray_0.13.10
## [23] rtracklayer_1.47.0 assertive.files_0.0-2
## [25] bookdown_0.18 scales_1.1.0
## [27] tRNA_1.5.2 askpass_1.1
## [29] rappdirs_0.3.1 stringr_1.4.0
## [31] digest_0.6.25 Rsamtools_2.3.7
## [33] rmarkdown_2.1 assertive.numbers_0.0-2
## [35] pkgconfig_2.0.3 htmltools_0.4.0
## [37] dbplyr_1.4.2 rlang_0.4.5
## [39] RSQLite_2.2.0 assertive_0.3-5
## [41] jsonlite_1.6.1 BiocParallel_1.21.2
## [43] dplyr_0.8.5 RCurl_1.98-1.1
## [45] magrittr_1.5 GenomeInfoDbData_1.2.2
## [47] Matrix_1.2-18 munsell_0.5.0
## [49] Rcpp_1.0.4 fansi_0.4.1
## [51] lifecycle_0.2.0 stringi_1.4.6
## [53] assertive.base_0.0-7 yaml_2.2.1
## [55] SummarizedExperiment_1.17.5 zlibbioc_1.33.1
## [57] BiocFileCache_1.11.4 grid_4.0.0
## [59] blob_1.2.1 crayon_1.3.4
## [61] lattice_0.20-41 tRNAdbImport_1.5.6
## [63] assertive.code_0.0-3 GenomicFeatures_1.39.7
## [65] hms_0.5.3 knitr_1.28
## [67] pillar_1.4.3 assertive.sets_0.0-3
## [69] codetools_0.2-16 biomaRt_2.43.5
## [71] XML_3.99-0.3 glue_1.4.0
## [73] evaluate_0.14 Structstrings_1.3.5
## [75] BiocManager_1.30.10 vctrs_0.2.4
## [77] gtable_0.3.0 openssl_1.4.1
## [79] purrr_0.3.3 assertive.strings_0.0-3
## [81] assertthat_0.2.1 ggplot2_3.3.0
## [83] xfun_0.12 assertive.types_0.0-3
## [85] assertive.data.uk_0.0-2 tibble_3.0.0
## [87] GenomicAlignments_1.23.2 memoise_1.1.0
## [89] assertive.matrices_0.0-2 ellipsis_0.3.0
## [91] assertive.data.us_0.0-2Jühling, Frank, Mario Mörl, Roland K. Hartmann, Mathias Sprinzl, Peter F. Stadler, and Joern Pütz. 2009. “TRNAdb 2009: Compilation of tRNA Sequences and tRNA Genes.” Nucleic Acids Research 37: D159–D162. https://doi.org/10.1093/nar/gkn772.
Sprinzl, Mathias, and Konstantin S. Vassilenko. 2005. “Compilation of tRNA Sequences and Sequences of tRNA Genes.” Nucleic Acids Research 33: D139–D140. https://doi.org/10.1093/nar/gki012.
Sun, Wen-Ju, Jun-Hao Li, Shun Liu, Jie Wu, Hui Zhou, Liang-Hu Qu, and Jian-Hua Yang. 2015. “RMBase: a resource for decoding the landscape of RNA modifications from high-throughput sequencing data.” Nucleic Acids Research 44 (D1): D259–D265. https://doi.org/10.1093/nar/gkv1036.
Xuan, Jia-Jia, Wen-Ju Sun, Peng-Hui Lin, Ke-Ren Zhou, Shun Liu, Ling-Ling Zheng, Liang-Hu Qu, and Jian-Hua Yang. 2017. “RMBase v2.0: deciphering the map of RNA modifications from epitranscriptome sequencing data.” Nucleic Acids Research 46 (D1): D327–D334. https://doi.org/10.1093/nar/gkx934.