0 Important announcement

NxtIRFcore’s full functionality (plus more) will be replaced by SpliceWiz in Bioconductor version 3.16 onwards!

1 Installation

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:

library(NxtIRFdata)

2 Obtaining the example NxtIRF genome and gene annotation files

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:

example_fasta = chrZ_genome()
example_gtf = chrZ_gtf()

3 Obtaining the example BAM file dataset for SpliceWiz

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:

bam_paths = example_bams(path = tempdir())

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.

4 Obtaining the Mappability Exclusion BED files

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"

5 Accessing NxtIRFdata via ExperimentHub

The data deposited in ExperimentHub can be accessed as follows:

library(ExperimentHub)
eh = ExperimentHub()
NxtIRF_hub = query(eh, "NxtIRF")

NxtIRF_hub

temp = eh[["EH6792"]]
temp

temp = eh[["EH6787"]]
temp

5 Information about the example BAM files

For more information about the example BAM files, refer to the NxtIRFdata package documentation:

?`NxtIRFdata-package`

SessionInfo

sessionInfo()
## 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