Answers.Body: <p>You can include that info in the design, though you'd have to decide whether to include it as a number of just a factor.  I've never heard of anyone including PCR cycles.  I've heard of people including RIN scores, though I've never done it.</p>
  <p>If the PCA looks alright, I wouldn't worry about.  I've got tons of examples of the lab people saying &quot;I did a few more cycles on these few samples&quot; and then you can't tell any difference in the PCA.  I think needing to do more cycles likely has more to do with variance in the prep than a real biological difference.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2020-09-30T20:05:45.560Z
Answers.Text: Thank you a lot for your reply!The problem is that PCA does not look good and I have observed that the samples with 15 cycles tend to cluster. I have tried to use removeBatchEffect from limma (I made a design matrix with cycle as a factor).then the pca got better. In the downstream analysis I used the new design matrix with cycle like in the question, and then the differentially expressed genes got dramatically reduced...(from 750 to 100).Any thoughts?
Answers.UserDisplayName: marilu
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2020-09-30T20:17:54.817Z
Answers.Text: It is possible to have some real DE genes even if your PCA looks discouraging.  100 is not a lot, but maybe that's realistic if your condition doesn't split nicely along your first few PCs. You might play around with making cycles a factor vs a number.  But if it's right that the number of cycles is really driving some difference, or representative of some real difference, accounting for it in the design will attenuate those differences.
Answers.UserDisplayName: swbarnes2
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2020-09-30T20:22:13.890Z
Answers.Text: Thank you so much for your reply, it really helps. One more question, I have observed that when I use vsd instead of rlog for PCA pot, PCA plot looks better. However, my samples are less than 30 (in the tutorial, rld is recommended especially for samples less than 30). Any ideas? Thank you in advance.
Answers.UserDisplayName: marilu
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2020-10-01T21:47:35.017Z
Answers.Text: I think the author recommends vst for everything now.
Answers.UserDisplayName: swbarnes2
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2020-09-30T19:56:17.207Z
Answers.LastActivityDate: 2020-09-30T19:56:17.207Z
Answers.UserDisplayName: swbarnes2
Body: <p>I am trying to understand library normalization in DESeq2. I would like to ask the following: I know that some samples have been run 15 cycles and some others 20, can I give this information to DESeq2, would it be useful? I mean, I know that DESEQ2 uses an algorithm to normalize library size, but since I have this piece of information, I wonder if it would be useful to provide it somehow to DESEQ2. If it is not, would it make sense to incorporate this information in the design matrix, something like: design(~cycle+condition)? Thank you a lot!</p>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-09-30T19:27:38.270Z
Comments.Text: What does "circle" mean in this context?  Do you mean some number of PCR cycles?
Comments.UserDisplayName: swbarnes2
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-09-30T19:49:51.117Z
Comments.Text: I am sorry, I meant cycles (I edited my question)!
Comments.UserDisplayName: marilu
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-09-30T19:51:19.097Z
Comments.Text: what is dds??? is it possible to ask 1 question
Comments.UserDisplayName: StupidWolf
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-09-30T19:54:23.943Z
Comments.Text: Ok, I have edited my question, only one question..I am sorry..
Comments.UserDisplayName: marilu
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-09-30T20:13:47.117Z
Comments.Text: dds is the name name of the DESeq2 object in pretty much every vignette and tutorial.
Comments.UserDisplayName: swbarnes2
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-09-30T21:37:50.940Z
Comments.Text: You can include it in the design matrix, but you need to be careful. If you use it as a continuous variable, then you are trying to find a linear relationship between counts and the cycle number. This may not be the case and in same cases you actually introduce more noise into the model
Comments.UserDisplayName: StupidWolf
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-09-30T21:38:20.617Z
Comments.Text: I think you are trying to make a model more complicated than it should be
Comments.UserDisplayName: StupidWolf
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-10-01T05:06:02.183Z
Comments.Text: Ok, so you recommend against..if I use it as a factor? Thank you!
Comments.UserDisplayName: marilu
Tags: <rna-seq><normalization><library>
Title: if I know the number of sequencing circles can I give this information to DESeq2?
ThreadText: # if I know the number of sequencing circles can I give this information to DESeq2?
  
  (Asked by: username_0 on 2020-09-30)
  
  I am trying to understand library normalization in DESeq2. I would like to ask the following: I know that some samples have been run 15 cycles and some others 20, can I give this information to DESeq2, would it be useful? I mean, I know that DESEQ2 uses an algorithm to normalize library size, but since I have this piece of information, I wonder if it would be useful to provide it somehow to DESEQ2. If it is not, would it make sense to incorporate this information in the design matrix, something like: design(~cycle+condition)? Thank you a lot!
  
  **Tags:** rna-seq, normalization, library
  
  > **Comments:**
  > **username_1**: What does "circle" mean in this context?  Do you mean some number of PCR cycles?  
  > **username_0**: I am sorry, I meant cycles (I edited my question)!  
  > **username_2**: what is dds??? is it possible to ask 1 question  
  > **username_0**: Ok, I have edited my question, only one question..I am sorry..  
  > **username_1**: dds is the name name of the DESeq2 object in pretty much every vignette and tutorial.  
  > **username_2**: You can include it in the design matrix, but you need to be careful. If you use it as a continuous variable, then you are trying to find a linear relationship between counts and the cycle number. This may not be the case and in same cases you actually introduce more noise into the model  
  > **username_2**: I think you are trying to make a model more complicated than it should be  
  > **username_0**: Ok, so you recommend against..if I use it as a factor? Thank you!  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_1 on 2020-09-30. Score: 2
  
  You can include that info in the design, though you'd have to decide whether to include it as a number of just a factor. I've never heard of anyone including PCR cycles. I've heard of people including RIN scores, though I've never done it.
  
  If the PCA looks alright, I wouldn't worry about. I've got tons of examples of the lab people saying "I did a few more cycles on these few samples" and then you can't tell any difference in the PCA. I think needing to do more cycles likely has more to do with variance in the prep than a real biological difference.
  
  > **Comments:**
  > **username_0**: Thank you a lot for your reply!The problem is that PCA does not look good and I have observed that the samples with 15 cycles tend to cluster. I have tried to use removeBatchEffect from limma (I made a design matrix with cycle as a factor).then the pca got better. In the downstream analysis I used the new design matrix with cycle like in the question, and then the differentially expressed genes got dramatically reduced...(from 750 to 100).Any thoughts?  
  > **username_1**: It is possible to have some real DE genes even if your PCA looks discouraging.  100 is not a lot, but maybe that's realistic if your condition doesn't split nicely along your first few PCs. You might play around with making cycles a factor vs a number.  But if it's right that the number of cycles is really driving some difference, or representative of some real difference, accounting for it in the design will attenuate those differences.  
  > **username_0**: Thank you so much for your reply, it really helps. One more question, I have observed that when I use vsd instead of rlog for PCA pot, PCA plot looks better. However, my samples are less than 30 (in the tutorial, rld is recommended especially for samples less than 30). Any ideas? Thank you in advance.  
  > **username_1**: I think the author recommends vst for everything now.  
  
  ---
  
  