SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with batch entrez Tsuyoshi Bioinformatics 5 01-24-2018 01:12 AM
Batch effect for RNAseq data sehrrot Bioinformatics 5 07-07-2014 03:20 PM
hello and BLAST batch question SophieP Introductions 10 02-18-2014 03:20 AM
Batch effect Amative Bioinformatics 3 04-28-2013 04:01 PM
Batch effect, SNPs YKay Genomic Resequencing 0 12-13-2011 05:42 AM

Reply
 
Thread Tools
Old 06-30-2014, 06:34 AM   #1
emolinari
Member
 
Location: New haven

Join Date: May 2013
Posts: 47
Default Batch Effect

Hi guys,

I am currently struggling over an issue which I am not really sure how to solve.
In brief, I have run 4 sessions of PE seq on human cells (3 cohorts, Young, Old and Control). I have followed the Tophat-Cufflinks-Cuffmerge-Cuffdiff pipeline and visualized the data on Cummerbund.
When running the MDS plot I see 3 different clusters, and I can clearly assess that the samples cluster according to the sequencing session rather than the proper cohort. I've checked the quality of the data several times, and all the logs look ok. The only "weird" behavior is the properly mate reads rate, that for one group is on average 80%, for the other 72% and the third 60%.

Could this thing alone determine such a strange clustering or is rather a "batch effect"? Any suggestions for solving it?
BTW, I also have HTSeq counts of this data...do you thing I should use that on a different program, such as DeSeq or EdgeR???

Please help!!!

Manu
Attached Files
File Type: pdf Quartz 2 [*].pdf (498.3 KB, 32 views)
emolinari is offline   Reply With Quote
Old 06-30-2014, 06:57 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

That's a pretty classic batch effect. Use the SVA package (specifically, ComBat()) with DESeq2/edgeR/limma and you'll get more meaningful results. As to why this occurred, who knows. I've seen prominent library creation date batch effects before, so if the libraries were made on different dates then that could certainly be the original source of the problem.
dpryan is offline   Reply With Quote
Old 06-30-2014, 07:34 AM   #3
emolinari
Member
 
Location: New haven

Join Date: May 2013
Posts: 47
Default

Quote:
Originally Posted by dpryan View Post
That's a pretty classic batch effect. Use the SVA package (specifically, ComBat()) with DESeq2/edgeR/limma and you'll get more meaningful results. As to why this occurred, who knows. I've seen prominent library creation date batch effects before, so if the libraries were made on different dates then that could certainly be the original source of the problem.
Thanks dpryan for your comment,
The libraries where indeed prepped on different dates, and cluster accordingly (with the only exception of one sample, that was prepped alone one day and clusters with the big group on the right). I will try to do what you suggest...still you don't think that the properly paired mates rate has an influence in determining the clustering? I am asking because maybe there is a way to try to fix it...

Thanks again!
emolinari is offline   Reply With Quote
Old 06-30-2014, 07:50 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

That could be the cause as well (my guess would be that it's not, but that's just a guess). Just subset the alignments to contain only properly paired alignments and then look at the PCA plot. If the clustering goes away then you know that's the cause and will have also solved the problem (though you'd be throwing information away, so you still might get slightly better results using ComBat()).
dpryan is offline   Reply With Quote
Old 06-30-2014, 08:04 AM   #5
emolinari
Member
 
Location: New haven

Join Date: May 2013
Posts: 47
Default

Quote:
Originally Posted by dpryan View Post
That could be the cause as well (my guess would be that it's not, but that's just a guess). Just subset the alignments to contain only properly paired alignments and then look at the PCA plot. If the clustering goes away then you know that's the cause and will have also solved the problem (though you'd be throwing information away, so you still might get slightly better results using ComBat()).
I see what you mean, I'll try to that! A quick question though: should I use FPKM values or HTSeq?
emolinari is offline   Reply With Quote
Old 06-30-2014, 08:09 AM   #6
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

HTSeq (or featureCounts, which is faster).
dpryan is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:56 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO