RNAmodR: creating classes for additional modification detection from high throughput sequencing. (original) (raw)
Introduction
For users interested in the general aspect of any RNAmodR based package please have a look at the main vignette of the package.
This vignette is aimed at developers and researchers, who want to use the functionality of the RNAmodR package to develop a new modification strategy based on high throughput sequencing data.
library(RNAmodR)Two classes have to be considered to establish a new analysis pipeline usingRNAmodR. These are the SequenceData and the Modifier class.
A new SequenceData class
First, the SequenceData class has to be considered. Several classes are already implemented, which are:
End5SequenceDataEnd3SequenceDataEndSequenceDataProtectedEndSequenceDataCoverageSequenceDataPileupSequenceDataNormEnd5SequenceDataNormEnd3SequenceData
If these cannot be reused, a new class can be implemented quite easily. First the DataFrame class, the Data class and a constructor has to defined. The only value, which has to be provided, is a default minQuality integer value and some basic information.
setClass(Class = "ExampleSequenceDataFrame",
contains = "SequenceDFrame")
ExampleSequenceDataFrame <- function(df, ranges, sequence, replicate,
condition, bamfiles, seqinfo){
RNAmodR:::.SequenceDataFrame("Example",df, ranges, sequence, replicate,
condition, bamfiles, seqinfo)
}
setClass(Class = "ExampleSequenceData",
contains = "SequenceData",
slots = c(unlistData = "ExampleSequenceDataFrame"),
prototype = list(unlistData = ExampleSequenceDataFrame(),
unlistType = "ExampleSequenceDataFrame",
minQuality = 5L,
dataDescription = "Example data"))
ExampleSequenceData <- function(bamfiles, annotation, sequences, seqinfo, ...){
RNAmodR:::SequenceData("Example", bamfiles = bamfiles,
annotation = annotation, sequences = sequences,
seqinfo = seqinfo, ...)
}Second, the getData function has to be implemented. This is used to load the data from a bam file and must return a named list IntegerList,NumericList or CompressedSplitDataFrameList per file.
setMethod("getData",
signature = c(x = "ExampleSequenceData",
bamfiles = "BamFileList",
grl = "GRangesList",
sequences = "XStringSet",
param = "ScanBamParam"),
definition = function(x, bamfiles, grl, sequences, param, args){
###
}
)Third, the aggregate function has to be implemented. This function is used to aggregate data over replicates for all or one of the conditions. The resulting data is passed on to the Modifier class.
setMethod("aggregateData",
signature = c(x = "ExampleSequenceData"),
function(x, condition = c("Both","Treated","Control")){
###
}
)A new Modifier class
A new Modifier class is probably the main class, which needs to be implemented. Three variable have to be set. mod must be a single element from the Modstrings::shortName(Modstrings::ModRNAString()). score is the default score, which is used for several function. A column with this name should be returned from the aggregate function. dataType defines the SequenceDataclass to be used. dataType can contain multiple names of a SequenceDataclass, which are then combined to form a SequenceDataSet.
setClass("ModExample",
contains = c("RNAModifier"),
prototype = list(mod = "X",
score = "score",
dataType = "ExampleSequenceData"))
ModExample <- function(x, annotation, sequences, seqinfo, ...){
RNAmodR:::Modifier("ModExample", x = x, annotation = annotation,
sequences = sequences, seqinfo = seqinfo, ...)
}dataType can also be a list of character vectors, which leads then to the creation of SequenceDataList. However, for now this is a hypothetical case and should only be used, if the detection of a modification requires bam files from two or more different methods to be used to detect one modification.
The settings<- function can be amended to save specifc settings (.norm_example_args must be defined seperatly to normalize input arguments in any way one sees fit).
setReplaceMethod(f = "settings",
signature = signature(x = "ModExample"),
definition = function(x, value){
x <- callNextMethod()
# validate special setting here
x@settings[names(value)] <- unname(.norm_example_args(value))
x
})The aggregateData function is used to take the aggregated data from theSequenceData object and to calculate the specific scores, which are then stored in the aggregate slot.
setMethod(f = "aggregateData",
signature = signature(x = "ModExample"),
definition =
function(x, force = FALSE){
# Some data with element per transcript
}
)The findMod function takes the aggregate data and searches for modifications, which are then returned as a GRanges object and stored in the modificationsslot.
setMethod("findMod",
signature = c(x = "ModExample"),
function(x){
# an element per modification found.
}
)A new ModifierSet class
The ModifierSet class is implemented very easily by defining the class and the constructor. The functionality is defined by the Modifier class.
setClass("ModSetExample",
contains = "ModifierSet",
prototype = list(elementType = "ModExample"))
ModSetExample <- function(x, annotation, sequences, seqinfo, ...){
RNAmodR:::ModifierSet("ModExample", x = x, annotation = annotation,
sequences = sequences, seqinfo = seqinfo, ...)
}Visualization functions
Additional functions, which need to be implemented, are getDataTrack for the new SequenceData and new Modifier classes andplotData/plotDataByCoord for the new Modifier and ModifierSetclasses. name defines a transcript name found in names(ranges(x)) andtype is the data type typically found as a column in the aggregate slot.
setMethod(
f = "getDataTrack",
signature = signature(x = "ExampleSequenceData"),
definition = function(x, name, ...) {
###
}
)
setMethod(
f = "getDataTrack",
signature = signature(x = "ModExample"),
definition = function(x, name, type, ...) {
}
)
setMethod(
f = "plotDataByCoord",
signature = signature(x = "ModExample", coord = "GRanges"),
definition = function(x, coord, type = "score", window.size = 15L, ...) {
}
)
setMethod(
f = "plotData",
signature = signature(x = "ModExample"),
definition = function(x, name, from, to, type = "score", ...) {
}
)
setMethod(
f = "plotDataByCoord",
signature = signature(x = "ModSetExample", coord = "GRanges"),
definition = function(x, coord, type = "score", window.size = 15L, ...) {
}
)
setMethod(
f = "plotData",
signature = signature(x = "ModSetExample"),
definition = function(x, name, from, to, type = "score", ...) {
}
)If unsure, how to modify these functions, have a look a the code in theModifier-Inosine-viz.R file of this package.
Summary
As suggested directly above, for a more detailed example have a look at theModInosine class source code found in the Modifier-Inosine-class.R andModifier-Inosine-viz.R files of this package.
Sessioninfo
sessionInfo()## R version 4.5.0 RC (2025-04-04 r88126)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 24.04.2 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.21-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.12.0 LAPACK version 3.12.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] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] RNAmodR_1.22.0 Modstrings_1.24.0 RNAmodR.Data_1.21.0
## [4] ExperimentHubData_1.34.0 AnnotationHubData_1.38.0 futile.logger_1.4.3
## [7] ExperimentHub_2.16.0 AnnotationHub_3.16.0 BiocFileCache_2.16.0
## [10] dbplyr_2.5.0 txdbmaker_1.4.0 GenomicFeatures_1.60.0
## [13] AnnotationDbi_1.70.0 Biobase_2.68.0 Rsamtools_2.24.0
## [16] Biostrings_2.76.0 XVector_0.48.0 rtracklayer_1.68.0
## [19] GenomicRanges_1.60.0 GenomeInfoDb_1.44.0 IRanges_2.42.0
## [22] S4Vectors_0.46.0 BiocGenerics_0.54.0 generics_0.1.3
## [25] BiocStyle_2.36.0
##
## loaded via a namespace (and not attached):
## [1] RColorBrewer_1.1-3 rstudioapi_0.17.1
## [3] jsonlite_2.0.0 magrittr_2.0.3
## [5] magick_2.8.6 farver_2.1.2
## [7] rmarkdown_2.29 BiocIO_1.18.0
## [9] vctrs_0.6.5 ROCR_1.0-11
## [11] memoise_2.0.1 RCurl_1.98-1.17
## [13] base64enc_0.1-3 tinytex_0.57
## [15] htmltools_0.5.8.1 S4Arrays_1.8.0
## [17] BiocBaseUtils_1.10.0 progress_1.2.3
## [19] lambda.r_1.2.4 curl_6.2.2
## [21] SparseArray_1.8.0 Formula_1.2-5
## [23] sass_0.4.10 bslib_0.9.0
## [25] htmlwidgets_1.6.4 plyr_1.8.9
## [27] Gviz_1.52.0 httr2_1.1.2
## [29] futile.options_1.0.1 cachem_1.1.0
## [31] GenomicAlignments_1.44.0 mime_0.13
## [33] lifecycle_1.0.4 pkgconfig_2.0.3
## [35] Matrix_1.7-3 R6_2.6.1
## [37] fastmap_1.2.0 GenomeInfoDbData_1.2.14
## [39] BiocCheck_1.44.0 MatrixGenerics_1.20.0
## [41] digest_0.6.37 colorspace_2.1-1
## [43] OrganismDbi_1.50.0 Hmisc_5.2-3
## [45] RSQLite_2.3.9 labeling_0.4.3
## [47] filelock_1.0.3 colorRamps_2.3.4
## [49] httr_1.4.7 abind_1.4-8
## [51] compiler_4.5.0 withr_3.0.2
## [53] bit64_4.6.0-1 backports_1.5.0
## [55] htmlTable_2.4.3 biocViews_1.76.0
## [57] BiocParallel_1.42.0 DBI_1.2.3
## [59] biomaRt_2.64.0 rappdirs_0.3.3
## [61] DelayedArray_0.34.0 rjson_0.2.23
## [63] tools_4.5.0 foreign_0.8-90
## [65] nnet_7.3-20 glue_1.8.0
## [67] restfulr_0.0.15 grid_4.5.0
## [69] stringdist_0.9.15 checkmate_2.3.2
## [71] reshape2_1.4.4 cluster_2.1.8.1
## [73] gtable_0.3.6 BSgenome_1.76.0
## [75] ensembldb_2.32.0 data.table_1.17.0
## [77] hms_1.1.3 xml2_1.3.8
## [79] BiocVersion_3.21.1 pillar_1.10.2
## [81] stringr_1.5.1 dplyr_1.1.4
## [83] lattice_0.22-7 deldir_2.0-4
## [85] bit_4.6.0 biovizBase_1.56.0
## [87] tidyselect_1.2.1 RBGL_1.84.0
## [89] knitr_1.50 gridExtra_2.3
## [91] bookdown_0.43 ProtGenerics_1.40.0
## [93] SummarizedExperiment_1.38.0 xfun_0.52
## [95] matrixStats_1.5.0 stringi_1.8.7
## [97] UCSC.utils_1.4.0 lazyeval_0.2.2
## [99] yaml_2.3.10 evaluate_1.0.3
## [101] codetools_0.2-20 interp_1.1-6
## [103] tibble_3.2.1 BiocManager_1.30.25
## [105] graph_1.86.0 cli_3.6.4
## [107] rpart_4.1.24 munsell_0.5.1
## [109] jquerylib_0.1.4 Rcpp_1.0.14
## [111] dichromat_2.0-0.1 png_0.1-8
## [113] XML_3.99-0.18 RUnit_0.4.33
## [115] parallel_4.5.0 ggplot2_3.5.2
## [117] blob_1.2.4 prettyunits_1.2.0
## [119] jpeg_0.1-11 latticeExtra_0.6-30
## [121] AnnotationFilter_1.32.0 AnnotationForge_1.50.0
## [123] bitops_1.0-9 VariantAnnotation_1.54.0
## [125] scales_1.3.0 purrr_1.0.4
## [127] crayon_1.5.3 rlang_1.1.6
## [129] KEGGREST_1.48.0 formatR_1.14