# Obtain and Display H3K27ac K562 track from the AnnotationHub

Paul Shannon

#### 2025-04-15

#### Package

igvR 1.28.0

# Contents

# 1 Overview

The Bioconductor AnnotationHub is a good source of genomic annotations of many different kinds.

H3K27ac is an epigenetic modification to the histone H3, an acetylation of the lysine residue at N-terminal position 27. H3K27ac is associated with active enhancers.

To the best of my knowledge, fetching data from the AnnotationHub does not support regions. The fetch is necessarily of the entire genomic resource - all chromosomes - and so may require time-consuming downloads. Subsetting by region takes place after the often time-consuming download.

Therefore, to run this vignette for the first time may take up to 20 minutes due to that download time.

Once downloaded, however, the resource is cached.

# 2 Display a genomic region of interest in igvR

    library(igvR)
    library(AnnotationHub)

    igv <- igvR()
    setBrowserWindowTitle(igv, "H3K27ac GATA2")
    setGenome(igv, "hg19")
    showGenomicRegion(igv, "GATA2")
    for(i in 1:4) zoomOut(igv)

# 3 Query the AnnotationHub

    aHub <- AnnotationHub()
    query.terms <- c("H3K27Ac", "k562")
    length(query(aHub, query.terms))  # found 7
    h3k27ac.entries <- query(aHub, query.terms)

The available data, key and title:

      AH23388 | wgEncodeBroadHistoneK562H3k27acStdPk.broadPeak.gz
      AH29788 | E123-H3K27ac.broadPeak.gz
      AH30836 | E123-H3K27ac.narrowPeak.gz
      AH31772 | E123-H3K27ac.gappedPeak.gz
      AH32958 | E123-H3K27ac.fc.signal.bigwig
      AH33990 | E123-H3K27ac.pval.signal.bigwig
      AH39539 | E123-H3K27ac.imputed.pval.signal.bigwig

# 4 Select Two Resources: boadPeaks and fc bigwig

If not in your cache, this step may take 20 minutes.

    x.broadPeak <- aHub[["AH23388"]]
    x.bigWig <- aHub[["AH32958"]]

The two resources are different data types, requiring different processing to render as tracks in igvR

- x.broadPeak is a GRanges object in memory
- x.bigWig is a bigwig file in your cache

# 5 broadPeaks: subset and display

The broadPeak data is a GRanges object already in memory. Subset to obtain only the 252 kb region in which we are interested.

    roi <- getGenomicRegion(igv)
    gr.broadpeak <- x.broadPeak[seqnames(x.broadpeak)==roi$chrom &
                                start(x.broadpeak) > roi$start &
                                end(x.broadpeak) < roi$end]

igvR’s GrangesQuantitativeTrack must have only one numeric column in the GRanges metadata. That column is used as the magnitudes the track will display.

    names(mcols(gr.broadpeak))
      #  "name"        "score"       "signalValue" "pValue"      "qValue"
    mcols(gr.broadpeak) <- gr.broadpeak$score
    track <- GRangesQuantitativeTrack("h3k27ac bp", gr.broadpeak, autoscale=TRUE, color="brown")
    displayTrack(igv, track)

# 6 bigWig: subset and display

We use the import function from the rtracklayer package to read in only a small portion of the large bigWig file. Note that, as read, there is only one numeric metadata colum, “score”, so no reduction of mcols is needed.

    file.bigWig <- resource(x.bigWig)[[1]]
    gr.roi <- with(roi, GRanges(seqnames=chrom, IRanges(start, end)))
    gr.bw <- import(file.bigWig, which=gr.roi, format="bigWig")
    track <- GRangesQuantitativeTrack("h3k27ac.bw", gr.bw, autoscale=TRUE, color="gray")
    displayTrack(igv, track)

# 7 Session Info

    sessionInfo()
    #> R version 4.5.0 RC (2025-04-04 r88126 ucrt)
    #> Platform: x86_64-w64-mingw32/x64
    #> Running under: Windows Server 2022 x64 (build 20348)
    #>
    #> Matrix products: default
    #>   LAPACK version 3.12.1
    #>
    #> locale:
    #> [1] LC_COLLATE=C                           LC_CTYPE=English_United States.utf8    LC_MONETARY=English_United States.utf8
    #> [4] LC_NUMERIC=C                           LC_TIME=English_United States.utf8
    #>
    #> 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.36.0
    #>
    #> loaded via a namespace (and not attached):
    #>  [1] digest_0.6.37       R6_2.6.1            bookdown_0.43       fastmap_1.2.0       xfun_0.52
    #>  [6] cachem_1.1.0        knitr_1.50          htmltools_0.5.8.1   rmarkdown_2.29      lifecycle_1.0.4
    #> [11] cli_3.6.4           sass_0.4.10         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_2.0.0
    #> [21] rlang_1.1.6