Prepare toy data

First, we import the sample data.

library(trackViewer)
extdata <- system.file("extdata", package="trackViewer",
                       mustWork=TRUE)
repA <- importScore(file.path(extdata, "cpsf160.repA_-.wig"),
                    file.path(extdata, "cpsf160.repA_+.wig"),
                    format="WIG")
## Because the wig file does not contain any strand info, 
## we need to set it manually.
strand(repA$dat) <- "-"
strand(repA$dat2) <- "+"

fox2 <- importScore(file.path(extdata, "fox2.bed"), format="BED",
                    ranges=GRanges("chr11", IRanges(122830799, 123116707)))
dat <- coverageGR(fox2$dat)
## We can split the data by strand into two different track channels
## Here, we set the dat2 slot to save the negative strand info. 
 
fox2$dat <- dat[strand(dat)=="+"]
fox2$dat2 <- dat[strand(dat)=="-"]

library(TxDb.Hsapiens.UCSC.hg19.knownGene)
library(org.Hs.eg.db)
gr <- GRanges("chr11", IRanges(122929275, 122930122), strand="-")
trs <- geneModelFromTxdb(TxDb.Hsapiens.UCSC.hg19.knownGene,
                         org.Hs.eg.db,
                         gr=gr)
optSty <- optimizeStyle(trackList(repA, fox2, trs))
trackList <- optSty$tracks
viewerStyle <- optSty$style

Using the browseTracks function as a helper

Since version 1.33.6, the interactive plot generated by browseTracks function will print the psuedocode about how to change the track styles at the bottom, which can be used as a quick reference.

browseTracks(trackList, gr=gr)



Adjust the x-axis or the X scale

In most cases, researchers are interested in the relative position of the peaks in the gene. Sometimes, margin needs to be adjusted to be able to show the entire gene model. The Figure below shows how to add an X scale (x-scale) and remove the x-axis using the setTrackXscaleParam and setTrackViewerStyleParam functions.

setTrackViewerStyleParam(viewerStyle, "xaxis", FALSE)
setTrackViewerStyleParam(viewerStyle, "margin", c(.01, .05, .01, .01))
setTrackXscaleParam(trackList[[1]], "draw", TRUE)
setTrackXscaleParam(trackList[[1]], "gp", list(cex=0.8))
viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
plot data with x-scale

plot data with x-scale

setTrackXscaleParam(trackList[[1]], attr="position", 
                    value=list(x=122929700, y=3, label=200))
viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
plot data with x-scale

plot data with x-scale

Adjust the y-axis

The y-axis can be put to the right side of the track by setting the main slot to FALSE in the y-axis slot of each track. In addition, the limit of y-axis (ylim) can be set by setTrackStyleParam.

setTrackViewerStyleParam(viewerStyle, "margin", c(.01, .05, .01, .05))
for(i in 1:2){
    setTrackYaxisParam(trackList[[i]], "main", FALSE)
}
## Adjust the limit of y-axis
setTrackStyleParam(trackList[[1]], "ylim", c(0, 25))
setTrackStyleParam(trackList[[2]], "ylim", c(-25, 0))
viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
plot data with y-axis in right side

plot data with y-axis in right side

The function setTrackYaxisParam can be sued to adjust the track color/size/etc of y-axis and the y-axis labels.

## change the 
setTrackYaxisParam(trackList[[1]], "gp", list(cex=.8, col="green"))
setTrackYaxisParam(trackList[[2]], "gp", list(cex=.8, col="blue"))
viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
plot data with adjusted track color

plot data with adjusted track color

Adjust the label of y-axis

The style of y-axis can be changed by setting the ylabgp slot in the style of each track.

setTrackStyleParam(trackList[[1]], "ylabgp", list(cex=.8, col="green"))
## set cex to avoid automatic adjust
setTrackStyleParam(trackList[[2]], "ylabgp", list(cex=.8, col="blue"))
setTrackStyleParam(trackList[[2]], "marginBottom", .2)
viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
plot data with adjusted color and size of y-axis label

plot data with adjusted color and size of y-axis label

The y-axis label can be put at the top or the bottom of each track.

setTrackStyleParam(trackList[[1]], "ylabpos", "bottomleft")
setTrackStyleParam(trackList[[2]], "ylabpos", "topright")
setTrackStyleParam(trackList[[2]], "marginTop", .2)
viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
plot data with adjusted y-axis label position

plot data with adjusted y-axis label position

For each transcript, the transcript name can be put on the upstream or downstream of the transcript.

trackListN <- trackList
setTrackStyleParam(trackListN[[3]], "ylabpos", "upstream")
setTrackStyleParam(trackListN[[4]], "ylabpos", "downstream")
## set cex to avoid automatic adjust
setTrackStyleParam(trackListN[[3]], "ylabgp", list(cex=.6))
setTrackStyleParam(trackListN[[4]], "ylabgp", list(cex=.6))
gr1 <- range(unname(unlist(GRangesList(sapply(trs, function(.ele) .ele$dat)))))
start(gr1) <- start(gr1) - 2000
end(gr1) <- end(gr1) + 2000
viewTracks(trackListN, gr=gr1, viewerStyle=viewerStyle)
plot data with adjusted transcripts name position

plot data with adjusted transcripts name position

Adjust the track color

The track color can be changed by setting the color slot in the style of each track. The first color is for the dat slot of track and the second color is for the dat2 slot.

setTrackStyleParam(trackList[[1]], "color", c("green", "black"))
setTrackStyleParam(trackList[[2]], "color", c("black", "blue"))
for(i in 3:length(trackList)) 
    setTrackStyleParam(trackList[[i]], "color", "black")
viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
plot data with adjusted track color

plot data with adjusted track color

Adjust the track height

The track height can be changed by setting the height slot in the style of each track. However, the total height for all the tracks should be 1.

trackListH <- trackList
setTrackStyleParam(trackListH[[1]], "height", .1)
setTrackStyleParam(trackListH[[2]], "height", .44)
for(i in 3:length(trackListH)){
    setTrackStyleParam(trackListH[[i]], "height", 
                       (1-(0.1+0.44))/(length(trackListH)-2))
}
viewTracks(trackListH, gr=gr, viewerStyle=viewerStyle)
plot data with adjusted track height

plot data with adjusted track height

Change the track names

The track names such as gene model names can be edited easily by changing the names of trackList.

names(trackList) <- c("cpsf160", "fox2", rep("HSPA8", 5))
viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
change the track names

change the track names

Show paired data in the same track

trackViewer can be used to show to-be-compared data in the same track side by side.

cpsf160 <- importScore(file.path(extdata, "cpsf160.repA_-.wig"),
                       file.path(extdata, "cpsf160.repB_-.wig"),
                       format="WIG")
strand(cpsf160$dat) <- strand(cpsf160$dat2) <- "-"
setTrackStyleParam(cpsf160, "color", c("black", "red"))
viewTracks(trackList(trs, cpsf160), gr=gr, viewerStyle=viewerStyle)
show two data in the same track

show two data in the same track

Flip the x-axis

The x-axis can be horizotally flipped for the genes in the negative strand.

viewerStyleF <- viewerStyle
setTrackViewerStyleParam(viewerStyleF, "flip", TRUE)
setTrackViewerStyleParam(viewerStyleF, "xaxis", TRUE)
setTrackViewerStyleParam(viewerStyleF, "margin", c(.1, .05, .01, .01))
vp <- viewTracks(trackList, gr=gr, viewerStyle=viewerStyleF)
addGuideLine(c(122929767, 122929969), vp=vp)
addArrowMark(list(x=122929650,
                  y=2),
             label="label",
             col="blue",
             vp=vp)
show data in the flipped track

show data in the flipped track

Optimize the theme

Currently, we support two themes: bw (black and white) and col (colored).

optSty <- optimizeStyle(trackList(repA, fox2, trs), theme="bw")
trackList <- optSty$tracks
viewerStyle <- optSty$style
vp <- viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
balck & white theme

balck & white theme

optSty <- optimizeStyle(trackList(repA, fox2, trs), theme="col")
trackList <- optSty$tracks
viewerStyle <- optSty$style
vp <- viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
colorful theme

colorful theme

optSty <- optimizeStyle(trackList(repA, fox2, trs), theme="safe")
trackList <- optSty$tracks
viewerStyle <- optSty$style
vp <- viewTracks(trackList, gr=gr, viewerStyle=viewerStyle)
safe theme

safe theme

Plot with breaks

We could plot the tracks with breaks by setting multiple genomic ranges.

gr.breaks <- GRanges("chr11", 
                     IRanges(c(122929275, 122929575, 122929775), 
                             c(122929555, 122929725, 122930122)), 
                     strand="-", percentage=c(.4, .2, .4))
vp <- viewTracks(trackList, gr=gr.breaks, viewerStyle=viewerStyle)
axis with breaks

axis with breaks

Session Info

sessionInfo()

R version 4.2.1 (2022-06-23) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu 20.04.5 LTS

