igvR 1.27.0
The UCSC Table Browser is a good source of genomic annotations of many different kinds. It has a clear, easily navigated user interface. It is a good complement to igvR.
The H3K4Me3 post-translational modification is frequently found in active promoters near transcription start sites. Here we obtain H3K4Me3 methylation marks in K562 cells in and around GATA2.
These are the steps involved:
All these steps are shown below.
library(igvR)
igv <- igvR()
setBrowserWindowTitle(igv, "H3K4Me3 GATA2")
setGenome(igv, "hg38")
showGenomicRegion(igv, "GATA2")
zoomOut(igv)
zoomOut(igv)
getGenomicRegion(igv)
$chrom
[1] "chr3"
$start
[1] 128454762
$end
[1] 128517865
$width
[1] 63104
$string
[1] "chr3:128,454,762-128,517,865"
With the output filename blank, the get output button shows you the selected data as text in your web browser:
Return to the previous UCSC Table Browser Screen, fill in a download filename, click get output
tbl <- read.table("~/drop/k562-h3k4me3-gata2.tsv", sep="\t", skip=1, as.is=TRUE, fill=TRUE)
colnames(tbl) <- c("chrom", "start", "end", "score")
Make sure the column classes are as expected:
lapply(tbl, class)
$chrom
[1] "character"
$start
[1] "integer"
$end
[1] "integer"
$score
[1] "numeric"
track <- DataFrameQuantitativeTrack("H3K4Me3 K562", tbl, autoscale=TRUE, color="darkGreen")
displayTrack(igv, track)
sessionInfo()
#> R Under development (unstable) (2025-01-20 r87609)
#> Platform: aarch64-apple-darwin20
#> Running under: macOS Ventura 13.7.1
#>
#> Matrix products: default
#> BLAS: /Library/Frameworks/R.framework/Versions/4.5-arm64/Resources/lib/libRblas.0.dylib
#> LAPACK: /Library/Frameworks/R.framework/Versions/4.5-arm64/Resources/lib/libRlapack.dylib; LAPACK version 3.12.0
#>
#> locale:
#> [1] C/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
#>
#> time zone: America/New_York
#> tzcode source: internal
#>
#> attached base packages:
#> [1] stats graphics grDevices utils datasets methods base
#>
#> other attached packages:
#> [1] BiocStyle_2.35.0
#>
#> loaded via a namespace (and not attached):
#> [1] digest_0.6.37 R6_2.5.1 bookdown_0.42 fastmap_1.2.0 xfun_0.50
#> [6] cachem_1.1.0 knitr_1.49 htmltools_0.5.8.1 rmarkdown_2.29 lifecycle_1.0.4
#> [11] cli_3.6.3 sass_0.4.9 jquerylib_0.1.4 compiler_4.5.0 tools_4.5.0
#> [16] evaluate_1.0.3 bslib_0.9.0 yaml_2.3.10 BiocManager_1.30.25 jsonlite_1.8.9
#> [21] rlang_1.1.5