Question: Analyzing a RNA-Seq dataset with paired samples and 2 batches
0
gravatar for ege.dedeoglu
21 days ago by
ege.dedeoglu0 wrote:

Hello everyone,

I have been tackling this issue for some time now. I have a dataset which has 20 patients and 40 samples (the samples are paired). There are two timepoints that I am checking differential expression for (eg. TP1 & TP2). However 11 of these patients are in one sequencing batch (B1) and the remaining 9 are in another batch (B2). In DESeq2 Vignette in the section "Model matrix not full rank" there is a part titled "Group-specific condition effects, individuals nested within groups" which explains how you can control for this type of situation. I have tried applying the design matrix as follows: (batch + batch:ba.pi + batch:condition) where ba.pi is the distinguisher of the individuals in the nested loop, but I got the same error. I have discarded 2 random samples from B1 and tried it again whicm made it work. I contemplated using ComBat to the normalized log2 converted DESeq2 counts seperately batch correct the two conditions and then do a t-test between the two groups but I am not sure if this is applicable. I would appreciate any input regarding this matter. The author's have seperately normalized the batches and worked with the common DEGs.

Regards,

Ege

sva deseq2 combat • 66 views
ADD COMMENTlink modified 21 days ago by Michael Love25k • written 21 days ago by ege.dedeoglu0
Answer: Analyzing a RNA-Seq dataset with paired samples and 2 batches
0
gravatar for Michael Love
21 days ago by
Michael Love25k
United States
Michael Love25k wrote:

If you have a design where you control patient baseline, and the paired samples per patient are all within batches, you don't need to control for batch effects with a term in the design. The patient baseline takes care of it.

ADD COMMENTlink written 21 days ago by Michael Love25k

I really could not understand what you are saying as English is not my first language. The batches contain 22 and 18 samples respectivly. Both the patients samples (TP1 & TP2) are in the same batch. My first design was :( ~patient + condition) but it gave me a distribution graph as the batches are seperated rather than the timepoints, however the time-points are seperated within the batches.

ADD REPLYlink modified 20 days ago • written 20 days ago by ege.dedeoglu0
Please log in to add an answer.

Help
Access

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