Matrix products: default BLAS: /home/biocbuild/bbs-3.16-bioc/R/lib/libRblas.so LAPACK: /home/biocbuild/bbs-3.16-bioc/R/lib/libRlapack.so

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

attached base packages: [1] grid stats4 stats graphics grDevices utils datasets [8] methods base

other attached packages: [1] httr_1.4.4
[2] VariantAnnotation_1.44.0
[3] Rsamtools_2.14.0
[4] Biostrings_2.66.0
[5] XVector_0.38.0
[6] SummarizedExperiment_1.28.0
[7] MatrixGenerics_1.10.0
[8] matrixStats_0.62.0
[9] org.Hs.eg.db_3.16.0
[10] TxDb.Hsapiens.UCSC.hg19.knownGene_3.2.2 [11] GenomicFeatures_1.50.0
[12] AnnotationDbi_1.60.0
[13] Biobase_2.58.0
[14] Gviz_1.42.0
[15] rtracklayer_1.58.0
[16] trackViewer_1.34.0
[17] Rcpp_1.0.9
[18] GenomicRanges_1.50.0
[19] GenomeInfoDb_1.34.0
[20] IRanges_2.32.0
[21] S4Vectors_0.36.0
[22] BiocGenerics_0.44.0

loaded via a namespace (and not attached): [1] colorspace_2.0-3 rjson_0.2.21 deldir_1.0-6
[4] ellipsis_0.3.2 biovizBase_1.46.0 htmlTable_2.4.1
[7] base64enc_0.1-3 dichromat_2.0-0.1 rstudioapi_0.14
[10] bit64_4.0.5 fansi_1.0.3 xml2_1.3.3
[13] codetools_0.2-18 splines_4.2.1 cachem_1.0.6
[16] knitr_1.40 Formula_1.2-4 jsonlite_1.8.3
[19] cluster_2.1.4 dbplyr_2.2.1 png_0.1-7
[22] grImport_0.9-5 graph_1.76.0 BiocManager_1.30.19
[25] compiler_4.2.1 backports_1.4.1 lazyeval_0.2.2
[28] assertthat_0.2.1 Matrix_1.5-1 fastmap_1.1.0
[31] cli_3.4.1 htmltools_0.5.3 prettyunits_1.1.1
[34] tools_4.2.1 gtable_0.3.1 glue_1.6.2
[37] GenomeInfoDbData_1.2.9 dplyr_1.0.10 rappdirs_0.3.3
[40] jquerylib_0.1.4 rhdf5filters_1.10.0 vctrs_0.5.0
[43] xfun_0.34 stringr_1.4.1 lifecycle_1.0.3
[46] ensembldb_2.22.0 restfulr_0.0.15 InteractionSet_1.26.0
[49] XML_3.99-0.12 zlibbioc_1.44.0 scales_1.2.1
[52] BiocStyle_2.26.0 BSgenome_1.66.0 ProtGenerics_1.30.0
[55] hms_1.1.2 parallel_4.2.1 rhdf5_2.42.0
[58] AnnotationFilter_1.22.0 RColorBrewer_1.1-3 yaml_2.3.6
[61] curl_4.3.3 memoise_2.0.1 gridExtra_2.3
[64] ggplot2_3.3.6 sass_0.4.2 biomaRt_2.54.0
[67] rpart_4.1.19 latticeExtra_0.6-30 stringi_1.7.8
[70] RSQLite_2.2.18 highr_0.9 BiocIO_1.8.0
[73] plotrix_3.8-2 checkmate_2.1.0 filelock_1.0.2
[76] BiocParallel_1.32.0 rlang_1.0.6 pkgconfig_2.0.3
[79] bitops_1.0-7 evaluate_0.17 lattice_0.20-45
[82] Rhdf5lib_1.20.0 GenomicAlignments_1.34.0 htmlwidgets_1.5.4
[85] bit_4.0.4 tidyselect_1.2.0 magrittr_2.0.3
[88] R6_2.5.1 generics_0.1.3 Hmisc_4.7-1
[91] DelayedArray_0.24.0 DBI_1.1.3 pillar_1.8.1
[94] foreign_0.8-83 survival_3.4-0 KEGGREST_1.38.0
[97] RCurl_1.98-1.9 nnet_7.3-18 tibble_3.1.8
[100] crayon_1.5.2 interp_1.1-3 utf8_1.2.2
[103] BiocFileCache_2.6.0 rmarkdown_2.17 jpeg_0.1-9
[106] progress_1.2.2 data.table_1.14.4 Rgraphviz_2.42.0
[109] blob_1.2.3 digest_0.6.30 munsell_0.5.0
[112] bslib_0.4.0