This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RNA sequencing data batch effect removal

Hi All,

Recently I am trying to analyze some RNA sequencing data and perform the differential expression analysis. Since the sequencing data were generated in different times (experiments), I am worrying about the potential batch effects and would like to find a way to figure out and remove it (remove is more important!)

I tried to process the data in two different ways:

  1. Obtain the RSEM values and perform the differential expression analysis, in this way, can anyone kindly provide a good tool to apply on RSEM values to remove the batch effect?
  2. Obtain the htseq raw counts and use Deseq2 to perform the differential expression analysis. I saw previously another discussion about about justifying the batch effect by adding the 'batch' in the 'design' [design ~ batch+treatment] parameter in Deseq2. Besides this, can anyone kindly provide another good way to remove the batch effect?

Thanks for the help in advance.

rna-seq batch-effect
• 13,866 views • link
updated by Ram Ram • written by Tenghui Chen • Tenghui Chen

When you cluster the samples, or do PCA analysis, do you see batch effects as confounding ? I am just curious.

• 0 views • link
written by GouthamAtla GouthamAtla

Yes, I did PCA, but it did not really show very obvious batch effect, may because I just have very few sample size in each batch (from 2 to 5).

• 0 views • link
written by Tenghui Chen • Tenghui Chen

1 answer

Here is an example using "sva" package:

http://www.bioconductor.org/help/workflows/rnaseqGene/#batch

Remove Batch Effect From RNAseq with SVAseq and Combat

Also,you can use EdgeR.(Example by Devon).Its actually similar to using DESeq.

edgeR: Correct pipeline for DE analysis with multiple conditions and batches

• 0 views • link
written by Ron Ron

Thank you. I do not actually quite understand how to achieve a best practice with SVAseq. Do you know how to decide the n.sv value in the function?

• 0 views • link
written by Tenghui Chen • Tenghui Chen

Log in to answer this question.

More posts like this