R/epiRomics_region_tracks.R

In Huising-Lab/epiRomics: Epigenomic Analysis Package Built for R (epiRomics)

#' Visualizes chromatin availability at custom regions
#'
#' @param epiRomics_region GRanges of region to visualize
#' @param epiRomics_dB epiRomics class database containing all data initially loaded
#' @param epiRomics_track_connection data frame containing bigwig track locations and their names
#' @return GViz plot
#' @export
epiRomics_region_tracks <-
  function(epiRomics_region,
           epiRomics_dB,
           epiRomics_track_connection) {
    Gviz::GeneRegionTrack(eval(parse(text = epiRomics_dB@txdb)))
    chr <-
      base::as.character(BiocGenerics::unique(GenomeInfoDb::seqnames(epiRomics_region)))
    gtrack <- Gviz::GenomeAxisTrack()
    txTr <-
      Gviz::GeneRegionTrack(
        base::eval(base::parse(text = epiRomics_dB@txdb)),
        start = BiocGenerics::start(epiRomics_region),
        end = BiocGenerics::end(epiRomics_region),
        chromosome = chr,
        collapseTranscripts = TRUE,
        symbol = epiRomics_region$SYMBOL,
        name = epiRomics_region$SYMBOL,
        transcriptAnnotation = "symbol",
        geneSymbol = TRUE,
        howId = TRUE,
        min.height = 10,
        size = 10
      )
    #  txTr <- Gviz::GeneRegionTrack(base::eval(base::parse(text = epiRomics_dB@txdb)), start = as.data.frame(epiRomics_region)[,"start"], end=as.data.frame(epiRomics_region)[,"end"], chromosome = chr,
    #     symbol = epiRomics_region$SYMBOL, name = epiRomics_region$SYMBOL, transcriptAnnotation = "symbol",
    #    geneSymbol = TRUE, showId = TRUE)
    itrack <- Gviz::IdeogramTrack(genome = "mm10", chromosome = chr)
    epiRomics_track_connection_chromatin <-
      epiRomics_track_connection[epiRomics_track_connection$type == "atac", ]
    range_max <-
      maxCovFiles(epiRomics_track_connection_chromatin[, 1], epiRomics_region)
    max_cov <-
      plyr::round_any((range_max$X + (range_max$X * 0.1)), accuracy = 10, f = ceiling)
    chromatin_tracks <- list()
    for (i in 1:base::dim(epiRomics_track_connection_chromatin)[1]) {
      chromatin_tracks <- c(
        chromatin_tracks,
        Gviz::DataTrack(
          size = 10,
          epiRomics_track_connection_chromatin[i, 1],
          col = epiRomics_track_connection_chromatin[i, 3],
          fill = epiRomics_track_connection_chromatin[i, 3],
          col.histogram = epiRomics_track_connection_chromatin[i, 3],
          ylim = c(0, max_cov),
          range = epiRomics_track_connection_chromatin[i, 1],
          genome = epiRomics_dB@genome,
          type = "hist",
          chromosome = chr,
          strand = epiRomics_region@strand,
          name = epiRomics_track_connection_chromatin[i, 2]
        )
      )
    }
    if (table(epiRomics_track_connection$type == "rna")[[TRUE]] > 0) {
      rna_tracks <- list()
      epiRomics_track_connection_rna <-
        epiRomics_track_connection[epiRomics_track_connection$type == "rna", ]
      range_max <-
        maxCovFiles(epiRomics_track_connection_rna[, 1], epiRomics_region)
      max_cov <-
        plyr::round_any((range_max$X + (range_max$X * 0.1)), accuracy = 10, f = ceiling)
      for (i in 1:base::dim(epiRomics_track_connection_rna)[1]) {
        rna_tracks <- c(
          rna_tracks,
          Gviz::DataTrack(
            size = 10,
            epiRomics_track_connection_rna[i, 1],
            col = epiRomics_track_connection_rna[i, 3],
            fill = epiRomics_track_connection_rna[i, 3],
            col.histogram = epiRomics_track_connection_rna[i, 3],
            ylim = c(0, max_cov),
            range = epiRomics_track_connection_rna[i, 1],
            genome = epiRomics_dB@genome,
            type = "hist",
            chromosome = chr,
            strand = epiRomics_region@strand,
            name = epiRomics_track_connection_rna[i, 2]
          )
        )
      }
    }
    Gviz::plotTracks(base::c(itrack, gtrack, txTr, chromatin_tracks, rna_tracks))
  }