Abstract
NxtIRFdata is a data package containing ready-made BED files of Mappability exclusion genomic regions. It also contains a fully-functioning example data set with a “mock” genome and genome annotation to demonstrate the functionalities of SpliceWiz, a powerful and interactive analysis and visualisation tool for alternative splicing.
NxtIRFcore’s full functionality (plus more) will be replaced by SpliceWiz in Bioconductor version 3.16 onwards!
To install this package, start R (version “4.2.1”) and enter:
if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("NxtIRFdata")Start using NxtIRFdata:
Examples in SpliceWiz are demonstrated using an artificial genome and gene annotation. A synthetic reference, with genome sequence (FASTA) and gene annotation (GTF) files are provided, based on the genes SRSF1, SRSF2, SRSF3, TRA2A, TRA2B, TP53 and NSUN5. These genes, each with an additional 100 flanking nucleotides, were used to construct an artificial “chromosome Z” (chrZ). Gene annotations, based on release-94 of Ensembl GRCh38 (hg38), were modified with genome coordinates corresponding to this artificial chromosome.
These files can be accessed as follows:
The set of 6 BAM files used in the SpliceWiz vignette / example code can be downloaded to a path of the user’s choice using the following function:
Note that this downloads BAM files and not their respective BAI (BAM file indices). This is because SpliceWiz reads BAM files natively and does not require RSamtools. BAI files are provided with BAM files in their respective ExperimentHub entries for users wishing to view these files using RSamtools.
NxtIRFdata retrieves the relevant records from AnnotationHub and makes a local copy of the BED file. This BED file is used to produce Mappability Exclusion information to SpliceWiz.
Note that this function is intended to be called internally by SpliceWiz. Users interested in the format or nature of the Mappability BED file can call this function to examine the contents of the BED file
# To get the MappabilityExclusion for hg38 as a GRanges object
gr = get_mappability_exclusion(genome_type = "hg38", as_type = "GRanges")
# To get the MappabilityExclusion for hg38 as a locally-copied gzipped BED file
bed_path = get_mappability_exclusion(genome_type = "hg38", as_type = "bed.gz",
path = tempdir())
# Other `genome_type` values include "hg19", "mm10", and "mm9"The data deposited in ExperimentHub can be accessed as follows:
For more information about the example BAM files, refer to the NxtIRFdata package documentation:
## R version 4.3.1 (2023-06-16)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.18-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] ExperimentHub_2.10.0 AnnotationHub_3.10.0 BiocFileCache_2.10.0
## [4] dbplyr_2.3.4 BiocGenerics_0.48.0 S4Vectors_0.40.0
## [7] IRanges_2.36.0 NxtIRFdata_1.8.0
##
## loaded via a namespace (and not attached):
## [1] tidyselect_1.2.0 dplyr_1.1.3
## [3] blob_1.2.4 filelock_1.0.2
## [5] R.utils_2.12.2 Biostrings_2.70.1
## [7] bitops_1.0-7 fastmap_1.1.1
## [9] RCurl_1.98-1.12 promises_1.2.1
## [11] GenomicAlignments_1.38.0 XML_3.99-0.14
## [13] digest_0.6.33 mime_0.12
## [15] lifecycle_1.0.3 ellipsis_0.3.2
## [17] KEGGREST_1.42.0 interactiveDisplayBase_1.40.0
## [19] RSQLite_2.3.1 magrittr_2.0.3
## [21] compiler_4.3.1 rlang_1.1.1
## [23] sass_0.4.7 tools_4.3.1
## [25] utf8_1.2.4 yaml_2.3.7
## [27] rtracklayer_1.62.0 knitr_1.44
## [29] S4Arrays_1.2.0 bit_4.0.5
## [31] curl_5.1.0 DelayedArray_0.28.0
## [33] abind_1.4-5 BiocParallel_1.36.0
## [35] withr_2.5.1 purrr_1.0.2
## [37] R.oo_1.25.0 grid_4.3.1
## [39] stats4_4.3.1 fansi_1.0.5
## [41] xtable_1.8-4 SummarizedExperiment_1.32.0
## [43] cli_3.6.1 rmarkdown_2.25
## [45] crayon_1.5.2 generics_0.1.3
## [47] httr_1.4.7 rjson_0.2.21
## [49] DBI_1.1.3 cachem_1.0.8
## [51] zlibbioc_1.48.0 parallel_4.3.1
## [53] AnnotationDbi_1.64.0 BiocManager_1.30.22
## [55] XVector_0.42.0 restfulr_0.0.15
## [57] matrixStats_1.0.0 vctrs_0.6.4
## [59] Matrix_1.6-1.1 jsonlite_1.8.7
## [61] bit64_4.0.5 jquerylib_0.1.4
## [63] glue_1.6.2 codetools_0.2-19
## [65] later_1.3.1 BiocVersion_3.18.0
## [67] GenomeInfoDb_1.38.0 BiocIO_1.12.0
## [69] GenomicRanges_1.54.0 tibble_3.2.1
## [71] pillar_1.9.0 rappdirs_0.3.3
## [73] htmltools_0.5.6.1 GenomeInfoDbData_1.2.11
## [75] R6_2.5.1 evaluate_0.22
## [77] shiny_1.7.5.1 Biobase_2.62.0
## [79] lattice_0.22-5 png_0.1-8
## [81] R.methodsS3_1.8.2 Rsamtools_2.18.0
## [83] memoise_2.0.1 httpuv_1.6.12
## [85] bslib_0.5.1 Rcpp_1.0.11
## [87] SparseArray_1.2.0 xfun_0.40
## [89] MatrixGenerics_1.14.0 pkgconfig_2.0.3