Gviz: Should the add53 parameter react to reverseStrand = TRUE ?
2
0
Entering edit mode
stianlagstad ▴ 90
@stianlagstad-9723
Last seen 4.8 years ago

The code below:

transcriptID <- "ENST00000395443"

mart <- biomaRt::useMart(
  biomart = "ENSEMBL_MART_ENSEMBL",
  dataset = "hsapiens_gene_ensembl",
  host = "dec2013.archive.ensembl.org")

# Create transcript track
biomartTrack <- Gviz::BiomartGeneRegionTrack(
  filters = list(
    ensembl_transcript_id = transcriptID),
  biomart = mart)

# Create axis track
axisTrack <- Gviz::GenomeAxisTrack()
displayPars(axisTrack) <- list(
  add53 = TRUE,
  labelPos = "below",
  cex = 1.5
)

# works as expected
Gviz::plotTracks(
  c(biomartTrack, axisTrack),
  chromosome = "chr3")

Produces this plot:

And if I add reverseStrand=TRUE, I get this:

# add53 in same position?
Gviz::plotTracks(
  c(biomartTrack, axisTrack),
  chromosome = "chr3",
  reverseStrand = TRUE)

As you can see, everything is reversed except for the 5' and 3' indicators. As the arrows on the axis change, shouldn't the add53 parameter react to reverseStrand also? Or should I react to that by using add35 instead?

gviz • 1.3k views
ADD COMMENT
2
Entering edit mode
@florianhahnenovartiscom-3784
Last seen 6.3 years ago
Switzerland

yes. 

And sorry for the late reply. I was a bit lazy recently with Gviz stuff. Will add it on the list of fixes.

Thanks

ADD COMMENT
0
Entering edit mode

The code in my original post above still produce the exact same plots (i.e. the 5' and 3' indicators are wrong with reverseStrand=TRUE). Using Gviz version 1.16. Could you check that the fix worked?

ADD REPLY
1
Entering edit mode
@florianhahnenovartiscom-3784
Last seen 6.3 years ago
Switzerland

Looks like I forgot to commit this fix. Will become available shortly with Gviz 1.16.1 and 1.17.1

ADD COMMENT

Login before adding your answer.

Traffic: 695 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6