[BioC] Gviz - Plot genes and data from - strand (3'-5') in 5'-3' direction

Steve Lianoglou lianoglou.steve at gene.com
Thu Sep 12 15:57:54 CEST 2013


Hi,

On Thu, Sep 12, 2013 at 6:14 AM, Robert Ivanek <robert.ivanek at unibas.ch> wrote:
> Hi Dominik,
>
> You can achieve that by changing the strand in corresponding track
> objects.
[snip]

While this might change the strand of a gene, I do not think this is
what the OP is really after -- as I have previously been keen to try
to implement the approach that (I think) the OP wants.

I believe the OP would like to plot "the mirror image" of the plot
that Gviz would "normally" show given the start/stop coordinates that
the plotting functionality infers (or is given).

The approach you suggest wouldn't properly "flip" my gene models and
the data associated with them -- it would only show them on the
positive strand. The 5'utr of the gene that was originally on the
negative strand would look as if it were the 3'utr of the gene on the
(now annotated) positive strand, and vice versa, and if I were (say)
plotting ChIP-seq data over this locus, the TF would appear to be
binding the 3'utr of a gene on the positive strand (or downstream of
it), when that's not really the case.

-steve

-- 
Steve Lianoglou
Computational Biologist
Bioinformatics and Computational Biology
Genentech



More information about the Bioconductor mailing list