plotMA function issue: can not circle genes

Question

I am trying to make MA-plot for my bulk-RNA-seq dataset and experiencing issues with it. I just copy-pasted the code that can be found here:

https://f1000research.com/articles/4-1070/v2

which is:

library("DESeq2")
plotMA(resLFC1, ylim=c(-5,5))
topGene <- rownames(resLFC1)[which.min(resLFC1$padj)]
with(resLFC1[topGene, ], {
  points(baseMean, log2FoldChange, col="dodgerblue", cex=2, lwd=2)
  text(baseMean, log2FoldChange, topGene, pos=2, col="dodgerblue")
})

However, I am not getting the right plot with the top gene circled, but the following one:

Maybe it is a bug in the plotMA function, not sure. How can I reproduce the figure from the publication?

Are you sure that resLFC1 is created as stated in the article ? There seems to be a lot more points than in the article. Also, do you get any error or message? — llrs, Apr 13 '18 at 10:23
Are you sure the weird double quotes (“ versus ") are not the problem? — benn, Apr 13 '18 at 15:22
I corrected them, of course. Rstudio complains about it right away. — Nikita Vlasenko, Apr 13 '18 at 18:39
Your plot shows all genes as DE (red points). That’s suspicious. — Konrad Rudolph, Apr 14 '18 at 17:44

score 4 · Accepted Answer · answered Apr 13 '18 at 21:07

It's likely that the fold-change of the topGene is greater than 5. Consequently, the circle is being drawn and the text produced, but they're outside the bounds of the plot (the bounds don't update if you add a circle and text label). Look at the log2FoldChange column of resLFC1[topGene, ] and increase the bounds to encompass that.

plotMA function issue: can not circle genes

1 Answers1