Compare different library
0
0
Entering edit mode
17 months ago
inminin • 0

Hello,

I have total RNA sequencing data (treatment 4, control 5), but we used different library sequencing for one treatment sample.

Can I do DESeq2 analysis with considering a batch effect for one sample?

sample treatment batch
s1 control A
s2 control A
s3 control A
s4 control A
s5 control A
s6 treatment A
s7 treatment A
s8 treatment A
s9 treatment A
s10 treatment B



dds <- DESeqDataSetFromMatrix(countData = countData,colData = colData,design = ~ batch + treatment)
dds$treatment <- factor(dds$treatment, levels=c("control", "treatment"))
dds$batch <- factor(dds$batch, levels=c("A", "B"))
dds <- DESeq(dds, full=design(dds), reduced = ~ batch)

Thanks,

RNA-Seq DESeq2 • 497 views
ADD COMMENT
1
Entering edit mode

In general, no, you cannot compare samples with different conditions in library prep.
RNA-direct? Like nanopore direct RNA-sequencing?

ADD REPLY
0
Entering edit mode

I agree with Wouter. You can but your results will be widely inaccurate because of all the compounding factors that arise when different sequencing methods are used. Can you repeat the controls using the same method as the treatment groups?

ADD REPLY
0
Entering edit mode

Thanks for your reply, I'll consider changing my method.

ADD REPLY
1
Entering edit mode

I added code markup to your post for increased readability. You can do this by selecting the text and clicking the 101010 button. When you compose or edit a post that button is in your toolbar, see image below:

101010 Button

ADD REPLY

Login before adding your answer.

Traffic: 2150 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6