Question: Is it appropriate to apply RUVseq on output from kallisto?
gravatar for unawaz
7 months ago by
unawaz40 wrote:


I've ran kallisto on my samples, and summarized the transcript level counts to gene level counts to perform a differential gene expression analysis. I'm also interested in differential transcript expression, which I will perform using Sleuth.

Exploratory analysis at the gene level revealed that a few of the samples were clustering by batch. Previously I've applied RUVg on my counts, but the counts were produced by featureCounts and the analysis was only at the gene level. I wanted to know if it's appropriate to use RUVSeq with counts produced from kallisto? I will very likely need to apply this at the transcript level as well. We would also have to create in silico empirical control genes (which are normally created by performing a DE analysis).

Otherwise are there any alternative methods I can use to remove unwanted variation for data produced by kallisto?

ADD COMMENTlink modified 7 months ago by kristoffer.vittingseerup2.2k • written 7 months ago by unawaz40
gravatar for kristoffer.vittingseerup
7 months ago by
European Union
kristoffer.vittingseerup2.2k wrote:

It sounds like you already know the batch effect? If that is the case you should not use RUVg (or SVA) instead:

For DE: You simply add the batch effect as a factor in the DE model - then it will be ignored when you test for the difference caused by your experimental conditions.

For EDA: you need to remove the batch effect and RUV, SVA::combat or limma::removeBatchEffects can all help you with this. All methods work on Kallisto data (although each need different input) and both on gene and transcript level.

Btw if you are interested in transcript analysis you can, with the data you already have, directly analyse fx isoform switches - something my R package IsoformSwitchAnalyzeR can help you with (and it also handles batch effects). You can find examples of what type of analysis you can do in this section of the vignette.

ADD COMMENTlink written 7 months ago by kristoffer.vittingseerup2.2k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1440 users visited in the last hour