# trackViewer Vignette: change the track styles

#### Jianhong Ou, Lihua Julie Zhu

#### 1 November 2022

Abstract

Visualize mapped reads along with annotation as track layers for NGS dataset such as ChIP-seq, RNA-seq, miRNA-seq, DNA-seq, SNPs and methylation data.

- Plot methylation/mutation/variant data by lollipop plot

- Plot high dense methylation/mutation/variant data by dandelion plot

- Plot genomic interaction data (HiC, cooler, etc) by viewTracks

## 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

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

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

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

## 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

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

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

## 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

## 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

## 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

### 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

## 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

## 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

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

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

## 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

## Session Info

    sessionInfo()

R version 4.2.1 (2022-06-23 ucrt) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows Server x64 (build 20348)

Matrix products: default

locale: [1] LC_COLLATE=C
[2] LC_CTYPE=English_United States.utf8
[3] LC_MONETARY=English_United States.utf8 [4] LC_NUMERIC=C
[5] LC_TIME=English_United States.utf8

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